利用“模擬—校正”法估算黃土高原半干旱區(qū)陸面蒸散發(fā)

楊 揚(yáng), 孫旭映, 張 良, 王麗娟, 胡 蝶

(中國氣象局 蘭州干旱氣象研究所 甘肅省(中國氣象局)干旱氣候變化與減災(zāi)重點(diǎn)實(shí)驗(yàn)室, 蘭州 730020)

陸地蒸散發(fā)是地表水平衡和能量平衡過程中的組成部分，對(duì)區(qū)域乃至全球的水循環(huán)和氣候均有重要的影響。一方面，蒸散發(fā)過程直接影響地表徑流、土壤濕度等水文變量；另一方面，潛熱通量是地表能量平衡的重要組成部分。蒸散發(fā)過程在發(fā)生相變時(shí)吸收和釋放能量會(huì)對(duì)周邊大氣運(yùn)動(dòng)產(chǎn)生顯著影響(液態(tài)水轉(zhuǎn)化為氣態(tài)水所需的能量是溫度升高1 K時(shí)所需能量的600倍，更是相同質(zhì)量空氣溫度升高1 K時(shí)所需能量的2 400倍[1])。研究表明，陸氣間的蒸散發(fā)過程消耗超過50%的地表凈輻射[2]，這些能量將約60%的地表降水通過蒸散發(fā)返回到大氣中[3]。不僅如此，蒸散發(fā)的變化趨勢(shì)在干旱監(jiān)測(cè)預(yù)警方面扮演著至關(guān)重要的角色[4]。例如，Wang等[5]發(fā)現(xiàn)驟發(fā)性干旱主要是由蒸散量劇增造成的。因而，準(zhǔn)確估算蒸散發(fā)對(duì)水資源管理和農(nóng)田灌溉具有重要意義[6]，尤其在水資源極其短缺的干旱和半干旱區(qū)，準(zhǔn)確估算實(shí)際蒸散發(fā)尤為重要。

目前估算實(shí)際蒸散發(fā)的方法主要有3類：第1類是經(jīng)驗(yàn)?zāi)Ｐ?，通常利用土壤水分、葉面積指數(shù)，輻射等估算蒸散發(fā)，這類模型[7-10]通常不考慮蒸散發(fā)的具體物理過程，僅僅利用相關(guān)關(guān)系，具有較大的經(jīng)驗(yàn)性，難以推廣到不同地區(qū)[11]；第2類是基于能量平衡的陸面蒸散發(fā)模型，如結(jié)合遙感和地表面觀測(cè)估算蒸散發(fā)的SEBS，SEBAL模型[12-13]，通過這類模型可以獲得像元尺度的蒸散發(fā)，但受限于衛(wèi)星遙感產(chǎn)品的反演精度，存在的問題是估算誤差較大；第3類詳細(xì)考慮了蒸散發(fā)實(shí)際物理過程的陸面過程模式，如CLM,NOAH等，但是這類模型通常在單點(diǎn)運(yùn)行模擬，難以估算區(qū)域或流域尺度的蒸散發(fā)。楊澤粟[14]利用實(shí)測(cè)資料分別驗(yàn)證了目前國際上主要的8個(gè)蒸散發(fā)估算模型，包括Wang經(jīng)驗(yàn)?zāi)Ｐ汀⑵搅鳌珊礎(chǔ)A模型、能量平衡SEBS模型、改進(jìn)的Priestley-Taylor(PT-JPL)模型、改進(jìn)的Penman-Monteith(RS-PM)模型、PT-JPL改進(jìn)模型(M1-PT-JPL，M2-PT-JPL模型)和陸面過程模式(CLM模式)，檢驗(yàn)了其在黃土高原的適用性。分析得出，各蒸散發(fā)估算模型模擬結(jié)果與觀測(cè)值均存在較大差異，而且不同模型也存在顯著差異，相比而言，PT-JPL，CLM模式估算誤差相對(duì)較小。因此，選用合理的蒸散發(fā)估算模型，提高蒸散發(fā)估算精度是一項(xiàng)重要的研究任務(wù)。

