The recent past and promising future for data integration methods to estimate species’ distributions

Emerging technologies in citizen science and potential for insect monitoring

Emerging technologies are increasingly employed in environmental citizen science projects. This integration offers benefits and opportunities for scientists and participants alike. Citizen science can support large-scale, long-term monitoring of species occurrences, behaviour and interactions. At the same time, technologies can foster participant engagement, regardless of pre-existing taxonomic expertise or experience, and permit new types of data to be collected. Yet, technologies may also create challenges by potentially increasing financial costs, necessitating technological expertise or demanding training of participants. Technology could also reduce people's direct involvement and engagement with nature. In this perspective, we discuss how current technologies have spurred an increase in citizen science projects and how the implementation of emerging technologies in citizen science may enhance scientific impact and public engagement. We show how technology can act as (i) a facilitator of current citizen science and monitoring efforts, (ii) an enabler of new research opportunities, and (iii) a transformer of science, policy and public participation, but could also become (iv) an inhibitor of participation, equity and scientific rigour. Technology is developing fast and promises to provide many exciting opportunities for citizen science and insect monitoring, but while we seize these opportunities, we must remain vigilant against potential risks. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Integrated distance sampling models for simple point counts

Point counts (PCs) are widely used in biodiversity surveys but, despite numerous advantages, simple PCs suffer from several problems: detectability, and therefore abundance, is unknown; systematic spatiotemporal variation in detectability yields biased inferences, and unknown survey area prevents formal density estimation and scaling-up to the landscape level. We introduce integrated distance sampling (IDS) models that combine distance sampling (DS) with simple PC or detection/nondetection (DND) data to capitalize on the strengths and mitigate the weaknesses of each data type. Key to IDS models is the view of simple PC and DND data as aggregations of latent DS surveys that observe the same underlying density process. This enables the estimation of separate detection functions, along with distinct covariate effects, for all data types. Additional information from repeat or time-removal surveys, or variable survey duration, enables the separate estimation of the availability and perceptibility components of detectability with DS and PC data. IDS models reconcile spatial and temporal mismatches among data sets and solve the above-mentioned problems of simple PC and DND data. To fit IDS models, we provide JAGS code and the new "IDS()" function in the R package unmarked. Extant citizen-science data generally lack the information necessary to adjust for detection biases, but IDS models address this shortcoming, thus greatly extending the utility and reach of these data. In addition, they enable formal density estimation in hybrid designs, which efficiently combine DS with distance-free, point-based PC or DND surveys. We believe that IDS models have considerable scope in ecology, management, and monitoring.

A new approach to geostatistical synthesis of historical records reveals capuchin spatial responses to climate and demographic change

Recent proliferation of GPS technology has transformed animal movement research. Yet, time-series data from this recent technology rarely span beyond a decade, constraining longitudinal research. Long-term field sites hold valuable historic animal location records, including hand-drawn maps and semantic descriptions. Here, we introduce a generalised workflow for converting such records into reliable location data to estimate home ranges, using 30 years of sleep-site data from 11 white-faced capuchin (Cebus imitator) groups in Costa Rica. Our findings illustrate that historic sleep locations can reliably recover home range size and geometry. We showcase the opportunity our approach presents to resolve open questions that can only be addressed with very long-term data, examining how home ranges are affected by climate cycles and demographic change. We urge researchers to translate historical records into usable movement data before this knowledge is lost; it is essential to understanding how animals are responding to our changing world.

Estimation and interpretation problems and solutions when using proportion covariates in linear regression models

Proportion variables, also known as compositional data, are very common in ecology. Unfortunately, few scientists are aware of how compositional data, when used as covariates, can adversely impact statistical analysis. We describe here how proportion covariates result in multicollinearity and parameter identifiability problems. Using simulated data on bird species richness as a function of land use, we show how these problems manifest when fitting a wide range of models in R, both in a frequentist and Bayesian framework. In particular, we show that similar models can often generate substantially different parameter estimates, leading to very different conclusions. Dropping a covariate or the intercept from the model can solve the multicollinearity and parameter identifiability problems. Unfortunately, these solutions do not fix the inherent challenges associated with interpreting parameter estimates. To this end, we propose focusing the interpretation on the difference of slope parameters to avoid the inherent unidentifiability of individual parameters. We also propose conditional plots with two x-axes and marginal plots as visualization techniques that can help users better interpret their modeling results. We illustrate these problems and proposed solutions using empirical data from the North American Breeding Bird Survey. The practical and straightforward approaches suggested in this article will help the fitting of linear models and interpretation of its results when some of the covariates are proportions.

Evaluating MedDRA-to-ICD terminology mappings

BACKGROUND

In this era of big data, data harmonization is an important step to ensure reproducible, scalable, and collaborative research. Thus, terminology mapping is a necessary step to harmonize heterogeneous data. Take the Medical Dictionary for Regulatory Activities (MedDRA) and International Classification of Diseases (ICD) for example, the mapping between them is essential for drug safety and pharmacovigilance research. Our main objective is to provide a quantitative and qualitative analysis of the mapping status between MedDRA and ICD. We focus on evaluating the current mapping status between MedDRA and ICD through the Unified Medical Language System (UMLS) and Observational Medical Outcomes Partnership Common Data Model (OMOP CDM). We summarized the current mapping statistics and evaluated the quality of the current MedDRA-ICD mapping; for unmapped terms, we used our self-developed algorithm to rank the best possible mapping candidates for additional mapping coverage.

