-
土壤密度是土壤的基本物理性质之一,对坡面土壤水分运动有重要影响[1]。受气候(如降雨)、土壤质地、地形和植被类型等差异的影响,土壤密度在坡面上有很大的空间异质性。准确探明坡面土壤密度空间变异及其影响因素对于深入理解坡面产流机制和土壤水分的坡面分配有重要意义[2]。
土壤密度空间变异研究涉及到了区域、流域、坡面和样地等不同的空间尺度[3-7]以及草地、灌丛和森林等不同的植被类型[8-11],但关于坡面土壤密度空间变异的研究,主要集中在草地[8]和灌丛[10]。现有的森林坡面研究多关注土壤密度的坡位差异[12],但通过几个坡位的简单比较还不能完整揭示坡面微地形、立地环境等各差异下的土壤密度空间变异特征,这是因为即使在特定坡位的土壤密度也存有较大空间变异[11]。因此,森林坡面土壤密度空间变异研究还需细化。影响土壤密度空间变异的因素在不同空间尺度存在差异,在坡面上,地形(如坡度)、植被类型和土壤质地(结构)等是影响坡面土壤密度空间差异的重要因素[13]。但在不同地区、植被类型下,其主导因素也有差别,如影响土壤密度空间变异的主导因素在河北坝上地区草地坡面为地形[8],在四川西缘山地坡面是植被[14],在喀斯特坡面是地质[10]。作为黄土高原的产流核心区-土石山林区,土壤密度异质性是影响坡面产流空间变化的重要因素,但坡面土壤密度空间变异的影响因素未有探索。另外,多数研究关注表层或特定深度土壤密度的空间变异[10, 15-16],但受不同土层石砾含量、根系生物量等差异影响[17-18],不同土层土壤密度空间变异可能存有差异,但目前关于土壤密度空间变异随土层深度变化还缺乏足够探索。
为克服坡面土壤密度空间变异的影响,往往需进行高密度的取样来确保坡面土壤密度的评估精度,但费时、费力[19]。若能量化土壤密度合理取样数及与空间变异的关系,将显著提高坡面土壤取样的效率和评估精度。目前有关土壤密度合理取样数的研究较少,刘帆等[11]虽量化了华北落叶松林表层土壤密度的合理取样数,但这仅是样地尺度的探索,且未探讨土壤密度空间变异与其合理取样数的关系。
华北落叶松(Larix gmelinii var. principis-rupprechtii (Mayr) Pilg.)是黄土高原土石山区主要的造林树种,对当地的水源涵养、保持水土有重要作用。理清华北落叶松林坡面土壤密度空间变异及合理取样数对于准确理解华北落叶松林坡面产流空间差异的影响机制具有重要价值。为此,本研究以六盘山半湿润区的华北落叶松林坡面为研究对象,结合经典统计学和地统计学方法,明确华北落叶松林坡面0~100 cm土层土壤密度的空间异质性及影响因素,并采用Monte Carlo模拟重抽样方法,量化各层土壤密度的合理取样数。
-
由表1可知,随土层加深,土壤密度逐渐增加,0~20、20~40、40~60、60~80和80~100 cm土层土壤密度平均值(变化范围)分别为0.88(0.73~1.08)、1.04(0.78~1.55)、1.23(0.91~1.66)、1.41(1.02~1.78)和1.49(1.09~1.73)g·cm−3,变异系数为0.09、0.16、0.17、0.15和0.11,土壤密度在0~20 cm土层为弱变异,在其它土层为中等变异。由图1可知,土壤密度空间变异系数与石砾含量的空间变异系数呈显著的(P<0.05)线性相关,表明土壤密度空间变异随土层深度的变化主要取决于石砾含量的空间异质性。
土层深度
Soil depth/
cm最小值
Minimum/
(g·cm−3)最大值
Maximum/
(g·cm−3)平均值
Mean/
(g·cm−3)方差
Variance标准差
Standard deviation变异系数
CV0~20 0.73 1.08 0.88 0.01 0.08 0.09 20~40 0.78 1.55 1.04 0.03 0.17 0.16 40~60 0.91 1.66 1.23 0.04 0.21 0.17 60~80 1.02 1.78 1.41 0.04 0.20 0.15 80~100 1.09 1.73 1.49 0.03 0.16 0.11 0~100 0.99 1.48 1.21 0.01 0.12 0.10 Table 1. Descriptive statistics characteristics of soil bulk density on the slope
-
由表2可知,土壤密度的半方差函数在0~20和20~40 cm土层为高斯模型,在40~60和80~100 cm土层为球状模型,在60~80 cm土层为指数模型。各层土壤密度块金值均为正值;土壤密度基台值沿土层深度的变化范围0.011~0.060,在0~20 cm土层最小,在40~60 cm土层最大。0~20、20~40、40~60、60~80和80~100 cm土层的块基比分别为0.295、0.348、0.300、0.495和0.498,表明各层土壤密度均属中等空间自相关性。各层土壤密度变程分别为257.70、126.20、179.30、215.42和294.50 m,表明随土层加深,空间自相关范围先减小后增大,在80~100 cm土层最大,在20~40 cm土层最小。
土层深度
Soil depth/cm块金值C0
Nugget基台值
Still块基比
Nugget/Still变程
Range/m决定系数
R2模型类型
Model type0~20 0.003 38 0.011 0.295 257.70 0.89 高斯模型 20~40 0.016 55 0.048 0.348 126.20 0.94 高斯模型 40~60 0.018 10 0.060 0.300 179.30 0.95 球状模型 60~80 0.027 34 0.055 0.495 215.42 0.65 指数模型 80~100 0.014 76 0.030 0.498 294.50 0.72 球状模型 Table 2. Semi-variogram theoretical models of soil bulk density in different soil layers and their related parameters
-
由图2可知,坡面各层土壤密度有明显的空间异质性,自坡顶至坡底,土壤密度总体上呈“低~高~低”的变化。各层土壤密度的最大值基本在坡中,最小值在坡下(0~20 cm除外,在坡上)。0~20 cm土层土壤密度跨度范围较小,斑块较大,空间分布较均匀;40~60 cm土层土壤密度跨度范围较大,条带密集,空间分布格局最复杂。随土层加深,土壤密度空间分布从斑块状逐渐变为条带状。
-
基于Monte Carlo模拟重抽样的方法确定了在95%的置信区间下,若控制坡面土壤密度在10%误差内的最少取样数。由图3可知,随取样数增加,土壤密度的平均值和置信区间逐渐向内收缩。在95%置信区间下,当取样数分别为4、10、10、8和5个时,得到的0~20、20~40、40~60、60~80和80~100 cm土层土壤密度可控制在10%误差内。由图4可知,土壤密度的最少取样数与其变异系数呈极显著的线性相关(P<0.01),即随土壤密度空间变异系数的增加,最少取样数逐渐增加。
Spatial Heterogeneity and Reasonable Sample Sizes of Soil Bulk Density on A Larch Plantation Hillslope in Liupan Mountains
- Received Date: 2021-12-08
- Accepted Date: 2022-01-11
- Available Online: 2022-08-20
Abstract: