辛強(qiáng)，李兆富，李瑞娟，郭泰，吳敏，潘劍君

（南京農(nóng)業(yè)大學(xué)資源與環(huán)境科學(xué)學(xué)院，江蘇 南京 210095）

基于溫度植被干旱指數(shù)的華東地區(qū)AMSR-E土壤水分?jǐn)?shù)據(jù)的空間降尺度研究

辛強(qiáng)，李兆富*，李瑞娟，郭泰，吳敏，潘劍君

（南京農(nóng)業(yè)大學(xué)資源與環(huán)境科學(xué)學(xué)院，江蘇 南京 210095）

土壤水分一直是土壤學(xué)領(lǐng)域中較為活躍的研究內(nèi)容，是陸面過程與水循環(huán)的重要影響因素。以AMSR-E為代表的被動微波遙感技術(shù)的發(fā)展為土壤水分的研究提供了方便，但是其粗糙的空間分辨率限制了其在中小尺度內(nèi)的應(yīng)用。因此，本研究利用MODIS溫度產(chǎn)品MOD11A2和歸一化植被指數(shù)產(chǎn)品MOD13A3構(gòu)建了月時間尺度下的溫度植被干旱指數(shù)（TVDI）；其次，利用溫度植被干旱指數(shù)TVDI和土壤水分之間的線性負(fù)相關(guān)關(guān)系，對AMSR-E三級土壤水分反演產(chǎn)品進(jìn)行空間降尺度研究，獲取2003年連續(xù)月時間尺度下空間分辨率為1 km的土壤水分反演結(jié)果，并利用地面實(shí)測土壤水分?jǐn)?shù)據(jù)對反演結(jié)果進(jìn)行驗(yàn)證。地面實(shí)測土壤水分值與降尺度反演結(jié)果顯著相關(guān)，每月的線性相關(guān)決定系數(shù)均在0.8以上，表明降尺度后的土壤水分反演結(jié)果具有較高的精度，能夠用來表示土壤水分的分布特征。

AMSR-E；TVDI；華東地區(qū)；土壤水分；降尺度

辛強(qiáng)， 李兆富， 李瑞娟， 郭泰， 吳敏， 潘劍君. 基于溫度植被干旱指數(shù)的華東地區(qū)AMSR-E土壤水分?jǐn)?shù)據(jù)的空間降尺度研究［J］.農(nóng)業(yè)現(xiàn)代化研究， 2016， 37（5）: 956-963.

Xin Q， Li Z F， Li R J， Guo T， Wu M， Pan J J. Downscaling AMSR-E soil moisture data based on temperature vegetation drought index in Eastern China［J］. Research of Agricultural Modernization， 2016， 37（5）: 956-963.

土壤水分是描述地表和大氣之間能量、水分交換的關(guān)鍵參數(shù)，對于水文、生態(tài)以及農(nóng)林牧業(yè)等研究具有非常重要的意義［1-2］。隨著微波遙感技術(shù)的發(fā)展，為區(qū)域尺度的土壤水分信息的獲取提供了有效手段［3］，極大促進(jìn)了土壤水分的研究進(jìn)程。國內(nèi)外許多學(xué)者開始使用遙感手段進(jìn)行土壤水分研究［4-7］，其中以AMSR-E為代表的被動微波遙感反演土壤水分已經(jīng)成為目前最主要的研究方法［8］。然而，被動微波輻射數(shù)據(jù)空間分辨率較低，大部分傳感器空間分辨率為25 km［9］，導(dǎo)致地表空間信息丟失［10］，難以滿足中小尺度對于數(shù)據(jù)的精度要求［11］以及在農(nóng)業(yè)應(yīng)用中的精度要求［12］。因此，基于被動微波輻射計的土壤水分?jǐn)?shù)據(jù)降尺度研究以獲取高分辨率土壤水分?jǐn)?shù)據(jù)具有現(xiàn)實(shí)意義［13-14］。

事實(shí)上，國內(nèi)外學(xué)者已經(jīng)嘗試使用遙感手段進(jìn)行AMSR-E土壤水分的降尺度反演研究［1，15-20］。其中，Sandholt等［21］提出的溫度植被干旱指數(shù)TVDI在土壤水分的研究中應(yīng)用較為廣泛。但是在現(xiàn)有的溫度植被干旱指數(shù)TVDI與被動微波輻射計土壤水分?jǐn)?shù)據(jù)的降尺度研究中，由于數(shù)據(jù)等方面的限制，大多數(shù)研究只集中于某一天或者某幾天的時間尺度，如曹永攀等［1］利用溫度植被干旱指數(shù)和MODIS數(shù)據(jù)對AMSR-E的降尺度研究，還有部分年際尺度上的研究，但是缺少對于連續(xù)月時間尺度下的AMSR-E土壤水分?jǐn)?shù)據(jù)的降尺度研究，而月尺度下的土壤水分遙感監(jiān)測對于土壤干旱的研究必不可少［20］。另外一方面，土壤水分的季節(jié)性變化規(guī)律愈來愈得到重視，而季節(jié)性變化又是以月尺度變化為基礎(chǔ)的［22］，因此，利用溫度植被干旱指數(shù)TVDI進(jìn)行AMSR-E土壤水分?jǐn)?shù)據(jù)的月尺度下的降尺度研究是非常必要的。

