基于氣候適宜度的水稻發(fā)育期預(yù)報(bào)模型

孫貴拓　楊若翰　楊柯　王鵬云　王輝　曾艷　胡春偉

摘要?利用昆明1994—2012年水稻發(fā)育期觀測(cè)資料和地面氣象觀測(cè)資料，計(jì)算各發(fā)育階段綜合氣候適宜度，構(gòu)建基于氣候適宜度的水稻發(fā)育期預(yù)報(bào)模型。結(jié)果表明，各發(fā)育階段間隔日數(shù)預(yù)報(bào)模型預(yù)測(cè)值與實(shí)測(cè)值平均絕對(duì)誤差均在3?d以內(nèi)，其中拔節(jié)—孕穗期平均絕對(duì)誤差為0;其次為三葉—返青、孕穗—抽穗、乳熟—成熟期平均絕對(duì)誤差均為1?d;返青—分蘗、分蘗—拔節(jié)期平均絕對(duì)誤差為2～3?d。全生育期預(yù)測(cè)值與實(shí)測(cè)值平均絕對(duì)誤差為1?d。預(yù)報(bào)模型模擬效果較好，可在水稻發(fā)育期預(yù)報(bào)業(yè)務(wù)中推廣應(yīng)用。

關(guān)鍵詞?氣候適宜度;水稻;發(fā)育期;預(yù)報(bào)模型

中圖分類號(hào)?S162文獻(xiàn)標(biāo)識(shí)碼?A

文章編號(hào)?0517-6611（2019）16-0231-04

doi：10.3969/j.issn.0517-6611.2019.16.066

開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Prediction?Model?of?Rice?Developmental?Phase?Based?on?Climatic?Suitability

SUN?Gui?tuo，YANG?Ruo?han，YANG?Ke?et?al?（Kunming?Agrometeorological?Station?of?Yunnan?Province，Kunming，Yunnan?650228）

Abstract?Based?on?the?observed?data?of?rice?developmental?phase?and?weather?data?from?1994?to?2012?in?Kunming，the?comprehensive?climatic?suitability?at?each?developmental?phase?was?calculated?to?construct?prediction?model?of?rice?developmental?phase.The?average?absolute?error?between?the?predicted?and?measured?values?of?the?prediction?model?of?the?interval?days?of?each?development?phase?was?within?3?days?all.The?average?absolute?error?of?the?joint?booting?stage?was?0，followed?by?the?three?leaves?returning，booting?heading，milking?maturity?mature?stage，the?average?absolute?error?was?1?d，and?the?returning?tillering，tilleringjointing?stage?were?2-3?d?last.The?average?absolute?error?between?the?predicted?and?measured?values?of?the?whole?growth?period?was?1?d.The?prediction?model?had?a?good?predicted?effect?and?can?be?popularized?and?applied?in?the?prediction?business?of?rice?developmental?stages.

Key?words?Climatic?suitability;Rice;Developmental?phase;Prediction?model

水稻是我國(guó)播種面積最大、總產(chǎn)量最多、單產(chǎn)量最高的糧食品種，在糧食生產(chǎn)和消費(fèi)中處于主導(dǎo)地位[1]。云南水稻播種面積約占全省農(nóng)作物播種面積的1/3，水稻是云南主要的糧食作物[2]。

云南地處低緯高原，地形復(fù)雜，形成了獨(dú)特的氣候，其在全國(guó)雜交水稻種植區(qū)屬最佳種植區(qū)域之列[3]，但與同緯度相鄰地區(qū)比較，水稻生長(zhǎng)期氣候條件的特點(diǎn)是氣溫偏低、水稻生育期長(zhǎng)、播種育秧期常有“倒春寒”、中期高溫不足[4]。昆明位于中國(guó)氣候的脆弱帶?，同時(shí)受到東亞、西南兩支季風(fēng)進(jìn)退異常和明顯年際變化的影響，氣候變率極不穩(wěn)定，氣象災(zāi)害頻發(fā)[5]。尤其是在全球變暖的氣候背景下，不同地區(qū)水稻種植的氣候適宜性也發(fā)生了變化，研究氣候變化對(duì)水稻生產(chǎn)的影響[6-7]已成為全球氣候變化研究中的一項(xiàng)重要課題，并且水稻發(fā)育期模型研究更是其中的關(guān)鍵環(huán)節(jié)[8]。加強(qiáng)水稻農(nóng)業(yè)氣象服務(wù)，對(duì)防御農(nóng)業(yè)氣象災(zāi)害，采用合理的栽培管理措施有十分重要的意義。然而，這些服務(wù)的前提是需要預(yù)知水稻未來所處的發(fā)育狀態(tài)，才能因時(shí)、因地、有針對(duì)性地進(jìn)行農(nóng)業(yè)氣象災(zāi)害預(yù)警和管理決策服務(wù)。因此，研究水稻發(fā)育期預(yù)報(bào)方法、準(zhǔn)確提供發(fā)育期預(yù)報(bào)十分必要。

