小時(shí)與日尺度PM公式的參照作物騰發(fā)量及其水稻單作物系數(shù)值差異

王可純 ， 蔡少杰， 衛(wèi) 琦,2， 呂玉平， 廖林仙， 徐俊增,2

(1.河海大學(xué) 南方地區(qū)高效灌排與農(nóng)業(yè)水土環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京 210098; 2.河海大學(xué) 水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京 210098)

1 研究背景

作物騰發(fā)量(crop evaportranspiration，ET)是制定作物灌溉制度與開展水資源配置的主要依據(jù)[1-4]。作物騰發(fā)量的準(zhǔn)確計(jì)算可在一定程度上節(jié)約灌溉用水量[5-6]，尤其是在農(nóng)業(yè)用水量較大、作物生長(zhǎng)對(duì)灌溉依賴程度較高的干旱和半干旱地區(qū)[7-10]。因此為了更好地管理作物灌溉用水量以及提高作物水分利用效率，必須準(zhǔn)確確定作物騰發(fā)量(ET)。

參考作物騰發(fā)量(reference crop evaportranspiration，ET0)是計(jì)算作物ET的關(guān)鍵參數(shù)，基于Penman-Monteith公式聯(lián)合國(guó)糧農(nóng)組織(FAO)和美國(guó)土木工程師協(xié)會(huì)(ASCE)給出了兩種計(jì)算ET0的標(biāo)準(zhǔn)：FAO定義的參考作物冠層阻力為70 s/m，反射率為0.23，并推薦改進(jìn)的PM公式計(jì)算ET0[11-12]；ASCE定義的矮型參考作物[13-15](0.12m高的修剪冷型草)冠層阻力rs按白晝與黑夜進(jìn)行了分別賦值，白天rs值為50 s/m，晚上rs值為200 s/m；而在以日為計(jì)算時(shí)段時(shí)，rs取70 s/m。此外，作物系數(shù)(crop coefficient，Kc)也是估算ET的關(guān)鍵因子[16-18]，其與作物的種類及生長(zhǎng)階段有關(guān)。不同作物上的研究表明，對(duì)Kc的區(qū)域性率定是準(zhǔn)確計(jì)算作物ET的關(guān)鍵[19]，隨著精準(zhǔn)灌溉的發(fā)展，小時(shí)尺度的作物騰發(fā)量估算越來(lái)越重要。以往對(duì)于日尺度作物騰發(fā)量的研究結(jié)果表明，兩種PM公式的計(jì)算結(jié)果關(guān)系密切，且ASCE PM公式的計(jì)算結(jié)果更為準(zhǔn)確[20-22]。而關(guān)于FAO56 PM和ASCE PM兩種公式在不同計(jì)算尺度下對(duì)日ET0計(jì)算結(jié)果的比較以及由于不同ET0計(jì)算對(duì)作物系數(shù)Kc取值的影響等方面的研究還較少。

因此，基于上述研究不足，本文擬以昆山2012-2013年常規(guī)小時(shí)氣象資料為基礎(chǔ)，以PM公式計(jì)算結(jié)果為對(duì)照，研究FAO56-PM和ASCE-PM兩種計(jì)算模型在不同時(shí)間尺度(小時(shí)、日)的ET0差異，分析不同日ET0計(jì)算策略對(duì)作物Kc值的影響，其研究結(jié)果對(duì)于華東地區(qū)節(jié)水農(nóng)業(yè)的灌溉制度制定具有重要的指導(dǎo)意義[23]。

2 材料與方法

2.1 試驗(yàn)區(qū)概況

