土壤碳库是陆地生态系统中最大而又最活跃的有机碳库，也是全球碳循环的重要组成部分，对全球气候变化有着重要的影响^[1]，其微小变化都将显著地影响着碳汇或碳源的大小和分布^[2]，增加土壤碳库储量可以减缓由于人类活动导致的大气CO₂的增加，相反，土壤中碳的释放将加速大气CO₂的升高^[3]，它是土壤质量评价和土地可持续利用管理中必须考虑的重要指标^[4]。

土地利用变化是影响陆地生态系统碳循环的主要因素之一，也是仅次于石化燃料燃烧而使大气CO₂浓度急剧增加的主要人为活动^[5]。土地利用变化发生后，由于影响土壤有机碳(SOC)的动态过程和相关性质发生了变化，使SOC的储量、分布等也发生了变化，这些变化又将进一步影响土壤向大气释放CO₂的强度与过程，因此土地利用变化对SOC储量、分布的影响，是对SOC影响的本质^[6]。使土壤碳库贮量积累的过程是碳汇，而使土壤碳贮量减少的过程是碳源。准确认识土地利用变化过程对土壤碳贮量的影响是增加还是减少，是直接影响还是间接影响，对于准确认识和评价陆地生态系统的碳源/汇功能具有重大意义^[7]。

SOC的时空变异性是指在土壤质地相同的区域内，SOC在同一时刻、各个空间位置上或在不同时刻、同一空间位置上量值的差异性，其变化情况不仅反映SOC变化趋势，印证了区域生态演变规律和生态建设的基本成果，也是国际全球气候变化问题研究的核心内容之一^[2]。目前对于SOC变化的研究多集中于同一时刻的不同土地利用类型之间的比较分析^[8]，而不同时刻同一地区土地利用类型转换后的土壤碳效应差异研究还有待深入。近年来，国内外学者的研究也表明合理的土地利用方式变化可以增加SOC含量并影响其时空动态变化^{[7, 8, 9]}，但在黄土高原小流域尺度上的相关试验研究较少。黄土高原地区地形破碎，土地利用类型复杂多样，也是我国SOC储量最低的地区之一^{[9, 10]}。20世纪90年代以来，该地区进行了大规模的生态环境修复工作，尤其是退耕还林工程，通过植被恢复引起了土地利用的变化，并改变了土壤内部及与植物群落的养分循环，必然引起土壤碳库变化^[11]。由于大量的SOC集中于20 cm以内的表层中，在不同的土地利用方式作用下，其空间变异较深层(20 cm以下)的影响更大^{[12, 13]}。因此，本文以黄土丘陵沟壑区的羊圈沟小流域为研究对象，结合碳“汇”、“源” 的概念，采用野外调查、室内分析并与地统计学相结合的方法，对小流域尺度在2006—2011年土地利用变化影响表层土(0—20 cm)有机碳密度、储量及其空间分布格局变化进行综合评价，以期为黄土高原植被恢复与生态建设提供参考，也可以为大尺度的研究提供依据。

研究区域为位于延安市宝塔区东北方向14 km处李渠镇的羊圈沟流域(36°42′N，109°31′E)，总面积为2.02 km²，为延河左岸的二级支沟，碾庄沟流域的一级支流。区域地貌类型为黄土梁和黄土沟，沟谷密度为2.74 km/km²，属典型黄土丘陵沟壑区。区域气候为半干旱大陆性季风气候，多年平均降水量535 mm，降水多集中在7—9月，年际变率大。流域内土壤以黄绵土为主，抗蚀性差，水土流失严重。近年来流域内实施了退耕还林工程，经过十几年的流域治理，农业景观发生了很大的变化。流域植被在区划上属于森林草原过渡带，流域内由于人为活动的干扰，自然植被破坏殆尽，多为人工种植而形成的次生植被，植物种类主要有刺槐(Robinia pseudoacacia)、柳树(Salix spp.)和杨树(Poplar spp.)等。

分别于2006年和2011年在黄土丘陵区羊圈沟流域选取6种具有代表性的土地利用方式进行采样：草地、灌木林地、乔木林地、梯田、坡耕地和坝地(水域和建设用地无采样点)。两次采样方法相同，均按土地利用类型和土壤断面，同时考虑微地貌布设采样点，采用地形剖面线法进行坡面采样，采用5点混合的方式采集混合样。采样点为114个，按照土地利用类型面积比例分配采样点个数，以保证样本对不同土地利用类型的代表性，面积比例较大的地类如乔木林地和灌木林地样点分布较多，而面积比例较小的地类，如坝地和坡耕地样点分布较少。取样器为荷兰Eijkekamp公司生产的直径为6 cm的半圆凿型土钻。采样深度为0—20 cm，土壤样点的具体分配和分布情况如表 1、图 1所示。野外采样时采用GPS定点，确定每个采样点的经纬度坐标，同时详细记录每个采样点的坡度、坡向及其土地利用类型等信息。

