雷凱文 盧宏瑋

摘要?基于氣象站點實測資料、MODIS植被產(chǎn)品數(shù)據(jù)集LAI（leaf area index）及世界土壤數(shù)據(jù)庫（HWSD）土壤數(shù)據(jù)，結(jié)合氣象插值軟件ANUSPLIN、土壤水分特征曲線擬合軟件RETC、空間信息處理軟件ArcGIS分別對氣象、土壤、植被數(shù)據(jù)進行計算，使用MIKESHE分布式水文模型對拉薩河流域內(nèi)2008—2017年3、7月及典型年（豐水年、中水年、枯水年）生長季（3—10月）內(nèi)的表層1 m土壤含水量進行模擬。結(jié)果表明，在空間分布上，拉薩河流域表層土壤含水量自西南向東北呈遞減趨勢，流域中部拉薩河中上游的高山地區(qū)土壤含水量相對較低，相較于中下游河谷及西部盆地地區(qū)，中部高山地區(qū)土壤水分變化更為劇烈。在時間分布上，典型年生長季內(nèi)不同月份土壤含水量變化量在33%～37%，總體而言，土壤含水量年內(nèi)變化大于年際變化。拉薩河谷平原、澎波盆地、流域西部念青唐古拉山脈以東的洪積寬谷、羊八井盆地及流域東北部麥地卡濕地植被覆蓋度和土壤水分含量較高、變化幅度小，農(nóng)業(yè)適宜性較強。

關(guān)鍵詞?土壤含水量，MODIS，MIKESHE，時空演變，拉薩河流域

中圖分類號?P334+.92文獻標(biāo)識碼?A

文章編號?0517-6611（2020）06-0050-05

Abstract?Based on meteorological site data，MODIS vegetation product dataset LAI （leaf area index） and soil data of the world soil database （HWSD），soil moisture content in the March and July during 2008-1017 in Lhasa River Basin was simulated and analyzed by the MIKESHE model，a distributed hydrology model.The results showed that in the spatial distribution，the surface soil water content across the Lhasa River Basin gradually decreased from the southwest to northeast，and the water content in the alpine region on the east and west sides of the middle stream was lower.Compared with the valley locate in the middlelower stream，the change of soil moisture in the central alpine area was more obvious.In terms of time distribution，the soil water content of the typical annual growth season changes in different months，varied between 33% and 37%，comparatively，the actual changes in soil moisture content were more obvious during the year.The basins and valley areas in the study area had higher vegetation coverage and higher soil moisture content，i.g.the Libo Basin and the Lhasa Valley Plain locate at middlelower stream of Lhasa River，Hongji Kuangu and Yangbajing Basins in the east of the Nyainqentanglha Mountains in the western part of the basin，and the Medica wetland in the northeastern part of the basin.The resulting agricultural suitability is stronger in these areas.

Key words?Soil moisture，MODIS，MIKESHE，Temporal and spatial variation，Lhasa River Basin

土壤含水量是氣候變化、水文循環(huán)、生態(tài)平衡過程中的一項重要參數(shù)，了解土壤含水量時空分布特征對區(qū)域旱澇預(yù)報[1-2]、農(nóng)牧業(yè)生產(chǎn)決策[3-4]與水資源調(diào)控[5-6]具有重要意義。目前常用的研究土壤含水量分布變化的方法主要包括實地測量法、遙感反演法、分布式水文模型模擬3種。劉雅莉等[7]通過在4種不同地形下設(shè)置采樣點并進行長期監(jiān)測與采樣，分析了微尺度下土壤水分的時空分布格局變化特征及其驅(qū)動因子，楊凱悅等[8]利用EM50數(shù)據(jù)采集器連續(xù)記錄監(jiān)測土壤體積含水量，對高寒沙區(qū)人工林在不同生長階段的土壤水分分布特征進行了研究。傳統(tǒng)的實地測量方法雖然可獲得精度更高的數(shù)據(jù)，但同時也存在數(shù)據(jù)獲取周期長、人工成本高的缺點，且站點網(wǎng)絡(luò)密度的限制使得傳統(tǒng)的站點監(jiān)測方法在區(qū)域尺度上的應(yīng)用被限制，難以用于大范圍區(qū)域土壤含水量的監(jiān)測[9]。隨著觀測技術(shù)與手段的不斷革新發(fā)展，基于衛(wèi)星遙感的土壤含水量反演技術(shù)獲得了很大發(fā)展，產(chǎn)生了以微波遙感法、熱慣量法、溫度-植被干旱指數(shù)法等土壤含水量的反演方法，Kolassa等[10]以ISMN站點監(jiān)測數(shù)據(jù)為基礎(chǔ)，對常用的土壤含水量反演產(chǎn)品SMOS、SMAP、MetopA及ESA CCI在不同區(qū)域及氣候條件下數(shù)據(jù)集的精度進行了評價。

