王 婷 ，王 鵬 ，李雨鴻 ，周 斌 ，李琳琳 ，胡春麗 ，劉東明

（1.遼寧省氣象科學(xué)研究所，遼寧沈陽(yáng)110166；2.遼寧省氣象服務(wù)中心，遼寧沈陽(yáng)110166）

水稻是我國(guó)重要的糧食作物之一，種植面積占糧食作物總播種面積的30%，產(chǎn)量占糧食總產(chǎn)量的40%以上[1]。遼寧省作為我國(guó)重要的水稻產(chǎn)區(qū)和商品水稻生產(chǎn)基地，水稻產(chǎn)量與該省乃至全國(guó)糧食安全密切相關(guān)[2]。移栽是水稻種植的關(guān)鍵環(huán)節(jié)，要根據(jù)適宜移栽秧齡大小，科學(xué)安排移栽期，水稻移栽期推遲，容易出現(xiàn)超秧齡的秧苗，使水稻有效分蘗期縮短，不利于前期發(fā)苗和保足穗、攻大穗，對(duì)水稻生長(zhǎng)發(fā)育和產(chǎn)量形成有較大影響[3]。因此，做好水稻移栽期預(yù)報(bào)，合理安排水稻移栽日期，及早做好茬口銜接，努力做到適期移栽是保證和提高水稻產(chǎn)量的基礎(chǔ)。做好水稻移栽期預(yù)報(bào)是現(xiàn)今農(nóng)業(yè)氣象預(yù)報(bào)的重要任務(wù)之一，對(duì)農(nóng)業(yè)生產(chǎn)管理有著重要的指導(dǎo)作用。影響水稻移栽期的氣象因子主要有氣溫、降水和日照，其中氣溫對(duì)水稻的影響最大。而有效積溫更能反映水稻生育期對(duì)熱量的要求[4-5]，因此使用有效積溫可以對(duì)水稻生育期進(jìn)行預(yù)測(cè)[6-9]。

研究表明，大氣環(huán)流對(duì)區(qū)域氣候要素的變化有著直接和重要的影響[10]，本文旨在篩選出對(duì)有效積溫影響顯著的大氣環(huán)流因子，建立有效積溫預(yù)報(bào)模型，并通過(guò)有效積溫和水稻移栽期建立相關(guān)性，建立水稻適宜移栽期的長(zhǎng)期預(yù)報(bào)模型，最終確定水稻適宜移栽的日期。

1 材料與方法

1.1 農(nóng)業(yè)資料

農(nóng)業(yè)資料來(lái)源于遼寧省5個(gè)農(nóng)業(yè)氣象觀測(cè)站（大洼、東港、新賓、開(kāi)原、遼中）的1981—2016年水稻生育期資料。氣象資料來(lái)源于遼寧省氣象局以上5個(gè)站1981—2016年日平均氣溫。大氣環(huán)流因子采用國(guó)家氣候中心提供的74項(xiàng)環(huán)流指數(shù)。

1.2 建模方法

1.2.1 有效積溫的計(jì)算。每一種植物都有其生長(zhǎng)的下限溫度。當(dāng)溫度高于下限溫度時(shí)，它才能夠發(fā)育和生長(zhǎng)。但僅僅溫度達(dá)到所需還不足以完成發(fā)育和生長(zhǎng)，因?yàn)檫€需要一定的時(shí)間，即需要一定的總熱量，稱為總積溫或者有效積溫。水稻的有效積溫是指對(duì)其生長(zhǎng)發(fā)育起有效作用的溫度累加值，即日平均氣溫（ti）與對(duì)水稻生長(zhǎng)無(wú)效的下限溫度（b0）之差的累加值，計(jì)算公式如下：

注：文中統(tǒng)計(jì)的是水稻播種至移栽期≥10℃的有效積溫，所以b0＝10℃，n為播種至移栽期天數(shù)。

1.2.2 水稻適宜移栽期模型建立。通過(guò)水稻適宜移栽期日序與水稻播種到移栽期有效積溫的高相關(guān)性建立適宜移栽期的回歸方程，然后利用國(guó)家氣候中心提供的74項(xiàng)環(huán)流指數(shù)，通過(guò)相關(guān)分析，得出有效積溫與前一年和當(dāng)年1—4月共74項(xiàng)環(huán)流指數(shù)的相關(guān)系數(shù)，并篩選出通過(guò)0.05顯著性檢驗(yàn)且絕對(duì)值大于0.4的指數(shù)，并且遵循因子獨(dú)立性、連續(xù)性及物理意義明確等原則，選出影響有效積溫的主要環(huán)流因子，建立預(yù)報(bào)有效積溫的方程，最終建立水稻適宜移栽期的預(yù)報(bào)模型。為了檢驗(yàn)預(yù)報(bào)模型的效果，對(duì)模型進(jìn)行歷史擬合和預(yù)報(bào)檢驗(yàn)。

