Microphysical process of precipitating hydrometeors from warm-front mid-level stratiform clouds revealed by ground-based lidar observations

Mid-level stratiform precipitations during the passage of warm front were detailedly observed on two occasions (light and 15 moderate rain) by a 355-nm polarization lidar and water-vapor Raman lidar, both equipped with waterproof transparent roof windows. The hours-long precipitation streaks shown in the lidar signal ( X ) and volume depolarization ratio ( δ v ) reveal some ubiquitous features of the microphysical process of precipitating hydrometeors. We find that for the light rain case, precipitation that reaches the surface begins as mixed-phase hydrometeors fall out of a shallow liquid cloud layer at altitudes above the 0 ° C isotherm level, and the depolarization ratio magnitude of falling hydrometeors increases from the liquid-20 water values ( δ v < 0.10) to the ice/snow values ( δ v > 0.25) during the first 100 − 200 m of their descent. Subsequently, the falling hydrometeors yield a dense layer with an ice/snow bright band occurring above and a liquid-water bright band occurring below (separated by a lidar dark band) as a result of crossing the 0 ° C level. The ice/snow bright band might be a manifestation of local hydrometeor accumulation. Most falling raindrops shrink or vanish in the liquid-water bright band due to evaporation, whereas a few large raindrops fall out of the layer. We also find that a prominent δ v peak (0.10 − 0.40) always 25 occurs at an

Abstract. Mid-level stratiform precipitations during the passage of warm fronts were detailedly observed on two occasions (light and moderate rain) by a 355 nm polarization lidar and water vapor Raman lidar, both equipped with waterproof transparent roof windows. The hours-long precipitation streaks shown in the lidar signal (X) and volume depolarization ratio (δ v ) reveal some ubiquitous features of the microphysical process of precipitating hydrometeors. We find that for the light-rain case precipitation that reaches the surface begins as ice-phase-dominant hydrometeors that fall out of a shallow liquid cloud layer at altitudes above the 0 • C isotherm level, and the depolarization ratio magnitude of falling hydrometeors increases from the liquid-water values (δ v < 0.09) to the ice/snow values (δ v > 0.20) during the first 100-200 m of their descent. Subsequently, the falling hydrometeors yield a dense layer with an ice/snow bright band occurring above and a liquid-water bright band occurring below (separated by a lidar dark band) as a result of crossing the 0 • C level. The ice/snow bright band might be a manifestation of local hydrometeor accumulation. Most falling raindrops shrink or vanish in the liquid-water bright band due to evaporation, whereas a few large raindrops fall out of the layer. We also find that a prominent δ v peak (0.10-0.40) always occurs at an altitude of approximately 0.6 km when precipitation reaches the surface, reflecting the collisioncoalescence growth of falling large raindrops and their subsequent spontaneous breakup. The microphysical process (at ice-bright-band altitudes and below) of moderate rain resembles that of the light-rain case, but more large-sized hydrometeors are involved.

Introduction
An observation-based understanding of the microphysical processes of precipitation is essential for weather/climate modeling and predictions. Such processes are difficult to observe since they involve a variety of hydrometeor sizes, shapes and phases at different altitudes, all of which are affected by cloud dynamics (Aggarwal et al., 2016;Pfitzenmaier et al., 2018). In situ aircraft observations deliver data on the sizes and numbers of hydrometeors only for small sampling volumes at single altitudes at any given time during preplanned case studies (Barrett et al., 2019). Although lidar and radar can measure the time-resolved vertical profiles of bulk backscattering quantities, retrievals of the microphysical properties of hydrometeors require numerous assumptions (e.g., the hydrometeor shape and size distributions). Furthermore, in most cases, ground-based lidar cannot penetrate high enough to sample complete precipitating hydrometeor layers (only profiling the lower part of a layer) due to signal attenuation. Thus, information about their source clouds is usually not available (Sassen et al., 2005;Di Girolamo et al., 2012;Mega et al., 2012). There is also a lack of systematically observed lidar data on precipitation processes, because most lidar systems are not protected from precipitation. Cloud/precipitation radars are insensitive to small raindrops and droplets in cloud layers. Therefore, the microphysical processes of precipitation formation are not well understood thus far.
Satellite observations have revealed that cold clouds are the major source of liquid precipitation over land (Mülmen-städt et al., 2015). Heterogeneous ice formation pertinent to cold clouds is believed to lead to the generation of rain (Field and Heymsfield, 2015;Bühl et al., 2016;Pfitzenmaier et al., 2018). The ice formation process has been studied extensively by observing liquid-layer-topped ice virgae, because ground-based lidar and radar can reliably sample the entire height ranges of ice virgae and their apparent source cloud bases de Boer et al., 2011;Bühl et al., 2016Bühl et al., , 2019. In stratiform cloud layers at temperatures above −20 • C, precipitating bulk ice particles (ice virga) occurred after bulk liquid phases had formed overhead de Boer et al., 2011). This suggests that the heterogeneous nucleation of ice proceeds via the freezing of supercooled droplets de Boer et al., 2011). Our polarization lidar observations have revealed the detailed vertical structures of falling virgae and their supercooled liquid source cloud layers, indicating that the depolarization ratio values of falling hydrometeors increase rapidly with decreasing altitude on the top of the virgae (Cheng and Yi, 2020).
To study the microphysical processes that occur at altitudes ranging from the apparent source cloud base down to the near-surface during surface precipitation, a 355 nm polarization lidar and a water vapor Raman lidar at the Wuhan University atmospheric observation site were equipped with waterproof transparent roof windows. According to an artificial water-splashing experiment, water accumulation on the lidar roof windows yielded nearly height-independent lidar signal (X, range-corrected signal) attenuation, whereas neither the X vertical structure nor the profile of the volume depolarization ratio δ v (the magnitude and vertical structure) were altered. In addition, water accumulation on the roof windows hardly impacted the lidar-observed subcloud water vapor mixing ratio (q v ) profiles. This allows us to systematically observe precipitation processes (light and moderate rains). Based on our lidar observations obtained on two warm-front occasions, a complete microphysical process is revealed for precipitating hydrometeors pertinent to warm-front-related mid-level stratiform precipitation (the ice-nucleating processes are not covered). This paper first depicts the relevant instrumentation and methodology. Section 3 presents two light and moderate warm-front precipitation cases observed at our lidar site. The summary and conclusions are given in Sect. 4.

Lidar
Precipitating hydrometeor observations were obtained with two newly developed lidars equipped with waterproof transparent roof windows at the Wuhan University atmospheric observatory (30.5 • N, 114.4 • E; 73 m above sea level). The roof windows were designed to project out from the sur-roundings, avoiding a heavy water accumulation on the window glass during rainfall. The two lidars can simultaneously deliver the sequential profiles of the range-corrected signal X, volume depolarization ratio δ v and water vapor mixing ratio q v . All the observation sessions started with clear-sky conditions and ended when heavy precipitation occurred. This allowed us to capture the evolving layer structures of light and moderate precipitation events as well as their precursor clouds present over our mid-latitude site.

