楊志輝,趙 軍,*,朱國鋒,2,劉家良,郭文镈,黃治化,王炎強(qiáng)

1 西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院, 蘭州 730070

2 中國科學(xué)院西北生態(tài)環(huán)境資源研究院,冰凍圈科學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 蘭州 730000

土壤水分(SM,soil moisture)是陸氣系統(tǒng)中的變量[1- 2],是水文循環(huán)和水熱平衡的組成部分[3],對地表植被生長也至關(guān)重要[4]。此外,SM對地表大氣溫度的演變[5]、降水的形成[6]和碳循環(huán)[7]有重要影響,對干旱事件的形成和發(fā)展也有重要的驅(qū)動作用[8]。因此土壤水分在氣候、水文、生態(tài)和農(nóng)業(yè)研究中起著重要的作用。

遙感技術(shù)的發(fā)展為大范圍、及時、準(zhǔn)確的SM遙感監(jiān)測提供了可能[9],成為檢測和監(jiān)測SM的主要手段[10]。早期SM光學(xué)遙感估算方法主要有光譜反射率法、短波紅外水分脅迫指數(shù)、植被狀況指數(shù)和溫度植被干燥指數(shù)等[11-13]。這些方法在估算SM時取得了較好的效果,但由于地表條件的差異,應(yīng)用于不同研究區(qū)也存在一定的差異。表觀熱慣量法適用于裸露土壤,溫度/植被指數(shù)法適用于植被覆蓋度較大區(qū)域[14],但由于不同植被覆蓋率的影響,降低了SM的估算精度。因此探索遙感估算SM的方法必須通過消除植被覆蓋影響來實(shí)現(xiàn)。消除混合像素中植被反射率貢獻(xiàn)的方法主要有：光譜混合分析法、光學(xué)植被覆蓋度法和梯形特征空間法[15]。但在實(shí)際應(yīng)用中光譜混合分析和光學(xué)植被覆蓋度法計算處理復(fù)雜,估算條件較難滿足。而梯形特征空間散點(diǎn)圖是由地表溫度/植被指數(shù)或反照率/植被覆蓋度的散點(diǎn)圖形成,對SM較為敏感且不依賴于地面觀測數(shù)據(jù),并可以與其他模型組合,所以梯形特征空間被廣泛應(yīng)用于估計植被蒸騰作用、土壤蒸發(fā)和土壤含水量[16]。

石羊河流域生態(tài)問題受廣大學(xué)者關(guān)注[17]。該流域農(nóng)業(yè)用水量比重較高,土壤水分的準(zhǔn)確估算對農(nóng)業(yè)灌溉用水的精確計算至關(guān)重要[18],因此,提高SM的估算精度、分析SM的空間格局是該流域目前亟待解決的問題之一?；诖吮疚囊肷Ⅻc(diǎn)圖梯形特征空間來消除植被的影響,提高SM估算精度,進(jìn)而分析空間格局,為土壤水分的有效利用提供依據(jù)。

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

1.1 研究區(qū)概況

石羊河流域位于甘肅河西走廊東部,祁連山北麓,介于101°22′—104°16′E,36°29′—39°27′N[19],屬河西走廊三大內(nèi)陸河流域之一,總面積4.16萬km2(圖1)。流域太陽輻射強(qiáng),溫差大,降水少,年降水量自南至北約為700—50mm、潛在蒸發(fā)量約為1200—2600mm,屬于典型的大陸溫帶干旱氣候,且氣候特征具明顯垂直分異。上游祁連山高寒半干旱濕潤區(qū),植被以高山和亞高山草甸為主,水土保持能力較強(qiáng)；中游走廊平原干旱區(qū),是重要的灌溉農(nóng)業(yè)區(qū)；下游溫暖干旱區(qū),主要為民勤綠洲和荒漠[20]。

圖1 研究區(qū)位置概況Fig.1 Location of the study area

1.2 研究數(shù)據(jù)

