刁家敏，孫思思，包云軒**，陸明紅，謝曉金

?

準雙周大氣低頻振蕩對江嶺稻區(qū)褐飛虱遷入的影響*

刁家敏1,2，孫思思1,2，包云軒1,2**，陸明紅3，謝曉金1,2

（1.南京信息工程大學氣象災害預報和評估協(xié)同創(chuàng)新中心，南京 210044；2.江蘇省農(nóng)業(yè)氣象重點實驗室/南京信息工程大學，南京 210044；3.農(nóng)業(yè)部全國農(nóng)業(yè)技術推廣與服務中心，北京 100125）

探明準雙周大氣低頻振蕩對中國褐飛虱遷入的影響，可為褐飛虱的延伸期預報提供重要的參考依據(jù)。根據(jù)1979-2015年中國105個植物保護站的逐日褐飛虱燈誘資料、同期NCEP/NCAR全球逐日氣象再分析資料，篩選出具有典型褐飛虱偏重以上發(fā)生特征的年份（2010年），分析江嶺稻區(qū)褐飛虱遷入量與準雙周低頻大氣振蕩之間的關系；利用HYSPLIT軌跡計算模型對遷入典型站點的褐飛虱種群作后向軌跡分析，探討大氣低頻流場對褐飛虱遷飛路徑和蟲源地的影響。結(jié)果表明：（1）2010年8-9月江嶺稻區(qū)褐飛虱遷入量存在顯著的10～20d振蕩周期；（2）同期850hPa和925hPa兩個高度的風場、位勢高度場及垂直速度場均存在10～20d的顯著低頻周期，且褐飛虱遷入量的低頻振蕩與這些大氣物理量的低頻振蕩呈顯著相關；（3）前期江嶺稻區(qū)主要受低頻西南風的影響，褐飛虱種群隨西南風由西南稻區(qū)向江嶺稻區(qū)輸送，后期部分褐飛虱種群隨東北風由江淮稻區(qū)向江嶺稻區(qū)輸送；（4）褐飛虱種群的后向軌跡模擬計算結(jié)果，進一步驗證了遷入江嶺稻區(qū)的褐飛虱的軌跡和蟲源地與低頻風場的盛行基本一致。

褐飛虱；小波分析；低頻流場；遷飛軌跡；蟲源地

褐飛虱[（St?l）]是水稻主要害蟲之一，其體積小、自主飛行速度慢，是風載遷飛性昆蟲，可隨風作遠距離遷飛，故有季節(jié)性、遠距離遷飛性等特點[1-2]。褐飛虱病蟲害在中國發(fā)生頻繁且危害嚴重，每年受害水稻面積有6×106～10×106hm2，約占總種植面積的50%或以上，年均稻谷損失達10×108～15×108kg，大發(fā)生年份造成的損失更大[3]。江嶺稻區(qū)為單雙季稻種植混合區(qū)，常年褐飛虱發(fā)生程度中等偏重，該地區(qū)褐飛虱在8月中下旬向江淮稻區(qū)輸送，9月上旬開始向南回遷，9月中旬出現(xiàn)由江淮稻區(qū)向江嶺稻區(qū)的回遷峰，9月下旬-10月上旬由長江中下游遷往南嶺稻區(qū)[4]。

褐飛虱災變性暴發(fā)需要在生理條件滿足的前提下，有較好的寄主環(huán)境和氣象條件[5]。在影響褐飛虱遷飛和災變性遷入的大氣背景（或氣象條件）中，不同時間尺度的大氣運動系統(tǒng)及其變化所產(chǎn)生的物理效應（包括動力效應、熱力效應、水分脅迫效應等）是關鍵因素。這些大氣運動系統(tǒng)及其變化按時間尺度由短至長依次可分為高頻變化、天氣變化、低頻變化、季節(jié)變化、甚低頻變化、年代際變化和地質(zhì)紀變化[6]。其中時間尺度在10d以上、一個季節(jié)以內(nèi)的大氣運動變化稱為大氣低頻變化，它包括了季節(jié)內(nèi)振蕩（30～60d）和準雙周振蕩（10～30d）[7-9]。

以往研究表明，包括氣壓場、風場、垂直速度場、溫度場、降水場等在內(nèi)的許多大氣物理量場[10-14]均存在顯著的大氣低頻振蕩，季風、副高、臺風等大氣環(huán)流系統(tǒng)均明顯受大氣低頻振蕩的影響[15-17]。包云軒等研究表明，褐飛虱遷飛受到了氣壓場（或位勢高度場）及其產(chǎn)生的大氣環(huán)流系統(tǒng)（如季風、副高、臺風等）的影響[18-21]；其遠距離遷飛的方向與高空風向基本一致[22]；強下沉氣流對降蟲的作用十分明顯，可使褐飛虱種群大量降落，降水有利于褐飛虱的降落[3,20]；褐飛虱遷飛的高度隨著氣溫的升降而升降，二者具有相對一致性[18,22]。故研究氣壓場、風場、垂直速度場等大氣物理量的低頻變化與褐飛虱遷入量低頻變化的關系可進一步探明大氣背景對褐飛虱遷入的影響，為褐飛虱遷入量的預報和遷入峰的預警提供重要的參考依據(jù)。

國內(nèi)外關于短期天氣過程和氣候變化對害蟲遷入影響的研究較多，但有關大氣低頻振蕩對害蟲遷飛影響的研究尚未見報道。準雙周大氣低頻振蕩（Quasi-biweekly Atmospheric Low Frequency Oscillation，QBWO）是以10～20d周期為主的振蕩[21]，參照前人的研究方法，本文對影響褐飛虱遷飛的部分大氣物理量作小波分析，根據(jù)其10～20d的振蕩周期規(guī)律進行低頻濾波處理，江嶺稻區(qū)褐飛虱遷入量有明顯10～30d的周期振蕩規(guī)律，故重點研究江嶺稻區(qū)典型遷飛季褐飛虱的10～20d的低頻遷入特征，進一步探討其與低頻大氣物理量場的關系，并通過HYSPLIT模式對褐飛虱的遷飛軌跡及蟲源地進行模擬驗證，旨在為褐飛虱的延伸期預報提供理論基礎。

