基于WRF-Flexpart的一次褐飛虱回遷過程模擬研究

劉 垚,包云軒,陸明紅,劉萬才

1 南京信息工程大學江蘇省農業氣象重點實驗室,南京 2100442 南京信息工程大學,氣象災害預報預警與評估協同創新中心,南京 2100443 中國氣象局旱區特色農業氣象災害監測預警與風險管理重點實驗室,銀川 7500024 農業部全國農業技術推廣服務中心,北京 100125

基于WRF-Flexpart的一次褐飛虱回遷過程模擬研究

劉 垚1,2,3,包云軒1,2,*,陸明紅4,劉萬才4

1 南京信息工程大學江蘇省農業氣象重點實驗室,南京 2100442 南京信息工程大學,氣象災害預報預警與評估協同創新中心,南京 2100443 中國氣象局旱區特色農業氣象災害監測預警與風險管理重點實驗室,銀川 7500024 農業部全國農業技術推廣服務中心,北京 100125

褐飛虱是影響亞洲地區水稻生產的遷飛性害蟲,其遠距離遷飛與種群分布會受到大氣環流和局地天氣條件的影響。采用WRF-Flexpart耦合模式、GIS空間分析和Vincenty方位角/距離估計等方法,選取2009年9月30日—10月7日發生在中國15個植保站的褐飛虱遷入過程,并根據模式輸出、GIS分析和軌跡計算結果探討了多種大氣物理脅迫共存下大氣動力場、降水、溫度場和濕度場對褐飛虱秋季回遷過程及其降落蟲量分布的影響。研究結果表明：(1)此次重大遷飛過程中,大氣動力場、降水、溫度場和濕度場在不同地區、不同時間段對褐飛虱種群的遷入和降落起了不同的作用,因而造成了遷入蟲量分布的時空差異。(2)從各蟲源地遷出的褐飛虱種群是向著溫暖而濕潤的地方遷飛的,遷飛方向主要受所經區域盛行風向的影響,遷飛距離由途徑區的風速所決定。(3)當遷飛種群途徑局地降水較多的區域時,初期遷入量并不多,經過短暫的一段時間后降蟲量會突然增加。(4)在一定的溫度范圍內,大氣濕度條件會影響褐飛虱的遷入蟲量：褐飛虱遷入蟲量在近地面溫度高于26℃、相對濕度(925 hPa)大于70%且下沉氣流較強的地區分布最多,而在近地面溫度低于24℃、相對濕度(925 hPa)小于50%、無明顯下沉氣流的地區未發現有褐飛虱的遷入。

褐飛虱；天氣條件；WRF-Flexpart模式；遷飛軌跡；遷入蟲量分布

褐飛虱,Nilaparvatalugens(St?l),具有遠距離遷飛的習性,是我國主要的水稻害蟲之一[1],其遷飛過程和遷入量與大氣環流形勢、天氣系統密切相關。通常,大尺度天氣形勢的變化為褐飛虱遠距離遷飛提供了一定的物理條件；中小尺度天氣系統的演變則決定了種群動態和分布特征[2]。

國內外對褐飛虱北遷過程已有大量的研究成果。當溫度在20℃左右時起飛頻繁,低于17℃時起飛受到抑制,小于10℃極少發生起飛現象[3]。褐飛虱最佳成層飛行的溫度在17℃左右,15℃等溫線限定了褐飛虱成層遠距離遷飛的上限高度,隨著空中溫度的升降,褐飛虱遷飛的高度也隨之升降[4]。氣流場對褐飛虱的起飛與降落有影響,決定著褐飛虱的遠距離遷飛特征[5- 7]。在我國,低空急流多盛行偏南風,對于主降區的褐飛虱有輸送作用,急流內存在適宜褐飛虱北遷的溫濕環境,并且其氣流內的下沉氣流決定了降蟲區的分布[4, 8- 9]。沈慧梅等[10]對云南稻飛虱北遷過程分析發現,低溫脅迫是影響稻飛虱集中降落的外在條件,溫度是稻飛虱遷入的重要因子。吳秋琳等[11]通過研究湖南白背飛虱前期遷入過程,得出降雨、低溫和下沉氣流是導致白背飛虱集中降落的直接原因。Turner[12]和Zhu等[13]采用BLAYER模型模擬褐飛虱從我國遷飛至韓國的遷飛軌跡,并分析了褐飛虱起飛、降落和飛行高度等參數以及不同飛行高度種群密度分布的狀況。Furuno等[14]利用GEARN模型對日本10個站的褐飛虱蟲情數據做后向模擬,分析了從中國遷飛至日本西部的褐飛虱的可能蟲源地及其與氣象因子的關系。

