近30年氣候異常對江蘇省褐飛虱災變性遷入的影響

包云軒，蔣 蓉，謝曉金，朱葉芹，楊榮明,朱 鳳

(1. 江蘇省農(nóng)業(yè)氣象重點實驗室，南京信息工程大學， 南京 210044； 2. 江蘇省植物保護站, 南京 210013)

近數(shù)十年來，全球氣候變化對農(nóng)業(yè)生態(tài)系統(tǒng)產(chǎn)生了重大的影響，對農(nóng)作物害蟲而言，也不例外。為此，研究氣候變化背景下害蟲的發(fā)生規(guī)律，及時有效地控制其危害，對減輕農(nóng)作物產(chǎn)量的損失，確保我國糧食安全是至關重要的。褐飛虱，Nilaparvatalugens(St?l)，是危害水稻的重要害蟲之一，具群聚性、遠距離遷飛性和災變突發(fā)性，它廣泛分布于中國、日本、朝鮮半島、東南亞、太平洋島嶼和澳大利亞等國家和地區(qū)[1]。近30多年來，受全球氣候變化和農(nóng)業(yè)產(chǎn)業(yè)結構調(diào)整等因素的影響，褐飛虱的遷入和致災趨于復雜化，出現(xiàn)了持續(xù)大發(fā)生的新特點，給水稻高產(chǎn)、穩(wěn)產(chǎn)帶來了嚴重威脅[2- 3]。

海溫異常變化導致海氣相互作用異常，繼而引發(fā)大氣環(huán)流異常，出現(xiàn)異常的天氣、氣候變化[4- 5]。ENSO是厄爾尼諾(El Nio)現(xiàn)象和南方濤動(Southern Oscillation)的總稱，它是反映赤道太平洋海溫和中、低緯海平面氣壓場變化的強信號，也是影響我國氣候異常的重要原因之一[6- 12]。在氣候變化研究中常用尼諾指數(shù)反映海溫場的異常，用南方濤動指數(shù)來反映海氣相互作用的異常。ENSO事件的發(fā)生可以通過改變大氣環(huán)流或天氣、氣候而影響害蟲種群的發(fā)生，而遷飛性害蟲的遷入和危害與天氣、氣候關系更為密切[13- 14]。

天氣、氣候條件是決定褐飛虱發(fā)生和種群發(fā)展的關鍵因子，由于異常的天氣、氣候條件常導致異常的褐飛虱發(fā)生時間、發(fā)生區(qū)域和危害程度。溫度是褐飛虱種群生命過程中的一個重要的環(huán)境因子，它可通過對發(fā)育速率、若蟲存活率、成蟲生殖力和空中遷飛層的影響而顯著改變種群的增長能力?？罩羞w飛的褐飛虱種群的降落主要靠降水和下沉氣流的動力作用。褐飛虱生長發(fā)育和繁殖的適宜溫度在20—30℃，最適溫度為26—28℃，空中遷飛的適宜溫度層為15—22℃；當溫度高于30℃或低于20℃，對它的生長發(fā)育和繁殖不利，溫度低于15℃或高于22℃的大氣層不利于褐飛虱的遠距離遷飛。褐飛虱生長發(fā)育、繁殖和遠距離遷飛的最適相對濕度為80%以上。“盛夏不熱，晚秋不冷，夏秋多雨”的氣候有利于褐飛虱的發(fā)生[15- 19]。

本文擬通過對近28年江蘇省褐飛虱的發(fā)生程度與海溫場、海平面氣壓場、海氣相互作用、氣候要素關系的分析，揭示ENSO事件和氣候異常對該省褐飛虱災變性遷入的影響，為江蘇省褐飛虱遷入量的中長期預報提供科學依據(jù)。

1 資料和方法

1.1 研究區(qū)概況

江蘇省位于亞洲大陸東岸的中緯度地帶，屬東亞季風氣候區(qū)，處在亞熱帶和暖溫帶氣候過渡地帶。全省地勢較為平坦，介于30° 46′—35° 07′ N，116° 22′—121° 55′ E之間，一般以淮河、蘇北灌溉總渠一線為界，以北地區(qū)屬暖溫帶濕潤、半濕潤季風氣候；以南地區(qū)屬亞熱帶濕潤季風氣候。江蘇擁有近1000 km長的海岸線，海洋對江蘇氣候有著顯著的影響。由于地處中緯度的海陸相過渡帶和氣候過渡帶，兼受中緯度西風帶和低緯度信風帶(東風帶)天氣系統(tǒng)的影響，在太陽輻射、大氣環(huán)流以及江蘇特定的地理位置、地貌特征的綜合影響下，江蘇省的基本氣候特點是：氣候溫和、季風顯著、春溫多變、夏秋多雨、雨量適中、光熱充沛。在這種良好的自然地理和氣候條件下，全省種植的糧食作物以水稻、小麥為主。而水稻種植制度經(jīng)過多年的優(yōu)勝劣汰，已穩(wěn)定在單季稻上，其主要生長季正與褐飛虱遷入長江中下游稻區(qū)的盛行期相遇，一直受到褐飛虱災變性遷入的威脅，歷年水稻產(chǎn)量也因此而產(chǎn)生明顯的增減波動。

1.2 資料來源及處理

1.2.1 褐飛虱蟲情資料

本文從江蘇省不同區(qū)域篩選出資料比較完整的1983—2010年高郵、通州和宜興3個測報站作為代表站，高郵站代表蘇北地區(qū)、通州站代表蘇中地區(qū)、宜興站代表蘇南地區(qū)，這些站的蟲情實測資料由江蘇省植物保護站提供。統(tǒng)計高郵、通州和宜興三站的歷年褐飛虱燈下累計蟲量時，從燈下始見日開始(通常在6月)至終見日(通常在10月)結束。

1.2.2 氣象資料

氣象資料來自國家氣象信息中心，包括1981—2010年全國各氣象臺站逐日氣溫和降水觀測資料。

1.2.3 ENSO指數(shù)數(shù)據(jù)

本文引用了兩個目前世界氣候變化研究領域最常用的ENSO指數(shù)：ONI指數(shù)和SOI指數(shù)，指數(shù)年限為1981—2010年。從熱力學角度判斷El Nio(厄爾尼諾)事件時，采用的是美國國家海洋大氣管理局(National Oceanic and Atmospheric Administration, NOAA)氣候預測中心(Climate Predicting Center, CPC)在業(yè)務上使用的海洋尼諾指數(shù)(Oceanic Nino Index，簡稱ONI)，它是根據(jù)NOAA擴展重建的海面溫度資料(Extended Reconstructed Sea Surface Temperatures，ERSST)計算的Nino3區(qū)和Nino4區(qū)海表溫度距平的3個月滑動平均值。ONI≥+0.5℃持續(xù)6個月以上時稱為1次El Nino(厄爾尼諾)事件，ONI≤-0.5℃持續(xù)6個月以上時稱為1次La Nina(拉尼娜)事件；ONI的變化也反映了海溫場和海氣相互作用對全球氣候變化的影響[20]。南方濤動指數(shù)(Southern Oscillation Index，SOI)，用來度量南太平洋副熱帶高壓與印度洋赤道低壓間氣壓變化的負相關關系，當南太平洋副熱帶高壓比常年增高(降低)時，印度洋赤道低壓就比常年降低(增高)，SOI為負值(正值)。通常用印度洋的塔希堤島(148° 5′ W，17° 53′ S)和南太平洋的達爾文港(130° 59′ E，12° 20′ S)之間的海平面氣壓差來表示南方濤動的振動和位相指數(shù)[21]，該指數(shù)值同樣來源于美國海洋大氣局的氣候預測中心(CPC)；當SOI<0持續(xù)6個月以上時稱為1次El Nino(厄爾尼諾)事件，當SOI>0持續(xù)6個月以上時稱為1次La Nina(拉尼娜)事件；SOI的變化從動力角度反映了海氣相互作用對全球氣候變化的影響。

