生態(tài)系統(tǒng)服務(wù)與人類活動協(xié)調(diào)性時空演變及預(yù)測

摘要：探究并預(yù)測生態(tài)系統(tǒng)服務(wù)與人類活動之間關(guān)系的演變特征對于推動建設(shè)人與自然和諧共生的生態(tài)新格局意義重大。以陜西省黃河流域為例，通過價值當(dāng)量法、人類活動強度模型、PLUS模型和多種統(tǒng)計分析方法，對生態(tài)系統(tǒng)服務(wù)與人類活動以及二者之間的協(xié)調(diào)性進行時空演變特征分析和多情景預(yù)測。結(jié)果表明：1980～2020年陜西省黃河流域生態(tài)系統(tǒng)服務(wù)價值總體增加54.33億元，與2020年相比，2030年慣性發(fā)展和經(jīng)濟優(yōu)先模式下生態(tài)系統(tǒng)服務(wù)價值分別減少7.81億元、20.11億元，耕地保護和生態(tài)保護模式下分別增加0.27億元、8.24億元。1980～2020年人類活動強度呈現(xiàn)出“N”形變化特征，與2020年相比，經(jīng)濟優(yōu)先和慣性發(fā)展模式下人類活動強度增強，耕地保護和生態(tài)保護模式下有所減弱。研究區(qū)生態(tài)系統(tǒng)服務(wù)與人類活動呈現(xiàn)逐漸增強的負(fù)相關(guān)關(guān)系，二者協(xié)調(diào)性較差，2030年生態(tài)保護模式下二者的負(fù)相關(guān)作用最弱，經(jīng)濟優(yōu)先模式下最強。綜合來看，未來經(jīng)濟優(yōu)先模式下的負(fù)面影響最大，陜西省黃河流域應(yīng)優(yōu)先選擇生態(tài)保護型的發(fā)展模式。

關(guān) 鍵 詞：生態(tài)系統(tǒng)服務(wù)；人類活動；多情景模擬；協(xié)調(diào)性；黃河流域

中圖法分類號：X24；F301.2；X171.1

文獻(xiàn)標(biāo)志碼：ADOI：10.16232/j.cnki.1001-4179.2024.10.014

0 引 言

黃河流域位于中國北部地區(qū)，近半個世紀(jì)以來，隨著人口增加和社會經(jīng)濟快速發(fā)展，流域地區(qū)城市化水平迅速提升，人類活動對黃河流域本就脆弱的生態(tài)環(huán)境影響頗深，如何平衡流域生態(tài)保護與人類活動一直以來是有關(guān)黃河流域的核心議題^［1］。習(xí)近平總書記多次對黃河流域做出重要指示批示，在2019年9月18日黃河流域生態(tài)保護和高質(zhì)量發(fā)展座談會上強調(diào)“黃河流域是我國重要的生態(tài)屏障和重要的經(jīng)濟地帶，是打贏脫貧攻堅戰(zhàn)的重要區(qū)域，在我國經(jīng)濟社會發(fā)展和生態(tài)安全方面具有十分重要的地位?！?，明確指出“黃河流域生態(tài)保護和高質(zhì)量發(fā)展，同京津冀協(xié)同發(fā)展、長江經(jīng)濟帶發(fā)展、粵港澳大灣區(qū)建設(shè)、長三角一體化發(fā)展一樣，是重大國家戰(zhàn)略”^［2］。黨中央、國務(wù)院及流域內(nèi)相關(guān)政府制定并落實了《黃河流域生態(tài)保護和高質(zhì)量發(fā)展規(guī)劃綱要》、省市《黃河流域生態(tài)保護和高質(zhì)量發(fā)展規(guī)劃》等政策規(guī)劃，關(guān)鍵目標(biāo)之一是維護區(qū)域內(nèi)生態(tài)系統(tǒng)功能，協(xié)調(diào)生態(tài)保護與人類活動，改善人民生態(tài)福祉。在此背景下，開展黃河流域生態(tài)系統(tǒng)服務(wù)價值與人類活動協(xié)調(diào)發(fā)展評價與預(yù)測的相關(guān)主題研究既符合黨中央、國務(wù)院對黃河流域生態(tài)保護和高質(zhì)量發(fā)展的戰(zhàn)略定位，也符合建設(shè)造福人民幸福河的群眾期盼，同時對于推動黃河流域可持續(xù)發(fā)展意義重大。