目前,對褐飛虱南向回遷過程的研究還較少,對氣象條件影響褐飛虱南遷的定量化研究還相對缺乏。秋季冷空氣南下,偏北風使得褐飛虱向著溫暖且濕潤的南方地區遷飛[15]。而冷鋒南下時,局部地區降溫或降水會促使褐飛虱向南回遷并且大量降落；地面冷鋒前溫度會降低3—4℃,迫使褐飛虱降低遷飛高度,且褐飛虱分布帶平行于冷鋒[16-17]。在褐飛虱秋季回遷的雷達觀測實驗中,觀測到了褐飛虱起飛大多在傍晚到黃昏之間,起飛后在近地層逆溫層頂以上,大氣邊界層頂以下的高度水平飛行,平均遷飛層頂高度為800—1100 m,很少到1500 m,并且褐飛虱密度的垂直分布出現明顯的成層現象[18-19]。

褐飛虱遷飛并不是由單一的因素所決定,是由多個因素共同影響的,其中一種或兩種因素起主要作用,造成褐飛虱遠距離的遷飛、不同種群動態的變化和分布特征。由于風媒昆蟲與氣溶膠粒子具有相似的性質,本文將利用WRF-Flexpart耦合模式,采用Pender[20]歸納的褐飛虱遷飛參數(起飛時間、飛行溫度閾值、飛行高度和飛行持續時間等)來研究褐飛虱南遷過程中,復雜大氣物理脅迫條件對褐飛虱南遷的影響,以期該研究能對中國褐飛虱的測報與防治、農業防災減災及糧食安全具有重要的參考價值。

1 資料與方法

1.1 研究數據及個例概況

褐飛虱蟲情資料選自我國南方稻區植保站的逐日褐飛虱燈誘數據匯總表,由農業部全國農業技術推廣服務中心提供。褐飛虱燈誘數據的當日記錄為前一日08:00時至當日08:00時(北京時間,下同)的實際遷入蟲量,模擬時間應從前一日開始。

2009年9月30日—10月7日,在我國南方地區出現了大范圍降水過程,期間,福建省沿海地區、江西南部和湖南省西南部地區降水較多且溫度較高,江西省中、北部降水量較少且溫度較低,而湖北省只有零星降水卻溫度較高。本文選取此次天氣過程對褐飛虱遷飛有明顯影響的重大遷飛事件,分析大氣動力場、降水、溫度場和濕度場對褐飛虱遷飛軌跡和降落蟲量分布的影響,旨在對比分析褐飛虱南遷過程中不同時間段、不同區域復雜大氣物理背景對褐飛虱遷飛和降落的影響。

1.2 方法

1.2.1 WRF模式及參數化方案

使用WRF 3.6雙向三重嵌套,模式初始場所用氣象資料來自美國國家環境預測中心(NCEP,National Centers for Environmental Prediction)的FNL氣象再分析資料,空間分辨率為1°×1°,時間間隔為6 h。模擬以(34.01°N,117.22°E)為中心點,垂直方向分27層,積分初始時間為2009年9月30日08:00,終止時間為10月8日07:00,積分時間步長是180 s,模式參數化方案見表1,每小時輸出1次模擬結果。

表1 WRF模式參數化方案

1.2.2 Flexpart模式及參數設置

Flexpart模式是挪威大氣研究所(Norwegian Institute for Air Research)研究開發的大氣擴散模式,該模式基于拉格朗日算法,適用于大范圍的大氣傳輸和擴散過程。通過后向算法,Flexpart模式可以用來分析諸如溫室氣體和顆粒灰塵等的潛在影響區域[21]。

