Global Relationships between Cropland Intensification and Summer Temperature Extremes over the Last 50 Years

Nathaniel D. Mueller Department of Earth System Science, University of California, Irvine, Irvine, California

Search for other papers by Nathaniel D. Mueller in
Current site
Google Scholar
PubMed
Close
,
Andrew Rhines Department of Atmospheric Sciences, University of Washington, Seattle, Washington

Search for other papers by Andrew Rhines in
Current site
Google Scholar
PubMed
Close
,
Ethan E. Butler Department of Forest Resources, University of Minnesota, Twin Cities, St. Paul, Minnesota

Search for other papers by Ethan E. Butler in
Current site
Google Scholar
PubMed
Close
,
Deepak K. Ray Institute on the Environment, University of Minnesota, Twin Cities, St. Paul, Minnesota

Search for other papers by Deepak K. Ray in
Current site
Google Scholar
PubMed
Close
,
Stefan Siebert Institute of Crop Science and Resource Conservation, University of Bonn, Bonn, Germany

Search for other papers by Stefan Siebert in
Current site
Google Scholar
PubMed
Close
,
N. Michele Holbrook Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts

Search for other papers by N. Michele Holbrook in
Current site
Google Scholar
PubMed
Close
, and
Peter Huybers Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts

Search for other papers by Peter Huybers in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

Conversion of native ecosystems to cropland and the use of irrigation are considered dominant pathways through which agricultural land-use change alters regional climate. Recent research proposes that increases in cropland productivity, or intensification, also influences climate through increasing evapotranspiration. Increases in evapotranspiration are expected to have the greatest temperature influence on extremely hot summer days with high vapor pressure deficits. Here, the generalizability and importance of such relationships are assessed by examining historical land-use and climate trends in seven regions across the globe, each containing a major temperate or subtropical cropping area. Trends in summer high-temperature extremes are sequentially compared against trends in cropland area, area equipped for irrigation, precipitation, and summer cropping intensity. Trends in temperature extremes are estimated using quantile regression of weather station observations, and land-use data are from agricultural inventories and remote sensing. Intensification is the best predictor of trends in extreme temperatures among the factors that are considered and is generally associated with trends that are 0.2°–0.4°C decade−1 cooler than in adjacent regions. Neither cropland area nor precipitation trends are systematically associated with extreme temperature trends across regions, although high temperatures are suppressed over those portions of central North America and East Asia experiencing growth in irrigation. Both the temperature trends associated with intensification and increased irrigation can be understood as a consequence of increased latent cooling. These results underscore that the weather experienced by crops is not entirely external but also depends on agricultural practices.

Denotes content that is immediately available upon publication as open access.

© 2017 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Nathaniel D. Mueller, nmueller@uci.edu

Abstract

Conversion of native ecosystems to cropland and the use of irrigation are considered dominant pathways through which agricultural land-use change alters regional climate. Recent research proposes that increases in cropland productivity, or intensification, also influences climate through increasing evapotranspiration. Increases in evapotranspiration are expected to have the greatest temperature influence on extremely hot summer days with high vapor pressure deficits. Here, the generalizability and importance of such relationships are assessed by examining historical land-use and climate trends in seven regions across the globe, each containing a major temperate or subtropical cropping area. Trends in summer high-temperature extremes are sequentially compared against trends in cropland area, area equipped for irrigation, precipitation, and summer cropping intensity. Trends in temperature extremes are estimated using quantile regression of weather station observations, and land-use data are from agricultural inventories and remote sensing. Intensification is the best predictor of trends in extreme temperatures among the factors that are considered and is generally associated with trends that are 0.2°–0.4°C decade−1 cooler than in adjacent regions. Neither cropland area nor precipitation trends are systematically associated with extreme temperature trends across regions, although high temperatures are suppressed over those portions of central North America and East Asia experiencing growth in irrigation. Both the temperature trends associated with intensification and increased irrigation can be understood as a consequence of increased latent cooling. These results underscore that the weather experienced by crops is not entirely external but also depends on agricultural practices.

Denotes content that is immediately available upon publication as open access.

© 2017 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Nathaniel D. Mueller, nmueller@uci.edu

1. Introduction

Climate is a central determinant of crop distribution and productivity, yet climate itself can be influenced by agricultural land use and land cover via biophysical changes to surface albedo, rates of evapotranspiration, and surface roughness (Foley et al. 2003; Brovkin et al. 2004; Feddema et al. 2005; Diffenbaugh 2009; Pielke et al. 2011). Conversion of native ecosystems to cropland and the use of irrigation have long been considered dominant pathways through which agricultural land use alters regional temperatures. In the United States, cropland expansion altered albedo and evapotranspiration patterns and is thought to have cooled growing season temperatures (Bonan 1999, 2001; Oleson et al. 2004; Twine et al. 2004; Diffenbaugh 2009). Irrigation increases evapotranspiration and decreases temperatures, a relationship that has been documented in the U.S. Great Plains (Adegoke et al. 2003; Mahmood et al. 2006; Bonfils and Lobell 2007; Lobell et al. 2008; Harding and Snyder 2012; Lu et al. 2015), the Central Valley of California (Bonfils and Lobell 2007), Sudan (Alter et al. 2015b), and Asia (Bonfils and Lobell 2007). More recently, other changes to cropland management have been shown to alter climate. Multiple-cropping practices influence the seasonality of evapotranspiration in the North China Plain (Jeong et al. 2014) and the Brazilian Cerrado region (Spera et al. 2016) and are associated with higher temperatures during the intercropping period (Jeong et al. 2014). No-till practices can increase postharvest albedo, and model simulations suggest that increased adoption of no till on winter season crops in western Europe could substantially cool summer extreme temperatures (Davin et al. 2014).

Another recently proposed pathway by which agricultural land use can influence climate is through the intensification of crop production on existing croplands and an associated increase in evapotranspiration. Mueller et al. (2016) demonstrated century-long cooling trends in the U.S. Midwest that were proportional to trends in intensification documented in crop survey data, where intensification was defined as a positive trend in local crop biomass production. Cooling was observed for both irrigated and rainfed croplands that have undergone intensification, but with the important caveat that temperatures revert to historically high magnitudes during drought conditions in rainfed regions. These results are broadly consistent with studies of climatic trends for cropland in the Canadian Prairies (Gameda et al. 2007; Betts et al. 2013), where it was found that summer maximum temperatures decreased over the past several decades. Gameda et al. (2007) and Betts et al. (2013) attributed this pattern to greater landscape productivity and evapotranspiration as a result of declines in summer fallow practices, although the U.S. Midwest findings (Mueller et al. 2016) suggest that increased productivity on planted areas also contributed to changes in evapotranspiration across the Canadian Prairies.

In addition to observational evidence from historical data, the expectation that higher-productivity landscapes exhibit greater evapotranspiration accords with a number of field-scale studies. Vegetation productivity is tightly coupled to rates of evapotranspiration, and vegetation mediates the relationship between surface energy fluxes and soil moisture (Williams and Torn 2015). High-nitrogen application has been shown to result in both a larger magnitude (Jones et al. 1986; Rudnick and Irmak 2014) and duration (Rudnick and Irmak 2014) of peak evapotranspiration in maize. Nitrogen stress can otherwise be an important control on evapotranspiration through inhibiting leaf area, stomatal conductance, and root development (Jones et al. 1986; Chapin et al. 1988) but is largely alleviated in high-intensity cropping systems. Some crops are now managed at much greater planting densities (Duvick 2005), a change that can also lead to greater rates of evapotranspiration (Jiang et al. 2014). Adoption of conservation tillage practices, common in the United States (Horowitz et al. 2010), suppresses soil evaporation early in the season and thus can conserve water for transpiration (Gallaher 1977). Changes in cultivars may also influence transpiration characteristics, as more recent cultivars tend to have higher rates of stomatal conductance and lower canopy temperatures (Fischer et al. 1998; Barker et al. 2005; Roche 2015).

Given that the pace of cropland expansion has been relatively slow since 1950 (Ramankutty and Foley 1999) and that widespread increases in crop productivity occurred during this time period due to the adoption of “green revolution” technologies and management practices (Tilman et al. 2002), intensification of existing croplands may now be a dominant mechanism through which agricultural practices change regional climate. However, this relationship has only been documented in the U.S. Midwest (Mueller et al. 2016), an area that exhibits the most pronounced peak summer vegetation growth of anywhere on the planet (Guanter et al. 2014; Mueller et al. 2016). It is unclear whether more modest increases in crop productivity would significantly influence high-temperature trends elsewhere, and variability in cropping practices, soils, and atmospheric conditions also raise questions about the geographic generalizability of the U.S. Midwest intensification–cooling relationship. Examination of other regions provides an opportunity to test whether intensification is systematically related to a suppression of high temperatures.

