• <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ì)
    免费大片18禁| 下体分泌物呈黄色| 妹子高潮喷水视频| 久久久国产欧美日韩av| 欧美激情 高清一区二区三区| 大香蕉久久网| 免费看光身美女| 不卡视频在线观看欧美| 国产免费福利视频在线观看| 久久99精品国语久久久| 一本色道久久久久久精品综合| 亚洲欧洲日产国产| 精品少妇久久久久久888优播| 午夜老司机福利剧场| 久久久久精品性色| 超碰97精品在线观看| 久久精品国产亚洲av天美| tube8黄色片| 欧美变态另类bdsm刘玥| 黑丝袜美女国产一区| 我的老师免费观看完整版| 日本与韩国留学比较| 亚洲av成人精品一二三区| 又粗又硬又长又爽又黄的视频| 亚洲欧洲精品一区二区精品久久久 | 成人亚洲欧美一区二区av| 国产女主播在线喷水免费视频网站| 伊人亚洲综合成人网| 亚洲av欧美aⅴ国产| 成人毛片a级毛片在线播放| 视频在线观看一区二区三区| 亚洲精品,欧美精品| 又大又黄又爽视频免费| 毛片一级片免费看久久久久| 免费观看的影片在线观看| av福利片在线| 大片电影免费在线观看免费| 狂野欧美激情性xxxx在线观看| 国产亚洲最大av| 午夜激情久久久久久久| 成人影院久久| 天美传媒精品一区二区| 好男人视频免费观看在线| 午夜福利视频精品| 久久精品夜色国产| 乱人伦中国视频| 亚洲成人一二三区av| av线在线观看网站| 欧美人与善性xxx| 成年美女黄网站色视频大全免费 | xxxhd国产人妻xxx| 一本色道久久久久久精品综合| 精品久久久久久久久av| 黄色毛片三级朝国网站| 在线 av 中文字幕| 建设人人有责人人尽责人人享有的| 黑人欧美特级aaaaaa片| 午夜免费鲁丝| 啦啦啦视频在线资源免费观看| 亚洲av电影在线观看一区二区三区| 亚洲精品自拍成人| 亚洲人成网站在线观看播放| 丰满少妇做爰视频| 精品久久久久久久久亚洲| 亚洲精品久久久久久婷婷小说| av网站免费在线观看视频| 国产熟女欧美一区二区| 成人免费观看视频高清| 人人妻人人爽人人添夜夜欢视频| 亚洲国产日韩一区二区| 插逼视频在线观看| 只有这里有精品99| 亚洲精品日韩av片在线观看| 五月玫瑰六月丁香| 日本爱情动作片www.在线观看| 999精品在线视频| 在线精品无人区一区二区三| 一区二区av电影网| 亚洲av电影在线观看一区二区三区| 大香蕉97超碰在线| 熟女人妻精品中文字幕| 2018国产大陆天天弄谢| av不卡在线播放| 丝袜喷水一区| 欧美日韩亚洲高清精品| 80岁老熟妇乱子伦牲交| 99视频精品全部免费 在线| 一区在线观看完整版| 99九九线精品视频在线观看视频| 97精品久久久久久久久久精品| 日本欧美视频一区| 看十八女毛片水多多多| av在线观看视频网站免费| 国产精品国产三级国产av玫瑰| 97超视频在线观看视频| 亚洲欧洲精品一区二区精品久久久 | 人人妻人人爽人人添夜夜欢视频| 国产精品久久久久久久电影| 欧美日韩视频精品一区| 一本—道久久a久久精品蜜桃钙片| 一级毛片电影观看| 久久国内精品自在自线图片| 男人爽女人下面视频在线观看| 久久久久久久大尺度免费视频| 两个人免费观看高清视频| 插逼视频在线观看| 久久久久久久久大av| 大片免费播放器 马上看| 熟妇人妻不卡中文字幕| 精品少妇黑人巨大在线播放| 国产永久视频网站| 日日爽夜夜爽网站| 欧美97在线视频| 一级,二级,三级黄色视频| 国产高清国产精品国产三级| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 成人影院久久| 纯流量卡能插随身wifi吗| 人体艺术视频欧美日本| 亚洲欧洲精品一区二区精品久久久 | 尾随美女入室| 亚洲国产av影院在线观看| 一级毛片我不卡| 午夜影院在线不卡| 91国产中文字幕| 少妇猛男粗大的猛烈进出视频| 丁香六月天网| freevideosex欧美| 美女福利国产在线| 91成人精品电影| 永久免费av网站大全| 多毛熟女@视频| 性色av一级| 日韩熟女老妇一区二区性免费视频| 女人精品久久久久毛片| 下体分泌物呈黄色| 少妇高潮的动态图| 国产亚洲午夜精品一区二区久久| 亚洲精品国产av蜜桃| 精品国产露脸久久av麻豆| 精品久久蜜臀av无| 久热久热在线精品观看| 国产高清不卡午夜福利| 啦啦啦在线观看免费高清www| 99热6这里只有精品| 黄色配什么色好看| 免费不卡的大黄色大毛片视频在线观看| 能在线免费看毛片的网站| 欧美一级a爱片免费观看看| 亚洲精品乱码久久久久久按摩| 18禁在线播放成人免费| 久久久a久久爽久久v久久| 亚洲av欧美aⅴ国产| 久久综合国产亚洲精品| 久久久国产欧美日韩av| 国产一区二区三区av在线| 美女福利国产在线| 久久狼人影院| 蜜臀久久99精品久久宅男| 国精品久久久久久国模美| 老女人水多毛片| 日韩一本色道免费dvd| 欧美性感艳星| 新久久久久国产一级毛片| 国产永久视频网站| 午夜91福利影院| 人妻 亚洲 视频| 国产无遮挡羞羞视频在线观看| 亚洲性久久影院| 日韩不卡一区二区三区视频在线| 狂野欧美白嫩少妇大欣赏| 国产成人freesex在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美中文字幕日韩二区| 欧美激情极品国产一区二区三区 | 国产免费一区二区三区四区乱码| 夜夜爽夜夜爽视频| 最新中文字幕久久久久| 国产黄色视频一区二区在线观看| 成人免费观看视频高清| 一级毛片电影观看| 久久久久久人妻| 大话2 男鬼变身卡| 精品少妇久久久久久888优播| 亚洲少妇的诱惑av| 久久久午夜欧美精品| 国产成人免费无遮挡视频| 久久国产精品男人的天堂亚洲 | 少妇的逼好多水| 国产日韩欧美视频二区| 精品国产露脸久久av麻豆| 欧美另类一区| 新久久久久国产一级毛片| 婷婷色麻豆天堂久久| 国产精品久久久久久av不卡| 久久av网站| 一级毛片 在线播放| 久久久久国产精品人妻一区二区| 国产av码专区亚洲av| 成人影院久久| 中文字幕久久专区| 国产成人91sexporn| 国产黄色视频一区二区在线观看| 日本免费在线观看一区| 色视频在线一区二区三区| 国产亚洲av片在线观看秒播厂| 国产精品国产三级专区第一集| 亚洲成人手机| 亚洲,欧美,日韩| 一本色道久久久久久精品综合| 欧美丝袜亚洲另类| 美女大奶头黄色视频| 男男h啪啪无遮挡| 午夜免费观看性视频| 亚洲av成人精品一二三区| 久久久久国产精品人妻一区二区| 免费观看在线日韩| 毛片一级片免费看久久久久| 亚洲在久久综合| 日韩中文字幕视频在线看片| 高清毛片免费看| 男女啪啪激烈高潮av片| 老司机亚洲免费影院| 不卡视频在线观看欧美| 搡老乐熟女国产| av在线app专区| 蜜臀久久99精品久久宅男| 一本久久精品| videosex国产| 久久久久久久精品精品| 亚洲人成网站在线观看播放| 亚洲av福利一区| 久久久久久久大尺度免费视频| 欧美精品人与动牲交sv欧美| 丰满少妇做爰视频| 大香蕉久久网| 中国国产av一级| 在线观看一区二区三区激情| 亚洲av成人精品一二三区| 99热6这里只有精品| 亚洲人成77777在线视频| 国产视频首页在线观看| 青春草视频在线免费观看| 欧美bdsm另类| 丰满饥渴人妻一区二区三| 中文字幕人妻丝袜制服| 久久 成人 亚洲| 赤兔流量卡办理| 国产精品国产三级国产av玫瑰| 飞空精品影院首页| 久久99蜜桃精品久久| 99久国产av精品国产电影| 国产成人freesex在线| 男女边吃奶边做爰视频| 亚洲第一av免费看| 少妇人妻久久综合中文| 久久久久久久大尺度免费视频| 男女边摸边吃奶| 国产免费一区二区三区四区乱码| 国产精品一区二区三区四区免费观看| 亚洲情色 制服丝袜| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 国产精品久久久久久精品古装| 在线看a的网站| 日日摸夜夜添夜夜爱| 赤兔流量卡办理| 一本色道久久久久久精品综合| 18禁动态无遮挡网站| 婷婷色麻豆天堂久久| 一本一本综合久久| 国产精品久久久久久久电影| 久久97久久精品| 考比视频在线观看| 男女边吃奶边做爰视频| 亚洲第一av免费看| 黑人高潮一二区| 成人二区视频| 如日韩欧美国产精品一区二区三区 | 老熟女久久久| 欧美bdsm另类| 亚洲欧美日韩另类电影网站| 美女大奶头黄色视频| 日韩免费高清中文字幕av| www.av在线官网国产| 亚洲国产精品一区三区| 久久久精品免费免费高清| 成人国产av品久久久| 欧美bdsm另类| 欧美 日韩 精品 国产| 最新的欧美精品一区二区| 国产亚洲精品久久久com| 久久久久久久久久久免费av| 日本爱情动作片www.在线观看| 久久久欧美国产精品| 各种免费的搞黄视频| 三上悠亚av全集在线观看| 超色免费av| 中文乱码字字幕精品一区二区三区| 久久久久久久久久人人人人人人| 男男h啪啪无遮挡| 欧美激情国产日韩精品一区| 9色porny在线观看| 久久亚洲国产成人精品v| 一级片'在线观看视频| 久久久国产欧美日韩av| 中国三级夫妇交换| 91成人精品电影| 视频区图区小说| 亚洲成人一二三区av| 这个男人来自地球电影免费观看 | 日本午夜av视频| a级毛色黄片| 制服丝袜香蕉在线| 国产精品嫩草影院av在线观看| 成人亚洲欧美一区二区av| 日日摸夜夜添夜夜爱| 99re6热这里在线精品视频| 激情五月婷婷亚洲| 国产不卡av网站在线观看| 一级片'在线观看视频| tube8黄色片| 国产成人免费观看mmmm| 国产乱人偷精品视频| 街头女战士在线观看网站| 久久精品夜色国产| 久久久午夜欧美精品| 熟女人妻精品中文字幕| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 人妻人人澡人人爽人人| 亚洲人成77777在线视频| 日日摸夜夜添夜夜爱| 欧美日本中文国产一区发布| 国产高清有码在线观看视频| 精品人妻在线不人妻| xxx大片免费视频| 精品亚洲成国产av| av又黄又爽大尺度在线免费看| 九九爱精品视频在线观看| 欧美日韩视频精品一区| 18禁动态无遮挡网站| 成人黄色视频免费在线看| 国产日韩欧美视频二区| 街头女战士在线观看网站| a级毛色黄片| 亚洲av国产av综合av卡| 人妻系列 视频| 一级,二级,三级黄色视频| 丝袜脚勾引网站| 精品国产一区二区三区久久久樱花| 精品久久久久久久久av| 色婷婷久久久亚洲欧美| 啦啦啦中文免费视频观看日本| 18禁动态无遮挡网站| 大话2 男鬼变身卡| 亚洲熟女精品中文字幕| 性色avwww在线观看| 国产精品人妻久久久久久| 美女主播在线视频| av卡一久久| av在线app专区| 亚洲av国产av综合av卡| 最后的刺客免费高清国语| 中国美白少妇内射xxxbb| 日韩制服骚丝袜av| 国产日韩欧美视频二区| 精品国产露脸久久av麻豆| 免费大片18禁| 精品一区二区三卡| 久久精品久久久久久噜噜老黄| 简卡轻食公司| 成人亚洲精品一区在线观看| 男女边吃奶边做爰视频| 亚洲精品乱码久久久v下载方式| 天天躁夜夜躁狠狠久久av| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 18禁在线无遮挡免费观看视频| 欧美最新免费一区二区三区| 高清欧美精品videossex| 日本-黄色视频高清免费观看| 中文字幕精品免费在线观看视频 | av在线观看视频网站免费| 99热这里只有是精品在线观看| 午夜日本视频在线| 天堂俺去俺来也www色官网| 欧美精品一区二区免费开放| 美女视频免费永久观看网站| 日韩精品有码人妻一区| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 一级毛片aaaaaa免费看小| 国产亚洲午夜精品一区二区久久| 熟女av电影| 天堂8中文在线网| 人人妻人人爽人人添夜夜欢视频| 一级片'在线观看视频| 婷婷色av中文字幕| 不卡视频在线观看欧美| 热99久久久久精品小说推荐| 久久久亚洲精品成人影院| 精品国产乱码久久久久久小说| a级毛片黄视频| 亚洲经典国产精华液单| 极品人妻少妇av视频| 99热6这里只有精品| 精品熟女少妇av免费看| 国产免费视频播放在线视频| 国产成人午夜福利电影在线观看| 日日撸夜夜添| 美女视频免费永久观看网站| 国产永久视频网站| 国产免费一区二区三区四区乱码| 麻豆成人av视频| 2018国产大陆天天弄谢| 999精品在线视频| 日韩人妻高清精品专区| 久久精品夜色国产| 国产午夜精品久久久久久一区二区三区| 香蕉精品网在线| 亚洲欧洲国产日韩| 一个人免费看片子| 色哟哟·www| 久久亚洲国产成人精品v| 最近中文字幕高清免费大全6| 国产欧美亚洲国产| 成人无遮挡网站| 高清午夜精品一区二区三区| av女优亚洲男人天堂| 成人毛片a级毛片在线播放| 久久人人爽人人片av| 亚洲色图 男人天堂 中文字幕 | 亚洲精品456在线播放app| av免费观看日本| 熟妇人妻不卡中文字幕| 精品午夜福利在线看| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 黄色配什么色好看| 99久久精品国产国产毛片| 久久热精品热| 亚洲精品中文字幕在线视频| 亚洲av成人精品一二三区| 又黄又爽又刺激的免费视频.| 国产亚洲午夜精品一区二区久久| 久久精品国产亚洲av涩爱| 美女内射精品一级片tv| 你懂的网址亚洲精品在线观看| 黄色视频在线播放观看不卡| 欧美精品人与动牲交sv欧美| 秋霞在线观看毛片| 成人黄色视频免费在线看| 亚洲av在线观看美女高潮| xxx大片免费视频| 18+在线观看网站| 国产av国产精品国产| 99精国产麻豆久久婷婷| 国产精品嫩草影院av在线观看| 婷婷色麻豆天堂久久| 18禁动态无遮挡网站| 成人国产av品久久久| 黑人猛操日本美女一级片| 国产日韩欧美在线精品| 久久久精品94久久精品| 99re6热这里在线精品视频| 精品卡一卡二卡四卡免费| 人妻夜夜爽99麻豆av| 麻豆乱淫一区二区| 成人国产麻豆网| 亚洲国产毛片av蜜桃av| 狂野欧美激情性xxxx在线观看| 九草在线视频观看| 性色avwww在线观看| 天天躁夜夜躁狠狠久久av| 婷婷色av中文字幕| 国产成人免费观看mmmm| 赤兔流量卡办理| 亚洲,一卡二卡三卡| 高清午夜精品一区二区三区| 在线免费观看不下载黄p国产| 日本午夜av视频| 久久久国产欧美日韩av| 成人影院久久| 久久女婷五月综合色啪小说| 简卡轻食公司| 婷婷色综合大香蕉| 亚洲国产毛片av蜜桃av| 亚洲av电影在线观看一区二区三区| 国产白丝娇喘喷水9色精品| 高清欧美精品videossex| 亚洲精品视频女| 亚洲国产成人一精品久久久| 777米奇影视久久| 色5月婷婷丁香| 能在线免费看毛片的网站| 男女边摸边吃奶| 亚洲人成网站在线播| 亚洲精品色激情综合| 国产伦理片在线播放av一区| 午夜老司机福利剧场| 亚洲一区二区三区欧美精品| 中文字幕免费在线视频6| av又黄又爽大尺度在线免费看| 亚洲国产精品999| av国产久精品久网站免费入址| 赤兔流量卡办理| 免费观看性生交大片5| 99视频精品全部免费 在线| 少妇人妻精品综合一区二区| 免费播放大片免费观看视频在线观看| 成人综合一区亚洲| 久久青草综合色| 高清午夜精品一区二区三区| 国产视频内射| 亚洲国产精品专区欧美| 久久毛片免费看一区二区三区| 美女福利国产在线| a级毛片在线看网站| 国产精品蜜桃在线观看| 精品午夜福利在线看| 久久久久久久久久成人| 中文字幕精品免费在线观看视频 | 黄色怎么调成土黄色| 全区人妻精品视频| 国产日韩欧美视频二区| 七月丁香在线播放| 国产不卡av网站在线观看| 精品国产国语对白av| 男女啪啪激烈高潮av片| 99久久人妻综合| 国产亚洲最大av| 成人亚洲欧美一区二区av| 少妇熟女欧美另类| 菩萨蛮人人尽说江南好唐韦庄| 赤兔流量卡办理| 精品卡一卡二卡四卡免费| 亚洲精品日韩在线中文字幕| av在线app专区| 亚洲av电影在线观看一区二区三区| 99国产综合亚洲精品| 国产乱来视频区| 亚洲av男天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩伦理黄色片| 大又大粗又爽又黄少妇毛片口| xxxhd国产人妻xxx| 欧美精品一区二区大全| 丝瓜视频免费看黄片| av一本久久久久| 成人亚洲欧美一区二区av| 建设人人有责人人尽责人人享有的| 成人国语在线视频| 免费看光身美女| 简卡轻食公司| 极品人妻少妇av视频| 日韩成人伦理影院| 国产精品熟女久久久久浪| 精品久久久久久久久亚洲| 国产淫语在线视频| 黄色毛片三级朝国网站| 亚洲精品乱码久久久久久按摩| 高清在线视频一区二区三区| 中文字幕最新亚洲高清| 欧美三级亚洲精品| 香蕉精品网在线| √禁漫天堂资源中文www| 国产精品.久久久| 99热6这里只有精品| 亚洲一级一片aⅴ在线观看| 伊人久久精品亚洲午夜| 精品人妻熟女av久视频| 99热网站在线观看| 亚洲无线观看免费| 久久久久网色| 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线观看播放| 亚洲综合精品二区| 国产一区二区三区av在线| 久久婷婷青草| 能在线免费看毛片的网站| 国产精品一国产av| 日韩一区二区视频免费看| 丝袜美足系列| 久久99蜜桃精品久久| 大又大粗又爽又黄少妇毛片口| 国精品久久久久久国模美| 视频在线观看一区二区三区| 欧美成人午夜免费资源| 欧美性感艳星| 人成视频在线观看免费观看| 丝瓜视频免费看黄片| 美女xxoo啪啪120秒动态图| 寂寞人妻少妇视频99o| 九九爱精品视频在线观看| 欧美性感艳星| 成人毛片a级毛片在线播放| 亚洲不卡免费看| 亚洲国产精品一区三区| 两个人的视频大全免费| 99久久人妻综合| 色视频在线一区二区三区| 亚洲av欧美aⅴ国产| 极品少妇高潮喷水抽搐| 午夜免费鲁丝| 免费观看在线日韩| 在线免费观看不下载黄p国产| 69精品国产乱码久久久| videos熟女内射| 精品久久国产蜜桃| 亚洲精华国产精华液的使用体验| 国产精品女同一区二区软件| 国产av一区二区精品久久| 亚洲国产精品一区三区|