陳 晗，黃津輝

(1.重慶交通大學河海學院，重慶 400074;2.南開大學環(huán)境科學與工程學院，天津 300350)

中國南方農(nóng)田蒸散量實測及其影響因素分析

陳 晗1，黃津輝2

(1.重慶交通大學河海學院，重慶 400074;2.南開大學環(huán)境科學與工程學院，天津 300350)

利用基于閉路QCLAS-EC激光分析儀的渦度相關(guān)法對湖南省岳陽市郊區(qū)的一片蔬菜地的實際蒸散發(fā)、水汽通量以及潛熱通量進行了連續(xù)兩年的野外觀測,并對原始觀測數(shù)據(jù)進行處理，計算蒸散量，以分析研究區(qū)域?qū)嶋H蒸散量的年際尺度周期變化趨勢和季節(jié)變化規(guī)律。結(jié)果表明，實驗區(qū)域全年蒸散量在730～803 mm之間；季節(jié)變化大，8月達到全年最高水平3.5 mm/d，而1月只有0.4 mm/d。整個研究區(qū)域?qū)崪y蒸散量變化規(guī)律是：春夏季較高，秋冬季較低；種植季較高，非種植季較低；每天的正午以及下午較高，夜晚較低。此外，利用高分辨率的實測工具觀察不同的農(nóng)業(yè)活動對實際蒸散的影響，發(fā)現(xiàn)種植引起蒸散發(fā)上升，收割引起蒸散發(fā)下降，并且收割使蒸散發(fā)下降的幅度要明顯高于種植使蒸散發(fā)上升的幅度；灌溉對蒸散發(fā)有促進作用但是作用較小，施肥對蒸散發(fā)的影響目前尚不明確。

蒸散發(fā)；南方農(nóng)田；蔬菜地；渦度相關(guān)法

人類對水循環(huán)特點及水資源的了解程度，直接影響流域水資源利用、開發(fā)、管理的合理性與科學性[1]。農(nóng)業(yè)是各個行業(yè)中第一用水大戶，而灌溉用水大部分通過蒸散發(fā)返回大氣中，因此研究農(nóng)田蒸散發(fā)(ET)，可以為農(nóng)業(yè)節(jié)水和高效水資源利用提供科學依據(jù)。

蒸散發(fā)是水分從地表傳輸?shù)酱髿庵械乃诌\移過程，是植物與外界環(huán)境水分交換的一種普遍現(xiàn)象。蒸散發(fā)不僅是生態(tài)圈水循環(huán)重要的組成部分，也是陸地以及大氣之間能量循環(huán)的重要組成部分。準確測定蒸散發(fā)量，對優(yōu)化區(qū)域水資源配置、水資源規(guī)劃與管理，以及對更加深入地了解陸面生態(tài)水文過程，正確評估氣候和人類活動對生態(tài)系統(tǒng)的影響具有重要的現(xiàn)實意義[2-5]。近年來隨著微氣象學理論以及量子激光技術(shù)的逐步成熟，利用微氣象學中渦度相關(guān)法高精度測定蒸散發(fā)成為可能。

湖南省是我國農(nóng)業(yè)大省，位于亞熱帶季風區(qū)，多年平均降水量為1 427 mm。從總量看，湖南省降水豐沛，但水資源分布不均，年內(nèi)各季節(jié)降水量相差較大，而且農(nóng)業(yè)用水消耗量大，用水效率較低[6-8]。農(nóng)業(yè)灌溉用水主要用于補給作物生長所需的水量，而作物需水量與土壤蒸發(fā)和作物蒸騰有著密切的關(guān)系，研究不同的農(nóng)業(yè)活動如何影響蒸散發(fā)至關(guān)重要[9]。

本研究通過長達兩年的對湖南省岳陽地區(qū)一年兩熟的蔬菜地進行高精度連續(xù)觀測，探索蒸散發(fā)變化規(guī)律和農(nóng)業(yè)活動對蒸散發(fā)的影響情況。

1 研究區(qū)概況

