東北地區不同熟性水稻適宜種植區障礙型冷害時空變化*

梁立江，武永峰**，劉 聰,2，邵 帥

東北地區不同熟性水稻適宜種植區障礙型冷害時空變化*

梁立江1，武永峰1**，劉 聰1,2，邵 帥3

（1．中國農業科學院農業環境與可持續發展研究所/農業部農業環境重點實驗室，北京 100081；2．東華理工大學地球科學學院，南昌 330000；3.杭州電子科技大學計算機學院，杭州 310000）

利用東北地區167個氣象站1986-2015年的逐日氣象資料和28個物候觀測站1991-2008年的逐年水稻生育期數據，參照水稻障礙型冷害國家標準，分三個時間段（1986-1995年、1996-2005年和2006-2015年）劃分東北地區不同熟性水稻適宜種植區，利用冷害發生頻率和站次比分析東北地區近30a（1986-2015年）不同熟性水稻適宜種植區障礙型冷害的時空規律。結果表明：研究期內，東北地區不適宜種植區域面積逐漸減小，不同熟性適宜種植區呈現北移東擴的趨勢，其中以晚熟區擴大范圍最為明顯。不同熟性水稻適宜種植區障礙型冷害具有波動發生的特點，不同熟性區年平均冷害站次比大小表現為早熟區>晚熟區>中熟區；2006-2015年為早熟區冷害多發時段，1986-1995年為中熟區冷害多發時段，1996-2005年為晚熟區冷害多發時段；在早熟區冷害多發時段內，重度冷害的發生頻率減少；2006年以后中熟區和晚熟區冷害總發生頻率減小，但部分區域的不同等級冷害發生頻率提高。可見，東北地區障礙型冷害發生仍具有不確定性，即使在整體冷害頻率下降的情況下，個別區域在某些年份仍可能發生不同程度的冷害。

水稻；熟性；活動積溫；障礙型冷害；東北地區

東北地區地處北溫帶季風氣候區，水熱條件適宜，糧食增產潛力較大，為保障中國糧食安全發揮了重要作用[1]。同時，東北地區所處緯度較高，階段性低溫現象時有發生，作物生長季內因低溫產生的冷害是對該地區水稻生產影響最大的氣象災害[2]。7、8月是水稻由營養生長轉為生殖生長的關鍵時期，若水稻在此階段遭受短期異常低溫天氣過程，其生長活動將受抑制，導致生殖器官生理機能遭到破壞，引起障礙型冷害[3-4]，造成稻米品質下降，水稻產量降低。因此，充分認識東北地區水稻障礙型冷害變化規律，對保障糧食生產安全具有重要的現實意義。

國內學者針對東北地區水稻障礙型冷害時空變化規律已開展相關研究。胡春麗等[5]以氣象站逐日氣溫數據為基礎，分析了1961-2010年東北地區水稻障礙型冷害變化與區域氣候增暖之間的關系，發現水稻障礙型冷害對關鍵發育期氣溫變化響應較為敏感，東北地區近年來關鍵發育期氣溫呈上升趨勢，水稻障礙型冷害事件總體呈減少趨勢，但局部地區受氣候變率增大的影響，水稻障礙型冷害事件增加。張軼等[6]利用格點氣象數據和農業氣象站點生育期數據分析了1981-2010年東北地區水稻障礙型冷害的多年風險，認為東北地區障礙型冷害風險雖然整體呈下降趨勢，但障礙型冷害風險范圍具有波動變化，進入20世紀初，雖然氣溫明顯增加，但障礙型冷害風險敏感區域較大，且從時間序列上看黑龍江省和吉林省在2000年以后都存在高風險冷害年份。馮喜媛等[7]以氣象站逐日氣溫數據為基礎資料，采用氣象影響模型分析1961-2010年東北地區水稻障礙型冷害時空變化規律，表明東北地區水稻障礙型冷害主要發生在吉林省東北部和黑龍江省東南部，進入21世紀后，東北地區水稻障礙型冷害事件明顯增多，黑龍江省和吉林省可能進入新的冷害多發階段。在全球氣候變暖背景下，東北地區仍存在水稻冷害風險，水稻障礙型冷害是否進入新的多發階段值得研究[8]。目前，國內學者均將東北地區水稻視為單一的熟性品種對其障礙型冷害時空規律進行研究，所參照農業氣象災害標準未考慮不同熟性水稻的耐冷性差異。為此，本研究在劃分不同熟性水稻適宜種植區的基礎上，利用中國國家標準化管理委員會2017年發布的《北方水稻低溫冷害等級》（GB/T34967-2017）國家標準（下文簡稱國家標準）[9]，對1986-2015年東北地區水稻障礙型冷害發生特點以及隨年代發展的變化趨勢進行探究，以期揭示東北地區水稻障礙型冷害隨熟性變化的時空規律，從而為水稻減災保產、優化種植布局提供參考依據。

1 資料與方法

1.1 數據及其來源

