基于137Cs的青藏高原高寒草甸土壤侵蚀及碳流失估算.pdf

目的定量分析青藏高原高寒草甸土壤侵蚀状况及其伴随的碳流失,为全面评估土壤侵蚀影响,实施有效水土保持措施提供参考。方法结合1 3 7C s示踪技术与前人研究,对青藏高原高寒草甸土壤的整体侵蚀水平及其土壤有机碳流失进行了估算。结果未受人为扰动的高寒草甸土壤自上而下表现出3个层次(A,B和C层)的理化性质特征,其1 3 7C s分布遵循显著指数递减模式。目前,高原草甸土壤年均侵蚀模数约为7 72 3 0t/k m2,推测其每年直接导致的土壤有机碳损失量平均不低于4.8 6t/k m2。结论青藏高原高寒草甸土壤侵蚀水平整体较弱,但因土壤侵蚀流失的有机碳不容忽视。在未来气候变化背景下,升温导致的土壤湿度下降对植被生长的限制,以及人类活动的影响,较大可能成为诱使青藏高原草甸土壤退化和有机碳流失的潜在因素。关键词:1 3 7C s示踪技术;高寒草甸;土壤侵蚀;有机碳;青藏高原文献标识码:A 文章编号:1 0 0 0-2 8 8 X(2 0 2 3)0 3-0 3 3 0-0 8 中图分类号:F 3 0 1.2 4,S 1 5 7.1,X 1 7 1.1文献参数:卢海涯,栗文佳,于世永.基于1 3 7C s的青藏高原高寒草甸土壤侵蚀及碳流失估算J.水土保持通报,2 0 2 3,4 3(3):3 3 0-3 3 7,3 4 8.D O I:1 0.1 3 9 6 1/j.c n k i.s t b c t b.2 0 2 3.0 3.0 3 7;L u H a i y a,L iW e n j i a,Y uS h i y o n g.1 3 7C s-b a s e de s t i m a t i o no fs o i le r o s i o na n do r g a n i cc a r b o nl o s si na l p i n em e a d o ws o i lo nT i b e t a nP l a t e a uJ.B u l l e t i no fS o i l a n dW a t e rC o n s e r v a t i o n,2 0 2 3,4 3(3):3 3 0-3 3 7,3 4 8.1 3 7C s-b a s e dE s t i m a t i o no fS o i lE r o s i o na n dO r g a n i cC a r b o nL o s s i nA l p i n eM e a d o wS o i l o nT i b e t a nP l a t e a uL uH a i y a1,L iW e n j i a2,Y uS h i y o n g1(1.S c h o o l o fG e o g r a p h y,G e o m a t i c sa n dP l a n n i n g,J i a n g s uN o r m a lU n i v e r s i t y,X u z h o u,J i a n g s u2 2 1 1 1 6,C h i n a;2.I n s t i t u t eo fT i b e t a nP l a t e a uR e s e a r c h,C h i n e s eA c a d e m yo fS c i e n c e s,S t a t eK e yL a b o r a t o r yo fT i b e t a nP l a t e a uE a r t hS y s t e m,B e i j i n g1 0 1 4 0 8,C h i n a)A b s t r a c t:O b j e c t i v eS o i l e r o s i o na n d i t sa s s o c i a t e dc a r b o n l o s s i na na l p i n em e a d o wo nt h eT i b e t a nP l a t e a uw e r eq u a n t i t a t i v e l ya n a l y z e d i no r d e r t op r o v i d ea r e f e r e n c e f o r c o m p r e h e n s i v e l ye v a l u a t i n gt h e i m p a c to f s o i le r o s i o na n d f o r i m p l e m e n t i n ge f f e c t i v es o i l a n dw a t e rc o n s e r v a t i o nm e a s u r e s.M e t h o d sW eu s e d1 3 7C st r a c i n gt e c h n o l o g yc o m b i n e dw i t ht h er e s u l t so fp r e v i o u sr e s e a r c hs t u d i e st od e t e r m i n ee r o s i o nr a t ea n dc o n s e q u e n tc a r b o nl o s sf o ra nu n d i s t u r b e da l p i n e m e a d o w s o i lo nt h e T i b e t a nP l a t e a u.R e s u l t sT h eu n d i s t u r b e dm e a d o ws o i l e x h i b i t e dt h r e ed i s t i n c t l a y e r s(A,B,a n dC)o f p h y s i o c h e m i c a l p r o p e r t i e s,a n d t h ed i s t r i b u t i o no f1 3 7C sd e c l i n e de x p o n e n t i a l l yw i t hi n c r e a s i n gs o i ld e p t h.T h er e c e n te r o s i o nr a t eo f t h ea l p i n em e a d o ws o i lw a sa b o u t 7 72 3 0t/(k m2y r),a n dw a s e s t i m a t e d t o r e s u l t i na na v e r a g e a n n u a l s o i l o r g a n i c c a r b o n l o s so fn o t l e s s t h a n4.8 6t/k m2.C o n c l u s i o nE v e n t h o u g ht h eo v e r a l l s o i l e r o s i o nr a t eo f t h e a l p i n em e a d o wo n t h eT i b e t a nP l a t e a uw a s l o w,i tso r g a n i cc a r b o nl o s sd u et oe r o s i o nc a n n o tb eu n d e r e s t i m a t e d.B o t hi n c r e a s e dh u m a na c t i v i t i e sa n dr e d u c e dm e a d o wg r o w t hd u e t od e c r e a s e ds o i lm o i s t u r eu n d e r f u t u r eg l o b a lw a r m i n ga r el i k e l yt or e s u l t i nf u r t h e rs o i l d e g r a d a t i o na n ds o i l o r g a n i cc a r b o nl o s s i na l p i n em e a d o w s.K e y w o r d s:1 3 7C s t e c h n o l o g y;a l p i n em e a d o w;s o i l e r o s i o n;o r g a n i cc a r b o n;T i b e t a nP l a t e a u 受自然环境变化及人类活动的影响,土壤侵蚀已成为了威胁着全球土壤安全的严重问题。研究显示,土壤侵蚀造成的土壤退化导致了大幅度的土壤可耕作深度缩减、持水能力下降及粮食作物减产1-2。同时,土壤作为陆地上最活跃、最庞大的碳库3,土壤侵蚀过程也将导致土壤内部的有机碳流失,影响全球碳平衡。因此,如何准确估计土壤的侵蚀状况及其导致的碳流失过程,对评价区域性生态系统服务效能变化及推动科学有效治理至关重要。当前,已有多种手段方法被开发应用于不同时空尺度的土壤侵蚀研究。如针对小尺度区域的侵蚀小区测量4或运用侵蚀图进行大范围估算。但上述方法通常存在反映时间尺度短,或当空间尺度不同时外推精度难以保证等问题。1 3 7C s作为一种人工核素,形成于2 0世纪5 0年代至7 0年代的大气核试验。大气中1 3 7C s主要随降水沉降,并迅速被土壤表层吸收,其后续再分配多取决于土壤的空间流动和垂直迁移,基本不受植物吸附作用或淋溶过程影响5。因此,通过对比分析研究区域土壤的1 3 7C s实际含量与背景值,即可获取样区的自1 3 7C s沉降以来的侵蚀或堆积信息。相较于需要建立长期观测的传统方法,1 3 7C s示踪技术能在短期采样中追溯土壤的侵蚀过程,因而更加便捷并被广泛应用6-8。作为“亚洲水塔”和“世界第三极”地区,青藏高原土壤温度较低,有机质分解缓慢9,在全球碳循环中发挥着重要的碳汇作用1 0。高寒草原(包括草甸和干草原)是青藏高原最主要的生态系统,覆盖了约2/3的面积。其中,高寒草甸土壤在草原类型中有机碳含量最高,对高原碳汇具有重要贡献。但高原的特殊环境使高寒草甸土壤易受风蚀、水蚀和冻融侵蚀的风险。长期侵蚀会造成土壤中C和N等元素流失,降低土壤质量。然而,由于缺乏直接测量,对高寒草甸土壤侵蚀导致的有机碳流失量化评估较少。另外,近几十年青藏高原气温以远高于全球平均的速率上升(0.3/1 0a)1 1。这预计会改变高寒草甸土壤侵蚀环境、加速土壤有机碳(S O C)分解,进而影响高原草甸土壤的碳汇功能1 2。因此,对高寒草甸土壤侵蚀和碳流失状况的估计,也可以为准确评估气候变化背景下的青藏高原土壤动态提供参考性信息。本研究目的在于根据高寒草甸土壤的1 3 7C s和S O C对其当前的土壤侵蚀状况及其土壤有机碳损失或累积状况进行评估,并尝试讨论在预测气候变化和人类活动影响背景下青藏高原高寒草甸土壤可能的侵蚀和碳流失潜力响应。1 研究区概况青藏高原位于中国西南部,平均海拔约43 8 0m,面积约2.6 01 06k m21 3,是世界上最大的高原。位于2 6 3 9 4 7 N,7 3 1 9 1 0 4 4 7 E之 间,包 括 了中国西藏的全部以及青海、新疆、甘肃、四川和云南的部分地区。海拔极高使高原内气候与周围迥异,形成了辐射强、日照多,但少积温、日较差大的独特高原季风 气候1 4。年 均 温 约 为1.6,年 降 水 约 为4 1 6mm1 5。高原内气候具有明显的东南西北向梯度变化,因此由南至北,高原上依次分布着森林、灌丛、草地和荒漠等植被生态1 6。高寒草原(包括草甸和干草原)为高原上最主要的植被类型,面积约占6 0%,其中高寒草甸约占2/5。本文研究样地位于青海省玉树市沱沱河南面约5.5k m的山地地区,海拔45 8 9m,平均坡度1 0。年均温约-3.8 6,年降水量4 4 0mm。研究区域气候为典型高原大陆性气候,冷暖和干湿两季分明。植被类型属 于 典 型 高 寒 草 甸,土 壤 为 草 甸 土。以 嵩 草(C a r e xm y o s u r o i d e s)、团垫黄耆(A s t r a g a l u s a r n o l-d i i)和 丛 生 钉 柱 委 陵 菜(P o t e n t i l l as a u n d e r s i a n av a r.c a e s p i t o s a)为群落优势种,伴随着白火绒草(L e-o n t o p o d i u mo c h r o l e u c u m)和矮生二裂委陵菜(P o-t e n t i l l ab i f u r c av a r.h u m i l i o r)等散生草本植物,植被覆盖度9 0%,裸地面积比例约占2%1 7。地表无凋落物、砾石、覆沙等,且无人类活动利用,但存在多处鼠、兔等啮齿类动物洞穴。2 材料和方法2.1 样品采集土壤样品于2 0 2 0年8月采集。在样地内沿顺坡方向设置1 0m1 0m样方。样方内沿对角和中心位置设置3个1m1m小样方,各小样方水平间隔5m左右。各小样方内在随机位置使用直径4 0mm土钻钻取一根钻孔,采集深度05 0c m。由于样地内裸露地段多分布有动物洞穴,样方内土壤样品钻取通常避开裸露位置(图1)。受实际地形和植被限制,样区内无法判断是否存在无侵蚀的平坦地段,未进行背景值采样,后续将通过模型计算并结合前人研究结果对样区背景值进行估算和验证。同时,样区内调查与记录植被组成和盖度等信息,并使用G P S定位仪器记录采样位置的地理坐标与海拔高度。133第3期 卢海涯等:基于1 3 7C s的青藏高原高寒草甸土壤侵蚀及碳流失估算 注:图a数据来源于青藏高原植被图(2 0 2 0)1 7、国家青藏高原科学数据中心,并基于自然资源部标准地图服务系统G S(2 0 2 0)4 6 3 4号标准地图校准制作。黑点为本文采样区域;橙点为其他研究的高寒草甸土壤1 3 7C s采样位置1 8-2 3。图1 青藏高原植被类型及高寒草甸样点位置F i g.1 V e g e t a t i o nt y p e sa n da l p i n em e a d o ws a m p l e s i t e so nT i b e t a nP l a t e a u2.2 试验测量与分析2.2.1 土壤理化性质与C,N测量 将钻孔样品按2c m为间隔进行分割,记录各层土样重量,并于5 0烘箱内烘干2 4h,再次称重,用于计算土壤含水量和土壤容重。土壤C,N使用EMA 3 0 0 0元素分析仪(T h e r m oF i s h e rS c i e n t i f i c,U S A)测量。样品烘干后 研磨过1 0 0目筛,取1/2上机测进行C,N测量,得到土壤总碳、总氮。另外1/2加入过量HC l,放置反应2 4h去除碳酸盐,加入蒸馏水不断离心直至样品呈中性(p H7)。于5 0烘干4 8h,研磨过1 0 0目筛网后,再次测量得到土壤有机碳。实验室标准重复分析典型标准差为0.1。2.2.2 1 3 7C s测量和计算 04 0c m土层的土样以24c m为间隔按深度混合土样,剩余4 05 0c m为组样,在5 0烘箱中放置约2 4h。烘干后,过2mm筛,去除植物残体和砾石。研磨过0.0 7 4mm筛网,混合均匀后填充入长7c m和直径1.5c m的聚氯乙烯标准管中压实进行1 3 7C s测量。采用美国Am e t e k公司的O R T E C高 纯 锗能 谱 测 量 系 统(H P G e,D e s p e c),在峰值分辨率(FWHM)为6 6 2k e V时检测土样的1 3 7C s活度,测试累积时间大于4 00 0 0秒,误差控制标准为5%。直接测量结果为单位质量上的活度(B q/k g),需要将 其 转 化 为 基 于 面 积 的 样 点1 3 7C s库 存 容 量(B q/m2),公式如下2 4:C P I=1 03ni=1CiBiDi(1)式中:C P I为样点的1 3 7C s库存容量(B q/m2);i为样土的土层编号;n为样土层数;Ci为i层的1 3 7C s活度(B q/k g);Bi为i层的土壤容重(g/c m3);Di为i层的土壤厚度(m)。所有试验于2 0 2 2年1月至6月开始和完成。试验数据运用E x c e l和O r i g i n2 0 1 8进行处理与分析,图形绘制使用A r c G I S1 0.2和O r i g i n2 0 1 8软件。2.3 1 3 7C s背景值确定确定1 3 7C s背景值是利用1 3 7C s技术准确计算样点土壤侵蚀速率的基础。理论上,应当选择样点周围地形平坦,地表植被覆盖充足,且无侵蚀和堆积作用的地段作为计算1 3 7C s背景值的区域2 5。但实际过程中,受复杂骤变的地形地势、局地小气候和植物演替影响,限制的空间范围内这种理想的1 3 7C s背景地块很难满足和确定。在缺少理想1 3 7C s背景值地段选择的情况下,模型计算是一种较好的替代性选择。前人研究表明,结合区域样点的1 3 7C s数据集和相关模型,通过对1 3 7C s分布模式函数进行模型拟合,即可得到区域的1 3 7C s背 景 值 及 其 剖 面 形 状 参 数。拟 合 函数2 6-2 7如下:A(x)=Ar e f(1-e-x/h0)(2)式中:x为质量深度(k g/m2);A(x)为质量深度x对应的1 3 7C s累计总量(B q/m2);Ar e f为样区的1 3 7C s背景值(B q/m2);h0为1 3 7C s的剖面形状参数。2.4 1 3 7C s土壤侵蚀转换模型利用1 3 7C s技术估算土壤样点的侵蚀速率,关键在于建立1 3 7C s损失量或累积量与土壤空间迁移和流233 水土保持通报 第4 3卷动间的联系。在没有或较少人为扰动土壤中,1 3 7C s的分布通常为随深度指数减小。因此,2 0世纪9 0年代提出了一种典型的剖面形状模型用以描述这种特征2 4,并被广泛应用于各地区未受干扰土地的土壤侵蚀评价。然而,1 3 7C s对土壤黏粒具有依附性,其后续重新分布还会受到土壤细粒垂直运动影响。由于上述模型没有考虑这一过程导致的1 3 7C s损失,可能存在土壤侵蚀的高估。基于此,Z h a n g等在W a l l i n g和H e的侵蚀模型2 8基 础 上,发 展 了 一 种 简 化 的1 3 7C s剖面迁移模型(S TM)2 9,有效减少了计算中的冗余和误差,提高了未扰动土地土壤侵蚀评价的精度。同时已有研究表明,剖面迁移模型在青藏高原和云贵高原山地地区具有较好适用性2 0-2 1,2 3,3 0。因此,本研究采用S TM转换模型计算土壤侵蚀速率,相关公式如下:K=1-Ti=1e r f cH2Di(3)式中:H为1 9 6 3年以来的年均侵蚀深度(c m/a);T为采样年份与1 9 6 3年的差值,本研究中值为5 7;D为1 3 7C s扩散系数;K为1 3 7C s损失比例;e r f c为补余误差函数。其中,1 3 7C s损失比例K值计算公式如下:K=Ar e f-AAr e f1 0 0%(4)式中:Ar e f为1 3 7C s背景值;A为样点的C P I值。扩散系数D通过1 3 7C s的深度分布数据集进行模型拟合得到,拟合公式2 9如下:C(z,T)=C(0,T)e-z2/4D T(5)式中:z为土层深度(c m);C(z,T)为T对应年份土壤z深度处的1 3 7C s活度含量(B q/k g);C(0,T)即T年份的土壤表层1 3 7C s活度含量;为半衰期常数(=0.0 2 3)。土壤年均侵蚀深度H通过E x c e l软件对e r f c函数(补余误差积分函数)进行单变量求解得出,通过以下公式转换为侵蚀模数单位Et/(k m2a):E=HB1 04(6)式中:B为样点的土壤容重(g/c m3);1 04为单位转换常量。3 结果与分析3.1 高寒草甸土壤理化特征高寒草甸样点剖面的土壤含水量、干密度、有机碳(TO C)和总碳(T N)分布见图2。结果显示,高寒草甸土壤剖面的物理和化学性质随深度呈现出规律变化,凸显了3个主要的土壤层次特征(A,B和C层)。具体上看,A层(01 0c m)土壤含水量最高,在68c m深度时出现了最大峰值,达5 3.1 0%。同时在A层,土壤干密度随着深度的增加而逐渐上升,至A,B层交界时达0.9 7g/c m3。T N和TO C含量在A层最高,均表现出随深度减少的规律。TO C在土壤表面最高,到A层6c m以下大幅降低(4 9.1 4%);T N最高峰出现在6c m处而非土壤表面,之后也随深度增加大幅减少(5 2.6 2%)。B层中4个方面的性质相对稳定,除土壤容重呈现了细微的波动增加外,土壤含水量、T N和TO C基本稳定。至C层,各深度土层的理化性质再次出现了显著的差异。土壤含水量随深度 不 断 减 小,至4 5c m处 减 小 到 剖 面 最 低 值1 3.9 7%,之后有所增加;土壤干密度波动上升,4 5c m左右最大。土壤T N和T O C在C层中呈现出了比较一致的特征,均随深度呈指数减少趋势。在4 2c m处T N和T O C含量最低,分别为0.0 7,1.8 0g/k g。图2 高寒草甸土壤剖面理化性质F i g.2 P h y s i c o c h e m i c a lp r o p e r t i e so f s o i l p r o f i l eo fa l p i n em e a d o ws o i l3.2 1 3 7C s背景值运用O r i g i n2 0 1 8软件,结合样点的1 3 7C s分布数据集,对公式2进行了拟合得到样点1 3 7C s背景值范围为1 7 1 3.4 71 0 8.2 6B q/m2。为了验证模拟的1 3 7C s背景值的可靠性,与前人研究在邻近的高寒草甸测量得到的1 3 7C s背景值进行了对比(表1)。整体上看,高寒草甸样点的1 3 7C s背景值符合K.T等3 1对亚洲3 0 4 0 N范围内的统计。邵全琴等2 1在玉树县样点(3 2 5 8 N,9 6 1 9 E)测 量 得 到 背 景 值 为16 1 6.7 6B q/m2,在玛多县(3 4 2 7 N,9 7 5 8 E)测得值 为19 2 6.4 5 B q/m2;严 平 等3 2在 五 道 梁 地 区(3 5 0 N,9 2 3 0 E)的测量背景值为14 0 0.9 0B q/m2。333第3期 卢海涯等:基于1 3 7C s的青藏高原高寒草甸土壤侵蚀及碳流失估算由于降水和纬度位置的影响,1 3 7C s参考值在青藏高原上的分布基本符合南部低,中北部高,且随降水呈自东向西减少的趋势3 3。鉴于本研究区在青藏高原中部,位于上述区域之间。本文认为将1 7 1 3.4 7B q/m2作为研究区的1 3 7C s背景值具有一定的可行性和合理性。表1 本文和其他研究高寒草甸的1 3 7C s背景值对比T a b l e1 C o m p a r i s o no f1 3 7C s r e f e r e n c ev a l u e sb e t w e e nt h i ss t u d ya n do t h e r s t u d i e so na l p i n em e a d o w样点位置 降水量/mm1 3 7C s背景值/(B qm-2)文献来源3 4 0 9 N,9 3 3 1 E4 4 017 1 3.4 71 0 8.2 6本文3 2 5 8 N,9 6 1 9 E4 7 916 1 6.7 62 13 4 2 7 N,9 7 5 8 E5 0 819 2 6.4 52 13 5 0 N,9 2 3 0 E2 6 514 0 0.9 03 23 0 4 0 N(亚洲)8 2 1.0 035 6 3.0 03 13.3 高寒草甸土壤1 3 7C s垂直分布1 3 7C s测量结果表明,研究样区剖面1 3 7C s含量具有显著差异(图3)。1 3 7C s含量最大值出现在土壤48c m,随着深度增加,各土层1 3 7C s含量呈逐渐减少趋势。到达2 0c m左右深度时,1 3 7C s含量接近检测下限0.7 7B q/k g(图3 a)。为了验证样区是否符合典型的未利用高寒草甸土壤,通过O r i g i n2 0 1 8软件对样地剖面各层次土壤的1 3 7C s进行指数模型拟合(参考公式5)。由于04c m土壤表层1 3 7C s活度含量较低,说明此深度土层已遭受侵蚀,无法表征原本的1 3 7C s特征,指数拟合时不参与拟合过程。结果表明,样区剖面各土层的1 3 7C s随深度呈现了极好的指数递减对应关系(图3 b)。这说明,该样区符合典型的未扰动高寒草甸区域,其剖面的1 3 7C s分布具有较好代表性。同时,拟合关系显示,在土壤最表层(04c m),1 3 7C s模拟活度值为5 3.6 3B q/k g,远大于其实测值5.7 3B q/k g,且样点的1 3 7C s总累计量为16 2 1.9 7B q/m2,略小于样区的背景值。这说明,该样区高寒草甸土壤存在着一定的侵蚀过程。图3 1 3 7C s剖面分布及指数模型拟合F i g.3 P r o f i l ed i s t r i b u t i o na n de x p o n e n t i a lm o d e l f i t t i n go f1 3 7C s3.4 青藏高原高寒草甸土壤侵蚀状况S TM模型计算结果表明,本研究的高寒草甸样区土壤平均侵蚀模数为7 7.1 42 7.5 0t/k m2/a。与前人的研究进行对比(图4),本文结果较其他研究区小。推测可能是由于研究区域的位置和地形、采样时间及模型选择差异导致。不同地区高寒草甸侵蚀因素规模和特征不同。该研究样区属于山地坡面高寒草甸区,风力与降雨形成的坡面侵蚀是主要土壤侵蚀因素,而前人研究如L iY.S.1 8,W a n gY i b o1 9和陈泽涛等2 2的样地中包括了河盆和沟谷部分,还受河流冲蚀或堆积的影响。同时,根据采样时间与对应侵蚀结果的对比发现,随着采样年份增加,高寒草甸土壤年均侵蚀速率整体具有减小趋势。主要有两个方面的原因:一是由于1 3 7C s示踪技术反映的是自1 3 7C s沉降以来的平均土壤再分配信息;二是近年来高原气候、生态因素的改善。另外,W a n gY i b o等采用的1 3 7C s转换模型为线性估计,在侵蚀计算中存在高估的可能性较大2 9。然而,虽然受限于上述因素影响,不同高寒草甸区的研究结果具有一定差异,但整体上看,青藏高原高寒草甸土壤类型的平均侵蚀模数都普遍较低。根据中国土壤侵蚀等级和分类标准,各高寒草甸地区均属于轻微侵蚀0.0 5)。因此,结合上述研究估计,青藏高原高寒草甸土壤S O C平均 含 量 在03 0,3 05 0c m内 分 别 为2 0.1 73 5.4 9,9.9 71 3.3 9g/k g。图5 高寒草甸不同土层S O C含量和研究结果对比F i g.5 S O Cc o n t r i b u t i o n i nd i f f e r e n t l a y e r so fa l p i n em e a d o ws o i l a n dr e s u l t s c o m p a r i s o n4 讨论与结论4.1 讨 论(1)青藏高原高寒草甸土壤侵蚀S O C流失。青藏高原覆盖了近5 8.1 2%的草地,草原生态的有机碳含量基本代表了高原土壤有机碳储量水平。高寒草甸区域降水较多且温度低,为草甸植物的生长繁殖提供了条件,促进了有机质的输入。同时低温环境抑制了土壤有机质分解和矿化,使得内部有机碳不断积累。因此,相较于其他草地类型,高寒草甸土壤拥有较高的有机质和有机碳,在高原碳汇中发挥着不容替代的作用。相较于深层土壤,土壤表层作为侵蚀的直接发生地和作用对象,长期侵蚀导致的土壤空间再分配过程,不仅会导致其土壤结构或性质的变化,改变土壤内的氧化还原过程,也会造成土壤内部有机物质和营养元素的流失或因暴露地表而被风化分解。通过高寒草甸土壤侵蚀模数及其S O C平均含量的乘积可以估算得出,青藏高原高寒草甸每年因土壤侵蚀直接造成的S O C流失量约为1.5 68.1 6t/k m2。据统计,青藏高原高寒草甸面积约5 3.4 71 04k m23 6,这533第3期 卢海涯等:基于1 3 7C s的青藏高原高寒草甸土壤侵蚀及碳流失估算意味着即使保持最低侵蚀水平,高原内高寒草甸土壤每年至少将流失8.3 41 05tS O C。作为全球重要碳库和敏感区域,大规模的S O C流失不仅严重影响着高原内高寒草甸生态的碳平衡过程,对全球气候变化过程也是极大的潜在风险。因此,对于高原草甸植被区域,采取针对性的保护举措、减缓或避免土壤侵蚀过程及S O C损失,对稳固和提高高寒草甸生态的碳汇能力、缓解全球C O2压力相当必要。(2)气候变化及人类活动对高原侵蚀和S O C流失的影响。青藏高原土壤浅薄,植物矮小且生长期短,拥有着仅次于极地的广泛冰川和冻土区域,使其成为对全球气候变化最敏感的地区之一3 8。近几十年,高原气候主要表现为表面温度增加、降水的空间格局变化和风力的年际变化。多数研究表明,这些变化对高原土壤侵蚀和固碳过程起着积极促进的作用。比如,近几十年来青藏高原表面风速整体减小,使高原土壤侵蚀强度以1 4t/k m2/a的速率逐渐减弱3 9;气候暖湿化促使高原草甸草原植物的返青期、生长时长和开花期提前4 0,提高了草原生态生长和有机质输入,间接削弱了土壤受到的侵蚀。虽然气候变化对青藏高原土壤侵蚀和有机碳储存带来许多可预期的好处。但升温导致的土壤蒸发量增加可能严重限制土壤水分有效性,反而增强土壤的可蚀性4 1,促使草原植被发生退化4 2。且相较于草原和灌丛植被,高寒草甸种群要求更高土壤湿度,土壤干化更易于导致草甸群落向中干草本群落转化。另外,放牧和耕作等人为活动也会导致高寒草甸退化,催使土壤侵蚀规模增加,促使S O C转化速率提升且储量减小4 3。同时,S O C转化和降解将向大气释放C O2,对高原气候变化形成正反馈4 4。相反,减少放牧、休耕还草和围栏护草等措施能直接增加S O C密度,推动草地植被恢复。因此,对高寒草甸区域进行持续侵蚀和碳流动监测、推广有效管理策略,对于维持和提高高寒草甸土壤固碳潜力及更好了解气候变化影响至关重要。4.2 结 论高寒草甸土壤是青藏高原重要土壤类型,在高原碳汇中发挥着重要作用。本文选择了青藏高原中部的典型未扰动的高寒草甸土壤,对其1 3 7C s和土壤有机碳测量进行了分析。结果表明,研究样区高寒草甸土的理化性质随土层有规律变化,1 3 7C s分布遵循显著指数递减模式,符合典型未扰动草甸土壤。通过对比本文和前人不同的地区研究发现,青藏高原草甸土壤整体侵蚀水平属于轻微侵蚀,估计其目前的年均侵蚀速率约为7 72 3 0t/(k m2a-1)。同时,根据土壤侵蚀速率及高寒草甸土壤上层S O C含量计算得出,青藏高原高寒草甸区由侵蚀直接导致的S O C损失量不容小觑,其年均流失量不少于8.3 21 05t。另外,在全球持续变暖趋势下,升温带来的土壤水分限制、放牧等人类活动可能是推动高寒草甸植被退化、削减土壤固碳能力,并形成正反馈机制的潜在威胁。对青藏高原高寒草甸进行有效监测和管理,防止其生态系统遭受侵蚀、退化,对保持高原碳汇功能和生态安全具有十分重要的意义。参考文献1 R o mm e n sT,V e r s t r a e t e nG,L a n gA,e t a l.S o i l e r o s i o na n ds e d i m e n td e p o s i t i o ni nt h eB e l g i a no e s sb e l td u r i n gt h e H o l o c e n e:E s t a b l i s h i n g as e d i m e n tb u d g e tf o ras m a l l a g r i c u l t u r a lc a t c h m e n tJ.T h eH o l o c e n e,2 0 0 5,1 5(7):1 0 3 2-1 0 4 3.2 B a k k e rM M,G o v e r sG,J o n e sRA,e t a l.T h e e f f e c t o fs o i l e r o s i o no nE u r o p esc r o py i e l d sJ.E c o s y s t e m s,2 0 0 7,1 0(7):1 2 0 9-1 2 1 9.3 Q i nY,Y i SH,L iNJ,e t a l.A d v a n c e i ns t u d i e so f c a r-b o nc y c l i n go na l p i n eg r a s s l a n d so ft h eQ i n g h a i-T i b e t a nP l a t e a uJ.A c t aP r a t a c u l t u r a eS i n i c a,2 0 1 2,2 1(6):2 7 5-2 8 5.4 B o a r d m a nJ,F a v i s-M o r t l o c k D,F o s t e rI.A 1 3-y e a rr e c o r do fe r o s i o no nb a d l a n ds i t e si nt h eK a r o o,S o u t hA f r i c aJ.E a r t h S u r f a c eP r o c e s s e sa n d L a n d f o r m s,2 0 1 5,4 0(1 4):1 9 6 4-1 9 8 1.5 V a n d e nB y g a a r tAJ,P r o t zR.B o m b-f a l l o u t1 3 7C sa sam a r k e ro fg e o m o r p h i cs t a b i l i t yi nd u n es a n d sa n d