三江平原寒地稻田蒸散量估算研究

摘要：基于非稱重式蒸滲儀測定的寒地稻田實際蒸散量，探討波文比能量平衡法和Penman-Monteith參考作物蒸散-作物系數法（PM-Kc模型）在寒地稻田蒸散量計算中的適用性。結果表明，波文比能量平衡法估算的蒸散量較實測蒸散量偏低9.1%，但總體有較好的相關性（R2=0.87）；校正后的PM-Kc模型低估了水稻生長中期和后期的累積ET，但其相對誤差均小于10%，而生長初期，PM-Kc模型得到的累積ET比實測的累積ETa高19.34%，以生長初期為例進一步分析PM-Kc模型與實際蒸散ETa的關系，展現了隨著間隔日數的增加，RMSE、MAE有降低的趨勢，一致性系數d有提高的趨勢，PM-Kc模型得到的7日平均模擬實際蒸散與實測蒸散數據最為接近。

關鍵詞：寒地稻田；蒸散；估算；波文比能量平衡法；Penman-Monteith參考作物蒸散-作物系數法

中圖分類號：S161.4 文獻標識碼：A 文章編號：0439-8114（2016）10-2525-04

DOI：10.14088/j.cnki.issn0439-8114.2016.10.019

Abstract： Based on the actual evapotranspiration of cold rice measured by the non-weighing lysimeter， the applicability of the Bowen ratio energy balance method and the Penman-Monteith reference evapotranspiration-crop coefficient method（PM-Kc model） in the estimation of the evapotranspiration of cold rice was discussed. The results showed that the evapotranspiration estimated by the Bowen ratio energy balance method was 9.1% lower than the measured actual evapotranspiration， but overall there was a good correlation between them， with the R2 of 0.87. The corrected PM-Kc model underestimated the cumulative ET of cold rice in the medium-term growth and the late-term growth stage， but the relative error was less than 10%， while the cumulative ET of cold rice in the stage of early growth estimated by the PM-Kc model was 19.34% higher than the measured actual cumulative ET. It was shown that RMSE and MAE decreased while the consistency coefficient increased with the increase of the number of interval days， through the further analysis of the relationship between the estimated evapotranspiration and the measured actual evapotranspiration of cold rice in the stage of early growth. And the 7 days’ average evapotranspiration estimated by PM-Kc model was closest to the measured actual evapotranspiration.

Key words： cold rice fields； evapotranspiration； estimation； bowen ratio energy balance method； Penman-Monteith reference evapotranspirationcrop coefficient model

作物的農田蒸散是農業用水的主要組成部分，是一定時段內作物蒸騰和土壤棵間蒸發之和，它與氣象因素、土壤的供水狀況以及作物種類和長勢有關[1]。準確估算和模擬植被蒸散可為提高水分利用效率、合理配置水資源及生態系統可持續經營管理提供科學依據[2-5]。蒸散觀測和蒸散量的估算也是農業節水研究的基礎，農業節水措施的制定關鍵是需要準確測定和估算農田的實際蒸散量。三江平原是中國的9大商品糧基地之一，是中國典型的灌溉農業區，如何有效地利用有限的水資源，滿足農作物正常需水是當地農業面臨的重大問題[6-8]。本研究以三江平原寒地稻田小區試驗資料為基礎，采用波文比能量平衡法和Penman-Monteith參考作物蒸散-作物系數法對三江平原寒地稻田蒸散量進行了估算，并與非稱重式蒸滲儀測定的稻田蒸散量進行比較分析，以期尋找適當的方法對三江平原寒地稻田蒸散量進行估算，為寒地稻田的適時灌溉及灌溉的管理和規劃設計提供科學依據。

1 試驗區概況與研究方法

1.1 試驗區概況與資料來源

試驗在中國科學院三江平原沼澤濕地生態試驗站農田綜合試驗場進行。該試驗場位于三江平原腹地。站區中心地理坐標為47.35°N，133.31°E，海拔高度55.6 m，具有溫帶大陸性季風氣候，年平均氣溫為1.9 ℃，7月平均氣溫為22.0 ℃，全年最高氣溫為37.0 ℃，最低溫度為-38.0 ℃，無霜期為125 d左右，年降水量為600 mm左右，且主要集中在7～9月，占全年降水量的60%[9]。供試水稻田是1995年由沼澤化草甸開墾而來的耕地，為全年一熟一冬休閑型，屬于寒地稻田，土壤理化性質詳見表1。

稻田蒸散通過非稱重式蒸滲儀（有底測坑）進行測定，測坑規格為2 m×2 m，深30 cm，有底測坑底部及四壁用0.06 mm的塑料膜完全圍隔，由于避免了滲漏的影響，測坑內的耗水量均是通過蒸散損失掉。每天記錄并測定有底測坑內的灌水量，降水量及水層或土壤水分的變化，由水量平衡法即可得到相應時段相應灌水處理的稻田實測蒸散量。

