This study presented airborne measurements of ice particle properties in three stratiform precipitating clouds over northern China. By using horizontal observations at selected altitudes, the distributions of ice water content (IWC), particle habits, and particle size spectrum parameters were investigated. The cloud cases were characterized by high IWC values due to the existence of embedded convective cells. Liquid water contents were rather low with the maxima of less than 0.3 g m−3 and the general values of less than 0.1 g m−3. The occurrence of large dendritic crystals as well as rimed capped columns and branched crystals indicated that ice seeding from the above cloud layer (6 km altitude or above) contributed significantly to both high ice crystal number concentrations and IWCs. Horizontal observations at selected levels suggested the general decreasing trend of IWC with decreasing temperature only in part of the cloud layers but not throughout the cold layer of the multilayered stratiform clouds. Both exponential and gamma functions were used to characterize the particle size spectrum parameters. The slope parameter values of exponential distributions were primarily in the range of 103–104 m−1. In comparison, slope values of the gamma distribution fits spanned more and a relationship was found between the dispersion and slope values.

1. Introduction

Ice particles in ice or mixed-phase clouds have been investigated in many previous studies [13], and they are complex in concentrations, habits, and relative growth modes, affecting precipitation processes and radiative transfer. Furthermore, the study of ice particle size distributions is critical for developing parameterizations for mesoscale and climate models [4].

The extensive use of airborne cloud probes has made it possible to quantitatively examine particle properties by measuring cloud particle sizes and images. Sekhon and Srivasta [5] concluded that the size distribution of snowflakes can be described using the exponential relation. By measuring the size distribution of precipitation particles in frontal clouds, Houze Jr. et al. [6] also found that particles with sizes larger than 1.5 mm generally follow an exponential size distribution. Since devising the spiral decent flight plan by Lo and Passarelli Jr. [7], it has been improved and used in studies of ice particle spectra evolution [810]. Heymsfield et al. [11] and Field et al. [12] presented parameterizations to estimate moments of snow size distributions that can be used in numerical models. Woods et al. [13] subgrouped size distributions of snow particles according to the habit composition from airborne imagery, which significantly improved correlations between the size spectrum parameters and temperature. Aircraft observations of particle size spectra were parameterized in cloud-resolving models in the 1980s in China. Case studies of cloud particle spectra in stratiform clouds over northern China [14, 15] have also been conducted in recent years, but there are still large uncertainties in quantifying cloud microphysical properties [16].

Apart from the above particle information, ice water content (IWC) is also an important parameter to characterize cloud microphysical properties. There are two ways to obtain ice water content from airborne data: one way is in situ measurements based on condensed- and gas-phase water [17] and the other way is by calculation from measured ice particle sizes and images [1820]. For airborne observations without direct IWC measurement, mass-dimensional relationships for ice particles are very useful to analyze microphysical structure and particle growth modes. Fleishauer et al. [21] investigated IWC profiles within six midlevel clouds by using the mass-dimensional expressions from Mitchell et al. [18]. Brown and Francis [19] introduced an improved method to estimate IWC from the mass-dimensional relationships. Gallagher et al. [22] discussed microphysical variations in tropical anvil cirrus outflow regions with the IWC calculated from Heymsfield et al. [23]. The vertical profiles of IWC in midlatitude mixed-phase clouds have also been studied with such method of IWC estimation [24, 25].

Since mixed-phase stratiform clouds are important in producing precipitation over northern China, numerous aircraft observations have been conducted since the 1980s [16, 26, 27]. However, horizontal traverses at a variety of altitudes within one cloud system did not occur frequently. The objective of this study is to investigate the characteristics of IWC distribution and ice particle size spectrum parameters at selected altitudes by using the aircraft data collected in 2009 and 2010. Ice crystal habits at corresponding altitudes are also compared to provide further insight into particle growth mechanisms.

2. Instrumentation and Data Processing

2.1. Instrumentation

