Integrating mechanistic and empirical model projections to assess climate impacts on tree species distributions in northwestern North America

Empirical and mechanistic models have both been used to assess the potential impacts of climate change on species distributions, and each modeling approach has its strengths and weaknesses. Here, we demonstrate an approach to projecting climate‐driven changes in species distributions that draws on both empirical and mechanistic models. We combined projections from a dynamic global vegetation model (DGVM) that simulates the distributions of biomes based on basic plant functional types with projections from empirical climatic niche models for six tree species in northwestern North America. These integrated model outputs incorporate important biological processes, such as competition, physiological responses of plants to changes in atmospheric CO₂ concentrations, and fire, as well as what are likely to be species‐specific climatic constraints. We compared the integrated projections to projections from the empirical climatic niche models alone. Overall, our integrated model outputs projected a greater climate‐driven loss of potentially suitable environmental space than did the empirical climatic niche model outputs alone for the majority of modeled species. Our results also show that refining species distributions with DGVM outputs had large effects on the geographic locations of suitable habitat. We demonstrate one approach to integrating the outputs of mechanistic and empirical niche models to produce bioclimatic projections. But perhaps more importantly, our study reveals the potential for empirical climatic niche models to over‐predict suitable environmental space under future climatic conditions.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Assessment of the relationships of geographic variation in species richness to climate and landscape variables within and among lineages of North American freshwater fishes

Aim Geographic variation in species richness is a well‐studied phenomenon. However, the unique response of individual lineages to environmental gradients in the context of general patterns of biodiversity across broad spatial scales has received limited attention. The focus of this research is to examine relationships between species richness and climate, topographic heterogeneity and stream channel characteristics within and among families of North American freshwater fishes. Location The United States and Canada. Methods Distribution maps of 828 native species of freshwater fishes were used to generate species richness estimates across the United States and Canada. Variation in species richness was predicted using spatially explicit models incorporating variation in climate, topography and/or stream channel length and stream channel diversity for all 828 species as well as for the seven largest families of freshwater fishes. Results The overall gradient of species richness in North American freshwater fishes is best predicted by a model incorporating variables describing climate and topography. However, the response of species richness to particular climate or landscape variables differed among families, with models possessing the highest predictive ability incorporating data on climate, topography and/or stream channel characteristics within a region. Main conclusions The correlations between species richness and abiotic variables suggest a strong influence of climate and physical habitat on the structuring of regional assemblages of North American freshwater fishes. However, the relationship between these variables and species richness varies among families, suggesting the importance of phylogenetic constraints on the regulation of geographic distributions of species.     

Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean‐temperate tree species

Populations occurring in areas of overlap between the current and future distribution of a species are particularly important because they can represent “refugia from climate change”. We coupled ecological and range‐wide genetic variation data to detect such areas and to evaluate the impacts of habitat suitability changes on the genetic diversity of the transitional Mediterranean‐temperate tree Fraxinus angustifolia. We sampled and genotyped 38 natural populations comprising 1006 individuals from across Europe. We found the highest genetic diversity in western and northern Mediterranean populations, as well as a significant west to east decline in genetic diversity. Areas of potential refugia that correspond to approximately 70% of the suitable habitat may support the persistence of more than 90% of the total number of alleles in the future. Moreover, based on correlations between Bayesian genetic assignment and climate, climate change may favour the westward spread of the Black Sea gene pool in the long term. Overall, our results suggest that the northerly core areas of the current distribution contain the most important part of the genetic variation for this species and may serve as in situ macrorefugia from ongoing climate change. However, rear‐edge populations of the southern Mediterranean may be exposed to a potential loss of unique genetic diversity owing to habitat suitability changes unless populations can persist in microrefugia that have facilitated such persistence in the past.     

The implication of irrigation in climate change impact assessment: a European‐wide study

