云南省褐飛虱早期遷入蟲源及其發生大氣背景分析

包云軒，尚 潔，孫思思，謝曉金，陸明紅，劉萬才

1 南京信息工程大學氣象災害預報和評估協同創新中心/南京信息工程大學，南京 210044 2 江蘇省農業氣象重點實驗室/南京信息工程大學，南京 210044 3 農業部全國農業技術推廣與服務中心，北京 100125

褐飛虱，Nilaparvatalugens(st?l)，是東亞和東南亞地區水稻生產上危害性最強的害蟲之一，通過直接吸取水稻的液態營養物質和傳播水稻草叢狀矮縮病、鋸齒矮縮病等病毒而導致水稻嚴重減產[1]，具有群聚性、遠距離遷飛性和災變突發性等特點[2]。褐飛虱的生存和繁殖需要合適的溫度和水稻寄主，為了尋找合適的棲息地，褐飛虱每年在東南亞和東亞地區進行季節性的南北往返遷飛[3- 5]。每年3月—5月褐飛虱從東南亞地區開始北遷進入中國境內的華南和西南稻區，繁殖一代后繼續向北遷飛擴散，從南向北每增加兩個地理緯度，褐飛虱的年發生減少一個世代，其中云南地區的褐飛虱可繁殖8—9代。從4月下旬到8月上旬，褐飛虱約發生5次重大北遷過程，繁殖的后代可覆蓋中國乃至日本和朝鮮半島的整個東亞水稻生長區。8月底9月初，褐飛虱開始向南回遷至溫度較高的地區，發生2次重大南遷過程；10月底，大多數褐飛虱種群遷飛至境外可以安全越冬的地域[3,6- 9]。

中國的溫帶、亞熱帶和熱帶大部分地區無法為褐飛虱提供合適的越冬條件，每年的初始褐飛虱遷入種群大多為3月底從中南半島北遷進來的境外蟲源。云南省地處低緯度高原地區，與越南、老撾、緬甸等地接壤(圖1)，錯綜復雜的地形地貌使得省內溫度的水平、垂直分布差異較大，也造就了云南水稻種植制度的多樣化。云南中部和東北部為一季中稻區，西南部和南部是雙季稻區，早中晚稻混栽，終年有水稻生長[10]，這為褐飛虱的北遷提供了棲息地，每年3月開始褐飛虱陸續遷入云南南部和東南部地區，趨集于當地的早稻和早栽中稻上繁殖為害，對云南本地的水稻生長造成嚴重危害，并為當年更北緯度的中國境內褐飛虱的爆發成災積累了初始蟲源[11]。云南因其特殊的地理位置，與廣西并列東南亞境外蟲源北遷進入中國的第一站，但是由于地形復雜、科研力量薄弱、測報數據缺乏等原因，對云南省褐飛虱遷飛路徑尚未有系統的研究，本文根據云南省現有的褐飛虱燈誘數據記錄，對云南省早期遷入種群的遷飛路徑進行分析，追溯境外蟲源地，可以根據蟲源地褐飛虱生物型和抗藥性實現褐飛虱的源頭治理，對當年中國褐飛虱的爆發預警具有重要意義[11]。

目前，對于國內外稻飛虱的軌跡模擬[11- 16]多采用HYSPLIT模式(the Hybrid Single Particle Lagrangian Integrated Trajectory Model)，該模式依賴天氣預報模式提供的風場，追溯遷飛昆蟲的遷出蟲源地、遷飛路徑和降落區。為了更全面地反應不同尺度天氣過程對褐飛虱起飛、空中飛行和降落的動力作用，本文采用WRF-Flexpart耦合模式對褐飛虱的遷飛軌跡進行數值模擬。包云軒等[14]對WRF- HYSPLIT和WRF-Flexpart兩種軌跡模式進行了比較研究，表明WRF-Flexpart模式在計算軌跡的過程中較WRF-HYSPLIT多考慮了對流參數、地表脅迫和各種地形參數，能更客觀地反應各軌跡參數的起伏變化；劉垚等[17- 18]利用WRF-Flexpart耦合模式對復雜地形條件下褐飛虱遷飛進行了數值模擬，表明該模式適用于褐飛虱遷飛軌跡的模擬，并且能較好地反應地表物理過程、大氣湍流結構和地形起伏對褐飛虱種群遷飛的影響。

