周括 冉令坤 蔡仁 屈濤 陳蕾

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

2 烏魯木齊市氣象局, 烏魯木齊830002

3 烏魯木齊市米東區(qū)氣象局, 烏魯木齊831499

4 貴州省氣象服務(wù)中心, 貴陽(yáng)550002

1 引言

新疆地處歐亞大陸中部，屬典型的大陸性干旱、半干旱氣候，但在有利的地形條件下，個(gè)別地區(qū)年降水量可以達(dá)到700 mm 以上，也會(huì)出現(xiàn)大到暴雨（謝澤明等, 2018）。近年來(lái)，西北干旱半干旱地區(qū)極端降水頻發(fā)（戴新剛等, 2007; Wang et al.,2013, 2017），新疆極端暴雨研究逐漸引起人們注意（馬淑紅和席元偉, 1997; 楊蓮梅等, 2011）。楊蓮梅和楊濤（2005）采用濕位渦對(duì)新疆暴雨分析發(fā)現(xiàn)高層干位渦能夠反映冷暖空氣活動(dòng)，而低層位渦與降水強(qiáng)弱有較好對(duì)應(yīng)關(guān)系。莊薇等（2006）利用雷達(dá)觀測(cè)資料研究了發(fā)生在烏魯木齊的強(qiáng)降水過程，并揭示出不同對(duì)流單體的三維風(fēng)場(chǎng)結(jié)構(gòu)。王江等（2015）利用鋒生函數(shù)在南疆暴雨中分析發(fā)現(xiàn)鋒生作用是造成降水的重要原因。曾勇等（2019）在新疆西部暴雨研究中提到天山迎風(fēng)坡持續(xù)時(shí)間較長(zhǎng)的輻合線是造成暴雨的主要系統(tǒng)。黃昕等（2021）在一次伊犁河谷暴雨研究中發(fā)現(xiàn)中亞低渦及地形輻合線激發(fā)了暴雨。新疆暴雨具有明顯的地理分布規(guī)律，大致是山區(qū)降水多于盆地，北疆降水多于南疆（馬淑紅和席元偉, 1997; 黃玉霞等, 2019）。21 世紀(jì)以來(lái)，南疆降雨頻次增加，降水量明顯上升（黃建平等, 2014）。整體而言，新疆強(qiáng)降水在水汽相對(duì)匱乏且地形復(fù)雜的條件下形成，其形成演變機(jī)制有待進(jìn)一步研究。

垂直運(yùn)動(dòng)是影響對(duì)流系統(tǒng)和暴雨發(fā)展的一個(gè)關(guān)鍵因素。為了研究影響垂直運(yùn)動(dòng)變化的原因，氣象學(xué)者發(fā)展了準(zhǔn)地轉(zhuǎn)ω方程（Holton, 1979），并用于診斷分析天氣尺度系統(tǒng)的發(fā)展演變（Strahl and Smith, 2001）。但對(duì)于次天氣尺度和中尺度系統(tǒng)，則需要考慮非地轉(zhuǎn)風(fēng)、非絕熱加熱過程和靜力穩(wěn)定度的水平變化等因素的影響。因此，更適用于中尺度系統(tǒng)研究的廣義垂直運(yùn)動(dòng)診斷方程逐漸發(fā)展（ Pauley and Nieman, 1992; R?is?nen, 1995,Rantanen et al., 2017）。Hoskins et al.（1978）由等壓坐標(biāo)系f平面近似的方程組得到準(zhǔn)地轉(zhuǎn)Q矢量散度驅(qū)動(dòng)的Omega 方程，其意義在于Q矢量散度作為垂直運(yùn)動(dòng)的唯一強(qiáng)迫項(xiàng)，解決了傳統(tǒng)垂直運(yùn)動(dòng)診斷方程中渦度平流的垂直梯度和溫度平流的拉普拉斯項(xiàng)之間的相互抵消作用。此后，國(guó)內(nèi)外學(xué)者進(jìn)一步豐富Q矢量的理論研究工作，相繼提出半地轉(zhuǎn)Q矢量（李柏和李國(guó)杰, 1997）、非地轉(zhuǎn)干Q矢量（Davies-Jones, 1991; 張興旺, 1999）、非地轉(zhuǎn)濕Q矢量（張興旺, 1998; 姚秀萍和于玉斌, 2000, 2001）等理論。Yang et al.（2007）將凝結(jié)概率函數(shù)導(dǎo)入Q矢量表達(dá)式中，拓展Q矢量的適用范圍。將Q矢量沿等位溫線（Keyser et al., 1992）和等高線（Jusem and Atlas, 1998）分解，所得分量對(duì)暴雨發(fā)展具有一定的指示意義。周小剛等（2011）研究指出Q矢量并不是有別于準(zhǔn)地轉(zhuǎn)ω方程的新理論，而是以一種新的形式去表達(dá)二級(jí)環(huán)流的強(qiáng)迫因子。

