摘要? 以寧晉縣為研究區(qū)，基于土地利用/覆被社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)資料等數(shù)據(jù)，評(píng)估不同時(shí)期的生態(tài)系統(tǒng)服務(wù)價(jià)值，分析其損益變化權(quán)衡關(guān)系的時(shí)空規(guī)律。結(jié)果表明：1995—2019年生態(tài)系統(tǒng)多功能性、生產(chǎn)功能和生態(tài)功能整體呈先增加后減少；氣候調(diào)節(jié)、養(yǎng)分固持、氣體調(diào)節(jié)功能呈增加趨勢(shì)；水源涵養(yǎng)功能明顯下降；污染物容納功能持續(xù)降低。耕地面積減少導(dǎo)致流失的生態(tài)系統(tǒng)服務(wù)最大。生態(tài)系統(tǒng)多功能性、生態(tài)功能、水源涵養(yǎng)功能與其他功能由協(xié)同關(guān)系變?yōu)闄?quán)衡關(guān)系。

關(guān)鍵詞? 生態(tài)系統(tǒng)服務(wù)價(jià)值；權(quán)衡分析；時(shí)空演變；土地利用

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

文章編號(hào)? 0517-6611（2024）04-0036-07

doi：10.3969/j.issn.0517-6611.2024.04.009

開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Temporal and Spatial Variation Trade.offs Analysis of Ecosystem Services in Ningjin County of Hebei Province

ZHOU Li.na

（Baoding University，Baoding，Hebei 071000）

Abstract? Taking Ningjin County as the research area，base on the LUCC and social and economic statistics，the value of ecosystem services in different periods were estimated，the temporal and spatial change rule of the trade.off relationship between profit and loss were analyzed.The results showed that during 1995-2019，multi.functionality of ecosystem，production function and ecological function showed increased first and then decreased;the functions of climate regulation，nutrient fixation and gas regulation were showing an increasing trend; the water conservation function had significantly decreased; the pollutant containment function continued to decrease.The loss of ecosystem services caused by the reduction of cultivated land area was the most.Multi.functionality，ecological function，water conservation function and other functions of the ecosystem had changed from synergy to trade.off.

Key words? Ecosystem service;Trade.offs analysis;Temporal and spatial change;Land use

基金項(xiàng)目? 河北省高等學(xué)?？茖W(xué)技術(shù)研究項(xiàng)目（ZC2021002）。

作者簡(jiǎn)介? 周麗娜（1983—），女，河北廊坊人，講師，碩士，從事農(nóng)業(yè)資源與環(huán)境研究。

收稿日期? 2023-03-21

生態(tài)系統(tǒng)服務(wù)是人類從生態(tài)系統(tǒng)中獲得的所有收益，是人類生存和發(fā)展的基礎(chǔ)［1］，然而隨著社會(huì)發(fā)展和人類急劇增長(zhǎng)對(duì)食物供給的要求不斷增加，對(duì)生態(tài)系統(tǒng)造成巨大壓力，生態(tài)系統(tǒng)服務(wù)能力正在受到嚴(yán)重削弱［2］，人類當(dāng)前面臨的多種生態(tài)問(wèn)題的本質(zhì)是由于生態(tài)系統(tǒng)功能與服務(wù)受到破壞與退化的后果［3-5］。國(guó)內(nèi)外土地生態(tài)系統(tǒng)服務(wù)研究多采用遙感影像解譯獲取 LUCC 數(shù)據(jù)，應(yīng)用 GIS 進(jìn)行時(shí)空變化分析［6-19］，使用 InVEST 模型模擬生態(tài)系統(tǒng)服務(wù)［20-22］，為生態(tài)補(bǔ)償［23-24］、生態(tài)安全格局構(gòu)建［25-27］、土地利用管理與規(guī)劃［28］提供基礎(chǔ)信息。

寧晉縣地處太行山東麓沖積平原，地勢(shì)低平，境內(nèi)有洨河、滏陽(yáng)河等河流，是國(guó)家重要的商品糧生產(chǎn)基地，也是冀中平原重要的生態(tài)屏障區(qū)。如何通過(guò)合理調(diào)控土地利用保障糧食安全、緩解生態(tài)壓力，實(shí)現(xiàn)不同類型生態(tài)系統(tǒng)功能協(xié)同發(fā)展是該區(qū)域面臨的重要問(wèn)題之一。筆者通過(guò)指標(biāo)量化，研究寧晉縣生態(tài)系統(tǒng)服務(wù)功能的價(jià)值變化，采用生態(tài)系統(tǒng)多功能的響應(yīng)比進(jìn)行權(quán)衡分析，為把握區(qū)域生態(tài)系統(tǒng)服務(wù)變化、構(gòu)建生態(tài)安全格局提供研究依據(jù)。