目前國內(nèi)關(guān)于土壤水分的研究，大多數(shù)選擇西部及西北部干旱以及半干旱地區(qū)作為研究區(qū)域［23，24］，而對于華東地區(qū)土壤水分的研究比較少見，而在邴飛龍等［25］的研究中指出，華東地區(qū)土壤水分蒸散量較高，土壤水分變化劇烈，針對華東地區(qū)展開土壤水分的研究具有重要的意義。因此，本研究選擇華東地區(qū)作為研究區(qū)域，首先利用MODIS溫度產(chǎn)品MOD11A2和歸一化植被指數(shù)產(chǎn)品MOD13A3計算2003年連續(xù)月時間尺度下的溫度植被指數(shù)TVDI，然后利用溫度植被干旱指數(shù)TVDI與土壤水分之間的負(fù)相關(guān)關(guān)系［21］，通過AMSR-E土壤水分反演結(jié)果計算TVDI與土壤水分間的負(fù)相關(guān)參數(shù)，進(jìn)而得到空間分辨率為1 km的土壤水分降尺度結(jié)果，并利用中國氣象數(shù)據(jù)共享網(wǎng)提供的地面實(shí)測土壤濕度數(shù)據(jù)對反演結(jié)果進(jìn)行驗(yàn)證。驗(yàn)證結(jié)果表明基于AMSR-E土壤水分?jǐn)?shù)據(jù)的空間降尺度反演結(jié)果具有較高的精度。

1　研究區(qū)概況

本研究選取華東地區(qū)為研究區(qū)域（圖1a），包括山東省、江蘇省、安徽省、浙江省、江西省、福建省和上海市，面積為79.83萬km2，經(jīng)緯度范圍24°27′36″—38°14′24″N，113°32′24″—122°42′36″E。氣候?qū)儆趤啛釒駶櫺约撅L(fēng)氣候和溫帶季風(fēng)氣候，以淮河為界限，淮河以北為溫帶季風(fēng)氣候，以南為亞熱帶季風(fēng)氣候。華東地區(qū)土地覆被類型多樣，以耕地和林地為主，兼有草地等用地類型（圖1b）。其中，耕地主要分布在山東、江蘇、安徽以及浙江北部，林地主要位于江西、福建以及浙江南部區(qū)域。地勢的總體特點(diǎn)是南北高、中部低，根據(jù)地形地貌特點(diǎn)大致可分為三部分：北部為黃淮平原，中部為長江中下游平原，南部為低山丘陵地區(qū)，最大高程為2 018 m（圖1c），地形地貌差異顯著，土壤水分含量差異較大。

2　研究方法

圖1　華東地區(qū)行政區(qū)劃及土壤水分測量站點(diǎn)（a）、土地覆被類型（b）與高程（c）Fig. 1　Administrative area and ground-based soil moisture （a）， land cover type （b）， and elevation of Eastern China （c）

2.1AMSR-E土壤水分降尺度方法

本研究選擇的土壤水分?jǐn)?shù)據(jù)AMSR-E空間分辨率為25 km，并且在現(xiàn)有大尺度的土壤水分研究中已經(jīng)得到廣泛應(yīng)用［2，10，26］，但是其粗糙的空間分辨率難以滿足中小尺度的精度要求。為了保證華東地區(qū)AMSR-E土壤水分?jǐn)?shù)據(jù)的精度，本研究利用溫度植被干旱指數(shù)TVDI與土壤水分之間的負(fù)相關(guān)關(guān)系［21］進(jìn)行AMSR-E土壤水分?jǐn)?shù)據(jù)的降尺度研究。溫度植被干旱指數(shù)最早是由Sandholt等［21］基于植被指數(shù)和地表溫度的關(guān)系提出的。TVDI越大，土壤濕度越低，反之，土壤濕度越高［27］，也就是TVDI與土壤水分之間存在良好的線性負(fù)相關(guān)關(guān)系。溫度植被干旱指數(shù)TVDI由植被指數(shù)和地表溫度計算得到，只依靠圖像數(shù)據(jù)，其表達(dá)式為：

式中：Ts為任意像元的地表溫度，Tsmin為最小的地表溫度，對應(yīng)的是濕邊；Tsmax=a+ b×NDVI，Tsmax為最大的地表溫度，對應(yīng)的是干邊，a和b為干邊擬合方程的系數(shù)。在干邊上TVDI為1，濕邊上TVDI為0，對于每個像元，利用NDVI確定Tsmax，進(jìn)而計算溫度植被干旱指數(shù)TVDI。TVDI越大，土壤濕度越低；TVDI越小，土壤濕度越高，也就是TVDI與土壤水分呈負(fù)相關(guān)關(guān)系［1］。在1 km空間分辨率尺度上，TVDI與AMSR-E土壤水分的負(fù)相關(guān)關(guān)系表述為：

