黃土塬區麥田能量平衡特征

苑 旭，朱元駿,，劉文兆,，韓曉陽

(1.西北農林科技大學水土保持研究所，黃土高原土壤侵蝕與旱地農業國家重點實驗室， 陜西 楊凌 712100；2.中國科學院水利部水土保持研究所， 陜西 楊凌 712100)

黃土塬區麥田能量平衡特征

苑 旭1，朱元駿1,2，劉文兆1,2，韓曉陽2

(1.西北農林科技大學水土保持研究所，黃土高原土壤侵蝕與旱地農業國家重點實驗室， 陜西 楊凌 712100；2.中國科學院水利部水土保持研究所， 陜西 楊凌 712100)

以黃土塬區冬小麥田為研究對象，基于渦度相關數據分析麥田能量平衡的日變化、季節變化和能量分配特征及其主控因子。結果表明，長武塬區麥田全年獲得的凈輻射(Rn)為2.56×103MJ·m-2·a-1，渦度相關系統的能量閉合度達到0.72。冬小麥生育期內，越冬期和灌漿期麥田主要能量支配項為感熱通量(H)，最大值出現在6月，為7.09 MJ·m-2·d-1；其他生育期和休閑期，主要能量支配項為潛熱通量(LE)，最大值出現在5月，為10.71 MJ·m-2·d-1。波文比(β)在生育期平均值為0.57，休閑期為0.46。土壤熱通量(G)年總量為-15.26 MJ·m-2·a-1，日總量最大值出現在6月，為1.85 MJ·m-2·d-1，10月至次年1月為負值，表現土壤釋放熱量。

冬小麥；渦度相關法；能量平衡；潛熱通量；感熱通量；黃土高原

黃土塬區是黃土高原重要的糧食產區，冬小麥是該區的主要作物，區域太陽輻射資源豐富而水資源相對短缺。然而，關于黃土塬區農田能量平衡方面的研究報道還比較少。本文利用麥田渦度相關系統[9-10]觀測數據，并結合同步微氣象資料進行系統分析，以期闡明麥田能量平衡分量的時間變化、分配特征及其主控因子，為進一步評價區域能量與水量平衡狀況及特征奠定基礎。

1 研究區概況

觀測樣地設在中國科學院長武黃土高原農業生態試驗站水熱通量觀測場，地處黃土高原南部陜甘交界處的陜西省長武縣洪家鎮王東村(北緯35°14′，東經107°4l′，海拔1 220 m)。該地區屬于暖溫帶半濕潤大陸性季風氣候，年均降水584 mm，年均氣溫9.1℃，無霜期171 d，地下水埋深50～80 m，是黃土高原典型旱作農業區。觀測場位于黃土塬上，土壤類型為黑壚土，質地均勻疏松。觀測場內種植冬小麥，其生長季為每年10月至次年6月，地面平整，構成大范圍均一的下墊面，滿足實驗要求盛行風向的風浪區長度[11]。

2 材料與方法

2.1 觀測方法

渦度相關系統安裝于1.86 m高度的伸展臂上，由三維超聲風速儀(1210R3，Gill Instruments Ltd.，UK)和開路CO2/H2O分析儀(Li-7500，Li-Cor Ltd.，USA)組成，原始采樣頻率為10 Hz，數據傳輸由數據采集器(CR5000，Campbell Scientific Ltd.，USA)進行存儲。冬小麥冠層上方還采用凈輻射傳感器(CNR-1，Kipp & Zonen，Netherlands)測量凈輻射，空氣溫濕儀(HMP-45C，Vaisala，Finland)測量空氣溫、濕梯度。將2個土壤熱通量板(HFP01，Hukse Flux，Netherlands)安置于土壤深5cm處用于測量土壤熱通量，同時進行5層土壤溫度(地面以下2、10、20、40、80 cm，105-T，Campbell Scientific Ltd.，USA)和5層土壤濕度(地面以下2、10、20、40、80 cm，CS616，Campbell Scientific Ltd.，USA)的測定，這些要素均按30 min計算平均值進行存儲。其他氣象因素由常規氣象觀測系統測定。觀測時間為2008年6月1日至2009年5月31日。

