陸面過(guò)程參數(shù)化對(duì)寧波地區(qū)雷暴過(guò)程模擬的影響

汪雅,苗峻峰,談?wù)苊?/p>

(1.南京信息工程大學(xué) 大氣科學(xué)學(xué)院,江蘇 南京 210044;2.上海市氣象局,上海 200030;

3.南京大學(xué) 中尺度災(zāi)害性天氣教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210093)

?

陸面過(guò)程參數(shù)化對(duì)寧波地區(qū)雷暴過(guò)程模擬的影響

汪雅1,2,苗峻峰1,談?wù)苊?

(1.南京信息工程大學(xué) 大氣科學(xué)學(xué)院,江蘇 南京 210044;2.上海市氣象局,上海 200030;

3.南京大學(xué) 中尺度災(zāi)害性天氣教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210093)

摘要：利用耦合復(fù)雜程度不同的陸面過(guò)程參數(shù)化方案(Noah、RUC)的新一代中尺度數(shù)值模式WRF,對(duì)2006年6月24日發(fā)生在寧波地區(qū)的一次典型的雷暴過(guò)程進(jìn)行了數(shù)值模擬試驗(yàn)。結(jié)果表明:對(duì)于雷暴發(fā)生前期近地面熱力、動(dòng)力場(chǎng)的特征,Noah方案的模擬較為逼真,RUC方案沒(méi)有反映出下墊面覆蓋的多樣性以及城市下墊面的影響,城鄉(xiāng)之間差異不明顯;Noah方案模擬的雷暴啟動(dòng)、發(fā)展過(guò)程與觀測(cè)較為一致,RUC方案較好地描述出了演變過(guò)程中的關(guān)鍵階段(3次合并過(guò)程);由于參數(shù)化所考慮的要素和物理過(guò)程存在一定差異,Noah方案在對(duì)降水的強(qiáng)度、降水中心位置的模擬方面具有一定優(yōu)勢(shì);雷暴的持續(xù)時(shí)間對(duì)陸面過(guò)程參數(shù)化方案的選擇比較敏感,兩個(gè)方案所模擬的雷暴過(guò)程持續(xù)時(shí)間不同程度地長(zhǎng)于實(shí)際雷暴持續(xù)時(shí)間;無(wú)論是哪種下墊面覆蓋類(lèi)型,白天Noah方案模擬的感熱通量均大于RUC方案,而Noah方案模擬的潛熱通量均小于RUC方案。

關(guān)鍵詞：雷暴過(guò)程;陸面過(guò)程參數(shù)化;Noah方案;RUC方案;數(shù)值敏感性

0引言

雷暴是在發(fā)展強(qiáng)烈的積雨云中產(chǎn)生的,伴有雷電、大風(fēng)和冰雹等多種天氣現(xiàn)象的中小尺度天氣過(guò)程,是我國(guó)夏季常見(jiàn)且嚴(yán)重的自然災(zāi)害之一。它發(fā)生發(fā)展的突發(fā)性和局地性強(qiáng),預(yù)報(bào)難度大,因而對(duì)雷暴的觀測(cè)和研究工作一直以來(lái)受到氣象學(xué)者們的廣泛關(guān)注。我國(guó)對(duì)雷暴的研究起步相對(duì)較晚,較多的工作為依賴(lài)于大尺度環(huán)境背景場(chǎng)和探空資料針對(duì)個(gè)例的天氣分析(李志楠和康玉霞,1995)或多個(gè)個(gè)例的氣候?qū)W統(tǒng)計(jì)(張敏鋒和馮霞,1998;陳思蓉等,2009),但由于時(shí)間和空間的分辨率較低,這些研究工作對(duì)于雷暴的認(rèn)識(shí)非常有限。隨著加密觀測(cè)網(wǎng)絡(luò)的布置和天氣雷達(dá)的不斷發(fā)展,雷暴在發(fā)生發(fā)展過(guò)程中精細(xì)特征的識(shí)別及其位置的精確定位逐漸成為可能(肖輝等,2002),但目前天氣雷達(dá)難以分辨到近地層的回波特征,并且很難深入地揭示出雷暴過(guò)程的影響機(jī)制(李昀英等,2008)。

近年來(lái),隨著計(jì)算機(jī)的快速發(fā)展和數(shù)值研究的逐漸深入,雷暴等強(qiáng)對(duì)流天氣過(guò)程的數(shù)值模擬研究工作也陸續(xù)開(kāi)展。何宏讓等(1999)采用一個(gè)二維云模式,對(duì)副熱帶高壓后部的強(qiáng)雷暴的產(chǎn)生及發(fā)展進(jìn)行了數(shù)值模擬,并探討了不同初始擾動(dòng)方案對(duì)模擬效果的影響。陳力強(qiáng)等(2004)運(yùn)用早期版本的MM5模式對(duì)東北冷渦誘發(fā)的強(qiáng)風(fēng)暴過(guò)程的垂直環(huán)境條件和不穩(wěn)定能量積累機(jī)制進(jìn)行了分析和討論。張春喜等(2005)在對(duì)山東地區(qū)的中尺度對(duì)流系統(tǒng)的大尺度環(huán)境進(jìn)行診斷分析的基礎(chǔ)上,用WRFV2通過(guò)對(duì)比試驗(yàn)的方法研究了初始對(duì)流的觸發(fā)機(jī)制。

上述數(shù)值模擬工作側(cè)重于中尺度動(dòng)力過(guò)程本身的研究,所使用的模式動(dòng)力學(xué)框架較為簡(jiǎn)單、精度不高,另外對(duì)于陸面(大氣下邊界)過(guò)程的影響被較少地考慮到。事實(shí)上,陸面過(guò)程勢(shì)必影響到中短期天氣的演變,因?yàn)樘?yáng)輻射、大氣逆輻射的能量直接由陸面過(guò)程分配為感熱和潛熱,并反饋到大氣(曾新民等,2012)。同時(shí),數(shù)值模式中對(duì)云以及大氣邊界層溫、濕度特征的準(zhǔn)確描述也主要依賴(lài)于下墊面的土壤、植被特征和物理生化過(guò)程(Xiu and Pleim,2001;Miao et al.,2008,2009)。因此,在對(duì)中短期天氣過(guò)程的數(shù)值模擬中,陸面方案的選取起著非常重要的作用。近年來(lái),不少研究表明中小尺度天氣過(guò)程模擬的結(jié)果對(duì)陸面參數(shù)化方案的選取具有依賴(lài)性,如,馬紅云等(2009)用WRF模式討論了不同陸面參數(shù)化方案對(duì)暴雨過(guò)程的影響,其研究結(jié)果表明:強(qiáng)降水中心位置、降水強(qiáng)度的分布對(duì)陸面物理過(guò)程是敏感的。曾新民等(2012)針對(duì)南方地區(qū)強(qiáng)降水事件也作了類(lèi)似的研究,并指出陸面方案通過(guò)影響地表蒸發(fā)量以及低層環(huán)流及水汽輻合,從而影響大范圍降水。然而,以上的研究區(qū)域主要集中于江淮之間和華南一帶,對(duì)浙江等地的研究較少;研究?jī)?nèi)容也多限于較強(qiáng)的暴雨事件,針對(duì)雷暴等強(qiáng)對(duì)流天氣的研究較少。