如能求得參數(shù)a和b，則可由公式（2）計算得到分辨率為1 km的土壤濕度。AMSR-E空間分辨率為25 km，由25×25個1 km的網(wǎng)格組成，本研究將像元采取平均計算，即25 km網(wǎng)格的土壤濕度是其中所有1 km網(wǎng)格的土壤濕度的平均值［1，10，11］，公式表示為：

公式（3）中f （tvdi1，tvdi2，…，tvdii，…）則為空間分辨率為25 km的溫度植被干旱指數(shù)TVDI在1 km尺度上的平均值［1，21］，即：

由此，則可以通過公式（4）計算參數(shù)a與b，從而利用公式（2）得到空間分辨率為1 km的土壤水分降尺度反演結(jié)果。

2.2降尺度結(jié)果驗(yàn)證方法

本研究使用單點(diǎn)驗(yàn)證方法［28］驗(yàn)證土壤水分降尺度反演結(jié)果。AMSR-E土壤濕度產(chǎn)品及反演結(jié)果的精度驗(yàn)證的常用方法有四種［28］：第一，長時間土壤濕度地面觀測站驗(yàn)證。地面站點(diǎn)測得的結(jié)果都是單點(diǎn)數(shù)據(jù)，不能完全代表像元尺度內(nèi)的整體土壤濕度，因此驗(yàn)證結(jié)果會存在一定程度的誤差。但是由于觀測的規(guī)律性和持續(xù)性，地面觀測的土壤濕度數(shù)據(jù)可以與AMSR-E土壤濕度反演值比較變化趨勢，來驗(yàn)證長時間序列基于AMSR-E的反演結(jié)果的精度。這種方法省事省力，在本研究中就是選擇的該方法。第二，短期野外試驗(yàn)驗(yàn)證。野外試驗(yàn)由于試驗(yàn)區(qū)域有限，只能在少數(shù)幾個像元內(nèi)做精度驗(yàn)證，而地球表面的異質(zhì)性限制了該方法的應(yīng)用。第三，再分析資料土壤濕度比較。通過模型得到的土壤濕度在精度上比野外試驗(yàn)獲得的精度低，驗(yàn)證結(jié)果較差。第四，同衛(wèi)星不同傳感器，或者不同衛(wèi)星傳感器的類似產(chǎn)品做比較。該方法屬于交叉驗(yàn)證法，但是很難保證土壤水分結(jié)果的實(shí)時性。

華東地區(qū)地面測量站點(diǎn)測得的土壤水分結(jié)果為土壤相對濕度，而降尺度反演得到的土壤水分?jǐn)?shù)據(jù)為體積含水量。土壤水分體積含水量θv與質(zhì)量含水量θm之間的換算關(guān)系如公式（5）所示：

式中：sc為土壤容重（g/cm3）。

土壤相對濕度q為土壤質(zhì)量含水量θm與土壤田間持水量fc之比，如公式（6）所示：

結(jié)合公式（5）和（6），土壤體積含水量θv為土壤相對濕度q與田間持水量fc以及土壤容重sc的乘積，即：

利用轉(zhuǎn)化后的地面實(shí)測土壤濕度數(shù)據(jù)對降尺度后的土壤濕度反演結(jié)果進(jìn)行驗(yàn)證，并選取以下幾個參數(shù)：相對誤差REVR，絕對誤差A(yù)BVR，均方根誤差RMSE以及最大偏差MaxE和最小偏差MinE，定義地面土壤水分為X，降尺度反演結(jié)果為X0，N是參與比較的樣本數(shù)，此處即地面土壤水分測量數(shù)目，公式定義如下：

3　數(shù)據(jù)來源與處理

3.1AMSR-E土壤水分產(chǎn)品

研究區(qū)土壤水分產(chǎn)品選擇AMSR-E三級地表土壤水分產(chǎn)品，該產(chǎn)品為基于Qp模型開發(fā)的雙通道反演的結(jié)果，其數(shù)值為土壤體積含水量，單位為m3/m3，每天的數(shù)據(jù)包括升軌和降軌數(shù)據(jù)，空間分辨率為25 km。本研究選擇以2003年為例，對研究年限內(nèi)每天的土壤水分升軌、降軌數(shù)據(jù)進(jìn)行拼接、裁剪處理并進(jìn)行算數(shù)平均［10］，得到華東地區(qū)每月的土壤水分?jǐn)?shù)據(jù)，進(jìn)而通過對2003年十二個月的土壤水分?jǐn)?shù)據(jù)進(jìn)行算數(shù)平均得到當(dāng)年的土壤水分?jǐn)?shù)據(jù)［10］。

3.2MODIS數(shù)據(jù)

本研究選擇MODIS提供的1 km空間分辨率的8日地表溫度產(chǎn)品MOD11A2以及以月為時間步長的歸一化植被指數(shù)產(chǎn)品MOD13A3來計算溫度植被干旱指數(shù)TVDI。由于MODIS地表溫度產(chǎn)品MOD11A2存在少量數(shù)據(jù)缺失的情況，為了保證數(shù)據(jù)的完整性，對地表溫度產(chǎn)品MOD11A2采用最大值合成每月的地表溫度產(chǎn)品，最大值合成法可有效剔除、替換非正常像元值，保留正常值［11］，同時在計算溫度植被干旱指數(shù)的過程中會進(jìn)一步剔除異常高溫值的影響，所以此處選擇最大值合成月時間尺度下的地表溫度。由此與歸一化植被指數(shù)MOD13A3時間周期統(tǒng)一，進(jìn)而計算月時間尺度下的溫度植被干旱指數(shù)TVDI。

