景海超,劉穎慧,賀 佩,張家琦,董婧怡,汪 怡

地表過程與資源生態(tài)國家重點實驗室,北京師范大學(xué)地理科學(xué)學(xué)部, 北京 100875

隨著1997年Costanza等[1]發(fā)表在Nature上的一篇關(guān)于核算全球生態(tài)系統(tǒng)服務(wù)價值的文章以及2001年聯(lián)合國千年生態(tài)系統(tǒng)評估項目[2]的開展,生態(tài)系統(tǒng)服務(wù)得到了越來越多的關(guān)注。生態(tài)系統(tǒng)服務(wù)作為連接自然和人類社會的橋梁和紐帶,是指生態(tài)系統(tǒng)與生態(tài)過程所形成及所維持的人類賴以生存的自然環(huán)境條件與效用,提供了人類賴以生存和發(fā)展的資源與環(huán)境基礎(chǔ)[3—4]。對生態(tài)系統(tǒng)服務(wù)的量化評估可以為生態(tài)安全、土地規(guī)劃等部門提供科學(xué)定量的生態(tài)環(huán)境信息,將有利于開展針對生態(tài)環(huán)境的管理和決策,這也使得近年來通過模型模擬的方式評估一區(qū)域生態(tài)系統(tǒng)服務(wù)并分析它們的特征與相互關(guān)系成為了熱點研究領(lǐng)域[5—6]。如Li等[6]以潮白河沿岸的城市帶為研究區(qū),模擬了當(dāng)?shù)毓烫?、土壤保持、防風(fēng)固沙、產(chǎn)水以及生境質(zhì)量五項生態(tài)系統(tǒng)服務(wù),并分析了五項服務(wù)之間的相互關(guān)系。Wang等[5]分析了京津冀地區(qū)產(chǎn)水、固碳、土壤侵蝕等服務(wù)的多年變化及權(quán)衡協(xié)同關(guān)系。

然而不同地區(qū)有著不同的環(huán)境及社會本底特征,各項生態(tài)系統(tǒng)服務(wù)的供給量必定有所區(qū)別,由此造成了生態(tài)系統(tǒng)服務(wù)在空間分布上的異質(zhì)性[7—8],識別生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要驅(qū)動因素成為了一項亟待解決的問題。目前對生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性驅(qū)動因子的研究多是通過建立生態(tài)系統(tǒng)服務(wù)與影響因子的相關(guān)關(guān)系、回歸模型或是通過情景模擬的方式分析驅(qū)動因子在不同情境下對生態(tài)系統(tǒng)服務(wù)的影響,如饒恩明等[9]和張玲玲等[10]用建立相關(guān)的方式分析了土壤保持、產(chǎn)水、生境質(zhì)量等服務(wù)與人口密度、GDP、氣候因子等的關(guān)系；侯文娟等[11]通過地理加權(quán)回歸分析了產(chǎn)流量與NDVI、景觀破碎化指數(shù)等因子的相互關(guān)系；Bai等[12]和Clerici等[13]則通過設(shè)置不同的情景,揭示氣候變化和土地利用類型對美國肯塔基州以及安第斯山脈東部兩個小流域生態(tài)系統(tǒng)服務(wù)時空格局的影響。上述研究雖然都在空間上分析了生態(tài)系統(tǒng)服務(wù)空間格局與他們影響因素間的關(guān)系,但沒有闡明具體哪些影響因素對該生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的形成起到主要作用。地理探測器是用來分析數(shù)據(jù)空間異質(zhì)性的有力工具[14],可用于揭示數(shù)據(jù)空間異質(zhì)性背后的主要驅(qū)動因子,并定量分析單個驅(qū)動因子及因子交互對因變量的解釋力度。

探索生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性并開展其主要影響因子的識別對地區(qū)生態(tài)系統(tǒng)服務(wù)功能調(diào)控和優(yōu)化具有重要意義。那曲市位于青藏高原的中心,是一個典型的生態(tài)脆弱地區(qū)[15]。同時那曲也是我國“兩屏三帶”[16]生態(tài)安全戰(zhàn)略格局中青藏高原生態(tài)安全屏障的核心區(qū)域,在水源涵養(yǎng)、氣候調(diào)節(jié)、生物多樣性維持等方面都具有重要的生態(tài)意義。隨著2018年5月那曲撤地設(shè)市,那曲市各層級、各項社會經(jīng)濟(jì)職能也逐漸的步入正軌,在未來經(jīng)濟(jì)建設(shè)過程中思考如何協(xié)調(diào)生態(tài)保護(hù),做到生態(tài)保護(hù)與民生改善同步提升意義重大。因此,本文以那曲市為例,詳細(xì)介紹了如何運用模型評估其多年的生態(tài)系統(tǒng)服務(wù)并分析生態(tài)系統(tǒng)服務(wù)的空間格局和相互關(guān)系,隨后引入地理探測器,選擇氣候、植被指數(shù)、土壤、地形、人類活動及土地利用等影響因子探究各生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要影響因素,試圖找出影響各生態(tài)系統(tǒng)服務(wù)空間分布的主要影響因子,研究結(jié)果可以為那曲市未來的生態(tài)系統(tǒng)功能優(yōu)化、生態(tài)安全格局構(gòu)建及生態(tài)文明建設(shè)等提供科學(xué)的決策依據(jù)。

1 研究區(qū)概況

那曲市(83°55′E—5°5′E,29°55′N—36°30′N)位于青藏高原生態(tài)安全屏障的核心區(qū)域,長江、怒江、瀾滄江、拉薩河等大江大河皆發(fā)源于那曲或支流流經(jīng)那曲[17],水資源極其豐富。那曲市湖泊濕地遍布且孕育著西藏最大面積的天然草場,在涵養(yǎng)水土、調(diào)節(jié)氣候的同時也為大量高原動植物物種提供了必要的生存場所,那曲市提供的多項生態(tài)系統(tǒng)服務(wù)功能必將會為整個高原乃至全國的生態(tài)安全保障做出積極貢獻(xiàn)。那曲市平均海拔約4500 m,整體地勢由西北向東南傾斜,西北部地勢高亢,海拔高度平均在5500 m左右,中西部地勢開闊,屬于高原丘陵地形,東部屬河谷地帶,多高山峽谷且山勢陡峭,海拔高度在3500—4500 m之間。那曲市屬亞寒帶氣候區(qū),基本特點是氣溫低,空氣稀薄,大風(fēng)天氣多,太陽輻射強(qiáng),同時其東南部地區(qū)易受夏季風(fēng)影響,降水與溫度條件明顯好于中西部地區(qū),那曲市的東南部地區(qū)也是藏北僅有的農(nóng)作物產(chǎn)區(qū),人口相對稠密,同時擁有少量的森林資源和灌木草場。那曲市的土地利用類型包括耕地、林地、草地、水域、建設(shè)用地、未利用地、永久性冰川積雪七大類,其中覆蓋面積最大的草地約占那曲市總面積的52%,林地、耕地、建設(shè)用地的占比較低(圖1)。

