森林在生态系统碳循环过程中发挥着十分重要的作用。森林碳储量是研究森林与大气之间碳交换的基本参数，也是计算碳排放和碳循环过程中的关键因子^[1-2]，因此研究森林碳储量的估算方法具有十分重要的意义。近年来，运用以森林蓄积量为基础的生物量转换因子法(BEF)来估测植被碳储量得到了广泛应用^[3]。方精云等建立了“蓄积量-生物量”回归方程，并测算了不同优势树种的计算参数^[4-6]，徐新良^[7]、曾伟生^[8]等又通过研究进一步提出了方程中参数a和b的优化取值。李海奎等以全国第六、七次森林资源连续清查资料为基础，研究比较了生物量转换因子法同其他方法在森林生物量和碳储量估算方面的异同^[9]，为提高生态系统植被碳储量估算精度打下了基础，实现了森林植被生物量和碳储量估算方法由样地调查向区域推算尺度的转换。

传统蓄积量的调查方法主要通过周期性测定样地内的树木生长指标来估算^[10]，这种方法精度较高，但野外工作量大，人力和资金成本较高，且很难满足大尺度以及实时的监测需求^[11]。随着遥感技术的发展，利用3S及其集成技术开展森林蓄积量的定量估测受到越来越多的关注^[12]。遥感数据能提供动态的植被覆盖信息，与森林清查数据相结合，能实现快速、连续、较准确的蓄积量估测及预报，极大提高了工作效率，诸多国内外学者在此基础上对森林蓄积量定量遥感估测开展研究^[13-14]。一些研究结果表明，对遥感影像的灰度值提取，其线性或非线性组合与森林蓄积量有较明显的相关性，以地面调查数据得到的蓄积量为因变量，以遥感数据若干波段的灰度值及其组合值作为自变量，同时综合考虑坡度、海拔等地理信息，可以建立起合适的森林蓄积量多元线性回归模型^[15-17]。

当前人类活动对陆地生态系统碳循环过程的影响越来越强烈，尤其是土地利用变化显著地影响了陆地生态系统的结构和功能，造成植被和土壤碳储量的变化，是对陆地生态系统碳循环影响的最大人为因素之一^[18-20]。研究表明，1850—1980年间，约180Pg碳通过土地利用变化从陆地生态系统中排放到大气中，而森林面积减少导致近100Pg碳排放到大气中^[21]。20世纪80年代以来，全球由于土地利用变化造成每年平均释放1.6Pg碳^[22]。因此，定量估算森林碳储量及其变化、评估碳储量对土地利用变化的响应不仅可以为森林资源经营管理和林业可持续发展提供重要科学依据，而且对碳排放及碳循环研究具有十分重要的意义。

本文主要对以下问题展开研究：(1) 基于RS与GIS技术，从Landsat5 TM影像中提取不同波段灰度值及其组合，结合海拔、坡度等因子作为自变量，以福州市森林资源二类调查中的单位面积蓄积量数据为因变量，采用逐步回归法建立福州市森林蓄积量多元线性估测模型。(2) 根据生物量转换因子法，对福州市森林植被碳储量和碳密度进行估算，分析2000年到2010年森林植被碳储量变化特征。(3) 建立2000年至2010年福州市土地利用转移矩阵，并分析土地利用变化影响下的森林碳储量变化特征。研究结果旨在为福州市利用遥感影像数据监测及预测连续清查间隔期内森林资源的动态变化提供科学依据。

1 研究区概况

福州市位于福建省中部东端、闽江下游，介于北纬25°15′—26°39′、东经118°08′—120°31′之间，总面积11968 km²，其中市区面积1786 km²。行政区划为五区七县(市)，包括市区(鼓楼区、台江区、仓山区、晋安区和马尾区)、闽清县、永泰县、闽侯县、连江县、罗源县、福清市、长乐市。

福州属亚热带季风气候，东临台湾海峡，温暖湿润，雨量充沛，四季常青。福州市横跨中国东部湿润森林区域中的两个植被带，即南亚热带雨林植被带和中亚热带照叶林植被带。福州天然植被包括常绿阔叶林、针叶林、针阔混交林、竹林和灌丛。常绿阔叶林为中亚热带地带性植被，植物种类以壳斗科栲属、石栎属为主，其它樟科、山茶科、杜鹃花科植物也占较多比例。针叶林为福州市分布最广、数量最多的植被类型，主要树种有马尾松、杉木、油杉、柳杉、黄山松、黑松等，其中马尾松、杉木覆盖面积最大，且多为纯林，基本上为人工营造。

2 材料与方法 2.1 数据来源及预处理

本研究采用Landsat5 TM遥感影像数据、福州市土地利用数据、福州市森林资源二类调查数据、DEM数据作为基本数据源。

