The Biodiversity and Climate Change Virtual Laboratory: Where ecology meets big data

Advances in computing power and infrastructure, increases in the number and size of ecological and environmental datasets, and the number and type of data collection methods, are revolutionizing the field of Ecology. To integrate these advances, virtual laboratories offer a unique tool to facilitate, expedite, and accelerate research into the impacts of climate change on biodiversity. We introduce the uniquely cloud-based Biodiversity and Climate Change Virtual Laboratory (BCCVL), which provides access to numerous species distribution modelling tools; a large and growing collection of biological, climate, and other environmental datasets; and a variety of experiment types to conduct research into the impact of climate change on biodiversity. Users can upload and share datasets, potentially increasing collaboration, cross-fertilisation of ideas, and innovation among the user community. Feedback confirms that the BCCVL's goals of lowering the technical requirements for species distribution modelling, and reducing time spent on such research, are being met

Recommender systems meet species distribution modelling

Publisher Copyright: © 2021 Copyright for this paper by its authors.Recommender systems techniques can naturally lend themselves to species distribution modelling if biological species are treated as items and places where they occur are treated as users. In this setting recommendation scores can reflect which habitats are suited for which species. Recommendation scores can also be used for reconstructing relative abundances of species, and analysing their rises and declines over millions of years in the past. Analysis of such predictions can shed light on the effects of changing environments on the biosphere now and in the past, as well as help to make predictions for the future. The major potential advantage of the recommender systems treatment over many existing solutions is the large spatial and temporal scale at which such analysis can be done within a single model. A single model makes predictions easier to compare globally in space and over time. While algorithmic application of recommender systems techniques to species distribution modelling is relatively straightforward, model selection and evaluation is particularly challenging, as there is no possibility for online tests or on-demand sampling, since the past worlds are long gone. Explainability is paramount in these tasks. Here we highlight the main challenges and promising directions of evaluation of such modelling, which is still in early stages of development. We show how aggregated prediction statistics and constraints may help for reliable model selection and evaluation. We illustrate the approaches on a case study of the mammalian fossil record from Europe around 8-17 millions of years ago.Peer reviewe

Incorporating knowledge uncertainty into species distribution modelling

Monitoring progress towards global goals and biodiversity targets require reliable descriptions of species distributions over time and space. Current gaps in accessible information on species distributions urges the need for integrating all available data and knowledge sources, and intensifying cooperations to more effectively support global environmental governance. For many areas and species groups, experts can constitute a valuable source of information to fill the gaps by offering their knowledge on species-environment interactions. However, expert knowledge is always subject to uncertainty, and incorporating that into species distribution mapping poses a challenge. We propose the use of the dempster–shafer theory of evidence (DST) as a novel approach in this field to extract expert knowledge, to incorporate the associated uncertainty into the procedure, and to produce reliable species distribution maps. We applied DST to model the distribution of two species of eagle in Spain. We invited experts to fill in an online questionnaire and express their beliefs on the habitat of the species by assigning probability values for given environmental variables, along with their confidence in expressing the beliefs. We then calculated evidential functions, and combined them using Dempster’s rules of combination to map the species distribution based on the experts’ knowledge. We evaluated the performances of our proposed approach using the atlas of Spanish breeding birds as an independent test dataset, and further compared the results with the outcome of an ensemble of conventional SDMs. Purely based on expert knowledge, the DST approach yielded similar results as the data driven SDMs ensemble. Our proposed approach offers a strong and practical alternative for species distribution modelling when species occurrence data are not accessible, or reliable, or both. The particular strengths of the proposed approach are that it explicitly accounts for and aggregates knowledge uncertainty, and it capitalizes on the range of data sources usually considered by an expert

Delimiting the geographical background in species distribution modelling