圖1 那曲市在青藏高原的位置、地形及2018年土地利用圖

2 數(shù)據(jù)來源及研究方法

2.1 數(shù)據(jù)來源

本研究使用的數(shù)據(jù)主要包括：土地利用數(shù)據(jù)、氣象數(shù)據(jù)、土壤數(shù)據(jù)、高程數(shù)據(jù)、歸一化植被指數(shù)(NDVI)、行政區(qū)劃數(shù)據(jù)、道路數(shù)據(jù)、人口密度數(shù)據(jù)等。其中2000、2010、2018年3期30 m分辨率的土地利用數(shù)據(jù)、土壤類型數(shù)據(jù)以及研究區(qū)縣級行政區(qū)劃數(shù)據(jù)均來自于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn/)。氣象數(shù)據(jù)中的氣溫和降水使用了國家氣象局CN05.1格點化觀測數(shù)據(jù)集[18],日照時數(shù)數(shù)據(jù)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://data.cma.cn/)。土壤質(zhì)地、土壤有機(jī)質(zhì)含量、土壤容重等數(shù)據(jù)通過聯(lián)合國糧農(nóng)組織(FAO)世界土壤數(shù)據(jù)庫以及國家青藏高原科學(xué)數(shù)據(jù)中心(http://data.tpdc.ac.cn)獲得[19]。DEM獲取自地理空間數(shù)據(jù)云(www.gscloud.cn/),并通過ArcGIS 10.3軟件提取研究區(qū)的坡度。NDVI使用了美國地質(zhì)調(diào)查局(USGS)陸地過程分布式數(shù)據(jù)檔案中心(http://lpdaac.usgs.gov/)的MOD13Q1數(shù)據(jù)。研究區(qū)的道路數(shù)據(jù)來自于國家基礎(chǔ)地理信息中心1：25萬全國基礎(chǔ)地理數(shù)據(jù)庫(https://www.webmap.cn/)。人口密度數(shù)據(jù)來自于WorldPop的人口密度數(shù)據(jù)集(https://www.worldpop.org/)。香農(nóng)多樣性指數(shù)(SHDI)借助Fragstats 4.2軟件計算得到。上述數(shù)據(jù)在預(yù)處理完成后均采用阿爾伯斯等積圓錐投影(Krasovsky_1940_Albers)并重采樣為1 km分辨率用于后續(xù)分析。由于獲取到的各類數(shù)據(jù)的年份集中在2000—2018年,且那曲市自然環(huán)境與社會發(fā)展在此期間均表現(xiàn)出一定變化,故本文選擇2000年、2010年、2018年3個年份的數(shù)據(jù)進(jìn)行后續(xù)研究。

2.2 研究方法

2.2.1生態(tài)系統(tǒng)服務(wù)估算

(1)產(chǎn)水服務(wù)

本研究使用InVEST模型的產(chǎn)水量模塊進(jìn)行了那曲市產(chǎn)水服務(wù)的模擬運算。產(chǎn)水模塊運用研究區(qū)的降水量、蒸散量、土壤根系深度、植物可利用含水量、土地利用等數(shù)據(jù),基于Budyko水熱耦合平衡假設(shè)[20]計算研究區(qū)每個柵格單元的年總產(chǎn)水量。公式如下：

(1)

(2)

式中,AETi表示柵格單元i上的年蒸散量(mm)；Pi表示柵格i上的年降水量(mm)；Yi為柵格i上全年的產(chǎn)水深度(mm)；PETi為柵格i上的潛在蒸散量(mm)；wi為經(jīng)驗參數(shù)。

PETi=Kc(li)×ET0i

(3)

(4)

PAWCi=54.509-0.132SAND%—0.003(SAND%)2-0.055SILT%-0.006(SILT%)2-0.738CLAY%+0.007(CLAY%)2-2.688OM%+0.501(OM%)2

(5)

式中,Kc(li)表示柵格i上的蒸散系數(shù),數(shù)值的確定參考了模型手冊[21]及文獻(xiàn)[22]；ET0i為柵格i的參考作物蒸散,通過Modified-Hargreaves法[23]計算得到；root.depth為根系深度參考了模型手冊[21]及文獻(xiàn)[24—25],Z為經(jīng)驗常數(shù)；PAWCi為植物可利用含水量[26],SAND%、SILT%、CLAY%、OM%分別表示土壤砂粒、粉粒、黏粒、有機(jī)質(zhì)的含量。

(2)土壤保持服務(wù)

InVEST土壤保持模塊通過氣象數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)以及DEM等數(shù)據(jù)計算研究區(qū)各柵格單元的潛在土壤侵蝕量和實際土壤侵蝕量,兩者之差則為研究區(qū)的土壤保持量,本文將采用此方法模擬那曲市各柵格單元上的土壤保持量。公式為：

RKLSi=Ri×Ki×LSi

(6)

USLEi=Ri×Ki×LSi×Ci×Pi

(7)

SDi=RKLSi-USLEi

(8)

式中,RKLSi表示柵格i上裸地的土壤侵蝕量(t/hm2)；USLEi為柵格單元i在植被覆蓋管理因子和水土保持措施因子作用下的土壤侵蝕量(t/hm2)；SDi為柵格單元i上的實際土壤保持量(t/hm2)；Ri是降雨侵蝕力因子(MJ mm hm-2h-1),采用Wischmeier的月尺度公式[27]計算；Ki是土壤可蝕性因子(t hm2h hm-2MJ-1mm-1),采用EPIC模型[28]估算,并通過修正公式[29]進(jìn)行本地化矯正；LSi為坡長坡度因子,參考InVEST模型手冊與已有研究[9, 21],模型中由DEM數(shù)據(jù)計算得到；Ci為植被覆蓋管理因子,參考了蔡崇法等[30—31]的研究,通過NDVI計算植被覆蓋度最終由植被覆蓋度計算得到各柵格單元的Ci值；Pi為水土保持措施因子,參考書籍[32]與前人研究[31, 33]獲得。

(3)生境質(zhì)量

