skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Understanding the non-stationary relationships between corn yields and meteorology via a spatiotemporally varying coefficient model.
The relationships between crop yields and meteorology are naturally non-stationary because of spatiotemporal heterogeneity. Many studies have examined spatial heterogeneity in the regression model, but only limited research has attempted to account for both spatial autocorrelation and temporal variation. In this article, we develop a novel spatiotemporally varying coefficient (STVC) model to understand non-stationary relationships between crop yields and meteorological variables. We compare the proposed model with variant models specialized for time or spatial, namely spatial varying coefficient (SVC) model and temporal varying coefficient (TVC) model. This study was conducted using the county-level corn yield and meteorological data, including seasonal Growing Degree Days (GDD), Killing Degree Days (KDD), Vapor Pressure Deficit (VPD), and precipitation (PCPN), from 1981 to 2018 in three Corn Belt states, including Illinois, Indiana, and Iowa. Allowing model coefficients varying in both temporal and spatial dimensions gives the best performance of STVC in simulating the corn yield responses toward various meteorological conditions. The STVC reduced the root-mean-square error to 10.64 Bu/Ac (0.72 Mg/ha) from 15.68 Bu/Ac (1.06 Mg/ha) for TVC and 16.48 Bu/Ac (1.11 Mg/ha) for SVC. Meanwhile, the STVC resulted in a higher R2 of 0.81 compared to 0.56 for SVC and 0.64 for TVC. The STVC showed better performance in handling spatial dependence of corn production, which tends to cluster estimation residuals when counties are close, with the lowest Moran’s I of 0.10. Considering the spatiotemporal non-stationarity, the proposed model significantly improves the power of the meteorological data in explaining the variations of corn yields.  more » « less
Award ID(s):
1922758
PAR ID:
10291051
Author(s) / Creator(s):
Date Published:
Journal Name:
Agricultural and forest meteorology
Volume:
301-302
Issue:
2021
ISSN:
0168-1923
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Accurate phenological information is essential for monitoring crop development, predicting crop yield, and enhancing resilience to cope with climate change. This study employed a curve-change-based dynamic threshold approach on NDVI (Normalized Differential Vegetation Index) time series to detect the planting and harvesting dates for corn and soybean in Kentucky, a typical climatic transition zone, from 2000 to 2018. We compared satellite-based estimates with ground observations and performed trend analyses of crop phenological stages over the study period to analyze their relationships with climate change and crop yields. Our results showed that corn and soybean planting dates were delayed by 0.01 and 0.07 days/year, respectively. Corn harvesting dates were also delayed at a rate of 0.67 days/year, while advanced soybean harvesting occurred at a rate of 0.05 days/year. The growing season length has increased considerably at a rate of 0.66 days/year for corn and was shortened by 0.12 days/year for soybean. Sensitivity analysis showed that planting dates were more sensitive to the early season temperature, while harvesting dates were significantly correlated with temperature over the entire growing season. In terms of the changing climatic factors, only the increased summer precipitation was statistically related to the delayed corn harvesting dates in Kentucky. Further analysis showed that the increased corn yield was significantly correlated with the delayed harvesting dates (1.37 Bu/acre per day) and extended growing season length (1.67 Bu/acre per day). Our results suggested that seasonal climate change (e.g., summer precipitation) was the main factor influencing crop phenological trends, particularly corn harvesting in Kentucky over the study period. We also highlighted the critical role of changing crop phenology in constraining crop production, which needs further efforts for optimizing crop management practices. 
    more » « less
  2. Abstract Understanding subfield crop yields and temporal stability is critical to better manage crops. Several algorithms have proposed to study within-field temporal variability but they were mostly limited to few fields. In this study, a large dataset composed of 5520 yield maps from 768 fields provided by farmers was used to investigate the influence of subfield yield distribution skewness on temporal variability. The data are used to test two intuitive algorithms for mapping stability: one based on standard deviation and the second based on pixel ranking and percentiles. The analysis of yield monitor data indicates that yield distribution is asymmetric, and it tends to be negatively skewed (p < 0.05) for all of the four crops analyzed, meaning that low yielding areas are lower in frequency but cover a larger range of low values. The mean yield difference between the pixels classified as high-and-stable and the pixels classified as low-and-stable was 1.04 Mg ha −1 for maize, 0.39 Mg ha −1 for cotton, 0.34 Mg ha −1 for soybean, and 0.59 Mg ha −1 for wheat. The yield of the unstable zones was similar to the pixels classified as low-and-stable by the standard deviation algorithm, whereas the two-way outlier algorithm did not exhibit this bias. Furthermore, the increase in the number years of yield maps available induced a modest but significant increase in the certainty of stability classifications, and the proportion of unstable pixels increased with the precipitation heterogeneity between the years comprising the yield maps. 
    more » « less
  3. Abstract With climate change threatening agricultural productivity and global food demand increasing, it is important to better understand which farm management practices will maximize crop yields in various climatic conditions. To assess the effectiveness of agricultural practices, researchers often turn to randomized field experiments, which are reliable for identifying causal effects but are often limited in scope and therefore lack external validity. Recently, researchers have also leveraged large observational datasets from satellites and other sources, which can lead to conclusions biased by confounding variables or systematic measurement errors. Because experimental and observational datasets have complementary strengths, in this paper we propose a method that uses a combination of experimental and observational data in the same analysis. As a case study, we focus on the causal effect of crop rotation on corn (maize) and soybean yields in the Midwestern United States. We find that, in terms of root mean squared error, our hybrid method performs 13% better than using experimental data alone and 26% better than using the observational data alone in the task of predicting the effect of rotation on corn yield at held-out experimental sites. Further, the causal estimates based on our method suggest that benefits of crop rotations on corn yield are lower in years and locations with high temperatures whereas the benefits of crop rotations on soybean yield are higher in years and locations with high temperatures. In particular, we estimated that the benefit of rotation on corn yields (and soybean yields) was 0.85 t ha−1(0.24 t ha−1) on average for the top quintile of temperatures, 1.03 t ha−1(0.21 t ha−1) on average for the whole dataset, and 1.19 t ha−1(0.16 t ha−1) on average for the bottom quintile of temperatures. This association between temperatures and rotation benefits is consistent with the hypothesis that the benefit of the corn-soybean rotation on soybean yield is largely driven by pest pressure reductions while the benefit of the corn-soybean rotation on corn yields is largely driven by nitrogen availability. 
    more » « less
  4. Abstract Crop N decision support tools are typically based on either empirical relationships that lack mechanistic underpinnings or simulation models that are too complex to use on farms with limited input data. We developed an N mineralization model for corn that lies between these endpoints; it includes a mechanistic model structure reflecting microbial and texture controls on N mineralization but requires just a few simple inputs: soil texture soil C and N concentration and cover crop N content and carbon to nitgrogen ratio (C/N). We evaluated a previous version of the model with an independent dataset to determine the accuracy in predictions of unfertilized corn (Zea maysL.) yield across a wider range of soil texture, cover crop, and growing season precipitation conditions. We tested three assumptions used in the original model: (1) soil C/N is equal to 10, (2) yield does not need to be adjusted for growing season precipitation, and (3) sand content controls humification efficiency (ε). The best new model used measured values for soil C/N, had a summertime precipitation adjustment, and included both sand and clay content as predictors ofε(root mean square error [RMSE] = 1.43 Mg ha−1;r= 0.69). In the new model, clay has a stronger influence than sand onε, corresponding to lower predicted mineralization rates on fine‐textured soils. The new model had a reasonable validation fit (RMSE = 1.71 Mg ha−1;r= 0.56) using an independent dataset. Our results indicate the new model is an improvement over the previous version because it predicts unfertilized corn yield for a wider range of conditions. 
    more » « less
  5. The possible influence of global climate changes on agricultural production is becoming increasingly significant, necessitating greater attention to improving agricultural production in response to temperature rises and precipitation variability. As one of the main winter wheat-producing areas in China, the temporal and spatial distribution characteristics of precipitation, accumulated temperature, and actual yield and climatic yield of winter wheat during the growing period in Shanxi Province were analysed in detail. With the utilisation of daily meteorological data collected from 12 meteorological stations in Shanxi Province in 1964–2018, our study analysed the change in winter wheat yield with climate change using GIS combined with wavelet analysis. The results show the following: (1) Accumulated temperature and precipitation are the two most important limiting factors among the main physical factors that impact yield. Based on the analysis of the ArcGIS geographical detector, the correlation between the actual yield of winter wheat and the precipitation during the growth period was the highest, reaching 0.469, and the meteorological yield and accumulated temperature during this period also reached its peak value of 0.376. (2) The regions with more suitable precipitation and accumulated temperature during the growth period of winter wheat in the study area had relatively high actual winter wheat yields. Overall, the average actual yield of the entire region showed a significant increasing trend over time, with an upward trend of 47.827 kg ha−1 yr−1. (3) The variation coefficient of winter wheat climatic yield was relatively stable in 2008–2018. In particular, there were many years of continuous reduction in winter wheat yields prior to 2006. Thereafter, the impact of climate change on winter wheat yields became smaller. This study expands our understanding of the complex interactions between climate variables and crop yield but also provides practical recommendations for enhancing agricultural practices in this region 
    more » « less