Modeling plant species distributions under future climates: how fine scale do climate projections need to be?

Recent studies suggest that species distribution models (SDMs) based on fine‐scale climate data may provide markedly different estimates of climate‐change impacts than coarse‐scale models. However, these studies disagree in their conclusions of how scale influences projected species distributions. In rugged terrain, coarse‐scale climate grids may not capture topographically controlled climate variation at the scale that constitutes microhabitat or refugia for some species. Although finer scale data are therefore considered to better reflect climatic conditions experienced by species, there have been few formal analyses of how modeled distributions differ with scale. We modeled distributions for 52 plant species endemic to the California Floristic Province of different life forms and range sizes under recent and future climate across a 2000‐fold range of spatial scales (0.008–16 km²). We produced unique current and future climate datasets by separately downscaling 4 km climate models to three finer resolutions based on 800, 270, and 90 m digital elevation models and deriving bioclimatic predictors from them. As climate‐data resolution became coarser, SDMs predicted larger habitat area with diminishing spatial congruence between fine‐ and coarse‐scale predictions. These trends were most pronounced at the coarsest resolutions and depended on climate scenario and species' range size. On average, SDMs projected onto 4 km climate data predicted 42% more stable habitat (the amount of spatial overlap between predicted current and future climatically suitable habitat) compared with 800 m data. We found only modest agreement between areas predicted to be stable by 90 m models generalized to 4 km grids compared with areas classified as stable based on 4 km models, suggesting that some climate refugia captured at finer scales may be missed using coarser scale data. These differences in projected locations of habitat change may have more serious implications than net habitat area when predictive maps form the basis of conservation decision making.     

The impact of climate change measured at relevant spatial scales: new hope for tropical lizards

Much attention has been given to recent predictions that widespread extinctions of tropical ectotherms, and tropical forest lizards in particular, will result from anthropogenic climate change. Most of these predictions, however, are based on environmental temperature data measured at a maximum resolution of 1 km², whereas individuals of most species experience thermal variation on a much finer scale. To address this disconnect, we combined thermal performance curves for five populations of Anolis lizard from the Bay Islands of Honduras with high‐resolution temperature distributions generated from physical models. Previous research has suggested that open‐habitat species are likely to invade forest habitat and drive forest species to extinction. We test this hypothesis, and compare the vulnerabilities of closely related, but allopatric, forest species. Our data suggest that the open‐habitat populations we studied will not invade forest habitat and may actually benefit from predicted warming for many decades. Conversely, one of the forest species we studied should experience reduced activity time as a result of warming, while two others are unlikely to experience a significant decline in performance. Our results suggest that global‐scale predictions generated using low‐resolution temperature data may overestimate the vulnerability of many tropical ectotherms to climate change.     

Simulating tick distributions over sub-Saharan Africa: the use of observed and simulated climate surfaces

