閆人華，熊黑鋼，張 芳

(1．新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，新疆 烏魯木齊 830046；2．北京聯(lián)合大學(xué)應(yīng)用文理學(xué)院城市系，北京 100083；3．教育部綠洲生態(tài)重點(diǎn)實驗室，新疆 烏魯木齊 830046)

陸面過程是連接近地層大氣圈、水圈、生物圈和土壤圈之間關(guān)系的紐帶,也是它們之間進(jìn)行物質(zhì)和能量交換的橋梁以及發(fā)生相互作用的必然環(huán)節(jié)[1-2]。由于其深刻影響著區(qū)域環(huán)境的形成、全球氣候變化等問題，陸面過程的研究日益成為大氣物理學(xué)、生態(tài)學(xué)、氣象氣候?qū)W、地理學(xué)等相關(guān)學(xué)科的重點(diǎn)研究領(lǐng)域[3-4]。以地表能量平衡為主的地-氣能量交換過程是陸面過程研究的主要內(nèi)容，它直接關(guān)系到陸面(土壤)溫度的形成和變化。而土壤溫度作為陸面過程中的重要物理特征，又反作用于地-氣間能量和物質(zhì)流的循環(huán)，是研究地氣相互作用的重要參數(shù)，對氣候變化起著非常重要的作用[5-7]。因此，國外開展的大型陸面過程實驗EFEDA[8]、NOPEX[9]、BOREAS[10]等均將陸面溫度和能量平衡過程列為主要觀測項目。

自20世紀(jì)80年代在黑河流域開展國內(nèi)第一個大型陸面過程觀測試驗(HEIFE)以來，我國已先后開展了西北河西走廊干旱區(qū)(NWC-ALIEX)[11-13]、淮河流域(HUPEX)[14]、西南云貴地區(qū)(JICA)[15]、青藏高原(TIPEX)[16]、內(nèi)蒙古草原(IMGRASS)[17]等具有典型自然生態(tài)性代表區(qū)域的重大陸面過程試驗。它們一方面較好地揭示了這些不同性質(zhì)下墊面的陸面過程特征與發(fā)生機(jī)理，特別是陸面溫度和能量平衡變化特征；另一方面也為我國奠定了陸面過程研究的理論和實踐基礎(chǔ)[18]。然而目前涉及非均勻下墊面的研究還不太完善，尤其是它在不同天氣條件下陸面溫度與能量平衡之間關(guān)系的差異研究尚需深入，以揭示其在不同環(huán)境條件下的響應(yīng)程度，為陸面過程的參數(shù)化提供依據(jù)。處在綠洲與荒漠交錯位置，性質(zhì)上介于兩者生態(tài)系統(tǒng)之間的綠洲荒漠過渡帶，是減輕荒漠對綠洲壓力的天然屏障和綠洲擴(kuò)展并改造荒漠的先鋒[19-20]。而其特殊的陸面溫度變化及其與地表能量之間的關(guān)系特征則為過渡帶的穩(wěn)定和發(fā)展提供能量支撐，因此，在該區(qū)開展上述相關(guān)研究是合理配置水熱資源，科學(xué)調(diào)控人類行為的基礎(chǔ)，關(guān)系著荒漠化的進(jìn)退。本文以2010年7-9月準(zhǔn)噶爾盆地東南部的綠洲荒漠過渡帶天然芨芨草地實測數(shù)據(jù)為基礎(chǔ)，分析了該區(qū)陸面溫度典型天氣條件下的日變化規(guī)律，采用相關(guān)分析和多元逐步回歸等方法定量研究了陸面溫度與地表能量通量之間的相互關(guān)系以及不同天氣背景對兩者關(guān)系的影響，并在此基礎(chǔ)上建立了預(yù)測陸面溫度的經(jīng)驗?zāi)Ｐ汀＿@些可深化對過渡帶陸面過程模式的研究，為進(jìn)一步提高區(qū)域氣候模型的模擬水平奠定基礎(chǔ)，有利于該地區(qū)水熱資源的優(yōu)化利用。

1 研究區(qū)概況

研究區(qū)是位于天山北麓，準(zhǔn)噶爾盆地東南緣的新疆奇臺縣西北部萬木春林場附近(89°26′E， 44°11′N，海拔611 m)的大面積天然芨芨草地，向北1 km與我國第二大沙漠——古爾班通古特沙漠鄰接，屬典型的綠洲-荒漠過渡帶。該地區(qū)年平均氣溫5 ℃，年平均風(fēng)速為2.9 m/s；多年平均降水量150 mm，而年平均蒸發(fā)潛力高達(dá)2 100 mm，屬溫帶大陸性干旱氣候。