2.2 數據處理

2.2.1 渦度相關數據處理 近地面層是大氣邊界層的最底層，在該氣層內空氣運動符合湍流交換規律，可以利用渦度相關法測定近地層大氣中熱量和水汽的垂直輸送通量。當下墊面均勻一致時，在近地層內潛熱通量、感熱通量計算公式如下[12]：

(1)

(2)

采用Edire軟件處理10Hz的原始數據，剔除野點，對風速進行二次坐標旋轉，計算30min的水熱通量均值，并對其進行頻率響應校正、AOA(Angleofattack)校正[13-14]、密度脈動和濕度校正[15]。受儀器故障、天氣等影響，渦度相關數據存在異常值或缺測[16]。缺失數據時間間隔小于或等于2h的采用線性內插；大于2h的采用平均晝夜變化法[17-18]。

2.2.2 能量閉合計算 能量閉合狀況可以用來檢驗渦度相關系統觀測結果的可靠性，即器測的湍流通量與可利用能量的平衡情況，方程為：

LE+H=Rn-G-S-Q

(3)

式中，Rn為凈輻射(W·m-2)，G為土壤熱通量(W·m-2)；S為地表到渦度相關系統觀測高度間的空氣與作物冠層儲熱(W·m-2)；Q為其他能量源和匯(W·m-2)。對裸地或植被低矮的生態系統，S和Q的數值通常很小，可忽略不計[19-20]。土壤熱通量G包括埋深5 cm的傳感器所測熱通量G0(W·m-2)和其上層土壤儲熱Gs(W·m-2)兩部分，其中Gs采用的計算公式如下[21-23]：

Gs=Zp×(ρwθCw+ρsCs)×

(4)

式中，Zp是土壤熱通量板安裝深度(5 cm)；ρw是水的密度(kg·m-3)；θ是土壤體積含水率(m3·m-3)；Cw是水的比熱容(J·kg-1·K-1)；ρs是土壤容重(1 300 kg·m-3)；Cs是干土比熱容(1 620 J·kg-1·K-1)；T1是2 cm土壤溫度(K)；ΔT是2 cm與Zp的土壤溫度差(K)；Δt是測定溫度的時間間隔(30 min)；5 cm土壤溫濕度由2 cm和10 cm土壤溫濕度線性內插得到。

本文采用湍流通量(LE+H)和有效能量(Rn-G)[24-25]的線性回歸關系來評價能量閉合狀況。理想情況下斜率為1，截距為0，但實際情況(LE+H)往往比(Rn-G)小[19,24,26-27]。

3 結果與分析

3.1 能量閉合度分析

