纳帕海湿地区表土土壤有机碳的活性水平及空间变异

　　全球陆地生态系统 0 ～100 cm 土层土壤有机碳(SOC)储量约为 1 462 ～1 548 Gt，分别是大气碳库和陆地生物量碳库的 2 ～3 倍。其中，湿地生态系统 SOC 库约占 20%以上，是极为重要的碳汇和碳源。高海拔或高纬寒区生态系统是生态响应敏感区，而寒区湿地生态系统 SOC 库更易受气候变化和日益增强的人类活动的影响而受关注。我国西部高原山地区发育有大面积的高寒湿地，因地处流域上游或源头区，高寒湿地 SOC的变化不仅影响自身的生态过程和功能，还影响其与大气和下游河湖系统等的碳交换。西部高寒湿地SOC 库及其活性组分的时空分异、动态变化及其与环境因子和人类活动等的关系研究越来越受到重视。

　　滇西北横断山中段纳帕海湿地是较独特的低纬高寒湿地，1984 年被列为省级自然保护区，2004 年被列入 Ｒamsar 名录。尽管保护受到重视，但在长期持续人类活动驱动下，纳帕海湿地水文情势、植被生态和景观格局等都发生了显着变化，湿地退化严重。针对纳帕海湿地土壤已开展不少研究，初步揭示了湿地区土壤的特征、退化状态及影响因子，发现微地貌控制下的水分和植被生态耦合分异对湿地区 SOC 以及氮磷等组分等空间分异产生显着影响。然而，高强度且极为频繁的干扰，如近万头大家畜的过载散养、日益增多的游客骑马观光等，对湿地区植被和土壤的影响在空间上是无序的。已有基于湿地区不同退化状态、植物群落或土地利用类型布设典型样带或样地的采样研究，难以反映湿地区整体的土壤特征，也不能反映这些无序干扰所产生的影响。

　　土壤有机碳活性组分又称活性有机碳(LOC –Labile OC)，是土壤碳库中最活跃、周转最快、易被微生物利用、且对环境变化和土地利用响应敏感的组分。根据相应的测试方法，LOC 被定义为轻组有机碳、颗粒有机碳、易氧化有机碳(EOC)、溶解有机碳(DOC)、微生物生物量碳等。尽管不同测试方法得到的 LOC 组分含量差异显着，但 SOC 与各活性组分以及活性组分之间都呈(极)显着正相关，因此，土壤 LOC 组分常被用来估算土壤 SOC 库的可矿化有机碳量、预测 SOC 库变化、评估气候变化和人类活动对 SOC 库的影响，以及利用 LOC/SOC(%)评价土壤有机碳库的活性水平等，表征不同的生态效应。本研究在纳帕海湿地区采用网格均匀布点采集表土(0 ～ 20 cm) 土样，测试其SOC、EOC、DOC 含量，应用地统计学等方法，揭示SOC、EOC、DOC、EOC / SOC ( % ) 、DOC / SOC ( % ) 空间变异的整体特征。

　　1 研究区概况与研究方法

　　1. 1 研究区概况

　　纳帕海湿地为典型的高原湖沼洪泛湿地。湿地区十分平坦，但微地貌分异仍较显着，总体上为东南部和南部向北部和西北部、两侧向中间逐渐降低。湿地区景观类型丰富多样性;土壤为典型亚高山土属。纳帕海湿地所在流域区为典型的季风气候区，季风气候、湿地区的微地貌分异共同决定着湿地区水文生态情势的空间分异及其显着的季节性和年际变化，并影响湿地区土壤相关特性的分异。此外，湿地区还存在多方面有序和无序的人为活动，如道路和水利工程建设、放牧、旅游等，也对湿地区的水文生态情势和土壤产生明显的影响。

　　1. 2 土壤采集和分析