Here, we examine the relationship between extremely hot maximum temperatures and summer cropland intensification, as well as the relative importance of intensification alongside changes in cropland area, irrigation growth, and precipitation, by analyzing land-use and extreme temperature trends for seven regions across the globe (Fig. 1). The management (Mueller et al. 2012; Mueller and Binder 2015; Siebert et al. 2015), productivity (Monfreda et al. 2008; Ray et al. 2012, 2013), and phenology (Sacks et al. 2010; Guanter et al. 2014) of crops vary widely across regions, providing a useful series of case studies to examine land–atmosphere connections with observational data. The analysis is restricted to subtropical and temperate regions because of greater availability of high-quality weather station records and the presence of well-defined seasonality in extreme temperatures and evaporative demand. We focus on summer as the season when evaporative demand is greatest and when temperature extremes generally have the greatest societal consequences, although crop damages from extreme heat will depend upon the specific timing of the exposure relative to sensitive periods of crop development (Gourdji et al. 2013; Butler and Huybers 2015). Consistent with Mueller et al. (2016), we examine the 95th percentile of summer daily maximum temperatures using quantile regression. Hot extremes exhibit unique trends relative to lower percentiles of the temperature distribution (McKinnon et al. 2016; Mueller et al. 2016) and are particularly sensitive to changes in evapotranspiration (Seneviratne et al. 2010; Mueller and Seneviratne 2012; Huybers et al. 2014; Mueller et al. 2016).

Fig. 1.
Fig. 1.

Regions examined for associations between agricultural land use, precipitation, and extreme temperatures are shown in orange boxes and include central North America, northern North America, western Europe, northern East Asia, southern East Asia, southern Australia, and southern South America. Within each region, a major cropping area is identified (in green), and these areas are used to characterize patterns of crop phenology within each region. Major cropping areas are defined as areas where the trend in our SCI index (defined in section 2d) is >1 gC m−2 summer−2, cropland area is >50% gridcell area, and gridcell centers are within the bounds identified by the dashed lines.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

2. Data and methods

The ability to document global-scale relationships between climatic trends and changes to summer cropping intensity, irrigation, and cropland area is only recently possible because of the release of several global historical land-use datasets used in coordination with weather station and satellite observations. Below, we detail our geographic areas of interest, the analysis of land-use trends, and the analysis of temperature and precipitation trends.

a. Regions and major cropping systems

Relationships between agricultural land-use and climate trends are examined across seven broad regions (orange lines in Fig. 1). We also identify grid cells comprising an intensified major cropping area in each region; these grid cells are utilized solely to characterize local crop phenology in a series of descriptive plots. To define these grid cells, we first delineate the most important continuous cropland regions (latitude and longitude boundaries are shown by the dashed lines in Fig. 1). Grid cells within these boundaries are then classified as a “major cropping area” if they contain greater than 50% cropland according to a circa-2000 dataset (Ramankutty et al. 2008) and exhibit positive trends in our summer cropping intensity index, defined below.

b. Cropland area trends

Historical cropland area is estimated from agricultural census records in combination with land-cover classifications from remote sensing Ramankutty and Foley (1999). The dataset has been recently updated (N. Ramankutty 2014, personal communication) and is now available at ½° resolution between 1961 and 2007. Trends are fit over this available interval using simple linear regression (Fig. 2a).

Fig. 2.
Fig. 2.

(a) Trends in cropland area for 1961–2007, (b) trends in AEI for 1961–2005, and (c) trends in summer precipitation for 1961–2014. Cropland area is from a historical dataset based on satellite and agricultural census data (Ramankutty and Foley 1999). AEI is determined from agricultural census and land-use records as recorded by Siebert et al. (2015). Precipitation data are from the GHCND weather station dataset, and dot sizes are scaled according to Voronoi polygons surrounding each station. Summer seasons are defined as JJA in the Northern Hemisphere and DJF in the Southern Hemisphere.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

c. Irrigated area trends

Data on area equipped for irrigation have been compiled by Siebert et al. (2015) into a gridded dataset at 5′ resolution covering the years 1900–2005, with maps available every 10 years from 1900 to 1980 and every 5 years after 1980 (available online at https://mygeohub.org/publications/8/2). This dataset is based on agricultural census information and detailed local land-use maps. We utilize the AEI–EarthStat–IR version of the GIS dataset that is constructed using the update to Ramankutty and Foley (1999) cropland areas and maximizes consistency with subnational irrigation statistics (denoted by “IR” in the version name). Trends are fit to gridcell area equipped for irrigation (AEI) values for 1961–2005 (Fig. 2b), where values for 1961 are calculated by linearly interpolating between 1960 and 1970 values in each grid cell. We fit trends at the native resolution of the irrigation dataset and all subsequent gridded data then upscale by averaging to ½° resolution so that all datasets are on a common grid.

d. Summer cropping intensity trends

To evaluate trends in summer cropping intensity (where a positive trend is considered “cropland intensification”), we develop an index of summer cropping intensity (SCI) that quantifies yearly summer crop biomass production across the landscape in units of grams of carbon per square meter produced over the summer growing season (gC m−2 summer−1). Yearly crop biomass production can be calculated from historical crop-specific harvested area and yield data, along with parameters that relate yields to total crop biomass. Unfortunately, these datasets do not also detail the seasonality of crop growth, a crucial consideration since changes to crop evapotranspiration characteristics only plausibly influence summer temperature extremes when crop growth occurs during the summer. Many temperate areas grow some crops during a “winter season,” when the crop is planted in the autumn and is harvested in the late spring or early summer, so a summer growing season cannot be assumed. In earlier work focused on the United States, Mueller et al. (2016) were able to isolate statistics for summer crop types, but this is not possible with the global crop datasets that we employ. To incorporate the seasonality of crop growth, we use remote sensing data to calculate the fraction of vegetation growth occurring during summer [vegetation summer fraction (VEGsf)]. We then utilize VEGsf as a fractional weight on crop biomass to convert annual cropping intensity to SCI. The crop datasets and calculations are described in greater detail below.

1) Calculating annual crop biomass production

To obtain trends in crop biomass production for six major crops, we first calculate the net primary productivity per harvested area (NPPha; in units of gC m−2 yr−1) of each crop from data on the yield (Y; converted to units of g m−2) of harvested crop products, as well as the dry fraction of the harvested product (DF; g g−1), the carbon content (C; gC g−1), the harvest index (HI; g g−1), and the aboveground fraction (AF; g g−1). Following Monfreda et al. (2008),
e1
where c is the crop type, y is the year, and i represents the index of each grid cell. We use gridded, crop-specific yield data spanning the years 1961–2008. Yield data for maize (grain, not silage), wheat, soybean, and rice are from Ray et al. (2012) and are generally resolved subnationally for major agricultural countries, although the temporal frequency of source data depends upon availability. Yield data for barley and rapeseed are from Monfreda et al. (2008) and are resolved subnationally for the year 2000. To obtain a historical time series, we scale these base maps to match the national-level average yield data from the Food and Agriculture Organization of the United Nations (FAO 2016), while preserving subnational spatial heterogeneity in yields from 2000. Values for DF, C, AF, and modern HI are directly from Monfreda et al. (2008). The harvest index of some crops has changed as a result of crop breeding, and historical values are reported in Table 1. In lieu of detailed data about the temporal evolution of HI, we assume a linear scaling between historical and modern values from 1910 to 1980, with modern values used for 1980 and all subsequent years. The use of historically varying HI values decreases the calculated intensification trend and works in opposition to the yield trends, but the latter are much larger and dominate the intensification trends. Our results are not sensitive to the use of historically varying harvest indices.
Table 1.

Historical and modern HI values by crop. All modern HI values are drawn from the compilation by Monfreda et al. (2008), and references for the historical values are listed.

Table 1.
Harvested area is relevant for considering the extent to which cropland evapotranspiration characteristics influence temperature. A large increase in evapotranspiration across a small field would, obviously, have limited influence on regional air temperatures. Thus, we multiply NPPha by the harvested area for each crop (HAc; in units of m2) relative to the total area within each grid cell (TA; m2), giving an area-normalized net primary productivity metric (NPPan),
e2

Harvested area data for our six crops are from the same sources (Monfreda et al. 2008; Ray et al. 2012; FAO 2016) as the yield data. The units for NPPan remain gC m−2 yr−1, although the area (m−2) is now relative to gridcell area and not harvested area. Trends in NPPan are fit for 1961–2008 (Fig. 3a) and provide a useful measure of cropland intensification for our six crops. However, as previously mentioned, these estimates do not indicate whether that intensification would have occurred during a summer growing season or at other portions of the year.