3.3地面土壤水分?jǐn)?shù)據(jù)

利用中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)（http://cdc. nmic.cn/）提供的土壤濕度數(shù)據(jù)驗(yàn)證空間降尺度結(jié)果的準(zhǔn)確性。該土壤濕度數(shù)據(jù)來源于以農(nóng)業(yè)氣象觀測為基礎(chǔ)的土壤相對濕度數(shù)據(jù)，數(shù)據(jù)集名稱為“中國農(nóng)作物生長發(fā)育狀況資料數(shù)據(jù)集”。華東地區(qū)共有125個土壤水分地面測量站點(diǎn)（圖1a），測量的土壤水分結(jié)果為土壤相對濕度，對每月三期的土壤濕度觀測結(jié)果進(jìn)行平均計算，以此來表示每月土壤濕度觀測結(jié)果，與AMSR-E時間步長統(tǒng)一。地面測得的土壤水分?jǐn)?shù)據(jù)為土壤相對濕度，在轉(zhuǎn)化為土壤體積含水量過程中所需要的土壤容重和田間持水量數(shù)據(jù)來源于Land-Atmosphere Interaction Research Group at Beijing Normal University （http://globalchange.bnu.edu.cn/）。

4　結(jié)果與分析

4.1華東地區(qū)AMSR-E土壤水分降尺度結(jié)果

表1統(tǒng)計了2003年十二個月內(nèi)溫度植被干旱指數(shù)TVDI與AMSR-E土壤水分?jǐn)?shù)據(jù)的線性回歸參數(shù)及兩者的決定系數(shù)R2。R2最大值為七月的0.91，最小值為十月的0.79。十二個月內(nèi)溫度植被干旱指數(shù)TVDI與AMSR-E土壤水分?jǐn)?shù)據(jù)均表現(xiàn)出較好的線性負(fù)相關(guān)關(guān)系，因此可以利用該線性參數(shù)對AMSR-E土壤濕度數(shù)據(jù)進(jìn)行降尺度運(yùn)算，得到空間分辨率為1 km的土壤水分降尺度反演結(jié)果（圖2）。

表1　AMSR-E與TVDI的線性回歸擬合參數(shù)Table 1　Linear fitting parameters between AMSR-E and TVDI

結(jié)合華東地區(qū)2003年AMSR-E土壤水分的降尺度反演結(jié)果，土壤水分含量及時空分布具有較大的變化。其中在二月份，山東北部地區(qū)土壤水分含量明顯偏高，造成該現(xiàn)象的原因是由于在研究中所使用的地表溫度數(shù)據(jù)存在部分缺失，這對于溫度植被干旱指數(shù)具有一定的影響，因此在反演結(jié)果上存在一定誤差，使得在山東區(qū)域內(nèi)土壤水分含量變化較為明顯。三、四、五月份山東、安徽等地土壤水分含量普遍偏低，浙江、福建土壤水分含量較前幾個月呈明顯增加的趨勢，主要原因在于山東、安徽等地屬于我國地理意義上的北方，此時進(jìn)入春旱時期，土壤水分含量明顯減少。從七月份開始，浙江、福建等南方地區(qū)開始進(jìn)入伏旱時期，尤其是在八月份，土壤水分含量與前段時期相比呈現(xiàn)出減少的變化趨勢。雖然土壤水分在時空分布上存在較大變化，但是從整個華東地區(qū)來說，以浙江、福建等林地覆蓋下土壤水分含量要高于山東、安徽等地耕地覆被下的土壤水分含量。

根據(jù)每月降尺度后的土壤水分反演結(jié)果，計算得到華東地區(qū)2003年土壤水分反演結(jié)果。降尺度前后AMSR-E土壤水分反演結(jié)果如圖3所示。圖3中左圖是經(jīng)過降尺度運(yùn)算后的結(jié)果，右圖是降尺度之前的原始AMSR-E土壤水分結(jié)果。經(jīng)過空間降尺度運(yùn)算后，將AMSR-E土壤水分?jǐn)?shù)據(jù)空間分辨率提高至1 km。與AMSR-E相比，降尺度后的土壤水分反演結(jié)果能夠滿足中小研究尺度范圍內(nèi)的數(shù)據(jù)精度要求，提高了研究精度，擴(kuò)大了AMSR-E土壤水分?jǐn)?shù)據(jù)的應(yīng)用范圍及領(lǐng)域。

圖2　華東地區(qū)2003年每月AMSR-E土壤水分降尺度反演結(jié)果Fig. 2　Retrieved soil moisture based on AMSR-E for each month of 2003 in Eastern China

圖3　AMSR-E降尺度前后對比圖Fig. 3　Comparison between downscaled result and AMSR-E soil moisture data

