Reducing sampling intensity in order to investigate spatial variability of soil pH,organic matter and available phosphorus using co-kriging techniques. A case study of acid soils in Eastern Croatia

Soil pH, hydrolytic acidity (HA), organic matter (OM) and plant available phosphorus (AP) are factors controlling the environmental-friendly soil management in agroecosystems. These parameters are highly variable in space. The objective of this work is to study spatial variability of pH, HA, OM and AP using several interpolation methods in Eastern Croatia. A total of 1004 (0–30 cm) soil samples were collected, and several univariate and multivariate interpolation performances were tested. The results showed that soils of the study area had high HA and AP, while pH and OM were low. The variogram analysis revealed different spatial structures among studied soil properties and demonstrate a need for variable-rate management. Soil pH and OM had lower spatial variability compared to AP and HA. Ordinary kriging was the most accurate method to estimate the studied variables. The incorporation of auxiliary variables increased the precision of the estimations for HA. Soil AP and OM showed different results for spatial prediction obtained by co-kriging. Overall, the incorporation of pH as auxiliary variable increased the prediction of the models. However, more co-variates should be incorporated in further models, in order to identify with more precision areas that need to be restored.