环境科学  2019, Vol. 40 Issue (2): 961-969   PDF    
复杂地形条件下根系对土壤有机碳的贡献
张彦军1,2, 郭胜利1     
1. 西北农林科技大学水土保持研究所, 杨凌 712100;
2. 宝鸡文理学院地理与环境学院, 陕西省灾害监测与机理模拟重点实验, 宝鸡 721013
摘要: 在地形条件复杂的地区,量化根系对土壤有机碳的贡献对科学评价水土流失区的土壤碳储量具有重要意义.本研究在黄土高原丘陵沟壑区的砖窑沟小流域内,基于地貌类型(梁峁坡、沟坡和沟谷)和植被措施(农田、林地和草地措施)两大因素采集土壤和根系样品,在流域尺度上研究根系密度(FRD)对土壤有机碳密度(SOCD)的贡献.在砖窑沟小流域内,地形、植被措施和土层厚度及其交互作用显著影响SOCD和FRD的空间分布.SOCD和FRD在不同地形部位下均呈现出沟谷>沟坡>梁峁坡的趋势,在不同植被措施下均呈现出林地措施>草地措施>农田措施的趋势,在不同土层厚度上均呈现出表层(0~20 cm)大于下层(20~100 cm)的趋势.此外,FRD对SOCD的影响显著(P < 0.05),SOCD随着FRD增加呈现出对数增加的趋势,且不同地形和植被措施下的根系-碳转化效率差异显著(P < 0.05).在农田措施下,沟谷(0.87)的根系-碳转化效率均是沟坡(0.43)和梁峁坡(0.43)的2.0倍;在草地措施下,沟坡(0.57)的根系-碳转化效率分别是沟谷(0.45)和梁峁坡(0.27)的1.3倍和2.1倍;在林地措施下,梁峁坡(0.56)的根系-碳转化效率是沟坡(0.44)的1.3倍.因此,在砖窑沟小流域内,从增加根系-碳转化效率的角度而言,沟谷适合进行农业生产,沟坡适宜进行退耕还草,而梁峁坡适合进行退耕还林.
关键词: 土壤有机碳      地貌类型      植被措施      根系-碳转化效率      黄土高原     
Contribution of Root Biomass to Soil Organic Carbon Under Complex Landforms Conditions
ZHANG Yan-jun1,2 , GUO Sheng-li1     
1. State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Northwest A & F University, Yangling 712100, China;
2. Shaanxi Key Laboratory of Disaster Monitoring and Mechanism Simulation, College of Geography and Environment Engineering, Baoji University of Arts and Sciences, Baoji 721013, China
Abstract: Quantifying the contribution of fine root density (FRD) to soil organic carbon density (SOCD) under extremely complex landforms is of great significance for the evaluation of carbon storage at an eroded small catchment in hilly regions of the Loess Plateau. Soil and root samples were collected from typical landforms (a ridge slope, gully slope, and valley bottom) and different vegetation types (cropland, grassland, and woodland) to investigate the contribution of FRD to SOCD at the Zhuan Yaogou watershed The spatial distribution of SOCD and FRD was influenced significantly by landforms, vegetation type, soil depth, or their interactions. SOCD and FRD tended to go in the following order:valley bottom > gully slope > ridge slope on different landforms, Woodland > grassland > cropland on different types of vegetation, and surface > lower depths at different soil depths. Additionally, FRD had a significant effect on SOCD (P < 0.05), SOCD increased logarithmically with FRD, and there was a significant difference in root-carbon conversion efficiency between different landforms and types of vegetation (P < 0.05). The root-carbon conversion efficiency on the valley bottom (0.87) was 2.0 times higher than those on ridge slopes (0.43) and gully slopes (0.43) on cropland; the root-carbon conversion efficiency on gully slopes (0.57) were 1.3 and 2.1 times greater than those in valley bottoms (0.45) and ridge slopes (0.27), respectively, on grassland; and the root-carbon conversion efficiency on ridge slopes (0.56) was 1.3 times greater than that on gully slopes (0.44). Therefore, from the perspective of increasing root-carbon conversion efficiency at the Zhuan Yaogou watershed in the hilly region of the Loess Plateau, valley bottoms are suitable for agriculture, gully slopes are suitable for returning farmland to grassland, and ridge slopes are suitable for returning farmland to woodland.
Key words: soil organic carbon      landforms      vegetation types      root-carbon conversion efficiency      Loess Plateau     

黄土高原地区地形破碎、沟壑纵横, 水土流失严重[1], 严重的水土流失导致其土壤碳库流失和土壤退化[2], 而水土流失治理可显著改善土壤理化性状和提高土壤碳储量[3, 4].而土壤碳库的恢复对维持地区的可持续发展和缓解气候变化具有重要意义[5, 6].因此, 在黄土高原地区, 研究地形和植被措施对土壤有机碳(SOC)时空变异性的影响, 对科学评价该区的土壤碳储量具有重要意义.

