福建省土壤有机碳储量估算的尺度效应研究*

陈中星 张 楠 张黎明† 袁 平 姚彩燕 邢世和邱龙霞 陈翰阅 范协裕

（1 福建农林大学资源与环境学院，福州 350002）

（2 土壤生态系统健康与调控福建省高校重点实验室，福州 350002）

土壤有机碳在生态系统中起着重要作用，其储量约为1 500 Pg（1 Pg=1015g），是大气碳库的3倍、植被碳库的2.5倍，其较小幅度变化将会对大气中二氧化碳（CO2）浓度产生剧烈影响［1］。我国的土地面积居世界第三位，土壤类型众多，在全球温室气候变化研究中占据着十分重要的地位［2］。近年来，国际社会不断寻找控制气候变化的最佳途径，土壤固碳效应被广泛认为是能够延缓全球气候变化的有效措施［1，3］，而精确估算土壤有机碳储量和明确其空间分布特征是制定固碳减排政策的前提。

目前，基于土壤普查资料并结合数字化土壤图进行有机碳储量估算的方法越来越普遍。国外，Eswaran［4］于1993年基于16 000 个土壤剖面和其他已经发表的文献数据，采用FAO-UNESCO土壤图估算出全球土壤有机碳的储量大约为1 576 Pg。Martín等［5］基于4 401个样点对西班牙有机碳储量进行估算，结果表明该地区表层土壤（0～20 cm）有机碳储量和平均密度分别为2.8 Pg和56.57 Mg hm-2。Veronesi等［6］使用英国和威尔士的土壤调查数据对米德兰兹郡西部面积约为13 948 km2区域的碳储量进行了估算，结果表明，该地区表层和深层土壤有机碳储量分别为162.91 ± 3.63 × 10-5Tg和66.71 ± 3.52 × 10-5Tg。Dorji等［7］基于数字化土壤图和186个土壤剖面对喜马拉雅东部的山地森林系统的有机碳储量进行估算，结果表明该区域有机碳储量为27.1 Mt。国内利用Shi等［8］提出的“PKB”（Pedological Knowledge Based method）法进行土壤有机碳储量估算的研究比较多。该方法主要根据土壤类型相一致与相似性、土壤成土母质相近或者相同、土壤剖面点位置与分布区域的一致或邻近等原则，将土壤数据库中每个土壤剖面的属性数据与空间数据库中相对应的土壤类型图斑单元连接起来。如，于东升等［9］利用1∶100万土壤数据库，并采用“PKB连接法”将7 292个剖面数据和土壤图中的94 000多个图斑连接，通过面积加权平均法得出我国的土壤有机碳储量为89.14 Pg，有机碳平均密度为9.60 kg m-2。Zhi 等［10］利用浙江省1∶5万土壤数据库和2 154个土壤剖面属性资料，采用“PKB连接法”估算出全省的土壤有机碳储量为831.49 Tg，土壤有机碳平均密度为8.69 kg m-2。李晓迪等［11］采用“PKB连接法”估算出苏北旱地地区在1∶5万土壤数据库基础上有机碳密度和储量分别为1.56 kg m-2和61.20 Tg （1 Tg=1012g）。

但是从以上的研究可以看出，我国的国家和区域土壤有机碳储量估算中大多数使用的是单一制图尺度土壤数据库，未考虑大、中、小系列比例尺土壤数据库对有机碳估算的影响。不同比例尺土壤数据库对有机碳空间分布的表达能力差异较大，土壤图图斑及其连接的属性会随着比例尺的改变而发生明显变化［12］。有统计表明，目前我国土壤有机碳储量估算的最小值为50 Pg，而最大值为185.69 Pg，二者相差近4倍，这说明使用不恰当的制图尺度将会导致估算结果产生极大的不确定性［13-14］；此外，国家尺度的土壤碳储量估算大多数是以小比例尺（1∶100万、1∶400万、1∶1 00 0万）为基础资料，而这些小比例尺下的图斑概化对有机碳储量估算精度的影响程度尚不清楚［15-16］。福建省森林覆盖率高达65.95%，居全国第一，森林生产力高，土壤输入有机质量多，有利于土壤有机碳积累［17-18］，明确该地区有机碳储量估算精度受目前常用制图比例尺土壤数据库的影响程度，对我国国家尺度和亚热带地区制定合理的固碳减排政策提供有力依据。基于此，本研究以整个地区均属于亚热带的福建省作为研究区，分别计算“PKB连接法”建立的目前我国常用6种制图尺度（1∶5万、1∶20万、1∶50万、1∶100万、1∶400万、1∶1 000万）表层和剖面土壤数据库下的有机碳密度和储量，并以目前土壤数据最为详尽的1∶5万尺度估算值为基准，定量评价其他制图尺度的误差，结果可为我国亚热带地区土壤有机碳研究中选择适宜的制图尺度和定量化不同制图尺度的估算精度提供理论依据。