使用遙感技術(shù)反演估算土壤含水量具有觀測快速、覆蓋區(qū)域廣的優(yōu)點，但同時觀測結(jié)果受云雨與復(fù)雜地表覆蓋類型的影響，通常難以獲得高精度、連續(xù)性強的土壤含水量[11]。由于遙感反演技術(shù)在高空間分辨率土壤含水率模擬方面的局限性，在利用衛(wèi)星遙感獲取的氣象、土地利用類型、LAI及DEM有效數(shù)據(jù)基礎(chǔ)上，使用分布式水文模型模擬目標(biāo)區(qū)域的土壤含水量在缺少觀測資料的地區(qū)具有廣泛的應(yīng)用前景。如Khan等[12]利用衛(wèi)星獲取非洲維多利亞湖某流域的氣象、土地利用類型等遙感數(shù)據(jù)，以此為基礎(chǔ)使用水文模型CREST（coupled routing and excess storage）對該流域的水文過程進行了模擬，Mason等[13]基于TerraSAR-X合成孔徑雷達(dá)（SAR）遙感數(shù)據(jù)，估測英國Tewkesbury城區(qū)的洪水淹沒區(qū)，研究認(rèn)為遙感獲取的洪水范圍能夠為水文模型的驗證提供較為可靠的信息。遙感技術(shù)為水文模型提供了大量基礎(chǔ)的地表水文信息，為缺乏實測資料區(qū)域的水文模擬提供了數(shù)據(jù)基礎(chǔ)，提高了定量描述流域地貌、水文過程的可能性?；谶b感數(shù)據(jù)的水文模擬在未來具有廣泛的應(yīng)用前景[14-15]。

該研究利用氣象站點監(jiān)測數(shù)據(jù)、HWSD土壤數(shù)據(jù)庫數(shù)據(jù)及MODIS遙感數(shù)據(jù)，結(jié)合分布式水文模型MIKESHE對拉薩河流域表層土壤含水量進行模擬，獲得2008—2017年流域表層土壤含水量，并對模擬結(jié)果與流域31個采樣點實測數(shù)據(jù)進行皮爾遜相關(guān)分析，選取各年枯水期與豐水期的空間分布結(jié)果，分析拉薩河流域表層土壤水分在不同情景下各年的空間分布狀況及其變化。

1?資料與方法

1.1?研究區(qū)概況

拉薩河流域地處雅魯藏布江流域中游左岸，面積33 823 km2，約占雅魯藏布江流域面積的13.5%。干流拉薩河全長約551 km，主要支流包括堆龍曲、墨竹馬曲、麥曲、拉曲、桑曲、和玉年曲等。拉薩河流域氣候以半干旱、半濕潤氣候為主，太陽輻射強烈，日照強烈，降水主要集中于夏季，氣候變化分明，從西南向東北大致分為溫和半干旱區(qū)、溫涼半干旱區(qū)、濕涼半濕潤區(qū)、冷涼半濕潤區(qū)4個氣候類型區(qū)[16]，流域內(nèi)旱澇并存，土壤水分空間分布存在較明顯差異。流域內(nèi)主要分布7種土類，分別為高山寒漠土、高山草甸土、亞高山草甸土、灌叢草原土、草甸土、潮土、新積土。中下游河谷區(qū)內(nèi)土壤類型主要為草甸土、草原土，土壤石礫含量較高。流域內(nèi)主要植被類型為高山沼澤草甸、高山灌叢草甸及少量高山灌叢[17-18]。

