趙寅成 孫雷 岳正波 吳楠

摘要 通過評價六安市生態(tài)系統(tǒng)為六安和周邊地區(qū)人們福祉和支撐經(jīng)濟(jì)社會發(fā)展所提供的產(chǎn)品和服務(wù)及其經(jīng)濟(jì)價值總和，基于指標(biāo)方法的規(guī)范和定期開展生態(tài)系統(tǒng)生產(chǎn)總值（GEP）核算的制度，核算了六安市2017年GEP。綜合分析核算結(jié)果可知，六安市2017年的GEP達(dá)3 794.48億元，其中生態(tài)系統(tǒng)調(diào)節(jié)功能總價值、生態(tài)系統(tǒng)產(chǎn)品提供總價值和生態(tài)系統(tǒng)文化功能總價值分別為3 140.8億、408.5億和245.2億元，分別占全市GEP的82.77%、10.77%和6.46%。六安市土地面積約15 461.16 km2，單位面積GEP為2 454.20萬元/km2。六安市GEP遠(yuǎn)遠(yuǎn)高于當(dāng)年的國內(nèi)生產(chǎn)總值（GDP）（1 218.7億元）?？芍彩蠫EP具有較好的發(fā)展?jié)摿?，可為促進(jìn)社會經(jīng)濟(jì)的可持續(xù)發(fā)展提供重要支持。

關(guān)鍵詞 生態(tài)系統(tǒng)生產(chǎn)總值;核算模型;安徽六安

中圖分類號 X171? 文獻(xiàn)標(biāo)識碼 A

文章編號 0517-6611（2021）04-0073-05

doi：10.3969/j.issn.0517-6611.2021.04.019

Study on Gross Ecosystem Product Accounting of Luan City，Anhui Province

ZHAO Yin-cheng1，2，SUN Lei3，YUE Zheng-bo1 et al （1.College of Resources and Environmental Engineering，Hefei University of Technology，Hefei，Anhui 230009;2.Wuhu Ecological Environment Bureau，Wuhu，Anhui 241000;3.Anhui Academy of Environmental Science，Hefei，Anhui 230071）

Abstract This study evaluated the products and services provided by the Luan ecosystem for the well-being of people in Luan and surrounding areas and supports the economic and social development and the sum of their economic value.Based on the standardization of the indicator method and the regular implementation of the system of gross ecosystem product （GEP） accounting，the 2017 GEP of Luan City was calculated.The comprehensive analysis and calculation results showed that the GEP of Luan City in 2017 reached 379.448 billion yuan，of which the total value of ecosystem regulation functions，the total value of ecosystem products and the total value of ecosystem cultural functions were 314.08 billion yuan，40.85 billion yuan and 24.52 billion yuan，accounting for 82.77%，10.77% and 6.46% of the citys GEP.The land area of Luan City was about 15 461.16 km2，and the GEP per unit area was 24.542 million yuan/km2.The GEP of Luan City was much higher than the GDP of that year （121.87 billion yuan）.It can be seen that GEP in Luan City has good development potential and can provide important support for promoting the sustainable development of social economy.

Key words Gross ecosystem product （GEP）;Accounting model;Luan City，Anhui Province

作者簡介 趙寅成（1991—），女，安徽桐城人，工程師，碩士，從事環(huán)境監(jiān)察工作。通信作者，高級工程師，博士，碩士生導(dǎo)師，從事環(huán)境規(guī)劃研究。

收稿日期 2020-07-07;修回日期 2020-07-16

