Geostatistical disaggregation of polygon maps of average crop yields by area-to-point kriging

D.J. Brus, H. Boogaard, T. Ceccarelli, T.G. Orton, S. Traore, M. Zhang

Research output: Contribution to journalArticleAcademicpeer-review

1 Citation (Scopus)

Abstract

Crop yield data are often available as statistics of areas, such as administrative units, generated by national agricultural surveys and censuses. This paper shows that such areal data can be used in area-to-point kriging (ATP kriging) to estimate the crop yield at the nodes of a fine grid that discretizes the study area, so that a more detailed map of the crop yield is obtained. The theory behind ATP kriging is explained, and illustrated with a one-dimensional simulation study and two real-world case studies. Vegetation, precipitation, temperature and soil data were used as potential covariates in the spatial trend part of the geostatistical model. ATP kriging requires the covariogram at point support, which can be recovered from the areal data by restricted maximum likelihood. The standard errors of the estimated variogram parameters can then be obtained by the Fisher information matrix. The average yields of only 17 administrative units in Shandong province (China) were not enough to obtain reliable estimates of the covariogram at point support. Also the ranges of the regional averages of the covariates were very narrow, so that the model must be extrapolated in the largest part of the study area. We were more confident about the covariogram parameters estimated from 45 provinces in Burkina Faso. We conclude that ATP kriging is an interesting method for disaggregation of spatially averaged crop yields. Contrary to other downscaling methods ATP kriging is founded on statistical theory, and consequently provides estimates of the precision of the disaggregated yields. Shortcomings are related to the uncertainty in the estimated covariogram parameters, as well as to the extrapolation of the model outside the range of the regional means of the covariates. Opportunities for future advancements are the use of modelled yields as covariates and the introduction of expert knowledge at different levels. For the latter a Bayesian approach to ATP kriging can be advantageous, introducing prior knowledge about the model parameters, as well as accounting for uncertainty about the model parameters.

LanguageEnglish
Pages48-59
JournalEuropean Journal of Agronomy
Volume97
DOIs
Publication statusPublished - 1 Jul 2018

Fingerprint

kriging
polygon
crop yield
model uncertainty
Burkina Faso
China
expert opinion
downscaling
variogram
census
uncertainty
statistics
parameter
case studies
vegetation
matrix
methodology
simulation
soil
temperature

Keywords

  • Aggregated data
  • Uncertainty
  • Yield gap

Cite this

@article{3ea6dd6593894f7284fae2295255389a,
title = "Geostatistical disaggregation of polygon maps of average crop yields by area-to-point kriging",
abstract = "Crop yield data are often available as statistics of areas, such as administrative units, generated by national agricultural surveys and censuses. This paper shows that such areal data can be used in area-to-point kriging (ATP kriging) to estimate the crop yield at the nodes of a fine grid that discretizes the study area, so that a more detailed map of the crop yield is obtained. The theory behind ATP kriging is explained, and illustrated with a one-dimensional simulation study and two real-world case studies. Vegetation, precipitation, temperature and soil data were used as potential covariates in the spatial trend part of the geostatistical model. ATP kriging requires the covariogram at point support, which can be recovered from the areal data by restricted maximum likelihood. The standard errors of the estimated variogram parameters can then be obtained by the Fisher information matrix. The average yields of only 17 administrative units in Shandong province (China) were not enough to obtain reliable estimates of the covariogram at point support. Also the ranges of the regional averages of the covariates were very narrow, so that the model must be extrapolated in the largest part of the study area. We were more confident about the covariogram parameters estimated from 45 provinces in Burkina Faso. We conclude that ATP kriging is an interesting method for disaggregation of spatially averaged crop yields. Contrary to other downscaling methods ATP kriging is founded on statistical theory, and consequently provides estimates of the precision of the disaggregated yields. Shortcomings are related to the uncertainty in the estimated covariogram parameters, as well as to the extrapolation of the model outside the range of the regional means of the covariates. Opportunities for future advancements are the use of modelled yields as covariates and the introduction of expert knowledge at different levels. For the latter a Bayesian approach to ATP kriging can be advantageous, introducing prior knowledge about the model parameters, as well as accounting for uncertainty about the model parameters.",
keywords = "Aggregated data, Uncertainty, Yield gap",
author = "D.J. Brus and H. Boogaard and T. Ceccarelli and T.G. Orton and S. Traore and M. Zhang",
year = "2018",
month = "7",
day = "1",
doi = "10.1016/j.eja.2018.05.003",
language = "English",
volume = "97",
pages = "48--59",
journal = "European Journal of Agronomy",
issn = "1161-0301",
publisher = "Elsevier",

}

