張 帥，汪 洋，安沙舟

(新疆農(nóng)業(yè)大學(xué)草業(yè)與環(huán)境科學(xué)學(xué)院/新疆草地資源與生態(tài)重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830052)

0 引 言

【研究意義】生態(tài)系統(tǒng)服務(wù)是人類通過生態(tài)系統(tǒng)結(jié)構(gòu)、過程和功能直接或間接獲取到與福祉關(guān)系密切相關(guān)的產(chǎn)品和服務(wù)[1]，是開展生態(tài)環(huán)境穩(wěn)定可持續(xù)發(fā)展的基礎(chǔ)[2]。全球生態(tài)系統(tǒng)服務(wù)價(jià)值(Ecosystem services values，ESV)評(píng)估體系，使得生態(tài)系統(tǒng)服務(wù)研究逐漸成為了生態(tài)學(xué)中的熱點(diǎn)[3]。探尋ESV對(duì)土地利用格局演變的響應(yīng)機(jī)制，對(duì)生態(tài)系統(tǒng)可持續(xù)發(fā)展與修復(fù)具有重要意義?！厩叭搜芯窟M(jìn)展】Dawson等[4]認(rèn)為是否考慮社會(huì)復(fù)雜性對(duì)于生態(tài)系統(tǒng)服務(wù)的影響，是有關(guān)生態(tài)系統(tǒng)服務(wù)的研究成果能否為實(shí)際社會(huì)中的生態(tài)環(huán)境管理和調(diào)配提供科學(xué)依據(jù)的關(guān)鍵。Hassan等[5]提出不同研究區(qū)生態(tài)系統(tǒng)存在空間異質(zhì)性，在不同研究區(qū)的評(píng)估中使用同一套生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，其研究結(jié)果將會(huì)存在爭議。謝高地等[6]基于Costanza等的評(píng)估模型，依據(jù)中國生態(tài)系統(tǒng)的實(shí)際情況改進(jìn)了生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，建立了一套可廣泛應(yīng)用于評(píng)估中國城市[7]、島嶼[8]及流域[9]等生態(tài)系統(tǒng)的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)。目前，當(dāng)量因子法以更為直觀，易于使用且所需數(shù)據(jù)少的優(yōu)點(diǎn)，被廣泛應(yīng)用在區(qū)域和全球范圍內(nèi)評(píng)估ESV中[10]。土地利用/覆被變化(Land Use and Land Cover Change，LUCC)是反映人類活動(dòng)與陸地生態(tài)系統(tǒng)之間耦合關(guān)系最根本的表現(xiàn)，是生態(tài)系統(tǒng)服務(wù)功能轉(zhuǎn)變的重要驅(qū)動(dòng)因素[11]。土地利用變化的過程中引起了生態(tài)系統(tǒng)服務(wù)功能、結(jié)構(gòu)和過程的變更，使ESV隨之變化[12]。然而，近年來人類活動(dòng)的加劇和社會(huì)經(jīng)濟(jì)發(fā)展的不斷提速[13]，不合理的土地覆被類型變動(dòng)致使生態(tài)系統(tǒng)遭到嚴(yán)重破壞，ESV也隨之改變，使得生態(tài)系統(tǒng)服務(wù)功能退化[14]?！颈狙芯壳腥朦c(diǎn)】塔里木河流域位于我國西北干旱區(qū)生態(tài)脆弱帶[15]流域?qū)儆诮z綢之路經(jīng)濟(jì)帶的核心區(qū)[16]，是我國最大的內(nèi)陸河流域和重要的棉花生產(chǎn)基地[17]。近年來，塔里木河流域土地利用格局演變劇烈[18]、天然草地退化[19]、冰川萎縮[20]、河流斷流[21]、水資源調(diào)配不合理和荒漠化進(jìn)程加速等[22]，使得流域生態(tài)環(huán)境愈加脆弱[23]。關(guān)注塔里木河流域生態(tài)系統(tǒng)可持續(xù)發(fā)展與修復(fù)[24]?！緮M解決的關(guān)鍵問題】研究提出1990～2018年間塔里木河流域ESV表現(xiàn)為損失趨勢的假設(shè)，選取塔里木河流域1990、2000、2010和2018年4期的土地利用數(shù)據(jù)，基于ESV估算模型與探索性空間分析手段，定量研究塔里木河流域土地利用變化對(duì)ESV的影響及ESV的時(shí)空動(dòng)態(tài)變化趨勢，為干旱區(qū)土地資源可持續(xù)利用和區(qū)域生態(tài)環(huán)境的管理決策提供科學(xué)依據(jù)。

