08.17.2020

Adaptation in U.S. Corn Belt increases resistance to soil carbon loss with climate change

By Yao Zhang et al.

Abstract

Increasing the amount of soil organic carbon (SOC) has agronomic benefits and the potential to mitigate climate change. Previous regional predictions of SOC trends under climate change often ignore or do not explicitly consider the effect of crop adaptation (i.e., changing planting dates and varieties). We used the DayCent biogeochemical model to examine the effect of adaptation on SOC for corn and soybean production in the U.S. Corn Belt using climate data from three models. Without adaptation, yields of both corn and soybean tended to decrease and the decomposition of SOC tended to increase leading to a loss of SOC with climate change compared to a baseline scenario with no climate change. With adaptation, the model predicted a substantially higher crop yield. The increase in yields and associated carbon input to the SOC pool counteracted the increased decomposition in the adaptation scenarios, leading to similar SOC stocks under different climate change scenarios. Consequently, we found that crop management adaptation to changing climatic conditions strengthen agroecosystem resistance to SOC loss. However, there are differences spatially in SOC trends. The northern part of the region is likely to gain SOC while the southern part of the region is predicted to lose SOC.

Introduction

Soil organic matter is essential for maintaining soil health and sustaining plant growth. Loss of soil organic matter often leads to degradation of soil quality1. It also constitutes the largest terrestrial organic carbon (C) pool (~ 2,400 Pg C in top 2 m soil). This soil organic C (SOC) pool is three times greater than the amount of C in the atmosphere2. An increase or decrease of C in soils by only a small percent represents a substantial C sink or source for atmospheric CO2. Studies have been conducted to predict SOC of croplands under climate change3,4,5. However, crop management changes with adaptation to future climate6,7,8,9 (such as changes in varieties and planting dates) were ignored in most studies. Thus, we carried out a regional study to assess the impact of projected climate change and elevated CO2 on SOC in agricultural systems with management adaptation.

A change in SOC is a result of the net effect of the changes in SOC decomposition rates (the main C outflow) and C inputs from plants10 (main C inflow). In a warmer climate, the higher temperature could increase the decomposition rate in both the short and long term11. Decomposition will also be sensitive to increases or decreases in precipitation10 predicted by climate models8. Carbon input is mainly from root and surface residue litter that is not removed from the system through harvest, grazing or burning. Changes in temperature, precipitation, and atmospheric CO2 all affect this input through plant growth and production6. These climatic changes are also likely to affect management practices in the future7. Although systemic management changes are possible that can limit the impact of climate change, such as moving from dryland to irrigated systems, smaller adjustments including changes in varieties and planting dates have less barriers for adoption. These adjustments to management have been found to affect both crop production (C input) and decomposition7,12, with potentially larger effects on crop production7. These interactions are complex, and the overall change to SOC pools in agricultural lands remains uncertain.

Here we present the results from a simulation of SOC dynamics for more than 54,000 locations, covering 34,600,000 ha of cropland in the U.S. Corn Belt (Supplementary Fig. 1a). These locations have historically been managed with a corn (Zea mays L.) and soybean (Glycine max L.) rotation, which is the most common crop rotation in this region. Corn and soybean yields in this region account for about 85% of US crop production13. In our simulation, the widely used biogeochemical model DayCent14,15 was driven by weather data from the Representative Concentration Pathway 4.5 (RCP 4.5) climate change scenario8. Predictions from three General Circulation Models (GCMs; GFDL-CM3, MIROC-ESM, and MRI-CGCM3) were downscaled to generate daily weather (32 km grid) and used to assess the uncertainty. The main goal was to quantify the change in surface soil C (0–20 cm) under management adaptation with the selection of alternative crop varieties and planting dates to maintain high levels of crop production in future climate scenarios from 2041 to 2071. Alternative crop varieties were based on those currently available in the United States without consideration of additional crop breeding or genetic modifications that could further enhance production in the future. A no climate change scenario (historical weather with current CO2 level) was used as the baseline for comparisons.

Crop production trends

We found that, without adaptation, the average corn yield during the 2041–2071 period in the Corn Belt would drop in all three GCM climate scenarios in comparison with the baseline with no climate change (Fig. 1). Yield decreased by 17% (GFDL-CM3), 34% (MIROC-ESM), and 2% (MRI-CGCM3) under the respective scenarios. The difference between GCMs can be attributed to differences in the projected temperature and precipitation (Supplementary Fig. 2). Similar yield decreases were seen for soybean as found for corn if there is no adaptation (decreased by 13% for GFDL-CM3, 28% for MIROC-ESM, and increased by 8% for MRI-CGCM3).

