包云軒,蒼薪竹,楊詩(shī)俊,陳 粲,謝曉金,陸明紅,劉萬(wàn)才

1 氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心/南京信息工程大學(xué), 南京 210044 2 江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室/南京信息工程大學(xué), 南京 210044 3 全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心/農(nóng)業(yè)農(nóng)村部, 北京 100125

稻飛虱是一類(lèi)危害極重的專(zhuān)食性水稻害蟲(chóng),褐飛虱(Nilaparvatalugens(St?l))是三種主要稻飛虱中產(chǎn)生危害最嚴(yán)重的一種,它們可借助高空氣流實(shí)現(xiàn)種群大規(guī)模遠(yuǎn)距離的遷飛,由此產(chǎn)生毀滅級(jí)的災(zāi)害[1- 2]。自20世紀(jì)60年代初期,東亞和東南亞熱帶和亞熱帶地區(qū)實(shí)行“綠色革命”以來(lái),在感蟲(chóng)性水稻品種種植面積快速擴(kuò)增、產(chǎn)量快速增長(zhǎng)的同時(shí),于20世紀(jì)60年代后期至90年代初出現(xiàn)了第一輪蟲(chóng)害爆發(fā)周期。之后這些地區(qū)濫用殺蟲(chóng)劑、過(guò)量施用氮肥等行為導(dǎo)致褐飛虱于21世紀(jì)00年代中期至今進(jìn)入第二輪大爆發(fā)[3]。常年從東南亞地區(qū)頻繁且大量遷入到中國(guó)境內(nèi)的褐飛虱嚴(yán)重地影響了中國(guó)水稻的品質(zhì)和產(chǎn)量[4-5],由此造成中國(guó)水稻年均減產(chǎn)可達(dá)數(shù)億噸[6],嚴(yán)重威脅了中國(guó)的糧食安全。

褐飛虱年內(nèi)初期遷入我國(guó)的情況往往可以決定整年蟲(chóng)害的發(fā)生程度[7]。華南稻區(qū)和西南稻區(qū)作為每年年初境外褐飛虱種群遷入我國(guó)的第一站,其褐飛虱的春季發(fā)生量與全國(guó)發(fā)生面積、發(fā)生程度有著至關(guān)重要的聯(lián)系[8],遷入期早、遷入峰期提前也加大了我國(guó)南方稻區(qū)早稻上褐飛虱種群爆發(fā)的風(fēng)險(xiǎn)[9]。

大氣溫度對(duì)褐飛虱遷入期的早遲有著十分重要的影響,這是因?yàn)榇髿鉁囟瓤梢韵拗坪诛w虱完整的遷飛行為[10]。例如,近地面溫度在13℃以下時(shí)大部分褐飛虱不會(huì)起飛,起飛后的褐飛虱種群也多集中在溫度范圍在16—18℃的對(duì)流層大氣中飛行[11- 12],遇低溫則出現(xiàn)集中降落的行為等。其中,促使褐飛虱發(fā)生降落行為的低溫屏障往往是南北方的冷暖氣團(tuán)交匯形成的低溫墻[13]。穩(wěn)定的低溫層結(jié)可迫使褐飛虱逐漸接近地面或向上遇冷停止振翅從而導(dǎo)致其直接降落。有研究表明,全球變暖引發(fā)的溫度普遍升高間接使褐飛虱年內(nèi)初期遷入我國(guó)的時(shí)間提前,其中一部分原因是蟲(chóng)源地溫度的升高更利于褐飛虱種群的遷出[14];而另一部分原因則是我國(guó)西南、華南兩稻區(qū)地處熱帶與亞熱帶的過(guò)渡地帶,初春時(shí)節(jié)氣溫變化復(fù)雜,極易產(chǎn)生低溫層結(jié),而溫度的升高極有可能弱化了低溫屏障的限制作用。因此,研究大氣低溫脅迫對(duì)境外褐飛虱年內(nèi)初始遷入中國(guó)的影響對(duì)我國(guó)褐飛虱種群早期入侵的預(yù)警具有重要的理論價(jià)值和實(shí)踐意義。但目前相關(guān)的研究多以春季遷入個(gè)例為研究對(duì)象[13, 15],幾乎沒(méi)有對(duì)多年褐飛虱初始遷入我國(guó)過(guò)程中的低溫脅迫發(fā)生特征進(jìn)行過(guò)系統(tǒng)性、規(guī)律性的總結(jié)與討論。

在現(xiàn)今全球變暖愈演愈烈、極端氣候事件頻發(fā)的大背景下,氣候異常增加了溫度對(duì)褐飛虱遷飛行為影響的不確定性,同時(shí),大氣溫度對(duì)褐飛虱遷飛過(guò)程的影響程度和作用機(jī)理也趨于復(fù)雜化。因此,系統(tǒng)性地了解歷史上褐飛虱年內(nèi)初始遷入我國(guó)時(shí)的低溫脅迫發(fā)生規(guī)律,明確當(dāng)前發(fā)生特征,對(duì)預(yù)測(cè)未來(lái)氣候變化情景下可能的發(fā)生趨勢(shì)顯得十分重要和迫切。本文擬利用中國(guó)多年褐飛虱遷入的蟲(chóng)情數(shù)據(jù),結(jié)合HYSPLIT軌跡模式,對(duì)2000—2017年境外褐飛虱年內(nèi)初始遷入中國(guó)華南和西南稻區(qū)的發(fā)生特征進(jìn)行統(tǒng)計(jì)分析和遷入過(guò)程模擬,以探明大氣低溫對(duì)褐飛虱遷入的影響,揭示低溫脅迫作用的時(shí)空差異。研究將有助于監(jiān)測(cè)境外褐飛虱初始遷入我國(guó)的動(dòng)態(tài),為全球變暖趨勢(shì)下褐飛虱遷入發(fā)生的準(zhǔn)確預(yù)測(cè)奠定基礎(chǔ)。