作物發(fā)育期的預(yù)報(bào)即是對(duì)作物從一個(gè)發(fā)育期到另一個(gè)相鄰發(fā)育期所經(jīng)歷時(shí)間間隔天數(shù)的預(yù)報(bào)，國(guó)內(nèi)外發(fā)展了很多作物發(fā)育期模擬模型，法國(guó)人?Reanmar早在1735年就創(chuàng)建基于積溫模擬的作物發(fā)育期模型（growing?degree?days，GDD），積溫模型隨后在國(guó)內(nèi)被廣泛應(yīng)用[9-10];高亮之等[11]提出適合于計(jì)算機(jī)應(yīng)用的動(dòng)態(tài)模擬模型——水稻鐘模型;戚昌瀚[12]、沈國(guó)權(quán)[13]建立了水稻發(fā)育溫度效應(yīng)的非線性模型。近年來，人們認(rèn)識(shí)到采用單一氣象因子的作物發(fā)育期模型已不能完全解釋作物生長(zhǎng)發(fā)育進(jìn)程，隨著作物發(fā)育期模型研究的深入，模型得到不斷改善，眾多學(xué)者提出了以生理發(fā)育時(shí)間（physiological?development?time，PDT）[14-15]為尺度，考慮了光、溫因子的作用，建立作物發(fā)育期預(yù)測(cè)模型，但模型中同樣沒有考慮水分條件對(duì)作物發(fā)育速率的影響。為綜合考慮光、溫、水等氣象條件對(duì)水稻生長(zhǎng)發(fā)育進(jìn)程的影響，任玉玉等[16]提出了作物氣候適宜度評(píng)價(jià)方法;孔令帥等[17]將氣候適宜度應(yīng)用于水稻產(chǎn)量及冬小麥發(fā)育期預(yù)報(bào)，取得了很好的效果。但采用氣候適宜度對(duì)水稻發(fā)育期模擬的研究尚不多見，筆者采用氣候適宜度對(duì)水稻發(fā)育期模型進(jìn)行模擬預(yù)報(bào)，為水稻生產(chǎn)氣象服務(wù)提供科學(xué)依據(jù)。

1?資料與方法

1.1?資料來源?選取昆明國(guó)家氣候基準(zhǔn)站1994—2014年逐日平均氣溫、日照時(shí)數(shù)、降水量資料，以及昆明國(guó)家農(nóng)業(yè)氣象試驗(yàn)站1994—2014年平行觀測(cè)的水稻發(fā)育期數(shù)據(jù)。?根據(jù)水稻生長(zhǎng)發(fā)育及昆明地區(qū)水稻種植特點(diǎn)，因播種—三葉期水稻為覆膜或搭棚育秧，其生長(zhǎng)發(fā)育受膜或棚內(nèi)小氣候條件影響，因此，該研究將水稻全生育期劃分為三葉—返青、返青—分蘗、分蘗—拔節(jié)、拔節(jié)—孕穗、孕穗—抽穗、抽穗—乳熟、乳熟—成熟7個(gè)生長(zhǎng)發(fā)育階段。

1.2?研究方法

1.2.1?溫度適宜度。