寧波位于我國(guó)東海沿岸,海陸風(fēng)現(xiàn)象較為顯著(宋潔慧等,2009)。海陸風(fēng)的發(fā)生發(fā)展會(huì)影響沿岸地區(qū)溫度場(chǎng)、濕度場(chǎng)、風(fēng)場(chǎng)的分布,并引起低層大氣層結(jié)狀況的改變,與沿岸積云對(duì)流系統(tǒng)的分布、移動(dòng)和發(fā)展有著密切的關(guān)系(Talbot et al.,2007;王玉國(guó)等,2009;汪雅等,2013;苗峻峰,2014)。統(tǒng)計(jì)顯示,寧波境內(nèi)全年雷暴日均在40 d以上,屬于雷暴高發(fā)地區(qū)。到目前為止,針對(duì)寧波地區(qū)雷暴等強(qiáng)對(duì)流天氣的模擬研究仍然很少。因此,本文以寧波及周邊地區(qū)為研究區(qū)域,在汪雅等(2013)的基礎(chǔ)上,利用新一代中尺度模式WRF分別耦合不同的陸面過(guò)程參數(shù)化方案對(duì)發(fā)生在該區(qū)域的一次強(qiáng)雷暴過(guò)程進(jìn)行了數(shù)值模擬,目的在于探討雷暴過(guò)程數(shù)值模擬對(duì)陸面過(guò)程參數(shù)化方案的敏感性,為天氣研究和預(yù)報(bào)業(yè)務(wù)中合理選擇和使用中尺度模式陸面過(guò)程參數(shù)化方案提供一些依據(jù)。

1雷暴個(gè)例與天氣過(guò)程簡(jiǎn)介

2006年6月24日14—23時(shí)(北京時(shí)間,下同),寧波市出現(xiàn)了強(qiáng)雷暴天氣,局部遭受8～10級(jí)的大風(fēng)和短時(shí)暴雨的襲擊。根據(jù)寧波市氣象局降水資料記錄,有11個(gè)站1 h雨量達(dá)到30 mm以上,屬短時(shí)極強(qiáng)暴雨。

1.1　天氣尺度背景

有利的大尺度環(huán)流形勢(shì)為本次雷暴的形成提供了良好的天氣背景條件:從該日08時(shí)由NCEP FNL 1°×1°資料分析的500和850 hPa等壓面環(huán)流形勢(shì)(圖1)可以看出,此時(shí)500 hPa高度上內(nèi)蒙古北側(cè)受低壓控制,副熱帶高壓脊線位于24°N附近,寧波位于高壓的頂后部,水汽從南海沿高壓西北側(cè)的西南氣流向浙江東部輸送。低層(850 hPa)上,寧波地區(qū)受西南氣流控制,無(wú)明顯切變線系統(tǒng)。地面上寧波位于低壓中心附近,風(fēng)場(chǎng)有氣旋式輻合(圖略)。24日上午,寧波云量較少,受短波輻射地表強(qiáng)烈加熱,下層暖濕上層干冷的垂直結(jié)構(gòu)形成強(qiáng)烈的不穩(wěn)定層結(jié),為對(duì)流的發(fā)展提供了充足的能量。圖1的陰影區(qū)分別給出了14時(shí)的地面抬升指數(shù)和對(duì)流有效位能,此時(shí)寧波地區(qū)具有較高的對(duì)流有效位能,地面抬升指數(shù)也為相對(duì)高的負(fù)值區(qū),有利于強(qiáng)對(duì)流的爆發(fā)。

圖1　由NCEP FNL 1°×1°資料分析的2006年6月24日08時(shí)風(fēng)場(chǎng)(箭矢;單位:m·s-1)和位勢(shì)高度場(chǎng)(等值線;單位:dagpm)　　a.500 hPa(陰影表示14時(shí)地面抬升指數(shù);單位:℃);b.850 hPa(陰影表示14時(shí)對(duì)流有效位能;單位:J·kg-1)Fig.1　Wind field(arrows;units:m·s-1) and geopotential height field(contours;units:dagpm) at 08:00 BST 24 June 2006 based on NCEP FNL 1°×1° reanalysis data　　a.500 hPa(shaded area represents the lift index at 14:00 BST;units:℃);b.850 hPa(shaded area represents convective available potential energy at 14:00 BST;units:J·kg-1)

1.2　雷達(dá)資料反映的對(duì)流發(fā)展特征

圖2給出了2006年6月24日午后由寧波雷達(dá)觀測(cè)到的對(duì)流在寧波及周邊地區(qū)的發(fā)展情況。13時(shí)之后,在寧波市沿海、慈溪余姚南部一線有輻合帶生成。14時(shí)以后,沿海生成的輻合帶上開(kāi)始有強(qiáng)對(duì)流回波生成(圖2a)。同時(shí),從杭州西北部生成的東移的強(qiáng)回波到達(dá)杭州灣頂,并逐漸影響余姚和慈溪(圖2a MCS1)。隨后,外來(lái)回波和本地生成回波一并向南發(fā)展傳播。16時(shí)在寧海北部和象山港(圖3b中的“Xiangshan Bay”)北側(cè)有強(qiáng)回波生成,向東偏北方向移動(dòng),與前面南移的局地回波在象山港北岸合并加強(qiáng)為MCS2(圖2b)。此后,合并的回波向南發(fā)展,到象山港上空與MCS1合并為MCS12(圖2c)。同時(shí),寧波南部有強(qiáng)回波(MCS3)向北移動(dòng),18時(shí)之后MCS3與MCS12在象山港山上空合并加強(qiáng)為MCS123并向南傳播,入海后強(qiáng)度減弱并逐漸消亡。另外,18時(shí)后寧波西南側(cè)有回波單體生成,并逐漸東移。由上可見(jiàn),本次雷暴過(guò)程是一次典型的外來(lái)回波與局地生成回波綜合作用而發(fā)展加強(qiáng)的強(qiáng)天氣過(guò)程。

圖2　2006年6月24日寧波雷達(dá)觀測(cè)到的雷達(dá)基本反射率(0.5°仰角;單位:dBz)a.14時(shí)51分;b.16時(shí)20分;c.17時(shí)40分;d.19時(shí)16分Fig.2　Radar basic reflectivity images observed by Ningbo Radar on 24 June 2006(elevation angle is 0.5°;units:dBz)　　a.14:51 BST;b.16:20 BST;c.17:40 BST;d.19:16 BST

