杜偉宏，吳天忠，霍艾迪，管文軻，韋 紅

(1.長(zhǎng)安大學(xué) 水利與環(huán)境學(xué)院，西安 710054； 2.長(zhǎng)安大學(xué) 旱區(qū)地下水文與生態(tài)效應(yīng)教育部重點(diǎn)實(shí)驗(yàn)室，西安 710054；3.新疆林業(yè)科學(xué)院，烏魯木齊 830000； 4.新疆塔里木河流域胡楊林生態(tài)系統(tǒng)國(guó)家定位觀測(cè)研究站，烏魯木齊 830000)

1 研究背景

土壤中水分不僅是陸地與大氣能量交換中的重要因素，同時(shí)也是表征土壤是否干旱的重要依據(jù)。在目前的生態(tài)環(huán)境研究中，大尺度的土壤濕度監(jiān)測(cè)已經(jīng)成為必不可少的研究?jī)?nèi)容，該項(xiàng)研究對(duì)一個(gè)區(qū)域的生態(tài)氣候、土壤水文等具有重大意義[1]。目前，國(guó)內(nèi)外專家學(xué)者主要通過3種方法獲取研究地區(qū)的土壤濕度，即田間實(shí)測(cè)法、土壤水分模型法、遙感監(jiān)測(cè)法，其中遙感監(jiān)測(cè)法監(jiān)測(cè)范圍廣、獲取數(shù)據(jù)快，已成為當(dāng)前監(jiān)測(cè)大尺度土壤濕度的重要手段[2-4]。

張仁華[5]以18 d重復(fù)一次的熱紅外資料為基礎(chǔ)，結(jié)合估產(chǎn)，通過改進(jìn)熱通量模型，提出了遙感信息的估產(chǎn)模式，奠定了我國(guó)遙感監(jiān)測(cè)土的基礎(chǔ)。楊鶴松等[6]以華北平原部分地區(qū)為例，利用2003—2006年的遙感數(shù)據(jù)、地區(qū)降水和土壤表層含水量數(shù)據(jù)進(jìn)行了歸一化指數(shù)和地表溫度的反演計(jì)算和驗(yàn)證。邢文淵[7]以新疆干旱地區(qū)為例，結(jié)合MODIS數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)，并以傳統(tǒng)的熱慣量模型為基礎(chǔ)，建立了適合新疆地區(qū)土壤濕度監(jiān)測(cè)的熱慣量法模型。魏珍[8]通過實(shí)測(cè)數(shù)據(jù)與Landsat7 ET，對(duì)陜北神木縣大柳塔煤炭開發(fā)區(qū)進(jìn)行了研究。楊東旭[9]利用MODIS衛(wèi)星遙感數(shù)據(jù)，結(jié)合實(shí)測(cè)土壤濕度數(shù)據(jù)對(duì)吉林中部地區(qū)的土壤濕度進(jìn)行了反演計(jì)算，結(jié)果表明表觀熱慣量與土壤濕度具有較好的相關(guān)性?；舭蟍10]通過提取MODIS數(shù)據(jù)中第6、第7波段的反射率，建立數(shù)值模型，驗(yàn)證了新疆干旱地區(qū)MODIS數(shù)據(jù)的第7波段反射率與土壤含水量呈負(fù)相關(guān)關(guān)系。

結(jié)合國(guó)內(nèi)外遙感反演土壤濕度的研究來看，大尺度的土壤濕度監(jiān)測(cè)主要利用TM與MODIS影像數(shù)據(jù)，且能較好地反映土壤濕度的變化。在塔里木河流域，土壤濕度的變化影響著流域區(qū)域的植被生長(zhǎng)與生態(tài)恢復(fù)?；舭蟍10]、韋紅等[11]、鐘家驊等[12]、杜偉宏等[13]通過研究發(fā)現(xiàn)，干旱地區(qū)土壤含水量與MODIS影響數(shù)據(jù)中的第7波段的反射率呈負(fù)相關(guān)；塔里木河干流胡楊林下土壤水鹽含量在空間上(地表至地下80 cm)不存在顯著差異；塔里木河流域的植被覆蓋度在19.696%～36.573%，存在空間異質(zhì)性。在此基礎(chǔ)上，本文通過實(shí)地測(cè)量的土壤濕度數(shù)據(jù)與MODIS衛(wèi)星遙感數(shù)據(jù)相結(jié)合的方法，驗(yàn)證了第7波段(Band 7)反射率可以用于該地區(qū)大尺度的土壤濕度監(jiān)測(cè)，建立了塔里木河流域土壤濕度反演模型，分析塔里木河流域近10 a來土壤濕度的變化情況，對(duì)生態(tài)恢復(fù)研究具有重要的研究意義。

