土壤有机碳（SOC）是重要的土壤属性之一，在陆地生态系统功能中起着重要作用^[1]，SOC含量大小对土壤性质以及质量起着至关重要的影响，是评价土壤肥力的重要指标^[2]。陆地生态系统中的碳储量大约是大气中CO₂的2倍多，陆地植被碳库的3倍^[3]。全球0.3 m土壤深度的碳储存量为684~724 Pg，1 m土壤深度碳储量为1 462~1 548 Pg^[4]。大量的有机碳被储存在土壤中，即使是土壤碳含量的微小变化也会显著改变大气中CO₂的浓度，并造成全球气候变化^[5]。因此，学者们在全球^[6]、国家^[7]和区域^[8]尺度上对SOC分布以及储量进行了大量的研究，然而这些研究由于数据来源的不同和SOC固有的空间变异性，导致结果的不确定性。SOC受到诸多因素的影响，包括气候，地形、土壤质地、土地利用等，因此在给定区域建立一个详细的数据库，应该考虑到以上因素。黄土高原地区地形破碎，土地利用类型复杂多样，在区域尺度上，形成了典型的黄土高原沟壑和黄土丘陵沟壑等流域侵蚀地貌特征，成为世界上水土流失最严重的地区之一^[9]。受水土流失的影响，SOC在黄土丘陵沟壑的梁峁顶、梁峁坡、沟谷地等呈现不同规律的迁移^[10-11]。经过长期的植被恢复与重建，黄土高原土地利用和植被发生了明显的变化。从另一个角度而言，地形支配着水、热资源的分配，从而影响植被和土地利用方式在空间上的配置，对SOC的输入产生作用^[12]。因此在黄土高原地区，地形和土地利用对SOC库的响应最为敏感^[13]，而在不同流域尺度SOC的研究^[14-15]中，地形与土地利用方式共同作用下有机碳变化研究的报道较少，并且已有的研究大多集中在传统统计和地理统计上^[16]，不同土层尤其是深层（> 20 cm）SOC受地形及土地利用因子共同作用的贡献率计算更是鲜见。

在本文中，选择在黄土高原典型流域-王茂沟流域开展相关工作，对研究区SOC的分布和储量进行研究。通过分析典型小流域0~100 cm土层SOC的分布特征，分析地形、土地利用对其分布的影响，基于方差探讨黄土高原小流域地形、植被及其交互作用对SOC的影响贡献率，从小流域尺度揭示各因子对SOC垂直分布的影响规律，为黄土高原地区植被恢复的土壤碳效应评估和生态规划提供理论依据，对科学认识黄土高原土壤在全球变化的作用与响应机制具有重要意义。

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

研究区选择陕西省绥德县王茂沟流域为研究区，其位于陕西省绥德县东南5 km处无定河流域（图 1a），是黄河中游黄土丘陵区第一副区具有代表性的流域之一（110°20′~110°22′E，37°34′~37°36′N）。王茂沟流域面积5.97 km²，主沟长3.75 km，沟道平均比降2.7，沟壑密度4.3 km·km^-2，海拔936~1 188 m，该流域的平均坡度为2.7%。流域内高地和沟壑土地面积为2.97 km²，占总面积的46.7%（图 1b）。流域属大陆性季风气候，年平均气温10.2 ℃，多年平均降雨量513.1 mm，主要集中在7—9月，汛期降水量占全年降水量的73.1%，且多以暴雨形式出现。土壤类型均为黄绵土，水土流失严重，流域年输沙量可达74.13 Mg·hm^-2·a^-1。作为黄土高原水土流失的典型流域，王茂沟流域主要以水力侵蚀为主。

