Large-scale environmental correlates of forest tree distributions in Catalonia (NE Spain)

Aim To explore the environmental correlates of tree species distributions in Catalonia according to the chorological status of the species. Location The study area is the region of Catalonia, in north‐eastern Spain. Methods We used presence‐absence data for 24 species, sampled in random plots distributed throughout forests of Catalonia. A climate model for the Catalonia region provided environmental variables. We used classification tree analysis to explore the environmental correlates of the realized niches of tree species. The predictive accuracy of the models was assessed using the ROC curve approach. Potential distribution maps of tree species were generated for the whole Catalonia region. Results Models were ranked from low to high accuracy for the 24 species. Differences in accuracy among species were related to the chorological status of species. Zonal species, or species at the core of their range (Mediterranean and Sub‐Mediterranean species), were generally well predicted, while extrazonal species, or species at the edge of their range, were predicted only moderately well. Mediterranean species distributions showed good correlations with extreme temperatures and annual precipitation. Main conclusions The above trends confirmed the difficulty of identifying the realized niche of species at the edges of their ranges. In contrast, Mediterranean and Sub‐Mediterranean species, which were at the core of their range, were well‐predicted, confirming the importance of extremes of temperature and annual precipitation as effective surrogates for variables having more direct physiological roles in limiting the ability of plants to survive and grow. Maps of potential tree distributions allowed us to define suitable habitats and to highlight areas where species have been planted outside their natural distribution.     

Can species distribution models be used to describe plant abundance patterns?

In recent years, there has been increasing interest in modelling of species abundance data in addition to presence data. In this study, we assessed the similarities and differences between presence‐absence distributions and abundance distributions along similar environmental gradients, derived, respectively, from presence‐absence and abundance data. Moreover, we examined the possibility of using presence‐absence distribution models to derive abundance distributions. For this purpose, we used Braun‐Blanquet abundance scores for 243 vascular species at 10 996 French forest sites. Species distribution models were used to analyse the link between the patterns of occurrence, low abundance and high abundance for each species with regard to mean annual temperature, June water balance, and soil pH. For each species, differences in the modelled distributions were characterised by the ecological optimum and ecological amplitude. A comparison of the presence‐absence and abundance distributions for all species revealed similar optima and different amplitudes along the three ecological factors. An abundant‐centre distribution was observed in environmental space, with species abundance being greatest at the optimal conditions and lower at less favourable conditions of the species occurrence response. Geographical habitat mapping also shows centred, high‐abundance suitability within the presence habitat of each species. We conclude that species distribution models derived from presence‐absence data provide useful information about the ecological optima of abundance distributions but overestimate the range of habitats suitable for high species abundance. This study demonstrates the utility of presence‐absence data for ecologist and conservation biologist when they are interested in the optimal conditions of high species abundance.     

How do steppe plants follow their optimal environmental conditions or persist under suboptimal conditions? The differing strategies of annuals and perennials

For a species to be able to respond to environmental change, it must either succeed in following its optimal environmental conditions or in persisting under suboptimal conditions, but we know very little about what controls these capacities. We parameterized species distribution models (SDMs) for 135 plant species from the Algerian steppes. We interpreted low false‐positive rates as reflecting a high capacity to follow optimal environmental conditions and high false‐negative rates as a high capacity to persist under suboptimal environmental conditions. We also measured functional traits in the field and built a unique plant trait database for the North‐African steppe. For both perennial and annual species, we explored how these two capacities can be explained by species traits and whether relevant trait values reflect species strategies or biases in SDMs. We found low false‐positive rates in species with small seeds, flowers attracting specialist pollinators, and specialized distributions (among annuals and perennials), low root:shoot ratios, wide root‐systems, and large leaves (perennials only) (R² = .52–58). We found high false‐negative rates in species with marginal environmental distribution (among annuals and perennials), small seeds, relatively deep roots, and specialized distributions (annuals) or large leaves, wide root‐systems, and monocarpic life cycle (perennials) (R² = .38 for annuals and 0.65 for perennials). Overall, relevant traits are rarely indicative of the possible biases of SDMs, but rather reflect the species' reproductive strategy, dispersal ability, stress tolerance, and pollination strategies. Our results suggest that wide undirected dispersal in annual species and efficient resource acquisition in perennial species favor both capacities, whereas short life spans in perennial species favor persistence in suboptimal environmental conditions and flowers attracting specialist pollinators in perennial and annual species favor following optimal environmental conditions. Species that neither follow nor persist will be at risk under future environmental change.     