生態(tài)系統(tǒng)生產(chǎn)總值（GEP）是用來評價與分析生態(tài)系統(tǒng)對人類經(jīng)濟(jì)社會發(fā)展支撐作用和對人類福祉的貢獻(xiàn)，是一套與GDP相對應(yīng)的、單獨的統(tǒng)計與核算體系[1-5]。生態(tài)系統(tǒng)的狀況和變化也可通過GEP核算認(rèn)識和了解[6]。GEP核算的重點是生態(tài)系統(tǒng)產(chǎn)品、調(diào)節(jié)服務(wù)和文化服務(wù)價值，非生態(tài)系統(tǒng)的經(jīng)濟(jì)價值是流量，而非存量，核算需限定時間和空間范圍[7-10]。該研究以六安市為對象，探討GEP核算的應(yīng)用方法，評價六安市生態(tài)系統(tǒng)為六安本地和周邊地區(qū)來源于人們福祉和支撐經(jīng)濟(jì)社會發(fā)展的產(chǎn)品和服務(wù)及其經(jīng)濟(jì)價值總和，基于指標(biāo)方法的規(guī)范和定期開展GEP核算的制度，定量反映六安市生態(tài)保護(hù)與建設(shè)進(jìn)展和成效，為全國其他地區(qū)“綠水青山”轉(zhuǎn)變?yōu)椤敖鹕姐y山”提供重要的示范引領(lǐng)作用。

1 資料與方法

1.1 研究區(qū)概況

六安市地處于安徽省西部、大別山北麓（115°20′～117°14′E、31°01′～32°40′N），東連合肥市，南接安慶市，西鄰河南信陽市，北與阜陽市隔河相望。東西寬176 km，南北長179 km，總面積15 459.2 km2（圖1）。六安市屬濕潤季風(fēng)氣候，地貌類型復(fù)雜多樣，河流隸屬淮河和長江流域。2017年，六安市能源消費總量538.21萬t標(biāo)準(zhǔn)煤，水資源總量為86.71億m3，生產(chǎn)總值（GDP）為1 218.7 億元，環(huán)境污染治理投資達(dá)19 332萬元。生態(tài)環(huán)境現(xiàn)狀方面，水環(huán)境質(zhì)量總體水質(zhì)為優(yōu)，化學(xué)需氧量下降趨勢明顯。城區(qū)環(huán)境空氣質(zhì)量達(dá)標(biāo)天數(shù)比例為80.5%，主要污染物為可吸入顆粒物和細(xì)顆粒物。城市和縣城生活污水集中處理率達(dá)92%。

1.2 六安市GEP核算方法研究

GEP的物質(zhì)量包含產(chǎn)品提供功能、調(diào)節(jié)功能和文化功能3個方面[11]。

1.2.1 生態(tài)系統(tǒng)產(chǎn)品價值核算。

具體核算方法參照下式[12]：

Vp=ni=1EPi×Pi（1）

式中，Vp為生態(tài)系統(tǒng)供給服務(wù)總值;EPi為第i類產(chǎn)品的產(chǎn)量;Pi為第i類產(chǎn)品的價格。

數(shù)據(jù)均來源于六安市統(tǒng)計年鑒。由于六安市不執(zhí)行峰谷電價，其電價按0.565 3元/（kW·h）計算。六安市自來水公司水費收費標(biāo)準(zhǔn)，2017年居民生活用自來水水價（居民階梯第二級）1.92元/t。

1.2.2 土壤保持服務(wù)功能價值核算。

具體評估方法參照下式[13-16]：

V1n=24%×Ac×C/ρ（2）

V1f=Ac×Ci×Pi/10 000（3）

式中，V1n為土壤保持的經(jīng)濟(jì)效益（元/a）;Ac為土壤保持量（t/a）;C為水庫清淤工程費用（元/m3）;ρ為土壤容重（t/m3）;V1f為減少面源污染經(jīng)濟(jì)效益（元/a）;Ci為土壤中氮、磷的純含量（kg/kg）;Pi為環(huán)境工程降解成本（元）。

1.2.3 洪水調(diào)蓄服務(wù)功能價值核算。具體評估方法參照下式[17-22]：

V5=c×Cr（4）

式中，V5為洪水調(diào)蓄價值;Cr為所有濕地（湖泊、水庫、沼澤）蓄水量;c為建設(shè)單位庫容的造價。

六安市內(nèi)各大中型湖庫2017年的蓄水動態(tài)數(shù)據(jù)來源于《2017年六安市水資源公報》，2017年六安市湖泊洪水調(diào)蓄能力17.19億m3，水庫洪水調(diào)蓄能力16.83億m3。由2017年遙感影像解譯得到的各縣區(qū)灘涂面積，按洪水期平均最大淹沒深度（1 m）推算沼澤地表滯水量，即沼澤洪水調(diào)蓄能力。

1.2.4 水源涵養(yǎng)服務(wù)功能價值核算。具體評估方法參照下式：

V4=TQ×c（5）

式中，V4為水源涵養(yǎng)總價值量（元/a），TQ為總水源涵養(yǎng)量（m3），c為建設(shè)成本?；诹彩兴Y源公報，計算得到六安市及各區(qū)縣2017年實際蒸散發(fā)量和水源涵養(yǎng)量。

1.2.5 大氣環(huán)境凈化服務(wù)功能價值核算。具體評估方法參照下式[23]：

V7a=2i=1ci×Qi（6）

式中，V7a為大氣環(huán)境凈化價值（元/a）;ci為治理大氣污染物的成本（元/t）;Qi為2017年大氣環(huán)境容量（t）。通過計算，得到 2017年六安市二氧化硫、氮氧化物和PM2.5的環(huán)境容量分別為5 354.70、3 449.92和5 086.97 t/a。

1.2.6 水質(zhì)凈化服務(wù)功能價值核算。具體評估方法參照下式[24]：

V7a=2i=1ci×Qi（7）

式中，V7a為生態(tài)系統(tǒng)水質(zhì)凈化的價值（元/a）;ci為治理水體污染物的成本;Qi為年水體污染物的凈化能力（t）。

根據(jù)《六安統(tǒng)計年鑒2018》，2017年六安市CODCr和氨氮的排放量分別達(dá)37 565.68和4 076.56 t/a。根據(jù)《2017年六安市水資源公報》，六安市CODCr和氨氮的入河排污量分別達(dá)9 779和1 827 t/a。由此得到2017年六安市水體污染物的凈化能力CODCr和氨氮分別為27 786.68和2 249.56 t/a。根據(jù)相關(guān)標(biāo)準(zhǔn)，CODCr征收標(biāo)準(zhǔn)為0.7元/1（CODCr當(dāng)量）×1 000=700元/t，氨氮征收標(biāo)準(zhǔn)為0.7元/0.8（氨氮當(dāng)量）×1 000=875元/t。

1.2.7 病蟲害控制服務(wù)功能價值核算。具體評估方法參照下式[25]：

Eb=NFa×（MFr-NFr）×Pb（8）

式中，Eb為功能價值（萬元），NFa為林業(yè)面積（km2）或農(nóng)業(yè)綜合防治面積（hm2），MFr為人工林或非綜合農(nóng)田防治病蟲害發(fā)生率（%），NFr為天然林或綜合防治農(nóng)田病蟲害發(fā)生率（%），Pb為單位面積病蟲害防治的費用（萬元/km2）。

通過統(tǒng)計年鑒可以得到六安市不同區(qū)域天然林面積和農(nóng)業(yè)病蟲害綜合防治面積，全市2017年天然林面積2 646.19 km2。通過對六安市相關(guān)專家的咨詢，六安市人工林和天然林的病蟲害發(fā)生率分別為10%和1%，防治成本約3.30萬元/km2。六安市農(nóng)作物耕作面積和農(nóng)業(yè)病蟲害綜合防治面積分別為6 677.16和1 335.43 km2。

1.2.8 固碳釋氧服務(wù)功能價值核算。具體評估方法參照下式[25]：

V8=（NPP/45%）×（1.62×CP+1.20×c）（9）