研究區包括黑龍江、吉林和遼寧三省，位于121°10'-135°05'E，53°33'-38°43'N，面積78.73萬km2。該地區地形以平原、山地、河流為主，平原區土壤肥沃，擁有良好的水稻種植條件，是中國重要的水稻主產基地[10-11]。同時，東北地區也是水稻障礙型冷害的高發區域，每隔2～3a就會發生一次規模較大的冷害，冷害年水稻減產20%以上[12]。

氣象觀測數據和水稻生育期數據來自國家氣象信息中心，包括東北地區167個氣象站（其中，黑龍江省72站，吉林省49站，遼寧省46站）1986-2015年的逐日氣溫觀測數據和28個物候觀測臺站（其中，黑龍江省10站，遼寧省8站，吉林省10站）1991-2008年的逐年觀測數據，氣象站點和物候觀測站點分布見圖1。

原始氣象觀測資料中存在部分異常值和缺失值，故對數據集進行預處理，包括使用EXCEL軟件對數據異常值進行剔除以及使用SPSS軟件中線性插補方法對缺失值進行替代[13-14]。東北地區水稻障礙型冷害災情資料源自《中國氣象災害大典》（遼寧卷、吉林卷和黑龍江卷）、《中國氣象災害年鑒》（2002-2016年）。

圖1 東北地區氣象站、物候觀測站分布

1.2 不同熟性水稻適宜種植區劃分標準

根據宛敏渭等[15]的研究，當其它生態條件滿足作物基本的生長發育要求時，溫度對作物的生長發育起主導性作用，此時作物生長發育所需的下限溫度及有效積溫將成為制約因素。有效積溫在相鄰年份中變化較少，每個品種在其生育期內所要求的有效積溫基本固定，所以用有效積溫作為水稻熟性劃分的熱量指標較為合理。因此，首先參照國家標準所確定的≥10℃活動積溫指標（表1），對東北地區的不同熟性水稻適宜種植區進行劃分。

表1 不同熟性水稻適宜種植區積溫劃分標準

1.3 不同熟性水稻障礙型冷害判斷標準

水稻發生障礙型冷害的關鍵時期位于孕穗期與開花期之間[3-4]，參照國家標準以水稻抽穗前20d為孕穗期，以水稻抽穗后10d為開花期[9]。在獲取的167個氣象站點觀測數據中，僅28個臺站點擁有完整的水稻關鍵發育期資料，曲輝輝等[16]在水稻障礙型冷害時空變化研究中，利用霍普金斯物候定律[17]推算東北地區無觀測數據站點的水稻多年平均生育期，其結果表明該定律在東北地區適用性較好。因此，為獲取所有站點的水稻發育期數據，采用霍普金斯物候定律推算無觀測記錄站點的水稻發育期。該物候定律主要內容為，在北美洲溫帶內，在其它條件相同的情況下，每向北移動1°，向東移動5°或海拔上升121.92m，植物的發育期在春季和初夏將延遲4d，在秋季則相反，提早4d。

胡春麗等[18]研究表明，水稻受災強度受日均氣溫與持續日數的雙重影響，故參照國家標準，以水稻關鍵生育期（孕穗期和開花期）日均氣溫和持續日數作為致災判斷因子（表2），對東北地區不同熟性水稻障礙型冷害等級進行評估。若某站點在任一水稻關鍵生育期內，有一日或連續幾日達到障礙型冷害的氣象評價指標，則認為該站點當年發生了一次障礙型冷害；若某站點發生冷害的時間間隔在一日以上，則認為該站點當年發生了多次障礙型冷害，以該站點在關鍵生育期內達到的最嚴重冷害等級為該站點當年的冷害等級。

表2 不同熟性種植區水稻障礙型冷害判斷標準

注：抽穗前20d為孕穗期，抽穗后10d為開花期，Ta為抽穗前20d內（孕穗期）逐日觀測平均氣溫，Tb為抽穗后10d內（抽穗期）逐日觀測平均氣溫；符號“-”代表持續保持，符號“+”代表超過某一持續天數。

Note: The booting stage is the 20 days before heading, the flowering stage is the 10 days after heading. Ta is the daily mean temperature in 20 days before heading; Tb is the daily mean temperature in 10 days after heading. The symbol "-" indicates continuous, the sign "+" indicates more than a certain duration.

1.4 水稻障礙型冷害時空格局分析

采用某一年份冷害發生站數與總站數的比值（站次比F）表示該年冷害出現范圍的大小[19]，即

式中，m為各熟性區內發生冷害的站點數，M為各熟性區內站點總數，各熟性區不包括不適宜種植區。同時，利用多元線性回歸方法分析冷害站次比隨研究年份的線性變化趨勢。

以10a為一個時段（1986-1995年、1996-2005年和2006-2015年），采用某一時段出現冷害的年數與該時段總年數的比值（Pi）表示某一站該時段冷害的發生頻率（%），即

式中，n為該站點指定時段內發生某種等級冷害的年數，N為指定時段內總年數，i為用于區別不同站點的代號。

2 結果與分析

2.1 不同熟性水稻適宜種植區劃分及其變化

