陳 晗，張翔宇

(1.重慶交通大學(xué) 河海學(xué)院，重慶 400074；2.華北水利水電大學(xué)，河南 鄭州 450045)

?

渦動相關(guān)法在蒸散發(fā)實測中的運用

陳 晗1，張翔宇2

(1.重慶交通大學(xué) 河海學(xué)院，重慶 400074；2.華北水利水電大學(xué)，河南 鄭州 450045)

在區(qū)域水循環(huán)各個環(huán)節(jié)中，蒸散發(fā)是占比最大最難估算的分量，如何對其進行長時間、高精度、高頻率的測定成為難題。近年來隨著微氣象學(xué)理論和量子激光技術(shù)的逐漸成熟，利用渦動相關(guān)法測定蒸散發(fā)成為可能。湖南省作為典型的中國南方地區(qū)，農(nóng)業(yè)主要以蔬菜等經(jīng)濟作物以及水稻為主，在總的用水量中，農(nóng)業(yè)用水占到了總用水量的70%，并且在農(nóng)業(yè)總失水量中蒸散發(fā)是很大的一部分，所以對蒸散發(fā)準確的測定是十分必要的。本文簡要介紹了渦動相關(guān)法的相關(guān)原理，利用EC-QCL分析儀和三維超聲風(fēng)速儀，以及其他一些輔助儀器對岳陽市廣興洲鎮(zhèn)的一個大面積包菜地進行了實地蒸散發(fā)的測定，并利用EddyPro軟件計算了2015年的蒸散量，分析了研究區(qū)域蒸散量在不同時間尺度的變化趨勢。本文采用的研究方法可為以后蒸散發(fā)的實測工作提供參考與借鑒。

蒸散發(fā)；渦動相關(guān)；實測；Eddypro軟件

蒸散發(fā)包括地表和植物表面的水分蒸發(fā)以及植物表面和植物體內(nèi)的水分蒸騰兩部分。蒸散發(fā)不僅是生態(tài)圈水循環(huán)重要的組成部分，同時也是陸地以及大氣之間能量循環(huán)的重要組成部分。到達地球表面的降水有70%通過蒸發(fā)或蒸騰作用回到大氣中，在干旱區(qū)這個數(shù)字可能達到90%。在當(dāng)前普遍面臨水資源短缺的大背景之下，準確測定蒸散量，對優(yōu)化區(qū)域水資源配置、水資源規(guī)劃與管理，以及更深入了解陸面過程，正確評估氣候和人類活動對生態(tài)系統(tǒng)的影響具有重要的現(xiàn)實與經(jīng)濟意義[1-3]。

1 渦動相關(guān)法簡介

渦動相關(guān)(eddy covariance technique，EC)是指某種物質(zhì)的垂直通量，即這種物質(zhì)的濃度與其垂直速度的協(xié)方差。渦動相關(guān)法提供了一種直接測量植被與大氣間CO2、水和熱通量的方法。渦動相關(guān)法的優(yōu)點在于能通過測量各種屬性的湍流脈動值來直接測量它們的通量與其他方法相比，它并不是建立在經(jīng)驗關(guān)系基礎(chǔ)之上，或從其他氣象參量推論而來，而是建立在所依據(jù)的物理原理之上，是一種直接測量亂流通量的方法，是各種實測方法中較精密和可靠的方法，并且該方法有利于蒸散發(fā)觀測工作的長期開展。

渦動相關(guān)法是根據(jù)高頻(一般10 Hz)測定的垂直風(fēng)速和氣體濃度來計算通量[4]。通量計算為：

(1)

式(1)中:Fc代表氣體的通量,g/(h·m2)；ρa代表空氣密度,g/m3；w′和c′分別代表垂直風(fēng)速(m/s)和氣體濃度的脈動，即w′為瞬時垂直風(fēng)速測量值偏離給定時間段內(nèi)風(fēng)速算數(shù)(或滑動)平均值的量，c′類推。上劃線代表計算平均通量的時間長度(一般為30 min，也有用1 h)。

