• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    中國(guó)稻縱卷葉螟發(fā)生特點(diǎn)及北遷的大氣背景

    2015-02-07 06:50:59包云軒謝曉金陸明紅王純枝劉萬(wàn)才
    生態(tài)學(xué)報(bào) 2015年11期
    關(guān)鍵詞:稻區(qū)卷葉螟水汽

    包云軒,曹 云,謝曉金,陸明紅,李 軒,王純枝,劉萬(wàn)才

    1 江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室,南京信息工程大學(xué),南京 210044 2 農(nóng)業(yè)部全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心,北京 100125 3 國(guó)家氣象中心,北京 100081

    中國(guó)稻縱卷葉螟發(fā)生特點(diǎn)及北遷的大氣背景

    包云軒1,*,曹 云1,謝曉金1,陸明紅2,李 軒3,王純枝3,劉萬(wàn)才2

    1 江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室,南京信息工程大學(xué),南京 210044 2 農(nóng)業(yè)部全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心,北京 100125 3 國(guó)家氣象中心,北京 100081

    稻縱卷葉螟因其發(fā)生面積大、暴發(fā)頻率高、致害程度重而成為我國(guó)水稻生產(chǎn)上的重要害蟲(chóng)。為了闡明大氣背景對(duì)其北遷重大過(guò)程的影響,利用2000—2012年中國(guó)稻縱卷葉螟燈誘數(shù)據(jù)分析了我國(guó)水稻主產(chǎn)區(qū)稻縱卷葉螟遷入的主要特征;選取2007年7月下旬的一次重大北遷過(guò)程作為典型個(gè)例,運(yùn)用中尺度天氣研究和預(yù)報(bào)模式WRF,結(jié)合NCEP氣象再分析資料,模擬了這一過(guò)程的大氣背景,推算了稻縱卷葉螟的遷飛軌跡,并分析了大氣背景對(duì)稻縱卷葉螟災(zāi)變性遷入的影響。研究表明:(1)這13年中,我國(guó)稻縱卷葉螟的遷入大多在3月初始見(jiàn),3—8月為北遷期,從南到北先后在華南、西南、江嶺、江淮稻區(qū)出現(xiàn)遷入峰;9—11月為南遷期,從北到南先后遷入江淮、江嶺、華南稻區(qū)并出現(xiàn)相應(yīng)的遷入峰,10月底至11月初為終見(jiàn)期。(2)北遷個(gè)例中各站逆推軌跡分析顯示:包括再遷飛在內(nèi)的不同時(shí)段蟲(chóng)源地基本上位于降蟲(chóng)區(qū)的西南方,遷飛高度都變化于550—850 m之間,但在沿海地區(qū)降落的稻縱卷葉螟遷飛軌跡是復(fù)雜的。(3)水平氣流是稻縱卷葉螟遠(yuǎn)距離北遷的主要運(yùn)載動(dòng)力,925 hPa上南方稻區(qū)一致的偏南氣流對(duì)稻縱卷葉螟北遷極為有利;三維流場(chǎng)的起伏、特別是垂直氣流的強(qiáng)弱變化對(duì)遷飛高度的變化起重要的作用。(4)下沉氣流和降水是稻縱卷葉螟降落的關(guān)鍵動(dòng)力因素,兩者都對(duì)降蟲(chóng)有明顯影響。(5)在這一北遷過(guò)程中,稻縱卷葉螟種群多降落在相對(duì)濕度大的區(qū)域,降蟲(chóng)區(qū)的相對(duì)濕度均在75%以上。

    稻縱卷葉螟;遷入;WRF模式;遷飛軌跡;大氣背景

    稻縱卷葉螟,Cnaphalocrocismedinalis,是廣泛分布在亞洲和東非水稻生長(zhǎng)區(qū)的一種遠(yuǎn)距離遷飛性害蟲(chóng)。我國(guó)作為世界上水稻生產(chǎn)和消費(fèi)的最大國(guó)家,自20世紀(jì)60年代以來(lái),一直遭受著其危害和困擾。稻縱卷葉螟在我國(guó)的分布較廣,向北可達(dá)吉林,向南可至海南、中國(guó)臺(tái)灣,東至沿海島嶼,西至云貴高原,均有觀測(cè)記載。每年3月份開(kāi)始,其種群從中南半島隨西南氣流遷入我國(guó)兩廣南部和云南南部地區(qū),隨后不斷向北深入,夏秋之交開(kāi)始南遷,10月下旬后陸續(xù)遷出我國(guó)[1]。我國(guó)稻縱卷葉螟的主要危害區(qū)在淮河以南的水稻主產(chǎn)區(qū),尤其以華南、江嶺和江淮稻區(qū)為甚[2]。本世紀(jì)以來(lái),稻縱卷葉螟在我國(guó)各稻區(qū)的發(fā)生蟲(chóng)量和危害面積明顯上升,自2003年來(lái)多年大發(fā)生,造成了我國(guó)水稻產(chǎn)量的巨大損失,已成為目前我國(guó)糧食安全的重大威脅。因此,揭示稻縱卷葉螟的發(fā)生規(guī)律,探明影響其發(fā)生和災(zāi)變的因素,提高稻縱卷葉螟發(fā)生的測(cè)報(bào)水平,對(duì)于保障我國(guó)糧食安全、維護(hù)社會(huì)穩(wěn)定有著重要的意義。

    稻縱卷葉螟屬鱗翅目螟蛾科,趨光、趨濕、趨嫩,具有遠(yuǎn)距離遷飛性,且具有較強(qiáng)的再遷飛能力,能經(jīng)過(guò)幾個(gè)夜晚的多次飛行,實(shí)現(xiàn)遠(yuǎn)距離轉(zhuǎn)移[3],這給監(jiān)測(cè)預(yù)報(bào)工作帶來(lái)了很大的困難。經(jīng)過(guò)國(guó)內(nèi)外學(xué)者多年的研究,已經(jīng)明確了影響遷飛性害蟲(chóng)災(zāi)變性遷入的重要因素有害蟲(chóng)自身的生理生態(tài)特性、寄主環(huán)境、大氣背景及地理?xiàng)l件等,其中它們的遠(yuǎn)距離遷飛和輾轉(zhuǎn)成災(zāi)只能借助適宜的氣象條件才能實(shí)現(xiàn)。因此,在全面了解害蟲(chóng)生理生態(tài)特征的基礎(chǔ)上,運(yùn)用現(xiàn)代氣象學(xué)理論和技術(shù),系統(tǒng)地研究和闡明大氣背景對(duì)稻縱卷葉螟災(zāi)變性遷入的影響機(jī)理顯得尤為重要。國(guó)內(nèi)外已有許多學(xué)者運(yùn)用氣象學(xué)理論和方法研究了影響遷飛性害蟲(chóng)災(zāi)變的大氣背景,如有研究指出不同類(lèi)型的天氣系統(tǒng)如鋒面、副熱帶高壓、大陸冷高壓等對(duì)褐飛虱降落有不同的影響[4],有研究發(fā)現(xiàn)氣象條件如風(fēng)場(chǎng)、垂直氣流、降水、溫度和濕度等是影響稻縱卷葉螟遷飛的重要因子[5]。Curtis R Wood等[6]曾用雷達(dá)觀測(cè)英國(guó)南部蛾類(lèi)夜間遷徙動(dòng)態(tài)時(shí)發(fā)現(xiàn),溫度、濕度以及風(fēng)速條件對(duì)蛾的起飛和成層遷飛有很大的影響。包云軒等[7]運(yùn)用MICAPS軟件繪制氣象背景圖,分析得出高空大氣環(huán)流形勢(shì)對(duì)稻縱卷葉螟的南、北遷降落起著重要的宏觀調(diào)控作用,高空風(fēng)場(chǎng)盛行偏南風(fēng)時(shí)對(duì)北遷有利,盛行偏北風(fēng)時(shí)對(duì)南遷有利,而降水條件則是稻縱卷葉螟集中遷入和降落的重要?jiǎng)恿γ{迫因子。王翠花等[2,8]分析了大氣動(dòng)力機(jī)制和水汽條件對(duì)2003年中國(guó)稻縱卷葉螟遷入成災(zāi)的影響,并指出高空的水汽輸送狀況和相對(duì)濕度分布對(duì)稻縱卷葉螟的遷入有著很強(qiáng)的指示意義。高月波等[9]在對(duì)稻縱卷葉螟的遷飛進(jìn)行多普勒雷達(dá)觀測(cè)的基礎(chǔ)上,推算了夏秋季發(fā)生在長(zhǎng)江中下游稻區(qū)的多個(gè)稻縱卷葉螟重大遷飛過(guò)程的后推軌跡,并簡(jiǎn)要分析了其遷入的大氣背景。前人主要從二維大氣背景場(chǎng)的水平分布上探討氣象條件對(duì)稻縱卷葉螟遠(yuǎn)距離遷飛的影響,極少結(jié)合多個(gè)大氣物理量場(chǎng)的垂直梯度變化從三維動(dòng)態(tài)上來(lái)研究大氣背景場(chǎng)對(duì)稻縱卷葉螟遷飛和再遷飛的影響,也鮮少運(yùn)用WRF模式輸出的精細(xì)化大氣背景場(chǎng)來(lái)系統(tǒng)分析各大氣物理量場(chǎng)對(duì)稻縱卷葉螟災(zāi)變性遷入的影響。

    本文針對(duì)過(guò)去研究中存在的問(wèn)題,在分析2000—2012年我國(guó)水稻主產(chǎn)區(qū)稻縱卷葉螟遷入主要特征的基礎(chǔ)上,選取2007年7月下旬一次典型北遷過(guò)程作為個(gè)例逆推稻縱卷葉螟的遷飛軌跡,運(yùn)用中尺度天氣研究和預(yù)報(bào)模式WRF,結(jié)合NCEP氣象再分析資料,精細(xì)化地模擬這一過(guò)程的三維大氣背景場(chǎng),分析其動(dòng)態(tài)變化對(duì)稻縱卷葉螟災(zāi)變性遷入的影響,旨在探明導(dǎo)致我國(guó)稻縱卷葉螟災(zāi)變性遷入的大氣物理機(jī)制,為我國(guó)遷飛性害蟲(chóng)的災(zāi)變監(jiān)測(cè)和預(yù)警提供科學(xué)依據(jù)。

    1 資料與方法

    1.1 資料

    1.1.1 蟲(chóng)情資料

    稻縱卷葉螟蟲(chóng)情資料是由農(nóng)業(yè)部全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心提供的全國(guó)168個(gè)病蟲(chóng)測(cè)報(bào)站的逐日燈誘數(shù)據(jù),資料時(shí)長(zhǎng)為2000—2012年,選取了資料較為完整且具有代表性的45個(gè)站點(diǎn)用于分析稻縱卷葉螟的發(fā)生特征。在此基礎(chǔ)上,選取了2007年7月下旬一次重大遷入過(guò)程(2007年7月22—24日)作為典型的災(zāi)變性北遷個(gè)例。

    1.1.2 氣象資料

    本文降水場(chǎng)實(shí)測(cè)資料由國(guó)家氣象信息中心提供。運(yùn)行WRF模式所用的初始場(chǎng)資料為美國(guó)國(guó)家環(huán)境預(yù)測(cè)中心(NCEP)提供的每6 h一次的水平分辨率為1°×1°的氣象再分析格點(diǎn)資料。

    1.1.3 基礎(chǔ)地理信息

    研究所用的基礎(chǔ)地理圖層為國(guó)家基礎(chǔ)地理信息中心提供的1∶4000000中國(guó)電子地圖。

    1.2 方法

    1.2.1 稻縱卷葉螟燈誘數(shù)據(jù)的處理

    首先,在Excel中整理出歷年各站點(diǎn)稻縱卷葉螟逐日燈誘資料,繪制出各站點(diǎn)稻縱卷葉螟歷年遷入蟲(chóng)量的時(shí)間變化曲線(xiàn),分析得出13年來(lái)各站遷入的始見(jiàn)期、南北遷高峰期、峰次數(shù)、終見(jiàn)期。因農(nóng)業(yè)部病蟲(chóng)測(cè)報(bào)規(guī)范中常以候?yàn)閱挝唬?,每?5日以前,每5d為1候,最后一候可以是3、4、5d或6d,本文在統(tǒng)計(jì)發(fā)生期時(shí)也應(yīng)用單位候。其次,在ArcGIS中繪制逐日稻縱卷葉螟燈誘蟲(chóng)量的空間分布圖,通過(guò)同期突增現(xiàn)象分析單站蟲(chóng)情遷入動(dòng)態(tài),并初步判斷出蟲(chóng)源地、遷飛路徑和降蟲(chóng)區(qū)。最后,從稻縱卷葉螟大發(fā)生的2007年燈誘資料中篩選出一個(gè)典型的重大遷飛過(guò)程2007年7月22—24日作為本文研究個(gè)例,詳細(xì)分析其遷飛的時(shí)空動(dòng)態(tài)。

    1.2.2 軌跡分析方法

    本文稻縱卷葉螟的遷飛軌跡計(jì)算是在NOAA官網(wǎng)(www.noaa.gov)上的HYSPLIT- 4軟件平臺(tái)上進(jìn)行的,采用后推示蹤法推算遷出蟲(chóng)源地和遷飛路徑,采用前推法分析主、副降蟲(chóng)區(qū)和遷飛路徑。采用UTC(標(biāo)準(zhǔn)世界時(shí)),可以每1 h輸出1組軌跡參數(shù),這些參數(shù)包括軌跡點(diǎn)所在的緯度、經(jīng)度、海拔高度等。模式中用于分析的氣象數(shù)據(jù)為GDAS格點(diǎn)值,NOAA官網(wǎng)同樣提供該類(lèi)數(shù)據(jù)的下載。選取了主、副降蟲(chóng)區(qū)的4個(gè)代表性站點(diǎn)(安徽太湖、江蘇張家港、江蘇丹陽(yáng)、浙江象山)為逆推軌跡起點(diǎn),進(jìn)行蟲(chóng)源地的逆推分析;考慮到盛夏7月稻縱卷葉螟的起飛峰時(shí)在18:00前后、降落峰時(shí)06:00前后,選取逆推時(shí)段為當(dāng)日06:00逆推至前一日18:00;根據(jù)昆蟲(chóng)雷達(dá)觀測(cè)結(jié)果[9- 10],稻縱卷葉螟盛夏的飛行高度變化于100—1500 m之間,聚集遷飛層在750 m上下。因此,逆推軌跡起點(diǎn)高度設(shè)為750 m,據(jù)此逆推出3個(gè)夜間12 h的遷飛軌跡線(xiàn),并在Excel中繪制出相應(yīng)的逆推軌跡高度,以反映稻縱卷葉螟的三維遷飛動(dòng)態(tài)。在分析軌跡及軌跡落點(diǎn)時(shí),考慮到該昆蟲(chóng)的生物學(xué)習(xí)性,軌跡終止落點(diǎn)必須在水稻種植區(qū)內(nèi),且該稻區(qū)正值稻縱卷葉螟遷出高峰[11]。

    1.2.3 個(gè)例大氣背景的數(shù)值模擬

    本文應(yīng)用WRF-ARW3.4模式對(duì)所選稻縱卷葉螟重大遷飛過(guò)程的大氣背景進(jìn)行數(shù)值模擬,采用NCEP 1°×1°氣象再分析資料作為初始場(chǎng),設(shè)置雙重嵌套區(qū)域,區(qū)域中心設(shè)為(116°E,28°N),外層網(wǎng)格水平方向?yàn)?22×200個(gè)格點(diǎn),格距為30 km;內(nèi)層網(wǎng)格水平方向?yàn)?80×229個(gè)格點(diǎn),格距為10 km;嵌套區(qū)域覆蓋我國(guó)稻縱卷葉螟各發(fā)生區(qū)域,模式積分步長(zhǎng)為60 s,考慮到稻縱卷葉螟的夜間遷飛特點(diǎn)及遷飛的時(shí)效性,模擬初始時(shí)刻選為2007年7月22日08:00(北京時(shí)),積分36 h,1 h輸出一次結(jié)果,在選擇微物理方案時(shí),內(nèi)外層微物理過(guò)程均采用WSM3簡(jiǎn)單冰方案,長(zhǎng)波輻射采用RRTM方案,短波輻射采用Dudhia方案,邊界層采用YSU方案,積云參數(shù)化采用Grell方案,模式輸出結(jié)果用氣象專(zhuān)用繪圖軟件Grads 2.0繪制各物理量場(chǎng)。

    2 結(jié)果與分析

    2.1 2000—2012年中國(guó)稻縱卷葉螟發(fā)生特點(diǎn)

    本文分析了2000—2012年中國(guó)稻縱卷葉螟燈誘資料,得出13年來(lái)平均的始見(jiàn)期、高峰期、終見(jiàn)期,并制作成圖1所示的GIS空間分布圖。

    稻縱卷葉螟的全國(guó)各地始見(jiàn)期分布圖(圖1),3月份,我國(guó)南方稻區(qū)西南季風(fēng)開(kāi)始盛行[12- 13],將境外的蟲(chóng)源不斷帶入我國(guó)境內(nèi)。這13年中,常年燈下見(jiàn)蟲(chóng)為3月5日—10日,最早的為3月1日,出現(xiàn)在2007年廣東陽(yáng)春,最遲的為3月25日,出現(xiàn)在2002年廣東陽(yáng)春、廣西合浦。在經(jīng)過(guò)平均時(shí)長(zhǎng)35候的順序北遷后,可蔓延至我國(guó)淮河以南的整個(gè)南方水稻主產(chǎn)區(qū),其間14—21候遷入兩廣南部;22—27候遷入兩廣北部、云南、貴州、湖南、福建和江西南部;28—33候遷入重慶、江西北部、湖北南部、安徽南部和浙江南部;34—39候遷入湖北北部、安徽北部、浙江北部及江蘇大部。

    稻縱卷葉螟的北遷高峰期分布圖(圖1),在西南季風(fēng)(特別是低空急流)的輸送下,伴隨著西太平洋副高的北抬西伸,從南到北我國(guó)各稻區(qū)依次出現(xiàn)北遷高峰期:25—33候,廣西東南部和廣東西南部首先出現(xiàn)北遷高峰期;緊接著34—39候,西南稻區(qū)進(jìn)入北遷高峰期;37—42候江嶺稻區(qū)經(jīng)歷北遷高峰期;43—47候江淮地區(qū)最后進(jìn)入北遷高峰期。

    南遷高峰期分布圖(圖1),和北遷相比,南遷具有遷入快、時(shí)間短、范圍小的特點(diǎn)。該時(shí)期,副高南撤,西南季風(fēng)和東南季風(fēng)開(kāi)始迅速衰退,偏北季風(fēng)開(kāi)始盛行,此時(shí)的運(yùn)載氣流有利于稻縱卷葉螟向南回遷,從北向南,江淮稻區(qū)南遷高峰期在50—51候;江嶺稻區(qū)在50—55候進(jìn)入南遷高峰期;華南稻區(qū)在56—59候進(jìn)入南遷高峰期;西南稻區(qū)只有貴州有回遷蟲(chóng)源出現(xiàn),其南遷高峰期在50—51候,云南和重慶地區(qū)由于海拔高、地形障礙和稻縱卷葉螟南遷高度比一般地形低等因素很難監(jiān)測(cè)到回遷蟲(chóng)源,而使得這兩個(gè)地區(qū)的稻縱卷葉螟遷入在8月底就已終見(jiàn)。

    稻縱卷葉螟的終見(jiàn)期分布圖(圖1),稻縱卷葉螟的遷入終見(jiàn)期最先于50—55候在江淮稻區(qū)出現(xiàn);其后是江嶺稻區(qū),發(fā)生在53—58候;59—62候華南稻區(qū)進(jìn)入終見(jiàn)期。

    本文對(duì)稻縱卷葉螟南北遷高峰次數(shù)進(jìn)行了一次統(tǒng)計(jì),北遷持續(xù)時(shí)間長(zhǎng),在各地峰次較多,峰期較長(zhǎng),從圖2看出,華南稻區(qū)常年平均有3—4個(gè)遷入峰,西南稻區(qū)一般有3—4個(gè)遷入峰,江嶺稻區(qū)一般有2—3個(gè)遷入峰,江淮稻區(qū)一般出現(xiàn)2—3個(gè)遷入峰;華南稻區(qū)是峰次發(fā)生較多的一個(gè)區(qū)域,華南稻區(qū)溫度、日照、降水適宜,稻田多分布在江河平原和丘陵谷底,多生長(zhǎng)以雙季稻為主的一年多熟水稻,豐富的食料和有利的地理?xiàng)l件使該地區(qū)成為稻縱卷葉螟北遷重要的輸送必經(jīng)之地,而江淮地區(qū)發(fā)生峰次則較少,該地區(qū)大部分水稻一季有余,兩季不足,雖峰次不多,但該地區(qū)有著地理優(yōu)勢(shì),尤其長(zhǎng)江中下游平原這一魚(yú)米之鄉(xiāng),在季風(fēng)影響之下,氣候溫潤(rùn),降水豐富,豐富的食料環(huán)境通常導(dǎo)致該地區(qū)稻縱卷葉螟暴發(fā)成災(zāi)。

    稻縱卷葉螟南遷的周期短于北遷周期,通常在9月初開(kāi)始回遷,至11月初結(jié)束,從圖2看出,江淮稻區(qū)常出現(xiàn)1—2個(gè)遷入峰,江嶺稻區(qū)一般有1—2個(gè)遷入峰,華南稻區(qū)常有2—3個(gè)遷入峰,西南稻區(qū)常年有1個(gè)遷入峰。

    圖1 2000—2012年我國(guó)稻縱卷葉螟遷入期的時(shí)空分布圖

    圖2 2000—2012年我國(guó)不同稻區(qū)稻縱卷葉螟遷入峰次圖

    圖3 2007年7月22—24日我國(guó)稻縱卷葉螟燈誘量空間分布圖

    2.2 2007年北遷個(gè)例分析

    2.2.1 稻縱卷葉螟蟲(chóng)情分析

    2007年是稻縱卷葉螟在我國(guó)稻區(qū)發(fā)生較重的一年,江淮稻區(qū)是主要受災(zāi)害的地區(qū),田間蟲(chóng)量之大,災(zāi)害持續(xù)時(shí)間之長(zhǎng),實(shí)屬罕見(jiàn),對(duì)水稻生產(chǎn)造成了嚴(yán)重的威脅[14]。通過(guò)在Excel中對(duì)2007年的稻縱卷葉螟燈誘數(shù)據(jù)進(jìn)行單站時(shí)間序列分析,可以知道2007年稻縱卷葉螟始見(jiàn)期早、終見(jiàn)期遲、在我國(guó)稻區(qū)為害時(shí)間長(zhǎng)的特點(diǎn)。燈誘始見(jiàn)期最早出現(xiàn)在廣東陽(yáng)春站,為3月1日,最晚出現(xiàn)在江蘇靖江,為7月4日;燈誘終見(jiàn)期最早出現(xiàn)在重慶墊江,為8月20日,最晚出現(xiàn)在廣西大多數(shù)站點(diǎn),均為11月10日。2007年7月23日,江淮稻區(qū)和江南稻區(qū)出現(xiàn)了一次同期突增現(xiàn)象,迎來(lái)了一次遷入高峰期,這一峰值時(shí)期也正與北遷高峰期統(tǒng)計(jì)結(jié)果一致,在此后的10 d內(nèi),這兩個(gè)稻區(qū)多個(gè)植保站的燈誘蟲(chóng)量居高不下。

    圖3是7月22日—24日稻縱卷葉螟燈誘蟲(chóng)量空間分布圖,可以看出,這幾天我國(guó)南方稻區(qū)都有蟲(chóng)源遷入,22日降蟲(chóng)區(qū)域主要分布在皖、湘、黔以及浙江、福建沿海地區(qū),降蟲(chóng)量相對(duì)較少,分布比較均勻,而23日降蟲(chóng)顯著增加,尤其以江、浙、皖稻區(qū)最明顯,浙江沿海站點(diǎn)日降蟲(chóng)量竟達(dá)到了上萬(wàn)頭,對(duì)比接下來(lái)幾天的降蟲(chóng)分布圖(2007-07-24),可以很清晰的看到此次同期突增持續(xù)時(shí)間較長(zhǎng)、范圍廣,主要降蟲(chóng)區(qū)域?yàn)檎憬睾5貐^(qū)和長(zhǎng)江中下游稻區(qū)的蘇、皖、鄂地區(qū),主要集中在30°N南北,副降蟲(chóng)區(qū)為湘、黔和閩東地區(qū),通過(guò)比較分析,可以確定7月23日是大氣背景對(duì)稻縱卷葉螟遷入產(chǎn)生作用的關(guān)鍵日。

    2.2.2 稻縱卷葉螟遷飛軌跡分析

    圖4 2007年7月23日06:00起、22日06:00起和21日06:00起太湖、象山、張家港和丹陽(yáng)站稻縱卷葉螟各夜夜間12h逆推軌跡的連線(xiàn)

    為了了解此次長(zhǎng)江中下游稻區(qū)稻縱卷葉螟同期突增的遷入過(guò)程,選取了代表性站點(diǎn)進(jìn)行高峰日逆推軌跡模擬,從圖4可以看出,23日蘇皖地區(qū)的蟲(chóng)源主要來(lái)源于其西南方向,逆推蟲(chóng)源地集中在江嶺稻區(qū)中部,說(shuō)明稻縱卷葉螟在這3個(gè)夜晚的遷飛中主要受西南氣流的影響,由湖南南部、江西中部輾轉(zhuǎn)至江蘇南部、安徽、浙江地區(qū),進(jìn)而到達(dá)長(zhǎng)江中下游地區(qū),而沿海稻區(qū)象山是此次重發(fā)生區(qū)域,其燈誘遷入量達(dá)到上萬(wàn)頭,其首次逆推路徑來(lái)自南路,逆推起點(diǎn)正好落在海上,無(wú)法判斷出蟲(chóng)源地,因此,在沿海地區(qū)降落的稻縱卷葉螟的遷飛路徑是復(fù)雜的。

    本文以750 m作為逆推起點(diǎn),從逆推軌跡高度變化來(lái)看(圖5),第一次逆推(2007-07-23),到達(dá)太湖站的蟲(chóng)群隨風(fēng)高度先降低后升高,丹陽(yáng)站的蟲(chóng)群則隨風(fēng)而上升,而到達(dá)張家港的蟲(chóng)群則一開(kāi)始就在較高的高度,此后的高度一致穩(wěn)定在750—800 m之間;第二次逆推(2007-07-22)、第三次逆推(2007-07-21)蟲(chóng)群遷飛高度也呈現(xiàn)一定的起伏變化,但遷飛高度大多集中在550—850 m之間。

    圖5 2007年7月23日、22日和21日分別于06:00起太湖、象山、張家港和丹陽(yáng)站稻縱卷葉螟各夜夜間12 h逆推軌跡高度曲線(xiàn)

    2.2.3 大氣環(huán)流形勢(shì)分析

    結(jié)合7月22日08:00 500 hPa天氣圖(圖6a)和7月23日08:00 500 hPa天氣形勢(shì)圖(圖略)可以發(fā)現(xiàn):500 hPa高度上,歐亞大陸高空維系著兩槽一脊的形勢(shì),兩低壓分別位于新西伯利亞西部和鄂霍次克海上空,高壓脊位于貝加爾湖上空附近,中高緯阻高處于減弱崩潰時(shí)期。在我國(guó)境內(nèi),河套地區(qū)有一切斷低壓,冷渦形勢(shì)不斷發(fā)展將會(huì)對(duì)我國(guó)北方地區(qū)造成降溫和不穩(wěn)定的雷陣雨天氣,在東部沿海地區(qū)西太平洋副熱帶高壓控制著粵東、福建大部和浙南地區(qū),并逐漸向北向西推移,副高脊線(xiàn)位于25°N—28°N之間,西南暖濕氣流和北方低壓后部南下的西北干冷氣流在副高西北側(cè)相遇,形成對(duì)流,有利于降水的發(fā)生。7月22日08:00(圖6b)以及23日08:00(圖略)700 hPa天氣圖上,可以清楚地看到江淮流域北部與朝鮮半島南部之間有一切變線(xiàn),切變線(xiàn)以南是我國(guó)水稻主產(chǎn)區(qū),盛行強(qiáng)勁的西南暖濕氣流,其中湘、鄂、皖、蘇、滬、浙等稻區(qū)的風(fēng)速已達(dá)到了低空急流的量級(jí),這對(duì)遷飛性害蟲(chóng)從蟲(chóng)源地大規(guī)模遷出十分有利。7月22日20:00 850 hPa(圖6c)天氣圖上顯示對(duì)流層低層高低壓配置與對(duì)流層中高層槽脊分布密切對(duì)應(yīng),我國(guó)長(zhǎng)江以南大部分地區(qū)受強(qiáng)盛的西南氣流控制,并與來(lái)自西北的冷性偏北風(fēng)形成切變線(xiàn),西南氣流源源不斷地將孟加拉灣和南海的水汽輸送至我國(guó)南方地區(qū),形成大片濕區(qū)范圍。

    圖6 不同高度位勢(shì)高度場(chǎng)(gpm)和風(fēng)場(chǎng)(m/s)

    2.2.4 稻縱卷葉螟遷飛的大氣動(dòng)力背景分析

    (1)水平流場(chǎng)分析

    7月22日20:00(圖7a)925 hPa風(fēng)場(chǎng)上,來(lái)自南海的偏南風(fēng)直入我國(guó)兩廣地區(qū),西南氣流在我國(guó)兩廣地區(qū)和湘、贛南部地區(qū)強(qiáng)勁而有力,深入內(nèi)陸,風(fēng)速減緩,海上黃海中部以南和陸上淮河以南的東部地區(qū)盛行較一致的強(qiáng)西南氣流,但此時(shí)我國(guó)陸地上的風(fēng)速并未達(dá)到低空急流的量級(jí),隨著整個(gè)環(huán)流形勢(shì)的發(fā)展,23日02:00(圖7b)925 hPa上風(fēng)場(chǎng)圖上,我國(guó)湘南地區(qū)出現(xiàn)了低空急流,且我國(guó)大部分地區(qū)均以西南氣流為主,風(fēng)速較大,對(duì)稻縱卷葉螟從南部稻區(qū)的北遷十分有利,而在淮河沿岸,西南氣流與東亞大槽后部南下的東北氣流交綏,形成對(duì)流,并產(chǎn)生降水。23日08:00(圖7c)925 hPa風(fēng)場(chǎng)上,兩廣南部有西南暖濕氣流不斷向東北方向輸送,湘南地區(qū)的西南低空急流向東北方向移動(dòng),已經(jīng)影響到贛北地區(qū)。而自東北方向南下的干冷氣流不斷加強(qiáng)南壓,逐漸影響到江淮稻區(qū)上空的風(fēng)場(chǎng)格局,皖南、蘇南及浙北地區(qū)氣流逐漸由西南向西偏,由于西南氣流與東北氣流的交綏導(dǎo)致蘇、皖、鄂形成降水帶,更有利于稻縱卷葉螟的降落成災(zāi)。圖8a是02:00沿30°N各層次水平風(fēng)場(chǎng)的緯向剖面圖,可以看出113°E—125°E之間風(fēng)向均為西南風(fēng),且925 hPa高度風(fēng)速達(dá)到10 m/s以上,高度越高風(fēng)速越大;圖8b是02:00沿120°E各層次水平風(fēng)場(chǎng)的經(jīng)向剖面圖,可以看出23°N—33°N之間風(fēng)向均為西南風(fēng),925 hPa高度風(fēng)速也達(dá)10 m/s以上,而這兩個(gè)范圍所形成的的集中區(qū)域正是稻縱卷葉螟降蟲(chóng)最頻繁的區(qū)域,主要遷飛層的風(fēng)向以西南風(fēng)為主,風(fēng)速大小也對(duì)稻縱卷葉螟向東北方向的遷飛極為有利。

    圖7 2007年7月22日20:00(a)、7月23日02:00(b)和08:00(c)925 hPa風(fēng)場(chǎng)(m/s)

    圖8 2007年7月23日02:00沿30°N緯向風(fēng)場(chǎng)(m/s)剖面圖(a)和沿120°E經(jīng)向風(fēng)場(chǎng)(m/s)剖面圖(b)

    (2)垂直速度場(chǎng)分析

    稻縱卷葉螟在空中順風(fēng)而行時(shí)如遇一定強(qiáng)度的下沉氣流或降水,會(huì)被迫降落地面,從7月23日02:00和08:00 925 hPa垂直速度分布圖(圖9)上可以看出,23日夜間下沉氣流分布在30°N以南地區(qū):湘贛南部地區(qū)、東部沿海地區(qū)和皖鄂南部地區(qū),從02:00到08:00,下沉氣流區(qū)域呈縮小趨勢(shì)。從圖10中可以看出,沿30°N緯向分布,02:00(圖10a),850 hPa以上的高層主要以上升氣流為主,尤其是113°E附近有強(qiáng)的上升氣流;850 hPa以下的低層區(qū)域主要以下沉氣流為主,尤其是在113°E—120°E之間,存在著兩個(gè)較為明顯的下沉氣流區(qū)域,分別位于鄂南、皖南地區(qū)和浙北地區(qū)。到了08:00(圖10b),高層113°E—116°E之間下沉氣流明顯增強(qiáng),低層下沉氣流區(qū)范圍縮小,集中至114°E—117°E之間,即范圍縮至鄂南、皖南地區(qū)。圖11是垂直速度的緯向分布圖,02:00(圖11a),有利于稻縱卷葉螟降落的下沉氣流主要集中在25°N—29°N之間,即集中在福建和浙江的沿海地區(qū),而至08:00(圖11b),下沉氣流主要集中在30°N以南,對(duì)閩浙兩地區(qū)稻縱卷葉螟的降落極為有利。根據(jù)以上的分析,我國(guó)東南部水稻主產(chǎn)區(qū)下沉氣流的分布與稻縱卷葉螟的集中降落區(qū)基本吻合。

    圖9 2007年7月23日02:00(a)和08:00 (b) 925 hPa垂直速度分布(10-2 m/s)

    圖10 2007年7月23日02:00 (a)和08:00 (b)沿30°N緯向垂直速度剖面圖(10-2 m/s)

    圖11 2007年7月23日02:00 (a)和08:00 (b)沿120°E經(jīng)向垂直速度剖面圖(10-2 m/s)

    此外,垂直氣流的強(qiáng)弱變化還對(duì)稻縱卷葉螟空中遷飛種群的遷飛高度有著至關(guān)重要的影響,如7月23日08:00從安徽太湖站(30.494°N,116.1698°E)和江蘇丹陽(yáng)站(31.9581°N,119.6336°E)逆推的軌跡在有強(qiáng)上升氣流存在的時(shí)刻(02:00)和地區(qū),其遷飛高度上升(圖4、圖9a、圖10a、圖11a、圖5),而從張家港站(31.9002°N,120.6294°E)逆推的軌跡在有強(qiáng)下沉氣流存在的時(shí)刻(02:00)和地區(qū),其遷飛高度下降(圖4、圖9a、圖10a、圖11a、圖5a)。同樣,在其它時(shí)次、軌跡經(jīng)過(guò)的其它地區(qū)也出現(xiàn)了相同的現(xiàn)象。

    (3)降水場(chǎng)分析

    對(duì)于稻縱卷葉螟的降落,在不考慮稻縱卷葉螟由于生理原因主動(dòng)降落的情況下,大氣動(dòng)力脅迫(如下沉氣流、降水)和地形脅迫都是造成稻縱卷葉螟集中降落的重要機(jī)制[7]。在了解下沉氣流對(duì)此次降蟲(chóng)影響的基礎(chǔ)上,為了更全面地闡明稻縱卷葉螟的降落機(jī)制,將7月22日08:00至23日08:00各站點(diǎn)的降水量按照氣象部門(mén)規(guī)定的24h降雨量級(jí)別劃分后繪制出如圖12所示的空間分布圖。從圖12中可以看出,大部分稻區(qū)有少量的降水,在30°N以北降水區(qū)主要分布在東部,以蘇、皖、鄂地區(qū)為主;在鄂、皖、豫交界處以及湘鄂交界處都出現(xiàn)了強(qiáng)降水,雨量達(dá)到了大到暴雨級(jí)別;江蘇北部以及西南大部分稻區(qū)雨量大多為大雨級(jí)別;蘇、浙交界處雨量為小到中雨級(jí)別。在30°N以南降水區(qū)集中在海拔高的高原和山區(qū),即云南、四川和重慶地區(qū),但因?yàn)樵摰貐^(qū)海拔高于稻縱卷葉螟的主要遷飛高度,大部分地區(qū)并未見(jiàn)蟲(chóng)。因此,此次蘇皖部分地區(qū)的降蟲(chóng)是由降水引起的,降水條件能導(dǎo)致稻縱卷葉螟降落,但降蟲(chóng)的主降區(qū)并不在降水區(qū)域最強(qiáng)的區(qū)域,而是分布在降水區(qū)附近。

    圖12 2007年7月22日08:00至7月23日08:00中國(guó)南方水稻主產(chǎn)區(qū)降水分布圖(mm)

    (4)相對(duì)濕度場(chǎng)和水汽通量場(chǎng)分析

    稻縱卷葉螟是趨濕性遷飛害蟲(chóng),通過(guò)對(duì)比分析稻縱卷葉螟遷入的空間分布,以及夜間降水情況,可以發(fā)現(xiàn):此次主降蟲(chóng)區(qū)主要在沿江、沿湖和沿海地區(qū),在福建沿海地區(qū)有著密集的降蟲(chóng)區(qū),說(shuō)明稻縱卷葉螟的遷飛和降落與水汽條件有關(guān)。在氣象學(xué)中,相對(duì)濕度、水汽通量等可以表征一個(gè)地區(qū)的水汽條件。從7月22日20:00 925 hPa的相對(duì)濕度場(chǎng)上來(lái)看(圖13a),相對(duì)濕度低數(shù)值區(qū)主要集中在浙江沿海、福建沿海和湖南、江西兩省的大部分地區(qū),其值大約在55%—70%之間。而高值區(qū)集中在云南、兩廣地區(qū),其值在75%—95%左右,經(jīng)過(guò)12 h的醞釀,降蟲(chóng)區(qū)925 hPa相對(duì)濕度均超過(guò)了75%(圖13b),對(duì)比可以發(fā)現(xiàn),浙江沿海稻區(qū)相對(duì)濕度上升了20%之多,而其余稻區(qū)的相對(duì)濕度也上升了10%之多,圖15是08:00相對(duì)濕度沿30°N和120°E的垂直剖面圖,鄂南—皖南—浙北一帶(113°E—122°E)和福建—浙江—江蘇一帶(28°N—31°N)正是此次稻縱卷葉螟集中降落區(qū)域,稻縱卷葉螟遷飛層的相對(duì)濕度達(dá)到了80%—95%,這為稻縱卷葉螟的遷入成災(zāi)提供了絕好的水汽條件。

    再看水汽通量的分布(圖14a),7月22日20:00水汽通量最大值位于廣西南部的沿海地區(qū),通量值達(dá)到18 g/s以上,南海的水汽直接在西南氣流的輸送下,已經(jīng)開(kāi)始影響到了兩廣和湘南地區(qū)空氣中的水汽含量,使這幾個(gè)地區(qū)的水汽通量值達(dá)到10 g/s以上,另一個(gè)水汽通量大值區(qū)位于蘇南和浙北沿海地區(qū),其值達(dá)到了16 g/s,水汽來(lái)源于西北太平洋,其余地區(qū)的水汽通量普遍低于6 g/s,但都為正值,表明有水汽輸入。隨著水汽不斷向內(nèi)陸的輸送,經(jīng)過(guò)一夜,各地區(qū)的水汽通量格局發(fā)生了變化,如圖14b 所示,除了降水多的地區(qū)外,大部分地區(qū)的水汽通量值高于10 g/s,蘇皖鄂的雨帶地區(qū)水汽通量低于6 g/s,這是由于該地區(qū)降水條件的發(fā)生,水汽含量已趨飽和,與其他地區(qū)的水汽梯度相對(duì)較小。而來(lái)自南海的水汽通過(guò)廣西—湘南—贛南有一水汽輸送帶,不斷向北輸送,這正與稻縱卷葉螟的北遷主徑相吻合,贛南地區(qū)的水汽通量達(dá)到了較高值,高于20 g/s。另外,來(lái)自西北太平洋的水汽也不斷向江浙皖稻區(qū)輸送,水汽通量普遍高于10 g/s,圖15是08:00水汽通量沿30°N和120°E垂直剖面圖,降蟲(chóng)區(qū)域水汽通量維持在10 g/s以上,降水區(qū)對(duì)流層中低層水汽通量維持在較高的水平,此次我國(guó)南方稻區(qū)水汽主要有兩個(gè)來(lái)源:一是來(lái)自南海和孟加拉灣,由西南方向向東北方向輸送;二是來(lái)自東海海面,由東南偏南至西北偏北方向向陸地輸送,水汽條件對(duì)稻縱卷葉螟的遷入起著很好的輔助作用。

    圖13 2007年7月22日20:00(a)和23日08:00(b)925 hPa相對(duì)濕度場(chǎng)(%)

    圖14 2007年7月22日20:00(a)和23日08:00(b)925 hPa水汽通量場(chǎng)(g/s)

    圖15 2007年7月23日08:00相對(duì)濕度(%,等值線(xiàn))和水汽通量(g/s,彩色陰影)沿30°N(a)和120°E(b)的垂直剖面圖

    3 討論

    稻縱卷葉螟南北往返的遷入和為害與我國(guó)的水稻種植制度、地理?xiàng)l件以及大氣背景都有著很重要的關(guān)系,本文利用2000—2012年中國(guó)稻縱卷葉螟逐日燈誘數(shù)據(jù)分析了最近10多年來(lái)稻縱卷葉螟在我國(guó)水稻主產(chǎn)區(qū)的遷入特征,發(fā)現(xiàn):稻縱卷葉螟在我國(guó)的遷入時(shí)空特征與以往基本相似[15],但最早始見(jiàn)期有所提前,最晚終見(jiàn)期有所推遲,特別是2003年,這與全球氣候變暖和氣候異常有關(guān)[16]。軌跡分析顯示:此次過(guò)程的稻縱卷葉螟蟲(chóng)源地基本上都位于降蟲(chóng)區(qū)的西南方,這與這一時(shí)期西南季風(fēng)的盛行和湖南、江西及安徽南部稻區(qū)的早稻黃熟有關(guān);遷飛高度的起伏變化范圍則由稻縱卷葉螟的生理生態(tài)特征和遷飛行為學(xué)參數(shù)及這一時(shí)期遷飛層的大氣物理結(jié)構(gòu)所決定的,也與以往的昆蟲(chóng)雷達(dá)觀測(cè)結(jié)果相吻合[9-10]。沿海降蟲(chóng)區(qū)復(fù)雜的軌跡走向是由大尺度大氣環(huán)流背景、中尺度天氣過(guò)程和局地海陸風(fēng)環(huán)流共同作用造成的,需要耦合WRF模式與海陸風(fēng)模式進(jìn)一步探明,本文限于篇幅,沒(méi)有展開(kāi)。925 hPa等壓面所處高度約為750 m,正好是稻縱卷葉螟的適遷層,該高度上的水平氣流速度遠(yuǎn)大于稻縱卷葉螟主動(dòng)遷飛的速度,因此,受氣流脅迫,種群只能隨風(fēng)而行,此次北遷過(guò)程的925 hPa上南方稻區(qū)盛行強(qiáng)而一致的西南氣流,這就決定了稻縱卷葉螟種群遷向東北。雖然垂直氣流的強(qiáng)弱與水平氣流相比,差了1—2個(gè)量級(jí),但兩者合成以后產(chǎn)生的慣性對(duì)北遷方向及遷飛高度的變化起了至關(guān)重要的作用。稻縱卷葉螟晝伏夜行,可以經(jīng)過(guò)多次輾轉(zhuǎn)飛行到達(dá)受災(zāi)地,在其異地轉(zhuǎn)移時(shí)需要借助有利的大氣動(dòng)力條件,高空大氣環(huán)流形勢(shì)、風(fēng)場(chǎng)、垂直速度場(chǎng)、降水無(wú)疑是影響稻縱卷葉螟遷飛降落的關(guān)鍵動(dòng)力因素,本文通過(guò)一次典型北遷個(gè)例的分析,探討大氣動(dòng)力場(chǎng)對(duì)稻縱卷葉螟遷入降落的影響,所得結(jié)論與以往的研究基本吻合[8],而三維流場(chǎng)特別是垂直氣流的強(qiáng)弱變化對(duì)稻縱卷葉螟遷飛高度的影響卻是初次探討,變化大氣動(dòng)力條件和水汽條件對(duì)稻縱卷葉螟的遷入成災(zāi)均有很好的指示作用:水平氣流是稻縱卷葉螟遠(yuǎn)距離遷飛的主要運(yùn)載動(dòng)力,下沉氣流和降水是稻縱卷葉螟降落的關(guān)鍵動(dòng)力因素,兩者都對(duì)降蟲(chóng)有明顯的影響。稻縱卷葉螟是趨濕性害蟲(chóng),適宜的水分條件有助于其遠(yuǎn)距離遷飛和降落,稻縱卷葉螟遷入種群多降落在濕度大的區(qū)域,降蟲(chóng)區(qū)地面的相對(duì)濕度均在75%以上,這與其對(duì)大氣水分環(huán)境的生態(tài)要求(70%以上相對(duì)濕度才能正常生存)是吻合的,與筆者以往的研究[2]不一樣的是:以往只間斷性地考慮了850、925 hPa和1000 hPa三個(gè)高度的相對(duì)濕度條件,而本文從緯向和經(jīng)向垂直剖面上連續(xù)性地考慮了稻縱卷葉螟遷飛發(fā)生區(qū)的相對(duì)濕度條件,這樣更科學(xué)、更全面。以往的研究發(fā)現(xiàn):稻縱卷葉螟的遷入與溫度的變化直接相關(guān),適宜稻縱卷葉螟遷入的溫度大于18 ℃[5],但盛夏稻縱卷葉螟北遷的溫場(chǎng)條件極易滿(mǎn)足,對(duì)其遷出、空中飛行和降落很難形成脅迫。而本文所選個(gè)例中的遷飛過(guò)程正好發(fā)生在盛夏時(shí)節(jié),稻縱卷葉螟遷入所需的熱量條件不成問(wèn)題,故沒(méi)有考慮溫度條件對(duì)其遷入的限制。

    HYSPLIT模型是用于計(jì)算和分析大氣污染物輸送、擴(kuò)散軌跡的專(zhuān)業(yè)模型,隨著模型的不斷完善,它開(kāi)始被廣泛應(yīng)用于昆蟲(chóng)遷飛軌跡研究及蟲(chóng)源地分析中[17- 18]。通過(guò)這一軌跡分析平臺(tái),可以很清楚的追溯稻縱卷葉螟的遷入蟲(chóng)源地分布并查看各時(shí)空點(diǎn)的遷飛高度。本文利用該模型推算軌跡時(shí),充分考慮了稻縱卷葉螟的較強(qiáng)再遷飛能力,同時(shí)根據(jù)前人的昆蟲(chóng)雷達(dá)觀測(cè)結(jié)果[9- 10],選取相關(guān)軌跡計(jì)算判斷依據(jù)(包括起始和終止時(shí)間、盛行遷飛層高度、軌跡落點(diǎn)區(qū)域等),與從實(shí)際蟲(chóng)情數(shù)據(jù)分析推斷的結(jié)果有較好的一致性。但軌跡起終點(diǎn)、遷飛路徑和遷飛距離與實(shí)際蟲(chóng)源區(qū)、遷飛路徑和降蟲(chóng)區(qū)的偏差還是難免的,原因是軌跡計(jì)算模型中沒(méi)有考慮稻縱卷葉螟在起飛、空中飛行和降落過(guò)程中的主動(dòng)性。尤其是在推算降落在沿海稻區(qū)的稻縱卷葉螟的遷飛路徑和蟲(chóng)源地時(shí),較容易出現(xiàn)軌跡落點(diǎn)出現(xiàn)在海上、導(dǎo)致逆推軌跡中斷的現(xiàn)象。后續(xù)研究將通過(guò)風(fēng)洞試驗(yàn)獲取其主動(dòng)遷飛速度與環(huán)境條件(如氣流、溫度、濕度等)關(guān)系的參數(shù),將其導(dǎo)入軌跡模型使之更為精確。WRF模式已廣泛用于世界上多數(shù)國(guó)家的天氣過(guò)程研究和數(shù)值預(yù)報(bào)業(yè)務(wù)中,其功能和效果得到國(guó)際氣象界的廣泛公認(rèn),因此,它對(duì)稻縱卷葉螟遷飛過(guò)程發(fā)生期間的大氣背景的模擬是準(zhǔn)確而可行的。在掌握大范圍基本蟲(chóng)源信息的前提下,可以利用它來(lái)對(duì)稻縱卷葉螟遷飛的軌跡和發(fā)生大氣背景進(jìn)行短期預(yù)測(cè),為了驗(yàn)證這種預(yù)測(cè)的可行性,需要通過(guò)大量的個(gè)例試驗(yàn)獲取可靠的準(zhǔn)確率和誤差數(shù)據(jù),只有達(dá)到一定的準(zhǔn)確率并將誤差控制在允許的范圍內(nèi)才能業(yè)務(wù)化應(yīng)用。

    [1] 王鳳英, 胡高, 陳曉, 沈慧梅, 羅善煜, 辛德育, 徐盛剛, 張孝羲, 翟保平. 近年來(lái)廣西南寧稻縱卷葉螟大發(fā)生原因分析. 中國(guó)水稻科學(xué), 2009, 23(5): 537- 545.

    [2] 王翠花, 包云軒, 王建強(qiáng), 翟保平, 向勇. 2003年稻縱卷葉螟大發(fā)生的水汽條件分析. 應(yīng)用生態(tài)學(xué)報(bào), 2006, 17(9): 1693- 1698.

    [3] 王鳳英, 張孝羲, 翟保平. 稻縱卷葉螟的飛行和再遷飛能力. 昆蟲(chóng)學(xué)報(bào), 2010, 53(11): 1265- 1272.

    [4] 談涵秋, 毛瑞曾, 程極益, 姚禾芬. 褐飛虱遠(yuǎn)距離遷飛中的降落和垂直氣流、降雨的關(guān)系. 南京農(nóng)學(xué)院學(xué)報(bào), 1984, (2): 18- 25.

    [5] 白先達(dá), 黃超艷, 唐廣田, 李國(guó)剛, 鄒麗霞, 賈方曌. 氣象條件對(duì)稻縱卷葉螟遷飛的影響分析. 中國(guó)農(nóng)學(xué)通報(bào), 2010, 26(21): 262- 267.

    [6] Wood C R, Chapman J W, Reynolds D R, Barlow J F, Smith A D, Woiwod L P. The influence of the atmospheric boundary layer on nocturnal layers of noctuids and other moths migrating over southern Britain. International Journal of Biometeorology, 2006, 50(4): 193- 204.

    [7] 包云軒, 王永平, 嚴(yán)明良, 袁成松, 于曉麗. 2003年我國(guó)稻縱卷葉螟發(fā)生特征及其災(zāi)變大氣背景的研究. 氣象科學(xué), 2008, 28(2): 184- 189.

    [8] 王翠花, 包云軒, 王建強(qiáng), 向勇, 翟保平. 2003年稻縱卷葉螟重大遷入過(guò)程的大氣動(dòng)力機(jī)制分析. 昆蟲(chóng)學(xué)報(bào), 2006, 49(4): 604- 612.

    [9] 高月波, 陳曉, 陳鐘榮, 包云軒, 楊榮明, 劉天龍, 翟保平. 稻縱卷葉螟(Cnaphalocrocismedinalis)遷飛的多普勒昆蟲(chóng)雷達(dá)觀測(cè)及動(dòng)態(tài). 生態(tài)學(xué)報(bào), 2008, 28(11): 5238- 5247.

    [10] Riley J R, Reynolds D R, Smith A D, Edwards A S, Zhang X X, Cheng X N, Wang H K, Cheng J Y, Zhai B P. Observations of the autumn migration of the rice leaf rollerCnaphalocrocismedinalis(Lepidoptera: Pyralidae) and other moths in eastern China. Bulletin of Entomological Research, 1995, 85(3): 397- 414.

    [11] 蔣春先, 齊會(huì)會(huì), 孫明陽(yáng), 武俊杰, 張?jiān)苹? 程登發(fā). 2010年廣西興安地區(qū)稻縱卷葉螟發(fā)生動(dòng)態(tài)及遷飛軌跡分析. 生態(tài)學(xué)報(bào), 2011, 31(21): 6495- 6504.

    [12] 包云軒, 謝杰, 向勇, 徐希燕, 劉宇. 低空急流對(duì)中國(guó)褐飛虱重大北遷過(guò)程的影響. 生態(tài)學(xué)報(bào), 2009, 29(11): 5773- 5782.

    [13] 包云軒, 黃金穎, 謝曉金, 陸明紅. 季風(fēng)進(jìn)退和轉(zhuǎn)換對(duì)中國(guó)褐飛虱遷飛的影響. 生態(tài)學(xué)報(bào), 2013, 33(16): 4864- 4877.

    [14] 劉宇,王建強(qiáng), 馮曉東, 蔣學(xué)輝.2007年全國(guó)稻縱卷葉螟發(fā)生實(shí)況分析與2008年發(fā)生趨勢(shì)預(yù)測(cè).中國(guó)植保導(dǎo)刊,2008,28(7):33-35.

    [15] 石尚柏, 姜瑞中. 農(nóng)作物病蟲(chóng)長(zhǎng)期運(yùn)動(dòng)規(guī)律與預(yù)測(cè). 北京: 中國(guó)農(nóng)業(yè)出版社, 1997.

    [16] 陳曉, 劉萬(wàn)才, 朱建生, 張孝羲, 翟保平. 2003年我國(guó)稻縱卷葉螟特大發(fā)生的氣候背景. 昆蟲(chóng)知識(shí), 2013, 50(3): 615- 621.

    [17] 齊國(guó)君, 蘆芳, 高燕, 梁居林, 藍(lán)日青, 張孝羲, 翟保平, 呂利華. 稻縱卷葉螟2010年的一次遷飛過(guò)程及其蟲(chóng)源分析. 昆蟲(chóng)學(xué)報(bào), 2011, 54(10): 1194- 1203.

    [18] 蘆芳, 齊國(guó)君, 陳曉, 董西貴, 郭玉人, 武向文, 倪秀紅, 陳時(shí)健, 張紀(jì)英, 張孝羲, 翟保平. 上海地區(qū)2007年褐飛虱的后期遷入和蟲(chóng)源地的個(gè)例分析. 生態(tài)學(xué)報(bào), 2010, 30(12): 3215- 3225.

    Migration pattern of rice leaf roller and impact of atmospheric conditions on a heavy migration event in China

    BAO Yunxuan1,*, CAO Yun1, XIE Xiaojin1, LU Minghong2, LI Xuan3, WANG Chunzhi3, LIU Wancai2

    1JiangsuKeyLaboratoryofAgriculturalMeteorology,NanjingUniversityofInformationScience&Technology,Nanjing210044,China2NationalAgriculturalTechnologyExtensionandServiceCenter,MinistryofAgricultural,Beijing100125,China3NationalMeteorologicalCenter,Bejing100081,China

    Rice leaf roller,CnaphalocrocismedinalisGuenee, is one of the most important crop pests in China due to its large area of occurrence, high frequency of outbreaks and serious damage to rice yields. In this paper, the lighting trap observational data of rice leaf rollers from 2000 to 2012 were analyzed to assess the impacts of atmospheric conditions on the migration and landfall of rice leaf roller. We have identified the main features of rice leaf roller′s migration in the main rice-growing areas of China and presented a detailed analysis on a typical migration event occurred in late July, 2007 and verified the assessment result. By means of the Hybrid Single Particle Lagrangian Integrated Trojectory Model for Version 4.0 (HYSPLIT- 4), a software of calculating particle trajectory in the air, the backward migration trajectories of the rice leaf roller landfalled the rice growing regions in the the Yangtze River Valley were calculated and the Weather Research and Forecast model (WRF) ′s simulations on the atmospheric backgrounds in this event combined with the reanalysis data of the National Ceter of Environmental Predicting (NCEP) were used to assess the impacts of weather conditions on the catastrophic migrations of rice leaf roller. Several important findings are obtained from this study. First, during the studying period from 2000 to 2012, the migration of rice leaf roller started at the beginning of March in a year, the period during March to August was the period of prevailing migration northward in which the populations immigrate into the south China, the southwest China, the regions between the Nanling Mountain and the Yangtze River and the valley regions between the Yangtze River and the Huaihe River in China successively and emerge a series of immigration peaks. During the period from September to November each year, the southward migration is the reversed processes in the above regions. The rice leaf roller migrated backward to the valley region between the Yangtze River and the Huaihe River, the regions between the Nanling Mountain and the Yangtze River and South China in turn. Second, the backward trajectories calculations of rice leaf roller′s migration northward included the beginning and re-emigrating populations in the selected event indicated that the sources of the rice leaf roller′s populations were situated at the southwest parts of the landfall regions and the migration altitudes varied from 550 m to 850 m above ground level. As a comparison, the migrating trajectories in the coastal regions were more complicated than the other regions. Third, the horizontal wind is the main driving force for the long-distance migration of rice leaf roller′s populations and the prevailing southerly at 925 hPa height is favorable for the migration northward. Three-dimensional windfield, especially the strong vertical airstream plays a dominant role in the up and down of the migration heights. Fourth, the downward flow and rainfall are the two key dynamic factors for the landfall of the rice leaf rollers. Fifth, as a kind of hygrotaxis pest, the suitable atmospheric moisture condition is beneficial to the long-distance migration and landfalling of rice leaf roller′s populations. In this event, most of the rice leaf rollers landfalled in the areas of the relative humidity higher than 75%.

    CnaphalocrocismedinalisGuenee; migration; WRF model; migration trajectory; atmospheric conditions

    國(guó)家自然科學(xué)基金面上項(xiàng)目(41075086);國(guó)家公益性行業(yè)(氣象)科研專(zhuān)項(xiàng)(GYHY201306053);江蘇省農(nóng)業(yè)科技自主創(chuàng)新項(xiàng)目(SCX(12)3058);江蘇省高校優(yōu)勢(shì)學(xué)科建設(shè)工程

    2013- 09- 26;

    2014- 06- 03

    10.5846/stxb201309262367

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

    包云軒,曹云,謝曉金,陸明紅,李軒,王純枝,劉萬(wàn)才.中國(guó)稻縱卷葉螟發(fā)生特點(diǎn)及北遷的大氣背景.生態(tài)學(xué)報(bào),2015,35(11):3519- 3533.

    Bao Y X, Cao Y, Xie X J, Lu M H, Li X, Wang C Z, Liu W C.Migration pattern of rice leaf roller and impact of atmospheric conditions on a heavy migration event in China.Acta Ecologica Sinica,2015,35(11):3519- 3533.

    猜你喜歡
    稻區(qū)卷葉螟水汽
    河南省粳稻新品種生態(tài)適應(yīng)性分析
    2023年全國(guó)水稻重大病蟲(chóng)害發(fā)生趨勢(shì)預(yù)報(bào)
    青藏高原上空平流層水汽的時(shí)空演變特征
    太倉(cāng)市近年稻縱卷葉螟發(fā)生情況及防治進(jìn)展
    棉大卷葉螟識(shí)別與為害癥狀
    縱卷葉螟絨繭蜂對(duì)稻縱卷葉螟的控制作用
    1979~2011年間平流層溫度及平流層水汽的演變趨勢(shì)
    深圳“5·11”特大暴雨過(guò)程的水汽輸送特征分析
    5%環(huán)蟲(chóng)酰肼防治稻縱卷葉螟田間藥效評(píng)價(jià)
    2015年全國(guó)水稻重大病蟲(chóng)害發(fā)生趨勢(shì)預(yù)報(bào)
    久久久久久伊人网av| 自线自在国产av| 日本wwww免费看| 成人毛片a级毛片在线播放| 国产成人精品福利久久| 两个人免费观看高清视频 | 国产精品国产三级国产专区5o| 五月天丁香电影| 国产在视频线精品| 午夜福利影视在线免费观看| 亚洲国产最新在线播放| 美女福利国产在线| 久久国产亚洲av麻豆专区| 日日摸夜夜添夜夜爱| 国产综合精华液| 中文字幕av电影在线播放| 国产免费一区二区三区四区乱码| 午夜91福利影院| 午夜激情久久久久久久| 一本色道久久久久久精品综合| 日韩精品有码人妻一区| 老司机影院毛片| 2018国产大陆天天弄谢| 少妇人妻一区二区三区视频| 日韩欧美 国产精品| 全区人妻精品视频| 久久人妻熟女aⅴ| 成人毛片60女人毛片免费| 伊人久久精品亚洲午夜| 亚洲无线观看免费| 欧美人与善性xxx| 亚洲国产精品999| www.av在线官网国产| 国产精品99久久久久久久久| 久久久久国产网址| 欧美老熟妇乱子伦牲交| 人人妻人人看人人澡| 久久久久久久久久成人| 国产日韩欧美在线精品| 欧美丝袜亚洲另类| a 毛片基地| 有码 亚洲区| 午夜激情久久久久久久| 一区二区av电影网| 99热国产这里只有精品6| 少妇人妻 视频| 女人精品久久久久毛片| 亚洲内射少妇av| 日韩强制内射视频| 亚洲国产精品一区三区| 黄色怎么调成土黄色| 精品久久久精品久久久| 久久久久精品性色| 久久 成人 亚洲| 99九九在线精品视频 | av视频免费观看在线观看| 亚洲精品一二三| 国产毛片在线视频| 国产精品.久久久| 久久国产精品大桥未久av | 一级毛片我不卡| 亚洲av免费高清在线观看| 熟女人妻精品中文字幕| 久久精品国产亚洲网站| 少妇猛男粗大的猛烈进出视频| av卡一久久| 国产精品久久久久成人av| 狂野欧美白嫩少妇大欣赏| 中文字幕精品免费在线观看视频 | 高清在线视频一区二区三区| 99久久精品国产国产毛片| 精品久久久精品久久久| 国产美女午夜福利| 欧美 亚洲 国产 日韩一| 欧美 亚洲 国产 日韩一| 2021少妇久久久久久久久久久| 啦啦啦视频在线资源免费观看| 我的女老师完整版在线观看| 大片电影免费在线观看免费| 777米奇影视久久| 最近中文字幕高清免费大全6| 亚洲欧美清纯卡通| 女的被弄到高潮叫床怎么办| 一级爰片在线观看| 久久久久国产网址| 欧美三级亚洲精品| 一级,二级,三级黄色视频| 大话2 男鬼变身卡| 菩萨蛮人人尽说江南好唐韦庄| 十分钟在线观看高清视频www | 久久综合国产亚洲精品| 99九九在线精品视频 | 肉色欧美久久久久久久蜜桃| 欧美三级亚洲精品| 这个男人来自地球电影免费观看 | 99精国产麻豆久久婷婷| 中文字幕人妻丝袜制服| 成人毛片a级毛片在线播放| 男的添女的下面高潮视频| 五月伊人婷婷丁香| 少妇丰满av| av又黄又爽大尺度在线免费看| 一级毛片 在线播放| 久久影院123| 青春草亚洲视频在线观看| 国产一区二区在线观看av| 国产精品欧美亚洲77777| 久久久精品94久久精品| 性色avwww在线观看| 国产69精品久久久久777片| 99久久精品一区二区三区| 久久女婷五月综合色啪小说| 欧美日本中文国产一区发布| av天堂中文字幕网| 男女边摸边吃奶| 黄片无遮挡物在线观看| 99久久精品热视频| 精品人妻熟女毛片av久久网站| 久热这里只有精品99| 在线播放无遮挡| 精品人妻偷拍中文字幕| 一区二区三区精品91| 97超视频在线观看视频| √禁漫天堂资源中文www| 国模一区二区三区四区视频| 亚洲四区av| 嫩草影院入口| 黄色怎么调成土黄色| 搡老乐熟女国产| 久久久久久久精品精品| 国产免费视频播放在线视频| 曰老女人黄片| 亚洲欧美日韩另类电影网站| 亚洲精品中文字幕在线视频 | 天堂俺去俺来也www色官网| 色视频www国产| 欧美3d第一页| 麻豆精品久久久久久蜜桃| 欧美xxxx性猛交bbbb| 国产日韩欧美视频二区| 日韩一区二区三区影片| 女的被弄到高潮叫床怎么办| 中文字幕久久专区| 国产精品.久久久| 天天躁夜夜躁狠狠久久av| 自拍偷自拍亚洲精品老妇| 亚洲欧美成人综合另类久久久| 一级a做视频免费观看| 啦啦啦中文免费视频观看日本| 亚洲国产最新在线播放| 日韩一区二区三区影片| 亚洲精品成人av观看孕妇| 99热6这里只有精品| 热re99久久精品国产66热6| 嘟嘟电影网在线观看| 亚洲精品中文字幕在线视频 | 久久久久久久久久久丰满| 自线自在国产av| 亚洲图色成人| 亚洲av综合色区一区| 99热这里只有是精品50| 麻豆成人av视频| 看免费成人av毛片| 成人免费观看视频高清| 成年美女黄网站色视频大全免费 | 亚洲欧美一区二区三区黑人 | 久久精品熟女亚洲av麻豆精品| 国产无遮挡羞羞视频在线观看| 亚洲不卡免费看| 夜夜爽夜夜爽视频| 国产精品人妻久久久久久| 国产色婷婷99| 日韩精品免费视频一区二区三区 | 最近中文字幕2019免费版| 3wmmmm亚洲av在线观看| 男女国产视频网站| 男男h啪啪无遮挡| 久久久久精品性色| 如何舔出高潮| 乱码一卡2卡4卡精品| 另类亚洲欧美激情| 精品少妇黑人巨大在线播放| 人妻夜夜爽99麻豆av| 五月伊人婷婷丁香| 久久人人爽人人爽人人片va| 精品少妇久久久久久888优播| 一区二区三区免费毛片| 国产精品国产三级国产专区5o| 国产精品久久久久久精品古装| 国产在线一区二区三区精| 久久国产精品男人的天堂亚洲 | av天堂中文字幕网| 欧美三级亚洲精品| 只有这里有精品99| 免费播放大片免费观看视频在线观看| 高清av免费在线| 久久综合国产亚洲精品| 高清不卡的av网站| 国产美女午夜福利| 在线观看www视频免费| 91精品国产九色| 新久久久久国产一级毛片| 99久久人妻综合| 国产一区有黄有色的免费视频| 欧美变态另类bdsm刘玥| av视频免费观看在线观看| 国产在线免费精品| 精品国产露脸久久av麻豆| 日韩人妻高清精品专区| av免费观看日本| 日日撸夜夜添| 午夜老司机福利剧场| 久久久久久久国产电影| 蜜桃在线观看..| 全区人妻精品视频| 一个人看视频在线观看www免费| 黄色一级大片看看| 国产永久视频网站| 大码成人一级视频| 一区二区三区免费毛片| 嘟嘟电影网在线观看| 国产成人精品一,二区| 精品人妻熟女毛片av久久网站| av网站免费在线观看视频| 最新中文字幕久久久久| 人人澡人人妻人| 人妻少妇偷人精品九色| 老熟女久久久| 国产av码专区亚洲av| 日韩制服骚丝袜av| 精品久久久久久久久av| 美女xxoo啪啪120秒动态图| 99久久综合免费| 秋霞在线观看毛片| 国产精品麻豆人妻色哟哟久久| 少妇高潮的动态图| 免费人妻精品一区二区三区视频| 国产亚洲欧美精品永久| 91精品伊人久久大香线蕉| 成人二区视频| 精品久久久噜噜| 夫妻性生交免费视频一级片| 亚洲精品自拍成人| 亚洲精品乱久久久久久| 人人澡人人妻人| 精品人妻熟女av久视频| 免费人成在线观看视频色| 亚洲成色77777| a级毛片在线看网站| 亚洲伊人久久精品综合| a级毛片免费高清观看在线播放| 久久青草综合色| 亚洲四区av| 欧美日韩亚洲高清精品| 免费黄网站久久成人精品| 精品一品国产午夜福利视频| 伦理电影大哥的女人| 最新的欧美精品一区二区| 亚洲av男天堂| 在线观看美女被高潮喷水网站| a 毛片基地| 精品久久久噜噜| 国产av精品麻豆| 日本vs欧美在线观看视频 | 免费黄网站久久成人精品| 人妻制服诱惑在线中文字幕| 国产午夜精品久久久久久一区二区三区| 日日爽夜夜爽网站| 国产视频首页在线观看| 亚洲精品国产av成人精品| 一级av片app| 国产精品99久久久久久久久| 91午夜精品亚洲一区二区三区| 啦啦啦在线观看免费高清www| 一本久久精品| 熟妇人妻不卡中文字幕| 国产午夜精品一二区理论片| 久久国产乱子免费精品| av国产久精品久网站免费入址| 免费黄频网站在线观看国产| 99视频精品全部免费 在线| 日本vs欧美在线观看视频 | 一级毛片aaaaaa免费看小| 成人免费观看视频高清| 国产 一区精品| 免费不卡的大黄色大毛片视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 人妻制服诱惑在线中文字幕| 国产黄色免费在线视频| 啦啦啦在线观看免费高清www| 一本久久精品| 亚洲第一av免费看| 精品国产一区二区三区久久久樱花| 欧美日韩精品成人综合77777| 久久人人爽av亚洲精品天堂| 午夜激情福利司机影院| 男人添女人高潮全过程视频| 国产av码专区亚洲av| 大码成人一级视频| 欧美日韩国产mv在线观看视频| 另类精品久久| 亚洲丝袜综合中文字幕| 亚洲电影在线观看av| 精品国产一区二区久久| 人人妻人人添人人爽欧美一区卜| 欧美日韩av久久| 亚洲国产毛片av蜜桃av| 成年美女黄网站色视频大全免费 | 夜夜看夜夜爽夜夜摸| 菩萨蛮人人尽说江南好唐韦庄| 观看美女的网站| 欧美3d第一页| 久热这里只有精品99| 亚洲国产日韩一区二区| 九九在线视频观看精品| av免费观看日本| 成人影院久久| 99久国产av精品国产电影| 成年女人在线观看亚洲视频| 成人毛片60女人毛片免费| 人人妻人人澡人人看| 一本大道久久a久久精品| 赤兔流量卡办理| 2022亚洲国产成人精品| 丰满乱子伦码专区| 国产成人aa在线观看| 久久国产精品男人的天堂亚洲 | 亚洲情色 制服丝袜| 一级黄片播放器| 久久久久久久久久久免费av| 桃花免费在线播放| 高清欧美精品videossex| 国产一区有黄有色的免费视频| 久久鲁丝午夜福利片| 精华霜和精华液先用哪个| 人人澡人人妻人| 日日啪夜夜爽| 亚洲av日韩在线播放| 中文精品一卡2卡3卡4更新| 国产美女午夜福利| av有码第一页| 免费久久久久久久精品成人欧美视频 | 日韩精品有码人妻一区| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲高清精品| 久久久久久久久久久免费av| 亚洲国产精品专区欧美| 国语对白做爰xxxⅹ性视频网站| 国产色爽女视频免费观看| 中文字幕人妻丝袜制服| 一级毛片久久久久久久久女| 免费看日本二区| 不卡视频在线观看欧美| 自拍偷自拍亚洲精品老妇| 久久久久久久久大av| 黄色毛片三级朝国网站 | 免费av不卡在线播放| 日本av免费视频播放| 一级二级三级毛片免费看| 国产乱来视频区| 国产又色又爽无遮挡免| 亚洲精品自拍成人| 国产免费视频播放在线视频| 老女人水多毛片| 亚洲精品自拍成人| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 久久久亚洲精品成人影院| 天美传媒精品一区二区| 蜜桃久久精品国产亚洲av| 波野结衣二区三区在线| 色吧在线观看| 国语对白做爰xxxⅹ性视频网站| 中文字幕免费在线视频6| a级片在线免费高清观看视频| 国产免费一级a男人的天堂| 一级a做视频免费观看| 精品少妇内射三级| 国产精品国产三级国产专区5o| 99re6热这里在线精品视频| 五月玫瑰六月丁香| 91精品国产九色| 插逼视频在线观看| 久久久精品94久久精品| 午夜精品国产一区二区电影| 国产精品一区二区三区四区免费观看| 欧美日韩在线观看h| 中国国产av一级| 性色av一级| 国产伦理片在线播放av一区| 精品一区二区免费观看| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 爱豆传媒免费全集在线观看| 啦啦啦啦在线视频资源| 高清午夜精品一区二区三区| 国产视频首页在线观看| 又爽又黄a免费视频| 人妻夜夜爽99麻豆av| 精品少妇久久久久久888优播| 观看免费一级毛片| av播播在线观看一区| 国产精品熟女久久久久浪| 亚洲欧美成人精品一区二区| 最黄视频免费看| 男男h啪啪无遮挡| 一个人免费看片子| freevideosex欧美| 最近手机中文字幕大全| 秋霞伦理黄片| 国产欧美日韩综合在线一区二区 | 夜夜骑夜夜射夜夜干| 免费少妇av软件| 国产精品福利在线免费观看| 91精品国产国语对白视频| 久久久精品免费免费高清| 国产精品成人在线| 99久久精品国产国产毛片| 观看免费一级毛片| 国产真实伦视频高清在线观看| 亚洲国产色片| av有码第一页| 午夜激情久久久久久久| 免费看av在线观看网站| 亚洲经典国产精华液单| 少妇的逼好多水| 国产欧美日韩一区二区三区在线 | 国产成人午夜福利电影在线观看| 日韩人妻高清精品专区| 美女xxoo啪啪120秒动态图| 少妇 在线观看| 人妻 亚洲 视频| 欧美bdsm另类| 国产片特级美女逼逼视频| 午夜免费男女啪啪视频观看| 美女主播在线视频| 王馨瑶露胸无遮挡在线观看| 99热这里只有是精品50| 在线观看免费日韩欧美大片 | 亚洲欧美一区二区三区黑人 | 国产男人的电影天堂91| 能在线免费看毛片的网站| 国产毛片在线视频| 人人妻人人澡人人看| 国模一区二区三区四区视频| 中文字幕久久专区| 日韩av不卡免费在线播放| 一级毛片久久久久久久久女| 免费黄频网站在线观看国产| av免费观看日本| 99久久精品国产国产毛片| 日本午夜av视频| 国产精品一区二区性色av| 丰满迷人的少妇在线观看| 免费黄网站久久成人精品| 欧美精品一区二区免费开放| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载| 一区二区三区四区激情视频| 大又大粗又爽又黄少妇毛片口| 观看av在线不卡| 国产中年淑女户外野战色| 国产av码专区亚洲av| av天堂久久9| 日本91视频免费播放| 成人黄色视频免费在线看| 狂野欧美激情性bbbbbb| 国产成人精品福利久久| av又黄又爽大尺度在线免费看| 国产亚洲一区二区精品| 两个人的视频大全免费| 亚洲高清免费不卡视频| 80岁老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| 国产美女午夜福利| 国产探花极品一区二区| 成人国产av品久久久| 国内少妇人妻偷人精品xxx网站| 人人澡人人妻人| videos熟女内射| 成人毛片a级毛片在线播放| 久久午夜综合久久蜜桃| 晚上一个人看的免费电影| 欧美最新免费一区二区三区| 又粗又硬又长又爽又黄的视频| 国产视频内射| 久久女婷五月综合色啪小说| 国产亚洲午夜精品一区二区久久| 18禁动态无遮挡网站| 精品国产露脸久久av麻豆| 久久久国产一区二区| 国产无遮挡羞羞视频在线观看| 黑丝袜美女国产一区| 成人特级av手机在线观看| 亚洲美女搞黄在线观看| 国产永久视频网站| 久久青草综合色| 黄色欧美视频在线观看| 欧美日韩视频精品一区| a 毛片基地| .国产精品久久| 亚洲欧洲精品一区二区精品久久久 | 国产一级毛片在线| 欧美bdsm另类| 国产精品一二三区在线看| 在线观看av片永久免费下载| 国产淫片久久久久久久久| 国产成人精品福利久久| 青春草亚洲视频在线观看| 80岁老熟妇乱子伦牲交| 97超视频在线观看视频| 国产欧美日韩一区二区三区在线 | 亚洲精品久久久久久婷婷小说| 久久国内精品自在自线图片| 久久久久国产网址| 久久久久久久精品精品| 午夜福利,免费看| 99九九线精品视频在线观看视频| 久热这里只有精品99| 黑人高潮一二区| 久久久久久久久久久丰满| av在线观看视频网站免费| 国产日韩欧美视频二区| 亚洲真实伦在线观看| 国产视频内射| 国产极品天堂在线| 亚洲欧美日韩另类电影网站| 久久久久久久久久久久大奶| 国产精品久久久久久av不卡| 午夜福利,免费看| 一级二级三级毛片免费看| 亚洲无线观看免费| 综合色丁香网| 大香蕉久久网| 在线观看免费高清a一片| 伦理电影大哥的女人| 一级毛片 在线播放| 国产亚洲一区二区精品| 免费高清在线观看视频在线观看| 韩国高清视频一区二区三区| freevideosex欧美| a 毛片基地| 国产精品国产三级专区第一集| 在线免费观看不下载黄p国产| 免费高清在线观看视频在线观看| 在线免费观看不下载黄p国产| 久久久久久久久久久久大奶| 国产精品久久久久成人av| 亚洲欧美日韩东京热| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| www.色视频.com| 国产爽快片一区二区三区| 国产极品粉嫩免费观看在线 | 久久 成人 亚洲| 国产白丝娇喘喷水9色精品| 我要看黄色一级片免费的| 男的添女的下面高潮视频| 水蜜桃什么品种好| 午夜福利在线观看免费完整高清在| 日韩强制内射视频| 秋霞伦理黄片| 有码 亚洲区| 欧美高清成人免费视频www| 亚洲欧洲日产国产| 亚洲丝袜综合中文字幕| 我的女老师完整版在线观看| 精品亚洲乱码少妇综合久久| 七月丁香在线播放| tube8黄色片| 女人精品久久久久毛片| 日日啪夜夜爽| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 18禁裸乳无遮挡动漫免费视频| 在线观看人妻少妇| 日韩大片免费观看网站| 国产日韩欧美亚洲二区| 桃花免费在线播放| 亚洲精品456在线播放app| 精品视频人人做人人爽| 久久99一区二区三区| 亚洲精品日本国产第一区| 国产男女超爽视频在线观看| 欧美日韩国产mv在线观看视频| 国产精品麻豆人妻色哟哟久久| 自拍偷自拍亚洲精品老妇| 51国产日韩欧美| 国产精品国产三级专区第一集| 欧美精品亚洲一区二区| 免费观看av网站的网址| 嫩草影院入口| 精品国产露脸久久av麻豆| 大片免费播放器 马上看| 国产男女超爽视频在线观看| 大片免费播放器 马上看| 国产高清不卡午夜福利| 高清不卡的av网站| 久久鲁丝午夜福利片| 人人妻人人添人人爽欧美一区卜| 日日摸夜夜添夜夜添av毛片| 九色成人免费人妻av| 国产色婷婷99| 一边亲一边摸免费视频| 久久久久久久亚洲中文字幕| 亚洲一级一片aⅴ在线观看| 又爽又黄a免费视频| 精品国产一区二区久久| 男人狂女人下面高潮的视频| 国产一区二区三区av在线| 国产日韩欧美视频二区| 国产永久视频网站| 男人和女人高潮做爰伦理| 国产熟女欧美一区二区| 免费观看无遮挡的男女|