式中，V8為生態(tài)系統(tǒng)固碳釋氧服務(wù)功能的價值;NPP為凈初級生產(chǎn)力;CP為國際碳交易價格;c為醫(yī)用氧氣市場價格。基于GLO-PEM模型進(jìn)行核算。固碳價格采用962.7元/t CO2。40 L規(guī)格儲氣量氧氣的價格為30元。

1.2.9 氣候調(diào)節(jié)服務(wù)功能價值核算。具體評估方法參照下式 [26]：

V8=Qe×q×p×εα×3 600+β×Qe×p×η（10）

式中，V8為生態(tài)系統(tǒng)氣候調(diào)節(jié)的價值;Qe為年蒸散量（m3）;q為標(biāo)準(zhǔn)大氣壓下水的汽化熱（2.26×106 J/kg）;p為電價;α為空調(diào)能效比，取值3;ε為空調(diào)運行系數(shù)，一般設(shè)為0.130;β為單位體積水蒸發(fā)耗電量（125 kW·h/m3）;η為加濕器運行系數(shù)，一般設(shè)為0.123。電價取0.565 3元/（kW·h）。

1.2.10 文化服務(wù)功能價值核算。具體評估方法參照下式 [27]：

UV=CC+CS（11）

CC=TC+TV（12）

TV=H×W（13）

CS=∫pm0f（x）dx（14）

式中，UV為文化服務(wù)功能價值，CC為消費者的總支出，CS為消費者的剩余值;TC為旅行費用，TV為旅行時間價值;H為游客在景點游覽的時間（h），W為工資率（元/h）;x為追加費用，f（x）為需求曲線;Pm為最大追加費用。

2 六安市GEP總體特征

2.1 生態(tài)系統(tǒng)產(chǎn)品提供價值特征 由表1和圖2可知，六安市生態(tài)系統(tǒng)產(chǎn)品提供總價值為408.50億元，占GEP總價

值的10.77%。其中，農(nóng)業(yè)產(chǎn)品價值160.56億元，占GEP總

價值的4.23%;林業(yè)產(chǎn)品價值38.07億元，占比1.00%;畜牧業(yè)產(chǎn)品價值124.01億元，占比3.27%;漁業(yè)產(chǎn)品價值35.11億元，占比0.93%;水資源價值45.18億元，占比1.19%;生態(tài)能源價值5.57億元，占比0.15%?？梢?，農(nóng)業(yè)產(chǎn)品價值占GEP總價值最高，最低占比是生態(tài)能源價值，前者約為后者30倍。

2.2 生態(tài)系統(tǒng)調(diào)節(jié)功能價值特征

從表2可以看出，六安市生態(tài)系統(tǒng)調(diào)節(jié)功能總經(jīng)濟(jì)價值為3 140.79億元，占GEP總價值的82.77%。六安市氣候調(diào)節(jié)功能價值為1 268.64億元，貢獻(xiàn)最大（33.43%）;其次為水源涵養(yǎng)功能價值，為668.41億元，占比17.62%;固碳釋氧功能價值為625.12億元，占比16.47%;洪水調(diào)蓄功能價值254.96億元，占比6.72%;空氣凈化功能價值249.20億元，占比6.57%;其余的土壤保持、水質(zhì)凈化以及病蟲害防治功能價值總計為74.46億元，僅占GEP總價值的1.96%?；谒科胶饽Ｐ涂芍彩猩鷳B(tài)系統(tǒng)蓄水服務(wù)價值為668.41億元。

依據(jù)通用水土流失方程（USLE）計算得到，因生態(tài)系統(tǒng)土壤保持功能而致使全市減輕水體中泥沙淤積的體積為1.14億m3。土壤中氮、磷含量來源于1∶100萬土壤數(shù)據(jù)庫，得到全市減少氮面源污染63.28萬t，減少磷面源污染31.07萬t。單位總氮和總磷的排污費征收標(biāo)準(zhǔn)分別為1 750和2 800元/t，土方的挖掘成本為47.3元/m3，可知六安生態(tài)系統(tǒng)減少泥沙淤積價值為54.07億元，減少氮源污染價值為11.07億元，減少磷污染價值為8.70億元，土壤保持功能總價值為73.84億元。