Polarization lidar
The polarization lidar has a configuration similar to our 532 nm system (Kong and Yi, 2015), but the transmitter employs a frequency-tripled Nd:YAG laser. It produces emissions of ∼ 150 mJ per pulse at 355 nm with a repetition rate of 30 Hz. A Brewster polarizer is added to improve the polarization purity of the transmitting laser (up to ∼ 10 000 : 1). After beam expansion, the beam with a divergence of 0.15 mrad is transmitted vertically into the atmosphere (to the zenith). The backscattered light is collected by a 20 cm Cassegrain telescope. The field of view (FOV) of the receiver is ∼ 1 mrad. After collimation, the elastically backscattered light passes an interference filter (with a 0.3 nm bandwidth centered at 355 nm) and is then incident on a polarization beam-splitter prism (PBS). To decrease the cross talk between the two orthogonal polarization channels, two additional polarizers are placed on the two output sides of the PBS. The light exiting from the two polarizers is focused onto two photomultiplier tubes (PMTs). The signals from the two PMTs are gathered by a PC-controlled twochannel transient digitizer (TR40-160, manufactured by Licel).
The raw lidar data are stored in both analog and photon counting modes with a range resolution of 3.75 m and a temporal resolution of 1 min. Based on a method originally proposed by Newsom et al. (2009) that was further developed by Zhang et al. (2014), the stored analog and photon-count data are glued to form a reasonable photon-count profile with a large dynamic range. For the cases in this study, the altitude range of signal gluing was ∼ 1.2-3.3 km. The range and temporal resolution of the processed photon-count profiles are 30 m and 1 min, respectively. The starting altitude of the lidar measurements is ∼ 0.3 km, determined based on the overlap of the laser and the field of view of the telescope. The altitude values referenced in this article are all relative to sea level.
The range-corrected lidar signal X is utilized to represent the backscattering intensity (returned laser power) of cloud particles and gravitationally falling hydrometeors (Ansmann et al., 2008). The volume depolarization ratio δ v , defined by the ratio of the perpendicular-to parallel-polarized backscatter coefficients, can be obtained from the two-channel lidar signals along with the relative gain of the parallel and perpendicular channels. The relative gain is determined in advance using a conventional method (Freudenthaler et al., 2009). The magnitude of the δ v value allows us to identify whether the dominant backscattering is attributed to ice crystals or water droplets in a given backscatter volume (Shupe, 2007). In general, liquid water droplets suspended in the atmosphere are nearly spherical and produce a very low depolarization ratio (close to zero) for single scattering at exactly 180 • , while ice crystals, which are usually nonspherical, generate a quite large depolarization ratio in the 180 • backscattering direction. For some mid-level stratiform precipitations, gravitationally falling hydrometeors form initially at altitudes above the 0 • C isotherm level. They fall often as icephase-dominant hydrometeors at subzero temperatures during their early descent. After the falling hydrometeors pass through the 0 • C isotherm level, the snowflake-to-raindrop (ice-to-raindrop) transition can yield a shallow layer of relatively smaller lidar echoes (a local X minimum) that is called "lidar dark band" (Sassen and Chen, 1995;Di Girolamo et al., 2012). The lidar dark band can be used to differentiate between the altitudinal regions with ice-containing particles above the dark band and pure liquid raindrops below the dark band.
It should be mentioned that the particle depolarization ratio δ p is conceptually a more suitable quantity in discriminating spherical and nonspherical particles (hydrometeors) in virga/cloud than the volume depolarization ratio δ v . But, the volume depolarization ratio δ v represents the more basic lidar measurement. In order to validly utilize the δ v magnitude in discriminating spherical and nonspherical depolarizations, we have examined the relationship between δ p and δ v . The δ p magnitude is a well-defined function of δ v , lidar backscattering ratio R and molecular depolarization ratio δ m (Cairo et al., 1999). The molecular depolarization ratio δ m has a value of ∼ 0.004 in terms of our lidar receiver bandwidth (0.3 nm) (Behrendt and Nakamura, 2002). Information about the R value range is available from a combined consideration of the earlier lidar measurements and our current observations on precipitation-related cloud/virga. The typical values of R for enhanced aerosol load are around 2 and for optically thin clouds up to around 10 (Lampert et al., 2010). The R values are ∼5-8 on the upper part of typical shallow (∼ 400 m thick) evaporating ice virgae (see Fig. 4 in Cheng and Yi, 2020). In this study, the R value should certainly be larger than 7 on the precipitation-related virga layer. Based on the analysis for the δ p expression (Cairo et al., 1999) for clouds and virgae, the particle depolarization ratio δ p has a quasi-linear dependence on the volume depolarization ratio δ v and a very weak dependence on lidar backscatter ratio R (when R ≥ 5). This favorable feature of the functional dependences allows us to utilize δ v in discriminating whether the dominant lidar backscattering is attributed to spherical or nonspherical particles in a given backscatter volume. If R min is the minimum of the R value range for the clouds/virgae of interest (e.g., R min = 7 for the precipitation-related virgae in this study), the discrimination criterion of spherical particles expressed by δ v (z) (equivalent to δ p < 0.1) takes the form (see Appendix A for mathematical derivation) The discrimination criterion of nonspherical particles expressed by δ v (z) (equivalent to δ p > 0.2) is given approximately by As noticed from the right-hand sides of Inequalities (1) and (2), the absolute differences between the discrimination threshold values expressed by δ p (0.1 and 0.2) and by δ v are small for clouds/virgae with R min > 7. The unambiguous cloud-phase discriminations based on the volume depolarization ratio δ v in earlier literature (Wang and Sassen, 2001;Intrieri et al., 2002;Shupe, 2007;Ansmann et al., 2009;Lampert et al., 2010) have confirmed the functional relationship between δ p and δ v mentioned above. This allows us to employ δ v with very little threshold-value change in discriminating whether the dominant lidar backscattering is attributed to spherical or nonspherical particles in a given backscatter volume. Specifically, at altitudes above the dark band, the δ v -based discrimination criteria are δ v < 0.09 for spherical water drops/droplets and δ v > 0.17 for ice crystals (based on the above discrimination criteria when R min = 7), while an enhanced depolarization ratio (δ v > 0.1) at altitudes below the dark band indicates the presence of large raindrops. We examined the multiple-scattering-induced depolarization ratio enhancements for an opaque cloud layer composed of dense spherical water droplets by putting a motorized iris on our polarization lidar system. It is indicated that for a receiver FOV of 1 mrad, the enhanced depolarization ratio δ v values due to multiple scattering increased from ∼ 0.03 at the X peak altitude to a maximum value of ∼ 0.27 at the weaksignal cutoff altitude with increasing penetration of laser light into the opaque water-droplet cloud layer. Note that for the same receiver FOV (∼ 1 mrad), the multiple-scatteringinduced depolarization ratio δ v values were all less than 0.04 within the laser light penetration range in a slightly dense water-droplet cloud layer (Hu et al., 2006). Combining the earlier multiple-FOV polarization lidar measurements (Hu et al., 2006) and our similar observations yields a suggestion that for the 1 mrad receiver FOV the multiple-scatteringinduced depolarization ratio values larger than 0.10 should result from an opaque water-droplet cloud layer (see Figs. 2 and 4 in Yi et al., 2021). In other words, for the 1 mrad receiver FOV, the vertical structure of hydrometeors and aerosols present above a dense water-droplet cloud layer with δ v values larger than 0.1 is undetectable by ground-based lidars. An artificial water-splashing experiment was performed on the lidar roof windows to examine the effects of water accumulation. A comparison of the lidar profiles with and without water accumulation on the lidar roof windows is given in  Enhanced lidar signal (X) and depolarization (δ v ) values at altitudes around 4.0 km resulted from an optically thick (opaque) water-droplet cloud layer, because there existed a high X value and near-zero δ v value (∼ 0.008) on the cloud base (∼ 3.9 km) (Wang and Sassen, 2001), and also there initially existed a monotonic rapid increase in both the values of X and δ v with increasing penetration of laser light into the layer. The cloud-related structures shown in both the X and δ v profiles were consistent before and after water splashing (particularly, cloud base altitudes). This comparison clearly shows that water accumulation on the lidar roof windows yielded nearly height-independent lidar signal (X) attenuation, and neither the cloud-related X vertical structure nor the profile of the volume depolarization ratio δ v (the magnitude and vertical structure) were altered. This result is physically reasonable.