从20世纪50年代开始，为改善黄土丘陵区严重的水土流失，黄河水利委员会绥德水土保持科学实验站将王茂沟流域作为试验性治理小流域之一。直至90年代，生物措施（退耕还林还草）及工程措施（梯田，淤地坝）在流域内大面积开展。其中修筑梯田112.47 hm²，植树造林199.96 hm²，恢复草地27.25 hm²，淤积坝地25.16 hm²（图 1c）。人工草被主要为紫花苜蓿（Medicago sativa L.）、赖草（Leymus secalinus T.）、大针茅（Stipa grandis P.）、阿尔泰狗娃花（Heteropappus altaicus N.）、长芒草（Stipa bungeana T.）等，人工林地主要为杏树（Armeniaca vulgaris L.）、刺槐（Robinia pseudoacacia L.）、枣树（Ziziphus jujube M.）等，坡耕地主要种植大豆（Glycine max M.）、谷子（Setaria italica）、马铃薯（Solanum tuberosum L.）等，灌木植被主要为柠条（Caragana Korshinskii M.）、沙棘（Hippophae rhamnoides L.）、酸枣（Ziziphus jujube M.）等，梯田主要植被类型为苹果（Malus domestica）、杏树（Armeniaca vulgaris L.）、侧柏（Platycladus orientalis L.）等，坝地主要种植玉米（Zea mays）。水土流失治理面积为367.81 hm²，侵蚀控制率为61.61%^[17]。

1.2 样品采集

基于小流域内地形和土地利用两大因素，采用“分层采样的方法”收集土壤样品。将地形（坡顶、坡上、坡中、坡下）作为一级层次，将不同地形部位土地利用方式（坡耕地、林地、草地、灌木地、梯田）作为二级层次，对土壤样品进行系统分析。在2014年7月，在研究区进行网格采样，网格间距为150 m，结合王茂沟流域实际地形后，样点分布情况见图 1c。使用直径为6.9 cm手钻进行土壤样品采集，共布设采样点位216个，每个点位采样深度为100 cm，分为5层（0~20 cm、20~40 cm、40~60 cm、60~80 cm、80~100 cm）。并对应在采样原位各层使用环刀采集土样，测量土壤容重，用于有机碳储量计算。在采集土壤样品同时，记录采样点土地利用类型、坡度和坡向等相关信息。其中，坡耕地土壤样品75个，林地55个，草地325个，灌木地90个，梯田225个，共得到分析样品770个，环刀样770个。

1.3 分析方法

有机碳采用德国耶拿公司生产的有机碳分析仪multi N/C^®3100测定，取研磨过100目筛的土样1.000 g，用1 mol·L^-1的盐酸溶液浸泡24 h，去除土壤中的无机碳，用TOC分析仪测定土壤中的有机碳含量；土壤容重采用环刀法测定（表 1）。

1.4 数据处理

经检验研究区的土壤有机碳数据满足平稳假设，因此本文采用ArcGIS 10.2中的普通Kriging插值方法进行有机碳储量的空间插值，重点考虑区域化变量数据的空间相关性，并以GS+软件拟合的半变异函数进行插值。估算某测量点z值的通用方程为^[18]：

$ Z_{0}=\sum_{i=1}^{s} Z_{X} W_{X} $

(1)

式中，Z₀为待估计值，Z_X为已知点的x值，W_X为与已知点关联的权重，s为用于估算的已知点的数目，权重由数据空间特性拟合的半变异函数得到。

有机碳储量可以通过以下等式计算^[19]：

$ D_{i}=\rho_{i} \cdot h_{i} \cdot c_{i} $

(2)

$ \mathrm{SOC}_{i}=D_{i} \cdot A_{i} $

(3)

式中，ρ_i为土壤容重（g·cm^-3），h_i为土层厚度（cm），c_i为有机碳含量（g·kg^-1），SOC_i为有机碳储量（kg），A_i代表每个土地利用面积（m²），D_i为有机碳密度（kg·m^-2）。

根据野外实测数据与实验室的分析结果，采用一般线性模型（GLM）分析了土地利用、地形、深度以及它们的交互作用对SOC的影响。利用SPSS18.0对SOC含量进行方差分析（ANOVA）。

2 结果 2.1 地形对剖面SOC分布特征的影响

研究区0~100 cm土层SOC含量在3.13~4.72 g·kg^-1之间，且以表层SOC含量最高。地形显著影响了剖面SOC的空间分布（表 2），表现为坡上（4.49 g·kg^-1） > 坡中（4.30 g·kg^-1） > 坡下（3.97 g·kg^-1） > 坡顶（3.34 g·kg^-1）。100 cm土层中，坡上、坡中、坡下SOC含量分别较坡顶高34.6%、28.9%和18.8%（P < 0.05）。坡顶、坡上和坡中位置，0~20 cm土层SOC含量均显著高于20~100 cm各土层SOC含量，同一坡位，20 cm以下土层SOC含量差异不显著；而在该流域坡下部位，包括表层在内，各层SOC含量之间均无显著性差异。

