森林质量指数变化对江西省遂川县多年降雨径流的影响

根据遂川县域1989—2018年的天然降雨、径流序列数据以及不同年代遥感影像图,采用小波分析、遥感监督解译及森林质量等级判定等方法,对该县域降雨径流与森林质量指数变化的响应进行研究,为定量监测和科学评价南方红壤丘陵区水土流失与生态环境状况提供依据.结果表明:(1)遂川县域降雨径流随着年际变化表现出缓慢而不显著的下降趋势;(2)森林质量指数随解析年份的增加整体呈上升-下降-上升趋势;(3)随着森林质量指数的提高,径流深逐渐降低,森林质量提升具有显著滞洪作用.     

基于MODIS卫星火点数据的浙江省林火时空动态变化特征

  目的  研究森林火灾时空动态变化特征,有助于掌握森林火灾发生规律及变化趋势,为火险期划分和防火资源的合理分配提供科学依据。  方法  本研究基于MODIS卫星火点数据,运用数据统计、Mann-Kendall趋势检验法和滑动t检验,对浙江省2001—2016年森林火灾时空分布特征、变化趋势和防火期的变化进行研究。  结果  （1）2001—2016年浙江省森林火灾整体呈上升趋势,春季、夏季和冬季是林火发生主要季节。其中春、夏两季森林火灾显著增长,秋、冬两季森林火灾呈先上升后下降的变化趋势。浙江省传统定义火险期的林火数量呈下降趋势,而更多的林火在夏季发生,火险期由11月至翌年4月变为1月19日至9月1日,向夏季偏移和延长。（2）浙江省森林火灾空间分布呈现出由西南向东北递减趋势,其中杭州市、衢州市、丽水市和温州市森林火灾发生频率较高,且在防火期、春季和夏季呈上升趋势。  结论  研究建议浙江省加强西南地区及春、夏两季森林火灾监测和防火宣传,适当调整火险期和防火资源时空分配,将有助于抑制浙江省森林火灾的增长趋势。      

天目山林区景观格局时空变化及驱动因素分析

  目的  研究南方集体林区的景观格局变化特征及驱动因素,可为森林可持续经营提供科学依据。  方法  采用景观类型面积转移矩阵和景观格局指数分析方法,对比研究了浙江天目山自然保护区国有林和周边社区集体林2000、2010和2017年的景观格局变化特征。  结果  生态林是保护区和社区的优势景观。保护区整体呈生态林面积增加,经济林面积减少趋势,而社区生态林面积逐年减少,经济林面积不断增加;保护区的竹林和山核桃Carya cathayensis林向生态林转化,社区呈现生态林向经济林、建筑、道路转化,经济林之间相互转化,经济林和耕地相互转化的特征。保护区的景观格局变化不大,社区景观的破碎度、异质性和均匀性逐年上升,生态林的优势地位逐年下降。政策因素、市场需求变化和道路是导致保护区和社区森林景观格局变化差异的主要驱动因素。  结论  社区景观格局的变化会影响森林的生态和文化功能,可以从政策、旅游业发展和道路建设等方面制定相关措施。图1表8参17     

氮添加对大果木姜子凋落叶分解及养分释放的影响

  目的  研究不同施氮水平对大果木姜子Cinnamomum migao凋落叶养分分解的影响。  方法  于2017年1月,以药用植物大果木姜子人工林凋落叶为研究对象,将凋落叶清洗风干后装入分解袋中,每袋10.00 g。试验期间放置在不同施氮处理[对照(ck,0 g·m?2·a?1)、N₁(5 g·m?2·a?1)、N₂(15 g·m?2·a?1)、N₃(30 g·m?2·a?1)]样地中,每处理3个重复,分别于试验的3、5、7、9、11月采集凋落叶样品,测定凋落叶质量及养分质量分数的变化,分析氮沉降对大果木姜子凋落叶养分释放动态影响。  结果  凋落叶分解试验结束时,各施氮组凋落叶质量损失率整体小于ck,凋落叶残留率整体大于ck,且N₂与ck差异不显著,N₂、N₃差异显著(P＜0.05),ck、N₁、N₂、N₃凋落叶分解95%时所需的时间分别为2.973、3.626、3.285、3.671 a;各处理凋落叶碳质量分数总体均呈下降趋势,全氮质量分数整体上均呈先增加后降低的趋势,全磷与全钾质量分数变化趋势相似,为分解初期下降,之后整体上升,最后趋于平稳;各处理凋落叶碳、全磷、全钾残留率总体呈降低趋势,全氮残留率整体上呈先上升后下降趋势。其中分解各时段各施氮处理碳残留率均显著高于ck (P＜0.05);随着时间推移,ck处理氮的残留率呈先上升后下降趋势,各施氮处理的氮残留率整体上呈先上升后下降趋势;整个分解过程中各施氮组碳氮比均小于ck,且分解前期与ck表现出显著差异(P＜0.05)。  结论  添加氮不利于大果木姜子凋落叶的分解及养分的释放,且施氮越多抑制分解作用更显著。图4表1参33     

