A virtual species set for robust and reproducible species distribution modelling tests

Predicting species potential and future distribution has become a relevant tool in biodiversity monitoring and conservation.In this data article we present the suitability map of a virtual species generated based on two bioclimatic variables, and a dataset containing more than 700,000 random observations at the extent of Europe. The dataset includes spatial attributes such as: distance to roads, protected areas, country codes, and the habitat suitability of two spatially clustered species (grassland and forest species) and a wide-spread species.

Modelling the influence of biotic factors on species distribution patterns

Biotic interactions can have large effects on species distributions yet their role in shaping species ranges is seldom explored due to historical difficulties in incorporating biotic factors into models without a priori knowledge on interspecific interactions. Improved SDMs, which account for biotic factors and do not require a priori knowledge on species interactions, are needed to fully understand species distributions. Here, we model the influence of abiotic and biotic factors on species distribution patterns and explore the robustness of distributions under future climate change. We fit hierarchical spatial models using Integrated Nested Laplace Approximation (INLA) for lagomorph species throughout Europe and test the predictive ability of models containing only abiotic factors against models containing abiotic and biotic factors. We account for residual spatial autocorrelation using a conditional autoregressive (CAR) model. Model outputs are used to estimate areas in which abiotic and biotic factors determine species’ ranges. INLA models containing both abiotic and biotic factors had substantially better predictive ability than models containing abiotic factors only, for all but one of the four species. In models containing abiotic and biotic factors, both appeared equally important as determinants of lagomorph ranges, but the influences were spatially heterogeneous. Parts of widespread lagomorph ranges highly influenced by biotic factors will be less robust to future changes in climate, whereas parts of more localised species ranges highly influenced by the environment may be less robust to future climate. SDMs that do not explicitly include biotic factors are potentially misleading and omit a very important source of variation. For the field of species distribution modelling to advance, biotic factors must be taken into account in order to improve the reliability of predicting species distribution patterns both presently and under future climate change.

An examination of the usefulness of single-species biomass models for the management of Lake Victoria fisheries

Nile perch (Lates niloticus), tilapia (Oreochromis spp), dagaa (Rastrineobola argentea, silver cyprinid), and haplochromines (Tribe Haplochromini) form the backbone of the commercial fishery on Lake Victoria. These fish stocks account for about 70% of the total catch in the three riparian states Uganda, Kenya, and Tanzania. The lake fisheries have been poorly managed, in part due to inadequate scientific analysis and management advice. The overall objective of this project was to model the stocks of the commercial fisheries of Lake Victoria with the view of determining reference points and current stock status. The Schaefer biomass model was fitted to available data for each stock (starting in the 1960s or later) in the form of landings, catch per unit effort, acoustic survey indices, and trawl survey indices. In most cases, the Schaefer model did not fit all data components very well, but attempts were made to find the best model for each stock. When the model was fitted to the Nile perch data starting from 1996, the estimated current biomass is 654 kt (95% CI 466–763); below the optimum of 692 kt and current harvest rate is 38% (33–73%), close to the optimum of 35%. At best, these can be used as tentative guidelines for the management of these fisheries. The results indicate that there have been strong multispecies interactions in the lake ecosystem. The findings from our study can be used as a baseline reference for future studies using more complex models, which could take these multispecies interactions into account.

Patterns of Cactaceae species distribution in a protected area in the semiarid caatinga biome of North-Eastern Brazil.

2016

Niche conservatism in non-native birds in Europe: niche unfilling rather than niche expansion

Aim Niche conservatism, or the extent to which niches are conserved across space and time, is of special concern for the study of non-native species as it underlies predictions of invasion risk. Based on the occurrence of 28 non-native birds in Europe, we assess to what extent Grinnellian realized niches are conserved during invasion, formulate hypotheses to explain the variation in observed niche changes and test how well species distribution models can predict non-native bird occurrence in Europe. Location Europe. Methods To quantify niche changes, a recent method that applies kernel smoothers to densities of species occurrence in a gridded environmental space was used. This corrects for differences in the availability of environments between study areas and allows discrimination between 'niche expansion' into environments new to the species and 'niche unfilling', whereby the species only partially fills its niche in the invaded range. Predictions of non-native bird distribution in Europe were generated using several distribution modelling techniques. Results Niche overlap between native and non-native bird populations is low, but niche changes are smaller for species having a higher propagule pressure and that were introduced longer ago. Non-native birds in Europe occupy a subset of the environments they inhabit in their native ranges. Niche expansion into novel environments is rare for most species, allowing species distribution models to accurately predict invasion risk. Main conclusions Because of the recent nature of most bird introductions, species occupy only part of the suitable environments available in the invaded range. This signals that apart from purely ecological factors, patterns of niche conservatism may also be contingent on population-specific historical factors. These results also suggest that many claims of niche differences may be due to a partial filling of the native niche in the invaded range and thus do not represent true niche changes.

