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

    南海北部內(nèi)孤立波生成對黑潮響應的初步研究

    2016-07-06 10:30:20郭大全中國科學院海洋研究所山東青島266071中國科學院大學北京100049中國科學院海洋環(huán)流與波動重點實驗室山東青島266071
    海洋科學 2016年2期
    關鍵詞:黑潮內(nèi)波潮流

    郭大全, 申 輝(1. 中國科學院海洋研究所, 山東 青島 266071; 2. 中國科學院大學, 北京 100049; 3. 中國科學院 海洋環(huán)流與波動重點實驗室, 山東 青島 266071)

    ?

    南海北部內(nèi)孤立波生成對黑潮響應的初步研究

    郭大全1, 2, 3, 申 輝1, 3
    (1. 中國科學院海洋研究所, 山東 青島 266071; 2. 中國科學院大學, 北京 100049; 3. 中國科學院 海洋環(huán)流與波動重點實驗室, 山東 青島 266071)

    為了研究南海北部黑潮入侵對內(nèi)孤立波生成造成的影響, 本研究在對合成孔徑雷達圖像處理分析的基礎上, 運用二維非靜力模式, 對南海北部內(nèi)孤立波生成以及其對黑潮入侵的響應做了初步的模擬和分析。對4 a(2007~2010年)的合成孔徑雷達(SAR)圖像統(tǒng)計分析發(fā)現(xiàn), 內(nèi)波出現(xiàn)頻率呈現(xiàn)明顯的季節(jié)變化。月平均數(shù)據(jù)表明, 內(nèi)波出現(xiàn)頻率在5月到8月期間較高并在6月達到峰值, 占全年總數(shù)的29%;而在11月到次年的2月份期間出現(xiàn)頻率較低并在12月和1月達到最低, 各占全年的0.5%。通過二維、非靜力近似數(shù)值模擬研究發(fā)現(xiàn): 黑潮入侵引起的平流對呂宋海峽區(qū)域內(nèi)孤立波的激發(fā)產(chǎn)生具有明顯的影響。具體表現(xiàn)為, 黑潮入侵引起的的西向平均流拉伸了內(nèi)孤立波波形并加強內(nèi)孤立波的能量的頻散,使得大振幅內(nèi)孤立波演變成多個分散的波動, 波長變長、振幅變小。敏感性實驗通過替換夏季溫鹽場為冬季溫鹽場, 對模擬結果并無實質(zhì)上的影響, 表明冬夏季溫鹽場的差異并非是SAR觀測內(nèi)波出現(xiàn)頻率呈季節(jié)變化的重要因素。

    內(nèi)波; 合成孔徑雷達圖像; 二維非靜力近似; MITgcm模式; 黑潮

    [Foundation: National Natural Science Foundation of China, No.41030855]

    大振幅海洋內(nèi)孤立波廣泛分布于南海北部, 以其非線性強、振幅大而受到廣泛關注。2001年的亞洲海洋聲學實驗在南海北部開展了較大規(guī)模的海洋內(nèi)孤立波現(xiàn)場觀測。衛(wèi)星遙感和數(shù)值模擬技術是現(xiàn)場觀測的重要補充, 在海洋內(nèi)孤立波時空分布特征研究和機制分析等方面具有獨到的優(yōu)勢。合成孔徑雷達被廣泛地應用于該海域的內(nèi)孤立波研究中, 如:直布羅陀海峽、蘇祿海和南中國海等[1-2]。Zheng等[3]在2007年利用收集到的1995~2001年的合成孔徑雷達(SAR)圖像統(tǒng)計分析得出, 中國南海北部內(nèi)孤立波的出現(xiàn)在時間和空間上呈現(xiàn)一定的規(guī)律性, 有著明顯的年際變化和季節(jié)變化, 并指出SAR觀測受海表面狀況如風速等影響很大, 因此海表面狀況成為修正SAR觀測內(nèi)波出現(xiàn)頻率的主要外部條件, 而黑潮作為呂宋海峽東部的強勁西邊界流是導致其呈現(xiàn)變化的另一個關鍵因素, 對此可能扮演著至關重要的角色。Cai等[4]在2002年用一個兩層模式, 解釋了中國南海的內(nèi)孤立波產(chǎn)生于一種類似于lee-wave的機制。Yuan等[5]在2006年開發(fā)了一個線性的模式用來分析黑潮的不穩(wěn)定性, 結果表明黑潮西分支不穩(wěn)定而東分支穩(wěn)定, 暗示了黑潮的西分支是呂宋海峽內(nèi)波的擾動源。Du等[6]在2008年用一個二維靜力模式研究了包括潮流和象征著黑潮平均流在內(nèi)的內(nèi)波產(chǎn)生機制, 表明當西邊界流在山脊處強于2 m/s時會明顯地加強海水的混合, 類似于Maxworthy[7], 這種強混合區(qū)的塌陷會導致西向傳播非線性內(nèi)孤立波的產(chǎn)生。Wang等[8]在2010年用一個二維非靜力模式模擬了呂宋海峽處內(nèi)孤立波的產(chǎn)生與傳播, 他們認為黑潮入侵南海對內(nèi)孤立波的產(chǎn)生并不十分重要。在加入了垂向變化的流動之后, 內(nèi)孤立波的產(chǎn)生形式?jīng)]有明顯的改變。

    對于南海北部內(nèi)孤立波生成的黑潮響應雖然已有以上的研究, 但黑潮分支究竟是如何影響內(nèi)孤立波的生成及其影響方式與機制尚無確定的結論。本文將首先通過分析大量海洋內(nèi)孤立波遙感影響獲取南海北部海洋內(nèi)孤立波季節(jié)變化特征, 進而通過數(shù)值模擬手段探討黑潮入侵引起的平流對通過呂宋海峽口激發(fā)生成的海洋內(nèi)孤立波的影響。

    1 基于SAR的南海北部內(nèi)孤立波季節(jié)變化特征

    目前對南海內(nèi)孤立波的研究主要集中在南海北部, 研究方法主要有現(xiàn)場觀測、衛(wèi)星遙感以及數(shù)值模擬研究。合成孔徑雷達(SAR)作為傳統(tǒng)的遙感觀測手段因其空間覆蓋范圍廣、空間分辨率高、不受惡劣天氣狀況影響以及資料獲取費用相對較低等優(yōu)點,對研究內(nèi)孤立波有得天獨厚的優(yōu)勢[9]。

    我們對香港衛(wèi)星地面站(CUHK)存檔的(2007~ 2010年)共計16 037幅SAR圖像開展了篩選與統(tǒng)計分析, 獲得捕捉到海洋內(nèi)孤立波信息的SAR圖像共計150幅。各年月的內(nèi)波出現(xiàn)頻率分布如圖1所示??梢钥闯? 2007年4、5、6月內(nèi)孤立波出現(xiàn)頻率較高, 6月最高占全年的21.5%, 12、1、2月出現(xiàn)較少, 2月份最少為1.3%; 2008年5、6、7月內(nèi)波出現(xiàn)頻率較高, 7月最高占全年的25%, 1、2、8、9、11、12月出現(xiàn)頻率都很低為2.8%; 2009年6、7、8月出現(xiàn)頻率較高, 6月最高為23%, 1、4、9、11、12月出現(xiàn)頻率都很低, 其中9月份出現(xiàn)頻率為0, 沒有一次觀測到的內(nèi)孤立波出現(xiàn); 2010年只有1月份到10月份的數(shù)據(jù), 6、7、8月出現(xiàn)內(nèi)波頻率較高, 6月最高為32%, 1、3、9、10月都很低, 其中10月份內(nèi)波出現(xiàn)頻率為0。

    圖1 CUHK 存檔的SAR觀測內(nèi)波出現(xiàn)頻率月份變化Fig. 1 Monthly SAR-observed internal solitary waves (ISWs) occurrence frequencies from CUHK

    通過對4 a(2007~2010年)平均的合成孔徑雷達(SAR)圖像統(tǒng)計分析發(fā)現(xiàn), 內(nèi)波出現(xiàn)頻率呈現(xiàn)明顯的季節(jié)變化。如圖2所示, 4 a(2007~2010年)月平均數(shù)據(jù)表明: 內(nèi)波出現(xiàn)頻率在5月到8月期間較高并在6月達到峰值, 占全年總數(shù)的29%; 而在11月到次年的2月份出現(xiàn)頻率較低并在十二月和一月達到最低,各占全年的0.5%。這一統(tǒng)計結果與Zheng等[3]對1995~2001共計7年的SAR圖像月平均統(tǒng)計結果相符, 表明內(nèi)波的出現(xiàn)頻率呈現(xiàn)出明顯的季節(jié)變化:夏季出現(xiàn)頻率最高, 春秋次之, 冬季最低。

    與此同時, 黑潮入侵南海存在多模態(tài)和季節(jié)變化特征。前人[10-11]研究表明, 呂宋海峽處的黑潮入侵也呈現(xiàn)明顯的季節(jié)變化特征, 主要表現(xiàn)為: 11月至來年3月為強盛期, 黑潮大規(guī)模入侵南海; 3月到6月為衰退期, 黑潮影響范圍逐漸北移直到退回呂宋海峽東側; 6月至9月, 黑潮入侵南海最弱, 沒有黑潮南海分支出現(xiàn); 9月到11月為成長期, 黑潮影響范圍逐漸向南移直到深入南海。

    圖2 CUHK 2007~2010年SAR觀測內(nèi)波出現(xiàn)頻率月平均變化Fig. 2 Monthly averaged SAR-observed ISWs occurrence frequencies from 2007 to 2010

    表1 黑潮入侵南海與SAR觀測內(nèi)波出現(xiàn)頻率季節(jié)變化對比Tab.1 Seasonal variability of SAR-observed IWs occurrence frequencies and intrusion of K-branch

    對表1的結果比較得出: 黑潮入侵南海與SAR觀測到的內(nèi)波出現(xiàn)頻率呈現(xiàn)極大的相關性, 當黑潮入侵南海最強盛時, SAR觀測內(nèi)波出現(xiàn)頻率最低; 而當黑潮入侵南海最少處于消亡期時, SAR觀測到的內(nèi)波出現(xiàn)頻率卻最高。

    盡管存在其他外部條件例如風速等對SAR觀測的影響, 但不可否認的是黑潮分支對南海內(nèi)孤立波的形成有著重要的影響。

    2 模式介紹與設置

    2.1 MITgcm模式

    MITgcm是美國麻省理工學院開發(fā)的基于原始方程的三維大氣-海洋通用環(huán)流數(shù)值模式[12]。模式可采用靜力近似、準靜力近似和非靜力近似等模塊, 可以模擬從小尺度(如對流過程)到大尺度(如全球環(huán)流)的各種現(xiàn)象; 采用有限體積法(finite volume method),使得方程的離散比較直觀; 采用正交曲線網(wǎng)格和網(wǎng)格切削技術, 可以很好地處理不規(guī)則地形和邊界。并且模式提供了成熟的并行計算方案, 可以廣泛的應用在各種平臺上。

    針對內(nèi)孤立波非靜力特性的考慮, 故采用該模式進行內(nèi)孤立波生成的實驗模擬。

    2.2 模擬區(qū)域地形設置

    呂宋海峽地形環(huán)境復雜, 在臺灣-呂宋海溝東側從北向南依次有巴坦群島、巴布延群島和許多小的海脊; 西側是相對高度較低的海脊。復雜的海底地形與潮流相互作用使其成為中國南海北部內(nèi)波生成的重要源地。

    在此區(qū)域產(chǎn)生的內(nèi)波傳播方向有兩種: 一種是沿著巴坦群島至東沙群島方向西向傳播; 另一種是自巴布延群島像西北方向傳播[6]。本文主要模擬西向傳播的內(nèi)孤立波, 數(shù)值計算區(qū)域如圖3所示: 區(qū)域范圍為117.9°~122.6°E, 20°~21°N, 采用的地形數(shù)據(jù)來自Etopo1 Global Relief Model, 分辨率為1’×1’。

    圖3 所選實驗海區(qū)地形圖Fig. 3 Topography of simulation domain

    如圖4所示, 為了更好地反映整個區(qū)域的地形走勢, 二維模擬將以上模擬區(qū)域做經(jīng)向平均, 藍色細線為20°~21°N間每隔1 000 m做一次的地形剖面,黑色實線為區(qū)域徑向平均地形; 此外, 本文關注的是內(nèi)孤立波的生成過程, 對平均后的地形做進一步的簡化, 將120.5°E以西的海域水深設置為–3 000 m, 將122.2°E以東的海域水深設置為–3 000 m, 兩海脊中間低于–3 000 m區(qū)域設為–3 000 m(圖中綠色實線標記)。模式采用的地形(黑色實線加綠色實線)最深深度為–3 000 m, 擁有A、B兩個海脊, A海脊最高處水深–490 m, B海脊最高處水深–1 650 m。

    圖4 二維模擬所用地形Fig. 4 Topography used in 2-D simulation

    2.3 溫鹽數(shù)據(jù)

    模式采用水平均一的海水分層, 溫度鹽度數(shù)據(jù)選自World Ocean Atlas, 數(shù)值模擬時將所得數(shù)據(jù)進一步插值到間距為10 m的格點上。圖5所示是模式采用的多點平均溫鹽數(shù)據(jù)和浮力頻率N的垂直分布。圖中實線描述的是夏季(6月份到8月份)情況, 虛線描述的是冬季(11月份到第二年1月份)情況。

    所選區(qū)域多點平均的溫鹽數(shù)據(jù)及計算出的N顯示:冬夏兩季海水溫鹽并無太大的差異, 在水深250 m以下幾乎一致。而差異只體現(xiàn)在上層海水, 具體體現(xiàn)為夏季海水的鹽度最小值比冬季約小0.3, 溫度最大值比冬季約大4度; 上層海水的溫鹽差異導致了浮力頻率的差異, 夏季海水浮力頻率最大值為0.017 8, 位于水深60 m左右; 冬季海水浮力頻率最大值為0.014 1, 位于水深80 m左右。

    圖5 模式所用溫度(T)、鹽度(S)和浮力頻率(N)垂直剖面圖Fig. 5 Profiles of temperature (T), salinity (S), and buoyancy frequency (N)

    2.4 網(wǎng)格及參數(shù)處理

    二維模擬x方向范圍為117.9°~122.6°E, 分為1 960個格點, 水平分辨率為250 m, y方向1 000 m,垂直方向共分155層, 1 500 m以上每10 m一層,1 500~3 000 m每300 m一層。為了防止正壓波動及斜壓波動在邊界處的反射, 在兩邊界處加了海綿區(qū)域, 使x方向的網(wǎng)格點間距從中間區(qū)域的250 m逐漸增大到邊界處的100 km。模式的時間步長設置為10 s以滿足CFL(Courant-Friedrichs-Lewy)條件。

    根據(jù)Vlasenko等[13]中的方法, 模式采用周期性的正壓潮流驅(qū)動, 選擇M2分潮, 水平流速振幅為U。具體體現(xiàn)為在緯向和經(jīng)向的動量平衡方程中分別加一個外力驅(qū)動項和, 其中對應水平流速振幅U所在位置的海水深度, H( x , y )為不同(x , y )點對應的水深, f為科氏力參數(shù), σ為對應分潮的角速率。此處選擇為區(qū)域平均緯度20.5°N, 因此計算得到科氏力參數(shù)f為5.0936×10–5。

    2.5 實驗設計

    針對黑潮入侵對內(nèi)孤立波的產(chǎn)生所造成的影響,同時考慮冬夏季海水溫鹽差異, 二維模擬設計了如下4組實驗(表2)。

    表2 二維數(shù)值模擬實驗設計Tab. 2 Design of 2-D numerical model

    標準實驗a選自夏季海水溫鹽條件(圖6所示為其初始場), M2分潮潮流最大流速U0為0.1 m/s, 模式初始時刻, 全場為靜止態(tài), 驅(qū)動潮流從0開始, 隨之為東向漸增的退潮潮流并作周期性往復。對比試驗b在實驗a的基礎上, 于東邊界疊加類似黑潮分支平均流(K-branch)。黑潮平均流速度(Uk)大小設置為0.2 m/s, 方向自東向西。對比實驗c將海水溫鹽場替換為冬季, 其他設置與實驗a相同, 用來模擬冬季海水層化條件下的內(nèi)波產(chǎn)生過程。對比實驗d在實驗c基礎上疊加類似黑潮分支平均流, 大小為0.2 m/s,方向自東向西。

    3 實驗結果及分析

    為模擬內(nèi)孤立波產(chǎn)生過程, 時間步長選擇10 s,運行9 360步, 模擬時間為26 h。如圖7a 中6幅圖像所示, 標準實驗a采用的是夏季海水分層, 不加黑潮分支平均流, 較好地模擬出了內(nèi)孤立波群的產(chǎn)生過程。

    圖6 實驗初始溫度場Fig. 6 Initial temperature field

    在實驗初始時刻, 東向漸增的退潮潮流使得海脊A西側潮流有向上的通量, 而在東側有向下的通量, 這會使得對應的海脊A西側的等溫線抬升而東側的等溫線下壓(圖7a, t =2 h); 隨后東向退潮潮流在3.1 h達到最大, 之后開始減小, 到6.2 h為0轉(zhuǎn)向為西向漸增的漲潮潮流, 之前海脊A西側抬升的等溫線由于沒有東向潮流能量的繼續(xù)供給, 反之還有西向的漲潮潮流的逐漸增大, 使得抬升的等溫線以整體的形式開始向下塌陷, 最大塌陷達到150 m, 并伴隨著西向漲潮潮流向西傳播(圖7a, t = 7, t =12 h); 隨后在12.42 h, 一個M2潮周期完成, 潮流重新變?yōu)闁|向漸增的退潮流, 巨大的塌陷由于慣性繼續(xù)向西傳播, 但同時因退潮流的阻礙, 塌陷開始變得平滑, 深度也有所減弱, 最大塌陷為70 m左右(如圖7a , t =17 h所示); 隨后在18.6 h潮流轉(zhuǎn)向, 需要指出的是, 海脊B的位置位于121°E左右(圖6紅色虛線標出), 由于海脊B的阻擋, 以及西向漸增的漲潮潮流的共同作用,使得之前較平滑的塌陷重新變得陡峭起來(圖7a, t =22 h); 最后西向傳播的塌陷經(jīng)過了地形的阻礙, 開始釋放變?yōu)榉€(wěn)定的內(nèi)孤立波列向西傳播(圖7a, t =26 h)。

    如圖7b, 6幅圖像所示, 實驗b在標準實驗a的基礎上疊加了黑潮分支平均流, 大小為0.2 m/s, 方向自東向西, 使得整個流場添加了西向的背景平均流場。在實驗初始時刻, 西向的平均流場使得海脊A東側等溫線抬升、西側等溫線下壓, 造成西側等溫線的塌陷(圖7b, t =2時), 但也正是由于平均流持續(xù)向西, 使得這種塌陷在并沒有獲得太多的能量聚集起來時便隨之向西傳播(圖7b, t = 7, 12 h), 盡管在121°E處, 由于海脊B的阻擋, 再次獲得能量的聚集,但同時由于海脊B較低及西向平均流均一、較快, 并無太明顯的塌陷加深, 并且隨著西向傳播而波形越來越小, 在實驗結束時, 初始形成的第一個波動已傳播到119°E以西, 振幅已經(jīng)減小到10 m左右; 在整個過程中, 雖然有西向平均流的干擾, 潮流變化所引起的塌陷仍然很明顯, 如圖7b, 2, 7, 12 h, 由于東向退潮流漸增再減小到0轉(zhuǎn)向, 再到西向漸增的西向漲潮潮流, 使得海脊A西側仍然形成類似實驗a中的等溫線抬升以及塌陷, 只不過由于背景平均流的作用, 將這種塌陷的能量通過類似頭波傳播的形式分散掉了, 導致塌陷在121°E處與海脊B相互作用時已經(jīng)沒有太過集中的能量與抖動, 最終形成的波動能量較小且不集中。

    圖7 內(nèi)孤立波形成過程Fig. 7 Generation of IWs

    在實驗結束時刻即26時, 一組內(nèi)孤立波群清晰地出現(xiàn)在120°E左右(圖8a左側黑框), 模擬出的波群與我們在2012年6月21日采集到的SAR圖像顯示的內(nèi)孤立波群有極大的吻合度(圖10黑色虛線框內(nèi)顯示): 一組自西向東傳播的下凹型內(nèi)孤立波; 頭波波長較長、振幅較大; 尾隨的波動波長與振幅依次減小。

    而對照實驗b同樣采用夏季海水分層, 相比實驗a增加了自東向西的類似黑潮分支平均流, 在實驗結束時刻即26 h, 并無清晰的內(nèi)孤立波群出現(xiàn),只是有振幅相對較小而波長較長的波動, 實際上,更像是a圖像波動的拉伸版本: 能量分散, 波長拉伸變長, 振幅減小(圖8左側黑框所示)。同時, 黑潮平均流也加劇了與海脊A的相互作用(圖8右側黑框所示), 相比實驗a, 實驗b在海脊A處等溫線的抖動更加劇烈復雜。

    圖8 a、b實驗結束時刻26 h溫度場對比圖Fig. 8 Comparison of temperature fields for cases a and b at t = 26 h

    圖9 Radarsat-2 采集合成孔徑雷達圖像Fig. 9 SAR image from Radarsat-2

    對比實驗c和實驗d所選為冬季海水溫鹽場, 其層化形式絕大部分都一致, 差異只體現(xiàn)在上150 m 處: 相比夏季溫鹽場, 冬季表層水最高溫度比夏季低, 最低鹽度比夏季高, 浮力頻率相對夏季也較低。由圖10可以看出溫鹽場的改變對內(nèi)孤立波群的形成影響不大, 其形成過程與實驗a過程類似, 最終的內(nèi)孤立波波形、波數(shù)、振幅等都與標準實驗a一致。而對比實驗d與實驗b相比同樣差別不大。由此表明溫鹽差異并非內(nèi)波出現(xiàn)頻率呈現(xiàn)季節(jié)變化的關鍵因素[13]。

    4 結論

    通過對合成孔徑雷達圖像統(tǒng)計分析發(fā)現(xiàn), SAR觀測內(nèi)波出現(xiàn)頻率呈現(xiàn)明顯的季節(jié)變化: 夏季出現(xiàn)頻率最高, 春秋次之, 冬季最低。

    二維、非靜力近似、半實際地形數(shù)值模擬較好地模擬出了內(nèi)孤立波群的產(chǎn)生過程, 并顯示: 在疊加了類似黑潮分支平均流后, 內(nèi)孤立波的形成受到了干擾, 能量被背景流場帶走分散掉, 以致最終產(chǎn)生的波動波長較長, 振幅較小, 能量非常不集中, 推測這一影響是導致SAR觀測內(nèi)波出現(xiàn)頻率呈現(xiàn)季節(jié)變化的關鍵原因。而對比實驗替換夏季溫鹽場為冬季溫鹽場, 對模擬結果并無實質(zhì)上的影響, 說明冬夏季溫鹽場的差異并非是SAR觀測內(nèi)波出現(xiàn)頻率呈季節(jié)變化的重要因素。

    圖10 冬季內(nèi)孤立波形成過程對比圖Fig. 10 Generation of IWs in winter

    由于二維近似, 忽略了模擬區(qū)域的三維地形及波波相互作用等因素, 模擬的準確性還有待提高。另外黑潮入侵南海的具體形式尚無明確的結論, 如入侵具體時間、位置、流速大小等。因此, 三維、實際地形、明確的黑潮入侵形式也成為今后進一步研究的目標。

    [1] Apel J R. A new analytical model for internal solitons in the ocean[J]. J Phys Oceanogr, 2003, 33: 2247-2269.

    [2] Helfrich K R, Melville W K. Long nonlinear internal waves[J]. Annu Rev Fluid Mech, 2006, 38: 395-425.

    [3] Zheng Q, Susanto R D, Ho C R , et al. Statistical and dynamical analyses of generation mechanisms of solitary internal waves in the northern South China Sea[J]. J Geophys Res, 2007, 112: C03021, doi: 10.1029/ 2006JC003551.

    [4] Cai Shuqun, Long X, Z Gan Z. A numerical study of the generation and propagation of internal solitary waves in the Luzon Strait[J]. Oceanologica Acta, 2002, 25: 51-60.

    [5] Yuan Y, Zheng Q, Dai D , et al. Mechanism of internal waves in the Luzon Strait[J]. J Geophys Res, 2006, 111: C11S17, doi: 10.1029/2005JC003198.

    [6] Du T, Tseng Y H, Yan X H. Impacts of tidal currents and Kuroshio intrusion on the generation of nonlinear internal waves in Luzon Strait[J]. J Geophys Res, 2008, 113: C08015.

    [7] Maxworthy T. A note on the internal solitary waves produced by tidal flow over a three-dimensional ridge[J]. J Geophys Res, 1979, 84: 338-346.

    [8] Wang G, Qiao F, Dai D. A 2-D numerical modeling of the generation and propagation of internal solitary waves in the Luzon Strait[J]. Acta Oceanologica Sinica, 2010, 29(6): 1-11.

    [9] 申輝. 海洋內(nèi)波的遙感與數(shù)值模擬研究 [D].北京:中國科學院研究生, 2005.

    Shen Hui. Remote sensing and numerical simulation for the internal waves in the ocean[D]. Beijing: University of Chinese Academy of Sciences, 2005.

    [10] 劉秦玉, 劉倬騰, 鄭世培, 等. 黑潮在呂宋海峽的形變及動力機制[J]. 青島海洋大學學報, 1996, 26(4): 413-420.

    Liu Qinyu, Liu Zhuoteng, Zheng Shipei, et al. Thedeformation mechanism of Kuroshio in Luzong Strait[J]. Journal of Ocean University of China, 1996, 26(4): 413-420.

    [11] Caruso M J, Gawarkiewicz G G, Beardsley R C. Interannual variability of the Kuroshio intrusion in the South China Sea[J]. Journal of Oceanography, 2006, 62(4), 559-575.

    [12] Marshall J, Adroft A, Hill C, et al. A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers[J]. Journal of Geophysical Research, 1997, 102(C3): 5753-5766.

    [13] Vlasenko V, Stashchuk N, Guo C , et al. Multimodal structure of baroclinic tides in the South China Sea[J]. Nonlin Processes Geophys, 2010, 17: 529-543.

    (本文編輯: 劉珊珊)

    Effect of the Kuroshio on the generation of internal waves in the Northern South China Sea

    GUO Da-quan1, 2, 3, SHEN Hui1, 3
    (1. Institute of Oceanology, the Chinese Academy of Sciences, Qingdao 266071, China; 2. University of the Chinese Academy of Sciences, Beijing 100049, China; 3. Key Laboratory of Ocean Circulation and Wave, the Chinese Academy of Sciences, Qingdao 266071, China)

    Mar., 25, 2013

    internal waves; synthetic aperture radar; 2 dimensional; MITgcm; the Kuroshio

    In this study, we analyzed the Kuroshio’s effect on the generation of ocean internal waves (IWs) in the northern South China Sea (NSCS). We statistically analyzed IW occurrences in satellite synthetic aperture radar (SAR) images from 2007 to 2010. The monthly SAR-observed IW occurrence frequencies show that the highest frequencies are distributed from May to August and reach a peak in June, with a maximum frequency of 29%. Lowest frequencies are distributed in winter from November through February, with a minimum frequency of 0.5% in both December and January. The intrusion of the Kuroshio into the SCS represents a seasonal variability that matches the IWs occurrence variability well. A 2-D non-hydrostatical model properly simulates the generation of IWs and, with the K-branch current added into the domain, the IWs appear to stretch and the energy of the depression appears to scatter.

    P731.24

    A

    1000-3096(2016)02-0136-09

    10.11759/hykx20130325003

    2013-03-25;

    2013-05-07

    國家自然科學基金重點項目(41030855)

    郭大全(1987-), 男, 山東青島人, 碩士研究生, 主要從事海洋遙感與數(shù)值模擬研究, 電話: 15963202590, E-mai: guodaquan129 @126.com; 申輝(1978-), 男, 研究員, 博士, 主要從事物理海洋與海洋遙感研究, 電話0532-82898783, E-mail: shenhui@qdio.ac.cn

    猜你喜歡
    黑潮內(nèi)波潮流
    孤立內(nèi)波對過渡海域聲場干涉結構的影響分析
    海洋通報(2022年6期)2023-01-07 02:49:06
    譜松弛動力降尺度方法及其在呂宋海峽黑潮模擬中的應用
    海洋通報(2022年6期)2023-01-07 02:48:56
    內(nèi)波與死水,連潛艇都怕的海浪
    電腦報(2021年20期)2021-08-10 02:28:41
    基于MODIS 遙感影像的安達曼海內(nèi)波特征參數(shù)分布及生成周期研究
    海洋學報(2020年9期)2020-10-09 08:50:08
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    潮流
    足球周刊(2016年14期)2016-11-02 11:47:59
    潮流
    足球周刊(2016年15期)2016-11-02 11:44:02
    潮流
    足球周刊(2016年10期)2016-10-08 18:50:29
    基于20 a衛(wèi)星高度計數(shù)據(jù)的黑潮變異特征
    海洋科學(2016年1期)2016-07-06 09:42:34
    東海冷渦對黑潮熱輸運年代際躍變的響應
    av专区在线播放| 在线看三级毛片| 国产在视频线在精品| 日韩精品中文字幕看吧| 一a级毛片在线观看| 国产精品爽爽va在线观看网站| 1024手机看黄色片| 国产大屁股一区二区在线视频| 午夜视频国产福利| 精品人妻偷拍中文字幕| 国产亚洲91精品色在线| 99久久精品国产国产毛片| 亚洲人成网站在线播| 欧美人与善性xxx| 亚洲精品在线观看二区| 偷拍熟女少妇极品色| 国产麻豆成人av免费视频| 欧美一级a爱片免费观看看| 婷婷精品国产亚洲av| 美女大奶头视频| 欧美性猛交╳xxx乱大交人| 国产精品一区www在线观看 | 最好的美女福利视频网| 日韩欧美免费精品| 老司机福利观看| 亚洲无线观看免费| 精华霜和精华液先用哪个| 美女xxoo啪啪120秒动态图| 精品人妻熟女av久视频| av在线亚洲专区| 久久欧美精品欧美久久欧美| 成人特级av手机在线观看| 最近在线观看免费完整版| 国产 一区 欧美 日韩| 亚洲精品一卡2卡三卡4卡5卡| 少妇被粗大猛烈的视频| 午夜福利在线观看吧| 99热这里只有精品一区| 国产精品亚洲一级av第二区| 韩国av一区二区三区四区| 不卡一级毛片| 天堂√8在线中文| 禁无遮挡网站| 大又大粗又爽又黄少妇毛片口| 亚洲第一区二区三区不卡| 香蕉av资源在线| 一边摸一边抽搐一进一小说| 乱码一卡2卡4卡精品| 在线a可以看的网站| 国产高清激情床上av| 欧美激情在线99| 色av中文字幕| 在线a可以看的网站| 国产精品久久久久久精品电影| 免费搜索国产男女视频| 乱码一卡2卡4卡精品| 亚洲人成网站在线播放欧美日韩| 欧美黑人欧美精品刺激| 可以在线观看的亚洲视频| 窝窝影院91人妻| 成人av一区二区三区在线看| 亚洲av日韩精品久久久久久密| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜添小说| 日本免费一区二区三区高清不卡| 听说在线观看完整版免费高清| 99热6这里只有精品| 麻豆一二三区av精品| 久久热精品热| 国产成人福利小说| 精品无人区乱码1区二区| 亚洲中文日韩欧美视频| 久久精品人妻少妇| 国产精品女同一区二区软件 | 精品国产三级普通话版| 亚洲精品久久国产高清桃花| 大又大粗又爽又黄少妇毛片口| 悠悠久久av| 麻豆国产97在线/欧美| 人妻少妇偷人精品九色| 伦精品一区二区三区| 人人妻人人澡欧美一区二区| 夜夜爽天天搞| 国产大屁股一区二区在线视频| 久久精品91蜜桃| 国产精品人妻久久久影院| 1024手机看黄色片| 成人永久免费在线观看视频| 精品久久久久久久久亚洲 | 不卡一级毛片| 99九九线精品视频在线观看视频| 日本黄色片子视频| 黄色日韩在线| 国产精品一及| 亚洲欧美日韩高清专用| 在线播放无遮挡| 午夜免费成人在线视频| 色播亚洲综合网| 日本三级黄在线观看| 麻豆成人av在线观看| 国产精品一及| av女优亚洲男人天堂| 国产精品野战在线观看| 亚洲乱码一区二区免费版| 亚洲av一区综合| 亚洲av不卡在线观看| 亚洲av五月六月丁香网| 91麻豆精品激情在线观看国产| 亚洲欧美清纯卡通| 午夜日韩欧美国产| 国内揄拍国产精品人妻在线| 欧美成人a在线观看| 午夜久久久久精精品| 成年女人看的毛片在线观看| 久久久精品欧美日韩精品| 久久精品国产99精品国产亚洲性色| 看免费成人av毛片| 国产成年人精品一区二区| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 九九久久精品国产亚洲av麻豆| 久久久久国内视频| 国产精品人妻久久久久久| 亚洲美女黄片视频| 一边摸一边抽搐一进一小说| 在线国产一区二区在线| 欧美日韩精品成人综合77777| 午夜福利在线观看吧| 亚洲av成人精品一区久久| 成人av一区二区三区在线看| 岛国在线免费视频观看| 国产又黄又爽又无遮挡在线| 日韩精品青青久久久久久| 国内毛片毛片毛片毛片毛片| 中文字幕免费在线视频6| 在线观看66精品国产| 国产高清视频在线观看网站| 免费看av在线观看网站| 国产精品99久久久久久久久| 日韩欧美一区二区三区在线观看| 毛片女人毛片| 日韩在线高清观看一区二区三区 | 99热6这里只有精品| 免费看光身美女| 日日撸夜夜添| 国产美女午夜福利| 精品福利观看| 变态另类丝袜制服| 色av中文字幕| 大型黄色视频在线免费观看| 久久午夜福利片| 免费观看精品视频网站| 18禁黄网站禁片免费观看直播| 中亚洲国语对白在线视频| 欧美一级a爱片免费观看看| 日韩精品青青久久久久久| 两个人视频免费观看高清| 午夜久久久久精精品| 精品人妻偷拍中文字幕| 亚洲av熟女| 国产aⅴ精品一区二区三区波| 我的老师免费观看完整版| 听说在线观看完整版免费高清| 欧美高清性xxxxhd video| 亚洲精品一卡2卡三卡4卡5卡| 十八禁网站免费在线| 久久精品国产亚洲网站| 九九在线视频观看精品| 欧美日韩乱码在线| 免费观看精品视频网站| 少妇人妻一区二区三区视频| 好男人在线观看高清免费视频| 伦精品一区二区三区| 色综合亚洲欧美另类图片| 麻豆国产av国片精品| 成人永久免费在线观看视频| 国产亚洲91精品色在线| 91麻豆av在线| 日韩人妻高清精品专区| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 国产国拍精品亚洲av在线观看| 国产日本99.免费观看| 国产熟女欧美一区二区| 欧美性感艳星| 高清在线国产一区| 午夜影院日韩av| 国产成人影院久久av| 小蜜桃在线观看免费完整版高清| 久久久成人免费电影| 国产真实乱freesex| 国产探花极品一区二区| 亚洲内射少妇av| 亚洲欧美精品综合久久99| 在线a可以看的网站| a级毛片免费高清观看在线播放| 老司机福利观看| 悠悠久久av| 69人妻影院| 国产成人a区在线观看| 国产三级在线视频| 欧美色视频一区免费| 国产精品乱码一区二三区的特点| 97超视频在线观看视频| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久久电影| 亚洲av熟女| 窝窝影院91人妻| 18禁黄网站禁片免费观看直播| 久久午夜亚洲精品久久| 特级一级黄色大片| АⅤ资源中文在线天堂| 色尼玛亚洲综合影院| 91午夜精品亚洲一区二区三区 | 国产极品精品免费视频能看的| 欧美性猛交╳xxx乱大交人| 国产综合懂色| 3wmmmm亚洲av在线观看| x7x7x7水蜜桃| 国产极品精品免费视频能看的| 观看美女的网站| 日本一本二区三区精品| 欧美bdsm另类| 美女高潮的动态| 亚洲第一区二区三区不卡| 一夜夜www| av在线老鸭窝| 最新中文字幕久久久久| 韩国av在线不卡| 欧美最新免费一区二区三区| 搡老妇女老女人老熟妇| 18+在线观看网站| 一进一出好大好爽视频| 人妻久久中文字幕网| 国产女主播在线喷水免费视频网站 | 1000部很黄的大片| 色精品久久人妻99蜜桃| 成人永久免费在线观看视频| 亚洲无线在线观看| 国产高清三级在线| www.色视频.com| 在线看三级毛片| 亚洲第一区二区三区不卡| 欧美成人性av电影在线观看| 亚洲精华国产精华液的使用体验 | 自拍偷自拍亚洲精品老妇| 国模一区二区三区四区视频| 欧美精品啪啪一区二区三区| a级毛片免费高清观看在线播放| 久久久久精品国产欧美久久久| 欧美又色又爽又黄视频| 国产真实乱freesex| .国产精品久久| 最近视频中文字幕2019在线8| 亚洲天堂国产精品一区在线| 一区福利在线观看| 国产v大片淫在线免费观看| 国产精品综合久久久久久久免费| 男女视频在线观看网站免费| 国内精品一区二区在线观看| 九九热线精品视视频播放| 亚洲人成网站在线播| 国产三级中文精品| 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 乱人视频在线观看| 中文字幕高清在线视频| 久久精品国产清高在天天线| 欧美日韩亚洲国产一区二区在线观看| 人人妻人人澡欧美一区二区| 久久午夜亚洲精品久久| 麻豆成人av在线观看| 一个人免费在线观看电影| 免费人成在线观看视频色| 黄色女人牲交| 成熟少妇高潮喷水视频| 2021天堂中文幕一二区在线观| 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 亚洲av熟女| 欧美日本视频| 国产欧美日韩精品亚洲av| 精品久久久久久久久久久久久| 亚洲欧美日韩无卡精品| 成年免费大片在线观看| 日本免费一区二区三区高清不卡| 久9热在线精品视频| 国产老妇女一区| 国产高清有码在线观看视频| 欧美人与善性xxx| 色综合站精品国产| 国产av麻豆久久久久久久| 国产私拍福利视频在线观看| 久久精品国产鲁丝片午夜精品 | 老女人水多毛片| 婷婷色综合大香蕉| 一级av片app| 嫩草影视91久久| 精品欧美国产一区二区三| 网址你懂的国产日韩在线| 国产伦人伦偷精品视频| 伦精品一区二区三区| 成人特级av手机在线观看| 日韩大尺度精品在线看网址| 一进一出好大好爽视频| 99久久成人亚洲精品观看| 97人妻精品一区二区三区麻豆| 美女 人体艺术 gogo| 亚洲中文字幕一区二区三区有码在线看| 久久亚洲真实| 色综合色国产| 91久久精品国产一区二区成人| 亚洲国产日韩欧美精品在线观看| 国产成人福利小说| 国产欧美日韩精品一区二区| 国产黄a三级三级三级人| 日本一本二区三区精品| 亚洲精品色激情综合| 久久精品国产亚洲av天美| 亚洲精品日韩av片在线观看| 老司机午夜福利在线观看视频| av中文乱码字幕在线| 精品人妻熟女av久视频| 日韩人妻高清精品专区| 少妇人妻精品综合一区二区 | 真人做人爱边吃奶动态| 亚洲精品粉嫩美女一区| 免费人成在线观看视频色| 亚洲av免费在线观看| 直男gayav资源| 少妇丰满av| 成人二区视频| 一级a爱片免费观看的视频| 在线看三级毛片| 亚洲av中文av极速乱 | 亚洲成人精品中文字幕电影| 夜夜夜夜夜久久久久| 亚洲精品亚洲一区二区| 久99久视频精品免费| 亚洲国产高清在线一区二区三| 精品久久久久久久久av| 两人在一起打扑克的视频| 网址你懂的国产日韩在线| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 免费搜索国产男女视频| av福利片在线观看| 夜夜爽天天搞| 伦理电影大哥的女人| 亚洲七黄色美女视频| 极品教师在线视频| 2021天堂中文幕一二区在线观| 亚洲人成伊人成综合网2020| 日本与韩国留学比较| 国内精品久久久久久久电影| 婷婷精品国产亚洲av在线| 51国产日韩欧美| ponron亚洲| 欧美色视频一区免费| 国内精品美女久久久久久| 桃红色精品国产亚洲av| 身体一侧抽搐| 国产淫片久久久久久久久| 国产一区二区三区视频了| 久久国内精品自在自线图片| 在线天堂最新版资源| 日韩欧美精品免费久久| 久久精品91蜜桃| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 免费观看的影片在线观看| 亚洲成人免费电影在线观看| 久久午夜亚洲精品久久| 久久久久久九九精品二区国产| 午夜爱爱视频在线播放| 国产精品久久久久久亚洲av鲁大| 色视频www国产| 亚洲成人精品中文字幕电影| 成人亚洲精品av一区二区| 国产不卡一卡二| 中文字幕免费在线视频6| 男女下面进入的视频免费午夜| 人妻久久中文字幕网| 黄片wwwwww| 成人鲁丝片一二三区免费| 深夜a级毛片| 又爽又黄无遮挡网站| 91av网一区二区| 99热这里只有是精品在线观看| 最好的美女福利视频网| 精品午夜福利视频在线观看一区| 久久精品夜夜夜夜夜久久蜜豆| 国产久久久一区二区三区| 亚洲综合色惰| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 99热精品在线国产| 少妇的逼好多水| 国产伦精品一区二区三区四那| 久久精品综合一区二区三区| 久久久久免费精品人妻一区二区| 国产黄色小视频在线观看| 国产精品乱码一区二三区的特点| 久久天躁狠狠躁夜夜2o2o| 国产老妇女一区| 99久久精品国产国产毛片| 麻豆久久精品国产亚洲av| 可以在线观看的亚洲视频| 欧美又色又爽又黄视频| 91精品国产九色| 日韩精品有码人妻一区| 国产在视频线在精品| 亚州av有码| 免费观看精品视频网站| 亚洲狠狠婷婷综合久久图片| 国产精品乱码一区二三区的特点| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 可以在线观看的亚洲视频| av天堂在线播放| 成人鲁丝片一二三区免费| eeuss影院久久| 欧美国产日韩亚洲一区| 日本精品一区二区三区蜜桃| 欧美潮喷喷水| 偷拍熟女少妇极品色| 久久久久国产精品人妻aⅴ院| 日本 欧美在线| 一区二区三区激情视频| 亚洲成人久久爱视频| 精品人妻1区二区| 男人狂女人下面高潮的视频| 国产免费av片在线观看野外av| 国内久久婷婷六月综合欲色啪| 国产精品,欧美在线| 欧美最黄视频在线播放免费| 欧美性猛交╳xxx乱大交人| 免费看光身美女| 日本 欧美在线| 免费观看在线日韩| 一级黄片播放器| 18禁黄网站禁片免费观看直播| 亚洲图色成人| 免费高清视频大片| 亚洲午夜理论影院| 99国产精品一区二区蜜桃av| 亚洲第一区二区三区不卡| 日韩 亚洲 欧美在线| 午夜福利18| 别揉我奶头 嗯啊视频| 久久久久久久久久黄片| 精品久久久久久久人妻蜜臀av| 一级毛片久久久久久久久女| 亚洲专区国产一区二区| 久9热在线精品视频| 精品人妻1区二区| 欧美极品一区二区三区四区| 亚洲av免费在线观看| 搡女人真爽免费视频火全软件 | 欧美一级a爱片免费观看看| 久久久久久大精品| 欧美日韩瑟瑟在线播放| 亚洲欧美清纯卡通| 日日啪夜夜撸| 亚洲久久久久久中文字幕| 亚洲精品色激情综合| 成人精品一区二区免费| 欧美激情在线99| 久久精品91蜜桃| 在线免费观看不下载黄p国产 | 九色国产91popny在线| 哪里可以看免费的av片| 亚洲性久久影院| 久久天躁狠狠躁夜夜2o2o| 悠悠久久av| 欧美3d第一页| 一进一出好大好爽视频| 赤兔流量卡办理| 岛国在线免费视频观看| 国产一区二区三区视频了| 免费在线观看成人毛片| 久久精品国产亚洲av香蕉五月| 免费观看精品视频网站| 天美传媒精品一区二区| av.在线天堂| 国产中年淑女户外野战色| 91av网一区二区| 最近最新中文字幕大全电影3| 国产黄a三级三级三级人| 国产免费男女视频| 男插女下体视频免费在线播放| 欧美+亚洲+日韩+国产| 色av中文字幕| 国产国拍精品亚洲av在线观看| 国产伦一二天堂av在线观看| 性色avwww在线观看| 91午夜精品亚洲一区二区三区 | 亚洲av成人av| 成人永久免费在线观看视频| 亚洲欧美清纯卡通| 人人妻人人澡欧美一区二区| 午夜久久久久精精品| 综合色av麻豆| 黄色日韩在线| 国产成人av教育| 一a级毛片在线观看| 中文字幕久久专区| 搡老妇女老女人老熟妇| 亚洲中文字幕一区二区三区有码在线看| 麻豆精品久久久久久蜜桃| 久久婷婷人人爽人人干人人爱| 99久久精品一区二区三区| 最近中文字幕高清免费大全6 | 中亚洲国语对白在线视频| 日本免费一区二区三区高清不卡| 成熟少妇高潮喷水视频| 床上黄色一级片| 不卡视频在线观看欧美| 国产精品嫩草影院av在线观看 | 在线免费观看不下载黄p国产 | a级一级毛片免费在线观看| 女人十人毛片免费观看3o分钟| 亚洲一区二区三区色噜噜| 国产69精品久久久久777片| 一进一出好大好爽视频| 欧美日本视频| 桃红色精品国产亚洲av| 亚洲成人精品中文字幕电影| 精品久久久久久久久久免费视频| 午夜亚洲福利在线播放| 三级国产精品欧美在线观看| 婷婷丁香在线五月| 九色成人免费人妻av| 亚洲国产精品合色在线| 成人av在线播放网站| 少妇丰满av| 日韩一区二区视频免费看| x7x7x7水蜜桃| 日日撸夜夜添| 91狼人影院| 免费观看的影片在线观看| 无遮挡黄片免费观看| 99热这里只有是精品50| 69人妻影院| 悠悠久久av| 精品99又大又爽又粗少妇毛片 | 久久久久国产精品人妻aⅴ院| 麻豆国产av国片精品| 色噜噜av男人的天堂激情| 日韩大尺度精品在线看网址| 亚洲av一区综合| 精品福利观看| 99久久精品国产国产毛片| 久久久久久九九精品二区国产| 午夜激情欧美在线| 日本色播在线视频| 别揉我奶头 嗯啊视频| 美女高潮喷水抽搐中文字幕| 99热精品在线国产| 中文字幕熟女人妻在线| 亚洲最大成人中文| 午夜福利在线观看吧| 身体一侧抽搐| 亚洲人成伊人成综合网2020| 搡老岳熟女国产| 99久久成人亚洲精品观看| 波多野结衣高清作品| 午夜日韩欧美国产| 欧美成人性av电影在线观看| 在线观看免费视频日本深夜| 午夜福利视频1000在线观看| 中文字幕av成人在线电影| 禁无遮挡网站| 久久国产精品人妻蜜桃| 亚洲国产欧美人成| 不卡视频在线观看欧美| 高清毛片免费观看视频网站| 欧美一区二区精品小视频在线| 美女高潮的动态| www.色视频.com| 赤兔流量卡办理| 欧美zozozo另类| 欧美区成人在线视频| 国产麻豆成人av免费视频| 3wmmmm亚洲av在线观看| 亚洲av成人精品一区久久| 中文字幕熟女人妻在线| 国产黄片美女视频| 亚洲欧美日韩高清专用| 亚洲精品亚洲一区二区| 九九爱精品视频在线观看| 国产亚洲91精品色在线| 搞女人的毛片| 禁无遮挡网站| 午夜福利高清视频| 成人二区视频| 日本精品一区二区三区蜜桃| 国产成人一区二区在线| a在线观看视频网站| 此物有八面人人有两片| 女的被弄到高潮叫床怎么办 | 国产91精品成人一区二区三区| 99久久精品一区二区三区| 亚洲七黄色美女视频| 日本-黄色视频高清免费观看| 51国产日韩欧美| 一区福利在线观看| 99热精品在线国产| 两个人视频免费观看高清| 国产亚洲精品久久久久久毛片| 亚洲人与动物交配视频| 97人妻精品一区二区三区麻豆| 亚洲第一电影网av| АⅤ资源中文在线天堂| 国产伦一二天堂av在线观看|