流域尺度水土保持生物措施减沙效益研究

1　引　　言黄土高原位于黄河中游 ,流经黄土高原的 9大支流进入河口镇至花园口的中游河段 ,总径流量为 187× 10 8ｍ3 ,输沙量占黄河流域总沙量的 90 % .建国以来以修建梯田坝地、造林种草、建谷坊、调整种植结构为主的水土保持工程、生物治理长年不懈 ,其目的在于减水减沙 ,改善流域生态环境 .但是 ,除部分中小尺度流域进入生态、经济效益良性循环的轨道外 ,大部分地区生态环境仍然十分脆弱 ,部分地区有恶化之势 .因此 ,定量评价水土保持工程、生物措施的治理效益是广为关注的研究课题 ,虽相关的研究成果不少 ,但是绝大部分工作是评价水土…     

赣江上游章水流域水沙变化及其归因分析

随着气候变化和人类活动的加剧,流域水沙情势发生变化,定量分析降水变化和人类活动的水沙效应是当前水科学研究的热点。本文基于赣江上游章水流域1961—2020年的水文、气象观测资料,采用累积量斜率变化率比较法、经验统计分析法定量分析了降水变化和人类活动对水沙变化的贡献率及人类活动的水沙效应对年降水量响应的临界值。结果表明：降水量和径流量均呈非显著下降趋势(P>0.05),输沙量呈极显著下降趋势(P<0.001);降水量和径流量的突变年份均为2003年,输沙量的突变年份为1998年;降水变化和人类活动对径流减少的贡献率分别为10.7%和89.3%,对输沙减少的贡献率分别为3.5%和96.5%;人类活动的减水、减沙效应对年降水量响应的临界值分别为1170 mm和858 mm。研究成果可为流域生态建设、水土资源优化管理提供科学依据。     

基于随机森林模型的茹河流域水沙变化及其影响因素研究

水沙关系不协调已成为制约黄河高质量发展的关键问题,揭示径流量与输沙量变化的主要影响因素,对适应生态经济发展条件下的水土流失治理工作具有指导意义。以茹河流域为研究对象,采用M-K检验和Pettitt检验分析茹河流域1989—2019年水沙变化趋势,运用双累积曲线法定量研究了降雨和人类活动对茹河流域径流量、输沙量变化的贡献率,通过随机森林回归分析探索了影响水沙变化的主控因子。结果表明：(1)茹河流域径流量和输沙量峰值对应情况良好,但在2012年之后径流量和输沙量呈现出明显分异,径流量逐渐增加,输沙量平稳减小。(2)茹河流域径流量突变点为1996年和2003年,输沙量突变点为1996年和2004年,影响期的两个阶段中,人类活动对径流量的贡献率分别为95.95%和77.52%,对输沙量的贡献率分别为78.31%和77.91%。(3)显著影响径流量的因子为农作物播种面积、耕地面积、草地面积和人口密度,显著影响输沙量的因子为植被覆盖度、水利环境投资占比、人口密度和工业总产值。(4)农作物播种面积和耕地面积的变化对径流量的影响均呈显著正相关,草地面积和人口密度对径流量响应关系为显著负相关;植被覆盖度...     

黄土高原丘陵沟壑区典型小流域水土流失治理模式

黄土高原是我国水土流失最为严重的区域,几十年的水土流失治理工作成果显著。当前生态文明建设已是我国国家战略之一,对黄土高原水土流失治理模式现状进行系统总结十分必要。对水土流失治理措施体系和流域土地利用情况系统分析,整理和总结了黄土高原6个典型流域的治理模式,包括治沟造地、防蚀固沙、生态农业发展、三大体系、梯田开发和水资源高效利用等6种模式。对比各治理模式发现,小流域综合治理模式是以环境治理为手段,从而达到改善人类生存和发展空间的终极目的,措施技术体系构建上体现了从坡面到沟道、工程措施植物措施相结合的特征;但不同的生态环境现状和社会经济条件下,治理模式有不同的表现形式。基于小流域水土流失治理模式的相似性和差异性,从自然地理条件和社会经济发展需求两方面探讨了流域水土流失治理模式,提出黄土高原水土流失治理模式形成机制框架,以期对当前黄土高原水土流失治理、黄河流域生态保护与高质量发展提供参考。     

近70年来不同地貌和植被类型区水沙特征演变及其对生态恢复的响应

