• <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
    東海冷渦對黑潮熱輸運年代際躍變的響應
    午夜免费男女啪啪视频观看| 少妇猛男粗大的猛烈进出视频 | 波多野结衣巨乳人妻| 国产又色又爽无遮挡免| 18+在线观看网站| 欧美人与善性xxx| 免费观看av网站的网址| 2021天堂中文幕一二区在线观| 97超视频在线观看视频| 久久久久久久久大av| 成人国产av品久久久| 两个人的视频大全免费| 亚洲欧美成人综合另类久久久| 激情五月婷婷亚洲| 国产在视频线精品| 五月天丁香电影| 久久精品久久久久久久性| 中文字幕人妻熟人妻熟丝袜美| av在线app专区| 亚洲精品一区蜜桃| 熟女av电影| 一个人看的www免费观看视频| 亚洲精品aⅴ在线观看| 欧美亚洲 丝袜 人妻 在线| av又黄又爽大尺度在线免费看| 久久97久久精品| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久精品古装| 国产日韩欧美亚洲二区| videos熟女内射| 熟女av电影| 人体艺术视频欧美日本| 久久久久国产网址| 久久女婷五月综合色啪小说 | 国内揄拍国产精品人妻在线| 精品午夜福利在线看| av在线天堂中文字幕| 一区二区三区免费毛片| 禁无遮挡网站| 韩国av在线不卡| 人妻系列 视频| 精品一区二区三区视频在线| 一级黄片播放器| 午夜精品国产一区二区电影 | 舔av片在线| 秋霞在线观看毛片| 国产伦理片在线播放av一区| 亚洲精品久久久久久婷婷小说| 夫妻午夜视频| 97精品久久久久久久久久精品| 久久这里有精品视频免费| 边亲边吃奶的免费视频| 搡老乐熟女国产| 99热全是精品| 人妻一区二区av| 一区二区三区精品91| 成人漫画全彩无遮挡| 久久99热6这里只有精品| 18+在线观看网站| 神马国产精品三级电影在线观看| 激情五月婷婷亚洲| 亚洲精品中文字幕在线视频 | 国产精品蜜桃在线观看| 亚洲成色77777| 一区二区av电影网| 在线观看免费高清a一片| 内地一区二区视频在线| 大陆偷拍与自拍| 草草在线视频免费看| 91精品伊人久久大香线蕉| 国产爽快片一区二区三区| 少妇高潮的动态图| 一本色道久久久久久精品综合| 亚洲综合精品二区| 观看免费一级毛片| 久久精品人妻少妇| 国产精品精品国产色婷婷| 天堂俺去俺来也www色官网| 91精品一卡2卡3卡4卡| 亚洲精品,欧美精品| 校园人妻丝袜中文字幕| 亚洲精品乱码久久久v下载方式| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 永久免费av网站大全| 麻豆国产97在线/欧美| 97精品久久久久久久久久精品| 97超碰精品成人国产| 一区二区三区精品91| 嫩草影院新地址| 欧美日韩亚洲高清精品| 国产精品久久久久久久电影| 国产亚洲91精品色在线| 免费在线观看成人毛片| 国产白丝娇喘喷水9色精品| 亚洲国产av新网站| 黄色视频在线播放观看不卡| 看黄色毛片网站| 一区二区三区精品91| 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 国产精品.久久久| 日本wwww免费看| 性色avwww在线观看| 免费黄频网站在线观看国产| 三级国产精品欧美在线观看| 久久久久久国产a免费观看| 一级黄片播放器| 国产精品99久久久久久久久| 99热国产这里只有精品6| 一个人观看的视频www高清免费观看| 久久精品国产鲁丝片午夜精品| 精品熟女少妇av免费看| 亚洲国产欧美在线一区| 国产精品熟女久久久久浪| 亚洲图色成人| 欧美人与善性xxx| 人妻夜夜爽99麻豆av| 一本一本综合久久| 美女视频免费永久观看网站| 成人欧美大片| 国产熟女欧美一区二区| av国产免费在线观看| 国产一区二区亚洲精品在线观看| 国模一区二区三区四区视频| 热re99久久精品国产66热6| 99热这里只有是精品50| 国产精品一区二区三区四区免费观看| freevideosex欧美| 亚洲精品色激情综合| 国产成人免费无遮挡视频| 成年免费大片在线观看| 日日摸夜夜添夜夜添av毛片| 婷婷色综合大香蕉| 最近手机中文字幕大全| 人人妻人人爽人人添夜夜欢视频 | 色综合色国产| 国国产精品蜜臀av免费| 欧美 日韩 精品 国产| 夫妻性生交免费视频一级片| 啦啦啦在线观看免费高清www| 水蜜桃什么品种好| 色综合色国产| 成人特级av手机在线观看| 嘟嘟电影网在线观看| 人体艺术视频欧美日本| 舔av片在线| 欧美高清性xxxxhd video| 国产精品久久久久久久电影| 丰满乱子伦码专区| 免费黄网站久久成人精品| 全区人妻精品视频| 欧美性感艳星| 欧美一级a爱片免费观看看| 汤姆久久久久久久影院中文字幕| 亚洲,一卡二卡三卡| 国产亚洲av片在线观看秒播厂| 成人毛片a级毛片在线播放| 免费大片18禁| 黄色欧美视频在线观看| 只有这里有精品99| 免费观看的影片在线观看| 精品久久久久久久人妻蜜臀av| av免费观看日本| 精品久久久精品久久久| 久久亚洲国产成人精品v| 26uuu在线亚洲综合色| 一区二区三区乱码不卡18| 美女被艹到高潮喷水动态| 亚洲熟女精品中文字幕| 国产精品久久久久久精品电影小说 | 日日摸夜夜添夜夜添av毛片| 亚洲欧洲国产日韩| 日韩一区二区三区影片| 国产精品.久久久| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看| 中文字幕av成人在线电影| 久久精品久久精品一区二区三区| 涩涩av久久男人的天堂| 97人妻精品一区二区三区麻豆| 亚洲欧美清纯卡通| 国产欧美亚洲国产| 久久久精品94久久精品| 别揉我奶头 嗯啊视频| 纵有疾风起免费观看全集完整版| 久久精品久久精品一区二区三区| 国产色婷婷99| 日本色播在线视频| 大香蕉97超碰在线| 亚洲色图av天堂| 精品久久国产蜜桃| 在线观看人妻少妇| 日本色播在线视频| 97精品久久久久久久久久精品| 成年版毛片免费区| 亚洲av电影在线观看一区二区三区 | 汤姆久久久久久久影院中文字幕| 久久午夜福利片| 80岁老熟妇乱子伦牲交| 97超视频在线观看视频| 亚洲人成网站在线播| 日本黄大片高清| 少妇裸体淫交视频免费看高清| 新久久久久国产一级毛片| 国产精品国产三级专区第一集| 欧美成人a在线观看| 波多野结衣巨乳人妻| 国产精品av视频在线免费观看| 国产精品久久久久久精品电影小说 | 全区人妻精品视频| 欧美zozozo另类| 久久精品熟女亚洲av麻豆精品| 亚洲天堂av无毛| 在现免费观看毛片| 亚洲色图av天堂| 成人特级av手机在线观看| 别揉我奶头 嗯啊视频| 国产黄色免费在线视频| 日本欧美国产在线视频| 国产黄片美女视频| 丰满乱子伦码专区| 欧美97在线视频| 精华霜和精华液先用哪个| 日日摸夜夜添夜夜爱| 国产探花在线观看一区二区| 欧美bdsm另类| 国产视频内射| 国产午夜福利久久久久久| 亚洲精品久久午夜乱码| 97超碰精品成人国产| 日韩大片免费观看网站| 午夜精品国产一区二区电影 | 大香蕉97超碰在线| 老司机影院毛片| 一级二级三级毛片免费看| 男女那种视频在线观看| 26uuu在线亚洲综合色| 国产黄色视频一区二区在线观看| 成人高潮视频无遮挡免费网站| 国产黄a三级三级三级人| 精品99又大又爽又粗少妇毛片| 久久热精品热| 老司机影院成人| 99久久九九国产精品国产免费| 激情五月婷婷亚洲| 精品视频人人做人人爽| 一级毛片 在线播放| 国产精品久久久久久精品古装| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 久久久欧美国产精品| 最近最新中文字幕免费大全7| 韩国av在线不卡| 精品熟女少妇av免费看| 丰满少妇做爰视频| 在线精品无人区一区二区三 | 久久精品夜色国产| 我的老师免费观看完整版| av卡一久久| 亚洲国产精品999| 一边亲一边摸免费视频| 国产一区有黄有色的免费视频| 精品久久久久久久久亚洲| 免费黄色在线免费观看| 国产av码专区亚洲av| 69人妻影院| a级一级毛片免费在线观看| 日日撸夜夜添| 99热这里只有精品一区| 亚洲欧美日韩无卡精品| 久久99热6这里只有精品| 两个人的视频大全免费| av国产久精品久网站免费入址| 91狼人影院| 午夜激情久久久久久久| 久久久久久久久久久丰满| 深夜a级毛片| 久久久久国产精品人妻一区二区| 亚洲欧美精品自产自拍| 最后的刺客免费高清国语| 精品少妇黑人巨大在线播放| 97在线视频观看| 亚洲人成网站高清观看| 日韩大片免费观看网站| videos熟女内射| 少妇丰满av| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线 | 性色av一级| 日韩欧美 国产精品| 一区二区三区免费毛片| 久久女婷五月综合色啪小说 | 欧美高清成人免费视频www| 一级毛片 在线播放| 国产爽快片一区二区三区| 国产精品成人在线| 韩国av在线不卡| 国产女主播在线喷水免费视频网站| 王馨瑶露胸无遮挡在线观看| 久久鲁丝午夜福利片| av一本久久久久| 亚洲在线观看片| 久久精品久久久久久噜噜老黄| 免费黄色在线免费观看| 在线看a的网站| 欧美成人a在线观看| 日韩,欧美,国产一区二区三区| 国产高清国产精品国产三级 | 亚洲国产高清在线一区二区三| 久久6这里有精品| 高清视频免费观看一区二区| 亚洲内射少妇av| 免费黄频网站在线观看国产| 亚洲综合色惰| 国产综合懂色| 亚洲av.av天堂| 狠狠精品人妻久久久久久综合| 91aial.com中文字幕在线观看| 狂野欧美激情性bbbbbb| 午夜福利网站1000一区二区三区| 日本欧美国产在线视频| 麻豆久久精品国产亚洲av| 亚洲精品国产av成人精品| 黄色日韩在线| 欧美成人a在线观看| 九草在线视频观看| 女人十人毛片免费观看3o分钟| 人人妻人人爽人人添夜夜欢视频 | 亚洲欧美精品专区久久| 韩国高清视频一区二区三区| 亚洲av国产av综合av卡| 亚洲精品第二区| 亚洲精品aⅴ在线观看| 毛片女人毛片| 午夜福利在线观看免费完整高清在| 久久6这里有精品| 国产成人精品一,二区| 久久国产乱子免费精品| 男人狂女人下面高潮的视频| 国产熟女欧美一区二区| 嫩草影院新地址| 欧美激情国产日韩精品一区| 午夜亚洲福利在线播放| 黄色配什么色好看| 少妇猛男粗大的猛烈进出视频 | 在现免费观看毛片| 午夜福利网站1000一区二区三区| 麻豆成人午夜福利视频| 日韩一本色道免费dvd| 国产淫片久久久久久久久| av.在线天堂| 午夜老司机福利剧场| 国产高清不卡午夜福利| 国产亚洲av嫩草精品影院| 99精国产麻豆久久婷婷| 大香蕉97超碰在线| 欧美精品一区二区大全| 午夜激情福利司机影院| av播播在线观看一区| 日本黄色片子视频| 成人黄色视频免费在线看| 成人鲁丝片一二三区免费| 免费高清在线观看视频在线观看| 麻豆乱淫一区二区| 亚洲欧美精品专区久久| 女的被弄到高潮叫床怎么办| 婷婷色综合大香蕉| 日本三级黄在线观看| 日韩制服骚丝袜av| 在线观看人妻少妇| 午夜激情久久久久久久| 另类亚洲欧美激情| 欧美性感艳星| 看免费成人av毛片| 大片免费播放器 马上看| 久久女婷五月综合色啪小说 | 亚洲国产日韩一区二区| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看| 久久午夜福利片| 欧美区成人在线视频| 一区二区三区乱码不卡18| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 亚洲第一区二区三区不卡| 毛片女人毛片| 精品人妻偷拍中文字幕| 国产精品一二三区在线看| 欧美亚洲 丝袜 人妻 在线| 天堂俺去俺来也www色官网| 黄片无遮挡物在线观看| 国产精品蜜桃在线观看| 亚洲伊人久久精品综合| 欧美一级a爱片免费观看看| 2018国产大陆天天弄谢| 免费看不卡的av| 99九九线精品视频在线观看视频| 国产成人freesex在线| 观看美女的网站| 伊人久久精品亚洲午夜| 大片免费播放器 马上看| 精品国产露脸久久av麻豆| 久久97久久精品| 亚洲欧美日韩东京热| 丝袜美腿在线中文| 全区人妻精品视频| 亚洲欧美中文字幕日韩二区| 汤姆久久久久久久影院中文字幕| 国产精品伦人一区二区| 少妇人妻 视频| 又爽又黄a免费视频| 18禁裸乳无遮挡动漫免费视频 | 午夜免费鲁丝| 一个人看的www免费观看视频| 一级毛片我不卡| 国产精品久久久久久久电影| 丝袜脚勾引网站| 国产极品天堂在线| 蜜桃亚洲精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 免费黄频网站在线观看国产| 国产亚洲午夜精品一区二区久久 | 欧美成人a在线观看| 熟女电影av网| 最近2019中文字幕mv第一页| 亚洲成人精品中文字幕电影| 久久精品国产亚洲av天美| 我的女老师完整版在线观看| 美女被艹到高潮喷水动态| 搡老乐熟女国产| 国产一区二区三区综合在线观看 | 新久久久久国产一级毛片| 菩萨蛮人人尽说江南好唐韦庄| av在线老鸭窝| 91久久精品国产一区二区三区| 中国美白少妇内射xxxbb| 天美传媒精品一区二区| 欧美日韩综合久久久久久| 大香蕉97超碰在线| 日日啪夜夜撸| 熟妇人妻不卡中文字幕| 18禁裸乳无遮挡动漫免费视频 | av在线播放精品| 三级经典国产精品| 国产老妇伦熟女老妇高清| 亚洲精品日韩在线中文字幕| 中文字幕制服av| 日日撸夜夜添| 80岁老熟妇乱子伦牲交| av国产精品久久久久影院| 大陆偷拍与自拍| 亚洲国产精品成人久久小说| 亚洲自拍偷在线| 高清毛片免费看| 亚洲精品色激情综合| 另类亚洲欧美激情| 免费观看无遮挡的男女| 国产成人精品久久久久久| 亚洲欧美日韩另类电影网站 | 精品午夜福利在线看| 精品久久久久久久人妻蜜臀av| 久久ye,这里只有精品| 亚洲最大成人中文| 最近手机中文字幕大全| 亚洲av中文av极速乱| 国产精品蜜桃在线观看| 各种免费的搞黄视频| 国产午夜精品一二区理论片| 亚洲aⅴ乱码一区二区在线播放| 全区人妻精品视频| 免费在线观看成人毛片| 一区二区av电影网| 亚洲丝袜综合中文字幕| 啦啦啦中文免费视频观看日本| av福利片在线观看| 国产美女午夜福利| 免费观看a级毛片全部| 亚洲国产精品国产精品| 久久久久久久午夜电影| 大香蕉久久网| 99久国产av精品国产电影| 亚洲一区二区三区欧美精品 | 国产成人精品婷婷| 亚洲国产精品999| 亚洲高清免费不卡视频| 日韩人妻高清精品专区| 欧美日韩亚洲高清精品| 777米奇影视久久| 亚洲国产精品999| 久久午夜福利片| 91久久精品电影网| 亚洲精品456在线播放app| 最近最新中文字幕免费大全7| 亚洲人成网站在线观看播放| 91久久精品国产一区二区成人| 亚洲欧美一区二区三区黑人 | 黄色视频在线播放观看不卡| 3wmmmm亚洲av在线观看| 夜夜爽夜夜爽视频| 精品少妇久久久久久888优播| 在线天堂最新版资源| 免费黄色在线免费观看| 草草在线视频免费看| 一级毛片aaaaaa免费看小| 亚洲av一区综合| 人体艺术视频欧美日本| 国产精品国产三级国产专区5o| 少妇人妻 视频| 日韩一区二区视频免费看| 国精品久久久久久国模美| 在线看a的网站| 男人狂女人下面高潮的视频| 五月玫瑰六月丁香| 搞女人的毛片| 老司机影院成人| 日本色播在线视频| 夜夜爽夜夜爽视频| 最近中文字幕高清免费大全6| 精品熟女少妇av免费看| 国产69精品久久久久777片| 久久久久久伊人网av| 欧美高清成人免费视频www| 亚洲婷婷狠狠爱综合网| 久久精品人妻少妇| 亚洲av在线观看美女高潮| 亚洲欧洲日产国产| 精品酒店卫生间| av女优亚洲男人天堂| 精华霜和精华液先用哪个| 日韩视频在线欧美| 久久6这里有精品| 纵有疾风起免费观看全集完整版| 成人亚洲精品av一区二区| 欧美精品国产亚洲| 人妻制服诱惑在线中文字幕| 在线免费十八禁| 少妇丰满av| 人妻 亚洲 视频| 国产免费视频播放在线视频| 日韩av不卡免费在线播放| 别揉我奶头 嗯啊视频| 精品一区在线观看国产| 欧美最新免费一区二区三区| 秋霞在线观看毛片| 最近最新中文字幕免费大全7| 国产伦在线观看视频一区| 国语对白做爰xxxⅹ性视频网站| 2021天堂中文幕一二区在线观| tube8黄色片| 亚洲精品久久午夜乱码| 两个人的视频大全免费| 亚洲人成网站在线播| 丰满乱子伦码专区| 99视频精品全部免费 在线| 性色av一级| xxx大片免费视频| 国产欧美亚洲国产| 亚洲欧美一区二区三区国产| 国产精品久久久久久精品电影| 中国美白少妇内射xxxbb| 国产亚洲精品久久久com| av免费观看日本| 国产欧美日韩精品一区二区| 国产精品久久久久久av不卡| 汤姆久久久久久久影院中文字幕| 最近最新中文字幕大全电影3| 少妇裸体淫交视频免费看高清| 国产精品一区二区性色av| 中文字幕亚洲精品专区| 丰满乱子伦码专区| 国产av码专区亚洲av| 亚洲精品中文字幕在线视频 | 国产黄色免费在线视频| 老女人水多毛片| 久久99热这里只频精品6学生| av在线蜜桃| 99热这里只有精品一区| 我要看日韩黄色一级片| av.在线天堂| 18禁动态无遮挡网站| 久久久精品94久久精品| 在线天堂最新版资源| 欧美成人精品欧美一级黄| 日韩av免费高清视频| 黄色日韩在线| 性色av一级| 高清av免费在线| av天堂中文字幕网| 超碰97精品在线观看| 一级a做视频免费观看| 色网站视频免费| 免费看日本二区| 欧美国产精品一级二级三级 | 白带黄色成豆腐渣| 亚洲欧美日韩另类电影网站 | 少妇裸体淫交视频免费看高清| 国产高潮美女av| 国产乱人偷精品视频| 偷拍熟女少妇极品色| 国产女主播在线喷水免费视频网站| 直男gayav资源| 观看美女的网站| 男女啪啪激烈高潮av片| 国产高清三级在线| 久久ye,这里只有精品| 成年版毛片免费区| 免费观看无遮挡的男女| 蜜桃亚洲精品一区二区三区| 亚洲真实伦在线观看| 成人亚洲精品av一区二区| 黑人高潮一二区| 人人妻人人看人人澡| 丝袜美腿在线中文| 国产免费福利视频在线观看|