1 材料與方法

1.1 資料來(lái)源

(1)蟲(chóng)情資料：稻飛虱蟲(chóng)情資料由全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心提供,文中主要使用了云南省、廣西壯族自治區(qū)、廣東省2000—2017年褐飛虱逐日燈誘數(shù)據(jù)。

(2)氣象資料：本研究所用的氣象資料為美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(National Center for Environmental Prediction,NCEP)發(fā)布的全球6 h一次(UTC 00、06、12和18時(shí)),空間分辨率為1°×1°的FNL氣象再分析資料。

(3)其他資料：全球行政區(qū)劃底圖(1∶4000000),來(lái)源于中國(guó)國(guó)家基礎(chǔ)地理信息中心;30 m全球地表覆蓋數(shù)據(jù),來(lái)源于中國(guó)地理信息資源目錄服務(wù)系統(tǒng)。

1.2 研究方法

1.2.1調(diào)查方法

褐飛虱僅在我國(guó)華南南部和西南南部的少數(shù)地區(qū)存在越冬蟲(chóng)源,而在每年早春遷入我國(guó)南方稻區(qū)的蟲(chóng)源大多是從中南半島等地越境遷入[16]。因此,選取緊鄰中南半島的西南稻區(qū)(以云南省為代表)、華南稻區(qū)(以廣東省、廣西壯族自治區(qū)為代表)為研究區(qū),使用稻區(qū)內(nèi)各植保站點(diǎn)的每日燈誘蟲(chóng)情數(shù)據(jù)進(jìn)行大氣低溫脅迫對(duì)褐飛虱年內(nèi)初始遷入行為影響的研究(圖1)。

圖1 西南、華南稻區(qū)褐飛虱監(jiān)測(cè)植保站點(diǎn)分布Fig.1 Distribution of N. lugens′ monitoring plant protection stations in the rice-growing regions of the southwest China and the south China

綜合考慮年內(nèi)初始遷入的蟲(chóng)情發(fā)生情況,并參考2009年發(fā)布的《稻飛虱測(cè)報(bào)調(diào)查規(guī)范》對(duì)始見(jiàn)期及高峰日的定義,將研究區(qū)年內(nèi)褐飛虱首次單燈單日有蟲(chóng)的日期定為始見(jiàn)期,年內(nèi)單燈單日蟲(chóng)量首次突增超過(guò)100頭為首次遷入峰發(fā)生日,并以此統(tǒng)計(jì)出2000—2017年兩個(gè)稻區(qū)的始見(jiàn)期和首次遷入峰日期的情況。

1.2.2軌跡分析方法

軌跡回推分析采用的是美國(guó)國(guó)家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)與澳大利亞氣象局(Bureau of Meteorology,BOM)開(kāi)發(fā)的HYSPLIT(Hybrid Single Particle Lagrangian Integrated Trajectory)模型。參數(shù)設(shè)置如下：(1)褐飛虱順風(fēng)遷飛,其主動(dòng)飛行能力是可忽略項(xiàng)[17-18];(2)軌跡回推到的蟲(chóng)源附近曾存在理論上可提供蟲(chóng)源取食的稻區(qū);(3)褐飛虱在日出前1 h和日落后1 h內(nèi)起飛[19-20]。西南、華南兩稻區(qū)3—4月份日出時(shí)刻在BST 07:00左右,日落時(shí)刻在BST 19:00左右。春季早晚氣溫較低,褐飛虱一般集中在較溫暖的下午飛行,屬于“日間單峰起飛”型[20]。故設(shè)定燈下見(jiàn)蟲(chóng)日當(dāng)晚20:00、燈下見(jiàn)蟲(chóng)第二日清晨06:00兩個(gè)起始時(shí)刻,分別回推24 h、10 h至前一日晚20:00。(4)褐飛虱的盛行遷飛高度位于海拔1500 m左右[21],特別說(shuō)明：在模型中設(shè)置的起始高度為離地高度。西南稻區(qū)地勢(shì)高且水稻種植區(qū)在海拔幾百米至幾千米的范圍內(nèi)均有所分布[22],在離地高度相同的條件下,大氣溫度遠(yuǎn)低于同緯度地區(qū)。而褐飛虱無(wú)法在環(huán)境溫度過(guò)低的條件下飛行,因此,模型中設(shè)置的西南稻區(qū)的離地高度相對(duì)較低?？梢?jiàn),褐飛虱的遷飛高度因受到溫度、地形等條件限制,需要根據(jù)季節(jié)及站點(diǎn)海拔的不同進(jìn)行相應(yīng)的調(diào)整,模型中兩稻區(qū)遷飛軌跡起始高度設(shè)定如表1。

表1 基于HYSPLIT模型的西南、華南稻區(qū)遷飛軌跡起始高度設(shè)定

1.2.3數(shù)據(jù)統(tǒng)計(jì)方法

