Bio‐ORACLE: a global environmental dataset for marine species distribution modelling

Aim The oceans harbour a great diversity of organisms whose distribution and ecological preferences are often poorly understood. Species distribution modelling (SDM) could improve our knowledge and inform marine ecosystem management and conservation. Although marine environmental data are available from various sources, there are currently no user‐friendly, high‐resolution global datasets designed for SDM applications. This study aims to fill this gap by assembling a comprehensive, uniform, high‐resolution and readily usable package of global environmental rasters. Location Global, marine. Methods We compiled global coverage data, e.g. satellite‐based and in situ measured data, representing various aspects of the marine environment relevant for species distributions. Rasters were assembled at a resolution of 5 arcmin (c. 9.2 km) and a uniform landmask was applied. The utility of the dataset was evaluated by maximum entropy SDM of the invasive seaweed Codium fragile ssp. fragile. Results We present Bio‐ORACLE (ocean rasters for analysis of climate and environment), a global dataset consisting of 23 geophysical, biotic and climate rasters. This user‐friendly data package for marine species distribution modelling is available for download at http://www.bio‐oracle.ugent.be . The high predictive power of the distribution model of C. fragile ssp. fragile clearly illustrates the potential of the data package for SDM of shallow‐water marine organisms. Main conclusions The availability of this global environmental data package has the potential to stimulate marine SDM. The high predictive success of the presence‐only model of a notorious invasive seaweed shows that the information contained in Bio‐ORACLE can be informative about marine distributions and permits building highly accurate species distribution models.     

Assessing the effectiveness of foraging radius models for seabird distributions using biotelemetry and survey data

Relatively simple foraging radius models have the potential to generate predictive distributions for a large number of species rapidly, thus providing a cost-effective alternative to large-scale surveys or complex modelling approaches. Their effectiveness, however, remains largely untested. Here we compare foraging radius distribution models for all breeding seabirds in Ireland, to distributions of empirical data collected from tracking studies and aerial surveys. At the local/colony level, we compared foraging radius distributions to GPS tracking data from seabirds with short (Atlantic puffin Fratercula arctica, and razorbill Alca torda) and long (Manx shearwater Puffinus puffinus, and European storm-petrel Hydrobates pelagicus) foraging ranges. At the regional/national level, we compared foraging radius distributions to extensive aerial surveys conducted over a two-year period. Foraging radius distributions were significantly positively correlated with tracking data for all species except Manx shearwater. Correlations between foraging radius distributions and aerial survey data were also significant, but generally weaker than those for tracking data. Correlations between foraging radius distributions and aerial survey data were benchmarked against generalised additive models (GAMs) of the aerial survey data that included a range of environmental covariates. While GAM distributions had slightly higher correlations with aerial survey data, the results highlight that the foraging radius approach can be a useful and pragmatic approach for assessing breeding distributions for many seabird species. The approach is likely to have acceptable utility in complex, temporally variable ecosystems and when logistic and financial resources are limited.     

Detecting outliers in species distribution data

Aim

Species distribution data play a pivotal role in the study of ecology, evolution, biogeography and biodiversity conservation. Although large amounts of location data are available and accessible from public databases, data quality remains problematic. Of the potential sources of error, positional errors are critical for spatial applications, particularly where these errors place observations beyond the environmental or geographical range of species. These outliers need to be identified, checked and removed to improve data quality and minimize the impact on subsequent analyses. Manually checking all species records within large multispecies datasets is prohibitively costly. This work investigates algorithms that may assist in the efficient vetting of outliers in such large datasets.

Location

We used real, spatially explicit environmental data derived from the western part of Victoria, Australia, and simulated species distributions within this same region.

Methods

By adapting species distribution modelling (SDM), we developed a pseudo‐SDM approach for detecting outliers in species distribution data, which was implemented with random forest (RF) and support vector machine (SVM) resulting in two new methods: RF_pdSDM and SVM_pdSDM. Using virtual species, we compared eight existing multivariate outlier detection methods with these two new methods under various conditions.