降尺度運(yùn)算后的土壤水分反演結(jié)果較原始影像相比，土壤水分最大值減小，并且部分區(qū)域土壤水分含量也較低。出現(xiàn)這種現(xiàn)象的原因主要是由于傳感器本身的影像：原始AMSR-E土壤水分影像是根據(jù)亮溫結(jié)果反演而來，每天亮溫數(shù)據(jù)的變化對于AMSR-E最終的土壤水分結(jié)果有較大影響，而降尺度反演結(jié)果在計算過程中采用算術(shù)平均計算，能夠減少由于溫度變化造成的土壤水分的影響。因此，降尺度反演結(jié)果比原始影像土壤水分值偏低。但是，經(jīng)過地面實(shí)測土壤水分的驗(yàn)證，降尺度反演的結(jié)果是符合精度要求的。

4.2AMSR-E土壤水分降尺度結(jié)果驗(yàn)證

計算降尺度后的土壤水分結(jié)果與地面實(shí)測值之間的決定系數(shù)R2，REVR，ABVR，RMSE，MaxE，MinE，統(tǒng)計結(jié)果如表2所示。土壤水分降尺度反演結(jié)果與地面實(shí)測土壤水分?jǐn)?shù)據(jù)間的決定系數(shù)R2均在0.80以上，表明兩者之間具有較好的線性關(guān)系。降尺度反演值與實(shí)測值最大差值為0.025 m3/m3，最小差值0.009 m3/m3，分別占降尺度土壤水分結(jié)果的17.86%和6.92%，誤差參數(shù)最大的是六月份，此時均方根誤差0.081，其次為七月、八月，主要原因是此時處于夏季，降雨量大，地面實(shí)測土壤水分含量較高，而降尺度的土壤水分則是月平均后的結(jié)果，因此存在部分誤差。其余月份均方根誤差均在0.070左右，最小值為十一月份的0.062。從絕對誤差和相對誤差值以及偏差值中也可以得出同樣的結(jié)論，即在夏季反演精度較低，春秋季反演精度較高。

表2　降尺度土壤水分反演結(jié)果誤差分析Table 2　Error analysis about retrieved soil moisture data

土壤水分反演結(jié)果與地面實(shí)測土壤水分?jǐn)?shù)據(jù)之間的驗(yàn)證結(jié)果證明：降尺度土壤水分結(jié)果具有較高的精度，線性回歸擬合決定系數(shù)R2均在0.8以上，表明降尺度土壤水分結(jié)果與地面土壤水分?jǐn)?shù)據(jù)具有較好的線性相關(guān)性，降尺度結(jié)果能夠描述土壤濕度的變化情況。雖然驗(yàn)證結(jié)果表明反演得到的土壤濕度具有較高的精度，但是仍然存在部分誤差，十二個月內(nèi)的地面土壤水分值均高于反演結(jié)果，綜合原因，主要有以下三個方面：第一，驗(yàn)證方法的選擇。如前所述，地面測量數(shù)據(jù)為點(diǎn)數(shù)據(jù)，用以驗(yàn)證面數(shù)據(jù)會存在一定程度的誤差。但是由于觀測的規(guī)律性和持續(xù)性，地面觀測的土壤濕度數(shù)據(jù)可以與AMSR-E土壤濕度反演值比較變化趨勢［29］。第二，測量深度的不同。地面實(shí)測土壤水分?jǐn)?shù)據(jù)測量深度在10 cm左右，而利用AMSR-E反演得到的土壤水分?jǐn)?shù)據(jù)的測量深度在0-5 cm土層，地表的土壤水分運(yùn)動較深土層劇烈，土壤水分散失嚴(yán)重，因此深土層土壤水分含量大于淺土層土壤水分含量。第三，土地覆被類型對結(jié)果的影響。華東地區(qū)土地覆被類型多樣，以耕地和林地為主（圖1）。地面測量站點(diǎn)分布于不同的土地覆被類型之上，不同用地類型的蓄水保水能力各不相同，總體來說，林地的蓄水保水能力要高于耕地和草地，這對地面測量值和反演值也有一定程度的影響。綜上所述，雖然土壤水分反演結(jié)果與地面實(shí)測土壤水分之間存在一定誤差，但是反演值仍然具有較高精度，能夠用來描述華東地區(qū)土壤水分的變化情況。

4.3華東地區(qū)土壤水分分布特征

