The evolution of cloud and aerosol microphysics at the summit of Mt. Tai, China

The influence of aerosols, both natural and anthropogenic, remains a major area of uncertainty when predicting the properties and the behaviours of clouds and their influence on climate. In an attempt to better understand the microphysical properties of cloud droplets, the simultaneous variations in aerosol microphysics and their potential interactions during cloud life cycles in the North China Plain, an intensive observation took place from 17 June to 30 July 2018 at the summit of Mt. Tai. Cloud microphysical parameters were monitored simultaneously with number concentrations of cloud condensation nuclei (NCCN) at different supersaturations, PM2.5 mass concentrations, particle size distributions and meteorological parameters. Number concentrations of cloud droplets (NC), liquid water content (LWC) and effective radius of cloud droplets (reff) show large variations among 40 cloud events observed during the campaign. The low values of reff and LWC observed at Mt. Tai are comparable with urban fog. Clouds on clean days are more susceptible to the change in concentrations of particle number (NP), while clouds formed on polluted days might be more sensitive to meteorological parameters, such as updraft velocity and cloud base height. Through studying the size distributions of aerosol particles and cloud droplets, we find that particles larger than 150 nm play important roles in forming cloud droplets with the size of 5–10 μm. In general, LWC consistently varies with reff. As NC increases, reff changes from a trimodal distribution to a unimodal distribution and shifts to smaller size mode. By assuming a constant cloud thickness and ignoring any lifetime effects, increase in NC and decrease in reff would increase cloud albedo, which may induce a cooling effect on the local climate system. Our results contribute valuable information to enhance the understanding of cloud and aerosol properties, along with their potential interactions on the North China plain. Published by Copernicus Publications on behalf of the European Geosciences Union. 13736 J. Li et al.: The evolution of cloud and aerosol microphysics

Abstract. The influence of aerosols, both natural and anthropogenic, remains a major area of uncertainty when predicting the properties and the behaviours of clouds and their influence on climate. In an attempt to better understand the microphysical properties of cloud droplets, the simultaneous variations in aerosol microphysics and their potential interactions during cloud life cycles in the North China Plain, an intensive observation took place from 17 June to 30 July 2018 at the summit of Mt. Tai. Cloud microphysical parameters were monitored simultaneously with number concentrations of cloud condensation nuclei (N CCN ) at different supersaturations, PM 2.5 mass concentrations, particle size distributions and meteorological parameters. Number concentrations of cloud droplets (N C ), liquid water content (LWC) and effective radius of cloud droplets (r eff ) show large variations among 40 cloud events observed during the campaign. The low values of r eff and LWC observed at Mt. Tai are comparable with urban fog. Clouds on clean days are more sus-ceptible to the change in concentrations of particle number (N P ), while clouds formed on polluted days might be more sensitive to meteorological parameters, such as updraft velocity and cloud base height. Through studying the size distributions of aerosol particles and cloud droplets, we find that particles larger than 150 nm play important roles in forming cloud droplets with the size of 5-10 µm. In general, LWC consistently varies with r eff . As N C increases, r eff changes from a trimodal distribution to a unimodal distribution and shifts to smaller size mode. By assuming a constant cloud thickness and ignoring any lifetime effects, increase in N C and decrease in r eff would increase cloud albedo, which may induce a cooling effect on the local climate system. Our results contribute valuable information to enhance the understanding of cloud and aerosol properties, along with their potential interactions on the North China plain.

