劉 婕, 王明霞, 毋兆鵬,2

(1.新疆師范大學(xué) 地理科學(xué)與旅游學(xué)院, 烏魯木齊 830054;2.新疆師范大學(xué) 干旱區(qū)湖泊環(huán)境與資源新疆維吾爾自治區(qū)重點(diǎn)實(shí)驗(yàn)室, 烏魯木齊 830054)

?

遙感模型支持下的精河流域綠洲表層土壤水分時(shí)空分布特征

劉 婕1, 王明霞1, 毋兆鵬1,2

(1.新疆師范大學(xué) 地理科學(xué)與旅游學(xué)院, 烏魯木齊 830054;2.新疆師范大學(xué) 干旱區(qū)湖泊環(huán)境與資源新疆維吾爾自治區(qū)重點(diǎn)實(shí)驗(yàn)室, 烏魯木齊 830054)

土壤水分對(duì)干旱區(qū)氣候變化、生態(tài)環(huán)境穩(wěn)定、綠洲農(nóng)業(yè)發(fā)展有著重要的影響作用。以新疆精河流域綠洲為研究區(qū)，基于Landsat TM和ETM+影像，借助修正后的溫度植被干旱指數(shù)(TVDIm)，進(jìn)而反演了研究區(qū)2002年和2011年春、夏、秋3季的表層土壤水分。結(jié)果表明：(1) 野外同步實(shí)測(cè)數(shù)據(jù)精度驗(yàn)證R2值為0.800 2，TVDIm能夠較好反映綠洲土壤水分分布狀況；(2) 研究區(qū)表層土壤水分的季節(jié)性變化明顯，即秋季>夏季>春季，近10年間土壤水分總體呈趨濕特征；(3) 研究區(qū)綠洲內(nèi)部的土壤水分狀況明顯優(yōu)于過(guò)渡帶及天然荒漠區(qū)；(4) 土地利用疊加分析指示，人類活動(dòng)對(duì)土壤水分的空間分布影響強(qiáng)烈。

遙感; 土壤水分; 時(shí)空分布; 精河流域

新疆光熱豐富，降水稀少，是典型的干旱區(qū)，形成了以沙漠、半沙漠、戈壁、草原、森林以及綠洲為主的自然景觀和綠洲經(jīng)濟(jì)。隨著全球變化研究的深入和發(fā)展，人類生存與發(fā)展對(duì)土地的開發(fā)利用以及土地覆蓋變化被認(rèn)為是全球環(huán)境變化的重要組成部分和主要原因[1-2]。區(qū)域土壤水分檢測(cè)不但對(duì)改善區(qū)域及全球氣候的研究具有重要意義，而且能夠有效地監(jiān)控土地退化和干旱狀況。大量學(xué)者的研究闡明了土壤水分變化對(duì)環(huán)境生態(tài)和氣候變化響應(yīng)的重要性，以及農(nóng)作物生長(zhǎng)發(fā)育和不同時(shí)期土壤水分狀況的關(guān)系，并且有學(xué)者建立了土壤干旱預(yù)警模型，更好地監(jiān)測(cè)土壤水分變化對(duì)干旱區(qū)農(nóng)業(yè)的影響[3-6]。本文通過(guò)對(duì)精河流域綠洲土壤水分狀況的研究，不僅可以了解預(yù)測(cè)研究區(qū)土壤水分的干濕程度，也可以有效地掌握地表水分的運(yùn)移狀態(tài)，對(duì)干旱綠洲的農(nóng)業(yè)活動(dòng)具有一定的指導(dǎo)作用，也為改善小尺度區(qū)域的氣候，以及改善精河流域綠洲生態(tài)環(huán)境打下基礎(chǔ)。

1　試驗(yàn)材料與方法

1.1研究區(qū)概況

研究區(qū)包含整個(gè)艾比湖湖區(qū)和精河流域綠洲，總面積達(dá)2 572.30 km2，地理位置為82°40′—83°10′E，44°20′—44°50′N，自南向北呈扇狀坡面，中部為沖積—洪積平原，北部為艾比湖(圖1)。研究區(qū)地處亞歐大陸腹地，遠(yuǎn)離海洋，干燥少雨，屬典型的溫帶大陸性干旱氣候，區(qū)內(nèi)降水量由西向東遞增[7]，年平均降水量約為251.9 mm，年平均氣溫7.8℃，年平均蒸發(fā)量2 447.1 mm[8]。