采用WRF輸出的氣象場來驅動Flexpart模式計算褐飛虱遷飛軌跡,采用后向模擬分析大氣物理脅迫對褐飛虱遷飛及其密度分布的影響。模擬設置15個植保站作為釋放點源(即降蟲地),采樣率為1200 s,同步間隔600 s。褐飛虱起飛大多發生在傍晚的日落前后,晚秋溫度較低時,褐飛虱遷出高峰多出現在下午溫度較高時刻[22],且褐飛虱夜間降落的情況多于白天。故本文的后向模擬設置以降蟲地為起點,釋放開始時間為每日18:30,釋放結束時間為前一日07:30,即假設當日18:30至前一日7:30之間均有褐飛虱降落。由于秋季褐飛虱遷飛的最長時長約為35 h左右,故本文設置的模擬最大時長為35 h(即后推35 h后褐飛虱的落點位置為軌跡的最遠起點),如18:30釋放的軌跡群,最長后推35 h,則褐飛虱種群應當是前一日7:30起飛的；而該軌跡群中包括了遷飛歷時分別為1—35 h不等的各種軌跡；以此類推,可通過釋放不同時間間隔的設置(如每隔3 h、1.5 h、1 h、0.5 h、10 min、5 min甚至1 s,本文中以1.5 h為間隔)計算出當日18:30至前一日7:30之間后推最大時長35 h或1—35 h不同時長的褐飛虱遷飛軌跡群,將這些軌跡群的后推終點位置導入ArcGIS中制成空間分布圖。分別對2009年10月2—7日中每日褐飛虱的遷入進行一次后向模擬,則得到遷飛過程No.1—No.6,將這些過程的軌跡群落點分布圖連續排列則得到如圖3所示的遷入蟲量模擬分布圖。

1.2.3 遷飛軌跡和蟲量分布的計算方法

由WRF-Flexpart模式得到褐飛虱遷飛的相關數據(經度、緯度、高度和該位置的氣象要素等)。有效數據判定：后向模擬釋放點(起點)為植保站測點,終止點為褐飛虱的可能蟲源地且有水稻生長的地區；空中飛行期間,其飛行高度不低于16℃溫度層的高度；褐飛虱遷飛過程中不能通過強降水區域或大范圍降水區域。

首先,分別統計不同釋放地點、不同時刻的褐飛虱位置數據,利用ArcGIS 10.0空間分析法[23]分別對各后向逆推時間段褐飛虱遷飛位置和路徑數據進行曲線擬合,即可得到不同時間段褐飛虱的遷飛軌跡(回推軌跡)。其次,把褐飛虱的位置數據轉化為點圖層,并且將模擬區域劃分成0.5°×0.5°網格,通過把褐飛虱點圖層與模擬區域的網格關聯來統計每個網格單元內的褐飛虱數量,從而得出模擬的褐飛虱蟲量密度分布。

1.2.4 遷飛方位角及遷飛距離計算方法

Thaddeus Vincenty基于地球是一個橢球體的假設,根據迭代算法計算已知的兩個球體表面點之間的方位角和橢球面距離(緯度以正北方向為正；經度以正東方向為正)[24-25]。這種算法比基于地球是球形的假設計算得到的大圓距離更加精確。本文中,褐飛虱遷飛方位角和遷飛距離度量的起點為后推模擬的褐飛虱起飛點的位置,終點為模擬的釋放點,即降蟲站的位置。

2 結果與分析

2.1 褐飛虱燈誘蟲量的分布特征

2009年10月1日—6日各植保站的褐飛虱燈誘蟲量分布圖(圖1),褐飛虱向西南方向遷入各地稻區,南方地區各站蟲量均有增減。期間,褐飛虱主要分布于福建東部沿海地區、江西南部和湖南西南部稻區,而10月1日—6日湖北東部、江西中部和北部各站褐飛虱遷入零星出現。2009年10月1日,江西南部(信豐縣)最多達到15000頭,福建褐飛虱遷入蟲量也較多,其東部沿海地區單站日遷入量可達到千頭以上,湖南省西南部地區遷入量僅次于福建稻區,單站日遷入量也達千頭以上。2日,福建各站褐飛虱遷入蟲量有所減少,均低于1000頭,浙江各站褐飛虱遷入蟲量增加。3日,浙江和福建中部褐飛虱蟲量開始減少,福建沿海地區蟲量較多。湖南西南部稻區褐飛虱蟲量增加。4—5日,福建、浙江和湖南各站褐飛虱蟲量均有不同程度的減少,福建沿海地區褐飛虱蟲量仍較多。6日,除湖南西南部站點的褐飛虱蟲量有所增加外,福建、浙江和湖南其他地區各站褐飛虱蟲量均繼續減少。

