• <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
    91成人精品电影| 在线观看午夜福利视频| 黑丝袜美女国产一区| 亚洲国产精品一区二区三区在线| a在线观看视频网站| cao死你这个sao货| 午夜福利欧美成人| 欧美精品一区二区免费开放| 欧美国产精品一级二级三级| 成在线人永久免费视频| 久久久国产欧美日韩av| 手机成人av网站| 国产男女内射视频| 免费女性裸体啪啪无遮挡网站| 国产免费av片在线观看野外av| 欧美日韩精品网址| 亚洲av成人不卡在线观看播放网| 精品国产美女av久久久久小说| 久久精品aⅴ一区二区三区四区| a级片在线免费高清观看视频| 岛国毛片在线播放| 亚洲精品中文字幕一二三四区| 日韩欧美一区视频在线观看| 日韩免费av在线播放| 国产高清国产精品国产三级| 成人国语在线视频| 欧美国产精品一级二级三级| 亚洲第一av免费看| 最新美女视频免费是黄的| 视频区图区小说| 黑丝袜美女国产一区| 两个人免费观看高清视频| 精品亚洲成a人片在线观看| 建设人人有责人人尽责人人享有的| 男人舔女人的私密视频| 国产主播在线观看一区二区| 国产欧美日韩精品亚洲av| 国产有黄有色有爽视频| 亚洲五月天丁香| 一级片'在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 欧美激情久久久久久爽电影 | av电影中文网址| 亚洲精品国产一区二区精华液| av超薄肉色丝袜交足视频| 国产亚洲精品久久久久5区| 黄色毛片三级朝国网站| 精品福利观看| 午夜精品国产一区二区电影| 18禁观看日本| av不卡在线播放| 国内久久婷婷六月综合欲色啪| 中国美女看黄片| 久久久久久久午夜电影 | 欧美久久黑人一区二区| 操美女的视频在线观看| 国产成人av教育| 丝袜人妻中文字幕| 精品福利永久在线观看| 在线观看免费午夜福利视频| 亚洲欧美一区二区三区久久| 91成人精品电影| 少妇裸体淫交视频免费看高清 | 热re99久久精品国产66热6| 乱人伦中国视频| 国产精品久久久人人做人人爽| 啦啦啦 在线观看视频| 亚洲av熟女| 成人手机av| 99久久综合精品五月天人人| 成人精品一区二区免费| 又紧又爽又黄一区二区| 国产成人系列免费观看| 女人久久www免费人成看片| www.999成人在线观看| 免费女性裸体啪啪无遮挡网站| 69精品国产乱码久久久| 国产精品久久电影中文字幕 | 欧美日韩乱码在线| 久久精品亚洲av国产电影网| av网站在线播放免费| 一级毛片精品| 黄色成人免费大全| 啦啦啦免费观看视频1| 免费看a级黄色片| 丝瓜视频免费看黄片| www.自偷自拍.com| 久久狼人影院| 两个人免费观看高清视频| 国产欧美日韩精品亚洲av| 69av精品久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 在线播放国产精品三级| 久久午夜综合久久蜜桃| 老司机午夜福利在线观看视频| 久久国产精品人妻蜜桃| 97人妻天天添夜夜摸| 伊人久久大香线蕉亚洲五| 一区二区三区激情视频| 成人18禁高潮啪啪吃奶动态图| 午夜91福利影院| 在线观看www视频免费| 精品福利观看| 不卡一级毛片| 精品欧美一区二区三区在线| 大香蕉久久网| 无遮挡黄片免费观看| 日本黄色视频三级网站网址 | 国产男女内射视频| 两性夫妻黄色片| 午夜亚洲福利在线播放| 国产午夜精品久久久久久| 欧美精品人与动牲交sv欧美| e午夜精品久久久久久久| 国产精品免费大片| 搡老熟女国产l中国老女人| 免费在线观看影片大全网站| 日韩欧美一区视频在线观看| 免费女性裸体啪啪无遮挡网站| 久久香蕉激情| 丰满人妻熟妇乱又伦精品不卡| 欧美乱码精品一区二区三区| 好男人电影高清在线观看| 91麻豆精品激情在线观看国产 | 91字幕亚洲| 国产欧美日韩一区二区三| 香蕉丝袜av| 免费在线观看黄色视频的| 成人18禁在线播放| 不卡av一区二区三区| 久久人人97超碰香蕉20202| 国产蜜桃级精品一区二区三区 | 国产成人欧美在线观看 | 美女 人体艺术 gogo| 久久国产精品男人的天堂亚洲| 国产精品电影一区二区三区 | 啦啦啦 在线观看视频| 超色免费av| 国产成人系列免费观看| 999精品在线视频| 久久久精品区二区三区| 亚洲精品国产区一区二| 欧美黑人精品巨大| 免费观看人在逋| 久久性视频一级片| 午夜久久久在线观看| 中文字幕精品免费在线观看视频| 热99国产精品久久久久久7| 老司机靠b影院| 中文字幕制服av| 亚洲第一青青草原| 一夜夜www| 国产极品粉嫩免费观看在线| 看片在线看免费视频| 一级黄色大片毛片| 日本欧美视频一区| 久久中文字幕人妻熟女| 午夜精品久久久久久毛片777| 国产精品 欧美亚洲| 欧美日韩瑟瑟在线播放| 美女午夜性视频免费| av有码第一页| 欧美乱妇无乱码| www.精华液| 欧美激情久久久久久爽电影 | 久久草成人影院| 亚洲情色 制服丝袜| 曰老女人黄片| 精品视频人人做人人爽| 一级片免费观看大全| 欧美丝袜亚洲另类 | 午夜福利在线免费观看网站| 香蕉国产在线看| svipshipincom国产片| 国产成+人综合+亚洲专区| 国产精品久久久久久精品古装| 精品国产一区二区三区四区第35| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲欧美98| 日日夜夜操网爽| 久热爱精品视频在线9| 国产男靠女视频免费网站| 久久久水蜜桃国产精品网| 免费人成视频x8x8入口观看| 亚洲第一欧美日韩一区二区三区| 国产高清国产精品国产三级| 国产精品国产高清国产av | 啪啪无遮挡十八禁网站| 亚洲成国产人片在线观看| 男男h啪啪无遮挡| 动漫黄色视频在线观看| 久久天堂一区二区三区四区| 这个男人来自地球电影免费观看| 亚洲第一欧美日韩一区二区三区| www日本在线高清视频| 免费黄频网站在线观看国产| 俄罗斯特黄特色一大片| 久久人人爽av亚洲精品天堂| 国产真人三级小视频在线观看| 男人操女人黄网站| 一级a爱视频在线免费观看| 捣出白浆h1v1| svipshipincom国产片| 亚洲 国产 在线| 高清黄色对白视频在线免费看| 淫妇啪啪啪对白视频| 涩涩av久久男人的天堂| 成人黄色视频免费在线看| 好看av亚洲va欧美ⅴa在| 在线观看免费视频网站a站| 精品免费久久久久久久清纯 | 黄色 视频免费看| 欧美大码av| 精品国内亚洲2022精品成人 | 啦啦啦 在线观看视频| 另类亚洲欧美激情| 老熟女久久久| 久久国产乱子伦精品免费另类| 日本欧美视频一区| 精品国产一区二区三区久久久樱花| 日本撒尿小便嘘嘘汇集6| 国产99久久九九免费精品| 亚洲欧美日韩另类电影网站| 乱人伦中国视频| 国产免费现黄频在线看| 久久久国产成人免费| 一级毛片女人18水好多| 建设人人有责人人尽责人人享有的| 中文亚洲av片在线观看爽 | 日韩免费av在线播放| 国产亚洲精品第一综合不卡| 日韩欧美免费精品| 国产成人欧美在线观看 | 国产激情欧美一区二区| 亚洲av片天天在线观看| 婷婷精品国产亚洲av在线 | 成在线人永久免费视频| 三级毛片av免费| 欧美中文综合在线视频| 欧美日韩亚洲国产一区二区在线观看 | 日本五十路高清| 在线播放国产精品三级| 国产激情久久老熟女| 国产aⅴ精品一区二区三区波| 我的亚洲天堂| 国产免费男女视频| 久久久国产欧美日韩av| 国产精品 国内视频| 搡老乐熟女国产| 身体一侧抽搐| 国产三级黄色录像| 女人被狂操c到高潮| 国产免费男女视频| 男人舔女人的私密视频| 国产一区二区三区在线臀色熟女 | 黄色视频,在线免费观看| 日本一区二区免费在线视频| 国产激情欧美一区二区| 久久人人97超碰香蕉20202| 国产精品一区二区在线观看99| 后天国语完整版免费观看| 窝窝影院91人妻| 国产视频一区二区在线看| 亚洲男人天堂网一区| 国产欧美日韩一区二区精品| 真人做人爱边吃奶动态| 久久国产乱子伦精品免费另类| 69av精品久久久久久| 久久ye,这里只有精品| 国产亚洲一区二区精品| 精品久久久久久久毛片微露脸| 国产单亲对白刺激| 一边摸一边做爽爽视频免费| 日韩三级视频一区二区三区| 久久久久久久久久久久大奶| 国产激情欧美一区二区| 成年版毛片免费区| 欧美日本中文国产一区发布| 如日韩欧美国产精品一区二区三区| 亚洲免费av在线视频| xxx96com| 黄色 视频免费看| 日韩精品免费视频一区二区三区| 最近最新免费中文字幕在线| 两个人看的免费小视频| 日韩制服丝袜自拍偷拍| 51午夜福利影视在线观看| 成年女人毛片免费观看观看9 | 中文字幕人妻丝袜一区二区| 欧美黄色淫秽网站| 精品少妇一区二区三区视频日本电影| 99re6热这里在线精品视频| 如日韩欧美国产精品一区二区三区| 啦啦啦视频在线资源免费观看| 新久久久久国产一级毛片| 亚洲av欧美aⅴ国产| 亚洲国产中文字幕在线视频| 午夜91福利影院| 国产精品香港三级国产av潘金莲| 国产成人影院久久av| av国产精品久久久久影院| 十八禁高潮呻吟视频| 大型av网站在线播放| 欧美日韩亚洲国产一区二区在线观看 | 久久人妻av系列| 久久久久久免费高清国产稀缺| 精品久久久久久久毛片微露脸| 91麻豆av在线| 久久久精品区二区三区| 91成年电影在线观看| 首页视频小说图片口味搜索| 水蜜桃什么品种好| 黑人巨大精品欧美一区二区蜜桃| 国产精品 国内视频| 免费观看精品视频网站| av天堂在线播放| 亚洲精品中文字幕一二三四区| 在线观看免费午夜福利视频| svipshipincom国产片| 国产一区二区三区视频了| 人人妻人人澡人人爽人人夜夜| 国产在线精品亚洲第一网站| 久久草成人影院| 久久ye,这里只有精品| 搡老岳熟女国产| 99久久国产精品久久久| 久久精品国产a三级三级三级| 久99久视频精品免费| 欧美日韩黄片免| 免费看十八禁软件| 亚洲欧美精品综合一区二区三区| 18禁美女被吸乳视频| 亚洲一区二区三区不卡视频| 久久国产精品大桥未久av| 国产1区2区3区精品| 日本黄色日本黄色录像| 国精品久久久久久国模美| 国产精品免费一区二区三区在线 | 国产精品影院久久| 两人在一起打扑克的视频| 12—13女人毛片做爰片一| 欧美日韩成人在线一区二区| 又黄又爽又免费观看的视频| 亚洲精品中文字幕一二三四区| 欧美乱色亚洲激情| 亚洲人成77777在线视频| 99国产精品99久久久久| 老司机影院毛片| 国产1区2区3区精品| 在线天堂中文资源库| 国产高清国产精品国产三级| 国产高清视频在线播放一区| 欧美日本中文国产一区发布| 国产精品亚洲av一区麻豆| 一本大道久久a久久精品| 黄色女人牲交| 波多野结衣一区麻豆| 王馨瑶露胸无遮挡在线观看| 建设人人有责人人尽责人人享有的| 捣出白浆h1v1| 人成视频在线观看免费观看| 久久久久久久久免费视频了| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 国产精品电影一区二区三区 | 欧美丝袜亚洲另类 | 色播在线永久视频| 丰满的人妻完整版| 免费在线观看亚洲国产| 中文亚洲av片在线观看爽 | 女人爽到高潮嗷嗷叫在线视频| 亚洲男人天堂网一区| 国产精品一区二区免费欧美| 人妻一区二区av| 国产aⅴ精品一区二区三区波| 天天躁狠狠躁夜夜躁狠狠躁| 女人久久www免费人成看片| 欧美黄色片欧美黄色片| 精品无人区乱码1区二区| 少妇 在线观看| 久久国产精品大桥未久av| 久久久久精品国产欧美久久久| 亚洲国产毛片av蜜桃av| 99香蕉大伊视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲人成77777在线视频| 国产片内射在线| 亚洲色图av天堂| 国产精品av久久久久免费| 亚洲中文日韩欧美视频| 热99久久久久精品小说推荐| 国产色视频综合| 午夜福利影视在线免费观看| 日韩欧美三级三区| 久久久久视频综合| ponron亚洲| 亚洲精华国产精华精| 成年人午夜在线观看视频| 男女下面插进去视频免费观看| 在线十欧美十亚洲十日本专区| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 午夜福利在线免费观看网站| 国产99久久九九免费精品| 一a级毛片在线观看| 免费女性裸体啪啪无遮挡网站| 成人免费观看视频高清| 午夜精品久久久久久毛片777| 美国免费a级毛片| 黄色女人牲交| 免费在线观看亚洲国产| 亚洲熟女毛片儿| 国产精品自产拍在线观看55亚洲 | 男人舔女人的私密视频| 精品无人区乱码1区二区| 午夜影院日韩av| 性少妇av在线| 国产男靠女视频免费网站| 国产国语露脸激情在线看| 99国产极品粉嫩在线观看| 国产97色在线日韩免费| 成年动漫av网址| 成人18禁在线播放| 久久久久精品人妻al黑| 亚洲成人手机| 精品免费久久久久久久清纯 | 亚洲熟女精品中文字幕| 婷婷精品国产亚洲av在线 | 国产精品免费一区二区三区在线 | 宅男免费午夜| 看免费av毛片| 国产有黄有色有爽视频| 国产精品98久久久久久宅男小说| 国产精品偷伦视频观看了| 亚洲精品中文字幕在线视频| 在线观看日韩欧美| 国产日韩欧美亚洲二区| 欧美日韩亚洲国产一区二区在线观看 | 久久人妻熟女aⅴ| 又黄又爽又免费观看的视频| 一边摸一边抽搐一进一小说 | 日本撒尿小便嘘嘘汇集6| 亚洲第一av免费看| 国产日韩一区二区三区精品不卡| 精品人妻熟女毛片av久久网站| 黄色毛片三级朝国网站| 18禁美女被吸乳视频| 国产精品久久电影中文字幕 | 国产单亲对白刺激| 丰满的人妻完整版| 啦啦啦在线免费观看视频4| 亚洲中文av在线| 久久精品亚洲精品国产色婷小说| 国产精品久久久久久人妻精品电影| 久久国产精品影院| 日韩有码中文字幕| 久久国产乱子伦精品免费另类| 久久久国产一区二区| 久久中文看片网| 国产1区2区3区精品| 午夜亚洲福利在线播放| 欧美成人午夜精品| 久久久精品国产亚洲av高清涩受| 搡老岳熟女国产| 久久ye,这里只有精品| 黄色视频不卡| 欧美最黄视频在线播放免费 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久久成人av| 欧美日韩中文字幕国产精品一区二区三区 | 999久久久精品免费观看国产| 国产成人精品久久二区二区91| 亚洲一区二区三区不卡视频| 99国产综合亚洲精品| 国产日韩一区二区三区精品不卡| 色在线成人网| 免费在线观看黄色视频的| 妹子高潮喷水视频| 国产精品久久视频播放| 9热在线视频观看99| 丝袜美腿诱惑在线| 国产成人精品久久二区二区免费| 欧美激情极品国产一区二区三区| 亚洲欧美日韩高清在线视频| 18禁国产床啪视频网站| 正在播放国产对白刺激| 色在线成人网| 精品久久蜜臀av无| 欧美黑人精品巨大| 啦啦啦视频在线资源免费观看| 久久人妻熟女aⅴ| 亚洲熟女毛片儿| 亚洲人成伊人成综合网2020| 80岁老熟妇乱子伦牲交| 男男h啪啪无遮挡| cao死你这个sao货| 国产一区二区三区视频了| 久久久久久亚洲精品国产蜜桃av| avwww免费| 国产精品久久电影中文字幕 | 人人妻人人澡人人看| 国产成人系列免费观看| 搡老乐熟女国产| 一个人免费在线观看的高清视频| 99国产精品一区二区三区| 亚洲av第一区精品v没综合| 不卡av一区二区三区| 高清毛片免费观看视频网站 | 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av高清一级| 久久久精品区二区三区| 欧美黄色淫秽网站| 午夜福利,免费看| 精品久久久精品久久久| 久久精品91无色码中文字幕| 国产欧美日韩一区二区三区在线| 高清av免费在线| 久久天堂一区二区三区四区| 精品国产国语对白av| 啪啪无遮挡十八禁网站| 亚洲专区中文字幕在线| 午夜成年电影在线免费观看| 女人久久www免费人成看片| 亚洲午夜精品一区,二区,三区| 国产熟女午夜一区二区三区| 亚洲国产精品sss在线观看 | 国产成+人综合+亚洲专区| 高清毛片免费观看视频网站 | 国产欧美日韩一区二区三| 美女福利国产在线| 国产亚洲精品久久久久5区| 亚洲成人免费电影在线观看| 最新在线观看一区二区三区| 人妻一区二区av| 无限看片的www在线观看| 18在线观看网站| 国产男靠女视频免费网站| 久久ye,这里只有精品| 亚洲性夜色夜夜综合| 亚洲专区国产一区二区| 日本vs欧美在线观看视频| 国产免费av片在线观看野外av| 丝袜美足系列| 亚洲人成电影观看| 在线视频色国产色| 人人澡人人妻人| 亚洲九九香蕉| 成人特级黄色片久久久久久久| 久久久久久亚洲精品国产蜜桃av| 美女高潮喷水抽搐中文字幕| 不卡av一区二区三区| 免费观看a级毛片全部| 国产不卡一卡二| 看片在线看免费视频| 精品福利永久在线观看| 丝瓜视频免费看黄片| 亚洲欧洲精品一区二区精品久久久| 久久国产精品大桥未久av| 黄网站色视频无遮挡免费观看| 99热网站在线观看| 高清视频免费观看一区二区| 欧美精品啪啪一区二区三区| 国产欧美日韩综合在线一区二区| 80岁老熟妇乱子伦牲交| 欧美日韩国产mv在线观看视频| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 亚洲 欧美一区二区三区| 夜夜爽天天搞| 热re99久久精品国产66热6| 日韩有码中文字幕| 新久久久久国产一级毛片| 美女高潮到喷水免费观看| 成熟少妇高潮喷水视频| 久久狼人影院| 久久热在线av| 国产精品二区激情视频| 91精品国产国语对白视频| 国产精品一区二区在线观看99| 午夜福利影视在线免费观看| 亚洲七黄色美女视频| 9热在线视频观看99| 桃红色精品国产亚洲av| 999久久久国产精品视频| 国产成人免费观看mmmm| 国产真人三级小视频在线观看| 国产精品永久免费网站| 精品人妻熟女毛片av久久网站| 亚洲国产欧美网| 高清黄色对白视频在线免费看| 女同久久另类99精品国产91| 久久午夜亚洲精品久久| 国产精品综合久久久久久久免费 | 曰老女人黄片| 国产精华一区二区三区| 国产欧美日韩一区二区三区在线| av片东京热男人的天堂| 久久国产精品男人的天堂亚洲| 制服人妻中文乱码| 久久中文字幕一级| 一进一出抽搐gif免费好疼 | 国产精品亚洲一级av第二区| 久久草成人影院| 久久久国产精品麻豆| 丰满饥渴人妻一区二区三| 国产97色在线日韩免费| 黄色视频不卡| 麻豆乱淫一区二区| 亚洲成a人片在线一区二区| 日韩欧美一区视频在线观看| 国产激情欧美一区二区| 亚洲av日韩在线播放| 丰满饥渴人妻一区二区三| 久久久久久久午夜电影 | 欧美日韩一级在线毛片|