Fig. 3.
Fig. 3.

(a) Trends in NPPan over 1961–2014, calculated using harvested area and yield records for six major crops: maize, wheat, rice, soybean, barley, and rapeseed. (b) The VEGsf, calculated using SIF from the GOME-2 satellite instrument. (c) Trends in the SCI index, calculated by multiplying NPPan trends and VEGsf.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

2) Weighting by the vegetation summer fraction to calculate SCI

The Global Ozone Monitoring Experiment-2 (GOME-2) satellite instrument record of sun-induced chlorophyll fluorescence (SIF) (Joiner et al. 2013) is our preferred source of data for calculating VEGsf. These data are available at monthly 0.5° resolution. Chlorophyll fluorescence has previously been shown to exhibit closer correspondence with cropland gross primary productivity (GPP) from eddy flux towers than reflectance-based indices (Guanter et al. 2014). However, the relatively coarse resolution implies that the fluorescence data capture photosynthesis from both native and managed vegetation. This limitation is more pronounced for heterogeneous landscapes (e.g., western Europe) as opposed to those that are comparatively dominated by crops (e.g., the North American Corn Belt).

Using the SIF data, we calculate VEGsf. Assuming a simple linear scaling between SIF and GPP, the units for VEGsf are grams of carbon per summer divided by grams of carbon per year [gC summer−1 (gC yr−1)−1]. Summer is defined as June–August (JJA) in the Northern Hemisphere and December–February (DJF) in the Southern Hemisphere. Thus, for the Northern Hemisphere,
e3
where m is the month. Any negative SIF values, which do arise because of measurement errors, are set to zero prior to calculating VEGsf. We use the average summer fraction during the recent years of 2007–12 (Fig. 3b), and we test whether this fraction has varied over time using NDVI data as described below. Summer fraction is not calculated for areas with insufficient signal, here specified as monthly average fluorescence less than 1/12 mW m−2 sr−1 nm−1 (these areas are shown as light gray in Fig. 3b).
Our final SCI index is constructed by using VEGsf to weight NPPan and is calculated for all locations in the extratropics,
e4

Trends in the SCI index (Fig. 3c) retain the prominent NPPan trends in summer cropping areas (e.g., the U.S. Corn Belt and the Canadian Prairies), while NPPan trends in predominantly winter-cropping areas are down weighted (e.g., in portions of the U.S. southern Great Plains and southern Australia).

3) VEGsf sensitivity analysis

An alternate source of data for calculating VEGsf is the Global Inventory Monitoring and Modeling System (GIMMS) normalized difference vegetation index (NDVI) record generated from the Advanced Very High Resolution Radiometer (AVHRR) (Tucker 2014). These data are available bimonthly at 5′ resolution and span 1982–2013. Despite the aforementioned drawbacks of reflectance-based indices, these NDVI data permit an alternate estimation of SCI for comparison against our standard SIF approach. To permit for direct comparison against the SIF estimate, NDVI seasonality is computed over the 2007–12 interval. SCI is calculated at the 5′ resolution permitted by the NDVI data and then averaged to 0.5° resolution. The long temporal record also allows us to examine the extent to which VEGsf has changed over time, a topic we return to in section 3h.

e. Crop calendar data

Additional information about the seasonal cycle of crop development can be determined from global crop calendar data, and we use these data as contextual information for interpreting our findings. Average regional planting and harvest dates by crop type, as well as typical ranges around those means, are taken from the Sacks et al. (2010) database. These data do not include information about trends in planting and harvest dates as influenced by management practices and climate trends (e.g., Kucharik 2006). Spatial averages across major cropping regions are calculated for each crop type, where averages are weighted according to gridcell crop harvested areas (Monfreda et al. 2008). Planting and harvest dates for summer rapeseed in Canada are from USDA (1994), because Sacks et al. (2010) only contains data on winter rapeseed. We also determine crop harvested areas (Monfreda et al. 2008) circa 2000 as fractions of the total land area within each major cropping system. These values are shown in planting and harvest date figures to indicate the relative importance of various crops in each region. Planting and harvest data are presented alongside seasonal cycles of SIF for further context on local phenology in each major cropping area.

f. Climate trends

Weather station data are from the Global Historical Climatology Network–Daily (GHCND) dataset (Menne et al. 2012). Observations with negative quality flags are removed. In the interest of achieving a relatively complete geographic sample, we include any station where a minimum of 60% of days (after quality filtering) report values of maximum temperature from 1961 to 2014. All regions have average coverage considerably above this baseline, as shown in Table 2.

Table 2.

The percent of summer station days reporting maximum temperature observations across all weather stations, listed by region and time period. Summer is defined as JJA in the Northern Hemisphere and DJF in the Southern Hemisphere.

Table 2.

Quantile regression (Koenker and Bassett 1978) is utilized to assess trends in temperature extremes, and we focus on trends in the 95th percentile of daily maximum temperature observations during the summer months of June–August in the Northern Hemisphere and December–February in the Southern Hemisphere (Fig. 4). Temperature data were originally recorded in Fahrenheit and Celsius at different levels of precision and then were rounded to standard increments of 0.1°C for inclusion in GHCND. This heterogeneity poses problems for understanding trends in extreme temperatures since quantile regression assumes continuously distributed data and is biased by rounding artifacts. We correct for the effects of rounding by adding an appropriate amount of jitter to each observation to approximately correct each temperature record to its unrounded distribution, where jitter amplitude is determined from the results of a precision-decoding algorithm (Rhines et al. 2015).

Fig. 4.
Fig. 4.

Quantile regression trends in 95th-percentile summer daily maximum temperatures from 1961 to 2014. Temperature data are from the GHCND weather station dataset, and dot sizes are scaled according to Voronoi polygons surrounding each station. Summer seasons are defined as JJA in the Northern Hemisphere and DJF in the Southern Hemisphere.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Although daily temperature observations are the most suitable record for directly examining large-scale changes in extreme temperatures, station data are subject to a number of uncertainties. Station moves, changes in the time of observation, and shifts in equipment can all influence temperature observations (Quayle et al. 1991; Pielke et al. 2007b; Menne and Williams 2010). Examination of trends in temperature extremes in North America, using the same GHCND data and quantile regression approach, shows consistency between neighboring stations as well as between stations and reanalyses (Rhines et al. 2017), indicating that the influence of inhomogeneities in the daily temperature data are minor relative to trends in extreme temperature. Furthermore, pairwise comparison of summer temperature trends calculated from GHCND and from nearby hourly stations sampled using a consistent time-of-day window indicate that time-of-observation biases are small compared with typical magnitudes of summer temperature trends (McKinnon et al. 2016). Within the United States, the widespread change in thermometers during the 1980s is thought to have introduced a cool bias to maximum temperatures of around 0.4°C (Quayle et al. 1991). We suggest that these inhomogeneities and uncertainties in the data, while important for understanding the absolute magnitude of temperature trends, will have less influence on our identification of land-use effects given our focus on spatial differences in temperature trends. Moreover, the extent to which results are consistent between countries with different weather station networks serves as an important check on the robustness of our results.

Trends in precipitation are analyzed for the same subset of stations used to examine temperature trends. Precipitation can influence extreme temperatures through the influence of soil moisture availability on evapotranspiration (Mueller and Seneviratne 2012) and can also be affected by land-use change (Pielke et al. 2007a; DeAngelis et al. 2010; Harding and Snyder 2012; Alter et al. 2015a,b; Mueller et al. 2016). The relationship between precipitation and evapotranspiration is modulated by the ability of vegetation to access stored soil moisture in the root zone, which generally acts to suppress the impacts of precipitation anomalies on evapotranspiration (Betts et al. 2014). Average precipitation per day is calculated by season and year, and from these averages, seasonal total precipitation is estimated for every year where at least 80% of daily observations are present. Trends are then calculated for seasonal total precipitation using simple linear regression for every station where at least 80% of the seasonal totals are present (Fig. 2c).

The land area most closely associated with each weather station is calculated using spherical Voronoi polygons (Renka 1997). For coastal stations that fall just outside of our coastal boundaries, a minimum area of 1 ha is associated with the station. Station area is used to calculate the widths of boxes in our box-and-whisker plot figures and to scale the dot sizes associated with weather station locations on figures showing temperature and precipitation trends.

g. Statistical analysis

