Comparison of soil texture maps synthetized from standard depth layers with directly compiled products

• Published in: Geoderma 2019-10, Vol.352, p.360-372
• Main Authors: Laborczi, Annamária, Szatmári, Gábor, Kaposi, András Dezső, Pásztor, László
• Format: Article
• Language: English
• ISSN: 0016-7061; 1872-6259
• DOI: 10.1016/j.geoderma.2018.01.020

There are increasing demands nowadays on spatial soil information in order to support environmental related and land use management decisions. Physical soil properties, especially particle size distribution play important role in this context. A few of the requirements can be satisfied by the sand-, silt-, and clay content maps compiled according to global standards such as GlobalSoilMap (GSM) or Soil Grids. Soil texture classes (e. g. according to USDA classification) can be derived from these three fraction data, in this way texture map can be compiled based on the proper separate maps. Soil texture class as well as fraction information represent direct input of crop-, meteorological- and hydrological models. The model inputs frequently require maps representing topsoil features, which refer most commonly to 0–30 cm depth. This is covered by three consecutive depth intervals according to standard specifications: 0–5 cm, 5–15 cm, 15–30 cm. Becoming GSM and SoilGrids the most detailed freely available spatial soil data sources, the common model users (e. g. meteorologists, agronomists, or hydrologists) would produce input map from (the thickness-weighted mean of) these three layers. However, if the basic soil data and proper knowledge is obtainable, a soil texture map targeting directly the 0–30 cm layer could be independently compiled. In our work we compared Hungary's soil texture maps compiled using the same reference and auxiliary data and inference methods but for differing layer distribution. We produced the 0–30 cm clay, silt and sand map as well as the maps for the three standard layers (0–5 cm, 5–15 cm, 15–30 cm). Maps of sand-, silt-, and clay content were computed through composite regression kriging applying Additive Log-Ratio (alr) transformation. In addition to the Hungarian Soil Information and Monitoring System as reference soil data, digital elevation model and its derived components, soil physical property maps, remotely sensed images, land use-, geological-, as well as meteorological data were applied as auxiliary variables. We compared the directly compiled and the synthetized clay-, sand content, and texture class maps by different tools. In addition to pairwise comparison of basic statistical features (histograms, scatter plots), we examined the spatial distribution of the differences. We quantified the taxonomical distances of the textural classes, in order to investigate the differences of the map-pairs. We concluded that the direct