Assimilation of Doppler Radar Data and Its Impact on Prediction of a Heavy Meiyu Frontal Rainfall Event

Operational Doppler radar observations have potential advantages over other above-surface observations when it comes to assimilation for mesoscale model simulations with high spatial and temporal resolution. To improve the forecast of a heavy frontal rainfall event that occurred in the Yangtze-Huaihe River Basin from 4 July to 5 July 2014 in China, operational radar observations are assimilated by the Local Analysis and Prediction System (LAPS). Radar reflectivity data are used primarily in the LAPS cloud analysis procedure, which retrieves the number of hydrometeors and adjusts the moisture and cloud fields. Radial velocity data are analyzed through the LAPS wind analysis-based successive correction method. A new correction method is developed to correct threedimensional radar reflectivity data based on hourly surface rain gauge observations. 'e performance of the correction method is demonstrated by assimilating radar reflectivity observations into LAPS. Experiments with different radar data assimilation are examined. Results show that the assimilation of radar data can effectively correct the background errors and improve the heavy rainfall forecast. 'e simulated intensity, pattern, and temporal evolution of the heavy rainfall event are better improved with radar reflectivity assimilation, especially when the correction method is implemented to correct radar observations.


Introduction
From mid-June to early July each year, accompanied with the northward advance of the East Asian summer monsoon, there always exists a period of seasonal rainfall maximum in the Yangtze-Huaihe River Basin between 28 and 34 °N to the east of Yichang in China.Rainfall during this period is called "Meiyu" in China and "Baiu" in Japan, since this period corresponds to the maturing season of plums.Flood disasters in the Yangtze-Huaihe River Basin triggered by rainstorms along the Meiyu front in the middle and lower reaches of the Yangtze River are among the most important meteorological disasters in China.Due to the special characteristics of synoptic background, mesoscale structure, multiscale system interaction, and triggering and maintenance mechanisms for the Meiyu frontal rainstorms, monitoring and forecasting the Meiyu frontal rainstorms is a difficult issue that attracts great attention of meteorologists in China.
Quantitative precipitation forecast (QPF) associated with Meiyu frontal rainstorm remains a key challenge for numerical models.Many studies have been conducted to explore various aspects of the Meiyu frontal storm, such as synoptic scale physical mechanisms for the storm [1], changes of the Meiyu system in the context of global warming [2], uncertainties and physical mechanisms for mesoscale convective systems to trigger heavy Meiyu frontal rainfall [3], impacts of urban expansion on precipitation during the Meiyu period [4], convergence of convective clouds during the heavy Meiyu frontal rainfall process [5], and simulations of Meiyu frontal rainstorms in the middle and lower reaches of the Yangtze River associated with mesoscale vortex [6].
With more and more Doppler radar deployment, realtime observational data become increasingly rich.In addition, Doppler radar can provide huge amounts of detailed information with high spatial and temporal resolution including three-dimensional clouds, wind and precipitation, which makes it possible to further improve the Meiyu frontal rainstorm forecast.Assimilating radar data into highresolution numerical models to improve short-term quantitative precipitation forecast has been a hot research topic worldwide in recent years [7][8][9][10][11][12].
Variables obtained from direct radar observations include radial velocity and reflectivity.Radar reflectivity reflects a variety of hydrometeors, which usually have a relatively short forecast valid period.Methods of assimilating radar reflectivity include direct, indirect, and transitional assimilation.With the indirect assimilation method, variables that can be adopted by numerical models are retrieved from radar observations first.For example, precipitation is first estimated by radar reflectivity, and condensational heating retrieved by precipitation is then forced to be adopted in numerical models using the Newton relaxation approximation method [13].Direct assimilation means that radar reflectivity is directly assimilated into numerical models as a model variable using variational assimilation method [14] or ensemble Kalman filter method [15][16][17].Transitional assimilation includes two types of assimilation.One is to accomplish the radar reflectivity assimilation in two steps; that is, a simple one-dimensional variational method is used first to retrieve a model forecast variable from radar observations, and a three-dimensional or four-dimensional variational assimilation system is then utilized to assimilate the retrieved variable [18,19].e other is to assimilate radar reflectivity by a cloud analysis procedure, such as the cloud analysis system in the Local Analysis and Prediction System (LAPS), to provide the initial condition that is not only dynamically and thermally consistent but also contains cloud microphysical variables to improve short-term quantitative precipitation forecast [11,20,21].
No matter what method is used for assimilation of radar reflectivity, the prerequisite is that radar reflectivity observational data must be accurate enough and the effective detection range should cover the model forecast domain.
e operational Doppler weather radars located in China have different wavebands in the east and west.e western and eastern regions are dominated by the C-band and S-band radars, whose wavelengths are 5 cm and 10 cm, respectively.Although radars produced by different manufacturers have the same band, there are differences in their spatial resolutions and data format.
is may impair the quality of radar observations and thus affect the assimilation in numerical models.Large amounts of rainfall observations collected at surface stations located within the region of China Mesoscale Observation Network can well compensate for possible missing information in the radar blind zone.e relationship between precipitation intensity I and radar reflectivity Z is the basis for quantitative radar measurement of precipitation.A new correction method based on the relationship between Z and I is developed to correct threedimensional radar reflectivity data using hourly surface rain gauge observations.Assimilation experiments on a mudslide rainstorm show that this method can improve the capability of radar reflectivity assimilation and thus yield better heavy rainfall forecasts [22].