A bootstrap test is utilized to assess the significance of 95th-percentile temperature trends for weather stations experiencing a given shift in precipitation or land use relative to stations experiencing little change in that explanatory variable. Groupings of stations by land use and precipitation are shown in subsequent box-and-whisker plots for each region. The test accounts for spatial autocorrelation by resampling all station observations identically and accounts for temporal autocorrelation by resampling 3-month seasonal blocks. For each bootstrap replicate (1000), 95th-percentile temperature trends are fit to the resampled data at each station using quantile regression. We then take the difference in the mean trend of stations experiencing a given shift in land use or precipitation and the mean trend of stations experiencing no change in that explanatory variable. This procedure generates a distribution of mean differences that is compared with zero to determine a two-sided p value. The test is similar to the approach taken in Mueller et al. (2016), although that analysis was with respect to whether temperature trends grouped by a given explanatory variable were significantly different than zero, whereas here we evaluate if temperature trends are significantly different from adjacent areas that have little change in the explanatory variable.

h. Case study

An example illustrating the temporal resolution of the land-use and climate data employed in this study is presented in Fig. 5 for Redwood County, Minnesota. Maize and soybean are the dominant crops in the area, and both crops exhibit increasing yields since 1960 (Figs. 5a,b). Increases in maize and soybean harvested area (Fig. 5a) have been at the expense of other crops, with total cropland area remaining roughly constant (Fig. 5e). Cropland area represents all land devoted to crops and therefore tends to be more stable than harvested areas of individual crops, which can be affected by changing market conditions and weather-induced crop failure (e.g., note the drop in maize harvested area during the flood of 1993). Area equipped for irrigation is negligible (Fig. 5e). Summer precipitation shows substantial interannual variability and a modest long-term trend of 7 mm decade−1 (Fig. 5e, regression line not shown).

Fig. 5.
Fig. 5.

An example showing local crop and land-use characteristics, weather data, and 95th-percentile maximum temperature trends for a weather station in Redwood County, Minnesota. (a) Crop harvested areas and (b) crop yields for all crops (of the six considered) where the maximum harvested area was greater than 1% of gridcell area. (c) The VEGsf as calculated using SIF and NDVI. (d) NPPan and SCI, calculated using crop harvested area, crop yield, and SIF-based VEGsf according to Eqs. (1)(4). (e) Cropland area, area equipped for irrigation, and summer (JJA) precipitation are also considered as predictors of changing extreme temperatures. (f) Daily summer maximum temperature observations, with the 95th-percentile quantile regression trend overlaid in maroon. The quantile regression trend is calculated after adding jitter to the observations to account for rounding artifacts. (g) A histogram of 95th-percentile maximum temperature trends derived from a block-bootstrap resampling of yearly observations. The trend line fit using all the data is shown in the thick maroon line, and dashed lines indicate the 95% confidence interval on the trend. All land-use data are extracted for the nearest grid cell to the weather station, and gridded data are used at the original resolution of each dataset (5′ for the crop harvested area and yield data, 5′ for the irrigation data, and ½° for the cropland area data).

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Yield and harvested area data are combined according to Eqs. (1) and (2) to calculate NPPan (Fig. 5d) and linearly scaled into SCI using SIF-determined summer fraction of photosynthesis (VEGsf) according to Eqs. (3) and (4). The approximately linear increase in SCI over time reflects increases in yield and greater land devoted to high-yielding maize and soybean crops (Fig. 5c). Variations in crop types, crop productivity, planting schedules, or weather could all cause the summer fraction of SIF to vary with year. Although disaggregating the reasons for variations in satellite-based estimates of VEGsf is beyond the scope of this paper, it is possible to examine the summer fraction as a function of year back to 2007 using SIF and 1981 using NDVI. Both products show interannual variability, but neither exhibit strong trends. VEGsf calculated using SIF is systematically higher than when calculated using NDVI, an expected pattern since SIF tracks GPP more closely than NDVI (Guanter et al. 2014).

The distribution of summer temperatures is indicated in Fig. 5f, where the size of dots indicate the frequency of temperature observations during the summer months, binned to the nearest 0.5°C for legibility. Quantile regression of the 95th-percentile temperature shows a decreasing temperature trend of −0.3°C decade−1. A block bootstrap of the daily temperature data is used to assess significance of the temporal trend. For each bootstrap replicate, years are sampled with replacement, and all summer temperature observations are used for every year sampled. Quantile regression trends are fit to the sampled data for 1000 bootstrap replicates. The distribution of 95th-percentile temperature trends from the bootstrap demonstrates that this trend significantly differs from zero at 95% confidence (Fig. 5g). Note that although this calculation of significance applies to the temporal trend for this individual station, our calculation of significance related to land-use and precipitation trends depends upon relative temperature trends between weather stations grouped according to various explanatory variables.

Four different predictor variables are considered candidates for explaining the observed trends in 95th-percentile temperatures: total cropland area, area equipped for irrigation, precipitation, and SCI. We consider the explanatory power of each of these variables by examining the regionwide associations between temperature trends and trends in each predictor variable. In Redwood County, we see that each variable other than SCI displays minor trends since 1960. When examining regionwide associations between the predictor variables and temperature trends, the Redwood County weather station would therefore be included in the control group of stations (see box-and-whisker plots below) for both trends in cropland area and trends in area equipped for irrigation. For the precipitation analysis, the station would be grouped with other stations with modest positive trends. SCI is the only predictor variable with a strong positive trend that co-occurs with the significant cooling in summer 95th-percentile temperatures.

3. Results and discussion

Trends in 95th-percentile summer maximum temperatures are systematically cooler over intensified croplands relative to neighboring areas. This relationship holds in every region where summer cropping is the dominant land use, including for central North America, northern North America, northern East Asia, southern East Asia, and southern South America. Median trends in 95th-percentile maximum temperatures are 0.2°–0.4°C decade−1 cooler in intensifying areas compared to adjacent areas not experiencing intensification. No relationship is found in western Europe and southern Australia, areas where winter cropping dominates. Consistent with earlier work (Mueller et al. 2016), cooling is found in rainfed areas, such as the Canadian Prairies and much of the North American Corn Belt, as well as in irrigated areas. Substantial irrigation growth has occurred in East Asia, helping facilitate increases in cropland productivity. Therefore, both irrigated area trends and summer intensification trends are related to cooler temperature extremes in these areas. Changes in cropland area and precipitation are generally weak predictors of trends in extreme temperatures.

In each region discussed below, the relationship between weather station 95th-percentile temperature trends and local trends in our candidate predictor variables is described, discussed in the context of the literature, and presented visually using a series of box-and-whisker plots. Candidate predictor variables are the local trends in cropland area, area equipped for irrigation, summer cropping intensity, and precipitation (from the same weather station). All trends in predictor variables are calculated using simple linear regression (sections 2bd and 2f). In each plot, weather stations are evenly binned into subsets of stations according to local trends in the predictor variables. Subsetting allows us to examine how temperature trends vary with trends in the predictors in a way that is independent of functional form and provides the basis for the aforementioned bootstrap test. Each box-and-whisker plot displays the full range of 95th-percentile temperature trends for a given subset of weather stations, with asterisks indicating the significance of the temperature trends.

a. Central North America

Cropland intensification is strongly associated with cooling in the central North America region (Fig. 6), which covers most of the continental United States and southeastern Canada. These results are consistent with earlier results identifying an association between intensification and cooling from 1910 to 2014 using USDA crop survey data of twelve summer crop types (Mueller et al. 2016). Trends in 95th-percentile temperatures (Fig. 4) tend to show cooling or absence of warming over intensified cropland areas, while much of the rest of the region shows warming of around 0.1°C decade−1.

Fig. 6.
Fig. 6.

Trends in central North American temperature extremes grouped according to candidate predictor variables: (a) cropland area, (b) AEI, (c) summer precipitation, and (d) SCI. Data points are from weather stations that have been associated with local (nearest ½° grid box) trends in land-use characteristics. Weather stations are evenly binned according to land-use or precipitation trends. Box-and-whisker plots display the full range of temperature trends across stations for each bin, with the boxes containing the interquartile range, whiskers extending up to 1.5 times the interquartile range, and crosses indicating outliers beyond this range. Asterisks indicate that 95th-percentile temperature trends for a given bin significantly differ from those in the control group (gray box, centered on zero trend) at the p < 0.05 level or p < 0.01 for double asterisks. The x-axis values are generally the midpoints of each bin, although edge bins include weather stations associated with outlier trends in each explanatory variable. Box widths are proportional to the area associated with the constituent weather stations except for the control bins that are narrowed by a factor of 5 for legibility. Box colors are consistent with the maps in Figs. 2 and 3.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