黄土高原水土流失治理和退耕还林(草)等生态环境建设,对不同植被类型区和治理度下的水沙过程均产生了深刻的影响。以北洛河流域为例,基于1952—2019年实测水文数据,采用水文统计法,对比分析了上、中、下游等不同地貌和植被类型区生态恢复对水沙数量、径流调节、影响因素贡献等的影响及其差异。结果显示：(1)上、中、下游年均径流深分别为30、23、21 mm,产沙量为0.6、0.05、-0.02亿t。中游年均产流量占流域总水量的51%,上游年均产沙量占流域总沙量的90%以上。(2)上、中、下游径流深极显著减少,减少率分别为-0.28、-0.16、-0.43 mm/a;上游产沙量以0.02亿t/a速度极显著降低,中、下游产沙量不显著减少。(3)上、中、下游的削洪、补枯程度逐年代增加,与20世纪60年代比,21世纪10年代洪水期流量(Q₅)分别减少了63%、45%、60%,上游枯水期流量(Q₉₅)增加了68%,中游20世纪80年代枯水期开始产流。(4)与20世纪60年代比,上、中、下游21世纪10年代人类活动对产流量减少的贡献率分别为94.7%、26.7%...     

无定河流域不同地貌区径流变化归因分析

分析了无定河流域干流与其支流(黄土丘陵区的大理河和风沙区的海流兔河)的年径流变化(1960—2012)及其成因,并预测了其径流的变化趋势。结果表明:无定河及其不同地貌区支流海流兔河和大理河流域1960—2012年径流量均显著下降,但年降水量未发生显著变化;无定河和海流兔河流域年蒸散量未发生显著变化,仅大理河流域年蒸散量在1990s年代后期显著增加。无定河流域径流变化突变点发生在1979年和1996年,海流兔河流域径流变化突变点在1971年和1990年,而大理河流域径流突变点发生在1971年。人类活动对大理河流域1972—2012年径流减少的贡献大约占50%,对海流兔河流域1972—1990年和1991—2012年两个时期径流减少的贡献率分别为44.4%和82.4%。未来,无定河及其支流年径流量均呈持续下降趋势。归因分析表明大规模的水土保持治理措施减少了大理河流域侵蚀产沙量的同时,也在一定程度上减少了大理河流域的径流量,而过度的农田灌溉引水是海流兔河径流量下降的主要原因。因此,未来在大理河流域要优化现有植被建设布局,减少流域蒸散发,减缓径流下降;在海流兔河流域要进一步退耕还林(草),适当控制农田灌溉面积,提高灌溉用水效率,在减少灌溉用水的同时提高流域水源涵养能力。     

榆神府矿区典型小流域侵蚀产沙对退耕还林(草)及煤矿开采的响应

在脆弱的自然环境和强烈的煤矿开采活动共同影响下,黄土高原以榆神府矿区为核心的生态脆弱带土壤侵蚀极其严重。虽然退耕还林(草)工程的持续实施有效缓解甚至遏制了土壤侵蚀,煤矿开采却引发了新的问题,改变了小流域侵蚀产沙特征。本研究在该区域选择代表性坝控小流域——老爷满小流域,运用复合指纹识别法判别淤地坝沉积泥沙来源,旨在揭示退耕还林(草)和煤矿开采对榆神府矿区典型小流域侵蚀产沙特征的影响;从已建立的坝地运行期间的沉积旋回时间序列中选取退耕还林(草)之前10年(1990—1999年)和之后10年(2000—2010年),对比分析源地对坝地沉积泥沙的贡献差异及其变化特征。结果表明: 2000—2010年沉积旋回的年均厚度由1990—1999年的(42.7±29.9) cm减小到(16.6±13.7) cm,坝地采样点单位面积年均沉积量由579.9 kg·m^-2减少到245.8 kg·m^-2。研究期间小流域首要泥沙源地均为沟壁(平均贡献为79.2%±10.1%),但与1990—1999年相比,2000—2010年间受退耕还林(草)工程影响较大的坡面对坝地沉积泥沙的贡献率呈现增加趋势,年均贡献增加7.4%。结合区域的降雨和煤矿开采情况,上述现象表明在退耕还林(草)和降雨侵蚀力减小的共同影响下,流域的侵蚀产沙强度总体明显减弱;煤矿开采对坡面侵蚀产沙有重要的影响,导致退耕还林(草)对坡面水土流失的防治效果受到明显的影响。由此可见,煤矿开采增加了坡面的侵蚀产沙风险,受煤矿开采影响的区域是坡面水土流失治理的重点区域。     

