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

    利用地基紅外高光譜發(fā)射率數(shù)據(jù)進行云參數(shù)反演(1): 云相態(tài)判別

    2016-06-05 14:58:29孫學(xué)金高太長
    光譜學(xué)與光譜分析 2016年12期
    關(guān)鍵詞:水云相態(tài)發(fā)射率

    劉 磊、孫學(xué)金、高太長

    解放軍理工大學(xué)氣象海洋學(xué)院、江蘇 南京 211101

    利用地基紅外高光譜發(fā)射率數(shù)據(jù)進行云參數(shù)反演(1): 云相態(tài)判別

    劉 磊*、孫學(xué)金、高太長

    解放軍理工大學(xué)氣象海洋學(xué)院、江蘇 南京 211101

    云相態(tài)是氣候模式中的重要參數(shù)、也是遙感反演過程中進行云滴有效半徑、云水含量等微物理參數(shù)反演的重要前提。在研究了云層有效發(fā)射率光譜對云相態(tài)敏感性的基礎(chǔ)上、提出了基于云層有效發(fā)射率光譜的云相態(tài)表達特征、包括800~900 cm-1區(qū)域的有效發(fā)射率斜率、900~1 000 cm-1區(qū)域的有效發(fā)射率斜率、上述兩個區(qū)域的有效發(fā)射率斜率之差、862.1與989.8 cm-1的有效發(fā)射率之比、862.1與989.8 cm-1的有效發(fā)射率之差、1 900.1與2 029.3 cm-1的有效發(fā)射率之比、遠紅外窗區(qū)有效發(fā)射率平均值與900 cm-1有效發(fā)射率之比等7個特征。建立了利用支持向量機進行云相態(tài)判別的方法、開展了模擬數(shù)據(jù)驗證試驗、并利用遺傳算法優(yōu)化了支持向量機的徑向基核函數(shù)參數(shù)和懲罰因子。將該方法用于處理ARM計劃中SGP站點的AERI儀器獲得的數(shù)據(jù)、得到的云相態(tài)判別結(jié)果與Shupe提出的多儀器綜合判別結(jié)果進行了比較。結(jié)果表明、利用紅外波段不同窗區(qū)的有效發(fā)射率光譜特征可以實現(xiàn)發(fā)射率低于0.95的云層的相態(tài)判別、建立的基于支持向量機的云相態(tài)判別方法與Shupe方法的總體判別結(jié)果較為一致、但有約30%的云層由于發(fā)射率較大而標(biāo)記為不透明云?;诩t外高光譜發(fā)射率數(shù)據(jù)的云相態(tài)判別技術(shù)充分考慮了光譜斜率、比值和差值等信息、是較為穩(wěn)定有效的薄云相態(tài)判別方法。

    云相態(tài); 紅外高光譜輻射; 云發(fā)射率; 支持向量機

    引 言

    云相態(tài)主要分水相、冰相和混合相、不同相態(tài)的云具有不同的散射和吸收特性。云相態(tài)是氣候模式的重要參數(shù)、對輻射通量計算有很大影響。此外、云相態(tài)又是進行云滴有效半徑、云水含量等微物理參數(shù)反演的重要前提、相態(tài)判別錯誤直接引起單次散射特性計算錯誤、而這種錯誤會導(dǎo)致20%~100%的云滴有效半徑和光學(xué)厚度的誤差[1-2]。

    目前已有很多有關(guān)云相態(tài)方面的研究。ISCCP[3]利用云體溫度進行水云和冰云的判別、但實際上云相態(tài)并不完全由溫度所決定。一些學(xué)者通過利用短波紅外輻射、可見光輻射觀測來推斷云相態(tài)[2,4]。該方法只能用于白天、且對太陽天頂角和儀器觀測角度比較敏感。當(dāng)水云粒子尺寸太大或者冰云粒子尺寸太小時、云相態(tài)的識別會變得不確定。此外、由于云內(nèi)不同相態(tài)的水成物粒子的形狀、大小、密度和空間取向的不同、對一定偏振狀態(tài)的電磁波有著不同的退偏振效應(yīng)、可以通過分析這些觀測量的差異反演出的云相態(tài)[5]。如利用偏振激光雷達獲得的線性退偏比、通過經(jīng)驗閾值進行云相態(tài)判斷。也有學(xué)者[6]研究采用微波雷達的退偏振因子進行降水云的相態(tài)識別、但這對雷達系統(tǒng)的靈敏度有很高的要求。Shupe[1]提出采用聯(lián)合激光雷達退偏振比和毫米波雷達的反射率、多普勒速度、譜寬以及云體溫度進行云相態(tài)判別、是目前較為實用的地基云相態(tài)判別方法。

    利用紅外波段進行云相態(tài)識別有較長時間的研究歷史?;舅悸肥歉鶕?jù)云內(nèi)水相粒子與冰相粒子對在某些紅外波段的吸收率(復(fù)折射指數(shù)的虛部)特征的差異(這種差異表現(xiàn)在發(fā)射率上)來判定云的相態(tài)。在8~10 μm、水和冰的吸收系數(shù)幾乎相同、而在10~13 μm、冰的吸收要大于水、這會造成云體亮溫、發(fā)射率等有所不相同。Ackerman等[7]首次提出了基于8、11和12 μm的三光譜云相態(tài)識別方法。在此基礎(chǔ)上、Strabala等[8]進一步建立了8.5、11和12 μm的云相態(tài)識別算法。兩者均采用了亮溫差分析方法。Dan Lubin[9]和Turner等[10]針對高光譜紅外輻射數(shù)據(jù)分別提出利用亮溫斜率和云發(fā)射率比的方法進行云相態(tài)識別、均取得了較好的效果。這些方法對水汽吸收比較敏感、濕度廓線的測量誤差容易造成相態(tài)識別的誤差。需要說明的是、利用這些方法識別的云相態(tài)、確切地說應(yīng)該是“輻射意義上的云相態(tài)”[11]、與實際云的熱力學(xué)相態(tài)是有差異的。

    本文利用地基紅外高光譜輻射儀AERI[12-13]的觀測數(shù)據(jù)進行云相態(tài)判別。該儀器可以測量530~3 050 cm-1(19.0~3.3 μm)波段的大氣向下高光譜紅外輻射、光譜分辨率為1 cm-1。由于大氣向下紅外高光譜輻射數(shù)據(jù)不僅與云層物理屬性有關(guān)、還受到氣象條件的影響、本文在只與云層物理屬性有關(guān)的有效發(fā)射率基礎(chǔ)上開展工作、進行了云相態(tài)敏感性試驗并建立了基于發(fā)射率光譜的云相態(tài)特征表述方法、提出了利用支持向量機進行云相態(tài)判別的方法、最后將該方法與Shupe發(fā)展的基于地基多傳感器的云相態(tài)判別結(jié)果進行了比較分析。

    1 物理基礎(chǔ)

    1.1 用于云相態(tài)判別的窗口

    采用高光譜紅外輻射數(shù)據(jù)進行云參數(shù)反演的重要優(yōu)勢在于可以在吸收線之間的多個小窗區(qū)進行計算、從而能夠顯著減少云下大氣對輻射的影響[14]。用于模擬輻射傳輸?shù)能浖門urner發(fā)展的LBLDIS[15]、該軟件采用LBLRTM結(jié)合DISORT的方式、可以模擬云、氣溶膠的紅外高光譜輻射特性。結(jié)合前人所做工作[16]、確定了用于研究云相態(tài)判別的窗口、如表1所示。

    在紅外波段、水云和冰云的折射率虛部mi隨波長有一定的變化、具體如圖1(a)所示。進一步按照式(1)計算了吸收系數(shù)κ[17]、如圖1(b)所示。

    (1)

    表1 用于云相態(tài)判別研究的窗口

    圖1中同時繪制出了表1中小窗口的位置(用虛線表示)。可見、在遠紅外波段(如圖示的小于600 cm-1的三個窗口)、水的吸收大于冰的吸收; 而在熱紅外波段直到1 000 cm-1左右、冰的吸收大于水的吸收、在1 100 cm-1附近冰與水的吸收幾乎相同; 對于中紅外波段(如圖1所示的1 900~2 200 cm-1的區(qū)域)、冰的吸收也總體大于水的吸收、并有個別窗口兩者吸收系數(shù)接近。傳統(tǒng)的三波段判斷云相態(tài)的方法主要選擇的是熱紅外波段。從上述分析來看、理論上遠紅外波段和中紅外波段的窗口也蘊含了云相態(tài)的信息。Rathke[18]提出了加入20 μm(500 cm-1)波段信息進行云相態(tài)判別的方案、這充分利用了水和冰的吸收系數(shù)出現(xiàn)交叉的優(yōu)勢。但在遠紅外波段、PWV超過30 mm后將由于達到飽和而無法使用[10]。

    圖1 水云和冰云的折射率虛部(a)及對應(yīng)的吸收系數(shù)(b)

    1.2 冰云、水云的發(fā)射率光譜特征分析

    根據(jù)上述分析、進一步模擬了不同光學(xué)厚度和云滴有效半徑條件下、水云和冰云的有效發(fā)射率分布情況。在忽略云上大氣的輻射以及云底反射的地面和云下大氣發(fā)射輻射的前提下、云層有效發(fā)射率εeff可以通過式(2)獲得[19]

    (2)

    冰云的云滴形狀取六角柱狀、有效半徑一般在20~100 μm、典型冰云粒子特征尺度在30 μm左右、圖2—圖4中繪出了20、30和50 μm的情況。水云的云滴有效半徑一般取3~15 μm、典型水云粒子特征尺度在7 μm左右、圖2—圖4中繪出了3、7和10 μm的情況。兩種相態(tài)云層的光學(xué)厚度均取0.2、0.5、1.0、3.0和8.0。

    圖2 熱紅外波段的水云和冰云有效發(fā)射率特征

    由圖2可見、在熱紅外波段、冰云的有效發(fā)射率隨波長的變化不大、而水云則有較大的變化。在800~1 000 cm-1區(qū)域、水云的有效發(fā)射率斜率比冰云的要大很多、而且800~900 cm-1區(qū)域的有效發(fā)射率斜率比900~1 000 cm-1區(qū)域的更陡峭??梢岳眠@兩個區(qū)域有效發(fā)射率斜率的差異進行冰云和水云的區(qū)分。在1 100 cm-1左右的區(qū)域、水云和冰云的有效發(fā)射率幾乎都不隨波長發(fā)生變化。不同光學(xué)厚度的冰云和水云、其有效發(fā)射率也有一定差異、光學(xué)厚度在0.5~3.0之間時、冰云和水云之間的差異較大; 較小的光學(xué)厚度(如圖2所示的0.2)和較大的光學(xué)厚度(如圖2所示的8.0)時、冰云和水云之間的差異不大、但仍有一定的區(qū)分度、尤其在較小光學(xué)厚度時。

    圖3 遠紅外波段的水云和冰云有效發(fā)射率特征

    由圖3可見、在遠紅外波段的三個窗口、由于水云在此區(qū)域有較強的吸收作用、其有效發(fā)射率較冰云的大。結(jié)合800 cm-1左右熱紅外窗區(qū)冰云和水云的有效發(fā)射率可知、可以利用550 cm-1附近和800 cm-1附近窗口的有效發(fā)射率之比進行冰云和水云的區(qū)分。

    由圖4可見、在中紅外波段、冰云和水云的有效發(fā)射率隨波長的變化整體上不如熱紅外波段顯著。比較1 900 cm-1附近和2 050 cm-1附近幾個窗口的有效發(fā)射率、可以看出、利用這兩個區(qū)域云的有效發(fā)射率之比可以進行冰云和水云的區(qū)分。2 100 cm-1以上的四個窗口的有效發(fā)射率隨波長變化不大。

    圖4 中紅外波段的水云和冰云有效發(fā)射率特征

    分析冰云和水云在遠紅外、熱紅外和中紅外波段有關(guān)窗口的有效發(fā)射率情況可知、可以利用適當(dāng)?shù)挠行Оl(fā)射率斜率或者不同窗口的比值等組合進行云相態(tài)判別。

    2 基于有效發(fā)射率光譜的云相態(tài)特征參數(shù)

    據(jù)上述的敏感性試驗、進一步建立了7個用于表達云相態(tài)的特征。這些特征大多數(shù)采用有效發(fā)射率斜率或者有效發(fā)射率之比進行表達、這樣可以最大程度地減小由于其他氣象條件測量不精確引起的特征表達誤差。具體名稱及物理意義如下:

    (1)800~900 cm-1區(qū)域的有效發(fā)射率斜率:主要利用862.1、873.6、892.9和900.0 cm-1四個波段的數(shù)據(jù)進行擬合、簡稱Slope1。若斜率為負、則云粒子較小、可判為水云;若斜率接近0、則認為粒子較大、可判為冰云、如圖5(a)所示。

    圖5 基于有效發(fā)射率光譜的云相態(tài)特征參數(shù)

    (2)900~1 000 cm-1區(qū)域的有效發(fā)射率斜率:主要利用934.9、961.9、989.8和1 080.0 cm-1四個波段的數(shù)據(jù)進行擬合、簡稱Slope2。物理意義同(1)、如圖5(b)所示。

    (3)上述兩個區(qū)域有效發(fā)射率斜率之差:簡稱Slope_Diff、Slope_Diff= Slope2-Slope1。冰和水在800~900 cm-1區(qū)域和900~1 000 cm-1區(qū)域的吸收系數(shù)變化不同、水云的Slope1更陡峭、會出現(xiàn)差異。Slope_Diff接近0、可判為冰云; Slope_Diff為正、可判為水云、如圖5(c)所示。

    (4)862.1與989.8 cm-1的有效發(fā)射率之比:簡稱Ratio1。在熱紅外區(qū)域、這兩個波段的冰和水吸收系數(shù)差異很大。862.1 cm-1波段的冰的吸收系數(shù)顯著高于水的吸收系數(shù)、而在989.8 cm-1波段、兩者吸收系數(shù)幾乎相同。這兩個波段的有效發(fā)射率之比若在1附近、則可判為冰云、否則可以判為水云。如圖5(d)所示。

    (5)862.1與989.8 cm-1的有效發(fā)射率之差:簡稱Diff_E。物理意義同(4)。若兩波段有效發(fā)射率之差接近0、則判為冰云、否則可以判為水云。

    (6)1 900.1與2 029.3 cm-1的有效發(fā)射率之比:簡稱Ratio2。選擇這兩個波段的物理意義與(4)基本相同。若比值在1附近、可判為冰云、否則可判為水云。如圖5(e)所示。

    (7)遠紅外窗區(qū)有效發(fā)射率平均值與900 cm-1有效發(fā)射率之比:簡稱Ratio3。遠紅外窗區(qū)有效發(fā)射率平均值主要指利用530.4、561.2和572.3 cm-1三個波段的平均值。該區(qū)域是為數(shù)不多的水比冰的吸收系數(shù)高的區(qū)域。若比值在1附近、則可判為冰云、否則可以判為水云。具體如圖5(f)。

    此外、由圖5可見、在高發(fā)射率區(qū)域、尤其是有效發(fā)射率大于0.95以后、冰云和水云幾乎難以區(qū)分。有研究表明、在900 cm-1處水云和冰云的有效發(fā)射率幾乎一致、這一點也可由圖2看出。因此、本研究將900 cm-1處發(fā)射率大于0.95的云視為不透明云、并不參與云相態(tài)的識別。

    3 基于支持向量機的云相態(tài)判別方法

    3.1 訓(xùn)練樣本集的建立

    利用LBLDIS模擬建立了冰云、水云和混合相云的訓(xùn)練樣本集。其中混合相云采用Sun等[20]的散射特性計算方案、冰云占比可定義為f=τi/(τw+τi)。模擬計算時、光學(xué)厚度均取0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、1.0、1.5、2.0、2.5、3.0、4.0、6.0和8.0。水云(取f=0)的云滴有效半徑取2、3、5、7、10和15 μm、冰云(取f=1)的云滴有效半徑取10、20、30、50、70和100 μm、混合云(取f=0.5)的云滴有效半徑組合取(reff,w=3 μm、reff,i=20 μm)、(reff,w=3 μm、reff,i=30 μm)、(reff,w=5 μm、reff,i=30 μm)、(reff,w=7 μm、reff,i=20 μm)、(reff,w=7 μm、reff,i=30 μm)、(reff,w=10 μm、reff,i=40 μm)、分別得到三種相態(tài)各96組訓(xùn)練樣本。

    需要注意的是、f可以按其他比例獲得混合相云特征、而且云滴有效半徑的混合方案也可以取更多種。但從訓(xùn)練數(shù)量均衡的角度出發(fā)、三種云相態(tài)初始訓(xùn)練樣本集的數(shù)目設(shè)為相同、其他數(shù)據(jù)作為測試數(shù)據(jù)。

    3.2 模擬驗證結(jié)果

    模擬驗證時、采用的支持向量機核函數(shù)為徑向基核函數(shù)。利用遺傳算法對徑向基核函數(shù)中的σ和懲罰因子C這兩個參數(shù)進行了優(yōu)選、最終確定的參數(shù)為:C取84.448 5、σ取3.031 4。利用3.1中方法模擬f為0.1~0.9(以0.1為步長)共9組混合相云的特征。模擬數(shù)據(jù)的驗證結(jié)果如表2所示。

    表2 混合相云的模擬分類結(jié)果

    可見、在f<0.3(水云占比超過70%)或者f>0.7(冰云占比超過70%)時、混合云的判別更偏向于占比較多的云類。這些情況下準(zhǔn)確判斷為混合相云的比例將低于60%。當(dāng)然、可以將所有不同冰云水云比例的模擬數(shù)據(jù)集加入訓(xùn)練集進行SVM訓(xùn)練、從而獲得更高的模擬驗證分類準(zhǔn)確性。但這一方面會使得訓(xùn)練樣本數(shù)失衡(混合相云的組合數(shù)量遠遠多于純水或純冰相); 另一方面考慮到這些某一相態(tài)占比較大的混合云的特征與該純相態(tài)的特征非常接近、加入這些訓(xùn)練樣本的結(jié)果會得到很大的懲罰因子、以期加大對錯分樣本的懲罰程度、從而使得SVM的泛化推廣性能大大降低。綜合考慮這些因素、加之各種比例的混合云模擬驗證效果較為合理、可以認為第三節(jié)中提出的特征以及采用優(yōu)化參數(shù)后的SVM進行云相態(tài)判別是可行的。

    4 實際判別結(jié)果及分析

    將以上提出的云相態(tài)判別算法得到的結(jié)果與Shupe提出的云相態(tài)地基綜合識別算法[1]的結(jié)果進行比較。Shupe云相態(tài)判別使用的數(shù)據(jù)包括激光雷達的退偏振比和后向散射強度、毫米波雷達的反射率、多普勒速度和譜寬以及溫度廓線等。其基本判別過程包括[1]:利用激光雷達數(shù)據(jù)進行云相態(tài)初判、利用毫米波測云雷達數(shù)據(jù)進行云相態(tài)修正、標(biāo)記降水?dāng)?shù)據(jù)、標(biāo)記激光雷達未能觀測到的數(shù)據(jù)、利用溫度廓線進行再分類、利用液態(tài)水路徑修正分類結(jié)果、均質(zhì)化處理等。上述過程涉及到的閾值大多是經(jīng)驗結(jié)果。

    采用的分析數(shù)據(jù)來自美國能源部大氣輻射測量計劃(ARM計劃)的SGP站點2010年3月份數(shù)據(jù)。首先、基于激光雷達的云底高數(shù)據(jù)將AERI觀測時次有云的數(shù)據(jù)選出、共3 374組; 然后、分別利用Shupe方法和本方法進行云相態(tài)判別。由于Shupe方法獲得的是云相態(tài)廓線、為了進行比較需要進一步進行統(tǒng)計處理。主要方法是:若云相態(tài)廓線存在冰云、水云及混合相云這三種相態(tài)的兩個及以上的相態(tài)且各相態(tài)的云點占該廓線總云點數(shù)均大于30%、則該時次的云層判為混合相云; 若存在雨、雪和毛毛雨、則判為有降水; 除此之外、將占比最大的相態(tài)作為該時次的云相態(tài)、若上述凝結(jié)物的類別均不存在、則判該時次為無云。最終Shupe方法統(tǒng)計后的云層可以為:水云、冰云、混合云、無云和降水等5類、本方法可能將云層判別為:水云、冰云、混合云和不透明云等4類。

    總體情況如表3所示??梢妰煞N方法對冰云和混合云的月總體判別結(jié)果較為一致、本方法得到的水云數(shù)量與Shupe方法差異較大。分析發(fā)現(xiàn)、出現(xiàn)差異的這些水云大多屬于發(fā)射率較大的云層、本方法將發(fā)射率超過0.95的云層判為不透明云、共1 115組、約占30%。這其中包含了183組降水情況。

    表3 2010年3月份SGP站點云相態(tài)分類結(jié)果比較表

    圖6 SGP站點2010年3月2日的激光雷達和毫米波云雷達數(shù)據(jù)

    選擇以SGP站點2010年3月2日的觀測個例做進一步分析。圖6是激光雷達和毫米波雷達觀測數(shù)據(jù)。其中(a)為激光雷達的后向散射強度、(b)為激光雷達的線性退偏振比、(c)為毫米波雷達的反射率、(d)為毫米波雷達的多普勒速度、(e)為毫米波雷達的譜寬。

    利用Shupe方法得到的云相態(tài)廓線如圖7(a)所示。00 h—03 h天空存在兩層云、其中上層為冰云、下層為水云; 03 h—08 h天空存在一層冰云; 11 h—22 h天空存在一層水云、其中11 h—15 h激光雷達觀測到云層而毫米波雷達沒能觀測到。已有研究表明[21]毫米波雷達對于光學(xué)厚度小于2 的薄卷云的漏測率超過了30%、且難以探測到光學(xué)厚度小于0.2 的薄卷云。

    利用AERI觀測的數(shù)據(jù)、基于本研究的支持向量機云相態(tài)判別方法獲得的結(jié)果如圖7(b)所示。該方法獲得的結(jié)果表達了整層云在輻射意義下的云相態(tài)、為了與Shupe方法獲得的云相態(tài)廓線進行比較、將判別結(jié)果繪制在該云層所在云底高度上。其中云底高度采用的是激光雷達反演的結(jié)果??梢?、兩種方法得到的云相態(tài)結(jié)果有很好的一致性。天空出現(xiàn)兩層云時、AERI的紅外輻射光譜只反映了下層水云特征。由于毫米波雷達、激光雷達與AERI儀器有一定的距離、且AERI觀測一次需要3 min、若天空云層不均勻且下層水云的縫隙足以露出上層冰云時、可能會出現(xiàn)以下情況:(1)AERI的原始觀測輻射光譜包含了水云和冰云的信息、此時Shupe的結(jié)果可能是上層冰云、下層水云、而利用AERI數(shù)據(jù)判別的結(jié)果為混合相云、如03 h數(shù)據(jù); (2)激光雷達提供的云底高度是低層云的云底高度、由于AERI和激光雷達有一定的距離、AERI數(shù)據(jù)可能測量的是上層冰云的輻射光譜、如01、04和06 h數(shù)據(jù)。這一情況也說明本研究采用發(fā)射率差值、斜率等特征進行云相態(tài)區(qū)分時、對云底高度的準(zhǔn)確性要求不高、判別結(jié)果更多依賴于發(fā)射率光譜的相對大小。分析圖7還發(fā)現(xiàn)、在01 h—02 h、18 h—20 h左右、存在AERI原始輻射光譜飽和導(dǎo)致無法判別云相態(tài)的數(shù)據(jù)、該段云層的發(fā)射率大多數(shù)大于0.95、標(biāo)記為不透明云。

    圖7 SGP站點2010年3月2日的云相態(tài)判別結(jié)果

    進一步、給出了部分時次的TSI全天空云圖、如圖8所示。由于TSI只能獲得白天的觀測數(shù)據(jù)、所以選擇了14:00、15:00、17:30和20:00的TSI圖像??梢园l(fā)現(xiàn)、14 h之后天空以層積云為主。18 h—22 h的云層較厚、部分觀測時次AERI的輻射光譜飽和、故難以進行云相態(tài)判別。

    圖8 SGP站點2010年3月2日部分時次的TSI圖像

    分析表3還發(fā)現(xiàn)、Shupe方法將492組判為無云。這些數(shù)據(jù)大多數(shù)是激光雷達探測到有云而毫米波云雷達沒有探測到云層、且該云層的線性退偏比和后向散射強度均較小。以2010年3月22日12 h—18 h的數(shù)據(jù)為例(如圖9—圖11所示)、查看TSI圖像可見、這段時間應(yīng)存在較高且很薄的卷云。這一誤判的原因是Shupe方法中將線性退偏比小于0.25的云層判為了氣溶膠層。

    圖9 SGP站點2010年3月22日的激光雷達和毫米波云雷達數(shù)據(jù)

    圖10 SGP站點2010年3月22日的云相態(tài)判別結(jié)果

    圖11 SGP站點2010年3月22日部分時次的TSI圖像

    5 結(jié) 論

    研究了云層有效發(fā)射率光譜對云相態(tài)的敏感性、提出了基于云層有效發(fā)射率光譜的云相態(tài)表達特征、建立了利用支持向量機進行云相態(tài)判別的方法、并將所提出的方法獲得的云相態(tài)判別結(jié)果與Shupe提出的多儀器綜合判別結(jié)果進行了比較。研究表明、對于發(fā)射率低于0.95的云層、利用紅外波段不同窗區(qū)的有效發(fā)射率斜率或比值可以實現(xiàn)云相態(tài)的判別。其中800~900 cm-1區(qū)域和900~1 000 cm-1區(qū)域的有效發(fā)射率斜率以及兩個區(qū)域有效發(fā)射率斜率之差、862.1與989.8 cm-1的有效發(fā)射率比值和差值、1 900.1與2 029.3 cm-1的有效發(fā)射率之比、遠紅外窗區(qū)有效發(fā)射率平均值與900 cm-1有效發(fā)射率之比等特征有效地反映了云的相態(tài)。建立的基于支持向量機的云相態(tài)判別方法與Shupe建立的地基綜合云相態(tài)廓線判別方法總體判別結(jié)果較為一致。約30%的云層發(fā)射率較大、本方法標(biāo)記為不透明云。Shupe方法由于閾值設(shè)置的原因、對約10%的線性退偏比和后向散射強度均較小的薄云層誤判為晴空。本研究基于紅外高光譜發(fā)射率數(shù)據(jù)的云相態(tài)判別技術(shù)充分考慮了光譜斜率、比值等信息、是較為穩(wěn)定有效的薄云相態(tài)判別方法。

    [1] Shupe M D. Geophysical Research Letters,2007,34:L22809.

    [2] Key J R,Intrieri J M. J. Appl. Meteor.,2000,39: 1797.

    [3] Rossow W B,Schiffer R A. Bull. Amer. Meteor. Soc.,1999,80: 2261.

    [4] Daniel J S,Solomon S,Portmann R W,et al. Journal of Geophysical Research,2002,107: 4599.

    [5] Sassen K. Bulletin of the American Meteorological Society,1991,72(12): 1848.

    [7] Ackerman S A,Smith W L,Revercomb H E,et al. Mon. Wea. Rev.,1990,118: 2377.

    [8] Strabala K I,Ackerman S A,Menzel W P. Journal of Applied Meteorology 1994,33: 212.

    [9] Lubin D. Journal of Geophysical Research,2004,109: D04204.

    [10] Turner D D,Ackerman S A,Baum B A,et al. Journal of Applied Meteorology,2003,42(6): 701.

    [11] Garrett T J,Zhao C. Atmospheric Measurement Techniques,2013,6(5): 1227.

    [12] Knuteson R O,Revercomb H E,Best F A,et al. J. Atmos. Oceanic Technol.,2004,21: 1763.

    [13] Knuteson R O,Revercomb H E,Best F A,et al. J. Atmos. Oceanic Technol.,2004,21: 1777.

    [14] Smith W L,Ma X L,Ackerman S A,et al. Journal of Atmospheric Sciences,1993,50: 1708.

    [15] Turner D D. Journal of Applied Meteorology,2005,44(4): 427.

    [16] DeSlover D H,Smith W L,Piironen P K,et al. J. Atmos. Oceanic Technol.,1999,16: 251.

    [17] Bohren C F,Huffman D R. Absorption and Scattering of Light by Small Particles. Germany: Wiley-VCH Verlag GmbH & Co. KGaA 2007.

    [18] Rathke C. Geophysical Research Letters,2002,29: 1209.

    [19] Linjun P,Daren L. Atmpspheric and Oceanic Science Letters,2012,5(5): 439.

    [20] Sun Z,Shine K P. J. Climate,1995、(8): 1874.

    [21] Borg L A,Holz R E,Turner D D. Geophysical Research Letters,2011,38(5): L05807.

    *Corresponding author

    Research on Cloud Phase Detemination Using Infrared Emissivity Spectrum Data (1): Cloud Phase Determination

    LIU Lei*,SUN Xue-jin,GAO Tai-chang

    College of Meteorology and Oceanography,PLA University of Science and Technology,Nanjing 211101,China

    As a key factor in the climate model,cloud phase is an important prerequisite to performing cloud property retrievals from remote sensor measurements. The ability to infer cloud phase using cloud emissivity spectra is investigated by numerical simulations. It is shown that for emissivity below 0.95,several spectral features such as the slopes,the ratios and the differences of the emissivity are consistent with the variation of cloud phase in some spectral regions. Specifically,these features include the slope of the cloud emissivity between 800 and 900 cm-1,the slope of the cloud emissivity between 900 and 1 000 cm-1,the difference in the mean emissivity between above-mentioned two regions,the ratio of the emissivity at 862.1 cm-1to the emissivity at 989.8 cm-1,the difference in the emissivity between 862.1 and 989.8 cm-1,the ratio of the emissivity at 1 900.1 cm-1to the emissivity at 2 029.3 cm-1,the ratio of the mean emissivity for far-infrared region to the emissivity at 900 cm-1. A cloud phase classifier is proposed based on support vector machines (SVM). A series of simulations including various cloud patterns are performed. The RBF kernel function parameters and the penalty factor of SVM are selected by using the genetic algorithm. The phase determination algorithm is applied for collecting data from the AERI at the SGP site. The results from the ground-based multisensor cloud phase classifier proposed by Shupe are used to validate the phase determination algorithm. It is found the two results are consistent in general. 30% clouds are indicated as opaque due to its high emissivity. The cloud with small lidar’s depolarization is misclassified as clear sky by the Shupe method. It can be concluded that the proposed algorithm considering the spectral information (spectral slopes,ratios and differences) is efficient for cloud phase determination of thin cloud.

    Cloud phase; Hyperspectral infrared radiance; Cloud emissivity; Support vector machine

    Dec. 23,2015; accepted Apr. 15,2016)

    2015-12-23、

    2016-04-15

    國家自然科學(xué)基金項目(41205125、41575024)資助

    劉 磊、1983年生、解放軍理工大學(xué)氣象海洋學(xué)院講師 e-mail: liuleidll@gmail.com *通訊聯(lián)系人

    P414.9

    A

    10.3964/j.issn.1000-0593(2016)12-3885-10

    猜你喜歡
    水云相態(tài)發(fā)射率
    微波水云間
    保健與生活(2023年6期)2023-03-17 08:39:54
    水云間
    氧氣A(O,O)波段氣輝體發(fā)射率和臨邊輻射強度模擬與分析
    古琴曲《瀟湘水云》的題解流變考
    胡忌先生讀《水云村稿》札記一則
    中華戲曲(2017年2期)2017-02-16 06:53:16
    SBS改性瀝青相態(tài)結(jié)構(gòu)的參數(shù)化表征方法
    低溫狀態(tài)下的材料法向發(fā)射率測量
    四川省降水相態(tài)識別判據(jù)研究
    PS/PLA共混物的相態(tài)結(jié)構(gòu)及其發(fā)泡行為研究
    中國塑料(2015年4期)2015-10-14 01:09:21
    塔克拉瑪干沙漠地表發(fā)射率及分布變化特征
    国产有黄有色有爽视频| 99精品在免费线老司机午夜| 亚洲精品成人av观看孕妇| 亚洲国产毛片av蜜桃av| 丝袜美腿诱惑在线| 国产精品秋霞免费鲁丝片| 国产高清国产精品国产三级| 手机成人av网站| 麻豆乱淫一区二区| 国产在线观看jvid| 满18在线观看网站| av又黄又爽大尺度在线免费看| 日日夜夜操网爽| 国产一区二区三区视频了| 精品一品国产午夜福利视频| 国产精品 国内视频| 日韩中文字幕视频在线看片| www日本在线高清视频| 精品亚洲成a人片在线观看| 大香蕉久久网| 久久久国产成人免费| 搡老乐熟女国产| 下体分泌物呈黄色| 国产精品香港三级国产av潘金莲| 久久亚洲真实| 国产精品欧美亚洲77777| 两性午夜刺激爽爽歪歪视频在线观看 | 色94色欧美一区二区| 国产av又大| 多毛熟女@视频| 五月天丁香电影| 搡老岳熟女国产| 俄罗斯特黄特色一大片| 成年女人毛片免费观看观看9 | 精品福利观看| 在线永久观看黄色视频| 十八禁人妻一区二区| www.熟女人妻精品国产| 精品国产亚洲在线| 国产欧美日韩精品亚洲av| 国产精品偷伦视频观看了| 黑人猛操日本美女一级片| www日本在线高清视频| 在线观看一区二区三区激情| 纵有疾风起免费观看全集完整版| 久久久精品94久久精品| 热re99久久精品国产66热6| 久久天堂一区二区三区四区| 亚洲av国产av综合av卡| 啦啦啦中文免费视频观看日本| 精品一区二区三区视频在线观看免费 | 精品国内亚洲2022精品成人 | 亚洲美女黄片视频| 国产男女超爽视频在线观看| 成人国语在线视频| 夜夜骑夜夜射夜夜干| 十八禁高潮呻吟视频| 女性生殖器流出的白浆| 日本vs欧美在线观看视频| 午夜91福利影院| 午夜免费成人在线视频| 久久久久久久精品吃奶| 亚洲精品国产一区二区精华液| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产av新网站| 在线观看免费高清a一片| 久久久久久久国产电影| 久久久国产精品麻豆| 国产成人一区二区三区免费视频网站| 精品国产国语对白av| 国产精品熟女久久久久浪| 男人舔女人的私密视频| 天堂动漫精品| 纵有疾风起免费观看全集完整版| 久久久久久久久免费视频了| 男人舔女人的私密视频| 狠狠狠狠99中文字幕| 亚洲第一欧美日韩一区二区三区 | 十八禁网站免费在线| 欧美亚洲 丝袜 人妻 在线| 精品第一国产精品| 国产无遮挡羞羞视频在线观看| 国产又色又爽无遮挡免费看| 少妇裸体淫交视频免费看高清 | 国产精品香港三级国产av潘金莲| 黄片大片在线免费观看| 少妇猛男粗大的猛烈进出视频| 老熟女久久久| cao死你这个sao货| 丰满迷人的少妇在线观看| 麻豆乱淫一区二区| 视频区欧美日本亚洲| 女警被强在线播放| 老司机深夜福利视频在线观看| 国产成人系列免费观看| 视频在线观看一区二区三区| 一级毛片电影观看| 亚洲精品中文字幕在线视频| 大陆偷拍与自拍| 亚洲免费av在线视频| 亚洲va日本ⅴa欧美va伊人久久| 丁香六月天网| 色精品久久人妻99蜜桃| 91麻豆av在线| 亚洲va日本ⅴa欧美va伊人久久| 一区福利在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲色图综合在线观看| 999久久久国产精品视频| 久久精品人人爽人人爽视色| 欧美一级毛片孕妇| 在线观看66精品国产| 国产精品影院久久| 日韩人妻精品一区2区三区| 正在播放国产对白刺激| 亚洲精品国产区一区二| 国产有黄有色有爽视频| 成人18禁高潮啪啪吃奶动态图| 美女视频免费永久观看网站| 中文字幕人妻丝袜制服| 法律面前人人平等表现在哪些方面| 女性被躁到高潮视频| av欧美777| 国产黄色免费在线视频| 少妇猛男粗大的猛烈进出视频| videos熟女内射| 狠狠狠狠99中文字幕| av一本久久久久| 妹子高潮喷水视频| 国产三级黄色录像| 国产成人免费观看mmmm| 老司机亚洲免费影院| 精品亚洲成国产av| 男女边摸边吃奶| 成人影院久久| 精品亚洲乱码少妇综合久久| 久久亚洲真实| 精品高清国产在线一区| 成年动漫av网址| 欧美亚洲 丝袜 人妻 在线| 窝窝影院91人妻| 亚洲成人免费av在线播放| 51午夜福利影视在线观看| 91麻豆av在线| 国产高清视频在线播放一区| 成人亚洲精品一区在线观看| 久久久久精品国产欧美久久久| 日韩中文字幕欧美一区二区| 成在线人永久免费视频| 午夜激情久久久久久久| 欧美日韩亚洲高清精品| 国产精品自产拍在线观看55亚洲 | 欧美 亚洲 国产 日韩一| h视频一区二区三区| 日韩人妻精品一区2区三区| 午夜免费鲁丝| 人人妻,人人澡人人爽秒播| 满18在线观看网站| 黄色毛片三级朝国网站| 亚洲国产毛片av蜜桃av| 久9热在线精品视频| 亚洲av成人一区二区三| 国产精品国产av在线观看| 国产av又大| 黑人操中国人逼视频| 波多野结衣av一区二区av| 久久久久久久久久久久大奶| 中文字幕av电影在线播放| 久久久国产欧美日韩av| 久久免费观看电影| 亚洲精品av麻豆狂野| 天天躁狠狠躁夜夜躁狠狠躁| 男人操女人黄网站| 9热在线视频观看99| 欧美黄色片欧美黄色片| av又黄又爽大尺度在线免费看| 中文字幕精品免费在线观看视频| 亚洲精品美女久久av网站| 亚洲国产av影院在线观看| www日本在线高清视频| 人人妻人人爽人人添夜夜欢视频| 中文字幕色久视频| 美女视频免费永久观看网站| 51午夜福利影视在线观看| 精品第一国产精品| 亚洲国产看品久久| 一区二区av电影网| 一二三四在线观看免费中文在| 国产欧美亚洲国产| 正在播放国产对白刺激| 99国产极品粉嫩在线观看| 男女之事视频高清在线观看| 久久99热这里只频精品6学生| 日本五十路高清| 欧美精品一区二区大全| 久久 成人 亚洲| 国产不卡一卡二| 男人舔女人的私密视频| 高清毛片免费观看视频网站 | 亚洲国产av新网站| 九色亚洲精品在线播放| 黑丝袜美女国产一区| 国产精品免费大片| 精品人妻1区二区| 中国美女看黄片| 亚洲av成人一区二区三| 亚洲第一青青草原| 国产极品粉嫩免费观看在线| 午夜激情av网站| 亚洲精品美女久久久久99蜜臀| 狂野欧美激情性xxxx| 变态另类成人亚洲欧美熟女 | 亚洲av成人一区二区三| 99re6热这里在线精品视频| 免费不卡黄色视频| 建设人人有责人人尽责人人享有的| 欧美日韩中文字幕国产精品一区二区三区 | av片东京热男人的天堂| 在线观看www视频免费| 香蕉丝袜av| 久久久久久久久免费视频了| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线| 国产亚洲一区二区精品| 大陆偷拍与自拍| 久久人妻av系列| 男女午夜视频在线观看| 国产av又大| 欧美激情久久久久久爽电影 | 免费av中文字幕在线| 黑人巨大精品欧美一区二区蜜桃| 丝袜人妻中文字幕| 欧美国产精品va在线观看不卡| 国产精品二区激情视频| 高清在线国产一区| 巨乳人妻的诱惑在线观看| 制服诱惑二区| 国产精品亚洲一级av第二区| cao死你这个sao货| 黑人操中国人逼视频| 久久av网站| 嫁个100分男人电影在线观看| 一边摸一边抽搐一进一小说 | 一本—道久久a久久精品蜜桃钙片| 午夜91福利影院| 亚洲全国av大片| 国产三级黄色录像| 男人操女人黄网站| 久久久久久久久免费视频了| 亚洲精品在线美女| 国产欧美日韩精品亚洲av| 久久久国产成人免费| 国产精品电影一区二区三区 | 在线亚洲精品国产二区图片欧美| 深夜精品福利| 日本一区二区免费在线视频| 少妇猛男粗大的猛烈进出视频| 午夜视频精品福利| 十八禁网站网址无遮挡| 久9热在线精品视频| 午夜91福利影院| 深夜精品福利| 欧美大码av| 免费黄频网站在线观看国产| 18禁黄网站禁片午夜丰满| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| 丰满少妇做爰视频| 午夜福利免费观看在线| 搡老乐熟女国产| 国产xxxxx性猛交| 国产精品美女特级片免费视频播放器 | 性色av乱码一区二区三区2| 黄色a级毛片大全视频| 精品国产一区二区三区四区第35| 极品人妻少妇av视频| 国产亚洲一区二区精品| a在线观看视频网站| 在线十欧美十亚洲十日本专区| 久久久国产精品麻豆| 淫妇啪啪啪对白视频| 午夜福利欧美成人| 一本—道久久a久久精品蜜桃钙片| 成人av一区二区三区在线看| 国产精品成人在线| 麻豆成人av在线观看| 在线看a的网站| netflix在线观看网站| 国产在线视频一区二区| 亚洲精品国产色婷婷电影| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 日韩制服丝袜自拍偷拍| 最黄视频免费看| 一级黄色大片毛片| 男女边摸边吃奶| 免费日韩欧美在线观看| 精品国产乱码久久久久久男人| 亚洲精品av麻豆狂野| 亚洲,欧美精品.| 午夜激情久久久久久久| 精品国产一区二区三区四区第35| 精品午夜福利视频在线观看一区 | 中国美女看黄片| 免费日韩欧美在线观看| 亚洲七黄色美女视频| 丰满少妇做爰视频| 99国产综合亚洲精品| 午夜福利视频精品| 美女高潮到喷水免费观看| 最新美女视频免费是黄的| 午夜福利乱码中文字幕| 人妻久久中文字幕网| 久久久国产欧美日韩av| 日韩大片免费观看网站| 丝袜喷水一区| 老司机午夜福利在线观看视频 | tocl精华| 精品国产超薄肉色丝袜足j| 午夜福利在线观看吧| 欧美性长视频在线观看| 久久99一区二区三区| 成人永久免费在线观看视频 | 777久久人妻少妇嫩草av网站| 国产亚洲精品久久久久5区| 丝袜美足系列| 日韩视频一区二区在线观看| 丰满少妇做爰视频| 亚洲av电影在线进入| 久久青草综合色| 高清视频免费观看一区二区| 黑人欧美特级aaaaaa片| a级片在线免费高清观看视频| 婷婷丁香在线五月| 少妇精品久久久久久久| 亚洲成国产人片在线观看| 国产在线一区二区三区精| 欧美日韩福利视频一区二区| 好男人电影高清在线观看| 欧美日韩一级在线毛片| 亚洲精品国产区一区二| 成年人黄色毛片网站| 亚洲三区欧美一区| tocl精华| 国产欧美日韩一区二区三区在线| 窝窝影院91人妻| 亚洲精品美女久久久久99蜜臀| 久久国产精品人妻蜜桃| 一区二区日韩欧美中文字幕| 操出白浆在线播放| 国产精品熟女久久久久浪| 日韩欧美一区视频在线观看| 亚洲少妇的诱惑av| 国产欧美日韩一区二区三区在线| 亚洲精品自拍成人| 老司机亚洲免费影院| 成年人免费黄色播放视频| 老司机福利观看| 亚洲 国产 在线| 高清毛片免费观看视频网站 | 欧美激情久久久久久爽电影 | 精品国产乱码久久久久久小说| 欧美日韩国产mv在线观看视频| 亚洲精品美女久久久久99蜜臀| 成人18禁在线播放| av天堂在线播放| 午夜福利免费观看在线| 黑人猛操日本美女一级片| av天堂在线播放| 日韩大码丰满熟妇| 脱女人内裤的视频| 午夜免费成人在线视频| 这个男人来自地球电影免费观看| 欧美乱码精品一区二区三区| 亚洲,欧美精品.| 中文字幕色久视频| 老汉色∧v一级毛片| 交换朋友夫妻互换小说| 欧美大码av| 一个人免费看片子| 波多野结衣一区麻豆| av不卡在线播放| 精品亚洲成a人片在线观看| 高清黄色对白视频在线免费看| 精品少妇久久久久久888优播| 日本欧美视频一区| 香蕉国产在线看| 国产成人av激情在线播放| 少妇粗大呻吟视频| 汤姆久久久久久久影院中文字幕| 欧美国产精品va在线观看不卡| 美女高潮到喷水免费观看| 美女福利国产在线| 国产亚洲精品久久久久5区| 久久精品亚洲熟妇少妇任你| av超薄肉色丝袜交足视频| 9热在线视频观看99| 一级,二级,三级黄色视频| 欧美精品亚洲一区二区| 成人手机av| 悠悠久久av| av福利片在线| 国产精品影院久久| 欧美一级毛片孕妇| 他把我摸到了高潮在线观看 | 免费高清在线观看日韩| 超碰成人久久| 1024视频免费在线观看| 999精品在线视频| 精品国产乱码久久久久久小说| www.自偷自拍.com| 亚洲精品乱久久久久久| 人妻 亚洲 视频| 制服人妻中文乱码| 国产欧美日韩精品亚洲av| 啦啦啦视频在线资源免费观看| 国产三级黄色录像| 夜夜骑夜夜射夜夜干| 国产福利在线免费观看视频| 久久国产精品影院| 夜夜爽天天搞| 午夜两性在线视频| 美女主播在线视频| 国产亚洲精品第一综合不卡| 在线av久久热| av欧美777| 欧美大码av| 国产一区二区三区综合在线观看| 中文字幕高清在线视频| 少妇裸体淫交视频免费看高清 | 欧美精品av麻豆av| 黄色视频不卡| 亚洲成人免费电影在线观看| 高清欧美精品videossex| 一进一出好大好爽视频| 中文字幕人妻熟女乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻,人人澡人人爽秒播| 欧美精品一区二区免费开放| 777久久人妻少妇嫩草av网站| 欧美激情高清一区二区三区| 亚洲精品久久午夜乱码| 亚洲欧美一区二区三区黑人| 99国产精品99久久久久| 亚洲九九香蕉| 欧美日韩av久久| 亚洲av欧美aⅴ国产| 亚洲三区欧美一区| 精品熟女少妇八av免费久了| 新久久久久国产一级毛片| 午夜激情久久久久久久| 1024香蕉在线观看| 我的亚洲天堂| 岛国毛片在线播放| netflix在线观看网站| 久久99一区二区三区| 国产精品亚洲一级av第二区| 国产熟女午夜一区二区三区| av线在线观看网站| 91麻豆av在线| 国产91精品成人一区二区三区 | 国产精品自产拍在线观看55亚洲 | 成年动漫av网址| 国产aⅴ精品一区二区三区波| 99国产精品免费福利视频| 大片免费播放器 马上看| 天堂中文最新版在线下载| 变态另类成人亚洲欧美熟女 | 男人舔女人的私密视频| 两人在一起打扑克的视频| 男人舔女人的私密视频| 一本色道久久久久久精品综合| 国产精品香港三级国产av潘金莲| 两个人看的免费小视频| 大香蕉久久网| 两性午夜刺激爽爽歪歪视频在线观看 | 久久亚洲精品不卡| 制服人妻中文乱码| 日韩熟女老妇一区二区性免费视频| 我要看黄色一级片免费的| 国产精品 国内视频| 老汉色∧v一级毛片| av片东京热男人的天堂| 久久久精品免费免费高清| 91老司机精品| 天天添夜夜摸| 水蜜桃什么品种好| 久久性视频一级片| 国产不卡av网站在线观看| 中文字幕av电影在线播放| 久久人妻熟女aⅴ| 下体分泌物呈黄色| 亚洲天堂av无毛| 免费黄频网站在线观看国产| 欧美一级毛片孕妇| 一本一本久久a久久精品综合妖精| 色播在线永久视频| 老鸭窝网址在线观看| 999久久久精品免费观看国产| 日本撒尿小便嘘嘘汇集6| 一区二区日韩欧美中文字幕| 啦啦啦中文免费视频观看日本| 少妇裸体淫交视频免费看高清 | 两个人看的免费小视频| 免费高清在线观看日韩| 九色亚洲精品在线播放| 欧美日韩福利视频一区二区| 99re在线观看精品视频| 人人妻人人添人人爽欧美一区卜| 久久国产精品大桥未久av| 国产精品1区2区在线观看. | 久久久国产一区二区| 99国产精品免费福利视频| 一本综合久久免费| a在线观看视频网站| av线在线观看网站| 色尼玛亚洲综合影院| 人人妻人人爽人人添夜夜欢视频| 搡老熟女国产l中国老女人| 久久天躁狠狠躁夜夜2o2o| 好男人电影高清在线观看| 日韩中文字幕视频在线看片| 亚洲国产看品久久| 免费av中文字幕在线| 国产午夜精品久久久久久| 欧美 亚洲 国产 日韩一| 91国产中文字幕| 免费看a级黄色片| 久久久精品94久久精品| 一个人免费在线观看的高清视频| 黄色成人免费大全| 久久这里只有精品19| 久久毛片免费看一区二区三区| 日日摸夜夜添夜夜添小说| 久久人妻熟女aⅴ| 国产无遮挡羞羞视频在线观看| 国产成人欧美在线观看 | 久久久国产成人免费| 国产成人系列免费观看| 精品少妇久久久久久888优播| 亚洲欧美一区二区三区黑人| 久久婷婷成人综合色麻豆| 可以免费在线观看a视频的电影网站| 人人妻人人爽人人添夜夜欢视频| 老司机深夜福利视频在线观看| 久久精品亚洲av国产电影网| 欧美精品高潮呻吟av久久| 国产一区二区 视频在线| 日韩欧美免费精品| 丰满饥渴人妻一区二区三| 极品人妻少妇av视频| 亚洲成人国产一区在线观看| 在线观看舔阴道视频| 国产精品偷伦视频观看了| 国产精品麻豆人妻色哟哟久久| 中文字幕最新亚洲高清| 精品国产乱码久久久久久男人| 热re99久久精品国产66热6| 成年女人毛片免费观看观看9 | 亚洲七黄色美女视频| 欧美 日韩 精品 国产| 亚洲avbb在线观看| 色94色欧美一区二区| 在线观看免费视频日本深夜| 日韩一卡2卡3卡4卡2021年| 亚洲第一青青草原| 亚洲va日本ⅴa欧美va伊人久久| 一区在线观看完整版| 欧美在线黄色| 日日夜夜操网爽| 久久天躁狠狠躁夜夜2o2o| 国产男靠女视频免费网站| 青青草视频在线视频观看| 高清黄色对白视频在线免费看| 日韩免费av在线播放| av一本久久久久| 日韩欧美免费精品| 欧美 日韩 精品 国产| 国产精品欧美亚洲77777| 亚洲性夜色夜夜综合| 不卡一级毛片| 欧美中文综合在线视频| 国产av精品麻豆| 精品一区二区三区视频在线观看免费 | 免费在线观看完整版高清| 亚洲va日本ⅴa欧美va伊人久久| 精品熟女少妇八av免费久了| 他把我摸到了高潮在线观看 | 黄片大片在线免费观看| 国产日韩欧美在线精品| 窝窝影院91人妻| 亚洲av国产av综合av卡| 国产一区有黄有色的免费视频| 精品少妇久久久久久888优播| 老鸭窝网址在线观看| 国产野战对白在线观看| 99精品在免费线老司机午夜| 中国美女看黄片| 欧美av亚洲av综合av国产av| 国产成人免费无遮挡视频| 1024视频免费在线观看| 欧美精品高潮呻吟av久久| 真人做人爱边吃奶动态| 亚洲av国产av综合av卡| 亚洲第一av免费看| 丝袜美足系列| 欧美成狂野欧美在线观看| 欧美黑人精品巨大| av电影中文网址| 欧美日韩福利视频一区二区| 国产91精品成人一区二区三区 | 国产欧美日韩一区二区精品| 久久精品aⅴ一区二区三区四区| 一本久久精品| 波多野结衣一区麻豆|