• <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| 亚洲欧美清纯卡通| 一区二区三区乱码不卡18| 久久青草综合色| 欧美少妇被猛烈插入视频| 亚洲国产av影院在线观看| 亚洲天堂av无毛| 午夜福利在线观看免费完整高清在| 大香蕉久久网| 国产伦理片在线播放av一区| 日韩电影二区| 秋霞伦理黄片| 久久精品国产亚洲av天美| 久久久a久久爽久久v久久| 久久久国产欧美日韩av| 水蜜桃什么品种好| 一本久久精品| 丝袜美足系列| 国产成人a∨麻豆精品| 丰满迷人的少妇在线观看| 精品一区在线观看国产| 欧美精品国产亚洲| 国产一区二区在线观看日韩| 一区二区三区精品91| 国产精品一区二区在线观看99| 久久精品国产亚洲av天美| 免费黄色在线免费观看| 欧美 亚洲 国产 日韩一| 亚洲婷婷狠狠爱综合网| 午夜福利视频精品| 国产精品久久久久久av不卡| www.熟女人妻精品国产 | 另类精品久久| 一级黄片播放器| 97精品久久久久久久久久精品| 免费日韩欧美在线观看| 春色校园在线视频观看| 18禁在线无遮挡免费观看视频| 国产精品久久久久久av不卡| 精品少妇内射三级| 国产1区2区3区精品| 热re99久久国产66热| 嫩草影院入口| 亚洲欧美日韩卡通动漫| 在线观看免费高清a一片| 男人爽女人下面视频在线观看| 观看美女的网站| 久久97久久精品| 1024视频免费在线观看| 成人午夜精彩视频在线观看| 一区二区日韩欧美中文字幕 | 亚洲性久久影院| 免费人成在线观看视频色| 久久精品aⅴ一区二区三区四区 | 丝袜脚勾引网站| 免费日韩欧美在线观看| 又黄又粗又硬又大视频| 免费观看性生交大片5| 久久久久久伊人网av| 五月伊人婷婷丁香| 欧美日韩成人在线一区二区| 国产熟女午夜一区二区三区| 欧美亚洲 丝袜 人妻 在线| 观看av在线不卡| 视频区图区小说| 亚洲精品乱码久久久久久按摩| 国产黄色免费在线视频| 91午夜精品亚洲一区二区三区| 插逼视频在线观看| 午夜视频国产福利| 在线观看免费视频网站a站| 久久毛片免费看一区二区三区| 国产成人精品一,二区| 秋霞在线观看毛片| 美女国产高潮福利片在线看| 国产精品久久久久成人av| av在线播放精品| 在线观看美女被高潮喷水网站| 一级毛片电影观看| √禁漫天堂资源中文www| 国产成人精品婷婷| 亚洲av.av天堂| 日韩av免费高清视频| xxx大片免费视频| 国产男女内射视频| 国产精品秋霞免费鲁丝片| 伊人亚洲综合成人网| 精品福利永久在线观看| 天天操日日干夜夜撸| 婷婷色麻豆天堂久久| 蜜臀久久99精品久久宅男| 国产精品偷伦视频观看了| 久久99精品国语久久久| 国产女主播在线喷水免费视频网站| 国产av码专区亚洲av| 99视频精品全部免费 在线| 久久女婷五月综合色啪小说| 免费女性裸体啪啪无遮挡网站| 天天躁夜夜躁狠狠久久av| 人体艺术视频欧美日本| a级毛片在线看网站| 超色免费av| 男人操女人黄网站| 日本午夜av视频| 欧美bdsm另类| 精品一区二区三区四区五区乱码 | 97超碰精品成人国产| 久久这里有精品视频免费| 精品一品国产午夜福利视频| 亚洲国产看品久久| kizo精华| 亚洲欧美成人精品一区二区| 伦精品一区二区三区| 久久精品国产亚洲av涩爱| 亚洲性久久影院| 国产视频首页在线观看| 婷婷色麻豆天堂久久| 欧美精品一区二区免费开放| 精品久久久久久电影网| 少妇的丰满在线观看| 大码成人一级视频| 日本猛色少妇xxxxx猛交久久| 2018国产大陆天天弄谢| 黑人猛操日本美女一级片| 日韩av免费高清视频| 黑人欧美特级aaaaaa片| 搡老乐熟女国产| 精品一区在线观看国产| 久久久久久久精品精品| 亚洲高清免费不卡视频| 熟女人妻精品中文字幕| 美女视频免费永久观看网站| 精品亚洲成国产av| 18禁观看日本| 国产乱人偷精品视频| 国产亚洲精品第一综合不卡 | 欧美精品一区二区大全| 午夜福利网站1000一区二区三区| 寂寞人妻少妇视频99o| 亚洲精品aⅴ在线观看| 26uuu在线亚洲综合色| 看非洲黑人一级黄片| 2022亚洲国产成人精品| 国产精品一国产av| 日本wwww免费看| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧洲日产国产| 国产综合精华液| 9191精品国产免费久久| av不卡在线播放| 欧美国产精品va在线观看不卡| 高清欧美精品videossex| 下体分泌物呈黄色| 丝袜脚勾引网站| 最近2019中文字幕mv第一页| 黄网站色视频无遮挡免费观看| 免费看av在线观看网站| 飞空精品影院首页| 国产午夜精品一二区理论片| 久久人妻熟女aⅴ| 日本vs欧美在线观看视频| 人人妻人人添人人爽欧美一区卜| 亚洲欧美日韩另类电影网站| 又黄又爽又刺激的免费视频.| videosex国产| 国产一区二区在线观看日韩| 2018国产大陆天天弄谢| av在线播放精品| 2021少妇久久久久久久久久久| 99久久人妻综合| 看免费成人av毛片| 亚洲精品日韩在线中文字幕| 亚洲美女视频黄频| 9色porny在线观看| 美国免费a级毛片| 国产精品国产三级国产专区5o| 精品亚洲乱码少妇综合久久| 久久久欧美国产精品| 欧美最新免费一区二区三区| av在线观看视频网站免费| 一级a做视频免费观看| 欧美日韩国产mv在线观看视频| 99久国产av精品国产电影| 少妇熟女欧美另类| 国产亚洲精品第一综合不卡 | 蜜臀久久99精品久久宅男| 一本色道久久久久久精品综合| 成人漫画全彩无遮挡| 看免费av毛片| 一级毛片我不卡| 欧美xxxx性猛交bbbb| 婷婷色av中文字幕| 在线观看人妻少妇| 国产淫语在线视频| 亚洲国产日韩一区二区| 熟女电影av网| 日本与韩国留学比较| 精品一区二区三区四区五区乱码 | 伊人亚洲综合成人网| 美女主播在线视频| 一级毛片 在线播放| 丰满乱子伦码专区| 午夜福利视频精品| 国产黄色视频一区二区在线观看| 在线观看人妻少妇| 女人久久www免费人成看片| 亚洲美女黄色视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品久久精品一区二区三区| 97人妻天天添夜夜摸| 涩涩av久久男人的天堂| 免费在线观看完整版高清| 狠狠婷婷综合久久久久久88av| 免费人妻精品一区二区三区视频| 99热国产这里只有精品6| 中国国产av一级| 我要看黄色一级片免费的| 99久久人妻综合| 欧美激情国产日韩精品一区| 天天躁夜夜躁狠狠久久av| 国产精品免费大片| www日本在线高清视频| 欧美日韩国产mv在线观看视频| 热99国产精品久久久久久7| 黄色 视频免费看| 中国美白少妇内射xxxbb| 国产女主播在线喷水免费视频网站| av免费在线看不卡| 久久久久久久精品精品| 99热6这里只有精品| 午夜免费男女啪啪视频观看| 在线看a的网站| 亚洲伊人久久精品综合| 国产高清国产精品国产三级| 久久国内精品自在自线图片| 久久精品国产鲁丝片午夜精品| 久久久久久久久久人人人人人人| 国产av码专区亚洲av| 国产精品一国产av| 国产男女内射视频| 黄网站色视频无遮挡免费观看| 欧美人与善性xxx| 亚洲av中文av极速乱| 国国产精品蜜臀av免费| 国产精品嫩草影院av在线观看| 伊人久久国产一区二区| av有码第一页| 如日韩欧美国产精品一区二区三区| 黑人猛操日本美女一级片| 国产片特级美女逼逼视频| 国产色婷婷99| 亚洲精品aⅴ在线观看| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| 亚洲图色成人| 伦理电影大哥的女人| 超色免费av| 欧美国产精品一级二级三级| 国产片特级美女逼逼视频| 国产亚洲精品第一综合不卡 | 少妇的丰满在线观看| 日本欧美国产在线视频| 国产日韩欧美亚洲二区| 精品熟女少妇av免费看| 极品人妻少妇av视频| 免费观看无遮挡的男女| 国产成人欧美| 欧美激情国产日韩精品一区| 亚洲精品乱码久久久久久按摩| 日韩中文字幕视频在线看片| 99九九在线精品视频| 99久久精品国产国产毛片| 秋霞在线观看毛片| 极品少妇高潮喷水抽搐| 97在线人人人人妻| 久久国产精品大桥未久av| 中国三级夫妇交换| 日本欧美国产在线视频| 久热这里只有精品99| av播播在线观看一区| 九色成人免费人妻av| 丝袜脚勾引网站| 丝袜喷水一区| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 9色porny在线观看| 国产精品一二三区在线看| 国产伦理片在线播放av一区| 亚洲精品自拍成人| 中文字幕最新亚洲高清| 青春草国产在线视频| 中国美白少妇内射xxxbb| 久久99蜜桃精品久久| 91久久精品国产一区二区三区| 最近的中文字幕免费完整| av视频免费观看在线观看| 人成视频在线观看免费观看| 又粗又硬又长又爽又黄的视频| 欧美国产精品va在线观看不卡| 亚洲丝袜综合中文字幕| 欧美精品一区二区大全| 欧美国产精品一级二级三级| 久久久久国产精品人妻一区二区| 精品视频人人做人人爽| 亚洲激情五月婷婷啪啪| 亚洲av电影在线观看一区二区三区| 青春草国产在线视频| 91精品三级在线观看| 亚洲国产精品专区欧美| 午夜激情久久久久久久| 久久久久网色| 午夜福利乱码中文字幕| 精品国产露脸久久av麻豆| 日韩三级伦理在线观看| 亚洲国产最新在线播放| 香蕉精品网在线| 伦理电影大哥的女人| 精品国产一区二区久久| 欧美国产精品va在线观看不卡| 国产色爽女视频免费观看| 日本-黄色视频高清免费观看| 777米奇影视久久| 亚洲综合精品二区| 日本91视频免费播放| 免费观看av网站的网址| 国产亚洲欧美精品永久| 欧美人与性动交α欧美软件 | 伊人亚洲综合成人网| 男女国产视频网站| 丰满饥渴人妻一区二区三| 深夜精品福利| 哪个播放器可以免费观看大片| 高清不卡的av网站| 欧美亚洲日本最大视频资源| 黄色配什么色好看| 亚洲欧美成人综合另类久久久| 丝袜美足系列| 日日啪夜夜爽| 久久女婷五月综合色啪小说| 少妇人妻久久综合中文| 亚洲 欧美一区二区三区| 九九爱精品视频在线观看| 99精国产麻豆久久婷婷| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产精品成人久久小说| 在线观看免费视频网站a站| 男女边摸边吃奶| av天堂久久9| 日本欧美国产在线视频| 巨乳人妻的诱惑在线观看| 2022亚洲国产成人精品| 国产淫语在线视频| 成人毛片a级毛片在线播放| 精品久久蜜臀av无| 国产成人aa在线观看| 尾随美女入室| 日韩伦理黄色片| 国产精品熟女久久久久浪| 欧美老熟妇乱子伦牲交| 美女主播在线视频| 欧美国产精品va在线观看不卡| 久久久精品免费免费高清| av在线观看视频网站免费| 一区二区三区四区激情视频| 亚洲av在线观看美女高潮| 一级毛片黄色毛片免费观看视频| 少妇的逼好多水| 热re99久久精品国产66热6| 久久精品人人爽人人爽视色| 少妇人妻久久综合中文| 美女国产高潮福利片在线看| 深夜精品福利| 国产欧美日韩综合在线一区二区| 咕卡用的链子| 日韩欧美一区视频在线观看| 亚洲激情五月婷婷啪啪| 久久午夜综合久久蜜桃| 内地一区二区视频在线| 一级毛片黄色毛片免费观看视频| 欧美国产精品一级二级三级| 久久99热6这里只有精品| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 久久韩国三级中文字幕| 国产亚洲精品第一综合不卡 | 久久久国产一区二区| a级毛片黄视频| 午夜福利,免费看| 日韩制服骚丝袜av| 亚洲综合精品二区| 少妇被粗大的猛进出69影院 | 亚洲情色 制服丝袜| 午夜精品国产一区二区电影| 性高湖久久久久久久久免费观看| 搡女人真爽免费视频火全软件| 99久久精品国产国产毛片| 精品一品国产午夜福利视频| 日韩av在线免费看完整版不卡| 日韩av免费高清视频| 啦啦啦啦在线视频资源| 色94色欧美一区二区| 中文精品一卡2卡3卡4更新| 午夜福利影视在线免费观看| 国产成人精品一,二区| 精品福利永久在线观看| 伦理电影免费视频| 香蕉精品网在线| 美女脱内裤让男人舔精品视频| 一本色道久久久久久精品综合| 亚洲五月色婷婷综合| 久久久精品免费免费高清| 亚洲美女视频黄频| 国产精品久久久久久精品电影小说| 欧美激情国产日韩精品一区| 七月丁香在线播放| 国产精品国产三级国产专区5o| 国产免费现黄频在线看| √禁漫天堂资源中文www| 美女视频免费永久观看网站| 欧美日韩综合久久久久久| xxx大片免费视频| 99热这里只有是精品在线观看| 亚洲av电影在线进入| 国产片内射在线| 国产在线一区二区三区精| 日韩av在线免费看完整版不卡| 久久韩国三级中文字幕| 亚洲三级黄色毛片| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 国产男女超爽视频在线观看| 亚洲情色 制服丝袜| 天天影视国产精品| 人人妻人人澡人人看| 九色亚洲精品在线播放| 看十八女毛片水多多多| 男女边吃奶边做爰视频| 国产欧美日韩一区二区三区在线| 天堂中文最新版在线下载| 免费大片18禁| 免费少妇av软件| 哪个播放器可以免费观看大片| 一本大道久久a久久精品| 欧美人与善性xxx| 99久久中文字幕三级久久日本| 国产精品久久久久久av不卡| 久久精品熟女亚洲av麻豆精品| 免费高清在线观看视频在线观看| 中文字幕精品免费在线观看视频 | 亚洲成色77777| 热99久久久久精品小说推荐| 国产精品无大码| 高清欧美精品videossex| 国产不卡av网站在线观看| 成人国产麻豆网| 人妻少妇偷人精品九色| 少妇高潮的动态图| 男女边吃奶边做爰视频| 自线自在国产av| 男女边摸边吃奶| av不卡在线播放| 久久精品人人爽人人爽视色| 亚洲av电影在线观看一区二区三区| 一个人免费看片子| 国产一区二区在线观看日韩| 人妻少妇偷人精品九色| 黑丝袜美女国产一区| 制服人妻中文乱码| 国产深夜福利视频在线观看| 一区二区三区精品91| 国产片特级美女逼逼视频| videosex国产| 亚洲国产精品一区三区| 欧美精品高潮呻吟av久久| 99热国产这里只有精品6| 日韩 亚洲 欧美在线| videosex国产| 在线观看国产h片| 成年女人在线观看亚洲视频| 国产精品久久久久成人av| 国产xxxxx性猛交| 国产精品一国产av| 人妻少妇偷人精品九色| 国产av精品麻豆| 欧美精品av麻豆av| 伊人久久国产一区二区| 国产极品天堂在线| 人妻人人澡人人爽人人| 色5月婷婷丁香| 爱豆传媒免费全集在线观看| 国产无遮挡羞羞视频在线观看| 精品熟女少妇av免费看| 午夜免费鲁丝| 日韩精品免费视频一区二区三区 | 久久午夜综合久久蜜桃| 成年女人在线观看亚洲视频| www.色视频.com| 热re99久久精品国产66热6| 国产免费一级a男人的天堂| 在线观看免费高清a一片| 亚洲在久久综合| 在线观看国产h片| h视频一区二区三区| 亚洲综合色惰| 丝瓜视频免费看黄片| 91国产中文字幕| 少妇 在线观看| 国产女主播在线喷水免费视频网站| 日韩欧美精品免费久久| 伦理电影免费视频| av黄色大香蕉| 欧美bdsm另类| 另类精品久久| 最新的欧美精品一区二区| 中文乱码字字幕精品一区二区三区| 亚洲人成77777在线视频| 亚洲,一卡二卡三卡| 高清视频免费观看一区二区| 国国产精品蜜臀av免费| 国产日韩欧美在线精品| 美女国产视频在线观看| 欧美成人午夜免费资源| 日韩中字成人| 国产一级毛片在线| 免费人妻精品一区二区三区视频| 一二三四在线观看免费中文在 | 成人手机av| 黄色视频在线播放观看不卡| 国产福利在线免费观看视频| 亚洲精品国产av成人精品| 久久久国产精品麻豆| 午夜福利视频精品| www.av在线官网国产| av天堂久久9| 国产午夜精品一二区理论片| 一级片'在线观看视频| 日本91视频免费播放| 十分钟在线观看高清视频www| 啦啦啦啦在线视频资源| av在线观看视频网站免费| 少妇 在线观看| 日韩一本色道免费dvd| 中文乱码字字幕精品一区二区三区| 午夜福利在线观看免费完整高清在| videossex国产| 女人久久www免费人成看片| 黄色毛片三级朝国网站| 亚洲国产日韩一区二区| 久久精品久久精品一区二区三区| 久久久久久久久久成人| 观看美女的网站| 国产欧美亚洲国产| 亚洲精品av麻豆狂野| 国产av精品麻豆| 女性被躁到高潮视频| 久久免费观看电影| 国产乱来视频区| 两性夫妻黄色片 | 国产午夜精品一二区理论片| 国产男女内射视频| 日韩电影二区| 色哟哟·www| 一区二区日韩欧美中文字幕 | 国产老妇伦熟女老妇高清| 午夜福利网站1000一区二区三区| 国产精品熟女久久久久浪| 午夜福利在线观看免费完整高清在| 又黄又粗又硬又大视频| 免费观看av网站的网址| av不卡在线播放| 亚洲一码二码三码区别大吗| 中文乱码字字幕精品一区二区三区| 永久网站在线| 女性被躁到高潮视频| 国产女主播在线喷水免费视频网站| 女性被躁到高潮视频| 国产女主播在线喷水免费视频网站| 成年人午夜在线观看视频| 母亲3免费完整高清在线观看 | 中文字幕最新亚洲高清| 亚洲精品美女久久久久99蜜臀 | 在线免费观看不下载黄p国产| 国产黄色免费在线视频| 亚洲国产欧美日韩在线播放| 国产成人91sexporn| 成年女人在线观看亚洲视频| 黑人欧美特级aaaaaa片| 欧美激情极品国产一区二区三区 | 黄色配什么色好看| videossex国产| 亚洲av日韩在线播放| 久久人人爽人人片av| 五月伊人婷婷丁香| 日日啪夜夜爽| 亚洲精品乱久久久久久| 蜜桃在线观看..| 91在线精品国自产拍蜜月| 九草在线视频观看| av国产久精品久网站免费入址| 欧美日韩综合久久久久久| 春色校园在线视频观看| 永久网站在线| 国产黄色视频一区二区在线观看| freevideosex欧美| 97精品久久久久久久久久精品| 免费高清在线观看视频在线观看| 国产无遮挡羞羞视频在线观看| 欧美日韩精品成人综合77777| 国产欧美日韩一区二区三区在线| 成人国产麻豆网| av在线老鸭窝| 欧美少妇被猛烈插入视频| 伦精品一区二区三区| 嫩草影院入口|