Water vapor Raman lidar
The configuration of the water vapor Raman lidar used in this study is similar to our 45 cm aperture Raman system (Wu and Yi, 2017), but the current Raman lidar shares the same transmitter with our 355 nm polarization lidar depicted above. It detects inelastic Raman backscatter from water vapor at 407 nm and nitrogen molecules at 387 nm as well as detecting elastic backscattered light by using a 20 cm receiver telescope. The water vapor mixing ratio q v , which is defined as the mass ratio between water vapor and dry air in a given volume, can be obtained from the Raman signals representing water vapor and nitrogen molecules (Whiteman et al., 1992). The Raman lidar system was calibrated by corresponding local radiosonde measurements. A comparison analysis showed that the lidar-derived water vapor mixing ratio profiles agree well with the coincident radiosonde data (the relative deviation is less than 10 % when the water vapor field is horizontally homogeneous on a scale of ∼ 20 km). During the daytime, the water vapor Raman signal is quite noisy at high altitudes due to strong sky background light, so the water vapor mixing ratio profiles are available only at altitudes below ∼ 2 km. A similar artificial water-splashing experiment to that described above was performed on the water vapor Raman lidar roof window. Water accumulation on the roof window hardly had an impact on the obtained subcloud q v profiles.

Radiosonde
The radiosondes were launched at 08:00 LT (00:00 UTC) and 20:00 LT (12:00 UTC) every day from the Wuhan weather station (∼ 23.4 km away from our lidar site). Profiles of the pressure, temperature, relative humidity, wind speed, and direction from the near surface up to a height of 20-30 km were measured. The obtained radiosonde profiles were used to quantitatively determine the meteorological conditions pertinent to the precipitation events and their precursor clouds. The temperature measurement error was less than 1 • C, and the uncertainty in the relative humidity measurement is less than 5 % when the temperature was higher than 10 • C (Nash et al., 2011).

All-sky camera and rain gauge
The cloud photographs are recorded every 2 min by a groundbased all-sky camera located at our lidar site. A tippingbucket rain gauge is used to measure the precipitation rate at the surface. It has a sampling interval of 1 min. For each 0.1 mm of precipitation, the bucket tips and empties, yielding an output signal. The cloud layer was mostly characterized by a mixed phase and had subcloud ice virgae during the later descent ( Fig. 2a and b). After the subcloud virgae reached an altitude (∼ 3.0 km) that was lower than the 0 • C level (∼ 3.6 km), as measured by a conventional radiosonde at approximately 20:00 LT on 27 December at the Wuhan weather station (∼ 23 km away from our lidar site), falling raindrops (precipitation streaks in the X and δ v contour plots) that reached the ground were frequently observed beneath the 3 km altitude until 05:38 LT on 28 December 2017 when the lidar operation terminated. Long survival time of falling ice crystals at altitudes below the 0 • C level might be ascribed to cooling of the surrounding air during their evaporation and melting. The associated water vapor mixing ratio, q v , increased steadily with the descent of the cloud layers (Fig. 2c).
In particular, a high-concentration moisture layer appeared in the subcloud region during the rainfall event. This moisture layer resulted mainly from the evaporation of snow/ice particles and raindrops. Corresponding photographs of the sky taken by a ground-based camera at our lidar site are given at the top of Fig. 2. The light rain lasted for ∼ 8 h and yielded an accumulated rainfall amount of 2.6 mm (rain gauge data obtained at our lidar site). Interestingly, a humid aerosol layer also moved downward from ∼ 4.2 km at ∼ 20:00 LT on 25 December to ∼ 2.3 km at 20:00 LT on 27 December 2017, which appeared to be associated with the warm front. Figure 3 presents the radiosonde profiles that are pertinent to the warm-front cloud at different stages and during precipitation, together with the 1 h mean lidar profiles obtained during the radiosonde launches. The temporally varying cloud properties (e.g., falling cloud base, increasing cloud thickness and variable cloud types) between 20:00 LT on 26 December and 20:00 LT on 27 December 2017 coincided with the classical picture of preceding upglide clouds of an advancing warmfront system. Accordingly, a downgoing moist layer was observed strengthening and broadening with time during this period ( Fig. 3b and c). At the cloud base (except cirrus), the relative humidity over ice had values close to the relative humidity threshold of 84 % that is conventionally used to determine the cloud base heights (Wang and Rossow, 1995;Zhang et al., 2018). Furthermore, the radiosonde data showed that the southwesterly wind mostly prevailed at the cloud altitudes (Fig. 3d, e and f), and the air pressure at altitudes of ∼ 0-5 km dropped continuously by ∼3-5 hPa in the period (not shown here), which did belong to the typical warm-front features.

Associated meteorological conditions
The radiosonde released at 08:00 LT on 28 December 2017 provided measurements of the meteorological conditions when precipitation reached the surface, although the lidar measurements had already terminated (at 05:38 LT) ∼ 2 h earlier. As seen from Fig. 3b (red), the relative humidity reached a maximum of 98 % with respect to water in an altitude range of ∼ 3-4 km, immediately above the tops of the liquid precipitation streaks (at ∼ 3 km; see Fig. 2a and b). Water vapor at altitudes of ∼ 3-9 km was advected from the southwest, as seen in the wind component profiles (Fig. 3f, red). The high water vapor mixing ratios observed at alti-tudes below ∼ 3 km came from the evaporation of falling raindrops.