Two Y-12E aircraft, from the Shanxi and Beijing Weather Modification Bureaus, equipped with cloud microphysical probes manufactured by droplet measurement technology (DMT) were used to conduct aircraft observations. The Shanxi aircraft was instrumented with a cloud droplet probe (CDP) with a size range of 2–50 m, a 2D cloud imaging probe (CIP) with a size range of 25–1550 m, and a precipitation imaging probe (PIP) with a size range of 100–6200 m. The CIP has a laser diode array of 64 elements, but particles shadowing an end diode are rejected. Therefore, up to 62 slices compose a particle image. In addition, the individual elements of the diode array must be 70% occulted to trigger the digital electronics for an accepted particle. The Beijing aircraft was equipped with a cloud, aerosol, and precipitation spectrometer (CAPS), including a cloud and aerosol spectrometer (CAS) with a size range of 0.6–50 m, and with a CIP and PIP with the same size ranges as those in the Shanxi aircraft. In addition to the above probes, the two aircraft were also equipped with cloud condensation nuclei counter (CCN) and LWC-100, but data from those probes were not considered in the study due to poor baselining before observations.

2.2. IWC Calculation

For cloud particle observations using imaging probes, the number of fragments in the images of shattered particles may reach several hundreds, so the shattering effect should be considered during data analysis [28, 29]. Field et al. [30] found that estimated IWC could be overestimated by 20–30% for narrow size distributions due to the shattering process. And they demonstrated that fragmented particles could be filtered by considering interarrival times. A large quantity of particles smaller than 100 m was found during some flights and concentrations for those small particles were not included in particle size distributions [31]. Due to lack of reliable statistics on shattering effects for aircraft observations in China, only particles larger than 100 m were considered valid for IWC and particle size parameter calculation in this study. The IWC estimation was based on the combined measurements with particles smaller than 1000 m measured by the CIP and particles larger than 1000 m measured by the PIP.

Ice water content can be calculated based on the mass-dimensional relationship of Heymsfield et al. [23] (H04): where (g) represents the mass of an ice particle and (cm) is the maximum dimension. This method produced a very good fit for , but it was also applied to larger particles.

Noh et al. [32] found that H04 values were typically 60–80% smaller than the Locatelli and Hobbs [33] values. The estimates from Mitchell [18] have discrepancies in a factor of 2 [21]. Since no direct IWC measurements were available during our field observations, we could not conclude which method had the smallest errors. Considering that the H04 method was based on the measured IWCs, we chose the estimation with the H04 relationship.

2.3. LWC Calculation

The forward scattering spectrometer probe (FSSP) can detect small particle size distributions from the amount of forward scattered light. In comparison, the CAS determines particle sizes from both forward and backward scattered light. And the CDP detects scattered light in all directions when particles pass through the laser beam. For mixed-phase conditions, droplet spectra measured from the FSSP were inevitably contaminated by ice particles and shattering. Cober et al. [34] found that the FSSP measured droplet spectra could be significantly biased by ice crystals for the size range above 35 m. Further evaluation by Cober et al. [35] demonstrated that the derived LWC agreed with the measured LWC from other probes and within the errors expected from such comparisons. The FSSP measurements of concentration and LWC were found to agree to within 34% and 38%, respectively. Therefore, LWC estimated from the FSSP was still used in previous studies [21, 24].

Compared with the FSSP, the improved measurement techniques used in the CAS minimized observation uncertainties [36], although shattering of ice crystals still occurred on the CAS inlet [37]. In comparison, the CDP had no inlet and, differed from the optical mask of the FSSP annulus detector, generally produced less degree of particle shattering and overcounting than the FSSP [38]. The estimated bias for LWC from the CDP was concentration dependent and it could be as low as −80% [39]. It should also be noted that both CAS and CDP have difficulty in distinguishing small ice particles from droplets during observations. However, contamination of small ice particles was not significant in mixed-phase clouds, as most ice particles grew rapidly to larger than 25 m in less than one minute [40]. Based on above research, it was considered reasonable to estimate LWC from the CAS and CDP probes.

3. Flight Paths

Eight horizontal flight traverses and two ascents within the cold layer associated with three precipitating clouds over northern China were selected in the study. All three cloud cases were multilayered clouds that consisted of primarily altostratus and stratocumulus. The corresponding horizontal flight information is listed in Table 1.

In total, four horizontal flight legs were conducted during 20 and 21 April 2010, including two flight legs in the morning and two other legs in the afternoon. It is worth noting that both legs 2 and 3 were conducted during 1645 and 1710 CST for the same cloud case, except that they were at different horizontal levels. The target region of another three flight legs in 2009 was about 350 km northeast of the observation region in 2010. Similarly, legs 5 and 6 were for the same cloud case except at different horizontal levels. After leg 5, the Beijing aircraft ascended to higher levels and conducted one shorter horizontal observation at 4.9 km. For the cloud case on 18 April 2009, the clouds were in the initialization stage, so no rain was on the surface yet by the time of aircraft observation. The two ascents with the first one from 3.7 to 5.8 km (−1 to −11°C) and the other one from 3.7 to 6.2 km (−1 to −12°C) for the 100420 cloud case were also used to complement horizontal measurements.