The dominant crops within the North American Corn Belt are maize (accounting for 25% of total area) and soybeans (24%), and their growth is centered on summer months (Fig. 7). Average values of summer SIF across the world’s midlatitudes are around 0.8 mW m−2 sr−1 nm−1, but in central North America they peak in July at values exceeding 3 mW m−2 sr−1 nm−1, the largest values found for any spatially extensive region on the globe (Guanter et al. 2014; Mueller et al. 2016). The anomalously high productivity of the region is reflected in values of SCI that are higher than any other major cropping area (Table 3). We infer that achieving these high rates of photosynthesis during the summer season has led to corresponding increases in evapotranspiration. This inference is supported by estimates of a positive evapotranspiration trend over the Mississippi River basin (Milly and Dunne 2001) and is consistent with trends toward greater specific and relative humidity during summer in regions of intensified crop growth (Sandstrom et al. 2004; Brown and DeGaetano 2013). Further, we note that climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) simulate temperature increases over the central United States in response to historical forcings (Kumar et al. 2013), further emphasizing the importance of mechanisms not included in the models to explain historical temperature trends.

Fig. 7.
Fig. 7.

Seasonal patterns of vegetative development for the major crop production areas of the central North American Corn Belt. (a) Median monthly SIF and the interquartile range of monthly values calculated across available years. (b) Average crop seasons—from planting to harvest—for major crops in the region from data compiled by Sacks et al. (2010). Ranges of typical planting and harvest dates are indicated with the dashed black lines. Harvested area of major crops (Monfreda et al. 2008) in each region are indicated next to crop names and are used to scale the width of the boxes devoted to each crop. Given that two seasons of wheat are present, bar area is divided equally between the two categories since crop harvested area data are not separated by season. Both SIF and crop season data are weighted spatial averages across those grid cells indicated for the central North America region in Fig. 1, where weights are cropland area from Ramankutty et al. (2008) for the SIF plot and individual crop harvested areas from Monfreda et al. (2008) for the crop season plot.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Table 3.

Average cropland area, AEI, VEGsf calculated from chlorophyll fluorescence data, SCI index, and summer precipitation for major cropping areas. The major cropping areas are defined by the green grid cells in Fig. 1. Each average is calculated over the full temporal range of the data, from 1961 to 2007 for cropland area, from 1961 to 2005 for AEI, from 1961 to 2008 for SCI, and from 1961 to 2014 for precipitation. VEGsf is calculated over the recent years of 2007–12. No precipitation data are shown for the Argentine Pampas owing to data limitations.

Table 3.

Extreme temperatures since 1961 have cooled most strongly over the western Corn Belt, an area of substantial land-use change and expanding commodity crop production (Lark et al. 2015). The stronger cooling over this area may arise from more influential land-use transitions or from the gradual reduction in aerosol forcing over eastern North America. The cooling influence of aerosols on temperatures is thought to have peaked during the 1970s–90s; therefore, reductions in forcing would contribute to a warming trend that may counteract the influence from intensification (Leibensperger et al. 2012a,b). Since the climate of the western Corn Belt was never as strongly influenced by aerosols, this may explain the stronger cooling observed in this area.

Increasing area equipped for irrigation is found across the Great Plains and in rice-growing areas adjacent to the Mississippi River. Those locations where area equipped for irrigation increased 2.5%–3.5% of gridcell area per decade show significant cooling of 95th-percentile summer temperatures (p < 0.05) relative to regions with near-constant irrigated area (Fig. 6b) and become yet more significant for decadal trends greater than 3.5%. However, the amount of cooling area associated with increasing irrigation is only 14% of that associated with intensification, emphasizing that increases in vegetation productivity influence evapotranspiration characteristics, whether in irrigated or rainfed areas. Area calculations are performed using Voronoi polygons associated with each weather station exhibiting negative 95th-percentile temperature trends and associated with either increases in irrigated area of at least 2.5% per decade (Fig. 6b) or intensification trends of at least 0.5 gC m−2 summer−2 (Fig. 6d).

Trends in cropland area are inconsistently related to 95th-percentile temperature trends (Fig. 6a). The appearance of significant cooling in relation to 2% decade−1 growth in cropland area may reflect greater evapotranspiration from cropland expansion but also may result from the fact that we test candidate mechanisms in isolation. The presence of extreme temperature trends primarily driven by changes in irrigation and intensification makes it more likely that a random subsetting of the region can contain temperature trends that are larger than that of the control group. In future work, a multifactor panel analysis would likely prove a better indicator of exact significance.

Weather stations with decreased precipitation have slightly higher extreme temperature trends than other stations, which would be consistent with the effects of lower soil moisture, decreased evapotranspiration, and greater sensible heating from the land surface (Fig. 6c). However, the warming relationship is not significant for all subsets of stations with decreasing precipitation, and stations with increasing precipitation do not exhibit significant cooling. In contrast, Mueller et al. (2016) found a significant relationship between precipitation increases and cooler temperatures in the midwestern United States in their study of trends over the last century. They noted that such trends may be partly due to cropland intensification (Mueller et al. 2016) or irrigation growth across the Great Plains (DeAngelis et al. 2010; Harding and Snyder 2012; Alter et al. 2015a), since precipitation in the region is strongly influenced by rates of evapotranspiration (Betts 2004). The present analysis focused on trends since 1961 shows some areas of increasing precipitation in the region (Fig. 2c) but no significant relationship between cooling and elevated precipitation.

b. Northern North America

Northern North America also shows cooling of 95th-percentile temperatures associated with cropland intensification (Fig. 8). Irrigation growth has been minimal and shows no strong relationship with the pattern of temperature trends. Crop phenology in the Canadian Prairies is strongly summer seasonal but with a shorter growing season than in the Corn Belt.

Fig. 8.
Fig. 8.

As in Figs. 6 and 7, but for northern North America. One outlier station where the 95th-percentile summer temperature trend was >2°C decade−1 has been removed from the box-and-whisker plots and statistical analysis. Phenology is shown in (e) and (f) for the major crop production areas of the Canadian Prairies.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Our findings align with earlier studies that identified a cooling of maximum temperatures and an increase in relative humidity during the period of peak crop growth in the Canadian Prairies (Gameda et al. 2007; Betts et al. 2013, 2016). This cooling was attributed to greater landscape evapotranspiration from declining cropland area left fallow during summer (Betts et al. 2013). Increased productivity on planted areas and declining summer fallow will both influence SCI, respectively, via changes to crop yields and harvested areas. The SCI trend in heavily cropped areas is 1.4 gC m−2 summer−2 (Fig. 2a). To distinguish harvested area and yield contributions to this trend, we recalculate SCI holding each fixed (Fig. 9). SCI trends from harvested area variations alone give a trend of 0.5 gC m−2 summer−2. Conversely, SCI trends are 0.9 gC m−2 summer−2 when only yields are allowed to vary.

Fig. 9.
Fig. 9.

Disaggregating contributions to SCI trends in the Canadian Prairies. (a) Trends in SCI calculated using yearly varying harvested area and average crop yields over 1961–2008. (b) Trends in SCI calculated using yearly varying yields and average harvested area. Note that the scale is truncated relative to Fig. 3 to better highlight differences between the calculations.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Insomuch as summer cooling is linearly proportional to SCI trends, which is far from clear but appears the simplest assumption, increasing productivity on planted areas (determined from the yield trends) is the dominant influence on cooling. However, we note that we find greater increases in SCI from changing harvested area in Alberta and Manitoba than in Saskatchewan, despite inventory data showing the greatest declines in fallow for Saskatchewan (Betts et al. 2013). This discrepancy may result from expansion of harvested area unrelated to declining summer fallow, crop types not included in our analysis, or local-scale changes that we do not resolve in our historical crop data. We are particularly limited in resolving spatial patterns of change for barley and rapeseed since our area and yield time series are generated by perturbing circa-2000 maps with national-level data. A more complete analysis of influences on temperature would be possible by utilizing higher-resolution data on agricultural practices and by running regional climate simulations with fallow and productivity scenarios.

Gameda et al. (2007) and Betts et al. (2013) found increases in precipitation associated with elevated evapotranspiration during peak crop growth, indicating greater precipitation recycling and increased potential for deep convection triggered by land management shifts (Raddatz 1998). We also find positive precipitation trends over the Canadian Prairies (Fig. 2c), lending support to this notion. However, the associations between 95th-percentile temperature trends and precipitation trends over the whole northern North America region are more ambiguous. Areas with greater precipitation do not systematically show significantly cooler temperatures. However, most stations experiencing drying trends do have significantly elevated warming trends relative to the control group, consistent with decreases in evapotranspiration and increases in sensible heating. Greater temperature sensitivity to decreases in precipitation than to increases in precipitation is consistent with the results of Betts et al. (2017) for the Canadian Prairies, where it was demonstrated that the diurnal temperature range in the region exhibits greater coupling with precipitation anomalies during dry conditions than during wet conditions.

c. Western Europe

