Combining two national‐scale datasets to map soil properties, the case of available magnesium in England and Wales

Abstract

Summary Given the costs of soil survey it is necessary to make the best use of available datasets, but data that differ with respect to some aspect of the sampling or analytical protocol cannot be combined simply. In this paper we consider a case where two datasets were available on the concentration of plant‐available magnesium in the topsoil. The datasets were the Representative Soil Sampling Scheme (RSSS) and the National Soil Inventory (NSI) of England and Wales. The variable was measured over the same depth interval and with the same laboratory method, but the sample supports were different and so the datasets differ in their variance. We used a multivariate geostatistical model, the linear model of coregionalization (LMCR), to model the joint spatial distribution of the two datasets. The model allowed us to elucidate the effects of the sample support on the two datasets, and to show that there was a strong correlation between the underlying variables. The LMCR allowed us to make spatial predictions of the variable on the RSSS support by cokriging the RSSS data with the NSI data. We used cross‐validation to test the validity of the LMCR and showed how incorporating the NSI data restricted the range of prediction error variances relative to univariate ordinary kriging predictions from the RSSS data alone. The standardized squared prediction errors were computed and the coverage of prediction intervals (i.e. the proportion of sites at which the prediction interval included the observed value of the variable). Both these statistics suggested that the prediction error variances were consistent for the cokriging predictions but not for the ordinary kriging predictions from the simple combination of the RSSS and NSI data, which might be proposed on the basis of their very similar mean values. The LMCR is therefore proposed as a general tool for the combined analysis of different datasets on soil properties. Highlights Differences in sample support mean that two datasets on a soil property cannot be combined simply. We showed how a multivariate geostatistical model can be used to elucidate the relationships between two such datasets. The same model allows soil properties to be mapped jointly from such data. This offers a general basis for combining soil datasets from diverse sources

Spatial prediction of the concentration of selenium (Se) in grain across part of Amhara Region, Ethiopia

“…Specifically we considered a multivariate version of the model ( Marchant and Lark, 2007 ) which allows us to combine collocated and non-collocated data. The empirical best linear unbiased prediction (E-BLUP), based on the fitted model, has an associated prediction error distribution, and on the basis of this we were able to quantify uncertainties in the prediction relative to threshold concentrations of interest, and to use strategies to communicate this uncertainty which have been used elsewhere ( Lark et al, 2014 , Lark et al, 2019 ; Mastrandrea et al, 2010 ).…”

Detailed global modelling of soil organic carbon in cropland, grassland and forest soils

“…This method introduces error in cases where belowground plant growth does not scale with aboveground yield. There is some evidence that belowground C varies less between crops than aboveground C. For example, Cagnarini et al [35] conclude that the use of constant factors overestimates belowground biomass when compared with field-measured belowground biomass, for maize and wheat in Switzerland. Taghizadeh-Toosi et al [93] used long-term experimental data to conclude that belowground biomass growth was independent of the crop.…”

“…SOC) to the input variables. The Rothamsted Carbon Model (RothC) [30] model is one of the most commonly used soil models today [30–35]. The set of inputs required by RothC is a crucial advantage when compared with other process-based soil models that require much larger datasets [33].…”

A reconnaissance survey of farmers’ awareness of hypomagnesaemic tetany in UK cattle and sheep farms

“…Index values are available for many soil properties, guiding pasture management decisions [25]. The majority of respondents in this study were from Wales and England (86%) where soil plant-available Mg concentrations are generally greater than the Index 2 (51 mg L -1 ) which is used as a Mg-fertiliser recommendation threshold [28, 29]. It is noteworthy that grassland soils are often not optimally managed for pH, with 53% being below the recommended value of 6.0 in a recent private sector data synthesis [29].…”