1 資料與方法

研究資料包括：（1）由農(nóng)業(yè)部全國農(nóng)業(yè)技術推廣服務中心提供，包括1979-2015年中國105個植保站的逐日褐飛虱燈誘資料；（2）由美國國家環(huán)境預測中心（NCEP）和國家大氣研究中心（NCAR）提供的全球逐日氣象再分析資料，分辨率為2.5°×2.5°；（3）NOAA同期驅(qū)動軌跡模式HYSPLIT的資料（ftp://arlftp.arlhq.noaa.gov/pub/archives/reanalysis/）。

根據(jù)農(nóng)業(yè)部全國農(nóng)業(yè)技術推廣服務中心提供的褐飛虱發(fā)生程度分級數(shù)據(jù)，2010年全國褐飛虱發(fā)生程度為5級，取24°N-30°N、109°E-121°E范圍覆蓋中國淮河以南的南方水稻主產(chǎn)區(qū)（特別是湖南、江西和福建三省稻區(qū)）作為研究區(qū)域，考慮資料的完整性和代表性，選取具體監(jiān)測站點19個（其中湖南12個，江西2個，福建5個），其分布如圖1所示。

在Fortran中采用Morlet小波分析法[23]分析褐飛虱遷入量的低頻周期規(guī)律，利用Butterworth帶通濾波[24]（即通過濾去其它頻率僅保留某一頻率范圍的方法）和合成分析[25]等氣象統(tǒng)計方法，研究8-9月江嶺稻區(qū)褐飛虱遷入量及大氣物理量的主要低頻振蕩特點，再運用位相合成法[25]在GrADS中分析低頻褐飛虱遷入蟲量與各低頻大氣物理量的關系，最后利用HYSPLIT4.9采用后向軌跡模式模擬遷出蟲源地及遷飛軌跡。

HYSPLIT（Hybrid Single Particle Lagrangian Integrated Trajectory Model）模式是由美國國家海洋和大氣管理局（NOAA）的空氣資源實驗室（ARL）與澳大利亞國家氣象局（BOM）聯(lián)合開發(fā)，是供質(zhì)點軌跡、擴散及沉降分析用的綜合模式系統(tǒng)[26]。該模式通常用來跟蹤氣流所攜帶的粒子或氣體的移動方向，可以實時研究質(zhì)點的運動軌跡，分為前向軌 跡模式和后向軌跡模式，分別研究氣團的去向和來向，對無數(shù)個運動軌跡可以進行聚類分型。本文利用其聚類分析某一時間段內(nèi)褐飛虱種群的遷飛軌跡，其結(jié)果既能定性分析褐飛虱遷出蟲源，又能定量分析其空中遷飛軌跡及各路徑的貢獻大小[27]。

2 結(jié)果與分析

2.1 2010年8-9月江嶺稻區(qū)典型站點褐飛虱遷入量的低頻特征

將研究區(qū)觀測站點的逐日褐飛虱燈誘蟲量數(shù)據(jù)作小波分析，部分振蕩周期代表性站點的小波分析結(jié)果如圖2所示。由圖可見，所選站點（福建霞浦、湖南洪江、湖南雙峰、江西婺源）的褐飛虱燈誘蟲量在8-9月具有顯著的10～30d的振蕩周期，福建霞浦和江西婺源站的10～20d蟲量對數(shù)的濾波振幅強度大于20～30d，20～30d的蟲量對數(shù)的濾波與蟲量對數(shù)趨勢較為一致，但10～20d的趨勢性更好。

注：灰色區(qū)域通過0.10水平的顯著性檢驗。下同

Note: The gray area is P<0.10. The same as below

2.2 江嶺稻區(qū)研究時段內(nèi)逐日褐飛虱遷入量的準雙周低頻振蕩特征

將江嶺稻區(qū)19個站點的逐日褐飛虱遷入量進行區(qū)域平均，得到區(qū)域平均逐日褐飛虱遷入量序列，其小波分析結(jié)果見圖3。由圖3（a）可見，江嶺稻區(qū)褐飛虱遷入量低頻振蕩的主要周期為10～20d和20～30d，其中10～20d較顯著，且褐飛虱大發(fā)生過程的蟲量的波動變化與10～20d低頻分量相對應。由圖3（b）可見，當兩個低頻分量正（負）位相疊加，褐飛虱遷入量的低頻分量增強（中斷）；若反位相疊加，則褐飛虱遷入量的低頻分量減弱。進一步 對研究區(qū)內(nèi)每個站點的逐日褐飛虱遷入量序列進行小波變換分析（圖略）可見，其以10～20d周期的低頻發(fā)生期主要出現(xiàn)在9月，均具有準雙周振蕩周期。其統(tǒng)計結(jié)果見表1，由表可見，10～30d的振蕩主要分為10～20d周期和20～30d周期，19個站點中既有10～20d周期的褐飛虱遷入量振蕩周期，又有20～30d的振蕩周期。其中18個站點的褐飛虱遷入量有10～20d的振蕩周期，13個站點有20～30d的振蕩周期，說明江嶺稻區(qū)的褐飛虱振蕩周期主要為10～20d。

2.3 褐飛虱低頻遷入量與低頻大氣物理量場的相關性分析