根據(jù)溫度三區(qū)間理論[18]，影響作物生長(zhǎng)發(fā)育的溫度指標(biāo)包括下限溫度、適宜溫度和上限溫度，作物各生長(zhǎng)發(fā)育階段三基點(diǎn)溫度各不相同。適宜溫度能夠促進(jìn)作物生長(zhǎng)，而極端溫度通過對(duì)作物組織的破壞作用影響干物質(zhì)累積，抑制作物生長(zhǎng)，溫度適宜度計(jì)算公式為：

S（T）=[（T-T?1）（T?h-T）B]/[（T?0-T?1）（T?h-T?0）B]（1）

B=（T?h-T?0）/（T?0-T?1）（2）

式中，S（T）為溫度適宜度，T為某日平均氣溫，T?0?、T?1和T?h分別為水稻各發(fā)育期所需的適宜溫度、下限溫度和上限溫度，水稻各生育期T?0、T?1、T?h不同（表1）。當(dāng)T≤T?1或T≥T?h時(shí)，S（T）=0;當(dāng)T=T?0?時(shí)，S（T）=1;當(dāng)T?1

1.2.2?日照適宜度。日照適宜度計(jì)算公式為：

S（S）=e??-[（S-S?0）]/b]2S<S?0

1S≥S?0（3）

式中，S（S）為水稻某發(fā)育階段某日日照適宜度，S為某日日照時(shí)數(shù)，S?0為臨界日照時(shí)數(shù)，即日照百分率達(dá)70%時(shí)的日照時(shí)數(shù)，b為經(jīng)驗(yàn)常數(shù)，水稻不同發(fā)育階段S?0和b的取值參照表2。當(dāng)實(shí)際日照時(shí)數(shù)大于臨界日照時(shí)數(shù)（最適宜日照時(shí)數(shù)）?時(shí)，S（S）?為?1，此時(shí)光照對(duì)水稻生長(zhǎng)發(fā)育最為適宜，S（S）也是在0～1取值的隸屬函數(shù)。

1.2.3?降水適宜度。降水適宜度計(jì)算公式如下：

S（R）=R/R?0R<R?0?R?0/RR≥R?0?（4）

式中，S（R）為降水適宜度，R為水稻某發(fā)育階段降水量，R?0為水稻某發(fā)育階段需水量。因水稻生長(zhǎng)期水分來源不僅是自然降水，即使降水不足，也可以通過人工灌溉補(bǔ)充水分，且昆明水稻種植區(qū)灌溉條件較好，在計(jì)算光溫水綜合適宜度時(shí)，默認(rèn)水分適宜度?S（R）為0.7[19]。

1.2.4?氣候適宜度。氣候適宜度能反映光、溫、水對(duì)水稻生育期的綜合影響，其計(jì)算式為：

S=3S（T）×S（R）×S（S）?（5）

式中，S為某日水稻氣候適宜度，S（T）、S（R）、S（S）為分別對(duì)應(yīng)的溫度、降水、日照適宜度。

1.2.5?發(fā)育日數(shù)與氣候適宜度。

當(dāng)氣候適宜度（S）計(jì)算值為1時(shí)，表明該日非常適宜水稻生長(zhǎng)，可將這一日確定為1個(gè)生長(zhǎng)發(fā)育日;當(dāng)氣候適宜度（S）計(jì)算值為0～1時(shí)，說明該日的氣象條件對(duì)水稻生長(zhǎng)發(fā)育有影響，可將這一日確定為相應(yīng)數(shù)值的生長(zhǎng)發(fā)育日數(shù);當(dāng)氣候適宜度（S）計(jì)算值為0時(shí)，說明該日不利于水稻生長(zhǎng)發(fā)育，水稻停止生長(zhǎng)。即某日的氣候適宜度值，就是水稻生長(zhǎng)發(fā)育日數(shù)，相鄰2個(gè)發(fā)育期間的發(fā)育日數(shù)由對(duì)應(yīng)時(shí)段內(nèi)的逐日氣候適宜度累加求得。