InVEST生境質(zhì)量模塊以土地利用數(shù)據(jù)和威脅源數(shù)據(jù)為基礎(chǔ),計算各柵格單元的生境質(zhì)量得分,得分在0—1之間,數(shù)值越大表示生境質(zhì)量越好[21]。本文結(jié)合研究區(qū)實際情況與相關(guān)研究[34—36],在選擇耕地、建設(shè)用地、未利用地、國道、鐵路為威脅源數(shù)據(jù)的同時構(gòu)建了生境質(zhì)量參數(shù)表,最終運行InVEST生境質(zhì)量模塊得到了那曲市各柵格單元的生境質(zhì)量得分。生境質(zhì)量計算公式如下：

(9)

(10)

式中,Qij為柵格i的生境質(zhì)量；Hj為土地利用類型j的生境適宜度；Dij為柵格i的棲息地退化程度；k為半飽和常數(shù)；Z為默認(rèn)參數(shù)2.5；R為威脅因子的數(shù)量；Pr是r威脅因子的總柵格數(shù)；Wr是威脅因子r的權(quán)重；rp是威脅強(qiáng)度；irip表示rp對i的威脅程度；βi表示i的可達(dá)性；Sjr土地利用類型j對r威脅因子的敏感性。

(4)固碳服務(wù)

本文運用CASA模型[37]模擬那曲市各柵格單元上植被對碳的吸收,即固碳服務(wù)。CASA模型是一個典型的光能利用率模型,用于模擬植被的凈初級生產(chǎn)力(NPP),該模型模擬植被NPP所需的參數(shù)相對較少,同時部分參量可直接通過遙感手段獲得,使之成為目前應(yīng)用廣泛且不斷發(fā)展的NPP模型之一[38]。運用CASA模型進(jìn)行NPP計算所需的數(shù)據(jù)包括土地利用數(shù)據(jù)、氣象數(shù)據(jù)、NDVI數(shù)據(jù)等,計算公式如下[39]：

NPPi=0.5×Rs×FPAR×f1×f2×W×emax

(11)

式中,NPPi為柵格單元i處的植被凈初級生產(chǎn)力(gC/m2)；Rs表示地表太陽輻射(MJ/m2),參照童成立等[40]的研究,基于日照時數(shù)數(shù)據(jù)算得；0.5為經(jīng)驗常數(shù),與Rs相乘近似表示植物可利用的地表太陽輻射；FPAR為植物轉(zhuǎn)化光合有效輻射的比例；f1和f2為低溫和高溫對植物光能利用率的影響；W為水分限制因子；emax植被的最大光能利用率,數(shù)值參考朱文泉等的研究成果[39]。

2.2.2生態(tài)系統(tǒng)服務(wù)綜合熱點區(qū)域

對研究區(qū)生態(tài)系統(tǒng)服務(wù)綜合熱點區(qū)域的識別可以為當(dāng)?shù)厣鷳B(tài)保護(hù)措施的制定與決策提供定性的依據(jù)。本研究在柵格尺度上將各項生態(tài)系統(tǒng)服務(wù)供給量超過各自均值的區(qū)域視為該項生態(tài)系統(tǒng)服務(wù)供給的熱點區(qū)域,通過疊加研究區(qū)各對應(yīng)年份的產(chǎn)水量、土壤保持、生境質(zhì)量、固碳4項服務(wù)的熱點區(qū)域,最終得到研究區(qū)各對應(yīng)年份的生態(tài)系統(tǒng)服務(wù)綜合熱點區(qū)域的分布。根據(jù)區(qū)域上單項生態(tài)系統(tǒng)服務(wù)熱點的數(shù)量,本研究將熱點區(qū)域劃分為極重要區(qū)域、高度重要區(qū)域、中度重要區(qū)域、一般重要區(qū)域、非重要區(qū)域五類(表1)。

表1 那曲市生態(tài)系統(tǒng)服務(wù)重要區(qū)域劃分

2.2.3生態(tài)系統(tǒng)服務(wù)間的相互關(guān)系

本研究基于各生態(tài)系統(tǒng)服務(wù)的模擬值進(jìn)行了生態(tài)系統(tǒng)服務(wù)間逐柵格點的相關(guān)分析并檢驗其顯著性水平,探究那曲市各項生態(tài)系統(tǒng)服務(wù)間的相互關(guān)系。依據(jù)spearman相關(guān)分析,得到4種生態(tài)系統(tǒng)服務(wù)之間的相關(guān)性。公式為：

(12)

相關(guān)性結(jié)果為正值時,表明兩種生態(tài)系統(tǒng)服務(wù)之間具有協(xié)同效應(yīng),即兩種生態(tài)系統(tǒng)服務(wù)在同一時間段具有同樣的上升或降低趨勢,一種服務(wù)的增加會對另一種服務(wù)產(chǎn)生一定的促進(jìn)和增幅作用,結(jié)果為負(fù)值時,表明兩種生態(tài)系統(tǒng)服務(wù)具有權(quán)衡關(guān)系,即某生態(tài)系統(tǒng)服務(wù)的增加會造成另一種生態(tài)服務(wù)供給量的下降。

2.2.4生態(tài)系統(tǒng)服務(wù)空間分異的驅(qū)動因子識別

地理現(xiàn)象具有空間異質(zhì)性,地理探測器是揭示生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性背后主要驅(qū)動因子的有力工具,它的核心假設(shè)為：若某自變量會對因變量造成一定影響,則它們在空間上的分布會存在一致性[42]。地理探測器可量化各驅(qū)動因子對地理數(shù)據(jù)空間異質(zhì)性的解釋度,結(jié)果用q統(tǒng)計量表示,q值的含義為自變量解釋了q×100%的因變量。為了探究各單項生態(tài)系統(tǒng)服務(wù)受哪些影響因素的驅(qū)動,本文以研究區(qū)各單項生態(tài)系統(tǒng)服務(wù)為因變量,以影響那曲市生態(tài)系統(tǒng)服務(wù)的各因素為自變量,選用地理探測器中的“因子探測器”和“交互探測器”分析各驅(qū)動因子及因子組合對研究區(qū)生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的解釋程度。因子探測公式如下[42]：

(13)

(14)

SST=Nσ2

(15)

在生態(tài)系統(tǒng)服務(wù)驅(qū)動因子的選擇方面,以那曲市生態(tài)環(huán)境與社會經(jīng)濟(jì)背景為特征,最終考慮了氣候、植被指數(shù)、土壤、地形、人類活動及土地利用六類共九個指標(biāo)。氣候要素中選擇了年降水、年均溫；植被指數(shù)要素中選擇了歸一化植被指數(shù)(NDVI)；土壤要素中選擇土壤類型為主要驅(qū)動因子；地形要素中選取了高程和坡度因子；人類活動要素則主要考慮了人口密度；土地利用方面的驅(qū)動因子選擇了土地利用類型及表示土地利用類型破碎化的香農(nóng)多樣性指數(shù)(SHDI)。