從陸氣相互作用角度而言，陸面過程模式估算蒸散發(fā)具有一定的優(yōu)勢(shì)，主要是陸面過程模式基于能量和水平衡建立，并詳細(xì)考慮了陸地蒸散發(fā)中植被截流降水蒸發(fā)、蒸騰和地表蒸發(fā)等各個(gè)重要的物理過程，尤其是第三代陸面模式(如CLM等)，包含了植被氣孔與光合作用之間的聯(lián)系等[15]，更加真實(shí)地刻畫了蒸散發(fā)過程，成為估算蒸散發(fā)重要的工具[16-19]。目前諸多計(jì)劃如全球土壤濕度計(jì)劃(GSWP)[20]、全球陸面數(shù)據(jù)同化系統(tǒng)(GLDAS)[16]、中國陸面數(shù)據(jù)同化系統(tǒng)(CLDAS)均利用不同的陸面模式估算了區(qū)域或全球陸面蒸散量，廣泛應(yīng)用于全球氣候和水文循環(huán)研究。但不可否認(rèn)的是，由于陸面過程模式中的參數(shù)化方案、參數(shù)選取等均存在較大的不確定性[21]，因而其模擬的蒸散發(fā)與觀測(cè)之間也存在較大的偏差。Haddeland等[17]比較了11個(gè)模式模擬的陸面蒸散發(fā)發(fā)現(xiàn)全球陸面蒸散量在415～586 mm/a之間變化。Xia等[18-19]表明NLDAS-2在區(qū)域和流域尺度上，SAC-SMA模式高估了蒸散量，而Noah和VIC模式則低估了蒸散量。崔園園等[22]基于第3次青藏高原大氣科學(xué)試驗(yàn)資料評(píng)估了CLDAS-V2.0和GLDAS-NOAH產(chǎn)品在青藏高原地區(qū)的適用性表明，CLDAS-V2.0土壤溫、濕度產(chǎn)品與觀測(cè)資料的相關(guān)性均優(yōu)于GLDAS-NOAH產(chǎn)品，他們之間的相關(guān)性在濕季大于干季。張強(qiáng)等[23]分析了FLUXNET和GLDAS的蒸散發(fā)資料在我國北方的適用性。結(jié)果表明，F(xiàn)LUXNET資料在我國北方地區(qū)的年平均蒸散量明顯比GLDAS資料的年平均蒸散量可靠。為此，在區(qū)域尺度上評(píng)估陸面過程模式模擬的蒸散發(fā)產(chǎn)品，并有效提高陸面產(chǎn)品的應(yīng)用，尤其是干旱監(jiān)測(cè)等方面的應(yīng)用值得深入研究。

為了提高陸面過程模式模擬的蒸散發(fā)的精度，一個(gè)重要的途徑是根據(jù)不同下墊面的重要物理過程，改進(jìn)參數(shù)化方案和優(yōu)化關(guān)鍵參數(shù)，如Sellers等[24]在模式中引入氣孔阻抗，更加準(zhǔn)確地描述了蒸散發(fā)過程，進(jìn)而提高了濕潤區(qū)的蒸散發(fā)模擬；Chen等[25]通過改進(jìn)熱力學(xué)粗糙度，改進(jìn)了對(duì)半干旱區(qū)NOAH模式模擬性能；孫菽芬等[26]發(fā)展了通用土壤水熱模型，改進(jìn)了土壤濕度模擬等；Yang等[27]利用粒子群優(yōu)化算法，改進(jìn)了半干旱區(qū)土壤濕度和蒸散發(fā)的模擬性能。近年來，Parr等[28]提出了提高蒸散發(fā)的另一種途徑，可稱之為“模擬—校正”方法，該方法利用某區(qū)域某一時(shí)段內(nèi)模式模擬的蒸散發(fā)與觀測(cè)的蒸散發(fā)建立函數(shù)關(guān)系，并假定這一關(guān)系不隨時(shí)間發(fā)生變化，并適用于該區(qū)域，因而可以用于校正模式模擬的蒸散發(fā)。Parr[28]、Wang[29]等通過研究發(fā)現(xiàn)，利用該方法校正模式的模擬的蒸散發(fā)后，誤差顯著減小。但是，該方法也存在一定的局限，因?yàn)樵趯?shí)際情形下氣候、下墊面狀況和時(shí)空尺度均會(huì)影響蒸散發(fā)的估算精度[7,30]。不僅如此，由于該方法連接了模擬和觀測(cè)，因而也與模式參數(shù)化方案密切相關(guān)。因此，如何在該方法中，引入影響區(qū)域蒸散發(fā)的關(guān)鍵因子，進(jìn)而提高陸面過程模式對(duì)蒸散發(fā)估算能力有重要的研究價(jià)值。