Potential source and sink locations for climate‐driven species range shifts in Europe since the Last Glacial Maximum

Aim To identify potential source and sink locations for climate‐driven species range shifts in Europe since the Last Glacial Maximum (LGM). Location Europe. Methods We developed a new approach combining past‐climate simulations with the concept of analogous climate space. Our index gives a continuous measure of the potential of a location to have acted as a source or a sink for species that have shifted their ranges since the LGM. High glacial source potential is indicated by LGM climatic conditions that are widespread now; high post‐glacial sink potential is indicated by current climatic conditions that were widespread at the LGM. The degree of isolation of source and sink areas was calculated as the median distance to areas with analogous climate conditions. Results We identified areas of high glacial source potential in the previously recognized refugial areas in the southern European peninsulas, but also in large areas in central‐western Europe. The most climatically isolated source areas were located in northern Spain, in north‐western Europe and in eastern Turkey. From here species would have had to cover substantial distances to find current climate conditions analogous to LGM conditions of these areas. Areas with high post‐glacial sink potential were mainly located in Fennoscandia and in central and south‐eastern Europe. Some of the most isolated sink areas were located in the Spanish highlands and around the Baltic Sea. Main conclusions Our species‐independent approach successfully identified previously recognized glacial refugial areas with high source potential for species range shifts in southern Europe and in addition highlighted other potential source areas in central Europe. This study offers new insights into how the distribution of past and current climatic conditions may have influenced past species range shifts and current large‐scale biodiversity patterns.     

Predicting potential distribution of poorly known species with small database: the case of four‐horned antelope Tetracerus quadricornis on the Indian subcontinent

Information gaps on the distribution of data deficient and rare species such as four‐horned antelope (FHA) in Nepal may impair their conservation. We aimed to empirically predict the distribution of FHA in Nepal with the help of data from the Indian subcontinent. Additionally, we wanted to identify core areas and gaps within the reported range limits and to assess the degree of isolation of known Nepalese populations from the main distribution areas in India. The tropical part of the Indian subcontinent (65°–90° eastern longitude, 5°–30° northern latitude), that is, the areas south of the Himalayan Mountains. Using MaxEnt and accounting for sampling bias, we developed predictive distribution models from environmental and topographical variables, and known presence locations of the study species in India and Nepal. We address and discuss the use of target group vs. random background. The prediction map reveals a disjunct distribution of FHA with core areas in the tropical parts of central to southern–western India. At the scale of the Indian subcontinent, suitable FHA habitat area in Nepal was small. The Indo‐Gangetic Plain isolates Nepalese from the Indian FHA populations, but the distribution area extends further south than proposed by the current IUCN map. A low to intermediate temperature seasonality as well as low precipitation during the dry and warm season contributed most to the prediction of FHA distribution. The predicted distribution maps confirm other FHA range maps but also indicate that suitable areas exist south of the known range. Results further highlight that small populations in the Nepalese Terai Arc are isolated from the Indian core distribution and therefore might be under high extinction risk.     

Reconstructing range dynamics and range fragmentation of European bison for the last 8000 years

Aim Understanding what constituted species’ ranges prior to large‐scale human influence, and how past climate and land use change have affected range dynamics, provides conservation planners with important insights into how species may respond to future environmental change. Our aim here was to reconstruct the Holocene range of European bison (Bison bonasus) by combining a time‐calibrated species distribution models (SDM) with a dynamic vegetation model. Location Europe. Method We used European bison occurrences from the Holocene in a maximum entropy model to assess bison range dynamics during the last 8000 years. As predictors, we used bioclimatic variables and vegetation reconstructions from the generalized dynamic vegetation model LPJ‐GUESS. We compared our range maps with maps of farmland and human population expansion to identify the main species range constraints. Results The Holocene distribution of European bison was mainly determined by vegetation patterns, with bison thriving in both broadleaved and coniferous forests, as well as by mean winter temperature. The heartland of European bison was in Central and Eastern Europe, whereas suitable habitat in Western Europe was scarce. While environmentally suitable regions were overall stable, the expansion of settlements and farming severely diminished available habitat. Main conclusions European bison habitat preferences may be wider than previously assumed, and our results suggest that the species had a more eastern and northern distribution than previously reported. Vegetation and climate transformation during the Holocene did not affect the bison’s range substantially. Conversely, human population growth and the spread of farming resulted in drastic bison habitat loss and fragmentation, likely reaching a tipping point during the last 1000 years. Combining SDM and dynamic vegetation models can improve range reconstructions and projections, and thus help to identify resilient conservation strategies for endangered species.     

Climatic stability in the Brazilian Cerrado: implications for biogeographical connections of South American savannas,species richness and conservation in a biodiversity hotspot

Aim To investigate the historical distribution of the Cerrado across Quaternary climatic fluctuations and to generate historical stability maps to test: (1) whether the ‘historical climate’ stability hypothesis explains squamate reptile richness in the Cerrado; and (2) the hypothesis of Pleistocene connections between savannas located north and south of Amazonia. Location The Cerrado, a savanna biome and a global biodiversity hotspot distributed mainly in central Brazil. Methods We generated occurrence datasets from 1000 presence points randomly selected from the entire distribution of the Cerrado, as determined by two spatial definitions. We modelled the potential Cerrado distribution by implementing a maximum‐entropy machine‐learning algorithm across four time projections: current, mid‐Holocene (6 ka), Last Glacial Maximum (LGM, 21 ka) and Last Interglacial (LIG, 120 ka). We generated historical stability maps (refugial areas) by overlapping presence/absence projections of all scenarios, and checked consistencies with qualitative comparisons with available fossil pollen records. We built a spatially explicit simultaneous autoregressive model to explore the relationship between current climate, climatic stability, and squamate species richness. Results Models predicted the LGM and LIG as the periods of narrowest and widest Cerrado distributions, respectively, and were largely corroborated by palynological evidence. We found evidence for two savanna corridors (eastern coastal during the LIG, and Andean during the LGM) and predicted a large refugial area in the north‐eastern Cerrado (Serra Geral de Goiás refugium). Variables related to climatic stability predicted squamate richness better than present climatic variables did. Main conclusions Our results indicate that Bolivian savannas should be included within the Cerrado range and that the Cerrado’s biogeographical counterparts are not Chaco and Caatinga but rather the disjunct savannas of the Guyana shield plateaus. Climatic stability is a good predictor of Cerrado squamate richness, and our stability maps could be used in future studies to test diversity patterns and genetic signatures of different taxonomic groups and as a higher‐order landscape biodiversity surrogate for conservation planning.     

Replication levels,false presences and the estimation of the presence/absence from eDNA metabarcoding data

Environmental DNA (eDNA) metabarcoding is increasingly used to study the present and past biodiversity. eDNA analyses often rely on amplification of very small quantities or degraded DNA. To avoid missing detection of taxa that are actually present (false negatives), multiple extractions and amplifications of the same samples are often performed. However, the level of replication needed for reliable estimates of the presence/absence patterns remains an unaddressed topic. Furthermore, degraded DNA and PCR/sequencing errors might produce false positives. We used simulations and empirical data to evaluate the level of replication required for accurate detection of targeted taxa in different contexts and to assess the performance of methods used to reduce the risk of false detections. Furthermore, we evaluated whether statistical approaches developed to estimate occupancy in the presence of observational errors can successfully estimate true prevalence, detection probability and false‐positive rates. Replications reduced the rate of false negatives; the optimal level of replication was strongly dependent on the detection probability of taxa. Occupancy models successfully estimated true prevalence, detection probability and false‐positive rates, but their performance increased with the number of replicates. At least eight PCR replicates should be performed if detection probability is not high, such as in ancient DNA studies. Multiple DNA extractions from the same sample yielded consistent results; in some cases, collecting multiple samples from the same locality allowed detecting more species. The optimal level of replication for accurate species detection strongly varies among studies and could be explicitly estimated to improve the reliability of results.     

Modelling chorotypes of invasive vertebrates in mainland Spain