1 材料与方法

1.1 研究区概况

福建省位于我国东南沿海（115°50′～120°40′E，23°33′～28°20′N），属于亚热带海洋性季风气候（图1）。根据全省66个气象站点的统计得出，年平均气温为14.6～21.3 ℃，年均降水量为1 037～2 051 mm，气候温暖，雨量充足。福建省土壤类型以红壤为主，达7.66 × 106hm2，占全省土壤总面积63.41%；水稻土为分布最广的耕作土壤，面积仅次于红壤，总面积1.07 ×106hm2，占全省土壤总面积8.71%；其他土壤类型（如滨海盐土、潮土、风砂土、山地草甸土、石灰土、新积土等）面积较小［19］。

1.2 数据来源

图1 福建省地理位置分布图Fig. 1 Geographical location map of Fujian Province

本研究所使用的制图尺度土壤数据库来自全国第二次土壤普查成果，由属性数据库和空间数据库两部分组成。其中，数据最为详细的1∶5万尺度土壤图空间数据库来自福建省8 4个县（市）第二次土壤普查纸质资料。借助ArcGIS工作平台，采用双标准纬线等积圆锥投影对第二次土壤普查纸质图件进行数字化、编辑、修改并建立拓扑关系，建立该地区的空间数据库，全省共计1 2个土类、2 3个亚类和8 3个土属。福建省第二次土壤普查的各县（市、区）《土种志》中共记录的剖面点为3 082个，其属性数据包含土壤类型名称、土壤肥力和土壤理化性质等指标（图2）。由于第二次土壤普查时期没有GPS（Global Positioning System）仪器进行剖面点的准确定位，本研究根据第二次土壤普查《土种志》中记录的剖面点详细采样位置，以村界为限定区域，找出与该剖面土壤类型一致的土壤图最大面积“图斑”，“图斑”的中心点即认为是第二次土壤普查剖面采集位置，从而实现了纸质记录剖面点位置的GIS（地理信息系统）空间化。空间数据与属性数据的连接采用Shi等［8］提出的“PKB”进行连接。在连接过程中，既要考虑剖面的土壤类型与图斑土壤类型I D值相同，也要兼顾剖面的空间位置，当出现多个剖面点与一个图斑相对应的情况时，取多个剖面点的属性平均值赋予该图斑，而当一个剖面点对应多个图斑时，则取该剖面点的属性值赋予相对应的多个图斑［16］。所有土壤图制图单元均采用中国土壤发生分类系统（Genetic Soil Classification of China，GSCC），随着制图比例尺的降低，制图单元所对应的土壤分类级别也随之升高，1∶5万～1∶1 000万这6个制图比例尺的制图单元分别为土种、土属、土属、土属、亚类和亚类。在制图尺度降低过程中，由于图斑概化导致的土壤类型归并，6个尺度土壤图的图斑数目依次为247 969、46 408、15 282、6 343、440和345个，分别利用了3 082、3 082、2 718、2 547、1 000和890个土壤剖面的理化属性数据（表1）。

图2 1982年福建省土壤剖面点分布图Fig. 2 Distribution of soil profiles （1982） in Fujian Province

表1 福建省1∶5万～1∶1 000万6个制图比例尺土壤图基本特征Table 1 Basic features of the six soil maps of Fujian different in mapping scale from 1∶50 000 ～ 1∶10 000 000

1.3 土壤有机碳密度及储量计算

本研究采用土壤类型法估算福建省土壤有机碳储量，以表层土壤0～20 cm和剖面土壤0～100 cm作为计算对象。

土壤有机碳密度指单位面积一定深度土层中的土壤有机碳储量。土壤有机碳密度的计算公式为［20］:

式中，SOCD为土壤剖面有机碳密度（kg m-2）；n为参与计算的土壤层次总数；äi%为第i层＞2 mm砾石含量（体积百分含量）；ρi为第i层土壤容重（g cm-3）；Ci为第i层土壤有机碳含量（g kg-1），由土壤有机质含量×0.58（Bemmelen转换系数）换算得到［21］；Ti为第i层土层厚度（cm）；100用于单位换算。由于1982年福建省土壤普查剖面点资料中缺少容重信息，本研究采用Song等［22］基于第二次土壤普查建立土壤容重与有机碳含量的回归模型：y=1.377×e-0.0048×SOC（R2=0.7870，p＜0.001，n=4 765），式中SOC代表土壤有机碳含量。

