王征

(中煤航測(cè)遙感集團(tuán)有限公司,陜西 西安 710018)

0 引 言

土壤濕度(Soil Moisture)是表示土壤干旱程度的重要指標(biāo),是陸-氣系統(tǒng)水熱過(guò)程的重要紐帶,同時(shí)也是植被生長(zhǎng)的主要水分來(lái)源[1]。目前主要通過(guò)兩種手段獲得土壤濕度信息[2],第一種是田間實(shí)測(cè)法,即利用野外觀測(cè)站或者采樣點(diǎn)直接測(cè)量不同深度的土壤濕度值,該方法雖然可以準(zhǔn)確地獲得土壤濕度信息,但是獲得的基于站點(diǎn)尺度數(shù)據(jù),不能夠刻畫大尺度、大范圍內(nèi)的土壤水分時(shí)空變化。第二種是采用遙感數(shù)據(jù)對(duì)土壤濕度信息進(jìn)行反演,是當(dāng)前較為常見的方法。通常結(jié)合熱慣量法、植被指數(shù)法等方法反演土壤濕度。馬春峰等人[3]以MODIS為數(shù)據(jù)源,采用熱慣量法對(duì)黑河流域土壤濕度信息進(jìn)行估測(cè),但該方法對(duì)于日較差數(shù)據(jù)精度有較高要求,無(wú)法滿足遙感定量研究的需要。2003年姚春生研究發(fā)現(xiàn)[4],土壤水分決定著植被冠層溫度,而冠層溫度間接決定土壤供水情況,基于此,提出溫度植被干旱指數(shù)(TVDI)。至此,國(guó)內(nèi)外學(xué)者基于TVDI監(jiān)測(cè)土壤濕度開展大量研究。Cao等人[5]利用AVHRR和MODIS數(shù)據(jù)構(gòu)建了Ts-NDVI特征空間監(jiān)測(cè)蒙古高原土壤濕度并分析干旱的時(shí)空變化特征,研究發(fā)現(xiàn)Ts-NDVI特征空間是監(jiān)測(cè)區(qū)域干旱的有效方法。Zhang等人[6]利用NOAA/AVHRR數(shù)據(jù)提取TVDI對(duì)2006年夏季四川盆地土壤濕度進(jìn)行監(jiān)測(cè),研究發(fā)現(xiàn)TVDI是監(jiān)測(cè)區(qū)域干旱的有效手段。程夢(mèng)園等人[7]基于LST和NDVI數(shù)據(jù)構(gòu)建Ts-NDVI特征空間對(duì)2005—2020年香日德-柴達(dá)木河流域的土壤濕度信息進(jìn)行提取,研究發(fā)現(xiàn)Ts-NDVI特征空間反演該流域土壤濕度具有較好的適用性,同時(shí)TVDI與海拔表現(xiàn)負(fù)相關(guān),而與平均氣溫表現(xiàn)正相關(guān)趨勢(shì)?？祱虻热薣8]基于MODIS數(shù)據(jù)中的NDVI和LST數(shù)據(jù)構(gòu)建NDVI-LST特征空間對(duì)2000—2019年蒙古高原干旱程度進(jìn)行監(jiān)測(cè),研究發(fā)現(xiàn)該地區(qū)的干旱情況小幅度加重,其增長(zhǎng)速率為 0.000 1/a。眾多研究表明應(yīng)用TVDI可以實(shí)現(xiàn)土壤濕度的反演,但已有研究在長(zhǎng)時(shí)序土壤濕度動(dòng)態(tài)監(jiān)測(cè)和驅(qū)動(dòng)因素方面仍然存在不足,理論情況下TVDI特征空間能夠達(dá)到裸土到全植被覆蓋的地表情況,但在實(shí)際應(yīng)用中,研究區(qū)域地表覆蓋情況難以滿足,得到的干濕邊實(shí)際上只是特征空間內(nèi)部存在的邊,并非理論的干濕邊。陜西省是重要的糧食作物產(chǎn)區(qū),而土壤濕度作為農(nóng)作物生產(chǎn)過(guò)程中的核心要素,需要長(zhǎng)時(shí)序動(dòng)態(tài)土壤濕度觀測(cè)數(shù)據(jù),但是陜西省氣象站點(diǎn)分布較少,難以獲取土壤濕度信息。