1? 資料與方法

1.1? 數(shù)據(jù)來(lái)源

基礎(chǔ)數(shù)據(jù)包括：①寧晉縣1995、2003、2011、2019年4期土地利用基礎(chǔ)數(shù)據(jù)，數(shù)據(jù)源為L(zhǎng)andsat7 ETM、Landsat8-9 OLI數(shù)據(jù)。通過(guò)非監(jiān)督分類和人工目視判讀相結(jié)合的方法進(jìn)行遙感解譯分為6類用地，空間分辨率為30 m，混合矩陣檢驗(yàn)Kappa系數(shù)在88%以上。②《河北省統(tǒng)計(jì)年鑒》《中國(guó)農(nóng)產(chǎn)品成本匯編》等。

1.2? 生態(tài)系統(tǒng)功能評(píng)價(jià)與研究方法

1.2.1? 生態(tài)系統(tǒng)功能評(píng)價(jià)方法。

根據(jù)《中國(guó)生物多樣性國(guó)情研究報(bào)告》［5］對(duì)生態(tài)系統(tǒng)功能進(jìn)行經(jīng)濟(jì)價(jià)值評(píng)估，構(gòu)建以下評(píng)估體系（表1）。森林和草地生態(tài)系統(tǒng)生態(tài)價(jià)值評(píng)估及相關(guān)系數(shù)計(jì)算參考《森林生態(tài)系統(tǒng)功能評(píng)估規(guī)范》［29］和相關(guān)研究［30-34］。

（1）生產(chǎn)功能。森林和淡水生產(chǎn)功能價(jià)值以林業(yè)和漁業(yè)總產(chǎn)值表示。農(nóng)田和草地生態(tài)系統(tǒng)生產(chǎn)功能以農(nóng)產(chǎn)品和牧草經(jīng)濟(jì)價(jià)值表示。計(jì)算公式如下：

EC=ni=1（gi×Yi）

Eg=Ac×Y草×H

式中：EC為農(nóng)田生產(chǎn)價(jià)值（元）；i為主要作物種類（小麥、玉米、大豆、花生、棉花和蔬菜）；Yi為某作物當(dāng)年產(chǎn)量（t/hm2）；gi為某農(nóng)作物當(dāng)年市場(chǎng)價(jià)格（元/t）。Eg為草地生產(chǎn)價(jià)值（元）；Ac為草地面積（hm2）；Y草為牧草產(chǎn)量（t/hm2），參考已有研究取2.53 t/hm2［30］；H為牧草價(jià)格（元/t），按照近年市場(chǎng)價(jià)800元/t計(jì)算。

（2）氣候調(diào)節(jié)。農(nóng)田生態(tài)系統(tǒng)氣候調(diào)節(jié)包括固碳和釋氧功能價(jià)值。按照作物每合成單位有機(jī)干物質(zhì)所吸收的碳量和釋氧量的方法［35］計(jì)算。計(jì)算公式如下：

MC=ni=1ci×Yi×（1-Wi）×T碳HIi

MO=ni=1ri×Yi×（1-Wi）×C制氧HIi

式中：MC 為農(nóng)田生態(tài)系統(tǒng)碳蓄積量（ t/a）；MO 為農(nóng)田生態(tài)系統(tǒng)氧氣的釋放量（ t/a ）；ci 為作物碳吸收率（%）；Yi 為作物經(jīng)濟(jì)產(chǎn)量（ t/a ）；Wi 為作物經(jīng)濟(jì)產(chǎn)品部分的含水量（%）；HIi為作物經(jīng)濟(jì)系數(shù)；ri為釋氧量（t/a）;T碳根據(jù)瑞典碳稅法，按 2018 年 8 月匯率計(jì)算為1 200元/t;C制氧為工業(yè)制氧價(jià)格1 000元/t。

（3）養(yǎng)分固持。農(nóng)田生態(tài)系統(tǒng)養(yǎng)分固持價(jià)值可通過(guò)秸稈還田價(jià)值表示，參考秸稈資源相關(guān)研究［29，36］，秸稈養(yǎng)分折算計(jì)算方法如下：

FCN=ni=1［yi×Si×N耕×C1/R1］

FCP=ni=1［yi×Si×P耕×C1/R2］

FCK=ni=1［yi×Si×K耕×C2/R3］

FC=FCN+FCP+FCK

式中：FC為農(nóng)田養(yǎng)分固持價(jià)值量（元）；Si為作物谷草比［37］；N耕、P耕、K耕分別為作物秸稈含N、P、K的百分比，根據(jù)相關(guān)研究計(jì)算作物養(yǎng)分含量［38］；R1、R2、R3分別為化肥磷酸二銨含氮量（%）、磷含量（%）、氯化鉀化肥鉀含量（%），按照14.00%、 15.01%、 50.00%計(jì)算;C1、C2采用 2007 年中國(guó)農(nóng)業(yè)信息網(wǎng)公示平均價(jià)格磷酸二銨為2 400 元/t、氯化鉀為2 200 元/t。