Niche dynamics in space and time.

Niche conservatism, the tendency of a species niche to remain unchanged over time, is often assumed when discussing, explaining or predicting biogeographical patterns. Unfortunately, there has been no basis for predicting niche dynamics over relevant timescales, from tens to a few hundreds of years. The recent application of species distribution models (SDMs) and phylogenetic methods to analysis of niche characteristics has provided insight to niche dynamics. Niche shifts and conservatism have both occurred within the last 100 years, with recent speciation events, and deep within clades of species. There is increasing evidence that coordinated application of these methods can help to identify species which likely fulfill one key assumption in the predictive application of SDMs: an unchanging niche. This will improve confidence in SDM-based predictions of the impacts of climate change and species invasions on species distributions and biodiversity.

Suitability, success and sinks: how do predictions of nesting distributions relate to fitness parameters in high arctic waders?

AimAlthough habitat suitability maps derived from species distribution models (SDMs) are often assumed to highlight locations that can sustain healthy populations over time, the relationship between suitability scores and fitness parameters has rarely been tested thoroughly. LocationZackenberg Valley, north-east Greenland. MethodsUsing 14years of data (1997-2010) representing three wader species (dunlin Calidris alpina, sanderling Calidris alba and ruddy turnstone Arenaria interpres), we tested the relationships between modelled suitability and fitness parameters at nesting locations. ResultsAmong the three species examined, only the ruddy turnstone exhibited significant relationships between suitability and nest success, but over time rather than space. During years with extensive snow cover in the landscape, the nesting sites of ruddy turnstone occurred in different habitats than were typically used across years. Moreover, in years with extensive snow cover, the ruddy turnstone initiated nests later and suffered from higher egg predation rates. Main conclusionOur results suggest that SDMs derived from species occurrences that include years of low reproductive success may over-estimate the potential suitable habitat in the landscape. Whenever possible, variation in reproductive success should be considered when building models to inform species' response to environmental change. species' response to environmental change.

Combining palaeodistribution modelling and phylogeographical approaches for identifying glacial refugia in Alpine Primula

Aim We investigated the late Quaternary history of two closely related and partly sympatric species of Primula from the south-western European Alps, P. latifolia Lapeyr. and P. marginata Curtis, by combining phylogeographical and palaeodistribution modelling approaches. In particular, we were interested in whether the two approaches were congruent and identified the same glacial refugia. Location South-western European Alps. Methods For the phylogeographical analysis we included 353 individuals from 28 populations of P. marginata and 172 individuals from 15 populations of P. latifolia and used amplified fragment length polymorphisms (AFLPs). For palaeodistribution modelling, species distribution models (SDMs) were based on extant species occurrences and then projected to climate models (CCSM, MIROC) of the Last Glacial Maximum (LGM), approximately 21 ka. Results The locations of the modelled LGM refugia were confirmed by various indices of genetic variation. The refugia of the two species were largely geographically isolated, overlapping only 6% to 11% of the species' total LGM distribution. This overlap decreased when the position of the glacial ice sheet and the differential elevational and edaphic distributions of the two species were considered. Main conclusions The combination of phylogeography and palaeodistribution modelling proved useful in locating putative glacial refugia of two alpine species of Primula. The phylogeographical data allowed us to identify those parts of the modelled LGM refugial area that were likely source areas for recolonization. The use of SDMs predicted LGM refugial areas substantially larger and geographically more divergent than could have been predicted by phylogeographical data alone

Developing methodologies to allow the robust prediction dynamic context dominated by perturbations