2.1.1 Landsat TM影像数据

TM遥感影像来源于2000年和2010年，均为植被生长季数据，单景影像覆盖全市范围，平均云量小于1%，噪音面积小于10%，影像质量较好，适合进行遥感反演。为削弱不利因素对预测模型精度的影响，预先对TM影像数据进行大气校正、正射校正、几何校正。校正过程在ENVI 5.0软件中进行。

2.1.2 土地利用数据

福州市土地覆被数据是利用遥感影像等资料，通过计算机辅助分类软件和人工目视解译相结合生成，共2000年、2010年两期数据。森林类型分类至二级生态系统，即常绿阔叶林、常绿针叶林、针阔混交林，并进行过野外核查修正，精度大于90%。

土地利用分类体系采用政府间气候变化和专门委员会(Intergovernmental Panel on Climate Change,IPCC)编著的《土地利用、土地利用变化和林业优良做法指南(GPG-LULUCF)》^[23]中使用的土地类别，即森林、耕地、草地、湿地、建设用地和未利用地，是目前公认的可以较全面、综合地评价土地利用对碳储量影响的土地分类体系。在此基础上，本文根据研究需要，增加灌木类别，共计七种土地利用类型。为和“灌木”区分，本研究中的“森林”土地利用类型特指乔木林。

2.1.3 地面调查数据

地面调查数据来源于2010年福州市森林资源二类调查数据(下文简称“二调数据”)，森林小班为统计的基本单元，以矢量地理信息数据的形式储存。在进行研究前对数据进行预处理，针对小班出现自相交的情况，在ArcGIS 10.2软件中使用Repair Geometry工具进行修复。

2.2 自变量因子选择及提取

从遥感数据和DEM数据中提取适用于构建蓄积量定量估测模型的因子为备选自变量，具体包括遥感因子变量、植被因子变量、地理因子变量。

遥感因子变量根据TM影像7个波段具有不同的波谱特征，将TM1-TM5、TM7六个单波段设置为遥感因子。TM6为热红外波段，本次研究不予考虑。除此之外，在借鉴前人研究的基础上，设置TM5/TM4、TM5/TM7、TM7/TM3、TM3/(TM1+TM2+TM3+TM4+TM5+TM7)、(TM4+TM5-TM2)/(TM4+TM5+TM2)、(TM5+TM7-TM2)/(TM5+TM7+TM2)、TM3×TM4/TM7七个组合波段。此外由于缨帽变换中的土壤亮度指数(BI)、植被绿度指数(GI)、湿度指数(WI)3个分量与地物类型有密切关系，故亦采用作为遥感因子变量^{[17, 24]}。

植被因子变量采用当前比较常用的、能反映植被覆盖信息的重要植被指数，主要有归一化植被指数NDVI、差值植被指数DVI、比值植被指数RVI。各植被指数在反映植被长势、生物量等方面具备各自优势。遥感因子变量和植被因子变量用ENVI 5.0软件进行提取。地理因子变量主要为海拔和坡度，从DEM数据中直接用ArcGIS 10.2软件提取。

因变量来源于二调数据。森林小班根据林业管理需要呈现出不规则的面积形态，一个森林小班中往往包含许多个遥感灰度值不同的像元，因此需要在不规则的小班面积覆盖范围内提取唯一的像元灰度值。将小班矢量图层和提取出的每个自变量因子栅格图层在ArcGIS 10.2软件中进行叠置，利用Spatial Statistics Tools做统计分析，可得到小班内某个自变量因子所有像元的平均值，因而将任意小班内的自变量因子转换成唯一值。

2.3 样本数据筛选

运用标准差分析法对小班异常离群数据进行剔除，公式为：|x_i-x|>2σ。式中x_i为实测小班单位面积蓄积量数据，x为平均值，σ为标准差。筛选过后的小班个数为常绿阔叶林1565个；常绿针叶林2722个；针阔混交林1045个。考虑到样本数据的数量需满足建模要求，最终每个森林类型随机选取400个小班数据作为建模样本，选取150个小班数据作为验证样本。各森林类型的建模数据和验证数据样本点分布如图 1所示。

2.4 逐步回归法

自变量个数及种类对模型估测精度的影响很大，需要克服备选变量间的共线性，筛选合适的自变量构建模型^[25]。在进行多元线性回归时，一般线性回归不论自变量对因变量的影响是否显著，都可以进入方程，这样会使方程估计的精度变低^[26]，因此本研究采用逐步回归法避免上述问题。