土壤有机碳储量指一定区域范围内土层中所包含的土壤有机碳总量。土壤有机碳储量的计算公式为［23］:

式中，SOCS为土壤剖面有机碳储量（kg）；SOCDi为第i个土壤图斑的有机碳密度（kg m-2）；Si为第i个土壤图斑的面积（m2）；n为图斑数。

1.4 不同制图比例尺下的相对偏差

本次研究采用相对偏差（%）的绝对值来表示不同制图尺度之间的土壤有机碳密度、储量和土壤面积差异大小：

式中，Y为相对偏差（%）；xs为1∶20万、1∶50万、1∶100万、1∶400万和1∶1 000万尺度计算的有机碳密度、储量和土壤面积；x0为土壤数据最详细的1∶5万尺度计算的有机碳密度、储量和土壤面积。

1.5 不同制图比例尺的差异性检验

应用统计软件SPSS 22.0对不同制图尺度下的土壤有机碳密度独立样本T检验［24］，分析1∶20万、1∶50万、1∶100万、1∶400万和1∶1 000万尺度的有机碳密度与土壤数据最为详尽的l∶5万尺度之间的差异显著性。

2 结果与讨论

2.1 不同制图尺度下福建省土壤有机碳密度与储量分布

通过面积加权平均方法得到福建省1∶5万～1∶1 000万6个制图尺度下表层（0～20 cm）土壤有机碳含量分别为：18.69 g kg-1、19.42 g kg-1、17.64 g kg-1、18.39 g kg-1、18.88 g kg-1和19.11 g kg-1；剖面（0～100 cm）土壤有机碳含量分别为：8.92 g kg-1、9.10 g kg-1、8.54 g kg-1、8.94 g kg-1、9.23 g kg-1和9.27 g kg-1。以数据最为详细的1∶5万制图尺度为基准，1∶50万尺度下表层和剖面土壤有机碳含量相对偏差最大，分别为5.59%和4.28%；1∶400万制图尺度下表层相对偏差最小为1.03%，1∶100万制图尺度下剖面相对偏差最小为0.20%。结果表明：福建省表层土壤有机碳含量远大于剖面土壤有机碳含量，这是由于表层地表枯落物较为丰富，有机质含量较高，而剖面土壤受枯落物的影响较低，有机质含量相对较小［25］；此外，受到人为活动的影响，表层活性有机碳会随土壤深度的增加而降低，从而导致有机碳含量降低［26］。

从土壤有机碳密度来看，在数据最为详细的1∶5万制图尺度下福建省土壤图斑约为247 969个，估算出全省土壤面积为12.08×106hm2，表层（0～20 cm）和剖面（0～100 cm）土壤有机碳密度分别为4.57 kg m-2和11.55 kg m-2，有机碳储量分别为552.23 Tg和1 396 Tg。从表2可以看出，福建省表层土壤有机碳密度高于全国平均水平［14］和临近的广东［27］、广西［28］两省；此外，剖面土壤有机碳密度也高于全国平均水平和临近的广东［27］、浙江［10］两省，主要是由于该省山脉较多、气候温和湿润、植物覆盖面积大等有利于土壤有机碳积累的特殊地理环境所致［18］。

表2 不同省份表层和剖面土壤有机碳密度与储量Table 2 SOC densities and SOCSs in the surface soil layers and profiles of different provinces

从图3a中可以看出，福建省表层土壤有机碳密度空间分布差异很大，最小值为0.09 kg m-2，最大值为13.62 kg m-2，二者之间相差接近158倍。土壤有机碳密度较大（＞10 kg m-2）地区主要分布于福建省东北部、北部以及西部地区，合计面积为0.41×106hm2，占研究区土壤总面积的3.36%。福建省中部和西部地区相对于东部的气温较低，降雨量较大，低温湿润的环境条件能促进植被生长，使得有机碳积累的更快。土壤有机碳密度较小（＜2.5 kg m-2）地区主要分布于福建省东南沿海地区，合计面积为2.19× 106hm2，占研究区土壤总面积的18.12%，这主要与该地区温度高，降雨量少不利于有机碳的积累有关。从图4a中可以看出，福建省剖面土壤有机碳密度空间分布差异也很大，最小值为0.10 kg m-2，最大值为50.42 kg m-2，二者之间相差超过500倍。剖面土壤中有机碳密度较高（＞20 kg m-2）的土壤分布于福建省东北部、北部以及西部地区，与土壤表层分布规律一致，面积为1.28×106hm2，占研究区土壤总面积的10.62%；土壤有机碳密度较小（＜5 kg m-2）的区域同样主要分布于福建省东南沿海地区。