（4）水源涵養(yǎng)。水源涵養(yǎng)價(jià)值包括水源調(diào)節(jié)價(jià)值和水分凈化價(jià)值。采用影子工程法和替代工程法進(jìn)行估算［29］。耕地水源涵養(yǎng)價(jià)值計(jì)算如下：

W耕調(diào)節(jié)=10×Ag×（P-E耕-θ耕）×L1

W耕凈化=10×Ag×（P- E耕-θ耕）×L2

式中：W耕調(diào)節(jié)為耕地水源調(diào)節(jié)價(jià)值（元）；W耕凈化為水源凈化價(jià)值（元）；P為降雨量（mm/a）；E耕為耕地水分的蒸散量（mm/a）；Ag為耕地面積（hm2）；θ耕為地表徑流量（mm/a）；L1為修建水庫(kù)的造價(jià)（元）；L2為水的凈化費(fèi)用（元）。依據(jù)文獻(xiàn)［39］可知，耕地土壤年蒸發(fā)量為196.58 mm，地表徑流系數(shù)以河北平原徑流系數(shù)的平均值計(jì)算為32%［40］。

水域水源涵養(yǎng)功能經(jīng)濟(jì)價(jià)值計(jì)算公式如下：

W水調(diào)節(jié)=10×Aw×P×L1

W水凈化=10×Aw×P×L2

式中：W水調(diào)節(jié)為水源涵養(yǎng)功能價(jià)值（元）；Aw為水域面積（hm2）；W水凈化為水的凈化價(jià)值（元）。其他指標(biāo)同前。

（5）氣體調(diào)節(jié)。

林地的氣體調(diào)節(jié)功能主要計(jì)算SO2的吸收和阻滯粉塵的功能價(jià)值。參考《中國(guó)生物多樣性國(guó)情研究報(bào)告》［5］取值。廢氣排污費(fèi)按照河北省2020年的排污費(fèi)收費(fèi)標(biāo)準(zhǔn)計(jì)算。

（6）污染物容納。

以農(nóng)田化肥流失損失費(fèi)用和農(nóng)藥、地膜污染造成產(chǎn)量經(jīng)濟(jì)損失來(lái)表征污染物容納能力價(jià)值量。化肥流失污染計(jì)算公式如下：

LF=-ni=1（1-q）×（LN×1814×0.70.8+LP×14262×0.70.05）

式中：q為化肥利用率（%）;LN為折純氮肥量（t）;LP為折純P2O5肥料量（t）;18/14為N元素與氨氮的摩爾系數(shù);0.8為氨氮的污染當(dāng)量值;0.7為每當(dāng)量污染物收費(fèi)額;142/62為P2O5和元素P的摩爾系數(shù);0.05為元素P的污染當(dāng)量值。

地膜污染造成損失以農(nóng)作物減產(chǎn)的經(jīng)濟(jì)損失表達(dá)，計(jì)算公式如下：

LG=-ni=1（Yi×di×Pi）

式中：Yi 為作物經(jīng)濟(jì)產(chǎn)量（t/a）;di為作物產(chǎn)量損失率（%）［41-42］;Pi為某作物的市場(chǎng)價(jià)格（元）。

農(nóng)藥污染損失的計(jì)算采用影子工程法，按農(nóng)藥主要成分氯氰菊酯參照工業(yè)水污染排污標(biāo)準(zhǔn)計(jì)算生產(chǎn)廢水，以其生化需氧量BOD5污染收費(fèi)來(lái)評(píng)估農(nóng)藥污染價(jià)值。計(jì)算公式如下：

LA=-Yi×T×LCOD×k×b

式中：Yi 為作物經(jīng)濟(jì)產(chǎn)量（t/a）;T為轉(zhuǎn)位廢水量（t/t）;LCOD為化學(xué)需氧量（mg/L）;k為COD與BOD5的轉(zhuǎn)換系數(shù)［43］;b為每當(dāng)量排污收費(fèi)額0.7元。

1.2.2? 權(quán)衡系數(shù)計(jì)算。

采用響應(yīng)比［44］計(jì)算權(quán)衡系數(shù)C，反映生產(chǎn)功能增減比例與生態(tài)功能增減比例，權(quán)衡評(píng)價(jià)生態(tài)系統(tǒng)功能性，計(jì)算公式如下：

Cc-o=ln（Ec/Eo）+ln（Yc/Yo ）

RE=ln（Ec/Eo）

RY=ln（Yc/Yo）

