地表反照率不同計(jì)算方法對(duì)干旱區(qū)流域蒸散反演結(jié)果的影響
——以新疆三工河流域?yàn)槔?/h1>

2019-05-31 00:51:32張振宇李小玉

生態(tài)學(xué)報(bào)訂閱 2019年8期 收藏

關(guān)鍵詞：反照率計(jì)算方法通量

張振宇,李小玉,孫 浩

1 浙江農(nóng)林大學(xué)林業(yè)與生物技術(shù)學(xué)院, 杭州 311300 2 中國(guó)科學(xué)院新疆生態(tài)與地理研究所荒漠與綠洲生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 烏魯木齊 830011 3 中國(guó)科學(xué)院大學(xué),北京 100049

地表蒸散(Evapotranspiration,簡(jiǎn)稱ET)是地球表面與大氣相互作用過(guò)程中的關(guān)鍵一環(huán),是地球水量平衡的重要組成部分,地表蒸散的過(guò)程伴隨著能量的吸收與釋放,也是維持地球能量平衡的重要環(huán)節(jié)。近年來(lái),針對(duì)全球氣候變化的研究不斷深入,對(duì)地表蒸散的研究也越來(lái)越多,在這個(gè)過(guò)程中產(chǎn)生了較多的ET估算模型,如SEBAL模型[1]、SEBS模型[2]、ETwatch模型[3]等,其中SEBAL模型是應(yīng)用最為廣泛的ET估算模型之一。SEBAL模型是由Bastiaanssen提出的一種地表蒸散估算模型,研究者利用SEBAL模型結(jié)合不同種類遙感數(shù)據(jù)進(jìn)行了不同區(qū)域的蒸散反演,并且取得了較好的反演結(jié)果。Santos[4]利用Landsat5 TM數(shù)據(jù)結(jié)合SEBAL模型進(jìn)行了巴西布吉達(dá)河流域蒸散反演,Rawat等[5]基于Landsat7 TM數(shù)據(jù)和SEBAL模型對(duì)印度哈里亞邦區(qū)域的小麥農(nóng)田進(jìn)行了蒸散反演,Alemayehu等[6]利用MODIS數(shù)據(jù)結(jié)合SEBAL模型進(jìn)行了東非瑪拉盆地的地表蒸散估算。曾麗紅、杜嘉等[7-9]利用Landsat TM和MODIS數(shù)據(jù)結(jié)合SEBAL模型對(duì)中國(guó)東北地區(qū)的扎龍濕地、三江平原以及松嫩平原地區(qū)進(jìn)行了蒸散量的反演與分析；張殿君等[10]利用SEBAL模型對(duì)中國(guó)黃土高原典型流域區(qū)進(jìn)行了蒸散反演,并討論了土地利用對(duì)蒸散的影響；張發(fā)耀等[11]利用MODIS數(shù)據(jù)結(jié)合SEBAL模型對(duì)浙江省蒸散量進(jìn)行了估算；王秋云[12]利用SEBAL模型進(jìn)行了北京平原造林區(qū)的蒸散研究；周玲等[13]基于Landsat TM/ETM+數(shù)據(jù),應(yīng)用SEBAL模型對(duì)中國(guó)西南地區(qū)漓江流域進(jìn)行了蒸散量變化分析；楊肖麗等[14]利用MODIS數(shù)據(jù)和SEBAL模型對(duì)中國(guó)半干旱區(qū)的老哈河流域進(jìn)行了蒸散研究；李寶富、陳亞寧等[15]利用TM數(shù)據(jù),結(jié)合SEBAL模型對(duì)中國(guó)干旱區(qū)的塔里木河干流區(qū)進(jìn)行了蒸散估算；張楠楠等[16]利用HJ- 1B數(shù)據(jù)結(jié)合SEBAL模型對(duì)淮河上游段進(jìn)行了蒸散估算。以上研究都是利用SEBAL模型實(shí)現(xiàn)了地表蒸散的反演,而對(duì)SEBAL模型的參數(shù)計(jì)算方法影響討論較少,只有少數(shù)研究討論了SEBAL模型反演結(jié)果受氣象條件[17]和地表類型的影響。地表反照率,是指地表物體各個(gè)方向上發(fā)射的太陽(yáng)輻射通量與到達(dá)該物體表面上的總輻射通量之比,是各個(gè)方向上反射率的積分。地表反照率是影響地球能量平衡的關(guān)鍵變量,相關(guān)研究[18-19]表明,地表反照率的增加促使地表凈輻射減少,造成區(qū)域降水和云量的降低,同時(shí),云量的減少使得地表凈輻射通量進(jìn)一步增大,促進(jìn)降水和云的形成,這種反饋?zhàn)饔镁S持著區(qū)域氣候與熱量的穩(wěn)定與平衡。在遙感反演研究中,針對(duì)同一物理量的多種求解方法可能導(dǎo)致衍生結(jié)果的差異。地表反照率作為影響地表能量平衡的關(guān)鍵變量,是SEBAL模型的重要參數(shù),分析不同地表反照率計(jì)算方法造成的SEBAL模型反演結(jié)果差異,對(duì)提高SEBAL模型反演精度具有重要意義。本文擬在控制SEBAL模型其他參數(shù)一致的基礎(chǔ)上,基于兩種地表反照率計(jì)算方法,對(duì)研究區(qū)進(jìn)行蒸散反演,通過(guò)反演值與實(shí)測(cè)值的擬合程度、偏離程度差異,實(shí)現(xiàn)不同地表反照率計(jì)算方法對(duì)SEBAL模型的影響分析。