圖1 2009年10月1日—10月6日中國南方水稻主產區各植保站褐飛虱燈誘蟲量分布Fig.1 Distribution of N. lugens′ light trap catches in the main rice-growing regions of south China during 1st — 6th Oct. in 2009

2.2 氣象要素場對褐飛虱遷飛距離和方位角的影響

褐飛虱的空中遷飛因其主動遷飛速度遠小于空中水平氣流的速度而只能順風遷飛,當風速較大時有利于褐飛虱的遠距離遷飛,當風速達到3—8 m/s時,褐飛虱順風遷飛很明顯,遷飛方向與風向一致[4]。遷飛過程No.1期間(9月30日—10月2日,圖2),冷空氣由我國東北地區和朝鮮半島侵入我國山東省境內。受冷空氣影響,浙江地區以東北風為主,江蘇及其附近地區地面氣溫較低約在21℃左右。河南、湖北和湖南地面氣溫約24℃左右,湖南地區風向以北風和東北風為主。福建、廣西和廣東地區以東北風為主,氣溫均高于24℃。隨著冷空氣繼續向西南方向推進,南方各地區氣溫繼續下降。此時,浙江南部比北部褐飛虱遷入蟲量多(圖1)。福建沿海地區風速較小,大量褐飛虱由江蘇和浙江北部向福建沿海地區遷入。湖南西部風速較大,以東北風為主,褐飛虱種群取東北—西南方向遷入湘西南。各站點褐飛虱遷飛距離和方向見表2,褐飛虱遷飛軌跡以偏北至偏南方向居多。

同樣,遷飛過程No.2—No.6期間(10月1日—7日),隨著冷空氣的南下,我國中部和南部地區大范圍降溫,福建和廣東以西北風-東北風為主,地面氣溫約在20—24℃(圖2)。褐飛虱種群主要從東北方向遷入且遷飛距離較遠,最遠可達600 km以上(表2)。而安徽、江蘇、浙江和江西東北部地面氣溫均低于22℃,浙江地區以西北風-東風為主,褐飛虱的主要遷入方向以偏西至偏南為主。湖南和湖北氣溫持續降低,約在22℃左右,以東北風為主(圖2)。遷入湖南湘陰、洪江和東安的褐飛虱種群的遷飛距離大于200.0 km,而湖南龍山因受山地地形阻擋,風速較小,褐飛虱遷飛距離較短,約為34.0—198.5 km(表2)。

2.3 氣象要素場對褐飛虱遷入蟲量分布的影響

2009年9月30日至10月7日褐飛虱歷次遷飛過程中的遷入蟲量密度分布如圖3。9月30日至10月3日,隨著冷鋒的推進,冷空氣影響我國長江下游地區,華東地區普遍開始降溫。長江下游地區盛行東北風-東風,盛行風促使褐飛虱種群從長江中下游地區逐漸向西南至東南方向遷入湖南、浙江和福建等稻區。湖南省的蟲源主要來自兩湖交界區,最遠可至湖北省東南部。浙江的蟲源主要來自江蘇和上海稻區。降落在福建稻區的褐飛虱種群主要來自浙江北部稻區。10月4日冷鋒后的冷高壓主體南下且逐漸減弱,我國南方地區降溫持續但降溫幅度逐漸減小,東南沿海地區以北到東北風為主且風速較大。因此,福建和浙江地區的褐飛虱遷入量有所增加。10月6日,又一股弱冷空氣從西北方向侵入我國東南部地區,該地區氣溫繼續降低,弱冷空氣對相對偏向西南的湖南等地的氣溫影響相對較小,但對東南部地區影響較大。東南沿海部分地區依舊盛行東南風,但湖南(因盛行風向轉變)、福建和浙江各地褐飛虱遷入量明顯減少。

圖2 地面氣溫場和925 hPa風場Fig.2 Temperature field on the surface and wind field on 925 hPa

