鄭德鳳，萬(wàn)巨影，白麗娜，呂樂(lè)婷

(遼寧師范大學(xué)地理科學(xué)學(xué)院，遼寧 大連 116029)

生態(tài)系統(tǒng)服務(wù)是生態(tài)系統(tǒng)供給人類的可使用的自然資源[1]，主要指人類通過(guò)直接或間接手段得到的可以支持生命發(fā)展的產(chǎn)品與服務(wù)，一般分為供給、調(diào)節(jié)、文化和支持4類典型生態(tài)系統(tǒng)服務(wù)[2]，可采用生態(tài)系統(tǒng)服務(wù)價(jià)值(ecosystem services value，ESV)來(lái)表征。不同時(shí)空變化過(guò)程中，各生態(tài)系統(tǒng)服務(wù)存在權(quán)衡/協(xié)同關(guān)系。權(quán)衡關(guān)系指地理要素空間分布存在某一項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值的增減，對(duì)其他生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生影響；協(xié)同關(guān)系指2種或2種以上生態(tài)系統(tǒng)服務(wù)價(jià)值同時(shí)增減。權(quán)衡/協(xié)同關(guān)系與人們對(duì)生態(tài)系統(tǒng)的需求相關(guān)[3-4]。國(guó)內(nèi)外學(xué)者針對(duì)濕地[5]、流域[6]、生態(tài)保護(hù)區(qū)[7]和農(nóng)耕區(qū)[8]，采用情景模擬法[9]、空間分析法[10]及相關(guān)統(tǒng)計(jì)學(xué)方法[11]等，圍繞生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同的理論基礎(chǔ)[12]、驅(qū)動(dòng)機(jī)制[13]及尺度效應(yīng)[14]展開了大量研究。盡管權(quán)衡/協(xié)同的研究?jī)?nèi)容較為廣泛，但研究尺度存在一定局限性，以往研究以單一尺度為側(cè)重點(diǎn)，分別對(duì)宏觀或微觀尺度網(wǎng)格單元[15-16]展開研究，而缺乏多尺度耦合探究。展開多空間尺度生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系研究，有利于促進(jìn)權(quán)衡/協(xié)同關(guān)系的轉(zhuǎn)換，通過(guò)時(shí)間尺度探究時(shí)間序列ESV演變規(guī)律，并結(jié)合空間尺度下網(wǎng)格單元、鄉(xiāng)鎮(zhèn)、縣域等微觀尺度梳理ESV權(quán)衡/協(xié)同空間分布形式，可以促進(jìn)區(qū)域生態(tài)保護(hù)措施落實(shí)，為生態(tài)系統(tǒng)可持續(xù)發(fā)展提供理論支撐。

燕山—太行山區(qū)作為國(guó)家14個(gè)集中連片特困區(qū)之一，跨越河北、山西和內(nèi)蒙古3個(gè)省份，承擔(dān)了如京津冀沙源治理、太行山綠化等重點(diǎn)生態(tài)建設(shè)任務(wù)。已有研究表明，生態(tài)系統(tǒng)服務(wù)的尺度關(guān)聯(lián)特征使同一組生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系在不同尺度存在差異[14-15]。關(guān)于燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系時(shí)空特征及動(dòng)態(tài)變化的研究較為薄弱，已有研究多集中于土地利用程度[17-18]、貧困度[19]對(duì)生態(tài)系統(tǒng)服務(wù)的影響，對(duì)ESV權(quán)衡/協(xié)同關(guān)系研究集中于河北省[20]、京津冀[21]等省域或城市群等大尺度范圍。在研究?jī)?nèi)容上，現(xiàn)有研究多側(cè)重生態(tài)補(bǔ)償與生態(tài)支撐區(qū)劃分[17]，對(duì)于生態(tài)系統(tǒng)服務(wù)關(guān)系研究較少；在研究區(qū)域方面，目前仍缺乏不同尺度生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系方面的相關(guān)研究。為此，根據(jù)燕山—太行山區(qū)區(qū)域特征對(duì)當(dāng)量因子進(jìn)行修正，估算生態(tài)系統(tǒng)服務(wù)價(jià)值并探討其空間格局變化情況，采用相關(guān)分析揭示燕山—太行山區(qū)時(shí)間序列各項(xiàng)生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系，并通過(guò)雙變量自相關(guān)法對(duì)比縣域、鄉(xiāng)鎮(zhèn)和網(wǎng)格單元3種研究尺度的生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系及空間分布特征，研究結(jié)果可為“后扶貧時(shí)代”燕山—太行山區(qū)生態(tài)系統(tǒng)維護(hù)與可持續(xù)發(fā)展提供依據(jù)。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況與數(shù)據(jù)來(lái)源

