Matthew C. Fitzpatrick

Climate change and the past, present, and future of biotic interactions.

Biotic interactions drive key ecological and evolutionary processes and mediate ecosystem responses to climate change. The direction, frequency, and intensity of biotic interactions can in turn be altered by climate change. Understanding the complex interplay between climate and biotic interactions is thus essential for fully anticipating how ecosystems will respond to the fast rates of current warming, which are unprecedented since the end of the last glacial period. We highlight episodes of climate change that have disrupted ecosystems and trophic interactions over time scales ranging from years to millennia by changing species’ relative abundances and geographic ranges, causing extinctions, and creating transient and novel communities dominated by generalist species and interactions. These patterns emerge repeatedly across disparate temporal and spatial scales, suggesting the possibility of similar underlying processes. Based on these findings, we identify knowledge gaps and fruitful areas for research that will further our understanding of the effects of climate change on ecosystems.

The projection of species distribution models and the problem of non-analog climate

By 2100, a quarter or more of the Earth’s land surface may experience climatic conditions that have no modern analog, with novel climates predicted to arise primarily in regions that currently support high levels of biodiversity (Williams et al. 2007). Further, global commerce will continue to transport species beyond long-standing dispersal barriers, potentially unleashing biological invaders into regions outside of those in which they evolved. Global climatic change and biological invasions will each have important and likely synergistic impacts on biodiversity. However, the emergence of non-analog climates (i.e., climatic conditions that do not presently exist) and the introduction of species to new biogeographical settings challenge our ability to anticipate these impacts because little information exists to predict how species may respond under novel environments.

Space can substitute for time in predicting climate-change effects on biodiversity

“Space-for-time” substitution is widely used in biodiversity modeling to infer past or future trajectories of ecological systems from contemporary spatial patterns. However, the foundational assumption—that drivers of spatial gradients of species composition also drive temporal changes in diversity—rarely is tested. Here, we empirically test the space-for-time assumption by constructing orthogonal datasets of compositional turnover of plant taxa and climatic dissimilarity through time and across space from Late Quaternary pollen records in eastern North America, then modeling climate-driven compositional turnover. Predictions relying on space-for-time substitution were ∼72% as accurate as “time-for-time” predictions. However, space-for-time substitution performed poorly during the Holocene when temporal variation in climate was small relative to spatial variation and required subsampling to match the extent of spatial and temporal climatic gradients. Despite this caution, our results generally support the judicious use of space-for-time substitution in modeling community responses to climate change.

Climate refugia: joint inference from fossil records, species distribution models and phylogeography

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial-interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia - fossil records, species distribution models and phylogeographic surveys - in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

Climatic drivers of hemispheric asymmetry in global patterns of ant species richness.

Although many taxa show a latitudinal gradient in richness, the relationship between latitude and species richness is often asymmetrical between the northern and southern hemispheres. Here we examine the latitudinal pattern of species richness across 1003 local ant assemblages. We find latitudinal asymmetry, with southern hemisphere sites being more diverse than northern hemisphere sites. Most of this asymmetry could be explained statistically by differences in contemporary climate. Local ant species richness was positively associated with temperature, but negatively (although weakly) associated with temperature range and precipitation. After contemporary climate was accounted for, a modest difference in diversity between hemispheres persisted, suggesting that factors other than contemporary climate contributed to the hemispherical asymmetry. The most parsimonious explanation for this remaining asymmetry is that greater climate change since the Eocene in the northern than in the southern hemisphere has led to more extinctions in the northern hemisphere with consequent effects on local ant species richness.

Observer bias and the detection of low-density populations

Monitoring programs increasingly are used to document the spread of invasive species in the hope of detecting and eradicating low-density infestations before they become established. However, interobserver variation in the detection and correct identification of low-density populations of invasive species remains largely unexplored. In this study, we compare the abilities of volunteer and experienced individuals to detect low-density populations of an actively spreading invasive species, and we explore how interobserver variation can bias estimates of the proportion of sites infested derived from occupancy models that allow for both false negative and false positive (misclassification) errors. We found that experienced individuals detected small infestations at sites where volunteers failed to find infestations. However, occupancy models erroneously suggested that experienced observers had a higher probability of falsely detecting the species as present than did volunteers. This unexpected finding is an artifact of the modeling framework and results from a failure of volunteers to detect low-density infestations rather than from false positive errors by experienced observers. Our findings reveal a potential issue with site occupancy models that can arise when volunteer and experienced observers are used together in surveys.

Ecological genomics meets community-level modelling of biodiversity: mapping the genomic landscape of current and future environmental adaptation