逐次将2.2节中提取的21个备选变量选入方程，比较前后变量对方程的重要性。若后引入的变量使先前变量失去了重要性，则将先前的变量淘汰，通过变量的引入和剔除，最终筛选出最适合建模的自变量^[27]。逐步回归中的每一步都进行F检验，以确保每次引入新的变量之前回归方程中只包含显著变量。以回归系数显著性检验中各自变量的F统计量的相伴概率值Sig≤0.05和F统计量的相伴概率值Sig≥0.10作为自变量是否引入模型或者从模型剔出的标准。当Sig≤0.05，认为该变量对因变量的影响是显著的，应被引入回归方程；当Sig≥0.10时，则认为该变量对因变量的影响是不显著的，应从回归方程剔除。

2.5 基于蓄积量的森林碳储量估算

本研究采用BEF法估算森林植被生物量和碳储量，该方法建立生物量与蓄积量之间的换算关系，可以实现由样地调查向区域推算的尺度转换，在国内城市森林的碳储量估算及碳汇评估中较为常见^[28]，其公式^[5]为：

式中，B为森林生物量(t)，V为森林蓄积量(m³)，a、b为参数，取值参照前人的研究成果(表 1)。

森林植被碳储量与碳密度计算采用公式：C=B×C_C;ρ=C/S。式中C为碳储量(t)；C_C为含碳率，采用国内学者常用的系数0.5^[30]；ρ为碳密度(t/hm²)；S为面积(hm²)。

3 结果与讨论 3.1 蓄积量定量估测模型

经过逐步回归法对自变量因子的反复引入和剔除，最终筛选出适合于常绿阔叶林、常绿针叶林、针阔混交林建模的自变量因子，建立多元线性回归方程(表 2)。

α=0.05的可靠性水平下，各模型均通过F检验和t检验。常绿阔叶林、常绿针叶林、针阔混交林的回归模型中的复相关系数分别为0.790、0.803、0.835，修正决定系数分别为0.599、0.679、0.694，能较好地反映出蓄积量与建模因子之间的线性关系，拟合效果较好。

3.2 模型适用性检验及精度验证

对所拟合的模型进行适用性检验，用本研究建立的模型所估测出的蓄积量作为自变量，以小班实测蓄积量数据为因变量，建立二者的线性回归方程。3种不同森林类型的检验方程及R²如图 2。

常绿阔叶林、常绿针叶林、针阔混交林的模型检验方程R²分别为0.881、0.969、0.962，结果表明，通过模型计算的林分蓄积量估测值与森林小班中的实测值无显著差异。

此外，对所拟合的模型进行精度验证。2010年福州市森林蓄积量实际统计值为2951万m³，本研究中通过模型的估测值为2593.95万m³，总体精度为2593.95/2951=87.90%，预测结果较为理想。

3.3 森林植被碳储量变化分析

采用2.2节同样方法提取2000年不同森林类型的建模因子。在ArcGIS软件中将参与蓄积量估测的因子变量栅格图分别用土地覆被矢量数据中的常绿阔叶林、常绿针叶林、针阔混交林掩膜裁剪，根据已建立的蓄积量估测模型和TM遥感影像反演值计算每个像元的蓄积量。根据2.5节的BEF法，得出2000年、2010年福州市森林植被碳储量总量分别为：12.499Tg、12.642Tg，碳密度分别为18.694 t/hm²、18.708 t/hm²，森林植被碳储量增加了1.44×10⁵t，平均碳密度为18.701 t/hm²。在ArcGIS软件支持下得到福州市2000年、2010年森林植被碳密度分布图(图 3)，并对十年间福州市各市县森林植被碳密度空间变化进行分析。

由图 3可以看出，从2000年到2010年，福州市森林植被碳密度总体呈上升的趋势。在东部沿海区域，植被碳密度大于25 t/hm²的像元总面积减少，在西部内陆地区，森林植被碳密度大于25 t/hm²的像元总面积增加。福州市的各市县中，闽清县、永泰县、闽侯县的森林植被碳密度常年保持较高水平，并呈现出增长趋势。市区(包括鼓楼区、仓山区、台江区、晋安区、马尾区)的植被碳储量平均值为18.91 t/hm²，变化较小，总体上保持稳定。罗源县、长乐市、连江县森林植被碳密度较低，并呈现下降趋势。闽清县、永泰县、闽侯县地处福州市西部，地广人稀，经济欠发达，多山脉多林场，并分布有多个自然保护区和森林公园，森林生态系统受到较好保护。该区域受人为负面干扰程度较小，森林植被碳储量常年保持较稳定增长的水平。处于沿海地区的各市县，如市区、长乐市、罗源县、福清市、连江县，随着社会经济的发展和建设用地对森林的侵占，森林植被碳密度有下降的趋势。

3.4 福州市土地利用变化分析 3.4.1 土地利用类型面积变化