由華東地區(qū)2003年每月的土壤水分反演結(jié)果進(jìn)行算數(shù)平均，得到2003年的土壤水分降尺度結(jié)果，如圖4所示。降尺度結(jié)果反映了土壤濕度的空間分布規(guī)律，其格局十分明顯，整體呈現(xiàn)北部低，南部高，中間過渡的分布形態(tài)，其中土壤水分含量最低的地區(qū)位于山東地區(qū)，土壤水分含量在0.113 m3/m3，其次是安徽、上海以及江蘇部分地區(qū)，土壤水分含量最大值的區(qū)域位于福建南部以及浙江部分地區(qū)，土壤水分含量最大值為0.146 m3/m3。這種分布規(guī)律與華東地區(qū)的經(jīng)度地帶性規(guī)律相符，同時也與研究區(qū)內(nèi)土地覆被類型緊密相關(guān)。因?yàn)槭芄庹蘸图撅L(fēng)影響，區(qū)域降水量自北向南逐漸增加，蒸發(fā)量則逐漸減少，而降水量和蒸發(fā)量基本上共同決定了該區(qū)域土壤水分狀況，同時由于在研究區(qū)域內(nèi)，耕地主要位于華東地區(qū)北部，以山東、江蘇、安徽、上海等地為主，而浙江、福建則以林地為主，林地的蓄水保水能力遠(yuǎn)遠(yuǎn)高于耕地，同時能夠減少水分的蒸發(fā)［30］，因此造成了土壤水分在華東地區(qū)內(nèi)的整體分布格局為南高北低。

圖4　2003年華東地區(qū)土壤水分分布Fig. 4　Spatial distribution of soil moisture in Eastern China in 2003

5　結(jié)論

本研究選擇MODIS地表溫度產(chǎn)品MOD11A2和歸一化植被指數(shù)產(chǎn)品MOD13A3，構(gòu)建溫度植被干旱指數(shù)TVDI，并利用土壤水分與溫度植被干旱指數(shù)之間的負(fù)相關(guān)關(guān)系，通過AMSR-E三級地表土壤水分產(chǎn)品計算得到TVDI與土壤水分之間的負(fù)相關(guān)線性參數(shù)，進(jìn)而對AMSR-E進(jìn)行降尺度計算，得到空間分辨率為1 km的土壤水分反演結(jié)果，并利用中國氣象數(shù)據(jù)共享服務(wù)網(wǎng)提供的地面實(shí)測土壤水分?jǐn)?shù)據(jù)對降尺度結(jié)果進(jìn)行驗(yàn)證。驗(yàn)證結(jié)果表明降尺度結(jié)果與地面實(shí)測土壤水分之間具有較好的線性關(guān)系，兩者的決定系數(shù)均大于0.8。不僅如此，還計算了土壤實(shí)測值與反演值之間的絕對誤差、相對誤差、均方根誤差以及最大差值和最小差值，得到的計算結(jié)果顯示了相同的結(jié)論：土壤水分實(shí)測值與降尺度反演結(jié)果具有很好的線性關(guān)系，反演結(jié)果能夠滿足精度要求。研究區(qū)域土壤水分時空分布差異較大，總體來說土壤水分含量較高的區(qū)域位于浙江、福建、江西等林地面積較大的地區(qū)，而山東、江蘇、安徽、上海等以耕地為主的地區(qū)，土壤水分含量較低。

6　討論

本研究利用Sandholt溫度植被干旱指數(shù)TVDI與土壤水分之間的線性負(fù)相關(guān)關(guān)系，對AMSR-E土壤數(shù)據(jù)進(jìn)行了空間降尺度運(yùn)算操作，獲取了連續(xù)月時間尺度下的空間分辨率為1 km的土壤水分反演結(jié)果。目前，在現(xiàn)有的研究中對連續(xù)時間尺度的AMSR-E土壤水分的降尺度研究還較為缺乏，本研究則是對該方面的一個補(bǔ)充。在后續(xù)的研究中，考慮如何進(jìn)一步減小中間參數(shù)的誤差是一個具有現(xiàn)實(shí)意義的研究方向。因?yàn)樵谀壳坝嬎銣囟戎脖桓珊抵笖?shù)的過程中，使用的都是MODIS提供的地表溫度產(chǎn)品和歸一化植被指數(shù)產(chǎn)品，而MODIS提供的日時間尺度和本研究應(yīng)用的8日時間尺度的地表溫度產(chǎn)品均存在不同程度的數(shù)據(jù)缺失。采用何種方法獲取日時間尺度下的溫度植被干旱指數(shù)，進(jìn)而得到連續(xù)日時間尺度下的高空間分辨率的土壤水分結(jié)果，是后續(xù)研究中比較重要的延伸方向。

致謝：感謝National Snow and Ice Data Center提供的AMSR-E土壤水分?jǐn)?shù)據(jù)，美國宇航局提供地表溫度和歸一化植被指數(shù)數(shù)據(jù)，中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供的地面土壤水分測量數(shù)據(jù)，Land-Atmosphere Interaction Research Group at Beijing Normal University提供的土壤容重和田間持水量數(shù)據(jù)。

［1］ 曹永攀， 晉銳， 韓旭軍， 等. 基于MODIS和AMSR-E遙感數(shù)據(jù)的土壤水分降尺度研究［J］. 遙感技術(shù)與應(yīng)用， 2011， 26（5）: 590-597.

Cao Y P， Jin R， Han X J， et al. A downscaling method for AMSR-E soil moisture using MODIS derived dryness index［J］. Remote Sensing Technology and Application， 2011， 26（5）: 590-597.

［2］ Chen X Z， Su Y X， LI Y， et al. Retrieving China’s surface soil moisture and land surface temperature using AMSR-E brightness temperatures［J］. Remote Sensing Letters， 2014， 5（7）: 662-671.