生態(tài)系統(tǒng)服務(wù)評估與人類活動管理的研究與實踐活動已在全球范圍內(nèi)得到了廣泛關(guān)注，眾多學(xué)者圍繞生態(tài)系統(tǒng)服務(wù)概念內(nèi)涵、量化方法以及時空特征、權(quán)衡協(xié)同、供需匹配及其相關(guān)的驅(qū)動機制^［3-8］等方面展開了大量討論，并取得了豐碩的研究成果。在此過程中，通過耦合人類活動以順應(yīng)生態(tài)管理實踐的研究議題逐漸成為了當(dāng)前人地關(guān)系研究的熱點^［9-10］，但多局限于針對“過去時”狀態(tài)的揭示，而對“將來時”生態(tài)系統(tǒng)服務(wù)與人類活動協(xié)調(diào)特征的預(yù)測研究仍有待豐富。因此，本研究以陜西省黃河流域為例，以提供生態(tài)系統(tǒng)服務(wù)和承載人類活動的土地為切入點，結(jié)合價值當(dāng)量法、人類活動強度模型、PLUS模型等方法和多種統(tǒng)計分析手段，揭示生態(tài)系統(tǒng)服務(wù)價值與人類活動強度的時空演變規(guī)律及未來發(fā)展趨勢，并分析和預(yù)測二者之間的協(xié)調(diào)特征，以期豐富生態(tài)系統(tǒng)服務(wù)與人類活動相互作用主題研究的內(nèi)容，為促進黃河流域生態(tài)保護與人類活動協(xié)調(diào)發(fā)展提供前瞻性的決策支撐。

1 研究區(qū)與數(shù)據(jù)來源

1.1 研究區(qū)概況

陜西省黃河流域位于黃河流域中部（圖1），覆蓋陜西省77個縣（市、區(qū)），面積12.9萬km²，約占黃河流域總面積的16.3%。陜西省黃河流域海拔介于316～3 742 m之間，大致呈現(xiàn)東南低西北高的特點，屬于半濕潤半干旱氣候區(qū)，多年平均降雨量為325～878 mm，表現(xiàn)出自南向北遞減的空間特征。該地區(qū)土地利用類型以耕地、草地和林地為主，土壤類型以黃土為主，土質(zhì)疏松、易被風(fēng)和水侵蝕，是中國水土流失較為嚴(yán)重的地區(qū)之一^［11］。陜西省黃河流域生態(tài)基礎(chǔ)薄弱，地形、氣候等要素具有明顯的空間異質(zhì)性，長期以來存在著干旱缺水、土地退化、水土流失等一系列典型生態(tài)問題^［12］。與此同時，在退耕還林還草工程實施、城市建設(shè)用地擴張等人類活動的交叉影響下，該地區(qū)社會-生態(tài)系統(tǒng)發(fā)生了顯著變化^［11］，有必要對其生態(tài)系統(tǒng)服務(wù)及人類活動進行量化和預(yù)測，并進一步揭示兩者協(xié)調(diào)演變規(guī)律，以更好地服務(wù)于黃河流域生態(tài)保護和高質(zhì)量發(fā)展。

1.2 數(shù)據(jù)基礎(chǔ)

本研究所需數(shù)據(jù)及其來源如表1所列。

2 研究方法

2.1 生態(tài)系統(tǒng)服務(wù)量化

基于謝高地等^［5］提出的價值當(dāng)量法計算陜西省黃河流域的生態(tài)系統(tǒng)服務(wù)價值。詳細(xì)過程如下：

（1）單位生態(tài)系統(tǒng)服務(wù)當(dāng)量經(jīng)濟價值的計算。結(jié)合公式（1）對陜西省黃河流域單位生態(tài)系統(tǒng)服務(wù)當(dāng)量的經(jīng)濟價值進行計算。

E=17∑Q×P×S∑S（1）

式中：E為單位生態(tài)系統(tǒng)服務(wù)當(dāng)量的經(jīng)濟價值；l為糧食作物種類；Q為l類糧食作物單產(chǎn)；P為l類糧食作物單價；S為l類糧食作物面積。

（2）生態(tài)系統(tǒng)服務(wù)價值系數(shù)的修訂。采用植被凈初級生產(chǎn)力（NPP）對生態(tài)系統(tǒng)服務(wù)價值系數(shù)進行本地化修正^［13］。計算公式為

VC=E×H×σ（3）

式中：σ為本地化修正系數(shù)；NPP、NPP分別為陜西省黃河流域和全國的NPP值；VC為i類用地的第f項生態(tài)系統(tǒng)服務(wù)價值系數(shù)；H為i類用地的第f項生態(tài)系統(tǒng)服務(wù)當(dāng)量。

（3）生態(tài)系統(tǒng)服務(wù)價值的量化。采用公式（4）計算陜西省黃河流域的生態(tài)系統(tǒng)服務(wù)價值。

ESV=∑∑S×VC（4）

式中：ESV為生態(tài)系統(tǒng)服務(wù)價值；S為i類用地類型的面積。由于研究區(qū)不同縣域單元面積的差異，為保證生態(tài)系統(tǒng)服務(wù)價值計算結(jié)果在不同縣域間的可比性，結(jié)合不同縣域面積，求取地均生態(tài)系統(tǒng)服務(wù)價值，用于相關(guān)分析中。

（4）生態(tài)系統(tǒng)服務(wù)價值結(jié)果的檢驗。以經(jīng)濟學(xué)中常用的彈性系數(shù)計算價值系數(shù)的敏感性，以確定陜西省黃河流域生態(tài)系統(tǒng)服務(wù)價值對價值系數(shù)變化調(diào)整的依賴程度，以檢驗結(jié)果的合理性。計算公式如下：

SC=ESV-ESV/ESVVC-VC/VC（5）