2 研究區(qū)域概況

塔里木河干流位于新疆維吾爾自治區(qū)，是我國(guó)最大的內(nèi)陸河，全長(zhǎng)1 321 km，由阿克蘇河、葉爾羌河、和田河3條河流匯合成[12-13]。塔里木河自3條河流匯入處肖夾克向東流向塔克拉瑪干沙漠，途經(jīng)阿克蘇、沙雅、庫爾勒等地，終到巴音郭楞蒙古自治州若羌縣臺(tái)特瑪湖。研究表明，塔里木河流域綜合徑流補(bǔ)給主要為冰川融水，占總徑流補(bǔ)給量的47.9%，雨、雪水占27.9%、地下水占24.2%。塔里木河流域內(nèi)氣候干燥，降雨較少，蒸發(fā)強(qiáng)烈，多年平均降雨量為17.4～42.8 mm，多年平均蒸發(fā)量1 125～1 600 mm，最大風(fēng)速35 m/s，多年平均氣溫為10.6～11.5 ℃，日照時(shí)數(shù)在3 000 h左右，年最大日較差一般>30 ℃，平均日較差14～16 ℃，極端最高溫43.6 ℃，極端最低溫-30.9 ℃。

3 數(shù)據(jù)來源與研究方法

3.1 野外調(diào)查與土樣采集

為了研究塔里木河流域土壤濕度的時(shí)空變化特征，2019年9月27日—10月6日，研究團(tuán)隊(duì)在塔里木河中下游布設(shè)土壤監(jiān)測(cè)點(diǎn)(見圖1)，并對(duì)表層土(0～10 cm)進(jìn)行采樣。在每個(gè)100 m×100 m采樣點(diǎn)中設(shè)置5個(gè)25 m×25 m采樣單元，然后利用隨機(jī)布點(diǎn)法在采樣單元隨機(jī)選取5個(gè)點(diǎn)，多點(diǎn)采樣，每個(gè)采樣點(diǎn)稱取約50～60 g的土壤裝入鋁盒中(見圖2)。

圖1 研究區(qū)域及樣方地布設(shè)Fig.1 Sampling sites in the study area

圖2 土樣采集Fig.2 Collection of soil samples

野外取土后，將鋁盒取回的樣品放入烘箱，在105 ℃的恒溫下烘8 h，待樣品在烘干箱內(nèi)自然冷卻至常溫，使用分析天平稱得烘干后的土壤質(zhì)量。按式(1)計(jì)算土壤水分含量W，即

(1)

式中：W1為土樣濕質(zhì)量；W2為土樣干質(zhì)量。

3.2 遙感影像的選擇及處理方法

3.2.1 遙感影像的選擇

本研究采用Terra MODIS產(chǎn)品中的MOD09 Al。MOD09 Al數(shù)據(jù)空間分辨率為500 m(實(shí)際上是463.3 m)分辨率的第1—第7波段8 d合成的地表反射率，去除了云及云影干擾，可以用來監(jiān)測(cè)地球植被的季節(jié)變化和年際變化。

3.2.2 基于遙感監(jiān)測(cè)土壤濕度的變化特征

烘干土壤后，通過式(1)計(jì)算得出采樣樣品的土壤濕度，通過ArcGIS提取MODIS影像數(shù)據(jù)中實(shí)測(cè)采樣點(diǎn)位置的第7波段反射率數(shù)值，利用線性擬合的方法驗(yàn)證實(shí)測(cè)土壤濕度與MODIS影像數(shù)據(jù)中第7波段反射率之間的關(guān)系。上述擬合驗(yàn)證成功后，通過獲取的2010—2019年MODIS衛(wèi)星遙感數(shù)據(jù)，提取MODIS影像中10 a的塔里木河流域第7波段反射率，分析10 a來塔里木河流域土壤濕度的變化，通過塔里木河上中下游的3個(gè)河段，分析土壤濕度的空間變化。