RESULTS

The identified MedDRA-ICD mapped pairs cover 27.23% of the overall MedDRA preferred terms (PT). The systematic quality analysis demonstrated that, among the mapped pairs provided by UMLS, only 51.44% are considered an exact match. For the 2400 sampled unmapped terms, 56 of the 2400 MedDRA Preferred Terms (PT) could have exact match terms from ICD.

CONCLUSION

Some of the mapped pairs between MedDRA and ICD are not exact matches due to differences in granularity and focus. For 72% of the unmapped PT terms, the identified exact match pairs illustrate the possibility of identifying additional mapped pairs. Referring to its own mapping standard, some of the unmapped terms should qualify for the expansion of MedDRA to ICD mapping in UMLS.

The method matters. A comparative study of biologging and camera traps as data sources with which to describe wildlife habitat selection

Habitat use is a virtually universal activity among animals and is highly relevant as regards designing wildlife management and conservation actions. This has led to the development of a great variety of methods to study it, of which resource selection functions combined with biologging-derived data (RSF) is the most widely used for this purpose. However this approach has some constraints, such as its invasiveness and high costs. Analytical approaches taking into consideration imperfect detection coupled with camera trap data (IDM) have, therefore, emerged as a non-invasive cost-effective alternative. However, despite the fact that both approaches (RSF and IDM) have been used in habitat selection studies, they should also be comparatively assessed. The objective of this work is consequently to assess them from two perspectives: explanatory and predictive. This has been done by analyzing data obtained from camera traps (60 sampling sites) and biologging (17 animals monitored: 7 red deer Cervus elaphus, 6 fallow deer Dama dama and 4 wild boar Sus scrofa) in the same periods using IDM and RSF, respectively, in Doñana National Park (southern Spain) in order to explain and predict habitat use patterns for three studied species. Our results showed discrepancies between the two approaches, as they identified different predictors as being the most relevant to determine species intensity of use, and they predicted spatial patterns of habitat use with a contrasted level of concordance, depending on species and scale. Given these results and the characteristics of each approach, we suggested that although partly comparable interpretations can be obtained with both approaches, they are not equivalent but rather complementary. The combination of data from biologging and camera traps would, therefore, appear to be suitable for the development of an analytical framework with which to describe and characterise the habitat use processes of wildlife.

An Integrated Multi-Omics and Artificial Intelligence Framework for Advance Plant Phenotyping in Horticulture

This review discusses the transformative potential of integrating multi-omics data and artificial intelligence (AI) in advancing horticultural research, specifically plant phenotyping. The traditional methods of plant phenotyping, while valuable, are limited in their ability to capture the complexity of plant biology. The advent of (meta-)genomics, (meta-)transcriptomics, proteomics, and metabolomics has provided an opportunity for a more comprehensive analysis. AI and machine learning (ML) techniques can effectively handle the complexity and volume of multi-omics data, providing meaningful interpretations and predictions. Reflecting the multidisciplinary nature of this area of research, in this review, readers will find a collection of state-of-the-art solutions that are key to the integration of multi-omics data and AI for phenotyping experiments in horticulture, including experimental design considerations with several technical and non-technical challenges, which are discussed along with potential solutions. The future prospects of this integration include precision horticulture, predictive breeding, improved disease and stress response management, sustainable crop management, and exploration of plant biodiversity. The integration of multi-omics and AI holds immense promise for revolutionizing horticultural research and applications, heralding a new era in plant phenotyping.

Combining passive acoustic data from a towed hydrophone array with visual line transect data to estimate abundance and availability bias of sperm whales (Physeter macrocephalus)

Visual line transect (VLT) surveys are central to the monitoring and study of marine mammals. However, for cryptic species such as deep diving cetaceans VLT surveys alone suffer from problems of low sample sizes and availability bias where animals below the surface are not available to be detected. The advent of passive acoustic monitoring (PAM) technology offers important opportunities to observe deep diving cetaceans but statistical challenges remain particularly when trying to integrate VLT and PAM data. Herein, we present a general framework to combine these data streams to estimate abundance when both surveys are conducted simultaneously. Secondarily, our approach can also be used to derive an estimate of availability bias. We outline three methods that vary in complexity and data requirements which are (1) a simple distance sampling (DS) method that treats the two datasets independently (DS-DS Method), (2) a fully integrated approach that applies a capture-mark recapture (CMR) analysis to the PAM data (CMR-DS Method) and (3) a hybrid approach that requires only a subset of the PAM CMR data (Hybrid Method). To evaluate their performance, we use simulations based on known diving and vocalizing behavior of sperm whales (Physeter macrocephalus). As a case study, we applied the Hybrid Method to data from a shipboard survey of sperm whales and compared estimates to a VLT only analysis. Simulation results demonstrated that the CMR-DS Method and Hybrid Method reduced bias by >90% for both abundance and availability bias in comparison to the simpler DS -DS Method. Overall, the CMR-DS Method was the least biased and most precise. For the case study, our application of the Hybrid Method to the sperm whale dataset produced estimates of abundance and availability bias that were comparable to estimates from the VLT only analysis but with considerably higher precision. Integrating multiple sources of data is an important goal with clear benefits. As a step towards that goal we have developed a novel framework. Results from this study are promising although challenges still remain. Future work may focus on applying this method to other deep-diving species and comparing the proposed method to other statistical approaches that aim to combine information from multiple data sources.