根據(jù)天氣預(yù)報(bào)產(chǎn)品，利用（1）～（5）式計(jì)算水稻某一發(fā)育階段逐日氣候適宜度，結(jié)合歷年水稻發(fā)育期觀測(cè)資料，統(tǒng)計(jì)水稻從一個(gè)發(fā)育期到進(jìn)入下一個(gè)相鄰發(fā)育期所需的累計(jì)氣候適宜度，計(jì)算該發(fā)育階段累積氣候適宜度的多年平均值，作為該發(fā)育階段的氣候適宜度生理發(fā)育指標(biāo)。各發(fā)育階段氣候適宜度生理發(fā)育指標(biāo)計(jì)算公式為：

Ijk=mi=1S?i（6）

I?k=1nnj=1Ijk（7）

式中，Ijk表示第k年第j個(gè)發(fā)育階段累計(jì)氣候適宜度，m為各發(fā)育階段持續(xù)日數(shù)，I?k為各發(fā)育階段生理發(fā)育指標(biāo)，n為發(fā)育期觀測(cè)資料年數(shù)。

2?結(jié)果與分析

2.1?發(fā)育期氣候適宜度預(yù)報(bào)模型

利用1994—2012年昆明水稻不同發(fā)育階段的發(fā)育日數(shù)（x）和相鄰2個(gè)發(fā)育期實(shí)際間隔日數(shù)（y）建立不同發(fā)育階段持續(xù)天數(shù)的預(yù)報(bào)模型y=ax+b，各發(fā)育階段模型參數(shù)a、b取值及相關(guān)系數(shù)見表3。結(jié)果顯示，除拔節(jié)—孕穗、乳熟—成熟2個(gè)發(fā)育階段模型相關(guān)系數(shù)稍差外，分別為0.753、0.768，其他生長(zhǎng)發(fā)育階段模擬模型相關(guān)性較好，均達(dá)到了0.800以上，三葉—返青、孕穗—抽穗的相關(guān)系數(shù)達(dá)到了0.900以上。

2.2?發(fā)育期動(dòng)態(tài)預(yù)報(bào)

采用1994—2014年水稻發(fā)育期觀測(cè)數(shù)據(jù)及氣象數(shù)據(jù)，計(jì)算并統(tǒng)計(jì)得到三葉—返青、返青—分蘗、分蘗—拔節(jié)、拔節(jié)—孕穗、孕穗—抽穗、抽穗—乳熟、乳熟—成熟各發(fā)育階段氣候適宜度生理發(fā)育指標(biāo)分別為16.2、9.8、15.1、6.9、4.8、13.8、12.0。當(dāng)某一發(fā)育期的具體日期確定時(shí)，根據(jù)最新天氣預(yù)報(bào)產(chǎn)品光、溫、水等氣象要素預(yù)測(cè)值，計(jì)算逐日氣候適宜度并進(jìn)行累加，當(dāng)累積氣候適宜度達(dá)到該發(fā)育階段的氣候適宜度生理發(fā)育指標(biāo)時(shí)，即說明水稻進(jìn)入下一發(fā)育期，該日即為下一相鄰發(fā)育期的具體預(yù)報(bào)日期。

2.3?檢驗(yàn)

2.3.1?模型檢驗(yàn)。

采用（1）～（5）式計(jì)算各發(fā)育階段逐日氣候適宜度，將各發(fā)育階段逐日氣候適宜度累加即為相應(yīng)發(fā)育階段間隔日數(shù)模擬預(yù)測(cè)值，比較各發(fā)育階段間隔日數(shù)實(shí)測(cè)值與模擬預(yù)測(cè)值。結(jié)果顯示（圖1），各發(fā)育階段間隔日數(shù)實(shí)測(cè)值與模擬預(yù)測(cè)值變化趨勢(shì)基本一致，實(shí)測(cè)值與預(yù)測(cè)值平均絕對(duì)誤差為1～4?d，其中孕穗—抽穗階段平均絕對(duì)誤差最小，為1?d;分蘗—拔節(jié)階段平均絕對(duì)誤差最大，為4?d;全生育期平均絕對(duì)誤差為3?d，模擬效果較好。

2.3.2?預(yù)報(bào)檢驗(yàn)。

利用2013和2014年昆明國(guó)家基準(zhǔn)氣候站地面觀測(cè)氣象數(shù)據(jù)及相應(yīng)時(shí)段昆明國(guó)家農(nóng)業(yè)氣象試驗(yàn)站水稻發(fā)育期大田觀測(cè)數(shù)據(jù)，對(duì)各發(fā)育階段發(fā)育期預(yù)報(bào)模型進(jìn)行

