南亚热带不同类型人工林的水土保持效应

孙 聃，刘素萍，申卫军，饶兴权，郭志峰，李永强

我国南方自20世纪80年代以来大规模植造人工林，以华南地区为例，山地丘陵占土地面积的70%以上，人工林种植面积为3 213公顷，水土流失面积达3 533 hm2，通过种植人工林改变地表径流和土壤侵蚀特征的研究显得尤为重要。森林影响径流的过程是一个十分复杂的过程，不同的树种和植被群落结构，不同的森林年龄在地区、季节都会产生不同效应。森林植被变化对径流的影响幅度相差较大，森林覆盖度的降低可以不同程度地增加流域的产水量，造林则导致流域产水量降低。但森林植被变化对流域产水量的影响却不近相同［1-3］，主要有3种不同结论:1)森林的存在会使径流量增加;2)森林的存在与径流量之间没有明显的关系;3)森林的存在会减少年径流量。

为了研究水土流失过程，许多国家和地区已经设立了小区径流场［4］，中国过去几十年中大量的水土流失研究也都使用了小区径流场［5-6］。尽管径流场的大小、保护措施和测量时长有所不同，但其数据对中国的水土流失提供了宝贵的信息，也揭示了不同植被类型与水土侵蚀之间存在的规律［7-12］。

受季风气候的强烈影响，我国南亚热带地区拥有丰富的森林动植物种类和资源，植物区系、植被类型以及稀有濒危植物的分布也存在着地带分异性。我国南亚热带也有大面积的各种人工林和次生林。人工林的主要造林树种以马尾松(Pinus massoniana)、杉木(Cunninghamia lanceolata)、马占相思(Castanopsis hystrix)、桉树(Eucalyptus spp)及一些乡土树种为主［13］。本文将结合9种植被类型和6年的观测数据来分析水土流失的动态过程，研究植被类型和林龄如何影响着水土保持，以期为华南地区人工林造林提供科学依据。

1 研究区概况

鹤山丘陵综合试验站是中国科学院华南植物研究所(园)与鹤山市林业科学研究所共建，2005年成为国家野外台站之一。该站位于广东省江门市鹤山市，地处南亚热带丘陵地区，平均海拔高度80 m，顶极森林群落是亚热带季风常绿阔叶林。年均温度22.5℃，最热为7月份，最冷为1月份。年均降雨量1 534 mm，干湿季交替明显。土壤为砂页岩发育而成的赤红壤。鹤山站先后建立成熟人工林和幼龄人工林2个研究样地。成熟人工林样地1984年建于桃源镇，位于 E 112°54'，N 22°41'，1986 年种植人工林种类包括:1)马占相思林AM(Acacia mangium);2)桉树林EU(Eucalyptus urophylla);3)荷木林SS(Schima superba);4)马尾松林PM(Pinus massoniana);幼龄人工林样地2005年建于共和镇，位于E 112°50'，N 22°34'各样地约 1 hm2。2005 年幼龄人工林样地包括:1)10种树种混交M1(10 species mixed forest);2)30种树种混交M2(30 species mixed forest);3)尾叶桉纯林EU1(Eucalyptus urophylla);4)厚荚相思纯林AC(Acacia crassicarpa);5)红锥纯林CH(Castanopsis hystrix);6)草坡GL本实验地立地条件相似，林下物种主要为芒萁(Dicranopteris pedata)。

2 材料与方法

2.1 数据采集

本研究采用小集水区径流场，把生态系统边界定义在小集水区的边界以内，我们分别在4个成熟人工林和5个幼龄人工林中建立了地表径流和含沙量监测系统。在每个林型内部用四周插入水泥板至20 cm深土层，插入后将原土回填，水泥板高出地面10 cm，使地表径流导入测流堰中，径流进入堰中通过三角堰口流出，也是径流的唯一出口。在坡底边设1个径流槽和1个仪器室，通过WGZ－1型自动水位计记录堰内水位高度从而换算出径流量。幼龄人工林试验区每个径流场投影面积大小为15 m×20 m，成熟人工林试验区每个径流场投影面积为10 m×20 m，通过测流堰用自动水位计记录径流量。每个试验区都配有1个自动雨量计用于测定大气降水。测流堰顶角为20°，流量计算公式［14］如下:

式中:Q0为水位自计曲线上某一水位高度的流量，m3/s，H为堰口水头高度，m。Q为某一降水过程的径流总量，Q0和Q'0为水位曲线上相临两点水位高度(H1和H2)的流量，m3/s。T为相临2点的时间差(s)，i=1，2，3，…，n，分别表示某一径流过程 1，2，3，…，n个相临2点的径流时段。最后将径流流量除以样地面积换算为径流深，mm。地表径流为连续监测，记录间隔为5 min 1次;推移质每年年底收集，在沉沙池底部称量泥沙总量，并取样品回实验室在烘箱中105℃干燥至恒重，计算出推移质干重。

2.2 分析方法

所有的统计分析均用SPSS(Version14.0)软件进行。不同林型之间径流量与系数和径流土壤产出的显著性差异均用(One-way ANOVA)单因素方差和(Repeated measure ANOVA)重复测定方差分析，使用LSD检验，当 P＜0.05时视为统计上差异显著。

3 结果与讨论

3.1 地表径流特征

图1显示成熟人工林试验区4种林型和幼龄人工林试验区5种林型的地表径流系数各月间的变化情况。4—9月之间各林型地表径流变化较大，从量上来说要高于10月到翌年3月份，显示出了明显的干湿季差异。从降雨量规律来说，6月份降雨量最大，5月份次之，2月份降雨量最小，1月份次之。其中，成熟人工林和幼龄人工林地表径流均远远小于草坡地表径流。成熟人工林6年月平均地表径流量5月份最大，6月份次之，1月份径流量最小，12月份次之。幼龄人工林月平均地表径流量为6月份最大，5月份次之，2月份径流量最小而1月份次之。各林型地表径流量5月份和6月份呈现峰值，1、2、11、12月处于谷底。这与申卫军等［15］等所得出的径流月度变化规律基本一致。

结合表1和2的降雨量，幼龄人工林实验区6年降雨平均值为1 695.6 mm，变异系数21.2%，最大降雨量2012年为2 217.3 mm，最小降雨量2011年为1 132.4 mm。成熟人工林实验区6年降雨平均值为1 754.7 mm，变异系数20.1%，最大降雨量在2012年，为2 196.8 mm，最小降雨量在2011年为1 194.5 mm。

表1中幼龄人工林试验区不同林型6年平均径流量为157.2 mm，径流量从小到大依次为:尾叶桉林EU1(109.4 mm)＜10物种混交林 M1(129.7 mm)＜厚荚相思林AC(169.1 mm)＜30物种混交林M2(182.8 mm)＜红锥纯林CH(194.8 mm)。变异系数为43.7%～70.1%。表2成熟人工林试验区不同林型6年平均径流量为57.1 mm，径流量从小到大依次为:荷木林SS(41.3 mm)＜窿橼桉林EU(42.0 mm)＜马占相思林AM(63.7 mm)＜松林PM(78.6 mm)。变异系数为26.3% ～33.3%。

表1 2008—2013年幼龄人工林降雨与径流动态Tab．1 Dynamics of rainfall and runoff in young plantation from 2008 to 2013

表2 2008—2013年成熟人工林降雨与径流动态Tab．2 Dynamics of rainfall and runoff in mature plantation from 2008 to 2013

从表1可以看出:2008—2010年期间，幼龄人工林各幼龄林型地表径流量随着时间推移径流系数有一个较为显著的降低过程，而在2010年到2013年的过程中趋于平稳，即使在2012年和2013年雨水相对充沛的年份里的地表径流系数也无显著增加，且各种林型也几乎趋于稳定;因此人工林在6林龄时基本处于拦蓄径流能力快速增加的末端，之后效果较为缓和，从森林类型上看，总体来说混交林径流要少于纯林。而在成熟人工林中(表2)，地表径流系数无显著变化，说明成熟林减流效果已经达到一个平稳较低的水平且荷木混交林相较于其他纯林减流效果最好。

3.2 成熟人工林与幼龄人工林径流差异

林龄和观测年份做两因素重复测量方差分析，表3的结果表明:成熟人工林和幼龄人工林在林龄间差异极为显著(P＜0.01)，说明成熟人工林相对于幼龄人工林更能有效减少径流;年际间的差异为极为显著(P＜0.01)，表明随着时间推移，径流有显著的降低趋势。