在地形条件复杂地区, 地形控制着地区的水土流失状况[1, 7], 支配着地区水热资源的空间分布[8, 9], 决定着地区植被群落的空间配置[10].因此, 植被措施、土壤微气候和土壤侵蚀均不同程度受到地形的制约, 从而导致有机碳储量在空间分布上的不均一[11, 12].随着植被的恢复, 不同地形部位的景观发生了巨大的改变[11, 13], 致使不同地形部位的SOC差异显著[14~16].除了地形因素以外, 植被措施/土地利用方式转化也是影响SOC的重要因素[17, 18], 且其主要通过改变进入土壤中有机物质的数量和质量影响SOC的时空变异性[19, 20].根系生物量占到陆地生态系统净初级生产力的70%以上[21], 而细根(d < 2 mm)的生物量占到总根系生物量的73%[22].此外, 大量的研究结果显示根系凋落物对土壤有机碳储量的贡献远远大于地表凋落物对土壤有机碳储量的贡献[19, 20, 23].但是, 目前大多数研究仅在相似或者单一地形条件下探讨不同植被条件下根系生物量对SOC的影响[3, 24~28], 但是同一植被措施下不同地形条件下根系生物量对SOC影响的研究很匮乏.

基于此, 本研究在黄土高原丘陵沟壑区的砖窑沟小流域内, 基于地貌类型(梁峁坡、沟坡和沟谷)和植被措施(农田、林地和草地措施)两大因素, 在流域尺度上采集土壤和根系样品, 探讨细根对土壤有机碳储量的贡献, 以期为合理估算地形复杂地区的碳储量提供数据支持和理论依据.

1 材料与方法 1.1 研究区域概况

砖窑沟流域(39°11′06″~39°13′47″N, 111°12′03″~111°19′28″E)位于晋西北河曲县沙坪乡, 地势东高西低, 海拔845~1 244 m, 主沟方向为东西向, 长约14.2 km, 支沟从南北两面汇入, 南北剖面呈凹形, 分水岭与沟底高差达150~200 m, 起伏剧烈.流域(黄土丘陵沟壑区)沟壑密度为6.24 km·km-2, 沟谷切割深度一般在50~100 m之间, 有的高达200 m左右.流域总面积29.11 km2, 地貌类型由梁峁坡、沟坡和沟谷组成, 三者在流域内所占的比例分别为47.7%、47.3%、5.0%.流域内坡度在0°~15°范围之间的土地面积占整个流域面积的40.40%; 15°~25°之间占29.81%; 25°~35°之间占12.44%; >35°的占17.35%.

流域属大陆性季风气候, 年平均降水量447.5 mm, 蒸发量1 913.7 mm, 年内降雨主要集中在6~9月, 其中7﹑8月占全年降水量的52.3%, 并多为大雨和暴雨. ≥10℃的积温天数为140 d, 无霜期为156 d, 全年日照达2 856 h.多年平均大风日数87.2 d.以春季最多, 占41.64%.流域属于黄土高原典型的水蚀和风蚀交错区, 地面割裂度为55.1%, 年土壤侵蚀模数>20 000 t·km-2, 属于剧烈侵蚀.流域的土壤以栗褐土为主, 土壤层次发育不明显, 多表现出母质性状, 质地偏轻, 有机质含量一般在0.50%左右, 全N﹑全P仅0.04%和0.05%, 全K含量2.00%左右.

流域植被措施多样, 主要农作物有糜子(Panicum miliaceum L.)﹑玉米(Zea mays L.)﹑谷子(Setaria italica)等; 林地包括乔木林和灌木林, 且主要为人工林, 乔木林地主要有刺槐(Robinia pseudoacacia)林、小杨树(Populus simonii Carr.)林、沙枣(Elaeagnus angustifolia L.)林等, 灌木林有柠条(Caragana korshinskii Kom)﹑沙棘(Hippophae rhamnoides L.); 草地分为人工草地和天然草地两种类型, 其中人工草地主要有苜蓿(Medicago sativa L.)、香蒲(Typha orientalis Presl)、沙打旺(Astragalus adsurgens Pall.)等, 且主要分布在梁峁坡上, 而天然草地主要有草木樨(Melilotus suaveolens Ledeb.)、冰草(Agropyron cristatum L. Gaertn.)等, 且主要分布在沟坡上.

1.2 测定项目及方法 1.2.1 土壤和根系样品的采集

基于小流域内的地貌类型和植被措施两大因素, 采用“分层次采样的方法” (0~20、20~40、40~60、60~80、80~100 cm, 共5层), 采集土壤和根系样品.将地貌类型[梁峁坡(对照)﹑沟坡和沟谷]作为一级层次, 将同一地貌类型条件下的不同植被措施[农田(对照)、林地措施﹑草地措施]作为二级层次.在流域内, 以横断面思路为依据, 采集土壤和根系样品, 并充分考虑空间异质性, 在沿途进行补点采样. 3种地貌类型在野外的判别依据为:沟缘线以上为梁峁坡(包括峁顶), 沟缘线以下到坡脚线以上为沟坡, 多个坡脚线之间相对平坦的地区为沟谷.结合砖窑沟流域的实际情况于2010年9月进行采样.采样时利用GPS定位, 并记录每个采样点的海拔、经纬度、地貌类型、植被措施等信息.