Results

The two new methods based on the pseudo‐SDM approach had higher true skill statistic (TSS) values than other approaches, with TSS values always exceeding 0. More than 70% of the true outliers in datasets for species with a low and intermediate prevalence can be identified by checking 10% of the data points with the highest outlier scores.

Main conclusions

Pseudo‐SDM‐based methods were more effective than other outlier detection methods. However, this outlier detection procedure can only be considered as a screening tool, and putative outliers must be examined by experts to determine whether they are actual errors or important records within an inherently biased set of data.     

A review of the occurrence of inter‐colony segregation of seabird foraging areas and the implications for marine environmental impact assessment

Understanding the determinants of species’ distributions is a fundamental aim in ecology and a prerequisite for conservation but is particularly challenging in the marine environment. Advances in bio‐logging technology have resulted in a rapid increase in studies of seabird movement and distribution in recent years. Multi‐colony studies examining the effects of intra‐ and inter‐colony competition on distribution have found that several species exhibit inter‐colony segregation of foraging areas, rather than overlapping distributions. These findings are timely given the increasing rate of human exploitation of marine resources and the need to make robust assessments of likely impacts of proposed marine developments on biodiversity. Here we review the occurrence of foraging area segregation reported by published tracking studies in relation to the density‐dependent hinterland (DDH) model, which predicts that segregation occurs in response to inter‐colony competition, itself a function of colony size, distance from the colony and prey distribution. We found that inter‐colony foraging area segregation occurred in 79% of 39 studies. The frequency of occurrence was similar across the four seabird orders for which data were available, and included species with both smaller (10–100 km) and larger (100–1000 km) foraging ranges. Many predictions of the DDH model were confirmed, with examples of segregation in response to high levels of inter‐colony competition related to colony size and proximity, and enclosed landform restricting the extent of available habitat. Moreover, as predicted by the DDH model, inter‐colony overlap tended to occur where birds aggregated in highly productive areas, often remote from all colonies. The apparent prevalence of inter‐colony foraging segregation has important implications for assessment of impacts of marine development on protected seabird colonies. If a development area is accessible from multiple colonies, it may impact those colonies much more asymmetrically than previously supposed. Current impact assessment approaches that do not consider spatial inter‐colony segregation will therefore be subject to error. We recommend the collection of tracking data from multiple colonies and modelling of inter‐colony interactions to predict colony‐specific distributions.     

Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Mapping from heterogeneous biodiversity monitoring data sources

Field monitoring can vary from simple volunteer opportunistic observations to professional standardised monitoring surveys, leading to a trade-off between data quality and data collection costs. Such variability in data quality may result in biased predictions obtained from species distribution models (SDMs). We aimed to identify the limitations of different monitoring data sources for developing species distribution maps and to evaluate their potential for spatial data integration in a conservation context. Using Maxent, SDMs were generated from three different bird data sources in Catalonia, which differ in the degree of standardisation and available sample size. In addition, an alternative approach for modelling species distributions was applied, which combined the three data sources at a large spatial scale, but then downscaling to the required resolution. Finally, SDM predictions were used to identify species richness and high quality areas (hotspots) from different treatments. Models were evaluated by using high quality Atlas information. We show that both sample size and survey methodology used to collect the data are important in delivering robust information on species distributions. Models based on standardized monitoring provided higher accuracy with a lower sample size, especially when modelling common species. Accuracy of models from opportunistic observations substantially increased when modelling uncommon species, giving similar accuracy to a more standardized survey. Although downscaling data through a SDM approach appears to be a useful tool in cases of data shortage or low data quality and heterogeneity, it will tend to overestimate species distributions. In order to identify distributions of species, data with different quality may be appropriate. However, to identify biodiversity hotspots high quality information is needed.     

Fine scale waterbody data improve prediction of waterbird occurrence despite coarse species data