黄土高原流域水沙变化研究进展

人类活动和气候变化是影响流域水文过程的两大驱动因素,径流输沙是流域水文过程的总体反映,变化环境下径流输沙的变化规律与成因分析是水文学和全球变化研究的热点问题。黄土高原是我国水土流失最严重的地区。20世纪50年代以来,黄土高原地区开展了大规模的生态环境建设和水土流失综合治理,显著改变了流域土地利用和植被覆盖。下垫面条件改变与气候变化综合作用,使得流域水沙情势发生剧变。围绕黄土高原流域水沙变化的时空尺度特征与驱动机制,总结了径流输沙和水沙关系变化特征的研究结果,归纳了径流输沙变化的归因分析方法与人类活动和气候变化影响的贡献分割结果,探讨了气候变化、植被恢复、水土保持工程措施以及流域景观格局对水沙变化的影响机制。未来应加强流域水沙演变的时空尺度特征特别是水沙关系非线性特征的定量研究,阐明极端事件对水沙动态的影响与贡献;开展水沙变化影响机制的多要素综合解析,发展耦合地表覆被动态特征和气候变化的降雨-径流-输沙模型,揭示生态恢复与水沙演变过程互馈机制;开展未来气候变化、社会经济发展和生态建设工程情景下水沙动态的趋势预测,为黄土高原生态综合治理和水资源管理与黄河水沙调控提供策略建议。     

黄土高原水土流失治理现状、问题及对策

黄土高原位于黄河中游地区,是世界上水土流失最为严重的区域之一;黄土高原水土流失治理问题历来受到国家高度关注,涌现了种类多样的水土流失治理模式,目前还缺乏对黄土高原区域尺度水土流失治理模式的总结和整理。本文总结和整理了建国以来黄土高原主要的四类水土流失治理模式,生物措施模式主要由退耕还林、荒山造林和封山育林工程组成,工程措施模式包括修建梯田和淤地坝,以及近期涌现的治沟造地工程;小流域综合治理模式主要体现在坡面、沟道系统整治,生物和工程措施相结合的特点,区域综合整理模式则强调对生态系统进行整体保护、系统修复和综合治理,达到生态、社会、经济可持续发展。基于黄土高原水土流失治理模式存在的问题,提出了黄土高原水土流失治理模式调整和优化建议,以期对黄土高原生态恢复建设和水土流失科学治理提供科学依据。     

黄土高原典型流域生态问题诊断与系统修复实践探讨——以山西汾河中上游试点项目为例

流域生态系统涵盖山、水、林、田、湖、草等自然生态各要素,对流域生态系统进行总体保护、系统修复和综合治理,是践行习近平总书记生态文明建设理念的重要抓手。黄土高原地处半干旱半湿润气候带,水土流失严重,生态环境脆弱,探索黄土高原典型流域生态保护修复的理论与实践,对于实现国家生态文明建设具有重要作用。从黄土高原典型流域生态系统特征及存在的生态问题出发,提出流域生态系统修复治理的四大核心理念：（1）宏观规划调控与局地修复调整相结合理念;（2）区域生态功能提升与社会经济发展相协调理念;（3）由“头痛治头、脚痛治脚”转向“整体把脉,系统治疗”理念;（4）由“开刀治病”工程治理向“健康管理”自然恢复逐步引导理念。以山西省汾河中上游为例,在梳理山水林田湖草生态保护修复工程项目实践设计思路的基础上,详细介绍了汾河中上游五大生态保护修复分区（管涔山汾河源头水源涵养及生物多样性保护区、汾河干流中上游水生态保护修复区、水土保持生态保护修复区、汾河上游水源涵养生态保护修复区和矿山生态环境治理修复区）及其修复任务。总结山水林田湖草生态保护修复基本经验,并对完善黄土高原典型流域生态保护修复体制建设与技术服务进行展望,以期为黄土高原乃至全国区域开展流域生态保护修复工作提供借鉴和参考。     

黄土高原典型坝系流域碳氮湿沉降与水体碳氮流失特征

