NPP估算CASE模型.docx

1、NPP估算CASE模型第三章 长江上游初级生产力评估植被既是重要的自然资源，又是自然条件（如地质、地貌、气候、土壤等）和人类开发利用资源状况的综合反映。随着全球变化的加剧及其对全球变化研究的不断深入，植被作为陆地生物圈的主体，在生态系统中的作用也日益受到重视，尤其是在全球物质循环、能量流动、调节全球碳平衡、减缓大气中CO2等温室气体浓度上升趋势以及维护全球气候稳定等方面具有不容忽视的意义。植被净初级生产力（Net Primary Productivity，NPP）是指绿色植物在单位面积、单位时间内所累积的有机物数量，表现为光合作用固定的有机碳中扣除植物本身呼吸消耗的部分，这一部分用于植被的生长

2、和生殖，也称净第一性生产力。NPP作为地表碳循环的重要组成部分，不仅直接反映了植被群落在自然环境条件下的生产能力，表征陆地生态系统的质量状况，而且是判定生态系统中碳源/碳汇和调节生态过程的主要因子，在全球变化及碳平衡中扮演着重要的作用。自20世纪60年代以来，各国学者对NPP的研究倍受重视，国际生物学计划（International Biological Programme，IBP，19651974）期间，曾进行了大量的植物NPP的测定，并以测定资料为基础联系气候环境因子建立模型对植被NPP的区域分布进行评估如Miami模型、Thornthwaite纪念模型、Chikugo模型等。建立于198

3、7年的国际地圈生物圈计划（International Geo-Biosphere Programme，IGBP）、全球变化与陆地生态系统（GCTE）和最近出台的京都协定（Kyoto Protocol）均把植被的NPP研究确定为核心内容之一（IGBP，1998）。长江上游地区面积广大，地形复杂，气候多样，植被类型丰富，是我国生态屏障关键区，也是重要的生态脆弱区和气候变化敏感区。同时，它还是我国生物多样性和自然保护集中区和我国重要的森林分布区、草地分布区。长江上游陆地生态系统特殊的地理区位、复杂的资源体系决定着它在我国生态系统中不可替代的战略地位。长江上游山地垂直气候多样，孕育了丰富的森林植被，是

4、我国第二大林区，其陆地生态系统生产力在我国陆地生态系统生产力中具有十分重要的作用。目前，长江上游地区正在大力进行封山育林、公益林和人工林基地建设和实施天然林保护工程、退耕还林工程。可以预计，随着森林面积的不断增加，长江上游地区陆地生态系统初级生产力在我国陆地生态系统初级生产力中的地位也会发生变化。因此开展长江上游植被NPP及其时空分布的研究，对于掌握该区域的生态系统质量状况和自然生产能力具有重要的意义。第一节 NPP研究概述人类研究陆地植被生产力及其地理分布有相当长的历史，有关学者从不同角度及学科对NPP的估算进行了深入细致的研究，取得了丰硕的成果。从空间尺度上来说，可分为NPP定位观测、区域

5、NPP模拟估算和全球NPP模拟估算三种尺度。对陆地植被初级生产力的研究方法主要有2种：野外收割实测法和利用模式模拟计算法。第一种方法就是在单位面积的土地上收割某种植物，晾干后称重计算该种植物的生产力。Wittaker等利用全球各种植被实测数据估算了全球陆地生态系统的净初级生产力为59 Pg/yr；第二种方法就是通过模式中的参数化方案来模拟各种植被的净初级生产力，该方法现已成为估算陆地生态系统生产力的主要手段。自Lieth发表了首张基于回归模型（Mi-ami模型）用计算机模拟的全球NPP分布图，全球陆地生态系统NPP的研究就成为全球碳循环研究的热点领域之一。Cramer对IGBP支持的17个模型

6、统NPP年平均值为54. 9 Pg，其中最高的NPP出现在水热条件良好的热带地区（1000gC/m2/yr），温带地区的NPP居中（500700 gC/m2/yr），NPP最低值出现在寒冷和干旱地区（20gC/m2/yr）。基于地面的NPP定位观测，只能收集到数公顷面积的不同生态系统类型的实测数据，然后根据各种生态系统类型，用以点代面的办法外推区域及全球NPP总量。由于这些估算是基于空间实测数据，迄今仍被用作全球NPP估算的参照。在区域或全球尺度上，人们无法直接和全面地测量NPP，因此利用模型估算NPP已成为一种重要而且被打多数学者广泛接受的研究方法。现有的NPP模型大体分为气候生产力模型、生

