Design-based Generalized Least Squares estimation of status and trend of soil properties from monitoring data

D.J. Brus, J.J. de Gruijter

Research output: Contribution to journalArticleAcademicpeer-review

15 Citations (Scopus)

Abstract

This paper introduces and demonstrates design-based Generalized Least Squares (GLS) estimation of spatial means at selected time points from data collected in repeated soil surveys with partial overlap, such as a rotating and a supplemented panel. The linear time trend of the spatial means can then be obtained as a linear combination of the estimated spatial means. The GLS estimator is the minimum variance linear unbiased estimator. Five space–time designs were compared under a first order autoregressive time series model for the spatial means, through the average and generalized sampling variances of the estimated spatial means, and through the sampling variance of the estimated linear trend. None of the designs scored best on all three quality measures. If the aim of soil monitoring is estimation of both status and trend, then these two conflicting aims must be prioritized in order to choose an efficient space–time design. The methodology is demonstrated by a case study on eutrophication and acidification of forest soils. The linear trends in the spatial means of pH and the ammonium and nitrate concentrations at three depths in the soil profile, as estimated from a rotational design with four sampling times at an interval of one year, were small and not significant. Exceptions were pH in the subsoil (- 0.06 pH units yr- 1) and ammonium in the middle soil horizon (- 0.086 mg N kg- 1 yr- 1). The linear trend is here defined as a linear combination of the true and unknown, but fixed spatial means. In quantifying the uncertainty of the estimated trend, only the sampling error in the estimated spatial means is accounted for. If there is a need to include uncertainty due to fluctuations of the true spatial means around a linear trend, then a super-population or time series model for the spatial means must be postulated which comprises a model error term. The linear trend is then defined as a model parameter, that can be estimated by Generalized Least Squares as in generalized linear modeling. --------------------------------------------------------------------------------
LanguageEnglish
Pages172-180
JournalGeoderma
Volume164
Issue number3-4
DOIs
Publication statusPublished - 2011

Fingerprint

least squares
soil properties
soil property
monitoring
time series analysis
uncertainty
sampling
soil surveys
soil horizons
subsoil
ammonium
forest soils
acidification
eutrophication
soil profiles
time series
nitrates
trend
monitoring data
soil survey

Cite this

@article{3d486b0d959c442580b0c5d3f0d0f42c,
title = "Design-based Generalized Least Squares estimation of status and trend of soil properties from monitoring data",
abstract = "This paper introduces and demonstrates design-based Generalized Least Squares (GLS) estimation of spatial means at selected time points from data collected in repeated soil surveys with partial overlap, such as a rotating and a supplemented panel. The linear time trend of the spatial means can then be obtained as a linear combination of the estimated spatial means. The GLS estimator is the minimum variance linear unbiased estimator. Five space–time designs were compared under a first order autoregressive time series model for the spatial means, through the average and generalized sampling variances of the estimated spatial means, and through the sampling variance of the estimated linear trend. None of the designs scored best on all three quality measures. If the aim of soil monitoring is estimation of both status and trend, then these two conflicting aims must be prioritized in order to choose an efficient space–time design. The methodology is demonstrated by a case study on eutrophication and acidification of forest soils. The linear trends in the spatial means of pH and the ammonium and nitrate concentrations at three depths in the soil profile, as estimated from a rotational design with four sampling times at an interval of one year, were small and not significant. Exceptions were pH in the subsoil (- 0.06 pH units yr- 1) and ammonium in the middle soil horizon (- 0.086 mg N kg- 1 yr- 1). The linear trend is here defined as a linear combination of the true and unknown, but fixed spatial means. In quantifying the uncertainty of the estimated trend, only the sampling error in the estimated spatial means is accounted for. If there is a need to include uncertainty due to fluctuations of the true spatial means around a linear trend, then a super-population or time series model for the spatial means must be postulated which comprises a model error term. The linear trend is then defined as a model parameter, that can be estimated by Generalized Least Squares as in generalized linear modeling. --------------------------------------------------------------------------------",
author = "D.J. Brus and {de Gruijter}, J.J.",
year = "2011",
doi = "10.1016/j.geoderma.2011.06.001",
language = "English",
volume = "164",
pages = "172--180",
journal = "Geoderma",
issn = "0016-7061",
publisher = "Elsevier",
number = "3-4",

}