[Aim]: The extent of the study area (geographical background, GB) can strongly affect the results of species distribution models (SDMs), but as yet we lack objective and practicable criteria for delimiting the appropriate GB. We propose an approach to this problem using trend surface analysis (TSA) and provide an assessment of the effects of varying GB extent on the performance of SDMs for four species. [Location]: Mainland Spain. [Methods]: Using data for four well known wild ungulate species and different GBs delimited with a TSA, we assessed the effects of GB extent on the predictive performance of SDMs: specifically on model calibration (Miller's statistic) and discrimination (area under the curve of the receiver operating characteristic plot, AUC; sensitivity and specificity), and on the tendency of the models to predict environmental potential when they are projected beyond their training area. [Results]: In the training area, discrimination significantly increased and calibration decreased as the GB was enlarged. In contrast, as GB was enlarged, both discriminatory power and calibration decreased when assessed in the core area of the species distributions. When models trained using small GBs were projected beyond their training area, they showed a tendency to predict higher environmental potential for the species than those models trained using large GBs. [Main conclusions]: By restricting GB extent using a geographical criterion, model performance in the core area of the species distribution can be significantly improved. Large GBs make models demonstrate high discriminatory power but are barely informative. By delimiting GB using a geographical criterion, the effect of historical events on model parameterization may be reduced. Thus purely environmental models are obtained that, when projected onto a new scenario, depict the potential distribution of the species. We therefore recommend the use of TSA in geographically delimiting the GB for use in SDMs.P.A. and A.J.-V. were supported by the Juan de la Cierva research program awarded by the Ministerio de Ciencia e Innovación – Fondo Social Europeo, and partly by the project CGL2009-11316/BOS – Fondos FEDER. P.A. is in Portugal thanks to a José Castillejo fellowship (2010–11) granted by the Ministerio de Ciencia e Innovación.Peer Reviewe

Species distribution modelling of Aloidendron dichotomum (quiver tree)

A variety of species distribution models (SDMs) were fit to data collected by a 15,000km road-side visual survey of Aloidendron dichotomum populations in the Northern Cape region of South Africa, and Namibia. We fit traditional presence/absence SDMs as well as SDMs on how proportions are distributed across three species stage classes (juvenile, adult, dead). Using five candidate machine learning methods and an ensemble model, we compared a number of approaches, including the role of balanced class (presence/absence) datasets in species distribution modelling. Secondary to this was whether or not the addition of species’ absences, generated where the species is known not to exist have an impact on findings. The goal of the analysis was to map the distribution of Aloidendron dichotomum under different scenarios. Precipitation-based variables were generally more deterministic of species presence or lack thereof. Visual interpretation of the estimated Aloidendron dichotomum population under current climate conditions, suggested a reasonably well fit model, having a large overlap with the sampled area. There however were some conditions estimated to be suitable for species incidence outside of the sampled range, where Aloidendron dichotomum are not known to occur. Habitat suitability for juvenile individuals was largely decreasing in concentration towards Windhoek. The largest proportion of dead individuals was estimated to be on the northern edge of the Riemvasmaak Conservancy, along the South African/Namibian boarder, reaching up to a 60% composition of the population. The adult stage class maintained overall proportional dominance. Under future climate scenarios, despite maintaining a bulk of the currently habitable conditions, a noticeable negative shift in habitat suitability for the species was observed. A temporal analysis of Aloidendron dichotomum’s latitudinal and longitudinal range revealed a potential south-easterly shift in suitable species conditions. Results were however met with some uncertainty as SDMs were uncovered to be extrapolating into a substantial amount of the study area. We found that balancing response class frequencies within the data proved not to be an effective error reduction technique overall, having no considerable impact on species detection accuracy. Balancing the classes however did improve the accuracy on the presence class, at the cost of accuracy of the observed absence class. Furthermore, overall model accuracy increased as more absences from outside the study area were added, only because these generated absences were predicted well. The resulting models had lower estimated suitability outside of the survey area and noticeably different suitability distributions within the survey area. This made the addition of the generated absences undesirable. Results highlighted the potential vulnerability of Aloidendron dichotomum given the pessimistic, yet likely future climate scenarios

Species distribution modelling to support marine conservation planning

This thesis explores some important practical considerations concerning the use of species distribution models in marine conservation planning. Using geo-referenced gorgonian distribution data, together with explanatory environmental variables, predictive models have been used to map the spatial distribution of suitable gorgonian (sea fan) habitat in two study sites; Hatton Bank, in the Northeast Atlantic, and Lyme Bay on the south coast of Devon. Generalized Linear Models (GLMs), Generalized Additive Models (GAMs) and a Maximum Entropy (Maxent) model have been used to support critical investigation into important model considerations that have received inadequate attention in the marine environment. The influence of environmental data resolution on model performance has been explored with specific reference to available datasets in the nearshore and offshore environments. The transferability of deep-sea models has been similarly appraised, with recommendations as to the appropriate use of transferred models. Investigating these practical issues will allow managers to make informed decisions with respect to the best and most appropriate use of existing data. This study has also used novel approaches and investigated their suitability for marine conservation planning, including the use of model classification error in the spatial prioritisation of monitoring sites, and the adaptation of an existing presence-only modelling method to include absence data. Together, these studies contribute both practical recommendations for marine conservation planning and novel applications within the wider species distribution modelling discipline, and consider the implications of these developments for managers, to ensure the ongoing improvement and development of models to support conservation planning

Joint species distribution modelling with the r-package Hmsc

Joint Species Distribution Modelling (JSDM) is becoming an increasingly popular statistical method for analysing data in community ecology. Hierarchical Modelling of Species Communities (HMSC) is a general and flexible framework for fitting JSDMs. HMSC allows the integration of community ecology data with data on environmental covariates, species traits, phylogenetic relationships and the spatio-temporal context of the study, providing predictive insights into community assembly processes from non-manipulative observational data of species communities. The full range of functionality of HMSC has remained restricted to Matlab users only. To make HMSC accessible to the wider community of ecologists, we introduce Hmsc 3.0, a user-friendly r implementation. We illustrate the use of the package by applying Hmsc 3.0 to a range of case studies on real and simulated data. The real data consist of bird counts in a spatio-temporally structured dataset, environmental covariates, species traits and phylogenetic relationships. Vignettes on simulated data involve single-species models, models of small communities, models of large species communities and models for large spatial data. We demonstrate the estimation of species responses to environmental covariates and how these depend on species traits, as well as the estimation of residual species associations. We demonstrate how to construct and fit models with different types of random effects, how to examine MCMC convergence, how to examine the explanatory and predictive powers of the models, how to assess parameter estimates and how to make predictions. We further demonstrate how Hmsc 3.0 can be applied to normally distributed data, count data and presence-absence data. The package, along with the extended vignettes, makes JSDM fitting and post-processing easily accessible to ecologists familiar with r.Peer reviewe

Joint species distribution modelling with the r-package Hmsc

Joint Species Distribution Modelling (JSDM) is becoming an increasingly popular statistical method for analysing data in community ecology. Hierarchical Modelling of Species Communities (HMSC) is a general and flexible framework for fitting JSDMs. HMSC allows the integration of community ecology data with data on environmental covariates, species traits, phylogenetic relationships and the spatio-temporal context of the study, providing predictive insights into community assembly processes from non-manipulative observational data of species communities. The full range of functionality of HMSC has remained restricted to Matlab users only. To make HMSC accessible to the wider community of ecologists, we introduce Hmsc 3.0, a user-friendly r implementation. We illustrate the use of the package by applying Hmsc 3.0 to a range of case studies on real and simulated data. The real data consist of bird counts in a spatio-temporally structured dataset, environmental covariates, species traits and phylogenetic relationships. Vignettes on simulated data involve single-species models, models of small communities, models of large species communities and models for large spatial data. We demonstrate the estimation of species responses to environmental covariates and how these depend on species traits, as well as the estimation of residual species associations. We demonstrate how to construct and fit models with different types of random effects, how to examine MCMC convergence, how to examine the explanatory and predictive powers of the models, how to assess parameter estimates and how to make predictions. We further demonstrate how Hmsc 3.0 can be applied to normally distributed data, count data and presence-absence data. The package, along with the extended vignettes, makes JSDM fitting and post-processing easily accessible to ecologists familiar with r.Peer reviewe