以往的垂直運(yùn)動(dòng)診斷方程主要建立在等壓坐標(biāo)系下，即滿足靜力平衡假定，而局地暴雨中對(duì)流活動(dòng)具有明顯的非靜力平衡特征，加之我國(guó)新疆地區(qū)交錯(cuò)分布的地形地貌，使得對(duì)流觸發(fā)過程更為復(fù)雜。基于以上考慮，本文發(fā)展了基于地形追隨坐標(biāo)的非地轉(zhuǎn)、非靜力的Q矢量散度驅(qū)動(dòng)的廣義垂直運(yùn)動(dòng)方程，針對(duì)2021 年6 月15～17 日南疆暴雨過程，診斷分析影響暴雨發(fā)生、發(fā)展的關(guān)鍵因素。

2 天氣過程

2021 年6 月15～17 日我國(guó)新疆南部地區(qū)發(fā)生一次暴雨天氣過程。和田國(guó)家站單日降水量達(dá)45.5 mm，打破當(dāng)?shù)赜袣庀笥涗浺詠?lái)單日降水記錄。南疆洛浦縣日降水量達(dá)86.2 mm，局部地區(qū)日降水量超過100 mm，達(dá)到大暴雨量級(jí)。

此次暴雨發(fā)生在有利的天氣背景下（圖1），2021 年6 月15 日00 時(shí)（協(xié)調(diào)世界時(shí)，下同），200 hPa 急流中心位于青海省以北地區(qū)，南疆處在急流入口區(qū)右側(cè)，盛行輻散氣流。500 hPa 中亞大槽向西加深并轉(zhuǎn)為橫槽，南疆處在槽前偏西氣流控制下。700 hPa 高壓系統(tǒng)在60°E 咸海以東發(fā)展加強(qiáng)，阿爾泰山東北部的低渦向其西南側(cè)伸展，高低值系統(tǒng)之間偏北風(fēng)增強(qiáng)，強(qiáng)烈的北風(fēng)沿著阿爾泰山和天山之間的縫隙進(jìn)入南疆地區(qū)，在昆侖山脈的阻擋下劇烈輻合。近地面層水汽通量顯示，暴雨中的水汽主要源自西伯利亞及其以北地區(qū)，由于昆侖山脈的阻擋作用，水汽在南疆地區(qū)堆積聚集。綜上所述，充足的水汽供應(yīng)和低層輻合、高層輻散的動(dòng)力配置為暴雨的發(fā)生、發(fā)展提供有利條件。

圖1 2021 年6 月15 日00 時(shí)（協(xié)調(diào)世界時(shí)，下同）（a）200 hPa 風(fēng)場(chǎng)（箭矢，單位：m s-1）和大于30 m s-1 風(fēng)速（填色，單位：m s-1）、（b）500 hPa 高度場(chǎng)（黑色等值線，單位：gpm）、（c）700 hPa 高度場(chǎng)（黑色等值線，單位：gpm）和風(fēng)場(chǎng)（箭矢，單位：m s-1）、（d）近地面水汽通量（箭矢，單位：g cm-1 hPa-1 s-1）及水汽通量散度（填色，單位：10-7 g cm-2 hPa-1 s-1）。圖c 中填色表示地形高度（單位：m）Fig. 1 (a) Wind field (arrows, units: m s-1) and wind speed more than 30 m s-1 (shadings, units: m s-1) at 200 hPa, (b) 500-hPa geopotential height(black contours, units: gpm), (c) geopotential height (black contours, units: gpm) and wind field (arrows, units: m s-1) at 700 hPa, (d) water vapor fluxes (arrows, units: g cm-1 hPa-1 s-1) and water vapor fluxes divergence (shadings, units: 10-7 g cm-2 hPa-1 s-1) at 0000 UTC on 15 June 2021. In Fig.c, the shaded areas denote the topographical height (units: m)

3 數(shù)值模擬

采用WRF 模式對(duì)此次暴雨天氣過程進(jìn)行數(shù)值模擬，模式背景場(chǎng)和側(cè)邊界來(lái)自NCEP（National Centers for Environmental Prediction）0.5°×0.5°全球預(yù)報(bào)系統(tǒng)（global forecast system，簡(jiǎn)稱GFS）分析場(chǎng)和預(yù)報(bào)場(chǎng)。模式水平分辨率為3 km，水平格點(diǎn)數(shù)為901×901，模式層氣壓頂為50 hPa。模擬起始時(shí)間為2021 年6 月13 日12 時(shí)，積分84 h，時(shí)間分辨率為1 h。采用WSM6云微物理方案（Hong and Lim, 2006）、RRTMG 長(zhǎng)波輻射和短波輻射方案（Iacono et al., 2008）、Noah 陸面模式（Tewari et al., 2004）、YSU 邊界層方案（Hong et al., 2006）。

本文觀測(cè)降水來(lái)自國(guó)家氣象信息中心日降水量產(chǎn)品（http://www.nmic.cn/data [2021-10-15]），2021年6 月14～15 日新疆中南部地區(qū)出現(xiàn)零星降水（圖2 左列）。15 日02 時(shí)降水系統(tǒng)逐漸在新疆南部的和田和喀什等地區(qū)加強(qiáng)。15 日09 時(shí)，降水系統(tǒng)發(fā)展旺盛，雨帶沿昆侖山脈呈東西走向。16 日00 時(shí)，降水系統(tǒng)向新疆西南部收縮并逐漸減弱。與觀測(cè)24 h 累計(jì)降水相比，16～17 日消亡階段，模擬降水強(qiáng)度偏強(qiáng)，但數(shù)值模式結(jié)果較好地還原了降水初生和發(fā)展階段的落區(qū)和強(qiáng)度，因此本文采用數(shù)值模式輸出數(shù)據(jù)對(duì)此次降水發(fā)生發(fā)展機(jī)制做進(jìn)一步分析。

圖2 2021 年6 月（a、b）15 日00 時(shí)、（c、d）16 日00 時(shí)、（e、f）17 日00 時(shí)實(shí)況（左列）與模擬（右列）的24 h 累計(jì)降水量（單位：mm）分布。紅色方框表示南疆暴雨區(qū)域Fig. 2 24-h accumulated precipitation (units: mm) observed (left column) and simulated (right column) at (a, b) 0000 UTC on 15 June 2021, (c, d)0000 UTC on 16 June 2021, (e, f) 0000 UTC on 17 June 2021. The red boxes denote rainstorm area in southern Xinjiang

4 地形追隨坐標(biāo)系垂直運(yùn)動(dòng)診斷方程

暴雨多發(fā)生在地形復(fù)雜地區(qū)，地形對(duì)暴雨有觸發(fā)和增幅作用，在暴雨動(dòng)力分析研究中需要考慮地形作用。此外，廣泛用于暴雨模擬和預(yù)報(bào)的中尺度WRF 模式采用地形追隨坐標(biāo)的非靜力平衡框架，避免了等壓面或等高面位于地形高度以下和垂直插值的誤差，那么利用WRF 模式面數(shù)據(jù)進(jìn)行計(jì)算和分析暴雨可能獲得一些新結(jié)果。垂直運(yùn)動(dòng)是影響暴雨發(fā)展的關(guān)鍵因素，也是暴雨機(jī)理研究重點(diǎn)方向之一。基于上述認(rèn)識(shí)，本文從WRF-ARW 模式質(zhì)量地形追隨坐標(biāo)控制方程出發(fā)，推導(dǎo)建立地形追隨坐標(biāo)的非靜力平衡廣義垂直運(yùn)動(dòng)方程，可以利用模式面數(shù)據(jù)診斷分析暴雨垂直運(yùn)動(dòng)發(fā)展機(jī)理。

WRF-ARW 模式質(zhì)量地形追隨坐標(biāo)定義為

采用Boussinesq 近似，經(jīng)整理得到質(zhì)量地形追隨坐標(biāo)非靜力平衡廣義垂直運(yùn)動(dòng)方程（具體推導(dǎo)過程詳見附錄A）：

將垂直運(yùn)動(dòng)方程中主要強(qiáng)迫項(xiàng)做如下分解：

表1 垂直運(yùn)動(dòng)方程右端主要強(qiáng)迫項(xiàng)物理意義Table 1 Physical meanings of the major forcing terms on right hand side of the vertical motion equation

需要說(shuō)明的是，廣義垂直運(yùn)動(dòng)方程（11）是垂直速度的診斷方程，描述哪些物理過程與當(dāng)前垂直運(yùn)動(dòng)狀態(tài)相關(guān)，在推導(dǎo)過程中采用水平動(dòng)量和熱力學(xué)方程，但沒用到垂直動(dòng)量方程。而WRF 模式垂直動(dòng)量方程描述的是哪些物理過程與垂直動(dòng)量個(gè)別變化有關(guān)。

5 暴雨垂直運(yùn)動(dòng)診斷分析

在2021 年6 月15 日02 時(shí)降水初生階段（圖3a），36.75°N 弱降水區(qū)上空為相當(dāng)位溫高值區(qū)，地形追隨坐標(biāo)下層結(jié)穩(wěn)定度表達(dá)式（Skamarock et al.,2021）為

其中， θe表示相當(dāng)位溫。弱降水區(qū)北側(cè)0.829 等η面以下為層結(jié)穩(wěn)定度負(fù)值區(qū)，即存在條件不穩(wěn)定層結(jié)，有利于垂直運(yùn)動(dòng)發(fā)展。此時(shí)弱降水區(qū)上空存在弱垂直運(yùn)動(dòng)，降水區(qū)北部偏北風(fēng)逐漸增強(qiáng)。受到地形影響，0.829～0.088 等η面之間的等相當(dāng)位溫線和風(fēng)場(chǎng)均呈現(xiàn)波狀形態(tài)。在15 日09 時(shí)降水成熟階段，37°N 強(qiáng)降水中心上空等相當(dāng)位溫線垂直伸展至0.626 等η面以上（圖3b），這種近乎中性層結(jié)分布特征促進(jìn)了垂直運(yùn)動(dòng)的增強(qiáng)。昆侖山北部偏北風(fēng)進(jìn)一步增強(qiáng)，受地形阻擋輻合上升。