Microphysical process of precipitating hydrometeors for the light warm-front rain
The X and δ v precipitation streaks were visible in the period between 23:51 LT on 27 December and 05:36 LT on 28 December 2017 ( Fig. 2a and b). The streaks extended from the starting height (∼ 0.3 km) of the lidar measurements to an altitude of ∼ 2.88 km when surface precipitation occurred. A lidar dark band (X minimum) appeared persistently on the top of the precipitation streaks at an ∼ 2.88 km altitude except when the dark band was concealed by a drifting small-scale cloud (at 2.2-2.6 km altitudes during 04:18-05:36 LT on 28 December). This is consistent with earlier lidar observations of stratiform precipitation (Sassen and Chen, 1995;Demoz et al., 2000;Roy and Bissonnette, 2001;Di Girolamo et al., 2012). An inapparent local depolarization (δ v ) minimum was also persistently present at an altitude of ∼ 2.76 km, lying just ∼ 100 m below the dark-band minimum (Fig. 2b). The local δ v minimum represented the completion of the melting process of most falling ice/snow particles. Note that the δ v value decreased as a whole from the ice/snow (including partially melted large particles) values (0.10-0.34) at altitudes above the lidar dark band to the small liquid drop level (≤ 0.04, far less than the δ v -based discrimination threshold value of spherical particles when the lidar backscatter ratio R ≥ 5) at an altitude ∼ 100 m below the dark-band minimum. The lidar dark band definitely differentiates the altitude regions of precipitating ice-containing hydrometeors occurring above and liquid raindrops occurring below. Although the rainfall-induced water accumulation on the roof window of the lidar varied with time, the precipitation streaks and dark band were steadily reasonably displayed in the X and δ v time-height plots ( Fig. 2a and b). This is consistent with the result of our water-splashing experiment.
To further clarify the microphysical process of precipitating hydrometeors, two sets of representative lidar profiles (X, δ v and q v ) for the period that precipitation reached the surface (in Fig. 2) are plotted in Figs. 4 and 5. Figure 4 gives three 1 min X and δ v profiles from 01:12 to 01:14 LT on 28 December 2017 and a 1 h averaged q v profile centered at 01:13 LT on the same day. The lidar dark band appeared at a 2.88 km altitude at approximately 01:13 LT, while the local δ v minimum (< 0.04, far less than the δ v -based discrimination threshold value of spherical particles when R ≥ 5) was located at a 2.76 km altitude. These altitudes represent a typical lidar signature of the snowflake-to-raindrop transition for a variety of stratiform precipitation events. An ice-containing bright band (ice bright band hereafter) with δ v values ranging from ∼ 0.13 to ∼ 0.39 was visible at altitudes ∼ 3.0-3.45 km, just above the lidar dark band (Fig. 4); these altitudes correspond to the "relative lidar bright band" in the literature water vapor mixing ratio q v measured by a water vapor Raman lidar on 26-28 December 2017, which exhibited the passage of a warm front and the resulting hours-long light rain. A sliding average of 60 min was applied to the Raman lidar data. The precipitation streaks surrounded by magenta lines are zoomed in to show their details. Shown on the top of the figure are the corresponding photographs of the sky taken by a ground-based camera at our lidar site, with the third photograph exhibiting the sky illuminated by a 532 nm laser beam during the onset of rainfall. (Sassen and Chen, 1995;Di Girolamo et al., 2012). The ice bright band peaked on its bottom (∼ 3.0 km). It showed a variable vertical structure and intensity (in both X and δ v ) on the timescale of minutes, representing the presence of small-scale fluctuations in the precipitating ice crystals and snowflakes. A liquid-water bright band appeared as a layer of relatively large particle backscatter values, located at ∼ 1.50-2.76 km altitudes, just below the lidar dark band (Fig. 4). It is called "weak lidar bright band" in the literature (Sassen and Chen, 1995;Di Girolamo et al., 2012). The δ v values in the water bright band were ∼ 0.03-0.06, indicating that the enhanced lidar backscattering therein was caused mainly by high-concentration quasi-spherical raindrops with diameters ≤ 1 mm. The water bright band actually represents a major precipitation-related lidar backscattering layer in the liquidphase stage of the light precipitation event. The water bright band appeared to have a larger vertical extent (∼ 1.26 km) than that of the lidar ice bright band.
The lidar q v profile (Fig. 4c) shows an enhanced water vapor mixing ratio at altitudes from ∼ 1.7-3.4 km, indicating the subcloud evaporation of precipitating hydrometeors. In particular, q v was maximized (5.95 g kg −1 ) around the water bright band center (at ∼ 2.34 km), suggesting that this altitude was a primary subcloud evaporation region for this light warm-front precipitation event. Furthermore, the q v values in the water bright band increased as precipitation continued (Fig. 2c). Combining the vertical structures of X, δ v and q v in the water bright band (Figs. 2 and 4) yields the suggestion that most falling small-sized raindrops shrunk or vanished in the water bright band due to evaporation, whereas a small portion of large-sized raindrops survived via collisioncoalescence processes and fell out of the water bright band.  At altitudes below the water bright band, the precipitationrelated lidar backscattering (X) apparently weakened (Fig. 4a, in which the enhanced X values at altitudes from 0.3-0.7 km resulted from the boundary layer aerosols), indicating low-density raindrops there, whereas δ v first increased with decreasing height and then decreased after reaching a maximum (0.13-0.16) at an altitude of approximately 0.6 km (Fig. 4b). Here the magnitude and altitude variation of the lidar depolarization ratio δ v values allow us to identify where large-sized raindrops form and break up. Falling small-sized raindrops (equivalent diameter ≤ 1.0 mm) are quasi-spherical (Pruppacher and Klett, 1997) and yield small δ v values (generally less than 0.1), whereas falling large-sized raindrops (equivalent diameter > 2.8 mm) become nonspherical (with flat or hollow bottom in falling direction) (Pruppacher and Klett, 1997) and lead to large δ v values (larger than 0.1). In fact, prominent δ v peaks (∼ 0.1-0.4) at altitudes of approximately 0.6 km are always observed in the δ v profiles related to reaching-surface precipitation in the present lightrain case (Fig. 2). The δ v maxima at an altitude of ∼ 0.6 km are much larger than the typical values (< ∼ 0.07) observed by our 355 nm polarization lidar at approximately the same altitude during rainless days. Here we can exclude a possibility that the δ v maxima (∼ 0.1-0.4) at ∼ 0.6 km altitude resulted from multiple scattering by dense droplets around this altitude. As mentioned above, for the 1 mrad receiver FOV, a dense water-droplet cloud layer with the multiplescattering-induced depolarization ratio δ v values larger than 0.1 is optically opaque. In contrast to this situation, in our case, when the prominent δ v peak (∼ 0.1-0.4) around 0.6 km altitude occurred, the vertical structure of the precipitation streaks at altitudes far above 0.6 km (e.g., ice bright band, lidar dark band and lidar water bright band) was unambiguously detected by our polarization lidar, indicating that the enhanced depolarization ratios around 0.6 km altitude cannot be caused by multiple scattering from dense spherical water droplets therein. Furthermore, since most falling raindrops evaporated and vanished in the liquid-water bright band as indicated by the enhanced water vapor mixing ratio therein and rapidly decreasing lidar signal on the bottom of the water bright band, small droplets at altitudes below the water bright band were hardly dense enough to generate a strong multiple scattering with δ v ≥ 0.1. Therefore, our observational results suggest that sparse large raindrops that fall out of the water bright band with higher fall velocities further grow in size by collecting smaller raindrops along their fall paths. They grow to sizes at which spontaneous breakup occurs at an altitude of approximately 0.6 km. In brief, our lidar observations reveal for the first time (to our knowledge) the collisioncoalescence growth and subsequent spontaneous breakup of falling raindrops that actually take place in the natural atmosphere. They represent the posterior microphysical processes necessary for the reaching-surface precipitation production. Interestingly, the size maximization of falling raindrops as shown by the strongest nonspherical shapes (maximum depolarization ratio values) always appeared at an altitude of ∼ 0.6 km for a variety of mid-level stratiform precipitations (in light of our observations). Obviously, the explanation to this ubiquitous feature needs further observational and modeling efforts. As seen in Fig. 2b, the boundary layer aerosols had little impact on the δ v precipitation streaks. In addition, at altitudes below 1.5 km, the q v values decreased with increasing altitude, reflecting a normal altitude distribution of the boundary layer water vapor.
Based on the radiosonde temperature data obtained at approximately 20:00 LT on 27 December 2017 (Fig. 3a, orange), the 0 • C isotherm level was at an altitude of ∼ 3.6 km, and a warm-front-related inversion layer appeared just below the 0 • C level with a local temperature maximum (2.2 • C) at 3.33 km and a local minimum (1.0 • C) at 2.84 km. The lidar dark band (at 2.88 km, with a temperature of ∼ 1.0 • C) was located ∼ 720 m below the 0 • C level. In comparison with the results reported in the literature (Sassen and Chen, 1995;Demoz et al., 2000;Sassen et al., 2005;Di Girolamo et al., 2012), the observed ∼ 720 m distance of the dark-band minimum to the 0 • C level and the low dark-band temperature (∼ 1.0 • C) are somewhat peculiar for light-precipitation cases. In the current case, the melting process might be delayed by the temperature structure (with a small lapse rate) of the inversion layer. However, it should be mentioned here that the radiosonde launching site was ∼ 23.4 km away from our lidar site. Figure 5 presents three 1 min lidar X and δ v profiles displaying the time span from 02:30 to 02:32 LT on 28 December 2017 and a 1 h averaged lidar q v profile centered at 02:31 LT on the same day, depicting the microphysical process of precipitating hydrometeors for slightly strong precipitation that reached the surface during the light rain event. Although the water bright band and aerosol backscatter layer below the dark band became evidently weak compared to those seen in Fig. 4a (due to precipitation attenuation), the altitude of the dark-band minimum (2.85 km) was very close to that (2.88 km) obtained from Fig. 4a. The magnitude (∼ 0.03) and occurrence altitude (2.76 km) of the local δ v minimum were consistent with the corresponding values (less than 0.04 and 2.76 km, respectively) observed in Fig. 4b. Furthermore, the depolarization maxima (∼ 0.17) associated with reaching-surface precipitation still appeared at an altitude of approximately 0.6 km, which was also similar to that seen in Fig. 4. The observational facts confirm the result gathered from our water-splashing experiment in which the thin liquid-water layer on the roof windows of the lidars caused nearly altitude-independent attenuation on the X profiles and had no effect on the δ v profiles. The profile characteristics shown in Fig. 5 are mostly similar to those mentioned above for Fig. 4, but some newly emerging features need to be illustrated. Figure 5 exhibits an ice bright band stronger than the concurrent water bright band. This result is different from our observations obtained at ∼ 01:12 LT (Fig. 4) but is consistent with earlier lidar observations (Sassen and Chen, 1995;Di Girolamo et al., 2012). The ice bright band observed at approximately 02:30 LT had the X maxima at its bottom (at an altitude of ∼ 3.0 km) and a small vertical extent (∼ 0.21 km, from 3.00 to ∼ 3.21 km due to precipitation attenuation). The X maxima corresponded to the local minima of the depolarization ratio ( Fig. 5a and b). Interestingly, this inverse relationship between the backscatter and depolarization values on the bottom of the ice bright band is nearly ubiquitous in the precipitation lidar profiles obtained in the present case. Since the depolarization δ v showed moderate minima (∼ 0.08-0.10) at an altitude of ∼ 3.0 km (Fig. 5b), the ice bright-band maxima observed at approximately 02:30 LT might reflect backscattering from high-concentration partially melted large particles therein. On the band's altitudinal extension (from 3.06-3.21 km), the markedly enhanced depolarization values (∼ 0.17-0.34) indicate the presence of ice crystals and large snowflakes (Sassen and Chen, 1995;Di Girolamo et al., 2012). The water vapor mixing ratio q v showed slight enhancements at altitudes of ∼ 1.5-3.0 km at approximately 02:30 LT compared with that measured at approximately 01:12 LT.
As seen from the X and δ v precipitation streaks at altitudes below ∼ 1.5 km ( Fig. 2a and b), precipitation that reached the surface was intermittent. During periods without reachingsurface precipitation, our lidars were able to sample both a complete virga (from the rain to the snow regions) and a shallow mixed-phase cloud layer immediately above the virga under weak optical attenuation conditions. Such an example is shown in Fig. 6. The lidar profiles above the dark band clearly exhibit the typical structure characteristics of a liquidtopped mixed-phase cloud (a shallow liquid cloud layer and ice virga below) (see Fig. 6 in Wang and Sassen, 2001). The mixed-phase cloud top layer (at altitudes of ∼ 4.6 km) was of high X values and very low δ v values (∼ 0.01), while lower part of the cloud was characterized by significantly lower X values and higher δ v values (with a maximum up to ∼ 0.33). Furthermore, the cloud top layer had a maximum water vapor mixing ratio q v and a temperature of ∼ −8.5 • C (based on radiosonde data at ∼ 20:00 LT on 27 December). Combining with the schematic representation of commonly observed mixed-phase cloud layers (see Fig. 1 in Bühl et al., 2016), the current observations suggest that the cloud top layer should mainly be composed of liquid droplets (that were not dense enough to yield detectable multiple scattering), and the lower part of the cloud was mainly precipitating ice crystals (falling ice virga). The liquid-topped mixedphase cloud (a liquid cloud layer and ice virga below) (Bühl et al., 2016) might be fundamental monomers that constitute mid-level precipitating stratiform clouds. According to the expressions on the right-hand sides of Inequalities (1) and (2), the δ v -based discrimination threshold values were, respectively, 0.09 for spherical particles and 0.17 for nonspherical particles when the lidar backscatter ratio R had a value of 7 (the minimum of the R value range) on the upper part of the precipitation-related virga (Lampert et al., 2010;Cheng and Yi, 2020). Thus, the δ v magnitude of the falling virga increased from the liquid-water values of ∼ 0.02-0.07 (< 0.09) at an altitude of 4.38 km to the ice/snow values of ∼ 0.21-0.33 (> 0.17) at an altitude of 4.02 km. The falling ice crystals yielded a very weak ice bright band at an altitude of ∼ 3.0 km and then melted into liquid drops at an altitude of ∼ 2.76 km (the local δ v minimum). During their further descent, the liquid drops fully vanished due to evaporation, leaving a lidar-detectable rain virga (water bright band) without reaching-surface precipitation. In contrast to the situation during precipitation that reached the surface, no clear-cut δ v enhancement occurred at an altitude of approximately 0.6 km when there were only virgae suspended in air. Similar results were discerned for other lidar profiles shown in Fig. 2, in which a complete mixed-phase cloud layer could be detected.
During the light warm-front rain event, since the reachingsurface precipitations and virgae occurred alternately on a small timescale from a few minutes to tens of minutes and since their precipitation streaks had nearly the same darkband structures ( Fig. 2a and b), both reaching-surface precipitation and virgae would come from the same source cloud (because a warm-front cloud system is generally widespread and slowly varying). Reaching-surface precipitation (drizzle) arose when the precipitation rate was high below the shallow water-droplet-dominated cloud layer (apparent source cloud), while virgae without reaching-surface precipitation took place when the subcloud precipitation rate was slightly low. Therefore, the current lidar observations reveal the microphysical process of precipitating hydrometeors related to light warm-front rain. Both reaching-surface rainfall and virgae suspended in air began as ice-phase-dominant hydrometeors fell out of a liquid apparent source cloud layer at altitudes above the 0 • C isotherm level. The depolarization ratio magnitude of falling hydrometeors increased from the liquidwater values (δ v < 0.09) to the ice/snow values (δ v > 0.20) during the first 100-200 m of their descent. Subsequently, the falling hydrometeors yielded a dense layer with an ice/snow bright band occurring above and a liquid-water bright band occurring below (separated by a lidar dark band) as a result of crossing the 0 • C level. In the ice/snow bright band, large particles would form via the cold rain processes (riming and aggregation), because the broad size distributions of the pristine hydrometeors falling out of the apparent source cloud base could lead to local accretion. The production efficiencies of large particles would depend on the magnitude of the rain rate below the apparent source cloud base and size distributions of the pristine falling hydrometeors. The local depolarization minimum (δ v ≤ 0.04, far less than the δ v -based discrimination threshold value of spherical particles when R ≥ 5) was persistently observed immediately beneath (∼ 100 m below) the lidar dark-band minimum (X minimum). This displayed that the completion of the melting process of most falling ice particles took place at altitudes (hundreds of meters) below the 0 • C isotherm level. The liquid-water bright band (with a geometrical thickness of ∼ 1 km) just below the lidar dark band was characterized by enhanced X values and small δ v values. There existed a high-concentration moisture (large q v values) in this bright band. These features indicate that the liquid-water bright band resulted from gravitationally falling, dense evaporating liquid drops. In terms of the lidar-measured profiles during reaching-surface precipitation, at altitudes below the water bright band, the precipitation-related lidar backscattering apparently weakened, while δ v first increased with decreasing altitude and then decreased after reaching a prominent maximum at an altitude of ∼ 0.6 km. The lidar profiles for the virgae showed narrower and weaker water bright bands than those observed when precipitation reached the surface. Moreover, during virga occurrence, there was no perceptible depolarization enhancement at an altitude of ∼ 0.6 km. By combining the abovementioned lidar observations, a picture on the microphysical processes of falling hydrometeors in liquid-phase stage emerged. After going through the dark band, most falling raindrops shrank or vanished in the water bright band due to evaporation, whereas a few large raindrops survived and fell out of the water bright band when the rain rate below the apparent source cloud base was high enough. The large raindrops might come from both the complete melting of large falling ice/snow particles and collision-coalescence formation in the dense water bright band. Sparse, large raindrops with high fall velocities further grew in size by collecting smaller raindrops along their fall paths. At an altitude of ∼ 0.6 km, the large raindrops grew to the sizes at which spontaneous breakup could occur, yielding reaching-surface precipitation. When the rain rate below the apparent source cloud base was low, nearly none of the large raindrops fell out of the water bright band. Consequently, there were only virgae suspended in air (without reachingsurface precipitation). Figure 7 shows an example of moderate warm-front precipitation that occurred on 4 March 2019. Both the descending precursor clouds and the X and δ v precipitation streaks are generally similar to those seen in the first example (Fig. 2). The precursor clouds are cirrus (photo I in Fig. 7), altostratus (photo II) and altocumulus (photo III). The reachingsurface precipitation started just after the subcloud ice virgae reached an altitude (∼ 2.7 km) slightly lower than the 0 • C level (∼ 3.0 km). The δ v precipitation streaks show the upper portion (ice bright band) containing ice/snow particles (mostly δ v > 0.3) and the lower portion (water bright band and below) being composed of liquid drops (δ v ≤∼ 0.12 except for the δ v maxima that occurred due to raindrop-size growth at an altitude of ∼ 0.6 km). The δ v values in both the ice bright band and water bright band (Fig. 7b) were generally larger than their counterparts in the light-rain example (Fig. 2b), indicating that more large ice/snow particles and raindrops were involved in the moderate precipitation than in the light-precipitation event. Partially melted, large falling particles sometimes concealed the lidar dark band produced by the melting effect of most relatively small-sized particles in precipitating hydrometeors, making the band somewhat fuzzy (Fig. 7a). Accordingly, the altitude of the local δ v minimum (on the lidar dark band) became somewhat unsteady (Fig. 7b). The δ v maxima at an altitude of approximately 0.6 km (Fig. 7b) were apparently larger than those shown in Fig. 2b, indicating that more breakupsize raindrops formed via collision-coalescence processes therein than in the light-rain case. Specifically, the δ v maxima at an altitude of ∼ 0.6 km were as high as ∼ 0.27-0.35 at ∼ 23:38 LT, which corresponded well to the large rainfall rate of 3.2 mm h −1 measured from our rain gauge on the ground. The (apparent) source cloud for this moderate rain event was invisible by the lidars due to strong optical attenuation. Therefore, the following analysis was limited to the ice bright band and below. A strong southerly wind prevailed at altitudes of 0-12 km in light of the radiosonde data obtained at 20:00 LT on 4 March 2019. A high-concentration moisture layer appeared in the subcloud region at altitudes from ∼ 0.5 to ∼ 3.0 km during the rainfall event, indicating the subcloud evaporation of precipitating hydrometeors. The moderate rainfall lasted for ∼ 14 h, yielding an accumulated rainfall amount of 23.9 mm on the ground.