本研究始于2014年10月，止于2016年12月。實驗地點位于湖南省岳陽市君山區(qū)廣興洲鎮(zhèn)六支渠村的一塊典型施肥蔬菜農(nóng)田。實驗地屬典型的亞熱帶季風氣候，多年平均氣溫16.5～17.0℃，冬季平均氣溫3.8～4.7℃；夏季平均氣溫29～30℃左右，年平均降水量1 250～1 450 mm[10]；秋冬季盛行風向為東北風，春夏季盛行風向為南風；土壤質(zhì)地為沙壤土，離地表0～1 m的平均土壤容重為1.38 g/cm3，土壤含水量為29%，土壤偏堿性，pH值在7.0到8.0之間。

研究區(qū)域蔬菜地總面積大約25 hm2，中間有很多溝和水渠把整個蔬菜地分為面積1.5～2.5 hm2的小斑塊，這些斑塊東西向長100 m，南北向?qū)?～15 m 不等，每塊地分別由不同農(nóng)戶種植，作物類型、種植時間、施肥時間、施肥量都由農(nóng)戶自行決定。在兩年的觀測時間里，每年都會有兩個種植季，春季種植季始于每年的冬季，大約在2月份；夏季種植季始于每年的夏季，大約在8月份。春季種植季的作物一般有包菜、辣椒、南瓜、甜瓜和西瓜等，夏季種植季一般為包菜。每個種植季之前農(nóng)戶們都會施用肥料(一般為復合肥和碳酸氫銨)，施肥量大約在300～400 kg/hm2，待作物成熟收割之后，會進行旋耕滅茬。當種植作物為包菜時，種植密度為2 700株/hm2，每株之間間距為50 cm左右，成熟期包菜高40 cm 左右，種植深度為4 cm。灌溉采用引水渠漫灌的方式，夏季大約每10～15 d灌溉一次，平均每次灌溉水量10～15 mm。

2 研究方法

2.1 微氣象學方法—閉路QCLAS-EC法

閉路QCLAS-EC法是基于渦度相關(guān)法對通量進行監(jiān)測的方法。渦度相關(guān)是指某種物質(zhì)的垂直通量，即這種物質(zhì)的濃度與其垂直速度的脈動協(xié)方差。渦度相關(guān)法的優(yōu)點在于能通過測量各種屬性的湍流脈動值來得到它們的通量。和其他方法相比，渦度相關(guān)法并不是建立在經(jīng)驗關(guān)系基礎之上，或從其他氣象參量推論而來，而是建立在微氣象學的物理原理之上，是一種直接測量湍流通量的方法，比其他各種實測方法更加精密和可靠，并且該方法有利于蒸散發(fā)觀測工作的長期開展[10-13]。傳統(tǒng)渦度相關(guān)法測量蒸散通量多采用基于開路式氣體分析儀(如LI-7500，LI-COR Bioscience，美國)，優(yōu)點是省電、維護簡便，缺點是容易受到天氣影響，降雨期間的數(shù)據(jù)無法使用。本研究采用閉路QCLAS-EC法測量農(nóng)田的H2O通量，可以最大程度增加數(shù)據(jù)有效性。測量的通量是觀測塔上風向斑塊蔬菜地(約100 m范圍內(nèi)，根據(jù)通量貢獻區(qū)模型估算得到)的平均通量。

2.2 實驗儀器

選取靠近整個蔬菜地中間的一個地塊，在地塊東西向中點位置搭建安放渦度相關(guān)系統(tǒng)氣體分析儀的臨時小屋。整個系統(tǒng)的主要設備包括：一臺三維超聲風速儀(CSAT3，Campbell Scientific Inc，美國)、一臺數(shù)據(jù)采集器(CR3000，Campbell Scientific Inc，美國)、一臺大流量真空泵(XDS 35i，Edwards，美國)和一臺雙光路量子級聯(lián)可調(diào)諧紅外激光差分吸收光譜儀(QC-TILDAS-DUAL，Aerodyne Research Inc，美國)。數(shù)據(jù)采集器安裝在不銹鋼制作的觀測塔上，超聲風速儀以垂直于實驗地主風向的方向(東西向)垂直安裝在2m的高度。采用的QCLAS分析儀安裝了兩個不同波段(2 256.3 cm-1和1 901 cm-1)的激光發(fā)射器，可以高頻(10 Hz)同時測量多種溫室氣體及H2O濃度。采樣口位于超聲風速儀一側(cè)20 cm處，樣氣在真空泵的抽力下，經(jīng)過一條長約10 m的特氟龍管路(內(nèi)徑6.4 mm)進入QCLAS分析儀，完成H2O濃度的實時檢測。10 Hz的三維風速和H2O濃度數(shù)據(jù)存儲在CR3000數(shù)據(jù)采集器內(nèi)。