Local adaptation is a central feature of most species occupying spatially heterogeneous environments, and may factor critically in responses to environmental change. However, most efforts to model the response of species to climate change ignore intraspecific variation due to local adaptation. Here, we present a new perspective on spatial modelling of organism-environment relationships that combines genomic data and community-level modelling to develop scenarios regarding the geographic distribution of genomic variation in response to environmental change. Rather than modelling species within communities, we use these techniques to model large numbers of loci across genomes. Using balsam poplar (Populus balsamifera) as a case study, we demonstrate how our framework can accommodate nonlinear responses of loci to environmental gradients. We identify a threshold response to temperature in the circadian clock gene GIGANTEA-5 (GI5), suggesting that this gene has experienced strong local adaptation to temperature. We also demonstrate how these methods can map ecological adaptation from genomic data, including the identification of predicted differences in the genetic composition of populations under current and future climates. Community-level modelling of genomic variation represents an important advance in landscape genomics and spatial modelling of biodiversity that moves beyond species-level assessments of climate change vulnerability.

MaxEnt versus MaxLike: empirical comparisons with ant species distributions

MaxEnt is one of the most widely used tools in ecology, biogeography, and evolution for modeling and mapping species distributions using presence-only occurrence records and associated environmental covariates. Despite its popularity, the exponential model implemented by MaxEnt does not directly estimate occurrence probability, the natural quantity of interest when modeling species distributions. Instead, MaxEnt generates an index of relative habitat suitability. MaxLike, a newly introduced maximum-likelihood technique, has been shown to overcome the problem of directly estimating the probability of occurrence using presence-only data. However, the performance and relative merits of MaxEnt and MaxLike remain largely untested, especially when modeling species with relatively few occurrence data that encompass only a portion of the geographic range of the species. Using geo-referenced occurrence records for six species of ants in New England, we provide comparisons of MaxEnt and MaxLike. We show that by most quantitative metrics, the performance of MaxLike exceeds that of MaxEnt, regardless of whether MaxEnt models account for sampling bias and include greater model complexity than implemented in MaxLike. More importantly, for most species, the relative suitability index estimated by MaxEnt often was poorly correlated with the probability of occurrence estimated by MaxLike, suggesting that the two methods are estimating different quantities. For species distribution modeling, MaxLike, and similar models that are based on an explicit sampling process and that directly estimate probability of occurrence, should be considered as important alternatives to the widely-used MaxEnt framework.

Modeling Range Dynamics In Heterogeneous Landscapes: Invasion Of The Hemlock Woolly Adelgid In Eastern North America

Range expansion by native and exotic species will continue to be a major component of global change. Anticipating the potential effects of changes in species distributions requires models capable of forecasting population spread across realistic, heterogeneous landscapes and subject to spatiotemporal variability in habitat suitability. Several decades of theory and model development, as well as increased computing power and availability of fine-resolution GIS data, now make such models possible. Still unanswered, however, is the question of how well this new generation of dynamic models will anticipate range expansion. Here we develop a spatially explicit stochastic model that combines dynamic dispersal and population processes with fine-resolution maps characterizing spatiotemporal heterogeneity in climate and habitat to model range expansion of the hemlock woolly adelgid (HWA; Adelges tsugae). We parameterize this model using multiyear data sets describing population and dispersal dynamics of HWA and apply it to eastern North America over a 57-year period (1951-2008). To evaluate the model, the observed pattern of spread of HWA during this same period was compared to model predictions. Our model predicts considerable heterogeneity in the risk of HWA invasion across space and through time, and it suggests that spatiotemporal variation in winter temperature, rather than hemlock abundance, exerts a primary control on the spread of HWA. Although the simulations generally matched the observed current extent of the invasion of HWA and patterns of anisotropic spread, it did not correctly predict when HWA was observed to arrive in different geographic regions. We attribute differences between the modeled and observed dynamics to an inability to capture the timing and direction of long-distance dispersal events that substantially affected the ensuing pattern of spread.

Environmental and historical imprints on beta diversity: insights from variation in rates of species turnover along gradients

A common approach for analysing geographical variation in biodiversity involves using linear models to determine the rate at which species similarity declines with geographical or environmental distance and comparing this rate among regions, taxa or communities. Implicit in this approach are weakly justified assumptions that the rate of species turnover remains constant along gradients and that this rate can therefore serve as a means to compare ecological systems. We use generalized dissimilarity modelling, a novel method that accommodates variation in rates of species turnover along gradients and between different gradients, to compare environmental and spatial controls on the floras of two regions with contrasting evolutionary and climatic histories: southwest Australia and northern Europe. We find stronger signals of climate history in the northern European flora and demonstrate that variation in rates of species turnover is persistent across regions, taxa and different gradients. Such variation may represent an important but often overlooked component of biodiversity that complicates comparisons of distance–decay relationships and underscores the importance of using methods that accommodate the curvilinear relationships expected when modelling beta diversity. Determining how rates of species turnover vary along and between gradients is relevant to understanding the sensitivity of ecological systems to environmental change.