為了解年內(nèi)初始遷入過(guò)程中低溫脅迫發(fā)生的整體情況,文中使用Excel等軟件統(tǒng)計(jì)2000—2017年中低溫屏障發(fā)生頻率、低溫屏障強(qiáng)度、低溫屏障幅度、低溫屏障出現(xiàn)時(shí)間分布等特征,并作相關(guān)描述。各特征具體判定方法如下：(1)低溫屏障：假設(shè)褐飛虱所在飛行高度層環(huán)境溫度低于16℃時(shí),褐飛虱即停止飛行[11];(2)假設(shè)降蟲(chóng)區(qū)整晚均有降蟲(chóng)。針對(duì)不同高度層,若降蟲(chóng)區(qū)在降蟲(chóng)當(dāng)晚最低氣溫小于16℃,即判定存在低溫屏障。18 a中出現(xiàn)低溫屏障的次數(shù)除以總年數(shù)記為低溫屏障發(fā)生頻率。僅針對(duì)受低溫脅迫的降蟲(chóng)過(guò)程計(jì)算低溫強(qiáng)度和降溫幅度。(3)低溫強(qiáng)度：針對(duì)不同高度層,在降蟲(chóng)區(qū)降蟲(chóng)當(dāng)晚20:00、次日02:00、次日08:00(記蟲(chóng))三個(gè)時(shí)刻中出現(xiàn)的最低溫度記為低溫強(qiáng)度;降溫幅度：將盛行遷飛高度層(西南稻區(qū)800 hPa,華南稻區(qū)850 hPa)上降蟲(chóng)區(qū)降蟲(chóng)當(dāng)晚的最大溫差記為降溫幅度。

2 結(jié)果與分析

2.1 初始遷入時(shí)期蟲(chóng)情統(tǒng)計(jì)結(jié)果

2.1.1始見(jiàn)期統(tǒng)計(jì)結(jié)果

2000—2017年間,褐飛虱在我國(guó)西南稻區(qū)最早于3月10日始見(jiàn),華南稻區(qū)最早于3月1日始見(jiàn)(表2)。在西南稻區(qū)內(nèi)所有站點(diǎn)的平均始見(jiàn)期為4月24日±2.15日(平均始見(jiàn)期平均值±標(biāo)準(zhǔn)誤差),華南稻區(qū)所有站點(diǎn)的平均始見(jiàn)期為4月10日±2.72日。在18 a間,褐飛虱始見(jiàn)于華南稻區(qū)的日期較同緯度的西南稻區(qū)整體早9—14 d。2001—2008年西南稻區(qū)內(nèi)最早始見(jiàn)期出現(xiàn)在4月10日前后,2009—2014年提前至4月初,到2015—2017年則穩(wěn)定在3月10日。2000—2003年華南稻區(qū)內(nèi)最早始見(jiàn)期在4月初,2004—2009年提前至3月上旬,到2010—2017年則穩(wěn)定在3月1日。如圖2中所示,在同一稻區(qū)的始見(jiàn)期,盡管年間存在波動(dòng),整體上仍呈顯著提前的趨勢(shì)(始見(jiàn)期與年份呈顯著負(fù)相關(guān),L1：r1=-0.803,P<0.01;L2：r2=-0.801,P<0.01)。相關(guān)研究表明,這種提前趨勢(shì)可能與近30 a的全球氣候變暖有關(guān)[14]。

表2 2000—2017年西南、華南稻區(qū)的褐飛虱遷入始見(jiàn)期

圖2 2000—2017年西南、華南稻區(qū)褐飛虱始見(jiàn)期Fig.2 Annual immigration beginning dates of N. lugens in the rice-growing regions of the southwest China and the south China during 2000 to 2017

2.1.2首次遷入峰統(tǒng)計(jì)結(jié)果

2000—2017年間,境外褐飛虱遷入峰首次出現(xiàn)在西南稻區(qū)的日期最早為3月21日,華南稻區(qū)最早為3月6日(表3)。在西南稻區(qū)內(nèi)全部站點(diǎn)的平均首遷峰日為5月29日±4.82日,華南稻區(qū)全部站點(diǎn)的平均首遷峰日為6月12日±3.96日,結(jié)果表明春季褐飛虱首次遷入西南稻區(qū)站點(diǎn)的時(shí)間在平均水平上反而早于華南稻區(qū)。為了排除因兩稻區(qū)站點(diǎn)分布差異對(duì)結(jié)果分析造成的影響,選取兩稻區(qū)緯度在22°N—24°N內(nèi)的所有站點(diǎn),調(diào)查其平均首遷峰日,所得結(jié)果相同。這顯示在每年的4—7月褐飛虱種群開(kāi)始大量遷入我國(guó)境內(nèi)時(shí),若以百頭蟲(chóng)量為準(zhǔn),則褐飛虱遷入西南稻區(qū)的時(shí)間更早。其原因一方面,在北半球南亞和東南亞低緯度地區(qū),西南季風(fēng)爆發(fā)早于南海季風(fēng),同時(shí),在陸上越靠近印度低壓的地區(qū)西南季風(fēng)越強(qiáng),就會(huì)越早出現(xiàn)有利于褐飛虱北遷的風(fēng)場(chǎng)[23-24]。因此,從地理經(jīng)度和緯度上講,西南稻區(qū)在褐飛虱初始遷入風(fēng)場(chǎng)上比華南稻區(qū)更有利;另一方面,在每年4—5月不僅從緬甸、泰國(guó)、越南等蟲(chóng)源區(qū)到中國(guó)的西南、華南稻區(qū)的空中溫場(chǎng)條件都能適應(yīng)于褐飛虱的北遷,且此時(shí)正為緬甸旱季稻的黃熟期,為西南稻區(qū)提供了大量蟲(chóng)源[13]。這說(shuō)明相較于大氣溫度,此時(shí)蟲(chóng)源量大小與是否存在適宜攜蟲(chóng)的風(fēng)場(chǎng)已成為影響褐飛虱遷入時(shí)間的主要因素。

