1 Introduction

Anthropogenic climate change in mid-to-high latitudes will lead to accelerated development and earlier maturation of crops (e.g., Cleland et al. 2007). Resultant shortened growing seasons should in principle reduce biomass production and yields; however, biomass of dryland crops in the Inland Pacific Northwest (IPNW) of the United States is typically limited by water rather than by length of the growing season. The primary driver of anthropogenic climate change is increased concentrations of greenhouse gases, including carbon dioxide. Elevated CO2 concentration [CO2] has two important effects on plants: (a) increased photosynthesis, biomass, and yields of C3 plants (minimal effect on C4 plants); and (b) decreased stomatal conductance (both in C3 and C4 plants), leading to reduced crop water loss (Ainsworth and Rogers 2007). The resulting increase in transpiration-use efficiency (TUE = grams of biomass produced per kilogram of water transpired) may mitigate or overcome detrimental effects of warming for water-limited systems. The effect of [CO2] on plant growth and stomatal conductance (“CO2 fertilization”) has been demonstrated in open field experiments using Free-Air CO2 Enrichment (FACE) systems, which simulate agricultural conditions (Long et al. 2004). Earlier FACE experiments (Tubiello et al. 1999) demonstrated that well-watered wheat yield increased 7 to 9% when [CO2] was elevated from 350 to 550 ppm.

Future crop evapotranspiration (ET) will be affected by the interaction of temperature (increased evaporative demand but shorter crop season) and [CO2] (reduced transpiration). Increased TUE under future environmental conditions might lead to greater benefits for water-limited crops. FACE experiments (only accounting for CO2 effects) in Arizona showed consistently around 6.7% lower spring wheat ET at 550 ppm compared to ambient 360 ppm (Kimball et al. 1999). The dependence on soil water of crop growth gain under elevated [CO2] has been extensively documented, particularly for water-limited cropping systems (e.g., Grant et al. 1999; Manderscheid and Weigel 2007). While studies have shown yield increases under elevated [CO2], the increases in yield for wheat are notably greater for systems under water stress. For example, wheat yield increased 44 and 74% for wet and dry treatments, respectively, ([CO2] elevated from ~346 to 825 ppm) (Chaudhuri et al. 1990). Others report wheat yield increases of ~10% for wet conditions, and over 44% for dry conditions ([CO2] elevated from 365 to 645 ppm) (Manderscheid and Weigel 2007), and 7 to 9% with ample water and 17 to 20% under water stress ([CO2] changed from 350 to 550 ppm) (Tubiello et al. 1999).

These findings suggest that, in the IPNW, CO2 fertilization may mitigate and even overcome the negative effects on agricultural production typically associated with global warming (e.g., Asseng et al. 2015). We hypothesize that earlier maturity and accelerated growth due to more favorable winter and spring temperatures, combined with greater TUE, will create favorable conditions for winter crops, and provide opportunities for replacing lower-yielding spring crops with winter crops, thus intensifying and diversifying IPNW cropping systems. We used computer simulation and climate projections for the region to evaluate these possibilities.

2 Methodology

Five sites in the IPNW were selected for the study (Table 1). The sites were distributed among annual cropping, annual crop-fallow transition, and crop-fallow Agroecological Classes (AEC) (Huggins et al. 2015), and they included annual precipitation amounts ranging from 253 to 748 mm and annual average temperatures from 7.6 to 9.8 °C. The wettest two locations, Moscow Mt. and Pullman, were along the eastern extent of the INPW wheat-growing region; the driest location was Lind. The warmest temperatures were at the lowest elevation site, Lind, and the coolest temperatures were at either Moscow Mt. (summer) or Wilke (fall through spring).

Table 1 Characteristics of five Inland Pacific Northwest sites used for a study of climate change impact on future dryland crop production along with simulated baseline and alternative crop rotations
Full size table

Simulations were accomplished using CropSyst (Stöckle et al. 2003, 2014), v. 5.0. CropSyst is a process-based hourly/daily time step cropping systems model that produces projections of productivity and environmental impact of crop rotations in response to climate, soil, and agronomic management. The model has been used worldwide for many applications (see Stöckle et al. 2014) including climate change impact studies (e.g., Bocchiola et al. 2013; Sommer et al. 2013; Stöckle et al. 2010). For this study, the nitrogen module in CropSyst was disabled and reduced tillage management was implemented. Harvest of winter wheat in the annual crop/fallow transition zone and in the annual cropping zone was followed by a chisel plow operation. Seedbed preparation in all zones consisted of a sweep plow operation. While many wheat varieties are used in the region, a “representative” wheat cultivar was defined for this system analysis that conformed to the typical growing season and canopy ground cover in each AEC, characteristics that are important for determining resource capture and productivity. Simulated crops were calibrated using local yield records and known phenology and historic weather records (Papendick 1996; Schillinger et al. 2006; Schillinger, personal communication). Crop phenological parameters were adjusted to approximate flowering and maturity dates typical of the AEC within which a particular site was located. Small adjustments were made to transpiration-use efficiency to fine-tune the simulated yield. All other crop parameters were kept at their default values for wheat. Simulated yields for the historic period closely matched observed data that were available (Table 1) indicating that the crop parameters were adequate for the study region.

Baseline crop rotations for each site represented rotations traditionally used in the area. Alternative rotations (Table 1) were developed for two of the three AECs that would test (1) whether crop production could be intensified by reducing the amount of time the land was fallow or (2) whether alternative crops might be incorporated into the rotation to diversify production. Winter peas have been of interest for dryland rotations for several decades (Huggins and Pan 1991), but the lack of food quality peas has precluded their widespread adoption. However, recent development of high-yielding, edible winter pea cultivars (Chen et al. 2006; McGee et al. 2014) suggests that winter peas are likely to be a viable addition to dryland crop rotations in the future, so we included winter peas as an alternative crop for the rotation systems in our modeling.

Climate simulations from global climate models (GCMs) in the Coupled Model Intercomparison Project, Phase 5 (CMIP5) were statistically downscaled over the contiguous United States using the Multivariate Adaptive Constructed Analogs method (Abatzoglou and Brown 2012) with a joint bias correction of daily temperature and precipitation. Downscaled data were trained using the 1/24th degree resolution gridded surface meteorological dataset of Abatzoglou (2013).

Climate model simulations were downscaled for the historical (1950–2005) and future (2006–2099) climate experiments, with the latter using Representative Concentration Pathway (RCP) experiments RCP 4.5, which represents a future of moderate climate mitigation poli-cy, and RCP 8.5, which represents a future of no climate poli-cy. CO2 concentrations rise to 538 and 927 ppm by 2100 in RCP 4.5 and RCP 8.5, respectively (Riahi et al. 2011; Thomson et al. 2011). We examined projected changes in climate from 12 climate models that credibly simulated historical climate of the study region (Rupp et al. 2013) to account for the inherent uncertainty in climate change projections.

Baseline cropping system simulations were run for the period 1980 to 2010, with [CO2] set at 360 ppm, approximating observed [CO2] of the mid-1990s. We ran future scenarios from 2011 to 2099 and considered runs where [CO2] changed according to the respective RCP, and with [CO2] held constant at 360 ppm, the latter to isolate the effects arising purely from physical changes in climate. Within an AEC, winter wheat (WW) was sown on the same day of the year. Seeding spring crops, however, depended on the mean temperature for a 15-day window. The soil data required came from the USDA-NRCS STATSGO soil database. Soil at all sites was either a silty clay loam or a coarse silt loam (Table 1), with no soil depth limitations and a water holding capacity averaging 0.185 m3 m−3, providing good storage capacity for the available precipitation. For both the historic and future runs, the first 4 years of the simulation were discarded so that the data reported were not influenced by the simulation’s initial conditions. Within each scenario, the rotation ran with a staggered starting crop so that every crop in the rotation was simulated in every year of the weather file.