While modelling habitat suitability and species distribution, ecologists must deal with issues related to the spatial resolution of species occurrence and environmental data. Indeed, given that the spatial resolution of species and environmental datasets range from centimeters to hundreds of kilometers, it underlines the importance of choosing the optimal combination of resolutions to achieve the highest possible modelling prediction accuracy. We evaluated how the spatial resolution of land cover/waterbody datasets (meters to 1 km) affect waterbird habitat suitability models based on atlas data (grid cell of 12 × 11 km). We hypothesized that the area, perimeter and number of waterbodies computed from high resolution datasets would explain distributions of waterbirds better because coarse resolution datasets omit small waterbodies affecting species occurrence. Specifically, we investigated which spatial resolution of waterbodies better explain the distribution of seven waterbirds nesting on ponds/lakes with areas ranging from 0.1 ha to hundreds of hectares. Our results show that the area and perimeter of waterbodies derived from high resolution datasets (raster data with 30 m resolution, vector data corresponding with map scale 1:10 000) explain the distribution of the waterbirds better than those calculated using less accurate datasets despite the coarse grain of the species data. Taking into account the spatial extent (global vs regional) of the datasets, we found the Global Inland Waterbody Dataset to be the most suitable for modelling distribution of waterbirds. In general, we recommend using land cover data of a resolution sufficient to capture the smallest patches of the habitat suitable for a given species’ presence for both fine and coarse grain habitat suitability and distribution modelling.     

Long-term monitoring data meet freshwater species distribution models: Lessons from an LTER-site

Long-term monitoring datasets provide a solid framework for ecological research. Such a dataset from the German long-term ecological research (LTER) site Rhine-Main-Observatory was used to set up a species distribution model (SDM) for the Kinzig catchment. The extensive knowledge on the monitoring data provided by the LTER-site framework allowed to calibrate a robust model for 175 taxa of stream macroinvertebrates and to project their distributions on the Kinzig River stream network using bioclimatic, topographical, hydrological, land use and geological predictors. On average, model performance was good, with a TSS of 0.83 (±0.09 SD) and a ROC of 0.95 (±0.03 SD). The model delivered valuable insights on three sources of bias that plague SDMs in general: (a) level of taxonomic identification of the modeled organisms, (b) the spatial arrangement of sampling sites, and (c) the sampling intensity at each sampling site. Taxonomic resolution did not affect SDM performance. The distribution of high predicted probabilities of occurrence in the stream network coincided with those segments in the stream network most densely and frequently sampled, indicating both a spatial and temporal sampling bias. Species richness curves confirmed the temporal sampling bias. Next to spatial bias, sampling frequency also plays an important role in data collection, affecting further analysis and modeling procedures. Results indicate an underrepresentation of low order streams, an important aspect that should be addressed by both monitoring schemes and modeling approaches.     

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.     

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.     

From topography to hydrology—The modifiable area unit problem impacts freshwater species distribution models

Species distribution models (SDMs) are statistical tools to identify potentially suitable habitats for species. For SDMs in river ecosystems, species occurrences and predictor data are often aggregated across subcatchments that serve as modeling units. The level of aggregation (i.e., model resolution) influences the statistical relationships between species occurrences and environmental predictors—a phenomenon known as the modifiable area unit problem (MAUP), making model outputs directly contingent on the model resolution. Here, we test how model performance, predictor importance, and the spatial congruence of species predictions depend on the model resolution (i.e., average subcatchment size) of SDMs. We modeled the potential habitat suitability of 50 native fish species in the upper Danube catchment at 10 different model resolutions. Model resolutions were derived using a 90‐m digital‐elevation model by using the GRASS‐GIS module r.watershed. Here, we decreased the average subcatchment size gradually from 632 to 2 km². We then ran ensemble SDMs based on five algorithms using topographical, climatic, hydrological, and land‐use predictors for each species and resolution. Model evaluation scores were consistently high, as sensitivity and True Skill Statistic values ranged from 86.1–93.2 and 0.61–0.73, respectively. The most contributing predictor changed from topography at coarse, to hydrology at fine resolutions. Climate predictors played an intermediate role for all resolutions, while land use was of little importance. Regarding the predicted habitat suitability, we identified a spatial filtering from coarse to intermediate resolutions. The predicted habitat suitability within a coarse resolution was not ported to all smaller, nested subcatchments, but only to a fraction that held the suitable environmental conditions. Across finer resolutions, the mapped predictions were spatially congruent without such filter effect. We show that freshwater SDM predictions can have consistently high evaluation scores while mapped predictions differ significantly and are highly contingent on the underlying subcatchment size. We encourage building freshwater SDMs across multiple catchment sizes, to assess model variability and uncertainties in model outcomes emerging from the MAUP.     