圖3 2021 年6 月15 日（a）02 時(shí)、（b）09 時(shí)層結(jié)穩(wěn)定度（填色，單位：K m-1）、相當(dāng)位溫（藍(lán)色等值線，單位：K）和風(fēng)場(chǎng)（箭矢，單位：m s-1）沿78.35°E 剖面。黑色填色表示地形高度（單位：km），紅色實(shí)線表示模擬1 h 累計(jì)降水量Fig. 3 Cross sections of stratification stability (shadings, units: K m-1), equivalent potential temperature (blue contours, units: K), and wind field(arrows, units: m s-1) along 78.35°E at (a) 0200 UTC on 15 June 2021; (b) 0900 UTC on 15 June 2021. The black shadings denote terrain height (units:km). The red solid lines denote the simulated 1-h accumulated precipitation

利用WRF 模式輸出的模式面數(shù)據(jù)計(jì)算地形追隨坐標(biāo)系垂直運(yùn)動(dòng)診斷方程（11）右端各強(qiáng)迫項(xiàng)，以此來(lái)分析暴雨中影響垂直運(yùn)動(dòng)發(fā)展的主要因素。在2021 年6 月15 日02 時(shí)降水觸發(fā)階段，低層垂直運(yùn)動(dòng)較弱，受地形影響，在低層的等η面上存在氣壓和位勢(shì)的水平梯度，因此沿著昆侖山脈地區(qū)存在Q矢量散度異常高值區(qū)（圖4a）。0.829 等η面以下的Q矢量散度主要分布在地形起伏較大的區(qū)域，高層0.626～0.208 等η面之間的Q矢量散度異常區(qū)和垂直運(yùn)動(dòng)分布存在對(duì)應(yīng)性（圖4c），垂直運(yùn)動(dòng)波狀特征明顯。在15 日09 時(shí)降水成熟階段，主雨帶沿高地形呈西北—東南走向，Q矢量散度異常增強(qiáng)，并與降水分布對(duì)應(yīng)（圖4b）。（38.5°N,79.6°E）附近存在另一個(gè)降水區(qū)和Q矢量散度的異常分布。在空間分布上，0.46 等η面以下的Q矢量散度較降水初生階段顯著增強(qiáng)（圖4d）。37°N處的降水高值區(qū)上空垂直運(yùn)動(dòng)發(fā)展旺盛，并存在Q矢量散度異常與之對(duì)應(yīng)。高層0.46 等η面以上Q矢量散度同樣有所增強(qiáng)，相應(yīng)波狀垂直運(yùn)動(dòng)增強(qiáng)。由此可見，地形追隨坐標(biāo)系下的Q矢量散度能夠反映空間垂直運(yùn)動(dòng)狀態(tài)。

圖4 2021 年6 月15 日02 時(shí)（左）、15 日09 時(shí)（右）0.863 等η 面上（a、b）Q 矢量散度（填色，單位：10-6 s-3）疊加1 h 累計(jì)降水量（黑色等值線，單位：mm）水平分布，（c、d）Q 矢量散度（填色，單位：10-6 s-3）疊加垂直速度（黑色等值線，單位：m s-1）沿78.35°E 剖面Fig. 4 (a, b) Horizontal distribution of Q vector divergence (shadings, units: 10-6 s-3) on 0.863 η surface and 1-h accumulated precipitation (black contours, units: mm) at 0200 UTC on 15 June 2021; (c, d) cross sections of Q vector divergence (shadings, units: 10-6 s-3) and vertical velocity (black contours, units: m s-1) along 78.35°E at 0900 UTC on 15 June 2021

對(duì)比方程（11）右端三個(gè)強(qiáng)迫項(xiàng)的強(qiáng)度變化（圖5）可知，經(jīng)向氣壓梯度力項(xiàng)是主要強(qiáng)迫項(xiàng)，其變化趨勢(shì)基本能夠反映總Q矢量散度的變化（圖5a）。緯向氣壓梯度力項(xiàng)是次要強(qiáng)迫項(xiàng)，強(qiáng)度始終弱于經(jīng)向氣壓梯度力強(qiáng)迫項(xiàng)，非地轉(zhuǎn)風(fēng)強(qiáng)迫項(xiàng)貢獻(xiàn)最弱。在經(jīng)向氣壓梯度力強(qiáng)迫項(xiàng)強(qiáng)度增強(qiáng)時(shí)段，小時(shí)降水量相應(yīng)增強(qiáng)，二者相關(guān)系數(shù)為0.63，高于另外兩強(qiáng)迫項(xiàng)與降水的相關(guān)系數(shù)。在0.863 等η面上，經(jīng)向氣壓梯度力強(qiáng)迫項(xiàng)仍然是總Q矢量散度的最大貢獻(xiàn)項(xiàng)（圖5b）。經(jīng)向氣壓梯度力強(qiáng)迫的空間分布和總Q矢量散度空間分布十分接近（圖6），因此進(jìn)一步討論經(jīng)向氣壓梯度力強(qiáng)迫項(xiàng)中包含的關(guān)鍵物理過程。