式中：SC為敏感性系數(shù)；ESV和ESV代表初始和調(diào)整后的生態(tài)系統(tǒng)服務(wù)價值；VC、VC表示i類用地類型的初始和調(diào)整后的價值系數(shù)。本研究將各用地類型的價值系數(shù)分別調(diào)整50%，來衡量生態(tài)系統(tǒng)服務(wù)價值的變化情況。若SC>1，說明生態(tài)系統(tǒng)服務(wù)價值對VC是富有彈性的，即可信度較低；若SC<1，則說明生態(tài)系統(tǒng)服務(wù)價值對VC是缺乏彈性的，計算結(jié)果是可信的^［13］。

經(jīng)計算，得到陜西省黃河流域不同土地利用類型的敏感性系數(shù)，如表2所列?？梢园l(fā)現(xiàn)，在1980～2020年以及2030年不同情景模式下所有土地利用類型的敏感性系數(shù)均遠(yuǎn)小于1，說明本研究生態(tài)系統(tǒng)服務(wù)價值的計算結(jié)果符合研究區(qū)實際。

2.2 人類活動強度量化

人類活動通過改變土地利用類型，改變生態(tài)系統(tǒng)的結(jié)構(gòu)和過程，進而與生態(tài)系統(tǒng)服務(wù)產(chǎn)生關(guān)聯(lián)。因此，本研究采用基于土地利用類型的人類活動強度模型^［14］來量化陜西省黃河流域的人類活動強度。公式如下：

HAI=∑S×H∑S（6）

式中：HAI為人類活動強度；S為i類用地的面積；H為i類用地的人類活動強度指數(shù)，以Delphi法、Lohani清單法和Leopold矩陣法的平均值表示^［15］（表3）。

2.3 不同情景模式模擬

（1）情景模式設(shè)置。本研究結(jié)合陜西省國土空間總體規(guī)劃成果、生態(tài)恢復(fù)背景以及區(qū)域城市化發(fā)展特征，共設(shè)置了以下4種情景模式，以模擬2030年不同情景模式下的土地利用格局。① 慣性發(fā)展模式：考慮城市化進程和現(xiàn)有的土地政策，延續(xù)歷史的土地利用發(fā)展趨勢，保持各土地利用類型轉(zhuǎn)移概率不變，模擬2030年的土地利用情況。② 經(jīng)濟優(yōu)先模式：陜西省黃河流域正處于城市化和經(jīng)濟快速發(fā)展和轉(zhuǎn)型階段，該模式強調(diào)以經(jīng)濟發(fā)展和城市建設(shè)為核心，依據(jù)城市建設(shè)用地擴張需求，基于慣性發(fā)展模式，通過限制建設(shè)用地向耕地、水域、林地的轉(zhuǎn)移概率，增加耕地、林地、草地、水域、未利用地向建設(shè)用地的轉(zhuǎn)移概率，實現(xiàn)對2030年土地利用格局的模擬。③ 耕地保護模式：耕地是糧食安全的載體，以耕地保護為首要目的，基于慣性發(fā)展模式，一方面增加未利用地向耕地的轉(zhuǎn)移概率，另一方面嚴(yán)格控制耕地向建設(shè)用地、草地、水域的轉(zhuǎn)移，實現(xiàn)對2030年土地利用格局的模擬。④ 生態(tài)保護模式：在當(dāng)前的國土空間規(guī)劃背景下，生態(tài)保護模式基于慣性發(fā)展模式，以生態(tài)環(huán)境保護為首要目的，限制城鎮(zhèn)化發(fā)展，突出嚴(yán)格管控生態(tài)用地的大規(guī)模開發(fā)，通過嚴(yán)格限制生態(tài)用地的轉(zhuǎn)出概率、增強未利用地向林地和草地的轉(zhuǎn)移，實現(xiàn)對2030年土地利用格局的模擬。

（2）模型模擬。PLUS模型克服了CA-Markov、FLUS等模型在轉(zhuǎn)化規(guī)則挖掘策略和景觀動態(tài)變化模擬策略等方面存在的限制，同時還克服了CLUE-S模型模擬區(qū)域受限、無法考慮政策驅(qū)動和引導(dǎo)作用的問題。因此，本研究采用PLUS模型對陜西省黃河流域2030年的土地利用格局進行模擬，模型詳細(xì)介紹見文獻(xiàn)［16］。為更準(zhǔn)確地模擬不同情景模式下陜西省黃河流域的土地利用情況，綜合前人研究^［16-17］與研究區(qū)實際情況，從自然環(huán)境、社會經(jīng)濟和可達(dá)性3個角度選取海拔、坡度、年平均降雨、年平均溫度、人口密度、經(jīng)濟密度、到居民點距離、到政府距離、到水域距離、到道路距離、到鐵路距離等11個驅(qū)動因子。

（3）模型結(jié)果驗證。采用PLUS模型對陜西省黃河流域2020年土地利用格局進行預(yù)測，并與2020年實際土地利用數(shù)據(jù)進行對比驗證。經(jīng)計算，Kappa系數(shù)為0.875，總體精度為0.913，模擬精度較高，結(jié)果可信，說明PLUS模型在陜西省黃河流域未來土地利用預(yù)測中具有良好的適用性。