2陸面過(guò)程參數(shù)化方案簡(jiǎn)介

本文數(shù)值模擬研究采用的模式是中尺度模式WRF,它重點(diǎn)考慮從云尺度到天氣尺度等重要天氣的預(yù)報(bào),在目前的短期天氣預(yù)報(bào)業(yè)務(wù)和研究中得到了廣泛應(yīng)用。與它的前身MM5中尺度模式相比,WRF模式考慮了較為詳細(xì)的陸面過(guò)程,可以描述不同下墊面的感熱、潛熱通量的傳輸過(guò)程,并將其作為下邊界條件耦合到大氣模式中,進(jìn)一步加強(qiáng)云、輻射、降水等過(guò)程的模擬。

本文選取了兩個(gè)比較復(fù)雜的并且被廣泛使用的陸面過(guò)程參數(shù)化方案來(lái)研究雷暴過(guò)程數(shù)值模擬對(duì)不同陸面過(guò)程參數(shù)化方案的依賴(lài)性:1)Noah方案,源于早期的俄勒岡州立大學(xué)(OSU)陸面模式(Pan and Mahrt,1987),其后多個(gè)機(jī)構(gòu)對(duì)其進(jìn)行了合作開(kāi)發(fā)并不斷完善逐漸發(fā)展為現(xiàn)在版本。從WRF V2.0之后,Noah方案考慮了城市下墊面的影響,許多與城市相關(guān)的參數(shù)(表面反照率、粗糙度、體積熱容量、土壤熱傳導(dǎo)率、植被成分)在Noah模式中被作了調(diào)整和改進(jìn)(Miao et al.,2007)。2)RUC方案,采用了美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)天氣業(yè)務(wù)預(yù)報(bào)系統(tǒng)RUC(Rapid Update Cycle)模式中的陸面物理參數(shù)化方案(Smirnova et al.,1997)。兩種陸面過(guò)程參數(shù)化方案的主要差異見(jiàn)表1和楊薇等(2014)。

3模式設(shè)置和試驗(yàn)方案

本研究采用的模式版本為WRF-ARW V3.2,模式的初始場(chǎng)和邊界條件由NCEP FNL 1°×1°再分析資料提供。模式使用四重嵌套方案,水平格距分別為27、9、3、1 km,垂直方向均具有不等間距的35個(gè)σ層,其中2 km以下設(shè)置24層,模式層頂氣壓設(shè)為100 hPa。圖3為模擬所用網(wǎng)格嵌套示意圖與最內(nèi)層網(wǎng)格區(qū)域的地形及主要自動(dòng)氣象站的分布。最外層區(qū)域?yàn)槲覈?guó)中東大部地區(qū),提供大的背景強(qiáng)迫。最內(nèi)層網(wǎng)格區(qū)域D4覆蓋了寧波及其周邊地區(qū),為

圖3　模擬區(qū)域示意圖　　a.四重嵌套區(qū)域;b.D4網(wǎng)格區(qū)域的地形(單位:m)及部分自動(dòng)氣象站的分布;c.D4區(qū)域陸面覆蓋類(lèi)型Fig.3　Modeling domain settings of the simulation　　a.Domains 1,2,3,and 4(denoted by D1,D2,D3,and D4);b.Model terrain(units:m) in D4 and the distribution of automatic weather stations;c.land cover category in D4

本文的研究區(qū)域(圖3b),格點(diǎn)數(shù)160×166。短波輻射采用Dudhia方案,長(zhǎng)波輻射采用RRTM方案,近地面層采用Monin-Obukhov方案,邊界層選用YSU方案,微物理過(guò)程采用Lin方案,積云對(duì)流采用Kain-Fritsch方案(D3、D4不采用)。陸面覆蓋采用Boston University/NCEP提供的MODIS 30″格點(diǎn)資料,其土地利用/土地覆蓋類(lèi)型分布如圖3c所示,其中主要類(lèi)型為:水體(藍(lán)色),混合林(綠色),農(nóng)田(黃色),城市及建筑用地(紅色)。模式結(jié)構(gòu)詳見(jiàn)汪雅等(2013)。

為了比較不同的陸面過(guò)程參數(shù)化方案對(duì)2006年6月24日的這次雷暴過(guò)程模擬的影響,本文選擇兩個(gè)比較復(fù)雜的并且被廣泛使用的Noah和RUC陸面參數(shù)化方案分別對(duì)本次雷暴過(guò)程作模擬試驗(yàn)。模式積分時(shí)間為2006年6月23日08時(shí)至25日08時(shí),模擬結(jié)果1 h輸出一次,各試驗(yàn)前16 h為模式積分調(diào)整(spin-up)時(shí)間。除了陸面過(guò)程參數(shù)化方案外,兩個(gè)試驗(yàn)所選的物理過(guò)程和模式配置完全一致,因此不同試驗(yàn)?zāi)M結(jié)果的差異可以看作是由不同的陸面過(guò)程參數(shù)化方案所引起的。如不做特殊說(shuō)明,文中對(duì)模擬結(jié)果的分析均是針對(duì)D4區(qū)域。

表1Noah與RUC陸面過(guò)程參數(shù)化方案的差異

Table 1Differences between Noah and RUC land surface process parameterization schemes

Noah方案RUC方案預(yù)報(bào)變量4個(gè)土壤層上的土壤溫度和濕度,地表溫度,冠層水分,積雪厚度和密度。6個(gè)土壤層的土壤溫度和濕度,地表溫度,冠層水分,積雪厚度和密度。植被參數(shù)化①采用隱式方案求解能量和水分平衡方程;②植被影響包括在其中,植被指數(shù)選項(xiàng)隨時(shí)間變化;③處理蒸散時(shí)使用土壤和植被類(lèi)型,植被參數(shù)(如:地面粗糙度和反照率)依賴(lài)于植被類(lèi)型,土壤熱力特性依賴(lài)于土壤類(lèi)型;④考慮了植被的蒸發(fā)作用,植物的蒸騰作用以及植被的截留作用⑤在地表交換系數(shù)中引入熱力粗糙度,潛在蒸發(fā)通過(guò)求解彭曼方程得到;⑥為行星邊界層提供熱量和水汽通量。①能量和水分平衡采用隱式方案求解;②植被影響包括在其中,地面粗糙度作為植被類(lèi)型的函數(shù)引入;③把植被與裸土的蒸發(fā)分開(kāi)考慮,植被的蒸騰作用較為簡(jiǎn)單,對(duì)參數(shù)化并不敏感;④各層能量和水汽收支近似平衡;⑤包括3個(gè)蒸散分量:裸土蒸發(fā)作用,植被蒸發(fā)作用,植被蒸騰作用;更多的長(zhǎng)波輻射冷卻;⑥為行星邊界層提供熱量和水汽通量。土壤參數(shù)化①4個(gè)土壤層,自上而下厚度依次為:10,30,60,100cm;②考慮熱力擴(kuò)散和水汽輸送,采用迭代法求解對(duì)角矩陣來(lái)更新土壤溫度。①6個(gè)土壤層,自上而下厚度依次為:2.5,2.5,15,20,120,140cm,頂層土壤較薄;②土壤濕度由Richards水分傳輸方程求解,土壤導(dǎo)水率由二階數(shù)值逼近法得出。