圖5 Q 矢量散度（divQ）及其分量（divQ1、divQ2、divQ3）強(qiáng)度（單位：10-6 s-3）在站點(diǎn)（37.24°N, 78.35°E）的（a）整層累加、（b）0.863 等η 面上的時(shí)間演變。左上角的數(shù)字代表各項(xiàng)與小時(shí)降水量的相關(guān)系數(shù)，綠色虛線表示1 h 累計(jì)降水量，下同F(xiàn)ig. 5 Time series of Q vector divergence (divQ) and its component (divQ1, divQ2, divQ3) intensity (units: 10-6 s-3) at (37.24°N, 78.35°E): (a) Full level accumulation; (b) 0.863-η level. The values in the top left-hand corner denote correlation coefficients between each term (divQ, divQ1, divQ2,divQ3) and hourly precipitation. The green dashed lines denote the 1-h accumulated precipitation (units: mm), the same below. divQ1, divQ2, divQ3 represent forcing terms of zonal pressure gradient, forcing term of meridional pressure gradient, ageostrophic wind, respectively

圖6 2021 年6 月15 日（a）02 時(shí)、（b）09 時(shí)divQ2 項(xiàng)（填色，單位：10-6 s-3）、垂直速度（黑色等值線，單位：m s-1）沿78.35°E 剖面。紅色實(shí)線表示模擬1 h 累計(jì)降水量（單位：mm）Fig. 6 Cross sections of divQ2 term (shadings, units: 10-6 s-3) and vertical velocity (black contours, units: m s-1) along 78.35°E at (a) 0200 UTC, (b)0900 UTC on 15 June 2021. The red solid lines denote the simulated 1-h accumulated precipitation (units: mm)

經(jīng)向氣壓梯度力強(qiáng)迫項(xiàng)divQ2中包含了氣壓梯度力平流項(xiàng)（17～19 式）、氣壓梯度力耦合水平風(fēng)切變作用項(xiàng)（20～22 式）、非絕熱加熱梯度項(xiàng)（23 式）和水物質(zhì)強(qiáng)迫耦合氣壓梯度力作用項(xiàng)（24 式）。分析發(fā)現(xiàn)，divQ2項(xiàng)中經(jīng)向氣壓梯度力耦合經(jīng)向風(fēng)切變項(xiàng)divQ2_5（23 式）是影響divQ2強(qiáng)度變化的主要貢獻(xiàn)項(xiàng)（圖7a）。垂直氣壓梯度力耦合緯向風(fēng)切變項(xiàng)divQ2_6（22 式）對(duì)divQ2強(qiáng)度變化的影響稍弱于divQ2_5項(xiàng)。在2021 年6 月15 日02 時(shí)至16 日06 時(shí)降水發(fā)展期間，非絕熱加熱梯度項(xiàng)divQ2_7對(duì)divQ2強(qiáng)度變化影響增強(qiáng)，divQ2_7與小時(shí)降水的相關(guān)系數(shù)達(dá)到0.62。其余各項(xiàng)對(duì)divQ2強(qiáng)度變化影響較小。divQ2各分量在0.863 等η面上的演變規(guī)律也有類似特點(diǎn)（圖7b）。

圖7 divQ2 項(xiàng)及其分量（divQ2_1、divQ2_2、……、divQ2_8）強(qiáng)度（單位：10-6 s-3）在（37.24°N, 78.35°E）的（a）整層累加、（b）0.863 等η 面上的時(shí)間演變Fig. 7 Time series of divQ2 term and its component (divQ2_1, divQ2_2,…, divQ2_8) intensity (units: 10-6 s-3) at (37.24°N, 78.35°E): (a) Full level accumulation; (b) 0.863-η level. divQ2_1, divQ2_2, …, and divQ2_8 represent vertical advection term of meridional pressure gradient force,meridional advection term of meridional pressure gradient force, zonal advection term of meridional pressure gradient force, zonal pressure gradient force coupled with zonal wind shear term, meridional pressure gradient force coupled with meridional wind shear term, vertical pressure gradient force coupled with zonal wind shear term, diabatic heating gradient term, and water mixing ratio coupled with pressure gradient term

進(jìn)一步分析發(fā)現(xiàn)經(jīng)向氣壓梯度力耦合經(jīng)向散度項(xiàng)divQ2_5_1（25 式）是divQ2_5項(xiàng)的主要組成項(xiàng)，二者的演變規(guī)律十分接近（圖8a）。對(duì)于divQ2_6項(xiàng)，兩個(gè)分量在15 日02 時(shí)前演變特征相近，在隨后的降水系統(tǒng)發(fā)展期間，垂直氣壓梯度力耦合緯向散度項(xiàng)divQ2_6_2（28 式）成為divQ2_6項(xiàng)的主要貢獻(xiàn)項(xiàng)（圖8b）。divQ2_7項(xiàng)中包含了非絕熱加熱經(jīng)向和垂直梯度的作用，其中非絕熱加熱經(jīng)向梯度項(xiàng)divQ2_7_1（29 式）在降水系統(tǒng)演變過程中起主導(dǎo)作用（圖8c）。