Predicting fundamental and realized distributions based on thermal niche: A case study of a freshwater turtle

Species distribution models (SDM) have been broadly used in ecology to address theoretical and practical problems. Currently, there are two main approaches to generate SDMs: (i) correlative, which is based on species occurrences and environmental predictor layers and (ii) process-based models, which are constructed based on species' functional traits and physiological tolerances. The distributions estimated by each approach are based on different components of species niche. Predictions of correlative models approach species realized niches, while predictions of process-based are more akin to species fundamental niche. Here, we integrated the predictions of fundamental and realized distributions of the freshwater turtle Trachemys dorbigni. Fundamental distribution was estimated using data of T. dorbigni's egg incubation temperature, and realized distribution was estimated using species occurrence records. Both types of distributions were estimated using the same regression approaches (logistic regression and support vector machines), both considering macroclimatic and microclimatic temperatures. The realized distribution of T. dorbigni was generally nested in its fundamental distribution reinforcing theoretical assumptions that the species' realized niche is a subset of its fundamental niche. Both modelling algorithms produced similar results but microtemperature generated better results than macrotemperature for the incubation model. Finally, our results reinforce the conclusion that species realized distributions are constrained by other factors other than just thermal tolerances.     

In search of relevant predictors for marine species distribution modelling using the MarineSPEED benchmark dataset

Aim

Ideally, datasets for species distribution modelling (SDM) contain evenly sampled records covering the entire distribution of the species, confirmed absences and auxiliary ecophysiological data allowing informed decisions on relevant predictors. Unfortunately, these criteria are rarely met for marine organisms for which distributions are too often only scantly characterized and absences generally not recorded. Here, we investigate predictor relevance as a function of modelling algorithms and settings for a global dataset of marine species.

Location

Global marine.

Methods

We selected well‐studied and identifiable species from all major marine taxonomic groups. Distribution records were compiled from public sources (e.g., OBIS, GBIF, Reef Life Survey) and linked to environmental data from Bio‐ORACLE and MARSPEC. Using this dataset, predictor relevance was analysed under different variations of modelling algorithms, numbers of predictor variables, cross‐validation strategies, sampling bias mitigation methods, evaluation methods and ranking methods. SDMs for all combinations of predictors from eight correlation groups were fitted and ranked, from which the top five predictors were selected as the most relevant.

Results

We collected two million distribution records from 514 species across 18 phyla. Mean sea surface temperature and calcite are, respectively, the most relevant and irrelevant predictors. A less clear pattern was derived from the other predictors. The biggest differences in predictor relevance were induced by varying the number of predictors, the modelling algorithm and the sample selection bias correction. The distribution data and associated environmental data are made available through the R package marinespeed and at http://marinespeed.org .

Main conclusions

While temperature is a relevant predictor of global marine species distributions, considerable variation in predictor relevance is linked to the SDM set‐up. We promote the usage of a standardized benchmark dataset (MarineSPEED) for methodological SDM studies.

     

Scaling of species distribution models across spatial resolutions and extents along a biogeographic gradient. The case of the Iberian mole Talpa occidentalis

