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

    海浪破碎對北太平洋海表面溫度模擬的影響

    2017-10-13 22:50:55劉子龍蔣國榮
    海洋科學(xué) 2017年3期
    關(guān)鍵詞:表面溫度浮標(biāo)環(huán)流

    劉子龍, 史 劍, 蔣國榮

    ?

    海浪破碎對北太平洋海表面溫度模擬的影響

    劉子龍, 史 劍, 蔣國榮

    (解放軍理工大學(xué)氣象海洋學(xué)院, 江蘇南京211101)

    基于海浪模式WAVEWATCH III模擬北太平洋海浪要素, 結(jié)合NDBC浮標(biāo)資料進(jìn)行驗(yàn)證, 發(fā)現(xiàn)模擬出的有效波高與浮標(biāo)測量值具有很好的一致性?;诟倪M(jìn)型白冠覆蓋率耗散模型, 利用海浪模式模擬出的有效波高、有效波周期和摩擦速度等海浪要素計(jì)算出單位面積水柱內(nèi)因海浪破碎產(chǎn)生的湍動(dòng)能通量。通過改變環(huán)流模式sbPOM湍動(dòng)能方程的上邊界條件, 引入海浪破碎產(chǎn)生的湍動(dòng)能通量, 并探究海浪破碎對北太平洋海表面溫度模擬的影響。研究表明, 由于海浪破碎的引入, 環(huán)流模式sbPOM對北太平洋海表面溫度模擬的準(zhǔn)確程度得到提升, 這為大氣模式提供一個(gè)準(zhǔn)確的北太平洋下邊界條件具有重要意義。

    海浪破碎; 海表面溫度; WAVEWATCH III; sbPOM

    近年來, 隨著對上層海洋和低層大氣物理過程耦合研究的不斷深入, 對海浪運(yùn)動(dòng)的研究也得到了前所未有的關(guān)注, 所關(guān)注的焦點(diǎn)并非海浪現(xiàn)象本身, 而是海浪對其他海洋動(dòng)力過程的影響[1]。海浪運(yùn)動(dòng)導(dǎo)致海洋上層混合的機(jī)制主要有: 波致雷諾應(yīng)力引發(fā)的湍流混合、波浪破碎誘導(dǎo)的湍流混合、Stocks漂與風(fēng)驅(qū)動(dòng)的剪切流相互作用形成的Langmuir環(huán)流混合等[2]。本文在前人研究的基礎(chǔ)上, 僅從海浪破碎誘導(dǎo)的湍流混合對海洋表面溫度模擬產(chǎn)生的影響方面展開研究。

    海浪是存在于海-氣界面上的一種重要的運(yùn)動(dòng)形式, 海浪破碎對海洋上混合層中的湍流生成及其混合有著直接的影響。大量實(shí)驗(yàn)觀測結(jié)果表明[3-5], 海浪破碎將促進(jìn)海洋上混合層中湍流生成, 加強(qiáng)湍流垂向混合, 其影響深度與海浪成長狀態(tài)有關(guān)。Toba等[6]根據(jù)觀測指出海浪破碎的影響深度可達(dá)5倍有效波高。Thomson等[27]在數(shù)值模式中引入波浪破碎導(dǎo)致的湍動(dòng)能通量, 發(fā)現(xiàn)波浪破碎的影響深度與垂直參考系相關(guān)。Alari等[20]基于耦合環(huán)流模式NEMO和海浪模式WAM引入海浪破碎作用, 發(fā)現(xiàn)海浪破碎對海表面溫度、溫度的垂直分布和上升流都產(chǎn)生較大影響。Breivik等[28]基于NEMO環(huán)流模式, 討論了海表應(yīng)力、波浪破碎導(dǎo)致的湍動(dòng)能通量和Stokes Coriolis應(yīng)力對混合層的影響, 發(fā)現(xiàn)3個(gè)物理過程對減小海表面溫度偏差起到了巨大作用。管長龍等[1]總結(jié)前人所得, 認(rèn)為波浪破碎為海表附近的湍流生成源, 在混合層上部形成一個(gè)影響磁層, 該層中的湍流耗散率分布不再遵從固壁定律所規(guī)定的結(jié)果, 波浪破碎的影響深度可達(dá)5~10倍有效波高。孫群等[8-10]為了研究海浪破碎對混合層中湍能量收支的影響, 采用垂向一維湍封閉混合模型, 通過改變湍動(dòng)能方程的上邊界條件, 引入海浪破碎產(chǎn)生的湍動(dòng)能通量, 發(fā)現(xiàn)海浪破碎在海表產(chǎn)生的湍動(dòng)能通量影響了海洋上混合層中的各項(xiàng)湍能量收支間的局部平衡關(guān)系。Sun等[11]利用M-Y一維湍封閉模型, 通過將波浪破碎作為海表湍流源引入, 發(fā)現(xiàn)考慮波浪破碎效應(yīng)后, 模式模擬的海表面溫度比未考慮波浪破碎效應(yīng)時(shí)降低, 同時(shí), 混合層深度增加, 模擬效果更顯著。

    為了探究海浪破碎效應(yīng)在較大范圍內(nèi)對海表面溫度的影響, 本文在采用先進(jìn)的并行版本sbPOM環(huán)流模式[12]的基礎(chǔ)上, 區(qū)別于孫群等[8]利用海浪模式WAM計(jì)算海浪破碎導(dǎo)致能量耗散率的方式, 利用Guan等[7]總結(jié)計(jì)算出的基于白冠破碎的能量耗散率表達(dá)式, 通過改變環(huán)流模式sbPOM湍動(dòng)能方程的上邊界條件, 引入海浪破碎產(chǎn)生的湍動(dòng)能通量, 探究海浪破碎對北太平洋海域不同季節(jié)的海表面溫度(Sea Surface Temperature, SST)的影響。

    1 研究方法

    1.1 WAVEWATCH III海浪模式簡介及設(shè)置

    WAVEWATCH III海浪模式(以下簡稱 WW3)是基于第三代海浪模式WAM發(fā)展起來的當(dāng)前國際上最為成熟的幾個(gè)海浪模式之一[13], 具有穩(wěn)定性好、計(jì)算精度高等特點(diǎn), 目前已成為美國海洋環(huán)境預(yù)報(bào)中心的業(yè)務(wù)化海浪預(yù)報(bào)模式。大量海浪數(shù)值模擬研究表明[14-16], WW3能夠模擬出精度較高的海浪要素。

    本文采用NCEP(National Centers for Environmental Prediction)再分析風(fēng)場數(shù)據(jù)驅(qū)動(dòng)WW3海浪模式(該數(shù)據(jù)空間分辨率為1.875°×1.9°, 時(shí)間間隔為6 h),模式地形數(shù)據(jù)由全球分辨率5′×5′的ETOPO5地形插值得到。模擬了2014年全年的北太平洋海域海浪過程。設(shè)置WW3模式考慮風(fēng)攝入波動(dòng)、非線性的波波相互作用、白帽耗散和底摩擦等物理過程, 相應(yīng)的檢驗(yàn)參數(shù)采用默認(rèn)值。WW3模擬的區(qū)域范圍為: 10°S~ 66°N, 100°~280°E, 空間分辨率為5°×5°。海浪譜頻率分布范圍為 0.041 8~0.41 Hz, 共分25個(gè)頻段, 譜空間上離散為24個(gè)波向, 方向角分辨率為15°。模式最大全局時(shí)間積分步長取為2 400 s,方向和-theta方向最大CFL時(shí)間積分步長設(shè)為1 200 s, 最小源函數(shù)項(xiàng)時(shí)間積分步長本文中設(shè)為300 s。

    1.2 sbPOM環(huán)流模式簡介及設(shè)置

    sbPOM(stony brook parallel ocean model)模式[12]是基于POM(princeton ocean model)模式[25]進(jìn)行改進(jìn)得到的一種可并行運(yùn)算的環(huán)流模式, 劉欣等[17]曾采用區(qū)域海-氣耦合模式WRF-sbPOM, 對理想熱帶氣旋和海洋暖渦間的海氣相互作用進(jìn)行模擬, 獲得較好的模擬效果。劉子龍等[26]曾在環(huán)流模式sbPOM中引入海浪攪拌混合作用, 并探究其對北太平洋海表面溫度模擬的影響, 發(fā)現(xiàn)海浪攪拌混合作用的引入能夠較大的改善原模式對北太平洋海表面溫度的模擬。sbPOM除將POM改為并行外, 物理模型較原來無較大變化, 垂向坐標(biāo)仍然采用坐標(biāo)。本文關(guān)注的湍動(dòng)能方程,

    (2)

    其中q為湍流垂向擴(kuò)散系數(shù);s、b和分別表示湍動(dòng)能的剪切生成項(xiàng)、浮力生成項(xiàng)和耗散率;為湍流速度尺度;1為常數(shù);為海水中的摩擦速度。

    在對北太平洋海表面溫度進(jìn)行模擬的過程中, sbPOM模式空間分辨率與上文WW3模式設(shè)置保持一致, 垂向采用40個(gè)層, 在海洋上層采用較高的分辨率, 而在中、深層分辨率較低。模式初始溫、鹽場取自SODA資料月平均溫度和鹽度, 全場運(yùn)動(dòng)速度取零, 地形同樣由全球分辨率為5′×5′的ETOPO5地形插值得到。風(fēng)場和熱通量取自空間分辨率為1.875°×1.9°的NCEP再分析資料中心, sbPOM模式外模時(shí)間步長為20 s, 內(nèi)模時(shí)間步長600 s, 積分2 a, 取第二年積分結(jié)果進(jìn)行分析。

    1.3 海浪破碎的引入

    海浪破碎能夠在海表面產(chǎn)生湍動(dòng)能通量, 進(jìn)而影響海洋上混合層湍能量收支的平衡關(guān)系。近年來, 大量實(shí)驗(yàn)觀測[18-19]也都證明海浪破碎對上混合層中湍能量收支有著重要的影響。因此海洋混合層模式中需要考慮海浪破碎的影響。

    根據(jù)文獻(xiàn)[21-22], 當(dāng)考慮波浪破碎對海表面溫度的影響時(shí), 將在海表形成一向下輸入的湍動(dòng)能通量, 此時(shí), 可將(2)式給定的海表邊界條件改為(3)式:

    結(jié)合Guan等[7]提出的改進(jìn)型白冠覆蓋率耗散模型,

    2 結(jié)果分析

    2.1 海浪模式數(shù)值試驗(yàn)結(jié)果分析

    選取模式模擬有效波高作為探究WW3海浪數(shù)值試驗(yàn)?zāi)M結(jié)果準(zhǔn)確性的驗(yàn)證要素?;诿绹鴩腋?biāo)資料中心NDBC(National Data Buoy Center), 選取如圖1所示的8個(gè)浮標(biāo)站數(shù)據(jù)進(jìn)行驗(yàn)證, 驗(yàn)證的月份分別為2014年2月、5月、8月、11月。

    限于篇幅, 圖2僅給出了46075號浮標(biāo)站點(diǎn)的4個(gè)月份的有效波高模擬結(jié)果和浮標(biāo)測量結(jié)果的對比圖。從圖中可以看出, WW3海浪模式模擬有效波高值和浮標(biāo)測量值具有很好的一致性, 但是在一些風(fēng)速較大、變化較劇烈的時(shí)間段, 模擬結(jié)果稍不理想, 如, 2月8日到2月15日期間, 模式模擬結(jié)果和浮標(biāo)測量結(jié)果存在一定偏差, 出現(xiàn)這種現(xiàn)象的原因可能是NCEP再分析風(fēng)場資料時(shí)空分辨率不夠高。

    對圖1所示8個(gè)浮標(biāo)站點(diǎn)的模式模擬有效波高和浮標(biāo)測量有效波高進(jìn)行統(tǒng)計(jì)分析, 分別計(jì)算出均方根相對誤差與相關(guān)系數(shù), 其表達(dá)式分別如(6)、(7)所示。對于浮標(biāo)因某些異常原因而未能進(jìn)行測量或測量數(shù)據(jù)較少的月份, 在表1中用“NaN”表示。

    (6)

    式中,o和f分別代表有效波高的浮標(biāo)測量值和模式模擬值,為樣本點(diǎn)的個(gè)數(shù)。

    表1 模擬有效波高值的誤差分析

    Tab.1 Error analysis of the simulation significant wave heights

    根據(jù)表1中8個(gè)浮標(biāo)4個(gè)月份的模式模擬有效波高和浮標(biāo)測量有效波高的均方根相對誤差和相關(guān)系數(shù)的統(tǒng)計(jì)結(jié)果, 4個(gè)月份的平均均方根相對誤差和平均相關(guān)系數(shù)差別不大, 其中2月份平均均方根相對誤差較其他月份稍大, 為22%, 平均相關(guān)系數(shù)相對稍小, 為0.83。整體上, 模式模擬有效波高和浮標(biāo)測量有效波高的均方根相對誤差基本控制在20%左右, 最小僅為16%; 關(guān)系數(shù)基本控制在0.80以上, 最大相關(guān)系數(shù)為0.92, 且相關(guān)系數(shù)在0.80以下的僅有兩個(gè)。

    基于對圖2和表1的分析結(jié)果, 可以得出, NCEP再分析風(fēng)場資料可以作為WW3海浪模式的驅(qū)動(dòng)風(fēng)場, 其模擬出來的北太平洋有效波高等海浪要素準(zhǔn)確度較高。

    2.2 環(huán)流模式數(shù)值試驗(yàn)結(jié)果分析

    圖3分別給出了4個(gè)月份的未考慮海浪破碎模擬月平均SST和ECMWF(european center for medium range weather forecasts)再分析月平均SST之差在北太平洋的空間分布圖。從圖中可以得出, 在未引入海浪破碎的情況下, 環(huán)流模式sbPOM已經(jīng)能夠較為準(zhǔn)確的模擬出北太平洋的海表面溫度。在空間尺度上, 環(huán)流模式sbPOM對北太平洋中心海區(qū)模擬效果較好, 誤差范圍基本控制在±1℃左右, 而對近大陸海區(qū)模擬效果較差, 誤差范圍在±3℃左右。赤道地區(qū)多呈現(xiàn)模擬SST較ECMWF再分析SST偏低的趨勢。在時(shí)間尺度上, 2月份和11月份北太平洋西部海區(qū), 亞洲東岸鄰近海區(qū)呈現(xiàn)模擬SST較ECMWF再分析SST偏高的趨勢, 5月份和8月份北太平洋東部海區(qū), 美洲西岸鄰近海區(qū)呈現(xiàn)模擬SST較ECMWF再分析SST偏低的趨勢, 且赤道地區(qū)西部海區(qū)常年呈現(xiàn)出模擬SST較ECMWF再分析SST偏高的趨勢, 赤道地區(qū)東部海區(qū)常年呈現(xiàn)出模擬SST較ECMWF再分析SST偏低的趨勢。

    為初步了解海浪破碎不同月份對北太平洋海表面的作用情況, 圖4給出了海浪破碎導(dǎo)致的湍動(dòng)能通量在北太平洋不同月份的月平均空間分布圖。從圖中可以看出, 在時(shí)間尺度上, 2月份和11月份海浪破碎導(dǎo)致的湍動(dòng)能通量較大, 而5月份和8月份海浪破碎導(dǎo)致的湍動(dòng)能通量較小; 在空間尺度上, 中高緯度地區(qū)海浪破碎導(dǎo)致的湍動(dòng)能通量較大, 而在低緯度地區(qū)導(dǎo)致的湍動(dòng)能通量較小。

    選取美國全球海洋數(shù)據(jù)同化實(shí)驗(yàn)室(USA Global Ocean Data Assimilation Experiment, USGODAE) 2014年2月27日如圖5所示的9個(gè)浮標(biāo)站點(diǎn)作為驗(yàn)證點(diǎn), 將Argo浮標(biāo)資料測量SST、未考慮海浪破碎模擬SST與考慮海浪破碎模擬SST的對比結(jié)果進(jìn)行統(tǒng)計(jì)分析, 主要比較表中給出的絕對誤差和相對誤差, 其定義分別為:

    (8)

    自沉浮式剖面探測浮標(biāo)又稱Argo浮標(biāo), 因首先應(yīng)用于國際Argo計(jì)劃而被稱為Argo浮標(biāo), 專門應(yīng)用于海洋次表層溫、鹽、深剖面測量。Argo浮標(biāo)每7 d或10 d浮出一次海面, 向衛(wèi)星傳輸經(jīng)度、緯度、深度、溫度和鹽度等資料, 然后再次潛入預(yù)定深度, 一般為2 000 m左右, 其經(jīng)度、緯度、鹽度和溫度資料的精確度都為0.001。本文浮標(biāo)測量SST數(shù)據(jù)取為Argo浮標(biāo)測量首層海溫?cái)?shù)據(jù), 其測量首層水深基本在5 m左右, 具體數(shù)據(jù)見表2。

    根據(jù)表2的統(tǒng)計(jì)結(jié)果, 在sbPOM模式未引入海浪破碎時(shí), 模式本身能夠較為準(zhǔn)確的模擬出北太平洋SST, 平均絕對誤差僅為0.855。在模式引入海浪破碎后, 其模擬SST的平均值相比未引入海浪破碎時(shí)得到一定提升, 同時(shí), 平均絕對誤差由原來的0.855降為0.753, 平均相對誤差由4.26%降到3.52%。綜上所述, sbPOM模式在引入海浪破碎后, 北太平洋海域模擬SST的準(zhǔn)確性將進(jìn)一步提升。

    3 討論

    為了能從整體上更加直觀地了解波浪破碎對北太平洋海域不同季節(jié)海表面溫度模擬的影響, 圖6給出了4個(gè)月份北太平洋未考慮海浪破碎模擬SST和考慮海浪破碎后模擬SST之差分布圖。

    表2 SST北太平洋模擬與浮標(biāo)實(shí)測的比較

    Tab.2 Comparisons of the simulation SST with the results obtained from buoy data in the northern Pacific

    注: “未考慮”表示未考慮海浪破碎; “考慮”表示考慮海浪破碎

    對圖6進(jìn)行分析可以得出, 空間尺度上, 海浪破碎的引入, 對中高緯度SST模擬改善效果較為顯著, 對低緯度海區(qū)SST模擬改善效果不明顯; 對沿海地區(qū)改善效果比大洋內(nèi)部及赤道地區(qū)效果明顯。時(shí)間尺度上, 2月份和11月份海浪破碎的引入使得北太平洋大部分海區(qū)模擬SST比未引入波浪破碎時(shí)模擬SST減小, 以北太平洋西北海域最為明顯; 而5月份和8月份北太平洋北部地區(qū)模擬SST主要呈現(xiàn)比未引入波浪破碎時(shí)模擬SST增大的趨勢, 這也正好與圖3得出的結(jié)論相對應(yīng): 未考慮海浪破碎時(shí)模擬SST比ECMWF再分析SST偏高地區(qū), 如2月和11月份亞洲東部海區(qū), 海浪破碎的引入, 使得模擬SST比未考慮海浪破碎時(shí)降低; 反之, 未考慮海浪破碎時(shí)模擬SST比ECMWF再分析SST偏低地區(qū), 如5月和8月份北太平洋北部海區(qū), 海浪破碎的引入, 使得模擬SST比未考慮海浪破碎時(shí)升高。結(jié)合圖4可以發(fā)現(xiàn), 海浪破碎導(dǎo)致的湍動(dòng)能通量愈強(qiáng)的地區(qū), 海浪破碎對北太平洋SST模擬的改善效果愈明顯, 主要表現(xiàn)為使得海區(qū)模擬SST較未引入波浪破碎時(shí)減小, 而對于某些模擬SST較未引入波浪破碎時(shí)增大的海區(qū), 如5月和8月份北太平洋北部海區(qū), 從圖4可以看出其海浪破碎導(dǎo)致的湍動(dòng)能通量并不強(qiáng), 其模擬SST的增大可能與海洋環(huán)流和海氣熱通量有關(guān)。一般情況下, 海洋環(huán)流模式在考慮海浪破碎后, 會(huì)增強(qiáng)海洋表層的混合, 通常使得SST降低, 但是海表面溫度的升高或者降低同時(shí)還受到了海洋環(huán)流和海氣熱通量的影響[23 24], 這兩個(gè)過程的作用都可能使得環(huán)流模式在加入海浪破碎后模擬的SST較原來偏大。通過對比可以發(fā)現(xiàn), 在sbPOM模式對北太平洋SST進(jìn)行模擬的過程中, 海浪破碎的引入使得海洋表面混合更加均勻, 模式模擬出來的SST更加準(zhǔn)確。

    4 結(jié)論

    本文在WW3-sbPOM聯(lián)合模型的基礎(chǔ)上, 利用NCEP再分析風(fēng)場驅(qū)動(dòng)WW3海浪模式模擬出的有效波高、有效波周期和摩擦速度等海浪要素, 結(jié)合Guan等[7]提出的改進(jìn)型白冠覆蓋率耗散模型, 計(jì)算出海浪破碎導(dǎo)致湍動(dòng)能通量。通過改變環(huán)流模式sbPOM湍動(dòng)能方程的上邊界條件, 引入海浪破碎產(chǎn)生的湍動(dòng)能通量, 探究海浪破碎對北太平洋SST模擬的影響。研究表明: (1)在未考慮海浪破碎的情況下, sbPOM模式已經(jīng)能夠較為準(zhǔn)確的模擬出北太平洋SST, 且對大洋中部模擬效果較好, 誤差在±1℃左右, 對沿岸地區(qū)以及赤道地區(qū)模擬效果較差, 最大絕對誤差可達(dá)±3℃左右。(2)從空間尺度上來說, 海浪破碎的引入, 對中高緯度SST模擬效果改善較大, 對低緯度SST模擬效果改善較小。從時(shí)間尺度上來說, 海浪破碎的引入, 使得2月份和11月份北太平洋大部分海區(qū)SST比未引入海浪破碎時(shí)模擬的SST減小, 而5月份和8月份北太平洋北部地區(qū)主要呈現(xiàn)出比未引入海浪破碎時(shí)模擬SST升高的趨勢。(3)sbPOM環(huán)流模式中海浪破碎的引入, 使得海表面混合更加均勻, 對北太平洋SST的模擬相對于未引入海浪破碎之前有較大幅度提高, 模擬效果更好。

    [1] 管長龍, 張文清, 朱冬琳, 等. 上層海洋中浪致混合研究評述——研究進(jìn)展及存在問題[J]. 中國海洋大學(xué)學(xué)報(bào)自然科學(xué)版, 2014, 44(10): 20-24. Guan Changlong, Zhang Wenqing, Zhu Donglin, et al. Review of research on surface wave induced mixing in upper ocean layer: progress and existing problems[J]. Periodical of Ocean University of China, 2014, 44(10): 20-24.

    [2] 張書文, 曹瑞雪, 朱風(fēng)芹. 波浪破碎湍流混合研究綜述[J]. 物理學(xué)報(bào), 2011, 60(11): 786-791. Zhang Shuwen, Cao Ruixue, Zhu Fengqin. A review of the turbulent mixing by wind wave breaking[J]. Acta Physica Sinica, 2011, 60(11): 786-791.

    [3] Kitaigorodskii S A, Lumley J L. Wave-turbulence interactions in the upper ocean. Part I: the energy balance of the interacting fields of surface wind waves and wind-induced three-dimensional turbulence[J]. Journal of Physical Oceanography, 1983, 13(11): 1977-1987.

    [4] Terray E A, Donelan M A, Agrawal Y C, et al. Estimates of kinetic energy dissipation under breaking waves[J]. Journal of Physical Oceanography, 1996, 26(5): 792-807.

    [5] Agrawal Y C, Terray E A, Donelan M A, et al. Enhanced dissipation of kinetic energy beneath surface waves[J]. Nature, 1992, 359(6392): 219-220.

    [6] Toba Y, Kawamura H. Wind-wave coupled downward-bursting boundary layer (DBBL) beneath the sea surface[J]. Journal of Oceanography, 1996, 52(4): 409- 419.

    [7] Guan Changlong, Hu Wei, Sun Jian, et al. The whitecap coverage model from breaking dissipation parametrizations of wind waves[J]. Journal of Geophysical Research, 2007, 112(112): 395-412.

    [8] 孫群. 海浪破碎對海洋上混合層影響的數(shù)值研究[D]. 青島: 中國海洋大學(xué), 2003. Sun Qun. Numerical investigation of the effect of wave breaking on ocean surface mixed layer[D]. Qingdao: Ocean University of China, 2003.

    [9] 孫群, 管長龍, 宋金寶. 海浪破碎對海洋上混合層中湍能量收支的影響[J]. 海洋與湖沼, 2006, 37(1): 69- 74. Sun Qun, Guan Changlong, Song Jinbao. Effect of wave breaking on turbulent energy budgets in ocean surface mixed layer[J]. Chinese Journal of Oceanology and Limnology, 2006, 37(1): 69-74.

    [10] Sun Qun, Guan Changlong, Song Jinbao. Wave breaking on turbulent energy budget in the ocean surface mixed layer[J]. Chinese Journal of Oceanology and Limnology, 2008, 26(1): 9-13.

    [11] Sun Qun, Song Jinbao, Guan Changlong. Simulation of the ocean surface mixed layer under the wave breaking[J]. Acta Oceanologica Sinica, 2005, 24(3): 9-15.

    [12] Jordi A, Wang D P. sbPOM: A parallel implementation of Princenton Ocean Model[J]. Environmental Modelling & Software, 2012, 38(6): 59-61.

    [13] Tolman H L. A thrid-generation model for wind waves on slowly varying, unsteady, and inhomogeneous depths and currents[J]. Journal of Physical Oceanography, 1991, 21(6): 782-797.

    [14] 張洪生, 辜俊波, 王海龍, 等. 利用WAVEWATCH和SWAN嵌套計(jì)算珠江口附近海域的風(fēng)浪場[J]. 熱帶海洋學(xué)報(bào), 2013, 32(1): 8-17. Zhang Hongsheng, Gu Junbo, Wang Hailong, et al. Simulating wind wave field near the Pearl River Estuary with SWAN nested in WAVEWATCH[J]. Journal of Tropical Oceanography, 2013, 32(1): 8-17.

    [15] Wang Jizhao, Zhang Jie, Yang Jungang. Numerical simulation and preliminary analysis on ocean waves during Typhoon Nesat in South China Sea and adjacent areas[J]. Chinese Journal of Oceanology & Limnology, 2014, 32(3): 665-680.

    [16] Cui Hong, He Hailun, Liu Xiaohui, et al. Effect of oceanic current on typhoon-wave modeling in the East China Sea[J]. Chinese Physics B, 2012, 21(10): 592-599.

    [17] 劉欣, 韋駿. 熱帶氣旋與海洋暖渦間的海-氣相互作用[J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 50(3): 456- 466. Liu Xin, Wei Jun. Air-Sea interaction between Tropical Cyclone and Ocean Warm Core Ring[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2014, 50(3): 456-466.

    [18] Agrawal Y C, Terray E A, Donelan M A, et al. Enhanced dissipation of kinetic energy beneath surface waves[J]. Nature, 1992, 359(6392): 219-220.

    [19] Terray E A, Donelan M A, Agrawal Y C, et al. Estimates of kinetic energy dissipation under breaking waves[J]. Journal of Physical Oceanography, 1996, 26(5): 792-807.

    [20] Alari V, Staneva J, ?yvind Breivik, et al. Surface wave effects on water temperature in the Baltic Sea: simulations with the coupled NEMO-WAM model[J]. Ocean Dynamics, 2016, 66(8): 917-930.

    [21] Mellor G L, Yamada T. A hierarchy of turbulence closure models for planetary boundary layers.[J]. Journal of the Atmospheric Sciences, 1974, 31(7): 1791-1806.

    [22] Mellor G L, Yamada T. Development of a turbulence closure model for geophysical fluid problems[J]. Reviews of Geophysics and Space Physics, 1982, 20(4): 851-875.

    [23] 宋振亞, 喬方利, 楊永增, 等. 波致混合對熱帶太平洋海氣耦合模式中冷舌模擬的改進(jìn)[J]. 自然科學(xué)進(jìn)展, 2006, 16(9): 1138-1145. Song Zhenya, Qiao Fangli, Yang Yongzeng, et al. The improvement of wave-induced mixing on the simulation of cold tongue in the tropical Pacific Ocean with sea-air coupled model[J]. Progress in Natural Science, 2006, 16(9): 1138-1145.

    [24] 宋振亞, 喬方利, 雷曉燕, 等. 大氣-海浪-海洋環(huán)流耦合數(shù)值模式的建立及北太平洋SST模擬[J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2007, 22(5): 543-548. Song Zhenya, Qiao Fangli, Lei Xiaoyan, et al. The establishment of an atmosphere-wave-ocean circulation coupled numerical model and its application in the North Pacific SST simulation[J]. Journal of Hydrodynamics, 2007, 22(5): 543-548.

    [25] Blumberg A F, Mellor G L. A description of a three dimensional coastal ocean circulation model[C]// Heaps N S. Three-Dimensional Coastal Ocean Models, vol 4. Washington DC: American Geophysical Union, 1987: 1-16.

    [26] 劉子龍, 史劍, 蔣國榮, 等. 海浪攪拌混合對北太平洋海表面溫度模擬的影響[J]. 海洋科學(xué), 2016, 40(12): 131-137. Liu Zilong, Shi Jian, Jiang Guorong, et al.Influence of wave-induced mixing on a sea surface temperature simulation of the North Pacific[J]. Marine Sciences, 2016, 40(12)131-137.

    [27] Thomson J, Schwendeman M S, Zippel S F, et al. Wave breaking turbulence in the ocean surface layer[J]. Journal of Physical Oceanography, 2016, 46(6): 1866- 1882.

    [28] Breivik ?, Mogensen K, Bidlot J, et al. Surface wave effects in the NEMO ocean model: Forced and coupled experiments[J]. Journal of Geophysical Research Oceans, 2015, 120(4): 135-143.

    (本文編輯: 劉珊珊)

    Effect of wave breaking on the sea surface temperature simulation in the North Pacific

    LIU Zi-long, SHI Jian, JIANG Guo-rong

    (College of Meteorology and Oceanography, PLA University of Science and Technology, Nanjing 211101, China)

    In this paper, wave parameters in the North Pacific were simulated by the WAVEWATCH III and were combined with NDBC buoy data verification. We found that the simulated significant wave height is highly precise. Based on the modified model of whitecap coverage dissipation, turbulence kinetic energy (TKE) fluxes due to wave breaking were calculated using the simulated significant wave height, wave period, and friction velocity. The influence of surface wave breaking could be introduced into ocean model sbPOM by modifying the existing surface boundary condition of the TKE equation and specifying its input. The effects of wave breaking on the SST simulation in the North Pacific were also investigated. The results indicated that wave breaking can improve the accuracy of the ocean model sbPOM on the North Pacific sea surface temperature simulation and it plays an important role in providing accurate lower boundary conditions in atmospheric models.

    wave breaking; sea surface temperature; WAVEWATCH III; sbPOM

    Sep.4, 2016

    [National Natural Science Foundation of China, No.41676014]

    P732.7

    A

    1000-3096(2017)03-0122-08

    10.11759/hykx20160904001

    2016-09-04;

    2017-01-18

    國家自然科學(xué)基金項(xiàng)目(41676014)

    劉子龍(1991-), 男, 湖南常德人, 碩士研究生, 研究方向: 海洋動(dòng)力學(xué)與數(shù)值模擬, 手機(jī): 18761683583, E-mail: 941117242@qq.com; 史劍(1981-), 通信作者, 男, 江蘇揚(yáng)州人, 副教授, 研究方向: 海洋動(dòng)力學(xué)與數(shù)值模擬, 手機(jī): 13813388338

    猜你喜歡
    表面溫度浮標(biāo)環(huán)流
    受了委屈的浮標(biāo)君
    受了委屈的浮標(biāo)君
    受了委屈的浮標(biāo)君
    家教世界(2023年7期)2023-03-22 12:11:24
    受了委屈的浮標(biāo)君
    家教世界(2023年4期)2023-03-04 07:31:28
    內(nèi)環(huán)流控溫技術(shù)應(yīng)用實(shí)踐與發(fā)展前景
    結(jié)合注意力機(jī)制的區(qū)域型海表面溫度預(yù)報(bào)算法
    熱鹽環(huán)流方程全局弱解的存在性
    謎底大揭秘
    熱電池新型隔熱結(jié)構(gòu)設(shè)計(jì)及表面溫度研究
    Kerr型中子星與黑洞表面溫度分布的研究
    色视频在线一区二区三区| 美女视频免费永久观看网站| 成人18禁高潮啪啪吃奶动态图 | 亚洲一级一片aⅴ在线观看| 亚洲第一区二区三区不卡| h视频一区二区三区| 国产精品一二三区在线看| 一区二区三区四区激情视频| 草草在线视频免费看| 永久网站在线| 欧美高清成人免费视频www| 极品教师在线视频| 97热精品久久久久久| 国产国拍精品亚洲av在线观看| 亚洲av中文av极速乱| 看免费成人av毛片| 久久精品国产鲁丝片午夜精品| 国产精品嫩草影院av在线观看| 国产视频首页在线观看| 美女主播在线视频| av.在线天堂| 另类亚洲欧美激情| 国产欧美日韩一区二区三区在线 | 女性生殖器流出的白浆| 国产精品麻豆人妻色哟哟久久| 狂野欧美激情性bbbbbb| 日本爱情动作片www.在线观看| 国产成人91sexporn| 国产爱豆传媒在线观看| 高清视频免费观看一区二区| 直男gayav资源| 国产人妻一区二区三区在| 少妇裸体淫交视频免费看高清| 国产成人91sexporn| 91午夜精品亚洲一区二区三区| 日韩,欧美,国产一区二区三区| 老师上课跳d突然被开到最大视频| 国产精品99久久99久久久不卡 | 亚洲aⅴ乱码一区二区在线播放| 亚洲精品色激情综合| 亚洲真实伦在线观看| av天堂中文字幕网| av视频免费观看在线观看| 久久久久久久国产电影| 亚洲av免费高清在线观看| 偷拍熟女少妇极品色| 亚洲国产精品专区欧美| 日韩中文字幕视频在线看片 | 久久久久久久久久久丰满| h视频一区二区三区| 亚洲av成人精品一区久久| 国产中年淑女户外野战色| av黄色大香蕉| 久久精品国产鲁丝片午夜精品| 韩国av在线不卡| 免费观看的影片在线观看| 久久99热这里只频精品6学生| 2018国产大陆天天弄谢| 成人国产麻豆网| 免费高清在线观看视频在线观看| 日本免费在线观看一区| 国产一区二区三区综合在线观看 | 成人18禁高潮啪啪吃奶动态图 | 国产免费视频播放在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 黄片wwwwww| 亚洲激情五月婷婷啪啪| 全区人妻精品视频| 久久鲁丝午夜福利片| av专区在线播放| 中文在线观看免费www的网站| 欧美3d第一页| 亚洲国产精品成人久久小说| 嫩草影院入口| 伦理电影免费视频| 亚洲精品色激情综合| 日韩精品有码人妻一区| av网站免费在线观看视频| 十八禁网站网址无遮挡 | 亚洲三级黄色毛片| 亚洲人成网站在线播| 亚洲精品成人av观看孕妇| 中文字幕精品免费在线观看视频 | 一级毛片aaaaaa免费看小| 国产精品一二三区在线看| 日本欧美国产在线视频| 一区二区三区免费毛片| 搡老乐熟女国产| 国内少妇人妻偷人精品xxx网站| 国产精品一区二区三区四区免费观看| tube8黄色片| 日韩 亚洲 欧美在线| 狂野欧美激情性xxxx在线观看| 大码成人一级视频| 国产精品一及| 久久99热这里只有精品18| av免费观看日本| 亚洲美女黄色视频免费看| 男女啪啪激烈高潮av片| 亚洲成人中文字幕在线播放| 又粗又硬又长又爽又黄的视频| 一级毛片久久久久久久久女| 国产在视频线精品| 插阴视频在线观看视频| 99久久精品国产国产毛片| 亚洲国产色片| 国产中年淑女户外野战色| 色婷婷久久久亚洲欧美| 我要看黄色一级片免费的| 久久精品夜色国产| 一区二区三区免费毛片| 亚洲精品国产成人久久av| 各种免费的搞黄视频| 日韩精品有码人妻一区| 国产黄频视频在线观看| 22中文网久久字幕| 亚洲第一区二区三区不卡| 久久久久久久久久久免费av| 国产成人一区二区在线| 亚洲丝袜综合中文字幕| tube8黄色片| 久久精品国产a三级三级三级| 欧美3d第一页| 99热6这里只有精品| 成年人午夜在线观看视频| 欧美xxxx黑人xx丫x性爽| av在线播放精品| 亚洲av电影在线观看一区二区三区| 一级片'在线观看视频| 我的女老师完整版在线观看| 黄片wwwwww| 国产 一区精品| 2022亚洲国产成人精品| 永久免费av网站大全| 国产有黄有色有爽视频| 亚洲成色77777| 九草在线视频观看| 精品少妇黑人巨大在线播放| 欧美日本视频| 国产午夜精品一二区理论片| 免费人成在线观看视频色| 少妇精品久久久久久久| 高清在线视频一区二区三区| 国产欧美日韩一区二区三区在线 | 天堂中文最新版在线下载| 亚洲成色77777| 久久鲁丝午夜福利片| 人人妻人人爽人人添夜夜欢视频 | 午夜视频国产福利| 在线观看人妻少妇| 精品人妻一区二区三区麻豆| 亚洲美女搞黄在线观看| 国产欧美日韩一区二区三区在线 | 自拍偷自拍亚洲精品老妇| 亚洲精品亚洲一区二区| 成人毛片a级毛片在线播放| 国产免费福利视频在线观看| 肉色欧美久久久久久久蜜桃| 欧美区成人在线视频| 国产男女超爽视频在线观看| h日本视频在线播放| 九草在线视频观看| 人妻制服诱惑在线中文字幕| 伊人久久精品亚洲午夜| 十八禁网站网址无遮挡 | 欧美变态另类bdsm刘玥| 80岁老熟妇乱子伦牲交| 久久久精品免费免费高清| av专区在线播放| 人妻 亚洲 视频| 99久久精品一区二区三区| 久久婷婷青草| 精品一品国产午夜福利视频| h视频一区二区三区| 人妻 亚洲 视频| 免费看不卡的av| 日日摸夜夜添夜夜添av毛片| 日韩成人伦理影院| 久久久久久久久久成人| 不卡视频在线观看欧美| 久久鲁丝午夜福利片| 大香蕉久久网| av天堂中文字幕网| 国产av码专区亚洲av| 天堂中文最新版在线下载| 国产成人精品久久久久久| 人妻夜夜爽99麻豆av| 国产爱豆传媒在线观看| 亚洲精品视频女| 自拍偷自拍亚洲精品老妇| 日本av免费视频播放| 日韩伦理黄色片| 国产精品av视频在线免费观看| 亚洲国产精品专区欧美| 一本色道久久久久久精品综合| 自拍欧美九色日韩亚洲蝌蚪91 | 国产欧美另类精品又又久久亚洲欧美| 黑人高潮一二区| 麻豆国产97在线/欧美| 亚洲国产欧美人成| 久久影院123| 精品午夜福利在线看| 久久精品久久久久久噜噜老黄| 欧美三级亚洲精品| 天堂中文最新版在线下载| 国产永久视频网站| 亚洲真实伦在线观看| 日本午夜av视频| 亚洲国产av新网站| 亚洲自偷自拍三级| 黑人猛操日本美女一级片| 一级二级三级毛片免费看| 亚洲,一卡二卡三卡| 最后的刺客免费高清国语| 久久久色成人| 韩国av在线不卡| 亚洲av不卡在线观看| 高清在线视频一区二区三区| 国产免费又黄又爽又色| 国产精品国产av在线观看| 久久精品国产亚洲av涩爱| 欧美最新免费一区二区三区| 久久人妻熟女aⅴ| 男人舔奶头视频| 日韩中字成人| 我的女老师完整版在线观看| 国产精品免费大片| 少妇裸体淫交视频免费看高清| 天天躁日日操中文字幕| 国产精品嫩草影院av在线观看| 免费黄色在线免费观看| 久热久热在线精品观看| 精品亚洲成a人片在线观看 | 色吧在线观看| 欧美另类一区| 高清不卡的av网站| 日本欧美国产在线视频| 99热这里只有是精品在线观看| 久久久久久久国产电影| 一级爰片在线观看| 91精品一卡2卡3卡4卡| 亚洲色图av天堂| 高清黄色对白视频在线免费看 | 精品人妻偷拍中文字幕| 一区在线观看完整版| 亚洲综合色惰| 男女国产视频网站| 伊人久久精品亚洲午夜| 五月开心婷婷网| 身体一侧抽搐| 亚洲真实伦在线观看| 99久久精品热视频| 国产亚洲午夜精品一区二区久久| 日本av手机在线免费观看| 亚洲综合精品二区| 久久久久性生活片| 久久精品国产鲁丝片午夜精品| 亚洲丝袜综合中文字幕| 免费看光身美女| 永久网站在线| 国产精品无大码| 各种免费的搞黄视频| 美女xxoo啪啪120秒动态图| av国产久精品久网站免费入址| 亚洲欧美一区二区三区国产| 国产爽快片一区二区三区| 成人一区二区视频在线观看| 亚洲最大成人中文| 国产精品一区二区性色av| 久久国产乱子免费精品| 大码成人一级视频| 黄色怎么调成土黄色| 国内精品宾馆在线| 国产精品久久久久久精品古装| 少妇精品久久久久久久| 麻豆成人av视频| av.在线天堂| 精品人妻视频免费看| 在线观看三级黄色| av在线播放精品| a 毛片基地| 丝瓜视频免费看黄片| 国产在线男女| 国产高清不卡午夜福利| 蜜臀久久99精品久久宅男| av免费在线看不卡| 成人漫画全彩无遮挡| 美女cb高潮喷水在线观看| 欧美3d第一页| 精品一区二区三卡| 国产欧美亚洲国产| 国内少妇人妻偷人精品xxx网站| 草草在线视频免费看| 噜噜噜噜噜久久久久久91| 最近手机中文字幕大全| 国产一级毛片在线| 秋霞伦理黄片| 亚洲成人av在线免费| av播播在线观看一区| 国产亚洲91精品色在线| 国产成人免费无遮挡视频| 在线免费观看不下载黄p国产| 精品亚洲成a人片在线观看 | 有码 亚洲区| 一区二区三区免费毛片| 国产69精品久久久久777片| 黄色欧美视频在线观看| 久久鲁丝午夜福利片| 日韩av不卡免费在线播放| 日韩制服骚丝袜av| 天堂8中文在线网| 水蜜桃什么品种好| 最近最新中文字幕大全电影3| 欧美激情国产日韩精品一区| 91久久精品电影网| 国产精品一及| 国产免费视频播放在线视频| 麻豆精品久久久久久蜜桃| 国产色婷婷99| 网址你懂的国产日韩在线| 在线观看免费日韩欧美大片 | 国产一区有黄有色的免费视频| 黄色配什么色好看| av在线观看视频网站免费| 一边亲一边摸免费视频| 午夜老司机福利剧场| 欧美日韩亚洲高清精品| 久久亚洲国产成人精品v| 亚洲第一区二区三区不卡| 搡女人真爽免费视频火全软件| 久久综合国产亚洲精品| 日韩 亚洲 欧美在线| 亚洲第一av免费看| 99精国产麻豆久久婷婷| 久久精品国产自在天天线| 欧美性感艳星| 51国产日韩欧美| 美女中出高潮动态图| 国产av码专区亚洲av| 精品国产乱码久久久久久小说| 美女xxoo啪啪120秒动态图| 国产精品99久久99久久久不卡 | 久久婷婷青草| 欧美xxⅹ黑人| 十分钟在线观看高清视频www | 精华霜和精华液先用哪个| 毛片一级片免费看久久久久| 中文字幕制服av| 性高湖久久久久久久久免费观看| 亚洲自偷自拍三级| 狂野欧美激情性bbbbbb| 大陆偷拍与自拍| 成人特级av手机在线观看| 午夜精品国产一区二区电影| 边亲边吃奶的免费视频| 久久99热这里只频精品6学生| 国产精品一区www在线观看| 久久久久性生活片| 亚洲国产成人一精品久久久| 国产精品.久久久| 亚洲国产av新网站| 精品一区二区三区视频在线| 国产成人a区在线观看| 日本黄大片高清| 精品久久国产蜜桃| 国产 精品1| 亚洲精华国产精华液的使用体验| 日产精品乱码卡一卡2卡三| 亚洲国产精品一区三区| 久久99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线| 偷拍熟女少妇极品色| 久久久a久久爽久久v久久| 99精国产麻豆久久婷婷| 天天躁夜夜躁狠狠久久av| 最黄视频免费看| 草草在线视频免费看| 国产亚洲91精品色在线| 亚洲欧美中文字幕日韩二区| av视频免费观看在线观看| 亚洲真实伦在线观看| 色视频在线一区二区三区| 在线亚洲精品国产二区图片欧美 | 丰满少妇做爰视频| 国产一区有黄有色的免费视频| 国产女主播在线喷水免费视频网站| 国产精品欧美亚洲77777| 三级国产精品片| 十八禁网站网址无遮挡 | 草草在线视频免费看| 亚洲综合精品二区| 成年女人在线观看亚洲视频| 久久久午夜欧美精品| 视频中文字幕在线观看| 欧美高清性xxxxhd video| 国产精品久久久久久久电影| 2018国产大陆天天弄谢| 免费在线观看成人毛片| 免费看光身美女| 日本vs欧美在线观看视频 | 亚洲精品日本国产第一区| 另类亚洲欧美激情| 欧美日韩在线观看h| 亚洲成人手机| 黄色欧美视频在线观看| 久久久午夜欧美精品| av在线播放精品| 国产大屁股一区二区在线视频| 身体一侧抽搐| 国产亚洲最大av| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 亚洲国产欧美人成| 美女视频免费永久观看网站| 国产男女超爽视频在线观看| 久久人人爽人人爽人人片va| a级毛色黄片| 欧美最新免费一区二区三区| 午夜激情久久久久久久| 日本黄色片子视频| 老司机影院毛片| 日本猛色少妇xxxxx猛交久久| 啦啦啦啦在线视频资源| av天堂中文字幕网| 精品久久国产蜜桃| 香蕉精品网在线| 久久人妻熟女aⅴ| 国产精品精品国产色婷婷| 欧美成人午夜免费资源| 亚洲av日韩在线播放| 久久综合国产亚洲精品| 又爽又黄a免费视频| 男女边摸边吃奶| 国产免费视频播放在线视频| 国产免费又黄又爽又色| 亚洲美女黄色视频免费看| 性色avwww在线观看| 另类亚洲欧美激情| 少妇人妻 视频| 水蜜桃什么品种好| www.av在线官网国产| 在线观看三级黄色| 亚洲精品中文字幕在线视频 | 成人美女网站在线观看视频| 国产精品福利在线免费观看| 97热精品久久久久久| 国产又色又爽无遮挡免| 在线观看一区二区三区| 在线免费十八禁| 免费大片黄手机在线观看| 国产亚洲91精品色在线| 久热久热在线精品观看| 欧美xxⅹ黑人| 日本爱情动作片www.在线观看| 亚洲精品第二区| 国产探花极品一区二区| 五月玫瑰六月丁香| 美女视频免费永久观看网站| 日本黄大片高清| 亚洲色图综合在线观看| 国产色婷婷99| 最黄视频免费看| 黄色视频在线播放观看不卡| 性高湖久久久久久久久免费观看| 麻豆国产97在线/欧美| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 一区二区三区免费毛片| 91精品一卡2卡3卡4卡| 丰满人妻一区二区三区视频av| 美女高潮的动态| 亚洲欧美一区二区三区黑人 | 久久久久久九九精品二区国产| 在线观看免费高清a一片| 国产午夜精品一二区理论片| 精品一区二区免费观看| 国产精品一二三区在线看| 天堂中文最新版在线下载| 最近中文字幕高清免费大全6| 日韩电影二区| 国产伦精品一区二区三区视频9| 在线免费观看不下载黄p国产| 美女脱内裤让男人舔精品视频| 国产综合精华液| 久久婷婷青草| 国产成人91sexporn| 免费看av在线观看网站| 边亲边吃奶的免费视频| 亚洲婷婷狠狠爱综合网| 91狼人影院| 九色成人免费人妻av| 国产精品一区二区性色av| 免费观看性生交大片5| h视频一区二区三区| www.色视频.com| 男人狂女人下面高潮的视频| 男的添女的下面高潮视频| 一区二区三区免费毛片| 高清视频免费观看一区二区| 色哟哟·www| 老师上课跳d突然被开到最大视频| 国产真实伦视频高清在线观看| 人妻少妇偷人精品九色| 熟女电影av网| 久久国产亚洲av麻豆专区| 日韩成人伦理影院| 九色成人免费人妻av| 中文字幕人妻熟人妻熟丝袜美| 婷婷色综合大香蕉| 久久精品国产亚洲网站| 1000部很黄的大片| 精品久久久精品久久久| www.av在线官网国产| 国产精品精品国产色婷婷| 精品视频人人做人人爽| 丰满人妻一区二区三区视频av| 亚洲精品国产av蜜桃| 国产精品成人在线| 久久精品夜色国产| av在线蜜桃| 午夜福利影视在线免费观看| 女人久久www免费人成看片| 舔av片在线| 国产av一区二区精品久久 | 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 成人漫画全彩无遮挡| a级一级毛片免费在线观看| 亚洲国产欧美在线一区| videos熟女内射| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线 | 久久久久久久国产电影| 国产午夜精品久久久久久一区二区三区| 国产黄色视频一区二区在线观看| 男女免费视频国产| 春色校园在线视频观看| 99精国产麻豆久久婷婷| 天天躁日日操中文字幕| 国产精品一区二区三区四区免费观看| 视频中文字幕在线观看| 尾随美女入室| 青春草亚洲视频在线观看| 少妇的逼水好多| 亚洲国产欧美人成| 免费人成在线观看视频色| xxx大片免费视频| 亚洲人成网站高清观看| 色吧在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产伦精品一区二区三区视频9| 亚洲精品国产av蜜桃| 亚洲一区二区三区欧美精品| 午夜老司机福利剧场| 不卡视频在线观看欧美| 欧美一区二区亚洲| 国产伦精品一区二区三区四那| 又粗又硬又长又爽又黄的视频| 日本欧美视频一区| 国产男人的电影天堂91| 精品久久久噜噜| 欧美少妇被猛烈插入视频| 中国美白少妇内射xxxbb| 校园人妻丝袜中文字幕| 久久久久国产网址| 卡戴珊不雅视频在线播放| 国产在线免费精品| videossex国产| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 边亲边吃奶的免费视频| 狂野欧美白嫩少妇大欣赏| 日韩一本色道免费dvd| 亚洲人与动物交配视频| 日韩人妻高清精品专区| 国产精品嫩草影院av在线观看| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 免费人成在线观看视频色| 日日撸夜夜添| 中文字幕免费在线视频6| 午夜福利在线在线| 国产淫片久久久久久久久| 精品国产一区二区三区久久久樱花 | 亚洲电影在线观看av| 免费观看无遮挡的男女| 中文资源天堂在线| 色哟哟·www| 亚洲第一av免费看| 日日撸夜夜添| 亚洲不卡免费看| 亚洲成人中文字幕在线播放| 青春草国产在线视频| 在线观看av片永久免费下载| 热99国产精品久久久久久7| 久久国产精品大桥未久av | 免费不卡的大黄色大毛片视频在线观看| 精品99又大又爽又粗少妇毛片| 欧美日韩在线观看h| 免费在线观看成人毛片| 在线观看av片永久免费下载| 色婷婷久久久亚洲欧美| 一级毛片黄色毛片免费观看视频| 国产精品久久久久久av不卡| 不卡视频在线观看欧美| 丰满迷人的少妇在线观看| 少妇的逼水好多| 亚洲欧洲日产国产| 国产淫语在线视频| 啦啦啦在线观看免费高清www| 一个人免费看片子| 一级毛片 在线播放| 日日啪夜夜爽| 精品久久久久久电影网| 亚洲美女搞黄在线观看|