抵消机制背景下企业森林碳汇需求价格模拟

  目的  从森林碳汇需求角度,测算与分析企业减排成本、森林碳汇需求价格及政策影响因素,以期更全面了解森林碳汇的市场需求潜力,推进森林碳汇市场的发展。  方法  以北京、上海、湖北、广东4个碳交易试点省(市)为案例区,选择火电、化工、钢铁3个碳排放密集型代表行业,通过方向性距离函数方法计算89家2 759个样本减排单位的碳边际减排成本,进而采用罗宾斯坦恩博弈模型测算森林碳汇的需求价格,运用云模型测度不同行业和地区企业的森林碳汇需求价格并模拟政策变化对价格的影响。  结果  样本结果显示:火电、钢铁、化工行业的森林碳汇需求价格均值分别为631、556和575元·t?1,上海市、北京市、广东省和湖北省企业的均值分别为305、456、877和715元·t?1;云模型模拟分析结果表明:随着政府允许碳汇抵消比例的提高和碳汇补贴额度的增加,各行业和各省(市)企业对森林碳汇的需求价格上升,随着碳税征收率的提高,需求价格呈现下降趋势。  结论  不同地区不同行业的二氧化碳边际减排成本存在较大差异且逐年上升,森林碳汇是未来的减排趋势;目前火电行业企业已经正式启动全国碳排放权交易市场,云模型模拟结果显示:化工行业比其他2个行业更易受政策影响,因此合理的允许抵消比例与补贴政策组合下,火电、化工行业将会是未来森林碳汇的重大需求者。森林碳汇市场的发展潜力是巨大的。图4表5参28     

青岛崂山次生植物群落物种分布及多样性特征

  目的  为明确山东崂山地区森林群落的组成结构和多样性特征及其与环境因子的关系,为生物多样性保护和森林管理提供理论依据。  方法  用物种累计曲线验证抽样是否充分,根据物种重要值矩阵将崂山森林群落划分为不同的群落类型,用典范对应分析法（CCA）研究不同群落类型、不同层次植被的分布与地形及土壤条件的关系,并通过局部加权回归散点平滑法（LOWESS）拟合植被物种多样性与地形及土壤条件的相关性回归曲线。  结果  50个标准样地抽样较为充分;通过WARD聚类法可将调查样地划分为10种群落类型;排序结果显示海拔、坡位、坡向和土壤pH值均与群落的分布具有显著相关性,从物种分布来看,只有海拔因子是显著的影响因子。在11个地形和土壤因子中,海拔和土壤速效磷含量与群落总体物种多样性具有显著相关性;随着海拔的升高,乔木、灌木、草本植物的Shannon-Wiener指数均表现为先增加后减小的趋势;乔木、灌木的Pielou指数随海拔的升高略微增加;当速效磷含量达到30 mg/kg时,Shannon-Wiener指数和Pielou指数均显著下降,随后基本保持平稳。  结论  地形和土壤条件对崂山群落类型的分布具有一定的影响,证明了生境过滤的作用。土壤速效磷含量和海拔是崂山森林群落物种多样性的重要限制因子。      

