王荣嘉, 高 鹏, 李 成, 刘潘伟, 孙鉴妮, 刘 晴

山东农业大学林学院,国家林业局泰山森林生态站,山东省土壤侵蚀与生态修复重点实验室, 泰安 271018

降雨能够引起产流,随着雨滴的击溅地表与地表径流的产生,造成土壤中氮磷等元素随地表径流流失[1],并产生面源污染。为了解径流产流与氮磷流失特征,减少氮磷流失,防治面源污染,国内外不少专家学者都进行了大量的研究。有些学者针对某一地区的不同土地利用类型的产流及氮磷流失特征进行研究,如Lewis等[2]针对热带湿林的氮磷流失进行研究,周林飞等[3]关于辽宁浑河流域不同土地类型地表径流和壤中流氮磷流失特征的研究,王国重等[4]在豫西南山区对不同土地类型氮磷流失进行研究,吕婷等[5]在长兴县合溪水库集水区根据不同土地利用方式对径流氮流失特征进行研究；也有学者根据不同水土保持措施对产流与氮磷流失影响进行研究,如王萍[6]研究拦截草带对坡耕地水土流失及氮磷输出的控制作用研究,王全九等[7]在黄土坡地研究不同植被覆盖对产流产沙及氮磷流失的影响,Choi等[8]关于植物过滤带(VFS)对氮磷污染的控制作用的研究。研究方法也以野外径流小区观测[9- 11]和人工模拟降雨[12- 14]居多。但目前的研究大多是对径流产流及氮磷流失特征进行研究,研究对象也以地表径流为主[15- 17],关于林地土壤中壤中流和地表径流流失比例及其氮素流失比例等问题尚缺少较为深入地研究。而已有研究发现地表径流与壤中流迁移特征差异较大,随壤中流迁移的氮磷元素对导致的水体富营养化的作用不容忽视[18]。因此,本研究选择鲁中南山区典型区域药乡小流域的麻栎林地及荒草地为研究对象,通过人工模拟降雨试验,研究麻栎林与荒草地的地表径流与壤中流及氮素流失特征,探讨麻栎林地产流及其氮素流失机制,为鲁中南山区径流和氮素流失规律研究及其有效控制与水土保持林的营造提供科学依据。

1 材料与方法

1.1 研究区概况

本研究于泰山森林生态系统国家定位观测研究站(117°04′—117°22′E, 36°17′—36°27′N)进行,该站位于山东省泰安市黄前水库药乡小流域,地处北方土石山区鲁中南山丘区西北部(图1),所属一级流域为黄河流域,支流为大汶河。气候类型属暖温带大陆性亚湿润季风气候,其多年平均气温为18.5 ℃,≥10 ℃的活动积温为3821 ℃,无霜期为196 d,多年平均降水量为727.9 mm,75%的降水集中于6—9月。土壤类型主要为棕壤,主要人工林植被有麻栎(QuercusacutissimaCarr.)、刺槐(RobiniapseudoacaciaL.)、赤松(PinusdensifloraSieb.)、日本落叶松(LarixkaempferiCarr.)等；灌木有胡枝子(LespedezabicolorTurcz.)、酸枣(ZiziphusjujubaMill. var.spinosa(Bunge) Hu)、三裂绣线菊(SpiraeatrilobataL.)以及南蛇藤(CelastrusorbiculatusThunb.)等；草本主要有黄背草(ThemedatriandraForsk. Var.)、野艾蒿(ArtemisialavandulaefoliaDC.)和白羊草(Bothriochloaischaemum(L.) Keng)等。

图1 研究区地理位置图Fig.1 Geographical location of the study area

1.2 试验设计

1.2.1野外取样

本研究是在泰山森林生态系统国家定位观测研究站进行的,该站位于山东省泰安市黄前水库药乡小流域,地处北方土石山区鲁中南山丘区的典型区域。麻栎林是该地区典型水土保持林树种,药乡小流域是代表鲁中南山丘区的典型区域,选择该流域麻栎林进行观测试验,具有一定的典型性和代表性。