3.2.3 MODIS數(shù)據(jù)波段監(jiān)測(cè)土壤濕度

前人通過研究發(fā)現(xiàn)，當(dāng)土壤的含水量較高時(shí)，水的吸收率曲線波長(zhǎng)值較低，且水的吸收率曲線在1.4、2.0 μm附近存在吸收峰[10]。MODIS衛(wèi)星影像數(shù)據(jù)中的第6波段的波長(zhǎng)為1.628～1.652 μm，第7波段的波長(zhǎng)為2.105～2.115 μm，和水的吸收峰接近。在荒漠化地區(qū)，土壤水分含量與MODIS影響數(shù)據(jù)中的第7波段反射率呈負(fù)相關(guān)，且通過對(duì)比研究發(fā)現(xiàn)，MODIS影像數(shù)據(jù)中第7波段的反射率比第6波段的反射率更接近荒漠地區(qū)土壤含水量的吸收峰[10]。故本次研究塔里木河流域土壤濕度遙感監(jiān)測(cè)采用Terra MODIS產(chǎn)品中數(shù)據(jù)空間分辨率500 m、8 d合成的MOD09A1第7波段來監(jiān)測(cè)土壤濕度。

4 結(jié)果與分析

4.1 塔里木河流域土壤濕度遙感監(jiān)測(cè)

利用MRT(MODIS Reprojection Tool)提取MODIS影像數(shù)據(jù)中第7波段，按照實(shí)測(cè)數(shù)據(jù)的位置在ArcGIS中輸入相同位置的經(jīng)緯度，提取該位置的第7波段反射率，并對(duì)比實(shí)測(cè)土壤濕度與MODIS第7波段反射率。通過對(duì)比分析發(fā)現(xiàn)，在塔里木河流域?qū)嶋H測(cè)得的土壤濕度與MODIS影像數(shù)據(jù)第7波段反射率呈一元線性回歸擬合關(guān)系(見圖3)，直線擬合公式為y=-0.45x+0.19,R2=0.48，即隨著反射率的增加，研究區(qū)域的土壤含水量逐漸減少。

圖3 2019年9月MODIS影像數(shù)據(jù)波段反射率與實(shí)測(cè)土壤濕度的擬合Fig.3 Fitted curve between band reflectance of MODISimage data and measured soil moisture inSeptember 2019

4.2 基于MODIS的土壤濕度時(shí)間變化特征

4.2.1 土壤濕度年內(nèi)動(dòng)態(tài)變化特征

通過對(duì)塔里木河流域2019年全年MODIS遙感數(shù)據(jù)的第7波段反射率提取分析，并在每個(gè)月中選擇同一時(shí)段內(nèi)的第7波段反射率代入上述一元線性擬合公式估算塔里木河上中下游及流域2019年各個(gè)月份的土壤濕度，得出2019年各個(gè)月份土壤濕度的動(dòng)態(tài)變化情況(見圖4)。

圖4 2019年塔里木河流域土壤濕度動(dòng)態(tài)變化Fig.4 Dynamic change of soil moisture in Tarim RiverBasin in 2019

由圖4可以看出，塔里木河流域土壤濕度年內(nèi)季節(jié)性變化較大，這與塔里木河流域的豐水期、枯水期相適應(yīng)。相關(guān)研究指出[14]，6—10月份塔里木河流域月徑流量占年徑流量的80.35%，1—5月份塔里木河的月徑流量?jī)H占年徑流量的7.74%。通過本次反演計(jì)算發(fā)現(xiàn)，塔里木河流域的土壤濕度在4.16%～6.45%之間，月平均土壤濕度最大為6.29%；2—5月份的土壤濕度均未超過5%，2019年最低土壤濕度僅為4.16%，出現(xiàn)在2月份，結(jié)合塔里木河的月徑流量來看，兩者情況是對(duì)應(yīng)的。

4.2.2 土壤濕度年際動(dòng)態(tài)特征

據(jù)了解，每年的6—10月份是塔里木河流域的汛期，本研究以塔里木河流域2010—2019年10 a的MODIS影響數(shù)據(jù)為基礎(chǔ)，選擇流域內(nèi)7、8月份及每年來水量較小的2、3月份進(jìn)行對(duì)比分析，利用在此期間的MODIS影響數(shù)據(jù)的第7波段反射率，反演計(jì)算出每年同期塔里木河流域的土壤濕度。