對8-9月逐日850hPa和925hPa大氣物理量序列的小波分析（圖略）表明，其均存在顯著的10～20d低頻變化周期，且在9月與對應的褐飛虱遷入量的低頻振蕩周期基本對應，因此，進一步對9月的位勢高度（hgt）、合成風（wnd）和垂直速度（omega）分別進行10～20d低頻濾波分析，并與相應時段10～20d褐飛虱遷入量的低頻濾波序列進行對比，結(jié)果見圖4。由圖可見，850hPa和925hPa的10～20d濾波曲線的峰、谷與10～20d低頻遷入量的峰（谷）、谷（峰）能較好地對應，低頻蟲量與同頻域的低頻850hPa和925hPa的大氣物理量場之間存在較好的正（負）相關關系，其相關系數(shù)見表2。

圖3 江嶺稻區(qū)區(qū)域平均褐飛虱逐日燈誘蟲量序列的小波分析（a）、逐日蟲量對數(shù)化序列及其10～20d和20～30d的濾波曲線（b）

Fig. 3 Wavelet analysis ofarea averaged BPH’s daily light trapped catches sequence in the region between the Nanling Mountain and the Yangtze River（a）, the logarithm sequence of BPH’s daily light trapped catches and its 10-20d and 20-30d filter curves（b）

表 1 江嶺稻區(qū)19個植保站褐飛虱逐日燈誘蟲量的低頻振蕩周期

從圖4和表2還可看出，（1）褐飛虱的低頻遷入量與前一位相850hPa和925hPa平均位勢高度場成顯著負相關，表明前一位相（前48h內(nèi)）蟲源區(qū)的負變高區(qū)（即上升氣流區(qū)）對后一位相（后48h）褐飛虱降蟲區(qū)的遷入量有正向促進作用，反之則呈反向削弱作用；（2）褐飛虱的低頻遷入量與前一位相上游區(qū)域850hPa和925hPa平均合成風場呈顯著正相關，表明前48h內(nèi)偏南風（合成風速為正值）越強，越有利于后48h內(nèi)褐飛虱的遷入；（3）褐飛虱的低頻遷入量與同一位相內(nèi)本地垂直速度場呈一定程度的正相關，表明48h內(nèi)本區(qū)域上空的下沉氣流對蟲源遷入有利，其中925hPa的強下沉氣流對降蟲十分有利。

注：虛線表示褐飛虱遷入量濾波序列的±0.4個標準差

Note: Dashed lines represent ±0.4SD(standard deviation) of the 10-20d filter frequency immigration sequence of BPH

表2 低頻位勢高度（hgt）、合成風（wnd）、垂直速度（omega）與低頻褐飛虱遷入量的相關系數(shù)

注：**表示P＜0.01，*表示P＜0.05。

Note:**is P<0.01,*is P<0.05.

2.4 低頻大氣物理量場對江嶺稻區(qū)褐飛虱遷入量的影響分析

2.4.1 低頻振蕩周期的位相設置

將一個振蕩周期分為8個位相，由圖5可見，位相3和位相7對應波列的波峰和波谷，分別稱為波動的活躍和中斷位相，位相1（位相5）表示波動由中斷位相（活躍位相）轉(zhuǎn)向活躍位相（中斷位相），位相2、4、6、8對應振幅達極值位相振幅的一半。一個波動周期循環(huán)結(jié)束則回到下一個周期的1號位相。

2.4.2 低頻位勢高度場影響的逐位相分析

圖6中，從未濾波的925hPa平均高度場上（等值線）可以看出，30°E以東180°E以西的歐亞地區(qū)位勢高度場一直呈現(xiàn)出兩槽一脊配置，槽脊始終徘徊于江嶺稻區(qū)，上升氣流和下沉氣流交替影響著褐飛虱在該地區(qū)的遷入。

圖6中，在低頻褐飛虱遷入過程開始前（位相1），925hPa低頻高度場上烏拉爾山、鄂霍次克海、西伯利亞以東和中國華北、東北、南海地區(qū)均表現(xiàn)為位勢高度正異常（即正變高），而伊朗高原至印度半島、江嶺地區(qū)、東南沿海地區(qū)均表現(xiàn)為負異常。位相2，江嶺地區(qū)的負異常開始減弱，同時西南地區(qū)的正異常也減弱。到了極端活躍位相（即位相3），伊朗高原至印度半島的負異常區(qū)分裂成兩個負異常區(qū)，日本海處的正異常與菲律賓群島的正異常匯合成一個大的正異常，在位相1-位相2中處于負異常的江嶺、東南沿海等地區(qū)此時表現(xiàn)為正異常，對應有下沉氣流，有利于褐飛虱的降落。位相4，江嶺地區(qū)負異常加強，正異常減弱，意味著上升氣流在增強，而下沉氣流在減弱，故褐飛虱的遷入量開始減少。位相5-位相8低頻位勢高度場分布與位相1-位相4基本相反，為低頻褐飛虱遷入量減少階段。特別是在中斷位相，中國大部分地區(qū)均表現(xiàn)為負異常（除中國西北地區(qū)），對應有上升氣流，所以很不利于褐飛虱的降落，此階段褐飛虱的遷入量達到一個低頻振蕩周期的最低值。位相1-位相3，江嶺稻區(qū)從負異常向正異常轉(zhuǎn)變，上升氣流轉(zhuǎn)變?yōu)橄鲁翚饬?，低頻遷入蟲量增加，在位相4正異常減弱，低頻遷入蟲量開始減少；位相5-位相8的變化則與之相反。

2.4.3 低頻高空風輸送影響的逐位相分析

圖7給出了8個位相925hPa 的低頻合成風場，從中可以看出大氣環(huán)流對褐飛虱種群遷飛的影響。褐飛虱低頻遷入過程開始前（位相1），東北半球30°E以東180°E以西的歐亞地區(qū)，偏西方向上，阿拉伯海的南部和伊朗西南沿海各有一西風帶氣旋，青藏高原西南側(cè)有一反氣旋，其北側(cè)為偏西氣流；東南方向上，有一較強的反氣旋盤踞于菲律賓群島，其西側(cè)有來自中國南海的較強偏南氣流向中國大陸輸送；偏東方向上，在西北太平洋上雖有瀕臨消失的熱帶氣旋存在，但離中國大陸較遠，故對褐飛虱種群遷飛影響程度較??；偏北方向上，有一冷高壓位于蒙古高原南部，其東側(cè)的東北和偏北氣流對遠在江淮和長江中下游稻區(qū)的褐飛虱種群遷出有微弱的影響，使這些稻區(qū)有少量的遷出褐飛虱種群得以向西南偏南方向遷飛。