綜合整個(gè)初始遷入期的統(tǒng)計(jì)結(jié)果來(lái)看,褐飛虱初始遷入我國(guó)華南稻區(qū)的日期整體上要早于西南稻區(qū)14 d左右。如圖3中所示,首次遷入峰日有提前趨勢(shì)但未形成明顯的時(shí)間差異性,主要是年度間的差異(首次遷入峰日早遲與年份負(fù)相關(guān)但不顯著,L3：r3=-0.698,P<0.01;L4：r4=-0.426,P>0.05)。例如,華南稻區(qū)2000和2001年、與2016和2017年幾乎相同;西南稻區(qū)2001和2002年與2011和2012年也幾乎相同(表3)。與首遷峰日相比,始見(jiàn)期在不同年份間變動(dòng)的幅度更小,出現(xiàn)的時(shí)間更穩(wěn)定。

表3 2000—2017年西南、華南稻區(qū)的褐飛虱首次遷入峰日

圖3 2000—2017年西南、華南稻區(qū)褐飛虱首次遷入峰日Fig.3 First immigration peak dates of N. lugens in the rice-growing regions of the southwest China and the south China during 2000 to 2017

2.2 年內(nèi)初始遷入種群的蟲(chóng)源地分布及其遷飛軌跡模擬

為探明多年來(lái)年內(nèi)初始遷入我國(guó)褐飛虱種群的境外蟲(chóng)源分布,本文對(duì)18 a中遷入中國(guó)西南、華南稻區(qū)的初始遷入過(guò)程進(jìn)行了軌跡回推。按不同稻區(qū)、不同遷入期將有效回推軌跡疊加到全球地表稻田分布圖上(30 m分辨率)(圖4)。以軌跡走向及蟲(chóng)源地分布為依據(jù),將各年回推軌跡進(jìn)行分類(lèi)統(tǒng)計(jì)(表4)。

總結(jié)多年回推軌跡的蟲(chóng)源地,得出年內(nèi)初始遷入我國(guó)西南稻區(qū)的蟲(chóng)源地主要分布在緬甸北部和中東部,華南稻區(qū)年內(nèi)初始遷入蟲(chóng)源主要來(lái)自老撾、越南北部和中國(guó)的海南南部。在褐飛虱年內(nèi)初始遷入西南稻區(qū)的過(guò)程中,共約有87.5%遷入過(guò)程的蟲(chóng)源來(lái)自緬甸北部和中東部稻區(qū),其余少部分蟲(chóng)源來(lái)自于老撾、越南。而在初始遷入華南稻區(qū)的過(guò)程中,共約有66.7%遷入過(guò)程的蟲(chóng)源來(lái)自于老撾和越南的中部或北部稻區(qū),22.2%遷入過(guò)程的蟲(chóng)源來(lái)自于緬甸、老撾交界邊境以及海南島。其余少量未歸類(lèi)回推軌跡的蟲(chóng)源地分布在雷州半島和廣西南部稻區(qū),為區(qū)內(nèi)的短途遷飛。如圖4,在整個(gè)初始遷入過(guò)程中回推軌跡多為東—西走向或西南—東北走向,只有個(gè)別過(guò)程取偏南至偏北走向。此外,褐飛虱在首次遷入峰時(shí)期較始見(jiàn)期相比遷入地范圍更大,蟲(chóng)源地更廣,軌跡聚合度更低。

表4 2000—2017年西南、華南稻區(qū)初始蟲(chóng)源地及其軌跡分類(lèi)

圖4 2000—2017年西南、華南稻區(qū)褐飛虱初始蟲(chóng)源地及其軌跡分布Fig.4 Distribution of source areas and trajectories for N. lugens′ initial immigration in the rice-growing regions of the southwest China and the south China during 2000 to 2017

2.3 低溫脅迫統(tǒng)計(jì)特征

使用1.2.3數(shù)據(jù)統(tǒng)計(jì)方法中介紹的方法得到低溫脅迫統(tǒng)計(jì)特征值如下表(表5、表6、表7)。為簡(jiǎn)明揭示溫度對(duì)褐飛虱整個(gè)遷飛過(guò)程的復(fù)雜影響,下面僅對(duì)褐飛虱在盛行遷飛層上的脅迫情況進(jìn)行討論。統(tǒng)計(jì)所有受低溫脅迫影響的降落過(guò)程發(fā)現(xiàn),在降蟲(chóng)當(dāng)晚,降溫過(guò)程有21.9%出現(xiàn)在20:00—次日02:00之間,25%出現(xiàn)在02:00—08:00之間,52.5%整晚降溫。在2000—2017年褐飛虱年內(nèi)初始遷入西南稻區(qū)的過(guò)程中,約有50.0%的降落過(guò)程受到了低溫影響,低溫過(guò)程平均強(qiáng)度為13.71℃(表7);在初始遷入華南稻區(qū)的過(guò)程中,約有58.3%的降落過(guò)程受到了低溫影響,低溫過(guò)程平均強(qiáng)度為13.18℃?？梢?jiàn),褐飛虱降落在華南稻區(qū)受低溫脅迫影響相對(duì)更大。相比于華南稻區(qū),西南稻區(qū)因緯度跨距大、海拔高度起伏明顯而溫度變化幅度大,但因云南稻區(qū)地形北高南低,受地形阻擋,一般冷空氣難以侵入西南南部稻區(qū)(即云南中南部稻區(qū)),導(dǎo)致褐飛虱種群降落對(duì)低溫脅迫的響應(yīng)概率小;華南稻區(qū)因緯度跨距小、海拔高度差異小而溫度變化幅度小,南嶺山系雖呈東西帶狀分布,但其五個(gè)山嶺之間是間斷的,且海拔高度明顯低于云南北部山系,對(duì)冷空氣的阻力也明顯小于西南南部稻區(qū),導(dǎo)致華南稻區(qū)受冷空影響范圍廣、發(fā)生概率高,褐飛虱種群降落對(duì)低溫脅迫的響應(yīng)概率和強(qiáng)度都大。綜上,春季西南稻區(qū)日均溫高,溫度日較差大,褐飛虱種群降落對(duì)低溫脅迫的響應(yīng)小;華南稻區(qū)日均溫低,溫度日較差小,褐飛虱降落對(duì)低溫脅迫的響應(yīng)大。