　　本次研究将湿地区划分成 1 km ×1 km 的网格(图 1)。基于各网格内景观特征(地貌、植物群落、土壤湿度等)的实地调查，根据典型性和代表性，在每个网格内布设 3 ～ 5 个采样点采集表层(0 ～ 20cm)土样，每个样点土样为 3 个土样的混合样，每个样点均利用 GPS 记录采样点位置信息。如网格内景观差异显着，适当增设采样点。北部主湖区因淹水难以通行而未采样。本研究共采集湿地区 141 个表土样品。【图1】

　　土样自然风干、去杂、磨碎，过 2 mm、1 mm、0. 5mm 和 0. 149 mm 筛，送中科院东北地理与农业生态研究所分析测试中心测试 SOC、EOC、DOC 等含量。SOC 采用重铬酸钾测定法 ( LY / T1237 － 1999 ) 。

　　EOC 采用高锰酸钾(333 mmol / l) 氧化法。DOC采用水溶萃取和总有机碳分析仪测定，即称取 10 克风干置于 100 ml 塑料瓶中，加入 50 ml 蒸馏水常温振荡 30 min，振荡后的样品以 4 000 转/min 离心 5min，取离心样品的上清液过 0. 45 um 滤膜，同时做三份空白，利用总有机碳分析仪测定滤液有机氮含量。
　　
　　1. 3 地统计分析
　　
　　土壤是时空连续的变异体，具有高度空间异质性，地统计分析与 GIS 结合是研究土壤空间异质性的常用手段。本研究利用 Arcgis 系统的 GeostatisticalAnalyst 模块的克里金插值功能(Kriging)来分析。

　　1． 3． 1 正态分布检验

　　利用偏度和峰度系数、K － S 检验的渐进型分布显着性水平(双尾)进行判识。当样本数据平均值与中值近似相等，偏度系数接近于 0 时、峰度系数绝对值 ＜ 3、K － S 渐进型分布检验显着性水平 ＞0. 05，则基本符合正态分布。若不符合正态分布，则进行数据变换(如对数转换)再检验。土壤是时空变异体，因土壤自身可能存在的显着变异、样地选择的人为性、采样时间的不同、实验测定条件的改变以及系统误差等，都可导致土壤测试数据存在一定数量的离群值，需剔除离群值开展正态分布检验。

　　本研究中，离群值利用“均值 ± 2 个标准差”进行判识，其中 SOC、EOC、DOC、EOC/SOC 分别剔除 12、12、1、4 个离群值，占总样本量比例分别为 8. 51% 、8. 51% 、0. 71% 、2. 84% ，不影响总样本的正态分布检验。基于 SPSS 频数分布统计和单变量 K － S 渐进分布的检验结果见表 1。【表1】

　　1． 3． 2 空间相关性分析

　　利用各变量半方差函数最优拟合模型计算得到块金值 C0、基台值 C0+ C、变程 a 等参数来检验。

　　土壤属性空间变异多表现为各向异性，在克里金插值前，首先选择高斯模型、指数模型等模型，并考虑各向异性，对各变量半方差函数进行拟合，获得相关特征参数及决定系数(Ｒ2)和残差(ＲSS)。Ｒ2越大、ＲSS 越小，相关模型拟合效果越好，进而确定各变量半方差函数最优拟合模型，结果见表 2。【表2】

　　半方差函数分析中，C0/(C0+ C)(基底效应系数)被用来判识变量空间相关性强弱。从表 2来看，DOC/SOC 的空间相关性强，其他变量为中等空间相关，可进行空间克里格插值。表 2 还给出了各变量基于样本数据的变程范围(Minor． a ～ Major． a)，它表征各变量空间相关性的各向异性的显着程度，当 Minor． a 和 Major． a 接近，变量空间相关性的各向异性不显着。表 2 所示，DOC、DOC / SOC 变程范围小，接近各向同性，其他变量各向异性特征较显着。尽管 DOC、DOC/SOC 接近于各向同性，但通常土壤属性的空间相关性表现为各向异性，所以后续克里金插值均考虑各向异性情形。

　　1． 3. 3 克里金插值

