• <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
    久久这里只有精品19| 人人妻人人添人人爽欧美一区卜| 看免费成人av毛片| 久久99一区二区三区| 欧美日韩一级在线毛片| 黑人欧美特级aaaaaa片| a级毛片黄视频| 高清在线视频一区二区三区| 老汉色av国产亚洲站长工具| 男女边摸边吃奶| 国产男人的电影天堂91| 欧美 亚洲 国产 日韩一| 三上悠亚av全集在线观看| 亚洲第一青青草原| 亚洲精品久久久久久婷婷小说| 日日撸夜夜添| 国产精品久久久久成人av| 国产深夜福利视频在线观看| 天天添夜夜摸| 久久久久精品久久久久真实原创| 亚洲国产欧美网| 中文字幕制服av| 亚洲久久久国产精品| 亚洲精品aⅴ在线观看| 爱豆传媒免费全集在线观看| 美女国产高潮福利片在线看| h视频一区二区三区| 亚洲av日韩在线播放| 色吧在线观看| 欧美国产精品va在线观看不卡| 精品人妻一区二区三区麻豆| 一级毛片黄色毛片免费观看视频| 亚洲综合色网址| 国产成人a∨麻豆精品| 欧美 亚洲 国产 日韩一| 美女国产高潮福利片在线看| 亚洲五月色婷婷综合| 91精品伊人久久大香线蕉| 2021少妇久久久久久久久久久| 一级毛片我不卡| 男人爽女人下面视频在线观看| 久久 成人 亚洲| 婷婷色综合www| 久久国产亚洲av麻豆专区| 人人妻人人添人人爽欧美一区卜| 99九九在线精品视频| 久热爱精品视频在线9| 午夜精品国产一区二区电影| 夫妻性生交免费视频一级片| 亚洲国产精品一区二区三区在线| 别揉我奶头~嗯~啊~动态视频 | 一区二区三区乱码不卡18| 欧美日韩av久久| 欧美成人精品欧美一级黄| 精品亚洲乱码少妇综合久久| 亚洲人成网站在线观看播放| 女性生殖器流出的白浆| 国产一区二区激情短视频 | 啦啦啦中文免费视频观看日本| 我要看黄色一级片免费的| 亚洲欧美日韩另类电影网站| 国产乱人偷精品视频| 国产欧美日韩一区二区三区在线| 亚洲av在线观看美女高潮| 新久久久久国产一级毛片| 青春草视频在线免费观看| 精品人妻在线不人妻| 亚洲欧美成人精品一区二区| 午夜福利免费观看在线| 亚洲激情五月婷婷啪啪| 99热网站在线观看| 色综合欧美亚洲国产小说| 丰满少妇做爰视频| 一个人免费看片子| 日本色播在线视频| 男男h啪啪无遮挡| 亚洲美女搞黄在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲午夜精品一区二区久久| 一边摸一边做爽爽视频免费| 九草在线视频观看| 亚洲一级一片aⅴ在线观看| 美国免费a级毛片| 精品亚洲成国产av| 午夜福利网站1000一区二区三区| 欧美少妇被猛烈插入视频| 国产精品久久久久久精品电影小说| h视频一区二区三区| 一级,二级,三级黄色视频| 成人三级做爰电影| svipshipincom国产片| 色婷婷久久久亚洲欧美| 自线自在国产av| 国产伦理片在线播放av一区| 国产精品国产三级专区第一集| 国产精品女同一区二区软件| 国产爽快片一区二区三区| 国产精品一区二区在线观看99| 人人妻人人澡人人看| 青青草视频在线视频观看| 日韩欧美精品免费久久| 亚洲五月色婷婷综合| 久久人人97超碰香蕉20202| 18禁动态无遮挡网站| 亚洲人成网站在线观看播放| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一国产av| 国产精品 欧美亚洲| 在线 av 中文字幕| 亚洲av成人不卡在线观看播放网 | 久久久久人妻精品一区果冻| 青青草视频在线视频观看| 99久久综合免费| 亚洲精品国产区一区二| 在线观看免费视频网站a站| 欧美日韩综合久久久久久| 午夜福利免费观看在线| 亚洲色图 男人天堂 中文字幕| av网站在线播放免费| 亚洲精品国产色婷婷电影| 亚洲精品国产av成人精品| 日韩一区二区三区影片| 国产成人系列免费观看| 五月开心婷婷网| 最新的欧美精品一区二区| 美女扒开内裤让男人捅视频| 9191精品国产免费久久| 性色av一级| 欧美在线黄色| 色网站视频免费| 在线观看免费午夜福利视频| 黑人猛操日本美女一级片| 王馨瑶露胸无遮挡在线观看| 久久久久精品国产欧美久久久 | 亚洲七黄色美女视频| 亚洲七黄色美女视频| 亚洲伊人久久精品综合| 亚洲精品久久午夜乱码| 少妇 在线观看| 亚洲久久久国产精品| 视频区图区小说| 夜夜骑夜夜射夜夜干| 精品人妻在线不人妻| 亚洲色图 男人天堂 中文字幕| 成人影院久久| 国产黄色视频一区二区在线观看| 男人舔女人的私密视频| 久久久欧美国产精品| 看十八女毛片水多多多| 亚洲人成电影观看| av网站在线播放免费| 一本久久精品| 日韩制服丝袜自拍偷拍| 捣出白浆h1v1| 99久久精品国产亚洲精品| 亚洲一区中文字幕在线| 国产在线视频一区二区| 国产男女内射视频| 亚洲美女搞黄在线观看| 日韩av在线免费看完整版不卡| 丰满迷人的少妇在线观看| 嫩草影视91久久| 亚洲七黄色美女视频| 久久影院123| 丝袜在线中文字幕| 人体艺术视频欧美日本| 91老司机精品| 日韩 欧美 亚洲 中文字幕| 又大又爽又粗| 人妻人人澡人人爽人人| 欧美变态另类bdsm刘玥| av片东京热男人的天堂| 激情视频va一区二区三区| 亚洲av日韩精品久久久久久密 | 国产成人欧美| 精品国产一区二区三区久久久樱花| www.精华液| 婷婷色av中文字幕| 国产女主播在线喷水免费视频网站| 精品第一国产精品| 精品一区二区三卡| 国产国语露脸激情在线看| 女性被躁到高潮视频| 亚洲男人天堂网一区| av片东京热男人的天堂| 美女视频免费永久观看网站| 性高湖久久久久久久久免费观看| 一本色道久久久久久精品综合| 九色亚洲精品在线播放| 91精品三级在线观看| 亚洲免费av在线视频| a级毛片黄视频| 老汉色av国产亚洲站长工具| www.熟女人妻精品国产| 自线自在国产av| 日韩人妻精品一区2区三区| 高清不卡的av网站| 午夜免费男女啪啪视频观看| 一本色道久久久久久精品综合| 亚洲三区欧美一区| 久久99一区二区三区| 在线观看国产h片| 精品人妻一区二区三区麻豆| 婷婷色综合大香蕉| 少妇 在线观看| 欧美最新免费一区二区三区| 日本av手机在线免费观看| 亚洲专区中文字幕在线 | 男女国产视频网站| 欧美日韩精品网址| 国产xxxxx性猛交| 中文字幕人妻丝袜制服| 菩萨蛮人人尽说江南好唐韦庄| 99久久综合免费| 亚洲少妇的诱惑av| 人体艺术视频欧美日本| 久久久久久久久免费视频了| 搡老岳熟女国产| 老鸭窝网址在线观看| 久久久久精品国产欧美久久久 | 欧美激情 高清一区二区三区| 视频在线观看一区二区三区| 国产精品 欧美亚洲| 亚洲欧美日韩另类电影网站| 男女床上黄色一级片免费看| 欧美日韩视频精品一区| 午夜福利视频精品| 国产精品久久久久久精品古装| 午夜福利,免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 91老司机精品| 蜜桃国产av成人99| www.av在线官网国产| 国产一区亚洲一区在线观看| 丝瓜视频免费看黄片| 日韩大码丰满熟妇| 成年人午夜在线观看视频| 日本黄色日本黄色录像| 久久热在线av| 国产午夜精品一二区理论片| 中文字幕色久视频| 亚洲精品一区蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费高清a一片| 天天躁狠狠躁夜夜躁狠狠躁| a 毛片基地| 欧美黑人精品巨大| 精品福利永久在线观看| 一级片'在线观看视频| 亚洲国产最新在线播放| tube8黄色片| 亚洲av欧美aⅴ国产| 99热网站在线观看| 午夜激情久久久久久久| 精品第一国产精品| 曰老女人黄片| 十八禁高潮呻吟视频| 大话2 男鬼变身卡| 国产精品亚洲av一区麻豆 | 久热这里只有精品99| a级毛片在线看网站| 亚洲av男天堂| 又大又爽又粗| 亚洲欧美成人综合另类久久久| 欧美国产精品一级二级三级| 亚洲精品成人av观看孕妇| 亚洲精品美女久久av网站| 少妇人妻久久综合中文| 男男h啪啪无遮挡| 晚上一个人看的免费电影| 女性生殖器流出的白浆| av免费观看日本| 看免费av毛片| 精品亚洲乱码少妇综合久久| 亚洲第一区二区三区不卡| 亚洲欧洲日产国产| 妹子高潮喷水视频| www.av在线官网国产| 一级毛片 在线播放| 日韩 亚洲 欧美在线| 高清视频免费观看一区二区| bbb黄色大片| 午夜激情av网站| 免费在线观看完整版高清| 人人妻人人澡人人爽人人夜夜| 人人妻,人人澡人人爽秒播 | 亚洲精品久久午夜乱码| 亚洲国产精品国产精品| 性高湖久久久久久久久免费观看| 一级毛片 在线播放| 80岁老熟妇乱子伦牲交| 亚洲国产日韩一区二区| 国产精品女同一区二区软件| 久久天堂一区二区三区四区| 丰满乱子伦码专区| 男人添女人高潮全过程视频| 777米奇影视久久| 国产在线一区二区三区精| 日本午夜av视频| 下体分泌物呈黄色| 亚洲成av片中文字幕在线观看| 成人黄色视频免费在线看| 日本爱情动作片www.在线观看| 久久人人97超碰香蕉20202| 夫妻性生交免费视频一级片| www.精华液| 丝袜在线中文字幕| 十八禁网站网址无遮挡| 国产一区二区在线观看av| 亚洲精品日本国产第一区| 宅男免费午夜| 亚洲欧洲国产日韩| 午夜激情av网站| 美女高潮到喷水免费观看| 成人18禁高潮啪啪吃奶动态图| 久久国产精品男人的天堂亚洲| 七月丁香在线播放| 亚洲少妇的诱惑av| 亚洲av男天堂| 无遮挡黄片免费观看| 久久这里只有精品19| 三上悠亚av全集在线观看| 日本av免费视频播放| 国产亚洲av高清不卡| 久久久亚洲精品成人影院| 人妻人人澡人人爽人人| 欧美日本中文国产一区发布| 久久97久久精品| 99re6热这里在线精品视频| 中文字幕最新亚洲高清| 欧美激情 高清一区二区三区| 男女高潮啪啪啪动态图| 亚洲精品一二三| av一本久久久久| 精品视频人人做人人爽| 国产亚洲午夜精品一区二区久久| 成人午夜精彩视频在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲第一av免费看| 久久久久视频综合| 巨乳人妻的诱惑在线观看| 亚洲国产av影院在线观看| 伊人久久大香线蕉亚洲五| 欧美日韩成人在线一区二区| 热re99久久精品国产66热6| 亚洲精品视频女| a 毛片基地| 国产精品久久久久久人妻精品电影 | 色婷婷av一区二区三区视频| 新久久久久国产一级毛片| 亚洲国产av新网站| 亚洲图色成人| 免费在线观看黄色视频的| 丰满饥渴人妻一区二区三| 亚洲一区二区三区欧美精品| 一区二区三区激情视频| 街头女战士在线观看网站| 久久久精品区二区三区| 免费观看a级毛片全部| 欧美亚洲日本最大视频资源| av福利片在线| 久久精品国产a三级三级三级| 在线观看免费午夜福利视频| 国产深夜福利视频在线观看| 一级毛片 在线播放| 精品亚洲成a人片在线观看| 精品一区二区免费观看| 欧美av亚洲av综合av国产av | 亚洲欧美色中文字幕在线| 色94色欧美一区二区| 国产欧美亚洲国产| 国产成人欧美在线观看 | 国产精品久久久久久人妻精品电影 | 欧美日韩视频精品一区| 久久韩国三级中文字幕| av视频免费观看在线观看| 午夜老司机福利片| 在线免费观看不下载黄p国产| 亚洲国产日韩一区二区| 国产精品 欧美亚洲| 91成人精品电影| 1024视频免费在线观看| 国产精品久久久久久精品电影小说| 纵有疾风起免费观看全集完整版| 下体分泌物呈黄色| 麻豆av在线久日| 日韩中文字幕欧美一区二区 | 黄色怎么调成土黄色| 男女床上黄色一级片免费看| 卡戴珊不雅视频在线播放| 国产在线免费精品| 一区二区三区乱码不卡18| av.在线天堂| www日本在线高清视频| 国产免费福利视频在线观看| 一二三四在线观看免费中文在| 亚洲精品日韩在线中文字幕| 久久久久久免费高清国产稀缺| 国产 一区精品| 纯流量卡能插随身wifi吗| 精品国产一区二区久久| 一级爰片在线观看| 一级毛片电影观看| 亚洲精品美女久久久久99蜜臀 | 999久久久国产精品视频| 亚洲美女视频黄频| 国产成人系列免费观看| 丰满乱子伦码专区| 深夜精品福利| 91精品伊人久久大香线蕉| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 最近最新中文字幕大全免费视频 | kizo精华| 亚洲色图综合在线观看| 成人亚洲欧美一区二区av| 久久国产亚洲av麻豆专区| 午夜免费观看性视频| 丝袜美腿诱惑在线| 亚洲国产欧美一区二区综合| 国产免费视频播放在线视频| 国产午夜精品一二区理论片| 亚洲成国产人片在线观看| 久久久久国产精品人妻一区二区| 久久久国产精品麻豆| 青青草视频在线视频观看| 80岁老熟妇乱子伦牲交| 777米奇影视久久| av线在线观看网站| 成人手机av| 老司机在亚洲福利影院| 最近中文字幕2019免费版| 免费在线观看视频国产中文字幕亚洲 | 高清欧美精品videossex| 国产免费福利视频在线观看| 韩国av在线不卡| 哪个播放器可以免费观看大片| 亚洲国产欧美日韩在线播放| 免费不卡黄色视频| 精品一品国产午夜福利视频| 日本猛色少妇xxxxx猛交久久| 亚洲人成网站在线观看播放| 亚洲 欧美一区二区三区| 黑丝袜美女国产一区| 亚洲精品日本国产第一区| 日韩大码丰满熟妇| 免费黄网站久久成人精品| 国产麻豆69| 欧美97在线视频| 久久国产亚洲av麻豆专区| kizo精华| 国产精品女同一区二区软件| 午夜福利在线免费观看网站| 亚洲国产欧美网| 欧美激情极品国产一区二区三区| 国产av国产精品国产| 久久国产精品大桥未久av| 制服人妻中文乱码| 宅男免费午夜| 观看美女的网站| 国产男女超爽视频在线观看| 国产 精品1| 如日韩欧美国产精品一区二区三区| 韩国精品一区二区三区| 国产男女内射视频| 日韩,欧美,国产一区二区三区| 99久久综合免费| 男女下面插进去视频免费观看| 青春草国产在线视频| 亚洲欧洲国产日韩| 成年动漫av网址| 色视频在线一区二区三区| 老司机影院毛片| a级毛片黄视频| 91国产中文字幕| 久久久久久人妻| 日韩中文字幕视频在线看片| 免费在线观看完整版高清| 大片免费播放器 马上看| 亚洲国产精品一区二区三区在线| 亚洲精品在线美女| 日韩一卡2卡3卡4卡2021年| 一区福利在线观看| 看免费av毛片| 国产成人精品久久久久久| 亚洲欧美日韩另类电影网站| 亚洲欧美成人综合另类久久久| 日韩视频在线欧美| 国产精品女同一区二区软件| 国产av精品麻豆| 亚洲欧美一区二区三区国产| 国产极品粉嫩免费观看在线| 美女中出高潮动态图| 久久热在线av| 欧美日韩亚洲综合一区二区三区_| 国产精品免费视频内射| av免费观看日本| 国产精品麻豆人妻色哟哟久久| 在线观看国产h片| 老司机影院毛片| 日韩 亚洲 欧美在线| 看非洲黑人一级黄片| 国产xxxxx性猛交| 人人澡人人妻人| 99国产综合亚洲精品| 侵犯人妻中文字幕一二三四区| 国产国语露脸激情在线看| 国产淫语在线视频| 亚洲色图综合在线观看| 成人毛片60女人毛片免费| 亚洲成人一二三区av| 国产成人欧美在线观看 | 午夜av观看不卡| 国产日韩欧美亚洲二区| 成人影院久久| 人体艺术视频欧美日本| 日韩一区二区三区影片| 精品一区二区三区四区五区乱码 | 女性生殖器流出的白浆| 欧美精品亚洲一区二区| 日韩精品有码人妻一区| 免费少妇av软件| 久久久久久久精品精品| 国产伦人伦偷精品视频| 蜜桃国产av成人99| 在线观看三级黄色| 在线天堂中文资源库| 中文乱码字字幕精品一区二区三区| 一区在线观看完整版| av在线老鸭窝| 观看美女的网站| 成年美女黄网站色视频大全免费| 亚洲精品中文字幕在线视频| 亚洲美女视频黄频| 免费看不卡的av| 国产免费一区二区三区四区乱码| 国产极品天堂在线| 丁香六月天网| 黄片播放在线免费| 天堂中文最新版在线下载| 婷婷成人精品国产| 高清av免费在线| 女人久久www免费人成看片| 国产精品嫩草影院av在线观看| av免费观看日本| 久久精品熟女亚洲av麻豆精品| 久久 成人 亚洲| 久久天躁狠狠躁夜夜2o2o | 伦理电影大哥的女人| 水蜜桃什么品种好| 婷婷色麻豆天堂久久| 亚洲国产欧美日韩在线播放| 在线观看免费午夜福利视频| av在线app专区| 美女大奶头黄色视频| 久久午夜综合久久蜜桃| 在线 av 中文字幕| 男女边摸边吃奶| 建设人人有责人人尽责人人享有的| 如日韩欧美国产精品一区二区三区| 亚洲专区中文字幕在线 | 麻豆精品久久久久久蜜桃| 日本一区二区免费在线视频| 丁香六月天网| 日韩av免费高清视频| 男女无遮挡免费网站观看| 欧美另类一区| 亚洲国产日韩一区二区| 2021少妇久久久久久久久久久| 汤姆久久久久久久影院中文字幕| 黄色 视频免费看| 极品人妻少妇av视频| 国产淫语在线视频| 国产精品女同一区二区软件| 丝袜美腿诱惑在线| 99久国产av精品国产电影| 亚洲av国产av综合av卡| 亚洲专区中文字幕在线 | 日韩一区二区视频免费看| 熟女少妇亚洲综合色aaa.| 亚洲,欧美,日韩| 久久精品国产综合久久久| 亚洲av欧美aⅴ国产| 老司机靠b影院| 大码成人一级视频| 在线精品无人区一区二区三| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 国产黄色视频一区二区在线观看| 别揉我奶头~嗯~啊~动态视频 | 久久精品久久精品一区二区三区| 成人亚洲欧美一区二区av| 色婷婷久久久亚洲欧美| 日本色播在线视频| 如何舔出高潮| 爱豆传媒免费全集在线观看| 人人妻人人爽人人添夜夜欢视频| 中文字幕高清在线视频| 母亲3免费完整高清在线观看| 国产精品蜜桃在线观看| 成人国产av品久久久| 亚洲国产日韩一区二区| 夜夜骑夜夜射夜夜干| 欧美亚洲 丝袜 人妻 在线| 亚洲美女搞黄在线观看| 日本猛色少妇xxxxx猛交久久| 久久久久精品人妻al黑| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美中文字幕日韩二区| 哪个播放器可以免费观看大片| 天天躁狠狠躁夜夜躁狠狠躁| 欧美久久黑人一区二区| 午夜福利在线免费观看网站| 狂野欧美激情性bbbbbb| 久久精品久久久久久噜噜老黄|