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

    不同微物理參數(shù)化方案對我國北方一次大范圍暴雪天氣過程的數(shù)值模擬研究

    2022-06-01 07:15:50王淑彩平凡孟雪峰李玉鵬
    大氣科學(xué) 2022年3期
    關(guān)鍵詞:實(shí)況暴雪對流

    王淑彩 平凡 孟雪峰 李玉鵬

    1 南京信息工程大學(xué)遙感與測繪工程學(xué)院, 南京 210044

    2 中國科學(xué)院大氣物理研究所云降水物理與強(qiáng)風(fēng)暴重點(diǎn)實(shí)驗(yàn)室(LACS), 北京 100029

    3 內(nèi)蒙古自治區(qū)氣象臺, 呼和浩特 010051

    4 吉林省氣象科學(xué)研究所, 長春 130062

    1 引言

    暴雪是我國北方冬季最常見的氣象災(zāi)害之一,常伴有大風(fēng)、寒潮以及冰凍等天氣,是主要的成災(zāi)因素。大范圍、長時間的暴雪事件對交通、電力、農(nóng)業(yè)等設(shè)施均能造成嚴(yán)重影響,局地性的大暴雪、特大暴雪事件常常危及人民生命財產(chǎn)安全,屬于典型的高影響天氣(滕方達(dá)等, 2020)。我國的暴雪天氣過程主要包括內(nèi)蒙暴雪、環(huán)渤黃海區(qū)域的回流暴雪、南方暴雪、高原暴雪以及東北暴雪,因此東北地區(qū)是我國暴雪多發(fā)的關(guān)鍵區(qū)之一,局地性的大雪和暴雪發(fā)生次數(shù)呈現(xiàn)波浪式上升趨勢(王煥毅等, 2017)。

    數(shù)值模式發(fā)展的早期,對暴雪的研究大部分還是基于天氣學(xué)分析方法對暴雪環(huán)流形勢、系統(tǒng)配置等進(jìn)行分析。秦華鋒和金榮花(2008)對2007 年東北地區(qū)的一次罕見暴雪過程成因進(jìn)行了分析,表明高低空明顯的垂直切變和高層輻散與低層輻合的配置導(dǎo)致強(qiáng)上升運(yùn)動;中低層深厚的正渦度產(chǎn)生和維持是強(qiáng)暴雪形成的主要動力機(jī)制;干冷空氣侵入是產(chǎn)生暴雪的主要觸發(fā)因素;南方暖濕水汽的北上是產(chǎn)生東北暴雪的重要原因之一。孟雪峰等(2012)對內(nèi)蒙古東北地區(qū)一次漏報的暴雪天氣進(jìn)行診斷分析,表明對流層中低層溫度平流隨高度的減小有利于對流層中低層不穩(wěn)定層結(jié)的建立;地面副冷鋒與氣旋合并加強(qiáng),850 hPa 中尺度低渦強(qiáng)烈發(fā)展,加強(qiáng)了對流層低層的輻合上升運(yùn)動,觸發(fā)不穩(wěn)定能量釋放,這是強(qiáng)降雪形成的主要原因,因此邊界層“冷墊”作用對強(qiáng)降雪有一定的增幅。前人研究表明,天氣學(xué)分析的方法僅能夠抓住降雪事件中的主要影響系統(tǒng),從而理解天氣過程的演變和發(fā)展規(guī)律,但是不能精細(xì)地描述降雪的落區(qū)和強(qiáng)度。

    隨著數(shù)值模式的發(fā)展,目前大部分學(xué)者的研究都傾向于采用天氣學(xué)分析方法和微物理過程分析方法相結(jié)合,來分析降雪天氣過程。數(shù)值模式中對云微物理過程的描述主要是通過微物理參數(shù)化方案實(shí)現(xiàn)的。根據(jù)對水成物的描述不同,云微物理參數(shù)化方案可分為兩類,一類是體積水方案(bulk),另一類是分檔方案(bin)。bulk 方案是用各類水凝物粒子對應(yīng)的經(jīng)驗(yàn)函數(shù)來描述云中各類水凝物粒子的總體粒子譜分布特征,該方案計(jì)算量小,因此廣泛應(yīng)用于實(shí)際的業(yè)務(wù)和理論研究者中。但是bulk方案只適合描述云中水凝物粒子的總體分布特征,不適合描述某個尺度段的粒子變化引起的譜演變(許煥斌和段英, 1999)。而bin 方案則是根據(jù)水成物的相態(tài)、粒子大小、形狀、密度等微物理特征,將水成物分成幾十或幾百檔,給出各檔粒子的預(yù)報方程及他們之間的相互轉(zhuǎn)化關(guān)系(史月琴和樓小鳳, 2006; Khain et al., 2015)。分檔方案將粒子群分檔處理, 描述了不同檔位粒子的相互作用,與自然圖像很接近。 但是bin 方案的計(jì)算量巨大,遠(yuǎn)不能滿足實(shí)際的業(yè)務(wù)需求,因此大多只用在理論研究中。

    云微物理參數(shù)化方案的選擇對數(shù)值模擬結(jié)果有重要影響,研究不同微物理參數(shù)化方案對天氣過程的影響是國內(nèi)外的研究熱點(diǎn)之一(Wu et al., 2013;Adams-Selin et al., 2013; 尹金方等, 2014; Khain et al., 2016; Qian et al., 2018),但是目前的大量研究都是基于颮線、冰雹、超級單體風(fēng)暴等強(qiáng)天氣過程進(jìn)行數(shù)值模擬研究(Fan et al., 2015, 2017; Yin et al., 2019)。Morrison and Milbrandt(2011)分別采用MY(Morrison-Milbrandt Yau)方案對一次超級單體風(fēng)暴進(jìn)行模擬,表明方案中對于霰、雹粒子的處理是導(dǎo)致模擬結(jié)果不同的重要原因。Yin et al.(2017)利用bin 方案和bulk 方案對我國東部的一次颮線過程進(jìn)行對比研究,表明bulk 方案模擬的雷達(dá)結(jié)構(gòu)較窄但是組織性良好,bin 方案模擬的結(jié)構(gòu)較松散但是強(qiáng)度更強(qiáng),更與實(shí)況一致,但是bulk 方案模擬的降水落區(qū)和強(qiáng)度更與實(shí)況一致。前人對暴雪等弱對流天氣過程的研究較少,但也有一些研究成果,Comin et al.(2018)利用WRF 模式,以FNL 和ERA-interim 數(shù)據(jù)為邊界和初始場條件,對安第斯山脈南部的復(fù)雜地形條件下的極端降雪事件進(jìn)行不同微物理參數(shù)化的研究,表明初始場對預(yù)報的結(jié)果至關(guān)重要,F(xiàn)NL 背景場條件下WSM6 模擬的降雪量效果最好,Morrison 和WDM5 方案模擬的過強(qiáng),效果最差;ERA-interim 的雪深模擬結(jié)果比FNL 更接近實(shí)況。Ghafarian(2021)利用WRF 模式,采用四種微物理方案(WSM6、Goddard、Morrison 和Thompson)和兩種行星邊界層方案 [Yonsei University(YSU)和Mellor-Yamada-Janjic(MYJ)] 得出了八種不同的組合,對2014年發(fā)生在里海西南海岸的湖泊效應(yīng)導(dǎo)致的一次降雪過程進(jìn)行數(shù)值模擬研究,結(jié)果表明Morrison-MYJ和Goddard-MYJ 是估算降水和降雪時空變化的最佳配置,MYJ 邊界層方案比YSU 邊界層方案更能模擬出降水的強(qiáng)度和分布。國內(nèi)也對暴雪天氣過程進(jìn)行了一些研究,林文實(shí)等(2009)在利用MM5模式模擬華北地區(qū)的一次暴雪天氣過程中,比較兩個純顯式冰相云微物理參數(shù)化方案的差異,兩個方案對降雪分布和強(qiáng)度預(yù)報差別較小,但云中的微物理過程有較大的不同。崔錦等(2014)研究了不同微物理參數(shù)化方案對我國東北降水相態(tài)的影響,表明不同微物理方案對降水落區(qū)和降水強(qiáng)度預(yù)報影響不明顯,而降水相態(tài)對微物理參數(shù)化方案較敏感,主要表現(xiàn)在對雨區(qū)和雨夾雪區(qū)預(yù)報影響顯著。滕方達(dá)等(2020)利用Lin(Purdue Lin)、WSM6(WRF Single-Moment 6-class scheme)、Thompson、WDM6(WRF Double-Moment 6-class scheme)四種微物理過程參數(shù)化方案對遼寧的一次暴雪過程進(jìn)行對比研究,表明Thompson 方案模擬出更多的雪粒子和最少的霰粒子,Lin 方案霰粒子南北范圍廣、伸展高度高,WSM6 和WDM6 兩種方案模擬出較少的霰粒子。綜上所述,不同的物理參數(shù)化方案對水成物的種類、密度、數(shù)目、譜分布的設(shè)置差別較大,對云微物理過程的物理描述和處理方法上也存在較大差異,因此不同微物理參數(shù)化方案模擬的暴雪天氣過程的熱動力過程、水凝物過程、雷達(dá)回波以及地面降水等也應(yīng)存在顯著區(qū)別。但是目前的大多數(shù)微物理方案的研究都是基于暴雨、颮線、冰雹等強(qiáng)對流天氣進(jìn)行研究的,對降雪天氣過程的研究很少。而降雪天氣過程的雷達(dá)回波很弱,對流發(fā)展也偏弱,不同微物理參數(shù)化模擬的動熱力、水凝物等是否存在很大差異還需要進(jìn)行大量研究。

    本文利用WRF 模式對我國北方的一次大范圍暴雪天氣過程進(jìn)行數(shù)值模擬,對比不同微物理參數(shù)化方案 [Thompson、Morrison、WDM6、NSSL、SBM fast(bin)] 模擬的雷達(dá)回波、降水量、動熱力以及水凝物等物理量特征,從而為暴雪天氣過程的預(yù)報提供一定的理論研究基礎(chǔ)。

    2 資料和方法

    2.1 個例介紹

    2020 年4 月19 日至20 日15 時(協(xié)調(diào)世界時,下同),我國東北地區(qū)受到一次大范圍暴雪天氣過程侵襲,東北三省大部分地區(qū)以及內(nèi)蒙古東部都不同程度的受到影響。內(nèi)蒙古呼倫貝爾、黑龍江齊齊哈爾和黑河局地積雪達(dá)15~25 cm。齊齊哈爾市為當(dāng)?shù)赜袣庀笥涗浺詠? 月下旬最大,其中,內(nèi)蒙古的扎蘭屯積雪深度最大,最深達(dá)47 cm,降雪分布范圍和強(qiáng)度如圖1 所示。本次大范圍降雪天氣過程造成了交通癱瘓,鐵路、公路以及機(jī)場關(guān)閉。且發(fā)生在春耕的關(guān)鍵時期,對農(nóng)業(yè)設(shè)施也造成了很大破壞,發(fā)生在疫情防控的關(guān)鍵時期,大范圍的降溫過程,增加了公眾健康安全風(fēng)險。

    圖1 2020 年4 月19 日12 時至20 日12 時(協(xié)調(diào)世界時,下同),24 小時累計(jì)降水量實(shí)況分布(單位:mm)。散點(diǎn)圖的顏色代表降水量的量級,紅色框?yàn)楸敬卧囼?yàn)主要關(guān)注的區(qū)域Fig. 1 Distribution of the 24-h accumulative precipitation from 1200 UTC 19 April to 1200 UTC 20 April 2020 (units: mm). The colors of the scatter plot represent the magnitude of precipitation. Red boxes are areas of major concern in this experiment

    2.2 天氣形勢分析

    為了分析此次暴雪天氣過程的天氣背景條件,利用歐洲中期天氣預(yù)報中心(European Centre For Medium Range Weather Forecasts,ECMWF)的全球再分析數(shù)據(jù)(ECMWF Reanalysis v5,ERA5),分別分析了此次大范圍暴雪天氣過程發(fā)生前的500 hPa、850 hPa 以及地面的環(huán)流場、溫度場和風(fēng)場。ERA5 的空間分辨率為0.25°×0.25°,時間分辨率為1 h。從500 hPa 高度場(圖2a)上可知,亞洲中高緯地區(qū)為“一槽一脊”形勢,蒙古國中西部地區(qū)以及中國東北華北地區(qū)被高壓脊所控制。我國東北的北部有一深厚的冷渦,中心位于貝加爾湖以東和庫頁島以西附近,冷空氣沿著槽線附近進(jìn)入我國內(nèi)蒙東部從而影響整個東北地區(qū),且我國大部分東北地區(qū)以及中東部地區(qū)受到急流的影響,對流不穩(wěn)定增強(qiáng)。從850 hPa(圖2b)低空環(huán)流和比濕的分布來看,低渦中心位于內(nèi)蒙東部地區(qū),低渦前部的西南氣流和偏南氣流從南方攜帶暖濕氣流,并向東北輸送。從地面天氣形勢(圖2c)可知,我國東北正好在低壓中心附近,受到低壓中心的控制,其中心氣壓低于1001 hPa,且低壓中心具有不斷加強(qiáng)的趨勢。北方的干冷氣流與南方的飽和暖濕氣流在東北地區(qū)匯集。這種低層輻合高層輻散的天氣系統(tǒng)加上低層大量水汽的輸送,導(dǎo)致這次大范圍暴雪天氣發(fā)生和維持。

    圖2 ERA5 再分析資料2020 年4 月19 日06 時天氣形勢圖:(a)500 hPa 位勢高度(藍(lán)色等值線,單位:10 gpm)、溫度(填色,單位:K)、風(fēng)場(箭矢,單位:m s-1)和急流(綠色陰影,風(fēng)速大于20 m s-1);(b)850 hPa 位勢高度(藍(lán)色等值線,單位:10 gpm)、風(fēng)場(箭矢,單位m s-1)和比濕(填色,單位:g kg-1);(c)地面海平面平均氣壓(藍(lán)色等值線,單位:hPa)、風(fēng)場(箭矢,單位:m s-1)、溫度(色階,單位:K)和溫度露點(diǎn)差(綠色陰影,<2 :單位:K)Fig. 2 (a) Geopotential height (blue contour, units: 10 gpm), temperature (shaded, units: K), wind (vector, units: m s-1), and wind jet (green shadow,wind speed over 20 m s-1) at 500 hPa; (b) geopotential height (blue contour, units: 10 gpm), wind (vector, units: m s-1), specific humidity (shaded,units: g kg-1) at 850 hPa; (c) mean sea level pressure (blue contour, units: hPa), wind (vector, units: m s-1), temperature (shaded, units: K) and temperature dew point difference (the green shadow is less than 2, units: K) on 1200 UTC 19 April 2020 based on the ERA5 reanalysis data

    2.3 模式設(shè)計(jì)和觀測資料

    本文采用WRF4.1.5 中尺度模式對本次暴雪天氣過程進(jìn)行數(shù)值模擬,模式為單層,模擬區(qū)域如圖3,空間分辨率為3 km,格點(diǎn)數(shù)為500×530,中心為(46.5°N,124°E),垂直方向?yàn)?1 層,模式層頂為50 hPa。采用分辨率為0.25°×0.25°的ERA5 再分析資料為初始和邊界條件,模擬的初始時刻為2020 年2 月19 日00 時,模擬時長為36 個小時,每隔1 小時輸出一次結(jié)果。為了研究模式不同微物理方案對這次大范圍暴雪天氣過程的影響,選 擇 了Thompson(8)、Morrison(10)、WDM6(16)、NSSL(18)、SBM fast (bin, 30) 五種云微物理過程方案進(jìn)行數(shù)值模擬,Thompson(8)、Morrison(10)、WDM6(16)、NSSL(18)這四個方案屬于bulk 方案,除云微物理方案不同外,均采用YSU 邊界層方案,Noah 陸面過程方案,不采用積云對流參數(shù)化方案,RRTM 長波輻射方案,Dudhia 短波輻射方案,具體的參數(shù)化方案設(shè)置如表1 所示,表2 為5 個微物理參數(shù)化方案對應(yīng)的混合比和數(shù)濃度。為了驗(yàn)證本次試驗(yàn)?zāi)M的效果,采用地面氣象站數(shù)據(jù)以及雷達(dá)基數(shù)據(jù)進(jìn)行驗(yàn)證。

    表1 模式設(shè)計(jì)Table 1 Mode parameterization scheme setting

    表2 微物理參數(shù)化方案的混合比和數(shù)濃度Table 2 Mixing ratio and number concentration of the five microphysical parameterization schemes

    圖3 模擬區(qū)域Fig. 3 Model domain for the numerical simulation

    3 模擬結(jié)果驗(yàn)證

    為了分析不同的微物理參數(shù)化方案對這次暴雪天氣過程模擬的效果,本文主要對不同微物理參數(shù)化方案模擬得到的雷達(dá)組合反射率、降水水平分布分布以及雷達(dá)反射率垂直分布與實(shí)況進(jìn)行對比。

    3.1 雷達(dá)反射率

    天氣雷達(dá)是利用云霧、雨、雪等降水粒子對電磁波的散射和吸收作用,從而可探測到降水的空間分布和鉛直結(jié)構(gòu),雷達(dá)觀測的反射率因子能很好地反映大氣中的水物質(zhì)含量,所以利用不同的云微物理參數(shù)化方案模擬的雷達(dá)反射率與雷達(dá)基數(shù)據(jù)進(jìn)行比較。此次大范圍暴雪天氣過程從2020 年4 月19日02 時左右開始生成并不斷向東發(fā)展,到19 日20 時左右發(fā)展成熟,并開始慢慢消散,持續(xù)時間大概36 個小時左右。本次試驗(yàn)重點(diǎn)關(guān)注的區(qū)域是降雪量比較大的地區(qū),其位置如圖1 所示,分別選取了該區(qū)域內(nèi)對流系統(tǒng)生成、發(fā)展和成熟三個階段進(jìn)行分析。

    總體上,5 種微物理參數(shù)化方案模擬的雷達(dá)反射率在時間和空間上分布基本一致。實(shí)測的雷達(dá)最強(qiáng)回波時間是19 日20 時(圖4a3),而模擬的最強(qiáng)回波時間是19 日21 時(圖4b3、c3、d3、e3、f3),較觀測晚了約1 小時。圖4 為不同的微物理方案模擬不同時次的雷達(dá)回波與實(shí)測的雷達(dá)回波對比圖。圖4a1、a2、a3 是雷達(dá)實(shí)況,由于雷達(dá)資料的不全,只有少數(shù)的幾個自動站,因此實(shí)況的雷達(dá)數(shù)據(jù)只覆蓋部分地區(qū)。第一列為主要關(guān)注區(qū)域的對流初生時刻,第二列為對流發(fā)展時刻,第三列為對流成熟時刻。5 種微物理方案模擬的雷達(dá)回波范圍明顯比實(shí)況偏大,但是雷達(dá)回波中心強(qiáng)度和落區(qū)基本與實(shí)況一致,強(qiáng)回波中心大于45 dBZ,在內(nèi)蒙古的東側(cè)有一個發(fā)展強(qiáng)烈且范圍廣的對流系統(tǒng),從地面天氣形勢圖(圖2c)可知,該地區(qū)受強(qiáng)低壓中心控制,低壓中心在該地區(qū)停留將近兩天(圖略),這與強(qiáng)對流系統(tǒng)持續(xù)時間一致。該對流系統(tǒng)導(dǎo)致扎蘭屯、齊齊哈爾、白城等區(qū)域大范圍、高強(qiáng)度、持續(xù)時間長的降雪過程。從整個過程看,bin方案模擬的強(qiáng)回波中心與實(shí)況相比偏弱(圖4f1、f2、f3),而Morrison 模擬的雷達(dá)回波中心偏強(qiáng)(圖4c1、c2、c3)。與bin 方案相比,bulk 方案模擬的回波結(jié)構(gòu)較松散,單體相對孤立,組織化程度較差。從整個時間周期上看(圖略),雷達(dá)回波在20 日12 時左右基本消散,而bulk 方案模擬的雷達(dá)回波中心還有30 dBZ左右,bin 方案的回波強(qiáng)度基本消散,因此bin 方案的發(fā)展周期比bulk 方案的短,更與實(shí)況一致。

    不同的微物理參數(shù)化方案模擬的雷達(dá)回波的水平分布差異較小,因此也對比了觀測和模擬的雷達(dá)反射率垂直結(jié)構(gòu)分布,各剖面的位置如圖4 紫色框?qū)蔷€位置。4 月19 日14 時(圖5a1),黑龍江省、內(nèi)蒙古與吉林的交界地帶有幾個初生的對流單體,回波中心強(qiáng)度在35 dBZ以下,懸垂在850~600 hPa 之間。5 種微物理方案模擬的雷達(dá)回波剖面(圖5b1、c1、d1、e1、f1)基本都與實(shí)況一致。但是NSSL 方案模擬的回波底偏低,WDM6 方案模擬的回波頂偏低。隨著對流系統(tǒng)的發(fā)展(圖5b1),對流單體的水平寬度以及垂直高度不斷擴(kuò)大,多對流中心發(fā)展并開始合并,強(qiáng)回波中心最高可達(dá)到45 dBZ,5 種微物理方案(圖5b2、c2、d2、e2、f2)均模擬出了多單體逐步發(fā)展合并的特征,且其形態(tài)基本與實(shí)況一致,但是WDM6(圖5d2)和NSSL(圖5e2)方案模擬的強(qiáng)回波中心與實(shí)況相比偏弱,在21 時(圖5a3),對流系統(tǒng)發(fā)展達(dá)到最強(qiáng),中心回波強(qiáng)度達(dá)到50 dBZ,Thompson(圖5b3)和Morrison(圖5c3)均模擬出50 dBZ的強(qiáng)回波中心。但是Morrison 模擬的回波形態(tài)比較寬廣,更與實(shí)況一致。從整個過程上看,Morrison模擬的回波強(qiáng)度、范圍和形態(tài)更與實(shí)況一致,bin與bulk 方案相比,模擬的單體更多,這可能這主要是因?yàn)閎in 方案將粒子群分檔處理,沒有捆綁不同粒子類型運(yùn)動,更能描述出不同檔位粒子相互作用。

    圖4 2022 年4 月19 日(a1)15:00、(a2)17:00、(a3)20:00 實(shí)測和(b1、c1、d1、e1、f1)14:00、(b2、c2、d2、e2、f2)16:00、(b3、c3、d3、e3、f3)21:00(b1-b3)Thompson 方案、(c1-c3)Morrison 方案、(d1-d3)WDM6 方案、(e1-e3)NSSL 方案及(f1-f3)bin 方案模擬的雷達(dá)組合反射率(單位:dBZ)。紅色對角線和方框分別表示下文剖面位置和區(qū)域平均的范圍Fig. 4 Observed (a1) 1500 UTC, (a2) 1700 UTC, (a3) 2000 UTC and simulated radar composite reflectivity (units: dBZ); (b1) 1400 UTC 19 April 2020, (b2) 1600 UTC 19 April 2020, (b3) 2100 UTC 19 April 2020) correspond to Thompson (the time of other schemes is the same below), (c1, c2,c3) is Morrison, (d1, d2, d3) is WDM6, (e1, e2, e3) is NSSL, and (f1, f2, f3) is bin. Red diagonal lines and red boxes indicate the slice lines of crosssections and the averaging area in the following analysis, respectively

    圖5 實(shí)測和微物理參數(shù)化方案模擬的雷達(dá)反射率垂直剖面(單位:dBZ)和風(fēng)場(箭矢:單位:m s-1)(其垂直速度風(fēng)矢擴(kuò)大了50 倍,下同),時間和剖面的位置如圖4 所示Fig. 5 Observed and simulated vertical cross-section of the radar composite reflectivity (units: dBZ) and winds (vector, units: m s-1; vertical vector speed is magnified by 50, the same below). The analysis time and locations of cross-sections are denoted in Fig. 4

    3.2 地面累計(jì)降水分布

    本文在對比雷達(dá)回波的同時,對降水的分布情況也進(jìn)行了對比。由于不同微物理參數(shù)化方案對水物質(zhì)的組成和分布描述不同,會導(dǎo)致模式對降水的模擬產(chǎn)生差異,5 種微物理參數(shù)化方案模擬得到的6 小時累計(jì)降水量與站點(diǎn)實(shí)況降水觀測的對比如圖6 所示,其中圖6a1、a2、a3 是地面自動站降水分布情況。在初期降水帶主要呈東北—西南走向(圖6a1),由地面天氣形勢圖(圖2c)可知,我國內(nèi)蒙古、黑龍江、吉林三省交界地帶受低壓中心控制,強(qiáng)低壓中心在此停滯不前,南部受強(qiáng)勁的西南風(fēng)影響,在西南風(fēng)不斷加強(qiáng)和推進(jìn)作用下,降水帶逐漸呈南北走向(圖6a2),最后呈西北—東南走向(圖6a3),最后移出我國區(qū)域,持續(xù)時間30 個小時左右。從圖6 可知5 種微物理參數(shù)化方案模擬的降水量分布、雨帶走勢、強(qiáng)度基本與實(shí)況一致,基本模擬出了強(qiáng)降水中心和發(fā)展的走勢,但是模擬的強(qiáng)度中心(>12.8 mm)和范圍明顯比實(shí)況要大。與此同時,也對比5 個方案模擬的24 小時累計(jì)降水量,從降水量分布情況(圖7)來看,5 個微物理參數(shù)化方案模擬的24 小時累計(jì)降水量的強(qiáng)度和落區(qū)與實(shí)況基本一致,WDM6 方案模擬的強(qiáng)降水(>25.6 mm)中心的強(qiáng)度和范圍明顯比實(shí)況偏大。

    圖6 觀測和模擬的6 小時累計(jì)降水量(單位:mm):(a1、b1、c1、d1、e1、f1)2020 年4 月19 日10:00~16:00;(a2、b2、c2、d2、e2、f2)2020 年4 月19 日20:00 至20 日02:00;(a3、b3、c3、d3、e3、f3)2020 年4 月20 日06:00~12:00Fig. 6 Observed and simulated of the 6-h accumulated rainfall (units: mm) during (a1, b1, c1, d1, e1, f1) 1000-1600 UTC 19, (a2, b2, c2, d2, e2, f2)2000 UTC 19-0200 UTC 20, and (a3, b3, c3, d3, e3, f3) 0600 UTC-1200 UTC 20

    為了定量評價不同微物理參數(shù)化方案對這次暴雪天氣過程的模擬結(jié)果,選取了東北258 個地面自動站的1 小時累計(jì)降水均值資料對模式的模擬能力進(jìn)行客觀定量評價。其結(jié)果如圖8 和表3 所示,4月19 日22 時左右,降水達(dá)到峰值,但是模擬的降水峰值與實(shí)況相比提前了1 小時左右。19 日23 時以前,模式模擬的降水量級與實(shí)況一致,但是23時以后出現(xiàn)了較大偏差,模擬降水的明顯比實(shí)況偏高。從表3 的相關(guān)系數(shù)和均方根誤差結(jié)果來看,相關(guān)系數(shù)最大的是Morrison 方案,達(dá)到0.9591,均方根誤差最小的是Thompson 方案,為0.0640。可以認(rèn)為本次過程中Thompson 和Morrison 方案整體優(yōu)于其他方案。從24 小時累計(jì)降水(圖7)主觀檢驗(yàn)上看,Morrison 方案的強(qiáng)降水中心(>25.6 mm)范圍比Thompson 方案的小,更與實(shí)況一致,因此綜合客觀定量評價和主觀檢驗(yàn),Morrison 方案對此次降雪過程的模擬更好。

    圖7 2020 年4 月19 日06 時至20 日06 時實(shí)況與不同微物理參數(shù)化模擬的24 小時累計(jì)降水量(單位:mm)Fig. 7 Observed and simulated of the 24 h accumulated rainfall during 0600 UTC 19 April 2020-0600 UTC 20 April 2020 for simulation (units:mm).

    表3 1 小時累計(jì)降水量統(tǒng)計(jì)評估Table 3 One-hour accumulated rainfall statistical assessment

    4 微物理參數(shù)化方案對云內(nèi)熱動力的影響

    云中的動力和熱力場能為云中的微物理過程提供背景流場,并影響云和降水粒子的數(shù)濃度、初始尺度譜以及相應(yīng)的物理性質(zhì),規(guī)范了微物理過程進(jìn)行的速率和持續(xù)的時間,以及最終降水物的尺度。反過來,云中的微物理過程對熱動力有很好的反饋?zhàn)饔?。通過凝結(jié)、蒸發(fā)、凝華等過程改變水汽分布,相變潛熱的釋放和吸收,為云動力過程提供了重要的熱源和熱匯。因此對云微物理結(jié)構(gòu)進(jìn)行分析需要對云內(nèi)熱動力結(jié)構(gòu)進(jìn)行分析。圖9 是5 個微物理參數(shù)化方案模擬的區(qū)域平均垂直速度場和風(fēng)場結(jié)構(gòu)的剖面。充足的水汽輸送和旺盛的上升運(yùn)動是降水產(chǎn)生的基本條件,上升運(yùn)動可使水汽從未飽和狀態(tài)達(dá)到飽和狀態(tài),從而有利于降水形成。從圖中的風(fēng)場結(jié)構(gòu)看,該地區(qū)主要吹偏南風(fēng),可以從南方帶來暖濕的水汽,有利于降水的形成。通常降雪過程的對流較弱,因此總體上模擬的垂直速度量級都在2 m s-1以下。5 個微物理參數(shù)化方案模擬垂直速度場都以上升運(yùn)動為主,300 hPa 以上垂直速度基本接近0。在整個對流過程中,對流發(fā)展階段(圖9a2、b2、c2、d2、e2)的上升運(yùn)動達(dá)到最強(qiáng),到了成熟階段(圖9a3、b3、c3、d3、e3),強(qiáng)對流已經(jīng)變?nèi)醪⒅鸩较?,但是Morrison 方案模擬的帶狀強(qiáng)上升運(yùn)動中心仍然比較強(qiáng),說明Morrison 方案模擬的強(qiáng)上升運(yùn)動持續(xù)時間比其他方案更長,這與Morrison 方案模擬的降水量(圖8)在19 日16~21 時這段時間內(nèi)比其他方案偏高是相對應(yīng)的。

    圖8 東北258 個地面自動站與不同微物理參數(shù)化模擬的平均1 小時累計(jì)降水量時間演變(單位:mm)Fig. 8 Temporal evolution of the average 1-hour accumulated rainfall at 258 automatic ground stations in Northeast China simulated by different microphysical parameterizations (units: mm)

    圖9 五種方案模擬的垂直速度 (填色,單位:m s-1)、風(fēng)場(箭矢,單位:m s-1),時間和剖面位置同圖4Fig. 9 Vertical cross-sections of the vertical velocity (shaded, units: m s-1) and winds (vectors, units: m s-1) simulated by different microphysical parameterizations. The analysis time and locations of cross-sections are denoted in Fig. 4

    圖10 為渦度場的垂直剖面,從渦度場可知,隨著對流系統(tǒng)的發(fā)展,正渦度不斷加強(qiáng),范圍也不斷向高層發(fā)展,說明輻合上升運(yùn)動不斷加強(qiáng)??傮w上低層渦度為正,為輻合區(qū),中層渦度為負(fù),為強(qiáng)的輻散區(qū),這種低層輻合,高層輻散的環(huán)境配置有利于強(qiáng)對流系統(tǒng)的形成發(fā)展,有利于降水形成(宗志平和劉文明, 2004)。劉寧微(2006)分析了2003 年3 月遼寧的暴雪過程,指出低空輻合和高空輻散,導(dǎo)致上升運(yùn)動的加強(qiáng)以及低層正渦度的產(chǎn)生和維持而產(chǎn)生的垂直方向上水汽凝結(jié)是該暴雪過程的形成機(jī)制,5 個微物理參數(shù)化方案模擬的渦度都體現(xiàn)了這個特征。Thompson 方案(圖10a3)和NSSL 方案(圖10d3)模擬的低層正渦度(>0.0008)延伸高度更高,說明輻合上升運(yùn)動更強(qiáng),這與Thompson 方案和NSSL 方案的雷達(dá)回波頂較高是一致的。在初始階段(圖10a1、b1、c1、d1、e1),正渦度中心在850 hPa 附近發(fā)展,隨著對流系統(tǒng)的發(fā)展(圖10a2、b2、c2、d2、e2),正渦度中心開始向高層發(fā)展,達(dá)到500 hPa 附近,正渦度中心伴隨著有負(fù)渦度中心,到了成熟階段(圖10a3、b3、c3、d3、e3),正渦度中心已經(jīng)逐漸消散。在對流發(fā)展階段,5 個微物理參數(shù)化方案在600 hPa以下均有波列分布,但是bin 方案的波列結(jié)構(gòu)更細(xì)小、數(shù)量更多、更明顯,這主要是因?yàn)閎in 方案將粒子群分檔處理,沒有捆綁不同粒子類型運(yùn)動,更能細(xì)致描述出不同粒子的下沉拖曳作用。

    圖10 微物理參數(shù)化方案模擬的渦度場(填色,單位:s-1)和風(fēng)場(箭矢,單位:m s-1)垂直剖面,時間和剖面位置同圖4Fig. 10 Vertical cross-sections of the vorticity field (shaded, units: s-1) and wind (vector, units: m s-1) simulation by different microphysical parameterizations. The analysis time and the locations of cross-sections are denoted in Fig. 4

    圖11 是擾動位溫和總水凝物混合比的垂直分布圖,擾動位溫是三個對流發(fā)展階段(初生、發(fā)展、成熟)分別與模式初始時刻(19 日00 時)的差值。從擾動位溫來看,5 種方案都模擬出了負(fù)—正—負(fù)的位溫結(jié)構(gòu),即在低層(700 hPa 以下)擾動位溫為負(fù)值,在中層(700~500 hPa)擾動位溫為正值,在中高層(500~300 hPa)擾動位溫為負(fù)值,5 種微物理參數(shù)化方案模擬的擾動位溫差異不明顯。隨著對流系統(tǒng)的發(fā)展,擾動位溫不斷加強(qiáng),Thompson方案模擬的成熟階段(圖11a3)的擾動位溫在低層發(fā)展最強(qiáng)烈,說明該方案模擬的冷池強(qiáng)度最強(qiáng),更有利于對流系統(tǒng)發(fā)展。水凝物主要分布在500 hPa以下的中低層區(qū),不同的微物理參數(shù)化方案模擬的總水凝物分布與垂直速度是對應(yīng)的,垂直上升運(yùn)動中心對應(yīng)著水凝物的混合比最大值中心,強(qiáng)上升運(yùn)動的延伸高度和寬度對應(yīng)著水凝物的延伸高度和寬度。

    圖11 不同方案模擬的擾動位溫(填色,單位:K)和總水凝物混合比(等值線,單位:g kg-1)的垂直剖面,紅色實(shí)線為0°C 等溫線,時間和剖面位置如圖4Fig. 11 Vertical cross-sections of the potential temperature perturbation (shaded, units: K) and the total mixing ratio (contours, units: g kg-1). Red solid lines represent the 0°C isotherms. The analysis time and locations of cross-sections are denoted in Fig. 4

    為了更細(xì)致地描述整個對流過程的垂直速度、渦度以及擾動位溫特征,進(jìn)行了最大最小垂直速度、渦度以及平均擾動位溫隨時間演變特征統(tǒng)計(jì)(圖12),不同微物理參數(shù)化方案所選取的區(qū)域位置如圖4(b3、c3、d3、e3、f3)所示。從圖中可知,各方案的垂直速度、渦度和擾動位溫走勢基本一致,渦度的走勢與垂直速度的走勢基本一致,從整體上看,19 日10 時左右,WDM6 方案模擬的最大垂直速度(圖12a1)和最大渦度(圖12b1)與其它方案相比明顯偏高,到了11 ~20 時左右,其最大垂直速度又明顯比其他方案偏小,這主要是因?yàn)樵?1~20 時這段時間,WDM6 方案模擬的降水(圖8)明顯比其他方案高,比較大的降水引起粒子的下沉拖曳作用,使得其垂直上升運(yùn)動明顯降低。本次暴雪過程主要發(fā)生在19 日16 時至20 日03 時,并不是發(fā)生在對流運(yùn)動最強(qiáng)烈的時段,而是發(fā)生在垂直速度極值和渦度極值都趨于變小的時段,說明相對靜穩(wěn)的環(huán)境,可能更有利于降雪發(fā)生。從平均擾動位溫分布走勢圖來看(圖12c),從19 日00 時開始,擾動位溫不斷增加,到19 日06 時,達(dá)到最大值,開始維持在一個穩(wěn)定的數(shù)值,到了19 日16 時,擾動位溫開始逐步下降,到20日03 時以后又開始回升,擾動位溫下降的時間段對應(yīng)著強(qiáng)降雪時段(圖8),所以氣溫的冷卻作用可能是此次暴雪天氣過程的形成的機(jī)制之一。從整個過程來看,bin 方案和Morrison 方案的擾動位溫最大,Thompson 方案的擾動位溫最小。

    圖12 不同微物理參數(shù)化方案模擬的(a1、a2)最大最小垂直速度、(b1、b2)渦度和(c)平均擾動位溫時間演變Fig. 12 Temporal evolution of the (a1, a2) max and min vertical velocity, (b1, b2) vorticity, and (c) potential temperature (PT) perturbation simulated by different microphysical parameterizations.

    5 微物理參數(shù)化方案對水凝物的影響

    5.1 水凝物的垂直廓線

    云微物理過程中,水成物的分布在很大程度上反映了微物理過程中水汽與降水粒子之間相互轉(zhuǎn)換的熱力動力過程(許廣等, 2017),對云微物理過程特別是水成物的分析有助于了解降水過程的機(jī)制。從圖13 可看出,5 種微物理參數(shù)化方案模擬的水凝物含量和分布存在較大差異,從整體的水凝物混合比可知,所有的微物理參數(shù)化方案中,占比最大的都是雪粒子,這與本次個例是暴雪天氣過程相對應(yīng)。冰晶主要集中在8 km 附近,3 km 為雪粒子最大層,2 km 為霰粒子最大層, 1.5 km 為云水最大層,雨水主要分布在2 km 以下的暖區(qū)。對于雪、霰、云水以及雨水的模擬,各方案模擬的廓線除了在量級上有區(qū)別外,其分布高度基本一致,而對冰晶的模擬,各方案模擬的廓線不管在強(qiáng)度和分布高度上都存在很大的差異,說明不同的微物理參數(shù)化方案對冰晶粒子的描述存在很大差異。

    圖13 區(qū)域平均的水凝物混合比垂直廓線(單位:g kg-1):(a1、a2、a3)冰晶;(b1、b2、b3)雪;(c1、c2、c3)霰;(d1、d2、d3)云水;(e1、e2、e3)雨水,時間和位置同圖4Fig. 13 Vertical profiles of the regionally average mixing ratio (units: g kg-1) of (a1, a2, a3) ice, (b1, b2, b3) snow, (c1, c2, c3) graupel, (d1, d2, d3)cloud water, and (e1, e2, e3) rainwater. The analysis time and the locations of cross-sections are denoted in Fig. 4

    通常雷達(dá)反射率與大氣中的雨、雪、霰等水凝物有關(guān),其中雨水粒子的反射率最強(qiáng),雪粒子的反射率最弱。5 種微物理參數(shù)化方案均模擬出了大氣低層的雨水粒子(圖13e1、e2、e3),但是WDM6方案模擬的雨水粒子明顯比其他方案的高,這是造成WDM6 方案模擬的反射率剖面在低層(850 hPa以下)比其他方案高且深厚的原因。NSSL 方案模擬的雨水粒子伸展高度最高,造成NSSL 方案模擬的雷達(dá)回波伸展高度也高。從雪粒子的整個時空分布來看,總體上bin 方案模擬的雪粒子量級在對流初生階段和發(fā)展階段比其他方案偏高,但是在成熟階段與其它方案基本一致。Thompson 方案模擬的雪粒子量級在成熟階段,其低層的雪粒子明顯比其他方偏高,說明Thompson 方案模擬的地面降雪量最大。而霰粒子的空間分布各方案模擬的正好與雪粒子相反,bin 和Morrison 方案模擬的雪粒子偏高,而模擬的霰粒子含量最低,WDM6 方案模擬的雪粒子偏低,而模擬的霰粒子含量偏高。云水粒子主要集中在2 km 高度附近,在對流初生階段(圖13d1)和對流發(fā)展階段(圖13d2),各方案模擬的云水粒子含量和高度出現(xiàn)較大偏差,但到了成熟階段(圖13d3),各方案模擬的云水粒子不管是在量級還是在分布高度上基本一致。從模擬的冰晶粒子分布來看,各方案存在很大的差異,WDM6 方案模擬的冰晶粒子含量最高,且其延伸范圍廣,從低層至高層9 km 附近均有冰晶粒子存在。Morrison方案和NSSL 方案模擬的冰晶粒子含量也較高,但是兩個方案模擬的冰晶粒子分布高度存在很大區(qū)別,Morrison 方案主要集中在8 km 的高層,而NSSL方案主要集中在3 km 的中低層,bin 和Thompson方案模擬的冰晶粒子含量最少。從整個時間演變上看,在對流系統(tǒng)初生和發(fā)展階段,總體上各方案模擬的水凝物在量級和分布高度上都存在較大的差距,到了成熟階段,各方案模擬的水凝物趨于一致。

    圖14 是數(shù)濃度的垂直廓線,由于WDM6 方案只預(yù)報出了雨水和云水?dāng)?shù)濃度,Thompson 方案只預(yù)報了冰晶和雨水?dāng)?shù)濃度,NSSL 方案沒有預(yù)報出云粒子,因此我們只給出了各個方案預(yù)報出的水凝物數(shù)濃度。從模擬的水凝物數(shù)濃度來看,各方案存在很大的差異。Morrison 模擬的冰晶、雪粒子以及霰粒子等固態(tài)水凝物的數(shù)濃度都很大,基本都大于其他方案,從云水的數(shù)濃度分布來看,WDM6 方案遠(yuǎn)大于bin 方案,NSSL 方案和WDM6 方案模擬的雨水遠(yuǎn)大于其他方案,bin 方案模擬的水凝物數(shù)濃度與其它方案相比都很弱。

    圖14 區(qū)域平均的水凝物數(shù)濃度垂直廓線(單位:105 kg-1):(a1、a2、a3)冰晶;(b1、b2、b3)雪;(c1、c2、c3)霰;(d1、d2、d3)云水;(e1、e2、e3)雨水。時間和位置如圖4Fig. 14 Vertical profiles of the regionally average number concentration of (a1, a2, a3) ice, (b1, b2, b3) snow, (c1, c2, c3) graupel, (d1, d2, d3) cloud water, and (e1, e2, e3) rainwater (units: 105 kg-1). The analysis time and locations of cross-sections are denoted in Fig. 4

    5.2 水凝物的時間—高度剖面

    不同的微物理參數(shù)化方案對水凝物的模擬在時間層上也存在較大差異(圖15),對于云粒子,各方案模擬的強(qiáng)度基本一致,Morrison 方案模擬的云頂高度比其他方案偏低。對于冰晶粒子和雪粒子,其分布范圍較廣,出現(xiàn)多個極值中心,從時間軸上看,NSSL 方案模擬的冰晶粒子出現(xiàn)時間比別的方案要晚了1 個小時,Thompson 方案在4 月19 日12~13 時沒有模擬出冰晶。Thompson 方案模擬的雪粒子在19 日18 時已經(jīng)開始延伸到地面,而其他方案模擬的雪粒子在19 時左右才開始延伸到地面,說明Thompson 方案模擬地面最早降雪時間比比其他方案模擬早了1 個小時。對于雨水粒子的模擬,雨水主要集中在上升氣流前側(cè)的低層,5 個方案模擬的雨水持續(xù)時間基本一致,開始降雨和結(jié)束時間也基本一致,但是bin 方案模擬的最強(qiáng)降雨時間比其他方案晚了1 小時左右。對于霰粒子的模擬,除了在強(qiáng)度上存在較大差異外,在時間周期上也存在較大差異,bulk 方案模擬的霰粒子在19 日07 時開始有霰粒子,而bin 方案在12 時才有霰粒子,晚了將近5 個小時;19 日23 時,霰粒子基本消失,但是WDM6 方案模擬的霰粒子在20 日02 時左右才全部小時,晚了將近3 小時,說明WDM6 方案模擬的霰粒子在強(qiáng)度上大于其它方案,時間周期上也比其它方案長。bin 方案和Morrison 方案模擬的霰粒子在強(qiáng)度上比其他方案弱,持續(xù)時間也比其他方案短。

    圖15 不同微物理參數(shù)化方案模擬的區(qū)域平均水凝物混合比時間—高度剖面分布:(a1、b1、c1、d1、e1)云水(Qc,填色)和冰晶(Qi,等值線);(a2、b2、c2、d2、e2)雨水(Qr,填色)和雪粒子(Qs,等值線);(a3、b3、c3、d3、e3)霰粒子(Qg,等值線),(單位:g kg-1,紅色實(shí)線為0°C 線)Fig. 15 Time-height cross-sections of the (a1, b1, c1, d1, e1) regionally average mixing ratio of cloud water (Qc, shaded) and cloud ice (Qi, contours),(a2, b2, c2, d2, e2) rainwater (Qr, shaded) and snow (Qs, contours), (a3, b3, c3, d3, e3) graupel (Qg, contours) of different microphysical parameterization schemes (units: g kg-1, red solid lines represent the 0°C isotherms)

    6 結(jié)論

    通過對2020 年4 月19~20 日發(fā)生在我國東北地區(qū)的一次大范圍暴雪天氣過程進(jìn)行不同微物理參數(shù)化方案(Thompson、Morrison、WDM6、NSSL、bin)的數(shù)值模擬,探究不同微物理過程對暴雪天氣過程的影響。比較了5 種微物理參數(shù)化方案模擬的雷達(dá)反射率、地面降水量,動熱力過程以及微物理場特征。進(jìn)而了解不同微物理參數(shù)化方案對我國東北地區(qū)的適應(yīng)性問題,以期為改進(jìn)微物理參數(shù)化方案提供一定的參考。本文得出的結(jié)論如下:

    (1)不同微物理參數(shù)化方案都較好的模擬出了這次暴雪天氣過程的形態(tài)和特征,但是Morrison方案模擬的效果更好。總體上bulk 方案模擬的雷達(dá)回波結(jié)構(gòu)較松散,單體相對孤立,組織化程度較差,bin 方案模擬的雷達(dá)回波中心與實(shí)況相比明顯偏弱。從模擬的雷達(dá)回波剖面上看,Morrison 模擬的回波強(qiáng)度、范圍和形態(tài)更與實(shí)況一致。綜合主觀檢驗(yàn)和相關(guān)系數(shù)、均方根誤差統(tǒng)計(jì)分析,Morrison方案模擬的降水量較其他參數(shù)化方案更好。因此,Morrison 方案更好的模擬出了本次降雪過程。

    (2)盡管不同方案模擬的降水量、雷達(dá)回波等差異不明顯,但是模擬的云內(nèi)動力及微物理結(jié)構(gòu)存在明顯差別,主要表現(xiàn)在不同云微物理方案模擬的垂直速度和水凝物的時空分布上。bin 方案在600 hPa 高度以下存在明顯的渦度波列,這主要是因?yàn)閎in 方案將粒子群分檔處理,沒有捆綁不同粒子類型運(yùn)動,更能細(xì)致描述出不同粒子的下沉拖曳作用,更與實(shí)況一致。在強(qiáng)降雪發(fā)生前,WDM6方案模擬的最大垂直速度明顯比其它方案強(qiáng),但到了強(qiáng)降雪階段,其最大垂直速度又明顯比其它方案偏小,這主要是因?yàn)閃DM6 方案模擬的降水明顯比其它方案高,比較大的降水引起粒子的下沉拖曳作用,使得其垂直上升運(yùn)動明顯降低。

    (3)不同方案模擬的雪、霰、云水以及雨水粒子較為接近,而對冰晶的模擬不管在量級還是在分布范圍上都存在很大的差異。WDM6 和NSSL方案,冰晶分布較其他方案更接近0°C 融化層及以下,甚至達(dá)至地面。從整個時間演變上看,在對流系統(tǒng)初生和發(fā)展階段,總體上各方案模擬的水凝物混合比在量級和分布高度上都存在較大的差距,到了成熟階段,各方案模擬的水凝物混合比趨于一致,這與垂直速度和渦度的演變特征是一致的。

    (4)bin 方案和bulk 方案的模擬差異主要體現(xiàn)在霰粒子和冰晶的模擬,其模擬的量級遠(yuǎn)小于bulk 方案,且bin 方案模擬的水凝物數(shù)濃度也遠(yuǎn)小于bulk 方案。從水凝物的時間演變特征來看,bin方案模擬的最強(qiáng)降雨時間比其它方案晚了將近1 小時,最早降霰時間也比bulk 方案晚了將近5 個小時,總體上bin 方案模擬的霰粒子在強(qiáng)度上比其它方案弱,持續(xù)時間也比其它方案短。理論上,bin方案對粒子變化的描述更與實(shí)況一致,與自然圖像更接近,但是在這次暴雪天氣過程中,bin 方案模擬的雷達(dá)回波和降水量不如bulk 方案,其優(yōu)勢并沒有體現(xiàn)出來,且在強(qiáng)度、空間分布上都不如bulk方案。

    本次研究是基于不同微物理方案對2020 年4月的一次大范圍暴雪天氣過程進(jìn)行的數(shù)值模擬結(jié)果,其研究結(jié)論還需要更多的暴雪天氣個例進(jìn)行驗(yàn)證和支撐,特別是局地性的暴雪天氣過程是否會得出同樣的結(jié)論還需要進(jìn)行進(jìn)一步的驗(yàn)證。由于本次研究沒有雙偏振雷達(dá)實(shí)際觀測資料,因此不能對云內(nèi)細(xì)致的動熱力和水凝物與實(shí)況進(jìn)行比較。未來的工作將基于偏振雷達(dá)的高精度觀測資料進(jìn)行細(xì)致地對比分析不同微物理方案模擬的云水粒子的時空分布及三維結(jié)構(gòu),探討云微物理影響降雪的可能機(jī)制。

    猜你喜歡
    實(shí)況暴雪對流
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    鄉(xiāng)村小學(xué)的愿望與現(xiàn)實(shí)——宜君鄉(xiāng)村教育實(shí)況
    天舟一號貨運(yùn)飛船發(fā)射實(shí)況掠影
    可愛潮咖們的獨(dú)門彩妝實(shí)況直播
    西藏西南部暴雪天氣分析
    西藏科技(2016年9期)2016-09-26 12:21:45
    基于ANSYS的自然對流換熱系數(shù)計(jì)算方法研究
    日喀則“118”大到暴雪天氣診斷分析
    西藏科技(2015年3期)2015-09-26 12:11:08
    二元驅(qū)油水界面Marangoni對流啟動殘余油機(jī)理
    暴雪之后街頭雪人秀
    “突降暴雪快出動,徹夜清除保暢通”等十二則
    自拍偷自拍亚洲精品老妇| 久久午夜福利片| 久久中文看片网| a级毛片a级免费在线| 久久久欧美国产精品| 特大巨黑吊av在线直播| 少妇猛男粗大的猛烈进出视频 | 啦啦啦韩国在线观看视频| 日韩 亚洲 欧美在线| 日本免费a在线| 国产伦精品一区二区三区视频9| 精品不卡国产一区二区三区| 亚洲av成人精品一区久久| 精品久久久久久久久久久久久| 日韩欧美精品免费久久| 国产高清有码在线观看视频| 美女国产视频在线观看| 一区二区三区高清视频在线| 熟妇人妻久久中文字幕3abv| 亚洲av电影不卡..在线观看| 久99久视频精品免费| 国产一区二区在线av高清观看| 精品久久国产蜜桃| 免费观看人在逋| 97人妻精品一区二区三区麻豆| 一个人看的www免费观看视频| 成人欧美大片| 日本与韩国留学比较| 国产亚洲精品久久久久久毛片| 伦理电影大哥的女人| 国产在线男女| 久久久久久久久久久丰满| 国产精品一区www在线观看| 成人二区视频| 如何舔出高潮| 人妻久久中文字幕网| 国产69精品久久久久777片| 麻豆乱淫一区二区| av在线观看视频网站免费| 亚洲国产精品合色在线| 午夜精品国产一区二区电影 | 国产免费男女视频| 久久久久久久亚洲中文字幕| 少妇熟女欧美另类| 人妻少妇偷人精品九色| 久久精品国产亚洲网站| 麻豆乱淫一区二区| 午夜福利成人在线免费观看| 亚洲不卡免费看| 精品欧美国产一区二区三| 久久久久网色| 国产人妻一区二区三区在| 免费观看的影片在线观看| 在现免费观看毛片| 久久这里只有精品中国| 在线观看一区二区三区| 色噜噜av男人的天堂激情| 国产三级中文精品| 亚洲人成网站在线播放欧美日韩| 一个人看的www免费观看视频| 啦啦啦观看免费观看视频高清| 欧美另类亚洲清纯唯美| 国产老妇伦熟女老妇高清| 91在线精品国自产拍蜜月| 97超视频在线观看视频| 97人妻精品一区二区三区麻豆| 插逼视频在线观看| 尾随美女入室| 国产 一区 欧美 日韩| 国产亚洲欧美98| 午夜免费激情av| 一个人免费在线观看电影| 亚洲一级一片aⅴ在线观看| 男女啪啪激烈高潮av片| 黄色配什么色好看| 欧美日本视频| 赤兔流量卡办理| 精品久久久久久久人妻蜜臀av| 国产高清三级在线| 午夜久久久久精精品| 狂野欧美白嫩少妇大欣赏| 国产精华一区二区三区| 午夜亚洲福利在线播放| 中国国产av一级| 国产精品人妻久久久影院| 一区福利在线观看| 五月伊人婷婷丁香| 最好的美女福利视频网| 99视频精品全部免费 在线| 人人妻人人看人人澡| 精品久久久噜噜| 直男gayav资源| a级毛色黄片| 国产老妇伦熟女老妇高清| 久久久久九九精品影院| 欧美一区二区国产精品久久精品| 久久久精品欧美日韩精品| 美女黄网站色视频| avwww免费| 免费av毛片视频| 亚洲18禁久久av| 91久久精品国产一区二区三区| 校园人妻丝袜中文字幕| 最近手机中文字幕大全| 国产美女午夜福利| 成人亚洲精品av一区二区| av黄色大香蕉| 国产伦精品一区二区三区四那| 小蜜桃在线观看免费完整版高清| 国产毛片a区久久久久| 色哟哟哟哟哟哟| 亚洲乱码一区二区免费版| 国产日本99.免费观看| 国产精品麻豆人妻色哟哟久久 | 中文字幕久久专区| 直男gayav资源| 女的被弄到高潮叫床怎么办| 一进一出抽搐动态| 99riav亚洲国产免费| 亚洲婷婷狠狠爱综合网| 日韩 亚洲 欧美在线| 久久精品国产亚洲网站| 国产精品蜜桃在线观看 | 级片在线观看| 国产在线精品亚洲第一网站| 天堂av国产一区二区熟女人妻| 男女那种视频在线观看| 老女人水多毛片| 久久久a久久爽久久v久久| 精品久久久噜噜| 亚洲精品乱码久久久v下载方式| 久久精品国产鲁丝片午夜精品| 亚洲美女搞黄在线观看| 成年免费大片在线观看| 日本黄色片子视频| 国产又黄又爽又无遮挡在线| 小蜜桃在线观看免费完整版高清| 不卡一级毛片| 成人国产麻豆网| www.色视频.com| www日本黄色视频网| 精品国产三级普通话版| 久久久成人免费电影| 波多野结衣高清作品| 久久99蜜桃精品久久| 26uuu在线亚洲综合色| 国产精品,欧美在线| 免费观看人在逋| 精品欧美国产一区二区三| 亚洲av不卡在线观看| 联通29元200g的流量卡| 成年版毛片免费区| 成人毛片a级毛片在线播放| 男女下面进入的视频免费午夜| 亚洲精品粉嫩美女一区| 色播亚洲综合网| 精品无人区乱码1区二区| 国产又黄又爽又无遮挡在线| 搞女人的毛片| 国产亚洲精品久久久久久毛片| 国产午夜福利久久久久久| 亚洲第一区二区三区不卡| 国产女主播在线喷水免费视频网站 | 韩国av在线不卡| 国产中年淑女户外野战色| 乱系列少妇在线播放| 午夜激情欧美在线| 三级毛片av免费| 免费在线观看成人毛片| 国产亚洲欧美98| 美女高潮的动态| 男人狂女人下面高潮的视频| 99久国产av精品| 女人被狂操c到高潮| 亚洲精品456在线播放app| 男人狂女人下面高潮的视频| 亚洲av不卡在线观看| 日本在线视频免费播放| 久久久国产成人免费| 国产v大片淫在线免费观看| 天堂av国产一区二区熟女人妻| 国语自产精品视频在线第100页| 国产一级毛片七仙女欲春2| 最好的美女福利视频网| 国产成人精品婷婷| 欧美zozozo另类| 在线观看午夜福利视频| 少妇人妻精品综合一区二区 | 一级黄片播放器| 欧美日韩国产亚洲二区| 亚洲自偷自拍三级| 日韩高清综合在线| 亚洲乱码一区二区免费版| 国产大屁股一区二区在线视频| 亚洲av一区综合| 九九在线视频观看精品| 又爽又黄无遮挡网站| 老熟妇乱子伦视频在线观看| 校园春色视频在线观看| 一本一本综合久久| 一级二级三级毛片免费看| 欧美日韩精品成人综合77777| 男人狂女人下面高潮的视频| 热99在线观看视频| 99国产精品一区二区蜜桃av| 人体艺术视频欧美日本| 久久国内精品自在自线图片| 色噜噜av男人的天堂激情| 亚洲国产精品国产精品| 51国产日韩欧美| 岛国毛片在线播放| 国产亚洲精品久久久久久毛片| 激情 狠狠 欧美| .国产精品久久| 99热这里只有是精品50| 日韩 亚洲 欧美在线| 高清日韩中文字幕在线| 伊人久久精品亚洲午夜| 国产91av在线免费观看| 天堂影院成人在线观看| 一个人免费在线观看电影| 国产一区二区在线av高清观看| 久久人人爽人人片av| 99热这里只有是精品50| 九九热线精品视视频播放| 亚洲成人精品中文字幕电影| 美女被艹到高潮喷水动态| 日本与韩国留学比较| 国产黄色视频一区二区在线观看 | 黄色视频,在线免费观看| 午夜免费男女啪啪视频观看| 亚洲无线在线观看| 欧美日韩在线观看h| 国产黄片美女视频| 蜜桃久久精品国产亚洲av| 日韩欧美一区二区三区在线观看| 全区人妻精品视频| 精品久久久久久久久久久久久| 两个人视频免费观看高清| 亚洲国产精品sss在线观看| 久久久久久久久中文| 成人亚洲精品av一区二区| 免费人成视频x8x8入口观看| 国产精品免费一区二区三区在线| 亚洲中文字幕日韩| 悠悠久久av| 看非洲黑人一级黄片| 国产精品三级大全| 久久久久九九精品影院| 色视频www国产| 国产av麻豆久久久久久久| 91av网一区二区| 久久精品久久久久久噜噜老黄 | 最近最新中文字幕大全电影3| 亚洲欧美日韩高清在线视频| 亚洲欧美日韩高清在线视频| 国产精品美女特级片免费视频播放器| 18禁在线播放成人免费| 亚洲欧美日韩东京热| 成人性生交大片免费视频hd| 18禁黄网站禁片免费观看直播| 人人妻人人看人人澡| 狠狠狠狠99中文字幕| 一级黄色大片毛片| 日韩,欧美,国产一区二区三区 | 成人高潮视频无遮挡免费网站| 国产片特级美女逼逼视频| 亚洲欧美成人综合另类久久久 | 免费av毛片视频| 啦啦啦观看免费观看视频高清| 亚洲欧美精品综合久久99| 亚洲最大成人手机在线| 亚洲aⅴ乱码一区二区在线播放| ponron亚洲| 特级一级黄色大片| 一级黄片播放器| 成人午夜高清在线视频| 日韩欧美 国产精品| 亚洲av电影不卡..在线观看| 国产av不卡久久| 特级一级黄色大片| 床上黄色一级片| 又爽又黄a免费视频| 一夜夜www| 91精品一卡2卡3卡4卡| 一本久久精品| 久久这里只有精品中国| 国产精品久久久久久亚洲av鲁大| 免费看日本二区| 亚洲精品日韩av片在线观看| 网址你懂的国产日韩在线| 日韩人妻高清精品专区| 一区二区三区高清视频在线| 亚洲av免费在线观看| 激情 狠狠 欧美| 亚洲av免费高清在线观看| 亚洲七黄色美女视频| 欧美激情在线99| 国产av一区在线观看免费| 搞女人的毛片| 91aial.com中文字幕在线观看| av国产免费在线观看| 在线观看av片永久免费下载| 亚洲欧美成人综合另类久久久 | .国产精品久久| 日本黄大片高清| 久久久成人免费电影| 久久午夜亚洲精品久久| 成人美女网站在线观看视频| 久久99蜜桃精品久久| 亚洲欧美清纯卡通| 国产精品99久久久久久久久| 少妇丰满av| 日韩一区二区三区影片| .国产精品久久| 两个人视频免费观看高清| 直男gayav资源| 日韩欧美精品v在线| av卡一久久| 久久久久久久久久成人| 26uuu在线亚洲综合色| 一本一本综合久久| 你懂的网址亚洲精品在线观看 | 少妇被粗大猛烈的视频| 亚洲成a人片在线一区二区| 中文字幕av在线有码专区| 婷婷六月久久综合丁香| 好男人在线观看高清免费视频| 99久久九九国产精品国产免费| 国产精品,欧美在线| 亚洲av不卡在线观看| 欧美日韩综合久久久久久| 少妇猛男粗大的猛烈进出视频 | 成人性生交大片免费视频hd| 免费看光身美女| 在线观看av片永久免费下载| 99国产精品一区二区蜜桃av| 高清在线视频一区二区三区 | 99国产极品粉嫩在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产伦一二天堂av在线观看| 国内揄拍国产精品人妻在线| 亚洲一区高清亚洲精品| 高清午夜精品一区二区三区 | 99精品在免费线老司机午夜| 亚洲av男天堂| 欧美丝袜亚洲另类| 精品不卡国产一区二区三区| 国产一区二区亚洲精品在线观看| 日韩大尺度精品在线看网址| 婷婷色综合大香蕉| 久久草成人影院| 国产极品精品免费视频能看的| 亚洲高清免费不卡视频| 看片在线看免费视频| 爱豆传媒免费全集在线观看| 啦啦啦啦在线视频资源| 村上凉子中文字幕在线| 国产精华一区二区三区| 看十八女毛片水多多多| 男女边吃奶边做爰视频| 在线国产一区二区在线| 日韩欧美精品免费久久| 国产69精品久久久久777片| 国内精品一区二区在线观看| av视频在线观看入口| 国产色爽女视频免费观看| 丰满的人妻完整版| 日本黄色片子视频| 寂寞人妻少妇视频99o| 久久九九热精品免费| 亚洲欧美成人精品一区二区| 噜噜噜噜噜久久久久久91| 九九爱精品视频在线观看| 小说图片视频综合网站| 午夜老司机福利剧场| 又黄又爽又刺激的免费视频.| 插逼视频在线观看| 激情 狠狠 欧美| 国产在视频线在精品| 成人国产麻豆网| 熟女人妻精品中文字幕| 国产午夜福利久久久久久| 偷拍熟女少妇极品色| 国产精品不卡视频一区二区| av又黄又爽大尺度在线免费看 | 久久人人爽人人爽人人片va| 久久精品国产亚洲网站| 国产精品乱码一区二三区的特点| 日韩强制内射视频| 99久久精品国产国产毛片| 亚洲国产精品合色在线| 欧美性猛交╳xxx乱大交人| 国产成人a区在线观看| 在线免费观看的www视频| 99热这里只有是精品在线观看| 九草在线视频观看| 少妇猛男粗大的猛烈进出视频 | 亚洲av二区三区四区| 哪里可以看免费的av片| 三级毛片av免费| 国产成人freesex在线| 久久这里只有精品中国| 观看免费一级毛片| 国产三级在线视频| 中出人妻视频一区二区| 国产三级中文精品| 男人狂女人下面高潮的视频| 日韩大尺度精品在线看网址| www日本黄色视频网| 精品久久国产蜜桃| 精品久久久久久久末码| 午夜久久久久精精品| av天堂在线播放| 免费看美女性在线毛片视频| 青春草国产在线视频 | 国产黄片美女视频| 精品熟女少妇av免费看| 欧美日韩精品成人综合77777| 99九九线精品视频在线观看视频| 99精品在免费线老司机午夜| 最好的美女福利视频网| 久久久久久久久久黄片| 欧美3d第一页| 伦精品一区二区三区| 边亲边吃奶的免费视频| 夜夜爽天天搞| 午夜精品一区二区三区免费看| 亚洲人成网站在线播放欧美日韩| 观看美女的网站| 特大巨黑吊av在线直播| 国产探花极品一区二区| 在线免费观看的www视频| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线播| 成熟少妇高潮喷水视频| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 国产黄色小视频在线观看| 国内久久婷婷六月综合欲色啪| 国产午夜精品论理片| 国产精品福利在线免费观看| 最近2019中文字幕mv第一页| 精品人妻一区二区三区麻豆| 日韩视频在线欧美| 中文字幕人妻熟人妻熟丝袜美| 国产午夜福利久久久久久| 亚洲国产精品成人综合色| 久久精品国产亚洲av天美| av在线观看视频网站免费| 深夜精品福利| 蜜桃亚洲精品一区二区三区| 国内精品一区二区在线观看| 2022亚洲国产成人精品| 一本久久精品| 十八禁国产超污无遮挡网站| 国产精品乱码一区二三区的特点| av国产免费在线观看| 亚洲国产精品成人综合色| 高清毛片免费观看视频网站| 国产伦理片在线播放av一区 | 精品久久久久久久久久久久久| 成人毛片60女人毛片免费| 午夜免费激情av| 赤兔流量卡办理| 成人三级黄色视频| 一级毛片久久久久久久久女| 色尼玛亚洲综合影院| 你懂的网址亚洲精品在线观看 | 亚洲精品乱码久久久久久按摩| 国产精品三级大全| 国产高清三级在线| 亚洲第一区二区三区不卡| 国产国拍精品亚洲av在线观看| 最近2019中文字幕mv第一页| 精品少妇黑人巨大在线播放 | 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 成人无遮挡网站| 亚洲美女搞黄在线观看| 亚洲av二区三区四区| 丰满人妻一区二区三区视频av| 国产黄片视频在线免费观看| 青春草国产在线视频 | 你懂的网址亚洲精品在线观看 | 亚洲国产日韩欧美精品在线观看| 日产精品乱码卡一卡2卡三| 欧美zozozo另类| 亚洲人成网站在线观看播放| 日产精品乱码卡一卡2卡三| 久久精品影院6| 亚洲精品国产av成人精品| 国语自产精品视频在线第100页| 亚洲成人中文字幕在线播放| 日韩一区二区三区影片| 成年版毛片免费区| 久久综合国产亚洲精品| 亚洲第一电影网av| 少妇人妻一区二区三区视频| 成熟少妇高潮喷水视频| 亚洲精品成人久久久久久| 日韩欧美国产在线观看| www.色视频.com| 午夜爱爱视频在线播放| 国产午夜福利久久久久久| 我要搜黄色片| 国内精品一区二区在线观看| АⅤ资源中文在线天堂| 亚洲无线观看免费| 亚洲精品粉嫩美女一区| www.av在线官网国产| 午夜亚洲福利在线播放| 日韩一区二区视频免费看| 国产一区二区在线观看日韩| а√天堂www在线а√下载| 最近的中文字幕免费完整| 精品一区二区三区视频在线| 亚洲图色成人| 三级经典国产精品| 日韩欧美一区二区三区在线观看| 午夜亚洲福利在线播放| 亚洲国产日韩欧美精品在线观看| 亚洲欧洲国产日韩| 免费搜索国产男女视频| av黄色大香蕉| 久久久成人免费电影| 亚洲人与动物交配视频| 久久久欧美国产精品| 亚洲精品乱码久久久久久按摩| 亚洲人与动物交配视频| 免费av观看视频| 少妇被粗大猛烈的视频| 晚上一个人看的免费电影| 成年版毛片免费区| 国产亚洲91精品色在线| 久久久精品大字幕| 精品久久久久久久久亚洲| 欧美丝袜亚洲另类| 伦精品一区二区三区| 亚洲综合色惰| 老师上课跳d突然被开到最大视频| 亚洲色图av天堂| 中国国产av一级| 免费人成视频x8x8入口观看| 美女内射精品一级片tv| 亚洲精品国产av成人精品| 国产激情偷乱视频一区二区| 91久久精品电影网| 欧美一级a爱片免费观看看| 日本一二三区视频观看| 亚洲精品成人久久久久久| 哪里可以看免费的av片| 午夜老司机福利剧场| 我要搜黄色片| 99久久精品国产国产毛片| 日韩中字成人| 人妻制服诱惑在线中文字幕| 亚洲精品自拍成人| 亚洲国产欧洲综合997久久,| 精品免费久久久久久久清纯| 国产黄色视频一区二区在线观看 | 欧美不卡视频在线免费观看| www日本黄色视频网| 亚洲aⅴ乱码一区二区在线播放| 最近最新中文字幕大全电影3| 少妇熟女aⅴ在线视频| 直男gayav资源| 亚洲欧美精品专区久久| 日韩亚洲欧美综合| 波多野结衣高清无吗| 国产黄色视频一区二区在线观看 | 成年女人看的毛片在线观看| 亚洲第一区二区三区不卡| 久久久a久久爽久久v久久| 日韩一本色道免费dvd| 久久久精品欧美日韩精品| 精品日产1卡2卡| 欧美+亚洲+日韩+国产| 嫩草影院新地址| 亚洲人成网站高清观看| 国产精品,欧美在线| 欧美丝袜亚洲另类| 大香蕉久久网| 97超视频在线观看视频| 女人十人毛片免费观看3o分钟| 最近2019中文字幕mv第一页| 中文亚洲av片在线观看爽| 亚洲无线在线观看| 成熟少妇高潮喷水视频| 国产精品综合久久久久久久免费| 深爱激情五月婷婷| 欧美+亚洲+日韩+国产| 午夜精品在线福利| 网址你懂的国产日韩在线| 国产片特级美女逼逼视频| 蜜臀久久99精品久久宅男| 亚洲成人久久爱视频| 亚洲国产精品sss在线观看| 午夜爱爱视频在线播放| 国产精品不卡视频一区二区| 99热网站在线观看| 精品国内亚洲2022精品成人| 高清日韩中文字幕在线| 欧美最黄视频在线播放免费| 深夜a级毛片| 久99久视频精品免费| 国产探花极品一区二区| 亚洲婷婷狠狠爱综合网| 亚洲精华国产精华液的使用体验 | 国产男人的电影天堂91| 嘟嘟电影网在线观看| www.色视频.com| 亚洲国产精品久久男人天堂| 人妻系列 视频|