1 材料與方法

1.1 材 料

塔里木河流域位于我國西北內(nèi)陸干旱區(qū)腹地，塔克拉瑪干沙漠坐落于流域中部，地處73°10′～94°05′ E、34°55′～43°08′ N，平均海拔800～2 000 m，具有荒漠、綠洲、山地、高原等地貌景觀，是我國最大的內(nèi)陸河流域；屬于暖溫帶大陸性干旱氣候，由于其獨(dú)特地理位置和特殊的地形，致使該流域氣候干燥(平均氣溫10.7℃)、少雨(年均降水量17.4～42.8 mm)[25]。該流域是葉爾羌河、和田河、阿克蘇河、開都-孔雀河4條主要源流及塔里木河干流組成的114條河流的總稱。行政區(qū)包括巴音郭楞蒙古自治州、阿克蘇地區(qū)、和田地區(qū)、喀什地區(qū)、克孜勒蘇柯爾克孜自治州5個(gè)地州的94個(gè)縣(市)，總面積1.02×106km2，約占全疆總面積的61.4%。圖1

注：1.巴音郭楞蒙古自治州；2.阿克蘇地區(qū)；3.和田地區(qū)；4.喀什地區(qū)；5.克孜勒蘇柯爾克孜自治州；地圖來源：該圖中底圖遙感數(shù)據(jù)來源于Google Earth(www.earth.goole.com)Note:1.Bayingol Mongolian Autonomous Prefecture;2.Aksu Prefecture;3.Hotan Prefecture;4.Kashgar Prefecture;5.Kizilsu Kirghiz Autonomous Prefecture;Map source:the remote sensing data of the base map in the map is from Google Earth(www.earth.goole.com)圖1 研究區(qū)概況Fig.1 Sketch map of study area