7、理生态过程模型、光能利用率模型和生态遥感耦合模型四类。NPP估算模型的优缺点及适用条件如表3.1所示。表3.1不同NPP估算模型优缺点比较表模型类型模型举例优点缺点适用条件气候生产力模型MiamiThornthwaiteChichugo(1)模型简单(2)气候参数易获取(1)生理生态机制不是很清楚(2)估算结果以点代面(3)估算误差较大，是一种潜在的NPP适用于区域潜在NPP的估算生理生态过程模型CENTURYTEMBIOME-BGC(1)生理生态机制清楚(2)可以模拟、预测全球变化对NPP的影响 (3)估算结果较准确(1)模型复杂(2)所需参数较多，而且难以获得 (3)区域尺度转换困难适用于

8、空间尺度较小、均质斑块上的NPP估算光能利用率模型CASA，GLO-PEMSDBM(1)区域尺度转换容易，适宜于向区域及全球推广(2)许多植被参数可由遥感获得(3)可以获得NPP的季节、年际动态(1)缺乏可靠的生理生态基础(2)无法实现NPP的模拟及预测(3)光能传递及转换的的过程中还存在许多不确定性适用于区域及全球尺度上的NPP估算生态遥感耦合模型BEPS，改进的PEM模型(1)遥感数据在获取NPP空间分布信息时得到了有效利用(2)具有模拟、预测功能，可以获得NPP的季节、年际动态(3)植被变化信息能立即反映在NPP估算上(1)BEPS模型比较复杂，所需参数较多，在参数确定上人为因素影响较大

9、(2)改进的PEM模型虽然对生理生态过程作了简化，但NPP估算精度受LAI影响较大。适用于小面积样区、区域及全球尺度上的NPP估算总体上看，气候生产力模型比较简单，在资料欠缺和技术落后的情况下，它的应用比较广泛。生理生态过程模型仍是当前陆地NPP估算研究的主要手段，而区域尺度转换则是它所面临的关键问题。近年来，随着遥感技术的发展，光能利用率模型已成为NPP估算的一种全新手段，它利用遥感所获得的全覆盖数据，使区域及全球尺度的NPP估算成为可能，但其生态学机理还有待于进一步研究。而生态遥感耦合模型则是综合了生理生态过程模型和光能利用率模型的优点，通过LAI将二者整合起来，可以反映区域及全球尺度的N

10、PP空间分布及动态，增强了陆地NPP估算的可靠性和可操作性，表现出巨大的发展潜力。近年来关于中国陆地植被净初级生产力时空特征分布的研究工作也陆续开展开，研究的方法主要采用静态统计模型或改进的光能利用率模型以及过程机理模型来估算我国陆地各种植被的净初级生产力，并且讨论了各种气候要素对于我国陆地植被净初级生产力的影响作用。利用光能利用率模型估算植被NPP有三大优点：1）模型比较简单，可直接利用遥感获得全覆盖数据，在实验和理论基础上适宜于向区域及全球推广；2）冠层绿叶所吸收的光合有效辐射的比例可以通过遥感手段获得；3）可以获得确切的NPP季节、年际动态。因此，近年来，光能利用率模型已成为NPP模型的

11、主要发展方向之一。本研究选择目前国际和国内上应用较多的CASA模型来评价长江上游地区植被净第一性生产力。第二节 研究方法与技术路线1. 研究方法1.1 模拟模型光能利用率模型估算植物的NPP主要是基于资源平衡的观点，即假定生态过程趋于调整植物特性以响应环境条件，认为植物的生长是资源可利用性的组合体，物种通过生态过程的排序和生理生化、形态过程的植物驯化，就趋向于所有资源对植物生长有平等的限制作用。在某些极端或者环境因子变化迅速的情况下，如果完全适应不可能，或者植物还来不及适应新的环境，NPP则受最紧缺资源的限制。由此认为，任何对植物生长起限制性的资源(如水、氮、光照等)均可用于NPP的估算，它们

12、之间可以通过一个转换因子联系起来，这一转换因子既可以是一个复杂的调节模型，也可以是一个简单的比率常数。NPP和限制性资源的关系可用公式表示如下：NPP=FcRu其中，Fc：转换因子，Ru：限制性资源。光合有效辐射（PAR）是植物光合作用的驱动力，对这部分光的截获和利用是生物圈起源、进化和持续存在的必要条件。光合有效辐射是植物NPP的一个决定性因子，而植物吸收的光合有效辐射（APAR）则尤为重要。著名的Monteith方程就是建立在此基础之上。NPP=APAR式中为植物光能利用率，它受水、温度、营养物质等的影响。 1.2 CASA模型与参数获取早在20世纪70年代Monteith就发现NPP和植

