基于遥感方法的生态资产评估模型20031015.doc

《基于遥感方法的生态资产评估模型20031015.doc》由会员分享，可在线阅读，更多相关《基于遥感方法的生态资产评估模型20031015.doc（36页珍藏版）》请在三一办公上搜索。

1、附件四国家高技术研究发展计划(863计划)子课题阶段验收报告书子课题名称： 区域生态资产遥感测量模型 课题名称：中国西部生态资产定量遥感测量技术体系与应用示范所属专题： 空间信息处理技术所属主题/重大专项： 信息获取与处理技术所属领域： 信 息 技 术 领 域课题责任人 (甲方)： 史 培 军 李 京子课题责任人（乙方）： 潘 耀 忠 起止年限： 2003 年 1 月 至2003 年 9 月北京师范大学资源科学研究所二三 年 九 月目录1 生态系统服务指标体系的建立11.1 生态系统服务及功能11.2 基于遥感手段的区域生态资产评估指标22 基于遥感方法的生态资产评估特点23 生态资产遥感测量

2、的概念模型与技术路线33.1 基于定量遥感测量的生态资产价值评估概念模型33.1.1 生态系统遥感分类43.1.2 生态参数遥感测量43.1.3 单位面积单项生态服务功能价值的定量计算53.1.4 测量结果实地检验53.2 生态资产生态遥感测量的技术路线54 区域生态资产估算的数学模型与方法74.1 植被净初级生产力(NPP)估算74.1.1年度NPP的算法74.1.2APAR的算法74.1.3光能转化率(x, t)的算法84.2生产有机物质的价值94.3吸收固定CO2的价值104.4释放O2的价值104.5N、P、K等营养物质循环的价值114.5.1每年吸收的氮元素（N）的价值114.5.2

3、每年吸收的磷元素（P）的价值114.5.3每年吸收的钾元素（K）的价值124.6涵养水源的价值124.6.1基于母岩性质的贮水量计算134.6.2基于降水转换率的贮水量计算134.7水土保持价值134.7.1减少的土壤侵蚀量134.7.2减少表土损失的价值134.7.3减少土壤养分损失的价值134.7.4减少淤积损失的价值135 数据准备135.1遥感与栅格化的气象数据135.2栅格化的土壤理化性质数据135.3统计数据136 区域生态资产评估结果（内蒙古）136.1基于MODIS数据的生态资产评估结果136.1.1主要的植被参数136.1.2 内蒙古生态资产评估结果136.2基于TM数据的生

4、态资产评估结果137 下一步工作设想137.1 文化与休闲娱乐价值估算137.2 生物多样性价值估算137.3 净化环境污染物的价值估算13参考文献131 生态系统服务指标体系的建立生态系统的功能是指其生境、生物学性质或生态系统过程。生态系统服务功能是指生态系统与生态过程所形成及所维持的人类赖以生存的自然环境条件与效用。生态系统服务功能可以概略地分为两大类，一类是生态系统产品，如为人类提供食物、工业原材料、药品等可以商品化的功能；第二类是支撑与维持人类赖以生存的环境，如生态系统对气候调节、水源涵养、水土保持、土壤肥力的更新与维持、营养物的循环、二氧化碳的固定等难以商品化的功能，从而表现为间接价

5、值。生态服务功能的间接价值虽不表现在国家的核算体制上，但它们的价值可能大大超过直接价值。而且直接价值常常源于间接价值。1.1 生态系统服务及功能Costanza等（1997）把生态系统服务分为17大类（表1-1），其中只包括可再生的服务，不包括不可再生的燃料和矿物质以及大气。表1-1生态系统服务及功能序号生态系统服务生态系统功能举例1气体调节大气化学成份调节CO2/O2平衡，O3防紫外线，SOx水平2气候调节全球温度、降水及其它由生物媒介的全球及地区性气候调节温室气体调节，影响云形成的DMS产物3干扰调节生态系统反应对环境波动的容量、衰减和综合风暴防止，洪水控制，干旱恢复等生境对主要受植被结构