地表主要植物為多年密叢生禾本科草本植物芨芨草(Achnatherumsplendens)，總蓋度30%左右，在100 m2內(nèi)約有30叢，其間稀疏分布有苦豆子(Sophoraalopecuroides)、豬毛菜(SalsolacollinaPall)等小型植物，土壤為沙質(zhì)土，地下水埋深3～5 m。

2 研究方法

2.1 數(shù)據(jù)觀測與處理

本研究主要觀測項目有輻射(包括總輻射、凈輻射)、空氣溫濕度、風(fēng)速風(fēng)向、土壤溫濕度、土壤熱通量、降水量。土壤溫濕度測量采用一體化土壤溫濕度傳感器TRIME-PICO-32/64(IMKO，Germany),分別埋設(shè)在地下5、10、15、20、40、60、80和100 cm等共8個深度；土壤熱通量板埋設(shè)在地表以下5 cm，降水量采用RE1雨量筒傳感器(Delta-T，England)測量。

總輻射、凈輻射分別采用ES2總輻射表(測量范圍：0～2 000 W；精度：±3%)和NR2凈輻射表(測量范圍：-500～2 000 W/m2；精度：±5%)觀測,均架設(shè)在距地1.5 m處?？諝鉁貪穸群惋L(fēng)速風(fēng)向均在芨芨草灌層以上0.5和2.5 m處進(jìn)行觀測。其中，空氣溫濕度利用AT2(測量范圍：-20～80 ℃；精度：±0.5%)和RHT2(測量范圍：0～100%；精度：±2%(5%95%))傳感器測量，風(fēng)速風(fēng)向采用AN3(測量范圍：0.15～75 m/s；精度：±0.1 m/s)和WD1(測量范圍：0～358°；精度:±2°(V>5 m/s))探頭觀測，上述儀器均為Delta-T公司(England)生產(chǎn)。以上各項數(shù)據(jù)由DT80自動數(shù)據(jù)采集器每30 min采集記錄一次，最后再由計算機(jī)導(dǎo)出。

利用SPSS統(tǒng)計分析軟件、DPS數(shù)據(jù)處理系統(tǒng)對野外所采集的數(shù)據(jù)進(jìn)行整理和分析。

2.2 地表能量計算

波文比能量平衡法是國際上廣泛采用的精度較高的能量估算方法之一。并且許多學(xué)者也對該方法適用的可靠性進(jìn)行了討論[21-24]，為此方法在觀測和計算的實際應(yīng)用中就如何更好地提高精度提供了許多可借鑒的經(jīng)驗。

波文比法基于能量平衡原理和梯度擴(kuò)散方程，主要是利用溫度梯度和水汽壓梯度之比將下墊面獲得的能量分為感熱通量和潛熱通量

Rn=LE+H+G

式中，Rn為凈輻射通量(W/m2)，由輻射儀直接測出；LE為潛熱通量(W/m2)；H為感熱通量(W/m2)；G為土壤熱通量(W/m2)，由土壤熱通量板測出；γ為濕度計常數(shù)(取值為0.646[25])；Δt為上下空氣溫度差(℃)；Δe為上下水汽壓差(hPa)。

3 結(jié)果與討論

3.1 不同天氣條件下陸面溫度日變化

地氣之間能量傳輸?shù)娜談討B(tài)變化必然會反映到陸面(土壤)溫度的變化上。選取典型晴天(總云量<2成，7月22日)、陰天(總云量>8成，7月25日)、降雨天(1.2 mm/h，8月11日)，對比分析它們的土壤溫度日變化特征，以揭示不同天氣條件下土壤溫度日變化的差異性。各種天氣條件下40、60、80和100 cm深土壤溫度幾乎沒有明顯的日變化特征，始終分別維持在28、26、23和22 ℃的低水平上，受大氣能量傳輸?shù)娜兆兓绊戄^小(圖1)。因此將只討論20 cm以上土壤溫度的日變化特征。