4. Horizontal Variation of IWC, LWC, and Particle Habits

The variation of IWC, LWC, and particles habits for the 100420 cloud case was examined in this section. The stratiform clouds on 20 April 2010 were associated with an upper-level trough and surface low-pressure center that moved from southwestern to eastern China. Light-to-moderate rain continued throughout 20 April 2010 and stopped by the early afternoon of the next day over Shanxi Province. Additional radar data in the higher ice region would be more helpful, but they were not available. Therefore, cloud development was shown by observations from a C-band doppler radar with the maximum range of 150 km at Taiyuan (longitude: 112.6°; latitude: 37.7°; ASL: 817.0 m).

Figure 1 shows three plan position indicator (PPI) images and one range height indicator (RHI) image observed by the radar at four selected times (China Standard Time, CST = UTC + 8). Figure 1(a) shows that, at 1031 CST 20 April, rather weak and scattered echo regions around Taiyuan were observed, with one relatively larger region to the southwest of Taiyuan. During 1000 and 1200 CST, light rain with 2 h accumulated rainfall of less than 0.5 mm was observed around Taiyuan.

By 1625 UTC 20 April (Figure 1(b)), the clouds had been more organized with the general reflectivity of 20–30 dBZ and the maximum values increasing to 35–41 dBZ. The 2 h accumulated precipitation during 1600 and 1800 CST increased to 1–6 mm. By 1036 CST 21 April (Figure 1(c)), the clouds had been in the dissipating stage, showing two small separated reflectivity regions with the maximum values of only 15–20 dBZ. The RHI image at 1704 CST 20 April (Figure 1(d)) suggests the typical radar bright band with values of about 30 dBZ at approximately 1.6 km above ground level, indicating weak updrafts inside clouds. The above radar information demonstrated that the 100420 case fell in the range of stratiform cloud precipitation.

Figure 2 shows the time series of IWC and LWC determined by the combined CIP and PIP across the horizontal flight leg 1 on 20 April 2010. The data have been averaged over 10 s to reduce spurious variability. From Figure 2, it can be seen that IWC fluctuated from the peak values of over 0.35 g m−3 to relatively lower values of less than 0.05 g m−3, suggesting spatial inhomogeneity of cloud development. In comparison, LWC also varied inside the cloud, including the maximum value of 0.23 g m−3 at about 1031 CST and extremely low values of less than 0.01 g m−3 at other times. High IWC did not always correspond to high LWC, as many other factors, such relative humidity, cloud temperature, and cloud dynamics, also had impacts on ice particle growth modes. For comparison, two short time periods, one characterized by both high IWC and LWC during 1031 and 1035 CST and the other with both lower IWC and LWC (<0.01 g m−3) during 1044 and 1046 CST, were chosen. The particle images observed during these two periods are shown in Figure 3.

For the first period (Figure 3(a)), there were needles, columns, and unidentifiable crystals in small sizes. It is notable that dendrites and aggregates of dendrites were dominant in this region, with sizes of many aggregates even larger than the maximum range of the CIP. Well-distinguished branches of dendrites and their assemblages from the PIP suggested that aggregation was the dominant growth mode in that region (about −2 to −1°C). According to past laboratory studies, dendrites should occur most frequently at −15 to −10°C. However, a significant number of dendritic particles were observed at −10 to 0°C level [41], indicating falling of ice particles from the above layer. Therefore, the large dendrites in the study were assumed to fall from colder temperature levels at around or above the 6 km altitude.

For the other region observed during 1044 and 1046 CST (Figure 3(b)), particles sizes were much smaller, with needles, columns, and irregular rimed crystals observed. As the two time periods were both in the environmental temperature of about −2 to −1°C, higher supersaturation with respect to ice aloft during 1031 and 1035 CST was considered as the main factor leading to the fast growth of large dendrites.

