Black swans in space: modeling spatiotemporal processes with extremes

Global determinants of insect mitochondrial genetic diversity

Understanding global patterns of genetic diversity is essential for describing, monitoring, and preserving life on Earth. To date, efforts to map macrogenetic patterns have been restricted to vertebrates, which comprise only a small fraction of Earth's biodiversity. Here, we construct a global map of predicted insect mitochondrial genetic diversity from cytochrome c oxidase subunit 1 sequences, derived from open data. We calculate the mitochondrial genetic diversity mean and genetic diversity evenness of insect assemblages across the globe, identify their environmental correlates, and make predictions of mitochondrial genetic diversity levels in unsampled areas based on environmental data. Using a large single-locus genetic dataset of over 2 million globally distributed and georeferenced mtDNA sequences, we find that mitochondrial genetic diversity evenness follows a quadratic latitudinal gradient peaking in the subtropics. Both mitochondrial genetic diversity mean and evenness positively correlate with seasonally hot temperatures, as well as climate stability since the last glacial maximum. Our models explain 27.9% and 24.0% of the observed variation in mitochondrial genetic diversity mean and evenness in insects, respectively, making an important step towards understanding global biodiversity patterns in the most diverse animal taxon.

Estimating phenology and phenological shifts with hierarchical modeling

Climate-driven changes to phenology are some of the most prevalent climate change impacts, yet there is no commonly accepted approach to modeling phenological shifts. Here, we present a hierarchical modeling framework for estimating intra-annual patterns in phenology (e.g., peak phenological expression) and analyzing interannual rates of change in peak phenology. Our approach allows for the estimation of multiple sources of uncertainty, including observation error (e.g., imperfect observations of intra-annual patterns in phenology like peak flowering date) and variation in phenological processes (e.g., uncertainty in the rate of change in annual peak phenological expression). Covariates may be included as predictors of annual peaks or interannual variability in phenological responses. We demonstrate the use of our hierarchical modeling framework in two migratory species-juvenile chum salmon and Swainson's thrush. We acknowledge that the complexity of hierarchical models can be difficult to implement from scratch and present an R package that can be used to model peak dates and range (number of days between 25th- and 75th-quartile dates), as well as a rate of change in peak phenology. Increasing precision, calculating uncertainty, and allowing for imperfect data sets when estimating phenological shifts should help ecologists understand how organisms respond to climate change.

Spatial correlation reverses the compound effect of multiple stressors on rocky shore biofilm

Understanding how multifactorial fluctuating environments affect species and communities remains one of the major challenges in ecology. The spatial configuration of the environment is known to generate complex patterns of correlation among multiple stressors. However, to what extent the spatial correlation between simultaneously fluctuating variables affects ecological assemblages in real-world conditions remains poorly understood. Here, we use field experiments and simulations to assess the influence of spatial correlation of two relevant climate variables - warming and sediment deposition following heavy precipitation - on the biomass and photosynthetic activity of rocky intertidal biofilm. First, we used a response-surface design experiment to establish the relation between biofilm, warming, and sediment deposition in the field. Second, we used the response surface to generate predictions of biofilm performance under different scenarios of warming and sediment correlation. Finally, we tested the predicted outcomes by manipulating the degree of correlation between the two climate variables in a second field experiment. Simulations stemming from the experimentally derived response surface showed how the degree and direction (positive or negative) of spatial correlation between warming and sediment deposition ultimately determined the nonlinear response of biofilm biomass (but not photosynthetic activity) to fluctuating levels of the two climate variables. Experimental results corroborated these predictions, probing the buffering effect of negative spatial correlation against extreme levels of warming and sediment deposition. Together, these results indicate that consideration of nonlinear response functions and local-scale patterns of correlation between climate drivers can improve our understanding and ability to predict ecological responses to multiple processes in heterogeneous environments.

The shadow model: how and why small choices in spatially explicit species distribution models affect predictions