2 材料與方法

2.1 實驗地點概況

本研究的實驗地點位于湖南省岳陽市君山區(qū)廣興洲鎮(zhèn)五一村，為中國科學(xué)院洞庭湖濕地生態(tài)系統(tǒng)觀測研究站(以下簡稱洞庭湖站)周邊的一塊典型施肥蔬菜農(nóng)田。洞庭湖站位于湖南省岳陽市郊區(qū)采桑湖南岸(北緯29°30′，東經(jīng)112°48′)，屬典型的亞熱帶季風(fēng)氣候，多年平均氣溫16.5～17.0 ℃，1月平均氣溫3.8～4.7 ℃；7月平均氣溫29 ℃左右，年平均降水量1250～1450 mm[5]。試驗點秋冬季盛行風(fēng)向為東北風(fēng)(圖1)。

圖1 試驗站2014年9月—2015年2月風(fēng)向頻率分布

實心圓表示渦動相關(guān)觀測系統(tǒng)(EC)和自動氣象站的位置，正方形為臨時搭建小屋的位置。 圖2 試驗點的樣地設(shè)置和觀測儀器位置示意

圖2為試驗點的樣地設(shè)置和觀測儀器位置示意圖。試驗點東西長102 m,南北寬56.3 m。渦動相關(guān)觀測系統(tǒng)(EC)位于試驗點的中央，小屋距離EC塔的距離為10 m，自動氣象站距離EC的距離為5 m。該地區(qū)最主要的農(nóng)作物是包菜、辣椒、南瓜、甜瓜和西瓜。通常一片地上每年8月到次年1月種植一茬作物，這一茬的作物大多為包菜；次年2月到7月種一茬作物，通常為辣椒、南瓜或包菜。

2.2 H2O通量的EC-QCL法觀測

在試驗樣地中主風(fēng)向的下游安裝EC-QCL觀測設(shè)備，包括：三維超聲風(fēng)速儀和數(shù)據(jù)采集器安裝在不銹鋼制作的觀測塔上，其中超聲風(fēng)速儀垂直安裝在2 m的高度；QCL氣體分析儀，為雙光路量子級聯(lián)可調(diào)諧紅外激光差分吸收光譜儀,該儀器及與之配套的真空抽氣泵均裝在距離觀測塔10 m且位于風(fēng)向頻率最小方位的臨時小屋內(nèi)，小屋內(nèi)配有空調(diào)用于恒溫除濕，同時配有大功率的不間斷電源，以免隨機性事故斷電對設(shè)備的影響；QCL氣體分析儀的樣氣采樣口安裝在位于三維超聲風(fēng)速儀的下方10 cm處，樣氣經(jīng)長15 m(內(nèi)徑3/4英寸)的特氟龍采樣管進入儀器，再由連接在儀器上的真空泵抽出。真空泵驅(qū)動使采樣口的樣氣以15 L/ min的流量經(jīng)QCL氣體分析儀，并使儀器樣品池(sample cell)中的氣體置換頻率高到能滿足以10 Hz的頻率測定有效H2O濃度要求。為除去樣氣中的顆粒物，采樣口下游20 cm處氣路上安裝粗過濾膜(直徑50 mm，孔徑5)，在QCL儀上游20 cm處安裝細過濾膜(直徑50 mm，孔徑0.45)，濾膜每4～5 d更換1次。為避免水汽在采樣管內(nèi)壁凝結(jié)和減小管壁對氣體的黏滯效應(yīng)，采樣管用不透光的保溫材料包裹，內(nèi)部用加熱帶進行加熱，通過自動控溫裝置讓管路的溫度保持30 ℃左右。

各個儀器的組成以及連接可參考下圖3。

圖3 渦動相關(guān)系統(tǒng)連接圖

2.3 H2O通量計算