在流域尺度上采集不同地貌类型、不同治理措施的土壤样本99个, 其中梁峁坡上共采集了55个土壤样本, 沟坡上共采集了32个土壤样本, 沟谷上共采集了12个土壤样本, 林地措施下共采集了19个土壤样本, 草地措施下共采集了36个土壤样本, 农田措施下共采集了44个土壤样本(图 1).每个样本3~5次重复(“S”型路线), 其中0~20 cm土层采集了60个土壤样本, 0~40 cm土层采集了20个土壤样本, 0~100 cm土层采集了19个土壤样本.共计767个土壤样品, 其中0~20 cm的土壤样品大约为60%, 20~40 cm的土壤样品大约为23%, 而40~100 cm的土壤样品大约为17%, 采样土钻直径为3 cm.

修改自文献[29] 图 1 砖窑沟小流域地貌类型和植被措施的样本分布示意 Fig. 1 Sampling sites at the Zhuanyaogou watershed based on vegetation type and landforms

在流域尺度上, 基于地貌类型和植被措施两大因素, 以横断面思路为依据, 采集根系样本, 每个样本5次重复, 采样根钻直径为9 cm.在林地措施下, 在距树干分别大约0.3、0.5、1、1.5和2 m处随机进行根系样品的采集, 在草地措施下, 按照对角线法进行根系样品的采集, 在农田措施下, 2个样品在植株间进行采集, 而3个样品在行间进行采集.在流域尺度上共采集根系样本17个, 其中梁峁坡上采集了8个样本, 沟坡上采集了5个样本, 沟谷上采集了4个样本, 林地措施下采集了5个样本, 草地措施下采集了5个样本, 农田措施下采集了7个样本(图 1), 共计416个根系样品.将采集好的根系样品带回实验室, 人为利用镊子挑选里面的根系, 并根据直径将其分为两类, 细根(d≤2 mm)和粗根(d≥2 mm), 并将其于烘箱中在70℃下烘72 h至恒重.

1.2.2 样品分析

新鲜样品混合均匀后, 在室内进行自然风干, 风干样品过0.15 mm筛后, 测定SOC (H2SO4-K2Cr2O7外加热法)含量(g·kg-1), 土壤容重的测定利用环刀法.

1.2.3 土壤有机碳密度和根系密度的估算

土壤有机碳密度(SOCD)利用下列方程进行计算:

(1)

式中, SOCD是土壤有机碳密度(kg·m-2), SOCi是第i层的土壤有机碳含量(g·kg-1), BDi是第i层的土壤容重(g·cm-3), Di是第i层的土层厚度(m), m是土层的数量.

总土壤有机碳储量(TSOC)由下列方程计算:

(2)

式中, ASOCDi是第i层的土壤有机碳密度(kg·m-2), Si是第i层的面积(m2), Di是第i层的土层厚度(m), m是土层的数量.

根系密度(FRD)由下列方程计算:

(3)

式中, FRDi是第i层的根系密度(g·m-2), Ri是第i层的根系生物量(mg·cm-3), Di是第i层的土层厚度(m), m是土层的数量.

1.3 统计分析

利用Excel软件进行数据的预处理, 利用Sigmplot软件进行有关基础图件的制作, 并用SAS软件的UNIVARIATE plot normal程序包对数据进行正态分布检验.利用SAS软件对不同地貌类型和不同植被措施条件下的SOCD和FRD进行了方差分析(GLM), 当F检验显著时再进行3种地貌类型(梁峁坡、沟坡、沟谷)和3种植被措施(农田措施、林地措施、草地措施)均值间的Duncan检验, 对同一地貌类型条件下不同植被措施以及同一植被措施条件下不同地貌类型进行方差分析(GLM), 当F检验显著时进行均值间的Duncan检验.此外, 利用SAS软件的PROC reg程序包对分析获得的SOCD (y)和FRD (x)进行回归分析(y=alnx+b), 拟合系数a值的大小可以间接反映根系-碳转化效率, 即每增加1单位的FRD其SOCD的增加量.

2 结果与分析 2.1 SOCD和FRD在流域尺度上的影响因素

在砖窑沟小流域内, SOCD (0~100 cm)变化于2.79~4.95 kg·m-2之间(表 1), 其均值为3.48 kg·m-2, 且数据接近正态分布(W=0.91). FRD (0~100 cm)的变化范围为387~1 303 g·m-2, 其均值为713 g·m-2, 数据类似正态分布(W=0.87).同时, 研究结果显示, 除了地貌类型和植被措施交互作用对FRD的影响不显著以外(P>0.05), 植被措施, 地貌类型、土层厚度以及他们之间彼此的交互作用均显著地影响SOCD和FRD的空间变异性(P < 0.05).

表 1 植被措施和土层厚度交互作用对SOCD和FRD的影响1) Table 1 Interaction of vegetation types, landforms, and soil depth with soil organic carbon density (SOCD) and fine root density (FRD) at the watershed scale

2.2 在流域尺度上植被措施对SOCD和FRD的影响

在同一地貌类型条件下, 植被措施显著地影响SOCD和FRD的剖面分布特征, 且二者均呈现出表层(0~20 cm)大于下层(20~100 cm)的趋势(图 2).

图 2 植被措施对SOCD和FRD的影响 Fig. 2 Effect of vegetation type on SOCD and FRD

