藏东南地区不同土地利用方式下土壤有机碳组分及 周转变化特征

徐梦，李晓亮，蔡晓布，李晓林，张旭博，张俊伶



藏东南地区不同土地利用方式下土壤有机碳组分及 周转变化特征

徐梦1,2，李晓亮3，蔡晓布4，李晓林2，张旭博1，张俊伶2

（1中国科学院地理科学与资源研究所生态网络观测与模拟重点实验室，北京 100101；2中国农业大学资源与环境学院/植物-土壤相互作用教育部重点实验室，北京 100193；3中国热带农业科学院热带作物品种资源研究所/农业部华南作物基因资源与种质创制重点开放实验室，海南儋州 571700；4西藏大学农牧学院，西藏林芝 860000）

【目的】自然植被转变为农业用地显著影响土壤有机碳储量。青藏高原东南部地区森林或草地转换为农田的面积逐年增加，但其对土壤有机碳组分及周转特征的影响尚不明确。因此，阐明藏东南地区不同土地利用方式对土壤有机碳储量的影响程度和作用机制，可为该地区农业土地资源合理利用提供科学依据。【方法】采集藏东南地区长期耕作的农田（50年以上）及毗邻的自然森林和草地土壤，采用物理-化学联合分组技术以及稳定性碳同位素测定，分析3种土地利用方式下土壤有机碳组分的数量、碳含量的差异，探究不同有机碳组分周转差异及其对农田耕作的响应规律。【结果】农田0—20 cm表层土壤有机碳储量为（39.4±2.0） Mg C·hm-2，比自然森林的（81.5±8.5）Mg C·hm-2和草地的（71.4±7.3）Mg C·hm-2分别降低了约52%和45%。农田耕作导致粗颗粒有机质（cPOM）数量相对于自然植被降低了63.4%—70.8%，微团聚体（μagg）和黏粉粒（dSilt+Clay）的数量分别增加了10.0%—25.9%和65.7%—86.2%。农田土壤的有机碳含量与森林和草地土壤相比降低了51.7%—58.1%，其中不稳定性、物理稳定性和生物化学稳定性有机碳库分别降低79.8%—86.3%、72.4%—73.1%、32.4%—39.8%，且与总有机碳的变化显著正相关，但化学稳定性有机碳库没有显著变化。土地利用方式不同导致不同有机碳组分的C/N值和13C值差异明显。农田土壤cPOM组分的C/N值（10.0±0.5）显著低于森林（13.5±0.4），而13C值（-21.6±0.5）‰则显著高于森林土壤（-23.6±0.4）‰。微团聚体保护的颗粒有机质（iPOM）和难酸解组分（NH-dSilt+Caly和NH-μSilt+Clay）具有较低的13C值（-25.3‰—-27.2‰），并且其C/N在农田土壤为8.4—9.4，显著低于森林土壤（13.5—15.9）。【结论】藏东南地区长期耕作的农田土壤有机碳储量相比于自然植被降低了约50%。农业耕作显著加速了不稳定颗粒有机质的周转，减少了稳定性有机碳组分如微团聚体保护的有机碳组分的形成，是导致土壤有机碳库明显下降的关键原因。因此，为有效降低农业耕作对土壤有机碳储量的负面影响，免耕和保护性耕作或可成为藏东南农耕区固碳增汇、维持该地区土地资源可持续利用的技术选择之一。

土壤有机碳；土地利用方式转换；农田；有机碳分组；稳定性碳同位素；青藏高原

0 引言

【研究意义】由于独特的高海拔地势和气候环境，青藏高原的土壤有机碳储量十分巨大[1-2]，并且对气候变化以及人类活动的干扰十分敏感[3]，在中国乃至全球碳源汇管理中起着举足轻重的作用。因此，明确青藏高原地区土壤有机碳储量对外界环境变化的响应规律尤为重要。土地利用方式的变化是影响土壤有机碳储量的重要因素之一[4]，由于森林或草地开垦为农田耕作而导致的土壤有机碳储量下降在全球各个区域均有报道[5-10]。青藏高原东南地区土地开发利用历史悠久、格局多样，森林资源十分丰富，同时也是西藏粮食重要产区[11]。近年来由于该区域经济发展和人口增长压力的需求，农业用地面积逐年增加，不仅对维持森林生态系统带来了极大的挑战，也对土壤有机碳储量产生潜在影响。因此，研究该区域不同土地利用类型（自然植被以及农田）对土壤有机碳储量的影响规律，不仅有助于深入认识青藏高原土壤有机碳稳定性，也能够为该区域农、林业土地资源的合理利用和可持续管理提供科学依据。【前人研究进展】土地利用方式的改变能够通过调控输入到土壤中有机碳的数量和质量[12]、影响土壤理化性质及过程[13-14]、改变土壤微生物群落结构及活性[15-16]，从而影响土壤有机碳的分解速率和稳定性[17]，最终导致土壤有机碳储量变化。近年来，研究人员对草地放牧[18]及草场管理方式[19]如何影响青藏高原土壤碳库储量及碳循环过程有了一些认识，但是对农田开垦对土壤有机碳储量的影响关注较少。并且以往研究中多以原状土有机碳总含量作为主要分析指标，从而无法深入辨析不同土地利用方式下土壤有机碳储量变化的特征及其机制。通过物理或化学分组法，土壤有机碳库可以被分解为许多在功能性或者生物组成上不同的组分[20]，而这些有机碳组分能够更加敏感地反映出土地利用方式或者管理措施的变化对土壤有机碳的影响并揭示其生物地球化学机制[21-23]。此外，利用稳定性碳同位素分析方法能够揭示土地利用方式变化对土壤有机碳动态周转的影响[24-25]，这是由于有机碳在周转过程中会发生同位素分馏，12C优先被分解成为CO2释放，从而使残留有机碳富集13C[26]。QIAO等[27]在青藏高原海北高寒草甸的研究发现，种植燕麦10年后其表层土壤13C值（-23.1‰）明显高于周边自然嵩草牧场（-25.6‰），农田种植加速土壤有机碳周转，导致土壤有机碳储量下降了10%。然而，土地利用方式变化如何影响各级土壤有机碳组分的数量和动态周转，尚缺乏深入研究。【本研究切入点】藏东南地区自然植被（森林和草地）转变为农田种植后，输入土壤的有机碳数量和质量变化、农田管理措施等因素能够影响土壤有机碳循环过程，加速土壤有机碳周转，进而降低土壤有机碳储量。然而，不同土地利用方式下土壤有机碳组成的差异性，尤其是明确哪些土壤有机碳组分的数量和周转对农田开垦的响应最为敏感，进而影响到土壤有机碳储量变化，亟需进行深入分析。【拟解决的关键问题】本研究针对藏东南地区耕作超过50年以上的农田以及自然森林和草地3种土地利用方式，利用物理-化学联合分组法对具有不同周转和稳定性机制的有机碳库组分进行分组，明确不同土地利用方式下土壤有机碳组成的差异性，并测定各个组分的碳含量以及稳定性碳同位素组成，以阐明土地利用方式的变化对土壤有机碳储量的影响程度及其作用机制，为藏东南地区农业土地资源合理利用提供科学依据。

