Characteristic Features of the Evolution of a Meiyu Frontal Rainstorm with Doppler Radar Data Assimilation

During the Meiyu period, floods are prone to occur in the middle and lower reaches of the Yangtze River because of the highly concentrated and heavy rainfall, which caused huge life and economic losses. Based on numerical simulation by assimilating Doppler radar, radiosonde, and surface meteorological observations, the evolution mechanism for the initiation, development, and decaying of a Meiyu frontal rainstorm that occurred from 4th to 5th July 2014 is analyzed. Results show that the numerical experiment can well reproduce the temporal variability of heavy precipitation and successfully simulate accumulative precipitation and its evolution over the key rainstorm area. (e simulated “rainbelt training” is consistent with observed “echo training” on both spatial structure and temporal evolution. (e convective cells in the mesoscale convective belt propagated from southwest to northeast across the key rainstorm area, leading to large accumulative precipitation in this area. (ere existed convective instability in lower levels above the key rainstorm area, while strong ascending motion developed during the rainstorm process. Combined with abundant water vapor supply, the above condition was favorable for the formation and development of heavy rainstorm. (e low-level jet (LLJ) provided sufficient energy for the rainstorm system, and the low-level convergence intensified, which was important for the maintenance of precipitation system and its eventual intensification to rainstorm. At its mature stage, the rainstorm system demonstrated vertically tilted structure with strong ascending motion in the key rainstorm area, which was favorable for the occurrence of rainstorm. In the decaying stage, unstable energy decreased and the rainstorm no longer had sufficient energy to sustain.(e rapid weakening of LLJ resulted in smaller energy supply to the convective system, and the stratification tended to be stable in the middle and lower levels.(e ascendingmotion weakened correspondingly, whichmade it hard for the convective system to maintain.


Introduction
When the East Asian summer monsoon advances northward from mid-June to early July of each year, there often exists a period of continuous rainy weather in the Jianghuai valley between 28 and 34 °N east of Yichang in China.Since this is the plum season, precipitation during this period is called "Meiyu" ("Baiu" in Japan and "Changma" in South Korea).
e Meiyu rain belt extends eastward by thousands of kilometers from the Yangtze River Basin in China to Japan [1][2][3].During the Meiyu period, floods are prone to occur in the middle and lower reaches of the Yangtze River due to the highly concentrated and heavy rainfall.e serious floods occurred in 1954,1991,1998,1999,2003, and 2016 caused huge economic losses [4].Moisture transport by the southwesterly flow along the western rim of the North Pacific Subtropical High (NPSH) feeds the rainband [5,6].
e Meiyu-Baiurainband can be identified by sharp gradients of specific humidity and equivalent potential temperature [7] and a horizontal wind shear line [5].Unlike moisture, the low-level temperature gradient is weak over China and modest around Japan. e Meiyu front is different from the traditional warm/cold front between warm and cold air masses.Despite its importance and socioeconomic impacts, the mechanism and characteristic of the Meiyu-Baiu rainstorm are poorly understood.
Meiyu as a weather phenomenon on meso-synoptic scales has been studied extensively.Many studies have investigated Meiyu precipitation from various perspectives like the physical mechanism for synoptic scale Meiyu frontal rainstorm [8], changes in the Meiyu system under global warming [9], uncertainties in the forecast of mesoscale convective systems that trigger Meiyu rainstorm [10], impacts of urbanization on Meiyu precipitation [11], and the merge of convective clouds during the rainstorm [12].TRMM observations have also been analyzed to explore features of clouds [13] and spatial and temporal characteristics of Meiyu precipitation [14][15][16].Sample and Xie [17] investigated the large-scale features and environmental forcing for the Meiyu rainbelt based on reanalysis products and suggested that the westerly jet stream is the major reason for the formation of Meiyu rainbelt.Sato et al. [18] analyzed the impact of the Kuroshio Current on mesoscale convective systems in the Meiyu front area and found that the intensified evaporation in the Kuroshio Current leads to increases in the moisture content of the air masses in the Meiyu front area, which is favorable for the maintenance of convective instability in the lower troposphere near the Meiyu front.
e results of subkilometer simulation of a torrential-rain-producing mesoscale convective system in Meiyu season show that the model could reproduce the evolution of the dissipating mesoscale convective system (MCS) [19].Based on numerical simulations, Shen et al. [20] studied the relationship between Meiyu rainstorm and mesoscale vortex disturbance originated from the planetary boundary layer (PBL).Nagata and Ogura [21] analyzed the interaction between rainstorm and low-level jet during the Meiyu season in Japan.
e above studies are conducted using reanalysis products, real-time observations, and numerical simulations that can realistically reproduce the observations.e physical understanding of internal convective feedback and external large-scale environment that leads to Meiyu frontal rainstorms may be identified through such methods.
To better understand the characteristics and mechanism of Meiyu frontal rainstorms, there is a need to carry out numerical simulation with high-resolution observed data assimilation.e present study investigates the characteristics and mechanism of a Meiyu frontal rainstorm occurred in the Yangtze River Basin with the local analysis and prediction system (LAPS) as the assimilation system and the weather research and forecast (WRF) as the forecast model.
LAPS, developed and operated by the National Oceanic and Atmospheric Administration's (NOAA) Earth System Research Laboratory in Boulder, Colorado, is a mesoscale meteorological data assimilation tool that employs a suite of observations (meteorological networks, radar, satellite, soundings, and aircraft) to generate a realistic, spatially distributed, time-evolving, three-dimensional representation of atmospheric features and processes [22][23][24][25][26]. LAPS can serve as a tool to initialize local-scale and mesoscale weather forecast models over local to regional domains [27][28][29][30].
e Institute of Heavy Rain (IHR) has been collaborating with NOAA/ESRL (Earth System Research Laboratory) to develop the application method of LAPS in China from 2006 to 2014 [31].e cloud fields have been improved greatly with assimilating Doppler radar data by LAPS [32].Li and Xu [33] developed a correction method to correct radar reflectivity based on hourly surface rain gauge observations, which helps to improve the rainstorm forecast.Gregow et al. [34] investigated the appropriateness of four different methods to produce precipitation accumulation fields using radar data alone or combined with precipitation gauge data.Numerical simulation of a Meiyu rainstorm using this correction method was conducted, and the heavy rainfall simulated is close to observations in terms of the location, intensity, and evolution [35].In the present study using the results of this numerical simulation, the evolution of the Meiyu rainstorm was analyzed to explore the physical mechanism for the genesis, development, and decaying of the rainstorm system.In Section 2, we describe the numerical experiment and compare the simulated precipitation with observations.In Section 3, synoptic background and evolution characteristics of the rainstorm system are analyzed.Section 4 presents the mechanism for the initiation and development of the rainstorm system.Section 5 presents the structure characteristics of the rainstorm system at its mature stage and the mechanism for its decaying.e summary and conclusions are given in Section 6.

Design of Experiments.
Figure 1 shows the model domain and locations of 83 Doppler radar sites, 79 radiosonde sites, and 2049 surface meteorological sites in this study.e observation variables of the operational Doppler radars in China include radar reflectivity, radial velocity, and wide spectrum.e variables of radiosonde observations include pressure, height, temperature, dew point, wind direction, and wind speed at significant levels at each site.e vertical levels include 1000 hPa, 925 hPa, 850 hPa, 700 hPa, 500 hPa, 400 hPa, 300 hPa, 250 hPa, 200 hPa, 150 hPa, and 100 hPa.
e variables of surface meteorological observations include pressure, temperature, dew point, wind direction, and wind speed.Surface rain gauge observations are interpolated to 0.1 × 0.1 °grid points via objective analysis and used for the verification of model simulation results.To ensure the quality of these observations, the single-radar measurements, such as radar reflectivity and radial velocity, have gone through strict quality control including calibration, nonmeteorological return filtering, and velocity dealiasing by the Severe Weather Automatically Nowcast System (SWAN) developed by China Meteorology Administration, which operates on multiple radars to yield real-time observations and products for severe weather nowcasting [36][37][38].
e new correction method [33] was adopted to correct three-dimensional radar reflectivity data based on hourly surface rain gauge observations.LAPS was used to assimilate the corrected radar reflectivity data, radiosonde observations, and surface meteorological observations [35].e

Case Overview and Comparative Analysis.
e observed 24-h accumulative precipitation during 00 UTC 4th to 00 UTC 5th July 2014 is presented in Figure 2(a), which shows a northeast-southwest-oriented heavy rain belt that is more than 700 km long and about 200 km wide.e rainstorm center is located in the central-southern Anhui Province (denoted by the white frame, hereafter key area of the rainstorm), which is about 150 km long and 100 km wide, extending along the east-west direction.
e maximum precipitation is up to 297.9 mm, reaching the extremely heavy rainstorm level.Numerical simulation of 24-h accumulative precipitation is consistent with observations with the maximum precipitation of 301.3 mm, and the simulated rain belt is also oriented northeast-southwest (Figure 2(b)).
Temporal evolution of hourly precipitation in the key area of the rainstorm (Figure 3(a)) shows that the simulated precipitation evolution and magnitude both agree well with the observation.Observations indicate that there existed three heavy precipitation periods, that is, 01-08 UTC, 13-17 UTC, and 18-23 UTC, and the three peak values of precipitation occurred at 06, 15, and 21 UTC, respectively.
e simulated heavy precipitation periods are during 01-11 and 13-23 and the peak values of precipitation are 7 mm•h −1 and 10 mm•h −1 , respectively, which are similar to those from observations.e precipitation intensity within the key area is classi ed into three categories (Figures 3(b)-3(d)) to compute hourly precipitation in each category, which can represent the temporal evolution of precipitation at di erent magnitudes.Here, the curves of Figures 3(b)-3(d) add up to those in Figure 3(a).Figure 3(b) shows that for observed precipitation that is less than 5 mm•h −1 , the temporal evolution of hourly precipitation is gentle with the  Advances in Meteorology magnitude less than 2 mm•h −1 ; however, the simulated hourly precipitation varies largely with a larger magnitude compared with observations.Large difference is found between the simulated and observed evolution of weak precipitation, and the model overestimates the weak precipitation.Figure 3(c) indicates that, for precipitation within the intensity range of 5-15 mm•h −1 , their temporal evolution of hourly precipitation is consistent between observations and simulations, but the model underestimates the precipitation within this intensity range.As shown in Figure 3(d), compared with simulations of precipitation in the previous two categories, the model performs the best for simulation of precipitation that is larger than 15 mm•h −1 , area-averaged hourly precipitation from simulation is consistent with observations in both magnitude and temporal evolution, except that the peak value is slightly biased.is result indicates that the model can well simulate the temporal evolution of heavy precipitation.
Although the numerical experiment successfully reproduces the accumulative precipitation and its temporal evolution over the key area of the rainstorm, there still exist certain deficiencies in the simulation.For example, the period during which the rainstorm maintains is slightly different from that shown by the echoes; the occurrence time of peak precipitation is also different.Despite these deficiencies, the model can well reproduce the three major periods of precipitation.e characteristics of the rainstorm and the mechanism for its evolution will be analyzed in the following sections based on the results of this numerical simulation.

Synoptic Circulation and Evolution
Characteristics of the Rainstorm System  e synoptic background in the lower troposphere demonstrated typical features of largescale background for Meiyu frontal rainstorm.e isoline of θ e 345 K roughly denoted the position of the Meiyu front [16].A quasi-west-east-extending band of dense contours of equivalent potential temperature θ e appeared near 30 °N with the value of θ e within 337-347 K. Signi cant horizontal wind shear, large humidity contrast, and relatively small temperature di erence existed between the northern and southern sides of the Meiyu front.e key area of the rain storm was located within the moist band in front of the Meiyu front, where vertically integrated precipitable water was larger than 60 mm, and a strong southwesterly jet existed in the same area (Figure 4(c)).At 850 hPa, large moisture ux appeared to the southwest of the key area of the rainstorm due to the strong water vapor transport by the southwesterly.e key area was located at the end of the strong southwest-northeast-oriented water vapor transport belt (Figure 4(d)).
e vertical cross section along 116.5 °E in the rainstorm area (Figure 5) indicates that the key area was located between 30 °N and 31 °N (denoted by the red triangle in Figure 5).High equivalent potential temperature controlled this region from surface up to 900 hPa, and the equivalent potential temperature near the surface exceeded 354 K, while a belt of dense, nearly vertically oriented contours of θ e appeared from 900 hPa to 600 hPa with large vertical gradient of θ e on both sides of the belt.θ e decreased with increasing height from surface to 800 hPa along 30 °N to the south of the dense θ e belt, indicating an instable atmospheric strati cation.Along 33 °N to the north of the dense θ e belt, θ e decreased with increasing height from surface to 700 hPa, suggesting that the atmospheric instability was even stronger than that to the south.From ground to the upper troposphere in the key area of the rainstorm, the atmospheric column is moist with the relative humidity exceeding 60% throughout the column, and the relative humidity is higher than 90% from surface to 400 hPa.e relative humidity is higher than 60% to the north of the key area and below 800 hPa to the south.A vertical belt of isothermal advection formed at the interface of the cold and warm air mass over the key area of the rainstorm, with strong cold advection above the region between 29 °N and 31 °N and warm advection to the northern and southern sides of this region.
e above analysis shows that at the beginning of the rainstorm, there existed su cient water vapor above the key area of the rainstorm while the atmospheric strati cation is unstable on both sides of it.Such a pattern of synoptic weather system provided unstable energy, su cient water vapor, Advances in Meteorology and dynamic lifting for subsequent formation of the rainstorm.

Evolution Characteristics of the Rainstorm System.
Radar echoes at 3 km (Figure 6) are implemented to analyze the evolution of the rainstorm system.From the daytime on 3 July to the morning of 4 July 2014, heavy precipitation largely occurred in southeastern Hubei.In central southern Anhui, heavy precipitation mainly occurred from the daytime of 4 July to the nighttime of 5 July 2014.Radar echoes show consecutive genesis of mesoscale convective system that propagated towards the key area from southwest to northeast.At 0300 UTC 4 July (Figure 6(a)), strong echo center was located at southeastern Hubei, while small-scale radar echo center was also found in the key area and maintained there in the subsequent three hours.At 0600 UTC (Figure 6(b)), radar echoes in southeastern Hubei weakened, and the strong echo center with the size of around 20-50 km was located at the key area of the rainstorm.Large areas of weak echoes were found to the north of the key area, suggesting that during 00-06 UTC 4 July, the rainstorm mainly comprised deep convective precipitation in the south and stratocumulus precipitation in the north.e strong echoes kept propagating northeastward in the subsequent three hours, and the mesoscale convective system moved to east of the key area of the rainstorm by 0900 UTC 4 July (Figure 6(c)).Meanwhile, new convective system formed in southwestern Hubei and intensi ed during its moving to the northeast.By 1200 UTC 4 July, the newly generated convective system moved to the west of the key rainstorm area (Figure 6(d)).At the same time, new convective system also formed to the southeast of the key rainstorm area.At 1500 UTC (Figure 6(e)), the convective system to the southeast developed vigorously, while weak echoes continuously formed to the west of the key area.At 1800 UTC (Figure 6(f)), 6 Advances in Meteorology the newly developed echoes to the west of the key area intensi ed and moved eastward.At 2100 UTC (Figure 6(g)), a few convective cells on the scale of around 20-50 km were found in the key area, and the echo intensity was greater than 45 dBZ.Since then, the convective system decayed and eventually disappeared.By 0000 UTC 5 July, only weak echoes were found in the key area of the rainstorm.e convective system in the key area of the rainstorm apparently experienced three stages.In the rst stage, convective cells were largely distributed to the west of the key area and maintained there for about 5-6 hours from 0000 UTC to 0600 UTC 4 July; in the second stage, new weak echoes kept developing and intensi ed during their eastward moving until they moved out of the key area.is stage lasted for about 10 hours from 0700 to 1700 UTC 4 July; in the third stage, strong convective cells continuously formed in central and southeastern portion of the key area.ese strong convective cells maintained for about ve hours from 1700 UTC to 2200 UTC 4 July.e convective cells that formed the mesoscale convective belt propagated from southwest to northeast across the key rainstorm area, resulting in large accumulative precipitation in this area.
is is the main reason for the rainstorm occurrence in this region.is phenomenon is often called echo training [20].
e third stage is the mature stage of the rainstorm, and the maximum precipitation occurred at about 1800 UTC.
Time-distance distribution of hourly precipitation in the key rainstorm area is shown in Figure 7, which indicates that during the entire rainstorm process, there existed three heavy rain belts that moved from southwest to northeast along the heavy rainfall axis EF (shown in Figure 2).e three heavy rain belts corresponded to the meso-β precipitation system with the horizontal scale of 20-250 km that moved from southwest to northeast.Its time scale is 1-8 hours.e rst rainstorm system occurred in the rst ve hours, and the rainstorm center was located within 1.5 °to the west of key rainstorm area with averaged hourly precipitation larger than 16 mm.e maximum precipitation occurred at the 4th hour, and the precipitation intensity was stronger than 20 mm•h −1 .e second rainstorm system occurred during 4-16 hours, when the precipitation was relatively weak but lasted longer.is system covered the entire key rainstorm area during its moving.e third rainstorm system occurred during 14-23 hours, when the precipitation intensity was the strongest and the system moved over a large area and lasted longer than the previous two systems.Precipitation during this period demonstrated the so-called "rainbelt training" phenomenon.e distribution and evolution of hourly precipitation were quite similar to that of radar echo observations (Figure 6), suggesting that the corrected radar re ectivity assimilated by LAPS can provide optimal initial condition for numerical models.e simulated "rainbelt training" is similar to the observed "echo training."Figure 8 shows that in the rst six hours, CAPE was small and decreased with time, while the maximum CAPE occurred at about 500 m in lower level.CAPE decreased from 370 J•kg −1 at 0300 UTC to 300 J•kg −1 at 0600 UTC 4 July, corresponding to weak precipitation in the rst stage (0-6 hours) shown in Figure 3. Since then, the maximum CAPE    Advances in Meteorology remained at about the height of 1 km and its value increased from 300 J•kg −1 to 650 J•kg −1 at 1500 UTC 4 July, and then decreased to 150 J•kg −1 at 0000UTC 5 July.CAPE in the lower levels was small during the entire process.Comparison of CAPE during the heavy rainfall period of 1300 UTC-2300 UTC and that at the maximum rainfall time of 1900 UTC 4 July; it is found that CAPE gradually increased before the heavy rainfall occurred.CAPE was released when the heaviest rainfall occurred and then gradually decreased.e distribution of CIN was consistent with that of CAPE.CIN remained large in lower levels and small above the height of 1 km, suggesting that CAPE was small in lower levels and the atmosphere was stable.Overall, the atmosphere was stable in lower levels with low CAPE but unstable in the upper level with continuously increasing CAPE during the rainstorm process.

Advances in Meteorology
Vertical cross sections of equivalent potential temperature, speci c humidity, and vertical velocity along the rainstorm center (∼31 °N) at various times are displayed in Figure 9, which shows a center of high equivalent potential temperature at 850 hPa over the rainstorm area of 116 °E- 118 °E since 1000 UTC 4 July.
e equivalent potential temperature decreased from lower levels upward to 850 hPa, corresponding to atmospheric instability below 850 hPa in this region.is pattern that is favorable for the maintenance of convective instability in lower levels maintained stably in the rainstorm area except that the center of high equivalent potential temperature slowly moved eastward and was out of the rainstorm area by 2200 UTC. e heaviest precipitation occurred at 1900 UTC (Figure 3(a)).In the prior 12 hours of 1900 UTC, contours of equivalent potential temperature were densely distributed below 900 hPa with large gradient above the western edge of the key rainstorm area along 116 °E, and ascending motion occurred from 900 hPa up to 350 hPa (Figure 9(a)).Convective rainfall started over the key rainstorm area.Five hours later, the strongest updraft moved to 117 °E and vertical velocity increased (Figure 9(b)).
e speci c humidity exceeded 16 g•kg −1 below 850 hPa during the entire rainfall process, suggesting favorable water vapor condition for rainfall.
e above analysis shows that over the rainstorm area, convective instability in lower levels, strong ascending motion during the heavy rainfall period, and abundant water vapor coexisted, leading to favorable condition for heavy rainfall.Convective instability is a kind of potential instability that needs external uplift forcing as a triggering mechanism.It is necessary to study how the low-level unstable energy was triggered and released that subsequently resulted in heavy rainfall during the rainstorm process.

Low-Level Jet.
LLJ is regarded as the most important mechanism that supplies moisture and energy for midlatitude rainstorms.LLJ can rapidly transport warm, moist air mass and lead to strong moisture convergence.Momentum, water vapor, and heat uxes often accumulate in the planetary boundary layer due to the e ect of LLJ before they are transported upwards.
Figure 4(c) indicates that, at the initial time, the southwesterly LLJ at 700 hPa was located to the southwest of the key rainstorm area with the jet axis in the area (29-30 °N, 113-115 °E). is southwesterly LLJ transported water vapor to the key rainstorm area downstream.Many previous studies have shown that the southwesterly LLJ is closely linked to Meiyu frontal rainstorm.Vertical cross sections of water vapor ux and horizontal winds along 29 °N at various times are shown in Figure 10, which is used to analyze changes in the LLJ to the southwest of the rainfall system during the rainstorm process.
e simulated heaviest rainfall in the key rainstorm area occurred at 1900 UTC 4 July (Figure 3(a)).Cross sections along 29 °N (Figure 10) show that LLJ maintained in the upstream to the southwest of the key rainstorm area, and the LLJ slowly moved eastward.Figure 10(a) displays low-level winds and water vapor ux below 3 km eight hours (1100 UTC) prior to the second stage of heavy rainfall.It is found that horizontal winds between 1.6 km and 3 km that exceeded 16 m•s −1 were largely distributed over 113 °E-117 °E, where high water vapor uxes larger than 90 g•kg −1 •m•s −1 also occurred.In the subsequent several hours, the LLJ kept intensifying and moving eastward, while the area of large water vapor ux expanded simultaneously to reach 120 °E and extended downward to the height of 600 m (Figure 10(b)).10

Advances in Meteorology
Temporal evolution of vertically integrated water vapor ux from surface to 850 hPa averaged over the LLJ domain to the south of the key rainstorm area (Figure 11) shows that the low-level water vapor decreased with time in the rst 10 hours from 623 g is is because the rst-stage heavy precipitation was decaying during this period.Since then, the LLJ intensi ed and water vapor ux increased, providing su cient water vapor condition for the development of the second-stage heavy precipitation.e water vapor transport gradually intensi ed during a 10-h period from 1000 UTC to 2000 UTC 4 July and the value increased to 565 g•kg −1 •m•s −1 . is period corresponded to the stage for the rainstorm to develop and mature in the key rainstorm area, suggesting that the LLJ provided enough energy for the rainstorm system.e intensi cation of lowlevel convergence was a critical reason for precipitation to reach its maximum in the key rainstorm area.

Steering Flow in the Mid-level.
e synoptic steering ow is analyzed based on changes of domain-averaged horizontal wind with height in the western half of the key rainstorm area (Figure 12) at various times during the development stage of the rainstorm system.During 1400-1800 UTC 4 July when the rainstorm was developing, easterly winds with the speed less than 10 m•s −1 prevailed in lower levels from 975 hPa to 925 hPa, while the meridional winds were relatively weak with wind speed within 1 to 7 m•s −1 .From 850 hPa to 200 hPa, westerly winds prevailed and the wind speed increased with height from 2 m•s −1 to 20 m•s −1 , and the meridional wind speed is about 10 m•s −1 .e wind speed at 700 hPa remained almost unchanged in ve hours, and the westerly wind speed was around 10 m•s −1 .e above results indicate that above the convective region, the rainstorm system moved northeastward following the southwesterly winds, which were the major steering ow in the midlevel.

Structure Characteristics of the Rainstorm
System at Its Mature Stage and the Mechanism for Its Decaying    12 Advances in Meteorology characteristic changes will be helpful to reveal the physical and dynamical processes in the rainstorm system.Horizontal distributions of vertical velocity w at various heights are presented in Figure 13, which shows that strong ascending motion was located in the key area during the mature stage of the rainstorm and extended to upper levels.Ascending motions in other areas were weak.At the height of 1 km (Figure 13(a)), vertical motion was weak over the entire system with the maximum temporal average being about 0.4 m•s −1 , and the ascending motion was largely found in central key area of the rainstorm.At the height of 3 km (Figure 13(b)), the ascending motion was relatively strong with the maximum temporal average of 1 m•s −1 .At the height of 5 km (Figure 13(c)), the ascending motion significantly intensified compared with that in the lower troposphere.e temporal average of vertical velocity exceeded 2 m•s −1 and strong vertical motions were largely distributed in the central key area.In the upper level of 10 km (Figure 13(d)), the ascending motion was slightly weaker than that in the middle level, and the maximum value was about 1.2 m•s −1 located to the northeast of the key area.Apparently, during the mature stage of the rainstorm system, ascending motion consistently occurred from low level to upper level above the key area with the largest vertical motion in the middle level.e area of ascending motion tilted northeastward with height, indicating a vertically tilted structure of the rainstorm system during its mature stage.e strong ascending motion above the key area was favorable for the formation of heavy precipitation.

Changes in Ambient Atmospheric
Stability.Vertical profiles of CAPE and CIN averaged over the key area during the decaying period of the rainstorm are displayed in Figure 14, which shows that following the decaying of the system, the maximum CAPE at the height of 1.5 km decreased from 350 J•kg −1 at 2000 UTC to 170 J•kg −1 at 2300 UTC. e atmospheric stratification in lower levels became more stable, but CIN below 500m remained large through the process.e stable atmospheric stratification in lower levels indicated that unstable energy decreased during the decaying stage and the rainstorm could not maintain without sufficient energy.

Evolution of the Rainstorm Structure.
Vertical cross sections of radar echo, equivalent potential temperature, and winds at various time along the axis of the heaviest rainfall (EF shown in Figure 2(a)) in the key rainstorm area are shown in Figure 15 for analysis of differences in echo structure, thermal structure, and circulation between the mature and decaying stages of the rainstorm system.Figure 15 shows that mesoscale convective systems on the scale of 20-50 km kept forming in the key rainstorm area, intensified during their northeastward moving, and eventually weakened and moved out of the area.Figure 15(a) shows the situation at 1400 UTC, which was five hours prior to the occurrence of the heaviest precipitation.Area-averaged rainfall intensity at this time was about 4 mm•h −1 and mainly comprised weak precipitation.θ e was lower between 3 and 6 km above the area and decreased with height below 4 km, demonstrating convective unstable structure.e flow field was relatively flat and a new convective cell formed to the southwest of the area.Two hours later, this convective cell moved slowly and intensified.Meanwhile, air mass with high θ e was transported to lower troposphere by southwesterly LLJ, leading to strong ascending motion below 5 km.
e convective cell intensified and the radar echo reached 40 dBZ (Figure 15(b)).At 1700 UTC, this convective cell further intensified and continued to move northeastward.At the same time, high θ e occurred throughout the entire atmospheric column, accompanied with strong ascending motion that could reach the height of 12 km and large echoes that exceeded 45 dBZ.e echo top increased from 7 km to 11 km.Low θ e centers were still found below 4 km along both sides of the convective cell, indicating that there still existed unstable stratification there.Following the development of LLJ and abundant water vapor transport to the key rainstorm area (Figure 15(c)), new convective cells continued to form while the previously formed convective cells intensified and enlarged.At 1800 UTC, the echo top reached 13 km and the ascending motion further intensified.An updraft occurred at lower levels on the southern side of the convective cell and vertically extended to the height of 13 km on the northern side, which effectively transported water vapor upward.Large amounts of water vapor condensed into rainwater during the ascent, preparing favorable conditions for the heaviest precipitation that occurred at 1900 UTC (Figure 15(d)).At 1900 UTC when the heaviest precipitation occurred, three convective cells simultaneously developed and sequentially moved northeastward, forming the so-called "rainbelt training."Meanwhile, continuous energy release during this process broke through unstable stratification in the middle and lower levels, and high θ e became dominant throughout the rainstorm system (Figure 15(e)).At this time, hourly precipitation was 10 mm•h −1 (shown in Figure 3(a)), area-average hourly precipitation for rainfall >15 mm•h −1 was 8 mm•h −1 (Figure 3(d)), indicating that precipitation at this time was largely convective precipitation.Since then, the convective system weakened in the key rainstorm area and moved northeastward.At 2100 UTC, the convective system was broken over the entire area and stratocumulus precipitation prevailed.θ e was lower in the middle and lower troposphere (Figure 15(f)).
Following the occurrence of the heaviest rainfall at 1900 UTC, the LLJ rapidly decreased and the water vapor transport by LLJ became much smaller.Changes in LLJ were highly correlated with the development and maintenance of rainstorm in this region.e rapid weakening of LLJ resulted in lower energy supply to this area, and the stratification in the middle and lower levels became stable, while ascending motion became weaker.As a result, the convective system could not maintain anymore and eventually disappeared.

Summary and Conclusions
Based on numerical simulations with radar data assimilation, the present study analyzes the characteristic features for Advances in Meteorology the initiation, development, and decaying of the Meiyu frontal rainstorm that occurred during 4-5 July 2014.
e results show that simulated temporal evolution of hourly precipitation agrees well with observations.Meanwhile, discrepancy is found between simulated and observed weak precipitation, and the model overestimates weak precipitation but underestimates precipitation of medium intensity.However, it can well reproduce the temporal variability of heavy precipitation and successfully simulate accumulative precipitation and precipitation evolution over the key rainstorm area.e simulated "rainbelt training" is consistent with observed "echo training" on both spatial structure and temporal evolution.Despite slight di erence between the simulated and observed occurrence time of peak precipitation and di erence between the maintenance periods of precipitation and echoes, the model can realistically simulate the three critical stages of precipitation.
e weather background in the lower troposphere demonstrated synoptic weather features typical for Meiyu frontal rainstorm.Signi cant horizontal wind shears, large humidity, and relatively weak temperature contrast could be found between northern and southern sides of the Meiyu front.e key rainstorm area was located in front of the Meiyu front, where the precipitable water was larger than 60 mm.In addition, huge amounts of water vapor were transported to this region by strong southwesterly LLJ.
Radar echoes showed that mesoscale convective systems kept forming and propagating towards the key rainstorm area from southwest to northeast.e convective cells in the mesoscale convective belt propagated from southwest to northeast across the key rainstorm area, leading to large accumulative precipitation and rainstorm in this area.
During this rainstorm process, the vertical atmospheric structure was characterized by low energy and stable strati cation in the lower levels overlapped by continuously intensifying high-energy air mass with unstable strati cation.
ere existed convective instability in lower levels above the key rainstorm area, while strong ascending motion developed during the period of heavy rainfall.Combined with abundant water vapor supply, the above condition was favorable for the formation and development of heavy rainfall.e LLJ provided su cient energy for the rainstorm system and intensi ed the low-level convergence, which was an important reason for the maintenance of precipitation system and its eventual intensi cation to rainstorm.e large-scale steering ow was southwesterly winds that prevailed in the middle and upper levels, which steered the rainfall system to move northeastward.
At its mature stage, the rainstorm system demonstrated vertically tilted structure with strong ascending motion in the key rainstorm area, which was favorable for the occurrence of heavy rainfall.Mesoscale convective systems on the scale of 20-50 km continuously formed and propagated northeastward in the key rainstorm area and eventually weakened and moved out of the area.At the time of the heaviest rainfall, three sequential convective cells formed simultaneously and moved northeastward continuously, leading to the so called "rainbelt training."Continuous energy release broke through the unstable strati cation in the middle and lower levels, and heavy convective precipitation dominated the entire rainstorm process.
In the decaying stage, unstable energy decreased, and the rainstorm no longer had su cient energy to sustain.Convective system weakened and moved northeastward, eventually became a stratocumulus precipitation system, and 14 Advances in Meteorology moved out of the key rainstorm area.e rapid weakening of LLJ resulted in smaller energy supply to the convective system, and the strati cation tended to be stable in the middle and lower levels.e ascending motion weakened correspondingly, which made it hard for the convective system to maintain.As a result, the convective system eventually disappeared.
Our conclusions here are mainly on a Meiyu frontal rainstorm case.Perhaps, the most encouraging aspect of this entire study is that the numerical experiment can well reproduce the temporal evolution of heavy precipitation over the key rainstorm area with corrected radar re ectivity data assimilation so that we can study the characteristic features of the evolution of the Meiyu frontal rainstorm.It has

Advances in Meteorology
further prompted us to vigorously pursue more Meiyu frontal rainstorm cases using real-time quantitative precipitation forecasting similar to that described here.

Figure 1 :
Figure 1: (a) Locations of Doppler radar sites in China (black solid circles), and (b) radiosonde (black hollow squares) and surface meteorological sites (×).

3. 1 .
Synoptic Circulation.Heavy precipitation already started at 0000 UTC 4th July 2014 over the key area in southeastern Hubei and central-southern Anhui, which is located in front of a deep trough.In the upper level, a northeast-southwest-oriented divergence zone occurred to the west of the trough (Figure4(a)), with the divergence larger than 15 × 10 −5 s −1 .In the mid-troposphere,

Figure 2 :
Figure 2: Observed (a) and simulated (b) 24 h accumulative precipitation during 00 UTC 4 -00 UTC 5 July 2014.e white frame indicates the key area of the rainstorm (unit: mm).

Figure 3 :
Figure 3: (a) Time series of hourly average precipitation over the key rainstorm area from observations (black solid line) and simulations (red dashed line) (unit: mm•h −1 ); (b), (c), and (d) are the same as (a) but for the rain rates less than 5 mm•h −1 , greater than 5 mm•h −1 and less than 15 mm•h −1 , and greater than 15 mm•h −1 , respectively.X axis indicates the forecasting time from 0100 UTC 4 July to 0000 UTC 5 July 2014.

Figure 4 :
Figure 4: Large-scale background at 0000 UTC 4 July 2014: (a) winds (wind barbs, each long barb represents 4 m•s −1 , the same hereafter), geopotential height (black solid contours, unit: 10 dgpm), temperature (red dashed lines, unit: °C), and horizontal divergence (shaded in grey, unit: 10 −5 •s −1 ) at 200 hPa, and the yellow frame denotes the key area of the rainstorm; (b) same as (a) but for 500 hPa, and the areas shaded in grey show temperature advection (unit: 10 −5 •K•s −1 ); (c) equivalent potential temperature black solid lines, unit: K; the red lines indicate the isolines of 345 K and vertically integrated precipitable water (unit: mm) at 700 hPa; (d) same as (a) but for 850 hPa, and the shaded areas indicate water vapor ux (unit: m•g•s −1 •kg −1 ).

4 .
Mechanism for the Initiation and Development of the Rainstorm System4.1.Atmospheric Stability.Domain-averaged convective available potential energy (CAPE) and convective inhibition (CIN) in the rainstorm area are computed at each level (Figure 8) to analyze the distribution of atmospheric energy.

Figure 5 :
Figure5: Vertical cross section along 116.5 °E of equivalent potential temperature black solid lines, unit: K the blue line is the isoline of 345•K, horizontal winds (wind barbs), temperature advection (red dashed lines: ≥0 °C, green dashed lines: <0 °C), relative humidity (shaded, unit: %). e red triangle denotes the key area of the rainstorm and the vertical coordinate is pressure (unit: hPa).

Figure 7 :Figure 8 :
Figure7: Temporal evolution of hourly precipitation (unit: mm•h −1 ) averaged over all grids along AB (shown in Figure2(a)) in the key rainstorm area (the white frame shown in Figure2(a)).X-axis indicates the distance (unit: 10 km) to the western boundary (i.e., AB) of the key rainstorm area along the BD direction (shown in Figure2(a)).Y-axis shows the forecasting time (unit: hour) from 0000 UTC 4 -0000 UTC 5 July.

Figure 15 :
Figure 15: Vertical cross sections along EF (shown in Figure 2(a)) of radar echo (shadings, unit: dBZ), equivalent potential temperature (solid lines, unit: K blue lines for ≤350 K and black lines for >350 K), and winds at various times (a) 1400 UTC, (b) 1600 UTC, (c) 1700 UTC, (d) 1800 UTC, (e) 1900 UTC, and (f ) 2100 UTC.Y-axis indicates the height (unit: km).X-axis shows the distance along the direction shown by BD in Figure 2(a) to the western boundary (AB) of the key rainstorm area (unit: km).
5 °E of equivalent potential temperature black solid lines, unit: K the blue line is the isoline of 345•K, horizontal winds (wind barbs), temperature advection (red dashed lines: ≥0 °C, green dashed lines: <0 °C), relative humidity (shaded, unit: %). e red triangle denotes the key area of the rainstorm and the vertical coordinate is pressure (unit: hPa).