2.4 統(tǒng)計分析方法

2.4.1 地理集中度

地理集中度可以在綜合考慮區(qū)域面積的前提下，刻畫要素的空間集聚狀態(tài)^［18］。因此，本研究采用地理集中度描述生態(tài)系統(tǒng)服務(wù)與人類活動強度的空間分布特征。公式為

R=ESV/∑ESVS/∑S（7）

R=HAI/∑HAIS/∑S（8）

式中：R與R分別為縣域j的生態(tài)地理集中度和人類活動地理集中度，其值越大說明該縣域生態(tài)系統(tǒng)服務(wù)價值或人類活動強度越集中，即生態(tài)系統(tǒng)服務(wù)價值或人類活動強度越高；ESV、HAI、S分別為縣域j的生態(tài)系統(tǒng)服務(wù)價值、人類活動強度和面積。

2.4.2 相關(guān)性分析

采用Spearman相關(guān)系數(shù)檢驗生態(tài)系統(tǒng)服務(wù)與人類活動之間的整體相互作用關(guān)系，公式如下：

ρ=1-6×∑kq=1d²k×（k²-1）（9）

式中：ρ為Spearman相關(guān)系數(shù)；d為要素之間同一組樣本的位次差；k為樣本數(shù)量。

2.4.3 空間關(guān)聯(lián)性分析

采用雙變量全局空間自相關(guān)分析揭示生態(tài)系統(tǒng)服務(wù)與人類活動在空間上的關(guān)聯(lián)特征，公式如下：

I=∑∑Wx-xy-yS²∑∑W（10）

式中：I為雙變量全局莫蘭指數(shù)，I<1表示空間負(fù)相關(guān)，I>1表示空間正相關(guān)，I=0表示不存在空間相關(guān)性；W為空間權(quán)重矩陣；x、y分別為地均生態(tài)系統(tǒng)服務(wù)價值、人類活動強度的觀測值；S²為所有樣本的方差。

2.4.4 協(xié)調(diào)性分析

采用一致性指數(shù)衡量生態(tài)系統(tǒng)服務(wù)與人類活動之間協(xié)調(diào)性的空間格局，可以用生態(tài)地理集中度與人類活動地理集中度的比值來表示^［18］。公式如下：

CPE=RR（11）

式中：CPE為一致性指數(shù)，其值越接近于1說明生態(tài)系統(tǒng)服務(wù)價值與人類活動強度的空間演變趨勢越一致，即生態(tài)系統(tǒng)服務(wù)與人類活動越協(xié)調(diào)，其值小于1說明生態(tài)系統(tǒng)服務(wù)價值演變滯后于人類活動強度，其值大于1說明生態(tài)系統(tǒng)服務(wù)價值演變超前于人類活動強度。

3 結(jié)果與分析

3.1 生態(tài)系統(tǒng)服務(wù)價值時空演變特征及預(yù)測

3.1.1 生態(tài)系統(tǒng)服務(wù)價值時間變化特征

陜西省黃河流域生態(tài)系統(tǒng)服務(wù)價值如圖2所示?？梢钥闯?，1980～2020年，研究區(qū)生態(tài)系統(tǒng)服務(wù)價值呈現(xiàn)出倒“N”形變化趨勢，先由1980年的2 987.64億元減少至1990年的2 983.93億元，隨后持續(xù)增加并于2010年達(dá)到研究期內(nèi)最大值3 053.46億元，最后又減少至2020年的3 041.97億元，研究期內(nèi)總體增加了54.33億元。2030年不同情景模式下的生態(tài)系統(tǒng)服務(wù)價值較2020年出現(xiàn)了不同方向的變化，其中慣性發(fā)展模式和經(jīng)濟優(yōu)先模式下生態(tài)系統(tǒng)服務(wù)價值出現(xiàn)減少，較2020年分別減少了7.81億元、20.11億元，而耕地保護模式和生態(tài)保護模式下生態(tài)系統(tǒng)服務(wù)價值出現(xiàn)增加，較2020年分別增加了0.27億元和8.24億元。總體來說，生態(tài)保護模式下的生態(tài)系統(tǒng)服務(wù)價值最高，經(jīng)濟優(yōu)先模式下的最低。

3.1.2 生態(tài)系統(tǒng)服務(wù)價值空間變化特征