于2017年7月选择该区具有代表性的麻栎林作为研究对象,建立60 m×80 m样方,调查林地基本情况。该麻栎林样地树龄为42年,林分密度825株/hm2,郁闭度0.65,样地地表枯落物层的厚度为3—4 cm。地表植被以麻栎幼苗和南蛇藤为主,样地灌木层盖度10%,草本层盖度41%,草本植物以黄背草与野艾蒿为主。该区土层平均厚度40 cm,土壤类型为棕壤,坡度为14°—18°。选择麻栎林与其对照荒草地的典型标准地,结合试验径流槽,每个标准地选取3个地形地貌、植被覆盖等条件相一致的取样样地(共6个样地),每个取样大小为1.5 m×2 m,用于平行降雨。根据该土石山区土层平均厚度约为40 cm,并考虑到标准地的土层厚度的实际情况,取土层深度为40 cm。试验原状土样采集时,参考了钱婧[19]、Mardamootoo[20]以及王国重[4]等人取土方法,尽可能保持地表植被和土层不被破坏,将地表植被、枯落物层和土层进行分开和分层(每10 cm土层取土一次并将土样按由上至下的顺序封装编号)取样。取样后,运至山东农业大学实验基地人工降雨大厅,每个土样的土体按照取样先后顺序放入试验径流槽内并静置堆放一定时间(30 d左右),并再次测定土壤容重等指标,使得径流槽土样的土壤容重接近原状土样的土壤容重值,然后进行人工模拟降雨产流试验。试验标准地基本情况见表1。

表1 麻栎林与荒草地标准地基本情况

粒径分级标准采用美国农部制；TN：全氮,Total nitrogen

1.2.2人工模拟降雨试验

人工模拟降雨试验采用南京南林电子科技有限公司研制的NLJY- 10人工模拟降雨系统,具体参数见表2。根据泰安当地近30年降雨资料中降雨强度水平,考虑到该区暴雨历时短、强度大的特点,并结合试验需要,模拟降雨试验采用一种代表性雨强,即重现期为2.5年的雨强75 mm/h。人工模拟降雨系统下部放有试验径流槽,内有供试土壤。考虑到该研究区标准地坡度为14°—18°,因此,径流槽倾斜坡度确定为16°,与标准地平均坡度一致且具有代表性。试验径流槽可同时收集地表径流和壤中流(图2)。为了保证试验准确性,模拟降雨时6个样地进行平行降雨实验,共进行6次降雨实验(麻栎林和荒草地各3次)。为保证降雨强度均匀稳定,降雨前,对降雨强度进行2次率定,2次率定误差不超过5%,于降雨结束后再次率定,取前后率定结果的平均值作为该试验的降雨强度,每次降雨历时40 min。降雨开始后,记录降雨开始时间、地表径流以及壤中流产流时间,径流产流后,每间隔5 min采集地表径流和壤中流各一次,并记录采样径流量,用于计算不同时刻的径流量。记录降雨停止时间、径流停止时间,采样直至降雨径流结束。

表2 人工模拟降雨系统参数表

图2 试验装置示意图Fig.2 Schematic diagram of experiment

1.3 样品分析方法及数据分析

水样采集后加浓硫酸至pH值为1左右进行保存并带到山东省土壤侵蚀与生态修复重点实验室进行实验。水样全氮(TN)测定采用过硫酸钾氧化-紫外分光光度法[21],水样中的全氮(TN)包括铵态氮、硝态氮、亚硝态氮和能被转换为铵盐的有机氮,本研究所采用的方法是将水样在过硫酸钾存在条件下,通过高压消煮,将亚硝酸根、铵、有机态氮均转换成硝酸根,通过计算出硝酸根的浓度从而计算出全氮(TN)的浓度；由于试验用水一致,故不考虑试验用水对结果的影响。土壤容重和土壤孔隙度测定采用环刀法。

1.4 数据处理与统计分析

文中数据、图表采用Excel 2007处理分析。

2 结果与讨论

2.1 麻栎林和荒草地产流过程特征