從中國科學(xué)院環(huán)境科學(xué)數(shù)據(jù)中心網(wǎng)站(http://www.resdc.cn/)選取1990、2000、2010和2018年4期的塔里木河流域土地利用現(xiàn)狀數(shù)據(jù)，空間分辨率30 m，分類精度達(dá)到90%以上[26]。依據(jù)全國土地利用分類體系，并根據(jù)研究區(qū)實(shí)際情況，將研究區(qū)內(nèi)的土地利用類型分為耕地、林地、高覆蓋度草地、中覆蓋度草地、低覆蓋度草地、水域、建設(shè)用地和未利用地；數(shù)據(jù)來源于1990～2018年《新疆統(tǒng)計(jì)年鑒》[27]、《新疆生產(chǎn)建設(shè)兵團(tuán)統(tǒng)計(jì)年鑒》[28]及新疆維吾爾自治區(qū)統(tǒng)計(jì)局網(wǎng)站(http://www.xjtj.gov.cn/)等。

1.2 方 法

1.2.1 土地利用動(dòng)態(tài)度

該指標(biāo)可定量描述研究區(qū)不同土地利用類型在一定時(shí)間范圍內(nèi)的變化速率和幅度[27]，反映單一土地利用類型受到的影響。依據(jù)塔里木河流域1990、2000、2010和2018年土地類型數(shù)據(jù)，分析塔里木河流域28 a的土地利用/覆被變化。

(1)

式中，K為j到i的時(shí)間段內(nèi)該土地利用類型的動(dòng)態(tài)度，Ui、Uj分別表示研究時(shí)段開始與結(jié)束時(shí)該土地利用類型的面積；t為研究時(shí)長。

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

研究參照謝高地等[28]提出的基于中國實(shí)際情況改進(jìn)的Costanza評(píng)價(jià)模型[3]，得出的中國生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量。將農(nóng)田與難利用地劃分為耕地和未利用地[29]，高、中、低覆蓋度草地統(tǒng)一按照一級(jí)土地利用類型草地計(jì)入，根據(jù)研究區(qū)實(shí)際情況進(jìn)行修正[30]，塔里木河流域生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值系數(shù)。根據(jù)服務(wù)價(jià)值系數(shù)和塔里木河流域各土地利用類型的面積，計(jì)算研究區(qū)的ESV，計(jì)算公式如下：

ESV=∑Aa×VCa.

ESVb=∑Aa×VCba.

(2)

式中，ESV為生態(tài)系統(tǒng)服務(wù)價(jià)值，ESVb為第b項(xiàng)的生態(tài)系統(tǒng)服務(wù)價(jià)值；Aa為研究區(qū)a類土地利用類型的面積，VCa表示a類土地利用類型的生態(tài)服務(wù)價(jià)值系數(shù)，VCba表示a類土地利用類型的第b項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值。表1

表1 塔里木河流域生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(元/hm2)Table 1 ESV coefficient of Tarim River Basin(yuan/hm2)

1.2.3 空間自相關(guān)性

基于Geoga模型中的Global Moran’s I和Univariate local Moran’s I探究ESV空間分布格局的集聚和異常[31]。公式為：

(3)

(4)

2 結(jié)果與分析

2.1 土地利用類型時(shí)空變化特征

研究表明，該流域未利用地占據(jù)研究區(qū)總面積的64%以上，是該研究區(qū)分布最廣、面積最大的土地利用類型，耕地是研究區(qū)面積增長最多的土地利用類型，增長了1.75%；低覆蓋度草地是研究區(qū)面積第二大的土地利用類型，占總面積的14%以上，研究時(shí)段內(nèi)表現(xiàn)為減少態(tài)勢，共減少了2.39%。林地、高覆蓋度草地、中覆蓋度草地、水域分別占總面積的1.25%、5.49%、6.12%及2.66%；建設(shè)用地面積最小，僅占總面積的0.29%。表2，圖2

表2 塔里木河流域1990～2018年土地利用類型結(jié)構(gòu)Table 2 Land use type in the Tarim River Basin from 1990 to 2018(%)

圖2 塔里木河流域1990～2018年土地利用現(xiàn)狀Fig.2 Status of land use in the Tarim River Basin from 1990 to 2018

研究表明，1990～2018年間，塔里木河流域各土地利用類型之間轉(zhuǎn)化頻繁，且變化特征存在顯著差異，主要以耕地、建設(shè)用地及未利用地的增加和林地、草地及水域的減少為主。耕地是各類土地利用類型中增長趨勢最顯著的，增長了226.91×104hm2，其中2010～2018年增長速率最快，動(dòng)態(tài)度為1.920，主要由中覆蓋度草地(56.37×104hm2)、低覆蓋度草地(73.84 ×104hm2)和未利用地(47.94×104hm2)轉(zhuǎn)化而來；建設(shè)用地1990～2018年間增長了25.32×104hm2，其中2010～2018年增長速率最快，動(dòng)態(tài)度為3.689，主要由耕地(12.91×104hm2)和未利用地(7.30×104hm2)轉(zhuǎn)化而來；未利用地1990年為6 634.59×104hm2，2018年為6 748.77×104hm2，增長速率最快的時(shí)間段為2010～2018年，動(dòng)態(tài)度為0.088，主要由中覆蓋度草地(928.85×104hm2)、低覆蓋度草地(928.85×104hm2)和水域(152.57×104hm2)轉(zhuǎn)化而來，由于近些來氣溫增長、水資源的不合理調(diào)配與人類活動(dòng)強(qiáng)度加重，使得荒漠化進(jìn)程加速；林地28年間共減少了5.27×104hm2，減少速率最快的時(shí)間段為2000～2010年間，動(dòng)態(tài)度為-0.45，主要轉(zhuǎn)化為高、中、低覆蓋度草地(23.24×104hm2、18.12×104hm2、14.94×104hm2)、未利用地(16.16×104hm2)及耕地(18.95×104hm2)；高覆蓋度草地在1990～2018年間增加了90.11×104hm2，中、低覆蓋度草地確呈現(xiàn)下降趨勢，分別減少了31.98×104hm2和243.24×104hm2，高覆蓋度草地在2010～2018年間增加最快，動(dòng)態(tài)度為0.957，中、低覆蓋度草地減少速度最快的時(shí)間段分別為2000～2010年、2010～2018年，動(dòng)態(tài)度分別為-0.157、-0.468，高覆蓋草地主要由中覆蓋度草地轉(zhuǎn)化而來，中覆蓋度草地主要轉(zhuǎn)化為了低覆蓋度草地，低覆蓋度草地主要轉(zhuǎn)化為未利用地(769.19×104hm2)和耕地(73.84×104hm2)。表3，表4

表3 塔里木河流域1990～2018年土地利用類型動(dòng)態(tài)度Table 3 Dynamic degree of land use types in the Tarim River Basin from 1990 to 2018(%)

表4 塔里木河流域1990～2018年土地利用轉(zhuǎn)移矩陣Table 4 Land use transfer matrix of the Tarim River Basin from 1990 to 2018(104 hm2)

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

2.2.1 生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)間維度的變化特征

