關鍵驅動力作用下的新疆生態系統服務時空格局分析*

李婧昕，許爾琪，張紅旗※

(1.中國科學院地理科學與資源研究所陸地表層格局與模擬重點實驗室，北京 100101； 2.中國科學院大學，北京 100049)

0 引言

生態系統服務是指人類能夠從生態系統中獲得的惠益[1-2]。隨著社會經濟的發展，人們對于生態系統中自然資源的多元需求與日俱增。聯合國千年生態系統評估(The Millennium Ecosystem Assessment，MA)報告顯示，全球已有63%的生態系統服務急劇退化，其中人為因素是主要驅動力[3]。宏觀上，人類活動對資源的開發及自然環境的改造直接體現在承擔著不同的生態功能的土地覆被及土地利用類型的變化上，進而對生態系統服務的供給與調節能力產生影響[4-6]，因此合理的土地利用方式是保障生態系統服務穩定有序輸出、實現可持續發展的重要途徑。

然而土地資源的有限性、多功能性及人類需求的多元性決定了土地利用方式的競爭性[7]，并使得相應生態系統服務之間也存在權衡與協同關系[8-10]，因此在土地利用狀況基礎上實現區域內重要生態系統服務功能的量化是實現對當地生態系統服務定量分析的基礎。目前，基于年鑒等統計數據和相關計量方法實現生態系統服務實際供給或潛在供給的數值分析是一種簡潔高效的評估手段[11-14]，優點是可以聯合大量指標對指定區域的生態系統進行分析，缺陷是這種方法通常建立在以最小行政區劃為單元的統計數據上，時間和空間分辨率受限，生態系統服務功能分布的空間異質性表達急需進行深入挖掘[15]； 更廣泛的方法則是將不同土地利用類型服務價值當量與已有的土地利用圖關聯，直接利用當量與相應土地利用類型面積的乘積來表示生態系統服務價值[16-19]，這一方法存在的問題為實例研究是否真正實現當量指數的本地化，此外以貨幣形式表示的絕對價值當量人為主觀性較強。因此，客觀性、全面性、本地化是構建生態系統服務評估體系的重要目標。

在評估生態系統服務的基礎上，明確生態系統服務的重要驅動因素及時空變化特征，是評估當地生態系統服務及變化的前提。但驅動因素之間、生態系統服務功能之間錯綜復雜的關系，以及驅動因素與服務功能之間可能存在的一對多、多對一及多對多的關系極大增加了生態系統服務驅動力分析的難度[15]，加之目前生態系統服務形成及傳輸理論尚不健全[1-2]，故對生態系統服務功能與驅動力之間關系的探索基本處于一對一的分析階段，因此首先需根據研究區實際情況識別重要的生態系統服務功能和關鍵驅動因素。如馬赫等[20]根據對瀘沽湖解譯得到的5種一級地類劃分出5種生態系統，通過不同指標變化矩陣的形式分別分析每種生態系統變化與指標變化的關系，最終針對生態系統發生的變化提出相應的規劃建議； 嚴恩萍等[21]將修訂后的生態服務價值當量匹配到通過解譯得到的不同土地利用類型中，分析自然因素、人類活動和政策驅動3個驅動因素與生態系統服務價值空間變化的關聯度，進而尋找出關鍵驅動因素；包翠榮[22]針對濕地生態系統，利用濕地資源調查數據并采用層次分析法評估每個指標對生態系統的影響，從而尋找到主要驅動因素。關鍵驅動因素有助于理解當地生態系統服務變化的原因，在一定情況下，可以通過合理可控的人為活動促進生態系統正向演變及生態系統內物質循環速率，從而得到更穩定高效的生態系統服務輸出[12]。

新疆地處我國西北部干旱區腹地，降水稀少、蒸發強烈，具有典型的溫帶大陸性氣候特征。近些年在全球氣候變化的作用下，暖濕化趨勢顯著[23-28]，是對當地生態環境產生影響的重要外部因素，而不斷增長的人口對當地有限資源的需求和不合理開發則形成激烈的內部矛盾[29-32]，以土地鹽堿化、沙漠化為代表的生態問題嚴重制約當地社會經濟發展[33-35]。作為我國重要的生態屏障區及“一帶一路”戰略核心區，新疆的生態環境狀況關系到國家的生態安全和發展戰略實施。文章基于這一背景，選取2000—2015年新疆全境為研究區，依據當地生態環境和土地利用實際情況，明確生態系統服務核心功能并建立評估體系； 選取當地生態系統服務的關鍵驅動因子，并探究關鍵因子驅動下的生態系統服務變化的熱點與冷點，以期為當地生態環境保護提供科學參考。

1 研究區概況及數據來源

1.1 研究區概況

圖1 新疆地形

新疆維吾爾自治區位于我國西北部，地處歐亞大陸中心，面積166.31萬km2，約占中國國土面積的1/6。新疆有獨特的三山兩盆地形特征，北部的阿爾泰山、南部的昆侖山及橫亙于中部的天山將新疆分為北疆和南疆(圖1)。北疆為溫帶大陸性干旱半干旱氣候，年均氣溫約8.5℃[36]，由于受到西風帶影響，全年降水量150～200mm； 而南疆則為溫帶大陸性干旱氣候，年均氣溫約12.8℃[36]，由于地形相對封閉，暖濕氣流無法到達，全年降水量僅為25～100mm。同時，南北疆經濟發展也較不均衡，全疆經濟核心區集中在北疆，南疆受資源匱乏、地形阻隔等因素影響，經濟發展相對滯后[37-38]。

1.2 數據來源

該文選取了2000年、2010年、2015年基于landsat8解譯得到的土地利用數據， 3期數據的土地利用類型一致，包括耕地、林地、草地、水域、居民地和未利用土地6個一級類型及24個二級類型，分辨率30m； 另外參考相同年份的年均降水量數據及歸一化植被指數數據，分辨率1km。以上數據均來源于中國科學院資源環境科學數據中心(http：//www.resdc.cn)。數據進行空間分析前均重采樣為250m。

2 研究方法

2.1 生態系統服務供給潛力評估

該文基于遙感目視解譯的土地利用類型對應的10個關鍵生態系統服務功能構建生態系統服務供給潛力評估體系。為盡可能符合研究區生態環境狀況，采取從一級地類至二級地類逐層分等定級原則使用非貨幣形式評價生態系統服務供給潛力。

Ling等根據前人研究[16， 39-40]并結合新疆干旱區的環境特點，基于6個一級地類對生態系統服務價值當量本地化[15]。為避免人為因素對貨幣形式生態系統服務價值產生干擾，該文對服務價值的相對大小進行排序，排序結果為各生態系統服務供給潛力數值的整數部分，當不同地類產生的相同生態系統服務價值相差懸殊時，小量級的生態系統服務可忽略不計。

對二級地類的生態系統服務的評價則基于Burkhard[1-2]在CORINE土地覆被分類體系框架下構建的生態系統服務供給能力矩陣。CORINE體系共分為44個土地覆被類型(https：//www.eea.europa.eu/publications/COR0-landcover)，通過對照兩種分類體系對各自土地覆被類型的定義將相同地類進行關聯，能夠得到每個一級地類下的二級地類生態系統服務大小排序，并作為各生態系統服務影響力數值的小數部分，計算方法為：

(1)

式(1)中，ESVi表示第i種二級地類的生態系統服務；ESVi1表示二級地類i所在的一級地類的生態系統服務整數值，即Inti1；ESVi2表示二級地類i所表示的生態系統服務小數值；N代表某一級地類下生態系統服務不為0的二級地類數量，n表示地類i在這些二級地類中的倒序排序。當一級地類對應的生態系統服務為0而二級地類對應同一服務不為0時，以0.1為間隔在整數0的基礎上進行累積； 反之當一級地類不為0而二級地類為0時，生態系統服務為0，最終結果見表1。

表1 生態系統服務供給潛力矩陣

地類編碼氣候調節水資源調節土壤形成及保持廢物凈化生物多樣性營養物質循環食物供給原料能源供給水源供給文化服務總值ⅠⅡⅠⅡFⅠⅡFⅠⅡFⅠⅡFⅠⅡFⅠⅡFⅠⅡFⅠⅡFⅠⅡFⅠⅡF111424.67323.67414.5313.5313.5616.5616.5000000212.535.341214.3313.3314.513.513.516.516.5000012.534.66221646.8535.75646.8545.8636.75343.8434.6646.8000545.852.92216.215.2516.215.216.2513.214.216.20015.247.92336.625.536.635.626.533.624.436.60035.6512426.435.7526.425.436.7523.444.826.40025.450.7331535.75434.75535.75424.67434.75434.75535.75000000333.7539.923225.524.525.514.3324.524.525.5000023.537.833315.2514.2515.250014.2514.2515.25000013.2531.75441323.5636.6000616.5515.5525.67414.5010.1616.33646.845.54223.526.4000015.515.3314.50016.3336.638.164323.526.4000015.515.3314.510.116.3336.638.264433.7546.800000000000026.6726.423.624613.2516.2000015.515.3300000016.226.48551000200000000000000200010.1000424.674.775200000000000012.510.10014.336.93530000000000000010.100000.1661010.1000000000200000100000000212.22.36210.1000000000000000012.22.36310.10000000010.111.5000022.44.16420.220.210.10012.520.200000042.866510.110.110.100000000000032.62.9660010.10000000000000032.62.76710.110.110.1000000000000000.3總值69.775.6551.844.571.2562.466526.525.6693.9 注:表中Ⅰ代表1級地類; Ⅱ代表2級地類; F為最終生態系統服務供給潛力值

2.2 生態系統服務主要驅動力分析

人口增加導致的人為活動加劇和氣候暖濕化趨勢共同作用于新疆的生態環境，同時伴隨著土地覆被類型和生態系統服務之間、不同生態系統服務之間的權衡和協同。研究表明西北干旱區蒸散發量自1950年到2010年整體無顯著波動且略呈下降趨勢，因此干旱區變濕趨勢主要由降水變化產生[41-44]。故該文以年均降水量代表當地的氣候變化情況； 人類活動干擾指數代表人口數量及空間分布對土地的利用強度的影響[13]； 歸一化植被指數表示當地不同土地利用類型的植被覆蓋情況。將3個因素作為研究區內生態系統服務重要驅動力，分析三者的時空分布及變化對于生態系統服務總量和時空分異特征的影響。其中改進后的人為活動干擾指數計算方法為：

(2)

式(2)中，HDIk表示研究區任意柵格k的認為活動干擾指數；Pki為柵格k對應的土地利用類型i的人為活動干擾強度參數，參考貢璐等[45]在新疆干旱區內陸河流域的研究結果，該研究區不同土地利用類型的取值見表2；d為柵格k到最近居民點或各類建設用地的歐式距離；IDW(d)是對歐式距離進行反距離加權，歐氏距離及反距離加權過程通過Arcgis10.2的空間分析模塊實現；Norm(IDW(d))則是對加權結果進行歸一化。

表2 人為活動干擾強度參數P

ID土地利用類型干擾強度ID土地利用類型干擾強度ID土地利用類型干擾強度11水田0.5433低覆蓋度草地0.2453其它建設用地0.9612旱地0.5441河渠0.1261沙地0.0821有林地0.1142湖泊0.1262戈壁0.0822灌木林0.1143水庫坑塘0.1263鹽堿地0.3823疏林地0.1144永久性冰川雪地0.0864沼澤地0.1524其它林地0.5446灘地0.1565裸土地0.0831高覆蓋度草地0.2451城鎮用地0.9666裸巖石質地0.0832中覆蓋度草地0.2452農村居民點0.9667其它0.08

考慮到不同降水量、干擾強度和植被覆被情況對當地10種主要生態系統服務功能有不同程度的供給和調節作用，即不同的驅動強度，因此先對3種驅動力進行驅動強度分級，其數值范圍為0～1，共分為5個等級[15]，每個等級對應的驅動強度值見表3。

表3 3種驅動因子對不同生態系統服務的驅動強度值

驅動因素強度分級數值范圍生態系統服務功能氣候調節水資源調節土壤形成及保持廢物凈化生物多樣性營養物質循環食物供給原料能源供給水源供給文化服務NDVI10～0.10.050.050.10.10.050.050.050.050.050.0520.1～0.30.40.250.40.40.250.250.250.40.250.2530.3～0.50.60.50.60.60.50.50.50.60.50.540.5～0.70.80.750.80.80.750.750.750.80.750.755>0.71111111111HDI1<0.010.050.10.050.10.050.050.050.050.050.120.01～0.03111110.750.510.5130.03～0.080.90.750.80.750.750.60.750.750.750.540.08～0.130.750.50.750.80.6110.610.255>0.130.50.250.50.750.40.80.850.250.850.1降水1>0.8111111111120.4～0.80.80.750.90.60.850.90.90.90.90.7530.2～0.40.60.50.80.40.750.80.550.550.550.540.05～0.20.40.250.70.20.40.70.250.250.250.255<0.050.20.10.60.10.10.60.10.10.10.1

單一驅動因子對當地生態系統產生的驅動力、所有因子產生的總驅動力及各驅動因子的驅動強度計算方法為：

(3)

(4)

(5)

式(3)～(5)中，DFj為第j個驅動因子的驅動力(j=1， 2， 3)αi為第i種土地利用類型對應的驅動強度數值，DF為所有納入分析的驅動因子產生的生態系統服務驅動力，DFIj表示驅動因子j對生態系統服務功能的驅動強度或生態系統服務變化對驅動因子的響應程度。

2.3 生態系統服務變化及熱點分析

生態系統服務結構及強度變化是多因素作用下生態環境變化的直接體現，分析研究區內生態系統服務時空變化及原因對當地生態環境保護、實現生態系統服務供需平衡具有重要意義。基于前文所述方法，分別計算關鍵因素驅動下不同時期生態系統服務供給潛力的差異，并對其進行熱點分析。該文參考Cai等對生態系統服務熱點地區的提取方法[46]，判斷如下：

H>X+δ

(6)

C

(7)

式(6)、(7)中，H和C分別代表熱點和冷點地區的生態系統服務潛力值，X和δ分別代表研究區內生態系統服務潛力平均值和標準差，即服務潛力值高于全區平均值與一倍標準差之和的地區均被視為生態系統服務的熱點地區。進一步將生態系統服務的熱點區與非熱點區范圍視為二值化圖，可以通過對不同時期熱點區重心位置的提取直觀反映生態系統服務高值區總體分布遷移特征； 而通過分析不同時期多因素驅動下的遷移特征有助于分析生態系統服務權衡形成的原因[47]。

3 結果與分析

3.1 生態系統服務供給潛力空間分析

3.1.1 土地利用狀況分析

新疆地區2000年、2010年和2015年3個時期土地利用圖級統計數據如圖2和表4所示。從全疆來看，由于中國面積最大的兩大沙漠塔克拉瑪干沙漠和古爾班通古特沙漠分別位于南疆和北疆地區，使得全疆未利用地占比超過60%，而3個時期沙地占比分別為22.52%,22.49%和22.47%，成為全疆面積占比最高的二級地類。15年來以耕地和建設用地為代表的高人為活動干擾區面積增加顯著，而林地、草地數量略有減少。北疆和南疆地區土地覆被類型均以草地和未利用地為主，但北疆地區草地面積多于未利用地， 3個時期平均面積占比分別為45.02%和38.30%； 南疆地區則呈現完全相反的景觀， 3個時期草地平均面積占比為25.30%，而未利用地面積高達67.94%； 未利用地中近95%為難以利用的沙地、戈壁和裸巖石質地，與北疆地區相比南疆可供開發利用的土地資源極為有限。

表4 2000—2015年全疆一級土地利用類型占比

%

圖2 3個時期新疆土地利用狀況分布

3.1.2 關鍵因子驅動下的生態系統服務空間分布

3個關鍵驅動因子作用下的生態系統服務潛力空間分布如圖3所示。可以看出北疆地區覆蓋了全疆大部分生態系統服務高值區域，其中天山地區為全疆生態系統服務最為集中的區域，另一個范圍相對較小的高值區位于阿爾泰山邊緣； 南疆高服務潛力區主要位于塔里木盆地西側及昆侖山北緣，雖然在2010年有局部提升，但15年來這兩個重要區域的生態系統服務均有顯著下降。

圖3 關鍵因子驅動下的生態系統服務空間分布

3.2 生態系統服務變化熱點及遷移分析

對生態系統服務變化顯著區域的識別有助于探索生態系統服務的權衡與協同關系。不同時期3個關鍵驅動因子共同作用下生態系統服務變化如圖4所示。其中2000—2010年總體呈增加趨勢，平均增效為2.27/柵格，且覆蓋范圍較廣，北疆主要體現在天山北緣的人工綠洲區和阿爾泰山南緣的天然綠洲區，變化主要由于這一區域10年來植被覆蓋增加； 南疆則主要體現在塔里木盆地以北、以南及西側邊緣，是降水與植被覆蓋增加同時驅動的效果。然而前10年生態系統服務增加顯著的區域卻在2010—2015年出現較大幅度的下降(圖5)，通過圖層間的相關性分析可知，南疆地區年降水量與生態系統服務下降的相關性為0.79，而北疆地區兩者的相關性高達0.83，由此可見當地脆弱的生態環境對以降水為代表的氣候變化極強的敏感性； 經過兩個不同變化趨勢的階段后， 15年來生態系統服務表現為總體上微弱的下降趨勢，平均減小幅度為1.45/柵格。減小幅度最大的地區出現在阿爾泰山北緣，使得降水對當地生態系統服務的驅動強度大幅降低。2010—2015年NDVI變化出現零散的低值斑塊，意味著植被覆蓋也出現較為破碎地退化； 生態系統服務增幅最大的地區位于天山北緣，另外南疆地區在塔里木盆地西側邊緣整體也呈現較明顯的生態系統服務增加趨勢。

圖4 生態系統服務變化空間格局

圖5 生態系統服務變化熱點冷點空間格局

不同時期的生態系統服務變化重心及其遷移路徑更直觀地體現出熱點與冷點的總體變化特征。由圖6可知，分區后的熱點遷移總體表現出由西向東的趨勢，而冷點遷移則基本呈現相反方向的變化。遷移路徑長度代表服務遷移強度。北疆地區熱點、冷點遷移路徑長度遠小于南疆地區，這雖然與北疆疆域較小有關，但相對封閉的小區域中的生態環境變化更容易打破當地原有的生態系統平衡。北疆冷點遷移強度遠小于熱點，說明北疆生態系統在面對外界干擾時表現出整體性，且抗性較強，同時在有適宜的物質或能量輸入介入生態系統時，易產生更高效的輸出，故熱點能夠形成較遠距離的遷移。而南疆的熱點與冷點遷移則呈現方向幾乎相反、距離相差不大的情況，一定程度上可以說明這一地區受到外界有利和不利因素影響后的響應程度相似，對不利因素的抗性不強。

圖6 生態系統服務變化熱點與冷點遷移特征

3.3 生態系統服務主要驅動力分析

3.3.1 人為干擾驅動分析

根據公式(4)可計算得到每個時期不同驅動因子的總驅動力，結果如圖7所示。2000—2015年間，新疆總人口數從1 846萬人增至2 360萬人[48-49]。人口數量的增加伴隨著建設用地面積的擴張。北疆地區由于經濟基礎和自然資源稟賦相對南疆受到的限制較小，故建設用地增長比例相對較高，土地利用類型的變化主要體現在北疆北部地區，南疆則主要體現在塔里木盆地北緣。2000年，2010年及2015年3個時期北疆人為活動干擾強度平均值分別為0.209、0.213和0.217，而南疆地區分別為0.145、0.146和0.148。在人為干擾強度增大的背景下，由其驅動的生態系統服務數值在北疆和南疆也均呈下降趨勢，北疆3個時期分別為7.108、7.106和7.104，南疆3個時期分別為6.721、6.696、6.650。

3.3.2 氣候變化驅動分析

3個時期北疆地區降水驅動的服務平均值分別為16.02，17.24和13.01，與之形成強烈反差的南疆地區則僅為6.78，7.86和5.49，由此可見，相比于人為活動對生態環境產生的干擾，自然環境對南疆地區生態系統服務提升的短板效應更顯著； 而從時間序列變化上看，無論北疆或南疆地區2015年較另外兩個時期生態系統服務均相對較低，這與2015年降水量少有直接關系。由此也可知，生態系統脆弱、生態系統服務低下的干旱半干旱地區對氣候因素的不穩定性會有極敏感的響應。

3.3.3 植被覆蓋驅動分析

與降水相似，總體上植被覆蓋度與生態系統服務供給、調節功能呈正相關(表3)。由于北疆地區草地覆蓋占比約為45%，遠高于南疆地區，故植被覆蓋對生態系統服務的驅動作用主要體現在北疆除沙地之外的大部分地區； 而對南疆地區的驅動則主要體現在天山以南及塔里木盆地狹長的邊緣地區。北疆地區3個時期NDVI驅動下的生態系統服務值分別為14.30，16.74和15.25，與之相比南疆地區僅為3.79，4.58和3.93。

3.3.4 關鍵因子驅動強度分析

由圖8結果可知， 3個因子在南北疆驅動強度差異較大。在社會經濟發展較好、草被覆蓋度高的北疆地區，人為活動及植被覆蓋對生態系統服務的影響已相對穩定，而本身不確定性較大的年降水量對生態系統服務的影響波動較大； 在南疆地區，由于植被覆蓋度極低，故從數值上分析NDVI所代表的植被覆蓋能夠提供的生態系統服務最低； 3個時期降水對當地生態系統服務驅動強度與年降水量大小變化也呈現相同趨勢。

圖7 3個時期不同驅動因子驅動力空間分布

圖8 南北疆不同因子驅動強度對比

4 討論與結論

4.1 討論

該文結合生態系統服務價值當量與生態系統服務功能供給矩陣兩種方法的優勢，與兩級土地利用類型逐級匹配，構建出符合研究區實際情況的生態系統服務供給潛力矩陣，結果表明該體系能夠很好展現出該區域的生態系統服務時空分布特征。耦合人類活動干擾指數、降水和歸一化植被指數作為干旱區生態系統服務及變化的關鍵驅動因子并分析三者對于環境背景差異較大的北疆和南疆的驅動強度，結果表明不同環境下關鍵因子的驅動強度差異顯著，這一方面直觀展現出不同環境中3個因子的重要性程度，另一方面則說明該文構建的生態系統服務供給潛力評估和關鍵驅動力組合對當地生態環境具有相當的敏感程度。通過關鍵驅動力作用下的生態系統服務變化熱點與冷點分析，可以明確造成當地生態環境變化的首要因素，進而在人為可控范圍內通過政策調整等方式指導當地人為活動向利于生態環境正向演化的方向進行。

該文在熱點冷點識別時采用了根據均值和方差的和與差來計算閾值的方法，除意義簡潔明確外，這一方法在研究區域較大時相對于Gi*統計法、核密度估計、莫蘭指數法等[50]有極高的計算效率，且基本不存在斑塊破碎化、只存在統計意義而無實際意義的情況。但這一方法只體現出生態系統服務數值上的熱點與冷點，對空間聚集情況分析不足，Gi*等方法能夠很好地對空間聚集情況加以計算。在未來的研究中，可以考慮將該文方法與Gi*等方法融合，結合熱點與冷點區變遷與研究區內重點保護區域的空間重疊度，分析不同的環境保護措施的保護成效，從而為當地相關政策的建立或調整提供依據。

4.2 結論

該文通過對兩級土地利用類型逐層分等定級原則評估研究區生態系統服務供給潛力，并耦合人類活動干擾指數、年降水量及歸一化植被指數3個驅動因素，評估三者的驅動強度并分析在其作用下新疆生態系統服務時空分布集變化情況，結果表明如下。

(1)新疆生態系統服務供給潛力的高值主要集中在天山周邊地區，最大值出現在阿爾泰山北緣，此外塔里木盆地以西也有分布狹長的高值區。北疆地區由于植被覆蓋度較高、氣候相對濕潤，生態系統服務潛力遠高于沙漠廣布、植被稀少的南疆地區。

(2)2010—2015年， 5年來生態系統服務的減少大于2000—2010年10年來生態系統服務增加，導致整個研究期間生態系統服務平均值呈小幅度下降。前10年生態系統服務熱點分布區大部分轉變為成為后5年的冷點區，其中降水變化驅動下的阿爾泰山北緣生態系統服務降低幅度最大。北疆生態系統在外界干擾下能夠體現出更強的抗性和整體性，也更易形成高效的生態系統服務輸出，而南疆地區對有利和不利因素的響應程度相似。

(3)3個時期北疆和南疆人為干擾強度均呈增加趨勢，由此驅動的生態系統服務相應呈減小趨勢。北疆地區3個時期由降水驅動的生態系統服務遠高于南疆地區；降水驅動的生態系統服務高于人類活動，但其不確定性造成的波動也較大。南北疆土地覆被類型的巨大差異使得植被覆蓋驅動下的生態系統服務懸殊。相比于人類活動的干擾，北疆地區以降水和植被覆蓋為代表的自然環境對生態系統服務驅動強度更高，而南疆人類活動和降水對生態系統服務的驅動則交替處于主導地位。