图 1 研究区域样点分布图 Fig.1 Distribution map of sample sites in study area

图选项

(1)粒径测定    取土壤样品过2 mm筛后测定土壤粒径，分析采用Mastersizer 2000激光粒度仪进行，测定前处理主要是去除土壤有机质和碳酸盐，上机前用六偏磷酸钠溶液分散土壤溶液，测量范围为0.02—2000 μm，重复测量误差小于±3%。

(2)SOC含量测定    取土壤样品过0.28 mm筛后，采用高温外热重铬酸钾氧化-容量法^[14]进行测定。

(3)土壤容重(D_b)测定 采用环刀法测定^[15]。

点击浏览公式

式中，DOC为SOC密度(g/m²)；SOC为SOC重量含量(g/kg)；γ为土壤容重(g/cm³)；H为土层厚度(cm)；δ_2mm为＞2 mm土壤粒径百分含量^[16]。

点击浏览公式

式中，POC为SOC储量(g)；S_i为第i种土地利用的面积(m²)，DOC_i为第i种土地利用类型的SOC重量密度(g/m²)^[5]。

土地利用变化对SOC的影响是一个长期过程，其研究方法包括间接方法和直接方法。间接方法包括相邻样地比较(Paired sites)和空间代替时间方法(Chronosequence)^{[17, 18]}，需要的时间较短、易于操作，但是不够客观，是预测的理想结果，忽略了土地利用变化中的实际演变过程；直接测定的方法为不同时间在同一土地利用方式下重复采样分析，这种方法虽然需要较长的时间和严格的控制条件，但是反映的是样地的实际变化特征，真实客观。本研究鉴于对研究地区的土地利用变化的长期观测、准确的时间记载和系统的研究，采用直接测定比较的方法。

首先在Mapgis软件中，将2006年和2011年通过GPS测得的采样点的经纬度坐标与羊圈沟流域的空间分布坐标相连接，将测得的数据导入羊圈沟流域的点位表中，添加到地形图中，得到研究区域采样点的分布图(图 1)。然后，根据采样点的数值，采用基于土地利用类型修正的普通克里金空间插值法^[19]生成2006年和2011年羊圈沟流域表层的SOC密度图，并提取出不同土地利用类型的平均有机碳密度，并据此分析2006—2011年羊圈沟流域表层SOC密度的变化情况。再对羊圈沟流域2006和2011年的土地利用遥感解译图进行空间叠加生成2006—2011年土地利用变化图，将此图与生成的表层SOC密度的变化图叠加后进行区域统计，即可得出土地利用变化与表层SOC密度的对应关系^[20]。

以羊圈沟流域2006和2011年的土地利用类型图(图 1)为底图得到的2006—2011年各土地利用转移矩阵如表 2所示。其中变化最为剧烈的主要有：草地转化为灌木林地，面积为27.680 hm²，占2006年草地总量的30.35%；草地转化为乔木林地的面积为30.927 hm²，占2006年草地总量的33.91%；坡耕地转为梯田，约为10.999 hm²，占2006年坡耕地总量的43.54%。这是由于这期间流域内加大了水土流失治理措施，调整了土地利用结构，使得流域内乔木林地和灌木林地所占比重增加，遍布整个小流域。由于植树造林，梯田转换为乔木林地和灌木林地的趋势较大，尤其是在流域东边和西边较为明显，大部分荒草地也转换为乔木林地、灌木林地和坡耕地，草地是最大的转出者，乔木林地和灌木林地是最大的转入者。

2006年羊圈沟流域表层土壤平均有机碳密度为10.27 g/m²，其中乔木林地>灌木林地>草地>坝地>梯田>坡耕地；2011年表层SOC密度较2006年略有所增加，平均值为11.01 g/m²，两年的SOC密度空间分布图见图 2，不同地类表层土壤碳密度两年间的变化情况见表 3。从羊圈沟流域SOC密度变化图(图 3)可知，西北地区表层SOC密度以增加为主，其转变用地类型主要为：转变为灌木林地、乔木林地的及保持不变的灌木林地、乔木林地等；而中部及东南地区的地区以减少为主，其转变用地类型主要为转变为及保持不变的坡耕地等。

图 2 2006年、2011年羊圈沟流域表层土壤有机碳密度图 Fig.2 Surface Soil Organic Carbon(SOC)density in Yangjuangou watershed in 2006 and 2011