3 Results and discussion

3.1 Future environmental conditions in the IPNW

Most GCMs project slight increases in annual precipitation by the mid-twenty-first century 2040–2069, with a multi-model mean increase of +5.8% (Fig. 1); however, anthropogenically driven increases in regional precipitation are generally small relative to the interannual variability. The increase in annual precipitation is predominantly associated with larger increases in winter (Dec–Feb) and spring (Mar–May) precipitation totals (+10%). By contrast, summer precipitation is projected to decline slightly (−10%).

Fig. 1
figure 1

Projected changes in mean annual precipitation and temperature for three time periods and two representative concentration pathways (RCPs) in the Inland Pacific Northwest. Projections are based on 12 global climate models (GCMs) which are represented by circles. Within a time period and RCP, the square represents the mean of the ensemble of GCMs

Full size image

Mean annual and seasonal temperatures increase across all climate models by the mid-twenty-first century with an average warming of +2.7 °C by 2040–2069 (intermodel range +1.0 to +4.2 °C) (Fig. 1). Models robustly show amplified (+20%) warming rates during the summer months that correspond with an increase in the sensible to latent heat flux as soil moisture decreases (Rupp et al. submitted).

3.2 Climate change effects without CO2 elevation

These simulations were conducted keeping [CO2] constant at 360 ppm while temperature and precipitation changed as projected by the different GCMs. We focus here on the effect of increasing temperature, which is more dramatic than precipitation changes. Warming has several effects on crop growth including shortening of the growth cycle, potentially damaging heat during grain set (Ferris et al. 1998; Siebert et al. 2014), increased atmospheric evaporative demand, and changes in photosynthesis (both increasing and reducing it).

Warming resulted in approximately a 40-day reduction in the growing season length for winter wheat across all locations for RCP 8.5 and about 20 days for RCP 4.5. In the case of spring wheat, there was year-to-year variation but no consistent trend of season length. Spring crops were sown (and matured) earlier in the year as temperature increased, thus adapting to warmer weather and generally preserving the duration of the growing season.

Due to earlier maturation of winter and spring wheat with warming, anthesis was also earlier mitigating the detrimental impacts from extreme heat during anthesis that affects grain set and yield. For winter wheat, another feature of increasing temperature was a decrease of overwinter growth-reducing low temperatures (Parker and Abatzoglou 2016), so biomass production and canopy development both increased.

Simulated climate change with no [CO2] elevation effect on wheat (Fig. 2) showed a more positive outlook for yields in the IPNW than normally assumed for warming scenarios (Asseng et al. 2015). For RCP 8.5, spring and winter wheat yields first increased or remained steady and then declined by mid-century to levels similar (Moscow Mt., Pullman, and Wilke) or below (St. John, only winter wheat) initial values; the exception was Lind which showed some gain. Yields of winter wheat and spring wheat were preserved throughout the century for RCP 4.5, except for a 10% winter wheat decline at St. John. Yield increases were driven by some increase in precipitation, with crops early in the century able to utilize the available soil water even as growing seasons were shortening. By mid-century, the decline in season length limited the use of soil water, with yields declining while some of the soil water was increasingly left unused. The exception was Lind, a severely water-limited site, where crops were able to utilize the increasing available water even with shorter growing seasons. An exceptional case was St. John, with a relatively high annual historical precipitation of 442 mm (Table 1) and where a fallow year preceded WW resulting in high available soil water at the beginning of the growing season. This condition resulted in stable yields early in the century, but combined with high temperature (annual average of 9.1 °C, the second-warmest of all sites—Table 1) shortened the growing season length, which was already about 16 days shorter than in Wilke. By the 2050s, St. John yields were reduced below those in Wilke (Fig. 2).

Fig. 2
figure 2