4模擬結(jié)果分析

4.1　陸面參數(shù)化方案對(duì)對(duì)流系統(tǒng)發(fā)生、發(fā)展過(guò)程的影響

由寧波雷達(dá)實(shí)際觀測(cè)結(jié)果(1.2節(jié))可以看出,本次雷暴過(guò)程中主要存在3個(gè)中尺度對(duì)流系統(tǒng)的發(fā)展演變過(guò)程。與圖2相對(duì)應(yīng),圖4給出了Noah方案和RUC方案模擬的對(duì)應(yīng)時(shí)間約1.5 km高度雷達(dá)回波。兩個(gè)方案清晰地再現(xiàn)了3個(gè)MCS初始、發(fā)展、合并與增強(qiáng)的過(guò)程。24日13時(shí),兩個(gè)方案均模擬出了局地生成回波,但RUC方案模擬的回波強(qiáng)度較小且位置分散,系統(tǒng)發(fā)展過(guò)程相對(duì)滯后。圖4b、f反映出兩個(gè)方案都模擬出了北來(lái)的回波系統(tǒng)已移動(dòng)到慈溪—余姚一帶,但是RUC方案模擬的MCS1回波強(qiáng)度和回波面積偏小。由圖4c、g可見(jiàn),RUC方案模擬出了MCS1與MCS2的合并過(guò)程,而Noah方案并沒(méi)有較好地展示出回波系統(tǒng)的合并過(guò)程。18時(shí)之后,兩個(gè)方案均展現(xiàn)出了3個(gè)系統(tǒng)的合并過(guò)程(圖4d、g),然而RUC方案所模擬的系統(tǒng)影響范圍較小,只有Noah的2/3,位置偏西南約20 km。此外,兩個(gè)方案均模擬出了18時(shí)之后于寧波西部生成的回波單體,該系統(tǒng)形成時(shí)間與實(shí)際相比晚了約1 h。縱觀整個(gè)雷暴演變過(guò)程,兩個(gè)方案分別在不同階段展現(xiàn)了各自的優(yōu)勢(shì):Noah方案對(duì)雷暴的啟動(dòng)與發(fā)展過(guò)程及其雷暴的強(qiáng)度和位置的模擬較為合理,而RUC方案較好地描述出了雷暴演變過(guò)程中的關(guān)鍵階段,如三次合并過(guò)程。

圖4　模擬的2006年6月24日模式D4區(qū)域約1.5 km高度的雷達(dá)回波(陰影;單位:dBz)和10 m風(fēng)場(chǎng)(單位:m·s-1)a.Noah方案13時(shí);b.Noah方案15時(shí);c.Noah方案18時(shí);d.Noah方案20時(shí);e.RUC方案13時(shí);f.RUC方案15時(shí);g.RUC方案18時(shí);h.RUC方案20時(shí)Fig.4　Simulated radar reflectivity(shaded areas;units:dBz) at 1.5 km level and wind field at 10 m(units:m·s-1) in D4 on 24 June 2006　　a.Noah scheme at 13:00 BST;b.Noah scheme at 15:00 BST;c.Noah scheme at 18:00 BST;d.Noah scheme at 20:00 BST;e.RUC scheme at 13:00 BST;f.RUC scheme at 15:00 BST;g.RUC scheme at 18:00 BST;h.RUC scheme at 20:00 BST

4.2　陸面參數(shù)化方案對(duì)近地面熱動(dòng)力結(jié)構(gòu)的影響

陸氣之間不斷進(jìn)行的熱量、動(dòng)量和水汽交換對(duì)大氣邊界層的熱動(dòng)力結(jié)構(gòu)產(chǎn)生重要影響,進(jìn)而可影響到雷暴等強(qiáng)對(duì)流天氣的發(fā)生和發(fā)展。蒙偉光等(2005)的研究表明陸地表面可通過(guò)對(duì)前期溫度場(chǎng)、氣壓場(chǎng)、流場(chǎng)的影響,對(duì)對(duì)流的啟動(dòng)和中尺度對(duì)流的形成起作用。

圖5a、b顯示了雷暴發(fā)生前期不同方案模擬的D4區(qū)域的近地面溫度場(chǎng)。由圖可見(jiàn),在陸地上,陸面覆蓋的復(fù)雜性引起了近地面氣溫分布的非均勻性,2種方案模擬近地面氣溫分布差別較大。Noah方案模擬的近地面氣溫分布較為合理,沿海溫度低,內(nèi)陸溫度高。同時(shí),Noah方案模擬出了白天寧波西部四明山地區(qū)由于谷地和山地的熱力對(duì)比形成溫度低于四周的事實(shí)。另外,Noah方案模擬出了城市熱島效應(yīng)(黃利萍等,2012,2013),即城區(qū)氣溫普遍高于周?chē)紖^(qū)的現(xiàn)象。而RUC方案模擬的近地面層氣溫呈現(xiàn)高溫分布,值普遍在34 ℃以上,沒(méi)有合理地反映出下墊面覆蓋的多樣性以及人為排熱對(duì)溫度分布的影響。

近地面氣溫的非均勻性引起了空氣密度和氣壓的區(qū)域性差異,并進(jìn)一步產(chǎn)生局地環(huán)流疊加在背景風(fēng)場(chǎng)之上。圖5c、d給出了24日12時(shí)的10 m風(fēng)場(chǎng)的分布,可見(jiàn),西部陸地上受偏西背景風(fēng)主導(dǎo),而東部海域上背景風(fēng)為偏南風(fēng)。在寧波市的西北岸和東北岸,2個(gè)方案都模擬出了比較明顯的海風(fēng)特征。海風(fēng)與背景風(fēng)在寧波北部沿岸一帶輻合。在輻合區(qū)上對(duì)應(yīng)著較強(qiáng)的抬升運(yùn)動(dòng)(圖略),配合一定的不穩(wěn)定層結(jié)和水汽條件,局地的弱對(duì)流系統(tǒng)逐漸發(fā)展起來(lái)。值得注意的是Noah方案不僅模擬出了海風(fēng)的特征,也模擬出了城鄉(xiāng)溫度差異引起的水平流場(chǎng),熱島流場(chǎng)與背景風(fēng)、海風(fēng)疊加并在城郊輻合,也進(jìn)一步促進(jìn)了抬升運(yùn)動(dòng)和對(duì)流的發(fā)展。另外,需要指出的是Noah方案在東南部海域上模擬的風(fēng)速明顯大于RUC方案,其原因有待進(jìn)一步研究。而RUC方案對(duì)城市下墊面特征的描述并不充分導(dǎo)致了近地面風(fēng)場(chǎng)的輻合抬升運(yùn)動(dòng)并沒(méi)有Noah方案強(qiáng)烈,這可能是RUC方案模擬的雷暴發(fā)展初期回波較弱的一個(gè)原因。