国内有研究分析认为森林对地表径流具有良好的调节功能，随着森林植被率的增加，地表径流的形成和土壤侵蚀明显减少。地表径流是引起流域水文变化的主要因子，是洪水流量的主要成分，同时也是水土流失、土壤侵蚀的一个重要因素［16］。

表3 成熟与幼龄人工林在2008年到2013年间的径流系数的重复测量方差分析Tab．3 P value of repeated measurement ANOVAs for young plantation and mature plantation during 2008－2013

从图2a可以看出:幼龄人工林不同林型间地表径流系数差异显著，尾叶桉纯林径流系数要显著小于其他林型，说明在幼龄人工林中，快速生长树种尾叶桉林能较早地起到减少径流的作用［17-18］，且其蒸腾作用较强［19-20］，减少了地表径流的流出量;红锥纯林除了30物种混交林外，径流系数显著大于其他林型，红锥树种在早期生长较为缓慢，属于后发树种。有研究发现其30年后能达到甚至超过桉林的蓄积量［21］，但早期由于个体较小，林冠不发达导致其无法有效截留降雨，导致林区径流量较大。

图2 2008—2013年成熟林年平均径流系数Fig．2 Average annual runoff coefficient of mature plantation from 2008 to 2013

从图2b可以看出，2008—2013年松树林径流系数显著高于其他林型，其次是马占相思林，最小的是蕯橼桉林和荷木林且显著低于其他林型。从结果来看阔叶林减少径流能力要比针叶林强，而同为速生树种的马占相思和蕯橼桉，在成熟林阶段对径流的减流能力有一定的差别别，桉林要好于马占相思林，这在幼龄林中也是呈现这样的趋势，作为乡土树种的荷木林减少径流的能力最强。

通过对不同林型间的产流对比发现，天然林、混交林优于人工林、纯林，复合层结构要好于单层，林龄与水土保持功能呈现正相关关系。这一结论与Narain 等［22］和 Sánchez等［23］等学者的结论一致，通过20多年恢复期的成熟林，年际间的地表径流林型间已经没有明显差异，从类型上看也是混交林地表径流要小于纯林。

3.3 成熟人工林与幼龄人工林土壤推移质特征

从林龄上来看(图3)，幼龄人工林(3～7林龄)中各林型推移质年平均值土壤侵蚀量是成熟人工林(24～28年林龄)的7.14倍。幼林人工林自2008—2012年土壤推移质量各林型平均值是0.5 t/hm2，而成熟人工林平均为0.07 t/hm2。从各自林型来看，成熟人工林中推移质量为:荷木林(SS)＜马占相思林(AM)＜蕯橼桉林(EU);幼龄人工林中表现为:厚荚相思林(AC)＜10物种混交林(M1)＜30物种混交林(M2)＜尾叶桉林(EU1)＜红锥纯林(CH)。

图3 幼龄人工林(a)和成熟人工林(b)推移质量Fig．3 Young plantation(a)and mature plantation(b)soil loss

3.4 成熟人工林与幼龄人工林土壤推移质差异

以林龄和观测年份对年际土壤推移质量做两因素重复测量方差分析，结果表明:成熟人工林和幼龄人工林在林龄间土壤流失量有差异但不显著(P=0.070)，而年际间的差异显著(P=0.034)，说明成熟人工林相对于幼龄人工林虽然在总的推移质量上差异不大，但年际间两者还是存在差异的。究其原因是因为后面几年的年际土壤流失量不同林龄之间差异越来越小，也越来越趋近于0，表明随着时间推移，年际推移质量逐年递减。

表4 成熟人工林与幼龄人工林在2008年到2012年间的推移质的重复测量方差分析Tab．4 P value of repeated measurement ANOVAs with young plantation and mature plantation

从图4可以看出，桉林和红锥林与其他林型土壤推移质量有显著性差异(P＜0.05)，说明在幼龄人工林中，桉林和红锥对土壤侵蚀的保护力比其他林型弱，原因可能是林冠不如其他林型茂密，穿透雨对地面侵蚀较大。从土壤流失和径流系数可以看出:虽然桉林对径流量的减流作用较好，但对土壤的保持能力却较弱;而厚荚相思虽然对减流能力比较弱，但是对土壤的保持能力却最强;红锥纯林相比其他树种对径流和土壤流失的阻滞能力最弱。