Modelling environmental DNA transport in rivers reveals highly resolved spatio-temporal biodiversity patterns

The ever-increasing threats to riverine ecosystems call for novel approaches for highly resolved biodiversity assessments across taxonomic groups and spatio-temporal scales. Recent advances in the joint use of environmental DNA (eDNA) data and eDNA transport models in rivers (e.g., eDITH) allow uncovering the full structure of riverine biodiversity, hence elucidating ecosystem processes and supporting conservation measures. We applied eDITH to a metabarcoding dataset covering three taxonomic groups (fish, invertebrates, bacteria) and three seasons for a catchment sampled for eDNA at 73 sites. We upscaled eDNA-based biodiversity predictions to approximately 1900 reaches, and assessed α- and β-diversity patterns across seasons and taxonomic groups over the whole network. Genus richness predicted by eDITH was generally higher than values from direct eDNA analysis. Both predicted α- and β-diversity varied depending on season and taxonomic group. Predicted fish α-diversity increased downstream in all seasons, while invertebrate and bacteria α-diversity either decreased downstream or were unrelated to network position. Spatial β-diversity mostly decreased downstream, especially for bacteria. The eDITH model yielded a more refined assessment of freshwater biodiversity as compared to raw eDNA data, both in terms of spatial coverage, diversity patterns and effect of covariates, thus providing a more complete picture of freshwater biodiversity.

Data integration reveals dynamic and systematic patterns of breeding habitat use by a threatened shorebird

Incorporating species distributions into conservation planning has traditionally involved long-term representations of habitat use where temporal variation is averaged to reveal habitats that are most suitable across time. Advances in remote sensing and analytical tools have allowed for the integration of dynamic processes into species distribution modeling. Our objective was to develop a spatiotemporal model of breeding habitat use for a federally threatened shorebird (piping plover, Charadrius melodus). Piping plovers are an ideal candidate species for dynamic habitat models because they depend on habitat created and maintained by variable hydrological processes and disturbance. We integrated a 20-year (2000-2019) nesting dataset with volunteer-collected sightings (eBird) using point process modeling. Our analysis incorporated spatiotemporal autocorrelation, differential observation processes within data streams, and dynamic environmental covariates. We evaluated the transferability of this model in space and time and the contribution of the eBird dataset. eBird data provided more complete spatial coverage in our study system than nest monitoring data. Patterns of observed breeding density depended on both dynamic (e.g., surface water levels) and long-term (e.g., proximity to permanent wetland basins) environmental processes. Our study provides a framework for quantifying dynamic spatiotemporal patterns of breeding density. This assessment can be iteratively updated with additional data to improve conservation and management efforts, because reducing temporal variability to average patterns of use may cause a loss in precision for such actions.

Are we telling the same story? Comparing inferences made from camera trap and telemetry data for wildlife monitoring

Estimating habitat and spatial associations for wildlife is common across ecological studies and it is well known that individual traits can drive population dynamics and vice versa. Thus, it is commonly assumed that individual- and population-level data should represent the same underlying processes, but few studies have directly compared contemporaneous data representing these different perspectives. We evaluated the circumstances under which data collected from Lagrangian (individual-level) and Eulerian (population-level) perspectives could yield comparable inference to understand how scalable information is from the individual to the population. We used Global Positioning System (GPS) collar (Lagrangian) and camera trap (Eulerian) data for seven species collected simultaneously in eastern Washington (2018-2020) to compare inferences made from different survey perspectives. We fit the respective data streams to resource selection functions (RSFs) and occupancy models and compared estimated habitat- and space-use patterns for each species. Although previous studies have considered whether individual- and population-level data generated comparable information, ours is the first to make this comparison for multiple species simultaneously and to specifically ask whether inferences from the two perspectives differed depending on the focal species. We found general agreement between the predicted spatial distributions for most paired analyses, although specific habitat relationships differed. We hypothesize the discrepancies arose due to differences in statistical power associated with camera and GPS-collar sampling, as well as spatial mismatches in the data. Our research suggests data collected from individual-based sampling methods can capture coarse population-wide patterns for a diversity of species, but results differ when interpreting specific wildlife-habitat relationships.

Forecasting the future of life in Antarctica

