张 骞, 张中华, 马 丽, 杨晓渊, 郭 婧, 徐文华, 周秉荣,邵新庆, 王 芳, 贾永忠, 贺有龙, 金 欣, 周华坤*
(1. 中国科学院西北高原生物研究所青海省寒区恢复生态学重点实验室, 青海 西宁 810001; 2. 青海省气象科学研究所,青海 西宁 810001; 3. 中国农业大学动物科技学院, 北京 100193; 4. 循化县积石镇政府农村经济综合服务中心,青海 海东 810600; 5. 青海省果洛州草原站, 青海 大武 814000; 6. 青海省气象灾害防御技术中心, 青海 西宁 810001)
三江源区位于青藏高原腹地,是长江、黄河、澜沧江的源头,其面积的70%以上为高寒草甸,具有重要的生态固持功能,是不可替代的生态屏障[1-3]。但是,由于自身生态环境较为敏感、脆弱,加之近年来气候变化、过度放牧等的影响,该区高寒草甸发生了极为严重的退化现象[4]。高寒草甸生态系统是在寒冷而湿润的气候条件下,由耐寒的多年生中生草本植物为建群种而形成的一种生态系统类型,广泛分布于青藏高原,是高原畜牧业生产的物质基础[5]。高寒草甸宝贵的牧草资源不仅对当地畜牧业的发展、民族经济发展具有不可代替的作用,也对高寒地区风沙、水蚀、盐碱、内涝、地下水位变化、土壤侵蚀、土地资源的流失等自然灾害的发生具有一定的屏障作用[6]。同时,牧草作为畜牧业发展的基础,对于草甸最大载畜率及生态承载力具有决定性的作用[7]。而高寒草甸退化导致牧草覆盖度及品质的下降,使得三江源区的畜牧业和民族经济发展受到巨大阻碍[8]。土壤在高寒草甸涵养水源、养分供给等功能中发挥着重要作用,高寒草甸的退化导致土壤生态功能下降,影响草甸生态系统的可持续发展[9-11]。因此,退化高寒草甸生态系统的恢复对于改善当前现状具有重大意义。
目前,对于三江源区高寒草甸的退化,学者们提出了一系列的恢复措施,如人工草地建植[12]、施肥[13]、围栏封育等,其中人工草地建植是恢复措施中效果最直观明显的[14]。但在人工草地建植过程中也存在许多的问题,如草种的选择、搭配。同时,建植后人工草地的土壤以及生物量状况缺乏研究。因此,为重度退化高寒草甸恢复提供适宜的牧草组合,以及明晰重建后草甸的土壤及生物量状况是十分必要的。垂穗披碱草(Elymusnutans)、早熟禾(Poapratensis)为青藏高原的本土禾本科植物,而且具有耐寒、利用价值高等优点,常被广泛作为建植种用于退化高寒草甸的恢复中[15-16]。苜蓿作为一种豆科植物,对于禾本科植物种子萌发生长具有促进作用[17]。有研究表明合理的豆禾比例与种类搭配能显著改善土壤养分的供给[18]。故本研究采用人工草地建植的方式,选取垂穗披碱草、草地早熟禾、呼伦贝尔苜蓿(Medicagofalcata)作为研究草种,通过分析不同牧草混播对退化高寒草甸土壤养分及生物量的影响,筛选对于土壤改良以及提高生物量有优势的牧草组合,以期为重度退化高寒草甸人工草地建植提供技术支撑。同时,也为退化高寒草甸治理模式的构建及高寒区生态恢复[19]提供理论依据。
原位试验在三江源高寒草地研究观测站军牧场试验点(青海省果洛州玛沁县军牧场,经纬度34°22′~34°20′ N,100°30′~100°29′ E,海拔4 100 m)进行。该区为典型高寒草甸,土壤类型为高山草甸土。2010-2018年平均年降水量为514.9 mm,生长期(5-8月)降水量上升12.0%。年平均气温为0.7℃,1月为2.1℃,7月为19.8℃,年日照小时为2 493.6 h。属高原大陆性气候,雨热同季。
本试验于2016年1月-2018年12月在一片相对平坦的重度退化的高寒草甸进行。试验草甸通过围栏进行保护,以排除放牧干扰,同时,协调当地牧民进行鼠害防治,尽可能的避免鼠兔等啮齿动物的干扰。采用随机区组试验设计,选择由青海省牧草良种繁殖场提供的高禾草垂穗披碱草、矮禾草草地早熟禾、豆科牧草呼伦贝尔苜蓿进行配置,共8个组合(表1),3次重复,每个处理小区为3 m×3 m,各小区以及重复之间距离均为1 m。2016年5月中旬进行播种,牧草播种前对试验样地进行翻耕处理,翻耕方式为农具翻耕,深度为20~22 cm。选取优质草种均匀混合进行撒播,当年苗期进行2次除杂,之后每年返青后进行一次除杂,除杂方式为人工除杂,剔除试验草种以外的杂草。对照组不作任何处理。播种量以青海建植人工草地的地方标准执行(垂穗披碱草播种量为3 g·m-2。草地早熟禾播种量为0.75 g·m-2。呼伦贝尔苜蓿播种量为1.5 g·m-2)。
表1 牧草混播处理及播种量Table 1 Pasture mixed treatments and seeded rates
注:对照组的处理方式为维持原始的植被及土壤状态,未做任何处理
Note:The control group was treated to maintain the original vegetation and soil status without any treatment
2018年9月中旬进行样品采集。试验样地处理组牧草长势良好,呼伦贝尔苜蓿存在少数死亡情况,但整体存活率较高。在每个处理小区选取能够代表小区牧草整体长势水平的0.5 m×0.5 m样方,刈割样方内地上部分牧草,装入大信封中,带回实验室测定地上生物量。在每个处理小区对角线上均匀选取2个点作为土壤取样点,使用内径5 cm的土钻分别取0~10 cm,10~20 cm,20~30 cm的土壤样品,并将每个处理的2个取样点同层样品混合,装入土样袋中,带回实验室测定土壤养分指标及地下生物量。
土壤养分指标:土壤全氮、全磷、铵态氮、硝态氮、速效磷采用全自动间断化学分析仪(CleverChem 380),土壤速效氮采用碱解扩散法,土壤有机碳采用油浴加热重铬酸钾氧化-容量法。土壤养分指标分别测定0~10 cm,10~20 cm,20~30 cm。本文使用三层土壤养分指标数据的加权平均值。
生物量:将采集的地上部分牧草,105℃杀青30 min,在65℃烘干48 h至恒重后称干重,作为地上生物量。将采集的土壤样品过筛后残留根系,105℃杀青30 min,在65℃烘干48 h至恒重后称干重,作为地下生物量。
采用2013 Excel对数据进行前期处理、隶属函数值的计算,SPSS 20对数据进行因子分析、相关性分析、显著性分析,Origin 2017进行作图。隶属函数值的计算方法:
式中:Zij表示土壤养分状况隶属函数平均值,Xij表示i处理j指标的测定值,Xjmin表示所有处理中j指标的最小值,Xjmax表示所有处理中j指标的最大值。
由表2可知,通过不同牧草混播,退化高寒草甸土壤有机碳含量与全氮含量呈极显著正相关(P<0.01),与速效氮和硝态氮含量呈显著正相关(P<0.05)。土壤全氮含量与速效氮含量呈极显著正相关(P<0.01),与硝态氮含量呈显著正相关(P<0.05)。土壤硝态氮含量与速效磷含量呈极显著正相关(P<0.01)。反映了氮含量显著影响着氮转换,土壤有机碳及速效磷含量与氮循环密切相关。
从表3主因子分析可以看出,3个主因子代表了7个土壤养分指标的94.150%的主要信息。故可选用这3个主因子来作为试验处理结果的分析基础。从表4主因子分析结果可以看出,主因子1的AP,NO 2个指标代表了总指标60.927%的信息。主因子 2的NH 1个指标代表了总指标18.199%的信息。主因子3的TP 1个指标代表了总指标15.024%的信息。因此,将AP,NO,NH,TP这4个土壤指标作为试验处理对土壤养分影响的主要分析指标。
表2 土壤养分指标相关性分析Table 2 Correlation analysis of soil nutrient indicators
注:**在 0.01 水平(双侧)上显著相关,*在 0.05 水平(双侧)上显著相关。表中字母SOC代表“有机碳”,TN代表“全氮”,TP代表“全磷”,AN代表“速效氮”,AP代表“速效磷”,NO代表“硝态氮”,NH代表“铵态氮”,下同
Note:**Significantly correlated at the 0.01 level (both sides),*significant correlation at the 0.05 level (both sides). The letters SOC in the table represent “Soil organic carbon”,TN stands for “Total nitrogen”,TP stands for “Total phosphorus”,AN stands for “Available nitrogen”,AP stands for “Available phosphorus”,NO stands for “Nitrate nitrogen”,and NH stands for “Ammonium nitrogen”,the same as below
表3 土壤养分指标主因子分析Table 3 Main factor analysis of soil nutrient indicators
注:类别中的编号代表主因子分析所提取的公因子,下同
Note:The number in the category represents the common factor extracted by the principal factor analysis,the same as below
表4 主因子分析结果Table 4 Main factor analysis results
不同牧草混播对退化高寒草甸土壤养分的影响显著。由图1A可知,EM,EP,PM处理组的土壤速效磷含量均高于CK(P<0.05),其中EM处理的土壤速效磷含量最高,且均高于其他各处理(P<0.05)。EPM,E,P,M处理组的土壤速效磷含量均高于CK(P<0.05),但处理间差异不显著。由图1B可知,各处理的土壤铵态氮含量波动较大,但均高于CK。其中M处理的土壤铵态氮含量最高,与CK相比提高了约1.2倍(P<0.05)。EPM,EM,P处理组间的土壤铵态氮含量无显著差异(P<0.05),但EPM处理的土壤铵态氮含量高于EM,P处理组。
由图1C可知,EP,EM,PM处理组的土壤硝态氮含量均高于CK,其中EP,EM处理组显著高于其他各处理(P<0.05)。EPM,E,P,M处理组的土壤硝态氮含量均低于CK(P<0.05),其中P处理的土壤硝态氮含量最低,与各处理差异显著(P<0.05)。由图1D可知,E处理的土壤全磷含量最高,且显著高于CK(P<0.05),其余各处理与CK相比均未产生显著差异。
综合图1来看,E处理的土壤全磷含量高于CK,但土壤速效磷含量低于CK(P<0.05)。M处理的土壤铵态氮含量高于CK,但土壤硝态氮含量低于CK(P<0.05)。
图1 各处理下土壤速效磷含量(A)、铵态氮含量(B)、硝态氮含量(C)、全磷含量(D)Fig.1 Soil available phosphorus content (A) and ammonium nitrogen content (B) and nitrate nitrogen content (C) and totalphosphorus content (D) under each treatment注:图中值为平均值+标准误,不同字母表示同一测定指标在P<0.05水平上差异显著。图中字母EPM代表“垂穗披碱草+草地早熟禾+呼伦贝尔苜蓿”,EP代表“垂穗披碱草+草地早熟禾”,EM代表“垂穗披碱草+呼伦贝尔苜蓿”,PM代表“草地早熟禾+呼伦贝尔苜蓿”,CK代表“对照组”,E代表“垂穗披碱草”,P代表“草地早熟禾”,M代表“呼伦贝尔苜蓿”,下同Note:The values in the figure are the mean + standard error,and the different letters indicate that the same measurement index is significantly different at the 0.05 level.The letter EPM stands for "Elymus nutans+Poa pratensis+Medicago falcata",EP stands for "Elymus nutans+Poa pratensis",EM stands for "Elymus nutans+Medicago falcata",PM stands for "Poa pratensis+Medicago falcata",CK stands for "control group",E stands for "Elymus nutans",P stands for "Poa pratensis",and M stands for "Medicago falcata",the same as below
选取不同牧草混播条件下退化高寒草甸土壤的AP,NO,NH,TP 4个主因子养分指标,利用隶属函数法对各指标的相对值进行计算,为不同牧草组合播种处理后土壤养分状况做出综合评价(表5),得到经不同牧草混播处理后土壤养分状况的排序为EM>EP>PM>E>CK>EPM>M>P。
表5 各处理下主效指标的综合评价值及排序Table 5 Comprehensive evaluation value and ranking of main effect indicators under each treatment
通过不同牧草混播处理,重度退化草甸的生物量发生了显著的变化。由图2A可知,各处理组的地上生物量均显著高于CK(P<0.05)。其中EM处理的地上生物量最高,较CK提高了约2.3倍。M处理的地上生物量最低,较CK提高了约1.5倍。如图2B可知,各处理组的地下生物量均高于CK(P<0.05),其中EM处理的地下生物量最高,较CK提高了约2.7倍。M处理的地下生物量最低,较CK提高了约1.2倍。由此可见,地上生物量与地下生物量的各处理组具有同步性,EM处理的效果最为显著。
图2 各处理下地上生物量(A)、地下生物量(B)Fig.2 Aboveground biomass (A) and Underground biomass (B) under each treatment
土壤养分状况的改善是退化草甸恢复的目标之一,它决定着草甸利用是否可持续[20]。养分条件的改变,将影响与之相关的土壤生态功能稳定性[21]。本研究中,不同牧草混播条件下退化高寒草甸土壤有机碳含量与全氮含量呈极显著正相关(P<0.01),与硝态氮含量呈显著正相关(P<0.05)。这是由于土壤氮素和有机碳之间的耦合效应,一方面土壤氮素含量的变化改变了土壤对CO2的吸收能力,影响着土壤碳的固定,另一方面土壤有机碳含量又对氮素的矿化、固氮及反硝化有重要作用[22-23]。不同牧草混播条件下退化高寒草甸土壤全氮含量与速效氮含量呈极显著正相关(P<0.01),与硝态氮含量呈显著正相关(P<0.05)。这是由于土壤全氮含量处于动态变化之中,它在土壤微生物的作用下影响着土壤净氮矿化速率及氮素有效性[24]。
从单个养分指标分析结果来看,土壤速效磷含量在EM,EP,PM处理下均高于CK,其中EM处理的土壤速效磷含量最高。这可能是由于混播导致凋落物量增加,而凋落物是形成腐殖质的重要原料,它会在铁、铝氧化物等胶体表面形成保护膜,减少对磷酸根的吸附,增加了土壤中速效磷的含量[25-26]。土壤全磷含量只有在E处理下高于CK(P<0.05),而其他各个处理均与CK无显著差异。这可能是由于垂穗披碱草的凋落物积累、根系分泌物使得土壤全磷产生了一定程度的富集,而其他各牧草混播的富集效果不明显[27]。土壤铵态氮和硝态氮主要来源于土壤有机氮的氨化和硝化,是植物能直接吸收利用的生物有效态氮[28]。本研究通过分析得出,土壤铵态氮含量在各牧草混播条件下的波动较大,但均高于CK。其中M处理的土壤铵态氮含量最高(P<0.05)。这可能是由于植物生长期间产生的可矿化氮,促进了氮素的转换过程。同时,由于果洛试验区的土壤大多受氮素转换的限制,而豆科植物的固氮作用可使土壤铵态氮含量提升,故呼伦贝尔苜蓿处理下的铵态氮含量最高[29-30]。土壤硝态氮含量在各牧草混播条件下既有升高又有减低,其中2种牧草混播条件下的土壤硝态氮含量均高于CK,而1种牧草和3种牧草混播条件下的土壤硝态氮含量均低于CK。这说明混播可以促进氮素的硝化作用,但适宜的牧草混播数为2种,其中豆禾混播效果更佳[31]。
从总体土壤养分状况的评价结果来看,2种牧草混播下的土壤养分状况均高于CK,而1种牧草和3种牧草混播条件下的土壤养分状况均低于CK。其中,EM处理下的土壤养分状况最佳。这可能是豆禾牧草混播下,不同种类牧草的生态位互补效应,能更有效地利用有限的资源,同时又有利于调节土壤环境条件,从而提高土壤养分状况[18,32]。而M处理的土壤养分状况最差,这是由于豆科单播对于土壤养分的利用较单一,不利于退化高寒草甸土壤养分状况的改善[33]。在豆禾牧草混播体系中豆科牧草将土壤中氮素固定后绝大部分储藏在自身体内,部分转移给禾本科植物,这增强了共生效应和相容性,有利于土壤稳定性及养分状况的提升[29,34]。
生物量的增加是恢复直观的判断标准,它决定着草地的承载力[35-36]。高寒地区的适宜牧草混播对草地生物量的积累和分配有较大的影响,这将有助于理解未来植被恢复和植物间的相互作用[37]。合理的牧草配置能够充分利用各种资源,促进种间有益作用,提高单位面积牧草产量有利于改善牧草品质[38]。本研究得出,不同牧草混播后,地上、地下生物量均高于CK,且多种牧草混播的地上、地下生物量普遍高于1种牧草的播种,其中EM处理的地上、地下生物量均最高。导致这一结果的原因可能为重度退化草甸自身的生产力低下,选取的恢复草种对高寒环境有较好的适应能力,其中垂穗披碱草为本土植物,其适应性在三江源大多区域都极强,呼伦贝尔苜蓿其自身的抗寒性及高寒地区适应性在青海省牧草良种繁殖场培育过程中得到了证明,本研究中其生物量高于对照组也说明了其适应于高寒地区的生长,这2种牧草的应用为重度退化高寒草甸的生态修复提供了优质的种子源[39],从而使得各牧草混播条件下的生物量均高于CK。EPM,EP,PM处理的地上、地下生物量均高于对照组,这可能是由于多种牧草的混播会使群落稳定性增加,初级生产力得到提升[40]。而1种牧草的播种虽然也高于对照组,但这可能是由于其自身的适应性导致。此外,豆禾牧草的混播由于其资源利用的分化,使得在退化草甸恢复过程中占有很大的优势[41],因而EM处理对于提升退化草甸生物量的效果最佳。
从本研究的研究结果来看,采用人工草地建植方式对三江源区重度退化高寒草甸进行恢复时,建议选择EM(垂穗披碱草+呼伦贝尔苜蓿)进行混播,这样不仅可以保证退化高寒草甸的土壤养分状况得到改善,同时也可为当地畜牧业提供优质的牧草资源。在保证草地资源的可持续利用的前提下,维持草地生态系统的稳定性。