不同遷入時(shí)期的降蟲(chóng)過(guò)程受低溫脅迫影響的差異在兩稻區(qū)表現(xiàn)不同,在華南稻區(qū)表現(xiàn)出明顯差異,在西南稻區(qū)表現(xiàn)為無(wú)差異。華南稻區(qū)褐飛虱始見(jiàn)期受到低溫脅迫影響的概率(88.9%)遠(yuǎn)高于西南稻區(qū)始見(jiàn)期(50%)。而在首次遷入峰期中,華南稻區(qū)受低溫脅迫影響的概率(27.8%)明顯低于西南稻區(qū)(50%)(表7)。對(duì)不同遷入時(shí)期的低溫脅迫發(fā)生情況統(tǒng)計(jì)得到,始見(jiàn)期約有70.6%的降落過(guò)程受低溫影響,低溫強(qiáng)度為12.53℃,而在首次遷入峰時(shí)期僅有38.2%的降落過(guò)程受到了低溫影響,低溫強(qiáng)度為14.36℃。統(tǒng)計(jì)發(fā)現(xiàn),一般情況下,華南稻區(qū)始見(jiàn)期較首次遷入峰日要早2—3個(gè)月,始見(jiàn)期日均溫低于16℃的頻率要遠(yuǎn)高于首次遷入峰時(shí)期,所以在始見(jiàn)期低溫脅迫對(duì)降蟲(chóng)過(guò)程產(chǎn)生的影響更顯著。但是,在西南稻區(qū)的統(tǒng)計(jì)卻顯示：這種低溫脅迫的時(shí)間差異性很小,一方面原因可能是北遷過(guò)程中地形阻擋因素會(huì)使部分褐飛虱遷入種群在遇到低溫之前就遇到了地形動(dòng)力脅迫而產(chǎn)生了降落,另一方面是因?yàn)樵颇洗杭練鉁刈兓刃?晝夜溫差大,故在時(shí)間維度上西南地區(qū)上空大氣溫度低于褐飛虱種群16℃飛行閾值的概率差異小,從而受到低溫脅迫影響降蟲(chóng)過(guò)程的時(shí)間差異減小。因此,不同遷入期的降蟲(chóng)過(guò)程受低溫脅迫影響的差異在西南稻區(qū)就很小。

綜上,褐飛虱初始遷入我國(guó)受低溫脅迫的降落過(guò)程在時(shí)間和空間分布上均存在差異,特別是在時(shí)間分布上的差異更加明顯。就始見(jiàn)期而言,華南稻區(qū)的降溫幅度小于西南稻區(qū),但低溫脅迫在華南稻區(qū)的影響頻次比西南稻區(qū)高、發(fā)生時(shí)間早、強(qiáng)度大;而在首次遷入峰上,低溫脅迫在西南稻區(qū)的影響頻次比華南稻區(qū)高、發(fā)生時(shí)間早、強(qiáng)度大、降溫幅度也大(表7)。

表5 2000—2017年西南稻區(qū)逐年800 hPa低溫脅迫特征值

表6 2000—2017年華南稻區(qū)逐年850 hPa低溫脅迫特征值

表7 2000—2017年西南、華南稻區(qū)不同高度層低溫屏障統(tǒng)計(jì)

2.4 遷入地溫度分析

為能更全面地了解褐飛虱遷入與低溫脅迫的關(guān)系,本文統(tǒng)計(jì)了在整個(gè)初始遷入時(shí)期降蟲(chóng)地盛行遷飛層的日均溫,遷入日用黑色橫粗線(xiàn)覆蓋。從圖5可見(jiàn),西南稻區(qū)整個(gè)初始遷入期日均溫低于16℃的頻次在整體趨勢(shì)上逐年減少,始見(jiàn)期、首次遷入峰逐年提前,這表明西南稻區(qū)褐飛虱提早遷入可能與降蟲(chóng)地盛行遷飛層的溫度升高有關(guān)。圖中顯示并非在溫度符合條件后(大于16℃溫度閾值)褐飛虱立即遷入降蟲(chóng)區(qū),而是遷入時(shí)的溫度基本都符合大于16℃的條件(如圖5中2010年)。這說(shuō)明溫度僅是限制攜蟲(chóng)氣流遷入某地的一個(gè)重要因素,而不是影響初始遷入時(shí)期早晚的決定性因素。