研究表明，1990～2018年塔里木河流域ESV整體呈現(xiàn)出下降趨勢，減少了1 057.19×108元，約12.11%。1990、2000和2010年研究區(qū)ESV分別為8 728.84、8 778.50和8 754.76×108元，變化趨勢較為平緩；2018年研究區(qū)ESV為7 671.65×108元，呈現(xiàn)明顯下降趨勢。依據(jù)各地類的生態(tài)系統(tǒng)服務(wù)價(jià)值變化可以發(fā)現(xiàn)水域和草地的ESV下降，28 a分別減少了1 026.32×108和261.65×108元，而耕地和未利用地的ESV呈現(xiàn)上升趨勢，分別增長了236.14×108和9.42×108元。表5

表5 1990～2018年塔里木河流域各土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值變化Table 5 Changes in ESV of various land use types in the Tarim River Basin from 1990 to 2018(108 yuan)

塔里木河流域的氣體調(diào)節(jié)、氣候調(diào)節(jié)、水源涵養(yǎng)和廢物處理服務(wù)功能整體上均呈現(xiàn)出下降趨勢，其中氣體調(diào)節(jié)、水源涵養(yǎng)和廢物處理服務(wù)功能在1990～2000年出現(xiàn)了輕微的回升，水源涵養(yǎng)和廢物處理功能在2010～2018年呈現(xiàn)劇烈下降的趨勢，分別下降了493.5、425.29×108元。調(diào)節(jié)服務(wù)功能主要受到耕地、林地、草地和水域的影響，由于研究區(qū)草地面積占總面積的比重較其他3個(gè)地類大，研究區(qū)調(diào)節(jié)服務(wù)的變化與草地面積的變化保持較高的相關(guān)關(guān)系。水源涵養(yǎng)和廢物處理服務(wù)功能呈現(xiàn)先上升后下降的趨勢，與草地、水域面積的變化一致。

塔里木河流域土壤形成與保護(hù)和生物多樣性保護(hù)2項(xiàng)服務(wù)功能并稱為支持服務(wù)，這2項(xiàng)功能在1900～2018年間始終保持平緩下降趨勢。支持服務(wù)主要受耕地、林地和草地的影響，林地的波動(dòng)性下降、耕地的持續(xù)性擴(kuò)增與草地的持續(xù)性下降，但在面積仍然占比很高的三者綜合作用下致使該項(xiàng)服務(wù)保持平緩下降的趨勢。

食物生產(chǎn)與原材料生產(chǎn)2項(xiàng)服務(wù)功能在研究區(qū)中生態(tài)系統(tǒng)服務(wù)價(jià)值所占比重較小。其中，供給功能主要受耕地、林地和草地的影響，耕地面積的擴(kuò)增與草地面積下降，使得農(nóng)作物和畜產(chǎn)品的產(chǎn)量增加，食物生產(chǎn)功能在1990～2010年呈現(xiàn)出平緩上升的趨勢，在2010～2018年出現(xiàn)了較為劇烈的增長；原材料生產(chǎn)服務(wù)功能在1990～2000年間呈現(xiàn)輕微的回升趨勢，在2000～2018年則保持持續(xù)平緩下降趨勢，這一變化與耕地和草地的面積變化趨勢相似，供給服務(wù)的變化主要由于耕地、林地和草地3個(gè)地類綜合作用，整體上保持上升趨勢。文化服務(wù)功能在1990～2000年呈現(xiàn)上升趨勢，在2000～2018年呈現(xiàn)出下降趨勢，其中2010～2018年的下降趨勢較為劇烈，該項(xiàng)服務(wù)功能主要受到林地、草地和水域的影響，其變化趨勢與林地和水域面積的變化趨勢一致。

塔里木河流域1990～2018年生態(tài)系統(tǒng)服務(wù)價(jià)值整體上保持下降趨勢，與調(diào)節(jié)服務(wù)變化趨勢一致，2000年以前呈現(xiàn)出輕微上升趨勢，2000年以后保持下降趨勢，2010～2018年下降趨勢劇烈。表6

2.2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值空間維度的變化特征