Antarctic ecosystems are under increasing anthropogenic pressure, but efforts to predict the responses of Antarctic biodiversity to environmental change are hindered by considerable data challenges. Here, we illustrate how novel data capture technologies provide exciting opportunities to sample Antarctic biodiversity at wider spatiotemporal scales. Data integration frameworks, such as point process and hierarchical models, can mitigate weaknesses in individual data sets, improving confidence in their predictions. Increasing process knowledge in models is imperative to achieving improved forecasts of Antarctic biodiversity, which can be attained for data-limited species using hybrid modelling frameworks. Leveraging these state-of-the-art tools will help to overcome many of the data scarcity challenges presented by the remoteness of Antarctica, enabling more robust forecasts both near- and long-term.

A dataset of bird distributions in zoogeographical regions of China

Background

China, the largest country in Asia, has a land area of approximately 9.6 million square kilometres. There are 1481 bird species (following the taxonomy of IOC World Bird List version 12.1) recorded in two zoogeographical realms, seven regions and 19 subregions in the country. From 1955 to 2017, six authoritative monographs were published, which recorded the distribution area for all bird species in China and were widely quoted by research papers and field guides. This massive amount of data could be used to address many hot topics in ornithology, biogeography and ecology. However, rapid changes in the taxonomic status and different schemes of zoogeographical regionalisation in these six monographs provided limits to the utilisation of these valuable data.

New information

By integrating the data from the six monographs, we presented an open-access dataset on the occurrences and residence types of all Chinese bird species in zoogeographical regions over the past 60 years. The taxonomic statuses for these species were determined following the IOC World Bird List version 12.1 and the zoogeographical regions were based on the updated scheme. These data provide valuable information for the research in bird ecology and conservation biology.

Joint analysis of structured and semi-structured community science data improves precision of relative abundance but not trends in birds

Estimating absolute and relative abundance of wildlife populations is critical to addressing ecological questions and conservation needs, yet obtaining reliable estimates can be challenging because surveys are often limited spatially or temporally. Community science (i.e., citizen science) provides opportunities for semi-structured data collected by the public (e.g., eBird) to improve capacity of relative abundance estimation by complementing structured survey data collected by trained observers (e.g., North American breeding bird survey [BBS]). We developed two state-space models to estimate relative abundance and population trends: one using BBS data and the other jointly analyzing BBS and eBird data. We applied these models to seven bird species with diverse life history characteristics. Joint analysis of eBird and BBS data improved precision of mean and year-specific relative abundance estimates for all species, but the BBS-only model produced more precise trend estimates compared to the joint model for most species. The relative abundance estimates of the joint model were particularly more precise than the BBS-only estimates in areas where species detectability was low resulting from either low BBS survey effort or low abundance. These results suggest that community science data can be a valuable resource for cost-effective improvement in wildlife abundance estimation.

Integration of mark-recapture and acoustic detections for unbiased population estimation in animal communities

Abundance estimation methods that combine several types of data are becoming increasingly common because they yield more accurate and precise parameter estimates and predictions than are possible from a single data source. These beneficial effects result from increasing sample size (through data pooling) and complementarity between different data types. Here, we test whether integrating mark-recapture data with passive acoustic detections into a joint likelihood improves estimates of population size in a multi-guild community. We compared the integrated model to a mark-recapture-only model using simulated data first and then using a data set of mist-net captures and acoustic recordings from an Afrotropical agroforest bird community. The integrated model with simulated data improved accuracy and precision of estimated population size and detection parameters. When applied to field data, the integrated model was able to produce, for each bird guild, ecologically plausible estimates of population size and detection parameters, with more precision compared with the mark-recapture model. Overall, our results show that adding acoustic data to mark-recapture analyses improves estimates of population size. With the increasing availability of acoustic recording devices, this data collection technique could readily be added to routine field protocols, leading to a cost-efficient improvement of traditional mark-recapture population estimation.

Integrating data types to estimate spatial patterns of avian migration across the Western Hemisphere

For many avian species, spatial migration patterns remain largely undescribed, especially across hemispheric extents. Recent advancements in tracking technologies and high-resolution species distribution models (i.e., eBird Status and Trends products) provide new insights into migratory bird movements and offer a promising opportunity for integrating independent data sources to describe avian migration. Here, we present a three-stage modeling framework for estimating spatial patterns of avian migration. First, we integrate tracking and band re-encounter data to quantify migratory connectivity, defined as the relative proportions of individuals migrating between breeding and nonbreeding regions. Next, we use estimated connectivity proportions along with eBird occurrence probabilities to produce probabilistic least-cost path (LCP) indices. In a final step, we use generalized additive mixed models (GAMMs) both to evaluate the ability of LCP indices to accurately predict (i.e., as a covariate) observed locations derived from tracking and band re-encounter data sets versus pseudo-absence locations during migratory periods and to create a fully integrated (i.e., eBird occurrence, LCP, and tracking/band re-encounter data) spatial prediction index for mapping species-specific seasonal migrations. To illustrate this approach, we apply this framework to describe seasonal migrations of 12 bird species across the Western Hemisphere during pre- and postbreeding migratory periods (i.e., spring and fall, respectively). We found that including LCP indices with eBird occurrence in GAMMs generally improved the ability to accurately predict observed migratory locations compared to models with eBird occurrence alone. Using three performance metrics, the eBird + LCP model demonstrated equivalent or superior fit relative to the eBird-only model for 22 of 24 species-season GAMMs. In particular, the integrated index filled in spatial gaps for species with over-water movements and those that migrated over land where there were few eBird sightings and, thus, low predictive ability of eBird occurrence probabilities (e.g., Amazonian rainforest in South America). This methodology of combining individual-based seasonal movement data with temporally dynamic species distribution models provides a comprehensive approach to integrating multiple data types to describe broad-scale spatial patterns of animal movement. Further development and customization of this approach will continue to advance knowledge about the full annual cycle and conservation of migratory birds.