對湍流通量與有效能量全年30min平均值進行閉合度分析(圖1a)，回歸直線斜率為0.65，截距是9.19W·m-2。如果強制通過原點，斜率增加到0.68。若考慮土壤熱儲量后再進行能量閉合分析(圖1b)，回歸直線斜率提高到0.71，截距減小為4.24W·m-2。與上一次強制通過原點相比，斜率變為0.72，增加了7%，但與理想曲線(斜率為1)仍相差較大。Wilson等[19]對全球通量網(FLUXNET，http://fluxnet.ornl.gov/)站點的觀測數據分析表明斜率變化范圍是0.53～0.99，平均值為0.79，截距變化范圍是-32.9～36.9W·m-2，平均值為3.7W·m-2。Li等[26]對中國通量網(ChinaFLUX，http://www.chinaflux.org/)站點的觀測數據分析指出斜率變化范圍是0.49～0.81，平均值為0.67，截距變化范圍是10.8～79.9W·m-2，平均值為28.9W·m-2。長武站通量觀測場能量閉合度接近上述報道的平均值。

注:a.未考慮土壤熱儲量的能量閉合;b.考慮土壤熱儲量的能量閉合Note:a.soilheatstoragewasnotconsideredforenergybalanceclosure;b.soilheatstoragewasconsideredforenergybalanceclosure.

圖1 麥田能量閉合狀況

Fig.1Energybalanceclosureinwinterwheatfield

3.2 能量平衡分量的日變化

麥田能量平衡分量在大部分生育階段均有顯著的日變化趨勢。以圖2所示各生育階段平均日變化為例，凈輻射、潛熱通量和感熱通量都呈早晚低、中午高的“單峰型”曲線特征，日變化峰值出現的時間在不同生育階段并無明顯差異，均在14∶00前后。其他地區農田研究[5-6]也發現潛熱通量、感熱通量和凈輻射日變化趨勢相似。凈輻射在日出后由負值逐漸轉為正值，即轉變為農田能量的收入項，午間達到最大值，其中抽穗期日變化峰值最大，為494.70W·m-2，而越冬期日變化峰值最小，為246.26W·m-2，日落前由正值逐漸轉為負值。一天中剩余大部分時間表現為長波輻射能的損失。

圖2 麥田能量平衡分量日變化

Fig.2Diurnalvariationoftheenergybalancecomponentsinwinterwheatfield

潛熱通量同樣在抽穗期日變化峰值最大，為234.83W·m-2，而越冬期日變化峰值最小，為32.79W·m-2。在夜間空氣濕度較高和風速較低的情況下，潛熱通量有時會表現為負值[30]。感熱通量在灌漿期日變化峰值最大，為201.62W·m-2，而出苗分蘗期日變化峰值最小，為52.58W·m-2。夜間出現逆溫層可使感熱通量為負值。白天潛熱和感熱通量遠遠大于夜間的潛熱和感熱通量，說明潛熱和感熱輸送主要受到太陽輻射強度的制約，與賈志軍等[6]對三江平原大豆田水熱通量的研究結論一致。

土壤熱通量在越冬期日變化平緩，而在其他生育階段都呈“單峰型”日變化。在白天，土壤熱通量表現為向下傳輸，即吸熱狀態。在夜間，有數小時表現為向上傳輸，即放熱狀態。吸熱與放熱狀態時間分配主要受到土壤溫度和植被生長狀況的影響。總之，農田氣候季節變化和小麥生長狀況是影響能量平衡分量日變化的主要因素。

3.3 能量平衡分量的季節變化

3.3.1 不同月份能量平衡分量的變化特征 凈輻射日總量在年內呈近似單峰變化(圖3)，其中雨季凈輻射的日間差異較大。全年凈輻射總量為2.56×103MJ·m-2·a-1，日總量峰值出現在7月份，為18.59MJ·m-2·d-1，但由于該月降水集中(降水天數17d，降水量134.8mm)，月平均值(10.90MJ·m-2·d-1)小于6月份(11.26MJ·m-2·d-1)。12月份平均值最小，為2.80MJ·m-2·d-1，分別小于1月份的3.40MJ·m-2·d-1與11月份的3.36MJ·m-2·d-1。土壤熱通量年總量為-15.26MJ·m-2·a-1，日總量最大值(1.85MJ·m-2·d-1)出現在6月。從10月至次年1月，土壤熱通量日總量始終為負值，即農田土壤釋放能量。說明在該時段對整個系統而言，土壤是熱源，能量通過土壤向大氣傳送。

圖3 麥田能量平衡分量季節變化

Fig.3Seasonalvariationoftheenergybalancecomponentsinwinterwheatfield

感熱通量日總量從4月中旬開始隨凈輻射增加而增大，最大值(7.09MJ·m-2·d-1)出現在凈輻射最強的6月份，隨后逐漸降低。11月又開始出現略微增加趨勢，受凈輻射持續降低的影響，最后趨于平穩。潛熱通量日總量從2月開始明顯增大，最大值(10.71MJ·m-2·d-1)出現在小麥生長旺盛的5月份，隨著小麥生長后期生理活動減弱，呈顯著降低趨勢。6月中旬以后，潛熱通量日總量又開始增大，至7月中旬達到第二個峰值，為9.08MJ·m-2·d-1。許多農田研究[5-6,31]中發現降水量是影響潛熱通量的重要因子，大氣降水是黃土塬區麥田水分的主要來源，降水增加導致土壤含水量升高，有助于作物根系吸水并促進土壤蒸發。因此，7月份麥田為裸地時潛熱通量日總量較高可能與降水量(7月份降水量比6月份增加51.4mm)增大、土壤含水量增高有關。

3月中旬，LE/Rn隨著小麥生長開始明顯增加(圖4)，最大值出現在5月份，而H/Rn變化平穩。小麥進入到生長后期，H/Rn超過LE/Rn，成為能量的主要支配項。隨后LE/Rn和H/Rn都進入到平穩變化階段，至12月中旬，H/Rn超過LE/Rn，成為能量的主要支配項。從2月中旬到9月，G/Rn一直呈平穩變化，隨著土壤溫度降低，G/Rn開始迅速下降，在冬季土壤熱通量大部分時間都是向系統釋放熱量。從全年來看，感熱通量占凈輻射25.5%，而潛熱通量占51.3%，這一比例小于Lei等[8]在華北平原冬小麥-玉米輪作灌溉農田的觀測結果(59%)和Suyker等[32]在美國中北部平原玉米-大豆輪作灌溉農田的觀測結果(60%)，原因主要在于前面說到的渦度相關系統能量不閉合的因素影響。

3.3.2 不同生育階段能量平衡分量的變化特征 表1給出麥田各生育階段能量平衡分量和主要環境因子的變化情況。一年中各生育階段數值差異較大。在越冬期潛熱值出現最低，平均值為0.84 MJ·m-2·d-1，大部分能量轉化成感熱，平均值為1.32 MJ·m-2·d-1，土壤熱通量值也較小，平均值為-0.37 MJ·m-2·d-1。波文比(β=H/LE，即感熱通量H與潛熱通量LE之比)在出苗分蘗期小于1，潛熱通量成為能量的主要支配項，越冬期β值增大到1.58，感熱通量變成能量的主要支配項。進入返青期后，隨著作物生長，β值減小至1以下，潛熱通量成為能量的主要支配項，直到抽穗期。灌漿期β值又增大至1.05，感熱通量超過潛熱通量，這可能與小麥生長后期蒸騰作用下降，潛熱通量相對減少有關。但成熟期后β值又開始小于1，并逐漸減小。整個生育期內波文比平均值為0.57，高于Hernandez-Ramirez等[3]在美國中西部玉米田的觀測結果(0.39)和大豆田的觀測結果(0.30)。從表1中可以看出，除了受農田物候變化的影響，氣候條件和土壤狀況也是β值季節變化的重要影響因子。

圖4 麥田能量分配比率季節變化

Fig.4 Seasonal variation of the ratio of energy partitioning in winter wheat field

表1 不同生育階段能量平衡分量與主要環境因子的平均值變化情況

注：Ts為5 cm土壤溫度；SWC為5 cm土壤體積含水量；Ta為空氣溫度；P為降水量。

Note:Tsis soil temperature at 5 cm depth;SWCis soil water content at 5 cm depth;Tais air temperature;Pis precipitation.

4 結 論

1) 潛熱通量和感熱通量之和比較于凈輻射和土壤熱通量之差，所給出的長武站麥田渦度相關系統能量閉合度為0.72。凈輻射、潛熱通量和感熱通量在不同生育階段日變化均呈“單峰型”， 日變化峰值出現的時間在不同生育階段并無明顯差異。受到太陽輻射強度與作物生理活動變化的制約，白天潛熱和感熱通量遠遠大于夜間的潛熱和感熱通量。土壤熱通量在越冬期日變化平緩，而在其他生育階段日變化模式與其他分量基本相同。農田氣候季節變化和冬小麥生長狀況是影響能量平衡分量日變化的主要因素。