晴天5、10和15 cm土層溫度日變化表現(xiàn)為準(zhǔn)“S”型曲線，峰值隨深度的增加而逐漸減小，其出現(xiàn)時間也逐步滯后(圖1a)。10 和15 cm土層溫度峰值依次比5 cm(40.56 ℃)低4.43 ℃、7.46 ℃，出現(xiàn)時間分別比5 cm晚約2h、5h。這表明土壤溫度對短波輻射能量加熱的響應(yīng)程度隨深度逐漸減小，主要是因為太陽輻射的直接照射，使地表溫度升溫迅速，而下層土壤則依賴于地表土壤熱量的傳遞，所以溫度上升緩慢，峰值時間延遲，數(shù)值減小。20 cm的土壤溫度則呈現(xiàn)出幅度較小的反相位變化，在中午14∶30達(dá)到最小值29.19 ℃，午夜0∶30左右出現(xiàn)最大值31.89 ℃。

雨天時，土壤溫度日變化與晴天規(guī)律相似，但振幅明顯減小，以5 cm土層為例，其變化極差為8 ℃，僅相當(dāng)于晴天極差的53%(圖1c)。而且受雨水的影響，峰值大大削弱，5～20 cm各層峰值分別比晴天減小9.40、4.27、2.94、3.01 ℃。此外，5和10 cm峰值出現(xiàn)時間分別比晴天滯后1 h、0.5 h。

陰天時，土壤溫度的日變化與晴天的日循環(huán)規(guī)律明顯不同，總體上，5、10、15 和20 cm土層溫度日內(nèi)均呈現(xiàn)出直線下降態(tài)勢。同時，這四層溫度隨深度的增加而增大(圖1b)。之所以出現(xiàn)該種情況，主要是與陰天土壤全天大部分時間在釋放熱量，自身能量虧損有關(guān)。這也說明了對于降水頻率少且強(qiáng)度小的該地區(qū)來說，陰天對土壤溫度變化的影響程度遠(yuǎn)大于雨天。

平均的土壤溫度日變化也和晴天保持較高的一致性(圖1d)。這是由于夏季干旱區(qū)以晴朗天氣為主，陰雨天較少的緣故。

圖1 不同天氣條件下各層土壤溫度日變化特征

土壤活動層5～20 cm的溫度梯度差在一定程度上反映了地氣之間能量傳遞的方向性。除陰天外，5～20 cm的溫度梯度差都很強(qiáng)，最大值達(dá)到10.88 ℃，但其峰值表現(xiàn)為晴天>平均>雨天(圖2)。平均情況下的溫度差日變化與晴天最為一致，雨天也較為相似，都表現(xiàn)出正弦曲線形態(tài)。晴天和平均在12∶00-24∶00之間溫度梯度>0，熱量向下傳輸，土壤貯藏?zé)崃浚?∶00-12∶00之間梯度<0，土壤放出熱量。雨天由于受降水低溫的影響，梯度由負(fù)值轉(zhuǎn)為正值較晴天滯后2 h、由正值轉(zhuǎn)為負(fù)值提前2.5 h。陰天土壤溫度差則在大部分時間里小于0，表明陰天由于太陽輻射的不足，土壤幾乎全天作為熱源釋放熱量來彌補(bǔ)地氣系統(tǒng)能量的不足。

圖2 不同天氣條件下5～20 cm土壤溫度差

3.2 陸面溫度與地表能量關(guān)系

3.2.1 不同深度土壤溫度與地表能量關(guān)系 為更準(zhǔn)確、深入地探討土壤溫度與地表能量之間的關(guān)系，本文將對土壤溫度進(jìn)行分層討論。因地表能量傳輸主要發(fā)生在日間，故將觀測期內(nèi)地表能量日間(8∶00-18∶00)數(shù)據(jù)與各層土壤溫度分別進(jìn)行相關(guān)性分析，并作顯著性水平檢驗(表1)。

各表層土壤溫度與地表能量分量的關(guān)系均呈現(xiàn)出:相關(guān)系數(shù)隨深度的增加而逐漸減小，由正相關(guān)轉(zhuǎn)為負(fù)相關(guān)，且負(fù)相關(guān)趨勢增加的規(guī)律(表1)。這主要是由于隨著土壤深度的增加，土溫峰值出現(xiàn)的時間逐漸滯后，從而導(dǎo)致它們與能量日變化曲線不同步性愈來愈強(qiáng)，相關(guān)系數(shù)變小并轉(zhuǎn)為負(fù)數(shù)。到20 cm時，土壤溫度日變化曲線呈反相位變化，與各能量變化趨勢相反，兩者負(fù)相關(guān)性達(dá)到最大。表層5 cm土壤溫度與各能量分量相關(guān)性最好，對其變化反應(yīng)靈敏性最強(qiáng)。其與土壤熱通量(G)的相關(guān)性分別是10、15 和20 cm、平均的1.61、6.81、2.72和1.54倍。因為土壤表層是SPAC(土壤—植被—大氣)系統(tǒng)的臨界面，能量交換首先從表層土壤開始。在下面的討論中，將只討論5 cm深土層溫度與各能量分量的關(guān)系。