為進(jìn)一步了解始見(jiàn)期出現(xiàn)早遲與近年來(lái)降蟲(chóng)區(qū)盛行遷飛層溫度變化間的關(guān)系,我們通過(guò)調(diào)查華南、西南稻區(qū)遷入頻率最高的六個(gè)站點(diǎn)(西南稻區(qū)：景洪JH、金平JP、屏邊PB,華南稻區(qū)：雷州LZ、武鳴WM、防城FC)2000—2017年3月份盛行遷飛層的月均溫,首先分析18年來(lái)月均溫的年度間變化得出：近年來(lái)西南稻區(qū)的站點(diǎn)溫度有升高的趨勢(shì)(圖6),而華南稻區(qū)的站點(diǎn)溫度幾乎沒(méi)有變化(圖6)。其次,分析月均溫與始見(jiàn)期出現(xiàn)時(shí)間的相關(guān)性得出：西南稻區(qū)3月份的站點(diǎn)溫度與始見(jiàn)期早遲具有較弱的負(fù)相關(guān)性,即盛行遷飛層溫度越高,始見(jiàn)期越早。這說(shuō)明西南稻區(qū)褐飛虱始見(jiàn)期提前可能與其盛行遷飛層上大氣溫度升高有關(guān),而華南稻區(qū)的始見(jiàn)期提前與其盛行遷飛層上的大氣溫度變化無(wú)關(guān)。

在華南稻區(qū)初始遷入期,特別是始見(jiàn)期(圖5)的遷入地盛行遷飛層溫度低于16℃的頻次逐年增多,初始遷入期卻仍有提前,甚至83.3%遷入日的日均溫都低于16℃,個(gè)別年份遠(yuǎn)低于16℃(最低溫9.75℃)。但眾所周知,褐飛虱在低于16℃的環(huán)境下是無(wú)法繼續(xù)飛行的。由此推測(cè),華南稻區(qū)始見(jiàn)期蟲(chóng)源的飛行高度一定低于1500 m(850 hPa),同時(shí)回推蟲(chóng)源地與降蟲(chóng)地的距離相對(duì)縮短,故可能為越南北部稻區(qū)或海南、華南南部稻區(qū)存在的少量越冬蟲(chóng)源隨偏南風(fēng)短途低空遷飛后降落本站的,而境內(nèi)越冬蟲(chóng)源的活動(dòng)期提前導(dǎo)致華南稻區(qū)褐飛虱燈下始見(jiàn)期提前與近年來(lái)氣候變暖使褐飛虱越冬蟲(chóng)源分布北界北移有關(guān)。

圖5 2000—2017年西南、華南稻區(qū)初始遷入期降蟲(chóng)地盛行遷飛高度層(西南稻區(qū)800 hPa/華南稻區(qū)850 hPa)日均溫(℃)Fig.5 Daily mean temperature on the prevailing migratory layers (on 800 hPa in the southwest China or 850 hPa in the south China) of N. lugens′ landing areas during the initial immigration dates in the rice-growing regions of the southwest China and the south China during 2000 to 2017

圖6 2000—2017年西南稻區(qū)站點(diǎn)(JH景洪、JP金平、PB屏邊),華南稻區(qū)站點(diǎn)(LZ雷州、WM武鳴、FC防城)盛行遷飛高度層(西南稻區(qū)800 hPa/華南稻區(qū)850 hPa)3月份月均溫(℃)與始見(jiàn)期Fig.6 Monthly mean temperature in March of the prevalent migration layers (on 800 hPa in the southwest China or 850 hPa in the south China) of the stations in the rice-growing regions of the southwest China (JH, JP, and PB) and the south China (LZ, WM, and FC), and N. lugens′ immigration beginning dates during 2000 to 2017

3 結(jié)論與討論

3.1 褐飛虱初始遷入我國(guó)時(shí)間提前,華南稻區(qū)始見(jiàn)期提早可能與越冬北界的北移有關(guān)

近年來(lái),我國(guó)褐飛虱年內(nèi)初始遷入期有所提前。有研究發(fā)現(xiàn)在氣候變暖背景下,廣西、云南等地稻飛虱的遷入始見(jiàn)期、首次遷入高峰期均呈逐年提早的趨勢(shì)[25-26, 9]。稻飛虱遷入始見(jiàn)期的出現(xiàn)早遲在我國(guó)東、西部間存在著差異,在相同緯度上我國(guó)東部始見(jiàn)期要早于西部地區(qū)[26]。本文統(tǒng)計(jì)多年褐飛虱遷入蟲(chóng)情得到,褐飛虱初始遷入華南稻區(qū)時(shí)間較西南稻區(qū)要早9—14 d。褐飛虱初始遷入我國(guó)東、西部稻區(qū)上的時(shí)間差異不僅與降蟲(chóng)區(qū)的地形、溫濕度條件、水稻種植等情況有關(guān),蟲(chóng)源地的蟲(chóng)源情況以及距離降蟲(chóng)地的遠(yuǎn)近也與之密切相關(guān)。華南稻區(qū)西連越南,南臨南海,盡管地貌類(lèi)型多樣,但整體海拔較低,與老撾、越南等蟲(chóng)源地距離近且無(wú)阻擋,能為兩廣稻區(qū)常年提供穩(wěn)定蟲(chóng)源,因此,褐飛虱年內(nèi)初始遷入期更早。而西南地區(qū)地處青藏高原東南側(cè),不僅受到東亞季風(fēng)和印度季風(fēng)及中高緯度天氣系統(tǒng)的交叉影響[27],而且西南地區(qū)海拔高、地形地勢(shì)多變,褐飛虱遷入受地理因素阻礙較大,褐飛虱遷入的降蟲(chóng)機(jī)制相比東部其他省份要更加復(fù)雜,多種因素綜合導(dǎo)致褐飛虱年內(nèi)遷入我國(guó)西部始見(jiàn)期要晚于東部地區(qū)。