神农架林区景观格局时空演变及其驱动力分析

  目的  20世纪90年代以来神农架林区森林保护工程陆续实施,使得神农架林区景观格局发生显著变化。需要分析这些变化及其驱动机制,以探究保护工作成效,同时为林区未来生态环境保护与土地合理利用对策的制定提供参考。  方法  在遥感(RS)和地理信息系统(GIS)技术支持下,以神农架林区为研究对象,选取1996、2007、2018年共3个不同时期的遥感影像为数据源,利用土地利用转移矩阵、景观格局分布重心模型与景观格局指数,定量分析了神农架林区1996?2018年景观格局的时空演变特征及驱动力。  结果  1996?2018年,林地与水域景观处于稳定向好状态,林地是研究区最主要的用地类型,占林区总面积的90%以上;草地与耕地面积总体呈减少趋势,其他土地用地类型向建设用地转化的趋势较明显。神农架林区总体上景观破碎度有所降低,而东北部区域与中南部区域的破碎化程度加剧,景观多样性也随之呈现扩散趋势。经济、人口、政治等是引起景观格局变化的主导因素。  结论  神农架林区退耕还林等生态保护与恢复工程初见成效,景观格局总体趋于整合,局部破碎化加重而景观多样性有所增加,社会经济空间与生态空间相互融合;近年来神农架大力开发的旅游产业造成了人工斑块对自然景观的侵占,对林区生态安全产生威胁。图4表5参32     

基于MODIS时间序列数据的竹林地上生物量估算

  目的  基于浙江省中分辨率成像光谱仪(MODIS)时间序列数据,对浙江省竹林地上生物量进行估算,为竹林碳汇遥感监测提供参考。  方法  以MODIS叶面积指数(LAI)、增强型植被指数(EVI)和比值指数(RVI)时间序列数据为变量,利用随机森林模型筛选变量,采用支持向量回归(SVR)模型估算研究区竹林地上生物量。  结果  随机森林模型共筛选出43个对竹林地上生物量影响最大的变量;基于43个变量,采用radial核函数构建的SVR模型预测能力最强,模型训练精度和测试精度分别为0.76和0.72,均方根误差分别为5.15和8.03 Mg·hm?2。浙江省全省竹林地上生物量均值为7.85 Mg·hm?2,总地上生物量为3.31×107 Mg;浙江省竹林地上生物量在各市具有明显的差异性,其中,湖州市、杭州市、金华市、绍兴市和宁波市的竹林地上生物量均值均大于全省均值,湖州市竹林地上生物量均值最大,为13.56 Mg·hm?2,舟山市地上生物量均值最小,为5.72 Mg·hm?2。  结论  耦合了MODIS LAI、EVI、RVI时间序列数据的SVR模型可实现浙江省竹林地上生物量较高精度的估算。图3表1参31     

基于Budyko理论定量分析窟野河流域植被变化对径流的影响

  目的  植被变化可影响流域下垫面条件,改变水文过程和水资源可持续利用量。揭示植被时空变化格局,明确其对径流变化的贡献是流域水资源管理规划的前提,也是植被恢复工程实施效果评价的依据。  方法  本研究以植被恢复的典型半干旱区窟野河流域为研究区,基于遥感植被指数NDVI和水文、气象数据,通过基于Budyko假设的弹性系数法,分析了1982—2015年植被变化对径流变化的影响。  结果  （1）径流突变点在1998年左右,突变点前后,径流和人类取耗水变化率分别达到?55.5%和92.3%,存在显著性差异（P < 0.05）;降水、干燥指数分别增加2.6%和1.9%,前后两个时段差异不显著（P > 0.05）。（2）流域植被指数NDVI在突变点前后差异显著,88.7%以上区域呈显著增加趋势（P < 0.001）,尤以下游明显,表明植被得到了极大恢复。Budyko弹性模型中表征流域属性信息的参数最优值为2.37,拟合出来的参数系数和常数项分别为20.407和?1.868（P < 0.05）,表明拟合效果较好,适用性较强。（3）植被变化对径流的贡献率达到84.9%,气候变化和人类活动影响分别为9.3%和5.8%。  结论  窟野河流域内植被得到了极大的恢复;植被变化成为决定径流变化的主导因素,其次为气候变化和人类活动。植被变化的长期生态水文效应有待进一步加强研究,本研究结果可为生态恢复工程实施和政策制定提供参考依据。      