6、控制的环境变化的反应4水调节水文流的调节为农业、工业和运输提供用水5供水水的贮存和保持向集水区、水库和含水岩层供水6控制侵蚀和保持沉积物生态系统内的土壤保持防止土壤被风、水侵蚀，把淤泥保存在湖泊和湿地中7土壤形成土壤形成过程岩石风化和有机质积累8养分循环养分的贮存、内循环和获取固氮，N、P和其它元素及养分循环9废物处理易流失养分的再获取、过多或外来养分、化合物的去除或降解废物处理，污染控制，解除毒性10传粉有花植物配子的运动提供传粉者以便植物种群繁殖11生物控制生物种群的营养动力学控制关键捕食者控制猎物种群，高级捕食者使食草动物减少12避难所为定居和迁徙种群提供生境育雏地，迁徙动物栖息地，当地

7、收获物种栖息地或越冬场所13食品生产总初级生产中可用为食物的部分通过渔、猎、采集和农耕收获的鱼、鸟兽、作物、坚果、水果等14原材料总初级生产中可用为原材料的部分木材、燃料和饲料产品15基因资源独一无二的生物和产品的资源医药、材料科学产品，用于农作物抗病和抗植物感染的基因，家养物种（宠物和植物栽培品种）16休闲娱乐提供休闲游乐活动机会生态旅游、钓鱼运动及其它户外游乐活动17文化提供非商业性用途的机会生态系统的美学、艺术、教育、精神及科学价值1.2 基于遥感手段的区域生态资产评估指标表1列出的生态系统服务和功能并不一定一一对应。在有些情况下，一种生态功能可提供两种或多种服务，而且许多生态系统功能是

8、相互依赖的。从宏观生态学角度，根据基于遥感手段的区域生态资产计算的特点，考虑数据获取的可能性和可靠性，以及我国在这方面的研究情况（欧阳志云，1999；黄兴文，1999；陈仲新，2000），最终确定的基于遥感手段的区域生态资产评估指标如表1-2所示：表1-2基于遥感手段的生态资产评估指标序号生态系统服务生态系统功能1太阳能的固定与食品生产生产有机物质2促进营养元素循环营养物质循环与贮存3水调节与供水涵养水源4控制侵蚀和保持沉积物保持土壤5气体调节吸纳二氧化碳6气体调节释放氧气2 基于遥感方法的生态资产评估特点（1）遥感模型用于生态资产评估可以克服传统方法“以点代面”的缺点，减少了统计模型用地面观

9、测值插值所引起的误差，提高了评估精度。（2）遥感数据具有不同的波谱、时空分辨率，通过多源数据的融合处理，可以较精确的进行各种生态系统类型的划分，并且可以反演出一些生态与非生态的模型参数（如生物量模型、NPP模型，蒸散量模型等）。（3）相对于传统的数据采集途径来说，遥感数据成本较低，而且获取容易，速度较快。（4）基于遥感方法的生态资产评估模型，可以较容易地搭建在GIS平台上，通过和地面统计数据的链接，利用计算机技术可以实现区域生态资产自动化评估。（5）受数据特点和自动化要求的限制，所选择的方法与模型会使一部分生态服务被忽略，但从工程性的角度来看，这样可以提高系统的可行性。对于被忽略的生态服务，如

10、果通过调查发现这部分价值不应该完全被忽略，则可以通过其它的途径进行计算，例如将某种被忽略的生态价值纳入到其它价值中进行计算，或者根据经验判断这部分价值在总价值中所占比重，这样可以通过计算其它价值来进行估算。在区域生态资产评估中，有些服务，譬如游憩价值和生物多样性的价值，比较难于定量，通常的方法都需要有人的互动环节，但在工程中，这显然不太适用。3 生态资产遥感测量的概念模型与技术路线3.1 基于定量遥感测量的生态资产价值评估概念模型一定区域内的生态资产的总量是一个随时间动态变化的量值，它是区域内所有生态系统类型提供的所有服务功能及其自然资源价值的总和，并随着区域内所含有的生态系统的类型、面积、质

11、量的变化而变化。Costanza（1997）等人在进行全球生态资产测量时，将全球生态系统划分为17个类型的研究方法，为基于定量遥感测量的生态资产价值评估模型的建立提供了一种可借鉴的方法。一定区域内的生态资产价值总量（V）可以表示为： （1）其中：c=1,2,n，表示生态系统的类型；Vc表示第C类生态系统生态资产价值。 （2）其中：i=1,2,n，表示第C生态系统的第i种生态服务功能；Vci表示第C生态系统的第i种生态服务功能类型的单位面积价值；j=1,2,m表示一定区域内Vci在空间上分布的斑块数，Sij表示各个斑块的面积大小，Rij表示Vci在不同斑块的调整系数，它是由生态系统的质量状况决定

12、的。 （3）其中：A1,A2,An是表征生态系统质量状况的生态参数，不同质量状况的同一生态系统类型的Vci是不同的，其差别由Rij来调整。根据区域生态资产价值计算的要求，基于定量遥感测量的价值评估模型主要包括以下四部分内容：生态系统遥感分类、生态参数遥感测量、单位面积（或单位质量）生态服务功能价值定量计算和测量结果经度检验，生态价值评估模型的技术流程如图3-1所示。多尺度遥感数据生态系统分类体系的确定生态系统类型测量生态服务功能的定义生态参数测量不同类型的空间分布单位面积价值量计算调整系数计算生态学/环境经济学估算方法/模型生态资产价值量图形输出统计输出图3-1 基于遥感测量的生态资产价值评估

13、模型流程图精度检验3.1.1 生态系统遥感分类不同的生态系统类型包含的自然资源价值和具有的生态服务功能是完全不同的，因而基于遥感的生态资产测量的第一步，就是根据不同的空间尺度确定适合于生态资产测量的生态系统分类体系，然后利用多尺度遥感技术获取这些类型的空间分布状况，在这一领域，遥感技术已经积累了相当丰富的经验。3.1.2 生态参数遥感测量不同空间分布的同一生态系统类型的质量状况是不同的，同时也是随时间变化的，因而其所具有的生态资产价值也是变化的，这种变化是由生态系统的质量状况来决定的。生态系统的质量状况可以由一系列的生态参数来表达，这些参数随时间、空间的变化的全覆盖数据只有通过遥感测量才能够得

14、到。根据前人的研究经验，本文认为：选取地上总生物量、第一性净生产力（反映有机质生产的两个重要指标），以及植被覆盖度（反映生态系统随时间、空间变化的重要指标）三个参数来表达生态系统的质量状况，基本符合生态资产遥感定量测量的要求。3.1.3 单位面积单项生态服务功能价值的定量计算包括两部分：自然资源单位面积价值和生态服务功能单位面积价值。同一生态系统同时具有多种生态服务功能，每一种生态服务功能的价值及其计算方法也是不同的，为使各种生态资产计算方法的结果具有可比性，建立单项服务功能的单位面积价值的标准计算方法和估算模型是必不可少的，这需要生态学和环境经济学为基础的大量的实地测量的基础性研究结果为依据

15、的，目前这方面已经取得了一定的成果。3.1.4 测量结果实地检验根据分层统计抽样的原理，根据不同测量尺度的要求，与遥感测量相匹配，按一定的标准和规范设计实地抽样测量样地和样点，按一定的指标体系和标准的调查测量规程，建立基于实地抽样测量的生态资产测量结果验证方法，并依此作为验证遥感测量结果精度的依据。3.2 生态资产生态遥感测量的技术路线生态资产生态遥感测量的技术体系如图3-2所示。模型的数据来源主要有三条途径：气象数据、遥感数据、地面观测与统计资料。遥感影像数据一方面用于生态系统类型的划分和面积估算，另一方面用于生态系统生物量、净初级生产力、植被覆盖度等生态参数的反演。生态系统有机物质生产、营

16、养物质循环、CO2吸收、O2释放这四项生态价值的评估是基于生态系统生物量和净初级生产力来计算的，即根据单位质量生物量所产生的价值进行估算。生态系统涵养水源的效益分两种情况估算。当下垫面为土壤时（如林地、灌丛、草地等），根据不同植被覆盖下的土壤储水能力来计算涵养水源价值；当下垫面为水时（如水域、水稻田、沼泽等），则根据降水转换率（降水贮存量占总降水量的百分比）来计算。生态系统保护土壤主要通过减少表土损失量，保护土壤肥力，减轻泥沙淤积灾害，减少风沙等灾害的4个相互联系的生态过程来实现其经济价值。在估算过程中，首先采用无植被覆盖的潜在土壤侵蚀量来估算各生态系统减少的土壤侵蚀量，在此基础之上再评价它们

17、对表土损失、肥力损失和减轻泥沙淤积灾害三方面的价值。所获得的各项量化的生态效益，最后通过影子工程法、市场价值法、机会成本法、替代价值法等转换成生态系统服务价值。市场价值法影子工程法机会成本法替代价值法碳税法工业制氧法能量固定替换法市场价值法影子工程法直接使用价值法遥感影像数据气象数据地面观测与统计资料生态系统类型及其面积各生态系统生物量各生态系统净初级生产力生产有机物质N、P、K等营养物质循环吸收固定CO2释放O2降水量涵养水源各类生态系统的土壤侵蚀量减少表土损失量保护土壤肥力：有机质、N、P、K减轻泥沙淤积水土保持生态系统服务价值图3-2生态资产生态遥感测量的技术路线4 区域生态资产估算的数

18、学模型与方法4.1 植被净初级生产力(NPP)估算4.1.1年度NPP的算法 （1）NPP(x)：象元x处全年净初级生产力（gC）；NPP(x, t)：t月份在象元x处的单位面积净初级生产力（gCm-2mon-1）；Sp：象元面积（m2）。 （2）APAR(x, t)：t月份在象元x处植被吸收的光合有效辐射（MJm-2）；(x, t)：t月份在象元x处的光能转化率（gCMJ-1）。4.1.2APAR的算法（3）SOL(x, t)：t月份在象元x处的太阳总辐射量（MJm-2）；FPAR(x, t)：t月份在象元x处植被层对入射的光合有效辐射的吸收比例；r：植被所能利用的太阳辐射占太阳总辐射的比例

19、，在此取值0.5。（4）（5）SR(x, t)：比值植被指数，SRmin取值为1.08，SRmax的大小与植被类型有关（各植被类型中的最大值），取值范围在4.146.17之间；NDVI(x, t)：归一化植被指数。4.1.3光能转化率(x, t)的算法（6）Te1(x, t)：t月份在象元x处的温度胁迫系数，反映在低温和高温时植物内在的生化作用对光合的限制而降低净第一性生产力（Potter, et al., 1993; Field, et al., 1995）Te2(x, t)：t月份在象元x处的温度胁迫系数，表示环境温度从最适温度Topt(x)向高温和低温变化时植物的光合转化率逐渐变小的趋势

20、（Potter, et al., 1993; Field, et al., 1995）We1(x, t)：t月份在象元x处的水分胁迫系数，反映水分条件的影响emax：理想条件下的最大光能转化率，在此取0.389 gCMJ-1（Potter, et al., 1993）4.1.3.1Te1(x, t) 的算法（7）（当某一月的平均气温MeanT(x, t)10时，Te1(x, t) = 0） Topt(x)：某一区域一年内NDVI达最大值时月份的平均气温（）。4.1.3.2 Te2(x, t)的算法（8）当T(x, t) Topt(x, t) + 10；或者T(x, t) Topt(x, t)

21、13时，用公式（9）计算：（9）T(x, t)：某月的平均温度（）。4.1.3.3We1(x, t) 的算法 (10)EET(x, t)：t月份在象元x处的区域估计蒸散量（mm）；PET(x, t)：t月份在象元x处的区域可能蒸散量(mm)；当该月降水量PPT(x, t)PET(x, t)时，EET(x, t)PET(x, t)，即We1(x, t)=1；当该月降水量PPT(x, t)PET(x, t)时，用以下公式计算EET(x, t)和PET(x, t)： (11)PPT(x, t)：t月份在象元x处的月降水量（mm）；Rn(x, t)：t月份在象元x处的月净辐射量(mmmon-1)； (

22、12)Ep(x, t)：局地可能蒸散量(mm) (13) (14) (15) (16)4.2生产有机物质的价值净初级生产力是反映有机物质生产的一个重要指标，它是指植物在某一时间段（如一年）所生产的有机物质量。Vn(x): 象元x处每年生产的有机物质价值（元）；NPP(x): 象元x处每年生产的有机物质（gC）；T(x)：象元x处有机物质的单位质量价值（元g-1C）；Vn：区域每年生产有机物质的价值（元）。4.3吸收固定CO2的价值在评估生态系统对CO2的吸收与固定作用时，以生态系统有机物质生产为基础，根据光合作用和呼吸作用的反应方程式推算，每形成1g干物质，需要1.62g CO2。Vr(x)：

23、象元x处每年吸收CO2的价值（元）；NPP(x): 象元x处每年生产的有机物质（gC）；R：碳税法中CO2的单位质量价值（元gCO2）；Vr：区域每年吸收CO2的价值（元）。4.4释放O2的价值在评估生态系统释放O2的价值时，以生态系统有机物质生产为基础，根据光合作用和呼吸作用的反应方程式推算，每形成1g干物质，可以释放1.2gO2。V(x)：象元x处每年释放O2的价值（元）NPP(x): 象元x处每年生产的有机物质（gC）R：工业制氧价格（元g-1O2）。V：区域释放O2的价值（元）。4.5N、P、K等营养物质循环的价值在评估生态系统营养物质循环的价值时，以生态系统的生物量和净初级生产力为基

24、础，根据各生态系统中N、P、K的质量分配率来计算(表5-1)。Va：区域每年吸收的营养物质价值（元）；Vna: 区域每年吸收的氮元素价值（元）；Vpa: 区域每年吸收的磷元素价值（元）；Vka: 区域每年吸收的钾元素价值（元）。4.5.1每年吸收的氮元素（N）的价值Vna(x)：象元x处每年吸收的氮元素价值（元）；NPP(x): 象元x处每年生产的有机物质（gC）；r1：不同生态系统中氮元素在有机物质中的分配率（gN/gC）；r2：纯N折算成化肥的比例，为79/14；M：氮肥的平均价格（元g-1）；Vna: 区域每年吸收的氮元素价值（元）。4.5.2每年吸收的磷元素（P）的价值Vpa(x)：象

25、元x处每年吸收的磷元素价值（元）；NPP(x): 象元x处每年生产的有机物质（gC）；r1：不同生态系统中磷元素在有机物质中的分配率（gN/gC）；r2：纯P折算成化肥的比例，为506/62；M：磷肥的平均价格（元g-1）；Vpa: 区域每年吸收的磷元素价值（元）。4.5.3每年吸收的钾元素（K）的价值Vka(x)：象元x处每年吸收的钾元素价值（元）；NPP(x): 象元x处每年生产的有机物质（gC）；r1：不同生态系统中钾元素在有机物质中的分配率（gK/gC）；r2：纯K折算成化肥的比例，为174/78；M：钾肥的平均价格（元g-1）；Vka: 区域每年吸收的钾元素价值（元）。4.6涵养水源

26、的价值V(x)：象元x处每年涵养水源的价值（元）；Q(x)：象元x处每年单位面积的降水贮水量(m3/m2)；p：每建设1m3库容的成本花费；Sp：象元面积（m2）；V：区域涵养水源的价值（元）。下垫面为土壤时（如林地、灌丛、草地等）采用4.6.1的方法计算涵养水源价值；下垫面为水时（如水域、水稻田、沼泽等）采用4.6.2的方法来计算。4.6.1基于母岩性质的贮水量计算Q(x)：象元x处每年单位面积的降水贮水量(m3/m2)；H(x)：象元x处的土壤深度（m）；P(x)：象元x处的土壤粗孔隙率。4.6.2基于降水转换率的贮水量计算Q(x)：象元x处每年单位面积的降水贮水量(m3/m2)；J(x)

