全 文 :澜沧江中游山地不同土地利用对土壤
有机碳的影响及预测∗
刘世梁1∗∗ 安南南1 杨珏婕1 董玉红2 王 聪1
( 1北京师范大学环境学院 /环境模拟与污染控制国家重点联合实验室,北京 100875; 2中国林业科学研究院林业研究所,北京
100091)
摘 要 基于 2011年实地采集的 210个表层(0~20 cm)土壤样本,研究了澜沧江中游山地 4
种典型土地利用类型土壤有机碳含量的空间分布差异.利用相关性分析研究土壤有机碳含量
与土地利用类型、地形因子、归一化植被指数之间的关系,采用多元线性回归模型和克里格插
值模型预测澜沧江中游山地有机碳含量分布.结果表明: 研究区 4 种土地利用类型表层土壤
有机碳含量的大小依次为:林地>灌丛>草地>农田.土地利用类型、坡向、坡度和平面曲率均是
影响土壤有机碳含量的重要因素.多元回归模型比地统计学模型对山地生态系统土壤有机碳
的预测效果更好.
关键词 土地利用; 土壤有机碳; 澜沧江; 多元线性回归; 克里格插值
文章编号 1001-9332(2015)04-0981-08 中图分类号 S153.6;U412.1 文献标识码 A
Effects of different land⁃use types on soil organic carbon and its prediction in the mountain⁃
ous areas in the middle reaches of Lancang River. LIU Shi⁃liang1, AN Nan⁃nan1, YANG Jue⁃
jie1, DONG Yu⁃hong2, WANG Cong1 ( 1 State Key Laboratory of Water Environment Simulation,
School of Environment, Beijing Normal University, Beijing 100875, China; 2Research Institute of
Forestry, Chinese Academy of Forestry, Beijing 100091, China) . ⁃Chin. J. Appl. Ecol., 2015, 26
(4): 981-988.
Abstract: Based on the analysis of 210 surface soil samples (0-20 cm) according to field investi⁃
gation in 2011, soil organic carbon (SOC) content and its distribution under four different land⁃use
types in the middle reaches of Lancang River were studied. Combining with the correlation analysis,
relationships among SOC content, land⁃use type, topographical factors and NDVI were revealed.
Furthermore, multiple linear regression and Kriging method were used to predict the SOC contents
distribution in the mountainous areas in the middle reaches of Lancang River. The results exhibited
that the rank order of SOC content under the four land⁃use types was forestland > scrubland > grass⁃
land > farmland. The important factors influencing the SOC content were land⁃use type, aspect,
slope and curvature. Comparing the accuracy of these two models, we found that the multiple regres⁃
sion model performed better than geo⁃statistics in prediction of the SOC content in the mountainous
region.
Key words: land⁃use type; soil organic carbon; Lancang River; multiple linear regression; Kriging
method.
∗环保公益项目(2012090290)资助.
∗∗通讯作者. E⁃mail: shiliangliu@ bnu.edu.cn
2014⁃05⁃26收稿,2014⁃12⁃24接受.
土壤有机碳(soil organic carbon,SOC)及其空间
分布特征是反映土壤质量或土壤健康的一个重要指
标,直接影响土壤肥力和生产力[1] .由于土壤有机碳
库的库容巨大[2],其微小变化就将影响大气 CO2浓
度,因而在地球碳循环过程中起着重要作用[3] .土壤
有机碳受到自然因子[4]及人为干扰等诸多因素的
影响,其中,土地利用变化是影响土壤有机碳动态的
关键因素[5] .土地利用方式的改变影响土壤有机质
的输入和输出方式,进而引起土壤有机碳含量的变
化[6] .流域内水电开发对土地利用 /土地覆盖的影响
是最为敏感的人类社会活动之一[7] .漫湾水电站位
应 用 生 态 学 报 2015年 4月 第 26卷 第 4期
Chinese Journal of Applied Ecology, Apr. 2015, 26(4): 981-988
于澜沧江流域的中游,该水电站于 1986 年开工,截
至 1996年底,漫湾电站共淹没耕地 414.97 hm2(其
中,水田 241.92 hm2、旱地 173.04 hm2)、林地 567.17
hm2(含经济林 69.57 hm2)、荒山荒地 765.54 hm2 [8] .
水库的蓄水,导致一系列现象出现,如水位升高、水
库水面面积增大、原有河谷景观消失[9]、大量移民
的产生、对土地利用进行改造、废弃原有耕地、开发
新的生产或生活用地等,淹没区的产生及其土地利
用方式的改变都导致土壤有机质含量发生变化.
Zhao等[10]以漫湾水电站为例分析了库区土地利用
的变化及各土地类型间的转化;杨少红[11]、李正
才[12]和陈亮中[1]分析了不同土地利用方式对土壤
有机碳的影响,认为不同土地利用类型的土壤碳含
量不同,并分析了土壤碳分布的影响因子.目前,对
于澜沧江中游山区这一特殊生境,研究不同植被类
型对土壤有机碳的影响,探讨以水电工程建设为主
的人为干扰对土壤有机碳含量影响的相关研究成果
较少.
澜沧江中游山区位于峡谷型山区,土壤有机层
薄、有机质含量低、土壤抗蚀性差而易发生水土流
失,山区生态系统具有极大的潜在脆弱性,岸坡陡
峻,自然环境较为复杂,再加上漫湾水电站的建设导
致人类活动频繁,改变了原有土地利用方式及土地
管理措施,使得该区域碳循环表现出极大的敏感性
和典型性,急需探索有效方法对此复杂的土壤有机
碳的空间分布格局进行精确预测.多元回归分析和
地统计学方法已经被广泛应用于定量分析土壤有机
碳的空间分布特征.如 Cheng 等[13]建立土壤有机碳
含量与土壤母质、地形因子之间多元回归方程,预测
了江西省兴国县土壤有机碳含量的空间分布.利用
多元线性回归和数字地形模型进行各种土壤属性的
空间分布预测在土壤学领域已经得到了广泛应
用[14-15] .基于地统计学的克里格方法是当前预测区
域土壤有机碳的主要方法[16-17] .如 McGrath 等[18]和
Zhang等[19]采用地统计学方法对爱尔兰草原土壤有
机碳含量的空间分布和时空变化进行预测.赵永存
等[20]结合地形因子,采用多元线性回归、泛克里格
和回归克里格方法,预测河北省土壤有机碳密度的
空间分布特征,认为回归克里格法的预测效果最好,
泛克里格法次之,而多元线性回归方法最差.张忠启
等[21]比较分析了普通克里格法和结合土壤类型信
息的克里格法 2 种方法对土壤有机碳的空间预测
精度.
水电站的建设会改变土地利用方式,进而影响
土壤碳含量及其分布.在水电站建设和土地利用变
化过程中,土壤有机碳库将会如何变化,不同的土地
利用类型对土壤有机碳有何影响,目前相关的研究
较少.因此,本文以澜沧江中游漫湾水电站库区为研
究对象,研究土地利用转化对土壤碳含量、分布及其
未来趋势的影响;利用多元线性回归和克里格插值
2 种方法对漫湾库区土壤有机碳进行空间预测,并
对比分析 2 种方法的预测精度,探讨不同预测模型
在漫湾库区土壤有机碳空间预测精度的可行性,以
期为该地区土壤碳储存现状与潜力的估算和评估奠
定基础.
1 研究地区与研究方法
1 1 研究区概况
澜沧江是一条国际河流,全长约 4500 km,在中
国境内长 2153 km,总落差约 5060 m,共规划有 4库
22级水电站,总装机容量约 31.89 GW[22] .漫湾水电
站(24°25′—24°40′ N、100°05′—100°25′ E)地处澜
沧江中游云南省境内,位于云县与景东县交界处的
漫湾镇上游附近.漫湾水电站是澜沧江流域梯级水
电开发建设的第一级电站,坝长 418 m、高 132 m,正
常蓄水位 994 m,总库容 10.6×108 m3,水库面积 23.6
km2,干流回水约 70 km,总装机容量为 150 × 104
kW[6] .本研究所指的库区范围包括电站水库所涉及
的澜沧江河段两岸分水岭以内的区域,上游至小湾
电站附近,下游至坝址所在山脊[23](图 1).研究区属
横断山系南部帚状山脉峡谷中山区,为河道型水
库[24] .库区河谷属于印度洋季风气候,年均气温
18~20 ℃,年降水量 1000 ~ 1150 mm[25],年降水量
的季节分配不均,有明显的干季和雨季,雨季多为 7
月或 8月,干季出现在 12 月或 1 月[26] .两岸植被覆
盖率高、类型多样,主要有河岸半落叶阔叶混交林、
山地针叶林、陡坡高草稀树林、山地常绿阔叶林、河
滩灌丛和荒地灌草林[23] .研究区域内,主要的土壤
群和土壤类型分别是红壤和沙土,其中,红壤的 pH
在 4.5~ 5.2、粘重、保水保肥力差、耕性较差,沙土土
壤土质疏松、透水透气性好、有机质含量在 2%左
右、pH在 6.2~7.3、耕种层超过 30 cm[27-28] .
1 2 样地的选择
2011年 6 月,在漫湾库区云县、凤庆、南涧和景
东境内共选取 210 个表层(0 ~ 20 cm)土壤样本(图
1).在漫湾库区采集不同土地利用方式(阔叶林、针
阔混交林、针叶林、人工林、灌丛、草地和农田)下面
积较广、生长良好而集中的典型植被下的土壤样本
289 应 用 生 态 学 报 26卷
图 1 研究区域及样线示意图
Fig.1 Sketch map of the study region and line transect dia⁃
gram.
作为 2006年遥感影像对应的土壤属性数据.漫湾库
区植被类型多样,阔叶林主要以旱冬瓜(Alnus ne⁃
palensis)、麻栎(Quercus acutissima)等为优势种,并
常与马尾松 (Pinus massoniana)、柏树 (Platycladus
orientalis)等混交成林;针叶林主要为云南松(Pinus
yunnanensis);人工林主要为思茅松(Pinus khasys)、
铁核桃(Juglans sigillata)和蓝桉(Blue eucalyptus);
灌木以虾子花 (Woodfordia fruticosa)、革叶算盘子
(Glochidion daltoni)和展毛野牡丹(Melastoma nor⁃
male)为主;草本以金发草 ( Pogonatherum panice⁃
um)、硬杆子草 ( Capillipedium assimile)和飞机草
(Chromolaena odoratum)为主.
1 3 样品采集和处理
在库区范围内,每种植被类型随机选取 5 个样
方,每个样方内随机采取 6 个重复的土壤表层样本
(0~20 cm),共 210个土壤样品.土层原状土壤用环
刀采集,重复 3次,烘干称量,获取库区表层土壤容
重.土样采回混合均匀并风干后,过 100 目筛,然后
用 H2SO4⁃K2Cr2O7外加热法测定土壤有机碳含量.用
1 ∶ 50000 数字高程模型 ( DEM) 图 (分辨率为
50 m ´50 m)提取地形因子.利用 2011年 TM影像计
算得出归一化植被指数(NDVI).
1 4 土壤有机碳含量预测模型
1 4 1多元统计分析方法 基于经典统计学中的多
元线性回归分析(multiple linear regression, MLR),
采用最小二次方的方法进行拟合,用一组自变量去
预测因变量,公式如下[29]:
y = a+∑
n
i = 1
bixi+c (1)
式中:a为截距;bi为回归系数;c为残差.
将漫湾水电站库区 210 个表层土壤碳数据,按
实际采样植被类型进行划分,对应的土壤有机碳平
均含量作为该植被类型的属性编码值.使用 ArcGIS
的空间分析模块在 DEM的基础上,提取出库区的海
拔(elevation)、坡度( slope)、坡向( aspect)、平面曲
率(curvature)和平面粗糙度( roughness),并根据提
取的地形因子计算复合地形指数(CTI).以植被类
型、地形因子、NDVI 为自变量,以土壤有机碳含量
为因变量,利用 SPSS 软件统计分析建立多元线性
回归模型.并按照多元线性回归的最优方程,将数字
类型图进行叠加运算,生成流域土壤有机碳含量空
间分布图.
1 4 2地统计学分析方法 地统计学模型是定量预
测土壤特性空间分异的常用方法之一.本研究用地
统计学的克里格插值法建立库区土壤有机碳预测模
型,并与前述的多元统计学模型进行比较,选取最适
合漫湾小流域的土壤碳预测模型.克里格插值法是
利用半方差图的模型进行观测点之间的最优内插,
是目前在地统计学中应用最广泛的一种线性无偏最
优估计方法,可以直观显示土壤属性的空间分布特
点,公式如下:
Z(x0) =∑
n
i = 1
λ iZ(xi) (2)
式中:Z(x0)为未观测点 x0上的内插估计值;Z( xi)
为在点 x0附近的若干观测点上获得的实测值;λ i为
考虑了半方差图中表示空间的权重.基于克里金插
值法,采用 GIS软件的地统计分析模块,绘制土壤有
机碳的预测图,并进行模型拟合精度的相关分析.
1 4 3预测模型的准确性 空间预测方法的优劣可
采用验证数据集进行验证,通过比较验证点位置上
土壤有机碳含量的预测值与实测值来评价预测的准
确性,采用平均标准差 ( mean standard deviation,
MSD) 和 均 方 根 误 差 ( root⁃mean square error,
RMSE) [30]进行评价,RMSE反映样点数据的估值灵
敏度和极值效应,其值越小,表明预测结果越准确.
1 5 数据处理
本研究用单因素方差分析(ANOVA)分析土地
利用类型对土壤有机碳含量和容重的影响,用 Excel
和 SPSS统计分析软件对结果进行数据统计分析.
3894期 刘世梁等: 澜沧江中游山地不同土地利用对土壤有机碳的影响及预测
2 结果与分析
2 1 不同土地利用类型的土壤有机碳含量
由表 1 可以看出,研究区不同土地利用类型的
土壤容重有所差异,其中,灌丛、草地和农田的土壤
容重差异不显著,而针叶林、针阔混交林和人工林的
土壤容重显著高于灌丛、草地和农田.在漫湾库区范
围的 7种土地利用类型表层土壤中,人工林的土壤
容重最高,为 1.41 g·cm-3,灌丛的土壤容重最低,
为 1.14 g·cm-3 .
不同土地利用类型的土壤有机碳含量也不同,
其中,灌丛、草地和农田土壤有机碳含量的差异不显
著,但天然林地(包括阔叶林、针叶林、针阔混交林)
的土壤有机碳含量显著高于其他 4 种植被类型(表
1).4种土地利用类型的土壤有机碳含量高低顺序
依次为林地>灌丛>草地>农田,农田土壤有机碳含
量明显低于其他 3 种土地利用类型.在 4 种不同林
地类型中,库区表层土壤有机碳含量最高的是阔叶
林 ,为42 .57 g·kg-1 ;最低的是人工林,为16 .44
表 1 不同土地利用类型的土壤容重及土壤有机碳含量
Table 1 Soil organic carbon (SOC) contents and soil bulk
densities for land⁃use types
土地利用类型
Land use
type
土壤容重
Soil bulk
density
(g·cm-3)
平均土壤
有机碳含量
Mean SOC
content
(g C· kg-1)
林地
Forestland
阔叶林
Broad⁃leaved forest
1.18±0.02a 42.57±19.70b
针阔混交林
Mixed coniferous and
broad⁃leaved forest
1.36±0.02b 40.10±19.63b
针叶林
Coniferous forest
1.27±0.03c 37.81±6.15b
人工林
Plantation
1.41±0.02d 16.44±3.24a
灌丛 Shrubland 1.14±0.06a 23.91±2.45a
草地 Grassland 1.17±0.03a 21.54±7.54a
农田 Farmland 1.14±0.03a 18.50±4.47a
F值 F value 74.21∗ 13.19∗
不同小写字母表示不同类型之间差异显著(P<0.05) Different small
letters meant significant difference among different types at 0. 05 level.
∗ P<0.05.
g·kg-1,说明自然林的固碳能力高于人工林.
2 2 土壤有机碳预测模型
云南省地形复杂,海拔高差悬殊较大,坡度大于
25°以上的坡耕地面积 74.87×104hm2,占耕地总面积
的 11.7%.水土流失严重,大面积滑坡和泥石流多
发,坡度 25°以上的坡耕地流失水土达 5542 × 104
t·km-2·a-1 .植被覆盖的减少会加剧水土流失,不
仅会导致流域内水土资源的重新分配,还会引起土
壤有机碳大量流失.因此,地形、NDVI 和植被类型都
是影响表层土壤有机碳空间分布的重要因素.本研
究分别采用多元统计学和地统计学方法建立漫湾库
区土壤有机碳含量的预测模型,分析地形因子、ND⁃
VI和植被类型对表层土壤有机碳空间分布的影响,
并对模型的拟合精度进行比较,选取最优模型.
在实际采样的 210 个剖面点中,随机抽取 158
个(75%)作为插值数据集进行土壤有机碳的空间
预测,其余 52个(25%)作为验证数据集,进行预测
的准确性检验.
2 2 1基于多元统计回归的土壤有机碳模型 由表
2 可知,研究区土壤有机碳与土地利用类型呈极显
著相关,与坡向、坡度和平面曲率显著相关,但与
NDVI和复合地形指数的相关性并不显著.说明在小
流域内,气候和土壤母质差异不大,土地利用类型决
定了土壤有机碳的输入,而坡度、坡向和平面曲率直
接决定了水分运动过程以及坡面物质的迁移和再分
配,因此,土地利用类型和地形因子是影响漫湾库区
土壤发生发育特别是土壤有机碳含量最普遍、最有
效的定量预测因子.
运用向后回归模型进一步筛选和建立研究区土
壤有机碳预测模型(表 3),结果表明,只用土地利用
类型预测土壤有机碳含量模型的 R2为 0.523,加入
复合地形指数 CTI、NDVI和坡向(aspect)后,模型的
参数相对少而且拟合程度最优,R2为 0.563,所以筛
选出如下回归方程:
SOC= 0 71CTI+0 019Aspect+6 225NDVI+
0 793SOC i-3.72 (3)
表 2 土壤有机碳(SOC)与地形因子、NDVI和植被类型的相关系数
Table 2 Correlation coefficient between SOC and terrain factor, NDVI and vegetation types
项目
Item
土地利用类型
Land use type
复合地形
因子 CTI
坡向
Aspect
平面曲率
Curvature
海拔
Elevation
粗糙度
Roughness
坡度
Slope
归一化植被指数
NDVI
Pearson 相关系数
Pearson correlation coefficient
0.73* * 0.21 0.32* -0.36* -0.27 0.06 -0.36* 0.27
P值 P value 0.00 0.19 0.05 0.03 0.09 0.72 0.03 0.09
∗ P<0.05; ∗∗ P<0.01.
489 应 用 生 态 学 报 26卷
表 3 基于向后回归模型的土壤有机碳预测模型
Table 3 Prediction models of SOC based on the backward
regression
模型参数
Model parameter
R2 标准估计的误差
Error of
standard estimate
常量、土地利用类型、复合地形指数、归一化植被
指数、海拔、平面曲率、坡度、粗糙度、坡向 Con⁃
stant, land⁃use types, CTI, NDVI, elevation, cur⁃
vature, slope, roughness, aspect
0.520 4.02
常量、土地利用类型、复合地形指数、归一化植被
指数、海拔、平面曲率、坡度、坡向 Constant, land⁃
use types, CTI, NDVI, elevation, curvature, slope,
aspect
0.536 3.95
常量、土地利用类型、复合地形指数、归一化植被
指数、海拔、平面曲率、坡向 Constant, land⁃use
types, CTI, NDVI, elevation, curvature, aspect
0.549 3.90
常量、土地利用类型、复合地形指数、归一化植被
指数、海拔、坡向 Constant, land⁃use types, CTI,
NDVI, elevation, aspect
0.562 3.85
常量、土地利用类型、复合地形指数、归一化植被
指数、坡向 Constant, land⁃use types, CTI, NDVI,
aspect
0.563 3.85
常量、土地利用类型、复合地形指数、归一化植被
指数 Constant, land⁃use types, CTI, NDVI
0.560 3.85
常量、土地利用类型
Constant, land⁃use types
0.523 4.01
表 4 研究区土壤有机碳模型预测精度的比较
Table 4 Comparison of accuracy for SOC prediction mod⁃
els in the study area (g·kg-1)
项目
Item
观测值
Observation
模型方法 Model method
多元回归法
Multiple
regression
method
克里格
插值法
Kriging
method
均值 Mean 3100 1800 2600
最大值 Maximum 13000 14500 13000
最小值 Minimum 450 330 450
标准方差 Standard deviation 1500 780 980
均方根误差 Root⁃mean⁃square eror 0.625 0.752
平均标准差 Mean standard deviation 0.511 0.675
均方根标准误差 Root⁃mean⁃square
standard error
1.069 1.418
2 2 2土壤有机碳模型的比较 采用平均标准差
(MSD)、均方根误差(RESD)作为精度检验的标准,
检验不同方法的预测效果.均方根误差和平均标准
差越接近于 0,均方根标准误差越接近于 1,说明模
型的精确度和稳定性越好.从表 4 可以看出,多元回
归模型比地统计学模型在漫湾库区土壤有机碳的预
测效果更好,适用性更强.
3 讨 论
3 1 不同土地利用类型对土壤有机碳含量的影响
土壤碳固定是缓解温室效应的有效方法之一,
逐渐成为研究热点.土壤有机碳含量主要取决于有
机物料的投入量、分解释放量及土壤性质和过程,而
土壤碳的主要输入源是地上植物的凋落物.土地利
用变化对土壤有机碳含量的影响主要是通过改变凋
落物的数量和质量及其环境条件的改变而影响有机
碳的储存量、组成和稳定性[31] .凋落物是影响土壤
有机碳空间分布的重要因子[32] .漫湾库区各土地利
用类型 0~20 cm土层的有机碳含量大小依次为:林
地>灌丛>草地>农田.相比其他土地利用类型,农田
土壤碳含量在 1979 和 2011 年都最低,主要原因在
于:1)收获农作物时地上部分被大量移走,大大减
少了植被对土壤的归还量;2)机械耕作使土壤温度
和湿度得到改善,为好氧微生物的活动创造了条件,
在一定程度上促进了土壤呼吸作用,加速了土壤有
机质的分解;3) 耕作等人类活动导致土壤结构发生
改变,降低了土壤的物理保护作用,造成土壤侵
蚀[33-34] .林地土壤碳含量远高于其他土地利用类
型,原因在于其表层覆盖有大量的枯枝落叶,加上该
地区气候湿润,导致腐殖质层较厚,引起土壤碳的积
累.人工林的土壤碳含量低于其他几种林地类型,主
要原因归结于人工林受人类活动的干扰较大,说明
人类活动对土壤碳储量影响的程度和速率远远超过
自然因素[35] .
3 2 土壤碳含量对土地利用变化的响应
土地利用发生变化后,影响土壤有机碳的动态
过程和相关因子发生变化,导致土壤有机碳的储量、
分布等也发生变化,进一步又影响土壤向大气释放
CO2 的强度与过程,因此土地利用变化对土壤碳储
量、分布的影响是温室效应的关键内容[12,30] .研究表
明,森林变成农田或草地后土壤有机碳在一定程度
上将会发生变化.李正才[12]认为,次生林转变成农
耕地后土壤有机碳平均含量下降了 28. 2%. Lug
等[36]研究结果表明,亚热带地区的森林被砍伐变成
农田后 10年,土壤有机碳含量下降 46%,80 年后下
降 70%.林地转为农地和草地一般会使土壤有机碳
降低的原因:一方面植被净生产力降低(或残留物
移除)减少了土壤有机碳输入量,另一方面是覆被
类型改变使地表土壤温度升高因而加快了土壤有机
碳分解,原来稳定的有机碳变得不稳定并加大了渗
透量[34] .而灌丛、农田和草地变成人工林后土壤有
机碳会大幅升高,说明林地能增加土壤有机碳的积
累, 提高土壤有机碳的有效性.该地区位于库区范
围内,林地砍伐变成农田或其他用地通常是毁灭性
的,增加当地土壤有机碳积累的不确定性.天然林地
砍伐变成人工林后,一般会引起表层土壤有机碳的
5894期 刘世梁等: 澜沧江中游山地不同土地利用对土壤有机碳的影响及预测
损失:人工林地的管理措施与天然林不同,人为活动
的加剧会导致土壤容重大于其他天然林地类型,水
土流失严重,这是导致人工林地表层土壤有机碳低
于其他林地类型的主要原因之一.有研究表明,人工
林幼林郁闭前,林地上通过凋落物和枯死细根的有
机碳归还较少[37],同时,人工林地的凋落物和细根
的分解速率也低于天然林[38],这些因素都会导致人
工林地的碳含量低于天然林地.本研究区域位于印
度洋季风气候区,湿润多雨,年降水量为 1000~1150
mm,加之本区山高多坡陡,自然生态环境具有潜在
的脆弱性,当自然林遭砍伐更替为人工林时,常造成
严重的土壤有机质流失.
3 3 土壤有机碳模型的比较
多元线性回归模型将自变量与因变量之间的关
系假设为是线性的,通过单方差因素分析表明,土壤
有机碳与植被类型呈极显著相关,与坡向、坡度和平
面曲率呈显著相关,通过向后回归模型,筛选出坡
向、CTI 和 NDVI这 3个自变量建立土壤有机碳的多
元线性回归模型.
通过对比两种模型的预测值和实际观测值,可
以评价模型预测制图的精度.本研究中多元线性回
归模型比克里格插值模型对土壤有机碳的预测效果
好,适用性更强.海拔、坡度、坡向可以综合体现环境
变量,在本研究中海拔、坡度和坡向都影响土壤有机
碳含量,一方面不同环境变量梯度植被类型及植被
的垂直地带性、水平地带性的不同进而影响土壤特
性[39-40],从而影响土壤有机物质的质和量;另一方
面又通过温度和水分等生态因子的变化影响土壤有
机质的分解和腐殖化过程,土壤地表枯落物腐殖化
过程和有机碳释放过程受诸多因素影响.研究区的
土壤有机碳含量由植被类型和地形特征所决定,而
地统计学模型是根据土壤样本变异与样本距离之间
的关系建立,忽略了土壤母质、植被类型和局部微气
候对土壤有机碳的影响.在漫湾库区,植被类型多
样、地形复杂、水热分配不均,土壤有机碳的空间分
布主要与植被类型和地形因子以及植被覆盖有关,
因此基于这几个因素建立的多元统计模型精确度和
稳定性更好.
4 结 论
利用 2011年的土壤剖面有机碳与土地利用图
的斑块连接和叠加,绘制土壤有机碳的空间分布差
值图,以比较土壤有机碳的空间变化.然后,分别用
多元统计学和地统计学的方法建立漫湾库区土壤有
机碳含量的预测模型,分析地形因子、NDVI 和植被
类型对表层土壤有机碳空间分布的影响,并对模型
的拟合精度进行比较,选取最优模型.结果表明:澜
沧江流域中游不同土地利用类型的土壤有机碳含量
不同,0~20 cm深度土壤有机碳含量的大小依次为:
林地>灌丛>草地>农田.天然林地(阔叶林、针阔混
交林和针叶林)转化为农田会导致土壤有机碳含量
分别降低 56.5%、53.9%和 51.1%,说明自然林地能
增加土壤有机碳的积累,提高土壤有机碳的有效性.
而退耕还林过程中,如果将农田转变为人工林并不
会增加土壤有机碳含量,反而会降低 12.5%.库区影
响土壤有机碳含量的显著因素主要是植被类型、坡
向、坡度和平面曲率.由此建立多元回归模型并与克
里金插值法模型进行比较,结果表明,多元回归模型
比地统计学模型在漫湾库区的预测效果更好,适用
性更强.
参考文献
[1] Chen L⁃Z (陈亮中). Studies on Soil Organic Carbon
under Main Forest Vegetation Types in the Three⁃Gorges
Reservoir Area. PhD Thesis. Beijing: Beijing Forestry
University, 2007 (in Chinese)
[2] Davidson EA, Trumbore SE, Amundson R. Soil war⁃
ming and organic carbon content. Nature, 2000, 408:
789-790
[3] Gao YN, Fang HJ, Yu GR, et al. Spatial distribution
and temporal dynamics of soil carbon removal caused by
water erosion in China. Journal of Resources and Ecolo⁃
gy, 2011, 2: 210-216
[4] Guo Y⁃F (郭月峰), Yao Y⁃F (姚云峰), Qin F⁃C (秦
富仓), et al. Impact of terrain factors on soil organic
carbon in Laoha River basin. Journal of Arid Land Re⁃
sources and Environment(干旱区资源与环境), 2014,
28(2): 156-161 (in Chinese)
[5] Zhang Y⁃G (张于光), Zhang X⁃Q (张小全), Xiao Y
(肖 烨). Effects of land use change on soil organic
carbon and microbial biomass carbon in Miyaluo forest
area. Chinese Journal of Applied Ecology(应用生态学
报), 2006, 17(11): 2029-2033 (in Chinese)
[6] Ding Y⁃K (丁越岿), Yang J (杨 劼), Song B⁃Y
(宋炳煜), et al. Effect of different vegetation types on
soil organic carbon in Mu Us desert. Acta Partaculturae
Sinica (草业学报), 2012, 21(2): 18- 25 ( in Chi⁃
nese)
[7] Zhou Q (周 庆), Ou X⁃K (欧晓昆), Zhang Z⁃M
(张志明), et al. Spatial temporal land use pattern
changes in Manwan hydropower station reservoir of Lan⁃
cang River, Yunnan, China. Journal of Mountain Sci⁃
ence (山地学报), 2008, 26(4): 481-489 ( in Chi⁃
nese)
[8] Zhou J⁃H (周嘉慧), Fu B⁃H (付保红). Cultivated
land use of the migrant areas of large⁃sized hydropower
689 应 用 生 态 学 报 26卷
stations in Yunnan Province: A case study of Manwan
Dam. Tropical Geography (热带地理), 2008, 28(6):
551-554 (in Chinese)
[9] Wei G⁃L (魏国良), Cui B⁃S (崔保山), Dong S⁃K
(董世魁), et al. Impact of hydropower development on
river ecosystem service: A case study from the Manwan
Hydropower Project. Acta Scientiae Circumstantiae(环境
科学学报), 2008, 28(2): 235-242(in Chinese)
[10] Zhao QH, Liu SL, Deng L, et al. Determining the influ⁃
encing distance of dam construction and reservoir im⁃
poundment on land use: A case study of Manwan Dam,
Lancang River. Ecological Engineering, 2013, 53:
235-242
[11] Yang S⁃H(杨少红). The Effect of Land Use Change on
Soil Organic Carbon and Soil Respiration. PhD Thesis.
Fuzhou: Fujian Agriculture & Forestry University, 2006
(in Chinese)
[12] Li Z⁃C(李正才). The Effects of Land⁃use Change on
the Soil Organic Carbon. PhD Thesis. Beijing: China
Academy of Forestry, 2006 (in Chinese)
[13] Cheng XF, Shi XZ, Yu DS, et al. Using GIS spatial
distribution to predict soil organic carbon in subtropical
China. Pedosphere, 2004, 14: 425-431
[14] Florinsky IV, Eilers RG, Manning GR, et al. Prediction
of soil properties by digital terrain modelling. Environ⁃
mental Modelling & Software, 2002, 17: 295-311
[15] Li QQ, Yue TX, Wang CQ, et al. Spatially distributed
modeling of soil organic matter across China: An appli⁃
cation of artificial neural network approach. Catena,
2013, 104: 210-218
[16] Zhao TC, Xu XH, Huang B, et al. Using robust kriging
and sequential Gaussian simulation to delineate the cop⁃
per⁃ and lead⁃contaminated areas of a rapidly industrial⁃
ized city in Yangtze River Delta, China. Environmental
Geology, 2007, 52: 1423-1433
[17] Kerry R, Oliver MA. Determining nugget: Sill ratios of
standardized variograms from aerial photographs to krige
sparse soil data. Precision Agriculture, 2008, 9: 33-56
[18] McGrath D, Zhang CS. Spatial distribution of soil organ⁃
ic carbon concentrations in grassland of Ireland. Applied
Geochemistry, 2003, 18: 1629-1639
[19] Zhang CS, McGrath D. Geostatistical and GIS analyses
on soil organic carbon concentrations in grassland of
southeastern Ireland from two different periods. Geoder⁃
ma, 2004, 119: 261-275
[20] Zhao Y⁃C (赵永存), Shi X⁃Z (史学正),Yu D⁃S (于
东升), et al. Different methods for prediction of spatial
patterns of soil organic carbon density in Heibei Prov⁃
ince, China. Acta Pedologica Sinica (土壤学报 ),
2005, 42(3): 379-385 (in Chinese)
[21] Zhang Z⁃Q (张忠启), Yu F⁃Z (于法展). Study on the
application of soil type information in spatial prediction
of soil organic carbon. Chinese Agricultural Science Bul⁃
letin(中国农学通报), 2013, 29(11): 139-144 ( in
Chinese)
[22] Liu Q (刘 琦),Liu S⁃L (刘世梁),Zhao Q⁃H (赵清
贺), et al. Gradient analysis of landscape pattern affect⁃
ed by hydroelectric station based on moving window
method. Journal of Mountain Science (山地学报),
2012, 30(5): 628-635 (in Chinese)
[23] Cui B⁃S (崔保山), Zhai H⁃J (翟红娟). Quality evalu⁃
ation of habitats disturbed by the Manwan Hydropower
Dam. Acta Scientiae Circumstantiae (环境科学学报),
2008, 28(2): 227-234 (in Chinese)
[24] Fu B⁃H (付保红), Chen L⁃H (陈丽晖), Zhu T (朱
彤). An analysis of ecological environment change in
the storage areas of Manwan hydropower station and its
managing countermeasures. Territory & Natural Re⁃
sources Study (国土与自然资源研究), 2005(1): 54-
55 (in Chinese)
[25] Wang P (王 平), Wang J⁃L (王金亮). Basic climat⁃
ic characteristics in the Manwan hydropower station res⁃
ervoir area of Lancang River, Yunnan Province, China.
Journal of Yunnan Normal University(Natural Science)
(云南师范大学学报:自然科学版), 2002, 22(6):
57-61 (in Chinese)
[26] Zhao Q⁃H (赵清贺), Liu S⁃L (刘世梁), Zhang Z⁃L
(张兆苓), et al. Effects of Manwan hydropower exploi⁃
tation on dynamic changes of landscape pattern in relat⁃
ed reservoir area. Chinese Journal of Ecology (生态学
杂志), 2011, 30(10): 2343-2350 (in Chinese)
[27] Zhao QH, Liu SL, Deng L, et al. Soil degradation asso⁃
ciated with water⁃level fluctuations in the Manwan Res⁃
ervoir, Lancang River Basin. Catena, 2014, 113:
226-235
[28] Zhou L⁃F (周乐福). Characteristics and zonality of soil
distribution in Yunnan Province. Mountain Research (山
地调查), 1983, 1(4): 31-38 (in Chinese)
[29] Florinsky IV, Eliers RG, Manning GR, et al. Prediction
of soil properties by digital terrain modeling. Environ⁃
mental Modelling & Software, 2002, 17: 295- 311
[30] Liu S (刘 莎), Ren H⁃Y (任红艳), Shi X⁃Z (史学
正), et al. Different methods for prediction of spatial
patterns of paddy soil organic carbon density in Chang⁃
xing County, Zhejiang Province. Journal of Geoinforma⁃
tion Science (地球信息科学学报),2010, 12 ( 2):
2180-2185 (in Chinese)
[31] Jiang G⁃Y (姜贵勇). The implementation of returning
farmland to forest project policy and technology in Jing⁃
dong County. Forest Inventory and Planning (林业调查
规划), 2006, 31(suppl.1): 214-217 (in Chinese)
[32] Wang G⁃Q (王国强), Mao Y⁃L (毛艳玲). The effect
of land use change on semitropics soil organic carbon.
China High Technology Enterprises (中国高新技术企
业), 2008(19): 123-126 (in Chinese)
[33] Qiao Y⁃M (乔有明), Wang Z⁃Q (王振群), Duan Z⁃H
(段中华). Effects of different land use types on soil
carbon and nitrogen contents in the northern region of
Qinghai Lake. Acta Prataculturae Sinica (草业学报),
2009, 18(6): 105-112 (in Chinese)
[34] Jiang Y (姜 勇), Zhuang Q⁃L (庄秋丽), Liang W⁃J
(梁文举). Soil organic carbon pool and its affecting
factors in farm land ecosystem. Chinese Journal of Ecolo⁃
gy (生态学杂志), 2007, 26(2): 278-285 ( in Chi⁃
7894期 刘世梁等: 澜沧江中游山地不同土地利用对土壤有机碳的影响及预测
nese)
[35] Li M⁃F (李明峰), Dong Y⁃S (董云社), Qi Y⁃C (齐
玉春), et al. Effect of land use change on the contents
of C & N in temperate grassland soils. Grassland of Chi⁃
na (中国草地), 2005, 27(1): 2-7 (in Chinese)
[36] Lugo A, Sanchez MJ. Land use and organic carbon con⁃
tent of some subtropical soils. Plant and Soil, 1986,
96: 185-196
[37] Chen Z (陈 朝), Lyu C⁃H (吕昌河), Fan L (范
兰), et al. Effects of land use change on soil organic
carbon: A review. Acta Ecologica Sinica (生态学报),
2011, 31(18): 5358-5371 (in Chinese)
[38] Yang Y⁃S (杨玉盛), Xie J⁃S (谢锦升), Sheng H (盛
浩), et al. The impact of land use / cover change on soil
organic carbon stocks and quality in mid⁃ subtropical
mountainous area of Southern China. Acta Geographica
Sinca (地理学报), 2007, 62(11): 1123- 1131 ( in
Chinese)
[39] Li X⁃H (李秀花), Shi Q⁃D (师庆东), Guo J (郭
娟), et al. The response of NDVI to climate variability
in northwest arid area of China from 1981 to 2001. Jour⁃
nal of Arid Land Resources and Environment (干旱区资
源与环境), 2009, 23(2):12-16 (in Chinese)
[40] Zhang G (张 国), Cao Z⁃P (曹志平), Hu C⁃J (胡
婵娟). Soil organic carbon fractionation methods and
their application in farmland ecosystem research: A re⁃
view. Chinese Journal of Applied Ecology (应用生态学
报), 2011, 22(7): 1921-1930 (in Chinese)
作者简介 刘世梁,男,1976年生,博士生导师,副教授.主要
从事景观生态学、GIS 与 RS 研究. E⁃mail: shiliangliu@ bnu.
edu.cn
责任编辑 杨 弘
封 面 说 明
图片由中国科学院生态环境研究中心城市与区域生态国家重点实验室于洋 2014 年 7 月拍摄
于甘肃省定西市龙滩流域,该流域属半干旱黄土丘陵沟壑区,海拔 1840 ~ 2260 m,流域面积
16.1 km2,多年平均降水量 386.3 mm,降水量季节分配不均,春季降水较为稀少,主要集中在 7—9
月,且多暴雨事件发生,水土流失较为严重.流域年均气温为 6.8 ℃,平均无霜期 152 d,潜在蒸发量
为 1649.0 mm,年均相对湿度 72%.流域土壤类型为黄绵土,流域内主要土地利用类型包括撂荒地、
灌木林地、乔木林地、农地、人工草地.同时分布着零星的天然草地.流域内的人工植被以柠条、侧
柏、山杏、油松、紫花苜蓿为主,天然植被以多年生草本为主,主要种类包括赖草、长芒草、阿尔泰狗
娃花等.多年来,中科院生态环境研究中心景观生态研究团队联合地方相关部门,在流域内开展了
长期系统的研究.伴随退耕还林工程的实施,流域内形成了不同的植被恢复模式,主要包括:退耕地
人工林草植被可持续经营模式、荒坡水土保持乔灌草配置模式、梁峁顶植被生态功能修复模式、侵
蚀沟道治理模式和农田安全与地埂保育模式.在多重植被恢复模式的基础上,流域内不仅构建了有
效的水土流失综合治理体系,也形成了丘陵沟壑交错的美丽景观.
889 应 用 生 态 学 报 26卷