由圖5(a)可知，2010—2014年、2017—2019年塔里木河流域2月份的土壤濕度存在下降的趨勢(shì)，2014—2017年存在上升的趨勢(shì)。2010—2019年中土壤濕度最大值為7.20%；在2014年，塔里木河流域的土壤濕度最小，僅為3.49%。變異系數(shù)CV是衡量年際變化規(guī)律的常用參數(shù)，一般情況下，若CV>15%，說明研究數(shù)據(jù)震蕩變化相對(duì)劇烈。經(jīng)計(jì)算，在塔里木河流域2月份土壤濕度的年際間變異系數(shù)為20.91%，上、中、下游的土壤濕度年際間變異系數(shù)分別為18.59%、24.31%、22.77%。

圖5 塔里木河流域不同月份土壤濕度年際變化Fig.5 Interannual variation of soil moisture in TarimRiver Basin in different months

由圖5(b)可知，近10 a來，塔里木河流域3月份的土壤濕度變化總體平穩(wěn)，存在一定的震蕩、整體來看略有增加。最大土壤濕度出現(xiàn)在2017年，為4.54%；最小土壤濕度出現(xiàn)在2015年，僅為3.39%。流域內(nèi)3月份的土壤濕度年際間變異系數(shù)為11.78%，上、中、下游土壤濕度年際間變異系數(shù)分別為9.71%、15.18%、18.41%。

由圖5(c)可知，近10 a來塔里木河流域7月份土壤濕度整體變化趨穩(wěn)定，略有震蕩。最大土壤濕度出現(xiàn)在2012年，為7.53%；最小土壤濕度出現(xiàn)在2015年，僅為5.76%。流域內(nèi)3月份的土壤濕度年際間變異系數(shù)為9.04%，上、中、下游土壤濕度年際間變異系數(shù)分別為9.70%、7.94%、20.40%。

由圖5(d)可知，近10 a來，2010—2014年塔里木河流域8月份土壤濕度變化略有震蕩，2014—2017年土壤濕度呈增長(zhǎng)趨勢(shì)，2017—2019年土壤濕度呈下降趨勢(shì)。最大壤濕度出現(xiàn)在2017年，為7.71%；最小土壤濕度出現(xiàn)在2019年，僅為5.10%。流域內(nèi)8月份的土壤濕度年際間變異系數(shù)為12.10%，上、中、下游土壤濕度年際間變異系數(shù)分別為8.38%、2.92%、28.32%。

為了直觀表示塔里木河流域土壤濕度的變化情況，本研究以2010年8月份、2019年8月份該地區(qū)MODIS影響數(shù)據(jù)的第7波段反射率的空間分布為例，如圖6所示，其余年份的第7波段反射率空間分布不再贅述。

圖6 2010年8月份和2019年8月份塔里木河流域波段反射率空間分布Fig.6 Spatial distribution of band reflectance inTarim River Basin in August 2010 and August 2019

4.3 基于MODIS的土壤濕度的空間分布特征

由于塔里木河是由阿克蘇河、和田河、葉爾羌河三條河匯流而成的，且塔里木河干流狹長(zhǎng)，致使塔里木河上、中、下游的水文生態(tài)環(huán)境有所差異，沿程的土壤濕度情況也會(huì)存在一定的差異。本研究分別選取了2019年的3月份、7月份、10月份、12月份，分別代表流域內(nèi)春、夏、秋、冬的典型時(shí)段，利用MODIS影響數(shù)據(jù)及相關(guān)關(guān)系反演計(jì)算塔里木河流域不同河段的土壤濕度，反演結(jié)果對(duì)比見圖7。由圖7可見，塔里木河流域上、中、下游土壤濕度在年內(nèi)空間上、時(shí)間上差異較大。具體表現(xiàn)在以下幾個(gè)方面。

圖7 2019年塔里木河流域不同河段土壤濕度分布Fig.7 Soil moisture distribution in different reaches ofTarim River Basin in 2019