該研究區(qū)域位于29°20′～31°15′N、90°05′～93°20′E，主要包括拉薩市、墨竹工卡縣、林周縣、當(dāng)雄縣、嘉黎縣等主要縣市。其中拉薩河河谷區(qū)位于流域拉薩河中下游地區(qū)，是青藏高原最重要的農(nóng)業(yè)產(chǎn)糧地，主要產(chǎn)出糧食作物包括青稞、小麥、油菜等，表層土壤水分是當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)的重要制約因素。因此，對拉薩河流域內(nèi)土壤水分分布的研究對該地區(qū)農(nóng)業(yè)發(fā)展具有重要意義。

1.2?數(shù)據(jù)來源

氣象降水?dāng)?shù)據(jù)來自青藏高原地區(qū)317個氣象站監(jiān)測數(shù)據(jù)，使用Anusplin 42插值得到，葉面積指數(shù)數(shù)據(jù)（LAI）來自于美國航天局（NASA）MODIS15產(chǎn)品，時間分辨率為8 d，空間分辨率為500 m，流域土壤性質(zhì)數(shù)據(jù)來自聯(lián)合國糧農(nóng)組織（FAO）的HWSD土壤數(shù)據(jù)庫，使用基于Van-Genuchten模型的土壤水分特征曲線擬合軟件RETC處理后得到，流域內(nèi)土地利用類型數(shù)據(jù)來自中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心的基于Landsat-TM/ETM影像的中國多時期土地覆被遙感數(shù)據(jù)監(jiān)測數(shù)據(jù)集（CNLUCC），DEM數(shù)據(jù)來自中國科學(xué)院資源環(huán)境數(shù)據(jù)中心的全國1 km DEM數(shù)據(jù)產(chǎn)品。

根據(jù)多年降雨量的值，采用國內(nèi)常用的分析枯水年、平水年、豐水年標(biāo)準(zhǔn)：①豐水年，Pi>+0.33W，其中Pi為第i年降水量，W為均方差，②枯水年，Pi<+0.33W，③Pi與+0.33W結(jié)果相近時為平水年。以此為依據(jù)計算得到2008年為豐水年、2015年為枯水年、2012年為平水年。

1.3?ANUSPLIN氣象數(shù)據(jù)插值

ANUSPLIN允許引入多元協(xié)變量線性子模型，擁有平穩(wěn)處理二維以上樣條的功能，可引入多個氣象空間要素作為協(xié)變量，為考慮多要素下氣象數(shù)據(jù)的空間插值提供了可能。薄盤光滑樣條原型中擴展了局部薄盤光滑樣條（partial thin plate smoothing splines）部分[19]，允許引入海拔、降水、海岸線關(guān)系等線性協(xié)變量因子，擴展了模型功能。這種模型的表述如下：

式（1）中，xi為獨立變量矢量，f為要估算的關(guān)于xi的未知平滑函數(shù)，yi為p維獨立協(xié)變量矢量，Zi為空間上第i點的因變量，b為yi的p維系數(shù)，N為觀測值個數(shù)。當(dāng)p=0時（即協(xié)變量不存在時）即為普通薄盤樣條模型，當(dāng)不存在獨立自變量f（xi）時，模型變?yōu)槎嘣€性回歸模型。

由最小二乘估計來確定函數(shù)f和系數(shù)b：

式（2）中，Jm（f）定義為函數(shù)的m階偏導(dǎo)數(shù)，ρ為光滑參數(shù)（正值），在數(shù)據(jù)保真度與曲面的粗糙度之間起平衡作用。

1.4?RETC擬合土壤水分特征曲線

土壤水分特征為影響土壤含水量的重要因素，同時，土壤水分特征曲線也是分析土壤水分遷移的重要基本資料[20]。對缺乏資料地區(qū)，利用數(shù)值模擬方法對土壤水分運動進行預(yù)測，對于干旱程度、農(nóng)牧業(yè)適宜性評價等領(lǐng)域具有重要意義[21-22]。在推求土壤水分運動過程及擬合土壤水分特征曲線的研究中，Van Genuchten、Gardner、GardnerRusso、BroodsCorey等土壤水分特征曲線擬合模型被廣泛用于描述水分在土壤中的入滲、遷移、蒸發(fā)等過程。其中Van Genuchten模型適用范圍廣、擬合效果好、采用最為廣泛，其水分特征曲線方程為：