27、：象元x处的年降水量（m）；R(x)：象元x处的降水贮存量占总降水量的百分比。4.7水土保持价值生态系统保护土壤主要通过减少表土损失量，保护土壤肥力，减轻泥沙淤积灾害，减少风沙等灾害的4个相互联系的生态过程来实现其经济价值。在估算过程中，首先采用无林地的土壤侵蚀量来估算各生态系统减少的土壤侵蚀量，然后再评价它们对表土损失、肥力损失和减轻泥沙淤积灾害三方面的价值。V(x)：在象元x处每年的水土保持价值（元）；V1(x)：在象元x处每年减少表土损失的价值（元）；V2(x)：在象元x处每年减少有机质损失的价值（元）；V3(x)：在象元x处每年减少氮肥损失的价值（元）；V4(x)：在象元x处每年减少磷

28、肥损失的价值（元）；V5(x)：在象元x处每年减少钾肥损失的价值（元）；V6(x)：在象元x处每年减少泥沙淤积损失的价值（元）。V：区域每年水土保持价值（元）；V1：区域每年减少表土损失的价值（元）；V2：区域每年减少有机质损失的价值（元）；V3：区域每年减少氮肥损失的价值（元）；V4：区域每年减少磷肥损失的价值（元）；V5：区域每年减少钾肥损失的价值（元）；V6：区域每年减少泥沙淤积损失的价值（元）。4.7.1减少的土壤侵蚀量Er(x)：象元x处减少的土壤侵蚀量（m3）；En(x)：象元x处在裸露条件下的潜在土壤侵蚀深度(m)；Ef(x)：象元x处在有树木、灌丛、草、农作物等植被覆盖情况下的

29、侵蚀深度(m)；Sp：象元面积（m2）；Fc(x)：象元x处的植被覆盖度。NDVI(x)：象元x处的归一化植被指数。NDVImin和NDVImax分别为植物生长季NDVI的最小值和最大值（西北地区分别取值0.005和0.92）。4.7.2减少表土损失的价值V1(x)：在象元x处每年减少表土损失的价值（元）；Er(x)：象元x处减少的土壤侵蚀量（m3）；D(x)：象元x处的土壤深度（m）；T(x)：象元x处的土壤生产机会成本（元/m2）（表52）。4.7.3减少土壤养分损失的价值减少养分损失包括减少有机质损失、减少N损失、减少P损失和减少K损失。4.7.3.1减少有机质损失的价值V2(x)：在象

30、元x处每年减少有机质损失的价值（元）；Er(x)：象元x处减少的土壤侵蚀量（m3）；De(x)：象元x处土壤的密度（t/m3）；Co(x)：象元x处土壤的有机质含量(kg/t)；R：薪柴转化为土壤有机质的比例，一般为2；C：薪柴的平均价格（元/kg）。4.7.3.2减少N损失的价值V3(x)：在象元x处每年减少氮肥损失的价值（元）；Er(x)：象元x处减少的土壤侵蚀量（m3）；De(x)：象元x处土壤的密度（t/m3）；Cn(x)：象元x处土壤的N含量(kg/t)；r：纯N折算成化肥的比例，为79/14；C：氮肥的平均价格（元/kg）。4.7.3.3减少P损失的价值V4(x)：在象元x处每年减

31、少磷肥损失的价值（元）；Er(x)：象元x处减少的土壤侵蚀量（m3）；De(x)：象元x处土壤的密度（t/m3）；Cp(x)：象元x处土壤的P含量(kg/t)；r：纯P折算成化肥的比例，为506/62；C：磷肥的平均价格（元/kg）。4.7.3.4减少K损失的价值V5(x)：在象元x处每年减少K肥损失的价值（元）；Er(x)：象元x处减少的土壤侵蚀量（m3）；De(x)：象元x处土壤的密度（t/m3）；Ck(x)：象元x处土壤的K含量(kg/t)；r：纯K折算成化肥的比例，为174/78；C：钾肥的平均价格（元/kg）。4.7.4减少淤积损失的价值土壤侵蚀造成的泥沙有24%淤积于水库、江河、湖