［3］ D’urso G， Minacapilli M. A semi-empirical approach for surface soil water content estimation from radar data without a-priori information on surface roughness［J］. Journal of Hydrology， 2006，321（1/2/3/4）: 297-310.

［4］ Notarnicola C， Ventura B， Pasolli L， et al. Towards an Operational Daily Soil Moisture Index Derived from Combination of MODIS，ASAR and AMSR-E Data［M］. Ieee International Geoscience and Remote Sensing Symposium， 2010: 816-819.

［5］ Brocca L， Hasenaure S， Lacava T， et al. Soil moisture estimation through ASCAT and AMSR-E sensors: An intercomparison and validation study across Europe［J］. Remote Sensing of Environment， 2011， 115（12）: 3390-3408.

［6］ Han N L， Chen S B， Wang Z J. Retrieval of Bare Surface Soil Moisture from AMSR-E Data［M］. 2nd Ieee International Conference on Advanced Computer Control， 2010: 67-72.

［7］ 楊濤， 宮輝力， 李小娟， 等. 土壤水分遙感監(jiān)測研究進(jìn)展［J］.生態(tài)學(xué)報， 2010， 30（22）: 6264-6277.

Yang T， Gong H L， Li X J， et al. Research progress of remote sensing monitoring soil moisture［J］. Acta Ecologica Sinica， 2010，30（22）: 6264-6277.

［8］ Sridhar V， Jaksa W T A， Fang B， et al. Evaluating bias-corrected AMSR-E soil moisture using in situ observations and model estimates［J］. Vadose Zone Journal， 2013， 12（3）: 299-314.

［9］ Song C Y， Jia L， Menenti M. Retrieving high-resolution surface soil moisture by downscaling AMSR-E brightness temperature using MODIS LST and NDVI data［J］. Ieee Journal of Selected Topics in Applied Earth Observations and Remote Sensing， 2014，7（3）: 935-942.

［10］ Feng H， Liu Y. Trajectory based detection of forest-change impacts on surface soil moisture at a basin scale Poyang Lake Basin， China［J］. Journal of Hydrology， 2014， 514: 337-346.

［11］ Zhao W， Li A N. A downscaling method for improving the spatial resolution of AMSR-E derived soil moisture product based on MSG-SEVIRI data［J］. Remote Sensing， 2013， 5（12）: 6790-6811.

［12］ Crow W T， Wood E F. The value of coarse-scale soil moisture observations for regional surface energy balance modeling［J］. Journal of Hydrometeorology， 2002， 3（4）: 467-482.

［13］ Crow W T， Wood E F， Dubayah R. Potential for downscaling soil moisture maps derived from spaceborne imaging radar data［J］. Journal of Geophysical Research: Atmospheres （1984-2012），2000， 105（D2）: 2203-2212.

［14］ Carlson T. An overview of the “triangle method” for estimating surface evapotranspiration and soil moisture from satellite imagery［J］. Sensors， 2007， 7（8）: 1612-1629.

［15］ Tuttle S E， Salvucci G D. A new approach for validating satellite estimates of soil moisture using large-scale precipitation: Comparing AMSR-E products ［J］. Remote Sensing of Environment， 2014， 142（4）: 207-222.

［16］ Pellarin T， Louvet S， Gruhier C， et al. A simple and effective method for correcting soil moisture and precipitation estimates using AMSR-E measurements［J］. Remote Sensing of Environment， 2013， 136（5）: 28-36.

［17］ Zeng J Y， Li Z， Chen Q， et al. A simplified physically-based algorithm for surface soil moisture retrieval using AMSR-E data［J］. Frontiers of Earth Science， 2014， 8（3）: 427-438.

［18］ 張顯峰， 趙杰鵬， 包慧漪， 等. 基于AMSR-E與MODIS數(shù)據(jù)的新疆土壤水分協(xié)同反演與驗(yàn)證［J］. 土壤學(xué)報， 2012， 49（2）: 205-211.

Zhang X F， Zhao J P， Bao H Y， et al. Co-inversion and validation of large-area soil moisture based on MODIS and AMSR-E data［J］. Acta Pedologica Sinica， 2012， 49（2）: 205-211.

［19］ 馮徽徽， 劉元波， 楊風(fēng)亭. 高植被覆蓋區(qū)AMSR-E土壤水分?jǐn)?shù)據(jù)校準(zhǔn)方法［J］. 土壤通報， 2014， 45（4）: 816-822.

Feng H H， Liu Y B， Yang F T. Calibration of the AMSR-E soil moisture data in dense vegetation area［J］. Chinese Journal of Soil Science， 2014， 45（4）: 816-822.

［20］ 李明星， 馬柱國. 基于模擬土壤濕度的中國干旱檢測及多時間尺度特征［J］. 中國科學(xué): 地球科學(xué)， 2015， 45（7）: 994-1010.

Li M X， Ma Z G. Soil moisture drought detection and multitemporal variability across China［J］. Science China: Earth Sciences， 2015， 45（7）: 994-1010.

［21］ 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（2/3）: 213-224.