2 結果與分析

2.1 江蘇省褐飛虱歷年發(fā)生情況

本文根據(jù)江蘇省植物保護站收集和統(tǒng)計的褐飛虱蟲情資料，結合農(nóng)業(yè)部1995年頒發(fā)的稻飛虱單站發(fā)生程度等級標準[22]，將江蘇省各縣級植保站褐飛虱發(fā)生危害程度分為5個等級：輕發(fā)生為1級，日燈誘蟲量小于1000 頭；偏輕發(fā)生為2級，日燈誘蟲量1000—2000 頭；中等發(fā)生為3級，日燈誘蟲量2001—3000 頭；偏重發(fā)生為4級，日燈誘蟲量3001—4000 頭；大發(fā)生為5級，日燈誘蟲量大于4000 頭。

從圖1中可以看出：高郵站1983、1987、1991、1995、2006、2007年是偏重發(fā)生(4級)或大發(fā)生年(5級)，通州站1983、1985、1987、1991、2006年是偏重發(fā)生或大發(fā)生年，宜興站1983、1985、1987、1988、1990、1991、2005、2006年是偏重發(fā)生或大發(fā)生年。3個站1983、1987、1991、2006年同為偏重發(fā)生或大發(fā)生年，這也是江蘇全省范圍內(nèi)褐飛虱發(fā)生較重的4a。

圖1 1983—2010年江蘇省高郵、通州和宜興三站褐飛虱的發(fā)生等級

表1為三站褐飛虱首次遷入峰日期，從中可以看出：褐飛虱的首次遷入峰期基本在6月和7月；高郵站1994、2001和2004年沒有遷入峰，1996、2002、2005和2009年的首次遷入峰期在8月和9月，沒有遷入峰或遷入峰晚的年份褐飛虱發(fā)生等級均為輕發(fā)生；通州站1984、1999、2001、2002和2004年沒有遷入峰，1992、1996、1998和2009年的首次遷入峰期在8月和9月，沒有遷入峰和遷入峰晚的年份褐飛虱發(fā)生等級也為輕發(fā)生和偏輕發(fā)生；宜興站2001、2002、2004和2009年沒有遷入峰，1984和2010年的首次遷入峰期在8月和9月，沒有遷入峰和遷入峰晚的年份褐飛虱發(fā)生等級也為輕發(fā)生和偏輕發(fā)生。

表1 褐飛虱首次遷入峰日期(月-日)

2.2 ENSO事件對褐飛虱遷入的影響

2.2.1海溫場(ONI指數(shù))對江蘇省褐飛虱遷入的影響

厄爾尼諾事件是指赤道中東太平洋每隔幾年發(fā)生的大規(guī)模表層海水持續(xù)(6個月以上)異常偏暖的現(xiàn)象，而把赤道中東太平洋表層海水大規(guī)模持續(xù)(6個月以上)異常偏冷的現(xiàn)象稱為反厄爾尼諾或拉尼娜事件。

本文參照相關文獻[20]和NOAA/CPC在業(yè)務上使用的ONI指數(shù)定義計算得到了ONI指數(shù)從1981年1月至2010年12月的時間變化曲線(見圖2)。選擇指數(shù)大于1.0(小于-1.0)的時間段為強厄爾尼諾(強拉尼娜)事件，大于0.5且小于1.0(大于-1.0且小于-0.5)的時間段為弱厄爾尼諾(弱拉尼娜)事件。從圖中可以看出：1983、1987、1988、1998、2003和2010年都為強厄爾尼諾年；1989、1999、2000和2008年為強拉尼娜年。

圖2 1981年1月—2010年12月ONI變化曲線，虛線為 ±0.5℃

表2為1981—2010年間發(fā)生的強、弱厄爾尼諾事件和強、弱拉尼娜事件的具體情況及其與3個站點褐飛虱發(fā)生程度的對比。結合圖1和表2可以看出：1983—2010年期間，高郵站褐飛虱遷入量較多(發(fā)生程度為4級以上，下同)的年份有6a，其中出現(xiàn)在強厄爾尼諾年的有4a，出現(xiàn)在弱厄爾尼諾年的有1a，還有1a是弱拉尼娜年(2006年)；通州站褐飛虱遷入量較多的年份有5a，其中出現(xiàn)在強厄爾尼諾年的有3a，出現(xiàn)在強拉尼娜年的有1a(1985年)，還有1a是弱拉尼娜年(2006年)；褐飛虱在宜興遷入量較多的年份有8a，其中出現(xiàn)在強厄爾尼諾年的有3a，出現(xiàn)在弱厄爾尼諾年的有1a，出現(xiàn)在強拉尼娜年的有2a，出現(xiàn)在弱拉尼娜年的有1a，還有1a是強厄爾尼諾事件向強拉尼娜事件的過渡年。因此，由海溫場反映出來的結果是：1983年以來褐飛虱偏重以上發(fā)生年(4級以上)出現(xiàn)在厄爾尼諾事件期間的概率為：高郵83.33%、通州60.0%和宜興50%。由此看來，厄爾尼諾事件可以作為褐飛虱偏重以上發(fā)生年的重要氣候背景信息，而ONI指數(shù)則可以作為海溫場預測褐飛虱發(fā)生程度的一個重要因子。

從表1和表2中還可以看出：從1982年以來，厄爾尼諾事件的發(fā)生在較大程度上引發(fā)了褐飛虱的大量遷入，1981—2010年的30a中共發(fā)生了9次厄爾尼諾事件，引發(fā)了江蘇省6個明顯的褐飛虱大量遷入年(4級以上)，褐飛虱遷入的偏重以上發(fā)生(以褐飛虱遷入江蘇的首次遷入峰日期為準)與厄爾尼諾事件開始期之間有1—14個月的滯后時間。

2.2.2 南方濤動對江蘇省褐飛虱遷入的影響