基于東北地區各氣象站點1986-2015年的逐日氣溫觀測數據，參照國家標準制定的80%保證率下≥10℃活動積溫指標，采用5日滑動平均法[16]判定≥10℃溫度的起止日期，即從一年中最長的一段≥10℃的5日滑動平均溫度序列中，選取第1個≥10℃的5日滑動平均組中的第1個≥10℃的日期作為起始日期，選取最后一個≥10℃的5日滑動平均組中的最后一個≥10℃的日期作為終止日期，≥10℃活動積溫為起止日期時段內≥10℃的日平均氣溫的總和，采用AUSPLIN氣象插值軟件對研究區內各站點有效積溫進行插值。參照國家標準劃定的不同熟性水稻適宜種植區活動積溫指標，采用ARCGIS軟件中重分類（Reclassify）工具對氣象插值結果進行分類，將東北地區分為早熟型、中熟型、晚熟型水稻適宜種植區以及不適宜水稻種植區4個類別，其對應的≥10℃積溫區間分別為1800～2400℃·d、2400～2850℃·d、≥2850℃·d和＜1800℃·d。

東北地區水稻熟性分區如圖2所示。由圖可見，1986-1995年（T1時段），不適宜種植區主要分布在黑龍江省西北部大興安嶺、黑河及伊春地區和吉林省東部白山、延邊地區，早熟區主要位于黑龍江省東部、吉林省中東部和遼寧省東北部，中熟區主要集中在黑龍江省西南部和吉林省中西部，晚熟區主要分布在遼寧省中西部；與T1時段相比，T2時段（1996-2005年）不適宜種植區范圍向西北部減小，減小部分主要來自黑龍江省西北部和吉林省東部不適宜種植區向早熟區的轉化，早熟區范圍減小明顯，減小部分主要來自黑龍江省東部早熟區向中熟區的轉化，中熟區范圍變化較大，位于黑龍江省東部的中熟區范圍向東北部增加，而位于吉林省中西部的中熟區范圍減小，晚熟區范圍增加明顯，增加部分主要來自吉林省中西部的中熟區向晚熟區的轉化；與T2時段相比，T3時段（2006-2015年）不適宜種植區和早熟區范圍持續向西北部減小，不適宜種植區縮至黑龍江省西北部大興安嶺地區，黑龍江省東北部早熟區持續向中熟區轉化，中熟區范圍因黑龍江省西北部早熟區的轉化而增加，晚熟區范圍則變化不明顯。

圖2 1986-2015年不同階段東北地區不同熟性水稻適宜種植區分布

可見，隨著時間的發展，東北地區水稻不適宜種植區范圍逐漸減小，不同熟性水稻適宜種植區范圍不斷北移東擴，其中以晚熟區擴大范圍最為明顯。

2.2 不同熟性水稻適宜種植區域障礙型冷害影響范圍及其變化

2.2.1 總冷害站次比

由圖3可見，研究期內歷年不同熟性水稻適宜種植區發生障礙型冷害的站次比波動很大，各熟性區存在差異。具體來看，30a間，早熟區有2a（1994、2000年）、中熟區有3a（1997、2008、2013年）、晚熟區有5a（1991、1994、2007、2008、2014年）未發生冷害（總冷害站次比為0），早熟區最少、晚熟區最多，具體年份也不相同。比較不同熟性區冷害影響站次比較大年份可見，早熟區冷害影響最大年份是2006年和2009年，有60%的站點受到冷害影響；中熟區冷害影響最大年份也是2006年，有71%的站點受到影響；晚熟區在2006年有37%的站點受到冷害影響，冷害影響范圍最大的是1989年，有55%的站點受到影響。從《氣象災害大典》和《中國氣象災害年鑒》的記載和姜麗霞等[20-22]的研究中可發現，1993、1994、2007、2008、2013和2014年為東北地區冷害較少年份，而在2006年東北地區遭遇了大范圍階段性低溫，當年黑龍江省有29.6萬hm2農田遭受低溫冷害，吉林省有23個縣（市）出現了障礙型冷害，遼寧省7月下旬北部開原、本溪、阜新等地出現階段性低溫冷害。可見，所選指標能較好地反映整個東北地區范圍內水稻障礙型低溫冷害的發生狀況?？傮w上看，在所有研究年份均有冷害發生，早熟區在T1-T2-T3時段內冷害發生年數變化不明顯，中熟區和晚熟區在T3時段進入冷害少發階段。早熟區和中熟區冷害影響范圍在T3時段呈增加趨勢，晚熟區冷害影響范圍在T3時段呈減小趨勢。