注：A表示反氣旋，C表示氣旋

Note: A is anticyclone, C is cyclone

至位相2，偏西方向上，原來位于伊朗西南沿海的氣旋東移至印度北部，位于阿拉伯海南部的氣旋東移至新馬群島北部并增強，位于青藏高原西南側(cè)的反氣旋卻減弱至無，而此時阿拉伯海有一氣旋新生。東南方向上，原來位于菲律賓群島的反氣旋西移北抬增強，其西南側(cè)的東南氣流與馬六甲海峽東側(cè)的西南氣流匯合增強，并向中國華南沿海和西南地區(qū)輸送，受其影響，華南大部、江嶺稻區(qū)南部盛行西南風，江嶺稻區(qū)北部盛行西北偏西風，這對華南稻區(qū)蟲源的北上和西南北部稻區(qū)蟲源的東輸有利。偏北方向上，位相1中蒙古高原南部的反氣旋減弱消失，華北地區(qū)開始盛行偏東風，江淮和長江中下游稻區(qū)偏北風分量減弱，導致褐飛虱種群的南遷量減少。因此，這一階段遷入江嶺稻區(qū)的褐飛虱種群仍主要來自于華南和西南稻區(qū)，少量來自于江淮稻區(qū)。

在極端活躍位相（位相3），偏西方向的印度氣旋繼續(xù)發(fā)展并向東北方向移動至阿富汗，逼近青藏高原，馬六甲海峽氣旋發(fā)展東移北抬，東南方向原位于菲律賓群島的反氣旋碎4減弱東移，但余威不減，受三者影響，中國華南、西南、江嶺南部盛行偏南風，江嶺北部盛行西風，有大量褐飛虱種群分別從華南稻區(qū)、西南稻區(qū)遷入江嶺稻區(qū)。

位相4，偏西方向上，原來位相3中位于青藏高原以西的氣旋跨越高原進入中國華中地區(qū)，使江嶺北部西南風分量增強；位于阿拉伯海的氣旋東移至印度半島，導致原來受偏南風控制的華南稻區(qū)、西南稻區(qū)和江嶺稻區(qū)南部的南風分量增強，同時在馬六甲海峽的氣旋東移至中國南海南部海面，使中國東南沿海的東南風分量增強，這三股褐飛虱運載氣流的增強，致使這一階段蟲源從華南稻區(qū)、西南稻區(qū)和東南沿海稻區(qū)源源不斷地遷入江嶺稻區(qū)。

位相5-位相8，受冷高壓東移南下影響，東南風逐漸減弱南退，西北風逐漸增強南壓，進入江嶺稻區(qū)的北遷褐飛虱蟲量逐漸減少，南遷蟲量逐漸增多。尤以低頻褐飛虱蟲量中斷位相（位相7），在偏西方向，在中國西南地區(qū)有一明顯的氣旋存在，西南稻區(qū)和江嶺稻區(qū)均受氣旋東北側(cè)東風的影響；在西北方向，中國西北部有一較強的冷式反氣旋南壓，其東側(cè)偏北風強勁南壓，其影響區(qū)域覆蓋江淮稻區(qū)大部和江嶺稻區(qū)北部，為以南遷為主的褐飛虱種群

進入江嶺稻區(qū)奠定了基礎；而江嶺稻區(qū)東南部仍受東南偏東風影響，部分蟲源從東南沿海向西北內(nèi)陸輸送。由于水稻生育期、地面溫場脅迫等迫使褐飛虱種群大規(guī)模南遷的因素還不成熟，故這一階段南遷降落江嶺稻區(qū)的蟲量并不多，這也是位相7開始，低頻遷入蟲量明顯減少的主要原因。

由位相1-位相8可知，此次褐飛虱低頻遷入過程，前期在位相1-位相4期間，江嶺稻區(qū)褐飛虱種群遷入受偏南盛行風影響，以北遷為主，東西向遷入為輔，偶有少量南遷；后期在位相5-位相8期間，受隨盛行風左轉(zhuǎn)影響，由東南向西北的北遷轉(zhuǎn)向由東北向西南的南遷，最后完全由南遷控制。另外，從上述分析還可以發(fā)現(xiàn)：（1）8月底-9月下旬，適宜遷飛的925hPa低頻風場受東亞季風轉(zhuǎn)換的影響存在一個明顯的盛行風向逆時針轉(zhuǎn)向的過程。（2）受其影響，這一盛行風向轉(zhuǎn)向期的褐飛虱種群遷飛不是單一的北遷或南遷盛行期，而是從北遷為主、南遷為輔的混合遷轉(zhuǎn)向以南遷為主、北遷為輔的混合遷，最終轉(zhuǎn)入南遷盛行期。

2.4.4 低頻垂直速度場影響的逐位相分析

由表2可知，在925hPa高度上褐飛虱遷入蟲量與垂直速度的相關系數(shù)大于0.5，且通過0.01水平的顯著性檢驗，說明二者呈極顯著的正相關關系。而在850hPa高度上，低頻垂直速度場與褐飛虱遷入量 的低頻分量相關性相對較小，加之，8月底以后，褐飛虱的適宜遷飛高度已開始明顯下降至1000m左右，更接近925hPa（750m高度）[28]，故本文不考慮850hPa低頻垂直速度場。

