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

    初始外圍尺度與加熱對臺風(fēng)次眼墻形成的影響

    2022-10-18 10:17:36林超然談?wù)苊?/span>
    氣象科學(xué) 2022年4期
    關(guān)鍵詞:雨帶渦度極大值

    林超然 談?wù)苊?/p>

    (南京大學(xué) 中尺度災(zāi)害性天氣教育部重點(diǎn)實(shí)驗(yàn)室/大氣科學(xué)學(xué)院,南京 210023)

    引 言

    臺風(fēng)發(fā)展過程伴隨著風(fēng)場、對流結(jié)構(gòu)變化[1-4],約50%的強(qiáng)臺風(fēng)會出現(xiàn)雙眼墻結(jié)構(gòu)、眼墻替換[5],引起臺風(fēng)強(qiáng)度的劇烈變化。臺風(fēng)雙眼墻具有兩個重要的結(jié)構(gòu)特征:一是典型的對流結(jié)構(gòu)特征,出現(xiàn)構(gòu)成主、次眼墻的對流環(huán)和二者之間以下沉運(yùn)動為主的moat區(qū)域;二是典型的風(fēng)場結(jié)構(gòu)特征:存在與主眼墻相對應(yīng)的切向風(fēng)主極大值和與次眼墻相對應(yīng)的低層切向風(fēng)次極大值。

    主眼墻和外圍雨帶的非絕熱加熱對臺風(fēng)強(qiáng)度、結(jié)構(gòu)變化至關(guān)重要。平衡動力學(xué)理論中,非絕熱加熱使臺風(fēng)在主環(huán)流上形成一個對流層低層徑向流入、眼墻處上升、對流層高層流出的次級環(huán)流結(jié)構(gòu)[6-8],其中低層入流向內(nèi)輸送絕對角動量,由于角動量在邊界層無摩擦?xí)r守恒,從而促使切向風(fēng)增強(qiáng)[9-10]。內(nèi)雨帶通常指3倍最大風(fēng)速半徑(Radius of Maximum Wind, RMW)以內(nèi)的雨帶,外雨帶指3倍RMW外的雨帶[4],雨帶分布不同導(dǎo)致的加熱率分布差異會對臺風(fēng)結(jié)構(gòu)和強(qiáng)度產(chǎn)生很大影響。

    同樣,非絕熱加熱對次眼墻發(fā)生、發(fā)展有著重要影響。WANG[4]指出增強(qiáng)外雨帶的非絕熱加熱有利于臺風(fēng)風(fēng)場擴(kuò)張和次眼墻的形成。Moon, et al[11]指出外雨帶對流和層云降水可以引發(fā)次切向風(fēng)極值出現(xiàn)。Rozoff, et al[12]指出臺風(fēng)外圍對流場加熱與風(fēng)場間存在正反饋?zhàn)饔?,即主眼墻增?qiáng)及外擴(kuò)風(fēng)場使慣性穩(wěn)定度增加,非絕熱加熱向動能的轉(zhuǎn)化率提高,促使切向風(fēng)次極值形成,角動量收支分析也表明臺風(fēng)外圍潛熱有利于次眼墻的形成。ZHU, et al[13]在模擬試驗(yàn)中抹去了臺風(fēng)外雨帶處2~3 km高度上的非絕熱加熱,發(fā)現(xiàn)次眼墻處的次切向風(fēng)極值無法形成,表明了外雨帶低層對流和加熱對次眼墻形成的重要性,而低層對流、切向風(fēng)加速和邊界層上層入流輻合、地表蒸發(fā)之間的正反饋能夠促進(jìn)外雨帶對流發(fā)展和次眼墻形成,且外雨帶加熱引起的內(nèi)側(cè)出流會導(dǎo)致主眼墻減弱。CHEN, et al[14]將臺風(fēng)內(nèi)核區(qū)主眼墻外的低層非絕熱冷卻增加30%,發(fā)現(xiàn)內(nèi)核區(qū)低層存在低熵空氣和輻散下沉運(yùn)動,而外圍產(chǎn)生顯著輻合和活躍的對流雨帶,誘發(fā)了內(nèi)核風(fēng)場外擴(kuò),有利于次眼墻的形成。

    WANG, et al[15]提出了外雨帶驅(qū)動雙眼墻形成的關(guān)鍵動力機(jī)制:當(dāng)臺風(fēng)初始外圍尺度較大時,對流結(jié)構(gòu)和風(fēng)場結(jié)構(gòu)才能通過外雨帶和邊界層的相互作用耦合,形成完整次眼墻;如果初始外圍尺度較小,臺風(fēng)增強(qiáng)過程中主要以內(nèi)雨帶為主,無法出現(xiàn)明顯的外雨帶結(jié)構(gòu),相應(yīng)只能在主眼墻外形成次對流環(huán),無法形成切向風(fēng)次極大值結(jié)構(gòu),即為“偽雙眼墻”;當(dāng)進(jìn)一步縮小初始外圍尺度時,主眼墻外無法形成外雨帶,主要以內(nèi)雨帶為主,相應(yīng)無法形成次對流環(huán)和切向風(fēng)次極大值結(jié)構(gòu)時,表明不存在雙眼墻。顯然,臺風(fēng)初始外圍尺度能夠影響外雨帶的建立,進(jìn)而控制著次眼墻形成。臺風(fēng)增強(qiáng)過程中,外雨帶出現(xiàn)時間和位置可由此決定。

    非絕熱加熱對臺風(fēng)次眼墻的形成中至關(guān)重要,但目前為止,關(guān)于臺風(fēng)初始風(fēng)場結(jié)構(gòu)如何影響或者控制次眼墻形成位置和出現(xiàn)時間的討論尚不深入。初始外圍尺度如何通過影響雨帶及其非絕熱加熱分布進(jìn)而對次眼墻的形成產(chǎn)生影響?這一問題仍需深入探討。

    本文基于一系列的臺風(fēng)理想數(shù)值模擬,探究不同初始外圍尺度的臺風(fēng)中次眼墻形成位置和出現(xiàn)時間特征,分析主眼墻、外雨帶加熱的分布及其對次眼墻形成的影響并探討其中的動力機(jī)制。

    1 模式簡介和試驗(yàn)設(shè)計

    本文采用熱帶氣旋數(shù)值模式TCM4[16]進(jìn)行臺風(fēng)雙眼墻的理想模擬試驗(yàn)。該模式使用多重嵌套可移動網(wǎng)格,是完全可壓、非靜力平衡的理想三維熱帶氣旋數(shù)值模式,在熱帶氣旋研究中有廣泛應(yīng)用[17-20]。模式層頂高度為38 km,無擾動地面氣壓為1 010 hPa。本文模擬均設(shè)定在北緯18°N的f平面下進(jìn)行,采用四層嵌套網(wǎng)格,隔層網(wǎng)格為雙向反饋,水平網(wǎng)格從外至內(nèi)的分辨率分別為54、18、6和2 km;相應(yīng)水平格點(diǎn)數(shù)分別為281×241、181×181、217×217、271×271。積分步長360 s,垂直方向?yàn)榉蔷鶆?2層,其中對流層低層分辨率較高。模擬試驗(yàn)均處于靜止環(huán)境場,海表面溫度為29 ℃,采用相同的參數(shù)化方案,具體設(shè)置如下:次網(wǎng)格垂直混合采用E-ε湍流閉合方案[21]、地面通量計算采用Monin-Obukhov通量方案[22]、微物理方案采用顯式混合相云微物理方案[23],除質(zhì)量守恒方程以外,其他方程均采用非線性四階水平擴(kuò)散方案,位溫方程使用牛頓冷卻項來模擬輻射冷卻作用[24]。

    模擬采用的初始風(fēng)廓線分布為[17]:

    (1)

    其中:r為模擬渦旋的半徑;R0為風(fēng)速為0的半徑;Vm為最大切向風(fēng)速半徑rm處的風(fēng)速;b為初始風(fēng)廓線形狀參數(shù),表征初始渦旋最大風(fēng)速半徑外切向風(fēng)的徑向衰減速率,b值越大代表渦旋外圍風(fēng)速隨半徑衰減越快,相應(yīng)的渦旋外圍尺度越小。本文模擬試驗(yàn)中僅通過改變外圍風(fēng)廓線參數(shù)b來控制模擬臺風(fēng)的初始渦旋外圍尺度大小。

    圖1為模擬試驗(yàn)中采用的部分初始渦旋廓線分布,其中展示的初始渦旋廓線的b值分別為0.2、0.27、0.34、0.41、0.48、0.55、0.6、0.7和0.8。實(shí)際模擬試驗(yàn)中b值改變范圍為0.1~1.0之間,模擬臺風(fēng)的總數(shù)為30個。初始最大切向風(fēng)速為30 m·s-1,初始最大風(fēng)速半徑為135 km。為此,b值分別設(shè)置為0.1、0.12、0.16、0.2、0.21、0.22、0.23、0.24、0.25、0.26、0.27、0.28、0.29、0.3、0.31、0.32、0.33、0.34、0.35、0.36、0.37、0.38、0.41、0.44、0.45、0.48、0.55、0.6、0.7和0.8,相應(yīng)的模擬臺風(fēng)試驗(yàn)名稱分別為C01、C012、C016、…、C07和C08。由于b值不同,模擬臺風(fēng)初始17 m·s-1風(fēng)圈半徑(R17)不同[25-26],從1 200 km逐漸下降到300 km左右,表征臺風(fēng)初始外圍尺度的減小。針對不同的初始最大切向風(fēng)速和最大風(fēng)速半徑的模擬試驗(yàn)均得到了類似結(jié)果。

    圖1 模擬渦旋初始切向風(fēng)的徑向廓線

    2 試驗(yàn)結(jié)果概述

    表1給出部分模擬臺風(fēng)演變過程中重要特征參數(shù)。其中,次眼墻形成時間定義為最大切向風(fēng)速次極大值大于等于最大切向風(fēng)速主極大值的時間[27]。次眼墻形成位置定義為最大切向風(fēng)速次極大值的最初出現(xiàn)位置。Moat區(qū)寬度定義為最大切向風(fēng)速主、次極大值之間的距離。按照WANG, et al[15],將臺風(fēng)演變過程中出現(xiàn)次對流環(huán)但沒有出現(xiàn)最大切向風(fēng)的次極大值的現(xiàn)象稱為“偽雙眼墻”。由于偽雙眼墻未形成與最大切向風(fēng)速主極大值分離的次極大值,因此,偽次眼墻生成時間和位置定義為為最大切向風(fēng)速半徑在外圍開始二次收縮時的時間和位置。

    在表1所示的模擬臺風(fēng)中,C01—C038試驗(yàn)發(fā)生完整雙眼墻,C04—C055試驗(yàn)發(fā)生偽雙眼墻,C06—C08試驗(yàn)沒有出現(xiàn)雙眼墻。顯然,臺風(fēng)初始外圍尺度較大時,有利于次眼墻的出現(xiàn)。

    表1 模擬臺風(fēng)的特征參數(shù)

    2.1 次眼墻形成的徑向位置和時間

    本節(jié)主要分析模擬臺風(fēng)中初始外圍尺度對次眼墻形成位置、時間的影響。圖2中r2表示切向風(fēng)次極大值建立時的徑向位置,r1表示此時切向風(fēng)主極大值的位置,r2-r1表示moat區(qū)寬度。隨著b值從0.1增加到0.55,臺風(fēng)初始外圍尺度減小,次眼墻建立位置更加靠近臺風(fēng)中心,從150 km下降到86 km,此時,主眼墻也逐步收縮。同時,moat區(qū)的寬度也呈現(xiàn)變小趨勢。當(dāng)初始外圍尺度減少一定大小時,無法出現(xiàn)完整雙眼墻。次眼墻出現(xiàn)位置和b值的線性擬合效果較好,二者的相關(guān)系數(shù)為-0.92,具有極強(qiáng)的負(fù)相關(guān)性。因此,外圍風(fēng)廓線衰減率參數(shù)b所決定的臺風(fēng)初始外圍尺度大小,對雙眼墻的切向風(fēng)次極大值能否形成、形成位置具有主導(dǎo)作用。

    圖2 主、次眼墻位置和b值的關(guān)系:r1為切向風(fēng)主極大值的徑向位置;r2為切向風(fēng)次極大值的徑向位置。紅色三角標(biāo)記偽雙眼墻試驗(yàn)

    圖3給出了雙眼墻出現(xiàn)時間隨參數(shù)b的變化。可知,次眼墻出現(xiàn)時間隨臺風(fēng)初始外圍尺度減小而推遲。雙眼墻出現(xiàn)時間與參數(shù)b之間線性擬合較好,相關(guān)系數(shù)為0.82,有較高的正相關(guān)性。因此,模擬臺風(fēng)初始外圍尺度越小,次眼墻出現(xiàn)就越晚。

    圖3 次眼墻形成時間和b值的關(guān)系(紅色三角標(biāo)記偽雙眼墻試驗(yàn))

    不同初始風(fēng)廓線的臺風(fēng)模擬試驗(yàn)結(jié)果表明,次眼墻形成的時間、位置與渦旋初始外圍尺度密切相關(guān):隨著初始外圍尺度遞減,模擬臺風(fēng)從能夠出現(xiàn)完整雙眼墻逐步過渡到偽雙眼墻、沒有雙眼墻。另外,對于出現(xiàn)次眼墻的臺風(fēng),其雙眼墻出現(xiàn)時間也隨之推遲且位置更加靠近臺風(fēng)中心。顯然,初始外圍風(fēng)場大小對臺風(fēng)次眼墻的出現(xiàn)、出現(xiàn)的徑向位置和時間具有主導(dǎo)作用。

    2.2 次眼墻的切向風(fēng)增強(qiáng)

    上述模擬結(jié)果表明初始風(fēng)場外圍尺度對最大切向風(fēng)次極大值、次眼墻形成的位置和時間具有主導(dǎo)作用。下面具體分析初始風(fēng)場結(jié)構(gòu)對次眼墻形成的影響機(jī)制,為此選取了C02(完整雙眼墻)和C055(偽雙眼墻)試驗(yàn)作為分析對象。

    圖4給出C02和C055試驗(yàn)降水率水平分布及軸對稱平均的加熱、切向風(fēng)垂直分布。C02試驗(yàn)在71.3 h存在內(nèi)外兩個對流環(huán)(圖4a),77 h主眼墻對流逐漸消散。軸對稱平均切向風(fēng)垂直分布上存在兩個極值中心,外側(cè)切向風(fēng)速大值中心超過50 m·s-1(圖4c),對應(yīng)次眼墻切向風(fēng)逐步建立的過程。74 h后次眼墻向內(nèi)收縮,主眼墻減弱,隨后發(fā)生眼墻替換。在C055試驗(yàn)中,外圍雨帶形成位置更加靠近主眼墻,moat區(qū)較狹窄,78 h左右主眼墻外的雨帶逐漸并入眼墻主雨帶(圖4b),主雨帶逐漸豐厚。此時與眼墻相關(guān)的切向風(fēng)極大值中心向外擴(kuò)張非常緩慢,無法形成兩個分離的大值中心(圖4d),在此意義上C055僅為偽雙眼墻。

    圖4 臺風(fēng)降水率分布(單位:mm·h-1)(距離中心每50 km以灰色圓圈標(biāo)示)以及非絕熱加熱率(填色;單位:K·h-1)和切向風(fēng)(等值線;單位:m·s-1)分布:(a、c)C02(71.3 h);(b、d)C055(78 h)

    圖5為C02、C055試驗(yàn)中軸對稱平均的切向風(fēng)、絕對渦度徑向輸送和非絕熱加熱隨時間的演變。加熱大值主要位于RMW外側(cè),而絕對渦度徑向輸送大值位于RMW內(nèi)側(cè)。C02試驗(yàn)中次眼墻生成前,切向風(fēng)、絕對渦度徑向輸送、非絕熱加熱在主眼墻處隨時間減弱,而在次眼墻形成區(qū)域不斷增強(qiáng),與WANG, et al[15]結(jié)果類似。與C02相比,C055有較大差別,盡管在主眼墻外側(cè)出現(xiàn)次加熱大值(圖5f),但絕對渦度徑向輸送大值始終維持在主眼墻及RMW附近(圖5d),導(dǎo)致主眼墻處的切向風(fēng)強(qiáng)度近乎不變,無法出現(xiàn)切向風(fēng)次極大中心(圖5b),所以C055試驗(yàn)中沒有出現(xiàn)完整雙眼墻,只出現(xiàn)了偽雙眼墻。

    圖5 模擬臺風(fēng)C02、C055軸對稱平均結(jié)構(gòu)的時間演變:(a、b)模式底層切向風(fēng)(填色; 單位:m·s-1);(c、d)模式底層絕對渦度徑向輸送(填色; 單位:10-3 m·s-2);(e、f)整層垂直積分的非絕熱加熱率(填色; 單位:K·h-1)。黑色實(shí)線為RMW

    為進(jìn)一步討論雙眼墻的切向風(fēng)次極大值形成,對切向風(fēng)傾向進(jìn)行診斷。軸對稱切向風(fēng)傾向診斷方程為[28]:

    (2)

    圖6表明切向風(fēng)傾向的模式結(jié)果(圖6a、b)和診斷結(jié)果(圖6c、d)基本一致。平均平流和摩擦項之和的分布(圖6e、f)表明了其對次眼墻形成區(qū)域切向風(fēng)增長的貢獻(xiàn)較大。

    圖6 C02在第66—72 h、C055在第72—78 h平均的(a、b)模式切向風(fēng)傾向,(c、d)切向風(fēng)診斷方程計算得到的切向風(fēng)傾向以及(e、f)平均平流與摩擦項之和(單位:m·s-1·h-1)

    圖7、8分別為C02、C055試驗(yàn)中切向風(fēng)診斷方程右邊各項分布。在C02次眼墻形成前6 h(第66—72 h),其切向風(fēng)收支分析表明平均絕對渦度徑向輸送對低層切向風(fēng)速增長、次眼墻形成具有重要貢獻(xiàn)(圖7a、g)。平均絕對渦度徑向輸送在邊界層內(nèi)是正貢獻(xiàn),主眼墻附近最大,外雨帶處對應(yīng)次極大值,而在邊界層上方由于出流,絕對渦度輸送使切向風(fēng)減速(圖7a)。動量垂直輸送項的貢獻(xiàn)則相反,眼墻處邊界層內(nèi)為負(fù),邊界層上為正,外雨帶處也有類似分布特征,但強(qiáng)度要弱(圖7b、e、h)。渦動相關(guān)項之和在邊界層主眼墻內(nèi)側(cè)為正,眼墻處為負(fù),moat區(qū)為正,moat區(qū)以外為負(fù)貢獻(xiàn)(圖7f)。摩擦耗散在邊界層主要為負(fù)貢獻(xiàn),邊界層上為正貢獻(xiàn)較小(圖7i)。平均平流項和摩擦項之和在眼墻外側(cè)為正貢獻(xiàn),moat區(qū)為負(fù)貢獻(xiàn),次眼墻處為正貢獻(xiàn)(圖6e)。

    圖7 模擬臺風(fēng)C02在第66—72 h平均的:(a)平均絕對渦度徑向輸送項;(b)平均動量垂直輸送項;(c)平均絕對渦度徑向、平均動量垂直輸送項之和;(d)渦動絕對渦度徑向輸送項;(e)渦動動量垂直輸送項;(f)渦動絕對渦度徑向、垂直輸送項之和,(g)平均、渦動絕對渦度徑向輸送項之和;(h)平均、渦動動量的垂直輸送項之和;(i)摩擦耗散項(單位:m·s-1h-1)

    切向風(fēng)診斷結(jié)果表明,次眼墻形成前低層絕對渦度徑向輸送對切向風(fēng)傾向的正貢獻(xiàn)較大,邊界層摩擦抵消這部分的貢獻(xiàn),二者之和基本決定了外雨帶處切向風(fēng)增長,這些結(jié)果與QIU, et al[29]、WANG, et al[15]類似。

    圖8為C055試驗(yàn)切向風(fēng)場向外圍擴(kuò)張前6 h(第72—78 h)的切向風(fēng)傾向方程中各項貢獻(xiàn)項分布。這一時間段中,外圍雨帶處形成的入流次極大值中心與主極大值很近,且由于平均絕對渦度徑向輸送次極大值未能與主眼墻的極大值中心分離(圖8a),而恰恰這個區(qū)域與摩擦耗散大值區(qū)域相重疊(圖8i),平均絕對渦度徑向輸送的正貢獻(xiàn)很大程度上被摩擦耗散的負(fù)貢獻(xiàn)所抵消,不利于次眼墻形成區(qū)域的切向風(fēng)增長,無法建立起切向風(fēng)的次極大值,相應(yīng)C055只能為偽雙眼墻。

    圖8 同圖7,但為模擬臺風(fēng)C055在第72—78 h平均的切向風(fēng)診斷結(jié)果

    綜上,平均絕對渦度徑向輸送和邊界層摩擦的位置、大小配置主要控制了雙眼墻形成區(qū)域的切向風(fēng)增強(qiáng)傾向。絕對渦度徑向輸送中渦度由切向風(fēng)場決定,入流由加熱和摩擦強(qiáng)迫得到,而非絕熱加熱發(fā)展較早、變化較快,對次眼墻形成和發(fā)展具有指示性作用。

    3 主、次眼墻的非絕熱加熱對次眼墻演變的影響

    前面分析可知,臺風(fēng)外圍尺度對次眼墻形成位置和時間具有主導(dǎo)作用。臺風(fēng)初始外圍尺度不同,導(dǎo)致主眼墻外側(cè)對流發(fā)展的位置和時間不同,相應(yīng)雙眼墻相關(guān)的次對流環(huán)、切向風(fēng)次極大值建立的位置和時間也不相同。臺風(fēng)外圍對流引起的非絕熱加熱可以通過影響邊界層入流進(jìn)而影響絕對渦度徑向輸送,進(jìn)一步促進(jìn)次眼墻形成區(qū)域的切向風(fēng)增長。切向風(fēng)次極大值形成與否、形成位置決定于絕對渦度徑向輸送是否被摩擦耗散抵消,因此,主眼墻和雨帶非絕熱加熱的分布對次眼墻的形成至關(guān)重要。

    到目前為止,雙眼墻形成過程中非絕熱加熱對低層徑向入流、切向風(fēng)增長的貢獻(xiàn)和影響機(jī)制尚不夠清晰,也并未細(xì)致區(qū)分不同風(fēng)場結(jié)構(gòu)中主、次眼墻加熱對次眼墻形成發(fā)展的貢獻(xiàn)以及主、次眼墻加熱之間的相互作用。為此,本節(jié)將重點(diǎn)探討不同初始風(fēng)場結(jié)構(gòu)中主、次眼墻的非絕熱加熱對次眼墻形成的影響,從加熱率分布的角度分析外圍尺度對臺風(fēng)次眼墻發(fā)生發(fā)展的重要性。為此采用平衡動力學(xué)方法來討論不同雨帶加熱對邊界層入流、絕對渦度徑向輸送以及雙眼墻的切向風(fēng)次極大值增強(qiáng)的影響作用。

    利用Sawyer-Eliassen(SE)方程診斷在不同加熱條件下臺風(fēng)的次級環(huán)流強(qiáng)度和切向風(fēng)傾向。SE方程如下[20]:

    (3)

    圖9表明SE方程診斷結(jié)果相比模式結(jié)果有一定的低估,但分布特征基本相似。從垂直速度來看(圖9a),動量強(qiáng)迫(紅色虛線)在主眼墻處相比其他位置有較大貢獻(xiàn),但總體貢獻(xiàn)比熱量強(qiáng)迫的貢獻(xiàn)要?。辉诖窝蹓μ?,熱量強(qiáng)迫有貢獻(xiàn)峰值,但動量強(qiáng)迫貢獻(xiàn)較小。從絕對渦度徑向輸送來看(圖9b),主眼墻處動量強(qiáng)迫貢獻(xiàn)較大,但在次眼墻形成區(qū)域,動量強(qiáng)迫與熱量強(qiáng)迫的貢獻(xiàn)相當(dāng)。從入流分布來看(圖9c),主眼墻處動量與熱量強(qiáng)迫的貢獻(xiàn)大致相當(dāng),但在次眼墻處,動量強(qiáng)迫貢獻(xiàn)要大于熱量強(qiáng)迫的貢獻(xiàn)。

    圖9 C02試驗(yàn)73 h,SE方程診斷出的物理量的徑向分布:(a)1.5 km高度垂直速度(單位:m·s-1);(b)0~1 km高度平均的絕對渦度徑向輸送(單位:10-3 m·s-2);(c)徑向風(fēng)速(單位:m·s-1)

    圖10給出了C02、C035、C055試驗(yàn)的SE方程診斷得到的次級環(huán)流分布。C02、C035均選取次眼墻建立時刻,分別為第72 h和79 h,而C055形成偽雙眼墻,因此選取主眼墻停止向內(nèi)收縮、風(fēng)場外擴(kuò)后的時間,具體為第81 h。圖10a—c為熱量和動量共同強(qiáng)迫的結(jié)果,圖d-f僅為熱量強(qiáng)迫,圖10g—i僅為動量強(qiáng)迫。顯然,熱量強(qiáng)迫對垂直速度分布起主導(dǎo)作用,在主、次眼墻處有明顯大值區(qū)。C055只存在一個入流極大值中心,約為-12 m·s-1,位于內(nèi)雨帶加熱大值的外側(cè)。動量強(qiáng)迫對主眼墻及外圍入流貢獻(xiàn)較大,大值中心均分布在主眼墻處,在3個試驗(yàn)中大小相似(圖10c、f、i),主眼墻處約為-6 m·s-1,次眼墻處約為-5~-6 m·s-1。

    圖10 SE方程診斷的垂直速度(填色;單位:m·s-1)與徑向風(fēng)(等值線;單位:m·s-1)的次級環(huán)流分布:(a-c)C02試驗(yàn)(73 h);(d-f)C035試驗(yàn)(79 h);(g-i)C055試驗(yàn)(81 h)

    為進(jìn)一步分析主、次眼墻加熱和動量強(qiáng)迫對邊界層入流的貢獻(xiàn),分離動量強(qiáng)迫和主、次眼墻的加熱后重新計算SE方程,其中C02主眼墻范圍為40~80 km,次眼墻范圍為110~260 km;C035主、次眼墻范圍分別為40~78 km和80~260 km,C055分別為40~78 km和78~260 km。圖11為3個試驗(yàn)的次級環(huán)流分布。若僅保留主眼墻加熱,引起的邊界層入流徑向范圍很窄,其大值中心主要位于主眼墻附近(圖11a、g、m)。若僅有次眼墻加熱,會在次眼墻形成區(qū)域強(qiáng)迫出較明顯的入流,形成大值中心,在次眼墻內(nèi)側(cè)強(qiáng)迫出出流(圖11d、j、p)。若保留主眼墻處的熱量和動量強(qiáng)迫,主眼墻外的入流大小增加約一倍,但徑向范圍擴(kuò)張不大(圖11b、h、n)。若保留主眼墻加熱和整個徑向范圍的動量強(qiáng)迫,則主眼墻垂直速度增加,入流伸展到臺風(fēng)外圍,主眼墻外側(cè)入流增加,主眼墻外圍垂直速度也存在較小正值(圖11c、i、o)。若保留次眼墻加熱和動量強(qiáng)迫,雨帶垂直速度增加,在次眼墻外側(cè)引起較大入流,比僅有次眼墻加熱時的入流中心約大4 m·s-1,而次眼墻內(nèi)側(cè)的出流中心不再出現(xiàn)(圖11e、k、q)。若保留次眼墻加熱和全場動量強(qiáng)迫,次眼墻外圍入流的中心位置不會改變,內(nèi)側(cè)的出流中心轉(zhuǎn)變?yōu)槿肓髦行?圖11f、l、r)。

    圖11 SE方程診斷的次級環(huán)流分布。強(qiáng)迫項分別為C02:試驗(yàn)(a)僅主眼墻加熱;(b)僅主眼墻處加熱和動量強(qiáng)迫;(c)主眼墻處加熱和全場動量強(qiáng)迫;(d)僅次眼墻處加熱;(e)僅次眼墻處加熱和動量強(qiáng)迫;(f)次眼墻處加熱和全場動量強(qiáng)迫。(g—l)C035試驗(yàn);(m—r)C055試驗(yàn)

    綜上,主眼墻加熱激發(fā)出的入流僅限于主眼墻外側(cè)較小范圍內(nèi),對次眼墻外圍入流的影響較小。次眼墻加熱越大,能強(qiáng)迫出的外圍入流越大,且入流中心位置和模式結(jié)果對應(yīng)。C02、C035的主、次眼墻加熱中心徑向距離較遠(yuǎn),能夠強(qiáng)迫出兩個徑向分離的低層入流中心。C055試驗(yàn)的內(nèi)雨帶加熱和主眼墻加熱距離很近,低層入流的極大值中心只有一個,位于內(nèi)雨帶加熱大值的外側(cè)。動量強(qiáng)迫引起的入流大值中心均僅有一個,位于主眼墻加熱率外側(cè)。由此可得,次眼墻處的加熱和動量強(qiáng)迫決定了次眼墻形成區(qū)域入流中心的位置和大小,對形成絕對渦度徑向輸送次極大值的貢獻(xiàn)較大。

    從圖12中可以看到,與C02 相比(圖12b),C055試驗(yàn)中偽次眼墻處熱量強(qiáng)迫出的入流峰值更靠近主眼墻(圖12e),主眼墻加熱對外圍絕對渦度徑向輸送的貢獻(xiàn)很小(圖12f),因而對外圍切向風(fēng)增長的影響較小,無法形成切向風(fēng)的次極大值。動量強(qiáng)迫對垂直速度的貢獻(xiàn)與加熱強(qiáng)迫相比要小(圖12a、d),但對入流和絕對渦度徑向輸送有較大影響(圖12b、c、e、f),然而由于動量強(qiáng)迫出的入流無法在次眼墻區(qū)域形成一個明顯的次極大值中心(圖10c、f、i),此時考慮到動量強(qiáng)迫下的入流次極大值中心仍緊靠主眼墻,因此,次眼墻形成區(qū)域的絕對渦度徑向輸送大小同時受次眼墻加熱和動量強(qiáng)迫控制,但次眼墻出現(xiàn)的位置主要由次眼墻加熱主導(dǎo)。

    圖12 C02試驗(yàn)(73 h),1.5 km高度(a)垂直速度(單位:m·s-1)、(b)徑向風(fēng)速(m·s-1)與(c)0~1 km高度平均的絕對渦度徑向輸送(單位:10-3 m·s-2)的分布(黑線:主眼墻加熱,綠線:主眼墻加熱和動量強(qiáng)迫,紅線:次眼墻加熱,黃線:次眼墻加熱和動量強(qiáng)迫;(d—f): 同a—c,但為C055試驗(yàn)(81 h))

    綜上,次眼墻處加熱對次眼墻處絕對渦度徑向輸送以及切向風(fēng)次極大值的出現(xiàn)位置具有主導(dǎo)影響,而切向風(fēng)次極大值大小與非絕熱加熱和動量強(qiáng)迫均有關(guān)。

    4 結(jié)論

    臺風(fēng)雙眼墻形成及其演變過程對臺風(fēng)強(qiáng)度及結(jié)構(gòu)變化有著重要影響作用,是臺風(fēng)預(yù)報的難點(diǎn)之一。臺風(fēng)次眼墻具有兩個基本結(jié)構(gòu)特征:一是主眼墻外側(cè)形成一個新閉合的對流環(huán),二是在低層形成與切向風(fēng)主極大值相分離的切向風(fēng)次極大值。對于一個臺風(fēng)出現(xiàn)完整的次眼墻,必須同時出現(xiàn)次對流環(huán)和切向風(fēng)次極大值,且兩者完全耦合。已有研究表明,臺風(fēng)外雨帶是臺風(fēng)次眼墻形成的最重要強(qiáng)迫,外雨帶通過其對流引起的非絕熱加熱影響臺風(fēng)邊界層徑向入流,進(jìn)而對絕對渦度徑向輸送產(chǎn)生影響,導(dǎo)致在主眼墻外側(cè)形成一個切向風(fēng)次極大值。本文對控制次眼墻形成位置和時間的因素進(jìn)行討論,重點(diǎn)討論了初始風(fēng)場外圍尺度及次眼墻加熱對次眼墻形成與發(fā)展影響作用。主要結(jié)論如下:

    (1)在一系列不同初始風(fēng)場結(jié)構(gòu)的臺風(fēng)理想模擬中,研究發(fā)現(xiàn)次眼墻形成的時間、位置與初始渦旋的外圍尺度顯著相關(guān):隨著初始風(fēng)場外圍尺度遞減,臺風(fēng)結(jié)構(gòu)從能夠形成完整雙眼墻向偽雙眼墻、沒有雙眼墻逐步過渡,且次眼墻出現(xiàn)時間逐步推遲、形成位置更加靠近臺風(fēng)中心。所以,渦旋初始外圍尺度大小對臺風(fēng)次眼墻形成位置和出現(xiàn)時間具有主導(dǎo)作用。

    (2)動力學(xué)診斷分析進(jìn)一步發(fā)現(xiàn),初始風(fēng)場結(jié)構(gòu)控制著臺風(fēng)外雨帶分布,雨帶的非絕熱加熱主導(dǎo)著主眼墻外側(cè)低層徑向入流和絕對渦度徑向輸送的分布和強(qiáng)度。而絕對渦度徑向輸送和摩擦耗散的相對大小及位置決定低層切向風(fēng)次極值(次眼墻)出現(xiàn)的位置和時間。此外,主眼墻加熱所引起的邊界層徑向入流其水平范圍較窄,對次眼墻形成區(qū)域的絕對渦度徑向輸送的貢獻(xiàn)較??;而次眼墻加熱主導(dǎo)了次眼墻外側(cè)的徑向入流和絕對渦度徑向輸送,二者大小受次眼墻加熱和動量強(qiáng)迫共同控制。非絕熱加熱對次眼墻形成區(qū)域的切向風(fēng)次極大值的出現(xiàn)位置具有主導(dǎo)影響,而非絕熱加熱和動量強(qiáng)迫共同貢獻(xiàn)切向風(fēng)次極大值大小。

    (3)綜上,初始外圍尺度大小決定臺風(fēng)發(fā)展階段的外雨帶徑向分布,進(jìn)而決定與之相關(guān)的非絕熱加熱徑向分布,進(jìn)一步控制了邊界層入流和絕對渦度徑向輸送,主導(dǎo)了切向風(fēng)增強(qiáng)以及切向風(fēng)大值區(qū)的出現(xiàn)位置和時間,從而進(jìn)一步對次眼墻的建立和發(fā)展產(chǎn)生重要影響。

    由于臺風(fēng)雙眼墻的建立與發(fā)展具有非線性、復(fù)雜性,除了初始外圍尺度會對臺風(fēng)強(qiáng)度和結(jié)構(gòu)演變、特別是對外雨帶發(fā)展產(chǎn)生影響,其他初始條件、環(huán)境條件也可對外雨帶等產(chǎn)生重要影響,在這些復(fù)雜條件下,臺風(fēng)雙眼墻形成的動力學(xué)機(jī)制需要進(jìn)一步的深入研究。

    猜你喜歡
    雨帶渦度極大值
    王金成:秋風(fēng)
    南海相對渦度的時空變化特征分析*
    斜壓渦度的變化與臺風(fēng)暴雨的關(guān)系研究
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    基于小波模極大值理論的勵磁涌流新判據(jù)研究
    基于經(jīng)驗(yàn)?zāi)B(tài)分解的自適應(yīng)模極大值去噪方法
    行人檢測中非極大值抑制算法的改進(jìn)
    各種各樣的雨
    一次東北南部特大暴雪的動力機(jī)制分析
    最大熵算法在氣象雨量預(yù)測中的應(yīng)用分析
    满18在线观看网站| 男女边摸边吃奶| 狂野欧美激情性xxxx| 91成年电影在线观看| 日韩 欧美 亚洲 中文字幕| 黄色成人免费大全| av国产精品久久久久影院| 亚洲国产av影院在线观看| 国产精品一区二区在线观看99| 九色亚洲精品在线播放| 黄色毛片三级朝国网站| 51午夜福利影视在线观看| 国产不卡一卡二| 我要看黄色一级片免费的| 亚洲av国产av综合av卡| 亚洲专区国产一区二区| 嫁个100分男人电影在线观看| 亚洲色图综合在线观看| 久9热在线精品视频| 亚洲精品美女久久久久99蜜臀| 国产精品1区2区在线观看. | 亚洲国产欧美网| √禁漫天堂资源中文www| www.999成人在线观看| 日韩制服丝袜自拍偷拍| 精品久久久精品久久久| 久久精品91无色码中文字幕| 国产亚洲欧美在线一区二区| 啦啦啦免费观看视频1| netflix在线观看网站| 久久亚洲真实| 日韩欧美一区视频在线观看| 中文字幕av电影在线播放| 青草久久国产| 黄网站色视频无遮挡免费观看| 黄色丝袜av网址大全| 国产精品偷伦视频观看了| 欧美亚洲日本最大视频资源| 欧美日韩国产mv在线观看视频| 51午夜福利影视在线观看| 久久狼人影院| 精品第一国产精品| 女性生殖器流出的白浆| 又紧又爽又黄一区二区| 国产麻豆69| 久久精品成人免费网站| 性高湖久久久久久久久免费观看| 国产主播在线观看一区二区| 纯流量卡能插随身wifi吗| 超碰成人久久| 青草久久国产| 大型黄色视频在线免费观看| 黄色丝袜av网址大全| av福利片在线| 一夜夜www| 日韩欧美国产一区二区入口| 2018国产大陆天天弄谢| 国产一区二区激情短视频| 精品亚洲乱码少妇综合久久| 水蜜桃什么品种好| 日日摸夜夜添夜夜添小说| 午夜老司机福利片| 男女边摸边吃奶| 成人免费观看视频高清| 亚洲国产毛片av蜜桃av| 日日夜夜操网爽| 久久精品91无色码中文字幕| 国内毛片毛片毛片毛片毛片| 欧美在线黄色| 老司机在亚洲福利影院| 首页视频小说图片口味搜索| 超色免费av| 欧美av亚洲av综合av国产av| 手机成人av网站| 满18在线观看网站| av视频免费观看在线观看| 欧美变态另类bdsm刘玥| 亚洲av成人不卡在线观看播放网| 女性被躁到高潮视频| netflix在线观看网站| 国产老妇伦熟女老妇高清| 免费一级毛片在线播放高清视频 | 国产黄频视频在线观看| 亚洲av美国av| 亚洲第一av免费看| 91九色精品人成在线观看| 精品国内亚洲2022精品成人 | 久久99热这里只频精品6学生| 成人黄色视频免费在线看| 日韩大码丰满熟妇| 国产精品免费一区二区三区在线 | 日韩视频在线欧美| 亚洲欧美激情在线| 美女国产高潮福利片在线看| 日日爽夜夜爽网站| 亚洲一区二区三区欧美精品| videosex国产| 中文字幕人妻丝袜一区二区| 亚洲国产成人一精品久久久| 精品久久久精品久久久| 丰满迷人的少妇在线观看| 国产亚洲欧美精品永久| 亚洲av第一区精品v没综合| 亚洲av电影在线进入| 十八禁网站网址无遮挡| 男女无遮挡免费网站观看| 女人高潮潮喷娇喘18禁视频| 亚洲精品国产区一区二| 一二三四在线观看免费中文在| 咕卡用的链子| 看免费av毛片| 一级毛片电影观看| 777米奇影视久久| 一级,二级,三级黄色视频| 久久久国产精品麻豆| 丁香欧美五月| 精品国产亚洲在线| 国产精品久久久久成人av| 91精品国产国语对白视频| 国产一区二区在线观看av| 中文字幕人妻熟女乱码| 久久久久国内视频| 国产伦人伦偷精品视频| 国产三级黄色录像| 日韩欧美三级三区| 亚洲美女黄片视频| 国产精品久久久久久精品古装| 国产精品熟女久久久久浪| xxxhd国产人妻xxx| 伦理电影免费视频| 大码成人一级视频| 免费人妻精品一区二区三区视频| 岛国毛片在线播放| 另类亚洲欧美激情| 国产在线观看jvid| 久久久久久久国产电影| 在线观看免费视频日本深夜| 久久ye,这里只有精品| 欧美亚洲 丝袜 人妻 在线| 国产精品av久久久久免费| 亚洲九九香蕉| 日韩视频一区二区在线观看| 亚洲成人免费av在线播放| 久久精品亚洲熟妇少妇任你| 午夜91福利影院| 韩国精品一区二区三区| 一级毛片精品| 久久ye,这里只有精品| 久久久精品94久久精品| 99久久精品国产亚洲精品| 亚洲国产精品一区二区三区在线| 久久毛片免费看一区二区三区| 一个人免费看片子| 国产精品麻豆人妻色哟哟久久| 欧美日韩亚洲高清精品| 国产精品久久久久久精品电影小说| 久久精品亚洲精品国产色婷小说| 国产aⅴ精品一区二区三区波| 免费久久久久久久精品成人欧美视频| 欧美一级毛片孕妇| 日本wwww免费看| 黄色丝袜av网址大全| 中文字幕高清在线视频| 成年人午夜在线观看视频| 成人特级黄色片久久久久久久 | 国产日韩欧美亚洲二区| 搡老岳熟女国产| 亚洲久久久国产精品| 日本a在线网址| 亚洲熟女精品中文字幕| 欧美在线一区亚洲| 国产亚洲精品久久久久5区| 黑丝袜美女国产一区| 国产一区有黄有色的免费视频| 中文字幕最新亚洲高清| 精品久久久久久电影网| 在线看a的网站| 久久久久久久久久久久大奶| 啪啪无遮挡十八禁网站| 久久久精品免费免费高清| 午夜福利一区二区在线看| 欧美变态另类bdsm刘玥| h视频一区二区三区| 丁香欧美五月| 日本五十路高清| 日韩欧美国产一区二区入口| 亚洲全国av大片| 嫩草影视91久久| 嫩草影视91久久| 欧美日韩亚洲高清精品| 在线 av 中文字幕| 9热在线视频观看99| 色婷婷久久久亚洲欧美| 久久影院123| 看免费av毛片| 丰满迷人的少妇在线观看| 亚洲精品乱久久久久久| 日本一区二区免费在线视频| 麻豆av在线久日| 最近最新免费中文字幕在线| 天天操日日干夜夜撸| 日韩视频一区二区在线观看| 亚洲人成77777在线视频| 成年人免费黄色播放视频| 大香蕉久久网| 国产免费现黄频在线看| 国产亚洲精品第一综合不卡| 亚洲欧美激情在线| 日韩免费av在线播放| 精品国内亚洲2022精品成人 | 国产成人免费无遮挡视频| 国产高清国产精品国产三级| 黄色丝袜av网址大全| 色视频在线一区二区三区| 9色porny在线观看| 国产99久久九九免费精品| 一级毛片精品| 不卡av一区二区三区| 亚洲午夜精品一区,二区,三区| 国产精品98久久久久久宅男小说| 丰满饥渴人妻一区二区三| 国产精品免费视频内射| 国产精品久久久久久精品古装| 99久久人妻综合| 精品亚洲乱码少妇综合久久| 精品人妻1区二区| h视频一区二区三区| kizo精华| av欧美777| 国产一区有黄有色的免费视频| 少妇被粗大的猛进出69影院| 亚洲第一av免费看| 国产一区二区三区综合在线观看| 精品久久久久久久毛片微露脸| 黄频高清免费视频| av网站免费在线观看视频| 欧美精品高潮呻吟av久久| 激情视频va一区二区三区| 侵犯人妻中文字幕一二三四区| 欧美日韩黄片免| 欧美 日韩 精品 国产| 免费黄频网站在线观看国产| 国产亚洲欧美在线一区二区| 久9热在线精品视频| 国产一区二区 视频在线| av天堂久久9| 俄罗斯特黄特色一大片| 久久中文字幕一级| 午夜福利视频精品| 麻豆国产av国片精品| 捣出白浆h1v1| 免费在线观看影片大全网站| 色综合婷婷激情| 一级毛片电影观看| 嫁个100分男人电影在线观看| 久久精品国产亚洲av高清一级| 亚洲熟女毛片儿| 亚洲欧洲日产国产| 这个男人来自地球电影免费观看| 黄色 视频免费看| 中文字幕另类日韩欧美亚洲嫩草| av视频免费观看在线观看| 一级a爱视频在线免费观看| 亚洲av成人一区二区三| 精品少妇内射三级| 最近最新中文字幕大全电影3 | 老熟女久久久| 成人精品一区二区免费| 久久久精品94久久精品| 久久免费观看电影| 视频在线观看一区二区三区| 亚洲一码二码三码区别大吗| avwww免费| 亚洲中文字幕日韩| 在线观看免费视频日本深夜| 国产亚洲av高清不卡| 日韩制服丝袜自拍偷拍| 欧美黄色片欧美黄色片| 波多野结衣av一区二区av| kizo精华| 国产成人精品久久二区二区91| 桃花免费在线播放| 久久精品亚洲熟妇少妇任你| 又黄又粗又硬又大视频| 狠狠精品人妻久久久久久综合| 国产av一区二区精品久久| 国产亚洲精品久久久久5区| 国产极品粉嫩免费观看在线| 成人18禁在线播放| 黑人猛操日本美女一级片| 黄片播放在线免费| 免费在线观看日本一区| 9191精品国产免费久久| 妹子高潮喷水视频| 亚洲国产成人一精品久久久| 亚洲熟妇熟女久久| 亚洲人成77777在线视频| 免费在线观看黄色视频的| 亚洲人成伊人成综合网2020| 男女高潮啪啪啪动态图| 男男h啪啪无遮挡| 亚洲精品av麻豆狂野| 又紧又爽又黄一区二区| 丁香六月天网| 视频区欧美日本亚洲| 欧美精品一区二区免费开放| 国产精品 欧美亚洲| 国产在视频线精品| 成人国产一区最新在线观看| 亚洲久久久国产精品| 交换朋友夫妻互换小说| 桃花免费在线播放| 国产精品欧美亚洲77777| 熟女少妇亚洲综合色aaa.| 亚洲avbb在线观看| 亚洲色图综合在线观看| 久久精品亚洲精品国产色婷小说| 午夜视频精品福利| 另类亚洲欧美激情| 亚洲国产精品一区二区三区在线| 国产欧美日韩综合在线一区二区| av一本久久久久| 欧美激情久久久久久爽电影 | netflix在线观看网站| 不卡一级毛片| 国产真人三级小视频在线观看| 日韩视频一区二区在线观看| 久久天躁狠狠躁夜夜2o2o| 日韩一卡2卡3卡4卡2021年| 中文字幕av电影在线播放| 久久精品亚洲熟妇少妇任你| 999久久久国产精品视频| 国产午夜精品久久久久久| 欧美精品人与动牲交sv欧美| aaaaa片日本免费| 午夜91福利影院| 韩国精品一区二区三区| 国产三级黄色录像| 一进一出抽搐动态| 亚洲国产欧美网| 伊人久久大香线蕉亚洲五| 国产不卡av网站在线观看| 久久人妻福利社区极品人妻图片| 国产成人一区二区三区免费视频网站| 国产av国产精品国产| 少妇粗大呻吟视频| 免费在线观看完整版高清| 极品教师在线免费播放| 无人区码免费观看不卡 | 男女午夜视频在线观看| 亚洲av电影在线进入| aaaaa片日本免费| 美国免费a级毛片| 精品亚洲成国产av| 久久国产精品大桥未久av| 亚洲伊人久久精品综合| 日日夜夜操网爽| 乱人伦中国视频| 91av网站免费观看| 日本av免费视频播放| 国产高清国产精品国产三级| 亚洲av成人一区二区三| 国产一区有黄有色的免费视频| 亚洲久久久国产精品| 在线观看免费午夜福利视频| 欧美国产精品一级二级三级| 亚洲精品一卡2卡三卡4卡5卡| 丝瓜视频免费看黄片| 久久国产精品大桥未久av| 成人手机av| av一本久久久久| 12—13女人毛片做爰片一| 啦啦啦 在线观看视频| 无人区码免费观看不卡 | 久久天堂一区二区三区四区| 1024视频免费在线观看| 中文字幕最新亚洲高清| 免费在线观看影片大全网站| 亚洲中文av在线| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美网| 国产在线观看jvid| 久久精品国产亚洲av高清一级| 菩萨蛮人人尽说江南好唐韦庄| 9191精品国产免费久久| 国产男靠女视频免费网站| 日韩熟女老妇一区二区性免费视频| 亚洲人成电影免费在线| 可以免费在线观看a视频的电影网站| 成人精品一区二区免费| 嫩草影视91久久| 精品国产乱码久久久久久男人| e午夜精品久久久久久久| 久久精品亚洲精品国产色婷小说| 国产又色又爽无遮挡免费看| a在线观看视频网站| 窝窝影院91人妻| 99国产极品粉嫩在线观看| 国产av精品麻豆| av片东京热男人的天堂| 国产精品国产高清国产av | 我要看黄色一级片免费的| 好男人电影高清在线观看| 欧美黑人欧美精品刺激| 亚洲一卡2卡3卡4卡5卡精品中文| 一级毛片精品| 精品人妻1区二区| 少妇裸体淫交视频免费看高清 | 极品少妇高潮喷水抽搐| 久久99一区二区三区| 欧美大码av| 999精品在线视频| 老汉色av国产亚洲站长工具| 欧美日韩福利视频一区二区| 高清欧美精品videossex| 国产主播在线观看一区二区| 国产精品久久久久久人妻精品电影 | 国产免费现黄频在线看| 十分钟在线观看高清视频www| 老熟妇乱子伦视频在线观看| 国产高清激情床上av| 亚洲熟女精品中文字幕| 精品高清国产在线一区| 波多野结衣av一区二区av| 电影成人av| 亚洲欧美日韩高清在线视频 | 如日韩欧美国产精品一区二区三区| 丁香六月欧美| h视频一区二区三区| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 成人精品一区二区免费| 丝瓜视频免费看黄片| 色94色欧美一区二区| 操出白浆在线播放| 50天的宝宝边吃奶边哭怎么回事| 婷婷成人精品国产| 天堂动漫精品| 成人黄色视频免费在线看| 窝窝影院91人妻| 韩国精品一区二区三区| 亚洲av第一区精品v没综合| 亚洲精品国产色婷婷电影| 亚洲精品成人av观看孕妇| 国产一区二区三区综合在线观看| 性少妇av在线| 激情在线观看视频在线高清 | 国产精品麻豆人妻色哟哟久久| 久久久久久久久久久久大奶| 亚洲 欧美一区二区三区| 无遮挡黄片免费观看| 一区二区三区精品91| 又黄又粗又硬又大视频| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 美女视频免费永久观看网站| 国产av精品麻豆| 成人亚洲精品一区在线观看| 69av精品久久久久久 | 国产在线一区二区三区精| 高清黄色对白视频在线免费看| 国产日韩一区二区三区精品不卡| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜制服| 国产精品一区二区免费欧美| 色在线成人网| 日韩大码丰满熟妇| 99久久99久久久精品蜜桃| 午夜福利乱码中文字幕| 一进一出抽搐动态| 国产精品免费视频内射| 亚洲av电影在线进入| 91大片在线观看| 欧美av亚洲av综合av国产av| 一区二区三区精品91| 国产精品久久久久成人av| 亚洲久久久国产精品| 亚洲综合色网址| 亚洲 国产 在线| 欧美亚洲 丝袜 人妻 在线| 香蕉久久夜色| 日本黄色日本黄色录像| e午夜精品久久久久久久| 成年人午夜在线观看视频| 久久中文看片网| 国产国语露脸激情在线看| 国产精品久久久久成人av| 国产高清激情床上av| 后天国语完整版免费观看| 欧美日韩国产mv在线观看视频| 成人国产一区最新在线观看| 国产黄色免费在线视频| 色视频在线一区二区三区| 涩涩av久久男人的天堂| 亚洲国产欧美一区二区综合| 人人澡人人妻人| 免费黄频网站在线观看国产| 国产成人av教育| 黄色视频在线播放观看不卡| 露出奶头的视频| 亚洲欧洲日产国产| 亚洲少妇的诱惑av| 国产欧美日韩一区二区精品| 手机成人av网站| 宅男免费午夜| av超薄肉色丝袜交足视频| 日韩免费高清中文字幕av| 日本精品一区二区三区蜜桃| 99国产精品99久久久久| 999精品在线视频| 天天操日日干夜夜撸| 两个人看的免费小视频| 天天躁夜夜躁狠狠躁躁| 久久精品国产亚洲av香蕉五月 | 新久久久久国产一级毛片| 一区二区三区精品91| 手机成人av网站| 国内毛片毛片毛片毛片毛片| 国产在线精品亚洲第一网站| 欧美黑人精品巨大| 精品国产一区二区久久| 国产无遮挡羞羞视频在线观看| 亚洲av第一区精品v没综合| 亚洲欧美日韩高清在线视频 | 国产色视频综合| 国产精品国产高清国产av | 国产一卡二卡三卡精品| 国产精品.久久久| svipshipincom国产片| 在线观看一区二区三区激情| 国产精品偷伦视频观看了| 高清欧美精品videossex| 国产欧美日韩精品亚洲av| 在线播放国产精品三级| 免费看a级黄色片| 99香蕉大伊视频| 亚洲精品在线美女| 欧美精品一区二区大全| 少妇猛男粗大的猛烈进出视频| 欧美午夜高清在线| 国产精品秋霞免费鲁丝片| 侵犯人妻中文字幕一二三四区| 国产亚洲精品久久久久5区| 中文字幕av电影在线播放| 欧美激情 高清一区二区三区| 一本色道久久久久久精品综合| 大型黄色视频在线免费观看| 国产精品偷伦视频观看了| 亚洲专区字幕在线| 国产精品九九99| 亚洲色图综合在线观看| 天堂动漫精品| 亚洲情色 制服丝袜| 久久精品aⅴ一区二区三区四区| 人人妻人人澡人人看| 国产无遮挡羞羞视频在线观看| 亚洲国产欧美一区二区综合| 国产精品熟女久久久久浪| 国产精品av久久久久免费| 色综合婷婷激情| 国产精品.久久久| 亚洲欧美日韩高清在线视频 | 日本vs欧美在线观看视频| 少妇被粗大的猛进出69影院| 免费看十八禁软件| 国产亚洲av高清不卡| 免费观看人在逋| 真人做人爱边吃奶动态| 无遮挡黄片免费观看| 69精品国产乱码久久久| a级毛片黄视频| 久久精品国产亚洲av高清一级| av天堂久久9| 久久精品aⅴ一区二区三区四区| 国产男女内射视频| 波多野结衣av一区二区av| 国产成人av教育| 三上悠亚av全集在线观看| 黄色片一级片一级黄色片| 国产精品影院久久| 亚洲人成伊人成综合网2020| 亚洲三区欧美一区| 国产黄色免费在线视频| 女同久久另类99精品国产91| 夜夜爽天天搞| 亚洲一区中文字幕在线| 18禁美女被吸乳视频| 一本综合久久免费| 午夜福利在线免费观看网站| 亚洲午夜理论影院| 俄罗斯特黄特色一大片| 一区二区三区乱码不卡18| 日韩视频一区二区在线观看| av天堂久久9| 久久精品熟女亚洲av麻豆精品| 777米奇影视久久| 搡老岳熟女国产| √禁漫天堂资源中文www| 久久精品国产亚洲av香蕉五月 | 欧美成狂野欧美在线观看| 另类精品久久| 国产熟女午夜一区二区三区| 黄色 视频免费看| 18禁黄网站禁片午夜丰满| 国产日韩欧美亚洲二区| 国产精品亚洲av一区麻豆| 极品少妇高潮喷水抽搐| 狂野欧美激情性xxxx| 精品视频人人做人人爽| 在线播放国产精品三级| a级片在线免费高清观看视频| av一本久久久久| 欧美黑人欧美精品刺激| 国产精品电影一区二区三区 | 日本精品一区二区三区蜜桃|