Determinants of habitat suitability models transferability across geographically disjunct populations: Insights from Vipera ursinii urs inii

Assessing future shifts in habitat suitability and connectivity to old-growth forests to support the conservation of the endangered giant noctule

Background

Suitable climate and availability of habitats for roosting, foraging, and dispersing are critical for the long-term persistence of bat species. The giant noctule (Nyctalus lasiopterus) represents one of the lesser-known European bats, especially regarding the environmental factors which shape its distribution.

Methodology

We integrated climate-based ecological niche models with information about topography and rivers' network to model weighted suitability for N. lasiopterus in the western Palearctic. The weighted suitability map was then used to estimate connectivity among the distinct occurrence localities of N. lasiopterus, as well as from these latter towards European old-growth forests, under current conditions and different combinations of future timeframes (2030, 2050, 2070) and shared socioeconomic pathways (SSPs 3.70 and 5.85).

Results

Current weighted suitability is highest in Andalusia, northern Iberia, southwestern France, peninsular Italy, coastal Balkans and Anatolia, with dispersed suitable patches elsewhere. A north-eastward shift of weighted suitability emerges in the considered future scenarios, especially under SSP 5.85. The major current ecological corridors for N. lasiopterus are predicted within a 'belt' connecting northern Spain and southwestern France, as well as in the Italian Alps. However, following changes in weighted suitability, connectivity would increase in central-eastern Europe in the future. The bioclimatic niche of the western N. lasiopterus populations does not overlap with those of the central and eastern ones, and it only overlaps with climatic conditions characterizing old-growth forests in western Europe.

Conclusions

The outcomes of our analyses would help in designing specific conservation measures for the distinct groups of giant noctule populations, favoring the possibility of range expansion and movement towards forested habitats.

Distribution model transferability for a wide-ranging species, the Gray Wolf

Using existing data can be a reliable and cost-effective way to predict species distributions, and particularly useful for recovering or expanding species. We developed a current gray wolf (Canis lupus) distribution model for the western Great Lakes region, USA, and evaluated the spatial transferability of single-state models to the region. This study is the first assessment of transferability in a wide-ranging carnivore, as well as one of few developed for large spatial extents. We collected 3500 wolf locations from winter surveys in Minnesota (2017–2019), Wisconsin (2019–2020), and Michigan (2017–2020). We included 10 variables: proportion of natural cover, pastures, and crops; distance to natural cover, agriculture, developed land, and water; major and minor road density; and snowfall (1-km res.). We created a regional ensemble distribution by weight-averaging eight models based on their performance. We also developed single-state models, and estimated spatial transferability using two approaches: state cross-validation and extrapolation. We assessed performance by quantifying correlations, receiver operating characteristic curves (ROC), sensitivities, and two niche similarity indices. The regional area estimated to be most suitable for wolves during winter (threshold = maximum sensitivity/specificity) was 106,465 km² (MN = 48,083 km², WI = 27,757 km², MI = 30,625 km²) and correctly predicted 88% of wolf locations analyzed. Increasing natural cover and distance to crops were consistently important for determining regional and single-state wolf distribution. Extrapolation (vs. cross-validation) produced results with the greatest performance metrics, and were most similar to the regional model, yet good internal performance was unrelated to greater extrapolation performance. Factors influencing species distributions are scale-dependent and can vary across areas due to behavioral plasticity. When extending inferences beyond the current occurrence of individuals, assessing variation in ecology such as habitat selection, as well as methodological factors including model performance, will be critical to avoid poor scientific interpretations and develop effective conservation applications. In particular, accurate distribution models for recovering or recovered carnivores can be used to develop plans for habitat management, quantify potential of unoccupied habitat, assess connectivity modeling, and mitigate conflict, facilitating long-term species persistence.

Worldclim 2.1 versus Worldclim 1.4: Climatic niche and grid resolution affect between-version mismatches in Habitat Suitability Models predictions across Europe