站點StationNo.1No.2No.3No.4No.5No.6距離Dis.方向Dir.距離Dis.方向Dir.距離Dis.方向Dir.距離Dis.方向Dir.距離Dis.方向Dir.距離Dis.方向Dir.福清75.8S79.3S479.4SSW518.1SW555.1SW531.4SSW同安74.2S107.8SW492.3SW577.8WSW524.5SW618.5SW霞浦61.1ESE74.9SE323.6SSW331.7WSW264.7SW419.5SW龍海49.0SW28.5SSE463.3SSW527.8SW384.7SW400.1SW新羅74.7SSE113.3SSW294.3SW396.7SW476.3SW451.9SW永安68.4SSE44.2SE216.1SW259.6WSW508.1SW374.2SW仙游73.7SW92.2SW104.3SSW479.6WSW402.9WSW571.7SW湘陰217.2SW237.8SW231.6SW198.0SW277.0SW225.7SW龍山87.8S198.5SW111.1SSW34.0SSW139.8S50.4WSW洪江220.0SW256.3SW249.3SW238.3SW344.0SW231.6WSW東安243.6WSW281.8SW295.8SW248.6SW383.3SW301.3WSW蕭山460.3SE303.2SSW303.0SSW213.8WSW246.0WSW241.6SW溫州124.8SSE274.8S390.4SSW389.6SW415.1SW397.5SW龍游387.4WSW274.0S324.2SSW242.9WSW352.0SW--宣州373.0WSW250.2SSW220.1SSW196.3SW176.9SW118.9WSW

圖3 不同遷飛期褐飛虱遷入蟲量的模擬分布Fig.3 Spatial distribution of N. lugens′ simulated immigrating populations in different migration periods

對褐飛虱遷飛過程No.1(9月30日—10月2日,如圖3)進行了詳細分析,從沿Line01方向的溫度、垂直速度和相對濕度剖面圖(圖4,圖5)上可以看出,福建沿海地區近地面至925 hPa氣溫約在22℃以上,850 hPa氣溫約在16—18℃。該地區有較適宜的溫度和相對濕度條件,尤其在高空925 hPa高度上相對濕度大于70%,更宜于褐飛虱向該地區的遷入。在(117.5°E,24.0°N)至(119.5°E,26.3°N)剖面附近,有較強的下沉氣流,非常有利于褐飛虱的降落,且福建沿海地區降水明顯,降水脅迫的動力拖曳作用也是此次福建地區褐飛虱集中降落的主要原因。而實際上這一時段福建沿海地區各植保站均有較多的褐飛虱降落。從沿Line02方向的溫度、垂直速度和相對濕度剖面圖(圖4,圖5)上可以看出,江西中部和北部稻區近地面氣溫低于24℃,925 hPa溫度低于20℃。該地區相對濕度較小,除中部略偏南的小范圍區域925 hPa相對濕度大于60%以外,其它區域均小于50%,溫度和相對濕度條件都不是很好,不適宜于褐飛虱的遷入。雖然江西中部16℃層(遷飛層頂)可達1750 m(即825 hPa),并有一片弱下沉氣流區,但由于其溫度和下沉氣流條件均不如Line01、Line03所經地區,因此,褐飛虱降落的熱力和動力條件不如福建沿海稻區和湘西南地區。從沿Line03的溫度、垂直速度和相對濕度剖面圖(圖4,圖5)上可以看出,湖南西南部稻區近地面至925 hPa氣溫均在20℃以上,925 hPa(最適遷飛層)氣溫約在20—22℃,且相對濕度約為50%—70%,溫濕條件非常適宜于褐飛虱的遷入；(111.0°E,26.1°N)附近剖面有較強的下沉氣流區。較之湖南東部地區,湖南西部近地面和925 hPa的氣溫更適宜,下沉氣流更強,更適宜于褐飛虱的降落。其它遷飛過程也能得到類似的結論(其溫度、垂直速度和相對濕度的剖面圖略)。

圖4 遷飛過程No.1的溫度和垂直速度剖面圖(9月30日—10月2日,剖線如圖3所示)Fig.4 Patterns of temperature and vertical velocity during 30th Sep. — 2nd Oct. (Section line as shown in Fig.3)

圖5 遷飛過程No.1的相對濕度剖面圖(9月30日—10月2日,剖線如圖3所示)Fig.5 Pattern of relative humidity during 30th Sep. — 2nd Oct. (Section line as shown in Fig.3)