試驗(yàn)區(qū)位于河海大學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室昆山試驗(yàn)研究基地(34°63′21″N，121°05′22″E)內(nèi)，該試驗(yàn)區(qū)屬于典型的亞熱帶季風(fēng)氣候區(qū)，多年平均氣溫15.5℃，多年平均蒸發(fā)量1 365.9 mm，多年平均降雨量1 097.1 mm，多年平均日照時(shí)數(shù)2 085.9 h，多年平均無(wú)霜期234 d。當(dāng)?shù)胤N植習(xí)慣為稻麥輪作，其中，水稻種植時(shí)間為6-10月份，小麥種植時(shí)間為11-次年5月，本試驗(yàn)主要在水稻生育期內(nèi)進(jìn)行。土壤屬潴育型黃泥土，土壤耕層質(zhì)地為重壤土，容重為1.30 g/cm3，pH值為7.4。

2.2 氣象因子觀測(cè)

2012-2013年常規(guī)氣象資料由自動(dòng)氣象站進(jìn)行觀測(cè)，主要包括空氣溫度(含最高溫度、最低溫度和日平均溫度)、氣壓、太陽(yáng)輻射、降雨量、相對(duì)濕度(含最高濕度、最低濕度和日平均濕度)、風(fēng)速和風(fēng)向等指標(biāo)，測(cè)定時(shí)間間隔為0.5 h。

2.3 實(shí)際騰發(fā)量ETa

實(shí)際騰發(fā)量ETa根據(jù)水量平衡原理計(jì)算而來(lái)，水量平衡公式如下：

h1+P+m-WC-d=h2

(1)

式中：h1為時(shí)段初田面水層深度；h2為時(shí)段末田面水層深度；P為時(shí)段內(nèi)降雨量；d為時(shí)段內(nèi)排水量；m為時(shí)段內(nèi)灌水量；WC為田間耗水量。以上單位均以mm計(jì)。實(shí)際騰發(fā)量ETa為田間耗水量WC與田間滲漏量之差。水層深度通過豎尺在固定觀測(cè)點(diǎn)觀測(cè)并記錄，灌溉水量直接由安裝在灌溉管道上的水表讀出，滲漏量則由翻斗式自動(dòng)測(cè)量裝置(精度為0.05 mm)測(cè)定。

2.4 ET0的標(biāo)準(zhǔn)化計(jì)算

ET0的標(biāo)準(zhǔn)化計(jì)算采用PM公式，PM公式在1998和2005年分別被FAO-56和ASCE-EWRI推薦使用，其具體計(jì)算公式如下：

(2)

式中：ET0為參考作物蒸騰蒸發(fā)量，mm/d或mm/h；G為土壤熱通量，MJ/(mm·d)或MJ/(mm·h)；Rn為輸入冠層的凈輻射，MJ/(m·d)或MJ/(m·h)；es、ea分別為飽和水氣壓與實(shí)際水汽壓，kPa；Δ為飽和水汽壓與氣溫關(guān)系曲線的斜率， kPa/℃；T為距地面2 m高處日平均氣溫，℃；u2為距地面2 m高處風(fēng)速，m/s；γ為濕度計(jì)常數(shù)，kPa/℃；Cn(K·mm·s3/(Mg·d))或(K·mm·s3/(Mg·h))和Cd(s/m)為由參考作物類型和計(jì)算時(shí)間步長(zhǎng)確定的常數(shù)(見表1)，0.408為計(jì)算常數(shù)，(m2·mm)/MJ。公式(2)中各參數(shù)的確定可參考文獻(xiàn)[2]和[7]。

小時(shí)尺度的參考作物騰發(fā)量記為ET0-h，由ASCE和FAO56 PM公式計(jì)算的小時(shí)ET0分別記為ET0-hA和ET0-hF，由ASCE和FAO56 PM公式計(jì)算的ET0-h累積得到的日參考作物騰發(fā)量分別記為ET0-dhA和ET0-dhF，由PM公式計(jì)算日尺度的參考作物騰發(fā)量記為ET0-d。

表1 Cn和Cd在FAO-56 PM 和ASCE PM中的取值

2.5 作物系數(shù)Kc的確定

節(jié)水灌溉條件下，一般需要考慮土壤水分脅迫系數(shù)Ks，可根據(jù)推薦的計(jì)算公式得到，公式如下[24]：

(3)

作物系數(shù)用下式求得：

(4)

