NPP增长驱动下的中国森林生态系统的碳汇及效率摘要：森林生态系统能够有效地吸收大气中的CO2，从而在一定程度上缓解全球变暖的压力生态系统固碳能力取决于两个关键因素：净初级生产力（NPP）的增长强度与生态系统碳周转时间本文基于校正后的遥感监测到的中国森林生态系统NPP增长趋势，结合森林样地实测数据得到的碳分配系数与周转时间，建立了中国森林生态系统碳周转模型，并模拟了1982〜1999年NPP增长驱动下的森林生态系统的碳汇量及其年际变化结果表明：在1982〜1999年间，我国森林生态系统的平均碳汇量为0.052PgCa-1，其中植物的碳汇量为0.034PgCa1，土壤的碳汇量为0.018PgCa-1;在不同森林类型中，常绿针叶林和常绿阔叶林的碳汇贡献最大，而落叶针叶林和针阔叶混交林贡献最小；进一步分析表明森林的固碳效率（碳汇量/NPP增长总量）显著地受到森林碳周转时间的控制关键词：NPP周转时间碳汇效率森林中国CapacityofforestcarbonsequestrationdrivenbyNPPincreasinginChinaAbstract:ForestecosystemcouldsignificantlysequestratesomeatmosphericCO2and,therefore,partlymitigatecurrentpressureonglobalwarming.ThecarbonsequestrationcapacityofforestecosystemisdeterminedbyboththeNPPincreasetrendandturnovertime.InordertoestimatethecapabilityofforestCsequestrationinChina,acarbonturnovermodel,whichbasesonNPPincreasetrendmonitoredbyremotesensingandcarbonturnovertimederivedfromforestobserveddata,wasdesigned.ModeledresultsillustratedthatChinaforestisanapparentcarbon-1-1sinkwithamagnitudeof0.052PgCa,inwhichabout0.034PgCainplanttissuesandtheother-1of0.018PgCainsoil.Thefurtheranalysisoncarbonsequestrationefficiency(CSE)indicatedthattheCSEvalueiscontrolledbycarbonturnovertime.Keyword:NPP,turnovertime,carbonsink,efficiency,forest,China1引言工业革命以来，化石燃料的使用及土地利用方式的改变向大气中释放了大量的温室气体CO2,使得全球气温持续升高，并进而引起一系列的环境问题（Solomonetal.,2007）。
如何有效地缓解大气CO2浓度的增长速度不仅是一个科学问题，而且是一个政治经济问题（Youngetal.,2003）由于森林的再生长与增强生长能够有效地吸收大气中的CO2，因而森林固碳受到多方面的广泛关注（方精云等，2007），它是缓解全球变暖的一项低成本高收益的有效举措（Andrasko,1990;Brownetal.,1996）随着《京都议定书》的签订和各个国家应对全球变暖政策的出台，森林碳汇的研究与应用受到了空前的重视（Cannell,1999;戴民汉等，2004）森林生态系统通过光合作用（碳吸收过程）与呼吸作用（碳排放过程）与大气进行碳交换，当生态系统吸收的碳量大于排放的碳量时，森林生态系统就成为大气C02的汇（简称碳汇），反之则为碳源影响生态系统的碳汇（或碳源）因素很多，如气候变化（Dai,1993）、C02施肥（Cramer,2001）、N沉降（Holland,1997）以及土地利用变化（Houghton,1999）等由于生态系统本身具有很大的空间异质性，加之影响因素各不相同，使得当前生态系统碳汇研究存在很大的不确定性（Denmanetal.,2007）从中国森林生态系统碳汇已有的研究看，无论是基于过程模型的模拟方法（Caoetal.,2003;Wangetal.,2007），还是基于森林资源的清查方法（刘国华，2000；方精云，2001，2007），它们在碳汇的估计上都存在很大的差异。
森林生态系统固碳能力大小取决于两个关键因素：NPP的增长强度与生态系统碳周转时间（Luoetal.,2003）森林再生长和增强生长导致更多的碳进入到生态系统中（表现为NPP增大），因而它是生态系统固碳的外部驱动因素但这部分碳能否有效固定在生态系统中，也就是说生态系统的固碳效率如何，则取决于生态系统的碳周转时间（Luoetal.,2003）当NPP增长驱动相同时，碳周转时间越长，森林的固碳能力越强，固碳效率也越大本研究基于NPP增长驱动下的森林生态系统碳周转模型，对1982〜1999年森林生态系统碳汇年际变化进行了模拟，估算了森林植被与土壤固碳能力，并探讨了碳周转时间对森林固碳效率的影响2方法与数据2.1模型结构NPP增长驱动下的碳周转模型结构如图1所示对于每一种森林生态系统，碳周转模型考虑乔木、灌丛、草本3个层校正后的遥感监测到的npp（Nppm（；a））按不同的分配系数（:t、:sh:-h）分配给乔木、灌丛与草本三个碳库，其中进入到乔木层的NPP又进一步分配给茎、叶、根三个子碳库（as、al、ar）进入到生态系统的有机碳一部分会保留下来并形成各自的生物量（qt、qsh、qh、q$、q），另一部分则通过植物的凋落过程脱离原来的碳库，其凋落速率取决于各碳库的周转时间（.t、.sh、.h、.s、.l、.J。
植物凋落物进入到凋落物与土壤有机碳库（qsoc）中，其中一部分有机碳保留下来形成土壤碳汇，而另一部分则以异养呼吸的形式释放回大气中，其速率取决于土壤碳周转时间（.soc）图1NPP增长驱动下的森林生态系统碳周转模型考虑到自20世纪70年代末开始，中国政府开始实施大规模的生态恢复工程(Shen,2000)，使森林生态系统的生产力大大增加同时Fang(2001)等的研究也表明中国森林固碳量增长开始于19世纪80年代因此，本文假设1982年为森林碳平衡点，并基于(1)〜(6)式对1982〜1999森林碳汇年际变化及累积量进行模拟如二tNPP—qt/tdtdqshdtdqh「shNPP-qsh/shdtdqsdtdql「tnpp-qs/sdt(1)(2)(3)(4)(5)t：rNPP7r/-r（6）dt2.2参数确定NPP增长趋势及强度是驱动生态系统碳汇的外部驱动因素，确定森林生态系统NPP的量及变化趋势是模拟生态系统碳汇的关键之一基于遥感数据的光能利用率模型（CASA模型）由于能较好地模拟区域尺度NPP的空间分布及变化，它们广泛应用于各种空间尺度的NPP监测（Potteretal.,1993;Piaoetal.,2005;朱文泉，2006）。
在CASA模型中，NPP是植物吸收的光合有效辐射（APAR）、最大光能利用率（；*）、温度胁迫系数（T；）及水分胁迫系数（W；）的函数，它可以由下式表示：NPP二fAPARPAR；*T；W；（7）其中fAPAR为光合有效辐射（PAR）被植物吸收的比例，它通过遥感数据NDVI来估算对CASA模型估算的NPP精度影响最大的参数是最大光能利用率；*（彭少麟，2000），它通常随植被类型的不同而有所差异（Ruimyetal.,1999;朱文泉，2006）因此，确定不同植被类型的最大光能利用率参数对于确定NPP增长的年际变化的绝对量（gCm-2a-1）非常重要然而，由于植被类型在空间分布上的年际变化很小，因此，Z的取值对NPP年际变化相对量（％）或趋势的影响不大（Piaoetal.,2005）正因为如此，Piao等（2005）禾9用CASA模型成功地模拟了最大光能利用率为常数（；*=0.405gCMJ-1）下的1982〜1999年中国不同生态系统的NPP增长趋势由于驱动森林碳汇变化的因素是NPP增长的绝对量而非相对量，因此，本文在利用Piao等（2005年）研究的NPP增长趋势的基础上，结合罗天祥等（1996年）1266块森林样地中NPP的观测值对CASA模型的最大光能利用率进行了校正（；；），并进而获得校正后的NPP（NPPm（；；））及年际变化量以驱动碳周转过程模型（图1）：(8)NPRNPPm(*)ttr*4*其中，E为最大光能利用率常数（0.405gCMJ-）,N®0为样地观测的NPP,NPPmQ）为最大光能利用率常数（E*）下模拟的多年平均NPP。
在本研究中，图1中各子库的碳分配系数和周转时间通过观测数据来获取其中碳分配系数由罗天祥的样地观测数据直接得到，而各子库的周转时间则基于罗天祥和王绍强等观测的植物与土壤的碳储量和碳通量数据进行估算（罗天祥，1996;Wangetal.,2003）3结果与讨论3.1最大光能利用率校正后的最大光能利用率（表1）表明，不同森林类型的最大光能利用率差异较大，其中常绿阔叶林与落叶阔叶林的最大光能利用率最大，分别达到了0.896和0.853,这与朱文泉等（2006）的研究结果比较一致（其值分别为0.985和0.692）这一结果也与彭少麟等（2000）的研究较为接近，其研究表明广东的最大光能利用率达到1.25gCMJ-1，大于CASA模型默认的光能利用率常数这表明校正后的最大光能利用率参数在保证NPP增长趋势与Piao等（2005）研究一致的同时，能显著降低CASA模型模拟的NPP增长量与观测结果的差异表1中国主要森林生态系统的最大光能利用率代码森林类型样本数土*；；a1常绿阔叶林2600.4050.8962落叶阔叶林3010.4050.8533常绿针叶林6130.4050.8464落叶针叶林480.4050.6955针阔混交林220.4050.6393.2森林碳汇量及年际变化中国森林生态系统碳汇模拟结果表明，NPP增长及趋势与森林碳汇的关系密切。
遥感监测到的NPP具有明显地增长趋势（图2a），使得进入到生态系统中的NPP的累积量显著增加（图2b），从而直接导致了森林生态系统碳汇量的年际变化（图2c）及累积量上的变化（图2d）因此，NPP增长是驱动中国森林生态系统碳汇变化的根本因素在不同森林类型中，常绿针叶林和常绿阔叶林由于同时具有最大的面积（分别占总面积46.5%和22.5%）和较大的NPP增长强度，因而其碳汇贡献最大，而落叶针叶林和针阔叶混交林贡献最小〔Ed*3360ISalMdcH.y94SIS.N-T0图2森林生态系统碳汇年际变化及累积量从1983~1999年，森林生态系统的累积的碳汇量为0.876PgC,年平均碳汇量为0.052(PgCa-1)，其中固定在植被中的碳汇量为0.034(PgCa-1)，固定在土壤中的碳汇量为0.018(PgCa-1)对于植被固碳来说，乔木的碳汇量最大(0.032PgCa-1)，灌丛次之(0.002PgCa。