Scaling is a key process in modelling approaches since it allows for translating information from one scale to another. However, the success of this procedure may depend on ‘source’ and ‘target’ scales, but also on the biogeographic/ecological context of the study area. We aimed to quantify the performance and success of scaling species distribution model (SDM) predictions across spatial resolution and extent along a biogeographic gradient using the Iberian mole as study case. We ran separate MaxEnt models at two extents (national and regional) using independent datasets (species locations and environmental predictors) collected at 10 km and 50 m resolutions respectively. Model performance and success of scaling SDMs were quantified on the basis of accuracy measures and spatial predictions. Complementarily, we calculated marginality and tolerance as indicators of habitat availability and niche truncation along the biogeographic gradient. Model performance increased with resolution and extent, as well as from north to south (mainly for high resolution models). When regional models were validated at different scales, their performance reduced severely, particularly in the case of coarse resolution models (some of them performed worse than random). However, when the 10 km‐national model was downscaled within regions, it performed better (AUC_test: 0.82, 0.85 and 0.55 respectively for Galicia, Madrid and Granada) than models specifically calibrated within each region at 10 km (0.47, 0.65, 0.44). Indeed, it also had a better accuracy when projected at 50 m (0.77, 0.91, 0.79) than models fitted at that resolution (0.62, 0.83, 0.96) in two of the three cases. The success of scaling model predictions decreased along the biogeographic gradient, being these differences associated to niche truncation. Models representing non‐truncated niches were more successfully scaled across resolutions and extents (particularly in areas not offering all possible habitats for species), which has important implications for SDM applications.     

National-scale conservation assessments at an appropriate resolution

Abstract.  Most national-scale conservation assessments are carried out at a resolution that is different from the actual size of protected areas in the study region. Coincidence between nature reserves and both hotspots (areas of high species richness) and complementary areas (sets of sites within which all species are represented) have been reported. However, the resolution (size of grid cells) of the species' distribution data upon which many of these studies are based is often close to an order of magnitude larger than the size of the reserves. Presumably, only a proportion of the species recorded in the coarse grid cells will actually be present on reserves. We use fine (2 × 2-km square grid cells) and coarse (10 × 10-km square grid cells) resolution data of national distributions for breeding birds throughout Great Britain, and presence data for avian species on Royal Society for the Protection of Birds (RSPB) nature reserves, to investigate the proportion of species in the local area (100 km²) that are actually present on reserves. RSPB reserves contain between 50% and 70% of species from the local area. These proportions are significantly higher than for randomly selected, non-reserve areas, indicating that RSPB reserves contain higher concentrations of bird species than the wider countryside. Furthermore, on RSPB reserves these proportions of threatened and non-threatened species are equal, whereas in non-reserve areas the proportions of non-threatened species are significantly higher than threatened species. Thus, reserves hold a higher proportion of threatened species than occurs in the wider countryside.     

Species distribution modelling as a macroecological tool: a case study using New World amphibians

Although species distribution modelling (SDM) is widely accepted among the scientific community and is increasingly used in ecology, conservation biology and biogeography, methodological limitations generate potential problems for its application in macroecology. Using amphibian species richness in North and South America, we compare species richness patterns derived from SDM maps and ‘expert’ maps to evaluate if: 1) richness patterns derived from SDM are biased toward climate‐based explanations for diversity when compared to expert maps, since SDM methods are typically based on climatic variables; and 2) SDM is a reliable tool for generating richness maps in hyperrich regions where point occurrence data are limited for many species. We found that although three widely used SDM methods overestimated amphibian species richness in grid cells when compared to expert richness maps in both North and South America due to systematic overestimation of range sizes, diversity gradients were reasonably robust at broad scales. Further, climatic variables statistically explained patterns of richness at similar levels among the different richness sources, although climatic relationships were stronger in the much better known North America than in South America. We conclude that in the face of the high deforestation rates coupled with incomplete data on species distributions, especially in the tropics, SDM represents a useful macroecological tool for investigating broad‐scale richness patterns and the dynamics between species richness and climate.     

Effects of ecohydrological variables on current and future ranges,local suitability patterns,and model accuracy in big sagebrush