本研究中發(fā)現(xiàn),華南稻區(qū)在盛行遷飛層上的大氣溫度無(wú)升高的趨勢(shì),但初始遷入期仍然提前了。特別是在華南稻區(qū)始見(jiàn)期,約有83.3%褐飛虱遷入過(guò)程的降蟲(chóng)區(qū)最適遷飛層的大氣溫度低于16℃,在這種情況下褐飛虱是無(wú)法飛行的,但實(shí)際情況卻是有蟲(chóng)遷入。這說(shuō)明褐飛虱在3月份遷入華南稻區(qū)時(shí)的遷飛高度遠(yuǎn)低于1500 m,同時(shí),近年來(lái)研究表明氣候變暖使褐飛虱越冬北界北移,在華南南部稻區(qū)存在少量的越冬蟲(chóng)源[28],在華南北部稻區(qū)則存在低空短途遷飛的現(xiàn)象[29]。Hu等[30]在研究2008年廣西北部灣稻區(qū)初始遷入的過(guò)程中模擬500 m高度下的蟲(chóng)源,結(jié)果顯示低空短途遷入華南的蟲(chóng)源可能來(lái)自于海南島。上述情況可以說(shuō)明存在少量越冬蟲(chóng)源在年內(nèi)早期低空短途遷入華南稻區(qū)的可能性。它的遷飛機(jī)制為褐飛虱主動(dòng)起飛升空后,一般如無(wú)直達(dá)16℃低溫層的強(qiáng)上升氣流,就會(huì)逐漸主動(dòng)適應(yīng)環(huán)境溫度在適溫層內(nèi)飛行;但若遇到直達(dá)16℃以下低溫層的強(qiáng)上升氣流或在隨氣流水平飛行時(shí)突遇高空大氣低溫層(即低于16℃的水平氣流層)時(shí),如其不能及時(shí)反應(yīng)并主動(dòng)調(diào)節(jié)翅振頻率和飛行高度,并快速進(jìn)入適溫區(qū),則會(huì)致死或驟停振翅,掉下地面。褐飛虱屬小型風(fēng)載遷飛性昆蟲(chóng),其在高空強(qiáng)氣流影響下對(duì)適宜溫度層主動(dòng)選擇的能力有限,而對(duì)流層卻有“邊界層頂效應(yīng)”[31],因此,在空中飛行時(shí),三維空間的氣流速度絕大多數(shù)情況下是遠(yuǎn)超其主動(dòng)遷飛的速度的,褐飛虱在這樣的環(huán)境流場(chǎng)和溫場(chǎng)條件下是被動(dòng)的。由此可以推測(cè),這可能是分布在越南北部稻區(qū)或海南、華南南部稻區(qū)的越冬蟲(chóng)源在遠(yuǎn)低于1500 m的飛行高度下短途遷入華南中部或北部稻區(qū)后而導(dǎo)致其他地區(qū)測(cè)報(bào)站在燈下見(jiàn)蟲(chóng)了。然而,盡管僅有10%左右的概率,仍不排除有蟲(chóng)源自越南中部跨海遷入的可能性[32]。

3.2 緬甸是西南稻區(qū)多年稻飛虱初始遷入過(guò)程的穩(wěn)定蟲(chóng)源區(qū),越南、老撾、海南是華南稻區(qū)多年稻飛虱初始遷入過(guò)程的穩(wěn)定蟲(chóng)源區(qū)。

云南、兩廣是地處我國(guó)南端邊陲的省份,也是稻飛虱遷入我國(guó)的第一站。關(guān)于云南以及兩廣地區(qū)的蟲(chóng)源問(wèn)題,許多科研工作者已經(jīng)做出了較為完整的研究。如沈慧梅[15]對(duì)2009年云南百背飛虱的蟲(chóng)源地分析得出云南早期蟲(chóng)源主要來(lái)自于緬甸和金三角地區(qū),而我國(guó)華南稻區(qū)的早期遷入蟲(chóng)源則主要來(lái)自越南中部。包云軒[13]等對(duì)云南省4、5月份早期褐飛虱遷入蟲(chóng)源地研究發(fā)現(xiàn),褐飛虱遷入的主要蟲(chóng)源區(qū)是緬甸,部分來(lái)自老撾和泰國(guó),少量來(lái)自越南,還有極少數(shù)來(lái)自孟加拉國(guó)。Hu等[30]等對(duì)2008年廣西早期遷入稻飛虱的蟲(chóng)源進(jìn)行分析后得出：廣西早期遷入蟲(chóng)源主要來(lái)自越南中部及老撾南部。王政等[33]對(duì)廣東白背飛虱多年早期遷入的蟲(chóng)源地進(jìn)行分析得出：蟲(chóng)源主要分布在越南中北部、老撾中南部、海南以及兩廣南部稻區(qū)。

