-
森林的持水性能是森林生态系统的重要功能之一,其凋落物层与土壤层调控着90%以上的降雨,是森林生态系统发挥水源涵养的主要作用层[1]。很多学者对不同林分类型的凋落物层和土壤层的持水性能开展了研究[2-4],如殷沙等[5]针对马尾松Pinus massoniana-木荷Schima superba不同比例混交林分析得出对马尾松纯林补植阔叶树种,有利于水源涵养功能的改善。吴迪等[6]对北京九龙山8种林分研究结果得出华北落叶松Larix principis-rupprechtii纯林和油松Pinus tabulaeformis-华北落叶松混交林水源涵养能力较好,樟子松Pinus sylvestris var. mongolica纯林则最差。石媛等[7]对河北雾灵山不同密度油松人工林研究得出密度为1 400株·hm-2的油松人工林水源涵养能力最强。陈莉莉等[8]针对不同密度油松人工林探究结果为油松人工林密度在2 617株·hm-2时水源涵养功能最好。高琛等[9]对北京松山不同密度丁香Syzygium aromaticum天然林分析得出丁香天然林在密度为1 210株·hm-2左右时水源涵养能力较好,而赵磊等[10]对江西大岗山不同密度杉木Cunninghamia lanceolata林研究得出,杉木林林分密度在3 500株·hm-2左右其水源涵养功能较好等。可以看出,不同学者研究的结论不尽相同,这可能主要是由于研究对象所在区域及林分特征不同所造成的。马尾松为针叶树种,具有抗旱、耐贫瘠、适应性强等特点,主要分布在中国的南方地区。由于受诸多历史因素的影响,至20世纪60年代,赣南是中国水土流失较严重的地区之一,尤其是该区的兴国县,水土流失范围广强度大,曾被称为“江南沙漠”。为改善生态环境,减少水土流失,20世纪60-90年代初赣南开展了大面积的马尾松飞播造林工作,飞播林占目前有林地面积的21.4%,而这些林分大部分以固土蓄水为主要经营目标。然而,由于飞播林特殊的造林方式,林分密度相差较大,密度的差异对林分凋落物及土壤持水性能的影响目前尚不清楚。因此,本研究以马尾松飞播造林具有代表性的兴国县为研究区,旨在通过对不同密度飞播马尾松林的凋落物和土壤的持水性能进行比较分析,为提高飞播马尾松林持水性能有效措施的制订提供科学依据。
HTML
-
兴国县(26°03′~26°42′N,115°01′~115°51′E)位居江西省的南部,赣州市的北部地区。海拔高为130 ~1 200 m,低山、丘陵为其主要地貌,土壤类型主要为红壤。气候是亚热带东南季风气候,年平均气温为18.9 ℃,年平均降水量为1 539.0 mm,无霜期为284 d。县域森林资源丰富,森林覆盖率为72.2%,植被类型主要有针叶林、常绿阔叶林等,其中飞播马尾松林面积达6.4万hm2,共占现有有林地面积的29.5%[11]。
-
由于研究区保存的飞播马尾松林多为中龄林,故以飞播马尾松中龄林为研究对象,依据其林木密度情况,分为Ⅰ组(900~1 500株·hm-2),Ⅱ组(1 500~2 100株·hm-2),Ⅲ组(2 100~2 700株·hm-2),Ⅳ组(2 700~3 300株·hm-2)等4种密度组,设立标准地3个·密度组-1,共12个标准地,标准地面积为0.04 hm2(20 m × 20 m)。标准地土壤类型均为红壤,成土母岩为花岗岩,坡向为半阳坡,坡位为中坡。草本植物主要是雀稗Paspalum thunbergii,铁芒萁Dicranopteris linearis和五节芒Miscanthus floridulus等,林下灌木主要是檵木Loropetalum chinense和胡枝子Lespedeza bicolor等。各密度组标准地基本概况见表 1。
密度组 平均林龄/a 平均胸径/cm 平均树高/m 林分平均密度/(株·hm-2) 郁闭度 土层厚度/cm 腐殖质层厚度/cm 坡度/(°) Ⅰ 26 10.8 9.3 1 250 0.4 78.5 3.3 22 Ⅱ 27 10.3 10.2 1 825 0.5 74.8 3.0 23 Ⅲ 24 9.1 9.5 2 475 0.6 81.7 3.2 18 Ⅳ 24 7.8 7.8 2 958 0.7 76.3 4.3 17 Table 1. Basic situation of different density groups
-
在各标准地的上部、中部、下部分别选取3个1 m × 1 m具有代表性的凋落物小样方,用收获法按未分解层(枯枝落叶保持原状,叶形完整,外表没有分解痕迹的凋落物层)、半分解层(大部分枯枝落叶已经粉碎,叶形不完整,分解成碎屑的凋落物层)分别收获凋落物并称鲜质量;将样品带回实验室内,置于80 ℃的烘箱内烘干测定干质量,计算其含水率,推算出1 hm2林地凋落物的储存量。选用室内浸泡的方法把烘干后的凋落物放入尼龙网袋中并搁置水中浸泡24 h,取出后称量,分别测定未分解层、半分解层的最大持水率、有效拦蓄率、最大持水量和有效拦蓄量等。计算公式为[5]M=k×m×10-2, kr=[(m1-m0)/m0]×100,krm=[(m2-m0)/m0]×100, krs=0.85krm-kr, W1=M×krm, W2=M×krs。其中:M为凋落物储存量(t·hm-2);k为换算系数,k =凋落物干质量(g)/凋落物鲜质量(g); m,m1,m0,m2为取样样方凋落物鲜质量、凋落物取样鲜质量、烘干质量、凋落物浸水24 h后的质量(g);kr,krm,krs为凋落物自然持水率、最大持水率、有效拦蓄率(%);W1,W2为凋落物最大持水量、有效拦蓄量(t·hm-2)。
-
在各标准地的上部、中部、下部各挖取1个典型土壤剖面,记录坡面的形态特征后,分别在0~10 cm,10~20 cm,20~40 cm,40~80 cm土层中的5 cm,15 cm,30 cm,60 cm处用环刀采集取样,取3次·层-1。将样品带回实验室内,测定其土壤含水率、容重、毛管孔隙度、非毛管孔隙度、总孔隙度等物理性质,并由此分层计算得出土壤的最大持水量、毛管持水量和有效持水量,计算公式分别为[5]Wt=10 000Pth, Wc=10 000Pch, Wo=10 000Poh。其中:Wt为土壤最大持水量(t·hm-2);Wo为土壤毛管持水量(t·hm-2);Wc为土壤非毛管持水量(t·hm-2);Pt为土壤总孔隙度(%);Po为土壤毛管孔隙度(%);Pc为土壤非毛管孔隙度(%);h为土层厚度(m)。
-
运用Excel 2003进行数据处理及绘图,运用SPSS 17.0软件进行单因素方差分析。
2.1. 标准地设置
2.2. 凋落物储存量及其水分含量测定
2.3. 土壤物理性质测定及其持水量计算
2.4. 数据分析
-
由图 1可知:所有标准地林分凋落物未分解层、半分解层及总储存量(未分解层+半分解层)的范围分别为0.47~1.34 t·hm-2,0.93~2.92 t·hm-2,1.40~4.26 t·hm-2。半分解层及总储存量表现为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ,未分解层储存量表现为密度组Ⅱ>密度组Ⅰ>密度组Ⅲ>密度组Ⅳ,各密度组半分解层储存量均大于未分解层。方差分析表明:未分解层和总凋落物储存量不同密度组间无显著差异(P>0.05),但密度组Ⅱ的总凋落物储存量明显高于其他密度组;半分解层除密度组Ⅱ与Ⅳ表现出显著差异外(P<0.05),其他各密度组间均无显著差异(P>0.05)。
-
由表 2可知:4种密度组凋落物的平均最大持水率为116.22%~267.21%,总最大持水量为1.80~11.55 t·hm-2,均为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ;最大持水率、最大持水量均为半分解层大于未分解层,且都表现为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ。方差分析发现,未分解层中,密度组Ⅱ和密度组Ⅲ的最大持水率显著高于其他密度组(P<0.05),密度组Ⅱ与密度组Ⅲ之间最大持水率无显著差异(P>0.05),而最大持水量则是密度组Ⅱ显著高于其他密度组(P<0.05);半分解层中,密度组Ⅱ除与密度组Ⅲ最大持水率无显著差异外(P>0.05),与其他各密度组均有显著差异(P<0.05),最大持水量也是密度组Ⅱ显著高于其他密度组(P<0.05);就整体凋落物层而言,平均最大持水率、总最大持水量均表现出密度组Ⅱ显著高于其他密度组(P<0.05)。
密度组 凋落物层 最大持水率/% 最大持水量/(t·hm-2) 有效拦蓄率/% 有效拦蓄量/(t·hm-2) 未分解层 141.53±21.13 b 1.10±0.58 b 120.06±17.96 b 0.94±0.49 b Ⅰ 半分解层 166.21±29.52 bc 2.39±0.58bc 141.00±25.08 bc 2.02±0.49 bc 平均或合计 153.87±25.84 b 3.49±1.09 b 130.53±22.28 b 2.96±0.93 b 未分解层 259.16±20.92 a 3.53±0.76 a 219.79±17.87 a 2.99±0.65 a Ⅱ 半分解层 275.26±8.92 a 8.02±0.53 a 233.42±7.55 a 6.80±0.45 a 平均或合计 267.21±12.27 a 11.55±1.13 a 226.61±13.27 a 9.79±0.96 a 未分解层 220.35±25.28 a 1.64±0.79 ab 186.94±31.72 a 1.40±0.67 ab Ⅲ 半分解层 241.26±37.32 ab 3.94±0.88 b 204.78±21.53ab 3.34±0.75 b 平均或合计 230.81±30.16 b 5.58±1.55 b 195.86±25.16 b 4.74±1.32 b 未分解层 98.77±6.55 b 0.47±0.18 b 83.80±5.55 b 0.40±0.15 b Ⅳ 半分解层 133.66±19.90 c 1.33±0.45 c 113.41±16.89 c 1.13±0.38 c 平均或合计 116.22±13.36 b 1.80±0.43 b 98.61±13.35 b 1.53±0.37 b 说明:平均值±标准误,“平均或合计”中最大持水率、有效拦蓄率为平均值,最大持水量、有效拦蓄量为合计;纵向不同小写字母表示相同分解层不同密度组均值间差异显著(P<0.05)。 Table 2. Water holding characteristics of litter in different density groups
由表 2可知:4种密度组凋落物平均有效拦蓄率、有效拦蓄量均以密度组Ⅱ最高,分别为226.21%和9.79 t·hm-2。经多重比较发现,凋落物有效拦蓄率和有效拦蓄量无论是不同分解层还是就整个凋落物层而言,其结果均与相应层次的凋落物最大持水率和最大持水量分析结果相一致。
综上所述,4种密度组的凋落物储存量、持水与蓄水功能均表现为密度组Ⅱ(1 500~2 100株·hm-2)优于其他密度组。这可能是由于密度过低其凋落物储存量较少不利于持水与蓄水,而密度过高,林分郁闭度增加,光照不足,温度降低,影响了凋落物的分解速度,使凋落物层的持水性能降低,不利于飞播马尾松林的水源涵养和水土保持。
-
由表 3可知:0~80 cm土壤容重平均值范围是1.29~1.41 g·cm-3,表现为密度组Ⅳ>密度组Ⅰ>密度组Ⅲ>密度组Ⅱ,土壤毛管孔隙度、非毛管孔隙度、总孔隙度平均值范围分别是29.45%~35.87%,4.57%~8.55%和47.92%~51.32%,表现为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ。方差分析发现,不同密度组间相同土层的土壤容重、毛管孔隙度、总孔隙度均无显著差异(P>0.05),非毛管孔隙度除40~80 cm土层密度组Ⅱ与密度组Ⅳ有显著差异外(P<0.05),其他相同土层不同密度组间无显著差异(P>0.05),0~80 cm土层土壤容重、毛管孔隙度均值无显著差异(P>0.05),非毛管孔隙度均值表现为密度组Ⅱ除与密度组Ⅲ无显著差异外(P>0.05),与其他各密度组均有显著差异(P<0.05),总孔隙度均值则表现为密度组Ⅱ与密度组Ⅳ有显著差异(P<0.05)。
密度组 土层厚度/cm 土壤容重/(g·cm-3) 毛管孔隙度/% 非毛管孔隙度/% 总孔隙度/% Ⅰ 0〜10 1.32 ± 0.02 a 34.54 ± 4.73 a 6.53 ± 1.42 a 50.19 ± 5.52 a 10〜20 1.35 ± 0.05 a 33.81 ± 5.06 a 6.06 ± 0.94 a 49.06 ± 5.45 a 20〜40 1.38 ± 0.01 a 32.27 ± 5.23 a 5.00 ± 0.76 a 47.92 ± 5.05 a 40〜80 1.41 ± 0.03 a 32.05 ± 3.73 a 4.17 ± 1.39 ab 46.79 ± 4.85 a 平均 1.36 ± 0.02 a 33.17 ± 2.04 a 5.44 ± 0.57 bc 48.49 ± 2.36 ab Ⅱ 0〜10 1.21 ± 0.02 a 36.08 ± 2.33 a 9.82 ± 3.81 a 54.34 ± 4.13 a 10〜20 1.25 ± 0.04 a 35.98 ± 0.49 a 9.67 ± 2.69 a 52.83 ± 3.08 a 20〜40 1.31 ± 0.08 a 35.85 ± 0.96 a 7.68 ± 1.95 a 50.57 ± 2.85 a 40〜80 1.39 ± 0.10 a 35.56 ± 2.72 a 7.01 ± 0.60 a 47.55 ± 2.36 a 平均 1.29 ± 0.03 a 35.87 ± 0.80 a 8.55 ± 1.15 a 51.32 ± 2.12 a Ⅲ 0〜10 1.28 ± 0.09 a 36.57 ± 5.34 a 11.38 ± 1.58 a 51.70 ± 3.84 a 10〜20 1.29 ± 0.06 a 35.76 ± 6.32 a 7.09 ± 2.65 a 51.32 ± 4.24 a 20〜40 1.33 ± 0.10 a 33.80 ± 6.66 a 6.37 ± 1.29 a 49.81 ± 5.72 a 40〜80 1.39 ± 0.02 a 32.32 ± 3.96 a 6.06 ± 1.19 ab 47.55 ± 3.47 a 平均 1.32 ± 0.03 a 34.61 ± 2.47 a 7.72 ± 0.99 ab 50.09 ± 2.55 ab Ⅳ 0〜10 1.32 ± 0.09 a 32.27 ± 1.49 a 5.34 ± 0.93 a 50.19 ± 2.32 a 10〜20 1.39 ± 0.04 a 31.97 ± 4.66 a 4.93 ± 0.53 a 47.55 ± 3.62 a 20〜40 1.40 ± 0.09 a 28.71 ± 3.13 a 4.49 ± 0.30 a 47.17 ± 2.83 a 40〜80 1.41 ± 0.08 a 24.84 ±10.87 a 3.50 ± 0.54 b 46.79 ± 2.84 a 平均 1.41 ± 0.08 a 29.45 ± 2.78 a 4.57 ± 0.33 c 47.92 ± 2.65 b 说明:平均值±标准误,“平均或合计”中最大持水率、有效拦蓄率为平均值,最大持水量、有效拦蓄量为合计;纵向不同小写字母表示相同分解层不同密度组均值间差异显著(P<0.05)。 Table 3. Soil bulk density and porosity of different density groups
-
由图 2可知:不同密度组0~80 cm土壤非毛管持水量、毛管持水量、最大持水量的范围分别是332.62~628.93,2 210.17~2 859.98,3 792.45~3 984.91 t·hm-2,均表现为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ。分析土壤不同密度组的非毛管持水量、毛管持水量和最大持水量发现,各密度组之间无显著差异(P>0.05)。
综上分析,4种不同密度组中,土壤物理性质和土壤蓄水能力均表现为在密度组Ⅱ时处于最大值,表明飞播马尾松林林木密度为1 500~2 100株·hm-2能较好地改善土壤结构,提高土壤持水能力。
3.1. 凋落物储存量及持水性能
3.1.1. 凋落物储存量
3.1.2. 凋落物持水性能
3.2. 土壤物理性质及持水性能
3.2.1. 土壤容重及孔隙度
3.2.2. 土壤持水性能
-
4种密度组凋落物总储存量为1.40~4.26 t·hm-2,表现为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ,以中等密度1 500~2 100株·hm-2的最大,且各密度组凋落物储存量均表现为半分解层大于未分解层。赵芳等[12]对飞播马尾松林不同林下植被类型凋落物及土壤水文效应研究以及张伟等[13]对冀北山地不同林分凋落物蓄积量研究等也表明:凋落物中半分解层大于未分解层,而刘宇等[14]对不同密度华北落叶松人工林研究结果也得出,中等密度林分的凋落物储存量最多。本研究得出凋落物总储存量随林分密度的增加表现为先增大后减小的趋势。究其原因,可能是凋落物最终的储存量是由凋落物总体的量和它分解的速度来决定的[12],由于林木密度的大小会对林木、林下植被的凋落物量产生影响,而且通过影响光照、温湿度等因子也会对凋落物的分解速度等产生影响[15-17],因此凋落物储存量并不是随着密度的增大而增大。
凋落物的持水能力与凋落物储存量有密切关系,储存量越大,其持水能力越强。本研究表明:凋落物总最大持水量、平均最大持水率的变化与凋落物储存量表现出同样的规律,即随林分密度的增加呈现先增大后减少的趋势,以密度组Ⅱ最大,而未分解层与半分解层的最大持水量和最大持水率也表现出相同的规律,且半分解层均大于未分解层。曾建军等[18]研究城市水源5种森林凋落物水文效应也得出类似的结果。一般情况下,凋落物的实际拦蓄能力是通过有效拦蓄量的大小来衡量的,飞播马尾松林凋落物的有效拦蓄量随密度的增加出现先增大后减小的变化趋势。鲁邵伟等[19]对北京山地不同密度侧柏Platycladus orientalis人工林凋落物变化以及杜阿朋等[20]对不同林分密度邓恩桉Eucalyptus dunnii凋落物持水特性的研究结果均得出同样的规律,但余林等[21]对不同密度毛竹Phyllostachys edulis林凋落物层水文特性研究结果则表明:最大密度林分(3 000株·hm-2)的凋落物持水特性最大。由此可见,虽然凋落物的持水功能与凋落物的储存量有密切关系,但随林分密度的变化规律不尽相同。这可能是由于林分类型、地形环境等因素差异所造成。
-
本研究中0~80 cm土壤容重平均值范围是1.29~1.41 g·cm-3,表现为密度组Ⅳ>密度组Ⅰ>密度组Ⅲ>密度组Ⅱ,且同一密度组的土壤容重随着土壤深度的增加而增加,而土壤毛管孔隙度、非毛管孔隙度、总孔隙度变化则与土壤容重的变化相反,这与LI等[22]对滇池流域4种森林生态系统水土保持能力、石媛等[7]对河北雾灵山油松Pinus tabuliformis人工林土壤物理性质以及殷沙等[5]对马尾松-木荷的不同比例混交林等研究所得出的结果类似,主要是因为随着土层深度的增加,土壤的有机物质含量、土壤团聚性降低,导致土壤的紧实性增加,也有可能是由于土壤表面层凋落物合成成分、凋落物储存量及其分解速率、土壤地下根系成长发育情况等差异造成的[23]。
土壤持水性能的大小主要是由土壤的非毛管孔隙度来决定,而孔隙度与土壤本身的透气性关系密切。本研究得出土壤非毛管持水量、毛管持水量以及最大持水量均表现为密度组Ⅱ>密度组Ⅲ>密度组Ⅰ>密度组Ⅳ,与土壤容重的变化规律相反,这也说明了土壤容重越小,土壤疏松多孔,其结构性就越好,从而增强土壤持水能力。土壤非毛管持水量、毛管持水量以及最大持水量各密度组间虽无显著差异(P>0.05),但均表现为密度组Ⅱ最大,这一结果与肖欣等[24]研究的同等密度下飞播马尾松土壤持水性能的结论相一致。
综上所述,飞播马尾松林密度为1 500~2 100株·hm-2的林分,枯枝落叶的储存量大,能够较好地改善土壤结构,使得其凋落物层和土壤层的持水性能都明显优于其他密度组。因此,在今后的飞播马尾松林的经营管理和水源涵养林建设中,可考虑对飞播马尾松林的密度进行控制,以改善林地土壤结构,提高其持水性能。由于条件所限,本研究没有对林冠层和林下植被层的水源涵养功能进行调查分析。在今后研究中有待将林冠层、林下植被层、凋落物层、土壤层等4方面结合起来,以较全面了解其水源涵养功能。