通過調(diào)查發(fā)現(xiàn)，六安市湖泊和水庫年調(diào)蓄洪水能力分別占總調(diào)蓄能力的47.94%和46.93%，依據(jù)洪水期沼澤的地表滯水的常量，計算得到沼澤年調(diào)蓄洪水能力1.84億m3，占比5.13%。水庫單位庫容的建造成本為7.11元/m3，因此六安濕地生態(tài)系統(tǒng)洪水調(diào)蓄功能的經(jīng)濟(jì)價值總量為254.96億元。

六安市大氣環(huán)境SO2容量5 354.7 t，氮氧化物3 449.9 t，細(xì)顆粒物5 087.0 t。硫化物和氮氧化物的排污費征收標(biāo)準(zhǔn)分別為1.26和0.63元/kg， PM2.5排放造成的健康危害經(jīng)濟(jì)損失為4 897.0元/kg。計算可知，2017年六安市生態(tài)系統(tǒng)凈化SO2、NOX以及細(xì)顆粒物的價值分別為0.07億、0.02億和249.11億元。

根據(jù)六安市各流域水體污染物排放量與水資源公報中的水污染物入河量差值，基于CODCr排污費征收標(biāo)準(zhǔn)（700元/t）和氨氮排污費征收標(biāo)準(zhǔn)（875元/t），計算得到2017年六安市生態(tài)系統(tǒng)凈化水污染物CODCr和氨氮的價值分別為0.19億和0.02億元。

根據(jù)光合作用和呼吸作用的方程式計算得到2017年全市生態(tài)系統(tǒng)釋放O2 1 216.03萬t，固定CO2 1 641.64萬t，依據(jù)2017年安徽醫(yī)用氧氣市場價格（3 841元/t）和瑞典碳稅價格（962.7元/t），可知六安市的生態(tài)系統(tǒng)釋氧價值為467.08億元，固碳價值為158.04億元，總價值為625.12億元。

根據(jù)擬合Budyko曲線的經(jīng)驗公式和不同生態(tài)系統(tǒng)類型上的經(jīng)驗參數(shù)，可知六安市森林生態(tài)系統(tǒng)蒸散發(fā)降溫增濕消耗能量817.17億kW·h，灌叢728.68億kW·h，草地564.64億kW·h，水面133.70億kW·h。基于安徽省居民電費為0.565 3元/（kW·h），可知全市生態(tài)系統(tǒng)氣候調(diào)節(jié)的總價值達(dá)到1 268.64億元。

六安市天然林面積和農(nóng)作物播種面積分別為2 646萬和6 677.16 km2，人工林病蟲害發(fā)生率約為10%，天然林病蟲害發(fā)生率約為1%，非綜合防治農(nóng)田病蟲害發(fā)生率為25%～35%，綜合防治農(nóng)田病蟲害發(fā)生率為10%～15%。2017年森林病蟲害平均防治費用約32 983.5元/km2，農(nóng)田病蟲害平均防治費用約149 925.0元/km2，可得六安市生態(tài)系統(tǒng)病蟲害控制價值為0.41億元。

2.3 生態(tài)系統(tǒng)文化功能價值特征

六安市年均接待海外游客和國內(nèi)游客分別為12.7萬和4 489.2萬人次，旅游總收入達(dá)350.27億;文化服務(wù)價值量達(dá)245.19億元，占GEP總價值的6.46%。

3 結(jié)論與討論

通過對六安市生態(tài)系統(tǒng)產(chǎn)品提供功能、調(diào)節(jié)功能和文化功能的指標(biāo)確定與模型構(gòu)建，研究了六安市GEP總體特征。具體結(jié)論如下：