圖1研究區(qū)及土壤水分樣點(diǎn)分布

1.2數(shù)據(jù)及預(yù)處理

本研究所用遙感數(shù)據(jù)主要來(lái)自美國(guó)陸地衛(wèi)星Landsat TM和ETM+影像共6景。其中，Landsat TM包括2002年5月7日(春季)、2002年7月10日(夏季)和2002年9月28日(秋季)；Landsat ETM+包括2011年5月24日(春季)、2011年7月11日(夏季)和2011年9月13日(秋季)。所有影像云量均小于5%，經(jīng)過(guò)輻射定標(biāo)、幾何糾正及大氣校正。此外，Landsat ETM+于1999年4月15日發(fā)射升空后，2003年因校正器故障，使得影像數(shù)據(jù)出現(xiàn)部分重疊和丟失。本文主要運(yùn)用ALR方法，即通過(guò)建立兩幅圖像灰度值的回歸方程，將填充影像的色調(diào)向待修復(fù)影像靠近，完成影像修復(fù)[9]，能夠滿足反演需求。

為了分析修正后的溫度植被干旱指數(shù)(TVDIm)與研究區(qū)域降水量、氣溫的關(guān)系，本文選取了2002年和2011年中4—10月的月平均溫度，氣象數(shù)據(jù)來(lái)源于精河氣象站，多年月平均降水量數(shù)據(jù)來(lái)自于相關(guān)文獻(xiàn)[8]。

1.3土壤水分反演算法

Sandholt等[10]在研究遙感反演土壤濕度的過(guò)程中，基于Ts-NDVI特征空間與土壤濕度的等值線關(guān)系，提出了可以表示土壤相對(duì)濕度狀況的溫度植被干旱指數(shù)(TVDI)。TVDI是由Ts-NDVI特征空間得到的比值，計(jì)算公式可以表示為：

(1)

式中：Tsmin——最小地表溫度，對(duì)應(yīng)的是濕邊；Ts——任意像元的地表溫度；Tsmax——指定像元NDVI對(duì)應(yīng)的最高地表溫度，對(duì)應(yīng)干邊；Tsmin，Tsmax的取值可以通過(guò)對(duì)Ts-NDVI特征空間的干邊及濕邊進(jìn)行模擬得到：

Tsmin=a2+b2NDVI

(2)

Tsmax=a1+b1NDVI

(3)

式中：a1，b1——干邊擬合方程的系數(shù)；a2，b2——濕邊擬合方程的系數(shù)。

由于NDVI在植被覆蓋度較高的區(qū)域容易飽和，導(dǎo)致對(duì)高植被覆蓋區(qū)敏感性較低[11]，運(yùn)用TVDI法反演濕度時(shí)會(huì)對(duì)結(jié)果造成一定的誤差。改進(jìn)型土壤調(diào)整植被指數(shù)(Modified Soil-Adjusted Vegetation Index，MSAVI)與歸一化植被指數(shù)(Normalized Difference Vegetation Index，NDVI)相比較，MSAVI修正了NDVI在植被密閉區(qū)或者低植被覆蓋區(qū)估測(cè)的誤差，增加了依據(jù)實(shí)際情況選定的土壤調(diào)節(jié)系數(shù)，可以較好地消除或減弱土壤背景對(duì)植被指數(shù)的影響，提高計(jì)算的精度。

(4)

式中：ρNIR——近紅外波段的反射率；ρR——可見光紅色波段的反射率。

因此本文采用改進(jìn)型土壤調(diào)整植被指數(shù)MSAVI替代NDVI參與計(jì)算，構(gòu)建Ts-MSAVI特征空間，減小因植被指數(shù)計(jì)算造成的誤差。結(jié)合公式(1)—(4)得到最終的TVDIm計(jì)算公式如下：

(5)