圖5　Noah(a,c)和RUC(b,d)方案模擬的寧波地區(qū)2006年6月24日12時(shí)的2 m溫度場(chǎng)(a,b;單位:℃)和10 m風(fēng)場(chǎng)(c,d;單位:m·s-1)(圖5a中多邊形區(qū)域?yàn)閷幉ㄊ械闹行某菂^(qū))Fig.5　(a,b)Air temperature(units:℃) at 2 m and (c,d)wind speed(units:m·s-1) at 10 m in Ningbo area at 12:00 BST 24 June 2006 simulated by (a,c)Noah and (b,d)RUC schemes(The polygon area in Fig.5a denotes the central city area of Ningbo city)

隨著海風(fēng)的移入,對(duì)流系統(tǒng)在輻合區(qū)上逐漸發(fā)展起來(lái)(圖4a、e),同時(shí)海風(fēng)從海面上攜帶了大量的水汽為后續(xù)對(duì)流的增強(qiáng)提供充沛水汽來(lái)源,這與Suresh(2007)的觀測(cè)結(jié)果相一致。之后,隨著北來(lái)的強(qiáng)對(duì)流系統(tǒng)逐漸登陸寧波,由于雷暴云中產(chǎn)生的下沉氣流到近地面時(shí)向四周散開(kāi),使得寧波北部近地面風(fēng)場(chǎng)結(jié)構(gòu)發(fā)生了較大的改變(圖4b、f)。風(fēng)場(chǎng)以強(qiáng)回波區(qū)為中心向四周輻散,輻散中心外圍風(fēng)速均大于10 m·s-1,甚至出現(xiàn)了18 m·s-1的大風(fēng),大風(fēng)和周?chē)谋尘帮L(fēng)疊加在回波附近形成更強(qiáng)烈的輻合,強(qiáng)烈的輻合引起的抬升促進(jìn)了回波的東南側(cè)新單體的生成,強(qiáng)度逐漸增強(qiáng)。同時(shí),雷暴大風(fēng)的覆蓋范圍也向東南擴(kuò)展。在移動(dòng)的過(guò)程中,強(qiáng)對(duì)流系統(tǒng)與東南岸的海風(fēng)相遇,使局地對(duì)流與移入對(duì)流在象山港西北側(cè)合并加強(qiáng)。由于RUC方案模擬的象山港北側(cè)的海風(fēng)強(qiáng)度較強(qiáng),RUC方案模擬的第三次合并過(guò)程早于Noah方案(圖4c、g)。由于太陽(yáng)輻射的消失引起海陸溫差的演變并且受雷暴系統(tǒng)的風(fēng)場(chǎng)掩蓋,18時(shí)之后海風(fēng)特征變得逐漸微弱,近地面主要由偏南風(fēng)主導(dǎo)。

為進(jìn)一步比較Noah方案和RUC方案模擬的風(fēng)場(chǎng)和溫度場(chǎng)的差異,分別對(duì)各時(shí)刻溫度場(chǎng)和風(fēng)場(chǎng)作區(qū)域內(nèi)站點(diǎn)平均并與觀測(cè)值對(duì)比(圖6a、b)。比較模擬與觀測(cè),2種方案均較為真實(shí)地模擬出了氣溫的日變化特征,但所模擬的最高氣溫出現(xiàn)的時(shí)間均提前了約1 h。不同陸面參數(shù)化方案對(duì)溫度量值的差異主要集中在白天:雷暴發(fā)生前期(09—12時(shí)),Noah方案所模擬的近地面氣溫與實(shí)際觀測(cè)較為吻合,而RUC方案模擬的近地面氣溫偏高(與圖5a、b的結(jié)論一致)。13時(shí)以后對(duì)流逐漸發(fā)展起來(lái),Noah方案所模擬的近地面氣溫偏低,這是由于Noah方案模擬的初始階段對(duì)流偏強(qiáng)形成更深厚的對(duì)流云遮蓋天空所致。雷暴發(fā)生前期(09—12時(shí)),模擬風(fēng)速與觀測(cè)風(fēng)速較為吻合,其差值約為0.3 m·s-1。此后(12—18時(shí)),兩個(gè)方案均模擬出由于雷暴的發(fā)生帶來(lái)的風(fēng)速增幅過(guò)程,然而模擬的風(fēng)速增大時(shí)間均較實(shí)況提前發(fā)生。從風(fēng)速日變化整個(gè)過(guò)程來(lái)看,Noah方案的模擬結(jié)果與觀測(cè)更為接近,而RUC方案與實(shí)際偏差稍微大一些。

圖6　不同方案模擬的D4區(qū)域內(nèi)站點(diǎn)平均值與實(shí)際觀測(cè)的站點(diǎn)平均值比較　　a.2 m溫度(單位:℃);b.10 m風(fēng)速(單位:m·s-1)Fig.6　Comparisons between the hourly site-averaged results simulated by different experiments and the observations in D4　　a.air temperature at 2 m(units:℃);b.wind speed at 10 m(units:m·s-1)

4.3　陸面參數(shù)化方案對(duì)降水時(shí)空分布特征的影響

圖7是24日雷暴過(guò)程降水最集中的12—20時(shí)實(shí)際觀測(cè)的累積降水量(圖7a)和模擬結(jié)果(圖7b、c)。實(shí)際觀測(cè)到的累積降水區(qū)呈南北向分布,降水大值中心位于(121.5°E,29.4°N)附近,中心值達(dá)到100 mm以上;在其北側(cè)有幾個(gè)雨量相對(duì)較小的中心,分別位于寧波東北沿岸及象山港上空。兩個(gè)方案模擬的降水強(qiáng)度與實(shí)際觀測(cè)相近,然而都高估了雷暴過(guò)程中50 mm以上累積降水量的范圍,且模擬的降水分布均較實(shí)際偏西。相比之下,Noah方案模擬的降水空間分布要優(yōu)于RUC方案:首先,Noah方案對(duì)降水中心位置和分布模態(tài)的表現(xiàn)較好,而RUC方案沒(méi)有模擬出象山港西岸的強(qiáng)降水中心,并且在象山港南北側(cè)模擬出了2個(gè)虛假的降水中心;另外,RUC方案低估了主城區(qū)降水。

