Spatio-temporal variability of some metal concentrations in the soil of eastern England, and implications for soil monitoring.

Abstract Previous workers have proposed the use of multivariate geostatistics for the problem of estimating temporal change in soil properties for soil monitoring, but this has yet to be evaluated. We present a case study of this approach from the Humber–Trent region in North East England. We extracted data from two sources on cobalt, nickel and vanadium concentrations in the topsoil on two dates. Auto-variograms were estimated for each metal on each date, and pseudo cross-variograms for each metal on the two dates. It was shown that robust estimators of the auto and pseudo cross-variograms were needed for the analysis of these data. A linear model of coregionalization was then fitted to describe the spatio-temporal variability of each metal. While the concentration of each metal in the soil showed pronounced spatial dependence that we know is driven by parent material, the change over time was only spatially structured for cobalt and vanadium. This shows that information on spatial variability from a single date may be a poor guide to the design of a monitoring scheme. We showed how the cokriging variance of the change in concentration of cobalt and vanadium depends on sampling effort and strategy. The change in these particular variables between two dates is best estimated by sampling with equal intensity at the same sites on both dates; and when resampling an existing baseline survey it is best to sample them at rather than between the original sites. The best strategy in any case depends on how the variable is coregionalized over time.

[1]  R. Olea Geostatistics for Natural Resources Evaluation By Pierre Goovaerts, Oxford University Press, Applied Geostatistics Series, 1997, 483 p., hardcover, $65 (U.S.), ISBN 0-19-511538-4 , 1999 .

[2]  N. Cressie,et al.  Robust estimation of the variogram: I , 1980 .

[3]  H. Bourennane,et al.  Short‐range spatial variability of metal contents in soil on a one hectare agricultural plot , 2000 .

[4]  R. M. Lark,et al.  A comparison of some robust estimators of the variogram for use in soil survey , 2000 .

[5]  R. M. Lark,et al.  Fitting a linear model of coregionalization for soil properties using simulated annealing , 2003 .

[6]  Hannes Flühler,et al.  Temporal change of spatially autocorrelated soil properties: optimal estimation by cokriging , 1994 .

[7]  A. Journel,et al.  Geostatistics for natural resources characterization , 1984 .

[8]  A. McBratney,et al.  Optimal interpolation and isarithmic mapping of soil properties: V. Co-regionalization and multiple sampling strategy , 1983 .

[9]  E. Bringmark,et al.  IMPROVED SOIL MONITORING BY USE OF SPATIAL PATTERNS , 1998 .

[10]  Shenglu Zhou,et al.  Evaluation of spatial and temporal changes of soil quality based on geostatistical analysis in the hill region of subtropical China , 2003 .

[11]  David McGrath,et al.  Geostatistical and GIS analyses on soil organic carbon concentrations in grassland of southeastern Ireland from two different periods , 2004 .

[12]  A. McBratney,et al.  Choosing functions for semi‐variograms of soil properties and fitting them to sampling estimates , 1986 .

[13]  T. Astatkie,et al.  Spatial and temporal trends in soil properties at a soil quality benchmark site in central Nova Scotia , 2000 .

[14]  R. Lark,et al.  Geostatistics for Environmental Scientists , 2001 .

[15]  Marc G. Genton,et al.  Highly Robust Variogram Estimation , 1998 .

[16]  T. R. Lister,et al.  The influence of parent material on topsoil geochemistry in eastern England , 2003 .

[17]  P.F.M. van Gaans,et al.  Future trends, detectable by soil monitoring networks? , 1998 .

[18]  Alex B. McBratney,et al.  Further Comparison of Spatial Methods for Predicting Soil pH , 1990 .

[19]  D. Myers Pseudo-cross variograms, positive-definiteness, and cokriging , 1991 .

[20]  W. V. Duijvenbooden Soil monitoring systems and their suitability for predicting delayed effects of diffuse pollutants , 1998 .

[21]  R. Webster,et al.  Optimal interpolation and isarithmic mapping of soil properties. II. Block kriging. , 1980 .

[22]  Richard Webster,et al.  Estimating temporal change in soil monitoring: I. Statistical theory , 1995 .

[23]  Peter A. Dowd,et al.  The Variogram and Kriging: Robust and Resistant Estimators , 1984 .

[24]  P. J. Loveland,et al.  The soil geochemical atlas of England and Wales. , 1992 .

[25]  M. Hubert,et al.  A Comparison of Some New Measures of Skewness , 2003 .

[26]  H. Künsch,et al.  Generalized cross-covariances and their estimation , 1996 .

[27]  R. Webster,et al.  Sampling To Monitor Soil In England And Wales , 1999 .

[28]  D. Hirst,et al.  Temporal changes in soil properties at an upland Scottish site between 1956 and 1997. , 2001, The Science of the total environment.

[29]  P. Rousseeuw,et al.  Alternatives to the Median Absolute Deviation , 1993 .

[30]  Richard Webster,et al.  Estimating temporal change in soil monitoring: II. Sampling from simulated fields , 1995 .

[31]  R. Webster,et al.  Geostatistical analysis of soil contamination in the Swiss Jura. , 1994, Environmental pollution.

[32]  Yadolah Dodge,et al.  L 1 -Statistical Analysis and Related Methods. , 1993 .

[33]  R. M. Lark,et al.  Modelling complex soil properties as contaminated regionalized variables , 2002 .

[34]  R. M. Lark,et al.  Robust estimation of the pseudo cross‐variogram for cokriging soil properties , 2002 .

[35]  D. J. Brus,et al.  Random sampling or geostatistical modelling? Choosing between design-based and model-based sampling strategies for soil (with discussion) , 1997 .

[36]  H. Künsch,et al.  On the pseudo cross-variogram , 1993 .

[37]  Gerard B. M. Heuvelink,et al.  Modelling soil variation: past, present, and future , 2001 .

[38]  J. Andrew Royle,et al.  Space: Time Dynamic Design of Environmental Monitoring Networks , 1999 .

[39]  P. Rousseeuw,et al.  Explicit scale estimators with high breakdown point , 1992 .

[40]  R. Reese Geostatistics for Environmental Scientists , 2001 .