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

    往返平飄式探空觀測系統(tǒng)對CMA-MESO的影響研究*

    2023-03-17 05:19:42王瑞文王金成陶亦為田偉紅
    氣象 2023年1期
    關(guān)鍵詞:探空降水數(shù)值

    王瑞文 王金成 王 丹 陶亦為 田偉紅

    1中國氣象局地球系統(tǒng)數(shù)值預(yù)報(bào)中心,北京 100081 2 災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室,北京 100081 3 國家氣象中心,北京 100081

    提 要: 往返平飄式探空觀測是我國研發(fā)的一種新型高空觀測技術(shù),除了具備與傳統(tǒng)探空觀測一致的上升段大氣垂直廓線觀測能力,同時(shí)還增加了平飄段和下降段的大氣探測,自動實(shí)現(xiàn)了探測廓線的時(shí)空加密。利用ERA5再分析資料作為“真值”,利用往返平飄式探空模擬仿真系統(tǒng)構(gòu)造了往返式探空模擬觀測,基于CMA-MESO區(qū)域模式和3D-Var同化系統(tǒng)進(jìn)行了觀測系統(tǒng)模擬試驗(yàn)(Observing System Simulation Experiments,OSSEs)。數(shù)值試驗(yàn)結(jié)果表明:相比傳統(tǒng)單次上升段探空觀測,往返平飄式探空在全國組網(wǎng)的情況下,其增加的下降段模擬探空觀測,能夠有效提高CMA-MESO的降水預(yù)報(bào)技巧,不同降水量級的ETS評分提高約2%~5%,同時(shí)改進(jìn)要素場(溫、濕場和風(fēng)場)的預(yù)報(bào),改進(jìn)率約為2%~5%。此外,典型天氣個(gè)例分析結(jié)果表明,增加往返平飄式探空觀測能夠改善模式初值偏差,從而更準(zhǔn)確地模擬降水分布。該文的研究結(jié)論為往返平飄式探空的未來科學(xué)布局和應(yīng)用提供了理論支撐。

    引 言

    數(shù)值天氣預(yù)報(bào)是基于數(shù)學(xué)物理方法客觀定量計(jì)算未來天氣演變的科學(xué),正經(jīng)歷了一場“靜悄悄的革命”(Bauer et al,2015),在天氣預(yù)報(bào)中處于決定性位置。數(shù)值預(yù)報(bào)實(shí)際上是一個(gè)初值問題(Bjerknes,1904),初值的準(zhǔn)確程度高度影響預(yù)報(bào)結(jié)果。資料同化是為數(shù)值天氣模式提供初值的理論技術(shù),而觀測資料是對大氣真實(shí)狀態(tài)的測量,是資料同化重要信息源。作為綜合氣象觀測系統(tǒng)的重要組成部分,探空觀測能夠提供完整的三維大氣的溫度、氣壓、濕度、風(fēng)等信息,是高空觀測體系中最為可信的基準(zhǔn)資料(Faccani et al, 2009),被廣泛應(yīng)用于天氣預(yù)報(bào)(葉篤正,1977a;1977b;王笑芳和丁一匯,1994;劉曉璐等,2014;劉超等,2017)、數(shù)值天氣預(yù)報(bào)領(lǐng)域(Radnóti et al,2012;郝民等,2014;Privé et al,2014;Ingleby et al,2016)。其質(zhì)量穩(wěn)定可靠、垂直層數(shù)多,常常作為同化分析中的基礎(chǔ)資料(于曉晶等,2018),也是衛(wèi)星輻射率和飛機(jī)報(bào)溫度等資料進(jìn)行偏差訂正的錨定資料(Eyre,2016; Cucurull et al,2014;王瑞文等,2017;吳瓊等,2020),在模式場檢驗(yàn)、分析中也具有重要作用。另有研究表明,探空資料時(shí)間觀測頻次加倍,將顯著提高數(shù)值天氣預(yù)報(bào)的預(yù)報(bào)技巧(Privé et al,2014;王丹等,2019)。但是由于受觀測成本的限制,目前我國只在00時(shí)和12時(shí)(世界時(shí))開展觀測,然而,每天兩次的觀測頻次并不能滿足中小尺度災(zāi)害性天氣的預(yù)測和預(yù)警要求(陳敏等,2011)。2011年起,我國探空系統(tǒng)也全面進(jìn)行更新?lián)Q代,由59型探空系統(tǒng)逐步升級為L波段探空系統(tǒng),它的采樣時(shí)間間隔為1.2 s,垂直間隔約為8 m,能提供更為密集的大氣垂直探測信息(郝民等,2014)。高分辨率的秒級探空資料可以更好地反映大氣層結(jié)的細(xì)致特征,有助于預(yù)報(bào)員進(jìn)行更為精細(xì)的天氣學(xué)分析(Allen and Vincent,1995;Wang and Geller,2003;王麗吉和楊程,2018),提高天氣預(yù)報(bào)的定性分析結(jié)果。在定量應(yīng)用研究方面,為了滿足數(shù)值預(yù)報(bào)的精細(xì)化需求,中國氣象局氣象探測中心近期提出了創(chuàng)新型的往返平飄式探空觀測技術(shù),即通過一次探空施放,同時(shí)實(shí)現(xiàn)“上升段-平飄段-下降段”三段觀測(曹曉鐘等,2019;錢媛,2019;王丹等,2020)。這個(gè)觀測手段不僅填補(bǔ)了平流層的資料空缺,同時(shí)通過平飄完成后的下降爆破階段進(jìn)行大氣狀態(tài)的再次測量,自動實(shí)現(xiàn)了大氣垂直探測的頻次加密。這種新型觀測技術(shù)還正處于研發(fā)初期,目前實(shí)現(xiàn)了小范圍組網(wǎng)布局(江淮區(qū)域),在大規(guī)模業(yè)務(wù)布網(wǎng)前,給出定性和定量的布局影響評估十分必要。該觀測技術(shù)的資料對于天氣預(yù)報(bào)尤其是數(shù)值天氣預(yù)報(bào)的影響如何?是否能被數(shù)值模式的同化系統(tǒng)更好應(yīng)用?這些科學(xué)問題都有待評估和驗(yàn)證。而觀測系統(tǒng)模擬試驗(yàn)(Observing System Simulation Experiments,OSSEs)是全面評估新型觀測資料對數(shù)值天氣預(yù)報(bào)潛在影響的一種重要手段。

    OSSEs方法是最先被用在氣象領(lǐng)域來評估前瞻性觀測(當(dāng)前不可大量獲得)對改進(jìn)數(shù)值預(yù)報(bào)潛在的影響(Lord et al,1997;Atlas et al,2003),同時(shí)還可以用來指導(dǎo)觀測站網(wǎng)的設(shè)計(jì)布局,測試新型觀測系統(tǒng)(Riishojgaard et al,2012),從而優(yōu)化對資料同化和預(yù)報(bào)系統(tǒng)的使用,改善天氣預(yù)報(bào),最終產(chǎn)生最大的社會和經(jīng)濟(jì)效益(Arnold and Dey,1986)。國內(nèi)在探空觀測的OSSEs方面也有一些嘗試,張宇等(2016)認(rèn)為通過增加同化青藏高原上構(gòu)造的探空偽觀測,能夠減少高原地區(qū)的初值不確定性,從而減小對下游預(yù)報(bào)的影響;王丹等(2019)就14時(shí)(北京時(shí))探空是否需要加密的觀測系統(tǒng)模擬試驗(yàn)對探空頻次加密的潛在意義進(jìn)行討論,結(jié)果表明14時(shí)(北京時(shí))加密探空對數(shù)值預(yù)報(bào)具有顯著改進(jìn)作用。由于往返平飄式探空是一種新的觀測技術(shù),以后大規(guī)模的業(yè)務(wù)化應(yīng)用是否可行,對中國數(shù)值預(yù)報(bào)的影響如何,這急需一個(gè)客觀定量的評估。

    本文利用ERA5再分析資料和往返平飄式軌跡預(yù)測系統(tǒng)(王金成等, 2021),構(gòu)造往返平飄式探空的模擬觀測,基于CMA-MESO,開展了往返平飄式探空的觀測系統(tǒng)模擬試驗(yàn),研究結(jié)果揭示了新型觀測資料對數(shù)值預(yù)報(bào)的影響,對往返平飄式探空的未來組網(wǎng)、布局具有一定指導(dǎo)意義。

    1 往返平飄式探空模擬觀測資料構(gòu)造

    采用OSSEs來分析評估往返平飄式探空在我國組網(wǎng)觀測情境下對數(shù)值天氣預(yù)報(bào)的影響。OSSEs通常包含3個(gè)部分: 替代“真實(shí)”大氣的“自然”場;由“自然”場生成的人造觀測,也叫偽觀測;用于評價(jià)“人造觀測”的資料同化和數(shù)值預(yù)報(bào)系統(tǒng)。由于ERA5再分析資料比較接近真實(shí)大氣狀態(tài),本文用ERA5再分析資料作為“自然”場,試驗(yàn)時(shí)間段是2019年6月1日06時(shí)至29日18時(shí)(世界時(shí),下同)。模擬觀測的構(gòu)造過程是:(1)將ERA5再分析資料插值到往返平飄式探空下降階段的經(jīng)緯度位置;(2)將該處構(gòu)造的觀測加上隨機(jī)誤差,隨機(jī)誤差是由探空資料該層的觀測誤差和由蒙特卡羅方法生成的隨機(jī)數(shù)相乘得到;(3)在OSSEs中,同化觀測資料的數(shù)值預(yù)報(bào)系統(tǒng)要和生成“自然”場的模式系統(tǒng)不同,這是為了避免“雙生”效應(yīng)(identical twin)造成資料影響被過高估計(jì)(Halem and Dlouhy,1984),因此,本文選取中國氣象局自主研發(fā)的區(qū)域高分辨率數(shù)值天氣預(yù)報(bào)系統(tǒng)CMA-MESO V5.0作為評價(jià)系統(tǒng)。

    由于往返平飄式探空從傳統(tǒng)探空站釋放觀測氣球,平飄一段時(shí)間(約4 h)后,其下降階段的觀測在空間位置上發(fā)生了變化,因此和傳統(tǒng)的探空相比,在垂向觀測上相當(dāng)于空間加密(和原觀測站位置較遠(yuǎn))和時(shí)間加密(00時(shí)釋放氣球后5 h開始下降,相當(dāng)于06時(shí))。下降階段的觀測位置在下降軌跡上,由軌跡預(yù)測系統(tǒng)得到(王金成等,2021),圖1是模擬的往返平飄式探空的三維軌跡。從圖中可以發(fā)現(xiàn):平飄探空明顯分為上升、平飄和下降三個(gè)階段,形成三個(gè)不同位置的觀測,和傳統(tǒng)探空相比,增加了平飄和下降階段的觀測,觀測空間得到加密。圖2a、2b分別是2019年6—7月和11—12月軌跡預(yù)測系統(tǒng)模擬的往返平飄式探空下降點(diǎn)位置平均分布。在夏季,探空氣球在上升和下降階段會受到一定的西風(fēng)帶影響,向東飄移,但是當(dāng)氣球繼續(xù)上升到平流層,主要受到風(fēng)速較大的平流層?xùn)|風(fēng)帶的影響,經(jīng)過平流層和對流層的不同風(fēng)帶的綜合影響,最后探空氣球的降落位置會在釋放的傳統(tǒng)探空氣球位置偏西。而在冬季,平流層有從西向東的西風(fēng)帶,最后模擬平飄探空降落的位置在傳統(tǒng)探空站以東。從圖2 可知,下降段探空起到了全網(wǎng)站點(diǎn)自動加密的作用。圖3給出了2019年6—7月構(gòu)造的模擬平飄探空觀測上升段和實(shí)際傳統(tǒng)探空的溫度、相對濕度、U風(fēng)和V風(fēng)的偏差和標(biāo)準(zhǔn)差,從對比結(jié)果可知,構(gòu)造的溫度模擬觀測和實(shí)際探空觀測的偏差在0.5 K 以內(nèi),標(biāo)準(zhǔn)差在2 K以內(nèi),相對濕度由于實(shí)際的觀測誤差比較大,因而構(gòu)造的模擬觀測和實(shí)際探空的相對濕度的偏差最大約為20%,而標(biāo)準(zhǔn)差約為12%,U風(fēng)和V風(fēng)的偏差在1 m·s-1以內(nèi),標(biāo)準(zhǔn)差在2~3 m·s-1。本文在傳統(tǒng)探空觀測資料的基礎(chǔ)上,研究觀測在時(shí)空加密的情況下對數(shù)值預(yù)報(bào)的影響。

    圖1 模擬平飄探空三維軌跡圖(灰白色代表軌跡,彩色陰影是地形高度)Fig.1 A three-dimensional trajectory map of simulated return sounding (gray and white: trajectory, colored shadow: height of terrain)

    圖2 2019年(a)夏季(6—7月)和(b)冬季(11—12月)00時(shí)傳統(tǒng)探空(藍(lán)色)和06時(shí)往返平飄式探空下降段平均位置(粉色)分布Fig.2 Distributions of classical radiosonde location (blue) at 00 UTC and return radiosonde descending location (pink) at 06 UTC in (a) summer (June-July) and (b) winter (November-December) of 2019

    圖3 2019年6—7月模擬觀測和實(shí)際傳統(tǒng)探空觀測的偏差(Bias)和標(biāo)準(zhǔn)差(Std)(a)溫度,(b)相對濕度,(c)U風(fēng),(d)V風(fēng)Fig.3 Bias and standard deviations (Std) of the simulated and actual classic sounding observations from June to July of 2019(a) temperature, (b) relative humidity, (c) U wind, (d) V wind

    2 數(shù)值試驗(yàn)

    2.1 方案設(shè)計(jì)

    我國自主研發(fā)的GRAPES是具有完全知識產(chǎn)權(quán)的數(shù)值天氣預(yù)報(bào)系統(tǒng)(薛紀(jì)善和陳德輝,2008),本文采用的GRAPES區(qū)域高分辨率(CMA-MESO V5.0)于2020年6月實(shí)現(xiàn)了業(yè)務(wù)化。CMA-MESO V5.0水平分辨率是0.03°×0.03°,模式覆蓋區(qū)域范圍是10°~60°N、70°~145°E,垂直層次為49層,模式層頂高為35 000 m。積分時(shí)間步長為30 s,采用3 h 循環(huán)同化,每12 h冷啟動一次,冷啟動背景場是由NCEP GFS (https:∥www.nco.ncep.noaa.gov/pmb/products/gfs/)的6 h預(yù)報(bào)場提供。

    本文采用CMA-MESO V5.0的三維變分(3D-Var)同化系統(tǒng),同化的資料有常規(guī)探空、地面報(bào)、船舶、飛機(jī)報(bào)等常規(guī)資料,還有GPS/PW(global position system/precipitable water)大氣可降水量、雷達(dá)VAD(velocity azimuth display)風(fēng)和GNSS/RO(global navigation satellite system/radio occultation)折射率等觀測信息。同化時(shí)間窗口為3 h,每天00時(shí)和12時(shí)冷啟動,每3 h更新循環(huán)同化,由于新增的模擬觀測在06時(shí)和18時(shí),所以將在06時(shí)和18時(shí)做36 h預(yù)報(bào),試驗(yàn)時(shí)段是2019年6月1日06時(shí)至29日18時(shí)。

    為了定量評估往返平飄式探空觀測對數(shù)值天氣預(yù)報(bào)的影響,設(shè)計(jì)兩組試驗(yàn):控制試驗(yàn)(Ctrl),同化上述觀測資料;往返平飄式探空觀測影響試驗(yàn)(Effect),除了上述資料外,還加入了構(gòu)造的往返平飄式探空模擬觀測資料(表1)。

    表1 試驗(yàn)方案Table 1 Experiment schemes

    2.2 對降水的影響分析

    本文用ETS(equitable treat score)評分方法對區(qū)域降水預(yù)報(bào)技巧進(jìn)行評價(jià)。圖4給出了06時(shí)起報(bào)的全月平均的逐6 h累計(jì)降水預(yù)報(bào)ETS評分,小雨、中雨、大雨、暴雨和大暴雨的降水閾值分別是0.1、4.0、13.0、25.0和60.0 mm,黃色代表影響試驗(yàn)結(jié)果,紫色代表沒有同化平飄探空的控制試驗(yàn)結(jié)果。圖中方框代表著兩個(gè)試驗(yàn)ETS評分差異的顯著性水平范圍,如果兩個(gè)數(shù)的差在方框之外,就說明通過了顯著性檢驗(yàn)。從圖中可以看出,同化構(gòu)造的往返平飄式探空下降段觀測(06時(shí))對于CMA-MESO模式降水的預(yù)報(bào)評分無論是小雨、中雨、大雨和暴雨都有提高,正效果可維持全預(yù)報(bào)時(shí)段(0~36 h預(yù)報(bào))。其中,小雨量級的ETS降水評分提高較大的是在6~12 h和24~30 h預(yù)報(bào)時(shí)段,分別提高約4%和5%,中雨和大雨預(yù)報(bào)ETS評分提高較大的在18~24 h預(yù)報(bào)時(shí)段,分別提高約8%和10%,暴雨預(yù)報(bào)ETS評分提高較大的在12~18 h預(yù)報(bào)時(shí)段,提高約20%。圖5給出了18時(shí)起報(bào)的全月平均的逐6 h累計(jì)降水預(yù)報(bào)ETS評分。從圖中結(jié)果可知,18時(shí)起報(bào)的影響試驗(yàn)的降水預(yù)報(bào)在各個(gè)降水量級的ETS評分上也均為正效果,其中,小雨量級預(yù)報(bào)降水ETS評分較大的在30~36 h預(yù)報(bào)時(shí)段,提高約8%。中雨量級ETS評分提高較大的在12~18 h和30~36 h預(yù)報(bào)時(shí)段,提高約10%。暴雨量級的ETS評分提高較明顯的是在18~24 h和30~36 h預(yù)報(bào)時(shí)段,提高約50%。綜合上述結(jié)果可知,增加06時(shí)和18時(shí)的往返平飄式探空下降階段的觀測,能夠使得降水預(yù)報(bào)的ETS評分有明顯提高。06時(shí)和18時(shí)比較而言,18時(shí)降水預(yù)報(bào)的ETS評分比06時(shí)更大,并且對于降水量級較大的暴雨的ETS評分提高較大。原因主要是在中國區(qū)域18時(shí)的觀測資料比06時(shí)減少較多,尤其是飛機(jī)資料減少較多,在資料相對少的情況下,加入平飄探空資料更能凸顯觀測的作用。

    圖4 2019年6月1—29日兩組試驗(yàn)06時(shí)起報(bào)的逐6 h累計(jì)降水預(yù)報(bào)月平均ETS評分(a)0~6 h,(b)6~12 h,(c)12~18 h,(d)18~24 h,(e)24~30 h,(f)30~36 h(紫色:控制試驗(yàn)Ctrl,黃色:影響試驗(yàn)Effect)Fig.4 Monthly average ETS scores of 6 h accumulated precipitation forecasts started at 06 UTC by the two experiments from 1 to 29 June 2019(a) 0-6 h, (b) 6-12 h, (c) 12-18 h, (d) 18-24 h, (e) 24-30 h and (f) 30-36 h (purple: control experiment, yellow: impact experiment)

    圖5 同圖4,但為18時(shí)起報(bào)的逐6 h累計(jì)降水預(yù)報(bào)ETS評分Fig.5 Same as Fig.4, but for ETS scores of 6 h accumulated precipitation started at 18 UTC

    2.3 對要素場預(yù)報(bào)的影響

    為了評估增加06時(shí)和18時(shí)的模擬往返探空觀測對要素預(yù)報(bào)的影響,本文選取1000、850、500和250 hPa四個(gè)典型層次,計(jì)算控制試驗(yàn)和影響試驗(yàn)中溫度、濕度和風(fēng)場觀測的0~36 h預(yù)報(bào)與參考場(ERA5再分析)的均方根誤差,利用改進(jìn)率(薛諶彬等,2013)來評估數(shù)值試驗(yàn)預(yù)報(bào)結(jié)果,其中改進(jìn)率I的計(jì)算如下:

    (1)

    式中:t為預(yù)報(bào)時(shí)效,nt為樣本數(shù)量,Ei和ei分別對應(yīng)控制試驗(yàn)和影響試驗(yàn)的均方根誤差。當(dāng)控制試驗(yàn)的均方根誤差大于影響試驗(yàn)的均方根誤差時(shí),It>0,此時(shí)說明影響試驗(yàn)比控制試驗(yàn)有改進(jìn);反之,說明有負(fù)效果。圖6分別是在1000,850,500和250 hPa的U風(fēng)逐6 h預(yù)報(bào)的改進(jìn)率,藍(lán)色代表06時(shí)起報(bào),紅色代表18時(shí)起報(bào)。從結(jié)果可知,改進(jìn)率均是正值,這說明同化06時(shí)和18時(shí)構(gòu)造的模擬觀測的U風(fēng)預(yù)報(bào)相比控制試驗(yàn)有所改進(jìn),并且在250 hPa要比其他層次要更明顯,最高達(dá)2%。圖7是V風(fēng)場的改進(jìn)率,和圖6類似,影響試驗(yàn)的預(yù)報(bào)也有提高,改進(jìn)幅度最高達(dá)2.5%。圖8是溫度場的改進(jìn)率,其改進(jìn)率整體上也是正值,并且在250 hPa 改進(jìn)明顯,改進(jìn)率最高達(dá)5%,但在模式低層850 hPa,18時(shí)起報(bào)之后的18~30 h,改進(jìn)率出現(xiàn)負(fù)值,約為-0.2%,整體負(fù)貢獻(xiàn)較小。比濕的改進(jìn)率和溫度場表現(xiàn)類同,在此不加以贅述。

    圖6 2019年6月1—29日同化模擬平飄探空觀測的U風(fēng)場預(yù)報(bào)改進(jìn)率(a)1000 hPa,(b)850 hPa,(c)500 hPa,(d)250 hPa(藍(lán)色是06時(shí)起報(bào),紅色是18時(shí)起報(bào))Fig.6 The forecast improvement rate of U wind field with assimilation of simulated observation of the return sounding from 1 to 29 June 2019 (a) 1000 hPa, (b) 850 hPa, (c) 500 hPa, (d) 250 hPa(Blue, red lines are beginning forecast times at 06 UTC and 18 UTC, respectively)

    圖7 同圖6,但為V風(fēng)場Fig.7 Same as Fig.6, but for V wind field

    圖8 同圖6,但為溫度場Fig.8 Same as Fig.6, but for temperature field

    2.4 典型個(gè)例預(yù)報(bào)影響分析

    為了進(jìn)一步評估往返平飄式探空模擬觀測對數(shù)值預(yù)報(bào)的影響,選取一個(gè)實(shí)況天氣過程個(gè)例來具體比較影響試驗(yàn)和控制試驗(yàn)。2019年6月4—5日,受低渦切變系統(tǒng)影響,西北地區(qū)東部、四川盆地中東部、山西南部等地出現(xiàn)大到暴雨,四川盆地東部局地出現(xiàn)大暴雨。圖9為2019年6月4日06時(shí)起報(bào)的24 h累計(jì)降水預(yù)報(bào)和累計(jì)降水實(shí)況觀測。從實(shí)況來看主要降水為三個(gè)區(qū)域,分別為西北地區(qū)東部和四川盆地、華南、華北和東北地區(qū),其中西北地區(qū)東部和四川盆地降水最強(qiáng),出現(xiàn)成片的暴雨,華南地區(qū)出現(xiàn)分散性暴雨。對比控制試驗(yàn)和影響試驗(yàn),兩種試驗(yàn)對于西南渦東移影響造成的西北地區(qū)東部至四川盆地一帶的強(qiáng)降水都有較好的體現(xiàn),但對于暴雨雨帶的位置和強(qiáng)度預(yù)報(bào)存在一定的偏差,整體來看對于四川盆地的暴雨雨帶的位置預(yù)報(bào)偏西北,對于甘肅南部至山西中部的暴雨雨帶位置預(yù)報(bào)偏西南,從850 hPa風(fēng)場來看影響試驗(yàn)(圖9b)黑框內(nèi)的偏南急流較控制試驗(yàn)偏弱一些,與實(shí)況更為接近,同時(shí)也可以看到,影響試驗(yàn)較控制試驗(yàn)對于陜西南部的暴雨雨帶預(yù)報(bào)略偏南,這也與實(shí)況雨帶更為接近,而控制試驗(yàn)急流偏強(qiáng)也造成暴雨雨帶略偏北。另外,對于華南的降水(暖區(qū)性質(zhì)降水),廣西南部和廣東中部出現(xiàn)分散性大到暴雨,兩種試驗(yàn)都有一定的體現(xiàn),但影響試驗(yàn)對于華南沿海的偏南氣流較控制試驗(yàn)偏強(qiáng)2 m·s-1(黑框內(nèi)),風(fēng)場的預(yù)報(bào)與實(shí)況更為一致,影響試驗(yàn)對于廣西南部和廣東中部大到暴雨的降水預(yù)報(bào)范圍略大于控制試驗(yàn),位置和強(qiáng)度與實(shí)況更為接近。

    圖9 2019年6月4日06時(shí)起報(bào)24 h累計(jì)降水預(yù)報(bào)(填色)與實(shí)況(圓點(diǎn))的比較(風(fēng)場為6月4日18 UTC)(a)控制試驗(yàn),(b)影響試驗(yàn)(黑色方框代表850 hPa風(fēng)場的偏南急流區(qū)域)Fig.9 Comparison of the 24 h accumulated precipitation forecast (colored) starting at 06 UTC 4 June 2019 with observation data (dot) (wind field: 18 UTC 4 June)(a) control experiment, (b) impact experiment(Black square represents the southerly jet stream zone of the 850 hPa wind field)

    3 結(jié)論與討論

    本文利用往返平飄式探空觀測系統(tǒng)的軌跡模擬系統(tǒng),結(jié)合ERA5再分析資料構(gòu)造往返探空下降段模擬觀測,基于CMA-MESO V5.0高分辨率數(shù)值預(yù)報(bào)系統(tǒng)進(jìn)行了一個(gè)月的觀測系統(tǒng)模擬試驗(yàn),研究了往返平飄式探空在全國組網(wǎng)觀測情況下對數(shù)值預(yù)報(bào)的影響,本文主要結(jié)論如下:

    (1)同化往返平飄式探空模擬觀測能提高模式降水預(yù)報(bào)技巧,表現(xiàn)為在各個(gè)降水量級上的ETS評分上有所提高,小雨、中雨和大雨的提高率在2%~5%,對暴雨的提高率更為明顯,最高可達(dá)50%。

    (2)在要素場的預(yù)報(bào)上,影響試驗(yàn)也比控制試驗(yàn)的預(yù)報(bào)有所改進(jìn),風(fēng)場、溫度場和比濕場基本上都是正效果,風(fēng)場和溫度場的整體改進(jìn)率分別為2%、5%,比濕的整體改進(jìn)率接近2%,其中在250 hPa的風(fēng)場和溫度場的改進(jìn)較其他層次更為明顯。

    (3)從定性分析來看,增加平飄探空下降階段的模擬觀測,能夠更準(zhǔn)確地捕捉降水的動力條件,提高雨帶分布的預(yù)報(bào)準(zhǔn)確度。

    本文通過觀測系統(tǒng)模擬試驗(yàn)論證了增加探空觀測頻次(06時(shí)、18時(shí))、觀測位置對數(shù)值預(yù)報(bào)的正面影響,表明新型往返平飄式探空大規(guī)模組網(wǎng)具備業(yè)務(wù)價(jià)值。文中控制試驗(yàn)的設(shè)置是包含常規(guī)資料和非常規(guī)資料的觀測,影響試驗(yàn)中僅對06時(shí)和18時(shí)的探空位置進(jìn)行模擬觀測的構(gòu)造,受全資料觀測的權(quán)重影響,因而導(dǎo)致要素場預(yù)報(bào)的改進(jìn)效果不是太明顯。此外,本文實(shí)況觀測的預(yù)報(bào)檢驗(yàn)僅是個(gè)例研究,對于本文的結(jié)論是否合理還需要結(jié)合更多的天氣過程進(jìn)行佐證,此外還需依靠將來的實(shí)際平飄探空的觀測數(shù)值試驗(yàn)來進(jìn)行綜合論證。

    致謝:中國氣象局地球系統(tǒng)數(shù)值預(yù)報(bào)中心鄧蓮堂、張進(jìn)在繪圖方面給予的幫助。

    猜你喜歡
    探空降水數(shù)值
    探空氣球?yàn)楦呖諝庀蟆鞍衙}”
    用固定數(shù)值計(jì)算
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    數(shù)值大小比較“招招鮮”
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對比分析
    降水現(xiàn)象儀模擬軟件設(shè)計(jì)與實(shí)現(xiàn)
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    淺談凈舉力對探空氣球升速及施放高度的影響
    河南科技(2014年14期)2014-02-27 14:12:08
    三上悠亚av全集在线观看 | 噜噜噜噜噜久久久久久91| 成人影院久久| 99久国产av精品国产电影| 在线观看www视频免费| 久久国产精品男人的天堂亚洲 | 亚洲丝袜综合中文字幕| 全区人妻精品视频| 亚洲天堂av无毛| 噜噜噜噜噜久久久久久91| 国产 精品1| 又爽又黄a免费视频| 九草在线视频观看| 日韩精品有码人妻一区| 亚洲成人手机| 精品人妻一区二区三区麻豆| 久久久久久久久久久免费av| 99热网站在线观看| 一区二区三区免费毛片| 欧美日韩亚洲高清精品| 亚洲在久久综合| 日日啪夜夜撸| 搡老乐熟女国产| 欧美精品亚洲一区二区| 80岁老熟妇乱子伦牲交| 亚洲真实伦在线观看| av一本久久久久| 国产精品久久久久久精品电影小说| 少妇精品久久久久久久| 国产伦理片在线播放av一区| 国产精品麻豆人妻色哟哟久久| 国产淫片久久久久久久久| 九色成人免费人妻av| 久久久精品94久久精品| 人人澡人人妻人| 能在线免费看毛片的网站| 乱系列少妇在线播放| 亚洲精品456在线播放app| 中文在线观看免费www的网站| 国产成人午夜福利电影在线观看| 少妇高潮的动态图| 人妻夜夜爽99麻豆av| av在线观看视频网站免费| 美女国产视频在线观看| 夫妻午夜视频| 精品卡一卡二卡四卡免费| 日韩 亚洲 欧美在线| 国产爽快片一区二区三区| 另类亚洲欧美激情| 欧美区成人在线视频| freevideosex欧美| 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂| 插阴视频在线观看视频| 9色porny在线观看| 亚洲国产日韩一区二区| 99re6热这里在线精品视频| 97超碰精品成人国产| 午夜视频国产福利| 亚洲精品一区蜜桃| 又黄又爽又刺激的免费视频.| 80岁老熟妇乱子伦牲交| h日本视频在线播放| 天天操日日干夜夜撸| 婷婷色综合大香蕉| 国产亚洲91精品色在线| 午夜免费鲁丝| 下体分泌物呈黄色| 成年av动漫网址| h视频一区二区三区| 亚洲人成网站在线观看播放| 午夜免费观看性视频| 免费看不卡的av| 精品熟女少妇av免费看| 一边亲一边摸免费视频| 免费观看的影片在线观看| 少妇精品久久久久久久| 亚洲精品乱码久久久久久按摩| 日本猛色少妇xxxxx猛交久久| 久久久久精品性色| 日韩在线高清观看一区二区三区| 最黄视频免费看| 午夜福利视频精品| 免费少妇av软件| 交换朋友夫妻互换小说| 校园人妻丝袜中文字幕| 成人18禁高潮啪啪吃奶动态图 | 国产精品人妻久久久影院| av免费在线看不卡| 免费大片18禁| 国产精品久久久久久精品电影小说| 99热6这里只有精品| 三上悠亚av全集在线观看 | 女人精品久久久久毛片| 欧美精品国产亚洲| 久久国产亚洲av麻豆专区| 国产精品熟女久久久久浪| 简卡轻食公司| 免费少妇av软件| 国产精品.久久久| 国内少妇人妻偷人精品xxx网站| 99久久中文字幕三级久久日本| 亚洲精品乱久久久久久| 久久国产亚洲av麻豆专区| av女优亚洲男人天堂| 久久精品久久久久久久性| 国产在线一区二区三区精| 亚洲丝袜综合中文字幕| 肉色欧美久久久久久久蜜桃| 一区在线观看完整版| 久久久国产一区二区| 多毛熟女@视频| 能在线免费看毛片的网站| 欧美3d第一页| 狠狠精品人妻久久久久久综合| 男人添女人高潮全过程视频| 久久这里有精品视频免费| 99re6热这里在线精品视频| 99re6热这里在线精品视频| 大香蕉久久网| av天堂中文字幕网| freevideosex欧美| 99久久精品热视频| 久久人人爽人人爽人人片va| 黄色怎么调成土黄色| 亚洲,一卡二卡三卡| 毛片一级片免费看久久久久| 久久韩国三级中文字幕| 99久久精品国产国产毛片| 丝袜在线中文字幕| 中文在线观看免费www的网站| 黄色一级大片看看| 午夜福利网站1000一区二区三区| 黑人高潮一二区| 国产精品一区二区在线不卡| 最后的刺客免费高清国语| 简卡轻食公司| 亚洲国产最新在线播放| 建设人人有责人人尽责人人享有的| 日韩视频在线欧美| 日韩欧美 国产精品| 黑人高潮一二区| 国模一区二区三区四区视频| 亚洲人与动物交配视频| 中文字幕人妻熟人妻熟丝袜美| 日韩精品免费视频一区二区三区 | 国产成人精品久久久久久| 在线看a的网站| 夜夜看夜夜爽夜夜摸| 亚洲欧美一区二区三区黑人 | 精品国产乱码久久久久久小说| 18禁在线播放成人免费| 国产精品一区www在线观看| 成人美女网站在线观看视频| 99国产精品免费福利视频| 成年女人在线观看亚洲视频| 久久久久久久亚洲中文字幕| 国产色婷婷99| 国产毛片在线视频| 91午夜精品亚洲一区二区三区| 日韩中文字幕视频在线看片| 国产精品久久久久久av不卡| 国产69精品久久久久777片| 久久久欧美国产精品| 久久99一区二区三区| 久久久久精品久久久久真实原创| 丰满乱子伦码专区| 久久韩国三级中文字幕| 另类精品久久| 国产成人午夜福利电影在线观看| 日韩一本色道免费dvd| 午夜久久久在线观看| 一本一本综合久久| 中文字幕人妻熟人妻熟丝袜美| 日韩熟女老妇一区二区性免费视频| 在线观看免费日韩欧美大片 | 国产成人91sexporn| 一区二区三区精品91| 精华霜和精华液先用哪个| 国产亚洲最大av| 亚洲国产色片| 国产免费又黄又爽又色| 精品国产露脸久久av麻豆| 极品人妻少妇av视频| 国产黄片视频在线免费观看| 日韩 亚洲 欧美在线| 国产免费又黄又爽又色| 肉色欧美久久久久久久蜜桃| 日日撸夜夜添| 精品少妇内射三级| 三上悠亚av全集在线观看 | 一级二级三级毛片免费看| 免费黄色在线免费观看| av视频免费观看在线观看| 久久免费观看电影| 国产男人的电影天堂91| 日韩免费高清中文字幕av| 久久99热6这里只有精品| 91精品国产九色| 久久久久久久亚洲中文字幕| 女性生殖器流出的白浆| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院入口| 18禁动态无遮挡网站| 亚洲,欧美,日韩| 亚洲一区二区三区欧美精品| 色哟哟·www| 国产美女午夜福利| 两个人免费观看高清视频 | 日本欧美国产在线视频| 在线观看人妻少妇| 九九在线视频观看精品| 七月丁香在线播放| 午夜老司机福利剧场| 亚洲精品乱码久久久v下载方式| 久久女婷五月综合色啪小说| 99热网站在线观看| 久久国内精品自在自线图片| 国产高清三级在线| 免费av不卡在线播放| 最新的欧美精品一区二区| 国产精品无大码| 日韩av免费高清视频| 丰满乱子伦码专区| 少妇熟女欧美另类| 午夜福利,免费看| 如日韩欧美国产精品一区二区三区 | 18禁动态无遮挡网站| 亚洲性久久影院| 久久久久久久久久人人人人人人| 纵有疾风起免费观看全集完整版| 亚洲综合色惰| 久久国产亚洲av麻豆专区| 99热这里只有精品一区| 波野结衣二区三区在线| 国产成人91sexporn| 日韩免费高清中文字幕av| 午夜影院在线不卡| 另类亚洲欧美激情| 欧美bdsm另类| 精品人妻熟女av久视频| 国产白丝娇喘喷水9色精品| 国产亚洲最大av| 高清在线视频一区二区三区| 在线观看免费日韩欧美大片 | 熟妇人妻不卡中文字幕| 中文欧美无线码| 午夜福利影视在线免费观看| 一个人看视频在线观看www免费| 精品午夜福利在线看| 久久精品国产亚洲网站| 国产精品一区二区性色av| 国产av国产精品国产| 大又大粗又爽又黄少妇毛片口| 国产一区二区在线观看日韩| 妹子高潮喷水视频| 这个男人来自地球电影免费观看 | 亚洲国产最新在线播放| 观看美女的网站| 纵有疾风起免费观看全集完整版| 啦啦啦在线观看免费高清www| 春色校园在线视频观看| 久久久久国产精品人妻一区二区| 久久久久久久国产电影| 国产欧美日韩一区二区三区在线 | 有码 亚洲区| 大香蕉97超碰在线| 精品久久久噜噜| 国产av码专区亚洲av| 亚洲伊人久久精品综合| 亚洲国产精品一区三区| 亚洲国产欧美日韩在线播放 | 美女内射精品一级片tv| 国产日韩欧美在线精品| 在线观看av片永久免费下载| 免费大片18禁| 建设人人有责人人尽责人人享有的| 人体艺术视频欧美日本| 日日爽夜夜爽网站| 久久久久久久久久成人| 伦理电影免费视频| 国产免费又黄又爽又色| 中国三级夫妇交换| 狠狠精品人妻久久久久久综合| 国产极品天堂在线| 五月天丁香电影| 亚洲经典国产精华液单| 国产高清有码在线观看视频| 伦精品一区二区三区| 黑人猛操日本美女一级片| 国产色爽女视频免费观看| 人人妻人人澡人人看| 三级国产精品片| 寂寞人妻少妇视频99o| 久久ye,这里只有精品| 男人添女人高潮全过程视频| 自拍欧美九色日韩亚洲蝌蚪91 | 九色成人免费人妻av| 99久久人妻综合| 国产真实伦视频高清在线观看| 欧美xxxx性猛交bbbb| 国产成人精品无人区| videossex国产| 亚洲欧美成人精品一区二区| 一级片'在线观看视频| 丁香六月天网| 亚洲av男天堂| 国产亚洲av片在线观看秒播厂| 3wmmmm亚洲av在线观看| 国产69精品久久久久777片| 亚洲激情五月婷婷啪啪| 免费黄色在线免费观看| 久久免费观看电影| 欧美日韩av久久| 免费观看av网站的网址| 99久久精品国产国产毛片| 人妻一区二区av| 国产亚洲午夜精品一区二区久久| 一级毛片aaaaaa免费看小| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久精品古装| 日韩一本色道免费dvd| 国产淫片久久久久久久久| 国产有黄有色有爽视频| 人妻人人澡人人爽人人| 99视频精品全部免费 在线| 嫩草影院入口| 色网站视频免费| 国产成人a∨麻豆精品| 精品少妇黑人巨大在线播放| 国产黄片视频在线免费观看| 黄色一级大片看看| 嫩草影院新地址| 成人综合一区亚洲| 日本猛色少妇xxxxx猛交久久| 日本欧美视频一区| 2018国产大陆天天弄谢| 精品视频人人做人人爽| 亚州av有码| 久久久a久久爽久久v久久| 精品一区二区免费观看| 久久精品国产鲁丝片午夜精品| av福利片在线| 久久久久久人妻| 精华霜和精华液先用哪个| 国产精品伦人一区二区| 久久精品国产自在天天线| 香蕉精品网在线| 一区在线观看完整版| 妹子高潮喷水视频| 日韩中字成人| 九九爱精品视频在线观看| 久久精品久久久久久噜噜老黄| 亚洲国产欧美日韩在线播放 | videos熟女内射| 国产欧美日韩综合在线一区二区 | 这个男人来自地球电影免费观看 | 99re6热这里在线精品视频| 偷拍熟女少妇极品色| 日韩成人伦理影院| 久久久久人妻精品一区果冻| 午夜福利,免费看| 18禁在线播放成人免费| 国产成人freesex在线| 久久久久久久精品精品| 性高湖久久久久久久久免费观看| 嫩草影院入口| 岛国毛片在线播放| 国产精品久久久久久av不卡| 精品国产露脸久久av麻豆| 国语对白做爰xxxⅹ性视频网站| 精品国产一区二区三区久久久樱花| 99久久综合免费| 成人免费观看视频高清| 一级毛片我不卡| 国产精品久久久久久久电影| 国内精品宾馆在线| 男的添女的下面高潮视频| 熟女人妻精品中文字幕| 亚洲精品日韩在线中文字幕| 亚洲av综合色区一区| 亚洲欧美精品专区久久| 亚洲国产欧美日韩在线播放 | 国产成人精品一,二区| 成人国产av品久久久| 国产精品人妻久久久久久| 噜噜噜噜噜久久久久久91| 久久亚洲国产成人精品v| 9色porny在线观看| 大香蕉久久网| 亚洲va在线va天堂va国产| 中文字幕人妻丝袜制服| 中文字幕精品免费在线观看视频 | 老司机影院毛片| 一级毛片我不卡| 免费观看av网站的网址| 欧美精品亚洲一区二区| 国产精品成人在线| 婷婷色av中文字幕| 亚洲美女视频黄频| 国产色爽女视频免费观看| 国产男人的电影天堂91| 国产午夜精品一二区理论片| 美女xxoo啪啪120秒动态图| 十八禁网站网址无遮挡 | 美女脱内裤让男人舔精品视频| 精品久久久精品久久久| 国产女主播在线喷水免费视频网站| 国产 精品1| 汤姆久久久久久久影院中文字幕| 晚上一个人看的免费电影| 亚洲欧美成人综合另类久久久| 蜜桃久久精品国产亚洲av| 桃花免费在线播放| 久久婷婷青草| 一级毛片久久久久久久久女| 少妇高潮的动态图| 美女cb高潮喷水在线观看| 午夜福利网站1000一区二区三区| 韩国av在线不卡| 最新的欧美精品一区二区| 亚洲四区av| 国产毛片在线视频| 99久久精品热视频| 亚洲成人一二三区av| 日韩欧美 国产精品| 国产精品欧美亚洲77777| av黄色大香蕉| 亚洲一区二区三区欧美精品| 大香蕉97超碰在线| 看十八女毛片水多多多| 成人18禁高潮啪啪吃奶动态图 | 亚洲av不卡在线观看| 亚洲一区二区三区欧美精品| 老司机影院成人| 久久女婷五月综合色啪小说| 人人妻人人添人人爽欧美一区卜| 久久狼人影院| 久久久久久伊人网av| 午夜精品国产一区二区电影| av播播在线观看一区| 午夜激情福利司机影院| 亚洲精品乱久久久久久| 国产精品熟女久久久久浪| 伊人亚洲综合成人网| 人妻系列 视频| 看非洲黑人一级黄片| videos熟女内射| 国产男女内射视频| 色5月婷婷丁香| 久久国产亚洲av麻豆专区| 国产有黄有色有爽视频| 狂野欧美激情性bbbbbb| 久久国产乱子免费精品| 久久久久久久久久久丰满| 777米奇影视久久| 午夜免费男女啪啪视频观看| 黄色毛片三级朝国网站 | a级毛色黄片| 99热国产这里只有精品6| 老女人水多毛片| 多毛熟女@视频| 久久精品国产亚洲av涩爱| 91精品伊人久久大香线蕉| 亚洲电影在线观看av| 成年人午夜在线观看视频| 只有这里有精品99| 九九爱精品视频在线观看| 日韩中文字幕视频在线看片| 日韩中文字幕视频在线看片| 国产精品不卡视频一区二区| 精品少妇内射三级| 91久久精品电影网| 另类精品久久| 我的女老师完整版在线观看| 亚洲人成网站在线观看播放| 91久久精品国产一区二区成人| 国产熟女欧美一区二区| 亚洲精品成人av观看孕妇| 91aial.com中文字幕在线观看| 熟妇人妻不卡中文字幕| 亚洲精品乱久久久久久| 日韩成人伦理影院| av国产久精品久网站免费入址| 永久免费av网站大全| 欧美3d第一页| 一本—道久久a久久精品蜜桃钙片| av专区在线播放| 天堂8中文在线网| 欧美日韩综合久久久久久| 简卡轻食公司| 看非洲黑人一级黄片| 午夜影院在线不卡| 久久久久人妻精品一区果冻| 国产成人一区二区在线| 天堂8中文在线网| 久久久久国产网址| 好男人视频免费观看在线| 下体分泌物呈黄色| 如何舔出高潮| 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区蜜桃 | 晚上一个人看的免费电影| 欧美bdsm另类| 亚洲不卡免费看| 亚洲综合色惰| 丝瓜视频免费看黄片| 欧美bdsm另类| 亚洲精品国产av蜜桃| 亚洲熟女精品中文字幕| 热re99久久精品国产66热6| 中文欧美无线码| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 国产乱来视频区| 国产一区二区三区av在线| 色视频www国产| 91午夜精品亚洲一区二区三区| 欧美 亚洲 国产 日韩一| 人妻夜夜爽99麻豆av| 99热6这里只有精品| 一区在线观看完整版| 久久久久久久久久久丰满| 国产成人免费观看mmmm| 精品人妻一区二区三区麻豆| 又大又黄又爽视频免费| 久久精品国产鲁丝片午夜精品| 国产成人免费观看mmmm| 99热这里只有是精品50| 亚洲国产欧美在线一区| 日韩电影二区| 国产国拍精品亚洲av在线观看| 中文在线观看免费www的网站| 色婷婷av一区二区三区视频| 黄色视频在线播放观看不卡| 日韩欧美 国产精品| 日本黄色片子视频| av免费观看日本| 国产爽快片一区二区三区| 九九久久精品国产亚洲av麻豆| 国产日韩欧美亚洲二区| 亚州av有码| 两个人免费观看高清视频 | 成人二区视频| 亚洲天堂av无毛| 中文字幕制服av| 国产有黄有色有爽视频| 久久精品夜色国产| 中文欧美无线码| 热re99久久精品国产66热6| 国产精品成人在线| 女的被弄到高潮叫床怎么办| 国内揄拍国产精品人妻在线| 99热这里只有是精品50| 91久久精品国产一区二区三区| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 国产黄片美女视频| 曰老女人黄片| 日本与韩国留学比较| 亚洲国产欧美日韩在线播放 | 极品少妇高潮喷水抽搐| 国产成人免费观看mmmm| 丁香六月天网| 激情五月婷婷亚洲| 三级经典国产精品| 六月丁香七月| 久久女婷五月综合色啪小说| 最后的刺客免费高清国语| 我的老师免费观看完整版| 亚洲欧美日韩另类电影网站| 国产熟女欧美一区二区| 免费观看的影片在线观看| 亚洲真实伦在线观看| 中文字幕亚洲精品专区| 亚洲欧美成人综合另类久久久| 亚洲av成人精品一区久久| 大片电影免费在线观看免费| 中文资源天堂在线| 国产免费视频播放在线视频| 亚洲成人av在线免费| 欧美 亚洲 国产 日韩一| av在线观看视频网站免费| 午夜精品国产一区二区电影| 国产成人免费观看mmmm| 日韩欧美一区视频在线观看 | 国产伦理片在线播放av一区| 男女国产视频网站| 人人妻人人爽人人添夜夜欢视频 | 精品视频人人做人人爽| 18禁在线播放成人免费| 插逼视频在线观看| 狂野欧美激情性bbbbbb| 久久精品国产自在天天线| 中文字幕人妻熟人妻熟丝袜美| 2018国产大陆天天弄谢| 熟女av电影| 日韩成人av中文字幕在线观看| av.在线天堂| 亚洲av综合色区一区| 视频区图区小说| 亚洲第一av免费看| 国产亚洲最大av| 免费看光身美女| 一区二区三区乱码不卡18| 寂寞人妻少妇视频99o| 国产乱来视频区| 亚洲综合色惰| 久久亚洲国产成人精品v| 嫩草影院入口| 国产永久视频网站| 中文字幕制服av| xxx大片免费视频| 亚洲,欧美,日韩| 国产精品一区二区在线观看99|