Abstract

Climate is the dominant factor determining the distribution and abundance of most insect species. In recent years, the issue of climatic changes caused by human activities and the effects on agriculture has raised concern. General circulation model scenarios were applied to a bioclimatic model of Melanoplus sanguinipes to assess the potential impact of global warming on its distribution and relative abundance. Native to North America and widely distributed, M. sanguinipes is one of the grasshopper species of the continent most responsible for economic damage to grain, oilseed, pulse, and forage crops. Compared to predicted range and distribution under current climate conditions, model results indicated that M. sanguinipes would have increased range and relative abundance under the three general circulation model scenarios in more northern regions of North America. Conversely, model output predicted that the range of this crop pest could contract in regions where climate conditions became limiting.

1. Introduction

Climate is the dominant factor determining the distribution and abundance of most insect species [1]. The issue of climatic changes caused by human activities and the effects on agriculture has raised concern in recent years. The overall global temperature has increased C over the last 100 years, with the 1990’s being the warmest decade on record [2]. Climate change scenarios using low greenhouse gas emissions suggest that temperatures will increase by 1– C over the next 100 years and temperatures have been predicted to increase by 3.5– C for scenarios with high gas emission [3]. However, Walther et al. [4] suggest that species respond to regional changes that are highly heterogeneous and not to approximated global averages. Many species have already responded to regional conditions that have occurred during the 20th century. In a study of 694 animal and plant species, Root et al. [5] investigated the change in timing of events over the past 50 years and reported that changes in timing of spring events (breeding, blooming) occurred 5.1 days earlier per decade. Warming conditions may impact grasshopper populations by extending the growing season, altering the timing of emergence from overwintering sites, increasing growth and development rates, shorting generation times, increasing the numbers of eggs laid, and changing their geographic distribution [6, 7].

Analogue scenarios which make use of existing climate data are useful to identify geographic regions that may be susceptible to establishment of insects, when comparing the results of climate change scenarios to those regions where the species in question is already established [8]. However, the magnitude of predicted temperature change associated with climate change is not within the historical experience of modern agriculture. Hence, it is unlikely that we can use historical data as analogues to predict the impact of climate change on pest species. As a result, simulation models have been used to assess impact and related system vulnerability due to climate change.

Bioclimate simulation models have been used successfully to predict the distribution and extent of insect establishment in new environments [912]. Bioclimatic modeling software, such as CLIMEX, enables the development of models that describe the potential distribution and relative abundance of a species based on climate [1, 13]. CLIMEX derives an Ecoclimatic Index (EI) which describes the suitability of specific locations for species survival and reproduction. Model parameters include temperature (TI), diapause (DI), light (LI), moisture (MI), heat stress (HS), cold stress (CS), wet stress (WS), and dry stress (DS). The EI values are obtained by combining a Growth Index (GI) with stress indices (dry, wet, cold, and hot) that describe conditions that are unfavourable for growth.

Native to North America and widely distributed, Melanoplus sanguinipes (Fabricius) (Orthoptera: Acrididae) is responsible for more economic damage to grain, oilseed, pulse, and forage crops than any other grasshopper species [1416]. A bioclimate model was developed to predict the potential distribution and relative abundance of M. sanguinipes, within Canada [17]. Ecological sensitivity analyses were then conducted using incremental scenarios for all combinations of temperature (0, +1, +2, +3, +4, +5, +6, and + C of climate normal temperature for each grid) and of precipitation ( 60%, 40%, 20%, 10%, 0%, 10%, 20%, 40%, 60% of climate normal precipitation for each grid). Compared to predicted range and distribution under current climate conditions, model results indicated that M. sanguinipes would have increased range and relative abundance for temperature increases between C and C. The model predicted that the range of this crop pest could be extended to regions that are not currently used for agricultural production in North America. Mika et al. [18] stated that at an ecosystem level, climatic variables will vary both spatially and temporally. Therefore, they suggested that the widely accepted and more commonly used general circulation models (GCMs) should be used in conjunction with bioclimate models, rather than incremental scenarios. Further, they encouraged the application of multiple GCMs due to the variability of climate projections between models.

The objective of this study was to use the bioclimate model for M. sanguinipes [17] to assess the impact of three general circulation models on population distribution and relative abundance across North America.

2. Methods