从图3b～图3f和图4b～图4f可以看出，随着制图尺度的减小，福建省整个表层和剖面的土壤有机碳密度空间分布变化差异明显。以1∶5万制图尺度为基准，表层和剖面土壤有机碳密度较大（＞10 kg m-2和 ＞20 kg m-2）区域面积发生明显变化，1∶20万和1∶50万尺度下该含量范围的面积急剧增加，相对偏差分别为85%～124%和11%～28%；而1∶100万～1∶1 000万尺度下明显减少，相对偏差分别为38%～44%和7%～49%。表层和剖面土壤有机碳密度较小（＜2 kg m-2和＜5 kg m-2）区域面积在不同制图尺度下也发生明显变化，相对偏差分别为1%～58%和7%～92%；其中1∶400万和1∶1 000万这两个制图尺度下有机碳储量和密度估算结果有较高的不确定性，这主要是小比例尺下土壤类型的概化过程较其他尺度更为剧烈所致。1∶400万制图尺度下滨海盐土、潮土、粗骨土、山地草甸土、新积土和紫色土这6个土类因图斑面积太小而被概化；1∶1 000万制图尺度下除了上述6个土类被概化外，石灰土和风砂土也同样被概化，这也说明当制图尺度发生改变时，一方面不同土壤有机碳密度含量的空间分布位置会发生变化，如随着制图尺度减小，土壤有机碳密度较高的区域从福建省东北部和北部逐渐转移至福建省西北部和中部地区，而土壤有机碳密度较小区域位置并未发生改变，仍为福建省东南沿海区域；另一方面制图尺度较大的土壤图中分布面积大、范围较为广泛的土壤类型可能是制图尺度较小的土壤图制图单元中重要组成部分，分布面积小，范围较窄的土壤类型则可能在图斑概化过程中被删除或者归化至其他土类中［24］。

图3 福建省不同制图尺度土壤表层（0～20 cm）有机碳密度空间分布Fig. 3 Spatial distribution of organic carbon density in the surface soil layer（0～20 cm）of Fujian Province relative to mapping scale of the soil map

以1∶5万尺度表层估算值为基准值，1∶100万尺度下的有机碳密度和储量相对偏差最小，分别为0.82% 和0.32%，而1∶20万尺度下的相对偏差最大，分别为15.57% 和15.34%，其他1∶50万、1∶400万和1∶1 000万尺度下的相对偏差分别为5.05% 和3.80%、12.64% 和11.13%、4.05%和3.79%。以1∶5万尺度剖面估算值为基准，同样1∶100万尺度下的有机碳密度和储量相对偏差最小，分别为0.26% 和0.10%，而1∶400万和1∶1 000万尺度下的相对偏差最大，分别为9.46% 和8.00%、9.98% 和9.72%，其他1∶20万、1∶50万和1∶400万尺度下的相对偏差分别为7.84% 和7.58%、4.22% 和5.35%。表层下1∶20万尺度下的相对偏差最大是因为该尺度下土壤有机碳密度较高（＞7.5 kg m-2）的区域面积为2.74 M hm2，占1∶20万研究区总面积的22.72%，相比于1∶5万土壤有机碳密度较高的区域面积（1.51 M hm2，占1∶5万研究区总面积的12.48%）增加近一倍，而在剖面下1∶400万和1∶1 000万尺度的相对偏差最大是由于图斑的概化，分别仅有1 000个和890个土壤剖面进入属性库，有机碳密度较小但面积较大的风砂土、滨海盐土、潮土等被归化为密度较大的水稻土和红壤，从而使得这两个尺度的土壤有机碳密度明显高于其他尺度。

上述研究表明：福建省土壤有机碳密度空间分布中大比例尺有机碳密度空间分布相对分散，具有一定的科学分布规律，而小比例尺有机碳密度在空间分布上趋于集中，这是由于小比例尺土壤数据库受概化过程的影响，图斑的单位面积增大、土壤剖面数量减少、土壤类型的相互归并所导致。随着制图尺度发生变化，某些土壤的空间分布发生改变，导致有机碳储量发生改变。有些密度较小的土壤如滨海盐土、粗骨土等在随着制图尺度发生变化时，这些土壤类型被概化为有机碳密度较大的红壤、黄壤等。同样某些密度较大的土壤如山地草甸土和新积土在制图尺度变小时，归并到水稻土、黄壤等有机碳相对密度较小的土壤类型中。由此可知，不同制图尺度所带来的估算结果差异性是由于空间分布所导致。

2.2 不同制图尺度下各土壤类型有机碳密度与储量分布

图4 福建省不同制图尺度土壤剖面（0～100 cm）有机碳密度空间分布Fig. 4 Spatial distribution of organic carbon density in the soil profiles （0～100 cm） of Fujian Province relative to mapping scale of the soil map