2.3 數(shù)據(jù)處理及通量計算

由于原始觀測數(shù)據(jù)有許多因為人為操作和儀器故障導致的無效值，所以必須進行數(shù)據(jù)預處理。數(shù)據(jù)的準確處理是渦動相關(guān)法能否正確測定蒸散量的關(guān)鍵。根據(jù)渦度相關(guān)系統(tǒng)得到的10 Hz湍流原始記錄，完整的通量計算過程如下：首先進行原始水汽濃度及風速數(shù)據(jù)的預處理，步驟如下：①查看超聲風速儀和激光分析儀的狀態(tài)信息報告，儀器出現(xiàn)異常狀態(tài)時要將此時間段的觀測數(shù)據(jù)刪除。②查看幾個觀測變量的原始數(shù)據(jù)，手動刪除明顯不合理的數(shù)據(jù)點(如明顯超出或者低于上下界閾值的數(shù)據(jù))，然后，進行幾項重要修正：①傾斜修正，即坐標旋轉(zhuǎn)，對于相對較為平坦均一的下墊面，一般采用2次坐標旋轉(zhuǎn)；而對于非平坦下墊面情況，采用“平面擬合”的方法。②對于每個測定的統(tǒng)計參數(shù)，如實測均值、實測方差、實測協(xié)方差，以及偏斜值、峰值等統(tǒng)計參數(shù)。依據(jù)數(shù)據(jù)獲取情況設置，必要時對水汽濃度做相應的校正或者轉(zhuǎn)換。③對切應力項，感熱通量以及潛熱通量和水汽通量等進行頻率響應修正。④校正各項通量，其中對顯熱通量做超聲虛溫修正，以用于去除濕度的影響。

3 結(jié)果與討論

3.1 蒸散量的季節(jié)變化趨勢

根據(jù)計算的0.5 h蒸散量，可得到研究區(qū)域日蒸散量(圖1)。從圖1中可以看出，蒸散量在2015年內(nèi)的日變化幅度較大，日蒸散量在0.4～6.7 mm之間，日與日之間變化無規(guī)律可循。初步猜測這與每日溫度、濕度、風速、降雨等氣象因子變化有關(guān)。

(a) 1月

(b) 2月

(c) 3月

(d) 4月

(e) 5月

(f) 6月

(g) 7月

(h) 8月

(i) 9月

(j) 10月

(k) 11月

(l) 12月

圖2為實驗地從2014年10月至2016年9月連續(xù)兩年的蒸散量實測數(shù)據(jù)。從圖2可以看出蒸散量在這兩年之內(nèi)具有相似的變化趨勢，變化規(guī)律大致可以總結(jié)為：蒸散量秋冬季較低，春夏季較高；蒸散量在種植季較高而在非種植季較低。具體來說，實測蒸散量從2014年10月至2015年1月逐漸下降，2015年1月蒸散量達到實驗期間較低水平，只有0.5 mm/d，蒸散量從2015年2月春季種植季開始逐漸升高，一直持續(xù)到2015年6月收割后蒸散量出現(xiàn)了明顯下降。造成這個現(xiàn)象的原因是每年的6月處于兩個種植季之間，土壤此時為裸土狀態(tài)，缺少了植物蒸騰導致。從2015年7月開始，隨著夏季的到來，氣溫和太陽輻射急劇升高，導致蒸散量也明顯增高；從2015年8月下半年種植季開始，蒸散量持續(xù)升高，在2015年8月達到全年最高水平，為3.5 mm/d；試驗地蒸散量從2015年9月開始持續(xù)降低，直到2016年1月的0.4 mm/d。造成蒸散量在冬季較低主要有兩個原因：①冬季溫度低，太陽輻射強度小于其他幾個季節(jié)；②冬季一般是非種植季，缺少了植物蒸騰。2016年的變化情況與2015年類似，不做贅述。經(jīng)計算，實驗地在2015年蒸散量為730 mm，每日平均蒸散量2.0 mm/d；2016年蒸散量為803 mm，每日平均蒸散量2.2 mm/d。