13、被吸收的光合有效辐射（APAR）之间存在着稳定的关系：当水分和肥料处在最适的条件下，农作物的NPP与APAR具有很强的线性相关。进一步的研究发现，NPP与APAR的时间序列积分有较好的相关性，但是不同的植被类型，或者同一植被类型在不同的生长条件下，所获得的经验模型存在着差异，这就意味着植被的NPP受植物本身及其生长环境的影响。目前，在区域及全球尺度的NPP估算模型当中，以CASA模型为代表的光能利用率模型得到了广泛应用。本研究所构建的NPP遥感估算模型将植被覆盖分类引入模型，并考虑植被覆盖分类精度对NPP估算的影响，由它们共同决定不同植被覆盖类型的NDVI最大值，由此获得各植被覆盖类型的比值植

14、被指数最大值，最后实现FPAR的估算。根据误差最小的原则，模拟出各植被类型的最大光能利用率，使之更符合长江上游的实际情况。利用气象数据（温度、降水、太阳净辐射），结合已有的区域蒸散模型来实现水分胁迫因子的估算，这样一方面可以保持模型原有的植物生理生态学基础，另一方面则在一定程度上对有关参数实行了简化，使其实际的可操作性得到加强。NPP估算模型的总体设计如图3.2所示。模型中所估算的NPP可以由植物吸收的光合有效辐射(APAR)和实际光能利用率()两个因子来表示，其估算公式如下： 式中： APAR(x,t)表示像元x在t月份吸收的光合有效辐射（单位：MJ/m2/月） (x,t)表示像元x在t月份

15、的实际光能利用率（单位：gC/MJ）1）APAR估算用遥感数据估算光合有效辐射(PAR)中被植物叶子吸收的部分(APAR)是基于植被对红外和近红外波段的反射特征实现的。植被吸收的光合有效辐射取决于太阳总辐射和植物本身的特征，可用以下公式进行计算： 式中： SOL(x,t)表示t月份在像元x处的太阳总辐射量(MJ/m2/月)；FPAR(x,t)为植被层对入射光合有效辐射的吸收比例（无单位）；常数0.5表示植被所能利用的太阳有效辐射（波长为0.380.71um）占太阳总辐射的比例。 对于FPAR的计算采用Potter等提出的如下计算公式:式中:SRmin取值为1.08，SRmax的大小与植被类型有

16、关，取值范围在4.14-6.17之间。SR(x，t)由NDVI(x，t)求得：2）光合利用率估算光能利用率是在一定时期单位面积上生产的干物质中所包含的化学潜能与同一时间投射到该面积上的光合有效辐射能之比。研究表明，不同植被的光能利用率差异很大。引起差异的主要原因是光能利用率受气温、水分、土壤、营养、疾病、个体发育、基因型差异和植被维持与生长的不同能量分配等因素的影响。所以，不同类型的植被以及同一植被在不同的环境条件下，其光能利用率不同。光能利用率的估算流程如图3.2所示： 图3.2 光能利用率估算流程Potter等认为在理想条件下植被具有最大光能利用率，而在现实条件下的最大光能利用率主要受温度

17、和水分的影响，其计算公式为：式中： T1(x,t)和T2(x,t)表示低温和高温对光能利用率的胁迫作用（无单位）；W(x,t)为水分胁迫影响系数（无单位），反映水分条件的影响；max是理想条件下的最大光能利用率（单位：gC/MJ）。T1(x,t)反映在低温和高温时植物内在的生化作用对光合的限制而降低净第一性生产力用下式计算：其中：Topt(x)为某一区域一年内NDVI值达到最高时的当月平均气温。已有许多研究表明，NDVI的大小及其变化可以反映植物的生长状况，NDVI达到最高时，植物生长最快，此时的气温可以在一定程度上代表植物生长的最适温度。T2(x,t)表示环境温度从最适温度Topt(x)向高