2000年到2010年福州市土地利用面积变化见表 3。福州市主要用地类型为森林，常年保持国土面积占比55%以上，其次为耕地和灌木。10年间，福州市土地利用面积发生较大变化，主要表现在：建设用地、森林、草地面积上升，其占福州市国土面积比例分别上升1.32%、0.60%、0.05%；灌木、耕地、湿地面积下降，占福州市国土面积比例分别下降1.25%、0.52%、1.05%。未利用地面积变化较小，变化仅0.01%。

3.4.2 土地利用类型转化

在ArcGIS软件进行数据融合和叠置分析，得到土地利用转化关系图，将数据导出后在Excel表格中计算土地利用转移矩阵，结果见表 4。

2000年到2010年，土地利用类型转化较明显。灌木和耕地是主要的转出者，灌木主要转化为森林，转化面积113.69 km²；耕地主要转化为建设用地和森林，转化面积分别为109.59 km²和61.38 km²。森林和建设用地是主要的转入者，森林主要由灌木和耕地转化，建设用地主要由耕地和湿地转化。森林主要向建设用地、耕地进行了转化，转换面积为14.67 km²和88.85 km²。可以看出经济发展和城市扩张导致对建设用地的诉求强烈，耕地转换为建设用地的强度增大，耕地安全需得到保障。森林面积的持续上升对森林资源的健康可持续发展有利。

在森林和其它土地利用类型的相互转化过程中，耕地和灌木转化最为剧烈。森林向耕地转化了88.75 km²，同时耕地向森林转化了61.38 km²，相互转化幅度较大；灌木向森林转化了113.69 km²，在所有转化过程中面积最大，而森林向灌木只转化了1.59 km²。草地、湿地和森林之间都发生面积转化，但变化幅度均较小。森林向建设用地转化14.67 km²，而建设用地只向森林转化了0.11 km²。

3.5 碳储量对土地利用变化的响应

根据土地利用类型转移矩阵及不同土地利用类型土壤碳密度和植被碳密度差异，估算福州市2000年至2010年森林碳储量变化对土地利用类型变化的响应。森林植被碳密度采用本文计算的平均值，其它不同土地利用类型碳密度来源于当地资料和相关文献^[31-34](表 5)。土地利用变化带来的森林碳储量变化见表 6。

由于土地利用变化，2000年到2010年福州市森林总碳储量减少了1.711×10⁴ t，其中土壤碳储量减少2.23×10³t，植被碳储量减少了1.489×10⁴ t，土地利用变化总体上造成了碳排放。森林向其他用地类型转化表现为碳排放效应，总碳储量减少1.49×10⁵ t。排放最多的为森林向耕地转化，总碳储量减少1.17×10⁵ t，占减少碳储量总量的79.07%；其次依次为森林向建设用地、草地、湿地、未利用地、灌木转化，碳储量减少量分别为2.465×10⁴、4.90×10³、8.6×10²、7.7×10²、6.1×10² t。其他用地类型向森林的转化表现为碳吸收效应，总碳储量增加1.32×10⁵t，吸收最多的为耕地向森林的转化，总碳储量增加8.10×10⁴t；其次依次为灌木、草地、湿地、建设用地向森林转化，碳储量增加量分别为4.32×10⁴、5.02×10³、2.39×10³、1.90×10²t。耕地和森林之间的土地利用转换造成的碳储量变化十分剧烈，建设用地对森林的侵占造成的碳储量减少应该引起相应的重视。

4 结论

(1) 不同的森林类型，即常绿阔叶林、常绿针叶林、针阔混交林蓄积量估测模型的复相关系数分别为0.599、0.679、0.694，模拟效果较好，具有较高的估测精度。

(2) 2000年、2010年福州市森林植被碳储量总量分别为：12.499、12.642Tg，碳密度分别为18.694、18.708 t/hm²，森林植被碳储量增加了1.44×10⁵ t。闽清县、永泰县、闽侯县的森林植被碳密度常年保持较高水平，并呈现出增长趋势。罗源县、长乐市、连江县森林植被碳密度较低，并呈现下降趋势。

(3) 2000年到2010年，福州市主要用地类型为森林，其次为耕地和灌木。土地利用面积变化和类型转化较明显。灌木和耕地是主要的转出者；森林和建设用地是主要的转入者。

(4) 由于土地利用变化，2000年到2010年福州市总碳储量减少了1.711×10⁴ t，其中土壤碳储量减少2.230×10³ t，植被碳储量减少了1.489×10⁴ t，土地利用变化总体上造成了碳排放。耕地和森林之间的土地利用转换造成的碳储量变化十分剧烈，建设用地对森林的侵占造成的碳储量减少应该引起相应的重视。