不同熟性區冷害平均影響站次比也不同，早熟區年平均冷害站次比為0.24，中熟區為0.13，晚熟區為0.15，可見，早熟區受冷害影響范圍最大，中熟區最小。采用線性回歸方法分析不同熟性區發生障礙型冷害站次比與研究年份的關系，得出早熟區、中熟區和晚熟區內兩者回歸系數的P值分別為0.86、0.41和0.68，均未通過0.05水平的顯著性檢驗，可見，不同熟性區障礙型冷害站次比隨研究年份的線性變化趨勢不顯著，但具有顯著的階段性特征。分階段計算不同熟性區冷害站次比年平均值可見（圖3），早熟區冷害站次比年平均值存在T1-T2-T3時段歷經先減小后增大的過程，在T3時段達到最大值；中熟區冷害站次比年平均值在T1-T2-T3時段逐漸增大，并在T3時段達到最大值；晚熟區冷害站次比年平均值在T1-T2-T3時段歷經先增大后減小的過程，最大值出現在T2時段。

圖3 1986-2015年東北地區不同熟性水稻適宜種植區障礙型冷害總站次比

注：1986-1995年為T1時段，1996-2005年為T2時段，2006-2015年為T3時段。下同。

Note: T1 stage is 1986-1995, T2 stage is 1996-2005, T3 stage is 2006-2015. The same as below.

2.2.2 不同強度等級冷害站次比

由圖4a可見，研究期內，早熟區發生輕度冷害的年份數最多，中度冷害次之，重度冷害最少，分別為25a、21a和17a。比較不同熟性區冷害影響站次比較大年份可見，早熟區輕度冷害影響最大的是2001年，31%的站點受到影響，中度冷害影響范圍最大的是2006年，27%的站點受到影響，重度冷害影響范圍最大的是2009年，33%的站點受到影響。早熟區各等級冷害年平均站次比也不同，輕度冷害年平均站次比為0.10，中度冷害為0.08，重度冷害為0.06，可見，早熟區輕度冷害的影響范圍最大，重度冷害的影響范圍最小。從早熟區各等級冷害站次比年平均值看出，輕度冷害站次比年平均值在T1-T2-T3時段逐漸增大；中度冷害站次比年平均值在T1-T2-T3時段經歷先減小后增大的過程，在T3時段達到最大值；重度冷害站次比年平均值在T1-T2-T3時段經歷先減小后增大的過程，最大值出現在T1時段。總體上看，早熟區發生冷害年份占所有研究年份的93%，輕度冷害和中度冷害在T1-T2-T3時段內發生年數變化不明顯，重度冷害在T3時段進入少發階段。進入T3時段，各冷害影響范圍呈增加趨勢。

圖4b顯示，過去30a中，中熟區發生輕度冷害年數最多，為25a，重度冷害發生年數最少，為15a，中度冷害年數為20a。中熟區各等級冷害在2006年影響范圍達到最大，24%的站點受到輕度冷害影響，27%的站點受到中度冷害影響，20%的站點受到重度冷害影響。中熟區輕度冷害年平均站次比最大，其次為中度冷害，重度冷害年平均站次比最小，分別為0.06、0.05和0.02。中熟區各等級冷害站次比年平均值顯示，輕度冷害站次比年平均值在T1-T2-T3時段經歷先增大后減小的過程，中度冷害和重度冷害站次比年平均值在T1-T2-T3時段逐漸增大。總體上看，中熟區發生冷害年份占到所有研究年份的90%，輕度冷害和中度冷害在T1-T2-T3時段內發生年數變化不明顯，重度冷害在T1-T2-T3時段內發生年數逐漸增加。進入T3時段，輕度冷害影響范圍減小，中度冷害和重度冷害影響范圍在研究期內逐漸擴大。

圖4c顯示，研究期內晚熟區各等級冷害發生年數差異較大，輕度冷害發生年數為24a，中度冷害發生年數為19a，重度冷害發生年數為13a，輕度冷害最多，重度冷害最少。晚熟區輕度冷害在1989年影響范圍達到最大，有48%的站點受到冷害影響，中度冷害和重度冷害影響范圍最大的是1996年，分別為29%和27%。晚熟區輕度冷害年平均站次比為0.12，影響范圍最大，重度冷害年平均站次比為0.03，影響范圍最小，中度冷害年平均站次比為0.08。分階段計算晚熟區各等級冷害站次比年平均值可見，輕度冷害站次比年平均值在T1-T2-T3時段逐漸減小，中度冷害和重度冷害站次比年平均值在T1-T2-T3時段經歷先增大后減小的過程?？傮w上看，晚熟區發生冷害年份占到所有研究年份的83%，T2時段為輕度冷害高發期，中度冷害在T1-T2-T3時段內發生年數變化不明顯，重度冷害在T3時段進入少發階段。進入T3時段，各等級冷害影響范圍呈減小趨勢。

圖4 1986-2015年東北地區早熟(a)、中熟(b)、晚熟區(c)水稻輕度(1)、中度(2)、重度(3)障礙型冷害站次比

2.3 不同熟性水稻適宜種植區障礙型冷害發生頻率及其變化

2.3.1 總冷害頻率

由圖5可見，整個分析期各站點障礙型冷害的頻率變化較大，不同熟性區存在差異。