由表3和图3可以看出,麻栎林与荒草地的地表径流均产生较早,麻栎林在3 min时开始产流,荒草地在2 min时开始产流,麻栎林的地表径流产流时间滞后于荒草地的,这是由于麻栎林地表植被覆盖度较高,起到拦截地表径流、延缓地表径流产流的作用。二者的地表径流量均呈现不断增加且趋于稳定的特征,这是由于在降雨初期,土壤含水量较低且未达到饱和,到达地表的降雨先进行土壤入渗,用于填充土壤孔隙,故初期的地表径流产流较少；随着降雨的不断进行,到达地表的水分也越来越多,土壤含水量逐渐增加达到饱和状态,并开始产生壤中流,降雨雨水下渗速度也开始接近于饱和土壤的稳定下渗速度,地表径流量到达最大值并保持稳定。陈晓冰等[22]在竹林地的研究也发现,降雨径流开始产生时径流量相对较小,随着降雨时间的持续,地表径流量逐渐增大,随着时间推移,地表径流量趋于稳定,这与本研究结果相一致。根据试验结果,麻栎林的地表径流量为9.7 L,荒草地的为15.8 L,麻栎林地表径流量是荒草地的61.4%,这是由于麻栎林的地表植被覆盖大于荒草地的,且有一层枯落物层,拦截了一部分到达地表的雨水,促进雨水下渗,对地表径流的有一定拦截作用,减少了地表径流量。陈安磊等[23]在桃源县的通过野外径流小区研究发现,自然林的径流量仅为草地区径流量的25.3%,自然林的径流量小于草地区的；吕婷等人[5]研究发现,在自然降雨条件下,林地(人工林)的径流量(360.3 m3/hm2)小于荒地的径流量(900.9 m3/hm2),这都与本研究结果一致。 Hosseini等[24]研究发现,植被覆盖对减少径流的起到了重要作用,而本研究中,麻栎林的植被覆盖大于荒草地的,其地表径流量也大于荒草地的,与Hosseini等[24]研究结果相一致。

图3 麻栎林和荒草地地表径流与壤中流产流特征Fig.3 Characteristics of surface flow and interflow inQ.acutissimaforest and grassland

样地 Sample area径流类型Runoff type样地1 Sample area 1样地2 Sample area 2降雨量Rainfall/mm产流时间Runoff time/min径流量 Runoff volume/L占比Ratio/%降雨量Rainfall/mm产流时间Runoff time/min径流量 Runoff volume/L占比Ratio/%麻栎林地地表径流-39.261.3-39.862.4Q. acutissima壤中流-85.837.7-75.937.6总产流50-15.0-50-15.7-荒草地Grassland地表径流-215.582.4-215.680.0壤中流-103.317.6-113.920.0总产流50-18.8-50-19.5-样地 Sample area径流类型Runoff type样地3 Sample area 3样地平均值 Sample average降雨量Rainfall/mm产流时间Runoff time/min径流量 Runoff volume/L占比Ratio/%降雨量Rainfall/mm产流时间Runoff time/min径流量 Runoff volume/L占比Ratio/%麻栎林地地表径流-210.161.6-39.761.8Q. acutissima壤中流-76.338.4-7637.2总产流50-16.4-50-15.7-荒草地Grassland地表径流-216.280.2-215.881.1壤中流-94.019.8-103.719.9总产流50-20.2-50-19.5-

由表3可以看出,麻栎林与荒草地的壤中流的产流时间分别为7 min和10 min时,均滞后于地表径流的,这是由于当土壤水分含量达到饱和才开始产生壤中流。由图3可以看出,麻栎林和荒草地的壤中流量波动比较小,产流过程保持相对稳定,这是因为壤中流产流是由于土壤中水分达到饱和时产生,下渗速度已基本稳定。本研究表明,麻栎林地的土壤孔隙度大于荒草地的,而麻栎林地的土壤容重则小于荒草地的(表1),这是因为与荒草地相比,麻栎林地中的地表植被较多,枯枝落叶物较多,加之地下部分树木根系的作用,较好地改良了土壤的结构,使得麻栎林地土壤孔隙度增加,土壤容重降低,有利于地表径流下渗,促进壤中流的形成与发展,本文中,麻栎林壤中流量为6 L,是荒草地壤中流量的162.2%。沈紫燕[25]通过模拟降雨试验研究发现,在相同土壤条件下,容重为1.05 g/cm3的土壤首先产流,容重为1.15、1.20 g/cm3的土壤相继产流,并且平均壤中流强度的大小和壤中流总量的大小关系为容重1.05>1.15>1.20,土壤容重越大,孔隙率越小,土壤越密实,越不利于壤中流的产生,容重越大壤中流越难产生。这与本研究中土壤容重越小越有利于壤中流产流的结论相一致。在杜祥运等[26]研究中,在同一雨强下,壤中流量大小排序依次为草灌>草本>灌木>裸地,地表覆盖度越高,壤中流量也越大,而本研究中,麻栎林的地表覆盖度大于荒草地的,壤中流量也大于荒草地的,与杜祥运等[26]研究相一致。