因此,本研究基于2000—2020年NDVI和LST數(shù)據(jù)構(gòu)建改進(jìn)的雙拋物線NDVI-LST特征空間,通過(guò)對(duì)TVDI指數(shù)進(jìn)行計(jì)算從而監(jiān)測(cè)陜西省土壤濕度情況,并結(jié)合氣候因素(降水、氣溫、日照時(shí)數(shù)和風(fēng)速)探討與其相關(guān)關(guān)系,以期為陜西省旱情監(jiān)測(cè)提供決策支持。

1 研究區(qū)概況

陜西省位于我國(guó)西北部地區(qū),地理范圍為105°29′～111°15′E,31°42′～39°35′N之間[9],總面積約為2.0×105km2(圖1)。陜西省區(qū)域內(nèi)地形地貌豐富多樣,地勢(shì)表現(xiàn)為南北高,中間低,平均海拔約為 1 100 m。

圖1 研究區(qū)行政區(qū)劃圖

陜西省地處濕潤(rùn)地區(qū)和干旱區(qū)的交界帶,陜南為濕潤(rùn)區(qū),主要有秦嶺和大巴山系,海拔高度在200～1 200 m之間,水熱條件相對(duì)較好,年均氣溫約為14℃,年降水量達(dá)到 800 mm;關(guān)中位于中間,屬于半濕潤(rùn)區(qū),主要是由渭河干流及兩岸支流沉積形成,地勢(shì)較低,海拔高度僅在300～800 m內(nèi),年均氣溫約為12℃,年降水量達(dá)到 600 mm;陜北則為半干旱區(qū),包括榆林和延安地區(qū),海拔高度在900～1 600 m之間,年均氣溫約為9℃,年降水量約為 500 mm。與此同時(shí),陜西省地處長(zhǎng)江流域中上游地區(qū),屬于國(guó)家生態(tài)環(huán)境重點(diǎn)建設(shè)區(qū)域,但是陜西省也是中國(guó)水土流失問(wèn)題較為嚴(yán)重的地區(qū)之一,因此,監(jiān)測(cè)區(qū)域土壤濕度時(shí)空變化對(duì)于陸地表面能量平衡具有重要作用[10]。

2 數(shù)據(jù)與方法

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

本研究采用的遙感數(shù)據(jù)為美國(guó)航天航空局的MODIS數(shù)據(jù)產(chǎn)品(https://modis.gsfc.nasa.gov/data/),獲取數(shù)據(jù)為2000—2020年,空間分辨率為 1 km×1 km,時(shí)間分辨率為 8 d的MOD11A2的合成產(chǎn)品和 16 d的MOD13A2合成產(chǎn)品。利用MRT軟件進(jìn)行一系列預(yù)處理操作,如拼接、投影轉(zhuǎn)換和格式轉(zhuǎn)換等[11],為消除云、霧以及大氣氣溶膠等外界因素的影響,本研究采用最大值合成法(MVC)計(jì)算得到月最大NDVI和月最大值LST,最后采用Python進(jìn)行裁剪處理。

氣象數(shù)據(jù)采用中國(guó)氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)2000—2020年陜西省及相鄰31個(gè)氣象站點(diǎn)的逐日氣溫、降水量、平均風(fēng)速和日照時(shí)數(shù)數(shù)據(jù),采用MATLAB軟件對(duì)該數(shù)據(jù)集進(jìn)行清洗和異常值處理,對(duì)于空缺數(shù)據(jù)采用雙線性內(nèi)插方法對(duì)其進(jìn)行填補(bǔ)[12],最后采用ANUSPLIN插值方法對(duì)其進(jìn)行空間插值,得到2000—2020年逐年 1 km分辨率的年均氣溫(Tmean)、總降水量(PRE)、年均平均風(fēng)速(WIN)和年日照總時(shí)數(shù)(SSD)數(shù)據(jù)。

