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

    黃海海浪天氣時(shí)間尺度變化的數(shù)值模擬研究

    2017-02-09 09:19:38劉子洲陳國(guó)光
    海洋科學(xué) 2017年10期
    關(guān)鍵詞:波高時(shí)間尺度黃海

    劉子洲, 陳國(guó)光, 陸 雪

    ?

    黃海海浪天氣時(shí)間尺度變化的數(shù)值模擬研究

    劉子洲1, 陳國(guó)光2, 陸 雪3

    (1. 中國(guó)海洋大學(xué) 海洋與大氣學(xué)院, 山東 青島 266100; 2. 國(guó)家海洋局北海環(huán)境監(jiān)測(cè)中心, 山東 青島 266033; 3. 青島市氣象局, 山東 青島 266003)

    利用第三代海浪數(shù)值模式(SWAN)系統(tǒng)研究了黃海海浪有效波高的天氣時(shí)間尺度變化的時(shí)空分布特征和相關(guān)動(dòng)力學(xué)過程。結(jié)果表明黃海海浪有效波高的天氣變化強(qiáng)度(W)具有顯著的時(shí)間變化特征和空間分布特征。其多年平均值在黃海的中東部存在由南向北延伸的高值區(qū), 同時(shí)向兩側(cè)近海區(qū)域逐漸減小。W在冬季最大, 夏季最小。從11月到翌年5月,W月氣候態(tài)平均值的空間分布與其年平均值的空間分布類似; 從6月到同年10月,W的月氣候態(tài)平均值在黃海與東海的分界處存在較強(qiáng)的由黃海到東海的空間梯度。進(jìn)一步分析表明黃海海域的W以風(fēng)浪占主, 涌浪的貢獻(xiàn)遠(yuǎn)小于風(fēng)浪貢獻(xiàn)。數(shù)值實(shí)驗(yàn)表明, 黃海海浪有效波高的天氣時(shí)間尺度變化主要是由大于天氣變化周期的海面風(fēng)強(qiáng)迫通過四波非線性相互作用產(chǎn)生的。

    黃海; 海浪; 有效波高; 天氣時(shí)間尺度變化; 第三代海浪數(shù)值模式(SWAN)

    黃海是位于太平洋西部、我國(guó)大陸和朝鮮半島之間的一個(gè)陸架淺海, 南以長(zhǎng)江口-濟(jì)州島一線為界與東海相鄰, 西以渤海海峽與渤海相接。由于黃海海槽的存在, 黃海的地形大體分布特征是: 中部較深, 向兩邊逐漸變淺(圖1)。該海域主要受東亞季風(fēng)系統(tǒng)的控制, 夏季盛行偏南的夏季風(fēng), 冬季盛行偏北的冬季風(fēng)。同時(shí)由于黃海暖流等的存在, 黃海水文狀況還受到黑潮的顯著影響[1]。

    在東亞季風(fēng)和黑潮等大氣和海洋環(huán)流的影響下, 黃海水動(dòng)力環(huán)境場(chǎng)的時(shí)間變化和空間分布均比較復(fù)雜。基于觀測(cè)數(shù)據(jù)和數(shù)值模式模擬結(jié)果等資料, 前人對(duì)黃海環(huán)流和溫鹽場(chǎng)的季節(jié)變化[2]和年際及年代際變化[3-4]等作了大量研究。但是前人對(duì)黃海水動(dòng)力環(huán)境場(chǎng)的另一方面——海表面波浪作的研究卻相對(duì)較少。在國(guó)內(nèi)外的基礎(chǔ)和應(yīng)用研究中, 海表面波浪和波能流密度等的空間分布、時(shí)間變化及其動(dòng)力學(xué)過程一直是物理海洋學(xué)家關(guān)注的重點(diǎn)[5]。海浪不僅對(duì)海洋環(huán)流和海氣相互作用等具有非常重要的影響[6],而且是海上工程威脅、海洋災(zāi)害和新型能源的重要來源[7]。

    到目前為止, 前人對(duì)黃海海浪有效波高時(shí)空變化的研究主要集中在其季節(jié)變化、年際變化和長(zhǎng)期趨勢(shì)等低頻變化。對(duì)于季節(jié)變化, 陳紅霞等[8]根據(jù)衛(wèi)星觀測(cè)資料指出黃海海域海浪的有效波高于冬季最大, 夏季最小。該觀測(cè)結(jié)果與李訓(xùn)強(qiáng)等[9]的數(shù)值模擬結(jié)果并不一致, 但卻與陳國(guó)光等[10]的數(shù)值模擬結(jié)果一致。對(duì)于年際及更長(zhǎng)周期的變化, 鄭崇偉等[11]基于WAVEWATCH-III的模擬結(jié)果認(rèn)為黃海海浪有效波高的年際變化主要與厄爾尼諾-南方濤動(dòng)有關(guān)。Zheng等[12]指出黃海海浪在過去30年間存在顯著地長(zhǎng)期變化趨勢(shì)。除此之外, 也有部分工作研究了黃海海浪由臺(tái)風(fēng)過程引起的短期變化。Cheon等[13]利用海洋數(shù)值模式研究了布拉萬臺(tái)風(fēng)經(jīng)過黃海海域時(shí)產(chǎn)生的大浪的時(shí)空特征和動(dòng)力學(xué)過程。

    盡管前人的研究工作極大增進(jìn)了我們對(duì)黃海海浪場(chǎng)時(shí)空變化特征的認(rèn)識(shí), 但是仍存在諸多不足。近期通過對(duì)黃海海浪場(chǎng)衛(wèi)星觀測(cè)資料的分析, 作者發(fā)現(xiàn)黃海海浪有效波高存在顯著的天氣時(shí)間尺度的變化。因此本文擬利用第三代海浪數(shù)值模式系統(tǒng)地研究黃海海浪有效波高在天氣時(shí)間尺度上的變化特征和動(dòng)力學(xué)過程。

    圖1 模式計(jì)算區(qū)域水深(m)分布圖

    1 資料與方法

    1.1 觀測(cè)數(shù)據(jù)

    本文利用的觀測(cè)數(shù)據(jù)是由衛(wèi)星高度計(jì)觀測(cè)的有效波高數(shù)據(jù), 該觀測(cè)數(shù)據(jù)融合了ERS-1/2, TOPEX, Jason-1/2, GEOSAT和ENVISAT等多個(gè)衛(wèi)星的觀測(cè)結(jié)果并處理成格點(diǎn)化資料, 空間分辨率為1.0°×1.0°, 時(shí)間分辨率是1 d。衛(wèi)星觀測(cè)結(jié)果在外海海域與海表面浮標(biāo)觀測(cè)結(jié)果具有較好的一致性[14], 但是在淺水區(qū)誤差仍較大[15]。觀測(cè)數(shù)據(jù)的下載網(wǎng)址是http: // www.aviso.altimetry.fr/en/home.html。

    1.2 模式配置

    為了研究海浪有效波高天氣時(shí)間尺度變化的動(dòng)力學(xué)過程, 本文采用第三代海浪數(shù)值模式SWAN (Simulating WAves Nearshore)[16]版本41.01進(jìn)行了后報(bào)模擬和多項(xiàng)數(shù)值實(shí)驗(yàn)。SWAN模式采用波作用量(,)來描述隨機(jī)海浪場(chǎng), 其中波作用量滿足如下方程

    其中,為海浪的相對(duì)頻率,為其傳播方向;CC分別是在方向和方向上的傳播速度, 而則是在譜空間(,)上的傳播速度。右側(cè)代表能量源匯, 其表達(dá)式為

    其中右側(cè)各項(xiàng)依次代表海面風(fēng)輸入的能量、四波非線性相互作用傳輸?shù)哪芰?、白冠耗散的能量、三波非線性相互作用傳輸?shù)哪芰?、底摩擦耗散的能量、深度誘導(dǎo)破碎的能量。后3項(xiàng)主要是在淺水中起作用, 因此可稱為淺水過程項(xiàng)。

    本文中SWAN的模擬區(qū)域?yàn)?17°~132°E, 26°~ 41°N(圖1), 空間分辨率為0.2°×0.2°, 時(shí)間步長(zhǎng)為1 h, 結(jié)果每4 h輸出一次。模式的地形數(shù)據(jù)來自全球陸地海洋1′高程數(shù)據(jù)[17], 海面10 m風(fēng)場(chǎng)數(shù)據(jù)來自歐洲中期天氣預(yù)報(bào)中心(European Center for Medium-Range Weather Forecasts) Interim大氣再分析資料(ERA- Interim)[18]。模擬的時(shí)間范圍是2009年1月1日—2015年6月30日, 本文選取2009年10月以后的模擬結(jié)果進(jìn)行分析。

    2 結(jié)果

    2.1 譜分析

    圖2首先比較了衛(wèi)星觀測(cè)與模式后報(bào)模擬的黃海平均的海浪有效波高的天平均時(shí)間序列。從圖中可以看出, 兩者吻合很好, 其同期相關(guān)系數(shù)為0.74 (高于95%的置信檢驗(yàn))。除了呈現(xiàn)出顯著的季節(jié)變化外[10],黃海海浪的有效波高還呈現(xiàn)出強(qiáng)烈的高頻變化。之后圖中分別給出了兩時(shí)間序列的功率譜密度。功率譜顯示衛(wèi)星觀測(cè)和模式后報(bào)模擬的有效波高均呈現(xiàn)出周期短于12 d的天氣時(shí)間尺度變化。為了查看該高頻變化的強(qiáng)度, 我們對(duì)有效波高的時(shí)間序列進(jìn)行了高通濾波, 提取出周期小于12 d的信號(hào), 記為有效波高天氣異常。圖2最后分別給出了衛(wèi)星觀測(cè)和后報(bào)模擬的有效波高天氣異常的標(biāo)準(zhǔn)差。從圖中可以看出, 兩者均從南往北逐漸減小。但同時(shí), 兩者之間也有較大差別。一方面, 衛(wèi)星觀測(cè)的結(jié)果比模式模擬的結(jié)果偏小, 這可能是因?yàn)槟J侥M所使用的海表面10 m風(fēng)速比實(shí)際風(fēng)速偏大[19]; 另一方面, 模式后報(bào)模擬結(jié)果顯示有效波高天氣異常的標(biāo)準(zhǔn)差從黃海中東部向沿岸區(qū)域逐漸減小, 而該現(xiàn)象在衛(wèi)星觀測(cè)結(jié)果中并不明顯, 這應(yīng)該與衛(wèi)星觀測(cè)結(jié)果在淺海區(qū)域誤差較大有關(guān)[15]。

    2.2 季節(jié)變化特征

    黃海海浪有效波高在天氣時(shí)間尺度上的變化有沒有季節(jié)差異呢?為此我們定義了有效波高在天氣時(shí)間尺度上的變化強(qiáng)度W(), 其定義如下: 對(duì)于某一有效波高時(shí)間序列, 首先進(jìn)行截?cái)嘀芷跒?2 d的高通濾波以獲得天氣時(shí)間尺度上的異常值, 然后以某一時(shí)刻為中心取前后共31 d的異常值時(shí)間序列計(jì)算其標(biāo)準(zhǔn)差作為有效波高在該時(shí)刻的天氣時(shí)間尺度變化強(qiáng)度。因此W()在統(tǒng)計(jì)意義上描述了有效波高在天氣時(shí)間尺度上的變化強(qiáng)弱。圖3首先分別給出了整個(gè)黃海、北黃海(37°N以北)和南黃海(37°N以南)分別平均的W時(shí)間序列。總體而言, 南黃海的W始終大于北黃海, 這與圖2最后描述的多年平均的W由南向北逐漸減小的空間分布特征一致。北黃海和南黃海分別平均的W與整個(gè)黃海平均的W之間的相關(guān)系數(shù)分別為0.92和0.99, 而它們的同期相關(guān)系數(shù)為0.87, 說明整個(gè)黃海海域的W呈現(xiàn)基本一致的時(shí)間變化。從圖中可以看出W的低頻變化主要以季節(jié)變化為主, 基本上在冬季最大, 而于夏季最小。同時(shí)我們也注意到這三個(gè)區(qū)域分別平均的W均在2011年8月、2012年9月和2014年8月等非冬季月份出現(xiàn)較大值, 從而表明W還具有其他時(shí)間尺度的變化。之后圖3進(jìn)一步給出了整個(gè)研究時(shí)間段內(nèi)W的標(biāo)準(zhǔn)差, 其空間分布與其平均值(圖2)類似, 均在黃海中東部最大, 并向沿岸區(qū)域逐漸減小。

    圖2 衛(wèi)星觀測(cè)和模式后報(bào)模擬得到的有效波高天氣異常

    隨后本文計(jì)算了W的月氣候態(tài)平均值, 結(jié)果如圖4所示, 其中粉紅色點(diǎn)和灰色點(diǎn)分別指示W(wǎng)的最大值與最小值發(fā)生的月份。在黃海內(nèi),W的最大值在大部分海域出現(xiàn)在12月, 而在少數(shù)沿岸海域出現(xiàn)在11月;W的最小值在大部分海域發(fā)生在6月, 而在江蘇—山東沿岸的部分海域則出現(xiàn)在7月。因?yàn)槿鄙儆^測(cè)資料的支持,W極大值發(fā)生的月份在部分近海海域與深海并不一致的現(xiàn)象還需要在將來進(jìn)一步的深入研究。W月氣候態(tài)平均值的空間分布也呈現(xiàn)出顯著地季節(jié)變化。從11月到翌年5月,W的月氣候態(tài)平均值均在黃海的中東部存在自南向北延伸的高值區(qū), 與其年平均值的空間分布類似(圖2)。從6月到同年的10月,W的月氣候態(tài)平均在黃海與東海的分界處存在較強(qiáng)的由黃海指向東海的空間梯度, 其等值線大致與長(zhǎng)江口-濟(jì)州島連線平行。在5月和6月, 南黃海的中東部區(qū)域和北部區(qū)域分別存在一個(gè)局地極大值區(qū)。

    圖3 黃海海浪SW的時(shí)間變化及其標(biāo)準(zhǔn)差的空間分布

    按照定義, 實(shí)際海洋中的波浪可以分為風(fēng)浪和涌浪兩部分, 那么風(fēng)浪和涌浪有效波高在天氣時(shí)間尺度上又有怎樣的變化呢?為此本文分別計(jì)算了風(fēng)浪和涌浪的W, 結(jié)果如圖5所示??傮w而言, 混合浪有效波高的天氣變化主要是由風(fēng)浪占主, 這與其季節(jié)變化的特征類似[10]。涌浪在南黃海的貢獻(xiàn)要大于在北黃海的貢獻(xiàn), 這應(yīng)該是由于南黃海受到來自東海及太平洋海域傳來的較強(qiáng)涌浪的影響。另外, 涌浪在夏天的部分月份也存在極大值, 如前面提到的2011年8月、2012年9月和2014年8月, 這可能是與臺(tái)風(fēng)浪有關(guān)[13]。

    2.3 動(dòng)力學(xué)過程探討

    海表面波浪都是由海面風(fēng)強(qiáng)迫產(chǎn)生的, 因此有必要探討W與海表面風(fēng)場(chǎng)的天氣時(shí)間尺度變化強(qiáng)度之間的關(guān)系。為此本文也按照類似的方法利用ERA-Interim風(fēng)場(chǎng)資料計(jì)算了海面10 m風(fēng)速的天氣時(shí)間尺度變化強(qiáng)度, 記為W。經(jīng)過比較,W的季節(jié)變化和空間分布與W的季節(jié)變化和空間分布存在很大差別, 如圖6所示(其中粉紅色點(diǎn)和灰色點(diǎn)分別指示W(wǎng)的最大值與最小值發(fā)生的月份), 因此海面風(fēng)的天氣時(shí)間尺度變化可能不是導(dǎo)致海浪有效波高天氣時(shí)間尺度變化的直接原因。

    圖4 黃海海浪SW的逐月氣候態(tài)平均

    為了探討黃海海浪有效波高天氣時(shí)間尺度變化的動(dòng)力學(xué)過程, 本文設(shè)計(jì)了6個(gè)數(shù)值實(shí)驗(yàn), 每個(gè)實(shí)驗(yàn)的模式配置如表1所示。圖7給出了利用不同數(shù)值實(shí)驗(yàn)?zāi)M結(jié)果計(jì)算得到的W時(shí)間序列。為了便于比較, 圖7中實(shí)驗(yàn)4的結(jié)果為真實(shí)結(jié)果的1/3。

    實(shí)驗(yàn)1和實(shí)驗(yàn)2比較了天氣時(shí)間尺度風(fēng)場(chǎng)(周期小于12 d)和低頻風(fēng)場(chǎng)(周期大于12 d)對(duì)W的影響。如圖7中灰色線和紅色線所示, 不管是在南黃海還是北黃海實(shí)驗(yàn)1的結(jié)果均遠(yuǎn)小于后報(bào)模擬結(jié)果, 而實(shí)驗(yàn)2的結(jié)果則與后報(bào)模擬結(jié)果非常接近。兩者得到的W與后報(bào)模擬得到的W之間的相關(guān)系數(shù)分別為0.5和0.8左右(表1), 這說明黃海海浪有效波高在天氣時(shí)間尺度上的變化主要是由海面風(fēng)場(chǎng)中較長(zhǎng)周期(大于天氣變化周期)的信號(hào)引起的, 而風(fēng)場(chǎng)中天氣時(shí)間尺度的信號(hào)則起次要作用。

    實(shí)驗(yàn)3—實(shí)驗(yàn)6進(jìn)一步探討了公式(2)中的不同物理過程對(duì)W的影響。實(shí)驗(yàn)3得到的W與后報(bào)模擬得到的W不僅大小一致, 而且同期相關(guān)系數(shù)非常高(≈0.98)。因此, 淺水過程項(xiàng)對(duì)南黃?;蛘弑秉S海區(qū)域平均的W均影響較小, 這是因?yàn)闇\水過程項(xiàng)主要在近海起作用[10], 而南黃?;蛘弑秉S海區(qū)域平均的W主要體現(xiàn)深水區(qū)域的有效波高天氣時(shí)間尺度變化強(qiáng)度特征。如果同時(shí)忽略四波非線性相互作用和白冠耗散等兩個(gè)深水過程, 那么實(shí)驗(yàn)4得到的W不僅在數(shù)值上較大于后報(bào)模擬結(jié)果(前者約為后者的3倍), 而且兩者的相關(guān)性也較低(≈0.1~0.4)。實(shí)驗(yàn)5的結(jié)果不僅在數(shù)值上遠(yuǎn)小于后報(bào)模擬結(jié)果, 而且與后者的同期相關(guān)性也非常低(≈0.2~0.3), 這說明四波非線性相互作用對(duì)有效波高的W非常重要。實(shí)驗(yàn)6得到的結(jié)果總體上比較接近后報(bào)模擬結(jié)果, 而且兩者之間的同期相關(guān)系數(shù)也相對(duì)較高(≈0.87), 這表明白冠耗散過程起次要作用。

    圖5 風(fēng)浪和涌浪的SW時(shí)間序列

    3 結(jié)論

    衛(wèi)星觀測(cè)表明黃海海浪的有效波高呈現(xiàn)出顯著地天氣時(shí)間尺度的變化(周期小于12 d), 為此本文利用第三代海浪數(shù)值模式SWAN系統(tǒng)研究了該變化的時(shí)空分布特征和相關(guān)動(dòng)力學(xué)過程, 具體結(jié)論如下:

    1) 黃海海浪有效波高的天氣變化強(qiáng)度W具有顯著的時(shí)間變化特征和空間分布特征。其多年平均值在黃海的中東部存在由南向北延伸的高值區(qū), 同時(shí)向兩側(cè)近海區(qū)域逐漸減小。W的低頻變化以季節(jié)變化為主, 基本是冬天最大, 夏天最小。W的空間分布也具有明顯的季節(jié)差異。從11月到翌年5月,W月氣候態(tài)平均值的空間分布與其年平均值的空間分布類似; 從6月到同年的10月,W的月氣候態(tài)平均值在黃海與東海的分界處存在較強(qiáng)的空間梯度:W由東海向黃海迅速減小, 其等值線大致與長(zhǎng)江口-濟(jì)州島連線平行。在5月和6月, 南黃海的中東部區(qū)域和北部區(qū)域分別存在一個(gè)極大值區(qū)。W以風(fēng)浪占主, 涌浪的貢獻(xiàn)遠(yuǎn)小于風(fēng)浪貢獻(xiàn)。

    2) 數(shù)值實(shí)驗(yàn)表明黃海海浪有效波高的天氣時(shí)間尺度變化并不是由海面風(fēng)場(chǎng)的天氣時(shí)間尺度變化導(dǎo)致的, 而是由較大變化周期的風(fēng)場(chǎng)主導(dǎo)的。進(jìn)一步的數(shù)值實(shí)驗(yàn)表明在影響深水區(qū)域波浪能量源匯的6個(gè)物理過程中, 四波非線性相互作用具有主導(dǎo)作用, 而白冠耗散起次要作用。因此, 黃海海浪有效波高的天氣時(shí)間尺度變化主要是由大于天氣變化周期的海面風(fēng)強(qiáng)迫通過四波非線性相互作用產(chǎn)生的。

    圖6 黃海海面10 m風(fēng)速WW的逐月氣候態(tài)平均

    表1 后報(bào)模擬和數(shù)值實(shí)驗(yàn)設(shè)置及結(jié)果

    圖7 后報(bào)模擬及數(shù)值實(shí)驗(yàn)得到的區(qū)域平均的有效波高的SW的時(shí)間序列

    [1] 臧家業(yè), 湯毓詳, 鄒娥梅, 等. 黃海環(huán)流的分析[J]. 科學(xué)通報(bào), 2001, 46(增): 7-15.Zang Jiaye, Tang Yuxiang, Zou Emei, et alAnalysis of Yellow Sea circulation[J]. Chinese Science Bulletin, 2001, 46(S): 7-15.

    [2] Pang Chongguang, Bai Xuezhi, BaiHong, et alSeasonal evolution of circulation and thermal structure in the Yellow Sea[J]. Chinese Journal of Oceanology and Limnology, 2005, 23(3): 269-274.

    [3] Han Guoqi, Huang Weigen. Pacific Decadal Oscillation and sea level variability in the Bohai, Yellow, and East China Seas[J]. Journal of Physical Oceanography, 2005, 38(12): 2772-2783.

    [4] Zhang Liping, Wu Lixin, Lin Xiaopei, et al. Modes and mechanisms of sea surface temperature low-frequency variations over the coastal China seas[J]. Journal of Geophysical Research-Oceans, 115: C08031.

    [5] 管長(zhǎng)龍. 我國(guó)海浪理論及預(yù)報(bào)研究的回顧與展望[J]. 青島海洋大學(xué)學(xué)報(bào), 2000, 30(4): 549-556. Guan Changlong. A review of history and prospect for study of sea wave theory and its forecast in China[J]. Journal of Ocean University of Qingdao, 2000, 30(4): 549-556.

    [6] Qiao Fangli, Ma Jian, Xia Changshui, et al. Influence of the surface wave-induced mixing and tidal mixing on the vertical temperature structure of the Yellow and East China Seas in summer[J]. Progress in Natural Science, 2006, 16(7): 739-746.

    [7] 許富祥, 余宙文. 中國(guó)近海及其鄰近海域?yàn)?zāi)害性海浪監(jiān)測(cè)和預(yù)報(bào)[J]. 海洋預(yù)報(bào), 1998, 15(3): 63-68. Xu Fuxiang, Yu Zhouwen. Monitoring and prediction of disastrous sea waves in the coastal waters of China and its adjacent areas[J]. Marine Forecasts, 1998, 15(3): 63-68.

    [8] 陳紅霞, 華鋒, 袁業(yè)立. 中國(guó)近海及鄰近海域海浪的季節(jié)特征及其時(shí)間變化[J]. 海洋科學(xué)進(jìn)展, 2006, 24(4): 407-415. Chen Hongxia, Hua Feng, Yuan Yeli. Seasonal characteristics and temporal variations of ocean wave in the Chinese offshore waters and adjacent sea areas[J]. Advances in Marine Science, 2006, 24(4): 407-415.

    [9] 李訓(xùn)強(qiáng), 鄭崇偉, 蘇勤, 等. 1988-2009年中國(guó)海波候、風(fēng)候統(tǒng)計(jì)分析[J]. 中國(guó)海洋大學(xué)學(xué)報(bào), 2012, 42(S1): 1-9. Li Xunqiang, Zheng Chongwei, Su Qin, et al. Wave climate and wind climate analysis in the China Sea from 1988 to 2009[J]. Periodical of Ocean University of China, 2012, 42: 1-9.

    [10] 陳國(guó)光, 翟方國(guó), 李培良, 等. 黃海海浪季節(jié)變化的數(shù)值模擬研究[J]. 海洋科學(xué), 2016, 40(11): 155-168. Chen Guoguang, Zhai Fangguo, Li Peiliang, et al. Numerical study of wave height seasonality in the Yellow Sea[J]. Marine Sciences, 2016, 40(11): 155-168

    [11] 鄭崇偉, 劉鐵軍, 錢粵海. 中國(guó)海海表風(fēng)場(chǎng)、海浪場(chǎng)與El Ni?o的相關(guān)性分析[J]. 云南大學(xué)學(xué)報(bào), 2014, 36(2): 214-223. Zheng Chongwei, Liu Tiejun, Qian Yuehai. The relationship between sea surface wind field, wave field and El Nino in the China sea[J]. Journal of Yunnan University, 2014, 36(2): 214-223.

    [12] Zheng Chongwei, Li Chongyin. Variation of the wave energy and significant wave height in the China Sea and adjacent waters[J]. Renewable and Sustainable Energy Reviews, 2015, 43: 381-387.

    [13] Cheon J K, Mu J W, Yong C J, et al. Simulation of the extreme waves generated by typhoon Bolaven (1215) in the East China Sea and Yellow Sea[J]. Acta Oceanol Sin, 2015, 34(12): 19-28.

    [14] Hwang P A, Bratos S M, Teague W J, et al. Wind and waves in the Yellow and East China Seas: A comparison of spaceborne altimeter measurements and model results[J]. Journal of Oceanography, 1999, 55(2): 307-325.

    [15] 楊樂. 衛(wèi)星雷達(dá)高度計(jì)在中國(guó)近海及高海況下遙感反演算法研究[D]. 南京: 南京理工大學(xué), 2009. Yang Le. Study on satellite radar altimeter retrieval algorithms over coastal seas and under high sea state events[D]. Nanjing: Nanjing University of Science and Technology, 2009.

    [16] Booij N, Ris R, Holthuijsen L. A third-generation wave model for coastal regions. Part I, Model description and validation[J]. Journal of Geophysical Research, 1999, 104(4): 7649-7666.

    [17] Smith W, Sandwell D. Global sea floor topography from satellite altimetry and ship depth soundings[J]. Science, 1997, 277: 1956-1962.

    [18] Dee D, Uppala S, Simmons A, et al. The ERA-Interim reanalysis: Configuration and performance of the data assimilation system[J]. Quarterly Journal of the Royal Meteorological Society, 2011, 137(656): 553-597.

    [19] Song Lina, Liu Zhiliang, Wang Fan. Comparison of wind data from ERA-Interim and buoys in the Yellow and East China Seas[J]. Chinese Journal of Oceanology and Limnology, 2015, 33(1): 282-288.

    (本文編輯: 劉珊珊)

    Numerical study on the variability of wave height on the weather time scale in the Yellow Sea

    LIU Zi-zhou1, CHEN Guo-guang2, LU Xue3

    (1. College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China; 2. North China Sea Environment Monitoring Center, State Oceanic Administration, Qingdao266033, China; 3. Qingdao Meteorological Administration, Qingdao 266003, China)

    Satellite observations indicate that the significant wave height (SWH) in the Yellow Sea (YS) exhibits significant variations in periods shorter than 12 days. The current study intends to investigate in detail the characteristics of the temporal and spatial variations in the weather with respect to the SWH in the YS as well as the associated dynamic processes based on Simulating WAves Nearshore (SWAN), a third-generation wave model. The results indicate that the strength of the weather variations in the SWH (denoted asW) shows strong characteristics in its temporal variations and spatial distributions. There is a high-value tongue ofWextending from the south to north in the eastern–central part of the YS. Meanwhile,Wdecreases shoreward. It is largest in winter and smallest in summer. From November to May in the following year, the monthly climatology ofWshows a spatial distribution similar to that of its annual mean. From June to October, however, there is a strong gradient in the monthly climatology ofWperpendicular to the border between the YS and the East China Sea. Further analysis indicates that theWin the YS is dominated by windsea, with secondary contributions from swell. Numerical experiments prove that the weather variations in the SWH are caused by wind forcing with longer periods via four-wave nonlinear interactions.

    the Yellow Sea; ocean wave; significant wave height; weather variation; Simulating WAves Nearshore (SWAN)

    Nov. 17, 2016

    P731

    A

    1000-3096(2017)10-0077-09

    10.11759/hykx20161117001

    2016-11-17;

    2016-12-20

    青島市博士后資助項(xiàng)目(86160504003)

    [Qingdao Postdoctoral Funding Project , No.86160504003]

    劉子洲(1987-), 男, 山東淄博人, 碩士, 實(shí)驗(yàn)師, 主要研究方向?yàn)榻>?xì)化預(yù)報(bào)及海洋調(diào)查, E-mail: lzz2013@ouc.edu.cn; 陳國(guó)光, 通信作者, 學(xué)士, E-mail: 13808960499@163.com

    猜你喜歡
    波高時(shí)間尺度黃海
    基于FHDI-GNWM 數(shù)據(jù)的全球超越概率波高宏觀分布特征分析
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對(duì)稱性與守恒量
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    非平整港池的多向不規(guī)則波試驗(yàn)研究
    東方濕地 黃海明珠
    黃海簡(jiǎn)介
    三角恒等變換去哪兒了
    大連市暴雨多時(shí)間尺度研究分析
    国产精品人妻久久久影院| 亚洲国产欧美人成| 一区二区三区高清视频在线| 麻豆一二三区av精品| 亚洲国产精品国产精品| 国产淫片久久久久久久久| 在线观看美女被高潮喷水网站| 国产一区二区激情短视频| 少妇的逼好多水| 欧美日韩精品成人综合77777| 久久精品人妻少妇| 高清午夜精品一区二区三区 | 日本成人三级电影网站| 免费黄网站久久成人精品| 成人亚洲精品av一区二区| 九九爱精品视频在线观看| 性欧美人与动物交配| 精品久久久久久久末码| 成人特级av手机在线观看| 国产精品国产三级国产av玫瑰| 精华霜和精华液先用哪个| 欧美一区二区精品小视频在线| 亚洲中文字幕一区二区三区有码在线看| 欧美区成人在线视频| 成年版毛片免费区| 亚洲欧美成人精品一区二区| 日韩大尺度精品在线看网址| 国产综合懂色| 成人亚洲精品av一区二区| 亚洲欧美中文字幕日韩二区| 一级毛片我不卡| 国产精品久久久久久av不卡| 国产黄片视频在线免费观看| 免费观看人在逋| 国产一区二区亚洲精品在线观看| 日韩欧美三级三区| 波多野结衣高清作品| 国产精品一及| 国产精品免费一区二区三区在线| 国产精品一及| 国产精品久久久久久精品电影| 搡女人真爽免费视频火全软件| 国内揄拍国产精品人妻在线| 美女大奶头视频| 日产精品乱码卡一卡2卡三| 天堂av国产一区二区熟女人妻| 一卡2卡三卡四卡精品乱码亚洲| 亚洲不卡免费看| 尤物成人国产欧美一区二区三区| 成熟少妇高潮喷水视频| 嘟嘟电影网在线观看| 免费看日本二区| 中文字幕制服av| 精品久久久久久成人av| 老司机福利观看| 人妻系列 视频| 久久午夜亚洲精品久久| 国模一区二区三区四区视频| 欧美zozozo另类| 久久午夜亚洲精品久久| 亚洲四区av| 午夜激情福利司机影院| 在线免费观看不下载黄p国产| 男人和女人高潮做爰伦理| av视频在线观看入口| 免费大片18禁| 色综合站精品国产| 狂野欧美激情性xxxx在线观看| 亚洲精品粉嫩美女一区| 天堂√8在线中文| av免费在线看不卡| 亚洲欧美成人综合另类久久久 | 亚洲人成网站在线播放欧美日韩| 最近手机中文字幕大全| 在线观看美女被高潮喷水网站| 久久久久久大精品| 久久久久久国产a免费观看| 九色成人免费人妻av| 日韩三级伦理在线观看| 亚洲一级一片aⅴ在线观看| 国产 一区精品| 3wmmmm亚洲av在线观看| 深夜a级毛片| 91av网一区二区| 日本五十路高清| av天堂在线播放| 国产极品天堂在线| 国产爱豆传媒在线观看| 91午夜精品亚洲一区二区三区| 91午夜精品亚洲一区二区三区| 日韩av不卡免费在线播放| 不卡一级毛片| 国产午夜福利久久久久久| 1000部很黄的大片| 欧美区成人在线视频| 身体一侧抽搐| 国产熟女欧美一区二区| 你懂的网址亚洲精品在线观看 | 精品人妻一区二区三区麻豆| 久久99热6这里只有精品| 一级黄片播放器| 国产麻豆成人av免费视频| 一本一本综合久久| 亚洲国产精品成人久久小说 | 91精品一卡2卡3卡4卡| 久久亚洲国产成人精品v| 久久久欧美国产精品| 欧美变态另类bdsm刘玥| 国产探花在线观看一区二区| 日韩欧美国产在线观看| 国产中年淑女户外野战色| 男人和女人高潮做爰伦理| 成人无遮挡网站| 日韩高清综合在线| 亚洲av熟女| 国模一区二区三区四区视频| 在线免费十八禁| 秋霞在线观看毛片| 性欧美人与动物交配| 亚洲国产精品久久男人天堂| 级片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 丝袜喷水一区| 联通29元200g的流量卡| 国产成人a区在线观看| 亚洲成av人片在线播放无| 亚洲人成网站在线观看播放| 午夜久久久久精精品| 久久精品国产亚洲av天美| 亚洲人成网站在线播| 国产色爽女视频免费观看| 亚洲天堂国产精品一区在线| 亚洲久久久久久中文字幕| 嘟嘟电影网在线观看| 老司机福利观看| 麻豆乱淫一区二区| 久久久国产成人精品二区| 成熟少妇高潮喷水视频| 中文字幕人妻熟人妻熟丝袜美| 中文字幕制服av| 国产色爽女视频免费观看| 啦啦啦观看免费观看视频高清| 日日摸夜夜添夜夜添av毛片| 国产免费男女视频| 小蜜桃在线观看免费完整版高清| 18+在线观看网站| 国产精品国产高清国产av| 伦理电影大哥的女人| 国语自产精品视频在线第100页| 搞女人的毛片| 夫妻性生交免费视频一级片| 白带黄色成豆腐渣| 波野结衣二区三区在线| 两个人的视频大全免费| 99国产极品粉嫩在线观看| 哪个播放器可以免费观看大片| 国产精品久久久久久精品电影| 狂野欧美激情性xxxx在线观看| 两个人视频免费观看高清| 亚洲国产精品成人综合色| 97热精品久久久久久| 午夜激情欧美在线| 老司机福利观看| 99久久人妻综合| 久久久久久久亚洲中文字幕| 久久久久免费精品人妻一区二区| 天天一区二区日本电影三级| 看免费成人av毛片| 人体艺术视频欧美日本| 亚洲国产高清在线一区二区三| 亚洲精品亚洲一区二区| 精品不卡国产一区二区三区| 99久久九九国产精品国产免费| 精品一区二区三区人妻视频| 久久热精品热| 午夜福利在线在线| 久久人人精品亚洲av| 黑人高潮一二区| 噜噜噜噜噜久久久久久91| 少妇的逼好多水| 亚洲七黄色美女视频| 免费搜索国产男女视频| 国产精品久久久久久精品电影小说 | 国产成人91sexporn| 亚州av有码| 乱人视频在线观看| 精品免费久久久久久久清纯| 日本一二三区视频观看| a级毛色黄片| 高清毛片免费看| 免费一级毛片在线播放高清视频| 天堂av国产一区二区熟女人妻| 欧美人与善性xxx| 日本成人三级电影网站| 男女啪啪激烈高潮av片| 日韩中字成人| 久久精品影院6| 少妇裸体淫交视频免费看高清| 一本久久精品| 精品一区二区三区视频在线| 99热这里只有是精品在线观看| 色尼玛亚洲综合影院| 在线播放无遮挡| 熟女电影av网| 91aial.com中文字幕在线观看| 99久久精品热视频| 听说在线观看完整版免费高清| 性色avwww在线观看| 久久精品国产亚洲网站| 日韩欧美国产在线观看| 亚洲精华国产精华液的使用体验 | 激情 狠狠 欧美| 91久久精品国产一区二区成人| 国产美女午夜福利| 干丝袜人妻中文字幕| 欧美丝袜亚洲另类| 婷婷色av中文字幕| 成人综合一区亚洲| 国产精品一区二区三区四区免费观看| 一进一出抽搐gif免费好疼| 一区福利在线观看| 欧美性猛交╳xxx乱大交人| 亚洲最大成人手机在线| 亚洲欧美清纯卡通| 欧美3d第一页| 成人综合一区亚洲| 99久久九九国产精品国产免费| 99国产极品粉嫩在线观看| 日日撸夜夜添| 亚洲国产欧洲综合997久久,| 日日干狠狠操夜夜爽| 可以在线观看的亚洲视频| 在线播放国产精品三级| 天堂影院成人在线观看| 久久精品国产亚洲av天美| 国产高清三级在线| 亚洲图色成人| 欧美最黄视频在线播放免费| av天堂中文字幕网| 日韩欧美精品v在线| 久久久久久国产a免费观看| 97热精品久久久久久| 嘟嘟电影网在线观看| 久久精品国产亚洲av涩爱 | 在线观看午夜福利视频| 成年免费大片在线观看| 观看免费一级毛片| 网址你懂的国产日韩在线| 只有这里有精品99| 波多野结衣高清作品| 色播亚洲综合网| 成年免费大片在线观看| 美女被艹到高潮喷水动态| 国产老妇女一区| 亚洲精品乱码久久久久久按摩| 99热只有精品国产| 久久韩国三级中文字幕| 麻豆av噜噜一区二区三区| 啦啦啦韩国在线观看视频| 精品不卡国产一区二区三区| av免费观看日本| 狂野欧美激情性xxxx在线观看| 在线观看午夜福利视频| 如何舔出高潮| 久久久欧美国产精品| 亚洲中文字幕日韩| 少妇被粗大猛烈的视频| 国产视频首页在线观看| 亚洲欧美日韩高清专用| 午夜久久久久精精品| 亚洲色图av天堂| 亚洲成人中文字幕在线播放| 淫秽高清视频在线观看| 免费观看的影片在线观看| 亚洲国产精品国产精品| 国产探花极品一区二区| 日产精品乱码卡一卡2卡三| 亚洲欧洲日产国产| 毛片女人毛片| 如何舔出高潮| 国语自产精品视频在线第100页| 国产乱人偷精品视频| a级一级毛片免费在线观看| 国产精品.久久久| 亚洲国产精品sss在线观看| 欧美性感艳星| 可以在线观看的亚洲视频| 精品无人区乱码1区二区| 日本黄色片子视频| 久久精品国产亚洲网站| 美女 人体艺术 gogo| 国产精品嫩草影院av在线观看| 中国美女看黄片| 嫩草影院新地址| 在线观看美女被高潮喷水网站| 国产黄a三级三级三级人| 国产成人精品婷婷| 国产国拍精品亚洲av在线观看| 亚洲18禁久久av| 亚洲内射少妇av| 亚洲av一区综合| 国产黄片视频在线免费观看| 99国产精品一区二区蜜桃av| videossex国产| 18禁裸乳无遮挡免费网站照片| 男女视频在线观看网站免费| 美女内射精品一级片tv| 插阴视频在线观看视频| 国内精品宾馆在线| 观看免费一级毛片| 国产极品精品免费视频能看的| 黑人高潮一二区| 久久午夜福利片| 中出人妻视频一区二区| 亚洲精华国产精华液的使用体验 | videossex国产| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看| 精品午夜福利在线看| 日本一二三区视频观看| 成熟少妇高潮喷水视频| 亚洲在久久综合| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 久久精品夜色国产| 在线国产一区二区在线| 国产乱人视频| 中国国产av一级| 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 日产精品乱码卡一卡2卡三| 丰满乱子伦码专区| 最近的中文字幕免费完整| 人人妻人人看人人澡| av免费观看日本| 国产亚洲av嫩草精品影院| 成人综合一区亚洲| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 99热这里只有精品一区| 在线天堂最新版资源| 色哟哟·www| 99热这里只有是精品在线观看| 麻豆成人av视频| 国产精品av视频在线免费观看| 18禁裸乳无遮挡免费网站照片| 天堂影院成人在线观看| 我的老师免费观看完整版| 黑人高潮一二区| 午夜激情欧美在线| 欧美变态另类bdsm刘玥| 1000部很黄的大片| 波多野结衣高清作品| 久久久久性生活片| 在线国产一区二区在线| 村上凉子中文字幕在线| 日韩成人av中文字幕在线观看| 欧美另类亚洲清纯唯美| 欧美日韩一区二区视频在线观看视频在线 | 免费看av在线观看网站| 国产伦理片在线播放av一区 | 丰满乱子伦码专区| 亚洲欧美中文字幕日韩二区| 免费看光身美女| 久久久久久久久大av| 精品一区二区三区视频在线| 欧美日本视频| 日韩欧美一区二区三区在线观看| 国产亚洲精品久久久com| av专区在线播放| 亚洲精品日韩av片在线观看| 亚洲精品影视一区二区三区av| 免费人成在线观看视频色| 国产一区二区三区av在线 | 伊人久久精品亚洲午夜| 久久婷婷人人爽人人干人人爱| 欧美色视频一区免费| 波多野结衣巨乳人妻| 日产精品乱码卡一卡2卡三| 国产高清有码在线观看视频| 精品久久久久久久久久久久久| 欧美精品一区二区大全| 综合色av麻豆| 午夜激情福利司机影院| 少妇熟女aⅴ在线视频| 亚洲精品日韩av片在线观看| 91久久精品国产一区二区成人| 在线播放国产精品三级| 亚洲人成网站在线播放欧美日韩| 欧美一区二区国产精品久久精品| 秋霞在线观看毛片| 一级黄色大片毛片| 久久久久久久午夜电影| 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 又粗又硬又长又爽又黄的视频 | 18+在线观看网站| 亚洲精品日韩av片在线观看| or卡值多少钱| 国产精品久久久久久精品电影小说 | 国产男人的电影天堂91| 日韩在线高清观看一区二区三区| 99久久精品一区二区三区| 九九爱精品视频在线观看| 国产中年淑女户外野战色| av女优亚洲男人天堂| 亚洲欧美中文字幕日韩二区| 日韩欧美精品免费久久| 少妇猛男粗大的猛烈进出视频 | 国产一区二区三区在线臀色熟女| 中文字幕制服av| 亚洲一区二区三区色噜噜| 18禁在线播放成人免费| 欧美高清性xxxxhd video| 免费观看在线日韩| 91久久精品国产一区二区三区| 国产91av在线免费观看| 小说图片视频综合网站| 久久久久久国产a免费观看| 亚洲人成网站在线播| 男人舔奶头视频| av卡一久久| 欧美潮喷喷水| 国产一区二区在线观看日韩| 国产亚洲av嫩草精品影院| 国产免费男女视频| 美女黄网站色视频| 久久人人爽人人片av| 欧美日韩综合久久久久久| 国产精品精品国产色婷婷| 在线播放国产精品三级| 日日啪夜夜撸| 最后的刺客免费高清国语| 欧美高清成人免费视频www| 精品一区二区免费观看| 精品无人区乱码1区二区| 国产黄片美女视频| 亚洲无线在线观看| 亚洲国产精品国产精品| 国产精品综合久久久久久久免费| 嫩草影院入口| 国产欧美日韩精品一区二区| 3wmmmm亚洲av在线观看| 日韩 亚洲 欧美在线| 在线观看66精品国产| 综合色av麻豆| 国产高清三级在线| 亚洲乱码一区二区免费版| 1024手机看黄色片| 黄色一级大片看看| 成人漫画全彩无遮挡| 99热这里只有是精品50| 看非洲黑人一级黄片| 麻豆精品久久久久久蜜桃| 亚洲人成网站在线播放欧美日韩| 天美传媒精品一区二区| 日韩一区二区视频免费看| 亚洲熟妇中文字幕五十中出| 日韩av不卡免费在线播放| 免费搜索国产男女视频| 亚洲av中文字字幕乱码综合| 精品国产三级普通话版| 亚洲四区av| 成年版毛片免费区| 久久久精品欧美日韩精品| 精品久久久久久久久久免费视频| 日本色播在线视频| 在线免费观看不下载黄p国产| 人体艺术视频欧美日本| 亚洲精品国产成人久久av| 又爽又黄a免费视频| 久久国产乱子免费精品| 亚洲欧美中文字幕日韩二区| 在线a可以看的网站| 一级毛片久久久久久久久女| 黄色欧美视频在线观看| 91av网一区二区| 熟女人妻精品中文字幕| 国语自产精品视频在线第100页| 日韩强制内射视频| 国产又黄又爽又无遮挡在线| 久久精品久久久久久噜噜老黄 | 久久6这里有精品| 联通29元200g的流量卡| 国产成人91sexporn| 女人被狂操c到高潮| 一区二区三区高清视频在线| 少妇猛男粗大的猛烈进出视频 | 内射极品少妇av片p| 日韩欧美一区二区三区在线观看| 性欧美人与动物交配| 99久久久亚洲精品蜜臀av| 国产精品久久久久久av不卡| 91狼人影院| 国国产精品蜜臀av免费| 亚洲第一电影网av| av免费观看日本| 国产精品一二三区在线看| 欧美精品一区二区大全| 精品久久久久久久末码| 赤兔流量卡办理| 丰满人妻一区二区三区视频av| 极品教师在线视频| 午夜福利高清视频| 日韩欧美精品免费久久| 成年版毛片免费区| 搡老妇女老女人老熟妇| 日韩欧美在线乱码| 国产成人福利小说| 国内精品久久久久精免费| 久久久久久久久大av| 亚洲人成网站在线播| 亚洲自拍偷在线| 免费在线观看成人毛片| 色哟哟哟哟哟哟| 有码 亚洲区| 国产在线精品亚洲第一网站| 精品一区二区免费观看| 免费黄网站久久成人精品| 亚洲成人中文字幕在线播放| 少妇熟女欧美另类| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 人人妻人人澡人人爽人人夜夜 | 性欧美人与动物交配| 男女啪啪激烈高潮av片| 亚洲成人久久爱视频| 亚洲中文字幕一区二区三区有码在线看| 嫩草影院入口| 欧美不卡视频在线免费观看| 欧美激情在线99| 精品少妇黑人巨大在线播放 | 亚洲熟妇中文字幕五十中出| 国产色爽女视频免费观看| 久久久成人免费电影| 国内久久婷婷六月综合欲色啪| 男女做爰动态图高潮gif福利片| 久久久国产成人免费| 国产免费一级a男人的天堂| 国产老妇伦熟女老妇高清| 在线天堂最新版资源| 日日干狠狠操夜夜爽| 看免费成人av毛片| 岛国在线免费视频观看| 久久鲁丝午夜福利片| 日本av手机在线免费观看| 日本熟妇午夜| 狠狠狠狠99中文字幕| 黄色配什么色好看| 69av精品久久久久久| 亚洲性久久影院| 18禁在线播放成人免费| 国产伦一二天堂av在线观看| 深爱激情五月婷婷| 国产成人午夜福利电影在线观看| 九九热线精品视视频播放| 18+在线观看网站| 久久久精品欧美日韩精品| 中文资源天堂在线| 中国国产av一级| 尾随美女入室| 在线观看美女被高潮喷水网站| 欧美人与善性xxx| 日日干狠狠操夜夜爽| 黄色日韩在线| videossex国产| 你懂的网址亚洲精品在线观看 | 99riav亚洲国产免费| 全区人妻精品视频| 春色校园在线视频观看| 插逼视频在线观看| 卡戴珊不雅视频在线播放| 久久久久久久久久久免费av| 精品久久久久久久末码| 亚洲中文字幕日韩| 久久九九热精品免费| av在线蜜桃| 精品人妻熟女av久视频| 欧美+亚洲+日韩+国产| 午夜爱爱视频在线播放| 1000部很黄的大片| 国产黄色视频一区二区在线观看 | 国产三级中文精品| 美女 人体艺术 gogo| av天堂在线播放| 91精品国产九色| 99久久精品热视频| av在线老鸭窝| 床上黄色一级片| 97超碰精品成人国产| 校园人妻丝袜中文字幕| 黄色配什么色好看| 亚洲精品久久国产高清桃花| 色哟哟哟哟哟哟| 久久人人爽人人片av| 国产一区二区三区av在线 | 亚洲精品成人久久久久久| 午夜精品国产一区二区电影 | 久久99蜜桃精品久久| 亚洲av不卡在线观看| 日日撸夜夜添| 一级毛片aaaaaa免费看小| 久久精品综合一区二区三区| 亚洲,欧美,日韩| 亚洲高清免费不卡视频| 国产成人精品久久久久久| 能在线免费观看的黄片| 亚洲激情五月婷婷啪啪| 久久草成人影院| 国产蜜桃级精品一区二区三区| 日本一本二区三区精品| 狂野欧美激情性xxxx在线观看| 国产伦一二天堂av在线观看| 毛片一级片免费看久久久久| 国产熟女欧美一区二区|