Forecasting of species and ecosystem responses to novel conditions, including climate change, is one of the major challenges facing ecologists at the start of the 21st century. Climate change studies based on species distribution models (SDMs) have been criticized because they extend correlational relationships beyond the observed data. Here, we compared conventional climate‐based SDMs against ecohydrological SDMs that include information from process‐based simulations of water balance. We examined the current and future distribution of Artemisia tridentata (big sagebrush) representing sagebrush ecosystems, which are widespread in semiarid western North America. For each approach, we calculated ensemble models from nine SDM methods and tested accuracy of each SDM with a null distribution. Climatic conditions included current conditions for 1970–1999 and two IPCC projections B1 and A2 for 2070–2099. Ecohydrological conditions were assessed by simulating soil water balance with SOILWAT, a daily time‐step, multiple layer, mechanistic, soil water model. Under current conditions, both climatic and ecohydrological SDM approaches produced comparable sagebrush distributions. Overall, sagebrush distribution is forecasted to decrease, with larger decreases under the A2 than under the B1 scenario and strong decreases in the southern part of the range. Increases were forecasted in the northern parts and at higher elevations. Both SDM approaches produced accurate predictions. However, the ecohydrological SDM approach was slightly less accurate than climatic SDMs (?1% in AUC, ?4% in Kappa and TSS) and predicted a higher number of habitat patches than observed in the input data. Future predictions of ecohydrological SDMs included an increased number of habitat patches whereas climatic SDMs predicted a decrease. This difference is important for understanding landscape‐scale patterns of sagebrush ecosystems and management of sagebrush obligate species for future conditions. Several mechanisms can explain the diverging forecasts; however, we need better insights into the consequences of different datasets for SDMs and how these affect our understanding of future trajectories.     

Decomposing the heterogeneity of species distributions into multiple scales: a hierarchical framework for large‐scale count surveys

We introduce a novel spatially explicit framework for decomposing species distributions into multiple scales from count data. These kinds of data are usually positively skewed, have non‐normal distributions and are spatially autocorrelated. To analyse such data, we propose a hierarchical model that takes into account the observation process and explicitly deals with spatial autocorrelation. The latent variable is the product of a positive trend representing the non‐constant mean of the species distribution and of a stationary positive spatial field representing the variance of the spatial density of the species distribution. Then, the different scales of emergent structures of the distribution of the population in space are modelled from the latent density of the species distribution using multi‐scale variogram models. Multi‐scale kriging is used to map the spatial patterns previously identified by the multi‐scale models. We show how our framework yields robust and precise estimates of the relevant scales both for spatial count data simulated from well‐defined models, and in a real case‐study based on seabird count data (the common guillemot Uria aalge) provided by large‐scale aerial surveys of the Bay of Biscay (France) performed over a winter. Our stochastic simulation study provides guidelines on the expected uncertainties of the scales estimates. Our results indicate that the spatial structure of the common guillemot can be modelled as a three‐level hierarchical system composed of a very broad‐scale pattern (～ 200 km) with a stable location over time that might be environmentally controlled, a broad‐scale pattern (～ 50 km) with a variable shape and location, that might be related to shifts in prey distribution, and a fine‐scale pattern (～ 10 km) with a rather stable shape and location, that might be controlled by behavioural processes. Our framework enables the development of robust, scale‐dependent hypotheses regarding the potential ecological processes that control species distributions.     

At-sea surveys confirm a North Atlantic biodiversity hotspot

Capsule A number of tracking studies have shown an area of the North Atlantic, south of the Charlie-Gibbs Fracture Zone, to be an important overwintering location for seabirds. We conducted seabird observations along a trans-Atlantic transect from Ireland to Canada in April 2014 to test the hypothesis that seabird species richness and abundance will peak in the area of the Mid-Atlantic Ridge. At-sea survey results agreed with previous tracking studies, highlighting the importance of the Mid-Atlantic Ridge area for seabirds.