陜西省黃河流域生態(tài)地理集中度的空間分布如圖3所示。從空間格局的時間變化上來說，不同年份以及不同情景模式下生態(tài)地理集中度存在著相似的空間分布特征，這很大程度上與生態(tài)系統(tǒng)的結(jié)構(gòu)和過程有關(guān)。一般來說，生態(tài)系統(tǒng)的結(jié)構(gòu)和過程決定著生態(tài)系統(tǒng)服務(wù)的供給，陜西省黃河流域不同時期同一縣域的生態(tài)系統(tǒng)結(jié)構(gòu)和過程較為相像，進而導(dǎo)致生態(tài)地理集中度隨時間表現(xiàn)出相似的空間格局。從空間分布上來看，研究期內(nèi)陜西省黃河流域最南部和中部地區(qū)的生態(tài)地理集中度較大，這里地形起伏較大，用地類型以林地和草地為主，生態(tài)環(huán)境相對較好，生態(tài)系統(tǒng)服務(wù)價值整體較高；相反，在研究區(qū)北部地區(qū)和中南部關(guān)中平原地區(qū)的生態(tài)地理集中度較小，其中北部地區(qū)未利用地分布廣泛，生態(tài)基礎(chǔ)較為薄弱，而關(guān)中平原地區(qū)地勢平坦，適宜耕種和居住，用地類型以耕地和建設(shè)用地為主，人類活動高強度干擾背景下的生態(tài)系統(tǒng)服務(wù)價值相對較低。

3.2 人類活動強度時空演變特征及預(yù)測

3.2.1 人類活動強度時間變化特征

圖4顯示了陜西省黃河流域人類活動強度變化情況。可以看出，1980～2020年，研究區(qū)人類活動強度呈現(xiàn)出“N”形變化特征，由1980年的0.324 5持續(xù)增加至2000年的0.327 4，隨后減小至2010年的0.319 5，最后又增加至2020年的0.322 0，以2000年人類活動強度最強，2020年人類活動強度最弱。在2030年的4種不同情景模式下，以經(jīng)濟優(yōu)先模式的人類活動強度最大（0.325 5），生態(tài)保護模式下的人類活動強度最?。?.321 2）。與2020年相比，經(jīng)濟優(yōu)先模式和慣性發(fā)展模式下的人類活動強度出現(xiàn)增強，而耕地保護模式和生態(tài)保護模式下的人類活動強度有所減弱。

3.2.2 人類活動強度空間變化特征

圖5顯示了陜西省黃河流域人類活動地理集中度的空間分布情況。從人類活動地理集中度空間格局的時間變化上來看，不同年份和不同情景模式下人類活動地理集中度在空間上存在著相似的分布特征。在空間分布上，陜西省黃河流域人類活動地理集中度的空間格局總體與生態(tài)地理集中度相反。研究區(qū)關(guān)中平原地區(qū)的人類活動地理集中度最高，這里人口密集，社會經(jīng)濟發(fā)展迅速，人類活動強度整體較高；相反，在研究區(qū)北部和中部地區(qū)的人類活動地理集中度較低，尤其是榆林市西北部地區(qū)的人類活動地理集中度最低，這里氣候干旱、溝壑縱橫，用地類型以未利用地為主，人口分布稀疏，人類活動強度最小。

3.3 生態(tài)系統(tǒng)服務(wù)價值與人類活動強度協(xié)調(diào)性及預(yù)測

3.3.1 相關(guān)性分析

用Spearman相關(guān)系數(shù)對地均生態(tài)系統(tǒng)服務(wù)價值與人類活動強度進行相關(guān)性分析，結(jié)果如表4所列。1980～2020年二者之間存在著高度的負(fù)相關(guān)關(guān)系，二者之間不協(xié)調(diào)性非常明顯，且隨時間呈現(xiàn)出進一步惡化趨勢。2030年不同情景模式下，耕地保護模式與生態(tài)保護模式下二者間的不協(xié)調(diào)性出現(xiàn)改善，而慣性發(fā)展模式和經(jīng)濟優(yōu)先模式下二者間的不協(xié)調(diào)性與2020年相比有所惡化，需進一步探討二者之間協(xié)調(diào)發(fā)展的時空特征。

3.3.2 空間關(guān)聯(lián)特征

利用GeoDa軟件得到了陜西省黃河流域1980～2020年以及不同情景模式下地均生態(tài)系統(tǒng)服務(wù)價值與人類活動強度的雙變量全局莫蘭指數(shù)，如表4所列?？梢园l(fā)現(xiàn)，地均生態(tài)系統(tǒng)服務(wù)價值與人類活動強度在研究期內(nèi)呈現(xiàn)出顯著的空間負(fù)關(guān)聯(lián)，且這種空間負(fù)關(guān)聯(lián)在1980～2020年隨時間呈現(xiàn)出增強趨勢。2030年，慣性發(fā)展模式和經(jīng)濟優(yōu)先模式下的空間負(fù)關(guān)聯(lián)較2020年有所增強，而耕地保護模式與生態(tài)保護模式下的空間負(fù)關(guān)聯(lián)與2020年相比有所減弱。

3.3.3 協(xié)調(diào)發(fā)展特征