［22］ 戴智慧. 貴州坡地黃壤水分時空變異特征研究［D］. 重慶: 西南大學(xué)， 2009.

Dai Z H. Soil moisture spatial and temporal variation characteristics research in Guizhou yellow［D］. Chongqing: Southwest University， 2009.

［23］ 李洪建， 王孟本， 柴寶峰. 黃土高原土壤水分變化的時空特征分析［J］. 應(yīng)用生態(tài)學(xué)報， 2003， 14（4）: 515-519.

Li H J， Wang M B， Chai B F. Spatial and temporal characteristics of soil moisture dynamics in Loess Plateau［J］. Chinese Journal of Applied Ecology， 2003， 14（4）: 515-519.

［24］ 李云鵬， 司瑤冰， 劉朋濤， 等. 基于空間信息的內(nèi)蒙古農(nóng)業(yè)干旱監(jiān)測研究［J］. 干旱區(qū)資源與環(huán)境， 2011， 25（11）: 125-131.

Li Y P， Si Y B， Liu P T， et al. Inner Mongolia agricultural drought monitoring research based on spatial information［J］， Resources and Environment in Arid Areas， 2011， 25（11）: 125-131.

［25］ 邴龍飛， 蘇紅波， 邵全琴， 等. 近30年來中國陸地蒸散量和土壤水分變化特征分析［J］. 地球信息科學(xué)學(xué)報， 2012， 14（1）: 1-13.

Bing L F， Su H B， Shao Q Q， et al. Changing characteristic of land surface evapotranspiration and soil moisture in China during the past 30 years［J］. Journal of Geo-information Science， 2012， 14（1）: 1-13.

［26］ Njoku E G， Chan S K. Vegetation and surface roughness effects on AMSR-E land observations［J］. Remote Sensing of Environment，2006， 100（2）: 190-199.

［27］ 王純枝， 毛留喜， 何延波， 等. 溫度植被干旱指數(shù)法（TVDI）在黃淮海平原土壤濕度反演中的應(yīng)用研究［J］. 土壤通報， 2009，40（5）: 998-1005.

Wang C Z， Mao L X， He Y B， et al. Application of temperaturevegetation dryness index （TVDI） in estimation of soil moisture in the HHH Plain［J］. Chinese Journal of Soil Science， 2009， 40（5）: 998-1005.

［28］ Xie X M， Xu J W， Zhao J F， et al. Validation of AMSR-E soil moisture retrievals over Huaihe River Basin， in China［J］. Applied Mechanics and Materials. 2014， 507: 855-858.

［29］ 陳潔. AMSR-E土壤濕度產(chǎn)品在我國西北地區(qū)的精度驗(yàn)證［D］.北京: 中國氣象科學(xué)研究院， 2010.

Chen J. Validation of AMSR-E soil moisture products in Northwest of China［D］. Beijing: The Chinese Academy of Meteorological Sciences， 2010.

［30］ Huang L， Shao Q， Liu J. Forest restoration to achieve both ecological and economic progress， Poyang Lake Basin， China［J］. Ecological Engineering， 2012， 44（3）: 53-60.

（責(zé)任編輯：王育花）

Downscaling AMSR-E soil moisture data based on temperature vegetation drought index in Eastern China

XIN Qiang， LI Zhao-fu， LI Rui-juan， GUO Tai， WU Min， PAN Jian-jun
（College of Resources and Environmental Sciences， Nanjing Agricultural University， Nanjing， Jiangsu 210095， China）

Soil moisture is not only an important factor in water cycle and land surface process， but also hot topic in soil science research. As a representative of passive microwave remote sensing technology， AMSR-E provides an avenue for soil moisture study. It is， however， restricted by the small and medium scales involved. Therefore， MOD11A2 and MOD13A3 were used to build the temperature vegetation drought index （TVDI） at a monthly time scale， the linear negative relationship between TVDI and soil moisture was then employed to downscale the spatial resolution of AMSR-E. This final objective was to obtain soil moisture at a spatial resolution of 1 km for each month in 2003. The results were well validated by the ground measured soil moisture data. Decision coefficients in all months were above 0.8. It indicated that the soil moisture results by downscaling were accurate enough and it could be used to characterize spatial distribution of soil moisture.

AMSR-E； TVDI； Eastern China； soil moisture； downscaling

Fundamental Research Funds for the Central Universities （KYZ201522）； National Natural Science Foundation of China （41571171，41171071）.

LI Zhao-fu， E-mail: lizhaofu@njau.edu.cn.

13 December， 2015；Accepted 5 April， 2016

S152.7

A

1000-0275（2016）05-0956-08

10.13872/j.1000-0275.2016.0099

中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)（KYZ201522）；國家自然科學(xué)基金項(xiàng)目（41571171、41171071）。

辛強(qiáng)（1988-），男，山東濟(jì)南人，碩士研究生，主要從事土壤濕度遙感反演，E-mail：xinqiang0711@126.com；通訊作者：李兆富（1977-），男，山東新泰人，博士，副教授，主要從事資源環(huán)境遙感研究，E-mail：lizhaofu@njau.edu.cn。

2015-12-13，接受日期：2016-04-05