圖2 2014年10月至2016年9月實測蒸散量變化趨勢

3.2 蒸散量在植物生長季的影響因子分析

觀測采用高分辨率的蒸散監(jiān)測儀器，最小時間尺度可觀測0.5 h蒸散量。相比于其他渦動相關(guān)監(jiān)測方法，蒸散監(jiān)測儀的分辨率和精度都更高，并且可以在雨天進行觀測。本研究采用了兩個指標來評估農(nóng)業(yè)活動對蒸散的影響：①農(nóng)業(yè)活動發(fā)生前后0.5 h蒸散變化量；②農(nóng)業(yè)活動發(fā)生前后10 d累積蒸散量變化量。通過監(jiān)測到蒸散量的變化來定量分析不同農(nóng)業(yè)活動對蒸散的影響。

3.2.1 種植和收割對蒸散的影響量化評估

量化評估結(jié)果見圖3～5，表1～2。

(a) 2014年下半年

(b) 2015年上半年

(c) 2015年下半年

(d) 2016年上半年

(a) 2015年上半年

(b) 2015年下半年

(c) 2014年下半年

(d) 2016年上半年

時間2014年下半年種植季2015年上半年種植季峰值累計量峰值累計量種菜前0.4516.80.1214.7種菜后0.6419.60.1917.9時間2015年下半年種植季2016年上半年種植季峰值累計量峰值累計量種菜前0.2817.70.2113.9種菜后0.2919.20.2117.5

表2收菜前后實測蒸散量變化情況 mm

時間2014年下半年種植季2015年上半年種植季峰值累計量峰值累計量收菜前0.4624.20.4818.6收菜后0.0518.60.0510.5時間2015年下半年種植季2016年上半年種植季峰值累計量峰值累計量收菜前0.6023.60.3920.2收菜后0.4220.10.1612.3

圖5 每個種植季種植后和收割后蒸散量的變化

通過蒸散發(fā)實測數(shù)據(jù)，可以清楚看出種菜和收割對農(nóng)田實際蒸散的影響。其中，種植使農(nóng)田實測蒸散升高，種植后的幾天之內(nèi)，4個施肥季0.5 h蒸散峰值分別升高了0.19、0.07、0.01和0 mm；而10 d累計蒸散量分別升高了2.8、3.2、1.5和3.6 mm。種植使蒸散升高的主要原因是由于原本裸露的土壤表面突然加入了植物產(chǎn)生了作物蒸騰，但是由于作物處于生長初期階段，作物高度和生物量都較小，因此并不會使實測蒸散發(fā)生明顯升高，10 d累積量只升高了2～3 mm左右。而收割后的幾天之內(nèi)，實測蒸散則發(fā)生了明顯下降，4個施肥季0.5h蒸散峰值分別下降了0.41、0.43、0.26和0.23 mm；10 d累計蒸散量分別下降了5.8、8.1、3.5和7.9 mm。需要注意的是，收割使蒸散發(fā)生變化的幅度要明顯高于種植，這是由于收割時作物已經(jīng)處于生長后期，生物量較大，作物蒸騰已經(jīng)在總蒸散中占據(jù)主導地位，所以收割后蒸散量的下降幅度更大。

3.2.2 施肥與灌溉對蒸散的影響量化評估

施肥與灌溉對蒸散的影響量化評估見圖6～8、表3～4。

(a) 2014年下半年

(b) 2015年上半年

(c) 2015年下半年

(d) 2016年上半年

時間2014年下半年種植季2015年上半年種植季峰值累計量峰值累計量施肥前0.4219.20.01716.4施肥后0.4620.10.01315.1時間2015年下半年種植季2016年上半年種植季峰值累計量峰值累計量施肥前0.5620.50.1314.9施肥后0.4519.90.1516.5

(a) 2014年下半年

(b) 2015年上半年

(c) 2015年下半年

(d) 2016年上半年