常規氣象資料選取2011年整個生長期的農田綜合試驗場的上下換位式波文比小氣候觀測系統30 min觀測資料。

1.2 研究方法

單作物系數法是計算作物系數的標準方法之一，也稱為作物分段單值平均法[11]。該方法把作物系數的變化過程概括為幾個階段，根據各階段植物蒸騰和土壤蒸發的變化規律，將描述實際作物表層與參考表層蒸騰與蒸發作用不同的系數Kcb和Ke歸結為一個系數Kc，即用一個時段平均值表示為該階段的作物系數，Kc=Kcb+Ke，從而簡化作物系數的計算過程，使得在灌溉規劃、設計和管理中的應用更為簡便。

1.2.3 蒸散量估算評估指標 R2描述了觀測指標的整個變異能夠被模型所解釋的比例，R2對異常值是非常敏感的，因此又選擇了另外3個指標對估算方法進行評估：均方根誤差（RMSE），絕對平均誤差（MAE），一致性指數（d）。

2 結果與分析

2.1 基于波文比能量平衡的三江平原寒地稻田蒸散量的估算

波文比法估算的蒸散量與稻田實測蒸散量的季節變化對比見圖1。由圖1可知，估算值和實測值季節變化規律基本相同，蒸散量波動比較大，特別在水稻生長前期，這主要與前期氣象因子波動較大有關。但與實測值相比，波文比法估算的蒸散量普遍偏低，波文比法估算的蒸散量在整個生長季的日均值為5.0 mm/d，較實測日均值（5.5 mm/d）偏低9.1%。原因可能是，波文比在測量過程中主要考慮的是能量平衡，主要受太陽凈輻射、氣溫、相對濕度等氣象因子的影響，忽略了空氣動力學方面的知識，并且沒有考慮到植株的生長情況，并且由于人為觀測原因，實測值本身與實際蒸散量也存在一定的誤差。

圖2為波文比能量平衡法估算的蒸散量（ETbowen）與實測蒸散量（ETa）的相關關系情況。估算值與觀測值有較好的相關性（R2=0.87），RMSE值為0.75 mm/d，MAE值為0.64 mm/d，與觀測值的一致性d為0.79。但估算值總體比實測值低，因此在使用波文比能量平衡法估算三江平原井灌區蒸散量時應該對估算值進行一定的修訂。

2.2 基于PM-Kc模型的三江平原寒地稻田蒸散量的估算

FAO56推薦的水稻生長初期、生長中期和后期的作物系數典型值Kc（Tab）分別為Kcini=1.1，Kcmid=1.2，Kcend=0.9，并要求根據當地的濕潤頻率、作物高度和氣候條件對其進行修正，修正的公式為：

Kc=Kc（Tab）+[0.04（u2-2）-0.004（RHmin-45）]（h/3）0.3

式中，RHmin為計算時段內每日最小相對濕度的平均值，20%

利用校正的Kc值與PM-Kc模型得到2011年模擬ET與實測ET比較的結果見表3。除了2011年生長初期以外，PM-Kc模型低估了不同生育期觀測的累積的ET為5.87%～9.18%，其差值小于10%。而2011年的生長初期，PM-Kc方法得到的累積ET比實測累積的ETa高19.34%，這可能與插秧后一周稻田未灌水而使土壤含水量較低蒸發降低，以及6月溫度較低產生低溫冷害造成蒸騰降低而使實際ETa較小有關，說明在水分受到限制并且受低溫影響的水田Penman-Monteith模型模擬的結果誤差較大，需要進一步調整模型需要的參數，以便更準確預測水分脅迫和冷害環境脅迫下的實際的蒸散。

以2011年生長初期為例，進一步分析PM-Kc模型與實際蒸散ETa的關系，展現了隨觀測平均間隔天數的增加模型估算值與觀測值更加接近的趨勢（表4）。2011年均以PM-Kc模型模擬的七日蒸散平均為最好，PM-Kc模型七日平均值與觀測值相關性最好（R2=0.71），與觀測值一致性最好（d=0.75），并且有最低的誤差（RMSE及MAE值分別為1.01 mm/d，1.05 mm/d）。PM-Kc模型五日平均值排名第二，PM-Kc模型3日平均值排名第三，最后是PM-Kc模型一日平均值。因此，PM-Kc模型七日平均蒸散估算效果最好，其次是PM-Kc模型五日平均，PM-Kc模型三日和一日平均值與實測值誤差較大，最差。

3 小結與討論

