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

    SST融合產(chǎn)品在渤黃海的質(zhì)量評(píng)估與修正研究

    2017-03-26 03:10:26曹凱翔靳熙芳孫偉富
    海洋科學(xué) 2017年9期
    關(guān)鍵詞:輻射計(jì)海溫黃海

    曹凱翔, 靳熙芳, 孫偉富, 王 進(jìn)

    ?

    SST融合產(chǎn)品在渤黃海的質(zhì)量評(píng)估與修正研究

    曹凱翔1, 靳熙芳2, 孫偉富1, 王 進(jìn)3

    (1. 國(guó)家海洋局第一海洋研究所, 山東 青島 266061; 2. 國(guó)家海洋局北海預(yù)報(bào)中心, 山東 青島 266061; 3. 青島大學(xué), 山東 青島 266003)

    為獲取高空間分辨率與空間覆蓋率的海表面溫度(Sea Surface Temperature, SST)產(chǎn)品, 基于最優(yōu)插值方法, 對(duì)微波輻射計(jì)WindSat、AMSR-E、ASMR2、HY-2 RM和紅外輻射計(jì)MODIS、AVHRR的SST觀測(cè)數(shù)據(jù)進(jìn)行融合, 生成了一種0.1°空間分辨率的每日SST融合產(chǎn)品, 利用浮標(biāo)數(shù)據(jù)在渤黃海區(qū)域進(jìn)行了精度評(píng)估和修正, 并分析了該區(qū)域的SST時(shí)空分布特征。結(jié)果表明: SST融合產(chǎn)品的在中國(guó)近海的精度為1.1℃, 利用浮標(biāo)數(shù)據(jù)修正后的精度略有改善; 利用修正的SST產(chǎn)品對(duì)渤黃海區(qū)域SST分布特征進(jìn)行了分析, 分析結(jié)果顯示, 渤黃海海域冬季海溫最為均勻, 春季在海水升溫過(guò)程中海溫不均勻性明顯。

    海表溫度; 數(shù)據(jù)融合; 渤黃海; 質(zhì)量評(píng)估; 時(shí)空分布特征

    海表面溫度(Sea Surface Temperature, SST)是海洋研究的基本參數(shù)之一, 其直接影響著大氣和海洋間的熱量、動(dòng)量和水汽交換, 是影響全球水循環(huán)的重要參數(shù), 影響著全球表面的能量收支平衡[1]。基于衛(wèi)星平臺(tái)的被動(dòng)遙感觀測(cè)技術(shù)是全球海表溫度監(jiān)測(cè)的主要手段。目前可進(jìn)行SST觀測(cè)的星載傳感器主要包括微波輻射計(jì)和紅外輻射計(jì)兩類。微波輻射計(jì)觀測(cè)可以穿透云霧, 受氣候影響小, 每天可覆蓋全球近90%的海域, 但是其空間分辨率較低, 為幾十千米。紅外輻射計(jì)工作頻段遠(yuǎn)高于微波輻射計(jì), 因此其空間分辨率可達(dá)千米量級(jí), 但紅外輻射計(jì)受天氣條件影響大, 不能穿透云霧, 單個(gè)紅外輻射計(jì)的每日有效觀測(cè)數(shù)據(jù)僅能覆蓋全球海洋的10%左右。采用一定的數(shù)據(jù)融合方法, 對(duì)星載微波輻射計(jì)和紅外輻射計(jì)SST觀測(cè)結(jié)果進(jìn)行融合, 可以有效地改善微波輻射計(jì)空間分辨率低、紅外輻射計(jì)空間覆蓋率差的問(wèn)題, 從而生成可覆蓋全球海洋的高分辨率SST多源融合產(chǎn)品。近年來(lái)相關(guān)研究受到各國(guó)研究機(jī)構(gòu)的普遍關(guān)注, 如美國(guó)國(guó)家海洋大氣管理局NOAA(National Oceanic and Atmospheric Administration)、NASA下屬噴氣推進(jìn)實(shí)驗(yàn)室JPL(Jet Propulsion Laboratory)、日本氣象廳JMA(Japan Meteorological Agency)、Met Office、遙感系統(tǒng)RSS(Remote Sensing System)均開(kāi)發(fā)了SST融合數(shù)據(jù)產(chǎn)品。這些產(chǎn)品空間分辨率為0.05°~0.25°, 時(shí)間分辨率一般為日均, 并普遍采用浮標(biāo)等實(shí)測(cè)數(shù)據(jù)進(jìn)行了精度驗(yàn)證[2]。

    SST是目前海洋學(xué)界的研究熱點(diǎn), 它直接影響著大氣和海洋間的熱量、動(dòng)量和水汽交換過(guò)程, 對(duì)全球表面的能量收支起重要作用。本文基于最優(yōu)插值方法, 對(duì)微波輻射計(jì)WindSat、AMSR-E、ASMR2、HY-2 RM和紅外輻射計(jì)MODIS、AVHRR數(shù)據(jù)進(jìn)行融合, 制作了0.1°空間分辨率的SST融合產(chǎn)品, 并基于實(shí)測(cè)數(shù)據(jù)在中國(guó)渤海和黃海北部開(kāi)展了SST產(chǎn)品的精度驗(yàn)證和修正研究。這為開(kāi)展海氣相互作用的研究提供了重要的科學(xué)依據(jù), 具有良好的應(yīng)用前景。

    1 SST融合產(chǎn)品介紹

    1.1 SST多源遙感數(shù)據(jù)

    SST融合產(chǎn)品所使用的遙感數(shù)據(jù)包括: WindSat、ASMR-E、AMSR2、HY-2 RM、MODIS(Aqua、Terra)和AVHRR等7個(gè)微波、紅外輻射計(jì)數(shù)據(jù)。

    微波輻射計(jì)AMSR-E、AMSR2和WindSat的每日SST數(shù)據(jù)由遙感系統(tǒng)RSS提供, 其空間分辨率為25 km, 數(shù)據(jù)格式為二進(jìn)制, 數(shù)據(jù)產(chǎn)品中包括每日升軌/降軌的SST、風(fēng)速、水汽和云液水產(chǎn)品[3]。ASMR-E的在軌時(shí)間為2002年3月至2011年10月, 目前已停止工作。AMSR2是AMSR-E的后繼星, 于2012年5月發(fā)射入軌, 至今正常在軌運(yùn)行[4]。WindSat是歷史上第一個(gè)星載全極化微波輻射計(jì), 于2003年1月發(fā)射, 至今在軌正常運(yùn)行。本文提取了該3顆星的SST數(shù)據(jù), 剔除了SST默認(rèn)值和異常值。WindSat和AMSR-E的數(shù)據(jù)樣例如圖1所示。

    圖1 WindSat(a)與AMSR-E(b)SST產(chǎn)品樣例

    HY-2 RM微波輻射計(jì)數(shù)據(jù)由國(guó)家衛(wèi)星海洋應(yīng)用中心提供, 為L(zhǎng)2級(jí)沿軌數(shù)據(jù), 格式為HDF, 其中包含HY-2微波輻射計(jì)反演的SST、水汽含量、云液水含量等數(shù)據(jù)[5]。本文提取了HY-2的原始分辨率SST數(shù)據(jù)并進(jìn)行了數(shù)據(jù)有效性檢查。HY-2 SST數(shù)據(jù)樣例及掃描時(shí)間如圖2所示。

    圖2 HY-2 SST網(wǎng)格數(shù)據(jù)

    MODIS Aqua/Terra數(shù)據(jù)由NASA提供, 空間分辨率為4 km, 數(shù)據(jù)格式為HDF[6]。本文讀取了SST數(shù)據(jù)和相應(yīng)的質(zhì)量標(biāo)記位, 并進(jìn)行了質(zhì)量控制。MODIS/Terra于1999年12月入軌, MODIS/Aqua于2002年3月入軌, 目前兩星皆在軌正常運(yùn)行。MODIS SST數(shù)據(jù)樣例如圖3所示。

    圖3 MODIS SST數(shù)據(jù)樣例

    AVHRR Pathfinder每日SST數(shù)據(jù)由NOAA下屬國(guó)家海洋數(shù)據(jù)中心NODC(National Ocean Data Center)提供, 空間分辨率為4 km, 數(shù)據(jù)格式為NetCDF[6]。本文讀取了該星的SST數(shù)據(jù)和相應(yīng)的質(zhì)量標(biāo)記位, 并進(jìn)行了質(zhì)量控制。AVHRR首星于1981年6月發(fā)射入軌, 目前在軌的是其后繼星。AVHRR SST數(shù)據(jù)樣例如圖4所示。

    圖4 AVHRR SST數(shù)據(jù)樣例

    1.2 SST融合方法

    圖5為融合SST產(chǎn)品的數(shù)據(jù)樣例。

    圖5 SST融合個(gè)例

    2 SST產(chǎn)品在渤黃海的質(zhì)量評(píng)估與修正

    2.1 數(shù)據(jù)時(shí)空匹配

    本文利用渤海和黃海北部(35°~41°N, 117°~ 127°E)6個(gè)浮標(biāo)的實(shí)測(cè)數(shù)據(jù), 對(duì)SST融合產(chǎn)品在該區(qū)域的精度進(jìn)行評(píng)估, 浮標(biāo)數(shù)據(jù)的時(shí)間范圍是2011—2014年。遙感-現(xiàn)場(chǎng)數(shù)據(jù)時(shí)空匹配過(guò)程中, 首先對(duì)浮標(biāo)數(shù)據(jù)進(jìn)行預(yù)處理, 剔除了浮標(biāo)觀測(cè)數(shù)據(jù)中的異常值; 由于浮標(biāo)對(duì)SST的觀測(cè)頻次很高, 因此對(duì)同一浮標(biāo)的觀測(cè)結(jié)果進(jìn)行了日平均; 進(jìn)而采用50 km、1 d的時(shí)空匹配窗口對(duì)浮標(biāo)和SST融合產(chǎn)品進(jìn)行了時(shí)空匹配, 共獲得匹配數(shù)據(jù)4 850組。

    2.2 SST產(chǎn)品精度評(píng)估

    基于以上時(shí)空匹配數(shù)據(jù), 本文計(jì)算了SST融合產(chǎn)品相對(duì)浮標(biāo)數(shù)據(jù)的平均偏差、均方根(root mean square, RMS)誤差和相關(guān)系數(shù), 作為SST產(chǎn)品精度評(píng)估的統(tǒng)計(jì)指標(biāo)。結(jié)果表明, SST融合產(chǎn)品的平均偏差為0.1℃, RMS誤差為1.1℃,相關(guān)系數(shù)為0.99。圖6是SST反演誤差的散點(diǎn)圖和誤差分布直方圖。

    圖6 SST產(chǎn)品與浮標(biāo)數(shù)據(jù)的比較(a)和誤差分布直方圖(b)

    一般而言, 星載微波輻射計(jì)和紅外輻射計(jì)在大洋區(qū)域的SST反演精度優(yōu)于1.0℃, HY-2微波輻射計(jì)的SST反演精度稍差, 約1.7℃[8]。本文生產(chǎn)的SST融合產(chǎn)品在中國(guó)近海的誤差略大于大洋區(qū)域, 其原因可能在于輻射計(jì)在近海區(qū)域的觀測(cè)數(shù)據(jù)容易受到陸地污染和無(wú)線電射頻干擾等因素的影響, 導(dǎo)致SST反演誤差增大[9]。

    2.3 SST融合產(chǎn)品修正

    本節(jié)利用SST融合產(chǎn)品和浮標(biāo)匹配數(shù)據(jù), 對(duì)SST融合產(chǎn)品進(jìn)行修正, 以改善其在中國(guó)近海的精度。首先將4 850組匹配數(shù)據(jù)隨機(jī)平均分為建模數(shù)據(jù)集和檢驗(yàn)數(shù)據(jù)集, 利用建模數(shù)據(jù)集建立SST融合產(chǎn)品與現(xiàn)場(chǎng)數(shù)據(jù)的經(jīng)驗(yàn)關(guān)系, 對(duì)融合產(chǎn)品進(jìn)行修正; 利用檢驗(yàn)數(shù)據(jù)集對(duì)修正后的SST產(chǎn)品進(jìn)行檢驗(yàn), 并與修正前的結(jié)果進(jìn)行比較。

    由圖6可知, SST融合數(shù)據(jù)和浮標(biāo)觀測(cè)數(shù)據(jù)在散點(diǎn)圖中基本都分布在1︰1線附近, 因此本文選擇線性方程對(duì)SST融合數(shù)據(jù)和浮標(biāo)數(shù)據(jù)進(jìn)行擬合, 得到線性擬合方程為:

    其中,Buoy為浮標(biāo)觀測(cè)數(shù)據(jù),Fusion為融合數(shù)據(jù)。

    利用以上回歸方程對(duì)檢驗(yàn)數(shù)據(jù)集中的SST融合結(jié)果進(jìn)行了修正, 并與浮標(biāo)數(shù)據(jù)進(jìn)行比較。結(jié)果表明, 修正后的SST融合產(chǎn)品的平均偏差為0.01℃, 明顯優(yōu)于修正前的結(jié)果, RMS誤差為1.0℃, 略優(yōu)于修正前的結(jié)果。下圖為修正后的SST產(chǎn)品與浮標(biāo)實(shí)測(cè)數(shù)據(jù)的比較。

    圖7 修正后的SST產(chǎn)品與浮標(biāo)的比較

    3 渤黃海區(qū)域SST分布特征分析

    本文利用以上制作的SST融合數(shù)據(jù), 以中國(guó)渤海和黃海北部(35°~41°N, 117°~127°E)為研究區(qū)域, 對(duì)該海域SST的時(shí)空分布特征進(jìn)行了分析, 計(jì)算了該海域2004—2014年各年4季的海表溫度平均值和標(biāo)準(zhǔn)差, 其中春季取3—5月, 夏季取6—8月, 秋季取9—11月, 冬季取12月—次年2月, 結(jié)果如圖8所示。

    圖8 2004~2014年海表溫度平均值和標(biāo)準(zhǔn)差

    圖中可見(jiàn), 渤黃海區(qū)域夏季海溫一般在22℃以上, 其中2011年夏季的平均海溫為20.9℃, 較其他年份同期海溫低1℃以上; 秋季海溫一般在18.5℃以上, 其中2011年秋季平均海溫18.1℃較往年略有偏低; 該海域春季平均海溫在2010—2011年也出現(xiàn)一個(gè)低值, 約7.1℃; 冬季海溫比較穩(wěn)定, 為3~3.5℃。海溫的標(biāo)準(zhǔn)差反映了海水溫度的均勻程度, 由圖8可知, 該海域春季海溫標(biāo)準(zhǔn)差最大, 為1.2℃以上, 表明海水在升溫過(guò)程中溫度梯度比較明顯; 冬季海溫標(biāo)準(zhǔn)差最小, 一般在0.6℃以下, 表明海溫此時(shí)最為均勻; 2011年夏季海溫的標(biāo)準(zhǔn)差明顯低于2010年和2012年, 表明上述2011年夏季平均海溫低于常年的現(xiàn)象是由于海溫整體偏低導(dǎo)致。

    4 結(jié)論

    本文利用國(guó)內(nèi)外多個(gè)微波、紅外輻射計(jì)的SST產(chǎn)品, 基于最優(yōu)插值方法制作了一種0.1°空間分辨率的全球每日SST產(chǎn)品; 并基于浮標(biāo)實(shí)測(cè)數(shù)據(jù)在我國(guó)渤黃海區(qū)域?qū)ST產(chǎn)品進(jìn)行了精度評(píng)估與修正, 最后對(duì)該海域11年的海溫分布特征進(jìn)行了分析, 得到結(jié)論如下: (1)相對(duì)于浮標(biāo)實(shí)測(cè)數(shù)據(jù), SST融合產(chǎn)品在渤黃海區(qū)域的偏差為0.1℃, RMS誤差為1.1℃; 由于輻射計(jì)在近海區(qū)域的數(shù)據(jù)容易受到陸地污染和無(wú)線電射頻干擾, 融合產(chǎn)品在近海的精度低于大洋區(qū)域。(2)本文利用浮標(biāo)實(shí)測(cè)SST數(shù)據(jù)對(duì)SST融合產(chǎn)品進(jìn)行了修正, 修正后SST融合產(chǎn)品的偏差為0.01℃, RMS誤差為1.0℃, 較修正前有所改善。(3)利用SST融合產(chǎn)品對(duì)渤黃海海域的SST分布特征進(jìn)行了分析, 結(jié)果表明, 渤黃海海域夏季水溫一般高于22℃, 冬季水溫一般為3~3.5℃; 該海域冬季海溫均勻, 標(biāo)準(zhǔn)差小于0.6℃; 春季海水升溫過(guò)程中海溫最不均勻, 海溫標(biāo)準(zhǔn)差大于1.2℃。

    [1] Dash P, Ignatov A, Martin M, et al. Group for high resolution sea surface temperature (GHRSST) analysis field sinter-comparisons—Part 2: Near real time web-based level 4 SST Quality Monitor (L4-SQUAM)[J]. Deep-Sea Research II, 2012, 80: 31-43.

    [2] Hosoda K, Qin Huiling. Algorithm for estimating sea surface temperatures based on aqua/MODIS global ocean data. 1. Developmentand validation of the algorithm [J]. Journal of Oceanography, 2011, 67 (1): 135-145.

    [3] H?yer J L, Karagali I, Dybkj?r G, et al. Multi sensor validation and error characteristics of arctic satellite sea surface temperature observations [J]. Remote Sensing of Environment, 2012, 121: 335-346.

    [4] Kurihara Y, Sakurai T, Kuragano T. Global daily sea surface temperature analysis using data from satellite microwave radiometer, satellite infrared radiometer and in-situ observations [J]. Weather Bull, 2006, 73: 1-18.

    [5] Martin M, Dash P, Ignatov A, et al. Group for high resolution sea surface temperature (GHRSST) analysis field sinter-comparisons. Part 1: A GHRSST multi-product ensemble (GMPE) [J]. Deep-Sea Research II, 2012, 80: 21-30.

    [6] Meissner T, Wentz F. Ocean Retrievals for WindSat[J]. IEEE Microrad, 2006, 7: 119-124.

    [7] Reynolds R, Smith T. Improved global sea surface temperature analyses using optimum interpolation[J]. Journal of Climate, 1994, 7(6): 929-948.

    [8] Reynolds W, Smith T, Liu C, et al. Daily high-resolu-tion-blended analyses for sea surface temperature[J]. Journal of Climate, 2007, 20(11): 5473- 5496.

    [9] Shaw A, Vennell R. A front-following algorithm for AVHRR SST imagery[J]. Remote Sensing of Environment, 2000, 72(3): 317-327.

    [10] Zabolotskikh E, Mitnik L, Reul N, et al. New possibilities for geophysical parameter retrieval opened by GCOM-W1 AMSR2[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 8(9): 4248-4261.

    [11] Zhao Yili, Zhu Jianhua, Lin Mingsen, et al. Assessment of the initial sea surface temperature product of the scanning microwave radiometer aboard on HY-2 Satellite[J]. Acta Oceanologica Sinica, 2014, 33(1): 109- 113.

    [12] 蔣興偉, 林明森, 宋清濤. 海洋二號(hào)衛(wèi)星主被動(dòng)微波遙感探測(cè)技術(shù)研究[J]. 中國(guó)工程科學(xué), 2013, 15(7): 4-11.Jiang Xingwei, Lin Mingsen, Song Qingtao. Active and passive microwave remote sensing technology of the HY-2A ocean satellite mission[J]. Engineering Sciences, 2013, 15(7): 4-11.

    [13] 孫廣輪, 關(guān)道明, 趙冬至, 等. 星載微波遙感觀測(cè)海表溫度的研究進(jìn)展[J]. 遙感技術(shù)與應(yīng)用, 2013, 28(4): 721-730. Sun Guanglun, Guan Daoming, Zhao Dongzhi, et al. Research on observing sea surface temperature(SST) based on microwave remote sensing by satellite[J]. Remote Sensing Technology and Application, 2013, 28(4): 721-730.

    (本文編輯: 劉珊珊)

    Quality assessment and correction of SST fusion product in the Bohai Sea and the Huanghai Seas

    CAO Kai-xiang1, JIN Xi-fang2, SUN Wei-fu1, WANG Jin3

    (1. The First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China; 2. North China Sea Marine Forecasting Center of State Oceanic Administration, Qingdao 266061, China; 3. Qingdao University, Qingdao 266003, China)

    Sea surface temperature (SST) is a key parameter used in physical ocean research and has an important influence on sea-air interaction and global climate change. Space-borne microwave and infrared radiometers are the main methods used in SST remote sensing. However, although microwave radiometers provide all-weather technology their spatial resolution is sparse, and although infrared radiometers have a good resolution they are affected by cloud. In this study, a 0.1° daily SST fusion product is analyzed using the optimum interpolation method. Data resources include WindSat, AMSR-E, ASMR2, HY-2 RM, MODIS, and AVHRR. Buoy data are used to validate and correct the fusion SST product and its SST characteristics in the Bohai and the Huanghai Seas are also studied. Results show that the RMS error of this SST product is 1.1℃, but this decreases slightly after correction. SST characteristics show the lowest standard deviation of SST is in winter, and SST tends to be inhomogeneous in spring when sea temperatures rise.

    SST(sea surface temperature); data fusion; the Bohai Sea and the Huanghai Sea; quality assessment; distribution characteristics

    Aug. 31, 2016

    曹凱翔(1994-), 男, 山東諸城人, 碩士, 研究方向?yàn)槲⒉ㄝ椛溆?jì)信息提取, 電話: 0532-58750699, E-mail: 676751961@ qq.com; 王進(jìn),通信作者, 博士, 主要從事被動(dòng)微波海洋遙感方向研究, 電話: 0532-85950680, E-mail: wangjin@qdu.edu.cn

    P731.11

    A

    1000-3096(2017)09-0050-06

    10.11759/hykx20160831001

    2016-08-31;

    2016-09-20

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFA0600102); 國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)(2013AA122803); 海洋公益性行業(yè)科研專項(xiàng)(2013418032)

    [The National Key Research and Development Program of China, No. 2016YFA0600102; National High-tech R&D Program of China (863 Program), No. 2013AA122803; Public Science and Technology Research Projects of Ocean, No. 2013418032]

    猜你喜歡
    輻射計(jì)海溫黃海
    基于深度學(xué)習(xí)的海溫觀測(cè)數(shù)據(jù)質(zhì)量控制應(yīng)用研究
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    風(fēng)云四號(hào)A星多通道掃描成像輻射計(jì)第一幅彩色合成圖像
    上海航天(2021年3期)2021-07-21 14:02:12
    東方濕地 黃海明珠
    基于CLEAN算法對(duì)一維綜合孔徑輻射計(jì)成像誤差的校正
    南印度洋偶極型海溫與中國(guó)西南地區(qū)初秋降水的關(guān)系
    黃海簡(jiǎn)介
    微波輻射計(jì)對(duì)昆明霧天的監(jiān)測(cè)及應(yīng)用
    三角恒等變換去哪兒了
    用于目標(biāo)輻射特性測(cè)量的W波段輻射計(jì)
    丝瓜视频免费看黄片| 久久精品国产自在天天线| av视频免费观看在线观看| 国产亚洲精品久久久com| 久久久午夜欧美精品| 少妇精品久久久久久久| 亚洲综合精品二区| 免费观看无遮挡的男女| 午夜福利视频精品| 久久精品久久精品一区二区三区| 天堂8中文在线网| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美 国产精品| 中文字幕制服av| 欧美精品国产亚洲| 美女主播在线视频| 精品酒店卫生间| 男人狂女人下面高潮的视频| 狂野欧美激情性xxxx在线观看| 男女免费视频国产| 国产精品蜜桃在线观看| 青青草视频在线视频观看| 免费人妻精品一区二区三区视频| 18禁在线无遮挡免费观看视频| 免费av不卡在线播放| 国产亚洲一区二区精品| 亚洲欧美成人综合另类久久久| 人妻 亚洲 视频| 亚洲精品国产色婷婷电影| 日韩人妻高清精品专区| 国产精品久久久久久精品古装| 日韩中字成人| 国产精品久久久久久久久免| 18禁在线播放成人免费| 少妇人妻一区二区三区视频| 内射极品少妇av片p| kizo精华| 久久久久久久大尺度免费视频| 亚洲国产欧美人成| 午夜福利在线在线| 精品一品国产午夜福利视频| 国产美女午夜福利| 高清av免费在线| 成人毛片a级毛片在线播放| 久热这里只有精品99| 日韩亚洲欧美综合| 国产男人的电影天堂91| 国产一区二区在线观看日韩| 成人高潮视频无遮挡免费网站| 午夜激情福利司机影院| 久久ye,这里只有精品| 街头女战士在线观看网站| 日韩国内少妇激情av| 纯流量卡能插随身wifi吗| 国语对白做爰xxxⅹ性视频网站| 美女国产视频在线观看| 亚洲欧美精品自产自拍| 联通29元200g的流量卡| 亚洲国产色片| 91久久精品国产一区二区三区| 亚洲真实伦在线观看| 精品酒店卫生间| 成年女人在线观看亚洲视频| 麻豆乱淫一区二区| 国产高清有码在线观看视频| 国产伦精品一区二区三区视频9| 免费观看在线日韩| 国产伦理片在线播放av一区| 超碰97精品在线观看| 精品少妇久久久久久888优播| 久久精品久久久久久噜噜老黄| 亚洲欧美一区二区三区黑人 | 国产免费又黄又爽又色| 亚洲四区av| 夫妻午夜视频| 亚洲欧美成人综合另类久久久| 纵有疾风起免费观看全集完整版| 亚洲人成网站高清观看| 水蜜桃什么品种好| 大片免费播放器 马上看| 精品午夜福利在线看| 日韩中字成人| 国产精品一区二区三区四区免费观看| 亚洲av在线观看美女高潮| 18禁裸乳无遮挡免费网站照片| 亚洲精品乱久久久久久| 国产永久视频网站| 91精品一卡2卡3卡4卡| 国产爽快片一区二区三区| 午夜福利影视在线免费观看| 亚洲av.av天堂| 国产亚洲av片在线观看秒播厂| 亚洲自偷自拍三级| 成人综合一区亚洲| 一级毛片电影观看| 欧美高清成人免费视频www| 伦理电影免费视频| 尾随美女入室| 久久6这里有精品| 久久久久国产精品人妻一区二区| av.在线天堂| 国产亚洲精品久久久com| 免费播放大片免费观看视频在线观看| 新久久久久国产一级毛片| 久久人人爽人人爽人人片va| 日韩国内少妇激情av| 国产亚洲精品久久久com| 日韩中字成人| 男女啪啪激烈高潮av片| 美女脱内裤让男人舔精品视频| 国产在线男女| 免费看av在线观看网站| 亚洲va在线va天堂va国产| 春色校园在线视频观看| 超碰av人人做人人爽久久| 精品国产三级普通话版| 欧美日韩综合久久久久久| 人妻 亚洲 视频| 亚洲精品一区蜜桃| 青春草国产在线视频| 国产精品嫩草影院av在线观看| 精品久久久久久久末码| 精品一区在线观看国产| 一本—道久久a久久精品蜜桃钙片| 观看av在线不卡| 91久久精品电影网| 欧美日韩综合久久久久久| 国产伦在线观看视频一区| 国产亚洲av片在线观看秒播厂| 性色av一级| 国产精品久久久久成人av| 久久精品久久精品一区二区三区| 中国三级夫妇交换| 国产成人精品久久久久久| 国产在线男女| 欧美日韩亚洲高清精品| 免费观看的影片在线观看| 人人妻人人添人人爽欧美一区卜 | 麻豆国产97在线/欧美| 日韩不卡一区二区三区视频在线| 午夜日本视频在线| 国产精品.久久久| 久久婷婷青草| 国产深夜福利视频在线观看| 亚洲精品成人av观看孕妇| 久久国产乱子免费精品| 国产深夜福利视频在线观看| 免费久久久久久久精品成人欧美视频 | 人妻系列 视频| 大片免费播放器 马上看| 久久午夜福利片| 一级二级三级毛片免费看| av又黄又爽大尺度在线免费看| 少妇 在线观看| 欧美日韩一区二区视频在线观看视频在线| 熟女人妻精品中文字幕| 一级a做视频免费观看| 自拍偷自拍亚洲精品老妇| 蜜桃在线观看..| 内地一区二区视频在线| 97精品久久久久久久久久精品| 国产精品久久久久久av不卡| 日韩av免费高清视频| 日本wwww免费看| 国产免费福利视频在线观看| 亚洲经典国产精华液单| 另类亚洲欧美激情| 亚洲美女搞黄在线观看| 国产无遮挡羞羞视频在线观看| 毛片女人毛片| 深爱激情五月婷婷| 免费不卡的大黄色大毛片视频在线观看| 免费观看av网站的网址| 成人一区二区视频在线观看| 人体艺术视频欧美日本| 国产免费又黄又爽又色| 又黄又爽又刺激的免费视频.| 国产一级毛片在线| 天美传媒精品一区二区| 亚洲自偷自拍三级| 七月丁香在线播放| 国产精品秋霞免费鲁丝片| 五月伊人婷婷丁香| 国产黄频视频在线观看| 制服丝袜香蕉在线| 免费av不卡在线播放| 久久亚洲国产成人精品v| 亚洲av中文字字幕乱码综合| 国产精品麻豆人妻色哟哟久久| 久久99热6这里只有精品| 亚洲av免费高清在线观看| 卡戴珊不雅视频在线播放| 大香蕉97超碰在线| 毛片一级片免费看久久久久| 成人18禁高潮啪啪吃奶动态图 | 只有这里有精品99| 黄色日韩在线| 日韩,欧美,国产一区二区三区| av视频免费观看在线观看| 亚洲图色成人| 少妇熟女欧美另类| 王馨瑶露胸无遮挡在线观看| 亚洲精品久久午夜乱码| 超碰av人人做人人爽久久| a级毛片免费高清观看在线播放| 亚洲欧洲国产日韩| 久久99蜜桃精品久久| 久久热精品热| 午夜福利在线观看免费完整高清在| 看免费成人av毛片| 午夜免费观看性视频| 久久久久久人妻| 国产精品久久久久久精品电影小说 | 国产视频首页在线观看| 国产精品国产三级国产av玫瑰| 免费播放大片免费观看视频在线观看| 一区在线观看完整版| 免费高清在线观看视频在线观看| 日韩人妻高清精品专区| 日本-黄色视频高清免费观看| 国产精品一区二区在线观看99| 免费av不卡在线播放| 99视频精品全部免费 在线| h日本视频在线播放| 十八禁网站网址无遮挡 | 18禁动态无遮挡网站| 免费黄网站久久成人精品| 国产久久久一区二区三区| 亚洲人成网站在线播| 亚洲精品一区蜜桃| kizo精华| 蜜桃在线观看..| 不卡视频在线观看欧美| 欧美+日韩+精品| 中文资源天堂在线| 日韩成人伦理影院| 色5月婷婷丁香| 亚洲久久久国产精品| 高清视频免费观看一区二区| 亚洲精品一二三| 久久亚洲国产成人精品v| 久久久欧美国产精品| h视频一区二区三区| 午夜视频国产福利| 人妻 亚洲 视频| 欧美日本视频| 日韩一区二区视频免费看| 欧美高清成人免费视频www| 午夜视频国产福利| 日韩在线高清观看一区二区三区| 久久这里有精品视频免费| av视频免费观看在线观看| 97超碰精品成人国产| 中文资源天堂在线| 成人高潮视频无遮挡免费网站| .国产精品久久| 免费少妇av软件| 国产黄色视频一区二区在线观看| 国产成人a∨麻豆精品| 亚洲av成人精品一二三区| 国产高清不卡午夜福利| 大又大粗又爽又黄少妇毛片口| 最近手机中文字幕大全| 18禁裸乳无遮挡动漫免费视频| 亚洲精品aⅴ在线观看| 成年美女黄网站色视频大全免费 | 亚洲av免费高清在线观看| 五月玫瑰六月丁香| av网站免费在线观看视频| 人妻一区二区av| 国产乱人偷精品视频| 亚洲丝袜综合中文字幕| 久久久欧美国产精品| 国产伦理片在线播放av一区| 国产成人freesex在线| 色网站视频免费| 22中文网久久字幕| 2022亚洲国产成人精品| 国产午夜精品久久久久久一区二区三区| 久久久久久久久大av| 国产精品人妻久久久影院| 插阴视频在线观看视频| 99热这里只有精品一区| 国产久久久一区二区三区| 乱码一卡2卡4卡精品| 成人一区二区视频在线观看| 国产精品.久久久| 国产永久视频网站| 这个男人来自地球电影免费观看 | 国产精品偷伦视频观看了| 免费观看性生交大片5| 国产精品人妻久久久久久| 九草在线视频观看| 内地一区二区视频在线| 熟女av电影| 精品久久久久久久末码| 麻豆国产97在线/欧美| 日韩欧美 国产精品| 成人亚洲精品一区在线观看 | 少妇 在线观看| 一本久久精品| 亚洲欧美日韩卡通动漫| 又粗又硬又长又爽又黄的视频| 免费人妻精品一区二区三区视频| freevideosex欧美| 午夜福利在线在线| 18禁在线无遮挡免费观看视频| 久久婷婷青草| freevideosex欧美| 王馨瑶露胸无遮挡在线观看| 国产欧美亚洲国产| 欧美一区二区亚洲| 尤物成人国产欧美一区二区三区| 日产精品乱码卡一卡2卡三| 一级毛片我不卡| 精品一区二区三卡| 高清黄色对白视频在线免费看 | 啦啦啦中文免费视频观看日本| 精品亚洲成国产av| 国产女主播在线喷水免费视频网站| 91精品伊人久久大香线蕉| 少妇精品久久久久久久| 五月伊人婷婷丁香| av不卡在线播放| 久久久久久久久久久免费av| 夫妻午夜视频| 乱码一卡2卡4卡精品| 成人毛片a级毛片在线播放| 国产高清有码在线观看视频| 日韩电影二区| 亚洲美女黄色视频免费看| 久久久久精品性色| 少妇高潮的动态图| 久久青草综合色| 国产大屁股一区二区在线视频| 人妻夜夜爽99麻豆av| 成人18禁高潮啪啪吃奶动态图 | 国产亚洲午夜精品一区二区久久| 91精品国产九色| 亚洲成色77777| 欧美性感艳星| 97在线人人人人妻| 国内少妇人妻偷人精品xxx网站| av国产久精品久网站免费入址| 成年女人在线观看亚洲视频| 欧美日韩精品成人综合77777| 欧美成人一区二区免费高清观看| 免费人成在线观看视频色| 国产精品秋霞免费鲁丝片| 亚洲四区av| 蜜桃在线观看..| 黄片wwwwww| 麻豆乱淫一区二区| 亚洲成人av在线免费| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区国产| 国产在线视频一区二区| 亚洲国产欧美在线一区| 国产视频内射| 成年美女黄网站色视频大全免费 | 男人爽女人下面视频在线观看| 尾随美女入室| 国国产精品蜜臀av免费| 狠狠精品人妻久久久久久综合| 肉色欧美久久久久久久蜜桃| 青春草亚洲视频在线观看| 久久婷婷青草| 九九久久精品国产亚洲av麻豆| 色婷婷久久久亚洲欧美| 亚洲最大成人中文| 亚洲国产精品一区三区| 国产大屁股一区二区在线视频| 一区二区三区精品91| 久久青草综合色| h视频一区二区三区| 精品人妻视频免费看| 夜夜骑夜夜射夜夜干| 国产精品一区二区在线观看99| 一级a做视频免费观看| 国产一级毛片在线| 国产黄片美女视频| 久久6这里有精品| 国产成人a区在线观看| 日产精品乱码卡一卡2卡三| 大片免费播放器 马上看| 成人一区二区视频在线观看| 日韩中文字幕视频在线看片 | 一级毛片黄色毛片免费观看视频| 搡老乐熟女国产| 国产精品国产av在线观看| 成人国产麻豆网| 乱系列少妇在线播放| kizo精华| 大码成人一级视频| 99热这里只有是精品50| 草草在线视频免费看| 在线免费十八禁| www.av在线官网国产| 高清av免费在线| 如何舔出高潮| 男女啪啪激烈高潮av片| 国语对白做爰xxxⅹ性视频网站| 免费av不卡在线播放| 男女国产视频网站| 乱码一卡2卡4卡精品| 中国三级夫妇交换| 免费看av在线观看网站| 亚洲第一区二区三区不卡| 夜夜骑夜夜射夜夜干| 人妻制服诱惑在线中文字幕| 免费观看无遮挡的男女| 美女中出高潮动态图| 国产亚洲午夜精品一区二区久久| 亚洲欧美成人精品一区二区| 国产精品人妻久久久影院| 久久久久久久亚洲中文字幕| 国产成人精品婷婷| 人人妻人人看人人澡| 亚洲精品国产成人久久av| 青春草国产在线视频| 在线天堂最新版资源| 国产在视频线精品| 三级经典国产精品| 最后的刺客免费高清国语| 久久亚洲国产成人精品v| 久久影院123| 亚洲第一av免费看| 欧美精品国产亚洲| 直男gayav资源| 丰满人妻一区二区三区视频av| 色视频www国产| 中文欧美无线码| 成年人午夜在线观看视频| 亚洲av日韩在线播放| 在线 av 中文字幕| 黄色日韩在线| 91久久精品国产一区二区三区| 国产 一区 欧美 日韩| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 亚洲国产av新网站| 蜜桃亚洲精品一区二区三区| av国产免费在线观看| 亚洲最大成人中文| 18禁裸乳无遮挡动漫免费视频| 99热全是精品| 天堂中文最新版在线下载| 午夜福利视频精品| 在线播放无遮挡| 色哟哟·www| 精品人妻一区二区三区麻豆| 高清黄色对白视频在线免费看 | 国产精品伦人一区二区| 精品少妇久久久久久888优播| 国内揄拍国产精品人妻在线| 色哟哟·www| 中国国产av一级| 日韩欧美精品免费久久| 三级国产精品片| 久久久久久久久久久丰满| 高清毛片免费看| 久热这里只有精品99| 伦精品一区二区三区| 成人毛片a级毛片在线播放| 最新中文字幕久久久久| 精品久久国产蜜桃| 久久青草综合色| 人妻少妇偷人精品九色| 国产深夜福利视频在线观看| 老师上课跳d突然被开到最大视频| 国产成人freesex在线| 亚洲av二区三区四区| 国产亚洲5aaaaa淫片| 色哟哟·www| 国产精品久久久久成人av| 99久国产av精品国产电影| 日日啪夜夜爽| 内射极品少妇av片p| 亚洲av中文av极速乱| av又黄又爽大尺度在线免费看| 视频中文字幕在线观看| 国产 精品1| 少妇精品久久久久久久| 网址你懂的国产日韩在线| 99热国产这里只有精品6| 美女xxoo啪啪120秒动态图| 狂野欧美白嫩少妇大欣赏| 国产色婷婷99| 欧美日韩一区二区视频在线观看视频在线| 一区在线观看完整版| 五月玫瑰六月丁香| 我要看日韩黄色一级片| 亚洲成人手机| 黄色视频在线播放观看不卡| 精品99又大又爽又粗少妇毛片| 欧美人与善性xxx| 自拍偷自拍亚洲精品老妇| 噜噜噜噜噜久久久久久91| 夜夜骑夜夜射夜夜干| 国产淫片久久久久久久久| 成年免费大片在线观看| 极品少妇高潮喷水抽搐| 亚洲最大成人中文| 交换朋友夫妻互换小说| 只有这里有精品99| 一本—道久久a久久精品蜜桃钙片| 国产av国产精品国产| av女优亚洲男人天堂| 国产精品麻豆人妻色哟哟久久| 人人妻人人澡人人爽人人夜夜| 亚洲精品亚洲一区二区| av免费在线看不卡| 一区二区三区免费毛片| 99久久综合免费| 18禁动态无遮挡网站| 国产探花极品一区二区| 久久精品人妻少妇| 日韩大片免费观看网站| 久久97久久精品| 久久99热6这里只有精品| 国产免费又黄又爽又色| 在线精品无人区一区二区三 | 欧美 日韩 精品 国产| 午夜视频国产福利| 亚洲高清免费不卡视频| 一本色道久久久久久精品综合| 蜜臀久久99精品久久宅男| 大话2 男鬼变身卡| 多毛熟女@视频| 人妻系列 视频| av福利片在线观看| 高清av免费在线| 成人二区视频| 亚洲精品aⅴ在线观看| 丰满乱子伦码专区| 超碰97精品在线观看| 亚洲av二区三区四区| 亚洲国产精品一区三区| 久久久久性生活片| 亚洲在久久综合| 日本av免费视频播放| 热99国产精品久久久久久7| tube8黄色片| 熟女电影av网| 这个男人来自地球电影免费观看 | 小蜜桃在线观看免费完整版高清| 我的老师免费观看完整版| 久久久久久久久大av| videossex国产| 中文在线观看免费www的网站| 欧美最新免费一区二区三区| av免费观看日本| 国产免费一区二区三区四区乱码| 国产91av在线免费观看| 熟女电影av网| 黑人猛操日本美女一级片| 欧美xxxx性猛交bbbb| 久久久精品免费免费高清| 国产黄色免费在线视频| 在现免费观看毛片| 直男gayav资源| 亚洲天堂av无毛| 最后的刺客免费高清国语| 欧美激情极品国产一区二区三区 | 免费av中文字幕在线| 成年av动漫网址| .国产精品久久| 最近最新中文字幕免费大全7| 亚洲精品乱久久久久久| 毛片女人毛片| 国产av一区二区精品久久 | 视频中文字幕在线观看| 在线免费十八禁| 亚洲精品一二三| 高清av免费在线| 成年免费大片在线观看| 国产男女超爽视频在线观看| 国产精品精品国产色婷婷| 麻豆乱淫一区二区| 国产精品一区二区三区四区免费观看| 91午夜精品亚洲一区二区三区| 亚洲va在线va天堂va国产| 国产成人免费观看mmmm| 性色avwww在线观看| 伦精品一区二区三区| 亚洲精品日韩在线中文字幕| 亚州av有码| 国产淫片久久久久久久久| 色网站视频免费| 美女国产视频在线观看| 亚洲av电影在线观看一区二区三区| 亚洲成人一二三区av| 国产淫语在线视频| 日本-黄色视频高清免费观看| 久久99热这里只频精品6学生| 九九爱精品视频在线观看| 寂寞人妻少妇视频99o| 丰满迷人的少妇在线观看| 在线观看人妻少妇| 蜜臀久久99精品久久宅男| 日日摸夜夜添夜夜爱| 国产色爽女视频免费观看| 免费观看a级毛片全部| 我要看黄色一级片免费的| 国产av一区二区精品久久 | 亚洲精品久久久久久婷婷小说| 欧美zozozo另类| 一个人免费看片子| 一级毛片 在线播放| 亚洲av综合色区一区| 人人妻人人澡人人爽人人夜夜| 欧美另类一区| 精品亚洲成a人片在线观看 | 搡女人真爽免费视频火全软件| 久久午夜福利片|