南方濤動主要指東南太平洋的副熱帶高壓與低緯印度洋和西太平洋的赤道低壓之間呈負相關的一種大尺度氣壓升降振蕩現(xiàn)象，其強弱用南方濤動指數(shù)SOI衡量。當SOI指數(shù)為負值時，表明兩站海平面氣壓差比常年偏大，如果SOI連續(xù)6個月小于0，則出現(xiàn)一次厄爾尼諾事件；反之，當SOI指數(shù)為正值時，表明兩站海平面氣壓差比常年偏小，如果SOI連續(xù)6個月大于0，則出現(xiàn)一次拉尼娜事件。圖3給出了SOI指數(shù)從1981年1月至2010年12月的時間變化曲線。選擇指數(shù)小于-2(大于2)的時間段為強厄爾尼諾(或強拉尼娜)事件。從圖中可以看出：1982、1987、1991、1994、1998和2004年都為強厄爾尼諾年；1989、1999和2008年為強拉尼娜年。

表3為由SOI確定的1981—2010年期間發(fā)生的強、弱厄爾尼諾事件和強、弱拉尼娜事件及與3個站點褐飛虱發(fā)生程度的對比。從中可以看出：SOI持續(xù)為負值的厄爾尼諾事件，會引發(fā)江蘇省褐飛虱遷入的偏重以上發(fā)生，且褐飛虱遷入的偏重以上發(fā)生(以褐飛虱遷入江蘇的首次遷入峰高峰日為準)與SOI變化反映出來的厄爾尼諾事件開始期之間有1—14個月的滯后時間。這與前述由海溫場分析的厄爾尼諾事件導致褐飛虱偏重以上發(fā)生的結論完全吻合。統(tǒng)計顯示：1983年以來褐飛虱偏重以上發(fā)生年(4級以上)出現(xiàn)在厄爾尼諾事件期間的概率也為：高郵83.33%、通州60%和宜興50%。

表2 高郵、通州和宜興三站褐飛虱發(fā)生程度與由ONI確定的厄爾尼諾(拉尼娜)事件的關系4

Table 2 Relationship between BPH′s occurrence grades of Gaoyou Station, Tongzhou Station and Yixing Station and El Nio & La Nia event defined by ONI from 1981 to 2010

表2 高郵、通州和宜興三站褐飛虱發(fā)生程度與由ONI確定的厄爾尼諾(拉尼娜)事件的關系4