通量用EddyPro軟件進行計算。

(1)去野點。在計算通量之前，根據(jù)Vickers和Mahrt(1997)[6]的方法對風(fēng)速和H2O濃度進行去野點處理。首先進行連續(xù)性檢驗，去除明顯不合物理意義的異常值。然后對每個變量的原始數(shù)據(jù)序列以n=10為窗口大小計算標準偏差σ，在每個窗口內(nèi)，若某個記錄值與窗口平均值的偏差大于5.5σ，則視為野點并刪除。若連續(xù)出現(xiàn)此類數(shù)據(jù)個數(shù)≥4，則視為有效數(shù)據(jù)予以保留。視為野點的數(shù)據(jù)被刪去，不再用于通量的計算。若半小時內(nèi)野點數(shù)據(jù)大于10%，則該半小時內(nèi)的通量被視為無效，當(dāng)缺測處理。

(2)坐標旋轉(zhuǎn)。EC法通量觀測有一個前提條件，即垂直風(fēng)速(w)的半小時平均值為零。但實際觀測中，由于地形、儀器安裝等因素而無法滿足這一要求。本研究采用二次坐標旋轉(zhuǎn)的方法對風(fēng)速進行修正，使其垂直風(fēng)速滿足平均值為零。

(3)信號延遲校正。因QCL氣體分析儀采用的是閉路測量原理，樣氣從采樣口經(jīng)過管路到達檢測器需要一定的時間。因此，測得的H2O濃度和垂直風(fēng)速的信號之間存在時間延遲，計算通量前，需要確定與垂直風(fēng)速同步的H2O濃度，這就需要按延遲時間對H2O濃度數(shù)據(jù)進行測量時間的修正?？赏ㄟ^計算H2O濃度和垂直風(fēng)速的協(xié)方差，用最大協(xié)方差出現(xiàn)的時刻來確定延遲時間[7]。

(4)信號去趨勢。通過信號去趨勢，建立垂直風(fēng)速和濃度的脈動時間序列。對于垂直風(fēng)速，采用了線性去趨勢(LDT)的方法。而對于H2O濃度信號，對它進行自相關(guān)滑動平均濾波處理，再用原始濃度減去滑動平均后的濃度，得到H2O的脈動時間序列。

(5)通量計算。根據(jù)式(1)計算半小時平均的H2O通量。

(6)高頻衰減修正。由于閉路氣體測量系統(tǒng)的管路衰減、儀器響應(yīng)能力不足、風(fēng)速和濃度傳感器的空間分離等原因，會造成氣體濃度湍流協(xié)譜(co-spectra)高頻部分的衰減，從而導(dǎo)致實測通量被系統(tǒng)性低估[8]，需要對這部分系統(tǒng)性低估進行修正。

(7)WPL修正。用EC技術(shù)觀測微量氣體成分的湍流通量，需要考慮因熱量或水汽通量的輸送而引起的微量氣體密度變化，因而要對氣體通量進行密度效應(yīng)修正[9]。

2.4 H2O通量數(shù)據(jù)的質(zhì)量控制處理

半小時通量數(shù)據(jù)質(zhì)量控制處理一般包括以下步驟。

(1)剔除有降水記錄或?qū)τ^測塔有操作的半小時通量數(shù)據(jù)。

(2)根據(jù)Foken和Wichura[10]的方法進行總體湍流特性檢驗，它用來判斷由于湍流混合不充分而被低估的通量，通過比較觀測到的和模擬(莫寧—奧布霍夫相似理論)到的總體湍流特性(歸一化的標準差)進行估計，如果差異超過100%，則半小時通量被剔除。

3 數(shù)據(jù)與結(jié)果

