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

    加密探空資料在華東區(qū)域業(yè)務(wù)模式中的應(yīng)用試驗(yàn)

    2016-10-25 08:35:07徐同王曉峰張蕾楊玉華李佳
    暴雨災(zāi)害 2016年4期
    關(guān)鍵詞:探空增量加密

    徐同,王曉峰,張蕾,楊玉華,李佳

    (1.中國(guó)氣象局上海臺(tái)風(fēng)研究所,上海200030;2.中國(guó)氣象局臺(tái)風(fēng)數(shù)值預(yù)報(bào)重點(diǎn)實(shí)驗(yàn)室,上海200030)

    加密探空資料在華東區(qū)域業(yè)務(wù)模式中的應(yīng)用試驗(yàn)

    徐同1,2,王曉峰1,2,張蕾1,2,楊玉華1,2,李佳1,2

    (1.中國(guó)氣象局上海臺(tái)風(fēng)研究所,上海200030;2.中國(guó)氣象局臺(tái)風(fēng)數(shù)值預(yù)報(bào)重點(diǎn)實(shí)驗(yàn)室,上海200030)

    采用中國(guó)氣象局2014年6月1日—30日14時(shí)加密探空資料,利用華東區(qū)域中尺度數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)比較同化加密探空觀測(cè)資料前后模式預(yù)報(bào)結(jié)果的差異。研究表明,同化加密探空資料后,對(duì)模式初始時(shí)刻不同高度的位勢(shì)高度、比濕、溫度、風(fēng)速等變量均有一定的影響;對(duì)位勢(shì)高度、溫度和風(fēng)場(chǎng)的影響在高層100—150 hPa比較顯著,而對(duì)比濕的影響主要體現(xiàn)在低層700—750 hPa。同化加密探空資料后模式初始場(chǎng)更接近實(shí)況。批量數(shù)值試驗(yàn)的統(tǒng)計(jì)檢驗(yàn)表明,同化加密探空觀測(cè)資料后對(duì)強(qiáng)降水及形勢(shì)場(chǎng)預(yù)報(bào)均有不同程度改進(jìn),24 h暴雨和大暴雨量級(jí)降水的預(yù)報(bào)技巧分別提高了2.5%和8.1%。

    加密探空觀測(cè);華東區(qū)域數(shù)值模式;同化試驗(yàn);數(shù)值模擬

    徐同,王曉峰,張蕾,等.加密探空資料在華東區(qū)域業(yè)務(wù)模式中的應(yīng)用試驗(yàn)[J].暴雨災(zāi)害,2016,35(4):306-314

    XU Tong,WANG Xiaofeng,ZHANG Lei,et al.The application test of intensive radiosonde observations in the Huadong regional numerical model system[J].Torrential Rain and Disaster,2016,35(4):306-314

    引言

    經(jīng)過(guò)一個(gè)世紀(jì)的數(shù)值天氣預(yù)報(bào)理論研究和半個(gè)世紀(jì)的業(yè)務(wù)化應(yīng)用實(shí)踐,數(shù)值天氣預(yù)報(bào)得到了迅速的發(fā)展[1]。多年來(lái),數(shù)值天氣預(yù)報(bào)效果的改善一直圍繞著模式本身和初值進(jìn)行。成功的數(shù)值預(yù)報(bào)離不開精確的模式初值和較好的數(shù)值模式,而目前模式的初值、數(shù)值模式本身都存在誤差(如較大的初值誤差常出現(xiàn)在對(duì)流不穩(wěn)定區(qū)域),因此通過(guò)同化高分辨率、高質(zhì)量的觀測(cè)資料提高模式初值精度就顯得非常重要[2]。全球范圍的探空觀測(cè)大致始于20世紀(jì)40年代[3]。探空資料可以反映大氣的垂直結(jié)構(gòu),在大氣污染擴(kuò)散研究、數(shù)值天氣預(yù)報(bào)、流體動(dòng)力學(xué)分析[4]、模式探空檢驗(yàn)[5]、再分析資料驗(yàn)證[6]及氣候變化研究等方面具有重要作用。

    盡管探空資料站的站點(diǎn)個(gè)數(shù)十分有限,但由于探空觀測(cè)資料質(zhì)量穩(wěn)定、可靠,垂直探測(cè)精度高,可以提供大氣三維結(jié)構(gòu)的完整描述,因此探空資料一直以來(lái)是數(shù)值模式預(yù)報(bào)最基本資料之一[7-8],也是數(shù)值模式預(yù)報(bào)效果和其它觀測(cè)資料可信度檢驗(yàn)的標(biāo)本,研究表明,同化探空觀測(cè)資料可以改進(jìn)大氣溫、濕、風(fēng)廓線的預(yù)報(bào)精度[9],提高對(duì)流天氣預(yù)報(bào)的準(zhǔn)確率[10-12]。然而,目前全球定常探空觀測(cè)一天僅有兩次(00 UTC和12 UTC),很難滿足模式的需求。因此對(duì)探空資料同化在時(shí)間上進(jìn)行加密并分析評(píng)估其對(duì)數(shù)值預(yù)報(bào)的影響是有必要的。

    目前,國(guó)內(nèi)外學(xué)者已開展了一系列加密探空資料在數(shù)值模式中的影響分析。如Faccani C等[13]將2006年非洲季風(fēng)多學(xué)科分析(AMMA)試驗(yàn)計(jì)劃中的加密探空觀測(cè)資料應(yīng)用到法國(guó)全球變分同化系統(tǒng)中,結(jié)果證明使用加密探空資料后,模式模擬的濕度場(chǎng)與觀測(cè)更為接近。李躍清等[14-16]將西南渦加密探空觀測(cè)資料應(yīng)用到數(shù)值預(yù)報(bào)業(yè)務(wù)模式中,研究表明,同化加密探空資料可以改進(jìn)模式初始場(chǎng)質(zhì)量,提高降水預(yù)報(bào)業(yè)務(wù)水平,尤其是區(qū)域性暴雨天氣預(yù)報(bào)。黃燕燕等[17]采用逐時(shí)同化預(yù)報(bào)系統(tǒng)CHAF比較了同化加密探空資料與僅同化常規(guī)探空資料的分析場(chǎng)和模式預(yù)報(bào),研究表明同化加密探空資料后所得分析場(chǎng)水汽的均方根誤差減小、水汽條件和大氣環(huán)流形勢(shì)與實(shí)況更為接近。陳朝平等[18]利用GRAPES_3DVAR系統(tǒng),分別對(duì)2011年發(fā)生在四川境內(nèi)的兩次強(qiáng)降水過(guò)程進(jìn)行加密探空和常規(guī)探空資料的同化對(duì)比試驗(yàn),結(jié)果表明,同化加密探空資料后降水預(yù)報(bào)有所改善,降水強(qiáng)度與降水落區(qū)預(yù)報(bào)與實(shí)況更為接近。盧萍等[19]將2010年在名山、九龍、稻城三地進(jìn)行的西南渦加密觀測(cè)試驗(yàn)取得的探空資料引入到AREM模式中,研究表明,在西南低渦存在的情況下,增加3個(gè)探空站資料后形成的初始場(chǎng)更有利于降水天氣的產(chǎn)生。魏蕾等[20]采用實(shí)時(shí)四維變分同化表明,加入加密探空觀測(cè)后,得到的雷達(dá)回波、高度場(chǎng)、溫度與實(shí)況比較接近,降水也更接近實(shí)況。然而上述研究工作大多基于個(gè)例分析,缺乏批量試驗(yàn)和實(shí)時(shí)業(yè)務(wù)系統(tǒng)的驗(yàn)證。2014年6月,中國(guó)氣象局組織了全國(guó)范圍內(nèi)逐日14時(shí)(北京時(shí),下同)的加密觀測(cè)試驗(yàn)。本文利用華東區(qū)域數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)開展加密探空對(duì)業(yè)務(wù)模式預(yù)報(bào)性能的影響評(píng)估。

    1 模式系統(tǒng)及試驗(yàn)方案設(shè)計(jì)

    華東區(qū)域業(yè)務(wù)模式系統(tǒng)基于中尺度數(shù)值預(yù)報(bào)模式WRF和ARPS Data Analysis System(ADAS)分析同化系統(tǒng)建立。模式中心位于(32°N,117°E),水平分辨率為9 km,水平格點(diǎn)數(shù)為265×265,垂直層次為35層。ADAS是由Oklahoma州立大學(xué)國(guó)家強(qiáng)風(fēng)暴試驗(yàn)室開發(fā)的一套氣象數(shù)據(jù)分析與同化系統(tǒng),從The Local Analysis and Prediction System(LAPS)發(fā)展而來(lái)。ADAS具有非絕熱初始技術(shù),主要是依靠復(fù)雜云分析模塊實(shí)現(xiàn)。

    同化背景場(chǎng)采用GFS模式初始場(chǎng),同化資料主要包括近地面資料(如2 m氣溫、10 m風(fēng)速)、雷達(dá)資料和其它常規(guī)觀測(cè)資料等。ADAS采用的同化分析方法主要為連續(xù)迭代方法(如Bratseth方法)。此外雷達(dá)資料主要使用云分析方法,該方法首先通過(guò)雷達(dá)觀測(cè)資料得到模式格點(diǎn)上的云量和云頂高度,然后根據(jù)三維雷達(dá)反射率資料通過(guò)Smith-Feddes模型,并考慮干空氣夾卷過(guò)程及凝結(jié)過(guò)程對(duì)云水量的損耗,得到云水、云冰混合比;而雨水、雪、冰雹等降水粒子場(chǎng)則是由模式格點(diǎn)上的溫度場(chǎng)和雷達(dá)反射率通過(guò)雷達(dá)反射率方程診斷分析獲得;最后根據(jù)與云水、云冰對(duì)應(yīng)的潛熱釋放,基于濕絕熱對(duì)云內(nèi)溫度場(chǎng)進(jìn)行調(diào)整。同化系統(tǒng)主要參數(shù)配置見表1。

    表1 ADAS同化系統(tǒng)主要參數(shù)配置Table 1 The main parameter configuration of ADAS

    為對(duì)比加密探空觀測(cè)資料對(duì)實(shí)時(shí)業(yè)務(wù)系統(tǒng)的影響,本文選取2014年6月1—30日做為試驗(yàn)樣本進(jìn)行統(tǒng)計(jì)分析,并對(duì)6月19—21日區(qū)域性的暴雨典型過(guò)程進(jìn)行個(gè)例分析。試驗(yàn)分為兩組,一組為控制試驗(yàn),即同化所有常規(guī)觀測(cè)資料和雷達(dá)觀測(cè)資料,另一組為敏感試驗(yàn),在其它條件一致的情況下,在每日14時(shí)增加了加密探空觀測(cè)資料的同化。圖1為模式對(duì)比試驗(yàn)的資料同化流程示意圖。圖2給出了模式區(qū)域范圍內(nèi)加密探空站點(diǎn)位置,增加的加密探空數(shù)量每次約60個(gè),探空數(shù)據(jù)的層次從40多層到60多層均有,平均為50層左右。本文檢驗(yàn)評(píng)估方法參考常見的檢驗(yàn)方法[21-23]。

    2 典型天氣個(gè)例試驗(yàn)分析

    2.1天氣實(shí)況

    2014年6月19—22日,受梅雨鋒和切變線影響,長(zhǎng)江以南地區(qū)發(fā)生一次大范圍降水過(guò)程,華東地區(qū)出現(xiàn)大范圍暴雨天氣。6月19—20日,安徽南部、浙江西北部、江西中部出現(xiàn)暴雨,江西中部和福建西北部局地有大暴雨(圖2a)。6月20—21日,暴雨區(qū)略有北抬東移,主要位于江西北部和皖、蘇、浙交界,大暴雨主要出現(xiàn)在江西北部(圖2b)。6月21—22日,暴雨區(qū)進(jìn)一步東移,主要位于江西西部和浙江境內(nèi),大暴雨主要出現(xiàn)在浙江西部(圖2c)。圖3顯示了本次個(gè)例的過(guò)程實(shí)況累計(jì)降水量分布情況,長(zhǎng)江以南地區(qū)出現(xiàn)大片的50 mm以上降水區(qū)域,其中從江西中西部至浙江西部有一條超過(guò)100 mm的強(qiáng)降水帶。

    2.2同化加密探空觀測(cè)前后模式初始場(chǎng)分析

    2.2.1增量和最大調(diào)整幅度

    為分析加密探空資料對(duì)模式初始場(chǎng)的影響程度,采用增量分析方法(即敏感試驗(yàn)與控制試驗(yàn)的差值)分析各同化變量的變化情況,其中增量為模式范圍內(nèi)的平均結(jié)果。分析發(fā)現(xiàn)從6月19日14時(shí)—6月21日14時(shí)幾個(gè)預(yù)報(bào)時(shí)次的擾動(dòng)振幅分布特征類似。以6月20日14時(shí)起報(bào)個(gè)例為例,可見,位勢(shì)高度的負(fù)增量主要位于500—900 hPa(圖4a),500 hPa向上均為正增量,100 hPa達(dá)到最大正增量值。比濕增量除100 hPa外,其余層次均為負(fù)增量,負(fù)增量最大值位于700 hPa(圖4b)。溫度增量從300—650 hPa為正增量(圖4c),其余層次均為負(fù)增量,正增量最大值在400 hPa,負(fù)增量最大值為925 hPa。u增量的正增量和負(fù)增量從低層到高層交替出現(xiàn),正增量最大值150 hPa,負(fù)增量最大值為500 hPa(圖4d)。

    與增量分析類似,本次個(gè)例同化變量增量的最大調(diào)整幅度(即增量最大值)各個(gè)時(shí)次具有相似的特點(diǎn),因此同樣選擇代表時(shí)次進(jìn)行分析,其它時(shí)次類似。從位勢(shì)高度的增量最大調(diào)整幅度可以看出(圖5a),隨著高度增高,調(diào)整幅度逐漸增加,在100 hPa達(dá)到最大值,總體而言,模式高層的調(diào)整幅度更為顯著。比濕最大調(diào)整幅度調(diào)整較大的層次主要位于1000—650hPa,最大值出現(xiàn)在750 hPa(圖5b),此后隨著高度的增加最大調(diào)整幅度逐漸減小。從溫度最大調(diào)整幅度分析可見(圖5c),調(diào)整幅度較大的層次分布在150 hPa以上和800—1 000 hPa之間,最大值出現(xiàn)在100 hPa。風(fēng)場(chǎng)的水平u分量最大調(diào)整幅度從低層到高層總體呈現(xiàn)增大的趨勢(shì)(圖5d),最大值出現(xiàn)在150 hPa。

    由以上分析可見,同化加密探空觀測(cè)資料后,對(duì)各模式變量在不同高度的初始場(chǎng)均有一定的影響。各模式變量最大調(diào)整幅度分析表明,對(duì)位勢(shì)高度、風(fēng)場(chǎng)和溫度場(chǎng)的影響在高層表現(xiàn)的更為顯著,對(duì)比濕的影響則在低層更為明顯。

    2.2.2增量場(chǎng)水平分布

    圖6給出模式各個(gè)變量增量場(chǎng)的水平分布特征。從圖中可以看到100 hPa位勢(shì)高度增量的水平分布情況(圖6a),模式同化加密探空資料后對(duì)位勢(shì)高度場(chǎng)的影響在整個(gè)積分區(qū)域都有體現(xiàn),增量場(chǎng)有正值和負(fù)值中心存在,正值中心的最大值達(dá)到5 dagpm,負(fù)值中心的最大值為-4 dagpm。同理,對(duì)850 hPa比濕場(chǎng)(圖6b)、100 hPa溫度場(chǎng)(圖6c)和150 hPa緯向風(fēng)(圖6d)的影響也涉及到模式整個(gè)區(qū)域,100 hPa溫度場(chǎng)的增量在模式區(qū)域內(nèi)多為負(fù)值,150 hPa緯向風(fēng)的增量在模式區(qū)域內(nèi)多為正值,且增量中心達(dá)到4 m·s-1。通過(guò)同化變量增量場(chǎng)水平分布的分析可見,同化了加密探空資料之后的影響不僅僅體現(xiàn)在局部地區(qū),而是面向整個(gè)積分區(qū)域的,并有明顯的正負(fù)值中心存在,這可能與探空站點(diǎn)在整個(gè)模式區(qū)域內(nèi)分布均勻有關(guān)。結(jié)合圖8和圖9,模式同化加密探空資料后,在24 h降水預(yù)報(bào)落區(qū)方面較控制試驗(yàn)結(jié)果更為合理,與實(shí)況落區(qū)更為吻合,預(yù)報(bào)技巧也有提高。并且從降水預(yù)報(bào)差值的分布來(lái)看(圖8d),整體是較控制試驗(yàn)增強(qiáng)的。以江西、浙江為例,敏感試驗(yàn)與控制試驗(yàn)的差值中心可達(dá)80mm·(24h)-1以上,對(duì)應(yīng)圖6b可見,敏感試驗(yàn)的850 hPa比濕增量場(chǎng)在相應(yīng)地區(qū)也存在大的正值區(qū)域,有較好的對(duì)應(yīng)。

    2.2.3初始場(chǎng)誤差分析

    將同化前后的分析場(chǎng)檢驗(yàn)變量分別插值到探空站進(jìn)行均方根誤差檢驗(yàn)。圖7為6月20日14時(shí)控制和敏感同化試驗(yàn)不同高度溫度、相對(duì)濕度、u風(fēng)和位勢(shì)高度同化前后的均方根誤差統(tǒng)計(jì)結(jié)果。圖中每一層次的均方根誤差為該層次所有站點(diǎn)均方根誤差平均值。從圖7分析可知敏感試驗(yàn)同化后各層次溫度的均方根誤差基本均小于控制試驗(yàn)(圖7a),相對(duì)濕度和u風(fēng)也有類似的結(jié)果,位勢(shì)高度則主要在500 hPa向上誤差有所減小。

    從以上分析可見,同化加密探空資料之后對(duì)模式初始場(chǎng)的影響從垂直和水平分布均較明顯,并減小了模式初始場(chǎng)誤差。這說(shuō)明加密探空資料中的各層次大氣風(fēng)場(chǎng)、水汽信息等加密觀測(cè)信息已經(jīng)很好地被吸入數(shù)值模式中,由于加密探空資料本身具有較高的垂直分辨率,對(duì)常規(guī)觀測(cè)資料的觀測(cè)盲區(qū)部分起到了很好的補(bǔ)充作用,因此在背景場(chǎng)中有效引入了更多的觀測(cè)信息,從而使得同化分析場(chǎng)更加接近于大氣的真實(shí)狀況。

    2.3降水預(yù)報(bào)對(duì)比分析

    圖8為6月20日14時(shí)降水實(shí)況及兩組試驗(yàn)?zāi)J侥M的24 h降水預(yù)報(bào)對(duì)比。如圖8a所示,6月20日14時(shí)—21日14時(shí),受梅雨鋒影響,華東地區(qū)從江西北部至上海形成一條西南-東北走向的暴雨帶,其中江西東北部、浙江西部和安徽南部交界處有大暴雨,局地特大暴雨。對(duì)比圖8b、8c可以看到,敏感試驗(yàn)和控制試驗(yàn)均預(yù)報(bào)了江西東北部有暴雨和大暴雨出現(xiàn),但控制試驗(yàn)中江西東北部大暴雨的范圍偏小,而敏感試驗(yàn)對(duì)大暴雨的影響范圍預(yù)報(bào)大于控制試驗(yàn),與實(shí)況更為接近,敏感試驗(yàn)和控制試驗(yàn)對(duì)暴雨量級(jí)以上降水均存在明顯的空?qǐng)?bào),特別是大暴雨量級(jí)。由兩試驗(yàn)的24 h· 48 h降水預(yù)報(bào)差值(圖8d)可見,敏感試驗(yàn)比控制試驗(yàn)在江西東北部降水增多,差值超過(guò)80 mm,對(duì)應(yīng)于強(qiáng)降水中心區(qū)域。

    圖9為不同試驗(yàn)?zāi)J?月20日14時(shí)起報(bào)的有無(wú)加密探空觀測(cè)資料24 h累計(jì)降水預(yù)報(bào)的各量級(jí)降水TS評(píng)分對(duì)比。如圖所示,增加加密探空資料同化后,模式對(duì)整個(gè)模式區(qū)域的24 h中雨、暴雨和大暴雨預(yù)報(bào)以及48 h各量級(jí)降水預(yù)報(bào)的TS評(píng)分比無(wú)加密探空資料高,其中對(duì)大暴雨量級(jí)的預(yù)報(bào)改進(jìn)更為明顯。說(shuō)明同化了加密探空資料之后,不僅改善了模式的初始場(chǎng),而且使得模式的降水預(yù)報(bào)在量級(jí)和空間分布上更加接近實(shí)況。

    圖10為不同試驗(yàn)?zāi)J?月20日14時(shí)起報(bào)的有無(wú)加密探空觀測(cè)資料6 h累計(jì)降水預(yù)報(bào)(預(yù)報(bào)至6月20日20時(shí))、18 h累計(jì)降水預(yù)報(bào)(預(yù)報(bào)至6月21日08時(shí))以及24 h累計(jì)降水預(yù)報(bào)(預(yù)報(bào)至4月21日14時(shí))的TS評(píng)分,如圖10a所示,增加加密探空觀測(cè)資料后對(duì)6 h累計(jì)降水預(yù)報(bào)的4 mm·(6 h)-1,13 mm·(6 h)-1和25 mm·(6 h)-1的TS評(píng)分都有提高,但對(duì)0.1 mm·(6 h)-1和1 mm·(6 h)-1預(yù)報(bào)評(píng)分是降低的。分析圖11(b)可知,增加加密探空資料后對(duì)18 h累計(jì)降水預(yù)報(bào)的8 mm·(18 h)-1預(yù)報(bào)技巧基本不變,對(duì)20 mm·(18 h)-1的預(yù)報(bào)技巧有提高,對(duì)0.1 mm·(18 h)-1,1 mm·(18 h)-1和40 mm·(18 h)-1的預(yù)報(bào)技巧比無(wú)加密探空資料低。24 h累計(jì)降水預(yù)報(bào)的10 mm·(24 h)-1,50 mm·(24 h)-1和100 mm·(24 h)-1預(yù)報(bào)技巧都有提高(圖10c),0.1 mm·(24 h)-1和25 mm·(24 h)-1預(yù)報(bào)技巧比無(wú)加密探空資料低。

    3 批量數(shù)值試驗(yàn)統(tǒng)計(jì)結(jié)果分析

    為進(jìn)一步對(duì)比加密探空觀測(cè)對(duì)數(shù)值預(yù)報(bào)實(shí)時(shí)業(yè)務(wù)系統(tǒng)的影響,選取2014年6月1—30日為期一個(gè)月時(shí)間進(jìn)行試驗(yàn)對(duì)比分析。

    3.1降水預(yù)報(bào)結(jié)果檢驗(yàn)

    從降水實(shí)況資料可知,2014年6月模式區(qū)域共有3次主要天氣過(guò)程,分別在6月1—2日,6月19—21日和6月25—26日,降水日共23 d。圖11為2014年6月1—30日為期一個(gè)月的降水預(yù)報(bào)統(tǒng)計(jì)結(jié)果。圖11a為24 h降水預(yù)報(bào)偏差值,可以看到同化加密探空資料后,各個(gè)量級(jí)降水預(yù)報(bào)的預(yù)報(bào)偏差均更加接近于1,進(jìn)一步提高了降水預(yù)報(bào)準(zhǔn)確率。圖11b為24 h累計(jì)降水預(yù)報(bào)TS評(píng)分,可見同化加密探空資料后,24 h暴雨、大暴雨以及48 h小雨、暴雨和大暴雨量級(jí)降水的預(yù)報(bào)技巧均有提高,其中24 h暴雨和大暴雨量級(jí)降水的預(yù)報(bào)技巧分別提高了2.5%和8.1%。

    為了比較同化08時(shí)的探空和14時(shí)的探空對(duì)于午后降水的預(yù)報(bào)能力的影響,對(duì)08時(shí)和14時(shí)起報(bào)的6 h累計(jì)降水(14—20時(shí))進(jìn)行了對(duì)比檢驗(yàn)(如圖12所示),可見同化14時(shí)探空后對(duì)各個(gè)量級(jí)的降水預(yù)報(bào)技巧評(píng)分均較08時(shí)有提高,說(shuō)明同化14時(shí)探空資料是非常有必要的,有助于提升區(qū)域模式對(duì)要素日變化的模擬性能。

    3.2高空預(yù)報(bào)要素檢驗(yàn)

    圖13a為不同試驗(yàn)的500 hPa位勢(shì)高度預(yù)報(bào)場(chǎng)誤差檢驗(yàn)統(tǒng)計(jì)結(jié)果??梢娂用芴娇召Y料后在大多數(shù)情況下均有改進(jìn),其中24 h、48 h位勢(shì)高度場(chǎng)預(yù)報(bào)均方根誤差分別為1.49和1.69;而控制試驗(yàn)的24 h、48 h預(yù)報(bào)均方根誤差分別為1.51和1.70。同樣,700 hPa相對(duì)濕度的預(yù)報(bào)誤差也有一定程度的減小(圖13b,同化加密探空觀測(cè)資料的24 h、48 h相對(duì)濕度預(yù)報(bào)均方根誤差分別為17.63和20.4,而控制試驗(yàn)的24 h、48 h預(yù)報(bào)均方根誤差分別為18.15和20.66。

    從為期一個(gè)月的測(cè)試樣本統(tǒng)計(jì)結(jié)果來(lái)看,同化加密探空觀測(cè)資料對(duì)形勢(shì)場(chǎng)和降水的預(yù)報(bào)都是有正貢獻(xiàn)的,其中位勢(shì)高度的24 h和48 h均方根誤差相對(duì)減小1.3%和0.6%,相對(duì)濕度的24 h和48 h均方根誤差相對(duì)減小2.8%和1.3%。

    4 結(jié)論和討論

    本文利用華東區(qū)域中尺度數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)分析了加密探空觀測(cè)資料對(duì)實(shí)時(shí)業(yè)務(wù)的影響,通過(guò)典型天氣個(gè)例分析和為期一個(gè)月的數(shù)值試驗(yàn),得到如下結(jié)論:

    (1)對(duì)模式初始場(chǎng)增量最大調(diào)整幅度分析結(jié)果表明,同化加密探空觀測(cè)資料后對(duì)模式初始場(chǎng)中不同高度上位勢(shì)高度、比濕、溫度、風(fēng)速均有一定的影響。對(duì)位勢(shì)高度、溫度和風(fēng)場(chǎng)的影響在高層100—150 hPa表現(xiàn)比較顯著,而對(duì)比濕的影響主要集中在低層700—750 hPa。

    (2)同化加密探空資料,可使模式分析場(chǎng)的位勢(shì)高度、相對(duì)濕度、溫度和風(fēng)速各層次的均方根誤差減小,使分析場(chǎng)更加接近于實(shí)況。模擬的暴雨和大暴雨量級(jí)的降水在水平空間分布和量級(jí)上也更接近于實(shí)況,預(yù)報(bào)技巧有所提高。

    (3)批量數(shù)值試驗(yàn)的分析表明,同化加密的探空觀測(cè)資料后對(duì)降水預(yù)報(bào)的提高主要在暴雨和大暴雨,而對(duì)暴雨以下量級(jí)的降水的改進(jìn)并不明顯,相關(guān)問題還有待進(jìn)一步研究。24 h暴雨和大暴雨量級(jí)降水的預(yù)報(bào)技巧分別提高了2.5%和8.1%,對(duì)形勢(shì)場(chǎng)預(yù)報(bào)也有所改善,位勢(shì)高度的24 h和48 h均方根誤差相對(duì)減小1.3%和0.6%,相對(duì)濕度的24 h和48 h均方根誤差相對(duì)減小2.8%和1.3%

    綜上所述,同化加密探空資料,有利于獲得更豐富的大氣觀測(cè)信息,合理利用這些加密觀測(cè)信息可進(jìn)一步改進(jìn)數(shù)值模式初始場(chǎng),在一定程度上提高數(shù)值模式的預(yù)報(bào)效果。由于受資料和試驗(yàn)條件所限,本文僅對(duì)2014年6月一個(gè)月的試驗(yàn)結(jié)果和個(gè)例進(jìn)行統(tǒng)計(jì)分析,今后將繼續(xù)深入開展相關(guān)的研究工作,針對(duì)分析同化加密探空資料對(duì)形勢(shì)場(chǎng)和降水預(yù)報(bào)的改進(jìn)是否具有廣泛適用性的問題作進(jìn)一步研究。

    [1]陳德輝,薛紀(jì)善.數(shù)值天氣預(yù)報(bào)業(yè)務(wù)模式現(xiàn)狀與展望[J].氣象學(xué)報(bào),2004,62(5):624-633

    [2]陳靜,薛紀(jì)善.一種新型的中尺度暴雨集合預(yù)報(bào)初值擾動(dòng)方法研究[J].大氣科學(xué),2005,29(5):717-726

    [3]阮新,朱艷峰,鞠曉慧.融合特性層觀測(cè)資料推算探空站海拔高度的方法[J].高原氣象,2011,30(2):532-537

    [4]劉清惓,戴偉,楊榮康,等.探空溫度測(cè)量太陽(yáng)輻射誤差的流體動(dòng)力學(xué)分析[J].高原氣象,2013,32(4):1157-1164

    [5]李耀孫,石春娥,楊軍,等.我國(guó)東部地區(qū)冬季模式邊界層探空效果評(píng)估[J].高原氣象,2012,31(6):1 690-1 703

    [6]支星,徐海明.三種再分析資料的高空溫度與中國(guó)探空溫度資料的季節(jié)平均特征對(duì)比分析[J].高原氣象,2012,32(1):97-109

    [7]伊里哈木,劉衛(wèi)平,胡義成,等.新疆探空儀器換型前后資料對(duì)比分析[J].綠洲與沙漠氣象,2013,7(2):61-67

    [8]姚作新,呂鳴,賀曉東,等.MP-3000A型地基微波輻射計(jì)在烏魯木齊探空站的輔助性探測(cè)研究[J].綠洲與沙漠氣象,2011,5(1):38-43

    [9]曲曉波,Julian Heming.下投式探空資料對(duì)數(shù)值預(yù)報(bào)初始場(chǎng)影響的個(gè)例分析[J].氣象與環(huán)境學(xué)報(bào),2006,22(2):34-37

    [10]李佳英,俞小鼎,王迎春.用探空資料檢驗(yàn)中尺度數(shù)值模式對(duì)強(qiáng)對(duì)流天氣的診斷分析能力[J].氣象,2006,32(7):13-17

    [11]廖曉農(nóng),俞小鼎,譚一洲.14時(shí)探空在改進(jìn)北京地區(qū)對(duì)流天氣潛勢(shì)預(yù)報(bào)中的作用[J].氣象,2007,33(3):28-32

    [12]魏東,孫繼松,雷蕾,等.三種探空資料在各類強(qiáng)對(duì)流天氣中的應(yīng)用對(duì)比分析[J].氣象,2011,37(4):412-422

    [13]Faccani C,Rabier F,F(xiàn)ourrie N.The impact of AMMA Radiosonde Data on the French Global Assimilation and Forecast system[J].Weather and forecasting,2009,24(5):1 268-1 286

    [14]李躍清,趙興炳,鄧波.2010年夏季西南渦加密觀測(cè)科學(xué)試驗(yàn)[J].高原山地氣象研究,2010,30(4):80-84

    [15]李躍清,趙興炳,張利紅.2011年夏季西南渦加密觀測(cè)科學(xué)試驗(yàn)[J].高原山地氣象研究,2011,31(4):7-11

    [16]李躍清,趙興炳,張利紅,等.2012年夏季西南渦加密觀測(cè)科學(xué)試驗(yàn)[J].高原山地氣象研究,2012,32(4):1-8

    [17]黃燕燕,萬(wàn)齊林,陳子通.加密探空資料在華南暴雨數(shù)值預(yù)報(bào)的應(yīng)用試驗(yàn)[J].熱帶氣象學(xué)報(bào),2011,27(2):179-188

    [18]陳朝平,龍柯吉,張利紅.加密探空資料在GRAPES同化系統(tǒng)中的試驗(yàn)結(jié)果分析[J].高原山地氣象研究,2012,32(4):15-20

    [19]盧萍,鄭偉鵬,趙興炳.川西西南渦加密探空資料分析及數(shù)值模擬實(shí)驗(yàn)[J].高原山地氣象研究,2012,32(1):1-7

    [20]魏蕾,蕾恒池.加密探空資料同化對(duì)一次降水預(yù)報(bào)能力改進(jìn)研究[J].氣候與環(huán)境研究,2012,17(6):809-820

    [21]陳超君,李俊,王明歡.2013年華中區(qū)域中尺度業(yè)務(wù)數(shù)值預(yù)報(bào)的客觀檢驗(yàn)[J].暴雨災(zāi)害,2014,33(2):187-192

    [22]劉靜,葉金印,張曉紅,等.淮河流域汛期面雨量多模式預(yù)報(bào)檢驗(yàn)評(píng)估[J].暴雨災(zāi)害,2014,33(1):58-64

    [23]于曉晶,于志翔,辛渝,等.沙漠綠洲隔壁區(qū)域同化預(yù)報(bào)系統(tǒng)降水預(yù)報(bào)檢驗(yàn)[J].暴雨災(zāi)害,2014,33(3):58-64

    (責(zé)任編輯:鄧雯)

    The application test of intensive radiosonde observations in the East China regional numerical model system

    XU Tong1,2,WANG Xiaofeng1,2,ZHANG Lei1,2,YANG Yuhua1,2,LI Jia1,2
    (1.Shanghai Typhoon Institute of China Meteorological Administration,Shanghai 200030;2.Key laboratory of Numerical Modeling for Tropical cyclone of China Meteorological Administration,Shanghai 200030)

    One-month data assimilation experiments were conducted to examine the impact of assimilating the intensive radiosonde observations in the East China regional numerical model system.Results indicated that assimilating intensive radiosonde observations had influence on geopotential height,specific humidity,temperature and wind speed at different heights in the model initial time.The influence on geopotential height and wind speed is relatively significant at high-level from 100 hPa to 150 hPa.The influence on specific humidity is relativey significant at low level from 700 hPa to 750 hPa.The model initial field obtained from including intensive radiosonde observations was more close to the observed results.The statistical results showed that assimilation of intensive radiosonde observations improves the forecast accuracy for heavy rainfall forecast and synoptic field forecast in some extent.The threat scores of 24 h rainstorm and heavy rain forecast were mproved about 2.5%and 8.1%,respectively.

    intensive radiosonde observation;regional numerical model of the East China;assimilation experiment;numerical simulation

    P435

    A

    10.3969/j.issn.1004-9045.2016.04.002

    2016-03-28;定稿日期:2016-07-01

    十二五”國(guó)家科技支撐計(jì)劃課題(2012BAC21B02);上海市科學(xué)技術(shù)委員會(huì)重點(diǎn)基金(13231203300)

    徐同,主要從事數(shù)值模式檢驗(yàn)和資料影響評(píng)估。E-mail:xut@mail.typhoon.gov.cn

    猜你喜歡
    探空增量加密
    探空氣球?yàn)楦呖諝庀蟆鞍衙}”
    提質(zhì)和增量之間的“辯證”
    “價(jià)增量減”型應(yīng)用題點(diǎn)撥
    一種基于熵的混沌加密小波變換水印算法
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對(duì)比分析
    基于均衡增量近鄰查詢的位置隱私保護(hù)方法
    認(rèn)證加密的研究進(jìn)展
    基于ECC加密的電子商務(wù)系統(tǒng)
    德州儀器(TI)發(fā)布了一對(duì)32位增量-累加模數(shù)轉(zhuǎn)換器(ADC):ADS1262和ADS126
    基于格的公鑰加密與證書基加密
    哪里可以看免费的av片| 成年人黄色毛片网站| 中文字幕精品亚洲无线码一区| 蜜桃久久精品国产亚洲av| 精品国内亚洲2022精品成人| 色综合色国产| 悠悠久久av| 一区福利在线观看| 欧美性猛交╳xxx乱大交人| 在线观看午夜福利视频| 女人被狂操c到高潮| 精品一区二区三区av网在线观看| 日韩,欧美,国产一区二区三区 | 色播亚洲综合网| 狂野欧美白嫩少妇大欣赏| 午夜久久久久精精品| 伊人久久精品亚洲午夜| 欧美又色又爽又黄视频| 欧美精品国产亚洲| 亚洲最大成人手机在线| 午夜爱爱视频在线播放| 别揉我奶头 嗯啊视频| 老司机深夜福利视频在线观看| 国产高潮美女av| 搡女人真爽免费视频火全软件 | 国内毛片毛片毛片毛片毛片| x7x7x7水蜜桃| 午夜视频国产福利| 国产爱豆传媒在线观看| 亚洲18禁久久av| 最近最新免费中文字幕在线| 成人亚洲精品av一区二区| 在线国产一区二区在线| 91久久精品电影网| 天堂√8在线中文| 毛片一级片免费看久久久久 | 动漫黄色视频在线观看| 日本 av在线| 欧美日韩乱码在线| 成年女人毛片免费观看观看9| 麻豆一二三区av精品| 少妇熟女aⅴ在线视频| 少妇熟女aⅴ在线视频| 在线看三级毛片| 午夜福利成人在线免费观看| 精品一区二区三区av网在线观看| 一个人看视频在线观看www免费| 国产一级毛片七仙女欲春2| 日本黄大片高清| 日韩一本色道免费dvd| 国产女主播在线喷水免费视频网站 | 国产高清不卡午夜福利| 女同久久另类99精品国产91| 色综合色国产| 免费观看在线日韩| 精品一区二区三区视频在线观看免费| 麻豆久久精品国产亚洲av| 成人国产一区最新在线观看| 男女边吃奶边做爰视频| 午夜激情欧美在线| 中国美白少妇内射xxxbb| 欧美+亚洲+日韩+国产| 国产av不卡久久| 亚洲经典国产精华液单| 欧美bdsm另类| 亚洲七黄色美女视频| 97超级碰碰碰精品色视频在线观看| 成人美女网站在线观看视频| 我的老师免费观看完整版| 亚洲av免费在线观看| 男人的好看免费观看在线视频| 别揉我奶头 嗯啊视频| 国产白丝娇喘喷水9色精品| eeuss影院久久| 69人妻影院| 蜜桃亚洲精品一区二区三区| 免费无遮挡裸体视频| 黄片wwwwww| 岛国在线免费视频观看| 精品日产1卡2卡| 亚洲国产精品合色在线| 久久精品久久久久久噜噜老黄 | 99久久九九国产精品国产免费| 欧美日韩精品成人综合77777| 天美传媒精品一区二区| 国内久久婷婷六月综合欲色啪| 国产乱人伦免费视频| 精品国产三级普通话版| 国产亚洲精品久久久com| av在线老鸭窝| 日本免费一区二区三区高清不卡| 日韩欧美免费精品| 91麻豆av在线| 亚洲人成网站在线播放欧美日韩| 十八禁网站免费在线| 在现免费观看毛片| 天堂av国产一区二区熟女人妻| 22中文网久久字幕| 一级毛片久久久久久久久女| 国产单亲对白刺激| 亚洲经典国产精华液单| 欧美精品国产亚洲| 少妇人妻精品综合一区二区 | 桃色一区二区三区在线观看| 欧美激情久久久久久爽电影| 色综合亚洲欧美另类图片| 国模一区二区三区四区视频| 啦啦啦观看免费观看视频高清| 欧美一区二区精品小视频在线| 亚洲精品粉嫩美女一区| 国产美女午夜福利| 日本熟妇午夜| 美女黄网站色视频| 可以在线观看的亚洲视频| 小蜜桃在线观看免费完整版高清| 亚洲欧美日韩无卡精品| а√天堂www在线а√下载| 欧美成人免费av一区二区三区| 欧美不卡视频在线免费观看| 久久中文看片网| 久久九九热精品免费| 啦啦啦观看免费观看视频高清| 国产乱人视频| 搡老妇女老女人老熟妇| 国产精品久久电影中文字幕| 在线a可以看的网站| 久久久久久久久大av| 免费观看的影片在线观看| 免费观看精品视频网站| 欧美+亚洲+日韩+国产| 我要搜黄色片| 久久精品国产亚洲av天美| 女同久久另类99精品国产91| a级毛片免费高清观看在线播放| 久久精品国产亚洲网站| 久久精品国产亚洲av涩爱 | av中文乱码字幕在线| 校园人妻丝袜中文字幕| 别揉我奶头 嗯啊视频| 精品一区二区三区视频在线观看免费| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 色吧在线观看| 非洲黑人性xxxx精品又粗又长| 在线观看舔阴道视频| 欧美bdsm另类| 久久九九热精品免费| 小蜜桃在线观看免费完整版高清| 国内揄拍国产精品人妻在线| 99在线视频只有这里精品首页| 天堂av国产一区二区熟女人妻| 国产午夜精品久久久久久一区二区三区 | 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| 久久香蕉精品热| 国产一区二区激情短视频| 国产日本99.免费观看| 伊人久久精品亚洲午夜| 精品午夜福利在线看| 非洲黑人性xxxx精品又粗又长| 亚洲va日本ⅴa欧美va伊人久久| 久久亚洲精品不卡| 国产美女午夜福利| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久| 成人二区视频| 国产午夜福利久久久久久| 成人特级av手机在线观看| 女人被狂操c到高潮| 变态另类成人亚洲欧美熟女| 91久久精品国产一区二区三区| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av| 免费不卡的大黄色大毛片视频在线观看 | 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 国产成人福利小说| 欧洲精品卡2卡3卡4卡5卡区| 12—13女人毛片做爰片一| 在线观看美女被高潮喷水网站| 九九爱精品视频在线观看| 不卡一级毛片| 俄罗斯特黄特色一大片| 最近最新中文字幕大全电影3| 久久久久久九九精品二区国产| 午夜亚洲福利在线播放| 国产在线精品亚洲第一网站| 男人舔奶头视频| 色播亚洲综合网| 成人综合一区亚洲| 午夜福利欧美成人| 乱系列少妇在线播放| 亚洲性久久影院| www.色视频.com| 18禁黄网站禁片午夜丰满| 亚洲国产精品合色在线| 国产男人的电影天堂91| 午夜亚洲福利在线播放| 欧美丝袜亚洲另类 | 国产精品日韩av在线免费观看| 亚洲精品久久国产高清桃花| 精品日产1卡2卡| 黄色配什么色好看| 91av网一区二区| 日韩av在线大香蕉| 亚洲第一电影网av| 国产熟女欧美一区二区| 有码 亚洲区| 无人区码免费观看不卡| a级一级毛片免费在线观看| 午夜福利在线在线| 欧美一区二区亚洲| 无人区码免费观看不卡| 国产伦精品一区二区三区四那| 国产男靠女视频免费网站| 中亚洲国语对白在线视频| 可以在线观看毛片的网站| 美女免费视频网站| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 欧美日韩中文字幕国产精品一区二区三区| 色视频www国产| 亚洲成人免费电影在线观看| 精品免费久久久久久久清纯| 亚洲在线观看片| 此物有八面人人有两片| www.www免费av| 小蜜桃在线观看免费完整版高清| 色吧在线观看| 99热这里只有精品一区| 人人妻,人人澡人人爽秒播| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添小说| 男女下面进入的视频免费午夜| 亚洲国产色片| 熟女电影av网| 一个人观看的视频www高清免费观看| 日韩大尺度精品在线看网址| 婷婷六月久久综合丁香| x7x7x7水蜜桃| 两个人的视频大全免费| 在线观看午夜福利视频| 国产毛片a区久久久久| АⅤ资源中文在线天堂| 久久久久久久久久久丰满 | 国产精华一区二区三区| 欧美日韩瑟瑟在线播放| 亚洲狠狠婷婷综合久久图片| 联通29元200g的流量卡| 97超视频在线观看视频| 色av中文字幕| 2021天堂中文幕一二区在线观| 午夜福利视频1000在线观看| 欧美绝顶高潮抽搐喷水| 久久精品综合一区二区三区| 又黄又爽又刺激的免费视频.| 美女大奶头视频| 欧美色视频一区免费| 亚洲av五月六月丁香网| 久久久久久大精品| 国产探花在线观看一区二区| 国产精品国产三级国产av玫瑰| 成人鲁丝片一二三区免费| 国产午夜福利久久久久久| 99国产精品一区二区蜜桃av| 亚洲性久久影院| 国产国拍精品亚洲av在线观看| 精品日产1卡2卡| 嫩草影视91久久| 一边摸一边抽搐一进一小说| 日韩精品中文字幕看吧| 国产黄色小视频在线观看| 人妻丰满熟妇av一区二区三区| 欧美日韩精品成人综合77777| 欧美xxxx性猛交bbbb| 国产熟女欧美一区二区| 色精品久久人妻99蜜桃| 成人高潮视频无遮挡免费网站| 国产av不卡久久| 高清在线国产一区| 日本 av在线| 色视频www国产| 一进一出抽搐gif免费好疼| 天天躁日日操中文字幕| 少妇高潮的动态图| 超碰av人人做人人爽久久| 大又大粗又爽又黄少妇毛片口| 日韩精品有码人妻一区| 久久草成人影院| 久久精品人妻少妇| 国内精品宾馆在线| 一区二区三区激情视频| 一区二区三区免费毛片| 国产精品一区二区三区四区免费观看 | 狂野欧美激情性xxxx在线观看| 国产综合懂色| 久久99热这里只有精品18| 夜夜夜夜夜久久久久| 国产精品一区二区性色av| 久久天躁狠狠躁夜夜2o2o| 国产伦在线观看视频一区| 国产黄a三级三级三级人| 欧美高清性xxxxhd video| av在线观看视频网站免费| 久久久国产成人免费| 一个人看视频在线观看www免费| 久久精品国产自在天天线| 亚洲成人久久爱视频| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 在线播放无遮挡| 亚洲内射少妇av| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 少妇人妻精品综合一区二区 | 欧美最黄视频在线播放免费| 国产高清视频在线观看网站| av在线老鸭窝| 亚洲精品色激情综合| 欧洲精品卡2卡3卡4卡5卡区| 日本 av在线| 99久国产av精品| 亚洲美女黄片视频| 在线看三级毛片| 亚洲 国产 在线| 国产淫片久久久久久久久| 在线观看免费视频日本深夜| 少妇猛男粗大的猛烈进出视频 | 精品国内亚洲2022精品成人| 国产av不卡久久| 亚洲人成网站高清观看| 丰满人妻一区二区三区视频av| 亚洲avbb在线观看| 日本色播在线视频| 美女高潮的动态| 99国产极品粉嫩在线观看| 国产成人aa在线观看| 国产午夜福利久久久久久| 成人无遮挡网站| 伦理电影大哥的女人| 国产伦人伦偷精品视频| 熟妇人妻久久中文字幕3abv| 校园人妻丝袜中文字幕| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看人在逋| 日韩欧美在线乱码| 听说在线观看完整版免费高清| 日日摸夜夜添夜夜添小说| АⅤ资源中文在线天堂| 国产激情偷乱视频一区二区| 色哟哟·www| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 波多野结衣高清作品| 国产高清三级在线| 天堂√8在线中文| 成人综合一区亚洲| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区久久| 99精品在免费线老司机午夜| 国产精品三级大全| 国产在线男女| or卡值多少钱| 国产成人一区二区在线| 搡老妇女老女人老熟妇| 婷婷亚洲欧美| 在线观看午夜福利视频| 午夜福利在线观看免费完整高清在 | 草草在线视频免费看| 内地一区二区视频在线| 久久久久久久久大av| 国产在视频线在精品| 亚洲av.av天堂| 别揉我奶头 嗯啊视频| 99久久精品热视频| 久久人妻av系列| 91麻豆精品激情在线观看国产| 级片在线观看| 亚洲熟妇中文字幕五十中出| 99热精品在线国产| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 欧美性猛交黑人性爽| 国产精品野战在线观看| 久久久精品欧美日韩精品| 国产探花在线观看一区二区| 亚洲精品国产成人久久av| 亚洲成av人片在线播放无| 日韩强制内射视频| 熟妇人妻久久中文字幕3abv| 免费人成视频x8x8入口观看| 欧美最新免费一区二区三区| 午夜福利在线观看吧| av天堂中文字幕网| 午夜精品在线福利| 国内少妇人妻偷人精品xxx网站| 熟女人妻精品中文字幕| 欧美又色又爽又黄视频| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| 国产久久久一区二区三区| 九九久久精品国产亚洲av麻豆| 日本三级黄在线观看| 一进一出好大好爽视频| 国产aⅴ精品一区二区三区波| 免费不卡的大黄色大毛片视频在线观看 | 国产在线男女| 免费观看人在逋| 欧美日韩瑟瑟在线播放| 久久久久国内视频| 1000部很黄的大片| 久久精品国产亚洲网站| 精品久久久久久久久亚洲 | 亚洲熟妇熟女久久| 亚洲国产日韩欧美精品在线观看| 两个人视频免费观看高清| 精品福利观看| 亚洲欧美清纯卡通| 国产乱人视频| 亚洲精品456在线播放app | 啦啦啦韩国在线观看视频| 国国产精品蜜臀av免费| 搞女人的毛片| 97热精品久久久久久| 日韩中文字幕欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 大型黄色视频在线免费观看| 色视频www国产| eeuss影院久久| 99久久成人亚洲精品观看| 亚洲一级一片aⅴ在线观看| 淫妇啪啪啪对白视频| 九色成人免费人妻av| 赤兔流量卡办理| 欧美性感艳星| 老熟妇乱子伦视频在线观看| 久久九九热精品免费| 大型黄色视频在线免费观看| 91精品国产九色| 男女做爰动态图高潮gif福利片| 久久人妻av系列| 国产成人福利小说| 特级一级黄色大片| 国产精品一及| 精品日产1卡2卡| 国产三级在线视频| 很黄的视频免费| 亚洲一区二区三区色噜噜| 少妇猛男粗大的猛烈进出视频 | 国产国拍精品亚洲av在线观看| 级片在线观看| 他把我摸到了高潮在线观看| 久久精品国产亚洲av天美| av国产免费在线观看| 国内精品久久久久久久电影| 久久热精品热| 日韩 亚洲 欧美在线| 亚洲一区二区三区色噜噜| 亚洲va在线va天堂va国产| 国产乱人伦免费视频| 欧美成人性av电影在线观看| h日本视频在线播放| 男人和女人高潮做爰伦理| 日日啪夜夜撸| 亚洲图色成人| 男女边吃奶边做爰视频| 亚洲国产精品合色在线| 久久国内精品自在自线图片| 97碰自拍视频| 亚洲av电影不卡..在线观看| 欧美精品国产亚洲| 熟妇人妻久久中文字幕3abv| 俺也久久电影网| 又粗又爽又猛毛片免费看| 成人二区视频| 97超视频在线观看视频| 中国美白少妇内射xxxbb| 久久亚洲真实| 一进一出抽搐gif免费好疼| 性欧美人与动物交配| 国内揄拍国产精品人妻在线| 99精品在免费线老司机午夜| 嫩草影院新地址| 在线a可以看的网站| 99久久精品国产国产毛片| 两个人的视频大全免费| 如何舔出高潮| 波多野结衣高清作品| 国产午夜精品论理片| 男女视频在线观看网站免费| 国产精品综合久久久久久久免费| 欧美中文日本在线观看视频| 91麻豆精品激情在线观看国产| 精品久久久久久久人妻蜜臀av| 成人午夜高清在线视频| 我的老师免费观看完整版| 欧美中文日本在线观看视频| 亚州av有码| 在线播放无遮挡| 亚洲精华国产精华精| 日本免费一区二区三区高清不卡| 校园春色视频在线观看| 亚洲精品粉嫩美女一区| 国产精品国产高清国产av| 久9热在线精品视频| 午夜福利18| 老女人水多毛片| 黄色视频,在线免费观看| 一级黄色大片毛片| 热99re8久久精品国产| 午夜a级毛片| 99精品在免费线老司机午夜| 女人十人毛片免费观看3o分钟| 午夜福利在线观看免费完整高清在 | 亚洲最大成人av| 亚洲经典国产精华液单| 美女 人体艺术 gogo| 一a级毛片在线观看| 亚洲欧美日韩高清在线视频| 免费观看人在逋| 在线观看66精品国产| 99热6这里只有精品| 国产主播在线观看一区二区| 成人精品一区二区免费| 啪啪无遮挡十八禁网站| 欧美+亚洲+日韩+国产| 欧美国产日韩亚洲一区| 久久久久久久午夜电影| 亚洲av五月六月丁香网| 国产成人av教育| 久久精品国产亚洲av香蕉五月| 小蜜桃在线观看免费完整版高清| 久久午夜福利片| 国产精品综合久久久久久久免费| 精品久久久久久成人av| 琪琪午夜伦伦电影理论片6080| 国产黄片美女视频| 亚洲人成伊人成综合网2020| 91麻豆精品激情在线观看国产| 亚洲国产日韩欧美精品在线观看| 男人的好看免费观看在线视频| 亚洲自拍偷在线| 十八禁国产超污无遮挡网站| 制服丝袜大香蕉在线| 日本黄大片高清| 亚洲av不卡在线观看| 天堂√8在线中文| 伦理电影大哥的女人| bbb黄色大片| 亚洲精华国产精华液的使用体验 | 亚洲七黄色美女视频| 国产女主播在线喷水免费视频网站 | 久久久久久久亚洲中文字幕| 免费大片18禁| 成年女人看的毛片在线观看| 欧美成人免费av一区二区三区| 国产黄a三级三级三级人| 国产白丝娇喘喷水9色精品| 热99re8久久精品国产| 99热网站在线观看| 搡老熟女国产l中国老女人| 给我免费播放毛片高清在线观看| 五月伊人婷婷丁香| 欧美潮喷喷水| 亚洲精华国产精华精| av中文乱码字幕在线| 99热这里只有是精品50| 女人十人毛片免费观看3o分钟| 久久热精品热| 天堂影院成人在线观看| 亚洲av一区综合| 99久久久亚洲精品蜜臀av| 亚洲久久久久久中文字幕| a级毛片a级免费在线| 69av精品久久久久久| 男女啪啪激烈高潮av片| 黄色丝袜av网址大全| 在线观看美女被高潮喷水网站| 人妻夜夜爽99麻豆av| 国产久久久一区二区三区| av在线天堂中文字幕| 亚洲第一区二区三区不卡| 亚洲av电影不卡..在线观看| a级毛片a级免费在线| 欧美三级亚洲精品| 校园春色视频在线观看| 中文字幕熟女人妻在线| 亚洲中文字幕一区二区三区有码在线看| 久久久久九九精品影院| 日韩欧美精品v在线| 欧美最黄视频在线播放免费| 亚洲精华国产精华精| 久久人人精品亚洲av| 99热精品在线国产| 国模一区二区三区四区视频| 少妇人妻一区二区三区视频| 成人一区二区视频在线观看| 中国美女看黄片| 久久人人精品亚洲av| 成人欧美大片| 日韩欧美精品免费久久| 18禁黄网站禁片免费观看直播| 国产av麻豆久久久久久久| 18禁裸乳无遮挡免费网站照片| 亚洲成人中文字幕在线播放| 亚洲中文日韩欧美视频| 亚洲最大成人手机在线| 久久精品久久久久久噜噜老黄 | 国产精品久久视频播放| 国产真实伦视频高清在线观看 | 少妇猛男粗大的猛烈进出视频 |