燕山—太行山區(qū)占地面積為9.3萬(wàn)hm2，地理位置為38°～43° N，113°～119° E，包括33個(gè)縣域。研究區(qū)內(nèi)降水分布不均，年均降水量為300～580 mm，區(qū)域土地利用類型以林地、草地和耕地為主(圖1)。

圖1 研究區(qū)位置與土地利用類型

2000、2005、2010、2015和2019年土地利用數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http:∥www.resdc.cn/)，分辨率為30 m。農(nóng)作物播種總面積、農(nóng)作物產(chǎn)量和糧食生產(chǎn)凈利潤(rùn)從2000—2019年《河北統(tǒng)計(jì)年鑒》《內(nèi)蒙古統(tǒng)計(jì)年鑒》《全國(guó)農(nóng)產(chǎn)品生產(chǎn)資料匯編》獲得。

土地利用數(shù)據(jù)分類標(biāo)準(zhǔn)與生態(tài)系統(tǒng)價(jià)值估算模型中的分類存在一定差別，需要對(duì)土地利用類型進(jìn)行重分類。根據(jù)謝高地等[22]提出的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表，將土地利用分為耕地、林地、草地、濕地、水域和荒漠6類。其中，燕山—太行山區(qū)未利用土地涵蓋沙地、鹽堿地、裸地、裸巖石質(zhì)地和沼澤，沙地和鹽堿地生態(tài)系統(tǒng)服務(wù)價(jià)值與荒漠一致，沼澤地生態(tài)系統(tǒng)服務(wù)價(jià)值則與濕地較為相近；裸土地、裸巖石質(zhì)地、城鄉(xiāng)工礦居民用地與裸地生態(tài)系統(tǒng)服務(wù)價(jià)值較為一致[8]，故將沙地、鹽堿地、裸土地、裸巖石質(zhì)地和城鄉(xiāng)工礦居民用地都?xì)w為荒漠[23-24]，沼澤地歸為濕地。

1.2 生態(tài)系統(tǒng)服務(wù)價(jià)值估算模型

生態(tài)系統(tǒng)為維持其內(nèi)部穩(wěn)定狀態(tài)，必須要消耗自身的一部分價(jià)值，參考謝高地等[22]提出的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估方法，并結(jié)合單位面積標(biāo)準(zhǔn)當(dāng)量對(duì)研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行修正，計(jì)算公式為

(1)

(2)

式(1)～(2)中，VES,z為z時(shí)期生態(tài)系統(tǒng)服務(wù)價(jià)值，元；Dz為單位標(biāo)準(zhǔn)當(dāng)量因子的價(jià)值量，元·hm-2；sk為k類糧食的播種面積，hm2；pk為凈利潤(rùn)，元·kg-1；qk為單位面積產(chǎn)量，kg·hm-2；s為糧食播種總面積，hm2；eij為土地利用類型i生態(tài)系統(tǒng)服務(wù)類型j的價(jià)值；Si為土地利用類型i的面積，hm2。單位標(biāo)準(zhǔn)當(dāng)量存在地區(qū)差異性，該研究確定2015年燕山—太行山區(qū)單位標(biāo)準(zhǔn)當(dāng)量為3 068.65元·hm-2，由此計(jì)算燕山—太行山區(qū)單位面積ESV，結(jié)果見表1。