Dynamic Species Distribution Modeling Reveals the Pivotal Role of Human-Mediated Long-Distance Dispersal in Plant Invasion

Plant invasions generate massive ecological and economic costs worldwide. Predicting their spatial dynamics is crucial to the design of effective management strategies and the prevention of invasions. Earlier studies highlighted the crucial role of long-distance dispersal in explaining the speed of many invasions. In addition, invasion speed depends highly on the duration of its lag phase, which may depend on the scaling of fecundity with age, especially for woody plants, even though empirical proof is still rare. Bayesian dynamic species distribution models enable the fitting of process-based models to partial and heterogeneous observations using a state-space modeling approach, thus offering a tool to test such hypotheses on past invasions over large spatial scales. We use such a model to explore the roles of long-distance dispersal and age-structured fecundity in the transient invasion dynamics of Plectranthus barbatus, a woody plant invader in South Africa. Our lattice-based model accounts for both short and human-mediated long-distance dispersal, as well as age-structured fecundity. We fitted our model on opportunistic occurrences, accounting for the spatio-temporal variations of the sampling effort and the variable detection rates across datasets. The Bayesian framework enables us to integrate a priori knowledge on demographic parameters and control identifiability issues. The model revealed a massive wave of spatial spread driven by human-mediated long-distance dispersal during the first decade and a subsequent drastic population growth, leading to a global equilibrium in the mid-1990s. Without long-distance dispersal, the maximum population would have been equivalent to 30% of the current equilibrium population. We further identified the reproductive maturity at three years old, which contributed to the lag phase before the final wave of population growth. Our results highlighted the importance of the early eradication of weedy horticultural alien plants around urban areas to hamper and delay the invasive spread.

Consistent trait-temperature interactions drive butterfly phenology in both incidental and survey data

Data availability limits phenological research at broad temporal and spatial extents. Butterflies are among the few taxa with broad-scale occurrence data, from both incidental reports and formal surveys. Incidental reports have biases that are challenging to address, but structured surveys are often limited seasonally and may not span full flight phenologies. Thus, how these data source compare in phenological analyses is unclear. We modeled butterfly phenology in relation to traits and climate using parallel analyses of incidental and survey data, to explore their shared utility and potential for analytical integration. One workflow aggregated “Pollard” surveys, where sites are visited multiple times per year; the other aggregated incidental data from online portals: iNaturalist and eButterfly. For 40 species, we estimated early (10%) and mid (50%) flight period metrics, and compared the spatiotemporal patterns and drivers of phenology across species and between datasets. For both datasets, inter-annual variability was best explained by temperature, and seasonal emergence was earlier for resident species overwintering at more advanced stages. Other traits related to habitat, feeding, dispersal, and voltinism had mixed or no impacts. Our results suggest that data integration can improve phenological research, and leveraging traits may predict phenology in poorly studied species.

Accounting for niche truncation to improve spatial and temporal predictions of species distributions

Species Distribution Models (SDMs) are essential tools for predicting climate change impact on species’ distributions and are commonly employed as an informative tool on which to base management and conservation actions. Focusing only on a part of the entire distribution of a species for fitting SDMs is a common approach. Yet, geographically restricting their range can result in considering only a subset of the species’ ecological niche (i.e., niche truncation) which could lead to biased spatial predictions of future climate change effects, particularly if future conditions belong to those parts of the species ecological niche that have been excluded for model fitting. The integration of large-scale distribution data encompassing the whole species range with more regional data can improve future predictions but comes along with challenges owing to the broader scale and/or lower quality usually associated with these data. Here, we compare future predictions obtained from a traditional SDM fitted on a regional dataset (Switzerland) to predictions obtained from data integration methods that combine regional and European datasets for several bird species breeding in Switzerland. Three models were fitted: a traditional SDM based only on regional data and thus not accounting for niche truncation, a data pooling model where the two datasets are merged without considering differences in extent or resolution, and a downscaling hierarchical approach that accounts for differences in extent and resolution. Results show that the traditional model leads to much larger predicted range changes (either positively or negatively) under climate change than both data integration methods. The traditional model also identified different variables as main drivers of species’ distribution compared to data-integration models. Differences between models regarding predicted range changes were larger for species where future conditions were outside the range of conditions existing in the regional dataset (i.e., when future conditions implied extrapolation). In conclusion, we showed that (i) models calibrated on a geographically restricted dataset provide markedly different predictions than data integration models and (ii) that these differences are at least partly explained by niche truncation. This suggests that using data integration methods could lead to more accurate predictions and more nuanced range changes than regional SDMs through a better characterization of species’ entire realized niches.