18、温和低温变化时植物光利用率逐渐变小的趋势，这是因为低温和高温时高的呼吸消耗必将会降低光利用率，生长在偏离最适温度的条件下，其光利用率也一定会降低。用下式计算：当某一月平均温度T(x,t)比最适温度Topt(x)高10或低13时，该月的T2(x,t)值等于月平均温度T(x,t)为最适温度Topt(x)时T2(x,t)值的一半。水分胁迫影响系数W(x,t)反映了植物所能利用的有效水分条件对光能利用率的影响，随着环境中有效水分的增加，W(x,t)逐渐增大，它的取值范围为0.5（在极端干旱条件下）到1（非常湿润条件下），由以下公式计算：其中：E(x,t)为区域实际蒸散量，可根据周广胜和张新时建立的区域

19、实际蒸散模型求取；Ep(x,t)为区域潜在蒸散量，可根据Boucher提出的互补关系)求取。 式中：P(x,t)为像元x在t月的降水量，Rn(x,t)为像元x在t月份的太阳净辐射量。 式中：Ep(x,t)(mm)为局地潜在蒸散量，可以由Thornthwaite的植被气候关系模型的计算方法求算。月最大光利用率max的取值因不同的植被类型而有所不同，由于全球最大光转化率的取值对NPP的估算结果影响很大，人们对它的大小一直存在争议，彭少麟等利用GIS和RS估算了广东植被光利用率，认为CASA模型中所使用的全球植被月最大光利用率（0.389gC/MJ）对广东植被来讲偏低。本文利用朱文泉等人根据生态生理

20、过程模型BIOME-BGG对10种植被类型所模拟的结果，它们的取值分别为：常绿针叶林0.389gC/MJ，常绿阔叶林0.985gC/MJ，落叶针叶林0.485gC/MJ，落叶阔叶林0.692gC/MJ，针阔混交林0.475gC/MJ，常绿、落叶阔叶混交林0.768 gC/MJ，灌丛0.429gC/MJ，草地0.542g C/MJ，矮林灌丛0.888g C/MJ，草地0.608gC/MJ，耕作植被0.604gC/MJ。其他如城市、水体等生态系统取CASA模型所估算的全球月平均最大光利用率0.542gC/MJ。2. 技术路线 技术路线如下图（图3.3）。文 献 调 研数据收集模型建立NPP计算方法

21、选取长江上游植被净初级生产力估算长江上游植被净初级生产力评估气象数据MODIS数据太阳总辐射量月平均温度月总降雨量NDVI植被类型植被模块APAR模块基于光能利用率的NPP估算模型图3.3 技术路线第三节 长江上游初级生产力估算1数据准备1.1遥感数据：由于NOAA/AVHRR数据具有全球覆盖、重访周期短的特点，在空间和时间上具有连续性，所以在植被、气候及地表长期变化研究中成为非常有效的数据源。本研究所采用的NOAA/AVHRRNDVI数据，来源于美国地球资源观测系统（Earth resources observation system：EROS）数据中心的探路者数据集（Pathfinder

22、AVHRR Land Data Set:PAL Data Set）。图像空间分辨率为1km，时间分辨率为月，时间序列为2006年1月2006年12月。1.2气象数据：本研究所用的基础数据来源于中国国家气象局，时间为2006年，数据内容为月降水量、月平均气温、月总太阳辐射，以及各气象站点的经度、纬度和海拔高度，共涉及全国726个气象站点（其中太阳辐射数据为112个站点）（图3.4）。对数据进行精度验证，剔除不可替代的错误数据后，计算植被的净初级生产力需要栅格化的气象数据，并从空间上与遥感数据相匹配。利用GIS的插值工具，根据各气象站点的经纬度信息，通过对气象数据进行Kriging插值，获取像元大

23、小与NDVI数据一致、投影相同的气象要素栅格图。图3.4 全国气象站点分布图1.3 植被分类数据： 植被分类图来源于欧盟联合研究中心(The Joint Research Centre:JRC)，原始分类图像由中国科学院遥感应用研究所编译，分类用的遥感数据为2000年的SPOTVGT 1km数据、以及潘耀忠等人基于NOAA/AVHRR和Holdridge的中国土地覆盖分类图，共分为22类。图3.5 长江上游植被分类图表3.2长江上游植被分类情况统计表代码植被类型像元数1落叶针叶林needleleaved deciduoud forest02常绿针叶林needleleaved evergreen

24、 forest2239363常绿阔叶林broadleaved evergreen forest942084落叶阔叶林broadleaved deciduoud forest140805灌丛bush1365126疏林sparse woods113287海边湿地seaside wet lands08高山、亚高山草甸alpine and sub_alpine meadow1913609坡面草地slope grassland1331210平原草地plain grassland2022411荒漠草地desert grassland153612草甸meadow013城市city51214河流river13