1 材料与方法

1.1 采样区域概况

本试验在青藏高原东南部林芝地区（N 26°52′—30°40′，E 92°09′—98°47′）进行。该区平均海拔约为3 000 m，年平均温度为7—16℃，＞10℃年积温为2 272℃，无霜期约为180 d，年平均降水量约为600—800 mm。林芝地区约有46%的面积为森林（2.64×106hm2），占西藏全部森林面积的80%。草地面积为2.91×105hm2，其中适合做放牧草场的面积约为5.25×104hm2。林芝地区是西藏粮食重要产区之一，该区的传统农业耕作为休闲和单作，化肥投入少，机械化程度低。适宜农业耕作的土地面积为3.00×104hm2，目前已被用作耕地的土地面积为1.87×104hm2。该区种植的粮食作物主要为小麦和青稞，近年来随着大棚技术的推广，蔬菜和水果的产量也逐渐增加。根据林芝地区统计局发布的《林芝地区经济和社会发展统计公报》，在2011—2013年间耕地面积以每年1.4%的速率增加，其中小麦和蔬菜种植面积增加速率分别为4%和1.6%。

1.2 土壤样品采集

本试验中土壤样品的采集于2012年7月进行。在林芝地区选取3个代表当地农田种植管理模式的典型农田（种植小麦）作为研究对象，农田耕作时间均为50年以上，施用少量有机肥，不使用除虫剂，每年作物种植前进行一次翻耕处理。在每个农田中，随机选取5个样方（5 m×5 m），在每个样方中随机采集3个表层（0—20 cm）土壤样品（20 cm×20 cm）并混合为一个样品重复。在每个农田附近（距离500 m以内）的自然植被下以同样的采样方法进行土壤样品的采集，一共采集45个土壤样品。所选取的森林是以川滇高山栎（）和林芝云杉（var.）为主要树种构成的温带针阔混交林次生林，林下小乔木和灌木包括花楸（）、小檗（）、杜鹃（L.）等。所选取的草地属于自然草地，轻度放牧，主要物种包括紫菀（Aster）、苔草（Carex）、禾本科（Poaceae）、蛇莓（）、米口袋（）等。将采集的土壤样品过2 mm筛，挑出可见的石块、土壤动物和植物残体。用于测定土壤有效氮（NH4+-N和NO3--N）含量的土壤样品则保存在4℃冰箱里，测定土壤其他指标的样品风干后室温下保存。

1.3 试验方法

1.3.1 土壤理化性质测定 土壤pH用1 mol·L-1KCl浸提（土水比为1﹕2.5）并用电位法测定。土壤总碳、总氮含量使用元素分析仪测定（EA1108, Carlo Erba, Turin, Italy）。土壤有效磷（Olsen-P）使用0.5 mol·L-1NaHCO3浸提-钼锑抗比色法测定。土壤有效氮含量的测定使用0.01 mol·L-1CaCl2溶液浸提（土水比1﹕10）并用流动分析仪（TRAACS 2000, Bran and Luebbe, Norderstedt, Germany）测定其中的NH4+-N和NO3--N含量，有效氮含量为两者之和。土壤容重使用环刀法测定。土壤团聚体分组使用湿筛法：首先将250、53和20 μm的筛子按次序放置于特制铁架上放入铁桶内，并将100 g风干土壤样品（过2 mm筛）均匀放置于最上层的250 μm筛子上。向铁桶内加入2 L去离子水并浸泡5 min后开始振荡，以振幅4 cm、频率35次/min持续振荡10 min。振荡结束后，将各级筛子上的土壤颗粒分别洗入铝盒中，置于105℃烘干，待铝盒内无明显积水后，再烘干24 h并称重。粒径为250—2000 μm的土壤颗粒为大团聚体（macroaggregate），粒径为20—53 μm和53—250 μm的土壤颗粒为微团聚体（microaggregate）。