交互探測用于評估不同驅(qū)動因子組合在一起對因變量的影響程度,兩因子的交互作用共有5種類型[42],若兩因子交互作用q值小于任意單因子q值,為非線性減弱；兩因子交互作用q值在兩個單因子q值之間,為單因子非線性減弱；兩因子交互作用q值大于任意單因子q值,為雙因子增強(qiáng)；兩因子交互作用q值等于兩個單因子q值之和,為獨立關(guān)系；兩因子交互作用q值大于兩個單因子q值之和,為非線性增強(qiáng)。

3 結(jié)果與分析

3.1 那曲市生態(tài)系統(tǒng)服務(wù)隨時間變化趨勢

2000—2018年那曲市產(chǎn)水量(圖2)總體呈明顯的下降趨勢,其平均產(chǎn)水深度由174.45 mm下降至113.19 mm,下降了35.1%。各縣區(qū)多年間的變化趨勢分為兩種類型,中部與東部縣區(qū)多表現(xiàn)為先降后增,西北部縣區(qū)多為逐年下降。那曲市各縣區(qū)2018年的平均產(chǎn)水深度普遍小于2000年的水平,尼瑪縣、雙湖縣、申扎縣的平均產(chǎn)水深度分別降低了76.9%、72.5%和68.5%。

那曲市土壤保持服務(wù)(圖2)在2000—2018年總體上表現(xiàn)為先降低后升高的趨勢,單位面積的土壤保持量由2000年的42.34 t/hm2下降至2010年的36.19 t/hm2后,隨后升高至2018年的40.32 t/hm2,2018年單位面積土壤保持量相較于2000年降低了4.8%,總體變化較小。各縣區(qū)多年間的變化趨勢與產(chǎn)水服務(wù)類似。2018年巴青縣、色尼區(qū)、聶榮縣、索縣單位面積的土壤保持量高于2000年,其余縣區(qū)的土壤保持量最終均低于2000年。

圖2 那曲市及各縣區(qū)2000、2010、2018年生態(tài)系統(tǒng)服務(wù)年均值

那曲市固碳服務(wù)(圖2)在2000—2018年總體呈上升趨勢,單位面積固碳量由2000年的219.07 gC/m2增長至2018年的231.40 gC/m2,增長了5.6%。那曲市各縣區(qū)2000—2018年的單位面積固碳量除申扎縣和班戈縣略微下降以外其余縣區(qū)均表現(xiàn)為增長,嘉黎縣、雙湖縣、安多縣分別增長了11.0%、10.7%和9.8%。

2000—2018年那曲市生境質(zhì)量(圖2)趨于穩(wěn)定,2000年、2010年和2018年那曲市生境質(zhì)量平均得分分別為0.364117、0.364131和0.365114,2018年得分相較于2000年提高了0.3%。那曲市各縣區(qū)的生境質(zhì)量得分在2000年至2018年間變化很小。

3.2 那曲市生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性

3.2.1單項生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性

那曲市各項生態(tài)系統(tǒng)服務(wù)的空間分布格局如圖3所示。各單項生態(tài)系統(tǒng)服務(wù)的分布都呈現(xiàn)出明顯的空間分異特征,同時各項服務(wù)的分布趨勢又存在一致性,即東南部地區(qū)的各單項生態(tài)系統(tǒng)服務(wù)普遍高于西北部地區(qū)。

圖3 2000—2018年那曲市生態(tài)系統(tǒng)服務(wù)時空格局

那曲市全市多年平均產(chǎn)水深度達(dá)133.79 mm,產(chǎn)水服務(wù)的空間分異明顯。那曲市產(chǎn)水服務(wù)的分布呈現(xiàn)出東南高西北低的特點,與該地區(qū)降雨量的分布特征基本吻合。產(chǎn)水量高值區(qū)主要分布在那曲市東部水熱條件較好的嘉黎縣、索縣、巴青縣、比如縣,該地區(qū)多高山峽谷且夏季受印度洋季風(fēng)的影響降水豐沛,四縣多年平均產(chǎn)水深度分別達(dá)到492.63 mm、486.11 mm、475.71 mm和470.66 mm,產(chǎn)水服務(wù)顯著高于西部各縣。產(chǎn)水量低值區(qū)主要分布在廣大的西部和北部地區(qū),該區(qū)地勢高亢,常年受西風(fēng)帶影響,降水稀少。尼瑪縣、雙湖縣多年平均產(chǎn)水深度僅為38.76 mm和45.23 mm。

那曲市單位面積土壤保持量的多年均值為39.62 t/hm2,土壤保持服務(wù)的空間分布格局與產(chǎn)水服務(wù)類似。土壤保持服務(wù)的高值區(qū)主要分布在東南部縣區(qū)的高山峽谷區(qū)域內(nèi),東南部區(qū)域豐沛的降水會使該區(qū)域受降雨侵蝕力的作用高于西北部,同時山區(qū)坡度變化大的區(qū)域會存蓄大量的流失土壤,東南部山地的山腰和山腳發(fā)育有那曲僅有的森林,保土能力強(qiáng)。索縣、嘉黎縣、比如縣、巴青縣多年單位面積土壤保持量分別為295.55 t/hm2、222.85 t/hm2、197.04 t/hm2、143.74 t/hm2,顯著高于全市平均水平。而那曲市西部及北部的雙湖縣、安多縣、尼瑪縣多年單位面積土壤保持量僅為6.02 t/hm2、10.49 t/hm2、20.28 t/hm2,明顯低于全市平均水平。

那曲市固碳服務(wù)的空間分布也具有明顯的分異性,全市單位面積固碳量多年均值為224.40 gC/m2。單位面積固碳量高的區(qū)域分布在那曲市中部和東部地區(qū),該地區(qū)氣候條件優(yōu)于西部與北部,同時該地區(qū)分布著大面積的高寒灌叢草甸與高寒草甸,西部與北部地區(qū)則主要為高寒草原。東部的比如縣、嘉黎縣、索縣的部分高山峽谷中還蘊(yùn)機(jī)著以冷杉、柏木、高山櫟等為主的森林。索縣、聶榮縣、色尼區(qū)、巴青縣單位面積固碳量多年均值分別達(dá)到489.16 gC/m2、478.12 gC/m2、436.44 gC/m2、411.75 gC/m2。低值區(qū)雙湖縣、尼瑪縣單位面積固碳量多年均值為135.98 gC/m2和174.39 gC/m2。

那曲市生境質(zhì)量得分的多年均值為0.36。生境質(zhì)量高值區(qū)主要分布在湖泊、河流及植被覆蓋狀況較好的區(qū)域,主要分布在那曲市的中部與東部地區(qū)。索縣、班戈縣、安多縣、比如縣的多年得分均值分別達(dá)到0.54、0.48、0.45、0.44。低值區(qū)即尼瑪縣、雙湖縣的多年得分均值僅為0.31和0.32。