图选项

图 3 2006年—2011年羊圈沟流域表层土壤有机碳密度空间格局变化 Fig.3 Spatial changes of the surface SOC density in Yangjuangou watershed from 2006 to 2011

图选项

SOC储量比较全面地反映了土地利用变化对SOC的影响。由于不同土地利用类型SOC密度存在着差异，因此SOC储量也存在差异^[5]。由保持利用类型的SOC密度及碳储量变化表(表 4)可知：2006—2011年保持不变的土地利用类型面积61.852 hm²，占整个流域总面积的30.62%。就SOC密度而言，虽然整个流域保持利用类型的SOC密度变化较小，但保持不变的草地、灌木林地、乔木林地、梯田及坝地的变化值均为正；而保持坡耕地的变化值均为负。整个流域保持用地类型表层SOC储量共增加211.64 kg，其中坡耕地减少量最多，为-88.25 kg，而乔木林地增加量最多，为177.53 kg。

转变土地利用类型SOC密度及碳储量变化情况如下：羊圈沟流域从2006到2011年，转变用地类型面积1.25×106 m²，占羊圈沟流域总面积的63.79%，转变土地利用类型的表层SOC密度的变化区间为-0.08—0.14 g/m²。从有机碳密度来看，转入为灌木林地、乔木林地及梯田的类型，有机碳密度多为增加，而转为坡耕地的均减少。

整个流域转变用地类型的SOC储量增加587.25 kg，其中草地转出类型增加的土壤表层有机碳储量最多，为441.64 kg；灌木林地转出类型减少的土壤表层有机碳储量最多，为-21.01 kg。

将小流域2006—2011年保持用地类型表层土壤碳密度的变化值与转变用地类型表层土壤的碳密度变化值进行比较分析发现(图 4),保持用地类型和转变用地类型的表层SOC密度均发生显著变化。草地-灌木林地、草地-乔木林地、坡耕地-草地、坡耕地-灌木林地、坡耕地-乔木林地、坡耕地-坝地、梯田-草地、梯田-灌木林地、梯田-乔木林地、梯田-坝地、坝地-草地、坝地-灌木林地、坝地-乔木林地等转换类型表现为表层土壤碳密度增加值高于保持用地类型碳密度的增加值，说明这些地类的转换有利于表层SOC储量的增加，即有利于表层土壤碳汇的形成；而其他地类之间的转换表现出转变用地类型的表层土壤碳密度增加值低于保持用地类型碳密度的增加值，说明在这些地类转换过程中造成了表层有机碳的释放，不利于土壤碳汇的形成。

图 4 2006年至2011年转变与保持用地类型有机碳密度变化值比较 Fig.4 Comparison of SOC density change between changed and maintained land use types from 2006 to 2011

图选项

研究表明，得益于国家退耕还林政策的实施，羊圈沟小流域的耕地面积明显减少，灌木林地和乔木林地面积显著增加，植被恢复效果明显。通过合理的土地利用方式转变，可以提高土壤有机物质的输入量，减少土壤有机物质的矿化分解，增加SOC汇^[21]。

林地植被覆盖度较高，林下凋落物较厚，拦蓄水土、阻挡表层土壤养分功能较强，且根系分泌物和凋落物输入较大，从而使表层SOC的密度高于其他用地，坡耕地和草地造林后SOC汇的功能强度都将增强^[13]。因而草地-灌木林地、草地-乔木林地、坡耕地-灌木林地、坡耕地-乔木林地、梯田-灌木林地、梯田-乔木林地、坝地-灌木林地、坝地-乔木林地等转变均有利于土壤表层碳汇的形成。造林作为干旱-半干旱区一种重要的生态恢复手段，已在世界范围引起重视和关注^[7]。

坝地可以拦截泥沙、流域地表的有机质及其它大量养分。土壤水分充足，则其透气性差，原有机碳分解速率小，外源有机残体在水分作用下易于腐烂降解成有机物质，从而有利于SOC密度的提高。坝地在流域的上、中和下游均有分布，黄土高原侵蚀剧烈，降雨侵蚀导致的水沙能在很短的时间内被坝地拦截，泥沙淤积迅速，所以SOC被氧化的可能性较小，并搬运到坝地中累积。在流域单元内，淤地坝工程的固碳效应是整个流域固碳的重要方面^[11]。坡耕地土壤水分不足，土壤孔隙度大，加大有机碳分解速率，不利于SOC的积累，因而，坡耕地-坝地的有机碳密度增加。