在梁峁坡上, 表层(0~20 cm)的SOCD呈现出林地措施(1.34 kg·m-2)>草地措施(1.09 kg·m-2)>农田(1.11 kg·m-2)的趋势(P < 0.05), 而表层农田和草地措施的SOCD差异不明显(P>0.05), 在20~100 cm土层, SOCD在农田、林地和草地措施之间的差异不显著(P>0.05);在沟坡上, 不同土层间(0~100 cm土层, 20 cm的间隔)的SOCD均呈现出林地措施(4.95 kg·m-2)>草地措施(3.04 kg·m-2)>农田(2.79 kg·m-2)的趋势; 在沟谷, 不同土层间(0~100 cm土层, 20 cm的间隔)的SOCD在农田和草地措施之间差异不显著, 但是仍呈现出农田大于草地措施的趋势(4.12 kg·m-2和3.19 kg·m-2).

在0~100 cm土层, 在梁峁坡上林地措施的FRD (1 037 g·m-2)分别是草地措施(619 g·m-2)和农田(387 g·m-2) FRD的1.7和2.7倍; 在沟坡上的FRD呈现出林地措施(1 303 g·m-2)>草地措施(756 g·m-2)>农田(472 g·m-2)的趋势; 在沟谷上的FRD也呈现出草地措施大于农田的趋势(571 g·m-2和559 g·m-2).

2.3 在流域尺度上地貌类型对SOCD和FRD的影响

在同一植被措施条件下, 地貌类型对SOCD和FRD的剖面分布特征影响显著(P < 0.05), 且二者均呈现出表层(0~20 cm)大于下层(20~100 cm)的趋势(图 3).

图 3 地貌类型对SOCD和FRD的影响 Fig. 3 Effect of landforms on SOCD and FRD

在农田措施中, 0~40 cm (20 cm间隔)土层, 沟谷的SOCD (2.34 kg·m-2)分别是梁峁坡(1.77 kg·m-2)和沟坡(1.59 kg·m-2)的1.3和1.5倍, 然而沟坡和梁峁坡的SOCD差异不明显(P>0.05), 在40~100 cm土层, 地貌类型对SOCD的影响不显著(P>0.05);在草地措施中, 表层(0~20 cm)的SOCD呈现出沟谷(1.34 kg·m-2)>沟坡(1.11 kg·m-2)>梁峁坡(1.09 kg·m-2)的趋势(P < 0.05), 然而梁峁坡和沟坡的SOCD没有差异(P>0.05), 同时在20~100 cm土层(20 cm间隔), 沟谷、梁峁坡和沟坡的SOCD差异不显著(P>0.05);在林地措施中, 0~60 cm土层(20 cm间隔), 沟坡的SOCD是梁峁坡的1.30倍(3.59 kg·m-2和2.76 kg·m-2), 在60~100 cm土层中, 二者的差异不明显(P>0.05).

在同一植被措施下, 地貌类型对FRD的影响显著(P < 0.05).在0~100 cm土层, 农田措施下的FRD呈现出沟谷(559 g·m-2)>沟坡(472 g·m-2)>梁峁坡(387 g·m-2)的趋势; 草地措施下沟坡(756 g·m-2)和梁峁坡(619 g·m-2)的FRD分别是沟谷(571 g·m-2)的1.3和1.1倍; 在林地措施中, 沟坡的FRD是梁峁坡的1.3倍(1 303 g·m-2和1 037 g·m-2).

2.4 在流域尺度上SOCD和FRD在剖面上的相对分布比例

在小流域尺度上, SOCD和FRD在不同土层厚度间(0~100 cm, 20 cm间隔)的相对比例差异显著(P < 0.05), 且二者均呈现出上层(0~40 cm)远远大于下层(40~100 cm)的趋势(图 4).

图 4 SOCD和FRD在剖面内的相对分布比例 Fig. 4 Distribution ratio of SOCD and FRD in the soil profile(0-100 cm)

在农田措施下, 梁峁坡0~40 cm土层的SOCD和FRD分别占整个剖面(0~100 cm)的57%和67%, 沟坡分别占57%和63%, 沟谷依次占57%和58%;在草地措施下, 梁峁坡0~40 cm土层的SOCD和FRD分别占整个剖面(0~100 cm)的55%和74%, 沟坡分别占58%和60%, 沟谷依次占67%和53%;在林地措施下, 梁峁坡上0~40 cm土层的SOCD和FRD分别占整个剖面(0~100 cm)的61%和62%, 沟谷依次占54%和60%.

在梁峁坡上, 农田措施下0~40 cm土层的SOCD和FRD分别占整个剖面(0~100 cm)的57%和67%, 草地措施下分别为55%和74%, 林地措施下依次为61%和62%;在沟坡上, 农田措施下0~40 cm土层的SOCD和FRD分别占整个剖面(0~100 cm)的57%和63%, 草地措施下分别为58%和60%, 林地措施下依次为54%和60%;在沟谷上, 农田措施下0~40 cm土层的SOCD和FRD分别占整个剖面(0~100 cm)的57%和58%, 草地措施下分别为67%和53%.

2.5 在流域尺度上FRD对SOCD的影响

在流域尺度上, FRD是影响SOCD的重要因素(表 2), 本研究结果显示SOCD (y)随着FRD (x)增加呈现出对数增加的趋势(y=a lnx+b), 且拟合系数a可以指示根系-碳转化效率, 即每增加1单位的FRD其SOCD的增加量.

表 2 不同植被类型和地形条件下的SOCD与FRD的对数回归分析 Table 2 Logarithmic regression analysis of SOCD and FRD for different landforms and vegetation types