大量研究表明，在半干旱區(qū)，植被覆蓋度和氣孔阻抗等是影響蒸散發(fā)關(guān)鍵因子，而在陸面模式中，這些因子對(duì)模擬結(jié)果也是敏感的。例如，植被覆蓋度決定了植被蒸騰和土壤蒸發(fā)的比例，因此對(duì)蒸散發(fā)的調(diào)控發(fā)揮著重要作用[31-32]，而氣孔導(dǎo)度與光合作用速率和土壤含水量緊密相連[33]，直接影響植被生長，而植被生長狀況的差異能夠良好地反映土壤濕度的不均勻性，植被越茂盛，其蒸騰作用越強(qiáng)，蒸散就越大[34]；在陸面模式中，植被覆蓋度和氣孔阻抗均與植被葉面積密切相關(guān)，為此將葉面積指數(shù)作為一個(gè)關(guān)鍵因子，可用于改善蒸散發(fā)估算。由于在實(shí)際應(yīng)用中，難以精確測(cè)量葉面積指數(shù)，因而通常用使用NDVI代替。為此，本文利用黃土高原半干旱區(qū)SACOL站的觀測(cè)資料，以NDVI為關(guān)鍵因子，結(jié)合CLM模式模擬結(jié)果，建立模擬、觀測(cè)與NDVI之間的函數(shù)關(guān)系，改進(jìn)“模擬—校正”方法，進(jìn)而提高蒸散發(fā)的估算精度，為該地區(qū)陸氣相互作用和區(qū)域氣候研究提供科學(xué)依據(jù)。

1 資料與方法

1.1 資 料

蘭州大學(xué)半干旱氣候與環(huán)境監(jiān)測(cè)站(Semi-Arid Climate Observatory and Laboratory of Lanzhou University,SACOL)位于中國黃土高原地區(qū)，地處溫帶半干旱氣候區(qū)，年降水量為381.8 mm，年平均氣溫6.7℃，代表了方圓幾百千米的半干旱區(qū)氣候狀況，儀器的描述及安裝情況可參考文獻(xiàn)[35—36]。本文使用SACOL站2006—2008年近地層觀測(cè)資料，包括空氣溫度、降水、風(fēng)速、相對(duì)濕度、大氣壓強(qiáng)；向下、向上的短波和長波輻射；感熱和潛熱通量，以及土壤溫度和含水量。本文使用的蒸散發(fā)觀測(cè)數(shù)據(jù)來自于渦動(dòng)相關(guān)系統(tǒng)測(cè)量，為保證觀測(cè)數(shù)據(jù)的可靠性，對(duì)30 min的原始資料進(jìn)行二次坐標(biāo)旋轉(zhuǎn)、超聲虛溫的側(cè)向風(fēng)校正和WPL修正。然后剔除由于降水、露水等天氣導(dǎo)致渦動(dòng)相關(guān)系統(tǒng)觀測(cè)的異常。最后對(duì)儀器故障、系統(tǒng)校正及降水等導(dǎo)致的資料丟失及異常值進(jìn)行插補(bǔ)。插補(bǔ)方法根據(jù)Falge等[37]的方法，對(duì)≤6 h的缺失數(shù)據(jù)進(jìn)行線性內(nèi)插；對(duì)于≥6 h的缺失數(shù)據(jù)用查表法插補(bǔ)。

1.2 模式簡(jiǎn)介