El principal objectiu del projecte era desenvolupar millores conceptuals i metodològiques que permetessin una millor predicció dels canvis en la distribució de les espècies (a una escala de paisatge) derivats de canvis ambientals en un context dominat per pertorbacions. En un primer estudi, vàrem comparar l'eficàcia de diferents models dinàmics per a predir la distribució de l'hortolà (Emberiza hortulana). Els nostres resultats indiquen que un model híbrid que combini canvis en la qualitat de l'hàbitat, derivats de canvis en el paisatge, amb un model poblacional espacialment explícit és una aproximació adequada per abordar canvis en la distribució d'espècies en contextos de dinàmica ambiental elevada i una capacitat de dispersió limitada de l'espècie objectiu. En un segon estudi abordarem la calibració mitjançant dades de seguiment de models de distribució dinàmics per a 12 espècies amb preferència per hàbitats oberts. Entre les conclusions extretes destaquem: (1) la necessitat de que les dades de seguiment abarquin aquelles àrees on es produeixen els canvis de qualitat; (2) el biaix que es produeix en la estimació dels paràmetres del model d'ocupació quan la hipòtesi de canvi de paisatge o el model de qualitat d'hàbitat són incorrectes. En el darrer treball estudiarem el possible impacte en 67 espècies d’ocells de diferents règims d’incendis, definits a partir de combinacions de nivells de canvi climàtic (portant a un augment esperat de la mida i freqüència d’incendis forestals), i eficiència d’extinció per part dels bombers. Segons els resultats dels nostres models, la combinació de factors antropogènics del regim d’incendis, tals com l’abandonament rural i l’extinció, poden ser més determinants per als canvis de distribució que els efectes derivats del canvi climàtic. Els productes generats inclouen tres publicacions científiques, una pàgina web amb resultats del projecte i una llibreria per a l'entorn estadístic R.

Understanding the ecological niche to elucidate spatial strategies of the southernmost Tupinambis lizards

Understanding factors that shape ranges of species is central in evolutionary biology. Species distribution models have become important tools to test biogeographical, ecological and evolutionary hypotheses. Moreover, from an ecological and evolutionary perspective, these models help to elucidate the spatial strategies of species at a regional scale. We modelled species distributions of two phylogenetically, geographically and ecologically close Tupinambis species (Teiidae) that occupy the southernmost area of the genus distribution in South America. We hypothesized that similarities between these species might have induced spatial strategies at the species level, such as niche differentiation and divergence of distribution patterns at a regional scale. Using logistic regression and MaxEnt we obtained species distribution models that revealed interspecific differences in habitat requirements, such as environmental temperature, precipitation and altitude. Moreover, the models obtained suggest that although the ecological niches of Tupinambis merianae and T. rufescens are different, these species might co-occur in a large contact zone. We propose that niche plasticity could be the mechanism enabling their co-occurrence. Therefore, the approach used here allowed us to understand the spatial strategies of two Tupinambis lizards at a regional scale.

Predicting present and future intra-specific genetic structure through niche hindcasting across 24 millennia.

Paleoclimatic reconstructions coupled with species distribution models and identification of extant spatial genetic structure have the potential to provide insights into the demographic events that shape the distribution of intra-specific genetic variation across time. Using the globeflower Trollius europaeus as a case-study, we combined (1) Amplified Fragment Length Polymorphisms, (2) suites of 1000-years stepwise hindcasted species distributions and (3) a model of diffusion through time over the last 24,000 years, to trace the spatial dynamics that most likely fits the species' current genetic structure. We show that the globeflower comprises four gene pools in Europe which, from the dry period preceding the Last Glacial Maximum, dispersed while tracking the conditions fitting its climatic niche. Among these four gene pools, two are predicted to experience drastic range retraction in the near future. Our interdisciplinary approach, applicable to virtually any taxon, is an advance in inferring how climate change impacts species' genetic structures.

Modelling ecological niches with support vector machines