在同一植被措施下, 不同地貌类型间的根系-碳转化效率差异显著(P < 0.05).在农田措施下, 其根系-碳转化效率呈现出沟谷(0.87)>沟坡(0.43)=梁峁坡(0.43)的趋势; 在草地措施下, 其根系-碳转化效率呈现出沟坡(0.57)>沟谷(0.45)>梁峁坡(0.27)的趋势; 在林地措施下, 其根系-碳转化效率呈现出梁峁坡(0.56)>沟坡(0.44)的趋势.

在同一地貌类型下, 根系-碳转化效率因不同植被措施差异显著(P < 0.05).在梁峁坡上, 其根系-碳转化效率呈现出林地(0.56)>农田(0.43)>草地(0.27)的趋势, 在沟坡上, 其根系-碳转化效率呈现出草地(0.57)>林地(0.44)>农田(0.43)的趋势, 在沟谷上, 其根系-碳转化效率呈现农田(0.87)>草地(0.45)的趋势.

3 讨论

在砖窑沟小流域的剖面中(0~100 cm), 其TSOC为3.48 kg·m-2, 这显著低于中国热带和亚热带东部的平均TSOC (21.4 kg·m-2)[30]、中国西南山区的平均TSOC (14.43 kg·m-2)[31]、整个中国的平均TSOC (9.6 kg·m-2)[32]、甚至也远远低于整个黄土高原地区的平均TSOC (7.7 kg·m-2)[33].砖窑沟小流域内其较低的TSOC可能归咎于严重的水土流失以及较少的降雨[1, 3, 4].但是, 砖窑沟小流域的TSOC却远远大于纸坊沟小流域的TSOC (3.48 kg·m-2和2.63 kg·m-2)[34], 这可能取决于地形条件, 因为在纸坊沟小流域内并没有考虑地形因素对SOC分布的影响, 而已有的研究结果显示地形显著地影响SOC的空间分布[33, 35].在本研究中, 植被措施和地貌类型对SOCD的影响显著(P < 0.05), 同一地貌类型条件下的SOCD呈现出林地措施>草地措施>农田措施的趋势, 其主要通过影响输入到土壤中有机物质(地上凋落物和地下根系)的数量和质量影响SOCD的时空变异性[19, 20].同一植被措施下的SOCD基本呈现出沟谷>沟坡>梁峁坡的趋势, 其主要通过侵蚀影响SOC的迁移与掩埋[2, 7, 36, 37]和控制土壤水分平衡, 从而间接影响SOCD的时空变异性[8, 38].

根系尤其是细根(d≤2 mm), 对土壤有机碳库的贡献高达25%~80%[39], 且根系在剖面内的分布状态影响着土壤剖面中碳的位置和数量[40].虽然根系凋落物是影响土壤有机碳储量的重要因素, 但是大多数研究仅在相似或者单一地形条件下探讨不同植被条件下根系生物量对SOC的影响[3, 25, 41], 但是同一植被措施下不同地形条件下根系生物量对SOC影响的研究很匮乏.本研究结果显示, SOCD和FRD具有较强的相关性(P < 0.05, 图 4表 2), 且均呈现出上层(0~40 cm)大于下层(40~100 cm)的趋势(图 4).在梁峁坡和沟坡中, 有57%的SOCD分布在上层, 然而却有超过64%的FRD分布在上层, 即植被的根系分布比SOC的分布浅, 这与已有的研究结果相类似[40, 42].但是, 在沟谷中, 有超过62%的SOCD分布在上层, 然而仅有56%的FRD分布在上层, 这可能因为遭受侵蚀且富含有机碳的表土在沟谷进行沉积所造成的[11, 36].进一步的研究结果显示, SOCD随着FRD的增加呈现出对数增加的趋势(表 2), 而在中国不同景观类型下的研究发现SOCD随着FRD的增加呈现出线性增加的趋势[41].除了地下根系凋落物以外, 土壤温度、土壤水分、地上凋落物、地形以及土壤微生物的群落结构等因素均影响SOCD的空间分布[7, 23, 43, 44], 因此对数关系模型更能客观公正地反映FRD对SOCD的贡献, 且拟合系数a值的大小, 能间接地反映根系-碳转化效率, 即每增加单位FRD时SOCD的增加量.

