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

    FY-2D云導(dǎo)風(fēng)在“鯰魚(yú)”臺(tái)風(fēng)模擬中的同化試驗(yàn)

    2011-06-29 06:14:30朱克云
    關(guān)鍵詞:冷啟動(dòng)鯰魚(yú)實(shí)況

    蘇 捷, 朱克云, 張 杰

    (1.成都信息工程學(xué)院大氣科學(xué)學(xué)院,四川成都610225;2.成都軍區(qū)空軍氣象中心,四川成都610041;3.中國(guó)人民解放軍96115部隊(duì),吉林通化134001)

    1 引言

    隨著近幾年大氣觀測(cè)技術(shù)的逐步提高,衛(wèi)星、雷達(dá)等一系列高分辨率的非常規(guī)觀測(cè)已在預(yù)報(bào)和數(shù)值模擬等方面取得了較為廣泛應(yīng)用,成效顯著,很大程度的彌補(bǔ)了常規(guī)觀測(cè)資料時(shí)空分布不足的缺點(diǎn)。對(duì)于觀測(cè)稀疏的海洋以及高原地區(qū),云導(dǎo)風(fēng)資料可以對(duì)數(shù)值預(yù)報(bào)的初始場(chǎng)做重要補(bǔ)充和改善。

    早在20世紀(jì)80年代的美國(guó)衛(wèi)星應(yīng)用會(huì)議上,Purdom et al.[1]就提出可以通過(guò)追蹤衛(wèi)星照片中云的移動(dòng)和發(fā)展,來(lái)觀測(cè)分析風(fēng)暴等中尺度系統(tǒng)的發(fā)展和演變。之后隨著新型號(hào)衛(wèi)星的發(fā)射及新探測(cè)器的出現(xiàn),云導(dǎo)風(fēng)在數(shù)值模擬中的使用和研究也越發(fā)廣泛,John F.Le Marshall et al.[2]通過(guò)對(duì)比了3種不同的資料同化方法的結(jié)果之后,得出:變分方法能更好的利用高時(shí)空分辨率的云導(dǎo)風(fēng)資料中所含的信息。R.C.Bhatia et al.[3]通過(guò)在有限區(qū)域模式中同化METEOSAT-5水汽通道的云導(dǎo)風(fēng)數(shù)據(jù),得出加入云導(dǎo)風(fēng)可以對(duì)熱帶氣旋的路徑預(yù)測(cè)產(chǎn)生了正效果。B.J.Soden et al.[4]通過(guò)在GFDL颶風(fēng)模式中同化GOES-8的云導(dǎo)風(fēng)資料減小了對(duì)颶風(fēng)路徑預(yù)報(bào)的偏西誤差,得出同化云導(dǎo)風(fēng)可以較好的優(yōu)化初始場(chǎng),減小預(yù)報(bào)誤差。云導(dǎo)風(fēng)對(duì)臺(tái)風(fēng)路徑的影響在國(guó)內(nèi)也有較多的研究和探索,張守峰等[5,6]使用GMS-5衛(wèi)星云導(dǎo)風(fēng)資料用最有插值法訂正客觀分析風(fēng)場(chǎng),發(fā)現(xiàn)新的分析場(chǎng)對(duì)臺(tái)風(fēng)路徑預(yù)報(bào)有明顯改進(jìn)。王棟梁等[7]通過(guò)一系列四維變分同化試驗(yàn)同化GMS-5水汽和紅外云導(dǎo)風(fēng)資料,結(jié)果表明同化之后對(duì)熱帶氣旋路徑預(yù)報(bào)有一定改善。莊照榮等[8]通過(guò)使用GRAPES同化GMS-5衛(wèi)星云導(dǎo)風(fēng)資料對(duì)“森拉克”臺(tái)風(fēng)做24小時(shí)預(yù)報(bào),得到了質(zhì)量更高的風(fēng)壓場(chǎng)和更準(zhǔn)確的預(yù)報(bào)路徑。何斌等[9]以同化試驗(yàn)和無(wú)同化控制試驗(yàn)的檢驗(yàn)對(duì)比,得出云導(dǎo)風(fēng)資料對(duì)u,v風(fēng)場(chǎng)的改進(jìn)主要集中在對(duì)流層高層,對(duì)中低層的影響不很明顯。我國(guó)自FY-2C靜止衛(wèi)星的發(fā)射至退役(現(xiàn)在被FY-2E衛(wèi)星的成功替代),其云導(dǎo)風(fēng)產(chǎn)品在數(shù)值模式中的使用大多用于對(duì)暴雨的模擬預(yù)報(bào)或數(shù)據(jù)的統(tǒng)計(jì)處理,而在臺(tái)風(fēng)預(yù)報(bào)方面并沒(méi)有進(jìn)行太多研究。如李華宏等[10]同化FY-2C云導(dǎo)風(fēng)資料發(fā)現(xiàn)在暴雨區(qū)附近初始風(fēng)場(chǎng)低層的西南氣流明顯加強(qiáng),能很好地改善24h暴雨預(yù)報(bào)的強(qiáng)度和落區(qū)。其他分析[11-13]也證實(shí)在暴雨預(yù)報(bào)業(yè)務(wù)和科研中,云導(dǎo)風(fēng)資料的使用對(duì)提高降水預(yù)報(bào)質(zhì)量是一個(gè)有效途徑。何志新等[14]在一個(gè)個(gè)例中同化了一個(gè)時(shí)次的FY-2C云導(dǎo)風(fēng)資料,得出其對(duì)臺(tái)風(fēng)的降水預(yù)報(bào)有較好改善。

    鑒于FY-2C衛(wèi)星退役之后,國(guó)內(nèi)對(duì)我國(guó)云導(dǎo)風(fēng)產(chǎn)品研究和應(yīng)用較為少缺,因此以不同于前人的研究思路和方法,對(duì)FY-2D導(dǎo)風(fēng)產(chǎn)品進(jìn)入WRF-3DVAR模塊的同化結(jié)果進(jìn)行探討,并分析不同通道的產(chǎn)品和不同同化方法對(duì)臺(tái)風(fēng)預(yù)報(bào)的可能產(chǎn)生的影響。

    2 云導(dǎo)風(fēng)資料來(lái)源及其預(yù)處理

    利用靜止氣象衛(wèi)星紅外云圖和水汽圖像資料,追蹤連續(xù)的幾幅云圖上的目標(biāo)云塊(示蹤云)的位移,估計(jì)出大氣的運(yùn)動(dòng)并利用球面三角公式推算出平均運(yùn)動(dòng)矢量,作為該區(qū)域示蹤云所在高度的風(fēng)的近似,這種運(yùn)動(dòng)矢量稱(chēng)為大氣運(yùn)動(dòng)矢量(Atmosphere Movtion Vector,AMV),即國(guó)內(nèi)所稱(chēng)的云導(dǎo)風(fēng)(云跡風(fēng))。

    所使用的云導(dǎo)風(fēng)資料來(lái)自國(guó)家衛(wèi)星中心的FY-2D大氣運(yùn)動(dòng)矢量產(chǎn)品,其等壓面高度是根據(jù)云塊的環(huán)境溫度推算得出。產(chǎn)品在生成過(guò)程中,國(guó)家衛(wèi)星中心產(chǎn)品開(kāi)發(fā)人員首先對(duì)接收到的數(shù)據(jù)進(jìn)行了嚴(yán)格的質(zhì)量檢驗(yàn)和容錯(cuò)處理,其后又應(yīng)用雙通道法對(duì)矢量高度進(jìn)行估計(jì),最后再經(jīng)過(guò)空間連續(xù)性,時(shí)間連續(xù)性,地物剔出等一系列質(zhì)量控制的調(diào)整[15]。因此,本文所用云導(dǎo)風(fēng)資料已經(jīng)通過(guò)質(zhì)量檢驗(yàn),數(shù)據(jù)真實(shí)可用。

    FY-2D 云導(dǎo)風(fēng)數(shù)據(jù)的計(jì)算范圍從 37°E-137°E,50°N-50°S,即 100°×100°,分辨率為 1°×1°。將下載的2010年10月18日到2010年10月23日的導(dǎo)風(fēng)數(shù)據(jù)解碼并提取,由于云導(dǎo)風(fēng)是依據(jù)云頂?shù)囊苿?dòng)確立,因此對(duì)該時(shí)間段內(nèi)所有導(dǎo)風(fēng)數(shù)據(jù)在高度上的分布情況進(jìn)行統(tǒng)計(jì),結(jié)果如表1所示。

    表1 云導(dǎo)風(fēng)產(chǎn)品的高度分布(單位:hPa)

    可見(jiàn)云導(dǎo)風(fēng)數(shù)據(jù)云跡風(fēng)垂直分布比較復(fù)雜,且大多不在標(biāo)準(zhǔn)的層次上,并且高層的導(dǎo)風(fēng)資料較多主要集中在500hPa以上,其中紅外通道70%,水汽通道95%,而對(duì)于中、低層資料較少,且水汽通道遠(yuǎn)少于紅外通道。

    為了避免或減小相鄰數(shù)據(jù)間的相關(guān)性,在進(jìn)行資料的同化之前需要將分布密集的數(shù)據(jù)稀疏化處理。處理方法:把處于標(biāo)準(zhǔn)氣壓層(包括:200hPa,500hPa,700hPa,850hPa)±50hPa范圍的所有數(shù)據(jù)的均值作為該標(biāo)準(zhǔn)氣壓層上的數(shù)據(jù)值。稀疏化處理之后的數(shù)據(jù)分布如圖1所示。

    3 試驗(yàn)的設(shè)計(jì)和方法

    3.1 臺(tái)風(fēng)個(gè)例的介紹

    1013號(hào)超強(qiáng)臺(tái)風(fēng)鲇魚(yú)是近20年來(lái),巔峰時(shí)刻中心氣壓低于890hPa的唯一一個(gè)臺(tái)風(fēng),“鲇魚(yú)”最顯著的特征就是中心附近的風(fēng)力很大,在10月17日晚上21點(diǎn)左右,美軍的探測(cè)飛機(jī)穿越臺(tái)風(fēng)鲇魚(yú)的中心進(jìn)行探測(cè)。落海儀器測(cè)得890hPa的海平面氣壓,和85m/s的一分鐘平均風(fēng)速。18日凌晨開(kāi)始,鲇魚(yú)進(jìn)入巔峰狀態(tài),并在18日中午12:25左右,在菲律賓呂宋島伊麗莎白省登陸,這是4年來(lái)全世界范圍內(nèi)登陸的最強(qiáng)臺(tái)風(fēng)。10月23日12時(shí)55分在福建省漳浦縣沿海登陸,登陸時(shí)中心附近最大風(fēng)力有13級(jí)(38m/s),中心最低氣壓為970hPa,成為2010年西太平洋和南海海域生成的最強(qiáng)臺(tái)風(fēng)。

    3.2 模式的介紹

    所用模式為WRFV3.2,微物理過(guò)程采用WSM3方案,對(duì)流參數(shù)化方案選擇修改過(guò)的Kain-Fritsch方案(能夠較好的同時(shí)描述深對(duì)流和淺對(duì)流過(guò)程),試驗(yàn)中初始場(chǎng)和側(cè)邊界條件選用NCEP 1°×1°預(yù)報(bào)資料,每6h更新一次側(cè)邊界條件,模式水平分辨率30km,垂直方向分28層,積分步長(zhǎng)為180秒,資料3小時(shí)輸出一次,云導(dǎo)風(fēng)資料的同化在每天的00時(shí)和12時(shí),同化時(shí)間窗口取±3小時(shí)。

    3.3 試驗(yàn)的目的和方案

    首先,通過(guò)對(duì)導(dǎo)風(fēng)數(shù)據(jù)的同化前和同化后的比較,包括對(duì)兩種不同通道數(shù)據(jù)同化之后效果的比較,不同同化方案的結(jié)果對(duì)比,以及對(duì)比分析資料稀疏化處理的效果,探討FY-2D云導(dǎo)風(fēng)產(chǎn)品對(duì)臺(tái)風(fēng)路徑和物理量場(chǎng)的改善作用。

    圖1 稀疏化處理后數(shù)據(jù)分布

    其次,以往在WRF模式的應(yīng)用中,模式的向前積分往往只限于48小時(shí)以?xún)?nèi),以避免模式中誤差的累積增大,然而如果在進(jìn)行長(zhǎng)時(shí)間連續(xù)積分的同時(shí)結(jié)合循環(huán)同化,效果會(huì)更好還是更差呢,將通過(guò)延長(zhǎng)積分時(shí)間(延長(zhǎng)為120小時(shí))來(lái)對(duì)此進(jìn)行試驗(yàn)和探討。

    研究設(shè)定了兩個(gè)同化方案,4組試驗(yàn)對(duì)比:

    同化方案一:初始時(shí)刻進(jìn)行模式的冷啟動(dòng),之后采用每12小時(shí)一次的同化循環(huán),總共積分120小時(shí)。

    同化方案二:在方案一基礎(chǔ)上,每積分24小時(shí)進(jìn)行一次冷啟動(dòng)。

    試驗(yàn)設(shè)計(jì)如表2所示:

    表2 試驗(yàn)設(shè)計(jì)

    4 試驗(yàn)結(jié)果的對(duì)比和分析

    4.1 臺(tái)風(fēng)路徑的模擬對(duì)比

    由于云導(dǎo)風(fēng)資料多分布于200-300hPa左右的對(duì)流層中上部,低層資料比較稀少,為了更好的比較同化之后的模擬效果,因此所做研究均取模擬后的200hPa高度為研究對(duì)象,取200hPa的結(jié)果畫(huà)路徑圖,為了區(qū)分圖中不同的試驗(yàn)結(jié)果,每隔6小時(shí)以不同符號(hào)進(jìn)行標(biāo)記,如圖2所示。

    圖2 路徑圖:其中實(shí)況(實(shí)心圓),控制試驗(yàn)(空心圓),IR1同化(空心方),IR3同化(實(shí)心方)

    圖2(a)中包含實(shí)況,控制試驗(yàn),對(duì)比試驗(yàn)1和對(duì)比試驗(yàn)2;圖2(b)中包含實(shí)況,控制試驗(yàn)和對(duì)比試驗(yàn)3;圖2(c)中包含實(shí)況,控制試驗(yàn)和對(duì)比試驗(yàn)4。分別對(duì)圖2的(a),(b),(c)進(jìn)行分析可以分別得出3組結(jié)論。

    從圖2(a)中可以看出:(1)沒(méi)有進(jìn)行任何同化的控制試驗(yàn)所模擬的臺(tái)風(fēng)路徑與實(shí)況路徑相比,在前50小時(shí)二者基本吻合,隨著模擬時(shí)間的延長(zhǎng),控制試驗(yàn)臺(tái)風(fēng)中心漸漸向東偏離,并且移速增大,登陸較早。(2)根據(jù)鯰魚(yú)臺(tái)風(fēng)的實(shí)況觀測(cè),臺(tái)風(fēng)在19日12點(diǎn)到14點(diǎn)開(kāi)始北折(即圖中路徑上第7個(gè)點(diǎn)),此時(shí)對(duì)比試驗(yàn)1和對(duì)比試驗(yàn)2的結(jié)果與實(shí)況的偏差開(kāi)始增大,均明顯偏東,比控制試驗(yàn)出現(xiàn)明顯偏差的時(shí)間偏早。尤其是對(duì)比試驗(yàn)2,負(fù)效果遠(yuǎn)大于另外兩組實(shí)驗(yàn),然而同化了IR1-AMV的對(duì)比試驗(yàn)1在22日00時(shí)(即第16個(gè)點(diǎn))以后,臺(tái)風(fēng)的路徑相比控制試驗(yàn)更加貼近于實(shí)況,移速也比其他兩組實(shí)驗(yàn)更接近實(shí)況,可見(jiàn)同化在此時(shí)產(chǎn)生了較明顯的正效果。(3)以模擬48小時(shí)為界(即圖中路徑上第9個(gè)點(diǎn)),在此時(shí)刻之前3條試驗(yàn)?zāi)M路徑近乎重合,與實(shí)況路徑之間的偏離也較小,同化的效果并不能明顯表現(xiàn)出來(lái),而此時(shí)刻之后,同化效果開(kāi)始凸現(xiàn),有正效果的(如IR1-AMV)會(huì)使路徑出現(xiàn)回正,減小與實(shí)況的偏差,有負(fù)效果的(如IR3-AMV)則會(huì)加劇偏差的增大。因此綜上所述:不論是控制試驗(yàn)還是兩組對(duì)比試驗(yàn),模擬的臺(tái)風(fēng)路徑都存在偏東的誤差;其中加入了同化過(guò)程的兩組對(duì)比試驗(yàn)出現(xiàn)明顯東偏的時(shí)間要早于控制試驗(yàn);從臺(tái)風(fēng)的路徑和移速來(lái)看,同化IR1-AMV的效果要優(yōu)于同化IR3-AMV的結(jié)果;并且同化效果的凸顯出現(xiàn)在48小時(shí)之后。

    圖2(b)采用同化方案二,每24小時(shí)一次的冷啟動(dòng)將由圖2(a)得出的同化后臺(tái)風(fēng)路徑普遍偏東的現(xiàn)象控制在了一定的范圍內(nèi),使得同化結(jié)果與實(shí)況路徑的偏差明顯減小,甚至二者都優(yōu)于控制試驗(yàn),幾乎與實(shí)況完全擬合,唯一不足之處是每次模式冷啟動(dòng)時(shí)都會(huì)出現(xiàn)離群點(diǎn)(即偏離其他點(diǎn)較遠(yuǎn)的點(diǎn))??梢?jiàn)是否進(jìn)行冷啟動(dòng)對(duì)臺(tái)風(fēng)路徑的預(yù)報(bào)的準(zhǔn)確率有很大影響,且每次冷啟動(dòng)后的積分時(shí)間均較短,因此不能如圖2(a)般很好的反映兩種云導(dǎo)風(fēng)數(shù)據(jù)的優(yōu)劣。

    對(duì)于沒(méi)有經(jīng)過(guò)稀疏處理的導(dǎo)風(fēng)數(shù)據(jù)同化結(jié)果如圖2(c)所示,從圖中可以看出兩次同化的路徑與控制試驗(yàn)的路徑相差無(wú)幾,不同于圖2(a)有明顯的正負(fù)同化效果,從整體上看,圖3同化后的兩條路徑只是在控制試驗(yàn)結(jié)果的基礎(chǔ)上有較小的變化,且路徑迂回蜿蜒,不符合臺(tái)風(fēng)的實(shí)際移動(dòng)特征,如圖中出現(xiàn)的:北移的臺(tái)風(fēng)突然向南折返。出現(xiàn)這種情況,分析其原因主要在于原始資料無(wú)規(guī)律的分散在垂直方向的各個(gè)高度上(經(jīng)統(tǒng)計(jì):原始云導(dǎo)風(fēng)數(shù)據(jù)中IR1通道處于200hPa的觀測(cè)有386個(gè),IR3通道中有562個(gè),而如圖1所示:稀疏后的IR1通道中處于200hPa的觀測(cè)有84760個(gè),IR3通道中有75126個(gè)),因此稀疏前能用于同化的200hPa的數(shù)據(jù)個(gè)數(shù)遠(yuǎn)遠(yuǎn)小于稀疏后的,使得該區(qū)域在該氣壓層只有幾個(gè)格點(diǎn)周?chē)杏^測(cè)值可以同化,使得同化效果不明顯并且影響了整體的平衡性,而造成了路徑回折的現(xiàn)象。

    4.2 物理量場(chǎng)的增量分析

    由于云導(dǎo)風(fēng)是衛(wèi)星反演的大氣運(yùn)動(dòng)矢量,體現(xiàn)的是高空氣流的流動(dòng),因此從同化對(duì)位勢(shì)高度場(chǎng)和風(fēng)場(chǎng)的改善作用進(jìn)行比較分析。由于增量圖的相減作用,去除模式自身的系統(tǒng)偏差和長(zhǎng)時(shí)間積分過(guò)程中模式誤差的積累,可以更好地比較兩種不同通道的導(dǎo)風(fēng)數(shù)據(jù)的優(yōu)劣。

    圖3 對(duì)比試驗(yàn)1,19日00時(shí)200hPa位勢(shì)高度增量圖(實(shí)線為正,虛線為負(fù),單位:gpm)

    圖4 對(duì)比試驗(yàn)2,19日00時(shí)200hPa位勢(shì)高度增量圖(實(shí)線為正,虛線為負(fù),單位:gpm)

    圖5 對(duì)比試驗(yàn)1,20日00時(shí)200hPa位勢(shì)高度增量圖(實(shí)線為正,虛線為負(fù),單位:gpm)

    圖6 對(duì)比試驗(yàn)2,20日00時(shí)200hPa位勢(shì)高度增量圖(實(shí)線為正,虛線為負(fù),單位:gpm)

    圖7 對(duì)比試驗(yàn)1,21日00時(shí)200hPa位勢(shì)高度增量圖(實(shí)線為正,虛線為負(fù),單位:gpm)

    圖8 對(duì)比試驗(yàn)2,21日00時(shí)200hPa位勢(shì)高度增量圖(實(shí)線為正,虛線為負(fù),單位:gpm)

    圖2-8為對(duì)比試驗(yàn)1和對(duì)比試驗(yàn)2的位勢(shì)高度增量圖(同化后減去同化前,各自取3個(gè)時(shí)次)。由圖可知,不論是同化哪個(gè)通道的云導(dǎo)風(fēng)資料,在臺(tái)風(fēng)以外的區(qū)域位勢(shì)高度的變化均在±10gpm左右,并且隨著積分時(shí)間的增長(zhǎng),其變化幅度遠(yuǎn)小于臺(tái)風(fēng)區(qū)域;且都有兩個(gè)強(qiáng)度中心,其中正值中心代表的是同化前的臺(tái)風(fēng)中心,負(fù)值中心代表同化后的臺(tái)風(fēng)中心,且負(fù)值中心均位于正值中心東面,符合同化后臺(tái)風(fēng)中心偏東的結(jié)論。以22日00時(shí)(即積分48小時(shí))為界,對(duì)比圖3、7和圖4、8,可知:在此時(shí)刻之前,兩組試驗(yàn)的正負(fù)中心位置均為西正東負(fù),中心強(qiáng)度為±40gpm;在此時(shí)刻之后的對(duì)比試驗(yàn)1中,正負(fù)中心位置由西東對(duì)稱(chēng)變成了北南對(duì)稱(chēng),中心強(qiáng)度的增幅為30-40gpm,對(duì)比試驗(yàn)2中的正負(fù)中心位置保持西正東負(fù)不變,正負(fù)中心強(qiáng)度卻增加到了±150gpm,增幅達(dá)100gpm以上。因此可以得出,當(dāng)模式連續(xù)積分48小時(shí)后,雖然兩個(gè)實(shí)驗(yàn)都存在誤差的累積增長(zhǎng),但不論是臺(tái)風(fēng)路徑還是中心強(qiáng)度,同化紅外通道導(dǎo)風(fēng)數(shù)據(jù)的效果都要好于同化水汽通道的數(shù)據(jù)。

    在對(duì)臺(tái)風(fēng)的研究中,臺(tái)風(fēng)的移動(dòng)主要取決于周?chē)h(huán)境中引導(dǎo)氣流的影響,而高空氣流的旋度和散度也是評(píng)測(cè)臺(tái)風(fēng)的重要因素,因此風(fēng)場(chǎng)的變化對(duì)臺(tái)風(fēng)預(yù)報(bào)有著較大的意義。雖然導(dǎo)風(fēng)數(shù)據(jù)的高度確定存在誤差,但是對(duì)于觀測(cè)稀少的洋面,云導(dǎo)風(fēng)可以近似看做高空的實(shí)況風(fēng)場(chǎng),因此從風(fēng)場(chǎng)的增量圖中可以清楚的看出同化之后云導(dǎo)風(fēng)對(duì)控制試驗(yàn)中風(fēng)向和風(fēng)速的改變值。

    圖9和圖10所使用的云導(dǎo)風(fēng)均為紅外通道,不同之處在于:前者取自對(duì)比試驗(yàn)1的熱啟動(dòng)同化,后者取自對(duì)比試驗(yàn)3的冷啟動(dòng)同化。如圖9,10所示,20日之前,兩組試驗(yàn)在臺(tái)風(fēng)以北區(qū)域都有較強(qiáng)的東風(fēng)分量的增量;20日以后,在臺(tái)風(fēng)以西區(qū)域出現(xiàn)較強(qiáng)的南風(fēng)氣流增量,北部東風(fēng)增量減弱。

    圖9 對(duì)比試驗(yàn)1的200hPa風(fēng)場(chǎng)增量圖(單位:m/s)

    圖10 對(duì)比試驗(yàn)3(IR1-AMV)的200hPa風(fēng)場(chǎng)增量圖(單位:m/s)

    此外對(duì)比兩組試驗(yàn)臺(tái)風(fēng)區(qū)域的旋度及高空散度可知:熱啟動(dòng)同化在高空會(huì)出現(xiàn)比冷啟動(dòng)同化更強(qiáng)的氣流輻散增量;冷啟動(dòng)同化則會(huì)反映出較熱啟動(dòng)同化更加明顯的反氣旋氣流增量??梢缘贸觯翰徽撌悄姆N同化方法,對(duì)環(huán)境風(fēng)場(chǎng)中的引導(dǎo)氣流都有明顯的增強(qiáng)作用;并且都可以反映出臺(tái)風(fēng)在移動(dòng)過(guò)程中在不斷減弱,但熱啟動(dòng)同化表現(xiàn)在氣流散度的改變,冷啟動(dòng)同化表現(xiàn)在旋度的改變。

    5 結(jié)論

    打破了以往在區(qū)域模式中進(jìn)行同化云導(dǎo)風(fēng)和模擬臺(tái)風(fēng)時(shí),以控制模式連續(xù)積分的時(shí)間(通常為48小時(shí)以?xún)?nèi))來(lái)減小誤差累積的慣常做法,將積分時(shí)間延長(zhǎng)至120小時(shí),并且結(jié)合了兩種同化方案:熱循環(huán)同化和冷啟動(dòng)同化與熱啟動(dòng)同化交替進(jìn)行。以此驗(yàn)證了對(duì)云導(dǎo)風(fēng)數(shù)據(jù)進(jìn)行稀疏處理的必要性,兩種同化方案和兩種通道的導(dǎo)風(fēng)資料的優(yōu)劣。

    對(duì)云導(dǎo)風(fēng)資料進(jìn)行垂直方向上的稀疏化,可以有效地增加其在水平面上的數(shù)據(jù)個(gè)數(shù),確保大部分格點(diǎn)有觀測(cè)可同化,避免同化之后分析場(chǎng)出現(xiàn)明顯的不平衡性(對(duì)比試驗(yàn)4)。雖然連續(xù)積分120小時(shí)使得臺(tái)風(fēng)的模擬路徑較實(shí)況路徑的偏差持續(xù)增大(控制試驗(yàn)),然而同化了紅外通道云導(dǎo)風(fēng)的模擬路徑卻在積分48小時(shí)之后,出現(xiàn)了明顯的正效果(對(duì)比試驗(yàn)1),水汽通道的則為明顯的負(fù)效果(對(duì)比試驗(yàn)2)。而將積分時(shí)間限制在48小時(shí)之內(nèi),雖然能保證臺(tái)風(fēng)模擬路徑最接近于實(shí)況,卻無(wú)法分辨兩個(gè)通道導(dǎo)風(fēng)的質(zhì)量好壞(對(duì)比試驗(yàn)3)。通過(guò)對(duì)位勢(shì)高度增量的分析,不論是對(duì)臺(tái)風(fēng)的路徑還是中心強(qiáng)度的模擬,紅外通道的導(dǎo)風(fēng)資料都要優(yōu)于水汽通道。從風(fēng)場(chǎng)的增量圖中,對(duì)比了熱啟動(dòng)同化和冷啟動(dòng)同化,發(fā)現(xiàn)兩種方法都可以較好的反映出臺(tái)風(fēng)在移動(dòng)過(guò)程中的不斷衰減,但熱啟動(dòng)同化表現(xiàn)在氣流散度的改變,冷啟動(dòng)同化表現(xiàn)在旋度的改變。

    [1]Purdom J F W,Vonder Haar T H,Stewart J K,Leary N E.Diagnosing the severe thunderstorm environment by mesoscale cloud tracking-A new approach and new information[A].IN:Conference on Satellite/Remote Sensing and Applications[C].Clearwater Beach.FL.Preprints.Boston.MA.American Meteorological Society.1984:100-105.

    [2]John F.Le Marshall,Lance M.Leslie,Connie Spinoso.The Generation and Assimilation of Cloud-Drift Winds in Numerical Weather Prediction[J].Meteor.Soc.Japan,1997,75(1):383-393.

    [3]R C Bhatia,P N Khanna,K Prasad.Use of meteosat-5 derived winds for analysis of two tropical cyclones affecting gujarat coast on 20 may,1999 and orissa coast on 29 october,1999[A].IN:5th International Winds Workshop[C].Eur.Organ,2000:135-142.

    [4]Brian J.Soden,Christopher S.Velden,Robert E.The Impact of Satellite Winds on Experimental GFDL Hurricane Model Forecasts[J].Mon Wea Rev,2001,129(4):835-852.

    [5]張守峰,王詩(shī)文.在臺(tái)風(fēng)業(yè)務(wù)系統(tǒng)中使用衛(wèi)星云導(dǎo)風(fēng)資料的試驗(yàn)[J].氣象.1999(8):22-25.

    [6]張守峰,王詩(shī)文.應(yīng)用衛(wèi)星云導(dǎo)風(fēng)進(jìn)行臺(tái)風(fēng)路徑預(yù)報(bào)試驗(yàn)[J].熱帶氣象學(xué)報(bào).1999,15(4):347-355.

    [7]王棟梁,梁旭東,端義宏.云跡風(fēng)在熱帶氣旋路徑數(shù)值預(yù)報(bào)中的應(yīng)用研究[J].氣象學(xué)報(bào),2005,63(3):351-358.

    [8]莊照榮,薛紀(jì)善.云跡風(fēng)資料的三維變分同化及對(duì)臺(tái)風(fēng)預(yù)報(bào)的影響試驗(yàn)[J].熱帶氣象學(xué)報(bào),2004,20(3):225-236.

    [9]何斌,翟國(guó)慶,高坤.云跡風(fēng)資料同化對(duì)東亞海域風(fēng)場(chǎng)數(shù)值模擬的影響[J];海洋學(xué)報(bào),2007,29(6):23-32.

    [10]李華宏,王曼,薛紀(jì)善.FY-2C云跡風(fēng)資料在中尺度數(shù)值模式中的應(yīng)用研究[J].氣象學(xué)報(bào).2008,66(1):50-58.

    [11]黃彥彬,雷恒池,王振會(huì).衛(wèi)星云跡風(fēng)資料對(duì)中尺度數(shù)值模式初始風(fēng)場(chǎng)改進(jìn)試驗(yàn)[J].南京氣象學(xué)院學(xué)報(bào),2003,26(5):668-676.

    [12]馮文,萬(wàn)齊林,陳子通.逐時(shí)云跡風(fēng)資料同化對(duì)暴雨預(yù)報(bào)的模擬試驗(yàn)[J].氣象學(xué)報(bào),2008,66(4):500-512.

    [13]周兵,徐海明,吳國(guó)雄.云跡風(fēng)資料同化對(duì)暴雨預(yù)報(bào)影響的數(shù)值模擬[J].氣象學(xué)報(bào),2002,60(3):309-317.

    [14]何志新,沈桐立,丁治英.云跡風(fēng)資料在臺(tái)風(fēng)數(shù)值模擬中的應(yīng)用[D].南京,南京信息工程大學(xué),2009.

    [15]許健民.風(fēng)云二號(hào)衛(wèi)星業(yè)務(wù)產(chǎn)品與衛(wèi)星數(shù)據(jù)格式實(shí)用手冊(cè)[M].北京:氣象出版社,2008:17-23.

    猜你喜歡
    冷啟動(dòng)鯰魚(yú)實(shí)況
    輕型汽油車(chē)實(shí)際行駛排放試驗(yàn)中冷啟動(dòng)排放的評(píng)估
    基于學(xué)習(xí)興趣的冷啟動(dòng)推薦模型
    客聯(lián)(2021年2期)2021-09-10 07:22:44
    鄉(xiāng)村小學(xué)的愿望與現(xiàn)實(shí)——宜君鄉(xiāng)村教育實(shí)況
    天舟一號(hào)貨運(yùn)飛船發(fā)射實(shí)況掠影
    可愛(ài)潮咖們的獨(dú)門(mén)彩妝實(shí)況直播
    讓鯰魚(yú)慌起來(lái)
    Space X會(huì)是攪局“鯰魚(yú)”?
    太空探索(2014年7期)2014-07-12 15:16:47
    陷 阱
    軍事技能“冷啟動(dòng)”式訓(xùn)練理念初探
    陷阱
    知識(shí)窗(2007年6期)2007-05-14 09:07:34
    日韩av不卡免费在线播放| 国产成人免费观看mmmm| 男女视频在线观看网站免费| 国产69精品久久久久777片| 真实男女啪啪啪动态图| 国产精品久久久久久久电影| 七月丁香在线播放| 2021少妇久久久久久久久久久| 一级av片app| 国产伦精品一区二区三区四那| 国产高清不卡午夜福利| 色吧在线观看| 大香蕉97超碰在线| 日韩电影二区| 午夜福利高清视频| 91久久精品国产一区二区三区| 男插女下体视频免费在线播放| 99久国产av精品| 亚洲无线观看免费| 五月伊人婷婷丁香| 男女视频在线观看网站免费| videos熟女内射| 黄色日韩在线| 精品一区二区免费观看| 18禁裸乳无遮挡免费网站照片| 丝瓜视频免费看黄片| 久久久久国产网址| 亚洲欧洲国产日韩| 精品久久久精品久久久| 夜夜爽夜夜爽视频| 91午夜精品亚洲一区二区三区| 高清av免费在线| 亚洲欧美清纯卡通| 亚洲三级黄色毛片| 中文天堂在线官网| 人人妻人人澡欧美一区二区| 国产综合精华液| 成年免费大片在线观看| www.av在线官网国产| 国产真实伦视频高清在线观看| 舔av片在线| 97精品久久久久久久久久精品| 免费人成在线观看视频色| 麻豆久久精品国产亚洲av| 色综合亚洲欧美另类图片| 久久韩国三级中文字幕| 中文在线观看免费www的网站| 国产白丝娇喘喷水9色精品| 精品久久久久久久久亚洲| av一本久久久久| 亚洲人与动物交配视频| 久久久欧美国产精品| 日韩人妻高清精品专区| 国产精品女同一区二区软件| 99热全是精品| 伊人久久精品亚洲午夜| 国产一区二区在线观看日韩| 亚洲国产最新在线播放| 黄色欧美视频在线观看| 2022亚洲国产成人精品| 久久综合国产亚洲精品| 欧美极品一区二区三区四区| 亚洲第一区二区三区不卡| 黑人高潮一二区| 高清欧美精品videossex| 男女那种视频在线观看| 国内少妇人妻偷人精品xxx网站| 又粗又硬又长又爽又黄的视频| 国产亚洲精品久久久com| 真实男女啪啪啪动态图| 亚洲av中文av极速乱| 狂野欧美激情性xxxx在线观看| 三级毛片av免费| 国产在线男女| 中文欧美无线码| 久久99热这里只频精品6学生| 国产一级毛片在线| 大陆偷拍与自拍| 精品一区二区免费观看| 成人性生交大片免费视频hd| 91狼人影院| 日韩欧美精品v在线| 亚洲av一区综合| 婷婷色av中文字幕| 欧美日韩综合久久久久久| 国产精品一及| 乱系列少妇在线播放| 欧美变态另类bdsm刘玥| 久久精品久久久久久噜噜老黄| 日韩av在线大香蕉| 免费少妇av软件| 午夜福利高清视频| 女人十人毛片免费观看3o分钟| 久久久久久久国产电影| 欧美xxⅹ黑人| 青春草亚洲视频在线观看| 午夜激情欧美在线| 中国美白少妇内射xxxbb| 看十八女毛片水多多多| 高清日韩中文字幕在线| 搡老妇女老女人老熟妇| 晚上一个人看的免费电影| 少妇高潮的动态图| 亚州av有码| 好男人在线观看高清免费视频| 国产淫语在线视频| 日本色播在线视频| 午夜精品在线福利| 欧美+日韩+精品| eeuss影院久久| 亚洲精品影视一区二区三区av| 又爽又黄a免费视频| 免费观看av网站的网址| 99九九线精品视频在线观看视频| 高清毛片免费看| 免费观看在线日韩| 久久综合国产亚洲精品| 99久久人妻综合| 日本免费在线观看一区| 最近视频中文字幕2019在线8| 超碰97精品在线观看| 国产毛片a区久久久久| 免费av毛片视频| 亚洲四区av| 小蜜桃在线观看免费完整版高清| 久久这里只有精品中国| 国产精品嫩草影院av在线观看| 秋霞伦理黄片| 黄片wwwwww| videossex国产| 免费少妇av软件| 亚洲av国产av综合av卡| videossex国产| 夜夜爽夜夜爽视频| 欧美日韩精品成人综合77777| 日韩欧美精品免费久久| 欧美潮喷喷水| 国产一区二区在线观看日韩| 国产老妇伦熟女老妇高清| 看黄色毛片网站| 亚洲欧美精品专区久久| 国产精品精品国产色婷婷| 国产精品国产三级专区第一集| 免费看a级黄色片| 色综合亚洲欧美另类图片| 国产精品女同一区二区软件| 麻豆成人午夜福利视频| 中文字幕免费在线视频6| 成人一区二区视频在线观看| 白带黄色成豆腐渣| 国产精品久久久久久久电影| 亚洲熟女精品中文字幕| 天堂√8在线中文| 欧美激情在线99| 乱系列少妇在线播放| 国产大屁股一区二区在线视频| 欧美+日韩+精品| 在线观看av片永久免费下载| 色综合亚洲欧美另类图片| 偷拍熟女少妇极品色| 亚洲欧美中文字幕日韩二区| 亚洲久久久久久中文字幕| 国产成人一区二区在线| 午夜日本视频在线| 丰满乱子伦码专区| 免费观看无遮挡的男女| 日本午夜av视频| 九色成人免费人妻av| 女人十人毛片免费观看3o分钟| av女优亚洲男人天堂| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 亚洲av福利一区| 六月丁香七月| 日本三级黄在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲欧美一区二区三区黑人 | 久久国内精品自在自线图片| 国产亚洲一区二区精品| 国产探花在线观看一区二区| 久久久久国产网址| 国产亚洲5aaaaa淫片| 一级片'在线观看视频| av卡一久久| 亚洲一区高清亚洲精品| 99久久精品一区二区三区| 久久久久久久久久人人人人人人| 免费观看a级毛片全部| 激情 狠狠 欧美| 男人爽女人下面视频在线观看| 最近最新中文字幕大全电影3| 三级经典国产精品| 久久精品久久久久久噜噜老黄| 国产成人a∨麻豆精品| 一级av片app| 欧美激情在线99| 18禁裸乳无遮挡免费网站照片| 国产精品美女特级片免费视频播放器| 成人美女网站在线观看视频| 少妇高潮的动态图| 人人妻人人澡欧美一区二区| 日韩av免费高清视频| 国语对白做爰xxxⅹ性视频网站| 18禁动态无遮挡网站| 国产成人午夜福利电影在线观看| 天堂中文最新版在线下载 | 中国美白少妇内射xxxbb| ponron亚洲| 最近最新中文字幕免费大全7| 亚洲精品影视一区二区三区av| 亚洲美女搞黄在线观看| 国内精品一区二区在线观看| 精品久久久久久久久亚洲| 老司机影院成人| 日韩国内少妇激情av| 爱豆传媒免费全集在线观看| 精品久久久久久成人av| 日韩精品青青久久久久久| 欧美xxxx黑人xx丫x性爽| 国产精品人妻久久久久久| www.色视频.com| 简卡轻食公司| 精品熟女少妇av免费看| 欧美成人a在线观看| 久久久欧美国产精品| 日日撸夜夜添| 久久久久久久国产电影| 亚洲自拍偷在线| 国产精品嫩草影院av在线观看| 黄色配什么色好看| 亚洲av电影在线观看一区二区三区 | 国产乱来视频区| 在线a可以看的网站| 搡老乐熟女国产| 中文字幕人妻熟人妻熟丝袜美| 久久热精品热| 禁无遮挡网站| 亚洲经典国产精华液单| 一区二区三区四区激情视频| 亚洲精华国产精华液的使用体验| 亚洲精品乱久久久久久| 亚洲精品一区蜜桃| 国产爱豆传媒在线观看| 亚洲av男天堂| 亚洲成人精品中文字幕电影| 色综合亚洲欧美另类图片| 好男人视频免费观看在线| 国产淫语在线视频| 国产 一区精品| 性插视频无遮挡在线免费观看| 免费人成在线观看视频色| 亚洲国产精品sss在线观看| 国产精品久久久久久av不卡| 国产日韩欧美在线精品| 国产午夜精品久久久久久一区二区三区| 毛片一级片免费看久久久久| 国产精品麻豆人妻色哟哟久久 | 国产精品国产三级国产av玫瑰| 免费观看的影片在线观看| 精品亚洲乱码少妇综合久久| 欧美日韩综合久久久久久| 又粗又硬又长又爽又黄的视频| 免费在线观看成人毛片| 亚洲精品久久午夜乱码| 久久久久九九精品影院| 亚洲av成人精品一区久久| 精品酒店卫生间| 男人舔女人下体高潮全视频| 嫩草影院入口| 99久国产av精品| 精品少妇黑人巨大在线播放| 亚洲精品,欧美精品| 国产综合精华液| 午夜精品一区二区三区免费看| 欧美日韩在线观看h| 一二三四中文在线观看免费高清| 亚洲天堂国产精品一区在线| 1000部很黄的大片| 听说在线观看完整版免费高清| 超碰97精品在线观看| 国产高清不卡午夜福利| 天堂√8在线中文| 青春草国产在线视频| 亚洲国产日韩欧美精品在线观看| 国产黄色免费在线视频| 淫秽高清视频在线观看| 简卡轻食公司| 97超碰精品成人国产| 免费高清在线观看视频在线观看| 国产视频首页在线观看| 51国产日韩欧美| 18+在线观看网站| 日本免费在线观看一区| 亚洲av成人av| 成人av在线播放网站| 色5月婷婷丁香| 日本与韩国留学比较| 一夜夜www| 狂野欧美激情性xxxx在线观看| 一级av片app| 日本熟妇午夜| 免费观看a级毛片全部| 91精品伊人久久大香线蕉| 十八禁国产超污无遮挡网站| 22中文网久久字幕| 午夜福利在线观看吧| 内射极品少妇av片p| 国产黄a三级三级三级人| 国产黄色视频一区二区在线观看| 纵有疾风起免费观看全集完整版 | 国产亚洲av片在线观看秒播厂 | 久久99热6这里只有精品| 91久久精品电影网| 精品一区二区三区视频在线| 国产女主播在线喷水免费视频网站 | 亚洲人成网站在线观看播放| 一边亲一边摸免费视频| 免费观看精品视频网站| 天堂√8在线中文| 国产色爽女视频免费观看| 99久久精品热视频| 别揉我奶头 嗯啊视频| 欧美bdsm另类| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 别揉我奶头 嗯啊视频| 日韩伦理黄色片| 蜜桃亚洲精品一区二区三区| 亚洲熟女精品中文字幕| 欧美成人精品欧美一级黄| 只有这里有精品99| 国产精品美女特级片免费视频播放器| 91午夜精品亚洲一区二区三区| 亚洲国产色片| 国产乱人视频| 少妇熟女欧美另类| 国产精品福利在线免费观看| 免费看日本二区| 亚洲成人精品中文字幕电影| 99久国产av精品| 国产真实伦视频高清在线观看| 久久99热6这里只有精品| 亚洲精品乱码久久久久久按摩| 在线免费十八禁| 国产一级毛片七仙女欲春2| 精品久久久精品久久久| 免费观看在线日韩| 精品一区二区三区视频在线| 在线观看av片永久免费下载| 久久人人爽人人片av| 22中文网久久字幕| 少妇熟女欧美另类| 欧美日韩在线观看h| 久久精品久久久久久久性| 午夜福利在线在线| 女人被狂操c到高潮| 国产三级在线视频| 久热久热在线精品观看| 免费播放大片免费观看视频在线观看| 成年av动漫网址| 精品国产三级普通话版| av线在线观看网站| 在线免费观看的www视频| 免费黄频网站在线观看国产| 午夜激情久久久久久久| 可以在线观看毛片的网站| 亚洲国产欧美在线一区| 国产免费又黄又爽又色| 精品久久久久久久久av| 97超碰精品成人国产| 嫩草影院精品99| 中文资源天堂在线| 视频中文字幕在线观看| 色综合色国产| 午夜老司机福利剧场| 免费av不卡在线播放| 午夜福利在线观看吧| 欧美日韩综合久久久久久| 欧美日韩国产mv在线观看视频 | 日本熟妇午夜| 嫩草影院入口| 亚洲第一区二区三区不卡| 午夜精品在线福利| 欧美精品一区二区大全| 97精品久久久久久久久久精品| 亚洲av中文av极速乱| 麻豆乱淫一区二区| 国产亚洲5aaaaa淫片| 成年免费大片在线观看| 精品久久久久久久久亚洲| 中文字幕免费在线视频6| 亚洲美女视频黄频| 男女那种视频在线观看| 国产av不卡久久| 精品国内亚洲2022精品成人| 国产亚洲最大av| 久久亚洲国产成人精品v| 国产综合精华液| 人妻制服诱惑在线中文字幕| 国模一区二区三区四区视频| 亚洲欧美一区二区三区黑人 | 五月伊人婷婷丁香| 伊人久久国产一区二区| 中文精品一卡2卡3卡4更新| 精品一区二区三卡| 国产高清国产精品国产三级 | 欧美97在线视频| 亚洲欧洲日产国产| 国产亚洲精品av在线| 超碰av人人做人人爽久久| 精品少妇黑人巨大在线播放| 六月丁香七月| 身体一侧抽搐| 国产乱来视频区| 日本色播在线视频| av黄色大香蕉| 一级片'在线观看视频| 亚洲av中文字字幕乱码综合| 日本免费a在线| 国产毛片a区久久久久| 午夜福利在线观看吧| av.在线天堂| 国产片特级美女逼逼视频| 在线观看av片永久免费下载| 又爽又黄无遮挡网站| 久久久久精品久久久久真实原创| 国产高清三级在线| 91久久精品国产一区二区三区| 伦精品一区二区三区| 久久99热这里只频精品6学生| 看十八女毛片水多多多| 一本一本综合久久| 日韩av在线免费看完整版不卡| 国产视频内射| 亚洲欧洲国产日韩| 国产免费又黄又爽又色| 亚洲av电影在线观看一区二区三区 | 老司机影院毛片| 亚洲四区av| 91午夜精品亚洲一区二区三区| 美女高潮的动态| 26uuu在线亚洲综合色| 又爽又黄a免费视频| 一级片'在线观看视频| 亚洲精品一二三| 69人妻影院| 老司机影院毛片| 在线 av 中文字幕| 街头女战士在线观看网站| 亚洲国产成人一精品久久久| 爱豆传媒免费全集在线观看| 国产又色又爽无遮挡免| 一区二区三区四区激情视频| 亚洲国产高清在线一区二区三| 乱系列少妇在线播放| 黄片无遮挡物在线观看| 国产乱人偷精品视频| 国产男人的电影天堂91| 91精品一卡2卡3卡4卡| 国精品久久久久久国模美| 夫妻性生交免费视频一级片| 在现免费观看毛片| 国产精品国产三级专区第一集| 久久久久久久大尺度免费视频| 国产 亚洲一区二区三区 | 国产日韩欧美在线精品| 听说在线观看完整版免费高清| 99久久九九国产精品国产免费| 免费观看的影片在线观看| 日韩精品青青久久久久久| 女人十人毛片免费观看3o分钟| 国产精品一区二区三区四区免费观看| 亚洲熟女精品中文字幕| 乱码一卡2卡4卡精品| 夫妻午夜视频| 国产毛片a区久久久久| 99久久精品国产国产毛片| 亚洲精品自拍成人| 精品99又大又爽又粗少妇毛片| 黑人高潮一二区| 97超碰精品成人国产| 日韩成人av中文字幕在线观看| 亚洲av二区三区四区| 国产片特级美女逼逼视频| 街头女战士在线观看网站| 直男gayav资源| 又粗又硬又长又爽又黄的视频| 国产成人一区二区在线| 亚洲va在线va天堂va国产| 男女啪啪激烈高潮av片| 99热这里只有是精品在线观看| 91aial.com中文字幕在线观看| 久久精品国产亚洲av天美| 国产亚洲91精品色在线| 日韩国内少妇激情av| 美女cb高潮喷水在线观看| 日韩av在线免费看完整版不卡| 亚洲精品国产av蜜桃| 国产男人的电影天堂91| 十八禁网站网址无遮挡 | 国产亚洲5aaaaa淫片| 国产久久久一区二区三区| 少妇的逼好多水| 日本熟妇午夜| 日韩亚洲欧美综合| 日本色播在线视频| 夜夜爽夜夜爽视频| 最近中文字幕高清免费大全6| 免费观看在线日韩| 超碰97精品在线观看| 丰满人妻一区二区三区视频av| 亚洲人成网站在线观看播放| 美女内射精品一级片tv| 在线观看av片永久免费下载| 内地一区二区视频在线| 国产亚洲午夜精品一区二区久久 | 九九久久精品国产亚洲av麻豆| 日韩视频在线欧美| av在线播放精品| 欧美xxⅹ黑人| 国产精品久久久久久久久免| 国产精品人妻久久久影院| 中文资源天堂在线| 日本一二三区视频观看| 2021天堂中文幕一二区在线观| 中文在线观看免费www的网站| or卡值多少钱| 人体艺术视频欧美日本| 少妇熟女欧美另类| 国产精品一及| 国产成人一区二区在线| 国产成人a∨麻豆精品| 国产av码专区亚洲av| 亚洲四区av| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 男女边吃奶边做爰视频| 中文资源天堂在线| 欧美zozozo另类| 日韩欧美国产在线观看| 九色成人免费人妻av| 91精品国产九色| 精华霜和精华液先用哪个| 亚洲不卡免费看| 免费不卡的大黄色大毛片视频在线观看 | 大又大粗又爽又黄少妇毛片口| 国产国拍精品亚洲av在线观看| 国产在线男女| 校园人妻丝袜中文字幕| 国产毛片a区久久久久| 亚洲精品影视一区二区三区av| 国产白丝娇喘喷水9色精品| 色播亚洲综合网| 色网站视频免费| 国产精品嫩草影院av在线观看| 亚洲激情五月婷婷啪啪| 日本黄色片子视频| 欧美日韩在线观看h| 精品久久久久久成人av| 国产 亚洲一区二区三区 | 国产在线一区二区三区精| 日日撸夜夜添| 成人美女网站在线观看视频| 天堂√8在线中文| 久久这里有精品视频免费| 99九九线精品视频在线观看视频| 欧美日本视频| 床上黄色一级片| 成人美女网站在线观看视频| 夫妻午夜视频| 天堂av国产一区二区熟女人妻| 最近最新中文字幕免费大全7| 在线观看人妻少妇| 亚洲国产欧美人成| 久久久久久久久久成人| 一区二区三区四区激情视频| 亚洲成人av在线免费| 国产亚洲精品久久久com| 亚洲精品亚洲一区二区| 一级av片app| 亚洲av.av天堂| 国产欧美日韩精品一区二区| 国产精品女同一区二区软件| 久久久久久久久久久免费av| 亚洲四区av| 免费高清在线观看视频在线观看| 亚洲av电影不卡..在线观看| 97超视频在线观看视频| 国产午夜精品一二区理论片| 91在线精品国自产拍蜜月| 三级男女做爰猛烈吃奶摸视频| 人妻少妇偷人精品九色| 精品一区二区免费观看| av在线亚洲专区| 亚洲国产色片| 久久久久久伊人网av| 免费看光身美女| 麻豆成人午夜福利视频| 国产成人91sexporn| 国产成人一区二区在线| 简卡轻食公司| 国产成人免费观看mmmm| 日韩,欧美,国产一区二区三区| 乱系列少妇在线播放| 99久久中文字幕三级久久日本| 久久鲁丝午夜福利片| 秋霞在线观看毛片| 爱豆传媒免费全集在线观看| 欧美不卡视频在线免费观看| 久久久久性生活片| 天堂√8在线中文| 精品国产一区二区三区久久久樱花 | 亚洲av成人精品一二三区| 久久久久免费精品人妻一区二区| av女优亚洲男人天堂|