厄爾尼諾事件/拉尼娜事件(El Ni?o event & La Ni?a event)事件Event出現(xiàn)時間(年-月)Beginning time(year-month)結束時間(年-月)(Ending time) (year-month)持續(xù)時間(月)(Lasting periods(month)強度(Strength)褐飛虱發(fā)生等級(BPH′s occurrence grades)年份Year植保站Plant protection station高郵通州宜興厄爾尼諾事件1982-051983-0614強1983555El Ni?o event1986-081988-0219強19875551991-051992-0614強19914551994-091995-037強19954311997-051998-0412強19981122002-052003-0210強20031212004-072005-017弱20051152006-092007-015弱20075332009-072010-0410強2010312拉尼娜事件1983-091984-015弱1984111La Ni?a event1984-101985-0912強19853551988-051989-0513強19902141995-091996-037弱19961111998-072001-0333強19991122005-112006-035弱20065552007-082008-0611強2008113

圖3 1981年1月—2010年12月SOI變化曲線

綜合考慮ONI和SOI反映的ENSO事件對江蘇省褐飛虱災變性遷入的影響，江蘇省褐飛虱遷入的偏重以上發(fā)生(以褐飛虱首次遷入峰高峰日為準)均出現(xiàn)在厄爾尼諾事件開始期的1—14個月后。

在ENSO事件與褐飛虱遷入偏重以上發(fā)生之間出現(xiàn)這種對應關系的主要原因是：一方面，由于在熱帶西太平洋及印度洋東部常出現(xiàn)多年平均海表溫度(SST)在28℃以上的暖海區(qū)，該海區(qū)被稱為大洋暖池(Warm Pool)，它東西跨越150個經(jīng)度，南北伸展約35個緯度。由于太陽輻射、熱量交換、自東向西信風吹送等的作用，大量暖水逐漸積蓄在暖池區(qū)，致使該區(qū)SST比赤道東太平洋海面高出3—9℃。因此，正常情況下，赤道太平洋SST西高東低，且溫差維持在一定的范圍內(nèi)。當發(fā)生厄爾尼諾事件時，西太平洋暖池的西部SST會出現(xiàn)降溫，但暖池范圍會增大，暖池中心及暖池東部邊界東移，使得赤道東太平洋SST升高；厄爾尼諾事件結束后的半年內(nèi)，暖池SST有較小升溫，暖池范圍逐漸縮小，但仍較常年偏大，暖池中心西移。由厄爾尼諾事件導致的暖池范圍擴大(雖然暖池區(qū)內(nèi)有降溫)其重要的影響結果之一是東亞和南亞地區(qū)近地表氣溫偏高，且這種偏高在冬、春反映尤其明顯。另一方面，在赤道太評洋上空，垂直方向盛行著沃克環(huán)流，它在低空吹偏東風，在西部暖池上方為上升氣流，到高空吹偏西風，東部轉(zhuǎn)為下沉氣流。在ENSO事件發(fā)生期間，沃克環(huán)流減弱，西太平洋哈德萊環(huán)流偏弱。哈德萊環(huán)流的偏弱又造成西北太平洋副高位置偏東且強度偏弱，同時副高西側和西北側的偏南氣流和西南氣流偏強且影響范圍廣。在這種形勢下，來自孟加拉灣和南海的水汽可以源源不斷地向東北或偏北方向輸送，而從西北方或北方高緯地區(qū)南下的冷空氣可以長驅(qū)直入在我國江淮流域或江南地區(qū)與暖濕空氣相遇，發(fā)生對流，從而使這些地區(qū)夏秋遷飛季內(nèi)降水偏多。

表3 高郵、通州和宜興三站褐飛虱發(fā)生程度與由SOI確定的厄爾尼諾(拉尼娜)事件的關系(1981—2010)

Table 3 Relationship between BPH′s occurrence grades of Gaoyou Station, Tongzhou Station and Yixing Station and El Nio & La Nia event defined by SOI from 1981 to 2010

表3 高郵、通州和宜興三站褐飛虱發(fā)生程度與由SOI確定的厄爾尼諾(拉尼娜)事件的關系(1981—2010)

厄爾尼諾事件/拉尼娜事件(El Ni?o event & La Ni?a event)事件Event出現(xiàn)時間(年-月)Beginning time(year-month)結束時間(年-月)(Ending time) (year-month)持續(xù)時間(月)(Lasting periods(month)強度(Strength)褐飛虱發(fā)生等級(BPH′s occurrence grades)年份Year植保站Plant protection station高郵通州宜興厄爾尼諾事件1982-051983-0412強1983555El Ni?o event1986-081988-0219強19875551991-031992-0414強19914551994-071995-028強19954311997-051998-0412強19981122002-052003-0614弱20031212004-062005-029強20051152006-072007-028強20075332009-082010-038弱2010312拉尼娜事件1983-051984-0210弱1984111La Ni?a event1984-111985-0911弱19853551988-071989-0713強19902141995-091996-0711弱19961111998-072001-0333強19991122005-092006-048強20065552007-102008-1013強2008113

褐飛虱是喜溫、喜濕性害蟲，除在海南、兩廣南部及云南南部冬季有少量蟲源存活外，在中國其它大部分地區(qū)常年均不能越冬，初始蟲源主要來源于中南半島[23- 24]。在厄爾尼諾年，中國及東南亞地區(qū)易出現(xiàn)暖冬。冬季偏暖為褐飛虱的生長發(fā)育提供了良好條件，不僅可使褐飛虱越冬范圍擴大、冬季繁殖或殘存數(shù)量增加、蟲源基數(shù)提高，而且還使越南等境外或我國稻飛虱越冬區(qū)的褐飛虱發(fā)育速率加快、遷入和始害期提前、為害期延長[24- 27]。此外，大量研究證實：ENSO發(fā)展年的夏季我國多雨帶位于江淮流域，ENSO恢復年的夏季我國長江及江南地區(qū)雨量偏多[28- 30]。一定強度的降水是脅迫褐飛虱降落的重要動力條件，降水多，常使得過境遷飛性害蟲在雨區(qū)內(nèi)降落。夏季是江蘇省褐飛虱集中遷入(降落)的盛期，夏季降水多，易導致遷入江蘇的蟲量多，繼而導致發(fā)生程度重。因為大氣的變化滯后于海水溫度的變化，所以ENSO事件對江蘇省天氣、氣候的影響有一定的滯后效應，而褐飛虱對天氣、氣候的反應又滯后于天氣、氣候的變化，因此，褐飛虱首次遷入峰出現(xiàn)時間明顯滯后于ENSO事件的開始期。

有研究證實，西太平洋副高的主體和大洋暖池的主體都位于西太平洋熱帶和副熱帶地區(qū)，二者關系密切，一般暖池范圍增大，副高面積會相應增大；暖池強度增強，副高強度也會增強[31]。因此，在厄爾尼諾發(fā)生期間，西太平洋副高的占據(jù)范圍擴大，中心和東部邊界東移，但西部邊界和覆蓋區(qū)少變，強度減弱，對我國大陸的影響減小，而此時西風帶環(huán)流增強，盛行范圍擴大，受其影響，在厄爾尼諾事件開始后的的當年或處于恢復期的次年，春、夏季西南季風開始早、強度大、影響范圍廣，導致當年或次年褐飛虱遷入江蘇的始見期早、遷入峰次多、蟲量大、發(fā)生程度重，如大發(fā)生的1987年和1991年。在厄爾尼諾事件結束后的半年內(nèi)，副高有較小程度的增強，雖然其覆蓋范圍逐漸縮小，但仍較常年偏大，中心西移，主體控制我國東部沿海的福建、浙江、上海和江蘇等省市，導致這些省市受強下沉氣流影響，不斷有褐飛虱遷入種群降落；由于這種增強和西移常是振蕩式(或非連續(xù)性)的。在其振蕩過程中，常在副高主體的西南側、西側和西北側陣發(fā)性地出現(xiàn)東南氣流、偏南氣流和西南氣流，并輸送蟲源北上，為褐飛虱種群在江蘇的降落提供動力，典型的年份有1983年、1995年和2005年。

有些年份雖然是拉尼娜年，但褐飛虱的遷入?yún)s是偏重以上的發(fā)生，如2006年江蘇全省褐飛虱普遍大發(fā)生，1985年江蘇中部和南部偏重以上發(fā)生，1990年江蘇蘇南地區(qū)大發(fā)生。據(jù)李崇銀[32]的研究，ENSO年(El Nino年)西北太平洋發(fā)生的熱帶氣旋總數(shù)偏少，在我國沿海登陸的熱帶氣旋次數(shù)也少；反ENSO年(La Nina年)則相反。楊淑瑞等[33]研究指出：ENSO年西北太平洋熱帶氣旋生成位置偏東偏南，強熱帶氣旋持續(xù)天數(shù)較多，首次熱帶氣旋出現(xiàn)日期較晚，而末次結束日期較早；反ENSO年則相反。已有學者就臺風(或強熱帶氣旋)對褐飛虱遷入影響開展過一系列研究[34- 35]，證實：在多強熱帶氣旋登陸的年份及其影響區(qū)，褐飛虱的遷入峰次多、遷入量大、發(fā)生程度重。對比歷年登陸并影響江蘇的強熱帶氣旋及其影響期間發(fā)生的褐飛虱遷入蟲情資料后發(fā)現(xiàn)：2006年有7個強熱帶氣旋登陸我國，其中5個影響到江蘇大部；1985年有9個強熱帶氣旋登陸我國，其中4個影響到江蘇中部和南部；1990年也有9個強熱帶氣旋登陸我國，其中3個影響到江蘇蘇南地區(qū)，而且這些登陸強熱帶氣旋都發(fā)生在7—9月，表明登陸熱帶氣旋的影響是導致江蘇省褐飛虱災變性遷入的重要因素之一。雖然其它拉尼娜年也有熱帶氣旋登陸并影響我國，但影響到江蘇的個數(shù)、范圍和程度都不如這3a。

還有個別年份如1988年，ONI指數(shù)絕對值很小，年均值為-0.2，但全年ONI值卻變化于-1.95—+1.63之間，是典型的強厄爾尼諾事件向強拉尼娜事件的過渡年。這一年蘇南地區(qū)褐飛虱的災變性遷入很明顯，這是因強厄爾尼諾事件的滯后效應而導致的。事實上，因強厄爾尼諾效應導致的1987—1988年的暖冬和1988年早春的西南季風偏早、偏強北上對越冬蟲源的高基數(shù)和春夏北遷起了關鍵作用，也為這一年蘇南地區(qū)褐飛虱大發(fā)生奠定了基礎。

2.3 氣候異常對褐飛虱遷入量的影響

反映某一地區(qū)氣候異常的關鍵要素是降水和氣溫，本文根據(jù)歷年高郵、通州、宜興3站逐日褐飛虱燈誘資料，結合當?shù)貧鉁?、降水資料，分析了6—10月的降水和氣溫對褐飛虱遷入的影響。

2.3.1 降水對褐飛虱遷入的影響

一個地方的基本降水特征可以從降水量和降水日數(shù)上反映出來。將1983—2010年間6—10月降水日數(shù)與這28a間高郵、通州和宜興三個站點褐飛虱遷入峰次對應起來分析后發(fā)現(xiàn)(圖4)：降水日數(shù)多的年份，褐飛虱遷入峰次多。高郵站褐飛虱遷入峰次達到8次的5a中有4a降水日數(shù)超過平均值，占80%；通州站遷入峰次達到8次的10a中有6a降水日數(shù)超過平均值，占60%；宜興遷入峰次達到8次的10a中有4a年降水日數(shù)達到或超過平均值，占40%。可見，褐飛虱的遷入峰次與降水日數(shù)有關，降水日數(shù)越多，遷入峰次就越多，降水對褐飛虱的遷入(降落)有明顯的動力迫降作用。但從三個站降水日數(shù)對遷入峰次的影響來看，降水日數(shù)多并不是遷入峰次多的唯一原因。實際上，無降水時的強下沉氣流也會對遷入峰次有較大的影響[36]。