2) 麥田凈輻射全年呈單峰變化，雨季凈輻射的日間差異較大主要受降水的影響；感熱通量和潛熱通量季節變化主要受凈輻射的影響，同時冬小麥生長狀況和降水也是影響潛熱通量季節變化的主要因素；土壤熱通量從10月至次年1月主要表現為土壤釋放熱量，其他時間主要表現為吸收熱量，其變化過程主要受土壤溫度梯度及植被覆蓋狀況影響。

3) 麥田全年主要以潛熱的形式進行能量交換。冬小麥生育期內，β值在越冬期平均為1.58，灌漿期略大于1，感熱通量為能量主要支配項；其他生育期和休閑期小于1，潛熱通量為能量主要支配項。除了受農田物候變化的影響，氣候條件和土壤狀況也是β值季節變化的重要影響因子。

致謝：本研究得到了中科院寒區旱區環境與工程研究所王介民研究員的指導，在此對王老師的幫助表示誠摯的感謝。

[1] Baldocchi D. A comparative study of mass and energy exchange over a closed C3(wheat) and an open C4(corn) canopy: I. The partitioning of available energy into latent and sensible heat exchange[J]. Agricultural and Forest Meteorology, 1994,67(3):191-220．

[2] Anderson D E, Verma S B, Rosenberg N J. Eddy correlation measurements of CO2, latent heat and sensible heat fluxes over a crop surface[J]. Boundary Layer Meteorology, 1984,29(3):263-272.