圖8 整層累加的（a）divQ2_5 項(xiàng)及其分量（divQ2_5_1、divQ2_5_2）強(qiáng) 度（單 位：10-6 s-3）、（b）divQ2_6 項(xiàng) 及 其 分 量（divQ2_6_1、divQ2_6_2）強(qiáng) 度（單 位：10-6 s-3）、（c）divQ2_7 項(xiàng) 及 其 分 量（divQ2_7_1、divQ2_7_2）強(qiáng) 度（單 位：10-6 s-3）在（37.24°N,78.35°E）的時(shí)間演變Fig. 8 Time series of full level cumulative intensity (units: 10-6 s-3) at(37.24°N, 78.35°E): (a) divQ2_5 term and its components (divQ2_5_1,divQ2_5_2); (b) divQ2_6 term and its components (divQ2_6_1, divQ2_6_2);(c) divQ2_7 term and its components (divQ2_7_1, divQ2_7_2). divQ2_5_1 and divQ2_5_2 represent meridional pressure gradient force coupled with meridional divergence term and meridional pressure gradient force coupled with vertical shear of meridional wind term, respectively.divQ2_6_1 and divQ2_6_2 represent vertical pressure gradient force coupled with meridional shear of zonal wind term and vertical pressure gradient force coupled with zonal divergence term, respectively.divQ2_7_1 and divQ2_7_2 represent meridional gradient term of diabatic heating and vertical gradient term of diabatic heating, respectively

受地形因素的影響，氣壓和位勢(shì)高度在等η面上水平梯度較大，并且十分穩(wěn)定。對(duì)于divQ2_5_1項(xiàng)而言，其變化趨勢(shì)主要由經(jīng)向散度和經(jīng)向氣壓梯度決定。如圖9a、c 所示，逐漸增強(qiáng)的偏北風(fēng)在37°N 昆侖山脈迎風(fēng)坡輻合，激發(fā)垂直上升運(yùn)動(dòng)。0.626 等η面以下地形陡峭區(qū)域存在顯著的經(jīng)向氣壓梯度，由于經(jīng)向氣壓梯度穩(wěn)定少變，divQ2_5_1項(xiàng)的變化主要由經(jīng)向散度引起，而經(jīng)向氣壓梯度增強(qiáng)了該項(xiàng)變化，因此divQ2_5_1項(xiàng)成為影響Q矢量散度變化的關(guān)鍵因素之一。在36.1°～36.8°N 之間，由于陡峭地形導(dǎo)致0.829 等η面以上產(chǎn)生波狀分布的經(jīng)向散度（圖9a），即存在地形慣性重力波活動(dòng)。重力波活動(dòng)導(dǎo)致昆山脈山區(qū)垂直運(yùn)動(dòng)發(fā)展，產(chǎn)生降水，并在37°N 山脈迎風(fēng)坡0.829 等η面以上產(chǎn)生高空輻散場(chǎng)，進(jìn)一步增強(qiáng)了主雨帶上空垂直上升運(yùn)動(dòng)。

緯向散度、垂直氣壓梯度和經(jīng)向位勢(shì)梯度是影響divQ2_6_2項(xiàng)的主要因素。由于氣壓和位勢(shì)梯度穩(wěn)定少變，緯向散度決定了divQ2_6_2的變化趨勢(shì)。偏東氣流的增強(qiáng)使得（38.5°N, 79.6°E）處的雨區(qū)東北部以及昆侖山脈以北部分區(qū)域存在輻合中心（圖9d），強(qiáng)迫垂直運(yùn)動(dòng)的發(fā)展。37.1°N 低層存在緯向散度的輻合區(qū)（圖9b），高層0.626 等η面以上存在波動(dòng)活動(dòng)導(dǎo)致的緯向散度輻散中心。相比經(jīng)向散度（圖9a），緯向散度強(qiáng)度較弱，對(duì)垂直運(yùn)動(dòng)的強(qiáng)迫作用弱于經(jīng)向散度。

圖9 2021 年6 月15 日09 時(shí)（a）經(jīng)向散度（填色，單位：10-4 s-1）、氣壓（紅色細(xì)實(shí)線，單位：hPa）、垂直速度（黑色等值線，單位：m s-1），（b）緯向散度（填色，單位：10-4 s-1）、位勢(shì)高度（藍(lán)色細(xì)實(shí)線，單位：gpm）、垂直速度（黑色等值線，單位：m s-1）沿78.35°E 剖面（紅色實(shí)線表示模擬1 h 累計(jì)降水量）；2021 年6 月15 日09 時(shí)0.981 等η 面上（c）經(jīng)向散度（填色，單位：10-4 s-1）、（d）緯向散度（填色，單位：10-4 s-1）疊加風(fēng)場(chǎng)（箭矢，單位：m s-1）、1 h 累計(jì)降水量（黑色等值線，單位：mm）水平分布Fig. 9 Cross sections of vertical velocity (black contours, units: m s-1) and (a) meridional divergence (shadings, units: 10-4 s-1) and pressure (red thin contours, units: hPa), (b) zonal divergence (shadings, units: 10-4 s-1) and geopotential height (blue thin contours, units: gpm) along 78.35°E at 0900 UTC on 15 June 2021. The red solid lines denote the simulated 1-h accumulated precipitation. Horizontal distributions of wind field (arrows, units: m s-1), 1-h accumulated precipitation (black contours, units: mm), and (c) meridional divergence (shadings, units: 10-4 s-1), (d) zonal divergence(shadings, units: 10-4 s-1) on 0.981 η surface at 0900 UTC on 15 June 2021