Intensification is not associated with cooling in western Europe (Fig. 10). The 95th-percentile temperature trends since 1961 show strong warming averaging 0.4°C decade−1 and have insignificant relationships with cropland area, irrigation, and SCI trends. Temperature trends appear to decline with increasing precipitation trends, but this relationship is insignificant and weak compared to the predictor relationships found elsewhere.

Fig. 10.
Fig. 10.

As in Figs. 6 and 7, but for western Europe. Phenology is shown in (e) and (f) for the major crop production areas of southern England and northwestern France.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

These negative findings appear to result from the dominance of winter cropping and the heterogeneity of the landscape. SIF peaks during May when the growing season for barley, rapeseed, and winter wheat varieties all coincide. Of the crops examined, only maize has a long summer season where peak transpiration and peak temperatures would align. Grain maize only covers 3% of the landscape within the heavily cropped areas of southern England and northwestern France. Silage maize for fodder is not included in our dataset, but judging from disaggregated maize area for France, including both would still only double this percentage (FAO 2016). For comparison, summer maize and soybean account for 49% of the total land area in the central North American Corn Belt (Fig. 6).

Moderate SIF values persist throughout the summer and give higher VEGsf values in western Europe compared with more homogenous winter cropping areas such as Kansas and southwestern Australia (Fig. 3b and Table 3). The resolution of the SIF input to VEGsf makes it difficult to separate this heterogeneous landscape into cropland and natural vegetation, leading to VEGsf values that are likely higher than would be observed on croplands alone.

The general warming in western Europe is therefore consistent with our hypothesis that intensification of summer crop production is associated with cooling. Given the low extent of summer cropping, the large majority of croplands are mature or harvested by late summer. The dominance of winter cropping systems affords the possibility of mitigating extremely hot temperatures by transitioning to no-till systems, which have increased summer surface albedo relative to tilled soil (Davin et al. 2014).

d. Northern East Asia

Intensification of summer crops coincides with suppressed extreme temperature trends in the northern East Asia region, which encompasses northern China; Mongolia; Hokkaido, Japan; and eastern Russia, with a southern boundary of 40°N, or roughly the latitude of Beijing, China. The major cropping area within this region is northeastern China, where summer cropping of maize, soybeans, and rice dominate the landscape. Warming of 95th-percentile temperatures at rates of around 0.2°C decade−1 is found in most of the region (Fig. 4), with the exception of an arc of near-zero warming extending north to south across northeastern China exhibiting strong trends in SCI and area equipped for irrigation (Fig. 11). The spatial patterns of the SCI trend and the irrigation trend are highly correlated because of the heavy reliance upon irrigation to facilitate increases in crop productivity and paddy rice production. Areas of northeastern China, where intensification and irrigation trends are strong, exhibit both increasing and decreasing area devoted to cropland. If we consider intensification and irrigation the primary drivers of cooling, this spatial overlap explains the counterintuitive finding that both increasing and decreasing cropland area trends are associated with cooler extreme temperature trends. Precipitation trends exhibit no consistent association with extreme temperature trends.

Fig. 11.
Fig. 11.

As in Figs. 6 and 7, but for northern East Asia. Phenology is shown in (e) and (f) for the major crop production areas of northeastern China.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Our results are consistent with several recent studies suggesting land use has cooled summer temperatures in northeastern China. Hu et al. (2010) compare surface temperature observations to reanalysis products that do not include land-use forcing—the “observation minus reanalysis” methodology—in order to estimate the influence of land-use change. Similar to our results, they find cooling in maximum temperatures in northeastern China relative to reanalysis. Cao et al. (2015) force a regional climate model with remotely sensed changes in biophysical land surface parameters, including increases in leaf area index and vegetated fraction, and find cooling in cropped areas between 2001 and 2010. Zhao et al. (2016) find cooling and wetting trends from 1960 to 2014 associated with cultivated land fraction, with May–September daily maximum temperature trends in heavily cultivated areas 0.10°C decade−1 cooler than areas with minimal cropland.

A major uncertainty is the climatic influence of aerosol emissions and tropospheric ozone across Asia (Liao et al. 2015). While black carbon emissions and tropospheric ozone lead to warming, other pollutants are expected to have a cooling effect on surface temperatures. One set of model simulations (Chang et al. 2009) indicates that, on net, these forcings have minimal influence on summer temperatures but cause cooling during the winter months. However, Du et al. (2017) use an observationally based attribution methodology to suggest suppression of average warm season air temperature trends in northeastern China are due to declines in surface solar radiation. Detailed modeling studies are needed to understand the relative contributions of land-use change, air pollution, and greenhouse gases on temperature trends.

e. Southern East Asia

Cropland intensification is associated with cooling in the southern East Asia region, which includes areas of China, the Korea Peninsula, and Japan from south of 40°N to the Tropic of Cancer. Warming in 95th-percentile temperatures of around 0.2°C is seen over most of the region, with the exception of cooling over the major cropping area of the North China Plain and an absence of major warming extending south from this region through central China (Fig. 4). The pattern of changes in temperature reflects that of SCI (Fig. 12). Areas of negative SCI trends in South Korea and Japan are associated with the greatest rates of warming, whereas intensified landscapes in the North China Plain exhibit the most cooling. Similar to northern East Asia, cropland intensification across much of this region is accompanied and supported by increases in irrigation, such that trends in the area equipped for irrigation are also significantly associated with reductions in 95th-percentile temperatures. Area equipped for irrigation is higher in the North China Plain than any other major cropping area examined (Table 3).

Fig. 12.
Fig. 12.

As in Figs. 6 and 7, but for southern East Asia. Phenology is shown in (e) and (f) for the major crop production areas of the North China Plain.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Our results for southern East Asia are consistent with the land-use influence identified in the analysis of observations and reanalysis by Hu et al. (2010) and the regional modeling of Cao et al. (2015). Bonfils and Lobell (2007) have also identified cooling of irrigated areas relative to surrounding unirrigated land in this region. Given that much of the heavily cultivated areas have experienced declines in cropland area while increasing productivity, decreases in cropland area are associated with reductions in extreme temperature trends. Precipitation trends appear unrelated to temperature trends. Aerosol emissions and tropospheric ozone are likely also important in this region. Although one modeling study indicates minimal net influence of pollutants on summer temperatures (Chang et al. 2009), other research points to a suppression of warm season air temperature trends in the North China Plain of over 0.1°C attributed to changes in surface solar radiation (Du et al. 2017). Since changes in evapotranspiration from land can also influence cloudiness, modeling studies exploring the interactions between pollution and land-use change are necessary.

The major cropping area in this region is the North China Plain, an area where much of the land is double cropped with winter wheat (Figs. 12e,f). The intercropping period is centered on June according to the SIF data, and a large peak in photosynthetic activity occurs during July and August, corresponding to growth of the second crop. These findings suggest that elevated evapotranspiration rates associated with intensification of the second crop are sufficient to contribute to a cooling of 95th-percentile temperatures over the 3-month summer season. Jeong et al. (2014) note that temperatures during the intercropping period in double-cropped areas are higher than in areas planted with a single crop as a result of lower rates of evapotranspiration.

f. Southern Australia

In extratropical Australia, no substantial correlation exists between any of our explanatory variables and patterns of warming (Fig. 13). The null result for intensification is to be expected given that winter wheat is dominant for the intensified production area in southwest Australia. Winter seasonality is clearly demonstrated in the annual cycle of SIF and in the planting and harvest data. As a result, no significant variation exists in SCI. It is possible we would find associations between extreme temperatures and intensification if we extended our analysis to the winter growing season, as previous work focused on the wheat lands of Western Australia found elevated latent heat fluxes during the winter growing season over cropped areas relative to neighboring natural vegetation (Ray et al. 2003).

Fig. 13.
Fig. 13.

As in Figs. 6 and 7, but for southern Australia. Phenology is shown in (e) and (f) for the major crop production areas of Western Australia.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

g. Southern South America

Data availability is limited in southern South America (Fig. 14); however, several stations overlap with intensified cropland area in the Argentine Pampas west of Buenos Aires, Argentina (Fig. 3). Consistent with expectations, those stations that have positive SCI trends all exhibit 95th-percentile temperature trends that are negative or indistinguishable from zero, while the average 95th-percentile temperature trend across all other areas is toward warming. Strong relationships are not observed between 95th-percentile temperature trends and other predictors. Precipitation records in this region have a high number of missing observations, limiting our ability to analyze associations between temperature and precipitation trends. Our results are consistent with those of Nuñez et al. (2008), who find cooling of maximum temperatures and diurnal temperature range over the Pampas using an observation minus reanalysis approach. These authors also analyze precipitation trends using a more complete network of stations, finding elevated precipitation cooccurring with areas of cooling. Crop phenology in the Argentine Pampas is a mix of winter wheat and summer crops. Soybeans are the most dominant crop, and the area planted to soybeans has expanded substantially in recent years (Nuñez et al. 2008).