32、泊，这部分泥沙直接造成了水库、江河、湖泊蓄水量的下降，在一定程度上增加了干旱、洪涝灾害发生的机会，因此可根据蓄水成本计算损失价值。V6(x)：在象元x处每年减少淤积损失的价值（元）；Er(x)：象元x处减少的土壤侵蚀量（m3）；R：淤积于水库、江河、湖泊中的泥沙占土壤总侵蚀量的比例，全国平均为24%。C：平均库容工程费（元/m3）。我国为0.67元/m3。5 数据准备5.1遥感与栅格化的气象数据1. 栅格化的全年每月的AVHRR/NDVI数据；2. 栅格化的全年每月的降水数据（mm）；3. 栅格化的全年每月的平均温度（）；4. 栅格化的全年每月的平均太阳总辐射数据（MJm-2）；5. 栅格化的

33、植被分类图。5.2栅格化的土壤理化性质数据1. 栅格化的土壤深度分布图（m）；2. 栅格化的土壤粗孔隙率分布图；3. 栅格化的裸露条件下的潜在土壤侵蚀深度分布图(m)；4. 栅格化的有植被覆盖情况下的土壤侵蚀深度分布图(m)；5. 栅格化的土壤的有机质含量分布图(kg/t)；6. 栅格化的土壤N含量分布图(kg/t)；7. 栅格化的土壤P含量分布图(kg/t)；8. 栅格化的土壤K含量分布图(kg/t)；9. 栅格化的土壤密度分布图（t/m3），土壤的平均密度为1.28 t/m3。10. 栅格化的土壤粘粒所占百分比；11. 栅格化的土壤砂粒所占百分比。5.3统计数据1. 有机物质的单位质量价值

34、（元/gC）；根据标煤的价值来替代，（标煤为354元t-1），其取值为2.3718104元/gC；2. 碳税法中CO2的单位质量价值（元/gCO2）。借用瑞典碳税率0.15美元/kg(碳)，换算成吸收CO2的税率为4.094105美元/g（CO2），如果按8.2元/美元的汇率计算，即3.36104元/g（CO2）；3. 工业制氧价格（元/gO2），工业制氧的现价为4104元/g(O2)；4. 各生态系统营养物质分配率（表51）；5. 氮肥的平均价格（元g-1），400元/t；6. 磷肥的平均价格（元g-1），350元/t；7. 钾肥的平均价格（元g-1），350元/t；8. 每建设1m3库容的

35、成本花费(元/m3)，中国是0.67元/m3；9. 下垫面为水时，降水贮存量占总降水量的百分比（由统计资料确定）；10. 土壤生产机会成本（元/m2）。其取值是根据不同的生态系统类型来确定（表52），如在19851990年期间林业、牧业每年生产的平均收益分别为263.58元/hm2与245.50元/hm2；11. 薪柴的平均价格（元/kg）；12. 淤积于水库、江河、湖泊中的泥沙占土壤总侵蚀量的比例。全国平均为24%。表51各生态系统营养物质分配率（g/g）项目林地灌丛草地耕地建筑难利用地水域氮0.0041800.0132940.0132890.0132880.0132030.0132730.

36、004204磷0.0008900.0000920.0000930.0000900.0000870.0000910.000901钾0.0018100.0089040.0089080.0089150.0088740.0089090.001802各生态系统营养物质质量分配率是根据中国生物多样性国情研究报告整理出来的表52各生态系统的土壤生产机会成本（元/ahm2）（基于内蒙古1999年的统计数据）项目林地灌丛草地耕地建筑难利用地水域机会成本160.3160.3335.02324.50006 区域生态资产评估结果（内蒙古）根据生态资产评估的技术路线和模型，以内蒙古为例，对其生态资产进行了评估。本项目研

37、究所利用的遥感数据为MODIS数据和TM数据。MODIS数据的图像空间分辨率为500m500m，时间分辨率为月，时间序列为2002年1月至2002年12月；TM数据的图像空间分辨率为30m30m，共包含内蒙古自治区内的86景图像（内蒙古东北的部分TM图像缺失）。气象数据来自内蒙古气象局，为2002年内蒙古108个气象站点每月的降水、平均温度及太阳总辐射数据。所利用的土地覆被图为本课题组基于TM影像的分类图（详见“生态资产遥感测量的地表分类”）、以及中国科学院遥感应用研究所2002年生产的全国土地覆被图（下载网址：http:/www.gvm.jrc.it/glc2000）。土壤类型图来自中国科学

