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

    一次梅雨鋒暴雨過程數(shù)值模擬的云微物理參數(shù)化敏感性研究

    2021-12-14 07:15:50周志敏崔春光胡揚(yáng)康兆萍
    大氣科學(xué) 2021年6期
    關(guān)鍵詞:云滴雨滴極值

    周志敏 崔春光 胡揚(yáng) 康兆萍

    1 中國氣象局武漢暴雨研究所暴雨監(jiān)測預(yù)警湖北省重點(diǎn)實(shí)驗(yàn)室,武漢 430205

    2 中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室,北京 100081

    1 引言

    梅雨鋒暴雨是我國長江中下游地區(qū)夏季汛期的主要?dú)庀鬄?zāi)害,其具有持續(xù)時(shí)間長和突發(fā)性強(qiáng)的特征,往往產(chǎn)生大范圍洪澇災(zāi)害,給該流域經(jīng)濟(jì)和人民生命財(cái)產(chǎn)安全帶來嚴(yán)重威脅和損失。由于對大氣環(huán)流的短時(shí)突變、強(qiáng)對流天氣系統(tǒng)發(fā)生、發(fā)展的機(jī)理缺乏深入了解,故針對暴雨的落區(qū)、突發(fā)時(shí)間、強(qiáng)度突變等預(yù)報(bào)難度較大。

    針對梅雨鋒暴雨的研究,過去幾十年的工作積累了不少成果。有些研究側(cè)重于動(dòng)力、熱力條件的資料分析和模擬研究,如環(huán)流特征,溫度場,風(fēng)場結(jié)構(gòu)、中尺度擾動(dòng)和水分及能量收支、不穩(wěn)定理論和位渦與輻射冷卻的診斷分析、高低空急流的耦合作用等(孫淑清, 1982; Yu et al., 1988; 祁秀香等,2015; 楊舒楠等, 2017);有些研究側(cè)重于梅雨鋒暴雨系統(tǒng)的結(jié)構(gòu)特征及形成與維持(王建捷和陶詩言, 2002; 鄭婧等, 2015; 諶偉等, 2017)和梅雨鋒上的渦旋系統(tǒng)(陸漢城等, 2002; 谷文龍等, 2010; 劉夢娟等, 2017)。越來越多的研究人員發(fā)現(xiàn),要提高梅雨鋒暴雨的預(yù)報(bào)水平,必須重視其中的云微物理過程(孫建華和趙思雄, 2000)。因?yàn)樵莆⑽锢磉^程對環(huán)境的動(dòng)力、熱力過程有著不可忽視的反饋?zhàn)饔?,尤其是水成物相態(tài)變化釋放(吸收)潛熱對熱力和動(dòng)力過程有著重要影響。

    目前為止,在沒有針對云系相關(guān)的非常規(guī)探測資料時(shí),除了降水之外,主要從宏觀熱力和動(dòng)力方面對模式結(jié)果進(jìn)行評估檢驗(yàn)。然而,這是不完善的,模式的云的微物理屬性同樣需要與觀測進(jìn)行比較。由于觀測的局限性,無論是利用飛機(jī)、衛(wèi)星,還是云雷達(dá)等非常規(guī)探測設(shè)備,都無法針對特定云系做完美的跟蹤觀測。因此,梅雨鋒暴雨系統(tǒng)中的云微物理過程無法追蹤。由于模式結(jié)果具有時(shí)間和空間上的連續(xù)性,就研究云的屬性及變化過程而言,數(shù)值模式可以作為一個(gè)重要的輔助工具。當(dāng)然,該工具也存在一定不足。因?yàn)樵莆⑽锢韰?shù)化是影響天氣過程模擬效果的重要因素之一,但云微物理參數(shù)化本身卻存在不確定性,且該不確定性至少由三個(gè)原因產(chǎn)生:

    首先,云微物理參數(shù)化的不確定性來自對水成物屬性的定義。由于冰相粒子形狀復(fù)雜,密度范圍較廣,對其增長有著重要影響(Morrison and Milbrandt, 2015)。根據(jù)冰相粒子的處理方式不同,云微物理參數(shù)化方案可以被分為兩類。一類是傳統(tǒng)方式,水成物被分為云滴、雨滴、冰晶、雪和霰(雹)等(Lin et al., 1983; Rutledge and Hobbs,1984),并且水成物的特征,如群密度、下落末速度、粒子形狀等都已經(jīng)在方案中預(yù)設(shè);云滴和雨滴均假設(shè)為球體,盡管雨滴在空中運(yùn)動(dòng)時(shí)存在形變,但大體上假設(shè)是合理的;冰相粒子則相對復(fù)雜,尤其是雪粒子,存在盤狀、針狀、多邊形等形態(tài)。有的方案將雪粒子視為球體(Morrison et al., 2005),有的方案將雪粒子視為球體和非球體進(jìn)行考慮(Milbrandt and Yao, 2005a, 2005b),有的方案對雪粒子分檔進(jìn)行處理(Thompson et al., 2008);同時(shí),由于云微物理參數(shù)化方案考慮不同冰相粒子之間的相互轉(zhuǎn)換,預(yù)設(shè)的冰相粒子屬性差異在此過程中,必將產(chǎn)生新的不確定性(Morrison and Milbrandt, 2015)。另一類云微物理參數(shù)化方案將冰相粒子合成一類進(jìn)行考慮(Lin and Colle, 2011;Morrison and Milbrandt, 2015),以期減少不同冰相粒子之間相互轉(zhuǎn)換帶來的不確定性。但是,該方法帶來的水成物聚合、凇附效率等方面的不確定性對數(shù)值模式的性能產(chǎn)生影響(Morrison and Milbrandt, 2015)。

    其次,云微物理參數(shù)化的不確定性來自不同方案對云微物理過程處理的差異。有的方案關(guān)注大尺度抬升和中尺度輻合造成的鋒面降水(Rutledge and Hobbs, 1984),有的方案側(cè)重于積雨云的降水過程(胡志晉和何觀芳, 1987),有的方案側(cè)重于冰雹云(許煥斌和王思微, 1985; 孔凡鈾等, 1990,郭學(xué)良等, 2001),有的方案側(cè)重于積層混合云(劉公波等, 1994; 樓小鳳, 2002),有的方案強(qiáng)調(diào)云凝結(jié)核和冰核的重要性及其對長短波輻射的影響(Morrison et al., 2009)。由于研究側(cè)重點(diǎn)不同,故參數(shù)化方案中對云微物理過程的處理方式也有所差異。為了比較不同云微物理參數(shù)化方案的結(jié)果,研究人員做了大量工作:孫晶和王鵬云(2003)認(rèn)為MM5 模式中的Reisner 霰方案對華南暴雨具有較好的模擬能力。樓小鳳(2002)等對比并詳細(xì)分析了Goddard、Reisner 和Schultz 方案的物理并指出,模擬臺風(fēng)最低氣壓和地面最大風(fēng)速時(shí),Reisner方案并不存在優(yōu)勢。閆之輝和鄧蓮堂(2007)對比分析了各種微物理參數(shù)化方案對降水預(yù)報(bào)的影響,認(rèn)為Lin 方案總體預(yù)報(bào)結(jié)果較好。尹金方(2014)調(diào)研了十余年來中尺度模式云微物理參數(shù)化方案在中國范圍內(nèi)敏感性試驗(yàn)研究成果后也認(rèn)為,WRF模式中Lin 方案的效果較好。而朱格利等(2014)利用WRF 模式中不同云微物理參數(shù)化方案對華南一次暴雨過程進(jìn)行數(shù)值模擬和性能分析后發(fā)現(xiàn),Lin 方案的模擬效果最差。因此,不同云微物理參數(shù)化方案的模擬性能隨著個(gè)例的不同而呈現(xiàn)出差異。

    最后,云微物理參數(shù)化的不確定性來自方案的設(shè)計(jì)思路。嚴(yán)格意義上來講,一個(gè)完美的數(shù)值模式除了準(zhǔn)確描述宏觀熱、動(dòng)力場之外,還必須完美描述微觀場的變化過程。許煥斌和尹金方(2017)提出了“人工影響天氣數(shù)值模式”的概念,因?yàn)槠渌鼣?shù)值模式(氣候、環(huán)境、天氣及單一云—降水微物理模式等)一般只關(guān)注云場、地面降水、潛熱加熱項(xiàng)等,而人工影響天氣模式關(guān)心云體內(nèi)部結(jié)構(gòu),需要顯式地描述云。云體結(jié)構(gòu)及其演變是宏微觀場相互作用的結(jié)果,參數(shù)化方法需具有動(dòng)態(tài)描云體演變的功能。本質(zhì)上,這個(gè)概念不僅僅適用于“人工影響天氣數(shù)值模式”,也應(yīng)該適用于其它模式。而事實(shí)是,現(xiàn)有模式的云微物理參數(shù)化方案往往是一個(gè)靜態(tài)的概念,云體演變過程中,參數(shù)化方法固定不變,如很多云微物理參數(shù)化方案將部分水成物轉(zhuǎn)換臨界值取為定常值,這種不足主要由于現(xiàn)階段研究人員對降水系統(tǒng)里面的云微物理過程認(rèn)識不足所致。隨著探測技術(shù)的逐步提升,對云微物理過程認(rèn)識的逐步深入,有些武斷的假設(shè)將參照大量觀測資料,重新進(jìn)行嚴(yán)謹(jǐn)?shù)貐?shù)化表述。

    盡管云微物理參數(shù)化方案存在諸多不確定性,但在現(xiàn)有條件下,為了做好云微物理參數(shù)化方案的本地化工作,提高梅雨鋒暴雨的預(yù)報(bào)水平,我們有必要深入了解不同方案中的云微物理參數(shù)化過程差異。目前,很多研究人員已經(jīng)通過水成物的源匯項(xiàng)收支比較了相關(guān)降水機(jī)制(Cui and Li, 2006; Wang et al., 2007; Cui, 2008; Huang et al., 2019),但相關(guān)研究主要聚焦于熱帶降水系統(tǒng),在梅雨鋒暴雨方面的分析相對較少。

    本文期望從水成物源匯項(xiàng)差異的角度出發(fā),通過一次梅雨鋒暴雨過程的數(shù)值模擬研究,分析不同云微物理參數(shù)化方案對梅雨鋒暴雨的影響,從而為改善云微物理方案做一些前期工作。本文將在第二部分介紹使用數(shù)據(jù)和方法,第三部分介紹本文所用個(gè)例的概況,第四部分對模擬結(jié)果進(jìn)行討論,最后在第五部分進(jìn)行總結(jié)和討論。

    2 數(shù)據(jù)和方法

    2.1 數(shù)據(jù)

    本文使用的降水?dāng)?shù)據(jù)是CMORPH,分辨率為0.1°×0.1°(王新華等, 2006),該數(shù)據(jù)已經(jīng)過質(zhì)量控制。歐洲中心的第五代再分析資料ERA5 作為模式的背景場,該數(shù)據(jù)水平分辨率為0.25°×0.25°,時(shí)間分辨率為1 h,垂直方向共37 層(從1000 hPa到10 hPa)。

    2.2 研究方法

    本文的數(shù)值模擬采用了中尺度數(shù)值模式WRF的V3.4.1 版本,雙層嵌套(見圖1),中心經(jīng)緯度為(31.5°N,110°E)。為盡可能減少模式的“spinup”現(xiàn)象,模式運(yùn)行從2018 年6 月29 日12:00(協(xié)調(diào)世界時(shí),下同)開始,止于6 月30 日11:00。并且為了對背景場進(jìn)行約束,對外層網(wǎng)格采用了“nudging-微調(diào)逼近”約束,約束的時(shí)間周期和模式運(yùn)行周期同步。模式采用了三個(gè)云微物理參數(shù)化方案(Morrison and Milbrandt, 2005; Milbrandt and Yao, 2005a, 2005b; Thompson et al., 2008)進(jìn)行比較和評估。模式結(jié)果1 小時(shí)輸出一次。模式設(shè)置如表1。

    圖1 模式模擬區(qū)域,紅色實(shí)心圓表示模式外層區(qū)域中心點(diǎn)Fig. 1 Geographical domain used in WRF model runs. The red solid circle indicates the central point of the outer domain

    表1 模式參數(shù)設(shè)置Table 1 The settings of model parameters

    2.3 分析方法

    本文選擇了梅雨鋒背景下的一次過境湖北省的強(qiáng)降水過程進(jìn)行分析,分析的側(cè)重點(diǎn)在于導(dǎo)致降水差異的不同云微物理方案中的云微物理過程差異。由于不同云微物理參數(shù)化方案的主要差異在于對冰相粒子的考慮(Morrison and Milbrandt, 2015),為減少冰相粒子之間相互轉(zhuǎn)換帶來的復(fù)雜性,本文首先將冰相粒子作為一個(gè)整體進(jìn)行考慮,比較不同方案的差異,然后再分析雪粒子和霰粒子與云滴、雨滴之間的碰并過程。根據(jù)分析需求,本文考慮的云微物理過程見表2。

    表2 本文所用的微物理過程Table 2 The microphysical source/sink terms

    3 概況(天氣形勢)

    本次過程始于2018 年6 月28 日22:00,時(shí)斷時(shí)續(xù),一直持續(xù)到7 月1 日。該過程在湖北境內(nèi)的連續(xù)時(shí)段為2018 年6 月29 日22:00 至30 日11:00。30 日12:00~16:00,雨帶南壓到湖南省北部,但是17:00 開始,湖北東部又開始發(fā)生明顯降水。為了持續(xù)跟蹤降水系統(tǒng),選擇2018 年6 月29 日22:00 至30 日11:00 進(jìn)行考慮。在該時(shí)段,最大累積降水強(qiáng)度為189 mm(基于CMORPH 數(shù)據(jù))。圖2 為2018 年6 月29 日22:00 至30 日11:00 基于ERA5 再分析場的不同高度層的風(fēng)場及水汽場等的平均分布情況,各層環(huán)流及其配置有利于湖北區(qū)域產(chǎn)生強(qiáng)降水過程。

    如圖2a 所示,200 hPa 高度上,湖北、河南一帶位于高空急流入口區(qū)右側(cè),處于次級環(huán)流的高空輻散區(qū),對應(yīng)次級環(huán)流的上升支。我國中東部地區(qū)低層則主要受到西南季風(fēng)的影響(圖2b 和c)。從圖2b 可知,700 hPa 高度上,貴州北部存在溫度大值區(qū),湖北南部及西部存在一條切變線,在暖平流的配合下,有利于上升運(yùn)動(dòng)的產(chǎn)生和不穩(wěn)定層結(jié)的建立。從圖2c 可知,850 hPa 高度上,湖北省受多個(gè)方向的水汽通量影響,主要來源是西南氣流輸送的水汽,其在湖北中部(31°N,112°~114°E)形成很強(qiáng)的輻合中心。湖北省北部的水汽來自東部,但對應(yīng)水汽的輻散區(qū)。

    圖2 2018 年6 月29 日22:00 至30 日11:00(協(xié)調(diào)世界時(shí),下同)基于ERA5 再分析場的相關(guān)物理量時(shí)間平均分布:(a)200 hPa 風(fēng)矢量(單位:m s-1)和緯向風(fēng)速(填色,單位:m s-1);(b)700 hPa 風(fēng)矢量(單位:m s-1)和溫度(填色,單位:°C);(c)850 hPa 水汽通量(單位:g s-1 hPa-1 cm-1)和水汽通量散度(填色,單位:10-7 g s-1 hPa-1 cm-2)Fig. 2 The time means (2200 UTC 29 June to 1100 UTC 30 June 2018) of (a) wind vectors (units: m s-1) and zonal wind velocity (shaded; units: m s-1) at 200 hPa, (b) wind vectors (units: m s-1) and air temperature (units: °C) at 700 hPa, and (c) moisture flux (units: g s-1 hPa-1 cm-1) and moisture flux divergence (shaded; units: 10-7 g s-1 hPa-1 cm-2) at 850 hPa based on the ERA5 reanalysis

    4 模擬結(jié)果分析

    4.1 天氣形勢模擬比較

    圖3 為2018 年6 月29 日22:00 至30 日11:00基于不同云微物理參數(shù)化方案、不同高度層的風(fēng)場及水汽場的平均分布情況。對比圖2 和圖3 可知,三個(gè)方案模擬的熱、動(dòng)力場,大體上和再分析場相近。從圖3a1、b1 和c1 可知,三個(gè)方案均模擬出了200 hPa 上的西風(fēng)急流,且湖北省均處于高空急流入口區(qū)右側(cè)。Thompson 方案中,在安徽和江蘇北部,西風(fēng)急流較再分析場更強(qiáng),也強(qiáng)于其它兩個(gè)方案。同時(shí),在湖北中部地區(qū)的風(fēng)場較再分析場及其它兩個(gè)方案略弱。在700 hPa 高度上(見圖3a2-c2),Thompson 方案的溫度梯度較大,溫度平流略強(qiáng),其次是Morrison 方案,MY 方案相對更弱。此外,三個(gè)方案模擬的850 hPa 水汽通量差別不大,水汽來源主要是西南季風(fēng)輸送;但水汽輻合區(qū)的位置則存在明顯差異。Morrison 方案中,最強(qiáng)的水汽輻合區(qū)呈“水平分布”,位于湖北中部偏南(30°N,112°~114°E)。Thompson 方案中的主要水汽輻合區(qū)大體呈“西北—東南”方向分布(30°~32°N,112°~114°E)。MY 方案中的水汽輻合區(qū)相對其它兩個(gè)方案來說,較為零散,在湖北區(qū)域沒有明顯的水汽輻合帶。

    圖3 同圖2,但為不同方案的模擬結(jié)果:(a1-a3)Morrison 方案;(b1-b3)Thompson 方案;(c1-c3)MY 方案Fig. 3 Same as Fig. 2 but for the model output based on the (a1-a3) Morrison, (b1-b3) Thompson, and (c1-c3) MY microphysical schemes

    綜上所述,三個(gè)方案均較好再現(xiàn)了本次過程的天氣尺度熱力、動(dòng)力背景場,因此,我們要分析的不同云微物理參數(shù)化方案對局地降水影響的可信度也較高。三個(gè)方案之所以能大致模擬出大尺度過程,原因之一可能是因?yàn)樵谠搮^(qū)域使用的是顯式云微物理參數(shù)化方案,而顯式方案對相對均勻的大尺度云系具有較好的描述能力(許煥斌和尹金方, 2017)。顯式方案對非均勻的中小尺度對流系統(tǒng)描述能力弱于隱式方案,當(dāng)系統(tǒng)中的非均勻性中小尺度系統(tǒng)較強(qiáng)時(shí),大尺度系統(tǒng)的描述就會相對較差。因此,三個(gè)方案的結(jié)果中,與大尺度場存在差異的區(qū)域,部分原因可能就是顯示方案的不足所致(嚴(yán)格來說,這是一個(gè)綜合的結(jié)果,還與熱、動(dòng)力場及其它物理過程對云微物理過程的反饋有關(guān))。

    圖4 為2018 年6 月29 日22:00 至30 日11:00不同云微物理參數(shù)化方案模擬的700 hPa 位勢高度場的標(biāo)準(zhǔn)差及垂直氣流。從圖4a1-c1 可知,Thompson方案中,700 hPa 位勢高度場的標(biāo)準(zhǔn)差在湖北中部存在一個(gè)大值中心,而其它兩個(gè)方案的結(jié)果則沒有明顯的中心,且在該時(shí)段內(nèi)湖北省內(nèi)的高度場的標(biāo)準(zhǔn)差基本小于10 gpm。這說明Thompson 方案中的局地環(huán)流變化(例如系統(tǒng)的發(fā)展、消亡或移動(dòng)等)更劇烈。相應(yīng)地,低層的垂直速度更大(圖4b2)。值得注意的是,對比圖4a2-c2 與圖3a3-c3 可知,垂直速度的大值區(qū)除了和位勢高度標(biāo)準(zhǔn)差的大值區(qū)相對應(yīng)之外,和水汽輻合區(qū)也存在較好的對應(yīng)關(guān)系。

    圖4 2018 年6 月29 日22:00 至30 日11:00(a1、a2)Morrison 方案、(b1、b2)Thompson 方案、(c1、c2)MY 方案模擬700 hPa 位勢高度場的標(biāo)準(zhǔn)差(第一行,單位:gpm)及垂直氣流(第二行,單位:m s-1)空間分布Fig. 4 Spatial distribution of the standard deviation of the geopotential height (top line, units: gpm) and vertical airflow (bottom line, units: m s-1) at 700 hPa during 2200 UTC 29 June to 1100 UTC 30 June 2018 based on the model output of (a1, a2) Morrison, (b1, b2) Thompson, and (c1, c2) MY microphysical schemes

    圖5 為湖北地區(qū)(選擇108.35°E~116.25°E 范圍)的平均經(jīng)向環(huán)流及西風(fēng)的垂直剖面圖。通過圖3 和圖5 可知,三個(gè)方案中,高空西風(fēng)急流的南側(cè)都存在著上升氣流,但高空急流下方上升氣流并不強(qiáng)。Thompson 中的垂直氣流在中低層最強(qiáng),且西風(fēng)氣流大值區(qū)的延伸范圍最大,向下伸展高度最低(圖5b)。Morrison 方案中的上升氣流和MY方案差別不大(圖5a 和c),前者略強(qiáng),低層西風(fēng)急流伸展范圍也略廣。

    圖5 2018 年6 月29 日22:00 至30 日11:00(a)Morrison 方案、(b)Thompson 方案、(c)MY 方案模擬湖北地區(qū)(108.35°E~116.25°E)的時(shí)間平均經(jīng)向環(huán)流(風(fēng)矢圖,單位:m s-1)及西風(fēng)(填色,單位:m s-1)的垂直剖面Fig. 5 Vertical section of the time-averaged meridional circulation (vectors, units: m s-1) and zonal wind (shaded, units: m s-1) based on the model output of (a) Morrison, (b) Thompson, and (c) MY microphysical schemes during 2200 UTC 29 June to 1100 UTC 30 June 2018

    4.2 降水及水成物分析

    圖6 為2018 年6 月29 日22:00 至30 日11:00的累積降水的觀測和模擬結(jié)果。紅色方框基本覆蓋整個(gè)湖北省,即為本文的重點(diǎn)研究范圍。從圖6 可知,三個(gè)方案均模擬出了降水的覆蓋區(qū)域,但同時(shí)覆蓋范圍更廣。在紅色方框范圍(下文分析均基于此區(qū)域開展),從累積降水大于75 mm 的范圍比較來看,Thompson 方案模擬的降水最強(qiáng),Morrison和MY 差別不太大。不同方案之間的降水極值區(qū)分布差異較大,Morrison 方案模擬得到的強(qiáng)降水區(qū)覆蓋范圍與觀測更為吻合。

    圖7 為模擬和觀測逐時(shí)降水強(qiáng)度(Rain Rate,簡稱RR)的概率分布圖。從圖7 可知,三個(gè)方案的模擬結(jié)果均高估了雨強(qiáng)。Thompson 方案的模擬結(jié)果在RR<60 mm h-1范圍基本最大,但超過該范圍后,MY 的降水模擬值最大。當(dāng)RR 處于30 到35 mm h-1范圍,MY 方案的模擬值小于Morrison方案。當(dāng)RR 在其它范圍,Morrison 則小于MY 方案。RR 的觀測值最大約為55 mm h-1,但三個(gè)方案的RR 模擬最大值幾乎為觀測值的一倍。

    圖8 為2018 年6 月29 日22:00 至30 日11:00降水的觀測和模擬的區(qū)域(圖6 所示紅色方框區(qū)域,下同)平均值時(shí)間變化。從圖8 可知,模擬和實(shí)測降水在大部分時(shí)段較為相似(如2018 年6 月29日22:00 至30 日04:00,6 月30 日09:00 至11:00),但模擬降水明顯強(qiáng)于觀測結(jié)果。雨滴含量與降水的時(shí)變曲線也大體相似,但其峰值出現(xiàn)時(shí)間并不一致,后者出現(xiàn)峰值的時(shí)間更早。從三個(gè)方案的比較結(jié)果來看,Thompson 方案模擬的降水最強(qiáng)(與圖6 保持一致),其次是Morrison 方案,MY 方案模擬的降水相對最弱,但其在部分時(shí)段強(qiáng)于Morrison方案,其變化相對劇烈。Thompson 方案模擬的雨滴含量在部分時(shí)段少于Morrison 方案,MY 方案模擬的雨滴含量在整個(gè)模擬時(shí)段均最小。雨強(qiáng)和雨滴含量之所以不完全匹配,原因之一是雨滴含量是高度上的積分值,而雨強(qiáng)僅僅是模式最底層的結(jié)果。另一個(gè)原因是,雨滴的分布還受動(dòng)力場影響,如平流項(xiàng)及湍流混合項(xiàng)均對其分布有著一定作用。盡管如此,兩者在大部分時(shí)段依然保持相似變化。

    圖8 2018 年6 月29 日22:00 至30 日11:00 降水的(a)小時(shí)區(qū)域(圖6 紅色方框所示區(qū)域)平均(單位:mm)及(b)雨滴含量(單位:g m-3)平均值的時(shí)間變化曲線Fig. 8 Time series of the (a) area-averaged (area indicated by red rectangle in Fig.6) hourly rainfall intensity (units: mm) and (b) rainwater content(units: g m-3) during 2200 UTC 29 June to 1100 UTC 30 June 2018

    從累計(jì)降水量值及其分布區(qū)域來說,Morrison方案的模擬效果相對最優(yōu);從小時(shí)降水來說,Morrison 和Thompson 方案在不同時(shí)段分別占優(yōu),但模擬效果均強(qiáng)于Thompson 方案。因此,綜合來看,Morrison 方案在本次模擬過程中效果相對最優(yōu)。

    圖9 為三個(gè)方案水成物含量(小時(shí)累積)區(qū)域平均值隨時(shí)間和高度的分布。盡管圖8 中三個(gè)方案的平均降水演變趨勢具有相似性,但水成物的演變卻存在明顯差別。從云滴分布來看(圖9a1-c1),三個(gè)方案不僅極值中心出現(xiàn)的時(shí)間有差異,在分布高度上也存在差異。Morrison 方案中,云滴的5 g m-3等值線超過-10°C 以上高度,極值為17.3 g m-3;Thompson 方案中,該等值線在-10°C 高度層以下,極值為15.65 g m-3,與Morrison 方案差別不大。而MY 方案中,云滴的5 g m-3等值線伸展到了300 hPa 高度附近(-30°C)。并且,極值為26.4 g m-3,明顯大于其它兩個(gè)方案。對三個(gè)方案模擬的雨滴來說,其5 g m-3等值線均在零度層以上,但三個(gè)方案中的雨滴含量演變及極值也存在較大差異,Morrison、Thompson、和MY 方案的雨滴含量極值分別為37.5 g m-3、57.0 g m-3和34.4 g m-3。

    冰相粒子的演變及垂直分布比液態(tài)粒子的差異更大。Thompson 方案中模擬的冰粒子最少(圖9b1),整個(gè)時(shí)段、所有高度上冰粒子平均含量均小于3 g m-3的區(qū)域;Morrison 方案中,存在小部分大于3 g m-3的區(qū)域;MY 中的冰粒子分布高度及時(shí)間跨度遠(yuǎn)大于其它兩個(gè)方案,且最大冰粒子平均含量為45.4 g m-3,位于-40°C 高度以上。在雪粒子的演變過程中,MY 方案中存在兩個(gè)極值中心(圖9b3),其它方案只有一個(gè)。然而,Morrison 和Thompson 方案中,極值中心的強(qiáng)度及分布高度區(qū)間存在較大差異(圖9b1 和b2)。Morrison 方案中的雪粒子平均含量極值所在高度位于0~-10°C 區(qū)間,且極值不超過40 g m-3,而Thompson 中的雪粒子平均含量極值所在高度位于-10°C~-20°C 區(qū)間,且極值超過100 g m-3。MY 方案中的雪粒子平均含量極值與Morrison 方案模擬結(jié)果相近。就霰粒子而言,三個(gè)方案的模擬結(jié)果差異也較大。Morrison 方案中,霰粒子平均含量在兩個(gè)時(shí)段存在極值,最大極值為37.8 g m-3;Thompson 方案中的霰粒子平均含量最少,僅為16.4 g m-3;MY 方案中的霰粒子平均含量略少于Morrison 方案,為36.8 g m-3。

    圖9 2018 年6 月29 日22:00 至30 日11:00 Morrison 方案(左列)、Thompson 方案(中間列)以及MY 方案(右列)中水成物的區(qū)域平均值隨時(shí)間演變:冰晶(第一行);雪粒子(第二行);霰粒子(第三行)。填色表示冰相粒子含量(單位:g m-3);黑色實(shí)線為平均溫度,等值線從0°C~-40°C,間隔10°C;紅色實(shí)線表示云滴含量,等值線從5~35 g m-3,間隔5 g m-3;藍(lán)色實(shí)線表示雨滴含量,等值線從5~55 g m-3,間隔5 g m-3。Fig. 9 Time evolution (during 2200 UTC 29 June to 1100 UTC 30 June 2018) of area-averaged vertical profiles of ice (top line), snow (second line),and graupel (bottom line) hydrometeors by Morrison (left column), Thompson (middle column), and MY (right column) scheme, respectively. The shaded area indicates the distribution of ice phase hydrometeors (units: g m-3). The black solid lines indicates the average temperature, and the contours are from 0°C to -40°C with 10°C of intervals. The red dashed lines indicate the cloud droplet content (units: g m-3), and the contours are from 5 g m-3 to 35 g m-3 with 5 g m-3 of intervals. The blue dashed lines indicate the raindrop content (units: g m-3), and the contours are from 5 g m-3 to 55 g m-3 with 5 g m-3 of intervals

    綜上所述,三個(gè)方案模擬的水成物存在較大差異,Thompson 方案中雪粒子含量明顯偏多,冰晶和霰相對較少;MY 方案中冰晶粒子含量最多,其次為Morrison 方案,兩者的霰粒子含量均多于雪粒子。液態(tài)粒子(尤其云滴)的分布及演變均存在較大差異,由于過冷水是影響凇附過程等的關(guān)鍵因素,因此,該差異是冰相粒子分布及演變存在不同的重要原因之一。

    4.3 液態(tài)粒子主要源匯項(xiàng)特征分析

    雨滴含量差異是導(dǎo)致降水差異的最直接原因。從圖8 可知,通過雨滴含量的分析,能大體上了解雨強(qiáng)的變化趨勢,為了分析雨滴含量的變化趨勢,特對雨滴的源匯項(xiàng)進(jìn)行了分析。圖10 表示雨滴的五個(gè)主要源匯項(xiàng)在空間平均的時(shí)變曲線。從圖10可知,冰相粒子融化項(xiàng)Mlt 及雨滴對云滴的搜集項(xiàng)CLcr 項(xiàng)是雨滴的核心源項(xiàng),云滴向雨滴的自動(dòng)轉(zhuǎn)化項(xiàng)Auto 相對于其它項(xiàng)來說基本忽略不計(jì)。并且,三個(gè)方案中的CLcr 項(xiàng)都大于Mlt 項(xiàng)。Morrison方案模擬的CLcr 項(xiàng)極值最大,其次是Thompson方案,MY 中的CLcr 項(xiàng)極值最小。Thompson 方案中的雨滴蒸發(fā)和凇附效應(yīng)最強(qiáng),其次是Morrison方案,MY 中的雨滴蒸發(fā)項(xiàng)Evap 最弱。除云滴自動(dòng)轉(zhuǎn)化外,其它四個(gè)源匯項(xiàng)隨時(shí)間演變大體上有著相似性,但三個(gè)方案也存在差異。Morrison 和MY 方案中,四個(gè)主要源匯項(xiàng)大致在2018 年6 月30 日07:00 達(dá)到極值,而Thompson 方案中,四個(gè)主要源匯項(xiàng)達(dá)到極值的時(shí)間為2018 年6 月30 日08:00。

    從圖10 可知,盡管不同方案的主要雨滴源匯項(xiàng)量值有些差異,但其隨時(shí)間的演變?nèi)跃哂幸欢ㄏ嗨菩浴D11 為不同方案中雨滴主要源匯項(xiàng)總量區(qū)域平均的垂直分布比較(這里考慮了源匯項(xiàng)的正負(fù),圖10 主要考慮量值比較)。由圖11 可知,三個(gè)方案的融化總量極值所在高度一致,Morrison 模擬的融化量極值最大,Thompson 方案略小,MY 方案最小。在冰相粒子對雨滴的搜集作用上,三個(gè)方案出現(xiàn)較大差異。Morrison 和MY 均認(rèn)為,在碰并過程中,雨滴是冰相粒子增長的源,從量值來看,Morrison 方案中雨滴貢獻(xiàn)明顯大于MY 方案,但其極值所在高度更低。而Thompson 方案則認(rèn)為,零度層以上,雨滴作為冰相粒子的增長源;而零度層以下,冰相粒子和雨滴的碰并將成為雨滴增長的源,即雨滴搜集到的冰相粒子將轉(zhuǎn)化為與雨滴。三個(gè)方案的雨滴搜集云滴過程形態(tài)也較為相似,但Thompson 方案中,該過程最強(qiáng),MY 方案中最弱。Morrison 和MY 方案中的雨滴蒸發(fā)項(xiàng)垂直分布較為相似,從零度層高度往下直到750 hPa 左右高度,蒸發(fā)逐步增強(qiáng)。750 hPa 高度以下,蒸發(fā)變化較小。Thompson 方案的雨滴蒸發(fā)過程則明顯強(qiáng)于其它兩個(gè)方案,且從零度層往下至近地層,蒸發(fā)量逐步增大,這與其它兩個(gè)方案存在顯著差異。

    圖10 2018 年6 月29 日22:00 至30 日11:00 的主要雨滴源匯項(xiàng)空間平均值(單位:g m-3)的時(shí)間變化:(a)Morrison 方案;(b)Thompson方案;(c)MY 方案Fig. 10 Time series of the volumetrically averaged content (units: g m-3) of the main source and sink terms of rain during 2200 UTC 29 June to 1100 UTC 30 June 2018 by (a) Morrison, (b) Thompson, and (c) MY schemes

    圖11 不同方案的雨滴主要源匯項(xiàng)總量(單位:g m-3)區(qū)域平均值的的垂直分布:(a)Morrison 方案;(b)Thompson 方案;(c)MY方案Fig. 11 Vertical profiles of area-averaged content (unit: g m-3) of the key sink and source termsof raindrop by (a) Morrison, (b) Thompson, and (c)MY schemes

    圖12 為不同方案的Mlt、CLcr 及Evap 項(xiàng)在湖北省區(qū)域的高度平均分布。對比圖6 可知,圖12中三個(gè)關(guān)鍵源匯項(xiàng)的大值區(qū)大體對應(yīng)了累積降水的大值區(qū)(尤其強(qiáng)降水區(qū))。盡管圖10 中,不同方案的Mlt、CLcr 及Evap 項(xiàng)變化趨勢較為相似,差別主要體現(xiàn)在量值上,但圖12 顯示,上述三個(gè)源匯項(xiàng)的空間分布在不同方案中也存在較大差異,CLcr 項(xiàng)對雨滴的增長貢獻(xiàn)明顯強(qiáng)于Mlt 項(xiàng),區(qū)域極值也更大。Morrison 方案中,CLcr 及Mlt 項(xiàng)呈水平分布,因此最強(qiáng)降水區(qū)及最強(qiáng)蒸發(fā)區(qū)也呈水平分布。Thompson 方案中,CLcr 和Mlt 項(xiàng)呈“西北—東南”走向,因此最強(qiáng)降水區(qū)及最強(qiáng)蒸發(fā)區(qū)也呈“西北—東南”分布。MY 方案中的雨帶分部相對復(fù)雜,一條從鄂西北開始延伸,另一條在湖北南部,CLcr、Mlt 和Evap 項(xiàng)分別與其有著良好對應(yīng)關(guān)系。這些差異導(dǎo)致不同方案中模擬降水的位置和強(qiáng)度差異。冰相粒子融化或雨滴搜集云滴的高估可能是造成區(qū)域平均降水模擬值偏強(qiáng)的重要原因之一。

    圖12 2018 年6 月29 日22:00 至30 日11:00 Morrison 方案(第一行)、Thompson 方案(第二行)以及MY 方案(第三行)中源項(xiàng)Mlt(左列)、CLcr(中間列)及匯項(xiàng)Evap(右列)的高度平均分布,單位:g m-3Fig. 12 Distribution of the height-averaged total amount of Mlt (left column), CLcr (middle column) source terms, and (c) Evap sink term (right column) of rain during 2200 UTC 29 June to 1100 UTC 30 June 2018 by Morrison (top line), (b) Thompson (second line), and MY (bottom line)schemes, units: g m-3

    需要注意的是,Evap 項(xiàng)的空間分布并非和Mlt 或CLcr 項(xiàng)一一對應(yīng),因?yàn)槠渲饕獙?yīng)的是雨滴含量,而影響雨滴含量分布的因子除了微物理過程項(xiàng)之外,還有平流項(xiàng)、湍流混合項(xiàng)等(主要針對小粒子)。

    云滴是影響降水的另一個(gè)重要因子。圖13 為云滴的三個(gè)關(guān)鍵源匯項(xiàng)區(qū)域平均的垂直分布圖。從圖13 可知,三個(gè)方案中的主要云滴源匯項(xiàng)的垂直分布存在相似之處,云滴凝結(jié)的峰值高度均分別位于在750 hPa 及550 hPa 左右高度,雨滴搜集云滴的峰值高度也大致相似,位于600 hPa 高度附近。但三者也存在明顯差異,Thompson 中(圖13b),水汽凝結(jié)明顯強(qiáng)于其它兩個(gè)方案,峰值為9.24 g m-3,次峰值為5.8 g m-3;而Morrison 方案中(圖13a),水汽凝結(jié)峰值為4.94 g m-3,次峰值為3.67 g m-3;MY 方案中(圖13c),盡管水汽凝結(jié)峰值略小于Morrison 方案,為4.82 g m-3,但次峰值明顯大于Morrison 方案,為4.16 g m-3。除了凝結(jié)總量之外,三個(gè)方案中水汽凝結(jié)在高層也存在明顯差異。MY方案中,300 hPa 以上水汽凝結(jié)明顯,最大凝結(jié)高度高于250 hPa,Morrison 方案中的最大凝結(jié)高度略低于250 hPa,而Thompson 方案中的最大凝結(jié)高度略高于300 hPa,且300 hPa 以上高度凝結(jié)量極少。該現(xiàn)象是圖9 中云滴能夠在400 hPa 高度以上仍存在大值中心的重要原因??傮w而言,MY 方案中的凝結(jié)效應(yīng)略強(qiáng)于Morrison 方案。

    圖13 不同方案的云滴主要源匯項(xiàng)總量(單位:g m-3))區(qū)域平均值的的垂直分布:(a)Morrison 方案;(b)Thompson 方案;(c)MY方案Fig. 13 Vertical profiles of area-averaged content (unit: g m-3) of the key sink and source termsof cloud droplet by (a) Morrison, (b) Thompson, and(c) MY schemes

    由于Thompson 方案中水汽凝結(jié)效應(yīng)最強(qiáng),因此,該方案中雨滴搜集云滴的總量顯著多于其它兩個(gè)方案,且峰值高度相近。盡管Morrison 方案中的水汽凝結(jié)略弱于MY 方案,但雨滴搜集云滴的量仍明顯高于MY 方案。對于冰相粒子搜集云滴來說(Col2 包含冰相粒子直接搜集云滴,與Col1類似,只是專門針對云滴而言),Thompson 方案依然最強(qiáng),MY 方案其次,Morrison 方案略微弱于Morrison 方案。三者的極值高度存在較大差異。Thompson 方案中,Col2 的極值高度最低,略高于500 hPa 高度;Morrison 方案中,Col2 極值高度在500 hPa 附近;MY 方案中,Col2 極值高度位置僅略低于400 hPa。

    由以上分析可知,Thompson 方案中,凝結(jié)效應(yīng)、雨滴搜集云滴增長及云滴對冰粒子增長的貢獻(xiàn)均最強(qiáng);MY 方案中的凝結(jié)效應(yīng)及云滴對冰相粒子增長的貢獻(xiàn)強(qiáng)于Morrison 方案,但后者雨滴搜集云滴增長強(qiáng)于前者,兩者凝結(jié)效應(yīng)及雨滴對云滴搜集增長的峰值高度相近,但云滴對冰相粒子貢獻(xiàn)的極值高度則存在明顯差異。

    4.4 冰相粒子主要源匯項(xiàng)特征分析

    圖14 為冰相粒子的主要源匯項(xiàng)平均值時(shí)間變化的比較。從前文可知,就本次過程而言,冰相粒子的融化是影響雨滴增長的關(guān)鍵因素之一。因此,首先將不同方案的冰相粒子都作為一個(gè)整體進(jìn)行比較。在輸出冰相粒子源匯項(xiàng)之后發(fā)現(xiàn),冰相粒子融化、水汽凝華及冰相粒子與雨滴和云滴的碰并過程是冰相粒子增長的重要來源。下文將分別對其進(jìn)行比較。就Mlt 項(xiàng)來說,Morrison 方案模擬含量最多,MY 和Thompson 方案差別不大(圖14a)。就VD和Col1 項(xiàng)而言,Thompson 方案模擬量最多,其次是Morrison 方案,MY 方案模擬相對最少(圖14b和c)。三個(gè)方案中的Col2 值大于Col1(圖14d),一個(gè)重要原因在于,水汽在較高處凝結(jié)成云滴(圖13),然后云滴和冰相粒子發(fā)生碰并過程。而雨滴與冰相粒子的碰并過程發(fā)生在相對低的高度,盡管雨滴含量相對較大,但參與凇附過程幾率相對較小,故三個(gè)方案中的Col2 項(xiàng)均強(qiáng)于Col1 項(xiàng)。尤其是MY 方案,Col2 值數(shù)倍于Col1 值,而Thompson方案中兩者差值相對較小。

    圖14 2018 年6 月29 日22:00 至30 日11:00 不同方案冰相粒子主要(a)匯項(xiàng)Mlt 及源項(xiàng)(b)VD、(c)Col1、(d)Col2 空間平均值(單位:g m-3)的時(shí)間變化Fig. 14 Time series of the volumetrically averaged content (units: g m-3) of the distribution of the height-averaged total amount of the (a) sink term Mlt and source term (b) VD, (c) Col1, (d) Col2 of rain during 2200 UTC 29 June to 1100 UTC 30 June 2018 by different schemes

    圖15 為Col2 和VD 項(xiàng)的垂直分布(冰相粒子融化項(xiàng)Mlt 及冰相粒子對雨滴的碰并搜集Col1 已在圖11 中討論,不再贅述)。從圖15a 和b 可知,Morrison 和Thompson 方案中,Col2 的極值高度差別不大,略低于MY 方案。Morrison 和MY 方案中的Col2 極值差別不大,約為Thompson 方案極值的一半。對VD 來說,Morrison 和Thompson 方案的形態(tài)相似,極值高度在350 hPa 附近,但前者數(shù)值約為后者一半。MY 方案中,VD 的分布形態(tài)表現(xiàn)出較大差異,廓線明顯窄于其它兩個(gè)方案,并且極值高度在250 hPa 附近,顯著高于其它兩個(gè)方案。需要說明的是,MY 方案中,凝華增長效應(yīng)在中底層為負(fù)值,這是因?yàn)樵摳叨葏^(qū)間的區(qū)域平均表現(xiàn)為升華。由圖15b 可知,Thompson 方案中的升華效應(yīng)最強(qiáng)。事實(shí)上,該方案的凝華效應(yīng)也表現(xiàn)為最強(qiáng)。Morrison 和Thompson 方案中,Col2 的極值高度低于VD 的極值高度,MY 方案則恰好相反。

    圖15 同圖11,但為云滴主要源匯項(xiàng)(CLcr, Col2, Cond)Fig. 15 Same as Fig.11 but for key sink and source of the cloud droplet (CLcr, Col2, and Cond)

    事實(shí)上,三個(gè)方案中Col2 和VD 的差異揭示了導(dǎo)致冰相粒子分布差異的重要原因之一(其它兩個(gè)原因在前面已經(jīng)描述)。下文將從不同冰相粒子的角度進(jìn)行分析。

    圖16 為冰晶、雪和霰粒子的凝華/升華效應(yīng)。從圖16 可知,三個(gè)方案的凝華/升華效應(yīng)存在顯著差異。Morrison 方案中,雪粒子凝華和升華效應(yīng)均強(qiáng)于霰粒子,冰晶粒子的凝華效應(yīng)較小,升華效應(yīng)幾乎可以忽略不計(jì)(圖16a)。Thompson 方案中,雪粒子的升華效應(yīng)占主導(dǎo)地位,其它兩類冰相粒子的貢獻(xiàn)與其相比微乎其微(圖16b)。MY 方案中,冰晶的凝華效應(yīng)略弱于雪粒子,從其廓線形態(tài)來看,冰晶粒子的凝華差異是導(dǎo)致圖16 中冰相粒子與其它兩個(gè)方案顯著差異的重要原因,也是圖9 中MY 方案冰相粒子明顯多于其它方案的關(guān)鍵原因。

    圖16 不同方案的冰晶、雪和霰粒子的凝華/升華項(xiàng)總量(單位:g m-3)區(qū)域平均值的的垂直分布:(a)Morrison 方案;(b)Thompson方案;(c)MY 方案Fig. 16 Vertical profiles of area-averaged content (units: g m-3) of the the deposition/sublimation of ice, snow, and graupel by (a) Morrison,(b) Thompson, and (c) MY schemes

    圖17 為不同云微物理參數(shù)化方案中冰相粒子同云滴、雨滴之間的凇附(碰并)過程。從圖14可知,不同云微物理參數(shù)化方案中冰相粒子與云滴和雨滴之間的碰并過程存在較大差異。圖17 試圖對產(chǎn)生該差異的原因進(jìn)行初步解釋。Morrison 和MY 中,CLcg 占主導(dǎo)地位(圖14a 和c),MY 方案中其它碰并過程相對CLcg 項(xiàng)來說幾乎可以忽略不計(jì),而Morrison 方案中,其它凇附過程的貢獻(xiàn)相對MY 方案來說,依然較為明顯。Thompson 方案與其它兩個(gè)方案有著顯著區(qū)別(圖14c):CLcs貢獻(xiàn)最大,并且,CLrg 強(qiáng)于CLrs 項(xiàng)。

    圖17 2018 年6 月29 日22:00 至30 日11:00(a)Morrison、(b)Thompson、(c)MY 方案主要凇附項(xiàng)空間平均值(單位:g m-3)的時(shí)間變化Fig. 17 Time series of the volumetrically averaged content (units: g m-3) of riming terms during 2200 UTC 29 June to 1100 UTC 30 June 2018 by (a)Morrison, (b) Thompson, (c) MY schemes

    圖18 為雪粒子(霰粒子)搜集云滴(雨滴)增長總量區(qū)域平均的垂直廓線。從圖18 所示的雪粒子(霰粒子)搜集雨滴(云滴)增長廓線圖可知,以上四個(gè)微物理過程在三個(gè)方案中存在明顯差異。Morrison 方案中(圖18a),CLcg 效應(yīng)最明顯,但CLrg 的極值略大于CLcg,CLcs 效應(yīng)強(qiáng)于CLrs,但前者極值小于后者。從圖18a 可知,霰粒子搜集液態(tài)粒子的增長大于雪粒子與液態(tài)粒子的碰并增長,這是圖9a2 和圖9a3 所示霰粒子方案中霰粒子含量大于雪粒子的重要原因之一。Thompson 方案中(圖18b),CLcs 效應(yīng)顯著強(qiáng)于其它幾個(gè)過程,并且CLrg 存在負(fù)值區(qū)。如前所述,這是因?yàn)門hompson 方案處理雨滴與冰相粒子的碰并時(shí)與其它方案存在差異,該方案在零度層以上,兩者碰并將有助于冰相粒子增長,而零度層以下,兩者碰并將貢獻(xiàn)于雨滴增長。從CLrg 的垂直廓線來看,霰粒子因搜集雨滴的增長不明顯,而零度層以下兩者的碰并主要貢獻(xiàn)于雨滴增長。這是圖18b3 中霰粒子較少的重要原因之一。MY 方案中(圖18c),CLcs、CLrs 和CLrg 的量值明顯小于CLCg(存在量級差別)。圖18c2 和圖18c3 中,霰粒子大于雪粒子的含量。霰粒子對云滴的搜集作用對其有著重要貢獻(xiàn)。

    圖18 同圖11,但為雪粒子(霰粒子)搜集雨滴(云滴)的增長(CLcs、CLrs、CLcg 和CLrg)Fig. 18 Same as Fig.11, but for the accretion (CLcs, CLrs, CLcg, and CLrg) of raindrops (cloud droplets) by the snow (graupel)

    5 結(jié)論與討論

    本文通過WRF(3.4.1)模式,針對2018 年6月29~30 日一次梅雨鋒背景下的暴雨過程進(jìn)行數(shù)值模擬,分別采用Morrison、Thompson 和MY 方案進(jìn)行對比分析。結(jié)果發(fā)現(xiàn),三個(gè)方案模擬的背景場在天氣尺度上都與再分析資料一致,能夠再現(xiàn)出有利于強(qiáng)降水發(fā)生的環(huán)流場。但不同云微物理過程對梅雨期暴雨的局地環(huán)流有著顯著影響,且不同方案存在顯著的差異。

    Thompson 方案產(chǎn)生更強(qiáng)的局地環(huán)流系統(tǒng)的變率和上升氣流。由于天氣尺度環(huán)流的一致性,這些差異直接與云物理參數(shù)化方案的差異有關(guān),同時(shí)也與因云微物理方案差異而導(dǎo)致的其它過程有關(guān)。這些來自多個(gè)方面的相互作用可能是非線性的,屬于云微物理參數(shù)化方案的間接影響。三個(gè)方案的模擬降水均有所夸大,小時(shí)降水率、小時(shí)累計(jì)降水量均始終大于觀測值。總體來看,Morrison 方案的模擬效果相對最優(yōu)。各方案模擬的水成物存在較大差異。Thompson 方案中雪粒子含量明顯偏多,冰晶和霰相對較少;MY 方案中冰晶粒子含量最多,其次為Morrison 方案,兩者的霰粒子含量均多于雪粒子。液態(tài)粒子(尤其云滴)的分布及演變均存在較大差異。零度層以上,MY 方案中的云滴含量最多,其次是Morrison 方案,Thompson 方案相對最少,而不同方案的雨滴在該高度區(qū)間的差異相對較小,即過冷水的差異(針對區(qū)域平均而言)主要體現(xiàn)在過冷云滴上。

    冰相粒子融化、雨滴對云滴搜集作用是雨滴的關(guān)鍵源項(xiàng),雨滴蒸發(fā)是最重要的匯項(xiàng),以上三個(gè)物理量在不同云微物理方案中存在空間上的差異,這是導(dǎo)致降水空間分布差異的直接原因??偟膩碚f,雨滴對云滴搜集效應(yīng)大于冰相粒子融化。冰相粒子對雨滴的搜集效應(yīng)相對較弱,但不同方案差異大。Thompson 方案中的雨滴蒸發(fā)項(xiàng)顯著大于其它兩個(gè)方案,在底層表現(xiàn)得最為明顯。Morrison 方案中,冰相粒子在零度層以下高度并不直接轉(zhuǎn)換為雨滴,而是用于提高冰相粒子融化率并增加參與碰并雨滴的脫落(shedding);MY 方案中,兩者碰并僅發(fā)生在零度層以上高度。凇附過程的主要影響因子包括碰并粒子的數(shù)濃度、直徑、下落末速度等。與此同時(shí),過冷水的時(shí)空分布差異也將對其產(chǎn)生影響。對不同方案來說,凇附相關(guān)的云微物理過程存在較大差異。Morrison 方案中雨滴貢獻(xiàn)明顯大于MY方案。而Thompson 方案中,零度層以上,僅少量雨滴作為冰相粒子的增長源;零度層以下,冰相粒子在碰并過程中直接轉(zhuǎn)換為雨滴。MY 方案中,過冷云滴能到達(dá)的最高處高于其它兩個(gè)方案,因此凇附過程發(fā)生在相對較高處(圖18)。并且,相對于其它淞附過程來說,零度層以上豐富的云滴和霰粒子之間的碰凍絕對占優(yōu)。而Thompson 方案中,雪粒子含量最多,其與過冷云滴之間的碰凍明顯大于其它淞附過程。由于凇附過程是冰相粒子的關(guān)鍵來源之一,故不同方案中過冷水的差異最終將影響冰相粒子的分布。同時(shí),這種差異將通過其它云微物理過程及熱動(dòng)力過程進(jìn)行反饋。

    Thompson 方案中,由于更多冰相粒子融化成雨滴,且零度層以下的水汽凝結(jié)效應(yīng)最強(qiáng),由于更多云滴參與碰并,因此雨滴搜集更多云滴,使得該方案模擬的雨滴含量最多,降水最強(qiáng)。該方案中水汽主要凝華為雪,且其在與過冷水碰并增長過程中占主導(dǎo)地位,因此該方案中模擬的雪最多。并且,從圖18 可知,零度層以下高度,該方案中被雨滴搜集的冰相粒子為霰,且零度層以上高度,極少量霰通過凇附方式增長。這也是Thompson 方案中霰粒子含量少于雪粒子的原因之一。盡管MY 方案中水汽凝結(jié)量多于Morrison 方案,但由于后者冰相粒子融化效應(yīng)更強(qiáng),生成的雨滴更多,導(dǎo)致雨滴搜集了更多云滴。Morrison 和MY 方案中,霰粒子含量均大于雪粒子,但兩者形成的關(guān)鍵原因并不同。Morrison 方案中,凝華形成的霰粒子含量大于雪粒子,且霰的凇附效應(yīng)強(qiáng)于雪。而MY 方案中,凝華形成的雪粒子含量大于霰,但霰粒子的凇附增長起主導(dǎo)地位,且大于凝華效應(yīng),從而使得霰粒子含量大于雪粒子。該方案中,水汽凝華在約300 hPa 以上高度形成大量冰晶,顯著強(qiáng)于其它方案,這是該方案中高層冰晶粒子明顯大于其它方案的重要原因之一。

    從以上討論可以發(fā)現(xiàn),云微物理過程通過不同方式影響降水。但是,對降水來說,最直接的影響還是雨滴源匯項(xiàng)本身。由于冰相粒子融化及雨滴搜集云滴增長的大值區(qū)與累積降水分布的大值區(qū)相匹配。因此,冰相粒子融化或雨滴搜集云滴效應(yīng)的高估,可能是造成區(qū)域平均降水模擬值偏強(qiáng)的重要原因之一。僅從降水比較來看,Morrison 方案在本次模擬過程中的結(jié)果相對最優(yōu)。由于缺少云微物理觀測資料,暫無法從微觀方面對三個(gè)方案進(jìn)行比較。

    需要說明的是,冰晶向雪的轉(zhuǎn)換也是重要的云微物理過程之一。但在凝華和凇附過程中,冰相粒子變化量數(shù)倍于該過程(圖略),并且,本文側(cè)重于將冰相粒子作為一個(gè)整體進(jìn)行考慮,故暫不討論冰相粒子之間轉(zhuǎn)換相關(guān)的云微物理量。

    不同云微物理參數(shù)化方案中,影響模擬結(jié)果的重要因素之一就是冰相粒子的處理(Morrison and Milbrandt, 2015; Zhou et al., 2020)。因?yàn)閷ε七^程而言,云滴和雨滴均假設(shè)為球體,不同方案之間差異并不大。而在冷云過程中,冰相粒子的定義(形狀、下落末速度、密度等)在不同方案中存在一定差異,這將導(dǎo)致冰相粒子的模擬結(jié)果存在差異。這種差異除了體現(xiàn)在地表累積降水中,更體現(xiàn)在冰相粒子的空間分布上。

    云微物理參數(shù)化方案的敏感性研究工作,除了了解不同方案之間的差異之外,最關(guān)鍵的是希望最終能夠在此基礎(chǔ)上進(jìn)行改進(jìn)。本文的模擬結(jié)果表明,三個(gè)方案的模擬降水相對于觀測來說有所夸大。僅從云微物理過程的影響來說,改進(jìn)參數(shù)化方案可以考慮從減小融化量、雨滴搜集云滴量、或者增加蒸發(fā)量等等方面開展工作。但是,熱動(dòng)力方面的影響同樣不可忽視。實(shí)際上,云微物理過程對背景場的熱動(dòng)力過程有著重要反饋?zhàn)饔茫▽O建華和趙思雄,2000),云微物理參數(shù)化差異會導(dǎo)致背景場的差異。加熱場的異??梢杂绊懢值匚粶u的分布(Ertel,1942; Hoskins et al., 1985);與不同微物理過程有關(guān)的異常加熱場的疊加更可以形成復(fù)雜的位渦分布,并且隨著氣流的上升和平流作用甚至可以影響高層和下游的環(huán)流變化(Joos and Wernli, 2012; Joos and Forbes, 2016; Crezee et al., 2017; Attinger et al.,2019),該變化又將對云微物理過程有著重要影響。

    因此,為了改進(jìn)云微物理參數(shù)化,除了深入了解不同云微物理參數(shù)化的差異,做更多敏感性工作之外,還必須跟蹤云微物理過程與熱動(dòng)力過程之間的相互作用,以便更精準(zhǔn)地調(diào)整相關(guān)云微物理過程參數(shù)化。同時(shí),外場試驗(yàn)的開展也非常有必要。對梅雨鋒暴雨而言,從宏觀和微觀角度對模式熱動(dòng)力場及微物理過程進(jìn)行驗(yàn)證、約束和修改有著重要意義。

    本文屬個(gè)例研究,更多的個(gè)例差異、更詳細(xì)的影響機(jī)制還有待于后續(xù)進(jìn)一步研究。

    猜你喜歡
    云滴雨滴極值
    極值點(diǎn)帶你去“漂移”
    福建九仙山一次佛光的云滴粒子尺度分析*
    極值點(diǎn)偏移攔路,三法可取
    小雨滴
    草原歌聲(2020年3期)2020-12-08 22:22:28
    可愛的小雨滴
    你不知道的雷雨知識
    一類“極值點(diǎn)偏移”問題的解法與反思
    小雨滴
    云微物理特性及云滴有效半徑參數(shù)化:一次降水層狀云的飛機(jī)觀測資料結(jié)果
    小雨滴
    兒童繪本(2015年7期)2015-05-25 17:55:30
    在线天堂中文资源库| 国产精品不卡视频一区二区| 晚上一个人看的免费电影| 最后的刺客免费高清国语| 亚洲国产精品一区三区| 最近最新中文字幕免费大全7| 香蕉国产在线看| 久久99一区二区三区| 欧美最新免费一区二区三区| 18禁国产床啪视频网站| 秋霞伦理黄片| 高清在线视频一区二区三区| 久久久精品区二区三区| 97人妻天天添夜夜摸| 高清欧美精品videossex| 卡戴珊不雅视频在线播放| 欧美另类一区| 日本猛色少妇xxxxx猛交久久| av片东京热男人的天堂| 五月伊人婷婷丁香| 国产黄色免费在线视频| 永久免费av网站大全| 国产成人精品在线电影| av有码第一页| 久久av网站| 亚洲国产精品999| 97超碰精品成人国产| 国产爽快片一区二区三区| 亚洲成人一二三区av| 韩国av在线不卡| 亚洲国产精品成人久久小说| 国产在视频线精品| 精品卡一卡二卡四卡免费| 国产男人的电影天堂91| 国产一区亚洲一区在线观看| 亚洲欧洲国产日韩| 国产熟女午夜一区二区三区| 大香蕉久久成人网| 亚洲国产精品一区二区三区在线| av.在线天堂| 成人亚洲欧美一区二区av| 黄色一级大片看看| 看免费成人av毛片| 深夜精品福利| 99re6热这里在线精品视频| 中国美白少妇内射xxxbb| 一区二区日韩欧美中文字幕 | 亚洲,欧美精品.| 18+在线观看网站| 国产av一区二区精品久久| 精品一品国产午夜福利视频| 国产高清三级在线| 免费黄网站久久成人精品| 在线观看www视频免费| av天堂久久9| 国国产精品蜜臀av免费| 自拍欧美九色日韩亚洲蝌蚪91| av.在线天堂| 精品少妇黑人巨大在线播放| 日韩制服丝袜自拍偷拍| 久久精品国产亚洲av天美| 久久久久久久大尺度免费视频| 亚洲五月色婷婷综合| 精品一区二区三区视频在线| 夜夜骑夜夜射夜夜干| 最后的刺客免费高清国语| 亚洲综合精品二区| 亚洲精品国产av成人精品| 美国免费a级毛片| 丰满少妇做爰视频| 咕卡用的链子| 亚洲综合精品二区| av卡一久久| 十分钟在线观看高清视频www| 亚洲欧美成人综合另类久久久| 中文字幕亚洲精品专区| 午夜影院在线不卡| 欧美人与善性xxx| 亚洲欧美精品自产自拍| 亚洲精品,欧美精品| 爱豆传媒免费全集在线观看| 亚洲情色 制服丝袜| 精品卡一卡二卡四卡免费| 欧美xxxx性猛交bbbb| 成年av动漫网址| 国产成人午夜福利电影在线观看| 久久人妻熟女aⅴ| 亚洲欧美精品自产自拍| 香蕉丝袜av| 国产av码专区亚洲av| 在线天堂最新版资源| 免费观看在线日韩| 在线天堂最新版资源| 高清黄色对白视频在线免费看| 免费观看在线日韩| 边亲边吃奶的免费视频| 男女啪啪激烈高潮av片| 婷婷色麻豆天堂久久| 欧美日韩视频精品一区| 欧美日韩视频精品一区| 国产一区二区在线观看日韩| 国产国语露脸激情在线看| 欧美日韩精品成人综合77777| 精品国产乱码久久久久久小说| 99香蕉大伊视频| 视频区图区小说| 女人精品久久久久毛片| 精品国产一区二区久久| 美女国产视频在线观看| av片东京热男人的天堂| 我的女老师完整版在线观看| 中文精品一卡2卡3卡4更新| 中文精品一卡2卡3卡4更新| 国产国拍精品亚洲av在线观看| 久久久精品94久久精品| 在线观看国产h片| 亚洲情色 制服丝袜| 啦啦啦视频在线资源免费观看| 国产成人a∨麻豆精品| 人妻人人澡人人爽人人| a级毛片黄视频| 最近的中文字幕免费完整| 在线观看www视频免费| 久久久久久久国产电影| 极品少妇高潮喷水抽搐| 高清毛片免费看| 新久久久久国产一级毛片| 巨乳人妻的诱惑在线观看| 精品一品国产午夜福利视频| 蜜桃在线观看..| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久久久精品精品| 丰满少妇做爰视频| 在线 av 中文字幕| 99热全是精品| 韩国av在线不卡| 青春草国产在线视频| 免费看av在线观看网站| 韩国精品一区二区三区 | 国产国语露脸激情在线看| 99热网站在线观看| 在线观看三级黄色| 久久影院123| 18在线观看网站| 18在线观看网站| 看非洲黑人一级黄片| 桃花免费在线播放| 亚洲欧洲日产国产| 精品第一国产精品| 一二三四中文在线观看免费高清| 国产免费一区二区三区四区乱码| 国产精品嫩草影院av在线观看| a 毛片基地| 国产精品久久久久久久久免| 亚洲少妇的诱惑av| 蜜桃在线观看..| 性高湖久久久久久久久免费观看| 最近中文字幕2019免费版| 永久免费av网站大全| 成人手机av| 成人国产麻豆网| 久久人人爽人人爽人人片va| 欧美性感艳星| 91精品伊人久久大香线蕉| 我要看黄色一级片免费的| 国产成人aa在线观看| 亚洲av综合色区一区| 亚洲av在线观看美女高潮| 巨乳人妻的诱惑在线观看| 久久久欧美国产精品| 亚洲久久久国产精品| 女人久久www免费人成看片| 精品国产一区二区三区久久久樱花| 男女免费视频国产| 三上悠亚av全集在线观看| 国产不卡av网站在线观看| 久久精品国产亚洲av涩爱| 99国产综合亚洲精品| 纵有疾风起免费观看全集完整版| 欧美精品一区二区免费开放| 成人国产麻豆网| 亚洲一级一片aⅴ在线观看| 国产精品女同一区二区软件| 好男人视频免费观看在线| 中文乱码字字幕精品一区二区三区| 亚洲欧美日韩另类电影网站| 在线观看www视频免费| 高清毛片免费看| 国产xxxxx性猛交| 欧美激情国产日韩精品一区| 在线天堂最新版资源| 这个男人来自地球电影免费观看 | 免费黄色在线免费观看| 黑人巨大精品欧美一区二区蜜桃 | av福利片在线| 大码成人一级视频| 久久97久久精品| 中文字幕精品免费在线观看视频 | 最近中文字幕2019免费版| 日韩人妻精品一区2区三区| 日韩不卡一区二区三区视频在线| 精品久久国产蜜桃| 看免费av毛片| 亚洲四区av| 国产成人欧美| 亚洲国产看品久久| 久久鲁丝午夜福利片| 欧美xxⅹ黑人| 日韩精品有码人妻一区| 亚洲内射少妇av| 久久韩国三级中文字幕| 老司机亚洲免费影院| 国产爽快片一区二区三区| 久久久欧美国产精品| 国产熟女欧美一区二区| 久久99精品国语久久久| 成人国语在线视频| 中国国产av一级| 久久精品国产综合久久久 | 久久久久久久久久人人人人人人| 人成视频在线观看免费观看| 日本午夜av视频| 中国三级夫妇交换| 黄片播放在线免费| 又黄又粗又硬又大视频| 国产一级毛片在线| 高清在线视频一区二区三区| 亚洲国产精品999| 国产av码专区亚洲av| 国产精品无大码| 久久久久久伊人网av| 欧美少妇被猛烈插入视频| 成人无遮挡网站| 久久人人爽人人片av| 三上悠亚av全集在线观看| 欧美变态另类bdsm刘玥| 亚洲成人av在线免费| 男的添女的下面高潮视频| 精品国产一区二区久久| 制服人妻中文乱码| 午夜老司机福利剧场| 欧美3d第一页| 国产xxxxx性猛交| 成人免费观看视频高清| a级毛片在线看网站| 少妇人妻久久综合中文| 成人漫画全彩无遮挡| 老司机亚洲免费影院| 18禁观看日本| 欧美激情 高清一区二区三区| 插逼视频在线观看| 91国产中文字幕| 亚洲综合色网址| 大片免费播放器 马上看| 亚洲国产精品成人久久小说| 久久精品国产鲁丝片午夜精品| 国产成人欧美| 乱码一卡2卡4卡精品| 国产深夜福利视频在线观看| 欧美变态另类bdsm刘玥| 九九在线视频观看精品| 日本与韩国留学比较| 2018国产大陆天天弄谢| av卡一久久| 秋霞在线观看毛片| 青春草视频在线免费观看| 亚洲经典国产精华液单| 中文欧美无线码| 久久久久视频综合| 国精品久久久久久国模美| 免费黄色在线免费观看| 国产精品麻豆人妻色哟哟久久| 国产探花极品一区二区| 久久99精品国语久久久| xxx大片免费视频| 人妻人人澡人人爽人人| 免费黄网站久久成人精品| kizo精华| 亚洲一区二区三区欧美精品| 一级毛片 在线播放| 看免费av毛片| 国产精品一二三区在线看| 一级爰片在线观看| 久久久国产精品麻豆| 久久99精品国语久久久| 国产亚洲午夜精品一区二区久久| 我要看黄色一级片免费的| 免费不卡的大黄色大毛片视频在线观看| 久久综合国产亚洲精品| 我要看黄色一级片免费的| 精品人妻一区二区三区麻豆| 久久久久国产网址| 一区二区三区精品91| 国产视频首页在线观看| av电影中文网址| 免费黄网站久久成人精品| 亚洲国产精品999| 1024视频免费在线观看| 国产国拍精品亚洲av在线观看| 久久 成人 亚洲| 久热久热在线精品观看| 一区二区日韩欧美中文字幕 | 国产一区二区三区av在线| 国产高清国产精品国产三级| 麻豆精品久久久久久蜜桃| 巨乳人妻的诱惑在线观看| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三区四区五区乱码 | 国产精品偷伦视频观看了| 欧美精品av麻豆av| 国产 精品1| 亚洲三级黄色毛片| 亚洲第一区二区三区不卡| 七月丁香在线播放| 亚洲国产最新在线播放| 韩国精品一区二区三区 | 久久精品久久精品一区二区三区| 51国产日韩欧美| 精品视频人人做人人爽| 男女下面插进去视频免费观看 | 国产精品久久久久久精品古装| 国产精品国产av在线观看| 寂寞人妻少妇视频99o| 一本色道久久久久久精品综合| 夜夜骑夜夜射夜夜干| 内地一区二区视频在线| 国产在线视频一区二区| 这个男人来自地球电影免费观看 | 久久婷婷青草| 中文字幕人妻丝袜制服| 色婷婷av一区二区三区视频| 国产1区2区3区精品| 国产欧美日韩一区二区三区在线| 人妻一区二区av| 国产精品 国内视频| 免费av中文字幕在线| 欧美 亚洲 国产 日韩一| 国产亚洲av片在线观看秒播厂| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到 | 精品第一国产精品| 狂野欧美激情性bbbbbb| 80岁老熟妇乱子伦牲交| 欧美精品一区二区大全| √禁漫天堂资源中文www| 成年人午夜在线观看视频| 国产淫语在线视频| 精品福利永久在线观看| 少妇的丰满在线观看| 大香蕉97超碰在线| 日韩电影二区| 亚洲精品视频女| 久久ye,这里只有精品| 侵犯人妻中文字幕一二三四区| 国产av一区二区精品久久| 日产精品乱码卡一卡2卡三| 国产日韩一区二区三区精品不卡| 韩国高清视频一区二区三区| 高清视频免费观看一区二区| 日本与韩国留学比较| av国产精品久久久久影院| 国产深夜福利视频在线观看| 久久女婷五月综合色啪小说| 亚洲av.av天堂| 2021少妇久久久久久久久久久| 亚洲第一av免费看| 伊人亚洲综合成人网| 黑人猛操日本美女一级片| 久久精品久久久久久久性| 高清视频免费观看一区二区| 日本av免费视频播放| 亚洲精品久久午夜乱码| 精品99又大又爽又粗少妇毛片| 成人免费观看视频高清| 精品国产乱码久久久久久小说| av卡一久久| 丁香六月天网| 午夜福利在线观看免费完整高清在| 高清av免费在线| 日韩 亚洲 欧美在线| 久久av网站| 中文字幕人妻熟女乱码| 国产成人a∨麻豆精品| av在线app专区| 午夜福利在线观看免费完整高清在| 日韩av免费高清视频| 在线观看www视频免费| 大香蕉久久网| 日韩大片免费观看网站| 97精品久久久久久久久久精品| 91精品三级在线观看| 久久精品久久久久久噜噜老黄| 成人无遮挡网站| 久久狼人影院| 国产精品偷伦视频观看了| 日日摸夜夜添夜夜爱| 新久久久久国产一级毛片| 亚洲欧美精品自产自拍| 一级a做视频免费观看| 热re99久久精品国产66热6| 亚洲中文av在线| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美清纯卡通| 成人综合一区亚洲| 这个男人来自地球电影免费观看 | 国产又色又爽无遮挡免| 亚洲色图 男人天堂 中文字幕 | 女人被躁到高潮嗷嗷叫费观| 人妻系列 视频| 久久精品国产亚洲av涩爱| 少妇被粗大猛烈的视频| 国产精品久久久久成人av| 亚洲精品久久成人aⅴ小说| 在线观看免费日韩欧美大片| 高清不卡的av网站| 丁香六月天网| 午夜福利网站1000一区二区三区| 国产亚洲最大av| 久久精品久久久久久久性| 日本黄色日本黄色录像| 国产极品粉嫩免费观看在线| 亚洲国产最新在线播放| 美女中出高潮动态图| a级毛片在线看网站| 亚洲 欧美一区二区三区| 黄网站色视频无遮挡免费观看| 国产成人精品福利久久| 美女国产视频在线观看| 水蜜桃什么品种好| 少妇被粗大的猛进出69影院 | 激情五月婷婷亚洲| 蜜桃国产av成人99| 精品国产露脸久久av麻豆| 亚洲内射少妇av| 国产精品久久久久久久久免| 在线天堂最新版资源| 国产男女内射视频| a级毛色黄片| 亚洲经典国产精华液单| 欧美另类一区| 新久久久久国产一级毛片| 22中文网久久字幕| 91精品伊人久久大香线蕉| 秋霞伦理黄片| av播播在线观看一区| 精品卡一卡二卡四卡免费| 一区二区日韩欧美中文字幕 | 搡老乐熟女国产| 老熟女久久久| 大话2 男鬼变身卡| 美女xxoo啪啪120秒动态图| 国产高清不卡午夜福利| 99热国产这里只有精品6| 国产精品久久久久成人av| 欧美+日韩+精品| 亚洲成国产人片在线观看| 韩国av在线不卡| 黑丝袜美女国产一区| 中国国产av一级| 精品人妻一区二区三区麻豆| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| av视频免费观看在线观看| 亚洲av成人精品一二三区| 免费黄色在线免费观看| 国产白丝娇喘喷水9色精品| 欧美国产精品一级二级三级| 久久精品国产a三级三级三级| 22中文网久久字幕| 亚洲高清免费不卡视频| 中国美白少妇内射xxxbb| 成人国产av品久久久| 精品卡一卡二卡四卡免费| 考比视频在线观看| 久久久久精品久久久久真实原创| 最近的中文字幕免费完整| 99精国产麻豆久久婷婷| 国产精品女同一区二区软件| 国产色婷婷99| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 国产av精品麻豆| 天堂8中文在线网| 巨乳人妻的诱惑在线观看| 成人国产麻豆网| 少妇 在线观看| av.在线天堂| 青青草视频在线视频观看| 青春草视频在线免费观看| 99热国产这里只有精品6| 欧美xxxx性猛交bbbb| 全区人妻精品视频| 欧美日韩视频高清一区二区三区二| 久久精品熟女亚洲av麻豆精品| 看非洲黑人一级黄片| 在线免费观看不下载黄p国产| 99精国产麻豆久久婷婷| av国产精品久久久久影院| h视频一区二区三区| 亚洲综合色惰| 国产成人欧美| 久久免费观看电影| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 亚洲综合色网址| 国产在线一区二区三区精| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| 亚洲经典国产精华液单| 国产精品女同一区二区软件| 精品国产一区二区三区四区第35| 少妇人妻久久综合中文| 啦啦啦啦在线视频资源| 日韩av在线免费看完整版不卡| 在线观看免费高清a一片| 18禁国产床啪视频网站| xxx大片免费视频| 国产乱人偷精品视频| 日本色播在线视频| 下体分泌物呈黄色| 精品熟女少妇av免费看| av国产精品久久久久影院| 亚洲国产精品999| 国产精品三级大全| 丁香六月天网| 久久精品久久久久久噜噜老黄| 在线免费观看不下载黄p国产| 国产成人精品无人区| 观看美女的网站| 国产又色又爽无遮挡免| 在线观看三级黄色| 亚洲少妇的诱惑av| 少妇人妻 视频| 免费看av在线观看网站| 1024视频免费在线观看| 老司机亚洲免费影院| 国产熟女欧美一区二区| 亚洲,一卡二卡三卡| 亚洲久久久国产精品| 中文字幕最新亚洲高清| 亚洲精品久久成人aⅴ小说| 免费大片黄手机在线观看| 欧美日韩一区二区视频在线观看视频在线| av在线播放精品| 99re6热这里在线精品视频| 亚洲av日韩在线播放| av卡一久久| 国产欧美亚洲国产| 日本爱情动作片www.在线观看| 国产成人一区二区在线| 久久久精品94久久精品| 亚洲三级黄色毛片| 日韩大片免费观看网站| 丁香六月天网| 极品人妻少妇av视频| 欧美成人午夜免费资源| 成人黄色视频免费在线看| 亚洲av综合色区一区| 亚洲美女黄色视频免费看| 97精品久久久久久久久久精品| 精品福利永久在线观看| 国产精品久久久久成人av| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩国产mv在线观看视频| 国产又色又爽无遮挡免| 91久久精品国产一区二区三区| 亚洲精品久久成人aⅴ小说| 一边亲一边摸免费视频| 精品一区二区三卡| 午夜激情久久久久久久| 一个人免费看片子| 日韩制服骚丝袜av| 免费黄网站久久成人精品| 日韩伦理黄色片| 中文字幕免费在线视频6| 欧美日韩精品成人综合77777| 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 丰满饥渴人妻一区二区三| 男女边摸边吃奶| 最近2019中文字幕mv第一页| 母亲3免费完整高清在线观看 | 国产熟女午夜一区二区三区| 黄色 视频免费看| av福利片在线| 亚洲精品国产av蜜桃| 精品一区二区三区四区五区乱码 | 午夜福利视频精品| 色94色欧美一区二区| av国产久精品久网站免费入址| 99视频精品全部免费 在线| 王馨瑶露胸无遮挡在线观看| 国产av精品麻豆| 亚洲欧美清纯卡通| 免费女性裸体啪啪无遮挡网站| 欧美人与性动交α欧美软件 | 久久韩国三级中文字幕| 亚洲成色77777| 日韩欧美精品免费久久| 日本免费在线观看一区| 亚洲四区av| 天天操日日干夜夜撸| 亚洲欧美一区二区三区国产| 亚洲一区二区三区欧美精品| 亚洲精品成人av观看孕妇| 蜜臀久久99精品久久宅男| 国产一级毛片在线| 欧美bdsm另类| 亚洲图色成人| 美女国产视频在线观看| 日本91视频免费播放| av又黄又爽大尺度在线免费看| 夜夜爽夜夜爽视频| 日本猛色少妇xxxxx猛交久久| 成人午夜精彩视频在线观看|