圖7　2006年6月24日12—20時(shí)的累積降水量(陰影;單位:mm)　　a.觀測(cè)結(jié)果;b.Noah方案模擬結(jié)果;c.RUC方案模擬結(jié)果Fig.7　Accumulated precipitation from 12:00 BST to 20:00 BST 24 June 2006(shaded areas;units:mm)　　a.observed results;b.simulated results by Noah scheme;c.simulated results by RUC scheme

圖8給出了不同試驗(yàn)?zāi)M的區(qū)域平均降水強(qiáng)度隨時(shí)間的演變情況,觀測(cè)的降水強(qiáng)度由自動(dòng)氣象站觀測(cè)值平均而得。總體來(lái)看,觀測(cè)降水時(shí)間集中發(fā)生在15—20時(shí),最大平均降水強(qiáng)度達(dá)到5.25 mm·h-1。而兩個(gè)試驗(yàn)?zāi)M的降水持續(xù)時(shí)間則比觀測(cè)時(shí)間長(zhǎng)約2 h,這可能與模式初始條件中仍缺少對(duì)環(huán)境場(chǎng)中尺度特征的詳細(xì)描述有關(guān)(蒙偉光等,2012)。與觀測(cè)相比,兩個(gè)方案模擬的最大平均降水強(qiáng)度較小:Noah方案所模擬的最大降水強(qiáng)度約5 mm·h-1,與觀測(cè)更為接近;RUC方案的模擬的最大降水強(qiáng)度最小,只有3.5 mm·h-1,比實(shí)際觀測(cè)低了近30%。同時(shí),Noah方案所模擬的平均降水強(qiáng)度演變特征與觀測(cè)特征均呈現(xiàn)單峰形變化,而RUC方案呈現(xiàn)出雙峰形的變化特征。

4.4　陸面參數(shù)化方案對(duì)地表通量及熱動(dòng)力不穩(wěn)定參數(shù)的影響

性質(zhì)復(fù)雜的下墊面所構(gòu)成的陸地表面與大氣之間不斷進(jìn)行著各種時(shí)空尺度的相互作用以及能量、水汽等物質(zhì)的交換。盧萍和宇如聰(2008)指出強(qiáng)降水對(duì)于陸氣通量有一定的依賴(lài)性,其中潛熱通量和感熱通量的作用最大。因此,地表通量特征是否被合理地描述對(duì)短時(shí)強(qiáng)降水事件的模擬效果有較大的影響。

為探討2個(gè)方案模擬的降水空間分布差異的原因,分別對(duì)2個(gè)方案D4區(qū)域主要覆蓋類(lèi)型(水體,混合林,農(nóng)田,城市)的表面通量和熱動(dòng)力不穩(wěn)定參數(shù)取平均,圖9顯示了2個(gè)參數(shù)化方案不同地表覆蓋類(lèi)型的平均感熱通量,潛熱通量,凈輻射,邊界層高度以及不穩(wěn)定參數(shù)的逐時(shí)變化。首先,可以看出地表通量和熱動(dòng)力不穩(wěn)定參數(shù)受下墊面影響顯著;其次,在白天,無(wú)論哪種陸面覆蓋類(lèi)型,Noah方案模擬的感熱(潛熱)通量均大(小)于RUC方案;此外,在陸地表面(混合林,農(nóng)田,城市),潛熱通量和凈輻射對(duì)陸面方案的選取最為敏感,其次是感熱通量,最不敏感是邊界層高度。最后,需要指出的是,在水體表面,由于下墊面性質(zhì)比較穩(wěn)定、均一,參數(shù)化方案的選取對(duì)模擬結(jié)果的影響較大;而在陸地表面,由于下墊面屬性的多樣性和復(fù)雜性,實(shí)際模擬的降水差異是由下墊面本身屬性的差異與參數(shù)化方案的差異共同作用的結(jié)果。

圖9　D4區(qū)域不同方案、不同陸地覆蓋表面在2006年6月24日的平均感熱通量(a;單位:W·m-2)、平均潛熱通量(b;單位:W·m-2)、平均凈輻射(c;單位:W·m-2)、平均邊界層高度(d;單位:km)、平均對(duì)流有效位能(e;單位:J·kg-1)、平均對(duì)流抑制能(f;單位:J·kg-1)的逐時(shí)變化(N:Noah方案,R:RUC方案;紅色代表城市,綠色代表混合林,黃色代表農(nóng)田,藍(lán)色代表水體)Fig.9　Hourly-averaged (a)sensible heat flux(units:W·m-2),(b)latent heat flux(units:W·m-2),(c)net radiation (units:W·m-2),(d)boundary layer height(units:km),(e)convective available potential energy(units:J·kg-1),and (f)convective inhibition energy(units:J·kg-1) simulated by different experiments with different land cover categories in D4 on 24 June 2006(N:Noah scheme,R:RUC scheme;the red,green,yellow and blue colors denote urban,mixed forest,croplands and water bodies,respectively)

5結(jié)論

本文利用中尺度模式WRF分別選用不同的陸面過(guò)程參數(shù)化方案(Noah、RUC)對(duì)寧波地區(qū)一次夏季強(qiáng)對(duì)流天氣(雷暴)過(guò)程進(jìn)行了數(shù)值模擬試驗(yàn),重點(diǎn)分析了2種不同復(fù)雜程度的陸面過(guò)程參數(shù)化方案對(duì)于雷暴系統(tǒng)發(fā)生、發(fā)展過(guò)程以及近地面熱動(dòng)力條件的影響,試驗(yàn)結(jié)果表明:

1)Noah方案不僅模擬出了雷暴前期海陸、山地平原的熱力對(duì)比引起的近地面氣溫的不均勻分布,而且還模擬出了城市熱島效應(yīng);然而RUC方案沒(méi)有明顯地反映出陸面覆蓋的多樣性以及人為排熱對(duì)溫度分布的影響。

2)Noah方案模擬的風(fēng)場(chǎng)輻合較強(qiáng),有利于抬升運(yùn)動(dòng)和對(duì)流的啟動(dòng);由于考慮到了城市下墊面對(duì)氣流的拖拽作用,Noah方案所模擬的城郊之間風(fēng)速的差異比RUC方案更為明顯。

3)Noah方案模擬的雷暴啟動(dòng)、發(fā)展過(guò)程與觀測(cè)較為一致,而RUC方案較好地描述出了雷暴演變過(guò)程中的關(guān)鍵階段(三次合并過(guò)程)。

4)不同陸面過(guò)程參數(shù)化方案所模擬的雷暴累積降水范圍均較實(shí)況偏大,但Noah方案在對(duì)降水的走向、強(qiáng)度、降水中心位置的模擬方面具有一定優(yōu)勢(shì)。

5)模擬的雷暴持續(xù)時(shí)間對(duì)陸面過(guò)程參數(shù)化方案有很強(qiáng)的依賴(lài)性,兩個(gè)方案所模擬的雷暴過(guò)程持續(xù)時(shí)間均長(zhǎng)于實(shí)際雷暴持續(xù)時(shí)間,但程度有差異。

