A Bayesian geostatistical approach to modeling global distributions of Lygodium microphyllum under projected climate warming

Comparison of new computational methods for spatial modelling of malaria

BACKGROUND

Geostatistical analysis of health data is increasingly used to model spatial variation in malaria prevalence, burden, and other metrics. Traditional inference methods for geostatistical modelling are notoriously computationally intensive, motivating the development of newer, approximate methods for geostatistical analysis or, more broadly, computational modelling of spatial processes. The appeal of faster methods is particularly great as the size of the region and number of spatial locations being modelled increases.

METHODS

This work presents an applied comparison of four proposed 'fast' computational methods for spatial modelling and the software provided to implement them-Integrated Nested Laplace Approximation (INLA), tree boosting with Gaussian processes and mixed effect models (GPBoost), Fixed Rank Kriging (FRK) and Spatial Random Forests (SpRF). The four methods are illustrated by estimating malaria prevalence on two different spatial scales-country and continent. The performance of the four methods is compared on these data in terms of accuracy, computation time, and ease of implementation.

RESULTS

Two of these methods-SpRF and GPBoost-do not scale well as the data size increases, and so are likely to be infeasible for larger-scale analysis problems. The two remaining methods-INLA and FRK-do scale well computationally, however the resulting model fits are very sensitive to the user's modelling assumptions and parameter choices. The binomial observation distribution commonly used for disease prevalence mapping with INLA fails to account for small-scale overdispersion present in the malaria prevalence data, which can lead to poor predictions. Selection of an appropriate alternative such as the Beta-binomial distribution is required to produce a reliable model fit. The small-scale random effect term in FRK overcomes this pitfall, but FRK model estimates are very reliant on providing a sufficient number and appropriate configuration of basis functions. Unfortunately the computation time for FRK increases rapidly with increasing basis resolution.

CONCLUSIONS

INLA and FRK both enable scalable geostatistical modelling of malaria prevalence data. However care must be taken when using both methods to assess the fit of the model to data and plausibility of predictions, in order to select appropriate model assumptions and parameters.

Habitat and climatic associations of climate-sensitive species along a southern range boundary

Climate change and habitat loss are recognized as important drivers of shifts in wildlife species' geographic distributions. While often considered independently, there is considerable overlap between these drivers, and understanding how they contribute to range shifts can predict future species assemblages and inform effective management. Our objective was to evaluate the impacts of habitat, climatic, and anthropogenic effects on the distributions of climate-sensitive vertebrates along a southern range boundary in Northern Michigan, USA. We combined multiple sources of occurrence data, including harvest and citizen-science data, then used hierarchical Bayesian spatial models to determine habitat and climatic associations for four climate-sensitive vertebrate species (American marten [Martes americana], snowshoe hare [Lepus americanus], ruffed grouse [Bonasa umbellus] and moose [Alces alces]). We used total basal area of at-risk forest types to represent habitat, and temperature and winter habitat indices to represent climate. Marten associated with upland spruce-fir and lowland riparian forest types, hares with lowland conifer and aspen-birch, grouse with lowland riparian hardwoods, and moose with upland spruce-fir. Species differed in climatic drivers with hares positively associated with cooler annual temperatures, moose with cooler summer temperatures and grouse with colder winter temperatures. Contrary to expectations, temperature variables outperformed winter habitat indices. Model performance varied greatly among species, as did predicted distributions along the southern edge of the Northwoods region. As multiple species were associated with lowland riparian and upland spruce-fir habitats, these results provide potential for efficient prioritization of habitat management. Both direct and indirect effects from climate change are likely to impact the distribution of climate-sensitive species in the future and the use of multiple data types and sources in the modelling of species distributions can result in more accurate predictions resulting in improved management at policy-relevant scales.

Vector Surveillance, Host Species Richness, and Demographic Factors as West Nile Disease Risk Indicators

West Nile virus (WNV) is the most common arthropod-borne virus (arbovirus) in the United States (US) and is the leading cause of viral encephalitis in the country. The virus has affected tens of thousands of US persons total since its 1999 North America introduction, with thousands of new infections reported annually. Approximately 1% of humans infected with WNV acquire neuroinvasive West Nile Disease (WND) with severe encephalitis and risk of death. Research describing WNV ecology is needed to improve public health surveillance, monitoring, and risk assessment. We applied Bayesian joint-spatiotemporal modeling to assess the association of vector surveillance data, host species richness, and a variety of other environmental and socioeconomic disease risk factors with neuroinvasive WND throughout the conterminous US. Our research revealed that an aging human population was the strongest disease indicator, but climatic and vector-host biotic interactions were also significant in determining risk of neuroinvasive WND. Our analysis also identified a geographic region of disproportionately high neuroinvasive WND disease risk that parallels the Continental Divide, and extends southward from the US-Canada border in the states of Montana, North Dakota, and Wisconsin to the US-Mexico border in western Texas. Our results aid in unraveling complex WNV ecology and can be applied to prioritize disease surveillance locations and risk assessment.