TVDIm=0表示土壤水分最大，地表最濕潤(rùn)狀況，即濕邊；TVDIm=1表示土壤水分最小，地表最干旱的狀況，即干邊。TVDIm的取值范圍為0～1，TVDIm值表示土壤的土壤水分狀況，TVDIm值越大，對(duì)應(yīng)的土壤水分越低。

2　結(jié)果與分析

2.1精河流域綠洲土壤濕度反演結(jié)果驗(yàn)證

為驗(yàn)證反演結(jié)果，在研究區(qū)內(nèi)選取典型26個(gè)樣點(diǎn)，與衛(wèi)星同步對(duì)土壤表層(0—10 cm)水分進(jìn)行實(shí)地監(jiān)測(cè)。監(jiān)測(cè)所用儀器為瑞士產(chǎn)MSR145土壤自動(dòng)記錄溫濕儀，可連續(xù)自動(dòng)記錄2年數(shù)據(jù)。本研究在衛(wèi)星過(guò)境前10 d埋入，衛(wèi)星過(guò)境后5 d取出。所有觀測(cè)樣點(diǎn)均采用GPS進(jìn)行定位，研究區(qū)土壤監(jiān)測(cè)點(diǎn)分布如圖1所示。本研究反演所用TM影像數(shù)據(jù)空間分辨率為30 m，而野外觀測(cè)點(diǎn)是根據(jù)土地利用類型進(jìn)行布設(shè)，二者雖在空間上難以對(duì)應(yīng)，但由于研究區(qū)地勢(shì)平坦，因此可以嘗試使用點(diǎn)上采集的土壤的含水量數(shù)據(jù)來(lái)驗(yàn)證區(qū)域面上的數(shù)據(jù)。TVDIm和土壤含水量，經(jīng)回歸分析t檢驗(yàn)(α=0.05)，得出線性回歸方程達(dá)到顯著。土壤含水量為橫坐標(biāo)，以TVDIm值為縱坐標(biāo)，對(duì)TVDIm和土壤含水量的相關(guān)性進(jìn)行線性擬合。TVDIm與實(shí)測(cè)土壤含水量的擬合結(jié)果，如圖2所示。TVDIm與表層(0—10 cm)土壤含水量之間表現(xiàn)出顯著的相關(guān)性，隨著土壤含水量增大，TVDIm呈現(xiàn)減小趨勢(shì)。由于整個(gè)研究區(qū)的土壤含水量較低，散點(diǎn)集中在土壤濕度低值區(qū)域。對(duì)TVDIm和土壤含水量的線性擬合結(jié)果，R2值為0.800 2，說(shuō)明TVDIm能夠反映地表土壤水分分布狀況，作為指示地表土壤濕度的指標(biāo)比較合理。

圖2土壤含水率與TVDIm的相關(guān)性

2.2精河流域綠洲土壤濕度分布特征分析

在遙感圖像處理軟件的支持下，根據(jù)上述TVDIm模型和方法，利用Landsat TM影像，得到2002年和2011年春、夏、秋3季的土壤水分狀況反演結(jié)果。參照TVDI指數(shù)反演土壤干旱等級(jí)劃分標(biāo)準(zhǔn)[12]，將研究區(qū)TVDIm劃分為5個(gè)等級(jí)：極濕潤(rùn)(0≤TVDIm<0.2)、濕潤(rùn)(0.2≤TVDIm<0.4)、無(wú)旱(0.4≤TVDIm<0.6)、干旱(0.6≤TVDIm<0.8)、極干旱(0.8≤TVDIm≤1.0)。在研究中通過(guò)TVDI土壤水分反演模型的劃分標(biāo)準(zhǔn)對(duì)研究區(qū)TVDIm進(jìn)行等級(jí)劃分，結(jié)果如圖3所示。

圖3精河流域土壤濕度反演結(jié)果分級(jí)