从表3可以得出，在数据最详细的1∶5 万尺度下，福建省表层和剖面土壤中有机碳密度最大的是山地草甸土，分别达到7.75 kg m-2和29.07 kg m-2，该土类主要分布于中山山顶的平缓部位，年平均气温低（18.8 ℃），海拔较高（1 430 m），降水量较大（1 667 mm），草甸植被生长茂密且分解缓慢，使得土壤有机质富集明显［29］；其次，黄壤的有机碳密度也较高，表层和剖面分别为6.48 kg m-2和16.18 kg m-2，该土类主要分布于中山地带，该区域温度低湿度高（18.2 ℃、1 709 mm），导致表层积累大量的有机质。水稻土受到人为施肥的影响，且在水耕熟化过程中有机碳普遍呈现出积累的趋势［30］，导致表层和剖面有机碳密度分别达到3.93 kg m-2和11.17 kg m-2。福建省有机碳密度最低的是风砂土和粗骨土，表层和剖面的有机碳密度分别为0.94 kg m-2、2.17 kg m-2和1.22 kg m-2、3.22 kg m-2，这主要是因为前者主要分布于东南沿海地区，该地区干旱少雨气候不利于有机碳积累；此外，该土类中砂粒含量较高，占颗粒组成的85%以上，有研究表明，有机碳密度与土壤砂粒的含量呈极显著负相关关系［31］；而后者主要分布于边缘山丘地区，有机质和土壤养分贫瘠导致有机碳密度低。其他土类表层和剖面土壤有机碳密度分别介于2.19～4.83 kg m-2和5.18～11.91 kg m-2之间。

随着制图尺度的变化，土壤有机碳密度受制图尺度影响最大的土类为粗骨土，在1∶5万尺度下表层和剖面的有机碳密度分别为1.22 kg m-2和3.22 kg m-2，而在1∶50 万尺度下高达4.55 kg m-2和11.01 kg m-2，相对偏差为273.24%和241.72%。土壤有机碳储量受制图尺度影响最大的也是粗骨土，在1∶5 万尺度下分别为0.001 Tg和0.0 035 Tg，而在1∶50 万尺度下分别为11.86 Tg和11.01 Tg，相对偏差分别达8.88×105%和8.13×105%。首先，该土类在1∶5万数据库下只有16个图斑，面积为1.10×102hm2，在1∶20万和1∶50万数据库下分别有107和378个图斑，面积分别为6.10×103hm2和2.61×105hm2，其次，在1∶20万和1∶50万数据库下密度相对较大的山地草甸土、黄壤、红壤、水稻土等由于制图尺度的降低，被概化或归并为粗骨土，导致该土类有机碳密度增大。其次，不同制图尺度下石灰土表层和剖面土壤有机碳密度差异也很大，1∶20万和1∶50万尺度下相对偏差均大于60%，而在1∶100万和1∶400万制图尺度下相对偏差均小于35%，这主要是由于在制图尺度降低过程中，该土类的图斑数量和面积均呈下降趋势。红壤和水稻土分别是福建省分布最广的自然和农业土壤，这两种土类储量占各制图尺度总储量的80%以上，这一方面由于二者的土壤有机碳密度相对较高，不同尺度下表层和剖面分别在4.83～5.76 kg m-2、3.92～4.89 kg m-2和11.12～13.28 kg m-2、11.09～14.08 kg m-2之间；另一方面由于他们所占的面积较大，合计在不同尺度土壤总面积的80%以上。但相对而言，红壤受不同制图尺度的影响相对较小，表层和剖面土壤相对偏差均在4%～19%之间；而水稻土受不同制图尺度的影响相对较大，表层土壤有机碳密度的相对偏差在4%～19%之间，剖面土壤有机碳密度的相对偏差在11%～54%之间。

2.3 不同制图尺度下各行政区域有机碳密度与储量分布

表4可以看出，在1∶5万制图尺度下，南平市和龙岩市表层和剖面土壤有机碳密度最高，分别为5.41 kg m-2和15.05 kg m-2、5.44 kg m-2和12.28 kg m-2，这一方面由于两市背靠武夷山脉，距海较远，年平均气温低且降雨量较大（18.0 ℃、1 722 mm和19.5 ℃、1 649 mm），低温多雨的气候条件有利于有机碳积累；另一方面由于两市耕地面积占福建省耕地资源比重最大，分别为16.69%、13.75%［32］，长期且持续肥料的施用以及生物积累能够有效增加土壤中有机碳含量。厦门市表层和剖面土壤有机碳密度最低，分别为3.03 kg m-2和7.84 kg m-2，由于距海较近，年平均气温高且降雨量较小（20.9 ℃、1 401 mm），且滨海地区土壤砂粒含量高，导致养分贫乏，从而有机碳含量低。