从上图(图4a)可以看出:2007—2009年土壤推移质量不同林型差异显著(P＜0.05)，尾叶桉林和红锥纯林最高，显著高于10物种混交林和厚荚相思林，而30物种混交林与其他林型差异不显著;2008—2010年土壤推移质量(图4b)中红锥纯林显著高于10物种混交林和厚荚相思林，而30物种混交林和尾叶桉林与其他林型差异不显著;2009—2011年(图4c)虽然红锥还大于其他林型，尽管差异未达到统计学上的显著水平，但其绝对值还是远高于其他林型。以上说明红锥林对土壤保持效果相对于混交林和桉树林效果较差，而桉林在早期也保持土壤的效果也较差，但随着时间推移，桉林的效果也逐渐得到了加强。

图4 幼林人工林不同时间段的年均推移质量Fig．4 Average annual bed load sediment of young plantation in different time

4 结论

研究期间，鹤山站降雨年内呈明显干湿季划分，年际间差异较大，是造成不同年份径流量和推移质差异的因素之一。研究结果表明:成熟人工林和幼龄人工林地表径流均远远小于草坡地表径流。成熟人工林的地表径流和推移质流失均小于幼龄人工林。其中成熟林地表径流相较于幼龄林减少了64.2%，推移质流失减少了86.0%。幼龄林在3～5年林龄期间，各林型水土流失有一个显著的下降趋势，7年林龄后达到一个稳定较低水平，而成熟林中泥沙流失趋势已经趋于平稳，年际间差异不大。从树种上看，纯林中的桉林在幼龄林阶段能有效减少地表径流，但推移质流失效果显示较慢。长期来看，桉林与荷木混交林均表现出良好的水土保持效果，得出这一结果可能是桉林早期的快速生长使得其蒸腾作用较强，降雨快速以蒸腾的方式回归大气，而后期桉树进入缓慢生长期时，由于其稀疏的冠层使得其他物质得到快速生长，从而径流也较少。而混交林由于不同物种的协同作用，使得径流和推移质侵蚀均较小。结果也表明，在华南地区气候水热条件充足的情况下，虽然幼龄不同林型水土保持能力有差异，只要保护得当，后期均能有不错的保持效果。

［1］ ZHOU Guoyi，WEI Xiaohua，CHEN Xiuzhi，et al．Global pattern for the effect of climate and land cover on water yield［J］．Nat Commun，2015，6(3):5918．

［2］ 李文华，何永涛，杨丽韫．森林对径流影响研究的回顾与展望［J］．自然资源学报，2001(5):3．LI Wenhua，HE Yongtao，YANGLiyun，et al．A summary and perspective of forest vegetation impacts on water yield［J］．Journal of Natural Resources，2001(5):3．

［3］ 石培礼，李文华．森林植被变化对水文过程和径流的影响效应［J］．自然资源学报，2001(5):86．SHI Peili，LI Wenhua．Influence of forest cover change on hydrological process and watershed runoff［J］．Journal of Natural Resources，2001(5):86．

［4］ GARCA-RUIZ JOS M，BEGUERA SANTIAGO，NADAL-ROMERO ESTELA，et al．A meta-analysis of soil erosion rates across the world［J］．Geomorphology，2015，239:16．

［5］ GUO Qiankun，HAO Yanfang，LIU Baoyuan．Rates of soil erosion in China:a study based on runoff plot data［J］．Catena，2015，124:68．

［6］ 郑粉莉，唐克丽，白红英．标准小区和大型坡面径流场径流泥沙监测方法分析［J］．人民黄河，1994(7):19．ZHENG Fenli，TANGKeli，BAI Hongying，et al．Analysis of monitoring methods of runoff and sediment in standard plot and large sloping runoff field［J］．Yellow River，1994(7):19．

［7］ 邓湘雯，康文星，田大伦，等．不同年龄阶段杉木人工林生态系统的径流规律［J］．林业科学，2007(6):1．DENG Xiangwen，KANG Wenxing，TIAN Dalun，et al．Runoff changes in Chinese fir plantations at different age classes，Huitong，Hunan Province［J］．Scientia Silvae Sinicae，2007(6):1．

［8］ 王斌，杨校生，李正才，等．不同气候区4种典型地带性植被土壤水文功能比较［J］．林业资源管理，2009(6):84．WANG Bin， YANG Xiaosheng，LI Zhengcai，et al．Comparison of hydrological functions of four typical zonal vegetation in different climatic zones［J］．Forest Resources Management，2009(6):84．

