株洲市森林生态系统服务功能价值评估

株洲市森林生态系统服务功能价值评估钟智群，陈少波，苏勇，邓文韬，廖常乐，王婧菲（株洲市林业局 湖南 株洲 412007）摘 要：采用株洲市“十二五”森林资源二类调查资料，根据国家林业局发布的森林生态系统服务功能评估规范)(LYT17212008)标准，对株洲市8种主要森林生态系统服务功能价值进行了评估。结果表明：8种森林生态系统服务功能总价值为787.57亿元/a，其中提供林产品67.72亿元/a，涵养水源416.52亿元/a，固土保肥43.83亿元/a，固碳释氧95.66亿元/a，积累营养物质5.35亿元/a，净化大气环境25.52亿元/a，生物多样性保护128.65亿元/a，森林游憩4.31亿元/a。关键词：森林生态；服务功能；价值评估中图分类号： 文献标识码：B 文章编号：森林是地球陆地最大的生态系统，其生态服务功能已成为近几十年来研究的热点15。森林生态系统服务功能是森林生态系统与生态过程所形成及维持的人类赖以生存的自然环境条件与效用。它主要包括森林在涵养水源、保育土壤、固氮释氧、积累营养物质、净化大气环境、森林保护、生物多样性保护和森林游憩等方面提供的生态服务功能13。对森林生态系统服务功能做出正确的价值评估，能有效指导森林资源的合理利用和可持续发展，为维护森林生态系统的健康发展提供有力保障。国外对生态系统服务价值的研究已有几十年历史，Constanza等人在世界上率先对全球生物圈生态系统服务功能价值进行了估算，他们估计全球生态系统的服务每年总价值为1654万亿美元，平均为33万亿美元8; Peter估计全球生态系统服务功能为3万亿美元/a9。我国自上世纪80年代开始森林生态系统服务评价工作，大多数研究是借鉴国外的方法和评估指标体系12。蒋延玲等(1999)、陈仲新等(2000)分别对我国森林生态系统和生态系统服务总价值进行了评估，其经济价值分别为117401亿美元/a和7783448亿元人民币a。韩素芸、田大伦、闫文德(2009)等对湖南省主要森林类型生态系统服务功能进行了研究；黄怀雄、赵红艳（2010）对长株潭地区森林生态系统固碳释氧进行了研究。而对株洲市森林生态系统综合服务功能的研究未见报道。本研究主要采用国家林业局2008年发布实施的林业行业标准森林生态系统服务功能评价规范( LY/T17212008)中的评估指标体系和评估公式，利用株洲市“十二五”森林资源二类调查成果资料，对株洲市森林生态系统服务功能进行价值评估，旨在为株洲市生态环境保护政策的制定，以及株洲市林业可持续发展提供参考依据，对合理经营森林资源有着十分重要的意义。1 研究区域概况1.1 自然地理概况株洲市位于北纬260305280107，东经11257301140715之间，地处罗霄山脉西麓，南岭山脉至江汉平原的倾斜地段上，市域总体地势东南高、西北低。北中部地形岭谷相间，盆地呈带状展布；东南部均为山地，山峦迭障，地势雄伟。水域637.27km2，占市域总面积的5.66%；平原1843.25 km2，占16.37%；低岗地1449.86 km2，占12.87%；高岗地738.74 km2，占6.56%；丘陵1916.61 km2，占17.02%；山地4676.47 km2，占41.52%。山地主要集中于市域东南部，岗地以市域中北部居多，平原沿湘江两岸分布。区内属亚热带季风性湿润气候，四季分明，雨量充沛、光热充足，风向冬季多西北风，夏季多正南风，年降水量1400-1700mm，年均气温17.2，年日照时数1677.1小时，无霜期286天以上。土壤属中亚热带常绿阔叶林红壤、黄壤地带，在低山、丘陵和低丘岗地主要分布着以板页岩、花岗岩和第四纪红色粘土层发育的红壤和山地红黄壤，部分低丘岗地分布着紫色砂页岩发育的紫色土，其它土壤类型较少。1.2 森林资源概况株洲市森林植被主要属中亚热带常绿阔叶林带，森林资源丰富，植物种类繁多，群落交错，分布混杂。自然分布和引种栽培的植物约有106科、296属、884种。其中珍稀乡土树种约有40余种。截止2013年底，全市林地总面积738376.7ha。其中有林地面积554722.3ha，占林地面积的75.13，森林覆盖率60.41%；疏林地面积6441.9ha，占林地面积的0.87%；灌木林面积129646.4ha，占17.56%；宜林地面积13522ha，占1.83%。全市森林活立木蓄积量2101.75万m3。其中杉木1036.47万m3，阔叶树728.18万m3，马尾松252.05万m3，湿地松66.44万m3，其它林木18.61万m3。全市经济林面积81252.8ha,楠竹林面积72974.9ha。从森林的年龄构成来看，幼龄林和中龄林比例大，分别占29.36和35.35； 近熟林、成熟林和过熟林分别占20.90、12.46和1.95。株洲市生态公益林面积277095.5ha，商品林面积461281.2ha。2 研究方法2.1 评价范围、指标及数据来源本研究以株洲市辖区范围自然生态系统为研究单元，对其中的主要森林生态系统（包括杉木林、阔叶林、马尾松林、国外松林、竹林、经济林、不含经济林的特灌林以及其他林分）服务功能产生的物质量及其经济价值进行研究，包括森林生态系统提供林产品以及涵养水源、土壤保持、营养物质循环、净化空气、固碳释氧和森林游憩等八项生态系统服务功能。数据来源为株洲市“十二五”森林资源二类调查成果资料以及有关专家的研究成果。2.2 采用的计算方法和公式 本研究对森林生态系统提供的实物产品采取产量和市场价值法进行评估，对其它生态服务功能的评价采用国家林业局发布的森林生态系统服务功能评估规范（LY/T17212008)中的公式计算。2.2.1 森林生态系统提供林产品服务功能指森林生态系统提供木材、食品、药材及其他工业原材料的功能。本研究采用市场价值法来评估其价值。具体计算方法是根据研究区域内生产各类林产品的单位面积产量乘以各种林分面积,得出年净生产量后,再乘以它们的市场价格,即可得出它们价值。对于木材来说,在得出各林分的年净生长量后, 以此净生长量直接与相应林分活立木的价格相乘即可得到活立木的价值。根据株洲市森林提供林产品的实际情况，我们选择木材、楠竹、油茶和松脂等主要林产品作为评估对象。2.2.1.1 森林提供木材产品功能的评价方法计算公式：U1=Q1P1R1/S1 （1）式中：U1：森林生态系统提供木材产品功能价值；Q1：森林净初级生产量（t/a）；P1：木材平均价格（元/m3）；R1：森林综合出材率（%）；S1：木材平均比重（0.48t/m3）。2.2.1.2 森林提供竹材产品功能的评价方法计算公式：U竹=A竹P竹R竹 （2）式中：U竹为竹林提供竹材价值（元/a）；A竹为竹林面积（ha）；P竹为竹材市场价格（元/根）；R竹为单位面积竹材产量。2.2.1.3 森林提供油茶产品功能价值评价方法计算公式：U油= R油P油 （3）式中：U油为森林生产油茶的价值（元/a）；R油为油茶年平均产量（t/a）；P油为油茶市场价格。2.2.1.4 森林提供松脂产品功能价值评价方法计算公式：U脂=A脂P脂R脂 （4）式中：U脂为松脂价值（元/a）；A脂为可采松树面积（ha）；R脂为单位面积可采松脂数量(t/ha)； P脂为松脂市场价格（元/t）。2.2.2 森林生态系统涵养水源功能价值评价方法森林对降水的截留、吸收和贮存，将地表水转为地表径流或地下水的作用13。主要功能表现在增加可利用水资源、净化水质和调节径流三个方面。本研究将涵养水源分解为调节水量和净化水质2个指标进行计算。1）调节水量价值：Um=10CvA(P-E-C) （5）式中：Um为林分年调节水量价值（元/a），Cv为水库建设单位库容投资价格（元/m 3),A为林分面积，P为降水量（mm/a），E为林分蒸散量（mm/a），C为地表径流量（mm/a）。2）净化水质价值：Uwq=10KA(P-E-C) （6）式中：Uwq为林分年净化水质价值（元/a），K为水的净化费用（元/t）。2.2.3 森林生态系统保育土壤功能价值评价方法森林中活地被物和凋落物截留降水，降低水滴对表土的冲击和地表径流的侵蚀作用；同时林木根系固持土壤，防止土壤崩塌泻溜，减少土壤肥力损失及改善土壤结构的功能13。保育土壤分为固土价值和保肥价值2个指标。1）固土价值：利用有林地与无林地的土壤侵蚀模数差值和泥沙输移比值计算土壤的流失量，再通过清淤工程费用计算固土减淤功能价值。计算公式：Usm=ACs(X2-X1)/p （7）式中：Usm为林分年固土价值（元/a)，Cs为挖取和运输单位体积土方所需费用（元/m 3),X2为无林地土壤侵蚀模数（t/haa）,X1为林地土壤侵蚀模数（t/haa），p为林地土壤密度（t/m 3)。2）保肥价值：按照有林地比无林地每年减少土壤侵蚀量中N、P、K的数量乘以N、P、K肥市场售价折算后带来的间接经济效益计算。计算公式：Uf=A(X2-X1)(NC1/R1+PsC1/R2+KsC2/R3+MsC3) （8） 式中：Uf为林分年保肥价值（元/a),N为林分土壤平均含氮量（%），C1为磷酸二铵化肥价格（元/t），R1为磷酸二铵化肥含氮量（%），Ps为林分土壤平均含磷量（%），R2为磷酸二铵化肥含磷量（%），Ks为林分土壤平均含钾量（%），C2为氯化钾化肥价格（元/t），R3为氯化钾化肥含磷量（%），M为林分土壤有机质平均含量（%），C3为有机质价格（元/t）。2.2.4 森林生态系统固氮释氧功能价值评价方法指森林生态系统通过森林植被、土壤动物和微生物固定碳素、释放氧气的功能13。森林固碳释氧包括固碳和释氧2个指标。固碳释氧的价值评估要先计算出森林固碳释氧的物质量，然后再计算其价值量。本研究采用光合作用方程式来计算固碳释氧量；固碳（CO2）价值量的计算方法采用人工固定CO2成本法；释氧价值量计算方法采用工业制氧法。计算公式如下：1）固碳指标：Uc=ACc(1.63RcBa+Fsc) （9）式中：Uc为林分年固碳价值（元/a）；Cc为固碳价格（元/t）；Rc为二氧化碳中含碳质量分数，为27.27%；Ba为林分净生产力（t/haa）；Fsc为单位面积林分土壤年固碳量（t/ haa）。2）释氧指标：Uo=1.19CoABa （10）式中：Uo为林分年释氧价值（元/a），Co为氧气价格（元/t）。2.2.5 森林生态系统积累营养物质功能价值评价方法森林植物通过生化反应，在大气、土壤和降水中吸收N、P、K等营养物质并贮存在体内各营养器官的功能13。森林积累营养物质功能对降低下游面源污染及水体富营养化有重要作用。计算公式：Un=ABa(NtC1/R1+PtC1/R2+KtC2/R3) （11）式中：Un为林分年营养物质积累价值（元/a），Nt为林木含氮质量分数（%），Pt为林木含磷质量分数（%），Kt为林木含钾质量分数（%）。2.2.6 森林生态系统净化大气环境功能价值评价方法森林生态系统对大气污染物（如二氧化硫、氟化物、氮氧化物、粉尘、重金属等）的吸收、过滤、阻隔和分解，以及降低噪音、提供负离子和萜烯类（如芬多精）物质等功能13。本研究主要对森林提供负离子，吸收二氧化硫、氟化物、氮氧化物和阻滞粉尘指标来进行价值计算。 计算公式：1）提供负离子指标：Uni=5.2561015AHKni (Qni-600) /L （12）式中：Uni为林分年提供负离子价值（元/a），H为林分高度（m)，Kni为每个负离子生产费用（元），Qni为林分负离子密度（个/cm3 ),L为负离子寿命（分钟）。2）吸收二氧化硫指标：Uso=KsoQsoA （13）式中：Uso为林分年吸收二氧化硫价值（元/a），Kso为二氧化硫治理费用（元/kg），Qso为单位面积林分年吸收二氧化硫量（kg/haa）。3）吸收氮氧化物指标：Uhf=KhfQhfA （14）式中：Uhf为林分年吸收氮氧化物价值（元/a），Khf为氮氧化物治理费用（元/kg），Qhf为单位面积林分年吸收氮氧化物量（kg/haa）。4）吸收氟化物指标：Uno = KnoQnoA （15）式中：Uno为林分年吸收氟化物价值（元/a），Kno为氟化物治理费用（元/kg），Qno为单位面积林分年吸收氟化物量（kg/haa）。5）滞尘指标：Udc=KdcQdcA （16）式中：Udc为林分年滞尘价值（元/a），Kdc为降尘清理费用（元/kg），Qdc为单位面积林分年滞尘量（kg/haa）。2.2.7森林生态系统生物多样性保护功能价值评价方法森林生态系统为生物物种提供生存与繁衍的场所，从而对其起到保护作用的功能13。主要为物种保育功能价值。计算公式：Us=SsA （17）式中：Us为林分年物种保育价值（元/a），Ss为单位面积年物种损失的机会成本（kg/haa）。2.2.8 森林生态系统森林游憩功能价值评价方法森林生态系统为人类提供休闲和娱乐的场所，使人消除疲劳、愉悦身心、有益健康的功能而产生的价值13。森林游憩价值包括直接价值和间接价值，其间接价值远远大于直接价值。影响森林游憩价值的因素很多，其评估方法也很多。广为流行的评估方法是旅行费用法和条件价值法。但不管哪种方法，操作起来都很复杂，评估方法之间的评估结果也相差很大。本研究将综合考虑，选择合适的评估方法对株洲森林游憩功能价值进行评估。 3 结果3.1 森林生态系统提供林产品功能价值经计算，株洲市森林生态系统提供林产品总价值为67.721亿元（详见表3-1）。3.1.1 木材价值森林净初级生产量Q1（不含竹林、经济林和特灌林）采用黄怀雄、赵红艳（2010）的研究，不同林分的净生产力（森林的净干重增长量t/haa）如下:阔叶树5.44, 杉木5.25, 松木4.27, 竹林11.98, 经济林1.98, 灌木林4.1316；森林综合出材率R1取50%；森林面积为“十二五”森林资源二类调查数据；木材价格P1综合2013-2014年全国主要城市木价格行情数据为1350元/m3。则计算得出木材价值为34.795亿元/a。3.1.2 竹材价值竹林面积A竹根据“十二五”二类调查数据，株洲市竹林面积为72963.9ha；竹材市场价格P竹取2014年湖南省和广东省毛竹市场平均价格：12元/根；单位面积竹材产量R竹按株洲市竹林平均生产150根/haa计算。则计算得出竹材价值为1.314亿元/a。3.1.3 油茶价值油茶年平均产量R油根据株洲市林业局统计年报，2011、2012、2013年产量分别为7987.8吨、9984.8吨、12931.6吨，三年平均为10301.4吨；油茶市场价格P油取近几年株洲市油茶平均价格80000元/t。则计算得出油茶价值为8.241亿元/a。3.1.4 松脂价值可采松树面积A脂，根据国家松脂采集规程规定，胸径20cm以上的松树才能采集松脂，胸径20cm以上的马尾松和湿地松都已达到近成熟。根据“十二五”森林资源二类调查资料，株洲市近熟林以上马尾松和湿地松面积分别为31595.1ha和4472.0ha。单位面积可采松脂数量R脂，据有关专家研究，可采松脂马尾松和湿地松单株年贮量为4.5kg19，按1200株/ha计算，单位面积松脂年贮量为5.4t/haa。P脂为松脂市场价格（元/t），据调查，2012-2013年松脂平均市场价格为12000元/t。则计算得出松脂价值为23.371亿元/a。表3-1 株洲市森林生态系统提供林产品功能价值林产品种类森林面积（ha）提供林产品数量价格价值（亿元/a）木材481755.4257.75万m31350元/ m334.795竹材72974.91094.62万根12元/根1.314油茶81252.810301.4t80000元/t8.241松脂36067.1194762t12000元/t23.371合计67.7213.2 森林生态系统涵养水源功能价值森林生态系统涵养水源功能价值包括调节水量价值和净化水质价值。经评估，株洲市森林生态系统涵养水源量为460753.38万t/a，价值为416.521亿元/a（详见表3-2）。年平均降水量P：查阅株洲市气象资料得出株洲市年平均总降雨量为1550mm。年平均蒸散量E：本研究取我国亚热带森林蒸散率O.481O.66020的平均值0.57计算，得出株洲市森林的平均蒸散量为883.5mm。水的净化费用K：取当前株洲市水价 2.93元/t。水库建设单位库容投资价格Cv：根据森林生态系统服务功能评估规范，单位库容造价为611元/t。地表径流量C：根据韩素芸等的研究资料计算21，湖南省森林平均地表径流量（C）为3.4mm/a。林分面积A：二类调查最新数据为694847.5ha。表3-2 株洲市森林生态系统涵养水源功能价值功能面积（ha）涵养水源量（万t/a）价格（元/t）价值（亿元/a）调节水量694847.5460753.386.11281.520净化水质694847.5460753.382.93135.001合计694847.5460753.38416.5213.3 森林生态系统固土保肥功能价值 森林的固土保肥功能分为固土和保肥2个指标。经计算，株洲市森林生态系统固土保肥总价值为43.831亿元/a。具体计算如下。3.3.1固土价值无林地土壤侵蚀模数X2采用康文星，田大伦的研究成果，为37.58 t/haa22。有林地土壤侵蚀模数X1：综合采用欧阳志云和康文星、田大伦的研究成果，阔叶林地0.73 t/haa、杉木林地为2.17 t/haa、马尾松林地1.85 t/haa、国外松林地2.15 t/haa、特灌林地0.52 t/haa22，竹林地和经济林地0.5 t/haa2，其它林地以柏树林地、杨树林地、桉树林地和三杉林地土壤侵蚀量的平均值2.31 t/haa作为估算依据。挖取和运输单位体积土方所需费用Cs：根据森林生态系统服务功能评估规范，人工挖土方I和土类每100m3需42个工时，按平均每工时80元计算，得挖取单位体积土方费用为33.6元/m3。土壤密度p：为1.2t/m321。计算结果见表3-3.表3-3 株洲市森林生态系统固土量和固土价值计算表森林类型林分面积（A）（X2）（X1）X2-X1固土量（t/a）土方费用（Cs）(P)固土价值（亿元/a）杉木林20038737.582.1735.417095703.6733.61.21.987阔叶树181550.337.580.7336.856690128.5633.61.21.873马尾松56949.437.581.8535.732034802.0633.61.20.570国外松27322.537.582.1535.43968036.1833.61.20.271其它15546.237.582.3135.27548314.4733.61.20.154竹林72974.937.580.537.082705909.2933.61.20.758经济林81252.837.580.537.083012853.8233.61.20.844特灌林58864.437.580.5237.062181514.6633.61.20.611合计694847.525237262.7133.61.27.068由表3-3得出，株洲市森林生态系统固土量估算结果为2523.73万t/a，固土价值为7.068亿元/a。3.3.2 保肥价值 林分土壤平均氮含量（N）为0.1775%、磷含量（Ps）为0.0569%、钾含量（Ks）为1.7204%、有机质含量（M）为3.8204%21。磷酸二铵化肥价格（C1）、氯化钾化肥价格（C2）、有机质价格（C3）：采用农业部中国农业信息网2014年秋季平均价格,分别为为2700元/t、2200元/t;600元/t (草炭土中含有机质90%),折合为有机质价格为667元/t)。市场相应产品磷酸二铵化肥含氮量（R1）为14%、磷酸二铵化肥含磷量（R2）15.01%、氯化钾化肥含钾量（R3）50%。表3-4 株洲市森林生态系统保肥量估算表森林类型森林面积（ha）固土量（t/a）保N量（t/a）保P量（t/a）保K量（t/a）保有机质量（t/a）杉木林2003877095703.6712594.874037.46122074.49271084.26阔叶树181550.36690128.611874.983806.68115096.97255589.67马尾松56949.42034802.13611.771157.8035006.7477737.58国外松27322.5968036.21718.26550.8116654.0936982.85其它15546.2548314.47973.26331.999433.2020947.81竹林72974.92705909.294802.991539.6646552.46103376.56经济林81252.83012853.825347.821714.3151833.14115103.07特灌林58864.42181514.663872.191241.2837530.7883342.59合计694847.525237262.8144796.1414379.99434181.87964164.39表3-5株洲市森林生态系统保肥价值估算表森林类型保N价值（万元/a）保P价值（万元/a）保K价值（万元/a）保有机质价值（万元/a）保肥价值（万元/a）杉木林24290.1067262.58653712.77618081.320103346.788阔叶树22901.756847.4650642.6717047.8397439.71马尾松6965.562082.6515402.975185.1029636.28国外松3313.79990.807327.802466.7614099.15其它1877.001597.1844150.6081397.2198022.012竹林9262.9092769.54220483.0826895.21739410.75经济林10313.6533083.70222806.5827677.37543881.312特灌林7467.802232.8216531.545558.9531791.11合计86392.56925866.744191058.02864309.771367627.112由表3-4和表3-5得知，株洲市森林生态系统保肥量及价值估算结果为：保N量44796.14万t/a，价值86392.569万元/a；保P量14379.99万t/a，价值25866.744万元/a；保K量434181.87万t/a，价值191058.028万元/a；保有机质量964164.39万t/a，价值64309.771万元/a。合计保肥价值36.763亿元/a。3.4 森林生态系统固碳释氧功能价值经计算，株洲市森林生态系统固碳释氧功能价值为95.664亿元/a。其中：固碳量为2390185.07t/a,价值为28.682亿元/a ；释氧量为4465434.77t/a,价值为66.982亿元/a 。3.4.1 固碳指标单位面积林分土壤年固碳量Fsc：本研究采用韩素芸、田大伦等研究成果，林分土壤年固碳量为0.6561.647t/hm，其中：阔叶树1.647、马尾松0.777、湿地松0.771，杉木、柏木、三杉和特灌林为0.65621，竹林和经济林取0.6561.647的中值1.1515。固碳价格Cc采用瑞典的碳税率150美元/t (折合人民币为12OO元/t )。林分净生产力Ba采用黄怀雄、赵红艳的研究结果：阔叶树5.44, 杉木5.25, 松木4.27, 竹林11.98, 经济林1.98, 灌木林4.1316表3-6 株洲市森林生态系统固碳量及价值计算表森林类型面积（ha）1.63RcBa+FscRcBaFsc固碳量（t/a）Cc固碳价值 (万元/a)杉木林2003872.989630250.27275.250.656599083.04120071889.964 阔叶树181550.34.065085440.27275.441.647738017.48120088562.098 马尾松56949.42.675019270.27274.270.777152340.74120018280.889 国外松27322.52.669019270.27274.270.77172924.2812008750.913 其它15546.22.794049810.27274.810.65643436.8612005212.423竹林72974.96.476621980.272711.981.1515472630.84120056715.701 经济林81252.82.031611980.27271.981.1515165074.16120019808.899特灌林58864.42.491789130.27274.130.656146677.67120017601.321合计694847.52390185.07286822.208 3.4.2 释氧指标：氧气价格Co采用2014年全国氧气平均价格：1500元/t。林分净生产力Ba采用黄怀雄、赵红艳的研究结果16。表3-7 株洲市森林生态系统释氧量及价值计算表森林类型面积（ha）Ba(t/haa)释氧量（t/a）Co(元/t)释氧价值(万元/a)杉木林2003875.251251861.561500187779.233阔叶树181550.35.441175284.021500176292.603 马尾松56949.44.27289376.99150043406.548国外松27322.54.27138833.82150020825.073 其它15546.24.8188984.89150013347.734竹林72974.911.981040344.771500156051.715 经济林81252.81.98191447.85150028717.177 特灌林*58864.44.13289300.87150043395.13合计694847.54465434.771500669815.2133.5 森林生态系统积累营养物质功能评价经计算，株洲市森林生态系统积累营养物质价值5.34亿元/a。其中：积累N营养物质21789.9t/a,价值为4.20亿元/年；积累P营养物质2276.6t/a,价值为0.41亿元/年；积累K营养物质16848.5t/a,价值为0.74亿元/年。具体计算如下。林木含氮质量分数Nt、含磷质量分数Pt、含钾质量分数Kt：采用莫江明、张晓东、项文化、郑荣、曹根群等人对不同森林生态系统的植物体营养元素含量的研究成果23，见表3-8：表3-8 林木含N、P、K质量分数表森林类型亚热带常绿阔叶林亚热带常绿针叶林亚热带竹林亚热带灌木林N(%)0.8260.4200.6510.420P(%)0.0350.0750.0790.031K(%)0.6330.2130.5500.649林分净生产力Ba：采用黄怀雄、赵红艳的研究成果16。肥料价格：采用农业部中国农业信息网2014年秋季平均价格,磷酸二铵C12700元/t，氯化钾C12200元/t。化肥中营养物质含量：N素R1为14.00%，P素R2为15.01%，K素R3为50.00%。计算结果见表3-9.表3-9 株洲市森林生态系统积累营养物质价值计算表树种积累N素（t/a）积累N素价值（万元/a）积累P素（t/a）积累P素价值（万元/a）积累K素（t/a）积累K素价值（万元/a）积累营养物质价值（万元/a）杉木林4418.58521.39789.01419.252240.8985.9510926.59阔叶树8157.915733.1345.7621.856251.72750.7519105.70马尾松1021.31969.65182.4328.1518.0227.922525.67国外松490.0945.0087.5157.40248.5109.341121.74其它林314.1605.7656.1100.91159.370.09776.76竹林5691.310976.08690.61242.254808.32115.6514333.98经济林675.71303.1449.989.761044.1459.401852.3特灌林*1021.11969.2675.4135.631577.8694.232799.12合计21789.942023.382276.64095.1516848.57413.3353441.863.6 森林生态系统净化环境功能价值经评估，株洲市森林生态系统净化环境功能总价值为25.52亿元/a。3.6.1 提供负氧离子价值林分高度H采用最新二类调查数据；负离子生产费用Kni采用韩素芸等的研究资料，为5819元10-18/个21。林分负离子密度Qni 由抽样实测获得。株洲市森林负氧离子平均浓度为：杉木林1067个/cm3、阔叶树1700个/cm3、松树林833个/cm3、竹林1300个/cm3、经济林1150个/cm3、灌木林1150个/cm3，其它林分按经济林计算。负离子寿命（L）：10分钟。 根据公式（12）计算得出株洲市森林生态系统产生负氧离子323.961014亿个/a,价值为1147万元/a。3.6.2 吸收二氧化硫价值二氧化硫治理费用Kso采用森林生态系统服务功能评估规范，为1.20元/kg13；森林类型吸收二氧化硫的能力Qso采用中国生物多样性国情研究报告资料和向会娟，曹明宏的研究成果，阔叶树对二氧化硫的吸收能力为88.65kg/ha，针叶林为215.60 kg/ha，柏类为41.16 kg/ha，杉类为117.6 kg/ha，松类为117.6 kg/ha24 25。针阔混交林取针叶林和阔叶林的平均值152.13（本研究中，其他、竹林、经济林、灌木林Qso都按针阔混交林来计算）。根据公式（13）计算得出株洲市森林生态系统吸收二氧化硫量为8435.31t/a,价值为1.01亿元/a 。3.6.3 吸收氟化物价值氮氧化物治理费用Khf采用森林生态系统服务功能评估规范，为0.63元/kg13。单位面积林分年吸收氮氧化物量Qhf采用向会娟，曹明宏的研究成果，单位面积森林吸收氮氧化物量为6.00kg/hma24。根据公式（13）计算得出株洲市森林生态系统吸收氮氧化物量为4169.09t/a,价值为263万元/a。3.6.4 吸收氟化物价值氟化物治理费用Kno采用森林生态系统服务功能评估规范，为0.69元/kg13；单位面积林分年吸收氟化物量Qno采用韩素芸、田大伦等的研究成果，阔叶树4.65、针叶树0.5021，其他、竹林、经济林、特灌林按平均值2.575计算。根据公式（14）计算得出株洲市森林生态系统吸收氟化物量为1575.28t/a,价值为109万元/a。3.6.5 滞尘价值降尘清理费用Kdc采用森林生态系统服务功能评估规范，为015元/kg13；单位面积林分年滞尘量Qdc采用中国生物多样性国情研究报告，针叶林的滞尘能力为33200kg/haa,阔叶林的滞尘能力为10110kg/haa25，其他、竹林、经济林、特灌林Qdc都按上述两种的平均值（21660kg/haa）计算。根据公式（15）计算得出株洲市森林生态系统滞尘量为16237311t/a,价值为24.356亿元/a。3.7 森林生态系统生物多样性保护功能价值 物种损失机会成本Ss按照Shannon-Wiener指数划分为6个等级：当指数1时，Ss为5000元/haa；当指数12时，Ss为10000元/haa；当指数23时，Ss为20000元/haa；当指数34时，Ss为30000元/haa；当指数45时，Ss为40000元/haa；当指数5时，Ss为50000元/haa13。本研究采用韩素芸等的研究成果，杉木林、马尾松林为10000元/haa，阔叶林为30000元/haa，其他、竹林、经济林和特灌林Ss都按其平均值20000元/haa21来计算。根据公式（16）计算得出株洲市森林生态系统生物多样性保护价值为128.65亿元/a。3.8 森林生态系统森林游憩功能价值全市现有森林公园7个，经营面积28169.5hm2。分别为炎陵县神农谷、茶陵县云阳、攸县攸州3个国家森林公园，炎陵县湘山、醴陵市仙岳山、石峰区九郎山3个省级森林公园和醴陵市章龙山1个县级森林公园。2014年全市森林公园接待游客500万人次，实现旅游收入19838万元,门票收入600万元。根据掌握的资料和结合株洲实际情况，本研究采用北京林业大学经济管理学院张颖教授的研究成果。他在“森林游憩资源价值评价模型研究”中，根据我国1990-2010年森林旅游的统计资料（GDP、林业总产值、营林基本建设投资、森林公园总面积、森林公园面积变化情况、旅游人次、旅游收入等），以森林游憩资源价值评价为研究对象，采用最优控制的方法，分析森林游憩的变化与经济发展的关系，建立森林游憩与经济发展的差分方程，进而求出了森林游憩资源价值核算的最优价格，为1.044万元/ha。2011-2014年，我国的森林旅游业的增长速度为10%。由此推算出，2014年我国森林游憩资源的最优价格为1.529万元/ha。以这个数据作为计算参数,与株洲市森林公园总面积28169.5ha相乘，得出2014年我市森林游憩功能的价值为4.307亿元。4. 结果与讨论4.1 株洲市森林生态系统生态服务功能实物量及其价值量根据以上评价方法,得出株洲市森林生态系统8个方面23个指标的生态服务功能的实物量和价值量（详见表4-1）。表4-1 株洲市森林生态系统服务功能实物量和价值量汇总表功能指标实物量价值（万元）小计（万元）林产品木材257.75万m3347950.00677210.00竹材1094.62万根13140.00油茶10301.4t82410.00松脂194762t233710.00涵养水源调节水量460753.38万t2815200.004165210净化水质460753.38万t1350010.00固土保肥固土2523.73万t70680.00438307.11保肥保N44796.14万t367627.112保P14379.99万t保K434181.87万t保有机质964164.39万t固碳释氧固碳2390185.07t286822.21956637.42释氧4465434.77t669815.21积累营养物质积累N21789.9t42023.3853531.86积累P2276.6t4095.15积累K16848.57413.33净化环境提供负离子323.961014亿个1147.32255200.71吸收二氧化硫8435.31t10122.37吸收氮氧化物4169.09t262.66吸收氟化物1575.28t108.69滞尘16237311t243559.67保护生物多样1286500.001286500.00森林游憩43071.1743071.17合计7875668.274.2 生态服务功能与国民经济指标对比按照株洲市“十二五”森林资源二类调查数据评估，株洲市森林生态系统服务功能总价值为787.567亿元/a。根据2013年株洲市国民经济和社会发展统计公报，2013年全市生产总值1948亿元，其中农林牧渔业产值154.7亿元。由此可知，株洲市森林生态系统服务功能价值占株洲市2013年国民经济生产总值的40.43%，是同期农业总产值的5.09倍。由于本研究只选取了森林生态系统服务功能中的8项功能，计算得出的功能价值量远低于实际功能价值量。因此，其在国民经济中的占有量还应有所提高。4.3 森林生态系统功能价值结构分析表4-2 株洲市森林生态系统服务各功能价值量百分比功能提供林产品涵养水源固土保肥固碳释氧积累营养物质净化大气环境生物多样性保护森林游憩价值（亿元/a）67.721416.52143.83195.6645.35325.520128.6504.307占比（%）8.6052.895.5712.150.683.2416.340.55排序41537628从表4-2得出株洲市森林生态系统各项功能价值在总功能价值量中所占的比例，其大小排列顺序是：涵养水源(52.89)生物多样性保护 (16.34)固碳释氧(12.15)提供林产品(8.60)固土保肥(5.57)净化大气环境(3.24)积累营养物质(0.68) 森林游憩(0.55)。这与王兵、李少宁、郭浩对江西省森林生态系统服务功能价值研究结果相近，江西省7项森林生态系统服务功能价值的贡献大小顺序依次均为:涵养水源（63.161%）固碳制氧（18. 523%）保护生物多样性（9. 586%）净化大气环境（5. 183%）保育土壤（2. 280%）林木营养积累（1. 264%）森林游憩与生态文化（0. 014%）17。株洲市森林生态系统各项服务功能中，涵养水源功能价值占总价值的52.84%，说明森林生态系统在水循环调节中占有举足轻重的地位，这对调节气候，降低水、旱灾害，保护农业生产和人民生活起到了重要作用。在大气调节功能(固碳释氧和净化大气环境)方面，其价值占总价值的15.39，其中吸收二氧化硫量（8435.31t/a）和氮氧化物量（4169.09t/a）分别为2013年全市二氧化硫排放4.43万t和氮氧化物排放量（3.65万t）的19.04%和11.42%，对改善株洲市大气环境起到积极作用。在提供林产品价值方面，其价值只占总价值的8.60，说明森林生态系统提供的间接经济价值远大于传统意义上森林生态系统提供木材和林副产品等直接经济价值，提供林产品不是森林生态系统服务价值的主要部分。因此，充分发挥株洲市森林生态系统其他的生态服务功能，才是发挥其生态系统服务价值的最佳途径。4.4 森林生态系统各林分类型服务功能分析表4-3 株洲市各森林类型服务功能价值核算表（不含森林游憩功能）森林类型面积（ha）单位面积价值（元/haa）单位面积价值排序总价值（万元/a）占总价值比重（%）总价值排序杉木林20038710211972046334.326.132阔叶林181550.312128832201992.928.111马尾松56949.41345031765985.329.785国外松27322.51090664297995.383.807其它15546.21088735169256.452.168竹林72974.91227552895803.1711.443经济林81252.81061666862630.211.014特灌林58864.41006568592507.547.566合计694847.51127237832505.31004-4 株洲市各森林类型纯生态服务功能价值核算表（不含森林游憩功能）森林类型面积（ha）单位面积价值（元/haa）单位面积价值排序总价值（万元/a）占总价值比重（%）总价值排序杉木林2003879473661898392.426.532阔叶林181550.311363822063106.928.831马尾松56949.4925487527052.737.376国外松27322.5924558252610.513.537其它15546.21021093158740.92.228竹林72974.91209551882667.6912.343经济林81252.896024578021910.904特灌林58864.41006564592507.548.285合计694847.51030197155297.6100由表4-3可知，株洲市森林生态系统单位面积服务功能价值大小排序为：马尾松林竹林阔叶林国外松林其它林分经济林杉木林特灌林。其中马尾松、竹林和阔叶林超过了单位面积生态服务功能价值112723元/haa的总平均值。马尾松排名第一，这主要是由于马尾松除提供纯生态功能外，还可提供木材和松脂产品；竹林、阔叶树、国外松、其它林分以及经济林、杉木林都提供一定量的林产品，只有特灌林没有提供林产品。各森林类型服务功能总价值大小顺序为：阔叶林杉木林竹林经济林马尾松特灌林国外松其他林分 。阔叶林和杉木林排名一至二位，与其面积在有林地总面积中的比例基本相符。由表4-4可知，株洲市森林生态系统单位面积纯生态服务功能价值大小排序为：竹林阔叶树其它林分特灌林经济林杉木林马尾松林国外松。其中竹林和阔叶林的单位面积纯生态服务功能价值超过了103019元/haa的全市单位面积纯生态服务功能价值总平均值。由此可以看出，竹林和阔叶林的单位面积纯生态服务功能价值大，应加大发展力度。株洲市森林生态系统纯生态服务功能总价值大小排序为：阔叶林杉木林竹林经济林特灌林马尾松国外松其它林分。这一方面基本与它们的面积大小相关，也与它们的单位面积纯生态服务功能价值大小排序相关。4.5 各县市区森林生态系统服务功能价值分析表4-5 株洲市各县市区森林生态系统服务功能价值比较县市区森林面积（ha）单位面积价值（元/haa）单位面积价值排序总价值（