（1）2017年六安市GEP為3 794.48億元，其中生態(tài)系統(tǒng)產(chǎn)品提供功能、調(diào)節(jié)功能和文化功能的總價值分別為408.50億、3 140.79億和245.19億元，分別占全市GEP的10.77%、82.77%和6.46%。六安市土地面積15 461.16 km2，單位面積GEP為2 454.20萬元/km2，2017年六安市國內(nèi)生產(chǎn)總值（GDP）總量為1 218.7億元，僅占GEP的32.12%。

（2）六安市生態(tài)系統(tǒng)產(chǎn)品提供總價值中，農(nóng)業(yè)產(chǎn)品價值160.56億元;林業(yè)產(chǎn)品價值38.07億元;畜牧業(yè)產(chǎn)品價值124.01億元;漁業(yè)產(chǎn)品價值35.11億元;水資源價值45.18億元;生態(tài)能源價值5.57億元。

（3）六安市生態(tài)系統(tǒng)調(diào)節(jié)功能總經(jīng)濟(jì)價值中氣候調(diào)節(jié)功能價值最高，為1 268.64億元，占GEP總價值的33.43%;其次為水源涵養(yǎng)功能價值為668.41億元，占比17.62%;固碳釋氧功能價值為625.12億元，占比16.47%;洪水調(diào)蓄功能價值254.96億元，占比6.72%;空氣凈化功能價值249.20億元，占比6.57%;其余的土壤保持、水質(zhì)凈化以及病蟲害防治功能價值總計為74.46億元，僅占GEP總價值的1.96%。

參考文獻(xiàn)

[1]CAMACHO-VALDEZ V，RUIZ-LUNA A，GHERMANDI A，et al.Effects of land use changes on the ecosystem service values of coastal wetlands[J].Environmental management，2014，54（4）：852-864.

[2] CHHINH N，LAWN P.A green measure of Cambodias GDP to assist in achieving sustainability[J].Environmental research journal，2009，3（4）：259.

[3] DURAIAPPAH A K.Ecosystem services and human well-being：Do global findings make any sense？[J].Bioscience，2011，61（1）：7-8.

[4] 李濤，甘德欣，楊知建，等.土地利用變化影響下洞庭湖地區(qū)生態(tài)系統(tǒng)服務(wù)價值的時空演變[J].應(yīng)用生態(tài)學(xué)報，2016，27（12）：3787-3796.

[5] 馬國霞，於方，王金南，等.中國2015年陸地生態(tài)系統(tǒng)生產(chǎn)總值核算研究[J].中國環(huán)境科學(xué)，2017，37（4）：1474-1482.

[6] TANG X G，ZHOU Y L，LI H P，et al.Remotely monitoring ecosystem respiration from various grasslands along a large-scale east-west transect across northern China[J].Carbon balance and management，2020，15（1）：1-14.

[7] ZHAO A Z，ZHANG A B，CAO S，et al.Spatiotemporal patterns of water use efficiency in China and responses to multi-scale drought[J].Theoretical and applied climatology，2020，140（1/2）：559-570.

[8] 馬國霞，趙學(xué)濤，吳瓊，等.生態(tài)系統(tǒng)生產(chǎn)總值核算概念界定和體系構(gòu)建[J].資源科學(xué)，2015，37（9）：1709-1715.

[9] COSTANZA R，DARGE R，DE GROOT R，et al.The value of the worlds ecosystem services and natural capital[J].Nature，1997，387：253-260.

[10] 于淼，金海珍，李強(qiáng)，等.呈貢區(qū)生態(tài)系統(tǒng)生產(chǎn)總值（GEP）核算研究[J].西部林業(yè)科學(xué)，2020，49（3）：41-48，55.

[11] 牟雪潔，饒勝，蔣洪強(qiáng)，等.北京市延慶區(qū)生態(tài)環(huán)境（綠色GDP）核算體系研究[J].環(huán)境生態(tài)學(xué)，2019，1（3）：1-9.