1.3.2 有机碳物理-化学联合分组 物理-化学综合法是根据SIX等[21]提出的有机碳稳定性模型，将物理和化学方法相结合的有机碳分组方法，具体步骤根据STEWART等[28]进行。第一步，将50 g风干土样均匀放置于250 μm+53 μm套筛上，同时放入50个直径4 mm的玻璃珠，用去离子水冲洗并振荡，残留在上层250 μm筛上的组分为不稳定粗颗粒有机质（cPOM）；在下层53 μm筛上的组分继续用去离子水冲洗，并在2 min内振荡50次，残留在筛上的组分为微团聚体（μagg）；套筛下大烧杯中收集的滤下组分为易分散的粉粒和黏粒（dSilt+Clay）。第二步是将第一步中收集的μagg组分进一步分组，首先用1.7 g·cm-3NaI溶液进行密度梯度离心，离心后的上层组分过0.45 μm滤膜，滤膜上的组分为不稳定细颗粒有机质（LF）；离心后下层剩余组分转移到53 μm筛上，用去离子水冲洗振荡，残留在筛上的组分为微团聚体保护的颗粒有机质（iPOM）；大烧杯中收集的滤下组分为微团聚体中的粉粒和黏粒（μSilt+Clay）。第三步是用酸解法对前两步中收集的粉粒和黏粒组分进行化学分组，分别将dSilt+Clay组分和μSilt+Clay组分置于玻璃试管中，用6 mol·L-1HCl在95℃下消煮16 h，之后用去离子水清洗残留组分，即难酸解组分（NH-dSilt+Clay和NH-μSilt+ Clay）；易酸解组分（H-dSilt+Clay和H-μSilt+Clay）为酸解前后组分之差。根据不同组分对有机碳保护机制的差异，将得到组分划分为不稳定碳库（unprotected C pool）、物理稳定性碳库（physically protected C pool）、化学稳定性碳库（chemically protected C pool）和生物化学稳定性碳库（biochemically protected C pool）[21]。不稳定碳库包括cPOM和LF组分；物理稳定性碳库包括μagg组分以及其中的iPOM组分；易酸解和难酸解组分分别为化学稳定性和生物化学稳定性碳库。

收集所有组分，于65℃烘干并称重，各组分碳氮含量和稳定性碳同位素组成值（13C）使用DeltaPlusXP质谱仪-元素分析仪（Thermo Scientific，Bremen，Germany）及元素分析仪（FlashEA 1112；CE Instruments，Wigan，UK）测定。元素分析仪的燃烧温度设定为1 020℃。测定样品的同位素组成（）用样品的同位素比值（R样品）与标准物质的同位素比值（R标准，对碳同位素使用V-PDB标准）的千分差表示。具体公式为：

其中R为碳元素重轻同位素丰度比（13C/12C）。碳同位素组成和碳含量测定的精确度分别为0.2‰和0.1%。

1.3.3 数据分析 各土壤样品的容重根据下式[4]进行估算：

BD=

其中，BD为土壤容重（Bulk density，g·cm-3），SOM为土壤有机质（%），利用土壤有机碳含量（SOC %）计算可得（SOM % = SOC %×1.724）。

土壤有机碳储量（Mg C·hm-2）根据下式进行计算：

SOC stock=BD×SOC(%)×depth

其中BD为估算出的土壤容重，depth为需要计算有机碳储量的各层剖面土壤深度（cm）。

使用单因素方差分析（ANOVA）比较土地利用方式对土壤理化性质及土壤有机碳库各组分的差异，并用Duncan’s multiple range test及Paired t-test检验各处理间差异的显著性（在＜0.05水平下）。使用Pearson相关分析土壤理化性质和土壤有机碳库特征（各组分碳含量及碳氮比）的相关性。所有统计分析都使用IBM SPSS Statistics（version 22.0.0.0）进行。

2 结果

2.1 土地利用方式对土壤理化性质的影响

不同土地利用方式下土壤化学性质有显著差异（表1）。与森林和草地土壤相比，农田土壤的pH分别升高了2.1和1.5个单位，土壤有机碳、总氮以及有效氮含量分别降低了51.7%—58.1%、38.8%—39.8%和63.0%—70.6%，土壤有效磷含量则增加了32.1%—52.1%。农田种植显著影响土壤物理性质，与森林和草地土壤相比，农田土壤容重增加了12.7%—13.6%，大团聚体数量降低了29.8%—36.3%，而微团聚体数量则增加了34.6%—49.3%。农田土壤粉粒和黏粒含量显著高于森林和草地土壤，而砂粒含量则显著降低。据估算，农田土壤表层20 cm的有机碳储量为（39.4±2.0）Mg C·hm-2，较森林的（81.5±8.5）Mg C·hm-2和草地的（71.4±7.3）Mg C·hm-2分别降低了51.6%和44.8%。

表1 3种土地利用方式下土壤理化性质及有机碳储量

数据表示为每种土地利用方式下的平均值±标准误差。数据后不同的小写字母表示各指标在不同土地利用方式下差异显著（＜0.05）。下同

Data are presented by mean value ± SE. Different lowercase letters indicate significant differences among three land use types by Duncan’s multiple range tests (＜0.05). The same as below

2.2 农田耕作对土壤有机碳组成的影响

与自然森林和草地相比，农田土壤中粗颗粒有机质cPOM的数量降低了63.4%—70.8%，黏粉粒dSilt+Clay的数量增加了65.7%—86.2%。农田土壤微团聚体μagg的数量与草地土壤相近，并显著高于森林土壤，但是微团聚体中的不稳定细颗粒有机质LF和稳定颗粒有机质iPOM的数量均显著低于森林和草地土壤，而微团聚体中的粉粒和黏粒μSilt+Clay的数量则分别是森林和草地土壤的4.9和3.6倍。各土壤有机碳组分的碳含量都表现出农田小于森林/草地的变化特征（表2）。

