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

    海洋運(yùn)動(dòng)對(duì)大尺度風(fēng)場(chǎng)響應(yīng)的解析研究*

    2012-09-06 08:22:40趙艷玲劉天軍
    海洋科學(xué)進(jìn)展 2012年3期
    關(guān)鍵詞:壓力梯度時(shí)間尺度風(fēng)場(chǎng)

    趙艷玲,鄧 冰,張 銘,劉天軍

    (1.中國(guó)人民解放軍61741部隊(duì),北京100094;2.北京應(yīng)用氣象研究所,北京100029;3.中國(guó)人民解放軍理工大學(xué)氣象學(xué)院 大氣環(huán)流與短期氣候預(yù)測(cè)實(shí)驗(yàn)室,江蘇 南京,211101)

    海洋運(yùn)動(dòng)對(duì)大尺度風(fēng)場(chǎng)響應(yīng)的解析研究*

    趙艷玲1,3,鄧 冰2,3,張 銘3,劉天軍2

    (1.中國(guó)人民解放軍61741部隊(duì),北京100094;2.北京應(yīng)用氣象研究所,北京100029;3.中國(guó)人民解放軍理工大學(xué)氣象學(xué)院 大氣環(huán)流與短期氣候預(yù)測(cè)實(shí)驗(yàn)室,江蘇 南京,211101)

    為探討海洋向大尺度風(fēng)應(yīng)力響應(yīng)的具體過程,建立了考慮海面風(fēng)應(yīng)力強(qiáng)迫的正壓原始方程海洋模型,并利用Green函數(shù)做了解析求解和討論。發(fā)現(xiàn)風(fēng)應(yīng)力引起的Ekman抽吸正比于風(fēng)應(yīng)力旋度。得到了海洋向風(fēng)應(yīng)力響應(yīng)的具體形式:即海洋流場(chǎng)和壓力場(chǎng)進(jìn)行著準(zhǔn)平衡的演變,并有海洋中壓力梯度與科氏力之間的準(zhǔn)地轉(zhuǎn)平衡以及海洋Ekman體積輸送的散度場(chǎng)與風(fēng)應(yīng)力旋度場(chǎng)之間的準(zhǔn)平衡;此時(shí)有短周期的重力慣性波被激發(fā),該波會(huì)較快彌散掉;而準(zhǔn)平衡的演變則始終持續(xù)著。海洋在大尺度風(fēng)場(chǎng)的強(qiáng)迫下的運(yùn)動(dòng)由重力慣性波、強(qiáng)迫慣性振蕩和緩慢演變的準(zhǔn)平衡過程三部分組成;前兩部分是快過程,第一部分更快,第三部分是慢過程。最終海洋運(yùn)動(dòng)表現(xiàn)為強(qiáng)迫慣性振蕩和準(zhǔn)平衡演變過程的疊加。

    大尺度風(fēng)場(chǎng);海洋模型;響應(yīng);Ekman抽吸

    風(fēng)應(yīng)力驅(qū)動(dòng)上層海洋的理論早在20世紀(jì)50年代就已出現(xiàn)[1],至90年代以來,海洋運(yùn)動(dòng)向風(fēng)場(chǎng)的響應(yīng)問題受到越來越多的關(guān)注。這種響應(yīng)是海氣相互作用的重要組成部分,在全球氣候變化中起到了重要作用。直到目前,這方面的研究仍受到廣泛重視。錢永甫等[2]利用POM模式(分辨率0.5°)研究了南海海流對(duì)冬季風(fēng)應(yīng)力的響應(yīng)特征,認(rèn)為南海淺表層海流和全流對(duì)定常風(fēng)應(yīng)力的響應(yīng)時(shí)間均為120d,深層海流的響應(yīng)時(shí)間較長(zhǎng)。蔡樹群等[3]利用一個(gè)分區(qū)性的正壓、斜壓銜接模式,考察了南海上層環(huán)流對(duì)季風(fēng)轉(zhuǎn)變的響應(yīng),發(fā)現(xiàn)在不同方式變化的季風(fēng)轉(zhuǎn)變、季風(fēng)過渡時(shí)期,南海流場(chǎng)的調(diào)整有較大的差別,但在過渡之后最終的流場(chǎng)結(jié)構(gòu)基本一致。對(duì)于渤、黃東海區(qū)域海流對(duì)冬季風(fēng)應(yīng)力的響應(yīng)研究,趙保仁[4]利用黃海和東海近岸3個(gè)測(cè)站的觀測(cè)資料作了初步探討,認(rèn)為上表層有離岸的Ekman輸送時(shí),趨岸的補(bǔ)償流大約從1/3水深處直到海底的較厚水層中向岸流去,表層流大約6~12h可建立與風(fēng)相適應(yīng)的風(fēng)生流動(dòng),底層流的適應(yīng)時(shí)間大約是15~20 h。Qiu[5]則研究了東太平洋風(fēng)應(yīng)力強(qiáng)迫的Rossby波。Cabanes等[6]研究了風(fēng)應(yīng)力驅(qū)動(dòng)的正壓Sverdrup平衡和一階斜壓Rossby波。Qiu和Chen[7]還研究了南太平洋海表面高度時(shí)空變化的動(dòng)力機(jī)制等。以上這些研究以診斷分析和數(shù)值模擬居多。在解析研究方面,Gill在文獻(xiàn)[1]中指出,在一定條件下,此時(shí)海洋運(yùn)動(dòng)可分解為兩部分:一部分是風(fēng)應(yīng)力作為直接強(qiáng)迫項(xiàng)驅(qū)動(dòng)的強(qiáng)迫慣性振蕩,在該部分壓力梯度不起作用;另一部分是風(fēng)應(yīng)力引起的Ekman抽吸,此時(shí)則存在該抽吸引起的垂直運(yùn)動(dòng),并通過其進(jìn)一步影響海洋運(yùn)動(dòng),該部分有壓力梯度。然而總的看來,這方面的研究尚不很多。本研究曾建立了一個(gè)水平二維正壓準(zhǔn)平衡海洋模型,并在考慮西海岸的情況下,對(duì)氣候尺度的準(zhǔn)定常和時(shí)變風(fēng)場(chǎng)(異常)強(qiáng)迫下的海洋流場(chǎng)響應(yīng)做了解析求解[8],并試圖解釋中緯度北太平洋實(shí)際流場(chǎng)異常的原因。因該模型做了準(zhǔn)平衡(準(zhǔn)無輻散)近似,重力慣性波被濾去了,故其僅能反映海洋對(duì)大尺度風(fēng)應(yīng)力響應(yīng)的最終狀態(tài),而不能給出其響應(yīng)的具體過程,這點(diǎn)為其局限所在。為此,本研究建立了一個(gè)考慮海面風(fēng)應(yīng)力強(qiáng)迫的線性化正壓原始方程海洋模型,在對(duì)該模型做了理論分析后,進(jìn)一步對(duì)該模型做了解析求解和討論,以便探討正壓海洋對(duì)大尺度風(fēng)應(yīng)力響應(yīng)的具體過程。

    1 數(shù)學(xué)模型

    海洋對(duì)風(fēng)應(yīng)力的線性化正壓響應(yīng)方程組[1]為:

    式中,u為x方向流速,v為y方向流速;Xs,Xb分別表示海面風(fēng)應(yīng)力和海底摩擦應(yīng)力的x方向分量,Ys,Yb則為y方向分量;η為受擾后流體厚度與未受擾流體厚度的偏差,H為未受擾流體厚度,ρ為流體密度,兩者均取為常數(shù)。令

    式中,UE,VE是海面與海底邊界層的Ekman體積輸送,其滿足方程

    而up,vp則滿足以下方程

    在此φ=gη,其為位勢(shì)偏差。在式(4)中,強(qiáng)迫項(xiàng)已被消去。將式(2)代入式(1)的連續(xù)方程中,再乘以g,則可得

    這里co2=gH,wE=?UE/?x+?VE/?y,其為由風(fēng)應(yīng)力造成的Ekman抽吸的垂直運(yùn)動(dòng)。當(dāng)海洋表面出現(xiàn)體積輻散時(shí),其下部的水要向上涌升,故有wE>0,反之亦然。

    由此可見,風(fēng)應(yīng)力造成的海洋運(yùn)動(dòng)能夠分解成兩部分。一部分是風(fēng)應(yīng)力作為直接強(qiáng)迫項(xiàng)驅(qū)動(dòng)的強(qiáng)迫慣性振蕩,其控制方程組由式(3)式?jīng)Q定,注意在該部分壓力梯度不起作用。另一部分是風(fēng)應(yīng)力引起的Ekman抽吸,此時(shí)存在Ekman抽吸引起的垂直運(yùn)動(dòng),并通過該垂直運(yùn)動(dòng)進(jìn)一步影響海洋運(yùn)動(dòng);其控制方程為式(4)和式(5),該部分有壓力梯度。

    對(duì)于第一部分即海洋的強(qiáng)迫慣性振蕩部分,文獻(xiàn)[1]已給出了直角坐標(biāo)系下的解,還討論了在極坐標(biāo)下熱帶氣旋造成的海洋強(qiáng)迫振蕩問題[9],故本研究就不再討論該問題了。下面重點(diǎn)對(duì)第二部分進(jìn)行討論。在討論之前,先來估計(jì)wE與風(fēng)應(yīng)力(Xs,Ys)之間的關(guān)系。

    若不考慮海底摩擦應(yīng)力,由式(3)可得,

    將以上UE、VE代入wE=?UE/?x+?VE/?y式中,則有,

    若假定海洋運(yùn)動(dòng)呈準(zhǔn)定?;蛳鄬?duì)慣性時(shí)間尺度f-1變化很緩慢,即其屬海洋大尺度運(yùn)動(dòng),此時(shí)在熱帶外海洋其運(yùn)動(dòng)的時(shí)間尺度大于1d,空間尺度大于100km(氣候問題中海洋運(yùn)動(dòng)顯然屬該大尺度范疇);本文章中的“海洋”均為以上所指的大尺度海洋,在此不討論海洋中小尺度的問題,以下不再贅述。這樣式(7)最后1個(gè)等號(hào)右邊的第二項(xiàng)相對(duì)于第一項(xiàng)是小項(xiàng),其可略,故此時(shí)有:

    即Ekman抽吸正比于風(fēng)應(yīng)力的旋度,且此時(shí)有海洋Ekman體積輸送的散度場(chǎng)和風(fēng)應(yīng)力旋度場(chǎng)之間的準(zhǔn)平衡,并能得到Ekman抽吸的估計(jì)式(8),用該式則可從風(fēng)應(yīng)力的分布診斷得到Ekman抽吸的垂直運(yùn)動(dòng)。

    2 方程的求解

    現(xiàn)設(shè)海洋運(yùn)動(dòng)屬海洋大尺度,此時(shí)方程(5)可寫為:

    式中,aφ=?YS/?x-?XS/?y為風(fēng)應(yīng)力的旋度。方程組(4)和(9)構(gòu)成了上述海洋運(yùn)動(dòng)第二部分的控制方程組。下面對(duì)該方程組進(jìn)行求解。

    首先用式(4)可分別求得渦度方程和散度方程:

    然后消去渦度,可得,

    這里算子Δ=?2/?x2+?2/?y2,為二維Laplace算子。最后將式(9)與式(11)聯(lián)立以消去散度,則可得一個(gè)與單一變量φ的傾向即?φ/?t有關(guān)的方程:

    在此算子

    其為正壓流體的波動(dòng)算子,而BΩ,Bφ,Bu,Bv則定義如下:

    由上式可見,BΩ,Bφ,Bu,Bv則均與風(fēng)應(yīng)力的旋度有關(guān)。

    對(duì)于全平面的情況,非齊次方程的解可由相應(yīng)齊次方程的解加上非齊次項(xiàng)導(dǎo)致的“推遲位勢(shì)”來表示。對(duì)式(13)~(15)來說,齊次方程組Λ(F)=0是典型的波動(dòng)方程(F 為?up/?t,?vp/?t,?φ/?t),存在有滿足初值F(x,y,0)=F(0)(x,y)和?F(x,y,0)/?t=F(1)(x,y)的唯一解;這樣對(duì)非齊次方程組Λ(F)=-A(x,y,t)而言(這里-A為非齊次項(xiàng),即式(13)~(15)中等號(hào)右端的部分,其均與風(fēng)應(yīng)力的旋度有關(guān)),則借鑒文獻(xiàn)[10]中的數(shù)學(xué)處理,可得到其解為:

    其中

    在此有

    而在式(18)~(20)中r的表達(dá)式則為

    注意,式(18)中的F 實(shí)為傾向,即?φ/?t、?up/?t和?vp/?t。將式(18)對(duì)時(shí)間t積分,此時(shí)要用到分部積分和積分號(hào)下求微商等運(yùn)算,最終可得[10]:

    其中,

    這樣就求出了非齊次方程組式(13)~(15)的解,該解也是控制方程組式(4)和(9)的解。以下將對(duì)該解做出分析討論。

    3 討 論

    以下仍然借鑒文獻(xiàn)[10]中的思路來進(jìn)行分析討論。當(dāng)存在風(fēng)應(yīng)力的外強(qiáng)迫,即Aφ≡0時(shí),則由解式(22)~(24)可見,其可分為兩部分(該兩部分在式(22)~(24)中分別用大括號(hào)括起)。第一部分(下標(biāo)為I的部分)為初始場(chǎng)的影響,其中包括初始風(fēng)應(yīng)力場(chǎng)的影響,但初始風(fēng)應(yīng)力場(chǎng)只直接影響φ(x,y,t),而不直接影響up(x,y,t),vp(x,y,t)(因這兩者無A(0)φ項(xiàng))。第二部分則為外強(qiáng)迫,即風(fēng)應(yīng)力強(qiáng)迫對(duì)解的影響,并給出了大尺度海洋向風(fēng)應(yīng)力響應(yīng)的具體形式。

    以上已指出,第一部分為受初始擾動(dòng)影響的φi,upI,vpI,注意到,若初始擾動(dòng)只局限在局部地區(qū),或有(?φ/?t)(0),(?up/?t)(0),(?vp/?t)(0)以及等均 ∈L2,即在全平面能量有限,則在較長(zhǎng)時(shí)間后,因重力慣性波的彌散,φI,upI,vpI將趨于平衡。這種平衡是海洋中壓力梯度與科氏力之間的地轉(zhuǎn)平衡。若(?φ/?t)(0)等不局限在有限地區(qū),但設(shè)初始時(shí)刻擾動(dòng)是準(zhǔn)地轉(zhuǎn)的(這表明初始風(fēng)應(yīng)力引起的海洋非地轉(zhuǎn)運(yùn)動(dòng)不明顯),則在任何時(shí)刻φI,upI,vpI始終是準(zhǔn)地轉(zhuǎn)的。

    以下分析第二部分。這部分決定外強(qiáng)迫即風(fēng)應(yīng)力對(duì)海洋適應(yīng)過程的影響(如果有適應(yīng)過程的話,由下文可見,對(duì)大尺度風(fēng)應(yīng)力的確有適應(yīng)過程)。從式(22)~(24)可見,BΩ,Bφ,Bu,Bv描述了風(fēng)應(yīng)力激發(fā)出的海洋擾動(dòng),其與 Green函數(shù)0(r,t-t′)和?0(r,t-t′)/?t′密切有關(guān)。這兩種 Green函數(shù)的圖形如圖1、圖2所示。圖1給出了Green函數(shù)隨時(shí)間的變化,其橫坐標(biāo)是無量綱時(shí)間τ=f·t,其中f為地轉(zhuǎn)參數(shù),t為時(shí)間,縱坐標(biāo)則是與Green函數(shù)0和?0/?t′有關(guān)的各種函數(shù)值κ(其中也包含這兩種Green函數(shù)的值),值κ也均為無量綱量,而各種函數(shù)的具體形式則見圖1右方數(shù)字1~5后的表達(dá)式。圖2給出了Green函數(shù)0(r,t-t′)/?t′和?0(r,t-t′)的空間分布,其橫坐標(biāo)為無量綱距離ξ=fr/c0→0)。由圖1可見,當(dāng)t→∞即τ=f·t→∞時(shí)則有~K0(r,t)→(fr/c0)(由該圖可見曲線2逼近直線3,在此K0是虛變量零階Bessel函數(shù)),即麥克唐納函數(shù))0(r,t)圍繞著平衡位置K0(fr/c0)上下擺動(dòng)(參見該圖上的曲線2);但當(dāng)r很小時(shí),即r<<c0/f,也即r遠(yuǎn)小于海洋運(yùn)動(dòng)的Rossby變形半徑c0/f時(shí),則0(r,t)恒為正;當(dāng)r較大時(shí)(fr/c0>1),因K0(fr/c0)很小,0(r,t)才有可能為負(fù),不過振幅很小。這說明在任一點(diǎn)上,其周圍不遠(yuǎn)地區(qū)的BΩ的貢獻(xiàn)是隨時(shí)間和隨地區(qū)積累的。若再設(shè)強(qiáng)迫海洋的風(fēng)是準(zhǔn)定常的或相對(duì)慣性時(shí)間尺度f-1變化得很緩慢(此時(shí)其屬大尺度大氣運(yùn)動(dòng),即其時(shí)間尺度大于1d,空間尺度T大于1d,空間k度大于2 000km,注意大尺度大氣的空間尺度典型值要遠(yuǎn)大于大尺度海洋),此時(shí)式(22)~(24)中第二個(gè)大括號(hào)中由BΩ所決定的這部分(該大括號(hào)中的第一項(xiàng))就沿著比較確定的方向緩慢地演變著,呈現(xiàn)準(zhǔn)平衡的狀態(tài),即是海洋中壓力梯度與科氏力之間的準(zhǔn)地轉(zhuǎn)平衡以及海洋Ekman體積輸送的散度場(chǎng)和風(fēng)應(yīng)力旋度場(chǎng)之間的準(zhǔn)平衡。再注意到,?0(r,t-t′)/?t′圍繞著0上下擺動(dòng)(參見圖1上曲線4),周期約略等于慣性振動(dòng)周期2π/f。這說明,各點(diǎn)周圍地區(qū)Bφ,Bu,Bu的貢獻(xiàn)在時(shí)間上和空間上是相互抵消的,由此可知,式(22)~(24)中第二個(gè)大括號(hào)中的這部分(該大括號(hào)中的第二項(xiàng))大體是沿時(shí)間振動(dòng)著的,而且周期短,其具有重力慣性波的特點(diǎn)。總之,外源強(qiáng)迫即大尺度風(fēng)應(yīng)力的作用是使流場(chǎng)和壓力場(chǎng)沿著某一方向進(jìn)行著上述準(zhǔn)平衡的演變,但同時(shí)又激發(fā)出環(huán)繞著該基本方向的短周期重力慣性波。

    圖1 Green函數(shù)隨時(shí)間的變化Fig.1 Variations of the Green function with time

    圖2 Green函數(shù)的分布Fig.2 Distribution of the Green function

    其中,

    注意到影響最大的是r=0附近的擾源,于是根據(jù)式(29)可知,若滿足以下的條件:的量級(jí)不小于的量級(jí),?/?y的量級(jí)不小于f的量級(jí),?/?x的量級(jí)不小于f的量級(jí),則在t>>f-1之后,從整體來說,在式(27)的兩項(xiàng)中以第一項(xiàng)占優(yōu),這種情況下外強(qiáng)迫即風(fēng)應(yīng)力所導(dǎo)致海洋響應(yīng)的總效果是比較緩慢的演變。此時(shí)海洋運(yùn)動(dòng)以準(zhǔn)平衡狀態(tài)占優(yōu)。流場(chǎng)的重力慣性波部分(式(27)中的第二項(xiàng))相對(duì)于流場(chǎng)的準(zhǔn)平衡部分(式(27)中的第一項(xiàng))而言,其處于次要的地位??紤]到該風(fēng)應(yīng)力激發(fā)出的重力慣性波還會(huì)向全平面彌散,而使其振幅衰減,且準(zhǔn)平衡的運(yùn)動(dòng)長(zhǎng)期保留,并隨時(shí)間緩慢演變,這樣海洋在該風(fēng)應(yīng)力強(qiáng)迫下的這部分運(yùn)動(dòng)就可分為“快過程”即重力慣性波和“慢過程”即準(zhǔn)平衡的演變過程兩部分。該風(fēng)應(yīng)力不斷激發(fā)出重力慣性波,同時(shí)使海洋沿著比較確定的方向比較緩慢地準(zhǔn)平衡演變,且激發(fā)出的重力慣性波會(huì)彌散掉,而準(zhǔn)平衡的演變則持續(xù)著。這就是海洋向該風(fēng)應(yīng)力的適應(yīng)以及其演變的情況。

    以上假設(shè)了風(fēng)是準(zhǔn)定常的或相對(duì)慣性時(shí)間尺度f-1變化很緩慢的(其屬大氣大尺度運(yùn)動(dòng)),而這種情況就是大尺度風(fēng)場(chǎng)強(qiáng)迫海洋的情況。在考慮到式(17)并注意到有后,則此時(shí)有

    式中,τ為大尺度大氣運(yùn)動(dòng)的時(shí)間尺度,即平流時(shí)間尺度(τ大于1d),并有τ>f-1,在此f-1是慣性振動(dòng)的時(shí)間尺度。此時(shí)有,

    注意到τ>f-1則有

    同樣,

    這里L(fēng)為大氣大尺度運(yùn)動(dòng)的空間尺度(L>2 000km)。此時(shí),

    故有,

    同理有,

    4 結(jié) 語

    本研究建立了一個(gè)考慮海面風(fēng)應(yīng)力強(qiáng)迫的正壓原始方程海洋模型,以探討海洋向風(fēng)應(yīng)力響應(yīng)的具體過程。將海洋運(yùn)動(dòng)對(duì)風(fēng)場(chǎng)的響應(yīng)分解成相對(duì)獨(dú)立的兩部分:第一部分是風(fēng)應(yīng)力作為直接強(qiáng)迫項(xiàng)驅(qū)動(dòng)的強(qiáng)迫慣性振蕩,在該部分壓力梯度不起作用;第二部分是風(fēng)應(yīng)力引起的Ekman抽吸,此時(shí)則有該抽吸引起的垂直運(yùn)動(dòng),并通過其進(jìn)一步影響海洋運(yùn)動(dòng),該部分有壓力梯度。因第一部分已有研究,故本文章只討論第二部分,并僅對(duì)大尺度風(fēng)場(chǎng)強(qiáng)迫海洋運(yùn)動(dòng)的情況做了解析研究。首先對(duì)風(fēng)應(yīng)力引起的Ekman抽吸做了估計(jì),得到在大尺度海洋Ekman抽吸正比于風(fēng)應(yīng)力旋度的結(jié)論。然后通過對(duì)該抽吸引起的該海洋模型控制方程組的求解,給出了海洋向大尺度風(fēng)場(chǎng)風(fēng)應(yīng)力響應(yīng)的具體形式:其是使流場(chǎng)和壓力場(chǎng)沿著某一方向進(jìn)行準(zhǔn)平衡的演變,并有海洋中壓力梯度與科氏力之間的準(zhǔn)地轉(zhuǎn)平衡和海洋Ekman體積輸送的散度場(chǎng)與風(fēng)應(yīng)力旋度場(chǎng)之間的準(zhǔn)平衡。與此同時(shí),有短周期的重力慣性波被激發(fā)出來。該重力慣性波會(huì)較快彌散掉,而準(zhǔn)平衡的演變則始終持續(xù)著。故海洋在大尺度風(fēng)場(chǎng)的強(qiáng)迫下其運(yùn)動(dòng)由以下三部分組成:第一部分是強(qiáng)迫慣性振蕩,其時(shí)間尺度是f-1,這是1個(gè)快過程;第二部分是比較緩慢演變著的準(zhǔn)平衡過程;第三部分是重力慣性波,這也是快過程,其時(shí)間尺度比慣性振蕩還要短。第二部分與第三部分相比,前者占優(yōu)且后者還要向全平面彌散,故最終海洋運(yùn)動(dòng)表現(xiàn)為強(qiáng)迫慣性振蕩和準(zhǔn)平衡演變過程。

    最后要指出的是,本研究只討論了海洋運(yùn)動(dòng)對(duì)大尺度風(fēng)場(chǎng)的響應(yīng),對(duì)于海洋運(yùn)動(dòng)對(duì)中小尺度風(fēng)場(chǎng)的響應(yīng),則在文獻(xiàn)[9]中做了介紹。

    參考文獻(xiàn)(References):

    [1] GILL A E.Atmosphere-ocean dynamics[M].Translating by ZHANG L Z,LE K T,ZHAO X Y.Beijing:Ocean press,1988.GILL A E.大氣-海洋動(dòng)力學(xué)[M].張立政,樂肯堂,趙徐懿,譯.北京:海洋出版社,1988.

    [2] QIAN Y F,WANG Q Q,ZHU B C.Responsive properties of oceanic currents in the south china sea to wind stress of winter monsoon[J].Scientia Meteorologica Sinica.2000,20(1):1-8.錢永甫,王謙謙,朱伯承.南海海流對(duì)冬季風(fēng)風(fēng)應(yīng)力的響應(yīng)特征[J].氣象科學(xué),2000,20(1):1-8.

    [3] CAI S Q,SU J L,GAN Z J.Response of south china sea upper circulation to monsoon transition[J].Journal of Troical Oceanography.2001,20(1):52-60.蔡樹群,蘇紀(jì)蘭,甘子鈞.南海上層環(huán)流對(duì)季風(fēng)轉(zhuǎn)變的響應(yīng)[J].熱帶海洋學(xué)報(bào),2001,20(1):52-60.

    [4] ZHAO B R.A study of the effects of local wind on currents in the off-shore areas of the east china sea and huanghai sea[J].Chinese Journal of Oceanology and Limnology,1982,13(6):479-490.趙保仁.局地風(fēng)對(duì)黃海和東海近岸淺海海流影響的研究[J].海洋與湖沼,1982,13(6):479-490.

    [5] QIU B.Kuroshio Extension variability and forcing of the Pacific decadal oscillations:responses and potential feedback[J].J.Phys.Oceanogr.,2003,33:2465-2482.

    [6] CABANES C,HUCK T,VERDEREE C D.Contributions of wind forcing and surface heating to interannual sea level variations in the Atlantic Ocean[J].J.Phys.Oceanogr.,2006,36:1739-1750.

    [7] QIU B,CHEN S.Decadal variability in the large-scale sea surface height field of the South Pacific Ocean:Observations and causes[J].J.Phys.Oceanogr.,2006,36:1751-1762.

    [8] ZHANG Y C,LU K C,ZHANG M.Barotropic quasi-equilibrium ocean model and its solution.PartⅠ:Mid-latitude forcing of large-scale wind field[J].Journal of the Meteorological Sciences Scientia Meteorologica Sinica,2011,31(1):11-16.張永垂,路凱程,張銘.正壓準(zhǔn)平衡海洋模型及其解Ⅰ-中緯度大尺度風(fēng)場(chǎng)強(qiáng)迫情況[J].氣象科學(xué),2011,31(1):11-16.

    [9] ZHAO Y L.Analytical study on response of ocean to wind stress and high resolution numerical research on sea wave and current in the Bohai Sea,the Yellow Sea and the East China Sea[D].Nanjing:Doctor degree thesis of PLA univ.of sci.and eng.,2008.5.趙艷玲.海洋對(duì)風(fēng)應(yīng)力響應(yīng)的解析研究及渤黃東海浪流的高分辨率數(shù)值研究[D].南京:中國(guó)人民解放軍理工大學(xué),2008.

    [10] ZENG Q C.Numerical physical basis of numerical weather prediction I[M].Beijing:Science press,1979.曾慶存.數(shù)值天氣預(yù)報(bào)的數(shù)值物理基礎(chǔ)(第一卷)[M].北京:科學(xué)出版社,1979.

    (杜素蘭 編輯)

    Analytic Research on the Response of Ocean Movement to Large-scale Wind Field

    ZHAO Yan-ling1,3,DENG Bing2,3,ZHANG Ming3,LIU Tian-jun2
    (1.Unit 61741,P.L.A,Beijing 100094;2.Beijing Institute of Applied Meteorology,Beijing 100029;3.Laboratory of Atmospheric Circulation and Short-range Climate Forecast,Meteorological College,P.L.A.University of Science and Technology,Nanjing 211101)

    A barotropic primitive ocean model forced by the sea surface wind stress is established for discussing the response processes of the ocean to the large-scale wind stress,and the Green function is used to solve the analytic solutions of the ocean equations.It has been found that the strength of the Ekman pumping induced by wind stress is positively proportional to the curl of wind stress,and that the response form of the ocean to the wind stress can be described as a quasi-equilibrium evolution between the flow field and the pressure field in the ocean,which includes the quasi-geostrophic equilibrium between the pressure gradient in the ocean and the Coriolis force and the quasi-equilibrium between the divergence of Ekman transportation and the wind stress curl.Meanwhile,inertia gravity waves with short period will be excited out,but disperse very soon.However,the quasi-equilibrium evolution is always sustained.The ocean movement forced by large-scale wind field consists of the following three parts:the inertia gravity waves,the forced inertial oscillation and the quasi-equilibrium evolution.The first two parts are quick processes,with the first one being quicker than the second,and the third one is a slow process.The final performance of the ocean movement is the superimpostion of the forced inertial oscillation and the quasi-equilibrium evolution.

    large-scale wind field;ocean model;response;Ekman pumping

    October 11,2011

    P732

    A

    1671-6647(2012)03-0328-10

    2011-10-11

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目——北太平洋副熱帶環(huán)流變異對(duì)大氣強(qiáng)迫的反饋?zhàn)饔茫?007CB411805)

    趙艷玲(1980-),女 ,河北唐山人,工程師,博士,主要從事物理海洋方面研究.E-mail:bbygg99@yahoo.com.cn

    猜你喜歡
    壓力梯度時(shí)間尺度風(fēng)場(chǎng)
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對(duì)稱性與守恒量
    基于FLUENT的下?lián)舯┝魅S風(fēng)場(chǎng)建模
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    “最美風(fēng)場(chǎng)”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    大連市暴雨多時(shí)間尺度研究分析
    壓力梯度在油田開發(fā)中的應(yīng)用探討
    側(cè)向風(fēng)場(chǎng)中無人機(jī)的飛行研究
    疊加原理不能求解含啟動(dòng)壓力梯度滲流方程
    致密砂巖啟動(dòng)壓力梯度數(shù)值的影響因素
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    久热久热在线精品观看| 国产爱豆传媒在线观看| 精品一区二区免费观看| 69av精品久久久久久| 亚洲精品成人av观看孕妇| 午夜精品一区二区三区免费看| 成人高潮视频无遮挡免费网站| 男人舔奶头视频| 久久久国产一区二区| 欧美成人一区二区免费高清观看| 中文在线观看免费www的网站| 白带黄色成豆腐渣| 免费观看a级毛片全部| 男女边吃奶边做爰视频| 韩国av在线不卡| 久久精品熟女亚洲av麻豆精品 | 3wmmmm亚洲av在线观看| 91在线精品国自产拍蜜月| 最近的中文字幕免费完整| 青春草视频在线免费观看| 日本免费在线观看一区| 亚洲18禁久久av| 18禁裸乳无遮挡免费网站照片| 中文字幕免费在线视频6| 欧美激情国产日韩精品一区| 欧美日韩精品成人综合77777| 国产综合精华液| 国产精品一区www在线观看| 国产一区有黄有色的免费视频 | 国产黄色小视频在线观看| 人妻系列 视频| 国产激情偷乱视频一区二区| 男女啪啪激烈高潮av片| 亚洲av电影在线观看一区二区三区 | 色播亚洲综合网| 免费看日本二区| 亚洲性久久影院| 日产精品乱码卡一卡2卡三| 精品久久国产蜜桃| 久久精品国产亚洲网站| 少妇熟女aⅴ在线视频| 国产午夜精品久久久久久一区二区三区| 激情五月婷婷亚洲| 建设人人有责人人尽责人人享有的 | 观看美女的网站| 麻豆乱淫一区二区| 亚洲在久久综合| 国产三级在线视频| 日日干狠狠操夜夜爽| 国产单亲对白刺激| 波野结衣二区三区在线| 国产精品一及| 99久久精品国产国产毛片| 亚洲国产精品国产精品| 全区人妻精品视频| 天堂√8在线中文| 国产精品熟女久久久久浪| 免费黄频网站在线观看国产| 成人毛片60女人毛片免费| 五月玫瑰六月丁香| 亚洲成人久久爱视频| 久久这里有精品视频免费| 22中文网久久字幕| 免费看美女性在线毛片视频| 精品久久久久久成人av| 全区人妻精品视频| 成人特级av手机在线观看| 亚洲成人久久爱视频| 日本三级黄在线观看| 搞女人的毛片| eeuss影院久久| av在线播放精品| 如何舔出高潮| 在线天堂最新版资源| 久久久久久伊人网av| 美女内射精品一级片tv| 亚洲av在线观看美女高潮| av福利片在线观看| 亚洲综合色惰| 麻豆成人午夜福利视频| 欧美潮喷喷水| 黄片无遮挡物在线观看| 亚洲精品影视一区二区三区av| 国产久久久一区二区三区| av福利片在线观看| 国产精品国产三级国产av玫瑰| 一本久久精品| 国产精品人妻久久久久久| 国产极品天堂在线| 免费看美女性在线毛片视频| 中文在线观看免费www的网站| 精品熟女少妇av免费看| 777米奇影视久久| 能在线免费观看的黄片| 人妻一区二区av| 亚洲国产精品sss在线观看| 国产黄片视频在线免费观看| 极品少妇高潮喷水抽搐| 国产爱豆传媒在线观看| 一区二区三区四区激情视频| 中文字幕人妻熟人妻熟丝袜美| 国产熟女欧美一区二区| 国产亚洲午夜精品一区二区久久 | 中文资源天堂在线| 热99在线观看视频| 国产人妻一区二区三区在| 成人午夜精彩视频在线观看| 神马国产精品三级电影在线观看| 偷拍熟女少妇极品色| 男人狂女人下面高潮的视频| 波野结衣二区三区在线| 夫妻午夜视频| 中文字幕免费在线视频6| 人体艺术视频欧美日本| 97在线视频观看| 人妻系列 视频| 国产一区二区三区综合在线观看 | 亚洲国产精品sss在线观看| 色播亚洲综合网| 久久久a久久爽久久v久久| 日本黄大片高清| 成人特级av手机在线观看| 欧美3d第一页| 99视频精品全部免费 在线| 丰满人妻一区二区三区视频av| 久久久色成人| 久久鲁丝午夜福利片| 久久这里只有精品中国| 熟女人妻精品中文字幕| 成年版毛片免费区| 日韩av免费高清视频| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av涩爱| 男人舔女人下体高潮全视频| 久久久久网色| 国产成人freesex在线| 亚洲欧洲国产日韩| 午夜免费激情av| 亚洲三级黄色毛片| 哪个播放器可以免费观看大片| 久久韩国三级中文字幕| 日韩电影二区| 亚洲欧洲国产日韩| 麻豆成人午夜福利视频| 99久久九九国产精品国产免费| 日韩大片免费观看网站| 午夜激情福利司机影院| 欧美不卡视频在线免费观看| 国产精品一区二区三区四区久久| 国产精品一二三区在线看| 女人久久www免费人成看片| 久久久久久久久大av| 亚洲精品日本国产第一区| 欧美高清成人免费视频www| eeuss影院久久| 午夜福利视频1000在线观看| 精品少妇黑人巨大在线播放| 国产精品国产三级国产av玫瑰| 国产美女午夜福利| 免费看光身美女| 亚洲伊人久久精品综合| 免费在线观看成人毛片| 夫妻性生交免费视频一级片| 夫妻午夜视频| 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 丰满少妇做爰视频| 日韩一区二区视频免费看| 2021天堂中文幕一二区在线观| 天天躁日日操中文字幕| 看黄色毛片网站| 日韩精品有码人妻一区| 熟女人妻精品中文字幕| 大陆偷拍与自拍| 色哟哟·www| 亚洲精品中文字幕在线视频 | 国产有黄有色有爽视频| av一本久久久久| 如何舔出高潮| 久久久久久久久久人人人人人人| 国产高清国产精品国产三级 | 久久这里只有精品中国| 国产伦在线观看视频一区| 精品一区二区三区人妻视频| 亚洲精品久久久久久婷婷小说| 成人国产麻豆网| 夫妻午夜视频| 久久97久久精品| 午夜激情欧美在线| 男女那种视频在线观看| 免费看光身美女| 午夜福利在线在线| 日韩欧美三级三区| 97热精品久久久久久| 国产91av在线免费观看| 免费看光身美女| 汤姆久久久久久久影院中文字幕 | 嫩草影院入口| 久久久久久久大尺度免费视频| 久久久欧美国产精品| 狠狠精品人妻久久久久久综合| 久久久久久久久久人人人人人人| 欧美日韩视频高清一区二区三区二| 国产在视频线在精品| 纵有疾风起免费观看全集完整版 | 男人爽女人下面视频在线观看| 少妇的逼水好多| 午夜福利网站1000一区二区三区| 偷拍熟女少妇极品色| 亚洲天堂国产精品一区在线| 美女大奶头视频| 亚洲欧美清纯卡通| 天天一区二区日本电影三级| 激情五月婷婷亚洲| 亚洲在线观看片| 亚洲自拍偷在线| 国语对白做爰xxxⅹ性视频网站| 日韩在线高清观看一区二区三区| 国产av不卡久久| 免费看日本二区| 三级国产精品欧美在线观看| 国精品久久久久久国模美| 亚洲人成网站高清观看| 久久这里只有精品中国| 久久草成人影院| 国产成人精品久久久久久| 国产精品女同一区二区软件| 午夜免费观看性视频| 热99在线观看视频| 干丝袜人妻中文字幕| av.在线天堂| 日本爱情动作片www.在线观看| 国产综合懂色| 国产亚洲5aaaaa淫片| 国语对白做爰xxxⅹ性视频网站| 人妻一区二区av| xxx大片免费视频| 欧美日韩亚洲高清精品| 国产爱豆传媒在线观看| 国内精品宾馆在线| 人妻制服诱惑在线中文字幕| 夫妻午夜视频| 亚洲人成网站高清观看| www.色视频.com| 免费在线观看成人毛片| 1000部很黄的大片| 午夜福利视频1000在线观看| 在线a可以看的网站| 中国国产av一级| 亚洲在线自拍视频| 久久久久久久久大av| 婷婷色av中文字幕| 亚洲精品自拍成人| 亚洲va在线va天堂va国产| 熟妇人妻久久中文字幕3abv| 国产高清国产精品国产三级 | 久久久久久久久久久丰满| 久久99精品国语久久久| 美女cb高潮喷水在线观看| 看免费成人av毛片| 欧美性感艳星| 18+在线观看网站| 亚洲av免费在线观看| 午夜精品国产一区二区电影 | 少妇丰满av| 99热全是精品| 午夜福利高清视频| 国内精品宾馆在线| 一级毛片黄色毛片免费观看视频| 久久韩国三级中文字幕| 少妇熟女aⅴ在线视频| 丰满乱子伦码专区| 国产色婷婷99| 日韩三级伦理在线观看| 中文在线观看免费www的网站| av免费观看日本| 久99久视频精品免费| 免费av不卡在线播放| 一级二级三级毛片免费看| 成人美女网站在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 成人亚洲欧美一区二区av| 国产黄色小视频在线观看| 波多野结衣巨乳人妻| 国产 亚洲一区二区三区 | 国产69精品久久久久777片| 1000部很黄的大片| 日韩亚洲欧美综合| 亚洲电影在线观看av| 丝袜美腿在线中文| 亚洲欧美日韩东京热| 丰满乱子伦码专区| 精品久久久久久久人妻蜜臀av| 七月丁香在线播放| 99久久人妻综合| 欧美日韩在线观看h| 能在线免费观看的黄片| 成人毛片a级毛片在线播放| 午夜免费激情av| 免费无遮挡裸体视频| 18禁在线播放成人免费| 久久久精品免费免费高清| 美女黄网站色视频| 国产久久久一区二区三区| 日本色播在线视频| 干丝袜人妻中文字幕| 午夜日本视频在线| 在线免费十八禁| 你懂的网址亚洲精品在线观看| av卡一久久| 久久鲁丝午夜福利片| 国产精品女同一区二区软件| 亚洲无线观看免费| 最近视频中文字幕2019在线8| 日韩欧美一区视频在线观看 | 高清视频免费观看一区二区 | 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 天天一区二区日本电影三级| 天美传媒精品一区二区| 尾随美女入室| 青春草亚洲视频在线观看| 久久精品熟女亚洲av麻豆精品 | 亚洲av在线观看美女高潮| 精品久久久久久久久亚洲| 韩国高清视频一区二区三区| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 中文字幕av成人在线电影| 亚洲最大成人手机在线| av在线蜜桃| 午夜老司机福利剧场| 国产精品.久久久| 成人午夜高清在线视频| 日韩,欧美,国产一区二区三区| 青春草亚洲视频在线观看| 亚洲人成网站在线观看播放| 舔av片在线| 少妇的逼水好多| 夫妻午夜视频| 国产精品人妻久久久久久| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 久久精品人妻少妇| 久久久久久久久久成人| 免费观看a级毛片全部| 亚洲精品久久午夜乱码| 国产乱人偷精品视频| 2022亚洲国产成人精品| 嫩草影院入口| 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 观看美女的网站| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 观看美女的网站| 国产 一区 欧美 日韩| 亚洲精品乱码久久久久久按摩| 日本爱情动作片www.在线观看| 国产精品麻豆人妻色哟哟久久 | 成人无遮挡网站| 成人午夜高清在线视频| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 日韩欧美精品免费久久| 精品不卡国产一区二区三区| 亚洲四区av| 亚洲电影在线观看av| 高清欧美精品videossex| 久久久久久久久久久免费av| 欧美性猛交╳xxx乱大交人| 精品久久久噜噜| 国产男人的电影天堂91| 精品99又大又爽又粗少妇毛片| 国产爱豆传媒在线观看| 久久精品国产亚洲网站| 美女大奶头视频| 特大巨黑吊av在线直播| 亚洲欧美成人综合另类久久久| 亚洲精品乱码久久久v下载方式| 亚洲欧美成人综合另类久久久| 亚洲精品乱码久久久v下载方式| 日韩制服骚丝袜av| 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| 青春草视频在线免费观看| 国产亚洲91精品色在线| 青春草视频在线免费观看| 亚洲内射少妇av| 精品久久久久久久人妻蜜臀av| 国产视频首页在线观看| 人妻一区二区av| 日本av手机在线免费观看| 免费看光身美女| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久久久久免费av| 一本久久精品| 男的添女的下面高潮视频| 三级毛片av免费| 久久久午夜欧美精品| 99热这里只有是精品在线观看| 亚洲成人精品中文字幕电影| 欧美成人精品欧美一级黄| 26uuu在线亚洲综合色| 狠狠精品人妻久久久久久综合| 日本wwww免费看| 直男gayav资源| 色综合亚洲欧美另类图片| 别揉我奶头 嗯啊视频| 一级毛片 在线播放| 99热这里只有是精品50| 麻豆国产97在线/欧美| 人人妻人人看人人澡| 嫩草影院精品99| 内射极品少妇av片p| 91久久精品国产一区二区成人| 国产男女超爽视频在线观看| 日韩亚洲欧美综合| freevideosex欧美| 女的被弄到高潮叫床怎么办| 九九久久精品国产亚洲av麻豆| 最近最新中文字幕大全电影3| 亚洲av福利一区| 午夜视频国产福利| 国产精品女同一区二区软件| 成人亚洲精品一区在线观看 | xxx大片免费视频| 天堂av国产一区二区熟女人妻| 国产亚洲最大av| 亚洲无线观看免费| 男女视频在线观看网站免费| 夜夜看夜夜爽夜夜摸| 成人性生交大片免费视频hd| av播播在线观看一区| 久久久久久久国产电影| 日韩制服骚丝袜av| 精品国内亚洲2022精品成人| 嫩草影院精品99| 久久鲁丝午夜福利片| av国产久精品久网站免费入址| 亚洲第一区二区三区不卡| 国产美女午夜福利| av免费观看日本| 精品亚洲乱码少妇综合久久| 国产成人福利小说| av在线亚洲专区| 欧美日韩视频高清一区二区三区二| 精品酒店卫生间| 国产在视频线精品| 亚洲激情五月婷婷啪啪| 18禁动态无遮挡网站| 亚洲av.av天堂| 天堂网av新在线| av一本久久久久| 日本熟妇午夜| 午夜福利视频精品| 国产午夜福利久久久久久| 国产精品一区www在线观看| 自拍偷自拍亚洲精品老妇| 免费看日本二区| 99热网站在线观看| 中文字幕久久专区| av在线天堂中文字幕| 熟妇人妻久久中文字幕3abv| 久久97久久精品| av.在线天堂| 亚洲精品日本国产第一区| 亚洲精品日韩在线中文字幕| 欧美 日韩 精品 国产| av免费在线看不卡| 成年版毛片免费区| 两个人视频免费观看高清| 国产av码专区亚洲av| av女优亚洲男人天堂| 亚洲欧美精品专区久久| 天堂影院成人在线观看| 免费电影在线观看免费观看| 免费av观看视频| 卡戴珊不雅视频在线播放| 国产av码专区亚洲av| 亚洲欧美一区二区三区国产| 婷婷色综合大香蕉| 成人美女网站在线观看视频| 超碰av人人做人人爽久久| 成年版毛片免费区| 简卡轻食公司| 少妇被粗大猛烈的视频| 亚洲av成人精品一区久久| 精品一区在线观看国产| h日本视频在线播放| 日本欧美国产在线视频| 国产亚洲5aaaaa淫片| 丝瓜视频免费看黄片| 久久久久久伊人网av| 欧美变态另类bdsm刘玥| 两个人视频免费观看高清| 亚洲第一区二区三区不卡| 99热这里只有是精品在线观看| 久久这里只有精品中国| 99久久精品一区二区三区| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2| 久久99热6这里只有精品| 国产精品久久久久久精品电影| 午夜免费观看性视频| 亚洲欧美日韩卡通动漫| 国产女主播在线喷水免费视频网站 | 精品酒店卫生间| 老司机影院成人| 免费看a级黄色片| 亚洲图色成人| 精品一区在线观看国产| .国产精品久久| 亚洲国产色片| 3wmmmm亚洲av在线观看| 久久精品熟女亚洲av麻豆精品 | 亚洲欧美成人精品一区二区| av国产久精品久网站免费入址| 久久精品人妻少妇| 别揉我奶头 嗯啊视频| 欧美另类一区| 夫妻性生交免费视频一级片| 午夜激情欧美在线| 亚洲国产av新网站| 人人妻人人看人人澡| 亚洲婷婷狠狠爱综合网| freevideosex欧美| 极品少妇高潮喷水抽搐| 国产精品精品国产色婷婷| 看十八女毛片水多多多| 97精品久久久久久久久久精品| 少妇丰满av| 秋霞伦理黄片| 91精品伊人久久大香线蕉| av在线蜜桃| 亚洲av不卡在线观看| 国产黄a三级三级三级人| 国产一区二区三区综合在线观看 | 亚洲av成人精品一二三区| 国产久久久一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 丝瓜视频免费看黄片| 日韩三级伦理在线观看| 不卡视频在线观看欧美| av女优亚洲男人天堂| 91av网一区二区| 天堂网av新在线| 一区二区三区四区激情视频| 国产精品综合久久久久久久免费| 内地一区二区视频在线| 淫秽高清视频在线观看| 亚洲成人av在线免费| 少妇人妻精品综合一区二区| 亚洲av男天堂| 国产精品99久久久久久久久| 高清欧美精品videossex| 国产 一区精品| 成人综合一区亚洲| 国产免费视频播放在线视频 | 美女高潮的动态| 爱豆传媒免费全集在线观看| 日本熟妇午夜| 亚洲成人久久爱视频| 热99在线观看视频| 日韩在线高清观看一区二区三区| 精品久久久精品久久久| 美女黄网站色视频| kizo精华| 搡老妇女老女人老熟妇| 黄色一级大片看看| 国内精品宾馆在线| 国产黄色免费在线视频| 内地一区二区视频在线| 激情 狠狠 欧美| 日本免费在线观看一区| 晚上一个人看的免费电影| 男女边吃奶边做爰视频| 亚洲成人av在线免费| 免费看av在线观看网站| 床上黄色一级片| 久久久久久久久久久免费av| 人妻系列 视频| 五月玫瑰六月丁香| 欧美不卡视频在线免费观看| 免费黄色在线免费观看| 亚洲精品乱久久久久久| 日韩,欧美,国产一区二区三区| 精品人妻熟女av久视频| 99久久精品热视频| 91精品一卡2卡3卡4卡| 亚洲va在线va天堂va国产| 免费观看a级毛片全部| 日韩欧美一区视频在线观看 | 久久精品人妻少妇| 国产大屁股一区二区在线视频| 青春草亚洲视频在线观看| 亚洲国产高清在线一区二区三| 久久草成人影院| 久久精品综合一区二区三区| 一级毛片黄色毛片免费观看视频| 特大巨黑吊av在线直播| 欧美另类一区| 亚洲第一区二区三区不卡| 久久久久久久国产电影| 听说在线观看完整版免费高清| av在线蜜桃| 久久精品国产鲁丝片午夜精品| 天美传媒精品一区二区| 男女啪啪激烈高潮av片| 国产精品麻豆人妻色哟哟久久 | 成人毛片a级毛片在线播放| 尤物成人国产欧美一区二区三区| 一夜夜www| 国产黄色免费在线视频| 日韩视频在线欧美|