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

    雷達在昆蟲學研究中的應用

    2011-08-27 03:59:42封洪強
    植物保護 2011年5期
    關鍵詞:昆蟲學定向昆蟲

    封洪強

    (河南省農(nóng)業(yè)科學院植物保護研究所,河南省農(nóng)作物病蟲害防治重點實驗室,鄭州 450002)

    雷達在昆蟲學研究中的應用

    封洪強*

    (河南省農(nóng)業(yè)科學院植物保護研究所,河南省農(nóng)作物病蟲害防治重點實驗室,鄭州 450002)

    昆蟲雷達因能無干擾地遠距離探測到自然狀態(tài)下昆蟲自由飛行行為,已被廣泛用于昆蟲高空遠距離遷飛和近地飛行行為研究,成為昆蟲學研究中一個卓越的研究工具。同時天氣雷達也越來越多地用于鳥和昆蟲的監(jiān)測,探地雷達也在探測堤壩蟻穴中起到了重要作用。本文分析了各種不同類型雷達的性能和近年來在昆蟲學研究中的應用,以昆蟲學家的視角審視了雷達能為昆蟲學研究做什么以及其中的挑戰(zhàn)。

    雷達; 昆蟲; 遷飛; 目標識別; 應用

    雷達(radar)是無線電檢測和測距(radio detection and ranging)的縮寫,是利用電磁波進行探測的電子設備,它發(fā)射電磁波對目標進行照射并接收其回波,從而測得目標的距離、距離變化率(徑向速度)、方位、高度及目標形狀等信息[1-3]。自從1949年被證明可以檢測到單個昆蟲以來,雷達不斷地在昆蟲學研究中得到應用,并逐漸演變出一系列專門設計的昆蟲雷達,國際上形成了雷達昆蟲學學科,還建立了雷達昆蟲學網(wǎng)站(http:∥www.pems.adfa.edu.au/~s9104004/trews/ww_re_hp.htm)。

    關于昆蟲雷達的原理、類型、發(fā)展歷史和對昆蟲行為研究的經(jīng)典貢獻已有較為詳細的總結[1,4-8],在此不再贅述。本文主要分析最近十年來用于昆蟲學研究的雷達的性能、對昆蟲學研究的貢獻以及面臨的挑戰(zhàn),并提出對我國雷達昆蟲學發(fā)展的建議,激發(fā)更多昆蟲學家的興趣。

    1 天氣雷達在昆蟲學研究中的應用

    1.1 天氣雷達鑒定昆蟲目標的意義

    敏感的微波多普勒天氣雷達經(jīng)常觀測到邊界層的晴空回波(clear-air echoes),晴空回波的產(chǎn)生原因自1940年代雷達應用之時開始就一直存在爭論[9]。晴空回波主要來自于兩類散射體或其混合體:(1)產(chǎn)生布拉格散射的折射率的梯度變化(折射率波動)或(2)產(chǎn)生粒子散射的空中粒子,包括空中生物(如昆蟲、鳥和蝙蝠等)、灰塵與其他粒子[9-12]。

    布拉格散射是大氣中折射率波動產(chǎn)生的散射,主要來自于濕度在雷達的半波長尺度上的變化或更小程度上溫度的變化[9],這是頻率在3GHz以下的雷達的主要目標。而鳥和昆蟲等粒子散射體是更高頻率的雷達的主要目標[11]。這兩種散射對正在運行的天氣雷達都普遍存在,如美國的3GHz的多普勒天氣預警雷達(weather-surveillance radar-1988,doppler,WSR-88D)[11]。通常飛行能力弱的非遷飛小型昆蟲在X波段雷達上的橫截面約為10-10~10-7m2,而大型昆蟲為10-5m2以上[11],因此典型的天氣雷達因其分辨率空間(resolution volume)大而觀測不到單個的昆蟲,但當昆蟲的密度足夠大時作為體積散射目標確實可以被這些天氣雷達“看”到[11]。

    在陸地上空從春到秋(甚至在冬天氣溫超過10℃時)昆蟲回波經(jīng)常出現(xiàn)在晴朗的天空,它們可以被看做不期望的雜波,也可以為被看做是需要的(盡管不是完美的)晴空示蹤者來測風[11,13]。使用昆蟲回波測風的可靠性取決于:(1)用一個正確的方法將昆蟲回波標記出來,從而最大限度地得到好的數(shù)據(jù)和減少雜波干擾;(2)如何估計測風的誤差,包括由昆蟲自主飛行產(chǎn)生的測風偏差[14]。氣象雷達學家近年來的研究主要集中在如何識別空中生物回波以及昆蟲回波是否可以用來測風上[13-14]。

    1.2 天氣雷達鑒定空中生物的原理

    一部多普勒測云雷達通常測量雷達基本反射率(base reflectivity)、徑向速度(radial velocity)和徑向速度的多普勒譜寬(width of spectrum)[15-16],如美國的 WSR-88D雷達網(wǎng)的二級數(shù)據(jù)[16]。用這些信息可以探測到空中的生物,如鳥、蝙蝠和昆蟲。

    徑向速度是用于測量目標朝向或背離雷達的運動分量[16],是根據(jù)多普勒效應計算得來的。所謂多普勒效應就是當聲音、光和無線電波等振動源與觀測者以相對速度V運動時,觀測者所收到的振動頻率與振動源所發(fā)出的頻率有所不同。由多普勒效應所形成的頻率變化叫做多普勒頻移,它與相對速度V成正比,與振動的頻率成反比[2]。速度方位顯示(velocity azimuth display,VAD)是(假設風場是線性的)對每一個距離檔上各個高度的徑向速度進行正弦曲線擬合得到的[13]。利用VAD與當?shù)氐娘L速進行比較可以鑒別回波是否為昆蟲[13]。使用該方法測到1.2~1.4km高度上的目標為昆蟲,它們的速度為18m/s,比同高度上的雨滴快3m/s[17]。

    多普勒譜寬用于測量脈量內目標的多普勒速度的變異,如果一個脈量內所有的目標都以相同徑向速度飛行,譜寬將達到最小值,如果一個脈量內有多種生物目標以不同的徑向速度飛行,譜寬就會很大[15-16]。大多數(shù)昆蟲的空中速率都較小而鳥的飛行速率較大,因而大氣中以昆蟲為主時譜寬較小,以鳥為主時譜寬較大[16]。

    雙極化(dual-polarization)雷達發(fā)射和接收水平(H)和垂直(V)極化的電磁波,這使得我們可以估計雷達目標的幾何形狀[18-20],因此雙極化雷達比單極化雷達能更有效地鑒定目標,而且可確定昆蟲是否具有定向行為以及定向是否與風向一致[13-14,21]。用于區(qū)分生物目標的雙極化變量包括差反射率因子(differential reflectivity)ZDR、微 分 相 位 (differential phase)δ、相互垂直的極化回波的相關系數(shù)[16,20,22]。

    差反射率因子ZDR指水平和垂直頻道的反射率比值并以對數(shù)表示,即ZDR=10log(Zh/Zv)(dB),其中Zh和Zv分別指水平和垂直極化的反射率因子[9-10]。對于M次脈沖每一分辨率空間ρhv就有M個譜系數(shù)的測量值[22]。令κ為一系列取值從1到M的譜系數(shù)的序號,κ可以轉化為相對應的徑向速度,然后,Sh(κ)和Sv(κ)分別為水平和垂直頻道的包含了信號與噪音在內的復雜譜系數(shù),那么微分相位δ(κ)=arg(〈Sh(κ)Sv(κ)〉)[20],其中,雙極化雷達ρhv可由水平和垂直頻道的復雜譜系數(shù)的三點平均求得[18,22]。

    由于昆蟲和鳥是似橢圓體的目標,鳥的寬與長之比為1∶2~1∶3,昆蟲的寬與長之比為1∶3~1∶10[9],垂直極化將給出一個各向同性的掃描,而水平極化掃描昆蟲側面時得到較強的信號,因此鳥和昆蟲的ZDR值比水滴小。一般來說,降雨的ZDR值為0~2dB[10],鳥的ZDR值為-1~3dB[20],昆蟲回波的ZDR值為5~10dB[9-10]。當昆蟲具有共同定向時,ZDR值可用于將昆蟲從其他目標中區(qū)分出來[13],而且比較ZDR值還可確定昆蟲是否在進行遷飛[9]。

    使用10cm波長的垂直分辨率為2.5m的高分辨率S波段調頻連續(xù)波雷達(FMCW)可以將小型昆蟲(<10mm)從布拉格散射體中分離出來[11]。這次FMCW觀測結果與毫米波雷達觀測到的結果[23-24]一致,但觀測到的平均反射率因子遠小于其他觀測[9-10],可能是由于其他觀測正好趕上昆蟲密度特別大的時期[11]。

    毫米波雷達(3mm和8.6mm或94GHz和35GHz)由于時間和空間分辨率高,使氣象學家可以探測到更小的冰晶來研究大氣的微觀結構[12],也同樣有利于觀測微小昆蟲。Luke等人用35GHz毫米波測云雷達(MMCRs)輔以氣球和風箏攜帶的遙控捕蟲裝置證明從地面至700m高度充滿了小型昆蟲(體長約2mm,翅長4~5mm,翅寬1~2mm),發(fā)明出一種以35GHz毫米波測云雷達MMCR記錄的多普勒譜、描述昆蟲譜的特征提取機和神經(jīng)網(wǎng)絡算法提取小型昆蟲回波的技術[12]。只用 MMCR多普勒譜作為輸入時這項技術能鑒定出92%以上的昆蟲回波。

    1.3 使用天氣雷達區(qū)分鳥與昆蟲

    大量的文獻表明鳥和昆蟲都是晴空回波的主要來源,氣象學家區(qū)分鳥與昆蟲是由于他們認為遷飛的鳥會對測風造成嚴重的偏差,而認為昆蟲造成的偏差較?。?5]??墒菦]有一種可靠的簡單方法把昆蟲與鳥區(qū)分開,鳥類學家用雷達時經(jīng)常假設他們的目標是鳥,而昆蟲學家則經(jīng)常假設為昆蟲[25]。使用天氣雷達在一定程度上可以區(qū)分鳥與昆蟲,要達到百分之百的鑒定昆蟲目標是不可能,因為昆蟲回波的所有參數(shù)都與其他類型目標的參數(shù)有重疊[13]。

    每一分辨率空間可以用極化雷達的一系列極化參數(shù)進行描述,根據(jù)這一空間內目標的種類建立各種參數(shù)的描述函數(shù),可以對不同目標進行分類[18]。根據(jù)反射率大小似乎可以區(qū)分鳥與昆蟲,但由于昆蟲在空中的密度會非常大,高密度的昆蟲可以產(chǎn)生很高的反射率,也就是說較低的反射率可以排除鳥的可能,但高的反射率卻不能排除昆蟲的可能性,因此單以反射率不能區(qū)分鳥與昆蟲[25]。

    在早期的利用極化雷達鑒定空中目標的研究中發(fā)現(xiàn),ZDR值可用于區(qū)分鳥與昆蟲[20,26],但是 Mueller認為鳥的ZDR值高于昆蟲,而后來的研究卻認為恰恰相反,即昆蟲的ZDR值大于鳥[18,20,22]。這是由于昆蟲和鳥的體型都有很多變化,他們的ZDR值的變化范圍相互重疊,很難說哪一類生物的體型更長[18,25]。同樣原因δ值也無法將昆蟲與鳥區(qū)分開[16]。如果兩種或更多種散射體同時存在于同一個分辨率空間,這個空間的極化變量取決于混合的比例、它們的排列、生物目標的定向、在這一空間內的位置(均勻填充、成層、成團或其他)。因此當多種類型的目標同時存在于一個分辨率空間時,極化變量的變異范圍就會變大,假設一種目標為主體得到的模糊分類標準就難以奏效[16,18]。

    最近建立的雙極化雷達的速度方位顯示譜圖(SVAD)與三維(方位、多普勒頻移和極化譜密度)方 位 譜 圖 (3-D spectral azimuthal histogram,3DASH)為區(qū)分目標種類提供了新的視覺分析技術[18,22]。SVAD是在VAD圖基礎上以不同灰度表示出多普勒系數(shù)以及極化變量(ZDR、δ、ρhv)的譜密度圖,這種譜圖可以在同一取樣空間內同時存在鳥與昆蟲時將兩者的回波區(qū)分開[22]。3DASH是一定距離檔上的方位—速度—δ分布的三維圖,并用灰度表示發(fā)生概率。它可以投影成方位—速度、方位—δ和速度—δ圖,通過設置不同的透明度可以除去干擾看清目標[18]。3DASH的螺旋結構印跡代表了不同目標的特性,甚至有望區(qū)分不同類型的昆蟲[18]。

    不同波段的雷達同時觀測,根據(jù)目標在不同波段的雷達上的反射率差異可以確定目標的種類并容易將昆蟲與鳥區(qū)分開。例如DOW3(3.198cm)、WSR-88D(10.0~11.1cm)和 UMASS(0.315 7cm)的聯(lián)合使用[25],S波段(10.7cm)、X 波段(3.2cm)和 C波段(5.4~5.5cm)雷達的聯(lián)合使用[9]。激光雷達的波長比雨滴和昆蟲小,對由高密度的微小雨滴(10μm)組成的云非常敏感,但一般探測不到低密度的昆蟲。因此英國人曾利用兩部毫米波(35GHz和94GHz)垂直雷達、一部厘米波垂直昆蟲雷達與激光雷達同時使用來區(qū)分昆蟲目標[15]。

    1.4 天氣雷達觀測到的空中昆蟲行為

    昆蟲,尤其是自主飛行能力很弱的微小昆蟲(通常體重為1mg[15]或體長小于10mm[12])白天借助由于地面受熱產(chǎn)生的上升氣流升空,其上限是對流邊界層(convective boundary layer,CBL)頂[13]。晚上地面冷卻形成逆溫層,行星邊界層比較穩(wěn)定,夜間遷飛的昆蟲個體一般較大(>40mg),有較強的自主飛行能力,用來測風誤差較大,因此大多數(shù)氣象雷達對昆蟲的觀測都集中在白天[13]。

    英國的一部毫米波垂直多普勒雷達觀測到小型昆蟲的上升速率可達4m/s,可上升到高出CBL 800m以上的高度[15]。大多數(shù)小型昆蟲(蚜蟲、飛虱)自主飛行的上升速率都小于0.2m/s,即使是大型昆蟲(蝗蟲和夜蛾)其上升速率也只有0.4~0.5m/s,那么觀測到的4m/s的上升速率肯定是昆蟲與上升氣流的速率之和[15]。飛行高度超出CBL表明昆蟲是一種主動飛行者而不是單純靠上升氣流飛行的[15]。

    還有研究表明盡管循環(huán)氣流可以將昆蟲帶到更高的地方,但昆蟲層頂卻在1 800mAGL,那兒的溫度為10~15℃,低于這個溫度蝗蟲是不飛行的[10]。如果高于這個高度,由于昆蟲對氣溫的反應,會通過折疊翅膀停止飛行甚至向下飛行主動下降使自己保持在細線(氣象雷達觀測到的昆蟲帶或層)內[10]。機載W波段(3mm)垂直雷達的觀測結果卻表明使CBL內昆蟲產(chǎn)生反應的是突然的垂直運動而不是溫度[24]。

    氣象雷達曾多次觀測到大量遷飛的昆蟲具有共同定向特性,而且頭向與風向不一致[10,13,19]。有趣的是英國用多普勒雷達觀測到昆蟲的定向在6月份更偏北一些,而在8月份更偏南一些,這與遷飛的季節(jié)性偏好方向一致[13]。

    2 昆蟲雷達在昆蟲遷飛研究中的應用

    由于用傳統(tǒng)的昆蟲學方法很難觀測高空中遷飛的昆蟲,20世紀60年代末開始使用專門設計建造的昆蟲雷達,這極大地豐富了我們對昆蟲遷飛行為的認識[1,5,7-8,27]。目前各國普遍采用的昆蟲監(jiān)測雷達(insect monitoring radar,IMR)為能長期自動化運行的旋轉極化的垂直昆蟲雷達(vertical looking radar,VLR)[28-31]。

    2.1 旋轉極化的垂直昆蟲雷達

    昆蟲監(jiān)測雷達或垂直昆蟲雷達(IMR或VLR)的天線不像傳統(tǒng)雷達那樣以各個仰角進行360°掃描并將結果在平面顯示器(plan position indicator,PPI)上顯示出來,而是產(chǎn)生一個垂直的波束,當昆蟲飛過波束時對昆蟲進行測量[28]。實際上1949年Crawford首次證明雷達能探測到昆蟲回波時就是這樣做的,這種新型昆蟲雷達的原理和技術在20世紀80年代初就確定了下來,90年代的“數(shù)字革命”給昆蟲雷達的發(fā)展帶來新生,迅速發(fā)展的廉價數(shù)字計算機以及電子設備使制造低成本的自動運行的昆蟲雷達成為現(xiàn)實[28]。美國的Beerwinkle首先研制了一個簡單的垂直雷達并無人值守地運行了兩年,英國和澳大利亞也緊跟著研制出更為先進的垂直雷達,并從上世紀90年代末開始了長期的自動觀測[28,32]。

    旋轉極化的垂直昆蟲雷達發(fā)射一個環(huán)形對稱垂直向上看的波束,并通過機械轉動波導饋電使波束內的線性極化平面旋轉。饋電與垂直方向有一個非常小的0.2°的偏角,這使得波束產(chǎn)生輕微的擺動。于是整個雷達波束就像陀螺旋轉變慢時一樣產(chǎn)生旋轉的錐形掃描[29]。由于饋電的垂直偏角非常小,能以300r/min的轉速高速轉動,因此當昆蟲穿越垂直的雷達波束時可以完成幾個完整的旋轉掃描[28]。由于波束的位置和極化在任何時候都是已知的,接收到的回波信號包含了目標的軌跡(相對于垂直方向的位置)和由極化平面旋轉帶來的背向散射截面的變化兩方面的信息[28]。因此這種掃描的效果相當于一個目標飛過頭頂時以不同角度去反復“切割”它,得到昆蟲不同截面的信息。

    這種昆蟲雷達的旋轉極化配置(多個“切面”)與天氣雷達的雙極化配置(兩個“切面”)相比,能采集到目標更豐富的體形信息,大大增強了對昆蟲的識別能力。這種雷達檢測到的回波能量可以用一個非常復雜的公式表示,包含了目標的水平速率、位移方向、身體的軸向、到波束中心線的最近距離、昆蟲的大小和兩個代表昆蟲體形的參數(shù)在內的7個參數(shù)[28-29]。由于這種雷達每一圈掃描都可以得到64個回波能量數(shù)據(jù),這些數(shù)據(jù)足夠用于擬合7個參數(shù),擬合的回歸系數(shù)可達0.9以上[28-29]。具體的數(shù)據(jù)擬合方法十分復雜[33-35],但隨著計算機技術的突飛猛進,能在現(xiàn)代微型計算機上運行的C++寫成的程序可以在4min內自動分析5萬個昆蟲信號[35]。將模擬的昆蟲回波信號與實測的昆蟲回波信號相減得到的殘差并不為0,也不是隨機分布的,而是出現(xiàn)有規(guī)律的大小相似的波動,每一個波峰就代表了一次振翅,根據(jù)這一原理可以將昆蟲的振翅頻率提取出來[36]。

    1.3觀察指標 根據(jù)線性視覺模擬評分標尺對患者舒適程度進行評估[3],分為無癥狀、輕度、中度、重度不舒適,分別為0分、1~4分、5~7分、8~10分。自行設計調查表,對患者進行護理滿意率調查。

    由于垂直雷達是垂直向上照射的,假設昆蟲在飛過雷達波束時不主動(或因氣流顛簸被動)地改變傾角、飛行偏角和滾動傾斜角,所有的昆蟲都以相同的姿勢飛過電磁波束,因此提取出的4個關于昆蟲本身特性的參數(shù)(昆蟲的個體大小、振翅頻率以及2個體形參數(shù))就可以有效地鑒別昆蟲的種類[37-39]。英國昆蟲學家則使用最大背向散射截面δxx與最小截面δyy(與δxx垂直)的比值來鑒定昆蟲種類[38-39]。較大的比值代表較細長的體形,如15∶1代表草蛉,較小的比值代表較短的體形,如5∶1代表步甲,而1∶1代表瓢蟲[29]。澳大利亞科學家用兩個體形參數(shù)(α2和α4)描述昆蟲的共極的線性極化類型(copolar-linear polarization pattern,CLPP),理論預測與實際觀測的結果一致[37]。用α2和α4描述的CLPP由于有豐富的圖形變化(2個瓣至4個瓣及瓣的不同形態(tài))可能代表了更為豐富的昆蟲體形信息,可用于昆蟲種類的詳細鑒定[37]。但實驗室測定一種昆蟲時得到的CLPP形態(tài)變化也較大,這就給CLPP鑒定種類帶來了困難,還需要進一步的研究。

    現(xiàn)行的IMR或VLR還不能連續(xù)地處理信號,即不能連續(xù)的觀測垂直空間,而是把垂直高度分為若干個距離門(range gate)即高度帶,兩個距離門之間有一定的間隔[29]。澳大利亞的IMR在垂直方向上有8個距離門,每個距離門的高度(即分辨率)為50m,距離門的間隔為150m,由于第一個距離門含有大量地物回波干擾成為無效的距離門,因此其有效的觀測高度為200~1 400m[28,34]。英國 VLR有15個有效的距離門,門高45m,間隔26m,有效觀測高度為150~1 166m[29]。因昆蟲的大小不同實際能探測到的高度(距離)也不同,1mg的昆蟲只能在200m以下被探測到(但由于地物干擾實際上是探測不到),2mg的昆蟲的探測高度可達350m,5mg的可達550m,15mg的可達1 200m[29]。

    隨著現(xiàn)代通訊技術的發(fā)展,遠程遙控的多機聯(lián)網(wǎng)自動觀測、甚至與現(xiàn)行的天氣雷達或雷達網(wǎng)進行聯(lián)合或聯(lián)網(wǎng)觀測也成為現(xiàn)實,監(jiān)測結果可通過Internet向大眾實時發(fā)布[40]。

    2.2 傳統(tǒng)掃描雷達的新生

    旋轉極化的垂直昆蟲雷達雖然在目標鑒定上有著無可比擬的優(yōu)勢,但它的取樣空間小而且200m以下無法觀測。傳統(tǒng)的掃描昆蟲雷達取樣空間大,但信號數(shù)字化和分析困難,這也是國際雷達昆蟲學界轉而開發(fā)垂直昆蟲雷達的主要原因。近年來隨著信息技術的迅速發(fā)展掃描雷達的數(shù)字化取得了長足的進步,使傳統(tǒng)掃描昆蟲雷達重獲新生[1,41-42]。

    傳統(tǒng)掃描昆蟲雷達依靠對雷達屏幕顯示的回波進行照相和攝像來記錄數(shù)據(jù),以供后期人工分析用。澳大利亞的Drake博士詳細描述了手動分析掃描昆蟲雷達數(shù)據(jù)的方法和過程[43]。1998年我國第2部昆蟲雷達(KC-2)建成,為了克服信號采集與分析的困難,在郭予元院士和吳孔明博士的大力支持下程登發(fā)博士開發(fā)出了掃描昆蟲雷達實時數(shù)據(jù)采集、分析系統(tǒng)[41-42]。

    KC-2昆蟲雷達的顯示單元提供了一個8針VGA信號接口,供連接其他計算機顯示器用。使用基于PCI總線的8bit高速度圖像采集卡可直接從雷達主機將PPI采集為分辨率為768×536的8位黑白圖像。使用VC++編程實現(xiàn)了PPI圖像的采集、存貯與再現(xiàn),并根據(jù)Drake的方法實現(xiàn)了昆蟲回波的自動識別和密度統(tǒng)計。該系統(tǒng)還可將連續(xù)采集的3張圖片疊加在一起分別以紅、綠、藍3種不同的顏色顯示回波點,使用戶可以清楚地分辨出地物回波(白色)以及昆蟲的位移軌跡,并能測量出位移速度[1,41-42]。將測到的昆蟲位移速度與同時實地測得的風速進行矢量差運算得到昆蟲的定向[1,44]。

    由于這套系統(tǒng)大大提高了掃描昆蟲雷達的工作效率,使我國昆蟲學家得以進行長期的雷達昆蟲學研究,先后觀測了甜菜夜蛾、草地螟、棉鈴蟲和黏蟲等重要害蟲以及黃蜻和毛婪步甲等天敵在我國華北地區(qū)的遷飛行為,獲得了重要的行為參數(shù)[1,44-54]。這部雷達以及數(shù)據(jù)采集與分析系統(tǒng)直到今天還在正常運行,是標志著我國獨立進行雷達昆蟲學研究的一個里程碑。

    無獨有偶,加拿大科學家最近開發(fā)了一個開源平臺radR,用于采集、存貯和分析航海雷達以監(jiān)測模式觀測到的生物目標[55]。底層數(shù)據(jù)處理是由C代碼實現(xiàn)的,用戶定義的函數(shù)是由R語言寫成的,可以在計算時按預設的步驟調用。輸出數(shù)據(jù)格式設計的可以在將來增加數(shù)據(jù)而不必改變C代碼。目前這個開源平臺支持兩種采集卡,可以在 Windows XP和Linux兩個操作系統(tǒng)下運行。

    2.3 遷飛昆蟲的定向與成層

    昆蟲的共同定向是指高空遷飛的昆蟲頭朝著相同的方向飛行的行為,是迄今為止昆蟲雷達觀測到的最神奇且神秘的行為現(xiàn)象[8,27,29,56]。自從 1968年Scheafer首次用雷達觀測到夜間遷飛昆蟲具有共同定向以來,這種行為在各個地區(qū)和各種夜間遷飛的昆蟲身上反復地被觀測到[27,56-59]。由于在夜間高空飛行,視覺的定向是不可能的,兩種最可能的機制是“風定向”和“羅盤定向”[8,57-59]。

    風定向機制認為昆蟲可以用觸角探測到風向,以便始終向下風方向定向[8],那么當不同高度上的風向發(fā)生轉變時,昆蟲的定向也應跟著轉變。有趣的是一個模擬模型預測的結果是不同高度上風向的轉變(Ekman螺旋)誤導昆蟲在北半球向右偏而在南半球向左偏[60]。英國觀測到的數(shù)據(jù)支持模型預測的結果,即北半球昆蟲順風定向時向右偏[60],但澳大利亞的觀測結果卻不支持模型的預測結果(王??郏瑐€人交流)。同樣處在北半球的中國觀測到的棉鈴蟲秋季回遷中當在下風方向大于225°(西南方向)時向左偏,而當下風方向小于225°時向右偏[47],這同樣也不支持模型的預測結果。如果昆蟲總是順風定向那么遷到高緯度地區(qū)繁殖的昆蟲的后代,很可能由于所謂的Pied Piper效應而滯留在高緯度地區(qū),無法安全回到低緯度地區(qū)越冬,然而近年觀測到昆蟲回 遷 的 事 例 越 來 越 多[44-45,47-48,51,61-62],這 表 明昆蟲不是盲目順風定向、順風遷飛的。

    羅盤定向機制認為昆蟲以某種羅盤識別地理方向,遷飛時總是向著自己偏好的方向定向。近年來越來越多的觀測表明夜間遷飛的昆蟲具有這種非常驚人的羅盤定向能力[44,47-48,51,62-64]。我國的草地螟春季遷飛時向東北定向,黏蟲秋季回遷時向西南定向,當風力較弱時這種主動定向可以對遷飛軌跡產(chǎn)生較大影響,但當風力較強時這種校正起到的效果甚微[48-49]。除了黃蜻飛行能力較強可以側頂風飛行外[51],總體上昆蟲的飛行能力較弱,選擇與下風方向呈銳角順風飛行。但這并非是被動順風飛行,棉鈴蟲秋季總是向著西南方向偏轉,而且選擇有適合的風向出現(xiàn)的夜晚或高度飛行[45,47]。昆蟲的定向是有季節(jié)適應性的,即在北半球春季向北而秋季向南[62],它的生態(tài)學意義是使昆蟲能更快地到達適生區(qū)繁殖或安全越冬[61]。

    夜間遷飛昆蟲的成層也是使用昆蟲雷達經(jīng)常觀測到的一個現(xiàn)象[8]。在高緯度地區(qū)的觀測發(fā)現(xiàn)昆蟲的成層與逆溫有關,由于昆蟲有最適飛行溫度,在溫度較低的環(huán)境里經(jīng)常聚集在逆溫層頂飛行[49,51,65-67]。當溫度較高時夜間遷飛的昆蟲選擇風速最大的高度成層,這使昆蟲以最小的能量實現(xiàn)最大位移[45,47,61-62,64],我國還發(fā)現(xiàn)昆蟲的成層與風向和風切變有關[45,47]。因此,昆蟲的成層可能與溫度、風向、風速和風切變等多個因素有關。

    3 諧波雷達與近地昆蟲行為研究

    3.1 諧波雷達原理

    昆蟲諧波雷達的原理是發(fā)射高能的微波照射粘在單個昆蟲身上的轉發(fā)器(transponder)并接收由轉發(fā)器轉發(fā)回來的二次調諧波來探測到昆蟲目標[7,68]。由于植物、地物和海洋不產(chǎn)生高次諧波而由它們產(chǎn)生的基波又被抑制,故諧波雷達具有很強的抗雜波能力,能有效地排除地物雜波的干擾[69]。根據(jù)這一原理現(xiàn)在已建成了用于研究昆蟲爬行行為的便攜式諧波雷達和用于研究昆蟲近地飛行行為的基站式諧波雷達[69]。

    使用諧波雷達需要在昆蟲身上粘一個轉發(fā)器,類似于標記回收試驗所用的標記物,也稱標簽(tag)。這種標簽由一個肖特基二極管和一段天線組成[68,70],肖特基二極管能夠吸收高頻電磁波的部分能量并反射出二次諧波,不需自帶電源,故可做得很小很輕[7,69],不會影響昆蟲本身的行為[71-73]。

    所有使用諧波雷達的昆蟲學家都要自制標簽,在標簽的探測性能與耐久性之間找到一個平衡。自制或修改標簽包括改變二極管、天線、二極管—天線配置和天線配置(直線形、圓形和卵圓形)[74]。天線長度理論上是雷達波長的一半,提高雷達的頻率可以有效地縮短標簽的長度。利用經(jīng)濟的航海雷達主機制造的昆蟲雷達系統(tǒng)由于使用的頻率是9.4GHz,因此可以使用重量小于3mg,長度約十幾毫米的標簽[7,68,73,75]。到目前為止最輕的標簽(0.4mg)是由Roland等制作的,盡管這種標簽在野外的壽命是非常有限的[69]。

    雖然標簽越小對昆蟲行為的影響越小,但使用較小的標簽就要提高雷達的頻率,頻率越高的電磁波穿透力就越弱,因此要根據(jù)昆蟲的行為特點及所要研究的特性來決定使用什么樣的諧波雷達和標簽。

    3.2 便攜式諧波雷達

    常用的便攜式RECCO收發(fā)系統(tǒng),重1.6kg,是尋找雪崩遇難者的雷達系統(tǒng)(www.recco.com),使用917MHz的初始電磁波,可穿透作物、雪和土壤,使用RECCO配發(fā)的標簽可以探測到土壤下1m處的目標[70,74]。諧波雷達的探測性能與標簽有關,使用O’Neal等的16cm標簽,RECCO只能探測到土下30cm處的目標[74]。由于商用的RECCO收發(fā)系統(tǒng)無法改變電磁波頻率,這就限制了使用更小的轉發(fā)器。如果研究的對象個體較大,運動速度不快且有隱藏于土壤或茂密的作物叢中的習性如步甲,使用商用的RECCO收發(fā)系統(tǒng)和十幾厘米長的標簽能得到令人滿意的結果[74,76-79]。2004年 Williams等人還在北京研究了光肩星天牛的飛行行為取得成功,其中最遠移動距離超過了90m[70]。

    為了研究馬鈴薯甲蟲的運動行為,加拿大的Colpitts研究小組研制了頻率為9.4GHz的便攜式諧波昆蟲雷達和相應的重3mg長16mm的電子標簽,并詳細描述了標簽的制作和標定方法,這種系統(tǒng)還可用于研究與馬鈴薯甲蟲相似的多種有益昆蟲的行為[68,75,80]。使用這套系統(tǒng)在馬鈴薯田探測113頭被電子標簽標記過的馬鈴薯甲蟲的發(fā)現(xiàn)率為87%,平均每頭用時(165±15)s,在距離蘋果樹3m處能發(fā)現(xiàn)78%的被標記的李象鼻蟲,當被標記的玉米根螢葉甲在正對雷達的玉米葉、雄穗和莖稈上時,在6~20m距離處的發(fā)現(xiàn)率可達80%~100%,但當玉米螢葉甲在玉米雌穗或莖稈背后時幾乎不能被發(fā)現(xiàn),在第二行玉米上的葉甲的發(fā)現(xiàn)率也明顯下降,但如果使用者可以繞過作物進行掃描,發(fā)現(xiàn)率會非常高[81]。

    昆蟲行為是否會因蟲體粘上了電子標簽而受到影響一直是廣受質疑的。Boiteau等的研究發(fā)現(xiàn),被標記的3種甲蟲的確會因增加了重量負擔或標簽的纏繞爬行速度受到一定的影響,影響因種類而異。但這種負面影響遠比環(huán)境本身對昆蟲行為的影響小,而且這輕微的影響與使用此技術獲得的昆蟲的重要生態(tài)學信息相比是微不足道的,便攜式諧波雷達不失為一種優(yōu)良的研究手段[72,82]。桂連友最近使用這種技術研究了饑餓馬鈴薯甲蟲的行為,發(fā)現(xiàn)饑餓能增加甲蟲的運動距離而且這一影響比試驗室研究的結果更強[83]。他從加拿大學習回國后也建造了相同的便攜式諧波昆蟲雷達,正在開展橘大實蠅行為的研究[84]。

    3.3 基站式諧波雷達

    便攜式諧波雷達從某種意義上講是一種“方向發(fā)現(xiàn)”裝置,即發(fā)現(xiàn)標記昆蟲的方向,它不能很快地繪出昆蟲飛行的軌跡,因為它的探測距離一般小于50m,操作者必須能夠跟得上昆蟲的爬行或飛行才行[7]。1995年英國的Riley領導的雷達昆蟲學研究團隊研制了世界首部真正的諧波昆蟲雷達。這部基站式諧波昆蟲雷達有兩個碟形的天線,下部的一個直徑1.5m的拋物面天線用于發(fā)射功率25kW、脈沖寬度0.1μs、波長3.2cm的電磁波,上部的第二個直徑0.7m的拋物面小天線用于接收反射回來的1.6cm波長的諧波,兩個天線都以20r/min的轉速旋轉,可以實時繪制出半徑約900m范圍內飛行高度在3m以下的標記昆蟲的飛行軌跡,空間解析度為5m[8,73,85]。

    由于社會性昆蟲蜜蜂(或其他的蜂)有“家”的概念,每次外出都要重新回家,用電子標簽標記蜜蜂一般不會丟失,同時蜜蜂的行為也充滿了不解之謎,所以成為基站式諧波昆蟲雷達研究的主要對象[8]。研究發(fā)現(xiàn),初次出巢的工蜂通過反復飛行訓練逐漸增加離巢距離從而保持回巢的能力,這種距離的增加是由于飛行速度的增加而不是飛行時間的增加,這種飛行行為的漸變與出巢的次數(shù)有關而與工蜂的年齡無關。熊蜂和蜜蜂都能正確地從覓食處徑直回到蜂巢,即使正在覓食的工蜂被人為移動了位置也會自動按既定的方向和路線飛行,飛行過程中可以克服風的飄移,當蜜蜂準確定位蜂巢或食物源時進行隨機尋找。研究還證明了von Frisch的蜜蜂舞蹈理論的正確性,即舞蹈交流包含了方向和距離兩方面的信息[8]。由于這種諧波雷達的出色表現(xiàn),最近德國也請洛桑研究院的雷達工程師為他們建造了世界第二部基站式諧波雷達,用于蜜蜂行為的研究,明確了蜜蜂使用一系列地面標志物作為測量距離的標尺[86]。

    盡管許多黃地老虎飛出雷達的視野,用于研究幾百米范圍內雄蛾尋找性信息素源的行為也取得了成功,發(fā)現(xiàn)當氣味濃度足夠大時,順風飛行的雄蛾轉為側風飛行而側風飛行的則轉為逆風飛行[87-88]。首次用諧波雷達研究5種蝴蝶近地面飛行行為發(fā)現(xiàn)了快速的直線型和慢速的包括循環(huán)覓食的非直線型的兩種飛行軌跡,展示了用于蝴蝶飛行研究的潛力[89]。慶網(wǎng)蛺蝶(Melitaeacinxia)的飛行能力不強,2h飛行距離的中位數(shù)是488m(最短9m,最長2 730m),而且喜歡在開闊地和林地的邊緣飛行,因此試驗在由一圈樹圍成的開闊草地上進行,成功證明了生活在片段化生境下的種群比生活在連續(xù)生境下的種群飛行擴散能力更強[90-91]。

    4 探地雷達在蟻巢探測中的應用

    探地雷達(ground penetrating radar,GPR)又稱地質雷達或透地雷達,是用無線電波來確定地下介質分布的一種方法,在考古、市政建設、建筑、鐵路、公路、水利、電力、采礦、航空等領域都有廣泛應用。中國科學家創(chuàng)造性地將其應用于堤壩白蟻巢的探測,在堤壩隱患排查中發(fā)揮著作用[92-94]。

    探地雷達由主機、收發(fā)天線、圖像處理顯示系統(tǒng)幾個部分組成。實際探測過程中天線沿地面移動,脈沖信號不斷地被發(fā)射和接收,顯示系統(tǒng)經(jīng)A/D轉換后得到的數(shù)據(jù)信號按一定方式進行編碼排列及處理,以二維形式給出連續(xù)的地下縱向剖面圖像[92]。電磁波在地下介質中傳播時遇到存在電性差異的界面時發(fā)生反射,根據(jù)接收到電磁波的波形、振幅強度和時間的變化特征推斷地下介質的空間位置、結構、形態(tài)和埋藏深度。根據(jù)探地雷達的影像特征可以區(qū)分出蟻巢、菌圃群、蟻道、樹根、紅火蟻巢、金屬物等不同的目標[92,95]。

    目前用于探測堤壩蟻巢的探地雷達都是由美國地球物理勘探設備公司(http:∥www.geophysical.com)生產(chǎn)的SIR(subsurface interface radar),該雷達可配備16~2 600MHz范圍內不同波段的天線。探測深度與土壤結構和含水量等因素有關,同樣條件下使用的電磁波頻率越高探測深度越淺[92,94]。使用SIR-3000GPR 400MHz的天線可探測到地下半米左右深度的直徑在10cm左右的蟻巢[94-95]。

    5 我國昆蟲雷達學的機遇與挑戰(zhàn)

    幾乎所有的雷達昆蟲學研究都是以控制重要經(jīng)濟害蟲為主要目的的研究項目的副產(chǎn)品[96]。我國的雷達昆蟲學研究也不例外,新中國成立后東北黏蟲的變幻莫測促使我國昆蟲學家建立了第一部掃描昆蟲雷達,棉鈴蟲的大暴發(fā)促生了第二部昆蟲雷達,進入新世紀草地螟和稻飛虱大暴發(fā),我國又投資建造了多部垂直昆蟲雷達、毫米波昆蟲雷達和多普勒調頻連續(xù)波垂直雷達[4,6,53,54,97-102]。經(jīng)過十幾年的迅速發(fā)展,我國已形成了以垂直昆蟲雷達為主要機型的昆蟲雷達網(wǎng)[6],同時雷達昆蟲學研究團隊也不斷壯大,觀測資料不斷積累。在今后的發(fā)展中應充分發(fā)揮多種雷達共同探測和聯(lián)網(wǎng)探測的優(yōu)勢,努力完善新型雷達軟件系統(tǒng),深入分析已有的觀測資料,不斷提高我國雷達昆蟲學研究和應用的水平。

    5.1 新型昆蟲雷達(網(wǎng))的建設

    出于不同的需求我國近年來研制了垂直昆蟲雷達(包括簡單的垂直雷達、旋轉極化的垂直雷達和收發(fā)分置的多普勒垂直雷達)。但這些雷達的信號采集與分析卻一直存在問題,如只能測到昆蟲的高度與密度[97-101],要得到昆蟲的更多信息必須進一步完善數(shù)字化采集與分析系統(tǒng)。最近發(fā)展起的開源radR平臺[55]有助于改進我國掃描雷達的信號采集與分析系統(tǒng)。我國首創(chuàng)了探地雷達探測蟻巢的方法,如果能進一步研究用來調查地下害蟲(或在地下的蟲態(tài))的密度與變化規(guī)律,將對地下害蟲以及害蟲地下越冬規(guī)律的研究起到不可估量的推動作用。

    我國目前雖然擁有世界上最龐大的昆蟲雷達網(wǎng),但由于沒有支撐其運行的軟件,雷達網(wǎng)尚未得到充分利用。即便是有軟件支持運行,由于各部雷達相距太遠[6],也無法形成有效的覆蓋。天氣雷達網(wǎng)經(jīng)常探測到的昆蟲,對氣象學研究是污染或干擾,但對昆蟲學家卻是非常有用的[103]。雖然現(xiàn)行的氣象雷達不能識別昆蟲的種類,但可以監(jiān)測到昆蟲是無定向的近距離擴散還是進行有定向的遠距離遷飛,單部雷達可監(jiān)測10~30km甚至50km遠和1~2km高[13-14],其覆蓋范圍也遠比昆蟲雷達網(wǎng)廣。

    早在1988年芬蘭人就開始用天氣雷達、誘蟲網(wǎng)和軌跡分析模型研究蚜蟲和小菜蛾向芬蘭的遷飛,并利用手機短信和計算機網(wǎng)絡向試驗用戶發(fā)送害蟲入侵警報[104,105]。事實上,美國的天氣預警雷達網(wǎng)也越來越多地用于昆蟲密度和昆蟲遷飛的常規(guī)觀測[24]。因此,可以設想將我國的天氣雷達網(wǎng)、昆蟲雷達網(wǎng)和誘蟲燈網(wǎng)相互結合形成能預測害蟲發(fā)生和遷飛的實用網(wǎng)絡[15,106],但如何協(xié)調各部門聯(lián)動是一個新的挑戰(zhàn)。

    5.2 毫米波與厘米波昆蟲雷達的協(xié)同觀測

    2007年我國第一部波長為8.0mm的掃描昆蟲雷達建成以來,不但觀測到了稻飛虱的空中密度、飛行高度、空中密度隨時間的變化等預期現(xiàn)象,而且還觀測到了啞鈴形回波[97-99]。但關于毫米波雷達觀測到的啞鈴形回波是不是由稻飛虱的共同定向產(chǎn)生的存在很大爭議。

    Mie散射能計算出任何物體在任何雷達波長下的回波強度,當粒子小于波長的1/10時近似于Rayleigh散射,即回波強度與球體直徑的6次方成正比[25]。對于傳統(tǒng)的3.2cm昆蟲雷達,一般昆蟲的頭尾方向的截面都落入Rayleigh散射區(qū),雷達回波會因反射截面的減小而驟減,因而當昆蟲具有定向時會在雷達屏幕上產(chǎn)生啞鈴形回波[7,56]。毫米波雷達的波長為8mm,一般的遷飛蛾類(如黏蟲、稻縱卷葉螟)不論身體的側面還是正面對準雷達,都不會落入Rayleigh散射區(qū),回波強度的變化不會差別很大,因此這類目標即使是有共同定向,在毫米波雷達上也不應表現(xiàn)出啞鈴形。稻飛虱由于體長與毫米波雷達的波長相似,其頭尾方向的截面落入Rayleigh散射區(qū),理論上可產(chǎn)生啞鈴形回波。有趣的是英國用多普勒雷達也觀測到微小昆蟲的共同定向,而且定向在6月份更偏北一些,在8月份更偏南一些,這與遷飛的季節(jié)性偏好方向一致[13]。因此幾乎可以斷定我國毫米波雷達觀測到的啞鈴形回波是稻飛虱產(chǎn)生的,即稻飛虱具有共同定向特性。

    根據(jù)不同波長雷達探測到的信號差來確定回波的種類在氣象雷達鑒定晴空回波的實踐中曾多次使用[9,25]。根據(jù)雷達原理厘米波昆蟲雷達“看”不到稻飛虱,而毫米波雷達卻可以“看”到稻飛虱,如果同時使用厘米波和毫米波雷達進行觀測,可以很容易地鑒定出啞鈴形分布是稻飛虱產(chǎn)生的還是遷飛的蛾類產(chǎn)生的。當然空中捕捉能得到更直接的證據(jù)。

    5.3 基站式和便攜式諧波雷達的聯(lián)合使用

    基站式掃描諧波雷達非常適合在開闊地帶研究昆蟲的近距離飛行行為,尤其是覓食行為,在蜜蜂神經(jīng)行為學、授粉昆蟲生態(tài)學、氣味介導的趨風性、最優(yōu)搜尋策略和近距離擴散研究中發(fā)揮了重要作用[8]。但當昆蟲飛入茂密的作物冠層或樹林里,或者飛出監(jiān)測范圍,就很難再發(fā)現(xiàn)和追蹤目標。由于諧波雷達不能區(qū)分不同的個體,研究時只能同時放飛1頭或幾頭昆蟲。盡管蜜蜂有“家”的概念,標記的蜜蜂也時有丟失,或者因為飛出監(jiān)測范圍長時間不歸,觀測者只能等待。這時如果繼續(xù)釋放新的標記昆蟲,前面標記過的昆蟲的意外歸來會干擾試驗,這種等待有時可持續(xù)半天或更長。便攜式諧波雷達由于體積小重量輕、可以移動,更適合于尋找進入作物冠層活動的昆蟲。由于兩種雷達可以使用同樣的標簽,如果在農(nóng)業(yè)昆蟲研究中能將兩種類型的雷達同時應用,即根據(jù)基站式雷達指示的方位使用便攜式雷達到較遠或茂密的作物叢中定位昆蟲,也許會得到意想不到的結果。

    5.4 現(xiàn)有觀測數(shù)據(jù)的分析與新的挑戰(zhàn)

    昆蟲定向與成層行為機理的研究在本世紀逐漸成為雷達昆蟲學研究的熱點,也取得了一些進展,但相互矛盾和不解之謎甚多。我國經(jīng)過十年的持續(xù)觀測積累了大量的夜間昆蟲遷飛的資料,而且都有同時測到的高空風溫場資料,這是一筆寶貴的資源。我們已初步建立了用風溫場模擬昆蟲成層的方法,現(xiàn)正在建立計算機分析程序,以期對歷史資料進行更深入的分析,明確我國遷飛昆蟲的成層機理。但昆蟲學工作者缺乏編程經(jīng)驗,而程序員不懂昆蟲,交叉學科人才不足給這項看似簡單的工作平添幾分困難。

    越來越多的觀測表明昆蟲具有羅盤定向能力,但它們是使用什么羅盤還不清楚。使用雷達雖然可以準確觀測昆蟲的定向,但難以控制幾百米高空的磁場、氣壓和風溫等環(huán)境條件。紅外三維軌跡捕捉及分析系統(tǒng)的嘗試也未取得好的結果,飛行模擬器可以在模擬環(huán)境下測出昆蟲的定向,但并不是所有的昆蟲都適合吊飛[107-108]。最近小型昆蟲的定向和成層也逐漸被認識[13,97],為昆蟲的成層和定向機理研究提出了新的課題和挑戰(zhàn)。

    [1]程登發(fā),封洪強,吳孔明.掃描昆蟲雷達與昆蟲遷飛監(jiān)測[M].北京:科學出版社,2005.

    [2]百度百科.http:∥baike.baidu.com,2011.

    [3]維基百科.雷達[EB/OL].http:∥zh.wikipedia.org/zh-cn/%E9%9B%B7%E9%81%94,2011.

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

    [5]翟保平.昆蟲雷達讓我們看到了什么?[J].昆蟲知識,2005,42(2):217-226.

    [6]封洪強.雷達昆蟲學40年研究的回顧與展望[J].河南農(nóng)業(yè)科學,2009(9):121-126.

    [7]Riley J R.Riding on the wind:a radar perspective of insect flight[M].Inaugural Lecture Series,the University of Greenwich,London,1999:31.

    [8]Chapman J W,Drake V A,Reynolds D R.Recent insights from radar studies of insect flight[J].Annu Rev Entomol,2011,56:337-356.

    [9]Wilson J W,Weckwerth T M,Vivekanandan J,et al.Boundary layer clear-air radar echoes:origin of echoes and accuracy of derived winds[J].Journal of Atmospheric and Oceanic Technology,1994,11:1184-1206.

    [10]Achtemeier G L.The use of insects as tracers for“clear-air”boundary-layer studies by Doppler radar[J].Journal of Atmospheric and Oceanic Technology,1991,8:746-765.

    [11]Contreras R F,F(xiàn)rasier S J.High-resolution observations of insects in the atmospheric boundary layer[J].Journal of Atmospheric and Oceanic Technology,2008,25:2176-2187.

    [12]Luke E P,Kollias P,Johnson K L,et al.A technique for the automatic detection of insect clutter in cloud radar[J].Journal of Atmospheric and Oceanic Technology,2008,25:1498-1513.

    [13]Rennie S J,Illingworth A J,Dance S L,et al.The accuracy of Doppler radar wind retrievals using insects as targets[J].Meteorological Applications,2010,17(4):419-432.

    [14]Rennie S J,Illingworth A J,Dance S L,et al.Utilization of Doppler radar wind measurements from insect returns[C]∥ERAD 2008-The Fifth European Conference on Radar in Meteorology and Hydrology,Helsinki,F(xiàn)inland,30June-4July 2008.Finnish Meteorological Institute:Helsinki,2008.

    [15]Wood C R,O’Connor E J,Hurley R A,et al.Cloud-radar observations of insects in the UK convective boundary layer[J].Meteorological Applications,2009,16(4):491-500.

    [16]Gauthreaux Jr S A,Livingston J W,Belser C G.Detection and discrimination of fauna in the aerosphere using Doppler weather surveillance radar[J].Integr Comp Biol,2008,48(1):12-23.

    [17]Leskinen M.Observed polarimetric signals of insects[C]∥ERAD 2008-The Fifth European Conference on Radar in Meteorology and Hydrology,Helsinki,F(xiàn)inland,30June-4July 2008.Finnish Meteorological Institute:Helsinki,2008.

    [18]Bachman S,Zrnic'D.Three-dimensional attributes of clear-air scatterers observed with the polarimetric weather radar[J].IEEE Geosciences and Remote Sensing Letters,2008,25:231-235.

    [19]Mueller E A,Larkin R P.Insects observed using dual-polarization radar[J].Journal of Atmospheric and Oceanic Technology,1985,2:49-54.

    [20]Zrnic D S,Ryzhkov A V.Observations of insects and birds with a polarimetric radar[J].IEEE Transactions on Geoscience and Remote Sensing,1998,36(2):661-668.

    [21]Rennie S,Illingworth A,Dance S.On differentiating ground clutter and insect echoes from Doppler weather radars using archived data[J].Atmospheric Measurement Techniques Discussions,2010,3:1843-1860.

    [22]Bachman S,Zrnic'D.Spectral density of polarimetric variables separating biological scatterers in the VAD display[J].Journal of Atmospheric and Oceanic Technology,2007,24:1186-1198.

    [23]Geerts B,Miao Q.The use of millimeter Doppler radar echoes to estimate vertical air velocities in the fair-weather convective boundary layer[J].Journal of Atmospheric and Oceanic Technology,2005,22:225-246.

    [24]Geerts B,Miao Q.Airborne radar observations of the flight behavior of small insects in the atmospheric convective boundary layer[J].Environ Entomol,2005,34:361-377.

    [25]Martin W J,Shapiro A.Discrimination of bird and insect radar echoes in clear air using high-resolution radars[J].Journal of Atmospheric and Oceanic Technology,2007,24:1215-1230.

    [26]Mueller E A.Differential reflectivity of birds and insects[C]∥Preprints,21st Conference on Radar Meteorology,Edmonton,AB,Canada:American Meteorological Society,1983:465-466.

    [27]Schaefer G W.Radar observations of insect flight[C]∥Symposium of the Royal Entomology Society,London,1976,7:157-197.

    [28]Drake V A.Automatically operating radars for monitoring insect pest migrations[J].Entomologia Sinica,2002,9(4):27-39.

    [29]Chapman J W,Reynolds D R,Smith A D.Vertical-looking radar:A new tool for monitoring high-altitude insect migration[J].BioScience,2003,53(5):503-511.

    [30]Chapman J W,Smith A D,Woiwod I P,et al.Development of vertical-looking radar technology for monitoring insect migration[J].Computers and Electronics in Agriculture,2002,35(2-3):95-110.

    [31]翟保平.昆蟲雷達:從研究型到實用型[J].遙感學報,2001,5(3):231-240.

    [32]Smith A D,Reynolds D R,Riley J R.The use of vertical-looking radar to continuously monitor the insect fauna flying at altitude over southern England[J].Bull Entomol Res,2000,90(3):265-277.

    [33]Smith A D,Riley J R,Gregory R D.A method for routine monitoring of the aerial migration of insects by using a verticallooking radar[J].Philosophical Transactions of the Royal Society of London.Series B:Biological Sciences,1993,340(1294):393-404.

    [34]Drake V A,Harman I T,Wang H K.Insect monitoring radar:stationary-beam operating mode[J].Computers and Electronics in Agriculture,2002,35(2-3):111-137.

    [35]Harman I T,Drake V A.Insect monitoring radar:analytical time-domain algorithm for retrieving trajectory and target parameters[J].Computers and Electronics in Agriculture,2004,43(1):23-41.

    [36]Wang H K,Drake V A.Insect monitoring radar:retrieval of wingbeat information from conical-scan observation data[J].Computers and Electronics in Agriculture,2004,43(3):209-222.

    [37]Dean T J,Drake V A.Monitoring insect migration with radar:the ventral-aspect polarization pattern and its potential for target identification[J].International Journal of Remote Sensing,2005,26(18):3957-3974.

    [38]Reynolds D R,Chapman J W,Brooks S J,et al.Seasonal variation in the migration strategies of the green lacewingChrysoperlacarneaspecies complex[J].Ecological Entomology,2006,31(4):378-388.

    [39]Chapman J W,Reynolds D R,Smith A D,et al.High-altitude migration of the diamondback mothPlutellaxylostellato the U.K.:a study using radar,aerial netting,and ground trapping[J].Ecological Entomology,2002,27(6):641-650.

    [40]Drake V A,Wang H K,Harman I T.Insect monitoring radar:remote and network operation[J].Computers and Electronics in Agriculture,2002,35(2-3):77-94.

    [41]Cheng D F,Wu K M,Tian Z,et al.Acquisition and analysis of migration data from the digitised display of a scanning entomological radar[J].Computers and Electronics in Agriculture,2002,35(2-3):63-75.

    [42]程登發(fā),吳孔明,田喆,等.掃描昆蟲雷達實時數(shù)據(jù)采集、分析系統(tǒng)[J].植物保護,2004,30(2):41-46.

    [43]Drake V A.Quantitative observation and analysis procedures for a manually operated entomological radar[R].CSIRO Australia Division of Entomology Technical Paper No.19,1981:1-41.

    [44]Feng H Q,Wu K M,Cheng D F,et al.Radar observations of the autumn migration of the beet armywormSpodopteraexigua(Lepidoptera:Noctuidae)and other moths in northern China[J].Bull Entomol Res,2003,93(2):115-124.

    [45]Feng H Q,Wu X F,Wu B,et al.Seasonal migration ofHelicoverpaarmigera(Lepidoptera:Noctuidae)over the Bohai Sea[J].J Econ Entomol,2009,102(1):95-104.

    [46]Feng H Q,Wu K M,Cheng D F,et al.Northward migration ofHelicoverpaarmigera(Lepidoptera:Noctuidae)and other moths in early summer observed with radar in northern China[J].J Econ Entomol,2004,97(6):1874-1883.

    [47]Feng H Q,Wu K M,Ni Y X,et al.Return migration ofHelicov-erpaarmigera(Lepidoptera:Noctuidae)during autumn in northern China[J].Bull Entomol Res,2005,95(4):361-370.

    [48]Feng H Q,Zhao X C,Wu X F,et al.Autumn migration ofMythimnaseparata(Lepidoptera:Noctuidae)over the Bohai Sea in northern China[J].Environ Entomol,2008,37(3):774-781.

    [49]Feng H,Wu K M,Cheng D F,et al.Spring migration and summer dispersal ofLoxostegesticticalis(Lepidoptera:Pyralidae)and other insects observed with radar in Northern China[J].Environ Entomol,2004,33(5):1253-1265.

    [50]Feng H Q,Wu K M,Ni Y X,et al.High-altitude windborne transport ofHelicoverpaarmigera(Lepidoptera:Noctuidae)and other moths in mid-summer in Northern China[J].Journal of Insect Behavior,2005,18(3):335-349.

    [51]Feng H Q,Wu K M,Ni Y X,et al.Nocturnal migration of dragonflies over the Bohai Sea in northern China[J].Ecological Entomology,2006,31(5):511-520.

    [52]Feng H Q,Zhang Y H,Wu K M,et al.Nocturnal windborne migration of ground beetles,particularlyPseudoophonusgriseus(Coleoptera:Carabidae),in China[J].Agricultural and Forest Entomology,2007,9(2):103-113.

    [53]吳孔明,翟保平,封洪強,等.華北北部地區(qū)二代棉鈴蟲成蟲遷飛行為的雷達觀測[J].植物保護學報,2006,33(2):163-167.

    [54]吳孔明,程登發(fā),徐廣,等.華北地區(qū)昆蟲秋季遷飛的雷達觀測[J].生態(tài)學報,2001,21(11):1833-1838.

    [55]Taylor P D,Brzustowski J M,Matkovich C,et al.radR:an open-source platform for acquiring and analysing data on biological targets observed by surveillance radar[J].BMC Ecol,2010,10:22.

    [56]Riley J R.Collective orientation in night-flying insects[J].Nature,1975,253(5487):113-114.

    [57]Riley J R,Reynolds D R.Orientation at night by high-flying insects[M]∥Danthanarayana W,ed.Insect Flight:Dispersal and Migration.Berlin:Springer-Verlag,1986:71-87.

    [58]Reynolds D R,Riley J R.The flight behaviour and migration of insect pests:radar studies in developing countries[R].Chatham (United Kingdom):Natural Resources Institute.NRI Bulletin no.71,1997.

    [59]Drake V A.Collective orientation by nocturnally migrating Australian plague locusts,Chortoicetesterminifera(Walker)(Orthoptera:Acrididae):a radar study[J].Bull Entomol Res,1983,73(04):679-692.

    [60]Reynolds A M,Reynolds D R,Smith A D,et al.A single wind-mediated mechanism explains high-altitude‘non-goal oriented’headings and layering of nocturnally migrating insects[J].Proceedings of the Royal Society B:Biological Sciences,2010,277(1682):765-772.

    [61]Chapman J W,Nesbit R L,Burgin L E,et al.Flight orientation behaviors promote optimal migration trajectories in highflying insects[J].Science,2010,327(5966):682-685.

    [62]Chapman J W,Reynolds D R,Hill J K,et al.A seasonal switch in compass orientation in a high-flying migrant moth[J].Curr Biol,2008,18(19):R908-909.

    [63]Chapman J W,Nesbit R L,Burgin L E,et al.Flight orientation behaviors promote optimal migration trajectories in highflying insects[J].Science,2010,327(5966):682-685.

    [64]Chapman J W,Reynolds D R,Mouritsen H,et al.Wind selection and drift compensation optimize migratory pathways in a high-flying moth[J].Curr Biol,2008,18(7):514-518.

    [65]Reynolds D R,Chapman J W,Edwards A S,et al.Radar studies of the vertical distribution of insects migrating over southern Britain:the influence of temperature inversions on nocturnal layer concentrations[J].Bull Entomol Res,2005,95(3):259-274.

    [66]Wood C R,Chapman J W,Reynolds D R,et al.The influence of the atmospheric boundary layer on nocturnal layers of noctuids and other moths migrating over southern Britain[J].Int J Biometeorol,2006,50(4):193-204.

    [67]Wood C R,Clarke S J,Barlow J F,et al.Layers of nocturnal insect migrants at high-altitude:the influence of atmospheric conditions on their formation[J].Agricultural and Forest Entomology,2010,12:113-121.

    [68]Colpitts B G,Boiteau G.Harmonic radar transceiver design:miniature tags for insect tracking[J].IEEE Transactions on Antennas and Propagation,2004,52(11):2825-2832.

    [69]桂連友,黃秀琴,李傳仁.昆蟲諧波雷達的發(fā)展和利用[J].應用昆蟲學報,2011,48(3):732-738.

    [70]Williams D W,Li G,Gao R.Tracking movements of individualAnoplophoraglabripennis(Coleoptera:Cerambycidae)adults:application of harmonic radar[J].Environmental Entomology,2004,33(3):644-649.

    [71]Riley J R,Smith A D.Design considerations for an harmonic radar to investigate the flight of insects at low altitude[J].Computers and Electronics in Agriculture,2002,35(2-3):151-169.

    [72]Boiteau G,Vincent C,Meloche F,et al.Harmonic radar:assessing the impact of tag weight on walking activity of Colorado potato beetle,Plumcurculio,and western corn rootworm[J].J Econ Entomol,2010,103(1):63-69.

    [73]Riley J R,Smith A D,Reynolds D R,et al.Tracking bees with harmonic radar[J].Nature,1996,379:29-30

    [74]O’Neal M E,Landis D A,Rothwell E,et al.Tracking insects with harmonic radar:a case study[J].American Entomologist,2004,50(4):212-218.

    [75]Boiteau G,Colpitts B.The potential of portable harmonic radar technology for the tracking of beneficial insects[J].International Journal of Pest Management,2004,50(3):233-242.

    [76]Mascanzoni D,Wallin H.The harmonic radar:a new method of tracing insects in the field[J].Ecological Entomology,1986,11(4):387-390.

    [77]Wallin H,Ekbom B S.Movements of carabid beetles(Coleoptera:Carabidae)inhabiting cereal fields:a field tracing study[J].Oecologia,1988,77(1):39-43.

    [78]Wallin H,Ekbom B.Influence of hunger level and prey densities on movement patterns in three species of pterostichus beetles(Coleoptera:Carabidae)[J].Environmental Entomology,1994,23(5):1171-1181.

    [79]Charrier S,Petit S,Burel F.Movements ofAbaxparallelepipedus(Coleoptera,Carabidae)in woody habitats of a hedgerow network landscape:a radio-tracing study[J].Agriculture,Ecosystems & Environment,1997,61(2-3):133-144.

    [80]Bliteau G,Colpitts B.Electronic tags for the tracking of insects in flight:effect of weight on flight performance of adult Colorado potato beetles[J].Entomologia Experimentalis et Applicata,2001,100:187-193.

    [81]Boiteau G,Vincent C,Meloche F,et al.Harmonic radar:efficacy at detecting and recovering insects on agricultural host plants[J].Pest Manag Sci,2011,67(2):213-9.

    [82]Boiteau G,Vincent C,Meloche F,et al.Evaluation of tag entanglement as a factor in harmonic radar studies of insect dispersal[J].Environ Entomol,2011,40(1):94-102.

    [83]Gui L Y,Boiteau G,Colpitts B G,et al.Random movement pattern of fed and unfed adult Colorado potato beetles in bareground habitat[J].Agricultural and Forest Entomology,2011.In press.

    [84]Gui L Y,Huang X Q,Li C R,et al.Validation of harmonic radar tags to study the movement of Chinese citrus fly[J].The Canadian Entomologist,2011:In press.

    [85]Carreck N L,Osborn J L,Capaldi E A,et al.Tracking bees with radar[J].Bee World,1999,80(3):124-131.

    [86]Menzel R,F(xiàn)uchs J,Nadler L,et al.Dominance of the odometer over serial landmark learning in honeybee navigation[J].Naturwissenschaften,2010,97(8):763-767.

    [87]Svensson G P,Valeur P G,Reynolds D R,et al.Mating disruption inAgrotissegetummonitored by harmonic radar[J].Entomologia Experimentalis et Applicata,2001,101(2):111-121.

    [88]Reynolds A M,Reynolds D R,Smith A D,et al.Appetitive flight patterns of maleAgrotissegetummoths over landscape scales[J].J Theor Biol,2007,245(1):141-149.

    [89]Cant E T,Smith A D,Reynolds D R,et al.Tracking butterfly flight paths across the landscape with harmonic radar[J].Proc Biol Sci,2005,272(1565):785-790.

    [90]Ovaskainen O,Smith A D,Osborne J L,et al.Tracking butterfly movements with harmonic radar reveals an effect of population age on movement distance[J].Proc Natl Acad Sci U S A,2008,105(49):19090-19095.

    [91]Wang R,Ovaskainen O,Cao Y D,et al.Dispersal in the Glanville fritillary butterfly in fragmented versus continuous landscapes:comparison between three methods[J].Ecological Entomology,2011,36:251-260.

    [92]徐興新,李棟.探地雷達探測堤壩白蟻巢研究[J].昆蟲學報,1996,39(1):46-52.

    [93]Henderson G.The termite menace in New Orleans:did they cause the floodwalls to tumble?[J].American Entomologist,2008,54(3):156-162.

    [94]陳家銀,李富長,司馬軍.白龜山水庫白蟻隱患的探地雷達檢測試驗[J].河南水利與南水北調,2007,2007(1):54-55.

    [95]Yang X H,Henderson G,Mao L X,et al.Application of ground penetrating radar in detecting the hazards and risks of termites and ants in soil levees[J].Environ Entomol,2009,38(4):1241-1249.

    [96]Drake V A,Reynolds D R.Radar entomology:observing insect flight and migration[M].CABI,2011:304.

    [97]齊會會.稻飛虱的雷達監(jiān)測與遷飛軌跡分析[D].北京:中國農(nóng)業(yè)科學院,2011.

    [98]齊會會,張云慧,程登發(fā),等.褐飛虱2009年秋季回遷的雷達監(jiān)測及軌跡分析[J].昆蟲學報,2010,53(11):1256-1264.

    [99]楊秀麗,陳林,程登發(fā),等.毫米波掃描昆蟲雷達空中昆蟲監(jiān)測的初步應用[J].植物保護,2008,34(2):31-36.

    [100]高月波,陳曉,陳鐘榮,等.稻縱卷葉螟(Cnaphalocrocismedinalis)遷飛的多普勒昆蟲雷達觀測及動態(tài)[J].生態(tài)學報,2008,28(11):5238-5247.

    [101]張云慧,陳林,程登發(fā),等.草地螟2007年越冬代成蟲遷飛行為研究與蟲源分析[J].昆蟲學報,2008,51(7):720-727.

    [102]張云慧,喬紅波,程登發(fā),等.垂直監(jiān)測昆蟲雷達空中昆蟲監(jiān)測的初步應用[J].植物保護,2007,33(3):23-26.

    [103]Zhang P F,Ryzhkov A,Zrnic D.Polarimetric prototype of the WSR-88Dradar observation of birds and insects[C]∥Preprints CD,12th Conference on Aviation,Range,and Aerospace Meteorology,Atlanta,GA,28January-3February 2006.American Meteorological Society:Boston,2006.

    [104]Leskinen M,Markkula I,Koistinen J,et al.Pest insect immigration warning by an atmospheric dispersion model,weather radars and traps[J].Journal of Applied Entomology,2010,135(1-2):55-67.

    [105]Nieminen M,Leskinen M,Helenius J.Doppler radar detection of exceptional mass-migration of aphids into Finland[J].Int J Biometeorol,2000,44(4):172-181.

    [106]Markkula I,Ojanen H,Tiilikkala K,et al.Insect migration case studies by polarimetric radar[C].In ERAD 2008-The Fifth European Conference on Radar in Meteorology and Hydrology,Helsinki,F(xiàn)inland,30June-4July 2008.Finnish Meteorological Institute:Helsinki,2008.

    [107]吳博,封洪強,趙奎軍,等.視頻軌跡分析技術在昆蟲行為學研究中的應用[J].植物保護,2009,35(2):1-6.

    [108]施翔宇,封洪強,劉中芳,等.實驗條件下黏蟲、棉鈴蟲和小地老虎的定向行為比較[J].植物保護,2010,36(4):60-63.

    Applications of radar in entomological research

    Feng Hongqiang
    (HenanKeyLaboratoryofCropPestControl,InstituteofPlantProtection,HenanAcademyofAgriculturalSciences,Zhengzhou450002,China)

    Insect radar has been used as an excellent tool in entomological research,such as flight behavior of insects during high-altitude and long-distance migration or just above ground,because of its fascinating ability in detecting flight behavior of insects in natural free flying from a long distance with little disturbance to them.Weather radar has increasingly been used to monitor birds and insects,and GPR (Ground Penetrating Radar)has also been used to detect nests of termites and ants in soil levees.This review evaluated the capability of different types of radars and their performance in entomological research,to tell what the radar can do and the challenges in using radar from the view of an entomologist.

    radar; insect; migration; target identification; application

    S 431.9

    A

    10.3969/j.issn.0529-1542.2011.05.001

    2011-09-05

    國家自然科學基金項目(30971873);河南省杰出青年基金項目(094100510024);河南省重大公益項目(081100911300)

    * 通信作者E-mail:feng_hq@163.com

    猜你喜歡
    昆蟲學定向昆蟲
    RFID昆蟲閱讀放大鏡
    玩具世界(2022年3期)2022-09-20 01:48:20
    森林昆蟲學經(jīng)典學術著作再版
    借昆蟲上課
    甘肅教育(2020年2期)2020-09-11 08:01:48
    法布爾
    我最喜歡的昆蟲——知了
    偏序集上的相對定向集及其應用
    昆蟲的冬天
    定向越野
    昆蟲學
    昆蟲學
    .国产精品久久| 卡戴珊不雅视频在线播放| 一区二区三区四区激情视频| 三级国产精品片| 欧美97在线视频| 国产精品一区二区三区四区免费观看| av在线观看视频网站免费| 2018国产大陆天天弄谢| 男人和女人高潮做爰伦理| 国产精品久久久久久久久免| 欧美人与善性xxx| 边亲边吃奶的免费视频| 久久人人爽人人爽人人片va| 国模一区二区三区四区视频| 中文在线观看免费www的网站| 日本熟妇午夜| 亚洲熟女精品中文字幕| 国产av国产精品国产| 青春草国产在线视频| 日韩一区二区视频免费看| 99视频精品全部免费 在线| 看非洲黑人一级黄片| 丰满乱子伦码专区| 久久久久性生活片| 色5月婷婷丁香| 国产成人精品久久久久久| 亚洲国产av新网站| 亚洲精品国产av成人精品| 久久久久久久久久久丰满| 久久鲁丝午夜福利片| 久久99蜜桃精品久久| 丰满乱子伦码专区| 国产精品国产三级国产专区5o| av专区在线播放| 色网站视频免费| 亚洲精品成人久久久久久| 国产精品女同一区二区软件| 干丝袜人妻中文字幕| kizo精华| 中文欧美无线码| 男女国产视频网站| 建设人人有责人人尽责人人享有的 | 又爽又黄无遮挡网站| 小蜜桃在线观看免费完整版高清| av在线蜜桃| 亚洲精品色激情综合| 欧美成人一区二区免费高清观看| 日本一本二区三区精品| av网站免费在线观看视频 | 成人鲁丝片一二三区免费| 亚洲欧美成人精品一区二区| 亚洲精品久久久久久婷婷小说| 麻豆国产97在线/欧美| 亚洲av免费高清在线观看| 国产一区有黄有色的免费视频 | 日韩成人av中文字幕在线观看| 亚洲人与动物交配视频| 人妻一区二区av| 国产av在哪里看| 国产精品爽爽va在线观看网站| 亚洲成人久久爱视频| 啦啦啦韩国在线观看视频| 寂寞人妻少妇视频99o| 丰满人妻一区二区三区视频av| 日韩欧美精品免费久久| 一级爰片在线观看| 亚洲精品aⅴ在线观看| 国产成年人精品一区二区| 久久久久九九精品影院| 国产探花在线观看一区二区| 亚洲精品aⅴ在线观看| 国产av不卡久久| 在线 av 中文字幕| 人人妻人人看人人澡| 国产高潮美女av| 国产精品一区二区三区四区久久| 国产美女午夜福利| 精品久久久久久久人妻蜜臀av| 久久久亚洲精品成人影院| 色哟哟·www| 国产精品久久久久久精品电影小说 | 好男人视频免费观看在线| 日本三级黄在线观看| 欧美最新免费一区二区三区| 黑人高潮一二区| 99久国产av精品国产电影| 99热这里只有是精品50| 亚洲激情五月婷婷啪啪| 18禁在线播放成人免费| 蜜桃久久精品国产亚洲av| .国产精品久久| 一级毛片 在线播放| 久久久久久久午夜电影| 一区二区三区乱码不卡18| 国产在线一区二区三区精| av天堂中文字幕网| 久久精品国产自在天天线| 国产精品人妻久久久影院| 男插女下体视频免费在线播放| 水蜜桃什么品种好| 少妇的逼水好多| 边亲边吃奶的免费视频| 偷拍熟女少妇极品色| 青春草国产在线视频| 高清视频免费观看一区二区 | 久久精品国产鲁丝片午夜精品| 人妻一区二区av| 国产精品一区二区在线观看99 | 亚洲成人久久爱视频| 亚洲av福利一区| 亚洲精品视频女| 日韩av免费高清视频| 一夜夜www| 91精品国产九色| 91狼人影院| 亚洲欧美精品专区久久| 中国国产av一级| 最近的中文字幕免费完整| 成人欧美大片| 美女脱内裤让男人舔精品视频| 免费观看av网站的网址| 99久久九九国产精品国产免费| 人妻夜夜爽99麻豆av| 欧美bdsm另类| 亚洲一级一片aⅴ在线观看| 好男人视频免费观看在线| 久久99热6这里只有精品| 国产午夜精品一二区理论片| 免费观看在线日韩| 精品一区二区免费观看| 午夜福利视频精品| 18禁动态无遮挡网站| 午夜视频国产福利| 一个人看的www免费观看视频| av.在线天堂| 美女主播在线视频| 亚洲国产欧美人成| 97热精品久久久久久| 欧美极品一区二区三区四区| 人妻一区二区av| 老司机影院毛片| 汤姆久久久久久久影院中文字幕 | 亚洲欧美日韩卡通动漫| 麻豆成人午夜福利视频| 国产久久久一区二区三区| 国产伦一二天堂av在线观看| 1000部很黄的大片| 日本色播在线视频| 国产精品一及| 中国国产av一级| 丝袜美腿在线中文| 亚洲av电影在线观看一区二区三区 | 亚洲成人av在线免费| 神马国产精品三级电影在线观看| 我的女老师完整版在线观看| 一夜夜www| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| 午夜福利在线观看免费完整高清在| 男的添女的下面高潮视频| 亚洲一级一片aⅴ在线观看| 春色校园在线视频观看| av在线老鸭窝| 一个人看的www免费观看视频| 汤姆久久久久久久影院中文字幕 | 久99久视频精品免费| 亚洲成人一二三区av| 久久久精品免费免费高清| 九九在线视频观看精品| 国产一区二区在线观看日韩| 国产亚洲最大av| 国产人妻一区二区三区在| 可以在线观看毛片的网站| 三级男女做爰猛烈吃奶摸视频| av又黄又爽大尺度在线免费看| 97精品久久久久久久久久精品| 特级一级黄色大片| 国产成人91sexporn| 久久久久精品性色| 99久久九九国产精品国产免费| 黑人高潮一二区| 久久精品国产亚洲网站| 黑人高潮一二区| 97精品久久久久久久久久精品| 十八禁网站网址无遮挡 | 真实男女啪啪啪动态图| 亚洲欧美日韩无卡精品| 国产精品熟女久久久久浪| 午夜视频国产福利| 亚洲av男天堂| 秋霞在线观看毛片| 精品国产露脸久久av麻豆 | 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99 | 高清在线视频一区二区三区| 亚洲精品日韩在线中文字幕| 嫩草影院入口| 久久亚洲国产成人精品v| 亚洲熟女精品中文字幕| 国产成人aa在线观看| 亚洲欧美日韩东京热| av免费在线看不卡| 嘟嘟电影网在线观看| 国产av码专区亚洲av| 亚洲欧洲日产国产| 国产在视频线在精品| 日韩精品青青久久久久久| 一个人观看的视频www高清免费观看| 精品久久久噜噜| 亚洲自拍偷在线| 国语对白做爰xxxⅹ性视频网站| 人人妻人人看人人澡| 亚洲av成人精品一二三区| 成人av在线播放网站| 亚洲最大成人中文| 欧美精品国产亚洲| 国产久久久一区二区三区| 中文字幕av成人在线电影| 久久6这里有精品| 午夜精品一区二区三区免费看| 亚洲精品,欧美精品| 一区二区三区免费毛片| 内射极品少妇av片p| 亚洲自偷自拍三级| 亚洲精品乱码久久久v下载方式| 国内少妇人妻偷人精品xxx网站| 欧美zozozo另类| 久久精品熟女亚洲av麻豆精品 | 亚洲欧美中文字幕日韩二区| 大陆偷拍与自拍| 亚洲综合色惰| 国产亚洲精品久久久com| av黄色大香蕉| 国产av码专区亚洲av| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区四区激情视频| 色尼玛亚洲综合影院| 少妇高潮的动态图| 久久久久久伊人网av| 精品久久久久久久末码| av在线老鸭窝| 高清午夜精品一区二区三区| 又粗又硬又长又爽又黄的视频| 久久久久九九精品影院| 免费人成在线观看视频色| 亚洲欧美一区二区三区黑人 | 久久99热这里只有精品18| 少妇猛男粗大的猛烈进出视频 | 国产成人精品一,二区| 寂寞人妻少妇视频99o| 免费观看无遮挡的男女| 人妻一区二区av| 久久人人爽人人爽人人片va| 丝袜喷水一区| 精品一区二区免费观看| 狠狠精品人妻久久久久久综合| 亚洲精品乱码久久久久久按摩| 日韩av在线大香蕉| 高清视频免费观看一区二区 | 十八禁网站网址无遮挡 | 黄色配什么色好看| 偷拍熟女少妇极品色| 联通29元200g的流量卡| 日韩大片免费观看网站| 亚洲不卡免费看| 国产三级在线视频| 99热全是精品| 我的老师免费观看完整版| 一个人观看的视频www高清免费观看| 欧美高清成人免费视频www| 在线免费观看的www视频| 91精品国产九色| 亚洲精品乱码久久久久久按摩| 亚洲va在线va天堂va国产| 国产成人freesex在线| 国国产精品蜜臀av免费| 免费观看精品视频网站| 亚洲18禁久久av| 精品久久久久久成人av| 日本色播在线视频| 超碰97精品在线观看| 国产永久视频网站| 在线 av 中文字幕| 欧美xxⅹ黑人| 大香蕉97超碰在线| 在线天堂最新版资源| 日日摸夜夜添夜夜爱| 97超碰精品成人国产| av又黄又爽大尺度在线免费看| 国产久久久一区二区三区| 深夜a级毛片| 乱系列少妇在线播放| 欧美区成人在线视频| 免费观看无遮挡的男女| 国内精品一区二区在线观看| 观看美女的网站| 久久久精品免费免费高清| 日韩大片免费观看网站| 国产在视频线在精品| 日韩伦理黄色片| 国产亚洲一区二区精品| 只有这里有精品99| 久久人人爽人人片av| 大话2 男鬼变身卡| 大香蕉97超碰在线| 国产午夜福利久久久久久| 午夜福利在线观看吧| 久久久久国产网址| 亚洲av一区综合| 亚洲av中文字字幕乱码综合| 九色成人免费人妻av| 天堂影院成人在线观看| 亚洲国产精品专区欧美| 精品久久久久久久久久久久久| 国产不卡一卡二| 女人被狂操c到高潮| 久久热精品热| 亚洲人与动物交配视频| 精品一区二区三区人妻视频| 午夜精品在线福利| 国产一级毛片七仙女欲春2| 午夜老司机福利剧场| 国产成人精品福利久久| 69av精品久久久久久| 欧美+日韩+精品| 国产精品爽爽va在线观看网站| 久久精品熟女亚洲av麻豆精品 | 高清毛片免费看| 亚洲熟女精品中文字幕| 国产单亲对白刺激| 精华霜和精华液先用哪个| 国产视频内射| 嫩草影院精品99| 亚洲av电影不卡..在线观看| 国产不卡一卡二| 久久草成人影院| 中文欧美无线码| 韩国av在线不卡| 亚洲精品亚洲一区二区| 秋霞在线观看毛片| 亚洲乱码一区二区免费版| 最近2019中文字幕mv第一页| 91午夜精品亚洲一区二区三区| 久久久久精品性色| 22中文网久久字幕| 国产麻豆成人av免费视频| 91久久精品电影网| 能在线免费观看的黄片| 亚洲精品亚洲一区二区| 久久精品国产亚洲网站| 91aial.com中文字幕在线观看| 韩国av在线不卡| 国产黄色免费在线视频| 麻豆精品久久久久久蜜桃| freevideosex欧美| 国产69精品久久久久777片| 汤姆久久久久久久影院中文字幕 | 免费观看在线日韩| 中文乱码字字幕精品一区二区三区 | 欧美一区二区亚洲| 日韩国内少妇激情av| 美女xxoo啪啪120秒动态图| 中文字幕av在线有码专区| 国产成人福利小说| 欧美变态另类bdsm刘玥| 听说在线观看完整版免费高清| 欧美高清成人免费视频www| 91久久精品国产一区二区三区| 国产激情偷乱视频一区二区| 男女国产视频网站| 亚洲无线观看免费| 国产黄色小视频在线观看| 国产伦在线观看视频一区| 深爱激情五月婷婷| 成人高潮视频无遮挡免费网站| 国产黄色小视频在线观看| 一个人看视频在线观看www免费| 午夜亚洲福利在线播放| 色综合色国产| 国模一区二区三区四区视频| 亚洲第一区二区三区不卡| 春色校园在线视频观看| 亚洲精品aⅴ在线观看| 简卡轻食公司| 国产黄色免费在线视频| 免费看av在线观看网站| eeuss影院久久| 成人亚洲精品av一区二区| 日韩伦理黄色片| 蜜桃久久精品国产亚洲av| 99久久精品热视频| 国产成人a区在线观看| 中文精品一卡2卡3卡4更新| 免费少妇av软件| 亚洲熟妇中文字幕五十中出| 中文字幕人妻熟人妻熟丝袜美| 日产精品乱码卡一卡2卡三| 午夜爱爱视频在线播放| 亚洲精品日韩在线中文字幕| 好男人在线观看高清免费视频| 黄色配什么色好看| 欧美日韩国产mv在线观看视频 | 精品酒店卫生间| 一个人观看的视频www高清免费观看| 黄色日韩在线| 毛片女人毛片| 身体一侧抽搐| 三级毛片av免费| 日产精品乱码卡一卡2卡三| 色视频www国产| 97人妻精品一区二区三区麻豆| 特大巨黑吊av在线直播| 日本wwww免费看| 最近最新中文字幕大全电影3| 亚洲自偷自拍三级| 午夜免费男女啪啪视频观看| 校园人妻丝袜中文字幕| 欧美最新免费一区二区三区| 国产亚洲精品久久久com| 午夜福利网站1000一区二区三区| 久久人人爽人人片av| 国产黄频视频在线观看| 亚洲国产成人一精品久久久| 久久久久久久久大av| 亚洲人成网站高清观看| av一本久久久久| 日韩成人伦理影院| 丝袜美腿在线中文| 亚洲国产高清在线一区二区三| 老司机影院成人| 91午夜精品亚洲一区二区三区| 国产视频内射| 3wmmmm亚洲av在线观看| 国产综合懂色| 成人特级av手机在线观看| 亚洲国产日韩欧美精品在线观看| 啦啦啦中文免费视频观看日本| 女人久久www免费人成看片| 一级毛片我不卡| 两个人的视频大全免费| 成年免费大片在线观看| 日韩不卡一区二区三区视频在线| 国产成人精品一,二区| 久久精品综合一区二区三区| 久久久久久久久久成人| 六月丁香七月| 日韩欧美 国产精品| a级毛色黄片| 亚洲精品成人久久久久久| 伊人久久国产一区二区| 国产成人福利小说| 内射极品少妇av片p| 六月丁香七月| 欧美精品一区二区大全| 免费看av在线观看网站| 99热网站在线观看| 国产欧美日韩精品一区二区| 国产成人a∨麻豆精品| 欧美最新免费一区二区三区| 男人爽女人下面视频在线观看| 午夜福利在线在线| 成年av动漫网址| 亚洲精品色激情综合| 波多野结衣巨乳人妻| 日本爱情动作片www.在线观看| 一个人看的www免费观看视频| 2021少妇久久久久久久久久久| 日韩av在线免费看完整版不卡| 亚洲成人一二三区av| 九九久久精品国产亚洲av麻豆| 国产乱人偷精品视频| 一边亲一边摸免费视频| 国产伦在线观看视频一区| 亚洲高清免费不卡视频| 国产免费一级a男人的天堂| 久久99热6这里只有精品| 精品久久久久久久久亚洲| 国产亚洲av片在线观看秒播厂 | 国产色爽女视频免费观看| 国产亚洲午夜精品一区二区久久 | 日本午夜av视频| 中文在线观看免费www的网站| 床上黄色一级片| 国产精品99久久久久久久久| 欧美日韩国产mv在线观看视频 | 夜夜看夜夜爽夜夜摸| 最近的中文字幕免费完整| 日本色播在线视频| 久久久久久久亚洲中文字幕| 中国国产av一级| 亚洲经典国产精华液单| 在线免费十八禁| 国产 一区精品| 尤物成人国产欧美一区二区三区| 色综合站精品国产| 欧美不卡视频在线免费观看| 好男人在线观看高清免费视频| 2021天堂中文幕一二区在线观| 内地一区二区视频在线| 一边亲一边摸免费视频| 久久鲁丝午夜福利片| 毛片女人毛片| 亚洲乱码一区二区免费版| 亚洲国产精品国产精品| 国产亚洲午夜精品一区二区久久 | 丝瓜视频免费看黄片| 2018国产大陆天天弄谢| 少妇熟女aⅴ在线视频| 久久国内精品自在自线图片| 亚洲乱码一区二区免费版| 国产在视频线在精品| 亚洲精品国产成人久久av| 亚洲色图av天堂| 亚洲欧美中文字幕日韩二区| 欧美日韩精品成人综合77777| 久久久久久久亚洲中文字幕| 欧美一区二区亚洲| 亚洲电影在线观看av| 国产视频内射| 免费电影在线观看免费观看| 日日摸夜夜添夜夜添av毛片| 国语对白做爰xxxⅹ性视频网站| 婷婷色麻豆天堂久久| 亚洲不卡免费看| 国产一区二区亚洲精品在线观看| 看十八女毛片水多多多| 成人性生交大片免费视频hd| 欧美bdsm另类| 午夜免费观看性视频| 亚洲aⅴ乱码一区二区在线播放| 身体一侧抽搐| 91精品国产九色| xxx大片免费视频| 国产探花极品一区二区| 久久久久网色| 欧美bdsm另类| 国产伦精品一区二区三区四那| 91久久精品国产一区二区成人| 欧美精品一区二区大全| 亚洲人成网站在线播| 丰满人妻一区二区三区视频av| 黄色欧美视频在线观看| 精品久久国产蜜桃| 最近2019中文字幕mv第一页| 国产精品久久久久久久久免| 日本色播在线视频| 欧美成人精品欧美一级黄| 国产不卡一卡二| 午夜久久久久精精品| 亚洲精品色激情综合| 欧美97在线视频| av在线天堂中文字幕| 天天一区二区日本电影三级| 性色avwww在线观看| 午夜福利高清视频| 成人鲁丝片一二三区免费| 亚洲怡红院男人天堂| 亚洲三级黄色毛片| 亚洲av一区综合| 一个人看的www免费观看视频| 国产精品熟女久久久久浪| 亚洲精品久久久久久婷婷小说| 国产单亲对白刺激| 99热网站在线观看| 久久精品熟女亚洲av麻豆精品 | 91精品伊人久久大香线蕉| 你懂的网址亚洲精品在线观看| 熟女人妻精品中文字幕| 中文欧美无线码| 久久久久久久国产电影| www.色视频.com| 美女被艹到高潮喷水动态| 久久精品综合一区二区三区| 联通29元200g的流量卡| 色哟哟·www| 在线播放无遮挡| 日日干狠狠操夜夜爽| 免费av不卡在线播放| 午夜激情福利司机影院| 免费电影在线观看免费观看| 十八禁网站网址无遮挡 | 精品久久久久久久久亚洲| 国产伦精品一区二区三区四那| 国产精品久久视频播放| av播播在线观看一区| 精品一区二区三卡| 91精品伊人久久大香线蕉| 国产极品天堂在线| 免费黄频网站在线观看国产| 搡老妇女老女人老熟妇| 久久久久久九九精品二区国产| 伊人久久国产一区二区| 少妇裸体淫交视频免费看高清| 亚洲av日韩在线播放| 高清毛片免费看| 18禁在线播放成人免费| 亚洲精品aⅴ在线观看| 久久精品国产鲁丝片午夜精品| 久久综合国产亚洲精品| 欧美精品一区二区大全| 又爽又黄a免费视频| 亚洲精品久久午夜乱码| 91午夜精品亚洲一区二区三区| 免费黄频网站在线观看国产| 久久97久久精品| 2022亚洲国产成人精品| 日韩视频在线欧美| 性插视频无遮挡在线免费观看| 午夜福利在线在线| 精品一区二区三卡| 97热精品久久久久久| 亚洲国产精品sss在线观看| 国产淫语在线视频| videossex国产| 日本色播在线视频|