[3] Hernandez-Ramirez G, Hatfield J L, Prueger J H, et al. Energy balance and turbulent flux partitioning in a corn-soybean rotation in the Midwestern US[J]. Theoretical and Applied Climatology, 2010,100(1-2):79-92.

[4] Suyker A E, Verma S B. Evapotranspiration of irrigated and rainfed maize-soybean cropping systems[J]. Agricultural and Forest Meteorology, 2009,149(3-4):443-452.

[6] 賈志軍,黃 耀.三江平原大豆田湍流熱通量研究[J].中國農業科學,2012,45(7):1274-1284.

[7] 郭家選,梅旭榮,林 琪,等.冬小麥農田暫時水分脅迫狀況下水、熱通量日變化[J].生態學報,2006,26(1):130-137.

[8] Lei H M, Yang D W. Interannual and seasonal variability in evapotranspiration and energy partitioning over an irrigated cropland in the North China Plain[J]. Agricultural and Forest Meteorology, 2010,150(4):581-589.

[9] 郭家選,梅旭榮,盧志光,等.測定農田蒸散的渦度相關技術[J].中國農業科學,2004,37(8):1172-1176.

[10] Liu S M, Xu Z W, Zhu Z L, et al. Measurements of evapotranspiration from eddy-covariance systems and large aperture scintillometers in the Hai River Basin, China[J]. Journal of Hydrology, 2013,487(9):24-38．

[11] 李雙江,劉文兆,高橋厚裕,等.黃土塬區麥田CO2通量季節變化[J].生態學報,2007,27(5):1987-1992.

[12] Baldocchi D D, Hincks B B, Meyers T P. Measuring biosphere-atmosphere exchanges of biologically related gases with micrometeorological methods[J]. Ecology, 1988,69(5):1331-1340.

[13] Nakai T, Shimoyama K. Ultrasonic anemometer angle of attack errors under turbulent conditions[J]. Agricultural and Forest Meteorology, 2012,162-163(3):14-26．

[14] Xu Z W, Liu S M, Li X, et al. Intercomparison of surface energy flux measurement systems used during the HiWATER-MUSOEXE[J]. Journal of Geophysical Research-atmospheres, 2013,118(23):13140-13157.

[15] 徐自為,劉紹民,宮麗娟,等.渦動相關儀觀測數據的處理與質量評價研究[J].地球科學進展,2008,23(4):357-370.

[16] Foken T, Wichura B. Tools for quality assessment of surface-based flux measurements[J]. Agricultural and Forest Meteorology, 1996,78(1):83-105.

[17] Falge E, Baldocchi D, Olson R, et al. Gap filling strategies for defensible annual sums of net ecosystem exchange[J]. Agricultural and Forest Meteorology, 2001,107(1):43-69.

[18] 徐自為,劉紹民,徐同仁,等.渦動相關儀觀測蒸散量的插補方法比較[J].地球科學進展,2009,24(4):372-382.

[19] Wilson K, Goldstein A, Falge E, et al. Energy balance closure at FLUXNET sites[J]. Agricultural and Forest Meteorology, 2002,113(1-4):223-243.

[20] Chávez J L, Howell T A, Copeland K S. Evaluating eddy covariance cotton ET measurements in an advective environment with large weighing lysimeters[J]. Irrigation Science, 2009,28(1):35-50.

[21] 陽 坤,王介民.一種基于土壤溫濕資料計算地表土壤熱通量的溫度預報校正法[J].中國科學(D輯:地球科學),2008,38(2):243-250.