農田蒸散量的實測難度很大，而在缺乏足夠氣象資料、土壤數據和作物生理學參數情況下，尋找簡便、可行的方法計算田間蒸散量，可以為制定合理灌溉制度提供可靠的定量化的理論支持。不同農田蒸散量估算方法在不同氣候區、不同農藝管理條件下的估算效果也是不盡相同的[10-15]。本研究探討了波文比能量平衡法和Penman-Monteith參考作物蒸散-作物系數法在寒地稻田蒸散量計算中的適用性。

波文比能量平衡法估算的蒸散量較實測蒸散量偏低9.1%，但總體有較好的相關性（R2=0.87），RMSE值為0.75 mm/d，MAE值為0.64 mm/d，與觀測值的一致性指數d為0.79。結合當地的土壤和大氣狀況對Penman-Monteith參考作物蒸散-作物系數法中的Kc值作了調整，利用校正的Kc值與Penman-Monteith模型得到的估算ET與實測ET比較的結果表明，PM-Kc模型低估了2011年水稻生長中期和后期的累積ET，但其相對誤差均小于10%，而2011年的生長初期，由于水分脅迫和冷害環境脅迫，PM-Kc模型得到的累積ET比實測的累積ETa高19.34%，這說明在某些天氣條件下，估算值與實測值有較大的偏差，這與陳四龍[16]的研究結果一致；以2011年生長初期為例進一步分析PM-Kc模型與實際蒸散ETa的關系，展現了隨著間隔日數的增加，RMSE、MAE有降低的趨勢，一致性系數有提高的趨勢，PM-Kc模型得到的七日平均模擬實際蒸散與實測蒸散數據最為接近，這表明可采用Penman-Monteith方程與校正的作物系數Kc得到的蒸散五日或七日平均值代表三江平原寒地稻田生育期內的實際蒸散量，這對于簡化寒地稻田水量平衡研究以及在水稻生產中評估和改善水分管理實踐具有重要意義。

參考文獻：

[1] 劉曉英，李玉中，郝衛平.華北主要作物需水量近50年變化趨勢及原因[J].農業工程學報，2005，21（10）：155-159.

[2] 尹海霞，張 勃，王亞敏，等.黑河流域中游地區近43年來農作物需水量的變化趨勢分析[J].資源科學，2012，34（3）：409-417.

[3] 王 瑤，趙傳燕，田風霞，等.黑河中游春小麥需水量空間分布[J].生態學報，2011，31（9）：2374-2382.

[4] GUNTER W，MARIJN VAN DER V，ALBERTO A，et al. Estimating irrigation water requirements in Europe[J]. Journal of Hydrology，2009，373（3）：527-544.

[5] 胡志橋，田霄鴻，張久東.石羊河流域主要作物的需水量及需水規律的研究[J].干旱地區農業研究，2011，29（3）：1-6.

[6] NIE X，WANG Y Y， LIU X T ，et al. Influence of intermittent irrigation on water consumption and yield of cold rice in Northeast China[J].Journal of Food，Agriculture Environment，2011，9（34）：315-320.

[7] 聶 曉，王毅勇，劉興土，等.節水灌溉對三江平原寒地水稻生理生態需水和產量的影響[J].華北農學報，2011，5（26）：163-168.

[8] 張世春，聶 曉，王毅勇.淺水灌溉對寒地井灌稻田溫度和產量構成因素的影響[J].東北林業大學學報，2011，39（12）：86-89.

[9] 陳剛起.三江平原沼澤研究[M].北京：中國科學出版社，1996.15-26.

[10] 王笑影，梁文舉，聞大中.北方稻田蒸散需水分析及其作物系數確定[J].應用生態學報，2005，16（1）：69-72.

[11] 劉 鈺，PEREIRA L S.對FAO推薦的作物系數計算方法的驗證[J].農業工程學報，2000，16（5）：26-30.

[12] KUSTAS W P，LI F，JACKSON T J，et al. Effects of remote sensing pixel resolution on modeled energy flux variability of croplands in Iowa[J]. Remote Sensing of Environment，2004， 92（4）：535-547.

[13] 賈 紅，胡繼超，張佳寶，等.應用Shuttleworth-Wallace模型對夏玉米農田蒸散的估計[J].灌溉排水學報，2008，27（4）：77-80.

[14] ALKHAIER F，SU Z B，FLERCHINGER G N. Reconnoitering the effect of shallow groundwater on land surface temperature and surface energy balance using MODIS and SEBS[J].Hydrology and Earth System Sciences Discussions，2012，16（7）：1833-1844.

[15] 趙春江，楊貴軍，薛緒掌，等.基于互補相關模型和IKONOS數據的農田蒸散時空特征分析[J].農業工程學報，2013，29（8）：115-124.

[16] 陳四龍.節水條件下冬小麥農田蒸散規律及其模擬研究[D].北京：中國科學院遺傳與發育生物學研究所，2006.