We investigated the existence of chorotypes – assemblages of species with similar geographical ranges – of invasive species in a host territory, and their potential use to advocate similar control or management strategies for species in the same chorotype. We analysed the distribution of 13 exotic terrestrial vertebrate species (six birds, six mammals, and one reptile) with well‐known distributions in mainland Spain. We used the presence/absence data on a grid of 10 km × 10 km UTM cells from the Atlases of terrestrial vertebrates of Spain. These data were aggregated to a grid of 50 km × 50 km UTM cells, because it entailed no loss of meaningful information and allowed dealing with a much lower number of cells. Using cluster analysis and a probabilistic assessment of the classification, we identified seven significant chorotypes: four multispecific and three monospecific. The compound chorotypes grouped together species that tended to share certain characteristics about their introduction, release cause, establishment, and spread. We modelled the chorotypes using a favourability function based on a generalized linear model and 31 variables related to spatial situation, topography, lithology, climatic stability, energy availability, water availability, disturbances, productivity, and human activity. Climatic factors affected the favourability for every chorotype, whereas human variables had a high influence in the distribution of three chorotypes involving eight species. On the basis of these variables, we identified favourable areas for all the chorotypes in mainland Spain. The favourability for a chorotype in an area may be a useful criterion for evaluating the local conservation concern due to the whole set of species. Favourable but unoccupied areas can be used to infer possible colonization areas for each chorotype. We recommend using chorotypes to optimize broad‐scale surveillance of invasive species.     

Strong genetic impoverishment from the centre of distribution in southern Europe to peripheral Baltic and isolated Scandinavian populations of the pearly heath butterfly

Aim Climatic changes and fluctuations in the past have strongly influenced the distribution of animal and plant species. Such fluctuations are also reflected in the patterns of genetic diversity on both local and global scales. The genetic pattern of the pearly heath butterfly, Coenonympha arcania, was used to evaluate the genetic differentiation of isolated (in north‐western Europe), peripheral (in north‐eastern Europe) and central (in southern Europe) populations in the context of post‐glacial distributional changes of the species. Location Europe (Sweden, Germany, the Baltic states, Italy, Slovenia, Hungary, Romania, Bulgaria). Thus, samples were collected from large parts of the species’ distribution representing the three categories mentioned above. Methods We analysed 18 loci of 569 individuals from 28 populations by allozyme electrophoresis. We used both individual‐based and population‐based analyses, including F‐statistics, various clustering methods and Markov chain Monte Carlo simulations. Results All loci, except Fum, were polymorphic. The mean F_ST for all samples was 0.18. The mean genetic distance among populations was 0.046. Two major genetic lineages were distinguished. Populations from the centre of the distributional range in southern Europe and the northern periphery of the distributional range differed significantly in their level of genetic variability. The central populations of south‐eastern Europe showed high levels of genetic diversity and no differentiation among populations. Main conclusions Most probably the two major genetic lineages evolved during glacial isolation in two disjunct Mediterranean refugia. The lack of genetic differentiation across south‐eastern Europe implies a continuous Würm ice age distribution in this area, thus supporting the functional existence of steppe forests throughout this region. The peripheral‐isolated populations in Sweden seem to have suffered from one or more severe bottlenecks, resulting in substantial genetic impoverishment. The peripheral‐connected eastern Baltic populations, on the other hand, are affected by post‐glacial and possibly recurrent gene flow from more central parts of the distribution.     

Estimation of climatic favourable areas for locust outbreaks in Spain: integrating species' presence records and spatial information on outbreaks

For millennia, locust swarms have recurrently devastated crop productivity across continents. In much of Europe, locust outbreaks have been considerably reduced by human pressure in recent decades, but important foci of outbreaks still exist in Spain. Distribution models are often used to derive spatial hypotheses and risk maps. Because insufficient information is available to include the extreme plasticity of the solitary and gregarious phases of locusts in large‐scale spatial models, we modelled the distribution either of Acrididae species or of outbreaks per se. Confirmed occurrences of Dociostaurus, Calliptamus and Chorthippus species were obtained from a field survey complemented by museum collection data and the published literature. The locations of confirmed or potential outbreaks covering two time periods of 20 years each were obtained from the literature and from Spanish autonomous community reports. Models were built with one topographic and eleven climatic predictors. We evaluated the ability of different models to predict outbreak recurrence and found that models based on Moroccan locust data or outbreak occurrence data performed the best. We generated a predictive map of the climatic favourability for locust outbreaks in Spain and found that the major foci of locust swarms were encompassed by those areas categorized by the models as areas of highest risk. Predictive maps of outbreak favourability can facilitate the more sustainable use of insecticides and more efficient integrated pest management.     

