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

    加速度遷移項(xiàng)位勢(shì)及其在高原低渦及臺(tái)風(fēng)系統(tǒng)分析中的應(yīng)用

    2013-09-22 06:43:20高守亭周玉淑鄧滌菲
    地球物理學(xué)報(bào) 2013年2期
    關(guān)鍵詞:東移位勢(shì)平流

    高守亭,周玉淑,鄧滌菲.2

    1 中國(guó)科學(xué)院大氣物理研究所 云降水物理與強(qiáng)風(fēng)暴實(shí)驗(yàn)室(LACS),北京 100029

    2 中國(guó)科學(xué)院大學(xué),北京 100049

    1 引 言

    大氣運(yùn)動(dòng)中,水平風(fēng)除了有旋轉(zhuǎn)和輻合輻散(氣象上用渦度及散度分別度量其旋轉(zhuǎn)及輻合輻散的強(qiáng)度)特性外,水平風(fēng)場(chǎng)還有一個(gè)顯著的特性,就是平流效應(yīng),即水平風(fēng)場(chǎng)自身的平流作用及其對(duì)各個(gè)物理量的水平輸送,如:溫度平流、渦度平流、水汽平流等,都被廣泛用于天氣系統(tǒng)發(fā)生發(fā)展的分析.周明煜等[1]對(duì)1998年夏季第2次青藏高原大氣科學(xué)試驗(yàn)的觀測(cè)資料的分析表明,在對(duì)流層中平均風(fēng)速垂直分布存在風(fēng)切變現(xiàn)象,水平平流作用對(duì)邊界層結(jié)構(gòu)有重要影響.在已往的研究結(jié)果中,風(fēng)場(chǎng)平流作用更多的是用于對(duì)渦度方程的分析,如姚秀萍等[2]在研究熱帶對(duì)流層上空東風(fēng)帶影響西太平洋副熱帶高壓活動(dòng)時(shí),對(duì)垂直渦度方程進(jìn)行了分析,指出水平風(fēng)場(chǎng)對(duì)渦度輸送作用的貢獻(xiàn)最大.周玉淑等[3]和Zhou等[4]則是考慮了風(fēng)場(chǎng)平流效應(yīng)的旋轉(zhuǎn),推導(dǎo)得到了包含水平風(fēng)平流旋轉(zhuǎn)形式的渦度方程,在2006年的bilis臺(tái)風(fēng)及2008年的Fung-Wong臺(tái)風(fēng)分析中得到了很好的應(yīng)用,但是該研究只討論了水平平流旋轉(zhuǎn)效應(yīng)對(duì)渦度變化的貢獻(xiàn),并沒有考慮到平流的輻合輻散作用.楊平章等[5]在研究水平動(dòng)量平流作用對(duì)熱帶氣旋運(yùn)動(dòng)的影響時(shí),分析得到動(dòng)量平流引起的水平動(dòng)量輻合或輻散對(duì)臺(tái)風(fēng)路徑異常變化有明顯的影響.最近,高守亭等[6]的研究提出加速度遷移項(xiàng)(反映水平風(fēng)場(chǎng)平流的輻合輻散)能代表氣塊受到的所有力的散度,即通過(guò)計(jì)算加速度遷移項(xiàng)就可以判斷氣塊所受力的情況,這為用加速度遷移項(xiàng)來(lái)研究大氣運(yùn)動(dòng)提供了新的思路.由于加速度遷移項(xiàng)(也即是水平平流散度)與風(fēng)場(chǎng)本身的分布一樣,是非均勻的,且空間及時(shí)間變化很大,也就是說(shuō),水平平流的分布具有波動(dòng)特征,除了具有旋轉(zhuǎn)效應(yīng)外,其輻合輻散效應(yīng)也是一個(gè)重要方面.以往對(duì)水平平流作用及其旋轉(zhuǎn)效應(yīng)的研究還有一些結(jié)果,但在其輻合輻散效應(yīng)方面的研究尚屬少見.而大氣中的很多劇烈天氣現(xiàn)象(尤其是中尺度災(zāi)害天氣)往往是由風(fēng)場(chǎng)及其水汽平流的不均勻產(chǎn)生的輻合輻散引起的.因此,研究加速度散度遷移項(xiàng)(其反映了水平風(fēng)場(chǎng)的平流作用)是中尺度天氣研究中的一個(gè)重要內(nèi)容.

    根據(jù)二維亥姆霍茲定律,水平風(fēng)場(chǎng)可以分解為旋轉(zhuǎn)風(fēng)分量和輻散風(fēng)分量,引入風(fēng)場(chǎng)的流函數(shù)和速度勢(shì)后,就可以直接用流函數(shù)和速度勢(shì)分別表示風(fēng)場(chǎng)的旋轉(zhuǎn)和輻合輻散性質(zhì),已在全球大氣和大洋環(huán)流的診斷和數(shù)值模擬中有重要應(yīng)用.依照這一思路,本文對(duì)加速度遷移項(xiàng)也進(jìn)行有旋和無(wú)旋分解,可得到對(duì)應(yīng)的加速度遷移項(xiàng)的流函數(shù)和位勢(shì),利用加速度遷移項(xiàng)位勢(shì)就可以反映出加速度遷移項(xiàng)的輻合輻散運(yùn)動(dòng),擴(kuò)展了水平風(fēng)場(chǎng)分解技術(shù)在中尺度天氣研究中的應(yīng)用.根據(jù)已有的研究,Chen和Kuo[7]提出的調(diào)和-余弦算法在有限區(qū)域風(fēng)場(chǎng)分解和重建的精度最高,而且物理意義也比其它方法明確,在臺(tái)風(fēng)等低渦系統(tǒng)的風(fēng)場(chǎng)分析中已有很好的應(yīng)用[8-9].因此,本文采用這一方法,對(duì)加速度遷移項(xiàng)也進(jìn)行調(diào)和-余弦的有旋和無(wú)旋分量的分解,引入加速度遷移項(xiàng)位勢(shì)概念,以便利用加速度遷移項(xiàng)位勢(shì)來(lái)研究東移高原低渦及臺(tái)風(fēng)系統(tǒng)演變過(guò)程.

    2 加速度遷移項(xiàng)的調(diào)和-余弦分解方法介紹

    引入加速度遷移項(xiàng)的流函數(shù)及位勢(shì)概念后,有

    則加速度遷移項(xiàng)的流函數(shù)ψA及加速度遷移項(xiàng)位勢(shì)A與其旋轉(zhuǎn)部分的垂直分量和散度分量之間有如下關(guān)系:

    即:

    以上Ar和Ad為加速度遷移項(xiàng)的旋轉(zhuǎn)分量和輻散分,Δ ,Δ量 ψA為加速度遷移項(xiàng)的流函數(shù)A為加速度遷移項(xiàng)的位勢(shì),ξA和DA則分別是加速度遷移項(xiàng)的旋度的垂直分量和水平散度,計(jì)算公式為

    為得到有限區(qū)域泊松方程(6)和(7)的解,在耦合邊界條件下用調(diào)和-余弦算法[7]進(jìn)行求解,求解基本算法如下:

    (1)通過(guò)已有水平風(fēng)分量u和v計(jì)算得到加速度遷移項(xiàng)的兩個(gè)分量Ax和Ay后,計(jì)算區(qū)域里的加速度遷移項(xiàng)項(xiàng)的渦度ξA和散度DA;

    (2)利用其次邊界條件,求解公式(10)所示的Possion方程,該方Δ程是公式(6)和(7)內(nèi)部部分,即:

    通過(guò)傅里葉變換得到Possion方程(10)的解為

    其中,F(xiàn)-1是傅里葉逆變換算子,Φ和X是針對(duì)加速度遷移項(xiàng)的流函數(shù)和位勢(shì)作了有限區(qū)域傅里葉余弦變換得到的雙余弦展開式,具體表達(dá)式為

    其中下標(biāo)i表示內(nèi)部部分(inner part),下標(biāo)c表示余弦函數(shù)展開.

    (3)由計(jì)算出的ψAic和Aic,利用二維亥姆霍茲原理得到加速度遷移項(xiàng)的內(nèi)部分量AxIc和AyIc為

    (4)計(jì)算加速度遷移項(xiàng)的外部分量AxEc和AyEc

    (5)利用(4)得到的耦合邊界條件:

    來(lái)求解Laplace方程,即在(17)的邊條件下求解方程(18):

    其中,下標(biāo)h代表調(diào)和函數(shù)部分(harmonic part).

    通過(guò)迭代求解(具體系數(shù)的計(jì)算參見文獻(xiàn)[7]),最終得到加速度遷移項(xiàng)的流函數(shù)和位勢(shì)的調(diào)和分量為

    其中,ψW(0),ψW(n),…,ψE(0),ψE(n),… 都是邊界上的加速度遷移項(xiàng)的流函數(shù)和位勢(shì)的譜系數(shù),下標(biāo)E,S,W,N分別表示東、南、西、北四個(gè)方位.

    (6)由方程(11)、(12)和(19)、(20),得到計(jì)算的有限區(qū)域的加速度遷移項(xiàng)的流函數(shù)和位勢(shì)分量為

    這樣求解得到的物理量具有很明確的物理意義:加速度遷移項(xiàng)的流函數(shù)和位勢(shì)的內(nèi)部分量ψAic和完全由有限區(qū)域內(nèi)部的加速度遷移項(xiàng)的旋轉(zhuǎn)垂直分量及其散度決定,邊界條件都是齊次的,所以邊條件對(duì)內(nèi)部分量沒有影響;而調(diào)和分量ψAhc和Ahc是通過(guò)求解Laplace方程得到的,完全由區(qū)域的邊界條件決定,因而極值只出現(xiàn)在邊界上.在這樣的求解過(guò)程中,只需要通過(guò)極少的迭代次數(shù)就能達(dá)到很好的收斂效果,計(jì)算效率和精度都大大提高.

    由于本文關(guān)心的主要是加速度遷移項(xiàng)的輻合輻散作用對(duì)天氣系統(tǒng)的影響,以下的分析就用加速度遷移項(xiàng)位勢(shì)部分,即分布來(lái)討論其在具體天氣系統(tǒng)演變中的示蹤作用.

    3 加速度遷移項(xiàng)位勢(shì)在低渦及臺(tái)風(fēng)系統(tǒng)分析中的應(yīng)用

    由于平流的作用在大氣運(yùn)動(dòng)中無(wú)處無(wú)時(shí)不在,其輻合輻散作用必然會(huì)對(duì)系統(tǒng)演變有所指示.由于低渦及臺(tái)風(fēng)系統(tǒng)都是具有明顯的平流旋轉(zhuǎn)和輻合的系統(tǒng),本節(jié)就以一個(gè)東移高原低渦及登陸我國(guó)的Bilis臺(tái)風(fēng)(2006年第4號(hào)臺(tái)風(fēng))系統(tǒng)為例,分析高原低渦東移及臺(tái)風(fēng)系統(tǒng)演變過(guò)程中加速度遷移項(xiàng)位勢(shì)的變化,以期為示蹤及反應(yīng)低渦及臺(tái)風(fēng)系統(tǒng)的演變提供一個(gè)新的物理變量.

    3.1 在東移高原低渦中的應(yīng)用

    高原低渦是引發(fā)高原地區(qū)暴雨天氣的重要系統(tǒng)之一,在適當(dāng)?shù)沫h(huán)流條件配合下可以向東移出高原,并在適當(dāng)條件下發(fā)展,往往會(huì)造成周邊及下游地區(qū)(如江淮流域等)出現(xiàn)暴雨等災(zāi)害性天氣[10-12].基于高原低渦研究的需要,本節(jié)首先將加速度遷移項(xiàng)位勢(shì)用于2008-07-20-08∶00到21日14∶00的一次青藏高原低渦東移的個(gè)例分析中[13],所用資料為NCEP/NCAR的分析資料,空間分辨率為110km,時(shí)間分辨率為6h.首先分析這次高原低渦東移出高原的天氣過(guò)程.

    在500hPa天氣圖上,圖1a中中高緯度影響中國(guó)的主要槽脊系統(tǒng)是青藏高原北部的低壓槽(簡(jiǎn)稱西部槽)和貝加爾湖以南的高壓脊(簡(jiǎn)稱東部脊),其中西部槽比較深厚,而東部脊相對(duì)范圍比較小,且東部脊的東南部河套平原附近存在一個(gè)閉合的低壓中心.在40°N以南,高原東部附近(97°E,35°N)上空存在一個(gè)明顯的閉合低壓中心,即為本節(jié)分析的高原低渦.黃海上空也存在一個(gè)小的閉合低壓系統(tǒng).6小時(shí)后(20日14∶00,圖1b),中高緯地區(qū)槽脊系統(tǒng)東移,東部脊進(jìn)一步加強(qiáng),高壓脊東南部河套附近低壓中心東移加強(qiáng),位勢(shì)高度最小值降至5775位勢(shì)米.此時(shí)青藏高原東部高原渦東移至(98°E,35°N)附近,5805等高線所圍面積減小.20日20∶00(圖1c),河套附近低壓中心進(jìn)一步東移加強(qiáng),5775等高線范圍擴(kuò)大,而高原低渦東移至高原東邊緣附近,強(qiáng)度變化不大,此時(shí)高原渦位于河套低壓西南邊緣.至21日02∶00(圖1d),河套低壓東移至朝鮮半島,強(qiáng)度進(jìn)一步加強(qiáng),等高線最小值達(dá)5745位勢(shì)米,此時(shí)高原渦不再有閉合等高線,而是并入東移的河套低壓區(qū),以小槽形式存在,同時(shí)高原東側(cè)邊緣的低槽較前一時(shí)刻向南發(fā)展,5850位勢(shì)米等值線南伸到23°N.21日08∶00(圖1e),原河套低壓繼續(xù)增強(qiáng)向東移動(dòng),高原東側(cè)低槽也進(jìn)一步加深東移,且槽線由橫槽逐漸轉(zhuǎn)為豎槽,位置在105°E附近.到21日14∶00(圖1f),高原東側(cè)低槽強(qiáng)度加強(qiáng),移出高原后向東北方向移動(dòng).

    圖1 500hPa位勢(shì)高度場(chǎng)(單位:位勢(shì)米)和風(fēng)場(chǎng)(單位:m/s)分布(a)2008-07-20-08∶00;(b)2008-07-20-14∶00;(c)2008-07-20-20∶00;(d)2008-07-21-02∶00;(e)2008-07-21-08∶00;(f)2008-07-21-14∶00.

    從公式(7)可知,用余弦譜展開后,加速度遷移項(xiàng)位勢(shì)的拉普拉斯符號(hào)與平流本身的符號(hào)相反,也就是說(shuō),當(dāng)加速度遷移項(xiàng)位勢(shì)為正時(shí),其散度為負(fù),表示加速度遷移項(xiàng)是輻合的.從加速度遷移項(xiàng)位勢(shì)的分布來(lái)看,對(duì)2008年7月20日08∶00,500hPa上的高原東部,也是在(97°E,35°N)附近有明顯的加速度遷移項(xiàng)位勢(shì)的正值區(qū)(圖2a),與位勢(shì)高度分析顯示的高原低渦的位置相當(dāng),表示低渦處有明顯的加速度遷移項(xiàng)的輻合.從加速度遷移項(xiàng)位勢(shì)的正值區(qū)及其中心分布來(lái)看,它比500hPa位勢(shì)高度表示的低渦系統(tǒng)更為清楚.而且,黃海上空的閉合低壓系統(tǒng)及槽區(qū),也有明顯的加速度遷移項(xiàng)位勢(shì)的正值中心及正值區(qū)對(duì)應(yīng),這說(shuō)明用加速度遷移項(xiàng)位勢(shì)確實(shí)可以示蹤和識(shí)別低渦系統(tǒng).6小時(shí)后,高原低渦東移到(100°E,35°N)附近(圖1b),對(duì)應(yīng)時(shí)刻的加速度遷移項(xiàng)位勢(shì)的正值區(qū)及正值中心也東移到同樣的位置(圖2b),說(shuō)明低渦系統(tǒng)中確實(shí)一直存在加速度遷移項(xiàng)的輻合作用,也正是由于水平風(fēng)的加速度遷移項(xiàng)的輻合,低渦才能得以維持和發(fā)展.到7月20日20∶00和21日02∶00,低渦移到高原東部斜坡地形附近,雖然風(fēng)場(chǎng)上還存在氣旋性環(huán)流,但在位勢(shì)高度分布上已經(jīng)沒有閉合的等壓線,而是一個(gè)槽區(qū)(圖1c和圖1d),但在同時(shí)刻的加速度遷移項(xiàng)位勢(shì)分布上(圖2c和圖2d),在低槽區(qū)位置,仍然有明顯的加速度遷移項(xiàng)位勢(shì)正值區(qū)和正值中心,對(duì)指示東移高原低渦的位置比位勢(shì)高度更為清晰.隨著低渦并入低槽并移出高原后,槽區(qū)范圍加強(qiáng)南伸并繼續(xù)東移,對(duì)應(yīng)的加速度遷移項(xiàng)位勢(shì)分布圖上(圖2e和圖2f)也有加速度遷移項(xiàng)位勢(shì)范圍的擴(kuò)展和加強(qiáng),與位勢(shì)高度反映的趨勢(shì)是一致的,說(shuō)明用加速度遷移項(xiàng)位勢(shì)確實(shí)可以示蹤高原低渦的東移演變過(guò)程.

    圖2 500hPa加速度遷移項(xiàng)位勢(shì)(單位:m2·s-2)分布(a)2008-07-20-08∶00;(b)2008-07-20-14時(shí);(c)2008-07-20-20∶00;(d)2008-07-21日02∶00;(e)2008-07-21日08∶00;(f)2008-07-21-14∶00.

    圖3 500hPa加速度遷移項(xiàng)位勢(shì)(實(shí)線,單位:m2·s-2),位勢(shì)高度(點(diǎn)點(diǎn)虛線,單位:gpm)及位渦(陰影區(qū)表示PV大于0.2的區(qū)域,單位:PVU)的經(jīng)度-時(shí)間剖面

    圖3 是500hPa面上位勢(shì)高度與加速度遷移項(xiàng)位勢(shì)及位渦在30°N—35°N之間平均的經(jīng)度-時(shí)間演變圖,從圖中可見,從7月20日08∶00開始,隨著低渦從高原中東部開始東移直至移出高原的過(guò)程中,實(shí)線表示的加速度遷移項(xiàng)位勢(shì)從高原到下游地區(qū)隨時(shí)間呈現(xiàn)出一條明顯的正值中心區(qū),與陰影區(qū)表示的位渦的正異常區(qū)域是一致的,說(shuō)明加速度遷移項(xiàng)位勢(shì)確實(shí)可以用以示蹤低渦的東移演變.不同的是,加速度遷移項(xiàng)位勢(shì)表示的輻合運(yùn)動(dòng)在高原上時(shí)較強(qiáng),而位勢(shì)高度及位渦則是隨后加強(qiáng)的,在時(shí)間上稍微滯后于加速度遷移項(xiàng)位勢(shì)的加強(qiáng).當(dāng)?shù)蜏u東移出高原并入低槽中加強(qiáng)后,加速度遷移項(xiàng)位勢(shì)是減弱的,此時(shí)位勢(shì)高度減小,位渦正異常也加大,表明低槽加深加強(qiáng).可以推測(cè),這個(gè)低槽加強(qiáng)可能與前期加速度遷移項(xiàng)位勢(shì)的輻合有關(guān),這與天氣學(xué)上用風(fēng)場(chǎng)輻合可導(dǎo)致低渦加強(qiáng)的動(dòng)力過(guò)程的描述也是一致的,也進(jìn)一步說(shuō)明加速度遷移項(xiàng)位勢(shì)所指示的輻合輻散對(duì)系統(tǒng)的強(qiáng)度變化有一定的動(dòng)力指示意義.

    從以上對(duì)東移高原低渦的分析可見,由于加速度遷移項(xiàng)位勢(shì)包含了水平風(fēng)場(chǎng)平流的輻合輻散動(dòng)力特征,因此,從加速度遷移項(xiàng)位勢(shì)的分布,我們就可以直接判斷低渦系統(tǒng)的水平風(fēng)場(chǎng)平流的輻合輻散動(dòng)力效應(yīng),而不用再分析其散度的分布.而且,隨著空間分辨率越來(lái)越高,計(jì)算的散度分布會(huì)越來(lái)越亂,而加速度遷移項(xiàng)位勢(shì)的分布相對(duì)于散度分布來(lái)說(shuō)要光滑很多.因此,作為一個(gè)新的動(dòng)力變量,加速度遷移項(xiàng)位勢(shì)的引入可用于示蹤東移高原低渦系統(tǒng)以及指示低渦系統(tǒng)的平流的輻合輻散效應(yīng).

    以上是加速度遷移項(xiàng)位勢(shì)在東移高原低渦系統(tǒng)分析中的應(yīng)用,以下分析其在登陸臺(tái)風(fēng)中的應(yīng)用.

    3.2 在登陸臺(tái)風(fēng)Bilis中的應(yīng)用

    登陸臺(tái)風(fēng)是造成我國(guó)暴雨的又一個(gè)主要系統(tǒng),2006年的4號(hào)臺(tái)風(fēng)Bilis在整個(gè)發(fā)展過(guò)程中,結(jié)構(gòu)松散,中心附近對(duì)流不強(qiáng),沒有眼壁和風(fēng)眼結(jié)構(gòu),沒有統(tǒng)一的臺(tái)風(fēng)中心,風(fēng)力不是特別強(qiáng),但含水量大,造成南方地區(qū)大面積洪澇,以及山體滑坡和泥石流等次生災(zāi)害,是近十年來(lái)臺(tái)風(fēng)災(zāi)害導(dǎo)致傷亡人數(shù)最多的一次,本小節(jié)以Bilis為例進(jìn)行加速度遷移項(xiàng)位勢(shì)的分析.所用資料為日本氣象廳區(qū)域譜模式(RSM)再分析資料,空間分辨率為20km,時(shí)間分辨率為6h.

    圖4 850hPa上加速度遷移項(xiàng)位勢(shì)(實(shí)線,單位:m2·s-2),無(wú)旋風(fēng)(矢量,單位:m·s)及散度(陰影區(qū),單位:1×10-4·s-1)分布(a)2006-07-12-20∶00;(b)2006-07-13-08∶00;(c)2006-07-14-02∶00;(d)2006-07-14-20∶00.圖中粗虛線為Bilis臺(tái)風(fēng)路徑,臺(tái)風(fēng)符號(hào)表示當(dāng)前時(shí)刻臺(tái)風(fēng)中心位置.

    從圖4可見,無(wú)論是Bilis臺(tái)風(fēng)登陸前(圖4a,b,c)還是登陸后(圖4d),與位勢(shì)高度類似(圖略),實(shí)線所示的加速度遷移項(xiàng)位勢(shì)都能較好地指示出臺(tái)風(fēng)系統(tǒng),加速度遷移項(xiàng)位勢(shì)的正值中心與臺(tái)風(fēng)中心位置比較接近,用加速度遷移項(xiàng)位勢(shì)的正異常中心可以大概判斷出臺(tái)風(fēng)中心位置.此外,從加速度遷移項(xiàng)位勢(shì)與水平風(fēng)散度的分布來(lái)看,加速度遷移項(xiàng)位勢(shì)等值線相對(duì)密集的區(qū)域,也就是加速度遷移項(xiàng)位勢(shì)梯度較大的區(qū)域,往往也是散度較大的區(qū)域(這個(gè)特征在圖4c中表現(xiàn)最顯著,環(huán)狀深色陰影區(qū)代表的輻合區(qū)與加速度遷移項(xiàng)位勢(shì)等值線分布較密集的區(qū)域基本是一致的).在加速度遷移項(xiàng)等值線曲率變化大的地方(這種區(qū)域加速度遷移項(xiàng)位勢(shì)的拉普拉斯算子最大),一般也對(duì)應(yīng)著散度大值區(qū),如圖4a中,在Bilis的東南部(122°E,21°N)附近,加速度遷移項(xiàng)位勢(shì)的400和600的等值線有明顯的彎曲,對(duì)應(yīng)著明顯的輻合發(fā)生(圖中陰影區(qū));圖4b中,臺(tái)風(fēng)西側(cè)的加速度遷移項(xiàng)位勢(shì)有小范圍彎曲的地方,也是輻合較強(qiáng)的區(qū)域;又如圖4c中,在臺(tái)風(fēng)中心南側(cè),從外到內(nèi)都有加速度遷移項(xiàng)位勢(shì)等值線發(fā)生明顯彎曲,也對(duì)應(yīng)著無(wú)旋風(fēng)的輻合大值區(qū)及散度大值區(qū).在圖4d中,由于臺(tái)風(fēng)登陸后,受復(fù)雜下墊面條件影響,臺(tái)風(fēng)結(jié)構(gòu)不再像洋面上時(shí)那樣規(guī)則,但是也可分析得到加速度遷移項(xiàng)位勢(shì)的等值線曲率大的區(qū)域,也基本對(duì)應(yīng)著較強(qiáng)的輻合運(yùn)動(dòng).圖4中陰影區(qū)只表示了輻合發(fā)生的區(qū)域,如果全標(biāo)出輻合輻散運(yùn)動(dòng)會(huì)發(fā)現(xiàn),整個(gè)散度場(chǎng)的分布比較凌亂,對(duì)分析臺(tái)風(fēng)動(dòng)力結(jié)構(gòu)并不清楚.但是從加速度遷移項(xiàng)位勢(shì)分布看,只要注意加速度遷移項(xiàng)位勢(shì)正值中心區(qū)及其等值線曲率變化較大的區(qū)域,就基本可以判斷出大的輻合區(qū)域.從加速度遷移項(xiàng)位勢(shì)的數(shù)值看,在洋面上時(shí),其數(shù)值較大,中心等值線為1400,最強(qiáng)時(shí)達(dá)到1600(圖略),隨著臺(tái)風(fēng)的移動(dòng)及登陸,數(shù)值從1600逐漸降低到300多,其數(shù)值演變趨勢(shì)與臺(tái)風(fēng)強(qiáng)度的變化也是一致的.因此,加速度遷移項(xiàng)位勢(shì)的強(qiáng)弱也可表示出臺(tái)風(fēng)系統(tǒng)的強(qiáng)度變化.

    從本節(jié)加速度遷移項(xiàng)位勢(shì)對(duì)Bilis臺(tái)風(fēng)的分析可見,它可指示臺(tái)風(fēng)系統(tǒng)位置及輻合大值區(qū)的位置以及臺(tái)風(fēng)強(qiáng)度的演變趨勢(shì).可見,加速度遷移項(xiàng)位勢(shì)可用于指示臺(tái)風(fēng)系統(tǒng)以及顯示臺(tái)風(fēng)系統(tǒng)里的輻合輻散運(yùn)動(dòng)特征,還可反映臺(tái)風(fēng)系統(tǒng)強(qiáng)度,可作為一個(gè)新的動(dòng)力診斷變量來(lái)示蹤臺(tái)風(fēng)的強(qiáng)度變化及其輻合輻散效應(yīng).

    4 結(jié) 論

    在高守亭等[6]對(duì)加速度遷移項(xiàng)散度作用研究的基礎(chǔ)上,即:加速度遷移項(xiàng)能代表氣塊受到的所有力的輻合或輻散效應(yīng),通過(guò)計(jì)算加速度遷移散度項(xiàng)就可以判斷氣塊受到的力的輻合或輻散情況,本文對(duì)加速度遷移項(xiàng)進(jìn)一步進(jìn)行調(diào)和-余弦的分解,并引入加速度遷移項(xiàng)位勢(shì)的概念,用加速度遷移項(xiàng)對(duì)應(yīng)的位勢(shì)部分(可反映平流的輻合輻散)對(duì)東移高原低渦及登陸臺(tái)風(fēng)Bilis的演變過(guò)程進(jìn)行了分析.結(jié)果表明加速度遷移項(xiàng)位勢(shì)對(duì)低渦系統(tǒng)有較好的描述作用,用于指示東移的高原低渦系統(tǒng)較常用的500hPa位勢(shì)高度場(chǎng)對(duì)低渦的分析更為清晰.此外,加速度遷移項(xiàng)位勢(shì)用于登陸臺(tái)風(fēng)Bilis的分析也表明,它可指示臺(tái)風(fēng)環(huán)流系統(tǒng),還可用以判斷洋面上臺(tái)風(fēng)的中心位置以及臺(tái)風(fēng)強(qiáng)度變化.由于加速度遷移項(xiàng)位勢(shì)可反映出水平平流的輻合輻散,可以作為一個(gè)新的動(dòng)力變量來(lái)診斷并示蹤低渦及臺(tái)風(fēng)等天氣系統(tǒng)的演變.

    (References)

    [1]周明煜,錢粉蘭,陳陟等.西藏高原斜壓對(duì)流邊界層風(fēng)、溫、濕廓線特征.地球物理學(xué)報(bào),2002,45(6):773-783.Zhou M Y,Qian F L,Chen Z,et al.The characteristics of the profiles for wind,temperature and humidity in the baroclinic convective boundary layer on Xizang plateau.Chinese J.Geophys.(in Chinese),2002,45(6):773-783.

    [2]姚秀萍,吳國(guó)雄,劉屹岷等.熱帶對(duì)流層上空東風(fēng)帶擾動(dòng)影響西太平洋副熱帶高壓的個(gè)例分析.氣象學(xué)報(bào),2007,65(2):198-207.Yao X P,Wu G X,Liu Y M,et al.Case study on the impact of the vortex in the easterlies over the tropical upper troposphere on the subtropical anticyclone over western Pacific Ocean.Acta Meteorologica Sinica (in Chinese),2007,65(2):198-207.

    [3]Zhao Y S,Liu L P.Analyses of wind structure of typhoon fung-wong(2008)and its relation to precipitation region.Journal of Tropical Meteorology,2009,15(2):192-203.

    [4]周玉淑,冉令坤.平流渦度方程及其在2006年Bilis臺(tái)風(fēng)分析中的應(yīng)用.物理學(xué)報(bào),2010,59(2):1367-1377.Zhou Y S,Ran L K.Advective vorticity equation and its application to the vorticity variation of typhoon Bills in 2006.Acta Physica Sinica (in Chinese),2010,59(2):1367-1377.

    [5]楊平章,賀海晏.非線性水平動(dòng)量平流作用對(duì)熱帶氣旋運(yùn)動(dòng)的影響.中山大學(xué)學(xué)報(bào) (自然科學(xué)版),1995,34(1):82-89.Yang P Z,He H Y.The effect of non-linear horizontal momentum advection on tropical cyclone motion.Acta Scientiarum Naturalium Universitatis Sunyatseni (in Chinese),1995,34(1):82-89.

    [6]高守亭,鄧滌菲,周玉淑.加速度遷移項(xiàng)散度在渦旋系統(tǒng)動(dòng)力識(shí)別上的初步應(yīng)用.地球物理學(xué)報(bào),2011,54(4):875-885.Gao S T,Deng D F,Zhou Y S.The effect of the divergence of relocation term in acceleration and its application to the dynamic identification of vortex system.Chinese J.Geophys.(in Chinese),2011,54(4):875-885.

    [7]Chen Q S,Kuo Y H.A consistency condition for wind-field reconstruction in a limited area and a harmonic-cosine series expansion.Mon.Wea.Rev.,1992,120(11):2653-2670.

    [8]周玉淑,曹潔,高守亭.有限區(qū)域風(fēng)場(chǎng)分解方法及其在臺(tái)風(fēng)SAOMEI研究中的應(yīng)用.物理學(xué)報(bào),2008,57(10):6654-6665.Zhao Y S,Cao J,Gao S T.The method of decomposing wind field in a limited area and its application to typhoon SAOMEI.Acta Physica Sinica (in Chinese),2008,57(10):6654-6665.

    [9]周玉淑,曹潔.有限區(qū)域風(fēng)場(chǎng)的分解和重建.物理學(xué)報(bào),2010,59(4):2898-2906.Zhou Y S,Cao J.Partitioning and reconstruction problem of the wind in a limited region.Acta Physica Sinica(in Chinese),2010,59(4):2898-2906.

    [10]高文良,郁淑華.高原低渦東移出高原的平均環(huán)流場(chǎng)分析.高原氣象,2007,26(1):206-212.Gao W L,Yu S H.Analyses on mean circulation field of the plateau low vortex moving out of Tibetan Plateau.Plateau Meteorology (in Chinese),2007,26(1):206-212.

    [11]郁淑華,高文良.高原低渦移出高原的觀測(cè)事實(shí)分析.氣象學(xué)報(bào),2006,64(3):392-399.Yu S H,Gao W L.Observational analysis on the movement of vortices before/after moving out the Tibetan Plateau.Acta Meteorologica Sinica (in Chinese),2006,64(3):392-399.

    [12]郁淑華,高文良.青藏高原低渦移出高原的大尺度條件.高原氣象,2008,27(6):1276-1287.Yu S H,Gao W L.The Large-Scale conditions of the vortex moving out of Qinghai-Xizang plateau.Plateau Meteorology(in Chinese),2008,27(6):1276-1287.

    [13]何光碧,高文良,屠妮妮.兩次高原低渦東移特征及發(fā)展機(jī)制動(dòng)力診斷.氣象學(xué)報(bào),2009,67(4):599-612.He G B,Gao W L,Tu N N.The dynamic diagnosis on easterwards moving characteristics and developing mechanism of two Tibetan Plateau vortex processes.Acta Meteorologica Sinica(in Chinese),2009,67(4):599-612.

    猜你喜歡
    東移位勢(shì)平流
    2007—2020年威海機(jī)場(chǎng)4—9月平流低云的變化特征及影響因子分析
    含Hardy位勢(shì)的非線性Schr?dinger-Poisson方程正規(guī)化解的多重性
    一類帶強(qiáng)制位勢(shì)的p-Laplace特征值問(wèn)題
    一次東移型西南低渦引發(fā)的強(qiáng)降水診斷分析
    基于衛(wèi)星和Argo觀測(cè)的阿拉伯海中北部海表鹽度季節(jié)和年際變化
    含變號(hào)位勢(shì)的ρ-Kirchhoff型方程組無(wú)窮多個(gè)高能量解的存在性
    含位勢(shì)的非線性雙調(diào)和方程解的存在性
    荊州市一次局地濃霧天氣特征分析
    浦東機(jī)場(chǎng)一次低云低能見度天氣氣象服務(wù)總結(jié)
    在玩中“悟”出逆推策略
    ——數(shù)學(xué)實(shí)驗(yàn)課《巧移汽車》教學(xué)設(shè)計(jì)
    www.av在线官网国产| 亚洲精品国产av成人精品| 国产一区二区三区av在线| 精品少妇黑人巨大在线播放| 亚洲精品乱码久久久v下载方式| 国产高清有码在线观看视频| 久久久国产一区二区| 欧美日韩精品成人综合77777| 欧美日韩综合久久久久久| 亚洲三级黄色毛片| 日本-黄色视频高清免费观看| 国产av在哪里看| 一级a做视频免费观看| 免费看美女性在线毛片视频| 最近2019中文字幕mv第一页| 国产亚洲av片在线观看秒播厂 | 日韩欧美一区视频在线观看 | 久久99精品国语久久久| 我要看日韩黄色一级片| 六月丁香七月| 又爽又黄a免费视频| 国产欧美日韩精品一区二区| 两个人的视频大全免费| 国产精品av视频在线免费观看| 搡老妇女老女人老熟妇| 午夜激情欧美在线| 国产伦精品一区二区三区四那| 免费观看在线日韩| 国产av不卡久久| 欧美成人a在线观看| 国产精品久久久久久久久免| 日韩av在线大香蕉| 色综合色国产| kizo精华| 看十八女毛片水多多多| www.av在线官网国产| 黄色一级大片看看| 欧美高清成人免费视频www| 久久久欧美国产精品| 日韩,欧美,国产一区二区三区| 色综合色国产| 亚洲人成网站高清观看| 18禁动态无遮挡网站| 观看美女的网站| 亚洲内射少妇av| 91在线精品国自产拍蜜月| 99久久中文字幕三级久久日本| 亚洲av中文av极速乱| 肉色欧美久久久久久久蜜桃 | 国产白丝娇喘喷水9色精品| 国产乱来视频区| 日本免费a在线| 一级片'在线观看视频| 国产高潮美女av| 欧美+日韩+精品| 日韩电影二区| 亚洲欧美日韩东京热| 亚洲不卡免费看| 51国产日韩欧美| 校园人妻丝袜中文字幕| 久久99热这里只频精品6学生| 97超视频在线观看视频| 国产乱人视频| a级一级毛片免费在线观看| 久久久精品免费免费高清| 免费黄频网站在线观看国产| 超碰97精品在线观看| 男人狂女人下面高潮的视频| 欧美3d第一页| 亚洲欧美中文字幕日韩二区| 在线播放无遮挡| 狠狠精品人妻久久久久久综合| 中文字幕亚洲精品专区| 一区二区三区乱码不卡18| 中文字幕制服av| 一级毛片我不卡| 国产亚洲精品久久久com| 女人被狂操c到高潮| 真实男女啪啪啪动态图| 精品少妇黑人巨大在线播放| 欧美区成人在线视频| 在线天堂最新版资源| 精华霜和精华液先用哪个| 99久久精品热视频| 国产精品日韩av在线免费观看| 中文字幕av成人在线电影| 少妇高潮的动态图| 国产一级毛片在线| 日日干狠狠操夜夜爽| 久久久久精品久久久久真实原创| 国产黄色免费在线视频| 九九在线视频观看精品| 欧美xxⅹ黑人| 纵有疾风起免费观看全集完整版 | 亚洲av电影不卡..在线观看| 大香蕉97超碰在线| 国产在线男女| 99热这里只有是精品在线观看| 最后的刺客免费高清国语| 18禁裸乳无遮挡免费网站照片| 亚洲综合精品二区| 日韩成人av中文字幕在线观看| 欧美xxxx性猛交bbbb| 国产精品伦人一区二区| 亚洲av不卡在线观看| 亚洲天堂国产精品一区在线| 色5月婷婷丁香| 尤物成人国产欧美一区二区三区| 一级黄片播放器| 在线观看免费高清a一片| 超碰av人人做人人爽久久| 国产乱人视频| 欧美三级亚洲精品| 亚洲经典国产精华液单| 伦理电影大哥的女人| 中文字幕av成人在线电影| 成人综合一区亚洲| 99久国产av精品| 51国产日韩欧美| 午夜老司机福利剧场| 日本色播在线视频| 亚洲aⅴ乱码一区二区在线播放| kizo精华| 久久久久久久大尺度免费视频| 亚洲人与动物交配视频| 久久久久久久国产电影| 久久精品国产亚洲av涩爱| 女人久久www免费人成看片| 永久免费av网站大全| 黑人高潮一二区| 日日摸夜夜添夜夜爱| 久久草成人影院| 亚洲av在线观看美女高潮| 亚洲一区高清亚洲精品| 亚洲欧美成人精品一区二区| 网址你懂的国产日韩在线| 国产毛片a区久久久久| 亚洲国产日韩欧美精品在线观看| 色视频www国产| 日韩国内少妇激情av| videos熟女内射| 能在线免费观看的黄片| 18禁在线无遮挡免费观看视频| 久久99热这里只频精品6学生| 免费看美女性在线毛片视频| 舔av片在线| 免费观看无遮挡的男女| 久久6这里有精品| 国产伦精品一区二区三区视频9| 日韩一本色道免费dvd| 免费高清在线观看视频在线观看| 国产高清不卡午夜福利| 网址你懂的国产日韩在线| 麻豆久久精品国产亚洲av| 1000部很黄的大片| 欧美高清性xxxxhd video| 日日撸夜夜添| 人妻夜夜爽99麻豆av| 亚洲精华国产精华液的使用体验| 国产 一区 欧美 日韩| 免费黄网站久久成人精品| 亚洲人成网站在线播| videossex国产| 美女cb高潮喷水在线观看| 午夜免费激情av| 国产老妇伦熟女老妇高清| 亚洲精品乱码久久久v下载方式| 如何舔出高潮| 中文字幕av成人在线电影| 成人毛片a级毛片在线播放| 蜜臀久久99精品久久宅男| 97人妻精品一区二区三区麻豆| 日韩不卡一区二区三区视频在线| 亚洲av免费在线观看| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久精品电影| 国产三级在线视频| 一级毛片黄色毛片免费观看视频| 日韩欧美国产在线观看| 91狼人影院| 亚洲精品一二三| 国产精品人妻久久久久久| 一级毛片aaaaaa免费看小| 亚洲精华国产精华液的使用体验| 日韩,欧美,国产一区二区三区| 国产精品一区二区在线观看99 | 免费av毛片视频| 日本-黄色视频高清免费观看| 日本av手机在线免费观看| 爱豆传媒免费全集在线观看| 99热网站在线观看| 丝袜美腿在线中文| 国产亚洲av嫩草精品影院| 久久久久性生活片| 国产免费又黄又爽又色| 不卡视频在线观看欧美| 高清视频免费观看一区二区 | 国产永久视频网站| 国产熟女欧美一区二区| 日韩人妻高清精品专区| 免费播放大片免费观看视频在线观看| 99久久中文字幕三级久久日本| 两个人视频免费观看高清| 久久久久九九精品影院| 久久久色成人| 中文字幕人妻熟人妻熟丝袜美| 国产成人午夜福利电影在线观看| 亚洲欧美精品自产自拍| 亚洲18禁久久av| 蜜桃亚洲精品一区二区三区| 国产精品蜜桃在线观看| 日韩成人伦理影院| 国产探花极品一区二区| 夜夜看夜夜爽夜夜摸| 美女国产视频在线观看| 日韩伦理黄色片| 搡女人真爽免费视频火全软件| 亚洲怡红院男人天堂| 国产 亚洲一区二区三区 | 国产亚洲5aaaaa淫片| 久久国产乱子免费精品| 久久97久久精品| 国产黄色视频一区二区在线观看| av黄色大香蕉| 免费观看av网站的网址| 国产精品国产三级专区第一集| 色综合亚洲欧美另类图片| 少妇的逼好多水| 舔av片在线| 国产精品久久久久久av不卡| 亚洲精品第二区| 久久久久久久亚洲中文字幕| 亚洲国产日韩欧美精品在线观看| 内射极品少妇av片p| 亚洲成色77777| 成人毛片60女人毛片免费| 一区二区三区高清视频在线| 精品人妻视频免费看| 婷婷色麻豆天堂久久| a级毛片免费高清观看在线播放| 亚洲欧美清纯卡通| 国产人妻一区二区三区在| 18+在线观看网站| 免费观看av网站的网址| 最近最新中文字幕免费大全7| 亚洲av一区综合| 麻豆成人午夜福利视频| 中文字幕制服av| 免费av不卡在线播放| 国内少妇人妻偷人精品xxx网站| 亚洲国产高清在线一区二区三| 两个人的视频大全免费| 蜜臀久久99精品久久宅男| 2018国产大陆天天弄谢| 婷婷色av中文字幕| 久久久久久久久久久丰满| 日韩成人伦理影院| 赤兔流量卡办理| 网址你懂的国产日韩在线| 婷婷六月久久综合丁香| 久久久精品94久久精品| 丰满乱子伦码专区| 狂野欧美激情性xxxx在线观看| 欧美日韩视频高清一区二区三区二| 亚洲高清免费不卡视频| 国产91av在线免费观看| 九九爱精品视频在线观看| 女的被弄到高潮叫床怎么办| 亚洲天堂国产精品一区在线| 天堂网av新在线| 男女那种视频在线观看| 一级毛片aaaaaa免费看小| 亚洲电影在线观看av| 亚洲国产高清在线一区二区三| av在线蜜桃| 纵有疾风起免费观看全集完整版 | 国产精品久久久久久久电影| 亚洲精品亚洲一区二区| 国产真实伦视频高清在线观看| 五月伊人婷婷丁香| 国产伦精品一区二区三区视频9| 我的女老师完整版在线观看| 草草在线视频免费看| 亚洲欧美成人综合另类久久久| 亚洲一级一片aⅴ在线观看| 亚州av有码| 免费播放大片免费观看视频在线观看| 乱人视频在线观看| 男女视频在线观看网站免费| 国产av国产精品国产| 国产高潮美女av| 久久久午夜欧美精品| 舔av片在线| 久久精品夜色国产| 亚洲人与动物交配视频| 久久久久久久久中文| av在线蜜桃| 九九久久精品国产亚洲av麻豆| 床上黄色一级片| 久久99精品国语久久久| 亚洲图色成人| 成人午夜精彩视频在线观看| 欧美三级亚洲精品| 欧美另类一区| videossex国产| 男人狂女人下面高潮的视频| av福利片在线观看| 汤姆久久久久久久影院中文字幕 | 欧美日本视频| 不卡视频在线观看欧美| 精品酒店卫生间| 国产熟女欧美一区二区| 2022亚洲国产成人精品| 男女国产视频网站| 亚洲国产最新在线播放| 国产精品嫩草影院av在线观看| 成年免费大片在线观看| 久久人人爽人人爽人人片va| 成人av在线播放网站| 久久这里有精品视频免费| 午夜日本视频在线| 乱人视频在线观看| 欧美成人午夜免费资源| 国产免费福利视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品亚洲一区二区| 国产精品久久久久久av不卡| 精品午夜福利在线看| 草草在线视频免费看| 国产片特级美女逼逼视频| 成年女人在线观看亚洲视频 | 特级一级黄色大片| 1000部很黄的大片| 国产av在哪里看| 啦啦啦中文免费视频观看日本| 日本黄色片子视频| 啦啦啦韩国在线观看视频| 久久99精品国语久久久| 国产男人的电影天堂91| 国模一区二区三区四区视频| 国产探花在线观看一区二区| 麻豆久久精品国产亚洲av| 日韩强制内射视频| 能在线免费看毛片的网站| 最近2019中文字幕mv第一页| 国产视频首页在线观看| 免费黄色在线免费观看| av黄色大香蕉| 精品国产露脸久久av麻豆 | 校园人妻丝袜中文字幕| 国产成人a区在线观看| 色吧在线观看| 女人被狂操c到高潮| 国产高潮美女av| 亚洲欧美日韩卡通动漫| 老师上课跳d突然被开到最大视频| 国产白丝娇喘喷水9色精品| 美女黄网站色视频| 综合色av麻豆| 少妇高潮的动态图| 美女被艹到高潮喷水动态| 嘟嘟电影网在线观看| 在线播放无遮挡| 天天躁夜夜躁狠狠久久av| 欧美精品国产亚洲| 久久精品夜色国产| 一个人看的www免费观看视频| 久久精品人妻少妇| 国产 一区精品| 2018国产大陆天天弄谢| 在线免费十八禁| 久久久久久久大尺度免费视频| 亚洲三级黄色毛片| 亚洲第一区二区三区不卡| 五月玫瑰六月丁香| 特大巨黑吊av在线直播| 亚洲经典国产精华液单| 亚洲电影在线观看av| 成人特级av手机在线观看| or卡值多少钱| 少妇高潮的动态图| 久久久欧美国产精品| 欧美精品一区二区大全| 国产亚洲最大av| 禁无遮挡网站| 午夜久久久久精精品| 最近中文字幕2019免费版| 97人妻精品一区二区三区麻豆| av又黄又爽大尺度在线免费看| 久久久色成人| 成人亚洲精品av一区二区| 国产成人freesex在线| 欧美精品一区二区大全| 国产精品福利在线免费观看| 国产成人精品婷婷| 天美传媒精品一区二区| 精品一区二区三卡| 免费看不卡的av| 国产极品天堂在线| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 纵有疾风起免费观看全集完整版 | 亚洲自偷自拍三级| 最近手机中文字幕大全| 国产成人精品一,二区| .国产精品久久| 九九久久精品国产亚洲av麻豆| 免费看av在线观看网站| 亚洲经典国产精华液单| 国产精品人妻久久久影院| 国产老妇伦熟女老妇高清| 国产 一区 欧美 日韩| 午夜福利视频精品| 亚洲熟女精品中文字幕| 久久这里只有精品中国| 久久99蜜桃精品久久| 免费av毛片视频| or卡值多少钱| 女人十人毛片免费观看3o分钟| 国产精品蜜桃在线观看| 欧美+日韩+精品| 日本黄大片高清| 国产成人免费观看mmmm| 欧美激情在线99| 欧美成人午夜免费资源| 精品少妇黑人巨大在线播放| av又黄又爽大尺度在线免费看| 毛片女人毛片| 国产精品麻豆人妻色哟哟久久 | 毛片一级片免费看久久久久| 伦理电影大哥的女人| 久久久久久久久久黄片| 水蜜桃什么品种好| 国产欧美日韩精品一区二区| 在线观看av片永久免费下载| 乱人视频在线观看| 亚洲不卡免费看| 欧美最新免费一区二区三区| www.色视频.com| 禁无遮挡网站| 欧美成人午夜免费资源| 午夜福利成人在线免费观看| 成人特级av手机在线观看| 九九久久精品国产亚洲av麻豆| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 国产精品综合久久久久久久免费| 免费观看性生交大片5| 夜夜爽夜夜爽视频| 日日啪夜夜撸| 高清日韩中文字幕在线| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 日韩,欧美,国产一区二区三区| 久久久久久国产a免费观看| 五月天丁香电影| videos熟女内射| 免费看光身美女| 中文字幕人妻熟人妻熟丝袜美| 国产69精品久久久久777片| 内地一区二区视频在线| 中文字幕av在线有码专区| 免费看不卡的av| 亚洲自拍偷在线| 久久99热这里只频精品6学生| 国产美女午夜福利| 久久这里有精品视频免费| 深夜a级毛片| 国产又色又爽无遮挡免| 亚洲国产日韩欧美精品在线观看| 在线观看美女被高潮喷水网站| 丝瓜视频免费看黄片| 亚洲av中文字字幕乱码综合| 日产精品乱码卡一卡2卡三| 我要看日韩黄色一级片| 中文字幕制服av| 亚洲av国产av综合av卡| 日本免费a在线| 免费看av在线观看网站| 日韩av不卡免费在线播放| 麻豆国产97在线/欧美| 亚洲av在线观看美女高潮| 女的被弄到高潮叫床怎么办| 男女那种视频在线观看| 国产高潮美女av| 五月玫瑰六月丁香| 内地一区二区视频在线| 91午夜精品亚洲一区二区三区| 亚洲va在线va天堂va国产| kizo精华| 免费av毛片视频| 极品教师在线视频| 99视频精品全部免费 在线| 最近中文字幕高清免费大全6| 午夜福利成人在线免费观看| 久久久久久久久久人人人人人人| 夫妻性生交免费视频一级片| 欧美日韩在线观看h| 亚洲aⅴ乱码一区二区在线播放| 成年av动漫网址| 日本黄大片高清| 激情五月婷婷亚洲| 久久精品久久久久久久性| 秋霞伦理黄片| 中文在线观看免费www的网站| 少妇的逼水好多| 亚洲自偷自拍三级| 成人鲁丝片一二三区免费| 一个人看的www免费观看视频| 午夜福利视频精品| 搞女人的毛片| 亚洲欧美日韩卡通动漫| 国产伦精品一区二区三区四那| 日本猛色少妇xxxxx猛交久久| 成人二区视频| 国产伦一二天堂av在线观看| 国产欧美另类精品又又久久亚洲欧美| 人妻系列 视频| 我的老师免费观看完整版| 卡戴珊不雅视频在线播放| 久久人人爽人人片av| 永久网站在线| 国产视频首页在线观看| 日韩一区二区三区影片| 国产精品一区二区性色av| 深爱激情五月婷婷| 老女人水多毛片| 精品一区在线观看国产| 波多野结衣巨乳人妻| 久久97久久精品| 啦啦啦韩国在线观看视频| 亚洲国产欧美人成| 18禁动态无遮挡网站| 日日撸夜夜添| 国产成人a∨麻豆精品| 岛国毛片在线播放| 日本猛色少妇xxxxx猛交久久| 国产久久久一区二区三区| 亚洲第一区二区三区不卡| 99久久九九国产精品国产免费| 久久久久久九九精品二区国产| 一夜夜www| 国产在视频线在精品| freevideosex欧美| 日本一二三区视频观看| 永久网站在线| 插阴视频在线观看视频| .国产精品久久| 亚洲av成人av| 伊人久久国产一区二区| 午夜激情福利司机影院| 亚洲丝袜综合中文字幕| av免费观看日本| 欧美激情在线99| 一级爰片在线观看| 在线观看av片永久免费下载| 亚洲欧美日韩卡通动漫| 久久这里有精品视频免费| 色综合亚洲欧美另类图片| 毛片一级片免费看久久久久| 简卡轻食公司| 少妇熟女aⅴ在线视频| 国产亚洲91精品色在线| 少妇熟女欧美另类| 午夜福利高清视频| 日日摸夜夜添夜夜爱| 久久久欧美国产精品| 草草在线视频免费看| 欧美成人午夜免费资源| 国产黄色免费在线视频| 99久国产av精品国产电影| 欧美变态另类bdsm刘玥| 国产成人精品婷婷| 国产精品.久久久| 午夜福利在线观看吧| 免费看a级黄色片| 国产单亲对白刺激| 久久精品人妻少妇| 亚洲成人一二三区av| 久久久亚洲精品成人影院| 又爽又黄a免费视频| 国产一区亚洲一区在线观看| 特级一级黄色大片| 久久97久久精品| 男女啪啪激烈高潮av片| 国产中年淑女户外野战色| 美女黄网站色视频| 乱码一卡2卡4卡精品| 日韩强制内射视频| 国产又色又爽无遮挡免| 午夜福利视频1000在线观看| 国产真实伦视频高清在线观看| 亚洲精品国产av蜜桃| 国国产精品蜜臀av免费| 国产又色又爽无遮挡免| 大香蕉97超碰在线| 亚洲欧美一区二区三区黑人 | 亚洲欧美精品专区久久| 精品亚洲乱码少妇综合久久| 少妇被粗大猛烈的视频| 免费看av在线观看网站| 久久99热这里只有精品18| 午夜精品一区二区三区免费看| 国产视频首页在线观看| 小蜜桃在线观看免费完整版高清| 乱人视频在线观看| av在线蜜桃| 99久国产av精品国产电影| 神马国产精品三级电影在线观看| 非洲黑人性xxxx精品又粗又长| 99久久中文字幕三级久久日本| 国产成人精品一,二区| 男女视频在线观看网站免费| 色网站视频免费| 蜜桃亚洲精品一区二区三区|