本文采用空間分辨率500m,時間分辨率8d的MODIS-MYD09A1數(shù)據(jù)的1—7波段計算地表反照率,從2018年11月至2019年10月共46期數(shù)據(jù),該數(shù)據(jù)能夠滿足生產(chǎn)活動對土壤濕度數(shù)據(jù)快速獲取和及時更新的要求。由空間分辨率為250m,時間分辨率為16d的MODIS-MYD13A1 NDVI產(chǎn)品數(shù)據(jù)計算植被覆蓋度(fractional vegetation cover,FVC),從2018年11月至2019年10月共24期。上述數(shù)據(jù)來源于NASA(https://ladsweb.modaps.eosdis.nasa.gov/)。

土地利用數(shù)據(jù)和土壤類型數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn)。

于2019年8月下旬用土鉆對石羊河流域0—10cm地表進(jìn)行土壤采樣,用手持GPS記錄樣點(diǎn)位置,樣點(diǎn)分布如圖1所示。將土壤樣本裝入鋁盒中并現(xiàn)場稱重,帶回實(shí)驗(yàn)室放入105℃的烤箱中持續(xù)烘烤8h,然后將其放置在干燥器中冷卻和稱重,持續(xù)操作至兩者差小于3mg,最后計算實(shí)際土壤含水量。實(shí)測SM數(shù)據(jù)用于構(gòu)建土壤水分估算模型和精度驗(yàn)證。

1.3 研究方法

1.3.1地表反照率

采用Liang Shunlin在2001年提出的MODIS表面反照率模型來計算地表反照率[21],該模型在反演華北土壤水分時也得到了充分應(yīng)用[22],公式如下：

a=0.16ρ1+0.29ρ2+0.243ρ3+0.116ρ4+0.112ρ5+0.08ρ7-0.0015

(1)

式中：a為全波段反照率,ρ1—ρ5、ρ7分別為MODIS相應(yīng)波段的地物反射率。

1.3.2建立梯形特征空間

梯形由asd,avd,avw和asw四個點(diǎn)形成(圖2)。其中asd和asw表示植被覆蓋率為零的最大和最小地表反照率。點(diǎn)avd和點(diǎn)avw表示植被覆蓋率為100%時最大和最小地表反照率。其中上邊界稱為“理論干燥邊”,表示在給定的植被覆蓋條件下土壤是非常干燥的。但在實(shí)際應(yīng)用中,“理論干燥邊”通常無法準(zhǔn)確獲得,取而代之的是“實(shí)際干邊”,表示植被覆蓋率變化時的實(shí)際表面反照率?！袄碚摑襁叀北硎驹谝欢ㄖ脖桓采w條件下土壤水分是充足的。同樣,在實(shí)際應(yīng)用中使用“實(shí)際濕邊”。位于“理論干邊”和“理論濕邊”之間的SM等值線便可以通過線性插值計算[23]。在SM等值線的基礎(chǔ)上可以容易計算裸土反照率。本文根據(jù)像元二分模型和歸一化植被指數(shù)計算研究區(qū)植被覆蓋度[24]。

圖2 表面反照率/植被覆蓋梯形特征空間示意圖Fig.2 Trapezoid space schematic diagram of surface albedo/vegetation coverage

1.3.3裸土反照率計算

研究區(qū)地表反照率主要取決于當(dāng)?shù)刂脖缓屯寥捞卣?。因?反照率是土壤和植被的加權(quán)總和[25],結(jié)合梯形特征空間得到裸土反照率計算公式,計算步驟詳見文獻(xiàn)[26]。

(2)

式中：asoil代表裸土反照率,a代表地表反照率。FVC表示植被覆蓋度。kd和kw分別表示“實(shí)際干邊”和“實(shí)際濕邊”的斜率。ai,max和ai,min表示植被覆蓋度較低時的實(shí)際地面最大和最小反照率。

1.3.4模型精度評估