Conceptualizing the 3D niche and vertical space use

Spatial partitioning in ecological communities has predominantly been described in two dimensions, yet habitat is complex and 3D. Complex space use mediates community structure and interaction strength by expanding spatial, temporal, and dietary dimensions. Vertical stratification of resources provides opportunities for novel specializations, creating a 3D niche. Competition and predation are mediated by 3D space use, as individuals use the vertical axis to access prey, flee predators, or avoid competitors. The 3D niche is important for long-term conservation strategies as species must navigate tradeoffs in habitat use between strata-specific threats and suboptimal habitat patches. Ultimately, elucidating the 3D niche has implications for protected area management and corridor design that directly influence species persistence and ecosystem function in a rapidly changing world.

Accounting for detection probability with overestimation by integrating double monitoring programs over 40 years

In the context of wildlife population declines, increasing computer power over the last 20 years allowed wildlife managers to apply advanced statistical techniques that has improved population size estimates. However, respecting the assumptions of the models that consider the probability of detection, such as N-mixture models, requires the implementation of a rigorous monitoring protocol with several replicate survey occasions and no double counting that are hardly adaptable to field conditions. When the logistical, economic and ecological constraints are too strong to meet model assumptions, it may be possible to combine data from independent surveys into the modelling framework in order to understand population dynamics more reliably. Here, we present a state-space model with an error process modelled on the log scale to evaluate wintering waterfowl numbers in the Camargue, southern France, while taking a conditional probability of detection into consideration. Conditional probability of detection corresponds to estimation of a detection probability index, which is not a true probability of detection, but rather conditional on the difference to a particular baseline. The large number of sites (wetlands within the Camargue delta) and years monitored (44) provide significant information to combine both terrestrial and aerial surveys (which constituted spatially and temporally replicated counts) to estimate a conditional probability of detection, while accounting for false-positive counting errors and changes in observers over the study period. The model estimates abundance indices of wintering Common Teal, Mallard and Common Coot, all species abundant in the area. We found that raw counts were underestimated compared to the predicted population size. The model-based data integration approach as described here seems like a promising solution that takes advantage of as much as possible of the data collected from several methods when the logistic constraints do not allow the implementation of a permanent monitoring and analysis protocol that takes into account the detectability of individuals.

Finding rare species and estimating the probability that all occupied sites have been found

Detecting occupied sites of rare species, and estimating the probability that all occupied sites are known within a given area, are desired outcomes for many ecological or conservation projects. Examples include managing all occupied sites of a threatened species or eradicating an emerging invader. Occupied sites may remain undetected because (1) sites where the species potentially occurs had not been searched, and (2) the species could have been overlooked in the searched sites. For rare species, available data are typically scant, making it difficult to predict sites where the species probably occurs or to estimate detection probability in the searched sites. Using the critically endangered Rose's mountain toadlet (Capensibufo rosei), known from only two localities, we outline an iterative process aimed at estimating the probability that any unknown occupied sites remain and maximizing the chance of finding them. This includes fitting a species distribution model to guide sampling effort, testing model accuracy and sampling efficacy using the occurrence of more common proxy species, and estimating detection probability using sites of known presence. The final estimate of the probability that all occupied sites were found incorporates the uncertainties of uneven distribution, relative area searched, and detection probability. Our results show that very few occupied sites of C. rosei are likely to remain undetected. We also show that the probability of an undetected occupied site remaining will always be high for large unsearched areas of potential occurrence, but can be low for smaller areas intended for targeted management interventions. Our approach is especially useful for assessing uncertainty in species occurrences, planning the required search effort needed to reduce probability of unknown occurrence to desired levels, and identifying priority areas for further searches or management interventions.

Estimating the movements of terrestrial animal populations using broad-scale occurrence data

As human and automated sensor networks collect increasingly massive volumes of animal observations, new opportunities have arisen to use these data to infer or track species movements. Sources of broad scale occurrence datasets include crowdsourced databases, such as eBird and iNaturalist, weather surveillance radars, and passive automated sensors including acoustic monitoring units and camera trap networks. Such data resources represent static observations, typically at the species level, at a given location. Nonetheless, by combining multiple observations across many locations and times it is possible to infer spatially continuous population-level movements. Population-level movement characterizes the aggregated movement of individuals comprising a population, such as range contractions, expansions, climate tracking, or migration, that can result from physical, behavioral, or demographic processes. A desire to model population movements from such forms of occurrence data has led to an evolving field that has created new analytical and statistical approaches that can account for spatial and temporal sampling bias in the observations. The insights generated from the growth of population-level movement research can complement the insights from focal tracking studies, and elucidate mechanisms driving changes in population distributions at potentially larger spatial and temporal scales. This review will summarize current broad-scale occurrence datasets, discuss the latest approaches for utilizing them in population-level movement analyses, and highlight studies where such analyses have provided ecological insights. We outline the conceptual approaches and common methodological steps to infer movements from spatially distributed occurrence data that currently exist for terrestrial animals, though similar approaches may be applicable to plants, freshwater, or marine organisms.