圖4 高郵、通州和宜興3站1983—2010年褐飛虱發(fā)生峰次與降水日數(shù)

分析1983—2010年間6—10月降水日數(shù)與這28a間3個站點褐飛虱發(fā)生等級情況可以發(fā)現(xiàn)(圖5)：高郵站褐飛虱偏重以上發(fā)生的6a中有4a降水日數(shù)超過了平均值，通州站褐飛虱偏重以上發(fā)生的5a所有年份降水日數(shù)都超過了平均值，宜興站偏重以上發(fā)生的8a中有4a降水日數(shù)超過了平均值。由此可見，降水日數(shù)的多少對褐飛虱的遷入量有明顯的影響，降水日數(shù)越多，遷入量就越多。但有些年份，褐飛虱遷入量多，降水日數(shù)不一定多，如2006年江蘇全省降水日數(shù)并不多，遷入量卻很多，這一年副高短期振蕩造成的間歇性、高頻次強下沉氣流給褐飛虱的遷入創(chuàng)造了良好的條件[36]。

從圖6給出的3個站點6—10月的降水量距平百分率的年際變化和圖1所提供的遷入蟲量來看，在高郵站褐飛虱偏重或大發(fā)生的6a中有4a降水量偏多，即1983、1987、1991和2007年，占66.67%；通州站5個褐飛虱偏重或大發(fā)生的年份中有3a降水量偏多，即1983、1987和1991年，占60%；宜興站偏重或大發(fā)生的8a中，1985、1987、1990和1991年降水量偏多，占50%。表明降水量多少對褐飛虱的降落有明顯的影響。

對一個植保站而言，某一年份的褐飛虱遷入發(fā)生程度重(4級以上)，有可能由兩種情況導致的：一是遷入量不大(須達到一定的量級)，但遷入峰次多，累積效應比較明顯，全年遷入蟲量達到了4級以上，這種情況發(fā)生概率最高，高郵站為4/6、通州站為4/5、宜興站為7/8；二是遷入峰次不多，但單次遷入量特大，也達到了全年發(fā)生程度4級以上，這種情況局地性很強，如高郵站的1995、2006年，通州站的2006年，宜興站的2006年，特別是2006年整個長江中下游稻區(qū)在8月底至9月初，因副高短期振蕩出現(xiàn)了兩次特大遷飛量和“飛虱雨”。當然也有遷入峰次很多，但單一峰次蟲量不是很大，總蟲量達不到4級以上的情況，如高郵站的1985、1990年。

圖5 高郵、通州和宜興3站1983—2010年褐飛虱發(fā)生等級與降水日數(shù)

圖6 高郵、通州、宜興三站1981—2010年降水距平百分率的年際變化圖

一個時間段的降水量累積效應可以從降水強度、降水持續(xù)性和降水頻次3個方面來體現(xiàn)，同樣的降水量值有可能由4種情況造成的：一是降水強度大，降水頻次多，但降水歷時短(如陣性降水，多由強對流性天氣導致)；二是降水頻次多，但降水強度小，降水歷時短(如間歇性降水)；三是降水歷時長，但降水強度不大，降水頻次也不多(如連續(xù)性降水)；四是降水強度特大，但降水歷時不長，降水頻次也少(如暴雨、大暴雨甚至特大暴雨)。第一、二、三種情況導致的降水持續(xù)效應，能使多日、多批次(可以是連續(xù)的，也可以是不連續(xù)的)遷飛途經(jīng)該降水區(qū)域的褐飛虱種群受到一定大小和力度的雨滴下沖及降水引起的拖曳下沉氣流的動力脅迫而降落地面，且不致于沖刷致死，能存活上燈，從而導致某一時間段或某一遷飛季的大遷入量和多遷入峰次。第四種情況對種群的遷入意義不大，即使有遷入，也會被沖刷至死，不能上燈計數(shù)；因此，降水對褐飛虱遷入的影響是有條件的，只有遷飛季內(nèi)降水在一定強度以下、頻次多、累積量大的年份才有利于褐飛虱偏重以上遷入的發(fā)生。

綜合降水量和降水頻次兩方面的影響可以看出：3站褐飛虱遷入偏重以上發(fā)生的年份中，高郵站6a中有4a 的6—10月遷飛季內(nèi)降水量大、降水頻次多，通州站5a中有4a降水量大、降水頻次多，宜興站8a中有7a降水量大、降水頻次多。由此可以看出；一定強度(中雨強度以下)以內(nèi)的降水，在蟲源輸送得到保證的前提下，降水區(qū)內(nèi)降水累積量越大、降水頻次越高，褐飛虱遷入累計量越大、遷入峰次也越多，最終導致其遷入和災變的發(fā)生程度越嚴重。

為了定量反映降水異常對褐飛虱遷入的影響，本文從氣候?qū)W角度將當月總降水量與當月降水日數(shù)的比值定義為月降水強度，分別以歷年6—10月降水量、降水日數(shù)和降水強度為預報因子，對褐飛虱遷入發(fā)生等級進行回歸分析。3個站點褐飛虱發(fā)生等級的回歸方程復相關系數(shù)(R)都較大，通過了F檢驗，置信度(P)達0.001。3個方程的回代檢驗正確率都在60%—80%之間，說明降水對遷入有明顯的影響。具體回歸方程為：

高郵站Y=1.577+0.006X1+0.205X2

式中，Y為褐飛虱遷入發(fā)生等級；X1為7月降水量，X2為6月降水強度；R=0.649；F=9.088；P=0.001

其中7月降水量對于褐飛虱的遷入影響比6月降水強度大。

通州站Y=0.121+0.622X1+0.005X2

Y為褐飛虱遷入發(fā)生等級；X1為8月降水強度，X2為7月降水量；R=0.639；F=8.622；P=0.001。

其中8月降水強度對褐飛虱遷入的影響比7月降水量大。

宜興站Y=-0.328+0.202X1+0.124X2+0.165X3+0.178X4

Y為褐飛虱遷入發(fā)生等級；X1為7月降水強度，X2為8月降水日數(shù)，X3為10月降水日數(shù)，X4為9月降水日數(shù)；R=0.742；F=7.040；P=0.001。

其中7月降水強度對褐飛虱遷入的影響最大，其它因子的作用大小依次為8月降水日數(shù)、10月降水日數(shù)和9月降水日數(shù)。

2.3.2 溫度對褐飛虱遷入的影響