計算SM估算值與實(shí)測土壤水分值的擬合度(R2),均方根誤差(RMSE)、平均絕對百分比誤差(MAPE)、不等式系數(shù)(TIC)、絕對誤差(AE)和相對誤差(RE)來評估土壤水分估算模型精度。此外,通過逐像元計算植被覆蓋度和土壤水分的Pearson′s相關(guān)系數(shù)來評估植被覆蓋對土壤水分空間分布的影響[27]。

1.3.5穩(wěn)定性計算

標(biāo)準(zhǔn)差是衡量數(shù)據(jù)分布程度的方法,可用于評估時間序列上SM的穩(wěn)定性[28]。公式如下：

(3)

1.3.6空間自相關(guān)

為分析整體尺度上SM在地理空間上是否存在集聚特征,采用全局Moran′sI(式(4))評估SM在空間上的自相關(guān)性。同時用局部空間相關(guān) Moran′sI(式(5))系數(shù)進(jìn)一步量化空間要素與周圍要素之間的差異程度及空間分布特征[29-30]。

(4)

(5)

1.3.7地理探測器模型

地理探測器主要包含以下4種: 因子探測,探測每個影響因子對 SM空間分異的影響大小,q值越高則對SM的影響力越大。風(fēng)險探測,探測因子對SM是否具有風(fēng)險性,指示因子在不同等級內(nèi)對SM分布的影響。交互作用探測,用于判斷不同影響因子對SM 的交互作用。生態(tài)探測,用于比較兩因子對SM空間分布的影響是否有顯著差異。分析過程詳見文獻(xiàn)[31]。q值是衡量一個因子對因變量的解釋力,表達(dá)式為：

(6)

2 結(jié)果與分析

2.1 生成梯形特征空間散點(diǎn)圖

Toby Carlson的研究已表明：沙漠可視為地表裸露的干土,反照率幾近于裸露的干土反照率；水體可視為地表裸露的濕土,反照率等于裸露的濕土反照率[32]?；诖瞬⒖紤]了不同植被覆蓋條件的地表反照率特征,在無覆蓋、少覆蓋、多覆蓋和全覆蓋四類植被覆蓋度區(qū)域內(nèi)劃定樣方(圖3)。獲得樣方內(nèi)對應(yīng)格網(wǎng)點(diǎn)的植被覆蓋率和地表反照率值,剔除異常值(FVC <0,或FVC> 1)后繪制反照率/植被覆蓋梯形特征空間散點(diǎn)圖(圖4)。

圖3 不同植被覆蓋度下樣方分布Fig.3 Sample distribution under different vegetation coverage

圖4 植被覆蓋度/表面反照率梯形特征空間散點(diǎn)圖Fig.4 Trapezoidal feature space scatter plot of vegetation coverage / surface albedo

2.2 土壤濕度估算模型的建立

劉培君等研究表明SM與地表反照率或裸土反照率呈自然對數(shù)關(guān)系[33]。因此,根據(jù)實(shí)測SM值和對應(yīng)樣點(diǎn)的地表反照率或裸土反照率值通過回歸分析建立估算模型和確定精度系數(shù)R2(表1),從R2可以看出,裸土反照率模型的準(zhǔn)確性要明顯高于地表反照率模型。

表1 地表反照率/裸土反照率的土壤水分回歸模型Table 1 Sregression model of surface albedo / bare soil albedo

2.3 模型精度分析

實(shí)測數(shù)據(jù)與SM估算模型之間的誤差分析(表2)表明：除01、02、11和12號點(diǎn)外,所有裸土反照率模型結(jié)果的AE和RE都小于地表反照率模型,RMSE、MAPE、TIC和R2值皆是裸土反照率模型小于地表反照率模型。此外,本文擬合度(R2=0.92)與同類研究(R2=0.57)和同流域研究(R2=0.79)比較精度都有所提高[34- 35]。結(jié)果表明裸土反照率模型能有效地消除植被覆蓋對地表土壤光譜的影響,明顯地提高了SM的估算精度。