本研究通过对黄土高原坝系流域的碳(C)、氮(N)湿沉降过程、降雨径流过程及其基流过程的动态监测,揭示C、N湿沉降对流域水体碳氮流失的贡献,探讨C、N流失负荷在降雨-径流中的分布,进一步评估黄土高原雨季水体碳氮流失状况。结果表明,研究区域碳、氮湿沉降通量分别为173.95 kg km~(-2)mon~(-1)和43.01 kg km~(-2)mon~(-1),而通过径流的碳、氮流失量为11.52 kg km~(-2)mon~(-1)和2.19 kg km~(-2)mon~(-1)。雨季C、N湿沉降对该地区水体C、N流失的贡献率分别为65.81%和100%,其中流水侵蚀引起的C、N流失贡献率为59.20%和56.16%。黄土高原主要以流水侵蚀为主,C、N主要集中在降雨径流后期,因此截留后期径流可以有效控制养分流失。     

黄土高原生态分区及概况

黄土高原地域广阔,水土流失区域差异显著。为了有效治理水土流失,评估水土流失治理技术和模式及生态恢复建设工程的成效性,需要对黄土高原进行区域划分。依据自然条件、水土流失治理技术和模式的区域性特征及差异,基于国家基础地理信息系统数据的县级行政界,对其进行合并,进行生态分区的划分,并分别统计其气候、地形地貌、植被特征及水土流失现状,以期为黄土高原水土流失治理技术和模式的改良优化提供依据。主要结论如下:(1)黄土高原分为黄土高塬沟壑区,黄土丘陵沟壑区,沙地和农灌区,土石山区及河谷平原区。其中黄土高塬沟壑区和黄土丘陵沟壑区分别划分为两个副区。(2)黄土高原的气候、植被、水土流失具有明显的分区差异。降水和植被覆盖度自东南向西北递减,二者在空间分布上具有很好的一致性,降水量大的分区,植被覆盖度也高。在年际变化方面,丘陵沟壑区B2副区降水量呈增加趋势,其他分区呈减小趋势,变化均不显著。80年代以来,黄土高原和各生态分区的植被覆盖度均逐渐增加,黄土丘陵沟壑区的增加量最大。各分区的面平均气温均呈非显著增加趋势,90年代以来增温明显。(3)1970年以来,黄土高原侵蚀产沙强度减弱趋势显著,至2002—2015年,多年平均输沙模数在0.13—3924 t km~(-2) a~(-1)之间,侵蚀强度最大为中度侵蚀(2500—5000 t km~(-2) a~(-1)),但面积较小,主要分布于第二高塬沟壑区的泾河流域。     

不同视角下基于InVEST模型的流域生态补偿标准核算——以渭河甘肃段为例

补偿标准是生态补偿项目研究的核心内容。为提高补偿标准的可操作性,以渭河干流甘肃段为研究区,基于流域生态补偿标准范围界定,在考虑耕地向多种林地转化情况下采用In VEST模型和加权法核算流域新增生态服务量,并以新增生态服务价值为补偿上限,农户机会成本为补偿下限,先从国家新一轮退耕还林政策导向视角,探讨15°—25°和25°以上坡耕地退耕还林的差异化补偿标准。再从利益公平分配视角,耦合机会成本投入与生态系统服务产出,探讨生态补偿净收益不同地区的差异化补偿标准。结果表明:(1)渭河干流甘肃段退耕还林,生态补偿标准范围为146.39—481.98万元km~(-2)a~(-1)。(2)重要水源地一级保护区和保留区迎水面15°—25°坡耕地全部退耕还林需补偿金额13.95—33.55亿元。25°以上坡耕地全部退耕还林需要补偿金额12.96—28.35亿元。(3)各县/区比较,退耕还林净收益麦积区最大,陇西县最小。其生态补偿范围分别为202.23—727.92、96.99—136.53万元km~(-2)a~(-1)。     

黄土高原浅沟空间分异特征