Aim A broad suit of climate data sets is becoming available for use in predictive species modelling. We compare the efficacy of using interpolated climate surfaces [Center for Resource and Environmental Studies (CRES) and Climate Research Unit (CRU)] or high‐resolution model‐derived climate data [Division of Atmospheric Research limited‐area model (DARLAM)] for predictive species modelling, using tick distributions from sub‐Saharan Africa. Location The analysis is restricted to sub‐Saharan Africa. The study area was subdivided into 3000 grids cells with a resolution of 60 × 60 km. Methods Species distributions were predicted using an established multivariate climate envelope modelling approach and three very different climate data sets. The recorded variance in the climate data sets was quantified by employing omnidirectional variograms. To further compare the interpolated tick distributions that flowed from using three climate data sets, we calculated true positive (TP) predictions, false negative (FN) predictions as well as the proportional overlaps between observed and modelled tick distributions. In addition, the effect of tick data set size on the performance of the climate data sets was evaluated by performing random draws of known tick distribution records without replacement. Results The predicted distributions were consistently wider ranging than the known records when using any of the three climate data sets. However, the proportional overlap between predicted and known distributions varied as follows: for Rhipicephalus appendiculatus Neumann (Acari: Ixodidae), these were 60%, 60% and 70%; for Rhipicephalus longus Neumann (Acari: Ixodidae) 60%, 57% and 75%; for Rhipicephalus zambeziensis Walker, Norval & Corwin (Acari: Ixodidae) 57%, 51% and 62%, and for Rhipicephalus capensis Koch (Acari: Ixodidae) 70%, 60% and 60% using the CRES, CRU and DARLAM climate data sets, respectively. All data sets were sensitive to data size but DARLAM performed better when using smaller species data sets. At a 20% data subsample level, DARLAM was able to capture more than 50% of the known records and captured more than 60% of known records at higher subsample levels. Main conclusions The use of data derived from high‐resolution nested climate models (e.g. DARLAM) provided equal or even better species distribution modelling performance. As the model is dynamic and process based, the output data are available at the modelled resolution, and are not hamstrung by the sampling intensity of observed climate data sets (c. one sample per 30,000 km² for Africa). In addition, when exploring the biodiversity consequences of climate change, these modelled outputs form a more useful basis for comparison with modelled future climate scenarios.     

Habitat distribution models for intertidal seaweeds: responses to climatic and non‐climatic drivers

Aim Because intertidal organisms often live close to their physiological tolerance limits, they are potentially sensitive indicators of climate‐driven changes in the environment. The goals of this study were to assess the effect of climatic and non‐climatic factors on the geographical distribution of intertidal macroalgae, and to predict future distributions under different climate‐warming scenarios. Location North‐western Iberian Peninsula, southern Europe. Methods We developed distribution models for six ecologically important intertidal seaweed species. Occurrence and microhabitat data were sampled at 1‐km² resolution and analysed with climate variables measured at larger spatial scales. We used generalized linear models and applied the deviance and Bayesian information criterion to model the relationship between environmental variables and the distribution of each target species. We also used hierarchical partitioning (HP) to identify predictor variables with higher independent explanatory power. Results The distributions of Himanthalia elongata and Bifurcaria bifurcata were correlated with measures of terrestrial and marine climate, although in opposite directions. Model projections under two warming scenarios indicated the extinction of the former at a faster rate in the Cantabrian Sea (northern Spain) than in the Atlantic (west). In contrast, these models predicted an increase in the occurrence of B. bifurcata in both areas. The occurrences of Ascophyllum nodosum and Pelvetia canaliculata, species showing rather static historical distributions, were related to specific non‐climatic environmental conditions and locations, such as the location of sheltered sites. At the southernmost distributional limit, these habitats may present favourable microclimatic conditions or provide refuges from competitors or natural enemies. Model performances for Fucus vesiculosus and F. serratus were similar and poor, but several climatic variables influenced the occurrence of the latter in the HP analyses. Main conclusions The correlation between species distributions and climate was evident for two species, whereas the distributions of the others were associated with non‐climatic predictors. We hypothesize that the distribution of F. serratus responds to diverse combinations of factors in different sections of the north‐west Iberian Peninsula. Our study shows how the response of species distributions to climatic and non‐climatic variables may be complex and vary geographically. Our analyses also highlight the difficulty of making predictions based solely on variation in climatic factors measured at coarse spatial scales.     

The contributions of topoclimate and land cover to species distributions and abundance: fine‐resolution tests for a mountain butterfly fauna