通過計算高郵、通州和宜興三站的整個夏季(6—8月)的平均氣溫距平值和秋季褐飛虱為害期(9—10月)的平均氣溫距平值，得到圖7所示的變化曲線。分析其特征發(fā)現(xiàn)：這3個站點氣溫變化趨勢比較一致，遷飛季內(nèi)氣溫波動都較大?？紤]到1983—2010年的28a中，1994年是全球變暖突變年，1983—1993年是冷期，1994年以來是暖期，以1994年為分界年，分別計算冷期、暖期兩個時段的6—8月和9—10月氣溫平均值，再結合圖1和圖7分析可以發(fā)現(xiàn)：在高郵站褐飛虱偏重或大發(fā)生的6a中，夏溫低的年份為1983、1987、1991年，正常的為1995、2007年，偏高的為2006年；秋季(9—10月)氣溫偏高的年份有6a。通州站褐飛虱偏重或大發(fā)生的5a中有4a是夏溫低的(2006年夏溫偏高)；秋溫高的則有4a(1991年秋溫偏低)。宜興站偏重或大發(fā)生的8a中，夏溫低的年份有4a(1983、1985、1987和1991年)，另外4a(1988、1990、2005和2006年)夏溫偏高；而秋溫高的年份有6a(1988、1990秋溫偏低)。

雖然2006年夏秋季長江中下游地區(qū)受副高控制，氣溫偏高，但副高短期振蕩明顯，當年8月底至9月初，中國大陸20—30°N之間的區(qū)域均為一致的偏南氣流，這對褐飛虱的北遷極為有利，而33°N附近已經(jīng)存在較強的西風帶，并一直向南移動、擴展，與偏南氣流在江蘇匯合，導致江蘇地區(qū)出現(xiàn)了歷史罕見的集中遷入[28]。通州站1991年秋季雖然溫度較低，但雨日和降水量均多，所以褐飛虱遷入量也較多。1988年9月22日有一個較強的臺風8817號臺風(風力在12級以上)登陸廣東，使得浙江南部、江西東部、福建和廣東東部等地出現(xiàn)了大暴雨，而安徽南部、浙江北部和江蘇南部等地區(qū)也有一定強度(小到中雨級別)的降水，所以宜興出現(xiàn)了褐飛虱的集中降落，而通州和高郵位于江蘇東部和中部，臺風對這兩個地區(qū)沒有影響，三維氣流場也不利于褐飛虱遷入，因此，這兩地也就沒有褐飛虱的明顯遷入。1990年7月11—15日和9月10—15日，宜興一直為雨日，而通州和高郵降水量和雨日都很少，所以宜興在這幾日內(nèi)有大量褐飛虱遷入，而通州和高郵遷入量較少。由此可見：在蟲源條件滿足的前提下，夏秋季三維流場、降水和溫度場的配置對褐飛虱的遷入有十分重要的影響，動力條件是關鍵，熱力條件只起輔助作用。

圖7 高郵、通州和宜興三站1983—2010年夏季和秋季平均溫度距平值年際變化

因此，秋溫高對褐飛虱的遷入和為害更有利，因為褐飛虱生長發(fā)育和繁殖的適宜溫度為20—30℃，尤其以26—28℃最為適宜，溫度過高或過低都不利于其生長發(fā)育。而且，秋季溫度在一定范圍內(nèi)的偏高(日平均氣溫不超過28℃)有利于褐飛虱的滯留，發(fā)生代數(shù)會比常年增加。而夏季氣溫偏低，必須有降水日數(shù)偏多來作為大遷入總量、多遷入頻次的支撐條件，才能引起夏季褐飛虱偏重以上程度的遷入和田間為害。

由此可見，“盛夏不熱，晚秋不冷，夏秋多雨”的氣候確實有利于江蘇稻區(qū)褐飛虱的遷入與為害。

3 結論與討論

本文根據(jù)1983—2010年江蘇省褐飛虱燈誘蟲情資料、1981—2010年ONI指數(shù)和SOI指數(shù)及1983—2010年江蘇省氣溫和降水數(shù)據(jù)，分析了江蘇省褐飛虱的發(fā)生特征，探討了氣候異常對江蘇省褐飛虱災變性遷入的影響。研究得到以下結論：

(1)赤道東太平洋海表水溫持續(xù)偏高的厄爾尼諾事件易引發(fā)江蘇省褐飛虱偏重以上程度的發(fā)生，且其首次遷入峰的出現(xiàn)時間與厄爾尼諾事件開始期之間有1—14個月的滯后期。

(2)南方濤動指數(shù)(SOI)持續(xù)出現(xiàn)負值的厄爾尼諾事件發(fā)生的當年或次年，褐飛虱為偏重以上發(fā)生程度，且其首次遷入峰的出現(xiàn)時間也滯后于SOI負值開始期1—14個月。

(3)登陸并影響江蘇的強熱帶氣旋偏多的拉尼娜年份會出現(xiàn)褐飛虱偏重以上的發(fā)生。

(4)降水對褐飛虱遷入的影響明顯，但有限制條件；降水量大、降水日數(shù)多的年份褐飛虱遷入量大、發(fā)生程度重；以6—10月降水量、降水日數(shù)和降水強度為預報因子，建立褐飛虱發(fā)生等級的回歸方程，方程擬合和試報效果較好。

(5)夏秋季6—10月氣溫的異常對江蘇省褐飛虱遷入的發(fā)生有一定的影響，夏溫低、秋溫高的年份會導致江蘇省褐飛虱偏重以上的發(fā)生，但秋溫高的作用比夏溫低更明顯。

海溫場以熱力作用為源通過海氣相互作用影響大氣環(huán)流，繼而驅(qū)動氣候變化，控制或影響天氣系統(tǒng)、天氣過程，最終影響褐飛虱的遷入和災變。南方濤動以東太平洋副熱帶高壓與低緯印度洋和西太平洋赤道低壓之間“蹺蹺板”式的負相關來體現(xiàn)南半球中、低緯氣壓場低頻振蕩，從動力作用上實現(xiàn)海氣相互作用對大氣環(huán)流調(diào)整的影響。兩者共同影響的結果是造成全球行星風帶、大氣環(huán)流及季風系統(tǒng)、天氣系統(tǒng)和天氣過程的時空變化，繼而影響褐飛虱蟲源地、遷飛區(qū)、降落區(qū)的天氣、氣候變化，最終影響到褐飛虱在我國的區(qū)域性遷入和災變。本文研究發(fā)現(xiàn)的褐飛虱遷入的偏重以上發(fā)生與ENSO事件、氣候異常之間存在的關系只是初步的，至于深層次的定量關系研究及海溫場—海氣相互作用—大氣環(huán)流—天氣、氣候—褐飛虱遷入量傳遞影響機制的進一步探討是下一步的重要研究目標。

：

[1] Zhu S X. Rice Planthoppers and Their Prevention. Shanghai: Shanghai Science and Techenology Press, 1984: 1- 110.

[2] Cheng J A, Zhu Z R. Analysis on the key factors causing the outbreak of brown planthopper in Yangtze Area, China in 2005. Plant Protection, 2006, 32(4): 1- 4.