随着制图比例尺变化，表层和剖面土壤有机碳密度受制图尺度影响最大的是福州市，以1∶5万尺度估算值为基准值，在1∶400万制图尺度下的相对偏差最大，分别达到33.99%和40.14%；此外，在1∶5万制图尺度下该市有16 922个图斑，包含新积土、石灰土、风砂土、潮土等11种土壤类型，而1∶400万制图尺度下只有54个图斑，仅有水稻土、红壤、黄壤、赤红壤4种土壤类型。表层土壤有机碳密度和储量受制图尺度影响最小的是泉州市，1∶20万～1∶1 000万5个尺度下的相对偏差分别在1.24%～18.32%和0.06%～17.12%

之间；而剖面土壤有机碳密度受制图尺度影响最小的是龙岩市，1∶20万～1∶1 000万5个尺度下的相对偏差分别在0.48%～8.69%和0.43%～7.84%之间。总体而言，1∶50万制图尺度下仅有厦门市剖面土壤有机碳密度相对偏差大于10%，其他行政区域剖面土壤有机碳密度和储量相对偏差均小于10%。很多学者研究表明，由于陆地系统是一个自然等级系统，在该系统内不同尺度上的主导过程是不同的，地表信息在时间上和空间上的分辨率均有很大的跨度，在某一尺度上人们观察到的性质、总结出来的原理和规律，在另一个尺度上可能仍然有效，可能有相似性，也可能需要修改［33］。从本研究的结果来看，在福建省不同地级市下1∶50万制图尺度空间有机碳密度及储量与目前大区域尺度数据最详尽的1∶5万数据库比较接近。因此，福建省各地级市在缺乏精确土壤数据库的情况下，可使用1∶50万制图尺度进行有机碳储量估算和农业管理措施的制定，这一方面可保证相对比较高的精度，另一方面数据容易获取，也可省时省力。

表3 不同制图尺度数据库下福建省各土壤类型有机碳储量（Tg）和密度（kg m-2）Table 3 SOC density（kg m-2） andSOC storage （Tg） in the soils of Fujian Province relative to soil type and mapping scale of the soil map

表4 不同制图尺度下福建省各行政区域土壤有机碳密度（kg m-2）和储量（Tg）Table 4 SOC density （kg m-2） and SOC storage （Tg） in different administrative areas of Fujian Province relative to mapping scale of the soil map

3 结 论

明确不同制图尺度土壤数据库是正确估算亚热带土壤有机碳储量，实施“固碳减排”国家战略的基础。从整个福建省来看，1∶100万尺度下估算的表层（0～20 cm）和剖面（0～100 cm）土壤有机碳密度及储量与目前大区域数据最详细的1∶5万尺度相对偏差比较小，分别为0.67% 和0.82%；而1∶20万尺度相对偏差比较大，分别达到15.57% 和15.34%。目前，1∶100万尺度是国家有机碳估算中使用最为广泛的比例尺，本研究一定程度上也证明了该尺度在国家级别土壤有机碳研究应用中的可靠性。从不同土壤类型来看，粗骨土受到制图比例尺影响最大，1∶50万尺度下表层和剖面土壤有机碳密度及储量的相对偏差分别达到273.24%、8.88×105%和241.72%、8.13×105%；面积最大的红壤有机碳密度相对偏差在表层土壤下1∶100万尺度最小，为4.28%，而在剖面土壤下1∶50万最低，仅为1.53%。从不同行政区来看，除厦门市外，其他地级市表层和剖面土壤有机碳密度的相对偏差均在1∶50万尺度下最低，这说明各地级市在缺乏精确土壤数据库的情况下，可使用该尺度进行有机碳储量估算和农业管理措施的制定。

［ 1 ］ Lal R. Soil carbon sequestration impacts on global climate change and food security. Science，2004，304（7）：1623—1627

［ 2 ］ 潘根兴，赵其国. 我国农田土壤碳库演变研究：全球变化和国家粮食安全. 地球科学进展，2005，20（4）：384—393 Pan G X，Zhao Q G. Study on evolution of organic carbon stock in agricultural soils of China：Facing the challenge of global change and food security （In Chinese）. Advances in Earth Science，2005，20（4）：384—393

［ 3 ］ Pan G X，Lu H F，Li L Q. et al. Soil carbon sequestration with bioactivity：A new emerging frontier for sustainable soil management. Advance in Earth Sciences，2015，30（8）：940—951