The bioclimatic model for M. sanguinipes, developed using CLIMEX 2.0 [19], has been previously described [17]. CLIMEX is a dynamic model that integrates the weekly responses of a population to climate using a series of annual indices. It uses an annual Growth Index to describe the potential for population growth as a function of soil moisture and temperature during favourable conditions, and Stress Indices (cold, wet, hot, and dry) to determine the effect of abiotic stress on survival in unfavourable conditions. The weekly Growth Index is a function of temperature (TI), diapause (DI), and moisture (MI). The growth and stress indices are calculated weekly and then combined into an overall annual index of climatic suitability, the Ecoclimatic Index (EI), that ranges from 0 for locations at which the species is not able to persist to 100 for locations that are optimal for the species [17]. Model parameter values are listed in Table 1. Initial parameter values were obtained from published papers. Model parameters were then adjusted to ensure that EI ≥ 30 in geographical regions historically affected by M. sanguinipes, indicating that climatic conditions were favorable for development of densities associated with crop loss. Historical grasshopper population data were used for model validation. Annual surveys of abundance of adult grasshoppers have been conducted in Saskatchewan since 1931 [20]. Relative abundance was validated by comparison with adult grasshopper survey data from Saskatchewan over the period of 1970 to 2004 [17]. The model was tested by comparing the occurrence of observed life history events against those predicted by the model.

Climate change projections were obtained from the Intergovernmental Panel on Climate Change [21] as monthly means for three GCMs, based on current climate, 30 yr average (1961–1990) dataset (A1B emission scenario) (CRU—Climate Research Unit, East Anglia, UK). The three GCMs selected were CSIRO Mark 3.0 (CSIRO, Australia), NCAR273 CCSM (National Centre for Atmospheric Research, USA), and MIROC-H (Centre for Climate Research, Japan). All three had relatively small horizontal grid spacing and the requisite climatic variables at a temporal resolution appropriate for CLIMEX. The data were pattern-scaled to develop individual change scenarios relative to the base climatology [22]. The three models cover a range of climate sensitivity, defined as the amount of global warming for a doubling of the atmospheric CO2 concentration compared with 1990 levels [23]. The respective sensitivities are: CSIRO Mark 3.0 ( C), NCAR-CCSM ( C), and MIROC-H ( C).

The resulting database was queried to analyze data at a regional scale. A geographic rectangle, latitude by longitude, was used to delineate a regional template. The defined region was approximately the size and shape of Colorado (270,000 km2) and, for each of the datasets, consists of 112 grid cells. Specific regions, based on latitude and longitude coordinates, were defined and output (averaged across the region) was generated for detailed analysis. The datasets permitted comparison of variables, both spatially and temporally (weekly intervals). Analyses were based on values centered on six locations including Lubbock, Texas ( N, W), Gillette, Wyoming ( N, W), Lincoln, Nebraska ( N, W), Saskatoon, Saskatchewan ( N, W), Peace River, Alberta ( N, W), and Fairbanks, Alaska ( N, W).

Contour maps were generated by importing EI values into geographic information system software, ArcView 8.1 [24]. Final EI values were displayed in the four categories defined above: “Unfavourable,” “Suitable;” “Favourable;” and “Very Favourable.”

3. Results and Discussion

Comparisons were made to determine if differences in baseline climate data would result in differences in output. The New et al. [25] climate data represents a splined world grid dataset. The EI output the baseline CRU data agreed with that produced using the New et al. [25] climate data set in Olfert and Weiss [17]. Initially, there appeared to be some differences in model output between the two approaches for Peace River and Saskatoon (Table 2). Olfert and Weiss [17] reported that the EI values for Peace River and Saskatoon were 24 and 30, respectively. This study showed that EI values for Peace River and Saskatoon were 14.7 and 26.7 (Table 2). These differences occurred because the original paper reported values for single grid cells. However, the current analysis was based on averages across large regions that are composed of 112 grid cells. When single grid cells for Peace River and Saskatoon were examined in the current study, it was found that EI values were indeed 24 and 30.

Results, based on the CRU data for current climate, indicated that M. sanguinipes would have highest EI values across most of the Great Plains of North America, extending from northern Texas to southern Saskatchewan (Figure 1). These results agreed with the distribution of M. sanguinipes as described by Riegert [20] and Pfadt [26]. Compared to these results, each of the three GCMs resulted in large differences for most model parameters, particularly EI (Figures 24; Table 2). Across North America, the overall mean EI values were 4.9 (CRU), 7.5 (CSIRO MARK 3.0), 7.9 (MIROC-H), and 7.3 (NCAR273 CCSM). Olfert and Weiss [17] grouped ecoclimatic indices into four categories: “Unfavourable” (EI = 0–5), “Suitable” (EI = 5–20), “Favourable” (EI = 20–30), and “Very Favourable” (EI ≥ 30). Unfavourable described regions where M. sanguinipes would be very rare or may not occur; “Suitable” defined areas were grasshoppers would occur, usually in low densities; “Favourable” defined areas were densities could be high enough to result in crop loss; “Very Favourable” defined areas where grasshoppers regularly occur in high enough densities that result in crop loss. Based on this study, the extent of the area predicted to be “Very Favourable” were 11.2% (CRU), 16.2% (CSIRO MARK 3.0), 16.2% (MIROC-H), and 18.1% (NCAR273 CCSM) of North America.