表1 燕山—太行山區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)

1.3 相關(guān)分析與雙變量空間自相關(guān)方法

相關(guān)分析法用于分析變量之間的線性關(guān)系[25]，利用相關(guān)系數(shù)可對(duì)變量間的相關(guān)性強(qiáng)弱進(jìn)行分析，若為正值，表明兩變量間存在同時(shí)增加或減少的協(xié)同關(guān)系；若為負(fù)值，則兩變量存在相互抑制作用。相關(guān)系數(shù)計(jì)算公式為

(3)

地理現(xiàn)象在空間相互作用和空間擴(kuò)散的影響下通常具有不同程度的空間相關(guān)性[26]，一般采用Moran′sI指數(shù)探討多變量之間的空間相關(guān)性。雙變量空間自相關(guān)結(jié)果表現(xiàn)為高-高集聚、低-低集聚空間分布特征，表示生態(tài)系統(tǒng)服務(wù)之間為協(xié)同作用，低-高集聚、高-低集聚則表示權(quán)衡作用[27]。該研究引用此方法進(jìn)行雙變量生態(tài)系統(tǒng)服務(wù)空間上的相關(guān)性研究，相關(guān)公式參照文獻(xiàn)[26]。

2 結(jié)果與分析

2.1 燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)變化

由表2可知，2000—2019年研究區(qū)ESV為增—減—增波動(dòng)態(tài)勢(shì)，不同服務(wù)類型變化各異。其中，2005—2010年調(diào)節(jié)服務(wù)增加1.59億元，供給、支持和文化服務(wù)3項(xiàng)表現(xiàn)為減少趨勢(shì)。由圖2可知，燕山—太行山區(qū)各類用地ESV排序由高到低為林地、草地、耕地、濕地、水域和荒漠。各類用地ESV變化趨勢(shì)不一致，林地可產(chǎn)生較高ESV，占比為54.48%。耕地和草地ESV在2000—2015年減少，2015—2019年增加，2000—2019年耕地和草地ESV分別增加14.08億和14.21億元。2005—2019年濕地服務(wù)價(jià)值逐年減少，共減少63.85億元。水域ESV變化較為顯著，2005—2010年增加2.70億元，2010—2019年減少0.89億元。

由表3可知，燕山—太行山區(qū)各縣域ESV較高地區(qū)以河北省為主，如圍場(chǎng)、豐寧和隆化等，較低縣域?yàn)橥己屯昕h。2000—2019年ESV明顯增加的縣域有圍場(chǎng)、隆化和尚義，明顯下降的有康保、沽源和張北，分別減少3.41億、4.66億和4.93億元。

表2 2000—2019年燕山—太行山區(qū)各類生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)

表3 2000—2019年燕山—太行山區(qū)縣域平均生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)

圖2 2000—2019年燕山—太行山區(qū)各類用地生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)

2.2 燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系

2.2.1燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)的相關(guān)性

研究區(qū)各類生態(tài)系統(tǒng)服務(wù)相關(guān)系數(shù)見表4。由表4可知，11類生態(tài)系統(tǒng)服務(wù)組成的121組相關(guān)系數(shù)值均通過(guò)0.01水平的顯著性檢驗(yàn)。供給服務(wù)中，食物生產(chǎn)與氣體調(diào)節(jié)、氣候調(diào)節(jié)、凈化環(huán)境、生物多樣性和美學(xué)景觀存在權(quán)衡關(guān)系；調(diào)節(jié)服務(wù)中，氣體調(diào)節(jié)與氣候調(diào)節(jié)、凈化環(huán)境、土壤保持、維持養(yǎng)分循環(huán)和生物多樣性存在協(xié)同關(guān)系(相關(guān)系數(shù)>0.9)；氣候調(diào)節(jié)與凈化環(huán)境和土壤保持為協(xié)同關(guān)系(相關(guān)系數(shù)>0.9)；支持服務(wù)中土壤保持與維持養(yǎng)分循環(huán)存在協(xié)同關(guān)系，生物多樣性與美學(xué)景觀存在協(xié)同關(guān)系(相關(guān)系數(shù)=0.98)?？傮w來(lái)看，協(xié)同關(guān)系為研究區(qū)生態(tài)系統(tǒng)服務(wù)的主導(dǎo)關(guān)系。