在Spearman相關(guān)性分析和雙變量空間自相關(guān)分析的基礎(chǔ)上，通過引入一致性指數(shù)進一步厘清陜西省黃河流域生態(tài)系統(tǒng)服務(wù)與人類活動在各縣域內(nèi)部的一致性特征及演變趨勢，以闡釋研究區(qū)生態(tài)系統(tǒng)服務(wù)與人類活動的協(xié)調(diào)性（圖6）。1980～2020年，陜西省黃河流域生態(tài)系統(tǒng)服務(wù)與人類活動的協(xié)調(diào)性較差，二者地理集中度的一致性指數(shù)在空間分布上多屬于大于1.5或小于0.5的狀態(tài)，說明研究區(qū)生態(tài)系統(tǒng)服務(wù)明顯超前或滯后于人類活動。從一致性指數(shù)空間格局的時間變化特征上來看，一致性指數(shù)在不同年份和不同情景模式下的空間格局較為相似，說明生態(tài)系統(tǒng)服務(wù)與人類活動的關(guān)系在陜西省黃河流域各縣域內(nèi)部隨時間比較穩(wěn)定。從一致性指數(shù)空間格局上來說，一致性指數(shù)在陜北地區(qū)明顯大于1，即生態(tài)系統(tǒng)服務(wù)價值演變超前于人類活動強度。這里生態(tài)基礎(chǔ)薄弱，是中國生態(tài)修復(fù)的核心區(qū)域之一，生態(tài)系統(tǒng)積極向好發(fā)展；與此同時，該地區(qū)屬于溝壑區(qū)，人口分布稀疏，經(jīng)濟活動水平較低，人類活動強度低，進而導(dǎo)致了生態(tài)系統(tǒng)服務(wù)演變超前于人類活動的局面。一致性指數(shù)明顯小于1的縣域主要位于關(guān)中地區(qū)，這里地勢平坦，人口聚集，經(jīng)濟活動活躍，人類活動強度高而生態(tài)系統(tǒng)服務(wù)價值低，生態(tài)系統(tǒng)服務(wù)明顯滯后于人類活動。

4 討 論

陜西省黃河流域作為全國水土流失最為嚴(yán)重的地區(qū)之一，自20世紀(jì)90年代以來進行了一系列生態(tài)建設(shè)活動^［19］，這也導(dǎo)致了1990年以后該地區(qū)生態(tài)系統(tǒng)服務(wù)價值出現(xiàn)明顯提升，特別是2000年以后，隨著退耕還林還草工程在全區(qū)范圍內(nèi)的大規(guī)模實施，該地區(qū)生態(tài)系統(tǒng)服務(wù)價值進一步提升。但在2010～2020年，該地區(qū)生態(tài)系統(tǒng)服務(wù)價值有所減少，究其原因，盡管該時段生態(tài)修復(fù)仍在進行，但城市建設(shè)用地擴張速度明顯快于生態(tài)空間增加幅度^［20］，生態(tài)修復(fù)所帶來的生態(tài)效益小于同期城市空間擴張帶來的損失。本研究還發(fā)現(xiàn)生態(tài)系統(tǒng)服務(wù)與人類活動之間的權(quán)衡作用，即人類活動增強會導(dǎo)致生態(tài)系統(tǒng)服務(wù)的退化，這與其他研究的結(jié)果一致^{［14-15，21］}。

黨的二十大報告強調(diào)：“必須牢固樹立和踐行綠水青山就是金山銀山的理念，站在人與自然和諧共生的高度謀劃發(fā)展?！标兾魇↑S河流域作為黃河流域生態(tài)保護和高質(zhì)量發(fā)展的重點區(qū)域，近40 a來社會經(jīng)濟快速發(fā)展，頻繁的人類活動所引發(fā)的土地利用結(jié)構(gòu)劇烈變化與生態(tài)系統(tǒng)服務(wù)退化等問題^{［19，22］}，是本研究發(fā)現(xiàn)生態(tài)系統(tǒng)服務(wù)與人類活動強度之間的負(fù)相關(guān)程度隨時間推移不斷加劇的原因。但是，不同情景模式下生態(tài)系統(tǒng)服務(wù)與人類活動強度以及二者之間時空關(guān)聯(lián)信息顯示，堅持生態(tài)優(yōu)先原則的人類活動可緩解二者之間的權(quán)衡強度。因此，結(jié)合研究結(jié)果提出以下建議：① 在人類活動規(guī)劃方面，應(yīng)嚴(yán)格防止城市盲目擴張侵吞自然生態(tài)空間，合理劃定城市發(fā)展邊界，避免因人類過度開發(fā)活動帶來的生態(tài)空間破碎，集約節(jié)約利用城市現(xiàn)有建設(shè)空間。② 在協(xié)調(diào)人與自然方面，對于生態(tài)系統(tǒng)服務(wù)價值演變超前于人類活動強度的縣域而言，未來應(yīng)繼續(xù)推進生態(tài)修復(fù)，不斷改善生態(tài)環(huán)境本底，在維護生態(tài)系統(tǒng)完整性的同時，可通過提升基本公共服務(wù)水平等以增強縣域人口吸引力。而對于生態(tài)系統(tǒng)服務(wù)價值水平明顯滯后于人類活動的縣域，人與自然矛盾突出，人類活動復(fù)雜且強度高，未來應(yīng)著重控制人口規(guī)模和不合理人類活動，界定人類活動影響閾值，大力倡導(dǎo)綠色低碳發(fā)展，同時可通過加快城市公園、綠地建設(shè)等以有效改善縣域生態(tài)環(huán)境。