This study evaluates the impacts of projected climate change on irrigation requirements and yields of six crops (winter wheat, winter barley, rapeseed, grain maize, potato, and sugar beet) in Europe. Furthermore, the uncertainty deriving from consideration of irrigation, CO₂ effects on crop growth and transpiration, and different climate change scenarios in climate change impact assessments is quantified. Net irrigation requirement (NIR) and yields of the six crops were simulated for a baseline (1982–2006) and three SRES scenarios (B1, B2 and A1B, 2040–2064) under rainfed and irrigated conditions, using a process‐based crop model, SIMPLACE . We found that projected climate change decreased NIR of the three winter crops in northern Europe (up to 81 mm), but increased NIR of all the six crops in the Mediterranean regions (up to 182 mm yr^?1). Climate change increased yields of the three winter crops and sugar beet in middle and northern regions (up to 36%), but decreased their yields in Mediterranean countries (up to 81%). Consideration of CO₂ effects can alter the direction of change in NIR for irrigated crops in the south and of yields for C3 crops in central and northern Europe. Constraining the model to rainfed conditions for spring crops led to a negative bias in simulating climate change impacts on yields (up to 44%), which was proportional to the irrigation ratio of the simulation unit. Impacts on NIR and yields were generally consistent across the three SRES scenarios for the majority of regions in Europe. We conclude that due to the magnitude of irrigation and CO₂ effects, they should both be considered in the simulation of climate change impacts on crop production and water availability, particularly for crops and regions with a high proportion of irrigated crop area.     

Mutualism influences species distribution predictions for a bromeliad‐breeding anuran under climate change

Ecological niche models, or species distribution models, have been widely used to identify potentially suitable areas for species in future climate change scenarios. However, there are inherent errors to these models due to their inability to evaluate species occurrence influenced by non‐climatic factors. With the intuit to improve the modelling predictions for a bromeliad‐breeding treefrog (Phyllodytes melanomystax, Hylidae), we investigate how the climatic suitability of bromeliads influences the distribution model for the treefrog in the context of baseline and 2050 climate change scenarios. We used point occurrence data on the frog and the bromeliad (Vriesea procera, Bromeliaceae) to generate their predicted distributions based on baseline and 2050 climates. Using a consensus of five algorithms, we compared the accuracy of the models and the geographic predictions for the frog generated from two modelling procedures: (i) a climate‐only model for P. melanomystax and V. procera; and (ii) a climate‐biotic model for P. melanomystax, in which the climatic suitability of the bromeliad was jointly considered with the climatic variables. Both modelling approaches generated strong and similar predictive power for P. melanomystax, yet climate‐biotic modelling generated more concise predictions, particularly for the year 2050. Specifically, because the predicted area of the bromeliad overlaps with the predictions for the treefrog in the baseline climate, both modelling approaches produce reasonable similar predicted areas for the anuran. Alternatively, due to the predicted loss of northern climatically suitable areas for the bromeliad by 2050, only the climate‐biotic models provide evidence that northern populations of P. melanomystax will likely be negatively affected by 2050.     

Multilocus phylogeography of a widespread savanna–woodland‐adapted rodent reveals the influence of Pleistocene geomorphology and climate change in Africa's Zambezi region

Understanding historical influences of climate and physiographic barriers in shaping patterns of biodiversity remains limited for many regions of the world. For mammals of continental Africa, phylogeographic studies, particularly for West African lineages, implicate both geographic barriers and climate oscillations in shaping small mammal diversity. In contrast, studies for southern African species have revealed conflicting phylogenetic patterns for how mammalian lineages respond to both climate change and geologic events such as river formation, especially during the Pleistocene. However, these studies were often biased by limited geographic sampling or exclusively focused on large‐bodied taxa. We exploited the broad southern African distribution of a savanna–woodland‐adapted African rodent, Gerbilliscus leucogaster (bushveld gerbil) and generated mitochondrial, autosomal and sex chromosome data to quantify regional signatures of climatic and vicariant biogeographic phenomena. Results indicate the most recent common ancestor for all G. leucogaster lineages occurred during the early Pleistocene. We documented six divergent mitochondrial lineages that diverged ~0.270–0.100 mya, each of which was geographically isolated during periods characterized by alterations to the course of the Zambezi River and its tributaries as well as regional ‘megadroughts’. Results demonstrate the presence of a widespread lineage exhibiting demographic expansion ~0.065–0.035 mya, a time that coincides with savanna–woodland expansion across southern Africa. A multilocus autosomal perspective revealed the influence of the Kafue River as a current barrier to gene flow and regions of secondary contact among divergent mitochondrial lineages. Our results demonstrate the importance of both climatic fluctuations and physiographic vicariance in shaping the distribution of southern African biodiversity.     

