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

    地基偏振拉曼激光雷達與CALIPS0星載激光雷達的對比觀測研究

    2016-11-28 02:20:08黃忠偉張瑾超辛亞男
    干旱氣象 2016年5期
    關(guān)鍵詞:沙塵偏振氣溶膠

    溫 春,黃忠偉,周 天,張瑾超,3,辛亞男,4

    (1.蘭州大學大氣科學學院,甘肅 蘭州 730000;2.新疆輪臺縣氣象局,新疆 輪臺 841600;3.浙江三門縣氣象局,浙江 三門 317100;4.甘肅省金昌市氣象局,甘肅 金昌 737100)

    ?

    地基偏振拉曼激光雷達與CALIPS0星載激光雷達的對比觀測研究

    溫 春1,2,黃忠偉1,周 天1,張瑾超1,3,辛亞男1,4

    (1.蘭州大學大氣科學學院,甘肅 蘭州 730000;2.新疆輪臺縣氣象局,新疆 輪臺 841600;3.浙江三門縣氣象局,浙江 三門 317100;4.甘肅省金昌市氣象局,甘肅 金昌 737100)

    云和氣溶膠具有顯著的氣候效應,且存在很大的不確定性,所以獲取云和氣溶膠的垂直結(jié)構(gòu)信息對于研究全球氣候變化具有重要意義。利用蘭州大學半干旱氣候變化教育部重點實驗室自行研制的地基偏振拉曼激光雷達,將2014年3月6日和22日的觀測資料與CALIPSO星載激光雷達的觀測資料進行對比研究。結(jié)果表明:晴天條件下,CALIPSO星載激光雷達和地基激光雷達的反演結(jié)果總體上有較好的一致性,而2 km以下的差異主要是因為地基激光雷達受到幾何重疊因子的影響。因此,在晴空條件下綜合地基和星載激光雷達觀測數(shù)據(jù),可計算出地基激光雷達的幾何重疊因子和盲區(qū);多云條件下,CALIPSO星載激光雷達和地基激光雷達在中高層的觀測結(jié)果具有較好的一致性,但低層存在較大差異,原因主要可能是由于兩者觀測方向的不同和云體分布不均。

    偏振拉曼激光雷達;CALIPSO激光雷達;云和氣溶膠垂直結(jié)構(gòu);對比觀測

    引 言

    云和氣溶膠在地球輻射收支平衡中扮演著重要角色。一方面,云通過吸收和發(fā)射長波輻射及反射和散射短波輻射來調(diào)節(jié)地氣系統(tǒng)的能量收支平衡,另一方面,云通過降水參與大氣中的水循環(huán),進而調(diào)節(jié)大氣溫度、濕度。云的輻射作用不僅取決于其微物理特性(如有效粒子半徑、數(shù)濃度等),還涉及云高、云量和云的輻射特性等多方面因素,存在很大的不確定性[1]。由于云對大氣的熱力和動力影響,對其研究已成為當今氣候變化研究的熱點,了解云的垂直結(jié)構(gòu)分布特征有助于更好地理解云是如何影響地氣系統(tǒng)的。由于氣溶膠的直接和間接輻射效應,它已經(jīng)成為全球氣候變化研究的最大不確定因子[2],而氣溶膠的垂直分布是造成氣溶膠直接輻射強迫不確定性的主要原因之一[3],并且近地面的氣溶膠對人類健康也有重要影響。而系統(tǒng)性的氣溶膠全球垂直分布廓線觀測的缺少,對氣溶膠長波輻射效應和有云條件下的短波輻射效應的評估仍有相當大的不確定性[4],此外云和氣溶膠相互作用所產(chǎn)生的間接氣候效應(對云的反照率和云的壽命的作用)也與氣溶膠濃度的垂直分布以及云的相對密度相關(guān)[2],因此對大氣氣溶膠的研究已成為氣候變化研究的熱點。

    激光雷達(Lidar)是一種光學主動遙感儀器。自20世紀60年代問世以來,激光雷達技術(shù)發(fā)展迅速,應用領(lǐng)域越來越廣泛。激光雷達在大氣遙感與環(huán)境監(jiān)測研究方面越來越受歡迎,已發(fā)展成為地基、車載[5]和星載多平臺形式。國外在利用激光雷達對云的光學特性、垂直結(jié)構(gòu)以及云和氣溶膠相互作用等方面研究較深入[6-15],國內(nèi)利用激光雷達進行了很多云和氣溶膠的探測試驗及研究[16-30]。隨著星載激光雷達技術(shù)的快速發(fā)展,其在大氣遙感中的應用也越來越廣泛[31-35]。

    實際觀測應用中,不論是地基激光雷達還是星載激光雷達都存在各自的優(yōu)缺點,地基激光雷達可以定點長期觀測,具有很好的時間連續(xù)性,但空間分辨率較差;星載激光雷達采用衛(wèi)星作為平臺,運行軌道高,探測視野廣,可以進行全球范圍的觀測,剛好彌補了地基激光雷達在探測空間范圍的不足。因此,地基和星載激光雷達相結(jié)合的探測方式已經(jīng)成為大氣遙感研究的發(fā)展趨勢。

    本文利用蘭州大學自行研制的地基偏振拉曼激光雷達的觀測資料與CALIPSO星載激光雷達觀測資料,對比研究地基激光雷達與星載激光雷達CALIPSO的探測結(jié)果,進一步評估星載激光雷達探測蘭州地區(qū)云和氣溶膠垂直結(jié)構(gòu)的精度。

    1 儀器與數(shù)據(jù)

    1.1 地基偏振拉曼激光雷達系統(tǒng)

    如圖1所示,所使用的激光雷達系統(tǒng)由激光發(fā)射﹑信號接收﹑數(shù)據(jù)采集3部分組成,其基本原理是:激光器產(chǎn)生的脈沖激光經(jīng)過擴束準直處理后向大氣中發(fā)射,激光在大氣中傳輸時被大氣分子﹑云和氣溶膠所散射,然后用望遠鏡接收這部分回波信號,接著對信號進行一系列處理后,將光信號轉(zhuǎn)換成電信號傳輸至電腦顯示與保存。

    圖1 地基偏振拉曼激光雷達系統(tǒng)結(jié)構(gòu)

    1.1.1 激光發(fā)射子系統(tǒng)

    激光發(fā)射系統(tǒng)主要由激光器﹑擴束器﹑反射鏡組成。該系統(tǒng)所用的激光器為Nd:YAG脈沖激光器(Continuum Surelite),基波波段為1 064 nm,經(jīng)二倍頻和三倍頻后可輸出532 nm和355 nm的激光,單脈沖能量可達430 MJ,脈沖重復頻率為20 Hz,脈沖寬度為8.7 ns,表1為激光雷達的主要技術(shù)參數(shù)。擴束器的作用是使激光光斑變大且發(fā)散角變小。被擴束、準直后的激光經(jīng)由3個反射鏡,以同軸方式垂直發(fā)射到大氣中。

    表1 地基偏振拉曼激光雷達系統(tǒng)的主要技術(shù)參數(shù)

    1.1.2 信號接收子系統(tǒng)

    信號的接收和處理系統(tǒng)主要由望遠鏡、可變光闌、光纖、接收光學單元、探測器組成。

    接收望遠鏡用于接收大氣反射的回波信號,本文所用的是卡塞格林反射式接收望遠鏡,是目前在激光雷達方面應用最廣泛的接收望遠鏡;可變光闌一般被置于望遠鏡的焦點處,通過調(diào)節(jié)可變光闌來調(diào)節(jié)激光雷達系統(tǒng)的接收視場角,以此來減小背景噪聲和多次散射的影響;利用透鏡將望遠鏡接收的信號光最大限度地輸入到光纖中去,光纖的另一端與成像光譜儀連接,將光信號輸入到光譜儀中;接收光學單元包括濾光片﹑分束鏡﹑準直透鏡及偏振晶體;在激光雷達的接收系統(tǒng)中,光電倍增管(PMT)和雪崩光電二極管(APD)是目前使用最廣泛的探測器,本文所采用的是光電倍增管(PMT),其作用是將望遠鏡接收到的光信號轉(zhuǎn)換為電信號。

    1.1.3 數(shù)據(jù)采集系統(tǒng)

    數(shù)據(jù)采集設備是專門用于激光雷達信號探測的Licel瞬態(tài)記錄儀,可以同時以模數(shù)轉(zhuǎn)換方式(A/D)和光子計數(shù)方式(PC)對PMT的高速信號進行采集,具有很高的動態(tài)范圍,且空間分辨率較高。它的主要部件包括可以實時進行信號平均的高速模數(shù)轉(zhuǎn)換器、單光子探測的甄別器、多通道記錄儀以及前置放大器。

    1.2 CALIPSO星載激光雷達

    為了研究云和氣溶膠在調(diào)節(jié)地球氣候中的作用及兩者之間的相互作用,1998年美國國家宇航局(NASA)與法國國家航天中心(CNES)合作實施“云—氣溶膠激光雷達和紅外探測者衛(wèi)星觀測”(Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations, CALIPSO)計劃,并于2008年4月28日成功發(fā)射了CALIPSO衛(wèi)星。CALIPSO衛(wèi)星主要載荷包括:正交偏振云—氣溶膠激光雷達(CALIOP)﹑寬視角相機(WFC)﹑紅外成像輻射計(IIR)。CALIOP采用了偏振檢測技術(shù),實現(xiàn)了全球覆蓋,它是世界上首個應用型的星載云和氣溶膠激光雷達,具有3個接收通道,一個用于接收波長為1 064 nm的后向散射信號,另外2個分別接收波長為532 nm的正交偏振后向散射信號。它的探測原理繼承了目前地基和空基激光雷達的探測原理,下面具體介紹CALIOP的發(fā)射和接收系統(tǒng),其參數(shù)見表2。

    表2 CALIOP系統(tǒng)主要參數(shù)

    CALIOP的發(fā)射系統(tǒng)有2套完全相同的Nd:YAG型激光器,單脈沖發(fā)射的最大能量為440 MJ,重復頻率為20.16 Hz。為了提高衛(wèi)星的可靠性和使用壽命,實際工作時單脈沖輸出的最大能量為220 MJ,再使用倍頻技術(shù)獲得波長為532 nm的激光輸出脈沖,最終輸出單脈沖能量變?yōu)?10 MJ的532 nm和1 064 nm激光,激光經(jīng)過擴束器并通過光束準直系統(tǒng)進入大氣。激光器的散射角為100 μrad,而視場校準結(jié)構(gòu)可以對發(fā)射系統(tǒng)視場和接收系統(tǒng)視場進行校準,其調(diào)節(jié)范圍為±1°,分辨率為16 μrad。

    CALIOP的接收系統(tǒng)包括望遠鏡、接收光學單元、探測器、信號放大器。CALIOP接收系統(tǒng)采用直徑為1 m、視場角為130 μrad的望遠鏡。望遠鏡系統(tǒng)接收的散射激光雷達脈沖經(jīng)過聚焦后通過濾鏡,然后準直入射到分光鏡上,532 nm的激光脈沖發(fā)射光,首先通過干涉濾光片和窄帶濾光片減小太陽背景光的影響,然后由偏振分光鏡分為532 nm平行光和垂直光;1 064 nm的激光脈沖發(fā)射光,經(jīng)過干涉分光鏡后直接由探測器探測。用于檢測532 nm的脈沖信號的探測器為光電倍增管(PMT)。由于PMT在1 064 nm處的量子效率較低,所以系統(tǒng)采用雪崩光電二極管(APD)來檢測1 064 nm的激光脈沖回波信號。

    1.3 數(shù) 據(jù)

    用于與CALIPSO星載激光雷達進行觀測對比的地基激光雷達系統(tǒng)(以下簡稱LZU- lidar)位于蘭州大學校園內(nèi)(36.05°N,103.85°E),此激光雷達系統(tǒng)可同時探測355 nm、532 nm垂直和平行通道的后向散射信號,觀測時間為2014年3月6日和22日。

    星載激光雷達的數(shù)據(jù)為CALIPSO提供的激光雷達Level 1B數(shù)據(jù)產(chǎn)品,是校準后的半軌道(白天或者夜間)激光雷達單射(兩單射點間距333 m)剖面數(shù)據(jù),包括532 nm和1 064 nm通道衰減后向散射數(shù)據(jù)和532 nm通道的退偏比數(shù)據(jù)等。使用的數(shù)據(jù)參數(shù)為:衰減后向散射剖面數(shù)據(jù)(包括532 nm總衰減后向散射系數(shù)、532 nm垂直通道衰減后向散射系數(shù)和1 064 nm衰減后向散射系數(shù))、方位及高度信息(包括經(jīng)緯度、激光雷達數(shù)據(jù)高度)。CALIPSO星載激光雷達的軌道周期為16 d,本文選取以LZU-lidar為中心,半徑為100 km范圍內(nèi)的CALIPSO星載激光雷達level 1數(shù)據(jù),其中2014年3月6日CALIPSO衛(wèi)星地面軌道上距LZU lidar最近點的距離為34.7107 km,該點緯度為36.1350°N;22日CALIPSO衛(wèi)星地面軌道上距LZU lidar最近點的距離為38.3531 km,該點的緯度為36.1431°N。這2條衛(wèi)星軌道過境時間均為03:23(北京時,下同),3月6日為多云天氣,22日為晴空有沙塵天氣。進行對比之前,先將LZU-lidar觀測數(shù)據(jù)在CALIPSO衛(wèi)星過境時間點前后各10 min的廓線進行平均,與之對應,將CALIPSO星載激光雷達以最近點為中心前后各100條廓線進行平均。

    2 方法介紹

    2.1 地基激光雷達幾何重疊因子的確定

    激光雷達系統(tǒng)的接收視場角理論上應略大于發(fā)射視場角,這樣回波信號才可能完全進入接收望遠鏡。然而,由于發(fā)射系統(tǒng)與接收系統(tǒng)不同軸的影響,導致回波信號在某段距離上有所損失,望遠鏡只能接收到部分信號,在近地距離尤為明顯,由此給觀測帶來誤差,必須對其進行訂正,而這一訂正因子就是幾何重疊因子。在激光雷達的實際操作中,即使調(diào)節(jié)時產(chǎn)生較小的誤差,也會給氣溶膠消光系數(shù)帶來很大誤差[36],因此幾何重疊因子的確定很重要。目前普遍用于幾何重疊因子的確定方法有實驗測量、幾何推導[37]以及其它新方法[38-39]。

    然而,對于本文所使用的激光雷達系統(tǒng)來說,采取實驗法求解幾何重疊因子比較困難,原因是該系統(tǒng)固定垂直發(fā)射,無法進行水平探測,以及該系統(tǒng)的幾何光學參數(shù)未知。

    對地基激光雷達來說,激光雷達方程為:

    式中:P0為t0時刻發(fā)射的能量,C為雷達系統(tǒng)參數(shù),β為r范圍內(nèi)的后向散射系數(shù),σ為高度r處的消光系數(shù),r為激光雷達到測量點的距離,r′是積分距離。

    對CALIPSO星載激光雷達來說,激光雷達方程為:

    式中:r為激光雷達到測量點的距離,P為接收到的測量信號,ξ為激光雷達系統(tǒng)參數(shù),E0為單束或多束平均激光能量,β為r范圍內(nèi)的后向散射系數(shù),σ為高度r處的消光系數(shù),r′是積分距離。

    將地基激光雷達和CALIPSO星載激光雷達的資料進行對比,計算地基激光雷達系統(tǒng)的幾何重疊因子。由于CALIPSO星載激光雷達和LZU-lidar觀測方向的區(qū)別,在近地面,CALIPSO星載激光雷達對信號基本100%接收,而地基激光雷達在遠地面對信號基本100%接收,根據(jù)這一特征,將地基激光雷達和星載激光雷達的雷達方程相比較,使用CALIPSO星載激光雷達和LZU-lidar的對比觀測資料,計算LZU-lidar的幾何重疊因子,計算方法如下:

    2.2 退偏比、色比的計算

    偏振激光雷達系統(tǒng)發(fā)射平行于激光電場強度振動方向的平行分量和垂直于該方向的垂直分量,根據(jù)米散射理論,對于球形粒子而言,其散射光仍然保持入射光的偏振態(tài)。對于非球形粒子,其散射光的振動方向相比于入射光發(fā)生轉(zhuǎn)動,即產(chǎn)生了退偏振。例如,對于硫酸鹽氣溶膠和水云等球形粒子的偏振垂直分量為0,而對于沙塵氣溶膠和冰云等非球形粒子的偏振垂直分量信號較強,因此退偏比定義為:

    式中,P⊥為偏振垂直分量,P‖為偏振平行分量。

    退偏比的大小不僅與顆粒物的形狀有關(guān),還與顆粒物的大小和濃度有關(guān)。此外,星載激光雷達系統(tǒng)可同時發(fā)射1 064 nm和532 nm波段的脈沖,接收這2個波段的回波信號可得到另一個觀測量,即色比,其定義為2個波段(1 064 nm和532 nm)的回波信號之比,用來表征顆粒物的大小,色比越大顆粒物越大。綜上所述,退偏比和色比可以提供更多有關(guān)氣溶膠和云的信息。

    3 結(jié)果分析

    3.1 晴天條件下觀測個例的廓線對比

    圖2是2014年3月22日CALIPSO星載激光雷達和LZU-lidar的廓線對比。從532 nm總后向散射衰減廓線(圖2a)看出,1.5 km以上,兩條廓線具有很好的一致性,1.5 km以下兩者的差異是受幾何重疊因子的影響。LZU-lidar和CALIPSO星載激光雷達在2 km內(nèi)均探測到沙塵層,但兩者觀測的沙塵層高度和強度有一定差異,2 km以下的差異可能是由于沙塵天氣時氣溶膠水平分布非均一和沙塵天氣的雙重影響。從退偏比(圖2b)廓線看,即使在晴天條件下,CALIPSO星載激光雷達也能探測到大氣背景氣溶膠的變化,而LZU-lidar則不能;對于2 km以下的沙塵氣溶膠,兩者的退偏比都要大于背景氣溶膠的退偏比,相比而言,CALIPSO得到的沙塵的退偏比更為合理。從色比廓線(圖2c)看,晴天條件下,LZU-lidar對于大氣背景氣溶膠顆粒的大小反映沒有CALIPSO星載激光雷達明顯,對于2 km以下的沙塵,CALIPSO有比較明顯的反映,雖然LZU-lidar也有反映,相比而言,CALIPSO對于沙塵的探測效果更為顯著。

    圖3給出CALIPSO星載激光雷達總后向散射衰減(532 nm)、退偏比(532 nm)和色比在不同經(jīng)緯度上的剖面,圖4是LZU-lidar總后向散射衰減(532 nm)、退偏比(532 nm)和色比的時間—高度剖面。對比圖3a和圖4a可以看出,即使晴天,CALIPSO星載激光雷達的532 nm總后向散射衰減在7 km以下有較明顯的變化,對應于圖2a的差異是合理的,特別是對2 km以下的沙塵,兩者都能探測到,但CALIPSO觀測到的強度明顯強于LZU-lidar,因此CALIPSO對于沙塵的觀測效果更好;7 km以上,由于大氣趨于干凈,兩者都無明顯變化,因此圖2a中7 km以上2條廓線具有很好的一致性。其次,對比圖3b和圖4b,CALIPSO的退偏比在晴天條件下對于大氣背景氣溶膠也有明顯反映,而LZU-lidar則基本無明顯變化,對于2 km以下的沙塵氣溶膠,CALIPSO的退偏比對沙塵的反映更為合理,因此,圖2b的廓線也是合理的。對比圖3c和圖4c,CALIPSO對2 km以下的沙塵反映比較明顯,對于2 km以上區(qū)域則變化不大,LZU-lidar也是如此,這與圖2c中的廓線基本一致。

    圖2 2014年3月22日CALIPSO星載激光雷達和LZU-lidar的廓線對比(a)532 nm總后向散射衰減,(b)532 nm退偏比,(c)色比

    圖4 2014年3月22日LZU-lidar 532 nm總后向散射衰減(a,單位:km-1·sr-1)、532 nm退偏比(b)和色比(c)的時間—高度剖面

    由圖2a可知,CALIPSO星載激光雷達和LZU-lidar對2 km以下的氣溶膠觀測差異較大,而造成LZU-lidar探測差異較大的原因,一方面由于高度越低,對應的幾何重疊因子越??;另一方面是由于近地面光信號和電子學系統(tǒng)干擾對最初幾個脈沖帶來較大的噪聲引起的[40]。可見幾何重疊因子對于氣溶膠觀測的影響還是很大的,特別是對于氣溶膠消光系數(shù)反演的影響。

    圖5是計算得到的LZU-lidar的幾何重疊因子,可以看出100~700 m之間<0.1,即說明在此探測范圍不到10%的回波信號被儀器所接收,從700 m開始幾何重疊因子快速增大,一直到1.38 km變?yōu)?,即激光雷達的發(fā)射場和接收場能夠很好地重疊。此外通過計算還得到LZU-lidar的盲區(qū)高度為120 m。

    圖5 計算得到的LZU-lidar幾何重疊因子Fig.5 The calculated geometric overlap factor of LZU-lidar

    3.2 有云條件下觀測個例的廓線對比

    圖6是2014年3月6日有云條件下CALIPSO星載激光雷達和LZU-lidar的廓線對比。從532 nm總后向散射衰減廓線(圖6a)看,兩者都可以探測到9 km左右的云,CALIPSO只顯示為一層云,而LZU-lidar卻探測到兩層云;5 km上下的云兩者都有反映,LZU-lidar顯示為多層云,而CALIPSO星載激光雷達對此反映卻不明顯;LZU-lidar探測到3 km處的云,而CALIPSO并未觀測到;在2 km內(nèi)的氣溶膠層,CALIPSO反映得更為明顯,而LZU-lidar在近地面由于受到幾何重疊因子的影響,對氣溶膠觀測效果不顯著。

    從CALIPSO星載激光雷達和LZU-lidar的退偏比(圖6b)看,2 km以上各對應位置,CALIPSO星載激光雷達的退偏比基本都大于LZU-lidar。在9 km高度有多層冰云的位置,LZU-lidar的退偏比在0.2左右,而CALIPSO在0.4左右;5 km以下有水云的位置,LZU-lidar的退偏比基本都在0.1左右,而CALIPSO在0.3左右。一般來說水云屬于球形粒子,而實際觀測結(jié)果是水云的退偏比較大,這是由于激光在云中的多次散射效應所致。CALIPSO星載激光雷達探測的水云和冰云的退偏比都很大,而LZU-lidar探測得到的兩者退偏比都比較小,因此星載和地基激光雷達在區(qū)分云相上有一定困難[41]。對于1~2 km的氣溶膠層,LZU-lidar探測的退偏比也要大于CALIPSO探測的值,但氣溶膠的退偏比要小于云的退偏比,由此可見退偏比的大小不僅與顆粒物的形狀有關(guān),還與顆粒物的大小和濃度有關(guān)。圖6c是CALIPSO星載激光雷達和LZU-lidar的色比對比,可以看出6 km以下,LZU-lidar的色比對于顆粒物大小的變化反映顯著,最大值可達2;而CALIPSO探測的6 km以下和9 km左右云的色比相差不大,說明僅根據(jù)云粒徑的大小區(qū)分云相也是有困難的。對于2 km以下區(qū)域的氣溶膠層,一般來說,對于同一套激光雷達系統(tǒng),云的色比要大于氣溶膠的,因為云粒徑大小遠大于氣溶膠粒徑,而圖6c的反映是合理的。

    圖7給出CALIPSO星載激光雷達探測的532 nm后向散射衰減、退偏比和色比在不同經(jīng)緯度上的剖面,圖8是LZU-lidar探測的532 nm后向散射衰減、退偏比和色比的時間—高度剖面。CALIPSO的觀測是在當?shù)睾0胃叨纫陨希虼藢Ρ葧r要減去當?shù)睾0胃叨取?/p>

    從圖7a看出,大約9 km高度有一層水平尺度達400 km的高云,云層厚度近2 km;5 km上下存在水平尺度約100 km的云。云頂高度以及云層厚度隨著CALIPSO衛(wèi)星地面軌道的改變而快速變化,并且LZU-lidar距離CALIPSO地面軌道最近也有36 km,所以圖6a中云層高度的差異也是合理的,衛(wèi)星過境時LZU-lidar觀測到多層云也是合理的,對應于圖8a也可看出上述結(jié)果。至于兩者觀測到云層厚度的差異,可能是由于地基和星載激光雷達探測方向的不同,地基激光雷達是垂直于地面向上進行觀測,而星載激光雷達是從高空向地面向下進行觀測,地基激光雷達在低空激光能量衰減小,因此接收到的信號要強,而CALIPSO星載激光雷達所發(fā)射的激光到達低空時已衰減很厲害,而且后向散射從低空到高空,要穿過多層云,能量進一步衰減,所以CALIPSO接收到的信號就比較弱。對比圖7b和圖8b,可以看出,LZU-lidar和CALIPSO星載激光雷達探測的退偏比在9 km和5 km以下有云區(qū)域差異不大,因此圖6b也是合理的,由此可見,只是根據(jù)顆粒物的退偏比很難區(qū)分云的相態(tài),還應結(jié)合顆粒物的大小及濃度來加以區(qū)分。

    對比圖7c和圖8c,可以看出,5 km以下云和氣溶膠的色比有明顯差異,對應于圖6c中LZU-lidar探測的色比在5 km以下的明顯變化。而在9 km以上,由于衛(wèi)星過境時,LZU-lidar觀測到的云快速變化,使得此處色比較小,對應于圖6c中高云的色比較?。粡膱D7c中可以看出色比在上述兩區(qū)域無明顯差異,圖6c也反映了CALIPSO在這兩區(qū)域探測的色比差異不大。

    圖6 2014年3月6日CALIPSO星載激光雷達和LZU-lidar的廓線對比(a)532 nm總后向散射衰減,(b)532 nm退偏比,(c)色比

    圖7 2014年3月6日CALIPSO星載激光雷達532 nm總后向散射衰減(a,單位:km-1·sr-1),532 nm退偏比(b)及色比(c)

    圖8 2014年3月6日LZU-lidar 的總后向散射衰減(532 nm)(a,單位:km-1·sr-1)、退偏比(532 nm)(b)和色比(c)的時間—高度分布

    4 結(jié)論與討論

    (1)晴天條件下,CALIPSO星載激光雷達和LZU-lidar的反演結(jié)果有較好的一致性,說明LZU-lidar的探測能力比較好,反演結(jié)果具有較高的可信度,而1.5 km以下兩者的差異主要是因為LZU-lidar受到幾何重疊因子的影響。

    (2)有云條件下,CALIPSO星載激光雷達和LZU-lidar在中高層的反演結(jié)果具有較好的一致性,而在低層,由于兩者觀測方向不同結(jié)果存在較大差異,說明在多云條件下,LZU-lidar的反演結(jié)果更為可信,CALIPSO星載激光雷達和LZU-lidar對于云相態(tài)的區(qū)分都有一定困難。

    (3)1.38 km以上,LZU-lidar系統(tǒng)的發(fā)射場和接收場能夠很好地重疊,LZU-lidar盲區(qū)高度為120 m。對于不能采用水平探測確定幾何重疊因子的激光雷達系統(tǒng),本文所使用的計算方法簡便實用、行之有效。

    為了更進一步評估CALIPSO星載激光雷達對蘭州地區(qū)云和氣溶膠的觀測效果和利用星載激光雷達來訂正地基激光雷達對氣溶膠的觀測,還應結(jié)合SACOL觀測站現(xiàn)有的可靠儀器來進行對比,此外為了更好地利用CALIPSO的觀測資料,應當對不同區(qū)域的觀測資料與可靠的地面儀器進行比較,并給出具有區(qū)域特色的激光雷達比。

    致謝:感謝蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)提供地基激光雷達的觀測數(shù)據(jù)及美國國家航空航天局CALIPSO科學組提供CALIPSO激光雷達觀測數(shù)據(jù)

    [1] Huang J P, Minnis P, Lin B, et al. Advanced Retrievals of Multilayered Cloud Properties using MuhispectralMeasurements[J]. Journalof Geophysical Research, 2005,110(D15), DOI:10.1029/2004JD005101.

    [2] 王標. 氣候變化研究中的大氣垂直結(jié)構(gòu)探測述評[J]. 大氣與環(huán)境光學學報,2011,6(1):2-10.

    [3] Meloni D, Sarra A D, Iorio T D, et al. Influence of the vertical profile of Saharan dust on the visible direct radiative forcing[J]. Journal of Quantitative Spectroscopyand Radiative Transfer, 2005,93(4):397-413.

    [4] Yu H, Kaufman Y J, Chin M, et al. A review of measurment-based assessments of the aerosol direct radiative effect and forcing[J]. Atmospheric Chemistry Physics, 2006,6(11):613-666.

    [5] 謝晨波,周軍,岳古明,等. 新型車載式拉曼激光雷達測量對流層水汽[J]. 光學學報,2006,26(9):1281-1286.

    [6] Charlson R J, Schwartz S E, Hales J M, et al. Climate forcing by anthropogenic aerosols[J]. Science, 1992,255(5043):423-430.

    [7] Shimian H, Sasano Y, Nakane H, et al. Large scale laser radar for measuring aerosol distribution over a wide area[J]. Applied Optics,1985, 24(5):617-626.

    [8] Qing P, Nakane H, Sasano Y, et al. Numerical simulation of the retrieval of aerosol size distribution from multiwavelength laser radar measurements[J]. Applied Optics, 1989,28(24):5259-5265.

    [9] Sasano Y, Nakane H. Significance of the extinction/backscatter ratio and the boundary value term in the solution for the two-component lidar equation[J]. Applied Optics, 1984,23(23):11-13.

    [10] Sasano Y, Browell E, Ismail S. Error caused by using a constant extintion/backscatteringratio in the lidar solution[J]. Applied Optics,1985,24(22):3929-3932.

    [11] Sasano Y, Nakane H. Quantitative analysis of RHI lidar data by an interative adjustment of the boundary condition term in the lidar solution[J]. Applied Optics, 1987,26(4):615-616.

    [12] Takamura T, Sasano Y. Ratio aerosol backscatter to extintion coefficients as determined from angular scattering measurements for use in atmospheric lidar application[J]. Optical Quantum Electronics,1987,19(5):293-302.

    [13] Takamura T, Sasano Y. Aerosol optical properties inferred from simultaneous lidar, aerosol-counter, and sunphotometer measurements[J]. Journal of the Meteorological Society of Japan,1990,68(6):729-739.

    [14] Takamura T, Sasano Y, Hayasaka T. Tropospheric aerosol optical properties derived from lidar, sun photometer, and optical particle counter measurements[J]. Applied Optics,1994,33(30):7132-7140.

    [15] Sasano Y. Tropospheric aerosol extinction coefficient profiles derived from scanning lidar measurements over Tsukuba, Japan, from 1990 to 1993[J]. Applied Optics, 1996,35(24):4941-4952.

    [16] 白宇波,石廣玉,田村耕一,等. 拉薩上空大氣氣溶膠光學特性的激光雷達探測[J]. 大氣科學,2000,24(4):559-567.

    [17] 黃建平,黃忠偉. 黃土高原大氣氣溶膠垂直結(jié)構(gòu)的微脈沖激光雷達測量[J]. 大氣和海洋科學通訊,2008(1):8-11.

    [18] 張文煜,宋嘉堯,胡交超,等. 激光雷達對半干旱地區(qū)云層的觀測[J]. 蘭州大學學報(自然科學版),2009,45(5):53-56.DOI:0455-2059(2009)05-0053-04

    [19] 鄧濤,張鐳,吳兌. 蘭州地區(qū)高云和氣溶膠光學特性及其輻射效應[J]. 高原氣象,2010,29(1):230-235.

    [20] 曹賢潔,張鐳,周碧,等. 利用激光雷達觀測蘭州沙塵氣溶膠輻射特性[J]. 高原氣象,2009,28(5):1115-1120.

    [21] 周碧,張鐳,曹賢潔,等. 利用激光雷達資料分析蘭州遠郊氣溶膠光學特性[J]. 高原氣象,2011,30(4):1011-1017.

    [22] 周天,黃忠偉,黃建平,等. 黃土高原地區(qū)云垂直結(jié)構(gòu)的激光雷達遙感研究[J]. 干旱氣象,2013,31(2):246-253.

    [23] 劉瑞金,張鐳,王宏斌,等. 半干旱地區(qū)卷云特征的激光雷達探測[J]. 大氣科學,2011,35(5):863-870.

    [24] 夏俊榮,張鐳. Mie散射激光雷達探測大氣氣溶膠的進展[J]. 干旱氣象,2006,24(4):68-81.

    [25] 張杰,唐從國. 干旱區(qū)一次春季沙塵過程的大氣氣溶膠垂直分布結(jié)構(gòu)及其特征[J]. 高原氣象,2012,31(1):156-166.

    [26] 胡向軍,張鐳,鄭飛,等. 基于激光雷達資料的氣溶膠輻射效應研究[J]. 干旱氣象,2014,32(2):248-255.

    [27] 李紅,馬媛媛,楊毅. 基于激光雷達資料的小波變換法反演邊界層高度的方法[J]. 干旱氣象,2015,33(1):78-97.

    [28] Huang Z W, Huang J P, Bi J R, et al. Dust aerosol vertical structure measurements using three MPL lidars during 2008 China-U.S. joint dust field experiment[J]. Journal of Geophysical Research Atomospheres, 2010,115(D7):1307-1314.

    [29] Huang Z W, Huang J P, Hayasaka T, et al. Short-cut transport path for Asian dust directly to the Arctic: a case study[J]. Environment Research Letters, 2015,10(11), DOI:10.1088/1748-9326/10/11/114018.

    [30] Zhou T, Huang J P, Huang Z W, et al. The depolarization-attenuated backscatter relationship for dust plumes[J]. Optics Express, 2013,21(13):15195-15204.

    [31] 丁曉東,黃建平,李積明,等. 基于主動衛(wèi)星遙感研究西北地區(qū)云層垂直結(jié)構(gòu)特征及其對人工增雨的影響[J]. 干旱氣象,2012,30(4):529-538.

    [32] 張寅超,李丹,陳炳龍,等. 基于CALIPSO和轉(zhuǎn)動拉曼—米散射激光雷達研究北京地區(qū)大氣氣溶膠特性[J]. 北京理工大學學報,2013,33(6):634-643.

    [33] 徐婷婷,秦艷,陳勇航,等. 環(huán)上海地區(qū)干霾氣溶膠垂直分布的季節(jié)變化特征[J]. 環(huán)境科學,2012,33(7):2165-2171.

    [34] 申莉莉,盛立芳,陳靜靜. 一次強沙塵暴過程中沙塵氣溶膠空間分布的初步分析[J]. 中國沙漠,2010,30(6):1483-1490.

    [35] 李積明,黃建平,衣育紅,等. 利用星載激光雷達資料研究東亞地區(qū)云垂直分布的統(tǒng)計特征[J]. 大氣科學,2009,33(4):698-707.

    [36] 張改霞,張寅超,陶宗明,等. 激光雷達幾何重疊因子及其對氣溶膠探測的影響[J]. 量子電子學報,2005,22(2):299-304.

    [37] 戴麗莉,丁紅星,薛國剛,等. 雙軸Mie散射激光雷達及其幾何重疊因子的計算[J]. 激光雜志,2007,28(5):49-50.

    [38] 曹開法,杜小勇,汪少林,等. 純轉(zhuǎn)動拉曼法確定激光雷達幾何因子[J]. 光學學報,2007,27(8):1349-1354.

    [39] 閆順生,李海花. 利用喇曼-米方法確定激光雷達幾何因子[J]. 激光技術(shù),2013,37(4):511-514.

    [40] 章正宇,眭曉林. 激光測距弱信號數(shù)字相關(guān)檢測技術(shù)的研究與仿真[J]. 中國激光,2002,29(7):661-665.

    [41] Hu Y, Vaughan M, Liu Z, et al. The Depolarization-attenuated backscatterrelation:CALIPSOlidarmeasurements vs.theory[J].Optics Express, 2007,15(9),5327-5332.

    Contrast Observation of Spaceborne CALIPSO Lidar and a Ground-based Polarized Raman Lidar

    WEN Chun1,2, HUANG Zhongwei1, ZHOU Tian1,ZHANG Jinchao1,3, XIN Yanan1,4

    (1.CollegeofAtmosphericScience,LanzhouUniversity,Lanzhou730000China;2.LuntaiMeteorologicalBureauofXinjiang,Luntai841600,China;3.SanmenMeteorologicalStationofZhejiangProvince,Sanmen317100,China;4.JinchangMeteorologicalBureauofGansuProvince,Jinchang737100,China)

    Clouds and aerosols have significant impact on globe climate change, and due to lack of information about clouds and aerosols vertical structure, there is large uncertainty to quantify their climatic effects. Lidar has been well known as a unique remote sensing for measuring the vertical structure of clouds and aerosols with high spatio-temporal resolution. To validate observations of spaceborne lidar, measurements from a developed ground-based polarized Raman lidar in Lanzhou on March 6 and March 22, 2014 were conducted. The results show that observations from CALIPSO lidar and ground-based lidar had a good consistency for clear day except for signals near ground surface because of overlap effect of ground-based lidar. So the geometric overlap factor of ground-based lidar could be calculated from the comparison. For cloudy condition, CALIPSO lidar observations had a good agreement with those of ground-based lidar for high cloud layer, but due to different direction of observation for ground-based and spaceborne lidar or inhomogeneous of cloud in the lower atmosphere, they were different.

    polarized raman lidar; CALIPSO lidar; vertical structure of clouds and aerosols; comparative observations

    10.11755/j.issn.1006-7639(2016)-05-0779

    2016-03-01;改回日期:2016-05-11

    國家自然科學基金項目(41375031、41305027、41575017)和中央高?;究蒲袠I(yè)務費(lzujbky-2015-k01)共同資助

    溫春(1990-),男,新疆庫爾勒市人,助理工程師,主要從事天氣預報工作. E-mail:wenchun622@163.com

    黃忠偉,博士,主要從事激光雷達大氣遙感與環(huán)境監(jiān)測的研究. E-mail:huangzhongwei@lzu.edu.cn

    1006-7639(2016)-05-0779-10 DOI:10.11755/j.issn.1006-7639(2016)-05-0779

    P407.5

    A

    溫 春,黃忠偉,周 天,等.地基偏振拉曼激光雷達與CALIPS0星載激光雷達的對比觀測研究[J].干旱氣象,2016,34(5):779-788, [WEN Chun, HUANG Zhongwei, ZHOU Tian, et al. Contrast Observation of Spaceborne CALIPSO Lidar and a Ground-based Polarized Raman Lidar[J]. Journal of Arid Meteorology, 2016, 34(5):779-788],

    猜你喜歡
    沙塵偏振氣溶膠
    什么是沙塵天氣?
    國土綠化(2024年3期)2024-04-17 01:02:22
    氣溶膠傳播之謎
    偏振糾纏雙光子態(tài)的糾纏特性分析
    電子制作(2019年12期)2019-07-16 08:45:20
    氣溶膠中210Po測定的不確定度評定
    基于LabVIEW的偏振調(diào)制激光測距測量控制系統(tǒng)
    可怕的沙塵天氣
    學生天地(2017年12期)2017-05-17 05:51:25
    偏振旋轉(zhuǎn)效應在全光緩存器中的應用與實現(xiàn)
    四川盆地秋季氣溶膠與云的相關(guān)分析
    大氣氣溶膠成核監(jiān)測
    光的偏振特性在機動車安全駕駛中的應用
    物理與工程(2011年6期)2011-03-25 10:04:13
    别揉我奶头 嗯啊视频| 91久久精品电影网| 美女黄网站色视频| 美女cb高潮喷水在线观看| 嘟嘟电影网在线观看| av专区在线播放| 永久网站在线| 国产亚洲91精品色在线| 久久久久网色| 丝袜喷水一区| 老司机影院毛片| 少妇丰满av| 一级二级三级毛片免费看| 亚洲人与动物交配视频| 亚洲精华国产精华液的使用体验| 国产av在哪里看| 国产高清视频在线观看网站| 日韩亚洲欧美综合| 亚洲美女搞黄在线观看| 白带黄色成豆腐渣| 99久久精品热视频| 能在线免费看毛片的网站| 久久久久久久国产电影| 亚洲最大成人中文| 久久精品夜夜夜夜夜久久蜜豆| 日本免费a在线| 久久精品夜色国产| 亚洲欧美成人精品一区二区| 一级毛片电影观看 | 日本一二三区视频观看| 国产成人a区在线观看| 精品国产露脸久久av麻豆 | 国产精品人妻久久久影院| 国模一区二区三区四区视频| 久久久a久久爽久久v久久| 欧美日韩综合久久久久久| 在线观看一区二区三区| 天天躁日日操中文字幕| 国产在视频线在精品| 小说图片视频综合网站| 亚洲精品乱久久久久久| 99在线视频只有这里精品首页| 观看免费一级毛片| 成人毛片a级毛片在线播放| 七月丁香在线播放| 亚洲国产精品成人久久小说| av专区在线播放| 国产爱豆传媒在线观看| 久久久精品欧美日韩精品| 最近2019中文字幕mv第一页| 精品久久久久久久人妻蜜臀av| 国内少妇人妻偷人精品xxx网站| 亚洲熟妇中文字幕五十中出| 乱系列少妇在线播放| 中文字幕久久专区| 国产三级中文精品| 激情 狠狠 欧美| 国产老妇伦熟女老妇高清| 夜夜爽夜夜爽视频| 亚洲精品亚洲一区二区| 少妇裸体淫交视频免费看高清| a级毛色黄片| 国产亚洲精品av在线| 波多野结衣高清无吗| 男女国产视频网站| 韩国av在线不卡| 91精品一卡2卡3卡4卡| 欧美又色又爽又黄视频| 国产高清有码在线观看视频| 亚洲精品一区蜜桃| 国内精品美女久久久久久| 国产男人的电影天堂91| 非洲黑人性xxxx精品又粗又长| 小蜜桃在线观看免费完整版高清| 久久国产乱子免费精品| 亚洲aⅴ乱码一区二区在线播放| 欧美极品一区二区三区四区| 国产老妇女一区| 国产精品一及| 精品人妻熟女av久视频| av专区在线播放| 又粗又爽又猛毛片免费看| 色综合站精品国产| 国产乱来视频区| 中国美白少妇内射xxxbb| 欧美日本亚洲视频在线播放| av在线播放精品| 欧美极品一区二区三区四区| 精品熟女少妇av免费看| 国产亚洲最大av| 建设人人有责人人尽责人人享有的 | 日本-黄色视频高清免费观看| 亚洲精品乱久久久久久| 久久久久性生活片| 国产精品.久久久| 日韩中字成人| 午夜福利视频1000在线观看| 亚洲人成网站在线观看播放| 精华霜和精华液先用哪个| 日日摸夜夜添夜夜爱| 亚洲国产精品久久男人天堂| 九九在线视频观看精品| 欧美一区二区国产精品久久精品| 久热久热在线精品观看| 欧美日韩一区二区视频在线观看视频在线 | 美女被艹到高潮喷水动态| 日韩人妻高清精品专区| 一区二区三区四区激情视频| 99热这里只有是精品50| 自拍偷自拍亚洲精品老妇| 91在线精品国自产拍蜜月| 午夜精品国产一区二区电影 | 国产 一区 欧美 日韩| 日本黄色视频三级网站网址| 91久久精品国产一区二区三区| 99热这里只有是精品在线观看| 久久精品国产亚洲av涩爱| 亚洲精品成人久久久久久| 99久久无色码亚洲精品果冻| 啦啦啦韩国在线观看视频| 午夜福利在线在线| 美女内射精品一级片tv| 亚洲av成人精品一区久久| 欧美精品一区二区大全| 国产高清国产精品国产三级 | 亚洲自偷自拍三级| 国产精品美女特级片免费视频播放器| 看免费成人av毛片| 天堂√8在线中文| 秋霞在线观看毛片| 波野结衣二区三区在线| 亚洲av成人av| 蜜桃亚洲精品一区二区三区| 欧美最新免费一区二区三区| 丰满乱子伦码专区| 内地一区二区视频在线| 中文字幕av成人在线电影| 国产精品女同一区二区软件| 精品人妻一区二区三区麻豆| 午夜福利在线在线| 亚洲综合精品二区| 伦精品一区二区三区| 搡女人真爽免费视频火全软件| 男插女下体视频免费在线播放| av专区在线播放| 中文欧美无线码| 欧美激情在线99| av在线老鸭窝| 国产真实伦视频高清在线观看| 男女视频在线观看网站免费| 免费观看精品视频网站| eeuss影院久久| 国产成人a∨麻豆精品| 午夜精品一区二区三区免费看| 91午夜精品亚洲一区二区三区| 免费大片18禁| 国产日韩欧美在线精品| 91av网一区二区| 最后的刺客免费高清国语| 国产探花在线观看一区二区| 亚洲精品色激情综合| 在线免费十八禁| 尾随美女入室| 国产成年人精品一区二区| 少妇人妻精品综合一区二区| 三级男女做爰猛烈吃奶摸视频| 国产av不卡久久| 日本黄色片子视频| 永久网站在线| av在线播放精品| 国产一区有黄有色的免费视频 | 成人国产麻豆网| 国产老妇伦熟女老妇高清| 日日摸夜夜添夜夜添av毛片| 亚洲一级一片aⅴ在线观看| 精品人妻熟女av久视频| 国产老妇女一区| 成人特级av手机在线观看| 中文字幕免费在线视频6| av.在线天堂| 日韩强制内射视频| 国产久久久一区二区三区| 大话2 男鬼变身卡| 91精品伊人久久大香线蕉| 欧美97在线视频| 亚洲精品自拍成人| 亚洲精品国产av成人精品| av专区在线播放| 看非洲黑人一级黄片| 国产亚洲精品久久久com| 国产伦一二天堂av在线观看| 亚洲av免费在线观看| 国产乱来视频区| 老司机影院毛片| 91久久精品国产一区二区三区| 亚洲美女搞黄在线观看| 久久99蜜桃精品久久| 一边摸一边抽搐一进一小说| 久久久久久久久久成人| 色哟哟·www| 老司机福利观看| 我的老师免费观看完整版| 精品免费久久久久久久清纯| 秋霞伦理黄片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女xxoo啪啪120秒动态图| 日韩欧美国产在线观看| 人人妻人人澡欧美一区二区| 99久久精品国产国产毛片| 亚洲性久久影院| 99热网站在线观看| 如何舔出高潮| 日日啪夜夜撸| 男人和女人高潮做爰伦理| 国产探花极品一区二区| 综合色av麻豆| 欧美一级a爱片免费观看看| 亚洲精品成人久久久久久| 欧美色视频一区免费| 成人高潮视频无遮挡免费网站| 黑人高潮一二区| 欧美3d第一页| 精品人妻偷拍中文字幕| 亚洲va在线va天堂va国产| 国产亚洲精品av在线| 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| 中文天堂在线官网| 美女cb高潮喷水在线观看| 久久精品国产亚洲av涩爱| 亚洲国产欧洲综合997久久,| 久久精品国产鲁丝片午夜精品| 美女脱内裤让男人舔精品视频| 欧美精品国产亚洲| 国产精品久久电影中文字幕| 精品一区二区三区人妻视频| 天美传媒精品一区二区| 国产在线男女| 日韩三级伦理在线观看| 亚洲欧美精品综合久久99| 一区二区三区高清视频在线| 日韩欧美精品免费久久| 人人妻人人看人人澡| 国产真实伦视频高清在线观看| 舔av片在线| 日韩大片免费观看网站 | 国产亚洲91精品色在线| 国产中年淑女户外野战色| 看十八女毛片水多多多| 日韩大片免费观看网站 | 国产亚洲精品av在线| 色噜噜av男人的天堂激情| 毛片女人毛片| 在线播放无遮挡| 亚洲欧洲国产日韩| 国产在视频线精品| 亚洲在线自拍视频| 久久99精品国语久久久| 亚洲五月天丁香| 美女被艹到高潮喷水动态| 欧美日韩在线观看h| 视频中文字幕在线观看| 久久草成人影院| 亚洲经典国产精华液单| 色综合亚洲欧美另类图片| 午夜免费激情av| 黄色欧美视频在线观看| 亚洲精品自拍成人| 免费大片18禁| 精品久久久久久成人av| 高清午夜精品一区二区三区| 99在线人妻在线中文字幕| a级毛色黄片| 看十八女毛片水多多多| 国产成人精品婷婷| 韩国av在线不卡| www.av在线官网国产| 美女黄网站色视频| 久久精品夜色国产| 亚洲在久久综合| 国产私拍福利视频在线观看| 久久久国产成人免费| 国产精品不卡视频一区二区| 日本色播在线视频| 一边摸一边抽搐一进一小说| 美女被艹到高潮喷水动态| 国产一区二区三区av在线| 日韩欧美精品免费久久| 国产伦理片在线播放av一区| 草草在线视频免费看| 男插女下体视频免费在线播放| 特级一级黄色大片| 观看美女的网站| 秋霞伦理黄片| 看非洲黑人一级黄片| 国产亚洲av嫩草精品影院| 不卡视频在线观看欧美| 免费电影在线观看免费观看| 国产视频内射| 国产av在哪里看| 国产精品美女特级片免费视频播放器| 亚洲电影在线观看av| 18禁在线无遮挡免费观看视频| 日本熟妇午夜| 69人妻影院| kizo精华| 欧美性猛交黑人性爽| 成人美女网站在线观看视频| av免费在线看不卡| 老女人水多毛片| 亚洲怡红院男人天堂| 中文字幕制服av| 村上凉子中文字幕在线| 青春草国产在线视频| 精品人妻偷拍中文字幕| 国产淫语在线视频| 久久久久久久久久黄片| 国产精品麻豆人妻色哟哟久久 | 色播亚洲综合网| 国产探花极品一区二区| 深爱激情五月婷婷| 久久久久免费精品人妻一区二区| 色哟哟·www| 啦啦啦观看免费观看视频高清| 国产精品国产三级国产av玫瑰| 亚洲图色成人| 九九久久精品国产亚洲av麻豆| 非洲黑人性xxxx精品又粗又长| 中文天堂在线官网| 国产黄色视频一区二区在线观看 | 亚洲一区高清亚洲精品| 国产真实伦视频高清在线观看| 人妻系列 视频| 22中文网久久字幕| 亚洲四区av| 国产精品久久久久久久电影| 久久国内精品自在自线图片| av在线蜜桃| 大香蕉97超碰在线| 大话2 男鬼变身卡| 午夜爱爱视频在线播放| 国内少妇人妻偷人精品xxx网站| 日韩亚洲欧美综合| 男的添女的下面高潮视频| 久久久色成人| 欧美日韩一区二区视频在线观看视频在线 | 夜夜看夜夜爽夜夜摸| 午夜福利在线观看吧| 少妇高潮的动态图| 免费搜索国产男女视频| 免费av不卡在线播放| 卡戴珊不雅视频在线播放| 黄色配什么色好看| 一个人观看的视频www高清免费观看| 免费看光身美女| 亚洲av一区综合| 国产精品人妻久久久久久| 成人特级av手机在线观看| a级一级毛片免费在线观看| 国产一区二区在线av高清观看| 免费人成在线观看视频色| 天堂√8在线中文| 联通29元200g的流量卡| 精品久久国产蜜桃| 少妇的逼水好多| 淫秽高清视频在线观看| av在线亚洲专区| 免费看美女性在线毛片视频| 亚洲国产日韩欧美精品在线观看| 亚洲av一区综合| 内射极品少妇av片p| 观看美女的网站| 国产探花极品一区二区| 久久韩国三级中文字幕| .国产精品久久| videossex国产| 99视频精品全部免费 在线| 一区二区三区免费毛片| 国产精品久久久久久久电影| 午夜精品在线福利| 久久久亚洲精品成人影院| 国产精品美女特级片免费视频播放器| 汤姆久久久久久久影院中文字幕 | 亚洲精品乱码久久久v下载方式| 精品国产三级普通话版| 日韩欧美国产在线观看| 亚洲熟妇中文字幕五十中出| 赤兔流量卡办理| 舔av片在线| 高清视频免费观看一区二区 | 欧美极品一区二区三区四区| 亚洲伊人久久精品综合 | 男女下面进入的视频免费午夜| 午夜福利成人在线免费观看| 久久久成人免费电影| 亚洲国产精品国产精品| 亚洲国产欧美在线一区| 久久久久久久久中文| 国产视频首页在线观看| 内射极品少妇av片p| 亚洲av熟女| 三级经典国产精品| 午夜福利视频1000在线观看| 久久久久久九九精品二区国产| 午夜激情福利司机影院| 美女大奶头视频| 亚洲人成网站在线播| 久久久成人免费电影| 久久人妻av系列| 中文精品一卡2卡3卡4更新| 可以在线观看毛片的网站| 国产亚洲最大av| 欧美日韩国产亚洲二区| 久久韩国三级中文字幕| 免费观看在线日韩| 成人av在线播放网站| 欧美又色又爽又黄视频| 女的被弄到高潮叫床怎么办| 国产精品一及| 99久久精品一区二区三区| 国产高清国产精品国产三级 | 2021少妇久久久久久久久久久| 日本-黄色视频高清免费观看| 麻豆成人午夜福利视频| av在线蜜桃| 在线免费观看不下载黄p国产| 2022亚洲国产成人精品| 久热久热在线精品观看| 久久久久久久久大av| 夜夜爽夜夜爽视频| 中国国产av一级| 内射极品少妇av片p| 91精品伊人久久大香线蕉| 国产熟女欧美一区二区| av女优亚洲男人天堂| 午夜福利在线在线| 精品久久久久久成人av| 亚洲18禁久久av| 国产激情偷乱视频一区二区| av在线天堂中文字幕| 久久鲁丝午夜福利片| av卡一久久| 亚洲中文字幕日韩| 少妇人妻精品综合一区二区| 3wmmmm亚洲av在线观看| 十八禁国产超污无遮挡网站| 亚洲欧洲日产国产| 欧美人与善性xxx| 毛片一级片免费看久久久久| 精品一区二区三区人妻视频| 韩国高清视频一区二区三区| 国产伦精品一区二区三区视频9| 亚洲人成网站高清观看| 噜噜噜噜噜久久久久久91| 免费看光身美女| av国产久精品久网站免费入址| 国产精品一区二区性色av| 韩国av在线不卡| 看非洲黑人一级黄片| 老女人水多毛片| 成年女人永久免费观看视频| 特级一级黄色大片| 国语对白做爰xxxⅹ性视频网站| 1000部很黄的大片| 成人二区视频| 亚洲国产精品专区欧美| 亚洲在线自拍视频| 亚洲av电影在线观看一区二区三区 | 干丝袜人妻中文字幕| 亚洲欧美成人综合另类久久久 | 青青草视频在线视频观看| 久久久久久久久中文| 成人午夜高清在线视频| 日本-黄色视频高清免费观看| 狂野欧美白嫩少妇大欣赏| 免费观看a级毛片全部| videos熟女内射| 久久精品综合一区二区三区| 伊人久久精品亚洲午夜| 亚洲欧美清纯卡通| 亚洲精品久久久久久婷婷小说 | 国产亚洲5aaaaa淫片| 欧美成人免费av一区二区三区| 亚洲一级一片aⅴ在线观看| 赤兔流量卡办理| 亚洲欧美日韩卡通动漫| av在线蜜桃| 干丝袜人妻中文字幕| 久久久久久久久大av| 精品久久久久久久人妻蜜臀av| 一区二区三区四区激情视频| 国内精品美女久久久久久| 麻豆av噜噜一区二区三区| 国内揄拍国产精品人妻在线| 可以在线观看毛片的网站| 嫩草影院新地址| 亚洲第一区二区三区不卡| 国产黄色视频一区二区在线观看 | 搞女人的毛片| 18禁裸乳无遮挡免费网站照片| 欧美一区二区精品小视频在线| 嘟嘟电影网在线观看| 97超视频在线观看视频| 在线播放无遮挡| 欧美xxxx黑人xx丫x性爽| 亚洲精品aⅴ在线观看| 超碰97精品在线观看| 久久久久久伊人网av| 人人妻人人看人人澡| 性插视频无遮挡在线免费观看| 麻豆久久精品国产亚洲av| 干丝袜人妻中文字幕| 丰满人妻一区二区三区视频av| 亚洲人成网站在线观看播放| 六月丁香七月| 久久欧美精品欧美久久欧美| 国产91av在线免费观看| 国产国拍精品亚洲av在线观看| 一个人看视频在线观看www免费| 午夜福利视频1000在线观看| 亚洲自拍偷在线| 男女那种视频在线观看| 99九九线精品视频在线观看视频| 欧美区成人在线视频| 99热全是精品| 免费av毛片视频| 人人妻人人看人人澡| 亚洲五月天丁香| 亚洲激情五月婷婷啪啪| 亚洲国产精品国产精品| 最近视频中文字幕2019在线8| 一级黄色大片毛片| 少妇人妻精品综合一区二区| 久久久成人免费电影| 日韩欧美在线乱码| 最后的刺客免费高清国语| 日产精品乱码卡一卡2卡三| 级片在线观看| 偷拍熟女少妇极品色| 欧美区成人在线视频| 建设人人有责人人尽责人人享有的 | 午夜老司机福利剧场| 日产精品乱码卡一卡2卡三| videos熟女内射| 国产精品女同一区二区软件| 永久网站在线| 亚洲国产成人一精品久久久| 床上黄色一级片| 伦理电影大哥的女人| 成人av在线播放网站| 亚洲欧洲国产日韩| 国产精品人妻久久久久久| 一区二区三区乱码不卡18| 中文欧美无线码| 精品久久久噜噜| 嫩草影院新地址| 边亲边吃奶的免费视频| 欧美极品一区二区三区四区| 久久久久久久午夜电影| 国产精品1区2区在线观看.| 韩国av在线不卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 韩国av在线不卡| 免费电影在线观看免费观看| 搞女人的毛片| 简卡轻食公司| www.色视频.com| 麻豆成人av视频| 久久久久网色| 草草在线视频免费看| 亚洲欧美成人精品一区二区| 免费观看的影片在线观看| 97人妻精品一区二区三区麻豆| 成人一区二区视频在线观看| 毛片一级片免费看久久久久| 蜜臀久久99精品久久宅男| av在线老鸭窝| 岛国在线免费视频观看| 亚洲精品自拍成人| 国产色爽女视频免费观看| 免费观看性生交大片5| 边亲边吃奶的免费视频| 亚洲精品一区蜜桃| 少妇裸体淫交视频免费看高清| 1000部很黄的大片| 直男gayav资源| 超碰av人人做人人爽久久| 我的老师免费观看完整版| 欧美成人午夜免费资源| 国产视频首页在线观看| 亚洲电影在线观看av| 午夜精品国产一区二区电影 | 亚洲自拍偷在线| 欧美潮喷喷水| 能在线免费看毛片的网站| 亚洲精品日韩av片在线观看| 国产欧美日韩精品一区二区| 国模一区二区三区四区视频| 亚洲精品日韩av片在线观看| 精品久久久久久久人妻蜜臀av| 亚洲精品国产av成人精品| 日韩中字成人| 免费av观看视频| 禁无遮挡网站| 最新中文字幕久久久久| 亚洲精华国产精华液的使用体验| 亚洲四区av| 婷婷六月久久综合丁香| 日日摸夜夜添夜夜添av毛片| 精品一区二区三区视频在线| 搡女人真爽免费视频火全软件| 亚洲,欧美,日韩| 最近2019中文字幕mv第一页| 日韩三级伦理在线观看| h日本视频在线播放|