1. The ecological niche is a fundamental biological concept. Modelling species' niches is central to numerous ecological applications, including predicting species invasions, identifying reservoirs for disease, nature reserve design and forecasting the effects of anthropogenic and natural climate change on species' ranges. 2. A computational analogue of Hutchinson's ecological niche concept (the multidimensional hyperspace of species' environmental requirements) is the support of the distribution of environments in which the species persist. Recently developed machine-learning algorithms can estimate the support of such high-dimensional distributions. We show how support vector machines can be used to map ecological niches using only observations of species presence to train distribution models for 106 species of woody plants and trees in a montane environment using up to nine environmental covariates. 3. We compared the accuracy of three methods that differ in their approaches to reducing model complexity. We tested models with independent observations of both species presence and species absence. We found that the simplest procedure, which uses all available variables and no pre-processing to reduce correlation, was best overall. Ecological niche models based on support vector machines are theoretically superior to models that rely on simulating pseudo-absence data and are comparable in empirical tests. 4. Synthesis and applications. Accurate species distribution models are crucial for effective environmental planning, management and conservation, and for unravelling the role of the environment in human health and welfare. Models based on distribution estimation rather than classification overcome theoretical and practical obstacles that pervade species distribution modelling. In particular, ecological niche models based on machine-learning algorithms for estimating the support of a statistical distribution provide a promising new approach to identifying species' potential distributions and to project changes in these distributions as a result of climate change, land use and landscape alteration.

Potential Impacts of Climate Change on Ecosystem Services in Europe: the Case of Pest Control by Vertebrates. BioScience.

Global environmental changes threaten ecosystems and cause significant alterations to the supply of ecosystem services that are vital for human well-being. We provide an assessment of the potential impacts of climate change on European diversity of vertebrates and their associated pest control services. We modeled the distributions of the species that provide this service using ensembles of forecasts from bioclimatic envelope models and then used their results to generate maps of potential species richness among vertebrate providers of pest control services. We assessed how potential richness of pest control providers would change according to different climate and greenhouse emissions scenarios. We found that potential richness of pest control providers was likely to face substantial reductions, especially in southern European countries that had economies highly dependent on agricultural yields. In much of central and northern Europe, where countries had their economies less dependent on agriculture, climate change was likely to benefit pest control providers

The future of terrestrial mammals in the Mediterranean basin under climate change.

The Mediterranean basin is considered a hotspot of biological diversity with a long history of modification of natural ecosystems by human activities, and is one of the regions that will face extensive changes in climate. For 181 terrestrial mammals (68% of all Mediterranean mammals), we used an ensemble forecasting approach to model the future (approx. 2100) potential distribution under climate change considering five climate change model outputs for two climate scenarios. Overall, a substantial number of Mediterranean mammals will be severely threatened by future climate change, particularly endemic species. Moreover, we found important changes in potential species richness owing to climate change, with some areas (e.g. montane region in central Italy) gaining species, while most of the region will be losing species (mainly Spain and North Africa). Existing protected areas (PAs) will probably be strongly influenced by climate change, with most PAs in Africa, the Middle East and Spain losing a substantial number of species, and those PAs gaining species (e.g. central Italy and southern France) will experience a substantial shift in species composition.

Does phylogeographic structure relate to climatic niche divergence? A test using maritime pine (Pinus pinaster Ait.)

Aim To disentangle the effects of environmental and geographical processes driving phylogenetic distances among clades of maritime pine (Pinus pinaster). To assess the implications for conservation management of combining molecular information with species distribution models (SDMs; which predict species distribution based on known occurrence records and on environmental variables). Location Western Mediterranean Basin and European Atlantic coast. Methods We undertook two cluster analyses for eight genetically defined pine clades based on climatic niche and genetic similarities. We assessed niche similarity by means of a principal component analysis and Schoener's D metric. To calculate genetic similarity, we used the unweighted pair group method with arithmetic mean based on Nei's distance using 266 single nucleotide polymorphisms. We then assessed the contribution of environmental and geographical distances to phylogenetic distance by means of Mantel regression with variance partitioning. Finally, we compared the projection obtained from SDMs fitted from the species level (SDMsp) and composed from the eight clade-level models (SDMcm). Results Genetically and environmentally defined clusters were identical. Environmental and geographical distances explained 12.6% of the phylogenetic distance variation and, overall, geographical and environmental overlap among clades was low. Large differences were detected between SDMsp and SDMcm (57.75% of disagreement in the areas predicted as suitable). Main conclusions The genetic structure within the maritime pine subspecies complex is primarily a consequence of its demographic history, as seen by the high proportion of unexplained variation in phylogenetic distances. Nevertheless, our results highlight the contribution of local environmental adaptation in shaping the lower-order, phylogeographical distribution patterns and spatial genetic structure of maritime pine: (1) genetically and environmentally defined clusters are consistent, and (2) environment, rather than geography, explained a higher proportion of variation in phylogenetic distance. SDMs, key tools in conservation management, better characterize the fundamental niche of the species when they include molecular information.