Differential success in northwards range expansion between ecotypes of the marble gallwasp Andricus kollari: a tale of two lifecycles

The Marble gallwasp Andricus kollari has a native range divided into two geographically separated lifecycles. In Eastern Europe and Turkey, the lifecycle involves a sexual generation on Turkey oak, Quercus cerris, while in Iberia and North Africa the sexual generation host is cork oak, Q. suber. Over the last 500 years, A. kollari has expanded its range into northern Europe, following human planting of Q. cerris from Italy and the Balkans. We ask: (i) what is the genetic relationship between eastern and western distributions of Andricus kollari? Can we determine which lifecycle is ancestral, and how long ago they diverged? (ii) To what extent have eastern and western native ranges contributed to northwards range expansion? (iii) Is there any evidence for hybridization between the two life cycle types? We present analyses of allozyme data for 13 polymorphic loci and of sequence variation for a 433 bp fragment of the mitochondrial cytochrome b gene. These show: (i) that four haplotype lineages (one in Spain, two in Hungary/Italy and one in Turkey) diverged more or less simultaneously between 1 and 2 million years ago, suggesting the existence of at least four refuges through recent ice age cycles. Our data cannot resolve which lifecycle type is ancestral. (ii) Populations north of putative refuges are divided into two sets. Populations in south‐west France are allied to Spain, while all remaining populations in northern Europe have been colonized from Italy and the Balkans. (iii) The transition from one race to another in south‐west France is marked by abrupt transitions in the frequency of refuge‐specific private alleles and corresponds closely to the northern limit of the distribution of cork oak. Although hybrids were detected in north‐west France, none were detected where the two lifecycles meet in south‐western France. The biology of oak gallwasps predicts that any hybrid zone will be narrow, and limited to regions where Q. cerris and Q. suber meet. Our data suggest that eastern and western A. kollari are effectively separate species.     

Spatially explicit estimation of occupancy,detection probability and survey effort needed to inform conservation planning

Aim It is increasingly recognized the importance of accounting for imperfect detection in species distribution modelling and conservation planning. However, the integration of detectability into a spatially explicit frame has received little attention. We aim (1) to show how to develop distribution maps of both detection probability and survey effort required to reliably determine a species presence/absence and (2) to increase awareness of the spatial variation of detection error inherent in studies of species occurrence. Location North‐western Spain. Methods  We registered the presence/absence of the endangered Egyptian vulture (Neophron percnopterus) in 213 surveys performed in 40 of 104 territories once known to be occupied. We model simultaneously both detection probability and occurrence, using site occupancy modelling. With the resulting regression equations, we developed distribution maps of both detection probability and required sampling effort throughout the area. Results Of the studied territories, 72.5% were detected as occupied, but after accounting for imperfect detection, the proportion of sites truly occupied was 79%. Detectability decreased in territories with higher topographical irregularity and increased with both the time of day of the survey and the progress of the season. Spatial distribution of detectability showed a mainly north–south gradient following the distribution of slope in the area. The likelihood of occupancy increased with rockier, less forested surface and less topographical irregularity within the territory. A minimum of five surveys, on average, are needed to assess, with 95% probability, the occupancy status of a site, ranging from ≤ 3 to > 24 visits/territory depending on survey‐ and site‐specific features. Main conclusions Accounting for detectability and its sources of variation allows us to elaborate distribution maps of detectability‐based survey effort. These maps are useful tools to reliably assess (e.g. with 95% probability) occupancy status throughout a landscape and provide guidance for species conservation planning.     

Atlas versus range maps: robustness of chorological relationships to distribution data types in European mammals