需要指出的是，本研究僅對陜西省黃河流域進行了案例分析，由于黃河流域地域遼闊，流域內(nèi)自然環(huán)境條件和社會經(jīng)濟發(fā)展水平等具有巨大的時空異質(zhì)性，未來應(yīng)加強黃河流域其他脆弱區(qū)如黃河源區(qū)、黃河三角洲等區(qū)域的相關(guān)研究，以更全面系統(tǒng)地掌握黃河流域的內(nèi)部差異。此外，本研究采用價值當(dāng)量法，并考慮了研究區(qū)生態(tài)實際對其進行了修正，但由于修正方法和參數(shù)的差異可能會導(dǎo)致生態(tài)系統(tǒng)服務(wù)價值量化結(jié)果與實際存有差異；而且，研究借助的價值當(dāng)量法中將建設(shè)用地的價值系數(shù)假定為0，實際上建設(shè)用地可能會有負(fù)的價值系數(shù)^［23］。因此，如何實現(xiàn)生態(tài)系統(tǒng)服務(wù)價值更為精準(zhǔn)的量化仍有深化空間。最后，本研究尚未深入探討生態(tài)系統(tǒng)服務(wù)與人類活動協(xié)調(diào)演變的驅(qū)動因素和機制，未來可加強背后影響因素的定量分析，以提供更具針對性和操作性的人與自然協(xié)調(diào)方案。

5 結(jié) 論

（1）1980～2020年陜西省黃河流域生態(tài)系統(tǒng)服務(wù)價值呈倒“N”形變化趨勢，總體增加54.33億元；空間上研究區(qū)最南部和中部地區(qū)的生態(tài)系統(tǒng)服務(wù)價值較高，而在研究區(qū)北部地區(qū)和中南部關(guān)中平原地區(qū)的價值較低；2030年，生態(tài)保護模式下的生態(tài)系統(tǒng)服務(wù)價值最高，而經(jīng)濟優(yōu)先模式下最低。

（2）1980～2020年陜西省黃河流域人類活動強度呈“N”形變化特征，2000年人類活動強度最強，2020年最弱；空間上關(guān)中平原地區(qū)的人類活動最強，而在研究區(qū)北部和中部地區(qū)的人類活動較弱，榆林市西北部地區(qū)的人類活動最弱；2030年，生態(tài)保護模式下的人類活動最弱，而經(jīng)濟優(yōu)先模式下最強。

（3）1980～2020年陜西省黃河流域生態(tài)系統(tǒng)服務(wù)與人類活動呈現(xiàn)逐漸增強的負(fù)相關(guān)關(guān)系，二者協(xié)調(diào)性較差；空間上陜北地區(qū)的生態(tài)系統(tǒng)服務(wù)價值演變超前于人類活動強度，而關(guān)中地區(qū)的生態(tài)系統(tǒng)服務(wù)明顯滯后于人類活動；2030年，生態(tài)保護模式下二者的負(fù)相關(guān)作用最弱，而經(jīng)濟優(yōu)先模式下最強。

參考文獻(xiàn)：

［1］ 傅伯杰.加強基礎(chǔ)研究的系統(tǒng)性和綜合性，助力黃河流域生態(tài)保護與高質(zhì)量發(fā)展［J］.中國科學(xué)基金，2021，35（4）：503.

［2］ 習(xí)近平.在黃河流域生態(tài)保護和高質(zhì)量發(fā)展座談會上的講話［J］.中國水利，2019（20）：1-3.

［3］ 李征遠(yuǎn)，李勝鵬，曹銀貴，等.生態(tài)系統(tǒng)服務(wù)供給與需求：基礎(chǔ)內(nèi)涵與實踐應(yīng)用［J］.農(nóng)業(yè)資源與環(huán)境學(xué)報，2022，39（3）：456-466.

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

［5］ 謝高地，甄霖，魯春霞，等.一個基于專家知識的生態(tài)系統(tǒng)服務(wù)價值化方法［J］.自然資源學(xué)報，2008（5）：911-919.

［6］ 傅伯杰，于丹丹.生態(tài)系統(tǒng)服務(wù)權(quán)衡與集成方法［J］.資源科學(xué)，2016，38（1）：1-9.

［7］ AI Z，YAO Y，DANG X，et al.Study on the trade-off/synergy between ecosystem services and human well-being based on land use in Weinan City，China［J］.Environmental Science and Pollution Research，2024，31：13549-13562.

［8］ 彭赤彬，錢湛，姜恒，等.沅江流域產(chǎn)水服務(wù)功能的時空變化及驅(qū)動力分析［J］.人民長江，2023，54（6）：95-102，125.

［9］ 馬恩樸，蔡建明，韓燕，等.人地系統(tǒng)遠(yuǎn)程耦合的研究進展與展望［J］.地理科學(xué)進展，2020，39（2）：310-326.

［10］HAN W，SU X，LU H，et al.Impacts of human activity intensity on ecosystem services for conservation in the Lhasa River Basin［J］.Ecosystem Health and Sustainability，2023，9：0088.

