Uncertainty assessment for mapping changes in soil organic matter using sparse legacy soil data and dense new-measured data in a typical black soil region of China

A total of 89 legacy soil organic matter (SOM) data (0–20 cm) sourced from the soil survey in 1980 in Hailun County of Northeastern China and 268 SOM data (0–20 cm) measured in the area in 2008 were collected in this study. The most straightforward mapping methods such as ordinary kriging (OK), inverse distance weighting interpolation (IDW), and the soil type-based method (ST) were used to derive the spatial patterns of changes in SOM; and the uncertainties associated with the mapping methods were quantitatively assessed through the uncertainty intervals of SOM realizations that were generated by sequential Gaussian simulation and/or Monte Carlo simulation. Results showed that the overall trend of changes in SOM over the last 28 years in the area decreased. The spatial patterns of changes in SOM derived by OK₂₀₀₈–OK₁₉₈₀(abbreviations before and after the dash represent the mapping method for SOM data in 2008 and 1980, respectively) and IDW₂₀₀₈–OK₁₉₈₀, OK₂₀₀₈–IDW₁₉₈₀and IDW₂₀₀₈–IDW₁₉₈₀, OK₂₀₀₈–ST₁₉₈₀and IDW₂₀₀₈–ST₁₉₈₀are similar. The change map obtained by the ST₂₀₀₈–ST₁₉₈₀method is also similar to those derived by OK₂₀₀₈–ST₁₉₈₀and IDW₂₀₀₈–ST₁₉₈₀methods, but local details of changes in SOM presented by the ST₂₀₀₈–ST₁₉₈₀method are coarser than those obtained by the OK₂₀₀₈–ST₁₉₈₀and IDW₂₀₀₈–ST₁₉₈₀methods. The OK₂₀₀₈–ST₁₉₈₀method, with an uncertainty tradeoff compared to the OK₂₀₀₈–OK₁₉₈₀and ST₂₀₀₈–ST₁₉₈₀methods, can be considered as the most suitable method for mapping changes in SOM in the area due to the limited legacy soil survey data and inaccurate sample location information records in the area.