Geostatistical interpolation and aggregation of crop growth model outputs

Research output: Contribution to journalArticleAcademicpeer-review

1 Citation (Scopus)

Abstract

Many crop growth models require daily meteorological data. Consequently, model simulations can be obtained only at a limited number of locations, i.e. at weather stations with long-term records of daily data. To estimate the potential crop production at country level, we present in this study a geostatistical approach for spatial interpolation and aggregation of crop growth model outputs. As case study, we interpolated, simulated and aggregated crop growth model outputs of sorghum and millet in West-Africa. We used crop growth model outputs to calibrate a linear regression model using environmental covariates as predictors. The spatial regression residuals were investigated for spatial correlation. The linear regression model and the spatial correlation of residuals together were used to predict theoretical crop yield at all locations using kriging with external drift. A spatial standard deviation comes along with this prediction, indicating the uncertainty of the prediction. In combination with land use data and country borders, we summed the crop yield predictions to determine an area total. With spatial stochastic simulation, we estimated the uncertainty of that total production potential as well as the spatial cumulative distribution function. We compared our results with the prevailing agro-ecological Climate Zones approach used for spatial aggregation. Linear regression could explain up to 70% of the spatial variation of the yield. In three out of four cases the regression residuals showed spatial correlation. The potential crop production per country according to the Climate Zones approach was in all countries and cases except one within the 95% prediction interval as obtained after yield aggregation. We concluded that the geostatistical approach can estimate a country's crop production, including a quantification of uncertainty. In addition, we stress the importance of the use of geostatistics to create tools for crop modelling scientists to explore relationships between yields and spatial environmental variables and to assist policy makers with tangible results on yield gaps at multiple levels of spatial aggregation.

LanguageEnglish
Pages111-121
JournalEuropean Journal of Agronomy
Volume77
DOIs
Publication statusPublished - 2016

Fingerprint

crop models
growth models
interpolation
crop production
crop
prediction
uncertainty
crop yield
cumulative distribution
environmental models
geostatistics
weather stations
kriging
millets
meteorological data
Western Africa
spatial variation
simulation models
land use
case studies

Keywords

  • Geostatistics
  • Spatial aggregation
  • Spatial prediction
  • Uncertainty
  • Yield gap
  • Yield potential

Cite this

@article{656b7a5e2a69431ea9398d8b4b23bf1a,
title = "Geostatistical interpolation and aggregation of crop growth model outputs",
abstract = "Many crop growth models require daily meteorological data. Consequently, model simulations can be obtained only at a limited number of locations, i.e. at weather stations with long-term records of daily data. To estimate the potential crop production at country level, we present in this study a geostatistical approach for spatial interpolation and aggregation of crop growth model outputs. As case study, we interpolated, simulated and aggregated crop growth model outputs of sorghum and millet in West-Africa. We used crop growth model outputs to calibrate a linear regression model using environmental covariates as predictors. The spatial regression residuals were investigated for spatial correlation. The linear regression model and the spatial correlation of residuals together were used to predict theoretical crop yield at all locations using kriging with external drift. A spatial standard deviation comes along with this prediction, indicating the uncertainty of the prediction. In combination with land use data and country borders, we summed the crop yield predictions to determine an area total. With spatial stochastic simulation, we estimated the uncertainty of that total production potential as well as the spatial cumulative distribution function. We compared our results with the prevailing agro-ecological Climate Zones approach used for spatial aggregation. Linear regression could explain up to 70{\%} of the spatial variation of the yield. In three out of four cases the regression residuals showed spatial correlation. The potential crop production per country according to the Climate Zones approach was in all countries and cases except one within the 95{\%} prediction interval as obtained after yield aggregation. We concluded that the geostatistical approach can estimate a country's crop production, including a quantification of uncertainty. In addition, we stress the importance of the use of geostatistics to create tools for crop modelling scientists to explore relationships between yields and spatial environmental variables and to assist policy makers with tangible results on yield gaps at multiple levels of spatial aggregation.",
keywords = "Geostatistics, Spatial aggregation, Spatial prediction, Uncertainty, Yield gap, Yield potential",
author = "Luc Steinbuch and Brus, {Dick J.} and {van Bussel}, {Lenny G.J.} and Heuvelink, {Gerard B.M.}",
year = "2016",
doi = "10.1016/j.eja.2016.03.007",
language = "English",
volume = "77",
pages = "111--121",
journal = "European Journal of Agronomy",
issn = "1161-0301",
publisher = "Elsevier",

}