Bridging the gap between commercial fisheries and survey data to model the spatiotemporal dynamics of marine species

Monitoring and assessment of natural resources often require inputs from multiple data sources. In fisheries science, for example, the inference of a species' abundance distribution relies on two main data sources, namely commercial fisheries and scientific survey data. Despite efforts to combine these data into an integrated statistical model, their coupling is frequently hampered due to differences in their sampling designs, which imposes distinct bias sources in the estimator of the abundance distribution. We developed a flexible species distribution model (SDM) that can integrate both data sources while filtering out their relative bias contributions. We applied the model on three different age groups of the western Baltic cod stock. For each age group, we tested the model on (1) survey data and (2) integrated data (survey + commercial) as a means to compare their differences and investigate how the inclusion of commercial fisheries data improved the spatiotemporal abundance estimator and parameter estimates. Moreover, we proposed a novel validation approach to evaluate whether the inclusion of commercial fisheries data in the integrated model is not in direct contradiction with the survey data. Following our approach, the results indicated that the use of commercial fisheries data is suitable for the integrated model. Across all age groups, our results demonstrated how commercial fisheries supplied additional information on cod's spatiotemporal abundance dynamics, highlighting sometimes abundance hot spots that were not detected by the survey model alone. Additionally, the integrated model provided a reduction of up to 20% and 10% in the uncertainty (SE) of the predicted abundance fields and fixed-effect parameters, respectively. The proposed model represents thus a valuable benchmark for evaluating spatiotemporal dynamics of fish, and strengthens the science-based advice for marine policymakers.

Open Data Practices among Users of Primary Biodiversity Data

Presence-only biodiversity data are increasingly relied on in biodiversity, ecology, and conservation research, driven by growing digital infrastructures that support open data sharing and reuse. Recent reviews of open biodiversity data have clearly documented the value of data sharing, but the extent to which the biodiversity research community has adopted open data practices remains unclear. We address this question by reviewing applications of presence-only primary biodiversity data, drawn from a variety of sources beyond open databases, in the indexed literature. We characterize how frequently researchers access open data relative to data from other sources, how often they share newly generated or collated data, and trends in metadata documentation and data citation. Our results indicate that biodiversity research commonly relies on presence-only data that are not openly available and neglects to make such data available. Improved data sharing and documentation will increase the value, reusability, and reproducibility of biodiversity research.

Estimating viral prevalence with data fusion for adaptive two-phase pooled sampling

The COVID-19 pandemic has highlighted the importance of efficient sampling strategies and statistical methods for monitoring infection prevalence, both in humans and in reservoir hosts. Pooled testing can be an efficient tool for learning pathogen prevalence in a population. Typically, pooled testing requires a second-phase retesting procedure to identify infected individuals, but when the goal is solely to learn prevalence in a population, such as a reservoir host, there are more efficient methods for allocating the second-phase samples.To estimate pathogen prevalence in a population, this manuscript presents an approach for data fusion with two-phased testing of pooled samples that allows more efficient estimation of prevalence with less samples than traditional methods. The first phase uses pooled samples to estimate the population prevalence and inform efficient strategies for the second phase. To combine information from both phases, we introduce a Bayesian data fusion procedure that combines pooled samples with individual samples for joint inferences about the population prevalence.Data fusion procedures result in more efficient estimation of prevalence than traditional procedures that only use individual samples or a single phase of pooled sampling.The manuscript presents guidance on implementing the first-phase and second-phase sampling plans using data fusion. Such methods can be used to assess the risk of pathogen spillover from reservoir hosts to humans, or to track pathogens such as SARS-CoV-2 in populations.

Long-term trends in the occupancy of ants revealed through use of multi-sourced datasets

We combined participatory science data and museum records to understand long-term changes in occupancy for 29 ant species in Denmark over 119 years. Bayesian occupancy modelling indicated change in occupancy for 15 species: five increased, four declined and six showed fluctuating trends. We consider how trends may have been influenced by life-history and habitat changes. Our results build on an emerging picture that biodiversity change in insects is more complex than implied by the simple insect decline narrative.

Data integration methods to account for spatial niche truncation effects in regional projections of species distribution

