Population dynamics affect the capacity of species distribution models to predict species abundance on a local scale

Using species distribution modeling to generate relative abundance information in socio-politically unstable territories: Conservation of Felidae in the central-western region of Mexico

The distribution range and population abundance of species provide fundamental information on the species-habitat relationship required for management and conservation. Abundance inherently provides more information about the ecology of species than do occurrence data. However, information on abundance is scarce for most species, mainly at large spatial scales. The objective of this work was, therefore, to provide information regarding the population status of six wild felids inhabiting territories in Mexico that are inaccessible or politically unstable. This was done using species distribution models derived from occurrence data. We used distribution data at a continental scale for the wild felids inhabiting Mexico: jaguar (Panthera onca), bobcat (Lynx rufus), ocelot (Leopardus pardalis), cougar (Puma concolor), margay (Leopardus wiedii), and jaguarundi (Herpailurus yagouaroundi) to predict environmental suitability (estimated by both Maxent and the distance to niche centroid, DNC). Suitability was then examined by relating to a capture rate-based index, in a well-monitored area in central western Mexico in order to assess their performance as proxies of relative abundance. Our results indicate that the environmental suitability patterns predicted by both algorithms were comparable. However, the strength of the relationship between the suitability and relative abundance of local populations differed across species and between algorithms, with the bobcat and DNC, respectively, having the best fit, although the relationship was not consistent in all the models. This paper presents the potential of implementing species distribution models in order to predict the relative abundance of wild felids in Mexico and offers guidance for the proper interpretation of the relationship between suitability and population abundance. The results obtained provide a robust information base on which to outline specific conservation actions and on which to examine the potential status of endangered species inhabiting remote or politically unstable territories in which on-field monitoring programs are not feasible.

Spatiotemporal connectivity dynamics in spatially structured populations

Connectivity is a fundamental concept linking dispersal to the emergent dynamics and persistence of spatially structured populations. Functional measures of connectivity typically seek to integrate aspects of landscape structure and animal movement to describe ecologically meaningful connectedness at the landscape and population scale. Despite this focus on function, traditional measures of landscape connectivity assume it is a static property of the landscape, hence abstracting out the underlying spatiotemporal population dynamics. Connectivity is, arguably, a dynamic property of landscapes, and is inherently related to the spatial distribution of individuals and populations across the landscape. Static representations of connectivity potentially overlook this variation and therefore adopting a dynamic approach should offer improved insights about connectivity and associated ecological processes. Using a large-scale, long-term time series of occupancy data from a metapopulation of water voles Arvicola amphibius, we tested competing hypotheses about how considering the dynamic nature of connectivity improves the ability of spatially explicit occupancy models to recover population dynamics. Iteratively relaxing standing assumptions of connectivity metrics, these models ranged from spatially and temporally fixed connectivity metrics that are widely applied, to the more flexible, but lesser used model that allowed temporally varying connectivity measures that incorporate spatiotemporally dynamic patch occupancy states. Our results provide empirical evidence that demographic weighting using patch occupancy dynamics and temporal variability in connectivity measures are important for describing metapopulation dynamics. We highlight the implications of commonly held assumption in connectivity modelling and demonstrate how they result in different and highly variable predictions of metapopulation capacity. Thus, we argue that the concept of connectivity and its potential applications would benefit from recognizing inherent spatiotemporal variation in connectivity that is explicitly linked to underlying ecological state variables.

Spatial epidemiology of Leptospira sp. exposure in bovines from Veracruz, México

Bovine leptospirosis is a bacterial disease that affects cattle herds, causing economic losses due to reproductive problems which require expensive treatments. The main source of transmission for cattle is still uncertain, but rodents and bats can play an important role in the transmission cycle by being maintenance hosts for the pathogenic species of the bacterium and spreading it through urine. In this study, we characterize possible risk areas for bovine leptospirosis exposure in the state of Veracruz, Mexico, based on the geographical distribution of flying (bats) and terrestrial (rodents and opossums) wild hosts of Leptospira sp. reported in Mexico, in addition to climate, geography, soil characteristics, land use and human activities (environmental variables). We used a generalized linear regression model to understand the association between the frequency of anti-Leptospira sp. antibodies (a proxy of exposure) in cattle herds exposed to Leptospira, the favourability of wild hosts of Leptospira as well as the environmental variables. The parameterized model explained 12.3% of the variance. The frequency of anti-Leptospira sp. antibodies exposure in cattle herds was associated with elevation, geographic longitude, pH of the soil surface and environmental favourability for the presence of rodents, opossums and bats. The variation in exposure was mainly explained by a longitudinal gradient (6.4% of the variance) and the favourability-based indices for wild hosts (9.6% of the variance). Describing the possible risks for exposure to Leptospira in an important and neglected livestock geographical region, we provide valuable information for the selection of areas for diagnosis and prevention of this relevant disease.

Tropical and subtropical Asia's valued tree species under threat