根据不同保护机制划分的4种土壤有机碳库中，不稳定性碳库、物理稳定性和化学稳定性碳库是主要组成部分，其碳含量总和超过土壤总有机碳库的90%（图1）。在自然植被土壤中，不稳定性碳库和物理稳定性碳库是土壤有机碳库的主要组成部分，这两个碳库的总和分别占到森林和草地土壤有机碳库的74.3%（22.9 g·kg-1）和65.1%（17.3 g·kg-1）。自然植被转变为农田后，有机碳库组成发生显著变化，不稳定性碳库和物理稳定性碳库占总有机碳的比例分别降低了49.6%—62.1%和27.9%—41.0%，而化学稳定性碳库则增加了0.9—1.6倍，成为农田土壤有机碳库的主要组分（6.8 g·kg-1，占有机碳总量60.4%）。在所有土地利用方式下，生物化学稳定性碳库只占土壤总有机碳库很小一部分（小于6%），其含碳量在农田土壤中显著低于森林或草地土壤。自然植被转变为农田后土壤总碳库的变化与不稳定性碳库（=0.949，＜0.001）、物理稳定性碳库（=0.924，＜0.001）及生物化学稳定性碳库（=0.792，＜0.001）呈现出极显著的正相关关系（图2）。土壤不稳定性、物理稳定性和生物化学稳定性碳库的碳含量与土壤总氮、C/N值、大团聚体和砂粒含量显著正相关，与土壤pH、容重、微团聚体及粉粒和黏粒含量显著负相关，而化学稳定性碳库的碳含量只与土壤总氮含量显著正相关（表3）。

2.3 农田耕作对土壤有机碳组分周转的影响

不同有机碳组分的碳氮比（C/N）存在显著的差异（图3-A）。在森林土壤中，NH-μSilt+Clay的C/N值最高（15.9±1.0），其次为cPOM（13.5±0.4）和NH-dSilt+Clay（13.5±0.7）、μagg（12.2±0.1）和iPOM（12.4±0.3），dSilt+Clay（11.4±0.2）和μSilt+Clay（11.7±0.2）的C/N值最低。草地土壤各个有机碳组分C/N值的变化规律与森林类似。而在农田土壤中，dSilt+Clay（11.0±0.2）和μSilt+Clay（10.8±0.3）的C/N值最高，其次为cPOM（10.0±0.5）和μagg（10.1±0.3），iPOM（8.4±0.2）、NH-dSilt+Clay（8.4±0.3）和NH- μSilt+Clay（9.4±0.2）的C/N值相对较低。除了黏粉粒组分（dSilt+Clay和μSilt+Clay），其他各组分的C/N值都表现出在森林土壤中显著高于草地和农田土壤。

表2 3种土地利用方式下各土壤有机碳组分的数量与碳含量

柱状图中每个柱子左边标注的数字表示各碳库的碳含量（g·kg-1），数据后的小写字母表示各碳库的碳含量在不同土地利用方式下差异显著（P＜0.05）。图例中不同的小写字母表示各碳库占总有机碳百分比在不同土地利用方式下差异显著。* P＜0.05; ** P＜0.01; *** P＜0.001

图2 各有机碳库组分与总有机碳库的相关性

表3 土壤总有机碳库及各有机碳库组分与土壤理化指标的相关性

加粗的数字表示相关性显著。*＜0.05; **＜0.01; ***＜0.001

Values presented in bold indicate significant correlations. *＜0.05; **＜0.01; ***＜0.001

与其他有机碳组分相比，cPOM和μagg相对富集13C，并且其在农田土壤的13C值（（-21.6±0.5）‰和（-22.4±0.3）‰）显著高于森林（（-23.6±0.4）‰和（-24.7±0.2）‰）和草地土壤（（-23.9±0.4）‰和（-24.0±0.5）‰）（图3-B）。黏粉粒组分（dSilt+Clay和μSilt+Clay）的13C值也相对较高，并且在农田（（-22.6±0.6）‰和（-23.9±0.3）‰）显著高于森林土壤（（-24.8±0.1）‰和（-25.2±0.2）‰）。微团聚保护的颗粒有机质iPOM的13C值在农田（-25.1±0.4）‰和森林（-25.8±0.5）‰显著低于草地土壤（-23.7±0.9）‰。难酸解组分（NH-dSilt+Clay和NH-μSilt+Clay）的13C值在3种土地利用方式下均明显低于其他组分，其变化范围为-25.3‰—-27.2‰。

图柱上的不同小写字母表示每种土地利用方式下不同有机碳组分的C/N值或δ13C值差异显著（P＜0.05）。图例中不同的小写字母表示各有机碳组分的C/N值或δ13C值在不同土地利用方式下有显著差异。* P＜0.05；** P＜0.01；*** P＜0.001

3 讨论

基于全球数据的Meta分析结果发现，森林转变为农田后其碳库平均损失42%，而放牧草原转变为农田后其碳库平均损失59%[4]。本研究中，相比于自然植被，农田种植导致土壤有机碳碳含量显著降低51.7%—58.1%（表1），农田土壤表层20 cm的有机碳储量（39.4 Mg C·hm-2）相比于森林（81.5 Mg C·hm-2）和草地（71.4 Mg C·hm-2）分别降低了51.6%和44.8%，与全球平均水平相近，高于中国黄土高原（34%）[29]，而低于热带雨林（78%）[6]。研究表明，不同土地利用方式对土壤有机碳储量的影响程度在全球各区域存在很大的变异，主要是由于各区域气候因子（例如年均温和降水量）的差异能够影响土壤有机碳周转速率[30-31]。另外，农田开垦的时间长度也能够影响土壤有机碳降低的程度[29,32]。QIAO等[27]在海北高寒草原的研究发现，种植燕麦10年后表层30 cm土壤有机碳储量相对于自然嵩草草原减少15 Mg C·hm-2（10.1%），显著低于本研究中农田耕作50年后土壤有机碳储量下降的程度，这说明青藏高原土壤有机碳储量的减少很可能随着农田开垦的时间延长而加剧。

农田开垦导致土壤有机碳储量的减少可进一步从不同有机碳组分的变化反映出来（图2）。不稳定性碳库和物理稳定性碳库是自然森林和草地土壤有机碳的主要组成组分（图1）。农田种植后土壤不稳定性碳库cPOM和LF组分的数量和碳含量均下降（表2），不稳定性碳库总量减少最高可达到86%。cPOM和LF组分主要由未被土壤矿物固持的半分解植物残体组成[21,33]。农业生态系统中，通过净初级生产力向土壤中输入的碳源（主要来自于作物根系及残茬）的数量通常少于自然生态系统（凋落物以及地下碳输入）[12]。并且，输入碳源中木质素、单宁酸以及其他植物次级化合物的含量往往相对低于自然生态系统的凋落物[22,34]，因而具有较快的分解速率[17]。本研究中，农田土壤cPOM的C/N值（10.0±0.5）显著低于森林（13.5±0.4），并且13C值（-21.6±0.5）‰显著高于森林（-23.6±0.4）‰（图3）。研究表明，土壤C/N值会随着微生物对有机碳分解程度的增加而降低[35]，而13C值则会相应增加[36]。因此，这一结果表明农田耕作显著促进了这部分有机碳组分的分解和周转，导致土壤有机碳储量减少。此外，cPOM和LF组分中易分解碳源和分解程度较低的植物残体是土壤团聚体形成的重要介质，对大团聚体的形成起到关键作用[21]。农田耕作显著加速了这两个有机碳组分的分解，导致土壤团聚体稳定性的下降，从而可能进一步加剧土壤有机碳储量的减少。

与不稳定性碳库的变化相似，土壤物理稳定性碳库在农田土壤中显著低于森林和草地土壤（图1），其总量减少约为73%。SIX等[37]提出，农业耕作能够破坏团聚体结构，加速有机碳矿化速率，使得富含有机碳的大团聚体数量减少，而有机碳含量较低的微团聚体数量则会增加。本研究也发现，自然植被转变为农田后，土壤大团聚体（> 250 μm）数量显著降低，微团聚体（以及μagg组分）数量明显增加但其碳含量却显著降低（表1和表2），并且这些团聚体结构的变化与土壤有机碳储量的减少显著相关（表3）。团聚体对土壤有机碳的保护作用是土壤碳库稳定的重要机制[38]，农业耕作导致的团聚体结构破碎而引起的的碳库损失可达到50%以上[10,29,39]，而免耕则有利于团聚体稳定性并增加土壤有机碳固持[40]。本研究中，微团聚体保护的颗粒有机质（iPOM组分）具有较低的13C值（图3-B），并且其碳含量在农田土壤中显著低于自然森林和草地土壤（表2）。前人研究也表明，微团聚体对土壤有机碳有很强的保护作用[28]，并且是反映土壤有机碳库变化响应土地利用变化、管理措施等重要指标[23]。此外，SIX等[37]提出的碳周转模型中认为，农业耕作能够加速不稳定性有机碳组分的周转，减少稳定性有机碳组分的形成，是导致土壤有机碳库的损失的主要原因[24]。本研究的结果正是对这一假说的科学验证：自然植被转变为农田种植后，粗颗粒有机质cPOM组分（不稳定性碳）周转加速，团聚体结构破碎导致稳定性有机碳组分（iPOM组分）减少（表2，图3），从而导致土壤有机碳储量显著下降。

化学稳定性碳库主要是黏粒和粉粒固持的有机碳组分（矿物结合态有机碳），其周转时间为几十年至百年[41]，对土壤固碳有十分重要的意义。矿物结合态有机碳是农田土壤有机碳的主要组分[35,42]。本研究中，化学稳定性碳库（H-dSilt+Clay和H-μSilt+Clay）占土壤总碳库的比例在农田土壤中显著高于森林和草地土壤，但是其总量在3种土地利用方式下没有显著的差异（图1）。与本研究类似，BAYER等[43]也发现农田耕作显著降低总有机碳和颗粒有机碳含量，但对矿物结合态有机碳没有显著影响。然而，也有一些研究表明农田耕作后土壤有机碳含量的变化归因于矿物结合态有机碳的减少[35,44]，但是两者的相关性随土壤类型而变化[42]。蔡岸东等[42]发现在农田、草地和林地3种土地利用方式下矿物结合态有机碳含量与粉黏粒含量显著正相关。而本研究发现，化学稳定性碳库的数量与土壤粉黏粒含量在森林（=0.640，=0.01）和草地（=0.604，＜0.05）具有正相关关系，但是在农田中具有负相关关系且不显著（=-0.302，=0.274）。SIX等[21]也观察到粉黏粒含量与其所固持碳含量的相关性随土地利用方式而变化，这很可能是由于农业耕作促使有机碳降解加速，而减少了被固持的数量。另外，一部分有机碳可能并未被矿物固持，只是由于物理作用被截留在粉粒尺径组分中[45]，且无法被现有有机碳分组方法分离出来，从而可能影响到对化学稳定性碳库的准确定量。

生物化学稳定性碳库主要是由化学结构复杂、难以降解的含碳化合物构成[21]，与CENTURY模型中的惰性碳库（passive C pool）十分类似，其含有的碳十分古老因而相对稳定[46]。与其他有机碳组分相比，生物化学稳定性碳库组分NH-dSilt+Clay和NH-μSilt+ Clay的13C值较低（图3-B），说明这两个组分的周转相对缓慢[26]，其稳定性强于其他有机碳组分。不过，13C值的变化也可能与碳源底物的分解性有关，因为容易分解的碳源，例如糖和纤维素等，要比难分解碳源例如木质素和脂质等富集13C[47]。本研究采样区域森林主要由川滇高山栎（）和林芝云杉（var.）构成，其凋落物中木质素和蜡质含量较高[48]，并且土壤中NH-dSilt+ Clay和NH-μSilt+Clay组分的C/N值较高（图3-A），也说明这部分有机碳组分的稳定性机制很可能是由于其化学组成的复杂性造成的。然而，农田土壤NH-dSilt+Clay和NH-μSilt+Clay的13C值与森林土壤没有显著差异，但是C/N值却显著低于森林土壤（图3）。这一结果说明这部分有机碳组分在农田土壤中的稳定性机制并不是其自身的难降解性，而很可能是由于黏粉粒的增加导致矿物固持作用的增强，或者是微团聚体的增加减少了与微生物分泌酶的接触，从而显著降低了周转速率[49]。尽管生物化学稳定性碳库与土壤有机碳储量的变化显著相关（图2），但由于其占总有机碳库的比例很小（小于6%）而且十分稳定，因此与其他3个有机碳库相比，不适宜作为表征土地利用方式转变对土壤有机碳储量影响的研究指标。

