The consequences of climate change at an avian influenza ‘hotspot’

V. L. Brown, Pejman Rohani

Abstract

Avian influenza viruses (AIVs) pose significant danger to human health. A key step in managing this threat is understanding the maintenance of AIVs in wild birds, their natural reservoir. Ruddy turnstones (Arenaria interpres) are an atypical bird species in this regard, annually experiencing high AIV prevalence in only one location—Delaware Bay, USA, during their spring migration. While there, they congregate on beaches, attracted by the super-abundance of horseshoe crab eggs. A relationship between ruddy turnstone and horseshoe crab (Limulus polyphemus) population sizes has been established, with a declining horseshoe crab population linked to a corresponding drop in ruddy turnstone population sizes. The effect of this interaction on AIV prevalence in ruddy turnstones has also been addressed. Here, we employ a transmission model to investigate how the interaction between these two species is likely to be altered by climate change. We explore the consequences of this modified interaction on both ruddy turnstone population size and AIV prevalence and show that, if climate change leads to a large enough mismatch in species phenology, AIV prevalence in ruddy turnstones will increase even as their population size decreases.

1. Introduction

Avian influenza viruses (AIVs) infect many bird species throughout the world, usually with relatively modest prevalence [1]. In some sites, however, the confluence of biotic and abiotic factors results in the proliferation of AIVs among multiple host species, leading to unusually high prevalence. Surveillance has identified Delaware Bay as such a ‘hotspot’, with virus isolation levels in shorebirds (ruddy turnstones in particular) exceeding cumulative figures across global surveillance sites by more than an order of magnitude [2]. Recent dissection of AIV epidemiology at Delaware Bay has identified the factors that generate this phenomenon [3], with food source (horseshoe crab eggs) availability during shorebird migration considered key. The interaction between these consumer (migratory shorebirds) and resource (horseshoe crabs) species is central to our understanding of both ruddy turnstone population biology and the contribution of these birds to the remarkably high AIV prevalence in Delaware Bay, a site of ‘hemispheric importance’ [2]. It is not surprising, therefore, that understanding and predicting the impact of climate change on horseshoe crab spawning times and bird migration timing [4], and their attendant consequences for AIV maintenance and spread, have been identified as both timely and important.

Horseshoe crabs have, for millions of years [5], existed in large numbers along the mid-Atlantic coast of North America [6]. Each year, adult crabs (more than 9 years old) travel to the beaches of Delaware Bay to spawn and lay their eggs [5], more than 99 per cent of which are lost through predation or exposure to the elements [7]. Declining populations of horseshoe crabs in recent years—owing to a combination of over-harvesting and habitat loss—are hypothesized [8] to account for declining numbers of ruddy turnstones [6] and other shorebirds [8].

Several studies have investigated the possible impact of climate change on bird populations [911]. Empirical evidence suggests that long-distance migrants with a seasonal food peak are more likely to experience population declines than ‘residents and short-distance migrants’ [11]. Ruddy turnstones are hypothesized to be impacted by climate-related changes in two distinct ways. Firstly, seasonal migration patterns may shift, with studies suggesting that migrating birds are now leaving their wintering grounds earlier, possibly in response to a changing climate [10]. Secondly, food resource (horseshoe crab eggs) availability may change, with evidence demonstrating that timing of horseshoe crab spawning is influenced by water temperature [12]. Moreover, because these changes in migration and resources are in response to environmental cues from different geographical regions, the predicted heterogeneous impact of climate change across the globe [13] means that their dynamical impacts are likely to be complex, rather than simply additive.

In this paper, we use a mathematical model to explore the consequences of an altered interaction between ruddy turnstones and horseshoe crabs in Delaware Bay as a result of climate change. Employing the most detailed parametrized model available of AIV dynamics in Delaware Bay—a model which includes three host species, seasonal forcing and multiple transmission routes [3]—we explore the consequences of two prominent potential scenarios: (I) climate change impacts migration in ruddy turnstones, altering their arrival date in Delaware Bay and (II) climate change affects horseshoe crab spawning, altering the timing of peak numbers of horseshoe crab's eggs on the Delaware Bay beaches. We show that a change in the availability of horseshoe crabs impacts the prevalence of AIV in ruddy turnstones and further demonstrate the impact of changing ruddy turnstone migration and horseshoe crab spawning on both AIV prevalence and population numbers of ruddy turnstones. We conclude that a small shift in either of these events will not qualitatively affect the dynamics of the system, but that the ecological and epidemiological consequences will be significant with greater mismatch in phenology.

2. Material and methods

We investigated the interaction between ruddy turnstones and horseshoe crabs using a model designed to capture the key interactions occurring in Delaware Bay. The model includes two duck species (one resident in Delaware Bay and the other a short-distance migrant that winters in Delaware Bay) as well as ruddy turnstones, who winter in South America, pause in Delaware Bay on their northbound migration for just under a month, then breed in the Arctic. The model includes seasonal breeding and transmission for all species and has a density-dependent birth rate in ruddy turnstones to reflect their dependence on horseshoe crabs (see below). An environmental component for the virus is also included, with transmission assumed to occur both directly and via the environment [14]. The model is fully described in the electronic supplementary material.