LAPS, developed and operated by Earth System Research
Laboratory of National Oceanic and Atmospheric Administration (NOAA) in Boulder, Colorado, combines a wide array of observed meteorological datasets into a unified atmospheric analysis with a time interval of an hour or less.It is also 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, threedimensional representation of atmospheric features and processes [20,[23][24][25][26]. LAPS can serve as a tool to initialize local-scale mesoscale weather forecast models over local to regional domains [27][28][29].Because it provides cloud ice, cloud water, rain, snow, and graupel analysis to initialize forecast models, LAPS can be used to provide a hot start for operational forecast models [20].
Since 2006, the Institute of Heavy Rain (IHR) has been collaborating with NOAA/ESRL to develop the application method of LAPS in China [30].e wind and humidity fields have been improved greatly with the input of wind profile radars and GPS observations into LAPS [31][32][33].Even though when sometimes the radar observation of China was assimilated by LAPS, there were no satisfactory effects because of the difference in the resolution of radar observations at different wavelengths.In this study, we attempt to adopt a new correction method to use the hourly surface rain gauge observations to correct the radar reflectivity factors, and study the impact of China radar observation assimilated by LAPS to a Meiyu frontal rainstorm forecast with WRF.It may improve the local rainstorm forecast if this method is applied in operation.We use LAPS to assimilate observations to a hot start of WRF.Numerical experiments are performed to simulate a heavy Meiyu frontal rainfall case and explore the sensitivity of radar data assimilation.e impact of the radar data assimilation with the correction method on the prediction of a heavy rainfall event is demonstrated.In Section 2, we describe the data and the correction method, the model and its configuration, the experiments used to assimilate radar observations, and the heavy rainfall case selected for the present study.Section 3 presents the results, and the summary and conclusions are given in Section 4.  e second letter means the different radar manufacturers in China.Both S-band and C-band operational Doppler weather radars are operating in two volume coverage patterns (VCPs) to sample precipitating systems.VCP11 provides 14 unique elevation scans every 5 minutes, while VCP21 provides 9 unique elevation scans every 6 minutes, for elevation angles of 0.5 through 19.5 °. e observation variables of all the operational Doppler radars are same including radar re ectivity, radial velocity, and wide spectrum, but there are di erences in spatial resolution and data format.

Data and Methodology
To ensure the quality of these observations, the singleradar measurements, such as radar re ectivity and radial velocity, have gone through strict quality control including calibration, nonmeteorological return ltering, and velocity dealiasing by the Severe Weather Automatically Nowcast System (SWAN) developed by China Meteorological Administration, which operates on multiple radars to yield real-time observations and products for severe weather nowcasting [34][35][36].en, the operational Doppler radar observation data of di erent types are uni ed to the same data format required by LAPS.Finally, a new correction method is used to correct the radar re ectivity factors.
In addition to Doppler radar data, other observations used in this study include data from 79 radiosonde sounding sites and 2049 surface meteorological sites (Figure 1(b)).e variables of radiosonde observations include pressure, height, temperature, dew point, wind direction, and wind speed at signi cant levels at each site.
e variables of surface meteorological observations include pressure, temperature, dew point, wind direction, and wind speed.Surface rain gauge observations are used for the veri cation of model simulation results.[22] introduced in detail a correction method to correct radar re ectivity using hourly surface rain gauge data.Figure 2 shows the diagram of the method with LAPS in this study.Observations of radar re ectivity are interpolated into LAPS three-dimensional (x, y, z) grids to obtain the radar re ectivity Z LAPS (x, y) at grids within the LAPS domain.Optimum coe cients in the dynamic relationship of Z − I are identi ed using the hourly surface rain gauge observations and radar re ectivity at each level.Suppose that each hourly rainfall I G (I, J) at the surface station (I, J) corresponds to the radar re ectivity dBZ, then dBZ can be used to derive rainfall intensity based on the relationship between radar re ectivity and rainfall intensity, that is, Z aI b .Hourly precipitation I R (I, J) estimated from radar observations can be obtained by time integration of the estimated rainfall intensity.To obtain the optimum coe cients a and b, the discriminant function CTF |I G − I R | is de ned.Initial values of a and b are speci ed rst, and these values then change within a certain range.When CTF reaches its minimum value, the corresponding values of a and b are the optimum values we seek.Based on the optimum coe cients of a and b at each level, the Z − I relationship, and hourly rainfall observations I G (I, J), we can retrieve the radar re ectivity Z G (I, J) at each level above the corresponding surface station.Based on Z G (I, J), Z LAPS (x, y) is corrected by using the correction method, and the three-dimensional radar re ectivity Z COR (x, y) can be obtained by iterative method.Finally, we use LAPS to assimilate Z COR (x, y) to perform the numerical experiments.

Case Overview and Design of Experiments.
A ected by the slow eastward movement of the low-level southwest vortex and the shear line over the Meiyu front, the rst heavy Meiyu frontal rainfall in the Yangtze-Huaihe River Basin ranges from 4 July to 5 July 2014.e observed 24 h accumulated precipitation is shown in Figure 3(a).e rain belt was oriented from southwest to northeast.e heavy rainfall area was mainly located in Hubei, Anhui, and Jiangsu provinces (the red box shown in Figure 3).Several obvious rainfall centers with 24 h accumulated precipitation amount greater than 100 mm were found along the rain belt.
e 24 h maximum of 297 mm was located at Yuexi in Anhui Province.e heavy rainfall event lasted for a long period and covered a large area, causing a severe ooding along the Yangtze River and tremendous economic losses.
e mesoscale nonhydrostatic Weather Research and Forecast (WRFV3.4.1) model is used in this study.LAPS is utilized to assimilate observations.We evaluate the performance of radar data assimilation mainly by examining the accuracy of the subsequent forecast.e primary goal of this work is to test the impact of assimilation of operational Doppler radar observations into initial conditions on the forecast of heavy rainfall events.ree experiments are conducted for the case study.In the control experiment (CON), radiosonde and surface meteorological observations are assimilated.Besides radiosonde and surface observations, operational Doppler radar observations are assimilated in the other two experiments (RAD and RADC).e di erence between these two experiments is whether the correction method is used to correct radar re ectivity.e method is implemented in RADC but not in RAD. e horizontal grid resolution of the WRF model is 9 km.ere are 45 levels in the vertical and the model top is 50 hPa.e model domain covers an area of 2000 km × 2000 km. e initial and lateral boundary elds of all experiments are extracted from the National Centers for Environmental Prediction (NCEP) Global Data Assimilation System (GDAS) Final Operational Global Analyses (1 °× 1 °).Important physics schemes used in the present study include the WSM6 microphysics scheme, the YSU planetary boundary layer scheme, the Kain-Fritsch cumulus parameterization scheme, the Rapid Radiative Transfer Model longwave radiation scheme, and the Dudhia shortwave radiation scheme.e model is integrated for 24 hours from 00 UTC 4 July to 00 UTC 5 July 2014.

Simulation of Precipitation.
ere existed a northeastsouthwest-oriented rainfall belt (black solid box) in this heavy rainfall event.e key rainstorm area with 24 h accumulated precipitation amount exceeding 100 mm was mainly located in eastern Hubei, and central and southern Anhui (denoted by the red solid box shown in Figure 3(a)).Comparison of the 24 h accumulated precipitation simulated by the three experiments with the rainfall observations shows that all the three experiments can realistically reproduce the location of the main rainfall belt, and the strongest precipitation area is located in the lower reaches of the Yangtze River.However, there are obvious di erences in rainfall intensity between the simulations and observations. .RAD can well simulate the heavy rainfall center located in central Anhui, but it creates larger area receiving rainfall than observations (Figure 3(c)).e rainstorm simulated in RADC experiment is similar to observations in terms of location and intensity, and RADC can successfully correct the overestimation of precipitation in Jiangsu found in the other two experiments (Figure 3(d)).
Regarding the temporal evolution of observed 6 h accumulated rainfall from 00 UTC 4 July to 00 UTC 5 July 2014, the rainfall amount of the rst 6 hours over the key rainstorm area located in eastern Hubei and Anhui is pretty large, and the largest rainfall amount at the rainstorm center is above 50 mm (Figure 4(a)).e accumulated rainfall in the second 6-hour interval is weak, and the rainfall area extends eastward (Figure 4(b)).e 6 h accumulated rainfall amount begins to increase again from 12 UTC to 18 UTC (Figure 4(c)).
e strongest precipitation over the key rainstorm area occurs in the last 6 hours (Figure 4(d)).Comparing the temporal evolution of 6 h accumulated precipitation simulated in the three experiments with observations, it is clear that the RADC simulation is the closest to observations.e rst 6 h rainfall amount over the key rainstorm area simulated in RADC is larger than that in the other two experiments, and the rainfall amount of the precipitation center is greater than 50 mm.e precipitation weakens six hours later and starts to increase with time 12 hours later (Figures 4(m)-4(p)).In CON, the 6 h rainfall amount keeps increasing with time (Figures 4(e)-4(h)).
e rainstorm center is located to the east of the observed center.e rainfall belt moves faster, and the rainfall amount is larger than observations.e evolution of 6 h precipitation simulated in RAD is better than that in CON.RAD reproduces the strong-weak-strong precipitation variation, which is similar to observations.However, the location of the rainstorm center simulated in RAD is located to the northeast of the observed center, which is worse than the simulation in RADC (Figures 4(i)-4(l)).
To compare the temporal evolutions of the rainfall belt between observations and simulations of the three experiments, we calculate hourly rainfall averaged over the area indicated by black solid box along line EF (Figure 3(a)) direction.Time series of hourly average rainfall along line FH direction (shown in Figure 3(a)) are presented in Figure 5.

Advances in Meteorology
Observations (Figure 5(a)) show clearly that the disperse heavy rainfall occurred in the southwestern area, and the rainfall in the northeastern area (red box) corresponds to temporal evolution of the heavy rainfall averaged over the key rainstorm area (red box in Figure 3(a)).Rainfall over this region is mainly concentrated over two periods: 0-8 hours and 14-24 hours.CON fails to reproduce the two precipitation periods, and the rainfall in the later period is heavier than observations (Figure 5(b)), while RAD and RADC with radar data assimilation can reproduce the two rainfall periods similar to observations.e simulated hourly precipitation amount over the later period in RAD is too high, and the evolution of rainfall simulated in RADC is closer to observation than that in RAD (Figures 5(c) and 5(d)).
Time series of hourly average precipitation over the key rainstorm area from observations and three experiments are displayed in Figure 6, which shows that all the three experiments can well simulate the two rainfall periods in this key area.
Here, OBS refers to the hourly surface rain gauge observations.e rainfall intensity in the rst period simulated in RAD and RADC with radar data assimilation is larger and that in CON is smaller than observation, while the rainfall intensity in the second period is better simulated in RAD and RADC than in CON.In particular, the evolution of hourly rainfall simulated in RADC is consistent with observations (Figure 6(a)).From the time series of hourly rain rate with various intensities, it can be found that the evolution of light rainfall (i.e., rain rate is less than 5 mm/hr) in the three experiments is similar to observations (Figure 6(b)).Simulations of moderate and heavy rainfall in RAD and RADC are better than that in CON, and the heavy rainfall simulation is the best in RADC, which yields comparable results with observations (Figures 6(c) and 6(d)).
e above results indicate that radar data assimilation can well improve the simulation of rainfall, while the assimilation of radar re ectivity corrected by the correction method can better improve the simulation of heavy rainfall.6

Advances in Meteorology
In summary, the simulated characteristics of precipitation distribution of this heavy rainfall event indicate that the three numerical experiments can realistically reproduce the southwestern-northeast precipitation belt.e simulations are close to observations, and the two rainfall periods over the key rainstorm area are well described, too.However, there exist certain di erences between the simulated and observed location and intensity of heavy rainfall center.It is found that the initial background eld and boundary conditions used in the experiments can re ect the large-scale background circulation during the heavy rainfall event, and the assimilation of observational data can adjust the initial condition of the model, which improves the simulation of the evolutional characteristics of the rainstorm system.e characteristics of rainfall distribution and evolution with time simulated in RAD and RADC with assimilation of radar data are better than that simulated in CON.Since the correction method is used to correct radar re ectivity based on hourly surface rain gauge observations in RADC, the location and intensity of heavy rainfall simulated in RADC are consistent with rainfall observations and better than that simulated in RAD.ese results show that the assimilation of radar data can e ectively correct the location errors in the background and improve the simulation of the heavy rainfall in the numerical model.e heavy rainfall simulated in RADC is the closest to observations in terms of the location, intensity, and evolution of rainfall with time.e overestimation of rainfall in the other two experiments is successfully corrected in RADC.

Analysis of Initial Conditions.
In order to explore the reasons for the improvement of the simulation of this rainfall event by radar data assimilation, the initial conditions of the three experiments are compared and analyzed.Figure 7 presents the initial elds of the three numerical experiments at 00 UTC 4 July 2014.e initial condition of wind and relative humidity elds at 700 hPa and geopotential height at 500 hPa for CON (Figure 7(a)) shows that the precipitation area (red box) is located in front of the upper trough, where relative humidity is high.Rainfall is a ected by an eastward-moving low vortex and low-level shears, and the rainfall area is underneath the region dominated by the southwesterly low-level jet stream (LLJ).Compared with the initial condition for CON (Figure 7(a)), the southwesterly LLJ is stronger in the initial conditions for RAD and RADC (Figures 7(b) and 7(c)).Radar data assimilation promotes the transport of water vapor to the heavy rainfall area, which is favorable for the maintenance and development of the rainstorm system.Comparing the wind elds in initial conditions for RAD and CON shown in Figure 7(d), it is obvious that the wind speed in the precipitation area (red box) is increased by 6 m/s and the southwesterly wind component is also increased by radar data assimilation.Di erence of speci c humidity in initial conditions between RADC and RAD shows that the speci c humidity is slightly increased by 1 g/kg with the correction of radar re ectivity.e above results suggest that the southwesterly low-level jet stream in the rainstorm area is enhanced by assimilating radial wind data, and the correction of radar re ectivity is helpful to increase humidity over the rainstorm area.e initial background elds are adjusted to be more realistic due to radar data assimilation, providing favorable environmental conditions for the development and maintenance of the rainstorm system.
e LAPS cloud analysis system has been updating microphysical variables decades [20].It also adjusts the humidity eld based on observational data.In this study, the increments of cloud microphysical variables are analyzed to evaluate the impact of radar data assimilation.Figure 8 shows cross sections of radar echo, humidity, and increments of cloud hydrometeors along line MN (shown in Figure 3(a)) at the initial time due to radar data assimilation.At 0000 UTC 4 July 2014, several strong radar echoes are distributed along the cross section shown in Figure 8. e strongest radar echo is located in the key rainstorm area (the red box in Figure 3(a)).From Figure 8(a), it can be seen that positive di erences in radar echo between initial conditions for RADC and RAD are mainly located in the key rainstorm area, and the di erence extends upward from surface to the middle troposphere (Figure 8(a)). is result indicates that radar echoes over the rainstorm area are enhanced by using Figure 6: Time series of hourly average precipitation over the key rainstorm area from observations (solid line, black) and simulations of three experiments (dashed lines) (a) (unit: mm/hr).(b), (c), and (d) are the same as (a) but for the rain rates less than 5 mm/hr, greater than 5 mm/hr and less than 15 mm/hr, and greater than 15 mm/hr, respectively.e horizontal axis is forecast hours from 0000 UTC 4 July to 0000 UTC 5 July 2014.
8 Advances in Meteorology the correction method to correct radar re ectivity.Since there are no hydrometers in the conventional observations and radar re ectivity is used to correct the cloud analysis, the analyzed hydrometeor increments are mainly concentrated in the area where radar echo increment occurs.e humidity analysis method developed by Birkenheuer [24] indirectly adjusts the humidity eld using radar re ectivity.is incloud humidity correction is shown in Figure 8(b) as positive humidity increases.It is mainly distributed below 850 hPa, and the maximum increment is more than 3.5 g/kg.e positive humidity disturbance may be bene cial for counteracting the negative buoyancy associated with the introduction of hydrometeors and evaporative cooling and for maintaining the development of rainstorm system.e positive increments of cloud water in the upper troposphere (Figure 8 e increment of the hydrometers in the lower troposphere is much larger than that in the upper layer.e increment of snow is the largest with a maximum value of more than 1.4 g/kg, followed by that of rainwater with a maximum value of about 0.7 g/kg.e increments of cloud water and cloud ice are the minimum, and their maximum increment is less than 0.1 g/kg. is may be due to the fact that the strong radar re ectivity at the middle and lower levels is greatly corrected by using the correction method, which subsequently a ects the analysis of low-level hydrometers in LAPS cloud analysis.

Simulation of Physical
Variables.In order to analyze in detail the forecast di erence between the three experiments, observed and simulated radar echoes at di erent time and levels are plotted.Figure 9 shows radar echoes at the height of 3 km at 0300 UTC 4 July 2014.e horizontal resolution of the simulation experiments is 9 km. is resolution probably is too coarse, and thus, small convection cells observed by radars are not simulated in the three experiments.However, the main convection cells can be reproduced.
ere are two main convective cells observed by radars.e one occurred in southeastern Hubei is named A, and the other occurred in southern central Anhui is named B (shown in Figure 9(a)).Compared with observations, the  e shape and position of the convective cell B simulated in RADC are better than that in RAD, and its position in RAD is far too west (Figures 9(c) and 9(d)).Similar situations can be found at other time.erefore, we can conclude that the initial model conditions are improved with radar data assimilation, and the development and evolution of convective cells in this heavy rainfall event are better simulated with radar data assimilation.RADC with the correction method yields the best result.
To quantitatively evaluate the precipitation forecast, 1 h precipitation forecast threat score (TS) [37] is calculated.e threat score (TS) is simply the ratio of successful event forecasts to the total number of event forecasts.Consider a set of forecasts that can have only two alternatives (e.g., yes and no).TS is de ned as where X denotes the number of positive (yes) forecasts that correspond to an occurrence of the event, Y denotes the number of events that occurred in conjunction with a negative forecast, and Z denotes the number of negative forecasts that did not have any associated events.Figure 10 shows 1 h TS for 0.1 mm, 5 mm, 10 mm, and 15 mm thresholds in the rainfall area (black box in Figure 3(a)).For the 1 h precipitation exceeding 0.1 mm, the TSs of CON, RAD, and RADC are similar (Figure 10(a)), which suggests that the rainfall location and evolution with time simulated by the three experiments are similar.For the heavier rainfall (Figures 10(b)-10(d)), the 1 h TSs of RAD and RADC with assimilation of radar data are greater than that of CON; especially in the rst eight hours, the scores of the two radar data assimilation experiments are far greater than that of CON.Eight hours later, the 1 h rainfall amount becomes larger and the 1 h TS is higher.e 1 h TS at most time is higher for RAD and RADC than for CON.Compared with RAD, the 1 h TS of RADC for heavier rainfall is higher than that of RAD.In summary, RADC with radar re ectivity correction has the best simulation result of rainfall during this event.
e intensity of updraft is an important factor of the rainstorm system.Figure 11 shows cross sections of vertical velocity and potential equivalent temperature (θ e ) over the rainstorm area.ree hours later, the low values of θ e simulated in the experiments are located at the height of 3-6 km above the rainstorm area.Convective unstable strati cation is found below 4 km, and θ e decreases with increasing height.In CON, the ascending motion over the heavy rainfall area is weak and the maximum radar echo is only 30 dBZ (Figure 11(a)).RAD and RADC with radar data assimilation create two convection cells with the scale of 10-20 km.
ere exist strong updrafts within the convection cells.
e ascending motion extends from the lower levels to the upper troposphere above 12 km (Figures 11(b) and 11(c)).Compared with that in RAD, the convective cells simulated in RADC are stronger and the maximum radar echo is over 45 dBZ.e simulated convection cells correspond to the locations of the two rainstorms at this time shown in Figure 5(a).is result indicates that, by assimilating radar re ectivity data, the convective cells in the heavy rainfall event can be well simulated.
e structure of the rainstorm system simulated in RADC using the correction method is a little better than that in RAD.

Summary and Conclusions
A heavy Meiyu frontal rainfall event occurred in the Yangtze-Huaihe River Basin from 4 to 5 July 2014 is studied with observed radar data assimilation.A new correction method is developed to correct three-dimensional radar reflectivity data based on hourly surface rain gauge observations.ree numerical experiments are carried out using LAPS and WRF to evaluate the impact of radar data assimilation on the simulation of this heavy rainfall event, and the performance of the correction method is evaluated by assimilating radar reflectivity observations into LAPS.e results show that the assimilation of radar data can effectively correct the location errors in the background and improve the simulation of the heavy rain.e heavy rainfall simulated with assimilation of radar data derived from hourly surface rain gauge observations is the closest to observations in terms of the location, intensity, and evolution of rainfall with time.e overestimation of rainfall in the other two experiments is successfully corrected due to the assimilation of radar data.
e southwesterly low-level jet stream in the rainstorm area is enhanced by assimilation of radial wind data, and the correction of radar reflectivity is helpful to increase the humidity over the rainstorm area.e initial background field adjusted by radar data assimilation provides a more realistic environmental condition for the development and maintenance of the rainstorm system.Radar echoes in the rainstorm area are enhanced by using the correction method, which corrects radar reflectivity.
e analyzed hydrometeors and humidity increments are mainly concentrated in the area where radar echoes are increased.e positive humidity disturbance might be able to counteract the negative buoyancy associated with the introduction of hydrometers and evaporative cooling and thus is favorable for the maintenance and development of rainstorm system.e convective cells in the heavy rainfall event can be well simulated with assimilation of radar reflectivity data, and the simulated structure of the rainstorm system by using the correction method is consistent with the evolution of observed rainfall.
Our conclusions here are mainly based on a single case study.Radar observations provide a promising source for better initialization in numerical studies of such kind of heavy rainfall events.More methods are needed to explore uncertain issues regarding how to better utilize radar observations and how to evaluate the impacts of radar data assimilation on numerical weather forecast.Our future work will focus on more case studies for the purpose to reach more reliable conclusions and further improve the ability of numerical models for the forecast of heavy Meiyu frontal rainfall in the Yangtze River Basin.

2. 1 .
Data Description.e types of operational Doppler radar in China include C-band and S-band, which are deployed in different provinces.

Figure 1 (
a) shows the model domain and locations of various types of Doppler radar sites.In the present study, 83 radars are used, including 7 SC radars, 13 CD radars, 5 CC radars, 6 CB radars, and 52 SA and SB radars.e first letter of the radar type acronyms means S-band or C-band.

Figure 1 :
Figure 1: (a) Locations of Doppler radar sites in China (black solid circle, hollow circle, solid square, hollow square, and × refer to SA/SB, CB, CD, CB, and SC, respectively), and (b) radiosonde (black hollow square) and surface meteorological sites (×).

Figure 2 :
Figure 2: Diagram of the correction method to correct radar re ectivity using hourly surface rain gauge observations.

Figure 3 :
Figure 3: Twenty-four-hour accumulated precipitation for the period 0000 UTC 4 July to 0000 UTC 5 July 2014 from observations (a) and simulations in (b) CON, (c) RAD, and (d) RADC (unit: mm).e black and red solid boxes represent the main rainfall belt and the key rainstorm area with 24 h accumulated precipitation amount greater than 100 mm, respectively.

Figure 5 :
Figure 5: Time-distance cross sections along line EF direction (shown in Figure 3(a)) of hourly rainfall averaged along line FH direction (Figure 3(a)) in the rainfall belt from 0000 UTC 4 July to 0000 UTC 5 July 2014: observations (a) and simulations in (b) CON, (c) RAD, and (d) RADC (unit: mm/hr).e red solid box represents temporal evolution of rainfall over the key rainstorm area (red box in Figure 3(a)).ehorizontal axis is the distance from line EF (Figure3(a)) (unit: 10 km).e vertical axis is the forecast hours.