Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

被引:173
作者
Kerry, R. [1 ]
Oliver, M. A.
机构
[1] Brigham Young Univ, Dept Geog, Provo, UT 84602 USA
[2] Univ Reading, Dept Soil Sci, Reading, Berks, England
关键词
variogram; residual maximum likelihood (REML); method of moments (MoM); prediction; sampling;
D O I
10.1016/j.geoderma.2007.04.019
中图分类号
S15 [土壤学];
学科分类号
0903 ; 090301 ;
摘要
It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites. (C) 2007 Elsevier B.V. All rights reserved.
引用
收藏
页码:383 / 396
页数:14
相关论文
共 25 条
[12]   MAXIMUM-LIKELIHOOD PARAMETER-ESTIMATION OF HYDROLOGIC SPATIAL PROCESSES BY THE GAUSS-NEWTON METHOD [J].
KITANIDIS, PK ;
LANE, RW .
JOURNAL OF HYDROLOGY, 1985, 79 (1-2) :53-71
[13]   Estimating variograms of soil properties by the method-of-moments and maximum likelihood [J].
Lark, RM .
EUROPEAN JOURNAL OF SOIL SCIENCE, 2000, 51 (04) :717-728
[14]  
MARDIA KV, 1984, BIOMETRIKA, V71, P135, DOI 10.2307/2336405
[15]  
Matheron G., 1965, Regionalized Variables and Their Estimation
[16]  
Matheron G., 1971, THEORY REGIONALIZED
[17]   THE DESIGN OF OPTIMAL SAMPLING SCHEMES FOR LOCAL ESTIMATION AND MAPPING OF REGIONALIZED VARIABLES .2. PROGRAM AND EXAMPLES [J].
MCBRATNEY, AB ;
WEBSTER, R .
COMPUTERS & GEOSCIENCES, 1981, 7 (04) :335-365
[18]   A SIMPLEX-METHOD FOR FUNCTION MINIMIZATION [J].
NELDER, JA ;
MEAD, R .
COMPUTER JOURNAL, 1965, 7 (04) :308-313
[19]   Inference of spatial indicator covariance parameters by maximum likelihood using MLREML [J].
Pardo-Iguzquiza, E .
COMPUTERS & GEOSCIENCES, 1998, 24 (05) :453-464
[20]   Maximum likelihood inference of spatial covariance parameters of soil properties [J].
Pardo-Iguzquiza, E ;
Dowd, PA .
SOIL SCIENCE, 1998, 163 (03) :212-219