Our focus here is the birth rate from the model for ruddy turnstones: b(t)(1 − θN(t)/E)N(t), where b(t) is the (seasonal) birth rate; θ is a shape parameter; N(t) is the (time-dependent) total ruddy turnstone population size and E is the number of horseshoe crabs. Details for calculating E are given in the electronic supplementary material. To incorporate the seasonality of horseshoe crab spawning (and hence egg availability), we calculated E from an exponential cosine function (Embedded Image, where Embedded Image is the midpoint of the ruddy turnstones' stay in Delaware Bay), with maximum E occurring as the ruddy turnstones are usually in Delaware Bay, and declining steeply away from that point. To address scenario I, we model a temporal shift in the ruddy turnstones' arrival in Delaware Bay. A single shift corresponds to a difference of three weeks and we consider two arrival times of either side of the current one. For scenario II, we shift the timing of peak horseshoe crab spawning by a month at a time and consider the two months prior to current spawning timing.

3. Results

We tested the two potential scenarios outlined in §1. Figure 1 presents the results from scenario I (i.e. changing the timing of ruddy turnstone migration) using our deterministic model. We calculated the number of horseshoe crabs available to the birds when they are in Delaware Bay from the exponential cosine curve given in §2 (for sensitivity to this function see the electronic supplementary material). The main panel of figure 1 shows the function for E, and the offset panels (figure 1ae) show the prevalence curves in all three hosts resulting from a change in migration timing (changed by three weeks each time).

Figure 1.

Exploring the consequences of potential scenario I. The assumed function for the number of horseshoe crabs, E, against time (main plot). The offset figures (ae) show prevalence throughout the year for each host species (for the final year of a 500 year simulation), with arrows showing when the ruddy turnstones were present in Delaware Bay.

As ruddy turnstones deviate from their current arrival time in Delaware Bay, AIV prevalence increases (figure 1d,e). Arrival earlier in the year (figure 1e) coincides with a higher AIV prevalence in the resident ducks (and hence higher virion concentration in the environmental reservoir). This, combined with two reasonable assumptions in the model—density-dependent direct transmission and seasonal breeding—leads to a greater number of ruddy turnstones present in Delaware Bay if they migrate a few weeks earlier. The reduction in horseshoe crab eggs, E, is not enough to significantly decrease the ruddy turnstone population size, so transmission is higher compared with later arrival. The latter increase in prevalence (figure 1d) occurs as their timing overlaps with the main peak in prevalence in resident ducks, sparking larger peak prevalences in both species than is otherwise seen. Current evidence from migrating bird species suggests their migration dates are advancing [10]. Applied to our results, this suggests AIV prevalence in ruddy turnstones may increase, all else being equal.

The results from scenario II (i.e. changing the timing of horseshoe crab spawning) are presented in figure 2. The offset figures show the AIV prevalence in ruddy turnstones and their population size in three different scenarios—the current situation and two earlier peaks in horseshoe crab spawning. As in figure 1, we find that for a small alteration in timing, prevalence (and population size) is not greatly affected, but for the earliest peak in horseshoe crab spawning both ruddy turnstone population size and AIV prevalence have declined hugely (figure 2c; for equivalent results for scenario I, see the electronic supplementary material).

Figure 2.

Three functions for peak horseshoe crab spawning, shifted through time (main plot). The offset figures show AIV prevalence and population size in ruddy turnstones for each of these functions (for the final year of a 500 year simulation). The blue bar shows when the tuddy turnstones are in Delaware Bay.

4. Discussion

We considered the consequences of climate change on both AIV prevalence and population sizes in ruddy turnstones. We compared two contrasting potential scenarios—that climate change modifies either ruddy turnstone migration timing or the horseshoe crab spawning season. Qualitatively, our results indicate that, as the mismatch between supply and demand for horseshoe crab eggs increases, so will the ruddy turnstone population decline. An unexpected finding, however, was that a change in the arrival time of ruddy turnstones in Delaware Bay is predicted to generate greater variability in outcomes. From stochastic (Markov process) realizations of the model, we observe a notable increase in AIV extinction frequency (table 1), which clearly dampens the potential for cross-species transmission. However, in instances where AIV persists, prevalence and, as a result, spillover potential to domestic birds, is substantially enhanced (electronic supplementary material, figure S3).

View this table:
Table 1.

Results of stochastic simulations (150 simulations for each time period). Parameters for stochastic model as in electronic supplementary material, table S1. Extinction probability is calculated as the number of simulations for which AIV does not persist in any of the three host species in the model after 500 years of simulation.

This research has exposed the potential reduction in population size faced by North American ruddy turnstones. In addition to current anthropogenic factors, such as harvesting of horseshoe crabs [8], which directly affect resource supply, ruddy turnstones are vulnerable to the effects of climate-driven shifts in phenological synchrony of horseshoe crab spawning and ruddy turnstone spring migration. These events, in turn, bear on AIV transmission in Delaware Bay, though the overall impact is predicted to be nuanced. Our analyses indicate a relatively small change in AIV prevalence among migratory ducks, but substantial effect on the phase and amplitude of AIV peak prevalence in resident ducks. Crucially, this affects Delaware Bay's role as a potential AIV gateway site: a change in ruddy turnstone migration biology may well to lead to an increase in peak prevalences among two of the host species in the system, and hence a possible increase in emergence threat.

Acknowledgements

This work was supported by the James S. McDonnell Foundation and the National Science Foundation (DEB-0917853). P.R. was also supported by the RAPIDD program of the Science and Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.

  • Received July 6, 2012.
  • Accepted August 2, 2012.

References

View Abstract