Design-based Generalized Least Squares estimation of status and trend of soil properties from monitoring data. / Brus, D.J.; de Gruijter, J.J.

In: Geoderma, Vol. 164, No. 3-4, 2011, p. 172-180.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Design-based Generalized Least Squares estimation of status and trend of soil properties from monitoring data

AU - Brus, D.J.

AU - de Gruijter, J.J.

PY - 2011

Y1 - 2011

N2 - This paper introduces and demonstrates design-based Generalized Least Squares (GLS) estimation of spatial means at selected time points from data collected in repeated soil surveys with partial overlap, such as a rotating and a supplemented panel. The linear time trend of the spatial means can then be obtained as a linear combination of the estimated spatial means. The GLS estimator is the minimum variance linear unbiased estimator. Five space–time designs were compared under a first order autoregressive time series model for the spatial means, through the average and generalized sampling variances of the estimated spatial means, and through the sampling variance of the estimated linear trend. None of the designs scored best on all three quality measures. If the aim of soil monitoring is estimation of both status and trend, then these two conflicting aims must be prioritized in order to choose an efficient space–time design. The methodology is demonstrated by a case study on eutrophication and acidification of forest soils. The linear trends in the spatial means of pH and the ammonium and nitrate concentrations at three depths in the soil profile, as estimated from a rotational design with four sampling times at an interval of one year, were small and not significant. Exceptions were pH in the subsoil (- 0.06 pH units yr- 1) and ammonium in the middle soil horizon (- 0.086 mg N kg- 1 yr- 1). The linear trend is here defined as a linear combination of the true and unknown, but fixed spatial means. In quantifying the uncertainty of the estimated trend, only the sampling error in the estimated spatial means is accounted for. If there is a need to include uncertainty due to fluctuations of the true spatial means around a linear trend, then a super-population or time series model for the spatial means must be postulated which comprises a model error term. The linear trend is then defined as a model parameter, that can be estimated by Generalized Least Squares as in generalized linear modeling. --------------------------------------------------------------------------------

AB - This paper introduces and demonstrates design-based Generalized Least Squares (GLS) estimation of spatial means at selected time points from data collected in repeated soil surveys with partial overlap, such as a rotating and a supplemented panel. The linear time trend of the spatial means can then be obtained as a linear combination of the estimated spatial means. The GLS estimator is the minimum variance linear unbiased estimator. Five space–time designs were compared under a first order autoregressive time series model for the spatial means, through the average and generalized sampling variances of the estimated spatial means, and through the sampling variance of the estimated linear trend. None of the designs scored best on all three quality measures. If the aim of soil monitoring is estimation of both status and trend, then these two conflicting aims must be prioritized in order to choose an efficient space–time design. The methodology is demonstrated by a case study on eutrophication and acidification of forest soils. The linear trends in the spatial means of pH and the ammonium and nitrate concentrations at three depths in the soil profile, as estimated from a rotational design with four sampling times at an interval of one year, were small and not significant. Exceptions were pH in the subsoil (- 0.06 pH units yr- 1) and ammonium in the middle soil horizon (- 0.086 mg N kg- 1 yr- 1). The linear trend is here defined as a linear combination of the true and unknown, but fixed spatial means. In quantifying the uncertainty of the estimated trend, only the sampling error in the estimated spatial means is accounted for. If there is a need to include uncertainty due to fluctuations of the true spatial means around a linear trend, then a super-population or time series model for the spatial means must be postulated which comprises a model error term. The linear trend is then defined as a model parameter, that can be estimated by Generalized Least Squares as in generalized linear modeling. --------------------------------------------------------------------------------

U2 - 10.1016/j.geoderma.2011.06.001

DO - 10.1016/j.geoderma.2011.06.001

M3 - Article

VL - 164

SP - 172

EP - 180

JO - Geoderma

T2 - Geoderma

JF - Geoderma

SN - 0016-7061

IS - 3-4

ER -