Aim Models relating species distributions to climate or habitat are widely used to predict the effects of global change on biodiversity. Most such approaches assume that climate governs coarse‐scale species ranges, whereas habitat limits fine‐scale distributions. We tested the influence of topoclimate and land cover on butterfly distributions and abundance in a mountain range, where climate may vary as markedly at a fine scale as land cover. Location Sierra de Guadarrama (Spain, southern Europe) Methods We sampled the butterfly fauna of 180 locations (89 in 2004, 91 in 2005) in a 10,800 km² region, and derived generalized linear models (GLMs) for species occurrence and abundance based on topoclimatic (elevation and insolation) or habitat (land cover, geology and hydrology) variables sampled at 100‐m resolution using GIS. Models for each year were tested against independent data from the alternate year, using the area under the receiver operating characteristic curve (AUC) (distribution) or Spearman's rank correlation coefficient (r_s) (abundance). Results In independent model tests, 74% of occurrence models achieved AUCs of > 0.7, and 85% of abundance models were significantly related to observed abundance. Topoclimatic models outperformed models based purely on land cover in 72% of occurrence models and 66% of abundance models. Including both types of variables often explained most variation in model calibration, but did not significantly improve model cross‐validation relative to topoclimatic models. Hierarchical partitioning analysis confirmed the overriding effect of topoclimatic factors on species distributions, with the exception of several species for which the importance of land cover was confirmed. Main conclusions Topoclimatic factors may dominate fine‐resolution species distributions in mountain ranges where climate conditions vary markedly over short distances and large areas of natural habitat remain. Climate change is likely to be a key driver of species distributions in such systems and could have important effects on biodiversity. However, continued habitat protection may be vital to facilitate range shifts in response to climate change.     

Anticipated climate and land‐cover changes reveal refuge areas for Borneo's orang‐utans

Habitat loss and climate change pose a double jeopardy for many threatened taxa, making the identification of optimal habitat for the future a conservation priority. Using a case study of the endangered Bornean orang‐utan, we identify environmental refuges by integrating bioclimatic models with projected deforestation and oil‐palm agriculture suitability from the 1950s to 2080s. We coupled a maximum entropy algorithm with information on habitat needs to predict suitable habitat for the present day and 1950s. We then projected to the 2020s, 2050s and 2080s in models incorporating only land‐cover change, climate change or both processes combined. For future climate, we incorporated projections from four model and emission scenario combinations. For future land cover, we developed spatial deforestation predictions from 10 years of satellite data. Refuges were delineated as suitable forested habitats identified by all models that were also unsuitable for oil palm – a major threat to tropical biodiversity. Our analyses indicate that in 2010 up to 260 000 km² of Borneo was suitable habitat within the core orang‐utan range; an 18–24% reduction since the 1950s. Land‐cover models predicted further decline of 15–30% by the 2080s. Although habitat extent under future climate conditions varied among projections, there was majority consensus, particularly in north‐eastern and western regions. Across projections habitat loss due to climate change alone averaged 63% by 2080, but 74% when also considering land‐cover change. Refuge areas amounted to 2000–42 000 km² depending on thresholds used, with 900–17 000 km² outside the current species range. We demonstrate that efforts to halt deforestation could mediate some orang‐utan habitat loss, but further decline of the most suitable areas is to be expected given projected changes to climate. Protected refuge areas could therefore become increasingly important for ongoing translocation efforts. We present an approach to help identify such areas for highly threatened species given environmental changes expected this century.     

Potential impacts of climate change on geographical distribution of three primary vectors of African Trypanosomiasis in Tanzania’s Maasai Steppe: G. m. morsitans,G. pallidipes and G. swynnertoni