　　考虑各变量空间相关性的各向异性，DOC/SOC选择普通克里金插值、其他变量选择对数正态克里金插值，结果图 2 所示。克里金插值交叉检验结果:

　　各变量平均标准误接近于均方根误差，两者比值变化于 1. 001 ～1. 116，且均方根标准误变化于 0. 966～ 1. 112，表明各变量克里金插值结果较理想。
　　
　　2 结果分析

　　2. 1 有机碳及活性组分空间变异特征

　　从表 1 来看，湿地区表土 SOC、EOC、DOC 含量均值分别为 51. 596 g/kg、4. 950 g/kg、155. 360 mg/kg。SOC、EOC 为强变异水平( CV ＞ 100% ) ，DOC 属中等变异水平(CV 介于 10% ～100%)，三个变量含量的极值比分别达到 82. 2、60. 5、56. 4。SOC ～EOC、SOC ～ DOC、EOC ～ DOC 的相关系数分别为0. 929、0. 741、0. 760 ( 通过 0. 001 显着性检验，双尾)。可见，湿地区表土有机碳和活性组分之间，以及不同的有机碳活性组分之间，均为极显着相关。

　　湿地区表土 EOC/SOC(%)、DOC/SOC(%)均值分别为 10. 388%、0. 403%，两者占土壤有机碳的比例并不高。两个变量的 CV 均介于 10% ～ 50%，属中等变异水平。但两个变量的极值比仍然较大，分别达到 16. 5、26. 2，可见湿地区表土有机碳活性组分的差异还是比较明显。

　　各变量变异系数和极值比体现了湿地区显着的水文生态分异和强的人为活动扰动等对土壤有机碳及活性组分的影响。各变量的基底效应系数 C0/(C0+ C)(%)表明，DOC/SOC(%)具有强空间相关性，其空间变异主要受湿地区结构因素影响;其他 4个变量均为中等空间相关性，其空间变异受到来自结构因素和随机因素的耦合影响。

　　2. 2 有机碳及活性组分空间分布格局

　　从图 2(a)、(b)、(c)来看，SOC 和 EOC 含量空间分异格局高度相似，均为由西南 － 南部向东北 －北部呈现高、低、高、低的近似西北 － 东南向的条带状分布，但 EOC 高值区分布范围与 SOC 不完全一致。DOC 含量空间分异与 SOC、EOC 整体上有一致性，但也有明显不同，东南部高，由东南向西北逐渐降低、西北向东北逐渐升高。【图2.略】

　　EOC / SOC( % ) 整体上为近似西南 － 东北向条带并由东南向西北的高、低、高、低变化，并呈现若干同心圆局地高值分布区(图 2(d))。EOC/SOC(%)条带走向(西南 － 东北)和 SOC 条带走向(东南 － 西北)表明两个变量具有一定的负相关特征，但 EOC/SOC(%)与 EOC 的空间分异态势不一致更明显。

　　EOC / SOC(%) 和 SOC、EOC 的相关系数分别为－ 0. 201(通过 0. 05 显着性检验，双尾)、－ 0. 007，也再现了彼此之间空间分异特征的异同。

　　DOC / SOC( % ) 空间分异格局近似为西北 － 东南向条带，但由西南向东北表现为低、高、低、高的变化态势(图 2(e))，整体上与 SOC 高低分异态势相反，而与 DOC 空间分异格局不同;DOC/SOC(%)和SOC、DOC 相关系数分别为 － 0. 499 ( 通过 0. 001 显着性检验，双尾)、－ 0. 023，表明 DOC/SOC(%)和SOC 为极显着负相关，而与 DOC 相关性未通过0. 05 水平显着性检验。

　　3 讨论

