Comparing Two Geostatistical Simulation Algorithms for Modelling the Spatial Uncertainty of Texture in Forest Soils

Uncertainty assessment is an essential part of modeling and mapping the spatial variability of key soil properties, such as texture. The study aimed to compare sequential Gaussian simulation (SGS) and turning bands simulation (TBS) for assessing the uncertainty in unknown values of the textural fractions accounting for their compositional nature. The study area was a forest catchment (1.39 km2) with soils classified as Typic Xerumbrepts and Ultic Haploxeralf. Samples were collected at 135 locations (0.20 m depth) according to a design developed using a spatial simulated annealing algorithm. Isometric log-ratio (ilr) was used to transform the three textural fractions into a two-dimensional real vector of coordinates ilr.1 and ilr.2, then 100 realizations were simulated using SGS and TBS. The realizations obtained by SGS and TBS showed a strong similarity in reproducing the distribution of ilr.1 and ilr.2 with minimal differences in average conditional variances of all grid nodes. The variograms of ilr.1 and ilr.2 coordinates were better reproduced by the realizations obtained by TBS. Similar results in reproducing the texture data statistics by both algorithms of simulation were obtained. The maps of expected values and standard deviations of the three soil textural fractions obtained by SGS and TBS showed no notable visual differences or visual artifacts. The realizations obtained by SGS and TBS showed a strong similarity in reproducing the distribution of isometric log-ratio coordinates (ilr.1 and ilr.2). Overall, their variograms and data were better reproduced by the realizations obtained by TBS.