3.2.2生態(tài)系統(tǒng)服務(wù)熱點區(qū)域的的空間異質(zhì)性

如圖4所示,那曲市生態(tài)系統(tǒng)服務(wù)綜合熱點區(qū)域的分布在2000—2018年間基本一致,各項生態(tài)系統(tǒng)服務(wù)高值區(qū)的重疊比例即高度重要區(qū)域與極重要區(qū)域的占比較低。其中,極重要區(qū)域在歷年來的占比保持在5.3%—5.4%之間,該區(qū)主要分布在那曲市中部與東部的高山峽谷與河谷地帶,植被與水熱條件明顯優(yōu)于其他地區(qū)。高度重要區(qū)域在歷年來的占比處于12.3%—13.8%之間,分布在那曲市中部與東部極重要區(qū)域的外圍。中度重要區(qū)域和一般重要區(qū)域的占比在多年間分別保持在16.2%—19.1%與32.2%—34.1%之間。非重要區(qū)域即沒有一項生態(tài)系統(tǒng)服務(wù)供給量超越其均值的區(qū)域占比處于28.6%—32.1%之間,該區(qū)域主要分布于那曲市西部與北部的尼瑪縣、雙湖縣及安多縣部分地區(qū),該區(qū)普遍海拔較高且干旱,溫度與降水條件明顯劣于那曲中部與東部區(qū)域,環(huán)境條件極為嚴(yán)苛,那曲北部區(qū)域也被劃歸在著名的“羌塘無人區(qū)”之中。

圖4 那曲市生態(tài)系統(tǒng)服務(wù)綜合熱點區(qū)域

3.2.3生態(tài)系統(tǒng)服務(wù)間的相互關(guān)系

在對那曲市四項生態(tài)系統(tǒng)服務(wù)進(jìn)行估算之后發(fā)現(xiàn)各生態(tài)系統(tǒng)服務(wù)的時空分布格局具有一定的一致性,各年份生態(tài)系統(tǒng)服務(wù)間的權(quán)衡與協(xié)同關(guān)系如表2所示,除固碳服務(wù)和生境質(zhì)量服務(wù)在2000年與2010年沒通過顯著性檢驗之外,其余年份各項生態(tài)系統(tǒng)服務(wù)之間的相關(guān)性均通過了0.01顯著性水平的檢驗。2000、2010、2018三年間各生態(tài)系統(tǒng)服務(wù)間的權(quán)衡與協(xié)同關(guān)系變化不大,穩(wěn)定在一定數(shù)值。其中產(chǎn)水服務(wù)、土壤保持服務(wù)與其余所有生態(tài)系統(tǒng)服務(wù)均表現(xiàn)為協(xié)同關(guān)系。產(chǎn)水量高的區(qū)域植被普遍長勢較好,保土能力強(qiáng),研究結(jié)果表明土壤保持服務(wù)和產(chǎn)水量服務(wù)之間有很強(qiáng)的協(xié)同關(guān)系,三個年份的相關(guān)系數(shù)保持在0.616—0.683之間。水和土壤是植被生長發(fā)育的基礎(chǔ),各年份中固碳服務(wù)與產(chǎn)水量服務(wù)、固碳服務(wù)與土壤保持服務(wù)、生境質(zhì)量與產(chǎn)水量服務(wù)間也表現(xiàn)出較強(qiáng)的協(xié)同關(guān)系。土壤保持與生境質(zhì)量服務(wù)之間雖然存在協(xié)同關(guān)系,但相關(guān)系數(shù)較低,土壤保持服務(wù)和生境質(zhì)量之間的相關(guān)系數(shù)多年間僅維持在0.11左右。

表2 那曲市生態(tài)系統(tǒng)服務(wù)間的相關(guān)性系數(shù)

3.3 那曲市生態(tài)系統(tǒng)服務(wù)空間分異的驅(qū)動因子識別

3.3.1生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的單因子歸因

那曲市各項生態(tài)系統(tǒng)服務(wù)均表現(xiàn)出了明顯的空間分異性,而這種分異性多是在自然因素和人類活動等的共同作用下形成的。本研究借助地理探測器的“因子探測”和 “交互探測”模塊得到了那曲市不同年份各驅(qū)動因子對不同生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的貢獻(xiàn)率(圖5)。由圖5可知,各驅(qū)動因子對那曲市生態(tài)系統(tǒng)服務(wù)的解釋力度在2000年、2010年、2018年間有較強(qiáng)的一致性,故本研究后續(xù)以2018年數(shù)據(jù)為例進(jìn)行那曲市生態(tài)系統(tǒng)服務(wù)影響因素的歸因分析。

圖5 2000—2018年各驅(qū)動因子對那曲市生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的影響程度

那曲市產(chǎn)水服務(wù)的空間分異受到多因子的共同作用,各因子按照q統(tǒng)計量的大小排序分別為降水(0.8785)>土壤類型(0.5910)>NDVI(0.5101)>溫度(0.2465)>坡度(0.1931)>土地利用類型(0.1828)>人口密度(0.1460)>高程(0.1212)>SHDI(0.0132)。因子探測結(jié)果表明,降水對那曲市產(chǎn)水服務(wù)空間分異的解釋力度最大,q值達(dá)到了0.8785,那曲市產(chǎn)水服務(wù)的空間分布也大致與降水的空間分布一直,表現(xiàn)出東南高、西北低的特點,年降水量處于749.89—880.93 mm之間的區(qū)域所帶來的產(chǎn)水量最高(圖6),平均產(chǎn)水量可達(dá)483.96 mm。同時,土壤類型、NDVI等自然因素對那曲市產(chǎn)水量服務(wù)的影響也極為顯著,那曲東部發(fā)育著大面積的草氈土和黑氈土而中西部和北部則多為寒鈣土,草氈土和黑氈土相較于寒鈣土有更強(qiáng)的肥力,且土層更厚,植被發(fā)育更好。草氈土和黑氈土覆蓋區(qū)的平均產(chǎn)水深度分別為314.92 mm、473.35 mm,平均產(chǎn)水深度最低的寒鈣土僅為33.92 mm。那曲市東部地區(qū)由于水熱狀況優(yōu)越,植被長勢好使得NDVI也明顯高于中西部與北部區(qū)域。NDVI處于0.73—0.99之間的區(qū)域產(chǎn)水量最大,平均產(chǎn)水深度為406.99 mm。那曲市人類活動因素對產(chǎn)水服務(wù)的影響有限,q值處在一個較低的水平。