Aim Chorological relationships describe the patterns of distributional overlap among species. In addition to revealing biogeographical structure, the resulting clusters of species with similar geographical distributions can serve as natural units in conservation planning. Here, we assess the extent to which temporal, methodological and taxonomical differences in the source of species’ distribution data can affect the relationships that are found. Location Western Europe. Methods We used two data sets – the Atlas of European mammals and polygon range maps from the IUCN Global Mammal Assessment – both as presence–absence data for UTM 50 km × 50 km squares. We performed pairwise comparisons among 156 species for each data set to build matrices of the similarity in distribution across species, using both Jaccard’s and Baroni‐Urbani & Buser’s indices. We then compared these similarity matrices (chorological relationships), as well as the species richness and occurrence patterns from the two data sets. Results As expected, range maps increased both the mean prevalence per species and mean species richness per grid cell in comparison to atlas data, reflecting the general view that these data types respectively over‐ and underestimate species occurrence. However, species richness and occurrence patterns in atlas and range map data were positively associated and, most importantly, the chorological relationships underlying the two data sets were highly similar. Main conclusions Despite many methodological, temporal and taxonomical differences between atlas data and range maps, the chorological relationships encountered between species were similar for both data sets. Chorological analyses can thus be robust to the data source used and provide a solid basis for analytical biogeographical studies, even over broad spatial scales.     

Invasive exotic aoudad (Ammotragus lervia) as a major threat to native Iberian ibex (Capra pyrenaica): a habitat suitability model approach

The introduction of alien species to new environments is one of the main threats to the conservation of biodiversity. One particularly problematic example is that of wild ungulates which are increasingly being established in regions outside their natural distribution range due to human hunting interests. Unfortunately, we know little of the effects these large herbivores may have on the host ecosystems. This study deals with a first comparative analysis of the habitat requirements of two ungulate species that may be facing competition for resources in the south of Europe: the native Iberian ibex (Capra pyrenaica) and the exotic aoudad (Ammotragus lervia). The aoudad is a North African caprid introduced in 1970 as a game species in south‐eastern Spain. It has adapted well, and populations have been freely expanding since then. Ecological Niche Factor Analysis is used to describe the realized niche of both species where their distribution ranges merge. Both species occupy marginal areas of rugged terrain in the region. Marginality is higher for the Iberian ibex, which also presents a higher tolerance of secondary environmental gradients than the aoudad. Highly suitable areas for each species are secondarily suitable for the other. Reclassified and cross‐tabulated habitat suitability maps showing the areas of potential spatial coexistence and differences in ecological traits between both species are provided. The results obtained do not allow inferring resource competition between these species. However, current aoudad expansion could result in it invading the favoured habitats of the ibex. Inadequate hunting policy and monitoring, and increasing climatic resemblance of the study region to the native aoudad areas, due to a strong desertification process, are facilitating a high rate of expansion. We strongly recommend to eradicate or, at least, monitor these exotic populations, and promote active conservation practices, if one wants to preserve the unique natural resources present in this European region.     

Downscaling European species atlas distributions to a finer resolution: implications for conservation planning

Aim One of the limitations to using species’ distribution atlases in conservation planning is their coarse resolution relative to the needs of local planners. In this study, a simple approach to downscale original species atlas distributions to a finer resolution is outlined. If such a procedure yielded accurate downscaled predictions, then it could be an aid to using available distribution atlases in real‐world local conservation decisions. Location Europe. Methods An iterative procedure based on generalized additive modelling is used to downscale original European 50 × 50 km distributions of 2189 plant and terrestrial vertebrate species to c. 10 × 10 km grid resolution. Models are trained on 70% of the original data and evaluated on the remaining 30%, using the receiver operating characteristic (ROC) procedure. Fitted models are then interpolated to a finer resolution. A British dataset comprising distributions of 81 passerine‐bird species in a 10 × 10 km grid is used as a test bed to assess the accuracy of the downscaled predictions. European‐wide, downscaled predictions are further evaluated in terms of their ability to reproduce: (1) spatial patterns of coincidence in species richness scores among different groups; and (2) spatial patterns of coincidence in richness, rarity and complementarity hotspots. Results There was a generally good agreement between downscaled and observed fine‐resolution distributions for passerine species in Britain (median Jaccard similarity = 70%; lower quartile = 36%; upper quartile = 88%). In contrast, the correlation between downscaled and observed passerine species richness was relatively low (rho = 0.31) indicating a pattern of error propagation through the process of overlaying downscaled distributions for many species. It was also found that measures of model accuracy in fitting original data (ROC) were a poor predictor of models’ ability to interpolate distributions at fine resolutions (rho = ?0.10). Although European hotspots were not fully coincident between observed and modelled coarse‐resolution data, or between modelled coarse resolution and modelled downscaled data, there was evidence that downscaled distributions were able to maintain original cross‐taxon coincidence of species‐richness scores, at least for terrestrial vertebrate groups. Downscaled distributions were also able to uncover important environmental gradients otherwise blurred by coarse‐resolution data. Main conclusions Despite uncertainties, downscaling procedures may prove useful to identify reserves that are more meaningfully related to local patterns of environmental variation. Potential errors arising from the presence of false positives may be reduced if downscaled‐distribution records projected to occur outside the range of original coarse‐resolution data are excluded. However, the usefulness of this procedure may be limited to data‐rich regions. If downscaling procedures are applied to data‐poor regions, then there is a need to undertake further research to understand the structure of error in models. In particular, it would be important to investigate which species are poorly modelled, where and why. Without such an assessment it is difficult to support unsupervised use of downscaled data in most real‐world situations.     