2.2.2燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)的多尺度空間集聚特征

為探究燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同空間分布特征，采用等間距方式將研究區(qū)劃分為3 908個(gè)5 km×5 km的網(wǎng)格單元，分別在網(wǎng)格單元、鄉(xiāng)鎮(zhèn)和縣域尺度計(jì)算雙變量空間自相關(guān)系數(shù)。由表5可知，4類生態(tài)系統(tǒng)服務(wù)空間均呈協(xié)同關(guān)系，各尺度雙變量空間自相關(guān)系數(shù)由高到低為縣域尺度、網(wǎng)格單元和鄉(xiāng)鎮(zhèn)尺度。

以2019年為例，繪制網(wǎng)格、鄉(xiāng)鎮(zhèn)和縣域3種尺度下4類典型生態(tài)系統(tǒng)服務(wù)LISA集聚圖，以探究生態(tài)系統(tǒng)服務(wù)在不同尺度的空間集聚特征(圖3～5)。由圖3可知，網(wǎng)格尺度生態(tài)系統(tǒng)服務(wù)關(guān)系存在空間差異，表現(xiàn)為協(xié)同關(guān)系的區(qū)域主要分布在研究區(qū)北側(cè)，且分布較為集中。表現(xiàn)為權(quán)衡關(guān)系的區(qū)域分布較為分散，其單元格數(shù)量不超過(guò)500個(gè)。但供給服務(wù)與文化服務(wù)表現(xiàn)為權(quán)衡關(guān)系的區(qū)域分布較為集中，單元格數(shù)量共552個(gè)。

由圖4可知，鄉(xiāng)鎮(zhèn)尺度由433個(gè)單元格組成，其雙變量空間自相關(guān)系數(shù)較低，4項(xiàng)服務(wù)表現(xiàn)為協(xié)同關(guān)系的區(qū)域集中在研究區(qū)北側(cè)。其中，調(diào)節(jié)-文化服務(wù)表現(xiàn)為協(xié)同關(guān)系的鄉(xiāng)鎮(zhèn)數(shù)量(62個(gè))最多。權(quán)衡關(guān)系中，供給服務(wù)與調(diào)節(jié)和支持服務(wù)表現(xiàn)明顯，單元格數(shù)量為39個(gè)，但沒(méi)有形成集聚效應(yīng)。

由圖5可知，縣域尺度雙變量空間自相關(guān)系數(shù)較高，4項(xiàng)生態(tài)系統(tǒng)服務(wù)都表現(xiàn)為協(xié)同關(guān)系。供給-調(diào)節(jié)服務(wù)表現(xiàn)為協(xié)同關(guān)系的區(qū)域分布在研究區(qū)北部豐寧、隆化和圍場(chǎng)，中部懷安、陽(yáng)高和陽(yáng)原，以及南部唐縣。

表4 燕山—太行山區(qū)各類生態(tài)系統(tǒng)服務(wù)的相關(guān)關(guān)系

表5 燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值雙變量自相關(guān)分析

3 討論

3.1 生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化分析