由表3还可以看出,麻栎林的总产流量小于于荒草地的,是荒草地的80.5%。麻栎林地表径流量占总产流量的61.8%,荒草地地表径流量占总产流量的81.1%,说明麻栎林和荒草地的产流都以地表径流为主。麻栎林和荒草地的壤中流量占总产流量比例分别是37.2%和18.9%,麻栎林的大于荒草地的,说明麻栎林与荒草地相比,可以促进降雨下渗,增加壤中流,减少地表径流的流失。由图3还可看出,壤中流的产流历时长于地表径流的,这是由于在降雨停止初期,地表径流产流停止,由于地表积水及土壤水分含量较大,壤中流依然保持着稳定的状态,随着降雨停止且缺乏后续降雨补充,壤中流的产流量逐渐减少。

通过对麻栎林和荒草地地表径流和壤中流随时间变化过程进行回归分析,得到关于地表径流量和壤中流量与降雨时间的函数关系方程,可用来模拟径流随时间的变化过程,其中地表径流量与降雨时间呈较好的对数函数关系,这与唐柄哲[27]研究结果相一致。

麻栎林地表径流量：y=839.63ln(x)+321.15,R2=0.781

荒草地地表径流量：y=1006.1ln(x)+625.45,R2=0.7071

麻栎林壤中流量：y=0.0514x3- 2.0952x2+22.968x+47.461,R2=0.9565

荒草地壤中流量：y=0.0496x3- 2.0126x2+22.659x- 9.0999,R2=0.9288

式中,x为降雨时间(min),y为径流量(mL/min)。

2.2 麻栎林和荒草地的径流氮素流失及其动态变化特征

2.2.1地表径流全氮迁移特征

图4 麻栎林与荒草地地表径流与壤中流全氮流失特征Fig.4 Characteristics of surface flow and interflow and nitrogen loss inQ.acutissimaforest and grassland

样地Sample area径流类型Runoff type样地1 Sample area 1样地2 Sample area 2径流量Runoff volume/L全氮流失量TN loss/mg占比 Ratio/%全氮平均浓度 TN average concentration/(mg/L)径流量Runoff volume/L全氮流失量TN loss/mg占比 Ratio/%全氮平均浓度 TN average concentration/(mg/L)麻栎林地地表径流9.2127.470.713.89.8127.270.413.0Q. acutissima壤中流5.852.929.39.15.953.529.69.1总产流15.0180.3-12.015.7180.7-11.5荒草地 Grassland地表径流15.5240.687.115.515.6240.187.115.4壤中流3.335.512.910.83.935.612.99.1总产流18.8276.1-14.719.5275.7-14.1样地Sample area径流类型Runoff type样地3 Sample area 3样地平均值Sample average径流量Runoff volume/L全氮流失量TN loss/mg占比 Ratio/%全氮平均浓度 TN average concentration/(mg/L)径流量Runoff volume/L全氮流失量TN loss/mg占比 Ratio/%全氮平均浓度 TN average concentration/(mg/L)麻栎林地地表径流10.1126.770.212.59.7127.170.413.1Q. acutissima壤中流6.353.829.88.5653.429.68.9总产流16.4180.5-11.015.7180.5-11.5荒草地Grassland地表径流16.2239.986.814.815.8240.287.015.2壤中流4.036.613.29.23.735.913.09.7总产流20.2276.5-13.719.5276.1-14.2

通过对麻栎林和荒草地地表径流全氮流失浓度随时间变化过程进行回归分析,发现地表径流全氮流失浓度与降雨时间呈较好的幂函数关系,即随时间变化的地表径流全氮流失浓度可以通过幂函数来模拟,这与王茹[30]研究结果相一致。

麻栎林：y=17.884x-0.165,R2=0.9221

荒草地：y=18.89x-0.108,R2=0.982

式中,x为降雨时间(min),y为全氮流失浓度(mg/L)。

2.2.2壤中流全氮迁移特征

