Biotic interactions mediate the expansion of black mangrove (Avicennia germinans) into salt marshes under climate change

Many species are expanding their distributions to higher latitudes due to global warming. Understanding the mechanisms underlying these distribution shifts is critical for better understanding the impacts of climate changes. The climate envelope approach is widely used to model and predict species distribution shifts with changing climates. Biotic interactions between species, however, may also influence species distributions, and a better understanding of biotic interactions could improve predictions based solely on climate envelope models. Along the northern Gulf of Mexico coast, USA, subtropical black mangrove (Avicennia germinans) at the northern limit of its distribution grows sympatrically with temperate salt marsh plants in Florida, Louisiana, and Texas. In recent decades, freeze‐free winters have led to an expansion of black mangrove into salt marshes. We examined how biotic interactions between black mangrove and salt marsh vegetation along the Texas coast varied across (i) a latitudinal gradient (associated with a winter‐temperature gradient); (ii) the elevational gradient within each marsh (which creates different marsh habitats); and (iii) different life history stages of black mangroves (seedlings vs. juvenile trees). Each of these variables affected the strength or nature of biotic interactions between black mangrove and salt marsh vegetation: (i) Salt marsh vegetation facilitated black mangrove seedlings at their high‐latitude distribution limit, but inhibited black mangrove seedlings at lower latitudes; (ii) mangroves performed well at intermediate elevations, but grew and survived poorly in high‐ and low‐marsh habitats; and (iii) the effect of salt marsh vegetation on black mangroves switched from negative to neutral as black mangroves grew from seedlings into juvenile trees. These results indicate that the expansion of black mangroves is mediated by complex biotic interactions. A better understanding of the impacts of climate change on ecological communities requires incorporating context‐dependent biotic interactions into species range models.     

Winter climate change and coastal wetland foundation species: salt marshes vs. mangrove forests in the southeastern United States

We live in an era of unprecedented ecological change in which ecologists and natural resource managers are increasingly challenged to anticipate and prepare for the ecological effects of future global change. In this study, we investigated the potential effect of winter climate change upon salt marsh and mangrove forest foundation species in the southeastern United States. Our research addresses the following three questions: (1) What is the relationship between winter climate and the presence and abundance of mangrove forests relative to salt marshes; (2) How vulnerable are salt marshes to winter climate change‐induced mangrove forest range expansion; and (3) What is the potential future distribution and relative abundance of mangrove forests under alternative winter climate change scenarios? We developed simple winter climate‐based models to predict mangrove forest distribution and relative abundance using observed winter temperature data (1970–2000) and mangrove forest and salt marsh habitat data. Our results identify winter climate thresholds for salt marsh–mangrove forest interactions and highlight coastal areas in the southeastern United States (e.g., Texas, Louisiana, and parts of Florida) where relatively small changes in the intensity and frequency of extreme winter events could cause relatively dramatic landscape‐scale ecosystem structural and functional change in the form of poleward mangrove forest migration and salt marsh displacement. The ecological implications of these marsh‐to‐mangrove forest conversions are poorly understood, but would likely include changes for associated fish and wildlife populations and for the supply of some ecosystem goods and services.     

Climate warming is associated with smaller body size and shorter lifespans in moose near their southern range limit

Despite the importance of body size for individual fitness, population dynamics and community dynamics, the influence of climate change on growth and body size is inadequately understood, particularly for long‐lived vertebrates. Although temporal trends in body size have been documented, it remains unclear whether these changes represent the adverse impact of climate change (environmental stress constraining phenotypes) or its mitigation (via phenotypic plasticity or evolution). Concerns have also been raised about whether climate change is indeed the causal agent of these phenotypic shifts, given the length of time‐series analysed and that studies often do not evaluate – and thereby sufficiently rule out – other potential causes. Here, we evaluate evidence for climate‐related changes in adult body size (indexed by skull size) over a 4–decade period for a population of moose (Alces alces) near the southern limit of their range whilst also considering changes in density, predation, and human activities. In particular, we document: (i) a trend of increasing winter temperatures and concurrent decline in skull size (decline of 19% for males and 13% for females) and (ii) evidence of a negative relationship between skull size and winter temperatures during the first year of life. These patterns could be plausibly interpreted as an adaptive phenotypic response to climate warming given that latitudinal/temperature clines are often accepted as evidence of adaptation to local climate. However, we also observed: (iii) that moose with smaller skulls had shorter lifespans, (iv) a reduction in lifespan over the 4‐decade study period, and (v) a negative relationship between lifespan and winter temperatures during the first year of life. Those observations indicate that this phenotypic change is not an adaptive response to climate change. However, this decline in lifespan was not accompanied by an obvious change in population dynamics, suggesting that climate change may affect population dynamics and life‐histories differently.     

