Predicting soil moisture conditions across a heterogeneous boreal catchment using terrain indices

Abstract

Soil moisture has important implications for drought and flooding forecasting, forest fire prediction and water supply management. However, mapping soil moisture has remained a scientific challenge due to forest canopy cover and small-scale variations in soil moisture conditions. When accurately scaled, terrain indices constitute a good candidate for modelling the spatial variation of soil moisture conditions in many landscapes. In this study, we evaluated seven different terrain indices at varying digital elevation model (DEM) resolutions and user-defined thresholds as well as two available soil moisture maps, using an extensive field dataset (398 plots) of soil moisture conditions registered in five classes from a survey covering a (68 km2) boreal landscape. We found that the variation in soil moisture conditions could be explained by terrain indices, and the best predictors within the studied landscape were the depth to water index (DTW) and a machine-learning-generated map. Furthermore, this study showed a large difference between terrain indices in the effects of changing DEM resolution and user-defined thresholds, which severely affected the performance of the predictions. For example, the commonly used topographic wetness index (TWI) performed best on a resolution of 16 m, while TWI calculated on DEM resolutions higher than 4 m gave inaccurate results. In contrast, depth to water (DTW) and elevation above stream (EAS) were more stable and performed best on 1–2 m DEM resolution. None of the terrain indices performed best on the highest DEM resolution of 0.5 m. In addition, this study highlights the challenges caused by heterogeneous soil types within the study area and shows the need of local knowledge when interpreting the modelled results. The results from this study clearly demonstrate that when using terrain indices to represent soil moisture conditions, modelled results need to be validated, as selecting an unsuitable DEM resolution or user-defined threshold can give ambiguous and even incorrect results.