表2 模型的精度分析與比較Table 2 Accuracy analysis and comparison of models

2.4 土壤濕度空間分布

參考WU Chengyong在青海對土壤水分的分級標(biāo)準(zhǔn)[36],結(jié)合研究區(qū)SM實(shí)際狀況,對石羊河流域全年平均土壤水分分級(圖5),結(jié)果顯示流域內(nèi)土壤水分整體較低,上游至下游土壤水分依次遞減,且各級區(qū)域依次相鄰分布。其中,SM大于25%的區(qū)域1254.53km2,占比最低為3.09%,主要分布在流域西南部祁連山山麓地帶；SM為15%—25%的區(qū)域面積占比8.68%；這兩級皆分布在海拔較高區(qū)域,中游僅紅崖山水庫區(qū)有小范圍分布；SM為7%—15%的區(qū)域占比最高為33.91%,主要分布于涼州農(nóng)灌區(qū)和民勤綠洲鄰近水庫部分；SM小于7%的區(qū)域則分布在民勤縣荒漠區(qū)。

圖5 土壤水分空間分布Fig.5 Spatial distribution of soil moisture

2.5 土壤水分穩(wěn)定性

用2018年10月—2019年11月的46期SM數(shù)據(jù),依據(jù)式(3)計算SM的穩(wěn)定性(圖6),反映土壤水分的年內(nèi)變化。良好的穩(wěn)定性代表SM幾乎保持不變,受季節(jié)變化和外部影響較??；差的穩(wěn)定性則變化較大,不利于SM的保持。圖6表明：穩(wěn)定性最好的區(qū)域主要分布于下游,因?yàn)樵搮^(qū)域常年裸露且SM值低。穩(wěn)定性較好的區(qū)域主要分布于中游和民勤綠洲,穩(wěn)定性差的區(qū)域主要分布于上游,但占比最小,主要受林地植被覆蓋隨季節(jié)變化和積雪的影響。SM不同水平穩(wěn)定性的面積占比如圖6所示。

圖6 土壤水分穩(wěn)定性特征Fig.6 Characteristics of soil moisture stability

2.6 SM空間自相關(guān)特征

利用ArcGIS和Geoda軟件分析研究區(qū)土壤水分的空間聚集特征,全局Moran′sI系數(shù)為0.88(Z-score=1852.94,P<0.01),表明石羊河流域SM存在強(qiáng)烈的全局空間自相關(guān)性； Moran′sI散點(diǎn)圖如圖7所示,樣點(diǎn)主要分布在第一(高-高,正相關(guān),占總散點(diǎn)41.95%)和第三象限(低-低,正相關(guān),占總散點(diǎn)57.97%),表明研究區(qū)SM值呈現(xiàn)高-高和低-低聚集格局。

圖7 土壤水分全局 Moran′s I 系數(shù)散點(diǎn)圖Fig.7 Scatter plot of soil moisture global Moran′s I coefficient

此外,基于局域Moran′sI指數(shù)分析(P<0.05)得到局部空間自相關(guān)特征(圖8)。SM聚類被分為四種類型：高-高(正空間自相關(guān),占21.95%),主要分布于石羊河上游的肅南、天祝、武威、永昌等地；低-低(正空間自相關(guān),占30.33%),主要集中分布于下游民勤荒漠區(qū)；低-高(負(fù)空間自相關(guān),占0.04%)和高-低(負(fù)空間自相關(guān),占0.01%)區(qū)域占比極低,且未出現(xiàn)聚集特征。中游區(qū)域則未表現(xiàn)出顯著的空間自相關(guān)特征。

圖8 土壤水分聚類分布圖Fig.8 Soil moisture cluster distribution

2.7 空間分異分析