在圖8中，位相1期間，華南和西南地區(qū)有強上升氣流，江嶺和江淮地區(qū)有強下沉氣流，有利于褐飛虱種群從華南稻區(qū)和西南稻區(qū)遷出，在江嶺稻區(qū)和江淮稻區(qū)降落。位相2期間，江嶺稻區(qū)仍維持強下沉氣流，江淮稻區(qū)的強下沉氣流區(qū)消失，華南和西南稻區(qū)的強上升氣流區(qū)東移減弱，垂直氣流對江嶺稻區(qū)降蟲作用仍十分明顯。位相3期間，江淮稻區(qū)強上升氣流區(qū)出現(xiàn)，對蟲源起飛南遷有利，江嶺稻區(qū)強下沉氣流區(qū)擴展并增強，有利于降蟲。位相4期間，江淮稻區(qū)強上升氣流區(qū)進一步擴展增強，垂直氣流場對北遷降蟲作用更弱，對南遷降蟲增強作用更明顯。位相5-位相8期間，江嶺稻區(qū)下沉氣流區(qū)強度逐漸減弱、范圍逐漸減小，導致降蟲量和降蟲范圍均減少，最終趨于消失（圖略）。

注：-為上升氣流、+為下沉氣流

Note：-is updraft airflow, + is downdraft airflow

2.5 研究期湖南雙峰站褐飛虱遷移規(guī)律驗證

選取湖南雙峰站（27.44°N，112.18°E），利用HYSPLIT后向軌跡模式進行褐飛虱后向軌跡模擬，以研究江嶺稻區(qū)8-9月的褐飛虱遷飛軌跡，驗證低頻風場褐飛虱遷入的規(guī)律。根據(jù)褐飛虱生物特性，首先假設：（1）褐飛虱順風遷移；（2）褐飛虱在日出前1h或日落后1h內(nèi)起飛。模擬時段為此次典型褐飛虱遷入蟲量的低頻振蕩的前半個周期（位相1-位相4），選取時間步長為12h，取通用協(xié)調(diào)時UTC時間為22：00（即BJT6：00）和10：00（即BJT18：00）為后推時間起點；選取1000m（925hPa）作為模擬高度；與上述周期分析對應，模擬空氣塊在位相2、3、4后向追蹤各48h的三維運動軌跡。然后利用HYSPLIT的聚類分析功能對褐飛虱的遷飛軌跡進行聚類，得到主要的遷飛軌跡及相應遷飛路徑的褐飛虱蟲量貢獻率（圖9）。

如圖9所示，位相1-位相2（圖9a），8-9月925hPa褐飛虱遷飛軌跡有5條，主要為4條偏南路徑和1條東北路徑。遷飛軌跡1和2 均起始于廣西稻區(qū)，一條起始于廣西中部地區(qū)，一條起始于廣西東南部地區(qū)，貢獻率分別為38%和13%，兩者貢獻率之和大于50%；遷飛軌跡3、4分別起始于廣東省的西南部和北部地區(qū)，貢獻率分別為17%和25%；最后一條軌跡5的起始點位于江西西北部地區(qū)，貢獻率為8%。由此可知，位相1-位相2的48h遷飛過程中，湖南雙峰站遷入的褐飛虱蟲源地主要是華南稻區(qū)，小部分褐飛虱種群從江西向西南遷入湖南雙峰站。

位相2-位相3（圖9b），8-9月925hPa褐飛虱遷飛軌跡有西南和東北兩條路徑。向后回推2d，西南路徑起始于廣西東北部，貢獻率達69%；東北路徑起始于安徽西南部，褐飛虱沿江西與湖北交界地區(qū)向湖南地區(qū)遷入，貢獻率達31%。由此可知，在此次遷飛過程中，仍以北遷的褐飛虱種群居多，南遷為少數(shù)。

注：圖中括號內(nèi)的百分數(shù)為該軌跡褐飛虱蟲量的貢獻率

Note: The percentage is the contribution rate of different trajectories in the total catches of BPH

位相3-位相4（圖9c），8-9月925hPa褐飛虱遷飛軌跡有兩條：一條主徑為偏南路徑，另一條次徑為偏東路徑。遷飛軌跡1起始點在福建西部地區(qū)，經(jīng)過江西北部向西遷飛。遷飛軌跡2起始于海南地區(qū)，途經(jīng)廣東東部、福建西部、江西南部，向湖南遷飛。遷飛軌跡1的貢獻率為69%，遷飛軌跡2的貢獻率為31%。由此可知，位相1-位相2的48h遷飛過程中，湖南雙峰站遷入的蟲源地主要是西南稻區(qū)，后隨風遷移至湖南地區(qū)。

將通過低頻位勢高度場影響的分析結(jié)果，與本節(jié)對湖南雙峰站水汽輸送路徑相比，HYSPLIT4.9后向軌跡模式基本驗證了圖7所得的結(jié)果，即位相1中江嶺稻區(qū)有從西南地區(qū)輸送來的暖濕氣流，也有較弱的東北風，至位相2江嶺稻區(qū)除了西南地區(qū)的暖濕氣流，華南地區(qū)也有氣流向北輸送，褐飛虱種群仍主要來自于華南和西南稻區(qū)，少量來自江淮稻區(qū)，與圖9a中位相1-位相2的后向軌跡相吻合。而位相3江嶺稻區(qū)盛行偏南風，至位相4改為盛行東南風，低頻風場的結(jié)果與HYSPLIT結(jié)果亦互相吻合，從而驗證了低頻位勢高度場分析的褐飛虱的遷飛路徑和蟲源地。由此可見，此次江嶺稻區(qū)的褐飛虱蟲源地主要是華南稻區(qū)，小部分種群從東北方向的沿江稻區(qū)向西南遷入江嶺稻區(qū)。

3 結(jié)論與討論

（1）2010年8-9月江嶺稻區(qū)19個站點褐飛虱遷入量的Morlet小波分析表明，江嶺稻區(qū)褐飛虱低頻振蕩的主要周期為10～20d和20～30d，其中10～20d振蕩更明顯；同期區(qū)域平均褐飛虱遷入量也存在10～20d和20～30d的振蕩周期，且10～20d的振蕩更明顯。