本文對(duì)西南、華南兩稻區(qū)多年褐飛虱遷入始見(jiàn)期、首次遷入峰蟲(chóng)源地分布進(jìn)行逆推后發(fā)現(xiàn)：多年褐飛虱初始遷入過(guò)程的蟲(chóng)源地分布總是穩(wěn)定在固定區(qū)域內(nèi)。2000—2017年西南稻區(qū)褐飛虱始見(jiàn)期蟲(chóng)源均來(lái)自于緬甸,首次遷入高峰期蟲(chóng)源區(qū)大部分在緬甸,少量分布在老撾和越南。這是由于近年來(lái)緬甸政府在農(nóng)業(yè)政策上作了的改革,增加了水稻的種植面積和產(chǎn)量。在每年的4月上中旬,緬甸中北部稻區(qū)所種植的水稻處于抽穗、灌漿期,5月初水稻收獲。同時(shí)4月份緬甸中北部的高空存在明顯的盛行西風(fēng),形成了褐飛虱蟲(chóng)源遠(yuǎn)距離遷飛的輸送通道,為我國(guó)云南省褐飛虱的早期遷入提供了蟲(chóng)源[34]。而華南稻區(qū)褐飛虱始見(jiàn)期蟲(chóng)源主要分布在越南中北部及老撾北部,少量蟲(chóng)源來(lái)自緬甸和海南。越南中北部的紅河三角洲平原以盛產(chǎn)水稻著稱(chēng),冬春稻在每年2月份播種,同年4月份水稻抽穗楊花,此時(shí)有大量褐飛虱種群開(kāi)始外遷。同樣由于老撾、越南等地在3、4月份盛行的西南風(fēng)、西南偏南風(fēng)為褐飛虱種群遷入我國(guó)華南稻區(qū)提供了可能[32]。由此可見(jiàn),華南稻區(qū)蟲(chóng)源地分布更為復(fù)雜。

3.3 始見(jiàn)期的低溫脅迫在華南稻區(qū)影響頻次比西南稻區(qū)高、發(fā)生時(shí)間早、強(qiáng)度大;首次遷入峰的低溫脅迫在西南稻區(qū)的影響頻次比華南稻區(qū)高、發(fā)生時(shí)間早、強(qiáng)度大、降溫幅度也大。

褐飛虱完整的遷飛過(guò)程并非像模型模擬的那樣,其中還存在著褐飛虱主動(dòng)選擇的行為。遷飛的褐飛虱個(gè)體在遇到降雨、下沉氣流、風(fēng)向幅合等外力作用時(shí)會(huì)集中迫降,或在飛行能源物質(zhì)耗盡或溫度驟變低于其飛行溫度閾值時(shí)也會(huì)選擇主動(dòng)降落[35]。在不考慮食料條件的影響下,一次完整的降落過(guò)程是降蟲(chóng)地地形及多種氣象因素交互作用的結(jié)果。盡管降水、下沉氣流是影響褐飛虱大規(guī)模降落的兩大主要?dú)庀笠蜃覽33, 36],但在初始遷入的過(guò)程中,時(shí)值早春,大氣溫度低,降雨少,低溫對(duì)稻飛虱降落的影響更加明顯[7]。如沈慧梅[15]在對(duì)2009年云南白背飛虱降落機(jī)制分析中指出,白背飛虱早期遷入云南時(shí),風(fēng)切變與低溫屏障幾乎存在于每一次降落過(guò)程中。本文重點(diǎn)研究了低溫對(duì)褐飛虱初始遷入我國(guó)的脅迫作用,且初始遷入蟲(chóng)量較少,故不考慮降水以及下沉氣流對(duì)褐飛虱降落的影響。徐蘭珍等[37]在研究滇東南白背飛虱早期遷入過(guò)程中得到,影響白背飛虱的降落因子中低溫屏障占56.5%,與本文中得到的西南稻區(qū)早期遷入過(guò)程中影響褐飛虱的降落因子中低溫屏障占50%相近。本文通過(guò)統(tǒng)計(jì)兩稻區(qū)18 a中共68次初始遷入過(guò)程發(fā)現(xiàn),約有54.4%的遷入過(guò)程受到了低溫影響,且在不同稻區(qū)、不同時(shí)期(始見(jiàn)期、首次遷入峰期)的發(fā)生情況均有差異。正是因?yàn)樵诳臻g尺度上不同地理位置獨(dú)有的地形特征與氣候特點(diǎn)和時(shí)間尺度上氣候變化的特點(diǎn)綜合作用,使華南稻區(qū)始見(jiàn)期的遷入過(guò)程受低溫影響更大。

盡管本文對(duì)有關(guān)褐飛虱年內(nèi)初始遷入過(guò)程中低溫脅迫的特征進(jìn)行了相關(guān)的統(tǒng)計(jì)分析,但研究中還存在一些不足,如兩稻區(qū)2005—2017年為每日燈誘蟲(chóng)情數(shù)據(jù),2000—2004年為每候燈誘蟲(chóng)情數(shù)據(jù)且開(kāi)燈日期較晚(開(kāi)燈日期：2000- 04- 06,2001- 03- 06,2002- 03- 10,2003- 04- 06,2004- 03- 06,2005—2017年均為3月1日),同時(shí)西南稻區(qū)缺少2000年和2005年的蟲(chóng)情數(shù)據(jù)。這些部分蟲(chóng)情數(shù)據(jù)不完整的問(wèn)題會(huì)對(duì)本文研究結(jié)果產(chǎn)生一定的影響,且本文僅對(duì)低溫屏障及其特征進(jìn)行了初步的定量統(tǒng)計(jì)分析,而沒(méi)有涉及到低溫是如何在飛行過(guò)程中脅迫褐飛虱產(chǎn)生降落行為的以及種群對(duì)低溫脅迫的響應(yīng)是否會(huì)在不同時(shí)空條件下存在差異,這些機(jī)理性問(wèn)題需要結(jié)合個(gè)例進(jìn)行更加詳細(xì)深入的闡釋與剖析。