表1 不同深度土壤溫度與各地表能量分量之間的相關(guān)系數(shù)1)

雖然Rn、LE、H、G與土壤溫度的線性相關(guān)性均達(dá)到極顯著性水平，但G最好，Rn和H居中，LE最差(表1)。這主要是由于G代表熱量進(jìn)出土壤的多少，直接影響著土壤溫度的變化。白天G為正，且持續(xù)增大時，熱流不斷由地表向土壤深層傳輸，土壤能量盈余，使得溫度不斷升高。而Rn、LE、H則是通過作用于G，對土壤溫度產(chǎn)生間接影響。這表明綠洲荒漠過渡帶土壤溫度的變化是以G變化為基礎(chǔ)，同時還受其它能量相互影響而發(fā)生的。

3.2.2 不同天氣條件下地表能量與土壤溫度關(guān)系 晴天G與土壤溫度關(guān)系最為密切，陰天略差，雨天最小(表2)。陰天因為一方面云量與土壤溫度呈負(fù)相關(guān)，削弱了G與土壤溫度關(guān)系；另一方面受云層的影響，凈輻射少，氣溫較低，致使G數(shù)值較小，土溫變化緩慢，造成兩者相關(guān)性R2低于晴天19%(表2b)。雨天則由于水的比熱容較大，降水(RF)對土壤溫度的冷卻作用明顯，與其呈顯著負(fù)相關(guān)，從而干擾了G和土壤溫度變化的同步性，使得兩者相關(guān)性僅為晴天的49%(表2c)。平均情況與晴天相當(dāng)都在0.61左右，表明在觀測期內(nèi)由于晴天較多，受云和降水影響較少，兩者關(guān)系與晴天較為一致(表2d)，這也是干旱區(qū)所特有的特征。土壤溫度對G變化的反應(yīng)強(qiáng)度(回歸系數(shù))也呈現(xiàn)為：晴天和平均相當(dāng)，陰天其次，雨天最小。

與G類似，Rn、LE與土壤溫度線性相關(guān)性均表現(xiàn)為：晴天>平均>陰天>雨天(表2)。這主要是因為晴天相比其它天氣，Rn、LE數(shù)值較大，引起土壤溫度變化較為明顯，致使與土壤溫度關(guān)系密切(表2a)。受云層和降水對Rn、LE的干擾，陰天和雨天相關(guān)性依次減弱，均低于晴天(表2b.c)。除晴天外，LE與土溫相關(guān)性均未達(dá)到顯著性水平，表明其對土壤溫度影響微弱。

H與土壤溫度線性相關(guān)性表現(xiàn)為：平均>雨天>晴天>陰天(表2)。土壤溫度對H變化的反應(yīng)強(qiáng)度(回歸系數(shù))卻表現(xiàn)出：平均>晴天>雨天>陰天。

表2 不同天氣條件下各地表能量分量與土壤溫度(T土-5cm)關(guān)系

3.3 陸面溫度預(yù)測模型的建立與檢驗

由于該區(qū)以晴天為主，陰雨天較少，因此從觀測期內(nèi)各選取典型晴天13天、陰天3天、雨天2天的日內(nèi)半小時值數(shù)據(jù)作為樣本，采用多元逐步回歸方法分別建立以土壤溫度為因變量，地表能量分量為自變量(陰天和雨天可能分別含有云量、降水量因子)的陸面溫度預(yù)測模型。不同天氣再另外各選取1 d相關(guān)數(shù)據(jù)作為樣本，進(jìn)行模型預(yù)測效果的檢驗。

從P值看，各回歸模型都達(dá)到了極顯著性水平，但晴天和平均建模效果較好，決定系數(shù)(R2)分別為0.809 3、0.947 1；陰天次之，R2為0.662 3；雨天建模效果不如前三者理想，R2為0.211 2(表3)。預(yù)測陸面溫度與實測陸面溫度的相關(guān)性檢驗與此相似。各模型的預(yù)測值與實測值相關(guān)性均達(dá)到極顯著性水平(P均小于0.01)，且斜率接近于1，平均偏差(MBE)較小，但從R2和斜率來看，仍表現(xiàn)為晴天最好，陰天其次，雨天相對略差(圖3)。