6)潛熱通量和凈輻射對(duì)陸面過(guò)程參數(shù)化方案的選擇最為敏感,其次是感熱通量,最不敏感的是邊界層高度。無(wú)論哪種陸面覆蓋類(lèi)型,白天Noah方案模擬的感熱通量均大于RUC方案,而Noah方案模擬的潛熱通量均小于RUC方案。

參考文獻(xiàn)(References)：

陳力強(qiáng),周小珊,楊森.2004.一次強(qiáng)風(fēng)暴的垂直環(huán)境特征數(shù)值模擬分析[J].氣象,30(9):3-8.Chen Liqiang,Zhou Xiaoshan,Yang Sen.2004.Numerical simulation and analysis of vertical environment character of a severe storm[J].Meteor Mon,30(9):3-8.(in Chinese).

陳思蓉,朱偉軍,周兵.2009.中國(guó)雷暴氣候分布特征及變化趨勢(shì)[J].大氣科學(xué)學(xué)報(bào),32(5):703-710.Chen Sirong,Zhu Weijun,Zhou Bing.2009.Climate characteristic and variation tendency of thunderstorm in China[J].Trans Atmos Sci,32(5):703-710.(in Chinese).

何宏讓,魏紹遠(yuǎn),曾文華.1999.江淮地區(qū)副高后部強(qiáng)對(duì)流的分析和數(shù)值模擬研究[J].氣象科學(xué),19(2):142-149.He Hongrang,Wei Shaoyuan,Zeng Wenhua.1999.Analysis and madeling study of severe thunderstorm behind a subtropical high over the Changjiang-Huaihe area[J].Scientia Meteorologica Sinca,19(2):142-149.(in Chinese).

黃利萍,苗峻峰,劉月琨.2012.天津城市熱島效應(yīng)的時(shí)空變化特征[J].大氣科學(xué)學(xué)報(bào),35(5):620-632.Huang Liping,Miao Junfeng,Liu Yuekun.2012.Spatial and temporal variation characteristics of urban heat island in Tianjin[J].Trans Atmos Sci,35(6):620-632.

黃利萍,苗峻峰,劉月琨,等.2013.天津地區(qū)夏季海陸風(fēng)對(duì)城市熱島日變化特征影響的觀測(cè)分析[J].大氣科學(xué)學(xué)報(bào),36(4):417-425.Huang Liping,Miao Junfeng,Liu Yuekun,et al.2013.Observational analysis of influence of sea-land breeze on diurnal characteristics of urban heat island in Tianjin during summer[J].Trans Atmos Sci,36(4):417-425.

李昀英,宇如聰,傅云飛,等.2008.一次熱對(duì)流降水成因的分析和模擬[J].氣象學(xué)報(bào),66(2):190-202.Li Yunying,Yu Rucong,Fu Yunfei,et al.2008.A case study on triggering of thermal convective precipitation[J].Acta Meteor Sinica,66(2):190-202.(in Chinese).

李志楠,康玉霞.1995.影響北京地區(qū)東北方路徑帶狀雷暴天氣分析[J].氣象,21(8):44-48.Li Zhinan,Kang Yuxia.1995.An analysis of northeastern route banded storm weather[J].Meteor Mon,21(8):44-48.(in Chinese).

盧萍,宇如聰.2008.地表潛熱通量對(duì)四川地區(qū)降水影響的數(shù)值分析[J].高原山地氣象研究,28(3):1-7.Lu Ping,Yu Rucong.2008.Numerical analysis on the impacts of surface latent heat flux transport on Sichuan rainfall process[J].Plateau and Mountain Meteorology Research,28(3):1-7.(in Chinese).

馬紅云,郭品文,宋潔.2009.耦合不同陸面方案的WRF模式對(duì)2007年7月江淮強(qiáng)降水過(guò)程的模擬[J].大氣科學(xué),33(3):557-567.Ma Hongyun,Guo Pinwen,Song Jie.2009.Simulation of “2007.7” heavy rainfall case in the Changjiang-Huaihe valley using the WRF model with different land surface schemes[J].Chinese J Atmos Sci,33(3):557-567.(in Chinese).

蒙偉光,李江南,王安宇,等.2005.凝結(jié)加熱和地表通量對(duì)華南中尺度對(duì)流系統(tǒng)(MCS)發(fā)生發(fā)展的影響[J].熱帶氣象學(xué)報(bào),21(4):368-376.Meng Weiguang,Li Jiangnan,Wang Anyu,et al.2005.Effects of condensation heating and surface fluxes on the development of a south China Mesoscale Convective System(MCS)[J].J Trop Meteor,21(4):368-376.(in Chinese).

蒙偉光,李昊睿,張艷霞,等.2012.珠三角城市環(huán)境對(duì)對(duì)流降水影響的模擬研究[J].大氣科學(xué),36(5):1063-1076.Meng Weiguang,Li Haorui,Zhang Yanxia,et al.2012.A modeling study of the impacts of Pearl River Delta urban environment on convective precipitation[J].Chinese J Atmos Sci,36(5):1063-1076.(in Chinese).

苗峻峰.2014.城市熱島和海風(fēng)環(huán)流相互作用的數(shù)值模擬研究進(jìn)展[J].大氣科學(xué)學(xué)報(bào),37(4):521-528.Miao Junfeng.2014.An overview of numerical studies of interaction of urban heat island and sea breeze circulations[J].Trans Atmos Sci,37(4):521-528.(in Chinese).

Miao J F,Chen D,Borne K.2007.Evaluation and comparison of Noah and Pleim-Xiu land surface models in MM5 using G?TE2001 data:Spatial and temporal variations in near-surface air temperature. J Appl Meteor Climatol,46:1587-1605.

Miao J F,Chen D,Wyser K,et al.2008.Evaluation of MM5 mesoscale model at local scale for air quality applications over the Swedish west coast:Influence of PBL and LSM parameterizations[J].Meteor Atmos Phys,99:77-103.

Miao J F,Wyser K,Chen D,et al.2009.Impacts of boundary layer turbulence and land surface process parameterizations on simulated sea breeze characteristics[J].Ann Geophys,27(6):2303-2320.

Pan H L,Mahrt L.1987.Interaction between soil hydrology and boundary-layer development[J].Bound-Lay Meteor,38:185-202.

Smirnova T G,Brown J M,Benjamin S G.1997.Performance of different soil model configurations in simulating ground surface temperature and surface fluxes[J].Mon Wea Rev,125(8):1870-1884.

宋潔慧,壽紹文,劉旭,等.2009.寧波一次典型夏季海陸風(fēng)過(guò)程觀測(cè)分析和數(shù)值模擬[J].熱帶氣象學(xué)報(bào),25(3):336-342.Song Jiehui,Shou Shaowen,Liu Xu,et al.2009.Observational studies and numerical simulation of a land and sea breezes case over eastern Ningbo[J].J Trop Meteor,25(3):336-342.(in Chinese).

