-
During 17–23 July 2021, a week-long heavy rainfall event occurred in Henan Province, China. Devastating flooding was caused in highly urbanized areas, leading to huge economic losses of over 100 billion RMB and more than 300 casualties; thus, it was a high-impact weather event that drew wide social concern. In particular, between 0800 and 0900 UTC (1600 and 1700 BJT, Beijing time) on 20 July 2021, an extreme rainfall rate of up to 201.9 mm h–1 (measured by rain gauge) was observed in Zhengzhou (hereafter, 21.7 extreme rainfall), breaking the hourly precipitation record of 198.5 mm h–1 set in August 1975 in China (Ding et al., 1978).
The occurrence of extreme rainfall usually depends on nonlinear interactions between multiscale processes, ranging from large-scale environmental conditions and mesoscale processes that help organize convective storms to storm-scale dynamics and cloud microphysical processes (Ding, 2015; Schumacher, 2017). According to radar observations and the Rapid Update Cycle Version 2 (RUC-2) analysis (Benjamin et al., 2004), Schumacher and Johnson (2005, 2006) found that of all precipitating weather systems, the mesoscale convective system (MCS) was the major producer of extreme rainfall which accounted for about two-thirds of the extreme rainfall events in the United States. Specifically, two common modes of convective organization were revealed for the extreme-rain-producing MCSs. The first type is termed “training line/adjoining stratiform” (TL/AS), which is a linear MCS with cell motion approximately parallel to the convective line. The second common extreme-rain-producing MCS pattern features “back-building/quasi-stationary” (BB) convective cells along a line. The TL/AS MCSs usually form in a very moist, unstable environment on the cool side of a pre-existing slow-moving surface boundary, while the BB MCSs depend more on the mesoscale/storm-scale processes. In both types of MCSs, extreme rainfall most often takes place in the region with repeated passage of convective cells generated upstream, a process commonly known as “echo training” (Doswell III et al., 1996).
There are also many observational and numerical studies of extreme-rain-producing MCSs in China, especially those which occurred in the warm sector of low-pressure systems (e.g., Sun et al., 2019). Previous studies have revealed the importance of a variety of mesoscale/storm-scale processes, such as cold pools, low-level jets, land-sea breezes, mesoscale convective vortices (MCV), boundary layer inertial oscillations, gravity waves, bores, orography, and urbanization, which have greatly advanced the understanding of extreme rainfall generation in China (e.g., Zhang and Zhang, 2012; Luo et al., 2014; Chen et al., 2016; Su and Zhai, 2017; Xu et al., 2017; Zhang et al., 2017; Liu et al., 2018; Zhang et al., 2019, 2020; Abulikemu et al., 2019; Du and Chen, 2019; Xia and Zhang, 2019; Zhang and Meng, 2019; Hua et al., 2020; Wang et al., 2021). For example, Wang et al. (2014) investigated the initiation, maintenance, and properties of convection in an extreme rainfall event observed in the coastal region of South China during the Southern China Monsoon Rainfall Experiment (SCMREX) (Luo et al., 2017). Convective cells were repeatedly triggered on the windward slope of coastal orography in the early stage of the MCS and, in the later stage, along the outflow boundary of the rainfall-induced cold pool. Two scales of training of convective elements were identified, i.e., echo training and band training, which ultimately caused the extreme rainfall. Zhang et al. (2013) examined the extreme rainfall event in Beijing on 21 July 2012. Their results showed that most of the extreme rainfall occurred in the warm sector rather than along the cold front, with the convective cells triggered repeatedly by local orography. Then the back-building convection moved in a manner of echo training along a quasi-stationary linear convective system into Beijing.
Similar to warm-sector heavy rainfall, the 21.7 record-breaking extreme rainfall also occurred in an environment of high equivalent potential temperature (θe) and convective instability (see the environmental conditions in section 3.1), with no pre-existing frontal boundary. However, back-building and subsequent training of deep convection were not found in this extreme rainfall event, different from many previous studies of warm-sector heavy rainfall in China. Instead, the extreme rainfall during 0800–0900 UTC was produced by a single convective storm nearly anchored to Zhengzhou, i.e., a quasi-stationary storm. The main purpose of this work is to reveal the mechanisms responsible for the stationarity of the extreme-rain-producing storm in Zhengzhou. A recent numerical study by Yin et al. (2022) proposed a possible dynamical mechanism that emphasizes the importance of a meso-γ-scale vortex surrounded by arc-shaped updrafts. They stated that the extreme-rain-producing storm was maintained by an arc-shaped convergence in the boundary layer caused by the southwesterly flow of the Huang-Huai cyclone west of Zhengzhou. Herein, we find that a northerly barrier jet (e.g., Overland, 1984; Li and Chen, 1998) causing northerly flows towards Zhengzhou also plays an important role in enhancing the low-level convergence in Zhengzhou. The dynamical balances between low-level flows from different directions and the outflow of the rain-induced cold pool are believed to be primarily responsible for the quasi-stationarity of the storm over Zhengzhou during the record-breaking rainfall hour. While none of the mesoscale processes involved is a new discovery, identifying their roles and contributions constitutes the main contribution of this study, as understanding key processes leading to record-breaking events is clearly important.
The remainder of the paper is structured as follows. In section 2, the observations and reanalysis datasets used in this paper are introduced, describing the numerical experiment setup and the Barnes-type filter (Barnes, 1973) employed for scale separation. The environmental conditions of the 21.7 extreme rainfall event are analyzed in section 3.1, while section 3.2 presents the evolution of the MCS according to both radar observation and numerical simulation. Section 3.3 examines the structure of the low-level mesoscale wind fields surrounding the extreme-rain-producing storm in Zhengzhou, focusing on the formation of a barrier jet along the eastern slopes of the Taihang Mountains (THM). A physical explanation is given in section 3.4 for the rainstorm remaining quasi-stationary in Zhengzhou. A discussion and summary are given in section 4.
-
The evolution of the 21.7 extreme-rain-producing convective storm was examined using the operational Doppler radar observations at the Zhengzhou station. The quality control processes of the radar reflectivity follow that of Huang et al. (2017) and Zhou et al. (2020). Non-meteorological echoes are first identified and eliminated according to a differential propagation phase shift (ΦDP) and specific differential phase (KDP). Then the attenuation of radar reflectivity is corrected based on the linear relationship between polarimetric parameters. The 0.01° × 0.01° gridded hourly precipitation data provided by the China Meteorological Data Service Center was adopted, which was produced via the merger of the Climate Precipitation Center Morphing technique (Joyce et al., 2004) and rain gauge observations (Shen et al., 2014). Moreover, the environmental conditions of the convective storm were studied utilizing a 6-hourly, 0.25° × 0.25° ERA5 reanalysis product (Hersbach et al., 2020).
-
A convection-permitting simulation was conducted using the Advanced Research Weather Research and Forecasting model (ARW-WRF v4.3, Skamarock et al., 2008). Indeed, it is one of the convection-permitting resolution ensemble forecasts that employed different initial and boundary conditions and used different physics parameterizations for the 21.7 extreme rainfall event (Zhu et al., 2022). While still not perfect, the current simulation best reproduced the overall evolution of the convective system and the extreme hourly rainfall under consideration (see section 3.2) among the ensemble members.
Figure 1a shows the coverage of the model’s single domain, which has a horizontal grid spacing of 4 km. There are 51 vertical levels, with the model top placed at 50 hPa. The vertical layer interval varies with height, from about 10 m near the surface, increasing to about 1000 m at the model top. The model’s initial conditions were derived from the ERA5 reanalysis, with the model boundary conditions updated hourly. This real-data simulation ran from 1200 UTC on 19 July 2021 and was integrated forward by 48 hours.
Figure 1. (a) Geographic location of the model domain. The triangle denotes the location of the radar in Zhengzhou. (b) Response functions of large-scale (blue), mesoscale (orange), and small-scale (green) Barnes-type filters.
For the model physics, the NSSL 1-moment 6-class microphysics scheme (Gilmore et al., 2004) was adopted, along with the Mellor–Yamada Nakanishi Niino (MYNN) 2.5-level planetary boundary layer scheme (Olson et al., 2019), the MYNN scheme for the surface layer, the Unified Noah land surface model (Tewari et al., 2004), and the Rapid Radiative Transfer Model for General circulation models (RRTMG) for longwave and shortwave radiation (Iacono et al., 2008). The parameterization of cumulus convection was turned off, given the 4 km horizontal resolution of the model.
-
Given the multiscale processes involved in the development of convective systems and associated severe weather (e.g., heavy rainfall, hail, and convective high winds; Xu et al., 2015a; Chen et al., 2017; Fu et al., 2018; Luo et al., 2018; Zeng et al., 2019), the Barnes filter was employed to decompose the meteorological fields (e.g., horizontal wind and temperature) into large-, small-, and meso-scale components. The former two were straightforwardly extracted from the model outputs by applying a low-pass and high-pass filter, respectively. Then the mesoscale component was obtained by subtracting the large- and small-scale output from the total. More details regarding implementing a Barnes filter can be found in the appendix of Xu et al. (2017). Figure 1b presents the response functions of the aforementioned Barnes filters. In general, the large-scale denotes the synoptic conditions with spatial scales greater than 1500 km; the peak response of mesoscale occurs at about 500 km, while the small-scale primarily consists of horizontal scales less than 200 km, particularly those smaller than 100 km.
-
Figure 2 shows the environmental conditions at 0000 UTC on 20 July 2021. The western Pacific subtropical high (WPSH) was mainly located over the Sea of Japan, with the center of the high near 40°N (Fig. 2a). Compared to its 30-yr climatology, the WPSH was displaced abnormally northward by more than 10° in latitude (Zhang et al., 2021). This northward shift of the WPSH was likely caused by Typhoon In-Fa (2021) over the subtropical ocean, which was approaching southeastern China at this time. Abundant water vapor was transported inland from the ocean by the easterlies between the WPSH and Typhoon In-Fa (2021), which turned to southeasterlies/southerlies along the southwestern periphery of the WPSH (Fig. 2b). Due to the blocking of low-level inflow by the THM to the northwest of Zhengzhou (see the topography in Figs. 2a and 2c), as well as the presence of a mesoscale vortex (Fig. 2c) in the vicinity and to the west of Zhengzhou (which was referred to as Huang-Huai cyclone in Yin et al., 2022), there was significant water vapor flux convergence over Zhengzhou, leading to high total column water (TCW) of more than 70 kg m–2 (equivalent to 70 mm of precipitable water) (Fig. 2b). This stable, large-scale water vapor supply, similar to what was present during the extreme rainfall that also occurred in Henan province in August 1975 (Ding, 2015), provided favorable moisture conditions for subsequent extreme rainfall (Doswell III et al., 1996).
Figure 2. Environmental conditions derived from the ERA5 reanalysis at 0000 UTC on 20 July 2021. (a) Geopotential height (contours, units: dagpm) and horizontal wind (vectors, units: m s−1) at 500 hPa, with the color shading denoting the topography (units: m). (b) Geopotential height (contours, units: dagpm), total column water (shading, units: kg m–2), and water vapor flux (vectors, units: 10−3 kg m s−1) at 700 hPa. (c) Geopotential height (contours, units: dagpm) and equivalent potential temperature (shading, units: K) at 850 hPa. Bold gray lines indicate 200 m terrain height. (d) Skew-T diagram for the inflow atmosphere averaged in an area of 1° × 1° centered at the green cross in (c). Red and blue lines denote temperature and dew-point temperature, respectively. The full (half) wind barb represents 4 m s−1 (2 m s−1).
Benefiting from the transport of warm, moist air from the ocean, Zhengzhou was in a region of high-θe (Fig. 2c), indicating sufficient convective instability. Figure 2d shows an atmospheric sounding in the inflow region (to the southeast of Zhengzhou by ~200 km; see the location of the green cross in Fig. 2c) extracted from the ERA5 reanalysis. The troposphere was nearly saturated, conducive to the production of heavy rainfall by preventing the entrainment of dry air (Schumacher and Johnson, 2005). A relatively dry and stable layer was present below about 900 hPa, which was likely a result of surface radiative cooling during the nighttime and early morning. A moderate convective available potential energy (CAPE) of about 1700 J kg–1 was found at 0000 UTC, which quickly increased to more than 3200 J kg–1 at 0600 UTC (i.e., two hours prior to the extreme rainfall in Zhengzhou) (not shown). The environment inflow veered with height from easterlies to southerlies, with a moderate vertical wind shear of 9.4 m s–1 (10.0 m s–1) between 1000–700 hPa (1000–500 hPa). This condition is consistent with previous studies that extreme rainfall often occurs in an environment of weak-to-moderate vertical wind shear throughout the cloud depth (Maddox et al., 1979).
-
Figures 3a–3d present the evolution of the convective system from radar observations. Also shown are the horizontal winds and the 850 hPa geopotential height from the ERA5 reanalysis. At 0200 UTC (Fig. 3a), there was an MCS in the Henan province, which can be traced back to at least 0600 UTC on 19 July 2021. Widespread convection was found on the eastern/northeastern flank of the low-level mesoscale vortex, which was a region preferable for convection initiation given the thermodynamic destabilization by vortex-induced isentropic lifting (e.g., Raymond and Jiang, 1990; Trier et al., 2000; Trier and Davis, 2002). In particular, a comma-head echo was apparent on the composite radar reflectivity, which wrapped cyclonically around the vortex center. It is essentially the well-recognized MCS-MCV morphology (e.g., Houze, 2004; Trier et al., 2006). As described below, it is just the southern part of the comma-head echo that caused the extreme rainfall in Zhengzhou during 0800–0900 UTC.
Figure 3. Observed radar composite reflectivity (shading, units: dBZ) and geopotential height (blue lines, units: gpm) and horizontal wind (vectors, units: m s−1) at 850 hPa from the ERA5 reanalysis at (a) 0200 UTC, (b) 0400 UTC, (c) 0600 UTC, and (d) 0800 UTC on 20 July 2021. The right column is the simulated composite reflectivity and geopotential height at (e) 0300 UTC, (f) 0530 UTC, (g) 0730 UTC, and (h) 0900UTC. The small-scale component of geopotential height was filtered in (e)−(f) to better reveal the low-level mesoscale vortex. Gray lines indicate 200 m terrain height.
Over time, the low-level mesoscale vortex moved northwestward as it was embedded within the synoptic southeasterly flow, and the MCS moved northwestward accordingly. Two hours later, at 0400 UTC (Fig. 3b), the convection of interest reached the Funiu Mountains (FNM) immediately to the west of Zhengzhou and strengthened, as indicated by the stronger radar reflectivity. The upslope lifting by the FNM played an important role in the organization of the convection (Yin et al., 2022). Furthermore, the insolation at noon remarkably reduced the convective inhibition (CIN) and enhanced the convective instability of the inflow that fed into the convection, as mentioned in section 3.1, which favored the development of intense convection. Notably, a new convective line developed to the west of the convection over the FNM, which was along the eastern/southeastern flank of the low-level mesoscale vortex, a preferable region for the initiation of secondary convection (e.g., Trier and Davis, 2002; Trier et al., 2006; Xu et al., 2017). This convective line intensified over time and became increasingly evident two hours later at 0600 UTC (Fig. 3c). During this time, the convection over the FNM moved further northward, eventually arriving in Zhengzhou at 0800 UTC (Fig. 3d). The convective storm then moved very slowly, producing an extraordinarily heavy rainfall event, exceeding 200 mm between 0800 and 0900 UTC.
Figures 3e–h are similar to Figs. 3a–d but for the evolution of the simulated MCS. In general, the convection-permitting simulation well reproduced the MCS, such as the well-defined comma-head reflectivity that wrapped around the low-level mesoscale vortex, with a long tail trailing southwestward (Fig. 3e). Compared to its observational counterpart, however, the simulated MCS was displaced by a few tens of kilometers to the west and was delayed by about 1 to 2 hours. These timing and positioning biases are common in real-data simulations of MCSs (e.g., Davis and Galarneau, 2009; Meng et al., 2012; Xu et al., 2015b; Abulikemu et al., 2016; Wang et al., 2016). For the simulated MCS, as a whole, it basically moved northwestward in association with the low-level mesoscale vortex. Yet, the southern part of the leading convective line (that eventually caused the extreme rainfall in Zhengzhou) moved northward/northeastward because of the southerly/southwesterly winds of the low-level mesoscale vortex. At 0530 UTC (Fig. 3f), the convection arrived in the FNM and intensified, whereas a secondary convective line began to develop to its west near 112°E. Both were qualitatively consistent with the radar observation. In the next few hours, the simulated convection over the FNM moved further northeastward (Fig. 3g) and reached Zhengzhou at 0900 UTC, i.e., one hour later than reality, with a southward displacement of about 50 km compared to observations (Fig. 3h).
Despite the timing and positioning biases mentioned above, the simulated convective storm also became quasi-stationary in Zhengzhou, producing extremely heavy rainfall of about 160 mm h–1, comparable to the CMORPH precipitation (about 170 mm h–1, Fig. 4). Therefore, it is justified to use this convection-permitting WRF simulation to study the mechanisms responsible for the stagnation of the extreme-rain-producing storm in Zhengzhou.
-
Figure 5a depicts the simulated horizontal winds at 300 m above the mean sea level (MSL) at 0900 UTC. Three airflows (i.e., easterlies, southeasterlies, and northeasterlies) were found to converge in Zhengzhou and the surrounding regions. The three-way converging flow formed a mesoscale convergence zone centered at the southern tip of the convective storm. (The location of the storm was indicated by the 45-dBZ isoline south of Zhengzhou). There were also a few convective-scale features of convergence embedded within this mesoscale region, which were one order greater in magnitude and mainly located on the southern-to-eastern flank of the storm (not shown). Similarly, Yin et al. (2022) also noticed an arc-shaped convergence zone along the south-to-east periphery of the convective storm (see their Fig. 14).
Figure 5. (a) Simulated horizontal wind (vectors, units: m s−1) and the convergence of mesoscale horizontal winds (shading, units: 104 s–1) at 300 m above MSL at 0900 UTC on 20 July 2021. Panels (b) and (c) are similar to (a) but for the large-scale and mesoscale components of the horizontal wind (vectors) and speed (shading, units: m s−1) and pressure (blue contours, units: hPa). Red contours represent the composite reflectivity of 35 dBZ. Colored lines in (a) are the 9-h backward trajectories for air parcels terminating at the core region of the barrier jet along the eastern slope of the Taihang Mountains in (c). Lines AB, CD, and EF in (c) indicate the locations of the vertical cross-section shown in Fig. 6. Red boxes in (c) denote the locations where the mesoscale pressure and northerly wind components are calculated in Fig. 7.
The converging flow pattern was induced by mesoscale winds rather than a large-scale flow, clearly shown after scale separation. The large-scale winds were predominantly southeasterlies in approximate geostrophic balance with the large-scale pressure field (Fig. 5b), indicating little convergence. (The horizontal winds 300 m above MSL were directed slightly toward low pressure due to surface friction.) From the mesoscale fields (Fig. 5c), the convective storm was in a region of obvious horizontal wind convergence, accompanied by negative pressure perturbations in a much broader area than the storm itself. (Hereafter, perturbations are meant as mesoscale components unless otherwise stated.) This meso-β-scale low pressure was caused by a low-level mesoscale vortex in the mountainous areas northwest of Zhengzhou at this time (Fig. 3h). The cyclonic circulation of the low-level vortex resulted in the mesoscale southwesterlies, southerlies, and southeasterlies to the south of the storm. However, to the north of the storm, there were predominantly mesoscale northerly/northeasterly winds, accompanied by positive pressure perturbations, poleward of about 35°N. The high-pressure center adhered to the eastern slope of the THM and gradually decayed to the east (i.e., away from the THM).
Figure 6a shows the wind fields and pressure perturbations in the west-east vertical plane along the high-pressure center. The maximum pressure perturbation occurred on the eastern slope (i.e., windward slope, given the prevailing easterly large-scale flow in Fig. 5b) of the THM, which suggested low-level flow blocking by the high terrain of the THM. According to Miranda and James (1992), mountain flow blocking can be measured by the non-dimensional Froude number, namely,
Figure 6. (a) Simulated in-plane wind fields (vectors, units: m s−1) and mesoscale pressure (shading, units: hPa) in the west-east cross-section along line AB in Fig. 5c at 0900 UTC on 20 July 2021. Vertical velocity (units: m s−1) is exaggerated by 10 times for clarity. Panels (b) and (c) are the same as (a), but for the vertical cross-sections along lines CD and EF in Fig. 5c. White contours in (a) denote the Fr, while they represent the northerly wind component (units: m s–1) in (b) and (c).
where U(z) is the horizontal wind speed at height z, N(z) is the buoyancy frequency, and Hm is the mountain amplitude. The smaller the Fr, the more the airflow is blocked by the orography. As shown in Fig. 6a, the Fr was less than unity below about 1 km above MSL, reaching as low as 0.25 near the surface. Therefore, significant low-level flow blocking took place east of the THM.
Indeed, low-level blocking emerged as early as 0000 UTC (Fig. 7a). In the early stage, the large-scale easterlies were approximately in geostrophic balance, with the southward-directed pressure gradient force (PGF) balanced by the northward-directed Coriolis force. On approaching the THM, the low-level incoming flow decelerated due to the blocking, leading to a reduction of Coriolis force which was thus overwhelmed by the PGF. Consequently, the quasi-geostrophic easterly flow was forced to turn left, producing an ageostrophic northerly wind component parallel to the THM. This pattern can be readily confirmed by the Lagrangian backward trajectories of parcels within the core region of the northerly flow shown in Fig. 5a. These parcels were released at 0900 UTC and traced backward for nine hours. The trajectories were calculated using the model output at 5-min intervals, based on the fourth-order Runge-Kutta scheme adopted in Xu et al. (2017) and Hua et al. (2020).
Figure 7. Time series of the simulated mesoscale (a) pressure (shading, units: hPa) and (b) northerly wind components (shading, units: m s−1) at low levels averaged in the top box and the bottom box in Fig. 5c, respectively.
Given the persistent airflow blocking and the resulting blocking high, the northerly winds accelerated with time (Fig. 7b), forming a barrier jet stronger than 10 m s–1 at 0900 UTC along the eastern slope of the THM (Fig. 6b). The jet core was located close to, but south of, the high-pressure (Figs. 5c and 6c). The barrier jet extended away from the THM by about 72 km (as outlined by the 7 m s–1 contour), which agrees very well with the Rossby radius of deformation, i.e., LR = U/f ≈ 71 km. Here, the large-scale horizontal wind speed and Coriolis parameter were U = 6.2 m s–1 and f = 8.8 × 10–5 s–1, respectively, as estimated at point B in Fig. 5c. This scale is often adopted to measure the width of airmass that could be trapped by a mountain (Overland and Bond, 1995). Moreover, the maximum wind speed of the barrier jet occurred at about 500 m above MSL, roughly midway up the mountain crest, also in good agreement with previous studies (e.g., Yang et al., 2017).
On approaching the storm in Zhengzhou, the northerly barrier jet gradually turned northeasterly, with even greater speed in the canyon region northwest of Zhengzhou (Fig. 5c). The further acceleration of the jet was due to the enhanced PGF, owing to the presence of negative pressure perturbations. This low pressure, as mentioned above, was associated with the low-level mesoscale vortex located to the west of Zhengzhou (Fig. 8a). Therefore, the mesoscale barrier jet was deflected westward upon accelerating.
Figure 8. (a) Simulated temperature (shading, units: K) at 300 m above MSL and mesoscale wind (vectors, units: m s−1) and pressure (contours, units: hPa) at 1500 m MSL at 0900 UTC on 20 July 2021. Red contours are the composite reflectivity of 35 dBZ. (b) Simulated potential temperature (shading, units: K), in-plane winds (vectors), and 45-dBZ reflectivity (black line) in the west-east cross-section along line GH in (a). Vertical velocity (units: m s−1) is exaggerated by 10 times for clarity. The yellow and green boxes denote the locations where the wind speed and easterly inflow are calculated in Fig. 9.
-
Radar observations showed that a single, quasi-stationary storm produced the extreme hourly rainfall in Zhengzhou, which was faithfully reproduced by the real-data, high-resolution simulations conducted herein, as well as in Yin et al. (2022). Notably, convective back-building and echo training, often observed in extreme rainfall events in China (e.g., Zhang et al., 2012; Wang et al., 2014; Liu et al., 2018), were not found in the present case.
Back building of convection depends largely on the mesoscale/storm-scale processes, in particular through the lifting by a cold pool outflow boundary (Schumacher and Johnson, 2005). The cold pool produced by the 21.7 rainstorm was fairly weak, only about 2–3 K colder than the pre-storm inflow to its southeast (Fig. 8a). A west-east cross-section through the cold pool center showed that the cold pool extended upward by no more than 1 km (Fig. 8b). This weak and shallow cold pool was a natural consequence of the very moist environment (Fig. 2d), indicating that the evaporative cooling of rainfall was suppressed substantially (e.g., Davis and Lee, 2012; Xu et al., 2017). Yin et al. (2022) reported no obvious cold pool in this case either. Nevertheless, the convective storm was still well organized and produced extremely heavy rainfall. As shown in Fig. 5, the storm was surrounded by low-level converging flows from different directions. This low-level convergence dynamically lifted the air parcels, promoted intense updrafts, and thus helped maintain the convective storm.
Interestingly, the northeastward-moving convective storm became quasi-stationary once it reached Zhengzhou. The northward movement of the storm was inhibited by the enhanced northerly flow to the north of the storm. The northerly flow originated from the barrier jet on the eastern slope of the THM, and it was accelerated further by the low-level mesoscale vortex west of Zhengzhou (Figs. 5c and 8a). Figure 9 depicts the evolution of the mean southerly wind component averaged between 0–1 km above MSL to the south of the storm and the mean northerly wind component to its north during 0900–1000 UTC. The low-level southerly wind was largely balanced by its northerly counterpart to the north of the storm, with the former about 2.8 m s–1 stronger than the latter. Consequently, the storm only experienced a slight northward movement during this period.
Figure 9. Evolution of the simulated cold pool propagation speed (solid black, units: m s−1) between 0900 UTC and 1000 UTC on 20 July 2021. Also shown are the simulated mean easterly (dashed black), southerly (red), and northerly (blue) wind components in the low levels (0–1 km), averaged in the green, bottom black, and top black boxes in Fig. 8a, respectively.
As for the west-east direction, the mean easterly wind component to the east of the storm was about 10.2 m s–1 (Fig. 9) due primarily to the large-scale flow (Fig. 5). This low-level easterly inflow was nearly completely countered and offset by the westerly outflow of the cold pool. The calculation of the cold-pool outflow speed followed that in previous studies (e.g., Weisman and Rotunno, 2004; Hua et al., 2020):
Here, H is the cold pool depth, set to 1.0 km, according to Fig. 8b. The buoyancy B is given by:
where θ is the potential temperature and qv, qc, and qr are the mixing ratios of water vapor, cloud water, and rainwater, respectively. The base-state variables indicated by overbars are defined as the areal mean in the near-inflow region (see the green box in Fig. 8a) and temporally averaged between 0900 UTC and 1000 UTC. The approximate balance between the low-level easterly inflow and the westerly outflow of the cold pool led to little eastward movement of the storm. These offsetting zonal flows, in conjunction with the analysis in the north-south direction, help to explain why the convective storm became quasi-stationary in Zhengzhou between 0900 UTC and 1000 UTC.
-
On 20 July 2021, heavy rainfall impacted Zhengzhou city in Henan Province, China, with extreme hourly precipitation of up to 201.9 mm occurring between 0800 and 0900 UTC. Based upon radar observations, ERA5 reanalysis, and a real-data WRF simulation with 4-km grid spacing, this work investigated the environmental conditions and evolution of the extreme-rain-producing convective storm, particularly focusing on the mechanisms responsible for the production of the extreme rainfall over the hour.
Results showed that this extreme rainfall event occurred in an environment similar to that of warm-sector heavy rainfall, featuring high temperature and moisture with no pre-existing synoptic frontal boundary. A large amount of warm moist air was transported from the ocean by the easterlies-to-southeasterlies due to Typhoon In-Fa (2021) and the WPSH, the latter of which was abnormally displaced northward. Due to blocking by the Taihang Mountains to the north and a low-level mesoscale vortex west of Zhengzhou, there was notable convergence of low-level momentum and water vapor flux in Zhengzhou and surrounding regions, yielding very favorable conditions for high precipitable water and heavy rainfall. While warm-sector heavy rainfall is often caused by convective back building along the outflow boundary of a rain-induced cold pool and subsequent echo training, interestingly, the 21.7 hourly extreme rainfall was produced by a single, quasi-stationary convective storm. This storm originated from the Funiu Mountains to the southwest of Zhengzhou as part of a convective line that developed on the eastern flank of the low-level mesoscale vortex. Initially, the convective storm moved northeastward due to the steering of the southwesterly flow of the low-level mesoscale vortex but became stagnant in Zhengzhou, where the extreme rainfall occurred.
The overall evolution of the extreme-rain-producing storm was well-captured by a convection-permitting WRF simulation, including the northeastward movement in the early stages and the later stagnation in Zhengzhou. The model results revealed that the northward movement of the convective storm was inhibited by the intensification of a low-level northerly wind component north of Zhengzhou, which was comparable to the southerly wind component associated with the low-level mesoscale vortex. According to scale separation and backward trajectory analyses, this mesoscale northerly flow formed as the larger scale easterly flow impinging on the east-southeast-facing slope of the Taihang Mountain range was blocked, leading to the formation of a southward, terrain-parallel barrier jet. Furthermore, a balance was found between the eastward outflow of the relatively weak precipitation-induced cold pool and the low-level easterly inflow, preventing the storm from moving further east. These low-level flows converged into Zhengzhou from three directions. The combination of low mountains on the west side and an abnormally rich water vapor supply established unusually favorable conditions for the organization and slow movement of sustained deep convection, leading to the record-breaking hourly rainfall over Zhengzhou. These processes are illustrated in a conceptual model given in Fig. 10.
Figure 10. Schematic diagram illustrating the multiscale dynamical processes responsible for maintaining the 21.7 Zhengzhou extreme rainfall. Synoptic-scale flows are denoted by the thick red arrows, whereas thin blue arrows indicate mesoscale flows. The extreme-rain-producing convective storm was fed by the moisture supply from the ocean (denoted by the red arrow labeled as “moisture conveyor belt”) and maintained by the dynamical lifting of low-level converging flows. The northerly flow originated from a mesoscale barrier jet on the eastern slope of the Taihang Mountains, which occurred as the synoptic easterly flow was blocked and formed a blocking high. The southerly flow was associated with a low-level, meso-β-scale vortex west of Zhengzhou. An overall balance was reached between the low-level northerly and southerly winds. Moreover, the easterly inflow prevented the eastward propagation of the rain-induced cold pool. Thus, the convective storm was nearly stagnant in Zhengzhou, resulting in record-breaking hourly precipitation.
Accumulated rainfall depends on both the rate and duration of rainfall. Heavy rainfall tends to occur in locations where high rainfall rates last for a long time (Doswell III et al., 1996). This work provides insights into how multiscale dynamical processes can work synergistically to keep a convective storm quasi-stationary while at the same time providing strong low-level convergent forcing and consequently produce recording-breaking hourly rainfall. Yet, more detailed studies are still needed to fully understand all processes contributing to the 21.7 hourly extreme rainfall. For instance, Yin et al. (2022) emphasized the significance of a meso-γ-scale vortex. As noted in Nielsen and Schumacher (2020), about half of the hourly extreme rainfall events in the United States were associated with low-level rotation, suggesting that the rotationally-induced, upward-directed dynamic perturbation pressure gradient force can promote an intense updraft and thus potentially enhance the short-term rainfall rate (Nielsen and Schumacher, 2018). However, this finding should be taken with caution. Both Trapp and Weisman (2003) and Xu et al. (2015a) noticed that the low-level rotation could also induce a downward-directed perturbation pressure gradient force that would aid in suppressing convection. Thus, the dynamic impacts of low-level rotation on the organization of such an extreme-rain-producing storm remain to be clarified. Moreover, microphysical processes can also affect quantitive precipitation. Using disdrometer and polarimetric radar observations, Chen et al. (2022) revealed unusually high concentrations of raindrops of all sizes (showing combined properties of both maritime and continental convection), implying that both ice-based and warm rain processes contributed importantly to the high rainfall rate in the present case. Detailed microphysical processes, as well as their possible coupling with dynamic processes that support very high rainfall rates, will be investigated in the future.
Acknowledgements. This work is jointly supported by the National Science Foundation of China (Grant No. 42122036), the Second Tibetan Plateau Scientific Expedition and Research (STEP) program (2019QZKK0105), the National Key R&D Programs of China (2018YFC1507300), the National Science Foundation of China (Grant No. 91837207), and the Beijing Climate Center (QHMS2021008).
On the Key Dynamical Processes Supporting the 21.7 Zhengzhou Record-breaking Hourly Rainfall in China
- Manuscript received: 2022-03-08
- Manuscript revised: 2022-05-13
- Manuscript accepted: 2022-05-24
Abstract: An extremely heavy rainfall event occurred in Zhengzhou, China, on 20 July 2021 and produced an hourly rainfall rate of 201.9 mm, which broke the station record for mainland China. Based on radar observations and a convection-permitting simulation using the WRF-ARW model, this paper investigates the multiscale processes, especially those at the mesoscale, that support the extreme observed hourly rainfall. Results show that the extreme rainfall occurred in an environment characteristic of warm-sector heavy rainfall, with abundant warm moist air transported from the ocean by an abnormally northward-displaced western Pacific subtropical high and Typhoon In-Fa (2021). However, rather than through back building and echo training of convective cells often found in warm-sector heavy rainfall events, this extreme hourly rainfall event was caused by a single, quasi-stationary storm in Zhengzhou. Scale separation analysis reveals that the extreme-rain-producing storm was supported and maintained by the dynamic lifting of low-level converging flows from the north, south, and east of the storm. The low-level northerly flow originated from a mesoscale barrier jet on the eastern slope of the Taihang Mountain due to terrain blocking of large-scale easterly flows, which reached an overall balance with the southerly winds in association with a low-level meso-β-scale vortex located to the west of Zhengzhou. The large-scale easterly inflows that fed the deep convection via transport of thermodynamically unstable air into the storm prevented the eastward propagation of the weak, shallow cold pool. As a result, the convective storm was nearly stationary over Zhengzhou, resulting in record-breaking hourly precipitation.
-
Keywords:
- extreme rainfall,
- multiscale processes,
- orography,
- barrier jet,
- low-level mesoscale vortex
摘要: 2021年7月20日,河南郑州发生一次极端降水,小时降水量高达201.9 mm,创造了中国大陆小时降水的新记录。基于天气雷达观测和WRF模式的对流可分辨数值模拟,本文详细研究了郑州极端小时降水的多尺度动力机制。结果表明,本次极端降水的发生环境与暖区暴雨类似,在异常北抬的西太平洋副热带高压和台风“烟花”的共同作用下,大量暖湿空气由海洋输送至郑州。然而,与暖区暴雨中常见的对流单体后向建立和列车效应不同,本次极端降水由单一的准静止风暴所导致,来自郑州北侧、南侧和东侧的低层气流在该地区辐合,通过动力抬升促进了风暴的维持。尺度分离结果表明,低层北风源自于大尺度东风受到太行山阻挡形成的中尺度障碍急流,与郑州西部中尺度低涡导致的南风达到平衡。大尺度东风在输送对流所需的热动力不稳定空气的同时,也有效阻止了较弱浅层冷池的向东移动。因此,对流风暴几乎在郑州停滞,产生了破纪录的极端小时降水。