Many species distribution models (SDMs) are built with precise but geographically restricted presence-absence data sets (e.g., a country) where only a subset of the environmental conditions experienced by a species across its range is considered (i.e., spatial niche truncation). This type of truncation is worrisome because it can lead to incorrect predictions e.g., when projecting to future climatic conditions belonging to the species niche but unavailable in the calibration area. Data from citizen-science programs, species range maps or atlases covering the full species range can be used to capture those parts of the species' niche that are missing regionally. However, these data usually are too coarse or too biased to support regional management. Here, we aim to (1) demonstrate how varying degrees of spatial niche truncation affect SDMs projections when calibrated with climatically truncated regional data sets and (2) test the performance of different methods to harness information from larger-scale data sets presenting different spatial resolutions to solve the spatial niche truncation problem. We used simulations to compare the performance of the different methods, and applied them to a real data set to predict the future distribution of a plant species (Potentilla aurea) in Switzerland. SDMs calibrated with geographically restricted data sets expectedly provided biased predictions when projected outside the calibration area or time period. Approaches integrating information from larger-scale data sets using hierarchical data integration methods usually reduced this bias. However, their performance varied depending on the level of spatial niche truncation and how data were combined. Interestingly, while some methods (e.g., data pooling, downscaling) performed well on both simulated and real data, others (e.g., those based on a Poisson point process) performed better on real data, indicating a dependency of model performance on the simulation process (e.g., shape of simulated response curves). Based on our results, we recommend to use different data integration methods and, whenever possible, to make a choice depending on model performance. In any case, an ensemble modeling approach can be used to account for uncertainty in how niche truncation is accounted for and identify areas where similarities/dissimilarities exist across methods.

Using single visits into integrated occupancy models to make the most of existing monitoring programs

A major challenge in statistical ecology consists of integrating knowledge from different data sets to produce robust ecological indicators. To estimate species distribution, occupancy models are a flexible framework that can accommodate several data sets obtained from different sampling methods. However, repeating visits at sampling sites is a prerequisite for using standard occupancy models. Occupancy models were recently developed to analyze detection/non-detection data collected during a single visit. To date, single-visit occupancy models have never been used to integrate several different data sets. Here, we showcase an approach that combines two data sets into an integrated single-visit occupancy model. As a case study, we estimated the distribution of common bottlenose dolphin (Tursiops truncatus) over the northwestern Mediterranean Sea by combining 24,624 km of aerial surveys and 21,464 km of at-sea monitoring. We compared the outputs of single- vs. repeated-visit occupancy models into integrated occupancy models. Integrated models allowed a better sampling coverage of the targeted population, which provided a better precision for occupancy estimates than occupancy models using data sets in isolation. Overall, single- and repeated-visit integrated occupancy models produced similar inference about the distribution of bottlenose dolphins. We suggest that single-visit occupancy models open promising perspectives for the use of existing ecological data sets.

Trends in bird abundance differ among protected forests but not bird guilds

Improved monitoring and associated inferential tools to efficiently identify declining bird populations, particularly of rare or sparsely distributed species, is key to informed conservation and management across large spatiotemporal regions. We assess abundance trends for 106 bird species in a network of eight forested national parks located within the northeast United States from 2006 to 2019 using a novel hierarchical model. We develop a multispecies, multiregion, removal-sampling model that shares information across species and parks to enable inference on rare species and sparsely sampled parks and to evaluate the effects of local forest structure. Trends in bird abundance over time varied widely across parks, but species showed similar trends within parks. Three parks (Acadia National Park and Marsh-Billings-Rockefeller and Morristown National Historical Parks [NHP]) decreased in bird abundance across all species, while three parks (Saratoga NHP and Roosevelt-Vanderbilt and Weir-Farm National Historic Sites) increased in abundance. Bird abundance peaked at medium levels of basal area and high levels of percent forest and forest regeneration, with percent forest having the largest effect. Variation in these effects across parks could be a result of differences in forest structural stage and diversity. By sharing information across both communities and parks, our novel hierarchical model enables uncertainty-quantified estimates of abundance across multiple geographical (i.e., network, park) and taxonomic (i.e., community, guild, species) levels over a large spatiotemporal region. We found large variation in abundance trends across parks but not across bird guilds, suggesting that local forest condition might have a broad and consistent effect on the entire bird community within a given park. Research should target the three parks with overall decreasing trends in bird abundance to further identify what specific factors are driving observed declines across the bird community. Understanding how bird communities respond to local forest structure and other stressors (e.g., pest outbreaks, climate change) is crucial for informed and lasting management.

A Framework for Estimating Human-Wildlife Conflict Probabilities Conditional on Species Occupancy

Managing human-wildlife conflicts (HWCs) is an important conservation objective for the many terrestrial landscapes dominated by humans. Forecasting where future conflicts are likely to occur and assessing risks to lives and livelihoods posed by wildlife are central to informing HWC management strategies. Existing assessments of the spatial occurrence patterns of HWC are based on either understanding spatial patterns of past conflicts or patterns of species distribution. In the former case, the absence of conflicts at a site cannot be attributed to the absence of the species. In the latter case, the presence of a species may not be an accurate measure of the probability of conflict occurrence. We present a Bayesian hierarchical modeling framework that integrates conflict reporting data and species distribution data, thus allowing the estimation of the probability that conflicts with a species are reported from a site, conditional on the species being present. In doing so, our model corrects for both false-positive and false-negative conflict reporting errors. We provide study design recommendations using simulations that explore the performance of the model under a range of conflict reporting probabilities. We applied the model to data on wild boar (Sus scrofa) space use and conflicts collected from the Central Terai Landscape (CTL), an important tiger conservation landscape in India. We found that tolerance for wildlife was a predictor of the probability with which farmers report conflict with wild boars from sites not used by the species. We also discuss useful extensions of the model when conflict data are verified for potential false-positive errors and when landscapes are monitored over multiple seasons.