研究發現,褐飛虱遷入蟲量在近地面溫度高于26℃、相對濕度(925 hPa)大于70%且下沉氣流較強的地區分布最多,其次是近地面溫度為24—26℃、相對濕度(925 hPa)介于70%和50%之間且下沉氣流較強的地區。而在近地面溫度低于24℃、相對濕度(925 hPa)小于50%、無明顯下沉氣流的地區未發現有褐飛虱種群的遷入。通過模擬2009年9月30日至10月7日的褐飛虱南遷過程,發現遷飛軌跡和遷入蟲量的模擬結果與同期實際褐飛虱燈誘蟲量的變化有較好的一致性(表3)。

表3 不同區域褐飛虱遷飛環境參數的對比

3 結論與討論

本研究選取2009年9月30日—10月7日在中國南方地區出現大范圍降水的一次褐飛虱重大南遷過程來研究氣象要素場對遷飛的影響。期間,福建省沿海地區、江西南部、湖南西南部降水較多且溫度較高,江西省中部和北部降水量較少且溫度較低,而湖北省只有零星降水卻溫度較高。利用WRF-Flexpart耦合模式模擬了此次天氣過程中褐飛虱遷飛軌跡及其遷入蟲量的空間分布。在褐飛虱遷飛過程No.1—No.6(9月30日—10月7日),湖北地區主要受大氣動力脅迫(無降水)和大氣熱力脅迫的共同作用,有一定量的褐飛虱種群向該地區遷飛和降落。福建沿海地區、江西南部地區和湖南西南部地區受大氣動力脅迫(包括降水)、大氣熱力脅迫和大氣水分脅迫(相對濕度)的共同作用,有大量的褐飛虱種群向這些地區遷飛和降落。而整個遷飛過程中,江西中部和北部地區有利于褐飛虱種群遷入和降落的大氣動力、大氣熱力和大氣水分條件均不是很適宜,雖然在遷飛過程No.1—No.2(9月30日—10月3日)中該地區有少量的降水,但在遷飛過程No.3—No.6(10月2日—10月7日)中該地區卻無降水,該地區的地面及低空(常規遷飛層)溫度一直較低、相對濕度一直較小,不利于褐飛虱種群向該地區的遷飛和降落。

研究結果表明,風和垂直氣流為此次褐飛虱遷飛過程提供了必要的動力條件,是褐飛虱南遷的主要動力因素,風的盛行方向決定了中國褐飛虱秋季回遷的方向,下沉氣流是秋季南遷降蟲的關鍵動力因素。其次是濕度條件,由于福建、江西中北部和湖北三地緯度相近,一般無降水情況下低空和近地氣層相對濕度相差不大,而此次降水過程卻使得3個地區大氣濕度分布相差明顯,導致降蟲分布差異明顯,表明大氣水分對褐飛虱南遷和降落有著重要的脅迫作用。溫度是在春秋過渡季影響褐飛虱遷入和降落的重要條件之一,無論是低空適宜遷飛層還是近地氣層,溫度在秋季回遷過程中的熱力脅迫作用不可忽視,本文的研究還表明：褐飛虱的秋季回遷總是向著溫暖而濕潤的地方遷飛的,溫濕條件共同制約或影響著褐飛虱的遷飛行為,這與江廣恒等[17]、程極益等[26]、胡高等[27]和程遐年等[28]的研究結論一致。本文只選取了一個典型的秋季回遷個例,分析了大氣物理脅迫對褐飛虱遷飛和降落的影響,普遍的規律、可業務化的氣象預警預報指標需要通過大量個例的綜合分析、實驗觀測的驗證等來總結、提取、歸納和推廣試用。

[1] 巫國瑞, 黃次偉, 陶林勇, 馮炳燦, 陳福云, 劉少華, 葉志長, 陸利全, 趙有良. 影響褐飛虱猖獗和為害的因素. 生態學報, 1984, 4(2): 157- 166.

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

[3] 李金建. 大氣背景對褐飛虱災變性遷入影響的個例研究. 南京: 南京信息工程大學, 2008.

[4] Rosenberg L J, Magor J I. Predicting windborne displacements of the brown planthopper,Nilaparvatalugensfrom synoptic weather data. 1. Long-distance displacements in the north-east monsoon. Journal of Animal Ecology, 1987, 56(1): 39- 51.

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