梯田主要分布在峁坡，其有效地控制了水土流失并承接了来自峁顶的汇水汇沙，因而使处于强烈侵蚀带的峁坡SOC密度升高。而坡耕地以果园为主，果园属于经济林，林种单一，凋谢物少，植物残体归还量少，导致碳输入量低，控制水土流失能力以及承接上方泥沙的能力远远低于梯田，因而坡耕地-梯田会使有机碳密度增加^{[13, 20]}。

草地主要以自然恢复天然草地为主，大面积分布在峁坡，植被覆盖度较高，且有零星的乔木和灌木植被，有效地减少了土壤侵蚀和水土流失。坡耕地转变为草地SOC蓄积量的增加主要来源于高生物量、凋落物留存和高根系生物量^[22]，因而坡耕地-草地的转变有利于有机碳密度的增加。正如Gebhart等^[23]发现美国中部平原的农地变成草地之后300 cm土壤中有机碳增加110.0 gC m^-2 a^-1，该值是最大的增长速率。Lugo等^[24]在亚热带湿润森林区的研究证明了农田被草地代替之后，平均积累速率为33.2 gC m^-2 a^-1等。然而草地拦截水土能力弱于林地，天然草地根系分泌物和凋落物输入高于农田和坡耕地(果园)但低于乔木林地和灌木林地^{[23, 24]}。

由于实施退耕还林工程，从2006年以来，流域内坡耕地基本消失，仅为0.23 hm²，到2011年，原有坡耕地逐渐转换为梯田及灌木林地，但是仍有少量处于坡面的荒草地被开发为耕地，于是出现了少量坡耕地。转变为坡耕地(果园和旱地)从而引起SOC损失的机制为：首先由于耕作的物理效应，即耕作破坏了土壤的团聚体结构，使土壤透气性、温度、水分条件得到改善，微生物呼吸作用增强，加速了SOC的分解，导致SOC的大量释放；其次，开垦使土壤的温度和湿度条件在耕作中得到了改善，在一定程度上促进了土壤的呼吸作用，加速了土壤有机质的分解；另外，收获作物时大量地上部分被移走，从而大大减少了植被对土壤碳的归还量，耕地的烧荒措施也使固定在植被中的碳全部释放到大气中^[22]。因而草地-坡耕地造成土壤碳密度的减少，应该降低土壤耕作强度，减少土壤呼吸强度，从而达到减少SOC消耗的目的^[5]。

本研究中2006—2011年小流域尺度土壤碳密度虽整体虽呈增加趋势，其变化幅度却很小。这是由于SOC的累积和恢复是一个长期过程，黄土丘陵区植被恢复应着眼于长期效应，频繁的土地利用变化可能会降低土壤碳截流效果，而长期保持更有利于本区域土壤固存大气中的CO₂^[25]。随着退耕还林的开展和坝地的增加，从生态系统流失的碳可以通过固碳重新被土壤和植被吸收，从而恢复退化的土壤，增加生态系统碳储量^[26]。由图 5可知，乔木林地与灌木林地虽均为碳汇^[16]，然而互相转换后均造成SOC密度的降低，即不利于土壤碳的保持，因而就林地来说，更适合保持不变，转换树种后，由于对土壤进行了扰动，反而造成了碳损失。增加生态系统的碳汇功能应从增加输入量、减少输出量和增加稳定性去实现^[7]。

(1) 2006年到2011年，羊圈沟流域土地利用类型发生较大变化，主要集中在乔木林地和灌木林地面积的增加以及荒草地及坡耕地面积的减少。黄土丘陵区植被恢复的土壤固碳效应表现出一定的时空变化规律，对评估本区域植被恢复的SOC截流量，以及预测未来SOC截流潜力有重要影响。

(2) 土地利用方式的变化会导致SOC密度及储量的变化，其中，转变为林地、梯田、坝地等类型的土地利用方式，有利于表层SOC储量的增加，应该保留；而转变为坡耕地的类型造成了表层土壤的碳排放，应该引起足够的重视。然而在研究区域内，居民点附近仍出现少量将荒草地开发为坡耕地的现象，因此应加强植被恢复的力度和管理强度，兼顾淤地坝、梯田的维护和建设。

(3) 土壤固碳应着眼于长期效应，频繁的土地利用类型转化可能会降低土壤碳截流效果，黄土丘陵区植被重建的长期利用和保持更有利于SOC的积累。在野外条件下开展关于土地利用变化后SOC库的研究，可以在最接近自然状况的条件下揭示SOC库与土地利用变化相互作用的规律，但会受到土壤异质性及人为频繁扰动的影响，虽然可以通过重复样地和多点采样进行弥补，仍会造成一些变异性和不确定性结果，或会掩盖某些规律性的变化，因此今后研究应辅助以可控性更强的室内模拟试验。