2000—2019年燕山—太行山區(qū)ESV總體保持增加—降低—增加趨勢(shì)，降低趨勢(shì)主要出現(xiàn)在2015年，各項(xiàng)服務(wù)價(jià)值均有所降低，降幅較大的為食物生產(chǎn)和水文調(diào)節(jié)，食物生產(chǎn)ESV降低主要是由于耕地面積變化所致。截至2015年，研究區(qū)耕地面積逐年減少，這與工業(yè)化發(fā)展、城鎮(zhèn)化推進(jìn)密切相關(guān)[27]。建設(shè)用地侵占耕地面積，而迫于社會(huì)發(fā)展的經(jīng)濟(jì)壓力，農(nóng)民對(duì)耕地需求量增加，出現(xiàn)水域、濕地轉(zhuǎn)為耕地的逆生態(tài)化過(guò)程，這為生態(tài)環(huán)境帶來(lái)巨大壓力[28]。水文調(diào)節(jié)ESV降低是由于提供調(diào)節(jié)服務(wù)的水域和濕地面積減少，林地面積逐年增加，植被蒸發(fā)量增加[29]，生態(tài)系統(tǒng)調(diào)節(jié)功能減弱。

伴隨燕山小流域治理、太行山綠化等生態(tài)恢復(fù)工程建設(shè)的實(shí)施[19]，研究區(qū)林地面積逐年增加，生態(tài)用地的質(zhì)量和數(shù)量得到提升[30]，為研究區(qū)生態(tài)恢復(fù)與環(huán)境改善奠定基礎(chǔ)。作為燕山—太行山綠化建設(shè)的責(zé)任主體，縣級(jí)及以上政府應(yīng)落實(shí)行政責(zé)任制，依據(jù)山區(qū)地形地貌特點(diǎn)，因地制宜采取相應(yīng)建設(shè)對(duì)策，在ESV高值區(qū)采取封山育林等政策維持現(xiàn)有林區(qū)狀況，在ESV較低的丘陵區(qū)采取人工造林、育林措施，形成分區(qū)治理、統(tǒng)籌協(xié)調(diào)的生態(tài)環(huán)境保護(hù)系統(tǒng)。燕山—太行山區(qū)連接河北、山西和內(nèi)蒙古3個(gè)省份，應(yīng)利用好區(qū)域優(yōu)勢(shì)，合理規(guī)劃生產(chǎn)、生活、生態(tài)的三生空間布局，堅(jiān)守生態(tài)與發(fā)展2條線并進(jìn)，提高人們的生態(tài)保護(hù)意識(shí)，促進(jìn)區(qū)域生態(tài)環(huán)境可持續(xù)發(fā)展。

圖例中括號(hào)內(nèi)數(shù)值為單元格數(shù)量。

圖例中括號(hào)內(nèi)數(shù)值為單元格數(shù)量。

圖例中括號(hào)內(nèi)數(shù)值為單元格數(shù)量。

3.2 生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系分析

在探討ESV權(quán)衡/協(xié)同關(guān)系的基礎(chǔ)上，增加不同研究尺度的空間相關(guān)分析步驟，采用時(shí)空分異特性尋求平衡貧困區(qū)經(jīng)濟(jì)發(fā)展與生態(tài)環(huán)境修復(fù)的標(biāo)準(zhǔn)[19]。相關(guān)分析可用于計(jì)算較長(zhǎng)時(shí)間間隔內(nèi)網(wǎng)格單元變量的相關(guān)程度，雙變量自相關(guān)分析可表征縣域、鄉(xiāng)鎮(zhèn)和網(wǎng)格單元3種不同尺度空間集聚性。相關(guān)分析中協(xié)同關(guān)系為主導(dǎo)關(guān)系，但供給服務(wù)與其他3項(xiàng)服務(wù)存在不同程度的權(quán)衡關(guān)系。食物生產(chǎn)與調(diào)節(jié)服務(wù)中氣體調(diào)節(jié)、氣候調(diào)節(jié)和凈化環(huán)境3項(xiàng)服務(wù)存在權(quán)衡關(guān)系，這與孫藝杰等[13]研究發(fā)現(xiàn)的調(diào)節(jié)因子中凈初級(jí)生產(chǎn)力與食物供給存在沖突關(guān)系表現(xiàn)一致。其次，食物生產(chǎn)與生物多樣性和美學(xué)景觀之間存在權(quán)衡關(guān)系，燕山—太行山區(qū)為食物供給提供來(lái)源的土地類型多為耕地，耕地存在美學(xué)景觀價(jià)值，但其價(jià)值小于林地、草地等用地類型。鞏杰等[27]對(duì)甘肅白龍江的研究表明，食物供給與生境質(zhì)量間存在協(xié)同關(guān)系。筆者的研究結(jié)果與其存在差別，其原因與權(quán)衡/協(xié)同關(guān)系在不同地區(qū)存在地域差異相關(guān)[6]。