25、56815湖泊lake198416沼泽swamp32017冰川glacier332818裸岩bare rocks787219砾石gravels51220荒漠dersert021耕地farmland19321622高山、亚高山草地alpine and sub-alpine plain grass57088合计984896注：植被类型根据2000年SPOT一VGT lkm的遥感植被分类图确定2. 结果与分析2.1长江上游初级生产力估算结果 利用所建立的模型和收集整理后的数据，计算长江上游地区2006年陆地植被总NPP，结果见图6。图3.6 长江上游地区2006年总NPP对长江上游地区2006年陆地

26、植被总NPP进行统计，得到2006年长江上游地区不同植被NPP总值和平均值，见表3.3。表3.3 2006长江上游不同植被NPP统计表代码植被类型像元数平均值总和1落叶针叶林00.000.002常绿针叶林223936355.5679.623常绿阔叶林94208949.8289.484落叶阔叶林14080569.188.015灌丛136512382.6852.246疏林11328456.105.177海边湿地00.000.008高山、亚高山草甸191360435.8983.419坡面草地13312608.708.1010平原草地20224316.196.3911荒漠草地153675.440.12

27、12草甸00.000.0013城市512170.490.0914河流13568401.395.4515湖泊1984353.690.7016沼泽320284.940.0917冰川3328216.120.7218裸岩7872211.851.6719砾石512153.120.0820荒漠00.000.0021耕地193216388.0374.9722高山、亚高山草地57088226.7012.94合计984896435.84429.26注：NPP平均值单位为gC/m2/yr，NPP总量单位为1012gC/yr计算结果表明，2006年长江上游植被净初级生产力为429.261012gC/yr，平均值为4

28、35.84gC/m2/yr。2.2 长江上游初级生产力空间分布通常人们把自然植被分布格局随经纬度和海拔高度的变化称之为“三向地带性”，正是这种“三向地带性”和复杂多变的环境、气候条件共同决定了陆地生态系统NPP的空间格局变化。长江上游东西部气温差异显著。四川盆地东部以及宜宾至宜昌长江干流沿岸一带，气温较高，年均温18以上。成都平原年均温16左右，盆周山区年均温1416。金沙江下游河谷区海拔1300m以下，年平均气温达到20以上，为长江上游气温最高地区。高原大部分年平均气温为08，其中石渠、色达以北及通天河一带，年平均气温低于0。江源区年温在-4左右，系长江上游气温最低地区。长江上游西部高原地区

29、年平均降水较少（200800mm），由东往西减少，源头地区五道梁仅270mm左右。长江上游东部广大地区在太平洋和印度洋暖湿气流控制下，降水丰富，年平均降水量一般在8001500mm，地区分布差异较大，具有由东向西较少的特性。但受地形的影响，在四川盆地周围山地及云贵高原东部分布着几个多雨中心。长江上游幅员辽阔，山地众多，地形复杂，气候多样。上游东部属北亚热带季风和中亚热带湿润季风气候，西北部为山地高原气候，横断山地属于亚热带高原季风气候。不仅具有从南到北热带、亚热带到寒温带的纬向温度变化以及从沿海到内陆的干湿变化，同时亦具有明显的沿海拔高度变化的垂直地带性。长江上游自然植被NPP的分布受水热条件

30、的限制相当明显，陆地生态系统植被NPP空间分布的基本特点是南高北低、东高西低，从西北向东南呈逐级递增趋势。从整个流域来看，青藏高原东部的山麓植被覆盖率较高，植被类型多样，因此净初级生产力也较高，NPP平均值超过了1000 gC/m2/yr；四川盆地以及重庆等丘陵地区主要为农田，NPP在500 gC/m2/yr左右；西北长江源头地区为草地及荒漠区，植被覆盖率极低，且降雨减少，有些地区甚至寸草不生，NPP最小，平均值在100 gC/m2/yr以下。植被净初级生产力受水分及温度的限制比较厉害，温度升高、降水增加都会减少这两个因素的限制作用，因此，NPP的相对增加幅度比较大。这与国际生物圈计划(IBP)的研究结果基本一致。2.3长江上游不同植被类型初级生产力长江上游植被类型复杂多样。按地域水平分异特征，包括不同地带性植被类型，如西部江源草甸草原与草甸、川西暗针叶林、川西南滇北亚热带偏干常绿阔叶林、亚热带偏湿常绿阔叶林、横断山北部温带落叶林以及四川盆地以农业植被为主的植被类型。同时，长江上游许多山地的植被