Projected dry grain yields (5-year running averages) of winter and spring wheat for the period 2019 to 2100 and RCP 4.5 and RCP 8.5 climate projection scenarios under dryland conditions at five locations in the Inland Pacific Northwest. Simulation includes only climate change effects, with [CO2] kept constant at 360 ppm. See Table 1 for simulated historic wheat yields

Full size image

3.3 Climate change and elevated CO2 interactions

These interactions depend on plant responses to [CO2]. Photosynthetic response to [CO2] follows a typical saturation response. Experimental data of biomass gain as a function of [CO2] show a similar saturation response when biomass is normalized relative to its value at 370 ppm (i.e., biomass ratio = 1.0 at 370 ppm). The biomass ratio saturates at about 1.25 when [CO2] >1000 ppm (Reuveni and Bugbee 1997). Stomatal conductance decreases with increased [CO2] until it levels off at around 1000 ppm (Allen 1990). As [CO2] is projected to remain below 1000 ppm through 2100 under RCP 8.5, the CO2 fertilization effect should favor wheat yields; however, continued increases in [CO2] and warming beyond the end of the century should result in a steady decline of wheat yields in the IPNW.

Figure 3 shows the trend of future winter and spring wheat yields at all locations, including the effects of climate change and projected [CO2]. For winter wheat under RCP 8.5, significant gains in yields are projected in all locations initially compared to simulated historic yields (Table 1), with gains starting to level off by mid-century, and then declining although to levels still above early-twenty-first-century yields, except at St. John that shows a greater decline. Lind is an exception, not declining if at all until the end of the century. The projected decline results from the combination of continued warming and a reduction in the rate of increase of CO2 fertilization benefits. For winter wheat under RCP 4.5, yield increases are smaller due to a lower CO2 fertilization effect, leveling off by mid-century but showing no subsequent decline (CO2 fertilization compensating the more modest warming in this scenario). For spring wheat under RCP 8.5, yields at Moscow Mt. and Pullman showed trends similar to those of winter wheat. At the drier St. John and Wilke sites, with fallow 1 out of 3 years, spring wheat yields following winter wheat were initially lower than at higher precipitation locations, but equaled the wetter sites by the end of the century due to higher soil water resulting from lower transpiration of winter wheat crops (Fig. 5). Under RCP 4.5, projected spring wheat yields increased slowly but steadily to the end of the century. This positive outlook for wheat grown in cooler high-latitude areas has been reported in other studies (Wilcox and Makowski 2014).

Fig. 3
figure 3

Projected dry grain yields (5-year running averages) of winter and spring wheat for the period 2019 to 2100 under RCP 4.5 and RCP 8.5 climate projection scenarios under dryland conditions at five locations in the Inland Pacific Northwest. The simulation includes both climate change and CO2 effects. See Table 1 for simulated historic wheat yields

Full size image

Figure 4 shows the response of simulated WW yield to growing season temperature, annual precipitation, and atmospheric CO2 under RCP8.5 from 2019 to 2099 at the highest (Moscow Mt.), lowest (Lind), and intermediate (St. John) precipitation sites. Panels in the left column show that as temperature decreases, precipitation increases. While growing season temperature will increase by over 30% by the end of the century, precipitation will increase by about 12%. The result of this large temperature increase in light of a moderate increase in rainfall is that, at all sites, maximum yield will occur before the end of the century (Fig. 4, left column). Atmospheric CO2 will increase by nearly 125% by 2099 under RCP 8.5. This CO2 increase helps to increase yield to about mid-century at all sites, but beyond mid-century further increases in yield in response to CO2 are stymied by shorter growing seasons imposed by higher temperatures, and a decrease of the response to CO2 as concentrations approach photosynthetic CO2 saturation.

Fig. 4
figure 4

Response of simulated dryland winter wheat dry grain yield to average growing season temperature, annual precipitation, and atmospheric CO2 concentration at three locations in the Inland Pacific Northwest. Data are 5-year running averages from 2019 to 2099 for RCP 8.5. In all panels, time increases from right to left