預(yù)報(bào)檢驗(yàn)（表4）。各生育期實(shí)測(cè)值與預(yù)測(cè)值平均絕對(duì)誤差均在3?d以內(nèi)，其中拔節(jié)—孕穗期平均絕對(duì)誤差為0，模擬預(yù)報(bào)效果最好;其次為三葉—返青、孕穗—抽穗、乳熟—成熟期，平均絕對(duì)誤差均為1?d;返青—分蘗、分蘗—拔節(jié)期預(yù)報(bào)效果稍差，平均絕對(duì)誤差為2～3?d。全生育期預(yù)測(cè)值與實(shí)測(cè)值平均絕對(duì)誤差為1?d，且各發(fā)育階段預(yù)測(cè)值與實(shí)測(cè)值散點(diǎn)圖比較集中（圖2），預(yù)報(bào)模型精確，能滿足業(yè)務(wù)應(yīng)用需求。

3?結(jié)論與討論

溫度、降水、日照等氣象條件對(duì)作物的生長(zhǎng)發(fā)育速度和產(chǎn)量形成有顯著影響[20]。該研究在考慮溫度、降水、日照因子綜合作用基礎(chǔ)上，基于氣候適宜度方法，建立了水稻三葉—成熟期各發(fā)育階段的發(fā)育期預(yù)報(bào)模型，并可鏈接天氣預(yù)報(bào)業(yè)務(wù)，實(shí)現(xiàn)水稻發(fā)育期動(dòng)態(tài)預(yù)報(bào)業(yè)務(wù)化應(yīng)用，相對(duì)于單一采用積溫或只考慮光、溫因子影響建立的預(yù)報(bào)模型有較大的改進(jìn)[21-22]。模型檢驗(yàn)結(jié)果表明，各發(fā)育階段間隔日數(shù)預(yù)測(cè)值與實(shí)測(cè)值一致性較好，平均絕對(duì)誤差均在3?d以內(nèi)，全生育期預(yù)報(bào)值與實(shí)測(cè)值平均絕對(duì)誤差為1?d，拔節(jié)—孕穗期模擬效果最好，平均絕對(duì)誤差為0;返青—分蘗、分蘗—拔節(jié)期預(yù)報(bào)效果稍差，平均絕對(duì)誤差為2～3?d。

由于水稻是一種灌溉作物，土壤水分的遷移和存儲(chǔ)過程十分復(fù)雜，而該研究的降水適宜度只考慮了大氣降水與作物需水的關(guān)系，未對(duì)土壤水分遷移和存儲(chǔ)做詳細(xì)探討。除此之外，作物的生長(zhǎng)發(fā)育進(jìn)程還與其生物學(xué)特性有關(guān)，吳玉潔等[23]利用不同積溫計(jì)算方法對(duì)作物發(fā)育期模擬效果對(duì)比時(shí)發(fā)現(xiàn)，分品種的模擬效果較不分品種的模擬效果更優(yōu)，該研究未對(duì)不同年份不同品種的水稻發(fā)育期模擬效果進(jìn)行分類比較，有關(guān)品種對(duì)水稻發(fā)育期預(yù)報(bào)模型模擬效果的影響機(jī)理有待進(jìn)一步完善。

參考文獻(xiàn)

[1] 段居琦.我國(guó)水稻種植分布及其對(duì)氣候變化的響應(yīng)[D].北京：中國(guó)氣象科學(xué)研究院，2012.

[2]?周以鴻.云南不同地區(qū)水稻產(chǎn)量與氣候因素關(guān)系的探討[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)，1990，5（1）：7-10.