figure1
Figure 1

With adaptation, both corn and soybean yields increased compared with no adaptation: an increase of 5% for corn and 19% for soybean (average of the three GCMs) compared with the baseline. The larger increase in soybean yield was due to the C3 crop being more responsive to the CO2 fertilization than C4 crops16. The standard deviation of yields across the three GCMs in the adaptation scenarios was lower than that of the no adaptation scenarios (15% and 39% lower for corn and soybean respectively). Simulated crop yields were more stable under climate change with adaptation management compared to without adaptation. The stability in yields is because, without crop adaptation associated with the selection of alternative crop varieties, an increase in temperature shortens the growing period of the crop17,18. Each GCM predicts a longer growing season with warmer temperatures. When a longer-season maturity variety was simulated as an adaptation pathway, the extended growing period allows the crop to use the full window of optimal solar radiation. This leads to similar amounts of production from year to year between GCM climate scenarios. The spatial pattern of the yield changes was very different for the two crops (Supplementary Fig. 3) due to the difference in crop response to temperature, day length (photoperiod), and elevated CO2. Our overall predictions of crop yield change were generally in agreement with other studies9,17,19,20.

Carbon input trends

Crop yields are good indicators of total net primary production and are found to be proportional to the total C input to soils in the U.S. Corn Belt21 (i.e., crop yield to total biomass does not vary much geographically). Our simulations predicted the average input to be 3.7 Mg C ha−1 year−1 under the no climate change scenario for the 2041–2070 period (Supplementary Fig. 4). With climate change but no adaptation, all counties (an administrative subdivision of a state in the U.S.) had lower C input (GCMs ensemble mean) compared with the baseline. However, with adaptation, more than half of the counties (most counties in the northern part of the region) had higher C input than those of the baseline. Compared with the no adaptation scenario, all counties in the Corn Belt maintained higher C input with adaptation. The average C input across climate scenarios in the Corn Belt region was predicted to be 3.0 and 3.9 Mg C ha−1 (a change of − 19% and 5% from the baseline) with no adaptation and adaptation scenarios, respectively. The standard deviation of the adaptation scenarios across the three GCMs was 47% lower than without adaptation, suggesting a counteracting effect of adaptation to climate change.

Decomposition factor trend

In contrast to C input, the decomposition factor, which reflected the relative change in decomposition rate due to temperature and moisture effect (ranges from 0 to 1, with higher rates associated with values closer to 1), increased in all counties regardless of the adaptation (Supplementary Fig. 4). Larger decomposition factors can be explained by an increase in soil temperature and wetter soil conditions associated with climate change projections. Wetter soil conditions were a result of a reduced transpiration rate under high levels of atmospheric CO2 and increased precipitation projected by the three GCMs. The average decomposition factor (across three GCMs) was predicted to be 0.41, 0.49, and 0.47 for baseline, without adaptation, and with adaptation scenarios. The standard deviations were 27% lower with adaptation for the decomposition factor across the three GCMs than without adaptation. If adaptation does not occur under the GCM scenarios, the growing period of the crops was reduced due to global warming and resulted in less total transpiration and wetter soil conditions, which increased decomposition. With adaptation, the longer-season maturity variety continued to grow over a longer period and consumed more water, reducing soil moisture to lower levels, thus lowering decomposition (average annual evapotranspiration was 2.0–5.1% higher in the adaptation scenario compared with no adaptation). The soil moisture differences among GCMs were lower with crop adaptation.

Soil organic C trends

Without climate change (baseline scenario), the predicted sub-region SOC in the top soil ranged from less than 30 Mg C ha−1 to more than 80 Mg C ha−1(Fig. 2a). The highest levels of SOC were found in the western part of the region where soil clay content is high (Supplementary Fig. 1). The low levels of SOC found in Michigan and northern Indiana can be attributed to the soils with high sand content. In addition, tillage intensity varies among the counties22, contributing to differences in SOC levels among counties. With climate change, there were losses of SOC in almost all counties if there was no adaptation (Fig. 2b). This was the net result of decreased C inputs and increased decomposition rates (Supplementary Fig. 4). With adaptation, the northern part of the Corn Belt tended to gain SOC (Fig. 2c) due to increased C inputs (Supplementary Fig. 4), while the other areas tended to lose SOC. Compared with no adaptation, the adaptation scenario resulted in higher SOC stocks in all counties of the Corn Belt, which was also found in a modeling study of a site in Michigan9.