2.2 土地利用方式对剖面SOC分布特征的影响

土地利用方式也显著影响剖面SOC含量空间分布特征（表 3）。100 cm土层中，坡耕地SOC含量最低，仅为3.47 g·kg^-1（P < 0.05）。林地（4.31 g·kg^-1）、梯田（4.25 g·kg^-1）、草地（4.12 g·kg^-1）和灌木（3.82 g·kg^-1）分别较坡耕地SOC含量增加24.2%、22.4%、18.7%和10.1%。0~20 cm土层林地SOC含量最高（P < 0.05）。SOC含量在60 cm出现临界层，在60 cm土层以上，坡耕地SOC含量最低，其他各土地利用类型之间SOC含量均没有显著性差异；在60 cm土层以下，仅梯田与坡耕地SOC含量表现出显著性差异（P < 0.05）。且对于同一土地利用而言，林地、草地、灌木表层（0~20 cm）SOC与其各自深层土壤SOC含量表现出显著差异。

2.3 地形与土地利用交互作用对剖面SOC分布特征的影响

地形和土地利用方式交互作用极显著（P < 0.01）影响土壤剖面SOC分布特征（图 2）。坡顶草地（3.19 g·kg^-1）、灌木（3.60 g·kg^-1）、梯田（3.33 g·kg^-1）SOC平均含量分别为坡耕地（2.93 g·kg^-1）的1.09倍、1.23倍、1.14倍，坡顶0~20 cm土层SOC含量均显著高于坡耕地，草地、灌木、梯田0~20 cm土层SOC含量均高于其他各土层SOC含量（P < 0.05）。除坡耕地外，其余各土地利用0~20 cm土层SOC含量均显著高于其他土层SOC含量（P < 0.05），20~40 cm土层林地、草地、灌木SOC含量均显著高于坡耕地，其余各土层间SOC含量分布较为均一。坡中位置，坡耕地在0~60 cm土层SOC含量最低，林地、草地、梯田0~20 cm土层SOC含量均显著高于各20~100 cm土层SOC含量（P < 0.05）。在坡下位置，坡耕地SOC含量在0~80 cm土层显著低于草地和灌木，草地和灌木0~20 cm土层SOC含量最高（P < 0.05），坡耕地各土层之间SOC含量差异不显著。同一地形条件下生态建设措施的不同，可以导致进入土壤的植物残体数量和性质各异，且由于人为干扰程度的不同，影响土壤碳素矿化、运输以及植物的吸收利用，最后导致土壤中SOC含量存在差异。在不同的地形条件下，坡耕地0~20 cm土层SOC含量均表现为最低，且各土层之间SOC含量变化也较为一致，与Zhao等^[20]研究结论一致。将坡耕地恢复为林地、草地、梯田或者灌木地可以显著增加SOC的含量，对于区域而言，进行生态恢复可增加碳截存量，从而加强土壤碳汇的能力及区域碳储量。

2.4 剖面SOC含量的影响因子

根据一般线性模型（GLM）的方差成分分析，对土地利用、地形、深度及其交互作用对SOC含量变异性的贡献进行了计算。各因子以及土地利用与地形的交互作用对流域剖面SOC含量空间分布有极显著影响（P < 0.01）。土地利用与深度的交互作用、地形与深度的交互作用以及土地利用、地形、深度三者的交互作用对剖面SOC含量分布影响不显著。7种不同角度的对比分析发现，地形对SOC含量影响最为显著，可解释32.50%的变异性。其次为土地利用与地形的交互作用，可解释7.40%的变异性，地形对其贡献率较其他影响因子高出3倍甚至更高。单因子的贡献率中，土层深度对SOC含量的影响最小，解释其5%的变异性。土地利用与深度对SOC的交互作用影响最小，仅解释其0.80%的变异性。由此可知，地形是本研究区SOC的重要影响因子。

2.5 区域SOC空间分布特征

