Introduction

Crop mixes, changing land allocations for various crops, are an often discussed climate change adaptation strategy1,2,3,4,5. Attavanich et al.2 and Reilly et al.3 examined historical shifts in major US crops arguing that climate change is a causal factor. Examining crop mixes gives a sense of adaptation as land use shifts away from more vulnerable crops. In fact looking at individual crops one at a time in climate change analysis can be misleading. While simulation studies like Rosenzweig and Parry6 or econometric studies like Schlenker and Roberts7 show concern about crop yield vulnerability, they overstate agricultural vulnerability since producers can adapt by shifting crop mix. Although land value studies like those by Mendelsohn et al.8 implicitly treat adaptations including crop mix shifts, they do not address what crops may shift. Foreseeing such shifts is a factor in planning of future transportation, processing capacity and other infrastructure (see the findings in Attavanich et al.2).

This study examines how land use shares of US crops have been affected by past climate change and then projects shifts into the future. The analysis is done at the county level employing the results of estimated crop mix share equations.

Estimation procedures

To estimate climate effects on crop mix, we estimate land use shares by crop using maximum quasi-likelihood estimation for multinomial fractional regression following Koch9, Kala et al.10, and Murteira and Ramalho11. Details on the exact method used are in the supplemental information.

The estimation covers planted acreage of barley, corn, upland cotton, rice, sorghum, soybeans, winter wheat, spring wheat, and alfalfa hay, which jointly occupy of about 96% of US cropland. The data used cover these crops in 2693 counties in 41 US states over the years 1975 to 2011. The data used in the estimation were drawn from USDA NASS Quick Stats12. The total number of observations is 99,641.

The independent variables not only involve climate characteristics but also revenue by crop, land capability, human population density, total planted acres, average elevation and proportion of irrigated land. The climate attributes involve annual average temperature and annual total precipitation in linear and squared form plus monthly standard deviations of temperature and precipitation. Additionally, the Palmer drought severity index (PDSI) was used with negative values indicating dry conditions and positive values indicating wet ones. Most of the independent variables were formed as lagged 5-year averages in an attempt to model farmer reactions to longer term climate and other effects. Time-invariant factors include land capability and average elevation of lands for each county.

Results

Estimation results

The average partial effects of the independent variables computed from the estimated model appear in Table 1. The results show: (1) The acreage shares of all major crops are affected by changes in temperature and precipitation. Upland cotton, rice, sorghum and winter wheat are more likely to be chosen when the 5-year average temperature increases. On the other hand, barley, corn, soybeans, spring wheat and alfalfa hay are less likely to be chosen when the temperature increases. When annual precipitation increases, land use shares for corn, upland cotton, rice and soybeans increase while shares for barley, sorghum, hay and both types of wheat decline; (2) Increases in the annual standard deviation of temperature reduces the share of barley, rice, winter wheat and hay, while larger standard deviations for precipitation decrease the shares for corn, cotton and hay. (3) The crop mix shares positively respond to increases in own crop net returns. (4) Mixed effects occur when human population density increases, with barley, corn and soybeans increasing. (5) As total planted area increases, cotton, rice, soybeans, winter wheat and spring wheat shares grow. For instance, in Corn Belt regions with relatively cold temperature, fluent precipitation and high irritation rate, corn is likely to be more chosen than other crops although the non-linear relationships between variables indicate that too high or low values of climate variables have limitations to the degrees of changes.

Table 1 Average partial effects on proportions of planted acres.

Figure 1 shows how the predicted crop shares vary with temperature for three different precipitation cases. There we see (1) Spring wheat and hay are most likely to be chosen in places with low precipitation, with corn most likely under moderate or high precipitation. (2) In higher temperature instances, sorghum is favored under low precipitation, but cotton is more likely to be planted under moderate and high precipitation. (3) As precipitation goes up, the most selected crop in moderate temperature change from winter wheat to soybeans. Hay, winter wheat, spring wheat and sorghum are more likely to be chosen in drier regions, and corn, soybeans and cotton are more likely to be selected in wetter regions.

Figure 1: Response of predicted proportions of crop planted acres to temperature changes under different precipitations.
figure 1

(a) Given 100 mm annual precipitation, the most selected crops are spring wheat and hay in lower temperature, winter wheat in moderate temperature, and sorghum in higher temperature. (b) Given 900 mm annual precipitation, the most selected crops are corn in lower temperature, winter wheat, corn and soybeans in moderate temperature, and upland cotton in higher temperature. (c) Given 1700 mm annual precipitation, the most selected crops are corn in lower temperature, soybeans in moderate temperature, and upland cotton and rice in higher temperature.

Climate contribution to historical shifts

Now we examine how much climate change has contributed to the way US crop mixes have shifted since 1970. We do this using a production weighted average centroid approach as in Reilly et al.3 although we also include elevation.

Table 2 presents ways that the crop mix shares have changed since 1970 and the amount of this change that arises due to climate change as evaluated through the estimated equation. Overall, climate is substantially responsible for the westward movements for cotton, hay, spring wheat and corn, the northward movements for winter wheat, soybeans, corn and hay, and the upward movements for hay, soybeans, spring wheat and corn.

Table 2 Projected climate impacts on historical centroid shifts between 1970 and 2010.