[6] 包云軒, 程極益, 程遐年, 儲長樹, 沈桐立. 中國盛夏褐飛虱北遷過程的動態數值模擬. 昆蟲學報, 2000, 43(2): 176- 183.

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

[8] 封傳紅, 翟保平, 張孝羲, 湯金儀. 我國低空急流的時空分布與稻飛虱北遷. 生態學報, 2002, 22(4): 559- 565.

[9] 包云軒, 謝杰, 向勇, 徐希燕, 劉宇. 低空急流對中國褐飛虱重大北遷過程的影響. 生態學報, 2009, 29(11): 5773- 5782.

[10] 沈慧梅, 陳曉, 胡高, 程遐年, 張孝羲, 翟保平. 2008年廣西北部灣稻區稻飛虱初遷入過程分析. 應用昆蟲學報, 2011, 48(5): 1268- 1277.

[11] 吳秋琳, 胡高, 陸明紅, 王標, 朱秀秀, 粟芳, 張仲剛, 翟保平. 湖南白背飛虱前期遷入種群中小尺度蟲源地及降落機制. 生態學報, 2015, 35(22): 7397- 7417.

[12] Turner R, Song Y H, Uhm K B. Numerical model simulations of brown planthopperNilaparvatalugensand white-backed planthopperSogatellafurcifera(Hemiptera: Delphacidae) migration. Bulletin of Entomological Research, 1999, 89(6): 557- 568.

[13] Zhu M, Song Y H, Uhm K B, Turner R W, Lee J H, Roderick G K. Simulation of the long range migration of brown planthopper,Nilaparvatalugens(St?l), by using boundary layer atmospheric model and the geographic information system. Journal of Asia-Pacific Entomology, 2000, 3(1): 25- 32.

[14] Furuno A, Chino M, Otuka A, Watanabe T, Matsumura M, Suzuki Y. Development of a numerical simulation model for long-range migration of rice planthoppers. Agricultural and Forest Meteorology, 2005, 133(1/4): 197- 209.

[15] Cheng J, Norton G A, Holt J. A systems analysis approach to brown planthopper control on rice in Zhejiang Province, China. II. Investigation of control strategies. Journal of Applied Ecology, 1990, 27(1):100- 112.

[16] Crummay F A, Atkinson B W. Atmospheric influences on light-trap catches of the brown planthopper rice pest. Agricultural and Forest Meteorology, 1997, 88(1/4): 181- 197.

[17] 江廣恒, 談涵秋, 沈婉貞, 程遐年, 陳若篪. 褐飛虱遠距離向南遷飛的氣象條件. 昆蟲學報, 1982, 25(2): 147- 155.

[18] 程遐年, 張孝羲, 程極益. 褐飛虱在中國東部秋季回遷的雷達觀察. 南京農業大學學報, 1994, 17(3): 24- 32.

[19] 翟保平. 追蹤天使——雷達昆蟲學30年. 昆蟲學報, 1999, 42(3): 315- 326.

[20] Pender J. Migration of the brown planthopper,Nilaparvatalugens(Stal.) with special reference to synoptic meteorology. Grana, 1994, 33(2): 112- 115.

[21] Stohl A, Forster C, Frank A, Seibert P, Wotawa G. Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2. Atmospheric Chemistry and Physics, 2005, 5: 2461- 2474.

[22] 陳若篪, 程遐年. 褐飛虱起飛行為與自身生物學節律、環境因素同步關系的初步研究. 南京農學院學報. 1980, (2): 42- 49.

[23] 湯國安, 楊昕. ArcGIS地理信息系統空間分析實驗教程 (第二版). 北京: 科學出版社, 2012.

[24] Vincenty T. Direct and inverse solutions of geodesics on the ellipsoid with application of nested equations. Survey Review, 1975, 23: 88- 93.

[25] Vincenty T. Geodetic Inverse Solution Between Antipodal Points (Technical Report). DMAAC Geodetic Survey Squadron, 1975.

[26] 程極益, 包云軒, 樊多琦, 王海扣, 程遐年. 褐飛虱遷飛軌跡研究. 南京農業大學學報, 1995, 18(3): 60- 67.

[27] 胡高, 包云軒, 王建強, 翟保平. 褐飛虱的降落機制. 生態學報, 2007, 27(12): 5068- 5075.

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