No‐analog climates and shifting realized niches during the late quaternary: implications for 21st‐century predictions by species distribution models

Empirically derived species distributions models (SDMs) are increasingly relied upon to forecast species vulnerabilities to future climate change. However, many of the assumptions of SDMs may be violated when they are used to project species distributions across significant climate change events. In particular, SDM's in theory assume stable fundamental niches, but in practice, they assume stable realized niches. The assumption of a fixed realized niche relative to climate variables remains unlikely for various reasons, particularly if novel future climates open up currently unavailable portions of species’ fundamental niches. To demonstrate this effect, we compare the climate distributions for fossil‐pollen data from 21 to 15 ka bp (relying on paleoclimate simulations) when communities and climates with no modern analog were common across North America to observed modern pollen assemblages. We test how well SDMs are able to project 20th century pollen‐based taxon distributions with models calibrated using data from 21 to 15 ka. We find that taxa which were abundant in areas with no‐analog late glacial climates, such as Fraxinus, Ostrya/Carpinus and Ulmus, substantially shifted their realized niches from the late glacial period to present. SDMs for these taxa had low predictive accuracy when projected to modern climates despite demonstrating high predictive accuracy for late glacial pollen distributions. For other taxa, e.g. Quercus, Picea, Pinus strobus, had relatively stable realized niches and models for these taxa tended to have higher predictive accuracy when projected to present. Our findings reinforce the point that a realized niche at any one time often represents only a subset of the climate conditions in which a taxon can persist. Projections from SDMs into future climate conditions that are based solely on contemporary realized distributions are potentially misleading for assessing the vulnerability of species to future climate change.     

Using species distribution models to select species resistant to climate change for ecological restoration of bowé in West Africa

Bowalization is a particular form of land degradation and leads to lateral expansion of ferricrete horizons. The process occurs only in tropical regions. In this study, the most adapted and resistant species towards climate change were identified on bowé. The 15 most common bowé species of the subhumid and semi‐arid climate zones of Benin were submitted together with significant environmental variables (elevation, current bioclimatic variables, soil types) to three ecological niche modelling programmes (Maxent, Domain and GARP). For future prediction (2050), IPCC4/CIAT and IPCC5/CMIP5 climate data were applied. Asparagus africanus, Andropogon pseudapricus and Combretum nigricans were identified as the most resistant species for ecological restoration of bowé in the semi‐arid climate zone and Asparagus africanus, Detarium microcarpum and Lannea microcarpa in the subhumid climate zone. The ‘Pull’ strategies were identified as appropriate for ecological restoration of bowé in Benin.     

Wrong place,wrong time: climate change‐induced range shift across fragmented habitat causes maladaptation and declined population size in a modelled bird species