（2）同期850hPa和925hPa兩個高度的位勢高度場、風場、垂直速度場均存在10～20d的顯著低頻周期，且這些大氣物理量的低頻振蕩與褐飛虱遷入量的低頻振蕩有著較顯著的相關性。由于褐飛虱是一種典型的風載遷飛性昆蟲，低頻風場的演變主要通過改變盛行風場的結(jié)構(gòu)從而決定褐飛虱種群在空中的遷飛路徑。垂直速度場通過上升氣流和下沉氣流的時空分布影響褐飛虱種群的遷出蟲源地和降蟲區(qū)。而位勢高度場的四維起伏則是生成不同時段、不同氣壓場分布型的根本原因，在水平方向上它決定了高低壓區(qū)和槽脊的分布，也決定了渦旋區(qū)、切變區(qū)和輻散輻合區(qū)的分布，在垂直方向上它則決定了上升氣流和下沉氣流的強弱和分布。因此，三者的低頻振蕩實際上從動力角度對褐飛虱遷入量的低頻振蕩起了至關重要的作用。

（3）從低頻盛行風場對褐飛虱遷飛路徑的影響看，前期江嶺稻區(qū)褐飛虱的遷入主要受盛行的低頻偏南風影響，褐飛虱種群以北遷為主；后期受盛行的低頻東北風影響，以南遷為主；在低頻垂直速度場上，垂直氣流對褐飛虱的起飛和降落有一定影響，上升氣流有利于褐飛虱起飛，不利于其降落；下沉氣流有利于褐飛虱的降落，不利于其起飛。這均與以往研究結(jié)論一致[18-20]。由于褐飛虱遷入量的低頻振蕩與降水量的低頻振蕩相關性不顯著，故未考慮降水對降蟲的低頻動力作用，實際上，降水強度、降水歷時、降水分布區(qū)等因素的復雜性有可能使這種相關性變?nèi)?，而且降水對降蟲的作用一定程度上與拖曳下沉氣流的強弱有關，有可能降水與垂直氣流的交互作用掩蓋了降水的作用信號。根據(jù)對2010年歷史資料的查驗，此低頻遷入時段內(nèi)江嶺稻區(qū)有較大的降水過程發(fā)生[29]，且降水受到前期大氣低頻振蕩的影響[30]，因此，這一過程中，降水對褐飛虱的低頻遷入具有一定的促進作用。

（4）對褐飛虱遷飛進行軌跡模擬時，因未考慮其生物學特性和生態(tài)行為參數(shù)，導致模擬精度不高。由基于MM5模式驅(qū)動軌跡計算模型到應用MM5模式和GEARN模式結(jié)合的稻飛虱實時預報系統(tǒng)[31]，再到利用HYSPLIT模式和FLEXPART模擬褐飛虱的飛行軌跡[32-33]，準確性和精度都大大提升。HYSPLIT軌跡計算模式對遷入典型站點的褐飛虱種群作后向軌跡分析，發(fā)現(xiàn)遷入江嶺稻區(qū)的褐飛虱的軌跡和蟲源地與低頻風場的盛行基本一致。

以往對褐飛虱遷飛的研究，多為短期預報或從氣候的角度進行討論，少有從延伸期預報的角度考慮。通過相關分析篩選出對褐飛虱種群和低頻遷入影響最顯著的大氣物理量低頻風場、低頻垂直速度場和低頻位勢高度場，而不是從原始的氣象場，這些物理量都可以相對直接地從動力角度來反映準雙周大氣低頻振蕩對褐飛虱遷入的影響。

[1]Kennedy J S A.Turning point in the study of insect migration[J].Nature,1961,189:785-791.

[2]Southwood T R E.Migration of terrestrial arthropods in relation to habitat[J].Biol.Rev,1962,37:171-214.

[3]胡高,包云軒,王建強,等.褐飛虱的降落機制[J].生態(tài)學報,2007,27(12):5068-5075.

Hu G,Bao Y X,Wang J Q,et al.Case studies on the landing mechanisms of the brown planthoppers(St?l)[J].Acta Ecologica Sinica,2007,27(12):5068-5075.(in Chinese)

[4]程遐年,陳若篪,習學,等.稻褐飛虱遷飛規(guī)律的研究[J].昆蟲學報,1979,22(1):1-21.

Cheng X N,Chen R C,Xi X,et al.Studies on the migrations of brown planthopper(St?l)[J].Acta Entomologica Sinica,1979,22(1):1-21.(in Chinese)

[5]包云軒,程極益,程遐年,等.中國盛夏褐飛虱北遷過程的動態(tài)數(shù)值模擬[J].昆蟲學報,2000,43(2):176-183.

Bao Y X,Cheng J Y,Cheng X N,et al.Dynamical and numerical simulations on the processes oflong-distance migration northward during the midsummer in China[J].Acta Entomologica Sinica,2000,43(2):176-183.(in Chinese)

[6]李崇銀.大氣低頻振蕩[M].北京:氣象出版社,1991:1-310.

Li C Y.Atmospheric low frequency oscillation[M]. Beijing: China Meteorological Press,1991:1-310.(in Chinese)

[7]Madden R D,Julian P.Detection of a 40-50 day oscillation in the zonal wind in the tropical Pacific[J]. Atmos. Sci, 1971, 28: 702-708.

[8]Madden R D,Julian P.Description of globe scale circulation cells in the tropics with 40-50 day period[J]. Atmos. Sci, 1972, 29:1109-1123.

[9]徐國強,朱乾根.大氣低頻振蕩研究回顧與概述[J].氣象科技,2003,31(4):193-200.

Xu G Q,Zhu Q G.Review of research on atmospheric low-frequency oscillation[J].Meteorological Science and Technology,2003,31(4):193-200.(in Chinese)