高山松地上生物量遥感估算的不确定性分析

  目的  采用遥感数据估算森林地上生物量仍存在一些不确定性问题,研究估算过程中的误差来源及其占比,对提高森林地上生物量的估测精度具有重要意义。  方法  从遥感影像提取因子,结合高山松Pinus densata外业调查数据,建立多元线性回归、梯度提升回归树、随机森林等3种地上生物量估测模型,对样地尺度与3种模型的不确定性进行分析和度量。  结果  ①高山松单株生物量模型不确定性为16.43%,样地尺度的不确定性为7.07%;②多元线性回归模型残差不确定性为34.86%,参数不确定性为21.30%,与样地不确定性合成后总不确定性为41.45%;③非参数模型中,梯度提升回归树估测高山松地上生物量的总不确定性为23.12%,随机森林为19.42%。  结论  3种遥感估算模型中,非参数模型的不确定性明显低于参数模型。相较于样地尺度,遥感估算模型的不确定性对地上生物量估算精度的影响较大。图3表3参26     

洪涝诱发灾害干扰下受损恢复林地土壤颗粒的组成变化特征及多重分形分析

  目的  土壤颗粒组成的分形维数可得到更细致的土壤颗粒分布信息,对准确了解洪涝诱发灾害干扰下受损林地恢复过程的土壤颗粒变化规律,加深灾害干扰前后林地自然恢复变化过程的认识等具有重要的理论意义。  方法  采用以时间代空间的研究方法,选取次生阔叶林、杉木林、毛竹林3种林地,采集不同受损状态（未受损、刚受损、受损恢复7年）不同土层的土壤样品进行粒径分析,基于分形理论计算各项多重分形参数。  结果  灾害干扰对林地的土壤颗粒分布存在显著影响;灾害干扰后不同恢复阶段、不同土层深度对3种林地均有影响,不同林地对干扰的响应存在差异。综合分析3种林地不同恢复状态的粒径组成和多重分形参数可知:与未受损林地比较,受损恢复林地的土壤质地较差。3种林地在不同受损状态下不同深度土壤颗粒组成的变化趋势不同。  结论  本研究从土壤颗粒组成变化的角度阐明了洪涝诱发灾害对次生阔叶林、杉木林、毛竹林的影响,可为灾害干扰前后土壤侵蚀的防治和恢复提供理论依据。      

基于高空间分辨率时间序列影像的对象级侵占林地图斑快速检测

【目的】为掌握森林资源动态变化情况,及时、快速、准确地发现侵占林地地块,并解决主流遥感变化检测方法对数据源和时相一致性要求高、人工干预多、过程繁琐等应用瓶颈,采用一种基于高空间分辨率时间序列影像的多尺度对象级分割和变化提取方法,对主流方法的分类和检测两个过程进行了融合和简化。【方法】以陕西省白水县为研究区,采用GF-1和ZY-3卫星数据源,将前后两期遥感影像波段拆分和重组形成时间序列影像,对时间序列影像进行多尺度面向对象的分割,通过分割结果的光谱变化值统计学抽样判断临界点并制定提取阈值,再利用NDVI变化值对结果进行优化。【结果】以人工目视解译结果作为参照,该方法的检测精度达86.2%。在成功检出的侵占林地图斑中,形状吻合较好或基本吻合的图斑占48.8%。【结论】该方法能够实现侵占林地图斑的快速检测,在检测效率、精度和适应性方面可满足大范围、多时相、混合数据源森林资源监测工作的实际应用需要。     

密云水库上游地区不同雨型对坡面径流特征的影响

  目的  探讨不同自然降雨雨型对径流过程的影响。  方法  根据降雨数量特征及雨峰集中时间对密云水库上游地区自然降雨进行二级分类,分别探究不同雨型及过程雨型下的径流特征。  结果  ①根据降雨量(P)、降雨历时(t)、最大30 min雨强(I₃₀)可将研究区降雨划分为3类:P₂t₃I₂雨型(中雨量、长历时、中雨强)、P₁t₂I₁雨型(小雨量、中历时、小雨强)、P₃t₁I₃雨型(大雨量、短历时、大雨强),其中,P₁t₂I₁雨型在研究区发生频次最高,占比81.01%。4种过程雨型的发生频次从大到小依次为均匀型、前峰型、后峰型、中峰型。②不同雨型下的径流特征差异明显,从单次降雨看,P₃t₁I₃雨型单次产流能力最强, P₂t₃I₂雨型在研究区的径流贡献率最大。③不同过程雨型对径流特征影响较大,降雨量集中在降雨前期最容易产流。④降雨量对径流深的影响最大,I₃₀对径流深、径流系数及径流深峰值的影响均比较大,且对径流深峰值的影响要大于对径流深的影响。  结论  坡面径流特征对降雨雨型及过程雨型敏感性很高,采用合理的雨型划分方法研究雨型—径流关系可提高水力侵蚀研究精度。图4表5参27     