[3]?朱勇.云南雜交水稻種植的氣候條件及區(qū)劃[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2000，15（1）：34-37.

[4]?鄒麗云，嚴(yán)華生.云南水稻氣候敏感性和栽培適應(yīng)性[J].南京氣象學(xué)院學(xué)報(bào)，1994，17（2）：249-254.

[5]?秦劍.氣候變化與昆明糧食生產(chǎn)的研究[J].南京氣象學(xué)院學(xué)報(bào)，1999，22（4）：671-679.

[6]?OSATO?T.The?assessment?of?the?rice?growing?stage?shift?under?climate?change?in?northern?part?of?Iwate，and?important?matters?to?obtain?for?stable?and?higher?rice?yields[J].Tohoku?journal?of?crop?science，2013，56：43-44.

[7]?MISHRA?A，SINGH?R，RAGHUWANSHI?N?S，et?al.Spatial?variability?of?climate?change?impacts?on?yield?of?rice?and?wheat?in?the?Indian?Ganga?Basin[J].The?science?of?the?total?environment，2013，468-469：132-138.

[8]?張帥，陶福祿.水稻發(fā)育期模型研究進(jìn)展[J].地理科學(xué)進(jìn)展，2012，31（11）：1485-1491.

[9]?錢拴，陳暉，王良宇.全國(guó)棉花發(fā)育期業(yè)務(wù)預(yù)報(bào)方法研究[J].應(yīng)用氣象學(xué)報(bào)，2007，18（4）：539-547.

[10]?張玉芳，王明田，游超.四川盆地小麥生育期預(yù)報(bào)方法研究與應(yīng)用[J].生物災(zāi)害科學(xué)，2013，36（4）：438-442.

[11]?高亮之，金之慶，黃耀，等.水稻計(jì)算機(jī)模擬模型及其應(yīng)用之一水稻鐘模型——水稻發(fā)育動(dòng)態(tài)的計(jì)算機(jī)模型[J].中國(guó)農(nóng)業(yè)氣象，1989，10（3）：3-10.

[12]?戚昌瀚.水稻生長(zhǎng)日歷模擬模型研究綜合報(bào)告[J].江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，1992，14（3）：218-223.

[13]?沈國(guó)權(quán).影響作物發(fā)育速度的非線性溫度模式[J].氣象，1980（6）：9-11.

[14]?張立禎，曹衛(wèi)星，張思平，等.基于生理發(fā)育時(shí)間的棉花生育期模擬模型[J].棉花學(xué)報(bào)，2003，15（2）：97-103.

[15]?王冀川，馬富裕，馮勝利，等.基于生理發(fā)育時(shí)間的加工番茄生育期模擬模型[J].應(yīng)用生態(tài)學(xué)報(bào)，2008，19（7）：1544-1550.

[16]?任玉玉，千懷遂.河南省棉花氣候適宜度變化趨勢(shì)分析[J].應(yīng)用氣象學(xué)報(bào)，2006，17（1）：87-92.

[17]?孔令帥，江勝國(guó)，張斗勝.基于氣候適宜度的一季稻產(chǎn)量預(yù)報(bào)模式研制：以安徽省桐城市為例[J].安徽農(nóng)業(yè)科學(xué)，2015，43（36）：286-288，291.

[18]?魏星，王品，張朝，等.溫度三區(qū)間理論評(píng)價(jià)氣候變化對(duì)作物產(chǎn)量影響[J].自然資源學(xué)報(bào)，2015，30（3）：470-479.

[19]?王輝，王鵬云，曾艷，等.昆明水稻生育期氣候因子定量評(píng)價(jià)[J].安徽農(nóng)業(yè)科學(xué)，2014，42（26）：9081-9082，9189.

[20]?王彥平，侯瓊，宋衛(wèi)士，等.氣候適宜度在內(nèi)蒙古東北部馬鈴薯發(fā)育期預(yù)報(bào)中的應(yīng)用[J].中國(guó)農(nóng)學(xué)通報(bào)，2015，31（3）：216-220.

[21]?趙峰，千懷遂，焦士興.農(nóng)作物氣候適宜度模型研究：以河南省冬小麥為例[J].資源科學(xué)，2003，25（6）：77-82.

[22]?李昊宇，王建林，鄭昌玲，等.氣候適宜度在華北冬小麥發(fā)育期預(yù)報(bào)中的應(yīng)用[J].氣象，2012，38（12）：1554-155.

[23]?吳玉潔，葉彩華，姜會(huì)飛，等.不同積溫計(jì)算方法作物發(fā)育期模擬效果比較[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2016，21（10）：117-126.