Introduction
Clouds are key factors in the atmospheric hydrological cycle, which play an important role in the atmospheric energy budget and significantly influence the global and regional climate (Chang et al., 2019b;Y. Zhang et al., 2004). Clouds can be physically described by their liquid water content (LWC), number concentrations of droplets (N C ) and effective radius of droplets (r eff ). These parameters may show small interannual variations for the same monitoring station (Möller et al., 1996), but they vary over a large range for different cloud types (Quante, 2004), for different cloud altitudes (Padmakumari et al., 2017;Zhao et al., 2018) and in different parts of a cloud (Deng et al., 2009).
The interactions between the clouds and the aerosols are complex. Clouds can efficiently remove aerosols by activating cloud condensation nuclei (CCN) to form cloud droplets (Croft et al., 2010;L. M. Zhang et al., 2004). The cloud processes can incorporate large amount of fine particulate mass (Heintzenberg et al., 1989), change their size distributions Schroder et al., 2015), and alter the CCN compositions through homogeneous and heterogeneous reactions (Roth et al., 2016). In addition, the variation of aerosol number concentrations and aerosol size distributions could alter the cloud microphysics. Through studying microphysical characteristics of cloud droplet residuals at Mt. Åreskutan, Noone et al. (1990) found that larger cloud droplets preferred to form on larger CCN. In addition to this, the aerosol-cloud interaction has been investigated for cloud processes formed under both clean and polluted conditions. Padmakumari et al. (2017) found that convective clouds over land were characterised by lower LWC and higher N C due to the increase in pollution aerosol. Ground-based observations by radiometers during the summers of the US studies in the mid-Atlantic region revealed that cloud events with smaller droplets (< 7 µm) were more frequently observed in the polluted years than in the clean years (S. . The influence of aerosols on the cloud microphysics is evident but varies for different regions and for different cloud types. For a given liquid water content, aerosol particles can act as CCN, lead to higher number concentrations of cloud droplets with smaller sizes and result in higher albedo (Twomey effect or first indirect effect, FIE) (Twomey, 1974). Based on the principle of the Twomey effect, calculations for evaluating the influence of aerosols on the cloud microphysics have been widely studied (Lohmann and Feichter, 2005;McComiskey et al., 2009;Twohy et al., 2005). However, arithmetic terms representing aerosol loading are different, such as using the number concentration of particles, the CCN concentration and the aerosol optical depth (AOD), which makes it difficult to compare the FIE from different studies.
The increase in the aerosol concentrations can result in a longer cloud lifetime, thus producing large cloud fractions (Koren et al., 2005;Albrecht, 1989), increasing cloud top height and increasing the cloud thickness (Fan et al., 2013). It would further influence the regional and global climate (Rosenfeld, 2006;Seinfeld et al., 2016), such as reducing the precipitation or drizzle (Andreae et al., 2004;Heikenfeld et al., 2019) and further delaying the hydrological cycle (Rosenfeld, 2006). Through model experiments with the Coupled Model Intercomparison Project phase 5 (CMIP5), Frey et al. (2017) also found that the addition of anthropogenic aerosols could increase the monthly mean cloud albedo of subtropical marine stratocumulus clouds.
In situ measurements of cloud microphysics by aircraft or on high-altitude monitoring sites have provided some additional information for insight into the cloud processes (Allan et al., 2008;Padmakumari et al., 2017;Reid et al., 1999). However, lacking knowledge of the size distributions of cloud droplets and aerosol particles makes it difficult to evaluate the cloud microphysics in small-scale regions (Fan et al., 2016;Khain et al., 2015;Sant et al., 2013). Discrepancy still exists between the widths of observed and simulated size distributions of cloud droplets (Grabowski and Wang, 2013). In addition to this, incomplete knowledge of the impact of cloud-aerosol interactions (Rosenfeld et al., 2014b), unresolved processes of cloud formation (Stevens and Bony, 2013) and the lack of research into the variation of cloud microphysical parameters at different cloud stages still hinder modelling studies.
In this study, in situ observations at the summit of Mt. Tai were presented to investigate the evolution of cloud microphysics coupled to simultaneous monitoring of aerosol size distributions, PM 2.5 mass and CCN concentrations within non-precipitating clouds. The summit of Mt. Tai is the highest point in the centre of the North China Plain (NCP). Sufficient moisture in summer and dramatic temperature differences between day and night make it ideal for in situ orographic cloud monitoring (J. . The summit of Mt. Tai is far away from anthropogenic emission sources on the ground. However, high concentrations of inorganic ions in PM 2.5 (Zhou et al., 2009), abundant bacterial communities (Zhu et al., 2018), and NH 3 and NO x emissions from biomass burning (Chang et al., 2019a) have been observed at the summit, and thus a strong anthropogenic influence must exist. Previous studies of cloud samples collected at the same position showed high inorganic ion concentrations (J. Wang et al., 2011), which can be attributed to the increase in anthropogenic aerosol. In this study, two typical cloud processes are discussed in detail to elucidate the relationship of N C , r eff , and LWC under clean or polluted conditions (indicated by N P and N CCN ) and during the cloud life cycle. This paper provides comprehensive information about cloud microphysical properties and their potential links to aerosol concentrations and size distribution. Implications of cloud and aerosol microphysics for cloud albedo and climate are discussed.

Observation duration and site
From 17 June to 30 July 2018, 40 cloud events in total were monitored at the Shandong Taishan Meteorological Station at summit of Mt. Tai (Tai'an,China;36 • 18 N,117 • 13 E; 1545 m a.s.l.; Fig. S1). Mt. Tai is the highest point in the central of North China Plain (NCP) and located within the transportation channel between the NCP and the Yangtze River Delta (Shen et al., 2019). The altitude of Mt. Tai is close to 1.6 km. This height is close to the top of the planetary boundary layer in central eastern China and is suitable for studying the characteristics of particles inputting to clouds (Hudson, 2007). Orographic clouds, which are mainly formed in the boundary layer via air approaching the ridge, being forced to rise up and then being cooled by adiabatic expansion (Choularton et al., 1997), frequently occur at the summit of Mt. Tai, especially in summer. Previous studies concentrated on cloud chemistry found that Mt. Tai is significantly influenced by anthropogenic emissions (J. Wang et al., 2011). In addition, fixed observation location is mainly applied to study the evolution of aerosol properties and cloud processing (Mertes et al., 2005;Roth et al., 2016). Thus, Mt. Tai is a good site for monitoring orographic clouds and simultaneously investigating aerosol and cloud microphysics. The arrangement of instruments is presented in Fig. S1c. As shown in Fig. S2, the prevailing wind direction during this summer campaign was easterly wind (23.3 %), southwesterly wind (22.8 %) and southerly wind (21.9 %), respectively. About 85.6 % of wind speed was less than 8 m s −1 . The monitored cloud events in this study were mainly influenced by southerly wind (34.7 %) and southwesterly wind (22 %).

Cloud microphysical parameters
A Fog Monitor (Model FM-120, Droplet Measurement Technologies Inc., USA), a forward-scattering optical spectrometer with a sampling flow of 1 m 3 min −1 , was applied in situ for real-time size distributions of cloud droplets and computing N C , LWC, median volume diameter (MVD), and effective diameter (ED) in the size range of 2 to 50 µm (Spiegel et al., 2012). The corresponding equations are where N i is the cloud number concentration at the ith bin, r i represents the radius at the ith bin and ρ w = 1 g cm −3 stands for the density of liquid water. Droplets are categorised into the manufacturer's 30 predefined size bins with a sampling resolution of 1 s. The size bin widths using this configuration were 1 µm for droplets < 15 and 2 µm for droplets >15 µm. The true air speed calibration and size distribution calibration of FM-120 were carried out by the manufacturer using borosilicate glass microspheres of various sizes (5.0, 8.0, 15.0. 30.0, 40.0 and 50.0 µm, Duke Scientific Corporation, USA). The difference in optical properties between the glass beads and water was taken into account during the calibration process. In this study, the sampling inlet nozzle faced the main wind direction and was horizontally set. Cloud events are defined by the universally accepted threshold values in N C and LWC, i.e. N C > 10 cm −3 and LWC > 0.001 g m −3 (Demoz et al., 1996). Too short cloud events with a duration < 15 min were excluded.

Aerosol size distribution
A Scanning Mobility Particle Sizer (SMPS, Model 3938, TSI Inc., USA) consisting of a Differential Mobility Analyser (DMA, Model 3082, TSI Inc., USA) and a Condensation Particle Counter (CPC, Model 3775, TSI Inc., USA) was applied to monitor the size distributions of dehumidified aerosols through a PM 10 inlet. The neutralised aerosols were classified by DMA to generate a monodisperse stream of known size according to their electrical mobility. The CPC placed downstream counts the particles and gives the number of particles with different sizes. In this study, each scan was fixed at 5 min for every loop with a flow rate of 1.5 L min −1 sizing particles in the range of 13.6-763.5 nm in 110 size bins.

CCN number concentration
The N CCN at certain supersaturations (SSs) were quantified by a Cloud Condensation Nuclei Counter (Model CCN-100, DMT Inc., USA). The CCN counter was set at five SS values sequentially for 10 min each at 0.2 %, 0.4 %, 0.6 %, 0.8 % and 1.0 % with a full scan time resolution of 50 min. Data collected during the first 5 min of each SS were excluded since the CCN counter needs time for temperature stabilisation after the change in SS. The ratio of sample flow to sheath flow was set at 1 : 10 with a total airflow of 500 ccm. The SSs of the CCN counter were calibrated before the campaign and checked at the end of the campaign with monodisperse ammonium sulfate particles of different sizes (Rose et al., 2008).

PM concentrations and meteorological parameters
The PM 2.5 mass concentration was measured using a beta attenuation and optical analyser (SHARP monitor, model 5030i, Thermo Scientific Inc., USA). Meteorological parameters, including the ambient temperature (T a , • ), relative humidity (RH), wind speed (WS, m s −1 ) and wind direction (WD, • ), were provided by Shandong Taishan Meteorological Station at the same observation point. The ground-level temperature (T g ), ground-level pressure (P g ), and dew point temperature (T gd ) were supported by National Meteorological Observatory -Tai'an Station (station number: 54 827; 36 • 9 N, 117 • 9 E; 128.6 m a.s.l.) (Fig. S1a), which is located on the southern plain of Mt. Tai.

Calculation of cloud base height
In this study, the estimated lifting condensation level (LCL) is applied to represent the cloud base height (CBH) due to the lack of corresponding instruments. The calculation of LCL depends on the meteorological parameters measured at Tai'an Station. The ground-level data of temperature, dew point temperature, and pressure were used as input parameters (Georgakakos and Bras, 1984): where p LCL is the LCL pressure, and T LCL is the LCL temperature.

Calculation of AIE
Aerosol indirect effect (AIE), which here simply represents approximations of the derivatives of the cloud microphysics (r eff and N C ) with respect to changes in aerosol concentrations (McComiskey et al., 2009;Feingold et al., 2001), is applied to study the influence of N P on cloud microphysics and calculated as follows: where N P is applied as an proxy of aerosol amount (Zhao et al., 2012;Zhao et al., 2018).

Calculation of cloud albedo
Cloud albedos can be calculated using the equations shown below (Seinfeld and Pandis, 2006). Assuming the cloud droplet size distribution can be approximated as monodisperse and the cloud is vertically uniform with respect to droplet size distribution (Stephens, 1978), the cloud optical thickness (τ c ) could be obtained by where h is the thickness of the cloud and ρ w is the density of cloud water. For the non-absorbing and horizontally homogeneous cloud, the two-stream approximation for the cloud albedo (R c ) gives the following equation (Lacis and Hansen, 1974): where g is the asymmetry factor. The radius of cloud droplets was much greater than the wavelength of visible light, and hence g is 0.85. The equation before becomes 3 Results and discussion

Overview of the cloud microphysics
From 17 June to 30 July 2018, 40 cloud events were captured at the summit of Mt. Tai. Large ranges of cloud microphysics were observed during the campaign. The averaged N C , LWC and r eff of the 40 cloud events at the summit of Mt. Tai varied over the ranges of 59-1519 cm −3 , 0.01-0.59 g m −3 and 2.6-7.4 µm, respectively (Table S1). The monitored number concentration of cloud droplets at Mt. Tai both in this study and in 2014 can reach 2000-3000 cm −3 (J. , which is much higher than those values (with a range of 10-700 cm −3 ) for city fog and convective and orographic clouds (Allan et al., 2008;Li et al., 2011;Padmakumari et al., 2017) ( Table 1). This suggests that clouds at Mt. Tai were characterised with high N C . The microphysics of different clouds and fog could generally be distinguished in a plot of r eff (or MVD) against LWC. As illustrated in Fig. 1, the LWC generally increased on the order of city fog, orographic clouds, and convective clouds, and Mt. Tai generally follows this rule. It was consistent with the study by Penner et al. (2004) that LWC within clouds increases linearly with altitude. For LWC values of clouds at Mt. Tai, both the high values, which were comparable with convective clouds, and the low values, which were similar to city fog (Fig. 1), were monitored. It indicated that clouds at Mt. Tai appeared to show a larger range of LWC values. The increase in LWC at Mt. Tai was determined by the increase in r eff and/or N C . However, sometimes only one factor played the dominant role. As illustrated in Table S1, N C , r eff and LWC in cloud event 20 (CE-20) were 1519 cm −3 , 5.2 µm and 0.54 g m −3 , respectively, while the corresponding values in cloud event 16 (CE-16) were 59 cm −3 , 9.8 µm and 0.14 g m −3 , respectively. Even though r eff of CE-20 was smaller when compared to CE-16, the higher N C determined the larger LWC of clouds in CE-20. In the following sections, the evolution of cloud and aerosol microphysical properties are presented. The influence of meteorological parameters  (such as updraft velocity and cloud base height) and aerosol particle on cloud microphysics are discussed.

Analysis of typical cloud processes
By assuming a density of particles ρ = 1.58 g cm −3 (Cross et al., 2007), the mass concentrations of particles, which were calculated from the aerosol number size distribution measured by SMPS and named PM 0.8 , were highly consistent with PM 2.5 , especially when PM 2.5 was less than 20 µg m −3 (Fig. 2c). Based on the mass concentration (PM 2.5 ) and the number concentration (N P , which represented the total number concentration of aerosol particles measured by SMPS) of aerosols, two typical cloud processes, "cloud process-1" (CP-1) and "cloud process-2" (CP-2), were selected and analysed with respect to their special characteristics. The variations in updraft velocity v up on cloud microphysics during CP-1 and CP-2 were ignored for simplicity (Table S2, Figs. S3 and S4). The sampling angle (θ s ) and v up for CP-1 and CP-2 were 11.9 • and 0.82 ± 0.29 m s −1 and 10.6 • and 0.92 ± 0.36 m s −1 , respectively (Table S2). According to the calculations provided by Spiegel et al. (2012), the aspiration efficiency and transmission efficiency of the fog monitor were all close to 1. In CP-1 (which only contained cloud event 19, CE-19), cloud droplets formed under a relatively stable (wind speed < 4 m s −1 ) and clean (PM 2.5 ≈ 10.9 µg m −3 and N P ≈ 1425 cm −3 ) conditions accompanied by a slow increase in T a (Figs. 2 and 3). During daytime, especially in the afternoon, the PM 2.5 mass concentration dramatically increased with few changes in wind speed and wind direction, meanwhile, N P reached to about 5000 cm −3 (Fig. 3). CP-1 persisted for 74 h, making it the longest cloud event during the presented campaign. Quite different from CP-1, CP-2 contained eight cloud events (CE-20 to cloud event 26, CE-26; see Fig. 3) and occurred periodically under high PM 2.5 (Fig. 2, 50.7 µg m −3 on average) and high N P (Fig. 3, 1694 cm −3 on average) conditions. Cloud events in CP-2 formed after sunset, with a sharp decrease in PM 2.5 and N P , and transitorily dissipated at noon accompanied with the increase in PM 2.5 , N P , T a and cloud base height (CBH). For cloud water samples collected during CP-1 and CP-2, the percentage of chemical compositions did not change a lot (Fig. S5). The total measured dominant ions (sulfate, nitrate and ammonia) were 93.39 % in CP-1 and 90.37 % in CP-2. The high concentration of secondary ions in the cloud water samples indicated that clouds at Mt. Tai were dramatically influenced by anthropogenic emissions.
Each cloud event of CP-2 was separated into an activation stage (S1), collision-coalescence stage (S2), stable stage (S3) and dissipation stage (S4) according to the regular changes in N C and LWC / N C (Fig. 3a). In S1, N C dramatically increased to its maximum value among the cloud events. In S2, N C declined sharply to a stable value, meanwhile LWC / N C reached the maximum value. In S3, N C was stable or slightly varied and LWC / N C started to decrease. In ; (e) size distribution of particles (13.6-763.5 nm) and corresponding geometric mean radius (GMr P ); (f) size distribution of cloud droplets (2-50 µm) and corresponding geometric mean radius (GMr C ); and (g) N C and LWC of cloud droplets. S4, both N C and LWC / N C decreased sharply and finally arrived at zero. Even though the two stages (S2 and S3) in cloud event 25 (CE-25) did not totally follow the division rules, the other six cloud events followed them closely. This indicated that the division was helpful for the study of the variations of cloud microphysical properties during CP-2. The newly formed cloud droplets during S1 were characterised by small size, high N C and low LWC / N C values (Figs. 2f and 3b). For example, about 2310 cm −3 of cloud droplets can quickly form in the first 2 h of CE-20. The r eff of these droplets was smaller than 4.1 µm, and LWC / N C was about 0.2 ng. In going from S2 to S3, the strong collision-coalescence between cloud droplets caused the increase in both r eff and LWC / N C . In S4, the increase in PM 2.5 , through evaporation of cloud droplets or lifting of CBH (Fig. 2), was observed to coincide with the vanishing of cloud events (Mazoyer et al., 2019;.

Relationships among N P , N CCN and N C
In this study, both consistent variation and inverse variation between N P and N C were observed. N P and N C showed consistent variation in CP-1. However, in CP-2, an obviously inverse variation was found between N P and N C in S1 and S4, while a consistent variation was found between N P and N C in S2 and S3 (Figs. 3a, 4b, and c). Some in situ observations (Lu et al., 2007;Mazoyer et al., 2019) and modelling studies (Heikenfeld et al., 2019;Zhang et al., 2014) supported the viewpoint that the increase in N P brings more CCN and further increases N C , which could cause the consistent variation between N P and N C . In contrast, some recent studies of fog also suggested that the increase in N P could decrease the ambient supersaturation and then decrease N C (Boutle et al., 2018;Mazoyer et al., 2019). Besides, Modini et al. (2015) found inverse variation between N C and the number of particles with diameters larger than 100 nm due to the reduction of supersaturation by coarse primary marine aerosol parti- cles. In general, the covariation between N p and N C could be affected by many factors, including competition of water vapour between aerosol particles and/or cloud droplets, the scavenging of particles by cloud droplets, and new particle formation through cloud processes. In this study, consistent variation between N P and N C was characterised with higher LWC / N C , while inverse variation between N P and N C appeared with lower LWC / N C . The average LWC / N C was 0.61 ng in CP-1 and was 0.15, 0.42, 0.39 and 0.16 ng in S1, S2, S3 and S4, respectively, in CP-2. The relatively higher LWC / N C values could indicate conditions where the amount of water vapour was not limiting cloud droplet formation. Once N P increased, part of the cloud water was taken away by the CCN in the particles to form new droplets, and the remaining amount of water was still sufficient to maintain the previous droplets in liquid state. Thus, both N P and N C simultaneously increased. On the other hand, relatively low LWC / N C values could, to some extent, limit the formation of new cloud droplets. The activated particles that grew at the beginning of the cloud cycle would lower the surrounding supersaturation and to some extent limit further aerosol activation (Ekman et al., 2011). The part of the water taken by the CCN in the particles was not enough to change all of them into new droplets, and the remaining amount of water was also insufficient to maintain all of the previous droplets in a liquid state. Following this, the N C would decrease and the higher the N P becomes, the sharper decrease the N C will be. Thus, the inverse variation between N P and N C was observed.
The ratio between N CCN and N P reflects the activation ratio of aerosol particles. As shown in Fig. S6, N CCN increased with the increase in SS. In addition, N CCN of CP-2 was higher than that of CP-1 at the same SS. In order to compare with previous studies as discussed below, SS = 0.2 % was chosen to calculate N CCN / N P , which represented the activation ratio of aerosol particles. As shown in Fig. 3b, N CCN  in CP-1 was smaller than that of 0.38 in CP-2, and values lower than 0.22 did not appear during CP-2. This indicated that the activation of aerosol particles in CP-2 was relatively easy. Both the size distribution and the chemical composition could impact the cloud-nucleating ability of aerosol particles (Dusek et al., 2006;Mazoyer et al., 2019). Figure S7 shows the relation between N CCN,0.2 / N P and GMr P during CP-1 and CP-2. As can be seen in Fig. S7, higher correlation of N CCN,0.2 / N P with GMr P was found during CP-1 than during CP-2, which suggested that the physical properties might have more influence on the activation of aerosols during CP-1. Aside from this, Asmi et al. (2012) found that higher N CCN / N P and a more concentrated plot of N CCN versus N P usually occurred during winter when a higher fraction of aged organics was observed during the observation programme at Puy-de-Dôme, France. In this study, the plot of N CCN,0.2 versus N P was found to be more scattered in CP-1 than in CP-2 (Fig. S8). Even though the chosen SS in this study (SS = 0.2 %) is different from that at Puy-de-Dôme (SS = 0.24 %), most of the data points of CP-1 and CP-2 were distributed between the two recommended dashed lines (the visually defined boundaries within most of the data are centred, Fig. S8) by Asmi et al. (2012). This suggested that the difference in aerosol organic chemical compositions during CP-1 and CP-2 might also explain the different activation ratio of aerosol particles during these two cloud processes.

The influence of N P , CBH and v up on cloud microphysics
No negative AIE r or AIE N were found in this study (Fig. 4). The positive AIE r and AIE N at Mt. Tai mean that the increase in N P is accompanied by decreased r eff and increased N C . However, in the studies of Yuan et al. (2008) and Tang et al. (2014), AOD was applied to represent aerosol loading and negative AIE r (indicating r eff increased with the increasing of AOD) near the coastlines of the Gulf of Mexico, the South China Sea and over eastern China. The reason proposed by Yuan et al. (2008) was the increasing slightly soluble organic (SSO) particles, which would increase the critical supersaturation and hinder the activation of the particles. Meanwhile, Tang et al. (2014) showed that the meteorological conditions, which favoured the transportation of pollutants and water vapour, led to simultaneous increases in both AOD and r eff . Different from the coastal area, the summit of Mt. Tai is relatively far from the sea (around 230 km from the Bohai Sea and Yellow Sea) (Guo et al., 2012), leading to less moisture in the air. This might hinder the growth of cloud droplets and cause the positive AIE r . In addition, the increase in LWC was found to covary with the decrease in AIE, especially at coastal sites (McComiskey et al., 2009;Zhao et al., 2012). However, non-obvious variation was found between AIE r and LWC at Mt. Tai (Fig. 4a). It might due to the high aerosol loading during cloud processes (Zhao et al., 2012). Although all positive AIE r and AIE N were found in CP-1 and CP-2, the specific values were different. According to the studies of AIE r and AIE N of CP-1 and CP-2, our results suggested that cloud droplet numbers were more sensitive to N P under conditions with lower aerosol concentrations. The calculation of AIE r was shown in Fig. S9 and summarised in Fig. 4. As shown in Fig. 4a, except for the out-of-bounds AIE r values calculated with insufficient data points when LWC was larger than 0.7 g m −3 , AIE r of 0.181-0.269 for CP-1 was always higher than that of 0.025-0.123 for CP-2 in corresponding narrow LWC ranges. When calculating AIE N , the number of cloud droplets may be underestimated during the activation and dissipation stages (Mazoyer et al., 2019), due to the limitations of the fog monitor, which could be a cause of the low R 2 of AIE N in CP-1. In CP-2, only the data of S2 and S3 were employed to calculate AIE N for excluding the points in S1 and S4, which may be an underestimation. As shown in Fig. 4b and c, both the slope (0.144) and R 2 (0.050) of CP-2 are lower than those (0.544 and 0.282, respectively) of CP-1. It was verified that cloud droplets in CP-2 were influenced little by aerosols. In the previous studies, both observation and modelling studies also found that AIE r was higher under smaller aerosol amount conditions. Twohy et al. (2005) measured the equivalent AIE r of 0.27 on the California coast, while Zhao et al. (2018) used satellite observations to attribute lower values of 0.10-0.19 for convective clouds over Hebei, a polluted region in China. Using an adiabatic cloud parcel model, Feingold (2003) found that AIE r increased from 0.199 to 0.301 when N P decreased to less than 1000 cm −3 . By using the Community Atmospheric Model version 5 (CAM5), Zhao et al. (2012) also found high AIE r values in the tropical western Pacific at Darwin (TWP) due to the low N P in December, January and February. Through studying the impact of ship-produced aerosols on the microstructure and albedo of warm marine stratocumulus clouds, Durkee et al. (2000) found that the clean and shallow boundary layers would be more readily perturbed by the addition of ship particle effluent. In this study, the higher values of AIE r and AIE N of CP-1 indicated that if the same amount of aerosol particles entered into the cloud, the microphysical parameters would be influenced more in CP-1 than in CP-2.
In addition, the meteorological conditions and the topography during the monitoring period would also affect the microphysical properties of clouds. During the observation period, CBH ranged from 460.3 to 3639.1 m with an average value of 1382.5 m. The observation station would be totally enveloped in clouds when cloud events occurred, and the corresponding distance between the observation point and CBH is represented in Fig. 2b. The sensitivity analysis of N C to CBH and v up was estimated by applying the equation as S(X i ) = ∂ ln N C /∂ ln X i , where X i represented CBH and v up . As shown in Table S2, CP-2 was more sensitive to the variation of meteorological parameters if compared with CP-1. This is consistent with the study done by McFiggans et al. (2006). They found that the sensitivity of N C to v up increased, while the sensitivity of N C to N P decreased when N P > 1000 cm −3 . The higher values of S(CBH) and S(v up ) of CP-2 indicated that CP-2 was more sensitive to the change in CBH and v up . This might cause the periodical variations of cloud microphysical properties during CP-2.

Size distribution of cloud droplets and particles
To illustrate the evolution of the aerosol particles and cloud droplets during the cloud processes, the size distributions of N P and N C during different cloud stages are plotted in Fig. 5. For each of the four size bins, ranging from 2 to 13 µm,  [150, 200) nm, and [200, 765) nm. (b) Five size ranges of N C and five size ranges of N P in SL1, SH1, SL2, SH2, and CP-2; (c) five size ranges of N C and five size ranges of N P in S1, S2, S3, S4, and NC ("NC" in c represents particle size distributions during the cloudless period). cloud number concentrations of SL1 and SL2 were lower than those of SH1 and SH2. However, in the size bin of 13-50 µm, N C of SL1 and SL2 were the largest (Fig. 5b). The size distributions of cloud droplets in SL1 and SL2 resulted in the larger r eff during the two stages, which was consistent with the result shown in Fig. 3b. During SH1 and SH2 in CP-1, the numbers of aerosol particles in all size bins increased. However, the increase in aerosol particles larger than 150 nm was the smallest, indicating that aerosols larger than 150 nm were more easily activated into cloud droplets. The activation of aerosol particles with sizes larger than 150 nm in this study dramatically increased N C of 5-10 µm and made N C of SH1 and SH2 in different size bins all comparable with those of CP-2 (Fig. 5b).
As shown in Fig. 5c, cloud droplets with D C ranging from 5 to 10 µm had high N C in each stage in CP-2, and cloud droplets with D C ranging from 13 to 50 µm had low N C in each stage if compared to CP-1. This caused the lower r eff in CP-2 than in CP-1. During CP-2, aerosol particles with diameters larger than 150 nm quickly decreased by activation when cloud events occurred, while the number of aerosol particles in the size of 50-150 nm were slightly influenced by cloud events (see the left-most graph of Fig. 5a). This was consistent with the study of Targino et al. (2007), who found aerosol size distributions of cloud residuals, which represented aerosol particles activated to cloud droplets, peaked at about 0.15 µm at Mt. Åreskutan. Mertes et al. (2005) also found that particles centred at d p = 200 nm could be efficiently activated to droplets, while most Aitken mode particles remained in the interstitial phase. Compared with other stages, S1 had the highest N C in two size bins of (2, 5) and [5,7) µm. This indicated that large numbers of cloud droplets with small sizes were formed in the beginning of cloud events in CP-2.
3.6 Relations among LWC, r eff and N C The 5 min averaged LWC for CP-1 and CP-2 is plotted against corresponding r eff in Fig. 6a. Large cloud droplets (r eff > 8 µm) were observed in CP-1, while the r eff for CP-2 varied narrowly in the range of 2.5-8 µm.
Cloud droplets with r eff > 8 µm only occurred in the two relatively clean stages, SL1 and SL2, during CP-1. It was due to the weaker competition among droplets at lower N CCN conditions. This was also observed in the US mid-Atlantic region where cloud droplets with larger sizes are more easily formed with lower N CCN (S. . At the same LWC level, the growth of cloud droplets during SH1 and SH2 was obviously limited if compared with SL1 and SL2, which is referred to as the "Twomey effect" (Twomey, 1977). This is consistent with the illustration in Fig. 3 that cloud droplets in SH1 and SH2 were smaller.
The variation r eff and/or N C can influence LWC, while the key factor may be different in different stages of the cloud. As shown in the lower row of Fig. 6a, CE-20 was taken as an example to discuss the relation among LWC, r eff and N C in different cloud stages. During S1, the existing numerous CCN (Fig. 3a) were quickly activated to form cloud droplets. The newly formed droplets were characterised with small sizes but large numbers. They will suppress the beginning of collision-coalescence processes (Rosenfeld et al., 2014a) and may further significantly delay raindrop formation (Qian et al., 2009). In S1, N C and r eff consistently varied. Both the increase in N C (from 1188 to 2940 cm −3 ) and the growth of r eff (from ∼ 3.5 µm to ∼ 4.5 µm) boosted the LWC in this stage. This is different from the results of Mazoyer et al. (2019), who found a clearly inverse relationship between the number and the size of droplets at the beginning of the first hour of fog events during an observation in suburban Paris. When compared with fog, cloud is usually formed under conditions with more condensible water vapour (Fig. 1). The limited growth of droplets in fog will not occur in cloud. This causes the positive relationship with cloud droplet number and droplet size. At the beginning of S2, N C reached the maximum. The high N C yielded a great coalescence rate between cloud droplets. Meanwhile, the coalescence process was self-accelerating (Freud and Rosenfeld, 2012) and thus caused the quick decrease in N C (Fig. 3a). This made cloud droplets in S2 that were characterised by larger sizes and lower number concentrations, whilst LWC simply varied in a relatively narrow range (Fig. 6a). During S3, N C was almost constant due to the formation, coagulation and evaporation of the cloud droplets reaching a balance. As shown in Fig. 6a, the relationship between r eff and LWC in this stage could be fitting as r eff = a × LWC 0.34±0.02 , which meant under the increase in LWC, the N C was almost unchanged. The variation of LWC values was mainly due to the changes in droplet sizes. At the dissipation stage of S4, the increase in CBH brought air with low RH and high N P to the summit of Mt. Tai and caused the dissipation of cloud events (Figs. 2c and 3a). The previously activated CCN returned to the interstitial aerosol phase due to the evaporation of the droplets (Verheggen et al., 2007). Both N C and r eff declined. It was also illustrated in Fig. 5c that all of the N C in the five size bins of cloud droplets decreased in S4.
In order to investigate the variation of r eff upon N C , the distribution of r eff was classified with different N C ranges in Fig. 6b. For N C < 1000 cm −3 , r eff displayed a trimodal distribution and concentrated on 3.25 µm (Peak-1), 4.86 µm (Peak-2) and 7.52 µm (Peak-3). Peak-1 corresponded to cloud droplets with low N C , LWC and r eff values, while the N CCN0.2 was very high (Fig. 6c). These points represented cloud droplets in the incipient stage or the dissipation stage of cloud events where large numbers of CCN exist in the atmosphere. Peak-2 and Peak-3 represented the matured stages for cloud events with different environmental conditions. Peak-3 represented cloud droplets formed under a relatively clean atmosphere. In this circumstance, CCN were efficiently activated and had a lower concentration remaining in the atmosphere (Fig. 6c). The sufficient ambient water vapour ac- celerated the growth of the formed droplets, which were characterised with low N C and LWC but large r eff . Peak-2 represented cloud droplets formed under relatively polluted conditions and was the only peak found for N C larger than 1000 cm −3 . With the increase in N C , the distribution of this peak narrowed and slightly moved to lower r eff mode.
The thickness of orographic cloud can be easily influenced by the specific topography and environmental conditions (Barros and Lettenmaier, 1994;Welch et al., 2008). When assuming the cloud thickness during CP-1 and CP-2 to be equal, albedo would depend on the values of LWC and N C , as described in Sect. 2.8. Cloud albedo during CP-2 was always higher than that during CP-1, especially when the cloud thickness was lower than about 2500 m (Fig. 6d).
Note that the increase in N C could enhance the evaporation and further reduce the lifetime of cloud, which was not taken into account when calculating the induced albedo. Through studying marine stratocumulus clouds in the north-eastern Pacific Ocean, Twohy et al. (2005) also found that the increase in N C by a factor of 2.8 would lead to 40 % increase in albedo, going from 0.325 to 0.458. This indicated that the higher N C would increase the cloud albedo when assuming no change to cloud thickness and when disregarding any cloud lifetime effects.

Conclusions
From 17 June to 30 July 2018, in situ observations of number concentrations and size distributions of aerosol particles and cloud droplets were used to study aerosol-cloud interactions at the summit of Mt. Tai. Large variations of the characteristic values in terms of N C , LWC and r eff were found during the observation period. Cloud events with small r eff and low LWC similar to urban fog were also observed at Mt. Tai.
Two typical cloud processes, CP-1 and CP-2, were applied to study the cloud-aerosol interactions based on the aerosol characteristics (especially N P and N CCN ) before cloud onsets. For the CP-1, which corresponded to relatively clean conditions, water content was sufficient, while N CCN was considered the limitation of cloud droplet formation according to the observation results. The newly formed cloud droplets were characterised with low N C , high LWC / N C and large r eff . With the increase in aerosol concentration, N C was found to dramatically increase. Large numbers of N CCN competed for the water content with the formed cloud droplets and, as a result, further dramatically decrease the LWC / N C and r eff values of cloud droplets. In CP-2, N P before the cloud onset was high and N CCN was sufficient. Water vapour was considered the limitation for cloud formation. Large numbers of small cloud droplets with low LWC / N C were observed in the incipient stage of cloud events. In addition, periodical changes in cloud microphysical properties were found. Both consistent variation and inverse variation between N P and N C were observed in this study, which were characterised with relatively high and low LWC / N C values, respectively.
Both positive AIE r and AIE N values observed at Mt. Tai indicate that the increase in N P will decrease r eff and increase N C of cloud droplets. The lower values of AIE r and AIE N with higher N P and N CCN suggest that the increase in N P will more strongly decrease the size and increase the number of cloud droplets under the conditions with smaller aerosol amounts. Through studying the size distributions of aerosol particles and cloud droplets, higher N C in the size bin of 13-50 µm resulted in the larger r eff during the two clean stages in CP-1. Particles larger than 150 nm were able to be efficiently activated to cloud droplets and made important contributions to the increase in cloud droplets in the size range of 5-10 µm.
The LWC of cloud depended on the change in r eff and N C . However, the decisive factor may differ at different stages of the cloud. In general, the r eff of cloud droplets consistently varied with LWC. However, in different N C ranges, the r eff of cloud droplets were observed with different distribution shapes. For N C < 1000 cm −3 , r eff displayed a trimodal distribution. Three peaks at 3.25, 4.86 and 7.52 µm were observed. With the increase in N C , a narrowed unimodal distribution of r eff appeared, and the peak value slightly moved towards lower r eff mode. For a constant cloud thickness and disregarding any cloud lifetime effects, the increased N C and decreased r eff dramatically increase the cloud albedo, which may further influence the regional climate in the North China Plain.
The local topography of the surrounding areas at Mt. Tai supplies potential access for aerosol transportation and can affect the measured cloud droplet distributions by increasing turbulence or causing orographic flows. Even though the summit of Mt. Tai is far away from the polluted sources, the transported CCN could change the cloud microphysical properties (e.g. during CP-1). The cloud microphysical parameters derived in our study characterised the cloud features in the North China Plain and provided valuable data for modelling studies of cloud microphysics in the future.
Data availability. All data used to support the conclusion are presented in this paper. Additional data are available upon request. Please contact the corresponding authors Jianmin Chen (jm-chen@fudan.edu.cn) and Hui Chen (hui_chen@fudan.edu.cn).
Author contributions. JC and HC conceived the study. JL and CZ performed the field experiments and sampled cloud water. JL analysed the data and wrote the main manuscript text. JC, HC, DZ, CZ and HH revised the initial manuscript. LX, XW and HL supported the meteorological data and PM 2.5 mass concentration. PL, JL, CZ, YM and WZ assisted in instrument maintenance. LZ, KL and ML contributed to the organisation and arrangement of the field observations. LZ provided the meteorological parameters of Tai'an City. All of the authors discussed the results and contributed to the final manuscript.
Competing interests. The authors declare that they have no conflict of interest.
Review statement. This paper was edited by Ilona Riipinen and reviewed by Kevin Noone and three anonymous referees.