［ 4 ］ Eswaran H. Organic carbon in soils of the world.Soil Science Society of America Journal，1993，57：192—194

［ 5 ］ Martín J A R，Álvaro-Fuentes J，Gonzalo J，et al.Assessment of the soil organic carbon stock in Spain.Geoderma，2016，264：117—125

［ 6 ］ Veronesi F，Corstanje R，Mayr T. Landscape scale estimation of soil carbon stock using 3D modelling.Science of the Total Environment，2014，487：578—586

［ 7 ］ Dorji T，Odeh I O A，Field D J，et al. Digital soil mapping of soil organic carbon stocks under different land use and land cover types in montane ecosystems，Eastern Himalayas. Forest Ecology and Management，2014，318：91—102

［ 8 ］ Shi X Z，Yu D S，Warner E D，et al. Cross-Reference system for translating between genetic soil classification of China and soil taxonomy. Soil Science Society of America Journal，2006，70：13—22

［ 9 ］ 于东升，史学正，孙维侠，等. 基于1︰100万土壤数据库的中国土壤有机碳密度及储量研究. 应用生态学报，2005，16（12）：2279—2283 Yu D S，Shi X Z，Sun W X，et al. Estimation of China soil organic carbon storage and density based on 1︰1 000 000 soil database （In Chinese）. Chinese Journal of Applied Ecology，2005，16（12）：2279—2283

［10］ Zhi J J，Jing C W，Lin S P，et al. Estimates of soil organic carbon stocks in Zhejiang Province of China based on 1︰50 000 soil database using the PKB method. Pedosphere，2015，25（1）：12—24

［11］ 李晓迪，王淑民，张黎明，等. 土壤数据源和制图比例尺对旱地土壤有机碳储量估算的影响. 土壤学报，2016，53（1）：58—71 Li X D，Wang S M，Zhang L M，et al. Impacts of source of soil data and scale of mapping on assessment of organic carbon storage in upland soil （In Chinese）. Acta Pedologica Sinica，2016，53（1）：58—71

［12］ 张勇. 滇黔桂地区土壤有机碳空间表征的尺度效应研究. 南京：中国科学院南京土壤研究所，2008 Zhang Y. Effect sacle of soil organic carbon spatial characterization in China area （In Chinese）.Nanjing：Institute of Soil Science，Chinese Academy of Sciences，2008

［13］ 梁二，蔡典雄，张丁辰，等. 中国陆地土壤有机碳储量估算及其不确定性分析. 中国土壤与肥料，2010（6）：75—79 Liang E，Cai D X，Zhang D C，et al. Terrestrial soil organic carbon storage in China：Estimates and uncertainty （In Chinese）. Soil and Fertilizer Sciences in China，2010（6）：75—79

［14］ 解宪丽，孙波，周慧珍，等. 中国土壤有机碳密度和储量的估算与空间分布分析. 土壤学报，2004，41（1）：35—43 Xie X L，Sun B，Zhou H Z，et al. Organic carbon denstiy and storage in soils of China and spatial analysis （In Chinese）. Acta Pedologica Sinica，2004，41（1）：35—43

［15］ 李加加. 苏北旱地土壤有机碳估算的尺度效应研究. 福州：福建农林大学，2013：14—22 Li J J. Scale effect of soil organic carbon estimation in the uplands of north Jiangsu Provinve，China （In Chinese）. Fuzhou：Fujian Agriculture and Forestry University，2013：14—22

［16］ 黄琳斌. 旱地有机碳情景分析模拟的尺度效应及其影响因素研究. 福州：福建农林大学，2016 Huang L B. Study on scale effect and its influencing factors in the scenario of simulating soil organic carbon dynamics in the upland （In Chinese）. Fuzhou：Fujian Agriculture and Forestry University，2016

［17］ 王义祥，翁伯琦. 福建省土壤有机碳密度和储量的估算. 福建农业学报，2005，20（1）：42—45 Wang Y X，Weng B Q. Estimate of soil organic carbon density and its stock in Fujian Province （In Chinese）. Fujian Journal of Agricultural Sciences，2005，20（1）：42—45

［18］ 福建省环境保护厅. 2015年福建省环境状况公报.2016-05. http：//www. fjepb. gov. cn/zwgk/kjjc/hjzl/qshjzkgb/ 201606/ P020160603649157645861. pdf Fujian Provincial Department of Environmental Protection. Bulletin of the environmental situation of Fujian Province in 2015 （In Chinese）. 2016-05.http：//www. fjepb. gov. cn/zwgk/kjjc/hjzl/qshjzkgb/201606/ P020160603649157645861. pdf