4 结论

西藏东南部自然森林、草地和农田3种不同土地利用方式显著影响了土壤理化性质。相比于森林和草地，农田土壤总碳、总氮、大团聚体和砂粒含量显著降低，而pH、容重、微团聚体、粉粒和黏粒含量显著增加。超过50年的农田耕作导致土壤有机碳储量降低约50%，其中土壤不稳定性碳库（cPOM组分）和物理稳定性碳库（μagg和iPOM组分）的响应最为剧烈，其有机碳总量降低可达到73%以上。各个土壤有机碳组分的C/N值和13C值有明显差异，并且受到土地利用方式的显著影响。综合分析有机碳分组以及稳定性碳同位素数据表明，农业耕作显著加速了不稳定颗粒有机质的周转，减少了稳定性有机碳组分如微团聚体保护的有机碳组分的形成，是导致农田土壤有机碳库显著低于自然植被的关键原因。西藏高海拔地区自然植被开垦为农田种植是区域经济发展和人口增长压力的必然选择，应加大力度逐步推广免耕和保护性耕作，有效降低农业耕作对土壤有机碳储量的负面影响，以维持该区域土地资源的可持续利用。

[1] 李克让, 王绍强, 曹明奎. 中国植被和土壤碳贮量. 中国科学, 2003, 33(1): 72-80.

LI K R, WANG S Q, CAO M K. Carbon storage of vegetation and soil in China., 2003, 33(1):72-80. (in Chinese)

[2] YANG Y, FANG J, TANG Y, JI C, ZHENG C, HE J, ZHU B. Storage, patterns and controls of soil organic carbon in the Tibetan grasslands., 2008, 14(7):1592-1599.

[3] 李娜, 王根绪, 高永恒, 籍长志. 青藏高原生态系统土壤有机碳研究进展. 土壤, 2009, 41(4): 512-519.

LI N, WANG G X, GAO Y H, JI C Z. On soil organic carbon of alpine ecosystem in Qinghai-Tibet Plateau., 2009, 41(4): 512-519. (in Chinese)

[4] GUO LB, GIFFORD R M. Soil carbon stocks and land use change: a meta analysis., 2002, 8(4): 345-360.

[5] MARTENS D A, REEDY T E, LEWIS D T. Soil organic carbon content and composition of 130-year crop, pasture and forest land-use managements., 2003, 10(1):65-78.

[6] YONEYAMA T, DACANAY E V, CASTELO O, KASAJIMA I, PARK T H. Estimation of soil organic carbon turnover using natural13C abundance in Asian tropics: a case study in the Philippines., 2004, 50(4): 599-602.

[7] BEHESHTI A, RAIESI F, GOLCHIN A. Soil properties, C fractions and their dynamics in land use conversion from native forests to croplands in northern Iran., 2012, 148(2): 121-133.

[8] DING F, HU Y L, LI L J, LI A, SHI S, LIAN P Y, ZENG D H. Changes in soil organic carbon and total nitrogen stocks after conversion of meadow to cropland in Northeast China., 2013, 373(1/2): 659-672.

[9] SINGH A K, RAI A, SINGH N. Effect of long term land use systems on fractions of glomalin and soil organic carbon in the Indo-Gangetic plain., 2016, 277: 41-50.

[10] ROLANDO J L, DUBEUX JR J C, PEREZ W, RAMIREZ D A, TURIN C, RUIZ-MORENO M, COMERFORD N B, MARES V, GARCIA S, QUIROZ R. Soil organic carbon stocks and fractionation under different land uses in the Peruvian high-Andean Puna., 2017, 307: 65-72.

[11] 关卫星, 焦国成, 刘启勇, 高小丽. 西藏一江两河地区农业耕作制度的现状与改革对策. 西藏农业科技, 2012, 34(4): 44-48.

GUAN W, JIAO G, LIU Q, GAO X. Present situation and reform countermeasure on agricultural tillage system in the “YLN” region of Tibet., 2012, 34(4): 44-48. (in Chinese)

[12] PAUL E A. The nature and dynamics of soil organic matter: Plant inputs, microbial transformations, and organic matter stabilization., 2016, 98: 109-126.

[13] CELIK I. Land-use effects on organic matter and physical properties of soil in a southern Mediterranean highland of Turkey., 2005, 83(2):270-277.

[14] WILSON B R, KOEN T B, BARNES P, GHOSH S, KING D. Soil carbon and related soil properties along a soil type and land-use intensity gradient, New South Wales, Australia., 2011, 27(4): 437-447.

[15] BISSETT A, RICHARDSON A E, BAKER G, THRALL P H. Long-term land use effects on soil microbial community structure and function., 2011, 51: 66-78.

[16] RODRIGUES J L M, PELLIZAZI V H, MUELLER R, BAEK K, JESUS E D, PAULA F S, MIRZA B, HAMAOUI G S, TSAI S M, FEIGL B, TIEDJE J M, BOHANNAN B J M, NUSSLEIN K. Conversion of the Amazon rainforest to agriculture results in biotic homogenization of soil bacterial communities., 2013, 110(3): 988-993.