The estimated IWC and LWC during the horizontal legs 2 and 3 are shown in Figure 4. IWC values at 4.3 km (around −3.5°C, Figure 4(a)) were significantly high, generally in the range of 1.0–3.0 g m−3 and the maximum value of 3.1 g m−3. LWC values at that height were only 0.02–0.06 g m−3. IWC values at lower 3.7 km level (Figure 4(b)) also exhibited high values of 1.1–4.2 g m−3. Overall, LWC values at this level were much lower, with some extremely low values of even less than 0.01 g m−3.

Figure 5 shows the ice particle images at selected times during legs 2 and 3. It can be seen in Figure 5(a) that ice particles were predominantly needles, small hollow columns, and combination of needles at 4.3 km at around 1653 CST. Few rimed crystals were also present. For the other region observed during 1700–1710 CST, significant riming was observed, including such rimed particles as branched crystals, radiating assemblage of plates, capped columns, and irregulars. For the 3.7 km level (Figure 5(b)), predominant needles were also present at 165209 CST, while combination of capped columns and branched crystals with obvious riming existed at the other region. On the whole, a significant number of large heavily rimed ice crystals were observed at both 4.3 and 3.7 km. The capped columns were assumed to form as columns at higher levels (−25–−20°C) fell through the planar-crystal region (−18–−12°C), suggesting that the initial columns appeared at least at 7.5 km altitude. In addition, the crystal number concentrations at 4.3 and 3.7 km were over 50 and 100 L−1, respectively, which were much higher than the values calculated from the ice nucleation formula. Therefore, ice seeding from above the cloud layer contributed to the formation of high IWCs.

As shown in Figures 3 and 5, needles dominated along parts of the flight tracks at temperatures . This was in agreement with previous observations [41] which found that needles occurred most frequently at temperatures between 0 and −5°C. Therefore, the many detected needles from legs 1 to 3 formed initially at that temperature range and later grow primarily through diffusional growth mechanism.

On the whole, the cloud case was characterized by high IWC values at temperatures above −5°C, with magnitudes from 0.1 to 5 g m−3. In contrast, LWC values were generally rather low with peak values of less than 0.3 g m−3 and a small magnitude of 0.01 g m−3. The existence of large dendritic crystals, capped columns, and branched crystals with some degree of riming indicated that ice seeding from the above cloud layer (6 km altitude or above) contributed significantly to both high ice crystal number concentrations and IWCs.

5. Vertical Distribution of IWC with Temperature

IWCs in cirrus generally decreased with decreasing temperature; however, they were considerably scattered, dependent on relative humidity, temperature, vertical velocity, and other particle characteristics [42, 43]. The IWC in single layered mixed-phase clouds maximized in the mid or lower parts of clouds, while higher IWC values might occur in the top of each cloud layer for multilayered clouds [21]. To show the variation of IWC with temperature, Figure 6 presents the IWCs and ice particle concentrations from all the eight horizontal legs listed in Table 1.

According to Figure 6(a), the observations from legs 2 and 3 of the 100420 cloud case suggested general decrease of IWC with decrease of temperature from −1 to −4°C. For the 090501 cloud case, maximum IWC values also decreased from about 2.0 g m−3 at −3°C to the magnitude of 0.1 g m−3 at −8°C. Comparison of IWCs from legs 5, 6, and 7 suggested that although IWC had decreased to less than 0.1 g m−3 at −8°C, higher values of between 0.2 and 0.8 g m−3 at even lower −10°C still could be observed, since the two levels belonged to different cloud layers. In addition, predominant low IWC values of less than 0.25 g m−3 were also observed at various levels due to different phase of cloud dynamics.

Ice particle concentrations (Figure 6(b)) varied from the small magnitude of 1 L−1 to the significant large values of up to 102 L−1 throughout the layers below −11°C. There was great variability even along the same horizontal flight leg. Ice particle concentrations of up to 102 L−1 were also observed in stratiform clouds over Northern China due to the existence of embedded convective cells and such process as ice multiplication [26]. Hobbs and Rangno [44] found the high ice particle concentrations of 40 L−1 within one stratocumulus layer as there was embedded convection within the cloud.

Therefore, IWC values varied significantly at the same temperature level due to different cloud development stages. The general decrease of IWC with increase of height could be seen within one cloud layer at one certain stage but not throughout the cold layer of multilayered clouds.