2.2 溫度植被干旱指數(shù)法(Ts-NDVI特征空間)

研究發(fā)現(xiàn)[13],當(dāng)研究區(qū)植被覆蓋達(dá)到從裸土到全植被覆蓋、土壤含水量達(dá)到從完全干旱到濕潤(rùn)的條件下,NDVI和地表溫度的分布圖表現(xiàn)出類四邊形形狀(如圖2所示)。四邊形四個(gè)頂點(diǎn)分別代表4種極限情況,AD為干邊,土壤含水量有效性最小,地表水分蒸發(fā)量也較小;BC為濕邊,所包含的土壤含水量最大,此時(shí)地表蒸散量屬于潛在蒸散量。從A到B點(diǎn),地表土壤水分蒸發(fā)急劇增大,從B到C點(diǎn),植被覆蓋度由裸土達(dá)到全植被覆蓋。

圖2 Ts-NDVI特征空間原理

一般來(lái)說(shuō),當(dāng)監(jiān)測(cè)區(qū)完全達(dá)到濕潤(rùn)狀態(tài)時(shí),地表土壤內(nèi)的水分也隨之偏大,單一時(shí)段下特征空間濕邊比理論值高;當(dāng)監(jiān)測(cè)區(qū)完全達(dá)到干旱狀態(tài)時(shí),地表土壤內(nèi)的水分也隨之偏小,其特征空間干邊比理論值低;為彌補(bǔ)這一問(wèn)題,本研究采用長(zhǎng)時(shí)序遙感數(shù)據(jù)對(duì)陜西省土壤濕度進(jìn)行監(jiān)測(cè),以此彌補(bǔ)單一時(shí)段下Ts-NDVI特征空間的弊端。序列數(shù)據(jù)越完整,其合成的特征空間就更符合理論特征空間。計(jì)算公式如下[14]:

(1)

式中:LST表示像元地表溫度;LSTmin和LSTmax分別表示相同NDVI值所對(duì)應(yīng)地表溫度的最小值和最大值;TVDI在0～1之間,其值越大表示越干旱。

本研究采用對(duì)干濕邊方程改進(jìn)后的特征空間來(lái)計(jì)算TVDI,其計(jì)算公式如下[15]:

LSTmin=a1+b1×NDVI+c1×NDVI2

(2)

LSTmin=a2+b2×NDVI+c2×NDVI2

(3)

其中:a1、b1、c1、a2、b2和c2分別為干濕邊方程的系數(shù)。

本文將TVDI作為衡量干旱程度的指標(biāo),為消除量綱,對(duì)NDVI和LST均歸一化處理,步長(zhǎng)為0.01,編寫代碼獲取同一NDVI下全部單元格對(duì)應(yīng)的LST最大值和最小值,從而確定Ts-NDVI特征空間圖,將不同時(shí)期特征空間的干濕邊進(jìn)行線性擬合,最終得到逐年TVDI值。

2.3 Theil-Sen Median趨勢(shì)分析和Mann-Kendall檢驗(yàn)

基于Theil-Sen Median趨勢(shì)分析結(jié)合Mann-Kendall趨勢(shì)檢驗(yàn)方法能夠分析長(zhǎng)時(shí)間序列數(shù)據(jù),該方法相較于傳統(tǒng)的一元線性回歸分析方法而言,具有無(wú)可比擬的優(yōu)點(diǎn),如無(wú)數(shù)據(jù)情況下能夠服從特定分布、可抵抗一定程度的數(shù)據(jù)誤差,其計(jì)算公式為[16]:

(4)

式中:Median為中值函數(shù);TVDIi、TVDIj分別為該像元在i和j年的TVDI數(shù)據(jù)。β>0,即TVDI呈上升狀態(tài);β=0,即TVDI呈穩(wěn)定狀態(tài);β>0,即TVDI呈下降狀態(tài)。

Mann-Kendall是一種非參數(shù)突變檢驗(yàn)方法,它能夠判斷趨勢(shì)的顯著程度,當(dāng)n>10時(shí),標(biāo)準(zhǔn)的正態(tài)系統(tǒng)變量通過(guò)下式計(jì)算[17]:

(5)

2.4 Hurst指數(shù)

Hurst指數(shù)是用于表征長(zhǎng)時(shí)序變量在未來(lái)變化的演變趨勢(shì),該方法已在水文、植被、氣象等領(lǐng)域得到廣泛應(yīng)用[18]。目前,Hurst指數(shù)常使用重標(biāo)極差(R/S),以此確定陜西省未來(lái)的干旱變化情況。有研究發(fā)現(xiàn)[6],采用R/S分析法得到的Hurst指數(shù)更真實(shí),計(jì)算公式為

均值序列:

(6)

累計(jì)離差:

(7)

極差:

(8)

標(biāo)準(zhǔn)差:

(9)

根據(jù)Hurst指數(shù)的H值能夠確定陜西省的TVDI時(shí)間序列的持續(xù)性。Hurst指數(shù)主要分為三種情況:當(dāng)0.5

3 結(jié)果分析

3.1 基于TVDI的陜西省干旱與氣候因素的時(shí)間變化特征

2000—2020年陜西省年均TVDI介于0.48～0.6之間,總體呈現(xiàn)波動(dòng)下降趨勢(shì),其增長(zhǎng)速率為-0.002/a(p<0.01),這表明陜西省土壤濕度整體呈現(xiàn)濕潤(rùn)狀態(tài)。其中,2010年TVDI均值最小,為0.492;2005年TVDI均值最大,為0.585。研究時(shí)段可以2010年為界限,劃分2000—2010年的由干往濕發(fā)展和2010—2020年的由濕往干發(fā)展。2000—2010年年均TVDI呈現(xiàn)下降趨勢(shì),其增長(zhǎng)速率為-0.004/a;而2010—2020年年均TVDI呈現(xiàn)上升趨勢(shì),其增長(zhǎng)速率為0.002/a。

對(duì)于氣候因素而言,平均氣溫、總降水、平均風(fēng)速和日照總時(shí)數(shù)總體呈現(xiàn)波動(dòng)上升的趨勢(shì),其增長(zhǎng)速率分別為0.004/a、1.732/a、0.005/a和4.527/a(圖3)。

圖3 TVDI與氣溫、降水、風(fēng)速和日照時(shí)數(shù)的時(shí)間變化趨勢(shì)

3.2 基于TVDI的陜西省干旱空間變化特征

由于陜西省地區(qū)氣候類型多樣,各地區(qū)植被覆蓋度也存在明顯不同,因此年均TVDI值不足以表示干旱空間變化特征。本研究結(jié)合Theil-Sen Median趨勢(shì)分析與Mann-Kendall檢驗(yàn)方法進(jìn)行空間疊加,對(duì)陜西省內(nèi)像元尺度下的旱情趨勢(shì)進(jìn)行分析(圖4左)。首先對(duì)所有像元的年均TVDI值進(jìn)行趨勢(shì)分析,獲取到趨勢(shì)值β的空間分布,再利用檢驗(yàn)方法得到TVDI變化顯著性空間分布,將計(jì)算得到趨勢(shì)值和顯著性進(jìn)行疊加分類,得到陜西省年均TVDI的5種類型變化趨勢(shì),如顯著改善、輕微改善、基本穩(wěn)定、輕微退化和嚴(yán)重退化。結(jié)果發(fā)現(xiàn),陜西省的干旱變化中,面積占比最大的類型是輕微改善(50.91%),占比最小的類型是嚴(yán)重退化(0.95%),這表明2000—2020年陜西省土壤濕度在空間上往濕潤(rùn)方向發(fā)展為主,其中,呈現(xiàn)持續(xù)改善的地區(qū)主要分布在延安市的吳起縣、安塞縣、延川縣、甘泉縣,渭南市的合陽(yáng)縣、蒲城縣,寶雞市的隴縣、千陽(yáng)縣,西安市的藍(lán)田縣、戶縣,漢中市的寧強(qiáng)縣等地區(qū),主要占陜西省總面積的16.31%。呈現(xiàn)退化趨勢(shì)的面積占陜西省總面積的18.13%,其中呈現(xiàn)嚴(yán)重退化的面積占比0.95%,主要分布在榆林市的定邊縣,寶雞市的太白縣等地區(qū)。