式（3）中，θ（h）為土壤含水量，h為壓力水頭（cm），θr為土壤剩余含水量，θs為土壤飽和含水率，α和n為經(jīng)驗擬合系數(shù)（或稱為土壤水分曲線形狀參數(shù)），m=1-1/n。

RETC軟件基于Broods-Corey模型和Van Genuchten模型對土壤水分特征曲線進行擬合，常用于非飽和帶土壤水力傳導(dǎo)特性和土壤水分的分析計算[23]。

1.5?MIKESHE分布式水文模型

MIKESHE（MIKE system hydrological european）模型是丹麥水利研究所在SHE模型基礎(chǔ)上研發(fā)而成。MIKESHE在濕地修復(fù)、洪水預(yù)報、水資源管理、水源地保護等流域具有廣泛的應(yīng)用。MIKESHE主要包括蒸散發(fā)（ET）、融雪（SM）、非飽和帶（UZ）、坡面流（OL）、河流湖泊（OC）、飽和帶（SZ）6個模塊，每個模塊模擬不同的水文過程，不同模塊間可以獨立作用也可相互作用[24]。MIKESHE模型可根據(jù)需要選擇單一模塊或耦合多個模塊進行建模[25]。該研究選用ET與UZ模塊對拉薩河流域地表土壤含水量進行模擬。近年來，對于缺乏資料地區(qū)的水文模擬逐漸成為國際水文水資源研究領(lǐng)域的難點與熱點[25]。這使得水文模型在缺乏資料地區(qū)的模擬應(yīng)用變?yōu)榭赡堋?/p>

MIKESHE蒸散發(fā)模塊及非飽和帶下滲模塊計算基于以下算法：

（1）Kristensen-Jensen模型。

采用Rutter模型模擬水循環(huán)中截留過程，該過程可計算蒸散發(fā)量、冠層實際蓄水量和到達(dá)地面的凈雨量。使用氣候與植被資料直接計算潛在蒸散發(fā)能力。Rutter模型的實質(zhì)功能是計算冠層蓄水量，其公式如下：

式（4）、（5）中，P是雨強（mm/s），S是冠層蓄水容量（mm），C是冠層實際含水量（mm），K、b是冠層蓄水參數(shù)，P1是地表植被覆蓋率，P2是總的葉片面積與植被覆蓋的地表面積之比， t是時間（s），Ep為土壤潛在蒸散發(fā)能力。

（2）植被實際蒸騰量。具體計算公式如下：

式（6）中，RDF為根系分布函數(shù)，f1（LAI）為葉片面積指數(shù)的函數(shù)，f2（θ）為土壤含水量函數(shù)。

對于發(fā)生在表層的土壤蒸發(fā)量，包括表層土壤蓄水量的蒸發(fā)量和基本蒸發(fā)量Epf3（θ）兩部分，其計算公式如下：

公式（8）、（9）中，θp為土壤田間持水量，θW為植物凋萎含水量，θM為土壤剩余含水量，C1、C2、C3為經(jīng)驗常數(shù)。

2?結(jié)果與分析

2.1?表層土壤含水量空間變化特征

拉薩河流域按地貌可分為和園區(qū)、高山盆地及河谷區(qū)3類，流域內(nèi)氣候類型屬于高原溫帶-寒溫帶半濕潤季風(fēng)氣候，年內(nèi)降水分布極不均勻，雨季一般為6月中旬之后開始，8月中旬左右結(jié)束，雨季降水量約占全年降水量的90%。利用拉薩河流域2008—2017年各年3月（枯水期）、7月（豐水期）降水、蒸散發(fā)、LAI數(shù)據(jù)，計算得到各年典型月（3、7月）表層土壤體積含水量及2008、2012、2015年的生長季3—9月表層土壤含水量。