式中：Cc-o表示不同年限生態(tài)系統(tǒng)功能的權(quán)衡系數(shù)；Ec/Eo表示不同年限生產(chǎn)功能的比例；Yc/Yo表示不同年限生態(tài)功能的比例；RE為生產(chǎn)功能權(quán)衡系數(shù)；RY為生態(tài)功能權(quán)衡系數(shù)。

2? 結(jié)果與分析

2.1? 土地利用變化分析

遙感影像解譯土地利用類型如圖1所示，土地利用面積占比統(tǒng)計(jì)（表2）表明，1995—2019年耕地面積明顯減少，由89.25%減少至60.55%，城鎮(zhèn)化發(fā)展使得建設(shè)用地面積由10.18%增加至33.29%。隨著生態(tài)環(huán)境保護(hù)力度增加和退耕還林、植樹(shù)造林政策的實(shí)施，林地面積增加了5.54百分點(diǎn)。草地和水域面積占比變化不大。

從轉(zhuǎn)出數(shù)據(jù)分析（表3），耕地轉(zhuǎn)出面積最多，建設(shè)用地的轉(zhuǎn)出面積第二，主要復(fù)墾為耕地，其次為建設(shè)用地、林地、水域。從轉(zhuǎn)入特征分析，城鎮(zhèn)擴(kuò)張導(dǎo)致建設(shè)用地轉(zhuǎn)入面積最多，其次為林地、耕地、水域、草地，未利用地轉(zhuǎn)入面積最少。1995—2019年寧晉縣土地利用表現(xiàn)為建設(shè)用地大量擴(kuò)張，耕地面積迅速減少，林地、水域和草地面積不斷增加的變化特征。農(nóng)田生態(tài)系統(tǒng)為研究區(qū)主要生態(tài)系統(tǒng)類型。

2.2? 生態(tài)系統(tǒng)功能價(jià)值評(píng)價(jià)

從表4可以看出，1995—2019年寧晉縣生態(tài)系統(tǒng)服務(wù)總價(jià)值、生產(chǎn)功能和生態(tài)功能整體上呈先增加后減少的趨勢(shì)，氣候調(diào)節(jié)功能增加明顯，養(yǎng)分固持功能穩(wěn)定增長(zhǎng)，水源涵養(yǎng)功能明顯下降；氣體調(diào)節(jié)功能呈快速增加趨勢(shì)，與林地面積的顯著增加有關(guān)；污染物容納功能逐年降低。各類型生態(tài)系統(tǒng)服務(wù)價(jià)值平均值從大到小依次為農(nóng)田>草地>森林>淡水（表5）。

2.3? 生態(tài)系統(tǒng)功能權(quán)衡/協(xié)同關(guān)系時(shí)空變化

權(quán)衡系數(shù)計(jì)算結(jié)果（表6）顯示，1995—2003年生態(tài)系統(tǒng)多功能性和各功能間為協(xié)同關(guān)系，其中污染物容納功能權(quán)衡系數(shù)大于0，農(nóng)田生態(tài)環(huán)境壓力增大。2003—2011年，生態(tài)系統(tǒng)多功能性、生態(tài)功能和水源涵養(yǎng)功能與生產(chǎn)功能、氣候調(diào)節(jié)、養(yǎng)分固持、氣體調(diào)節(jié)功能之間呈權(quán)衡關(guān)系；生產(chǎn)功能與生態(tài)功能呈此消彼長(zhǎng)的權(quán)衡關(guān)系；水源涵養(yǎng)功能大幅降低。2011—2019年，生態(tài)系統(tǒng)多功能性、生產(chǎn)功能、生態(tài)功能、養(yǎng)分固持和水源涵養(yǎng)功能與其他功能呈權(quán)衡關(guān)系；環(huán)境污染程度增加，農(nóng)田生態(tài)環(huán)境壓力持續(xù)增大；生態(tài)功能與生產(chǎn)功能呈協(xié)同下降趨勢(shì)，生產(chǎn)功能急劇銳減；生態(tài)功能下降速度較上一時(shí)期緩慢，說(shuō)明近年退耕還林、大面積植樹(shù)造林等生態(tài)建設(shè)為遏制生態(tài)功能的下降發(fā)揮了一定作用。

利用ArcGIS根據(jù)土地利用和生態(tài)系統(tǒng)價(jià)值數(shù)據(jù)繪制生態(tài)系統(tǒng)功能損益圖，如圖2～3所示。從圖2～3可以看出，1995—2003年，生態(tài)系統(tǒng)多功能性以及各功能增益和減損面積變化較小。2003—2011年，氣候調(diào)節(jié)功能增益面積顯著增加，減損面積較少；其他功能增益減損面積較少。2011—2019年，生態(tài)功能增益和減損面積變化最為顯著，氣體調(diào)節(jié)功能損益面積變化較小，其他功能增益和減損面積變化一致。