浅沟侵蚀是黄土高原重要的侵蚀类型之一,在区域尺度回答浅沟空间分布特征是浅沟治理的关键。基于系统抽样方法在黄土高原选取137个调查单元,综合运用亚米级影像目视解译方法和GNSS RTK野外高精度实测,并通过GIS空间分析方法,获取浅沟(特指瓦背状地形底部,下同)长度、密度、临界坡长、所在坡面土地利用等信息,阐明了目前黄土高原浅沟空间分布,分析了仍分布在耕地上的浅沟空间格局。结果表明：(1)基于亚米级影像目视解译的浅沟调查方法与野外实测方法得到的浅沟长度无显著性差异,前者可以用于区域尺度浅沟的快速调查;(2)黄土高原33.6%的调查单元存在浅沟,平均浅沟密度为3.41 km/km~2,最大为21.92 km/km~2,平均沟数密度约为65条/km~2;平均沟长63.31 m;发生浅沟的坡面临界坡长主要为40—60 m,平均为56.20 m;延安市西北部吴起、志丹、榆林市靖边一带浅沟密度最大;在不同治理分区内,黄土丘陵沟壑区平均浅沟密度最大,为5.32 km/km~2,浅沟侵蚀强度极强烈;黄土高塬沟壑区平均浅沟密度为2.07 km/km~2,浅沟侵蚀强度为中度。(3)目前38.3%的浅沟分布...     

DEM栅格分辨率和子流域划分对杏子河流域水文模拟的影响

以黄土丘陵沟壑区杏子河流域为研究区,采用数字滤波法分割了杏子河招安水文站1958—1974年的基流量,评价了SWAT模型在该流域水文模拟的适用性,并分析了DEM栅格分辨率和子流域划分对SWAT水文模拟的影响。结果表明,SWAT模型适用于该流域年河川径流、地表径流、基流及产沙量的模拟。当DEM栅格分辨率在20—150 m之间时,SWAT能有效地模拟年河川径流、地表径流、基流及产沙量,各水文要素模拟结果的R2和NSE分别在0.93和0.51以上,RSR在0.43以下;而当栅格分辨率大于150 m时,各水文要素的模拟效果存在差异。子流域划分对流域产流模拟影响较小,而对产沙模拟影响较大。当子流域提取阈值在12—100 km2之间时,不同的子流域划分对产沙量几乎没有影响,若超出该阈值范围,模型会低估产沙量。因此,可针对不同的水文要素选择合适的DEM和子流域提取阈值,以提高模拟精度和运行效率。     

基于RUSLE模型的中国土壤水蚀时空规律研究

RUSLE模型是计算土壤水蚀的经典模型,在大尺度研究时参数率定比较困难。基于气候、土地覆盖、地形特征等空间分异特征,对RUSLE模型的降雨侵蚀力(R)、植被覆盖与管理因子(C)、水土保持措施因子(P)进行了率定,估算了2000、2005、2010、2015年的中国的土壤侵蚀量。结果表明:(1)土壤侵蚀强度较大的地区集中在中国长江中下游平原区、云贵高原、黄土高原区、昆仑山山麓区域,占统计总面积的9.65%。(2)土壤侵蚀明显增大的区域面积达10.36×10~4km~2,分布于新疆农田区、四川盆地、云贵高原东南部、长江中下游平原和东北平原。(3)土壤侵蚀显著改善的区域分布于黄土高原南部、秦岭地区和东南沿海地区,面积约13.6×10~4km~2。通过对RUSLE模型参数的率定,阐明了全国尺度土壤水蚀的时空分异规律,可对不同地区制定相应的土壤修复措施提供依据。     

青藏高原植被生态系统垂直分布变化的情景模拟

如何模拟和揭示青藏高原植被生态系统垂直分布在全球气候变化驱动下的时空变化情景,对定量解析青藏高原陆地生态系统对气候变化响应效应具有重要意义。该论文基于Holdridge life zone （HLZ）模型,结合数字高程模型（DEM）数据,改变模型输入参数模式,发展了改进型HLZ生态系统模型。结合1981-2010（T0）时段的气候观测数据和IPCC CMIP5 RCP2.6、RCP4.5、RCP8.5三种情景2011-2040（T1）、2041-2070（T2）、2071-2100（T3）三个时段气候情景数据,实现了青藏高原植被生态系统垂直分布的时空变化情景模拟。引入生态系统平均中心时空偏移趋势模型和生态多样性指数模型,定量揭示了青藏高原植被生态系统在不同垂直带上的时空变化情景。结果显示：青藏高原共有16种植被生态系统类型;冰雪/冰原、高山潮湿苔原和亚高山湿润森林为青藏高原主要的植被生态系统类型,其面积之和占到了青藏高原总面积的56.26%;高山干苔原、亚高山潮湿森林、山地灌丛、山地湿润森林和荒漠等对气候变化的敏感性总体上高于其它类型;在T0-T3期间,青藏高原的高山湿润苔原、高山干苔原、荒漠呈持续减少趋势,平均每10年将分别减少1.96×10⁴km²、0.15×10⁴km²和1.58×10⁴km²;亚高山潮湿森林、山地湿润森林和山地灌丛呈持续增加趋势,平均每10年将分别增加3.42×10⁴km²、2.98×10⁴km²和1.19×10⁴km²;RCP8.5情景下青藏高原的植被生态系统平均中心的偏移幅度最大,RCP4.5情景下的偏移幅度次之,而RCP2.6情景下的偏移幅度最小。另外,在三种气候变化情景驱动下,青藏高原植被生态系统的生态多样性呈减少趋势。总之,未来不同情景的气候变化将直接影响青藏高原植被生态系统的时空分布格局及其生态多样性,气候变化强度越高,影响就越大,而且气候变化对青藏高原植被生态系统的影响呈现出从低海拔到高海拔递增的影响效应。     