本文使用公用陸面過程模式CLM (Community Land Surface Model,Version 4.0)來計(jì)算蒸散發(fā)。CLM模式所需的地表參數(shù)主要包括植被功能類型，莖、葉面積指數(shù)，植株高度，土壤顏色、質(zhì)地、有機(jī)質(zhì)含量等，與植被功能類型有關(guān)的數(shù)據(jù)來自于MODIS的數(shù)據(jù)，而土壤質(zhì)地來自于IGBP數(shù)據(jù)[33]，土壤有機(jī)質(zhì)含量來自于Lawrence等[38]。CLM模式地表溫度及各層土壤溫度的初始值取為274 K，土壤含水量初始化為0.3，而最下面的5層初始化為0.0。

1.3 模擬-校正方法

Parr等[28]指出利用某區(qū)域某一時(shí)段內(nèi)模式模擬的蒸散發(fā)與觀測(cè)的蒸散發(fā)建立函數(shù)關(guān)系，并假定這一關(guān)系不隨時(shí)間發(fā)生變化，并適用于該區(qū)域，因而可以用于校正模式模擬的蒸散發(fā)，如下式所示：

lEest=k×lEm

(1)

式中：lEest是修訂后的通量；lEm是模式模擬值；k為某段時(shí)間內(nèi)觀測(cè)和模擬值的線性擬合系數(shù)。

本研究中，考慮到NDVI對(duì)蒸散發(fā)的重要影響，改進(jìn)了上述方法，建立了NDVI，觀測(cè)值和CLM模擬值與修正通量之間的如下關(guān)系：

lEest=f(NDVI,lEm)

(2)

利用上述的資料和方法，本文首先利用2006年9—12月的資料進(jìn)行了spin-up，然后使用2007—2012年的模擬結(jié)果進(jìn)行對(duì)比分析。

為評(píng)估模式的模擬性能，給出模擬值與觀測(cè)值的平均偏差(Mean Bias Error,MBE)、均方根偏差(Root Mean Bias Error,RMBE)和相關(guān)系數(shù)(Correlation Coefficient,R)定量評(píng)估模擬值與觀測(cè)值的偏差。

2 結(jié)果與分析

2.1 CLM模擬的蒸散發(fā)分析

圖1給出了CLM模擬的2007—2012年SACOL站蒸散發(fā)(潛熱通量)與實(shí)測(cè)值日均值的比較?？梢?，CLM能較好地模擬出蒸散發(fā)的變化趨勢(shì)，但模擬值較觀測(cè)值偏大，在非生長季，模式的模擬值較觀測(cè)值的偏差較大；在生長季，模擬的蒸散發(fā)與觀測(cè)值的偏差相對(duì)較小。表1為SACOL站2007—2012年地表能量及環(huán)境參數(shù)的變化。6 a的平均降水量(371.8 mm)略低于常年(1961—2008年)的年均降水量(386.5 mm)，但SACOL站的年降水量有很大的年際變率，2008年的年降水量為386.9 mm，與常年平均降水量相差不多，認(rèn)為是正常年；2007年年降水量達(dá)到了555.5 mm，比常年平均降水量偏高44%，2012年的偏高2%，屬于濕潤年；而2009—2011年的年均降水量比常年平均降水量低22%，14%，32%，屬于干旱年。從圖1還可以看出，CLM對(duì)濕潤年蒸散發(fā)的模擬偏差較小。從偏差的統(tǒng)計(jì)結(jié)果來看，濕潤年兩者的RMBE(27.95 W/m2)小于干旱年的36.81 W/m2，而濕潤年的相關(guān)系數(shù)(0.68)大于干旱年的0.59，表明模式在濕潤年對(duì)蒸散發(fā)的模擬偏差較小，且相關(guān)性較高。年降水量決定著蒸散發(fā)能量和水分的可利用性。降水作為CLM的大氣強(qiáng)迫項(xiàng)，通過改變土壤含水量和能量分配影響陸面蒸散發(fā)。在年降水量較多的年份，CLM更準(zhǔn)確地捕捉了土壤含水量的變化趨勢(shì)，因此其對(duì)蒸散發(fā)的模擬效果較好。

圖1 2007-2012年CLM模式模擬的SACOL站潛熱通量與觀測(cè)值的日均值比較

表1 2007-2012年SACOL站地表能量及環(huán)境參數(shù)的變化

注：Rn為凈輻射，H為感熱通量，lE為潛熱通量，Ta為日平均溫度，P為年降水量。