2.2.1精河流域綠洲土壤濕度時(shí)間分布特征由圖3可以看出，精河流域綠洲的土壤水分季節(jié)特征明顯，呈單峰型變化。從春季開始，TVDIm值開始明顯升高，表明土壤水分持續(xù)下降，連續(xù)最大的TVDIm值出現(xiàn)在夏季(6—8月)，反映出夏季土壤水分相對(duì)較低，進(jìn)入秋季TVDIm值則出現(xiàn)迅速減少趨勢(shì)，說(shuō)明研究區(qū)土壤水分呈現(xiàn)回升趨勢(shì)。從空間分布特征來(lái)看，春季，艾比湖湖區(qū)南部及西南部有部分呈斑塊狀區(qū)域土壤水分相對(duì)較高，其余區(qū)域土壤水分狀況均為干旱和極干旱。夏季，土壤水分較低的干旱區(qū)域呈片狀向四周擴(kuò)展，研究區(qū)除精河沿河區(qū)、內(nèi)沖積—洪積平原及綠洲耕作區(qū)土壤水分狀況為無(wú)旱外，其余地區(qū)為不同程度的干旱、極干旱。秋季，研究區(qū)除常年較干的南部、東南部外，整體土壤水分狀況開始大范圍好轉(zhuǎn)。

對(duì)精河流域過(guò)去近50年的年內(nèi)降水研究表明精河流域的年內(nèi)降水分配極不均勻，具有明顯的集中性[8]。將本研究的TVDIm數(shù)值與李加強(qiáng)等[8]多年平均降水量的月份分配數(shù)值疊加分析得出，研究區(qū)綠洲的土壤水分狀況的季節(jié)性變化與多年降水量呈明顯相關(guān)性(圖4)，這似乎與“降水量增長(zhǎng)則TVDIm值應(yīng)降低，即土壤水分升高”的常規(guī)性結(jié)論相悖。事實(shí)上這是由于研究區(qū)的特殊自然環(huán)境特點(diǎn)造成的。精河流域多年平均蒸發(fā)量2 447.1 mm，而多年平均降水量只有251.9 mm，相差近10倍，由此也表明，盡管降水量對(duì)土壤水分補(bǔ)充有一定影響，但不如其他氣象因子對(duì)土壤水分變化的影響強(qiáng)烈[13]。

圖4TVDIm與月平均降水量的關(guān)系曲線

進(jìn)一步將TVDIm數(shù)值與年內(nèi)逐月氣溫疊加分析，二者的趨勢(shì)一致性同樣十分明顯(圖5)。這說(shuō)明較高的氣溫會(huì)大大增加土壤蒸散發(fā)從而抵消降水對(duì)土壤水分的補(bǔ)給[14]。因此，對(duì)精河流域綠洲而言，在干旱區(qū)的大背景下，氣溫升高導(dǎo)致蒸發(fā)量增大的作用明顯大于降水量對(duì)土壤水分的補(bǔ)給作用，其影響起主導(dǎo)作用[7,13]。

圖5TVDIm與月平均氣溫的關(guān)系曲線

隨著全球氣候變暖，將引起蒸散發(fā)的增加，導(dǎo)致土壤干旱的加劇[8]。然而在這樣的大背景下，研究區(qū)2011年土壤水分狀況與2002年相比較好。這是因?yàn)榘群竦?007年由自治區(qū)級(jí)自然保護(hù)區(qū)晉升為國(guó)家級(jí)自然保護(hù)區(qū)。研究區(qū)內(nèi)退耕還林、退耕還濕、引水治沙、引水治草等的保護(hù)措施逐步實(shí)施，保護(hù)力度的逐步加強(qiáng)，使得土壤水分相對(duì)較高。

2.2.2精河流域綠洲土壤濕度空間分布特征盡管如此，由于人文作用會(huì)影響地表能量平衡，改變地表結(jié)構(gòu)及土壤水分的分布，因此，研究區(qū)人文活動(dòng)影響的疊加尤其是土地利用狀態(tài)，仍使地表土壤水分的狀況出現(xiàn)了較大的空間差異。本研究利用TM影像數(shù)據(jù)提取了2002年和2011年的土地利用數(shù)據(jù)(圖6)，并將TVDIm的月均值圖像與土地利用分類圖進(jìn)行疊加分析，得到不同土地利用條件下的土壤水分分布狀況(圖7)。

結(jié)果表明：耕地月平均TVDIm為0.65～0.85，而草地和林地月平均TVDIm為0.6～0.75，城鎮(zhèn)建筑用地月平均TVDIm為0.75～0.8，未利用地月平均TVDIm為0.8～0.9，土壤水分與土地利用之間關(guān)系密切(圖7)。