退耕还林还草工程生态效应的地域分异特征

退耕还林还草工程作为我国投资最大、涉及面最广的一项重大生态工程,其生态经济社会效应是行业部门及学术界关注的焦点。选择退耕还林还草面积、植被覆盖度、土壤侵蚀量作为指标,依据遥感反演和模型模拟得到的结果,分析了2000—2015年县域退耕还林还草时空差异,退耕还林还草区域植被覆盖度与土壤侵蚀变化及其对县域生态状况变化的贡献,并基于规划目标评估了退耕还林还草工程的宏观生态效应及其地域分异特征。结果表明:(1)近15年,工程区耕地转林地面积12.75万km~2,耕地转林地区域植被覆盖度年增加0.32%、土壤水蚀和风蚀模数分别年减少0.43 t/hm~2和0.21 t/hm~2。(2)耕地转草地面积9.43万km~2,耕地转草地区域植被覆盖度年增加0.43%,土壤水蚀和风蚀模数分别年减少0.55 t/hm~2和0.94 t/hm~2。(3)工程区全县域植被覆盖度年增加0.17%,平均土壤水蚀和风蚀模数分别年减少0.13 t/hm~2和0.68 t/hm~2。(4)遥感估算结果与工程规划目标相比,退耕还林的面积完成率达到87%,工程区林草覆盖率增加4.8%—6.5%,工程区平均土壤水蚀和风蚀模数逐年减少,大部分县域15度及以上坡耕地退耕比例超过50%。工程在黄土丘陵沟壑、东北山地及沙地、云贵高原等区域凸显了提高植被覆盖度与土壤抗蚀效应的正面作用。     

黄土高原人工苜蓿草地固碳效应评估

苜蓿草地作为黄土高原草地的重要组成部分,在黄土高原草地生态系统碳循环占有重要地位。同时由于苜蓿具有固氮和固碳的双重作用,使它在低碳农业的发展中具有重要意义。本文以一年生（建植初期）、三年生（盛产期）、五年生（稳产期）和十年生（衰败期）的人工苜蓿草地为研究对象,利用地上地下生物量和土壤碳贮量,定量分析了人工苜蓿草地不同发育阶段的固碳能力;结合种植面积,对黄土高原区人工苜蓿草地的总固碳效应进行了评价。研究结果表明：苜蓿草地地上生物固碳能力以三年生和五年生最强,并显著高于其他年限;地下生物固碳和土壤固碳量均以十年苜蓿草地最高,且显著高于其他年限。总体而言,人工苜蓿草地以土壤固碳为主,并集中在0~80cm土壤,各年限0-80cm土层土壤固碳量分别占到总固碳量的41.02%、39.43%、41.56%、39.59%;而地下生物固碳主要发生在0~20cm土层,各年限该层生物固碳量分别占到地下生物总固碳量的45.74%、 55.68%、53.12%和 43.28%。苜蓿草地在整个生命周期中总固碳量随生长年限逐年增长,但增长速度不明显,各年限单位面积总固碳量分别为12.69kg/m2、13.49 kg/m2、13.31 kg/m2和14.83 kg/m2。估算得到黄土高原的人工苜蓿草地年总固碳量为1.430?1011 Kg（143.0 Tg）,其中地上、地下、土壤分别为5.43Tg、1.15Tg和136.4Tg。本研究结果为正确分析区域碳循环和发展低碳大农业提供了重要理论参考。