為了進(jìn)一步分析CLM對(duì)干濕不同條件下蒸散發(fā)的模擬性能，本文選取2008年生長季(4—10月)的數(shù)據(jù)進(jìn)行研究。圖2給出了2008年SACOL站日平均土壤含水量和降水的季節(jié)變化。如圖2所示，土壤含水量的變化主要受控于降水。2008年生長季降水時(shí)段主要集中在8月上旬—10月中旬，5 cm的土壤含水量(SWC)均大于0.13 m3/m3，SWC的平均值為0.19 m3/m3；而在6月中旬—8月上旬，SWC均小于0.12 m3/m3，平均值為0.10 m3/m3。兩個(gè)時(shí)段的土壤含水量幾乎相差1倍，因此，將2008年6月23日—8月6日選取為干旱時(shí)段，8月7日—10月15日選取為濕潤時(shí)段來分析CLM對(duì)蒸散發(fā)的模擬性能。圖3給出了干旱、濕潤條件下模式模擬值(lEclm)和觀測(cè)值(lEobs)的比較。如圖3所示，CLM能夠模擬出干旱、濕潤時(shí)段蒸散發(fā)的變化趨勢(shì)，但干旱時(shí)段的偏差相對(duì)要大于濕潤時(shí)段的。另外從不同時(shí)段的散點(diǎn)圖來看，濕潤時(shí)段兩者的散點(diǎn)均勻地分布在擬合線兩側(cè)，而干旱時(shí)段特別是蒸散發(fā)較大時(shí)的散點(diǎn)分布較分散，濕潤時(shí)段的相關(guān)系數(shù)0.77大于干旱時(shí)段的0.62，表明濕潤時(shí)段的相關(guān)性更高，因此模式在濕潤時(shí)段對(duì)蒸散發(fā)的模擬性能優(yōu)于干旱時(shí)段的。楊澤粟[14]研究CLM在黃土高原不同氣候區(qū)對(duì)陸面蒸散發(fā)的模擬效果發(fā)現(xiàn)模式在半濕潤區(qū)站點(diǎn)對(duì)蒸散發(fā)的模擬效果要優(yōu)于半干旱區(qū)的。

2.2 模擬結(jié)果的改進(jìn)

針對(duì)CLM模擬的蒸散發(fā)(潛熱通量)偏差過大的問題，考慮到植被對(duì)模擬的影響可以用NDVI表征，并結(jié)合Parr等[28]的研究思路，通過擬合NDVI和CLM模擬值與觀測(cè)值間的相對(duì)偏差，提出一種新的“模擬—校正”蒸散發(fā)的新方法。首先，利用2008年的觀測(cè)的蒸散發(fā)建立模擬—校正函數(shù)；然后，利用其余時(shí)段的CLM模擬值，驗(yàn)證方法的可行性。如圖4所示，分別給出了干旱、濕潤條件下模擬蒸散發(fā)的相對(duì)誤差隨NDVI的變化。從圖中可以看出，干旱條件下相對(duì)誤差與NDVI成正比，隨NDVI的增大而增大；濕潤條件下，相對(duì)誤差與NDVI之間滿足二次曲線變化規(guī)律，并擬合得到如下的函數(shù)關(guān)系式用于校正蒸散發(fā)。

圖2 2008年SACOL站半干旱草地日平均土壤含水量(5 cm)及降水的變化

(3)

(4)

假設(shè)CLM模擬的黃土高原區(qū)干濕季蒸散發(fā)與觀測(cè)值和NDVI之間的關(guān)系始終保持不變，那么利用CLM模擬的值，可以得到校正值可以表示為：

(5)

(6)

圖3 CLM模式模擬的SACOL站干旱、濕潤條件潛熱通量與觀測(cè)值的比較日變化、散點(diǎn)圖

圖4 CLM模式模擬的SACOL站蒸散發(fā)相對(duì)誤差隨NDVI的變化