Influence of northern limit range on genetic diversity and structure in a widespread North American tree,sugar maple (Acer saccharum Marshall)

Due to climate change, the ranges of many North American tree species are expected to shift northward. Sugar maple (Acer saccharum Marshall) reaches its northern continuous distributional limit in northeastern North America at the transition between boreal mixed‐wood and temperate deciduous forests. We hypothesized that marginal fragmented northern populations from the boreal mixed wood would have a distinct pattern of genetic structure and diversity. We analyzed variation at 18 microsatellite loci from 23 populations distributed along three latitudinal transects (west, central, and east) that encompass the continuous–discontinuous species range. Each transect was divided into two zones, continuous (temperate deciduous) and discontinuous (boreal mixed wood), based on sugar maple stand abundance. Respective positive and negative relationships were found between the distance of each population to the northern limit (D_north), and allelic richness (A_R) and population differentiation (F_ST). These relations were tested for each transect separately; the pattern (discontinuous–continuous) remained significant only for the western transect. structure analysis revealed the presence of four clusters. The most northern populations of each transect were assigned to a distinct group. Asymmetrical gene flow occurred from the southern into the four northernmost populations. Southern populations in Québec may have originated from two different postglacial migration routes. No evidence was found to validate the hypothesis that northern populations were remnants of a larger population that had migrated further north of the species range after the retreat of the ice sheet. The northernmost sugar maple populations possibly originated from long‐distance dispersal.     

Warming and provenance limit tree recruitment across and beyond the elevation range of subalpine forest

Climate niche models project that subalpine forest ranges will extend upslope with climate warming. These projections assume that the climate suitable for adult trees will be adequate for forest regeneration, ignoring climate requirements for seedling recruitment, a potential demographic bottleneck. Moreover, local genetic adaptation is expected to facilitate range expansion, with tree populations at the upper forest edge providing the seed best adapted to the alpine. Here, we test these expectations using a novel combination of common gardens, seeded with two widely distributed subalpine conifers, and climate manipulations replicated at three elevations. Infrared heaters raised temperatures in heated plots, but raised temperatures more in the forest than at or above treeline because strong winds at high elevation reduced heating efficiency. Watering increased season‐average soil moisture similarly across sites. Contrary to expectations, warming reduced Engelmann spruce recruitment at and above treeline, as well as in the forest. Warming reduced limber pine first‐year recruitment in the forest, but had no net effect on fourth‐year recruitment at any site. Watering during the snow‐free season alleviated some negative effects of warming, indicating that warming exacerbated water limitations. Contrary to expectations of local adaptation, low‐elevation seeds of both species initially recruited more strongly than high‐elevation seeds across the elevation gradient, although the low‐provenance advantage diminished by the fourth year for Engelmann spruce, likely due to small sample sizes. High‐ and low‐elevation provenances responded similarly to warming across sites for Engelmann spruce, but differently for limber pine. In the context of increasing tree mortality, lower recruitment at all elevations with warming, combined with lower quality, high‐provenance seed being most available for colonizing the alpine, portends range contraction for Engelmann spruce. The lower sensitivity of limber pine to warming indicates a potential for this species to become more important in subalpine forest communities in the coming centuries.     

Variation in subsurface thermal characteristics of microrefuges used by range core and peripheral populations of the American pika (Ochotona princeps)

