包云軒,蒼薪竹,楊詩俊,陳 粲,謝曉金,陸明紅,劉萬才

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

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

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

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

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

1 材料與方法

1.1 資料來源

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

(2)氣象資料：本研究所用的氣象資料為美國國家環(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),來源于中國國家基礎(chǔ)地理信息中心;30 m全球地表覆蓋數(shù)據(jù),來源于中國地理信息資源目錄服務(wù)系統(tǒng)。

1.2 研究方法

1.2.1調(diào)查方法

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

圖1 西南、華南稻區(qū)褐飛虱監(jiān)測植保站點(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)初始遷入的蟲情發(fā)生情況,并參考2009年發(fā)布的《稻飛虱測報(bào)調(diào)查規(guī)范》對始見期及高峰日的定義,將研究區(qū)年內(nèi)褐飛虱首次單燈單日有蟲的日期定為始見期,年內(nèi)單燈單日蟲量首次突增超過100頭為首次遷入峰發(fā)生日,并以此統(tǒng)計(jì)出2000—2017年兩個(gè)稻區(qū)的始見期和首次遷入峰日期的情況。

1.2.2軌跡分析方法

軌跡回推分析采用的是美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)與澳大利亞氣象局(Bureau of Meteorology,BOM)開發(fā)的HYSPLIT(Hybrid Single Particle Lagrangian Integrated Trajectory)模型。參數(shù)設(shè)置如下：(1)褐飛虱順風(fēng)遷飛,其主動(dòng)飛行能力是可忽略項(xiàng)[17-18];(2)軌跡回推到的蟲源附近曾存在理論上可提供蟲源取食的稻區(qū);(3)褐飛虱在日出前1 h和日落后1 h內(nèi)起飛[19-20]。西南、華南兩稻區(qū)3—4月份日出時(shí)刻在BST 07:00左右,日落時(shí)刻在BST 19:00左右。春季早晚氣溫較低,褐飛虱一般集中在較溫暖的下午飛行,屬于“日間單峰起飛”型[20]。故設(shè)定燈下見蟲日當(dāng)晚20:00、燈下見蟲第二日清晨06:00兩個(gè)起始時(shí)刻,分別回推24 h、10 h至前一日晚20:00。(4)褐飛虱的盛行遷飛高度位于海拔1500 m左右[21],特別說明：在模型中設(shè)置的起始高度為離地高度。西南稻區(qū)地勢高且水稻種植區(qū)在海拔幾百米至幾千米的范圍內(nèi)均有所分布[22],在離地高度相同的條件下,大氣溫度遠(yuǎn)低于同緯度地區(qū)。而褐飛虱無法在環(huán)境溫度過低的條件下飛行,因此,模型中設(shè)置的西南稻區(qū)的離地高度相對較低?？梢?褐飛虱的遷飛高度因受到溫度、地形等條件限制,需要根據(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)初始遷入過程中低溫脅迫發(fā)生的整體情況,文中使用Excel等軟件統(tǒng)計(jì)2000—2017年中低溫屏障發(fā)生頻率、低溫屏障強(qiáng)度、低溫屏障幅度、低溫屏障出現(xiàn)時(shí)間分布等特征,并作相關(guān)描述。各特征具體判定方法如下：(1)低溫屏障：假設(shè)褐飛虱所在飛行高度層環(huán)境溫度低于16℃時(shí),褐飛虱即停止飛行[11];(2)假設(shè)降蟲區(qū)整晚均有降蟲。針對不同高度層,若降蟲區(qū)在降蟲當(dāng)晚最低氣溫小于16℃,即判定存在低溫屏障。18 a中出現(xiàn)低溫屏障的次數(shù)除以總年數(shù)記為低溫屏障發(fā)生頻率。僅針對受低溫脅迫的降蟲過程計(jì)算低溫強(qiáng)度和降溫幅度。(3)低溫強(qiáng)度：針對不同高度層,在降蟲區(qū)降蟲當(dāng)晚20:00、次日02:00、次日08:00(記蟲)三個(gè)時(shí)刻中出現(xiàn)的最低溫度記為低溫強(qiáng)度;降溫幅度：將盛行遷飛高度層(西南稻區(qū)800 hPa,華南稻區(qū)850 hPa)上降蟲區(qū)降蟲當(dāng)晚的最大溫差記為降溫幅度。

2 結(jié)果與分析

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

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

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

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

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

綜合整個(gè)初始遷入期的統(tǒng)計(jì)結(jié)果來看,褐飛虱初始遷入我國華南稻區(qū)的日期整體上要早于西南稻區(qū)14 d左右。如圖3中所示,首次遷入峰日有提前趨勢但未形成明顯的時(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)。與首遷峰日相比,始見期在不同年份間變動(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)初始遷入種群的蟲源地分布及其遷飛軌跡模擬

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

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

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

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

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

綜上,褐飛虱初始遷入我國受低溫脅迫的降落過程在時(shí)間和空間分布上均存在差異,特別是在時(shí)間分布上的差異更加明顯。就始見期而言,華南稻區(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í)期降蟲地盛行遷飛層的日均溫,遷入日用黑色橫粗線覆蓋。從圖5可見,西南稻區(qū)整個(gè)初始遷入期日均溫低于16℃的頻次在整體趨勢上逐年減少,始見期、首次遷入峰逐年提前,這表明西南稻區(qū)褐飛虱提早遷入可能與降蟲地盛行遷飛層的溫度升高有關(guān)。圖中顯示并非在溫度符合條件后(大于16℃溫度閾值)褐飛虱立即遷入降蟲區(qū),而是遷入時(shí)的溫度基本都符合大于16℃的條件(如圖5中2010年)。這說明溫度僅是限制攜蟲氣流遷入某地的一個(gè)重要因素,而不是影響初始遷入時(shí)期早晚的決定性因素。

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

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

圖5 2000—2017年西南、華南稻區(qū)初始遷入期降蟲地盛行遷飛高度層(西南稻區(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月份月均溫(℃)與始見期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 褐飛虱初始遷入我國時(shí)間提前,華南稻區(qū)始見期提早可能與越冬北界的北移有關(guān)

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

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

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

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

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

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

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

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