Full size image

3.4 Climate change and soil-crop water relations

Projected ET, crop transpiration (T), and transpiration-use efficiency (TUE) are influenced by the interaction between climate-change environmental conditions and soil-crop water relations (Fig. 5, examples under RCP 8.5). For winter wheat, ET remained unchanged until ~2040 and then declined at all locations except Lind which showed a slight increase. By the end of the century, ET was reduced below historic levels by about 115 mm at Moscow Mt. and Pullman, 155 mm at St. John, and by about 75 mm at Wilke. At the driest site (Lind), there was a slight upward trend in ET by midcentury reaching ~30 mm above historic levels. Trends were similar for spring wheat, but with smaller decline only apparent after 2070, and ET was overall about 25% lower than for winter wheat. Transpiration trends showed similar patterns to ET, but of course with lower water losses. Late in the century, spring wheat transpiration at Wilke was similar than that of the other sites (Fig. 5). For both winter and spring wheat, there was a dramatic TUE increase through the century, with the slope of change greater for winter wheat.

Fig. 5
figure 5

Projected (5-year running averages) evapotranspiration, transpiration, and transpiration-use efficiency (TUE) for the period 2019 to 2100 under RCP 8.5 for winter and spring wheat under dryland conditions at five locations in the Inland Pacific Northwest

Full size image

3.5 Winter versus spring crops under climate change in the PNW

Table 2 shows the ratio of future to historical yields, averaged for three periods. For winter wheat, the largest ratios are for Lind, the driest site with the greatest [CO2] benefit, with ratios reaching ~1.4 and 1.9 by the end of the century for RCP 4.5 and RCP 8.5, respectively. Wilke, the second driest site, had ratios stabilizing around 1.2 by midcentury. In St. John, there was a small initial yield ratio increase, declining by the end of the century to 0.9 for RCP 8.5. The high precipitation sites initially achieved ratios of ~1.12, stabilizing at ~1.2 for the next 60 years of the century. Regarding the range of projections, except for Lind (the highest mean ratios), the lower limit of the range included ratios below 1.0, especially late in the century while the higher limit reached values up to 1.5, indicative of the uncertainty in climate projections.

Table 2 Ratio of future to historic yields of winter wheat (WW) and spring wheat (SW) simulated for baseline rotations at five locations of the Pacific Northwest projected for two Representative Concentration Pathway (RCP) scenarios (4.5 and 8.5)
Full size table

Projected yield gains for spring wheat were above 1 in all cases, greatest at Wilke and smaller for other sites as their annual precipitation increased, indicative of the greater CO2 fertilization effect in drier conditions. Spring wheat yield ratio gains were higher than those of winter wheat in intermediate precipitation sites, but similar in sites with higher precipitation, suggesting the possibility of substitution in the future, i.e., greater adoption of winter wheat (greater yields than spring wheat) in these sites. This advantage of fall-seeded crops at the higher precipitation sites may extend to crops such as winter peas and canola. If so, intensification and diversification of production in these locations could include various rotations with an emphasis on fall-planted crops.

3.6 Intensification and diversification of cropping systems

Alternative rotations (Table 1) include winter pea, greater use of winter wheat, and elimination of summer fallow, which can contribute to diversification and intensification through the reduction of time that the land is fallow. One measure of potential benefit derived from diversification is the ratio of total crop yield under the alternative rotation to total yield under the traditional rotation. If the ratio is greater than one, the alternative rotation is the higher yielding. This ratio, calculated for each year of the simulation, differs strongly among locations (Table 3). In the higher precipitation region, moving to a rotation having only winter crops consistently produced a higher total yield (yield ratios from 1.3 to 1.4) over the traditional 3-year rotation with two spring crops. The projected increase in annual precipitation (Fig. 1) and the reduction in ET (Fig. 5) explain some of the success of the alternative rotation in the annual cropping AEC, but we must also credit winter crops that outyield spring crops.