figure2
Figure 2

Over the period from 2041 to 2070, the average SOC in the Corn Belt reached a new equilibrium state in the baseline scenario with no climate change (Fig. 3). Without adaptation in the climate change scenarios, SOC decreased over time. Rapid loss of C was predicted for the 2041–2050 period with the rate decreasing gradually over the 2051–2070 period. With adaptation, two GCMs (GFDL-CM3 and MRI-CGCM3) predicted similar changes over time as the baseline, while the SOC values for the other GCM (MIROC-ESM) were slightly lower than the baseline.

figure3
Figure 3

These results show that adaptation with the selection of longer-season varieties can lead to a similar SOC level as the baseline and GCM scenarios under a moderate climate change projection (RCP 4.5). Although the predicted climate was very different for the three GCMs, the SOC levels with adaptation were similar and not very different from the baseline with no climate change. In contrast, without adaptation, SOC storage levels were farther apart from each other and the baseline. This indicates a strong resistance to the effects of climate change in agricultural systems in the Corn Belt region if there is the selection for alternative varieties and planting dates that are better adapted to the changing climatic conditions.

We found the variation (CV 2.89%) of the predicted total SOC stock across the GCMs in the adaptation scenarios was much smaller than the variation (CV 5.88%) in the no adaptation scenarios. This corresponds to reduced variation in C input and decomposition factors. Because the increased decomposition rate was compensated by the higher overall C input, the SOC stock maintained similar levels as the baseline without climate change.

Limitations

In this study, we did not evaluate the possibility of new technologies that could be developed in the future and influence production, decomposition and other variables influencing SOC levels, as these changes are difficult to predict. For example, new technologies associated with crop breeding and other genetic improvements, pest control, and other developments could increase C input and result in higher SOC stocks than our predictions. Management practices such as adding cover crops, which increases C input, may also further enhance SOC stocks23. In contrast, large areas of removal of corn stover for biofuel production (not considered in our simulation) could reduce the total SOC stock24. Although most SOC is concentrated in the surface soil, subsoil C has been found to respond to warming climates and affects SOC stocks25. Future research should also address subsoil C dynamics.

Conclusions

The Corn Belt is the core area of U.S. crop production, and we found that management adaptation, which is often ignored in SOC studies, results in strong resistance to SOC losses that would otherwise occur with climate change. This resistance could enhance food security through its impact on soil health, fertility, and water holding capacity26. It also indicates a very low risk of losing SOC in this region, which will limit future carbon emissions to the atmosphere and further carbon debt from anthropogenic management of these lands27. However, there are counties in the southern part of the Corn Belt that are predicted to lose carbon even with adaptation. With future technology advancements and other management changes that are not considered in our study, it may be possible to reduce losses of SOC throughout the entire Corn Belt and potentially create a net carbon sink for atmospheric CO2 in the future28.

Methods

The region of our simulation included 12 states in the U.S. Corn Belt (IA, IL, IN, KS, MI, MN, MO, ND, NE, OH, SD, and WI) (Supplementary Fig. 1).

Climate change scenarios and climate data

Future weather data were generated using three General Circulation Models (GCMs) simulating the Representative Concentration Pathway (RCP) 4.5 greenhouse gas concentration trajectory8. Historical weather data used in the baseline scenario were derived from the 1980–2010 AgMERRA historical daily weather data29 (resampled at 32 km spatial resolution). GCM projections were downscaled to the AgMERRA spatial resolution and daily temporal resolution using the statistical method referred as “delta approach”30. The three GCMs used are MRI-CGM331, MIROC-ESM32, and GFDL-CM333. Supplementary Fig. 2 showed the characteristics of the three GCMs of the 2041–2070 period.

DayCent model, data and model implementation