圖6精河流域綠洲土地利用狀況

耕地區(qū)和草地區(qū)的TVDIm月均值在4—7月逐漸升高，8月穩(wěn)定在0.8左右，表明植物處于生長(zhǎng)需水量較大的時(shí)期，9—10月TVDIm減小，土壤水分增加，表明植物已過(guò)生長(zhǎng)需水期或農(nóng)作物已收割、植物已枯萎，耗水量下降。

圖7不同土地利用類型TVDIm變化

林地區(qū)的TVDIm月均值，4—8月為0.65～0.82，與草地區(qū)TVDIm值相比增加趨勢(shì)平穩(wěn)，則是因?yàn)橐荒晟荼九c多年生木本需水量變化情況不同。

城鎮(zhèn)建筑用地的TVDIm月均值為0.75～0.8，未利用地的TVDIm月均值為0.8～0.9，且后者TVDIm值整體比前者高，表明土壤水分極少，但二者時(shí)間序列的土壤水分分布趨勢(shì)基本一致。

從研究區(qū)整體來(lái)看，耕地、草地、林地和城鎮(zhèn)用地的TVDIm值明顯低于未利用地，這也表明人類活動(dòng)對(duì)土壤水分的分布情況有很大的影響。

3　結(jié) 論

(1) 利用Ts-MSAVI模型能夠反映地表土壤水分分布狀況，可以作為有效指示地表土壤水分狀況的指標(biāo)。

(2) 精河流域綠洲的土壤水分季節(jié)特征變化明顯。春季氣溫升高，土壤水分降低，隨著夏季的到來(lái)，溫度持續(xù)升高，蒸發(fā)量大于降水量，研究區(qū)土壤水分普遍下降。秋季降水量減少，氣溫降低，蒸發(fā)量減少，土壤干旱狀況得到緩解。在干旱區(qū)的大背景下，相對(duì)于降水量而言，氣溫升高導(dǎo)致的蒸發(fā)量增大是影響土壤水分狀況的主導(dǎo)因素。

(3) 雖然耕地、草地和林地這3種土地利用類型的土壤水分在年內(nèi)的變化趨勢(shì)具有較好的一致性，但不同土地利用類型的土壤在不同的時(shí)期，由于人類活動(dòng)的影響，其TVDIm存在差異。對(duì)于生態(tài)環(huán)境極其脆弱的干旱區(qū)綠洲而言，在人口增長(zhǎng)帶來(lái)對(duì)耕地需求增加的背景下，如若不提高農(nóng)田管理質(zhì)量而導(dǎo)致大面積撂荒，就會(huì)極易引起土地鹽漬化或荒漠化，對(duì)當(dāng)?shù)厝嗣竦纳婧桶l(fā)展造成潛在威脅，繼而帶來(lái)一系列生產(chǎn)經(jīng)濟(jì)問(wèn)題。因此，在合理規(guī)劃土地利用，注重經(jīng)濟(jì)效益的同時(shí)，應(yīng)繼續(xù)推行退耕還林、退耕還草政策。

[1]李秀彬.全球環(huán)境變化研究的核心領(lǐng)域:土地利用/土地覆被變化的國(guó)際研究動(dòng)向[J].地理學(xué)報(bào),1996,51(6):553-558.

[2]Nunes C, Auge J L. Land-use and Land-cover Change(LUCC):Implementation strategy[M].1999.