　　在土壤有机碳及活性组分研究中，SOC 测试方法很成熟，不同研究获得的 SOC 可进行比较。但对于 LOC 组分而言，测试方法不同意味不同的组分类型。而且，到目前为止，即便是针对特定活性组分，如 EOC、DOC、MBC(Microbial Biomass C)等，还未形成广泛接受的分析方法。以 DOC 为例，有直接采集土壤水测试、采集湿土或风干土水溶提取上清液测试等方法，尚未有研究分析不同测试方法得到的土壤 DOC 含量的差异。因此，在进行不同研究的活性有机碳组分含量比较时，需要考虑测量方法的相对一致性。

　　3. 1 活性有机碳组分含量水平

　　湿地区表土 EOC 含量均值水平为若尔盖湿地区沼泽土易氧化有机碳含量一半，但湿地区的泥炭土、沼泽土和部分典型沼泽化草甸土(15 个点位)的表土 EOC 含量变化于 8. 0 ～30. 0 g/kg，与若尔盖湿地区的沼泽土、泥炭土相当。考虑纳帕海湿地区土壤在整体上为亚高山草甸土，我们选择了两个案例研究进行比较。吴建国等对宁夏六盘山林区牧草地 0 ～10 cm、10 ～20 cm 层土壤采样分析(高锰酸钾氧化比色法)，0 ～ 10 cm、10 ～ 20 cm 层土壤EOC 含量均值分别为 1. 156、1. 959 g / kg。蔡晓布等对藏北典型草地区的正常、轻度退化和严重退化草地 0 ～ 10 cm、10 ～ 20 cm 层土壤采样分析其EOC 含量( 重铬酸钾 － 硫酸氧化法)，三类草地土壤上下两层 EOC 含量均值分别变化于 0. 777 ～ 0. 802g / kg、0. 953 ～ 1. 312 g / kg。相较于海拔更高和纬度偏北的藏北高寒草地、纬度偏北且由半湿润向半干旱区过渡的六盘山林区牧草地而言，纳帕海湿地区的亚高山草甸土表层土壤 EOC 含量均值较高。

　　纳帕海湿地区表土 DOC 含量均值为 155. 360mg / kg，远低于王忠梅等研究得到的三江平原湿地、张金波等研究得到的三江平原小叶章湿地的表土 DOC 含量(均采用新鲜土 － 水溶法)，但纳帕海湿地区的泥炭土、沼泽土、典型沼泽化湿草甸土表层土壤 DOC 含量变化于 307. 14 ～550. 14 mg/kg，接近三江平原岛状林和小叶章湿地表土 DOC 含量均值。从亚高山草甸土整体来看，纳帕海湿地区表土DOC 含量均值明显高于藏西北典型草地区的三类草地土壤 0 ～10 cm、10 ～20 cm 层的 DOC 均值含量(2. 710 ～ 2. 970 mg / kg、3. 520 ～ 4. 520 mg / kg)。

　　3. 2 有机碳活性水平及其指示意义
　　
　　纳帕 海 湿 地 区 表 土 EOC/SOC (%) 均 值 为10. 388。对比有关研究发现，不同土地类型表土SOC、EOC 含量差异虽然较大，但 EOC / SOC( % ) 均值较接近，如若尔盖湿地区沼泽和泥炭表土 EOC/SOC(%)均值为 9. 22%、8. 31%;藏西北典型草地区的三类草地 0 ～ 10 cm、10 ～ 20 cm 土壤 EOC/SOC( % ) 均值分别变化于 10. 82 ～ 11. 13% 、10. 87～ 11. 58%;六盘山 8 类土地类型 0 ～10 cm、10 ～20 cm 层土壤 EOC / SOC( % ) 均值变化于 10. 0% ～12. 5% 、9. 8% ～ 12. 9%。从亚高山草甸土来看，较草地和森林土壤，纳帕海湿地区表土 EOC含量较高，EOC/SOC(%)略低;而相较典型泥炭土和沼泽土则正好相反。