[10]高守亭,朱文妹,董敏.大氣低頻變異中的波流相互作用(阻塞形勢)[J].氣象學報,1998,56(6):665-680.

Gao S T,Zhu W M,Dong M.On the wave-flow interaction in the low frequency atmospheric variation:blocking pattern[J]. Acta Meteorologica Sinica,1998,56(6):665-680.(in Chinese)

[11]Sooraj K P,Kim D,Kug J S,et al.Effects of the low-frequency zonal wind variation on the high frequency atmospheric variability over the tropics[J].Clim.Dyn.,2009,33:495-507.

[12]崔雪鋒,黃榮輝.熱帶強對流區(qū)域的位置和水平尺度對大氣低頻振蕩的影響[J].中國科學院研究生院學報, 2002, 19(1): 49-58.

Cui X F,Huang R H.Impact of position and horizontal-scale of strong convection area on the low-frequency oscillation in the tropical atmosphere[J].Institute of Atmospheric Physics of the Chinese Academy of Sciences,2002,19(1):49-58.(in Chinese)

[13]董佩明,劉健文,董婕.用奇異譜分析診斷地面氣溫場低頻振蕩周期[J].氣象,1997,23(11):6-10.

Dong P M,Liu J W,Dong J.Diagnosis of low-frequency oscillation period of the surface temperature field in the Southeast China[J].Meteorological Monthly, 1997, 23(11): 6-10.(in Chinese)

[14]陸爾,丁一匯.1991年江淮特大暴雨與東亞大氣低頻振蕩[J].氣象學報,1996,54(6):730-736.

Lu E,Ding Y H.Low frequency oscillation in east Asia during the 1991 excessively heavy rain over Changjiang-Huaihe river basin[J].Acta Meteorologica Sinica,1996, 54(6):730- 736.(in Chinese)

[15]Wang B,Ding Y H.An overview of the Madden-Julian oscillation and its relation to monsoon and mid-latitude circulation[J].Advances in Atmospheric Sciences, 1992, 9(1): 93-111.

[16]Zhang J,Yu S.A diagnostic study on the relationship between the assembling of low frequency waves in the Pacific Ocean and the abnormality of the subtropical high[J].Advances in Atmospheric Sciences,1998,15(2):247-257.

[17]胡堅,王玉清.夏半年大氣低頻振蕩與西北太平洋臺風路徑[J].氣象學報,1992,50(4):420-428.

Hu J,Wang Y Q.Atmosphere and ocean low-frequency oscillation and their effects on tropical cyclone tracks over the Northwestern Pacific[J]. Acta Meteorologica Sinica, 1992,50(4):420-428.(in Chinese)

[18]包云軒,程極益,程遐年.盛夏褐飛虱北遷大發(fā)生的氣象背景:個例分析[J].南京農(nóng)業(yè)大學學報,1999,22(4):35-40.

Bao Y X,Cheng J Y,Cheng X N.Analysis on the meteorological backgrounds for the great events of’ long-distance migration northward during midsummers in China[J].Journal of Nanjing Agricultural University, 1999,22(4): 35-40.(in Chinese)

[19]包云軒,黃金穎,謝曉金,等.季風進退和轉(zhuǎn)換對中國褐飛虱遷飛的影響[J].生態(tài)學報,2013,33(16):4864-4877.

Bao Y X,Huang J Y,Xie X J,et al.Influence of monsoon's advancing,retreating and conversion on migrations of(St?l) in China[J].Acta Ecologica Sinica, 2013,33(16):4864-4877.(in Chinese)

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

Bao Y X,Li J J,Miao Q L,et al.Simulation of atmospheric dynamical background for a great migration event of the brown planthopper(St?l):a case study[J].Chinese Journal of Agrometeorolog, 2008,29(3): 347-352.(in Chinese)

[21]Bao Y X,Ding W W,Xie X J,et al.Statistical analysis on the influence of the landfalling strong tropical cyclones in the catastrophic immigrations of(st?l) in China[J].Journal of Tropical Meteorology,2014,20,(1):8-16.

[22]侯婷婷,霍治國,李世奎,等.影響稻飛虱遷飛規(guī)律的氣象環(huán)境成因[J].自然災害學報,2003,12(3):142-148.

Hou T T,Huo Z G,Li S K,et al.Causes of meteorological environment influencing on migration of rice planthopper [J].Journal of Natural Disaster, 2003,12(3): 142-148. (in Chinese)

[23]馬飛,許曉風,翟保平,等.害蟲災害研究的復雜性理論框架[J].昆蟲知識,2003,40(4):307-312.

Ma F,Xu X F,Zhai B P,et al.The framework of complexity theory in insect pest disaster research[J].Entomological Knowledge,2003,40(4):307-312.(in Chinese)

[24]吳洪寶,吳蕾.氣候變率診斷和預測方法[M].北京:氣象出版社,2005:208-225.

Wu H B,Wu L.Diagnosis and prediction of climate variability[M]. Beijing:China Meteorological Press, 2005: 208-225.(in Chinese)

[25]魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預測技術[M].北京:氣象出版社,2009:1-296.

Wei F Y.Modern weather forecasting techniques and statistical diagnosis[M].Beijing:China Meteorological Press, 2009:1-296.(in Chinese)

[26]Draxler R R,Hess G D.An overciew of HYSPLIT_4 modeling system for trajectories dispersion and deposition [J].Australian Meteorological Magazine, 1998, 47: 295-308.

[27]唐廣田,包云軒.廣西褐飛虱發(fā)生特點及其遷飛路徑分析[J].中國農(nóng)業(yè)氣象,2015,36(1):74-82.

Tang G T,Bao Y X.Analysis on occurrence pattern and migratory paths of brown plant-hopper in Guangxi province[J].Chinese Journal of Agrometeorology, 2015, 36(1):74-82.(in Chinese)