表3 不同天氣條件下陸面溫度預(yù)測模型的建立

圖3 陸面溫度模型預(yù)測值與實測值的相關(guān)性檢驗

經(jīng)過檢驗的預(yù)測模型表明，無論何種天氣條件下土壤溫度均主要受干熱性能量因子(Rn、G、H)的影響，而與生態(tài)系統(tǒng)因水份蒸發(fā)而引起的潛熱(LE)能量交換關(guān)系較小(表3)。這充分說明在降水稀少、水分短缺的干旱區(qū)，綠洲荒漠過渡帶SPAC系統(tǒng)(土壤-植被-大氣連續(xù)體)實際蒸散微弱，對土壤溫度影響較小。而在陰天和雨天，云量、降水因子則由于已通過地表能量受其影響后的大小變化來反映出其對陸面溫度的作用，因此在逐步回歸模型中未被引入。

4 結(jié) 論

1)晴天表層土壤(5、10、15 cm)溫度變化表現(xiàn)為標(biāo)準(zhǔn)的“S”曲線，峰值的大小和時間因土壤深度的增加而分別減小、滯后。而20 cm土壤溫度卻出現(xiàn)反相位變化。陰天由于較多的土壤熱量向大氣釋放，總體上，5～20 cm土壤溫度都大致表現(xiàn)出隨時間變化的直線遞減特征。雨天土壤溫度日變化形態(tài)與晴天相似，但振幅明顯偏小，5 cm土溫變化極差僅為晴天的53%。平均土壤溫度日變化也與晴天保持較高的一致性。無論何種天氣條件下，40 cm及其以下深層土壤溫度幾乎沒有日變化。

2)隨著深度的增加，表層(0～15 cm)土溫峰值出現(xiàn)的時間以2.5h/5cm速度推延，使得它們與能量日變化曲線不同步性愈來愈強(qiáng)。這些造成各能量分量與土壤溫度關(guān)系均表現(xiàn)為相關(guān)系數(shù)隨土壤深度的增加而逐漸減小，由正相關(guān)轉(zhuǎn)為負(fù)相關(guān)，到20 cm負(fù)相關(guān)性達(dá)到最大的規(guī)律。

3)各地表能量分量與土壤溫度相關(guān)性雖都達(dá)到極顯著性水平，但表現(xiàn)為G>Rn>H>LE。這主要是由于G直接代表進(jìn)出土壤熱量的多少，而其它能量則是通過G對土壤溫度進(jìn)行間接影響。

4)不同天氣條件下，G、Rn、LE與土壤溫度相關(guān)性均表現(xiàn)出晴天>平均>陰天>雨天，土壤溫度對其能量變化的反應(yīng)強(qiáng)度(回歸系數(shù))也均表現(xiàn)出：晴天和平均相當(dāng)，陰天其次，雨天最差。雨天降水對土壤溫度具有冷卻作用，與其呈顯著負(fù)相關(guān)(R=-0.338)。陰天云量也與土壤溫度呈負(fù)相關(guān)(R=-0.530)。

5)適合不同天氣條件的陸面溫度預(yù)測模型表明：由于干旱區(qū)SPAC系統(tǒng)實際蒸散微弱，綠洲荒漠過渡帶土壤溫度主要受干熱性能量因子(Rn、G、H)的影響，與濕熱性能量因子(LE)關(guān)系較小。

參考文獻(xiàn)：

[1]劉宏誼，楊興國，張強(qiáng)，等.敦煌戈壁冬夏季地表輻射與能量平衡特征對比研究[J]．中國沙漠，2009，29(3)：558-565．

[2]曾劍，張強(qiáng)，王勝.中國北方不同氣候區(qū)晴天陸面過程區(qū)域特征差異[J]．大氣科學(xué)，2011，35(3)： 483-494．

[3]周余華，葉伯生，胡和平.土壤凍融條件下的陸面過程研究綜述[J]．水科學(xué)進(jìn)展，2005，16(6)：887-891．

[4]丁張巍，溫之平，黎偉標(biāo)，等.廣東龍門地區(qū)草地下墊面地表輻射和能量平衡觀測和分析[J]．中山大學(xué)學(xué)報：自然科學(xué)版，2011，50(3)：123-129．

[5]孫昭萱，張強(qiáng).黃土高原半干旱區(qū)陸面溫度和能量的氣候特征分析[J]．中國沙漠，2011，31(5)：1302-1308．