Fig. 14.
Fig. 14.

As in Figs. 6 and 7, but for southern South America. Phenology is shown in (d) and (e) for the major crop production areas of the Argentine Pampas.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

h. Vegetation seasonality from NDVI data

Global patterns of vegetation seasonality remain similar when calculating VEGsf using the GIMMS NDVI data instead of GOME-2 SIF data for the years 2007–12; however, the magnitudes of NDVI-based VEGsf tend to be slightly lower (Fig. 15a) than the SIF-based values since reflectance-based indices do not track the seasonality of vegetation growth as tightly as SIF (Guanter et al. 2014). Consistent associations are seen between SCI, calculated using NDVI-based VEGsf (SCI-NDVI), and summer temperature trends (Fig. 16).

Fig. 15.
Fig. 15.

(a) VEGsf calculated using the GIMMS NDVI data over 2007–12, consistent with the calculation for SIF. (b) The decadal trend in VEGsf calculated using GIMMS NDVI data over 1982–2013. Areas where VEGsf was not calculated using the SIF data are masked.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Fig. 16.
Fig. 16.

(a) The SCI index calculated using GIMMS NDVI data instead of SIF to calculate the VEGsf (SCI-NDVI). Associations between SCI-NDVI and 95th-percentile summer temperature trends for (b) central North America, (c) northern North America, (d) northern East Asia, (e) southern East Asia, (f) and southern South America.

Citation: Journal of Climate 30, 18; 10.1175/JCLI-D-17-0096.1

Trends in VEGsf using NDVI over 1982–2013 (Fig. 15b) show weak, but positive, trends over the western Corn Belt, the Canadian Prairies, and the Argentine Pampas. Positive trends presumably reflect cropland intensification, soybean expansion in Argentina, and declining summer fallow in Canada. Negative trends in the North China Plain could be the result of increased double cropping (Ray and Foley 2013; Gray et al. 2014a; Jeong et al. 2014). If SCI could be calculated with yearly varying VEGsf over the full record, the VEGsf trend analysis suggests that the magnitudes of SCI would be slightly higher in many cropped regions, with the exception of the North China Plain. However, the spatial patterns of intensified (high SCI trend) versus nonintensified (zero or low SCI trend) areas would likely be minimally affected, suggesting little bearing on our conclusions.

4. Conclusions

A significant relationship between intensification and cooler temperature extremes is found across all regions with substantial trends toward intensified summer cropping. Intensification is consistently the strongest land-use predictor of extreme temperature trends and is associated with cooling in both rainfed and irrigated cropping systems. In portions of central North America and East Asia, growth in area equipped for irrigation is also closely related to cooling. Median 95th-percentile temperature trends in intensified areas are systematically 0.2°–0.4°C decade−1 lower than in neighboring areas not experiencing intensification. Cooling associated with both intensification and increased irrigation can be understood as a consequence of increased latent cooling associated with elevated rates of evapotranspiration. Regional cooling can thus be added to the list of impacts associated with cropland intensification, alongside land demand (Matson and Vitousek 2006; Burney et al. 2010), nutrient application (Vitousek et al. 2009), the seasonality of atmospheric carbon dioxide (Gray et al. 2014b), water use (Siebert and Döll 2010), and water quality (Diaz and Rosenberg 2008).

Because extreme high temperatures are associated with crop damages, their amelioration by enhanced evapotranspiration raises the interesting question of how much of the agricultural intensification that we estimate, which is largely driven by improvement in yield, can be characterized as a positive feedback. There are, however, a number of extenuating circumstances regarding the operation of such a feedback. Cooling from evapotranspiration in rainfed areas is lost during drought conditions, leading to greater temperature shocks when soil moisture is depleted (Mueller et al. 2016). Increased soil water consumption could also increase crop exposure to dry spells unless water is recycled through increased rainfall. Further, although extreme temperatures may be reduced over the summer months, temperature trends during key early season reproductive periods are often toward warming (Gourdji et al. 2013). Higher atmospheric CO2 concentrations increase plant water use efficiency (Leakey et al. 2009), a change that may offset some of the otherwise expected increases in evapotranspiration. Also of note is that increased humidity levels may lead to little net change in heat index extremes for local human populations despite cooler air temperatures (Lobell et al. 2008).

Suppression of extreme temperatures by high-intensity croplands can be considered a climate regulation service (West et al. 2011), but the total climatic influence of any ecosystem is a function of both biophysical and biogeochemical climate forcings. On an annual basis, the modeling and accounting performed by Anderson-Teixeira et al. (2012) indicate U.S. croplands and grasslands have similar climate regulation values, driven by high rates of evapotranspiration in cropland and high carbon storage in grasslands.

Further analyses are needed to understand the contribution of intensification-driven amelioration of temperature extremes on historical and future crop productivity. Crop yield models typically treat temperatures as an exogenous driver of productivity, although crop development and productivity play an important role in modifying surface energy fluxes (Williams and Torn 2015) and temperature extremes (Mueller et al. 2016). Moreover, the cooling effect of evapotranspiration on crop canopy temperature is much larger than the cooling effect on air temperature measured at standard weather stations (Siebert et al. 2014), and only recently has systematic modeling of canopy temperature been introduced into crop models to better reflect the impact of transpiration-driven cooling on crop heat stress (Webber et al. 2017). The degree to which management practices alter local weather and climate may have first-order implications for future yield trends.

Acknowledgments

Funding provided by USDA AFRI Fellowship 2016-67012-25208 to NDM and NSF Hydrological Sciences Grant 1521210. We thank Alan Betts and two anonymous reviewers for helpful comments and feedback, and we thank Marena Lin, Karen McKinnon, and Martin Tingley for helpful conversations.

