森林生态系统净初级生产力模拟研究
0引言
人类社会发展、科技进步的历史同时也是一部干预和破坏环境、改变全球气候的历史。自工业革命以来，随着现代工业的迅猛发展和矿物燃料的广泛利用，以及陆地植被破坏等人类活动引起的地球大气中的“温室气体”以前所未有的速度增加，特别是CO2浓度的急剧增加等带来了全球变暖等严重问题[1]，已经引起了全球的广泛关注，而森林作为陆地生态系统的主体，与其他陆地生态系统类型相比，森林植被具有最广泛的分布面积、最高的生产力和最大的生物量累积，其在全球碳平衡中的巨大的贡献也日益被人们重视。森林生态系统的主体地位主要体现在地圈、生物圈的生物地球化学过程中起着重要的“缓冲器”和“阀”的功能，在陆地碳循环中发挥着重要的作用。
自20世纪80年代以来，随着对森林生态系统结构、功能和生态过程认识的不断深入以及遥感、地理信息系统（GIS）和计算机技术的发展，森林生态系统碳循环模型研究渐渐成为森林碳循环研究中的主流方向之一，并已在区域或全球尺度上对森林生态系统当前及未来的碳存储模式做出了一些预测。其森林碳循环模型按空间尺度可分为斑块尺度的森林碳循环模型和区域尺度的陆地碳循环模型[2]。本文在对此两大类的模型对比分析的基础上，利用英国爱丁堡大学Williams等[2]提出的斑块尺度SPA（Soil-Plant-Atmosphere）森林碳循环模型，以栅格数据为主，其他数据资料为辅等多源数据将该模型从空间尺度上进行了扩展，并利用机群图像并行技术来解决扩展后模型计算量大、计算时间过长的缺点，基于上述条件建立了新的森林生态系统碳循环RSPA模型（RegionSoil-Plant-Atmosphere），通过该模型对2004年三明市的森林生态系统多种森林植被的生产力做出了评估、预测与模型结果验证。
1 SPA模型简介
SPA模型是植被冠层尺度生理生态学的过程模型[3-6]，它是以光合作用-气孔导度-蒸腾作用耦合模型为基础，涉及生理、生化和物理等机理，结合了生态学、植物生理学、气象学和水文学等多种自然科学方法来模拟森林的呼吸、光合、水量平衡、碳的分配之间的关系。在该模型中，植被的冠层被分为若干层次，并且规定植被之间的能量传输必须通过冠层，其冠层层次是具有垂直结构的，不同垂直高度上的植被生理生态学特性都不相同。它是通过逐层计算各通量，最后累加为冠层水平通量。该模型根据植被的生理、物理过程分为两个大的子模块。物理子模块主要包括冠层辐射传输、叶片边界层传导、地表能量平衡、土壤水热传输、根系水吸收等子模块；生理模块包括呼吸作用、光合作用、气孔导度、植被水力等子模块。SPA模型通过对各子模块模拟循环耦合进行森林生态系统生产力的模拟与预测，该循环的关键过程如碳循环、水循环和氮循环过程的生物地球化学过程，主要包括冠层截流、蒸发、植被蒸腾、光合作用、生长和维持呼吸作用、碳在地上与地下器官中的分配、凋落、分解及氮
硝化过程。还可以分析斑块研究区内森林生态系统的总初级生产力（GrossPrimary Productivity，GPP）、土壤水分蒸发蒸腾损失总量（Evapotranspiration，ET）的季节和年际变化，从而研究气候变化对森林生态系统的影响。目前，这类植被冠层尺度生理生态学过程模型已成为森林生态系统碳循环模型的核心之一。
2研究区域介绍
福建省三明市位于25°29′N~27°07′N，116°22′E~118°39′E。该市森林资源丰富，是中国南方集体林区综合改革试验区，享有福建“绿色宝库”的美誉，林木品种繁多，森林覆盖率达76.8%，活立木蓄积量1.15亿m3（为全福建省的1/3），毛竹储量3.8亿株，植物种类多达1 300多种，其中有秃杉、南方红豆杉等20多种国家重点保护植物和9种省级重点保护植物。森林植被类型主要有常绿阔叶林、常绿针叶林、针阔混交林、竹林、经济林、灌木林等。
3森林碳循环模型研究实验方法
3.1试验主要数据
本文中对所有的空间图像数据的坐标系统和空间分辨率大小进行统一，主要是利用MODIS产品的250 m×250 m分辨率，数据的坐标系统参数如表1所示。
3.2 SPA模型驱动参数变量空间化
SPA模型是一个只能计算斑块尺度的单一森林植被类型的森林生态系统碳循环模型，无法用于大区域森林生态系统碳循环研究，而在本研究中却需要该模型空间化以面状展开，因此必须将SPA模型各种参数空间化，将原来SPA模型的输入数据从点数据到面数据进行转换，这就需要空间数据的收集与处理，此工作是本文整个研究工作的基础，数据的质量将直接影响模拟结果的精度。本研究利用的空间数据有土地利用数据、植被类型数据、经纬度数据、风速数据、降水数据、温度数据、动力效应模型（Dynamic Effect Model，DEM）数据（图1~图4和表2~表5，表3中的面积百分数是指各土壤类型占整个三明市地区土地的面积百分比）等，以此作为空间化数据的接口，其中DEM数据为NASA的航天飞机雷达拓扑测绘（Shuttle Radar Topography Mission，SRTM）数据；土壤类型数据、叶面积指数（Leaf Area Index，LAI）数据、森林植被数据由福州大学空间数据挖掘与信息共享教育部重点实验室提供。本文利用空间插值方法对风速数据、降水数据、温度数据插值出每天24 h的各种所需数据，将这些数据作为SPA模型的空间输入数据源之一。位时间内所固定的有机物质总量，是由GPP中扣除自养呼吸量Ra后的剩余部分[6]NPP=GPP·Ra（1）当计算完GPP后，需要建立植被自养呼吸模型来求得森林生态系统自养呼吸部分，以求得NPP。通常规定森林植被的Ra为维持性呼吸Rm和生长性呼吸Rg之和，即Ra=Rm+Rg（2）对于Rg，实际计算时主要采取先主后分的方法：先考虑总的生长呼吸占总光合量的比例，然后再考虑森林植被各组份的生长呼吸在总生长呼吸量中所占的比例。即Rg=γ·GPP（3）Rg,j=Rg·ra,j（4）其中，γ为植被总生长呼吸占总光合量的比例，ra,j为植物各组份的生长呼吸比例。对于Rm，
