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

    微波真實(shí)光譜響應(yīng)對(duì)快速輻射模擬影響分析

    2023-10-09 10:23:06谷松巖康志明
    光譜學(xué)與光譜分析 2023年10期
    關(guān)鍵詞:響應(yīng)函數(shù)天氣預(yù)報(bào)光學(xué)

    陳 昊,王 皓,韓 威,谷松巖,張 鵬,康志明

    1. 江蘇省氣象臺(tái),江蘇 南京 210008 2. 中國(guó)氣象局交通氣象重點(diǎn)開(kāi)放實(shí)驗(yàn)室,江蘇 南京 210008 3. 中國(guó)氣象局地球系統(tǒng)數(shù)值預(yù)報(bào)中心,北京 100081 4. 國(guó)家衛(wèi)星氣象中心,北京 100081

    引 言

    微波在云雨等復(fù)雜氣象條件下具有一定的穿透性,因而星載微波遙感器具有全天候多時(shí)相工作的特點(diǎn),其被廣泛應(yīng)用于災(zāi)害性天氣預(yù)報(bào)、防災(zāi)減災(zāi)監(jiān)測(cè)等領(lǐng)域,同時(shí)也在全球數(shù)值天氣預(yù)報(bào)中發(fā)揮著重要作用。

    風(fēng)云三號(hào)D星(FengYun-3D,FY-3D)是我國(guó)最新一顆極軌氣象衛(wèi)星,其上搭載有諸多微波遙感器,包含微波成像儀(MicroWave Radiation Imager,MWRI)、微波溫度探測(cè)儀(MicroWave Temperature Sounder-2,MWTS-2)、微波濕度探測(cè)儀(MicroWave Humidity Sounder-2,MWHS-2)等[1]。其中MWTS-2主要用于對(duì)大氣溫度狀態(tài)的三維探測(cè),其觀測(cè)資料可被應(yīng)用于短時(shí)臨近災(zāi)害性天氣監(jiān)測(cè)與預(yù)報(bào)及全球/區(qū)域數(shù)值天氣預(yù)報(bào)系統(tǒng)。

    對(duì)于全球數(shù)值天氣預(yù)報(bào)系統(tǒng)而言,為得到更加精確的初始場(chǎng),常常需要同化大量觀測(cè)資料。在諸多觀測(cè)資料中,衛(wèi)星遙感觀測(cè)資料由于其具有覆蓋面廣、觀測(cè)頻次高、觀測(cè)要素豐富等特點(diǎn),在全球數(shù)值天氣預(yù)報(bào)模式中發(fā)揮著重要作用,衛(wèi)星遙感觀測(cè)資料在眾多被數(shù)值天氣預(yù)報(bào)系統(tǒng)同化的觀測(cè)資料中占比也最高。而在眾多衛(wèi)星遙感資料中,微波溫度計(jì)對(duì)于全球數(shù)值天氣預(yù)報(bào)準(zhǔn)確率貢獻(xiàn)最大[2]。

    在全球數(shù)值天氣預(yù)報(bào)系統(tǒng)中同化微波溫度計(jì)資料,一般直接同化微波溫度計(jì)觀測(cè)到的微波輻射。在同化流程中,核心是快速準(zhǔn)確地模擬微波溫度計(jì)觀測(cè)到的微波輻射。目前,國(guó)際上主流的快速輻射傳輸模型主要有RTTOV(radiative transfer for TOVS)、CRTM(community radiative transfer model)等[3]。在主流的全球數(shù)值天氣預(yù)報(bào)系統(tǒng)中,歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts,ECMWF)研發(fā)的ECMWF-IFS(ECMWF-Integrated Forecasting System)全球數(shù)值天氣預(yù)報(bào)系統(tǒng)的同化模塊采用RTTOV作為快速輻射傳輸模型,美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(National Centers for Environmental Prediction,NCEP)研發(fā)的NCEP-GFS(NCEP-Global Forecast System)全球數(shù)值天氣預(yù)報(bào)系統(tǒng)的同化模塊中采用CRTM作為快速輻射傳輸模型。

    快速輻射傳輸模型被用于快速模擬微波溫度計(jì)的過(guò)程中,需要依賴(lài)于微波溫度計(jì)在整個(gè)設(shè)計(jì)通帶上的光譜響應(yīng)函數(shù)和大氣吸收模型來(lái)預(yù)先計(jì)算基準(zhǔn)大氣狀態(tài)下測(cè)量到的大氣微波光學(xué)厚度,而后計(jì)算不同大氣狀態(tài)相比于基準(zhǔn)大氣狀態(tài)的微波光學(xué)厚度變化量,進(jìn)而模擬逐層大氣微波輻射及消光作用,最終模擬出微波溫度計(jì)觀測(cè)到的大氣整層微波輻射[4]。光譜響應(yīng)函數(shù)準(zhǔn)確與否,直接關(guān)系到微波輻射模擬的精確度。

    光譜響應(yīng)函數(shù)是微波遙感器的重要硬件參數(shù)特征,通常由中心頻點(diǎn)、頻帶寬度等參數(shù)定義,反映了微波遙感器在設(shè)計(jì)通道內(nèi)對(duì)電磁波的響應(yīng)能力。通常情況下,微波遙感器在設(shè)計(jì)過(guò)程中光譜響應(yīng)被設(shè)計(jì)成方波響應(yīng),即在設(shè)計(jì)通帶范圍內(nèi),各頻點(diǎn)對(duì)電磁波能量的測(cè)量響應(yīng)能力均勻一致。目前國(guó)際上主要的數(shù)值天氣預(yù)報(bào)系統(tǒng)ECMWF-IFS和NCEP-GFS等的同化模塊中,均采用理想方波響應(yīng)作為包括微波溫度計(jì)在內(nèi)的微波遙感器的光譜響應(yīng)函數(shù)[4]。

    近些年隨著對(duì)微波遙感器響應(yīng)測(cè)量及特征分析的研究深入,已有部分針對(duì)微波遙感真實(shí)光譜響應(yīng)的研究工作開(kāi)展,例如搭載在S-NPP(Suomi-National Polar-orbiting Partnership)衛(wèi)星上的與MWTS-2同類(lèi)型儀器ATMS(Advanced Technology Microwave Sounder)其微波光譜響應(yīng)特征對(duì)微波輻射模擬影響已經(jīng)得到部分分析[5],但針對(duì)風(fēng)云三號(hào)的微波遙感器的光譜響應(yīng)分析的詳細(xì)研究還較為少見(jiàn),筆者針對(duì)大氣完全吸收通道已開(kāi)展部分分析工作[4]。但是,針對(duì)同時(shí)反映地表和大氣吸收作用的混合通道分析工作少見(jiàn)開(kāi)展,真實(shí)光譜響應(yīng)對(duì)微波輻射模擬誤差的海陸分布特征如何等問(wèn)題還缺乏相關(guān)研究。本文以同時(shí)反映地表和大氣輻射特征的FY-3D MWTS-2第四通道為例,對(duì)其光譜響應(yīng)特征在國(guó)產(chǎn)數(shù)值天氣預(yù)報(bào)系統(tǒng)中對(duì)微波輻射模擬影響開(kāi)展相關(guān)研究,揭示地表和大氣混合微波輻射通道的真實(shí)光譜響應(yīng)函數(shù)在數(shù)值天氣預(yù)報(bào)系統(tǒng)中對(duì)微波輻射模擬的影響,并分析這種影響的海陸差異,從而為相關(guān)通道輻射資料的同化工作提供參考。

    1 實(shí)驗(yàn)部分

    中國(guó)氣象局全球預(yù)報(bào)系統(tǒng)(China Meteorological Administration-Global Forecast System,CMA-GFS)是中國(guó)氣象局自主研發(fā)的全球數(shù)值天氣預(yù)報(bào)系統(tǒng),其自帶有基于四維變分(4DVar)技術(shù)的同化模塊。本研究采用CMA-GFS作為模擬試驗(yàn)系統(tǒng)。

    CMA-GFS與ECMWF-IFS類(lèi)似,均采用RTTOV作為快速輻射傳輸模擬內(nèi)核。RTTOV模擬誤差大小決定了CMA-GFS同化資料的能力。在RTTOV中采用預(yù)先訓(xùn)練好的光學(xué)厚度擬合系數(shù),結(jié)合大氣實(shí)際狀態(tài)曲線和光學(xué)厚度擬合預(yù)測(cè)函數(shù)集,來(lái)實(shí)現(xiàn)微波輻射的快速模擬。光學(xué)厚度擬合預(yù)測(cè)函數(shù)自RTTOV v7開(kāi)始一直使用固定函數(shù)集,本研究中沿用自RTTOV v7以來(lái)一直使用的擬合預(yù)測(cè)函數(shù)集。預(yù)測(cè)函數(shù)集是一組以大氣狀態(tài)參數(shù)為自變量的函數(shù)組合。

    光學(xué)厚度擬合系數(shù)反映了需要模擬的遙感器的光譜響應(yīng)特征,其計(jì)算流程如圖1所示。本研究中,分別用真實(shí)光譜響應(yīng)函數(shù)和理想光譜響應(yīng)函數(shù)分別輸入進(jìn)圖1所示的流程中,計(jì)算真實(shí)光譜響應(yīng)的光學(xué)厚度擬合系數(shù)和理想光譜響應(yīng)的光學(xué)厚度擬合系數(shù)。通過(guò)這兩組不同的光學(xué)厚度擬合系數(shù)再輸入進(jìn)CMA-GFS的同化系統(tǒng)中,可以分別計(jì)算得到真實(shí)光譜響應(yīng)和理想光譜響應(yīng)在相同時(shí)空區(qū)間內(nèi)的觀測(cè)背景差(observation minus background,OMB)。OMB通常在數(shù)值天氣預(yù)報(bào)系統(tǒng)的同化模塊中,用于對(duì)同化模塊中的快速輻射傳輸模型模擬微波輻射能力的評(píng)估。OMB的絕對(duì)值越大,說(shuō)明快速輻射傳輸模型的模擬值相比衛(wèi)星遙感觀測(cè)值的偏移量越大,模擬精度較差;而當(dāng)OMB的絕對(duì)值越小,則說(shuō)明快速輻射傳輸模型模擬的輻射值越接近于衛(wèi)星觀測(cè)值,則快速輻射傳輸模型的準(zhǔn)確性越高。

    圖1 光學(xué)厚度擬合系數(shù)計(jì)算流程Fig.1 Calculation process of optical depth fitting parameters

    大氣吸收模型主要采用毫米波傳播模型(millimeter propagation model,MPM)來(lái)計(jì)算大氣中各類(lèi)氣態(tài)分子成分對(duì)微波的吸收作用,大氣廓線集合采用ECMWF的54層大氣廓線訓(xùn)練集,該訓(xùn)練集可有效表征全球各類(lèi)典型的大氣狀態(tài),確保訓(xùn)練獲得的光學(xué)厚度擬合系數(shù)具有全球區(qū)域代表性,能夠適用于各類(lèi)典型的大氣狀態(tài)下微波輻射的快速計(jì)算[6-7]。擬合方法采用最小二乘法對(duì)預(yù)測(cè)函數(shù)值集和典型氣壓層的大氣光學(xué)厚度進(jìn)行擬合,并獲得擬合系數(shù)(圖1)。

    在CMA-GFS采用的RTTOV快速微波輻射模擬模塊中,分別使用真實(shí)光譜響應(yīng)和理想光譜響應(yīng)獲得的擬合系數(shù),使用NCEP GFS分析場(chǎng)作為本研究中采用FY-3D MWTS-2的真實(shí)光譜響應(yīng)函數(shù)作為輸入,計(jì)算典型氣壓層的光學(xué)厚度值。預(yù)測(cè)函數(shù)采用RTTOV中使用的針對(duì)微波遙感器的預(yù)測(cè)函數(shù)。訓(xùn)練集采用ECMWF的54層訓(xùn)練集。

    FY-3D MWTS-2是典型的星載微波溫度探測(cè)儀,其具有13個(gè)探測(cè)通道,工作頻率、工作帶寬、極化方式等如表1所示。

    表1 FY3D MWTS-2 通道參數(shù)Table 1 Channel parameters of FY3D MWTS-2

    根據(jù)FY-3D MWTS-2各通道權(quán)重函數(shù)模擬,其第4通道所接收到的微波輻射既有大氣層的微波輻射也有地表發(fā)射的微波輻射貢獻(xiàn),本研究以FY-3D MWTS-2第4通道為例,開(kāi)展真實(shí)光譜響應(yīng)對(duì)CMA-GFS中微波輻射模擬影響分析,并分析在海上和陸地上的影響特征。本文所采用的試驗(yàn)數(shù)據(jù)為2019年FY-3D MWTS-2觀測(cè)數(shù)據(jù)。

    如圖2所示,為FY-3D MWTS-2第4通道的光譜響應(yīng)函數(shù),紅色實(shí)線為表1中通道參數(shù)所對(duì)應(yīng)的理想光譜響應(yīng),即MWTS-2第4通道設(shè)計(jì)的目標(biāo)狀態(tài)。藍(lán)色實(shí)線為測(cè)量?jī)x器實(shí)測(cè)到的MWTS-2第4通道真實(shí)的光譜響應(yīng)。對(duì)比紅色和藍(lán)色實(shí)線可以發(fā)現(xiàn),真實(shí)光譜響應(yīng)在低頻端具有一定的帶寬溢出,同時(shí)在通帶內(nèi),響應(yīng)并不呈現(xiàn)完全平坦的狀態(tài),在低頻段有一定的衰減,而在接近通帶的高頻端出現(xiàn)迅速衰減。

    圖2 FY-3D MWTS-2第4通道理想光譜響應(yīng)(方波響應(yīng))與真實(shí)光譜響應(yīng)函數(shù)Fig.2 Ideal spectral response and real spectral response function of the forth channel of FY-3D MWTS-2

    2 結(jié)果與討論

    利用真實(shí)光譜響應(yīng)函數(shù)和理想光譜響應(yīng)函數(shù)分別模擬出的RTTOV快速微波輻射模擬模型,分別嵌入進(jìn)CMA-GFS的四維變分同化系統(tǒng),計(jì)算兩種模擬模型的OMB。

    (1)

    圖3 FY-3D MWTS-2第4通道理想光譜響應(yīng)微波輻射模擬誤差分析Fig.3 Microwave radiance simulation error analysis of the forth channel of FY-3D MWTS-2 by using ideal spectral response

    圖4 FY-3D MWTS-2第4通道真實(shí)光譜響應(yīng)微波輻射模擬誤差分析Fig.4 Microwave radiance simulation error analysis of the forth channel of FY-3D MWTS-2 by using real spectral response

    通過(guò)分析真實(shí)光譜響應(yīng)函數(shù)模擬得到的OMB與理想光譜響應(yīng)函數(shù)模擬得到的OMB的正態(tài)分布特征可見(jiàn),真實(shí)光譜響應(yīng)函數(shù)模擬得到的OMB其正態(tài)分布的數(shù)學(xué)期望和方差相比理想光譜響應(yīng)函數(shù)模擬結(jié)果更小,顯示真實(shí)光譜響應(yīng)函數(shù)的應(yīng)用可以有效降低模擬輻射亮溫與真實(shí)值之間的差距(如圖5所示)。通過(guò)分析兩種光譜響應(yīng)情況下,模擬得到的OMB值隨緯度變化特征可見(jiàn),除了在南半球靠近極地的極少部分地區(qū)外,真實(shí)光譜響應(yīng)函數(shù)模擬OMB結(jié)果均小于理想光譜響應(yīng)函數(shù)模擬結(jié)果(見(jiàn)圖6),幾乎所有的真實(shí)光譜響應(yīng)與理想光譜響應(yīng)模擬OMB差都小于0,樣本點(diǎn)絕大多數(shù)分布于零值虛線的下方。

    圖5 FY-3D MWTS-2第4通道理想光譜響應(yīng)(方波響應(yīng))與真實(shí)光譜響應(yīng)模擬OMB正態(tài)分布擬合Fig.5 Normal distribution fitting of OMB by ideal and real spectral responses of the forth channel of FY-3D MWTS-2

    圖6 FY-3D MWTS-2第4通道真實(shí)光譜響應(yīng)與理想光譜響應(yīng)(方波響應(yīng))模擬微波輻射OMB差緯向分布特征Fig.6 Latitudinal distribution of microwave radiance OMB of the forth channel of FY-3D MWTS-2 by ideal and real spectral response

    由于FY3D MWTS第四通道同時(shí)反映大氣和地表輻射的特征,因此,針對(duì)樣本點(diǎn)模擬得到的真實(shí)光譜響應(yīng)函數(shù)模擬得到的OMB與理想光譜響應(yīng)函數(shù)模擬得到的OMB,按照全球(表2)、海洋(表3)、陸地(表4)及海陸交界地區(qū)(表5)等不同地表情況分類(lèi),統(tǒng)計(jì)分析不同海陸類(lèi)型下OMB分布的均方根誤差(root mean square error,RMSE)[式(2)]、標(biāo)準(zhǔn)差(standard deviation,STDV)[式(3)]、平均絕對(duì)誤差(mean absolute error,MAE)[式(4)]。定義如式(2)—式(4)所示

    表2 全球真實(shí)光譜響應(yīng)影響OMB統(tǒng)計(jì)Table 2 Global statistical characteristics of OMB by ideal and real spectral responses

    表3 陸地地區(qū)真實(shí)光譜響應(yīng)影響OMB統(tǒng)計(jì)Table 3 Terrestrial region statistical characteristics of OMB by ideal and real spectral responses

    表4 海洋地區(qū)真實(shí)光譜響應(yīng)影響OMB統(tǒng)計(jì)Table 4 Oceanic region statistical characteristics of OMB by ideal and real spectral responses

    表5 海陸交界地區(qū)真實(shí)光譜響應(yīng)影響OMB統(tǒng)計(jì)Table 5 Land-sea border region statistical characteristics of OMB by ideal and real spectral responses

    (2)

    (3)

    (4)

    其中N為總樣本數(shù)量,其余定義同式(1)。

    針對(duì)全球總體而言,采用真實(shí)光譜響應(yīng)函數(shù)模擬得到的OMB相對(duì)于理想光譜響應(yīng),無(wú)論是RMSE、STDV還是MAE均有所改善,其中RMSE減小約2.17%,MAE減小約4.6%,STDV減小約0.75%,可見(jiàn)在數(shù)值天氣預(yù)報(bào)系統(tǒng)中采用真實(shí)光譜響應(yīng)函數(shù)代替理想光譜響應(yīng)函數(shù)來(lái)計(jì)算微波光學(xué)厚度擬合系數(shù),可以改善快速微波輻射傳輸模型對(duì)輻射亮溫的模擬能力。

    在陸地地區(qū),采用真實(shí)光譜響應(yīng)函數(shù)來(lái)計(jì)算快速微波輻射傳輸模型的觀測(cè)算子后,RMSE誤差減小約1.80%,STDV減小約0.95%,MAE減小約1.91%;在海洋地區(qū),采用真實(shí)光譜響應(yīng)函數(shù)來(lái)計(jì)算快速微波輻射傳輸模型的觀測(cè)算子后,RMSE誤差減小約2.6%,STDV減小約0.54%,MAE減小約3.64%;在海陸交界地區(qū),采用真實(shí)光譜響應(yīng)函數(shù)來(lái)計(jì)算快速微波輻射傳輸模型的觀測(cè)算子后,RMSE誤差減小約4.72%,STDV減小約0.79%,MAE減小約6.25%??梢?jiàn),采用真實(shí)光譜響應(yīng)函數(shù)后對(duì)海陸交界地區(qū)的微波輻射模擬誤差改善效果最為明顯,海洋地區(qū)其次,陸地地區(qū)最小,但總體均有改善效果。

    3 結(jié) 論

    真實(shí)光譜響應(yīng)作為微波遙感器的重要探測(cè)特征,相比于傳統(tǒng)的中心頻率和工作帶寬等理想描述參數(shù),更能反映微波遙感器在整個(gè)響應(yīng)通帶上的連續(xù)探測(cè)性能。在國(guó)產(chǎn)數(shù)值天氣預(yù)報(bào)系統(tǒng)中,真實(shí)光譜響應(yīng)尚未被業(yè)務(wù)應(yīng)用,相關(guān)研究也較少。筆者通過(guò)在中國(guó)氣象局自主研發(fā)的數(shù)值天氣預(yù)報(bào)系統(tǒng)CMA-GFS中,采用國(guó)產(chǎn)FY-3D MWTS-2真實(shí)光譜響應(yīng)代替?zhèn)鹘y(tǒng)的理想光譜響應(yīng),計(jì)算新的代表微波遙感器探測(cè)特征的微波光學(xué)后續(xù)擬合系數(shù),經(jīng)過(guò)對(duì)比試驗(yàn),發(fā)現(xiàn)在CMA-GFS中,基于真實(shí)光譜響應(yīng)的快速微波輻射模擬模塊對(duì)于全球微波輻射的模擬能力有明顯提升,其中在海陸交界地區(qū)提升最明顯,海洋地區(qū)次之。研究結(jié)果顯示真實(shí)光譜響應(yīng)函數(shù)可改善數(shù)值天氣預(yù)報(bào)初始場(chǎng)精度。未來(lái)在CMA-GFS中采用真實(shí)光譜響應(yīng)對(duì)預(yù)報(bào)能力改善的試驗(yàn)工作有待進(jìn)一步開(kāi)展。

    猜你喜歡
    響應(yīng)函數(shù)天氣預(yù)報(bào)光學(xué)
    天氣預(yù)報(bào)員
    不同探源距下241Am的α譜模擬與解析
    一類(lèi)具有Beddington-DeAngelis響應(yīng)函數(shù)的階段結(jié)構(gòu)捕食模型的穩(wěn)定性
    滑輪組的裝配
    光學(xué)常見(jiàn)考題逐個(gè)擊破
    相機(jī)響應(yīng)函數(shù)定標(biāo)的正則化方法
    天氣預(yù)報(bào)的前世今生
    中期天氣預(yù)報(bào)
    克服動(dòng)態(tài)問(wèn)題影響的相機(jī)響應(yīng)函數(shù)標(biāo)定
    小小天氣預(yù)報(bào)員
    免费观看性生交大片5| 老司机在亚洲福利影院| www.自偷自拍.com| 啦啦啦 在线观看视频| 亚洲精品日韩在线中文字幕| 国产激情久久老熟女| 国语对白做爰xxxⅹ性视频网站| 久久精品久久久久久久性| 考比视频在线观看| 老鸭窝网址在线观看| 桃花免费在线播放| 亚洲国产最新在线播放| 涩涩av久久男人的天堂| 色网站视频免费| 日韩大码丰满熟妇| 久久综合国产亚洲精品| 国产精品国产三级国产专区5o| 久久婷婷青草| 国产免费福利视频在线观看| 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 国产深夜福利视频在线观看| 五月开心婷婷网| 国产 一区精品| 在线观看国产h片| 午夜激情av网站| 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 在现免费观看毛片| 欧美精品高潮呻吟av久久| 亚洲av在线观看美女高潮| 国产男人的电影天堂91| 在线观看三级黄色| 亚洲天堂av无毛| 又大又黄又爽视频免费| 久久久久久久久久久免费av| 国产成人精品在线电影| 欧美在线一区亚洲| 亚洲精品日韩在线中文字幕| 69精品国产乱码久久久| 国产麻豆69| 熟女少妇亚洲综合色aaa.| 欧美精品av麻豆av| 综合色丁香网| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区精品视频观看| 妹子高潮喷水视频| 午夜精品国产一区二区电影| 美女中出高潮动态图| 亚洲,欧美,日韩| 美女脱内裤让男人舔精品视频| 欧美激情高清一区二区三区 | 又黄又粗又硬又大视频| 咕卡用的链子| 亚洲五月色婷婷综合| 日韩一卡2卡3卡4卡2021年| 色吧在线观看| 天美传媒精品一区二区| 国产片特级美女逼逼视频| 婷婷色av中文字幕| 亚洲欧美色中文字幕在线| www.av在线官网国产| 老鸭窝网址在线观看| 亚洲熟女精品中文字幕| 亚洲国产欧美在线一区| 黄片播放在线免费| 美女福利国产在线| 最近2019中文字幕mv第一页| 亚洲国产欧美网| 看免费成人av毛片| 亚洲国产看品久久| 一个人免费看片子| 久久久国产欧美日韩av| 国产精品麻豆人妻色哟哟久久| 18禁国产床啪视频网站| 女人久久www免费人成看片| 欧美日韩亚洲综合一区二区三区_| 男女边摸边吃奶| 精品一区二区三区四区五区乱码 | 亚洲欧美一区二区三区黑人| 1024视频免费在线观看| 美女主播在线视频| 免费久久久久久久精品成人欧美视频| 精品少妇内射三级| 最近2019中文字幕mv第一页| 国产日韩欧美亚洲二区| 国产免费现黄频在线看| 日韩一卡2卡3卡4卡2021年| 精品久久久精品久久久| 午夜福利乱码中文字幕| 欧美精品一区二区免费开放| 色播在线永久视频| 欧美97在线视频| 一级片免费观看大全| 国产激情久久老熟女| 亚洲国产最新在线播放| 观看美女的网站| 无限看片的www在线观看| 女的被弄到高潮叫床怎么办| 精品国产露脸久久av麻豆| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 国产精品二区激情视频| 性少妇av在线| 亚洲欧美一区二区三区国产| 欧美激情极品国产一区二区三区| 精品久久久久久电影网| 一边摸一边抽搐一进一出视频| 两个人免费观看高清视频| 十八禁人妻一区二区| 我的亚洲天堂| 97在线人人人人妻| 人成视频在线观看免费观看| 中国国产av一级| 国产精品久久久久久精品古装| 久久精品国产亚洲av涩爱| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 国产精品国产av在线观看| 大香蕉久久成人网| 国产成人a∨麻豆精品| 无遮挡黄片免费观看| 午夜福利网站1000一区二区三区| 欧美日韩视频高清一区二区三区二| 悠悠久久av| videosex国产| 高清av免费在线| 欧美成人午夜精品| 伦理电影大哥的女人| 桃花免费在线播放| 亚洲五月色婷婷综合| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区二区精品视频观看| 亚洲av日韩精品久久久久久密 | 日本猛色少妇xxxxx猛交久久| 少妇人妻精品综合一区二区| 精品酒店卫生间| 桃花免费在线播放| 免费黄色在线免费观看| 90打野战视频偷拍视频| 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 国产成人精品福利久久| 午夜激情久久久久久久| 丝袜在线中文字幕| 国产成人精品久久久久久| 丰满乱子伦码专区| 亚洲av成人不卡在线观看播放网 | 90打野战视频偷拍视频| 国产精品欧美亚洲77777| 丝袜在线中文字幕| 国产成人精品无人区| 亚洲情色 制服丝袜| 国产深夜福利视频在线观看| 丝袜人妻中文字幕| 午夜91福利影院| 精品一品国产午夜福利视频| bbb黄色大片| 国产精品熟女久久久久浪| 亚洲视频免费观看视频| 十分钟在线观看高清视频www| 亚洲欧美一区二区三区国产| 捣出白浆h1v1| 欧美人与性动交α欧美精品济南到| 中文字幕色久视频| 国产成人精品福利久久| 日韩av在线免费看完整版不卡| 各种免费的搞黄视频| 中文字幕高清在线视频| 亚洲欧美精品综合一区二区三区| 三上悠亚av全集在线观看| 精品一区二区三卡| 久久精品人人爽人人爽视色| 男女边摸边吃奶| 亚洲av男天堂| 精品卡一卡二卡四卡免费| 热re99久久精品国产66热6| 亚洲人成网站在线观看播放| 免费观看av网站的网址| 国产成人免费观看mmmm| 国产1区2区3区精品| 一级黄片播放器| 久久久久精品久久久久真实原创| 久久久久精品性色| 波野结衣二区三区在线| 欧美精品av麻豆av| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看黄色视频的| 色婷婷久久久亚洲欧美| 美女高潮到喷水免费观看| 2018国产大陆天天弄谢| 九九爱精品视频在线观看| 天天操日日干夜夜撸| 激情五月婷婷亚洲| 最近的中文字幕免费完整| av在线观看视频网站免费| 99精国产麻豆久久婷婷| 只有这里有精品99| 观看美女的网站| 精品国产国语对白av| 国产精品久久久久成人av| 9热在线视频观看99| 一区二区三区激情视频| 纯流量卡能插随身wifi吗| 极品少妇高潮喷水抽搐| 午夜免费观看性视频| 交换朋友夫妻互换小说| 一本色道久久久久久精品综合| 欧美日韩国产mv在线观看视频| 女人高潮潮喷娇喘18禁视频| 国产精品av久久久久免费| 成人国产av品久久久| av一本久久久久| 久久久久国产精品人妻一区二区| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| xxxhd国产人妻xxx| 纯流量卡能插随身wifi吗| 亚洲成色77777| e午夜精品久久久久久久| 国产精品嫩草影院av在线观看| 考比视频在线观看| av在线观看视频网站免费| 观看美女的网站| 在线亚洲精品国产二区图片欧美| 夫妻性生交免费视频一级片| 最近最新中文字幕免费大全7| 久久久久网色| 两个人免费观看高清视频| 一区二区三区精品91| 最近中文字幕高清免费大全6| 无遮挡黄片免费观看| 制服丝袜香蕉在线| 中国三级夫妇交换| 国产成人a∨麻豆精品| 十八禁高潮呻吟视频| 欧美激情高清一区二区三区 | 悠悠久久av| 成人免费观看视频高清| 婷婷色综合大香蕉| 青春草亚洲视频在线观看| 777久久人妻少妇嫩草av网站| 搡老岳熟女国产| 亚洲第一区二区三区不卡| 日韩制服骚丝袜av| 亚洲精品av麻豆狂野| 国产乱来视频区| 人妻一区二区av| 黄色视频不卡| 又大又黄又爽视频免费| 亚洲 欧美一区二区三区| 伊人久久大香线蕉亚洲五| 性高湖久久久久久久久免费观看| 少妇被粗大猛烈的视频| 丁香六月欧美| 欧美黑人精品巨大| 一本—道久久a久久精品蜜桃钙片| 热re99久久精品国产66热6| 亚洲精品aⅴ在线观看| 国产免费视频播放在线视频| 久久久久精品人妻al黑| 青春草亚洲视频在线观看| 日本欧美视频一区| 亚洲一级一片aⅴ在线观看| 如何舔出高潮| 免费人妻精品一区二区三区视频| 热re99久久国产66热| 欧美97在线视频| 最近中文字幕2019免费版| 伊人久久国产一区二区| 久久久欧美国产精品| 天天操日日干夜夜撸| 各种免费的搞黄视频| 侵犯人妻中文字幕一二三四区| www.av在线官网国产| 美女扒开内裤让男人捅视频| 晚上一个人看的免费电影| 51午夜福利影视在线观看| 久久久精品94久久精品| 欧美精品亚洲一区二区| 国产探花极品一区二区| 亚洲av成人精品一二三区| 欧美 亚洲 国产 日韩一| 成年动漫av网址| 国产高清不卡午夜福利| 99热网站在线观看| 国产片特级美女逼逼视频| 天天影视国产精品| 欧美精品一区二区免费开放| 日本av手机在线免费观看| 国产精品久久久久久久久免| 国产探花极品一区二区| 波多野结衣一区麻豆| 久久热在线av| 亚洲精品中文字幕在线视频| 岛国毛片在线播放| 性少妇av在线| 超碰成人久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产欧美一区二区综合| 制服人妻中文乱码| 欧美日韩亚洲高清精品| 久久国产精品大桥未久av| 日本猛色少妇xxxxx猛交久久| 精品亚洲成a人片在线观看| 男人操女人黄网站| 国产亚洲av高清不卡| 亚洲欧美精品自产自拍| 2021少妇久久久久久久久久久| 伦理电影大哥的女人| 91老司机精品| 国产伦理片在线播放av一区| 中国国产av一级| 秋霞伦理黄片| 国产男女内射视频| 久久国产精品男人的天堂亚洲| 亚洲婷婷狠狠爱综合网| 午夜福利,免费看| 激情五月婷婷亚洲| 免费看不卡的av| 国产亚洲欧美精品永久| 亚洲国产看品久久| 亚洲国产中文字幕在线视频| 午夜福利乱码中文字幕| 亚洲国产中文字幕在线视频| 一区二区日韩欧美中文字幕| 日本猛色少妇xxxxx猛交久久| 天天影视国产精品| 又粗又硬又长又爽又黄的视频| 午夜福利视频在线观看免费| 久久久久视频综合| 国产xxxxx性猛交| 日韩大片免费观看网站| 亚洲国产毛片av蜜桃av| 最近最新中文字幕免费大全7| 国产深夜福利视频在线观看| 日韩一卡2卡3卡4卡2021年| 18禁观看日本| 免费在线观看黄色视频的| 如何舔出高潮| 国产国语露脸激情在线看| 18禁动态无遮挡网站| 街头女战士在线观看网站| 麻豆精品久久久久久蜜桃| 天天躁夜夜躁狠狠久久av| 亚洲欧美中文字幕日韩二区| 日韩人妻精品一区2区三区| 在现免费观看毛片| 中文天堂在线官网| 97精品久久久久久久久久精品| 伊人亚洲综合成人网| 在线 av 中文字幕| 久久久久久久久久久久大奶| 精品少妇一区二区三区视频日本电影 | 在线观看一区二区三区激情| 久久久久国产一级毛片高清牌| 在线观看人妻少妇| 久久久久国产一级毛片高清牌| 国产精品麻豆人妻色哟哟久久| 18禁裸乳无遮挡动漫免费视频| 91成人精品电影| 亚洲免费av在线视频| 亚洲伊人色综图| 亚洲国产av影院在线观看| 黑人巨大精品欧美一区二区蜜桃| e午夜精品久久久久久久| 极品少妇高潮喷水抽搐| 久久久久视频综合| 国产又爽黄色视频| 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 久久ye,这里只有精品| 国产av一区二区精品久久| 男女高潮啪啪啪动态图| 亚洲男人天堂网一区| 青青草视频在线视频观看| 在线精品无人区一区二区三| 一二三四中文在线观看免费高清| 成人亚洲精品一区在线观看| 午夜福利影视在线免费观看| 国产av码专区亚洲av| 亚洲欧美清纯卡通| 人成视频在线观看免费观看| 秋霞在线观看毛片| 波多野结衣一区麻豆| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 高清在线视频一区二区三区| 精品久久蜜臀av无| 人妻 亚洲 视频| 欧美另类一区| 中文字幕人妻丝袜一区二区 | 欧美日韩一区二区视频在线观看视频在线| 免费日韩欧美在线观看| 男女边摸边吃奶| 国产乱来视频区| 亚洲美女黄色视频免费看| 老汉色∧v一级毛片| av免费观看日本| 免费黄网站久久成人精品| 日本vs欧美在线观看视频| 久久久国产欧美日韩av| 午夜久久久在线观看| 国产av一区二区精品久久| 国产精品一区二区在线不卡| 中文字幕人妻丝袜一区二区 | 久久精品国产亚洲av高清一级| 欧美国产精品va在线观看不卡| 大话2 男鬼变身卡| 亚洲精品成人av观看孕妇| 久久精品亚洲av国产电影网| 精品国产超薄肉色丝袜足j| 哪个播放器可以免费观看大片| 成年av动漫网址| 亚洲av男天堂| h视频一区二区三区| 搡老乐熟女国产| 亚洲第一青青草原| 高清不卡的av网站| 免费观看av网站的网址| 最近的中文字幕免费完整| 91国产中文字幕| 亚洲精品中文字幕在线视频| 国产一区二区在线观看av| 亚洲国产毛片av蜜桃av| 激情五月婷婷亚洲| 亚洲精品美女久久av网站| 99久久综合免费| 精品国产露脸久久av麻豆| 亚洲精品日本国产第一区| 日韩av在线免费看完整版不卡| 丰满饥渴人妻一区二区三| 日本午夜av视频| 丝瓜视频免费看黄片| 日韩大片免费观看网站| av视频免费观看在线观看| 久久人人爽av亚洲精品天堂| 免费观看a级毛片全部| 女的被弄到高潮叫床怎么办| 中国三级夫妇交换| 精品人妻在线不人妻| 成人漫画全彩无遮挡| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美色中文字幕在线| 午夜久久久在线观看| 18禁国产床啪视频网站| 在线观看免费高清a一片| 中文字幕人妻丝袜制服| 亚洲五月色婷婷综合| 男人操女人黄网站| videos熟女内射| av网站在线播放免费| 亚洲精华国产精华液的使用体验| 高清不卡的av网站| 精品免费久久久久久久清纯 | 亚洲第一青青草原| 中文字幕高清在线视频| 成人毛片60女人毛片免费| 青春草国产在线视频| 日本一区二区免费在线视频| 欧美精品高潮呻吟av久久| www.熟女人妻精品国产| 在线天堂中文资源库| 午夜福利影视在线免费观看| 亚洲av中文av极速乱| 如何舔出高潮| 可以免费在线观看a视频的电影网站 | 最新在线观看一区二区三区 | 久久久久久久大尺度免费视频| 久久亚洲国产成人精品v| 久久99一区二区三区| 久久精品国产综合久久久| 啦啦啦 在线观看视频| 伦理电影免费视频| 日韩一区二区三区影片| 男的添女的下面高潮视频| 天天躁日日躁夜夜躁夜夜| 成年人免费黄色播放视频| 在线观看国产h片| 美女中出高潮动态图| 国产精品亚洲av一区麻豆 | 国产精品 国内视频| 国产片特级美女逼逼视频| 欧美xxⅹ黑人| 中文字幕精品免费在线观看视频| 国产男女超爽视频在线观看| 王馨瑶露胸无遮挡在线观看| 精品酒店卫生间| 高清不卡的av网站| 欧美激情极品国产一区二区三区| 啦啦啦在线免费观看视频4| 在线精品无人区一区二区三| 欧美日韩视频高清一区二区三区二| 国产毛片在线视频| 国产一区二区 视频在线| 久久精品国产综合久久久| videosex国产| 下体分泌物呈黄色| 看非洲黑人一级黄片| 巨乳人妻的诱惑在线观看| 黄片播放在线免费| 狠狠婷婷综合久久久久久88av| 一级黄片播放器| 日本wwww免费看| 人妻 亚洲 视频| 卡戴珊不雅视频在线播放| 国产女主播在线喷水免费视频网站| 婷婷色av中文字幕| 18禁裸乳无遮挡动漫免费视频| 18禁国产床啪视频网站| 咕卡用的链子| 汤姆久久久久久久影院中文字幕| 99久久人妻综合| 国产精品三级大全| 成人免费观看视频高清| 国产黄色视频一区二区在线观看| 精品国产乱码久久久久久小说| 亚洲国产看品久久| 一区二区三区四区激情视频| 国产成人精品无人区| 悠悠久久av| 亚洲综合色网址| 老汉色av国产亚洲站长工具| 女性被躁到高潮视频| av在线观看视频网站免费| 极品少妇高潮喷水抽搐| 精品第一国产精品| 肉色欧美久久久久久久蜜桃| 免费看av在线观看网站| xxx大片免费视频| 欧美黑人欧美精品刺激| av不卡在线播放| 成人手机av| 欧美精品高潮呻吟av久久| 天堂8中文在线网| 男女之事视频高清在线观看 | 涩涩av久久男人的天堂| 久久热在线av| 男女之事视频高清在线观看 | 国产精品麻豆人妻色哟哟久久| 人人澡人人妻人| 国产有黄有色有爽视频| 亚洲av综合色区一区| 精品午夜福利在线看| 亚洲精品美女久久av网站| a级毛片黄视频| 欧美日韩一级在线毛片| 一区福利在线观看| 中文字幕高清在线视频| 日韩大码丰满熟妇| 一个人免费看片子| 日本爱情动作片www.在线观看| 成人毛片60女人毛片免费| 三上悠亚av全集在线观看| 多毛熟女@视频| 免费高清在线观看日韩| 丰满饥渴人妻一区二区三| 成年女人毛片免费观看观看9 | 国产伦人伦偷精品视频| 欧美亚洲日本最大视频资源| 一边亲一边摸免费视频| 人人妻人人澡人人爽人人夜夜| 少妇人妻精品综合一区二区| 最新在线观看一区二区三区 | 天堂中文最新版在线下载| 亚洲国产欧美日韩在线播放| 亚洲综合色网址| 伦理电影大哥的女人| 国产亚洲av高清不卡| 国产精品亚洲av一区麻豆 | 免费在线观看视频国产中文字幕亚洲 | 国产伦理片在线播放av一区| 热99国产精品久久久久久7| 看非洲黑人一级黄片| 国产午夜精品一二区理论片| 亚洲精品国产av成人精品| 哪个播放器可以免费观看大片| h视频一区二区三区| 精品国产国语对白av| 黄片无遮挡物在线观看| 欧美在线黄色| 丝袜美足系列| 超色免费av| 日韩视频在线欧美| 嫩草影院入口| 两性夫妻黄色片| h视频一区二区三区| 久久精品国产综合久久久| 制服诱惑二区| 精品人妻在线不人妻| 一级毛片我不卡| 国产 一区精品| 成人午夜精彩视频在线观看| 日韩免费高清中文字幕av| 99精品久久久久人妻精品| 热99久久久久精品小说推荐| 免费在线观看视频国产中文字幕亚洲 | 99国产精品免费福利视频| a级毛片黄视频| av片东京热男人的天堂| 亚洲精品美女久久av网站| 免费黄色在线免费观看| 国产一区亚洲一区在线观看| 久久精品人人爽人人爽视色| 激情视频va一区二区三区| 亚洲精品国产一区二区精华液| 天天添夜夜摸| 老司机靠b影院| 欧美少妇被猛烈插入视频| 久久人妻熟女aⅴ| www.自偷自拍.com| 国产成人av激情在线播放|