In recent years there has been an increased demand for digital soil mapping (DSM) products. DSM has become
the ultimate soil spatial representation framework due to its quantitative results, replicability, and uncertainty
analysis. The present study aimed to map the probability distribution of the derived soil profiles (DSPs) of soil
typological units (STUs). DSPs are statistical representation of the properties of the soil profiles belonging to
STUs. STUs aggregate individual profiles into a group. The criteria used for grouping were homogeneity for
World Reference Base (WRB) reference soil group (WRB-RSG), qualifiers (WRB-qu), and Soil Taxonomy particle
size for the family classification (USDA-PS), and the belonging to a specific Soil Region. The European Soil
Regions have been suggested as the primary grouping criteria for soil mapping at the European continental scale
since they define continental-scale soilscapes, distinguished mainly by their climate and geology. To map DSPs,
we firstly mapped STUs. The grouping criteria of STUs were mapped at 500 m spatial resolution, using a Neural
Network trained on 18,707 georeferenced and analyzed soil profiles selected from the Italian national soil
database. A 10% of the soil profiles were randomly sampled using a stratified sampling approach for validation.
In particular, the procedure consisted of: i) mapping the grouping criteria WRB-RSG, WRB-qu, and USDA-PS, on a
500 m national grid, through Neural Network; ii) grouping soil profiles on the base of the combinations of
grouping criteria (WRB-RSG, WRB-qu, USDA-PS, and Soil Regions) as mapped with the first step at each grid
node, to produce a map of Soil Typological Units (STUs); iii) calculating statistics for the soil parameters of the
groups of soil profiles created, to produce a map of Derived Soil Profiles (DSPs). DSPs statistics (average, standard
deviation, and sample numerosity) were elaborated for the following parameters: soil rooting depth, pH (in
water), soil organic carbon, clay, silt, sand, coarse fragments, and cation exchange capacity. The maps obtained
were validated against the test set. The same test set was used for the comparison with the National benchmark
map (Soils Map of Italy 1:1,000,000) and with the global scale SoilGrids at 250 m spatial resolution. The overall
accuracy was 45.98% for the WRB-RSG map compared with the 30.74% of WRB-RSG as mapped with the Soil
Map of Italy, and 28.79% as mapped with SoilGrids; 33.07% for WRB-qu compared with the 15.69% of WRB-qu
as mapped with the Soil Map of Italy, and 12.45% as mapped with SoilGrids, and 45.48% for USDA-PS, not
comparable with the National and Global benchmarks. Tau statistics showed a higher accuracy Kappa of our
approach than in others, due to the unbalanced classes numerosity. The predictive ability in the validation of
DSPs parameters resulted in a R2 of 0.35 for clay (0.16 with SoilGrids), 0.28 for sand (0.08 with SoilGrids), 0.18
for pH in water (0.21 with SoilGrids). The proposed approach produced harmonized soil type maps with higher
accuracy than the previous generation of conventional field-based soil maps for the national benchmark and the
calculation of the uncertainty. The STUs express variability of soil properties between groups so their knowledge
might improve our understanding of the soil distribution, the planning of their management, monitoring, and the
decisions for further surveys. A future challenge will be including more dynamic parame
FANTAPPIÈ Maria;
L'ABATE Giovanni;
SCHILLACI Calogero;
COSTANTINI Edoardo;
2024-09-10
ELSEVIER
JRC130005
2352-0094 (online),
https://doi.org/10.1016/j.geodrs.2023.e00619,
https://publications.jrc.ec.europa.eu/repository/handle/JRC130005,
10.1016/j.geodrs.2023.e00619 (online),