TY - JOUR

T1 - Geostatistical interpolation and aggregation of crop growth model outputs

AU - Steinbuch, Luc

AU - Brus, Dick J.

AU - van Bussel, Lenny G.J.

AU - Heuvelink, Gerard B.M.

PY - 2016

Y1 - 2016

N2 - Many crop growth models require daily meteorological data. Consequently, model simulations can be obtained only at a limited number of locations, i.e. at weather stations with long-term records of daily data. To estimate the potential crop production at country level, we present in this study a geostatistical approach for spatial interpolation and aggregation of crop growth model outputs. As case study, we interpolated, simulated and aggregated crop growth model outputs of sorghum and millet in West-Africa. We used crop growth model outputs to calibrate a linear regression model using environmental covariates as predictors. The spatial regression residuals were investigated for spatial correlation. The linear regression model and the spatial correlation of residuals together were used to predict theoretical crop yield at all locations using kriging with external drift. A spatial standard deviation comes along with this prediction, indicating the uncertainty of the prediction. In combination with land use data and country borders, we summed the crop yield predictions to determine an area total. With spatial stochastic simulation, we estimated the uncertainty of that total production potential as well as the spatial cumulative distribution function. We compared our results with the prevailing agro-ecological Climate Zones approach used for spatial aggregation. Linear regression could explain up to 70% of the spatial variation of the yield. In three out of four cases the regression residuals showed spatial correlation. The potential crop production per country according to the Climate Zones approach was in all countries and cases except one within the 95% prediction interval as obtained after yield aggregation. We concluded that the geostatistical approach can estimate a country's crop production, including a quantification of uncertainty. In addition, we stress the importance of the use of geostatistics to create tools for crop modelling scientists to explore relationships between yields and spatial environmental variables and to assist policy makers with tangible results on yield gaps at multiple levels of spatial aggregation.

AB - Many crop growth models require daily meteorological data. Consequently, model simulations can be obtained only at a limited number of locations, i.e. at weather stations with long-term records of daily data. To estimate the potential crop production at country level, we present in this study a geostatistical approach for spatial interpolation and aggregation of crop growth model outputs. As case study, we interpolated, simulated and aggregated crop growth model outputs of sorghum and millet in West-Africa. We used crop growth model outputs to calibrate a linear regression model using environmental covariates as predictors. The spatial regression residuals were investigated for spatial correlation. The linear regression model and the spatial correlation of residuals together were used to predict theoretical crop yield at all locations using kriging with external drift. A spatial standard deviation comes along with this prediction, indicating the uncertainty of the prediction. In combination with land use data and country borders, we summed the crop yield predictions to determine an area total. With spatial stochastic simulation, we estimated the uncertainty of that total production potential as well as the spatial cumulative distribution function. We compared our results with the prevailing agro-ecological Climate Zones approach used for spatial aggregation. Linear regression could explain up to 70% of the spatial variation of the yield. In three out of four cases the regression residuals showed spatial correlation. The potential crop production per country according to the Climate Zones approach was in all countries and cases except one within the 95% prediction interval as obtained after yield aggregation. We concluded that the geostatistical approach can estimate a country's crop production, including a quantification of uncertainty. In addition, we stress the importance of the use of geostatistics to create tools for crop modelling scientists to explore relationships between yields and spatial environmental variables and to assist policy makers with tangible results on yield gaps at multiple levels of spatial aggregation.

KW - Geostatistics

KW - Spatial aggregation

KW - Spatial prediction

KW - Uncertainty

KW - Yield gap

KW - Yield potential

U2 - 10.1016/j.eja.2016.03.007

DO - 10.1016/j.eja.2016.03.007

M3 - Article

VL - 77

SP - 111

EP - 121

JO - European Journal of Agronomy

T2 - European Journal of Agronomy

JF - European Journal of Agronomy

SN - 1161-0301

ER -