圖6 那曲市各項生態(tài)系統(tǒng)服務(wù)隨驅(qū)動因子的變化趨勢

各影響因子對那曲市土壤保持服務(wù)空間分異的解釋力度如圖5所示,按照q統(tǒng)計量排序分別為坡度(0.5954)>土壤類型(0.2813)>降水(0.2634)>NDVI(0.1585)>高程(0.0748)>溫度(0.0720)>土地利用類型(0.0653)>人口密度(0.0395)>SHDI(0.0059)。坡度是對那曲市土壤保持服務(wù)影響最為顯著的因子,q統(tǒng)計量達(dá)到0.5954,土壤保持服務(wù)功能較強(qiáng)的區(qū)域主要集中在那曲市東南部的高山峽谷與河谷等坡度變化較大的區(qū)域(圖3),隨著坡度的增大,土壤保持量也呈現(xiàn)出一個增長的趨勢(圖6)。土壤類型、降水以及NDVI對土壤保持服務(wù)的空間分布也起到了較為重要的作用。不同土壤類型中,黑氈土覆蓋區(qū)域的土壤保持能力最強(qiáng),單位面積土壤保持量為356.22 t/hm2。土壤保持服務(wù)高值區(qū)與降水、NDVI的高值區(qū)也基本吻合,降水量處在749.89—880.93 mm之間、NDVI處于0.73—0.99之間的區(qū)域土壤保持服務(wù)也比較強(qiáng),單位面積土壤保持量分別為213.27 t/hm2、215.99 t/hm2。人類活動及土地利用因素對那曲市土壤保持服務(wù)空間分布的影響較小(圖5)。

各驅(qū)動因子對那曲市固碳服務(wù)空間異質(zhì)性的解釋力度比較高(圖5),按照q統(tǒng)計量排序分別為NDVI(0.6999)>土壤類型(0.4599)>降水(0.3992)>人口密度(0.2441)>溫度(0.2288)>高程(0.1586)>坡度(0.0558)>土地利用類型(0.0390)>SHDI(0.0234)。分析結(jié)果表明,NDVI對那曲市固碳服務(wù)空間分布的貢獻(xiàn)率最高,q統(tǒng)計量達(dá)到了0.6999。NDVI是指歸一化植被指數(shù),其高值區(qū)普遍為植被覆蓋水平較高的區(qū)域,故固碳量高,NDVI處于0.73—0.99之間的區(qū)域固碳量最大,平均固碳量達(dá)到了585.36 gC/m2。土壤類型因子方面,那曲市固碳服務(wù)高的區(qū)域發(fā)育有大面積的黑氈土、沼澤土、草氈土,其中黑氈土上發(fā)育著那曲市僅有的林地,沼澤土、草氈土上發(fā)育有大面積的優(yōu)質(zhì)牧草。那曲市黑氈土覆蓋區(qū)域單位面積的固碳量最高(圖6),可達(dá)到530.42 gC/m2。那曲市降水量高的地區(qū)固碳服務(wù)也比較強(qiáng),研究區(qū)降雨量處于642.67—749.89 mm之間的地區(qū)(圖6),其單位面積的固碳量最大,達(dá)到了477.01 gC/m2。

由圖5可知,各因子對那曲市生境質(zhì)量空間分異的解釋度按照q統(tǒng)計量排序分別為土地利用類型(0.3356)>NDVI(0.1757)>土壤類型(0.1350)>人口密度(0.1002)>SHDI(0.0374)>降水(0.0233)>高程(0.0168)>溫度(0.0083)>地形(0.0014)。土地利用類型是對研究區(qū)生境質(zhì)量空間分布貢獻(xiàn)率最大的影響因子,其次為NDVI和土壤類型。那曲市生境質(zhì)量的空間異質(zhì)性明顯,主要區(qū)域即湖泊、河流及林草地覆蓋區(qū)域的生境質(zhì)量要明顯優(yōu)于其他區(qū)域,生境質(zhì)量得分較高,而人口最密集、生產(chǎn)建設(shè)活動劇烈的建設(shè)用地區(qū)域生境質(zhì)量低。NDVI值越高的地區(qū)生境質(zhì)量的平均得分也越高,NDVI值處于0.73—0.99之間的區(qū)域,其生境質(zhì)量平均得分達(dá)到了0.54。除水域外,那曲市生境質(zhì)量最高的區(qū)域發(fā)育的土壤主要為黑氈土(圖6),黑氈土覆蓋區(qū)域的生境質(zhì)量平均得分為0.62。土地利用類型、NDVI、土壤類型等因子的解釋力度明顯高于高程和地形因子,地形與海拔對那曲市生境質(zhì)量的影響有限。

3.3.2生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性影響因子的交互作用探測

交互探測的結(jié)果用于表明任意兩因子的組合是否會增強(qiáng)或減弱它們各自對生態(tài)系統(tǒng)服務(wù)空間分布的解釋力度。各生態(tài)系統(tǒng)服務(wù)的因子交互探測結(jié)果顯示(表3),任意兩個因子的交互作用均大于單一因子對各生態(tài)系統(tǒng)服務(wù)空間分布的解釋力度,因子間兩兩交互的類型主要表現(xiàn)為“非線性增強(qiáng)”和“交互增強(qiáng)”,這樣的交互作用可以表明那曲市各項生態(tài)系統(tǒng)服務(wù)的空間格局是由多個因子共同作用的結(jié)果。

由表3可知,在產(chǎn)水服務(wù)中降水量與其他各因子的交互作用對那曲市產(chǎn)水服務(wù)空間格局的解釋力度最大,任意兩因子組合的q值均大于單一降水因子的q值(0.8785),其中降水∩土地利用類型(降水因子與土地利用類型因子的交互作用)的q值更是達(dá)到了0.9551,說明在同一降水水平下,土地利用類型的不同會對產(chǎn)水服務(wù)的空間格局造成很大的影響。NDVI∩其他因子以及土壤類型∩其他因子對產(chǎn)水服務(wù)空間格局的影響力也比較強(qiáng),q值處在0.5101—0.8984之間。

對土壤保持服務(wù)空間格局的交互探測結(jié)果表明(表3),地形因子中的坡度與其他因子的交互作用對土壤保持服務(wù)的空間格局影響最大,其中坡度∩NDVI、土壤類型等自然因子對土壤保持服務(wù)解釋力度達(dá)到了0.7以上,而坡度∩人口密度的q值僅比單一坡度的q值高了0.0082,也可說明人為因素對那曲市土壤保持服務(wù)空間格局的影響程度有限。降水、土壤類型與其他因子交互對土壤保持服務(wù)空間格局的解釋力度也比較大,q值處在0.2716—0.7248之間。

