• <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ì)算法求體積本征值
    伦理电影大哥的女人| 日韩一本色道免费dvd| 日本在线视频免费播放| 夜夜爽天天搞| 精品人妻视频免费看| 久久久久国产精品人妻aⅴ院| 免费黄网站久久成人精品| 日本爱情动作片www.在线观看 | 日本一二三区视频观看| 国产精品久久久久久精品电影| 在线观看免费视频日本深夜| 亚洲av成人精品一区久久| 麻豆成人午夜福利视频| 丰满乱子伦码专区| 亚洲成人精品中文字幕电影| 国产亚洲精品综合一区在线观看| 国产69精品久久久久777片| 一个人免费在线观看电影| 精品久久久久久久久亚洲| 日本成人三级电影网站| 日本一本二区三区精品| 国产精品嫩草影院av在线观看| 色哟哟·www| 村上凉子中文字幕在线| 欧美不卡视频在线免费观看| 天堂影院成人在线观看| 欧美色视频一区免费| 国产精品精品国产色婷婷| 欧美精品国产亚洲| 亚洲一区二区三区色噜噜| 午夜视频国产福利| 成人亚洲精品av一区二区| 国产av一区在线观看免费| 看片在线看免费视频| 干丝袜人妻中文字幕| 1024手机看黄色片| av在线播放精品| АⅤ资源中文在线天堂| 91在线精品国自产拍蜜月| 在现免费观看毛片| 男人狂女人下面高潮的视频| 免费人成在线观看视频色| 欧美色欧美亚洲另类二区| 韩国av在线不卡| 在线观看午夜福利视频| 亚洲内射少妇av| 18禁黄网站禁片免费观看直播| 国产免费一级a男人的天堂| 精品福利观看| 尾随美女入室| 观看免费一级毛片| 国产精品伦人一区二区| 99久国产av精品国产电影| av卡一久久| 我要看日韩黄色一级片| 亚洲av免费在线观看| 日本色播在线视频| 精品一区二区三区人妻视频| 人妻制服诱惑在线中文字幕| 中文字幕久久专区| 高清毛片免费看| 成年免费大片在线观看| 欧美色视频一区免费| 18禁在线无遮挡免费观看视频 | 亚洲丝袜综合中文字幕| 天堂影院成人在线观看| 亚洲av不卡在线观看| 能在线免费观看的黄片| 欧美成人a在线观看| 免费大片18禁| 亚洲内射少妇av| 99热只有精品国产| 日韩精品有码人妻一区| 亚洲精品一区av在线观看| 永久网站在线| 一级黄色大片毛片| 又黄又爽又免费观看的视频| 久久久欧美国产精品| 国内久久婷婷六月综合欲色啪| 国产亚洲欧美98| 成人av在线播放网站| 毛片一级片免费看久久久久| 黄色视频,在线免费观看| 国产一区亚洲一区在线观看| 久久99热这里只有精品18| 3wmmmm亚洲av在线观看| 最好的美女福利视频网| 欧美三级亚洲精品| 伊人久久精品亚洲午夜| 搡老熟女国产l中国老女人| 中文字幕久久专区| 最近最新中文字幕大全电影3| 国产精品久久电影中文字幕| 国产单亲对白刺激| 久久久久性生活片| 午夜爱爱视频在线播放| 久久九九热精品免费| 国产麻豆成人av免费视频| 女生性感内裤真人,穿戴方法视频| 国产三级中文精品| 日韩一区二区视频免费看| 男女那种视频在线观看| 成年av动漫网址| 校园春色视频在线观看| 国产精品伦人一区二区| 亚洲一区高清亚洲精品| 久久99热6这里只有精品| 少妇熟女欧美另类| 51国产日韩欧美| 99在线视频只有这里精品首页| 一边摸一边抽搐一进一小说| 亚洲va在线va天堂va国产| 国内精品宾馆在线| 欧美成人免费av一区二区三区| 日本黄色片子视频| 桃色一区二区三区在线观看| 美女 人体艺术 gogo| 最近的中文字幕免费完整| 乱系列少妇在线播放| 99久国产av精品国产电影| 国产精华一区二区三区| 亚洲五月天丁香| 成人高潮视频无遮挡免费网站| 色尼玛亚洲综合影院| 欧美极品一区二区三区四区| 国产v大片淫在线免费观看| 亚洲精品456在线播放app| 女人被狂操c到高潮| 秋霞在线观看毛片| aaaaa片日本免费| 夜夜夜夜夜久久久久| 蜜桃亚洲精品一区二区三区| 久久精品国产99精品国产亚洲性色| 日韩 亚洲 欧美在线| 乱系列少妇在线播放| 亚洲成人久久爱视频| 欧美xxxx性猛交bbbb| 国产中年淑女户外野战色| 久久精品国产亚洲网站| 桃色一区二区三区在线观看| 成人二区视频| 免费看日本二区| 午夜福利成人在线免费观看| 成年av动漫网址| 国产久久久一区二区三区| 欧美高清性xxxxhd video| 亚洲av电影不卡..在线观看| 亚洲精品乱码久久久v下载方式| 色尼玛亚洲综合影院| 国产 一区 欧美 日韩| 亚洲在线观看片| 国产伦精品一区二区三区四那| 日韩一区二区视频免费看| 人人妻人人看人人澡| 国产亚洲精品综合一区在线观看| 欧美日本视频| 久久久久久久久久成人| 色综合色国产| 国产精品久久久久久久电影| 日本黄大片高清| 午夜福利高清视频| 中文字幕av成人在线电影| av天堂在线播放| 欧美日韩在线观看h| 欧美日韩精品成人综合77777| 久久久成人免费电影| а√天堂www在线а√下载| 日本免费a在线| 国产成人福利小说| 久久久久久久久久黄片| 亚州av有码| 长腿黑丝高跟| 女生性感内裤真人,穿戴方法视频| 亚洲熟妇熟女久久| 精品久久国产蜜桃| 成人特级av手机在线观看| 色综合站精品国产| 97在线视频观看| 我要看日韩黄色一级片| 俄罗斯特黄特色一大片| 校园人妻丝袜中文字幕| 最近中文字幕高清免费大全6| 成人漫画全彩无遮挡| 一个人免费在线观看电影| 久久6这里有精品| 亚洲高清免费不卡视频| 最后的刺客免费高清国语| 国产亚洲av嫩草精品影院| 亚洲va在线va天堂va国产| 久久午夜亚洲精品久久| 麻豆成人午夜福利视频| 一进一出抽搐动态| 在线免费观看不下载黄p国产| 非洲黑人性xxxx精品又粗又长| 久久人人爽人人片av| 久久久国产成人精品二区| 午夜福利成人在线免费观看| 亚洲国产精品久久男人天堂| 最近最新中文字幕大全电影3| 国产探花在线观看一区二区| 久久久久久久久中文| 亚洲,欧美,日韩| 久久久久国产精品人妻aⅴ院| 国产黄片美女视频| 在线天堂最新版资源| 99视频精品全部免费 在线| 悠悠久久av| av在线观看视频网站免费| 精华霜和精华液先用哪个| 菩萨蛮人人尽说江南好唐韦庄 | 日韩强制内射视频| 欧美日韩乱码在线| 国产精品久久视频播放| 俺也久久电影网| 亚洲人与动物交配视频| 99久国产av精品| 中文字幕人妻熟人妻熟丝袜美| 黄色欧美视频在线观看| 亚洲国产日韩欧美精品在线观看| 国产精品乱码一区二三区的特点| 国产熟女欧美一区二区| 一区福利在线观看| 热99在线观看视频| 深夜精品福利| 久久综合国产亚洲精品| 久久久久国内视频| 在线免费观看不下载黄p国产| 国产69精品久久久久777片| 亚洲成av人片在线播放无| 国内少妇人妻偷人精品xxx网站| 亚洲精品色激情综合| 最近中文字幕高清免费大全6| 最近2019中文字幕mv第一页| 老女人水多毛片| av在线亚洲专区| 99久久久亚洲精品蜜臀av| 亚洲国产欧洲综合997久久,| 一区二区三区免费毛片| 国产日本99.免费观看| 看黄色毛片网站| 悠悠久久av| 久久99热这里只有精品18| 久久久精品94久久精品| 一夜夜www| 99热6这里只有精品| 成人综合一区亚洲| 搞女人的毛片| 嫩草影视91久久| 精品乱码久久久久久99久播| 国产精品爽爽va在线观看网站| 亚洲精品久久国产高清桃花| 老司机福利观看| 免费看光身美女| 日本与韩国留学比较| 一级毛片电影观看 | 欧美成人精品欧美一级黄| 91在线观看av| 午夜福利在线观看免费完整高清在 | 亚洲av免费高清在线观看| 免费看日本二区| 亚洲,欧美,日韩| 午夜福利高清视频| 日韩欧美一区二区三区在线观看| 乱人视频在线观看| 国产黄色视频一区二区在线观看 | 国产精品久久久久久精品电影| 亚洲性夜色夜夜综合| 非洲黑人性xxxx精品又粗又长| a级一级毛片免费在线观看| 精品一区二区三区视频在线观看免费| 日本色播在线视频| 国产精品日韩av在线免费观看| 亚洲欧美清纯卡通| 国产久久久一区二区三区| 国产亚洲精品综合一区在线观看| 身体一侧抽搐| 一进一出好大好爽视频| 丰满乱子伦码专区| 亚洲最大成人手机在线| 亚洲美女搞黄在线观看 | av卡一久久| 国产69精品久久久久777片| 又黄又爽又免费观看的视频| 99久久中文字幕三级久久日本| 午夜爱爱视频在线播放| 波多野结衣巨乳人妻| 2021天堂中文幕一二区在线观| 免费人成视频x8x8入口观看| 在线免费观看不下载黄p国产| 午夜久久久久精精品| 黄色配什么色好看| 永久网站在线| 一区福利在线观看| 十八禁网站免费在线| 一区二区三区四区激情视频 | 九色成人免费人妻av| 久久久国产成人精品二区| 尾随美女入室| 日韩中字成人| 春色校园在线视频观看| 国产欧美日韩精品一区二区| 少妇人妻一区二区三区视频| 亚洲av第一区精品v没综合| 国产av一区在线观看免费| 夜夜看夜夜爽夜夜摸| 女生性感内裤真人,穿戴方法视频| 一区二区三区免费毛片| 秋霞在线观看毛片| 国产精品一区二区免费欧美| 伦精品一区二区三区| 欧美一区二区精品小视频在线| 最新中文字幕久久久久| 色综合站精品国产| 性色avwww在线观看| 啦啦啦啦在线视频资源| 亚洲人成网站在线观看播放| 亚洲国产高清在线一区二区三| 色综合站精品国产| 午夜精品在线福利| 尤物成人国产欧美一区二区三区| 亚洲无线观看免费| 国产精品久久久久久亚洲av鲁大| 哪里可以看免费的av片| 99久久精品国产国产毛片| 国产精品久久电影中文字幕| 欧美激情国产日韩精品一区| 嫩草影视91久久| 色哟哟·www| 午夜激情福利司机影院| 午夜精品一区二区三区免费看| 身体一侧抽搐| 亚洲国产欧洲综合997久久,| 日韩欧美精品免费久久| 久久精品国产鲁丝片午夜精品| 亚洲av五月六月丁香网| 日韩中字成人| 99国产精品一区二区蜜桃av| 精品午夜福利视频在线观看一区| ponron亚洲| 免费电影在线观看免费观看| 男女视频在线观看网站免费| 老司机午夜福利在线观看视频| 国产精品不卡视频一区二区| 日韩人妻高清精品专区| 91麻豆精品激情在线观看国产| 夜夜爽天天搞| 校园春色视频在线观看| 搡老熟女国产l中国老女人| 婷婷精品国产亚洲av在线| 国产一区二区激情短视频| 亚洲av免费高清在线观看| 淫妇啪啪啪对白视频| 日本-黄色视频高清免费观看| 黑人高潮一二区| 久久九九热精品免费| 亚洲欧美日韩高清在线视频| 人妻夜夜爽99麻豆av| 亚洲内射少妇av| 男人和女人高潮做爰伦理| 大又大粗又爽又黄少妇毛片口| av.在线天堂| 亚洲国产欧洲综合997久久,| 日韩高清综合在线| 国产aⅴ精品一区二区三区波| 在线天堂最新版资源| 丰满人妻一区二区三区视频av| 一本久久中文字幕| 成人美女网站在线观看视频| 国产成人精品久久久久久| 亚洲熟妇熟女久久| 在线免费十八禁| a级一级毛片免费在线观看| 此物有八面人人有两片| 我要看日韩黄色一级片| 干丝袜人妻中文字幕| 人妻制服诱惑在线中文字幕| 男人舔女人下体高潮全视频| 两个人视频免费观看高清| 九九爱精品视频在线观看| 在线观看66精品国产| 久久中文看片网| 国产欧美日韩一区二区精品| 99精品在免费线老司机午夜| 久久热精品热| or卡值多少钱| 干丝袜人妻中文字幕| 亚洲国产精品合色在线| 婷婷六月久久综合丁香| 亚洲精品国产av成人精品 | 热99在线观看视频| 性色avwww在线观看| 亚洲三级黄色毛片| 亚洲精品亚洲一区二区| 欧美一区二区亚洲| 精品99又大又爽又粗少妇毛片| 成人亚洲欧美一区二区av| 我的女老师完整版在线观看| 久久午夜亚洲精品久久| 小蜜桃在线观看免费完整版高清| 久久人人精品亚洲av| 欧美在线一区亚洲| 最近最新中文字幕大全电影3| 免费在线观看影片大全网站| 狂野欧美激情性xxxx在线观看| 亚洲精品456在线播放app| 变态另类成人亚洲欧美熟女| 十八禁网站免费在线| 秋霞在线观看毛片| 内射极品少妇av片p| 男女做爰动态图高潮gif福利片| 高清毛片免费观看视频网站| 国产精品电影一区二区三区| 自拍偷自拍亚洲精品老妇| 少妇猛男粗大的猛烈进出视频 | 国产精品一区二区三区四区久久| 欧美xxxx性猛交bbbb| 亚洲最大成人中文| 国产精品一及| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费| 久久婷婷人人爽人人干人人爱| 亚洲av免费高清在线观看| 97超碰精品成人国产| 日韩精品有码人妻一区| 国产精品嫩草影院av在线观看| 国产亚洲精品综合一区在线观看| eeuss影院久久| 国产精华一区二区三区| 18+在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美日韩一区二区精品| 成人永久免费在线观看视频| 国产精品久久久久久精品电影| 麻豆成人午夜福利视频| 欧美日韩国产亚洲二区| 日韩欧美一区二区三区在线观看| 午夜爱爱视频在线播放| 免费看光身美女| 国产精品人妻久久久影院| 桃色一区二区三区在线观看| 国产精品亚洲一级av第二区| 成人特级av手机在线观看| 精品99又大又爽又粗少妇毛片| av在线蜜桃| а√天堂www在线а√下载| 亚洲在线自拍视频| 淫秽高清视频在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲精品一卡2卡三卡4卡5卡| 亚洲无线在线观看| 99热这里只有精品一区| 亚洲成人久久性| 岛国在线免费视频观看| a级毛片免费高清观看在线播放| 如何舔出高潮| 国产免费一级a男人的天堂| 我要搜黄色片| 久久久精品大字幕| 日本色播在线视频| 男人和女人高潮做爰伦理| 亚洲,欧美,日韩| 亚洲自偷自拍三级| 91精品国产九色| 两个人的视频大全免费| 成人午夜高清在线视频| 精品国内亚洲2022精品成人| 亚洲精品成人久久久久久| 国产片特级美女逼逼视频| 成人高潮视频无遮挡免费网站| 亚洲av成人av| 成人亚洲精品av一区二区| 国产高清视频在线播放一区| 丝袜美腿在线中文| 国产激情偷乱视频一区二区| 男女做爰动态图高潮gif福利片| 51国产日韩欧美| 99视频精品全部免费 在线| 欧美性猛交╳xxx乱大交人| 亚洲天堂国产精品一区在线| 日本与韩国留学比较| 午夜精品在线福利| 亚洲国产精品合色在线| 精品久久国产蜜桃| 亚洲第一电影网av| 亚洲乱码一区二区免费版| 99久国产av精品国产电影| 国产精品一区二区三区四区免费观看 | 舔av片在线| 久久人妻av系列| 3wmmmm亚洲av在线观看| 在线免费观看不下载黄p国产| 欧美区成人在线视频| 国产黄色小视频在线观看| 啦啦啦观看免费观看视频高清| 日韩精品有码人妻一区| 免费在线观看成人毛片| 女同久久另类99精品国产91| 俄罗斯特黄特色一大片| 久久久精品大字幕| 一a级毛片在线观看| 午夜激情欧美在线| 熟女电影av网| 日本a在线网址| 可以在线观看毛片的网站| 久久精品国产亚洲av香蕉五月| aaaaa片日本免费| 亚洲第一区二区三区不卡| 色播亚洲综合网| 国产伦一二天堂av在线观看| 日韩欧美三级三区| 一卡2卡三卡四卡精品乱码亚洲| 在线观看av片永久免费下载| 性色avwww在线观看| 久久久久国产网址| a级毛片a级免费在线| 久久久久久久亚洲中文字幕| 日韩中字成人| 麻豆精品久久久久久蜜桃| 成人av一区二区三区在线看| 午夜福利在线观看吧| 天天一区二区日本电影三级| 亚洲精品国产成人久久av| 九九在线视频观看精品| 岛国在线免费视频观看| 国产亚洲精品av在线| 国产国拍精品亚洲av在线观看| 精品99又大又爽又粗少妇毛片| 婷婷精品国产亚洲av在线| 给我免费播放毛片高清在线观看| 精品福利观看| а√天堂www在线а√下载| 久久亚洲国产成人精品v| 精品久久久久久久久av| 国产色爽女视频免费观看| 久久精品国产亚洲av香蕉五月| 一级av片app| 国产av在哪里看| 国产美女午夜福利| 麻豆久久精品国产亚洲av| 午夜免费激情av| 变态另类成人亚洲欧美熟女| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站在线观看播放| 深夜精品福利| 午夜久久久久精精品| 久久这里只有精品中国| 老司机午夜福利在线观看视频| 亚洲真实伦在线观看| 51国产日韩欧美| 一进一出抽搐gif免费好疼| 亚洲最大成人手机在线| 精品不卡国产一区二区三区| 亚洲精品影视一区二区三区av| 黄色日韩在线| 精品欧美国产一区二区三| 亚洲精品日韩在线中文字幕 | 1024手机看黄色片| 精品熟女少妇av免费看| 亚洲欧美清纯卡通| 午夜福利在线观看免费完整高清在 | 亚洲av中文av极速乱| 人妻少妇偷人精品九色| 国产高潮美女av| 亚洲最大成人手机在线| 中文亚洲av片在线观看爽| 禁无遮挡网站| 在线观看av片永久免费下载| 桃色一区二区三区在线观看| 人人妻人人澡欧美一区二区| 欧美xxxx性猛交bbbb| 国内少妇人妻偷人精品xxx网站| 国产精品嫩草影院av在线观看| 不卡一级毛片| 免费观看的影片在线观看| 熟女电影av网| 国内精品宾馆在线| 国产视频一区二区在线看| 久久午夜福利片| 亚洲成人av在线免费| 成人无遮挡网站| 不卡一级毛片| 可以在线观看的亚洲视频| 12—13女人毛片做爰片一| 能在线免费观看的黄片| 久久久午夜欧美精品| 久久人人爽人人片av| 国产精品久久久久久亚洲av鲁大| 秋霞在线观看毛片| 国产精品福利在线免费观看| 在线看三级毛片| 国产高清视频在线播放一区| 最近手机中文字幕大全| .国产精品久久| 丝袜美腿在线中文| 一个人看视频在线观看www免费| 激情 狠狠 欧美| 午夜视频国产福利| 色哟哟·www| 在线国产一区二区在线| 久久久久久久午夜电影| 大又大粗又爽又黄少妇毛片口| 免费不卡的大黄色大毛片视频在线观看 | 最好的美女福利视频网| 秋霞在线观看毛片| 在线观看av片永久免费下载| 国产片特级美女逼逼视频| 天天躁夜夜躁狠狠久久av| 精品无人区乱码1区二区| 日本与韩国留学比较| 欧美极品一区二区三区四区| 国产成人91sexporn| 成人综合一区亚洲| 级片在线观看| 国产亚洲91精品色在线| 赤兔流量卡办理|