MARCH 2016
CHEN ET AL.
621
Raindrop Size Distribution in a Midlatitude Continental Squall Line Measured by Thies Optical Disdrometers over East China BAOJUN CHEN Key Laboratory of Mesoscale Severe Weather/MOE, and School of Atmospheric Sciences, Nanjing University, Nanjing, China
JUN WANG AND DIANLI GONG Weather Modification Office of Shandong Province, Jinan, China (Manuscript received 8 May 2015, in final form 6 October 2015) ABSTRACT Disdrometer data measured by ground-based optical disdrometers during a midlatitude continental squall line event on 18 August 2012 in Shandong Province, eastern China, are analyzed to study characteristics of raindrop size distribution (DSD). Four disdrometers simultaneously performed continuous measurements during the passage of the convective line. The convective line was partitioned into three regions: the convective center, leading edge, and trailing edge. Results show distinct differences in DSDs and integral rainfall parameters between the convective-center and the edge regions. The convective center has higher drop concentrations, larger mean diameters, and wider size distributions when compared with the edge regions. The leading and trailing edges have similar drop concentrations, but the latter has larger mean diameters and wider size distributions. The shape of DSD for the convective center is convex down, whereas it is convex upward in tropical continental squall lines, as reported in the literature. There is also spatial variability of the DSD and its integral rainfall parameters in the along-convective-line direction.
1. Introduction Squall lines are linear or quasi-linear mesoscale convective systems that can produce heavy rains, hail, strong winds, and possibly tornadoes. According to the Glossary of Meteorology of the American Meteorological Society (Glickman 2000), the squall line is defined as ‘‘a line of active thunderstorms, either continuous or with breaks, including contiguous precipitation areas resulting from the existence of the thunderstorms.’’ Squall lines have been frequently observed in both the tropics and midlatitudes. There exist three distinct modes of squall lines (Parker and Johnson 2000): trailing stratiform (TS), leading stratiform (LS), and parallel stratiform (PS). Observations indicate that about 60% of midlatitude
Denotes Open Access content.
Corresponding author address: Baojun Chen, School of Atmospheric Sciences, Nanjing University, Xianlin Campus, 163 Xianlin Avenue, Nanjing 210023, China. E-mail:
[email protected] DOI: 10.1175/JAMC-D-15-0127.1 Ó 2016 American Meteorological Society
squall lines are of the TS type, while the remaining 40% are of the LS and PS types (Parker and Johnson 2000; Meng et al. 2013). The squall line is an ideal precipitation system for studying microphysical properties of different precipitation types, since it combines two distinctly different precipitation types (convective and stratiform) in one single event. As the most fundamental aspect of precipitation microphysics, knowledge of the raindrop size distribution (DSD) is essential for understanding precipitation processes, estimating rainfall, and improving microphysics parameterizations in numerical cloud models. There have been several studies of DSD measurements in squall lines. Maki et al. (2001) analyzed the characteristics of the DSD in 15 tropical continental squall lines observed in Darwin, Australia, with the Joss–Waldvogel (JW) impact disdrometer. They showed distinct differences in DSD patterns between the convective-center and stratiform regions of the squall line. Nzeukou et al. (2004) and Moumouni et al. (2008) also studied the DSD characteristics of tropical squall lines using disdrometer data collected at Dakar (Senegal) and Benin in West Africa with the JW and
622
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
optical disdrometers, respectively. Jung et al. (2012) reported microphysical properties of DSD in a subtropical maritime squall line in Taiwan, using Palomar Observatory Sky Survey (POSS) observation data. Relatively few published studies have focused on the DSD of midlatitude continental squall lines (Uijlenhoet et al. 2003). Although squall lines have been frequently observed in midlatitude eastern China, and their general features have also been investigated (Meng et al. 2013), not much is known about the DSD characteristics in squall lines occurring in this region because of the lack of disdrometer measurements. Although studies mentioned above have addressed the DSD of squall lines, they mainly focused on the difference in two rain types (convective and stratiform) or three precipitation regions (convective line, stratiform, and transition). Many recent studies have shown the small-scale spatial variability of DSDs and their integral rainfall parameters, especially for convective rainfall (e.g., Tokay and Bashor 2010; Jaffrain and Berne 2012; Jameson et al. 2015a). Since a squall line is an organized line of convective storms, is there any spatial difference in DSD along the convective line direction? In this paper, we report the characteristics of the DSD in a midlatitude continental PS squall line observed in eastern China. The data observed during the passage of convective lines with four Thies optical disdrometers are analyzed to study the characteristics of the DSD and its spatial variability.
2. Data and methods a. Thies disdrometer The DSD data used in this study were collected with the Thies disdrometer during a squall line that passed eastward over the central part of Shandong Province, eastern China, on 18 August 2012. During the overpass of the squall line, four Thies disdrometers simultaneously performed continuous measurements at a time resolution of 60 s. Figure 1 shows a horizontal radar image of the squall line at 1.58 elevation angle observed by the China New Generation Weather Radar (CINRAD)-SA Doppler weather radar located at Jinan (36.808N, 116.788E). The disdrometer observation sites are indicated in the figure: Zhangqiu (ZHQ), Baiquan (BQ), Duozhuang (DZH), and Sanchacun (SCC). The radar image corresponds to a time around 1842 LST, when the storm had developed into the mature stage, showing a typical PS squall line structure as defined in Parker and Johnson (2000). Since all of four
VOLUME 55
FIG. 1. Horizontal radar image at 1.58 elevation angle of the 18 Aug 2012 squall line at 1842 LST, observed by the CINRAD-SA Doppler weather radar located at Jinan (36.808N, 116.788E). Four Thies disdrometers were collocated at ZHQ, BQ, DZH, and SCC, respectively.
disdrometers were located in the convective precipitation region of the squall line system, we have not studied the characteristics of the DSD from the stratiform precipitation region because of the lack of disdrometer data. Thies is a laser-optical disdrometer, which can simultaneously count and measure the size and fall velocity of hydrometeors (Bloemink and Lanzinger 2005). Thies instruments were assessed in detail by Frasson et al. (2011) and Sarkar et al. (2015) and have been used in many other studies (e.g., Brawn and Upton 2008; Fernández-Raga et al. 2010; Jameson et al. 2015a,b). The core element of the instrument is an optical sensor that generates a parallel horizontal light beam of 0.75-mm thickness with a measuring area of 45.6 cm2 (228 mm long, 20 mm wide). When the precipitation particles fall through this beam, the receiving signal is reduced. The amplitude of the reduction is related to the size of the particles, and the duration of the reduction is related to the fall speed. Particle sizes are subdivided into 22 classes, and the measurement range is 0.2–8 mm in diameter. The size class width is varied from 0.125 to 0.5 mm. Fall speeds are subdivided into 20 classes, and the measurement range is 0.2– 10 m s21. The velocity class width is also varied from 0.2 to 1.0 m s21. For a detailed description of the instrument, please refer to the instrument’s manual (Thies Clima 2007).
MARCH 2016
623
CHEN ET AL.
FIG. 2. Time series of the 1-min rain rate derived from the disdrometer at ZHQ and measured by a collocated tipping-bucket rain gauge from 1830 to 1910 LST. The green line indicates the tipping-bucket rain gauge’s measurements, while blue and red lines indicate the rain rate computed from the quality-controlled and non-quality-controlled disdrometer data, respectively.
The measurement accuracy of optical-like disdrometers might be affected by strong winds and heavy rainfall, resulting in spurious particles with unrealistic fall speeds and diameters to be observed in the data (Friedrich et al. 2013). In addition, if particles are not totally within the light beam, they may be registered as small particles that fall faster than other particles observed at that size. These particles are called as ‘‘margin fallers’’ (Yuter et al. 2006). To minimize the measurement errors caused by strong winds, splashing, or margin fallers, particles outside the 660% of the fall speed– diameter relationship for rain (Atlas et al. 1973) are removed from the observed data following Friedrich et al. (2013). The data containing fewer than 10 drops or with a rain rate less than 0.1 mm h21 have been also disregarded as noise (Tokay and Bashor 2010). Comparing the disdrometer-derived rain rate (Figs. 2, 3), one can see that the quality-controlled data gave a closer agreement with rain gauge measurements, especially during heavy rainfall. Non-quality-controlled data tended to overestimate significantly the rain rate; the same as reported by Lanzinger et al. (2006). The quality control substantially decreased bias between the disdrometer and rain gauge, indicating that the qualitycontrol scheme is reliable. While the bias after quality control is likely caused by miscalculation of the disdrometer’s sensing area as pointed out by Frasson et al. (2011).
FIG. 3. Scatterplots of the 1-min rain rate calculated by the disdrometer at ZHQ and measured by the collocated tipping-bucket rain gauge. Rain rates derived from the quality-controlled and nonquality-controlled disdrometer data are indicated by blue and red circles, respectively.
b. Raindrop size distribution and rainfall integral parameters The number concentration of raindrops per unit volume per unit size interval at the discrete time interval has been calculated from the Thies disdrometer counts using the following equation: 20
N(Di ) 5
n
å ADtVijDD ,
j51
j
(1)
i
where N(Di) is the drop size distribution in size class i (m23 mm21), Di is the mean diameter of class i (mm), DDi is the width of size class i (mm), nij is the number of drops within the size class i and velocity class j, A is the measuring area (m2), Dt is the time interval in seconds, and Vj is the fall speed of velocity class j (m s21). Measured DSDs can be used to calculate rainfall microphysical parameters. Here we selected four integral parameters to characterize the observed DSD, including the total number concentration NT (m23), rainwater content W (g m23), rain rate R (mm h21), and radar reflectivity factor Z (mm6 m23), which were computed as follows: 22
NT 5 W5
20
n
ij , å å ADtV
i51 j51
(2)
j
22 p 3 1023 rw å 6 i51
20
n
ij , å D3i ADtV
j51
j
(3)
624
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
VOLUME 55
TABLE 1. Summary of the disdrometer observations. Times are LST. The mean and peak rain rate (mm h21) for each region are shown in parentheses. Site
No. of 1-min spectra
Convective center
Leading edge
Trailing edge
ZHQ BQ DZH SCC
77 77 72 66
1837–1854 (70.4, 191.4) 1835–1853 (68.2, 154.5) 1829–1909 (47.4, 95.1) 1833–1853 (61.8, 95.1)
1823–1836 (2.5, 7.4) 1822–1834 (1.8, 9.1) 1813–1828 (5.6, 12.6) 1812–1832 (1.7, 6.7)
1855–1940 (3.6, 14.7) 1854–1938 (4.7, 18.0) 1910–1925 (5.1, 17.8) 1854–1917 (1.3, 9.7)
22
R 5 6p 3 1024 å
20
n
ij , å D3i ADt
and
(4)
i51 j51
22
Z5
20
n
ij , å å D6i ADtV
i51 j51
(5)
j
where rw is the density of water (1 g cm23). For the convenience of presentation, we also used the decibel units of radar reflectivity, dBZ (510 log10Z), instead of the radar reflectivity factor. The gamma size distribution proposed by Ulbrich (1983) is used to parameterize the observed DSD: N(D) 5 N0 Dm exp(2LD) ,
(6)
where D is the raindrop diameter (mm), N0 is the intercept parameter (mm212m m23), m is the shape parameter, and L is the slope parameter (mm21). The three parameters were estimated from the truncated-moment method using the second, fourth, and sixth moments of the observed spectra, respectively. Detailed descriptions on the truncated-moment method have been addressed by Ulbrich and Atlas (1998) and Vivekanandan et al. (2004) and are not presented herein. In addition, the mass-weighted mean diameter Dm (mm) and the generalized intercept parameter NW (mm21 m23) were used to illustrate the observed DSDs. The NW can be derived as follows (Bringi et al. 2003): 44 NW 5 prw
! 103 W . D4m
(7)
c. Partitioning of convective lines According to Maki et al. (2001), DSDs in tropical continental squall lines show distinct differences between the convective-center and the edge regions of convective lines. The convective center was defined by Maki et al. (2001) as the strongest part of the convective line, which had a rain rate larger than 20 mm h21. The leading and the trailing edges were defined as the front and rear sides of the convective center, respectively. As a primary purpose of this study is to characterize the DSDs of different regions for the convective line in the observed continental squall line and to further compare
the DSDs with those of tropical continental squall lines reported in Maki et al. (2001), the method for partitioning convective lines proposed by Maki et al. (2001) is adopted in this work. It can be seen from Table 1 that the three regions and the four observation locations were diverse in terms of duration and intensity of rainfall. Convective centers have higher rain rates than the edge regions. Precipitation originating in the convective center was observed 21%–58% of the time (34% for four-site averaging), but 84%–95% of the total rainfall (90% for four-site averaging) was provided by convective centers. One should keep in mind that the current study used 20 mm h21 as the threshold for separating convective centers. Sensitivity tests showed that the changing threshold [e.g., 10 mm h21 as in Testud et al. (2001)] had an impact on the DSD and integral rainfall parameters. However, there are insignificant differences between the two results because of a few samples observed in the range of 10–20 mm h21 (three for ZHQ, six for BQ, two for DZH, and zero for SCC).
3. Results and discussion a. Time series of raindrop size distribution Figure 4 shows the time evolution of the drop size distribution N(D) during the passage of the squall line, as derived from the 1-min disdrometer observations at the four sites, respectively. Rain rate during the convective center phase was also plotted to separate the three phases of the convective line. All four DSDs exhibit a similar behavior, namely, low number concentrations and narrow size spectra during the leading and trailing edges where the maximum drop diameter was less than 4.75 mm, and high concentrations and wide size spectra during the convective center phase in which the maximum drop diameter was more than 6.25 mm. This figure clearly demonstrates that the convective center tends to produce wide spectra, whereas the leading and trailing edges exhibit narrow spectra. We note from Fig. 4 that nearly all of the raindrops with concentration exceeding 1000 m23 mm21 were concentrated entirely in the size range of D , 1 mm, and they contributed to about 90% of the total number concentration. Further, it is found that the raindrops
MARCH 2016
CHEN ET AL.
FIG. 4. Temporal evolution of the raindrop size distribution N(D) (m23 mm21) measured at (a) ZHQ, (b) BQ, (c) DZH, and (d) SCC during the passage of the squall line. The scale for the drop concentration is given in (d). The red lines indicate the rainfall rate during the convective center.
with concentration larger than 10 000 m23 mm21 were below 0.6 mm in diameter, occurring in the convective center phase and contributing to more than 95% of the total concentration in this period. It is also noted that the highest concentration of drops measured by disdrometers at all four sites was in the first two size classes, which correspond to an average diameter of 0.188 and 0.312 mm. Such a large number of small drops, however, contributed to over 50% of the total concentration, though they contributed only 0.7% of the rain rate. Comparing Figs. 4a and 4b, there were similar structure in the time series of DSD between ZHQ and BQ, where their separation distance was about 4 km. The maximum raindrop diameter in the spectra at the two sites was 8 mm. The overall coefficient of drop variation (CVD) for the entire time period, here defined as the standard deviation of the maximum diameter divided by the mean value of the maximum diameter to measure the fluctuating magnitude of DSD’s width, was 0.40 and 0.42 for ZHQ and BQ, respectively. When considering
625
the convective center phase, the derived value of CVD from this period was 0.23 at both sites. The differences between the two time series were found in the number concentration. The highest total number concentration was 30 978 m23 in ZHQ, occurring at 1846 LST, and 26 690 m23 in BQ at 1844 LST. The time series of DSD during the convective line’s passage in DZH showed a substantial difference from the previous two sites (Fig. 4c). The event’s highest total number concentration, 21 615 m23, was observed at 1831 LST. The largest drop at DZH was 6.25 mm in diameter, and CVD was 0.25 for the entire period and 0.11 for the convective center phase, both smaller than those in the previous two sites. During the convective line’s passage in SCC (Fig. 4d), the highest total number concentration was 21 795 m23 occurring at 1849 LST, and the maximum drop diameter was 8 mm. The values of CVD were 0.47 and 0.12, respectively, for the entire period and the convective center phase. The results showed temporal and spatial variability of the DSD in convective lines. It can be seen from Fig. 4 that there is an abrupt change in DSDs during the transition period from the leading edge to the convective center, and then to the trailing edge. To illustrate this variation, we plotted the size distributions for the last minute of the leading edge (LE_E), the first and last minutes of the convective center (CC_S and CC_E, respectively), and the first minute of the trailing edge (TE_S) at each observation site (Fig. 5). The standard deviation of the mass spectrum sM 5 Dm/(4 1 m)1/2, as defined by Ulbrich and Atlas (1998), was also presented to quantify the spectrum width. The CC_S had a wider spectrum breadth than LE_E, but there was no significant difference between CC_E and TE_S. The drop concentration in the CC_S (CC_E) was higher than in LE_E (TE_S) at all size ranges. This shows that there is a sudden increase of the drop size and concentration from the leading edge to convective center and a decrease of the drop size and concentration from the convective center to the trailing edge, thereby resulting in a dramatic increase or decrease in rain rates. Furthermore, the difference in DSDs between LE_E and CC_S was noticeably larger than that between CC_E and TE_S.
b. Averaged raindrop spectra The averaged raindrop spectra and their integral rainfall parameters for the leading edge (LE), convective center (CC), and trailing edge (TE) of the convective line at each observation site are shown in Fig. 6 and Table 2, respectively. Here we used the empirical fall speed–diameter relationship of Atlas et al. (1973) to compute rain rate from the average spectra. In addition, the gamma distribution fits on average spectra is also
626
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
VOLUME 55
FIG. 5. The 1-min drop size distributions for the end of the leading edge (LE_E), the start and end of the convective center (CC_S, CC_E), and the start of the trailing edge (TE_S), at four observation sites, as indicated. Values in parentheses show the standard deviation of the mass spectrum, as defined by Ulbrich and Atlas (1998).
plotted in the figure, and the corresponding parameter values are listed in Table 3. As shown in Fig. 6 and Table 3, the DSDs in the convective center at all four sites are convex down, with the shape parameter m in the range from 21.631 to 21.272, whereas the DSDs in the leading edge and trailing edge have convex shapes downward or upward, with the parameter m in the ranges from 20.319 to 1.229 and from 21.705 to 0.170 for the leading edge and trailing edge, respectively. The maximum raindrop diameters exhibit a similar tendency at four sites, all in the order of CC . TE . LE. The slope parameter L for the convective center is in the 0.689–1.001 range, while for the leading edge and the trailing edge they are in the 2.148–3.462 and 1.018–2.195 ranges, respectively. The slope parameters at all sites are in the order of CC , TE , LE. Overall, Table 3 shows that among the three
precipitation regions the convective center had the highest N0, the lowest m, and the smallest L, while all these three parameters had higher values for the leading edge than for the trailing edge. Compared to the leading edge or trailing edge for each observation site, the convective center had higher drop concentrations at all size regimes. This is the reason why the convective center has the highest total number concentration. Meanwhile, it is also one of the reasons why the convective center has the highest rain rate and rainwater content and the highest reflectivity. It is found from Table 2 that the trailing edge at ZHQ and BQ had larger values for all integral rainfall parameters than the leading edge, since the DSD of the trailing edge had higher concentration at all size ranges (Figs. 6a,b). However, the distribution patterns of the DSD for the two edge regions differed in some size intervals at
MARCH 2016
CHEN ET AL.
627
FIG. 6. Average raindrop size spectra for the LE, CC, and TE at the (a) ZHQ, (b) BQ, (c) DZH, and (d) SCC observation sites. Fitting curves of gamma distribution functions using the truncated-moment method are also shown in solid lines.
the other two sites. Comparing the trailing edge with the leading edge for DZH (Fig. 6c), we found a decrease in the number of small and intermediate drops (D , 3 mm) and an increase in the number of large drops (D . 3 mm). This characteristic DSD is the reason why the trailing edge has a lower concentration, less rainwater content and rainfall rate, and a higher reflectivity than the leading edge. For SCC (Fig. 6d), we noted from comparing the trailing edge with the leading edge 1) an increase in the number of small drops (D , 1 mm), 2) a decrease in the number of midsize drops (1 , D , 3 mm), and 3) an increase in the number of large drops (D . 3 mm). Such DSD behavior may explain why there are no significant differences in the integral rainfall parameters at this site between the leading edge and the trailing edge.
The small drops dominated the number concentration. However, they contributed less than 10% of the rain rate on average. The rainfall intensity is mainly contributed by the midsize and large drops, accounting for 55%–70% and 20%–35% of the total rainfall over the four sites, respectively. The contributions to the rain rate from small, midsize, and large drops are different in three precipitation regions. For the leading edge the small and midsize drops contributed 11%– 28% and 70%–85% of the rain rate, respectively, while the contribution of the large drops was very insignificant. For the trailing edge the contributions from the small, midsize, and large drops was 8%–13%, 47%– 74%, and 18%–44%, respectively. Very large drops (D . 5 mm) contributed negligibly to the rain rate since they were occasionally observed during the trailing
628
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
TABLE 2. Rainfall integral parameters derived from the averaged drop size spectra for the LE, CC, and TE. Terms NT, W, R, and dBZ represent total drop concentration (m23), rainwater content (g m23), rain rate (mm h21), and radar reflectivity, respectively. Mean mass diameter Dm (mm) and maximum diameter Dmax (mm) are also given for each case. Site ZHQ
BQ
DZH
SCC
All
LE CC TE LE CC TE LE CC TE LE CC TE LE CC TE
NT
W
R
dBZ
Dm
Dmax
387.4 13489.5 546.8 329.9 11309.0 812.2 1245.7 12751.2 823.9 344.1 11187.8 462.5 576.1 12264.6 659.6
0.146 3.705 0.193 0.097 3.619 0.309 0.347 2.555 0.253 0.098 3.407 0.087 0.171 3.143 0.224
2.69 90.7 4.12 1.90 95.2 7.20 7.32 60.7 6.09 1.91 86.8 1.82 3.42 78.2 5.07
31.4 55.1 36.2 31.1 56.6 40.3 38.0 51.5 40.8 31.7 55.0 35.0 34.2 54.3 38.6
1.47 2.76 1.90 1.62 3.23 2.23 1.83 2.50 2.44 1.64 2.95 2.03 1.70 2.82 2.14
3.25 8.00 4.75 3.75 8.00 4.75 4.75 6.25 5.25 4.25 8.00 4.75 4.75 8.00 5.25
edge. For the convective center the contributions of the small, midsize, and large drops to the rain rate were 4%–6%, 35%–49%, and 45%–59%, respectively, while the contribution from very large drops achieved 3%–12%. It is noted from Table 2 that there are differences in the integral rainfall parameters derived from average size spectra between the four different sites. For example, the site-to-site relative difference of the total drop concentration varies from 4.2% to 175.5% and from 1.4% to 61% for the leading and trailing edges, respectively, while it varies from 1.1% to 18.8% for the convective center. Larger variability is found in the rain rate, ranging from 0.5% to 179.6% and from 16.8% to 185.1% for the leading and trailing edges, respectively, and from 4.4% to 138.9% for the convective center. To further examine the difference in the DSDs between different observation sites, the averaged raindrop spectra for LE, CC, and TE regions at four observation sites are compared in Fig. 7. The standard deviation of the mass spectrum or the spectrum width was also presented. As shown in Table 3 and Fig. 7, both the LE and TE regions show substantial differences in their distribution shapes (m) and spectrum widths (sM) of the average DSDs between different observation sites (Figs. 7a,c), for the CC regions the spectrum widths are distinctly different and the distribution shapes are similar (Fig. 7b). Ulbrich and Atlas (1998) used the normalized standard deviation of the mass spectrum sM/DM to measure the spectral shape. The normalized standard deviation of the average DSD in Fig. 7 was 0.43–0.52 for LE, 0.49–0.66 for TE, and 0.61–0.65 for CC. This
VOLUME 55
TABLE 3. Gamma distribution parameters derived from the average spectra using the truncated-moment method. Site ZHQ
BQ
DZH
SCC
All
LE CC TE LE CC TE LE CC TE LE CC TE LE CC TE
N0 (mm212m m23)
m
L (mm21)
5505.9 3996.0 1349.0 1273.5 2520.2 727.6 2256.6 3836.9 379.2 950.9 3689.6 138.2 2694.4 3606.1 734.6
1.229 21.616 0.170 0.285 21.631 20.548 20.319 21.455 21.099 20.112 21.272 21.705 0.704 21.456 20.432
3.462 0.845 2.195 2.579 0.689 1.451 2.148 1.001 1.095 2.317 0.939 1.018 2.795 0.908 1.647
further illustrates the difference (similarity) in the distribution shapes in the edge (convective center) regions between different observation sites.
c. Gamma distribution parameters for the convective center regions Here only the convective center regions have been studied since the convective center dominated the precipitation of the convective line. Since the units of the intercept parameter N0 depend on the value of the shape parameter m, this makes it very inconvenient to discuss the variability of N0. The generalized intercept parameter NW is independent of the shape of the DSD, which means that whatever the shape of an observed DSD is, the corresponding NW is the intercept parameter of the equivalent exponential DSD [i.e., Eq. (6) with m 5 0] with the same W and Dm. Hence, we used NW instead of N0. The temporal evolution of the parameters NW, m, and L during the convective center phases is shown in Fig. 8. All three parameters exhibited an intraphase variation during the passage of the convective center over the four observation sites, with the coefficient of variation from 0.47 to 0.55 for NW, from 0.13 to 0.37 for m, and from 0.31 to 0.41 for L. The NW had a larger variation as compared to the other two gamma parameters. For the whole convective-center dataset the NW was in the 900–14 700 mm21 m23 range, while the m and L ranges were from 22.379 to 20.322 and from 0.272 to 1.519, respectively. Note that all m of convective centers in this study are negative, while they are mostly positive in tropical continental squall lines observed in Darwin (Maki et al. 2001). The disdrometer data used in Maki et al. (2001) were collected with the JW disdrometer, which allowed the measurement of drops ranging from 0.3 to 5.5 mm in diameter (Joss and Waldvogel 1967). Tokay et al.
MARCH 2016
CHEN ET AL.
629
FIG. 8. Temporal evolution of the generalized intercept parameter NW (mm21 m23), slope parameter L (mm21), and shape parameter m during the convective center phases.
FIG. 7. Average raindrop size spectra for (a) LE, (b) CC, and (c) TE regions of convective lines. The standard deviation of the mass spectrum is shown in parentheses.
(2001) found that the JW disdrometer severely underestimated the number of small drops (D , 1.5 mm), thus producing a pronounced effect on the shape parameter of gamma-fitted DSDs. The size ranges of the Thies disdrometer are 0.2–8.0 mm in diameter. Large concentrations of small drops that were measured by Thies are likely the main reason why the DSD in the convective center tends to be convex down and to have negative values of m. To illustrate this, the concentration of small drops for diameters less than 1.5 mm in the average size spectra for the convective center region at ZHQ was decreased by 50% and 90%, respectively. The resulting value of m was 20.937 and 0.861, respectively, indicating that the shape parameter increases with decreasing small-drop concentrations. Since the Thies measured very large drops that were beyond the JW disdrometer maximum size limit, the effect of very large drops on the shape parameter for the convective center was also examined by removing raindrops with diameters larger than 5.5 mm. The corresponding value of m was 21.893, which was very close to the value of 21.616 calculated from the untruncated spectra. Therefore, the differences in the DSD patterns for the convective center between this study and Maki et al. (2001) are mainly due to the concentration of small drops measured by the Thies and JW disdrometers.
630
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
VOLUME 55
For each observation site, the time series of m and L showed a rather similar behavior and a very strong correlation. Brandes et al. (2003) and Zhang et al. (2003) proposed an empirical polynomial m–L relation for convective precipitation with rain rates greater than 5 mm h21 as follows: L 5 0:0365m2 1 0:735m 1 1:935.
(8)
This m–L relation was obtained from disdrometer observations in Florida during the summer of 1998, and it should be tuned using local disdrometer data collected at various geographical locations or climatological regimes and in various seasons (Zhang et al. 2003; Vivekanandan et al. 2004). Recent studies have demonstrated the variability in the m–L relation at different locations (e.g., Brawn and Upton 2008; Chu and Su 2008; Chang et al. 2009; Kumar et al. 2011; Chen et al. 2013). The characteristics of m and L in squall lines have been studied as well (e.g., Maki et al. 2001; Uijlenhoet et al. 2003); nevertheless, not much is known about the m–L relation of squall lines. The scatterplots of m–L values for convective-center regions at the four sites are shown in Fig. 9. The corresponding m–L relations derived here using a polynomial least squares fit are shown as solid lines. The specific equations are given as follows: L(ZHQ) 5 0:0860m2 1 0:845m 1 1:959,
(9)
L(BQ) 5 0:0460m 1 0:704m 1 1:697,
(10)
2
L(DZH) 5 0:0209m 1 0:838m 1 2:169, and 2
(11)
L(SCC) 5 0:1462m 1 1:333m 1 2:524,
(12)
2
where acronyms in parentheses indicate the observation sites. All four fitted curves in Fig. 9 show a similar behavior, namely, the smaller the value of L (higher rain rates), the smaller the negative value of m. This suggests that the DSD tends to be more convex down with increasing rainfall intensity and to have a lower fraction of midsize drops and a higher fraction of small and large drops, reflecting less evaporation of small drops or more active drop breakup and coalescence mechanisms at higher rain rates. The fitted m–L relations exhibited a rather large difference among the four different sites and therefore have different m values given the same L value. This suggests that the m–L relation obtained from single-point observations cannot reliably represent the convective line. The m–L relations obtained in this work differ from that in Florida given in Eq. (8). The differences are likely attributed in part to the natural variability of rainfall and in part to dataset processing, because the threshold values of rain rate to determine convective spectra are
FIG. 9. Scatterplots of the m–L values for convective center phases. The solid lines indicate least squares fits of the polynomial m–L relationship.
different. To eliminate the errors from different criteria, the thresholds for drop counts .1000 and rain rate .5 mm h21 used in Brandes et al. (2003) and Zhang et al. (2003) were adopted to separate convective spectra. The revised plot is shown in Fig. 10 and the corresponding m–L relation is given in Eq. (13). Here the disdrometer data with the above-mentioned threshold collected at all four sites were put together to obtain an overall reliable formula. The figure contains 108 data points and their m values are still negative: L 5 0:0585m2 1 0:812m 1 1:934.
(13)
The revised m–L relation agrees well with Florida’s, and the difference is found in the squared term and linear term. Comparing Eqs. (8) and (13), we see that the coefficients of the squared term and linear term in Eq. (8) are smaller than those in Eq. (13); hence, for a given m value Eq. (8) gives a higher L value than Eq. (13), as shown in Fig. 10.
d. Distribution of Dm and NW The mass-weighted mean diameter Dm and the generalized intercept parameter NW can be used to characterize an overall feature of the DSD. Observation studies have shown that both the NW and Dm vary with climatic regimes, rain types, and rainfall intensity (e.g., Bringi et al. 2003; Ulbrich and Atlas 2007; Sharma et al. 2009; Marzano et al. 2010; Thurai et al. 2010; Chen et al. 2013). To examine the dependence of these two
MARCH 2016
CHEN ET AL.
631
FIG. 10. Scatterplots of the m–L values for DSDs with rain rate .5 mm h21 and total drop counts .1000. The red line represents the fitting of the m–L relationship. Equation (8) is overlaid for comparison.
parameters on rainfall intensity, Fig. 11 shows log10NW and Dm versus rain rate R for the four observation sites. In addition, the fitted power-law relationships using a least squares method are also provided in the figure to obtain a quantitative description of the two parameters with respect to rain rate. Both the NW and Dm values tend to increase with increasing rain rates, and the exponents in the NW–R and Dm–R relations are positive. Note that both the NW–R and Dm–R relations exhibit slight differences between different sites. This indicates that NW and Dm are likely different at various locations even though they are measured within the same one precipitation system and at the same rain rate. Previous studies (e.g., Bringi et al. 2003; Ulbrich and Atlas 2007) showed that the behavior of the DSD parameters NW and Dm was distinctly different for the maritime and continental convective storms. For maritime-like spectra, log10(hNWi) ’ 4–4.5 with hDmi ’ 1.5–1.75 mm (angle brackets denote averages), that is, a higher concentration of smaller-sized drops. For continental-like spectra, log10(hNWi) ’ 3–3.5 with hDmi ’ 2–2.75 mm, that is, a higher concentration of larger-sized drops. Figure 12 shows the average Dm and log10NW for convective spectra in the squall line. It should be noted that the separation of convective spectra in this section is different from the abovementioned method. Herein we adopted the classification scheme of Bringi et al. (2003) to compare our data
FIG. 11. (a) Term log10NW vs R and (b) Dm vs R for the four observation sites. The fitted power-law relationships using a least squares method are also provided in each panel.
with those of their study. The scheme is based on the standard deviation of rain rate over five consecutive 2-min samples. If the standard deviation is .1.5 mm h21 and rain rate is $5.0 mm h21, then it is classified as convective. The classification procedure used here is similar to Bringi et al. (2003), except for the standard deviation over 11 consecutive 1-min samples. It can be seen from Fig. 12 that all of the data points are centered near the continental convective clusters reported by Bringi et al. (2003). The log10NW ranges from 3.43 to 3.68
632
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
VOLUME 55
4. Summary and conclusions
FIG. 12. The average value of log10NW vs average Dm for various locations for convective rains based on the classification of Bringi et al. (2003). The two outlined rectangles correspond to the maritime and continental convective clusters reported by Bringi et al. (2003) and the dotted line represents their stratiform cases.
and Dm ranges from 2.25 to 2.9 mm. The four-site average is 3.55 for log10NW and 2.53 mm for Dm. The NW– Dm pair matches the continental-like cluster well, which further demonstrates the continental characteristics of the squall line in the present study. The NW–Dm pair for convective rain in Darwin matches the maritime-like cluster, as reported by Bringi et al. (2003) and implied in Thurai et al. (2010), which is distinctly different from ours. This is likely associated with local atmospheric aerosols and/or moisture. The observation study by May et al. (2011) showed that tropical thunderstorms in Darwin have a lower NW and larger median volume diameter D0 (which is closely related to Dm for gamma DSD) under the high-aerosol conditions when compared with low-aerosol conditions. The sounding in Zhangqiu at 0000 UTC on 18 August 2012 shows that the squall line forms in a warm, moist atmospheric conditions with a warm layer of ;5 km MSL and a water vapor mixing ratio of ;19 g kg21. The thermodynamical conditions are similar to that in Darwin. However, in eastern China, aerosol concentrations are markedly high (Streets et al. 2008). With increased aerosol loading, the warm-rain processes in deep convective clouds can be suppressed, resulting in lower raindrop concentrations; meanwhile, the ice-phase processes can be enhanced, thus producing more graupel/ hail and resulting in larger raindrops via the melting. Further research is needed to validate that hypothesis.
Raindrop size distributions representative of a parallel stratiform squall line that occurred on 18 August 2012 in Shandong Province in eastern China were analyzed. Four Thies disdrometers with different separation distances simultaneously performed continuous measurements during the passage of the convective line. The convective line was further partitioned into the convective center, leading edge, and trailing edge using a threshold rain rate of 20 mm h21 following Maki et al. (2001). The continental characteristic of the squall line was identified well in the NW–Dm space following Bringi et al. (2003). There are distinct differences in DSDs between the convective-center and edge regions of the convective line. The convective center has higher drop concentrations, larger mean diameters, and wider size distributions as compared to the leading and trailing edges. The leading and trailing edges have similar drop concentrations, but the latter has larger mean diameters and wider size distributions. The shape of DSD for the convective center is convex down, but it is not monodisperse for the leading and trailing edges. The DSD shape of the convective center distinctly differs for this study and Maki et al. (2001) in which it is convex upward. There are also differences in DSDs and the integral rainfall parameters in the convective line direction. The leading and trailing edges show substantial differences in their distribution shapes and spectrum widths of the average DSDs between different observation sites, for the convective center the spectrum widths are distinctly different and the distribution shapes are similar. The m–L relations are varied in different observation sites. The results reveal distinct characteristics of the DSD in continental squall lines between midlatitude eastern China and tropical regions, particularly the DSD shape for the convective center. The convective-center DSD observed in this study is characterized by a downward convex shape with an extremely high concentration of small drops in the presence of very large drops (.6 mm in diameter), likely reflecting the microphysical mechanisms of rain formation in the squall line through the melting of graupel and hailstones followed by drop coalescence and breakup. In future work we plan to use cloud-resolving models with two-moment bulk microphysics scheme to investigate the impact of gamma shape parameters on the structure and precipitation of the squall line to improve model representative. Cloud modeling with spectral bin microphysics would be valuable for understanding aerosol effects on DSDs in this region. The findings of this
MARCH 2016
CHEN ET AL.
paper were limited to a case study of PS squall lines in eastern China. Different types of squall lines (i.e., TS and LS types) should be also considered in future work. Acknowledgments. The authors thank Dr. Remko Uijlenhoet and two anonymous reviewers for their constructive comments, which have helped to improve the quality of the paper. This work was supported by the National Natural Science Foundation of China (41175118, 41275044) and the National Basic Research Program of China (2013CB430105). Author JW acknowledges funding from Shandong Meteorological Bureau (2012SDQX12). REFERENCES Atlas, D., R. C. Srivastava, and R. S. Sekhon, 1973: Doppler radar characteristics of precipitation at vertical incidence. Rev. Geophys. Space Phys., 11, 1–35, doi:10.1029/RG011i001p00001. Bloemink, H. I., and E. Lanzinger, 2005: Precipitation type from the Thies disdrometer. WMO Technical Conf. on Meteorological and Environmental Instruments and Methods of Observation 2005, Bucharest, Romania, WMO, 7 pp. [Available online at https://www.wmo.int/pages/prog/www/IMOP/publications/ IOM-82-TECO_2005/Papers/3(11)_Netherlands_4_Bloemink.pdf.] Brandes, E. A., G. Zhang, and J. Vivekanandan, 2003: An evaluation of a drop distribution–based rainfall estimator. J. Appl. Meteor., 42, 652–660, doi:10.1175/1520-0450(2003)042,0652: AEOADD.2.0.CO;2. Brawn, D., and G. Upton, 2008: On the measurement of atmospheric gamma drop-size distributions. Atmos. Sci. Lett., 9, 245–247, doi:10.1002/asl.198. Bringi, V. N., V. Chandrasekar, J. Hubbert, E. Gorgucci, W. L. Randeu, and M. Schoenhuber, 2003: Raindrop size distribution in different climatic regimes from disdrometer and dualpolarized radar analysis. J. Atmos. Sci., 60, 354–365, doi:10.1175/1520-0469(2003)060,0354:RSDIDC.2.0.CO;2. Chang, W.-Y., T.-C. C. Wang, and P.-L. Lin, 2009: Characteristics of the raindrop size distribution and drop shape relation in typhoon systems in the western Pacific from the 2D video disdrometer and NCU C-band polarimetric radar. J. Atmos. Oceanic Technol., 26, 1973–1993, doi:10.1175/ 2009JTECHA1236.1. Chen, B., J. Yang, and J. Pu, 2013: Statistical characteristics of raindrop size distribution in the Meiyu season observed in eastern China. J. Meteor. Soc. Japan, 91, 215–227, doi:10.2151/ jmsj.2013-208. Chu, Y.-H., and C.-L. Su, 2008: An investigation of the slope-shape relation for gamma raindrop size distribution. J. Appl. Meteor. Climatol., 47, 2531–2544, doi:10.1175/2008JAMC1755.1. Fernández-Raga, M., and Coauthors, 2010: The kinetic energy of rain measured with an optical disdrometer: An application to splash erosion. Atmos. Res., 96, 225–240, doi:10.1016/ j.atmosres.2009.07.013. Frasson, R. P. M., L. K. da Cunha, and W. F. Krajewski, 2011: Assessment of the Thies optical disdrometer performance. Atmos. Res., 101, 237–255, doi:10.1016/j.atmosres.2011.02.014. Friedrich, K., S. Higgins, F. J. Masters, and C. R. Lopez, 2013: Articulating and stationary PARSIVEL disdrometer measurements in conditions with strong winds and heavy rainfall.
633
J. Atmos. Oceanic Technol., 30, 2063–2080, doi:10.1175/ JTECH-D-12-00254.1. Glickman, T. S., Ed., 2000: Glossary of Meteorology. 2nd ed. Amer. Meteor. Soc., 855 pp. [Available online at http:// glossary.ametsoc.org/.] Jaffrain, J., and A. Berne, 2012: Quantification of the small-scale spatial structure of the raindrop size distribution from a network of disdrometers. J. Appl. Meteor. Climatol, 51, 941–953, doi:10.1175/JAMC-D-11-0136.1. Jameson, A. R., M. L. Larsen, and A. B. Kostinski, 2015a: Disdrometer network observations of finescale spatial–temporal clustering in rain. J. Atmos. Sci., 72, 1648–1666, doi:10.1175/ JAS-D-14-0136.1. ——, ——, and ——, 2015b: On the variability of drop size distributions over areas. J. Atmos. Sci., 72, 1386–1397, doi:10.1175/ JAS-D-14-0258.1. Joss, J., and A. Waldvogel, 1967: Ein Spektrograph für Niederschlagestrophen mit automatisher Auswerting (A spectrograph for rain drops with automatic analysis). Pure Appl. Geophys., 68, 240–246, doi:10.1007/BF00874898. Jung, S. A., D. I. Lee, B. J. D. Jou, and H. Uyeda, 2012: Microphysical properties of maritime squall line observed on June 2, 2008 in Taiwan. J. Meteor. Soc. Japan, 90, 833–850, doi:10.2151/jmsj.2012-516. Kumar, L. S., Y. H. Lee, and J. T. Ong, 2011: Two-parameter gamma drop size distribution models for Singapore. IEEE Trans. Geosci. Remote Sens., 49, 3371–3380, doi:10.1109/ TGRS.2011.2124464. Lanzinger, E., M. Theel, and H. Windolph, 2006: Rainfall amount and intensity measured by the Thies laser precipitation monitor. WMO Technical Conf. on Meteorological and Environmental Instruments and Methods of Observation 2006, Geneva, Switzerland, WMO, 9 pp. [Available online at https://www.wmo.int/ pages/prog/www/IMOP/publications/IOM-94-TECO2006/ 3(3)_Lanzinger_Germany.pdf.] Maki, M., T. D. Keenan, Y. Sasaki, and K. Nakamura, 2001: Characteristics of the raindrop size distribution in tropical continental squall lines observed in Darwin, Australia. J. Appl. Meteor., 40, 1393–1412, doi:10.1175/1520-0450(2001)040,1393: COTRSD.2.0.CO;2. Marzano, F. S., D. Cimini, and M. Montopoli, 2010: Investigating precipitation microphysics using ground-based microwave remote sensors and disdrometer data. Atmos. Res., 97, 583– 600, doi:10.1016/j.atmosres.2010.03.019. May, P. T., V. N. Bringi, and M. Thurai, 2011: Do we observe aerosol impacts on DSDs in strongly forced tropical thunderstorms? J. Atmos. Sci., 68, 1902–1910, doi:10.1175/ 2011JAS3617.1. Meng, Z., D. Yan, and Y. Zhang, 2013: General features of squall lines in east China. Mon. Wea. Rev., 141, 1629–1647, doi:10.1175/ MWR-D-12-00208.1. Moumouni, S., M. Gosset, and E. Houngninou, 2008: Main features of rain drop size distributions observed in Benin, West Africa, with optical disdrometers. Geophys. Res. Lett., 35, L23807, doi:10.1029/2008GL035755. Nzeukou, A., H. Sauvageot, A. D. Ochou, and C. M. F. Kebe, 2004: Raindrop size distribution and radar parameters at Cape Verde. J. Appl. Meteor., 43, 90–105, doi:10.1175/ 1520-0450(2004)043,0090:RSDARP.2.0.CO;2. Parker, M. D., and R. H. Johnson, 2000: Organizational modes of midlatitude mesoscale convective systems. Mon. Wea. Rev., 128, 3413–3436, doi:10.1175/1520-0493(2001)129,3413: OMOMMC.2.0.CO;2.
634
JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY
Sarkar, T., S. Das, and A. Maitra, 2015: Assessment of different raindrop size measuring techniques: Inter-comparison of Doppler radar, impact and optical disdrometer. Atmos. Res., 160, 15–27, doi:10.1016/j.atmosres.2015.03.001. Sharma, S., M. Konwar, D. K. Sarma, M. C. R. Kalapureddy, and A. R. Jain, 2009: Characteristics of rain integral parameters during tropical convective, transition, and stratiform rain at Gadanki and its application in rain retrieval. J. Appl. Meteor. Climatol., 48, 1245–1266, doi:10.1175/2008JAMC1948.1. Streets, D. G., C. Yu, Y. Wu, M. Chin, Z. Zhao, T. Hayasaka, and G. Shi, 2008: Aerosol trends over China, 1980–2000. Atmos. Res., 88, 174–182, doi:10.1016/j.atmosres.2007.10.016. Testud, J., S. Oury, P. Amayenc, and R. A. Black, 2001: The concept of ‘‘normalized’’ distributions to describe raindrop spectra: A tool for cloud physics and cloud remote sensing. J. Appl. Meteor., 40, 1118–1140, doi:10.1175/1520-0450(2001)040,1118: TCONDT.2.0.CO;2. Thies Clima, 2007: Instructions for use: Laser Precipitation Monitor 5.4110.xx.x00 V2.4x STD. Adolph Thies GmbH and Co., 64 pp. Thurai, M., V. N. Bringi, and P. T. May, 2010: CPOL radar-derived drop size distribution statistics of stratiform and convective rain for two regimes in Darwin, Australia. J. Atmos. Oceanic Technol., 27, 932–942, doi:10.1175/2010JTECHA1349.1. Tokay, A., and P. G. Bashor, 2010: An experimental study of smallscale variability of raindrop size distribution. J. Appl. Meteor. Climatol., 49, 2348–2365, doi:10.1175/2010JAMC2269.1. ——, A. Kruger, and W. F. Krajewski, 2001: Comparison of drop size distribution measurements by impact and optical
VOLUME 55
disdrometers. J. Appl. Meteor., 40, 2083–2097, doi:10.1175/ 1520-0450(2001)040,2083:CODSDM.2.0.CO;2. Uijlenhoet, R., M. Steiner, and J. A. Smith, 2003: Variability of raindrop size distributions in a squall line and implications for radar rainfall estimation. J. Hydrometeor., 4, 43–61, doi:10.1175/1525-7541(2003)004,0043:VORSDI.2.0.CO;2. Ulbrich, C. W., 1983: Natural variations in the analytical form of the raindrop size distribution. J. Climate Appl. Meteor., 22, 1764–1775, doi:10.1175/1520-0450(1983)022,1764:NVITAF.2.0.CO;2. ——, and D. Atlas, 1998: Rainfall microphysics and radar properties: Analysis methods for drop size spectra. J. Appl. Meteor., 37, 912–923, doi:10.1175/1520-0450(1998)037,0912: RMARPA.2.0.CO;2. ——, and ——, 2007: Microphysics of raindrop size spectra: Tropical continental and maritime storms. J. Appl. Meteor. Climatol., 46, 1777–1791, doi:10.1175/2007JAMC1649.1. Vivekanandan, J., G. Zhang, and E. Brandes, 2004: Polarimetric radar estimators based on a constrained gamma drop size distribution model. J. Appl. Meteor., 43, 217–230, doi:10.1175/ 1520-0450(2004)043,0217:PREBOA.2.0.CO;2. Yuter, S. E., D. E. Kingsmill, L. B. Nance, and M. Löffler-Mang, 2006: Observations of precipitation size and fall speed characteristics within coexisting rain and wet snow. J. Appl. Meteor., 45, 1450–1464, doi:10.1175/JAM2406.1. Zhang, G., J. Vivekanandan, E. A. Brandes, R. Meneghini, and T. Kozu, 2003: The shape–slope relation in observed gamma raindrop size distributions: Statistical error or useful information? J. Atmos. Oceanic Technol., 20, 1106–1119, doi:10.1175/1520-0426(2003)020,1106:TSRIOG.2.0.CO;2.