從生態(tài)系統(tǒng)功能時(shí)間權(quán)衡與協(xié)同變化趨勢(shì)（圖4）可以看出，1995—2019年生態(tài)系統(tǒng)多功能性和生產(chǎn)功能的增長(zhǎng)速度明顯減慢，根據(jù)發(fā)展趨勢(shì)未來(lái)可能出現(xiàn)減損情況；污染物容納功能持續(xù)減損，生態(tài)功能總體減損有所減緩；林地的增加減緩了生態(tài)功能減損速度，水源涵養(yǎng)功能下降是導(dǎo)致區(qū)域生態(tài)功能下降的主要原因。

由生態(tài)系統(tǒng)功能損益圖（圖5）可知：多功能性、生產(chǎn)功能、生態(tài)功能、氣候調(diào)節(jié)、養(yǎng)分固持功能損益的空間變化規(guī)律一致，隨時(shí)間減損面積有較大增加；早期水源涵養(yǎng)功能大面積增加，近期呈大面積減損；氣體調(diào)節(jié)功能保持變化較?。晃廴疚锶菁{功能減損和增益均有較大面積，增益空間為耕地轉(zhuǎn)出面積。1995—2019年區(qū)域生態(tài)功能呈大面積增益而權(quán)衡系數(shù)為負(fù)，表明生態(tài)功能時(shí)空權(quán)衡關(guān)系不一致；農(nóng)田生態(tài)系統(tǒng)生態(tài)功能單位價(jià)值有所提升且面積占比最大，但該區(qū)域生態(tài)功能價(jià)值卻整體下降，進(jìn)一步說(shuō)明耕地大面積轉(zhuǎn)入建設(shè)用地導(dǎo)致生態(tài)系統(tǒng)價(jià)值流失。農(nóng)田生態(tài)系統(tǒng)多功能性呈現(xiàn)快速增長(zhǎng)—緩慢增長(zhǎng)—逐漸下降的變化過(guò)程；生產(chǎn)功能增速變緩，與生態(tài)功能由協(xié)同關(guān)系轉(zhuǎn)變?yōu)闄?quán)衡關(guān)系。草地生態(tài)系統(tǒng)

多功能性逐漸增加，其生產(chǎn)與生態(tài)功能呈協(xié)同關(guān)系。淡水生態(tài)系統(tǒng)多功能性逐漸增加，其生態(tài)功能與水源涵養(yǎng)為協(xié)同關(guān)系。森林生態(tài)系統(tǒng)各功能間保持協(xié)同增長(zhǎng)。

3? 討論

土地利用變化對(duì)生態(tài)系統(tǒng)價(jià)值的影響十分顯著，這與多數(shù)研究結(jié)果一致［17，45-47］。研究區(qū)耕地大量流向建設(shè)用地導(dǎo)致生態(tài)系統(tǒng)多功能性降低。增加林地類型結(jié)構(gòu)有利于緩解區(qū)域生態(tài)環(huán)境壓力和提升整體生態(tài)功能。有研究顯示，林地面積的增長(zhǎng)對(duì)生態(tài)系統(tǒng)的積極影響一定程度抵消了建設(shè)用地面積的加劇帶來(lái)的消極影響［10］。有研究表明氣候?qū)烫挤?wù)貢獻(xiàn)較大［48］。

對(duì)生態(tài)系統(tǒng)服務(wù)與功能價(jià)值評(píng)估采用計(jì)算經(jīng)濟(jì)價(jià)值的方法，結(jié)合文獻(xiàn)查閱采用普查數(shù)據(jù)、系數(shù)修正等方法評(píng)估。此外經(jīng)典的評(píng)估方法如謝高地等［49-51］核算了生態(tài)系統(tǒng)服務(wù)單位面積價(jià)值量，部分文獻(xiàn)通過(guò)計(jì)算植物凈初級(jí)生產(chǎn)力評(píng)估生態(tài)系統(tǒng)功能［52-53］。對(duì)于消納污染物的能力評(píng)估，有研究通

過(guò)評(píng)估農(nóng)藥地膜等對(duì)農(nóng)田生態(tài)系統(tǒng)的污染造成的損失來(lái)表示［54］。結(jié)合相關(guān)研究［42］引入減產(chǎn)系數(shù)，對(duì)地膜污染造成的生態(tài)系統(tǒng)功能價(jià)值損失進(jìn)行評(píng)估，此外有文獻(xiàn)采用意愿調(diào)查法［55］進(jìn)行評(píng)估?；柿魇Ш娃r(nóng)藥使用的污染負(fù)價(jià)值評(píng)價(jià)根據(jù)污染的生產(chǎn)排污收費(fèi)計(jì)算。近年來(lái)對(duì)區(qū)域生態(tài)系統(tǒng)服務(wù)與功能權(quán)衡分析認(rèn)為生態(tài)系統(tǒng)功能權(quán)衡分析是多重的，應(yīng)考慮時(shí)間和空間尺度效應(yīng)［56-59］。