同一地貌类型条件下的根系-碳转化效率的差异与有机物质的数量和性质关系密切, 这与已有的研究结果基本吻合[3, 25, 41].但是, 同一植被措施下, 不同地形条件下根系-碳转化效率的差异可能是地形、人类活动、根系凋落物归还量以及土壤侵蚀综合作用的结果[4, 7, 50, 51].在农田措施下, 其根系-碳转化效率呈现出沟谷(0.87)>沟坡(0.43)=梁峁坡(0.43), 其可能是土壤侵蚀和人类活动综合作用的结果, 一方面随侵蚀迁移和搬运的泥沙和径流分别在沟谷和沟坡进行沉积, 造成该部位的有机碳的富集[4, 45], 另一方面由于沟谷地形平坦, 人类活动(耕作、施肥、灌溉)的加剧间接导致农田措施下的FRD呈现出沟谷(559 g·m-2)>沟坡(472 g·m-2)>梁峁坡(387 g·m-2)的趋势, 进而影响根系-碳转化效率; 在草地措施下, 其根系-碳转化效率呈现出沟坡(0.57)>沟谷(0.45)>梁峁坡(0.27)的趋势, 其可能与植被措施类型有关, 沟坡由于其较大的坡度不利于人类活动, 导致沟坡上的草地以自然恢复的草地为主(本土植被), 而本土植被对环境具有较强的适应性以及抵抗外界干扰的能力[46, 47], 可能不同程度地增加了其根系-碳转化效率; 在林地措施下, 其根系-碳转化效率呈现出梁峁坡(0.56)>沟坡(0.44)的趋势, 这可能与地形限制人类活动和人类对植被措施的选择等因素有关, 而植被措施和人类活动反过来通过影响土壤的物理、化学和生物属性而间接影响SOC的空间分布和储量[33], 进而影响根系-碳转化效率.因此, 从增加根系-碳转化效率的角度而言, 根据对数关系模型结果的进一步分析得到:梁峁坡是恢复林地措施的最佳地形部位, 沟坡是恢复草地措施的最佳地形部位, 而沟谷是农作物生长的最佳地形部位(平坦的地形条件和较高的水分以及SOC储量).

4 结论

在地形条件复杂的地区, 量化根系对SOC的贡献对科学评价水土流失区的土壤碳储量具有重要意义.地形、植被措施和土层厚度及其交互作用显著地影响SOCD和FRD在流域内的空间分布.在流域尺度上, FRD对SOCD的影响显著(P < 0.05), 且SOCD (y)随着FRD (x)增加呈现出对数增加的趋势(y=a lnx+b).在农田措施下, 其根系-碳转化效率呈现出沟谷>沟坡=梁峁坡的趋势; 在草地措施下, 其根系-碳转化效率呈现出沟坡>沟谷>梁峁坡的趋势; 在林地措施下, 其根系-碳转化效率呈现出梁峁坡>沟坡的趋势.因此, 在砖窑沟小流域内, 从增加根系-碳的转化效率角度而言, 沟谷适合进行农业生产, 沟坡适宜进行退耕还草, 而梁峁坡适合进行退耕还林.