The influence of climate on the distribution of taxa has been extensively investigated in the last two decades through Habitat Suitability Models (HSMs). In this context, the Worldclim database represents an invaluable data source as it provides worldwide climate surfaces for both historical and future time horizons. Thousands of HSMs-based papers have been published taking advantage of Worldclim 1.4, the first online version of this repository. In 2017, Worldclim 2.1 was released. Here, we evaluated spatially explicit prediction mismatch at continental scale, focusing on Europe, between HSMs fitted using climate surfaces from the two Worldclim versions (between-version differences). To this aim, we simulated occurrence probability and presence-absence across Europe of four virtual species (VS) with differing climate-occurrence relationships. For each VS, we fitted HSMs upon uncorrelated bioclimatic variables derived from each Worldclim version at three grid resolutions. For each factor combination, HSMs attaining sufficient discrimination performance on spatially independent test data were projected across Europe under current conditions and various future scenarios, and importance scores of the single variables were computed. HSMs failed in accurately retrieving the simulated climate-occurrence relationships for the climate-tolerant VS and the one occurring under a narrow combination of climatic conditions. Under current climate, noticeable between-version prediction mismatch emerged across most of Europe for these two VSs, whose simulated suitability mainly depended upon diurnal or yearly variability in temperature; differently, between-version differences were more clustered toward areas showing extreme values, like mountainous massifs or southern regions, for VSs responding to average temperature and precipitation trends. Under future climate, the chosen emission scenarios and Global Climate Models did not evidently influence between-version prediction discrepancies, while grid resolution synergistically interacted with VSs' niche characteristics in determining extent of such differences. Our findings could help in re-evaluating previous biodiversity-related works relying on geographical predictions from Worldclim-based HSMs.

Fine-Tuned Ecological Niche Models Unveil Climatic Suitability and Association with Vegetation Groups for Selected Chaetocnema Species in South Africa (Coleoptera: Chrysomelidae)

Despite beetles (Coleoptera) representing most existing animal species, the ecological and biogeographical factors shaping their distribution are still unclear in many regions. We implemented state-of-the-art ecological niche models (ENMs) and niche overlap analysis to investigate climate–occurrence patterns for five flea beetle species of the genus Chaetocnema in South Africa (C. brincki, C. danielssoni, C. darwini, C. gahani, and C. natalensis). ENMs were fitted through Maxent and Random Forests, testing various parameterizations. For each species, tuned ENMs attaining good discrimination on spatially independent test data were selected to predict suitability across the study region and individuate its main climatic drivers. Percentage coverage of climatically suitable areas by seventeen Afrotropical vegetation formations was also computed. Predicted suitable areas do not extend far away from known presence localities, except for C. brincki and C. gahani in north-eastern South Africa. Temperate grasslands and shrublands cover most of suitable areas for C. brincki and C. gahani, along with warm temperate forests, as well as for C. danielssoni, in this case being followed by tropical flooded and swamp forests. Climatic suitability for C. darwini mainly relates to the Mediterranean grasslands and scrublands of the southern coastal region, while suitable areas for C. natalensis encompass various vegetation formations, coherently with its wide distribution. The environmental niche of C. danielssoni significantly overlaps with those of the wide-ranging C. darwini and C. natalensis, suggesting that historical factors, rather than low climatic tolerance, has determined its restricted distribution in the Western Cape Province. Maxent and Random Forests were confirmed to be of great help in disentangling the environment–occurrence relationships and in predicting suitability for the target species outside their known range, but they need to be properly tuned to perform at their best.

Climate Change Increases the Expansion Risk of Helicoverpa zea in China According to Potential Geographical Distribution Estimation

Helicoverpa zea, a well-documented and endemic pest throughout most of the Americas, affecting more than 100 species of host plants. It is a quarantine pest according to the Asia and Pacific Plant Protection Commission (APPPC) and the catalog of quarantine pests for plants imported to the People's Republic of China. Based on 1781 global distribution records of H. zea and eight bioclimatic variables, the potential geographical distributions (PGDs) of H. zea were predicted by using a calibrated MaxEnt model. The contribution rate of bioclimatic variables and the jackknife method were integrated to assess the significant variables governing the PGDs. The response curves of bioclimatic variables were quantitatively determined to predict the PGDs of H. zea under climate change. The results showed that: (1) four out of the eight variables contributed the most to the model performance, namely, mean diurnal range (bio2), precipitation seasonality (bio15), precipitation of the driest quarter (bio17) and precipitation of the warmest quarter (bio18); (2) PGDs of H. zea under the current climate covered 418.15 × 104 km2, and were large in China; and (3) future climate change will facilitate the expansion of PGDs for H. zea under shared socioeconomic pathways (SSP) 1-2.6, SSP2-4.5, and SSP5-8.5 in both the 2030s and 2050s. The conversion of unsuitable to low suitability habitat and moderately to high suitability habitat increased by 8.43% and 2.35%, respectively. From the present day to the 2030s, under SSP1-2.6, SSP2-4.5 and SSP5-8.5, the centroid of the suitable habitats of H. zea showed a general tendency to move eastward; from 2030s to the 2050s, under SSP1-2.6 and SSP5-8.5, it moved southward, and it moved slightly northward under SSP2-4.5. According to bioclimatic conditions, H. zea has a high capacity for colonization by introduced individuals in China. Customs ports should pay attention to host plants and containers of H. zea and should exchange information to strengthen plant quarantine and pest monitoring, thus enhancing target management.