How future climate and tree distribution changes shape the biodiversity of macrofungi across Europe

Aim

Climate change is affecting biodiversity at an accelerating rate. Despite the importance of fungi in ecosystems in general, and in the global carbon and nitrogen cycle in particular, there is little research on the response of fungi to climate change compared with plants and animals. Earlier studies show that climatic factors and tree species are key determinants of macrofungal diversity and distribution at large spatial scales. However, our knowledge of how climate change will affect macrofungal diversity and distribution in the future remains poorly understood.

Location

Europe.

Methods

Using openly available occurrence data of 1845 macrofungal species from eight European countries (i.e. Norway, Sweden, Finland, Denmark, Netherlands, Germany, France and Spain), we built ensemble species distribution models to predict macrofungal response to climate change alone and combined climate and tree distribution change under the IPCC special report on 2080 emissions scenarios (SRES A2 and B2).

Results

Considering climate change alone, we predict that about 77% (74.1%–80.7%) of the modelled species will expand their distribution range, and around 57% (56.1%–58.4%) of the modelled area will have an increase in macrofungal species richness. However, when considering the combined climate and tree species distribution change, only 50% (50%–50.9%) of the species are predicted to expand their distribution range and 49% (47.4%–51.1%) of the modelled area will experience an increase in macrofungal species richness.

Main Conclusions

Overall, our models projected that large areas would exhibit increased macrofungal species richness under future climate change. However, tree species distribution might play a restrictive role in the future distributional shifts of macrofungi. In addition, macrofungal responses appear heterogeneous, varying among species and regions. Our findings highlight the importance of including tree species in the projection of climate change impacts on the macrofungal diversity and distribution on a continental scale.     

Forest composition in Mediterranean mountains is projected to shift along the entire elevational gradient under climate change

Aim Species distribution models have been used frequently to assess the effects of climate change on mountain biodiversity. However, the value and accuracy of these assessments have been hampered by the use of low‐resolution data for species distributions and climatic conditions. Herein we assess potential changes in the distribution and community composition of tree species in two mountainous regions of Spain under specific scenarios of climate change using data with a high spatial resolution. We also describe potential changes in species distributions and tree communities along the entire elevational gradient. Location Two mountain ranges in southern Europe: the Central Mountain Range (central west of the Iberian Peninsula), and the Iberian Mountain Range (central east). Methods We modelled current and future distributions of 15 tree species (Eurosiberian, sub‐Mediterranean and Mediterranean species) as functions of climate, lithology and availability of soil water using generalized linear models (logistic regression) and machine learning models (gradient boosting). Using multivariate ordination of a matrix of presence/absence of tree species obtained under two Intergovernmental Panel on Climate Change (IPCC) scenarios (A2 and B2) for two different periods in the future (2041–70 and 2071–2100), we assessed the predicted changes in the composition of tree communities. Results The models predicted an upward migration of communities of Mediterranean trees to higher elevations and an associated decline in communities of temperate or cold‐adapted trees during the 21st century. It was predicted that 80–99% of the area that shows a climate suitable for cold–wet‐optimum Eurosiberian coniferous and broad‐leaved species will be lost. The largest overall changes were predicted for Mediterranean species found currently at low elevations, such as Pinus halepensis, Pinus pinaster, Quercus ilex ssp. ballota and Juniperus oxycedrus, with sharp increases in their range of 350%. Main conclusions It is likely that areas with climatic conditions suitable for cold‐adapted species will decrease significantly under climate warming. Large changes in species ranges and forest communities might occur, not only at high elevations within Mediterranean mountains but also along the entire elevational gradient throughout this region, particularly at low and mid‐elevations. Mediterranean mountains might lose their key role as refugia for cold‐adapted species and thus an important part of their genetic heritage.     