In the Maasai Steppe, public health and economy are threatened by African Trypanosomiasis, a debilitating and fatal disease to livestock (African Animal Trypanosomiasis -AAT) and humans (Human African Trypanosomiasis—HAT), if not treated. The tsetse fly is the primary vector for both HAT and AAT and climate is an important predictor of their occurrence and the parasites they carry. While understanding tsetse fly distribution is essential for informing vector and disease control strategies, existing distribution maps are old and were based on coarse spatial resolution data, consequently, inaccurately representing vector and disease dynamics necessary to design and implement fit-for-purpose mitigation strategies. Also, the assertion that climate change is altering tsetse fly distribution in Tanzania lacks empirical evidence. Despite tsetse flies posing public health risks and economic hardship, no study has modelled their distributions at a scale needed for local planning. This study used MaxEnt species distribution modelling (SDM) and ecological niche modeling tools to predict potential distribution of three tsetse fly species in Tanzania’s Maasai Steppe from current climate information, and project their distributions to midcentury climatic conditions under representative concentration pathways (RCP) 4.5 scenarios. Current climate results predicted that G. m. morsitans, G. pallidipes and G swynnertoni cover 19,225 km², 7,113 km² and 32,335 km² and future prediction indicated that by the year 2050, the habitable area may decrease by up to 23.13%, 12.9% and 22.8% of current habitable area, respectively. This information can serve as a useful predictor of potential HAT and AAT hotspots and inform surveillance strategies. Distribution maps generated by this study can be useful in guiding tsetse fly control managers, and health, livestock and wildlife officers when setting surveys and surveillance programs. The maps can also inform protected area managers of potential encroachment into the protected areas (PAs) due to shrinkage of tsetse fly habitats outside PAs.     

Spatial scale,topography and thermoregulatory behaviour interact when modelling species’ thermal niches

The spatial scale at which climate and species’ occupancy data are gathered, and the resolution at which ecological models are run, can strongly influence predictions of species performance and distributions. Running model simulations at coarse rather than fine spatial resolutions, for example, can determine if a model accurately predicts the distribution of a species. The impacts of spatial scale on a model's accuracy are particularly pronounced across mountainous terrain. Understanding how these discrepancies arise requires a modelling approach in which the underlying processes that determine a species’ distribution are explicitly described. Here we use a process‐based model to explore how spatial resolution, topography and behaviour alter predictions of a species thermal niche, which in turn constrains its survival and geographic distribution. The model incorporates biophysical equations to predict the operative temperature (T_e), thermal‐dependent performance and survival of a typical insect, with a complex life‐cycle, in its microclimate. We run this model with geographic data from a mountainous terrain in South Africa using climate data at three spatial resolutions. We also explore how behavioural thermoregulation affects predictions of a species performance and survival by allowing the animal to select the optimum thermal location within each coarse‐grid cell. At the regional level, coarse‐resolution models predicted lower T_e at low elevations and higher T_e at high elevations than models run at fine‐resolutions. These differences were more prominent on steep, north‐facing slopes. The discrepancies in T_e in turn affected estimates of the species thermal niche. The modelling framework revealed how spatial resolution and topography influence predictions of species distribution models, including the potential impacts of climate change. These systematic biases must be accounted for when interpreting the outputs of future modelling studies, particularly when species distributions are predicted to shift from uniform to topographically heterogeneous landscapes.     

A consistent occupancy–climate relationship across birds and mammals of the Americas

At broad spatial scales, species richness is strongly related to climate. Yet, few ecological studies attempt to identify regularities in the individual species distributions that make up this pattern. Models used to describe species distributions typically model very complex responses to climate. Here, we test whether the variability in the distributions of birds and mammals of the Americas relates to mean annual temperature and precipitation in a simple, consistent way. Specifically, we test if simple mathematical models can predict, as a first approximation, the geographical variation in individual species’ probability of occupancy for 3277 non‐migratory bird and 1659 mammal species. We find a Gaussian model, where the probability of occupancy of a 10⁴ km² quadrat decreases symmetrically and gradually around a species ‘optimal’ temperature and precipitation, was generally the best model, explaining an average of 35% of the deviance in probability of occupancy. The inclusion of additional terms had very small and idiosyncratic effects across species. The Gaussian occupancy–climate relationship appears general among species and taxa and explains nearly as much deviance as complex models including many more parameters. Therefore, we propose that hypotheses aiming to explain the broad‐scale distribution of species or species richness must also predict generally Gaussian occupancy–climate relationships. Synthesis Science aims to identify regularities in a complex natural world. General patterns should be identified before one searches for potential mechanisms and contingencies. However, species geographic distributions are often modelled as complex (sometimes black box), species‐specific, functions of their environment. We asked whether a simple model could account for as much of the geographic variation in a species' probability of occupancy, and be widely applicable across thousands of species. As a first approximation, we found that a simple Gaussian occupancy‐climate relationship is very common in Nature, whether it be causal or not.     