［9］ 江淼华，谢锦升，王维明．闽北不同土地利用方式水土流失与降雨历时的关系研究［J］．南昌工程学院学报，2011(4):6．JIANG Miaohua，XIE Jinsheng，WANG Weiming，et al．Relationship between rainfall duration and loss of soil and water in diverse land-use types of northern Fujian［J］．Journal of Nanchang Institute of Technology，2011(4):6．

［10］黄承标，刘运华，秦武明，等．3种植被恢复过程中水土流失规律的研究［J］．南京林业大学学报(自然科学版)，2010(2):59．HUANG Chengbiao，LIU Yunhua，QIN Wuming，et al．Study on the water and soil erosion rules during the three types of vegetation recovery process［J］．Journal of Nanjing Forestry University(Natural Sciences Edition)，2010(2):59．

［11］张喜，薛建辉，许效天，等．黔中喀斯特山地不同森林类型的地表径流及影响因素［J］．热带亚热带植物学报，2007(6):527．ZHANG Xi，XUE Jianhui，XU Xiaotian，et al．Forest surface runoff and its influence factors in Karst mountainous area in center of Guizhou Province，China［J］．Journal of Tropical and Subtropical Botany，2007(6):527．

［12］罗兴录，樊吴静，杨鑫，等．不同植被下水土流失研究［J］．中国农学通报，2013(29):162．LUO Xinglu，FAN Wujing，YANG Xin，et al．The study on the water and soil loss in different vegetations［J］．Chinese Agricultural Science Bulletin，2013(29):162．

［13］广东森林编辑委员会．广东森林［M］．北京:中国林业出版社，1990:45．Guangdong Forest Editorial Board．Guangdong forest［M］．Beijing:China Forestry Publishing House，1990．

［14］周国逸．生态系统水热原理及其应用［M］．北京:中国气象出版社，1997:147．ZHOU Guoyi．Ecosystem hydrothermal principle and appplication［M］．Beijing:China Meteorological Press，1997:147．

［15］申卫军，周国逸，彭少麟，等．南亚热带鹤山5种生态系统的地表径流［J］．热带亚热带植物学报，1999(4):273．SHEN Weijun，ZHOU Guoyi，PENGShaolin，et al．Surface runoff in five ecosystems of Heshan subtropical hilly land［J］．Journal of Tropical and Subtropical Botany，1999(4):273．

［16］刘世荣，温远光，王兵等．中国森林生态系统水文生态功能规律［M］．北京:中国林业出版社，1996:201．LIU Shirong，WEN Yuanguang，WANG Bin，et al．Hydrological laws functions laws of major forset ecosystems in China［M］．Beijing:China Forestry Publishing House，1996:201．

［17］ LIU Zhanfeng，WU Jianping，ZHOU Lixia，et al．Effect of understory fern(Dicranopteris dichotoma)removal on substrate utilization patterns of culturable soil bacterial communities in subtropical Eucalyptus plantations［J］．Pedobiologia-International Journal of Soil Biology，2011，55(1):7．

［18］ DUAN Wenjun，REN Hai，FU Shenglei，et al．Community comparison and determinant analysis of understory vegetation in six plantations in South China［J］．Restoration Ecology，2010，18(2):206．

［19］ POORE M E D，FRIESC．The Ecological effects of Eucalyptus［M］．FAO Paper No．59:4．

［20］ BAI J Y，GAN S M．Social，economic and ecological problems of Eucalyptus plantations［J］．Word Forestry Research，1996(2):63．

［21］李松海，蒋军林，吴敏，等．29年生红锥人工林经济效益分析［J］．南方农业学报，2016(6):984．LI Songhai，JIANG Junlin，WU Min，et al．Economic benefit evaluation of 29-year-old Castanopsis hstrix plantation［J］．Journal of Southern Agriculture，2016(6):984．

［22］ NARAIN P，SINGH R K，SINDHWAL N S，et al．Agroforestry for soil and water conservation in the western Himalayan Valley Region of India 2．Crop and tree production［J］．Agroforestry Systems，1997，39(2):175．

［23］ SNCHEZ L A，ATAROFF M，LPEZ R．Soil erosion under different vegetation covers in the Venezuelan Andes［J］．The Environmentalist，2002，22(2):161．