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

    浙南沿岸海表面溫度分布特征

    2011-12-28 08:17:40韓小燕潘曉東馬林芳
    海洋通報(bào) 2011年6期
    關(guān)鍵詞:浙南時(shí)間尺度測站

    韓小燕,潘曉東,馬林芳

    (溫州海洋環(huán)境監(jiān)測中心站,浙江 溫州 325013)

    浙南沿岸海表面溫度分布特征

    韓小燕,潘曉東,馬林芳

    (溫州海洋環(huán)境監(jiān)測中心站,浙江 溫州 325013)

    文章根據(jù)浙南沿岸3個(gè)海洋觀測站大陳、坎門、南麂站51a(1960-2010年)SST觀測資料,進(jìn)行統(tǒng)計(jì)分析,得出SST時(shí)空分布特征。結(jié)果表明,浙南沿岸的SST總體呈現(xiàn)北低南高的趨勢,冬季南北差異最小,夏季南北差異最大;51 a年浙南沿岸SST總體呈現(xiàn)上升趨勢,上升速率為0.017 ℃/a;SST變化存在3.5 a、5~6 a、2~7 a及11 a左右的振蕩周期;SST主要由年循環(huán)分量決定,諧波擬合效果良好。

    浙南沿岸;SST;均一性;功率譜;諧波

    1 概 述

    海溫,特別是表層海溫(SST)是物理海洋的重要要素,是進(jìn)入國際氣象觀測網(wǎng)定時(shí)觀測的少有的海洋要素之一,也是對物理海洋研究以及全球氣候異常研究最有參考價(jià)值的海洋要素,與厄爾尼諾、拉尼娜現(xiàn)象、海平面上升現(xiàn)象都有著密切聯(lián)系;對海洋開發(fā)而言,SST的觀測和研究直接服務(wù)于水產(chǎn)養(yǎng)殖和海洋漁業(yè),服務(wù)于赤潮和海冰災(zāi)害防災(zāi)減災(zāi)業(yè)務(wù)。浙南海域東接西北太平洋,經(jīng)常遭受臺風(fēng)、赤潮、海水入侵等海洋災(zāi)害的影響,尤其隨著海洋開發(fā)活動(dòng)的日趨頻繁,影響愈發(fā)嚴(yán)重。因此,研究浙南沿岸海表溫度(SST)的變化規(guī)律十分必要,對于防災(zāi)減災(zāi)、更好地促進(jìn)和保障海洋漁業(yè)的發(fā)展具有現(xiàn)實(shí)意義。

    2 資料收集及預(yù)處理

    為較全面地反映浙南沿岸SST的總體變化特征,自北向南收集了大陳、坎門、南麂3個(gè)海洋站的SST觀測資料進(jìn)行研究,資料序列長度大陳站46 a(1960-2001年、2007-2010年),坎門站50 a(1960-2001年、2003-2010年),南麂站46 a(1960-2001年、2007-2010年),資料選用月、年平均值。所有觀測資料均按照《海濱觀測規(guī)范》進(jìn)行采集與整理,測站位置見表1。

    由于資料年限較長,期間經(jīng)歷了站點(diǎn)變遷、觀測儀器更換、觀測規(guī)范改版等,觀測數(shù)據(jù)存在不連續(xù)的可能。本文采用U檢驗(yàn)法[1,2]對觀測資料進(jìn)行歸一性檢驗(yàn),以此對SST數(shù)據(jù)有無顯著差異進(jìn)行簡單判別。

    表 1 測站經(jīng)緯度Tab.1 Latitude and longitude of stations

    SST資料收集過程中,有個(gè)別年份的數(shù)據(jù)由于儀器原因缺測,考慮到SST變化的連續(xù)性,采用三次樣條插值的方法進(jìn)行數(shù)據(jù)訂正,該方法不僅收斂性、穩(wěn)定性強(qiáng),而且方法簡單實(shí)用,計(jì)算量小,其方法為;假設(shè)在區(qū)間[a,b]上取n+1個(gè)節(jié)點(diǎn),三次樣條函數(shù)s(x)存在,mi為s(x)在點(diǎn)xi處的微商值,則每個(gè)小區(qū)間上

    同時(shí)由于坎門站資料序列相對較長,考慮利用大陳、南麂兩站與坎門站SST的相關(guān)性對缺測數(shù)據(jù)進(jìn)行訂正。采用相關(guān)系數(shù)的取值在-1.0~+1.0之間。當(dāng)r越接近+1.0時(shí),表明兩序列正相關(guān)越顯著;反之,r越接近-1.0,兩序列負(fù)相關(guān)越顯著。計(jì)算得R南麂-坎門=0.89,R大陳-坎門=0.95,相關(guān)性顯著。然后剔除重大天氣系統(tǒng)影響下的SST,根據(jù)南麂和大陳2站與坎門站逐時(shí)、逐日、逐月及逐年的SST差值對上述插值結(jié)果進(jìn)行比對訂正,最終得到完整的SST時(shí)間序列值。

    3 浙南沿岸SST分布特征

    3.1 SST空間分布特征

    浙南地區(qū)瀕臨東海。東海是西北太平洋的邊緣海,SST的分布特征主要取決于太陽輻射、暖流、海—?dú)鉄峤粨Q及入海徑流、降水等因素的影響。此外,海區(qū)地理環(huán)境、海岸形態(tài)等因素的影響,使得浙南沿岸SST空間分布較為復(fù)雜。

    年平均SST空間分布3站總體呈南高北低;南麂站最高,為19.0 ℃;坎門站次之,為18.5 ℃;大陳站最低,為17.9 ℃。SST空間分布不僅體現(xiàn)在年平均變化,季節(jié)變化也存在差異。根據(jù)海洋季節(jié)劃分[3],冬季為1-3月,春季為4-6月,夏季為7-9月,秋季為10-12月,以相應(yīng)月份的SST累年平均分析浙南沿岸SST空間分布特征(表2)。

    表 2 SST累年季平均、年平均統(tǒng)計(jì)(溫度/℃)Tab.2 Statistics of seasonal and annual average

    冬季太陽輻射為一年中最弱,主要受到北方陸地冷空氣南下影響,同時(shí)還受到黑潮與臺灣暖流的影響[4],為浙南沿岸SST最低的季節(jié)。與年平均SST變化規(guī)律略有不同,南麂站最高為10.6 ℃,坎門站SST最低為9.4 ℃;冬季平均SST各站相差最大為1.2 ℃,為一年中浙南地區(qū)南北差異最小的季節(jié)。

    春季是冬季向夏季過渡的季節(jié),隨著海域東北季風(fēng)的減弱和太陽輻射的加強(qiáng),浙南沿岸SST逐步升高。坎門站由于受基巖、水動(dòng)力環(huán)境的影響升溫較快,升幅為9.9 ℃;大陳、南麂兩站升溫較慢,且基本持平,其中大陳站8.4℃,南麂站8.5℃;春季平均SST坎門站最高為19.3 ℃,大陳站最低為18.0 ℃;各站相差最大為1.3 ℃。

    夏季太陽輻射為一年中最強(qiáng),強(qiáng)烈的太陽輻射使得浙南沿岸水溫進(jìn)一步大幅上升;坎門站升溫8.4 ℃,大陳、南麂兩站升溫基本持平,大陳站7.9℃,南麂站8.0℃;空間分布與春季相似,坎門站最高為27.7 ℃,大陳站最低為25.9 ℃;各站相差最大為1.8 ℃,為一年中浙南地區(qū)南北差異最大的季節(jié)。

    秋季是海域溫度由夏季型分布向冬季型分布過渡的季節(jié),夏季所存在的某些水文學(xué)現(xiàn)象基本消失[5],同時(shí)隨著太陽輻射的逐漸減弱,冷空氣影響的逐漸增強(qiáng),浙南沿岸水溫呈現(xiàn)下降趨勢;坎門站降溫幅度最大為10.4 ℃,大陳站降溫最小為7.8 ℃;空間分布與冬季相似,南麂站最高為19.0 ℃,坎門站最低為17.3 ℃;各站相差最大為1.7 ℃。

    3.2 SST年平均變化特征

    浙南沿岸年平均SST長期變化不僅受到高空氣候系統(tǒng)(如西太副高和東亞大槽的面積及強(qiáng)度)的影響,還受到厄爾尼諾、拉尼娜現(xiàn)象等因素的作用。通過對浙南地區(qū)3個(gè)測站年平均SST的趨勢分析可以看出(圖1),50 a年SST總體呈現(xiàn)上升趨勢,上升速率為0.017 ℃/a。其中南麂站上升速率為0.014 ℃/a,坎門、大陳兩站上升速率相近,分別為0.017 ℃/a 和0.019 ℃/a(表3)。

    浙南地區(qū)年平均SST除了明顯的長期趨勢變化,還具有階段性變化特征。文章采用6階多項(xiàng)式進(jìn)行趨勢擬合,擬合結(jié)果見圖1。具體表現(xiàn)為:3個(gè)測站變化趨勢基本接近,20世紀(jì)60年代中后期之前,SST呈現(xiàn)下降現(xiàn)象;60年代中后期到70年代中期,SST呈現(xiàn)緩慢上升;70年代中期到80年代中期,SST開始又一輪下降過程;之后到2005年,SST一直呈現(xiàn)上升,且上升幅度明顯比前個(gè)上升過程大;2005年前后至2010年,呈現(xiàn)下降。

    不同的是,坎門、南麂2站上世紀(jì)60年代中后期到80年代中期SST的變化趨勢明顯比大陳站緩慢;且坎門站80年代中期開始的SST上升事件中,擬合線曲率明顯比大陳、南麂兩站小,升溫緩慢,導(dǎo)致隨后的下降趨勢略有延遲。經(jīng)過對3個(gè)測站地理位置、觀測儀器使用情況的了解,產(chǎn)生上述現(xiàn)象的原因估計(jì)與各站不同觀測儀器使用的時(shí)間差、儀器布放海域的開闊度等因素有關(guān)。

    表 3 年均SST變化速率統(tǒng)計(jì)(速率/℃·a-1)Tab.3 Statistics of annual mean SST variation rate

    圖 1 年平均SST變化曲線Fig.1 Variation curve of annual mean SST

    圖2為浙南沿岸年平均SST距平變化曲線,可以看出上世紀(jì)80年代末期之前,SST變化基本處于常態(tài)(△T≤∣±0.2 ℃∣)或偏冷(△T≤-0.5 ℃)、略偏冷(-0.5 ℃<△T<-0.2 ℃)狀態(tài),浙南沿岸各測站年平均SST距平概率統(tǒng)計(jì)見表4。

    80年代后期開始,浙南沿岸SST距平逐步上升,偏暖(△T≥0.5 ℃)年份明顯增加,2000-2010年間,南麂、坎門、大陳3站偏暖年數(shù)分別為6 a、8 a和7 a。通過分析3站年平均SST及年平均氣溫之間相互關(guān)系,發(fā)現(xiàn)兩者變化趨勢非常接近,且相關(guān)性很高,相關(guān)系數(shù)分別為0.77(大陳)、0.92(坎門)、0.81(南麂),均通過顯著性檢驗(yàn)。因此一般情況下,人們常把氣候變暖作為引起SST長序列線性升高的主要因素。

    3.3 年平均SST變化周期分析

    為具體了解浙南沿岸SST變化的周期分布特征,文章采用功率譜分析方法[6]對SST的主要振動(dòng)周期進(jìn)行分析。首先,考慮到SST序列中包含多種時(shí)間尺度變化,采用9階遞歸濾波器對SST觀測數(shù)據(jù)進(jìn)行5年濾波處理;計(jì)算自相關(guān)系數(shù)R(i)(i=0,1,2,…,m),考慮到功率譜估計(jì)的真實(shí)性,最大滯后長度m取51/3;采用Hanning系數(shù)對得到的功率譜估計(jì)進(jìn)行平滑處理,得到連續(xù)的譜值;對譜估計(jì)進(jìn)行顯著性檢驗(yàn),由于自相關(guān)系數(shù)R(1)為較大正值,表明譜估計(jì)序列具有持續(xù)性,因此采用紅噪聲標(biāo)準(zhǔn)譜進(jìn)行顯著性檢驗(yàn),根據(jù)顯著性水平0.05的紅噪聲標(biāo)準(zhǔn)譜限值檢驗(yàn)得出譜估計(jì)的顯著周期(表5)。其中,3個(gè)測站SST同時(shí)具有15 a 、6 a、5 a、4.3 a和3.8 a時(shí)間尺度的變化周期。除此以外,大陳、南麂2個(gè)測站具有30 a時(shí)間尺度的變化,坎門、南麂2個(gè)測站具有7.5 a時(shí)間尺度的變化,南麂站具有10 a時(shí)間尺度的變化。

    圖 2 年平均SST距平變化曲線Fig.2 Anomaly variation curve of annual mean SST

    表 4 年平均SST距平概率統(tǒng)計(jì)Tab.4 Anomaly probability statistics of annual mean SST

    可以看出,浙南沿岸3個(gè)測站的周期變化與郭偉其等[7]、湯明義等[8]、閻俊岳等[9]、丁非等[10]對東海海域SST變化分析結(jié)果相近,表明浙南沿岸SST的周期變化主要與厄爾尼諾、南方濤動(dòng)現(xiàn)象及太陽活動(dòng)等密切相關(guān)。

    大陳、南麂2站存在30 a的周期振蕩,與黃海仁等[11]提出的Nino區(qū)SST存在30 a以上時(shí)間尺度變化結(jié)論一致,分析原因可能與大陳、南麂2站均布設(shè)在海島,比坎門站更靠近外海,受大陸影響較小有關(guān)。

    湯明義等[8]分析結(jié)果表明,黑潮途徑海域SST存在10~11 a、6~7 a周期變化,南麂站的10 a周期變化及坎門、南麂2站的7.5 a周期變化與此接近,可能與黑潮流經(jīng)區(qū)域有關(guān)。

    3.4 累年月平均SST分布特征

    一年之中,浙南沿岸SST變化受到太陽輻射、臺灣暖流、浙江沿岸流及海陸性質(zhì)的影響,趨勢接近余弦曲線(圖3),呈一峰一谷的年分布特征;夏季月平均SST坎門站>南麂站>大陳站,估計(jì)與各站不同的下墊面環(huán)境有關(guān)。具體為2月為峰谷,此時(shí)浙南沿岸SST最低,隨著緯度的變化,總體表現(xiàn)為南高北低,其中南麂站SST為9.8 ℃,但由于站位布設(shè)受陸地氣候的影響及海水熱容量的不同,坎門站低于大陳站,為8.4 ℃, SST南北相差最大為1.4 ℃。從3月份開始,浙南沿岸SST進(jìn)入快速升溫期,到8月份SST達(dá)到最高,此時(shí)月平均SST坎門站最高,為28.8 ℃,大陳站最低,為26.7 ℃, SST南北相差最大為2.1 ℃。

    表 5 顯著周期統(tǒng)計(jì)Tab.5 Statistics of remarkable cycles

    通過對3個(gè)測站累年月平均SST進(jìn)行諧波分析,得到較為顯著的年循環(huán)及半年循環(huán),其振幅及初始位相見表6,其中,A1、φ1分別表示年循環(huán)的振幅和位相,A2、φ2分別表示半年循環(huán)的振幅和位相。諧波分析中,振幅決定大小,位相決定時(shí)間??梢钥闯?,年循環(huán)分量A1坎門站最高,半年循環(huán)分量A2南麂站最高;浙南沿岸SST變化主要由年循環(huán)分量決定。

    由于初始位相決定SST極值出現(xiàn)時(shí)間,浙南沿岸φ1變化范圍為233°~241°,根據(jù)各月年循環(huán)相角大小初步判斷,一年中浙南沿岸SST極值出現(xiàn)時(shí)間為8月中旬(最高)和 2月中旬(最低)左右,各站出現(xiàn)時(shí)間前后相差約7天,分析結(jié)果與SST實(shí)測數(shù)據(jù)吻合。

    表 6 SST主要諧波分析統(tǒng)計(jì)Tab.6 Statistics of SST major harmonic analysis

    圖 3 累年月平均SST變化曲線Fig.3 Curve of monthly mean SST

    圖 4 南麂站諧波擬合曲線Fig.4 Harmonic fitting curve of Nanji station

    以南麂觀測站為例,分別對年循環(huán)、年循環(huán)+半年循環(huán)的諧波分析結(jié)果與觀測數(shù)據(jù)進(jìn)行擬合(圖4),可以看出,擬合效果良好。 但是由于半年循環(huán)與年循環(huán)初始位相相差達(dá)111°,對年循環(huán)振幅有削弱的作用,使得年循環(huán)+半年循環(huán)擬合曲線極值減小。

    4 小 結(jié)

    (1)浙南沿岸SST總體受太陽輻射影響,呈現(xiàn)北低南高的地理分布特征。各站布設(shè)的地理環(huán)境、海水熱容量、不同流系等因素的影響,使得浙南沿岸SST分布各有特點(diǎn);

    (2)浙南沿岸SST總體呈現(xiàn)上升趨勢,上升速率為0.017℃/a;浙南沿岸北部測站上升速率大于南部;3個(gè)測站變化趨勢基本接近,20世紀(jì)60年代中后期之前,SST呈現(xiàn)下降趨勢;60年代中后期到70年代中期,SST呈現(xiàn)緩慢上升;70年代中期到80年代中期,SST開始又一輪下降過程;之后到2005年,SST一直呈現(xiàn)上升,且上升幅度明顯比前個(gè)上升過程大;2005年前后至2010年,呈現(xiàn)下降。

    (3)浙南沿岸SST前30 a基本處于偏冷期,由北往南各站偏冷年份比例分別為64%、60%和61%;2000-2010年間,大陳、坎門、南麂3站偏暖年數(shù)大幅上升,分別為7 a、8 a和6 a。

    (4)功率譜分析結(jié)果表明,浙南沿岸SST同時(shí)具有15 a 、6 a、5 a、4.3 a和3.8 a時(shí)間尺度的變化,主要與厄爾尼諾、南方濤動(dòng)現(xiàn)象及太陽活動(dòng)等密切相關(guān);大陳、南麂2個(gè)測站具有30 a時(shí)間尺度的變化,坎門、南麂2個(gè)測站具有7.5 a時(shí)間尺度的變化,南麂站具有10 a時(shí)間尺度的變化,可能與測站的選址、黑潮的流經(jīng)區(qū)域有關(guān)。

    (5)浙南沿岸SST年變化趨勢近似余弦曲線,其變化主要受年循環(huán)分量決定;年循環(huán)分量振幅坎門站最高,半年循環(huán)分量振幅南麂站最高;年循環(huán)、半年循環(huán)諧波擬合效果良好;南麂站半年循環(huán)分量對年循環(huán)分量有削弱的作用。

    致謝:感謝國家海洋局東海分局預(yù)報(bào)減災(zāi)處沙偉高級工程師對本文的大力指導(dǎo)!

    [1]余建民.慈溪市氣象站遷站前后氣象資料序列的均一性檢驗(yàn) [J].浙江氣象科技, 1999, 20(1): 54-55.

    [2]李慶祥, 劉小寧, 張洪政, 等.定點(diǎn)觀測氣候序列的均一性研究[J].氣象科技, 2003, 30(1): 3-22.

    [3]陳上及, 何維煥, 姚湜予, 等.中國近海海洋水文氣候季節(jié)的劃分 [J].海洋學(xué)報(bào), 1992, 14(6): 1-11.

    [4]鮑獻(xiàn)文, 萬修全, 高郭平, 等.渤海、黃海、東海AVHRR海表溫度場的季節(jié)變化特征 [J].海洋學(xué)報(bào), 2002, 24(5): 125-133.

    [5]湯明義, 劉宇中, 李洪海, 等.渤、黃海及東海北部平均表面水溫場季節(jié)變化特征及其成因的初步分析 [J].海洋學(xué)報(bào), 1989, 11(5): 544-553.

    [6]魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測技術(shù) [M].北京: 氣象出版社, 2007.

    [7]郭偉其, 沙偉, 沈紅梅, 等.東海沿岸海水表層溫度特征分析[J].海洋學(xué)報(bào), 2005, 27(5): 1-8.

    [8]湯明義, 尹杰, 曲如美, 等.東海黑潮途徑海域島嶼測站表面海溫年變化及年際變化特征 [J].海洋學(xué)報(bào), 1994,16(3): 12-20.

    [9]閻俊岳, 李江龍.東海及鄰近地區(qū)百年來的溫度變化 [J].海洋學(xué)報(bào), 1997,19(6): 121-126.

    [10]于非, 許一.東海海表面溫度長期變化趨勢研究 [J].海洋科學(xué)進(jìn)展, 2003, 21(4): 477-481.

    [11]黃海仁, 孫衛(wèi)國, 程炳巖.Nino海區(qū)SSTA的時(shí)頻結(jié)構(gòu)特征 [J].氣象科技, 2005, 33(3): 218-226.

    Sea surface temperature distribution characteristics along the southern coastal region of Zhejiang Province

    HAN Xiao-yan, PAN Xiao-dong, MA Lin-fang

    ( Wenzhou Oceanic Situation Monitoring Station, Wenzhou 325013, China)

    The SST temporal and spatial distribution in southern Zhejiang coast is presented as a result of statistical analysis based on 51 years observation data from Dachen, Kanmen and Nanji gauge stations.It shows that the overall SST along the coast is higher in the southern part and lower in the northern part, and the discrepancy between south and north parts becomes smallest in winter and largest in summer.Along the southern coast of Zhejiang, SST as a whole generally has a uptrend at rising rate of 0.017 ℃/a in last 51 years.In addition, there are oscillation periods of 3.5 a, 5-6 a and 11 a stay with SST changes.Mainly, the annual cyclical variation determines SST with a good match of harmonic fitting effect.

    the southern coastal region of Zhejiang; SST; homogeneous; power spectrum; harmonic

    P731.11

    A

    1001-6932(2011)06-0619-06

    2010-10-28;

    2011-06-08

    韓小燕(1977-),女,工程師,研究生在讀,主要從事海洋預(yù)報(bào)及觀測數(shù)據(jù)分析,電子郵箱:xiaoyong0728@163.com。

    猜你喜歡
    浙南時(shí)間尺度測站
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對稱性與守恒量
    GNSS鐘差估計(jì)中的兩種測站選取策略分析
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    浙南采風(fēng),余偶微恙,未克成行
    全球GPS測站垂向周年變化統(tǒng)計(jì)改正模型的建立
    測站分布對GPS解算ERP的影響分析
    城市勘測(2018年1期)2018-03-15 03:37:02
    大連市暴雨多時(shí)間尺度研究分析
    從血脈到肌理看盡浙南鳳凰部落
    浙南廊橋上的鳳凰部落
    欧美极品一区二区三区四区| 亚洲国产毛片av蜜桃av| 午夜免费观看性视频| 亚洲自偷自拍三级| 久久人人爽人人爽人人片va| 一边亲一边摸免费视频| 一二三四中文在线观看免费高清| 少妇 在线观看| 伊人久久国产一区二区| 久久久久国产精品人妻一区二区| 免费观看无遮挡的男女| 免费播放大片免费观看视频在线观看| 欧美精品人与动牲交sv欧美| 一个人看的www免费观看视频| av视频免费观看在线观看| 国产乱人偷精品视频| 人妻系列 视频| 亚洲四区av| 国产一区二区在线观看日韩| 欧美精品一区二区免费开放| 中文欧美无线码| 国产午夜精品久久久久久一区二区三区| 欧美三级亚洲精品| 日韩欧美一区视频在线观看 | 中文字幕av成人在线电影| 国产成人freesex在线| 看免费成人av毛片| 青春草视频在线免费观看| 51国产日韩欧美| 毛片女人毛片| 欧美xxxx黑人xx丫x性爽| 99久久精品国产国产毛片| 久久97久久精品| 亚洲一级一片aⅴ在线观看| 国产亚洲5aaaaa淫片| 国产高清不卡午夜福利| 亚洲精品成人av观看孕妇| 丝袜脚勾引网站| 99久久综合免费| 美女xxoo啪啪120秒动态图| 欧美精品国产亚洲| 亚洲欧美精品专区久久| av国产免费在线观看| 亚洲欧美精品专区久久| 国产亚洲5aaaaa淫片| 亚洲精品乱码久久久久久按摩| 国产精品一二三区在线看| 久久人人爽av亚洲精品天堂 | 亚洲av福利一区| 精品国产三级普通话版| 日韩三级伦理在线观看| 亚洲成人中文字幕在线播放| 欧美亚洲 丝袜 人妻 在线| 身体一侧抽搐| 男人舔奶头视频| 国产一区二区三区av在线| 十分钟在线观看高清视频www | 欧美成人午夜免费资源| 国产无遮挡羞羞视频在线观看| 秋霞在线观看毛片| 成人美女网站在线观看视频| 精品亚洲乱码少妇综合久久| 欧美 日韩 精品 国产| 99久久人妻综合| 尤物成人国产欧美一区二区三区| 国产精品蜜桃在线观看| 日韩电影二区| 蜜桃在线观看..| 尾随美女入室| 久久久精品免费免费高清| av线在线观看网站| 国产淫语在线视频| 国产在线一区二区三区精| 特大巨黑吊av在线直播| 亚洲av不卡在线观看| 麻豆乱淫一区二区| 久久97久久精品| 18禁在线播放成人免费| 黄色怎么调成土黄色| 99九九线精品视频在线观看视频| 夫妻午夜视频| 麻豆国产97在线/欧美| 亚洲国产精品专区欧美| 在线观看一区二区三区| 99热网站在线观看| 观看av在线不卡| 街头女战士在线观看网站| 欧美亚洲 丝袜 人妻 在线| 国产一区二区在线观看日韩| 亚洲美女视频黄频| 国产成人精品婷婷| 久久久久网色| 人人妻人人爽人人添夜夜欢视频 | 老熟女久久久| 少妇的逼水好多| 这个男人来自地球电影免费观看 | 18禁裸乳无遮挡免费网站照片| 男人爽女人下面视频在线观看| 久久精品夜色国产| h视频一区二区三区| 久久99热这里只频精品6学生| 国内揄拍国产精品人妻在线| 午夜日本视频在线| 亚洲国产精品专区欧美| 天堂8中文在线网| 嫩草影院新地址| 免费黄网站久久成人精品| 国产精品国产av在线观看| 亚洲经典国产精华液单| 一边亲一边摸免费视频| 国产美女午夜福利| 联通29元200g的流量卡| 黄片无遮挡物在线观看| 精品久久久久久久久av| 天天躁日日操中文字幕| 国产老妇伦熟女老妇高清| 美女视频免费永久观看网站| 赤兔流量卡办理| 高清午夜精品一区二区三区| 日韩av免费高清视频| 一区二区三区免费毛片| 在线亚洲精品国产二区图片欧美 | 精品久久久久久久久亚洲| 国产av国产精品国产| 亚洲欧美一区二区三区国产| 亚洲成人一二三区av| 国产成人精品福利久久| 久久99蜜桃精品久久| 91aial.com中文字幕在线观看| 91在线精品国自产拍蜜月| 少妇 在线观看| 中文在线观看免费www的网站| 在线看a的网站| 99久久精品热视频| 欧美zozozo另类| 伊人久久国产一区二区| 欧美最新免费一区二区三区| 哪个播放器可以免费观看大片| 激情 狠狠 欧美| 久久久色成人| 秋霞伦理黄片| 亚洲欧美精品专区久久| 国产精品一区二区在线不卡| 纯流量卡能插随身wifi吗| 日韩成人伦理影院| 国产真实伦视频高清在线观看| 91久久精品国产一区二区成人| 韩国高清视频一区二区三区| 国产毛片在线视频| 久久久久久久久大av| av在线播放精品| 国产精品国产三级国产专区5o| 欧美激情极品国产一区二区三区 | 18禁在线播放成人免费| 男人狂女人下面高潮的视频| 国产在线免费精品| 午夜激情久久久久久久| 午夜福利视频精品| 亚洲精品国产av蜜桃| 熟女电影av网| av国产免费在线观看| 麻豆国产97在线/欧美| 国产伦精品一区二区三区视频9| 日本av免费视频播放| 日本色播在线视频| 久久国产乱子免费精品| 国产精品无大码| 亚洲美女黄色视频免费看| 五月伊人婷婷丁香| 亚洲欧美成人综合另类久久久| 男女国产视频网站| 视频中文字幕在线观看| 免费久久久久久久精品成人欧美视频 | av免费在线看不卡| 午夜激情久久久久久久| 午夜免费鲁丝| 国产亚洲精品久久久com| 久热久热在线精品观看| 2021少妇久久久久久久久久久| 午夜福利网站1000一区二区三区| 两个人的视频大全免费| 亚洲第一区二区三区不卡| 亚洲色图av天堂| 亚洲国产高清在线一区二区三| 国产高潮美女av| 国产免费一级a男人的天堂| 人体艺术视频欧美日本| 成人黄色视频免费在线看| 97超视频在线观看视频| 国国产精品蜜臀av免费| 精品少妇久久久久久888优播| 亚洲综合精品二区| 国产久久久一区二区三区| 日韩一区二区视频免费看| 国产成人精品福利久久| 欧美精品一区二区大全| 日韩大片免费观看网站| 国产欧美日韩一区二区三区在线 | 能在线免费看毛片的网站| 欧美高清性xxxxhd video| 欧美一级a爱片免费观看看| 亚洲精品,欧美精品| 一级二级三级毛片免费看| 色婷婷久久久亚洲欧美| 人人妻人人看人人澡| 在线观看免费日韩欧美大片 | 午夜福利高清视频| 嫩草影院入口| 一级爰片在线观看| 高清黄色对白视频在线免费看 | 久久韩国三级中文字幕| 日日啪夜夜爽| 国产精品无大码| 天美传媒精品一区二区| 国产免费又黄又爽又色| 国产乱人视频| 成人漫画全彩无遮挡| 精品亚洲乱码少妇综合久久| av在线蜜桃| 最近中文字幕高清免费大全6| 亚洲av不卡在线观看| 五月开心婷婷网| 啦啦啦在线观看免费高清www| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 精品久久国产蜜桃| 国产精品久久久久久久久免| 久久精品国产亚洲av涩爱| 成人国产麻豆网| 日韩欧美精品免费久久| 麻豆成人av视频| 精华霜和精华液先用哪个| 综合色丁香网| 精品人妻熟女av久视频| 国产在线男女| 免费av不卡在线播放| 欧美性感艳星| 内地一区二区视频在线| 国产午夜精品久久久久久一区二区三区| 亚洲精品久久午夜乱码| 亚洲精品456在线播放app| 欧美精品一区二区大全| 黄色怎么调成土黄色| 亚洲丝袜综合中文字幕| 成人黄色视频免费在线看| 免费观看无遮挡的男女| 日韩欧美一区视频在线观看 | 大话2 男鬼变身卡| 免费观看无遮挡的男女| 免费少妇av软件| 哪个播放器可以免费观看大片| 大片免费播放器 马上看| 国产高清三级在线| 少妇人妻 视频| 亚洲精品国产av成人精品| 久久久久久久国产电影| 亚洲精品视频女| 久久影院123| 青春草国产在线视频| 精品熟女少妇av免费看| 中国国产av一级| av在线老鸭窝| 日本av免费视频播放| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 少妇 在线观看| 黑人高潮一二区| 一区二区三区精品91| 免费播放大片免费观看视频在线观看| 激情 狠狠 欧美| 日本欧美国产在线视频| 九九爱精品视频在线观看| 国产精品久久久久久久久免| 女人久久www免费人成看片| 精品熟女少妇av免费看| 国产精品国产三级专区第一集| 国产免费一级a男人的天堂| 99热网站在线观看| av在线播放精品| 在线观看免费高清a一片| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| 少妇的逼好多水| 极品少妇高潮喷水抽搐| 天天躁夜夜躁狠狠久久av| 精品久久国产蜜桃| 欧美97在线视频| 国产亚洲5aaaaa淫片| 欧美成人一区二区免费高清观看| 久久99热6这里只有精品| 十八禁网站网址无遮挡 | 草草在线视频免费看| 男人和女人高潮做爰伦理| 欧美变态另类bdsm刘玥| 日韩av免费高清视频| a级毛色黄片| 亚洲性久久影院| 97在线视频观看| 国产v大片淫在线免费观看| 久久国产精品大桥未久av | 91久久精品国产一区二区成人| 尤物成人国产欧美一区二区三区| 亚洲av国产av综合av卡| 精品一区二区三区视频在线| 岛国毛片在线播放| 日韩三级伦理在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产av新网站| 深夜a级毛片| xxx大片免费视频| 大片电影免费在线观看免费| 国产免费一区二区三区四区乱码| 香蕉精品网在线| 国产成人午夜福利电影在线观看| 少妇人妻一区二区三区视频| 久久久精品94久久精品| 国产片特级美女逼逼视频| av在线老鸭窝| 久久女婷五月综合色啪小说| 欧美精品一区二区免费开放| 国产免费视频播放在线视频| 欧美3d第一页| 18+在线观看网站| 亚洲av国产av综合av卡| 国产精品免费大片| 少妇裸体淫交视频免费看高清| 在线亚洲精品国产二区图片欧美 | 久久精品熟女亚洲av麻豆精品| 欧美丝袜亚洲另类| 三级国产精品片| 老熟女久久久| 80岁老熟妇乱子伦牲交| 日韩亚洲欧美综合| 超碰97精品在线观看| 久久久久久久国产电影| 全区人妻精品视频| 在线精品无人区一区二区三 | 亚洲精品成人av观看孕妇| a级毛色黄片| 久久这里有精品视频免费| av播播在线观看一区| 国产精品三级大全| 五月天丁香电影| 男人添女人高潮全过程视频| 中文欧美无线码| 久久影院123| 韩国av在线不卡| 黄色欧美视频在线观看| 寂寞人妻少妇视频99o| 我的女老师完整版在线观看| 高清不卡的av网站| 国产精品蜜桃在线观看| 女性生殖器流出的白浆| 黑人猛操日本美女一级片| 亚洲欧美中文字幕日韩二区| 国产午夜精品一二区理论片| 午夜福利在线在线| 国产精品99久久99久久久不卡 | 99九九线精品视频在线观看视频| 国产精品国产av在线观看| 视频区图区小说| 国产 精品1| 国产精品一区二区性色av| 免费观看性生交大片5| 精品亚洲成a人片在线观看 | 欧美丝袜亚洲另类| 成人亚洲精品一区在线观看 | 干丝袜人妻中文字幕| 波野结衣二区三区在线| 亚洲aⅴ乱码一区二区在线播放| 国产在线男女| 色婷婷av一区二区三区视频| 国产无遮挡羞羞视频在线观看| 国产91av在线免费观看| 秋霞伦理黄片| 天堂俺去俺来也www色官网| 一级爰片在线观看| 女性生殖器流出的白浆| 久久久午夜欧美精品| av免费在线看不卡| 我要看黄色一级片免费的| 新久久久久国产一级毛片| 91精品一卡2卡3卡4卡| 视频区图区小说| 亚洲高清免费不卡视频| 精品人妻熟女av久视频| 精品一品国产午夜福利视频| 亚洲一区二区三区欧美精品| 一区二区三区免费毛片| 国产日韩欧美亚洲二区| 最近中文字幕2019免费版| xxx大片免费视频| 嘟嘟电影网在线观看| 十八禁网站网址无遮挡 | 欧美丝袜亚洲另类| 国产熟女欧美一区二区| 视频中文字幕在线观看| 欧美激情极品国产一区二区三区 | 国产精品嫩草影院av在线观看| 欧美三级亚洲精品| 毛片一级片免费看久久久久| h日本视频在线播放| 久久影院123| 蜜桃久久精品国产亚洲av| 亚洲av国产av综合av卡| 国产 精品1| 亚洲精品日韩在线中文字幕| 国产亚洲精品久久久com| 成人毛片60女人毛片免费| videossex国产| 亚洲人成网站在线观看播放| 综合色丁香网| 免费看不卡的av| 欧美日韩视频高清一区二区三区二| 国产成人a区在线观看| 日本色播在线视频| 中文字幕亚洲精品专区| 国产精品成人在线| 国产欧美日韩一区二区三区在线 | 欧美日韩亚洲高清精品| 亚洲真实伦在线观看| 99久久精品国产国产毛片| 边亲边吃奶的免费视频| 亚洲精品国产av成人精品| 又黄又爽又刺激的免费视频.| 色吧在线观看| 老熟女久久久| 又粗又硬又长又爽又黄的视频| 亚洲图色成人| 在线观看av片永久免费下载| 婷婷色av中文字幕| 亚洲精品久久久久久婷婷小说| 久久99热6这里只有精品| 国产女主播在线喷水免费视频网站| 免费观看性生交大片5| 波野结衣二区三区在线| av卡一久久| 精品国产露脸久久av麻豆| 国产av国产精品国产| 人妻 亚洲 视频| 综合色丁香网| 熟女av电影| 久久青草综合色| 80岁老熟妇乱子伦牲交| 97超视频在线观看视频| 免费看不卡的av| 一级片'在线观看视频| 午夜免费观看性视频| 亚洲av福利一区| 蜜桃久久精品国产亚洲av| 亚洲一级一片aⅴ在线观看| 成年av动漫网址| 免费不卡的大黄色大毛片视频在线观看| 中文字幕免费在线视频6| 这个男人来自地球电影免费观看 | 免费观看无遮挡的男女| 国产伦在线观看视频一区| 成人国产麻豆网| 亚洲精品日韩在线中文字幕| 五月天丁香电影| 精品国产一区二区三区久久久樱花 | 乱码一卡2卡4卡精品| 美女高潮的动态| 亚洲最大成人中文| 七月丁香在线播放| 简卡轻食公司| 国产在线男女| 国产成人免费观看mmmm| 国产精品久久久久成人av| 国产精品国产av在线观看| av国产免费在线观看| 99久久综合免费| 成人高潮视频无遮挡免费网站| 国产成人免费观看mmmm| 免费在线观看成人毛片| 亚洲国产日韩一区二区| 我的女老师完整版在线观看| 欧美高清性xxxxhd video| 在线看a的网站| 成人漫画全彩无遮挡| 婷婷色av中文字幕| 国产亚洲5aaaaa淫片| 六月丁香七月| 国产精品国产av在线观看| 国内揄拍国产精品人妻在线| 国精品久久久久久国模美| 91狼人影院| av又黄又爽大尺度在线免费看| 欧美高清性xxxxhd video| 国产成人精品婷婷| 亚洲丝袜综合中文字幕| 精品久久久久久电影网| 蜜桃久久精品国产亚洲av| 亚洲精品乱码久久久久久按摩| 18禁动态无遮挡网站| 一区在线观看完整版| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 国产亚洲午夜精品一区二区久久| 成人高潮视频无遮挡免费网站| 欧美日韩视频高清一区二区三区二| 人妻 亚洲 视频| 免费观看无遮挡的男女| 色综合色国产| 亚洲精品日本国产第一区| 国产真实伦视频高清在线观看| 精品亚洲成国产av| 成人影院久久| 欧美日韩综合久久久久久| 最近2019中文字幕mv第一页| 亚洲欧洲国产日韩| 国产av一区二区精品久久 | 久久精品久久久久久久性| 久久久久久久久大av| 国产成人午夜福利电影在线观看| 美女内射精品一级片tv| 亚洲精品亚洲一区二区| 亚洲av在线观看美女高潮| 成人特级av手机在线观看| 九草在线视频观看| 丝袜脚勾引网站| 日韩伦理黄色片| 成年女人在线观看亚洲视频| 欧美zozozo另类| 成人毛片a级毛片在线播放| 女的被弄到高潮叫床怎么办| 高清av免费在线| 黄色配什么色好看| 简卡轻食公司| 少妇人妻 视频| 欧美变态另类bdsm刘玥| videossex国产| 国产精品99久久久久久久久| 亚洲av免费高清在线观看| 在线看a的网站| 晚上一个人看的免费电影| 中文字幕av成人在线电影| 91精品国产国语对白视频| 网址你懂的国产日韩在线| 成人影院久久| 乱码一卡2卡4卡精品| 亚洲精品亚洲一区二区| 亚洲高清免费不卡视频| 午夜免费男女啪啪视频观看| 亚洲伊人久久精品综合| 国产精品国产三级国产专区5o| 网址你懂的国产日韩在线| 97在线视频观看| 欧美高清成人免费视频www| 亚洲av.av天堂| 干丝袜人妻中文字幕| 天堂中文最新版在线下载| 亚洲四区av| 久久人妻熟女aⅴ| av女优亚洲男人天堂| 三级国产精品欧美在线观看| 亚洲国产最新在线播放| av天堂中文字幕网| 建设人人有责人人尽责人人享有的 | h日本视频在线播放| 黄色怎么调成土黄色| 国产女主播在线喷水免费视频网站| 乱系列少妇在线播放| 精品久久久久久电影网| 久久久a久久爽久久v久久| 久久久久久久精品精品| 亚洲欧美日韩东京热| 老司机影院成人| 国产永久视频网站| 久久精品国产亚洲av天美| 高清日韩中文字幕在线| 亚洲精品456在线播放app| 精品国产乱码久久久久久小说| 18禁动态无遮挡网站| 777米奇影视久久| 亚洲三级黄色毛片| 啦啦啦视频在线资源免费观看| 欧美激情极品国产一区二区三区 | kizo精华| 久久综合国产亚洲精品| 午夜福利在线观看免费完整高清在| 国产精品99久久99久久久不卡 | 国产成人精品一,二区| 97超视频在线观看视频| 大又大粗又爽又黄少妇毛片口| 久久久久久久精品精品| 亚洲最大成人中文| 最近2019中文字幕mv第一页| 久久99蜜桃精品久久| 美女中出高潮动态图| 伦理电影免费视频| 午夜福利在线观看免费完整高清在| 简卡轻食公司| 亚洲av.av天堂| 国产色婷婷99| 国产精品久久久久久av不卡| 老司机影院毛片| 国产精品久久久久久精品古装| 精华霜和精华液先用哪个| 一级爰片在线观看| 日本wwww免费看| 国模一区二区三区四区视频| 欧美xxxx黑人xx丫x性爽| 国产精品不卡视频一区二区| 男女国产视频网站| 九九爱精品视频在线观看| 国产精品无大码| 亚洲四区av| 波野结衣二区三区在线| 一二三四中文在线观看免费高清| 亚洲美女搞黄在线观看| 这个男人来自地球电影免费观看 | 卡戴珊不雅视频在线播放| 一级毛片黄色毛片免费观看视频| 国产成人aa在线观看|