2 結(jié)果和分析

2.1 水稻適宜移栽期與播種到移栽期有效積溫的相關(guān)分析

從1981—2015年遼寧省水稻適宜移栽期的時(shí)間序列和有效積溫對(duì)比來(lái)看，兩者的年際變化基本同步，呈現(xiàn)極顯著的正相關(guān)（P＜0.01）。由此可見(jiàn)，利用有效積溫對(duì)適宜移栽期進(jìn)行預(yù)報(bào)是可行的，但由于遼寧省跨度較大，區(qū)域間存在較大差異，因此需要分區(qū)建立水稻移栽期的預(yù)報(bào)模型。

圖1 遼寧省水稻移栽期和播種到移栽期有效積溫的時(shí)間變化曲線

2.2 基于環(huán)流指數(shù)的水稻適宜移栽期預(yù)報(bào)模型建立

通過(guò)水稻適宜移栽期與播種到移栽期的有效積溫的高相關(guān)性建立一元回歸方程（表1），然后通過(guò)播種到移栽期有效積溫與大氣環(huán)流指數(shù)74項(xiàng)特征值的相關(guān)分析，篩選出與遼寧水稻有效積溫發(fā)生有顯著關(guān)系的大氣環(huán)流指數(shù)因子，采用逐步回歸方法建立有效積溫預(yù)報(bào)模型（表2）。影響大洼有效積溫的環(huán)流因子：前一年6月北半球極渦中心位置（H11），前一年 8 月亞洲經(jīng)向環(huán)流指數(shù)（H12），前一年10月東亞槽強(qiáng)度（H13），前一年3月南方濤動(dòng)指數(shù)（H14）；影響東港有效積溫的環(huán)流因子：前一年4月西太平洋副高脊線（H21），前一年10月西太平洋副高北界（H22），當(dāng)年 3 月北半球極渦中心位置（H23），當(dāng)年4月東亞槽位置（H24）；影響新賓有效積溫的環(huán)流因子：前一年7月西太平洋副高面積指數(shù)（H31），前一年10月東亞槽位置（H32），前一年10月東亞槽強(qiáng)度（H33），當(dāng)年 2 月西藏高原（H34）；影響開(kāi)原有效積溫的環(huán)流因子：前一年3月北半球極渦中心位置（H41），前一年9月大西洋歐洲環(huán)流型W（H42），前一年3月大西洋歐洲環(huán)流型C（H43），前一年9月東亞槽位置（H44），前一年 2 月東亞槽強(qiáng)度（H45），當(dāng)年 3月西太平洋副高面積指數(shù)（H46），當(dāng)年1月西太平洋副高北界（H47），當(dāng)年 2月亞洲緯向環(huán)流指數(shù)（H48）；影響遼中有效積溫的環(huán)流因子：前一年4月亞洲區(qū)極渦強(qiáng)度指數(shù)（H51），前一年 7 月西藏高原（H52），當(dāng)年4月東亞槽位置（H53）。

表1 遼寧水稻適宜移栽期與播種到移栽有效積溫的回歸方程

表2 利用大氣環(huán)流指數(shù)因子建立遼寧省水稻播種到移栽期有效積溫的預(yù)報(bào)模型

2.3 水稻適宜移栽期和有效積溫預(yù)報(bào)模型的效果檢驗(yàn)

為了檢驗(yàn)表1和表2中預(yù)報(bào)模型的擬合效果，利用表2中的預(yù)報(bào)模型可由前期大氣環(huán)流指數(shù)因子預(yù)報(bào)出當(dāng)年的播種到移栽期有效積溫，然后利用表1的回歸方程預(yù)報(bào)當(dāng)年的適宜移栽期日序，將實(shí)際值和模擬值進(jìn)行對(duì)比（圖2），可以看到模擬出的有效積溫和適宜移栽期都非常接近實(shí)際值。大洼、東港、新賓、開(kāi)原、遼中5個(gè)站點(diǎn)的適宜移栽期預(yù)報(bào)準(zhǔn)確率分別為98.0%、98.0%、98.4%、98.7%、97.1%。移栽期的模擬偏差主要在4 d以內(nèi)，模型通過(guò)0.01的顯著性水平檢驗(yàn)，擬合效果較好。利用2016年資料進(jìn)行試報(bào)檢驗(yàn)，表3為移栽期預(yù)報(bào)結(jié)果和實(shí)際移栽期對(duì)比，大洼、東港、新賓、遼中的誤差是在1～2 d，開(kāi)原移栽期預(yù)報(bào)比實(shí)際值偏早6 d。