Many species are locally adapted to decreased habitat quality at their range margins, and therefore show genetic differences throughout their ranges. Under contemporary climate change, range shifts may affect evolutionary processes at the expanding range margin due to founder events. In addition, populations that are affected by such founder events will, in the course of time, become located in the range centre. Recent studies investigated evolutionary changes at the expanding range margin, but have not assessed eventual effects across the species' range. We explored the possible influence of range shift on the level of adaptation throughout the species' total range. For this we used a spatially explicit, individual‐based simulation model of a woodland bird, parameterized after the middle spotted woodpecker ( Dendrocopos medius) in fragmented habitat. We simulated its range under climate change, and incorporated genetic differences at a single locus that determined the individual's degree of adaptation to optimal temperature conditions. Generalist individuals had a large thermal tolerance, but relatively low overall fitness, whereas climate specialists had high fitness combined with a small thermal tolerance. In equilibrium, the populations in the range centre were comprised of the specialists, whereas the generalists dominated the margins. In contrast, under temperature increase, the generalist numbers increased at the expanding margin and eventually also occupied the centre of the shifting range, whereas the specialists were located in the retracting margins. This was caused by founder events and led to overall maladaptation of the species, which resulted in a reduced metapopulation size and thus impeded the species' persistence. We therefore found no evidence for a complementary effect of local adaptation and range shifts on species' survival. Instead, we showed that founder events can cause local maladaptation which can amplify throughout the species' range, and, as such, hamper the species' persistence under climate change.     

Upward elevation and northwest range shifts for alpine Meconopsis species in the Himalaya–Hengduan Mountains region

Climate change may impact the distribution of species by shifting their ranges to higher elevations or higher latitudes. The impacts on alpine plant species may be particularly profound due to a potential lack of availability of future suitable habitat. To identify how alpine species have responded to climate change during the past century as well as to predict how they may react to possible global climate change scenarios in the future, we investigate the climatic responses of seven species of Meconopsis, a representative genus endemic in the alpine meadow and subnival region of the Himalaya–Hengduan Mountains. We analyzed past elevational shifts, as well as projected shifts in longitude, latitude, elevation, and range size using historical specimen records and species distribution modeling under optimistic (RCP 4.5) and pessimistic (RCP 8.5) scenarios across three general circulation models for 2070. Our results indicate that across all seven species, there has been an upward shift in mean elevation of 302.3 m between the pre‐1970s (1922–1969) and the post‐1970s (1970–2016). The model predictions suggest that the future suitable climate space will continue to shift upwards in elevation (as well as northwards and westwards) by 2070. While for most of the analyzed species, the area of suitable climate space is predicted to expand under the optimistic emission scenario, the area contracts, or, at best, shows little change under the pessimistic scenario. Species such as M. punicea, which already occupy high latitudes, are consistently predicted to experience a contraction of suitable climate space across all the models by 2070 and may consequently deserve particular attention by conservation strategies. Collectively, our results suggest that the alpine high‐latitude species analyzed here have already been significantly impacted by climate change and that these trends may continue over the coming decades.     

Contrasting responses to Pleistocene climate changes: a case study of two sister species Allium cyathophorum and A. spicata (Amaryllidaceae) distributed in the eastern and western Qinghai–Tibet Plateau

It has been hypothesized that species occurring in the eastern and the western Qinghai–Tibet Plateau (QTP) responded differently to climate changes during the Pleistocene. Here, we test this hypothesis by phylogeographic analysis of two sister species, Allium cyathophorum and A. spicata. We sequenced two chloroplast DNA (cpDNA) fragments (accD‐psaI and the rpl16 intron) of 150 individuals, and the nuclear (ITS) region of 114 individuals, from 19 populations throughout the distributional ranges of these species. The divergence between the two species was dated at 779 ‐ 714 thousand years before the present and was likely initiated by the most major glaciation in the QTP. Analysis of chlorotype diversity showed that A. spicata, the species occurring in the western QTP, contains much lower genetic diversity (0.25) than A. cyathophorum (0.93), which is distributed in the eastern QTP. Moreover, multiple independent tests suggested that the A. spicata population had expanded recently, while no such expansion was detected in A. cyathophorum, indicating a contrasting pattern of responses to Pleistocene climate changes. These findings highlight the importance of geographical topography in determining how species responded to the climate changes that took place in the QTP during the Pleistocene.