近年來，由于全球氣候變暖、水稻耕作制度變化和褐飛虱自身抗藥性増加，褐飛虱遷飛規律和區域性降落特征已發生明顯的變化，其遷出蟲源地、遷飛路徑、降落區均變得更為復雜[19]。本文根據云南省褐飛虱燈誘逐日數據，通過分析選取合適的年份、站點，利用WRF-Flexpart耦合模式對選取的個例燈下高峰日進行軌跡模擬，結合遷飛期間的大氣環流形勢和天氣背景分析，明確近年來云南早期褐飛虱遷入種群的蟲源地分布狀況，探討種群遷入和降落的大氣背景，以期為褐飛虱的跨境監測、源頭治理提供科學依據和技術支撐。

1 材料與方法

1.1 資料來源

(1)蟲情資料：褐飛虱蟲情數據由全國農業技術推廣服務中心提供，整理可得到云南省褐飛虱逐日燈誘數據。

(2)氣象資料：為美國國家環境預報中心(NCEP)提供的FNL氣象再分析資料，空間分辨率為1°×1°，時間間隔為6 h。

(3)地理信息圖層：STRM90 m高程數據來源于中國科學院計算機網絡信息中心地理空間數據云平臺(http://www.gscloud.cn)。

1.2 研究方法

1.2.1 個例選取

圖1 研究區域地形概況及站點分布圖

褐飛虱在云南省南部極少數地區存在越冬蟲源[20]，這些蟲源不足以在廣大稻區內爆發成災，每年7月下旬主害代的發生情況取決于4—5月跨境遷入代的蟲源基數，前期遷入峰早、遷入峰多、遷入量大往往會導致褐飛虱的大發生[21- 22]。

云南省地處中國西南地區，因地勢復雜、科研力量薄弱等原因導致云南省褐飛虱燈誘測報數據缺乏，從1999開始才陸續有站點記錄數據，2009年有記錄數據的站點增加至15個，2012年最多達到22個站點，且有站點缺測情況，所涉及的26個觀測站點的分布情況如圖1所示。考慮到數據在時間上的完整性和空間上的代表性，根據2009年發布的《稻飛虱測報調查規范》[23]對燈下高峰期及高峰日的定義，結合實際蟲情發生情況，以4—5月份期間褐飛虱單燈單日蟲量過100頭、連續幾年有有效數據、站點基本覆蓋云南全省為標準，篩選出了勐海、芒市、江城、師宗、麻栗坡、廣南6個代表性站點。結合圖1可見，空間上，選取的站點分布在云南南部、東部以及西部，位于北部的站點每年的有效燈誘數據少于10頭，可以認為對本省褐飛虱分布及其危害影響不大而不予考慮，因此所選站點基本覆蓋發生區每個方位，具有一定的代表性；時間上，因為2014年云南省褐飛虱早期蟲量整體偏少，故選取2013年、2015年和2016年3年的數據，詳細的站點情況及個例蟲情參見表1。

表1 用于云南省蟲源地追溯的6個代表性站點的蟲峰簡況

1.2.2 WRF模式及參數化方案

WRF(Weather Research and Forecasting Model)模式是由美國環境預測中心(National Center of Environmental Prediction in USA, NCEP)和美國國家大氣研究中心(National Center of Atmospheric Research in USA, NCAR)等科研機構聯合開發的一種統一的中尺度天氣預報與研究模式[24]。該模式為完全可壓縮及非靜力模式，采用Fortran90語言編寫，模式水平方向采用Arakawa C網格點(重點考慮1—10 km)，垂直方向采用地形跟隨質量坐標，在時間積分方面采用三階或者四階的Runge-Kutta算法，是改進從云尺度到各種不同天氣尺度的重要天氣特征預報精度的工具。本文應用WRF V3.7 模式版本，設置雙向三重嵌套，模擬采用NECP提供的1°×1°分辨率的FNL全球氣象再分析格點數據作為初始場，積分時間步長為180 s，結果每1 h輸出1次，以最內層嵌套區域輸出的物理量場進行分析計算，模式配置及參數化方案見表2。

表2 WRF模式配置及方案選擇

1.2.3 Flexpart模式及參數設置

對于大氣流場的變化，大多數數值模式采用歐拉法來計算，但是歐拉法側重于整個流場的變化，忽視了其中流點的個體變化和移動，因此只能不斷向前模擬，難以探究流點的起源。此外，當研究的對象沒有足夠密集、足夠細小到能稱之為微氣塊時(如：褐飛虱)，本身就是針對連續流體而設計的歐拉法更難勝任。Flexpart模式是一種拉格朗日粒子擴散模式，由挪威大氣研究所(Norwegian Institute for Air Research，NIAR)開發，適用于大范圍的大氣傳輸和擴散過程[25- 27]，可以通過后向運算來確定對于測定排放源有潛在影響的分布區域。不同于歐拉計算方法，拉格朗日法以研究單個流體質點運動過程作為基礎，綜合所有質點的運動，再去構建整個流體的變化。Flexpart模式之所以有比較好的模擬效果，是因為它并沒有采用復雜的化學過程，而是使用數值模式中常用的參數化方案思想，結合拉格朗日法的動力學特點。進而實現了對示蹤物準確、高效的模擬。最初的Flexpart模式采用歐洲中期天氣預報中心(European Centre for Medium-Range Weather Forecasts，ECMWF)預報數據作為輸入場來構建大氣環流背景，但是由于數據的分辨率較低。很多時候對中尺度甚至小尺度范圍內的研究并不適用，于是就有了WRF-Flexpart耦合模式。

WRF-Flexpart耦合模式使用WRF模式輸出的氣象背景場驅動Flexpart模式，設置褐飛虱降蟲地為起點，采用后向模擬，可以獲得不同站點釋放褐飛虱的連續位置，將輸出的聚類軌跡導入ArcGIS中，利用GIS空間分析法對不同站點的各時刻褐飛虱位置點進行擬合，即可得到褐飛虱的遷飛軌跡，追溯到褐飛虱的蟲源地[21]。本文采用Flexpart-WRF v3.7模式[28]，計算云南春季褐飛虱回推12、24 h和36 h的軌跡，模式釋放的起止時間、釋放點的經緯度、釋放種群數量等參數由站點具體情況設定，模式參數設置詳見表3 。褐飛虱春、夏季遷出高峰通常出現在日出前或日落后，呈現“晨昏雙峰型”[3,29- 30]。經查詢，2013、2015和2016年6個站點在4—5月份的日出時刻介于06:13—7:08之間，日落時刻介于19:27—20:05之間，因此模式假設：(1)褐飛虱在日出前1 h內和日落后1 h內起飛，設置釋放時間為當日05:00至前一日21:00之間，飛行時間持續36 h[29,31]；(2)褐飛虱遷飛高度因季節不同，夏季飛行高度集中在1500—2000 m[31]，最高高度不超過高空中15℃層所在高度，各站點的釋放高度結合該站點的海拔高度確定；(3)軌跡終止點必須在水稻種植區，且該稻區內水稻生長正處于生長后期；(4)回推軌跡落點區的褐飛虱正值遷出高峰。

表3 WRF-Flexpart 模式耦合計算參數設置

2 結果與分析

2.1 高峰日回推軌跡

根據云南省2013、2015和2016年3年的褐飛虱逐日燈誘數據，結合褐飛虱相關生物學和行為生態學特性，選取軌跡計算相關參數，如起飛、降落時間和遷飛高度等，利用WRF-Flexpart耦合模式模擬了6個站點高峰日不同高度層的遷飛軌跡(圖2)。圖2中用不同顏色代表不同年份、不同釋放高度的遷飛軌跡，為了便于分析，每個站點同一年同一釋放高度線條顏色保持一致；根據各站點的海拔高度，芒市設置1000、1500 m和2000 m 3個釋放高度，師宗設置2000 m釋放高度，其余站點均設置1500 m和2000 m釋放高度；軌跡中的點從站點開始依次每小時回推，累積回推出12、24 h和36 h前的蟲源地，據此可得到各站點褐飛虱發生高峰日回推12、24、36 h的蟲源地。6個代表性站點的高峰日蟲源地基本可以代表當年云南早期褐飛虱的來源，結合相應時段的天氣背景，可以分析褐飛虱的降落機制。

圖2 2013年、2015年、2016年各站點高峰日褐飛虱回推36 h遷飛軌跡

芒市位于云南省西部，3年早期遷入的褐飛虱均來自緬甸西部，不同高度層從降蟲區到蟲源地之間軌跡走向基本重合，說明褐飛虱遷飛基本上受氣流脅迫作用所致：12 h和24 h回推軌跡落點分別處于緬甸中部、西部，36 h的回推軌跡略有不同，其中2013年和2015年36 h軌跡回推至緬甸西部，2016年36 h軌跡回推落點在孟加拉國，從某種程度上孟加拉國存在直接為我國提供蟲源的可能；勐海位于云南的南部，3年褐飛虱高峰日回推軌跡落點緬甸東部到西部均有，12 h的軌跡落點在緬甸東部，24 h和36 h蟲源地顯示為緬甸中部、中部偏西，蟲源軌跡基本一致；江城位于云南南部，褐飛虱軌跡隨釋放高度的變化基本不發生偏移，2013年12、24 h和36 h 3個時次的蟲源地依次為緬甸東部、中部和中西部，2015年為緬甸東部、東南部和西南部，2016年分別為老撾北部、老撾西北部和泰國西北部，蟲源地隨時間延長向東南方向發生偏移。和江城類似，位于云南東部的師宗3年高峰日回推12、24 h和36 h 3個時次的蟲源地各不相同，2013年依次為云南中南部、緬甸東部偏北和緬甸中部，2015年為云南中南部、緬甸東部和和緬甸西南部，2016年分別為越南北部、老撾東北部和泰國東北部。麻栗坡蟲源地隨時間延長向西南方向發生偏移：2013年蟲源來自越南，2015年12 h和24 h蟲源來自越南北部、36 h蟲源地落點為云南西南部，2016年12、24 h和36 h回推落點依次為越南北部、云南南部和緬甸東部。廣南2013年24、36 h蟲源分別來自緬甸東部和中部，2015年36 h蟲源落在緬甸東部，其余時段的蟲源均來自云南省內。綜上所述，近年來云南地區早期褐飛虱境外蟲源主要來自緬甸，部分來自老撾和泰國，少數來自越南，還有極少量來自孟加拉國。

云南地形起伏變化明顯，地表脅迫作用對褐飛虱遷飛高度有一定的影響，WRF-Flexpart模式考慮了低層大氣的湍流結構、地表物理過程和地形環境對不同高度三維流場的影響。本文根據站點的海拔高度，分別模擬了1000、1500、2000 m 3個不同釋放高度的褐飛虱回推軌跡，并對遷飛高度的變化進行了擬合(圖3)。可以看出，遷飛高度隨釋放高度的增加而增加，在釋放高度分別為1500 m和2000 m時，褐飛虱遷飛高度的起伏趨勢表現出較好的一致性，這表明在1500 m至2000 m這個高度層內，系統性垂直氣流對不同褐飛虱種群遷飛的影響得到了較一致的反映。值得注意的是，在芒市1000 m釋放高度的時候，遷飛高度的起伏變化較大，這說明模式不僅較好地反映了系統性垂直氣流的影響，而且更客觀真實地反映了下墊面起伏對其上層三維流場和種群遷飛的影響。此外，經過計算分析，當釋放高度為1500 m時，褐飛虱遷飛的平均高度和盛行高度分別為2167.1 m和2216.3 m；而當釋放高度為2000 m時，這兩個值分別為2454.5 m和2489 m，說明WRF-Flexpart模式在計算遷飛軌跡的時候考慮了不同釋放高度的差異(圖3)。

2.2 降蟲期大氣背景場分析

褐飛虱體積微小，自身飛行能力只有0.3 m/s左右，通過適應和利用大氣環境來完成種群的遷移和繁殖[32]。在合適的大氣環流形勢和天氣條件下，褐飛虱能夠主動升空、在高空隨氣流集群遠距離遷飛，其遷飛活動受到氣象條件的顯著影響。云南省地處青藏高原南麓，向南延伸為中南半島，西臨孟加拉灣，東望中國南海，特殊的海陸分布以及復雜的地貌形成了云南特殊的立體氣候背景，既受到東亞季風的影響，也受到西南季風的影響，還不時受到來自南海和孟加拉灣熱帶氣旋的侵擾。

圖3 各站點高峰日當日09:00開始逆推36 h的不同時次褐飛虱遷飛高度的變化曲線

通過對2013、2015和2016年云南6個站點高峰日降蟲的大氣背景進行分析，發現云南省早期褐飛虱遷入發生的降蟲期的大氣背景大體相似，故本文選取2013年5月22日至26日的一次遷飛過程進行大氣背景分析，背景場包括位勢高度場、風場、垂直速度場、溫度場和相對濕度場。因云南屬高原地形，平均海拔在1000 m以上，地面氣壓接近850 hPa，因此，主要對850 hPa高度層進行天氣背景分析。

2.2.1 位勢高度場

從大尺度大氣環流來看，中國地處亞洲季風區，南亞季風和東亞季風的南北進退為褐飛虱通過長距離遷飛進出中國和東亞其他地區提供了便利條件。結合2013年5月21日至26日的位勢高度場(圖4)可以看出，5月份東南亞地區主要受印緬低壓控制或影響，蟲源區位于低壓系統的東部或東側，從地面到低層大氣均盛行強上升氣流，有利于褐飛虱種群的起飛遷出；同時高空有西南氣流盛行，有利于起飛后的褐飛虱種群從蟲源區向東北方向的中國西南稻區輸送。

圖4 2013年5月21—26日850 hPa天氣形勢

2.2.2 風溫場

風和溫度是影響褐飛虱遷飛活動的兩個主要因素，褐飛虱因其主動遷飛速度遠小于空中水平氣流的速度，因而只能順風遷飛，遷飛方向與風向保持一致。溫度是決定褐飛虱能否起飛的關鍵因子，一般來說，褐飛虱在溫度處于15—25℃左右最有可能起飛，上升到一定的溫度層(18℃層所在髙度)之后開始進入平飛階段，空中飛行的最低氣溫閾值為15℃，低于此溫度便可能迫降或降低飛行速度，甚至直接降落[31]。

圖5給出了2013年5月22日—26日850 hPa云南地區20：00的風溫場分布，從風向上來看，850 hPa高度上東南亞地區主要盛行西南風和東南風，緬甸地區主要盛行西南風，可以為蟲源輸送至云南南部提供動力，這與蟲源地回推軌跡(圖2)的模擬結果相吻合。從溫度條件上看，22日、24日、26日主遷峰期間云南全省溫度高于20℃，而西南地區的四川西南部、貴州大部等地受低溫脅迫，形成水平方向以東北或偏東風為主、垂直方向以下沉氣流為主的“低溫屏障墻”，使得空中遷飛種群遇“墻”降落或“反彈”降落。因此，此次遷飛過程的主要遷入地為云南地區。26日之后，隨著低溫區逐漸向南、向西加強擴展，云南東南部和中部等地區逐漸受到低溫脅迫作用，不太有利于褐飛虱的遷入，本次遷飛過程基本結束(圖5)。

圖5 2013年5月21—26日20時云南850 hPa風溫場

2.2.3 氣溫、垂直速度和相對濕度的垂直分布

褐飛虱的降落行為是主觀和客觀因素綜合作用的結果，如果能尋找到適合取食和繁衍的水稻田，褐飛虱將主動降落；如果蟲匯區遷飛層至地表有下沉氣流存在，或者有降雨的沖刷，或者種群飛行時間過長且蟲體脂肪能消耗將盡，或者種群飛行所經區域的溫度低于15℃等褐飛虱則容易產生迫降[29]。

從圖6、圖7中可以看出，850 hPa(約1500 m)褐飛虱最適宜遷飛層的溫度都不存在低溫脅迫。3個遷入高峰期遷飛過程發生區內，種群能遷飛的溫度層上限(15℃層)變化于940—2500 m，平均高度為2000 m；最適宜遷飛層(18℃層)變化于950—2080 m，平均高度為1462 m；遷飛層下限(21℃層)變化于150—1800 m，平均高度為987 m。從蟲源區到降蟲區。雖然遷飛層頂、最適宜遷飛層和遷飛層底高度在中途有一些起伏變化，但總體趨勢是遷飛高度在逐漸由高走低并最終降至300 m以下，直至迫降地面。

從圖6、圖7中還可以看出，垂直速度場上，在(16—22°N，94—102°E)范圍內，先后存在著不同強度的上升氣流(圖中虛線)；在(23—25°N，100—105°E)范圍內，先后存在著不同強度的下沉氣流(圖中實線)。而前者對應著蟲源區，后者對應著降蟲區，這表明蟲源區的上升氣流促使了蟲源的遷出、而降蟲區的下沉氣流促使了種群的降落，垂直氣流在此次遷飛過程中起了明顯的正向促進作用。

圖6 沿25°N氣溫(℃)和垂直速度(cm/s)的緯向剖面圖

圖7 沿105°E氣溫(℃)和垂直速度(cm/s)的經向剖面圖

有研究表明：當夏季稻田相對濕度在75%以上或春秋季稻田相對濕度在40%以上時，有利于褐飛虱種群的增長[3]；而空中夏季相對濕度在60%以上或春秋季相對濕度在40%以上時，均有利于褐飛虱種群的遷飛[33]。從圖8、圖9所示的相對濕度場上來看，此次遷飛過程中，以溫度場15℃層作為遷飛頂層向下至地表的各時次各區域(圖6、圖7)三維空間內，相對濕度均大于50%，因此，相對濕度對此次過程褐飛虱種群的遷飛不構成水分脅迫。

圖8 沿25°N相對濕度(%)的緯向剖面圖

2.3 蟲源地水稻生育期狀況

云南地處青藏高原的東南側，西部接壤緬甸地區，東部和南部緊靠老撾、越南，常年受東亞季風和西南季風的雙重影響，獨特的地理位置和氣候條件，使得云南早期遷入蟲源主要來自緬甸和金三角地區，緬甸是云南春季褐飛虱遷入的主要蟲源區。

圖9 沿105°E相對濕度(%)的經向剖面圖

素有“稻米之國”美譽的緬甸，是世界上第7大水稻生產國，土地肥沃，地形平坦，大部分地區屬熱帶熱帶季風氣候。就氣候條件來說，緬甸水稻可周年生長，能栽植三熟稻，但因雨水分布和水利設施等原因，多數地區每年只栽植一季稻，在灌溉水源良好或雨季后有積水存在的地方栽植二季稻，有較好人工灌溉設施的地方才栽植三季稻。歷史上，緬甸水稻主要種植一季，生產季節從5月至11月，稱為“雨季稻”，自1992年以來，緬甸政府采取興修水利、減輕田稅等多種措施，大力推廣夏季稻，以充分利用12月至次年4月較好的天氣進行旱季稻的栽種。通過多年來緬甸政府對土地政策的調整和大力投資水利水庫建設，緬甸的旱季稻實現從無到有的質的飛躍，每年1月份播種，5月份收割，2007年旱季稻的種植面積達到23萬hm2[34]。

每年4—5月為緬甸旱季稻的黃熟期，為了尋找新的取食環境和棲息地，褐飛虱種群會大規模地起飛遷出，而此時正是云南一季中稻的移栽時段，該時段的水稻可以為褐飛虱的北遷提供食料和棲息地。加上有利的大氣環流形勢和溫度條件，導致了每年4、5月份東南亞境外蟲源北遷進入云南為害，并為當年褐飛虱的發生積累了初始蟲源。

3 結論與討論

云南省毗鄰緬甸、老撾、越南，與水稻種植大國泰國距離也很近，是境外褐飛虱進入中國西南稻區的第一站，研究云南褐飛虱的蟲源地具有重要意義。近年來有些學者對云南稻飛虱發生規律和遷飛軌跡進行了研究[9- 10,12,21- 22,35]，但研究多根據氣象再分析場資料，利用HYSPLIT模式進行軌跡模擬，且研究對象大多為稻飛虱中的白背飛虱。隨著氣候變暖、水稻種植制度的改變和抗藥性的上升，褐飛虱的蟲源分布和降落結構相發生了復雜的變化，為了探明這些變化，本文采用WRF-Flexpart耦合模式模擬得到褐飛虱遷飛的相關參數，包括軌跡點所在經緯度、高度和該點所處的大氣背景場等，將模式輸出的聚類軌跡進行整合，再利用ArcGIS空間分析功能對各后向逆推時間段褐飛虱的三維遷飛坐標和路徑數據進行曲線擬合，得到了不同時間段褐飛虱的聚類遷飛軌跡及不同釋放高度上褐飛虱遷飛高度的變化特征，并根據褐飛虱的遷飛軌跡判斷出褐飛虱的境外蟲源地。WRF-Flexpart耦合模式在計算遷飛軌跡的時候，考慮了地表物理過程、大氣湍流結構和地形起伏對褐飛虱種群遷飛的影響，通過一系列模擬和分析得到如下結論：

(1)云南省4—5月份早期褐飛虱遷入的主要蟲源區是緬甸，部分來自老撾和泰國，少數來自越南，還有極少量來自孟加拉國。沈慧梅等[12,34]分析了2009年云南省白背飛虱早期遷入種群的蟲源地范圍，得到緬甸是云南春季白背飛虱遷入的蟲源地的結論，這與本文中春季褐飛虱遷入主要蟲源地為緬甸的結論是一致的。褐飛虱和白背飛虱均屬同翅目飛虱科，但是兩者在發生規律方面存在顯著差異。白背飛虱初代蟲源由南方熱帶稻區隨氣流逐代逐區遷入，發生時間要早于褐飛虱；白背飛虱在水稻分蘗盛期、孕穗、抽穗期發生較嚴重，而褐飛虱在抽穗揚花期發生最嚴重；最重要的是，褐飛虱的耐寒能力弱于白背飛虱，白背飛虱在我國海南南部和云南最南部地區終年繁殖，而褐飛虱在我國只有海南省南部可越冬，因此4—5月份云南省早期遷入褐飛虱大多數為境外直接進入蟲源。

(2)這一時期，當釋放高度分別為1500、2000 m時，褐飛虱的遷飛高度分別集中在2216 m和2489 m，平均遷飛高度分別是2167、2454 m，不同釋放高度褐飛虱的回推遷飛高度的起伏趨勢具有較好的一致性。但在釋放高度較低的時候，遷飛高度的起伏變化較大，這表明模式不僅較好地反映了系統性垂直氣流的影響，而且客觀真實地反映了下墊面起伏對其上層三維流場和種群遷飛的影響。這與我們前期對國內褐飛虱遷飛過程研究所得的結論相吻合[14]。

(3)風場切變、低溫屏障、氣流垂直擾動和降水等因素是造成褐飛虱集中降落的直接原因。本文選取2013年5月22—26日發生在云南勐海、麻栗坡和廣南的一次典型褐飛虱遷入過程，分析了大氣背景場對遷飛的影響，結果表明：受印緬低壓控制或影響，高空有從境外蟲源區至云南降蟲區有西南水平氣流作為種群輸送動力；溫場在降蟲區東北側有“低溫屏障墻”存在；垂直速度場上蟲源區有上升氣流促使種群起飛遷出，降蟲區有下沉氣流促使種群降落；弱降水形成的拖曳下沉氣流對降蟲十分有利；相對濕度場對此次遷飛過程不形成任何脅迫。因此，此次過程的大氣動力、熱力和水分條件都有利于云南稻區褐飛虱蟲源的遷入和降落。

近年來隨著東南亞地區水稻灌溉條件的改善，老撾、緬甸等國家北部的一些地區在1—4月份的旱季也可以順利種植水稻，而5月份中下旬云南水稻處于一季中稻移栽階段和分蘗初期，為褐飛虱的北遷提供了豐富的食料和適合的棲息地，加上合適的大氣背景，為褐飛虱的早期遷入提供了有利條件。2013、2015和2016年云南地區早期(4—5月)褐飛虱遷入的蟲源地與2009年相比發生了一些變化[31]，為了有效防治褐飛虱遷入在西南稻區第一站云南的危害，必須根據蟲源地褐飛虱發生規律和種群遷出動態，開展行之有效的早期預警和短期異地預測，合理防控，以減輕褐飛虱種群遷入下游稻區為害水稻造成的減產。