[3]蔣沖,王飛,穆興民,等.黃土高原農(nóng)田土壤濕度演變及其與氣候變化的響應(yīng)關(guān)系[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(3):234-243.

[4]王曉賢,張學(xué)培.晉西黃土高原土壤水分垂直變化規(guī)律及生態(tài)影響[J].水土保持研究,2009,16(2):200-204.

[5]張哲,焦峰,梁寧霞,等.延河流域不同植被帶土地利用方式的土壤水分分析評(píng)價(jià)[J].水土保持研究,2014,21(4):74-79.

[6]蒲金涌,姚曉紅,辛昌業(yè),等.甘肅黃土高原土壤水分預(yù)測(cè)及旱情預(yù)警模型研究[J].干旱地區(qū)農(nóng)業(yè)研究,2010,28(3):254-258.

[7]王范霞,毋兆鵬.近40 a來(lái)精河流域綠洲土地利用/土地覆被時(shí)空動(dòng)態(tài)演變[J].干旱區(qū)資源與環(huán)境,2013,27(2):150-155.

[8]李加強(qiáng),陳亞寧,李衛(wèi)紅,等.天山北坡中小河川降水與徑流變化特征:以精河為例[J].干旱區(qū)地理,2010,33(4):615-622.

[9]壽敬文,陳雪,馬建文,等.采用ALR算法對(duì)Landsat-7圖像缺行修復(fù)的應(yīng)用研究[J].光電子·激光,2006,17(3):368-371.

[10]Sandholt I, Rasmussen K, Andersen J. A simple interpretation of the surface temperature vegetation index space for assessment of surface moisture status[J].Remote Sensing of Environment,2002,79(23):213-224.

[11]Qi J, Huete A R, Moran M S,et al. Interpretation of vegetation indices derived from multi-temporal SPOT images[J].Remote Sensing of Environment,1993,89(44):89-90.

[12]姚春生,張?jiān)鱿?汪瀟.使用溫度植被干旱指數(shù)法(TVDI)反演新疆土壤濕度[J].遙感技術(shù)與應(yīng)用,2004,19(6):473-478.

[13]王磊,文軍,韋志剛,等.中國(guó)西北區(qū)西部土壤濕度及其氣候響應(yīng)[J].高原氣象,2008,27(6):1257-1266.

[14]蘇明峰,王會(huì)軍.全球變暖背景下中國(guó)夏季表面氣溫與土壤濕度的年代際共變率[J].科學(xué)通報(bào),2007,52(8):965-971.

Spatiotemporal Distribution Characteristics of Soil Moisture in Surface Layer of Oasis in Jinghe Basin Under the Support of Remote Sensing Model

LIU Jie1, WANG Mingxia1, WU Zhaopeng1,2

(1.CollegeofGeographyandTourismScience,XinjiangNormalUniversity,Urumqi830054,China;2.MunicipalKeyLaboratoryofAridLakeEnvironmentandResources,XinjiangNormalUniversity,Urumqi830054,China)

Soil moisture has an important influence on the change of climate, the stability of ecological environment and the development of oasis farming in the arid area. Based on the Landsat TM and ETM+ images, we selected the Oasis of Jinghe Basin as the study site and investigated the moisture changes of the topsoil in spring, summer and autumn during the period from 2002 and 2011 using the modified Temperature-Vegetation Dryness Index (TVDIm). The findings are listed as follows： (1) By comparing and measuring, we found that theR2value of precision is 0.800 2, which means that TVDIm can better reflect the distribution of soil moisture in oasis； (2) The seasonal changes of moisture in the topsoil of the study area are significant, that is to say the moisture decreased in the order: autumn>summer>spring, and the water content is becoming higher in recent 10 years； (3) The condition of soil moisture in the inland of the study area is better than that in the transition belts and natural deserts； (4) The overlay analysis of land use indicated that human activities had the strong influence on the spatial distribution of soil moisture.

remote sensing; soil moisture; spatiotemporal distribution; Jinghe Basin

2015-06-02

2015-06-18

新疆維吾爾自治區(qū)自然科學(xué)資助項(xiàng)目(2013211A022)；新疆師范大學(xué)研究生科技創(chuàng)新項(xiàng)目(XSY201502006);新疆師范大學(xué)地理學(xué)博士點(diǎn)支撐學(xué)科(XJNU-DL-201515)；新疆維吾爾自治區(qū)重點(diǎn)實(shí)驗(yàn)室“新疆干旱區(qū)湖泊環(huán)境與資源實(shí)驗(yàn)室”基金(XJDX0909201206)

劉婕(1992—),女,河南駐馬店人,碩士研究生,研究方向?yàn)橘Y源環(huán)境地理信息系統(tǒng)。E-mail:liujie_ixj@163.com

毋兆鵬(1976—),男,陜西西安人,博士,教授,主要從事3S技術(shù)與綠洲資源開發(fā)研究。E-mail:wuzhaopengxj@sina.com

S152.7

A

1005-3409(2016)03-0095-05