REFERENCES

  • Adegoke, J. O., R. A. Pielke Sr., J. Eastman, R. Mahmood, and K. G. Hubbard, 2003: Impact of irrigation on midsummer surface fluxes and temperature under dry synoptic conditions: A regional atmospheric model study of the U.S. high plains. Mon. Wea. Rev., 131, 556564, doi:10.1175/1520-0493(2003)131<0556:IOIOMS>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Alter, R. E., Y. Fan, B. R. Lintner, and C. P. Weaver, 2015a: Observational evidence that Great Plains irrigation has enhanced summer precipitation intensity and totals in the midwestern United States. J. Hydrometeor., 16, 17171735, doi:10.1175/JHM-D-14-0115.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Alter, R. E., E.-S. Im, and E. A. B. Eltahir, 2015b: Rainfall consistently enhanced around the Gezira scheme in East Africa due to irrigation. Nat. Geosci., 8, 763767, doi:10.1038/ngeo2514.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Anderson-Teixeira, K. J., P. K. Snyder, T. E. Twine, S. V. Cuadra, M. H. Costa, and E. H. DeLucia, 2012: Climate-regulation services of natural and agricultural ecoregions of the Americas. Nat. Climate Change, 2, 177181, doi:10.1038/nclimate1346.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Barker, T., and Coauthors, 2005: Improving drought tolerance in maize. Plant Breeding Reviews, J. Janick, Ed., Vol. 25, John Wiley and Sons, 173–253, doi:10.1002/9780470650301.ch7.

    • Crossref
    • Export Citation
  • Betts, A. K., 2004: Understanding hydrometeorology using global models. Bull. Amer. Meteor. Soc., 85, 16731688, doi:10.1175/BAMS-85-11-1673.

  • Betts, A. K., R. Desjardins, D. Worth, and D. Cerkowniak, 2013: Impact of land use change on the diurnal cycle climate of the Canadian prairies. J. Geophys. Res. Atmos., 118, 11 99612 011, doi:10.1002/2013JD020717.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Betts, A. K., R. Desjardins, D. Worth, and B. Beckage, 2014: Climate coupling between temperature, humidity, precipitation, and cloud cover over the Canadian prairies. J. Geophys. Res. Atmos., 119, 13 30513 326, doi:10.1002/2014JD022511.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Betts, A. K., R. Desjardins, and D. Worth, 2016: The impact of clouds, land use and snow cover on climate in the Canadian prairies. Adv. Sci. Res., 13, 3742, doi:10.5194/asr-13-37-2016.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Betts, A. K., A. B. Tawfik, and R. L. Desjardins, 2017: Revisiting hydrometeorology using cloud and climate observations. J. Hydrometeor., 18, 939955, doi:10.1175/JHM-D-16-0203.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bonan, G. B., 1999: Frost followed the plow: Impacts of deforestation on the climate of the United States. Ecol. Appl., 9, 13051315, doi:10.1890/1051-0761(1999)009[1305:FFTPIO]2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bonan, G. B., 2001: Observational evidence for reduction of daily maximum temperature by croplands in the midwest United States. J. Climate, 14, 24302442, doi:10.1175/1520-0442(2001)014<2430:OEFROD>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bonfils, C., and D. Lobell, 2007: Empirical evidence for a recent slowdown in irrigation-induced cooling. Proc. Natl. Acad. Sci. USA, 104, 13 58213 587, doi:10.1073/pnas.0700144104.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brovkin, V., S. Sitch, W. von Bloh, M. Claussen, E. Bauer, and W. Cramer, 2004: Role of land cover changes for atmospheric CO2 increase and climate change during the last 150 years. Global Change Biol., 10, 12531266, doi:10.1111/j.1365-2486.2004.00812.x.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Brown, P. J., and A. T. DeGaetano, 2013: Trends in U.S. surface humidity, 1930–2010. J. Appl. Meteor. Climatol., 52, 147163, doi:10.1175/JAMC-D-12-035.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Burney, J., S. J. Davis, and D. B. Lobell, 2010: Greenhouse gas mitigation by agricultural intensification. Proc. Natl. Acad. Sci. USA, 107, 12 05212057, doi:10.1073/pnas.0914216107.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Butler, E. E., and P. Huybers, 2015: Variations in the sensitivity of US maize yield to extreme temperatures by region and growth phase. Environ. Res. Lett., 10, 034009, doi:10.1088/1748-9326/10/3/034009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Cao, Q., D. Yu, M. Georgescu, Z. Han, and J. Wu, 2015: Impacts of land use and land cover change on regional climate: A case study in the agro-pastoral transitional zone of China. Environ. Res. Lett., 10, 124025, doi:10.1088/1748-9326/10/12/124025.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chang, W., H. Liao, and H. Wang, 2009: Climate responses to direct radiative forcing of anthropogenic aerosols, tropospheric ozone, and long-lived greenhouse gases in eastern China over 1951–2000. Adv. Atmos. Sci., 26, 748762, doi:10.1007/s00376-009-9032-4.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chapin, F. S., III, C. H. Walter, and D. T. Clarkson, 1988: Growth response of barley and tomato to nitrogen stress and its control by abscisic acid, water relations and photosynthesis. Planta, 173, 352366, doi:10.1007/BF00401022.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Davin, E. L., S. I. Seneviratne, P. Ciais, A. Olioso, and T. Want, 2014: Preferential cooling of hot extremes from cropland albedo management. Proc. Natl. Acad. Sci. USA, 111, 97579761, doi:10.1073/pnas.1317323111.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • DeAngelis, A., F. Dominguez, Y. Fan, A. Robock, M. D. Kustu, and D. Robinson, 2010: Evidence of enhanced precipitation due to irrigation over the Great Plains of the United States. J. Geophys. Res., 115, D15115, doi:10.1029/2010JD013892.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Diaz, R. J., and R. Rosenberg, 2008: Spreading dead zones and consequences for marine ecosystems. Science, 321, 926929, doi:10.1126/science.1156401.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Diffenbaugh, N. S., 2009: Influence of modern land cover on the climate of the United States. Climate Dyn., 33, 945958, doi:10.1007/s00382-009-0566-z.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Du, J., K. Wang, J. Wang, and Q. Ma, 2017: Contributions of surface solar radiation and precipitation to the spatiotemporal patterns of surface and air warming in China from 1960 to 2003. Atmos. Chem. Phys., 17, 49314944, doi:10.5194/acp-17-4931-2017.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Duvick, D. N., 2005: Genetic progress in yield of United States maize (Zea mays L.). Maydica, 50, 193202.

  • FAO, 2016: FAOSTAT. Food and Agriculture Organization of the United Nations. [Available online at http://www.fao.org/faostat/en/#data/QC.]

  • Feddema, J. J., K. W. Oleson, G. B. Bonan, L. O. Mearns, L. E. Buja, G. A. Meehl, and W. M. Washington, 2005: The importance of land-cover change in simulating future climates. Science, 310, 16741678, doi:10.1126/science.1118160.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fischer, R. A., D. Rees, K. D. Sayre, Z.-M. Lu, A. G. Condon, and A. L. Saavedra, 1998: Wheat yield progress associated with higher stomatal conductance and photosynthetic rate, and cooler canopies. Crop Sci., 38, 14671475, doi:10.2135/cropsci1998.0011183X003800060011x.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Foley, J. A., M. Costa, C. Delire, N. Ramankutty, and P. K. Snyder, 2003: Green surprise? How terrestrial ecosystems could affect earth’s climate. Front. Ecol. Environ., 1, 3844, doi:10.1890/1540-9295(2003)001[0038:GSHTEC]2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Gallaher, R. N., 1977: Soil moisture conservation and yield of crops no-till planted in rye. Soil Sci. Soc. Amer. J., 41, 145147, doi:10.2136/sssaj1977.03615995004100010040x.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gameda, S., B. Qian, C. A. Campbell, and R. L. Desjardins, 2007: Climatic trends associated with summerfallow in the Canadian prairies. Agric. For. Meteor., 142, 170185, doi:10.1016/j.agrformet.2006.03.026.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gourdji, S. M., A. M. Sibley, and D. B. Lobell, 2013: Global crop exposure to critical high temperatures in the reproductive period: Historical trends and future projections. Environ. Res. Lett., 8, 024041, doi:10.1088/1748-9326/8/2/024041.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gray, J., M. Friedl, S. Frolking, N. Ramankutty, A. Nelson, and M. K. Gumma, 2014a: Mapping Asian cropping intensity with MODIS. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens., 7, 33733379, doi:10.1109/JSTARS.2014.2344630.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gray, J., S. Frolking, E. A. Kort, D. K. Ray, C. J. Kucharik, N. Ramankutty, and M. A. Friedl, 2014b: Direct human influence on atmospheric CO2 seasonality from increased cropland productivity. Nature, 515, 398401, doi:10.1038/nature13957.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Guanter, L., and Coauthors, 2014: Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence. Proc. Natl. Acad. Sci. USA, 111, E1327E1333, doi:10.1073/pnas.1320008111.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Harding, K. J., and P. K. Snyder, 2012: Modeling the atmospheric response to irrigation in the Great Plains. Part I: General impacts on precipitation and the energy budget. J. Hydrometeor., 13, 16671686, doi:10.1175/JHM-D-11-098.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hay, R., 1995: Harvest index: A review of its use in plant breeding and crop physiology. Ann. Appl. Biol., 126, 197216, doi:10.1111/j.1744-7348.1995.tb05015.x.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Horowitz, J., R. Ebel, and K. Ueda, 2010: “No-till” farming is a growing practice. USDA Tech. Rep., 22 pp.

  • Hu, Y., W. Dong, and Y. He, 2010: Impact of land surface forcings on mean and extreme temperature in eastern China. J. Geophys. Res., 115, D19117, doi:10.1029/2009JD013368.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Huybers, P., K. A. McKinnon, A. Rhines, and M. Tingley, 2014: U.S. daily temperatures: The meaning of extremes in the context of nonnormality. J. Climate, 27, 73687384, doi:10.1175/JCLI-D-14-00216.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Jeong, S.-J., C.-H. Ho, S. Piao, J. Kim, P. Ciais, Y.-B. Lee, J.-G. Jhun, and S. K. Park, 2014: Effects of double cropping on summer climate of the North China Plain and neighbouring regions. Nat. Climate Change, 4, 615619, doi:10.1038/nclimate2266.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Jiang, X., S. Kang, L. Tong, F. Li, D. Li, R. Ding, and R. Qiu, 2014: Crop coefficient and evapotranspiration of grain maize modified by planting density in an arid region of northwest China. Agric. Water Manage., 142, 135143, doi:10.1016/j.agwat.2014.05.006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Joiner, J., and Coauthors, 2013: Global monitoring of terrestrial chlorophyll fluorescence from moderate spectral resolution near-infrared satellite measurements: Methodology, simulations, and application to GOME-2. Atmos. Meas. Tech., 6, 28032823, doi:10.5194/amt-6-2803-2013.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Jones, J. W., B. Zur, and J. M. Bennett, 1986: Interactive effects of water and nitrogen stresses on carbon and water vapor exchange of corn canopies. Agric. For. Meteor., 38, 113126, doi:10.1016/0168-1923(86)90053-5.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Koenker, R., and G. Bassett Jr., 1978: Regression quantiles. Econometrica, 46, 3350, doi:10.2307/1913643.

  • Kucharik, C. J., 2006: A multidecadal trend of earlier corn planting in the central USA. Agron. J., 98, 15441547, doi:10.2134/agronj2006.0156.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kumar, S., J.