非絕熱加熱經(jīng)向梯度項(xiàng)在降水發(fā)展期間逐漸增強(qiáng)，并且與小時(shí)降水間存在較好相關(guān)性（圖8c）。由（29）式可知，非絕熱加熱經(jīng)向梯度項(xiàng)中的氣壓強(qiáng)迫項(xiàng)Fp主要受微物理過程水汽相變導(dǎo)致的潛熱加熱FΘ和 水汽源匯項(xiàng)Fqv的影響。降水過程中，水汽相變以及水汽輻合增強(qiáng)，使得Fqv增強(qiáng)（圖10b）。微物理過程釋放潛熱加熱大氣（圖10a），造成大氣熱力不均勻分布狀態(tài)，促進(jìn)垂直上升運(yùn)動(dòng)的發(fā)展。水汽以及潛熱釋放增強(qiáng)導(dǎo)致非絕熱加熱經(jīng)向梯度項(xiàng)（divQ2_7_1）增大。

圖10 2021 年6 月15 日09 時(shí)（a）微物理過程釋放的潛熱（填色，單位：K s-1）、（b）水汽混合比傾向（填色，單位：g g-1 s-1）疊加垂直速度（黑色等值線，單位：m s-1）沿78.35°E 剖面。紅色實(shí)線表示模擬1 h 累計(jì)降水量Fig. 10 Cross sections of vertical velocity (black contours, units: m s-1) and (a) latent heat (shadings, units: K s-1) from microphysical processes, (b)tendency (shadings, units: g g-1 s-1) of water vapor mixing ratio along 78.35°E at 0900 UTC on 15 June 2021. The red solid lines denote the simulated 1-h accumulated precipitation

水汽和大氣熱力狀態(tài)是影響divQ2_7_1項(xiàng)的主要因素，而廣義位溫能夠綜合表現(xiàn)大氣中水汽和熱力狀態(tài)，結(jié)合經(jīng)向散度與緯向散度的時(shí)間演變（圖11）可以發(fā)現(xiàn)，在15 日02 時(shí)之前，經(jīng)向散度以輻合作用為主，促進(jìn)垂直運(yùn)動(dòng)發(fā)展，緯向散度輻合強(qiáng)度較弱，并且廣義位溫較低，因此經(jīng)向輻合增強(qiáng)是觸發(fā)降水的主要因素。在15 日02 時(shí)至16日06 時(shí)，廣義位溫升高，表明降水區(qū)水汽和熱力集中的特點(diǎn)，經(jīng)向散度和緯向散度在此期間輻合作用增強(qiáng)，水平氣流輻合與非絕熱加熱過程共同促進(jìn)垂直運(yùn)動(dòng)發(fā)展，進(jìn)而增強(qiáng)降水。

圖11 0.981 等η 面上經(jīng)向散度（紅色實(shí)線，單位：10-4 s-1）、緯向散度（藍(lán)色實(shí)線，單位：10-4 s-1）和廣義位溫（黑色虛線，單位：K）在（37.24°N, 78.35°E）的時(shí)間演變Fig. 11 Time series of meridional divergence (red solid line, units: 10-4 s-1), zonal divergence (blue solid line, units: 10-4 s-1), and generalized potential temperature (black dashed line, units: K) on 0.981-η level at point (37.24°N, 78.35°E)

結(jié)合以上分析得到此次南疆暴雨發(fā)生、發(fā)展機(jī)制的概念圖（圖12）。60°E 咸海以東高壓系統(tǒng)和阿爾泰山東北部的低渦發(fā)展增強(qiáng)，高低值系統(tǒng)之間的偏北風(fēng)沿地形狹口灌入南疆地區(qū)。偏北風(fēng)在昆侖山脈大地形的阻擋下輻合，經(jīng)向散度輻合增強(qiáng)，配合經(jīng)向氣壓梯度，造成經(jīng)向氣壓梯度力耦合經(jīng)向散度項(xiàng)（divQ2_5_1）異常增強(qiáng)，激發(fā)了垂直上升運(yùn)動(dòng)，進(jìn)而導(dǎo)致降水發(fā)生。降水發(fā)展過程中，經(jīng)向輻合作用進(jìn)一步增強(qiáng)，同時(shí)緯向散度輻合增強(qiáng)，造成垂直氣壓梯度力耦合緯向散度項(xiàng)（divQ2_6_2）異常增強(qiáng)，強(qiáng)迫垂直運(yùn)動(dòng)增強(qiáng)，促進(jìn)降水發(fā)展。水汽的輻合作用以及云微物理過程釋放的熱量導(dǎo)致非絕熱加熱經(jīng)向梯度項(xiàng)（divQ2_7_1）強(qiáng)度增加，這一過程在降水發(fā)展及成熟期均維持一定強(qiáng)度。與此同時(shí)，偏西過山氣流在大地形的阻擋下激發(fā)出慣性重力波，有利于山區(qū)對(duì)流系統(tǒng)發(fā)展。重力波在下游山脈迎風(fēng)坡造成高層輻散，有利于垂直運(yùn)動(dòng)的增強(qiáng)和維持。在經(jīng)向和緯向輻合增強(qiáng)、非絕熱加熱過程以及地形重力波等多個(gè)因素作用下促成了這次罕見的南疆暴雨。

