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

    基于高分五號DPC氧氣A吸收波段的云頂壓強(qiáng)反演

    2021-06-08 03:54:16伽麗麗李正強(qiáng)朱思峰謝一凇洪津代海山麻金繼
    關(guān)鍵詞:云頂頂角波段

    喬 瑞,伽麗麗,許 華,李正強(qiáng)?,朱思峰,謝一凇,洪津,代海山,麻金繼

    (1中國科學(xué)院空天信息創(chuàng)新研究院,國家環(huán)境保護(hù)衛(wèi)星遙感重點(diǎn)實(shí)驗(yàn)室,北京 100101;2中國科學(xué)院大學(xué),北京 100049;3中國科學(xué)院合肥物質(zhì)科學(xué)研究院安徽光學(xué)精密機(jī)械研究所,中國科學(xué)院通用光學(xué)定標(biāo)與表征技術(shù)重點(diǎn)實(shí)驗(yàn)室,安徽 合肥 230031;4上海衛(wèi)星工程研究所,上海 201109;5安徽師范大學(xué)地理與旅游學(xué)院,安徽 蕪湖 241000)

    0 引言

    云大約覆蓋了海洋面積的71.7%和陸地面積的58.4%[1],是影響全球輻射收支的關(guān)鍵因素,也是氣候變化中最大的不確定性源[2]。云頂壓強(qiáng)與云頂高度具有較穩(wěn)定的對應(yīng)關(guān)系,是重要的云物理參數(shù)之一。準(zhǔn)確獲取云頂壓強(qiáng)對于天氣預(yù)報(bào)尤其是強(qiáng)對流天氣的監(jiān)測預(yù)警有著重要的作用[3]。此外,利用云頂壓強(qiáng)對一些較難識別的云(如薄云、高亮地表(包括冰雪)上空的云)有較好的識別效果。云頂壓強(qiáng)還可作為一些高度特征明顯的云(如卷云、深對流云等)分類的重要判據(jù)。

    利用衛(wèi)星數(shù)據(jù)反演云頂壓強(qiáng),主要有成像幾何和通道輻射特性兩類方法[4]。成像幾何方法主要是利用不同衛(wèi)星或同一衛(wèi)星的不同角度對于同一目標(biāo)的視差和衛(wèi)星的位置來反演云頂高度,并轉(zhuǎn)化為云頂壓強(qiáng)[5]。而通道輻射特性方法主要有單通道或多通道亮溫法[6]、分裂窗法[7,8]、CO2切片法[9,10]、偏振反射率法[11,12]和氧氣A吸收帶法等。

    氧氣A吸收帶位于758~778 nm,此波長范圍除氧氣外幾乎沒有其他氣體吸收,而氧氣在大氣中的含量比例穩(wěn)定,所以該吸收帶是一個(gè)反演云頂壓強(qiáng)的理想通道[13]。利用氧氣A吸收帶反演云頂壓強(qiáng)的想法最早由Yamamoto和Wark提出[14]。Kuze和Chance[15]利用SCIAMACHY的觀測值和模擬計(jì)算的氧氣吸收率通過最小二乘擬合獲取了云頂高度。Br′eon等[16]利用POLDER傳感器在氧氣A吸收帶的寬窄兩個(gè)波段進(jìn)行了表面壓強(qiáng)反演的理論分析,其利用輻射傳輸模型進(jìn)行正向模擬,得到兩波段比值與壓強(qiáng)的經(jīng)驗(yàn)關(guān)系式,在POLDER云頂壓強(qiáng)產(chǎn)品中得到了應(yīng)用[17]。Preusker等[18]利用印度遙感3號衛(wèi)星MOS傳感器761 nm和765 nm波段的輻亮度比值計(jì)算了云頂高度。張巖等[19]對SCIAMACHY的兩種經(jīng)典云頂高度算法SACULA和FRESCO+的反演結(jié)果進(jìn)行了比較,并用地基雷達(dá)和探空數(shù)據(jù)進(jìn)行了驗(yàn)證。

    高分五號(GF-5)衛(wèi)星搭載的多角度偏振探測儀(Directional polarimetric camera,DPC),是中國第一個(gè)以氣溶膠探測為目標(biāo)的多角度偏振星載傳感器[20]。GF-5位于太陽同步軌道,過境時(shí)間為當(dāng)?shù)貢r(shí)間13:30,重訪周期為2天。DPC的波段設(shè)置和POLDER[21]相近,包括可見光至近紅外的共8個(gè)探測波段(443,490,565,670,763,765,865,910 nm),其中490、670、865 nm波段為偏振測量波段,763 nm和765 nm波段位于氧氣A吸收帶上,帶寬分別為10 nm和40 nm。DPC通過在衛(wèi)星飛行過程中沿軌快速框幅式攝影,實(shí)現(xiàn)每個(gè)地面像元至少9個(gè)角度的觀測。

    本文在POLDER反演云頂壓強(qiáng)[16]方法的基礎(chǔ)上,利用DPC的763 nm和765 nm兩個(gè)氧氣A吸收波段,研究DPC云頂壓強(qiáng)反演算法。主要考慮了DPC兩氧氣A吸收波段的儀器光譜響應(yīng)函數(shù),擴(kuò)展了輻射傳輸模擬樣本的海拔覆蓋范圍,改進(jìn)了兩氧氣A波段反射率比值與表面壓強(qiáng)經(jīng)驗(yàn)關(guān)系的擬合公式,并對云頂壓強(qiáng)反演進(jìn)行了初步驗(yàn)證。文中,首先介紹了氧氣A吸收帶反演云頂壓強(qiáng)的原理;然后,通過輻射傳輸模擬,得到DPC兩個(gè)氧氣A波段表觀輻亮度比值反演云頂壓強(qiáng)的擬合公式,校正了擬合誤差,并簡要分析了影響反演結(jié)果的其他誤差因素,如大氣氣溶膠、大氣廓線等;進(jìn)一步,利用晴空地表對反演方法進(jìn)行了穩(wěn)定性分析和原理性驗(yàn)證(對比DEM),并比較了DPC云頂壓強(qiáng)反演結(jié)果與MODIS云頂壓強(qiáng)產(chǎn)品的一致性。最終,得到了分辨率和精度均較高的DPC云頂壓強(qiáng)反演結(jié)果。

    1 氧氣A吸收帶反演云頂壓強(qiáng)的原理

    氧氣A吸收帶波長范圍為758~778 nm,是氧氣吸收最強(qiáng)的一個(gè)吸收帶,且在該吸收帶內(nèi)幾乎沒有其他大氣吸收,比較“純凈”。GF-5搭載的DPC在氧氣A吸收帶設(shè)置了兩個(gè)探測波段:一個(gè)是中心波長為763 nm的窄波段,覆蓋了大部分較強(qiáng)的氧氣吸收譜線;另一個(gè)是中心波長為765 nm的寬波段,覆蓋氧氣A吸收帶的同時(shí)還覆蓋了部分非吸收波段。因此,在傳感器實(shí)際觀測中,765 nm波段測得的表觀輻亮度高于763 nm波段。圖1為兩個(gè)波段的光譜響應(yīng)函數(shù)以及6S模擬的該波段整層大氣氧氣吸收透過率。

    圖1 DPC/GF-5 763 nm和765 nm波段傳感器光譜響應(yīng)函數(shù)及6S模擬的整層大氣氧氣吸收透過率Fig.1 Sensor spectral response and the total atmosphere oxygen transmittance simulated by 6S in DPC/GF-5 763nm and 765nm bands

    某一波段氧氣的吸收系數(shù)可以利用逐線積分法計(jì)算得到,假定波段內(nèi)的每條譜線都是Lorentz型[22],則該波段的吸收系數(shù)k的計(jì)算公式為

    式中:N為該波段范圍包含的吸收譜線個(gè)數(shù);Si為第i條吸收譜線的線強(qiáng),是溫度T的函數(shù);ai為第i條吸收譜線的半寬,是氣壓P和溫度T的函數(shù);ν和ν0,i分別是該波段中心波數(shù)和第i條譜線的中心波數(shù)。

    單一吸收譜線的線強(qiáng)Si和線寬ai的計(jì)算公式為

    式中:S0為標(biāo)準(zhǔn)溫度T0的線強(qiáng),E為吸收光譜對應(yīng)的電子躍遷低能級能量,a0為標(biāo)準(zhǔn)氣壓P0和標(biāo)準(zhǔn)溫度T0下的半寬,這三個(gè)參數(shù)均為常數(shù),對應(yīng)的標(biāo)準(zhǔn)氣壓P0=1013 hPa、標(biāo)準(zhǔn)溫度T0=273 K。

    從上述公式可以看出,對于特定的波段,氧氣的吸收透過率(可以由吸收系數(shù)計(jì)算得到)可以看作大氣壓強(qiáng)P和溫度T的函數(shù),而對于一個(gè)確定的大氣溫度廓線,氧氣的透過率只和大氣壓強(qiáng)P有關(guān)。這樣,可以利用氧氣吸收影響權(quán)重不同的763 nm和765 nm波段觀測值之間的對比關(guān)系,反演云頂(或晴空地表)大氣壓強(qiáng)P。

    2 云頂壓強(qiáng)反演方法

    2.1 基于輻射傳輸模擬確定反演擬合公式

    對于不透明的云,大氣頂表觀反射率可以認(rèn)為大致由瑞利散射、氣溶膠散射和云頂反射三部分構(gòu)成,即

    式中:tatm是云上大氣透過率,在DPC的763 nm和765 nm波段內(nèi),該透過率主要受氧氣吸收影響,其他大氣分子的吸收如臭氧、水汽影響非常小(見表1)。

    表1 6S模擬的DPC/GF-5 763 nm和765 nm波段吸收氣體整層大氣吸收透過率t763和t765Table 1 The total atmosphere transmittance t763and t765of absorption gases simulated by 6S in DPC/GF-5 763nm and 765nm bands

    瑞利散射貢獻(xiàn)Rrayleigh可以精確計(jì)算,其大小依賴于高度,在反演前是未知量,但由于這兩個(gè)波段瑞利散射光學(xué)厚度很小(約0.02),該項(xiàng)對RTOA貢獻(xiàn)也較小。氣溶膠散射貢獻(xiàn)Raerosol值依賴于云上氣溶膠含量,在云層較高的情況下,該項(xiàng)對于RTOA的影響較小。云頂反射Rcloud難以精確計(jì)算,因?yàn)橛?jì)算云層透過及內(nèi)部多次散射的部分比較困難。對于本研究,由于兩個(gè)波段非常接近,可以認(rèn)為Rrayleigh、Raerosol和Rcloud在兩個(gè)波段是相等的。

    利用以上特性,定義比例系數(shù)X為763 nm與765 nm波段表觀反射率的比值,其表達(dá)式為

    根據(jù)上述分析,X≈tO2(763nm)/tO2(765nm),而不同高度大氣中氧氣比例恒定約為21%,氧氣透過率tO2可以等價(jià)于云頂壓強(qiáng)P,所以X與P存在函數(shù)關(guān)系。

    輻射傳輸模擬時(shí)還需要考慮光程影響,定義大氣光學(xué)質(zhì)量m=1/μs+1/μν表示不同角度入射和反射光路的光程差異,μs為太陽天頂角的余弦,μν為觀測天頂角的余弦。

    為了得到X與P的函數(shù)關(guān)系,不考慮云層透過及內(nèi)部多次散射的影響,利用輻射傳輸模型對比例系數(shù)X進(jìn)行比較全面的模擬,壓強(qiáng)P范圍為1013.2~75.6 hPa(對應(yīng)海拔高度0~18 km,步長1 km);大氣光學(xué)質(zhì)量m范圍為2~5.84(對應(yīng)太陽天頂角0~70°,觀測天頂角0~70°,步長10°);云頂反射率設(shè)置為0.7;氣溶膠光學(xué)厚度(Aerosol optical depth,AOD)值設(shè)置為0.1;大氣模式設(shè)置為熱帶大氣。

    模擬結(jié)果如圖2所示。由圖2(a)可以看出mP與X有趨勢上的關(guān)系,但比較離散,無法擬合出關(guān)系曲線??紤]到一定高度范圍的氧氣吸收光學(xué)厚度與氣壓存在近似的平方關(guān)系[23],故可以擬合mP2與X的關(guān)系曲線,如圖2(b)所示。為了減少擬合偏差,嘗試多種擬合模型后,選擇用4次多項(xiàng)式進(jìn)行擬合。擬合公式為

    圖2 763 nm和765 nm波段表觀反射率比例系數(shù)X與大氣光學(xué)質(zhì)量m和大氣壓強(qiáng)P的關(guān)系。(a)X與mP關(guān)系;(b)X與mP2關(guān)系Fig.2 The function between X(ratio of 763 nm and 765 nm TOA)and m(air mass),P(atmospheric pressure).(a)Function of X and mP,(b)function of X and mP2

    擬合相關(guān)系數(shù)R2=0.9995,說明X與mP2有非常強(qiáng)的相關(guān)性。擬合參數(shù)見表2。

    表2 擬合公式(7)各項(xiàng)參數(shù)值Table 2 Fitting formula(7)parameters

    2.2 擬合誤差分析與校正

    由于X與mP2的數(shù)量級相差較大,擬合誤差的影響不可忽視,下面進(jìn)行擬合誤差分析和校正。

    以太陽和觀測天頂角均為0°、60°、70°為例,對于不同的X,用擬合公式得到mP2,計(jì)算壓強(qiáng)誤差相對值ΔP/P,得到誤差趨勢圖(圖3),其他天頂角的誤差趨勢與圖3一致,其值介于0°和70°的樣本曲線之間(類似圖 3 中 60°樣本)。

    圖3 不同太陽和觀測天頂角(0°,60°,70°)壓強(qiáng)擬合相對誤差Fig.3 Relative error of pressure fitting at different solar zenith angles and view zenith angles(0°,60°,70°)

    從圖3可以看出,擬合誤差主要和兩個(gè)因素有關(guān),一個(gè)是X的大小,另一個(gè)是太陽與觀測天頂角。對于X,在X<0.9時(shí),壓強(qiáng)的擬合誤差都較小,在±2%以內(nèi);而0.90.95時(shí),誤差均為正,且快速增大。對于角度,天頂角較小時(shí)擬合誤差較大,而天頂角較大時(shí)擬合誤差顯著變小。這主要是因?yàn)榻嵌鹊脑龃笤黾恿斯獬?限制了X的增大。對反演來說,由擬合產(chǎn)生的誤差主要會(huì)對云頂壓強(qiáng)較小的情況(小于182 hPa,對應(yīng)海拔高度約13 km以上)產(chǎn)生影響,而對于云頂壓強(qiáng)較大的情形,則對反演結(jié)果影響較小。

    因此,擬合誤差的校正要同時(shí)考慮X和太陽、觀測角度。對X>0.9的情況進(jìn)行校正,角度通過大氣光學(xué)質(zhì)量m反映。為了減少多次計(jì)算引起的誤差放大,直接對擬合項(xiàng)f(X)進(jìn)行校正。定義誤差校正量為Ecor,計(jì)算公式為

    式中,m0和m70分別對應(yīng)太陽和觀測天頂角均為0°或70°時(shí)的大氣光學(xué)質(zhì)量(其值分別為2和5.8476),中間角度m的值介于m0和m70之間。擬合誤差校正后P?X公式為

    重新計(jì)算得到新的誤差趨勢圖(圖4)??梢钥闯?經(jīng)過擬合誤差校正以后,擬合誤差在X>0.9時(shí)得到了顯著降低,基本都小于1.5%,提高了反演的準(zhǔn)確性。

    2.3 其他因素誤差分析

    實(shí)際反演時(shí),大氣和地表情況各異,與擬合公式的模擬條件有差異,進(jìn)一步分析由大氣氣溶膠、大氣廓線和地表反射幾方面因素帶來的反演誤差。

    1)大氣氣溶膠

    氣溶膠主要分布在海拔3 km以下的大氣中[24],對于具有一定高度且不透光的云來說,云上氣溶膠影響一般很小。氣溶膠對表觀反射率的影響比較復(fù)雜,有散射、吸收、多次散射等效應(yīng)??紤]到復(fù)雜性,利用輻射傳輸模擬計(jì)算其導(dǎo)致誤差最大的一種情況,即考慮整層大氣的影響。在實(shí)際反演中,誤差都小于該情況。分析不同AOD時(shí),利用式(11)反演云頂壓強(qiáng)的誤差,結(jié)果如圖5所示。

    圖5 不同AOD(0.01,0.3,0.5)時(shí)模擬反演的云頂壓強(qiáng)相對誤差Fig.5 Relative error of retrieved pressure for different AOD(0.01,0.3,0.5)

    圖5中,縱坐標(biāo)為不同AOD值(0.01,0.3,0.5)時(shí),基于擬合公式(11)(AOD值0.1)反演的云頂壓強(qiáng)相對誤差??梢钥闯?AOD大于0.1時(shí),云頂壓強(qiáng)反演結(jié)果偏大;AOD小于0.1時(shí),云頂壓強(qiáng)反演結(jié)果偏小;且相對誤差均不超過±0.2%。因此,氣溶膠含量對云頂壓強(qiáng)反演結(jié)果的影響可以忽略。

    2)大氣溫壓廓線

    利用擬合公式(11)模擬時(shí)選取的大氣模式為熱帶大氣,而地球上不同緯度和季節(jié)的大氣廓線會(huì)有一些差異,例如極地附近的對流層高度比熱帶低的多,冬季的大氣一般比夏季干燥等,尤其不同大氣模式的溫度壓強(qiáng)廓線也有一定的差別。本研究采用的763 nm和765 nm波段,除氧氣吸收外,其他氣體(如臭氧和水汽)的吸收都很小,所以其他氣體含量差異不會(huì)帶來很大的誤差。模擬五種不同的標(biāo)準(zhǔn)大氣模式:熱帶大氣(Tro)、中緯度夏季(MLS)、中緯度冬季(MLW)、亞北極區(qū)夏季(SAS)、亞北極區(qū)冬季(SAW)時(shí)763 nm和765 nm波段表觀反射率,分析不同大氣模式導(dǎo)致的云頂壓強(qiáng)反演誤差。

    圖6為其它四種大氣模式時(shí),基于擬合公式(11)(大氣模式為Tro)反演的壓強(qiáng)的相對誤差。可以看出,大氣模式差異會(huì)使壓強(qiáng)的反演結(jié)果偏大或偏小,如大氣模式為MLW時(shí),在太陽和觀測天頂角都為70°(m=5.8476)時(shí)壓強(qiáng)的相對誤差約為1.25%;大氣模式為SAW時(shí),在太陽和觀測天頂角都為0°(m=2)時(shí)壓強(qiáng)的相對誤差約為?1.55%??傮w來看,不同的大氣模式所導(dǎo)致的壓強(qiáng)反演誤差約在±1.5%以內(nèi),影響較小。

    圖6 不同大氣模式(MLS,MLW,SAS,SAW)時(shí)模擬反演的云頂壓強(qiáng)相對誤差Fig.6 Relative error of retrieved pressure for different atmosphere profiles(MLS,MLW,SAS,SAW)

    3)地表反射

    對于光學(xué)厚度較小的云,一部分光會(huì)穿過云層經(jīng)地表反射再被傳感器接收。這樣,反演得到的壓強(qiáng)會(huì)介于云頂壓強(qiáng)和地表壓強(qiáng)之間,造成結(jié)果偏大。一方面地表類型復(fù)雜多樣,另一方面薄云的性質(zhì)也難以描述,很難定量研究地表反射對于反演誤差的影響。因此,該算法對于不透光的厚云適用性更好。

    3 結(jié)果與分析

    3.1 利用晴空地表驗(yàn)證算法原理

    本方法利用氧氣A吸收帶寬窄兩個(gè)波段的比值來反演云頂壓強(qiáng)(或表面壓強(qiáng)),理論上可以用于任何下墊面(包括云頂或晴空地表)壓強(qiáng)的反演[16]。因此,本小節(jié)采用晴空地表的觀測對云頂壓強(qiáng)(表面壓強(qiáng))反演算法進(jìn)行原理性驗(yàn)證。

    3.1.1 反演穩(wěn)定性測試

    由于云具有三維結(jié)構(gòu)復(fù)雜且變化迅速的特點(diǎn),為驗(yàn)證本反演算法在不同太陽、觀測角度下的穩(wěn)定性,選擇DPC/GF-5同一晴空區(qū)域、不同觀測日期的多角度觀測數(shù)據(jù),反演得到晴空地表的表面壓強(qiáng),并評估反演結(jié)果的差異大小。選擇一個(gè)適合進(jìn)行穩(wěn)定性測試的區(qū)域:蘇萊曼山脈東部(緯度范圍30°N~31°N,經(jīng)度范圍69°E~70°E,海拔范圍933~2352 m)。該區(qū)域?yàn)樯降?地形復(fù)雜,能夠較好地檢驗(yàn)不同角度和日期情況下反演的穩(wěn)定性。

    圖7為2018年5月27日和2018年11月29日兩個(gè)不同日期對應(yīng)的9個(gè)不同觀測角度反演結(jié)果的差異系數(shù)(標(biāo)準(zhǔn)差與樣本均值的比值,CV),橫坐標(biāo)表示不同樣本??梢钥闯?大部分樣本點(diǎn)不同觀測角度反演結(jié)果的CV都小于4%。綜合兩個(gè)不同日期的平均CV為3.5%。個(gè)別樣本點(diǎn)CV稍大,超過5%,最大達(dá)7.6%,但樣本量較少,僅占總樣本數(shù)的3.6%。這可能是山地復(fù)雜的地形結(jié)構(gòu)導(dǎo)致DPC不同觀測角度觀測的地面目標(biāo)有差異而造成的。

    圖7 DPC不同日期表面壓強(qiáng)反演結(jié)果差異系數(shù)Fig.7 The coefficient of variation(CV)for DPC retrieved pressure from different dates

    3.1.2 基于DEM的表面壓強(qiáng)反演驗(yàn)證

    忽略不同氣壓帶和氣象因素對于氣壓的影響,海拔高度與大氣壓強(qiáng)有著穩(wěn)定的對應(yīng)關(guān)系。因此,可以利用該反演方法得到陸地晴空區(qū)域的地表壓強(qiáng),與利用DEM估算的地表壓強(qiáng)進(jìn)行對比,驗(yàn)證算法反演結(jié)果的可靠性。其中,DEM數(shù)據(jù)可以從DPC 1級數(shù)據(jù)中獲取,空間分辨率與DPC觀測數(shù)據(jù)一致。利用DEM數(shù)據(jù)估算地表壓強(qiáng)的公式為

    選取一個(gè)海拔高度變化范圍較大的地理區(qū)域(塔里木盆地和青藏高原交界處,緯度范圍36°N~37°N,經(jīng)度范圍82°E~83°E,海拔范圍1445~6286 m)為驗(yàn)證區(qū)域。數(shù)據(jù)觀測日期為2018年5月27日。圖8中橫坐標(biāo)為利用DEM通過高度-氣壓估算公式(公式12)估算的地表壓強(qiáng),縱坐標(biāo)為DPC氧氣A波段反演得到的地表壓強(qiáng)。

    圖8 DPC反演的表面壓強(qiáng)與DEM估算壓強(qiáng)對比結(jié)果Fig.8 The comparison between DPC retrieved atmospheric pressure and DEM estimated atmospheric pressure

    兩個(gè)地表壓強(qiáng)的線性擬合斜率k=1.115,截距b=?90.044,相關(guān)系數(shù)R=0.9421,說明DPC反演的地表壓強(qiáng)與DEM估算的地表壓強(qiáng)具有很好的一致性。經(jīng)統(tǒng)計(jì),二者均方根誤差(RMSE)約為47.6 hPa,最大偏差不超過100 hPa。該偏差大于氧A波段反演理論誤差的原因主要與DPC像元空間分辨率(3.3 km)較低有關(guān),在這個(gè)較大范圍內(nèi)地表復(fù)雜的地形結(jié)構(gòu)導(dǎo)致氧A帶與DEM估算對應(yīng)性之間存在較大的隨機(jī)誤差。但應(yīng)該強(qiáng)調(diào),兩者結(jié)果很好的線性擬合相關(guān)系數(shù)說明了氧A帶方法的有效性。

    3.2 全球反演結(jié)果

    圖9展示了DPC/GF-5 2018年5月27日全球云頂壓強(qiáng)反演結(jié)果。從圖中可以看出:云覆蓋了地球尤其是海洋的大部分區(qū)域,晴空多出現(xiàn)在干旱和內(nèi)陸地區(qū)。低緯度和中緯度的云形狀大多較碎;中高緯度的云多呈東西向的條帶狀,和地球風(fēng)帶有關(guān);高緯度的云多呈塊狀。云頂壓強(qiáng)較低即較高的云一般出現(xiàn)在低緯度的海陸交界區(qū)域附近,這是因?yàn)榈途暥葘α鲗痈叨容^高,而中緯度和高緯度要低的多,限制了云頂高度。

    圖9 2018年5月27日DPC云頂壓強(qiáng)全球反演結(jié)果樣例。(a)DPC全球圖像;(b)云頂壓強(qiáng)反演結(jié)果Fig.9 The sample of DPC global retrieved cloud top pressure on May 27,2018.(a)DPC global image,(b)retrieved cloud top pressure

    3.3 與MODIS云頂壓強(qiáng)產(chǎn)品的對比

    為了分析DPC云壓強(qiáng)(基于氧氣A帶法)空間分布結(jié)果的反演性能,將DPC的云頂壓強(qiáng)反演結(jié)果與MODIS云頂壓強(qiáng)產(chǎn)品(基于CO2切片法)對比,比較二者的差異。

    為去除DPC觀測數(shù)據(jù)中的晴空、太陽耀光和冰雪等區(qū)域,在反演過程中增加一些簡單的約束條件:1)要求海洋區(qū)域765 nm波段反射率大于0.08(大于該波段瑞利散射貢獻(xiàn)反射率的最大強(qiáng)度);2)要求陸地和海陸交界區(qū)域443、565、670 nm反射率差值小于0.01(利用云像元可見光各波段反射率基本相等的特性);3)剔除南半球海洋區(qū)域反演壓強(qiáng)大于950 hPa的像元(太陽耀光區(qū)域);4)剔除PDEM估算-P反演小于100 hPa的像元(冰雪覆蓋區(qū)域)??紤]云的復(fù)雜三維結(jié)構(gòu),這里選擇DPC觀測天頂角最小的觀測數(shù)據(jù)用于云頂壓強(qiáng)反演。

    選取MODIS/Aqua云頂壓強(qiáng)產(chǎn)品(MYD06)作為對比數(shù)據(jù)。MODIS/Aqua與DPC/GF-5過境時(shí)間相近,均為當(dāng)?shù)貢r(shí)間13:30,可以保證云的位置、形狀變化較小。MYD06云頂壓強(qiáng)產(chǎn)品利用MODIS紅外通道中第33~36波段通過CO2切片法反演得到,空間分辨率為5 km,對于卷云和中層水云其誤差小于50 hPa,其他云誤差小于200 hPa[25]。圖10展示了DPC和MODIS在2018年5月27日的云頂壓強(qiáng)對比結(jié)果。圖中從上到下每行依次對應(yīng)下墊面為:海洋(大西洋北部海域,緯度范圍32°N~46°N,經(jīng)度范圍55°W~45°W)、沙漠(阿拉伯半島附近,緯度范圍18°N~32°N,經(jīng)度范圍50°E~60°E)、植被(美國與加拿大交界處,緯度范圍40°N~54°N,經(jīng)度范圍110°W~100°W)。從圖10可以看出,DPC云頂壓強(qiáng)反演結(jié)果與MODIS產(chǎn)品的空間分布趨勢具有較好的一致性。且DPC的云頂壓強(qiáng)結(jié)果相比MODIS空間分布更為平滑,像元分辨率更高。此外,發(fā)現(xiàn)MODIS產(chǎn)品存在一些晴空像元的誤反演現(xiàn)象[如圖10(b)黑色方框所示,該區(qū)域明顯為晴空海面,卻反演出了“云頂壓強(qiáng)”]和一些條帶噪聲現(xiàn)象[如圖10(b)白色方框所示,這可能是由傳感器探測元正反掃描響應(yīng)差異、傳感器機(jī)械運(yùn)動(dòng)等造成的[26]]。而DPC云頂壓強(qiáng)結(jié)果則存在一些薄云和云邊緣區(qū)域的漏反演現(xiàn)象,這是由于氧氣A帶算法對厚云適用性更好,在數(shù)據(jù)處理時(shí)采用了相對嚴(yán)格的云識別判據(jù)。

    圖10 2018年5月27日DPC和MODIS的云頂壓強(qiáng)對比結(jié)果。(a)海洋區(qū)域;(b)沙漠區(qū)域;(c)植被區(qū)域Fig.10 Comparison of cloud top pressure between DPC and MODIS on May 27,2018.(a)Ocean region,(b)desert region,(c)vegetation region

    進(jìn)一步將DPC云頂壓強(qiáng)反演結(jié)果與MODIS/Aqua云頂壓強(qiáng)產(chǎn)品(MYD06)進(jìn)行定量對比。由于MODIS云頂壓強(qiáng)產(chǎn)品的空間分辨率為5 km,而DPC的空間分辨率3 km,為了盡量降低尺度因素的影響,同時(shí)避免邊緣效應(yīng),選取大面積云目標(biāo)的中間部分為樣本區(qū)域,并對DPC與MODIS數(shù)據(jù)進(jìn)行像元匹配。同時(shí),選取的樣本覆蓋高云、中云和低云,以增加對比結(jié)果的代表性。

    對比結(jié)果如圖11所示,圖中橫坐標(biāo)為MODIS產(chǎn)品的云頂壓強(qiáng),縱坐標(biāo)為利用DPC氧氣A波段觀測反演得到的云頂壓強(qiáng)。從圖中樣本橫坐標(biāo)可以看出,MODIS產(chǎn)品的壓強(qiáng)值不連續(xù),這是因?yàn)镸ODIS產(chǎn)品對壓強(qiáng)做了近似處理(最小單位為5 hPa)。DPC和MODIS的云頂壓強(qiáng)線性擬合斜率k=1.040,截距b=39.556,相關(guān)系數(shù)R=0.9567。這說明DPC反演結(jié)果與MODIS產(chǎn)品有很好的一致性。二者均方根誤差(RMSE)約為98.6 hPa,其中大部分的DPC反演結(jié)果較MODIS偏高。此外,在云頂壓強(qiáng)較低即云層較高時(shí),二者偏差較大??赡苁怯刹煌惴?氧氣A帶算法和CO2切片算法)在低壓強(qiáng)時(shí)敏感性不同,反演誤差增大造成的。

    圖11 DPC反演的云頂壓強(qiáng)與MODIS云頂壓強(qiáng)產(chǎn)品逐像元對比結(jié)果Fig.11 Pixel-by-pixel comparison between DPC retrieved cloud top pressure and MODIS cloud top pressure product

    4 結(jié)論

    基于DPC/GF-5的兩個(gè)氧氣A吸收波段(763 nm和765 nm,帶寬分別為10 nm和40 nm),研究了云頂壓強(qiáng)反演算法。首先,通過輻射傳輸模擬和多項(xiàng)式擬合得到了云頂壓強(qiáng)和兩波段表觀反射率比值的經(jīng)驗(yàn)關(guān)系,并對擬合公式進(jìn)行了誤差分析和校正,擬合公式整體誤差小于1.5%。然后,分析了大氣氣溶膠、大氣溫壓廓線等地氣條件對反演結(jié)果的影響,結(jié)果表明這些因素對云頂壓強(qiáng)反演結(jié)果影響均較小,可忽略。進(jìn)一步利用DPC晴空地表觀測數(shù)據(jù),對反演方法進(jìn)行了穩(wěn)定性測試和原理性驗(yàn)證。穩(wěn)定性測試結(jié)果表明該算法在不同太陽、觀測角度條件下,反演結(jié)果差異在5%以內(nèi),具有較好的穩(wěn)定性。晴空地表DPC反演的表面壓強(qiáng)與DEM估算的壓強(qiáng)對比的平均偏差約為47.6 hPa。最后,將DPC反演的云頂壓強(qiáng)與MODIS云頂壓強(qiáng)產(chǎn)品(MYD06)對比,結(jié)果表明二者有很好的一致性,云頂壓強(qiáng)平均偏差約為98.6 hPa。

    本反演方法理論上對于云和晴空地表均適用,云頂壓強(qiáng)結(jié)果不但可用于云高估計(jì),還可用于云識別和云分類。由于忽略云三維形狀、云的透過及云層間多次散射的影響,本反演對于不透光且充滿整個(gè)像元的厚云誤差較小,約為幾十hPa;而對于薄云或未充滿整個(gè)像元的云邊緣區(qū)和碎云則誤差較大,反演結(jié)果相對真實(shí)值偏大。未來可通過增加氧A帶光譜分辨率、傳感器像元空間分辨率等途徑進(jìn)一步改進(jìn)云頂壓強(qiáng)產(chǎn)品的精度和適用性。

    猜你喜歡
    云頂頂角波段
    春日暖陽
    王麗敏作品
    大眾文藝(2022年17期)2022-09-27 16:05:38
    一般三棱鏡最大頂角與折射率的關(guān)系
    FY-3D MERSI-II云頂產(chǎn)品算法及精度檢驗(yàn)
    涼亭中的數(shù)學(xué)
    頂角為100°的等腰三角形性質(zhì)的應(yīng)用
    M87的多波段輻射過程及其能譜擬合
    日常維護(hù)對L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    L波段雷達(dá)磁控管的使用與維護(hù)
    河南科技(2014年18期)2014-02-27 14:14:53
    雙元計(jì)算法求體積本征值
    美女黄网站色视频| 午夜久久久久精精品| 免费观看精品视频网站| 老司机影院毛片| av在线观看视频网站免费| 久久久久久久久久黄片| 日韩高清综合在线| 日韩欧美国产在线观看| 在线观看66精品国产| 五月伊人婷婷丁香| 中国国产av一级| 国产亚洲精品久久久com| 久久草成人影院| 色尼玛亚洲综合影院| 亚洲av成人精品一区久久| 夜夜爽夜夜爽视频| 尾随美女入室| 黄色配什么色好看| 免费av不卡在线播放| 两个人的视频大全免费| 国产 一区精品| 男人舔女人下体高潮全视频| 国产精品av视频在线免费观看| 国产成人91sexporn| 国产av不卡久久| 精品久久久久久久久av| or卡值多少钱| 亚洲av一区综合| 免费大片18禁| 成人午夜高清在线视频| 久久6这里有精品| 久久久久九九精品影院| 男人的好看免费观看在线视频| 99久久精品国产国产毛片| 最近2019中文字幕mv第一页| 久久久国产成人精品二区| 高清毛片免费看| 久久久精品大字幕| 午夜免费男女啪啪视频观看| 国产激情偷乱视频一区二区| 亚洲丝袜综合中文字幕| 少妇的逼水好多| 最近最新中文字幕大全电影3| 国产在视频线在精品| 成人二区视频| 色尼玛亚洲综合影院| 十八禁国产超污无遮挡网站| 99久久成人亚洲精品观看| 夜夜爽夜夜爽视频| 国产成人一区二区在线| 亚洲经典国产精华液单| 国产免费一级a男人的天堂| 欧美日韩在线观看h| 国产亚洲最大av| 国产黄色视频一区二区在线观看 | 美女国产视频在线观看| 国产亚洲av嫩草精品影院| 国产在视频线在精品| 最近中文字幕2019免费版| 精品一区二区免费观看| 日韩av在线大香蕉| 国产色爽女视频免费观看| 日本免费在线观看一区| 高清av免费在线| 亚洲欧美日韩高清专用| 中文字幕av在线有码专区| 日韩一本色道免费dvd| 精品久久国产蜜桃| 色综合站精品国产| a级毛色黄片| 国产伦一二天堂av在线观看| 夜夜看夜夜爽夜夜摸| 中文天堂在线官网| 人人妻人人澡欧美一区二区| 热99在线观看视频| 村上凉子中文字幕在线| 日韩欧美精品v在线| 日日撸夜夜添| 亚洲精品国产成人久久av| 五月玫瑰六月丁香| 丰满乱子伦码专区| 久久久久九九精品影院| 三级毛片av免费| av线在线观看网站| 小蜜桃在线观看免费完整版高清| 成人高潮视频无遮挡免费网站| 亚洲成人av在线免费| 国内少妇人妻偷人精品xxx网站| 草草在线视频免费看| 亚洲精品影视一区二区三区av| 欧美成人午夜免费资源| 大话2 男鬼变身卡| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看 | 日本爱情动作片www.在线观看| 又爽又黄无遮挡网站| 少妇熟女aⅴ在线视频| 色噜噜av男人的天堂激情| 国产成人a∨麻豆精品| 色尼玛亚洲综合影院| 国产av码专区亚洲av| 黑人高潮一二区| 99在线视频只有这里精品首页| 久久这里有精品视频免费| 日日啪夜夜撸| 一个人观看的视频www高清免费观看| 久久精品国产亚洲av天美| 欧美日韩综合久久久久久| 免费电影在线观看免费观看| 美女cb高潮喷水在线观看| 国产日韩欧美在线精品| av在线观看视频网站免费| 久久99热这里只频精品6学生 | 亚洲av中文av极速乱| 青青草视频在线视频观看| 久久草成人影院| 国产精品一区二区在线观看99 | 国产视频首页在线观看| 免费黄网站久久成人精品| 三级男女做爰猛烈吃奶摸视频| 国产精品一及| 亚洲欧美一区二区三区国产| 少妇高潮的动态图| 国产男人的电影天堂91| 99热这里只有是精品在线观看| 国产伦在线观看视频一区| 男女下面进入的视频免费午夜| 97热精品久久久久久| 久久久久久久亚洲中文字幕| 亚洲自拍偷在线| 色播亚洲综合网| 亚洲天堂国产精品一区在线| 两个人的视频大全免费| 嫩草影院精品99| 国产精品久久久久久久久免| 亚洲国产高清在线一区二区三| 少妇的逼好多水| 少妇熟女欧美另类| 国产精品国产高清国产av| 久久久久精品久久久久真实原创| 国产精品女同一区二区软件| 亚洲最大成人手机在线| 国产精品蜜桃在线观看| 观看免费一级毛片| 不卡视频在线观看欧美| 插逼视频在线观看| www日本黄色视频网| 免费电影在线观看免费观看| 我的女老师完整版在线观看| 国产黄片视频在线免费观看| 久久久久性生活片| 欧美成人免费av一区二区三区| 欧美潮喷喷水| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久久电影| 国产女主播在线喷水免费视频网站 | 亚洲国产日韩欧美精品在线观看| 欧美日韩综合久久久久久| 中文欧美无线码| 中文精品一卡2卡3卡4更新| 卡戴珊不雅视频在线播放| 在线免费观看不下载黄p国产| 男女那种视频在线观看| 波野结衣二区三区在线| 自拍偷自拍亚洲精品老妇| 变态另类丝袜制服| 91久久精品国产一区二区成人| 久久国内精品自在自线图片| 国产精品伦人一区二区| 成人美女网站在线观看视频| 插逼视频在线观看| 久久亚洲精品不卡| 乱人视频在线观看| 久99久视频精品免费| 亚洲在线自拍视频| 一个人看的www免费观看视频| 国产极品天堂在线| 美女黄网站色视频| 天堂√8在线中文| 精品欧美国产一区二区三| 老司机福利观看| 亚洲性久久影院| 国产精品一区二区在线观看99 | 一个人看视频在线观看www免费| 亚洲精品自拍成人| 九九热线精品视视频播放| 欧美丝袜亚洲另类| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 一级二级三级毛片免费看| 亚洲美女搞黄在线观看| 中文天堂在线官网| 欧美激情在线99| 国产单亲对白刺激| kizo精华| 人妻制服诱惑在线中文字幕| 亚洲经典国产精华液单| 国产午夜精品久久久久久一区二区三区| 高清日韩中文字幕在线| 性插视频无遮挡在线免费观看| 亚洲av电影不卡..在线观看| 国产午夜福利久久久久久| 97在线视频观看| 国产91av在线免费观看| 日本黄色视频三级网站网址| 91久久精品国产一区二区成人| 啦啦啦韩国在线观看视频| 精品人妻视频免费看| 午夜免费激情av| 日韩在线高清观看一区二区三区| 精品久久久久久电影网 | eeuss影院久久| 国产精品嫩草影院av在线观看| 内射极品少妇av片p| 免费不卡的大黄色大毛片视频在线观看 | 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 中文亚洲av片在线观看爽| 精品久久久久久久久亚洲| 欧美一区二区国产精品久久精品| 久99久视频精品免费| 久久人妻av系列| 一个人看的www免费观看视频| 中文精品一卡2卡3卡4更新| 国产v大片淫在线免费观看| av在线蜜桃| 国产黄片视频在线免费观看| 久久久久久久久中文| 婷婷六月久久综合丁香| 如何舔出高潮| 2021天堂中文幕一二区在线观| 男人的好看免费观看在线视频| 中文字幕久久专区| 国产免费一级a男人的天堂| 九九久久精品国产亚洲av麻豆| 嫩草影院新地址| 国产美女午夜福利| 欧美变态另类bdsm刘玥| 汤姆久久久久久久影院中文字幕 | 91久久精品国产一区二区成人| av专区在线播放| 免费播放大片免费观看视频在线观看 | 免费观看精品视频网站| 亚洲av.av天堂| 超碰97精品在线观看| 国产免费福利视频在线观看| 好男人视频免费观看在线| 中国美白少妇内射xxxbb| 在线免费观看的www视频| 亚洲精品456在线播放app| 91久久精品电影网| 亚洲va在线va天堂va国产| 久久久久免费精品人妻一区二区| 一级毛片电影观看 | 亚洲精品成人久久久久久| 精品一区二区三区视频在线| 国产精品.久久久| 九色成人免费人妻av| 51国产日韩欧美| 欧美一级a爱片免费观看看| 亚洲成色77777| 自拍偷自拍亚洲精品老妇| 一边亲一边摸免费视频| 非洲黑人性xxxx精品又粗又长| 久久久久国产网址| 我要看日韩黄色一级片| 啦啦啦观看免费观看视频高清| 男人舔女人下体高潮全视频| 在线免费观看不下载黄p国产| 国产伦理片在线播放av一区| av国产免费在线观看| 免费观看的影片在线观看| 午夜免费男女啪啪视频观看| 黄片wwwwww| 六月丁香七月| 波多野结衣巨乳人妻| www.色视频.com| 青春草国产在线视频| 国产伦理片在线播放av一区| 日韩大片免费观看网站 | 黄色配什么色好看| 久久久久久伊人网av| 哪个播放器可以免费观看大片| 99久久精品国产国产毛片| av在线亚洲专区| 免费av毛片视频| 中文乱码字字幕精品一区二区三区 | 欧美3d第一页| 国产黄片美女视频| av又黄又爽大尺度在线免费看 | 日韩国内少妇激情av| 最近最新中文字幕免费大全7| 亚洲四区av| 国产色婷婷99| 日本黄大片高清| 99在线人妻在线中文字幕| 三级毛片av免费| 国产精品电影一区二区三区| 婷婷色综合大香蕉| 一区二区三区乱码不卡18| 国产成人freesex在线| 丰满乱子伦码专区| 少妇的逼水好多| 国产国拍精品亚洲av在线观看| 亚洲精品一区蜜桃| 久久这里只有精品中国| 欧美另类亚洲清纯唯美| 五月伊人婷婷丁香| 亚洲国产成人一精品久久久| 日本午夜av视频| 国产亚洲一区二区精品| 亚洲欧美精品专区久久| 干丝袜人妻中文字幕| 久久精品人妻少妇| 51国产日韩欧美| 国产女主播在线喷水免费视频网站 | 99在线视频只有这里精品首页| 最近中文字幕2019免费版| 99国产精品一区二区蜜桃av| 国产精品乱码一区二三区的特点| 91在线精品国自产拍蜜月| 国产亚洲精品av在线| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 少妇的逼好多水| 国产一区有黄有色的免费视频 | 久久久精品94久久精品| 波多野结衣巨乳人妻| 99在线人妻在线中文字幕| 国产一区有黄有色的免费视频 | 国产一区有黄有色的免费视频 | 久久久久久大精品| 午夜视频国产福利| 亚洲国产欧美人成| 观看免费一级毛片| 国产一区亚洲一区在线观看| 特大巨黑吊av在线直播| 午夜a级毛片| 老司机福利观看| 久久国产乱子免费精品| 18禁动态无遮挡网站| 久久99蜜桃精品久久| 99热全是精品| av播播在线观看一区| 日韩一区二区三区影片| 亚洲成人精品中文字幕电影| 免费看日本二区| 亚洲欧美日韩无卡精品| 久久草成人影院| 美女大奶头视频| 成年免费大片在线观看| av在线观看视频网站免费| 非洲黑人性xxxx精品又粗又长| 久久久亚洲精品成人影院| 欧美色视频一区免费| 国产极品精品免费视频能看的| 成人无遮挡网站| 三级经典国产精品| 日韩欧美精品免费久久| 日本黄色视频三级网站网址| 1024手机看黄色片| 中文天堂在线官网| 少妇熟女欧美另类| 日韩欧美精品v在线| 女人被狂操c到高潮| 国产精品福利在线免费观看| 久热久热在线精品观看| 免费播放大片免费观看视频在线观看 | 免费看光身美女| 午夜亚洲福利在线播放| 天天一区二区日本电影三级| 日本爱情动作片www.在线观看| 亚洲在久久综合| 日本av手机在线免费观看| 特级一级黄色大片| 日本av手机在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩国内少妇激情av| 免费观看性生交大片5| 日本熟妇午夜| 18禁在线无遮挡免费观看视频| 亚洲成色77777| 好男人视频免费观看在线| 一级黄色大片毛片| 能在线免费看毛片的网站| 熟女电影av网| 在线免费十八禁| 亚洲中文字幕一区二区三区有码在线看| 大香蕉97超碰在线| 国产女主播在线喷水免费视频网站 | 尤物成人国产欧美一区二区三区| 成年免费大片在线观看| 在线观看66精品国产| 91久久精品国产一区二区三区| 听说在线观看完整版免费高清| 国产精品久久电影中文字幕| 日本与韩国留学比较| 久久久国产成人免费| 久久精品91蜜桃| 国产黄a三级三级三级人| 国产综合懂色| 国产精品av视频在线免费观看| 成人性生交大片免费视频hd| 看免费成人av毛片| 一级毛片电影观看 | 99久久九九国产精品国产免费| 丰满乱子伦码专区| 亚洲精品一区蜜桃| 日本免费一区二区三区高清不卡| 亚洲va在线va天堂va国产| av福利片在线观看| 91精品伊人久久大香线蕉| 国产毛片a区久久久久| 我要搜黄色片| 欧美日本视频| 欧美一区二区国产精品久久精品| 啦啦啦观看免费观看视频高清| av免费在线看不卡| 免费黄网站久久成人精品| 3wmmmm亚洲av在线观看| 国产精品一二三区在线看| 男女啪啪激烈高潮av片| 成人性生交大片免费视频hd| 精品熟女少妇av免费看| 性色avwww在线观看| av天堂中文字幕网| 一二三四中文在线观看免费高清| 特级一级黄色大片| 在线观看美女被高潮喷水网站| 亚洲成av人片在线播放无| 精品久久国产蜜桃| 别揉我奶头 嗯啊视频| 久久99精品国语久久久| 在现免费观看毛片| 在线观看66精品国产| 在线观看美女被高潮喷水网站| 国产精品综合久久久久久久免费| 偷拍熟女少妇极品色| 国产午夜精品一二区理论片| 五月玫瑰六月丁香| 最近最新中文字幕大全电影3| 99久久精品热视频| 精品99又大又爽又粗少妇毛片| 久久精品国产自在天天线| 婷婷色av中文字幕| 国产乱人视频| 久久久色成人| 在线a可以看的网站| 免费观看性生交大片5| 欧美日本视频| 久久精品影院6| 久久人妻av系列| 99热网站在线观看| 国产精品伦人一区二区| 亚洲精品日韩在线中文字幕| 日韩欧美精品免费久久| 女人久久www免费人成看片 | 久久久精品大字幕| 97人妻精品一区二区三区麻豆| 岛国在线免费视频观看| 99热精品在线国产| 久久久久久久久大av| 久99久视频精品免费| www.av在线官网国产| 成人漫画全彩无遮挡| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 婷婷色综合大香蕉| 国产中年淑女户外野战色| 99视频精品全部免费 在线| 亚洲国产精品成人综合色| 尾随美女入室| 亚洲欧洲日产国产| 日韩视频在线欧美| 国产精品久久久久久精品电影| 国产精品久久久久久久久免| 国产三级中文精品| 午夜精品在线福利| 中文亚洲av片在线观看爽| 欧美zozozo另类| 欧美人与善性xxx| 我要看日韩黄色一级片| 午夜视频国产福利| 高清av免费在线| 夜夜爽夜夜爽视频| 美女黄网站色视频| 日本三级黄在线观看| 亚洲欧美精品专区久久| 欧美xxxx性猛交bbbb| 国产女主播在线喷水免费视频网站 | 少妇被粗大猛烈的视频| 国产乱来视频区| 亚洲成色77777| 亚洲五月天丁香| 欧美+日韩+精品| 99九九线精品视频在线观看视频| 精品久久久噜噜| 久久久成人免费电影| 高清毛片免费看| 日韩精品有码人妻一区| 国产69精品久久久久777片| 男人舔奶头视频| 男人和女人高潮做爰伦理| 麻豆久久精品国产亚洲av| 久99久视频精品免费| 国产爱豆传媒在线观看| 51国产日韩欧美| 国产av在哪里看| 国产精品伦人一区二区| 男人舔女人下体高潮全视频| 亚洲国产精品久久男人天堂| 国产精品,欧美在线| 成人综合一区亚洲| a级毛片免费高清观看在线播放| 国产精品人妻久久久久久| 18禁在线播放成人免费| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看性生交大片5| 人妻制服诱惑在线中文字幕| av在线播放精品| 成年版毛片免费区| 美女内射精品一级片tv| 在线a可以看的网站| 一级毛片我不卡| 国产av码专区亚洲av| 免费观看的影片在线观看| av在线蜜桃| 国产黄色视频一区二区在线观看 | 99久国产av精品国产电影| 99视频精品全部免费 在线| 亚洲国产最新在线播放| 欧美日韩综合久久久久久| 精品久久久久久久末码| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 伊人久久精品亚洲午夜| 综合色丁香网| 国产伦一二天堂av在线观看| 综合色丁香网| 日韩在线高清观看一区二区三区| 一级毛片久久久久久久久女| 黄色配什么色好看| 欧美xxxx黑人xx丫x性爽| 色网站视频免费| 亚洲精品国产av成人精品| 日本免费在线观看一区| 久久久色成人| 国产精品一区二区性色av| 99在线视频只有这里精品首页| 纵有疾风起免费观看全集完整版 | 精品久久国产蜜桃| 麻豆一二三区av精品| 丰满乱子伦码专区| 久久亚洲国产成人精品v| 青青草视频在线视频观看| 亚洲欧美日韩东京热| 亚洲高清免费不卡视频| 三级男女做爰猛烈吃奶摸视频| 日韩高清综合在线| 国产午夜福利久久久久久| 国产精品嫩草影院av在线观看| 波野结衣二区三区在线| 一本久久精品| 日本黄色视频三级网站网址| 中文乱码字字幕精品一区二区三区 | 天堂网av新在线| 麻豆国产97在线/欧美| 高清视频免费观看一区二区 | 日本一二三区视频观看| 国产乱人视频| 国产一级毛片在线| 色噜噜av男人的天堂激情| 级片在线观看| 欧美性猛交黑人性爽| 国内精品一区二区在线观看| 午夜福利视频1000在线观看| 久久精品人妻少妇| 国产精品久久久久久久久免| 99热网站在线观看| 可以在线观看毛片的网站| 国产国拍精品亚洲av在线观看| 成人美女网站在线观看视频| 欧美性猛交黑人性爽| 欧美人与善性xxx| 国模一区二区三区四区视频| 国产精品一区二区性色av| 精品一区二区三区视频在线| 精品一区二区三区人妻视频| 日韩亚洲欧美综合| 国产精品一区二区性色av| 国内精品宾馆在线| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| 亚洲中文字幕日韩| 午夜福利视频1000在线观看| 欧美三级亚洲精品| 一区二区三区乱码不卡18| 熟女电影av网| 亚洲av免费高清在线观看| 久久热精品热| 国产中年淑女户外野战色| 亚洲av中文av极速乱| 日韩一区二区三区影片| 中文欧美无线码| 欧美成人精品欧美一级黄| 少妇被粗大猛烈的视频| 免费黄网站久久成人精品| 搡女人真爽免费视频火全软件| 观看免费一级毛片| 国产精品爽爽va在线观看网站| 亚洲欧美日韩东京热| 能在线免费看毛片的网站| 精品一区二区三区人妻视频|