参考文献
[1] 唐克丽. 黄土高原水蚀风蚀交错区治理的重要性与紧迫性[J]. 中国水土保持, 2000(11): 11-12, 17.
Tang K L. Importance and urgency of harnessing the interlocked area with both water and wind erosion in the Loess Plateau[J]. Soil and Water Conservation in China, 2000(11): 11-12, 17. DOI:10.3969/j.issn.1000-0941.2000.11.005
[2] Lal R. Soil erosion and the global carbon budget[J]. Environment International, 2003, 29(4): 437-450. DOI:10.1016/S0160-4120(02)00192-7
[3] Wang Z Q, Guo S L, Sun Q Q, et al. Soil organic carbon sequestration potential of artificial and natural vegetation in the hilly regions of Loess Plateau[J]. Ecological Engineering, 2015, 82: 547-554. DOI:10.1016/j.ecoleng.2015.05.031
[4] 张彦军, 郭胜利, 南雅芳, 等. 水土流失治理措施对小流域土壤有机碳和全氮的影响[J]. 生态学报, 2012, 32(18): 5777-5785.
Zhang Y J, Guo S L, Nan Y F, et al. Effects of soil erosion control measures on soil organic carbon and total nitrogen in a small watershed[J]. Acta Ecologica Sinica, 2012, 32(18): 5777-5785.
[5] Lal R. Soil carbon sequestration to mitigate climate change[J]. Geoderma, 2004, 123(1-2): 1-22. DOI:10.1016/j.geoderma.2004.01.032
[6] Pan G X, Smith P, Pan W N. The role of soil organic matter in maintaining the productivity and yield stability of cereals in China[J]. Agriculture, Ecosystems & Environment, 2009, 129(1-3): 344-348.
[7] Dialynas Y G, Bastola S, Bras R L, et al. Topographic variability and the influence of soil erosion on the carbon cycle[J]. Global Biogeochemical Cycles, 2016, 30(5): 644-660. DOI:10.1002/2015GB005302
[8] Griffiths R P, Madritch M D, Swanson A K. The effects of topography on forest soil characteristics in the Oregon Cascade Mountains (USA):implications for the effects of climate change on soil properties[J]. Forest Ecology and Management, 2009, 257(1): 1-7. DOI:10.1016/j.foreco.2008.08.010
[9] McNab W H. A topographic index to quantify the effect of mesoscale landform on site productivity[J]. Canadian Journal of Forest Research, 1993, 23(6): 1100-1107. DOI:10.1139/x93-140
[10] Kobayashi Y, Koike F. Separating the effects of land-use history and topography on the distribution of woody plant populations in a traditional rural landscape in Japan[J]. Landscape and Urban Planning, 2010, 95(1-2): 34-45. DOI:10.1016/j.landurbplan.2009.11.003
[11] Wang Y F, Fu B J, Lü Y H, et al. Effects of vegetation restoration on soil organic carbon sequestration at multiple scales in semi-arid Loess Plateau, China[J]. CATENA, 2011, 85(1): 58-66. DOI:10.1016/j.catena.2010.12.003
[12] Wang Z, Liu G B, Xu M X, et al. Temporal and spatial variations in soil organic carbon sequestration following revegetation in the hilly Loess Plateau, China[J]. CATENA, 2012, 99: 26-33. DOI:10.1016/j.catena.2012.07.003
[13] Gong J, Chen L, Fu B, et al. Effect of land use on soil nutrients in the loess hilly area of the Loess Plateau, China[J]. Land Degradation & Development, 2006, 17(5): 453-465.
[14] Chaplot V, Bouahom B, Valentin C. Soil organic carbon stocks in Laos:spatial variations and controlling factors[J]. Global Change Biology, 2010, 16(4): 1380-1393. DOI:10.1111/gcb.2010.16.issue-4
[15] Fang X, Xue Z J, Li B C, et al. Soil organic carbon distribution in relation to land use and its storage in a small watershed of the Loess Plateau, China[J]. CATENA, 2012, 88(1): 6-13. DOI:10.1016/j.catena.2011.07.012
[16] Oueslati I, Allamano P, Bonifacio E, et al. Vegetation and topographic control on spatial variability of soil organic carbon[J]. Pedosphere, 2013, 23(1): 48-58. DOI:10.1016/S1002-0160(12)60079-4
[17] Deng L, Liu G B, Shangguan Z P. Land-use conversion and changing soil carbon stocks in China's 'Grain-for-Green' Program:a synthesis[J]. Global Change Biology, 2014, 20(11): 3544-3556. DOI:10.1111/gcb.12508
[18] Song X Z, Peng C H, Zhou G M, et al. Chinese Grain for Green Program led to highly increased soil organic carbon levels:a meta-analysis[J]. Scientific Reports, 2014, 4: 4460.
[19] Kätterer T, Bolinder M A, Andrén O, et al. Roots contribute more to refractory soil organic matter than above-ground crop residues, as revealed by a long-term field experiment[J]. Agriculture, Ecosystems & Environment, 2011, 141(1-2): 184-192.
[20] Rasse D P, Rumpel C, Dignac M F. Is soil carbon mostly root carbon? Mechanisms for a specific stabilisation[J]. Plant and Soil, 2005, 269(1-2): 341-356. DOI:10.1007/s11104-004-0907-y
[21] De Graaff M A, Schadt C W, Rula K, et al. Elevated CO2 and plant species diversity interact to slow root decomposition[J]. Soil Biology and Biochemistry, 2011, 43(11): 2347-2354. DOI:10.1016/j.soilbio.2011.07.006
[22] Yuan Z Y, Chen H Y H. Fine root biomass, production, turnover rates, and nutrient contents in boreal forest ecosystems in relation to species, climate, fertility, and stand age:literature review and meta-analyses[J]. Critical Reviews in Plant Sciences, 2010, 29(4): 204-221. DOI:10.1080/07352689.2010.483579
[23] Hu Y L, Zeng D H, Ma X Q, et al. Root rather than leaf litter input drives soil carbon sequestration after afforestation on a marginal cropland[J]. Forest Ecology and Management, 2016, 362: 38-45. DOI:10.1016/j.foreco.2015.11.048
[24] Wang Y, Ji H F, Wang R, et al. Impact of root diversity upon coupling between soil C and N accumulation and bacterial community dynamics and activity:result of a 30 year rotation experiment[J]. Geoderma, 2017, 292: 87-95. DOI:10.1016/j.geoderma.2017.01.014
[25] 王彦丽. 不同植被恢复措施下剖面根系与SOC的分布特征[J]. 水土保持研究, 2013, 20(6): 19-23, 43.
Wang Y L. Profile distribution characteristics of root biomass and soil organic carbon under different vegetation restoration measures[J]. Research of Soil and Water Conservation, 2013, 20(6): 19-23, 43.
[26] 韩畅, 宋敏, 杜虎, 等. 广西不同林龄杉木、马尾松人工林根系生物量及碳储量特征[J]. 生态学报, 2017, 37(7): 2282-2289.
Han C, Song M, Du H, et al. Biomass and carbon storage in roots of Cunninghamia lanceolata and Pinus massoniana plantations at different stand ages in Guangxi[J]. Acta Ecologica Sinica, 2017, 37(7): 2282-2289.
[27] 匡文浓, 钱建强, 马群, 等. 五种荒漠灌木群落土壤有机碳垂直分布及其与根系分布的关系[J]. 生态学杂志, 2016, 35(2): 275-281.
Kuang W N, Qian J Q, Ma Q, et al. Vertical distribution of soil organic carbon and its relation to root distribution in five desert shrub communities[J]. Chinese Journal of Ecology, 2016, 35(2): 275-281.
[28] 罗永清, 赵学勇, 王涛, 等. 沙地植物根系特征及其与土壤有机碳和总氮的关系[J]. 草业学报, 2017, 26(8): 200-206.
Luo Y Q, Zhao X Y, Wang T, et al. Characteristics of the plant-root system and its relationships with soil organic carbon and total nitrogen in a degraded sandy grassland[J]. Acta Prataculturae Sinica, 2017, 26(8): 200-206.
[29] 中国科学院水利部西北水土保持研究所. 黄土高原综合治理试验示范区:专题地图集[M]. 北京: 测绘出版社, 1991.
[30] Li Z, Jiang X, Pan X Z, et al. Organic carbon storage in soils of tropical and subtropical China[J]. Water, Air, and Soil Pollution, 2001, 129(1-4): 45-60.
[31] Zhang Y, Zhao Y C, Shi X Z, et al. Variation of soil organic carbon estimates in mountain regions:a case study from Southwest China[J]. Geoderma, 2008, 146(3-4): 449-456. DOI:10.1016/j.geoderma.2008.06.015
[32] 于东升, 史学正, 孙维侠, 等. 基于1:100万土壤数据库的中国土壤有机碳密度及储量研究[J]. 应用生态学报, 2005, 16(12): 2279-2283.
Yu D S, Shi X Z, Sun W X, et al. Estimation of China soil organic carbon storage and density based on 1:1000000 soil database[J]. Chinese Journal of Applied Ecology, 2005, 16(12): 2279-2283. DOI:10.3321/j.issn:1001-9332.2005.12.011
[33] Liu Z P, Shao M A, Wang Y Q. Effect of environmental factors on regional soil organic carbon stocks across the Loess Plateau region, China[J]. Agriculture, Ecosystems & Environment, 2011, 142(3-4): 184-194.
[34] Li M M, Zhang X C, Pang G W, et al. The estimation of soil organic carbon distribution and storage in a small catchment area of the Loess Plateau[J]. CATENA, 2013, 101: 11-16. DOI:10.1016/j.catena.2012.09.012
[35] Zhu H H, Wu J S, Guo S L, et al. Land use and topographic position control soil organic C and N accumulation in eroded hilly watershed of the Loess Plateau[J]. CATENA, 2014, 120: 64-72. DOI:10.1016/j.catena.2014.04.007
[36] 杜兰兰, 王志齐, 王蕊, 等. 模拟条件下侵蚀-沉积部位土壤CO2通量变化及其影响因素[J]. 环境科学, 2016, 37(9): 3616-3624.
Du L L, Wang Z Q, Wang R, et al. Variation of soil CO2 flux and environmental factors across erosion-deposition sites under simulation experiment[J]. Environmental Science, 2016, 37(9): 3616-3624.
[37] 李如剑, 张彦军, 赵慢, 等. 坡度和降雨影响土壤CO2通量和有机碳流失的模拟研究[J]. 环境科学学报, 2016, 36(4): 1336-1342.
Li R J, Zhang Y J, Zhao M, et al. Simulation on the effects of slope and rainfall on soil CO2flux and SOC loss[J]. Acta Scientiae Circumstantiae, 2016, 36(4): 1336-1342.
[38] Zhao X, Wu P, Gao X, et al. Soil quality indicators in relation to land use and topography in a small catchment on the Loess Plateau of China[J]. Land Degradation & Development, 2015, 26(1): 54-61.
[39] Usman S, Singh S P, Rawat Y S, et al. Fine root decomposition and nitrogen mineralisation patterns in Quercus leucotrichophora and Pinus roxburghii forests in central Himalaya[J]. Forest Ecology and Management, 2000, 131(1-3): 191-199. DOI:10.1016/S0378-1127(99)00213-3
[40] Jobbágy E G, Jackson R B. The vertical distribution of soil organic carbon and its relation to climate and vegetation[J]. Ecological Applications, 2000, 10(2): 423-436. DOI:10.1890/1051-0761(2000)010[0423:TVDOSO]2.0.CO;2
[41] Wang Y G, Li Y, Ye X H, et al. Profile storage of organic/inorganic carbon in soil:from forest to desert[J]. Science of the Total Environment, 2010, 408(8): 1925-1931. DOI:10.1016/j.scitotenv.2010.01.015
[42] Olupot G, Daniel H, Lockwood P, et al. Root contributions to long-term storage of soil organic carbon: theories, mechanisms and gaps[A]. In: Proceedings of the 19th World Congress of Soil Science, Soil Solutions for a Changing World[C]. Brisbane, Australia: DVD, 2010. 112-115.
[43] Dai X Q, Wang H M, Fu X L. Soil microbial community composition and its role in carbon mineralization in long-term fertilization paddy soils[J]. Science of the Total Environment, 2017, 580: 556-563. DOI:10.1016/j.scitotenv.2016.11.212
[44] 邹俊亮, 郭胜利, 李泽, 等. 小流域土壤有机碳的分布和积累及土壤水分的影响[J]. 自然资源学报, 2012, 27(3): 430-439.
Zou J L, Guo S L, Li Z, et al. Soil organic carbon accumulation and distribution in a small watershed and the effects of soil water[J]. Jounnal of Natural Resources, 2012, 27(3): 430-439.
[45] Khan F, Hayat Z, Ahmad W, et al. Effect of slope position on physico-chemical properties of eroded soil[J]. Soil & Environment, 2013, 32(1): 22-28.
[46] Hyvönen T. Impact of temperature and germination time on the success of a C4 weed in a C3 crop:Amaranthus retroflexus and spring barley[J]. Agricultural and Food Science, 2008, 20(2): 183-190.
[47] 南志标, 王锁民, 王彦荣, 等. 我国北方草地6种乡土植物抗逆机理与应用[J]. 科学通报, 2016, 61(2): 239-249.
Nan Z B, Wang S M, Wang Y R, et al. Stress tolerance mechanisms of 6 native plant species growing in China's northern grassland and their utilization[J]. Chinese Science Bulletin, 2016, 61(2): 239-249.