Associated meteorological conditions
The conventional radiosonde profiles associated with the moderate warm-front precipitation and its precursor clouds and the 1 h mean lidar profiles obtained during the radiosonde launches are plotted in Fig. 8. At ∼ 08:00 LT on 4 March 2019, the sky was nearly cloudless (Fig. 8d, blue), and high relative humidity occurred only at altitudes below 1.2 km (Fig. 8b, blue), while the northwesterly wind prevailed at altitudes from 1.7-10.5 km. This indicated that the warm front had not yet reached our lidar site. At ∼ 20:00 LT on 4 March, a moist layer occurred at altitudes ranging from ∼ 4.8 to 8.0 km with increased relative humidity over water of 80 %-95 % (Fig. 8b, green). An evaporating ice virga was observed at altitudes from ∼ 3.6-4.6 km ( Fig. 8d and e, green), just below the moisture layer peak. The apparent source cloud of the virga was invisible by lidars due to strong optical attenuation. A potential occurrence region for the apparent source cloud ranged in altitude from 4.8-6.0 km, where the relative humidity was larger than 90 % (Fig. 8b, green). The southerly wind prevailed at altitudes from 0-12 km (Fig. 8f, green), indicating that the moisture layer and altocumulus (photo III in Fig. 7) were precursors of the warm-front precipitation event. The radiosonde profiles obtained at 08:00 LT on 5 March 2019 showed the meteorological conditions during the moderate warm-front precipitation event after the lidar measurements had already terminated (at 00:51 LT on 5 March). As shown in Fig. 8b (orange), the relative humidity over water had values of 97 %-98 % at altitudes from 0-5.65 km, corresponding to a precipitation rate of ∼ 1.8 mm h −1 (rain gauge record) at approximately 08:00 LT on 5 March. The air pressure from the radiosonde data at Wuhan showed a persistent decrease (by ∼ 2-4 hPa at altitudes of ∼ 0-5 km) during the observational period from the precursor clouds to precipitation (between 08:00 LT on 4 March to 08:00 LT on 5 March) that reflected the warmfront passage.

Microphysical process of precipitating
hydrometeors for the moderate warm-front rain identical dark-band locations (the X minima is located at ∼ 2.04 km, and the local δ v minima is located at ∼ 1.96 km). The dark-band minima appeared ∼ 960 m below the 0 • C level at a radiosonde temperature of ∼ 6.0 • C. Such a long survival time of falling ice crystals at altitudes below the 0 • C level was due to cooling of the surrounding air during their evaporation and melting. At 22:22 LT, a weak X peak occurred at the dark-band altitudes with δ v values ranging from ∼ 0.21-0.29, indicating that partially melted large particles passed through the dark band. As seen from Fig. 9, the X and δ v precipitation streaks had complicated vertical structures at altitudes below the dark band and showed strong variations on the timescale of minutes. In particular, enhanced depolarization (0.07-0.12) occurred within the water bright band. These profile details confirm that large-sized particles sometimes fell out of the ice bright band during the moder- The radiosonde profiles quantitatively present the meteorological conditions pertinent to the warm-front clouds and precipitation.
ate warm-front precipitation event, concealing the lidar dark band produced by the melting effect of most relatively smallsized particles in precipitating hydrometeors. This effect appears to explain why the lidar dark band became fuzzy for the moderate precipitation event (Fig. 7). Note that the δ v maxima (∼ 0.2), which occurred at an altitude of approximately 0.6 km, were slightly larger in the moderate warm-front rainfall than those observed in the light-rainfall warm-front example. This suggests a larger concentration of raindrops of spontaneous breakup sizes around this altitude. The water vapor mixing ratio q v had values ranging from 3.4-4.4 g kg −1 at altitudes from 1.0-3.0 km (from the bottom of the water bright band to the 0 • C level). Figure 10 gives three 1 min lidar X and δ v profiles representing the period from 23:37 to 23:39 LT on 4 March 2019 and a 1 h averaged lidar q v profile centered at 23:38 LT on the same day; these profiles exhibit the vertical structures of the X and δ v precipitation streaks as well as the water vapor mixing ratio observed when the surface precipitation rate was highest (3.2 mm h −1 ) during the moderate warm-front precipitation event (yielding thick liquid-water accumulation on the roof windows of the lidars; see photo V in Fig. 7). The ice bright band, dark band and water bright band were roughly discernible in the three 1 min X profiles despite the considerable fluctuations that occurred on the timescale of minutes. Large ice/snow particles occurred on the ice bright band (at an altitude of approximately 2.5 km), because the δ v values were larger than 0.3 therein. The dark band located ∼ 700 m below the 0 • C level (3.0 km) had δ v values ranging from 0.13-0.19 and a temperature of 4.3 • C, reflecting that there were partially melted large particles present in the dark band. In the height range of the water bright band, the depolarization ratio increased from ∼ 0.04-0.06 at an altitude of approximately 2.09 km to ∼ 0.12-0.15 at an alti- tude of 0.9 km, indicating that more large raindrops formed via collision-coalescence processes therein than in the lightrainfall warm-front example (Figs. 4 and 5). The δ v maxima observed at an altitude of ∼ 0.6 km were as high as ∼ 0.27-0.35 corresponding well to the high measured rainfall rate of 3.2 mm h −1 (rain gauge record on the ground). As mentioned above, for the 1 mrad receiver FOV, if such large δ v values (∼ 0.27-0.35) came from the multiple scattering by a dense water-droplet cloud layer around 0.6 km altitude, the cloud layer would be optically opaque. It would conceal the vertical structure of the precipitation streaks at altitudes above 0.6 km. In contrast to this situation, as seen from Fig. 10, Figure 10. Three 1 min lidar X and δ v profiles covering the period from 23:37 to 23:39 LT on 4 March 2019 and a 1 h averaged lidar q v profile centered at 23:38 LT on the same day, exhibiting the vertical structure of the X and δ v precipitation streaks as well as the water vapor mixing ratio when the surface precipitation rate was highest (3.2 mm h −1 ) during the studied moderate warm-front precipitation event (yielding thick liquid-water accumulation on the roof windows of the lidars; see photo V in Fig. 7). the vertical structure of the precipitation streaks at altitudes above 0.6 km was clearly discerned by our ground-based polarization lidar, indicating that the enhanced depolarization ratios around 0.6 km altitude cannot be caused by multiple scattering from dense spherical water droplets therein. Furthermore, since most falling raindrops evaporated and vanished in the liquid-water bright band as indicated by the enhanced water vapor mixing ratio therein and rapidly decreasing lidar signal on the bottom of the water bright band, small droplets at altitudes below the water bright band were hardly dense enough to generate a strong multiple scattering with δ v ≥ 0.1. Therefore, it is suggested that the prominent δ v peak at an altitude of approximately 0.6 km reflected the collision-coalescence growth of falling large raindrops and their subsequent spontaneous breakup. The q v values at altitudes from ∼ 0.7-3.0 km ranged from 5.3-7.3 g kg −1 (Fig. 10c), indicating overall moisture enhancement compared to those values measured at the onset of the moderate warm-front precipitation event (Fig. 9c).