［11］王帥，傅伯杰，武旭同，等.黃土高原社會-生態(tài)系統(tǒng)變化及其可持續(xù)性［J］.資源科學(xué)，2020，42（1）：96-103.

［12］耿甜偉，陳海，劉迪，等.縣域尺度下生態(tài)系統(tǒng)服務(wù)供需匹配及生態(tài)建設(shè)：以陜西省為例［J］.地域研究與開發(fā)，2021，40（2）：140-144，150.

［13］LIU M，JIA Y，ZHAO J，et al.Revegetation projects significantly improved ecosystem service values in the agro-pastoral ecotone of northern China in recent 20 years［J］.Science of the Total Environment，2021，788：147756.

［14］朱琳娜，趙牡丹，李云緋，等.西安都市圈生態(tài)系統(tǒng)服務(wù)價值與人類活動強度的時空關(guān)系［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報，2024，40（3）：325-334.

［15］LI G，WU Z，CHENG G，et al.Research on the spatiotemporal dynamic relationship between human activity intensity and ecosystem service value in the Three Gorges Reservoir Area［J］.Sustainability，2023，15（21）：15322

［16］王孟超，潘洪義，高成，等.基于政策導(dǎo)向的岷江流域土地利用模擬及生態(tài)系統(tǒng)服務(wù)價值評估［J］.長江流域資源與環(huán)境，2024，33（2）：322-335.

［17］金一諾，黃銀洲，尤鳳，等.基于Markov-PLUS模型的石羊河流域土地利用及生態(tài)系統(tǒng)服務(wù)價值模擬［J］.干旱區(qū)資源與環(huán)境，2024，38（5）：130-139.

［18］冉端，李江風(fēng).中國人口集聚與生態(tài)集聚的空間協(xié)調(diào)度研究［J］.統(tǒng)計與決策，2020，36（18）：54-57.

［19］王嵐，曹巍，黃麟.中國重大生態(tài)工程近40年生態(tài)成效整合分析［J］.生態(tài)學(xué)報，2024（7）：2673-2687.

［20］WEN X，WANG J，HAN X，et al.Where are the trade-offs in multiple ecosystem services in the process of ecological restoration？ a case study on the vegetation restoration area in the Loess Plateau，northern Shaanxi［J］.Land，2024，13（1）：70.

［21］肖芳，魏文穎，周斌.黃河流域生態(tài)系統(tǒng)服務(wù)與城市化時空交互作用分析：以黃河寧夏段為例［J］.人民長江，2023，54（8）：93-100.

［22］孫澤興，李汶怡，劉嘉敏，等.陜西省生態(tài)恢復(fù)綜合效益評估［J］.生態(tài)學(xué)報，2022，42（7）：2718-2729.

［23］BRYAN B A，YE Y，ZHANG J，et al.Land-use change impacts on ecosystem services value：Incorporating the scarcity effects of supply and demand dynamics［J］.Ecosystem Services，2018，32：144-157.

（編輯：黃文晉）

Spatiotemporal evolution and prediction of coordination between ecosystem

service and human activityLIU Tao¹，ZHANG Haipeng²，YANG Menghao³

（1.Henan Institute of Geological Sciences Co.Ltd.，Zhengzhou 450007，China； 2.Key Research Institute of Yellow River Civilization and Sustainable Development & Yellow River Civilization by Provincial and Ministerial Co-construction of Collaborative Innovation Center，Henan University，Kaifeng 475001，China； 3.College of Surveying and Geo-Informatics，North China University of Water Resources and Electric Power，Zhengzhou 450046，China）

Abstract： Exploring and predicting the evolution characteristics of relationship between ecosystem service and human activity is of great significance for promoting construction for the ecological security pattern of harmonious coexistence between human and nature.Taking the Yellow River Basin in Shaanxi Province as an example，this paper analyzed the spatiotemporal evolution characteristics and multi-scenario prediction of ecosystem services and human activities，as well as their coordination，through value equivalence method，human activity intensity model，PLUS model，and various statistical analysis methods.The results showed that the overall ecosystem service value of the Shaanxi Yellow River Basin increased by 5.433 billion yuan from 1980 to 2020.Compared with the result in 2020，the ecosystem service value under the inertia development mode and economic priority mode would decrease by 781 million yuan and 2.011 billion yuan respectively in 2030.The ecosystem service value under the cultivated land protection mode and ecological protection mode would increase by 27 million yuan and 824 million yuan respectively in 2030.From 1980 to 2020，the human activity intensity showed a "N" shaped change pattern.Compared with 2020，the human activity intensity increased under the economic priority mode and inertia development mode，while weakened under the cultivated land protection mode and ecological protection mode in 2030.In the study area，the negative correlation between ecosystem service and human activity was gradually increasing，and such coordination was poor.In 2030，the negative correlation between ecosystem service and human activity would be the weakest in the ecological protection mode，while the strongest in the economic priority model.Overall，the negative impact of economic priority mode is the greatest in the future，and the Yellow River Basin in Shaanxi Province should give priority to the ecological protection mode.

Key words： ecosystem service；human activity；multi-scenario simulation；coordination；Yellow River Basin