Spatial pattern and heterogeneity of soil properties in sand dunes under grazing and restoration in Horqin Sandy Land, Northern China

Applying a combination of classical and geostatistical methods, we identified soil properties and their spatial variation in a 5-year grazed sand dune (GSD5) and a 20-year recovered sand dune (RSD20) in Horqin Sandy Land, northern China. The paper assesses the effect of grazing, topography and vegetation restoration on spatial heterogeneity of soil properties. The results showed that soil organic carbon, total nitrogen, very fine sand (0.1–0.05 mm) content and their coefficients of variation were lower in GSD5 than in RSD20, while soil water contents (0–20 cm and 20–40 cm depths) were higher in GSD5 than in RSD20. Geostatistical analysis revealed that the spatial structured variance accounted for the largest proportion of total sample variance in soil properties at the measured scale under grazing and restoration. The spatial autocorrelation ranges were 66.30 m for soil organic carbon and 50.80 m for total nitrogen in GSD5 less than those in RSD20 (70.00 m and 76.10 m, respectively), while the spatial autocorrelation ranges of soil particle size fractions and soil water contents in RSD20 were less than those in GSD5. Kriging-interpolated maps also showed that the heterogeneity of soil organic carbon and total nitrogen and their degree of patch fragmentation were higher in GSD5 than in RSD20. These results suggested that continuous grazing resulted in an increase in spatial variability of soil nutrient and a decrease in spatial variability of soil particle size fractions and soil water content. Soil organic carbon and total nitrogen of sand dunes are associated closely with soil particle size fractions, relative height of sampling site and vegetation cover. Spatial patterns of soil properties are most strongly related to grazing, topography and plant-induced heterogeneity in sand dune ecosystems prone to wind erosion.