In addition to above horizontal flights, two vertical profiles were also obtained for the 100420 cloud case, including one ascent from 3.7 to 5.8 km during 1057 and 1122 CST and the other ascent from 3.7 to 6.2 km during 1723 and 1750 CST. The profiles of IWC, ice particle concentrations from the CIP, and cloud droplet concentrations from the CDP are shown in Figure 7. All the values were averaged over 100 m in the vertical. From Figure 7(a), it can be seen that IWC values in the morning were mostly between 0.1 and 0.6 g m−3, corresponding to the general ice particle concentrations of 10–40 L−1 (as shown in Figure 7(b)). The peak IWC value at −8°C corresponded to the high ice particle concentrations of up to 64 L−1. By the afternoon, IWC values at all levels generally had increased significantly to over 0.8 g m−3 at most levels. During that stage, ice particle concentrations also increased to over 30 L−1 or even up to 151 L−1 at −5.5°C.

The overall IWC profiles at these two stages suggested little correlation between IWC and temperature, so IWC does not always decrease with decreasing temperature and was linked more to cloud dynamics which produced different ice particle concentrations. Cloud droplet concentrations at both stages had the extremely low values, with the magnitude of 0.1 cm−3 except one peak value of 27 cm−3 at −9.1°C.

To show the ice particle characteristics more clearly, Figure 8 presents the 2D images from the CIP during those two vertical ascents. It should be noted that apart from the large particles shown in Figure 8, large numbers of small irregulars also existed in both stages. During the first ascent (Figure 8(a)), large particles were primarily hexagonal plates between the −9 and −11°C level and rimed dendrites below the −9°C level. By the time of the second ascent (Figure 8(b)), particles not only become larger in size, but also changed a lot in habits. Many heavily rimed capped columns, graupel-like snow, and plates appeared at various levels. In addition, aggregates of capped columns and plates also occurred. Based on particle fall-speed measurements by Locatelli and Hobbs [33], densely rimed dendrites, 500–1000 m in diameter, are assumed to have fall velocities of 2.2 to 2.8 m s−1. In comparison, graupel-like snow of lump type with the same diameter should have fall velocities of between 3.3 to 4.0 m s−1 and densely rimed columns have much higher velocity of around 10 m s−1. The updrafts in stratiform clouds are generally in the magnitude of 0.1 m s−1. Therefore, the capped columns and graupel-like snow in the second ascent were considered to have higher net velocities than those in the first ascent, which contributed to higher IWCs at various temperature levels.

6. Particle Size Distribution and Spectrum Parameters

Ice particle size distributions within clouds over northern China have been fitted in previous studies following an exponential form [14]. However, the data were obtained in the 1980s using OAP-2D-C and OAP-2D-P. The use of new instrumentation in our aircraft measurements is better in providing particle size information. Both exponential and gamma functions were used in the study to parameterize ice particle size distributions. The equation that describes the gamma distribution is where , , and are intercept, dispersion, and slope parameters to be derived. When , the gamma fit becomes an exponential function:

As particle size distributions are highly variable from one cloud region to another, Figure 9 firstly shows the average particle size distributions for different legs. Size distributions in Figure 9(a) included particles below 1000 m detected by the CIP and above 1000 to 6200 m detected by the PIP. It can be seen that measurements from the two probes overlapped quite well at 1000 m and most particles were smaller than 2000 m. For particles larger than 500 to 1000 m, the concentration (plotted on logarithmic axes) decreased linearly with size (plotted on linear axes), suggesting that above these sizes the particle size distributions were of exponential shape. In contrast, for particles smaller than 500 to 1000 m, a linear decrease is not so obvious when plotted on the logarithmic axes. Average particle concentrations measured by the CDP during legs 3, 4, and 6 are plotted in Figure 9(b), which varied significantly for different flight legs. The CDP data were not combined with the CIP measurement for two reasons. On the one hand, the Beijing aircraft was equipped with the CAS, not the CDP. On the other hand, particles in the size range 25–100 m from the CIP detection were not used in the study, so there was no overlap between 50 (the maximum size bin for the CDP) and 100 m.

Figure 10 shows the versus and versus points for the various horizontal legs. Each data point corresponded to a CIP and PIP combined size distribution averaged over 10 s. For the exponential fit (Figure 10(a)), values of generally had the magnitude of 107 to 109 m−4 and values spanned over one order of magnitude. The data were clustered into two groups, including legs 2, 3, 5, 6, and 7 with higher and lower values, as well as legs 1, 4, and 8 with lower and relatively higher values. The variation indicated that ice particles during legs 2, 3, 5, 6, and 7 were larger and more concentrated in size, which was in agreement with the occurrence of many rimed and aggregated particles. In comparison, the values (Figure 10(b)) also centered around 103 m−1 but spanned to over two orders of magnitude.