以上文数据和分析为基础，通过GS+软件对王茂沟流域不同土层SOC进行半方差函数分析（表 4），研究土壤特性空间分布的方向性和相关距离，根据软件得出的决定系数和残差平方和选择最优插值模型并对流域未采集部位SOC进行最优内插估值和成图。对Kriging插值模型结果进行正交检验，参考最优模型评价标准：标准平均值最接近于0，均方根预测误差最小，平均标准误差最接近于均方根预测误差，标准均方根预测误差最接近于1，插值结果精度越高。因此，本文采用普通Kriging插值模型对王茂沟流域土壤有机碳含量的插值结果满足精度要求。

将流域内各层SOC含量通过K-S正态检验，通过比较检测P值，均与正态性没有显著差异（P > 0.05），流域SOC符合正态分布，因此根据半方差分析所得出的参数，对各土层的SOC进行普通Kriging插值（图 3），空间分布基本呈现相似的格局。总体而言，王茂沟流域SOC集中呈斑块分布。随着深度的增加，流域SOC的空间分布向均一性发展，这是因为土壤表层与外界环境因子相互作用更加密切，异质性更强。SOC高的地区对应于林地、灌木地、草地和梯田，与实际调查相符。这些地区由于管理良好，土壤流失减小，同时林地、灌木地、梯田凋落物等外源出入碳聚集在土壤表层，可为微生物活动提供充足的能量和物质来源，从而会促进土壤中生物的活动。草地地上部分的枯枝落叶可以为土壤提供大量的碳源。因此，这些地区具有良好的保肥保水作用。

3 讨论 3.1 地形和土地利用对SOC的影响

复杂的地形因子对土壤的物质、能量过程、理化性质、土地生产力、土壤水分状况、微气候等均会产生影响^[21]。坡位、坡向、海拔等相关地形因素控制着小尺度区域尤其是坡面土壤的水、肥、气、热等。千沟万壑的黄土高原，地形破碎、复杂多变，起伏的地形地貌景观对SOC有重要影响。坡位不同，土壤有机质空间分布特征及凋落物腐殖化、矿化过程也不尽相同，与本研究结果相同。在黄土高原腹地延河流域，Sigua等^[22]也发现沟坡（即坡面）SOC含量较峁顶高。坡（峁）顶土地利用类型相对单一，农作物和果树是其主要的植被类型，由于人为翻耕，土壤的通透性增加，土壤质地松软，水蚀作用强烈，水土流失严重，土壤中SOC也随之流失；作物收获后其植物体内的碳不能归还至土壤中，导致坡耕地不利于SOC的积累和保护^[23]。坡上、坡中位置植被类型最为完整且具有多样性，凋落物等外源输入碳聚集在土壤表层，同时为微生物活动提供充足的能量和物质来源，从而促进了土壤表层的生物活动^[24]，表层较多的SOC输入，以及发达的植被细根分布，使其控制土壤侵蚀和水土流失能力较强^[25]。尽管坡下有来自其他坡位侵蚀泥沙的碳输入，但经过坡上和坡中植被的拦截作用，到达坡下的泥沙所携带的SOC含量较低，因此坡下SOC输入较小，使得坡下SOC含量相对较低。正是由于不同地形部位对SOC含量影响的不同，使得地形因子对SOC变异性的解释度最高。

针对土地利用方式的结构调整，不仅可以明显控制和减缓土壤侵蚀的发生，而且是提高陆地生态系统SOC的有效措施。在耕作过程中，土壤在耕层内上下扰动，导致本研究区坡耕地SOC含量在垂直方向上未呈现出表层高于底层的现象。与坡耕地相比，梯田种植物种多样性较高，枯落物层、根系发达程度均较坡耕地高，导致梯田整体SOC含量高于坡耕地。0~60 cm土层是SOC在土壤剖面中分布变化最为活跃的层次，尤其表现在退耕还林（草）之后的土地利用类型中；60~100 cm土层是SOC分布相对稳定的区域^[26]，这与本文得出的流域SOC空间分布结果一致。林地凋落物在土壤表层大量积累和分解，致使0~20 cm土层SOC急剧增加。研究结果显示将坡耕地恢复为灌木、林地、草地、梯田后，均可以增加0~60 cm土层SOC含量。总体来看坡耕地在退耕还林或撂荒后，由于枯落物输入、微生物活动等影响，表层土壤有机碳含量明显增加。而梯田作为我国较常用的水土保持耕作措施，改善了农业生产条件，增加土壤通透性，提高了水分入渗量和蓄存量，从而影响土壤水、肥、气、热状况^[27]。