[28]齊會會,張云慧,程登發(fā),等.褐飛虱2009年秋季回遷的雷達監(jiān)測及軌跡分析[J].昆蟲學報,2010,53(11):1256-1264.

Qi H H,Zhang Y H,Cheng D F,et al.Radar observation and trajectory analysis on the autumn return migration of the brown planthopper,(Homoptera: Delphacidae) in 2009 in China[J].Acta Entomologica Sinica,2010,53(11):1256 -1264.(in Chinese)

[29]黃奕武,趙偉.2010年9月大氣環(huán)流和天氣分析[J].氣象,2010,36(12):116-121.

Huang Y W,Zhao W.Analysis of the September 2010 atmospheric circulation and weather[J].Meteorological Monthly, 2010,36(12):116-121.(in Chinese)

[30]李麗平,許冠宇,柳艷菊.2010年華南前汛期低頻水汽輸送對低頻降水的影響[J].熱帶氣象學報,2014,30(3):423-431.

Li L P,Xu G Y,Liu Y J.Influences of low-frequency moisture transportation on low frequency precipitation anomalies in the annually first rain season of south China in 2010[J].Journal of Tropical Meteorology, 2014,30 (3): 423-431.(in Chinese)

[31]Akira O,Tomonari W,Yoshito S,et al.Real?time prediction system for migration of rice planthoppers(horváth) and(st?l) (:)[J].Japanese Society of Applied Entomology and Zoology,2005,40(2):221-229.

[32]齊會會,張云慧,蔣春先,等.桂東北稻區(qū)第七代褐飛虱遷飛規(guī)律及蟲源分析[J].生態(tài)學報,2014,34(8):2039-2049.

Qi H H,Zhang Y H,Jiang C X,et al.Migration of the 7th generation of brown planthopper in northeastern GuangxiZhuang Autonomous Region,and analysis of source areas[J].Acta Ecologica Sinica,2014,34(8):2039-2049.(in Chinese)

[33]劉垚,包云軒,魏巍,等.復雜地形下褐飛虱遷飛的數(shù)值模擬:個例研究[J].生態(tài)學報,2016,36(16):5263-5275.

Liu Y,Bao Y X,Wei W,et al.Simulations of the migration of the brown planthopper,(st?l) ,across the complex terrain of China:a case study[J].Acta Ecologica Sinica,2016,36(16):5263-5275.(in Chinese)

Impacts of Quasi-biweekly Atmospheric Low Frequency Oscillation on Immigration of（St?l）in Cropping Area of Jiangling

DIAO Jia-min1,2, SUN Si-si1,2, BAO Yun-xuan1,2, LU Ming-hong3, XIE Xiao-jin1,2

（1.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China; 2.Key Laboratory of Agricultural Meteorology in Jiangsu Province/Nanjing University of Information Science & Technology, Nanjing 210044; 3.National Agricultural Technology Extension and Service Center, Ministry of Agricultural, Beijing 100081）

Understanding the impacts of the atmospheric oscillation with quasi-biweekly low frequency on the immigration of brown planthopper (BPH),(St?l) could provide an important reference to the extended-period forecast of its occurrence. The daily BPH’s light trap data from 105 plant protection stations in China, the daily reanalyzed meteorological data from the National Center of Environment Prediction (NCEP), and the data from the National Center of Atmospheric Research (NCAR) from 1979 to 2015 were collected. It took 2010 as a typical year with occurrence characteristic of BPH, and the relationship between the amount of BPH migrated in rice-growing region from Nanling Mountain to the Yangtze River and the atmospheric oscillation with quasi-biweekly low frequency was analyzed. HYSPLIT (hybrid single particle Lagrangian integrated trajectory) model was used to analyze the backward trajectory of the brown planthopper population at the typical site. The impacts of the atmospheric oscillation on insect sources and migratory paths of BPH was discussed. The results showed that, （1）10-20d oscillation cycle of BPH amount existed significantly in the rice-growing region from the Nanling Mountain to the Yangtze River during August and September in 2010. （2）Wind fields, geo-potential height fields and vertical velocity fields on the isobaric surfaces of 850hPa and 925hPa appeared significant low-frequency cycles of 10-20d simultaneously. And the low-frequency oscillation of BPH’s immigrating amount was correlated significantly with these atmospheric physical variables. （3）Under the impacts of the low-frequency southwestern wind, the populations of BPH were transported from the rice-growing region of southwestern China to the rice-growing region, from the Nanling Mountain to the Yangtze River. In the later period, the part of BPH’s populations immigrated from the rice-growing region between the Huaihe River and the Yangtze River to the rice-growing region between the Nanling Mountain and the Yangtze River under the impacts of the low-frequency northeastern wind.（4）Results of backward trajectory analysis of Brown Planthopper Population was provided to confirm that the insect sources and migration paths of BPH’s populations that immigrated into the rice-growing region between the Nanling Mountain and the Yangtze River were consistent with the low frequency prevailing wind fields.

(St?l); Wavelet analysis; Low-frequency atmospheric airflow field; Migration trajectory; Region of insect source

10.3969/j.issn.1000-6362.2017.03.006

2016-08-25

通訊作者：。E-mail：baoyunxuan@163.com

國家自然科學基金（41475106；41075086）；江蘇省農(nóng)業(yè)科技自主創(chuàng)新基金項目（CX(12)3056）；江蘇省高校自然科學研究項目（14KJA170003）；江蘇高校優(yōu)勢學科建設工程資助項目（IRT1147）

刁家敏（1991-），碩士生，研究方向為農(nóng)業(yè)氣象、農(nóng)業(yè)病蟲測報學。E-mail：1121967219@qq.com

刁家敏,孫思思,包云軒,等.準雙周大氣低頻振蕩對江嶺稻區(qū)褐飛虱遷入的影響[J].中國農(nóng)業(yè)氣象,2017,38(3):181-195