The biogeochemical DayCent model was used in our study. We simulated 7 scenarios of corn/soybean rotation for the period of 2011–2070, which were combinations of three climate change scenarios (three GCMs of RCP 4.5) and two adaptation scenarios. Additionally, we simulated a baseline with no climate change. In the adaptation scenario, crop planting dates and crop maturity groups for each location were changed every decade to optimize production given the climate change conditions from the GCMs. In the baseline scenario, we used historic weather data and a stable CO2 concentration of 389 ppm. In the adaptation scenario, CO2 concentration increased from 389 in 2010 to 524 by 20708. The DayCent model was used to project SOC dynamics for 54,912 survey point locations that were in a corn/soybean rotation across the Corn Belt between 1990 and 2010. U.S. county boundaries were used to represent counties and all the data were summarized at the sub-region scale. Soil information was derived from the Soil Survey Geographic Database (SSURGO)34. More detailed information about the model and the regional simulation can be found in the Supplementary Materials.

Management data were derived from various sources or predicted. U.S. county-level tillage data (grouped into full, reduced, and no-tillage) were from the 2008 Crop Residue Management Survey conducted by the Conservation Technology Information Center (CTIC)22. Nitrogen fertilization rates were calculated using the predicted crop yield potential from DayCent for each scenario and recommended application rates to meet the targeted yields.

Planting dates for the baseline and no adaptation scenario were from the U.S. Department of Agriculture National Agricultural Statistics Service (NASS)35. To simulate the varieties of crops, we divided the crops into maturity groups (Supplemental Table 1). The spatial distributions of the maturity groups were derived from information provided by seed companies. The planting dates and maturity group distributions for each crop in the adaptation scenario were calculated based on a regression with temperature and estimated frost dates. Grain harvests were scheduled at the end of the growing season after maturity for corn and soybeans.

To establish a base SOC level (initial values for SOC pools) for each location, we conducted a model intialization36 using cropping histories from a variety of literature and historical databases developed for the US national greenhouse gas inventory37. In the pre-2011 historical simulation period, the radiation use efficiency (RUE) parameter was adjusted every 10 years from 1950 through 2010 to fit the simulated regional yield to the NASS reported historical yield trends for corn and soybean to represent crop breeding effects on production. Fertilizer rate was assumed to linearly increase from 1950 to 1980. Reduced till and no-till management was assumed to start in 1981 and only full tillage was used before 1981 in the historical simulations. The same initial values for SOC were used in each of the projections.

Model parameterization and verification

The crop parameters used for each maturity group are shown in Supplementary Table S1. The parameters were either from published data or model calibration using the experimental sites (Supplementary Table 2). We evaluated the DayCent simulated crop yields with NASS reported data13,38. More information about model parameterization and verification can be found in the Supplementary Materials.

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on request. However, the NRI survey data and location information are not publicly accessible because the data contain confidential information as mandated by law, 7 USC 2276, and interpretive policy delineated in NRCS General Manual Title 290, Part 400.11, B(4).