[6]劉輝志，涂鋼，董文杰，等.半干旱地區(qū)地氣界面水汽和二氧化碳通量的日變化及季節(jié)變化[J]．大氣科學(xué)，2006，30 (1)：108-118．

[7]代成穎，高志球，王琳琳，等.兩種土壤溫度算法的對比分析[J]．大氣科學(xué)，2009，33 (1)：135-144．

[8]BOLLE H J，ANDERE J C，ARRUE，et al.EFEDA-European field experiment in a desertification-threatened area[J]．Annual Geophysics，1993，11(2)：173-189．

[9]HALLDIN S，GOTTSCHALK L，GRIEND A A，et al.NOPEX-Anorthern hemisphere climate processes land surface experiment[J]．Journal of Hydrology，1998，212(1)：172-187．

[10]SELLERS P，HALL F，RANSON K J，et al.The boreal ecosystem atmosphere study (BOREAS)：An overview and early results from the 1994 field year[J]．Bulletin of the American Meteorological Society，1995，76(9)：1549-1577．

[11]張強(qiáng)，黃榮輝，王勝，等.西北干旱區(qū)陸—?dú)庀嗷プ饔迷囼?(NWC-ALIEX)及其研究進(jìn)展[J]．地球科學(xué)進(jìn)展，2005，20(4)：427-441．

[12]曹曉彥，張強(qiáng).西北干旱區(qū)荒漠戈壁陸面過程的數(shù)值模擬[J]．氣象學(xué)報，2003，61(2)： 219-225．

[13]張強(qiáng)，曹曉彥.敦煌地區(qū)荒漠戈壁地表熱量和輻射平衡特征[J]．大氣科學(xué)，2003，27(2)：245-254．

[14]周小剛，羅云峰.“九五”基金重大項目“淮河流域能量與水分循環(huán)試驗和研究”介紹與回顧[J]．中國科學(xué)基金，2004，18(4)：226-228．

[15]施曉暉，徐祥德，程興宏.2008年雪災(zāi)過程高原上游關(guān)鍵區(qū)水汽輸送機(jī)制及其前兆性“強(qiáng)信號”特征[J]．氣象學(xué)報，67(3)：478-487．

[16]徐祥德，陳聯(lián)壽.青藏高原大氣科學(xué)試驗研究進(jìn)展[J]．應(yīng)用氣象學(xué)報，2006，17(6)：756-772．

[17]呂達(dá)仁，陳佐忠，陳家宜，等.內(nèi)蒙古半干旱草原土壤—植被—大氣相互作用(IMGRASS)綜合研究[J]．地學(xué)前緣，2002，9(2)：295-306．

[18]張強(qiáng)，胡向軍，王勝，等.黃土高原陸面過程試驗研究(LOPEX)有關(guān)科學(xué)問題[J]．地球科學(xué)進(jìn)展，2009，24(4)：363-371．

[19]張凱，王潤元，張強(qiáng)，等.綠洲荒漠過渡帶夏季晴天地表輻射和能量平衡及小氣候特征[J]．中國沙漠，2007，27(6)：1055-1061．

[20]張勃，張華，張凱，等.黑河中游綠洲及綠洲—荒漠生態(tài)脆弱帶土壤含水量空間分異研究[J]．地理研究，2007，26(2)：321-327．

[21]JENSEN M E，BURMAN R D，ALLEN R G.Evapotranspiration and irrigation water requirements.ASCE Manuals and Reports on Engineering Practices No.70 [ M]．ASCE，New York，NY，1990．

[22]TODD R W，EVETT S R，HOWELL T A.The Bown ratio energy balance method for estimating latent heat flux of irrigated alfalfa evaluated in a semiarid，advective environment [J]．Agricultural and Forest Meteorology，2000，103：335 -348．

[23]GAVILAN P，BERENGENA J.Accuracy of the Bowen ratio-energy balance method for measuring latent heat flux in a semiarid advective environment[J]．Irrigation Science，2007，25(2)：127-140．

[24]XING Z H， CHOW L E， MENG F R，et al.Validating evapotranspiraiton equations using Bowen Ratio in New Brunswick，Maritime，Canada[J]．Sensors，2008，8：412-428．

[25]ALLEN R G，PEREIRA L S，RAES D，et a1.FAO irrigation and drainage paper No．56，crop evapotranspiration( guidelines for computing crop water requirements) [M]．FAO，Rome，1998．