QCL儀器在運行的過程中，一般每隔15 d左右下載數(shù)采的數(shù)據(jù)，將下載的原始數(shù)據(jù)先用LoggerNet轉(zhuǎn)化成TOA5格式后，分割成半小時的小文件，然后用EddyPro軟件進行通量的計算，最后將每次計算的結(jié)果存到一個文件中。實驗起始于2014年8月，終止于2016年10月。其中2015年的數(shù)據(jù)較為完整，數(shù)據(jù)質(zhì)量較高，我們選取2015年內(nèi)的數(shù)據(jù)進行蒸散量在不同時間尺度的變化趨勢分析。

圖4 水汽通量日變化趨勢(2015年8月15日)

圖4為水汽通量在2015年8月15日的變化趨勢圖，可以看到水汽通量具有明顯的日變化趨勢，在夜晚水汽通量接近于0，在正午以及下午時達到最高，日最大值出現(xiàn)在下午16:30，達到了5.6 mmol/(s·m2)，這是由于白天太陽輻射和氣溫較高，蒸散量較大，而在夜晚太陽輻射為零，蒸散量也接近于零。

圖5表示了不同天氣下的潛熱和水汽通量日變化趨勢，從圖中可以清楚地看到，潛熱通量和水汽通量具有完全相同的變化趨勢，這是由于他們只是表征蒸散量的不同單位。在不同的天氣，潛熱通量和水汽通量差異明顯，其中，8月6日—9日下雨，8月10日—14日為多云，8月15日—17日為晴天，相對應(yīng)的，潛熱通量和水汽通量排放峰值在雨天較低，在晴天較高。8月10日正午時的通量分別只有40.41 w/m2和0.93 mmol/(s·m2)，而8月15日正午時通量分別達到了453.54 w/m2和10.37 mmol/(s·m2)，這是由于不同天氣下氣溫和太陽短波輻射差異導(dǎo)致，在晴天太陽短波輻射和氣溫較高，相應(yīng)的蒸散量也較高，而在陰天短波輻射和氣溫較低，對應(yīng)的蒸散量也處于一個較低的水平。

2015年，實際蒸騰的年內(nèi)變化趨勢見圖6。在冬季1月維持較低水平，只有0.38 mm/d；從2月

圖5 潛熱通量和水汽通量在不同天氣變化趨勢(2015.8.6—2015.8.17)

圖6 實際蒸散量在2015年內(nèi)變化趨勢

份種植季的開始，由于植物蒸騰的參與，實際蒸散量開始逐漸增大，3、4、5月持續(xù)上升，直到6月蒸散量明顯降低，6月平均蒸散量只有1.3 mm/d。這是因為6月份實驗地開始收菜，缺少了植物蒸騰導(dǎo)致(這也從側(cè)面印證了植物蒸騰在總蒸散量中的主導(dǎo)地位)。下半年從7月開始逐漸增高，8月份由于種植季的開始持續(xù)升高，直到10月份達到全年最高值，為2.6 mm/d，而到了11月和12月又呈現(xiàn)下降趨勢。 經(jīng)計算，實驗地2015年總蒸散量521 mm，日平均蒸散量1.4 mm/d。

4 結(jié) 論

在眾多蒸散發(fā)的實測方法中，渦動相關(guān)法是可以直接測量亂流通量的方法，它的優(yōu)點是采樣頻率高，可以獲得長時間高分辨率的蒸散量；缺點是對下墊面的要求較高，需要均一的下墊面條件；其次就是在夜晚由于湍流發(fā)展不充分，在夜晚的數(shù)據(jù)有很大誤差。但總體來說，渦動相關(guān)法是蒸散發(fā)各種實測方法中較精密和可靠的方法，因此常常將渦動相關(guān)技術(shù)作為一種標準值。本文通過對岳陽市廣興洲鎮(zhèn)的一大塊包菜地進行蒸散發(fā)的實地測量，證明了這種方法的適用性。結(jié)果表明，實測蒸散量在不同的時間尺度具有明顯的變化趨勢。以年為尺度，實際蒸散量在1月份較小，從2月份種植季開始逐漸增大，到5月份達到最高，6月份開始減小，下半年種植季從8月份開始蒸散量逐漸增高，10月份達到最大，11月份開始逐漸減小；以月為尺度，蒸散量在晴天較高，在陰天和雨天較低；以天為尺度，蒸散量在正午以及下午時處于每天的高峰期，而在其他時間較低。