[3] Qi G J, Lu F, Hu G, Wang F Y, Cheng X N, Shen H M, Huang S S, Zhang X X, Zhai B P. Dynamics and population analysis of the brown planthopper Nilaparvata lugens (St?l) in the early rice field in Guangxi Municipality, 2007. Acta Ecologica Sinica, 2010, 30(2): 462- 472.

[4] Wu G X, Li J P, Zhou T J, Lu R Y, Yu Y Q, Zhu J, Mu M, Duan A M, Ren R C, Ding Y H, Li W J, He J H, Wang F, Yu W D, Qiao F L, Yuan D L, Qi Y Q. The key region affecting the short-term climate variations in China: the joining area of Asia and Indian-Pacific Ocean. Advances in Earth Science, 2006, 21(11): 1109- 1118.

[5] Xu Z Q, Fan K. Possible process for influences of winter and spring Indian Ocean SST anomalies interannual variability mode on summer rainfall over eastern China. Chinese Journal of Atmospheric Sciences, 2012, 36(5): 879- 888.

[6] Huang R H, Wu Y F. The influence of ENSO on the summer climate change in China and its mechanism. Advances in Atmospheric Sciences, 1989, 6(1): 21- 32.

[7] Liu Y Q, Ding Y H. Reappraisal of the influence of ENSO events on seasonal precipitation and temperature in China. Scientia Atmospherica Sinica, 1995, 19(2): 200- 208.

[8] Li S S, Shou S W. Equatorial Eastern Pacific SST and analysis on causes of summer floods/droughts in the Changjiang and Huaihe River Basin. Quarterly Journal of Applied Meteorology, 2000, 11(3): 331- 338.

[9] Sun S Q, Ma S J. Analysis and numerical experiment on the relationship between the 1998 summer monsoon activities and SSTA in tropical regions. Chinese Journal of Atmospheric Science, 2003, 27(1): 36- 51.

[10] Xu W C, Ma J S, Wang W. A review of studys on the influence of ENSO events on the climate in China. Journal of the Meteorological Sciences, 2005, 25(2): 212- 220

[11] Xin Z B, Xie Z R. The impact of ENSO events on Huaihe River Basin′s preciptation. Scientia Meteorologica Sinica, 2005, 25(4): 346- 354.

[12] Sun Z B, Chen H S, Tan G R. Basis on Short-Term Climate Prediction. Beijing: Meteorological Press, 2010: 126- 138.

[13] Jiang G H, Tan H Q, Shen W Z, Cheng X N, Chen R C. The relation between long-distance northward migration of the brown planthopper (NilaparvatalugensSt?l) and synoptic weather conditions. Acta Entomologica Sinica, 1981, 24(3): 251- 261.

[14] Drake V A, Farrow R A. The influence of atmospheric structure and motions on insect migration. Annual Review of Entomology, 1988, 33(1): 183- 210.

[15] Pu M H, Chen J M. A preliminary study on statitical forecast of brown planthopper ocurrence grade. Plant Protection, 1979, 5(5): 1- 9.

[16] Chen R C, Qi L Z, Cheng X N. Studies on the population dynamics of brown planthopper,Nilaparvatalugens(St?l), Ⅰ. effects of temperature and diet conditions on the growth of experimental population. Journal of Nanjing Agricultural University, 1986, 9(3): 23- 33.

[17] Zhao S J. Analysis on the meteorological cause of rice planthoppers outbreak in 1987. Plant Protection, 1988, 14(2): 2- 5.

[18] Bao Y X, Zhai B P, Cheng X N. Numerical simulation of the migration parameters of the brown planthopper, Nilaparvata lugens (St?l). Acta Ecologica Sinica, 2005, 25(5): 1107- 1114.

[19] Huo Z G, Li M S, Li N, Wang L, Huang D P, Wang C Y. Impacts of seasonal climate warming on crop diseases and pests in China. Scientia Agricultura Sinica, 2012, 45(11): 2168- 2179.

[20] He X C, Ding Y H, He J H. Response characteristics of the East Asian winter monsoon to ENSO events. Chinese Journal of Atmospheric Science, 2008, 32(2): 335- 344.

[21] Ropelewski C F, Halpert M S. Global and regional scale precipitation patterns associated with the El Nio/Southern Oscillation. Monthly Weather Review, 1987, 115(8): 1606- 1626.

[22] Crop Diseases and Pests Monitoring and Forecasting Station of Ministry of Agriculture in China. The Observation and Investigation Standard of Rice Planthoppers (GB/T 15794- 1995). Beijing: China Agriculture Press, 1995.

[23] Cheng X N, Chen R C, Xi X, Yang L M, Zhu Z L, Wu J C, Qian R G, Yang J S. Studies on the migrations of brown planthopper,Nilaparvatalugens(St?l). Acta Entomologica Sinica, 1979, 22(1): 1- 21.

[24] Huo Z G, Chen L, Ye C L, Liu L. Effect of climate on outbreak of China rice planthopper. Journal of Natural Disasters, 2002, 11(1): 97- 102.

[25] Li S H. The influence of Greenhouse effect on rice insect pests. Journal of Catastrophology, 1995, 10(2): 43- 47.

[26] Tang J Y, Hu B H, Wang J Q. Outbreak analysis of rice migratory pests in China and management strategies recommended. Acta Ecologica Sinica, 1996, 16(2): 167- 173.

[27] Lu X F, Huo Z G, Shen S H, Huang D P, Wang L, Xiao J J, Yu C X. Effects of climate warming on the northern distribution boundary of brown planthopper (Nilaparvatalugens(St?l)) overwintering in China. Chinese Journal of Ecology, 2012, 31(8): 1977- 1983.

[28] Guo Y J, Ni Y Q. Effects of the tropical Pacific convective activities on precipitation and temperature over China in summer. Meteorological Monthly, 1998, 24(6): 11- 16.

[29] Jin Z H, Tao S Y. A study on the relationships between ENSO cycle and rainfalls during summer and winter in Eastern China. Chinese Journal of Atmospheric Science, 1999, 23(6): 663- 672.

[30] Huang R H, Wu Y F. The influence of ENSO on the summer climat e change in China and its mechanism. Advances in Atmospheric Science, 1989, 6(1): 21- 30.

[31] Li W B, Zhou C P. The relationship between tropical Western Pacific warm pool and subtropical high. Acta Meteorologica Sinica, 1998, 56(5): 619- 626.

[32] Li C Y. A study on the influence of El Nino upon typhoon action over Western Pacific. Acta Meteorologica Sinica, 1987, 45(2): 229- 236.

[33] Yang S R, Zhan S Y, Zhao Q Y. A preliminary analysis of the relationship between the El Nino events and the typhoon in Western North Pacific. Marine Forecasts, 1987, 4(4): 10- 17.

[34] Yao D H, Chen X F, Yao S T, Chen Q H. Analysis of effects of typhoon Kanu(2005) on the sudden increase of Nilaparvata lugens. Chinese Journal of Agrometeorology, 2007, 28(3): 347- 349, 353- 353.