由2008—2017年模擬結(jié)果可知，拉薩河流域表層土壤含水量自西南向東北呈逐漸減少趨勢，流域內(nèi)盆地地區(qū)平均土壤表層含水量相對較高，包括位于拉薩河中游林周縣境內(nèi)的澎波盆地、下游拉薩市、墨竹工卡縣中部境內(nèi)拉薩河谷平原地區(qū)，其含水量在雨季為36%～40%，旱季為34%～36%，流域西部當(dāng)雄縣境內(nèi)，念青唐古拉山脈以東的中部洪積寬谷盆地、南部羊八井盆地，其含水量在雨季為38%～40%，旱季為34%～36%，流域東北地區(qū)嘉黎縣北部麥地卡濕地，其含水量在雨季。流域中部自拉薩河上游兩岸向東西兩側(cè)的高山地區(qū)平均土壤含水量相對較低，其含水量在雨季為34%～36%，旱季為30%～34%。各典型年旱季和雨季表層土壤平均含水量見圖1。

從拉薩河流域不同年份雨季與旱季的表層土壤含水量變化量空間分布（圖2）可以看出，在各不同年份雨季全流域土壤表層含水量呈增加趨勢，旱季除東北部少量河源濕地地區(qū)外全流域土壤表層含水量呈下降趨勢。在雨季，流域內(nèi)除西部念青唐古拉山脈地區(qū)外，全流域土壤呈變濕潤態(tài)勢，其中南部河谷地區(qū)及西部高山盆地區(qū)域內(nèi)表層土壤含水量上升了2%～4%，流域東北部河源濕地地區(qū)表層土壤含水量上升了1%～3%，流域中部山地地區(qū)土壤含水量上升了3%～5%，相較而言流域中部山地地區(qū)濕潤化程度更為明顯，在旱季，除流域東北部少數(shù)區(qū)域外，全流域呈變干旱趨勢，其中南部河谷地區(qū)土壤含水量下降了1%～2%，東北部河源濕地地區(qū)下降基本小于1%，西部高山盆地地區(qū)下降了1%～3%，流域中部山地地區(qū)下降了1%～4%?？傮w而言，流域中部山地地區(qū)干旱化程度更為顯著。

2.2?表層土壤含水量時間變化特征

計算得到2008—2017年各年平均含水量及各年生長季內(nèi)土壤含水量變化量如圖3所示，結(jié)果顯示，各年雨季降水量充沛，雨季、旱季表層土壤含水量變化趨勢與年際間降水變化趨勢大致相同，在植物生長期內(nèi)各年7月末8月初平均表層土壤含水量最大，各年3、4月表層土壤含水量最小，在枯水年7月由于植被蒸騰作用茂盛而降雨量稀少，導(dǎo)致水分支出失衡，故表層土壤呈干旱化趨勢。

從圖4可以看出，典型年內(nèi)3—10月流域內(nèi)表層土壤含水量變化趨勢隨枯水年、平水年、豐水年基本呈上升趨勢。流域內(nèi)3—5月為旱季，降水較少，此時土壤水分變化較小，為33%～34%，7—8月為雨季，降水充沛，土壤含水量主要依靠降水補給，變化較大，為35%～37%。總體而言，土壤含水量年內(nèi)變化大于年際變化。

2.3?土壤含水量與氣象要素及LAI的相關(guān)性分析

土壤含水量受多種影響因子如氣象、植被覆蓋類型、高程、土壤成分、土層厚度等共同影響，作用機理復(fù)雜，各影響因子間亦存在相互促進或削弱的作用。例如降水、蒸發(fā)、氣溫等氣象因子直接影響流域內(nèi)土壤的補給水量，地表植被覆蓋類型的不同影響土壤蒸發(fā)量、植物蒸騰量、冠層截留水量及根系截留水量，高程及坡度影響流域內(nèi)壤中流的走向[22]，壤中流在重力作用下向海拔較低的地區(qū)匯流，故往往海拔較低的地區(qū)如流域內(nèi)的河谷、盆地地區(qū)的土壤水分含量相對較高，這也與該研究所得結(jié)果相吻合，土壤類型、土壤成分（黏土、礫石組分）、土層厚度等因素則會影響土壤的蓄水能力與導(dǎo)水能力。土壤水分的分布結(jié)果是各類影響因子共同作用的結(jié)果，該研究主要探索氣象、植被覆蓋對土壤水分含量分布的影響，利用SPSS 18統(tǒng)計分析軟件對拉薩河流域2008—2017年土壤水分含量與降水、蒸發(fā)、葉面積指數(shù)（LAI）進行皮爾遜相關(guān)分析（Pearson），結(jié)果表明（表1），拉薩河流域內(nèi)土壤水分含量與降水、LAI呈正向關(guān)關(guān)系，置信度較高，故可認(rèn)為降水及植被覆蓋度是流域內(nèi)土壤水分含量的重要作用因子。

