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

    華東沿海地形對登陸熱帶氣旋運動影響的理想數(shù)值研究

    2018-09-10 16:52:38陳小宇吳立廣吳玉琴
    大氣科學(xué)學(xué)報 2018年4期

    陳小宇 吳立廣 吳玉琴

    摘要利用中尺度數(shù)值模式設(shè)計一組高分辨率理想試驗,采用位渦趨勢方法定量診斷分析熱帶氣旋在登陸我國華東沿海地形時,其運動發(fā)生的精細(xì)化變化以及不同因子的貢獻(xiàn)。結(jié)果表明,平地的存在使得登陸熱帶氣旋移速相對更快,當(dāng)華東沿海地形存在時,熱帶氣旋移速顯著增大,這種增速現(xiàn)象主要是由于平地和地形所引起的非對稱氣流以及相應(yīng)的引導(dǎo)氣流變化所致,這很可能是導(dǎo)致預(yù)報路徑誤差的一個重要原因。平地試驗中,陸地在熱帶氣旋低層激發(fā)出中小尺度的非對稱氣流,與之不同的是,實際地形的加入激發(fā)出更大尺度并且更強(qiáng)的非對稱偏南氣流。位渦趨勢方法的診斷結(jié)果表明,非引導(dǎo)效應(yīng)總體而言對熱帶氣旋運動貢獻(xiàn)較小,這是因為這些因子相互抵消,但在不同的垂直層次上,不同的非引導(dǎo)因子貢獻(xiàn)存在明顯的差異。

    關(guān)鍵詞登陸熱帶氣旋運動;地形;非對稱氣流;位渦趨勢方法

    熱帶氣旋(TC)的路徑和運動一直以來都是熱帶氣旋預(yù)報和研究中最為重要的方面之一,路徑預(yù)報的準(zhǔn)確度直接關(guān)系著強(qiáng)度、降水以及災(zāi)害等方面的預(yù)報結(jié)果。早期的觀點大多認(rèn)為熱帶氣旋的運動主要受大尺度環(huán)境氣流的引導(dǎo)效應(yīng)所影響(Adem and Lezama,1960;George and Gray,1976;Neumanm,1979;Chan and Gray,1982)。但許多學(xué)者在研究中發(fā)現(xiàn),觀測中熱帶氣旋的運動方向和環(huán)境氣流的方向之間始終存在一定的偏差,并且這種偏差并不是由于引導(dǎo)氣流計算方式和范圍的差別所引起的(George and Gray,1976;DeMaria,1985;Chan and Williams,1987;Fiorino and Elsberry,1989)。早在20世紀(jì)80年代初,Holland(1983)就提出了β漂移(βdrift)理論,對這一現(xiàn)象的機(jī)制進(jìn)行了解釋,發(fā)現(xiàn)行星渦度的經(jīng)向梯度在熱帶氣旋內(nèi)部激發(fā)出一對渦旋對,渦旋對之間的向極向西的通風(fēng)流是導(dǎo)致熱帶氣旋運動與環(huán)境引導(dǎo)氣流存在偏差的原因。這一理論得到后續(xù)大量的觀測和數(shù)值研究的驗證(Holland,1983;Chan and Williams,1987;Fiorino and Elsberry,1989;Willoughby,1990,1992;Peng and Williams,1990;Li and Wang,1994;Wang et al.,1997)。

    除了環(huán)境導(dǎo)引和β漂移之外,Wu and Wang(2000)發(fā)現(xiàn)熱帶氣旋運動還受到非引導(dǎo)氣流效應(yīng)的影響。他們提出了一種新的診斷方法,把熱帶氣旋看作具有較強(qiáng)正位勢渦度(PV)異常的系統(tǒng),因此將熱帶氣旋運動與位勢渦度的趨勢變化聯(lián)系起來,這種方法被稱為位渦趨勢(PVT)方法。該方法不僅考慮傳統(tǒng)的引導(dǎo)氣流和β漂移的作用,還包含了垂直運動、非絕熱加熱以及摩擦等非引導(dǎo)效應(yīng)對于熱帶氣旋運動的貢獻(xiàn)。Wu and Wang(2001a,2001b)利用這一方法進(jìn)一步分析了絕熱加熱和非絕熱加熱情形下,引導(dǎo)和非引導(dǎo)效應(yīng)對于熱帶氣旋運動的貢獻(xiàn)。Chan et al.(2002)采用該方法診斷觀測中熱帶氣旋路徑變化時,發(fā)現(xiàn)不同的因子對于各種類型的熱帶氣旋路徑貢獻(xiàn)有所差別。此外,一些研究在分析熱帶氣旋登陸的理想試驗和個例模擬中利用該診斷方法區(qū)別不同物理過程的作用(Wong and Chan,2006;吳玉琴等,2015)。這些結(jié)果表明,位渦趨勢方法相比于傳統(tǒng)的引導(dǎo)氣流方法,可以更好地應(yīng)用于熱帶氣旋運動的診斷分析中。

    登陸熱帶氣旋給我國沿海地區(qū)的人民生命財產(chǎn)安全造成了巨大的損失(Zhang et al.,2009;張嬌艷等,2011;曹劍等,2012;李肖雅等,2014),因此對于登陸熱帶氣旋路徑的研究直接關(guān)系著災(zāi)害損失和防臺減災(zāi)工作的展開。Tuleya and Kurihara(1978)最早設(shè)計了理想試驗以研究登陸熱帶氣旋的特征變化。在這之后,許多研究針對島嶼和山脈地形(例如臺灣島)進(jìn)行了大量的觀測和數(shù)值研究(Chang,1982;Bender et al.,1987;Yeh and Elsberry,1993a,1993b;Lin et al.,2005;Jian and Wu,2008;Yang et al.,2008;Huang et al.,2011;Wang et al.,2012;Xie and Zhang,2012)。與之相對應(yīng)的是,較少有研究涉及大陸地形對熱帶氣旋路徑的影響。針對這一問題,Wong and Chan(2006)設(shè)計了幾組理想試驗,研究在沒有背景氣流和地形的影響下,熱帶氣旋登陸f平面的平坦陸地過程中運動的變化。他們發(fā)現(xiàn)陸地的存在激發(fā)熱帶氣旋環(huán)流產(chǎn)生不對稱的結(jié)構(gòu),這種不對稱性逐漸發(fā)展,最終使得熱帶氣旋產(chǎn)生向陸地的運動趨勢,運動速度約為1 m/s。這些結(jié)果表明,陸地使得熱帶氣旋的不對稱結(jié)構(gòu)發(fā)展從而可能導(dǎo)致登陸過程中移速增加。后續(xù)的數(shù)值研究進(jìn)一步驗證,在β平面和三角洲地形影響下這種運動加速現(xiàn)象的存在(AuYeung and Chan,2010;Szeto and Chan,2010),但是這種研究中同樣沒有考慮環(huán)境氣流的影響。針對這一問題,Li et al.(2013)在試驗中加入了均勻的環(huán)境氣流,發(fā)現(xiàn)了類似的現(xiàn)象,但是他們在研究中著重關(guān)注熱帶氣旋登陸過程中降水的變化,缺乏對于運動變化的定量性分析。在熱帶氣旋登陸大陸地形的過程中,較少有研究考慮環(huán)境氣流的影響,尤其當(dāng)環(huán)境氣流和地形高度同時存在時,會對登陸熱帶氣旋運動產(chǎn)生何種影響,仍然缺乏更深入的認(rèn)識。

    近三十年間,登陸我國的熱帶氣旋陸上持續(xù)時間和移動速度都出現(xiàn)了明顯的變化(Chen et al.,2011),一些研究結(jié)果也表明,登陸我國大陸與海南島的熱帶氣旋造成了重大的直接經(jīng)濟(jì)損失和人員傷亡(Zhang et al.,2009;張嬌艷等,2011),但是對于登陸我國大陸沿海地區(qū)的熱帶氣旋的運動研究仍然偏少。因此,本文設(shè)計高分辨率的理想試驗?zāi)M熱帶氣旋登陸我國華東沿海地區(qū),通過對比試驗結(jié)果定量地診斷背景氣流和華東地區(qū)實際地形對于登陸熱帶氣旋運動的影響,以及不同的物理過程對于運動變化的貢獻(xiàn)。

    1資料與試驗設(shè)計

    使用WRFARW V 221模式,試驗設(shè)計采用了五層網(wǎng)格雙向嵌套設(shè)置,最外層網(wǎng)格為水平分辨率27 km,格點數(shù)為230×210,網(wǎng)格中心點為1325°E、300°N。四層內(nèi)網(wǎng)格的水平分辨率分別為9、3、1、1/3 km,格點數(shù)分別為432×399、333×333、501×501以及720×720。模式的垂直層次為40層,模式頂氣壓為50 hPa。最外兩層網(wǎng)格為固定網(wǎng)格,其范圍設(shè)置足夠包含了西北太平洋海域的主要系統(tǒng)的發(fā)展和演變;最內(nèi)三層網(wǎng)格為移動網(wǎng)格,伴隨著模擬的熱帶氣旋運動而移動,主要覆蓋了逐個時次上的熱帶氣旋影響范圍。

    選取2005年典型熱帶氣旋“麥莎”活動并登陸的2005年8月5—9日這一時段,利用Lanczos時間濾波方法(Duchon,1979)得出該時段的20 d以上低頻背景場作為試驗的背景場(圖1)。8月5號00時(世界時,下同),“麥莎”位于臺灣東北側(cè)洋面之上,此時該熱帶氣旋的近地面風(fēng)速約為45 m·s-1,在后續(xù)的3 d中,“麥莎”沿著環(huán)境場氣流向西北方向移動并于8月5日19時40分在我國東南沿海地區(qū)登陸,登陸之后強(qiáng)度開始下降并逐漸消散。初始時刻,該低頻背景場呈現(xiàn)出典型的西北太平洋夏季的環(huán)流形式,主要表現(xiàn)為較為明顯的季風(fēng)槽和副高的分布特征。需要指出的是由于環(huán)境場已經(jīng)過濾波處理(已去除20 d以下頻段的系統(tǒng)),因此后續(xù)登陸試驗中的TC路徑與強(qiáng)度均與臺風(fēng)“麥莎”存在明顯區(qū)別。

    最外層網(wǎng)格的微物理參數(shù)化方案和積云對流方案分別采用WSM3簡單冰方案和KainFritropical cycloneh方案(Dudhia,1989;Kain and Fritch,1993),第二至第五層網(wǎng)格采用WSM6微物理方案并關(guān)閉積云對流方案(Hong and Lim,2006)。邊界層過程選用的是Yonsei University 邊界層方案(Noh et al.,2003),輻射過程選擇了Dudhia 短波方案(Dudhia,1989)、Rapid Radiative Transfer Model(RRTM)長波方案(Mlawer et al.,1997),陸面過程采用的是Noah陸面過程模式。模式初始化選用美國國家環(huán)境預(yù)報中心(National Centers for Environmental Prediction,NCEP)發(fā)布的分辨率為1°×1°的FNL全球分析資料(Final Operational Global Analysis),該資料時間間隔為6 h。試驗中對最外層網(wǎng)格進(jìn)行張弛逼近(Nudging),以便積分過程的背景場更為穩(wěn)定。最外層網(wǎng)格的積分時間步長為90 s,總積分時間為96 h。

    在該低頻背景場中,放入人工構(gòu)造的理想軸對稱渦旋。為了進(jìn)一步改進(jìn)人造渦旋的結(jié)構(gòu)和風(fēng)壓關(guān)系,本研究中,預(yù)先將構(gòu)造的理想軸對稱渦旋放入一個f平面無環(huán)境場的全海洋試驗中,讓渦旋在該環(huán)境中進(jìn)行初旋(Spinup),使其結(jié)構(gòu)和風(fēng)壓關(guān)系得到模式動力與熱力過程的調(diào)整和改善。進(jìn)行初旋之前,構(gòu)造的理想渦旋的海平面最低氣壓為994 hPa,近地面最大風(fēng)速約30 m·s-1。經(jīng)過12 h的初旋后,該渦旋結(jié)構(gòu)依然保持對稱,海平面最低氣壓下降至965 hPa,近地面最大風(fēng)速達(dá)到40 m·s-1,已達(dá)到了臺風(fēng)強(qiáng)度(圖2)。之后從初旋試驗結(jié)果中計算得出該渦旋的風(fēng)場、溫度場、濕度場以及氣壓場等物理變量,加入圖1所示低頻背景場之中與熱帶氣旋“麥莎”(8月5日00時)相近的位置(1226°E,251°N),使其沿著季風(fēng)槽與副高之間的東南氣流向西北方向運動,之后登陸我國華東沿海地區(qū)。

    在低頻背景場和理想渦旋的基礎(chǔ)上,設(shè)計一組理想試驗,以此對比分析在高分辨率的情況下,熱帶氣旋在登陸過程中運動的精細(xì)化變化。三個試驗分別為全海洋試驗(OCN),平地試驗(FLT)和地形試驗(TOP)。全海洋試驗中,將整個區(qū)域的下墊面和表層屬性全部設(shè)置成海洋屬性;平地試驗中,去除所有的地形高度,將其全部設(shè)置為0;而在地形試驗中,包括真實的地形高度和分布。同時,在所有試驗中,海平面溫度統(tǒng)一設(shè)置為29 ℃,在平地和地形試驗中,地表溫度也同樣設(shè)置為29 ℃,這么設(shè)置是為了排除海平面和地表溫度的差異所可能造成的熱帶氣旋登陸過程中的差異。同時,在平地和地形試驗中,陸表屬性也設(shè)置為均一的類型,本研究中選擇的陸表類型是我國華中華北地區(qū)常見的灌溉型農(nóng)田與牧場(irrigated cropland and pasture),該類型的地表摩擦層厚度為007 m。

    2模擬結(jié)果

    圖3給出了模擬的海洋試驗、平地試驗以及地形試驗中熱帶氣旋的積分前48 h移動路徑,其中熱帶氣旋中心定義為最大方位角平均的切向風(fēng)所對應(yīng)的中心(Wu et al.,2006)。為了剔除地形對低層氣旋環(huán)流的影響,三組實驗中將700 hPa氣旋中心的移動軌跡作為氣旋整體的移動路徑。從圖3可以看出,三組試驗中熱帶氣旋的運動方向一致但是移速存在區(qū)別,其中海洋試驗、平地實驗、地形試驗中熱帶氣旋分別在積分42 h、39 h、24 h到達(dá)海岸線,這表明在平地試驗、地形試驗中熱帶氣旋的移動速度較海洋試驗中更快。需要指出的是,在海洋試驗中去除了陸地和地形,因此在該試驗的分析中將熱帶氣旋移動至原陸地位置的時刻定義為該試驗的登陸時刻。

    從強(qiáng)度模擬結(jié)果(圖4)看,平地實驗、地形試驗中熱帶氣旋登陸后強(qiáng)度迅速減弱,48 h后,兩個試驗中熱帶氣旋近地面風(fēng)速降為熱帶風(fēng)暴級別。在海洋試驗中,渦旋在積分前12 h強(qiáng)度增加,并維持在40 m/s。與海洋試驗類似,平地試驗和地形試驗中氣旋在海洋中移動強(qiáng)度增加,在登陸前3 h氣旋強(qiáng)度開始減弱。熱帶氣旋登陸時,平地試驗、地形試驗中熱帶氣旋強(qiáng)度為310、313 m/s,登陸9 h后兩個試驗中強(qiáng)度減弱為194、176 m/s。積分48 h后,三個試驗中熱帶氣旋強(qiáng)度分別為:383、191、148 m/s。海平面氣壓場的變化表明,在登陸前24 h,三個試驗中熱帶氣旋的海平面最低氣壓接近,但這之后,海洋試驗熱帶氣旋繼續(xù)增強(qiáng),氣壓持續(xù)下降。與之相反的是,平地和地形試驗中的熱帶氣旋在登陸前氣壓即開始上升,登陸后12 h,海洋試驗、平地試驗和地形試驗中的熱帶氣旋海平面最低氣壓分別為950、976、989 hPa,同樣表明了地形試驗中熱帶氣旋強(qiáng)度下降的速率更大。

    吳玉琴等(2015)利用WRF模式模擬熱帶氣旋“莫拉克”(2008),通過考慮熱帶氣旋在移動過程中自身的PV一波發(fā)展改進(jìn)了原有的PVT診斷方法,并通過對比PVT估算速度、大尺度環(huán)境引導(dǎo)氣流與模擬的熱帶氣旋移速,發(fā)現(xiàn)改進(jìn)后的PVT估算速度較引導(dǎo)氣流能更好地表征熱帶氣旋的運動,特別是在熱帶氣旋運動出現(xiàn)振蕩時刻。因此,本文利用改進(jìn)后的PVT方法診斷理想試驗結(jié)果,重點分析熱帶氣旋登陸過程中引導(dǎo)效應(yīng)與非對稱結(jié)構(gòu)變化之間的關(guān)系。

    為了評估大尺度環(huán)境引導(dǎo)氣流與PVT估算速度在表征熱帶氣旋移動上的差異,在此分別計算了海洋、平地、地形數(shù)值試驗中熱帶氣旋移速、PVT估算速度、引導(dǎo)氣流的對比(圖5)。參考前人計算方法(Wu and Wang,2000;吳玉琴等,2015),計算PVT各項估算速度和引導(dǎo)氣流時采用第二層網(wǎng)格(D02)的固定坐標(biāo)系數(shù)據(jù),熱帶氣旋移速采用第五層高分辨率網(wǎng)格中的氣旋中心計算得出。此外,三個數(shù)值試驗中各垂直層上PVT估算速度均為以熱帶氣旋中心為中心,直徑為540 km的圓形區(qū)域計算,垂直方向上采用850~300 hPa多層平均,各層引導(dǎo)氣流采用相同的計算范圍和垂直層次所得,熱帶氣旋移速為前后兩小時位置差計算所得。

    32熱帶氣旋移速的差異

    三個數(shù)值試驗中,熱帶氣旋的移速出現(xiàn)明顯的差異。海洋、平地和地形試驗中氣旋在登陸前后的37 h平均移速分別為267、313、462 m/s,同時段平均引導(dǎo)氣流大小分別為251、271、423 m/s。這表明在平地試驗中即出現(xiàn)了熱帶氣旋登陸增速的現(xiàn)象,與Wong and Chan (2006)結(jié)果類似,但僅約05 m/s的增速幅度明顯小于他們所發(fā)現(xiàn)的1 m/s的增幅,這可能是由于本研究中所采用的地表摩擦層厚度僅為007 m,小于Wong and Chan (2006)的試驗設(shè)置。更值得注意的是,在地形試驗中,熱帶氣旋移速相對于海洋試驗中的結(jié)果增加了約2 m/s,表明實際地形的存在使得熱帶氣旋在登陸過程中移速顯著增大,并且這種現(xiàn)象很可能是由于引導(dǎo)氣流的增大所導(dǎo)致。楊瓊瓊和吳立廣(2015)在分析臺風(fēng)路徑預(yù)報誤差時發(fā)現(xiàn)大多數(shù)的預(yù)報誤差來源于臺風(fēng)移速的預(yù)報偏差而非移動方向的偏差,因此在環(huán)境氣流和地形共同影響下,熱帶氣旋登陸時運動所出現(xiàn)的變化很可能是導(dǎo)致路徑預(yù)報誤差的一個重要原因。

    分析試驗中熱帶氣旋移速的經(jīng)向和緯向分量后發(fā)現(xiàn),PVT估算速度較引導(dǎo)氣流更接近熱帶氣旋的移速(圖5),特別是在熱帶氣旋移動過程中出現(xiàn)速度突變情況下,例如海洋試驗和平地試驗中熱帶氣旋登陸前后6 h經(jīng)向和緯向速度均出現(xiàn)短時振蕩特征(圖3a—d),PVT估算速度很好地表征了熱帶氣旋移動中的速度變化特征,而傳統(tǒng)的引導(dǎo)氣流計算難以表征這一現(xiàn)象。熱帶氣旋登陸前24 h到登陸后12 h內(nèi),平地試驗中平均經(jīng)向、緯向速度分別為-123、286 m/s,平均PVT估算速度為-116、262 m/s,相應(yīng)的平均引導(dǎo)氣流為-057、263 m/s。PVT估算速度與熱帶氣旋移速之間均方根誤差分別為013、029 m/s,引導(dǎo)氣流與熱帶氣旋移速之間均方根誤差為066、04 m/s,從計算的平均移速、均方根誤差數(shù)值上都可以看出PVT估算速度相較于引導(dǎo)氣流能更好地表征熱帶氣旋的移動,并且在海洋試驗和地形試驗中存在類似的結(jié)果。這是由于PVT診斷方法不僅包含了大尺度環(huán)境引導(dǎo)氣流對熱帶氣旋運動的作用,還包含了熱帶氣旋結(jié)構(gòu)變化、非絕熱加熱、PV垂直平流等對熱帶氣旋運動的作用。需特別指出的是,在地形試驗中的熱帶氣旋經(jīng)向速度分析中,PVT速度與熱帶氣旋移速的均方根誤差略大于引導(dǎo)氣流的偏差,這主要是由于該試驗中熱帶氣旋強(qiáng)度的迅速下降所導(dǎo)致的,但PVT方法仍然可以更好地再現(xiàn)熱帶氣旋在登陸前后的緯向速度以及經(jīng)向速度的振蕩。結(jié)合上述結(jié)果,可利用PVT診斷方法診斷三個數(shù)值試驗中熱帶氣旋運動的差異。

    Wu and Wang(2000)指出局地PV一波變化由PV的水平平流、PV垂直平流、非絕熱加熱、摩擦四項組成,為了評估各物理過程對熱帶氣旋運動的作用,在此分別計算三組數(shù)值試驗中PV的水平平流(簡稱HA)、PV垂直平流(簡稱VA)、非絕熱加熱(簡稱DH)、摩擦(簡稱FR)對應(yīng)的估算速度與PVT估算速度的對比(圖6)。

    從圖6可以看出,三組數(shù)值試驗中PV水平平流對應(yīng)的估算速度最接近PVT估算速度,在海洋數(shù)值試驗中,PVT、HA、VA、DH、FR各項所對應(yīng)的緯向估算速度為-087、-021、-07、005、-001 m/s,經(jīng)向速度為242、261、038、-043、-014 m/s,相應(yīng)的引導(dǎo)氣流為-052、243 m/s。從數(shù)值上可以看出,HA對應(yīng)的估算速度接近引導(dǎo)氣流的大小,從長時間平均來看,HA對熱帶氣旋運動的作用要遠(yuǎn)大于其他三項的作用總和,此外,VA與DH對熱帶氣旋運動的作用是相互抵消的,這也導(dǎo)致了熱帶氣旋的移動基本由HA決定。與海洋試驗相同,平地試驗與地形試驗中HA對PVT的貢獻(xiàn)最大,但HA與引導(dǎo)氣流之間的差異相比海洋試驗增加,這可能是兩個試驗登陸過程中由熱帶氣旋非對稱結(jié)構(gòu)的發(fā)展引起的非引導(dǎo)效應(yīng)所導(dǎo)致的??傮w而言,PVT方法診斷的結(jié)果表明,HA項對于三個試驗熱帶氣旋的運動起著最重要的貢獻(xiàn),同時可以表征三個熱帶氣旋移速的明顯差異(圖6)。Wu and Wang (2001a)指出,HA項的貢獻(xiàn)主要表示了熱帶氣旋中的引導(dǎo)效應(yīng),但與傳統(tǒng)意義的引導(dǎo)氣流有所區(qū)別,不僅包括了大尺度環(huán)境引導(dǎo)氣流和β漂移的作用,還包含了熱帶氣旋內(nèi)垂直切變、垂直運動、絕熱和非絕熱加熱等過程所造成的非對稱氣流對熱帶氣旋運動的直接貢獻(xiàn)。

    33登陸熱帶氣旋的非對稱結(jié)構(gòu)

    為了進(jìn)一步分析產(chǎn)生登陸熱帶氣旋運動變化的原因,分別計算三組試驗中熱帶氣旋登陸時刻降水率與700 hPa的風(fēng)場進(jìn)行對比(圖7)。從圖7可以看出,海洋試驗中熱帶氣旋登陸時刻降水率與700 hPa風(fēng)場基本呈對稱結(jié)構(gòu),降水主要集中在距離熱帶氣旋中心120 km區(qū)域內(nèi),且熱帶氣旋眼墻北側(cè)區(qū)域內(nèi)降水最大(圖7a,b)。平地實驗中,登陸時刻熱帶氣旋的風(fēng)場、降水速率結(jié)構(gòu)上非對稱性增強(qiáng),大值區(qū)均在氣旋東北側(cè)(圖7c,d),地形試驗中風(fēng)場、降水分布與平地試驗存在相似結(jié)構(gòu),但強(qiáng)度上地形試驗的非對稱性要略弱于平地試驗。三個試驗中的降水大值均分布在順切變左側(cè)區(qū)域,與之前的研究結(jié)果一致,表明強(qiáng)對流更傾向于發(fā)生在順切變左側(cè)區(qū)域(Frank and Ritchie,2001;Chen and Yau,2003)。

    為了進(jìn)一步了解地形對熱帶氣旋非對稱結(jié)構(gòu)的影響,分別計算平地試驗和海洋試驗中不同時間段內(nèi)中低層熱帶氣旋的一波非對稱風(fēng)場差值以及降水率的差值(圖8a—d)。為了更深入地研究登陸前熱帶氣旋結(jié)構(gòu)的變化,在此將熱帶氣旋登陸前的24 h分為兩個時間段:第一個為登陸前24~13 h(b1時間段),第二個為登陸前12 h(b2時間段),下同。從圖中可以看出,平地試驗中b1時間段內(nèi),700 hPa一波非對稱風(fēng)場熱帶氣旋中心處有弱西北氣流,同時在熱帶氣旋西北側(cè)(陸地的區(qū)域)出現(xiàn)了氣旋式的風(fēng)場異常,平地試驗和海洋試驗中的降水速率分布差異是這種低層地形導(dǎo)致輻合氣流的結(jié)果。熱帶氣旋內(nèi)核區(qū)域的西南向的通風(fēng)流與兩個試驗中水平速度的差異相當(dāng)。b2時間段內(nèi),低層非對稱氣流的增強(qiáng)與地形試驗中熱帶氣旋在登陸之前的加速度是一致的,700 hPa上出現(xiàn)了較為明顯的南風(fēng)氣流,氣流左側(cè)的氣旋性環(huán)流異常與陸地的分布相對應(yīng)(圖6d)。但是這種非對稱結(jié)構(gòu)在中層難以被識別出,500 hPa上的南風(fēng)非對稱氣流特征明顯偏弱,且在熱帶氣旋內(nèi)核區(qū)域出現(xiàn)了北風(fēng)異常,與熱帶氣旋的移速增速趨勢相反。從中低層對比可以看出,熱帶氣旋在平地試驗的登陸過程中的增速現(xiàn)象主要來源于700 hPa上非對稱南風(fēng)氣流的貢獻(xiàn),說明平地試驗中地形所激發(fā)的非對稱結(jié)構(gòu)主要表現(xiàn)在低層。

    在地形試驗中,b1時間段內(nèi),相較于海洋試驗和平地試驗,地形導(dǎo)致的非對稱氣流強(qiáng)度更強(qiáng);非對稱的南風(fēng)氣流在低層到高層上都比平地試驗中的結(jié)果更顯著(圖9a—c)。另一個顯著特征是,非對稱氣流的半徑超過熱帶氣旋的半徑,達(dá)到約1 000 km半徑的范圍,雨帶的分布與平地試驗類似。在b2時間段內(nèi),700 hPa上非對稱氣流最顯著,這可能是因為熱帶氣旋低層環(huán)流受地形的影響最為顯著。結(jié)合上述分析:在地形試驗中,非對稱南風(fēng)氣流與平地試驗中的結(jié)果相比明顯增強(qiáng),非對稱南風(fēng)氣流不僅出現(xiàn)在中低層,在高層也有明顯的分布特征,同時氣流強(qiáng)度和覆蓋范圍也明顯增大,因此導(dǎo)致地形試驗中熱帶氣旋的移速最快。此外值得一提的是,非對稱氣流的半徑明顯大于熱帶氣旋的尺度,這一結(jié)果表明,地形不僅導(dǎo)致了熱帶氣旋出現(xiàn)了非對稱結(jié)構(gòu),還顯著地影響了環(huán)境風(fēng)場。這一現(xiàn)象體現(xiàn)了地形與平地對于登陸熱帶氣旋運動的不同影響程度,平地只能對中低層的熱帶氣旋環(huán)流產(chǎn)生影響,而地形可以對熱帶氣旋和環(huán)境氣流同時造成非常重要的影響,使得氣旋在登陸過程中產(chǎn)生明顯的增速現(xiàn)象。

    34非引導(dǎo)效應(yīng)的作用

    吳玉琴等(2015)曾指出非引導(dǎo)效應(yīng)對于臺風(fēng)“莫拉克”的移速振蕩有著非常重要的作用。那么在登陸過程中,其余非引導(dǎo)效應(yīng)起著何種作用?為什么熱帶氣旋運動主要受到引導(dǎo)氣流的影響?針對這些問題,為了進(jìn)一步定量地分析非引導(dǎo)效應(yīng)的作用,分別計算三個試驗中熱帶氣旋在b1時間段內(nèi)、b2時間段內(nèi)各非引導(dǎo)效應(yīng)的平均貢獻(xiàn),結(jié)果如圖10所示。

    Wu and Wang(2001a)進(jìn)一步將PV水平平流項分為兩部分:

    ∧1(-V·P)≈-V1·Ps-Vs·P1。

    其中:腳標(biāo)s和1表示對稱分量和一波分量;等式右端第一項為非對稱氣流引起的對稱PV分量的輸送,簡稱ASAF,包含大尺度引導(dǎo)氣流和通風(fēng)流對熱帶氣旋運動的影響,與一波風(fēng)場的區(qū)域平均引導(dǎo)氣流對應(yīng)較好;第二項為對稱氣流引起的非對稱PV分量的平流,簡稱AASF,包含了β渦對、垂直切變及非對稱加熱產(chǎn)生的非對稱氣流等物理過程對熱帶氣旋運動的作用。本研究中將引導(dǎo)氣流以外的各物理過程統(tǒng)稱為非引導(dǎo)效應(yīng),包括了AASF項,PV垂直輸送項(VA),非絕熱加熱項(DH)以及摩擦項(FR)。

    圖10海洋試驗(a,d)、平地試驗(b,e)、地形試驗(c,f)中登陸前24~13 h(a,b,c)和12~1 h(d,e,f)平均的多層平均位渦趨勢各非引導(dǎo)項估算移速(VA項估算速度,藍(lán)色;DH項估算速度,紅色;FR項目估算速度,灰色;AASF項估算速度,紫色;上述各項的矢量和,棕色)

    Fig.10Averaged nonsteering terms including vertical advection(VA,blue),diabatic heating(DH,red),friction(FR,gray),advection of the asymmetric PV structure by the symmetric flow(AASF,purple) and vector sum of nonsteering terms(brown) during (a,b,c)24—13 h and (d,e,f)12—1 h before landfall in (a,d)OCN,(b,e)FLT and (c,f)TOP experiments

    由圖10可以看出,三組試驗中,b1、b2時間段內(nèi),去除引導(dǎo)氣流后PVT各分項對熱帶氣旋運動的作用不一致,PVT各分量的矢量和小于引導(dǎo)氣流作用。海洋試驗中b1、b2時間段內(nèi)非引導(dǎo)效應(yīng)總貢獻(xiàn)的矢量和分別為052、022 m/s,平地試驗中兩個時間段矢量和分別為069、022 m/s(圖10中棕色矢量),而在兩組試驗中熱帶氣旋移動的平均速度大約3 m/s,因此海洋試驗、平地試驗中登陸時段內(nèi)非引導(dǎo)效應(yīng)對熱帶氣旋運動的貢獻(xiàn)遠(yuǎn)小于引導(dǎo)氣流的貢獻(xiàn)。在地形試驗中,b1、b2時間段內(nèi)非引導(dǎo)效應(yīng)的總貢獻(xiàn)分別為116、138 m/s,明顯大于海洋和平地試驗的結(jié)果,但在該試驗中熱帶氣旋移速明顯增大,因此非引導(dǎo)效應(yīng)對于移速的整體貢獻(xiàn)依然偏小。值得注意的是,由圖10可以看出:VA、DH項各自的估算速度數(shù)值上均大于或者接近非引導(dǎo)項的總貢獻(xiàn),但矢量方向上存在一定程度的反位相分布特征,這說明在登陸過程中,非引導(dǎo)效應(yīng)各項之間存在一定程度相互抵消,因此對熱帶氣旋運動的作用整體貢獻(xiàn)偏小,熱帶氣旋在登陸過程中總體仍受到引導(dǎo)效應(yīng)的作用。此外,非引導(dǎo)效應(yīng)導(dǎo)致熱帶氣旋有向西運動的趨勢,同時在地形試驗中,由于地形的存在使得熱帶氣旋登陸過程中FR項加大。

    圖10的結(jié)果表明垂直方向上多層平均的非引導(dǎo)效應(yīng)對熱帶氣旋運動的總體作用較小,但在各垂直層上是否有同樣的作用呢?Wu and Wang(2001b)在數(shù)值研究中發(fā)現(xiàn),加熱場在熱帶氣旋高低層呈現(xiàn)相反的作用。為了進(jìn)一步分析各物理因子在不同垂直層次上的不同貢獻(xiàn),分別診斷了三個試驗中,熱帶氣旋登陸前12 h內(nèi)高中低層PVT各分項對熱帶氣旋運動的貢獻(xiàn)(圖11)。從圖11可以看出:熱帶氣旋登陸前12 h內(nèi),三個試驗中PVT方程各項的高低層貢獻(xiàn)存在較大差異,例如HA項和ASAF項高低層(700 hPa與300 hPa)貢獻(xiàn)相反,DH項與ASAF項在同一層次貢獻(xiàn)相反,這些發(fā)現(xiàn)與Wu and Wang(2001b)的結(jié)論一致。需要指出的是,海洋和平地試驗中ASAF項在中層貢獻(xiàn)較小,而在高層對熱帶氣旋運動的作用均為負(fù)貢獻(xiàn),但在地形試驗的ASAF項中高層均為較為顯著的正貢獻(xiàn),這一結(jié)果與該試驗中引導(dǎo)氣流的顯著增大對應(yīng)。此外,兩個登陸試驗的低層DH項與海洋試驗結(jié)果存在明顯差異,這說明熱帶氣旋在登陸過程中受陸地和地形的影響,低層非絕熱加熱以及相應(yīng)的非對稱對流、降水等結(jié)構(gòu)發(fā)生了相應(yīng)的變化,例如圖7中的降水率差別所示。以上結(jié)果表明,雖然非引導(dǎo)效應(yīng)對于熱帶氣旋運動的整體貢獻(xiàn)相對于引導(dǎo)效應(yīng)偏小,但在不同的垂直層次上各項的作用存在明顯差異,尤其在登陸試驗中,DH項發(fā)生了明顯的變化。

    圖11海洋試驗(a,b,c)、平地試驗(d,e,f)、地形試驗(h,i,j)中登陸前12~1 h平均的300 hPa(a,d,h)、500 hPa(b,e,i)和700 hPa(c,f,j)位渦趨勢各項估算移速(VA項估算速度,藍(lán)色;DH項估算速度,紅色;FR項目估算速度,灰色;AASF項估算速度,紫色;ASAF估算速度,黑色)

    Fig.11Averaged PVT terms including vertical advection(VA,blue),diabatic heating(DH,red),friction(FR,gray),advection of the asymmetric PV by the symmetric flow(AASF,purple) and advection of the symmetric PV by the asymmetric flow(ASAF,black) at (a,d,h)300 hPa,(b,e,i)500 hPa and (c,f,j)700 hPa during 12—1h before landfall in (a,b,c)OCN,(d,e,f)FLT and (h,i,j)TOP experiments

    4結(jié)論與討論

    對于熱帶氣旋登陸過程中的運動變化一直以來都是備受關(guān)注的一個問題,但是長期以來,更多的研究關(guān)注于特殊的島嶼地形(例如臺灣島)對于路徑的影響,僅有少量研究涉及大陸地形的作用(Wong and Chan,2006;AuYeung and Chan,2010;Szeto and Chan,2010;Li et al.,2013)。然而這些研究中,通常又沒有考慮地形高度和環(huán)境氣流的影響。因此,本研究針對這一問題,設(shè)計了一組高分辨率的數(shù)值理想試驗,通過試驗結(jié)果對比以及位渦趨勢方法的診斷結(jié)果,分析研究登陸我國華東沿海地區(qū)的熱帶氣旋在登陸過程中的運動變化,以及不同的物理因子對于這種變化的不同貢獻(xiàn)。

    分析結(jié)果表明:當(dāng)?shù)湫偷牡皖l背景場存在時,在熱帶氣旋登陸前24 h至登陸后12 h這一段時間內(nèi),熱帶氣旋在平地試驗中的移速相比在海洋試驗中的移速增加了約05 m/s,在具有實際地形高度的試驗中,熱帶氣旋移速存在約2 m/s的增幅。進(jìn)一步的分析發(fā)現(xiàn),這些移速增加的現(xiàn)象主要是由地形所引起的非對稱氣流以及相應(yīng)的引導(dǎo)效應(yīng)變化所引起的。在熱帶氣旋登陸前,平地下墊面可以激發(fā)出中低層的非對稱氣流,這種熱帶氣旋尺度的非對稱氣流導(dǎo)致了引導(dǎo)氣流的增大,從而產(chǎn)生移速的加快。當(dāng)?shù)匦胃叨却嬖跁r,地形可以激發(fā)出更大尺度的非對稱南風(fēng)氣流,同時這些非對稱氣流更強(qiáng),從熱帶氣旋的低層發(fā)展到高層,因此造成了地形試驗中顯著增大的熱帶氣旋移速。需要指出的是,這種地形所造成的非對稱氣流同時包括了地形對于熱帶氣旋環(huán)流和對于環(huán)境氣流的分別貢獻(xiàn),如何定量化地區(qū)別兩者的不同貢獻(xiàn),有待今后的研究中進(jìn)一步的分析。

    除了引導(dǎo)效應(yīng)之外,非引導(dǎo)效應(yīng)對于登陸過程中的熱帶氣旋運動的整體貢獻(xiàn)相對較小,但是其作用不可忽視。位渦趨勢方法的診斷結(jié)果表明,這些非引導(dǎo)因子整體存在相互抵消的現(xiàn)象,因此對熱帶氣旋運動總貢獻(xiàn)偏小,這與前人研究結(jié)果一致。此外,在不同的垂直層次上,不同的非引導(dǎo)效應(yīng)的貢獻(xiàn)存在著明顯的差異,這些差異與熱帶氣旋登陸過程中強(qiáng)度以及結(jié)構(gòu)等方面的變化有著密切聯(lián)系,表明熱帶氣旋強(qiáng)度和結(jié)構(gòu)的變化同樣可能在一定程度會對其運動產(chǎn)生影響。

    綜合以上結(jié)論,可以看出我國大陸華東沿海地區(qū)的地形可以使得熱帶氣旋在登陸該區(qū)域時移速發(fā)生顯著的增大,增幅明顯大于前人在平地試驗中發(fā)現(xiàn)的結(jié)果,這種現(xiàn)象在前人研究中未有涉及,可能是導(dǎo)致臺風(fēng)路徑預(yù)報誤差的一個原因,值得更深入的研究和探討。

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

    Adem J,Lezama P,1960.On the motion of cyclone embedded in a uniform flow[J].Tellus,12:255258.

    AuYeung A Y M,Chan J C L,2010.The effect of a river delta and coastal roughness variation on a landfalling tropical cyclone[J].J Geophys Res,115(D19).

    Bender M A,Tuleya R E,Kurihara Y,1987.A numerical study of the island terrain on tropical cyclones[J].Mon Wea Rev,115:130155.

    曹劍,吳立廣,潘維玉,20122006年7—9月西北太平洋熱帶氣旋季節(jié)活動的數(shù)值模擬[J].大氣科學(xué)學(xué)報,35(2):148162.Cao J,Wu L G,Pan W Y,2012.Simulated seasonal activity of tropical cyclone over the western North Pacific during JulySeptember 2006[J].Trans Atmos Sci,35(2):148162.(in Chinese).

    Chan J C L,Gray W M,1982.Tropical cyclone motion and surrounding flow relationship[J].Mon Wea Rev,110:13541374.

    Chan J C L,Williams R T,1987.Analytical and numerical studies of betaeffect in tropical cyclone motion.Part Ⅰ:zero mean flow[J].J Atmos Sci,44:12571265.

    Chan J C L,Ko F M F,Lei Y M,2002.Relationship between potential vorticty tendency and tropical cyclone motion[J].J Atmos Sci,59:13171336.

    Chang S W J,1982.The orographic effects induced by an island mountain range on propagating tropical cyclones[J].Mon Wea Rev,110:12551268.

    Chen X,Wu L,Zhang J,2011.Increasing duration of tropical cyclones over China[J].Geophys Res Lett,38:L02708.doi:101029/2010GL046137.

    Chen Y,Yau M K,2003.Asymmetric structures in a simulated landfalling hurricane[J].J Atmos Sci,60(18):22942312.

    DeMaria M,1985.Tropical cyclone motion in a nondivergent barotropic model[J].Mon Wea Rev,113:11991210.

    Duchon C E,1979.Lanczos filtering in one and two dimensions[J].J Appl Meteor,18:10161022.

    Dudhia J,1989.Numerical study of convection observed during the winter monsoon experiment using a mesoscale twodimensional model[J].J Atmos Sci,46:30773107.

    Fiorino M J,Elsberry R L,1989.Some aspects of vortex structure related to tropical cyclone motion[J].J Atmos Sci,46:975990.

    Frank W M,Ritchie E A,2001.Effects of vertical wind shear on the intensity and structure of numerically simulated hurricanes[J].Mon Wea Rev,129:22492269.

    George J E,Gray W M,1976.Tropical cyclone motion and surrounding parameter relationships[J].J Appl Meteor,15:12521264.

    Holland G J,1983.Tropical cyclone motion:environmental interaction plus a beta effect[J].J Atmos Sci,40:328342.

    Hong S Y,Lim J O J,2006.The WRF singlemoment 6class microphysics scheme (WSM6)[J].J Korean Meteor Soc,42:129151.

    Huang Y H,Wu C C,Wang Y,2011.The influence of island topography on typhoon track deflection[J].Mon Wea Rev,139:17081727.

    Jian G J,Wu C C,2008.A numerical study of the track deflection of super typhoon Haitang (2005) prior to landfall in Taiwan[J].Mon Wea Rev,136:598615.

    Kain J S,F(xiàn)ritch J M,1993.Convective parameterization for mesoscale models:the KainFritch scheme.The representation of cumulus convection in numerical models Meteor Monogr[J].Amer Meteor Soc,46:165170.

    Li X,Wang B,1994.Acceleration of the hurricane beta drift by shear strain rate of an environmental flow[J].J Atmos Sci,53:327334.

    李肖雅,吳立廣,宗慧君,2014.季風(fēng)渦旋影響西北太平洋臺風(fēng)生成初步分析[J].大氣科學(xué)學(xué)報,37(5):653664.Li X Y,Wu L G,Zong H J,2014.Analysis of influence of monsoon gyres on tropical cyclogenesis over the western North Pacific[J].Trans Atmos Sci,37(5):653664.(in Chinese).

    Li Y,Cheung K K W,Chan J C L,2013.Numerical study on the development of asymmetric convection and vertical wind shear during tropical cyclone landfall[J].Quart J Roy Meteor Soc,140:18661877.

    Lin Y L,Chen S Y,Hill C M,et al.,2005.Control parameters for the influence of a mesoscale mountain range on cyclone track continuity and deflection[J].J Atmos Sci,62:18491866.

    Mlawer E J,Taubman S J,Brown P D,1997.Radiative transfer for inhomogeneous atmosphere:RRTM,a validated correlatedk model for the longwave[J].J Geophys Res,102(D14):1666316682.

    Neumann C J,1979.A guide to Atlantic and Eastern Pacific models for the prediction of tropical cyclone motion[C]//NOAA Tech Memo NWS NHC11.

    Noh Y,Cheon W G,Hong S Y,2003.Improvement of the Kprofile model for the planetary boundary layer based on largeeddy simulation data[J].Bound Layer Meteor,107:401427.

    Peng M S,William R T,1990.Dynamics of vortex asymmetries and their influence on vortex motion on a betaplane[J].J Atmos Sci,47:19872003.

    Szeto K C,Chan J C L,2010.Structural changes of a tropical cyclone during landfall:βplane simulations[J].Adv Atmos Sci,27:11431150.

    Tuleya R E,Kurihara Y,1978.A numerical simulation of the landfall of tropical cyclones[J].J Atmos Sci,35:242257.

    Wang B,Li X,Wu L,1997.Hurricane beta drift direction in horizontally sheared flows[J].J Atmos Sci,54:14621471.

    Wang C C,Kuo H C,Chen Y H,et al.,2012.Effects of asymmetric latent heating on typhoon movement crossing Taiwan:the case of Morakot (2009) with extreme rainfall[J].J Atmos Sci,69:31723196.

    Willoughby H E,1990.Temporal changes of the primary circulation in tropical cyclones[J].J Atmos Sci,47:242264.

    Willoughby H E,1992.Linear motion of a shallow water barotropic vortex as an initialvalue problem[J].J Atmos Sci,47:242264.

    Wong M L M,Chan J C L,2006.Modeling the effects of landsea roughness contrast on tropical cyclone winds[J].J Atmos Sci,64:32493264.

    Wu L,Wang B,2000.A potential vorticity tendency diagnostic approach for tropical cyclone motion[J].Mon Wea Rev,128:18991911.

    Wu L,Wang B,2001a.Movement and vertical coupling of adiabatic baroclinic tropical cyclones[J].J Atmos Sci,58:18011814.

    Wu L,Wang B,2001b.Effects of convective heating on movement and vertical coupling of tropical cyclones:a numerical study[J].J Atmos Sci,58:36393649.

    Wu L,Braun S A,Halverson J,et al.,2006.A numerical study of Hurricane Erin (2001).Part I:model verification and storm evolution[J].J Atmos Sci,63:6586.

    吳玉琴,吳立廣,梁佳,2015.利用位渦趨勢法診斷臺風(fēng)“莫拉克”(2009)的運動[J].氣象學(xué)報,73(2):236246.Wu Y Q,Wu L G,Liang J,2015.Diagnosi of the movement of Typhoon Morakot (2009) using the potential vorticity tendency approach[J].Acta Meteor Sinica,73(2):236246.(in Chinese).

    Xie B,Zhang F,2012.Impacts of typhoon track and island topography on the heavy rainfalls in Taiwan associated with Morakot (2009)[J].Mon Wea Rev,140:33793394.

    Yang M J,Zhang D L,Huang H L,2008.A modeling study of typhoon Nari (2001) at landfall.Part I:topography effects[J].J Atmos Sci,65:30953115.

    楊瓊瓊,吳立廣,2015.中國24 h臺風(fēng)路徑預(yù)報難點及其大尺度環(huán)流分析[J].氣象科學(xué),35(4):389397.Yang Q Q,Wu L G,2015.Analysis on difficulty of 24 h typhoon track forecast in China and associated largescale circulation[J].J Meteor Sci,35(4):389397.(in Chinese).

    Yeh T C,Elsberry R L,1993a.Interaction of typhoons with the orography.Part I:upstream track deflections[J].Mon Wea Rev,121:31933212.

    Yeh T C,Elsberry R L,1993b.Interaction of typhoons with the orography.Part II:continuous and discontinuous tracks across the island[J].Mon Wea Rev,121:32133233.

    張嬌艷,吳立廣,張強(qiáng),2011.全球變暖背景下我國熱帶氣旋災(zāi)害趨勢分析[J].熱帶氣象學(xué)報,27(4):442454.Zhang J Y,Wu L G,Zhang Q,2011.Tropical cyclone damages in China under the background of global warming[J].J Tropical Meteor,27(4):442454.(in Chinese).

    Zhang Q,Liu Q,Wu L,2009.Tropical cyclone damages in China 1983—2006[J].Bull Amer Meteor Soc,90:489495.

    Idealized numerical study on influence of topography along East China coast on landfall tropical cyclone motion

    日韩免费高清中文字幕av| 侵犯人妻中文字幕一二三四区| 国产视频一区二区在线看| 久久久国产欧美日韩av| 啦啦啦 在线观看视频| 久久久亚洲精品成人影院| 首页视频小说图片口味搜索 | 蜜桃在线观看..| 午夜福利在线免费观看网站| 99久久精品国产亚洲精品| 亚洲成av片中文字幕在线观看| 国产亚洲精品久久久久5区| 国产成人av教育| 脱女人内裤的视频| 一区二区三区激情视频| 高清视频免费观看一区二区| 老鸭窝网址在线观看| 久久久久久亚洲精品国产蜜桃av| 王馨瑶露胸无遮挡在线观看| 成人国产一区最新在线观看 | 啦啦啦视频在线资源免费观看| 咕卡用的链子| 9色porny在线观看| 国产精品成人在线| 日本色播在线视频| 精品一区二区三区四区五区乱码 | 91精品国产国语对白视频| 亚洲熟女精品中文字幕| av福利片在线| √禁漫天堂资源中文www| 黄片播放在线免费| 婷婷色综合www| 咕卡用的链子| 在现免费观看毛片| 一区二区三区激情视频| 极品少妇高潮喷水抽搐| 亚洲人成网站在线观看播放| 黄色一级大片看看| 老鸭窝网址在线观看| 亚洲激情五月婷婷啪啪| 国产免费视频播放在线视频| 又大又黄又爽视频免费| 极品人妻少妇av视频| bbb黄色大片| 欧美日韩一级在线毛片| 国产亚洲欧美精品永久| 国产欧美日韩综合在线一区二区| 男人操女人黄网站| 成年人黄色毛片网站| 少妇人妻久久综合中文| 日本色播在线视频| 亚洲精品久久久久久婷婷小说| 一本—道久久a久久精品蜜桃钙片| 最新在线观看一区二区三区 | 亚洲欧美一区二区三区久久| 超碰97精品在线观看| 人人妻,人人澡人人爽秒播 | 狠狠婷婷综合久久久久久88av| 亚洲九九香蕉| 免费人妻精品一区二区三区视频| 免费看十八禁软件| 少妇人妻 视频| 人人妻人人澡人人看| 亚洲精品乱久久久久久| 麻豆国产av国片精品| 婷婷丁香在线五月| 校园人妻丝袜中文字幕| 黑人欧美特级aaaaaa片| 人人妻人人添人人爽欧美一区卜| 午夜91福利影院| 热99国产精品久久久久久7| 免费在线观看视频国产中文字幕亚洲 | 亚洲九九香蕉| 最近中文字幕2019免费版| 男的添女的下面高潮视频| 交换朋友夫妻互换小说| 亚洲av国产av综合av卡| 最黄视频免费看| 女警被强在线播放| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美一区二区三区久久| 一区二区三区激情视频| 搡老岳熟女国产| 国产精品麻豆人妻色哟哟久久| 日韩欧美一区视频在线观看| 欧美人与善性xxx| 叶爱在线成人免费视频播放| 最近最新中文字幕大全免费视频 | www日本在线高清视频| 麻豆乱淫一区二区| 欧美亚洲日本最大视频资源| 亚洲av在线观看美女高潮| 午夜福利乱码中文字幕| 少妇被粗大的猛进出69影院| 永久免费av网站大全| 精品熟女少妇八av免费久了| 免费久久久久久久精品成人欧美视频| 少妇的丰满在线观看| 精品高清国产在线一区| 一级a爱视频在线免费观看| 久久精品国产a三级三级三级| 丝袜喷水一区| 亚洲男人天堂网一区| 国产黄色免费在线视频| 久久精品久久久久久噜噜老黄| 欧美中文综合在线视频| 黄色怎么调成土黄色| 在现免费观看毛片| 男人操女人黄网站| 亚洲av成人精品一二三区| 日韩av免费高清视频| 免费一级毛片在线播放高清视频 | 国产有黄有色有爽视频| 久久99精品国语久久久| 色婷婷av一区二区三区视频| 国产欧美日韩综合在线一区二区| 亚洲精品国产色婷婷电影| 人人妻人人添人人爽欧美一区卜| 午夜福利视频精品| 下体分泌物呈黄色| www.自偷自拍.com| 美女中出高潮动态图| 男女无遮挡免费网站观看| 老熟女久久久| 青春草视频在线免费观看| 国产爽快片一区二区三区| 99久久人妻综合| 天天影视国产精品| 老司机影院毛片| 嫁个100分男人电影在线观看 | 性色av一级| 韩国精品一区二区三区| 少妇被粗大的猛进出69影院| 大陆偷拍与自拍| 大香蕉久久成人网| 七月丁香在线播放| xxxhd国产人妻xxx| 久久久久网色| 999久久久国产精品视频| 少妇被粗大的猛进出69影院| 国产伦人伦偷精品视频| 欧美日韩一级在线毛片| 中国国产av一级| a 毛片基地| 在线亚洲精品国产二区图片欧美| 亚洲av日韩在线播放| av一本久久久久| 999精品在线视频| 日韩熟女老妇一区二区性免费视频| 日韩一卡2卡3卡4卡2021年| 美国免费a级毛片| 成年动漫av网址| 国产色视频综合| av国产久精品久网站免费入址| 美女午夜性视频免费| 国产精品国产三级专区第一集| 亚洲,一卡二卡三卡| 亚洲精品国产av蜜桃| 欧美日韩亚洲综合一区二区三区_| svipshipincom国产片| 婷婷成人精品国产| 精品国产乱码久久久久久男人| 日本猛色少妇xxxxx猛交久久| 韩国高清视频一区二区三区| 成人亚洲精品一区在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩精品网址| 波野结衣二区三区在线| 久久国产精品人妻蜜桃| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 国产精品久久久久成人av| 国产日韩欧美在线精品| 99精品久久久久人妻精品| 亚洲欧美清纯卡通| 99国产精品免费福利视频| 女警被强在线播放| 一二三四在线观看免费中文在| 看免费成人av毛片| 可以免费在线观看a视频的电影网站| 久久国产精品影院| 狠狠婷婷综合久久久久久88av| 亚洲国产欧美在线一区| 男人舔女人的私密视频| 丝袜人妻中文字幕| 老司机影院成人| 亚洲久久久国产精品| 精品国产一区二区久久| 丝袜脚勾引网站| 欧美黄色片欧美黄色片| 精品欧美一区二区三区在线| 91字幕亚洲| 国产成人精品久久二区二区91| 男女高潮啪啪啪动态图| 悠悠久久av| 男女免费视频国产| 黄色a级毛片大全视频| 极品少妇高潮喷水抽搐| 69精品国产乱码久久久| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜一区二区| 国产精品三级大全| 精品人妻一区二区三区麻豆| 麻豆国产av国片精品| 99国产精品一区二区三区| 午夜福利,免费看| 精品国产一区二区久久| 国产黄色视频一区二区在线观看| 一级毛片黄色毛片免费观看视频| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 在线观看免费日韩欧美大片| 欧美国产精品va在线观看不卡| 91精品三级在线观看| 最新在线观看一区二区三区 | 国产精品久久久久久人妻精品电影 | 免费高清在线观看视频在线观看| 在现免费观看毛片| 欧美日韩综合久久久久久| 伦理电影免费视频| 91成人精品电影| 这个男人来自地球电影免费观看| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一二三| 七月丁香在线播放| 久久人人97超碰香蕉20202| 天天操日日干夜夜撸| 91老司机精品| 亚洲国产看品久久| 七月丁香在线播放| 欧美+亚洲+日韩+国产| 久久久久久久久久久久大奶| 好男人视频免费观看在线| kizo精华| 可以免费在线观看a视频的电影网站| 亚洲精品中文字幕在线视频| 男的添女的下面高潮视频| 91精品伊人久久大香线蕉| 亚洲精品久久成人aⅴ小说| 黄色片一级片一级黄色片| 一级黄色大片毛片| 最近最新中文字幕大全免费视频 | 久久99一区二区三区| 国产欧美日韩一区二区三 | 亚洲精品国产一区二区精华液| 欧美成人精品欧美一级黄| 纵有疾风起免费观看全集完整版| 国产1区2区3区精品| 人妻 亚洲 视频| 69精品国产乱码久久久| 老司机在亚洲福利影院| 国产免费福利视频在线观看| 在线观看免费高清a一片| 成人午夜精彩视频在线观看| 国产成人91sexporn| 激情视频va一区二区三区| 高清黄色对白视频在线免费看| 午夜免费鲁丝| 精品国产一区二区三区四区第35| 欧美黄色淫秽网站| 丝瓜视频免费看黄片| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区国产| 亚洲七黄色美女视频| 国产1区2区3区精品| 宅男免费午夜| 丰满人妻熟妇乱又伦精品不卡| 在线观看人妻少妇| 99久久精品国产亚洲精品| 丝袜在线中文字幕| 亚洲一区二区三区欧美精品| 欧美国产精品va在线观看不卡| 亚洲av日韩精品久久久久久密 | 亚洲国产看品久久| 国产一区亚洲一区在线观看| 欧美+亚洲+日韩+国产| 成人手机av| 夜夜骑夜夜射夜夜干| 高清视频免费观看一区二区| √禁漫天堂资源中文www| 亚洲图色成人| 国产欧美日韩综合在线一区二区| 在线观看www视频免费| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜爱| 亚洲,欧美精品.| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 精品一区二区三区四区五区乱码 | 色精品久久人妻99蜜桃| 成年人午夜在线观看视频| 九草在线视频观看| 脱女人内裤的视频| 久久精品国产亚洲av涩爱| 激情五月婷婷亚洲| 精品人妻在线不人妻| 久久久精品国产亚洲av高清涩受| 天天躁夜夜躁狠狠躁躁| 欧美日韩一级在线毛片| 午夜视频精品福利| 亚洲天堂av无毛| 9热在线视频观看99| 国产有黄有色有爽视频| 女人爽到高潮嗷嗷叫在线视频| 日韩,欧美,国产一区二区三区| 免费av中文字幕在线| 极品少妇高潮喷水抽搐| 99久久精品国产亚洲精品| 婷婷成人精品国产| 国产黄色免费在线视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲色图 男人天堂 中文字幕| 777久久人妻少妇嫩草av网站| 国产精品av久久久久免费| av国产精品久久久久影院| 国产主播在线观看一区二区 | 国产伦理片在线播放av一区| 青春草亚洲视频在线观看| avwww免费| 国产精品人妻久久久影院| av在线app专区| 久久性视频一级片| 婷婷色av中文字幕| 男人操女人黄网站| 性色av一级| av网站免费在线观看视频| 午夜91福利影院| 在线精品无人区一区二区三| 亚洲国产精品一区三区| 久久天堂一区二区三区四区| 国产视频首页在线观看| 蜜桃在线观看..| 啦啦啦在线观看免费高清www| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频 | 日本黄色日本黄色录像| 亚洲国产精品一区三区| 国产又色又爽无遮挡免| 国产高清国产精品国产三级| 另类亚洲欧美激情| 男女高潮啪啪啪动态图| 80岁老熟妇乱子伦牲交| 欧美另类一区| 亚洲视频免费观看视频| 国产高清国产精品国产三级| 国产午夜精品一二区理论片| 久久午夜综合久久蜜桃| 久久久欧美国产精品| svipshipincom国产片| av欧美777| av网站在线播放免费| 欧美在线黄色| √禁漫天堂资源中文www| 搡老岳熟女国产| 狂野欧美激情性xxxx| 无遮挡黄片免费观看| 久久久久久亚洲精品国产蜜桃av| 国产精品麻豆人妻色哟哟久久| 日本wwww免费看| 午夜精品国产一区二区电影| 欧美日韩福利视频一区二区| 丝瓜视频免费看黄片| 亚洲av成人精品一二三区| 亚洲 国产 在线| 18禁黄网站禁片午夜丰满| 亚洲国产成人一精品久久久| 18禁黄网站禁片午夜丰满| 伊人久久大香线蕉亚洲五| 久久性视频一级片| 午夜视频精品福利| 欧美人与性动交α欧美精品济南到| 美女视频免费永久观看网站| 久久九九热精品免费| 亚洲中文字幕日韩| 婷婷色综合大香蕉| 国产在线一区二区三区精| 国产精品秋霞免费鲁丝片| 黄色 视频免费看| 69精品国产乱码久久久| 欧美日韩亚洲高清精品| 日本wwww免费看| 亚洲精品一区蜜桃| 可以免费在线观看a视频的电影网站| 国产一区有黄有色的免费视频| 国产精品秋霞免费鲁丝片| 中文字幕亚洲精品专区| 亚洲av成人精品一二三区| 亚洲精品一区蜜桃| av又黄又爽大尺度在线免费看| 久久国产精品男人的天堂亚洲| xxx大片免费视频| 男女无遮挡免费网站观看| 色综合欧美亚洲国产小说| 精品国产一区二区三区四区第35| 一级毛片电影观看| 亚洲,欧美,日韩| 免费在线观看日本一区| 波野结衣二区三区在线| 欧美 亚洲 国产 日韩一| 色网站视频免费| 无限看片的www在线观看| 国产成人精品久久久久久| 美女国产高潮福利片在线看| 欧美精品一区二区大全| 亚洲成人免费电影在线观看 | 91成人精品电影| av片东京热男人的天堂| 精品国产一区二区三区四区第35| 国产高清videossex| 久久亚洲国产成人精品v| 91老司机精品| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 国产三级黄色录像| 天天躁日日躁夜夜躁夜夜| 欧美xxⅹ黑人| 亚洲精品第二区| 欧美精品一区二区免费开放| 老司机靠b影院| 国产国语露脸激情在线看| 亚洲专区中文字幕在线| 久久久久久亚洲精品国产蜜桃av| 免费少妇av软件| 一边亲一边摸免费视频| 蜜桃在线观看..| 亚洲国产精品一区三区| 久久久久精品人妻al黑| 亚洲国产最新在线播放| 国产成人a∨麻豆精品| 岛国毛片在线播放| 日韩免费高清中文字幕av| 在线精品无人区一区二区三| 啦啦啦中文免费视频观看日本| 看十八女毛片水多多多| 久久精品亚洲熟妇少妇任你| 男女边摸边吃奶| 国产无遮挡羞羞视频在线观看| 制服诱惑二区| 国产高清视频在线播放一区 | netflix在线观看网站| 欧美激情高清一区二区三区| 一边摸一边做爽爽视频免费| 国产三级黄色录像| 丰满少妇做爰视频| 国产精品久久久久久人妻精品电影 | 久久久久精品国产欧美久久久 | 中文字幕人妻丝袜制服| 欧美国产精品va在线观看不卡| 日韩av不卡免费在线播放| 免费看av在线观看网站| 日韩制服骚丝袜av| 美女午夜性视频免费| 欧美精品av麻豆av| 亚洲中文av在线| 色播在线永久视频| 满18在线观看网站| 侵犯人妻中文字幕一二三四区| 男人操女人黄网站| 99久久人妻综合| 中文精品一卡2卡3卡4更新| 久久天堂一区二区三区四区| 在线观看免费高清a一片| 成年美女黄网站色视频大全免费| 国产视频一区二区在线看| av福利片在线| 91成人精品电影| 又粗又硬又长又爽又黄的视频| 午夜免费男女啪啪视频观看| 啦啦啦在线免费观看视频4| 99国产综合亚洲精品| 一本大道久久a久久精品| 高清黄色对白视频在线免费看| 国产黄色视频一区二区在线观看| 中文字幕人妻熟女乱码| 天天影视国产精品| 成人影院久久| 看免费成人av毛片| 免费日韩欧美在线观看| 交换朋友夫妻互换小说| 一区二区三区激情视频| 99久久精品国产亚洲精品| 久久亚洲精品不卡| 啦啦啦在线免费观看视频4| 肉色欧美久久久久久久蜜桃| 色婷婷久久久亚洲欧美| 18在线观看网站| 悠悠久久av| 各种免费的搞黄视频| 欧美av亚洲av综合av国产av| 久久国产精品人妻蜜桃| 人人妻人人澡人人爽人人夜夜| 亚洲国产欧美网| 亚洲黑人精品在线| 国产成人欧美| 中文字幕色久视频| 午夜福利,免费看| 国产1区2区3区精品| 在线看a的网站| 亚洲欧美精品综合一区二区三区| 日本一区二区免费在线视频| 精品国产一区二区三区四区第35| 黄色怎么调成土黄色| 看十八女毛片水多多多| a 毛片基地| 欧美精品一区二区免费开放| 人人妻,人人澡人人爽秒播 | 国产真人三级小视频在线观看| 亚洲情色 制服丝袜| 欧美激情极品国产一区二区三区| 亚洲精品日本国产第一区| 亚洲美女黄色视频免费看| 丝瓜视频免费看黄片| 免费久久久久久久精品成人欧美视频| 日本黄色日本黄色录像| 婷婷色综合www| 只有这里有精品99| 91成人精品电影| 黑人猛操日本美女一级片| 老汉色av国产亚洲站长工具| 国产精品久久久久成人av| 电影成人av| 国产爽快片一区二区三区| 久久久久国产一级毛片高清牌| 人人妻人人爽人人添夜夜欢视频| 9191精品国产免费久久| 国产黄色免费在线视频| 亚洲少妇的诱惑av| av一本久久久久| 亚洲成色77777| 亚洲欧美一区二区三区国产| 精品国产一区二区久久| 日韩精品免费视频一区二区三区| 国产av国产精品国产| 欧美在线一区亚洲| 999久久久国产精品视频| 日本猛色少妇xxxxx猛交久久| av有码第一页| 亚洲午夜精品一区,二区,三区| 中文字幕人妻熟女乱码| h视频一区二区三区| 亚洲精品久久午夜乱码| 久久99一区二区三区| 久久久久久免费高清国产稀缺| 日韩大片免费观看网站| 亚洲久久久国产精品| 国语对白做爰xxxⅹ性视频网站| 久久99精品国语久久久| 亚洲国产毛片av蜜桃av| 欧美久久黑人一区二区| 亚洲欧美一区二区三区黑人| 晚上一个人看的免费电影| 熟女少妇亚洲综合色aaa.| 国产精品.久久久| 一级黄色大片毛片| 桃花免费在线播放| 亚洲人成网站在线观看播放| 一级毛片我不卡| 欧美精品一区二区免费开放| 亚洲少妇的诱惑av| 欧美另类一区| 人体艺术视频欧美日本| 妹子高潮喷水视频| 18禁裸乳无遮挡动漫免费视频| 国产一卡二卡三卡精品| 老司机深夜福利视频在线观看 | 一区二区三区四区激情视频| 欧美精品人与动牲交sv欧美| 在线亚洲精品国产二区图片欧美| 国产精品 欧美亚洲| 亚洲av电影在线进入| 国产日韩欧美视频二区| 日本一区二区免费在线视频| 午夜免费男女啪啪视频观看| 午夜福利影视在线免费观看| 欧美亚洲 丝袜 人妻 在线| 亚洲精品一区蜜桃| 黄片播放在线免费| 肉色欧美久久久久久久蜜桃| 在线观看www视频免费| 日韩大片免费观看网站| 国产一区二区三区综合在线观看| 欧美日韩精品网址| av片东京热男人的天堂| 国产欧美日韩精品亚洲av| 老司机在亚洲福利影院| 国产成人91sexporn| 国产精品二区激情视频| 欧美日韩精品网址| 国产成人免费观看mmmm| 国产在线视频一区二区| 美女国产高潮福利片在线看| 精品国产乱码久久久久久小说| avwww免费| 欧美日韩国产mv在线观看视频| 这个男人来自地球电影免费观看| 波多野结衣av一区二区av| 国精品久久久久久国模美| 久久天躁狠狠躁夜夜2o2o | 青春草亚洲视频在线观看| 久久精品国产亚洲av高清一级| 伊人久久大香线蕉亚洲五| 久久国产精品人妻蜜桃| 国产成人精品无人区| 亚洲精品国产av蜜桃| 在线观看国产h片| 国产精品av久久久久免费| 每晚都被弄得嗷嗷叫到高潮| av一本久久久久| 精品免费久久久久久久清纯 | 亚洲精品国产一区二区精华液| 国产精品一区二区在线观看99| 亚洲精品第二区| 国产精品人妻久久久影院| 国产免费现黄频在线看| 又紧又爽又黄一区二区|