圖5a為T1時段水稻障礙型冷害頻率的空間分布。由圖可見，T1時段冷害發生頻率>30%的站點主要位于早熟區黑龍江省西北部孫吳、五營地區、中東部嘉蔭、海倫地區和吉林省東部羅子溝地區，中熟區在黑龍江省東北部龍江、富裕地區，晚熟區位于遼寧省南部熊岳和凌源地區，其它地區冷害較少發生。T1時段為早熟區冷害高發期，冷害發生頻率>30%的站點占29%，其中冷害發生頻率>50%的站點位于早熟區黑龍江省西北部孫吳、五營地區；由圖5b可見，與T1時段相比，T2時段早熟區冷害發生程度減輕，早熟區冷害頻率>30%的站點占16%，主要位于黑龍江省西北部呼瑪地區，中熟區和晚熟區冷害發生程度加重，冷害頻率>30%的站點占比增加，其中，中熟區吉林省東部和龍地區冷害發生頻率>50%；圖5c顯示，T3時段早熟區黑龍江省西北部呼瑪、五營地區發生冷害頻率較高，為50%以上，占站點總數的20%，除此之外，T3時段發生冷害頻率>30%的站點主要位于中熟區黑龍江省東北部北安、樺川地區，晚熟區，遼寧省南部葫蘆島、蓋州地區，其它地區冷害較少發生。

圖5 1986-1995年(a)、1996-2005年(b)和2006-2015(c)年東北地區水稻障礙型冷害發生頻率

2.3.2 不同強度等級冷害頻率

（1）輕度等級冷害。由圖6a可見，T1時段早熟區黑龍江省西北部克山、依安地區，晚熟區遼寧省中部雙遼、法庫地區的站點發生輕度冷害頻率>20%，其它地區冷害較少發生，T1時段為晚熟區輕度冷害發生頻率最高的時段，發生輕度冷害頻率>20%的站點占9%（圖6a1）；T2時段早熟區輕度冷害比T1時段發生程度加重明顯，有8%的站點發生輕度冷害的頻率>30%，T2時段輕度冷害發生頻率>20%的站點主要位于早熟區吉林省東北部安圖地區，中熟區黑龍江省東北部鶴崗區和吉林省東部通化、和龍地區，晚熟區遼寧省東南部岫巖地區，其它地區冷害較少發生（圖6a2）；T3時段早熟區黑龍江省西北部呼瑪地區和晚熟區遼寧省東南部大石橋地區的站點發生輕度冷害頻率>20%，而其它地區發生冷害頻率較低，其中中熟區大部發生輕度冷害頻率在10%～20%（圖6a3）。

（2）中度等級冷害。圖6b顯示，T1時段發生中度冷害頻率>20%的站點主要集中在早熟區黑龍江省西北部五大連池地區和東北部嘉蔭地區，晚熟區遼寧省南部喀左、撫順地區，其它地區冷害發生較少（圖6b1）；T2時段為早熟區中度冷害發生頻率最高的時段，有8%的站點發生中度冷害的頻率>30%，主要分布在吉林省東部安圖地區，發生中度冷害頻率>20%的站點主要位于黑龍江省東北部寶清地區和吉林省東部汪清地區，晚熟區遼寧省南部熊岳和海城地區（圖6b2）；進入T3時段后，整個研究區內發生中度冷害的頻率整體減少，但中熟區發生中度冷害頻率>20%站點數相對T2時段增加，其余大部地區中度冷害較少發生（圖6b3）。

（3）重度等級冷害。由圖6c可見，T1時段早熟區重度冷害發生頻率較高，10%的站點發生重度冷害頻率>20%，其中黑龍江省西北部孫吳、伊春地區的站點發生中度冷害頻率>30%，而中熟區和晚熟區大部分站點發生冷害頻率均≤20%（圖6c1）；T2時段與T1時段相比，各熟性區發生重度冷害程度減輕，僅早熟區吉林省東部和龍地區的站點發生重度冷害頻率在20%～30%，其余地區發生冷害頻率均≤20%（圖6c2）；T3時段早熟區發生重度冷害程度繼續減輕，各熟性區發生重度冷害頻率均≤20%（圖6c3）。

圖6 1986-1995年(1)、1996-2005年(2)和2006-2015年(3)東北地區水稻輕度(a)、中度(b)和重度(c)障礙型冷害發生頻率

3 結論與討論

3.1 結論

（1）1986-2015年東北地區不適宜水稻種植的區域逐漸減小，不同熟性適宜種植區范圍呈現北移東擴的趨勢，適宜水稻種植區域范圍增加。研究期內東北水稻障礙型冷害在不同熟性種植區的時空分布表現為，早熟區冷害年平均站次比最高，晚熟區次之，中熟區最低。

（2）2006-2015年為早熟區冷害多發時段，1986-1995年為中熟區冷害多發時段，1996-2005年為晚熟區冷害多發時段；在2006-2015年，早熟區重度冷害發生頻率最低，中熟區中度冷害發生頻率增加明顯，晚熟區輕度冷害發生頻率顯著增加。

（3）2006年以后，中熟區和晚熟區冷害發生頻率整體下降，但是個別區域在某些年份仍出現不同程度的冷害，表明東北地區冷害發生具有不確定性。不同熟性區冷害多發時間段并不一致，不同等級冷害的多發時間段也不相同。