Potential distribution patterns of scorpions in north‐eastern Brazil under scenarios of future climate change

Current predictions about the responses of species to climate change strongly rely on projecting altered environmental conditions on their distributions. In this study, we investigated the effects of future climate change scenarios on the potential distribution of 10 species of scorpions in north‐eastern Brazil in the context of their degree of specialisation to closed (Atlantic and Amazon Forests) and open (Caatinga and Cerrado) habitats. Scorpion species were classified as habitat specialists or generalists according to the IndVal index, and present and future species distribution models were prepared using minimum volume ellipsoids. According to IndVal, four species were classified as closed‐forest specialists (Ananteris mauryi, Tityus brazilae, Tityus pusillus and Tityus neglectus), four as open‐forest specialists (Jaguajir agamemnon, Jaguajir rochae, Physoctonus debilis and Bothriurus rochai), and two as generalists (Tityus stigmurus and Bothriurus asper). All species presented a drastic reduction in potential distribution, ranging from 44% to 72%, when compared with their current distribution. In addition, we found a reduction in scorpion species richness under future climate change scenarios. This finding has implications for scorpion conservation. Further, the results show that climate change may impact the composition of scorpion assemblages in north‐eastern Brazil, revealing important implications for human–scorpion interactions.     

Inclusion of local environmental conditions alters high-latitude vegetation change predictions based on bioclimatic models

Current predictions of how species will respond to climate change are typically based on coarse-grained climate surfaces utilizing bioclimate envelope modelling. However, the suitability of environmental conditions for a given species might result from a variety of factors including some unrelated to climate. To address this issue, we investigated whether the inclusion of topographical and soil information in bioclimatic envelope models would significantly alter predictions of climate change—induced fine-scale tree and shrub species range size changes at the tree-limit in subarctic Europe. Using generalized additive models and data on current climate and species distributions and three different climate scenarios for the period 2040–2069, we developed predictions of the currently suitable area and potential range size changes of seven tree and shrub species in an area of 1,100 km² at a resolution of 1-ha. The inclusion of topography and soil information increased the predictive accuracy of climate-only models for all studied species. The predicted changes in species distribution volumes were contradictory, and the predicted occurrences varied greatly depending on the model used. Our results therefore support the arguments that vegetation responses to climate change can be influenced by local environmental conditions and that attention should be paid to the combined effects of these factors. We conclude that disregarding local topography and soil conditions in bioclimatic models may result in biased projections of range expansions and the associated colonization, extinction and turnover assessments.     

Multi‐model comparison highlights consistency in predicted effect of warming on a semi‐arid shrub

A number of modeling approaches have been developed to predict the impacts of climate change on species distributions, performance, and abundance. The stronger the agreement from models that represent different processes and are based on distinct and independent sources of information, the greater the confidence we can have in their predictions. Evaluating the level of confidence is particularly important when predictions are used to guide conservation or restoration decisions. We used a multi‐model approach to predict climate change impacts on big sagebrush (Artemisia tridentata), the dominant plant species on roughly 43 million hectares in the western United States and a key resource for many endemic wildlife species. To evaluate the climate sensitivity of A. tridentata, we developed four predictive models, two based on empirically derived spatial and temporal relationships, and two that applied mechanistic approaches to simulate sagebrush recruitment and growth. This approach enabled us to produce an aggregate index of climate change vulnerability and uncertainty based on the level of agreement between models. Despite large differences in model structure, predictions of sagebrush response to climate change were largely consistent. Performance, as measured by change in cover, growth, or recruitment, was predicted to decrease at the warmest sites, but increase throughout the cooler portions of sagebrush's range. A sensitivity analysis indicated that sagebrush performance responds more strongly to changes in temperature than precipitation. Most of the uncertainty in model predictions reflected variation among the ecological models, raising questions about the reliability of forecasts based on a single modeling approach. Our results highlight the value of a multi‐model approach in forecasting climate change impacts and uncertainties and should help land managers to maximize the value of conservation investments.     

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.     