圖2 遼寧省水稻適宜移栽期和播種到移栽期有效積溫預(yù)報(bào)模型擬合檢驗(yàn)

表3 2016年遼寧水稻適宜移栽期預(yù)報(bào)模型的試報(bào)檢驗(yàn)

3 結(jié)論與討論

本文根據(jù)有效積溫和環(huán)流因子的相關(guān)關(guān)系，篩選出對(duì)有效積溫影響顯著的大氣環(huán)流因子，通過(guò)0.05的顯著性檢驗(yàn)，建立有效積溫預(yù)報(bào)模型，并通過(guò)有效積溫和水稻移栽期之間的相關(guān)性，建立水稻適宜移栽期的長(zhǎng)期預(yù)報(bào)模型。結(jié)果表明：

1)遵循因子連續(xù)性和獨(dú)立性原則，篩選出對(duì)有效積溫影響較大的環(huán)流因子作為預(yù)報(bào)因子，基于多元線性回歸建立有效積溫預(yù)報(bào)模型，考慮到遼寧跨度較大，農(nóng)業(yè)氣候存在差異，因此分區(qū)選取代表站建立了模型。

2)遼寧省水稻適宜移栽期與播種到移栽≥10℃有效積溫存在極為顯著的線性正相關(guān)，可以利用有效積溫建立適宜移栽期的預(yù)報(bào)模型。預(yù)報(bào)模型歷史擬合效果較好，準(zhǔn)確率基本在95%以上，預(yù)報(bào)檢驗(yàn)誤差基本在1～2 d，預(yù)報(bào)準(zhǔn)確率較高，可以在業(yè)務(wù)中使用，預(yù)報(bào)結(jié)果對(duì)水稻種植具有重要的實(shí)用價(jià)值，可以提升農(nóng)業(yè)氣象的專業(yè)服務(wù)能力。

水稻的適宜移栽期除了受氣溫直接影響外，還應(yīng)考慮品種、降水、日照、土壤、社會(huì)因素等影響，本文僅通過(guò)有效積溫確定移栽期，還需兼顧其他因素繼續(xù)深入研究，以進(jìn)一步提高預(yù)報(bào)模型的準(zhǔn)確性。

參考文獻(xiàn):

[1]王 健.東北地區(qū)粳稻低溫冷害風(fēng)險(xiǎn)評(píng)估[D].北京:中國(guó)農(nóng)業(yè)科學(xué)院,2013.

[2]王晾晾,楊曉強(qiáng),李 帥,等.東北地區(qū)水稻霜凍災(zāi)害風(fēng)險(xiǎn)評(píng)估與區(qū)劃[J].氣象與環(huán)境學(xué)報(bào),2012,28(5):40-45.

[3]周 靜,孟桂元,馬國(guó)輝,等.超級(jí)雜交晚稻豐源優(yōu)299合理施氮量研究[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)),2013,28(4):536-539.

[4]薛大偉,方茂庭,錢 前.有效積溫在水稻生產(chǎn)中的應(yīng)用[J].中國(guó)稻米,2004,10(4):47-48.

[5]朱海霞,李秀芬,王 萍,等.黑龍江省水稻生長(zhǎng)季積溫計(jì)算方法[J].應(yīng)用氣象學(xué)報(bào),2017,28(2):247-256.

[6]沈國(guó)權(quán).當(dāng)量積溫及其應(yīng)用[J].氣象,1981(7):23-25.

[7]黃中雄,歐釗榮.改進(jìn)的當(dāng)量積溫法做水稻生育期預(yù)測(cè)[J].廣西農(nóng)學(xué)報(bào),1999(2):44-47.

[8]高 蘋,徐 敏,徐經(jīng)緯,等.基于大尺度海溫因子的江蘇省水稻適宜移栽期預(yù)報(bào)模型研究 [J].氣象,2015,41(12):1538-1544.

[9]高 蘋,張 佩,謝小萍,等.基于海溫和環(huán)流特征量的江蘇省小麥適播期預(yù)測(cè)[J].氣象,2012,38(12):1572-1578.

[10]龔道溢,王紹武.大氣環(huán)流因子對(duì)北半球氣溫變化影響的研究[J].地理研究,1999,18(1):31-38.