References

  1. 1.Liu, X., Herbert, S. J., Hashemi, A. M., Zhang, X. & Ding, G. Effects of agricultural management on soil organic matter and carbon transformation: A review. Plant Soil Environ. 52, 531–543. https://doi.org/10.17221/3544-pse (2006).CAS  Article  Google Scholar 
  2. 2.Batjes, N. H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 47, 151–163. https://doi.org/10.1111/j.1365-2389.1996.tb01386.x (1996).CAS  Article  Google Scholar 
  3. 3.Wan, Y. F., Lin, E., Xiong, W., Li, Y. & Guo, L. P. Modeling the impact of climate change on soil organic carbon stock in upland soils in the 21st century in China. Agric. Ecosyst. Environ. 141, 23–31. https://doi.org/10.1016/j.agee.2011.02.004 (2011).Article  Google Scholar 
  4. 4.Alvaro-Fuentes, J., Easter, M. & Paustian, K. Climate change effects on organic carbon storage in agricultural soils of northeastern Spain. Agric. Ecosyst. Environ. 155, 87–94. https://doi.org/10.1016/j.agee.2012.04.001 (2012).CAS  Article  Google Scholar 
  5. 5.Wiesmeier, M. et al. Projected loss of soil organic carbon in temperate agricultural soils in the 21st century: Effects of climate change and carbon input trends. Sci. Rep. https://doi.org/10.1038/srep32525 (2016).Article  PubMed  PubMed Central  Google Scholar 
  6. 6.Adams, R. M. et al. Global climate change and US agriculture. Nature 345, 219–224 (1990).ADS  Article  Google Scholar 
  7. 7.Challinor, A. J. et al. A meta-analysis of crop yield under climate change and adaptation. Nat. Clim. Change 4, 287–291. https://doi.org/10.1038/nclimate2153 (2014).ADS  Article  Google Scholar 
  8. 8.IPCC. Fifth Assessment Report. (Cambridge University Press, Cambridge, 2014).
  9. 9.Liu, L. & Basso, B. Impacts of climate variability and adaptation strategies on crop yields and soil organic carbon in the US Midwest. PLoS ONE 15, e0225433. https://doi.org/10.1371/journal.pone.0225433 (2020).CAS  Article  PubMed  PubMed Central  Google Scholar 
  10. 10.Parton, W. J., Schimel, D. S., Cole, C. V. & Ojima, D. S. Analysis of factors controlling soil organic matter levels in Great Plains grasslands. Soil Sci. Soc. Am. J. 51, 1173–1179 (1987).ADS  CAS  Article  Google Scholar 
  11. 11.Knorr, W., Prentice, I. C., House, J. I. & Holland, E. A. Long-term sensitivity of soil carbon turnover to warming. Nature 433, 298–301. https://doi.org/10.1038/nature03226 (2005).ADS  CAS  Article  PubMed  Google Scholar 
  12. 12.Ogle, S. M., Breidt, F. J. & Paustian, K. Agricultural management impacts on soil organic carbon storage under moist and dry climatic conditions of temperate and tropical regions. Biogeochemistry 72, 87–121. https://doi.org/10.1007/s10533-004-0360-2 (2005).Article  Google Scholar 
  13. 13.NASS. USDA-National Agricultural Statistics Service (NASS) Quick stats. https://quickstats.nass.usda.gov/. (2017).
  14. 14.Parton, W. J., Hartman, M., Ojima, D. & Schimel, D. DAYCENT and its land surface submodel: Description and testing. Glob. Planet. Change 19, 35–48. https://doi.org/10.1016/s0921-8181(98)00040-x (1998).ADS  Article  Google Scholar 
  15. 15.Zhang, Y., Suyker, A. & Paustian, K. Improved crop canopy and water balance dynamics for agroecosystem modeling using DayCent. Agron. J. 110, 511–524. https://doi.org/10.2134/agronj2017.06.0328 (2018).CAS  Article  Google Scholar 
  16. 16.Rogers, H. H., Bingham, G. E., Cure, J. D., Smith, J. M. & Surano, K. A. Responses of selected plant species to elevated carbon dioxide in the field. J. Environ. Qual. 12, 569–574. https://doi.org/10.2134/jeq1983.00472425001200040028x (1983).CAS  Article  Google Scholar 
  17. 17.Rose, G., Osborne, T., Greatrex, H. & Wheeler, T. Impact of progressive global warming on the global-scale yield of maize and soybean. Clim. Change 134, 417–428. https://doi.org/10.1007/s10584-016-1601-9 (2016).ADS  Article  Google Scholar 
  18. 18.Kent, J. & Paustian, K. In AGU Fall Meeting Abstracts.
  19. 19.Deryng, D., Sacks, W. J., Barford, C. C. & Ramankutty, N. Simulating the effects of climate and agricultural management practices on global crop yield. Glob. Biogeochem. Cycles. https://doi.org/10.1029/2009gb003765 (2011).Article  Google Scholar 
  20. 20.Osborne, T., Rose, G. & Wheeler, T. Variation in the global-scale impacts of climate change on crop productivity due to climate model uncertainty and adaptation. Agric. For. Meteorol. 170, 183–194. https://doi.org/10.1016/j.agrformet.2012.07.006 (2013).ADS  Article  Google Scholar 
  21. 21.Prince, S. D., Haskett, J., Steininger, M., Strand, H. & Wright, R. Net primary production of U.S. Midwest croplands from agricultural harvest yield data. Ecol. Appl. 11, 1194–1205. https://doi.org/10.1890/1051-0761(2001)011[1194:nppous]2.0.co;2 (2001).Article  Google Scholar 
  22. 22.CTIC. Crop Residue Management Survey. (Conservation Technology Information Center, West Lafayette, Indiana, USA, 2008).
  23. 23.Kuo, S., Sainju, U. M. & Jellum, E. J. Winter cover crop effects on soil organic carbon and carbohydrate in soil. Soil Sci. Soc. Am. J. 61, 145–152. https://doi.org/10.2136/sssaj1997.03615995006100010022x (1997).ADS  CAS  Article  Google Scholar 
  24. 24.Wilhelm, W. W., Johnson, J. M. E., Karlen, D. L. & Lightle, D. T. Corn stover to sustain soil organic carbon further constrains Biomass supply. Agron. J. 99, 1665–1667. https://doi.org/10.2134/agronj2007.0150 (2007).CAS  Article  Google Scholar 
  25. 25.Pries, C. E. H., Castanha, C., Porras, R. C. & Torn, M. S. The whole-soil carbon flux in response to warming. Science 355, 1420–1422. https://doi.org/10.1126/science.aal1319 (2017).ADS  CAS  Article  Google Scholar 
  26. 26.Zimnicki, T. et al. On quantifying water quality benefits of healthy soils. BioScience 70, 343–352 (2020).Article  Google Scholar 
  27. 27.Sanderman, J., Hengl, T. & Fiske, G. J. Soil carbon debt of 12,000 years of human land use. Proc. Natl. Acad. Sci. U.S.A. 114, 9575–9580. https://doi.org/10.1073/pnas.1706103114 (2017).CAS  Article  PubMed  PubMed Central  Google Scholar 
  28. 28.Paustian, K. et al. Climate-smart soils. Nature 532, 49–57. https://doi.org/10.1038/nature17174 (2016).ADS  CAS  Article  PubMed  Google Scholar 
  29. 29.Ruane, A. C., Goldberg, R. & Chryssanthacopoulos, J. Climate forcing datasets for agricultural modeling: Merged products for gap-filling and historical climate series estimation. Agric. For. Meteorol. 200, 233–248. https://doi.org/10.1016/j.agrformet.2014.09.016 (2015).ADS  Article  Google Scholar 
  30. 30.Fowler, H. J., Blenkinsop, S. & Tebaldi, C. Linking climate change modelling to impacts studies: Recent advances in downscaling techniques for hydrological modelling. Int. J. Climatol. 27, 1547–1578. https://doi.org/10.1002/joc.1556 (2007).Article  Google Scholar 
  31. 31.Yukimoto, S. et al. Meteorological Research Institute-Earth System Model Version 1 (MRI-ESM1)—model description. Technical Report of the Meteorological Research Institute 64. (2011).
  32. 32.Watanabe, S. et al. MIROC-ESM 2010: Model description and basic results of CMIP5-20c3m experiments. Geosci. Model Dev. 4, 845–872. https://doi.org/10.5194/gmd-4-845-2011 (2011).ADS  Article  Google Scholar 
  33. 33.Griffies, S. M. et al. The GFDL CM3 Coupled Climate Model: Characteristics of the ocean and sea ice simulations. J. Clim. 24, 3520–3544. https://doi.org/10.1175/2011jcli3964.1 (2011).ADS  Article  Google Scholar 
  34. 34.NRCS. Soil Survey Geographic (SSURGO) Database. USDA Natural Resources Conservation Service. (2009).
  35. 35.NASS. Field Crops: Usual Planting and Harvesting Dates. (USDA, National Agricultural Statistics Service, 2010).
  36. 36.Basso, B. et al. Procedures for initializing soil organic carbon pools in the DSSAT-CENTURY model for agricultural systems. Soil Sci. Soc. Am. J. 75, 69–78. https://doi.org/10.2136/sssaj2010.0115 (2011).ADS  CAS  Article  Google Scholar 
  37. 37.USEPA. Inventory of U.S. Greenhouse Gas Emissions and Sinks: 1990–2014. (2015).
  38. 38.Zhang, Y. et al. DayCent model predictions of NPP and grain yields for agricultural lands in the contiguous U.S. J. Geophys. Res. Biogeosci. 125, e2020JG005750. https://doi.org/10.1029/2020JG005750 (2020).Article  Google Scholar 
  39. 39.Kahle, D. & Wickham, H. ggmap: Spatial visualization with ggplot2. R Journal 5, 144–161 (2013).Article  Google Scholar 

Download references

Acknowledgements

This study was supported by Grants from USDA-NIFA (2012-67003-19904). The authors thank Paul Adler for critical discussions.

Author information

Affiliations

  1. Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO, USAYao Zhang, Ernie Marx, Stephen Williams, Ram Gurung, Stephen Ogle & Keith Paustian
  2. Department of Ecosystem Science and Sustainability, Colorado State University, Fort Collins, CO, USAStephen Ogle
  3. Center for Climate Systems Research, Columbia University, New York, NY, USARadley Horton & Daniel Bader
  4. Department of Soil and Crop Sciences, Colorado State University, Fort Collins, CO, USAKeith Paustian
]]>

Join our effort to build a greener future.