CN115689395B - 一种生态系统调节服务价值核算方法 - Google Patents
一种生态系统调节服务价值核算方法 Download PDFInfo
- Publication number
- CN115689395B CN115689395B CN202211706450.3A CN202211706450A CN115689395B CN 115689395 B CN115689395 B CN 115689395B CN 202211706450 A CN202211706450 A CN 202211706450A CN 115689395 B CN115689395 B CN 115689395B
- Authority
- CN
- China
- Prior art keywords
- service
- unit
- value
- formula
- soil
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 120
- 230000033228 biological regulation Effects 0.000 title claims abstract description 48
- 239000002689 soil Substances 0.000 claims abstract description 92
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 90
- 238000000746 purification Methods 0.000 claims abstract description 35
- 241000607479 Yersinia pestis Species 0.000 claims abstract description 30
- 239000004576 sand Substances 0.000 claims abstract description 30
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 28
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 28
- 239000001301 oxygen Substances 0.000 claims abstract description 28
- 230000001105 regulatory effect Effects 0.000 claims abstract description 19
- 238000006467 substitution reaction Methods 0.000 claims abstract description 19
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 18
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 18
- 201000010099 disease Diseases 0.000 claims abstract description 17
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims abstract description 17
- 230000002265 prevention Effects 0.000 claims abstract description 11
- 238000012423 maintenance Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 57
- 230000003628 erosive effect Effects 0.000 claims description 52
- 239000003344 environmental pollutant Substances 0.000 claims description 25
- 231100000719 pollutant Toxicity 0.000 claims description 25
- 239000003337 fertilizer Substances 0.000 claims description 23
- 229910052698 phosphorus Inorganic materials 0.000 claims description 18
- 238000001704 evaporation Methods 0.000 claims description 17
- 229910052700 potassium Inorganic materials 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 16
- 229910052757 nitrogen Inorganic materials 0.000 claims description 15
- WCUXLLCKKVVCTQ-UHFFFAOYSA-M Potassium chloride Chemical compound [Cl-].[K+] WCUXLLCKKVVCTQ-UHFFFAOYSA-M 0.000 claims description 14
- 230000008020 evaporation Effects 0.000 claims description 14
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 12
- 235000015097 nutrients Nutrition 0.000 claims description 12
- 230000005068 transpiration Effects 0.000 claims description 12
- 241000238631 Hexapoda Species 0.000 claims description 11
- 230000005611 electricity Effects 0.000 claims description 10
- 238000004162 soil erosion Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 8
- 238000012407 engineering method Methods 0.000 claims description 8
- 230000014759 maintenance of location Effects 0.000 claims description 8
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 claims description 7
- 239000004202 carbamide Substances 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 7
- 239000001103 potassium chloride Substances 0.000 claims description 7
- 235000011164 potassium chloride Nutrition 0.000 claims description 7
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 6
- 239000001569 carbon dioxide Substances 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- YYRMJZQKEFZXMX-UHFFFAOYSA-N calcium;phosphoric acid Chemical compound [Ca+2].OP(O)(O)=O.OP(O)(O)=O YYRMJZQKEFZXMX-UHFFFAOYSA-N 0.000 claims description 5
- 230000001276 controlling effect Effects 0.000 claims description 5
- 230000035558 fertility Effects 0.000 claims description 5
- 239000002426 superphosphate Substances 0.000 claims description 5
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000011835 investigation Methods 0.000 claims description 4
- 238000003973 irrigation Methods 0.000 claims description 4
- 230000002262 irrigation Effects 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 239000002420 orchard Substances 0.000 claims description 4
- 239000013049 sediment Substances 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 239000004927 clay Substances 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 claims description 3
- GVVPGTZRZFNKDS-JXMROGBWSA-N geranyl diphosphate Chemical compound CC(C)=CCC\C(C)=C\CO[P@](O)(=O)OP(O)(O)=O GVVPGTZRZFNKDS-JXMROGBWSA-N 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000003746 surface roughness Effects 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 229910000019 calcium carbonate Inorganic materials 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000012216 screening Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 102000003712 Complement factor B Human genes 0.000 claims 1
- 108090000056 Complement factor B Proteins 0.000 claims 1
- 239000004016 soil organic matter Substances 0.000 claims 1
- 241000196324 Embryophyta Species 0.000 description 11
- 238000011156 evaluation Methods 0.000 description 10
- 238000011160 research Methods 0.000 description 9
- 239000000047 product Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 3
- 239000011574 phosphorus Substances 0.000 description 3
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 2
- YYRMJZQKEFZXMX-UHFFFAOYSA-L calcium bis(dihydrogenphosphate) Chemical compound [Ca+2].OP(O)([O-])=O.OP(O)([O-])=O YYRMJZQKEFZXMX-UHFFFAOYSA-L 0.000 description 2
- 229910000389 calcium phosphate Inorganic materials 0.000 description 2
- 235000013877 carbamide Nutrition 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 2
- 235000019691 monocalcium phosphate Nutrition 0.000 description 2
- 230000029553 photosynthesis Effects 0.000 description 2
- 238000010672 photosynthesis Methods 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- MGWGWNFMUOTEHG-UHFFFAOYSA-N 4-(3,5-dimethylphenyl)-1,3-thiazol-2-amine Chemical compound CC1=CC(C)=CC(C=2N=C(N)SC=2)=C1 MGWGWNFMUOTEHG-UHFFFAOYSA-N 0.000 description 1
- 244000025254 Cannabis sativa Species 0.000 description 1
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 1
- 241001272720 Medialuna californiensis Species 0.000 description 1
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004887 air purification Methods 0.000 description 1
- 239000003513 alkali Substances 0.000 description 1
- 229910021529 ammonia Inorganic materials 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003181 biological factor Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 229910001882 dioxygen Inorganic materials 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 235000008216 herbs Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000008235 industrial water Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000011591 potassium Substances 0.000 description 1
- 239000008213 purified water Substances 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011155 quantitative monitoring Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 239000003403 water pollutant Substances 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Cultivation Of Plants (AREA)
Abstract
本发明公开了一种生态系统调节服务价值核算方法，涉及生态系统服务价值技术领域，包括以下步骤：收集相应生态系统调节服务的基础数据，生态系统调节服务包括水源涵养服务、土壤保持服务、洪水调蓄服务、防风固沙服务、固碳释氧服务、大气净化服务、水质净化服务、气候调节服务和病虫害防治服务；根据基础数据的数据类型采用不同的数据预处理方法获得相应生态系统调节服务的实物量数据，再分别对生态系统调节服务中的各个服务的价值进行核算；对每个服务价值的核算结果进行求和。本发明通过采用各种价值化方法，如市场价值法、替代成本法、防护费用法等，实现了各类生态系统调节服务的价值化。
Description
技术领域
本发明属于生态系统服务价值技术领域，尤其涉及一种生态系统调节服务价值核算方法。
背景技术
生态系统服务是指通过生态系统的结构、过程和功能直接或间接得到的生命支持产品和服务。目前，国内外就生态系统服务价值的评估方法开展了大量的研究工作，但尚未形成一套统一的评估体系。已有研究在生态系统服务价值指标体系的设置上，一般采用生态系统产品供给、调节服务、文化服务三大类一级指标，区别主要在细化的二级指标上。生态系统服务价值的核算方法大致可以分为两类，即基于单位服务功能价格的方法（功能价值法）和基于单位面积价值当量因子的方法（当量因子法）。功能价值法首先分别计算不同生态系统服务功能的实物量，如土壤保持量、水源涵养量、固碳释氧量等，实物量通过水平衡方程、土壤流失方程、风力侵蚀模型、植物净化模型等可靠的自然-生态模型计算获得；然后基于实物量核算结果，采用市场法、影子工程法等方法将各类服务功能价值化。当量因子法是在区分不同种类生态系统服务功能的基础上，基于可量化的标准构建不同类型生态系统服务功能价值当量，然后结合生态系统的分布面积进行评估；该方法不需计算实物量，直接计算价值量，已经比较成熟。
通过对比功能价值法和当量因子法，可以发现功能价值法的主要优点在于可以分别采用专门方法计算各类生态系统服务实物量，如土壤保持功能的核算首先采用通用土壤流失方程计算土壤保持量，而后采用替代成本法计算土壤保持功能价值量，核算准确度相对较高；但功能价值法所需数据较为繁杂，一般只适用于中小尺度的区域生态系统服务价值核算。当量因子法在构建标准的当量因子表后，只需使用研究区域不同类型土地面积乘以对应的服务价值当量，即可得到相应的生态系统服务价值，优点在于所需数据简单，适用于大区域的生态系统服务价值核算；但当量因子法假设同种土地利用类型的单位面积生态系统服务价值相同，在核算准确性上低于功能价值法。
公开号为CN109544021A的中国专利公开了一种区域水生态服务价值综合评价方法，包括：将区域水生态服务功能分为提供产品、调节、支持和文化服务功能；确定各功能的评价指标，构建区域水生态服务价值评价指标体系；根据指标体系，计算各功能价值；根据计算价值，计算区域水生态服务功能的总价值；其中，提供产品功能的评价指标包括居民生活用水、产业用水、渔业产品、水电蓄能和水源地温；调节功能的评价指标包括地表水资源调蓄、地下水资源调蓄与补给、水质净化功能、气候调节功能、洪水灌溉和净化空气；支持功能的评价指标包括固碳、释氧、营养物质循环、生物多样性保护、生活质量改善和预防地面沉降；文化服务功能的评价指标包括旅游休闲娱乐、水景观功能和水文化传承。但是该专利并未考虑到病虫害防治、防风固沙、洪水调蓄、土壤保持等方面的服务价值。
因此，需要设计一种考虑到多方面的服务价值且准确度较高的核算方法。
发明内容
本发明针对现有技术存在的问题，提供了一种生态系统调节服务价值核算方法。
为实现上述目的，本发明采用的技术方案如下：
一种生态系统调节服务价值核算方法，包括以下步骤：
步骤S1：收集相应生态系统调节服务的基础数据，生态系统调节服务包括水源涵养服务、土壤保持服务、洪水调蓄服务、防风固沙服务、固碳释氧服务、大气净化服务、水质净化服务、气候调节服务和病虫害防治服务；
步骤S2：根据基础数据的数据类型采用不同的数据预处理方法获得相应生态系统调节服务的实物量数据，再分别对生态系统调节服务中的各个服务的价值进行核算；
步骤S3：对每个服务价值的核算结果进行求和。
基于上述技术方案，更进一步地，所述水源涵养服务价值的核算方法，包括以下步骤：
采用替代工程法从水量平衡的角度计算水源涵养服务价值，计算公式为：
公式一：；
公式二：；
公式三：；
式中，为生态系统涵养水源的价值，单位为元/年；/>为涵养水源服务的功能量，单位为m3；/>为水库单位库容的建设成本，单位为元/m3；/>为降雨量，单位为mm；为地表径流量，单位为mm；/>为蒸散发量，单位为mm；/>为不同土地利用类型的面积，单位为m2；/>为不同土地利用类型的平均地表径流系数。
基于上述技术方案，更进一步地，所述土壤保持服务价值的核算方法，包括以下步骤：
采用替代成本法进行计算，计算公式为：
公式四：；
公式五：；
公式六：；
公式七：；
式中，为生态系统土壤保持服务价值，单位为元/年；/>为保肥价值，单位为元/年；/>为减少泥沙淤积价值；/>为土壤保持量，单位为吨/年；/>为土壤中N、P、K营养元素含量，取表层土N、P、K含量的平均值分别为0.17%，0.08%，2.17%；/>为化肥的价格，单位为元，N、P、K分别转换成尿素、过磷酸钙、氯化钾三种肥料；/>为化肥中i种营养元素的比重。
基于上述技术方案，更进一步地，所述洪水调蓄服务价值的核算方法，包括以下步骤：
按照水库总库容的35%计算水库洪水调蓄功能，按照替代工程法计算价值量，计算公式如下：
公式八：；
公式九：；
式中，为水库调蓄洪水的价值，单位为元/年；/>为水库调蓄洪水能力，单位为m3；/>为水库总库容，单位为m3；/>为人工的洪水调蓄价格，单位为元/(m3·年)。
基于上述技术方案，更进一步地，所述防风固沙服务价值的核算方法，包括以下步骤：
采用修正风蚀方程计算风蚀量，并利用替代成本法，将治理这些沙化土壤的成本和减少风蚀土壤肥力作为防风固沙功能的价值，计算公式为：
公式十：；
式中，为防风固沙价值，单位为元/年；/>为固沙量，单位为吨/年；/>为砂砾堆积密度，单位为吨/m3；/>为土壤沙化标准覆盖沙厚度，单位为m；/>为治沙工程的平均成本，单位为元/m2；/>为土壤中N、P、K营养元素含量（%），N、P、K分别转换成尿素、过磷酸钙、氯化钾三种肥料。
基于上述技术方案，更进一步地，所述固碳释氧服务价值的核算方法，包括以下步骤：
通过工业制氧成本法计算生态系统释氧价值，计算公式为：
公式十一：；
公式十二：；
公式十三：；
式中，为生态系统固碳释氧价值，单位为元/年；/>为生态系统固碳量，单位为吨/年；/>为生态系统释放氧气量，单位为吨/年；/>为二氧化碳价格，单位为元/吨；/>为氧气价格，单位为元/吨；NPP为植被净初级生产力，单位为吨/年。
基于上述技术方案，更进一步地，所述大气净化服务价值的核算方法，包括以下步骤：
采用替代成本法，以各类大气污染物治理成本计算价值量，计算公式为：
公式十四：；
公式十五：；
式中，为大气净化价值，单位为元/年；/>为大气污染物净化量，单位为吨；Qvj为第v 类生态系统第 j 种大气污染物的单位面积净化量，v为生态系统类型，无量纲，j为大气污染物类别，无量纲；Av为第v类生态系统类型面积，单位为km2，cj为不同种类大气污染物处理成本，单位为元/吨。
基于上述技术方案，更进一步地，所述水质净化服务价值的核算方法，包括以下步骤：
采用替代成本法，根据达到同样净化效果所需的人工净化成本计算水质净化价值，计算公式为：
公式十六：；
公式十七：；
式中，为湿地生态系统水质净化价值，单位为元/年；/>为治理i种污染物的成本，单位为元/吨；/>为湿地生态系统对第i种污染物的净化量，单位为吨/年；A2为水体面积，单位为km2；/>为单位面积湿地生态系统对第i种污染物的净化量，单位为吨/（km2·年）。
基于上述技术方案，更进一步地，所述气候调节服务价值的核算方法，包括以下步骤：
采用替代成本法，根据达到同样效果的用电量和电价计算相应的价值量；
公式十八：；
式中，V为气候调节价值，单位为元/年；Q为生态系统蒸腾蒸发消耗能量，单位为（kw·h)/年，p为电价，单位为元/(kw·h)。
基于上述技术方案，更进一步地，所述生态系统蒸腾蒸发消耗能量的计算公式为：
公式十九：；
式中：Q为生态系统蒸腾蒸发消耗能量，单位为kw·h；GPP为不同生态系统类型单位面积蒸腾消耗热量；Av为第v类生态系统类型面积，单位为km2；R1为空调能效比；d为空调开放天数的，单位为天；i为研究区不同生态系统类型；EQ为水面蒸发量，单位为m3；γ为将1m3水转化为蒸汽的耗电量，单位为kw·h。
基于上述技术方案，更进一步地，所述病虫害防治服务价值的核算方法，包括以下步骤：
采用防护费用法，以人工林病虫害高出天然林造成的损失计算森林病虫害控制功能及价值，计算公式为：
公式二十：；
公式二十一：；
式中，为病虫害控制功能价值，单位为元/年；/>为病虫害控制功能量，单位为km2；/>为单位面积病虫害防治的费用，单位为元/（km2·年）；/>为天然林面积，单位为km2；/>为人工林病虫害发生率，单位为%；/>为天然林病虫害发生率。
相对于现有技术，本发明具有以下有益效果：
本发明通过采用各种价值化方法，如市场价值法、替代成本法、防护费用法等，将水源涵养服务、土壤保持服务、洪水调蓄服务、防风固沙服务、固碳释氧服务、大气净化服务、水质净化服务、气候调节服务和病虫害防治服务等各类生态系统调节服务进行整合，实现了整体价值化；并在核算过程中，设定了相应的修正方程，进一步提高了价值核算结果的准确度。
具体实施方式
实施例1
一种生态系统调节服务价值核算方法，包括以下步骤：
步骤S1：收集相应生态系统调节服务的基础数据，生态系统调节服务包括水源涵养服务、土壤保持服务、洪水调蓄服务、防风固沙服务、固碳释氧服务、大气净化服务、水质净化服务、气候调节服务和病虫害防治服务；
步骤S2：根据基础数据的数据类型采用不同的数据预处理方法得到相应生态系统调节服务的实物量数据，再分别对生态系统调节服务中的各个服务的价值量进行核算；例如统计数据采用数据对比与筛选，矢量数据采用GIS克里金插值，栅格数据采用GIS将栅格数据矢量化，进而获得实物量数据等方式；
具体的，水源涵养指生态系统对降水进行截留、渗透、蓄积，并通过蒸散发实现对水流、水循环的调控，是陆地生态系统重要的服务功能之一，其变化直接影响区域气候、水文、植被和土壤等状况，是区域生态系统状况的重要指示器。本项目参考已有研究，采用替代工程法计算水源涵养服务的经济价值。涵养水源服务的实物量采用水量平衡方程进行计算。目前，主要的功能量计算方法包括降水储存量法、多因子回归法、综合蓄水能力法、水量平衡法、林冠截留剩余量法以及森林水文模型法等。其中水量平衡法应用范围最广，从水量平衡的角度计算，结合遥感和气象资料，应用方便，结果可信度较高，该水源涵养服务价值的核算方法，计算公式为：
公式一：；
公式二：；
公式三：；
式中，为生态系统涵养水源的价值，单位为元/年；/>为涵养水源服务的功能量，单位为m3；/>为水库单位库容的建设成本，单位为元/m3；/>为降雨量，单位为mm；为地表径流量，单位为mm；/>为蒸散发量，单位为mm；/>为不同土地利用类型的面积，单位为m2；/>为不同土地利用类型的平均地表径流系数。
土壤保持服务也叫水土保持服务，是生态系统提供的一项重要调节服务，指生态系统防止土壤流失的侵蚀调控能力以及对泥沙的储积保持能力。水土保持是指防治水土流失、保护、改良与合理利用山区、丘陵区和风沙区水土资源，维护和提高土地生产力，以利于充分发挥水土资源的经济效益和社会效益，建立良好生态环境的综合性科学技术。土壤侵蚀的危害主要表现在以下方面：一是造成局地水土流失、土壤肥力下降，影响粮食安全；二是增加河流输沙量，影响下游地区的灌溉、污水处理和水库使用；三是造成水体中营养元素超标污染环境。因此本报告从保持土壤肥力、减少泥沙淤积等两个方面计算生态系统土壤保持功能价值。目前水土保持由三大类措施组成：水土保持农业技术措施、水土保持林草措施和水土保持工程措施。土壤侵蚀的定量评价与监测对于水土保持能力评估有着重要的意义，土壤侵蚀模型可以分为经验统计模型和物理过程模型。生态系统土壤保持实物量采用修正通用土壤流失方程（RUSLE），通过潜在土壤侵蚀量与实际土壤侵蚀量的差值来计算。
其中，地表覆被、土壤可蚀性、地形条件、气候条件等都会对水土保持效果有影响，因此水土保持服务能力指数为：
；
式中：K1为土壤可蚀性因子（无量纲）；S为坡度因子（无量纲）；NPP为植被净初级生产力，各因子需要归一化处理。
土壤流失方程模型是在原来USLE模型上的进一步改进，使其可以更好的适用于国内实际情况。CSLE可以用于计算坡面上的土壤流失量，并且形式较为简单，可以在不同地区进行推广，方程基本形式为：；式中：A为土壤侵蚀模数，单位为t·hm-2·a-1，即吨/(hm2·年)，主要指由降雨和径流引起的坡面细沟或细沟间侵蚀的年均土壤流失量；R为降雨侵蚀力因子，单位为MJ·mm·hm-2·h-1·a-1，即兆焦耳·mm/(hm2·小时·年)，它反映降雨引起土壤流失的潜在能力；K1为土壤可蚀性因子；L为坡长因子；S为坡度因子；B为植被覆盖与生物措施因子；T为耕作措施因子。
降雨侵蚀力因子R计算：
；
；
；
式中：取1，2……，24，指将一年划分为24个半月；/>是第g个半月的降雨侵蚀力；N1指所计算的时间序列，按五年顺序顺延更新；j取0，1，……，m；m是第i年第g个半月内侵蚀性降雨日的数量（侵蚀性降雨日指日雨量大于等于10mm）；/>是第i年第g个半月第j个降雨侵蚀性降雨量（mm）；如果某年某个半月内没有侵蚀性降雨量，即j=0，则/>=0；/>为参数，暖季（2-9月）取0.3937，冷季（10-12月，1-4月）取0.3101；/>为第g个半月平均降雨侵蚀力占多年平均降雨侵蚀力的比例。
土壤可蚀性因子K计算：
；
；式中：SAN为土壤砂砾含量（%）；SIL为土壤粉粒含量（%）；CLA为土壤粘粒含量（%）；C1为土壤有机质含量（%）。
坡子和坡度因子计算：
坡度的计算利用ArcGIS坡度计算工具，从ASTGTM2 30m DEM获取研究区域的坡度数据。
坡长计算公式如下：
；
式中：slope为坡度，单位为°；为坡长；
坡度因子计算公式如下：
；
坡长因子计算公式如下：
；
；
式中：L为坡长因子。
获得24个半月的植被覆盖度，结合24个半月降雨侵蚀力因子和土地利用类型计算B。具体计算方法为：
；
式中：为第i个半月降雨侵蚀力占全年降雨侵蚀力的比例，取值范围为0-1；为第i个半月园地、林地和草地的土壤流失比例，无量纲，取值范围为0-1，计算公式如下：
果园、其它果园、有林地和其他林地计算公式：
；
草地计算公式：
；
式中：是基于NDVI计算的植被覆盖程度，取值范围0-1；/>是乔木林的林下盖度，取值范围0-1，包括除乔木林冠层以外的所有植被（灌木、草本和枯落物）构成的林下盖度，按实地调查或经验取值，本方法选取0.55。
；
式中：NDVImax，NDVImin为NDVI最大值与最小值；k2为转换系数，取该像元内NDVI平均值。
非园地、林地和草地的B因子赋值如表1所示：
表1 B因子赋值表
其中，表1中的其他土地包括盐碱地、沙地、裸土地、裸岩石砾地等。
所述土壤保持服务价值的核算方法，包括以下步骤：
采用替代成本法进行计算，计算公式为：
公式四：；
公式五：；
公式六：；
公式七：；
式中，为生态系统土壤保持服务价值，单位为元/年；/>为保肥价值，单位为元/年；/>为减少泥沙淤积价值；/>为土壤保持量，单位为吨/年；/>为土壤中N、P、K营养元素含量，取表层土N、P、K含量的平均值分别为0.17%，0.08%，2.17%；/>为化肥的价格，单位为元，N、P、K分别转换成尿素、过磷酸钙、氯化钾三种肥料。
洪水调蓄主要体现在湿地生态系统减轻洪水威胁的经济价值。沼泽湿地是蓄水防洪的天然海绵，湖泊是抵御洪水灾害的天然屏障，水库是现代防洪工程体系的重要组成部分。本项目参考已有研究，运用替代工程法，通过建设水库的费用成本计算坑塘水面、沼泽、水库等湿地生态系统的洪水调蓄价值。其中，坑塘水面是人工或天然形成的需水量小于等于10万立方米的坑塘。本项目参照湖泊调蓄能力计算，公式如下：
；
；
式中，为坑塘水面洪水调蓄价值，单位为元；/>为坑塘洪水调蓄能力，单位为万m3；/>为水库单位库容建造成本。
所述洪水调蓄服务价值的核算方法，包括以下步骤：
按照水库总库容的35%计算水库洪水调蓄功能，按照替代工程法计算价值量，计算公式如下：
公式八：；
公式九：；
式中，为水库调蓄洪水的价值，单位为元/年；/>为水库调蓄洪水能力，单位为m3；/>为水库总库容，单位为m3；/>为人工的洪水调蓄价格，单位为元/(m3·年)。
沼泽洪水调蓄价值主要计算内陆滩涂在调蓄洪水的价值，与其面积及洪水期淹没深度有关。而确定沼泽洪水调蓄功能价值计算公式如下：
；
；
式中，为沼泽洪水调蓄价值，单位为元；/>为洪水调蓄功能量，单位为m3；为沼泽面积，单位为m2。
防风固沙是指生态系统通过其结构与过程，降低因风蚀导致的地表土壤裸露，增强地表粗糙程度，降低风速，减少风沙输沙量，减少土地沙化的功能。所述防风固沙服务价值的核算方法，包括以下步骤：
采用修正风蚀方程计算风蚀量，并利用替代成本法，将治理这些沙化土壤的成本和减少风蚀土壤肥力作为防风固沙功能的价值，计算公式为：
公式十：；
式中，为防风固沙价值，单位为元/年；/>为固沙量，单位为吨/年；/>为砂砾堆积密度，单位为吨/m3；/>为土壤沙化标准覆盖沙厚度，单位为m；/>为治沙工程的平均成本，单位为元/m2；/>为土壤中N、P、K营养元素含量（%），N、P、K分别转换成尿素、过磷酸钙、氯化钾三种肥料。例如，三种化肥的价格参考全国化肥网及陕西省市场价格；/>为化肥的价格，单位为元；/>为化肥中i种营养元素的比重（%），分别为尿素中氮含量46%、过磷酸钙中磷含量17%、氯化钾中钾含量57%。
其中，结合两个修正风蚀方程，对固沙量进行计算，计算公式如下：
；
式中，为固沙量，单位为吨/年；/>为栅格i的单位面积固沙量，单位为吨/(km2·年)；n为研究区域的栅格数量；A1为栅格面积，单位为km2。
采用修正风蚀方程一，计算得到实际风力侵蚀量SL：
；
；
；
采用修正风蚀方程二，计算得到潜在风力侵蚀量：
；
；
；
其中，；
；
；
；
；
；
；
；
对于每个栅格，单位面积固沙量的计算公式为：
；/>
式中，SR为单位面积固沙量，单位为吨/(km2·年)；为潜在风力侵蚀量，单位为吨/(km2·年)；SL为实际风力侵蚀量，单位为吨/(km2·年) ；/>为最大转移量，单位为kg/m；Z为最大风蚀出现距离，取50m；WF为气候因子，单位为kg/m；/>为地表糙度因子；EF为土壤可蚀因子；SCF为土壤结皮因子；C为植被覆盖因子；Wf为各风力因子，单位为m/s3；g为重力加速度，取9.8m/s2；/>为空气密度，取1.293kg/m3；/>为2m处起沙风速，取5m/s；为气象站2m处月均风速值，单位为m/s；/>为各月2m处日均风速大于5 m/s的天数；SW为土壤湿度因子，无量纲；/>为潜在蒸发量，单位为mm，采用彭曼系数计算；P为降雨量，单位为mm；I为灌溉量，单位为mm（本次取0）；/>为降雨天数；N为计算周期天数； SAN、SIL、CLA、SOM、CaCO3分别为土壤中砂砾、粉粒、粘粒、有机质及碳酸钙含量； SC为植被覆盖度（%）；/>为坡度。
固碳释氧是生态系统中的植物通过光合作用吸收空气中的二氧化碳，同时向空气中释放氧气的服务。根据光合作用方程式，植物生产1g干物质需要1.63g二氧化碳，并释放1.2g氧气。本项目采用植物净初级生产力计算生态系统干物质量，参考已有研究通过市场价值法，采用二氧化碳市场价格计算生态系统固碳价值，通过工业制氧成本法计算生态系统释氧价值。所述固碳释氧服务价值的核算方法，计算公式为：
公式十一：；
公式十二：；
公式十三：；
式中，为生态系统固碳释氧价值，单位为元/年；/>为生态系统固碳量，单位为吨/年；/>为生态系统释放氧气量，单位为吨/年；/>为二氧化碳价格，单位为元/吨；/>为氧气价格，单位为元/吨；NPP为植被净初级生产力，单位为吨/年。
生态系统具有一定的自净能力，能够通过一系列的物理、化学和生物因素的共同作用，使环境介质中人类生产生活所产生的污染物浓度降低。本项目选择滞尘量、二氧化硫和氮氧化物的吸收量来衡量生态系统大气净化服务价值。采用单位面积生态系统对颗粒物（PM2.5、PM10）、二氧化硫、二氧化氮的吸收量与各类生态系统面积的乘积作为生态系统大气环境净化的功能量，采用核算大气净化功能价值常用的替代成本法，以各类大气污染物治理成本计算价值量，所述大气净化服务价值的核算方法，计算公式为：
公式十四：；
公式十五：；
式中，为大气净化价值，单位为元/年；/>为大气污染物净化量，单位为吨；Qvj为第v类生态系统第 j 种大气污染物的单位面积净化量，v为生态系统类型，无量纲，j为大气污染物类别，无量纲；Av为第v类生态系统类型面积，单位为km2，cj为不同种类大气污染物处理成本，单位为元/吨。
净化水质指生态系统吸收、处理、过滤水体污染物，净化水质的功能。本项目采用植物净化模型核算湿地生态系统净化水质功能主要包括对COD、总氨、总磷的净化量，并参照已有研究，使用替代成本法根据达到同样净化效果所需的人工净化成本计算水质净化价值，所述水质净化服务价值的核算方法，计算公式为：
公式十六：；
公式十七：；
式中，为湿地生态系统水质净化价值，单位为元/年；/>为治理i种污染物的成本，单位为元/吨；/>为湿地生态系统对第i种污染物的净化量，单位为吨/年；A2为水体面积，单位为km2；/>为单位面积湿地生态系统对第i种污染物的净化量，单位为吨/km2·年；根据相关文献数据，确定湿地生态系统单位面积COD、总氮、总磷平均吸收量分别为110.43t/km2、22.81t/km2、5.97t/km2。
生态系统气候调节功能的价值指生态系统蒸腾、蒸发消耗能量从而节约用电量产生的价值，包括植被蒸腾和水面蒸发两方面。
公式十九：；
式中，Q为生态系统蒸腾蒸发消耗能量，单位为kw h；GPP为不同生态系统类型单位面积蒸腾消耗热量；Av为第v类生态系统类型面积，单位为km2；R1为空调能效比；d为空调开放天数的，单位为天；i为研究区不同生态系统类型；EQ为水面蒸发量，单位为m3；γ为将1m3水转化为蒸汽的耗电量，单位为kw h，取125kw h/m3。
所述气候调节服务价值的核算方法，包括以下步骤：
采用替代成本法，根据达到同样效果的用电量和电价计算相应的价值量；
公式十八：；
式中，V为气候调节价值，单位为元/年；Q为生态系统蒸腾蒸发消耗能量，单位为kwh/年，p为电价，单位为元/kw h。
生态系统通过增加物种多样性水平可以提高生态系统的稳定性，同时能够增加天敌而降低植食性昆虫的种群数量，从而起到控制病虫害的作用。本项目参考已有相关研究，采用防护费用法，以人工林病虫害高出天然林造成的损失计算森林病虫害控制功能及价值，所述病虫害防治服务价值的核算方法，计算公式为：
公式二十：；
公式二十一：；
式中，为病虫害控制功能价值，单位为元/年；/>为病虫害控制功能量，单位为km2；/>为单位面积病虫害防治的费用，单位为元/（km2·年）；/>为天然林面积，单位为km2；/>为人工林病虫害发生率，单位为%；/>为天然林病虫害发生率（%）。/>
步骤S3：对每个服务价值的核算结果进行求和。
具体的，水源涵养服务价值计算为公式一：，单位为元/年；
土壤保持服务价值计算为公式四：，单位为元/年；
洪水调蓄服务价值计算为公式八：，单位为元/年；
防风固沙服务价值计算为公式十：
，单位为元/年；
固碳释氧服务价值计算为公式十一：
，单位为元/年；
大气净化服务价值计算为公式十五：
，单位为元/年；
水质净化服务价值计算为公式十六：
，单位为元/年；
气候调节服务价值计算为公式十八：，单位为元/年；
病虫害防治服务价值计算为公式二十：，单位为元/年；
将公式一、公式四、公式八、公式十、公式十一、公式十五、公式十六、公式十八、公式二十相加计算的结果即为生态系统调节服务价值，单位为元/年。
最后应当说明的是，以上内容仅用以说明本发明的技术方案，而非对本发明保护范围的限制，本领域的普通技术人员对本发明的技术方案进行的简单修改或者等同替换，均不脱离本发明技术方案的实质和范围。
Claims (8)
1.一种生态系统调节服务价值核算方法，其特征在于，包括以下步骤：
步骤S1：收集相应生态系统调节服务的基础数据，生态系统调节服务包括水源涵养服务、土壤保持服务、洪水调蓄服务、防风固沙服务、固碳释氧服务、大气净化服务、水质净化服务、气候调节服务和病虫害防治服务；
步骤S2：根据基础数据的数据类型采用不同的数据预处理方法获得相应生态系统调节服务的实物量数据，再分别对生态系统调节服务中的各个服务的价值进行核算；其中，统计数据采用数据对比与筛选，矢量数据采用GIS克里金插值，栅格数据采用GIS将栅格数据矢量化，进而获得实物量数据；涵养水源服务的实物量采用水量平衡方程进行计算；
其中，水源涵养服务价值的核算方法，包括以下步骤：
采用替代工程法从水量平衡的角度计算水源涵养服务价值，计算公式为：
公式二：；
公式三：；
式中，为涵养水源服务的功能量，单位为m3；/>为降雨量，单位为mm；/>为地表径流量，单位为mm；/>为蒸散发量，单位为mm；/>为不同土地利用类型的面积，单位为m2；为不同土地利用类型的平均地表径流系数；
其中，土壤保持服务设有水土保持服务能力指数：
；
式中：K1为土壤可蚀性因子；S为坡度因子；NPP为植被净初级生产力；
土壤保持服务包括用于计算坡面上的土壤流失量的CSLE，方程为：；式中：A为土壤侵蚀模数，单位为t·hm-2·a-1，指由降雨和径流引起的坡面细沟或细沟间侵蚀的年均土壤流失量；R为降雨侵蚀力因子，单位为MJ·mm·hm-2·h-1·a-1，反映降雨引起土壤流失的潜在能力；L为坡长因子；B为植被覆盖与生物措施因子；T为耕作措施因子；
降雨侵蚀力因子R计算：
；
；
；
式中：g取1，2……，24，指将一年划分为24个半月；是第g个半月的降雨侵蚀力；N1指所计算的时间序列，按五年顺序顺延更新；j取0，1，……，m；m是第i年第g个半月内侵蚀性降雨日的数量；/>是第i年第g个半月第j个降雨侵蚀性降雨量；/>为第g个半月平均降雨侵蚀力占多年平均降雨侵蚀力的比例；
土壤可蚀性因子K1计算：
；
；式中：SAN为土壤砂砾含量；SIL为土壤粉粒含量；C1LA为土壤粘粒含量；C1为土壤有机质含量；
坡子和坡度因子计算：
坡度的计算利用ArcGIS坡度计算工具；
坡长计算公式：
；
式中：slope为坡度，单位为°；为坡长；
坡度因子计算公式：
；
坡长因子计算公式：
；
；
式中：L为坡长因子；
获得24个半月的植被覆盖度，结合24个半月降雨侵蚀力因子和土地利用类型计算植被覆盖与生物措施因子B，计算方法为：
；
式中：为第i个半月降雨侵蚀力占全年降雨侵蚀力的比例，取值范围为0-1；/>为第i个半月园地、林地和草地的土壤流失比例，无量纲，取值范围为0-1，计算公式为：
果园、其它果园、有林地和其他林地计算公式：
；
草地计算公式：
；
式中：是基于NDVI计算的植被覆盖程度，取值范围0-1；/>是乔木林的林下盖度，取值范围0-1；
；
式中：NDVImax，NDVImin为NDVI最大值与最小值；k2为转换系数；
其中，防风固沙服务价值的核算方法，包括以下步骤：
采用修正风蚀方程计算风蚀量，并利用替代成本法，将治理沙化土壤的成本和减少风蚀土壤肥力作为防风固沙功能的价值，计算公式为：
为固沙量，单位为吨/年；/>为砂砾堆积密度，单位为吨/m3；/>为土壤沙化标准覆盖沙厚度，单位为m；/>为治沙工程的平均成本；/>为土壤中N、P、K营养元素含量，N、P、K分别转换成尿素、过磷酸钙、氯化钾三种肥料；NPi为化肥中i种营养元素的比重；
其中，结合两个修正风蚀方程，对固沙量进行计算，计算公式：
；
式中，为栅格i的单位面积固沙量，单位为吨/(km2·年)；n为研究区域的栅格数量；A1为栅格面积，单位为km2；
采用修正风蚀方程一，计算得到实际风力侵蚀量SL：
；
；
；
采用修正风蚀方程二，计算得到潜在风力侵蚀量：
；
；
；
其中，；
；
；
；
；
；
；
；
对于每个栅格，单位面积固沙量的计算公式为：
；
式中，SR为单位面积固沙量，单位为吨/(km2·年)；为潜在风力侵蚀量，单位为吨/(km2·年)；SL为实际风力侵蚀量，单位为吨/(km2·年)；/>为最大转移量，单位为kg/m；Z为最大风蚀出现距离；WF为气候因子，单位为kg/m；/>为地表糙度因子；EF为土壤可蚀因子；SCF为土壤结皮因子；C为植被覆盖因子；Wf为各风力因子，单位为m/s3；g为重力加速度；为空气密度；/>为2m处起沙风速；/>为气象站2m处月均风速值，单位为m/s；/>为各月2m处日均风速大于5m/s的天数；SW为土壤湿度因子，无量纲；/>为潜在蒸发量，单位为mm；P为降雨量，单位为mm；I为灌溉量，单位为mm；/>为降雨天数；N为计算周期天数；SAN、SIL、CLA、CaCO3分别为土壤中砂砾、粉粒、粘粒及碳酸钙含量；SC为植被覆盖度；/>为坡度；
步骤S3：对每个服务价值的结果进行求和计算：
，式中，/>为土壤保持量，单位为吨/年；/>为土壤中N、P、K营养元素含量，取表层土N、P、K含量的平均值分别为0.17%，0.08%，2.17%；/>为化肥的价格，N、P、K分别转换成尿素、过磷酸钙、氯化钾三种肥料；/>为化肥中i种营养元素的比重/>为生态系统固碳量，单位为吨/年；/>为生态系统释放氧气量；/>为二氧化碳价格；/>为氧气价格；/>为水库总库容,单位为m3；/>为人工的洪水调蓄价格；/>为大气污染物净化量，单位为吨，cj为不同种类大气污染物处理成本；/>为湿地生态系统对第i种污染物的净化量；/>为治理i种污染物的成本；p为电价；/>为单位面积病虫害防治的费用；/>为水库单位库容的建设成本；Q为生态系统蒸腾蒸发消耗能量，单位为kw h；/>为病虫害控制功能量，单位为km2。
2.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，土壤保持服务价值的核算方法，包括以下步骤：
采用替代成本法进行计算，计算公式为：
公式五：；
公式六：；
公式七：；
式中，为保肥价值；/>为减少泥沙淤积价值。
3.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，洪水调蓄服务价值的核算方法，包括以下步骤：
按照水库总库容的35%计算水库洪水调蓄功能，按照替代工程法计算价值量，计算公式如下：
公式九：；
式中，为水库调蓄洪水能力，单位为m3；/>为水库总库容，单位为m3。
4.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，固碳释氧服务价值的核算方法，包括以下步骤：
通过工业制氧成本法计算生态系统释氧价值，计算公式为：
公式十二：；
公式十三：；
式中，为生态系统固碳量，单位为吨/年；/>为生态系统释放氧气量，单位为吨/年。
5.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，大气净化服务价值的核算方法，包括以下步骤：
采用替代成本法，以各类大气污染物治理成本计算价值量，计算公式为：
公式十四：；
式中，Qvj为第v类生态系统第j种大气污染物的单位面积净化量，v为生态系统类型，无量纲；j为大气污染物类别，无量纲；Av为第v类生态系统类型面积，单位为km2。
6.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，水质净化服务价值的核算方法，包括以下步骤：
采用替代成本法，根据达到同样净化效果所需的人工净化成本计算水质净化价值，计算公式为：
公式十七：；
式中，A2为水体面积，单位为km2；为单位面积湿地生态系统对第i种污染物的净化量，单位为吨/（km2·年）。
7.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，气候调节服务价值的核算方法，包括以下步骤：
采用替代成本法，根据达到同样效果的用电量和电价计算相应的价值量；
其中，生态系统蒸腾蒸发消耗能量的计算公式为：
公式十九：；
式中：GPP为不同生态系统类型单位面积蒸腾消耗热量；Av为第v类生态系统类型面积，单位为km2；R1为空调能效比；d为空调开放天数的，单位为天；EQ为水面蒸发量，单位为m3；γ为将1m3水转化为蒸汽的耗电量，单位为kw h。
8.根据权利要求1所述的一种生态系统调节服务价值核算方法，其特征在于，病虫害防治服务价值的核算方法，包括以下步骤：
采用防护费用法，以人工林病虫害高出天然林造成的损失计算森林病虫害控制功能及价值，计算公式为：
公式二十一：；
式中，为天然林面积，单位为km2；/>为人工林病虫害发生率，单位为%；/>为天然林病虫害发生率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211706450.3A CN115689395B (zh) | 2022-12-29 | 2022-12-29 | 一种生态系统调节服务价值核算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211706450.3A CN115689395B (zh) | 2022-12-29 | 2022-12-29 | 一种生态系统调节服务价值核算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115689395A CN115689395A (zh) | 2023-02-03 |
| CN115689395B true CN115689395B (zh) | 2023-09-05 |
Family
ID=85056052
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211706450.3A Active CN115689395B (zh) | 2022-12-29 | 2022-12-29 | 一种生态系统调节服务价值核算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115689395B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118195176A (zh) * | 2024-05-15 | 2024-06-14 | 四川省生态环境科学研究院 | 一种生态系统调节服务价值确定方法、装置、设备以及介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2243554C1 (ru) * | 2004-05-17 | 2004-12-27 | Некоммерческое Партнерство "Экспертно консультативный центр по сертификации материалов и эффективных технологий" | Способ комплексной оценки экологической обстановки и эффективности экологического менеджмента в регионе |
| KR20180000619A (ko) * | 2016-06-23 | 2018-01-03 | (주)해동기술개발공사 | Ｇｉｓ 기반 토사유실 평가 방법 |
| CN108805466A (zh) * | 2018-06-27 | 2018-11-13 | 南京林业大学 | 一种滨海湿地生态系统服务价值评价方法 |
| CN110458388A (zh) * | 2019-06-27 | 2019-11-15 | 中国水利水电科学研究院 | 基于自然理念的流域水生态服务价值的测算方法及装置 |
| AU2020103139A4 (en) * | 2020-10-30 | 2021-01-07 | Gengyuan LIU | Marine Ecological Monitoring And Its Assessment Methods |
-
2022
- 2022-12-29 CN CN202211706450.3A patent/CN115689395B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2243554C1 (ru) * | 2004-05-17 | 2004-12-27 | Некоммерческое Партнерство "Экспертно консультативный центр по сертификации материалов и эффективных технологий" | Способ комплексной оценки экологической обстановки и эффективности экологического менеджмента в регионе |
| KR20180000619A (ko) * | 2016-06-23 | 2018-01-03 | (주)해동기술개발공사 | Ｇｉｓ 기반 토사유실 평가 방법 |
| CN108805466A (zh) * | 2018-06-27 | 2018-11-13 | 南京林业大学 | 一种滨海湿地生态系统服务价值评价方法 |
| CN110458388A (zh) * | 2019-06-27 | 2019-11-15 | 中国水利水电科学研究院 | 基于自然理念的流域水生态服务价值的测算方法及装置 |
| AU2020103139A4 (en) * | 2020-10-30 | 2021-01-07 | Gengyuan LIU | Marine Ecological Monitoring And Its Assessment Methods |
Non-Patent Citations (1)
| Title |
|---|
| 茹河流域固原二期工程水土保持生态价值评价;王秀琴;;水土保持通报(第06期);240-244 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115689395A (zh) | 2023-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Baron et al. | Nitrogen fluxes in a high elevation Colorado Rocky Mountain basin | |
| Guan et al. | Assessing vegetation response to climatic variations and human activities: spatiotemporal NDVI variations in the Hexi Corridor and surrounding areas from 2000 to 2010 | |
| Qiu et al. | Water use efficiency and evapotranspiration of winter wheat and its response to irrigation regime in the north China plain | |
| Rosenberry et al. | Comparison of 13 equations for determining evapotranspiration from a prairie wetland, Cottonwood Lake area, North Dakota, USA | |
| Krogh | Field and village nutrient balances in millet cultivation in northern Burkina Faso: a village case study | |
| Guo et al. | Effects of water stress on water use efficiency and water balance components of Hippophae rhamnoides and Caragana intermedia in the soil–plant–atmosphere continuum | |
| CN115689395B (zh) | 一种生态系统调节服务价值核算方法 | |
| Olejnik et al. | A model for heat and water balance estimation and its application to land use and climate variation | |
| CN114936339B (zh) | 一种提升农业面源污染防治效果的模拟方法 | |
| Wang et al. | Spatial-temporal patterns of organic carbon sequestration capacity after long-term coastal wetland reclamation | |
| Wisniewski et al. | Carbon dioxide sequestration in terrestrial ecosystems | |
| CN113642842A (zh) | 一种基于部门报表的生态系统生产价值业务化核算方法 | |
| Ritschard et al. | SPATIAL AND TEMPORAL ANALYSIS OF AGRICULTURAL WATER REQUIREMENTS IN THE GULF COAST OF THE UNITED STATES 1 | |
| FOLEGA et al. | Riparian and stream forests carbon sequestration in the context of high anthropogenic disturbance in Togo | |
| McCoy et al. | Wind erosion and intensive prehistoric agriculture: A case study from the Kalaupapa field system, Moloka'i Island, Hawai'i | |
| Botkin | Global Warming: What It Is What Is Controversial about It and What We Might Do in Response to It | |
| Gunn et al. | Global energy balance and regional hydrology: A Burgundian case study | |
| de Medeiros et al. | Climate Change in Tacaimbó-PE, Brazil | |
| Allah | Research progress of grassland ecosystem: A case study of China | |
| Yang et al. | Land use impact on nitrogen discharge by stream: a case study in subtropical hilly region of China | |
| Zhang et al. | Wind mediates the responses of net ecosystem carbon balance to climatic change in a temperate semiarid steppe of Northern China | |
| Abbadie | Nitrogen inputs to and outputs from the soil-plant system | |
| Omasa et al. | Climatic change scenarios for monsoon Asia based on 2× CO 2-GCM experiments | |
| Chen et al. | RESEARCH ON CHANGES OF SOIL PROPERTIES AND QUALITY OF CULTIVATED LAND IN DIFFERENT UTILIZATION YEARS. | |
| Wang et al. | Characteristics and sustainability of land salinization and utilization in the Yellow River Delta |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |