Microphysical Characteristics of Three Convective Events with Intense Rainfall Observed by Polarimetric Radar and Disdrometer in Eastern China

Polarimetric radar and disdrometer observations obtained during the 2014 Observation, Prediction, and Analysis of Severe Convection of China (OPACC) field campaign are used in this study to investigate the microphysical characteristics of three primary types of organized intense rainfall events (meiyu rainband, typhoon outer rainband, and squall line) in eastern China. Drop size distributions (DSDs) of these three events on the ground are derived from measurements of a surface disdrometer, while the corresponding three-dimensional microphysical structures are obtained from the Nanjing University C-band polarimetric radar (NJU-CPOL). Although the environmental moisture and instability conditions are different, all three events possess relatively high freezing level favorable for warm-rain processes where the high medium to small raindrop concentration at low levels is consistent with the high surface rainfall rates. Convection is tallest in the squall line where abundant ice-phase processes generate large amounts of rimed particles (graupel and hail) above the freezing level and the largest surface raindrops are present among these three events. The storm tops of both the typhoon and meiyu rainbands are lower than that in the squall line, composed of less active ice processes above the freezing level. The typhoon rainrate is more intense than that of meiyu, enhanced by higher coalescence efficiency. A revised generalized intercept parameter versus mass-weighted mean diameter (Nw-Dm) space diagram is constructed to describe the DSD distributions over the three events and illustrate the relative DSD positions for heavy precipitation. DSDs of these intense rainfall convections observed in this midlatitude region of eastern Asia somewhat represent the typical DSD characteristics in low latitudes, suggesting that the parameterization of microphysical characteristics in eastern China in numerical models needs to be further investigated to improve rain fall forecasts in these heavy rainfall events.


Introduction
Rainfall from convective systems accounts for a significant portion of the total precipitation during the monsoon season in eastern China [1,2].Flash floods associated with these extreme precipitation events occur frequently in eastern China, causing severe damage to properties and loss of lives.To improve quantitative precipitation forecast (QPF) skills for these convection systems, better understanding of precipitation structures and microphysical characteristics is essential to improve physical and microphysical representations in numerical weather prediction models.For example, different microphysics schemes in the Weather Research and Forecasting (WRF) model perform differently in certain rainfall events [3][4][5], which is mainly related to the different representations of microphysical processes.Therefore, observations of the microphysical characteristics in precipitation systems are crucial toward refining microphysics schemes in numerical models and ultimately improving QPF.
Rainfall drop size distributions (DSDs) can be obtained from disdrometer observations [6][7][8].The three-dimensional precipitation structures, particularly the vertical structures associated with the evolution of microphysical processes and hydrometeor distributions, are mainly inferred through polarimetric radar observations [9][10][11].Observations from disdrometers in conjunction with polarimetric radars have been commonly used for investigating DSDs and precipitation structures in different climate regimes [12][13][14][15].The combined analyses have also provided insights into the microphysical structures of certain precipitation cases (e.g., squall lines, typhoon rainbands, etc.) in eastern China [16,17].However, the quantitative comparison of observed microphysical characteristics among these summertime convection systems which commonly bring heavy precipitation amounts have yet to be discussed in detail.
The microphysical characteristics of convective clusters were first classified into maritime and continental types [12].Subsequent studies revealed that the variation of environmental conditions, in addition to geographic locations, also affects the microphysical properties of convective systems.The precipitation in warmer and more humid maritime (drier continental) environments is mainly composed of raindrops with higher (lower) number concentration and smaller (bigger) mean drop size [2,14,[18][19][20].Recently, Dolan et al. [15] (hereafter, DO18) classified the global DSD characteristics into six groups based on the principal component analysis to twelve disdrometer datasets (Supplementary Materials) across three latitude bands on Earth.In their study, a typical clustering of DSDs (group 5) as the signature of robust warm-rain processes (high rain rate, large number concentration, and moderate mean size) was shown in the low latitudes but rarely emerged in the midlatitudes.On the other hand, another distinct group of DSDs (group 6) associated with ice-based convection (large mean drop size but small number concentration), was frequently emerged in the midlatitudes.It indicated that the distribution of DSDs vary in different latitudes, and the two DSD groups above seem to be helpful to distinguish the DSD characteristics of the midlatitudes from the low latitudes.
During the 2014 Observation, Prediction, and Analysis of Severe Convection of China (OPACC) field campaign in eastern China (in the midlatitude of Asian monsoon region) [16], several precipitation events were captured by multiple instruments.Three heavy precipitation events [meiyu rainband (meiyu), typhoon outer rainband (typhoon) over land, and squall line], were observed and selected for the study.The cases belong to three primary types of convective systems which produce the majority of extreme rainfall events in China [21].Chen et al. [22] has documented the radar quantitative precipitation estimation (QPE) performances of the three events, the effects of microphysical features on the radar QPE results near surface has also been discussed.This study expands their work to compare the microphysical features among the three precipitation events, for a better understanding of the similarity and differences of microphysical characteristics in these three types of intense convective systems in eastern China.The variation of DSD characteristics for these heavy precipitation events occurring is also summarized.
This paper is organized as follows: data and methodology are introduced in Section 2. The results are described in Section 3. The microphysical characteristics are discussed in Section 4, and the conclusions are given in Section 5.

The NJU-CPOL Radar and Disdrometer
The radar data were collected from the Nanjing University C-band polarimetric radar (NJU-CPOL) during the OPACC in summer 2014.The data quality control (QC) procedures for the polarimetric radar variables are the same as in Chen et al. [22].Firstly the system biases for the radar reflectivity (Z H ) and differential reflectivity (Z DR ) are calibrated using a metal sphere experiment and the vertically pointing data in light rain, respectively.After eliminating the non-meteorological echoes, an enhanced linear programming method [23] is used for the recalculation of the specific differential phase (K DP ).Then the attenuation of Z H (Z DR ) are corrected based on the linear relation between the specific horizontal attenuation A H (the specific differential attenuation A DP ) and K DP [22].After QC, the radar data in each volume coverage pattern (VCP), consisting 14-elevation Plan Position Indicator (PPI) scans, are interpolated onto constant altitude PPI (CAPPI) in a Cartesian grid using the NCAR REORDER software package with 1.0 km horizontal and 0.5 km vertical grid spacing starting from 0.5 km above ground level (AGL).
A second-generation OTT PARSIVEL disdrometer (OTT-2), located at ZhuangMu (ZM) station ~20 km north of the NJU-CPOL, is used to collect DSDs.The OTT-2 measures 32 bins of diameter from 0 to 25 mm and 32 bins of fall speeds from 0 to 22.4 m s −1 [24].The first two size bins (0.062 and 0.187 mm) are not used because of the low signal-to-noise ratios and the smallest size starts at the third bin with 0.312 mm in diameter.Spurious drops larger than 8 mm or 40% above or below the fall velocity-diameter relationship are also eliminated, to remove the impact of wind or splash on the disdrometer measurements [25][26][27].A size correction proposed by Battaglia et al. [28] has also been applied to the OTT-2 observations to minimize measurement errors.Raindrops below 1-mm equivalent sphere diameter are assumed to be spherical, where the ratio between the height and width (hereafter referring to axis ratio) is 1.While the axis ratios are assumed to change linearly from 1 to 0.7 for the raindrops with diameters between 1 and 5 mm.But for drops with diameters above 5 mm, the axis ratio is constant and set to 0.7.Previous studies have found that the OTT PARSIVEL disdrometer would underestimate small raindrops and resulted in relatively higher D m and lower N w [29,30], which should be treated carefully when compared with the DSD results observed with other instruments, like two-dimensional video disdrometer (2DVD).

Raindrop Size Distribution Parameters
The integral rainfall parameters, including rain rate R mm h −1 and liquid water content (LWC, g m −3 ), are computed from the disdrometer observed raindrop size distribution: where L is the total number of bins (30) for the OTT-2, D i (mm) is the equivalent spherical raindrop diameter, ∆D I is the corresponding diameter interval (mm), V i m s −1 is the fall speed obtained with the theoretical formula from Brandes et al. [31], and N(I) mm −1 m −3 represents the corresponding number concentration of raindrops for bin I, respectively.The nth-order moment of the DSDs can be expressed as: where D (mm) represents the equivalent diameter and N(D) mm −1 m −3 is the number concentration of raindrops in a unit volume of air and unit size interval.The mass-weighted mean diameter D m (mm) equals the ratio of the fourth to the third moment of the size distribution: and the generalized intercept parameter N w mm −1 m −3 is computed as Bringi et al. [12]: where ρ w is the density of water (1.0 g cm −3 ).

Hydrometeor Identification Algorithm (HID)
The HID is a fuzzy-logic algorithm based on the construction of membership functions for each polarimetric variable to identify the hydrometeor type that dominates the radar resolution volume [32][33][34][35].The NJU-CPOL HID follows Dolan et al. [34] that considers the unique scattering characteristics of C-band polarimetric radar, such as the C-band resonance effects when drops with diameters of 5.5-7 mm exist in the radar volume.Five variables including temperature, Z H , Z DR , correlation coefficient (ρ hv ) and K DP are used to identify the hydrometeor types, temperature data is obtained from the Jiangning radiosonde located ~170 km east of the NJU-CPOL.Membership beta functions are the same as in Dolan et al. [34].Ten hydrometeor types are identified in the HID including: drizzle, rain, big drops, ice crystals, vertical aligned ice, aggregates, wet snow, high-density (HD) graupel, low-density (LD) graupel and hail.Limitations of fuzzy-logic based HID include: (1) it can only identify the dominant rather than all the hydrometeor types in one radar sample volume, and (2) the overlap of polarimetric boundaries among different hydrometeor types sometimes turns to erroneously identified particle species.Radar data quality can also limit the accuracy of the HID results [9].However, it remains the most reliable method for the identification of hydrometeor types by polarimetric radars.

Rain and Ice Water Content Estimation
The method for estimating rain water content using polarimetric radar observations varies at different heights, which is similar to Chang et al. [36].In the presence of mixed-phase precipitation, the difference reflectivity (Z DP , dB) method [37] Then M w and M i , which represent the total mass of liquid and solid hydrometeors in the radar resolution volume, are estimated using the reflectivity-water content (Z-M) relationships adopted from previous studies [16,36,37]: and where Z rain h is the linear form of Z rain H , and Z ice h is calculated as follows: Values of negative Z ice h can be obtained when the values of Z H and Z DR are high, it happens when the mean raindrop size is high and treated as pure rain.When only ice exists, the ice mass (M i ) is directly estimated via (8).When in pure rain regions, the liquid water content (M w ) is calculated from the C-G model retrieval method [38].Discontinuity may exist between the C-G method and Z-M relationship, while more accurate liquid water content can be obtained by using two parameters (Z H and Z DR ) from the C-G method [36].

Event Description and Environmental Conditions
Three precipitation events (meiyu, typhoon, and squall line) observed by the NJU-CPOL and an OTT-2 during the OPACC field campaign in 2014 in eastern China are selected for the study.Time periods during the passages of the three convection events over the OTT-2 and corresponding radar scanning sequences are provided in Table 1.The selected NJU-CPOL radar reflectivity CAPPIs at 2-km AGL are shown in Figure 1.The black arrows represent the moving directions for the three precipitation events, it indicates that the ZM OTT basically captured the main part of the three precipitation systems during their passage through this site.Values of Z H in the mature convective cells for the meiyu system are over 45   Table 1.Time periods for the analysis of the three weather systems.The environmental conditions of these events are depicted by the skew T-logP diagrams from the Jiangning radiosonde (Figure 2).Several parameters to characterize the environment are presented in Table 2.The convective available potential energy (CAPE) for the squall line case is ~3428 Jkg −1 , suggesting strong instability favorable for the development of deep convection [39].In contrast, the CAPEs are modest in the meiyu (940 Jkg −1 ) and typhoon rainband (1404 Jkg −1 ) events, limiting the maximum updraft and vertical development of convection.The 0-6 km vertical wind shears are moderate for the three events (13.9, 13.5, and 11.6 m s −1 ) (Table 2), which are consistent with the organization of multicellular convection results in longer lifetime and larger rainfall accumulation when compared with single-cell convection [40,41].As for the environment moisture, the air is almost saturated below 400 hPa (the red line in Figure 2b) before the passage of the typhoon rainband and is accompanied by a mean relative humidity ∼91% between 0.5 and 3.5 km.The environmental air is also moist at low levels in the meiyu event with a mean relative humidity between 0.5-3.5 km, ~84%.For the squall line event, the relative humidity is ∼59% in the same vertical layer, suggesting the driest environment at low levels among these three events.Precipitable water (PW) from the surface to 300 hPa is equivalent to the vertical accumulation of water vapor content, a good indicator of the environmental moisture field [42].Values of PW (Table 2) indicate that the moist environments of the three events are suitable for the intense rainfall, especially for the typhoon case (71.1 mm).On the other hand, all three events possess low LCLs (<1 km) and high freezing levels (>5 km) (Table 2).The presence of thick warm cloud layers (>4 km) are in favor of warm-rain processes [43].The environmental conditions of these events are depicted by the skew T-logP diagrams from the Jiangning radiosonde (Figure 2).Several parameters to characterize the environment are presented in Table 2.The convective available potential energy (CAPE) for the squall line case is ~3428 Jkg , suggesting strong instability favorable for the development of deep convection [39].In contrast, the CAPEs are modest in the meiyu (940 Jkg ) and typhoon rainband (1404 Jkg ) events, limiting the maximum updraft and vertical development of convection.The 0-6 km vertical wind shears are moderate for the three events (13.9, 13.5, and 11.6 m s-1) (Table 2), which are consistent with the organization of multicellular convection results in longer lifetime and larger rainfall accumulation when compared with single-cell convection [40,41].As for the environment moisture, the air is almost saturated below 400 hPa (the red line in Figure 2b) before the passage of the typhoon rainband and is accompanied by a mean relative humidity ~91% between 0.5 and 3.5 km.The environmental air is also moist at low levels in the meiyu event with a mean relative humidity between 0.5-3.5 km, ~84%.For the squall line event, the relative humidity is ~59% in the same vertical layer, suggesting the driest environment at low levels among these three events.Precipitable water (PW) from the surface to 300 hPa is equivalent to the vertical accumulation of water vapor content, a good indicator of the environmental moisture field [42].Values of PW (Table 2) indicate that the moist environments of the three events are suitable for the intense rainfall, especially for the    The microphysical characteristics during the passage of the three convective systems over the OTT-2 are investigated using the DSD time series and the instantaneous rain rate derived from the 1-min OTT-2 observations (Figure 3).The convective and stratiform types of samples are classified with the method used in Bringi et al. [12] and illustrated by the red and black bars on the top of the diagram, respectively.For the typhoon and squall line (Figure 3b,c The microphysical characteristics during the passage of the three convective systems over the OTT-2 are investigated using the DSD time series and the instantaneous rain rate derived from the 1min OTT-2 observations (Figure 3).The convective and stratiform types of samples are classified with the method used in Bringi et al. [12] and illustrated by the red and black bars on the top of the diagram, respectively.For the typhoon and squall line (Figure 3b    For further comparison of the surface DSDs among the three events, the DSD time series during the passage of precipitation systems are shown in Figure 4.For the meiyu event (Figure 4a), the instantaneous rain rate increases to the first peak at 11:25 and reaches the second peak at 11:55 (Figure 4a3) while D m increases from 1 mm to 2.3 mm (Figure 4a1) and the maximum value of LWC reaches 2 g m −3 (Figure 4a2).Values of log 10 N w ~4 and D m ~2.2 mm (Figure 4a1) when at the peak rain rate ( R ∼ 50 mmh −1 ), and their distributions fall between the maritime and continental convective clusters (gray rectangles) reported by Bringi et al. [12].As the instantaneous rain rate decreases, LWC reduced to near zero and the distributions of D m and N w move toward the bottom left and pass through the black dashed stratiform rain types fitting line [12], into the stratiform region.

Vertical Evolution at the OTT-2 Site
These three convective systems were also observed by the NJU-CPOL in conjunction with the ZM OTT-2.The combined observations are helpful in understanding the microphysical structures and of the precipitation in addition to the DSDs on the ground.The time-height plots of radar variables at the OTT-2 site are illustrated in Figure 5.There is one peak of Z and K near the surface in the meiyu event (Figure 5a1, a3), with the values ~ 45 dBZ and 0.9° km , respectively.And the maximum Z is less than 1.2 dB during the event.On the other hand, maximum radar reflectivity exceeds 50 dBZ below the 3km level in the typhoon system (Figure 5b1), with the corresponding Z less than 1.4 dB and the maximum K > 3.3° km .High (low) values of As for the typhoon event (Figure 4b), the maximum instantaneous rain rate (>100 mmh −1 ) around 17:50 (Figure 4b3) corresponds to a LWC of 5 g m −3 (Figure 4b2), the number concentration log 10 N w of 4.3 and a D m ∼2.0 mm (Figure 4b1), close to the characteristics of maritime convection.At the second rainfall peak after 18:00 in (Figure 4b3), the microphysical characteristics move toward the continental convective clusters with log 10 N w decreasing to 4 and D m exceeding 2.5 mm (Figure 4b1).Therefore, high precipitation intensity in the typhoon system can be accompanied by either high value of N w or D m .The rain rate and LWC on the trailing part of the typhoon rainband decrease rapidly, along with the further decrease of N w and D m .For the squall line, the initial distribution of N w versus D m locates within the stratiform region (Figure 4c1).Within the convective part of the system, D m increases rapidly to 3.5 mm while N w decrease slightly, showing a typical DSD characteristic of continental convection.When D m decreases and log 10 N w reaches 4 around the time at 20:12, the instantaneous rain rate exceeds 100 mm h −1 (Figure 4c3) and LWC is also peaked (Figure 4c2).The corresponding distribution of N w versus D m locates between the maritime and continental convective clusters, similar to the meiyu and typhoon systems, but the mean D m is much higher than the other two systems.Beyond the heavy precipitation region, D m become smaller and N w becomes lower, resulting in a lower LWC.
Results obtained from observations of the OTT-2 indicate that the squall line generally possesses the largest mean size of raindrops during the passage of the rainfall system.And the higher concentration of raindrops in the typhoon and squall line systems correspond to more intense rainfall rate when compared with the meiyu event.On the other hand, the distribution of N w versus D m can vary widely during the passage of convective systems, sometimes changing the DSD characteristics from continental convective clusters to maritime convective regions (the squall line event).And during the periods with peak instantaneous rain rate and LWC, the DSD characteristics are mainly represented as high N w and medium D m , indicating that heavy precipitation is resulted from a combination of high raindrop concentration and medium raindrops.However, due to the instrument limitation, observed N w (D m ) from the OTT-2 are lower (higher) than the true values, but the variation and comparison of the DSD characteristics above are basically reliable.

Vertical Evolution at the OTT-2 Site
These three convective systems were also observed by the NJU-CPOL in conjunction with the ZM OTT-2.The combined observations are helpful in understanding the microphysical structures and of the precipitation in addition to the DSDs on the ground.The time-height plots of radar variables at the OTT-2 site are illustrated in Figure 5.There is one peak of Z H and K DP near the surface in the meiyu event (Figure 5a1, a3), with the values ~45 dBZ and 0.9 • km −1 , respectively.And the maximum Z DR is less than 1.2 dB during the event.On the other hand, maximum radar reflectivity exceeds 50 dBZ below the 3km level in the typhoon system (Figure 5b1), with the corresponding Z DR less than 1.4 dB and the maximum K DP > 3.3 • km −1 .High (low) values of K DP (Z DR ) reflect that the typhoon event is mainly composed of small to median raindrops in this layer, while the rainfall rate and LWC are high, consistent with the OTT-2 observations on the ground.As for the squall line event, the maximum Z H in the convective core exceeds 55 dBZ, indicating the most intense convection among these three events.High Z DR and K DP are also observed, suggesting the existence of large drops and high rain rate.The maximum Z DR ∼3 dB occurs first around 20:05 at the 3 km level, and decrease rapidly at 20:14, then the second peak value exists at surface around 20:30, which is consistent with the variation of D m in Figure 4c1.In addition, the maximum K DP >4 • km −1 exists at 20:14 at the 4km height, then the K DP core descends to the ground, corresponding to the high rainfall rate at surface (Figure 4c3).

Composite Vertical Structure of Radar Variables
To further compare the vertical structures and investigate microphysical processes over the three systems, the contoured frequency by altitude diagrams (CFADs) [44] and mean profiles of Z H , Z DR , and K DP in the convective regions (based on Steiner et al. [45]) around the OTT-2 site are portrayed in Figure 6.The selected domain of datasets is a 20 km × 20 km box (gray box in Figure 1) centered on the OTT-2 site, to investigate the precipitation characteristics around the disdrometer.For the convection in the meiyu event, the major distribution (defined as the contoured frequency over 30%) of Z H is between 35 and 47 dBZ near the surface (Figure 6a1).Distributions of Z DR and K DP are relatively narrow with maximum Z DR < 1.5 dB and K DP < than 2 • km −1 located below the freezing level.The 30 dBZ echo top is a proxy of the convective intensity [46,47], while it is basically below the −10 • C level for the convection in the meiyu system.The downward increase of reflectivity and decrease of Z DR above the −10 • C level are evident in the meiyu event, representing the aggregation of ice crystals [9,48].Meanwhile, the increase in Z H and Z DR is also clear around the freezing level, suggesting the melting of aggregates [9,49,50].On the other hand, the major distribution of Z H is between 35 and 51 dBZ near the ground for the typhoon event (Figure 6b1), the maximum K DP >3 • km −1 and Z DR ~2 dB, which reflects more intense convection in the typhoon system when compared to meiyu.In addition, Z H and K DP increase more rapidly toward the ground below the freezing level in typhoon, indicating the higher efficiency for the growth of LWC through warm-rain processes compared to the meiyu system.The major distribution for the 30 dBZ echo top is around the −10 • C level and the processes of aggregation and melting are also evident in the typhoon system, similar to those in the meiyu event.2), respectively.

Composite Vertical Structure of Radar Variables
To further compare the vertical structures and investigate microphysical processes over the three systems, the contoured frequency by altitude diagrams (CFADs) [44] and mean profiles of Z , Z , and K in the convective regions (based on Steiner et al. [45]) around the OTT-2 site are portrayed in Figure 6.The selected domain of datasets is a 20 km × 20 km box (gray box in Figure 1) centered on the OTT-2 site, to investigate the precipitation characteristics around the disdrometer.For the convection in the meiyu event, the major distribution (defined as the contoured frequency over 30%) of Z is between 35 and 47 dBZ near the surface (Figure 6a1).Distributions of Z and K are relatively narrow with maximum Z < 1.5 dB and K < than 2° km located below the freezing level.The 30 dBZ echo top is a proxy of the convective intensity [46,47], while it is basically below the −10 C level for the convection in the meiyu system.The downward increase of reflectivity and decrease of Z above the −10 C level are evident in the meiyu event, representing the aggregation of ice crystals [9,48].Meanwhile, the increase in Z and Z is also clear around the freezing level, suggesting the melting of aggregates [9,49,50].On the other hand, the major distribution of Z is between 35 and 51 dBZ near the ground for the typhoon event (Figure 6b1), the maximum K >3 ° km and Z ~2 dB, which reflects more intense convection in the typhoon system when compared to meiyu.In addition, Z and K increase more rapidly toward the ground below the freezing level in typhoon, indicating the higher efficiency for the growth of LWC through warm-rain processes compared to the meiyu system.The major distribution for the 30 dBZ echo top is around the −10 C level and the processes of aggregation and melting are also evident in the typhoon  2), respectively.
For the squall line event, the major distribution of Z H is between 33 and 54 dBZ at low levels, the maximum frequency near the surface is ~50 dBZ, much higher than the other two systems.Maximum values of Z DR > 3 dB and K DP > 4 • km −1 below the freezing level imply both the existence of large-sized particles and high LWC of raindrops.The major distribution for the 30 dBZ echo top extends 9 km height and above the −20 • C level, which reflects that convection develops much deeper than the other two events.And large amount of high Z H (>30 dBZ) and low Z DR (~0 dB) exist between the −20 • C and 0 • C levels, indicating the existence of graupel formed through riming process [9,51].
The vertical profiles of mean Z H , Z DR , and K DP are used for the quantitative comparison of the vertical microphysical characteristics among the three systems (Figure 6d1,d2,d3).Firstly, mean Z H increase while Z DR decrease downwardly above the −10 • C level for meiyu and typhoon events.This is consistent with the growth of dendrites and the aggregation of ice crystals at these levels [50].As for the squall line, convection develops deeper due to a stronger updraft where more water vapor, liquid water and rimed aggregates can be transported into higher levels, producing large irregular ice particles with high (low) value of Z H (Z DR ).Z H and Z DR both increase rapidly downward between the −10 • C. and 0 • C levels in all three systems, where riming, aggregation and melting can coexist in this layer and affect the growth of particles [10,16].Below the freezing level, Z H , Z DR and K DP increase slowly downward to the surface in the meiyu event, indicating moderate growth of raindrops.In the typhoon event, Z H , Z DR and K DP increase toward the surface with the mean K DP close to 1 • km −1 near the surface, almost doubling that in the meiyu event.The rapid downward increase of these radar variables in conjunction with a high freezing level reflects both the growth of raindrop and LWC through warm-rain processes, which may be a result of the high efficiency of raindrop growth by collision-coalescence and the collection of cloud water in such a moist environment (Figure 2b).In contrast, Z H , Z DR and K DP of the squall line rapidly increase downward to the 4 km level due to the melting of large ice particles like hail or graupel.Then Z DR continues to increase while Z H and K DP remain constant down to the surface.It may indicate more complex warm-rain processes like the coexistence of size sorting, evaporation and coalescence and needs for further discussion.

HID Results
The derived HID can be used for further understanding the differences of microphysical characteristics and ice-phase processes like aggregation and riming represented by certain species of

HID Results
The derived HID can be used for further understanding the differences of microphysical characteristics and ice-phase processes like aggregation and riming represented by certain species of ice particles [9,10].The normalized vertical distributions of HID results in all three systems are portrayed in Figure 7. Ice crystals dominate the convective region above the −20 • C level in the meiyu event (Figure 7a).A large proportion of aggregates are identified between the −20 • C and −10 • C levels, indicating increased aggregation of ice crystals.Graupel reaches the −10 • C level and the fraction increases downward to the freezing level, indicating active riming process in this layer.Wet snow is an indicator of the melting process [9], and its dominant fraction below the 0 • C level suggesting the presence of a melting layer.Pure rain is classified at lower heights due to the melting of ice particles.Vertical distribution of hydrometeors in the typhoon event is similar to that of the meiyu event (Figure 7b) except for the higher fraction of aggregates in the upper levels and the presence of graupel at higher altitudes, reflecting more intense convection with a stronger updraft (as shown in Figure 6a1, b1).On the other hand, the presence of hail at the mid-levels of squall line reflects much more active ice-phase processes.Graupel presents above the 12 km level, where more water vapor, supercooled water and ice particles are transported there by strong updraft.The riming process plays an important role between the −20 • C and 0 • C levels with the fraction of graupel exceeding 40%.Hail melts completely before it reaches the ground (no hail was reported) and big drops are identified with small quantities at low levels, indicating the melting of large ice particles like hail and graupel.

Inferred Microphysical Processes
To quantitatively analyze the contributions of ice-phase and warm-rain processes to the surface precipitation among the three systems, separate profiles of ice and liquid water contents are shown in Figure 8. Profiles of the ice water content (IWC) in the meiyu and typhoon events are similar where the maximum values are around 0.25 g m −3 at about 6 km height, similar to that in the North American monsoon convection [19,48] in the mid-latitude.Both meiyu and typhoon events possess maximum LWC (>1.5 g m −3 ) near the surface five times more than the IWCs above the freezing level, indicating the importance of warm-rain processes below the melting layer in these two events.The low-level LWC in the typhoon event (∼2.0 g m −3 ) is higher than that in the meiyu event (LWC∼1.5 g m −3 ), and the faster downwardly increase of LWC represents more efficient raindrop growth from warm-rain processes.A three-fold higher IWC exists above the freezing level in the squall line event when compared with the other two systems, with the maximum value over 0.8 g m −3 around 6 km.It indicates more active ice-phase processes and higher production of graupel and hail above the freezing level in the squall line (Figure 7c).On the other hand, convection in the squall line also possesses high LWC (∼2.1 g m −3 ) at low levels, the growth of liquid particles is also related to the warm-rain processes in the deep warm cloud layer in this system.Compared to the meiyu and squall line, rainfall at surface from the typhoon event is more dependent on the growth of LWC through efficient warm-rain processes.And regardless of the different environmental conditions (e.g., CAPE and moisture) among the three events, the deep warm cloud layers are always favorable to the development of warm-rain processes.
Since Z H represents the intensity of precipitation and Z DR indicates raindrop mean size, the changes in Z H and Z DR over a 3-km layer (0.5-3.5 km AGL) of pure rain is able to reflect the variation of number concentration and size of raindrops, and they are used to infer the dominant warm-rain processes like coalescence, breakup, evaporation and size sorting in recent studies [20,52].Scatters of Z H and Z DR changes from the 3-km pure rain layer over the three systems are plotted in Figure 9, and each dot represents the result derived from one radar volume scan.It shows that the major samples from the typhoon event fall in the coalescence-dominant quadrant (positive ∆Z H and positive ∆Z DR ), indicating the importance of the coalescence process below the melting layer.Samples in the meiyu system are mainly close to the origin with near zero positive ∆Z H and ∆Z DR , indicating the warm-rain processes are less active.As for the squall line event, samples with high ∆Z DR and low ∆Z H fall in the coalescence-dominant quadrant, suggesting the slow growth of large-sized raindrops and the depletion of smaller particles, since value of Z H is sensitive to the variation of big raindrops, but it is less sensitive to the change of small drops when compared with Z DR .Several samples are also located in the top left quadrant, indicating the non-negligible role of evaporation and size sorting in this convection system, which develops in a relatively dry environment (Table 2).On the other hand, the red (blue) colors of these symbols, which represent the increase (decrease) of the mean LWC values from the higher layer to the lower layer, are also shown in Figure 9. Red symbols are almost located in the coalescence-dominant quadrant, indicating the increase of LWC, which may be related to the raindrop growth by collection of cloud water.While the blue symbols are mainly in the left quadrants, suggesting the decrease of LWC due to the breakup, evaporation and size sorting of the raindrops.The results are consistent with the former analysis using the changes in Z H and Z DR over the 3-km layer.
Remote Sens. 2019, 11, x FOR PEER REVIEW 13 of 21 of ice particles.Vertical distribution of hydrometeors in the typhoon event is similar to that of the meiyu event (Figure 7b) except for the higher fraction of aggregates in the upper levels and the presence of graupel at higher altitudes, reflecting more intense convection with a stronger updraft (as shown in Figure 6a1, b1).On the other hand, the presence of hail at the mid-levels of squall line reflects much more active ice-phase processes.Graupel presents above the 12 km level, where more water vapor, supercooled water and ice particles are transported there by strong updraft.The riming process plays an important role between the −20 C and 0 C levels with the fraction of graupel exceeding 40%.Hail melts completely before it reaches the ground (no hail was reported) and big drops are identified with small quantities at low levels, indicating the melting of large ice particles like hail and graupel.

Inferred Microphysical Processes
To quantitatively analyze the contributions of ice-phase and warm-rain processes to the surface precipitation among the three systems, separate profiles of ice and liquid water contents are shown freezing level in the squall line (Figure 7c).On the other hand, convection in the squall line also possesses high LWC (~2.1 g m at low levels, the growth of liquid particles is also related to the warm-rain processes in the deep warm cloud layer in this system.Compared to the meiyu and squall line, rainfall at surface from the typhoon event is more dependent on the growth of LWC through efficient warm-rain processes.And regardless of the different environmental conditions (e.g., CAPE and moisture) among the three events, the deep warm cloud layers are always favorable to the development of warm-rain processes.Since Z represents the intensity of precipitation and Z indicates raindrop mean size, the changes in Z and Z over a 3-km layer (0.5-3.5 km AGL) of pure rain is able to reflect the variation of number concentration and size of raindrops, and they are used to infer the dominant warm-rain processes like coalescence, breakup, evaporation and size sorting in recent studies [20,52].Scatters of Z and Z changes from the 3-km pure rain layer over the three systems are plotted in Figure 9, and each dot represents the result derived from one radar volume scan.It shows that the major samples from the typhoon event fall in the coalescence-dominant quadrant (positive ∆Z and positive ∆Z ), indicating the importance of the coalescence process below the melting layer.Samples in the meiyu system are mainly close to the origin with near zero positive ∆Z and ∆Z , indicating the warm-rain processes are less active.As for the squall line event, samples with high ∆Z and low ∆Z fall in the coalescence-dominant quadrant, suggesting the slow growth of largesized raindrops and the depletion of smaller particles, since value of Z is sensitive to the variation of big raindrops, but it is less sensitive to the change of small drops when compared with Z .Several samples are also located in the top left quadrant, indicating the non-negligible role of evaporation and size sorting in this convection system, which develops in a relatively dry environment (Table 2).On the other hand, the red (blue) colors of these symbols, which represent the increase (decrease) of the mean LWC values from the higher layer to the lower layer, are also shown in Figure 9. Red symbols are almost located in the coalescence-dominant quadrant, indicating the increase of LWC, which may be related to the raindrop growth by collection of cloud water.While the blue symbols are mainly in the left quadrants, suggesting the decrease of LWC due to the breakup, evaporation and size sorting of the raindrops.The results are consistent with the former analysis using the changes in Z and Z over the 3-km layer.

Discussion
Based on the microphysical analyses above, Figure 10 presents the variation of DSD

Discussion
Based on the microphysical analyses above, Figure 10 presents the variation of DSD characteristics of the three precipitation events in the Nw-Dm space.The shift of D m and N w values due to the limitation of OTT-2 have been adjusted according to Wen et al. [30].The colored ellipses depict the distribution of DSDs from the convective regions in the three cases as well as the distribution of DSDs related to heavy precipitation.The surface DSDs are composed of raindrops with relatively high N w but small D m (blue ellipse) in the convective region of the meiyu event, indicating the overall high number concentration and small mean size of raindrops there.As for the typhoon case, raindrops are generally with slightly larger mean size and even higher number concentration at surface, and the DSD distributions finally shift from the blue ellipse to the green one.On the other hand, the wide distribution of D m and N w in the squall line event (yellow ellipse) suggesting the large variation of DSD characteristics at surface.DSDs distributions of meiyu and typhoon are more close to the maritime convective clusters reported by Bringi et al. [12].But for the squall line case, it generally locates between the maritime and continental convective clusters.And the surface DSDs with intense rainfall (pink ellipse) are mainly shown the combination of high N w and medium D m , indicating that heavy precipitation is resulted from a combination of high raindrop concentration and medium raindrops.
Our results show that the overall surface DSD distribution of the three precipitation events are closer to the characteristics of warm-rain convection type (group 5 in DO18, dashed orange contour) rather than the ice-based convection type (group 6 in DO18, dashed purple contour).Since the warm-rain (ice-based) convection type of DSD is distinct in the low latitudes (midlatitudes) according to the statistical results of DO18.It indicates that these intense rainfall convections initiated in this midlatitude region of eastern Asia somewhat represent the typical DSD characteristics in low latitudes, suggesting the unique microphysical characteristics here.On the other hand, the ice-based convective DSD features can also be shown in certain periods of the squall line case.
We know that the results of surface DSD characteristics in these precipitation systems should be linked with the vertical microphysical processes and environment conditions.Convections from the meiyu event are in the relatively moist and stable environment, which are relatively shallow but with active warm-rain processes like coalescence.While for the typhoon case, environment is equally stable but with more low-level moisture, convections are intensified at warm cloud layers with higher Z H and K DP due to the growth of raindrop mean size and number concentration via collision-coalescence.As for the squall line event in a drier and more unstable environment, convection develops much deeper, ice-phase processes like aggregation and riming are more active above the freezing level, resulting in higher fraction of graupel and hail.The melting of the large ice particles leads to large raindrops at upper warm cloud layer, while the evaporation and even size-sorting processes can cause low number concentration and high mean size of raindrops at surface in certain periods.On the other hand, the three summertime precipitation events are all with equally deep warm cloud layers, warm-rain processes are essential in determining the potential for intense rainfall.

Conclusions
To describe the microphysical characteristics over different summertime convection systems with intense rainfall in eastern China, data from the combined observations of polarimetric radar and disdrometer during the OPACC field campaign is used in this study.The precipitation structures and DSD characteristics of three precipitation events belonging to different types of convective systems (meiyu rainband, typhoon outer rainband and squall line) are compared.The variation of DSD characteristics in different microphysical processes and environmental conditions are discussed as well.
The characteristics of the raindrop size distribution (DSDs) of these three precipitation systems are compared using data from an OTT-2 disdrometer, while the environmental conditions are presented from the observation of a radiosonde.Results show that the mean raindrop size in the squall line event is larger than the other two cases, indicating the unstable air and a relatively dry environment of this event promotes larger drops at the surface.Nevertheless, short-term extreme rainfall in this convection system is also produced by a large number of smaller raindrops.As for the meiyu and typhoon events in relatively moist environmental conditions, the microphysical characteristics are closer to the maritime precipitation, the number concentrations of surface rainfall

Conclusions
To describe the microphysical characteristics over different summertime convection systems with intense rainfall in eastern China, data from the combined observations of polarimetric radar and disdrometer during the OPACC field campaign is used in this study.The precipitation structures and DSD characteristics of three precipitation events belonging to different types of convective systems (meiyu rainband, typhoon outer rainband and squall line) are compared.The variation of DSD characteristics in different microphysical processes and environmental conditions are discussed as well.
The characteristics of the raindrop size distribution (DSDs) of these three precipitation systems are compared using data from an OTT-2 disdrometer, while the environmental conditions are presented from the observation of a radiosonde.Results show that the mean raindrop size in the squall line event is larger than the other two cases, indicating the unstable air and a relatively dry environment of this event promotes larger drops at the surface.Nevertheless, short-term extreme rainfall in this convection system is also produced by a large number of smaller raindrops.As for the meiyu and typhoon events in relatively moist environmental conditions, the microphysical characteristics are closer to the maritime precipitation, the number concentrations of surface rainfall DSDs are high while

21 3. 1 . 1 .
dBZ when almost reaching the OTT-2 site around 11:50 local standard time (LST) (LST = UTC + 8) on 12 July, 2014.The maximum Z H for the typhoon rainband exceeds 50 dBZ from 16:50 to 17:50 LST on 24 July, 2014, while it moved westward.As for the squall line, maximum radar reflectivity also exceeds 50 dBZ when passing through the OTT-2 site at 20:30 LST on 30 July, 2014.Remote Sens. 2019, 11, x FOR PEER REVIEW 5 of Event Description and Environmental Conditions Three precipitation events (meiyu, typhoon, and squall line) observed by the NJU-CPOL and an OTT-2 during the OPACC field campaign in 2014 in eastern China are selected for the study.Time periods during the passages of the three convection events over the OTT-2 and corresponding radar scanning sequences are provided in Table 1.The selected NJU-CPOL radar reflectivity CAPPIs at 2km AGL are shown in Figure 1.The black arrows represent the moving directions for the three precipitation events, it indicates that the ZM OTT basically captured the main part of the three precipitation systems during their passage through this site.Values of Z in the mature convective cells for the meiyu system are over 45 dBZ when almost reaching the OTT-2 site around 11:50 local standard time (LST) (LST = UTC + 8) on 12 July, 2014.The maximum Z for the typhoon rainband exceeds 50 dBZ from 16:50 to 17:50 LST on 24 July, 2014, while it moved westward.As for the squall line, maximum radar reflectivity also exceeds 50 dBZ when passing through the OTT-2 site at 20:30 LST on 30 July, 2014.

Figure 1 .
Figure 1.The selected constant altitude Plan Position Indicator (CAPPI) scans of radar reflectivity at 2-km AGL from the Nanjing University C-band polarimetric radar (NJU-CPOL) radar in (a) meiyu, (b) typhoon, and (c) squall line, respectively.The solid circle represents the location of the ZhuangMu (ZM) disdrometer, which is about 20 km to the NJU-CPOL (solid triangle).The moving directions of the three precipitation events are shown with black arrows.And the 20 km × 20 km rectangle centered round the disdrometer is shown with the gray box.

Figure 1 .
Figure 1.The selected constant altitude Plan Position Indicator (CAPPI) scans of radar reflectivity at 2-km AGL from the Nanjing University C-band polarimetric radar (NJU-CPOL) radar in (a) meiyu, (b) typhoon, and (c) squall line, respectively.The solid circle represents the location of the ZhuangMu (ZM) disdrometer, which is about 20 km to the NJU-CPOL (solid triangle).The moving directions of the three precipitation events are shown with black arrows.And the 20 km × 20 km rectangle centered round the disdrometer is shown with the gray box.

Figure 2 .
Figure 2. The skew T-LogP diagrams and vertical wind profiles of (a) meiyu, (b) typhoon, and (c) squall line from the Jiangning radiosonde.The black curve is the ascending path of a surface-based parcel, while the blue and red lines are temperature and dew point temperature profiles, respectively.

Figure 2 .Table 2 .
Figure 2. The skew T-LogP diagrams and vertical wind profiles of (a) meiyu, (b) typhoon, and (c) squall line from the Jiangning radiosonde.The black curve is the ascending path of a surface-based parcel, while the blue and red lines are temperature and dew point temperature profiles, respectively.Table 2. Environment parameters obtained from radiosonde of the three systems.
), high concentrations [ N(d) ∼ 10 4 m −3 mm −1 ] of small drops (d < 1 mm) are observed, while the concentration of medium drops ( d ∼ 1-3 mm) reached 10 3 m −3 mm −1 with a few large raindrops (d > 3 mm) recorded.The maximum instantaneous rain rates (pink lines) in these two systems exceeded 100 mm h −1 , consistent with the high concentrations of small raindrops.Between these two systems, more large drops (d > 3 mm) are present in the squall line.On the other hand, fewer small and medium drops are observed in meiyu with almost no large raindrops measured, resulting in the lowest maximum instantaneous rain rate ( R ∼ 50 mmh −1 ) among the three systems.Remote Sens. 2019, 11, x FOR PEER REVIEW 7 of 21 and 3c), high concentrations [N(d)~10 m mm ] of small drops (d < 1 mm) are observed, while the concentration of medium drops (d~1-3 mm) reached 10 m mm with a few large raindrops (d > 3 mm) recorded.The maximum instantaneous rain rates (pink lines) in these two systems exceeded 100 mm h , consistent with the high concentrations of small raindrops.Between these two systems, more large drops (d > 3 mm) are present in the squall line.On the other hand, fewer small and medium drops are observed in meiyu with almost no large raindrops measured, resulting in the lowest maximum instantaneous rain rate (R ~ 50 mmh ) among the three systems.

Figure 3 .
Figure 3.Time series of drop size distributions (DSDs) and rain rate calculated from ZM disdrometer in (a) meiyu, (b) typhoon, and (c) squall line.The color shading represents the DSD in logarithmic units of mm m and the y axis on the left indicates the equivalent volume diameter (mm) of raindrops, while the pink line represents the instantaneous rain rate (mm h ) is consistent with y axis on right.The classified convective and stratiform samples are illustrated by the red and black

Figure 3 .
Figure 3.Time series of drop size distributions (DSDs) and rain rate calculated from ZM disdrometer in (a) meiyu, (b) typhoon, and (c) squall line.The color shading represents the DSD in logarithmic units of mm −1 m −3 and the y axis on the left indicates the equivalent volume diameter (mm) of raindrops, Remote Sens. 2019, 11, x FOR PEER REVIEW 9 of 21

Figure 4 .
Figure 4. Time series of (a1, b1, and c1)  ( ), (a2, b2, and c2) LWC( ), and (a3, b3, and c3) rain rate from: (a) meiyu, (b) typhoon, and (c) squall line, where time increases as the colored markers move from blue to red.Observed samples are shown with four kinds of shapes in different thresholds of rain rate.The two gray rectangles in (a1, b1, and c1) represent the maritime and continental convective clusters, while the black dash line represents the stratiform rain type fitting line, as reported by Bringi et al. [12].

Figure 4 .
Figure 4. Time series of (a1, b1, and c1) N w (D m ), (a2, b2, and c2) LWC(D m ), and (a3, b3, and c3) rain rate from: (a) meiyu, (b) typhoon, and (c) squall line, where time increases as the colored markers move from blue to red.Observed samples are shown with four kinds of shapes in different thresholds of rain rate.The two gray rectangles in (a1, b1, and c1) represent the maritime and continental convective clusters, while the black dash line represents the stratiform rain type fitting line, as reported by Bringi et al. [12].

Figure 6 .
Figure 6.CFADs of convective-only (a1, b1, and c1) Z , (a2, b2, and c2) Z , and (a3, b3, and c3) K data in (a) meiyu, (b) typhoon and (c) squall line.The shaded colors represent the normalized frequency of occurrence relative to the maximum absolute frequency from to 100% and the outer black closed curve represent the frequency of 1%.The gray dashed lines indicate the levels of 0 C, −10 C and −20 C as in Figure 5.The average profiles of the radar variables in the meiyu, typhoon and squall line cases are shown in d1, d2 and d3.The CFAD analysis is computed over the domain display in Figure 1 (gray box).

Figure 6 .
Figure 6.CFADs of convective-only (a1, b1, and c1) Z H , (a2, b2, and c2) Z DR , and (a3, b3, and c3) K DP data in (a) meiyu, (b) typhoon and (c) squall line.The shaded colors represent the normalized frequency of occurrence relative to the maximum absolute frequency from 5% to 100% and the outer black closed curve represent the frequency of 1%.The gray dashed lines indicate the levels of 0 • C, −10 • C and −20 • C as in Figure 5.The average profiles of the radar variables in the meiyu, typhoon and squall line cases are shown in d1, d2 and d3.The CFAD analysis is computed over the domain display in Figure 1 (gray box).