38、院地理研究所资源与环境信息系统国家重点实验室生产的1:400万的失量图（中国科学院南京土壤研究所，1978），土壤理化性质数据均是在土壤分类图的基础上参考相关的文献资料更新得到（熊毅，李庆逵，1987；沈善敏，1998）。6.1基于MODIS数据的生态资产评估结果6.1.1主要的植被参数1.栅格化的全年净初级生产力（gC/m2）分布图；图61 内蒙古NPP分布图（gC/m2）2. 栅格化的植被分类图；图62 内蒙古植被分类图3. 栅格化的植被覆盖度分布图；图63 内蒙古植被覆盖度分布图6.1.2 内蒙古生态资产评估结果1.栅格化的每年生产的有机物质价值（元/hm2）分布图；图64 每年生产的有

39、机物质的价值（元/hm2）2.栅格化的每年吸收CO2的价值（元/hm2）分布图；图65 每年吸收的CO2的价值（元/hm2）3.栅格化的每年释放O2的价值（元/hm2）分布图；图66 每年释放O2的价值（元/hm2）4. 栅格化的每年吸收的氮元素的价值（元/hm2）分布图；图67 每年吸收氮元素的价值（元/hm2）5. 栅格化的每年吸收的磷元素的价值（元/hm2）分布图；图68 每年吸收的磷元素的价值（元/hm2）6. 栅格化的每年吸收的钾元素的价值（元/hm2）分布图；图69 每年吸收的钾元素的价值（元/hm2）7. 栅格化的每年涵养水源的价值（元/hm2）分布图；图610 每年涵养水源的价

40、值（元/hm2）8. 栅格化的每年的水土保持价值（元/hm2）分布图；图611 水土保持的价值（元/hm2）9. 栅格化的每年减少表土损失的价值（元/hm2）分布图； 图612 每年减少表土损失的价值（元/hm2）10. 栅格化的每年减少有机质损失的价值（元/hm2）分布图；图613 每年减少有机质损失的价值（元/hm2）11. 栅格化的每年减少氮肥损失的价值（元/hm2）分布图；图614 每年减少氮肥损失的价值（元/hm2）12. 栅格化的每年减少磷肥损失的价值（元/hm2）分布图；图615 每年减少磷肥损失的价值（元/hm2）13. 栅格化的每年减少钾肥损失的价值（元/hm2）分布图； 图

41、616 每年减少钾肥损失的价值 （元/hm2） 14. 栅格化的每年减少泥沙淤积损失的价值（元/hm2）分布图；图617 每年减少泥沙淤积损失的价值（元/hm2）15. 栅格化的总生态价值（元/hm2）分布图；图618 总生态价值分布图（元/hm2）16. 七种生态系统类型各项生态服务功能的基本情况表61七种生态系统类型各项生态服务功能统计表项目森林灌木草地耕地建筑用地难利用地水体合计面积(106 hm2)13.093.7955.8813.710.0001527.582.58116.62NPP (106 t/a)61.064.4993.5333.270.0013 5.359.52206.95C

42、O2年固定量(106 t/a)98.97.26151.553.90.002 8.6615.4335.6O2 释放量(106 t)73.285.38112.2539.930.0002 6.4511.43248.72氮固定量(06 t/a)1.440.347.012.49 0 0.400.2311.92磷固定量(106 t/a)0.440.0030.070.02 0 0.0040.070.61钾固定量(106 t/a) 0.220.081.670.59 0 0.090.032.6917. 七种生态系统类型各项生态价值分布情况表62七种生态系统类型各项生态价值统计表项目森林灌木草地耕地建筑用地难利用地水体合计生产有机物(109 元)14.481.0622.187.890.00003 1.272.2649.14固定CO2(109 元)33.242.4450.9118.11 0 2.915.18112.79释放O2 (109 元)29.312.1544.9015.97 0 2.584.5799.48营养物质循环 (109 元)0.810.173.411.21 0 0.190.135.92涵养水源(109 元)11.801.4531.897.85 0 2.972.1858.14水土保持(109 元)20.202.5043.20 0 0 0065.90总计(109 元)109