根據(jù)上述關(guān)系式，利用2007年、2009—2012年黃土高原半干旱區(qū)SACOL站的資料，對(duì)CLM模擬結(jié)果進(jìn)行了校正。如圖5所示，校正后在濕潤、干旱條件計(jì)算的蒸散發(fā)與觀測(cè)值的變化趨勢(shì)一致，兩者之間的誤差很小。圖6給出了研究時(shí)段內(nèi)極濕年(2007年)和極干年(2011年)估算值與觀測(cè)值的散點(diǎn)圖。從圖6中可以看出，估算值與觀測(cè)值的擬合效果更好，兩者的擬合系數(shù)更接近于1，相關(guān)性較好，“模擬—校正”法對(duì)濕潤、干旱條件蒸散發(fā)的估算效果有顯著改進(jìn)。從偏差的統(tǒng)計(jì)結(jié)果(表2)來看，濕潤條件下，“模擬—校正”法計(jì)算的蒸散發(fā)的平均偏差由CLM計(jì)算的4.36 W/m2減小到2.17 W/m2，均方根偏差由27.95 W/m2減小到9.48 W/m2；干旱條件下，新方法計(jì)算的蒸散發(fā)的平均偏差由CLM計(jì)算的13.27 W/m2減小到3.09 W/m2，均方根偏差由36.81 W/m2減小到7.61 W/m2，相關(guān)系數(shù)由0.59增加到0.94?？傮w來看，利用模擬蒸散發(fā)相對(duì)誤差及NDVI的關(guān)系構(gòu)建的“模擬—校正”法對(duì)半干旱區(qū)蒸散量的估算精度有顯著提高，更準(zhǔn)確地捕捉了不同環(huán)境條件下蒸散發(fā)的變化特征。

3 結(jié) 論

CLM能較好地模擬出蒸散發(fā)的變化趨勢(shì)，降水量是影響CLM計(jì)算蒸散發(fā)精度的重要因素，在濕潤環(huán)境條件下模式對(duì)蒸散發(fā)的模擬性能較優(yōu)。引入歸一化植被指數(shù)研究植被對(duì)陸面模式模擬蒸散發(fā)的影響表明，干旱時(shí)段模擬蒸散發(fā)的相對(duì)誤差與NDVI成正比；濕潤時(shí)段相對(duì)誤差與NDVI之間滿足二次曲線變化規(guī)律。基于以上關(guān)系構(gòu)建“模擬—校正”法估算蒸散發(fā)表明，新方法在濕潤、干旱條件下計(jì)算的蒸散發(fā)與觀測(cè)值的變化趨勢(shì)一致，兩者的誤差很小，擬合效果更好。

圖5 “模擬-校正”法計(jì)算的SACOL站濕潤(2007年、2012年)、干旱條件(2009-2011年)下生長季蒸散發(fā)與觀測(cè)值的對(duì)比

圖6 “模擬-校正”法計(jì)算的SACOL站濕潤(2007年)、干旱(2011年)條件下生長季蒸散發(fā)與觀測(cè)值的散點(diǎn)圖

表2 CLM模擬與“模擬—校正”法估算的蒸散發(fā)與觀測(cè)值的統(tǒng)計(jì)參數(shù)

利用模擬蒸散發(fā)相對(duì)誤差及NDVI在干濕不同環(huán)境條件下的關(guān)系構(gòu)建的“模擬—校正”法估算蒸散發(fā)方法顯著提高了半干旱區(qū)蒸散發(fā)的估算精度，能更準(zhǔn)確地捕捉不同環(huán)境條件下蒸散發(fā)的變化特征，在一定程度上彌補(bǔ)了陸面過程模式植被生理過程參數(shù)化方案不完善引起的不確定性?；诮?jīng)驗(yàn)建立的該方法更趨向于是計(jì)算蒸散發(fā)的技術(shù)，在模擬資料中加入遙感信息，利用反映地表非均勻性的遙感信息改進(jìn)模式參數(shù)化方案不確定性引起的誤差，可以得到更為準(zhǔn)確的區(qū)域蒸散量。但是，“模擬—校正”法僅針對(duì)黃土高原特殊的氣候環(huán)境提出，對(duì)陸面蒸散發(fā)的計(jì)算效果還需要更多測(cè)站的資料來驗(yàn)證。

致 謝：感謝蘭州大學(xué)干旱氣候與環(huán)境觀測(cè)站(SACOL站)為本文提供數(shù)據(jù)支持。