3.2 討論

（1）本研究在關鍵生育期時間確定上，利用霍普金斯物候定律[17]推算東北地區無觀測數據站點的水稻關鍵生育期時間，與實際生育期時間存在一定偏差。因此，今后應對研究區各站點生育期資料進行完善，使指標能更好地體現冷害對不同熟性區水稻的影響。

（2）由于冷害依據的分級標準、分析方法等不盡相同，本研究所分析的水稻冷害分布區域、頻率與其它文獻[23-24]有差別，主要體現在輕度和中度冷害分布區，而重度冷害的分布區域、頻率與余會康等[25-27]的研究結果大體相同。

（3）東北地區不同熟性水稻適宜種植區遭受冷害的頻率差異較大、波動性明顯，而近50a（1961-2010年）極端低溫事件呈增加趨勢[28]，其波動性導致冷害發生不確定性增加，因此，障礙型冷害仍是東北地區水稻生產中應著重關注的氣象災害。在氣候變化背景下，應結合障礙型冷害發生的表現特征，合理安排種植制度和提前準備防災救災措施[29]，以減小障礙型冷害對水稻產量的影響。

[1] 程葉青.東北地區中低產田改造的區域模式與對策措施[J].干旱區資源與環境,2010,24(11):120-124.

Cheng Y Q, Regional modes and countermeasure on the improvement of medium-productivity farmland in northeastern China[J].Journal of Arid Land Resources and Environment,2010,24(11):120-124.(in Chinese)

[2] 方修琦,王媛,朱曉禧.氣候變暖的適應行為與黑龍江省夏季低溫冷害的變化[J].地理研究,2005,(5):664-672.

Fang X Q, Wang Y, Zhu X X. Change of cool summer hazard under an adaptation behavior to the climate warming in Heilongjiang Province, Northeast China[J].Geographical Research,2005,(5):664-672.(in Chinese)

[3] 蔡志歡,張桂蓮.水稻低溫冷害研究進展[J].作物研究,2018,32(3):249-255.

Cai Z H, Zhang G L. Research progress of low temperature and chilling injury of rice[J].Crop Research,2018,32(3): 249-255.(in Chinese)

[4] 王彤彤,王連敏.我國寒地水稻障礙型冷害研究進展[J].自然災害學報,2013,22(4):167-174.

Wang T T, Wang L M, Advance in research of sterile-style cooling injury of rice in cold region of China[J].Journal of Natural Disasters,2013,22(4):167-174.(in Chinese)

[5] 胡春麗,李輯,焦敏,等.東北地區水稻障礙型低溫冷害變化對區域氣候增暖的響應[J].氣象科技,2015,43(4):744-749.

Hu C L, Li J, Jiao M, et al. Characteristics of severe cold disasters and its responses to regional climate warming in Northeast China[J].Meteorological Science and Technology, 2015,43(4):744-749.(in Chinese)

[6] 張軼.種植范圍變化下的東北水稻低溫冷害風險及其指數保險研究[D].北京:中國農業科學院,2018:26-32.

Zhang Y. The risk of rice chilling injury and index insurance in Northeast under the variation of planting area[D].Beijing: CAAS,2018:26-32.(in Chinese)

[7] 馮喜媛,郭春明,陳長勝,等.基于氣象模型分析東北三省近50年水稻孕穗期障礙型低溫冷害時空變化特征[J].中國農業氣象,2013,34(4):462-467.

Feng X Y, Guo C M, Chen C S, et al. Spatial-temporal variation of sterile-type chilling damages at rice booting stage in northeast China in last 50 years based on meteorological model[J].Chinese Journal of Agrometeorology,2013,34(4): 462-467.(in Chinese)

[8] 王品,魏星,張朝,等.氣候變暖背景下水稻低溫冷害和高溫熱害的研究進展[J].資源科學,2014,36(11):2316-2326.

Wang P, Wei X, Zhang C, et al. A review of cold injury and heat damage to rice growth under global warming[J]. Resources Science,2014,36(11):2316-2326.(in Chinese)

[9] 馬樹慶,馬力文,襲祝香,等.GB/T34967-2017北方水稻低溫冷害等級[S].北京:中華人民共和國國家質量監督檢驗檢疫總局, 中華人民共和國國家標準化管理委員會,2017.

Ma S Q, Ma L W, Xi Z X, et al.GB/T34967-2017 grade of rice cold damage in northern China[S].Beijing: General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China,2017.

[10] 唐亮,吳東立,苗微,等.東北地區食物安全可持續發展戰略研究[J].中國工程科學,2019,21(5):19-27.

Tang L, Wu D L, Miao W, et al. Sustainable development of food security in Northeast China [J].Strategic Study of CAE,2019,21(5):19-27.(in Chinese)

[11] 陳浩,李正國,唐鵬欽,等.氣候變化背景下東北水稻的時空分布特征[J].應用生態學報,2016,27(8):2571-2579.