Summary and conclusions
Observations of precipitation and associated precursor clouds were made with two co-located lidars (a 355 nm polarization lidar and water vapor Raman lidar) equipped with waterproof transparent roof windows at the Wuhan University atmospheric observatory (30.5 • N, 114.4 • E; 73 m above sea level). The lidar observations obtained during reachingsurface precipitation events indicate that the rainfall-induced liquid-water accumulation on the roof windows of the li-dars yielded a nearly height-independent lidar signal (rangecorrected signal X) attenuation, whereas neither the X vertical structure nor the magnitude or vertical structure of the volume depolarization ratio (δ v ) were altered. Furthermore, the liquid-water accumulation on the roof windows of the lidars also had nearly no effect on the obtained subcloud profiles of the water vapor mixing ratio measured by the Raman lidar. These observations are consistent with the results of our artificial water-splashing experiment on the roof windows.
Warm-front precipitation events and their precursor cloud evolution were reported in this paper based on two case studies corresponding to light and moderate rainfall occurring at the Earth's surface. The lidar-observed precursor clouds showed a systematic descent for each case. The descending clouds changed gradually from cirrus and altocumulus to altostratus before rainfall occurred, with gradually increasing moisture, and the southwesterly wind prevailed over most altitude ranges of the cloud layers. These features indicate that, in each case, a warm front was approaching our lidar site. The precursor clouds had underlying ice virgae in their later descent phases. When the subcloud virgae reached an altitude slightly below the 0 • C level, rainfall at the surface began. The hours-long precipitation streaks shown in the lidar signal (X) and volume depolarization ratio (δ v ) profiles reveal some ubiquitous features of the microphysical processes of precipitating hydrometeors.
For the light warm-front rain event, since the reachingsurface precipitations and virgae occurred alternately over a short timescale from a few minutes to tens of minutes and since their respective precipitation streaks had nearly the same dark-band structures, both reaching-surface precipitations and virgae originate from the same source cloud (because a warm-front cloud system is generally widespread and slowly varying). Through an analysis combining the lidar profiles of reaching-surface precipitations and virgae, we find that the reaching-surface precipitation began as icephase-dominant hydrometeors fell out of a liquid apparentsource-cloud layer at altitudes above the 0 • C isotherm level. The depolarization ratio magnitude of falling hydrometeors increased from the liquid-water values (δ v < 0.09) to the ice/snow values (δ v > 0.20) during the first 100-200 m of their descent. Subsequently, the falling hydrometeors yielded a dense layer with an ice/snow bright band occurring above and a liquid-water bright band occurring below (separated by a lidar dark band) as a result of crossing the 0 • C level. In the ice/snow bright band, larger particles formed by riming and/or aggregation, because the broad size distributions of the pristine hydrometeors falling out of their apparentsource-cloud base could lead to local accretion. The completion of the melting process of most falling ice particles appeared at altitudes (hundreds of meters) below the 0 • C isotherm level, as indicated by the local depolarization minimum located immediately beneath (∼ 100 m) the observed lidar dark-band minimum. After going through the dark band, most falling raindrops shrunk or vanished in the water bright band due to evaporation, whereas a few large raindrops survived and fell out of the water bright band when the rainfall rate below the liquid apparent-source-cloud base was high enough. Large raindrops might originate from both the complete melting of falling large ice/snow particles and collision-coalescence formation in the dense water bright band. We also find that a prominent depolarization δ v peak (0.10-0.40) always occurred at an altitude of approximately 0.6 km when precipitation reached the surface, reflecting the collision-coalescence growth of large falling raindrops (sparse large raindrops with high fall velocities further grew in size by collecting smaller raindrops along their fall paths) and subsequent spontaneous breakup. The δ v peak observed at an altitude of ∼ 0.6 km provides an indicator in advance (∼ 2 min) of precipitation that reached the surface.
For the moderate warm-front rain event, although the apparent source cloud was undetected, owing to strong attenuation, the lidar-detectable microphysical process (at the altitudes of the ice bright band and below) was similar to that observed in the light-rain case. However, the δ v values in both the ice bright band and water bright band were generally larger than their counterparts in the light-rainfall case, indicating that more large ice/snow particles and raindrops were involved in moderate precipitation. Furthermore, the X and δ v precipitation streaks had complicated vertical structures at altitudes around and below the dark band and showed strong variations on the timescale of minutes. These profile details suggest that large particles sometimes fell out of the ice bright band during moderate precipitation, concealing the lidar dark band produced by the melting effect of most relatively small particles in precipitating hydrometeors. Thus, the lidar dark band became fuzzy. The δ v maxima observed at an altitude of approximately 0.6 km were also larger than those observed in the light warm-front rain case. This suggests larger concentrations of raindrops with spontaneous breakup sizes around this altitude.
has a value range of [0.2− 0.24 R min +0.2 , 0.2). Since δ v,threshold (R) is a slowly varying function of R as seen from Eq. (A7) (particularly when R min ≥ 5), the discrimination criterion of nonspherical particles expressed by δ v (z) can be written approximately as When R min = 7, the discrimination criterion of nonspherical particles is given by δ v (z) > 0.167, which is equivalent to δ p (z) > 0.2 approximately. In conclusion, the particle depolarization ratio δ p has a quasi-linear dependence on the volume depolarization ratio δ v and a very weak dependence on lidar backscatter ratio R (when R ≥ 5). This favorable functional dependence allows us to utilize δ v in discriminating whether the dominant lidar backscattering is attributed to spherical or nonspherical particles in a given backscatter volume. If R min is the minimum of the R value range for interested clouds/virgae (e.g., R min = 7 for the precipitation-related virgae), the discrimination criterion of spherical particles expressed by δ v (z) is given by Eq. (A5), while the discrimination criterion of nonspherical particles expressed by δ v (z) is given approximately by Eq. (A8).
Data availability. Lidar data used to generate the results in this work are available from the corresponding author with permission (E-mail: yf@whu.edu.cn).
Author contributions. YY performed the lidar measurements, made the data analysis and wrote the initial article. FY conceived the project, led the study and finalized the article. FL, YZ and CY built the lidar systems for precipitation observations. YH participated in scientific discussions and suggested analysis. All authors discussed the results and commented on the article.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.