同時(shí)期不同河段對(duì)比：2019年3月份，塔里木河流域上游河段的土壤濕度分別是同期中、下游土壤濕度的1.85倍、2.75倍。隨著季節(jié)的變化，塔里木河流域來水量不斷增大，流域土壤濕度隨之增大，2019年10月份，塔里木河上、中、下游的土壤濕度變大，流域上游的土壤濕度為9.72%，分別是同期中、下游土壤濕度的1.75倍、2.94倍。通過對(duì)2019年塔里木河不同河段的土壤濕度分析比較，塔里木河流域土壤濕度表現(xiàn)為上游>中游>下游。

同河段不同時(shí)期對(duì)比：通過2019年3、7、10、12月份的土壤濕度反演計(jì)算，3—10月份，上游的土壤濕度增加了2.84%，中游土壤濕度增加了1.82%，下游土壤濕度增加了0.81%。通過同河段不同時(shí)期對(duì)比分析發(fā)現(xiàn)，塔里木河土壤濕度增加量表現(xiàn)為上游>中游>下游。

5 討 論

5.1 基于MODIS影像反演土壤濕度的差異性

土壤濕度是農(nóng)業(yè)干旱信息最重要的表征因子，是地球環(huán)境生態(tài)循環(huán)中的重要組成部分[15-16]。本文通過MODIS影像數(shù)據(jù)中的MOD09A1驗(yàn)證了塔里木河流域土壤濕度的反演計(jì)算模型，結(jié)果顯示實(shí)際測(cè)得的土壤濕度的與MOD09A1中第7波段反射率呈負(fù)相關(guān)關(guān)系，一元線性回歸擬合公式為y=-0.45x+0.19，R2=0.48，實(shí)測(cè)的土壤濕度與MODIS影像數(shù)據(jù)MOD09A1的第7波段反射率擬合相關(guān)系數(shù)較小。這主要有以下幾方面原因：

(1)衛(wèi)星影像數(shù)據(jù)與實(shí)際測(cè)得數(shù)據(jù)在時(shí)間上存在一定的差異。即使在本次實(shí)際取樣前15 d，研究區(qū)域沒有降雨，但MOD09A1是8 d合成的遙感影像數(shù)據(jù)，短短幾分鐘，衛(wèi)星就可能將研究區(qū)域過境，而現(xiàn)場(chǎng)取樣的土壤是幾天內(nèi)完成的，從而使得衛(wèi)星影像數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)時(shí)間不同步，可能造成一定誤差。

(2)衛(wèi)星影像數(shù)據(jù)的面平均數(shù)據(jù)與實(shí)際采樣點(diǎn)的點(diǎn)數(shù)據(jù)存在差距。本研究中衛(wèi)星遙感數(shù)據(jù)采用的是分辨率為500 m的MOD09A1，通過ArcGIS提取實(shí)測(cè)位置的第7波段反射率代表的是一個(gè)象元平均值，而實(shí)地測(cè)得的土壤濕度是在100 m×100 m樣方地內(nèi)。雖然樣方地內(nèi)多點(diǎn)取樣，測(cè)得土壤濕度后，計(jì)算樣方內(nèi)土壤濕度的平均值，單點(diǎn)數(shù)據(jù)的平均值也很難代表一個(gè)象元的平均土壤濕度。

(3)遙感數(shù)據(jù)反演計(jì)算與實(shí)測(cè)深度的土壤濕度數(shù)據(jù)關(guān)系不明確。本研究實(shí)測(cè)土壤濕度為地面10 cm 以內(nèi)，而MODIS影像數(shù)據(jù)的第7波段反映的是地面表層土壤濕度，二者的關(guān)系目前尚未明確，會(huì)產(chǎn)生相應(yīng)的誤差。

(4)不能精確地去除大氣影響。改進(jìn)的方法主要是盡量做到時(shí)間同步，提高點(diǎn)到面數(shù)據(jù)的代表性和建立更高精度的擬合模型。

基于以上，在荒漠化地區(qū)，MODIS反演土壤從點(diǎn)上和小尺度擬合效果并不是很理想。李薇[17]通過研究發(fā)現(xiàn)MODIS數(shù)據(jù)可用于監(jiān)測(cè)沙漠化地區(qū)大面積、多時(shí)相土壤含水量，可以反映出土壤含水量的變化。2010—2019年塔里木河流域土壤濕度的變化與塔里木河來水量及生態(tài)輸水有著密切的聯(lián)系，塔里木河流域土壤濕度在地表水的作用下發(fā)生變化。2018年5月，據(jù)塔里木河流域管理局統(tǒng)計(jì)，2017年為特豐年，上游阿拉爾來水量達(dá)68億m3，為近年來干流最大來水量，與反演土壤濕度年際變化中2017年土壤濕度最大的結(jié)果一致。綜上所述，MODIS影響數(shù)據(jù)中MOD09A1第7波段反射率可以用于反演、監(jiān)測(cè)塔里木河流域的土壤濕度。