Tree diversity in Asia's tropical and subtropical forests is central to nature-based solutions. Species vulnerability to multiple threats, which affect provision of ecosystem services, is poorly understood. We conducted a region-wide, spatially explicit assessment of the vulnerability of 63 socioeconomically important tree species to overexploitation, fire, overgrazing, habitat conversion, and climate change. Trees were selected for assessment from national priority lists, and selections were validated by an expert network representing 20 countries. We used Maxent suitability modeling to predict species distribution ranges, freely accessible spatial data sets to map threat exposures, and functional traits to estimate threat sensitivities. Species-specific vulnerability maps were created as the product of exposure maps and sensitivity estimates. Based on vulnerability to current threats and climate change, we identified priority areas for conservation and restoration. Overall, 74% of the most important areas for conservation of these trees fell outside protected areas, and all species were severely threatened across an average of 47% of their native ranges. The most imminent threats were overexploitation and habitat conversion; populations were severely threatened by these factors in an average of 24% and 16% of their ranges, respectively. Our model predicted limited overall climate change impacts, although some study species were likely to lose over 15% of their habitat by 2050 due to climate change. We pinpointed specific natural areas in Borneo rain forests as hotspots for in situ conservation of forest genetic resources, more than 82% of which fell outside designated protected areas. We also identified degraded areas in Western Ghats, Indochina dry forests, and Sumatran rain forests as hotspots for restoration, where planting or assisted natural regeneration will help conserve these species, and croplands in southern India and Thailand as potentially important agroforestry options. Our results highlight the need for regionally coordinated action for effective conservation and restoration.

Species distribution models for invasive Eurasian watermilfoil highlight the importance of data quality and limitations of discrimination accuracy metrics

AIM

Availability of uniformly collected presence, absence, and abundance data remains a key challenge in species distribution modeling (SDM). For invasive species, abundance and impacts are highly variable across landscapes, and quality occurrence and abundance data are critical for predicting locations at high risk for invasion and impacts, respectively. We leverage a large aquatic vegetation dataset comprising point-level survey data that includes information on the invasive plant Myriophyllum spicatum (Eurasian watermilfoil) to: (a) develop SDMs to predict invasion and impact from environmental variables based on presence-absence, presence-only, and abundance data, and (b) compare evaluation metrics based on functional and discrimination accuracy for presence-absence and presence-only SDMs.

LOCATION

Minnesota, USA.

METHODS

Eurasian watermilfoil presence-absence and abundance information were gathered from 468 surveyed lakes, and 801 unsurveyed lakes were leveraged as pseudoabsences for presence-only models. A Random Forest algorithm was used to model the distribution and abundance of Eurasian watermilfoil as a function of lake-specific predictors, both with and without a spatial autocovariate. Occurrence-based SDMs were evaluated using conventional discrimination accuracy metrics and functional accuracy metrics assessing correlation between predicted suitability and observed abundance.

RESULTS

Water temperature degree days and maximum lake depth were two leading predictors influencing both invasion risk and abundance, but they were relatively less important for predicting abundance than other water quality measures. Road density was a strong predictor of Eurasian watermilfoil invasion risk but not abundance. Model evaluations highlighted significant differences: Presence-absence models had high functional accuracy despite low discrimination accuracy, whereas presence-only models showed the opposite pattern.

MAIN CONCLUSION

Complementing presence-absence data with abundance information offers a richer understanding of invasive Eurasian watermilfoil's ecological niche and enables evaluation of the model's functional accuracy. Conventional discrimination accuracy measures were misleading when models were developed using pseudoabsences. We thus caution against the overuse of presence-only models and suggest directing more effort toward systematic monitoring programs that yield high-quality data.

Earth and field observations underpin metapopulation dynamics in complex landscapes: Near-term study on carabids

Understanding risks to biodiversity requires predictions of the spatial distribution of species adapting to changing ecosystems and, to that end, Earth observations integrating field surveys prove essential as they provide key numbers for assessing landscape-wide biodiversity scenarios. Here, we develop, and apply to a relevant case study, a method suited to merge Earth/field observations with spatially explicit stochastic metapopulation models to study the near-term ecological dynamics of target species in complex terrains. Our framework incorporates the use of species distribution models for a reasoned estimation of the initial presence of the target species and accounts for imperfect and incomplete detection of the species presence in the study area. It also uses a metapopulation fitness function derived from Earth observation data subsuming the ecological niche of the target species. This framework is applied to contrast occupancy of two species of carabids (Pterostichus flavofemoratus, Carabus depressus) observed in the context of a large ecological monitoring program carried out within the Gran Paradiso National Park (GPNP, Italy). Results suggest that the proposed framework may indeed exploit the hallmarks of spatially explicit ecological approaches and of remote Earth observations. The model reproduces well the observed in situ data. Moreover, it projects in the near term the two species' presence both in space and in time, highlighting the features of the metapopulation dynamics of colonization and extinction, and their expected trends within verifiable timeframes.