Predicting the distribution of four species of raptors (Aves: Accipitridae) in southern Spain: statistical models work better than existing maps

Aim To test the effectiveness of statistical models based on explanatory environmental variables vs. existing distribution information (maps and breeding atlas), for predicting the distribution of four species of raptors (family Accipitridae): common buzzard Buteo buteo (Linnaeus, 1758), short‐toed eagle Circaetus gallicus (Gmelin, 1788), booted eagle Hieraaetus pennatus (Gmelin, 1788) and black kite Milvus migrans (Boddaert, 1783). Location Andalusia, southern Spain. Methods Generalized linear models of 10 × 10 km squares surveyed for the presence/absence of the species by road census. Statistical models use as predictors variables derived from topography, vegetation and land‐use, and the geographical coordinates (to take account of possible spatial trends). Predictions from the models are compared with current distribution maps from the national breeding atlas and leading reference works. Results The maps derived from statistical models for all four species were more predictive than the previously published range maps and the recent national breeding atlas. The best models incorporated both topographic and vegetation and land‐use variables. Further, in three of the four species the inclusion of spatial coordinates to account for neighbourhood effects improved these models. Models for the common buzzard and black kite were highly predictive and easy to interpret from an ecological point of view, while models for short‐toed eagle and, particularly, booted eagle were not so easy to interpret, but still predicted better than previous distribution information. Main conclusions It is possible to build accurate predictive models for raptor distribution with a limited field survey using as predictors environmental variables derived from digital maps. These models integrated in a geographical information system produced distribution maps that were more accurate than previously published ones for the study species in the study area. Our study is an example of a methodology that could be used for many taxa and areas to improve unreliable distribution information.     

Climate‐driven shifts in the distribution of koala‐browse species from the Last Interglacial to the near future

The koala's Phascolarctos cinereus distribution is currently restricted to eastern and south‐eastern Australia. However, fossil records dating from 70 ± 4 ka (ka = 10³ yr) from south‐western Australia and the Nullarbor Plain are evidence of subpopulation extinctions in the southwest at least after the Last Interglacial (~128–116 ka). We hypothesize that koala sub‐population extinctions resulted from the eastward retraction of the koala's main browse species in response to unsuitable climatic conditions. We further posit a general reduction in the distribution of main koala‐browse trees in the near future in response climate change. We modelled 60 koala‐browse species and constructed a set of correlative species distribution models for five time periods: Last Interglacial (~128–116 ka), Last Glacial Maximum (~23–19 ka), Mid‐Holocene (~7–5 ka), present (interpolations of observed data, representative of 1960–1990), and 2070. We based our projections on five hindcasts and one forecast of climatic variables extracted from WorldClim based on two general circulation models (considering the most pessimistic scenario of high greenhouse‐gas emissions) and topsoil clay fraction. We used 17 dates of koala fossil specimens identified as reliable from 70 (± 4) to 535 (± 49) ka, with the last appearance of koalas at 70 ka in the southwest. The main simulated koala‐browse species were at their greatest modelled extent of suitability during the Last Glacial Maximum, with the greatest loss of koala habitat occurring between the Mid‐Holocene and the present. We predict a similar habitat loss between the present and 2070. The spatial patterns of habitat change support our hypothesis that koala extinctions in the southwest, Nullarbor Plain and central South Australia resulted from the eastward retraction of the dominant koala‐browse species in response to long‐term climate changes. Future climate patterns will likely increase the extinction risk of koalas in their remaining eastern ranges.