基于不同地下水位环境分解过程探究碳、氮、磷释放和中木质素降解的关联性植物分解是有机养分矿化的关键环节,是维持生态系统能量流动、物质循环和养分平衡的重要生态过程[1]。湿地植物的分解会影响湿地生态系统的养分循环周期和湿地土壤肥力的可持续供给,进而通过影响植物的生长速率、发育水平、物种组成和分布格局等来影响植物种群间的关系和群落的构建[2,3]。由于木质素在植物中的生物学作用主要是提高细胞壁的完整性和抵御病原体的攻击[4],通常被认为是阻碍有机物进行生物过程分解的难分解物质[5],其保护下的水溶性成分等易分解组分是制约分解速率的主要因素[6],进而控制着植物分解过程[7],因此,木质素的分解速率在很大程度上可以调节生态系统的碳通量变化和养分循环过程[8]。最近的研究则发现,与陆地生态系统相比,湿地系统具有更强的木质素可分解性,在0—60天内木质素分解速率即可达到最大值[9,10]。然而鲜有研究关注湿地植物分解过程中,木质素的降解对碳、氮、磷等营养元素释放的影响。鄱阳湖是中国最大的淡水湖泊,同时也是世界上重要的季节性洪泛平原湿地[11]。水位波动变化是鄱阳湖湿地生态系统结构和功能的重要调控因子[12],不仅影响了湿地植物的生长和空间分布[13],还会改变湿地植物分解、养分循环等生态过程,进而对湖泊水体环境产生深刻作用[14,15]。鄱阳湖受长江和“五河”的双重影响,10月份后逐渐退水,洲滩出露[16],为以灰化薹草为代表的典型优势植物提供了适宜的生存环境,经过2个多月的秋季生长过程,薹草一般在12月开始枯萎,进入分解状态[17],为鄱阳湖洲滩湿地提供了稳定的营养来源。作者前期对鄱阳湖湿地分解过程与地下水位变化的关系及其环境效应展开了系列研究[18,19,20],而对分解过程中大分子物质降解和元素的释放关系研究尚有不足。此外,湿地洲滩植物衰亡凋落后的分解过程释放的营养元素会通过地表径流和地下水输送到水体,对湖泊水体产生明显的污染效应[20],本研究的开展能够为湖区营养盐管控和水体富营养化监视监测提供重要的参考。本文以鄱阳湖典型碟形子湖白沙湖洲滩湿地优势植物灰化薹草的分解过程为研究对象,分析了木质素分解对碳、氮、磷元素释放的影响,以期为深入理解湿地植物分解过程及湿地生物地球化学过程进而强化对鄱阳湖湿地生态系统的有效保护提供科学认识。1、材料与方法1.1研究区概况本研究选择位于鄱阳湖南矶山湿地国家级自然保护区内代表性碟形湖泊湿地—白沙湖为试验地(图1)。该保护区位于鄱阳湖南部,地处赣江北支、中支和南支入湖冲击形成的三角洲前缘,海拔在12—16m,属亚热带暖湿型季风气候,夏季炎热多雨,冬季低温少雨,年平均降雨量为1387—1795mm,年最大降雨量为2452.8mm(1954年),年最小降雨量为1082.6mm(1078年),年平均蒸发量800—1200mm,多年平均气温17.6℃,最冷月平均气温5.1℃,最热月平均气温29.5℃。受鄱阳湖季节性水文节律的影响,研究区湿地有明显的丰、枯水期,呈现出典型的水陆相交错变化:丰水期(4—9月)除南山岛和矶山岛外其他洲滩均被洪水淹没;枯水期(10月—次年3月)湖水消退归入河道和碟形洼地,不同高程洲滩相继出露,呈现河、湖、洲交错景观,发育有大量肥沃土壤和草洲,其中灰化薹草(Carexcinerascens)集中成片分布,遍及整个湿地草洲,盖度为85%—100%。图1研究区位置与样点布设1.2试验设计与样品分析研究样地面积为200m×300m,由于灰化薹草分布面积最广,优势度最大,因此选择灰化薹草为分解实验材料。为了确保研究结果的典型性和科学性,根据地下水位将研究区样地划分为4个梯度带,由湖岸到湖心的方向平行分布,依次为GT-A(-50—-25cm),GT-B(-25—-15cm),GT-C(-15—-5cm),GT-D(-5—5cm)(图1),相邻梯度带间距约100m。每个地下水位梯度带内根据微地形差异设置3—6个监测样点,共计20个分解试验监测点,每个样点插入内径为2cm的PVC管(绕管壁等距离钻取4个内径约3.5mm的小孔,埋入地下约70cm,用铅垂线保证管壁与水平面垂直)用于测定地下水位,地下水位测量时间为2016年10月15日,11月15日,12月15日和2017年1月15日,每日9:00,13:00和16:00各测量一次,测量方法为用钢卷尺测量PVC管中水面与地表间的距离。具体方法参见文献[21]。分解实验采用尼龙网袋法。为避免分解袋中植物残体的非分解损失,同时保证不限制分解作用,选择了100目(0.15mm)、规格为15cm×15cm的网孔分解袋。2016年10月10日在远离水体的高地草洲上采集灰化薹草的成熟叶片,用去离子水冲洗后,剪成10cm长的小段,置于60℃烘箱中烘干至恒重,每个分解袋中放入5.00g烘干处理后的分解材料。将分解袋用竹竿固定在样点地表。分解试验开始于2016年10月15日,样品回收时间分别为试验开始后的第15、30、60、90、120天。分解袋带回实验室后,清除表面杂物并用去离子水冲洗干净后置于60℃烘箱中烘干至恒重,球磨仪研磨后测量其干物质质量、木质素含量以及碳、氮、磷元素含量。木质素的测定方法见文献[22],碳、氮元素含量用元素分析仪(ElementarVarioMaxCN;Hanau,德国)测定,磷元素含量用钼锑抗比色法测定。由于枯水期中期湿地洲滩不同地下水位土壤理化性质差异最为明显,而分解第60天植物分解过程受分解环境的影响最为显著[18,19],因此分别于2016年12月即枯水期的中期和分解第60天左右进行分解试验监测点处的土壤样品采集。利用土钻(内径5cm)对每个样点进行五点法混合取样,采集0—20cm土壤样品约500g装入聚乙烯自封袋内,编号后置于保温箱内运回实验室进行后续分析。土壤含水量采用烘干法测定,土壤容重(BD)采用环刀法测定,土壤pH采用水土比为2.5∶1的pH计测定(Sartorius,德国),土壤有机碳(TOC)采用重铬酸钾外加热法测定,土壤总氮(TN)采用元素分析仪测定(ElementarVarioMaxCN;Hanau,德国),土壤全磷(TP)采用钼锑抗比色法测定,土壤微生物量碳(MBC)和微生物量氮(MBN)采用氯仿熏蒸-K2SO4浸提法测定。1.3数据处理与统计灰化薹草分解过程中木质素的分解过程参数计算指标与方法为:实测瞬时残留率(Rt)计算公式[23]:Rt=Mt×LtM0×L0×100%式中,Rt表示t时刻残留率,Mt表示t时刻重量,Lt表示t时刻木质素的浓度,M0表示初始重量,L0表示木质素的初始浓度,t表示分解时间。瞬时分解速率(k)用Olson负指数衰减模型[24]计算:Mt=M0e-kt式中,k表示t时刻瞬时分解速率,k值越大分解速度越快。Mt表示t时刻重量,M0表示初始重量,t表示分解时间。碳、氮、磷元素在分解过程中的相对养分归还指数(RRI,Relativereturnindex)计算方法为[23]:RRIt=M0×C0−Mt×CtM0×C0×100%式中,Ct为t时刻某元素的浓度,C0为某元素的初始浓度,Mt表示t时刻物质重量,M0表示初始物质重量。用CRRI、NRRI和PRRI分别表示碳的相对归还指数、氮的相对归还指数和磷的相对归还指数。采用SPSS20.0进行数据统计分析,采用单因素方差分析不同地下水位环境梯度下的木质素残留率,碳、氮、磷元素相对归还指数及土壤理化性质的变化,多重比较方法选择用Dunnett′sT3比较检验,采用线性回归分析研究木质素分解速率与碳、氮、磷元素相对归还指数的关系,以及土壤木质素酶与碳、氮、磷含量的关系,采用CANOCO(CanocoforWindows5.0)进行冗余分析(RDA)研究土壤环境因子对木质素分解速率与碳、氮、磷元素相对归还指数的影响。制图选用OrigninlabPro8.0。2、结果与分析2.1不同地下水位环境木质素分解过程不同地下水位梯度木质素残存率随时间变化的方差分析表明(图2),分解第15天,木质素残留率随地下水位梯度的升高而增大,分解第30天,不同地下水位梯度木质素分解残留率差异性显著(P>0.05),而分解30d后,木质素残留率随地下水位的升高而显著降低(P<0.05)。不同地下水位梯度木质素分解速率随时间的变化均呈现先升高后降低的趋势(图2)。低地下水位梯度(GT-A和GT-B)木质素分解速率的最大值出现在15—30d,高地下水位梯度(GT-C和GT-D)木质素分解速率的最大值出现在30—60d。分解初期(20d之前),木质素分解速率随着地下水位的升高而降低,而分解30d以后,木质素分解速率随着地下水位的升高而升高。图2灰化薹草分解过程中木质素残存率及分解速率动态变化不同小写字母表示在0.05水平上具有显著差异2.2不同地下水位环境碳、氮、磷释放过程在整个分解过程中灰化薹草的碳、氮、磷元素都表现为净释放(RRIs>0;图3)。不同地下水位梯度间,CRRI在分解60d之前和120d均无显著性差异(P>0.05),而60d—90d表现为GT-D显著高于其他梯度(P<0.05);NRRI随时间的变化与CRRI趋势相似,在分解60d之前无显著性差异(P>0.05),60d—90d为GT-D显著高于其他梯度(P<0.05),而分解第120天则又表现为GT-A显著高于GT-C和GT-D(P<0.05),而GT-D与其他地下水位梯度差异性不显著(P>0.05);PRRI则在整个分解过程中均为随着地下水位的升高而显著增大(P<0.05)。随着分解时间的延长,CRRI呈先升高后降低的趋势,最高值出现在分解第60天;NRRI呈上下波动变化;PRRI则随着时间持续升高,60d后趋于稳定。图3灰化薹草分解过程中碳、氮、磷元素相对归还指数(CRRI、NRRI、PRRI)动态变化2.3木质素分解与碳、氮、磷释放的关系碳、氮、磷元素相对归还指数与木质素分解速率的线性回归分析表明(图4),分解0—60d,CRRI和NRRI与木质素分解速率的正相关性逐渐增强,分解60d后相关性减弱,分解第60天相关性显著程度最高(P<0.01)。PRRI在分解第3天随着木质素分解速率的增大而显著降低(P<0.01),分解60d以后则随着木质素分解速率的增大而显著升高(P<0.05)。在分解第60天对应的洲滩出露中期,随着地下水位的升高,土壤TOC、TN和土壤容重均逐渐降低(P<0.05),土壤TP表现为GT-D显著低于其他地下水位梯度(P<0.05),而其他梯度间差异不显著(P>0.05),土壤pH、含水量和土壤微生物量碳氮均随地下水位升高显著增大(P<0.05);微生物量碳、氮在GT-D显著高于其他地下水位梯度(P<0.05)。表1不同地下水位环境土壤参数对分解第60天灰化薹草分解过程中所有地下水位带木质素分解速率及碳、氮、磷的相对归还指数进行冗余分析发现,前2个排序轴累积贡献值为94.8%,其中第1排序轴解释了99.3%的变异(图5)。分解过程中,土壤环境因子对木质素分解速率及碳、氮、磷元素相对归还指数的作用一致。在所有环境因子中,MonteCarlo检验表明,木质素分解速率及碳、氮、磷元素的相对归还指数与土壤含水量,pH和微生物量碳、氮含量呈显著正相关性(P<0.05),而与土壤有机碳、全氮和土壤容重呈显著的负相关性(P<0.05)。总体上,在分解第60天,随着地下水位环境梯度升高,土壤含水量、pH和土壤微生物量升高,土壤有机碳、全氮、容重含量降低,灰化薹草木质素分解速率降低,碳、氮、磷元素相对归还指数减少。3、讨论地下水位环境梯度和时间是影响植物分解过程的主要因素。在整个分解时间段,不同地下水位环境梯度灰化薹草木质素残存率、分解速率以及磷元素的相对归还指数均具有显著差异,且在干湿交替环境下的高地下水位梯度条件下木质素的分解程度和磷元素的归还程度最高。这说明,较高的地下水位和间歇性的淹没会通过促进微生物代谢水平和物质淋溶过程来增加植物的分解速率[25]。Foulquier等[26]曾提出,较高的地下水位梯度能够为分解过程提供合适的微生态环境,因为适宜的土壤pH,土壤质地以及土壤微生物群落会促进有机碳的可利用性和周转速率,进而直接影响有机质的分解过程[27]。土壤微生物群落结构可能是影响植物分解过程中养分迁移转化和改变土壤生物、物理、化学进程的重要因素[28,29]。图5不同环境梯度灰化薹草分解参数与土壤环境因子的RDA排序不同地下水位环境梯度间,有机碳和总氮的相对归还指数差异不显著,这主要是因为在元素水平上,微生物在分解有机质的同时也会同化固定外界环境中的无机物,进而在分解残体中表现为元素的内外双向流动,而不是和木质素一样的单调减少[30]。但是整体上在灰化薹草的分解过程中,木质素的分解速率与碳、氮、磷元素表现出一定的相关性,且相关程度与分解时间密切相关,在分解第60—90天相关性最强。相关研究表明,相对湿润环境下的分解过程中微生物作用在30—90天迅速增强[31],并分泌出大量分解木质素的胞外酶[32],说明在灰化薹草分解第60—90天内,木质素在强烈的微生物作用下快速分解为小分子的化合物[4],这一时间段内有机碳的释放和归还主要来自于木质素的降解。在分解0—30d,由于有机碳和氮的释放主要靠糖类、低分子酚类等易分解物质[33],因此碳、氮元素的相对归还指数与木质素分解速率相关性较弱。在分解90天后,由于在微生物的同化作用下,外部碳、氮元素富集程度升高,而内部碳、氮的释放速率减慢[34],所以碳、氮元素相对归还指数与木质素相关性显著降低。Connin等[35]通过分析分解过程中δ13C的变化发现,早期分解过程以碳的释放为主,而分解后期外部碳元素开始逐渐富集。Oson等[36]研究发现,分解过程中磷元素的净释放量与木质素的净损失量具有显著的相关性,并且磷元素对木质素分解速率的敏感性比碳和氮更强,与本研究的结果一致,然而木质素降解对磷元素释放的影响机制至今尚不明确。总体上木质素分解速率与碳、氮、磷释放归还量变化的关系在一定程度上能够反映湿地植物分解阶段性差异。本研究冗余分析表明分解过程中木质素分解速率与碳、氮、磷元素的相对归还指数对土壤环境因子的响应特征相同,进一步说明了木质素分解与碳、氮、磷元素的释放过程具有协同性。地下水位变化是影响土壤理化性质分异的主要原因[37],湿地土壤在干湿交替条件下会产生土壤的膨胀和收缩,使土壤团聚体中有机质暴露,激发土壤微生物代谢活性,促进有机质的分解[38]。本研究中不同地下水位梯度土壤在洲滩出露中期差异性显著。高地下水位梯度环境下,水体波动引起的冲刷和筛选及淋溶作用导致土壤胶体形态变化,引起土壤pH、容重等性质的变化,使原来无法分解的有机质因团聚体崩溃而加速分解,因此土壤有机碳和全氮含量随着地下水位的升高而降低。洲滩出露中期,土壤微生物量碳氮均随着地下水位升高而显著增大,主要是因为土壤pH和含水量综合作用的结果[37]。较高的地下水位梯度能够为分解过程提供合适的微生态环境,因为适宜的土壤含水量、pH及土壤微生物活性的增强都会提高有机碳的可利用性和周转速率[22],进而促进了木质素的分解和碳、氮、磷等元素的释放。已有研究表明,木质素能通过提升土壤中水溶性磷等有效磷的含量来减少土壤对磷元素的固定[39],并且土壤中木质素的降解与土壤碳、氮含量显著相关。本研究从湿地植物分解过程入手,揭示了木质素分解速率与碳、氮、磷元素释放的关系,尤其是分解15d后,木质素的快速分解对磷元素的释放具有显著的促进作用,说明在土壤-植物系统中,木质素的含量及其分解过程对碳、氮、磷元素的迁移转化和生物地球化学循环具有重要的作用。然而,木质素如何通过影响元素的形态和转化来改变其在土壤及植物分解残体中的含量有待进步深入研究。4、结论地下水位环境梯度对灰化薹草分解过程中木质素分解速率受分解时间的影响,在分解第15天,木质素分解速率随地下水位升高而降低,15d后则随地下水位升高而显著升高,直至90—120d后趋于稳定;地下水位升高显著加强了灰化薹草分解过程中碳、磷元素的相对归还指数,而氮元素相对归还指数在不同地下水位梯度间的差异不显著;在分解第60—90d,碳、氮、磷元素的相对归还指数随着木质素分解速率的增大而显著升高,在分解第30天,碳、磷元素的相对归还指数也与木质素分解速率具有显著正相关性。总体上碳、氮、磷元素相对归还指数与木质素分解速率的相关性随着分解时间先增强,后减弱;较高的地下水位梯度能够为分解过程提供了适宜的pH、含水量及土壤微生物生物量等微生态环境条件,进而促进了木质素的降解和碳、氮、磷元素的释放和归还。参考文献：[14]胡茂林,吴志强,刘引兰.鄱阳湖湖口水位特性及其对水环境的影响.水生态学杂志,2010,3(1):1-6.[15]丁庆章,刘学勤,张晓可.水位波动对长江中下游湖泊湖滨带底质环境的影响.湖泊科学,2014,26(3):340-348.[16]胡振鹏,葛刚,刘成林,陈伏生,李述.鄱阳湖湿地植物生态系统结构及湖水位对其影响研究.长江流域资源与环境,2010,19(6):597-605.[17]吴琴,尧波,朱丽丽,幸瑞新,胡启武.鄱阳湖典型苔草湿地生物量季节变化及固碳功能评价.长江流域资源与环境,2012,21(2):215-219.[20]张广帅,于秀波,刘宇,张欢,张全军,李雅,段后浪.鄱阳湖碟形湖泊植物分解和水位变化对水体碳、氮浓度的叠加效应.湖泊科学,2018,30(3):668-679.[21]李丽,雷光春,高俊琴,吕偲,周延,贾亦飞,杨萌,索郎夺尔基.地下水位和土壤含水量对若尔盖木里苔草沼泽甲烷排放通量的影响.湿地科学,2011,9(2):173-178.[37]张广帅,于秀波,张全军,李雅,刘宇,段后浪.鄱阳湖湿地土壤微生物群落结构沿地下水位梯度分异特征.生态学报,2018,38(11):3825-3837.[38]孟伟庆,莫训强,胡蓓蓓,贺梦璇,李洪远.模拟干湿交替对湿地土壤呼吸及有机碳含量的影响.土壤通报,2015,46(4):910-915.[39]陈倩,穆环珍,黄衍初,杨问波.木质素对土壤磷素形态转化及对磷有效性的影响.农业环境科学学报,2003,22(6):745-748.张广帅,张全军,于秀波,闫吉顺,赵全民,蔡悦荫,宫玮.不同地下水位灰化薹草分解过程中木质素降解与碳、氮、磷释放的关系[J].生态学报,2020,40(09):3006-3015.基金：国家自然科学基金项目(41471088);国家海洋环境监测中心博士启动基金(2019-A-03).