A modeling study ofNilaparvatalugens′ return migration by using WRF-Flexpart

LIU Yao1,2,3, BAO Yunxuan1,2,*, LU Minghong4, LIU Wancai4

1 Key Laboratory of Agricultural Meteorology in Jiangsu Province, Nanjing University of Information Science and Technology, Nanjing 210044, China2 Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044, China3 Key Laboratory of Characteristic Agrometeorological Disaster Monitoring and Early Warning and Risk Management in Arid Regions, CMA, Yinchuan 750002, China4 National Agricultural Technology Extension and Service Center, Ministry of Agricultural, Beijing 100125, China

Nilaparvatalugens(N.Lugens) is one of the pests which have a large detrimental impact on rice production in Asia. Long-distance migration and population distribution are highly dependent on large-scale atmospheric circulation and local weather condition. The observational data of light-trapped pests from 1979 to 2015 were analyzed to quantify the impact of large-scale atmospheric circulation on southward migration ofN.Lugens. A typical heavy event of southward migration occurring in the major rice-growing regions of the southern China from September 30th to October 7th, 2009 was presented to discuss the impact of different weather conditions on immigration population ofN.Lugensin different rice-growing areas in South or Southeast China. The event was characterized by large-range precipitation in South China. A large amount of precipitation and high temperature were observed in the coastal regions of Fujian Province, South Jiangxi Province, and Southwest Hunan Province. In contrast, the rainfall was much less and temperature was relatively lower in the middle and north parts of Jiangxi Province compared to other regions and the rainfall was sporadic and the temperature was high in Hubei Province. As a result, the amount of light-trappedN.Lugenswas the largest in the coastal regions of Fujian Province, South Jiangxi Province, and Southwest Hunan Province, and the least in the middle and north parts of Jiangxi Province. The trapped amount in Hubei Province lied in between the two previous group regions. In this study, an online-coupling system of the Weather Research and Forecast Model (WRF) with FLEXible PARTicle dispersion model (FLEXPART) was utilized to simulate the meteorological conditions and to calculate the trajectories and the parameters associated withN.lugens′ migration. ArcGIS (a geographic information system software) was used to perform detailed spatial analysis and the Vincenty′s formula was applied to the calculation of azimuth and distance of migration. Several important results were identified from the study. First, meteorological factors including wind fields, rainfall, temperature, humidity conditions played vital and different roles in the immigration and land-falling ofN.lugensin different areas during different periods and contributed the discrepancy of temporal and spatial distributions ofN.lugens′ migration population. Second,N.lugenstended to immigrate toward the warm and humid areas. Migrating direction and migrating distance were determined by wind direction and wind speed of prevailing winds, respectively. Third, the immigrating amount ofN.lugenswas small at the beginning but can be accumulated to a significant level in a short time period. Fourth, within the certain range of air temperature, air humidity is the key parameters in determining migration population and temporal-spatial distribution. ImmigratingN.lugenswere concentrated in the regions with 2 m temperature higher than 26℃ and the relative humidity higher than 70% at the 925 hPa level whereas very few immigratingN.lugenscan be observed in the region with 2 m temperature lower than 24℃ and the relative humidity less than 50% at 925 hPa level.

Nilaparvatalugens(St?l); weather condition; WRF-FLEXPART model; migration trajectory; distribution of immigrating heads

國家自然科學基金(41475106,41075086)；江蘇省農業科技自主創新基金項目(CX(12)3056)；江蘇省高校自然科學研究項目(14KJA170003)；江蘇省普通高校研究生科研創新計劃項目(CXLX13_496)；江蘇高校優勢學科建設工程資助項目(IRT1147)

2015- 12- 31; 網絡出版日期:2017- 02- 23

10.5846/stxb201512312615

*通訊作者Corresponding author.E-mail: baoyx@nuist.edu.cn

劉垚,包云軒,陸明紅,劉萬才.基于WRF-Flexpart的一次褐飛虱回遷過程模擬研究.生態學報,2017,37(13):4466- 4475.

Liu Y, Bao Y X, Lu M H, Liu W C.A modeling study ofNilaparvatalugens′ return migration by using WRF-Flexpart.Acta Ecologica Sinica,2017,37(13):4466- 4475.