3.2 各因子及其交互作用对SOC变异的贡献

分布在0~100 cm深度的SOC，极易受人为干扰和气候变化的影响，同时土地利用^[28]、地形^[29]、土壤类型、成土母质、管理措施^[30]等多种因子也会使得SOC分布及空间异质性发生改变。我国黄土高原地区千沟万壑，支离破碎，土地利用类型复杂多样，形成了典型的沟壑侵蚀地貌（梁峁顶、梁峁坡及沟谷地）。因此该区有机碳的分布也往往是多种因子综合作用的结果，针对某一具体地区而言，SOC的主控因子也就会因时间、空间和尺度不同而不同。赵明月等^[31]对陕西安塞县典型小流域的有机碳研究发现，该区地形因子对有机碳含量和密度的影响高于土地利用；而在黄河腹地延河流域^[32]，植被类型对SOC的影响大于地形因子。在本研究流域，地形因子对SOC含量影响最大，其中又以坡位对SOC含量的影响最大^[33]。黄土高原梁峁顶土地利用方式相对单一，由于常年遭受暴雨冲蚀，且植被盖度较低，侵蚀剧烈，在坡位较高、单一土地利用方式两种因素的共同影响下，其SOC含量最低。峁坡上广泛分布梯田、林地、草地及果园等，植被盖度较高，有效减缓了水土-养分流失的发生，使得峁坡面SOC含量高于坡顶。因此，与其他因子的交互作用相比，土地利用与地形的交互作用对本研究流域SOC含量变异性贡献最大（7.4%）。此外，分布在坡顶的坡耕地表层土壤结构疏松，植被盖度较低，且遭受雨滴击溅、冲刷程度更大，因而较其他坡位、土地利用更易发生侵蚀，伴随土壤侵蚀的表层有机碳也随之流失，从而表现出坡顶表层SOC含量最低，随着土层变深，SOC含量逐渐增加的规律。林草地有机碳的输入多以表层枯落物为主，加之梯田肥料施入，且侵蚀相对坡顶较弱，因此坡上、坡中、坡下各土地利用方式下均以土壤表层SOC含量最高。较深层土壤水分匮缺，通气性较好，微生物相对活跃，促进了SOC的分解，使得深层土壤SOC含量较低。不同土地利用方式下SOC含量的垂直变化规律较为一致，因而土地利用与土层深度的交互作用对SOC含量变异的解释度最低。除梯田外，其他土地利用方式深层（> 20 cm）SOC含量显著小于表层（0~20 cm）土壤，说明人为活动对土壤干扰程度明显，经过人为精耕细作与水土保持措施的有效发挥，SOC可以有效被固存。

综上所述，为了增加研究区碳储量，不仅要合理规划不同地貌部位的土地利用，还应减少坡耕地面积及人为干扰等不利因素，重视水土保持措施生态、社会、经济效益，降低区域内SOC的流失，增加区域土壤“碳汇”能力。

4 结论

黄土丘陵第一副区王茂沟小流域在0~100 cm土层中，SOC平均含量坡上、坡中含量最高，其次为坡下，坡顶SOC含量最低；坡耕地SOC含量最低。林地、梯田、草地和灌木分别较坡耕地SOC含量增加24.2%、22.4%、18.7%和10.1%。表层SOC更易受到环境因子的影响，梯田等水土保持措施可明显固存深层SOC。方差成分估计表明，土地利用、地形、深度以及土地利用与地形的交互作用对流域剖面SOC含量空间分布有着极显著的影响（P < 0.01），其中地形对SOC含量的贡献率最高，达到32.50%。土地利用与地形的交互作用在各因子的交互作用中表现出对SOC含量变异解释度最高。流域SOC在空间上呈斑块分布，随着深度的增加，流域SOC的空间分布向均一性发展。