Jointly estimating spatial sampling effort and habitat suitability for multiple species from opportunistic presence‐only data

Christophe Botella (1) , Alexis Joly (2) , Pierre Bonnet (1) , François Munoz (3) , Pascal Monestiez (4, 5)

1 UMR AMAP - Botanique et Modélisation de l'Architecture des Plantes et des Végétations
2 ZENITH - Scientific Data Management
3 LECA - Laboratoire d'Ecologie Alpine
4 INRAE - Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement
5 BioSP - Biostatistique et Processus Spatiaux

1. Building reliable species distribution models (SDMs) from presence-only information requires a good understanding of the spatial variation in the sampling effort. However, in most cases, the sampling effort is unknown, leading to biases in SDMs. This study proposes a method to jointly estimate the parameters of sampling effort and species densities to avoid such biases. The method is particularly suited to the analysis of massive but highly heterogeneous presence-only data. 2. The proposed method is based on estimating the variation in sampling effort over units of a spatial mesh in parallel with the environmental density of multiple species using a marked Poisson process model. Based on simulations with realistic settings, we examined the performance and robustness of parameter estimations. We also analysed a large-scale citizen science dataset with highly heterogeneous sampling (Pl@ntNet), including around 300,000 occurrences of 150 plant species. 3. We found that sampling effort was correctly estimated when the true sampling effort was constant within the cells of a spatial mesh. Estimation bias arose when sampling effort and environmental drivers strongly covaried within cells. Otherwise, the inference was correct and robust to sampling variation within cells. Running the model on real occurrences of 150 plant species provided an estimated map of relative sampling effort for 15% of French territory. We also found that the density estimated for an exotic invasive plant was consistent with prior data. 4. This is the first method jointly estimating species densities depending on environment, and sampling effort as an explicit spatial function, from occurrence data of multiple species. An asset of the method is that a few frequently observed species greatly contribute to correctly estimate sampling effort, thereby improving density estimation of all other species. This approach can thus provide reliable SDM for This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Domaines: Sciences du Vivant [q-bio]; Sciences de l'environnement; Mathématiques [math]; Informatique [cs]

Christophe Botella, Alexis Joly, Pierre Bonnet, François Munoz, Pascal Monestiez. Jointly estimating spatial sampling effort and habitat suitability for multiple species from opportunistic presence‐only data. Methods in Ecology and Evolution, 2021, 12 (5), pp.933-945. ⟨10.1111/2041-210X.13565⟩. ⟨hal-03150701⟩