Species are more vulnerable to variations in temperature and precipitation when located near the outer limits of their geographic range than when located in the core area of the range. Sutherst et al. [19] defined a core area as a region with high EI values and little or no stress. Populations near the outer limits of the core area spend a greater amount of time in climates that are marginally suitable (exposed to climatic stress), while populations near the core experience a greater amount of time in favourable conditions (minimal exposure to climatic stress). In this study, EI values tended to increase in a northwestern direction and decrease for southern locations when the three GCMs were applied to the bioclimate model for M. sanguinipes. The percent of area (on a regional basis) with EI ≥ 20 varied across North America. For example, under current climate conditions (CRU), the model predicted that 0% of the Fairbanks region had EI ≥ 20 (Table 3). This value increased to as much as 57% of the area under conditions predicted by NCAR273 CCSM. As a result, the increase in the biological suitability of Fairbanks, AK, due to climate change was predicted to be similar to that of Lincoln, NE, under current climate conditions (CRU). In turn, the model predicted that the area surrounding Lincoln, NE, where EI ≥ 20 would decrease to 6.3% (NCAR273 CCSM), compared to 59.8% under current climate conditions (CRU).

As indicated, there were regional differences across North America in output of the bioclimate model for M. sanguinipes when the three different GCMs were applied (Figures 24). The application of CSIRO MARK 3.0 climate data resulted in a northward shift of areas predicted to have reduced suitability for grasshopper populations within the southern Great Plains, relative to current climate conditions (CRU). There was a significant reduction in EI values in states such as Colorado, Wyoming, and Missouri (Figure 2). In northwest Texas, the EI values were predicted to decrease to less than 10. In more northern regions, however, EI values were predicted to be higher in Alaska, northern Alberta, and Saskatchewan, relative to current climate conditions (CRU). Output based on the MIROC-H dataset resulted in a northwest shift of regions with EI ≥ 20 (Figure 3). Compared to current climate data (CRU), the MIROC-H GCM predicted that the overall area suitable for M. sanguinipes in the USA would be less than under current climate conditions. However, the suitable areas along the Rocky Mountains were observed to increase somewhat. The MIROC-H dataset predicted large EI increases across most of the Canadian prairies, and extending northwest to include a continuous area northwest to Peace River, Alberta. Of the three GCMs, MIROC-H output resulted in the largest, continuous areas with EI ≥ 20. Unlike CSIRO MARK 3.0 and MIROC-H, NCAR273 CCSM predicted a reduction in EI values for eastern North America. This GCM also predicted increased EI values in the interior of British Columbia.

In order to assess the potential impact of climate change in a more regional context, the resulting database was queried to analyze data at six regional locations between Lubbock, Texas, and Fairbanks, Alaska. Overall, the largest differences in EI values were observed at northern and southern regions of North America. The shifts in EI values were less in central locations. Compared to current climate (CRU), EI values derived from GCMs resulted in increased EI for the areas surrounding the three northern regions Saskatoon, Peace River, and Fairbanks (Figures 14, Table 2). The magnitude of the increase in EI values, based on regional means, was 252%, 77%, and 33% greater for Fairbanks, Peace River, and Saskatoon, respectively, than those under current climate conditions. As a result, warming conditions were predicted to result in increased potential for M. sanguinipes outbreaks in these three regions. Outbreaks of M. sanguinipes have been recently reported in northern areas of North America. This species has been reported to be a sporadic, potentially damaging grasshopper pest of small grain crops in Alaska [27] and recent outbreaks of grasshoppers have been reported in the Peace River region of Alberta [28]. The three southern locations (Lincoln, Gillette, and Lubbock) had lower EI values when GCMs were used as inputs into the model. Relative to the EI values under current climate, the regional mean EI values for Gillette, Lincoln, and Lubbock were predicted to be 23%, 29%, and 69% less, respectively. The regional responses to model input varied for the three GCMs (Table 3). The MIROC-H GCM resulted in the largest increase in EI for the Peace River and Saskatoon regions, while the NCAR273 CCSM model resulted in the largest increase for the Fairbanks region. Of the three more southern regions, Lubbock exhibited the largest decrease in EI values for MIROC-H.