4? 結(jié)論

1995—2019年，研究區(qū)生態(tài)系統(tǒng)多功能性、生產(chǎn)功能和生態(tài)功能整體呈現(xiàn)先增加后減少趨勢(shì)；氣候調(diào)節(jié)、養(yǎng)分固持功能穩(wěn)定增加，氣體調(diào)節(jié)功能增加迅速；水源涵養(yǎng)功能明顯下降；污染物容納功能持續(xù)降低。農(nóng)田生態(tài)系統(tǒng)是研究區(qū)生態(tài)系統(tǒng)的主體，由耕地面積減少導(dǎo)致流失的生態(tài)系統(tǒng)服務(wù)最大，該地類生態(tài)系統(tǒng)價(jià)值主要流向建設(shè)用地。森林、草地生態(tài)系統(tǒng)對(duì)提升區(qū)域生態(tài)的系統(tǒng)價(jià)值具有重要意義。土地利用對(duì)生態(tài)系統(tǒng)多功能性有重要影響。研究期內(nèi)區(qū)域生態(tài)功能、水源涵養(yǎng)功能、生態(tài)系統(tǒng)多功能性和其他功能由協(xié)同關(guān)系變?yōu)闄?quán)衡關(guān)系。農(nóng)田生態(tài)功能的協(xié)同/權(quán)衡關(guān)系時(shí)空表現(xiàn)不一致。

參考文獻(xiàn)

［1］王曉萌，潘佩佩，王曉旭，等.基于土地利用的河北省生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系研究［J］.地理與地理信息科學(xué)，2021，37（1）：80-88.

［2］鞏杰，燕玲玲，徐彩仙，等.1997—2018 年生態(tài)系統(tǒng)服務(wù)研究熱點(diǎn)變化與動(dòng)向［J］.應(yīng)用生態(tài)學(xué)報(bào)，2019，30（10）：3265-3276.

［3］殷楠，王帥，劉焱序.生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估：研究進(jìn)展與展望［J］.生態(tài)學(xué)雜志，2021，40（1）：233-244.

［4］COSTANZA R，DARGE R，DE GROOT R，et al.The value of the worlds ecosystem services and natural capital［J］.Nature，1997，387（6630）：253-260.

［5］《中國(guó)生物多樣性國(guó)情研究報(bào)告》編寫組.中國(guó)生物多樣性國(guó)情研究報(bào)告［M］.北京：中國(guó)環(huán)境科學(xué)出版社，1998.

［6］吳攀，楊木壯，劉偉聰.2007—2015年?yáng)|莞市土地利用及其生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空演變［J］.水土保持通報(bào)，2020，40（1）：250-255.

［7］李曼，李園園，劉煥才.阿克蘇河流域中下游生態(tài)系統(tǒng)服務(wù)價(jià)值對(duì)土地利用變化的響應(yīng)［J］.內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，41（3）：33-37.

［8］耿紀(jì)文，帥紅，張英，等.洞庭湖生態(tài)經(jīng)濟(jì)區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空變化［J］.中南林業(yè)科技大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2020，14（4）：53-61，76.

［9］喬斌，祝存兄，曹曉云，等.格網(wǎng)尺度下青海瑪多縣土地利用及生態(tài)系統(tǒng)服務(wù)價(jià)值空間自相關(guān)分析［J］.應(yīng)用生態(tài)學(xué)報(bào)，2020，31（5）：1660-1672.

［10］封建民，郭玲霞，李曉華.漢中市土地利用時(shí)空變化及其對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.水土保持研究，2020，27（1）：275-282.

［11］李龍，吳大放，劉艷艷.國(guó)內(nèi)外土地利用與生態(tài)系統(tǒng)服務(wù)研究熱點(diǎn)與趨勢(shì)：基于Citespace計(jì)量分析［J］.水土保持研究，2020，27（5）：396-404.

［12］歐陽(yáng)曉，賀清云，朱翔.多情景下模擬城市群土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響：以長(zhǎng)株潭城市群為例［J］.經(jīng)濟(jì)地理，2020，40（1）：93-102.

［13］包妲，李杰，劉鋒，等.1995—2015年緬甸土地利用/覆被變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.生態(tài)學(xué)報(bào)，2021，41（17）：6960-6969.

［14］郭亞紅，阿布都熱合曼·哈力克，魏天寶，等.基于土地利用變化的和田地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值分析［J］.生態(tài)學(xué)報(bào)，2021，41（16）：6363-6372.

［15］黃金亭，曹艷萍，秦奮.基于土地利用/覆蓋變化模擬的黃河流域生態(tài)環(huán)境質(zhì)量分析［J］.河南大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，50（2）：127-138.