借助ArcGIS 10.2軟件的漁網(wǎng)工具將研究區(qū)劃分為大小10 km×10 km[32]的10 641個(gè)網(wǎng)格，分別計(jì)算網(wǎng)格內(nèi)各地類面積，采用自然斷點(diǎn)法對(duì)研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值與生態(tài)系統(tǒng)服務(wù)價(jià)值的損益進(jìn)行空間化展示。

研究表明，除了塔里木河流域東北部地區(qū)以外，ESV存在中部地區(qū)低、四周高的空間分布特征，這與研究區(qū)獨(dú)特的地貌類型分布特征高度契合，即盆地中部分布著大面積的沙漠，屬于未利用地，ESV較低；四周的平原地區(qū)則環(huán)繞著沙漠分布著以耕地為主的綠洲地帶以及分布著大面積草地的山地地區(qū)，ESV較高；ESV的高值區(qū)主要分布在研究區(qū)的北部、西部和西南部地區(qū)的草地和水域。對(duì)比1990和2018年ESV空間分布變化，高值區(qū)呈現(xiàn)出縮減態(tài)勢，低值區(qū)的面積逐年擴(kuò)大，整體上來看，研究區(qū)ESV呈現(xiàn)出顯著下降趨勢。

在過去的28 a間，塔里木河流域ESV整體上減損較大，主要是由于人類活動(dòng)和城市經(jīng)濟(jì)發(fā)展進(jìn)程加速造成流域內(nèi)ESV較高的用地類型(如草地、水域)大量轉(zhuǎn)化為ESV較低的用地類型(如未利用地)。而在研究區(qū)的東南部地區(qū)出現(xiàn)板塊狀分布的增益區(qū)域，主要是由于水域和草地面積的增長；研究區(qū)東北部地區(qū)還存在一些以點(diǎn)狀分布的增益區(qū)，主要得益于生態(tài)修復(fù)后局部的草地、水域面積的增長；研究區(qū)西部地區(qū)的ESV增益區(qū)則主要得益于水資源開發(fā)使得部分區(qū)域草地覆蓋度增加。圖3

表6 1990～2018年塔里木河流域各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化Table 6 Changes in ESV in the Tarim River Basin from 1990 to 2018(108 yuan)

圖3 塔里木河流域ESV空間變化分布特征Fig.3 Spatial distribution characteristics of ESV in the Tarim River Basin

2.3 生態(tài)系統(tǒng)服務(wù)價(jià)值空間相關(guān)性

2.3.1 空間自相關(guān)性

研究表明，并通過Local Moran’s I指數(shù)分析出研究區(qū)ESV的空間集聚結(jié)構(gòu)。研究區(qū)的ESV分布存在非常顯著的空間自相關(guān)性，即ESV值高的區(qū)域空間分布集聚，反之亦然。Global Moran’s I指數(shù)呈現(xiàn)出逐年下降趨勢，1990年為最大值，0.661，研究區(qū)ESV的空間相關(guān)性逐年下降，這一變化體現(xiàn)出研究區(qū)草地和水域ESV價(jià)值較高的區(qū)域面積逐年下降。ESV動(dòng)態(tài)變化相關(guān)性為先降低再升高的變化趨勢，且在2010～2018年達(dá)到最大值，為0.598，這一現(xiàn)象主要是由于2010～2018年間土地利用類型變化劇烈。

ESV高值區(qū)主要分布在研究區(qū)的北部、西部和南部地區(qū)，該區(qū)域主要是分布著草地、耕地和水域資源ESV價(jià)值較高的地類資源豐富，是研究區(qū)社會(huì)經(jīng)濟(jì)發(fā)展水平較高的區(qū)域；ESV低值區(qū)主要分布在研究區(qū)的西部和東南部地區(qū)，該類型區(qū)域分布主要地類為未利用地和部分草地，ESV價(jià)值高的地類資源較少；研究區(qū)中部和東北部地區(qū)ESV分布較為離散，主要地類為ESV價(jià)值低的未利用地，且生態(tài)環(huán)境脆弱，多為難開發(fā)利用區(qū)域。表7，圖4

表7 1990～2018年塔里木河流域全局莫蘭指數(shù)Table 7 Changes in the Global Moran’s I in the Tarim River Basin from 1990 to 2018

圖4 1990～2018年塔里木河流域全局莫蘭指數(shù)Fig.4 Changes in the Global Moran's I in the Tarim River Basin from 1990 to 2018