The use of species distribution models (SDMs) has rapidly increased over the last decade, driven largely by increasing observational evidence of distributional shifts of terrestrial and aquatic populations. These models permit, for example, the quantification of range shifts, the estimation of species co-occurrence, and the association of habitat to species distribution and abundance. The increasing complexity of contemporary SDMs presents new challenges-as the choices among modeling options increase, it is essential to understand how these choices affect model outcomes. Using a combination of original analysis and literature review, we synthesize the effects of three common model choices in semi-parametric predictive process species distribution modeling: model structure, spatial extent of the data, and spatial scale of predictions. To illustrate the effects of these choices, we develop a case study centered around sablefish (Anoplopoma fimbria) distribution on the west coast of the USA. The three modeling choices represent decisions necessary in virtually all ecological applications of these methods, and are important because the consequences of these choices impact derived quantities of interest (e.g., estimates of population size and their management implications). Truncating the spatial extent of data near the observed range edge, or using a model that is misspecified in terms of covariates and spatial and spatiotemporal fields, led to bias in population biomass trends and mean distribution compared to estimates from models using the full dataset and appropriate model structure. In some cases, these suboptimal modeling decisions may be unavoidable, but understanding the tradeoffs of these choices and impacts on predictions is critical. We illustrate how seemingly small model choices, often made out of necessity or simplicity, can affect scientific advice informing management decisions-potentially leading to erroneous conclusions about changes in abundance or distribution and the precision of such estimates. For example, we show how incorrect decisions could cause overestimation of abundance, which could result in management advice resulting in overfishing. Based on these findings and literature gaps, we outline important frontiers in SDM development.

What is winter? Modeling spatial variation in bat host traits and hibernation and their implications for overwintering energetics

White-nose syndrome (WNS) has decimated hibernating bat populations across eastern and central North America for over a decade. Disease severity is driven by the interaction between bat characteristics, the cold-loving fungal agent, and the hibernation environment. While we further improve hibernation energetics models, we have yet to examine how spatial heterogeneity in host traits is linked to survival in this disease system. Here, we develop predictive spatial models of body mass for the little brown myotis (Myotis lucifugus) and reassess previous definitions of the duration of hibernation of this species. Using data from published literature, public databases, local experts, and our own fieldwork, we fit a series of generalized linear models with hypothesized abiotic drivers to create distribution-wide predictions of prehibernation body fat and hibernation duration. Our results provide improved estimations of hibernation duration and identify a scaling relationship between body mass and body fat; this relationship allows for the first continuous estimates of prehibernation body mass and fat across the species' distribution. We used these results to inform a hibernation energetic model to create spatially varying fat use estimates for M. lucifugus. These results predict WNS mortality of M. lucifugus populations in western North America may be comparable to the substantial die-off observed in eastern and central populations.

Effects of stochasticity on the length and behaviour of ecological transients

There is a growing recognition that ecological systems can spend extended periods of time far away from an asymptotic state, and that ecological understanding will therefore require a deeper appreciation for how long ecological transients arise. Recent work has defined classes of deterministic mechanisms that can lead to long transients. Given the ubiquity of stochasticity in ecological systems, a similar systematic treatment of transients that includes the influence of stochasticity is important. Stochasticity can of course promote the appearance of transient dynamics by preventing systems from settling permanently near their asymptotic state, but stochasticity also interacts with deterministic features to create qualitatively new dynamics. As such, stochasticity may shorten, extend or fundamentally change a system's transient dynamics. Here, we describe a general framework that is developing for understanding the range of possible outcomes when random processes impact the dynamics of ecological systems over realistic time scales. We emphasize that we can understand the ways in which stochasticity can either extend or reduce the lifetime of transients by studying the interactions between the stochastic and deterministic processes present, and we summarize both the current state of knowledge and avenues for future advances.

Predicting eastern equine encephalitis spread in North America: An ecological study

Eastern equine encephalitis (EEE) is a rare but lethal mosquito-borne zoonotic disease. Recent years have seen incursion into new areas of the USA, and in 2019 the highest number of human cases in decades. Due to the low detection rate of EEE, previous studies were unable to quantify large-scale and recent EEE ecological dynamics. We used Bayesian spatial generalized-linear mixed model to quantify the spatiotemporal dynamics of human EEE incidence in the northeastern USA. In addition, we assessed whether equine EEE incidence has predictive power for human cases, independently from other environmental variables. The predictors of the model were selected based on variable importance. Human incidence increased with temperature seasonality, but decreased with summer temperature, summer, fall, and winter precipitation. We also found EEE transmission in equines strongly associated with human infection (OR: 1.57; 95% CI: 1.52–1.60) and latitudes above 41.9°N after 2018. The study designed for sparse dataset described new and known relationships between human and animal EEE and environmental factors, including geographical directionality. Future models must include equine cases as a risk factor when predicting human EEE risks. Future work is still necessary to ascertain the establishment of EEE in northern latitudes and the robustness of the available data.

•

We collected EEE infections in humans and equines in the northeastern USA (2006–2019), at the county level.

•

We used reliably interpolated weather data from The PRISM Climate Group.

•

The first use of horse cases to predict human cases, controlling for weather and spatial effect.

•

Human risk was correlated with equine infection rates, year, latitude, temperature, and precipitation.

•

Cases increased in 2019 and above 41.9 degrees latitude, more studies are needed to confirm a northward shift.