［16］何楠.基于土地利用變化的湖北省生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估研究［D］.武漢：華中師范大學(xué)，2020.

［17］雷金睿，陳宗鑄，陳小花，等.1980—2018 年海南島土地利用與生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化［J］.生態(tài)學(xué)報(bào)，2020，40（14）：4760-4773.

［18］郭年冬，陳召亞，李恒哲，等.基于土地利用變化下的縣域生態(tài)敏感性及灰色預(yù)測(cè)：以河北省平山縣為例［J］.水土保持研究2016，23（5）：229-234.

［19］丁夢(mèng)婷，裴鳳松，胡引翠，等.基于 LUCC 的江蘇省生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化特征研究［J］.生態(tài)學(xué)報(bào)，2020，40（19）：6801-6811.

［20］王川，劉春芳，烏亞汗，等.黃土丘陵區(qū)生態(tài)系統(tǒng)服務(wù)空間格局及權(quán)衡與協(xié)同關(guān)系：以榆中縣為例［J］.生態(tài)學(xué)雜志，2019，38（2）：521－531.

［21］李瑾璞.基于InVEST模型的土地利用變化與生態(tài)系統(tǒng)碳儲(chǔ)量研究：以雄安新區(qū)為例［D］.保定：河北農(nóng)業(yè)大學(xué)，2020.

［22］陳昳.基于InVEST模型的舟山群島生態(tài)系統(tǒng)服務(wù)能值研究［D］.舟山：浙江海洋大學(xué)，2020.

［23］嚴(yán)有龍，王軍，王金滿.基于生態(tài)系統(tǒng)服務(wù)的閩江流域生態(tài)補(bǔ)償閾值研究［J］.中國(guó)土地科學(xué)，2021，35（3）：97-106.

［24］宋昌素，歐陽(yáng)志云.面向生態(tài)效益評(píng)估的生態(tài)系統(tǒng)生產(chǎn)總值 GEP 核算研究：以青海省為例［J］.生態(tài)學(xué)報(bào)，2020，40（10）：3207-3217.

［25］梁友嘉，劉麗珺.生態(tài)系統(tǒng)服務(wù)與景觀格局集成研究綜述［J］.生態(tài)學(xué)報(bào)，2018，38（20）：7159-7167.

［26］劉玉斌，王曉利，侯西勇，等.四個(gè)時(shí)期和模擬情景下2025年黃河三角洲的土地利用格局和生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估［J］.濕地科學(xué)，2020，18（4）：424-436.

［27］楊肖杰.河南省沿黃地區(qū)土地利用變化與生態(tài)安全研究［D］.開(kāi)封：河南大學(xué)，2020.

［28］文雯，彭哲，吳新星，等.基于耕地多功能的農(nóng)區(qū)全域土地綜合整治研究［J］.安徽農(nóng)業(yè)科學(xué)，2020，48（19）：262-263，268.

［29］國(guó)家林業(yè)局.森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范：LY/T 1721—2008［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2008.

［30］劉治國(guó)，劉玉華，于清軍，等.河北省草地生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估［J］.河北師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2021，45（3）：304-313.

［31］唐麗清，邱爾發(fā)，韓玉麗，等.石家莊市森林生態(tài)服務(wù)價(jià)值評(píng)估及預(yù)測(cè)［J］.中國(guó)城市林業(yè)，2014，12（2）：5-9.

［32］趙同謙，歐陽(yáng)志云，賈良清，等.中國(guó)草地生態(tài)系統(tǒng)服務(wù)功能間接價(jià)值評(píng)價(jià)［J］.生態(tài)學(xué)報(bào)，2004，24（6）：1101-1110.

［33］劉洋洋，任涵玉，周榮磊，等.中國(guó)草地生態(tài)系統(tǒng)服務(wù)價(jià)值估算及其動(dòng)態(tài)分析［J］.草地學(xué)報(bào)，2021，29（7）：1522-1532.

［34］蘇琪琪.太行山低山區(qū) 6 種林分類型水源涵養(yǎng)功能評(píng)價(jià)［D］.保定：河北農(nóng)業(yè)大學(xué)，2018.

［35］田志會(huì)，劉瑞涵.北京市農(nóng)田生態(tài)系統(tǒng)碳匯及釋氧功能年際變化研究［J］.生態(tài)經(jīng)濟(jì)，2016，32（1）：68-71.

［36］戴志剛，魯劍巍，周先竹，等.中國(guó)農(nóng)作物秸稈養(yǎng)分資源現(xiàn)狀及利用方式［J］.湖北農(nóng)業(yè)科學(xué)，2013，52（1）：27-29.

［37］彭立群，張強(qiáng)，賀克斌.基于調(diào)查的中國(guó)秸稈露天焚燒污染物排放清單［J］.環(huán)境科學(xué)研究，2016，29（8）：1109-1118.

