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

    渤海夏季第一斜壓羅斯貝變形半徑的計算與分析

    2015-11-28 11:08:24張曉璐熊學(xué)軍
    海洋科學(xué)進(jìn)展 2015年3期
    關(guān)鍵詞:斜壓浮性渤海

    張曉璐,熊學(xué)軍*

    (1.國家海洋局 第一海洋研究所,山東 青島266061;2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點(diǎn)實(shí)驗(yàn)室,山東 青島266061)

    渤海,是一個深入中國大陸的半封閉淺海,北、西、南三面被陸地包圍,僅通過東面的渤海海峽與黃海溝通。注入渤海的沿岸河流眾多,北有遼河、大遼河、大凌河,西有海河、灤河,南有黃河、小清河,因此形成了渤海獨(dú)特的水文特征及動力結(jié)構(gòu)。相較于冬季,渤海夏季的風(fēng)速小,主要為熱鹽流。近年來,大量的數(shù)值模擬工作及調(diào)查研究表明,夏季渤海多存在一些局地范圍的環(huán)流,如渤海中部的順時針環(huán)流,長興島以西的遼東灣海域的逆時針環(huán)流,以及在渤海南部的萊州灣存在弱的順時針環(huán)流等[1]。

    在海洋環(huán)流的理論研究和數(shù)值模擬中,羅斯貝變形半徑是一個極其重要的水平尺度,它描述了地球旋轉(zhuǎn)對流體運(yùn)動的影響。這個概念,最早是由Carl-Gustav Rossby于1938年在研究地轉(zhuǎn)調(diào)整中提出的。若海面高度變化在調(diào)整過程中的尺度與羅斯貝變形半徑相比為小量,則壓強(qiáng)梯度非常大,重力支配運(yùn)動的狀態(tài);但是,當(dāng)海面高度的變化發(fā)展到具有羅斯貝變形半徑尺度的范圍時,科氏力項(xiàng)變得同壓強(qiáng)梯度項(xiàng)同等重要[2]。簡言之,即在羅斯貝變形半徑這個距離尺度上,科氏力使自由面變形的趨勢與壓強(qiáng)梯度力使自由面復(fù)原的趨勢相平衡。在海洋中,斜壓羅斯貝變形半徑是一種常與邊界現(xiàn)象(諸如邊界流、鋒面)以及渦旋有聯(lián)系的自然尺度,是設(shè)置數(shù)值模式參數(shù)的重要參考。有研究表明,在模擬中尺度運(yùn)動時,網(wǎng)格大小應(yīng)當(dāng)是羅斯貝變形半徑的一半,或者是其三分之一[3]。

    國外對斜壓羅斯貝變形半徑的研究已有不少。Emery等[4]在1984年用氣候態(tài)的溫鹽剖面數(shù)據(jù)計算了北太平洋和北大西洋的斜壓羅斯貝變形半徑;Houry等[5]在1987年用NODC的溫鹽數(shù)據(jù)計算了南大西洋的斜壓羅斯貝變形半徑;特別地,Chelton等[6]在1996年利用NODC的1°×1°的氣候態(tài)的溫鹽數(shù)據(jù)計算了全球大洋的第一斜壓羅斯貝變形半徑以及對應(yīng)的全球范圍的第一斜壓重力波相速度,對后來學(xué)者的研究工作起了重要的啟示和參考作用。除了大洋中,在內(nèi)海也有針對斜壓羅斯貝變形半徑的專門研究,如Fennel等[7]、Alenius等[8]以及Osinsik等[9]曾分別利用航測溫鹽剖面數(shù)據(jù)計算了波羅的海、芬蘭灣的斜壓羅斯貝變形半徑。但在國內(nèi),有關(guān)斜壓羅斯貝變形半徑的研究工作不多,甘子鈞和蔡樹群[9],蔡樹群等[10],曾用Levitus和NODC標(biāo)準(zhǔn)層上的氣候態(tài)的溫鹽資料計算過南海斜壓羅斯貝變形半徑的地理及季節(jié)變化以及相應(yīng)的第一斜壓重力波相速度。

    在渤海進(jìn)行第一斜壓羅斯貝變形半徑的計算和研究尚未發(fā)現(xiàn)。本文在觀測數(shù)據(jù)基礎(chǔ)上,分別利用WKB近似和數(shù)值解法,對渤海夏季第一斜壓羅斯貝變形半徑進(jìn)行了計算和分析,希望有助于渤海環(huán)流、中尺度渦旋,以及與此相關(guān)的其他海洋現(xiàn)象的研究。

    1 資料

    本文所用資料為渤海2006年夏季的溫度和鹽度大面觀測資料,觀測儀器為日本ALEC公司的AAQ1183系列的CTD,觀測時間為2006-07-18—08-12。為方便對下文計算結(jié)果分析對比,將整個渤海劃分為6個區(qū)域,分別為渤海東部海區(qū)、遼東灣海域、北戴河附近海域、天津附近海域、黃河口三角洲海域以及萊州灣海域。觀測站位分布,渤海區(qū)劃及渤海水深如圖1所示。

    2 計算方法

    羅斯貝變形半徑在數(shù)學(xué)表達(dá)上,可以理解為:在地球旋轉(zhuǎn)特征周期這一時間尺度上,長重力波傳播的特征距離[11]。在緯度大于5°的海域,羅斯貝變形半徑可以寫為:

    式中,f=2Ωsin(θ)是科氏參數(shù);Ω為地球自轉(zhuǎn)角速度;θ為地理緯度。當(dāng)n=0時,c0為海洋中正壓波的相速度,對應(yīng)的R0為正壓羅斯貝變形半徑;當(dāng)n=1時,c1為第一斜壓重力波相速度,對應(yīng)的R1為第一斜壓羅斯貝變形半徑。第一斜壓羅斯貝變形半徑在海洋環(huán)流理論方面起重要作用,是中尺度渦、沿岸射流和赤道流等的水平尺度[10],故本文主要研究第一斜壓重力波相速度以及對應(yīng)的第一斜壓羅斯貝變形半徑。通常如不特別說明,斜壓羅斯貝變形半徑均指第一斜壓羅斯貝變形半徑。

    在不考慮底摩擦和風(fēng)應(yīng)力,引入平底、海面剛蓋近似,及Boussinesq近似的假設(shè)下,穩(wěn)定層化海洋水平大尺度波動的垂直結(jié)構(gòu)滿足:

    式中,H為水深;φ(z)為垂向速度分離變量后只與z有關(guān)的垂向函數(shù);N(z)為Brunt-V?is?l?頻率,亦稱浮性頻率,忽略壓縮性的影響,其表達(dá)式為:

    式中,ρ表示密度;g表示重力加速度。

    方程(2)及邊界條件(3)組成了典型的Sturm-Liouville問題,求解羅斯貝變形半徑,即轉(zhuǎn)換成求解Sturm-Liouville本征值的問題[2]。對于連續(xù)層化的海洋,方程(2)對應(yīng)的本征值是一組有序、非負(fù)的數(shù)列:。對應(yīng)的本征函數(shù)為φ1(z),φ2(z),φ3(z)…,下標(biāo)表示模態(tài)數(shù)(mode number)。cn表征第n模態(tài)斜壓重力波相速度,由c1和科氏參數(shù)f便可確定第一斜壓羅斯貝變形半徑R1。

    2.1 WKB近似

    在求解由方程(2)及邊界條件(3)組成的Sturm-Liouville問題中,WKB近似是一種快捷有效的方法[6]。其基本思想是把φ(z)寫成下面的小參數(shù)冪級數(shù)的指數(shù)形式:

    根據(jù)本文的情況,在小參數(shù)冪級數(shù)中引入虛數(shù)單位i。將式(5)代入方程(2)中,令ε的同次冪系數(shù)相等,可得遞推方程組:

    通過遞推方程組求解χ0(z),χ1(z),代入到式(5)中,只保留冪級數(shù)顯著的前兩項(xiàng),即可得出φ(z)的WKB近似解:

    式中,A,B為任意常系數(shù),代入邊界條件(3)得:

    當(dāng)n=1時,通過式(8)和式(9)便可求解出第一斜壓重力波相速度c1和第一斜壓羅斯貝變形半徑R1。

    2.2 數(shù)值解法

    數(shù)值解法是將本征值問題離散化后求解的過程。在相鄰離散深度zk、zk+1上,利用中間密度梯度法(neutral density gradient method)[6],將現(xiàn)場密度絕熱調(diào)整到中間深度ζk=zk+1/2(圖2),然后估計浮性頻率的平方N2(z):

    圖2 數(shù)值解法中在離散深度上求解浮性頻率的示意圖Fig.2 Schematic diagram of vertical discretization in numerical method

    據(jù)Chelton等[6]的建議,為了求羅斯貝變形半徑,如果得到的N2(ζk)為負(fù)值,用相鄰較淺深度ζk-1上的N2(ζk-1)代替;若N2(ζ1)為負(fù)值,令N2(ζ1)為10-8。

    為了得到數(shù)值解,在深度ζk上離散本征值方程。首先,φ(z)在z=ζk(k=1,…,n-1)上進(jìn)行泰勒展開:

    通過式(11)求取φ(ζk-1)、φ(ζk+1),消除φ′(ζk)項(xiàng)后有:

    其中,

    把式(11)代入到方程(2)中,得到n-1個離散方程,移項(xiàng)整理得:

    式中,

    根據(jù)本文的情況,定義邊界條件為:

    至此,將求解Sturm-Liouville本征值問題轉(zhuǎn)換成求n-1階矩陣特征值的問題[6]:

    式中,

    由式(18)可解出n-1個特征值μ1,…,μn-1,對應(yīng)第一個特征值μ1,通過式(16)解出第一斜壓重力波相速度c1和第一斜壓羅斯貝變形半徑R1。

    3 結(jié)果分析

    圖3和圖4分別給出了通過WKB近似和數(shù)值解法求出的渤海夏季第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑的分布。兩種算法下得到的渤海夏季第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑的分布格局相似。斜壓重力波相速度與斜壓羅斯貝變形半徑依賴于浮性頻率的積分深度,即與實(shí)際水深有密切關(guān)系,故它們的分布與渤海水深分布大體一致,等值線呈耳狀分布。在北戴河和長興島附近海域,斜壓重力波相速度與斜壓羅斯貝變形半徑最大,然后向岸逐漸遞減。在渤海中部等值線發(fā)生明顯彎曲,在靠近渤海海峽內(nèi)側(cè)斜壓重力波相速度與斜壓羅斯貝變形半徑出現(xiàn)一個小值,而與之對應(yīng)在渤海中部以西出現(xiàn)一個大值,導(dǎo)致等值線呈NE-SW向傾斜的“S”狀。

    圖3 WKB近似結(jié)果Fig.3 Result of WKB method

    圖4 數(shù)值方法計算結(jié)果Fig 4 Result of numerical method

    表1 WKB近似計算得出的第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑Table 1 The first baroclinic gravity wave phase speed and Rossby deformation radius calculated by WKB method

    表2 數(shù)值方法計算得出的第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑Table 2 The first baroclinic gravity wave phase speed and Rossby deformation radius calculated by numerical method

    在WKB近似下,渤海夏季第一斜壓羅斯貝變形半徑的最大值出現(xiàn)在長興島附近海域,為3.13km,對應(yīng)的斜壓重力波相速度最大為0.28m/s;在萊州灣海域斜壓羅斯貝變形半徑和對應(yīng)的斜壓重力波相速度最大值僅為0.95km和0.08m/s;在遼東灣海域、北戴河附近海域、天津附近海域和黃河口三角洲海域,斜壓羅斯貝變形半徑和斜壓重力波相速度的最大值分別為1.58,2.56,2.02和1.60km,以及0.15,0.24,0.18和0.14m/s。

    在數(shù)值解法下,第一斜壓羅斯貝變形半徑與第一斜壓重力波相速度的值均比在 WKB近似下的值大。最大值出現(xiàn)在北戴河附近海域,為4.72km和0.44m/s;次最大值出現(xiàn)在長興島附近海域,為4.51km和0.42m/s;在萊州灣海域斜壓羅斯貝變形半徑和斜壓重力波相速度的最大值較小,分別為1.76km和0.16 m/s。在遼東灣海域、天津附近海域和黃河口三角洲海域的最大值分別為2.97,2.37,2.55km和0.28,0.22,0.23m/s。兩種算法得出的第一斜壓羅斯貝變形半徑的比較見圖5。

    圖5 兩種算法得出的第一斜壓羅斯貝變形半徑的比較Fig.5 Comparison between the WKB approximation and numerical estimate R1 of the first baroclinic Rossby deformation radius

    從圖5可以看出,渤海夏季第一斜壓羅斯貝變形半徑在WKB近似下的結(jié)果總體低于數(shù)值解法求出的結(jié)果,同樣的情況也出現(xiàn)在其它海區(qū)的相關(guān)研究中。Fennel等[7]和Osinski等[8]在波羅的海的斜壓羅斯貝變形半徑的計算中均表明,WKB近似的結(jié)果偏小于數(shù)值解,F(xiàn)ennel指出偏小10%~30%[7]。這主要是因?yàn)椋琖KB近似適用于浮性頻率垂向變化較弱的情況。由式(9)和式(10)可以看出,在 WKB近似中,第一斜壓重力波的相速度和第一斜壓羅斯貝變形半徑依賴于浮性頻率隨水深的積分,而在實(shí)際求解的過程中,將連續(xù)的海洋分為若干層處理,浮性頻率一般取所在層的平均值,忽略了浮性頻率的變化。故當(dāng)浮性頻率垂向變化較大,WKB近似計算誤差較大。但在深海大洋中計算羅斯貝變形半徑,WKB近似的結(jié)果偏大[6],在此不作詳細(xì)討論。

    圖3b和圖4b中標(biāo)出了斜壓羅斯貝變形半徑出現(xiàn)極值的區(qū)域,其中在北戴河和長興島附近海域的極大值區(qū)有13個站點(diǎn),在渤海中部以西海域的極大值區(qū)有7個站點(diǎn),而在渤海中部以東海域的極小值區(qū)有9個站點(diǎn)。根據(jù)這些站點(diǎn)的溫鹽數(shù)據(jù),計算了極值區(qū)域的浮性頻率的垂向分布,如圖6所示。

    圖6 渤海三極值區(qū)浮性頻率垂向分布圖Fig.6 Profile of Brunt-V?is?l?frequency in the three extreme areas

    由圖6可見,在北戴河和長興島附近海域,以及渤海中部以西海域,浮性頻率的垂向分布較大,導(dǎo)致該海域較強(qiáng)的斜壓性,因此計算得出的第一斜壓羅斯貝變形半徑在該海域出現(xiàn)極大值。而在渤海中部以東海域,浮性頻率的垂向分布整體較小,計算得出的第一斜壓羅斯貝變形半徑在該海域出現(xiàn)極小值。

    圖7繪制了整個渤海夏季浮性頻率的垂向分布??梢娫谒?dbar附近浮性頻率達(dá)到最大,達(dá)0.13s-1,這一深度也是渤海夏季躍層發(fā)生的深度。隨著深度變大,浮性頻率逐漸減小,到水深20dbar的附近,浮性頻率降至約0.04s-1,因此可以看出,渤海夏季浮性頻率垂向上變化較大,故在計算第一斜壓羅斯貝變形半徑的過程中,WKB近似的誤差較大,建議采用斜壓羅斯貝變形半徑的數(shù)值解。

    圖7 整個渤海夏季浮性頻率垂向分布圖Fig.7 Vertical distribution of Brunt-V?is?l?frequency in the whole Bohai Sea

    4 結(jié)論

    本文利用渤海2006年夏季的溫度和鹽度大面觀測資料,運(yùn)用WKB近似和數(shù)值方法計算了渤海夏季第一斜壓羅斯貝變形半徑和第一斜壓重力波相速度。結(jié)果顯示,兩種不同算法得出的斜壓羅斯貝變形半徑和斜壓重力波相速度的分布格局相似,與渤海等深線分布大體一致。數(shù)值解法下,在北戴河和長興島附近海域,斜壓羅斯貝變形半徑較大,最大為4.72km,對應(yīng)的斜壓重力波相速度最大為0.44m/s;在萊州灣海域,斜壓羅斯貝變形半徑較小,最大僅為1.76km,對應(yīng)的斜壓重力波相速度最大為0.16m/s;在渤海中部斜壓羅斯貝變形半徑和斜壓重力波相速度的等值線發(fā)生明顯彎曲,在靠近渤海海峽內(nèi)側(cè)出現(xiàn)一個小值,而與之對應(yīng)在渤海中部以西出現(xiàn)一個大值,導(dǎo)致等值線呈NE-SW向傾斜的“S”狀。在遼東灣海域、天津附近海域和黃河口三角洲海域,斜壓羅斯貝變形半徑和斜壓重力波相速度的最大值分別為2.97,2.37和2.55km,以及0.28,0.22和0.23m/s。

    渤海夏季浮性頻率的垂向變化較大,在水深8dbar附近浮性頻率達(dá)到最大。在WKB近似計算渤海夏季第一斜壓羅斯貝變形半徑的過程中,因其依賴于浮性頻率對水深的積分,得出的結(jié)果誤差較大,其值低于數(shù)值解法求出的結(jié)果。

    (References):

    [1]XIONG X J.The circulation structure and mechanism studies on the China Seas[D].Qingdao:Ocean University of China,2013.熊學(xué)軍.中國近海環(huán)流及其發(fā)生機(jī)制研究[D].青島:中國海洋大學(xué),2013.

    [2]GILL A E.Atmosphere-Ocean dynamics[M].New York:Academic Press,1982.

    [3]ALENIUS P,NEKRASOV A,MYRBERG K.Variability of the baroclinic Rossby radius in the Gulf of Finland[J].Continental Shelf Research,2003,23(6):563-573.

    [4]EMERY W J,LEE W G,MAGAARD L.Geographic and deasonal distributions of Brunt-V?is?l?frequency and Rossby radii in the North Pacific and North Atlantic[J].Journal of Physical Oceanography,1984,14(2):294-317.

    [5]HOURY S,DOMBROWSKY E,DEMEY P,et al.Brunt-V?is?l?frequency and Rossby Radii in the South Atlantic[J].Journal of Physical Oceanography,1987,17(10):1619-1626.

    [6]CHELTON D B,DESZOEKE R A,SCHLAX M G.Geographical variability of the first baroclinic Rossby radius of deformation[J].Journal of Physical Oceanography,1998,28(3):433-460.

    [7]FENNEL W,SEIFERT T,KAYSER B.Rossby radii and phase speeds in the Baltic Sea[J].Continental Shelf Research,1991,11(1):23-36.

    [8]OSINSKI R,RAK D,WALCZOWSKI W,et al.Baroclinic Rossby radius of deformation in the southern Baltic Sea[J].Oceanologia,2010,52(3):417-429.

    [9]GAN Z J,CAI S Q.Geographical and seasonal variability of Rossby radii in South China Sea[J].Journal of Tropical Oceanography,2001,20(1):1-9.甘子鈞,蔡樹群.南海羅斯貝變形半徑的地理及季節(jié)變化[J].熱帶海洋學(xué)報,2001,20(1):1-9.

    [10]CAI S Q,LONG X M,WU R H,et al.Geographical and monthly variability of the first baroclinic Rossby radius of deformation in the South China Sea[J].Journal of Marine Systems,2008,74(1-2):711-720.

    [11]YU Z H,YANG D S,HE H Y,et al.Geophysical fluid dynamics[M].Beijing:China Meteorological Press,1996.余志豪,楊大升,賀海晏,等.地球物理流體動力學(xué)[M].北京:氣象出版社,1996.

    猜你喜歡
    斜壓浮性渤海
    渤海大學(xué)作品精選
    聲屏世界(2022年17期)2022-12-18 01:42:44
    淺析船舶壓載水處理技術(shù)及處理系統(tǒng)選型
    倒扣駁船扳正過程的仿真模擬
    航海(2020年5期)2020-11-02 09:35:16
    Gemini表面活性劑體系下鈣質(zhì)磷礦中白云石的可浮性研究
    “烈焰”,在渤海邊燃燒
    斜壓渦度的變化與臺風(fēng)暴雨的關(guān)系研究
    基于Visual FoxPro的船舶靜水力計算及曲線繪制
    價值工程(2017年3期)2017-05-14 02:11:40
    遼代“斡魯朵”內(nèi)的渤海人移民
    Rossby波的線性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運(yùn)的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    久久久国产成人免费| 色哟哟·www| 美女黄网站色视频| 丰满人妻一区二区三区视频av| 91久久精品国产一区二区三区| 国产视频首页在线观看| 嫩草影院入口| 青春草国产在线视频 | 亚洲成人精品中文字幕电影| 日韩精品有码人妻一区| 国产视频首页在线观看| 一级二级三级毛片免费看| 又粗又爽又猛毛片免费看| 自拍偷自拍亚洲精品老妇| 人妻系列 视频| av免费观看日本| 啦啦啦啦在线视频资源| 久久婷婷人人爽人人干人人爱| 九草在线视频观看| 免费一级毛片在线播放高清视频| 最近手机中文字幕大全| 国产精品日韩av在线免费观看| 婷婷六月久久综合丁香| 此物有八面人人有两片| 免费观看精品视频网站| 日韩一区二区三区影片| 岛国在线免费视频观看| 国产成人福利小说| 国内揄拍国产精品人妻在线| 国产精品综合久久久久久久免费| 精品不卡国产一区二区三区| 午夜久久久久精精品| 最近2019中文字幕mv第一页| 精品99又大又爽又粗少妇毛片| 尤物成人国产欧美一区二区三区| 男女那种视频在线观看| 午夜视频国产福利| 久久精品国产亚洲av涩爱 | 国产亚洲5aaaaa淫片| 又粗又爽又猛毛片免费看| 欧美日韩精品成人综合77777| 中文精品一卡2卡3卡4更新| 成人国产麻豆网| 波多野结衣高清无吗| 99国产极品粉嫩在线观看| 成人午夜高清在线视频| avwww免费| 成人国产麻豆网| 国产精品久久久久久精品电影| 久久久久久九九精品二区国产| 亚洲精品乱码久久久久久按摩| 白带黄色成豆腐渣| 亚洲中文字幕一区二区三区有码在线看| 91aial.com中文字幕在线观看| 亚洲精品日韩在线中文字幕 | 国产在线精品亚洲第一网站| ponron亚洲| 亚洲成a人片在线一区二区| а√天堂www在线а√下载| 免费电影在线观看免费观看| 18禁黄网站禁片免费观看直播| 成人毛片a级毛片在线播放| 一级黄色大片毛片| 伊人久久精品亚洲午夜| 最近最新中文字幕大全电影3| 一级av片app| 最近的中文字幕免费完整| 插阴视频在线观看视频| 色吧在线观看| 狠狠狠狠99中文字幕| 亚洲精品久久国产高清桃花| 人妻系列 视频| 国产精品永久免费网站| 成人永久免费在线观看视频| 99视频精品全部免费 在线| 亚洲精品日韩在线中文字幕 | 日韩强制内射视频| 日本av手机在线免费观看| 精品一区二区三区视频在线| 老师上课跳d突然被开到最大视频| 99热这里只有是精品50| 欧美日韩综合久久久久久| 亚洲国产精品sss在线观看| 亚洲一级一片aⅴ在线观看| 18禁在线播放成人免费| 秋霞在线观看毛片| 国产精品电影一区二区三区| 国产成人精品久久久久久| 亚洲乱码一区二区免费版| 97超视频在线观看视频| 亚洲精品影视一区二区三区av| 一本久久中文字幕| 别揉我奶头 嗯啊视频| 97超碰精品成人国产| av福利片在线观看| 久久久a久久爽久久v久久| 欧洲精品卡2卡3卡4卡5卡区| 日韩国内少妇激情av| 国产精品综合久久久久久久免费| 99久久精品热视频| 日韩三级伦理在线观看| 中文精品一卡2卡3卡4更新| 91久久精品国产一区二区成人| 丝袜美腿在线中文| 极品教师在线视频| 精品久久久久久久人妻蜜臀av| .国产精品久久| 美女脱内裤让男人舔精品视频 | 三级男女做爰猛烈吃奶摸视频| 日韩精品青青久久久久久| 国产亚洲av嫩草精品影院| 国产成人aa在线观看| 高清毛片免费看| 亚洲av中文字字幕乱码综合| 在线观看免费视频日本深夜| 亚洲综合色惰| 中国国产av一级| 亚洲va在线va天堂va国产| 天天躁日日操中文字幕| 国产精品伦人一区二区| 久久精品国产亚洲av涩爱 | .国产精品久久| 成人一区二区视频在线观看| 国内少妇人妻偷人精品xxx网站| 久久精品综合一区二区三区| 欧美激情在线99| 欧美色欧美亚洲另类二区| 国产成人a区在线观看| 亚洲在线观看片| 一个人看视频在线观看www免费| 亚洲中文字幕一区二区三区有码在线看| 五月伊人婷婷丁香| 久久久久久国产a免费观看| or卡值多少钱| 99热只有精品国产| 国产精品久久久久久av不卡| 亚洲欧美中文字幕日韩二区| 日本熟妇午夜| 国产精品,欧美在线| 99久国产av精品国产电影| 桃色一区二区三区在线观看| 国产成人freesex在线| 日韩欧美一区二区三区在线观看| 乱系列少妇在线播放| 欧美日本视频| 99热网站在线观看| 午夜激情欧美在线| 变态另类成人亚洲欧美熟女| 成年女人永久免费观看视频| 亚洲精品456在线播放app| 亚洲av二区三区四区| 国产精品一区二区三区四区久久| 亚洲一区高清亚洲精品| 国产精品电影一区二区三区| av.在线天堂| 亚洲一区二区三区色噜噜| 久久久久久久久大av| 搡老妇女老女人老熟妇| 国产单亲对白刺激| 国产亚洲精品av在线| 中文字幕av成人在线电影| 亚洲av男天堂| 成年女人永久免费观看视频| 精品欧美国产一区二区三| 色视频www国产| 国产黄片视频在线免费观看| 给我免费播放毛片高清在线观看| 久久精品国产亚洲av天美| 国产精品野战在线观看| 亚洲人成网站高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 黄色配什么色好看| 一级av片app| av女优亚洲男人天堂| 国产黄片视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 日日摸夜夜添夜夜添av毛片| 成人高潮视频无遮挡免费网站| 欧美潮喷喷水| 天堂av国产一区二区熟女人妻| 久久久久国产网址| 国产熟女欧美一区二区| 日韩制服骚丝袜av| 精品久久久久久久末码| 国产在视频线在精品| 国产毛片a区久久久久| 国产三级在线视频| 欧美性感艳星| 一进一出抽搐动态| 久久久a久久爽久久v久久| 国语自产精品视频在线第100页| 中文欧美无线码| avwww免费| 国产高清有码在线观看视频| 91精品国产九色| 国产精品不卡视频一区二区| 男人舔女人下体高潮全视频| 成年女人永久免费观看视频| 久久精品久久久久久久性| 看黄色毛片网站| 国产三级中文精品| 美女xxoo啪啪120秒动态图| 午夜爱爱视频在线播放| 国产精品野战在线观看| 此物有八面人人有两片| 亚洲av一区综合| 麻豆乱淫一区二区| 精品一区二区三区视频在线| 99热精品在线国产| 日韩精品青青久久久久久| 国产成人91sexporn| 欧美性猛交黑人性爽| 精品欧美国产一区二区三| 久久久精品大字幕| 国产女主播在线喷水免费视频网站 | 别揉我奶头 嗯啊视频| a级毛片a级免费在线| 欧美成人a在线观看| 校园春色视频在线观看| 白带黄色成豆腐渣| 欧美精品国产亚洲| 国产麻豆成人av免费视频| 小蜜桃在线观看免费完整版高清| 国产亚洲5aaaaa淫片| 久久综合国产亚洲精品| 高清毛片免费看| 在线a可以看的网站| 国产 一区精品| 日韩强制内射视频| 一级黄片播放器| 国产精品人妻久久久久久| 久久亚洲国产成人精品v| 成人亚洲欧美一区二区av| 99久久中文字幕三级久久日本| 91久久精品国产一区二区成人| 国产久久久一区二区三区| 在线观看66精品国产| 国内揄拍国产精品人妻在线| 91麻豆精品激情在线观看国产| 亚州av有码| 18禁黄网站禁片免费观看直播| 国产一区二区在线观看日韩| 天堂中文最新版在线下载 | 精品国产三级普通话版| av免费在线看不卡| 一夜夜www| 欧美激情在线99| 亚洲va在线va天堂va国产| 麻豆成人午夜福利视频| 亚洲av不卡在线观看| 国产亚洲av片在线观看秒播厂 | 99久久精品国产国产毛片| 午夜精品在线福利| 啦啦啦观看免费观看视频高清| 欧美三级亚洲精品| 搡女人真爽免费视频火全软件| 精品久久久久久久末码| 丝袜喷水一区| 亚洲精品自拍成人| av在线天堂中文字幕| 国产v大片淫在线免费观看| 欧美激情久久久久久爽电影| 亚洲美女视频黄频| 久久久久九九精品影院| 久久久精品大字幕| 又粗又爽又猛毛片免费看| 偷拍熟女少妇极品色| 国产精品伦人一区二区| 在线天堂最新版资源| 国产国拍精品亚洲av在线观看| 精品人妻熟女av久视频| 熟妇人妻久久中文字幕3abv| 3wmmmm亚洲av在线观看| 能在线免费看毛片的网站| 老女人水多毛片| 大香蕉久久网| 99热只有精品国产| 国产大屁股一区二区在线视频| av国产免费在线观看| 国产国拍精品亚洲av在线观看| 丝袜美腿在线中文| 日韩,欧美,国产一区二区三区 | av在线天堂中文字幕| 亚洲三级黄色毛片| 亚洲一级一片aⅴ在线观看| 午夜福利成人在线免费观看| av免费观看日本| 美女 人体艺术 gogo| 亚洲欧美精品自产自拍| 69av精品久久久久久| 久久久久免费精品人妻一区二区| 国产一区二区亚洲精品在线观看| 国产av一区在线观看免费| 亚洲乱码一区二区免费版| 久99久视频精品免费| 婷婷精品国产亚洲av| 中国国产av一级| 真实男女啪啪啪动态图| 国产一区二区在线观看日韩| 99久国产av精品| 91久久精品电影网| 亚洲国产精品成人久久小说 | 国产人妻一区二区三区在| 91狼人影院| 一边亲一边摸免费视频| 国产久久久一区二区三区| 久久国产乱子免费精品| 夜夜夜夜夜久久久久| 两个人的视频大全免费| 午夜精品在线福利| 亚洲人成网站在线播放欧美日韩| 97热精品久久久久久| 搞女人的毛片| 精品久久国产蜜桃| 国内精品久久久久精免费| 久久久久久伊人网av| 成人鲁丝片一二三区免费| 波多野结衣高清无吗| av女优亚洲男人天堂| 免费观看在线日韩| 一本久久中文字幕| 日本黄色视频三级网站网址| 国产精华一区二区三区| 在线观看免费视频日本深夜| 美女内射精品一级片tv| 国产伦理片在线播放av一区 | 亚洲欧美精品综合久久99| 自拍偷自拍亚洲精品老妇| 99久国产av精品国产电影| 国产精品一区二区在线观看99 | 高清午夜精品一区二区三区 | 欧美一区二区精品小视频在线| 成人无遮挡网站| 少妇熟女aⅴ在线视频| 久久久精品欧美日韩精品| 亚洲av中文字字幕乱码综合| 久久久久久久午夜电影| 美女被艹到高潮喷水动态| 中文精品一卡2卡3卡4更新| 亚洲国产高清在线一区二区三| 欧美在线一区亚洲| 欧美激情久久久久久爽电影| 亚洲欧美日韩卡通动漫| 国产黄片美女视频| 夜夜爽天天搞| 久久久久久久久久久免费av| 两个人视频免费观看高清| 久久久久免费精品人妻一区二区| 我要搜黄色片| 亚洲欧美日韩东京热| 黄色视频,在线免费观看| 精品久久久久久久人妻蜜臀av| 国产一区二区在线av高清观看| 色尼玛亚洲综合影院| 国产一区二区在线观看日韩| 国产精品久久久久久久久免| 日本色播在线视频| 嘟嘟电影网在线观看| 亚洲欧美日韩无卡精品| 一级av片app| 亚洲国产精品sss在线观看| 热99re8久久精品国产| 欧美一区二区精品小视频在线| 精品一区二区三区视频在线| 日韩人妻高清精品专区| 久久精品夜色国产| 国产精品久久久久久亚洲av鲁大| 秋霞在线观看毛片| 亚洲自拍偷在线| 久久精品夜色国产| 2021天堂中文幕一二区在线观| 最近的中文字幕免费完整| 男女啪啪激烈高潮av片| 夜夜夜夜夜久久久久| 村上凉子中文字幕在线| 男人狂女人下面高潮的视频| 少妇丰满av| 特大巨黑吊av在线直播| 最新中文字幕久久久久| 日韩欧美国产在线观看| 国产精品免费一区二区三区在线| 成人综合一区亚洲| 色综合站精品国产| 亚洲精华国产精华液的使用体验 | 黄色配什么色好看| 好男人视频免费观看在线| 看十八女毛片水多多多| 欧美xxxx黑人xx丫x性爽| 人人妻人人澡人人爽人人夜夜 | 亚洲精品乱码久久久久久按摩| 免费观看的影片在线观看| 少妇熟女aⅴ在线视频| 丰满的人妻完整版| 中出人妻视频一区二区| 久久人人精品亚洲av| 国产麻豆成人av免费视频| 国产高潮美女av| 久久国内精品自在自线图片| 婷婷精品国产亚洲av| 少妇裸体淫交视频免费看高清| 亚洲va在线va天堂va国产| 久久久久久久亚洲中文字幕| 永久网站在线| 免费大片18禁| 直男gayav资源| 免费看光身美女| 久久人人爽人人爽人人片va| 少妇熟女欧美另类| 精品久久久久久久人妻蜜臀av| 国产精品一及| 国产精品一二三区在线看| 欧美三级亚洲精品| av在线播放精品| 欧美高清成人免费视频www| 黑人高潮一二区| 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| 亚洲av免费高清在线观看| 看十八女毛片水多多多| 深爱激情五月婷婷| 可以在线观看毛片的网站| 日日撸夜夜添| 欧美成人精品欧美一级黄| 12—13女人毛片做爰片一| 亚洲五月天丁香| av在线亚洲专区| 超碰av人人做人人爽久久| 亚洲,欧美,日韩| 日日摸夜夜添夜夜爱| 网址你懂的国产日韩在线| 高清毛片免费观看视频网站| 悠悠久久av| 一级黄片播放器| 久久国产乱子免费精品| 国产av一区在线观看免费| av又黄又爽大尺度在线免费看 | 久久精品国产亚洲av涩爱 | 国产高潮美女av| 秋霞在线观看毛片| 国产一区二区在线观看日韩| 久久久a久久爽久久v久久| 国产av不卡久久| 亚洲av中文字字幕乱码综合| 亚洲av第一区精品v没综合| 人人妻人人澡人人爽人人夜夜 | АⅤ资源中文在线天堂| 在线播放无遮挡| 日韩成人av中文字幕在线观看| 九草在线视频观看| 天堂影院成人在线观看| 波多野结衣高清作品| 日韩大尺度精品在线看网址| 国产视频首页在线观看| www.av在线官网国产| 观看美女的网站| 婷婷精品国产亚洲av| 亚洲四区av| 欧美一区二区精品小视频在线| 99热这里只有是精品在线观看| 此物有八面人人有两片| 久久精品久久久久久久性| 日韩精品有码人妻一区| 黄色欧美视频在线观看| 久久久精品欧美日韩精品| 国内精品一区二区在线观看| 亚洲精品影视一区二区三区av| 国产极品精品免费视频能看的| 天天躁日日操中文字幕| 黄色一级大片看看| 波多野结衣高清作品| 国产男人的电影天堂91| 男女那种视频在线观看| 国产成年人精品一区二区| 在线播放国产精品三级| 久久精品夜夜夜夜夜久久蜜豆| 又爽又黄无遮挡网站| 边亲边吃奶的免费视频| 春色校园在线视频观看| 国产乱人偷精品视频| 97人妻精品一区二区三区麻豆| 久久久久性生活片| 欧美最新免费一区二区三区| 国产精品综合久久久久久久免费| 一区二区三区四区激情视频 | av国产免费在线观看| 午夜激情福利司机影院| 国产一级毛片七仙女欲春2| 一边摸一边抽搐一进一小说| 国语自产精品视频在线第100页| 女人十人毛片免费观看3o分钟| 伊人久久精品亚洲午夜| 日本免费a在线| 在线播放无遮挡| or卡值多少钱| 在线播放无遮挡| 激情 狠狠 欧美| 午夜视频国产福利| 狂野欧美激情性xxxx在线观看| 看片在线看免费视频| 少妇的逼好多水| 三级经典国产精品| 日本爱情动作片www.在线观看| 欧美+日韩+精品| 综合色av麻豆| 18+在线观看网站| 亚洲精品国产成人久久av| 美女被艹到高潮喷水动态| АⅤ资源中文在线天堂| 国产日韩欧美在线精品| av又黄又爽大尺度在线免费看 | 中文字幕免费在线视频6| 男女下面进入的视频免费午夜| 亚洲第一电影网av| 九九久久精品国产亚洲av麻豆| 欧美潮喷喷水| 亚洲国产精品成人综合色| 亚洲欧美精品自产自拍| 国产精品伦人一区二区| 五月玫瑰六月丁香| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲最大成人手机在线| 久久国产乱子免费精品| 精品久久久噜噜| 校园人妻丝袜中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 日韩欧美国产在线观看| 舔av片在线| 精品不卡国产一区二区三区| 欧美色欧美亚洲另类二区| 丰满的人妻完整版| 91精品一卡2卡3卡4卡| 在线观看午夜福利视频| 最好的美女福利视频网| 国产高清三级在线| 一个人免费在线观看电影| 天堂中文最新版在线下载 | 亚洲av中文av极速乱| 精品久久久久久成人av| 久久久久久久午夜电影| 夫妻性生交免费视频一级片| 99热全是精品| 久久中文看片网| 亚洲,欧美,日韩| 亚洲最大成人手机在线| 18禁黄网站禁片免费观看直播| 日韩大尺度精品在线看网址| 女人十人毛片免费观看3o分钟| 一进一出抽搐gif免费好疼| 如何舔出高潮| 国产一区二区在线av高清观看| 老熟妇乱子伦视频在线观看| 最近的中文字幕免费完整| 久久久成人免费电影| 久久久精品大字幕| 九九热线精品视视频播放| 亚洲成人av在线免费| 能在线免费观看的黄片| 久久精品国产亚洲av香蕉五月| 国产激情偷乱视频一区二区| 免费搜索国产男女视频| 国产高潮美女av| 99热这里只有是精品在线观看| 婷婷色综合大香蕉| 日本一本二区三区精品| 天美传媒精品一区二区| 三级国产精品欧美在线观看| 亚洲久久久久久中文字幕| 男女边吃奶边做爰视频| 国产黄片美女视频| 日本撒尿小便嘘嘘汇集6| 99热只有精品国产| 久久鲁丝午夜福利片| 麻豆国产av国片精品| 亚洲欧美成人综合另类久久久 | 蜜臀久久99精品久久宅男| 精品人妻一区二区三区麻豆| 在线观看美女被高潮喷水网站| 亚洲第一区二区三区不卡| 精品日产1卡2卡| 久久久精品94久久精品| 一级黄色大片毛片| 干丝袜人妻中文字幕| 少妇的逼好多水| 夫妻性生交免费视频一级片| 啦啦啦啦在线视频资源| 国产av不卡久久| 国产精品1区2区在线观看.| 国产视频内射| 亚洲av免费在线观看| 国产成人freesex在线| 国产黄色视频一区二区在线观看 | 毛片女人毛片| 蜜臀久久99精品久久宅男| 午夜福利成人在线免费观看| 国产精品久久久久久久电影| 日韩中字成人| 亚洲不卡免费看| 国产av不卡久久| 伊人久久精品亚洲午夜| 国产亚洲av嫩草精品影院| 久久欧美精品欧美久久欧美| 亚洲av不卡在线观看| 亚洲av.av天堂| 丝袜喷水一区| 长腿黑丝高跟| 色5月婷婷丁香| 少妇猛男粗大的猛烈进出视频 | 自拍偷自拍亚洲精品老妇| 国产单亲对白刺激| 一区福利在线观看| 别揉我奶头 嗯啊视频| 久久人人爽人人爽人人片va| 国产亚洲精品av在线| 亚洲高清免费不卡视频| 日本免费一区二区三区高清不卡|