Species interactions weakly modify climate‐induced tree co‐occurrence patterns

Aims

Species distributions are hypothesized to be underlain by a complex association of processes that span multiple spatial scales including biotic interactions, dispersal limitation, fine‐scale resource gradients and climate. Species disequilibrium with climate may reflect the effects of non‐climatic processes on species distributions, yet distribution models have rarely directly considered non‐climatic processes. Here, we use a Joint Species Distribution Model (JSDM) to investigate the influence of non‐climatic factors on species co‐occurrence patterns and to directly quantify the relative influences of climate and alternative processes that may generate correlated responses in species distributions, such as species interactions, on tree co‐occurrence patterns.

Location

US Rocky Mountains.

Methods

We apply a Bayesian JSDM to simultaneously model the co‐occurrence patterns of ten dominant tree species across the Rocky Mountains, and evaluate climatic and residual correlations from the fitted model to determine the relative contribution of each component to observed co‐occurrence patterns. We also evaluate predictions generated from the fitted model relative to a single‐species modelling approach.

Results

For most species, correlation due to climate covariates exceeded residual correlation, indicating an overriding influence of broad‐scale climate on co‐occurrence patterns. Accounting for covariance among species did not significantly improve predictions relative to a single‐species approach, providing limited evidence for a strong independent influence of species interactions on distribution patterns.

Conclusions

Overall, our findings indicate that climate is an important driver of regional biodiversity patterns and that interactions between dominant tree species contribute little to explain species co‐occurrence patterns among Rocky Mountain trees.     

Ecological and human impacts on stand density and distribution of tamarind (Tamarindus indica L.) in Senegal

Indigenous fruit tree species such as tamarind (Tamarindus indica L.) in African sub‐Saharan traditionally act to build resilience into the farming system in terms of food security, income generation and ecosystem stability. Therefore, increasing our knowledge on their ecology and distribution is a priority. Tamarind is mainly grown for the fruits but is also a valuable timber species. The fruit pulp has a high content of vitamin B and is eaten fresh or made into jam, chutney, juice or sweets. Flowers, leaves and seeds are also edible and used in a variety of dishes. The main objective of this study is to evaluate actual density of tamarind in Senegal and the climate change effects on its distribution for better conservation strategies. Tamarind's distribution and density around villages were recorded and modelled in different agro‐ecological zones in Senegal using transect method and under current and future climates. Distribution under two future climate scenarios were modelled using four climate models and three time slices (2020, 2050 and 2080). Results show a decreasing gradient in tree density (from 7 to 1 trees km^?2) from the Sudano agro‐ecological zone (in the south) to the Sahel (in the north). Future climate predictions show that although tamarind distribution will increase in the north‐west and south of the country in 2020; by 2050, the area identified as suitable for its growth will be greatly reduced. Areas in the north‐west basin appear to be an important refugia for the species under future climate conditions. However, density around villages in this area was found to be relatively low indicating that this could lead to problems of poor connectivity and inbreeding depression. This region should therefore be highlighted as important conservative management and protection strategies of tamarind in this region.     

Spatial conservation planning under climate change: Using species distribution modeling to assess priority for adaptive management of Fagus crenata in Japan