生態(tài)系統(tǒng)服務(wù)雙變量權(quán)衡/協(xié)同關(guān)系空間集聚特征存在尺度差異，縣域尺度雙變量相關(guān)系數(shù)與空間集聚程度明顯高于網(wǎng)格單元和鄉(xiāng)鎮(zhèn)尺度，這與劉華超等[7]、喬斌等[15]得出的生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同關(guān)系對(duì)研究尺度存在依賴性的結(jié)論相符。對(duì)網(wǎng)格單元、鄉(xiāng)鎮(zhèn)和縣域3種尺度的研究結(jié)果進(jìn)行對(duì)比分析，可以增加生態(tài)系統(tǒng)服務(wù)評(píng)估結(jié)果的可信程度。結(jié)合《河北省人大常委會(huì)關(guān)于加強(qiáng)太行山燕山綠化建設(shè)的決定》提及的將市、縣、鄉(xiāng)鎮(zhèn)分級(jí)規(guī)劃，各級(jí)單位配合落實(shí)區(qū)域生態(tài)建設(shè)[31]，從側(cè)面反映開展縣域、鄉(xiāng)鎮(zhèn)尺度研究的必要性，這有利于研究區(qū)結(jié)合實(shí)際情況，因地制宜開展生態(tài)環(huán)境保護(hù)相關(guān)工作。

3.3 不足與展望

基于當(dāng)量因子法的ESV估算模型可以運(yùn)用貨幣價(jià)值表達(dá)ESV，方法便于操作，但該模型也存在一定局限性[20]，如忽略了同種土地利用類型因所處環(huán)境不同而導(dǎo)致的結(jié)果差異；其次，由于計(jì)算過(guò)程中未涉及城市用地生態(tài)系統(tǒng)服務(wù)價(jià)值[18]，忽略了城市綠化和生態(tài)公園提供的ESV。未來(lái)應(yīng)對(duì)生態(tài)系統(tǒng)服務(wù)內(nèi)部驅(qū)動(dòng)因素展開進(jìn)一步研究，考慮不同情景下生態(tài)系統(tǒng)服務(wù)相互影響的動(dòng)態(tài)變化情況。

4 結(jié)論

(1)2000—2019年燕山—太行山區(qū)ESV表現(xiàn)為增加—降低—增加的變化趨勢(shì)。受綠化政策影響，林地ESV逐年增加，需要重點(diǎn)關(guān)注研究區(qū)林地、耕地和草地的轉(zhuǎn)換關(guān)系。

(2)燕山—太行山區(qū)生態(tài)系統(tǒng)服務(wù)以協(xié)同關(guān)系為主，供給服務(wù)中食物生產(chǎn)與氣體調(diào)節(jié)、氣候調(diào)節(jié)、凈化環(huán)境、生物多樣性和美學(xué)景觀存在權(quán)衡關(guān)系。

(3)不同研究尺度生態(tài)系統(tǒng)服務(wù)存在空間差異，網(wǎng)格單元與鄉(xiāng)鎮(zhèn)尺度的生態(tài)系統(tǒng)服務(wù)權(quán)衡關(guān)系分布分散，協(xié)同關(guān)系集中分布在研究區(qū)北部。研究區(qū)縣域尺度生態(tài)系統(tǒng)服務(wù)表現(xiàn)為協(xié)同關(guān)系，主要分布在荒山造林、荒地造林等工程覆蓋的縣域，如豐寧、隆化、圍場(chǎng)、懷安、陽(yáng)原、陽(yáng)高和唐縣等。