圖4 2000—2020年陜西省TVDI時(shí)空變化趨勢(shì)及未來(lái)趨勢(shì)分析

為揭示陜西省土壤濕度未來(lái)變化的持續(xù)性特征,本研究進(jìn)一步將趨勢(shì)值與Hurst指數(shù)進(jìn)行耦合,得到陜西省年均TVDI的5種類型未來(lái)變化趨勢(shì),包括持續(xù)下降、下降、上升、持續(xù)上升、穩(wěn)定不變(圖4右)。陜西省年均TVDI的Hurst指數(shù)值為0.441,大于0.5的像元數(shù)所占比例為25.70%,小于0.5的像元數(shù)所占比例為74.30%。結(jié)果表明陜西省大部分地區(qū)TVDI未來(lái)變化趨勢(shì)與過(guò)去相反。在空間分布上,陜西省有52.37%的區(qū)域未來(lái)干旱情況有好轉(zhuǎn)趨勢(shì)。陜北地區(qū)植被以低矮灌木為主,生態(tài)環(huán)境相對(duì)較為脆弱,所以干旱的未來(lái)趨勢(shì)無(wú)法確定。蒸發(fā)量、降水量和CO2濃度等因素都是影響植被生理效應(yīng)的主要因素,但是這些因素在未來(lái)存在較強(qiáng)的不確定性,這就使得未來(lái)干旱演變軌跡很難監(jiān)測(cè)。但是本文通過(guò)對(duì)TVDI的趨勢(shì)值和Hurst指數(shù)進(jìn)行耦合綜合分析,實(shí)現(xiàn)了陜西省未來(lái)干旱情況的監(jiān)測(cè)。

3.3 陜西省TVDI與氣候因素相關(guān)性分析

氣溫、降水、日照時(shí)數(shù)和風(fēng)速變化都是影響土壤濕度的重要影響因子,因此,本研究基于陜西省的氣象數(shù)據(jù)進(jìn)一步結(jié)合陜西省TVDI時(shí)序數(shù)據(jù)分析其與降水、氣溫、日照時(shí)數(shù)和風(fēng)速的相關(guān)性(圖5)。結(jié)果發(fā)現(xiàn),陜西省TVDI與總降水之間以負(fù)相關(guān)關(guān)系為主(-0.109),這表面在降水量增加的同時(shí),TVDI值變小,土壤濕度呈現(xiàn)增加趨勢(shì)。根據(jù)相關(guān)系數(shù)可看出,負(fù)相關(guān)占陜西省總面積的70.86%,主要分布在榆林市的橫山縣、神木縣、佳縣、米脂縣,漢中市的城固縣、寧強(qiáng)縣等地區(qū),降水是土壤水分的重要補(bǔ)給,能夠緩解干旱的加重,當(dāng)降水量逐漸減少時(shí),土壤中的水分含量也在逐漸減少,特別是淺層土壤。正相關(guān)占陜西省總面積的29.14%,主要分布在咸陽(yáng)市的乾縣、彬縣、永壽縣、禮泉縣,渭南市的白水縣等地區(qū),隨著降水量的增加,植被蒸騰作用導(dǎo)致植物需水量大幅度增加,也間接性導(dǎo)致這些地區(qū)土壤水分消耗加快,使得TVDI與降水出現(xiàn)了正相關(guān)關(guān)系。

圖5 2000—2020年陜西省TVDI與降水、氣溫、風(fēng)速和日照時(shí)數(shù)之間的空間相關(guān)性