Protected areas are the basis of modern conservation systems, but current climate change causes gaps between protected areas and the species distribution ranges. To mitigate the impact of climate change on species distribution ranges, revision of protected areas are necessary. Alternatively, active management such as excluding competitive species or transplanting target species would be effective. In this study, we assessed optimal actions (revision of protected areas or active management) in each geographical region to establish an effective spatial conservation plan in Japan. Gaps between the protected areas and future potential habitats were assessed using species distribution models and 20 future climate simulations. Fagus crenata, an endemic and dominant species in Japan, was used as a target species. Potential habitats within the protected areas were predicted to decrease from 22,122 km² at present to 12,309 km² under future climate conditions. Sustainable potential habitats (consistent potential habitats both at present and in future) without the protected areas extended to 13,208 km², and were mainly found in northeast Japan. These results suggest that, in northeast Japan, revisions to protected areas would be effective in preserving sustainable potential habitats under future climate change. However, the potential habitats of southwestern Japan, in which populations were genetically different from northeastern populations, were predicted to virtually disappear both within and outside of protected areas. Active management is thus necessary in southwestern Japan to ensure intraspecific genetic diversity under future climate change.     

Geological habitat template overrides late Quaternary climate change as a determinant of range dynamics and phylogeography in some habitat‐specialist water beetles

Aim We investigated the roles of lithology and climate in constraining the ranges of four co‐distributed species of Iberian saline‐habitat specialist water beetles (Ochthebius glaber, Ochthebius notabilis, Enochrus falcarius and Nebrioporus baeticus) across the late Quaternary and in shaping their geographical genetic structure. The aim was to improve our understanding of the effects of past climate changes on the biota of arid Mediterranean environments and of the relative importance of history and landscape on phylogeographical patterns. Location Iberian Peninsula, Mediterranean. Methods We combined species distribution modelling (SDM) and comparative phylogeography. We used a multi‐model inference and model‐averaging approach both for assessment of range determinants (climate and lithology) and for provision of spatially explicit estimates of the species current and Last Glacial Maximum (LGM) potential ranges. Potential LGM distributions were then contrasted with the phylogeographical and population expansion patterns as assessed using mitochondrial DNA sequence data. We also evaluated the relative importance of geographical distance, habitat resistance and historical isolation for genetic structure in a causal modelling framework. Results Lithology poses a strong constraint on the distribution of Iberian saline‐habitat specialist water beetles, with a variable, but generally moderate, additional influence by climate. The degree to which potential LGM distributions were reduced and fragmented decreased with increasing importance of lithology. These SDM‐based suitability predictions were mostly congruent with phylogeographical and population genetic patterns across the study species, with stronger geographical structure in the genetic diversity of the more temperature‐sensitive species (O. glaber and E. falcarius). Furthermore, while historical isolation was the only factor explaining genetic structure in the more temperature‐sensitive species, lithology‐controlled landscape configuration also played an important role for those species with more lithology‐determined ranges (O. notabilis and N. baeticus). Main conclusions Our data show that lithology is an important constraint on the distribution and range dynamics of endemic Iberian saline‐habitat water beetles, in interaction with climate and long‐term climate change, and overrides the latter in importance for some species. Hence, geological landscape structure and long‐term history may codetermine the overall range and the distribution of genetic lineages in endemic species with specialized edaphic requirements.     

Predicting the potential distribution of a riparian invasive plant: the effects of changing climate,flood regimes and land‐use patterns

Climate change is likely to affect plants in multiple ways, but predicting the consequences for habitat suitability requires a process‐based understanding of the interactions. This is at odds with existing approaches that are mostly phenomenological and largely restricted to predicting the effects of changing temperature and rainfall on species distributions at a coarse spatial scale. We examine the multiple effects of climate change, including predicting the effects of altered flood regimes and land‐use change, on the potential distribution of the invasive riparian species lippia (Phyla canescens) across a 26 000 km² catchment in eastern Australia. We determined habitat suitability for lippia by combining process‐understanding of experts and an eco‐physiological bioclimatic model within a Bayesian belief network. The bioclimatic model predicted substantial changes in habitat suitability by 2070 under both a wetter (Echam Mark 3) and drier (Hadley Centre Mark 2) climate change scenario, but only the more likely drier scenario reduced suitability in our test region. The area suitable for lippia was predicted to increase at least threefold with increased flooding under a wet climate scenario, although this would be partially negated by land‐use change to cultivation. The region would become unsuitable to lippia with reduced flooding under a drier scenario irrespective of land‐use changes, although existing populations would persist if grazing persisted. Independent field validation verified model structure and parameterization, and therefore the opinion of experts, but identified site‐scale deficiencies in the available environmental data layers. Model predictions suggest that adaptation options for managing lippia will be greatly reduced under a drying scenario, but identify potential restoration opportunities under either scenario. This work highlights the value of predictive models that incorporate process‐understanding at sufficiently fine spatial resolution to capture the important processes underpinning habitat suitability.     

