Soil aggregate stability within morphologically diverse areas

Knowledge of spatial distribution of soil aggregate stability as an indicator of soil degradation vulnerability and its possible prediction are required for many scientific and practical environmental studies. The goal of our study was to provide a model for predicting soil aggregate stability within morphologically diverse areas, where soil properties have been affected by soil material redistribution due to erosion. The study was performed on a study site (6ha area) in the loess region of Southern Moravia, Czech Republic. Haplic Chernozem, which is an original dominant soil unit, has been transformed into different soil units (eroded phases of Chernozem, Regosol, colluvial Chernozem and Colluvial soil). 36 sampling spots were selected in order to represent diverse soils. The following soil properties were measured: oxidable organic carbon content (Cox), CaCO3 content, pHH2O, pHKCl, soil particle density (ρs), bulk density (ρd), porosity (P), actual field soil-water content (θfield), content of iron and manganese (in ammonium oxalate extract, Feo and Mno, and dithionite–citrate extract, Fed and Mnd) and mass specific magnetic susceptibility (χlf and χhf). The aggregate stability was assessed using various tests to study different disruption mechanisms. Terrain attributes were derived from a digital elevation model.In general, the lowest soil aggregate stability was observed on steep slopes, which were highly impacted by soil erosion. The highest aggregate stability was measured on soils sampled at relatively flat upper parts, which were less influenced by erosion processes. Higher stability was also obtained on toe slopes, where the sedimentation of previously eroded soil material occurred. The simple correlations revealed that characteristics resulting from the tests studying aggregate slaking due to the compression of the entrapped air (Water Stable Aggregate index and coefficient of vulnerability from fast wetting test) were positively impacted by the Cox, P, Feo, Mno, Fed, Mnd, χlf and χhf values, and negatively by the ρd value. The soil aggregate stability was also negatively influenced by the plan and total terrain curvatures, i.e. larger aggregate stability was measured at concave parts in comparison with that at convex parts. Almost no statistically significant relationships were found in the case of the tests evaluating either aggregate disintegration caused by the micro-cracking due to the different swelling, or by the physico-chemical dispersion due to the osmotic stress or the mechanical aggregate breakdown. The multiple linear regressions resulted in the model for estimating the WSA index using the Cox content, total terrain curvature and actual field soil-water content (θfield). In this model the Cox content positively and the total terrain curvature and θfield value negatively influenced the value of the WSA index. Since Cox was positively related with iron content and thus also with the magnetic susceptibility, the alternative model was proposed for less costly and time consuming WSA estimation. The WSA index may be predicted by combining the mass specific magnetic susceptibility (χlf and χhf), total terrain curvature and actual field soil-water content (θfield).