Figure 7 .
Figure 7. Normalized frequency of occurrence of hydrometeor classification by height (km) in (a) meiyu, (b) typhoon and (c) squall line.Horizontal gray lines indicate the 0 C, −10 C and −20 C levels, as in Figure 5.

Figure 8 .
Figure 8. Mean values of liquid (solid lines) water and ice (dashed lines) in the convective region of three precipitation systems.

Figure 8 .
Figure 8. Mean values of liquid (solid lines) water and ice (dashed lines) in the convective region of three precipitation systems.

Figure 9 .
Figure 9. Change in Z as a function of the change in Z over the 3-km rain shaft in the meiyu (triangles), typhoon (squares) and squall line (circles).The microphysical processes represented by each quadrant are annotated.The red (blue) colors of these symbols represent the increase (decrease) of the mean LWC values from the higher layer to the lower layer, respectively.

Figure 9 .
Figure 9. Change in Z DR as a function of the change in Z H over the 3-km rain shaft in the meiyu (triangles), typhoon (squares) and squall line (circles).The microphysical processes represented by

Figure 10 .
Figure 10.DSD characteristics of the three precipitation events in the (D ) space.The blue, green and yellow ellipses depict the distribution of DSDs from the convective regions in the meiyu, typhoon and squall line cases, while the pink one shows the distribution of DSDs related to heavy precipitation.Dashed orange and purple contours represents the DSDs from warm-rain convections (group 5) and Ice-based convections (group 6) in DO18, respectively.The two gray rectangles represent the maritime and continental convective clusters in Bringi et al. (2003), while the black dash line represents the stratiform rain type fitting line, the same as in Figure 4.The gray area shows the distribution of DSDs in the stratiform region.

Figure 10 .
Figure 10.DSD characteristics of the three precipitation events in the N w (D m ) space.The blue, green and yellow ellipses depict the distribution of DSDs from the convective regions in the meiyu, typhoon and squall line cases, while the pink one shows the distribution of DSDs related to heavy precipitation.Dashed orange and purple contours represents the DSDs from warm-rain convections (group 5) and Ice-based convections (group 6) in DO18, respectively.The two gray rectangles represent the maritime and continental convective clusters in Bringi et al. (2003), while the black dash line represents the stratiform rain type fitting line, the same as in Figure 4.The gray area shows the distribution of DSDs in the stratiform region.

Table 1 .
Time periods for the analysis of the three weather systems.

Table 2 .
Environment parameters obtained from radiosonde of the three systems.