浙江松阳县生态公益林群落分类排序及优势种种间关联分析

  目的  对浙江松阳县生态公益林群落进行分类和排序,探讨多样性变化规律及种间关系,有助于揭示植物群落和多样性在复杂环境梯度上的分布规律。  方法  以松阳县重点公益林128个固定监测样地为研究对象,利用双向指示种分析(TWINSPAN)、典范对应分析(CCA)、物种多样性指数及种间联结分析等对该县公益林群落特征展开研究。  结果  TWINSPAN结果显示:松阳公益林群落可分为6种群丛类型。CCA结果显示:松阳县公益林群落分布和类型在环境梯度上呈现明显的分异规律,其主要受到海拔、土壤类型和坡向的影响。物种多样性结果表明:针阔混交林和阔叶林为主导的群丛Shannon-Wiener指数和Simpson指数显著高于其他群丛,以毛竹Phyllostachys edulis为主导的群丛多样性显著低于其他群丛。种间联结分析显示:样地内大多数乔木树种呈现不显著联结,正负关联比为1.9,说明群落整体稳定。  结论  松阳县公益林应进行差别化管理,积极推广针叶林混交化和阔叶化。可考虑对不同立地条件不同林分的公益林进行不同的补植措施,对毛竹林进行边缘控制,对一些能适应各类环境的群落组合进行适当推广。图4表3参37     

昕水河流域“水文法”基准期降水输沙相关关系统计方法研究

  目的  基于昕水河流域日尺度降雨和输沙数据,分析1958—2015年昕水河流域降水与输沙变化趋势及相关关系,比较年降雨量、汛期降雨量、侵蚀性降雨量、轻度产沙降雨量、中度产沙降雨量、强度产沙降雨量等6种统计标准和降水空间分布对“水文法”基准期分析降水输沙相关关系的影响,筛选最佳的降水输沙相关关系统计方法,提高“水文法”的计算精度。  方法  利用Mann-Kendall非参数检验方法分析降水量与输沙量变化趋势,双累积曲线法判定“水文法”基准期,Pearson相关分析法比较水沙关系相关性结果。  结果  昕水河流域年降水量无显著变化趋势,年输沙量呈显著下降趋势;6种统计标准下的水沙关系相关性由高至低依次为:强度产沙降雨量、年降雨量、侵蚀性降雨量、中度产沙降雨量、轻度产沙降雨量、汛期降雨量;低降雨强度但高降雨量降雨事件对研究区产沙贡献小,剔除产沙无效降雨能提高水沙关系相关性;午城、黄土、蒲县雨量站水沙关系相关性较好,大宁、桑峨雨量站水沙关系相关性较差。  结论  昕水河流域年输沙量变化与三北防护林及淤地坝建设等地表人类活动变化直接相关;研究区的产沙主要发生在夏秋季的强降雨事件中,强度产沙降雨量用于降水输沙相关关系研究结果最优,汛期降雨量不适用于降水输沙相关关系研究,产沙无效降雨会增加水沙关系研究的偏差,在缺乏降雨强度准确数据时,应选择降雨历时短、降雨量大的产沙有效降雨进行水沙关系分析;在流域尺度,降水空间分布对降水输沙相关关系有影响,距流域出口站距离较近雨量站降水数据的水沙相关关系较好。      

基于无人机激光雷达的天然林林隙识别

  目的  自然干扰引起的森林冠层林隙是天然林更新动态的主要驱动力，林隙的分布、形状和范围可以影响光照和土壤水分等生态因子。林隙的识别和特征描述对于理解森林的动态变化有着重要的意义。  方法  以云南省普洱太阳河保护区无人机飞行区域为研究区，根据无人机激光雷达点云数据提取冠层高度模型；然后使用固定阈值法、相对高度阈值法和面向对象分类法对冠层高度模型数据进行林隙识别，通过图像目视解释获得独立验证样本进行精度评估；最后精度选取最优方法提取的林隙描述其空间特征。  结果  固定阈值法的总体精度为92.00%，高于相对高度阈值法（66.00%）和面向对象分类法（88.00%）。研究区内林隙主要以中小林隙为主，干扰事件较少；研究区林隙的形状指数均值为1.97，多为形状指数较小、边缘效应不太明显的林隙，并且林隙的空间分布为聚集分布。  结论  利用无人机激光雷达数据和固定阈值法可以准确绘制出小范围亚热带天然林的林隙空间分布特征。      

