海南省生态产品总值（GEP）核算技术规范TechnicalspecificationforAccountingGrossEcosystemProduct(GEP)HainanProvince（试行）海南省生态环境厅2023年3月I目次前言............................................................................I引言...........................................................................II1范围................................................................................12规范性引用文件......................................................................13术语和定义..........................................................................14总则................................................................................34.1核算区域与核算单元..............................................................34.2核算原则........................................................................34.3核算程序........................................................................35评估指标体系........................................................................46调查与评估使用数据要求..............................................................67陆地生态产品总值核算方法............................................................67.1物质产品供给服务................................................................67.2调节服务........................................................................77.3文化服务.......................................................................158海洋和滨海湿地生态产品总值核算方法.................................................158.1物质产品供给服务...............................................................158.2调节服务.......................................................................178.3文化服务.......................................................................239附则...............................................................................239.1质量控制.......................................................................239.2成果汇总.......................................................................239.3成果应用.......................................................................24附录A（规范性）生态产品清单...................................................25附录B（资料性）生态产品实物量和价值量核算关键参数.............................26B.1水源涵养服务核算主要参数.......................................................26B.2土壤保持服务核算主要参数.......................................................29B.3洪水调蓄服务核算主要参数.......................................................31B.4固定二氧化碳服务核算参数.......................................................32B.5局部气候调节服务核算主要参数...................................................35B.6空气净化服务核算主要参数.......................................................36B.7水体净化服务核算主要参数.......................................................36B.8海洋废弃物处理服务主要参数.....................................................37B.9物种保育更新服务核算参数.......................................................38B.10价格参数......................................................................39附录C（资料性）生态产品总值核算基础数据要求....................................41II附录D（资料性）生态产品总值核算报告编写大纲...................................46附录E（资料性）生态系统划分...................................................47参考文献.......................................................................49I前言本文件按照GB/T1.1—2020《标准化工作导则第1部分：标准化文件的结构和起草规则》的规定起草。本文件由海南省生态环境厅组织编制并归口。本文件起草单位：海南省环境科学研究院（国家生态文明试验区（海南）研究中心）。本文件主要起草人：胡小飞、王晨野、穆晓东、邢巧、王敏英、吴晓晨、许积层、马字伟、邬乐雅、杨安富、覃茂运、侯德佳、蓝才鑫、左永令、王艺臻II引言为贯彻落实《关于建立健全生态产品价值实现机制的意见》《国家生态文明试验区（海南）实施方案》等文件精神，按照《海南省建立健全生态产品价值实现机制实施方案》要求，持续推动生态产品价值核算工作向纵深发展，需构建符合海南省区域特征的生态产品价值核算指标体系和方法，指导和规范海南省生态产品总值核算工作，提高核算的科学性、规范性、可操作性，推动核算结果可追溯、可核查、可比较，制定本技术规范。本技术规范基于海南省陆海区域特征，构建适用于海南省陆地生态系统和海洋生态系统的生态产品总值核算指标体系，界定本地化关键参数，规范生态产品总值核算方法和核算数据要求。本文件可为政府决策、绩效考核评价，以及生态保护补偿、生态环境损害赔偿、经营开发融资、生态资源权益交易、生态银行等生态产品价值实现领域提供参考支撑。1海南省生态产品总值（GEP）核算技术规范1范围本文件规定了海南省陆地、滨海湿地与海洋生态系统生态产品总值核算技术流程、指标体系、功能量核算方法、定价方法和价值量核算方法等内容。本文件适用于海南省省、市、县（区）、镇（乡）、村等各级行政区域的生态产品总值核算，其他功能相对完整的生态自然地理区域以及由不同生态系统类型组合而成的地域单元，可参照本文件开展生态产品总值核算。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB3838—2002地表水环境质量标准GB/T38582-2020森林生态系统服务功能评估规范GB/T28058海洋生态资本评估技术导则HY/T0305-2021养殖大型藻类和双壳贝类碳汇计量方法碳储量变化法HJ1169全国生态状况调查评估技术规范——湿地生态系统野外观测3术语和定义下列术语和定义适用于本文件。3.1陆地生态系统terrestrialecosystem指陆地表面一定空间范围内植物、动物和微生物群落及与非生物环境相互作用形成的功能整体，包括森林、草地、农田、湿地、城市（镇）等生态系统类型。3.2海洋生态系统marineecosystem指由海洋中的生物群落与其非生物环境相互作用形成的功能整体。3.3滨海湿地生态系统coastalwetlandsecosystem指分布在低潮时海平面以下6米处至大潮线之上，包括与内河流域相连的淡水或半咸水湖沼以及海水上溯未能抵达的入海河的河段范围内的生物群落及其无机环境相互作用形成的有机整体。3.4森林生态系统forestecosystem指以乔木、灌丛和其中的动物为主体的生物群落与其非生物环境相互作用形成的功能整体。3.5农田生态系统farmlandsecosystem2指以农作物为主体的生物群落与其非生物环境相互作用形成的功能整体。3.6草地生态系统grasslandsecosystem指以草本植物和食草动物为主体的生物群落与其非生物环境相互作用形成的功能整体。3.7内陆湿地生态系统inlandwetlandsecosystem指由陆地或水域相互作用区域内的各种生物与其非生物环境相互作用形成的兼顾水域和陆地生态系统特征的功能整体，包括沼泽湿地、湖泊湿地、河流湿地等自然湿地，以及库塘、盐田及其他重点保护野生动物栖息地或者重点保护野生植物原生地等人工湿地等。3.8城市生态系统urbanecosystem指由城市居民、生活在其中的动植物与其非生物环境相互作用而形成的功能整体。3.9生态系统服务ecosystemservices人类从生态系统中得到的惠益，包括物质供给、调节服务、文化服务以及支持服务。支持服务属于中间过程服务，不属于最终服务，为了避免重复计算，在生态系统生产总值核算中不考虑支持服务。3.10生态产品1)ecologicalproduct生态系统为经济活动和其他人类活动提供且被使用的货物与服务贡献，包括物质产品供给、调节服务和文化服务三类。3.11物质产品供给服务materialprovisioningservices生态系统在一定时间内为人类提供并被使用的各种物质产品，如粮食、油料、蔬菜、水果、木材、生物质能、水产品、中草药、牧草、花卉、种质资源等。3.12调节服务regulatingservices指生态系统为维持或改善人类生存环境提供的惠益，如水源涵养、土壤保持、防风固沙、海岸带防护、洪水调蓄、空气净化、水体净化、固碳、局部气候调节、噪声消减等。3.13文化服务culturalservices生态系统为提高人类生活质量提供的非物质惠益，如精神感受、灵感激发、旅游观光、休闲娱乐和美学体验等。3.14实物量bio-physicalvalue生态系统产品与服务的物理量，如粮食产量、洪水调蓄量、土壤保持量、固碳量与景点旅游人数等。1)“生态产品”是具有中国特色的概念。2010年12月国务院发布的《全国主体功能区划》（国发（2010）46号）中提出“生态产品指维系生态安全、保障生态调节功能、提供良好人居环境的自然要素，包括清新的空气、清洁的水源和宜人的气候等”和“生态功能区提供生态产品的主体功能主要体现在：吸收二氧化碳、制造氧气、涵养水源、保持水土、净化水质、防风固沙、调节气候、清洁空气、减少噪音、吸附粉尘、保护生物多样性、减轻自然灾害等”。这一界定主要强调生态系统的调节和维持服务。随着时间推移，将物质供给和文化服务纳入生态产品范畴已成为主流共识。国际上与之最接近的一个概念是联合国等国际组织最新制定的《环境经济核算体系——生态系统核算》（SEEAEA）中的“生态系统最终服务”，指生态系统为经济活动和其他人类活动提供且被使用的最终产品，包括物质供给、调节服务、文化服务。33.15价值量monetaryvalue生态产品的货币价值。3.16生态产品总值grossecosystemproduct,GEP一定空间范围内各类生态系统在核算期内提供的所有生态产品的货币价值之和。4总则4.1核算区域与核算单元陆地生态系统以自然资源和规划主管部门确定的海岸线以上陆域作为核算区域；海洋生态系统以海南省管辖的近岸海域作为核算区域2)。核算单元的大小根据核算区域的空间范围确定。省级行政区以30m×30m的空间网格作为基本核算单元，地级市、县（市、区）级行政区建议以10m×10m的空间网格作为基本核算单元。4.2核算原则4.2.1客观性原则核算时应优先使用实测数据，现有监测体系资料不能满足核算需求的，相关部门单位应积极完善监测体系，开展补充性调查监测，力求做到基础资料可靠、核算方法科学、核算结果准确可比。4.2.2基于交换价值原则采用与国民账户体系一致的估价方法，基于交换价值测算生态产品价值，以支持与标准国民账户中的信息整合。4.2.3实用性原则根据海南省生态系统特点、核算目的和数据可获得性，确定适合海南省区域特征的核算内容、指标和方法。4.2.4本地化原则基于热带雨林国家公园、近岸海域以及市县等生态产品总值核算结果，结合海南省区域特征，建立适合海南省本地化参数表，实现核算参数本地化。4.2.5开放性原则根据生态产品价值核算最新研究成果，发展和完善价值核算的指标与方法。4.3核算程序生态产品总值核算主要工作程序包括：a)确定核算地域范围。根据核算目的，确定生态产品总值核算的区域范围。可按下列方式确定：行政区域单元、功能相对完整的生态系统地域单元、由不同生态系统类型组合而成的特定地域单元等。2)海南岛周边近岸管理海域即海南岛海岸线至领海外缘线之间的海域，琼州海峡以广东省和海南省间海域行政区域界线为界。4b)明确生态系统类型与分布。调查分析核算区域内生态系统类型、面积与分布，绘制生态系统空间分布图。c)编制生态产品目录清单。根据生态系统类型与核算用途，调查分析核算范围内生态产品种类，确定生态产品总值核算的重点，编制生态产品目录清单，确定核算基准时间。d)收集资料与补充调查。收集生态产品总值核算所需相关文献资料、各类生态系统监测与统计数据以及基础地理与地形图件，开展必要的实地观测调查，进行数据预处理以及参数本地化。e)核算生态产品实物量。按照本规范规定的技术方法与核算模型，核算各类生态产品实物量。当区域参数无基础数据支撑时，可参考本规范提供的技术参数。f)确定生态产品单位价值。按照本规范规定的技术方法与核算模型，运用市场价值法、替代成本法等价值核算方法，采用当年价确定每一类生态产品的参考价格。涉及多年对比时可以采用基准年不变价。g)核算生态产品价值量。根据生态产品实物量及其单位价值，对生态系统提供的供给服务、调节服务和文化服务价值进行核算。h)核算生态产品总值。将核算区范围内各类生态产品价值加总，得到生态产品总值，核算按式（1）。𝐺𝐸𝑃=𝐸𝑃𝑉+𝐸𝑆𝑉+𝐸𝐶𝑉+𝐸𝑃𝑉𝑆+𝐸𝑆𝑉𝑆+𝐸𝐶𝑉𝑆··········································(1)式中：𝐺𝐸𝑃——生态产品总值（万元/a）；𝐸𝑃𝑉、𝐸𝑃𝑉𝑆——分别是陆地、海洋和滨海湿地物质产品供给服务价值（万元·a-1）；𝐸𝑆𝑉、𝐸𝑆𝑉𝑆——分别是陆地、海洋和滨海湿地调节服务价值（万元·a-1）；𝐸𝐶𝑉、𝐸𝐶𝑉𝑆——分别是陆地、海洋和滨海湿地文化服务价值（万元·a-1）。i)开展核算结果分析。对核算区生态产品总值的空间分布、区域特征及单位面积或人均生态产品总值等指标进行分析，完成核算报告。5评估指标体系海南省生态产品总值核算指标体系由物质产品供给服务、调节服务和文化服务3大类服务构成。其中，物质产品供给服务主要包括农业产品、林业产品、畜牧业产品、渔业产品、淡水资源、种质资源；调节服务主要包括水源涵养、土壤保持、固碳、氧气释放、局部气候调节、洪水调蓄、海岸带防护、物种保育更新、空气净化、水体净化、废弃物处理、负氧离子释放；文化服务主要包括旅游康养。生态产品总值核算指标体系与核算方法见表1，各类生态系统类型生态产品总值核算指标见表2。5表1生态产品总值核算指标体系与核算方法核算区域类型一级指标二级指标实物量核算方法价值量核算方法陆地生态系统物质产品供给服务农林牧渔（淡水）产品统计法增加值法或市场价值法淡水资源统计法市场价值法生态能源统计法市场价值法种质资源统计法市场价值法调节服务水源涵养水量平衡法影子工程法土壤保持修正通用土壤流失方程替代成本法洪水调蓄水量储存模型影子工程法固碳机理模型市场价值法氧气释放机理模型替代成本法空气净化污染物净化模型替代成本法水体净化污染物净化模型替代成本法局部气候调节实测法、蒸散模型替代成本法负氧离子释放机理模型替代成本法海岸带防护统计调查法替代成本法物种保育更新Shannnon-Weine指数替代成本法文化服务旅游康养统计调查法统计调查法、旅行费用法海洋和滨海湿地生态系统物质产品供给服务养殖生产统计法市场价值法捕捞生产统计法市场价值法种质资源统计法市场价值法调节服务固碳机理模型市场价值法氧气释放机理模型替代成本法废弃物处理统计监测法治理成本法局部气候调节实测法、蒸散模型替代成本法海岸带防护统计调查法影子工程法物种保育更新统计法保育价值法文化服务旅游康养统计法统计调查法、旅行费用法注：为可选项，根据实际情况选择是否核算。6表2各类生态系统生态产品总值主要核算指标类别核算指标森林草地农田内陆湿地城市滨海湿地海洋物质产品供给服务物质产品供给√√√√√√√调节服务水源涵养√√×√×××土壤保持√√√√×××洪水调蓄√√×√×××空气净化√√√√√√×水体净化×××√√√×固碳√√√√√√√氧气释放√√×√√√×局部气候调节√√×√√√√废弃物处理××××××√负氧离子提供√√×√√√×海岸带防护√×××√√×物种保育更新√√√√×√√文化服务旅游康养√√√√√√√注1：×为不评估；√为评估。注2：城市生态系统中的林地、草地、农用地、湿地等产生的生态服务功能按照本标准提供的方法核算后，按城市边界（城市建成区范围）进行汇总，得到城市生态产品总值。注3：为可选项，根据实际情况选择是否核算。6调查与评估使用数据要求生态产品总值评估使用的地理空间数据、监测数据、统计数据、实地调查数据以及参考文献数据等数据来源、单位、频次等应符合附录C的规定。7陆地生态产品总值核算方法7.1物质产品供给服务7.1.1实物量核算主要对农业、林业、畜牧业、渔业（淡水）、种子产品、（淡）水资源产品等直接利用物质产品以及生态能源等转化利用产品产量进行统计，不包括野生动物、野生植物或小规模自产自用的其它生态产品。其中一些以干货进行统计的产品如食用菌、茶叶等应按照经验系数换算为湿重核算。直接利用物质产品总产量按式（2）进行求和；生态能源等转化利用供给产品按式（3）。7Epro=∑𝐸𝑖𝑛𝑖=1··········································································(2)式中：Epro——物质产品总产量（t·a-1）；𝑖——核算区物质产品种类，𝑖=1，2，3，……，n；𝐸𝑖——第𝑖种物质产品的产量（t·a-1）。生态能源等转化利用供给产品按式（3）。𝐸𝑒𝑒=∑𝐸𝑒𝑒𝑖𝑛𝑖=1·········································································(3)式中：𝐸𝑒𝑒——生态能源总产量或使用量（kW·h·a-1）；𝑖——核算区生态能源的类型，𝑖=1，2，3，……，n；𝐸𝑒𝑒𝑖——𝑖类生态能源类型的产量（kW·h·a-1）。7.1.2价值量核算核算地区可获得本项产品增加值时，推荐采取增加值进行核算，没有增加值数据时，用市场价格法核算产品产值进行代替。物质产品供给服务价值量按式（4）。𝐸𝑃𝑉=（∑𝐸𝑖×𝑃𝑖𝑛𝑖=1+∑𝐸𝑒𝑒𝑖×𝑃𝑒𝑖𝑛𝑖=1）×10−4···········································(4)式中：𝐸𝑃𝑉——物质产品供给的总价值（万元·a-1）；𝑃𝑖——第𝑖类直接利用产品的价格（元·kg-1），单位视具体产品而定；𝑃𝑒𝑖——第𝑖类转化利用产品的价格（元·（kW·h）-1）7.2调节服务7.2.1水源涵养服务7.2.1.1实物量核算利用水量平衡模型计算水源涵养量（𝑄𝑤𝑟），核算按式（5）。𝑄𝑤𝑟=[∑𝐴𝑗𝑛𝑗=1×[∑∑（12𝑚=1𝑃𝑖,𝑚−𝑄𝐹𝑖,𝑚−𝐴𝐸𝑇𝑖,𝑚）𝑛𝑖=1]×10−3···································(5)式中：𝑄𝑤𝑟——水源涵养量（m3·a-1）；𝐴𝑗——表示j类生态系统面积（m2）；𝑃𝑖,𝑚——核算区内𝑖像元的m月多年平均降雨量（mm·a-1），建议采用近5年移动平均值；𝑄𝐹𝑖,𝑚——核算区内的𝑖像元m月地表径流量（mm·a-1），建议采用近5年移动平均值，或参考附录B.1；𝐴𝐸𝑇𝑖,𝑚——核算区内𝑖像元m月实际蒸散发量（mm·a-1），核算见附录B.1；𝑖——核算区内的第𝑖个像元；𝑛——核算区内的总像元个数。7.2.1.2价值量核算8采用影子工程法，以建设蓄水量与生态系统水源涵养量相当的水利设施所需成本作为水源涵养价值。核算按式（6）𝑉𝑤𝑟=𝑄𝑤𝑟×𝐶𝑤𝑒×10−4··································································(6)式中：𝑉𝑤𝑟——蓄水保水价值（万元·a-1）；𝐶𝑤𝑒——水库单位库容的工程造价（元·m-3），见附录B.10。7.2.2土壤保持服务7.2.2.1实物量核算采用通用土壤流失方程（RUSLE）核算土壤保持服务实物量。核算按式（7）。𝑄𝑠𝑟=∑[𝑅×𝐾×𝐿×𝑆×（1−𝐶×𝑃）×𝐴𝑗𝑛𝑗=1]············································(7)式中：𝑄𝑠𝑟——核算区土壤保持量（t·a-1）；𝐴𝑗——j类生态系统的面积（m2）；𝑅——降雨侵蚀力因子（MJ·mm·hm-2·h-1·a-1），用多年平均年降雨侵蚀力指数表示，核算见附录B.2；𝐾——土壤可蚀性因子（t·h·MJ-1·mm-1），用标准样方上单位降雨侵蚀力所引起的土壤流失量来表示，核算见附录B.2；𝐿——坡长因子（无量纲），核算见附录B.2；𝑆——坡度因子（无量纲），核算见附录B.2；𝐶——盖度和管理因子（无量纲），核算见附录B.2；𝑃——水土保持措施因子（无量纲），核算见附录B.2；n——核算区生态系统类型总数。7.2.2.2价值量核算采用替代成本法进行土壤保持价值量核算，按式（8）-（11）。Vsr=Vsd+V𝑑𝑝ℎ+Vsf····································································(8)Vsd=λ×(Qsr𝜌)×C×10−4·······························································(9)V𝑑𝑝ℎ=∑Qsr×Ci×Pini=1×10−4·························································(10)Vsf=∑Qsr×Ki×Ptini=1×10−4+Qsr×Cc×Pc···········································(11)式中：Vsr——土壤保持总价值（万元·a-1）；ρ——土壤容重（t·m-3），参考附录B.2；λ——泥沙淤积系数，参考附录B.2；Vsd——减少泥沙淤积价值（万元·a-1）；V𝑑𝑝ℎ——减少面源污染价值量（万元·a-1）；Vsf——保持土壤肥力价值量（万元·a-1）；9C——单位水库清淤工程费用（元·m-3），可参考附录B.10；i——土壤中氮、磷等营养物质数量，i=1，2，…，n；Qsr——土壤保持量（t·a-1）；Ci——土壤中氮、磷等营养物质的纯含量（%），参考附录B.2；Pi——处理氮、磷废水的成本（元·t-1）；Ki——土壤保持的养分（氮、磷、钾）折算为尿素、过磷酸钙、氯化钾的折算系数；Pti——化肥（尿素、过磷酸钙和氯化钾）价格（元·t-1）；Cc——土壤中有机质平均含量（%）；Pc——有机质价格（元·t-1）。7.2.3洪水调蓄服务7.2.3.1实物量核算用可调蓄洪水量作为核算指标，核算按式（12）。𝐶𝑓𝑚=𝐶𝑓𝑐+𝐶𝑟𝑐+𝐶𝑚𝑐································································(12)式中：𝐶𝑓𝑚——洪水调蓄量（m3·a-1）；𝐶𝑓𝑐——森林、灌丛、草地洪水调蓄总量（m3·a-1），核算见附录B.3；𝐶𝑟𝑐——水库洪水调蓄量（m3·a-1），核算见附录B.3；𝐶𝑚𝑐——沼泽湿地洪水调蓄量（m3·a-1），核算见附录B.3。7.2.3.2价值量核算运用替代成本法核算区域内洪水调蓄价值，核算按式（13）。𝑉𝑓𝑚=（𝐶𝑓𝑚×𝐶𝑤𝑒）×10−4··························································(13)式中：𝑉𝑓𝑚——生态系统洪水调蓄价值（万元·a-1）；𝐶𝑤𝑒——水库单位库容的工程造价（元·m-3），参考附录B.10。7.2.4空气净化服务7.2.4.1实物量核算依据二氧化硫、氮氧化物等污染物浓度是否超过环境空气功能区质量标准，选用污染物排放量或空气净化能力，作为生态系统空气净化实物量。方法一：如果污染物浓度未超过环境空气功能区质量标准，大气污染物净化量为污染物排放量，核算按式（14）。𝑄𝑎𝑝=∑𝑄𝑖𝑛𝑖=1·······································································(14)式中：𝑄𝑎𝑝——大气污染物排放总量（kg·a-1）；i——大气污染物类别（二氧化硫、氮氧化物），（无量纲）；10n——大气污染物类别的数量（无量纲）；j——大气污染物类别（二氧化硫、氮氧化物），（无量纲）；𝑄𝑖——第i类大气污染物排放量（kg.a-1）方法二：如果污染物浓度超过环境空气功能区质量标准，大气污染物净化量为生态系统自净能力，核算按式（15）。𝑄𝑎𝑝=∑∑𝑄𝑖𝑗×𝐴𝑖𝑛𝑗=1𝑚𝑖=1·································································(15)式中：𝑄𝑎𝑝——生态系统空气净化能力（kg·a-1）；m——生态系统类型的数量（无量纲）；i——生态系统类型（无量纲）；n——大气污染物类别的数量（无量纲）；j——大气污染物类别（二氧化硫、氮氧化物）（无量纲）；𝑄𝑖𝑗——核算区内第i类生态系统第j种大气污染物的净化能力（kg·km-2·a-1），参考附录B.6；𝐴𝑖——第i类生态系统类型面积（km2）。7.2.4.2价值量核算采用替代成本法进行空气净化价值量核算，通过工业治理大气污染物成本评估生态系统空气净化价值，按式（16）。𝑉𝑎𝑝=∑𝑄𝑖𝑛𝑖=1×𝐶𝑖×10−4·····························································(16)式中：𝑉𝑎𝑝——生态系统空气净化价值（万元·a-1）；Ci——第i类大气污染物的治理成本（元·t-1），参考附录B.10。7.2.5水体净化服务7.2.5.1实物量核算采用水体污染物净化量作为核算指标，核算方法有以下两种，优先选用方法一。方法一：采用水体污染物净化量作为核算指标，按照GB3838—2002对水环境质量应控制项目的规定，选取COD、氨氮等指标，以水环境质量Ⅲ类为环境容量，核算水域湿地水体净化实物量，核算按式（17）。𝑄𝑤𝑝=（∑𝑄𝑤𝑖𝑛𝑖=1×𝑆𝑅）×10−9······················································(17)式中：𝑄𝑤𝑝——水体净化量（t·a-1）；𝑄𝑤𝑖——为水域湿地对第i类污染物的Ⅲ类水质标准（mg·L-1），参考附录B.7；i——污染物类别（无量纲）；𝑆𝑅——内陆湿地生态系统径流量（m3）。方法二：采用水域湿地生态系统对污染物（COD、氨氮）的自净能力计算水体净化实物量，核算按式（18）。11𝑄𝑤𝑝=∑𝑄𝑖𝑛𝑖=1×𝐴···································································(18)式中：𝑄𝑤𝑝——水体净化总量（t·a-1）；n——水体污染物类别的数量（无量纲）；i——污染物类别（无量纲）；𝑄𝑖——水域湿地对第i类水体污染物的年净化能力（t·km-2·a-1），可参考附录B.7；A——水域面积（km2）。7.2.5.2价值量核算采用替代成本法进行水体净化价值量计算，通过水体污染物治理成本进行核算，核算按式（19）。𝑉𝑤𝑝=∑𝑄𝑤𝑝𝑖𝑛𝑖=1×𝐶𝑖×10−4··························································(19)式中：𝑉𝑤𝑝——水体净化价值（万元·a-1）；𝑄𝑤𝑝𝑖——第i类水体污染物净化总量（t·a-1）；Ci——第i类污染物的处理成本（元·t-1），参考附录B.10。7.2.6固碳服务7.2.6.1实物量核算采用生态系统固定二氧化碳量作为核算指标，森林、灌丛和草地等生态系统植被固碳量采用净生态系统生产力法（NEP法），湿地与土壤固碳实物量采用固碳速率法核算。核算按式（20）-（23）。𝑄𝑇=𝑄𝐸+𝑄𝑆+𝑄𝑤··································································(20)𝑄𝐸=𝑀𝑐𝑜2𝑀𝐶×𝑁𝐸𝑃····································································(21)𝑄𝑤=𝑀𝑐𝑜2𝑀𝐶×∑𝐵𝑗×𝑊𝑗𝑚𝑗······························································(22)𝑄𝑆=𝑀𝑐𝑜2𝑀𝐶×∑𝐴𝑖×𝑆𝑖𝑛𝑖································································(23)式中：𝑄𝑇——陆地每年生态系统总固碳量（t·CO2·a-1）；𝑄𝐸——森林、灌丛和草地植被每年二氧化碳固定量（t·CO2·a-1）；𝑄𝑆——土壤每年二氧化碳固定量（t·CO2·a-1）；𝑄𝑤——湿地每年二氧化碳固定量（t·CO2·a-1）；𝑀𝑐𝑜2𝑀𝐶——C转化为CO2的系数，44/12；NEP——净生态系统生产力（t·C·a-1），核算见附录B.4；𝐵𝑗——j类湿地的面积（hm2）；𝑊𝑗——j类湿地的固碳速率（gC·m-2·a-1），参考附录B.4；𝐴𝑖——不同生态系统的面积（hm2）；12𝑆𝑖——为不同生态系统实测土壤固碳速率（t·CO2·hm-2·a），或参考附录B.4。7.2.6.2价值量核算生态系统固碳价值采用碳交易法进行核算，固定二氧化碳价值量按式（24）核算。𝑉𝑐𝑓=𝑄𝑇×𝐶𝑐×10−4································································(24)式中：𝑉𝑐𝑓——固碳价值（万元·a-1）；Cc——碳交易价格（元·t-1），参考附录B.10。7.2.7氧气释放服务7.2.7.1实物量核算氧气释放实物量核算按式（25）。𝑄𝑜𝑝=𝑁𝑃𝑃×1.19···································································(25)式中：𝑄𝑜𝑝——生态系统释氧量（t·a-1）；𝑁𝑃𝑃——为生态系统净初级生产力（t）。7.2.7.2价值量核算𝑉𝑜𝑝=𝑄𝑜𝑝×𝐶𝑐𝑜×10−4·······························································(26)式中：Vop——生态系统释放氧气价值（万元·a-1）；Cco——工业制氧成本（元·t-1），参考附录B.10。7.2.8局部气候调节服务7.2.8.1实物量核算局部气候调节实物量核算有两种方法，优先选用方法一，在无实测数据情况下采用方法二。方法一：采用实际测量生态系统内外温差吸收的大气热量作为生态系统气候调节的实物量，核算按式（27）。𝐸𝑐𝑟=∑∆𝑇𝑖×𝜌𝑐×𝑉𝑛𝑖=1·······························································(27)式中：𝐸𝑐𝑟——生态系统内外温差吸收的总能量（kW·h·a-1）；∆𝑇𝑖——第i天生态系统内外实测温差（℃）；𝜌𝑐——空气的比热容J/m3·℃；V——生态系统内空气的体积（m3）；n——年内日最高温超过26℃的总天数。13方法二：采用生态系统蒸散发过程中消耗的能量作为局部气候调节的实物量，核算按式（28）。𝐸𝑐𝑟=𝐸𝑝𝑝+𝐸𝑤𝑒·····································································(28)式中：𝐸𝑐𝑟——生态系统蒸腾蒸发消耗的总能量（kW·h·a-1）；𝐸𝑝𝑝——生态系统植被蒸腾消耗的能量（kW·h·a-1），核算见附录B.5；𝐸𝑤𝑒——生态系统水面蒸发消耗的能量（kW·h·a-1）,核算见附录B.5。7.2.8.2价值量核算𝑉𝑐𝑟=（𝐸𝑐𝑟×𝑝）/γ×10−4···························································(29)式中：𝑉𝑐𝑟——生态系统气候调节价值（万元·a-1）；p——电价（元/kW·h），参考附录B.10。γ——空调能效比，取3.25。7.2.9负氧离子释放服务7.2.9.1实物量核算按照GB/T38582-2020中负氧离子功能量核算方法，按式（30）。𝐺𝑛𝑜𝑖=5.256×1015×𝑄𝑛𝑜𝑖×𝐴×𝐻𝐿·····························································(30)式中：𝐺𝑛𝑜𝑖——森林中提供负氧离子个数（个·a-1）；𝑄𝑛𝑜𝑖——森林负氧离子浓度（个·cm-3）；A——森林面积（hm2）；H——实测森林树木高度（m）；L——负氧离子寿命（min）。7.2.9.2价值量核算价值量核算按式（31）。𝑉𝑛𝑜𝑖=5.256×1015×（𝑄𝑛𝑜𝑖−600）×𝐴×𝐻𝐿×𝐾𝑛𝑜𝑖×10−4···············································(31)式中：𝑉𝑛𝑜𝑖——森林提供负氧离子价值（万元）；𝐾𝑛𝑜𝑖——负氧离子生产费用（元·个-1），参考附录B.10。7.2.10海岸带防护服务7.2.10.1实物量核算根据实际调查获取的海防林表征陆地生态系统海岸带防护服务的实物量。147.2.10.2价值量核算海岸带防护价值量核算按式（32）。𝑉𝐷=𝐷×𝐶·········································································(32)式中：𝑉𝐷——海岸带防护服务价值（万元）；C——人工岸线建设成本与年均维护成本（万元·km-1）。7.2.11物种保育更新7.2.11.1实物量核算以动植物及其生境为主要保护对象的各类自然保护地、生态保护红线的面积为实物量单元，并把濒危动植物、特有动植物和古树名木的数量纳入物种保育实物量计算，体现生物多样性保护等级信息。核算按式（33）。𝐺𝑏𝑖𝑜=𝐴×（1+0.1×∑𝐸𝑚+0.1×∑𝐵𝑛+0.1×∑𝑂𝑟）𝑧𝑟=1𝑦𝑛=1𝑥𝑚=1································(33)式中：𝐺𝑏𝑖𝑜——物种保育的实物量（hm2）；𝐸𝑚——区域内物种m的濒危分值，核算见附录B.9；𝐵𝑛——区域内物种n的特有值，核算见附录B.9；𝑂𝑟——区域内物种r的古树年龄指数，核算见附录B.9；x——为计算濒危指数物种数量；y——为计算特有种指数物种数量；z——为计算古树年龄指数物种数量；A——自然保护区面积（hm2）。7.2.11.2价值量核算物种保育价值量核算按式（34）。𝑉𝑏𝑖𝑜=𝐺𝑏𝑖𝑜×𝑆𝑐×10−4·······························································(34)式中：𝑉𝑏𝑖𝑜——物种保育价值（万元·a-1）；𝑆𝑐——为单位面积物种保育价值（元·hm-2·a-1），根据Shannnon-Weiner指数进行单位面积物种保育价值确定，具体参考附录B.10。7.2.12陆地调节服务总价值调节服务价值总量核算按式（35）。𝐸𝑆𝑉=𝑉𝑤𝑟+𝑉𝑠𝑑+𝑉𝑑𝑝ℎ+𝑉𝑓𝑚+𝑉𝑎𝑝+𝑉𝑤𝑝+𝑉𝑐𝑓+𝑉𝑐𝑟+𝑉𝑛𝑜𝑖+𝑉𝐷+𝑉𝑏𝑖𝑜···················(35)式中：𝐸𝑆𝑉——调节服务价值总量（万元）。157.3文化服务7.3.1旅游康养7.3.1.1实物量核算采用区域内统计年鉴中接待游客总人数作为核算指标。𝑁𝑡=∑𝑁𝑡𝑖𝑛𝑖=1········································································(36)式中𝑁𝑡——游客总人数（人次）；𝑁𝑡𝑖——第𝑖个旅游景观的游览人数（人次）；𝑛——旅游景观个数，𝑖=1，2，…，n。7.3.1.2价值量核算方法一：采用区域统计年鉴自然景观区域内旅游总收入作为旅游康养价值量。方法二：采用支付意愿问卷调查法。旅游康养价值核算按式（37）-（39）。𝑉𝑐=∑×𝑁𝑡×𝑇𝐶𝑗𝑛𝑗=1·································································(37)𝑇𝐶𝑗=𝑇𝑗×𝑊𝑗+𝐶𝑗···································································(38)𝐶𝑗=𝐶𝑡𝑐,𝑗+𝐶𝑙𝑓,𝑗+𝐶𝑒𝑓,𝑗······························································(39)式中：𝑉𝑐——休闲旅游总价值（万元·a-1）；𝑁𝑡——从j地区到核算区受调查游客人次（万人次）；𝑇𝐶𝑗——来自j地区的每人次游客的平均旅游消费水平（元·人次-1）；𝑇𝑗——表示来自j地的游客用于旅途和核算旅游地点的平均时间（天·人-1）；𝑊𝑗——表示来自j地的游客的当地平均工资（元·人-1·天-1）；𝐶𝑗——表示来自j地的游客花费的平均直接旅行费用（元·人-1），其中包括游客从j地到核算区域的交通费用𝐶𝑡𝑐,𝑗（元/人）、食宿花费𝐶𝑙𝑓,𝑗（元·人-1）和门票费用𝐶𝑒𝑓,𝑗；𝑗——来核算地区旅游的游客所在区域个数，j=1,2,3，……，n。7.3.2文化服务总价值𝐸𝐶𝑉=𝑉𝑐×10−4····································································(40)式中：𝐸𝐶𝑉——文化服务价值总量（亿元·a-1）。8海洋和滨海湿地生态产品总值核算方法8.1物质产品供给服务8.1.1养殖生产8.1.1.1实物量核算16养殖生产的实物量应采用评估海域的主要类别养殖水产品的年产量进行评估，分鱼类、甲壳类、贝类、藻类、其它等五类。𝑌𝑓=∑𝑌𝑓𝑖𝑛𝑖=1········································································(41)式中：𝑌𝑓——养殖生产总产量（kg·a-1）；𝑌𝑓𝑖——𝑖类供给产品产量（kg·a-1）；𝑛——核算海域养殖生产产品的类别数。8.1.1.2价值量核算核算地区能获取本项产品增加值时，推荐采取增加值进行核算，没有增加值数据时，采用市场价格法进行评估。核算按式（42）。𝑀𝑆𝑀=∑(𝑄𝑆𝑀𝑖×𝑃𝑀𝑖)×10−4𝑛𝑖=1·······················································(42)式中：𝑀𝑆𝑀——养殖生产价值（万元·a-1）；𝑄𝑆𝑀𝑖—第𝑖类养殖水产品产量（kg·a-1），𝑖=1，2，3，4，5分别代表鱼类、甲壳类、贝类、藻类和其它；𝑃𝑀𝑖——第𝑖类养殖水产品的平均市场价格（元·kg-1）。8.1.2捕捞生产8.1.2.1实物量核算捕捞生产的实物量采用捕捞年产量进行评估，数据来自统计年鉴。海南省捕捞生产的统计包括远洋捕捞和近海域捕捞，需扣除远洋捕捞，只对近海域捕捞的实物量进行核算。如评估海域存在商业捕捞，则捕捞生产的实物量应采用捕捞年产量进行评估。如评估海域存在商业捕捞或者非商业捕捞活动，但是没有捕捞产量统计数据，捕捞生产的实物量应根据邻近海域同类功能区主要品种的捕捞量与资源量的比例推算。如缺少评估海域渔业资源现存量数据，可采用临近海域同类功能区单位面积海域渔业资源现存量数据推算。𝑌𝑐=∑𝑌𝑐𝑖𝑛𝑖=1········································································(43)式中：𝑌𝑐——捕捞生产总产量（kg·a-1）；𝑌𝑐𝑖——第𝑖类供给产品产量（kg·a-1）；𝑛——核算海域捕捞生产产品的类别数。8.1.2.2价值量核算核算地区能获取本项产品增加值时，推荐采取增加值进行核算，没有增加值数据时，采用市场价格法进行评估。计算公式为：𝑉𝑆𝐶=∑(𝑄𝑆𝐶𝑖×𝑃𝐶𝑖)×10−4𝑛𝑖=1·························································(44)式中：17𝑉𝑆𝐶——捕捞生产价值（万元·a-1）；𝑄𝑆𝐶𝑖——第𝑖类捕捞水产品的产量（kg/a），𝑖=1，2，3，4，5，6分别代表鱼类、甲壳类、贝类、藻类、头足类和其它；𝑃𝐶𝑖——第𝑖类捕捞水产品的平均市场价格（元·kg-1）。捕捞水产品的平均市场价格应采用评估海域临近海产品批发市场的同类海产品批发价格进行计算。应根据评估海域实际情况进行调整。8.1.3种质资源8.1.3.1实物量核算种质资源的实物量采用评估海域主要培育的水产品育苗的年产量进行评估，分海水鱼苗、虾类育苗、贝类育苗、藻类育苗和海参育苗等5类。8.1.3.2价值量核算种质资源的价值量采用市场价格法进行评估。𝑉𝐺𝑀=∑(𝑄𝐺𝑀𝑖×𝑃𝐺𝑀𝑖)×10−4𝑛𝑖=1······················································(45)式中：𝑉𝐺𝑀——种质资源价值（万元·a-1）：𝑄𝐺𝑀𝑖——第i类种质资源的产量（尾或个）；i=1,2,3,4,5分别代表海水鱼苗、虾类育苗、贝类育苗、藻类苗和海参育苗；𝑃𝐺𝑀𝑖——第i类种质资源产品的市场价格（元·尾-1或个-1）。8.1.4物质产品供给服务价值海洋生态系统物质产品供给服务总价值核算按式（46）。𝐸𝑃𝑉𝑆=𝑀𝑆𝑀+𝑉𝑆𝐶+𝑉𝐺𝑀·······························································(46)式中：𝐸𝑃𝑉𝑆——物质产品供给服务总价值（万元·a-1）。8.2调节服务8.2.1固碳服务8.2.1.1实物量核算选用固定二氧化碳量作为滨海湿地与海洋生态系统固碳服务实物量的评价指标。采用红树林、海草床、浮游植物、大型藻类和滤食性贝类计量方法核算海洋生态系统固碳服务。核算按式（47）。𝑄𝐶𝑂2=𝑄𝑚+𝑄𝑝+𝑄𝑎+𝑄𝑠····························································(47)式中：𝑄𝑚——红树林、海草床年固碳量（t·a-1）；𝑄𝑝——浮游植物的年固碳量（t·a-1）；𝑄𝑎——大型藻类的年固碳量（t·a-1）；𝑄𝑠——滤食性贝类年固碳量（t·a-1）。18——红树林、海草床等滨海湿地固碳实物量核算。红树林、海草床固碳实物量计算方法主要有三种：根据数据可获得性，建议优先选择生物量法，其次选择固碳速率法，最后选择净生态系统生产力法（NEP法）。方法一：固碳速率法，核算按式（48）。𝑄𝑚=∑𝑀𝑐𝑜2𝑀𝐶×𝐵𝑗×𝑊𝑗𝑚𝑗=1·····························································(48)式中：𝐵𝑗——红树林、海草床的面积（hm2）；𝑊𝑗——红树林、海草床的固碳速率（t·CO2·hm-2·a-1）。方法二：生物量法，核算按式（49）-（52）。𝑄𝑚=（𝑄𝑎𝑏+𝑄𝑢𝑛+𝑄𝑠）×𝑀𝐶𝑂2/𝑀𝐶····················································(49)𝑄𝑎𝑏=𝐺𝑎𝑏.𝑡∗𝐶𝐹𝑎𝑏.𝑡−𝐺𝑎𝑏.𝑡−1∗𝐶𝐹𝑎𝑏.𝑡−1··················································(50)𝑄𝑢𝑛=𝐺𝑢𝑛.𝑡∗𝐶𝐹𝑢𝑛.𝑡−𝐺𝑢𝑛.𝑡−1∗𝐶𝐹𝑢𝑛.𝑡−1··················································(51)𝑄𝑠=𝐺𝑠.𝑡∗𝐶𝐹𝑠.𝑡−𝐺𝑠.𝑡−1∗𝐶𝐹𝑠.𝑡−1························································(52)式中：𝑄𝑚——红树林、海草床年固碳量（t·a-1）；𝑄𝑎𝑏——红树林、海草床地上年固碳量（t·a-1）；𝑄𝑢𝑛——红树林、海草床地下年固碳量（t·a-1）；𝑄𝑠——红树林、海草床土壤年固碳量（t·a-1）；𝐺𝑎𝑏.𝑡和𝐺𝑎𝑏.𝑡−1——分别是核算年和核算前一年红树林、海草床的地上生物量（t·a-1），参考HJ1169监测获得；𝐶𝐹𝑎𝑏.𝑡和𝐶𝐹𝑎𝑏.𝑡−1——分别是核算年和核算前一年红树林、海草床的地上部分含碳量（%）；𝐺𝑢𝑛.𝑡和𝐺𝑢𝑛.𝑡−1——分别是核算年和核算前一年红树林、海草床的地下生物量（t·a-1），参考HJ1169监测获得；𝐶𝐹𝑢𝑛.𝑡和𝐶𝐹𝑢𝑛.𝑡−1——分别是核算年和核算前一年红树林、海草床的地上部分含碳量（%）；𝐺𝑠.𝑡和𝐺𝑠.𝑡−1——分别是核算年和核算前一年红树林、海草床的土壤质量（t·a-1）；𝐶𝐹𝑠.𝑡和𝐶𝐹𝑠.𝑡−1——分别是核算年和核算前一年红树林、海草床的土壤含碳量（%）。方法三：采用净生态系统生产力估算方法，核算按式（53）。𝑄𝑚=𝑀𝐶𝑂2/𝑀𝐶×𝑁𝐸𝑃·······························································(53)式中：𝑄𝑚——红树林、海草床等生态系统二氧化碳固定总量（t·CO2·a-1）；𝑀𝐶𝑂2/𝑀𝐶——CO2与C的分子量之比，44/12；𝑁𝐸𝑃——净生态系统生产力（t·C·a-1）。——海洋浮游植物、大型藻类、贝类等固碳实物量评估方法有两种，根据数据可获得性选取评估方法。如两种方法都适用，选取最能反映潜力的作为最终实物量。方法一：基于海洋浮游植物、大型藻类和贝类固定二氧化碳的原理计算，适用于所有海域评估，核算按式（54）-（57）。19𝑄𝑐=𝑄𝑝+𝑄𝑎+𝑄𝑠···································································(54)𝑄𝑝=3.67×𝑄𝑝𝑝×𝑆×365×10−3×𝑀𝐶𝑂2/𝑀𝐶··········································(55)𝑄𝑎=1.63×𝑄𝐴×𝑀𝐶𝑂2/𝑀𝐶···························································(56)𝑄𝑠=𝑃×𝑀𝐶𝑂2/𝑀𝐶(𝑅𝑠𝑡∙𝑊𝑐+𝑅𝑠∙𝑊𝑐𝑠)·················································(57)式中：𝑄𝑃𝑃——浮游植物的初级生产力（mg·m-2.d-1），参考附录B4.2；𝑆——评估海域的水域面积（km2）；𝑄𝐴——大型藻类的干重（t·a-1）；𝑃——贝类产量（t·a-1）；𝑅𝑠𝑡——软体组织干质量比例（%）；𝑊𝑐——软体组织含碳量（%）𝑅𝑠——贝壳干质量比例（%）；𝑊𝑐𝑠——贝壳含碳量（%）。方法二：基于海洋吸收大气二氧化碳的原理计算，适用于有海气界面二氧化碳通量监测数据的大面积海域评估，核算按式（58）。𝑄𝑐=𝑄𝐶𝑂2×𝑆·······································································(58)式中：𝑄𝑐——海洋生态系统固定二氧化碳总量（t·a-1）;𝑄𝐶𝑂2——单位面积海域每年吸收二氧化碳的量（t·CO2/（a·km2）），可参考GB/T28058，南海海域每年吸收二氧化碳的量4.76t·CO2·a-1·km-2）;𝑆——评估海域的水域面积（km2）。8.2.1.2价值量核算固碳价值量应采用市场价格法进行评估。核算按式（59）。𝑉𝐶𝑂2=𝑄𝐶𝑂2×𝑃𝐶𝑂2×10−4····························································(59)式中：𝑉𝐶𝑂2——固碳价值（万元·a-1）；𝑄𝐶𝑂2——固碳实物量（t·a-1）；𝑃𝐶𝑂2——二氧化碳排放权的市场交易价格（元·t-1），采用海南省碳交易市场价格或全国当年碳交易市场均价，见附录B.10。8.2.2氧气释放服务8.2.2.1实物量核算用二氧化碳固定量换算成氧气释放量作为评估指标。核算按式（60）。𝑄𝑚𝑜𝑝=𝑎×𝑄𝑐𝑜2·····································································(60)式中：20𝑄𝑚𝑜𝑝——近岸海域生态系统氧气生产量（t）；𝑄𝑐𝑜2——近岸海域生态系统中海岸带红树林、海草床等植被、浮游植物和大型藻类的固碳量（t）；𝑎——近岸海域生态系统中固碳量换算成氧气释放量的系数，取32/44。8.2.2.2价值量核算释氧价值量应采用市场价格法进行评估。核算按式（61）。𝑉𝑜𝑟=𝑄𝑜𝑟×𝐶𝑜×10−4································································(61)式中：𝑉𝑜𝑟——近岸海域生态系统氧气生产价值（万元·a-1）；𝑄𝑜𝑟——近岸海域生态系统氧气释放量（t·a-1）；𝐶𝑜——工业制氧成本（元/t），见附录B.10。8.2.3废弃物处理服务8.2.3.1实物量核算废弃物处理的实物量评估有三个方法可以选用：方法一：对于已知环境容量的海域，宜采用环境容量值进行评估，废弃物处理的实物量按COD、无机氮、活性磷酸盐的容纳量计算。依据GB3097，选取化学需氧量、无机氮、活性磷酸盐等为主要污染物不同水质标准浓度限值，按照《海南省近岸海域环境功能区划（2010年修编）》中划定的纳污海洋功能区的面积，进行废弃物处理实物量核算。核算按式（62）。𝑄𝑆𝑊𝑇=∑𝑄𝑖𝑗×3𝑗=1𝐴𝑗×𝐻𝑗×10−6·····················································(62)式中：𝑄𝑆𝑊𝑇——为维护目标水质要求的第i项污染物的环境容纳量（t·a-1）；𝑄𝑖𝑗——第i项污染物第j种水质的浓度限值（mg·L-1），参考附录B8；I——化学需氧量、无机氮、活性磷酸盐3项主要污染物；j——二类、三类和四类水质类型；𝐴𝑗——j类近岸海域纳污海洋功能区面积（km2），参考附录B8；𝐻𝑗——j类近岸海域纳污海洋功能区平均水深（m）。方法二：基于海洋生态系统对陆源污染物无机氮、活性磷酸盐的净化及排海化学需氧量的处理能力，选用各类水体污染物净化量，作为海洋生态系统废弃物处理实物量的评价指标。核算按式（63）。𝑄𝑆𝑊𝑇=𝑄𝑡𝐶𝑂2×16106+𝑄𝑡𝐶𝑂2×1106+𝑃𝐶𝑂𝐷················································(63)式中：𝑄𝑡𝐶𝑂2——海洋生态系统固碳量（t·a-1）；16/106——根据浮游植物对营养盐的吸收总体上遵循Redfield比值（C:N:P=106:16:1）的规律取得；1/106——根据浮游植物对营养盐的吸收总体上遵循Redfield比值（C:N:P=106:16:1）的规律取得；𝑃𝐶𝑂𝐷——排海化学需氧量（t·a-1）。方法三：采用排海废弃物的数量进行评估，排废弃物主要考虑废水、COD、无机氮、活性磷酸盐等。核算按式（64）。21𝑄𝑆𝑊𝑇=𝑄𝑊𝑊−𝑄𝑊𝑊×𝑤×20%+𝑄𝑊𝑇················································(64)式中：𝑄𝑆𝑊𝑇——废弃物处理（考虑排海废水）的实物量（t·a-1）;𝑄𝑊𝑊——工业和生活废水产生量（t·a-1）;w——工业和生活废水所含污染物的质量分数（%）;𝑄𝑊𝑊×𝑤——表示工业和生活废水所含的污染物总量（t·a-1）;𝑄𝑊𝑊×𝑤×20%——表示废水通过河流、沟渠入海过程中滞留在途中的污染物量（t·a-1），按20%的滞留率计算。8.2.3.2价值量核算废弃物处理的价值量应采用替代成本法进行评估。核算按式（65）。𝑉𝑆𝑊=∑𝑄𝑆𝑊𝑇𝑚𝑖=1×𝑃𝑖×10−4··························································(65)式中：𝑉𝑆𝑊——废弃物处理的价值量（万元·a-1）；𝑄𝑆𝑊𝑇——废弃物处理的实物量（t·a-1）；𝑃𝑖——第𝑖项污染物的处理成本（元·t-1），见附录B.10。8.2.4局部气候调节8.2.4.1实物量核算选用生态系统蒸腾蒸发消耗的能量作为局部气候调节实物量评估指标，也可使用实际测量的生态系统内外温差转化为生态系统吸收的大气热量评估实物量，优先选用实际测量方法。方法同7.2.8。8.2.4.2价值量核算方法同7.2.8。8.2.5海岸带防护8.2.5.1实物量核算运用红树林、珊瑚礁等生态系统防护或替代海堤等防护工程的长度估算海岸带防护实物量。核算按式（66）。𝐷𝑐𝑙=∑𝐷𝑐𝑙𝑖𝑛𝑖=1·······································································(66)式中：𝐷𝑐𝑙——海岸带生态系统防护的海岸线总长度（km·a-1）；𝐷𝑐𝑙𝑖——i类生态系统防护的海岸线长度（km·a-1），i=1，2，…，n；𝑛——核算海域海岸带防护生态系统类型的数量。8.2.5.2价值量核算通过替代成本法核算海岸防护的经济价值，通过统计生态系统防护范围内由于减轻风暴潮造成海水淹没海岸带的经济损失和减轻海浪造成损毁岸堤的经济损失得到。核算按式（67）：22𝑉𝑐𝑝=∑𝐷𝑐𝑙𝑖×𝑊𝑐𝑙𝑖×𝑉𝑐𝑙𝑎×𝐷𝑟𝑐𝑙+𝐶𝑐𝑙×10−4𝑛𝑖=1·········································(67)式中：𝑉𝑐𝑝——海岸带防护总价值（万元·a-1）；𝐷𝑐𝑙𝑖——𝑖类生态系统防护的海岸线长度（km），𝑖=1，2，…，n；𝑊𝑐𝑙𝑖——𝑖类生态系统海岸防护的权重，即单位长度𝑖类生态系统海岸防护的生态效益与单位长度人工岸线防护的生态效益的比值，无量纲；𝑉𝑐𝑙𝑎——单位长度人工岸线的造价成本（元·km-1·a-1），包括减轻风暴潮造成海水淹没海岸带的经济损失、减轻海浪造成损毁岸堤的经济损失等，人工岸线寿命按照核算地域所涉工程建设标准或设计规范计算，默认值取50年；𝐶𝑐𝑙——单位长度人工岸线年维护成本（元·km-1·a-1）；𝐷𝑟𝑐𝑙——人工岸线工程年折旧率；𝑛——核算海域海岸带防护生态系统类型的数量。8.2.6物种保育更新8.2.6.1实物量核算方法一：采用对评估海域内分布的国家级、省级的自然保护区、海洋特别保护区和水产种质资源保护区数量来进行评估，采用区域保护区面积进行核算。核算按式（68）𝐺𝑏𝑖𝑜𝑝=𝑆···········································································(68)式中：Gbiop——物种保育实物量（hm2）；𝑆——保护区面积（hm2）。方法二：如区域未划定自然保护区，物种保育更新实物量核算包含采用评估海域内分布的海洋保护物种数（国家级、省级）以及在当地有重要价值（科学的、文化的、宗教的、经济的）的海洋物种数来进行评估。核算按式（69）：𝐺𝑏𝑖𝑜𝑝=𝐴×（1+0.1×∑𝐸𝑚𝑥𝑚=1+0.1×∑𝐵𝑛𝑦𝑛=1）·······································(69)式中：𝐺𝑏𝑖𝑜𝑝——物种保育实物量；𝐸𝑚——区域内保护物种的濒危分值，见附录B.9；𝐵𝑛——区域内物种n的特有值，见附录B.9；𝑥——计算濒危指数物种数量；𝑦——计算特有种指数物种数量；𝐴——群落面积（hm2）。8.2.6.2价值量核算物种保育更新价值量推荐采用保护区保护价值法。核算按式（70）：𝑉𝑏𝑖𝑜𝑝=𝐺𝑏𝑖𝑜𝑝×𝑆𝑐×10−4·····························································(70)式中：𝑉𝑏𝑖𝑜𝑝——物种保育价值（万元·a-1）;23𝐺𝑏𝑖𝑜𝑝——物种保育实物量；𝑆𝑐——为单位面积保育成本（元·hm-2·a-1），见附录B.10。8.2.7调节服务总价值调节服务价值总量核算按式（71）。𝐸𝑆𝑉𝑆=(𝑉𝐶𝑂2+𝑉𝑜𝑟+𝑉𝑆𝑊+𝑉𝑡𝑡+𝑉𝑐𝑝+𝑉𝑏𝑖𝑜𝑝)×10−4····································(71)式中：𝐸𝑆𝑉𝑆——调节服务价值总量（亿元）。8.3文化服务8.3.1旅游康养8.3.1.1实物量核算旅游康养服务评估主要评估海域以自然海洋景观为主体的海洋旅游景区。休闲娱乐的实物量采用海洋旅游景区的年旅游人数评估。若旅游人数很少，可不进行该项评估。方法同7.3.1。8.3.1.2价值量核算方法同7.3.1。8.3.2文化服务总价值𝐸𝐶𝑉𝑆=𝑉𝐶×10−4···································································(72)式中：𝐸𝐶𝑉𝑆——海洋和滨海湿地文化服务价值总量（亿元·a-1）。9附则9.1质量控制9.1.1核算管理核算工作应按照第7章与第8章的规定开展陆地与海洋生态产品的实物量和价值量核算，并实行自检、互检、技术负责人检查和专家验收“三检一验”制度，检验内容包括每个核算科目的数据来源、实物量和价值量核算结果。9.1.2数据要求生态产品总值核算所需数据应符合附录C中的规定。数据质量应保证可获得数据的正确性、完整性、逻辑性、一致性。9.2成果汇总9.2.1成果类型24核算成果可汇总为下列类型：——图件成果：1)生态系统分类标准产品：在自然资源调查数据或遥感专题信息自动提取成果的基础上，经过转换形成满足生态产品总值核算的生态系统分类标准产品。2)专题图件：根据生态产品总值核算技术规范，制作相应的实物量和价值量专题图件。——表册成果：包括打印签字版或盖章扫描版的数据来源确认表、问卷调查表、核算结果汇总表等。——文本成果：区域生态产品总值核算报告。9.2.2数据格式除文字和表册外，其他空间数据需转换为GEOTIFF格式（栅格数据格式）或SHAPEFILE格式（矢量数据格式）。9.2.3成果安全生态产品总值核算过程中涉及国家秘密的资料和数据，应按保密规定进行管理，不得发生失密、泄密问题。9.3成果应用9.3.1决策应用核算成果可应用于绿色发展财政奖补、国土空间管控、生态系统保护修复、环境治理评估、自然资源资产负债表编制等领域，调节服务产品核算结果可作为制定生态保护补偿标准、生态环境损害赔偿标准等制度的依据。9.3.2规划应用核算成果反映出的生态产品总值结构特征和变化趋势可应用于促进绿色高质量发展的战略规划、行动计划、生态保护修复规划、实施方案的编制，评估生态产品价值影响，确保生态产品保值增值。9.3.3项目应用核算成果可应用于国家级、省级等重大生态保护修复工程建设效益评价。核算成果可应用于土地出让、生态产品开发经营与市场化交易，推动生态惠民富民。9.3.4考核应用核算成果可应用于市县高质量发展综合考核、党政领导班子绩效考核、领导干部自然资源资产离任审计等工作中，推动完善绿色政绩考核制度。25AA附录A（规范性）生态产品清单表A.1为生态产品清单。表A.1生态产品清单一级指标二级指标说明物质产品供给服务农业产品从农业生态系统中获得的初级产品，如谷物、豆类、油料、糖料；热带水果、蔬菜、薯类、坚果、饮料和香料作物、胡椒；中药材等。林业产品林木产品、林产品以及与森林资源相关的初级产品，如木材、竹材、槟榔、橡胶、林药、茶叶、菌类、花卉、苗木等。畜牧业产品用圈养、散养或者两者结合的方式，饲养禽畜以取得动物产品，如牛、羊、猪、家禽年底出栏数；奶类；禽蛋等。渔业产品人类利用水域中生物的物质转化功能，通过捕捞、养殖等方式取得的水产品。（淡）水资源优质的生态环境为本地或下游提供高品质水资源。种质资源生态系统提供农、林、牧、渔业等种子产品。生态能源生态系统中的生物物质及其所含的能力，如水能、沼气、秸秆等。调节服务水源涵养生态系统通过其结构和过程拦截滞蓄降水，增强土壤下渗，有效涵养土壤水分和补充地下水、调节河川流量，增加可利用水资源的功能。土壤保持生态系统通过其结构与过程降低雨水的侵蚀能力，减少土壤流失的功能。洪水调蓄森林、灌丛、草地、湿地等生态系统通过调节暴雨径流、削减洪峰，减轻洪水危害的功能。陆地固碳陆地生态系统通过光合作用将大气中的CO2转化为碳水化合物，并以有机碳的形式固定在植物体内或土壤中，降低大气中CO2浓度的功能。海洋与滨海湿地固碳海洋与滨海湿地生态系统通过红树林、浮游植物、大型藻类等光合作用固定CO2，和滤食性贝类通过摄食活动直接吸收海水中的碳酸氢根，并最终以有机碳或其他形式输送到深层海水，并埋藏到海底，降低大气中CO2浓度的功能。氧气释放生态系统通过光合作用释放氧气，减缓大气中CO2浓度升高，维持大气氧气浓度稳定的功能。空气净化生态系统吸收、阻滤和分解大气中的污染物，如SO2、NOx、粉尘等，有效净化空气，改善大气环境的生态效应。水体净化水环境通过一系列物理和生化过程对进入其中的污染物进行吸附、转化以及生物吸收等，使水体得到净化的生态效应。废弃物处理海洋生态系统对人类产生的各种排海污染物如COD、氨氮、总磷、总氮、石油烃等的降解、吸收和转换服务。局部气候调节生态系统通过植被蒸腾作用和水面蒸发过程吸收能力，调节温湿度的功能。负氧离子释放生态系统通过物理和生物过程，释放负氧离子，产生空气净化和健康疗效的功能。海岸带防护生态系统降低海浪，避免或减小海堤或海岸侵蚀的功能。26一级指标二级指标说明物种保育更新生态系统为动植物物种提供生存与繁衍的场所，从而对其起到保育作用的功能。文化服务旅游康养为人类提供美学价值、灵感、教育价值等非物质惠益的自然景观，其承载的价值对社会具有重大的意义。BB附录B（资料性）生态产品实物量和价值量核算关键参数B.1水源涵养服务核算主要参数B.1.1地表径流有实测径流系数时，按式（B.1）核算地表径流量：𝑄=𝑃×𝑅········································································(B.1)无实测径流系数时，采用SCS模型进行地表径流量计算。SCS概念模型如下：𝑄={(𝑃−0.2𝑆)2𝑃+0.8𝑆,𝑃≥0.2𝑆0,𝑃<0.2𝑆····························································(B.2)𝑆=25400𝐶𝑁−254···································································(B.3)式中：𝑄——日径流量（径流深度，mm）；𝑃——日降雨量（mm）；𝑅——径流系数（无量纲）；𝑆——土壤的潜在最大蓄水载荷（mm）；𝐶𝑁——径流曲线系数（无量纲），取值见表B.1.2。𝑆𝐶𝑆概念模型细化到每场降雨和每种地类，公式如下：𝑄𝑖,𝑗={(𝑃𝑖−0.2𝑆𝑖,𝑗)2𝑃𝑖+0.8𝑆𝑖,𝑗,𝑃≥0.2𝑆𝑖,𝑗0,𝑃<0.2𝑆𝑖,𝑗······················································(B.4)𝑆𝑖,𝑗=25400𝐶𝑁𝑖,𝑗−254·································································(B.5)式中：𝑖和𝑗分别表示降雨的频次和生态系统类型。在年尺度上，对每种生态系统类型所有降雨场次产生的径流量进行累计：𝑄𝑗=∑Q𝑖,𝑗𝑛𝑖=1·····································································(B.6)27式中：𝑄𝑗——地类𝑗每年的径流总量。𝐶𝑁是反映降雨前下垫面地表特征的一个综合指标，𝐶𝑁值在0-100之间，𝐶𝑁值越大，𝑆值越小，越容易产生净流；反之，则相反。𝐶𝑁需要根据土壤前期湿润程度（Antecedentmoisturecondition，简称AMC）、坡度、土壤类型和土地利用现状等综合推算。𝐶𝑁值的确定有以下步骤：a)雨量站控制单元的划分。基于雨量站空间位置信息，在ArcGIS软件下划分泰森多边形，确定每个雨量站的控制单元。b)土壤水文类型的划分。以土壤类型图为输入数据，根据表B.1划分土壤水文类型（A、B、C、D四类）。表B.1土壤水文类型划分标准类别土壤水文性质稳定下渗率（mm/h）A厚层沙、厚层黄土、团粒化粉砂土7.26-11.73B薄层黄土、砂壤土3.81-7.26C粘壤土、薄层砂壤土、粘质含量高的土壤1.27-3.81D吸水后显著膨胀的土壤、塑性的粘土、某些盐渍土0-1.27c)降雨场次的选取。为了充分表达雨洪的概念，对降雨强度为大雨及其以上（日降雨量大于25mm）的降雨场次进行径流模拟。以此为标准，逐雨量站进行降雨场次的选取。d)生态系统类型的划分。根据土地利用类型将生态系统进行划分归类，尽可能的进行细分，如低地雨林、山地雨林、季雨林、橡胶林、经济果木林、用材林、灌木林、湿地、农田、城市用地、裸地等。e)CN参数的确定。根据土壤水文类型和生态系统类型，按照上文叙述中确定的赋值标准，逐斑块确定CN参数。可以参考的CN值如表B.2所示。表B.2海南省部分生态系统类型的CN值参考参数生态系统类型CN_ACN_BCN_CCN_D低地雨林30304148山地雨林30304148季雨林30304148橡胶林46667783经济果木林32587279用材林46667783灌木林49687984草地49697984湿地98989898农田67788589不透水地面91919191裸地7786919428B.1.2实际蒸散发实际蒸散发可用如下公式表示：𝐴𝐸𝑇𝑖,𝑚=min(𝑃𝐸𝑇𝑖,𝑚;𝑃𝑖,𝑚−𝑄𝐹𝑖,𝑚+𝛼𝑚𝛽𝑖𝐿𝑠𝑢𝑚.𝑎𝑣𝑎𝑖𝑙,𝑖)·································(B.7)𝑃𝐸𝑇𝑖,𝑚=𝐾𝑐,𝑖,𝑚×𝐸𝑇0,𝑖,𝑚···························································(B.8)𝐸𝑇0=0.408∆(𝑅𝑛−𝐺)+𝛾900𝑇+273𝜇2(𝑒𝑠−𝑒𝑎)∆+𝛾(1+0.34𝑢2)····················································(B.9)𝐿𝑠𝑢𝑚.𝑎𝑣𝑎𝑖𝑙,𝑖=∑𝑝𝑖𝑗𝑗∈{排水到像元𝑖的邻近栅格}×(𝐿𝑎𝑣𝑎𝑖𝑙,𝑗+𝐿𝑠𝑢𝑚.𝑎𝑣𝑎𝑖𝑙,𝑗)······························(B.10)𝐿𝑎𝑣𝑎𝑖𝑙,𝑖=𝑚𝑖𝑛(𝜕𝐿𝑖,𝐿𝑖)····························································(B.11)式中：𝐾𝑐,𝑖,𝑚——像元𝑖上生态系统m月的作物系数；通过查阅文献获得，可参考表B.3；𝐸𝑇0,𝑖,𝑚——像元𝑖上生态系统的m月的参考蒸散发，单位mm；𝑅𝑛——净辐射，单位MJ/（m2·d）；𝐺——土壤热通量，单位MJ/（m2·d）；𝛾——干湿表常数，取4；𝑇——日均温度（℃）；通过气象监测站获取；𝑢2——2m风速（m/s）；通过气象监测站获取；𝑒𝑠和𝑒𝑎——分别为饱和水汽压和实际水汽压（kPa）；∆——饱和水汽压曲线斜率（kPa/℃）；𝛼𝑚——第m月可用的上游年可用补给的部分；𝛽𝑖——可用于下游蒸散的上游补给的部分；𝜕——可用于下游栅格的栅格重新补给部分；𝑅𝑛、𝑒𝑠、𝑒𝑎和Δ——具体计算过程参考手册“FAOIrrigationandDrainagePaper，No.56，CropEvapotranspiration”。表B.3海南省不同生态系统类型不同月份作物蒸散发参考值生态系统类型Kc_1Kc_2Kc_3Kc_4Kc_5Kc_6Kc_7Kc_8Kc_9Kc_10Kc_11Kc_12低地雨林0.8370.8910.8330.8870.8810.8560.8460.8130.8050.7830.7980.899山地雨林0.8370.8910.8330.8870.8810.8560.8460.8130.8050.7830.7980.899季雨林0.8370.8910.8330.8870.8810.8560.8460.8130.8050.7830.7980.899灌木林0.6640.6580.6820.8440.7880.5610.4020.340.3790.4110.6110.59草地0.6880.6760.7110.8350.8450.7040.5850.4880.5030.5450.6170.67湿地1.031.0281.0311.0991.091.0350.9720.9240.9350.9541.0121.016经济果0.7070.7230.7110.8490.8610.7570.6060.4850.4990.5480.6750.70629木林橡胶林0.950.950.950.950.950.950.950.950.950.950.950.95农田0.6680.6750.6780.8360.8370.7070.5430.4280.4560.510.650.665城镇0.760.7280.7470.9240.8860.7480.5910.5010.5240.5610.730.704裸地0.30.30.30.30.30.30.30.30.30.30.30.3B.2土壤保持服务核算主要参数B.2.1降雨侵蚀力因子降雨侵蚀力因子R核算方法如下：基于月降雨量的计算模型，核算按式（B.12）：𝑅=∑(−1.5527+0.1792𝑃𝑖)12𝑖=1×17.02············································(B.12)式中：𝑅——降雨侵蚀力值（MJ·mm·hm-2·h-1·a-1）；𝑖——月份；𝑃𝑖——月降雨量（mm）。B.2.2土壤可蚀性因子土壤可蚀性因子K核算按式（B.13）：𝐾=(−0.01383+0.51515×𝐾𝐸𝑃𝐼𝐶)×0.1317······································(B.13)𝐾𝐸𝑃𝐼𝐶={0.2+0.3𝑒[−0.0256𝑚𝑠×(1−𝑚𝑠𝑖𝑙𝑡)]}×[𝑚𝑠𝑖𝑙𝑡𝑚𝑐+𝑚𝑠𝑖𝑙𝑡]0.3×{1−0.25𝑜𝑟𝑔𝐶[𝑜𝑟𝑔𝐶+𝑒(3.72−2.95×𝑜𝑟𝑔𝐶)]}×{1−0.7×(1−𝑚𝑠100)(1−𝑚𝑠100)+𝑒[−5.51+22.9(1−𝑚𝑠100)]}·············································································(B.14)式中：𝐾——土壤可蚀性值（t·hm2·h·hm-2·MJ-1·mm-1）；𝑜𝑟𝑔𝐶——土壤有机质含量百分比（%）；𝑚𝑠——粘粒（<0.002mm）含量百分比（%）；𝑚𝑠𝑖𝑙𝑡——粉粒（0.002mm～0.05mm）含量百分比（%）；𝑚𝑐——砂粒（0.05mm～2mm）含量百分比（%）。各土壤类型质地与有机质含量可参考表B.4。表B.4海南省土壤质地和有机质表土壤类型砂砾（%）粉粒（%）粘粒（%）有机质（%）砖红壤85.1113.061.831.905赤红壤48.37537.8813.7451.905黄壤39.8939.7820.331.905水稻土35.7745.3618.871.559燥红土38.9944.8816.130.69330紫色土58.1532.649.211.732火山灰土43.9942.8813.132.771石质土61.5531.217.244.330新积土69.7126.144.151.039盐土31.4854.0614.464.560注：数据来源于《海南岛土壤地球化学手册》。B.2.3坡长和坡度因子坡长因子L核算按式（B.15）、式（B.16）和式（B.17），坡度因子S的核算按式（B.18）：𝐿=(𝜆22.13)𝑚····································································(B.15)𝑚=𝛽1+𝛽·······································································(B.16)𝛽=(𝑠𝑖𝑛𝜃0.0896)(3.0×𝑠𝑖𝑛𝜃0.8+0.56)·····························································(B.17)𝑆={10.8×𝑠𝑖𝑛(𝜃)+0.03,𝜃<9%16.8×𝑠𝑖𝑛(𝜃)−0.5,9%≤𝜃≤18%21.91×𝑠𝑖𝑛(𝜃)−0.96,𝜃>18%·····································(B.18)式中：𝐿——坡长因子；𝜆——坡长（m）；𝑚——无量纲常数，取决于坡度百分比值；𝑆——坡度因子；𝜃——坡度（弧度）。B.2.4盖度和管理因子盖度和管理因子C是指在一定的覆盖度和管理措施下，一定面积土地上的土壤流失量与采取连续清耕、休闲处理的相同面积土地上的流失量的比值，为无量纲数，介于0-1之间。采用植被覆盖度计算得到，核算按式（B.19）。植被覆盖度f基于植被指数NDVI数据计算得到，核算按式（B.20）。𝐶={1𝑓=00.6508−0.3436𝑙𝑔𝑓0<𝑓≤78.3%𝟎𝑓>78.3%········································(B.19)𝑓=(𝑁𝐷𝑉𝐼−𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙)(𝑁𝐷𝑉𝐼𝑚𝑎𝑥−𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙)····························································(B.20)式中：𝐶——植被覆盖因子；𝑓——植被覆盖度（%）；𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙——纯裸土象元的𝑁𝐷𝑉𝐼值；31𝑁𝐷𝑉𝐼𝑚𝑎𝑥——纯植被象元的𝑁𝐷𝑉𝐼值。B.2.5水土保持措施因子水土保持措施因子P反映水土保持措施对土壤侵蚀的抑制作用，是指采取水土保持措施的土壤侵蚀量与采用顺坡种植时土壤侵蚀量的比值。P值的取值范围为0-1，值越小，表示水土保持能力越好，土壤越不易被侵蚀。在不能通过实验测定的情况下，可以通过土地利用类型的典型P值替代。海南省水田多位于地势平坦区域，旱地种植多采用等高耕作，因此水田的P值取0.15，旱地取值0.35；乔灌果园部分采取挖排水沟、防水壕等简易水土保持措施，因此园地的P值取0.6，其余土地利用类型基本没有采取水土保持措施。不同土地利用类型取值见表B.5。表B.5不同土地利用类型的水土保持措施因子P值土地利用类型水田旱地园地林地草地裸地不透水地面水体P值0.150.350.611100B.2.6泥沙淤积系数泥沙淤积系数取值24%。B.2.7土壤容重与土壤氮、磷、有机质含量根据《海南岛生态地球化学调查》与农业部门土壤调查成果，剔除最高值、最低值与异常值后，整理得到海南省11种主要土壤类型的土壤容重与养分元素均值，可作为土壤保持功能核算的参考。表B.6海南省土壤容重与土壤中氮、磷含量参考值序号土壤类型平均容重(g/cm3)有机质(%)全N（％）全磷（％）全钾（％）1滨海砂土1.670.5200.0370.0211.2452砖红壤1.521.9050.0900.0402.0743赤红壤1.531.9050.0950.0383.0704黄壤1.441.9050.0990.0362.9045水稻土1.511.5590.0750.0431.0796燥红土1.480.6930.0460.0222.2407紫色土1.061.7320.0840.0312.2408火山灰土1.072.7710.1380.1470.3329石质土1.134.3300.2070.1700.49810新积土1.041.0390.0630.0452.32311盐土1.214.5600.1200.3102.240B.3洪水调蓄服务核算主要参数B.3.1林灌草洪水调蓄能力核算32Cfc=∑(𝑃𝑖−𝑅𝑓𝑖)×𝐴𝑖𝑛𝑖=1×1000·················································(B.21)式中：𝑃𝑖——研究区域内暴雨降雨量量（mm）；𝑅𝑓𝑖——研究区域内暴雨径流量（mm）。暴雨径流量指暴雨期（一般以日降雨50mm为标准）的地表径流量。核算方式见附录B.1。B.3.2水库洪水调蓄能力根据水库防洪库容与水库泄洪次数构建模型，由此估算水库总防洪库容：𝐶𝑟𝑐=𝐶𝑡×0.29·································································(B.22)式中：𝐶𝑟𝑐——水库总防洪库容（亿m3）；0.29——库容转换为防洪库容的系数，取值参考国家发改委国家统计局印发的《生态产品总值核算规范（试行）》中东部平原区防洪库容系数；𝐶𝑡——水库库容（亿m3）。B.3.3沼泽湿地洪水调蓄能力湿地土壤具有特殊的水文物理性质，湿地土壤的草根层和泥炭层孔隙度达72%-93%，饱和持水量达830%-1030%，在汛期大量水资源被储存于沼泽湿地土壤中或者以表面水的形式保存在湿地中，湿地就像一个巨大的天然蓄水库容纳水分，直接减少了洪水径流，同时湿地植被可拦截径流减缓洪水流速，因此削减和滞后了洪峰，减低了下游洪峰的水位，并使之缓慢下泄，这样有效地对洪水形成缓冲滞纳作用，度过汛期，进而减少洪灾发生。通过构建沼泽土壤蓄水量和地表滞水量模型计算沼泽湿地洪水调蓄能力。𝐶𝑚𝑐=𝑆×ℎ×𝜌×(𝐹−𝐸)×106×𝜌𝑤+𝑆×𝐻×106·····························(B.23)式中：𝐶𝑚𝑐——沼泽湿地调蓄洪水的能力（万m³）；𝐻——沼泽湿地地表滞水高度（m﹒a-1）；ℎ——沼泽湿地土壤蓄水深度（m﹒a-1）；𝜌——沼泽湿地土壤容重（t﹒m-3）；𝐹——沼泽湿地土壤饱和含水率（无量纲）；𝐸——沼泽湿地洪水淹没前的土壤自然含水率（无量纲）；𝜌𝑤——水的密度（t﹒m-3）；𝑆——沼泽湿总面积（hm2）。表B.7海南省沼泽湿地洪水调蓄能力关键参数参考值土壤含水率差值土壤蓄水深度地表滞水高度0.2540680.4m0.3mB.4固定二氧化碳服务核算参数B.4.1陆地生态系统固碳核算参数33（1）植被固碳采用净生态系统生产力NEP核算方法如下：NEP根据海南省的NEP和NPP的转换系数计算得到，核算按式（B.24）。𝑁𝐸𝑃=𝛼×𝑁𝑃𝑃×𝑀𝐶6/𝑀𝐶6𝐻10𝑂5··················································(B.24)式中：𝑁𝐸𝑃——净生态系统生产力（t·C﹒a-1）；𝛼——𝑁𝐸𝑃和𝑁𝑃𝑃的转换系数；𝑁𝑃𝑃——净初级生产力（t·干物质﹒a-1）；𝑀𝐶6/𝑀𝐶6𝐻10𝑂5=72/162——干物质转化为C的系数。根据中国植被区划，海南省北部属于热带东部偏湿性季雨林区域，南部属于中热带季雨林、湿润雨林区域。根据海南尖峰岭森林生态系统国家野外科学观测研究站多年连续观测研究基础数据（2000-2020年）、海南东寨港红树林生态系统定位研究站数据以及相关文献整理得出海南省典型植被类型NEP-NPP转换系数参考值，见表B.8。表B.8海南典型植被类型NEP-NPP转换系数植被类型低地雨林山地雨林季雨林橡胶林用材林经济果木林灌木林草地α系数0.19090.35190.17050.25030.21930.14390.13790.3069（2）固碳速率法各类陆地生态系统面积、化学氮肥NF、复合肥施用量CF、作物j在当年的产量CYj等数据来自自然资源、林业、农业、水利和统计等部门的遥感数据、统计数据、实地调查或相关文献数据。森林及灌丛固碳速率、草地土壤固碳速率、湿地固碳速率来自于实测数据或者参考表B.9和表B.10。土壤固碳速率由海南尖峰岭森林生态系统国家野外科学观测研究站多年连续观测研究基础数据（2000-2020年）整理而得，见表B.9。表B.9海南省土壤固碳速率（（t·C·hm-2·a-1））生态系统类型低地雨林山地雨林季雨林橡胶林经济果木林用材林灌木林草地土壤固碳速率1.21.552.14.50.063.550.1140.02湿地生态系统固碳速率参考国家发展改革委国家统计局关于印发《生态产品总值核算规范（试行）》中给的参考值以及海南东寨港红树林生态系统定位研究站数据，见表B.10。表B.10海南省湿地生态系统固碳速率湿地类型沿海滩涂内陆滩涂红树林地沼泽地湖泊固碳速率（gC·m-2·a-1）235.6267.11444.2767.1156.67（3）农田土壤固碳速率𝐶𝑆𝐶𝑆=(𝐵𝑆𝑆+𝑆𝐶𝑆𝑅𝑁+𝑃𝑅×𝑆𝐶𝑆𝑅𝑆)×𝑆𝐶········································(B.25)34式中：𝐶𝑆𝐶𝑆——农田固碳量（tC/a）；𝐵𝑆𝑆——无固碳措施条件下的农田土壤固碳速率（tC·hm-2·a-1）；𝑆𝐶𝑆𝑅𝑁——施用化学氮肥和复合肥的农田土壤固碳速率（tC·hm-2·a-1）；𝑆𝐶𝑆𝑅𝑆——秸秆全部还田的农田土壤固碳速率（tC·hm-2·a-1）；𝑃𝑅——农田秸秆还田推广施行率（%）；𝑆𝐶——农田面积（hm2）。——无固碳措施条件下的农田土壤固碳速率：𝐵𝑆𝑆=𝑁𝑆𝐶×𝐵𝐷×𝐻×0.1······················································(B.26)式中：𝑁𝑆𝐶——无化学肥料和有机肥料施用的情况下，农田土壤有机碳的变化（g·kg-1·a-1）；𝐵𝐷——土壤容重（g·cm-3）；𝐻——土壤厚度（cm）。——施用化学氮肥、复合肥和秸秆还田的土壤固碳速率：施用化学氮肥的固碳速率：𝑆𝐶𝑆𝑅𝑁=1.5339×𝑇𝑁𝐹−266.7··················································(B.27)𝑇𝑁𝐹=(𝑁𝐹+𝐶𝐹×0.3)/𝑆𝑃······················································(B.28)式中：𝑇𝑁𝐹——单位面积耕地化学氮肥、复合肥总施用量（kgN·hm-2·a-1）；𝑁𝐹——化学氮肥施用量（t）；𝐶𝐹——复合肥施用量（t）；𝑆𝑃——耕地面积（hm2）。秸秆还田的固碳速率：𝑆𝐶𝑆𝑅𝑆=43.548×𝑆+375.1·····················································(B.29)𝑆=∑𝐶𝑌j𝑛𝑗=1×𝑆𝐺𝑅𝑗/𝑆𝑃··························································(B.30)式中：𝑆——单位耕地面积秸秆还田量（t·hm-2·a-1）；𝐶𝑌j——作物j在当年的产量（t）；𝑆𝑃——耕地面积（hm2）；𝑆𝐺𝑅𝑗——作物j的草谷比；𝑗——作物类别，j=1，2，…，n；𝑛——作物的种类。农田秸秆还田推广施行率、土壤容重、无化学肥料和有机肥料施用的情况下农田土壤有机碳的变化、土壤厚度、作物j的草谷比来自于实测数据或者参考文献，推荐参考国家发展改革委国家统计局关于印发《生态产品总值核算规范（试行）》的取值。B.4.2浮游植物初级生产力核算参数海洋浮游植物初级生产力的估算采用叶绿素a法，按联合国教科文组织（UNESCO）推荐的下列公式估算：35𝑃=𝐶ℎ𝑙𝑎∙𝑄∙𝐷∙𝐸2···································································(B.31)式中：𝑃——初级生产力（mg•C·m-2•d-1）；𝐶ℎ𝑙𝑎——真光层内平均叶绿素a含量（mg·m-3）；𝑄——不同层次同化指数算术平均值，取3.71；𝐷——昼长时间（h），全省取12小时；𝐸——真光层深度（m），取透明度（m）×2.71。B.4.3藻类与双壳贝类固碳核算参数参考《HY/T0305-2021养殖大型藻类和双壳贝类碳汇计量方法碳储量变化法》。表B.11海南常见养殖藻类成藻含碳率种类江篱麒麟菜石花菜坛紫菜其他藻类碳含量推荐取值28.40%30.36%26.37%41.96%30.36%表B.12海南常见双壳贝类含碳率种类软体部含碳率推荐值贝壳含碳率推荐值牡蛎44.90%11.52%蚶45.86%11.29%贻贝45.98%12.68%扇贝43.87%11.44%蛤44.90%11.52%其他42.82%11.45%B.5局部气候调节服务核算主要参数B.5.1植被蒸腾消耗能量森林、灌丛、草地生态系统植被蒸腾消耗的能量。𝑄𝑃𝑃=∑𝐸𝑃𝑃𝑖×𝑆𝑖×𝐷×106𝑛𝑖=1/(3600×𝛾)·······································(B.32)式中：𝑄𝑃𝑃——生态系统植被蒸腾消耗的能量（kW·h·a-1）；𝐸𝑃𝑃𝑖——第i生态系统类型单位面积蒸腾消耗热量（kJ·m-2·d-1）；𝑆𝑖——第i种生态系统类型面积（km2）；𝛾——空调能效比；𝐷——空调开放天数（天）；𝑖——研究区不同生态系统类型（森林、灌丛、草地）。36数据优先采用实际测量数据。如无实测数据，可参考相关文献获得。单位面积森林、灌丛、草地蒸腾吸热量参考值见下表B.13。表B.13生态系统单位面积蒸腾吸热量（kJm-2d-1）参考值生态系统类型森林灌丛草地蒸腾吸热2837.271300.95969.83B.5.2水面蒸发消耗能量水面蒸发降温增湿消耗的能量。𝑄𝑤𝑒=𝐸𝑤×𝑞×1033600×𝛾+𝐸𝑤ℎ×𝑦·················································(B.33)式中：𝑄𝑤𝑒——水面蒸发消耗能量（kW·h）；𝐸𝑤——开放空调降温期间水面蒸发量（m3）；𝐸𝑤ℎ——开放加湿器增湿期间水面蒸发量（m3）；𝑞——挥发潜热，即蒸发1克水所需要的热量（J·g-1），取2443.9；𝛾——空调能效比；𝑦——加湿器将1m3水转化为蒸汽的耗电量（kW·h·m-3），取120，仅计算湿度小于45%时的增湿功能。小型蒸发器观测的蒸发量与自然水体蒸发量的比值，海南为0.685。数据来源于国家发展改革委国家统计局关于印发《生态产品总值核算规范（试行）》中给的水面蒸发折算系数。B.6空气净化服务核算主要参数根据文献李意德等著的《海南生态公益林生态服务功能价值评估研究》，单位生态系统面积主要污染物净化量见表B.14。表B.14海南省单位生态系统面积主要污染物净化量参考值（t·km-2·a-1）低地雨林山地雨林季雨林橡胶林经济果木林用材林灌木林草地SO2净化能力4.5005.2205.0007.0746.5157.0748.9641.13氮氧化物净化能力0.850.911560.850.744460.634110.80570.6000.061滞尘能力1.5001.9551.5001.9120.5051.9120.0870.03B.7水体净化服务核算主要参数根据《地表水环境质量标准》（GB3838-2002）中对水环境质量应控制项目的规定，地表水各类水环境质量污染物浓度限值见表B.15。表B.15地表水污染物浓度限值单位：mg·L-137污染物Ⅰ类Ⅱ类Ⅲ类Ⅳ类Ⅴ类化学需氧量≤1515203040氨氮（NH3-N）≤0.150.51.01.52.0总氮≤0.20.51.01.52.0总磷≤0.02（湖、库0.01）0.1（湖、库0.025）0.2（湖、库0.05）0.3（湖、库0.1）0.4（湖、库0.2）表B.16单位面积湿地对各类水体污染物的净化量COD的净化能力（t·km-2·a-1）氨氮净化能力（t·km-2·a-1）总磷净化能力（t·km-2·a-1）湿地110.4308.5640.540B.8海洋废弃物处理服务主要参数根据《海南省近岸海域环境功能区划（2010年修编）》，海南省近岸海域共划分了106个功能区，各环境功能区面积与水质目标见表B.17。表B.17海南省近岸海域环境功能区统计表（km2）市县名称一类水质目标二类水质目标三类水质目标四类水质目标个数面积个数面积个数面积个数面积昌江164.17249.3746.46澄迈256250.63儋州284.614191.16345.52256.8东方259.78396.87243.47210.1海口4209.3823.1326.53乐东170.59113.5312.5300临高182.33173.6237.1400陵水3183.64332.6610.6400琼海2137.4511.6326.2700三亚3166.617140.85521.3812.7万宁5281.5524924.500文昌6307.794110.02324.8712.7儋州临高0013.6000038东方、昌江00121.230000海口、澄迈00143.730000三亚、乐东198.09000000南海海区224000000000西沙群岛1196.47000000合计3025733.05391192.6327195.9110129.46表B.18海水水质标准（mg/L）分类第二类第三类第四类化学需氧量345无机氮0.30.40.5活性磷酸盐0.030.030.045B.9物种保育更新服务核算参数表B.19物种濒危指数体系濒危指数濒危等级4极危3濒危2易危1近危注：物种种类参见《IUCN受胁物种红色名录》《国家重点保护动物名录》《国家重点保护野生植物名录》；濒危等级参考《IUCN受胁物种红色名录》，以及《国家重点保护动物名录》《国家重点保护野生植物名录》保护等级。将保护等级为“国家一级”对应“极危”、国家二级对应“濒危”，结合IUCN等级，采用就高不就低原则确定物种濒危等级。表B.20特有种指数体系特有种指数分布范围4仅限于范围不大的山峰或特殊的自然地理环境下分布3仅限于某些较大的自然地理环境下分布的类群，如仅分布于较大的海岛（岛屿）、高原、若干条山脉等2易仅限于某个大陆分布的分布群1至少在2个大陆都有分布的分布群0世界广布的分类群表B.21古树年龄指数体系39指数等级古树年龄1100-299年2300-499年3≥500年表B.22Shannon-Wiener指数等级划分及其价值量等级Shannon-Wiener指数单价/（元.hm-2·a-1）I指数≥666760II5≤指数≤653410III4≤指数≤540060IV3≤指数≤426700V2≤指数≤313350VI1≤指数≤26680B.10价格参数表B.23价值量核算参数参考表类别参数单位2020年参考价定价依据水资源生活取用水元/m³0.15《海南省物价局省财政厅省水务厅关于调整水资源费征收标准的通知》中2015年水资源费征收标准城镇公共供水价格地表水与地下水收费均值。可根据省最新水资源费征收标准调整。工业取用水元/m³0.20《海南省物价局省财政厅省水务厅关于调整水资源费征收标准的通知》中2015年水资源费征收标准工业和其他行业取用水与工商业取用水均值。可根据省最新水资源费征收标准调整。农业取用水元/m³0.085《海南省物价局省财政厅省水务厅关于调整水资源费征收标准的通知》中2015年水资源费征收标准农业取用水价格的地表水与地下水均值。可根据省最新水资源费征收标准调整。水源涵养、洪水调蓄水库单位库容造价元/m³6.11参考LY/T1721—2008中第5章给出的推荐使用价格。建议从当地水利局获取该参数。若无法获取，可根据PPI指数，调整到核算年份。土壤保持水库清淤工程费用元/m345参考LY/T1721—2008中第5章给出的推荐使用价格。建议从当地水利局获取该参数。若无法获取，可根据PPI指数，调整到核算年份。磷酸二铵价格元/t3895海南市场调研价格氯化钾价格元/t3300海南市场调研价格有机肥价格元/t2180海南市场调研价格40类别参数单位2020年参考价定价依据含氮废水处理成本元/t3500《海南省人民代表大会常务委员会关于海南省大气污染物和水污染物环境保护税适用税额的决定》含磷废水处理成本元/t11200空气净化净化二氧化硫元/t2526.3净化氮氧化物元/t2526.3净化工业粉尘元/t600水体净化净化COD元/t2800净化氨氮元/t3500净化总磷元/t11200固碳碳交易价格元/t28全国碳交易年平均成交价（核算年）释放氧气工业制氧价格元/t8002020年海南市场调研工业制氧生产均价调节气候普通居民合表用户电价元/kW·h0.63海南电网有限责任公司提供普通居民合表用户电价均价负氧离子负氧离子生产费用元/1018个9.5根据负氧离子发生器价格结合平均电价折算估算注：如需多年比较，可采用海南省多年PPI数据将定价参数折算到其他核算年的可比价。41附录C（资料性）生态产品总值核算基础数据要求C.1为生态产品总值核算基础数据要求。表C.1价值量核算参数参考表服务类别核算科目核算指标单位数据来源数据要求供给服务农业产品谷物谷物产量、单价、增加值统计部门、农业农村部门核算年度、乡镇尺度薯类薯类油料油料豆类豆类糖料糖料蔬菜蔬菜水果热带水果…...…...（根据核算样地实际产品更新）林业产品竹木采伐木材统计部门、林业部门竹材林产品槟榔橡胶林药茶叶菌类苗木…...…...（根据核算样地实际产品更新）畜牧业产品畜禽产量牛肉统计部门、农业农村部门羊肉猪肉禽肉猪出栏数牛出栏数羊出栏数42服务类别核算科目核算指标单位数据来源数据要求家禽出栏数…...…...（根据核算样地实际产品更新）奶类牛奶……禽蛋禽蛋小计渔业产品水产品淡水养殖海水养殖捕捞……种质资源种子产品农产品种子苗木水产种苗…...（根据核算样地实际产品更新）水资源生活用水量统计部门、水务部门工业用水量灌溉用水量小计生态能源生物质能源沼气统计部门、农业农村部门秸秆…...水能水能发电量统计部门、发改部门…...…...43服务类别核算科目核算指标单位数据来源数据要求（根据核算样地实际产品更新）小计调节服务调节服务土地利用数据（矢量图）-自然资源和规划部门核算年度、乡镇尺度森林资源清查数据-海岸带面积hm2海岸线数据（含岸线属性）km海防林面积hm2海防林长度km单位长度海岸防护工程建设成本近五年均值单位长度海岸防护工程年维护成本海洋保护动植物种类类每10年提供一次、空间分布土壤容重t/m3每10年提供一次土壤中氮的纯含量%土壤中磷的纯含量%年降雨量mm气象部门核算年度、乡镇尺度年蒸发量mm年均夏季气温℃年26℃气温天数天多年均降雨量mm每年提供20年平均、乡镇尺度多年均蒸发量mm多年月均温度℃多年月均湿度%水库单位库容建设成本元/m3水务部门近五年均值水库单位库容年维护成本元/m3水资源价格元/m3核算年度、市县尺度水库总库容万m3核算年度、分水库、分站点水库防洪库容万m3河流水文站点日径流量m3/s44服务类别核算科目核算指标单位数据来源数据要求河流水文站点泥沙含量kg/m3大中型水库出、入境流量m3/s生物量亿t林业部门核算年度、小班尺度小班面积hm2郁闭度%平均树高m古树名木棵核算年度、市县尺度造林成本元/亩濒危物种、特有物种数量种、只/个核算年度、保护区尺度自然保护区空间范围自然保护区保护成本万元/年各类污染物排放量元/t生态环境部门核算年度、市县尺度COD治理成本元/t二氧化硫治理成本元/t氮氧化物治理成本元/t工业粉尘治理成本元/t氨氮治理成本元/t入海河流污染物排放量t/a近岸海域浮游植物初级生产力mg/m2•d近岸海域表层海水叶绿素a含量mg•C/（m2•d）电价元/kW.h发改委、物价局核算年度、市县尺度单位长度海岸防护工程造价元/m氮肥、磷肥、钾肥、有机肥价格元/t负氧离子浓度个/cm3文化和旅游部门核算年度、站点数据文化服务旅游康养旅游人数万人/年核算年度、景区尺度旅游收入万元/年旅游景区清单-旅游调查样表-提供一次45表C.2生态产品总值核算涉及的空间专题图制作要求项目要求坐标系统采用2000国家大地坐标系。图面要素利用提取获得的应用专题产品作为主要图面要素，叠加省界、地区界、县界等行政区界要素，并根据实际需要添加各类名称注记，例如城市名称、行政区域注记、图名、核算时段、注记、比例尺等内容。图件产品主要图件包括要素详图和各个专题产品详图；其他图件方面可根据实际需要制作。46附录D（资料性）生态产品总值核算报告编写大纲XXX地区生态产品总值核算报告前言介绍核算背景、目的意义、任务来源等。A.1区域概况介绍核算区域地理范围、自然环境状况、经济社会状况、生态环境保护状况等基本情况。A.2核算目标与原则介绍核算目标、核算原则、核算依据、核算基准年。A.3核算方法与数据介绍主要核算思路、方法、数据来源与数据处理方法。A.4生态产品总值实物量介绍生态产品总值实物量的计算过程与结果，包括物质产品实物量、调节服务功能量和文化服务功能量。A.5生态产品总值价值量介绍生态产品总值价值量的计算过程与结果，包括物质产品价值量、调节服务价值量和文化服务价值量。A.6生态产品总值空间格局及变化驱动因素A.7结论与建议介绍生态产品总值核算的结论，提出相关政策建议。A.8附件包括生态产品总值核算过程中相关的技术资料及附表、附图等。47附录E（资料性）生态系统划分生态系统类型划分主要以自然资源部门第三次全国国土调查分类土地利用调查结果为基础，参考《国土空间调查、规划、用途管制分类指南（试行）》，结合森林资源清查以及实地植被调查划分。见表E.1。表E.1生态系统划分一级生态系统类型土地利用类型二级生态系统类型森林生态系统乔木林地低地雨林山地雨林季雨林用材林灌木林地灌木林竹林地用材林草地生态系统草地草地内陆湿地生态系统河流水面内陆湿地生态系统水库水面湖泊水面坑塘水面内陆滩涂水工建筑用地沟渠沼泽地森林沼泽灌丛沼泽海洋生态系统珊瑚礁滨海湿地生态系统海草床红树林地盐田沿海滩涂近海近海海岛海岛农田生态系统果园（干果园、油料园）经济果木林橡胶园橡胶林果园（干果园、油料园之外的果园）农田生态系统茶园其他园地水田水浇地旱地48表E.1生态系统划分（续）一级生态系统类型土地利用类型二级生态系统类型农田生态系统设施农用地农田生态系统田坎城市生态系统（建成区）商业服务业用地城市生态系统工矿用地住宅用地公共管理与公共服务用地特殊用地城市水体交通运输用地其他生态系统盐碱地其他生态系统沙地裸土地裸盐石砾地注：根据海岛覆被划分生态系统类型。49参考文献[1]GB3095—2012环境空气质量标准[2]GB3097海水水质标准[3]GB/T12343国家基本比例尺地图编绘规范[4]GB/T13923基础地理信息要素分类与代码[5]GB/T28592降水量等级[6]GB17378.7-2007海洋监测规范第7部分：近海污染生态调查和生物监测[7]CH/T9005基础地理信息数据库基本规定[8]LY∕T2899湿地生态系统服务评估规范[9]LY/T1812林地分类[10]LY/T2586空气负（氧）离子浓度观测技术规范[11]NY/T2997草地分类[12]DB33/T2274生态系统生产总值（GEP）核算技术规范陆域生态系统[13]DB45/T1230-2015红树林湿地生态系统固碳能力评估技术规程[14]生态环境部.地生态系统生产总值（GEP）核算技术指南（2020）[15]国家发展和改革委员会国家统计局.《生态产品总值核算规范（试行）》[M].人民出版社,2022[16]国家林业和草原局农业农村部.《国家重点保护野生动物名录（2021）》[17]国家林业和草原局农业农村部.《国家重点保护野生植物名录（2020）》[18]环境经济核算体系——生态系统核算（SystemofEnvironmentalAccounting-EcosystemAccounting，SEEAEA）[19]李意德，杨众养，陈德祥，等.海南生态公益林生态服务功能及其价值研究[M].中国林业出版社,2015.[20]何玉生，张固成、薛桂澄，等.海南岛土壤地球化学手册[M].中国地质大学出版社,2021.