Suresh R.2007.Observation of sea breeze front and its induced convection over Chennai in Southern Peninsular India using Doppler weather radar[J].Pure Appl Geophys.164:1511-1525.

Talbot C,Augustin P,Leroy C,et al.2007.Impact of a sea breeze on the boundary-layer dynamics and the atmospheric stratification in a coastal area of the North Sea[J].Bound-Lay Meteor,125(1):133-154.

汪雅,苗峻峰,談?wù)苊?2013.寧波地區(qū)海—陸下墊面差異對(duì)雷暴過(guò)程影響的數(shù)值模擬[J].氣象學(xué)報(bào),71(6):1146-1159.Wang Ya,Miao Junfeng,Tan Zhemin.2013.Numerical study of the impact of differences between sea and land underlying surface on thunderstorm in Ningbo area[J].Acta Meteor Sinica,71(6):1146-1159.(in Chinese).

王玉國(guó),吳增茂,常志清,等.2009.遼東灣西岸三維海陸風(fēng)特征的數(shù)值模擬[J].海洋預(yù)報(bào),26(3):12-21.Wang Yuguo,Wu Zengmao,Chang Zhiqing,et al.2009.Numerical simulation of three dimensional structure of sea-land breeze at west coast of Liao Dong Bay[J].Marine Forecasts,26(3):12-21.(in Chinese).

肖輝,吳玉霞,胡朝霞,等.2002.旬邑地區(qū)冰雹云的早期識(shí)別及數(shù)值模擬[J].高原氣象,21(2):159-166.Xiao Hui,Wu Yuxia,Hu Zhaoxia,et al.2002.Earlier identification and numerical simulation of hail stoms occurring in Xunyi region[J].Plateau Meteorology,21(2):159-166.(in Chinese).

Xiu A,Pleim J E.2001.Development of a land surface model.Part I:Application in a mesoscale meteorological model[J].J Appl Meteor,40:192-209.

楊薇,苗峻峰,談?wù)苊?2014.陸面過(guò)程參數(shù)化對(duì)太湖地區(qū)雷暴過(guò)程模擬的影響[J].大氣科學(xué)學(xué)報(bào),37(6):804-814.Yang Wei,Miao Junfeng,Tan Zhemin.2014.Impact of land surface parameterizations on simulation of a thunderstorm over Tai Lake region[J].Trans Atmos Sci,37(6):804-814.(in Chinese).

曾新民,吳志皇,宋帥,等.2012.WRF模式不同陸面方案對(duì)一次暴雨事件模擬的影響[J].地球物理學(xué)報(bào),55(1):16-28.Zeng Xinmin,Wu Zhihuang,Song Shuai,et al.2012.Effects of land surface schemes on the simulation of a heavy rainfall event by WRF[J].Chinese Journal of Geophysics,55(1):16-28.(in Chinese).

張春喜,朱佩君,鄭永光,等.2005.一次春季暴雨不穩(wěn)定條件和對(duì)流觸發(fā)機(jī)制的數(shù)值模擬研究[J].北京大學(xué)學(xué)報(bào):自然科學(xué)版,41(5):746-753.Zhang Chunxi,Zhu Peijun,Zheng Yongguang,et al.2005.A study to the trigger mechanism and unstable conditions of spring heavy rainfall with numerical simulation[J].Acta Scientiarum Naturalium Universitatis Pekinensis,41(5):746-753.(in Chinese).

張敏鋒,馮霞.1998.我國(guó)雷暴天氣的氣候特征[J].熱帶氣象學(xué)報(bào),14(2):156-162.Zhang Minfeng,Feng Xia.1998.A study on climatic features and anomalies of the thunderstorm in China[J].J Trop Meteor,14(2):156-162.(in Chinese).

(責(zé)任編輯：張福穎)

Impact of land surface parameterizations on simulated thunderstorm process over the Ningbo area

WANG Ya1,2,MIAO Jun-feng1,TAN Zhe-min3

(1.School of Atmospheric Sciences,NUIST,Nanjing 210044,China;2.Shanghai Meteorological Bureau,Shanghai 200030,China; 3.Key Laboratory of Mesoscale Severe Weather of Ministry of Education,Nanjing University,Nanjing 210093,China)

Abstract:The new generation mesoscale model WRF is used to simulate a typical thunderstorm process in Ningbo area on 24 June 2006,and the impacts of land surface parameterization schemes(Noah and RUC) on numerical simulation of the thunderstorm process are tested.The results show that,for surface thermodynamic field before the thunderstorm happened,the simulated results with Noah scheme are more close to the observations.However,with RUC scheme,the simulated results do not reflect the diversity of underlying surface coverage and the impacts of city underlying surface,and the difference between rural and urban areas is not obvious.The start and development processes simulated by Noah scheme are more reasonable,and the results simulated by RUC scheme describe the key stages in the evolution process of thunderstorm well,such as three merging processes.Because the considered factors and physical process parameterizations are different,the strength and center of accumulated precipitation simulated by Noah scheme are more consistent with the observations than those simulated by RUC scheme.The duration of thunderstorm is strongly sensitive to land surface process parameterization,and both durations simulated by the two schemes are longer than actual duration of the thunderstorm to different degree.Whatever the land cover category in Ningbo area,the sensible heat flux(latent heat flux) simulated by Noah scheme is greater than(less than) that simulated by RUC scheme during the day.

Key words:thunderstorm process;land surface process parameterization;Noah scheme;RUC scheme;numerical sensitivity

doi：10.13878/j.cnki.dqkxxb.20130310001

文章編號(hào)：1674-7097(2015)03-0299-11

中圖分類(lèi)號(hào)：P4

文獻(xiàn)標(biāo)志碼：A

通信作者：苗峻峰,博士,教授,博士生導(dǎo)師,研究方向?yàn)橹谐叨葦?shù)值模擬,miaoj@nuist.edu.cn.

基金項(xiàng)目:公益性行業(yè)(氣象)科研專(zhuān)項(xiàng)經(jīng)費(fèi)項(xiàng)目(GYHY201006004)

收稿日期：2013-03-10；改回日期：2013-05-04

汪雅,苗峻峰,談?wù)苊?2015.陸面過(guò)程參數(shù)化對(duì)寧波地區(qū)雷暴過(guò)程模擬的影響[J].大氣科學(xué)學(xué)報(bào),38(3):299-309.doi：10.13878/j.cnki.dqkxxb.20130310001.

Wang Ya,Miao Jun-feng,Tan Zhe-min.2015.Impact of land surface parameterizations on simulated thunderstorm process over the Ningbo area[J].Trans Atmos Sci,38(3):299-309.(in Chinese).