The weekly temperature index (TI) describes the weekly response of M. sanguinipes to the daily temperature cycles that occur during the growing season. Melanoplus sanguinipes overwinters in the egg stage. The timing and duration of spring hatch is influenced by the level of embryonic development going into winter, natural enemies, and soil temperature and moisture [29, 30]. In northerly regions, M. sanguinipes produces only one generation per year; in more southerly areas, a small proportion of the eggs oviposited do not enter diapause and may result in a lesser second generation. This species prefers warm, dry weather conditions. Warm temperatures early in spring favour nymphal development and in turn the timing of adulthood. Conversely, cool and wet conditions in spring results in increased nymphal mortality and delayed development. Crop loss due to feeding damage can occur throughout the growing season. Newly emerged seedlings in spring are most vulnerable, however, gradual plant defoliation may also contribute to decreased crop yield and quality [14, 16]. Later in the growing season, an extended, warm fall influences the longevity of adults, allowing them to continue reproducing until freeze-up [29, 30]. As a result, economic infestations are often associated with a prolonged period of consecutive seasons with above-normal temperatures [29]. Intermittent warm seasons tend to result in fluctuating populations [31, 32]. The GCM datasets, associated with temperatures that are warmer than CRU values, resulted in increased TI values for northern regions and reduced TI for southern regions. Olfert and Weiss [17] reported that incremental scenarios of C and C resulted in increased TI values and increases in both EI and the potential area of Canada that would potentially be exposed to grasshopper outbreaks. At Fairbanks, TI values increased from 12.0 (CRU) to 20.1 (NCAR273 CCSM), resulting in more favourable temperatures during the growing season. Changes in central North America were less dramatic. Temperature indices in the Saskatoon region were predicted to increase from 30.7 (CRU) to 37.8 (NCAR273 CCSM). Excessively warm temperatures have been shown to hinder grasshopper populations [29, 33]. Output indicates that increased temperatures would result in higher heat stress (HS) values in northern Texas and Nebraska.

The growth index (GI) is a weekly thermo hydrological index that describes conditions that are favourable for growth. CLIMEX outputs the number of weeks where the growth index is nonzero, effectively determining the length of the growing season. Growing season length and cold stress accumulation are two factors that limit the potential for population growth in the Fairbanks region. Increased temperatures were predicted to not only decrease the rate of cold stress accumulation, but to also increase both the diapause index (DI) and the length of the growing season from 14 weeks to 17-18 weeks. The growing season in the Peace River region was predicted to increase from 19 weeks to 22 weeks and would result in a growing season that is similar to the current growing season in the Saskatoon region. Mills [34] predicted that regions north of N and west of W have soils that are suitable for agricultural production and that climate change could positively impact small grain production in the area. This would suggest that M. sanguinipes populations could become established in these new agricultural areas in the event that they become accessible in the future. In southern regions, however, mean GI values and the number of weeks where GI values were positive decreased. Output indicated that prolonged periods of warm temperatures during the growing season could limit potential for grasshopper population growth. Extreme heat and drought tends to reduce crop growth while increase grasshopper feeding activity. Mukerji et al. [32] reported that increased competition for food can also result in population decline due to high mortality because of starvation.

In conclusion, bioclimatic models have proven useful for studies investigating the potential impact of climate on insect populations. However, some cautions have been expressed regarding the utilization of this approach including: (i) biotic interactions may not remain the same over time (adaptation can, and is likely to, occur); (ii) genetic and phenotypic composition of populations may change over time and space; (iii) most species have some limitation to dispersal [35, 36]. In the instance of M. sanguinipes, the impact of biotic factors such as natural enemies (e.g., diseases, parasites) must also be considered. For example, termination of several grasshopper outbreaks in Canada were attributed to cool, wet weather and epizootics of Entomophthora grylli Fres. [20, 37]. Even though conditions may be predicted to be conducive to grasshopper populations under climate change, diseases could result in population decline. In these instances, bioclimate and GCMs may not account for changes in population, and may overestimate populations. To address these naturally occurring phenomena, bioclimate modeling of grasshopper populations will benefit from a multitrophic approach (host plants—grasshoppers—natural enemies).

Acknowledgment

D. Giffen is gratefully acknowledged for preparation of the maps.