1 研究區(qū)概況

三工河流域(圖1),位于新疆阜康市境內(nèi),區(qū)域主要由綠洲與部分荒漠構(gòu)成,地勢(shì)南高北低。平原區(qū)北部為古爾班通古特沙漠,南部為天山山脈,境內(nèi)自西向東分布有水磨河、三工河、四工河3條河流,區(qū)域內(nèi)干旱少雨,年降水量150—240 mm,水資源主要依賴高山冰川和積雪融水。研究區(qū)氣候類型為典型的溫帶大陸性氣候,夏季日照時(shí)間長(zhǎng),熱量高,冬季較早,寒冷漫長(zhǎng),流域內(nèi)主要農(nóng)作物為小麥、玉米等經(jīng)濟(jì)作物。

圖1 三工河流域位置示意圖Fig.1 Location of Sangong river basin

2 研究數(shù)據(jù)及預(yù)處理

2.1 遙感和其他空間數(shù)據(jù)獲取與預(yù)處理

①Landsat 8 OLI/TIRS影像數(shù)據(jù)來(lái)源于美國(guó)地質(zhì)調(diào)查局網(wǎng)站(https://glovis.usgs.gov/),共六期,日期分別為2016年5月25日、2016年6月26日、2016年7月12日、2016年7月28日、2016年8月29日和2016年9月30日。影像數(shù)據(jù)整體云量小于0.6%,其中6月26日和7月12日影像云量較高,分別為1.1%和7.45%。對(duì)6期Landsat 8 OLI/TIRS影像數(shù)據(jù)進(jìn)行輻射定標(biāo)、大氣校正獲取糾正后數(shù)據(jù)。

②ASTER高程數(shù)據(jù)來(lái)源于馬里蘭大學(xué)全球土地覆蓋數(shù)據(jù)庫(kù)(http://www.landcover.org/data/),空間分辨率為30 m,對(duì)高程數(shù)據(jù)進(jìn)行拼接與裁剪得到所需研究區(qū)高程數(shù)據(jù)。

③三工河流域矢量邊界由中國(guó)科學(xué)院新疆遙感與地理信息系統(tǒng)重點(diǎn)實(shí)驗(yàn)室提供。

2.2 波文比數(shù)據(jù)

波文比數(shù)據(jù)主要來(lái)自于中國(guó)科學(xué)院新疆遙感與地理信息系統(tǒng)重點(diǎn)實(shí)驗(yàn)室設(shè)立的2臺(tái)波文比觀測(cè)站,經(jīng)緯度坐標(biāo)分別為(44.33°N,87.89°E)和(44.36°N,87.87°E),波文比觀測(cè)站放置于均一的農(nóng)田內(nèi)部,基于波文比-能量平衡法(BREB)[20]獲得瞬時(shí)蒸散,進(jìn)一步獲得每小時(shí)蒸散和日蒸散。大量研究證明[21-22]利用波文比觀測(cè)儀器測(cè)量蒸散精度較高。波文比數(shù)據(jù)監(jiān)測(cè)間隔為1h。

3 模型及參數(shù)

SEBAL模型[23]的理論基礎(chǔ)是地表能量平衡方程,即地表的凈輻射通量由從下墊面到大氣的感熱通量、從下墊面到大氣的潛熱通量和土壤熱通量組成：

Rn=G+H+λET

(1)

式中,Rn代表凈輻射通量(W/m2),G代表土壤熱通量(W/m2),H代表感熱通量(W/m2),λET代表潛熱通量(W/m2),λ代表水的汽化潛熱,通常為2.49×106W m-2mm-1；ET表示蒸散量,kg m-1s-1。

3.1 凈輻射通量(Rn)計(jì)算

凈輻射通量計(jì)算公式為：

Rn=(1-a)Rsd+Rld-Rlu-(1-k)Rld

(2)

式中,a為地表反照率,Rsd為太陽(yáng)下行短波輻射(W/m2),Rld為太陽(yáng)下行長(zhǎng)波輻射(W/m2),Rlu為地表反射的長(zhǎng)波輻射(W/m2),k為地表比輻射率(發(fā)射率)。

3.2 土壤熱通量(G)的計(jì)算

土壤熱通量是單位時(shí)間、單位面積上的土壤熱交換量,一般存儲(chǔ)在土壤和植被中,這部分能量在整個(gè)熱量平衡過(guò)程中占比重較小,而且計(jì)算較為困難,這里使用曾麗紅[24]推導(dǎo)的公式。

(3)

式中,Ts為地表溫度,NDVI為歸一化植被指數(shù),c表示衛(wèi)星過(guò)境時(shí)間對(duì)G的影響,過(guò)境時(shí)間在地方時(shí)12：00點(diǎn)以前取0.9；在12：00點(diǎn)到14：00點(diǎn)之間取1.0；在14：00點(diǎn)到16：00點(diǎn)之間取1.1。

3.3 感熱通量(H)的計(jì)算

感熱通量表征的是下墊面與大氣間的湍流形式的熱量交換,與大氣穩(wěn)定度、風(fēng)速、表面粗糙度有關(guān)。通用計(jì)算公式如下：

(4)

式中,Pair為空氣密度(kg/m3)；Cp是空氣定壓比熱(1004 J kg-1K-1)；dT是距離地面高Z1和Z2處的溫差(一般Z1為0.1 m,Z2為2 m),rah為空氣動(dòng)力學(xué)阻抗。

3.4 小時(shí)蒸散量與日蒸散量計(jì)算

小時(shí)蒸散量的計(jì)算公式[25]為：

(5)

λ=2.501-0.002361×(Ts-273.15)

(6)

式中,λ為水的汽化潛熱,ET為蒸散量,ETinst為小時(shí)蒸散量。

通過(guò)小時(shí)蒸散量結(jié)合日照時(shí)長(zhǎng)計(jì)算日蒸散量ETday：

(7)

式中,N為區(qū)域內(nèi)平均日照時(shí)長(zhǎng),t為日出到衛(wèi)星過(guò)境時(shí)間間隔(h)。

4 地表反照率計(jì)算

4.1 基于各波段反射率的大氣表觀反射率計(jì)算

求取大氣表觀反射率,這一過(guò)程在輻射定標(biāo)中即可完成：

(8)

式中,ρ代表表觀反射率,DN為影像灰度值,μ、φ分別為landsat8 OLI/TIRS影像的反射率增量和修正值,θ為太陽(yáng)高度角。

4.2 求取地表反照率

①Liang地表反照率計(jì)算方法(下稱Liang方法)[26]:

(9)

式中,ρ2、ρ4、ρ5、ρ6、ρ7為各波段大氣表觀反射率,ap為大氣程輻射值影響,一般取0.03,tsw為大氣單向透過(guò)率。

大氣單向透過(guò)率計(jì)算公式如下[27]：

tsw=0.75+2×10-5×Z

(10)

式中,Z為海拔高度(m),從DEM數(shù)據(jù)中獲取。

②Smith地表反照率計(jì)算方法(下稱Smith方法)[28]:

(11)

式中,ρi為波段反射率(對(duì)于Landsat 8 OLI/TIRS數(shù)據(jù)使用2、4、5、6、7波段),Δλi為各波段范圍,Si為1—7波段的平均太陽(yáng)輻照度(表1),ap為大氣程輻射值影響,一般取0.03,tsw為大氣單向透過(guò)率。

表1 Landsat 8 OLI數(shù)據(jù)波段平均太陽(yáng)輻照度[28]

Si,平均太陽(yáng)輻照度,Average solar irradiance

5 結(jié)果與分析

5.1 日蒸散量比較

基于SEBAL模型得到兩種地表反照率計(jì)算方法下反演的地表日均蒸散量(表2)及其空間分布圖(圖2和圖3),對(duì)比6個(gè)研究日區(qū)域日蒸散量的分布,結(jié)果顯示兩種地表反照率方法下最終得到的日蒸散量空間和時(shí)間分布規(guī)律保持一致。時(shí)間上,日蒸散量在6月—7月末明顯高于其他月份；空間上,流域綠洲的日蒸散量明顯高于其他類型區(qū)域(水體除外)；另一方面,兩種地表反照率方法下基于SEBAL模型得到的日蒸散量分布曲線(圖4)則有較大不同。5月25日,兩種地表反照率方法下的日蒸散量分布曲線在形態(tài)上保持一致,但是Liang方法下的日蒸散量分布曲線在6 mm/d處開(kāi)始接近橫坐標(biāo)軸,而Smith方法下日蒸散量分布曲線在日蒸散量7 mm/d處才開(kāi)始接近橫坐標(biāo)軸。6月26日、7月12日、7月28日三天同一地表反照率計(jì)算方法下的日蒸散量分布曲線形態(tài)上較為一致,但是不同地表反照率計(jì)算方法下得到的曲線則差異較大：Liang方法下得到的3期日蒸散量分布曲線波峰出現(xiàn)在日蒸散量4 mm/d附近,在日蒸散量6 mm/d處開(kāi)始接近橫坐標(biāo)軸,而Smith方法下的3期日蒸散量分布曲線均在日蒸散量6—6.5 mm/d處達(dá)到峰值,在日蒸散量7.5 mm/d處才開(kāi)始接近橫坐標(biāo)軸。8月29日和9月30日兩期日蒸散量分布曲線則呈現(xiàn)出與前4期較為不同的特征：曲線波峰開(kāi)始向低日蒸散量區(qū)間移動(dòng),Liang方法下得到的日蒸散量分布曲線波峰分布在日蒸散量1—1.5 mm/d附近,而Smith方法下得到的日蒸散量分布曲線在日蒸散量1.5—2.5 mm/d處達(dá)到峰值。兩種方法下最終得到的多條日蒸散量分布曲線在日蒸散量區(qū)間末端均出現(xiàn)了一個(gè)較為明顯的波峰,是由區(qū)域內(nèi)的水庫(kù)水體所引起。總的來(lái)說(shuō),Smith方法下基于SEBAL模型反演的區(qū)域日蒸散量分布曲線峰值出現(xiàn)更晚。

此外,兩種地表反照率計(jì)算方法下最終求得的日均蒸散量也有較大差異,Liang方法下最終得到的區(qū)域日均蒸散量在所選取的6個(gè)研究日內(nèi)均小于Smith方法下最終得到的結(jié)果,其中,春季(5月25日)、秋季(8月29日和9月30日)的日均蒸散量差距較小,平均差異0.2 mm/d,而夏季(6月26日、7月12日、7月28日)的日均蒸散量差距較大,平均差異0.64 mm/d。

表2 兩種地表反照率計(jì)算方法下得到的區(qū)域日均蒸散量/(mm/d)

圖2 Liang地表反照率計(jì)算方法下日蒸散量分布Fig.2 Daily evapotranspiration distribution using Liang surface albedo formula

圖3 Smith地表反照率計(jì)算方法下日蒸散量分布Fig.3 Daily evapotranspiration distribution using Smith surface albedo formula

圖4 Liang反照率計(jì)算方法和Smith反照率計(jì)算方法下的日蒸散量分布曲線Fig.4 Distribution curve of daily evapotranspiration using Liang and Smith surface albedo method

5.2 兩種地表反照率計(jì)算方法下的日蒸散反演精度比較

將SEBAL模型反演值與實(shí)測(cè)值進(jìn)行擬合,并使用決定系數(shù)R2、均方根誤差RMSE、平均偏差BIAS和平均絕對(duì)偏差MAE[29]對(duì)兩種地表反照率方法下基于SEBAL模型反演的日蒸散量進(jìn)行比較與評(píng)價(jià)。評(píng)價(jià)指標(biāo)計(jì)算公式如下:

(10)

(11)

(12)

式中,N為觀測(cè)個(gè)數(shù),Pi為反演值,Qi為實(shí)測(cè)值。得到的結(jié)果如表3所示。

表3 兩種地表反照率計(jì)算方法下SEBAL模型反演值精度比較

R2,決定系數(shù),Coefficient of determination；RMSE,均方根誤差,Root mean square error；BIAS,平均偏差,Bias；MAE,平均絕對(duì)偏差,Mean absolute error

兩種地表反照率計(jì)算方法下基于SEBAL模型反演的日蒸散量與實(shí)測(cè)值散點(diǎn)圖及擬合方程(圖5)表明,Smith方法下最終得到的日蒸散量與實(shí)測(cè)值的線性擬合程度更高,決定系數(shù)達(dá)到了0.85,相比之下,Liang方法下得到的模型反演日蒸散量與實(shí)測(cè)值的線性擬合決定系數(shù)僅為0.76；F檢驗(yàn)顯示,Smith方法下最終得到反演值與實(shí)測(cè)值線性擬合P值<0.001,而Liang反照率計(jì)算方法下得到的擬合直線P值為0.0016,大于0.001,可見(jiàn)相同時(shí)間空間條件下,Smith方法下基于SEBAL模型得到的日蒸散量與實(shí)測(cè)值更為接近,可信程度更高。另外通過(guò)RMSE、BIAS、MAE 3種指標(biāo)對(duì)反演值與實(shí)測(cè)值進(jìn)行偏離程度分析,Smith方法下最終得到的日蒸散量與實(shí)測(cè)值的偏離程度小于使用Liang方法下得到的結(jié)果,多項(xiàng)結(jié)果表明：在應(yīng)用SEBAL模型對(duì)干旱區(qū)地表進(jìn)行蒸散反演過(guò)程中,使用Smith方法結(jié)果作為SEBAL模型地表反照率參數(shù)輸入更優(yōu)。

圖5 Liang反照率計(jì)算方法和Smith反照率計(jì)算方法下的日蒸散量擬合Fig.5 Daily evapotranspiration fitting using Liang and Smith surface albedo method

6 討論

6.1 兩種方法計(jì)算下的地表反照率比較

2016年6個(gè)研究日的流域地表反照率分布曲線(圖6)表明,兩種方法得到的地表反照率分布曲線類似,都呈較明顯的正態(tài)分布,然而不同方法下的曲線波峰處反照率存在明顯差異,基于Liang方法得到的地表反照率分布曲線波峰處反照率在6個(gè)時(shí)期均高于Smith方法得到的曲線波峰反照率,波峰處地表反照率相差0.02—0.06,盡管這種差異較小,但是這種差異使得整個(gè)區(qū)域的反照率分布出現(xiàn)較大不同,Smith方法下得到的地表反照率更低,兩種方法計(jì)算出的地表反照率均值(表4)可進(jìn)一步印證這一點(diǎn),地表反照率的差異使得由此計(jì)算出的流域凈輻射通量、顯熱通量產(chǎn)生了更大的差異,進(jìn)一步的造成了兩種地表反照率計(jì)算方法下基于SEBAL模型反演結(jié)果在分布上的差異。

圖6 兩種方法計(jì)算下反照率像元分布曲線Fig.6 Albedo pixel distribution curves under two methods

地表反照率計(jì)算方法Surface albedo calculation method5月25日May 25th6月26日J(rèn)une 26th7月12日J(rèn)uly 12th7月28日J(rèn)uly 28th8月29日August 29th9月30日September 30thLiang0.36380.34690.32850.31830.30130.2916Smith0.29440.26840.25520.25580.25200.2520

6.2 兩種方法計(jì)算下地表反照率與MCD43A3產(chǎn)品的比較

圖7 地表反照率檢測(cè)點(diǎn)位置Fig.7 Location of surface albedo detection points

為進(jìn)一步比較兩種方法計(jì)算出的地表反照率在研究區(qū)的可靠性,在缺乏地表反照率測(cè)量?jī)x器的情況下,參考He等[30]的方法,使用MODIS地表反照率產(chǎn)品MCD43A3與兩種方法計(jì)算出的地表反照率進(jìn)行比較,由于研究區(qū)所選Landsat 8數(shù)據(jù)云量均較少(晴空條件),因此提取MCD43A3數(shù)據(jù)中的短波白空反照率(WSA)作為參考標(biāo)準(zhǔn)。在研究區(qū)內(nèi)部隨機(jī)布置了36個(gè)監(jiān)測(cè)點(diǎn)(圖7),通過(guò)6個(gè)研究日共得到216個(gè)數(shù)據(jù)。數(shù)據(jù)比對(duì)結(jié)果(圖8)表明：①兩種方法下得到的地表反照率均偏高于MCD43A3 WSA,Liang方法計(jì)算出的地表反照率與MCD43A3 WSA擬合程度(R2=0.35)略高于Smith方法計(jì)算結(jié)果(R2=0.2),說(shuō)明Liang方法計(jì)算出的地表反照率一致性更好。②RMSE比較上,Smith方法計(jì)算出的地表反照率與MCD43A3 WSA誤差相對(duì)較小(Smith方法RMSE=0.075；Liang方法RMSE=0.13)。

圖8 兩種方法計(jì)算下地表反照率與MCD43A3短波白空反照率對(duì)比Fig.8 Comparison between surface albedo calculated by two methods and MCD43A3 shortwave white sky albedo

同時(shí)由于Landsat 8數(shù)據(jù)部分參數(shù)官方尚未明確公布,計(jì)算Landsat 8數(shù)據(jù)地表反照率時(shí)仍采用針對(duì)上一代傳感器的方法,這也可能會(huì)對(duì)結(jié)果造成一定影響；基于此,Smith提出了針對(duì)Landsat 8數(shù)據(jù)地表反照率的計(jì)算方法,并在阿拉斯加、佛羅里達(dá)等地區(qū)進(jìn)行了實(shí)驗(yàn),近期部分學(xué)者[30]提出利用大氣頂反照率和站點(diǎn)測(cè)量數(shù)據(jù)直接估計(jì)Landsat 8地表反照率,實(shí)驗(yàn)效果較好,但由于中國(guó)西北地區(qū)缺少地表反照率國(guó)際測(cè)量站點(diǎn),因此,各種方法在中國(guó)西北地區(qū)的適用性仍有待進(jìn)一步檢驗(yàn)。

另一方面,由于USGS(美國(guó)地質(zhì)調(diào)查局)至今未公布Landsat 8數(shù)據(jù)各波段平均太陽(yáng)輻照度,因此在使用Smith方法過(guò)程中使用了其提供的各波段平均太陽(yáng)輻照度。在理想狀態(tài)下,平均太陽(yáng)輻照度可通過(guò)大氣光譜曲線和各波段響應(yīng)函數(shù)獲得,分別使用WRC太陽(yáng)光譜曲線和Thuillier 太陽(yáng)光譜曲線結(jié)合Landsat 8數(shù)據(jù)1—7波段響應(yīng)函數(shù)進(jìn)行平均太陽(yáng)輻照度估計(jì)[31-32],結(jié)果(表5)顯示,兩種譜線對(duì)應(yīng)的波段平均太陽(yáng)輻照度存在差異,其中1和4波段差異較大,其他波段結(jié)果相近,同時(shí)Smith方法的1、2、3波段平均太陽(yáng)輻照度均小于兩種譜線對(duì)應(yīng)估計(jì)值,而在第6波段(短波紅外)上高于譜線估計(jì)值,這使得Smith方法在計(jì)算地表反照率時(shí)更側(cè)重于短波紅外波段,而由此得出的地表反照率與相關(guān)產(chǎn)品的偏差較小,這說(shuō)明Landsat 8地表反照率計(jì)算中短波紅外波段的權(quán)重增加可能改善了研究區(qū)地表反照率反演質(zhì)量,部分研究[33]也表達(dá)了相似的結(jié)果。

表5 利用Thuillier和WRC太陽(yáng)光譜曲線估算的Landsat 8 OLI數(shù)據(jù)各波段平均太陽(yáng)輻照度

ThuillierSi,利用Thuillier太陽(yáng)光譜曲線估算的波段平均太陽(yáng)輻照度,Band average solar irradiance based on Thuillier Solar spectral curve；WRCSi,利用WRC太陽(yáng)光譜曲線估算的波段平均太陽(yáng)輻照度,Band average solar irradiance based on WRC Solar spectral curve

7 結(jié)論

研究基于兩種地表反照率計(jì)算方法結(jié)合SEBAL模型,對(duì)西北干旱區(qū)典型流域日蒸散進(jìn)行反演并進(jìn)行比較,結(jié)果表明：(1)Simth方法下最終得到的日均蒸散量均高于Liang方法下最終得到的日均蒸散量,這種差異在夏季最明顯。同時(shí),Smith方法下基于SEBAL模型得到的反演結(jié)果與實(shí)測(cè)值的擬合程度和偏離程度略優(yōu)于Liang方法下得到的結(jié)果。(2)兩種方法計(jì)算的地表反照率差異是導(dǎo)致SEBAL模型反演日蒸散量差異的主要原因,相同日期內(nèi)Liang方法得到的地表反照率均高于Smtih方法結(jié)果(高約0.04—0.08)。(3)與MCD43A3產(chǎn)品的比較結(jié)果顯示,Liang方法得到的地表反照率與MCD43A3產(chǎn)品一致性較好,但是其偏離程度也較高。綜合考慮之下,Smith方法在研究區(qū)的適用性略好。

猜你喜歡

反照率計(jì)算方法通量

基于藍(lán)天空反照率的氣溶膠輻射強(qiáng)迫模擬

中國(guó)環(huán)境科學(xué)(2023年9期)2023-09-25 01:36:46

浮力計(jì)算方法匯集

中學(xué)生數(shù)理化·八年級(jí)物理人教版(2022年5期)2022-06-05 06:57:38

冬小麥田N2O通量研究

農(nóng)業(yè)災(zāi)害研究(2022年1期)2022-05-07 01:31:04

薩吾爾山木斯島冰川反照率時(shí)空變化特征研究

冰川凍土(2022年6期)2022-02-12 08:31:06

長(zhǎng)江三角洲地區(qū)大氣氣溶膠柱單次散射反照率特性研究

南京信息工程大學(xué)學(xué)報(bào)(2021年2期)2021-05-22 01:22:56

隨機(jī)振動(dòng)試驗(yàn)包絡(luò)計(jì)算方法

現(xiàn)代工業(yè)經(jīng)濟(jì)和信息化(2016年2期)2016-05-17 05:34:16

不同應(yīng)變率比值計(jì)算方法在甲狀腺惡性腫瘤診斷中的應(yīng)用

中國(guó)醫(yī)學(xué)影像學(xué)雜志(2015年9期)2015-12-15 11:03:26

緩釋型固體二氧化氯的制備及其釋放通量的影響因素

化工進(jìn)展(2015年6期)2015-11-13 00:26:29

一種伺服機(jī)構(gòu)剛度計(jì)算方法

導(dǎo)航定位與授時(shí)(2014年2期)2014-04-27 13:41:09

春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量

中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版)(2014年12期)2014-02-28 12:22:06

生態(tài)學(xué)報(bào)2019年8期

生態(tài)學(xué)報(bào)的其它文章

中國(guó)城市空間形態(tài)的環(huán)境效應(yīng)

近30年淀山湖地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值對(duì)土地利用變化的響應(yīng)

杭州市西湖景區(qū)春季空氣細(xì)菌多樣性特征

不同抗旱性冬小麥根系時(shí)空分布與產(chǎn)量的關(guān)系

低溫脅迫下6種木蘭科植物的生理響應(yīng)及抗寒相關(guān)基因差異表達(dá)

外源氮持續(xù)輸入對(duì)相應(yīng)氮梯度下堿蓬殘?bào)w分解及硫養(yǎng)分釋放的影響