陜西省TVDI與平均氣溫主要呈正相關(guān)(0.003),這表明當(dāng)氣溫升高時(shí),TVDI值變大,則土壤濕度逐漸降低。根據(jù)相關(guān)系數(shù)結(jié)果可知,正相關(guān)占陜西省總面積的50.85%,主要分布在榆林市的子洲縣、橫山縣、定邊縣,延安市的甘泉縣、安塞縣、子長(zhǎng)縣、延長(zhǎng)縣,渭南市的大荔縣、澄城縣等地區(qū),當(dāng)溫度高于植被適宜值時(shí),植被葉面蒸騰作用以及土壤水分蒸散量加快,使得該地區(qū)具有干旱化趨勢(shì)。負(fù)相關(guān)占陜西省總面積的49.15%,主要分布在漢中市的佛坪縣,西安市的周至縣、戶縣等地區(qū),適宜的溫度使得地表土壤水分不易于蒸發(fā),植被鎖水能力增強(qiáng),使得該地區(qū)愈發(fā)濕潤(rùn)。

陜西省TVDI與平均風(fēng)速之間主要呈正相關(guān)(0.006),這表明當(dāng)風(fēng)速增加時(shí),TVDI值變大,則土壤濕度逐漸降低。根據(jù)相關(guān)系數(shù)結(jié)果可知,正相關(guān)占陜西省總面積的52.92%,主要分布于榆林市的子洲縣、神木縣、定邊縣,延安市的甘泉縣、延長(zhǎng)縣、子長(zhǎng)縣,渭南市的大荔縣、澄城縣等地區(qū)。大風(fēng)會(huì)帶走土壤中的水分,使得該地區(qū)具有干旱化趨勢(shì)。負(fù)相關(guān)占陜西省總面積的47.08%,主要分布于漢中市佛坪縣,西安市戶縣等地區(qū),適當(dāng)?shù)娘L(fēng)速不僅帶來(lái)植被生長(zhǎng)所需要的二氧化碳,而且加快了植被葉面蒸騰促進(jìn)了植物根部吸收土壤水分。

陜西省TVDI與日照總時(shí)數(shù)主要呈正相關(guān)(0.118),這說(shuō)明隨著日照時(shí)數(shù)的增加,TVDI值變大,則土壤濕度逐漸降低。根據(jù)相關(guān)系數(shù)結(jié)果可知,正相關(guān)占陜西省總面積的71.58%,主要分布于榆林市的橫山縣、定邊縣、子洲縣,延安市和渭南市的大部分地區(qū)。隨著日照時(shí)間的增長(zhǎng),加快土壤水分的增加,再加之這些地區(qū)本身干旱缺水。負(fù)相關(guān)占陜西省總面積的28.42%,主要分布于漢中市的佛坪縣,西安市的戶縣等地區(qū),適宜的日照可以促進(jìn)植被光合作用,有利于植被的生長(zhǎng),植被鎖水能力增強(qiáng),使得該地區(qū)愈發(fā)濕潤(rùn)。

4 結(jié) 論

(1)從時(shí)間演變特征上分析1,陜西省21年的年均TVDI在0.48～0.6波動(dòng),增速為-0.002/a(p<0.01),整體呈現(xiàn)波動(dòng)下降趨勢(shì);由Theil-Sen Median趨勢(shì)分析與Mann-Kendall檢驗(yàn)疊加分析表明,空間上陜西省主要以輕微改善為主,說(shuō)明土壤濕度在空間上往濕潤(rùn)方向發(fā)展為主。

(2)在未來(lái)演變特征上,陜西省年均TVDI表現(xiàn)出較強(qiáng)的反向持續(xù)性,其中Hurst指數(shù)小于0.5的地區(qū)占74.30%,表明陜西省大部分地區(qū)TVDI具有相反的變化趨勢(shì)。通過(guò)采用Theil-Sen Median趨勢(shì)分析與Hurst指數(shù)結(jié)果疊加分析表明,陜西省大部分地區(qū)未來(lái)干旱情況總體呈現(xiàn)好轉(zhuǎn)趨勢(shì)。

(3)從TVDI與氣候因素相關(guān)性分析結(jié)果來(lái)看,TVDI與降水呈負(fù)相關(guān),與氣溫、風(fēng)速和日照時(shí)數(shù)呈正相關(guān),并且TVDI與降水、氣溫、風(fēng)速的相關(guān)系數(shù)分別為-0.109、0.003、0.006和0.118。