Chen H, Li Z G, Tang P Q, et al. Rice area change in Northeast China and its correlation with climate change[J]. Chinese Journal of Applied Ecology, 2016,27(8): 2571-2579.(in Chinese)

[12] 孫曉巍,蔣大凱,王明華,等.遼寧省水稻低溫冷害發生規律及天氣學成因[J].現代農業科技,2016,(14):228-229.

Sun X W, Jiang D K, Wang M H, et al. Distribution of rice chilling damages and its meteorological factors in Liaoning Province[J].Modern Agricultural Science and Technology, 2016,(14):228-229.(in Chinese)

[13] 閆麗莉,溫少妍,高文晶,等.整點氣溫缺測的插補方法研究及其初步應用[J].震災防御技術,2019,14(2):446-455.

Yan L L, Wen S Y, Gao W J, et al. Interpolating method for missing data of integral point temperature and its preliminary application[J].Technology for Earthquake Disaster Prevention,2019,14(2):446-455.(in Chinese)

[14] 王海軍,涂詩玉,陳正洪.日氣溫數據缺測的插補方法試驗與誤差分析[J].氣象,2008,(7):83-91.

Wang H J, Tu S Y, Chen Z H. Interpolating method for missing data of daily air temperature and its error analysis[J]. Meteorological Monthly,2008,(7):83-91.(in Chinese)

[15] 宛敏渭,劉秀珍.物候觀測方法[M].北京:科學出版社, 1979:1-17.

Wan M W, Liu X Z. Phenology observation method[M]. Beijing: Science Press,1979:1-17.(in Chinese)

[16] 曲輝輝,姜麗霞,王冬冬.氣候變化對黑龍江省水稻障礙型冷害的影響[J].生態學報,2016,36(3):769-777.

Qu H H, Jiang L X, Wang D D. Influence of climate change on sterile-type cooling injury in rice in Heilongjiang Province, China[J].Acta Ecological Sinica,2016,36(3): 769-777.(in Chinese)

[17] Hopkins A D. Periodic events and natural laws as guides to agricultural research and practice[J].U.S. Dept. Agr. Monthly Weather Rev,1918,9(Suppl) :10-12.

[18] 胡春麗,李輯,林蓉,等.東北水稻障礙型低溫冷害變化特征及其與關鍵生育期溫度的關系[J].中國農業氣象, 2014,35(3):323-329.

Hu C L, Li J, Lin R, et al. Characters of sterile-type chilling damage in rice and its relationship with temperature at key growth stage in Northeast China[J].Chinese Journal of Agrometeorology,2014,35(3):323-329.(in Chinese)

[19] 紀瑞鵬,于文穎,馮銳,等.寒地水稻障礙型冷害指數構建及應用:以遼寧省為例[J].地理科學進展,2017,36(4):437-445.

Ji P R, Yu W Y, Feng R, et al. Construction and application of Rice Sterile-type Cold Damage Index in cold area: a case study of Liaoning Province[J].Progress in Geography,2017, 36(4):437-445.(in Chinese)

[20] 姜麗霞,李帥,閆平,等.黑龍江水稻孕穗期障礙型冷害及其對產量的影響[J].中國農業氣象,2009,30(3):463-468.

Jiang L X, Li S, Yan P, et al. Sterile type cool injury of rice during booting stage and its impacts on rice yield in Heilongjiang province[J]. Chinese Journal of Agrometeorology,2009,30(3):463-468. (in Chinese).

[21] 張海娜,李菲,戴欣童.1951-2010年東北地區玉米及水稻低溫冷害統計分析[J].現代農業科技,2015,(3):264-266.

Zhang H N, Li F, Dai X T. Chilling injury statistical analysis of corn and rice in Northeast China from 1951 to 2010[J].Modern Agricultural Science and Technology, 2015,(3):264-266.(in Chinese)

[22] 紀瑞鵬,于文穎,武晉雯,等.氣候平均值變化對遼寧水稻延遲型冷害評估結果的影響[J].應用生態學報,2015,26(6): 1695-1703.

Ji P R, Yu W Y, Wu J W, et al. Effect of climatic mean value change on the evaluation result of rice delayed cold damage in Liaoning Province, Northeast China[J].Chinese Journal of Applied Ecology,2015,26(6):1695-1703.(in Chinese)

[23] 紀仰慧,王晾晾,姜麗霞,等.黑龍江省2009年水稻障礙型冷害評估[J].氣象科技,2011,39(03):374-378.

Ji Y H, Wang L L, Jiang L X, et al. Evaluation of sterile-type rice cold injury in Heilongjiang Province in 2009[J]. Meteorological Science and Technology,2011, 39(03):374-378.

[24] 王小一.東北地區水稻低溫冷害時空分布特征與災變規律分析[D].沈陽農業大學,2019.

Wang X Y. Temporal and spatial distribution characteristics and catastrophic laws of rice chilling in Northeast China[D]. Shenyang Agricultural University,2019.

[25] 余會康,郭建平.氣候變化下東北水稻冷害時空分布變化[J].中國生態農業學報,2014,22(5):594-601.