3?結(jié)論與討論

（1）拉薩河流域表層土壤含水量自西南向東北呈逐漸減少趨勢，流域內(nèi)盆地地區(qū)平均土壤表層含水量相對較高，流域南部河谷、西部盆地及東北部濕地地區(qū)土壤水分含量高于流域中部山地地區(qū)，各年旱季中，中部山地地區(qū)干旱化趨勢強于流域中河谷、盆地地區(qū)，各年雨季中，流域中部山地地區(qū)濕潤化趨勢強于流域中河谷、盆地地區(qū)。拉薩河流域河谷及盆地地區(qū)表層土壤砂石含量相對中部山地較高，植被覆蓋度更高。因此造成的土壤蓄水、導(dǎo)水能力的不同，形成了土壤水分空間變化的差異?？傮w而言，高海拔的流域中部地區(qū)土壤濕度相對較低，更適宜安排耐旱作物或品種。

（2）拉薩河流域各年雨季土壤平均含水量在33.5%～35.5%。干旱年份雨季土壤水分含量增加較少。各典型年內(nèi)3—4月土壤含水量相對較低，變化幅度較小，5—6月降水充沛、植被相對茂盛的年份其土壤補給水量充沛、植被截留水量增加，相應(yīng)的土壤水分含量升高，在7月末至8月初為植被最為茂盛的階段，由于蒸騰作用劇烈導(dǎo)致土壤水分略有下降，8月后降水減少，氣溫下降、植被覆蓋度降低，蒸騰減弱，故流域內(nèi)土壤水分含量略有上升后持續(xù)降低?？傮w而言，土壤水分含量的年內(nèi)變化相較于年際變化更為劇烈，可將流域內(nèi)表層土壤水分的旱澇差異作為安排旱澇作物的依據(jù)。

（3）拉薩河流域尺度上土壤水分含量與降水、LAI呈正相關(guān)關(guān)系（P<0.05），與蒸發(fā)呈負(fù)相關(guān)關(guān)系（P>0.1），可能是因為降水及植物覆被直接影響土壤水分補給量，關(guān)系較為簡單，而蒸發(fā)受到地表植被蒸騰作用及土壤貯水、導(dǎo)水能力的影響，作用更為復(fù)雜，故而無明顯相關(guān)關(guān)系。之后的研究中可進一步加強對土壤水分及其主要影響因子的長時間序列監(jiān)測，以揭示更穩(wěn)定的土壤水分變化規(guī)律，這對調(diào)整布局農(nóng)作物生產(chǎn)具有積極意義。

參考文獻

[1] 郭燾，于紅博，馬梓策，等.基于MODIS的土壤含水量時空變化及干旱化程度分析[J].水土保持研究，2019，26（4）：185-189.

[2] BRONSTERT A，CREUTZFELDT B，GRAEFF T，et al.Potentials and constraints of different types of soil moisture observations for flood simulations in headwater catchments[J].Natural hazards，2012，60（3）：879-914.

[3] 張淑芳，柴守璽，常磊，等.冬小麥不同播期對土壤水分及產(chǎn)量的影響[J].中國農(nóng)學(xué)通報，2018，34（29）：13-19.

[4] 連玉珍，曹麗花，劉合滿，等.西藏農(nóng)牧交錯帶農(nóng)田與相鄰草地土壤含水量空間分布特征[J].西南農(nóng)業(yè)學(xué)報，2019，32（5）：1092-1097.

[5] 高海峰，白軍紅，王慶改，等.霍林河下游典型洪泛區(qū)濕地土壤pH值和土壤含水量分布特征[J].水土保持研究，2011，18（1）：268-271.

[6] BROCCA L，TULLO T，MELONE F，et al.Catchment scale soil moisture spatialtemporal variability[J].Journal of hydrology，2012，422/423：63-75.