[17] GRÜZWEIG J M, SPARROW S D, CHAPIN III F S. Impact of forest conversion to agriculture on carbon and nitrogen mineralization in subarctic Alaska., 2003, 64(2): 271-296.

[18] HAFNER S, UNTEREGELSBACHER S, SEEBER E, LENA B, XU X L, LI X G, GUGGENBERGER G, MIEHE G, KUZYAKOV Y. Effect of grazing on carbon stocks and assimilate partitioning in a Tibetan montane pasture revealed by13CO2pulse labeling., 2012, 18(2): 528-538.

[19] CHANG X F, ZHU X X, WANG S P, CUI S J, LUO C Y, ZHANG Z H, WILKES A. Impacts of management practices on soil organic carbon in degraded alpine meadows on the Tibetan Plateau., 2014, 11(13): 3495-3503.

[20] VON LÜTZOW M, KÖGEL-KNABNER I, EKSCHMITT K, FLESSA H, GUGGENBERGER G, MATZNER E, MARSCHNER B. SOM fractionation methods: relevance to functional pools and to stabilization mechanisms., 2007, 39(9): 2183-2207.

[21] SIX J, CONANT R T, PAUL E A, PAUSTIAN K. Stabilization mechanisms of soil organic matter: Implications for C-saturation of soils., 2002, 241(2): 155-176.

[22] SOLOMON D, LEHMANN J, KINYANGI J, AMELUNG W, LOBE I, PELL A, RIHA S, NGOZE S, VERCHOT L, MBUGUA D, SKJEMSTAD J, SCHAFER T. Long-term impacts of anthropogenic perturbations on dynamics and speciation of organic carbon in tropical forest and subtropical grassland ecosystems., 2007, 13(2): 511-530.

[23] SIX J, PAUSTIAN K. Aggregate-associated soil organic matter as an ecosystem property and a measurement tool., 2014, 68(1): A4-A9.

[24] DEL GALDO I, SIX J, PERESSOTTI A, COTRUFO M F. Assessing the impact of land-use change on soil C sequestration in agricultural soils by means of organic matter fractionation and stable C isotopes., 2003, 9(8): 1204-1213.

[25] GUILLAUME T, MUHAMMAD D, KUZYAKOV Y. Losses of soil carbon by converting tropical forest to plantations: erosion and decomposition estimated by13C., 2015, 21(9): 3548-3560.

[26] WERTH M, KUZYAKOV Y.13C fractionation in transformations at the interface between roots, microorganisms, and soil: a review and synthesis., 2010, 42(9): 1372-1384.

[27] QIAO N, XU X L, CAO G M, OUYANG H, KUZYAKOV Y. Land use change decreases soil carbon stocks in Tibetan grasslands., 2015, 395(1-2):231-241.

[28] STEWART C E, PAUSTIAN K, CONANT R T, PLANTE A F, SIX J. Soil C saturation: implications for measurable C pool dynamics in long-term incubations., 2009, 41(2): 357-366.

[29] WEI X, SHAO M, GALE W J, ZHANG X, LI L. Dynamics of aggregate-associated organic carbon following conversion of forest to cropland.2013, 57(3): 876-883.

[30] DAVIDSON E A, JANSSENS I A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change., 2006, 440(7081): 165-173.

[31] ZHANG K, DANG H, ZANG Q, CHENG X. Soil carbon dynamics following land-use change varied with temperature and precipitation gradients: evidence from stable isotopes.2015, 21(7): 2762-2772.

[32] SCHWEIZER S A, FISCHER H, HÄRING V, STAHR K. Soil structure breakdown following land use change from forest to maize in Northwest Vietnam., 2017, 166: 10-17.

[33] ROVIRA P, VALLEJO V R. Labile, recalcitrant, and inert organic matter in Mediterranean forest soils., 2007, 39(1): 202-215.

[34] SOLOMON D, LEHMANN J, ZECH W. Land use effects on soil organic matter properties of chromic Luvisols in semi-arid northern Tanzania: carbon, nitrogen, lignin and carbohydrates., 2000, 78(3): 203-213.

[35] PLAZA C, COURTIER-MURIAS D, FERNÁNDEZ J M, POLO A, SIMPSON A J. Physical, chemical, and biochemical mechanisms of soil organic matter stabilization under conservation tillage systems: a central role for microbes and microbial by-products in C sequestration., 2013 57(2): 124-134.

[36] GUNINA A, KUZYAKOV Y. Pathways of litter C by formation of aggregates and SOM density fractions: Implications from13C natural abundance., 2014, 71(6): 95-104.

[37] SIX J, ELLIOTT E T, PAUSTIAN K. Soil macroaggregate turnover and microaggregate formation: a mechanism for C sequestration under no-tillage agriculture., 2000, 32(14): 2099-2103.

[38] SIX J, BOSSUYT H, DEGRYZE S, DENEF K. A history of research on the link between (micro)aggregates, soil biota, and soil organic matter dynamics., 2004, 79(1): 7-31.

[39] KALHORO S A, XU X, CHEN W, HUA R, RAZA S, DING K.

Effects of different land-use systems on soil aggregates: a case study of the Loess Plateau (Northern China)., 2017, 9(8): 1349.

[40] GARCIA-FRANCO N, ALBALADEJO J, ALMAGRO M, MARTINEZ-MENA M. Beneficial effects of reduced tillage and green manure on soil aggregation and stabilization of organic carbon in a Mediterranean agroecosystem., 2015, 153: 66-75.

[41] TORN M S, KLEBER M, ZAVALETA E S, ZHU B, FIELD C B, TRUMBORE S E. A dual isotope approach to isolate carbon pools of different turnover times., 2013, 10(12): 8067-8081.