　　纳帕 海 湿 地 区 表 土 DOC/SOC (%) 均 值 为0. 403%(0. 036% ～ 0. 943%)，所有点位均低于1% 。对比相关研究案例发现，利用风干土或湿土(field moist soil)经常温水溶测试得到的若干类型土壤的表土 DOC/SOC(%) 均值基本低于1. 5% 。相对 EOC 而言，土壤 DOC 对水分的分异和变化响应更敏感。将 DOC/SOC(%) 由低而高排序，将 20 个低值和 20 个高值(对应的 DOC)点位挑选出来;同样将 DOC 由高而低排序，将 20 个高值和20 个低值( 对应的 DOC / SOC( % )) 点位挑选出来。两者的对应关系如图 3 所示。【图3】

　　结合野外调查记录，以 DOC/SOC(%)排序筛选(图3(a))为例分析。在 DOC/SOC(%)值低的20个点位中，尽管有 6 个点位 DOC 含量也较低，但这20 个点分别为沼泽或泥炭地(9 个) 、湿草甸(10) 或位于水域边缘(1 个)。DOC/SOC(%)高的 20 个点位中，DOC 含量基本偏低，这些点中有 10 个为中生草甸，或为局地高地，或靠近农田/公路，其余 10 个为湿草甸。图 3(b)也显示出一定的对应关系，特别是 DOC 高值的 7 个点全部为泥炭土、沼泽土或典型的湿草甸土，对应的 DOC/SOC (%) 明显较低。

　　Zhang 等和 Ahn 等研究也指出，湿地土壤表土DOC 含量要较森林、农田、弃耕地以及其他类型土壤表土 DOC 含量高，但 DOC/SOC(%)明显要低。

　　综合本研究和其他案例研究来看，自然发育土壤表土 LOC 组分含量绝对值和相对值(占 SOC 比例)应表现为一定程度上的负相关特征;典型湿地土壤和其他类型土壤(如森林、草地等)的两者对应关系表明水分条件是最重要的决定因子之一;相较LOC 组分含量绝对值而言，LOC / SOC ( % ) ，特别是DOC / SOC(%)，可能对外界环境(特别是水分)的分异及其变化的响应更敏感。但外来干扰可能会改变这一对应关系。以本研究 DOC 为例，141 个点位DOC 和 DOC / SOC( % ) 相关系数仅为 － 0. 023，未通过 0. 05 显着性检验，可能与其他因素对土壤的影响有关。

　　以上分析也意味着:气候变干或人为的疏水排干会导致湿地向中生草甸的演替，从而带来湿地区土壤 LOC(特别是 DOC)组分含量下降，但土壤有机碳活性水平反而增加，会加剧向下游水体的碳输送。

　　湿地土壤 DOC 组分含量及其向下游水体的输送量的变化已经备受关注，因为 DOC 是水域生态系统重要的降解资源，且能吸附污染物、衰减 UV－ B 辐射，促进水生生物生长，进而影响水体水质。

　　高纬或高寒湿地土壤 DOC 含量一般较高，是向河流和海洋输送有机碳极为重要的碳库，气候变化和人类活动驱动下的湿地土壤 DOC 含量及输送通量的变化应予以加强。

　　3. 3 空间变异和分异的影响因素

　　半方差函数分析中，影响土壤属性空间变异的因素被分为结构因素和随机因素。结构因素主要有土壤母质、植被类型、地貌、水文情势等，自然状态下发育的土壤受结构因素为主的控制，土壤属性具有强的空间相关性。随机因素主要来自人为活动、土壤采样的随机性和土样测试误差等，削弱土壤属性的空间相关性。在纳帕海湿地区，微地貌分异整体上控制着湿地区水文和植被生态的分异，是决定湿地区表土 SOC 及活性组分空间变异和分异最基本和最重要的因子。家畜过载散养、道路和水利工程建设、骑马观光旅游等人为活动必然削弱湿地区土壤属性的空间相关性。

　　采样过程中，部分地段因淹水较深而未采样，也会带来一定的影响。由表 2 可知，湿地区表土 SOC、EOC、DOC 均具有中等空间相关性，受到结构因素和随机因素的耦合影响，但其基底效应系数均接近25%，意味着结构因素的影响仍起着主导作用。