［38］宋大利，侯勝鵬，王秀斌，等.中國(guó)秸稈養(yǎng)分資源數(shù)量及替代化肥潛力［J］.植物營(yíng)養(yǎng)與肥料學(xué)報(bào)，2018，24（1）：1-21.

［39］顧旭彪，張晉國(guó)，蘇磊峰，等.小麥免耕播種土壤水分蒸發(fā)量試驗(yàn)［J］.河北農(nóng)機(jī)，2017（8）：28-29.

［40］封春華，李娜，劉連梅.城鎮(zhèn)雨排水與農(nóng)田排澇流量計(jì)算方法對(duì)比分析［J］.水科學(xué)與工程技術(shù)，2016（2）：68-70.

［41］周傳豹.農(nóng)戶地膜使用行為與農(nóng)業(yè)面源污染問(wèn)題研究：以新疆地區(qū)為例［D］.烏魯木齊：新疆農(nóng)業(yè)大學(xué)，2014.

［42］黃星炯，陳仲清，劉香春.地膜殘留土壤對(duì)花生生育影響的研究［J］.中國(guó)油料，1993，5（3）：45-48.

［43］王亞潔.環(huán)境監(jiān)測(cè)中生化需氧量與化學(xué)需氧量的相關(guān)性研究［J］.河南科技，2022，41（6）：125-128.

［44］GONG S X，HODGSON J A，TSCHARNTKE T，et al.Biodiversity and yield trade.offs for organic farming［J］.Ecology letters，2022，25（7）：1699-1710.

［45］王軍，頓耀龍.土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響研究綜述［J］.長(zhǎng)江流域資源與環(huán)境，2015，24（5）：798-808.

［46］王雅，蒙吉軍.黑河中游土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響［J］.干旱區(qū)研究，2017，34（1）：200-207.

［47］王金偉.雙鴨山市土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響研究［D］.哈爾濱：東北農(nóng)業(yè)大學(xué)，2015.

［48］劉閔.1992—2015 年氣候變化和人類活動(dòng)對(duì)中國(guó)生態(tài)系統(tǒng)服務(wù)功能權(quán)衡及協(xié)同關(guān)系的影響［D］.貴陽(yáng)：貴州師范大學(xué)，2022.

［49］謝高地，張彩霞，張昌順，等.中國(guó)生態(tài)系統(tǒng)服務(wù)的價(jià)值［J］.資源科學(xué)，2015，37（9）：1740-1746.

［50］謝高地，張彩霞，張雷明，等.基于單位面積價(jià)值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法改進(jìn)［J］.自然資源學(xué)報(bào)，2015，30（8）：1243-1254.

［51］謝高地，魯春霞，成升魁.全球生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估研究進(jìn)展［J］.資源科學(xué)，2001，23（6）：5-9.

［52］劉苗苗.京津冀地區(qū)耕地生態(tài)系統(tǒng)服務(wù)研究［D］.石家莊：河北師范大學(xué)，2022.

［53］白琳紅，王衛(wèi)，張玉.河北省縣域農(nóng)田生態(tài)系統(tǒng)供給功能的健康評(píng)價(jià)［J］.生態(tài)學(xué)報(bào)，2012，32（7）：2033-2039．

［54］元媛，劉金銅，靳占忠.欒城縣農(nóng)田生態(tài)系統(tǒng)服務(wù)功能正負(fù)效應(yīng)綜合評(píng)價(jià)［J］.生態(tài)學(xué)雜志，2011，30（2）：2809-2814.

［55］楊志新.北京郊區(qū)農(nóng)田生態(tài)系統(tǒng)正負(fù)效應(yīng)價(jià)值的綜合評(píng)價(jià)研究［D］.北京：中國(guó)農(nóng)業(yè)大學(xué)，2006.

［56］戴路煒，唐海萍，張欽，等.北方農(nóng)牧交錯(cuò)帶多倫縣生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系研究［J］.生態(tài)學(xué)報(bào)，2020，40（9）：2863-2876.

［57］卓志清，興安，孫忠祥，等.東北旱作區(qū)農(nóng)業(yè)生態(tài)系統(tǒng)協(xié)同發(fā)展與權(quán)衡分析［J］.中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2018，26（6）：892-902.

［58］孫藝杰，任志遠(yuǎn)，郝夢(mèng)雅，等.黃土高原生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同時(shí)空變化及影響因素：以延安市為例［J］.生態(tài)學(xué)報(bào)，2019，39（10）：3443-3454.

［59］祝萍，劉鑫，鄭瑜晗，等.北方重點(diǎn)生態(tài)功能區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同［J］.生態(tài)學(xué)報(bào)，2020，40（23）：8694-8706.