由表3可知,土地利用類型與其他因子交互作用對那曲市生境質(zhì)量空間格局的解釋力度最大,其次是NDVI∩其他因子以及土壤類型∩其他因子。土地利用類型∩NDVI以及土地利用類型∩土壤類型的q值分別為0.3895和0.3783,其次為土地利用類型∩人口密度和土地利用類型∩SHDI。

在固碳服務(wù)的因子交互作用中(表3),NDVI與其他因子的交互作用對那曲市固碳服務(wù)空間格局的影響程度最大,q值均在0.7以上,其次為土壤類型∩其他因子以及降水量∩其他因子。單因子的探測結(jié)果顯示土地利用類型的q值要遠(yuǎn)小于土壤類型和降水,然而交互作用中,NDVI∩土地利用類型的q值高于NDVI∩降水與NDVI∩土壤類型,達(dá)到了0.7915,表明在主導(dǎo)因素NDVI一致的情況下,土地利用類型的變化成為影響那曲市固碳服務(wù)空間格局的重要因素。

表3 驅(qū)動因子交互作用對那曲市生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的影響程度

4 討論與結(jié)論

4.1 討論

量化區(qū)域生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性并揭示異質(zhì)性形成的主要驅(qū)動因子可幫助決策者有針對性的施行生態(tài)優(yōu)化和調(diào)控措施。然而近些年的相關(guān)研究多是以回歸和相關(guān)為手段,沒能明確指出哪些驅(qū)動因子及因子組合會對生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性造成主要影響,而地理探測器在這一方面擁有獨特優(yōu)勢。本文嘗試在模型模擬生態(tài)系統(tǒng)服務(wù)時空格局的基礎(chǔ)上,用地理探測器探究那曲市各項生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要驅(qū)動因子,較好的揭示了那曲市各項生態(tài)系統(tǒng)服務(wù)的主控因子和因子組合,可為類似研究提供參考。本研究的模型模擬結(jié)果中那曲市各生態(tài)系統(tǒng)服務(wù)的值與臨近研究區(qū)已有的研究結(jié)果類似[36, 44—47]。如,潘韜等[45]的研究表明三江源2000年的平均產(chǎn)水深度為247.1 mm；楊蕾[44]在三江源的研究表明三江源多年平均土壤保持量為72.03 t/hm2；陳童堯等[46]在祁連山保護(hù)區(qū)的研究顯示祁連山區(qū)2015年的單位面積土壤保持量為212.63 t/hm2,該區(qū)域環(huán)境本底特征與那曲東南部縣區(qū)類似；鄒德富[47]的研究表明青藏高原2003—2010年的NPP均值為298 gC/m2。同時,上述研究中各項生態(tài)系統(tǒng)服務(wù)的空間分布趨勢也與本研究的模擬結(jié)果存在連續(xù)性,故模型模擬結(jié)果可信,能夠用于生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性歸因分析研究。

生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性是在眾多影響因子的共同作用下形成的。降水、土壤類型、NDVI這類自然因子對那曲市各項生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性都起到了重要作用,該結(jié)論與類似區(qū)域已有的研究結(jié)果一致[43, 48—49]。在這些主要影響因子的共同作用下,那曲市各項生態(tài)系統(tǒng)服務(wù)的空間格局存在相似性。除了上述不易受人為因素直接影響的自然因子外,土地利用類型可以在很大程度上解釋生境質(zhì)量的空間格局,從交互探測的結(jié)果來看,土地利用類型分別對那曲市固碳服務(wù)和產(chǎn)水服務(wù)空間異質(zhì)性起到重要影響,Hopping等[50]和Tang等[51]的研究也表明青藏高原土地利用及土地利用政策變化會對地區(qū)生態(tài)系統(tǒng)服務(wù)產(chǎn)生重大影響。在國家對口援藏政策[52]的持續(xù)支持下,近年來那曲市城鎮(zhèn)聚集效應(yīng)明顯提升,城鎮(zhèn)基礎(chǔ)設(shè)施建設(shè)也在逐步完善,伴隨而來的將是大量的其他用地類型向建設(shè)用地的轉(zhuǎn)變。2018年那曲市建設(shè)用地相較于2000年增長了近95%,一段時間內(nèi)那曲市城市、農(nóng)村居民點及工礦用地等建設(shè)用地面積還將持續(xù)增長,而這勢必會對當(dāng)?shù)氐纳鷳B(tài)系統(tǒng)服務(wù)產(chǎn)生影響。對此,建議那曲市應(yīng)當(dāng)：1)在未來城市建設(shè)與經(jīng)濟(jì)發(fā)展過程中繼續(xù)將生態(tài)保護(hù)擺在優(yōu)先位置,繼續(xù)推進(jìn)生態(tài)安全屏障建設(shè)；2)嚴(yán)格落實生態(tài)保護(hù)紅線,對土地利用進(jìn)行合理規(guī)劃,合理統(tǒng)籌生態(tài)用地、城鎮(zhèn)建設(shè)用地、工礦用地、農(nóng)業(yè)用地等功能空間；3)提升對全市各區(qū)域天然林、草地、濕地的保護(hù)力度,同時對已受建設(shè)及工礦活動影響的生態(tài)脆弱區(qū)采取劃區(qū)輪牧、退牧育草、大面積封禁管護(hù)等措施進(jìn)行恢復(fù)重建。

