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

    往返平飄式探空在CMA-MESO三維變分中的同化及對模式預(yù)報的影響*

    2024-03-15 07:35:44王金成王瑞文
    氣象 2024年2期
    關(guān)鍵詞:下降段探空降水

    王金成 王 丹 王瑞文 譚 娟 容 娜

    1 中國氣象局地球系統(tǒng)數(shù)值預(yù)報中心,北京100081

    2 災(zāi)害天氣國家重點實驗室,北京100081

    3 中國氣象局氣象發(fā)展與規(guī)劃院,北京100081

    4 新疆克拉瑪依氣象局,克拉瑪依834000

    提 要:往返平飄式探空通過一次探空氣球施放實現(xiàn)“上升段-平飄段-下降段”三段觀測,其下降段能實現(xiàn)在06時(世界時,下同)和18時自動垂直加密觀測大氣,具備提升區(qū)域高分辨率快速同化循環(huán)預(yù)報系統(tǒng)在06時和18時的預(yù)報技巧潛力。為了實現(xiàn)往返平飄式探空在區(qū)域高分辨率模式中的同化,分析其對預(yù)報的影響,初步提出了“選取模式層最接近觀測”的垂直稀疏化方法來預(yù)處理資料,深入分析了稀疏化對同化效果的影響,論證了資料垂直稀疏化對于同化應(yīng)用的必要性;在此基礎(chǔ)上,開展了為期1個月的批量同化影響試驗,著重分析了往返平飄式探空在長江中下游區(qū)域的組網(wǎng)觀測對CMA-MESO模式預(yù)報技巧的影響。稀疏化敏感性試驗結(jié)果表明,同化不稀疏化的往返平飄式探空相比同化傳統(tǒng)業(yè)務(wù)探空,分析和預(yù)報誤差顯著增加,降水預(yù)報評分也顯著降低,相反,“選取最接近模式層”數(shù)據(jù)的垂直稀疏化方案能提高模式的分析和預(yù)報技巧,表明往返平飄式探空同化前必須進行垂直稀疏化。批量同化試驗結(jié)果表明,在冷啟動時刻(00時和12時,為常規(guī)探空釋放時刻),同化往返平飄式探空(上升段)相對同化傳統(tǒng)業(yè)務(wù)探空,分析誤差和預(yù)報誤差變化較小。但在暖啟動時刻(03、06、09、15、18、21時,無常規(guī)探空釋放時刻),增加往返平飄式探空下降段數(shù)據(jù),相比控制試驗,分析場精度提高了約0.4%。此外,0~12 h 累計降水預(yù)報的ETS評分變化較小,但12~24 h累計降水預(yù)報在0.1、1.0、5.0、10.0、25.0 mm量級降水ETS評分提高了約0.5%,在50.0 mm量級的降水ETS評分提高了約2.3%。總體而言,同化往返平飄式探空對于區(qū)域高分辨率快速同化循環(huán)預(yù)報系統(tǒng)在暖啟動時刻的降水預(yù)報技巧有正貢獻。

    引 言

    探空觀測是目前全球唯一能夠同時直接觀測高空大氣質(zhì)量場(氣壓、溫度、濕度)和運動場(風(fēng)場)的觀測系統(tǒng),并且其資料質(zhì)量穩(wěn)定可靠、垂直分辨率高(于曉晶等,2018),因此成為了數(shù)值同化應(yīng)用的基礎(chǔ)核心觀測資料(Radnóti et al,2012;郝民等,2014;Ingleby et al,2016)。在資料同化中,探空資料還是衛(wèi)星輻射率和飛機報溫度等資料偏差訂正的錨定資料(Cucurull et al,2014;王瑞文等,2017;吳瓊等,2020)。但受限于人力物力,全球約三分之二的探空仍采用每日定常00時和12時(世界時,下同)兩次的傳統(tǒng)觀測模式(Choi et al,2015),其時間分辨率低,且遠不能滿足中小尺度災(zāi)害天氣系統(tǒng)的監(jiān)測和預(yù)測需求。為此,近年來我國在汛期或重大天氣過程中啟用了06時探空加密觀測模式,并且實現(xiàn)了探空加密觀測在區(qū)域高分辨率模式中的同化應(yīng)用(高篤鳴等,2018;王丹等,2019;李秋陽等,2022),顯著提高了區(qū)域模式在06時起報的強降水預(yù)報技巧。但是,06時探空加密觀測成本巨大,不能維持長久觀測。為了在不顯著增加成本的情況下,實現(xiàn)對06時和18時探空的連續(xù)觀測,滿足數(shù)值天氣預(yù)報對探空觀測的需求,提高數(shù)值天氣預(yù)報對中小尺度災(zāi)害性天氣的預(yù)報技巧,中國氣象局氣象探測中心創(chuàng)新性地提出了往返平飄式探空探測技術(shù)和觀測模式,即通過一次探空氣球施放,獲取“上升段1 h-平飄段4 h-下降段1 h”三段觀測(曹曉鐘等,2019),實現(xiàn)了探空觀測的間接時空加密,獲得的探測垂直廓線數(shù)量加倍,同時能獲取平流層大氣高頻次連續(xù)觀測資料。為了推動往返平飄式探空觀測系統(tǒng)的發(fā)展,開展了模擬觀測系統(tǒng)影響試驗(Observing System Simulation Experiments,OSSEs),試驗結(jié)果表明,往返平飄式探空在全國組網(wǎng)的情況下,其增加的下降段模擬探空觀測,能夠有效提高CMA-MESO在06時和18時起報的降水預(yù)報技巧,不同降水量級的ETS評分提高了約2%~5%(王瑞文等,2023),這一結(jié)論夯實了往返平飄式探空觀測站網(wǎng)建設(shè)的重要性。

    自2018年以來,中國氣象局氣象探測中心在長江中下游地區(qū)建設(shè)了包含6個站的往返平飄式探空區(qū)域觀測網(wǎng),并于2018年、2019年和2021年開展了不同時間長度的區(qū)域組網(wǎng)觀測試驗,獲取了大量的觀測資料。王丹等(2020)針對往返平飄式探空觀測特點,設(shè)計了往返平飄式探空質(zhì)量控制方案,并利用2018年試驗觀測數(shù)據(jù)評估了質(zhì)量控制的效果,分析了往返平飄式探空各變量各階段的不確定性,試驗數(shù)據(jù)的不確定性分析表明往返平飄式探空探測精度達到了世界氣象組織規(guī)定的突破目標,部分探測要素甚至實現(xiàn)了理想目標,探測資料具有良好的可用性,具備很好的應(yīng)用前景。以上研究為往返平飄式探空在數(shù)值天氣預(yù)報中的同化應(yīng)用奠定了基礎(chǔ)。

    然而,往返平飄式探空獲得的高垂直分辨率秒級觀測數(shù)據(jù)的同化技術(shù)和實際試驗觀測數(shù)據(jù)同化對數(shù)值天氣預(yù)報的影響和效益究竟如何?目前而言,這方面的研究相對較少。張旭鵬等(2021)利用WRF(weather research and forecast)模式及WRFDA(WRF data assimilation)同化系統(tǒng)針對長江中下游一次暴雨個例開展了同化往返平飄式探空下降段數(shù)據(jù)試驗,結(jié)果表明同化往返平飄式探空下降段觀測數(shù)據(jù)能夠提高暴雨的預(yù)報技巧。但該研究仍然采用了傳統(tǒng)探空觀測數(shù)據(jù)的稀疏化方案(李慶雷等,2018),并未結(jié)合模式本身的垂直分層特點。為了深入認識往返平飄式探空下降段數(shù)據(jù)對數(shù)值預(yù)報的影響,并實現(xiàn)往返平飄式探空在數(shù)值天氣預(yù)報中的業(yè)務(wù)應(yīng)用,還需要更多的數(shù)值試驗來評估和驗證往返平飄式探空下降段觀測數(shù)據(jù)同化對數(shù)值預(yù)報技巧的影響,且在同化過程中考慮資料的處理策略,特別是結(jié)合模式分層的垂直稀疏化方案。因此,本文首先結(jié)合理論分析和數(shù)值試驗探討了往返平飄式探空同化前開展垂直稀疏化的必要性,提出了一種“選取模式層最接近觀測”的稀疏化方案。進而,利用CMA-MESO快速同化循環(huán)預(yù)報系統(tǒng)開展了為期1個月的批量試驗,系統(tǒng)性評估了往返平飄式探空下降段資料同化對區(qū)域高分辨率數(shù)值天氣預(yù)報的影響。

    1 觀測資料與模式介紹

    選取2021年7月1—31日的往返平飄式探空試驗觀測數(shù)據(jù),共有6個觀測站,分別是安慶、長沙、武漢、贛州、南昌和宜昌,其位置分布如圖1所示。2021年7月,每個站共計開展了62次往返平飄式探空觀測,共成功開展下降段觀測243次,平均單站成功實施觀測約41次。但每次下降段的有效觀測數(shù)據(jù)相對上升段明顯偏少,主要原因是在現(xiàn)有區(qū)域組網(wǎng)條件下,接收站密度不足,加之山地等因素影響導(dǎo)致了下降段信號丟失的高度比較高,最終導(dǎo)致從地面到對流層中層(500 hPa)的下降段觀測較少。

    圖1 往返平飄式探空組網(wǎng)觀測現(xiàn)狀分布(紅點)與傳統(tǒng)業(yè)務(wù)探空站點分布(黑點和紅點)Fig.1 Distribution of round-trip horizontal drift radiosonde network (red dots) and operational radiosonde network (both red and black dots)

    利用CMA-MESO V5.1區(qū)域快速同化循環(huán)預(yù)報系統(tǒng)開展往返平飄式探空同化試驗,采用的同化方法是三維變分同化技術(shù)(3DVar)。該模式目前可同化的觀測資料包括常規(guī)探空、地面、飛機報、船舶、云導(dǎo)風(fēng)、雷達徑向風(fēng)、風(fēng)廓線雷達資料以及GPS/PW觀測。CMA-MESO V5.1版本的同化循環(huán)間隔為3 h,空間分辨率為3 km,垂直層次為50層,模式層頂高為33 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ù)報場。模式技術(shù)細節(jié)可參考黃麗萍等(2022)。

    2 質(zhì)量控制、觀測誤差和垂直稀疏化

    2.1 質(zhì)量控制

    往返探空的質(zhì)量控制方案包含兩部分:一是采用王丹等(2020)研發(fā)的面向往返平飄式探空資料本身的質(zhì)量控制方法,目的是在觀測資料端對試驗資料進行嚴格的質(zhì)量控制,剔除異常數(shù)據(jù)和不滿足同化要求的數(shù)據(jù);二是在同化系統(tǒng)內(nèi)采用背景場檢查方案對質(zhì)量做進一步約束,具體方案如下:

    式中:y是觀測變量,在本文中為溫度、氣壓、濕度和風(fēng);xb是分析變量的背景值;H是觀測算子,表示模式格點向觀測點的水平和垂直插值算法;σa是觀測誤差均方根;σb是觀測變量的背景誤差均方根;α表示閾值參數(shù)。[y-H(xb)]稱為新息增量,一般也用OMB表示。當前,CMA-MESO V5.1 3DVar中的背景檢查質(zhì)量控制中沒有考慮σb的影響,其中,風(fēng)場的閾值參數(shù)α=4.0,溫度的閾值參數(shù)α=2.5,濕度的閾值參數(shù)α=3.5。

    2.2 觀測誤差設(shè)置

    觀測誤差是變分資料同化中的核心關(guān)鍵參數(shù)之一,觀測誤差設(shè)置對模式分析場的精度具有很大的影響。由于目前往返平飄式探空只有6個站的區(qū)域組網(wǎng)觀測資料,尚無法采用精度較高的H-L方法(王金成等,2015)估計其觀測誤差,但已有研究表明往返平飄式探空上升段和下降段各變量的不確定性與傳統(tǒng)探空資料相當(王丹等,2020)。因此本研究仍然采用CMA-MESO 3DVar系統(tǒng)中業(yè)務(wù)探空資料的觀測誤差設(shè)置,觀測誤差如圖2所示。

    圖2 CMA-MESO V5.1 3DVar同化系統(tǒng)中探空觀測資料的觀測誤差(a)溫度,(b)濕度,(c)u風(fēng),(d)v風(fēng)Fig.2 Observation error of radiosonde in 3DVar of CMA-MESO V5.1(a) temperature, (b) humidity, (c) u, (d) v

    2.3 垂直稀疏化

    往返平飄式探空同化前需要進行垂直稀疏化。一方面,往返平飄式探空垂直分辨率與模式分辨率差異較大。往返平飄式探空采用衛(wèi)星導(dǎo)航定位系統(tǒng)(北斗+GPS)定位探空所在空間位置,其位置觀測精度大幅度提高,可以提供秒級的探空定位和觀測數(shù)據(jù),其上升段和下降段垂直分辨率平均約為15 m。目前而言,數(shù)值模式垂直分辨率相對較低,以CMA-MESO V5.1版本為例,其垂直共有50層,模式層頂高度為33 000 m,模式層厚度從地面向高空逐漸增大(如圖3所示),近地面垂直層厚度約為50~100 m,模式層頂附近模式層厚度增大到6.7 km 附近。然而,往返平飄式探空觀測的垂直分辨率(~15 m)明顯高于模式的垂直分辨率,二者相差兩個數(shù)量級。如果不對觀測進行垂直稀疏化,觀測中相對較高分辨率的垂直信息在模式中將被視為噪音,與模式不協(xié)調(diào),造成模式變量間的不協(xié)調(diào),影響預(yù)報效果,也可能引起模式積分不穩(wěn)定。另一方面,往返平飄式探空高垂直分辨率觀測的觀測誤差存在明顯垂直相關(guān)性,而當前主流業(yè)務(wù)同化系統(tǒng)中并未考慮觀測誤差垂直相關(guān)的影響。因此如果往返平飄式探空不進行垂直稀疏化,將在同化系統(tǒng)中過高估計往返平飄式探空的作用,相當于顯著減小觀測誤差,這將導(dǎo)致背景誤差與實際觀測誤差不協(xié)調(diào),分析增量中觀測權(quán)重過高,引起分析誤差增大,也會加重分析場與模式不協(xié)調(diào),引起預(yù)報誤差顯著增加。

    注:模式垂直坐標為高度地形追隨坐標。

    圖4是宜昌、武漢和長沙站稀疏化后業(yè)務(wù)探空資料(黑線)、原始往返平飄式探空(紅線)和稀疏化后往返平飄式探空(藍線)溫度與模式背景差(OMB)的垂直廓線。從圖中可發(fā)現(xiàn),往返平飄式探空原始資料的OMB相比業(yè)務(wù)探空和稀疏化后的往返平飄式探空OMB,高頻噪音明顯較大,特別是高層;此外,往返平飄式探空原始資料OMB在對流層內(nèi)有模式不可分辨的信號,如圖4c中700~500 hPa,未進行稀疏化的資料(紅線)明顯存在很多高頻信號,這些信號并不是觀測噪音,而是觀測中存在模式背景場無法分辨的高頻信號。無論是從模式的理論特征分析,還是結(jié)合實際分析,結(jié)論均認為同化往返平飄式探空前必須進行垂直稀疏化。濕度和風(fēng)場情況與溫度相同(圖略)。

    注:黑線為業(yè)務(wù)探空資料,紅線為往返平飄式探空原始觀測資料,藍線為往返平飄式探空稀疏化后得的資料。

    為了更深入地理解垂直稀疏化的必要性和對預(yù)報的影響,采用數(shù)值試驗方法來進一步定量化驗證垂直稀疏化的重要作用。在模式物理過程和其他參數(shù)不變的情況下,初始場的精度和協(xié)調(diào)性對降水預(yù)報具有決定性作用,因此,稀疏化作用可以通過分析場誤差和降水預(yù)報技巧來體現(xiàn)。為此設(shè)計了三組個例試驗,具體試驗方案如表1所示。本文借鑒掩星觀測(GNSS RO)資料垂直稀疏化的常用方法,即根據(jù)模式層高度,“選取最接近每個模式層”的往返平飄式觀測數(shù)據(jù),因為該稀疏化方法能避免模式層間觀測資料高頻信息對模式初始場造成的不協(xié)調(diào)問題,因此,在GNSS RO折射率和彎角同化中得到廣泛采用(Liu and Xue,2014;Wang et al,2020)。

    表1 數(shù)值試驗方案設(shè)計Table 1 Numerical test schemes

    利用歐洲中期天氣預(yù)報中心的ERA5再分析資料作為參考,計算三組試驗在往返平飄式探空的組網(wǎng)周邊區(qū)域(20°~37°N、105°~125°E)的分析場的均方根誤差,比較三組試驗分析誤差的差異。圖5 是2021年7月2日00時所選區(qū)域內(nèi),試驗ALL與THIN的溫度分析場的均方根誤差與試驗CTRL的相對差異。從圖5中可以發(fā)現(xiàn),試驗ALL溫度分析誤差在750~450 hPa和250 hPa以上明顯大于試驗CTRL,相對增大2%左右,而試驗THIN溫度分析誤差與控制試驗CTRL的分析基本相當。結(jié)果表明,同化原始高垂直分辨率的往返平飄式探空,分析誤差會有不同程度的增大;然而,對往返平飄式探空進行稀疏化后,溫度分析誤差會有明顯減小。其他變量結(jié)果與溫度變量結(jié)果相同,在此不再贅述。上述結(jié)果進一步說明往返平飄式探空高垂直分辨率資料同化前需要進行垂直稀疏化。

    注:黑色虛線是零參考線。

    采用全國10 796個國家級自動觀測站的降水觀測數(shù)據(jù),對上述三組試驗進行降水評分檢驗,降水檢驗區(qū)域與誤差分析的區(qū)域相同。圖6是三組試驗在2021年7月2日00時起報的24 h累計降水ETS評分。從圖中可以發(fā)現(xiàn),試驗ALL預(yù)報的24 h 累計降水ETS評分明顯低于試驗CTRL和THIN,而試驗THIN與CTRL預(yù)報的24 h累計降水ETS評分基本相當。該結(jié)果表明,往返平飄式探空不進行垂直稀疏化會明顯降低模式降水評分,進一步佐證了理論分析結(jié)果,說明往返平飄式探空這類高垂直分辨率觀測資料同化前需要進行垂直稀疏化。

    注:模式起報時間為2021年7月2日00時。

    此外,“選取與模式層接近”的往返平飄式探空觀測資料用于同化的稀疏化方案效果明顯,可以用于往返平飄式探空上升段和下降段資料同化?;诖?在CMA-MESO V5.1 3DVar同化系統(tǒng)中采用該稀疏化方案,對往返平飄式探空開展批量同化試驗,評估往返平飄式探空下降段資料對數(shù)值天氣預(yù)報的影響,為實現(xiàn)其在業(yè)務(wù)中的應(yīng)用提供量化支撐。

    3 批量試驗結(jié)果分析

    為了明確往返探空資料同化為數(shù)值天氣預(yù)報帶來的效益,采用第2節(jié)中的質(zhì)量控制、稀疏化和觀測誤差,利用CMA-MESO V5.1快速循環(huán)同化預(yù)報系統(tǒng),進行了2021年7月1—31日為期1個月的批量試驗。具體試驗設(shè)計分別為表1中的控制試驗CTRL和影響試驗THIN。

    對比分析00時和12時,同化往返平飄式探空和同化業(yè)務(wù)探空的分析和預(yù)報結(jié)果非常接近,兩者差異在0.1%以內(nèi),這從側(cè)面表明往返平飄式探空上升段資料(00時和12時對應(yīng)上升段觀測)的觀測精度和業(yè)務(wù)探空相當。由于往返平飄式探空相比業(yè)務(wù)探空增加了下降段廓線觀測,因此本研究中著重分析下降段往返平飄式探空對分析和預(yù)報的影響。

    圖7是CMA-MESO V5.1暖啟動時刻(03、06、09、15、18、21時)同化的往返式探空觀測數(shù)據(jù)量,分別對應(yīng)3 h同化時間窗內(nèi)往返平飄式探空下降段數(shù)據(jù)量。從圖7中可以發(fā)現(xiàn),經(jīng)過稀疏化和質(zhì)量控制后,能夠同化的往返平飄段下降段數(shù)據(jù)十分有限,每3 h平均同化數(shù)據(jù)量僅有64個數(shù)據(jù),也即相比業(yè)務(wù)探空,每天增加的數(shù)據(jù)量平均為384個數(shù)據(jù),基本相當于2條完整探空站的同化觀測數(shù)量。經(jīng)過分析,這是由兩個方面原因造成的:一是7月觀測試驗中6個站成功下降的觀測僅有243次,也就是成功下降觀測的比率為65.3%,二是下降段觀測的數(shù)據(jù)結(jié)束高度較高,平均僅到500 hPa附近高度。目前CMA-MESO V5.1的模式垂直層次設(shè)計呈現(xiàn)“低層密集、高層稀疏”的特點,低層探測數(shù)據(jù)不足加上稀疏化策略的影響,最終導(dǎo)致經(jīng)過稀疏化后的觀測數(shù)據(jù)量就變得十分有限。

    圖7 暖啟動時刻(03、06、09、15、18、21時)每3 h同化的往返平飄式探空數(shù)據(jù)量(包括溫度、濕度、u風(fēng)和v風(fēng))Fig.7 Assimilated observation number of round-trip horizontal drift sounding data (temperature, humidity, u and v) every 3 hours in 3DVar at warm start time (03 UTC, 06 UTC, 09 UTC, 15 UTC, 18 UTC, 21 UTC)

    3.1 往返平飄式探空同化對高空預(yù)報的影響

    圖8是06時試驗THIN高空場分析和預(yù)報均方根誤差相對試驗CTRL的分析和預(yù)報均方根誤差的減小率,負值表示THIN的均方根誤差相比CTRL減小,正值表示THIN的均方根誤差相比CTRL增大。從圖中可以看出,試驗THIN各變量分析和預(yù)報誤差與試驗CTRL相比差異非常小,差異僅在-0.5%~0.5%。并且試驗THIN相比CTRL分析和預(yù)報誤差差異較大的層次在600 hPa以上,這和往返平飄式探空下降段數(shù)據(jù)主要在500 hPa 以上有關(guān)。其他時刻結(jié)果與06時的結(jié)果十分相似,此處圖略。結(jié)果表明,往返平飄式探空下降段對CMA-MESO高空場分析和預(yù)報的影響相對較小,整體誤差減小率僅為0.1%。

    注:負值表示誤差減小,正值表示誤差增大,黑線是分析時刻,彩線代表不同預(yù)報時次,此處為06時起報結(jié)果。

    3.2 往返平飄式探空同化對降水預(yù)報的影響

    圖9是暖啟動時刻(03、06、09、15、18、21時)起報的結(jié)果,計算了試驗THIN預(yù)報的12 h累計降水ETS評分相比試驗CTRL的改進率,正值表示THIN相比CTRL試驗降水預(yù)報技巧有提高,負值則為降低。從圖9a中可以發(fā)現(xiàn),暖啟動時刻起報的0~12 h累計降水,試驗THIN在量級為0.1、1.0、50.0 mm的降水預(yù)報ETS評分雖有改善,但是變化幅度平均小于0.04%,在5.0、10.0、25.0 mm的降水預(yù)報ETS評分則略有降低,變化幅度僅為0.07%??傮w而言,同化往返平飄式探空6個站的試驗數(shù)據(jù)對0~12 h的降水預(yù)報影響非常小。查看暖啟動時刻起報的12~24 h累計降水評分結(jié)果,從圖9b結(jié)果來看:同化往返平飄式探空6個站的試驗數(shù)據(jù)對12~24 h累計降水預(yù)報在所有量級降水預(yù)報ETS評分有所增加,ETS評分平均提高了0.7%,特別是在50.0 mm降水量級,12~24 h降水的ETS評分提高了2.2%。整體而言,往返探空資料同化的正貢獻在模式暖啟動的12~24 h預(yù)報時段更為明顯。

    注:正值表示降水預(yù)報技巧有提高,負值則為降低。

    4 結(jié)論與討論

    通過理論分析和個例試驗,本文論證了往返平飄式探空同化前必須進行垂直稀疏化,并初步提出了“選取模式層最近觀測數(shù)據(jù)”的稀疏化方法。利用CMA-MESO V5.1 快速循環(huán)同化預(yù)報系統(tǒng),開展了2021年7月1—31日為期1個月的同化預(yù)報影響試驗,重點分析了同化2021年7月長江中下游6個試驗站網(wǎng)的下降段觀測數(shù)據(jù)對CMA-MESO V5.1暖啟動時刻分析和預(yù)報的影響。主要結(jié)論如下:

    (1)無論從理論分析還是實際同化角度,往返平飄式探空同化前均需要進行垂直稀疏化?!斑x取距離模式層最接近的觀測”是一種相對簡潔、有效的稀疏化方案。

    (2)由于試驗數(shù)據(jù)在下降段相對較少,同化增加的觀測數(shù)量相對有限。具體表現(xiàn)為:2021年7月長江中下游地區(qū)往返平飄式探空試驗觀測數(shù)據(jù)下降段成功觀測率約為65%,下降段觀測的結(jié)束高度較高,經(jīng)過稀疏化和質(zhì)量控制,相比業(yè)務(wù)探空,平均每天同化數(shù)據(jù)量增加384個數(shù)據(jù),僅相當于2條上升段廓線數(shù)量。

    (3)在同化往返平飄式探空資料量相對稀少的下降段觀測數(shù)據(jù)條件下,對暖啟動時刻高空場的分析和預(yù)報影響相對較小,誤差變化率在-0.5%~0.5%。

    (4)在同化往返平飄式探空資料量相對稀少的下降段觀測數(shù)據(jù)條件下,暖啟動時刻0~12 h累計降水ETS評分變化較小(略減小0.07%),12~24 h不同量級累計降水ETS評分均有提高,平均提高了0.7%,50.0 mm以上降水ETS評分提高了2.2%。

    上述研究結(jié)果表明,在目前觀測試驗下降段資料量有限和當前稀疏化策略下,盡管同化往返平飄式探空的下降段觀測數(shù)據(jù)量有限,但是還能在一定程度上提高高分辨率區(qū)域模式的降水預(yù)報技巧。結(jié)合往返平飄式探空觀測系統(tǒng)模擬試驗(OSSEs)結(jié)果(王瑞文等,2023),一旦實現(xiàn)往返平飄式探空組網(wǎng)觀測,CMA-MESO快速循環(huán)同化預(yù)報系統(tǒng)暖啟動時刻的降水預(yù)報技巧還能更為進步。

    目前而言,往返平飄式探空還未完全發(fā)揮在數(shù)值天氣預(yù)報方面的效益,在觀測和同化技術(shù)角度還有很多方面需要改進,主要包括三個方面:

    (1) 提高往返平飄式探空下降段觀測成功率,布設(shè)更為完備的信號接收網(wǎng),盡可能地接收到下降段探空從高空到地面完整的廓線觀測數(shù)據(jù)。

    (2) 研發(fā)更優(yōu)的垂直稀疏化方案,最大限度地提高同化往返平飄式數(shù)據(jù)量,從而同化更多的下降段數(shù)據(jù)。在導(dǎo)航衛(wèi)星掩星觀測同化中,已有學(xué)者提出了較好的稀疏化方案,這為往返平飄式探空稀疏化提供了很多參考,如馬旭林等(2019)針對GPS掩星反演資料同化提出了一種適應(yīng)數(shù)值模式垂直坐標的稀疏化方案,取得了很好的效果,后續(xù)將在往返平飄式探空同化中嘗試該方法并分析評估該方法的效果。

    (3)研究往返平飄式探空觀測誤差高精度估計方法,提高往返平飄式探空同化中的觀測誤差的精度,進一步提高資料的同化效果。

    猜你喜歡
    下降段探空降水
    探空氣球為高空氣象“把脈”
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    月球探測器動力下降段最優(yōu)軌跡參數(shù)化方法
    “阿波羅”是如何返回地球的?
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    球載式下投國產(chǎn)北斗探空儀測風(fēng)性能評估
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對比分析
    降水現(xiàn)象儀模擬軟件設(shè)計與實現(xiàn)
    循環(huán)流化床中顆粒振蕩循環(huán)現(xiàn)象的實驗研究
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    国产欧美亚洲国产| 在线精品无人区一区二区三 | 我要看黄色一级片免费的| 亚洲高清免费不卡视频| 日韩大片免费观看网站| 最黄视频免费看| 亚洲美女黄色视频免费看| 身体一侧抽搐| 成人亚洲精品一区在线观看 | 日韩中字成人| 人妻少妇偷人精品九色| 亚洲国产精品成人久久小说| 亚洲精品国产av成人精品| 多毛熟女@视频| 一边亲一边摸免费视频| 我的女老师完整版在线观看| 日本-黄色视频高清免费观看| 韩国高清视频一区二区三区| 黄色怎么调成土黄色| 亚洲精品国产成人久久av| 欧美精品一区二区大全| 一区二区三区免费毛片| 免费观看无遮挡的男女| 久久久久人妻精品一区果冻| 大码成人一级视频| 成人特级av手机在线观看| 久久人人爽人人爽人人片va| 人人妻人人添人人爽欧美一区卜 | 久久久久性生活片| 在线免费观看不下载黄p国产| 女人十人毛片免费观看3o分钟| 18禁在线无遮挡免费观看视频| 最近中文字幕高清免费大全6| 免费大片18禁| 91精品伊人久久大香线蕉| 欧美成人a在线观看| 毛片女人毛片| 成人综合一区亚洲| 大话2 男鬼变身卡| 国产男女内射视频| 简卡轻食公司| 国产成人午夜福利电影在线观看| 国产成人精品福利久久| 成人毛片a级毛片在线播放| 26uuu在线亚洲综合色| 日韩伦理黄色片| 18禁在线播放成人免费| 国产69精品久久久久777片| a级一级毛片免费在线观看| 一级毛片黄色毛片免费观看视频| 亚洲自偷自拍三级| 亚洲四区av| 中文字幕人妻熟人妻熟丝袜美| 性高湖久久久久久久久免费观看| 亚洲伊人久久精品综合| 国语对白做爰xxxⅹ性视频网站| 精品99又大又爽又粗少妇毛片| 中国美白少妇内射xxxbb| 久久精品人妻少妇| 亚洲精品国产色婷婷电影| 亚洲国产色片| 亚洲三级黄色毛片| 欧美精品人与动牲交sv欧美| 欧美少妇被猛烈插入视频| 男人爽女人下面视频在线观看| 一个人免费看片子| 成人毛片60女人毛片免费| 国产欧美另类精品又又久久亚洲欧美| 2021少妇久久久久久久久久久| 国产 精品1| 成人二区视频| 简卡轻食公司| 亚洲精品第二区| 免费在线观看成人毛片| 三级国产精品片| 国产精品久久久久久精品古装| 97在线人人人人妻| 欧美日韩综合久久久久久| 少妇人妻 视频| 久久国产亚洲av麻豆专区| 啦啦啦啦在线视频资源| 少妇高潮的动态图| 精品一区二区免费观看| 中文字幕精品免费在线观看视频 | 国产成人91sexporn| 美女福利国产在线 | 国产精品国产三级专区第一集| 男女国产视频网站| 日本黄色片子视频| 日韩在线高清观看一区二区三区| 久久国产乱子免费精品| 少妇猛男粗大的猛烈进出视频| 国产精品蜜桃在线观看| 99国产精品免费福利视频| 亚洲美女黄色视频免费看| 观看美女的网站| 免费看日本二区| 国产成人freesex在线| 免费观看a级毛片全部| 免费观看性生交大片5| 久久久久久久大尺度免费视频| 久久久久久久大尺度免费视频| 青春草国产在线视频| 大片免费播放器 马上看| 日韩一区二区三区影片| 一个人免费看片子| 亚洲在久久综合| 亚洲av中文字字幕乱码综合| 赤兔流量卡办理| 80岁老熟妇乱子伦牲交| 久久亚洲国产成人精品v| 中文资源天堂在线| 少妇高潮的动态图| 久久久久久九九精品二区国产| 精华霜和精华液先用哪个| 亚洲天堂av无毛| 久热这里只有精品99| 日日摸夜夜添夜夜添av毛片| 亚洲精品中文字幕在线视频 | 黑人猛操日本美女一级片| 女人十人毛片免费观看3o分钟| 2021少妇久久久久久久久久久| 日日啪夜夜爽| 色综合色国产| 亚洲av日韩在线播放| 国产成人91sexporn| 亚洲伊人久久精品综合| 久久99热这里只有精品18| 啦啦啦在线观看免费高清www| 国产精品三级大全| av在线播放精品| 国产在视频线精品| 三级国产精品片| av在线app专区| 午夜免费观看性视频| 日韩欧美 国产精品| 精品午夜福利在线看| 简卡轻食公司| 不卡视频在线观看欧美| 久久久久视频综合| 少妇 在线观看| 人妻一区二区av| 看非洲黑人一级黄片| 99久久综合免费| av在线老鸭窝| av播播在线观看一区| 免费大片18禁| 涩涩av久久男人的天堂| 男人舔奶头视频| 观看美女的网站| 国产深夜福利视频在线观看| av网站免费在线观看视频| 激情 狠狠 欧美| 国产精品蜜桃在线观看| 精品熟女少妇av免费看| 久久综合国产亚洲精品| 中文在线观看免费www的网站| 1000部很黄的大片| 亚洲色图综合在线观看| 少妇精品久久久久久久| 免费高清在线观看视频在线观看| 97热精品久久久久久| 免费少妇av软件| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 欧美一级a爱片免费观看看| 国产乱人视频| 国产一区亚洲一区在线观看| 下体分泌物呈黄色| 成人免费观看视频高清| 一级毛片电影观看| 特大巨黑吊av在线直播| 国内少妇人妻偷人精品xxx网站| 亚洲av电影在线观看一区二区三区| 亚洲精品国产av蜜桃| 精品国产一区二区三区久久久樱花 | 国产一区二区在线观看日韩| 婷婷色综合www| 99久久精品一区二区三区| 国产 一区精品| 97热精品久久久久久| 日韩欧美一区视频在线观看 | 天堂8中文在线网| 欧美极品一区二区三区四区| 22中文网久久字幕| 99久久综合免费| h视频一区二区三区| 深夜a级毛片| 能在线免费看毛片的网站| 成人午夜精彩视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩卡通动漫| 91精品国产国语对白视频| 欧美日韩在线观看h| 99热国产这里只有精品6| av一本久久久久| 国模一区二区三区四区视频| 亚洲婷婷狠狠爱综合网| 男人狂女人下面高潮的视频| 不卡视频在线观看欧美| 色哟哟·www| 香蕉精品网在线| 精品熟女少妇av免费看| 一区二区三区四区激情视频| 午夜视频国产福利| 亚洲欧洲国产日韩| 简卡轻食公司| 国产有黄有色有爽视频| 啦啦啦中文免费视频观看日本| 99热这里只有是精品50| 欧美bdsm另类| 少妇人妻精品综合一区二区| 欧美三级亚洲精品| 在线观看免费高清a一片| 联通29元200g的流量卡| 亚洲av.av天堂| 麻豆乱淫一区二区| 大话2 男鬼变身卡| 欧美xxxx黑人xx丫x性爽| 蜜桃在线观看..| 亚洲欧洲国产日韩| 啦啦啦啦在线视频资源| 在线观看av片永久免费下载| 亚洲va在线va天堂va国产| 亚洲无线观看免费| 一边亲一边摸免费视频| 又爽又黄a免费视频| 日韩人妻高清精品专区| 日本欧美国产在线视频| 波野结衣二区三区在线| 丰满迷人的少妇在线观看| 中文字幕人妻熟人妻熟丝袜美| 欧美亚洲 丝袜 人妻 在线| 欧美成人一区二区免费高清观看| 亚洲,欧美,日韩| 777米奇影视久久| 国产亚洲一区二区精品| 亚洲av欧美aⅴ国产| 最近中文字幕高清免费大全6| 青春草视频在线免费观看| 日韩伦理黄色片| 精品人妻熟女av久视频| 日本免费在线观看一区| 毛片一级片免费看久久久久| 不卡视频在线观看欧美| 久久影院123| 国产色爽女视频免费观看| 亚洲欧美日韩无卡精品| 亚洲国产精品成人久久小说| 你懂的网址亚洲精品在线观看| 偷拍熟女少妇极品色| 久久久a久久爽久久v久久| 中文欧美无线码| 日本欧美国产在线视频| 波野结衣二区三区在线| 亚洲,一卡二卡三卡| 日韩一区二区视频免费看| 亚洲av成人精品一二三区| 少妇熟女欧美另类| 交换朋友夫妻互换小说| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 水蜜桃什么品种好| 大话2 男鬼变身卡| 亚洲精品国产av成人精品| 亚洲国产日韩一区二区| 久久人人爽av亚洲精品天堂 | 亚洲欧美日韩无卡精品| 韩国av在线不卡| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久av不卡| 成人高潮视频无遮挡免费网站| 日本猛色少妇xxxxx猛交久久| 九九爱精品视频在线观看| 日日啪夜夜撸| 丝袜脚勾引网站| 国产高清不卡午夜福利| 卡戴珊不雅视频在线播放| 欧美老熟妇乱子伦牲交| 亚洲精品国产成人久久av| tube8黄色片| 99热国产这里只有精品6| 日韩一区二区三区影片| 天天躁夜夜躁狠狠久久av| 丰满少妇做爰视频| h日本视频在线播放| 99九九线精品视频在线观看视频| 国产在线男女| 亚洲国产欧美人成| 男女无遮挡免费网站观看| 成人国产麻豆网| 麻豆精品久久久久久蜜桃| 三级国产精品欧美在线观看| 2021少妇久久久久久久久久久| 亚洲色图av天堂| 国产综合精华液| 秋霞伦理黄片| 少妇的逼好多水| 尤物成人国产欧美一区二区三区| 看非洲黑人一级黄片| 国产视频内射| 午夜福利在线在线| a 毛片基地| 免费观看av网站的网址| 亚洲欧美清纯卡通| 国产熟女欧美一区二区| videossex国产| 国产一区二区三区综合在线观看 | 黄色配什么色好看| 成人漫画全彩无遮挡| 日韩电影二区| 日韩av免费高清视频| 一个人免费看片子| 一本色道久久久久久精品综合| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 少妇熟女欧美另类| 日产精品乱码卡一卡2卡三| 18禁在线无遮挡免费观看视频| 国产日韩欧美在线精品| 久久久成人免费电影| 欧美 日韩 精品 国产| 欧美亚洲 丝袜 人妻 在线| 国产人妻一区二区三区在| 伊人久久精品亚洲午夜| 久久久久视频综合| 国产亚洲av片在线观看秒播厂| 日产精品乱码卡一卡2卡三| 免费观看无遮挡的男女| 日韩成人伦理影院| 亚洲av电影在线观看一区二区三区| 免费人成在线观看视频色| 97超视频在线观看视频| 九九爱精品视频在线观看| 一区在线观看完整版| 男人和女人高潮做爰伦理| 成人高潮视频无遮挡免费网站| av一本久久久久| 国产淫语在线视频| 国产亚洲欧美精品永久| 国产淫语在线视频| 亚洲欧洲日产国产| 久热这里只有精品99| 97超视频在线观看视频| 国产一区有黄有色的免费视频| 美女国产视频在线观看| 涩涩av久久男人的天堂| 丰满人妻一区二区三区视频av| av线在线观看网站| 亚洲av男天堂| 国产亚洲av片在线观看秒播厂| 日韩制服骚丝袜av| 老师上课跳d突然被开到最大视频| 成人国产麻豆网| 纯流量卡能插随身wifi吗| 日本av免费视频播放| 亚洲va在线va天堂va国产| kizo精华| 亚洲国产av新网站| 国产亚洲最大av| 91精品国产九色| 久久国产精品大桥未久av | 在线观看三级黄色| 国产精品一区www在线观看| 麻豆成人av视频| 九九久久精品国产亚洲av麻豆| 熟女人妻精品中文字幕| 国产在线视频一区二区| 18禁裸乳无遮挡动漫免费视频| 99热6这里只有精品| 国产黄片美女视频| 欧美国产精品一级二级三级 | 免费黄网站久久成人精品| 国产精品一区二区性色av| 国产av国产精品国产| 观看免费一级毛片| 久久热精品热| 国产成人免费无遮挡视频| 免费人成在线观看视频色| 草草在线视频免费看| a级毛色黄片| 成人毛片60女人毛片免费| 卡戴珊不雅视频在线播放| 亚洲,欧美,日韩| 国产一区亚洲一区在线观看| 日本一二三区视频观看| 欧美zozozo另类| 国产亚洲av片在线观看秒播厂| 亚洲精品日韩在线中文字幕| 国产成人免费无遮挡视频| 国产成人一区二区在线| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 亚洲天堂av无毛| 免费在线观看成人毛片| 我的老师免费观看完整版| 一区二区三区四区激情视频| 精品一区二区三卡| 最近最新中文字幕免费大全7| av女优亚洲男人天堂| 成人亚洲欧美一区二区av| 久久精品久久精品一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 中文字幕免费在线视频6| 亚洲精品久久久久久婷婷小说| 亚洲av免费高清在线观看| 国产高清不卡午夜福利| 成人特级av手机在线观看| 一级爰片在线观看| 日韩精品有码人妻一区| 秋霞在线观看毛片| 亚洲精品aⅴ在线观看| 毛片女人毛片| 亚洲经典国产精华液单| 国产精品秋霞免费鲁丝片| 国产精品嫩草影院av在线观看| 免费看不卡的av| 最近最新中文字幕免费大全7| 国产乱来视频区| 精品熟女少妇av免费看| 精品久久久久久久久亚洲| videossex国产| 亚洲,欧美,日韩| 中文精品一卡2卡3卡4更新| 欧美3d第一页| 午夜福利在线观看免费完整高清在| 丰满迷人的少妇在线观看| 欧美97在线视频| 少妇人妻 视频| 韩国av在线不卡| 天堂中文最新版在线下载| 久久97久久精品| 国产在线男女| 一本—道久久a久久精品蜜桃钙片| 一级av片app| 亚洲精品久久久久久婷婷小说| 在线免费十八禁| 全区人妻精品视频| 亚洲成人手机| www.av在线官网国产| av在线播放精品| 亚洲国产日韩一区二区| 欧美性感艳星| 国产淫语在线视频| 爱豆传媒免费全集在线观看| 内地一区二区视频在线| 最近最新中文字幕免费大全7| 亚洲三级黄色毛片| 丰满迷人的少妇在线观看| 看免费成人av毛片| 高清欧美精品videossex| 国产女主播在线喷水免费视频网站| 少妇人妻精品综合一区二区| 观看av在线不卡| 日日啪夜夜爽| 亚洲国产欧美人成| 免费看光身美女| 国模一区二区三区四区视频| 国产女主播在线喷水免费视频网站| 秋霞在线观看毛片| 日韩精品有码人妻一区| 免费少妇av软件| 日韩强制内射视频| 国产国拍精品亚洲av在线观看| 偷拍熟女少妇极品色| 亚洲中文av在线| 精品久久久精品久久久| 天天躁夜夜躁狠狠久久av| 欧美丝袜亚洲另类| 一区二区三区精品91| 国产欧美日韩一区二区三区在线 | 天天躁夜夜躁狠狠久久av| 日本黄色片子视频| 久久久久精品久久久久真实原创| 一个人免费看片子| 麻豆国产97在线/欧美| 高清视频免费观看一区二区| 大片免费播放器 马上看| 人妻夜夜爽99麻豆av| 免费在线观看成人毛片| 一区二区av电影网| 男女无遮挡免费网站观看| 欧美一区二区亚洲| 国产av精品麻豆| a 毛片基地| 深夜a级毛片| 国产日韩欧美亚洲二区| 夫妻午夜视频| 久久鲁丝午夜福利片| 综合色丁香网| 日韩欧美一区视频在线观看 | 高清欧美精品videossex| 亚洲久久久国产精品| 日韩av免费高清视频| 丝袜喷水一区| 国产成人91sexporn| 丝瓜视频免费看黄片| 免费播放大片免费观看视频在线观看| 女的被弄到高潮叫床怎么办| 亚洲精品久久久久久婷婷小说| 国产黄色视频一区二区在线观看| 大香蕉久久网| av在线播放精品| 欧美激情极品国产一区二区三区 | 国产成人aa在线观看| 纵有疾风起免费观看全集完整版| 五月天丁香电影| 1000部很黄的大片| 久久久久网色| 两个人的视频大全免费| 热99国产精品久久久久久7| 人妻制服诱惑在线中文字幕| 超碰av人人做人人爽久久| 色5月婷婷丁香| 99九九线精品视频在线观看视频| 人人妻人人爽人人添夜夜欢视频 | 久久久久久伊人网av| 伊人久久国产一区二区| 精品一区二区免费观看| 日日啪夜夜爽| 色婷婷av一区二区三区视频| 中文字幕人妻熟人妻熟丝袜美| 丰满乱子伦码专区| www.av在线官网国产| 在线 av 中文字幕| 亚洲精品日韩av片在线观看| 久久久a久久爽久久v久久| 99九九线精品视频在线观看视频| 色网站视频免费| 99热6这里只有精品| 国产精品秋霞免费鲁丝片| 80岁老熟妇乱子伦牲交| 国产成人a∨麻豆精品| 久久久久久九九精品二区国产| 女的被弄到高潮叫床怎么办| 国产精品免费大片| 岛国毛片在线播放| 麻豆成人午夜福利视频| 观看av在线不卡| 夜夜爽夜夜爽视频| 日本av手机在线免费观看| 肉色欧美久久久久久久蜜桃| 欧美丝袜亚洲另类| 久久久久久伊人网av| 蜜桃在线观看..| 能在线免费看毛片的网站| 午夜福利在线在线| 精品亚洲成国产av| 国国产精品蜜臀av免费| 亚洲成色77777| 久久精品夜色国产| 亚洲伊人久久精品综合| 97超碰精品成人国产| av卡一久久| 在线观看av片永久免费下载| 晚上一个人看的免费电影| 爱豆传媒免费全集在线观看| 99国产精品免费福利视频| 热99国产精品久久久久久7| 亚洲天堂av无毛| 男女无遮挡免费网站观看| 高清不卡的av网站| 久久热精品热| 久久精品国产a三级三级三级| 日韩强制内射视频| 九九在线视频观看精品| 成年免费大片在线观看| 97超视频在线观看视频| 日韩,欧美,国产一区二区三区| 亚洲四区av| 久久精品国产自在天天线| 久久这里有精品视频免费| 一区二区三区免费毛片| 久久99精品国语久久久| 免费少妇av软件| 亚洲色图综合在线观看| 精品人妻偷拍中文字幕| 尤物成人国产欧美一区二区三区| 国产精品三级大全| 在线精品无人区一区二区三 | 亚洲国产欧美在线一区| 精品少妇黑人巨大在线播放| 国产淫片久久久久久久久| 亚洲欧洲国产日韩| 亚洲欧美成人精品一区二区| 中文字幕精品免费在线观看视频 | 欧美成人午夜免费资源| 在线观看国产h片| 国产有黄有色有爽视频| 国产精品伦人一区二区| 久久久国产一区二区| 久久久久久久久久久免费av| 久久精品国产自在天天线| 少妇的逼水好多| 日韩免费高清中文字幕av| 观看av在线不卡| 人妻 亚洲 视频| 亚洲精品国产av成人精品| 亚洲精品日韩av片在线观看| a级一级毛片免费在线观看| 国模一区二区三区四区视频| 深夜a级毛片| 国产欧美日韩精品一区二区| 亚洲国产最新在线播放| 自拍偷自拍亚洲精品老妇| 欧美97在线视频| 99久久综合免费| av播播在线观看一区| 成年女人在线观看亚洲视频| 国产免费一区二区三区四区乱码| 国产午夜精品久久久久久一区二区三区| 黄色欧美视频在线观看| 亚洲高清免费不卡视频| 国产黄频视频在线观看| 一级毛片aaaaaa免费看小|