[35] Wang C H, Zhai B P, Bao Y X. Effects of typhoon "Haitang" airflow field on the northward migration route of rice brown planthopper. Chinese Journal of Applied Ecology, 2009, 20(10): 2506- 2512.

[36] Bao Y X, Li J J, Miao Q L, Shen S H, Wang J Q. Simulation of atmospheric dynamical background for a great migration event of the brown planthopper (NilaparvatalugensSt?l): a case study. Chinese Journal of Agrometeorology, 2008, 29(3): 347- 352.

參考文獻:

[1] 朱紹先. 稻飛虱及其防治. 上海: 上?？茖W技術出版社, 1984: 1- 110.

[2] 程家安, 祝增榮. 2005年長江流域稻區(qū)褐飛虱暴發(fā)成災原因分析. 植物保護, 2006, 32(4): 1- 4.

[3] 齊國君, 蘆芳, 胡高, 王鳳英, 程遐年, 沈慧梅, 黃所生, 張孝羲, 翟保平. 2007年廣西早稻田褐飛虱發(fā)生動態(tài)及蟲源分析. 生態(tài)學報, 2010, 30(2): 462- 472.

[4] 吳國雄, 李建平, 周天軍, 陸日宇, 俞永強, 朱江, 穆穆, 段安民, 任榮彩, 丁一匯, 李維京, 何金海, 王凡, 于衛(wèi)東, 喬方利, 袁東亮, 齊義泉. 影響我國短期氣候異常的關鍵區(qū): 亞印太交匯區(qū). 地球科學進展, 2006, 21(11): 1109- 1118.

[5] 徐志清, 范可. 冬季和春季印度洋海溫異常年際變率模態(tài)對中國東部夏季降水的可能影響過程. 大氣科學, 2012, 36(5): 879- 888.

[7] 劉永強, 丁一匯. ENSO事件對我國季節(jié)降水和溫度的影響. 大氣科學, 1995, 19(2): 200- 208.

[8] 勵申申, 壽紹文. 赤道東太平洋海溫與我國江淮流域夏季旱澇的成因分析. 應用氣象學報, 2000, 11(3): 331- 338.

[9] 孫淑清, 馬淑杰. 海溫異常對東亞夏季風及長江流域降水影響的分析及數(shù)值試驗. 大氣科學, 2003, 27(1): 36- 51.

[10] 許武成, 馬勁松, 王文. 關于ENSO事件及其對中國氣候影響研究的綜述. 氣象科學, 2005, 25(2): 212- 220.

[11] 信忠保, 謝志仁. ENSO事件對淮河流域降水的影響. 氣象科學, 2005, 25(4): 346- 354.

[12] 孫照渤, 陳海山, 譚桂容. 短期氣候預測基礎. 北京: 氣象出版社, 2010: 126 - 138.

[13] 江廣恒, 談涵秋, 沈婉貞, 程遐年, 陳若篪. 褐飛虱遠距離向北遷飛的氣象條件. 昆蟲學報, 1981, 24(3): 251- 261.

[15] 浦茂華, 陳潔明. 褐稻虱發(fā)生程度數(shù)理預報的初步研究. 植物保護, 1979, 5(5): 1- 9.

[16] 陳若篪, 綦立正, 程遐年. 褐飛虱種群動態(tài)的研究—Ⅰ. 溫度、食料條件對種群增長的影響. 南京農(nóng)業(yè)大學學報, 1986, 9(3): 23- 33.

[17] 趙圣菊. 從氣象因素分析1987年稻飛虱大發(fā)生的原因. 植物保護, 1988, 14(2): 2- 5.

[18] 包云軒, 翟保平, 程遐年. 褐飛虱遷飛參數(shù)的數(shù)值模擬. 生態(tài)學報, 2005, 25(5): 1107- 1114.

[19] 霍治國, 李茂松, 李娜, 王麗, 黃大鵬, 王春艷. 季節(jié)性變暖對中國農(nóng)作物病蟲害的影響. 中國農(nóng)業(yè)科學, 2012, 45(11): 2168- 2179.

[20] 何溪澄, 丁一匯, 何金海. 東亞冬季風對ENSO事件的響應特征. 大氣科學, 2008, 32(2): 335- 344.

[22] 農(nóng)業(yè)部農(nóng)作物病蟲測報站. 稻飛虱測報調(diào)查規(guī)范(GB/T 15794- 1995). 北京: 中國農(nóng)業(yè)出版社, 1995.

[23] 程遐年, 陳若篪, 習學, 楊聯(lián)民, 朱子龍, 吳進才, 錢仁貴, 楊金生. 稻褐飛虱遷飛規(guī)律的研究. 昆蟲學報, 1979, 22(1): 1- 21.

[24] 霍治國, 陳林, 葉彩鈴, 劉玲. 氣候條件對中國水稻稻飛虱為害規(guī)律的影響. 自然災害學報, 2002, 11(1): 97- 102.

[25] 李淑華. 氣候變化對水稻蟲害的影響. 災害學, 1995, 10(2): 43- 47.

[26] 湯金儀, 胡伯海, 王建強. 我國水稻遷飛性害蟲猖獗成因及其治理對策建議. 生態(tài)學報, 1996, 16(2): 167- 173.

[27] 盧小鳳, 霍治國, 申雙和, 黃大鵬, 王麗, 肖晶晶, 于彩霞. 氣候變暖對中國褐飛虱越冬北界的影響. 生態(tài)學雜志, 2012, 31(8): 1977- 1983.

[28] 郭艷君, 倪允琪. 赤道太平洋對流活動異常對我國夏季降水和氣溫的影響. 氣象, 1998, 24(6): 11- 16.

[29] 金祖輝, 陶詩言. ENSO循環(huán)與中國東部地區(qū)夏季和冬季降水關系的研究. 大氣科學, 1999, 23(6): 663- 672.

[31] 李萬彪, 周春平. 熱帶西太平洋暖池和副熱帶高壓之間的關系. 氣象學報, 1998, 56(5): 619- 626.

[32] 李崇銀. 厄爾尼諾影響西太平洋臺風活動的研究. 氣象學報, 1987, 45(2): 229- 236.

[33] 楊淑瑞, 戰(zhàn)淑蕓, 趙黔源. 埃爾-尼諾與西北太平洋臺風關系的初步分析. 海洋預報, 1987, 4(4): 10- 17.

[34] 姚德宏, 陳雄飛, 姚士桐, 陳勤海. “卡努”(0515)臺風在褐飛虱突增中的作用探析. 中國農(nóng)業(yè)氣象, 2007, 28(3): 347- 349, 353- 353.

[35] 王翠花, 翟保平, 包云軒. “海棠”臺風氣流場對褐飛虱北遷路徑的影響. 應用生態(tài)學報, 2009, 20(10): 2506- 2512.

[36] 包云軒, 李金建, 繆啟龍, 申雙和, 王建強. 2006年褐飛虱重大遷入過程的大氣動力背景模擬——個例研究. 中國農(nóng)業(yè)氣象, 2008, 29(3): 347- 352.