［19］ 福建省土壤普查办公室. 福建土壤. 福州：福建科学技术出版社，1991：287—289 Soil Survey Office of Fujian Province.Fujian soil （In Chinese）. Fuzhou：Fujian Science & Technology Publishing House，1991：287—289

［20］ Pan G X，Li L Q，Zhang Q，et al. Organic carbon stock in topsoil of Jiangsu Province，China，and the recent trend of carbon sequestration. Journal of Environmental Sciences，2005，17（1）：1—7

［21］ 袁芳，赵小敏，乐丽红，等. 江西省表层土壤有机碳库储量估算与空间分布特征. 生态环境，2008，17（1）：268—272 Yuan F，Zhao X M，Le L H，et al. Organic carbon reserves of tops oil and spatial distribution in Jiangxi province （In Chinese）. Ecology and Environment，2008，17（1）：268—272

［22］ Song G H，Pan G X，Zhang Q. Topsoil SOC storage of China agricultural soils and its loss by cultivation.Biogeochemistry，2005，74（1）：47—62

［23］ Wang S Q，Tian H Q，Liu J Y，et al. Pattern and change of soil organic carbon storage in China：1960s—1980s. Tellus Series B：Chemical and Physical Meteorology，2003，55（2）：416—427

［24］ 李昕，张明明. SPSS 22. 0 统计分析从入门到精通. 北京：电子工业出版社，2015 Li X，Zhang M M. SPSS 22. 0 Statistical analysis from entry to mastery （In Chinese）. Beijing：Publishing House of Electronics Industry，2015

［25］ 渠开跃，冯慧敏，代力民，等. 辽东山区不同林型土壤有机碳剖面分布特征及碳储量研究. 土壤通报，2009，40（6）：1316—1320 Qu K Y，Feng H M，Dai L M，et al. Profile distribution and storage of soil organic carbon of main forest types in eastern mountainous region of Liaoning（In Chinese）. Chinese journal of soil science，2009，40（6）：1316—1320.

［26］ 王彩霞，岳西杰，葛玺祖，等. 不同耕作措施对（土娄）土土壤有机碳形态及活性的影响. 干旱地区农业研究，2010，28（6）：58—63 Wang C X，Yue X J，Ge X Z，et al. Effect of different cultivation measures on activity and bound forms of organic carbon in Lou Soil （In Chinese）.Agricultural Research in the Arid Area，2010，28（6）：58—63.

［27］ 黄玉麟，郭治兴，吴志峰，等. 土壤有机碳储量结果精度分析——以广东省为例. 中国农学通报，2014，30（17）：300—307 Huang Y L，Guo Z X，Wu Z F，et al. Accuracy analysis of soil organic carbon—A case study of Guangdong Province，China （In Chinese）. Chinese Agricultural Science Bulletin，2014，30（17）：300—307

［28］ 陈曦. 广西土壤有机碳储量估算及与全国部分省区的比较研究. 地理科学，2014，34（10）：1247—1253 Chen X. Estimation of soil organic carbon reserves in Guangxi and comparion study with some provinces in China （In Chinese）. Scientia Geographica Sinica，2014，30（10）：1247—1253

［29］ 陈丽娜. 武夷山山地草甸土成土因素分析. 辽宁师专学报（自然科学版），2005，7（3）：91—92 Chen L N. Analysis of soil factors of mountain meadow soil in Wuyishan （In Chinese）. Journal of Liaoning Teachers College（Natural Science Edition），2005，7（3）：91—92

［30］ 郭乃嘉，史学正，赵永存，等. 人为与环境因子对农田土壤有机质影响的比较研究——以典型黑土区和水稻土区为例. 土壤学报，2016，53（5）：1097—1106 Guo N J，Shi X Z，Zhao Y C，et al. Comparative study on impacts of anthropogenic and environment factors on soil organic matter：A case study of typical black soil region and paddy soil region （In Chinese）.Acta Pedologica Sinica，2016，53（5）：1097—1106

［31］ 曹樱子，王小丹. 藏北高寒草原样带土壤有机碳分布及其影响因素. 生态环境学报，2012，21（2）：213—219 Cao Y Z，Wang X D. Spatial distribution of soil organic carbon and influencing factors in the alpine steppe of northern Tibet （In Chinese）. Ecology and Environment Sciences，2012，21（2）：213—219

［32］ 邢世和，张黎明，周碧青. 福建农用地利用区划. 北京：中国农业科学技术出版社，2012 Xing S H，Zhang L M，Zhou B Q. Regionalization of agricultural land use in Fujian （In Chinese）.Beijing：China Agricultural Science and Technology Press，2012

［33］ Raffy M. Change of scale theory： A capital challenge for space observation of earth. International Journal of Remote Sensing，1994，15（12）： 2353—2357