[7] 劉雅莉，杜劍卿，李鋒，等.微尺度下城市公園人造綠地土壤水分的時空分異格局及其驅(qū)動機制[J].生態(tài)學(xué)報，2019，39（18）：6794-6802.

[8] 楊凱悅，賈志清，張立恒，等.高寒沙區(qū)典型人工林土壤水分空間分布特征研究[J].干旱區(qū)資源與環(huán)境，2019，33（5）：88-94.

[9] 汪倩倩，汪權(quán)方，王新生，等.地面資料稀缺區(qū)域的農(nóng)田土壤水分微波與光學(xué)遙感協(xié)同反演方法研究[J].中國農(nóng)學(xué)通報，2018，34（36）：117-123.

[10] KOLASSA J，REICHLE R H，LIU Q，et al.Estimating surface soil moisture from SMAP observations using a Neural Network technique[J].Remote sensing of environment，2017，204：43-59.

[11] WANG L L，QU J J.Satellite remote sensing applications for surface soil moisture monitoring：A review[J].Frontiers of earth science in China，2009，3（2）：237-247.

[12] KHAN S I，HONG ?Y，WANG J H，et al.Satellite remote sensing and hydrologic modeling for flood inundation mapping in Lake Victoria Basin：Implications for hydrologic prediction in ungauged basins[J].IEEE Transactions on Geoscience & Remote Sensing，2010，49（1）：85-95.

[13] MASON D C，SPECK R，DEVEREUX B，et al.Flood detection in urban areas using terraSAR-X[J].IEEE Transactions on Geoscience & Remote Sensing，2010，48（2）：882-894.

[14] GASHAW T，TULU T，ARGAW M，et al.Modeling the hydrological impacts of land use/land cover changes in the Andassa watershed，Blue Nile Basin，Ethiopia[J].Science of the total environment，2018，619/620：1394-1408.

[15] 徐宗學(xué)，程濤，洪思揚，等.遙感技術(shù)在城市洪澇模擬中的應(yīng)用進展[J].科學(xué)通報，2018，63（21）：2156-2166.

[16] GUO D L，F(xiàn)AN J，MI M X.Scaledependency of spatial variability of surface soil moisture under different land use types in Heihe Oasis，China[J].Chinese journal of applied ecology，2013，24（5）：1199-1208.

[17] 李筱金，徐琳瑜.拉薩河流域水電開發(fā)帶來的河岸帶土壤特征變化研究[J].環(huán)境科學(xué)與技術(shù)，2015，38（5）：148-156.

[18] 王慶海，李翠，龐卓，等.中國草地主要有毒植物及其防控技術(shù)[J].草地學(xué)報，2013，21（5）：831-841.

[19] HONG Y，NIX H A，HUTCHINSONB M F，et al.Spatial interpolation of monthly mean climate data for china[J].International journal of climatology，2005，25（10）：1369-1379.

[20] 劉貫群，孫蓓蓓，朱良超，等.內(nèi)蒙孿井灌區(qū)包氣帶土壤水分特征研究[J].中國海洋大學(xué)學(xué)報（自然科學(xué)版），2011，41（Z1）：107-112.

[21] 王全九，邵明安，鄭紀(jì)勇.土壤中水分運動與溶質(zhì)遷移[M].北京：水利水電出版社，2007.

[22] 陳雪，宋婭麗，王克勤，等.基于Van Genuchten模型的等高反坡階下土壤水分特征[J].水土保持研究，2019，26（5）：45-52.

[23] MA L，HE C G，BIAN H F，et al.MIKE SHE modeling of ecohydrological processes：Merits，applications，and challenges[J].Ecological engineering，2016，96：137-149.

[24] QIU Y，F(xiàn)U B J，WANG J，et al.Spatial prediction of soil moisture content using multiple-linear regressions in a gully catchment of the Loess Plateau，China[J].Journal of arid environments，2010，74（2）：208-220.

[25] 胡相明，趙艷云，程積民，等.云霧山自然保護區(qū)環(huán)境因素對土壤水分空間分布的影響[J].生態(tài)學(xué)報，2008，28（7）：2964-2971.