時間2014年下半年種植季2015年上半年種植季峰值累計量峰值累計量灌溉前0.0719.80.0815.9灌溉后0.1421.50.2716.6時間2015年下半年種植季2016年上半年種植季峰值累計量峰值累計量灌溉前0.2420.20.0116.3灌溉后0.6322.80.4416.9

圖8 每個種植季施肥和灌溉后蒸散量的變化

4個種植季施肥后的幾天之內(nèi)，蒸散量雖然發(fā)生了變化，但是這種變化并沒有明顯的規(guī)律。2014年下半年種植季和2016年上半年種植季，施肥后的10 d累計蒸散量分別上升了0.9 mm和1.6 mm，而在另外兩個種植季，2015年上半年種植季和2015年下半年種植季，施肥后的10 d累計蒸散量分別下降了1.3 mm和0.6 mm。因此，農(nóng)田施肥對于實際蒸散是促進還是抑制作用從本次試驗中并不能得出定性的結(jié)論，這是由于每次施肥對土壤注入了氮磷鉀等元素，并沒有直接的水分加入。但是灌溉對蒸散的影響卻比較顯著，從圖8可以看出，4個種植季每次灌溉之后0.5 h蒸散峰值分別升高了0.07、0.19、0.39和0.43 mm，灌溉前后10 d累計蒸散量分別升高了1.7、0.7、2.3和0.6 mm。灌溉對農(nóng)田實際蒸散的促進作用主要是由于兩個原因?qū)е拢孩倜看喂喔戎蠖紩鹜寥篮康纳?，而土壤蒸發(fā)的來源是土壤里的水分，因此灌溉會增大土壤蒸發(fā)量。②植物生長是一個需水過程，灌溉水的注入會促進作物生長，進而加大了植物蒸騰量。

4 結(jié) 論

蒸散量與農(nóng)業(yè)種植規(guī)律有很強的相關(guān)性，同時也表現(xiàn)出很強的季節(jié)差異性：實際蒸散量在1月份較小，從2月份種植季開始逐漸增大，到5月份達到最高；6月份由于蔬菜收割，農(nóng)地變成裸地，蒸散量顯著減小；下半年種植季從9月份開始蒸散量逐漸增高，10月份達到最大，11月份開始逐漸減小?？傮w來說，蒸散量在春夏季較高，秋冬季較低；在種植季較高，非種植季較低。此外，人類的種植活動增加了田間蒸散量，而收割則使蒸散量發(fā)生明顯下降；灌溉對總蒸散有促進作用，但是這種促進作用幅度比較小，而施肥對于蒸散的影響卻并不明顯。

[1]宋小寧.基于植被蒸散法的區(qū)域缺水遙感監(jiān)測方法研究[D].北京：中國科學院遙感應用研究所,2004.

[2]辛曉洲,田國良,柳欽火.地表蒸散定量遙感的研究進展[J].遙感學報,2003,7(3):233-240.(XIN Xiaozhou,TIAN Guoliang,LIU Qinhuo.A review of researches on remote sensing of land surface evapotranspiration[J].Journal of Remote Sensing,2003,7(3):233-240.(in Chinese))

[3] 趙傳燕.甘肅省祖厲河流域潛在生態(tài)條件的GIS輔助模擬研究[D].蘭州：蘭州大學,2003.

[4] 童瑞,楊肖麗,任立良,等. 黃河流域1961—2012年蒸散發(fā)時空變化特征及影響因素分析[J]. 水資源保護，2015,31(3):16-21.(TONG Rui,YANG Xiaoli,REN Liliang，et al. Temporal and spatial characteristics of evapotranspiration in the Yellow River Basin during 1961—2012 and analysis of its influence factors[J].Water Resources Protection,2015,31(3):16-21.(in Chinese))

[5] 張守平,蒲強,李麗琴，等.基于可控蒸散發(fā)的狹義水資源配置[J].水資源保護，2012,28(5):13-18.(ZHANG Shouping,PU Qiang,LI Liqin，et al.Special water resources allocation based on controllable ET[J].Water Resources Protection,2012,28(5):13-18.(in Chinese))