[22] Oncley S P, Foken T, Vogt R, et al. The Energy Balance Experiment EBEX-2000. Part I: overview and energy balance[J]. Boundary-Layer Meteorology, 2007,123(1):1-28.

[23] Liebethal C, Foken T. Evaluation of six parameterization approaches for the ground heat flux[J]. Theoretical and Applied Climatology, 2007,88(1):43-56.

[24] Foken T. The energy balance closure problem: an overview[J]. Ecological Applications, 2008,18(6):1351-1367.

[25] 黃 輝,孟 平,張勁松,等.華北低丘山地人工林蒸散的季節變化及環境影響要素[J].生態學報,2011,31(13):3569-3580.

[26] Li Z Q, Yu G R, Wen X F, et al. Energy balance closure at ChinaFLUX sites[J]. Science in China(Series D:Earth Sciences), 2005,48(S1):51-62.

[27] Liu S M, Xu Z W, Wang W Z, et al. A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem[J]. Hydrology and Earth System Sciences, 2011,15(4):1291-1306.

[28] Foken T, Mauder M, Liebethal C, et al. Energy balance closure for the LITFASS-2003 experiment[J]. Theoretical and Applied Climatology, 2010,101(1):149-160.

[29] Foken T, Aubinet M, Finnigan J J, et al. Results of a panel discussion about the energy balance closure correction for trace gases[J]. Bulletin of the American Meteorological Society, 2011,92(4):ES13-ES18.

[30] Liu W T. Moisture and latent heat flux variabilities in the tropical Pacific derived from satellite data[J]. Journal of Geophysical Research, 1988,93(C6):6749-6760.

[31] 陽伏林,張 強,王文玉,等.黃土高原春小麥農田蒸散及其影響因素[J].生態學報,2014,34(9):2323-2328．

[32] Suyker A E, Verma S B. Interannual water vapor and energy exchange in an irrigated maize-based agroecosystem[J]. Agricultural and Forest Meteorology, 2008,148(3):417-427.

Characteristics of energy balance in winter wheat field on the Loess Tableland

YUAN Xu1, ZHU Yuan-jun1,2, LIU Wen-zhao1,2, HAN Xiao-yang2

(1.StateKeyLaboratoryofSoilErosionandDrylandFarmingontheLoessPlateau,InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China；2.InstituteofSoilandWaterConservation,ChineseAcademyofSciencesandMinistryofWaterResources,Yangling,Shaanxi712100,China)

In this paper, eddy covariance data from a winter wheat field on the Loess Tableland were used to investigate the temporal variation of energy components, as well as the characteristics of energy partitioning and its main controlling factors. Results showed that the total net radiation (Rn) throughout the year was 2.56×103MJ·m-2·a-1, and the energy balance closure degree was 0.72. In the wintering and grain filling periods, the energy components were dominated by upward sensible heat flux (H), with a maximum value of 7.09 MJ·m-2·d-1occurring in June. However, in other growing stages and the fallow period, the energy components were dominated by upward latent heat flux (LE), with a maximum value of 10.71 MJ·m-2·d-1occurring in May. The average of Bowen ratio (β) was 0.57 in growing season and 0.46 in the fallow period. The annual total amount of soil heat flux (G) was -15.26 MJ·m-2·a-1. The daily total amount of soil heat flux (G) reached a maximum value of 1.85 MJ·m-2·d-1in June, and it was negative from October to the next January, which showed the release of heat from soil.

winter wheat; eddy covariance method; energy balance; latent heat flux; sensible heat flux; Loess Plateau

1000-7601(2017)02-0227-07

10.7606/j.issn.1000-7601.2017.02.37

2016-01-11基金項目：國家自然科學基金項目(41571036)；公益性行業(氣象)科研專項(重大專項)(GYHY201506001)

苑 旭(1988—)，男，吉林松原人，碩士，主要從事流域管理研究。E-mail:yuanxu833@126.com。

劉文兆(1960—)，男，陜西乾縣人，研究員，博士生導師，主要從事水文生態與流域管理研究。 E-mail:wzliu@ms.iswc.ac.cn。

P422.4

A