利用地理探測器模型對SM空間分異特征歸因,按照不同因子對SM空間分異特征的解釋能力排序:植被覆蓋度(q=0.42,P<0.001)>土壤類型(q=0.37,P<0.001)>高程(q=0.22,P<0.001)>土地利用(q=0.09,P<0.001)。表明FVC是該地區(qū)SM空間分布的主導(dǎo)因素,土壤類型次之。風(fēng)險探測分析(圖9)和交互探測器分析(表3)表明：不同土地利用對SM的影響存在差異,其中林地SM均值最大,耕地次之,未利用地最小。不同土壤類型對SM的影響也存在較大差異。且SM往往不受單因子作用,而是多因子協(xié)同交互影響。表3得出,交互作用最強(qiáng)的因子組合為土壤類型和植被覆蓋度(q=0.503),說明土壤類型的不同顯著增加了FVC對 SM的解釋力。因子組合皆產(chǎn)生了較高的q值,同時因子交互作用表現(xiàn)為較強(qiáng)的非線性增強(qiáng)交互作用,且兩因子對SM空間分布的影響有顯著差異。

圖9 風(fēng)險探測器結(jié)果Fig.9 Results of risk detector

表3 交互探測器和生態(tài)探測器結(jié)果Table 3 Results of interactive detector and ecological_detector

3 討論

空間統(tǒng)計不同土壤類型的SM特征(表4)表明：石羊河流域內(nèi)半淋溶土面積最小,主要分布在海拔較高的山谷和陰坡,其SM均值最高為24.28%。鈣層土和干旱土依次與半淋溶土鄰近分布,其SM均值也相對較高。漠土和初育土分布廣泛,覆蓋中下游廣大區(qū)域,該土壤類型區(qū)內(nèi)無河流分支,植被稀少,蒸散發(fā)強(qiáng)烈,故SM均值較低。半水成土主要分布在武威和民勤綠洲邊緣區(qū),處在初育土和人為土的過渡帶,鹽堿土零星分布在初育土中,SM值皆較低。而人為土分布在中游,主要由灌溉農(nóng)田組成。高山土分布在流域內(nèi)海拔最高的西北部高山區(qū),該區(qū)域常年氣溫較低,林地植被郁閉度較高,土壤腐殖質(zhì)層厚,蒸發(fā)量小,益于SM的保持。這種現(xiàn)象與其自身的理化性質(zhì)、地形因素、局地氣候和植被郁閉度等有關(guān),這些因素對SM的有效保持和空間分布都產(chǎn)生了重要的影響[37]。

表4 石羊河流域主要土壤類型的土壤水分特征Table 4 Soil moisture characteristics of the main soil types in the Shiyang River Basin

為探討不同土地利用類型的土壤水分特征,依據(jù)土地利用圖(圖10)空間統(tǒng)計各土地利用的SM狀況(表5)。草地大部分處于中等以下水平,其中較干區(qū)域占29.79%,中等濕度區(qū)域占53.48%。耕地SM大部分處于中等水平,其中中等濕度區(qū)域占69.92%。結(jié)合圖6,發(fā)現(xiàn)草地與耕地的SM穩(wěn)定性較好。林地SM整體較高,中等濕度及以上部分占80.97%。未利用地大部分的SM較低,其中干土壤44.7%,較干土壤49.63%。不同的土地利用類型對SM有重要影響,林地在涵養(yǎng)土壤水分方面優(yōu)于耕地和草地[38]。

圖10 土地利用類型Fig.10 Types of land use

表5 不同土地利用的不同土壤水分百分比Table 5 Different soil moisture percentages for different land uses