[42] 蔡岸冬, 徐香茹, 张旭博, 徐明岗, 张文菊. 不同利用方式下土壤矿物结合态有机碳特征与容量分析. 中国农业科学, 2014, 47(21): 4291-4299.

CAI A D, XU X R, ZHANG X B, XU M G, ZHANG W J. Capacity and characteristics of mineral associated soil organic carbon under various land uses.2014, 47(21): 4291-4299. (in Chinese)

[43] BAYER C, MIELNICZUK J, GIASSON E, MARTIN-NETO L, PAVINATO A. Tillage effects on particulate and mineral‐associated organic matter in two tropical brazilian soils.2006, 37(3/4): 389-400.

[44] DENEF K, SIX J, MERCKX R, PAUSTIAN K. Carbon sequestration in microaggregates of no-tillage soils with different clay mineralogy., 2004, 68(6): 1935-1944.

[45] VIRTO I, BARRÉ P, CHENU C. Microaggregation and organic matter storage at the silt-size scale., 2008, 146(1/2): 326-335.

[46] PAUL E A, COLLINS H P, LEAVITT S W. Dynamics of resistant soil carbon of Midwestern agricultural soils measured by naturally occurring14C abundance., 2001, 104(3/4): 239-256.

[47] HOBBIE E A, WERNER R A. Intramolecular, compound-specific, and bulk carbon isotope patterns in C-3 and C-4 plants: a review and synthesis., 2004, 161(2): 371-385.

[48] XU M, LI X, CAI X, GAI J, LI X, CHRISTIE P, ZHANG J. Soil microbial community structure and activity along a montane elevational gradient on the Tibetan Plateau., 2014, 64(6): 6-14.

[49] DUNGAIT J A J, HOPKINS D W, GREGORY A S, WHITMORE A P. Soil organic matter turnover is governed by accessibility not recalcitrance., 2012, 18(6): 1781-1796.

（责任编辑 李云霞）

Impact of Land Use Type on Soil Organic Carbon Fractionation and Turnover in Southeastern Tibet

XU Meng1, 2, LI XiaoLiang3, CAI XiaoBu4, LI XiaoLin2, ZHANG XuBo1, ZHANG JunLing2

(1;2;3;4)

【Objective】Land use conversion from native vegetation to cropland has a great impact on soil organic carbon (SOC) storage. Recently, the area of croplands that converted from native forests or grasslands increases gradually in the southeastern part of Tibetan Plateau, yet its influence on SOC fractionation and turnover remains unknown. It is therefore in great need to understand the extent and mechanisms that difference in land use type has exerted on SOC storage, which will provide scientific basis for sustainable management of agricultural soils in southeastern Tibet. 【Method】In the present study, soil samples were collected from cropland that had been cultivated for more than 50 years, as well as adjacent native forest and grassland. A combination of physical and chemical method was conducted to partition the SOC into different fractions. For each SOC fractions, the present in soil mass, SOC content and stable carbon isotope composition (13C) were analyzed to investigate the differences in SOC fractionation and turnover under different land use types. 【Result】The SOC that stored in top 20 cm soil of cropland was (39.4±2.0) Mg C·hm-2, which was 52% and 45% lower than native forest ((81.5±8.5) Mg C·hm-2) and grassland ((71.4±7.3) Mg C·hm-2), respectively. Compared to native forests and grasslands, long-term cultivation led to a decrease of 63.4%-70.8% in the mass of coarse particulate organic matter (cPOM), whereas the mass of microaggregate (μagg) and easily dispersed silt and clay (dSilt+Clay) increased 10.0%-25.9% and 65.7%-86.2%, respectively. The C content in each SOC fraction was significantly lower in cropland soils than that in forest and grassland soils. SOC content in cropland soils was 51.7%-58.1% lower than that in forest and grassland soils. The SOC contents of unprotected C, physically protected C and biochemically protected C pool reduced 79.8%-86.3%, 72.4%-73.1% and 32.4%-39.8% in cropland soils, respectively, and were positively correlated with changes in total SOC content following land use conversion. The SOC content in chemically protected C pool, however, was not affected by land use conversion. The C/N ratio and13C differed among different SOC fractions and three land use types. The C/N ratio of cPOM in cropland soils (10.0±0.5) was significantly lower than forest soils (13.5±0.4), whereas its13C value ((-21.6±0.5) ‰) was significantly higher than forest soils (-23.6±0.4) ‰. The13C values of microaggregate-protected particulate organic matter (iPOM) and non-hydrolyzable fraction (NH-dSilt+Clay and NH-μSilt+Clay) were comparably lower (-25.3‰- -27.2‰) than other SOC fractions, and had significantly lower C/N ratio in cropland soils (8.4-9.4) compared to forest soils (13.5-15.9). 【Conclusion】Results of the present study indicated that long-term cultivation had resulted in c.a. 50% reduce in SOC stock compared to native vegetation in southeastern Tibet. Agricultural cultivation strongly promoted the turnover of unprotected particulate organic matter, and suppressed the formation of more stabilized SOC, such as microaggregate-protected SOC fraction, which contributed great proportions to the dramatic decrease in SOC storage. Therefore, implement of no-tillage or other protective management would be necessary to reduce the negative influence of agricultural land use on SOC storage and to maintain the sustainable utilization of natural soil sources in southeastern Tibet.

soil organic carbon; land use conversion; croplands; soil organic carbon fractionation; stable carbon isotope; Tibetan Plateau

2018-04-02；

2018-05-25

国家重点基础研究发展计划（2015CB150500）、国家自然科学基金项目（41701333，41461054）

徐梦，Tel：15201029279；E-mail：mengx@cau.edu.cn。 通信作者张旭博，Tel：15210464308；E-mail：zhangxb@igsnrr.ac.cn。通信作者张俊伶，Tel：010-62732574；E-mail：junlingz@cau.edu.cn

10.3864/j.issn.0578-1752.2018.19.009