Microrefuges provide microclimates decoupled from inhospitable regional climate regimes that enable range‐peripheral populations to persist and are important to cold‐adapted species in an era of accelerated climate change. However, identifying and describing the thermal characteristics of microrefuge habitats is challenging, particularly for mobile organisms in cryptic, patchy habitats. We examined variation in subsurface thermal conditions of microrefuge habitats among different rock substrate types used by the American pika (Ochotona princeps ), a climate‐sensitive, rock‐dwelling Lagomorph. We compared subsurface temperatures in talus and lava substrates in pika survey sites in two US national park units; one park study area on the range periphery and the other in the range core. We deployed paired sensors to examine within‐site temperature variation. We hypothesized that subsurface temperatures within occupied sites and structurally complex substrates would be cooler in summer and warmer in winter than unoccupied and less complex sites. Although within‐site variability was high, with correlations between paired sensors as low as 47%, we found compelling evidence that pikas occupy microrefuge habitats where subsurface conditions provide more thermal stability than in unoccupied microhabitats. The percentage of days in which microhabitat temperatures were between ?2.5 and 25.5°C was significantly higher in occupied sites. Interestingly, thermal conditions were substantially more stable (p  < .05) in the lava substrate type identified to be preferentially used by pikas (pahoehoe vs. a'a) in a previous study. Our study and others suggest that thermal stability appears to be the defining characteristic of subsurface microrefuges used by American pikas and is a likely explanation for enigmatic population persistence at the range periphery. Our study exemplifies an integrated approach for studying complex microhabitat conditions, paired with site use surveys and contextualized with information about gene flow provided by complementary studies.     

A hierarchical Bayesian Beta regression approach to study the effects of geographical genetic structure and spatial autocorrelation on species distribution range shifts

Global climate change (GCC) may be causing distribution range shifts in many organisms worldwide. Multiple efforts are currently focused on the development of models to better predict distribution range shifts due to GCC. We addressed this issue by including intraspecific genetic structure and spatial autocorrelation (SAC) of data in distribution range models. Both factors reflect the joint effect of ecoevolutionary processes on the geographical heterogeneity of populations. We used a collection of 301 georeferenced accessions of the annual plant Arabidopsis thaliana in its Iberian Peninsula range, where the species shows strong geographical genetic structure. We developed spatial and nonspatial hierarchical Bayesian models (HBMs) to depict current and future distribution ranges for the four genetic clusters detected. We also compared the performance of HBMs with Maxent (a presence‐only model). Maxent and nonspatial HBMs presented some shortcomings, such as the loss of accessions with high genetic admixture in the case of Maxent and the presence of residual SAC for both. As spatial HBMs removed residual SAC, these models showed higher accuracy than nonspatial HBMs and handled the spatial effect on model outcomes. The ease of modelling and the consistency among model outputs for each genetic cluster was conditioned by the sparseness of the populations across the distribution range. Our HBMs enrich the toolbox of software available to evaluate GCC‐induced distribution range shifts by considering both genetic heterogeneity and SAC, two inherent properties of any organism that should not be overlooked.     

Effects of extreme weather events and legume presence on mycorrhization of Plantago lanceolata and Holcus lanatus in the field

Little is known about direct and indirect effects of extreme weather events on arbuscular mycorrhizal fungi (AMF) under field conditions. In a field experiment, we investigated the response of mycorrhization to drought and heavy rain in grassland communities. We quantified AMF biomass in soil, mycorrhization of roots of the grass Holcus lanatus and the forb Plantago lanceolata, as well as plant performance. Plants were grown in four‐species communities with or without a legume. We hypothesised that drought increases and heavy rain decreases mycorrhization, and that higher mycorrhization will be linked to improved stress resistance and higher biomass production. Soil AMF biomass increased under both weather extremes. Heavy rain generally benefitted plants and increased arbuscules in P. lanceolata. Drought neither reduced plant performance nor root mycorrhization. Arbuscules increased in H. lanatus several weeks after drought, and in P. lanceolata several weeks after heavy rain spells. These long‐lasting effects of weather events on mycorrhization highlight the indirect influence of climate on AMF via their host plant. Legume presence increased plant community biomass, but had only minor effects on mycorrhization. Arbuscule colonisation was negatively correlated with senescence during the dry summer. Mycorrhization and biomass production in P. lanceolata were positively related. However, increased mycorrhization was related to less biomass in the grass. AMF mycelium in soil might generally increase under extreme events, root colonisation, however, is host species specific. This might amplify community shifts in grassland under climate change by further increasing stress resistance of species that already benefit from changed precipitation.     

Climate‐driven range shifts explain the distribution of extant gene pools and predict future loss of unique lineages in a marine brown alga

The climate‐driven dynamics of species ranges is a critical research question in evolutionary ecology. We ask whether present intraspecific diversity is determined by the imprint of past climate. This is an ongoing debate requiring interdisciplinary examination of population genetic pools and persistence patterns across global ranges. Previously, contrasting inferences and predictions have resulted from distinct genomic coverage and/or geographical information. We aim to describe and explain the causes of geographical contrasts in genetic diversity and their consequences for the future baseline of the global genetic pool, by comparing present geographical distribution of genetic diversity and differentiation with predictive species distribution modelling (SDM) during past extremes, present time and future climate scenarios for a brown alga, Fucus vesiculosus. SDM showed that both atmospheric and oceanic variables shape the global distribution of intertidal species, revealing regions of persistence, extinction and expansion during glacial and postglacial periods. These explained the distribution and structure of present genetic diversity, consisting of differentiated genetic pools with maximal diversity in areas of long‐term persistence. Most of the present species range comprises postglacial expansion zones and, in contrast to highly dispersive marine organisms, expansions involved only local fronts, leaving distinct genetic pools at rear edges. Besides unravelling a complex phylogeographical history and showing congruence between genetic diversity and persistent distribution zones, supporting the hypothesis of niche conservatism, range shifts and loss of unique genetic diversity at the rear edge were predicted for future climate scenarios, impoverishing the global gene pool.     

Climate and host plant availability impact the future distribution of the bean leaf beetle (Cerotoma trifurcata)

The bean leaf beetle, Cerotoma trifurcata, has become a major pest of soybean throughout its North American range. With a changing climate, there is the potential for this pest to further expand its distribution and become an increasingly severe pest in certain regions. To examine this possibility, we developed bioclimatic envelope models for both the bean leaf beetle, and its most important agronomic host plant, soybean (Glycine max). These two models were combined to examine the potential future pest status of the beetle using climate change projections from multiple general circulation models (GCMs) and climate change scenarios. Despite the broad tolerances of soybean, incorporation of host plant availability substantially decreased the suitable and favourable areas for the bean leaf beetle as compared to an evaluation based solely on the climate envelope of the beetle, demonstrating the importance of incorporating biotic interactions in these predictions. The use of multiple GCM–scenario combinations also revealed differences in predictions depending on the choice of GCM, with scenario choice having less of an impact. While the Norwegian model predicted little northward expansion of the beetle from its current northern range limit of southern Ontario and overall decreases in suitable and favourable areas over time, the Canadian and Russian models predict that much of Ontario and Quebec will become suitable for the beetle in the future, as well as Manitoba under the Russian model. The Russian model also predicts expansion of the suitable and favourable areas for the beetle over time. Two predictions that do not depend on our choice of GCM include a decrease in suitability of the Mississippi Delta region and continued favourability of the southeastern United States.     

The effects of elevated CO2 and eutrophication on surface elevation gain in a European salt marsh

Salt marshes can play a vital role in mitigating the effects of global environmental change by dissipating incident storm wave energy and, through accretion, tracking increasing water depths consequent upon sea level rise. Atmospheric CO₂ concentrations and nutrient availability are two key variables that can affect the biological processes that contribute to marsh surface elevation gain. We measured the effects of CO₂ concentrations and nutrient availability on surface elevation change in intact mixed‐species blocks of UK salt marsh using six open‐top chambers receiving CO₂‐enriched (800 ppm) or ambient (400 ppm) air. We found more rapid surface elevation gain in elevated CO₂ conditions: an average increase of 3.4 mm over the growing season relative to ambient CO₂. Boosted regression analysis to determine the relative influence of different parameters on elevation change identified that a 10% reduction in microbial activity in elevated CO₂‐grown blocks had a positive influence on elevation. The biomass of Puccinellia maritima also had a positive influence on elevation, while other salt marsh species (e.g. Suaeda maritima) had no influence or a negative impact on elevation. Reduced rates of water use by the vegetation in the high CO₂ treatment could be contributing to elevation gain, either directly through reduced soil shrinkage or indirectly by decreasing microbial respiration rates due to lower redox levels in the soil. Eutrophication did not influence elevation change in either CO₂ treatment despite doubling aboveground biomass. The role of belowground processes (transpiration, root growth and decomposition) in the vertical adjustment of European salt marshes, which are primarily minerogenic in composition, could increase as atmospheric CO₂ concentrations rise and should be considered in future wetland models for the region. Elevated CO₂ conditions could enhance resilience in vulnerable systems such as those with low mineral sediment supply or where migration upwards within the tidal frame is constrained.     

Winter climate change and the poleward range expansion of a tropical invasive tree (Brazilian pepper—Schinus terebinthifolius)

Winter climate change is expected to lead to the tropicalization of temperate ecosystems, where tropical species expand poleward in response to a decrease in the intensity and duration of winter temperature extremes (i.e., freeze events). In the southeastern United States, freezing temperatures control the northern range limits of many invasive nonnative species. Here, we examine the influence of freezing temperatures and winter climate change on the northern range limits of an invasive nonnative tree—Schinus terebinthifolius (Brazilian pepper). Since introduction in the 1800s, Brazilian pepper has invaded ecosystems throughout south and central Florida to become the state's most widespread nonnative plant species. Although Brazilian pepper is sensitive to freezing temperatures, temperature controls on its northern distribution have not been adequately quantified. We used temperature and plant occurrence data to quantify the sensitivity of Brazilian pepper to freezing temperatures. Then, we examined the potential for range expansion under three alternative future climate scenarios (+2°C, +4°C, and +6°C). Our analyses identify a strong nonlinear sigmoidal relationship between minimum temperature and Brazilian pepper presence, with a discrete threshold temperature occurring near ?11°C. Our future scenario analyses indicate that, in response to warming winter temperatures, Brazilian pepper is expected to expand northward and transform ecosystems in north Florida and across much of the Gulf of Mexico and south Atlantic coasts of the United States. These results underscore the importance of early detection and rapid response efforts to identify and manage the northward invasion of Brazilian pepper in response to climate change. Looking more broadly, our work highlights the need to anticipate and prepare for the tropicalization of temperate ecosystems by tropical invasive species.     

Unravelling the effects of contemporary and historical range expansion on the distribution of genetic diversity in the damselfly Coenagrion scitulum

Although genetic diversity provides the basic substrate for evolution, there are a limited number of studies that assess the impact of recent climate change on intraspecific genetic variation. This study aims to unravel the degree to which historical and contemporary factors shape genetic diversity and structure across a large part of the range of the range‐expanding damselfly Coenagrion scitulum (Rambur, 1842). A total of 525 individuals from 31 populations were genotyped at nine microsatellites, and a subset was sequenced at two mitochondrial genes. We inferred the importance of geography, environmental factors, and recent range expansion on genetic diversity and structure. Genetic diversity decreased going westwards, suggesting a signature of historical post‐glacial expansion from east to west and the presence of eastern refugia. Although genetic differentiation decreased going northwards, it increased in the northern edge populations, suggesting a role of contemporary range expansion on the genetic make‐up of populations. The phylogeographical context was proven to be essential in understanding and identifying the genetic signatures of local contemporary processes. Within this framework, our results highlight that recent range expansion of a good disperser can decrease genetic diversity and increase genetic differentiation which should be considered when devising suitable conservation strategies.     

Pleistocene climate change and phylogeographic structure of the Gymnocarpos przewalskii (Caryophyllaceae) in the northwest China: Evidence from plastid DNA,ITS sequences,and Microsatellite

Northwestern China has a wealth of endemic species, which has been hypothesized to be affected by the complex paleoclimatic and paleogeographic history during Quaternary. In this paper, we used Gymnocarpos przewalskii as a model to address the evolutionary history and current population genetic structure of species in northwestern China. We employed two chloroplast DNA fragments (rps16 and psbB‐psbI), one nuclear DNA fragment (ITS), and simple sequence repeat (SSRs) to investigate the spatial genetic pattern of G. przewalskii. High genetic diversity (cpDNA: h_S = 0.330, h_T = 0.866; ITS: h_S = 0.458, h_T = 0.872) was identified in almost all populations, and most of the population have private haplotypes. Moreover, multimodal mismatch distributions were observed and estimates of Tajima's D and Fu's FS tests did not identify significantly departures from neutrality, indicating that recent expansion of G. przewalskii was rejected. Thus, we inferred that G. przewalskii survived generally in northwestern China during the Pleistocene. All data together support the genotypes of G. przewalskii into three groups, consistent with their respective geographical distributions in the western regions—Tarim Basin, the central regions—Hami Basin and Hexi Corridor, and the eastern regions—Alxa Desert and Wulate Prairie. Divergence among most lineages of G. przewalskii occurred in the Pleistocene, and the range of potential distributions is associated with glacial cycles. We concluded that climate oscillation during Pleistocene significantly affected the distribution of the species.     

Contrasting patterns of postglacial range shifts between the northern and southern hemisphere in Herbertus (Herbertaceae,Marchantiophyta)

The leafy liverwort genus Herbertus exhibits considerably variable morphology and widely disjunct distributions in both hemispheres. Here, the biogeographic history of the genus and its phylogenetic relationships with the focus on the taxonomically difficult, northern hemispheric disjunct species, were investigated. We conducted a time-calibrated, molecular-based phylogenetic analysis using psbA, trnL-F, and ITS1-2 loci and different approaches for ancestral range inference of the genus. Herbertus is inferred to have originated in the Cenozoic era about 51 million years ago, in an ancestral area including southern South America, the Neotropics, Oceania, and South-east Asia. The current distribution of the genus is supported to have resulted from long-term in situ persistence, short and long distance dispersals, extinctions and recolonizations. Postglacial range shifts between the southern and northern hemisphere in the genus show distinct patterns. In the southern hemisphere, Herbertus is characterized by in situ persistence, not showing further dispersal until the uplift of the Andean Cordillera. Species of the northern hemisphere showed wide range expansions and repeated recolonizations, including north- and southward dispersals, recolonizations and extinctions. Our results support that the ancestor of South-east Asiatic Herbertus had a Gondwanan origin and arrived in Asia via Indian Plate migration. The uplift of high mountains must have had a strong influence in the diversification and dispersal of the genus. Our results further suggest that climate changes must have had a profound effect on the evolution and biogeography of the species of Herbertus in the northern hemisphere, and might also have influenced the reproductive strategies of the genus. Few genetic differentiations amongst currently recognized species H. aduncus, H. dicranus, H. hutchinsiae, H. stramineus, H. delavayi, and H. kurzii, and amongst H. sendtneri, H. armitanus, and H. circinatus were shown, suggesting that the morphological characters that are currently used for delimiting species should be re-evaluated.     

Population structure of the South American Estuarine crab, Chasmagnathus granulatus (Brachyura: Varunidae) near the southern limit of its geographical distribution: comparison with northern populations

Crabs are among the most conspicuous and ecologically important invertebrates of the large intertidal zones that characterize estuarine and protected coastal areas in temperate regions. The habitat, population structure and breeding cycle of Chasmagnathus granulatus (Brachyura: Varunidae), a semiterrestrial burrowing crab endemic to the warm temperate coasts of the Southwestern Atlantic, were studied in San Antonio Bay (Argentina), near the southern limit of its range. San Antonio Bay has no freshwater input, winter is relatively colder, and summer warmer, than northern habitats of this species. Crabs lived both in vegetated and unvegetated zones, but density and sex ratio varied among dates and zones. The maximum observed density was 136 crabs/m² , the maximum carapace width (CW) was 32 mm (males) and 29.8 mm (females), ovigerous females were found only in November and January, and the smallest ovigerous female measured 17 mm CW. The population structure, spatial distribution, and recruitment pattern of C. granulatusdid not differ between San Antonio Bay and northern habitats. The higher density, smaller maximum size and shorter reproductive cycle observed in San Antonio cannot be atributed to changes associated with a latitudinal cline and other factors, such as thermal amplitude and food availability, need to be studied.     

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end‐glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land‐bridge route deriving from a “Carpathian” glacial refugium and one via a north‐eastern route from an “Eastern” glacial refugium near the Ural Mountains. Clustering of genome‐wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red‐backed voles (Myodes rutilus) took place in Fennoscandia just after end‐glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.