由图4和表4看出,在整个壤中流产流过程中,麻栎林和荒草地的壤中流全氮输出浓度一直处于一个相对稳定的状态,荒草地壤中流全氮输出浓度在10 mg/L上下波动,麻栎林的在7—10 mg/L之间波动,波动起伏变化比较小。造成这种现象的原因主要是因为麻栎林和荒草地的壤中流流速比较缓慢,流速波动较小,使雨水在下渗过程中所产生的壤中流与土壤中的氮元素交换作用达到一种平衡,氮元素的流失量相对比较稳定。周林飞等[3]在辽宁浑河流域研究发现,耕地与草地壤中流中总氮的浓度在整个径流过程中呈波浪状输出,总体趋势相对稳定,无较大的起伏变化,与本研究结果相似。麻栎林的壤中流全氮输出浓度平均值为8.9 mg/L,而荒草地的为9.7 mg/L,麻栎林的比荒草地的低8.2%,且在整个产流过程中,荒草地的壤中流全氮输出浓度始终高于麻栎林的,这是由于麻栎林土壤中植物根系和枯落物较多,可以有效吸收拦截随壤中流迁移的氮元素,减少麻栎林中氮元素的流失浓度。常龙飞等[31]在巢湖低山丘陵地区对6种主要土地利用类型(弃耕地、尾矿裸地、灌木林地、荒草地、马尾松林地和人工恢复林地) 进行研究,壤中流全氮流失浓度按大小排序依次为尾矿裸地>荒草地>次生马尾松林地>灌木林地>弃耕地>人工恢复林地,发现荒草地的壤中流全氮流失浓度大于马尾松林地和人工恢复林地的,这与本研究中荒草地的壤中流全氮输出浓度高于麻栎林的的结论相类似。

通过对麻栎林和荒草地壤中流全氮流失浓度随时间变化过程进行回归分析,得到关于壤中流全氮流失浓度与降雨时间的函数关系方程,可用来模拟壤中流全氮流失浓度随时间的变化过程：

麻栎林：y=0.0002x3-0.0136x2+0.2939x+7.1068,R2=0.9262

荒草地：y=-0.0002x4+0.0096x3-0.1746x2+1.118x+8.0246,R2=0.746

式中：x为降雨时间(min),y为全氮流失浓度(mg/L)。

2.2.3地表径流和壤中流全氮迁移比例分析

图5 麻栎林与荒草地全氮流失比例 Fig.5 The ratio of total nitrogen loss inQ.acutissimaforest and grassland

由表4和图5看出,麻栎林和荒草地的地表径流全氮输出浓度均高于壤中流的,并且麻栎林的全氮总流失浓度比荒草地的低19.0%,其中,麻栎林的地表径流全氮流失量占全氮总流失量的70.4%,壤中流占29.6%；荒草地的地表径流全氮流失量占全氮总流失量的87.0%,壤中流占13.0%；荒草地径流全氮流失量是麻栎林的1.9倍。说明与壤中流相比,荒草地的地表径流迁移而造成的全氮流失是主要的流失方式。同时,地表径流全氮输出浓度始终高于壤中流的,也说明土壤垂直方向上,全氮的流失浓度由上到下减少,说明全氮的流失受土壤虑减机制的影响,减少通过下渗进入地下水,减少全氮的流失量,避免地下水被氮元素污染。可见,通过植被恢复,促进地表径流下渗,减少地表径流量,能够有效减少氮素的流失。

3 结论

(1)麻栎林的总产流量、地表径流量、壤中流量分别是荒草地的80.5%、61.4%、162.2%,麻栎林的总产流量与地表径流量均小于荒草地,但壤中流量高于荒草地的。地表径流产流过程呈现不断增加且趋于稳定的特征；壤中流产流时间滞后于地表径流的,产流过程流量波动比较小,保持相对稳定。地表径流量和壤中流量随降雨时间的变化过程可分别通过对数函数和多项式函数进行模拟。

(2)麻栎林的全氮总流失浓度、地表径流全氮流失浓度、壤中流全氮流失浓度为11.5、13.1、8.9 mg/L,分别比荒草地低19.0%、13.8%、8.2%。地表径流全氮流失浓度一般前期较大,而后递减,至后期趋于稳定；壤中流全氮流失浓度在整个产流过程中保持相对稳定。地表径流和壤中流全氮流失浓度随降雨时间的变化过程可分别通过幂函数和多项式函数进行模拟。

(3)在整个降雨产流过程中,麻栎地表径流量占总产流量的61.8%,荒草地地表径流量占总产流量的81.1%,麻栎地表径流全氮流失量占总流失量的70.4%,荒草地地表径流全氮流失量占总流失量的87.0%,径流、氮素的流失都以地表径流为主,土壤滤减机制明显。可见,通过植被恢复,促进地表径流下渗,减少地表径流的流失量,能够有效减少氮素的流失。与荒草地相比,麻栎林具有明显增加壤中流、减少地表径流和氮素流失效果。