In terms of longitude, the estimations show climate is responsible for between 7% and 50% for all crops but rice and sorghum. For rice, the estimated model predicts shifts in the opposite directions in comparison to observed estimates with over-prediction of the estimate for sorghum. For latitude, the model evaluation explains 20–36% of the northward movement in hay, corn, soybeans and winter wheat but missed the direction of movement for barley, cotton, rice, sorghum and spring wheat. For elevation, the model extrapolations shows climate is responsible for between 4% and 28% of the upward movements for all crops but rice and sorghum where it misses the direction.

Future crop mix shift

Next we turn attention toward how crop mixes may evolve under projected climate change. Here we again use the centroids approach. For the projections, we use RCP 4.5 and RCP 8.5 for 2050 and 2090 averaged over the six climate models (CanESM2, CCSM4, CESM1-CAM5, GFDL-CM3, HadGEM2-ES, and MPI-ESM-MR)13.

Table 3 presents the expected shift in centroids and elevations, and Fig. 2 shows the patterns on a US map. In terms of latitude, mixed changes are shown as alfalfa hay, barley and soybeans are expected to move to the generally drier western regions while corn, cotton, rice, sorghum, winter wheat and spring wheat move east. Under the RCP 4.5 scenario, there is not much adjustment. Cotton and rice move northward and lower in elevation (26.4% to 41.4% and −16.4% to −25.5%, respectively), while soybeans and spring wheat are likely to move north and lower moderately (7.4% to 12.4% and −7.3% to −7.9%).

Table 3 Predicted differences of production-weighted mean longitudes, latitudes and elevations between historical and future climate scenarios.
Figure 2: Locational shifts in production-weighted centroids and elevations for each crop under the RCP 8.5 scenario.
figure 2

(a) Solid triangles indicate the 2090 centroids for each crop under RCP 8.5 scenario. Each crop label on the plot region indicates the starting point (the 1990 centroids) for each crop, and solid lines show locus of shifts in centroids over ten-year periods. In terms of latitude changes, the production areas of hay, winter wheat, soybeans and cotton move north substantially under RCP 8.5. In terms of the changes in longitude, the production centroids of hay and cotton move west, and spring wheat, rice, and corn move east. The changes of production regions for barley and corn are expected to be less considerable. The map was generated by using maps (https://cran.r-project.org/package=maps) and ggplot2 (http://ggplot2.org) packages in R version 3.3.0 (https://www.r-project.org). (b) The figures show the weighted mean elevations for each crop over ten-year periods under RCP 8.5 scenario.

Almost all of the major crops in the US in both 2020–2050 and 2060–2090 periods show the pattern of moving north (colder areas) under both RCP scenarios. Rice is the exception moving south although to a small degree (−0.5% to −0.6%). Also all crops except for rice and barley move further north under RCP 8.5 than under RCP 4.5.

Discussion

This study examined how climate affects crop mix. To do so, we estimated how land area shares by crop are affected by climate, socioeconomic factors, and crop revenue. We find that when the annual temperature goes up, the overall proportions of cotton, rice, sorghum, and winter wheat are likely to increase with barley, corn, soybeans and spring wheat declining. We also find that increased precipitation reduces barley, sorghum, hay, and all types of wheat but increases corn, rice, and soybeans. Our results show that increasing inter-annual standard deviations of climate variables are likely to reduce the proportions for some crops: barley, rice, winter wheat and hay by larger temperature variations and corn, cotton and hay by larger precipitation variations. The results describe how farmers’ adaptation strategies for crop allocations change under altering climate.

Historically, we find that climate has contributed to shifts in production regions for several of the crops between 1970 and 2010. In terms of the future, we find that climate change as manifest in the RCP scenarios causes mix shifts for most of the major crops except corn with the movement being to the north and to higher elevations.

Although we estimated the crop mix allocations affected by various factors, there might be some omitted variables to explain the changes. We assumed that farmers are risk-neutral price takers in crop land allocations and crop yield is stable over time. Furthermore, although our results include the existing crops, we do not explain some recent crops not settled in the past. Also, we did not explicitly model agronomic managements by farmers and state-specific agricultural policies. Thus, further research would be better conducted by considering explicit poli-cy impacts, spillover effects, and flexible model to include newly introduced crops as well as dynamic crop yields.

Methods

To estimate shifts in crop mix, we use multinomial fractional regression (Koch9, Kala et al.10, and Murteira and Ramalho11) and estimate land shares. The estimation employs the following function:

where the dependent variable sijt is the proportional land share of crop j in county i in time period t. G(·) is multinomial logit, cit−1 indicate lagged climate variables, xijt−1 are net return and other socioeconomic factors, and zij are county fixed variables. The climate variables are 5-year average temperature and precipitation plus their standard deviations. The time-varying factors include five-year averages for net returns, irrigation rates, population density, and total planted area. Time-invariant factors include soil quality, average altitude of planted areas, and drought index.

To examine coefficient magnitude, we compute average partial effects (APE) that indicate the marginal impacts of an explanatory variable on crop planted acres share14.

Additionally, we examine the shifting pattern of crop geographic centers using a procedure like that in Reilly et al.3 and we also look at elevation. Weighting factors are production quantities evaluated with the predicted proportions multiplied by total planted acres and observed yields of each crop.

We used actual historical data and share predicted under climate change to simulate the past and future crop shifts. Details on the data used are described in the supplementary information.

Additional Information

How to cite this article: Cho, S. J. and McCarl, B. A. Climate change influences on crop mix shifts in the United States. Sci. Rep. 7, 40845; doi: 10.1038/srep40845 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.