Geostatistical disaggregation of polygon maps of average crop yields by area-to-point kriging. / Brus, D.J.; Boogaard, H.; Ceccarelli, T.; Orton, T.G.; Traore, S.; Zhang, M.

In: European Journal of Agronomy, Vol. 97, 01.07.2018, p. 48-59.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Geostatistical disaggregation of polygon maps of average crop yields by area-to-point kriging

AU - Brus, D.J.

AU - Boogaard, H.

AU - Ceccarelli, T.

AU - Orton, T.G.

AU - Traore, S.

AU - Zhang, M.

PY - 2018/7/1

Y1 - 2018/7/1

N2 - Crop yield data are often available as statistics of areas, such as administrative units, generated by national agricultural surveys and censuses. This paper shows that such areal data can be used in area-to-point kriging (ATP kriging) to estimate the crop yield at the nodes of a fine grid that discretizes the study area, so that a more detailed map of the crop yield is obtained. The theory behind ATP kriging is explained, and illustrated with a one-dimensional simulation study and two real-world case studies. Vegetation, precipitation, temperature and soil data were used as potential covariates in the spatial trend part of the geostatistical model. ATP kriging requires the covariogram at point support, which can be recovered from the areal data by restricted maximum likelihood. The standard errors of the estimated variogram parameters can then be obtained by the Fisher information matrix. The average yields of only 17 administrative units in Shandong province (China) were not enough to obtain reliable estimates of the covariogram at point support. Also the ranges of the regional averages of the covariates were very narrow, so that the model must be extrapolated in the largest part of the study area. We were more confident about the covariogram parameters estimated from 45 provinces in Burkina Faso. We conclude that ATP kriging is an interesting method for disaggregation of spatially averaged crop yields. Contrary to other downscaling methods ATP kriging is founded on statistical theory, and consequently provides estimates of the precision of the disaggregated yields. Shortcomings are related to the uncertainty in the estimated covariogram parameters, as well as to the extrapolation of the model outside the range of the regional means of the covariates. Opportunities for future advancements are the use of modelled yields as covariates and the introduction of expert knowledge at different levels. For the latter a Bayesian approach to ATP kriging can be advantageous, introducing prior knowledge about the model parameters, as well as accounting for uncertainty about the model parameters.

AB - Crop yield data are often available as statistics of areas, such as administrative units, generated by national agricultural surveys and censuses. This paper shows that such areal data can be used in area-to-point kriging (ATP kriging) to estimate the crop yield at the nodes of a fine grid that discretizes the study area, so that a more detailed map of the crop yield is obtained. The theory behind ATP kriging is explained, and illustrated with a one-dimensional simulation study and two real-world case studies. Vegetation, precipitation, temperature and soil data were used as potential covariates in the spatial trend part of the geostatistical model. ATP kriging requires the covariogram at point support, which can be recovered from the areal data by restricted maximum likelihood. The standard errors of the estimated variogram parameters can then be obtained by the Fisher information matrix. The average yields of only 17 administrative units in Shandong province (China) were not enough to obtain reliable estimates of the covariogram at point support. Also the ranges of the regional averages of the covariates were very narrow, so that the model must be extrapolated in the largest part of the study area. We were more confident about the covariogram parameters estimated from 45 provinces in Burkina Faso. We conclude that ATP kriging is an interesting method for disaggregation of spatially averaged crop yields. Contrary to other downscaling methods ATP kriging is founded on statistical theory, and consequently provides estimates of the precision of the disaggregated yields. Shortcomings are related to the uncertainty in the estimated covariogram parameters, as well as to the extrapolation of the model outside the range of the regional means of the covariates. Opportunities for future advancements are the use of modelled yields as covariates and the introduction of expert knowledge at different levels. For the latter a Bayesian approach to ATP kriging can be advantageous, introducing prior knowledge about the model parameters, as well as accounting for uncertainty about the model parameters.

KW - Aggregated data

KW - Uncertainty

KW - Yield gap

U2 - 10.1016/j.eja.2018.05.003

DO - 10.1016/j.eja.2018.05.003

M3 - Article

VL - 97

SP - 48

EP - 59

JO - European Journal of Agronomy

T2 - European Journal of Agronomy

JF - European Journal of Agronomy

SN - 1161-0301

ER -