式中:Kc為特定作物生長(zhǎng)階段和作物管理?xiàng)l件下的當(dāng)?shù)氐淖魑锵禂?shù)；ETa為實(shí)際作物騰發(fā)量，mm/d或mm/h；ET0為參照作物騰發(fā)量，mm/d或mm/h；Ks為水分脅迫系數(shù)，由ET0-d、ET0-dhA與ET0-dhF求得的水稻單作物系數(shù)分別記為Kc-d、Kc-dhA與Kc-dhF。

3 結(jié)果與討論

3.1 小時(shí)尺度參考作物騰發(fā)量ET0-h的結(jié)果分析

圖1為2012和2013年不同PM公式計(jì)算的小時(shí)尺度的參考作物騰發(fā)量(ET0-h)的結(jié)果。由圖1可以看出，兩種PM公式計(jì)算得到的ET0-h變化規(guī)律較為相似，但由于計(jì)算參數(shù)取值的不同，導(dǎo)致其在數(shù)值上有差異?？傮w上，由ASCE PM公式計(jì)算的ET0-h值范圍為-0.05～0.92 mm/h，由FAO56 PM公式計(jì)算的ET0-h值范圍為-0.05～0.89 mm/h。且當(dāng)ET0-h值越大時(shí)，二者之間的差異也越大，但總體差異維持在-0.02～0.5 mm/h范圍內(nèi)。其中在白天，由于Cd取值的不同，導(dǎo)致ASCE PMET0-h值大于FAO56 PMET0-h值，而在夜間，由于Cd取值以及能量項(xiàng)和空氣動(dòng)力學(xué)項(xiàng)的差異導(dǎo)致兩種計(jì)算公式ET0-h值大小關(guān)系并沒有一致的規(guī)律。此外，若2 m高處風(fēng)速u2接近或?yàn)?，不同的Cd的取值將不影響計(jì)算結(jié)果，兩種計(jì)算標(biāo)準(zhǔn)求得的ET0-h值基本相同。

進(jìn)一步分析ASCE PM公式和FAO56 PM公式在不同生育期的ET0-h日變化趨勢(shì)(圖2)，可以發(fā)現(xiàn)，在不同生育期內(nèi)，ET0-h變化均呈拋物狀分布，且ET0-h值在中午時(shí)刻(12：00左右)達(dá)到最大，而夜間則維持在較低水平，這一特征與太陽(yáng)凈輻射的分布較為相似，說(shuō)明太陽(yáng)凈輻射對(duì)ET0-h影響較大。兩種計(jì)算公式下ET0-h值在同一年內(nèi)差異較小，但在不同年際間有所差異。

3.2 日尺度參考作物騰發(fā)量ET0-d的結(jié)果分析

日參考作物騰發(fā)量由兩種PM公式計(jì)算得到，計(jì)算值與小時(shí)累積值對(duì)比見圖3。根據(jù)圖3可以看出，整體上，ET0-dhA、ET0-dhF與ET0-d之間具有良好的線性關(guān)系，其相關(guān)系數(shù)R2達(dá)0.97以上，且其大小關(guān)系表現(xiàn)為ET0-dhA>ET0-d>ET0-dhF，其中2012年的ET0-dhA、ET0-dhF與ET0-d之間的誤差分別為1.69%、1.81%，而2013年的誤差則略有所減小(1.19%和1.63%)。進(jìn)一步比較不同計(jì)算策略下日參考作物騰發(fā)量在不同生長(zhǎng)時(shí)期的日ET0均值(表2)，并結(jié)合作物不同生長(zhǎng)時(shí)期內(nèi)的氣象因子變化，可以發(fā)現(xiàn)，2013年全生育期內(nèi)的平均溫度(27.24℃)和平均凈輻射(126.67 W/m2)分別較2012年的值提高了7.6%和5.5%，對(duì)比2012和2013年不同生長(zhǎng)時(shí)期的日ET0值，可發(fā)現(xiàn)除生長(zhǎng)后期外，2013年日ET0均大于2012年，由此推測(cè)出造成年際間日ET0值發(fā)生變化的原因可能是由太陽(yáng)凈輻射Rn和溫度T在年際間的變化引起的。

3.3 不同時(shí)間尺度計(jì)算ET0對(duì)確定水稻單作物系數(shù)的差異

通過分析不同計(jì)算策略下日ET0對(duì)作物系數(shù)的影響(圖4)，可以看出，作物系數(shù)Kc值與日ET0之間具有良好的線性關(guān)系，其相關(guān)系數(shù)R2達(dá)0.90以上，且大小關(guān)系表現(xiàn)為Kc-dhF>Kc-d>Kc-dhA，其中2012年Kc-dhA、Kc-dhF與Kc-d之間的誤差為1.13%、4.61%，2013年誤差為2.37%、0.52%，整體誤差較小(0～5%)。進(jìn)一步分析作物不同生長(zhǎng)時(shí)期Kc的均值(表3)可以發(fā)現(xiàn)，總體上，Kc在生長(zhǎng)初期較小，而隨著作物的生育期的進(jìn)行，其Kc值逐漸增大，并在生長(zhǎng)中期達(dá)到最大，此后又隨作物的枯萎而逐漸減小。此外，通過進(jìn)一步分析由實(shí)測(cè)日騰發(fā)量ETa求得的Kc值和FAO推薦值之間的關(guān)系，可以發(fā)現(xiàn)，由水稻實(shí)測(cè)日騰發(fā)量ETa求得的Kcini、Kcmid和Kcend值2012年分別為1.20、1.60和0.76，2013年分別為1.14、1.63和0.96，均略大于由FAO推薦的單作物系數(shù)修正法計(jì)算得到的Kcini、Kcmid和Kcend的值(1.05、1.12、0.85)。由此可以看出，單作物系數(shù)Kc值的區(qū)域性率定對(duì)于作物騰發(fā)量ETa的準(zhǔn)確計(jì)算至關(guān)重要。

圖1 ASCE PM公式與FAO56 PM公式計(jì)算的小時(shí)作物騰發(fā)量ET0之間的對(duì)比

圖2 小時(shí)ET0日間變化趨勢(shì)

圖3 兩種PM公式計(jì)算的ET0-d值與小時(shí)累積值對(duì)比

表2 不同生長(zhǎng)期及全生育期日騰發(fā)量ET0均值 mm/d

表3 水稻單作物系數(shù)不同生長(zhǎng)期的均值及差異分析 mm/d

圖4 兩種PM公式計(jì)算的Kc-d值與小時(shí)累積值對(duì)比

4 結(jié) 論

以昆山試驗(yàn)基地2012-2013年自動(dòng)氣象站的小時(shí)氣象資料為依據(jù)，研究了FAO56 PM和ASCE PM公式在小時(shí)尺度和日尺度的ET0變化特征，揭示了由小時(shí)尺度ET0值累積得到的日ET0值和由日尺度計(jì)算得到的日ET0值之間的差異，評(píng)價(jià)了不同計(jì)算策略對(duì)作物系數(shù)(Kc)的影響。得出的主要結(jié)論如下：

(1)在小時(shí)尺度上，由ASCE PM公式和FAO56 PM公式計(jì)算得到的ET0-h值變化規(guī)律較為一致，但由于Cd取值的不同，導(dǎo)致其在數(shù)值上略有差異，且整體上白天的ASCE PMET0-h值大于FAO56 PMET0-h值，而夜間兩種公式計(jì)算的ET0-h大小之間沒有一致性規(guī)律。

(2)在日尺度上，不同計(jì)算策略下的日ET0計(jì)算結(jié)果之間表現(xiàn)出較好的線形關(guān)系，其大小關(guān)系總體表現(xiàn)為ET0-dhA>ET0-d>ET0-dhF。

(3)用不同日ET0計(jì)算策略下得到的節(jié)水灌溉稻田的Kc值，其大小關(guān)系表現(xiàn)為：Kc-dhA