[1] 汪明霞,陳曉飛,王鐵梁,等.騰發(fā)量的測定和計算方法研究綜述[J].中國農(nóng)村水利水電,2006(12):9-12.

[2] 劉晨峰,張志強,查同剛,等.渦度相關(guān)法研究土壤水分狀況沙地楊樹人工林生態(tài)系統(tǒng)能量分配和蒸散日變化的影響[J].生態(tài)學(xué)報,2006,26(8):2549-2557.

[3] 辛?xí)灾?田國良,柳欽火.地表蒸散定量遙感的研究進展[J].遙感學(xué)報,2003,7 (3):233-240.

[4] Aubinet M, Vesala T, Papale D. Eddy Covariance: A Practical Guide to Measurement and Data Analysis[J].Springer, 2012:365-376.

[5] 侯志勇, 謝永宏, 陳心勝,等. 洞庭湖濕地的外來入侵植物研究[J]. 農(nóng)業(yè)現(xiàn)代化研究, 2011, 32(6):744-747.

[6] V Seabury J, Mclouth L, Blodgett P. Quality control and flux sampling problems for tower and aircraft data[J]. Journal of Atmospheric & Oceanic Techno.log.y, 1997, 14(3):514-526.

[7] Wienhold F G, Frahm H, Harris G W. Measurements of N 2 O fluxes from fertilized grassland u.sin.g a fast response tunable diode laser spectrometer[J]. Journal of Geophysical Research Atmospheres, 1994, 99(D8):557-16.

[8] Moore C J. Frequency response corrections for eddy correlation systems[J]. Boundary-Layer Meteoro.log.y, 1986, 37(1):17-35.

[9] Webb E K, Pearman G I, Leuning R. Correction of flux measurements for density effects due to heat and water vapour transfer[J]. Quarterly Journal of the Royal Meteoro.log.ical Society, 1980, 106(447):85-100.

[10] Foken T.Wichura B.Tools for quality assessment of surface-based flux measurements[J].Agritulcural and Forest Meteoro.log.y.1996(78)：83-105.

Application of eddy covariance method inevapotranspiration measurement

CHEN Han1,ZHANG Xiangyu2

(1.ChongqingJiaotongUniversityHehaiInstitute,Chongqing400074，China;2.NorthChinaUniversityofwaterresourcesandhydropower,Zhengzhou450045,China)

The evapotranspiration is the most important and the most difficult to estimate the component in the regional water cycle, how to measure and estimate accurately is the focus of research on evapotranspiration. With the micro-meteorological theory and quantum laser technology gradually mature in recent years, the use of eddy covariance method for the determination of evapotranspiration become possible. Hunan Province as a typical southern region of China, agriculture mainly vegetables and other economic crops and rice-based, agricultural water accounted for over 70% of the total water consumption, and evapotranspiration accounts for over half of the total water loss, so it is very necessary to measure and estimate accurately. This paper introduce the relevant theory of eddy covariance method, we use EC-QCL analyzer, three-dimensional ultrasonic anemometer and other auxiliary instruments measure the actual evapotranspiration in a large vegetable farmland over two years observation in the Guangxingzhou, Yueyang. Then we use EddyPro to calculate the actual evapotranspiration, and analysis the change trend of actual evapotranspiration in different time scales. The research method used in this paper can provide reference for the future measurement of evapotranspiration.

evapotranspiration; eddy covariance; actual measurement; EddyPro software

陳 晗(1992-)， 男， 碩士研究生， 研究方向為水文與水資源。

P339；S161.4

A

2096-0506(2016)11-0018-05