[6]鄒君.湖南省水資源可持續(xù)利用綜合評價研究[J].節(jié)水灌溉,2007(2):18-21.(ZOU Jun.Study on comprehensive evaluation of water resources in Hunan Province[J].Water Saving Irrigation,2007(2):18-21.(in Chinese))

[7]湖南省國土委員會編.湖南國土資源[M].長沙:湖南科技出版社,1985.

[8]劉培亮,毛德華,周慧,等.1990-2013年湖南四水入洞庭湖汛期徑流量的變化規(guī)律[J].水資源保護,2015,31(4):52-61.(LIU Peiliang,MAO Dehua,ZHOU Hui,et al.Variation law of runoff in flood seasons into Dongting Lake from Four Rivers in Hunan Province during 1990-2013[J].Water Resources Protection,2015,31(4):52-61.(in Chinese))

[9]張強,孫鵬,程辰,等.ENSO影響下安徽省旱澇災害及農(nóng)業(yè)生產(chǎn)損失時空變化特征[J].水資源保護,2016,32(6):6-18.(ZHANG Qiang,SUN Peng,CHENG Chen,et al.Temporal and spatial variations of ENSO-induced flood and drought disasters and loss in agricultural production in Anhui Province[J].Water Resources Protection,2016,32(6):6-18.(in Chinese))

[10] 侯志勇,謝永宏,陳心勝,等.洞庭湖濕地的外來入侵植物研究[J].農(nóng)業(yè)現(xiàn)代化研究,2011,32(6):744-747.(HOU Zhiyong,XIE Yonghong,CHEN Xinsheng,et al.Study on invasive plants in dongting lake wetlands[J].Research of Agricultural Modernization,2011,32(6):744-747.(in Chinese))

[11]MASSMAN W J.A simple method for estimating frequency response corrections for eddy covariance systems[J].Agricultural amp; Forest Meteorology,2000,104(3):185-198.

[12]SWINBANK W C.The measurement of vertical transfer of heat and water vapor by eddies in the lower atmosphere [J].Journal of Meteor,1951,8:135-145.

[13]SUN X M,ZHU Z L,WEN X F,et al.The impact of averaging period on eddy fluxes observed at China FLUX sites[J].Agricultural amp; Forest Meteorology,2006,137(3):188-193.

Studyonin-situevapotranspirationmeasurementanditsinfluentialfactorsinfarmlandinSouthernChina

CHENHan1,HUANGJinhui2

(1.HehaiCollegeofChongqingJiaotongUniversity,Chongqing400074,China; 2.CollegeofEnvironmentalScienceandEngineering,NankaiUniversity,Tianjin300350,China)

Field observations for two consecutive years were carried out on the actual evapotranspiration, water vapor flux and latent heat flux (LE) for two consecutive years in the vegetable land in the suburb of Yueyang City, Hunan Province by using eddy covariance technique based on closed-path QCLAS-EC Analyzer. The original observation data were processed and the evapotranspiration was calculated so as to analyze the annual periodic change trend and seasonal variation law of the regional evapotranspiration. The results show that annual evapotranspiration in the study area were between 730 and 803 mm; the seasonal variation was substantial, reaching the highest level in August, 3.5 mm/d and in January, only 0.4mm/d. The actual evapotranspiration in the whole study area can be summarized as follows: higher in spring and summer, lower in autumn and winter; higher in planting season, lower in non-planting season; higher at midday and afternoon, lower at night. In addition, the effects of different agricultural activities on the actual evapotranspiration were observed by using high-resolution measured data. The findings are that the evapotranspiration was increased by planting and was decreased by harvesting, and the reduction of evapotranspiration by harvesting was significantly higher than that by planting. While irrigation also contributed to the promotion of evapotranspiration but was less pronounced, the effect of final fertilization on evapotranspiration was not yet clear.

evapotranspiration; southern farmland; vegetable land; eddy covariance method

10.3880/j.issn.1004-6933.2017.06.12

國家自然科學基金(41561124015)

陳晗(1992—)男，碩士研究生，研究方向為生態(tài)水文。E-mail:15620691073@163.com

黃津輝，教授。E-mail:huangj@nankai.edu.cn

P332.2

A

1004-6933(2017)06-0079-08

2017-04-19 編輯：彭桃英)