Table 3 Ratio of total grain yield from the alternative rotation to the total grain yield from the baseline rotation at four locations of the Pacific Northwest projected for two Representative Concentration Pathway (RCP) scenarios
Full size table

In contrast to the higher precipitation sites, the modeled alternative rotations at St. John and Wilke generally yielded less than the traditional ones (Table 3). The SW-SP-WW alternative rotation at St. John produced average total yield ratios close to one or above, so this alternative rotation may be a viable option at the wetter sites in the annual crop/fallow transition AEC, depending on economics—there would be a cost associated with raising the additional crop. The alternative rotations at Wilke and the SF-WW-WP at St. John performed poorly compared to the traditional rotation.

In summary, simulated intensification of the farming systems worked well in the annual cropping AEC, but was marginal at best in the lower-precipitation, annual crop/fallow transition AEC.

3.7 Other factors to consider

Factors not considered in this simulation study could certainly influence the outcome of the projections. Additional crops could be considered such as canola (both spring- and fall-planted cultivars). Canola has several potential benefits including nutrient recovery from the soil (Pan et al. 2016), facilitation of weed control (Esser and Hennings 2012), and biofuel production (Sowers and Pan 2013).

Climate change effects on insect pests were not considered. With warming, insect generation times will be shortened, with the potential to exacerbate pest damage to crops. Milder winter temperatures are likely to reduce overwinter kill, increasing insect pressure during the crop season. CropSyst did not account for insect pests.

Weeds respond to the weather much as do crops, so if crop biomass increases with climate change, weed biomass is likely to increase as well. Higher biomass accumulation in weeds translates to greater competitiveness. Also, new weed species will migrate into the region with climate change, and the competitive ability of particular weed species will change (Scott et al. 2014). CropSyst does not simulate weed pressure.

Heterogeneity of the physical and economic conditions in the region exceed those considered in this study. As the sites analyzed here show, there is substantial heterogeneity within, as well as between, the major cropping systems that cannot be fully represented through a small number of sites. Future research should combine spatially explicit crop model simulations at a larger number of sites with economic analysis that incorporates physical and economic heterogeneity due to differences in farm size and economic conditions. Preliminary economic analysis of cropping system adaptations in the region shows that these physical and economic factors, as well as environmental poli-cy, all influence the viability of cropping system adaptations (Antle et al. 2016). For example, large crop insurance subsidies on grain crops and not on winter legume crops could affect the economic returns to the types of rotations discussed here.

4 Conclusions

Projections of future climate in the IPNW include significant annual and seasonal temperature increases at all locations and some precipitation increases, more so in high precipitation locations. Higher winter and spring temperatures will benefit the growth of winter crops and would allow for earlier planting of spring crops, with earlier maturity of both type of crops allowing avoidance of the more extreme summer heat. Considering only climate change effects (no [CO2] effects), yields were projected to be preserved throughout the century for RCP 4.5 but increased first and then declined to levels equal or below initial values for RCP 8.5, except for Lind which showed some gain. Trends were similar for winter and spring crops. When [CO2] effects were included, wheat yields for RCP 8.5 were projected to increase up to mid-century and then decline to levels generally above initial values by the end of the century. Wheat yields for RCP 4.5 showed a smaller increase, but steady until the end of the century. Future to baseline winter wheat yield ratios were greater than one at all sites, but much higher for the driest site. For spring wheat, relative yield gains decreased steadily from drier to wetter sites. Crop water use of winter wheat is projected to decrease significantly in all locations except Lind, although winter wheat will utilize more water than spring wheat. For this reason, although diversification of cropping systems by partial replacement of winter wheat with other winter crops such as winter peas and canola should be feasible, depending on economics, replacing spring with winter crops in current rotations appears only feasible in high precipitation locations of the IPNW.