圖12 南疆暴雨發(fā)展演變概念圖Fig. 12 Schematic of the evolution of southern Xinjiang rainstorm

6 結(jié)論與討論

2021 年6 月15～17 日新疆南部地區(qū)發(fā)生了一次暴雨天氣過程，多地打破降水歷史極值。利用中尺度WRF 模式對(duì)此次過程進(jìn)行高分辨率數(shù)值模擬，模式結(jié)果較好地還原了暴雨的發(fā)展演變過程。暴雨主要發(fā)生在昆侖山脈北坡，為考慮地形對(duì)降水發(fā)生、發(fā)展的作用，本文從地形追隨坐標(biāo)控制方程出發(fā)推導(dǎo)建立非靜力平衡廣義垂直運(yùn)動(dòng)方程，結(jié)合模式面數(shù)據(jù)診斷暴雨中垂直運(yùn)動(dòng)發(fā)展的機(jī)理，得到如下初步結(jié)論：

（1）采用WRF 模式質(zhì)量地形追隨坐標(biāo)系控制方程，并取Boussinesq 近似，推導(dǎo)了非靜力平衡廣義垂直運(yùn)動(dòng)方程。方程左端為包含垂直速度的算子項(xiàng)，右端不顯含垂直速度，右端三項(xiàng)分別表示緯向氣壓梯度力強(qiáng)迫、經(jīng)向氣壓梯度力強(qiáng)迫和非地轉(zhuǎn)風(fēng)的強(qiáng)迫作用。右端強(qiáng)迫項(xiàng)中除了隱含了傳統(tǒng)Q矢量散度外，還能夠體現(xiàn)地形動(dòng)力強(qiáng)迫作用，因此，方程能夠用于南疆復(fù)雜地形下暴雨垂直運(yùn)動(dòng)的診斷分析。

（2）在降水觸發(fā)階段，低層Q矢量散度主要分布在地形起伏較大地區(qū)，高層Q矢量散度異常和垂直運(yùn)動(dòng)對(duì)應(yīng)較好。在降水成熟階段，低層垂直上升運(yùn)動(dòng)顯著增強(qiáng)，與之對(duì)應(yīng)Q矢量散度異常增強(qiáng)，高層Q矢量散度分布和波狀結(jié)構(gòu)的垂直運(yùn)動(dòng)對(duì)應(yīng)。地形追隨坐標(biāo)下的Q矢量散度能夠反應(yīng)垂直運(yùn)動(dòng)的狀態(tài)。

（3）經(jīng)向氣壓梯度力耦合經(jīng)向散度項(xiàng)（divQ2_5_1）、垂直氣壓梯度力耦合緯向散度項(xiàng)（divQ2_6_2）和非絕熱加熱經(jīng)向梯度項(xiàng)（divQ2_7_1）是影響暴雨垂直運(yùn)動(dòng)發(fā)展演變的主要強(qiáng)迫項(xiàng)。divQ2_5_1 項(xiàng)體現(xiàn)了增強(qiáng)的北風(fēng)在昆侖山脈大地形影響下的輻合作用，地形強(qiáng)迫抬升激發(fā)了對(duì)流活動(dòng)。與此同時(shí)，南疆地區(qū)偏東風(fēng)的增強(qiáng)使得緯向輻合增強(qiáng)，促進(jìn)垂直運(yùn)動(dòng)發(fā)展，其強(qiáng)度弱于經(jīng)向輻合的作用。在等η面上，緯向和垂直方向氣壓梯度增強(qiáng)了divQ2_5_1項(xiàng)和divQ2_6_2 項(xiàng)的強(qiáng)度變化，但地形作用下的氣壓梯度較為穩(wěn)定，兩個(gè)強(qiáng)迫項(xiàng)變化主要由水平散度決定。隨著降水發(fā)展，水汽輻合增強(qiáng)以及云微物理過程釋放的熱量使得非絕熱加熱經(jīng)向梯度增加，進(jìn)一步促進(jìn)垂直運(yùn)動(dòng)的發(fā)展。此外，偏西過山氣流在地形的作用下激發(fā)形成慣性重力波，主雨帶上空在波動(dòng)的影響下存在輻散氣流，對(duì)流活動(dòng)在經(jīng)向和緯向上輻合增強(qiáng)、非絕熱加熱過程以及地形重力波等多個(gè)因素作用下劇烈發(fā)展，造成了此次南疆極端暴雨。

本文利用質(zhì)量地形追隨坐標(biāo)的垂直運(yùn)動(dòng)方程診斷了南疆暴雨過程中影響垂直運(yùn)動(dòng)發(fā)展的因素，方程能夠反映復(fù)雜地形影響下垂直運(yùn)動(dòng)狀態(tài)，垂直運(yùn)動(dòng)方程建立在Boussinesq 近似假設(shè)下，其診斷效果還需要更多個(gè)例的檢驗(yàn)。

附錄A