1.Key Laboratory of Mesoscale Severe Weather/Ministry of Education of China, and School of Atmospheric Sciences, Nanjing University, Nanjing 210093, China 2.Institute of Meteorology and Oceanography, National University of Defense Technology, Nanjing 211101, China Manuscript received: 2017-12-19 Manuscript revised: 2018-05-16 Manuscript accepted: 2018-07-31 Abstract:In this study, a three-dimensional mesoscale model was used to numerically simulate the well-known "98.7" heavy rainfall event that affected the Yangtze Valley in July 1998. Two experiments were conducted to analyze the impact of moist processes on the development of meso-β scale vortices (MβV) and their triggering by mesoscale wind perturbation (MWP). In the experiment in which the latent heat feedback (LHF) scheme was switched off, a stable low-level col field (i.e., saddle field——a region between two lows and two highs in the isobaric surface) formed, and the MWP triggered a weak MβV. However, when the LHF scheme was switched on as the MWP was introduced into the model, the MβV developed quickly and intense rainfall and a mesoscale low-level jet (mLLJ) were generated. The thickness of the air column and average temperature between 400 and 700 hPa decreased without the feedback of latent heat, whereas they increased quickly when the LHF scheme was switched on, with the air pressure falling at low levels but rising at upper levels. A schematic representation of the positive feedbacks among the mesoscale vortex, rainfall, and mLLJ shows that in the initial stage of the MβV, the MWP triggers light rainfall and the latent heat occurs at low levels, which leads to weak convergence and ageostrophic winds. In the mature stage of the MβV, convection extends to the middle-to-upper levels, resulting in an increase in the average temperature and a stretching of the air column. A low-level cyclonic circulation forms under the effect of Coriolis torque, and the mLLJ forms to the southeast of the MβV. Keywords: moist process, latent heat feedback, mesoscale vortex, col field 摘要:中纬度地区鞍型场常常产生中尺度涡旋,造成暴雨甚至特大暴雨。本文利用三维中尺度模式对“98.7”鄂东特大暴雨进行了数值模拟。设计两个试验分析了湿物理过程对中尺度风场扰动激发β中尺度涡旋以及β中尺度涡旋发展的影响。 在关闭凝结潜热加热反馈方案试验中,鄂东地区对流层低层形成一个稳定的鞍型场,中尺度风场扰动在鞍型场中激发出一个弱的β中尺度涡旋。在另一个试验中,当加入中尺度风场扰动后,打开凝结潜热加热反馈方案,β中尺度涡旋迅速发展,产生强降水和中尺度低空急流。400hPa和700hPa之间气柱的厚度和平均温度在没有凝结潜热加热反馈时随积分减小,而在凝结潜热加热反馈作用下,却迅速加大,伴有低层气压下降,高层气压升高。概括了中尺度涡旋、降水和中尺度低空急流之间正反馈作用的概念模型:β中尺度涡旋形成初始阶段,中尺度风场扰动激发弱降水,凝结潜热加热出现在低层,产生弱辐合和非地转风。在β中尺度涡旋成熟阶段,对流发展到中高层,导致气柱平均温度升高,气柱变厚,在科氏力作用下,低层形成气旋性环流,在β中尺度涡旋东南侧形成中尺度低空急流。 关键词:湿物理过程, 凝结潜热加热反馈, 中尺度涡旋, 鞍型场
HTML
--> --> --> -->
3.1. Description of the "98.7" heavy rainfall event in China
The "98.7" heavy rainfall event is well known all over China. During 21-22 July 1998, heavy rainfall fell along the Yangtze River Valley in eastern Hubei Province. The accumulated precipitation was 633 mm at the Parrot Island station in Wuhan, and 501.8 mm of rain fell in Huangshi over the two days, leading to serious flooding across the Yangtze basin. It has been suggested that this heavy rainfall was directly triggered by a meso-β scale weather system (Bei et al., 2002; Jiang and Wang, 2012). Figure 1 shows the 500-hPa synoptic map and 700-hPa streamlines at 0000 UTC 21 July 1998 based on the 0.5°× 0.5° resolution GAME (GEWEX Asian Monsoon Experiment; GEWEX: Global Energy and Water Cycle Experiment) reanalysis data provided by the Japan Meteorological Agency and the Earth Observation Research Center/National Space of Development Agency of Japan. The precipitation dataset was obtained from the National Meteorological Center of China. The southwest vortex (SWV), which was located over (26°-32°N, 102°-112°E), was so strong that it is very clear at 500 hPa. Wuhan was located to the east of the SWV, where the winds were relatively light. A col field is evident between the SWV and the vortex in Anhui Province at 700 hPa, and Wuhan is close to the col point. The southwest flows to the south of the col field brought warm and moist air to the region of heavy rainfall. Figure1. Synoptic maps from (a) 500 hPa and (b) 700 hPa at 0000 UTC 21 July 1998. The solid lines denote geopotential height (contoured every 20 gpm); the dashed lines denote temperature (contoured every 2°C); full and half barbs represent 4 and 2 m s-1, respectively; the solid lines with arrows denote streamlines; the shading shows the 48-h accumulated precipitation (mm) from 0000 UTC 21 to 0000 UTC 23 July 1998; "▲" marks the location of Wuhan in Hubei Province; and the ellipse frame denotes the 700-hPa col field.
2 3.2. EXP1 -->
3.2. EXP1
Figure 2 shows the evolution of the 700-hPa streamlines and vorticity in EXP1 and EXP2. In EXP1, the MWPs induce a mesoscale positive vertical relative vorticity (hereafter, vorticity) region with a maximum vorticity greater than 16× 10-5 s-1near the col point at t=17 h. The maximum vorticity then falls rapidly to about 8× 10-5 s-1 at t=19 h, and to about 4× 10-5 s-1 at t=21 h. The vorticity is then maintained at the same intensity until t=27 h, but the maximum value falls to less than 4× 10-5 s-1 at t=30 h. Under the influence of MWPs, a dynamic MβV forms near the col point at t=21 h. The MβV is so weak that it persists for less than 2 h. This suggests that MWPs are able to induce a weak MβV. However, the MβV is short-lived without the LHF. Figure2. 700-hPa streamlines and vorticity (units: 10-5 s-1) in EXP1 at (a) t=17 h, (b) t=19 h and (c) t=21 h, and in EXP2 at (d) t=19 h, (e) t=21 h and (f) t=23 h. The shading denotes vorticity greater than 4× 10-5 s-1.
2 3.3. EXP2 -->
3.3. EXP2
As in EXP1, the positive vorticity also decreases quickly as the perturbation is introduced into the model (Figs. 2d-f). The maximum vorticity fell to about 8× 10-5 s-1 at t=19 h because of the weak rainfall and associated weak LHR before t=19 h. Figure 3 shows the evolution of precipitation intensity within (30°-31°N, 113.5°-114.5°E) in EXP1 and EXP2. The region covers the 48-h accumulated precipitation center near Wuhan in EXP2. The precipitation amount is extremely small in EXP1 (without LHF), whereas it increases rapidly after the LHF scheme is switched on in EXP2. With the intensification of the rainfall in EXP2, the vorticity increases after t=21 h. The MβV forms at t=19 h, which is earlier than in EXP1 by about 2 h. When the rainfall intensifies, the MβV also develops, with a maximum vorticity greater than 12× 10-5 s-1 at t=21 h and 20× 10-5 s-1 at t=23 h. Our simulation results show that the LHF must play an important role in the development of the MβV. Figure3. 700-hPa vorticity (shaded; contoured every 3× 10-5 s-1), divergence (black lines; contoured every 2× 10-5 s-1), wind vectors (green arrows; units: m s-1), and 2.5-min accumulated precipitation (red lines; every 0.03 mm), at (a) 0 min, (b) 5 min, (c) 10 min and (d) 30 min, from t=17 h, in EXP2.
Figure 4 shows the evolution of the 700-hPa winds, vorticity, divergence, and 2.5-min accumulated precipitation in EXP2. The southwesterly MWPs result in the formation of the vorticity and a divergence couplet. The vorticity maintains its intensity until t=30 min, with a maximum vorticity greater than 12× 10-5 s-1. This vorticity is the source of the weak MβV formation. Although the divergence couplet decreases rapidly, it triggers the deep mesoscale vertical circulation (not shown). There are two inverse circulations above and below 800 hPa. The 2.5-min accumulated precipitation patterns show that the rainfall occurs mainly in the 700-hPa divergence region, but not in the convergence region. The water vapors transported by the southwesterly flows condense after being lifted by the updrafts associated with the shallow vertical circulation below 800 hPa. Figure4. Evolution of average precipitation intensity (mm h-1) in region of 30°-31°N, 113.5°-114.5°E in EXP1 (gray line) and EXP2 (black line).
2 3.4. Comparative analysis -->
3.4. Comparative analysis
Figure 5 shows the differences in the 700-hPa wind speed and geopotential height between EXP2 and EXP1. After 1 h of integration (i.e., t=18 h) as the MWP is introduced into the model, the wind speed near the perturbation in EXP2 is slightly greater than that in EXP1. The maximum difference is about 0.6 m s-1. However, the geopotential height near the perturbation in EXP2 is greater than that in EXP1, with a maximum difference greater than 4 gpm. With the intensification of rainfall in EXP2, the latent heat increases. An MβV forms to the northwest of the perturbation at t=19 h. The wind speeds decrease near the vortex center and increase to the northeast and southwest of the vortex. The geopotential heights in EXP2 also decrease quickly near the vortex center, with a difference of less than -8 gpm at t=20 h and less than -14 gpm at t=24 h. The wind speeds to the south and north of the vortex in EXP2 increase with the integration, respectively. The maximum wind speed differences are greater than 2.5 m s-1 at t=24 h and greater than 4 m s-1 at t=30 h. As the scale of the strong wind region is about 100 to 200 km, we follow (Chen et al., 1998) and define it as an mLLJ. This mLLJ can be regarded as the mesoscale disturbance associated with the rainstorms in their study. Our simulation results suggest that the low-level pressure drops are caused by the LHR, and then the mLLJ intensifies and the MβV develops. The southwesterly and northeasterly mLLJs to the south and north of the MβV play important roles in the development of the MβV. Figure5. Differences in 700-hPa wind speed (shaded according to scale; units: m s-1) and geopotential height (lines; contoured every 2 gpm) between EXP2 and EXP1 at (a) t=18 h, (b) t=20 h, (c) t=24 h and (d) t=30 h.
Figure 6 shows the differences in the 500-hPa wind speed and geopotential height between EXP2 and EXP1. The wind speeds above the perturbation in EXP2 increase clearly from t=18 h to t=20 h, while the wind speeds to the southwest of the perturbation decrease. A weak wind speed difference center forms above Wuhan at t=24 h, along with a weak cyclonic circulation. At the same time, the wind speeds to the north and south of the weak vortex in EXP2 increase and decrease, respectively. The evolution of the 500-hPa geopotential height is markedly different from that at 700 hPa. Under the influence of the latent heat in EXP2, the 500-hPa geopotential height rises near the perturbation, leading to a difference greater than 10 gpm by t=20 h. Figure6. As in Fig. 5 but at 500 hPa.
According to hydrostatic balance theory, the thickness of the air column is proportional to average temperature, \begin{equation} \Delta z=\frac{R\overline{T}}{g_0}\ln\frac{p_{\rm b}}{p_{\rm t}} , \ \ (2)\end{equation} where p b and p t are the pressure of isobaric surfaces; ? z and $\overline T$ are the thickness and average temperature of the air column between p b and p t, respectively; R is the gas constant; and g0 is the average of the gravitational acceleration. We define the thickness of the air column and average temperature between 400 and 700 hPa as DH and TA, respectively: \begin{eqnarray} {\rm DH}&=&{\rm {H}}_{400}-{\rm {H}}_{700} ;\ \ (3)\\ {\rm TA}&=&\frac{1}{{\rm log}p_{\rm 6}-{\rm log}p_1}\sum_{n=2}^{5}T_n({\rm log}p_{n+1}-{\rm log}p_{n-1}) . \ \ (4)\end{eqnarray} Here, H is the geopotential height; Tn (n=2...5) is the temperature at the 400, 500, 600 and 700 hPa levels; and pn(n=1...6) is the pressure at 300, 400, 500, 600, 700 and 800 hPa. Next, we define the differences in DH and TA between EXP2 and EXP1 as DDH and DTA: \begin{eqnarray} {\rm DDH}&=&{\rm DH}_{{\rm EXP2}}-{\rm DH}_{{\rm EXP1}} ;\ \ (5)\\ {\rm DTA}&=&{\rm TA}_{{\rm EXP2}}-{\rm TA}_{{\rm EXP1}} . \ \ (6)\end{eqnarray} Figure 7 shows the evolution of average DH and average TA within (30°-31°N, 113.5°-114.5°E). The evolution of average DH is similar to that of average TA, because DH is directly proportional to TA under the hydrostatic balance constraint. In EXP1, average DH decreases from 4504.3 gpm at t=17 h to 4484.2 gpm at t=24 h, and varies little after that; whereas, in EXP2, average DH increases rapidly from about 4495.8 gpm at t=19 h to about 4514.0 gpm at t=21 h. During this period, the upper-level pressure rises and the low-level pressure drops, resulting in the increase in DH. The value of average TA in EXP2 rises by 2.7°C at t=21 h compared with that in EXP1. Figure7. Evolution of (a) average DH (units: gpm) and (b) average TA (units: °C) within (30°-31°N, 113.5°-114.5°E), where the rainfall was relatively intense, in EXP2. The gray line indicates EXP1 and the black line EXP2.
The horizontal cross sections of DDH and DTA are shown in Fig. 8. DDH is negative near the perturbation at t=18 h, which corresponds to a weakly negative DTA center. The DDH intensifies quickly with the integration. The maximum DDH is greater than 21 gpm in association with a positive DTA center with a maximum DTA greater than 2.7°C at t=20 h. The regions of positive DDH and DTA extend, whereas the intensities weaken with the integration before t=24 h. By t=30 h, the maximum values of DDH and DTA decrease to about 18 gpm and 1.6°C, respectively. Figure8. Horizontal cross section of DDH (lines; contoured every 3 gpm) and DTA (shaded according to scale; units: °C) at (a) t=18 h, (b) t=20 h, (c) t=24 h and (d) t=30 h.
Figure 9 shows the horizontal cross section of the 700-hPa difference vorticity and wind between EXP2 and EXP1. Both the difference in vorticity and wind is very small at t=18 h. The difference in vorticity increases in association with the intensification of the difference in wind. The maximum value of the difference in vorticity is greater than 3× 10-5 s-1 near the perturbation at t=20, and it increases rapidly after t=21 h, with the maximum value being greater than 16× 10-5 s-1 at t=24 h. A mesoscale cyclonic circulation forms in the difference in the wind field at t=21 h (not shown). By t=30 h, the difference in vorticity is rather large, with the positive-value region and the cyclonic circulation extending to the meso-α scale. We also see that the regions of positive difference in vorticity all correspond to the cyclonic circulations, suggesting that the latent heating is responsible for the intensification of the MβV. Figure9. Horizontal cross section of the difference in vorticity (shaded according to scale; units: × 10-5 s-1) and wind (arrows; units: m s-1) at 700-hPa between EXP2 and EXP1 at (a) t=18 h, (b) t=20 h, (c) t=24 h and (d) t=30 h. The rectangle denotes the main region of the vortex circulation.
To analyze the energy difference between the two simulations, we defined the difference in the total energy (DTE) per unit mass (Zhang et al., 2003a; Tan et al., 2004) and total kinetic energy per unit mass. The DTE is expressed as \begin{equation} {\rm DTE}=\frac{1}{2}\sum(U'^2_{i,j,k} +V'^2_{i,j,k}+\kappa T'^2_{i,j,k}) , \ \ (7)\end{equation} where U\prime, V\prime and T\prime are the differences in the wind components and the difference in temperature between EXP2 and EXP1; \(\kappa=C_p/R\); i and j run over x and y grid points over the region (20°-50°N, 80°-120°E); and k runs over the 11 levels of 1000, 925, 850, 775, 700, 600, 500, 400, 300, 200 and 100 hPa. The total kinetic energy (TKE) per unit mass is defined as \begin{equation} \label{eq7} {\rm TKE}=\frac{1}{2}\sum(U_{i,j,k}^2+V_{i,j,k}^2+W_{i,j,k}^2) , \ \ (8)\end{equation} where U and V are the x- and y-components of wind speed, respectively, and W is the vertical velocity in the z coordinate. The evolution of the DTE is shown in Fig. 10a. After the MWP is introduced into the model, DTE increases linearly with the integration. However, both TKEs over the same region in EXP1 and EXP2 decrease with the integration (shown in Fig. 10b). This indicates that the decrease of the total energy is compensated by the latent heat energy in EXP2. As some of the latent heat energy is converted into the kinetic energy, the TKE in EXP2 decreases more slowly than that in EXP1. Figure10. (a) Evolution of DTE (units: m2 s-2). Evolution of TKE (units: m2 s-2) in EXP1 and EXP2 over (b) the same region as the DTE computation region, and (c) the Mβ V circulation field as shown in Fig. 9d.
To reveal the development of the MβV, we calculated the 700-hPa TKE over the MβV circulation field (Fig. 10c). The TKE varies little in EXP1, whereas in EXP2 the TKE increases linearly after t=18 h. Consequently, the LHF plays an important role in the evolution of the dynamic fields and the development of the MβV.