人類活動本身會對生態(tài)系統(tǒng)服務(wù)的空間分布帶來很大影響[53],但截至目前人口密度因子對各生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的解釋力度并不高,Li等[54]研究同樣表明在1990—2010年間青藏高原整體的生態(tài)環(huán)境受人類活動的影響較低。那曲市幅員遼闊但人口稀少,雖然近年來那曲市城鎮(zhèn)集聚效應(yīng)在逐步提升但城鎮(zhèn)化整體仍處在起步階段,建設(shè)用地占比較低。同時,隨著近年來陸續(xù)施行的限牧禁牧、退牧還草、生態(tài)補(bǔ)獎等[55—56]政策,放牧活動對草地的破壞也得到了一定的限制。這使得在那曲全域尺度上人類活動的規(guī)模和強(qiáng)度仍相對較小,對生態(tài)系統(tǒng)服務(wù)供給的影響尚處在較低水平。那曲大多數(shù)縣區(qū)剛剛實現(xiàn)脫貧摘帽,但大量居民尤其是牧民的生活水平仍然處在貧困線附近。今后,那曲的經(jīng)濟(jì)發(fā)展及城鎮(zhèn)化建設(shè)力度仍會增強(qiáng),生態(tài)保護(hù)和經(jīng)濟(jì)發(fā)展、民生改善之間的矛盾可能會更加突出。對此,建議那曲市應(yīng)當(dāng)：1)在帶領(lǐng)藏區(qū)群眾走向全面小康的同時,仍需加大力度落實生態(tài)扶貧及生態(tài)補(bǔ)償政策,在生態(tài)系統(tǒng)服務(wù)重要區(qū)域內(nèi)按照草蓄平衡原則,實施劃區(qū)輪牧、階段性禁牧和季節(jié)性休牧,在高海拔生態(tài)脆弱區(qū)積極推進(jìn)易地生態(tài)搬遷,在提高民生的同時將人類活動對生態(tài)環(huán)境的影響降至最低程度,努力做到生態(tài)環(huán)境與民生改善的協(xié)調(diào)發(fā)展；2)增強(qiáng)群眾的環(huán)保意識,對全市居民及外地游客廣泛宣傳生態(tài)環(huán)境相關(guān)的法律法規(guī)以及環(huán)保知識,呼吁群眾愛護(hù)自己賴以生存的生態(tài)環(huán)境；3)加強(qiáng)對各區(qū)域城鎮(zhèn)及農(nóng)村環(huán)境的綜合治理,實現(xiàn)城鎮(zhèn)及廣大農(nóng)村區(qū)域生活垃圾的收運處置,減少各類生活垃圾流通到環(huán)境中造成的環(huán)境污染；4)打造高質(zhì)量生態(tài)經(jīng)濟(jì)產(chǎn)業(yè),推動生態(tài)旅游業(yè)等特色產(chǎn)業(yè)的發(fā)展,在保護(hù)生態(tài)環(huán)境的同時獲取經(jīng)濟(jì)利益,提升附近居民生活水平。

本文所采用的研究方法中存在有一定的不確定性。InVEST模型和CASA模型憑借其操作簡單、參數(shù)輸入相對較少、采用空間數(shù)據(jù)作為模型的輸入與輸出并且模型的輸出結(jié)果可直接應(yīng)用于分析生態(tài)系統(tǒng)服務(wù)的時空分布格局的優(yōu)勢,獲得了越來越多科研人員的使用和掌握。本研究運用InVEST模型和CASA模型進(jìn)行生態(tài)系統(tǒng)服務(wù)的模擬評估,在各模塊的生物物理參數(shù)的計算與選擇過程中均結(jié)合研究區(qū)特點參考了已有的計算方法和相關(guān)文獻(xiàn)資料,但在參數(shù)的確定中仍然存在著不確定性,影響模擬精度。原因在于各參數(shù)(如產(chǎn)水模塊中的植被根系深度、土壤保持模塊中的水土保持措施因子、生境質(zhì)量模塊中各地類對威脅源的敏感性等)在基于研究區(qū)特點進(jìn)行本地化修訂時還需要針對研究區(qū)特點結(jié)合大量的野外觀測調(diào)研及實驗分析,從而驗證參數(shù)并確定參數(shù)應(yīng)用的合理性,進(jìn)而增強(qiáng)模型的模擬精度。本文使用地理探測器分析生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要驅(qū)動因子,探測結(jié)果可以定量化的表現(xiàn)各驅(qū)動因子和因子組合對生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的解釋力度,但此方法使用過程中在對連續(xù)的自變量因子進(jìn)行離散化處理時沒有明確的分類標(biāo)準(zhǔn),離散化方式的不同會對最終的探測結(jié)果產(chǎn)生一定影響,日后應(yīng)當(dāng)針對此方面加強(qiáng)分析,進(jìn)一步提升地理探測器方法的可靠性。在生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的評估過程中,本研究目前僅在全域尺度上分析了那曲市各生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要影響因素。然而在不同的空間尺度上同一生態(tài)系統(tǒng)服務(wù)的主要影響因子可能會存在差異,今后的研究中應(yīng)當(dāng)著力在更小的空間尺度上進(jìn)行生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的歸因分析,為不同區(qū)域制定區(qū)別化的生態(tài)保護(hù)政策提供依據(jù)。

4.2 結(jié)論

本文以模型模擬的方式評估了那曲市2000—2018年4項生態(tài)系統(tǒng)服務(wù),分析了它們的時空格局及相互關(guān)系,隨后借助地理探測器揭示了生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要影響因子及因子組合,得到如下結(jié)論：1)2000—2018年間,那曲市產(chǎn)水服務(wù)總體呈下降趨勢,土壤保持服務(wù)表現(xiàn)為先降低后增高,生境質(zhì)量服務(wù)趨于穩(wěn)定,固碳服務(wù)呈現(xiàn)上升趨勢。那曲市西北部各縣區(qū)生態(tài)系統(tǒng)服務(wù)年際變化趨勢相近,東南部各縣區(qū)生態(tài)系統(tǒng)服務(wù)的變化也存在一致性。2)2000—2018年間產(chǎn)水、土壤保持、固碳服務(wù)的高值區(qū)主要集中在那曲市東部與中部水熱與植被條件好的區(qū)域,尤其是高山峽谷與河谷地區(qū),生境質(zhì)量的高值區(qū)分布在湖泊河流以及植被覆蓋率高的地區(qū)。那曲市生態(tài)系統(tǒng)服務(wù)極重要區(qū)域與高度重要區(qū)域即各生態(tài)系統(tǒng)服務(wù)高值重疊區(qū)域則主要位于那曲市東南部高山峽谷與河谷區(qū)及其臨近地區(qū),為那曲市整體生態(tài)環(huán)境最優(yōu)的區(qū)域。3)2000—2018年間那曲市生態(tài)系統(tǒng)服務(wù)除固碳服務(wù)和生境質(zhì)量服務(wù)之間不存在明顯關(guān)系以外,其余組合均表現(xiàn)為相互促進(jìn)的協(xié)同關(guān)系。4)降水、土壤類型、NDVI等自然因子是對產(chǎn)水服務(wù)和固碳服務(wù)空間異質(zhì)性解釋力最強(qiáng)的驅(qū)動因子,坡度、土壤類型、降水等因子對土壤保持服務(wù)空間異質(zhì)性的解釋力度強(qiáng),生境質(zhì)量服務(wù)的空間異質(zhì)性主要受土地利用類型、NDVI、土壤類型等的影響,各生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的主要驅(qū)動因子存在相似性。此外,在交互作用中土地利用類型的重要性凸顯,土地利用類型∩NDVI對固碳服務(wù)和生境質(zhì)量有最大的解釋力度,土地利用∩降水可在最大程度上解釋產(chǎn)水服務(wù)的空間分異。截止目前代表人類活動的人口密度因子對各生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的解釋力度尚處在一個較低的水平。