　　表土有机碳及活性组分高值区主要分布在保存下来的泥炭地、沼泽地、沼泽化草甸、湿草甸和易洪泛的低洼湿地区。但在东南部靠近农田和公路的区域也是高值区，调查发现，该区原本是低洼且在湿季难以通行的沼泽地或湿草甸，后因发展农牧业、朗举刷河和奶子河裁弯取直(排水)等，逐渐演变成狼毒群落(Com． Euphorbia jolkinii Boiss—Stellera chamae-jasme) 集中成片分布的中生草甸。但从局地地貌来看，该区地势相对较低，在雨季能形成较大面积的滞水或浅积水环境，狼毒群落会显着增加群落斑块下表层土壤有机质，大型家畜极少进入狼毒群落区觅食，这些都利于该区有机碳及活性组分的积累，从而形成表土有机碳及活性组分的高值分布区。

　　湿地区的沟渠堤化、环湿公路建设和近年来日益增多的骑马观光(路线相对固定)等，显着改变各自沿线的水文 － 生态情势，必然对土壤有机碳及其活性组分的空间变异和分异产生显着影响。特别是河渠堤化直接改变了湿地区的洪泛过程和区域，对SOC、EOC 和 DOC 近似东南 － 西北条带走向分布产生明显的影响。

　　湿地区近万头大型家畜过载散养，长期践踏、翻拱和啃食等对湿地区水文、植物群落和土壤产生显着影响，如在湿地区西北部靠近纳帕村一带，大型家畜长期散养导致该区地表植被稀疏，且形成大面积的破坏裸地，形成表土有机碳及其活性组分的低值区。而且，大型家畜的散养对湿地区土壤的影响在空间上是无序的，这必然增大湿地区土壤有机碳及其活性组分的空间变异和分异，削弱湿地区土壤属性变异的空间相关性，也会带来土壤碳素的流失。

　　4 结 论

　　1. 纳帕海湿地区土壤为典型的亚高山草甸土，湿地区表土 SOC、EOC、DOC、EOC/SOC(%)、DOC/SOC(%) 均值分别为 51. 596 g/kg、4. 950 g/kg、155. 360 g / kg、10. 388% 、0. 403% ;各变量的空间变异水平分别为强、强、中、中、中。

　　2. 纳帕海湿地区表土有机碳和两类活性组分之间为极显着正相关;EOC/SOC(%)和 SOC 为显着负相关，与 EOC 为负相关但不显着;DOC/SOC(%)和 SOC 为极显着负相关，与 DOC 为负相关但不显着。各变量的空间分异格局的一致性和差异性也验证了彼此之间的相关显着性。

　　3. 纳帕海湿地区表土 DOC / SOC( % ) 空间变异具有强的空间相关性，受结构因素为主的控制;其他均为中等空间相关性，受结构因素和随机因素的耦合影响。SOC、EOC、EOC/SOC(%)空间相关性的各向异性特征显着，而 DOC、DOC/SOC(%)空间相关性的各向异性特征不显着。

　　4. 在纳帕海湿地区，泥炭和沼泽表土有机碳及活性组分含量均值要明显高于农田和中生草甸土，但其有机碳活性水平均值低于农田和中生草甸土。作为亚高山草甸土，纳帕海湿地区表土有机碳及活性组分含量均值高于毗邻青藏高原区典型高寒草地表土、低于典型泥炭沼泽土表土，而活性水平正好相反。这意味着不同类型土壤的表土活性组分含量绝对值均值和相对值均值(占 SOC 比例)为一定程度负相关，水分可能是影响这一特征的重要决定因子。

　　5. 相较于 EOC、DOC 含量绝对值而言，EOC /SOC( % ) 和 DOC / SOC( % ) 能更好地反映环境以及人为干扰的空间分异对表土有机碳活性组分含量及其活性水平空间分异的影响。