土壤水分和植被覆蓋度空間相關(guān)性分析(圖11)表明：SM和FVC具有較強(qiáng)的相關(guān)性,植被覆蓋度對土壤水分產(chǎn)生了重要的影響,涼州區(qū)和民勤縣低海拔區(qū)表現(xiàn)為負(fù)向強(qiáng)相關(guān)(R=-0.57,P<0.01),其中草地占66.51%,耕地占25.96%。主要因?yàn)樵搮^(qū)域降水稀少、蒸散發(fā)強(qiáng)烈,淺根系植被水分吸收以地表土壤水為主,不利于土壤水的保持,但耕地長期耕作土壤質(zhì)地松軟且有灌溉水源,因此耕地占比較小。FVC與SM存在負(fù)相關(guān)(R=-0.37,P<0.01)區(qū)域則主要分散在下游,該區(qū)域靠近巴丹吉林沙漠和騰格里沙漠,沙漠化效應(yīng)影響較大。結(jié)論與柳菲等研究結(jié)論基本一致[39]。上游高海拔區(qū)域SM與FVC則為正向強(qiáng)相關(guān)(R=0.49,P<0.01),其中林地占80.16%,海拔較高區(qū)域的FVC與SM為正相關(guān)(R=0.13,P<0.01),其中林地占38.25%,因此林地對SM表現(xiàn)出了正向影響。

圖11 土壤水分與植被覆蓋相關(guān)性Fig.11 Correlation between soil moisture and vegetation cover

李浙華等[40]研究顯示,土壤水分通常受到氣候、水文、地形、植被覆蓋度以及人類活動等因素的影響。但在小尺度上氣象和水文要素基本一致,地形、植被和土地利用的變化則成為影響SM保持和變化的重要因素。植被覆蓋度高的區(qū)域地表蒸散發(fā)較小,有利于SM保持。同時,不同高程也有影響作用,上游高海拔區(qū)域溫度較低且多為林地,使地表蒸散發(fā)更小,且植被多為高大灌木,根系發(fā)達(dá),水分吸收多源于深層土壤。而下游植被則多為草本和農(nóng)田,冠層密度低且根系較淺,主要吸收地表土壤水。從而不同的高程和土地利用影響了植被對土壤水分的影響程度。

限于作者水平,高山地區(qū)冰川積雪對土壤水分遙感估算的影響有待于后續(xù)深入研究,以期進(jìn)一步提高土壤水分的遙感估算精度。

4 結(jié)論

(1)通過地表反照率和植被覆蓋度構(gòu)建梯形特征空間計算裸土反照率,進(jìn)而估算地表土壤水分,可明顯地消除植被覆蓋的影響,提高遙感土壤水分估算精度,為流域尺度的SM計算提供了方法參考。該模型適用于石羊河流域,結(jié)果符合研究區(qū)實(shí)際情況。

(2)石羊河流域土壤水分表現(xiàn)出西南高東北低的態(tài)勢,上游至下游SM值逐漸降低,SM值大于15%的區(qū)域面積為11.77%,主要分布在上游高海拔區(qū)域,SM值7%—15%的區(qū)域面積為33.91%,主要分布在武威綠洲和民勤綠洲鄰近水庫部分,SM值小于7%的區(qū)域面積為54.32%,主要分布在下游民勤境內(nèi)的荒漠區(qū)。

(3)石羊河流域SM空間自相關(guān)特征顯著,高-高聚集主要分布于石羊河上游林地,占研究區(qū)21.95%,低-低聚集主要集中分布于下游民勤的荒漠,占研究區(qū)30.33%。上游林地區(qū)SM呈高高聚集,下游荒漠區(qū)SM呈低低聚集。土壤水分和植被覆蓋度存在顯著相關(guān)性且不同土地利用影響程度存在差異,上游林地與SM呈顯著正相關(guān)(R=0.49,P<0.01),中游草地與SM呈顯著負(fù)相關(guān)(R=-0.57,P<0.01)。

(4)石羊河流域植被覆蓋因素對SM的空間差異解釋力最強(qiáng),各因子對SM空間分異的解釋能力存在顯著差異,且因子間交互作用均增強(qiáng)了其對SM空間分異的解釋力。

(5)石羊河流域土壤水分穩(wěn)定性整體較好,中下游皆表現(xiàn)為穩(wěn)定。不同土地利用類型的SM差異較大,未利用地大部分SM小于7%,大部分草地和耕地SM處于中等水平,SM值為7%—15%,林地SM值大于25%,水平整體較高。