Whole range and regional‐based ecological niche models predict differing exposure to 21st century climate change in the key cool temperate rainforest tree southern beech (Nothofagus cunninghamii)

The warmer and drier climates projected for the mid‐ to late‐21st century may have particularly adverse impacts on the cool temperate rainforests of southeastern Australia. Southern beech (Nothofagus cunninghamii; Nothofagaceae), a dominant tree species in these forests, may be vulnerable to minor changes in its climate envelope, especially at the edge of the species range, with Holocene fossil evidence showing local extinction of populations in response to small climate changes. We modelled the stability of this species climate envelope using the maximum entropy algorithm implemented in Maxent and two thresholds of presence/absence by projecting the modern climate envelope onto four Global Circulation Models forecasted for two time periods (2050s and 2070s). The climate envelope, as estimated from the species present climatic range, is predicted to shrink by up to 49% by the 2050s and up to 64% by the 2070s. The greatest predicted reduction is in Victoria with 91–100% of its current range being climatically unsuitable by the 2070s. Climatically similar areas to the species present range are predicted to remain in mountainous areas of western Tasmania, the Northeast Highlands of Tasmania, and the Baw Baw Plateau in the Central Highlands of Victoria. However, region‐specific modelling approaches made very different predictions from the whole‐range based models, especially in the severity of the predicted decline for Victorian populations of N. cunninghamii which occur in much warmer climates than the rest of the species geographical range. This shows that, for widespread species that span a range of climate zones, the exposure of current populations to climate change may be better modelled using a regional based approach. How the species responds to climate change will depend on the species ability to respond to drier and warmer climates and the concomitant increase in fire intensity.     

Scale dependency of metapopulation models used to predict climate change impacts on small mammals

To investigate potential range shifts in a changing climate it is becoming increasingly common to develop models that account for demographic processes. Metapopulation models incorporate the spatial configuration of occupied habitat (i.e. arrangement, size and quality), population demographics, and inter‐patch dispersal making them suitable for investigating potential threats to small mammal range and abundance. However, the spatial scale (resolution) used to represent species–environment dynamics may affect estimates of range shift and population resilience by failing to realistically represent the spatial configuration of suitable habitat, including stepping stones and refugia. We aimed to determine whether relatively fine‐scale environmental information influenced predictions of metapopulation persistence and range shift. Species distribution models were constructed for four small terrestrial mammals from southern Australia using environmental predictors measured at 0.1 × 0.1 km (0.01 km²) or 1.0 × 1.0 km (1 km²) resolution, and combined with demographic information to parameterise coupled niche‐population models. These models were used to simulate population dynamics projected over 40‐yr under a stable and changing climate. Initial estimates of the area of available habitat were similar at both spatial scales. However, at the fine‐scale, habitat configuration comprised a greater number of patches (ca 12 times), that were more irregular in shape (ca 8 times the perimeter:area), and separated by a tenth of the distance than at the coarse‐scale. While small patches were not more prone to extinction, populations generally declined at a higher rate and were associated with a lower expected minimum abundance. Despite increased species vulnerability at the fine‐scale, greater range shifts were measured at the coarse‐scale (for species illustrating a shift at both scales). These results highlight the potential for range shifts and species vulnerability information to be misrepresented if advanced modelling techniques incorporating species demographics and dispersal inadequately represent the scale at which these processes occur.