Comparing abundance distributions and range maps in spatial conservation planning for migratory species

Most spatial conservation planning for wide-ranging or migratory species is constrained by poor knowledge of species' spatiotemporal dynamics and is only based on static species' ranges. However, species have substantial variation in abundance across their range and migratory species have important spatiotemporal population dynamics. With growing ecological data and advancing analytics, both of these can be estimated and incorporated into spatial conservation planning. However, there is limited information on the degree to which including this information affects conservation planning. We compared the performance of systematic conservation prioritizations for different scenarios based on varying the input species' distributions by ecological metric (abundance distributions versus range maps) and temporal sampling resolution (weekly, monthly, or quarterly). We used the example of a community of 41 species of migratory shorebirds that breed in North America, and we used eBird data to produce weekly estimates of species' abundances and ranges. Abundance distributions at a monthly or weekly resolution led to prioritizations that most efficiently protected species throughout the full annual cycle. Conversely, spatial prioritizations based on species' ranges required more sites and left most species insufficiently protected for at least part of their annual cycle. Prioritizations with only quarterly species ranges were very inefficient as they needed to target 40% of species' ranges to include 10% of populations. We highlight the high value of abundance information for spatial conservation planning, which leads to more efficient and effective spatial prioritization for conservation. Overall, we provide evidence that spatial conservation planning for wide-ranging migratory species is most robust and efficient when informed by species' abundance information from the full annual cycle.

Relationships between population densities and niche-centroid distances in North American birds

Correlational ecological niche models have seen intensive use and exploration as a means of estimating the limits of actual and potential geographic distributions of species, yet their application to explaining geographic abundance patterns has been debated. We developed a detailed test of this latter possibility based on the North American Breeding Bird Survey. Correlations between abundances and niche-centroid distances were mostly negative, as per expectations of niche theory and the abundant niche-centre relationship. The negative relationships were not distributed randomly among species: terrestrial, non-migratory, small-bodied, small-niche-breadth and restricted-range species had the strongest negative associations. Distances to niche centroids as estimated from correlational analyses of presence-only data thus offer a unique means by which to infer geographic abundance patterns, which otherwise are enormously difficult to characterise.

Machine learning of large-scale spatial distributions of wild turkeys with high-dimensional environmental data

Species distribution modeling often involves high-dimensional environmental data. Large amounts of data and multicollinearity among covariates impose challenges to statistical models in variable selection for reliable inferences of the effects of environmental factors on the spatial distribution of species. Few studies have evaluated and compared the performance of multiple machine learning (ML) models in handling multicollinearity. Here, we assessed the effectiveness of removal of correlated covariates and regularization to cope with multicollinearity in ML models for habitat suitability. Three machine learning algorithms maximum entropy (MaxEnt), random forests (RFs), and support vector machines (SVMs) were applied to the original data (OD) of 27 landscape variables, reduced data (RD) with 14 highly correlated covariates being removed, and 15 principal components (PC) of the OD accounting for 90% of the original variability. The performance of the three ML models was measured with the area under the curve and continuous Boyce index. We collected 663 nonduplicated presence locations of Eastern wild turkeys (Meleagris gallopavo silvestris) across the state of Mississippi, United States. Of the total locations, 453 locations separated by a distance of ≥2 km were used to train the three ML algorithms on the OD, RD, and PC data, respectively. The remaining 210 locations were used to validate the trained ML models to measure ML performance. Three ML models had excellent performance on the RD and PC data. MaxEnt and SVMs had good performance on the OD data, indicating the adequacy of regularization of the default setting for multicollinearity. Weak learning of RFs through bagging appeared to alleviate multicollinearity and resulted in excellent performance on the OD data. Regularization of ML algorithms may help exploratory studies of the effects of environmental factors on the spatial distribution and habitat suitability of wildlife.

Environmental suitability models predict population density, performance and body condition for microendemic salamanders

Species can show strong variation of local abundance across their ranges. Recent analyses suggested that variation in abundance can be related to environmental suitability, as the highest abundances are often observed in populations living in the most suitable areas. However, there is limited information on the mechanisms through which variation in environmental suitability determines abundance. We analysed populations of the microendemic salamander Hydromantes flavus, and tested several hypotheses on potential relationships linking environmental suitability to population parameters. For multiple populations across the whole species range, we assessed suitability using species distribution models, and measured density, activity level, food intake and body condition index. In high-suitability sites, the density of salamanders was up to 30-times higher than in the least suitable ones. Variation in activity levels and population performance can explain such variation of abundance. In high-suitability sites, salamanders were active close to the surface, and showed a low frequency of empty stomachs. Furthermore, when taking into account seasonal variation, body condition was better in the most suitable sites. Our results show that the strong relationship between environmental suitability and population abundance can be mediated by the variation of parameters strongly linked to individual performance and fitness.