[12] 馬國霞，於方，王金南，等.中國2015年陸地生態(tài)系統(tǒng)生產(chǎn)總值核算研究[J].中國環(huán)境科學(xué)，2017，37（4）：1474-1482.

[13] 陳媛媛，朱記偉，周蓓.基于系統(tǒng)動力學(xué)的西安市復(fù)合生態(tài)系統(tǒng)情景分析[J].水資源與水工程學(xué)報，2018，29（6）：31-40.

[14] 歐陽志云，鄭華，謝高地，等.生態(tài)資產(chǎn)、生態(tài)補(bǔ)償及生態(tài)文明科技貢獻(xiàn)核算理論與技術(shù)[J].生態(tài)學(xué)報，2016，36（22）：7136-7139.

[15] WANG D G，SHI A Q，WANG X B，et al.Assessment of island land ecosystem services value and their spatial distribution[J].Nature environment & pollution technology，2014，13（4）：743-748.

[16] WEIR K.Rising global value of ecosystem services dwarfs GDP[J].Pacific ecologist，2015（23）：57-63.

[17] 王金南，於方，曹東.中國綠色國民經(jīng)濟(jì)核算研究報告2004[J].中國人口·資源與環(huán)境，2006，16（6）：11-17.

[18] BUTTERFIELD Z，BUERMANN W，KEPPEL-ALEKS G.Satellite observations reveal seasonal redistribution of northern ecosystem productivity in response to interannual climate variability[J/OL].Remote sensing of environment，2020，242[2020-03-05].https/doi.org/10.1016/j.rse.2020.111755.

[19] CAVALLETTI B，DI FABIO C，LAGOMARSINO E，et al.Ecosystem accounting for marine protected areas：A proposed framework[J/OL].Ecological economics，2020，173[2020-03-05].https：//doi.org/10.1016/j.ecolecon.2020.106623.

[20] 於方，蔣洪強(qiáng)，曹東，等.中國綠色國民經(jīng)濟(jì)核算技術(shù)體系與方法概論[J].環(huán)境保護(hù)，2006，34（18）：30-39.

[21] BANERJEE O，BAGSTAD K J，CICOWIEZ M，et al.Economic，land use，and ecosystem services impacts of Rwandas Green Growth Strategy：An application of the IEEM +ESM platform[J/OL].Science of the total environment，2020，729[2020-03-05].https：//doi.org/10.1016/j.scitotenv.2020.138779.

[22] CIVEIRA G，LADO LIARES M，VIDAL VAZQUEZ E，et al.Ecosystem services and economic assessment of land uses in urban and periurban areas[J].Environmental management，2020，65（3）：355-368.

[23] 曾杰，李江風(fēng)，姚小薇.武漢城市圈生態(tài)系統(tǒng)服務(wù)價值時空變化特征[J].應(yīng)用生態(tài)學(xué)報，2014，25（3）：883-891.

[24] 孫秀美.基于GIS的山東沂沭泗流域水土保持生態(tài)功能評價及區(qū)劃研究[D].濟(jì)南：山東師范大學(xué)，2005.

[25] SMITH S V，HOLLIBAUGH J T.Annual cycle and interannual variability of ecosystem metabolism in a temperate climate embayment[J].Ecological monographs，1997，67（4）：509-533.

[26] SUN G，CALDWELL P，NOORMETS A，et al.Upscaling key ecosystem functions across the conterminous United States by a water-centric ecosystem model[J].Journal of geophysical research：Biogeosciences，2011，116：1-16.

[27] ZHAO M S，RUNNING S W，NEMANI R R.Sensitivity of Moderate Resolution Imaging Spectroradiometer （MODIS） terrestrial primary production to the accuracy of meteorological reanalyses[J].Journal of geophysical research：Biogeosciences，2006，111：1-13.