5.2 塔里木河流域土壤濕度的時(shí)空分布

塔里木河流域土壤濕度年內(nèi)季節(jié)性變化較大，其中2019年6—10月份的土壤濕度超過6%，最大土壤濕度為6.29%；2—5月份的土壤濕度均未超過5%，全年最低土壤濕度僅為4.16%。近10 a來，2017年塔里木河流域土壤濕度最大，塔里木河2月份土壤濕度的數(shù)據(jù)變化較大，變異系數(shù)CV均超過15%。2010—2019年間，塔里木河下游的土壤濕度數(shù)據(jù)變化較大，CV均超過15%。

塔里木河流域上、中、下游土壤濕度在年內(nèi)空間上、時(shí)間上差異較大。通過同期不同河段的對(duì)比分析發(fā)現(xiàn)，塔里木河流域土壤濕度表現(xiàn)為上游>中游>下游。通過同河段不同時(shí)期對(duì)比分析發(fā)現(xiàn)，塔里木河土壤濕度增加量表現(xiàn)為上游>中游>下游。

2010—2019年塔里木河流域土壤濕度的變化與塔里木河流域來水量及生態(tài)輸水有著密切的聯(lián)系。近幾年來，新疆維吾爾族自治區(qū)啟動(dòng)了塔里木河流域胡楊林生態(tài)保護(hù)重大專項(xiàng)行動(dòng)，據(jù)塔里木河流域管理局統(tǒng)計(jì)，截至2016年秋，新疆已經(jīng)連續(xù)17次向塔河下游進(jìn)行生態(tài)輸水，累計(jì)輸送生態(tài)用水57.8億m3。在地表水的作用下，流域土壤濕度發(fā)生變化。從以上分析可以看出，塔里木河流域土壤濕度空間分異較大，且下游土壤濕度數(shù)據(jù)變化較大，雖然近幾年來塔里木河來水量較大，但如果僅靠天然輸水，塔里木河流域的土壤濕度分異難以調(diào)控，河道兩岸的生態(tài)環(huán)境難以恢復(fù)。因此，需要結(jié)合生態(tài)輸水，加強(qiáng)水資源人工調(diào)控，充分利用河岸生態(tài)閘縱向分水以及下游輸水，來增大流域的土壤濕度，滿足塔里木河流域兩岸植被的需水，對(duì)塔里木河流域生態(tài)恢復(fù)具有重要意義。

6 結(jié) 論

本研究表明，基于MOD09A1土壤濕度反演的驗(yàn)證，土壤濕度與Band 7波段反射率呈一元線性回歸擬合，擬合公式為y=-0.45x+0.19,R2=0.48。塔里木河流域土壤(0～10 cm)濕度年內(nèi)季節(jié)性變化較大，塔里木河下游的土壤濕度數(shù)據(jù)變化較大。塔里木河流域上、中、下游土壤(0～10 cm)濕度在年內(nèi)空間上、時(shí)間上差異較大，塔里木河流域土壤濕度表現(xiàn)為上游>中游>下游，土壤濕度增加量表現(xiàn)為上游>中游>下游。

本文采用的是野外實(shí)地調(diào)查、遙感定量反演和統(tǒng)計(jì)分析相結(jié)合的方法，該方法存在著一些不足，如本研究尚未考慮土壤取樣點(diǎn)的植被覆蓋度情況，且取樣深度為10 cm。我國(guó)于 2001 年開始接收 MODIS 影像數(shù)據(jù)。在今后的研究中，結(jié)合土壤取樣點(diǎn)的土壤濕度及植被覆蓋度情況，將取樣深度增加為20 cm，將時(shí)間序列適當(dāng)延展，分析塔里木河流域土壤濕度的變化特征，作出更全面的分析與預(yù)測(cè)。

致謝：感謝新疆林業(yè)科學(xué)院、塔里木河流域干流管理局和胡楊保護(hù)區(qū)管理局對(duì)野外調(diào)查的大力支持。