Yu H K, Guo J P. Variation in spatial and temporal distribution of chilling injury of rice under climate change in Northeast China[J].Chinese Journal of Eco-Agriculture, 2014,22(5):594-601.(in Chinese)

[26] 姜麗霞,李帥,申雙和,等.近46a黑龍江水稻障礙型冷害及其與氣候生產力的關系[J].大氣科學學報,2010,33(3): 315-320.

Jiang L X, Li S, Shen S H, et al. Sterile-type cool injury on rice and its relationship with climate productivity in Heilongjiang Province from 1961 to 2006[J].Transactions of Atmospheric Sciences,2010,33(3):315-320.(in Chinese)

[27] 李雨鴻,張微瑋,李琳琳,等.遼寧水稻低溫冷害分布特征及危險性評估[J].大麥與谷類科學,2018,35(5):47-52.

Li Y H, Zhang W W, Li L L, et al. Identification of the distribution patterns of chilling injury in rice in Liaoning Province and its risk assessment[J].Barley and Cereal Sciences,2018,35(5):47-52.(in Chinese)

[28] 王艷華,任傳友,韓亞東,等.東北地區活動積溫和極端持續低溫的時空分布特征及其對糧食產量的影響[J].農業環境科學學報,2011,30(9):1742-1748.

Wang Y H, Ren C Y, Han D Y, et al. The tempo-spatial patterns of active accumulated and consecutive extreme low temperature and their impacts on grain crop yield in Northeast China[J].Journal of Agro-Environment Science, 2011,30(9):1742-1748.(in Chinese)

[29] 王媛,方修琦,徐錟,等.氣候變暖與東北地區水稻種植的適應行為[J].資源科學,2005,(1):121-127.

Wang Y, Fang X Q, Xu Y, et al. Impact of climate warming and adaptation activities of rice plantation in Northeast China[J].Resources Science,2005,(1):121-127.(in Chinese)

Temporal and Spatial Distribution of Sterile-type Chilling Damage in Rice Growth- suited Areas with Different Maturity in Northeast China

LIANG Li-jiang1, WU Yong-feng1, LIU Cong1,2, SHAO Shuai3

(1. Key Laboratory of Agricultural Environment, Ministry of Agriculture/Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081, China; 2. College of Earth Sciences, East China University of Technology, Nanchang 330000; 3.Department of Computer Science, Hangzhou Dianzi University, Hangzhou 310000)

Northeast China is the main rice producing area in China, and chilling damage is the most serious meteorological disaster affecting rice production in this area. Study on the temporal and spatial distribution of sterile-type chilling damage in the rice growth-suited areas with different maturity is helpful to provide the scientific basis for the rational production and management of rice in Northeast China. Based on the daily mean air temperature data at 167 meteorological stations and rice growth stage record data at 28 phenological observation stations, Chilling Damage Grade of Rice Northern China (GB / t 34967-2017) was selected to evaluate sterile-type chilling damage in the rice growth-suited areas with different maturity in Northeast China. First, the rice growth-suited areas with different maturity at three stages 1986-1995, 1996-2005 and 2006-2015 were determined based on the 80% guarantee rate of active accumulated temperature. Then, the ratio (F) of the meteorological stations with sterile-type chilling damage and frequency (P) of chilling damage were calculated to analyze temporal and spatial distribution of the damage from 1986 to 2015. The results showed that the area unsuitable for rice growth decreased gradually and the rice growth-suited area with different maturity showed an expanding trend towards north and east. Among them, the late-maturating rice area showed the most significant increase. The ratio (F) of the stations with sterile-type chilling damage was the largest in the early-maturating rice area, followed by the late-maturating rice area and midseason-maturating rice area. The occurrence frequency (P) of sterile-type chilling damage varied greatly among stations and rice growth-suited areas with different maturity. During 1986-1995, 1996-2005 and 2006-2015, obvious sterile-type chilling damage mainly occurred in the midseason-maturating, late-maturating and early-maturating rice area, respectively. During 2006-2015, the frequency (P) of sterile-type chilling damage with Sever grade reduced in the early-maturating rice area, while the frequency (P) with Middle grade in the midseason-maturating rice area and Slight grade in the late-maturating rice area increased. It should be noted that sterile-type chilling damage in Northeast China was still uncertain despite the overall decline in frequency, and probably occurred in any area and any year. The fluctuation of sterile-type chilling damage led to the increase of its uncertainty. In the context of climate change, the planting system, disaster prevention and relief measures should be reasonably arranged in combination with the characteristics of sterile-type chilling damage occurrences, so as to reduce its impact on rice yield.

Rice; Maturity; Active accumulated temperature; Sterile-type chilling damage of rice; Northeast China

10.3969/j.issn.1000-6362.2020.05.005

梁立江,武永峰,劉聰,等.東北地區不同熟性水稻適宜種植區障礙型冷害時空變化[J].中國農業氣象,2020,41(5):308-319

2019?11?25

武永峰，E-mail: wuyongfeng@caas.cn

十三五國家重點研發計劃子課題（2017YFD0300402-2）

梁立江，E-mail：llj1943@qq.com