Climate change refugia and habitat connectivity promote species persistence
Climate Change Responses volume 4, Article number: 8 (2017)
Climate change refugia, areas buffered from climate change relative to their surroundings, are of increasing interest as natural resource managers seek to prioritize climate adaptation actions. However, evidence that refugia buffer the effects of anthropogenic climate change is largely missing.
Focusing on the climate-sensitive Belding’s ground squirrel (Urocitellus beldingi), we predicted that highly connected Sierra Nevada meadows that had warmed less or shown less precipitation change over the last century would have greater population persistence, as measured by short-term occupancy, fewer extirpations over the twentieth century, and long-term persistence measured through genetic diversity.
Across California, U. beldingi were more likely to persist over the last century in meadows with high connectivity that were defined as refugial based on a suite of temperature and precipitation factors. In Yosemite National Park, highly connected refugial meadows were more likely to be occupied by U. beldingi. More broadly, populations inhabiting Sierra Nevada meadows with colder mean winter temperatures had higher values of allelic richness at microsatellite loci, consistent with higher population persistence in temperature-buffered sites. Furthermore, both allelic richness and gene flow were higher in meadows that had higher landscape connectivity, indicating the importance of metapopulation processes. Conversely, anthropogenic refugia, sites where populations appeared to persist due to food or water supplementation, had lower connectivity, genetic diversity, and gene flow, and thus might act as ecological traps. This study provides evidence that validates the climate change refugia concept in a contemporary context and illustrates how to integrate field observations and genetic analyses to test the effectiveness of climate change refugia and connectivity.
Climate change refugia will be important for conserving populations as well as genetic diversity and evolutionary potential. Our study shows that in-depth modeling paired with rigorous fieldwork can identify functioning climate change refugia for conservation.
Attention on climate change refugia, areas relatively buffered from contemporary climate change over time that enable persistence of valued physical, ecological, and socio-cultural resources , has increased as scientists attempt to explain anomalies in the expected patterns of population responses to climate change [2,3,4,5]. Although managing climate change refugia is frequently identified as a key climate change adaptation strategy [1, 6], spatially explicit predictions of refugia and their effects on biodiversity are rare [5, 7, 8]. Most studies skip the key step of using independent data to evaluate whether these hypothesized refugia increase species persistence (Although see [9, 10]), which is necessary to confirm that investing limited resources into refugia is warranted .
We test whether highly connected montane meadows that have experienced little recent climate change  increase persistence of an associated species and thus act as refugia in the face of anthropogenic climate change. Montane meadows are key to alpine and subalpine animal communities, hydrological function, and human recreation, agriculture, and economy [11,12,13]. Recent work  used high resolution, dynamically downscaled climate models to map refugial meadow networks to identify where climate-sensitive species could persist, and restoration and conservation actions could be focused. Here we use data on site occupancy, site persistence, and population persistence as indicated by genetic diversity for the Belding’s ground squirrel. Urocitellus beldingi is a 250 g social, semi-fossorial mammal that feeds primarily on grasses and is commonly found in mid to high elevation meadows in western North America [14, 15]. Although ecologically important as key prey items and ecosystem engineers, U. beldingi in California have recently undergone large climate-related declines in occupancy, except where “anthropogenic refugia” were created by supplements with water or food through high human activity (such as campgrounds and municipal parks) . We predicted that populations located in highly connected refugial meadows identified by Maher et al.  on the basis of temperature and precipitation changes over the twentieth century would have higher contemporary occupancy, persistence, and, as a result, higher genetic diversity than sites that were not identified as climate change refugia. Our results validate the use of climate refugia as a key aspect of change climate adaptation.
Climate change refugia
Monthly climate data, specifically temperature, precipitation, and snow water equivalent (SWE), were derived using monthly PRISM data statistically downscaled to 800 m  and then, based on a hydrologic process model, downscaled further to 270 m . The following thresholds were used to define climate change refugia, based on comparing the historical period (1910–1939) to the recent period (1970–1999) [4 for more details]: (1) ≤1 °C changes in temperature (mean annual temperature, annual minimum temperature, annual maximum temperature, and annual mean temperature of the coldest quarter); (2) ≤10% change in relative annual mean precipitation and spring snowpack (yearly estimates of SWE); and (3) ≤2 months per year on average exceeding the extreme maximum (95th percentile) historical temperature and precipitation and the extreme minimum (95th percentile) historical precipitation variation.
Maher et al.  generated a set of friction surfaces to represent how landscape features may act as barriers or facilitators to species’ dispersal among Sierra Nevada meadows based upon four hypotheses: (1) isolation by distance; (2) isolation by topography, tested with a raster layer of distance weighted by elevation (upslope or downslope) between a layer of delineated meadows  with no dispersal permitted if the slope was greater than 45°; (3) isolation by watercourses in two forms: a) rivers and streams weighted with a high resistance value (100) compared to the surrounding matrix (“river_pres”), and b) distance from rivers and streams weighted as an increasing barrier to dispersal (“river_dist”); and (4) isolation by roads, considering distance from primary or secondary roads, which also acted as a proxy for human activity and presences. From these friction surfaces, maps of potential connectivity between meadow polygons buffered by 150 m were created using Circuitscape (ver 3.5.8; ). The resulting maps were continuous values of flow through a given pixel and represent the expected frequency of movement across a pixel. The approach was agnostic to species, such that a breadth of surfaces were examined, some of which are not applicable to the natural history of U. beldingi. In addition to these surfaces, we use two friction surfaces representing spatial variation in climatic water deficit (CWD) in the historical (“histCWD”) and modern (“modCWD”) periods, which allows a test of whether less variation in moisture availability facilitated dispersal.
Estimating site occupancy and persistence
A survey of U. beldingi occurrence in 38 meadows in Yosemite National Park was conducted in the summer of 2011 (Fig. 1, inset). Observers walked meadow transects while visually scanning for U. beldingi and listening for alarm calls; a site was considered occupied if a U. beldingi was seen or heard anywhere in the meadow. Meadows were buffered by 2 km using the R package rgeos  and mean climate change values were extracted from rasters using the package raster . We used a series of two-sample Wilcoxon tests to determine whether climate change refugia (as defined by the thresholds described above) and connectivity predicted increased probability of occupancy.
Next, as a measure of persistence, we examined results from resurveys (conducted 2003–2011) in or near 74 short-grass sites (either human-modified (Fig. 2a-b) or natural meadows (Fig. 2c)) in the Sierra Nevada, southern Cascades, and Modoc Plateau of California, USA (Fig. 1), where U. beldingi were historically (1902 to 1966) observed and/or trapped (based on detailed field notes from the University of California Museum of Vertebrate Zoology (MVZ)). As reported elsewhere , U. beldingi detectability per site averaged 0.96 and they were found to have been extirpated from 42% of the resurvey sites. To test whether climate change refugia were related to site persistence over the twentieth century, we used binomial tests based on the mean value of select climatic variables for each site with a 2 km radius around its centerpoint. Campgrounds, municipal parks, and agricultural fields were grouped as human-modified sites (the effect of connectivity on persistence of populations in them was analyzed separately).
There is a temporal mismatch between the modern climate window that Maher et al. 2017 used to define meadow refugia and the Belding’s ground squirrel resurvey data we used to test them. We believe this is justified because the measures of change in conditions need to precede the survey window in order to measure a potential impact of those changes, especially when trying to measure long-term persistence or extirpation that may show lag effects. If our climate-related variables shifted in magnitude and/or direction after 1999 and prior to the resurvey, then we would expect to find a poor relationship between the environmental data and resurvey data, assuming that climate-related effects are the predominant factor in persistence. Thus at the very least, datasets from Maher et al. (2017) represent a conservative estimate of climate changes relative to the resurvey.
Analysis of genetic diversity and gene flow
We sampled 244 U. beldingi trapped at 15 sites spread across its California range (Fig. 2) from 2003 to 2011 using Sherman or Tomahawk Live Traps, the majority in 2010 and 2011. Genetic samples were collected from 187 adults from an ear snip prior to release or from a liver biopsy in the case of MVZ voucher specimens (archived samples described in Additional file 1: Table S1), extracted using a salt extraction protocol , then stored at 4 °C until amplification.
In brief (see Additional file 2: Supplemental Methods for more details), 14 polymorphic microsatellite loci were amplified by PCR using fluorescently labelled primers and conditions optimized from protocols developed for other sciurids (see Additional file 3: Table S2). After adjusting for null alleles and linkage disequilibrium, we calculated the fixation index FST [24, 25], where lower values indicate higher gene flow, to test predictions of meadow connectivity and sample size-corrected allelic richness (Ar) to test the effects of both connectivity and refugia. We focused on allelic richness as a measure of genetic diversity; in both theory and practice, allelic richness is more sensitive to population reduction than is heterozygosity [26, 27].
structure 2.3.4  was used to examine genetic structuring and identify distinct genetic populations based on the 15 sampled sites (Fig. 2; Additional file 4: Table S3). Each run was replicated ten times with K set to run from 1 to 8 populations with an initial burn-in of 100,000 iterations followed by 106 Markov chain Monte Carlo repetitions. We summarized across runs using Structure Harvester , the statistic ΔK to choose the best K , Clumpp 1.1.2  to cluster results from structure run repetitions, and Distruct 1.1  to visualize population structure.
We examined the pattern of genetic diversity for 11 populations (n = 124) of U. beldingi in the central Sierra Nevada including Yosemite National Park (hereafter “Central Sierra”, Fig. 2) and the western Great Basin using linear regressions to test for correlations between allelic richness and minimal change in climate factors, as well as mean minimum temperature in the recent period.
To test correlations with the connectivity predictions, we calculated the minimum Euclidean distance between populations using the package sp.  in R (for isolation by distance) and determined the effective resistances among populations using Circuitscape (for the remaining connectivity hypotheses). We used the Central Sierra populations as focal nodes and the friction surfaces as resistance (topography, presence of watercourses, environmental variation) and conductance (distance from watercourses, distance from roads). The Multiple Matrix Regression with Randomization method  and backward selection were used to test whether a combination of distance matrices could explain the FST patterns. The advantage of this approach, as opposed to Mantel tests, is that it allowed the comparison of several distance matrices simultaneously, as well as the amount of variation explained by the given model. We compared connectivity models to patterns in FST among the Central Sierra sites. We also used linear regressions to test for correlations between mean connectivity and allelic richness.
Finally, BayesAss  was used to test for directional gene flow between central Sierra populations. It estimates the recent proportion of migrants among proximal populations and thus gives a different, and more precise, measure of gene flow than FST . The number of iterations was 108, with 106 as burn-in period, and a sampling frequency of 2000. Three randomly seeded replicates were run with Δp (allele freq) of 0.80, Δm (mig rate) of 0.80, and ΔF (inbreeding) of 0.80.
Principal component analysis
To summarize overall trends in climate change across the suite of environmental variables and since univariate refugia may be unlikely, we also conducted a principal components analysis (PCA) with data across all of our sites (n = 102). Because variance was unequal among variables, we used a correlation matrix, centering and normalizing all variables before calculations. Eigenvalues of the first three components of the PCA were greater than 1 and together explained 85.4% of the variation in climate change (Additional file 5: Table S4 and Additional file 6: Table S5), with sites separated out in general by increases in precipitation (PC1), increases in SWE vs. warming conditions (PC2), and wetter conditions with more months of extreme maximum temperatures (positive) vs. increases in minimum temperature (PC3).
For these first three components, we compared how well the summarized change explained patterns in occupancy, persistence, and genetic diversity. For occupancy and persistence, we used a series of logistic regressions to test the ability of each set of component scores to explain our observations. For genetic diversity, we compared allelic richness to each component score using linear models.
The 20 meadows occupied by U. beldingi in Yosemite National Park in 2011 experienced less change over the last century than the 18 unoccupied meadows based on the following five climate change indices (Fig. 3, Additional file 7: Table S6): change in annual minimum temperature (occupied meadow average = +1.30 °C, unoccupied = +1.66 °C; P < 0.001), change in mean annual precipitation (occupied = +8.38 mm, unoccupied = +50.06 mm; P < 0.001), number of months warmer than the historical extreme minimum temperature (occupied = 39.87, unoccupied = 62.94; P < 0.001), number of months wetter than the historical precipitation extreme (occupied = 25.45, unoccupied = 27.74; P = 0.007), and change in April 1st SWE (occupied = +11.68, unoccupied = +69.12; P < 0.001). According to these indices, U. beldingi were more likely to be found in climate change refugia. Occupied sites within Yosemite National Park also had higher values of connectivity than unoccupied sites in relation to specific predicted barriers (topography: P < 0.001; rivers: P river_pres = 0.007, W = 250, P river_dist = 0.021; roads: P road < 0.001; and climate-water deficit: WhistCWD = 310, P < 0.001, WmodCWD = 282, P = 0.001).
In addition, logistic regression for individual principal components of the PCA, which accounted for correlation among climate variables, strongly explained patterns in occupancy (all P < 0.05; estimated PC1 = 0.9271, estimated PC2 = −3.6038, estimated PC3 = 2.1670). Considering all components using a multiple logistic regression followed by stepwise removal leads to essentially equal support of patterns representing PC1 and PC3 (estimated PC1 = 0.6676, P = 0.054 and estimated PC3 = 1.5326, P = 0.055 when considered together). Logistic regressions of individual components indicated higher occupancy at cooler, drier sites (P < 0.05 in each case; estimated PC1 = 0.9271, estimated PC2 = −3.6038, estimated PC3 = 2.1670; Fig. 3, Additional File 8: Figure S1).
U. beldingi populations in California were also more likely to persist throughout the twentieth century in sites that were defined as climate change refugia from a suite of temperature and precipitation variables. Specifically, binomial tests showed that persistence was significantly correlated with refugial status (P < 0.05) for 7 of 8 climatic thresholds, with the exception of change in SWE (Fig. 4). Based on logistic regression, only PC1 predicted persistence (est. = 0.3598, P < 0.01).
Connectivity did not predict persistence of U. beldingi across all sampled sites over the last century in the Central Sierra, but did affect persistence when the 21 anthropogenic refugia were excluded from the analysis (Wtopo = 61, P < 0.0001). As with occupancy, each connectivity layer contributed to this effect (Wriver_pres = 65, P < 0.0001; Wriver_dist = 63, P < 0.0001; Wroad_road = 64, P < 0.0001; WhistCWD = 56, P < 0.0001; WmodCWD = 51, P < 0.0001; Fig. 5). In contrast, populations that persisted in anthropogenic refugia had the lowest connectivity values.
California populations of U. beldingi fell into four distinct genetic clusters, although the distinction between the northern and central sites could be due to lack of samples from intermediate populations. Mean allelic richness when adjusted for sample size was 3.3 ± 0.3 alleles per locus for the Central Sierra sites (range 2.28–3.85; Additional file 7: Table S3). Most sites where there were enough U. beldingi to trap for a robust genetic analysis ranked as refugia; thus a comparison of the effect of refugial status on genetic diversity was not possible for most climate indices. However, Central Sierra sites with lower current mean minimum temperature had higher values of allelic richness (β = −0.24, t9 = −4.48, R2 adj = 0.66, P = 0.002; Fig. 6) at temperature-buffered sites. As predicted, higher allelic richness of Central Sierra meadows was also correlated with higher meadow connectivity (β = 0.48, t9 = 2.92, R2 adj = 0.43, P = 0.02; Fig. 6), but a multiple regression suggested that mean minimum temperature contributed more strongly (βa = −0.28, t7 = −4.83, P = 0.002, βConn = 2.3e-06, t7 = 1.52, P > 0.1, βInt = 4.0e-07, t7 = 1.62; R2 adj = 0.82, F3,7; = 10.39, P = 0.006). There was no significant effect of either meadow area or change in environmental conditions (i.e., environmental change PC axes) on allelic richness per site. Meadow area also was not predictive of persistence at these sites (Morelli et al. 2012).
Gene flow (INTENDED AS A HEADER) across the entire region, pairwise FST values (mean = 0.092) ranged from 0.012 to 0.23. The highest FST values were between sites in the north and those in the central and southern region (mean FST = 0.18) and the lowest values were among the central sites (mean FST = 0.092). Analysis of FST across all sites indicated that higher connectivity, based on absence of a river barrier, was correlated with higher gene flow (βriver_pres = −0.011, βriver_dist = 133.8, F = 35.78, P = 0.001, R2 = 0.68). A more parameter-rich model supporting multiple connectivity hypotheses was found (βpathdist = −1.04e-05, βriver_pres = 20.76, βriver_dist = −0.027, βmodCWD = −1.20e-06, F = 61.69, P = 0.001, R2 = 0.89), but the extra coefficients, while significant, were close to unity. Focusing only on the Central Sierra sites, there was again support for watercourses as barriers to gene flow (βriver_pres = 54.99, F = 53.32, P = 0.001, R2 = 0.34). See Fig. 7 for an example that considered two equidistant pairs of sites, with higher connectivity predictions correlated with lower values of FST and thus higher levels of gene flow. This result was corroborated by the migrant rates estimated by BayesAss, with the lowest rates of migration into and out of anthropogenic refugia (Table 1).
This study reveals clear support that climate change refugia, here identified as areas that have changed little over the twentieth century in terms of temperature and precipitation, can predict persistence, with a focus on the climate-sensitive U. beldingi. Yosemite National Park meadows that were predicted to be refugial based on several climate change indices were more likely to contain U. beldingi. Across California, persistence over the last century was higher in sites that were more stable climatically, with the exception of sites that were artificially supplemented with food or water. Allelic richness, an index of effective metapopulation size and population persistence, was highest for the coolest sites across the central Sierra Nevada. Although sample size precluded a direct test of the correlation between change in temperature and precipitation with genetic diversity across sites, this indicates that refugial sites that warm less in the future will better maintain U. beldingi populations. Thus, we show that evidence from occupancy surveys can be combined with simple genetic data to test and validate the buffering qualities of climate change refugia.
We considered several different delineations of climate change refugia in this analysis in order to explore the idea of using a single species as a test of whether an area is a climate refugium. However, species undoubtedly differ in the climate factors that are most relevant for their persistence . Our previous research  indicated that climate change has resulted in increased body size in U. beldingi, perhaps related to longer growing seasons; thus less access to food does not appear to be driving population extirpations. We found that minimum temperature was a good predictor of persistence based on occupancy and genetic measures. Other studies have shown that mammals living in temperate montane systems could be negatively affected by warming in the winter [16, 39]. It is hypothesized that homeothermic species such as American pika (Ochotona princeps)  could freeze to death if the insulating snow layer disappears due to higher minimum temperatures in winter. This would be particularly relevant for hibernating species like U. beldingi that can lose necessary fat stores trying to maintain body temperature over winter. Shifts in winter timing also appear to be associated with declines in fitness of Columbian ground squirrels .
Conversely, it’s possible that our refugia measures were correlated with other variables and thus are just indirect predictors of persistence in the face of climate change. For example, predators, parasites, or habitat may be directly related to the change in climate, and U. beldingi may actually be responding to these other factors. Further study, and experimental manipulations, would be required to partition this variation.
Our study also provides support for biologically relevant habitat connectivity as a predictor of functioning climate change refugia. Among the connectivity hypotheses examined, the impact of watercourses appeared to have the most explanatory power for U. beldingi. As expected, we found that the small, semi-fossorial U. beldingi finds a river to be a substantial barrier to movement, whereas even the largest roads crossing the mid-elevations of the Sierra Nevada are not complete impediments to their dispersal. Genetic diversity and gene flow were also predicted by the habitat connectivity models.
Previous work showed that anthropogenic refugia had the highest densities of U. beldingi ; however, our study found they had the lowest gene flow and habitat connectivity values. For example, the Mono Lake County Park population was genetically distinct from the nearest Central Sierra Nevada populations (Fig. 2b); other known historical populations of U. beldingi had all been extirpated from the area . Thus, this research indicates that Mono Lake County Park may be acting as a refugium for a previously more ubiquitous, genetically distinct U. beldingi population in an area that is particularly low and warm for the species. Though there was almost no immigration into these sites, they may be acting as a source for nearby populations. However, this complex example highlights the need to distinguish anthropogenic refugia from ecological traps .
Climate change refugia will be important for conserving populations as well as genetic diversity and evolutionary potential . As during past epochs, they may act as safe havens of genetic diversity over the forthcoming decades and centuries of drastic change across the landscape, which is consistent with our finding across the full set of study sites that recent environmental shifts and isolation were negatively associated with persistence. As such, they are key targets with potential for management action [1, 44, 45], particularly while more long-term, sustainable climate change adaptation strategies are developed. Our study shows that in-depth modeling paired with rigorous fieldwork can identify functioning climate change refugia for conservation.
- FST :
Principal components analysis
Polymerase chain reaction
Morelli TL, Daly C, Dobrowski SZ, Dulen DM, Ebersole JL, Jackson ST, et al. Managing climate change Refugia for climate adaptation. PLoS One. 2016;11(8):e0159909. doi:10.1371/journal.pone.0159909.
Ashcroft MB. Identifying refugia from climate change. J Biogeogr. 2010;37(8):1407–13. doi:10.1111/j.1365-2699.2010.02300.x.
Lenoir J, Hattab T, Pierre G. Climatic microrefugia under anthropogenic climate change: implications for species redistribution. Ecography. 2017;40(2):253–66. doi:10.1111/ecog.02788.
Maher SP, Morelli TL, Hershey M, Flint AL, Flint LE, Moritz C, et al. Erosion of refugia in the Sierra Nevada meadows network with climate change. Ecosphere. 2017;8(4):e01673. doi:10.1002/ecs2.1673.
Isaak DJ, Young MK, Nagel DE, Horan DL, Groce MC. The cold-water climate shield: delineating refugia for preserving salmonid fishes through the 21st century. Glob Chang Biol. 2015;21(7):2540–53. doi:10.1111/gcb.12879.
Keppel G, Mokany K, Wardell-Johnson GW, Phillips BL, Welbergen JA, Reside AE. The capacity of refugia for conservation planning under climate change. Front Ecol Environ. 2015;13(2):106–12. doi:10.1890/140055.
McLaughlin BC, Zavaleta ES. Predicting species responses to climate change: demography and climate microrefugia in California valley oak (Quercus Lobata). Glob Chang Biol. 2012;18(7):2301–12. doi:10.1111/j.1365-2486.2011.02630.x. PubMed PMID: WOS:000304820300019
Franklin J, Davis FW, Ikegami M, Syphard AD, Flint LE, Flint AL, et al. Modeling plant species distributions under future climates: how fine scale do climate projections need to be? Glob Chang Biol. 2013;19(2):473–83. doi:10.1111/gcb.12051. PubMed PMID: WOS:000314219200013
Lourenço CR, Zardi GI, McQuaid CD, Serrão EA, Pearson GA, Jacinto R, et al. Upwelling areas as climate change refugia for the distribution and genetic diversity of a marine macroalga. J Biogeogr. 2016;43(8):1595–607. doi:10.1111/jbi.12744.
Maclean IMD, Hopkins JJ, Bennie J, Lawson CR, Wilson RJ. Microclimates buffer the responses of plant communities to climate change. Glob Ecol Biogeogr. 2015;24(11):1340–50. doi:10.1111/geb.12359.
Hammersmark CT, Dobrowski SZ, Rains MC, Mount JF. Simulated effects of stream restoration on the distribution of wet-meadow vegetation. Restor Ecol. 2010;18(6):882–93. doi:10.1111/j.1526-100X.2009.00519.x. PubMed PMID: WOS:000283718000012
Kuhn TJ, Safford HD, Jones BE, Tate KW. Aspen (Populus Tremuloides) stands and their contribution to plant diversity in a semiarid coniferous landscape. Plant Ecol. 2011;212(9):1451–63. doi:10.1007/s11258-011-9920-4. PubMed PMID: WOS:000293344100005
Roche LM, Latimer AM, Eastburn DJ, Tate KW. Cattle grazing and conservation of a meadow-dependent amphibian species in the sierra Nevada. PLoS One. 2012;7(4) doi:10.1371/journal.pone.0035734. PubMed PMID: WOS:000305345200072
Grinnell J, Dixon J. Natural history of the ground squirrels of California. Sacramento Cal State Comm Horticulture Monthly Bul. 1919;7:597–708. PubMed PMID: ZOOREC:ZOOR05600003484
Sherman PW, Morton ML. Demography of Belding's ground squirrels. Ecology. 1984;65(5):1617–28.
Morelli TL, Smith AB, Kastely CR, Mastroserio I, Moritz C, Beissinger SR. Anthropogenic refugia ameliorate the severe climate-related decline of a montane mammal along its trailing edge. Proceedings of the Royal Society B-Biological Sciences. 2012;279(1745):4279–86. doi:10.1098/rspb.2012.1301. PubMed PMID: WOS:000308748000024
Daly C, Halbleib M, Smith JI, Gibson WP, Doggett MK, Taylor GH, et al. Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States. Int J Climatol. 2008;28(15):2031–64. doi:10.1002/joc.1688.
Flint AL, Flint LE. Downscaling future climate scenarios to fine scales for hydrologic and ecological modeling and analysis. Ecol Process. 2012;1:2.
Sierra Nevada Multi-Source Meadow Polygons Compilation (v 1.0) [Internet]. Center for Watershed Sciences, UC Davis. 2012 [cited December 2012]. Available from: https://meadows.ucdavis.edu/.
McRae BH, Dickson BG, Keitt TH, Shah VB. Using circuit theory to model connectivity in ecology, evolution, and conservation. Ecology. 2008;89(10):2712–24. doi:10.1890/07-1861.1.
Bivand R, Rundel C. rgeos: interface to geometry engine - open source (GEOS). 0.3–1 ed2013.
Hijmans RJ, van Etten J. raster: Geographic data analysis and modeling. R package version 2.0–31 ed2012.
Aljanabi SM, Martinez I. Universal and rapid saltextraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 1997;25:4692–3.
Hedrick PW. A standardized genetic differentiation measure. Evolution. 2005;59(8):1633–8. doi:10.1111/j.0014-3820.2005.tb01814.x. PubMed PMID: WOS:000231658900003
Weir BS, Cockerham CC. Estimating F-Statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70. doi:10.2307/2408641. PubMed PMID: WOS:A1984TY40400017
Nei M, Maruyama T, Chakraborty R. The bottleneck effect and genetic variability of populations. Evolution. 1975;29:1–10.
Rubidge EM, Patton JL, Lim M, Burton AC, Brashares JS, Moritz C. Climate-induced range contraction drives genetic erosion in an alpine mammal. Nat Clim Chang. 2012;2(4):285–8. doi:10.1038/nclimate1415.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59. PubMed PMID: WOS:000087475100039
Earl DA, vonHoldt BM. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61. doi:10.1007/s12686-011-9548-7. PubMed PMID: WOS:000303536400036
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20. doi:10.1111/j.1365-294X.2005.02553.x. PubMed PMID: WOS:000229961500029
Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23(14):1801–6. doi:10.1093/bioinformatics/btm233. PubMed PMID: WOS:000249248300012
Rosenberg N. DISTRUCT: A program for the graphical display of population structure. Mol Ecol Notes. 2004;4:137–8.
Pebesma EJ, Bivand RS. Classes and methods for spatial data in R. R News 5 (2)2005.
Wang IJ. Examining the full effects of landscape heterogeneity on spatial genetic variation: a multiple matrix regression approach for quantifying geographic and ecological isolation. Evolution. 2013;67(12):3403–11. doi:10.1111/evo.12134.
Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163(3):1177–91. PubMed PMID: WOS:000182046900029
Whitlock MC, McCauley DE. Indirect measures of gene flow and migration: FST not equal to 1/(4Nm + 1). (0018-067X (Print)).
Rapacciuolo G, Maher S, Schneider A, Hammond T, Jabis M, Walsh R, et al. Beyond a global warming fingerprint of climate change: heterogeneous biogeographic responses in California. Glob Chang Biol. 2014;20(9):2841–55.
Eastman LM, Morelli TL, Rowe KC, Conroy CJ, Moritz C. Size increase in high elevation ground squirrels over the last century. Glob Chang Biol. 2012;18(5):1499–508. doi:10.1111/j.1365-2486.2012.02644.x. PubMed PMID: WOS:000302543500003
Williams CM, HHA L, Sinclair BJ. Cold truths: how winter drives responses of terrestrial organisms to climate change. Biol Rev. 2015;90:214–35.
Beever EA, Ray C, Mote PW, Wilkening JL. Testing alternative models of climate-mediated extirpations. Ecol Appl. 2010;20(1):164–78.
Lane JE, Kruuk LEB, Charmantier A, Murie JO, Dobson FS. Delayed phenology and reduced fitness associated with climate change in a wild hibernator. Nature. 2012;489:554–8.
Schlaepfer MA, Runge MC, Sherman PW. Ecological and evolutionary traps. Trends in Ecology & Evolution. 2002;17(10):474–80. PubMed PMID: ISI:000178075900011
Moritz C, Agudo R. The future of species under climate change: resilience or decline? Science. 2013;341(6145):504–8. doi:10.1126/science.1237190.
Shoo LP, Hoffmann AA, Garnett S, Pressey RL, Williams YM, Taylor M, et al. Making decisions to conserve species under climate change. Clim Chang. 2013;119(2):239–46. doi:10.1007/s10584-013-0699-2. PubMed PMID: WOS:000321955100002
Keppel G, Mokany K, Wardell-Johnson GW, Phillips BL, Welbergen JA, Reside AE. The capacity of refugia for conservation planning under climate change. Front Ecol Environ. 2015; doi:10.1890/140055.
We are grateful to J Patton, K Rowe, A Smith, and other members of the Grinnell Resurvey Project for access to samples. We thank the faculty, staff, students, and friends who helped with data collection and analysis. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.
This work was financially supported by a grant from the California Landscape Conservation Cooperative (to TLM, SRB, and CM), an NSF Bioinformatics Postdoctoral Research Fellowship (to TLM), and the DOI Northeast Climate Science Center.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its Additional files. Maps of refugial meadows with connectivity rankings are available through https://databasin.org.
Ethics approval and consent to participate
Research permits from California Fish and Game (SC-10774) and Yosemite National Park Service (YOSE-2011-SCI-0085) and IACUC approval R278–0412 facilitated this research.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of sampled museum specimens. (PDF 42 kb)
Supplemental Methods. Details of genetic analysis. (PDF 42 kb)
Information on primers used to analyze Urocitellus beldingi genetic diversity and gene flow. (PDF 15 kb)
Summary of genetic information for each population corresponding with Fig. 2. (PDF 15 kb)
Eigenvalues of the first three principal components. (PDF 10 kb)
Loadings from the principal components analysis. (PDF 10 kb)
Climate variables defining climate change refugia predicted the probability of U. beldingi occupancy in Yosemite National Park. (PDF 14 kb)
Results of principal components analysis for Yosemite National Park. (PDF 99 kb)
About this article
Cite this article
Morelli, T.L., Maher, S.P., Lim, M.C.W. et al. Climate change refugia and habitat connectivity promote species persistence. Clim Chang Responses 4, 8 (2017). https://doi.org/10.1186/s40665-017-0036-5