Plots of versus for various legs are shown in Figure 11. The values of were between −3 and 4, including evident fluctuations. A relationship similar to that of Heymsfield [4] between and was fitted, with the values varying from positive values at large to negative values at small . However, due to significant variability in particle size distributions, it was difficult to find an accurate relationship between and . To compare the fits from exponential and gamma functions, Figure 12 presents correlation coefficients () versus . The correlation coefficients were mostly between 0.6 and 0.9 for both exponential and gamma distributions and also included some scattered lower values. The value averaged over all exponential fits was 0.68. Similarly, the value for gamma fits was 0.66. Therefore, the two functions had similar degree of accuracy in fitting particle size distributions.

In addition, plots of particle size spectrum parameters versus cloud temperature and IWC were also examined, but no accurate relationship was found. It is speculated that the data only in a small temperature range of 0 to −5°C from three cloud cases is not enough. To better fit particle size distributions and other related properties, more data of higher quality are still needed.

7. Discussion and Conclusions

This study presented the measurements of ice particle properties observed from horizontal flight legs at selected altitudes as well as from ascent legs for three stratiform precipitating cloud cases. The distributions of IWC, particle habits, and particle size spectrum parameters were examined.

IWC distribution was first analyzed by comparing IWC, LWC, and particle habits from horizontal observations for the 100420 cloud case. Spatial inhomogeneity of IWC was very common in stratiform clouds. High IWC values varying with magnitudes from 0.1 to 5 g m−3 were detected below the −5°C level, while LWCs were less than 0.3 g m−3. Needles and hollow columns dominated along parts of flight tracks, indicating that the diffusional growth mechanism was important for the production of relatively smaller ice crystals. At the same time, the occurrence of large dendritic aggregates and rimed particles such as capped columns, branched crystals, and graupel-like snow in the other regions suggested ice seeding from the above layer at around 6-7 km. Those large rimed particles were also observed during the two ascents, whose relatively higher fall velocities contributed to the appearance of high IWCs at lower levels.

In addition, high IWC also attributed to cloud dynamics. It was speculated that the embedded convection cells contributed to the generation of such high IWC (over 0.5 g m−3), since scattered convective regions with the reflectivity of up to 45 dBZ for the 090501 case and the small reflectivity regions of around 40 dBZ for the 100420 case were observed. By examining layer clouds around the world, Ryan [45] found that deeper stratiform clouds were generally glaciated ones with embedded convection; however, away from regions of embedded convection the cloud water contents were small.

The stratiform precipitating clouds over northern China are usually multilayered clouds with cirrus, altostratus, and stratocumulus. Vertical variability of IWCs suggested the general decrease trend of IWC maxima with decreasing temperature only in part of the cloud layers, such as the maximum value of more than 4.0 g m−3 at −1°C decreasing to around 3.0 g m−3 at −4°C in the 100420 case and the maximum 2.1 g m−3 at −3.5°C decreasing to less than 0.1 g m−3 at −8.5°C in the 090501 case. However, relatively lower IWC values of less than 0.25 g m−3 were also observed throughout the cold layer.

On the whole, values of varied primarily between 107 and 1010 m−4 and values were generally in the range of 103–104 m−1 from exponential distributions. In comparison, both and spanned more for gamma distributions. Further calculation of correlation coefficients suggested that the two functions had the similar degree of accuracy. By analyzing the slope and intercept parameters for 77 horizontal flight legs, Woods et al. [13] found that values of varied over three orders of magnitudes and values spanned less than a full order of magnitude. Heymsfield et al. [11] found the monotonic relationship between and after using exponential curves to fit particle size distributions for tropical clouds. The small sample size in this study did not allow us to conduct further investigation of the relationship between particle size spectrum parameters and temperature, but the initial analyses were still useful to show the overall characteristics of stratiform clouds over northern China.

Only three cloud cases were examined in this study. To obtain the statistical characteristics of cloud microphysical parameters, a larger sample size is needed. Future work such as analyzing errors in observational data and investigating key parameters affecting aggregation and other microphysical processes should also be conducted.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work was supported by the National Basic Research Program of China (973 Program, Grant no. 2013CB430105), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA05100300), the National Natural Science Foundation of China (Grant no. 41105095), and the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant no. KZCX2-EW-203).