戴云山南坡不同海拔森林土壤优势细菌群落特征及影响因素

  目的  探讨土壤细菌群落在戴云山保护区不同海拔梯度（900 ~ 1 500 m）分布特征,为理解海拔影响森林土壤结构和功能提供理论依据。  方法  基于高通量测序探讨不同海拔土壤细菌群落组成及多样性,并分析环境因子对土壤优势细菌群落结构的影响过程。  结果  （1）随海拔升高,土壤全磷含量总体呈逐渐递减;土壤有效磷含量整体呈单峰模式;土壤全碳和全氮含量呈双峰变化趋势。（2）不同海拔中土壤细菌优势菌门为变形菌门、酸杆菌门和放线菌门（相对丰度 > 10%）。（3）不同海拔梯度土壤细菌多样性指数如物种数、Chao1指数、Shannon-Wiener多样性指数和ACE指数沿海拔梯度呈先上升后下降的趋势,均在1 100 m处达到峰值,且达到显著水平（P < 0.05）。（4）微生物共现网络分析表明戴云山不同海拔土壤优势细菌群落具有明显模块化结构,关键类群包括变形菌门、酸杆菌门、放线菌门、拟杆菌门和疣微菌门的部分属,其中变形菌门的细菌关键类群最多。  结论  海拔、坡度、pH值、土壤全氮、水解氮和土壤有效磷是驱动不同海拔森林土壤优势细菌群落结构及多样性的主要因子。      

基于随机森林法的农作物高光谱遥感识别

  目的  不同农作物种类光谱差异小,通过探测众多窄波段范围的细微差别,提取区分不同农作物的特征波段,是目前实现农作物高光谱遥感识别的重要途径。如何提取区分不同农作物的特征波段,进而实现农作物的精确识别是一个挑战。近来出现的随机森林方法在多变量目标的分类识别方法展现了优势,为解决这一难题提供了一个新手段。  方法  利用随机森林法与传统方法分析杭州地区8种典型农作物的反射光谱,提取特征波段并进行分类,对比不同方法的识别效果。  结果  不同作物的反射光谱及其一阶微分、二阶微分、倒数的对数、去包络线法所提取的特征波段只能区分部分作物;随机森林法无需对反射光谱预处理,直接对全波段反射光谱数据处理,不仅筛选出了区分不同作物的特征波段,且运用所选择的波段对作物进行随机森林分类的效果也是最优的。  结论  随机森林法选择的波段（550、2 490、370、770、560、380、540、530、570、350 nm）不仅能区分不同作物,还能反映农作物生化属性的不同,使得用于分类的波段及分类方法体现了不同作物间物化性质的不同,在展现高光谱遥感识别农作物优势的同时,也为大面积农作物遥感精细分类提供借鉴。     

山西太岳山油松林土壤质量与有机碳稳定性随林龄的演变特征

  目的  通过分析土壤质量和有机碳稳定性随林龄的变化特征,为合理评估森林生态系统功能的恢复进程提供理论支撑。  方法  选取山西太岳山区域4种林龄（20年、40年、80年和110年生）的油松林,分析其表层土壤的有机碳及其组分、黏粒含量、根系生物量和土壤酶活性等理化指标,以有机碳与黏粒含量的比例表征土壤的物理质量,以惰性有机碳与活性有机碳含量的比值表征土壤有机碳的稳定性。  结果  土壤的水分含量、土壤酶活性、铁铝氧化物含量都随着林龄的增长呈现增大趋势;土壤质量和有机碳稳定性都随着林龄的增长而增大,在0 ~ 10 cm土层,以40年生林地的土壤质量最低,为0.12,110年生林地的最高,为0.40,有机碳稳定也是在40年生林地最低,为2.69,110年生林地最高,为6.72,土壤质量与有机碳稳定性间存在极显著正相关关系。  结论  络合态铁、铝氧化物含量对土壤有机碳稳定性的提高发挥了积极的促进作用,而土壤水分、根系生物量和土壤酶活性则对土壤质量的改善发挥的作用更强一些。