3 討 論

塔里木河流域在發(fā)展經(jīng)濟(jì)的同時(shí)，應(yīng)重視草地面積的大幅減少對(duì)脆弱生態(tài)環(huán)境的沖擊，以及塔里木河上游耕地面積增加和水利設(shè)施增多造成水資源浪費(fèi)和下游生態(tài)供水量短缺等生態(tài)問題[33]。

通過對(duì)比研究區(qū)與新疆其他地區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化特征，發(fā)現(xiàn)變化趨勢不盡相同。1994～2016年塔里木河流域北緣生態(tài)系統(tǒng)服務(wù)價(jià)值表現(xiàn)為持續(xù)下降態(tài)勢[34]；以全疆為尺度，2007～2016年生態(tài)系統(tǒng)服務(wù)價(jià)值總體上表現(xiàn)為增長[35]。各地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值差異主要是由土地利用結(jié)構(gòu)差異性、地形、時(shí)間尺度和研究區(qū)尺度的不同而導(dǎo)致的。研究認(rèn)為塔里木河流域1990～2000年生態(tài)系統(tǒng)服務(wù)價(jià)值呈現(xiàn)出上升趨勢，2000～2018年則呈現(xiàn)出下降趨勢。白元等[36]基于貨幣量的價(jià)值評(píng)價(jià)法估算出2005～2010年塔里木河干流生態(tài)系統(tǒng)服務(wù)價(jià)值表現(xiàn)為上升的趨勢；黃湘等[37]基于謝高地等制定的中國陸地生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子估算了塔里木河流域生態(tài)系統(tǒng)服務(wù)價(jià)值，結(jié)果表明，流域1994～2005年的整體生態(tài)系統(tǒng)服務(wù)價(jià)值呈現(xiàn)出上升趨勢。在當(dāng)量因子法評(píng)估生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)，當(dāng)量因子是否根據(jù)研究區(qū)的實(shí)際情況進(jìn)行修正以及當(dāng)量因子修正方法的選取與參數(shù)選擇的不同是造成使用相同評(píng)估方法但評(píng)估結(jié)果存在差異的主要原因之一，且研究時(shí)段的選取也會(huì)造成生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估結(jié)果呈現(xiàn)出差異。

通過基于網(wǎng)格單元的空間自相關(guān)分析方法能夠從空間角度較好的表達(dá)生態(tài)系統(tǒng)服務(wù)價(jià)值對(duì)土地利用變化的響應(yīng)，利用生態(tài)系統(tǒng)服務(wù)價(jià)值LISA聚類圖，ESV分布變化特征，為探究土地利用演變對(duì)生態(tài)系統(tǒng)服務(wù)的影響提供一種可行的方式。但對(duì)造成流域土地利用變化與生態(tài)系統(tǒng)服務(wù)價(jià)值的變化背后的自然和社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)因素尚未開展有效的工作，且空間關(guān)聯(lián)性不明確。

4 結(jié) 論

4.1未利用地、草地和耕地是研究區(qū)面積占比較大的土地利用類型，三者之和占比均在85%以上。流域內(nèi)耕地面積呈現(xiàn)為持續(xù)增長趨勢，未利用地和建設(shè)用地呈現(xiàn)出波動(dòng)式增長趨勢；中、低覆蓋度草地表現(xiàn)為持續(xù)下降，水域表現(xiàn)為波動(dòng)式減少趨勢。其中建設(shè)用地的漲幅最高，達(dá)到74%，耕地次之，未利用地漲幅最低；減幅最大的土地利用類型為水域，未利用地和林地次之。

4.2研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值1990～2000年表現(xiàn)為輕微上升趨勢，2000～2018年表現(xiàn)為下降趨勢，總體流域生態(tài)系統(tǒng)服務(wù)價(jià)值呈現(xiàn)出下降趨勢。其中2010～2018年間ESV下降趨勢劇烈，流域內(nèi)草地面積下降趨勢加劇，ESV高值土地利用類型轉(zhuǎn)向低值土地利用類型頻繁。

4.3研究區(qū)1990～2018年ESV的增值區(qū)主要集中分布在研究區(qū)的東南部地區(qū)；ESV的損失區(qū)主要分布在研究區(qū)北部和西部地區(qū)，塔里木河流域生態(tài)系統(tǒng)服務(wù)價(jià)值處于損失狀態(tài)。