2015, Kupiszewski, P., Weingartner, E., Vochezer, P., Schnaiter, M., Bigi, A., Gysel, M., Rosati, B., Toprak, E., Mertes, S., Baltensperger, U.

Climate predictions are affected by high uncertainties partially due to an insufficient knowledge of aerosol–cloud interactions. One of the poorly understood processes is formation of mixed-phase clouds (MPCs) via heterogeneous ice nucleation. Field measurements of the atmospheric ice phase in MPCs are challenging due to the presence of much more numerous liquid droplets. The Ice Selective Inlet (ISI), presented in this paper, is a novel inlet designed to selectively sample pristine ice crystals in mixed-phase clouds and extract the ice residual particles contained within the crystals for physical and chemical characterization. Using a modular setup composed of a cyclone impactor, droplet evaporation unit and pumped counterflow virtual impactor (PCVI), the ISI segregates particles based on their inertia and phase, exclusively extracting small ice particles between 5 and 20 μm in diameter. The setup also includes optical particle spectrometers for analysis of the number size distribution and shape of the sampled hydrometeors. The novelty of the ISI is a droplet evaporation unit, which separates liquid droplets and ice crystals in the airborne state, thus avoiding physical impaction of the hydrometeors and limiting potential artefacts. The design and validation of the droplet evaporation unit is based on modelling studies of droplet evaporation rates and computational fluid dynamics simulations of gas and particle flows through the unit. Prior to deployment in the field, an inter-comparison of the optical particle size spectrometers and a characterization of the transmission efficiency of the PCVI was conducted in the laboratory. The ISI was subsequently deployed during the Cloud and Aerosol Characterization Experiment (CLACE) 2013 and 2014 – two extensive international field campaigns encompassing comprehensive measurements of cloud microphysics, as well as bulk aerosol, ice residual and ice nuclei properties. The campaigns provided an important opportunity for a proof of concept of the inlet design. In this work we present the setup of the ISI, including the modelling and laboratory characterization of its components, as well as field measurements demonstrating the ISI performance and validating the working principle of the inlet. Finally, measurements of biological aerosol during a Saharan dust event (SDE) are presented, showing a first indication of enrichment of bio-material in sub-2 μm ice residuals.

2017, Schladitz, A., Müller, T., Kaaden, N., Massling, A., Kandler, K., Ebert, M., Weinbruch, S., Deutscher, C., Wiedensohler, A.

In situ measurements of optical and physical properties of mineral dust were performed at the outskirts of the Saharan Desert in the framework of the Saharan Mineral Dust Experiment part 1 (SAMUM-1). Goals of the field study were to achieve information on the extent and composition of the dust particle size distribution and the optical properties of dust at the ground. For the particle number size distribution, measured with a DMPS/APS, a size dependent dynamic shape factor was considered. The mean refractive index of the particles in this field study is 1.53–4.1 × 10-3i at 537 nm wavelength and 1.53–3.1 × 10-3i at 637 nm wavelength derived from measurements of scattering and absorption coefficients, as well as the particle size distribution. Whereas the real part of the refractive index is rather constant, the imaginary part varies depending on the mineral dust concentrations. For high dust concentration the single scattering albedo is primarily influenced by iron oxide and is 0.96 ± 0.02 and 0.98 ± 0.01 at 537 nm and 637 nm wavelength, respectively. During low dust concentration the single scattering albedo is more influenced by a soot-type absorber and is 0.89 ± 0.02 and 0.93 ± 0.01 for the same wavelengths.

2021, Bousiotis, Dimitrios, Pope, Francis D., Beddows, David C. S., Dall'Osto, Manuel, Massling, Andreas, Nøjgaard, Jakob Klenø, Nordstrøm, Claus, Niemi, Jarkko V., Portin, Harri, Petäjä, Tuukka, Perez, Noemi, Alastuey, Andrés, Querol, Xavier, Kouvarakis, Giorgos, Mihalopoulos, Nikos, Vratolis, Stergios, Eleftheriadis, Konstantinos, Wiedensohler, Alfred, Weinhold, Kay, Merkel, Maik, Tuch, Thomas, Harrison, Roy M.

New particle formation (NPF) events occur almost everywhere in the world and can play an important role as a particle source. The frequency and characteristics of NPF events vary spatially, and this variability is yet to be fully understood. In the present study, long-term particle size distribution datasets (minimum of 3 years) from 13 sites of various land uses and climates from across Europe were studied, and NPF events, deriving from secondary formation and not traffic-related nucleation, were extracted and analysed. The frequency of NPF events was consistently found to be higher at rural background sites, while the growth and formation rates of newly formed particles were higher at roadsides (though in many cases differences between the sites were small), underlining the importance of the abundance of condensable compounds of anthropogenic origin found there. The growth rate was higher in summer at all rural background sites studied. The urban background sites presented the highest uncertainty due to greater variability compared to the other two types of site. The origin of incoming air masses and the specific conditions associated with them greatly affect the characteristics of NPF events. In general, cleaner air masses present higher probability for NPF events, while the more polluted ones show higher growth rates. However, different patterns of NPF events were found, even at sites in close proximity (<ĝ€¯200ĝ€¯km), due to the different local conditions at each site. Region-wide events were also studied and were found to be associated with the same conditions as local events, although some variability was found which was associated with the different seasonality of the events at two neighbouring sites. NPF events were responsible for an increase in the number concentration of ultrafine particles of more than 400ĝ€¯% at rural background sites on the day of their occurrence. The degree of enhancement was less at urban sites due to the increased contribution of other sources within the urban environment. It is evident that, while some variables (such as solar radiation intensity, relative humidity, or the concentrations of specific pollutants) appear to have a similar influence on NPF events across all sites, it is impossible to predict the characteristics of NPF events at a site using just these variables, due to the crucial role of local conditions. © Author(s) 2021.

2018, DeMott, Paul J., Möhler, Ottmar, Cziczo, Daniel J., Hiranuma, Naruki, Petters, Markus D., Petters, Sarah S., Belosi, Franco, Bingemer, Heinz G., Brooks, Sarah D., Budke, Carsten, Burkert-Kohn, Monika, Collier, Kristen N., Danielczok, Anja, Eppers, Oliver, Felgitsch, Laura, Garimella, Sarvesh, Grothe, Hinrich, Herenz, Paul, Hill, Thomas C. J., Höhler, Kristina, Kanji, Zamin A., Kiselev, Alexei, Koop, Thomas, Kristensen, Thomas B., Krüger, Konstantin, Kulkarni, Gourihar, Levin, Ezra J. T., Murray, Benjamin J., Nicosia, Alessia, O'Sullivan, Daniel, Peckhaus, Andreas, Polen, Michael J., Price, Hannah C., Reicher, Naama, Rothenberg, Daniel A., Rudich, Yinon, Santachiara, Gianni, Schiebel, Thea, Schrod, Jann, Seifried, Teresa M., Stratmann, Frank, Sullivan, Ryan C., Suski, Kaitlyn J., Szakáll, Miklós, Taylor, Hans P., Ullrich, Romy, Vergara-Temprado, Jesus, Wagner, Robert, Whale, Thomas F., Weber, Daniel, Welti, André, Wilson, Theodore W., Wolf, Martin J., Zenker, Jake

The second phase of the Fifth International Ice Nucleation Workshop (FIN-02) involved the gathering of a large number of researchers at the Karlsruhe Institute of Technology's Aerosol Interactions and Dynamics of the Atmosphere (AIDA) facility to promote characterization and understanding of ice nucleation measurements made by a variety of methods used worldwide. Compared to the previous workshop in 2007, participation was doubled, reflecting a vibrant research area. Experimental methods involved sampling of aerosol particles by direct processing ice nucleation measuring systems from the same volume of air in separate experiments using different ice nucleating particle (INP) types, and collections of aerosol particle samples onto filters or into liquid for sharing amongst measurement techniques that post-process these samples. In this manner, any errors introduced by differences in generation methods when samples are shared across laboratories were mitigated. Furthermore, as much as possible, aerosol particle size distribution was controlled so that the size limitations of different methods were minimized. The results presented here use data from the workshop to assess the comparability of immersion freezing measurement methods activating INPs in bulk suspensions, methods that activate INPs in condensation and/or immersion freezing modes as single particles on a substrate, continuous flow diffusion chambers (CFDCs) directly sampling and processing particles well above water saturation to maximize immersion and subsequent freezing of aerosol particles, and expansion cloud chamber simulations in which liquid cloud droplets were first activated on aerosol particles prior to freezing. The AIDA expansion chamber measurements are expected to be the closest representation to INP activation in atmospheric cloud parcels in these comparisons, due to exposing particles freely to adiabatic cooling. The different particle types used as INPs included the minerals illite NX and potassium feldspar (K-feldspar), two natural soil dusts representative of arable sandy loam (Argentina) and highly erodible sandy dryland (Tunisia) soils, respectively, and a bacterial INP (Snomax®). Considered together, the agreement among post-processed immersion freezing measurements of the numbers and fractions of particles active at different temperatures following bulk collection of particles into liquid was excellent, with possible temperature uncertainties inferred to be a key factor in determining INP uncertainties. Collection onto filters for rinsing versus directly into liquid in impingers made little difference. For methods that activated collected single particles on a substrate at a controlled humidity at or above water saturation, agreement with immersion freezing methods was good in most cases, but was biased low in a few others for reasons that have not been resolved, but could relate to water vapor competition effects. Amongst CFDC-style instruments, various factors requiring (variable) higher supersaturations to achieve equivalent immersion freezing activation dominate the uncertainty between these measurements, and for comparison with bulk immersion freezing methods. When operated above water saturation to include assessment of immersion freezing, CFDC measurements often measured at or above the upper bound of immersion freezing device measurements, but often underestimated INP concentration in comparison to an immersion freezing method that first activates all particles into liquid droplets prior to cooling (the PIMCA-PINC device, or Portable Immersion Mode Cooling chAmber-Portable Ice Nucleation Chamber), and typically slightly underestimated INP number concentrations in comparison to cloud parcel expansions in the AIDA chamber; this can be largely mitigated when it is possible to raise the relative humidity to sufficiently high values in the CFDCs, although this is not always possible operationally. Correspondence of measurements of INPs among direct sampling and post-processing systems varied depending on the INP type. Agreement was best for Snomax® particles in the temperature regime colder than -10°C, where their ice nucleation activity is nearly maximized and changes very little with temperature. At temperatures warmer than -10°C, Snomax® INP measurements (all via freezing of suspensions) demonstrated discrepancies consistent with previous reports of the instability of its protein aggregates that appear to make it less suitable as a calibration INP at these temperatures. For Argentinian soil dust particles, there was excellent agreement across all measurement methods; measures ranged within 1 order of magnitude for INP number concentrations, active fractions and calculated active site densities over a 25 to 30°C range and 5 to 8 orders of corresponding magnitude change in number concentrations. This was also the case for all temperatures warmer than -25°C in Tunisian dust experiments. In contrast, discrepancies in measurements of INP concentrations or active site densities that exceeded 2 orders of magnitude across a broad range of temperature measurements found at temperatures warmer than -25°C in a previous study were replicated for illite NX. Discrepancies also exceeded 2 orders of magnitude at temperatures of -20 to -25°C for potassium feldspar (K-feldspar), but these coincided with the range of temperatures at which INP concentrations increase rapidly at approximately an order of magnitude per 2°C cooling for K-feldspar. These few discrepancies did not outweigh the overall positive outcomes of the workshop activity, nor the future utility of this data set or future similar efforts for resolving remaining measurement issues. Measurements of the same materials were repeatable over the time of the workshop and demonstrated strong consistency with prior studies, as reflected by agreement of data broadly with parameterizations of different specific or general (e.g., soil dust) aerosol types. The divergent measurements of the INP activity of illite NX by direct versus post-processing methods were not repeated for other particle types, and the Snomax° data demonstrated that, at least for a biological INP type, there is no expected measurement bias between bulk collection and direct immediately processed freezing methods to as warm as -10°C. Since particle size ranges were limited for this workshop, it can be expected that for atmospheric populations of INPs, measurement discrepancies will appear due to the different capabilities of methods for sampling the full aerosol size distribution, or due to limitations on achieving sufficient water supersaturations to fully capture immersion freezing in direct processing instruments. Overall, this workshop presents an improved picture of present capabilities for measuring INPs than in past workshops, and provides direction toward addressing remaining measurement issues.

2016, Brock, Charles A., Wagner, Nicholas L., Anderson, Bruce E., Attwood, Alexis R., Beyersdorf, Andreas, Campuzano-Jost, Pedro, Carlton, Annmarie G., Day, Douglas A., Diskin, Glenn S., Gordon, Timothy D., Jimenez, Jose L., Lack, Daniel A., Liao, Jin, Markovic, Milos Z., Middlebrook, Ann M., Ng, Nga L., Perring, Anne E., Richardson, Matthews S., Schwarz, Joshua P., Washenfelder, Rebecca A., Welti, Andre, Xu, Lu, Ziemba, Luke D., Murphy, Daniel M.

Aircraft observations of meteorological, trace gas, and aerosol properties were made between May and September 2013 in the southeastern United States (US). Regionally representative aggregate vertical profiles of median and interdecile ranges of the measured parameters were constructed from 37 individual aircraft profiles made in the afternoon when a well-mixed boundary layer with typical fair-weather cumulus was present (Wagner et al., 2015). We use these 0–4 km aggregate profiles and a simple model to calculate the sensitivity of aerosol optical depth (AOD) to changes in dry aerosol mass, relative humidity, mixed-layer height, the central diameter and width of the particle size distribution, hygroscopicity, and dry and wet refractive index, while holding the other parameters constant. The calculated sensitivity is a result of both the intrinsic sensitivity and the observed range of variation in these parameters. These observationally based sensitivity studies indicate that the relationship between AOD and dry aerosol mass in these conditions in the southeastern US can be highly variable and is especially sensitive to relative humidity (RH). For example, calculated AOD ranged from 0.137 to 0.305 as the RH was varied between the 10th and 90th percentile profiles with dry aerosol mass held constant. Calculated AOD was somewhat less sensitive to aerosol hygroscopicity, mean size, and geometric standard deviation, σg. However, some chemistry–climate models prescribe values of σg substantially larger than we or others observe, leading to potential high biases in model-calculated AOD of  ∼  25 %. Finally, AOD was least sensitive to observed variations in dry and wet aerosol refractive index and to changes in the height of the well-mixed surface layer. We expect these findings to be applicable to other moderately polluted and background continental air masses in which an accumulation mode between 0.1–0.5 µm diameter dominates aerosol extinction.

2012, Wiedensohler, A., Birmili, W., Nowak, A., Sonntag, A., Weinhold, K., Merkel, M., Wehner, B., Tuch, T., Pfeifer, S., Fiebig, M., Fjäraa, A.M., Asmi, E., Sellegri, K., Depuy, R., Venzac, H., Villani, P., Laj, P., Aalto, P., Ogren, J.A., Swietlick, E., Williams, P., Roldin, P., Quincey, P., Hüglin, C., Fierz-Schmidhauser, R., Gysel, M., Weingartner, E., Riccobono, F., Santos, S., Grüning, C., Faloon, K., Beddows, D., Harrison, R., Monahan, C., Jennings, S.G., O'Dowd, C.D., Marinoni, A., Horn, H.-G., Keck, L., Jiang, J., Scheckman, J., McMurry, P.H., Deng, Z., Zhao, C.S., Moerman, M., Henzing, B., de Leeuw, G., Löschau, G., Bastian, S.

Mobility particle size spectrometers often referred to as DMPS (Differential Mobility Particle Sizers) or SMPS (Scanning Mobility Particle Sizers) have found a wide range of applications in atmospheric aerosol research. However, comparability of measurements conducted world-wide is hampered by lack of generally accepted technical standards and guidelines with respect to the instrumental set-up, measurement mode, data evaluation as well as quality control. Technical standards were developed for a minimum requirement of mobility size spectrometry to perform long-term atmospheric aerosol measurements. Technical recommendations include continuous monitoring of flow rates, temperature, pressure, and relative humidity for the sheath and sample air in the differential mobility analyzer. We compared commercial and custom-made inversion routines to calculate the particle number size distributions from the measured electrical mobility distribution. All inversion routines are comparable within few per cent uncertainty for a given set of raw data. Furthermore, this work summarizes the results from several instrument intercomparison workshops conducted within the European infrastructure project EUSAAR (European Supersites for Atmospheric Aerosol Research) and ACTRIS (Aerosols, Clouds, and Trace gases Research InfraStructure Network) to determine present uncertainties especially of custom-built mobility particle size spectrometers. Under controlled laboratory conditions, the particle number size distributions from 20 to 200 nm determined by mobility particle size spectrometers of different design are within an uncertainty range of around ±10% after correcting internal particle losses, while below and above this size range the discrepancies increased. For particles larger than 200 nm, the uncertainty range increased to 30%, which could not be explained. The network reference mobility spectrometers with identical design agreed within ±4% in the peak particle number concentration when all settings were done carefully. The consistency of these reference instruments to the total particle number concentration was demonstrated to be less than 5%. Additionally, a new data structure for particle number size distributions was introduced to store and disseminate the data at EMEP (European Monitoring and Evaluation Program). This structure contains three levels: raw data, processed data, and final particle size distributions. Importantly, we recommend reporting raw measurements including all relevant instrument parameters as well as a complete documentation on all data transformation and correction steps. These technical and data structure standards aim to enhance the quality of long-term size distribution measurements, their comparability between different networks and sites, and their transparency and traceability back to raw data.

2011, Sorribas, M., de la Morena, B.A., Wehner, B., López, J.F., Prats, N., Mogo, S., Wiedensohler, A., Cachorro, V.E.

This study focuses on the analysis of the sub-micron aerosol characteristics at El Arenosillo Station, a rural and coastal environment in South-western Spain between 1 August 2004 and 31 July 2006 (594 days). The mean total concentration (NT) was 8660 cm−3 and the mean concentrations in the nucleation (NNUC), Aitken (NAIT) and accumulation (NACC) particle size ranges were 2830 cm−3, 4110 cm−3 and 1720 cm−3, respectively. Median size distribution was characterised by a single-modal fit, with a geometric diameter, median number concentration and geometric standard deviation of 60 nm, 5390 cm−3 and 2.31, respectively. Characterisation of primary emissions, secondary particle formation, changes to meteorology and long-term transport has been necessary to understand the seasonal and annual variability of the total and modal particle concentration. Number concentrations exhibited a diurnal pattern with maximum concentrations around noon. This was governed by the concentrations of the nucleation and Aitken modes during the warm seasons and only by the nucleation mode during the cold seasons. Similar monthly mean total concentrations were observed throughout the year due to a clear inverse variation between the monthly mean NNUC and NACC. It was related to the impact of desert dust and continental air masses on the monthly mean particle levels. These air masses were associated with high values of NACC which suppressed the new particle formation (decreasing NNUC). Each day was classified according to a land breeze flow or a synoptic pattern influence. The median size distribution for desert dust and continental aerosol was dominated by the Aitken and accumulation modes, and marine air masses were dominated by the nucleation and Aitken modes. Particles moved offshore due to the land breeze and had an impact on the particle burden at noon, especially when the wind was blowing from the NW sector in the morning during summer time. This increased NNUC and NAIT by factors of 3.1 and 2.4, respectively. Nucleation events with the typical "banana" shape were characterised by a mean particle nucleation rate of 0.74 cm−3 s−1, a mean growth rate of 1.96 nm h−1 and a mean total duration of 9.25 h (starting at 10:55 GMT and ending at 20:10 GMT). They were observed for 48 days. Other nucleation events were identified as those produced by the emissions from the industrial areas located at a distance of 35 km. They were observed for 42 days. Both nucleation events were strongly linked to the marine air mass origin.

2014, Pfeifer, S., Birmili, W., Schladitz, A., Müller, T., Nowak, A., Wiedensohler, A.

Multiple-charge inversion is an essential procedure to convert the raw mobility distributions recorded by mobility particle size spectrometers, such as the DMPS or SMPS (differential or scanning mobility particle sizers), into true particle number size distributions. In this work, we present a fast and easy-to-implement multiple-charge inversion algorithm with sufficient precision for atmospheric conditions, but extended functionality. The algorithm can incorporate size distribution information from sensors that measure beyond the upper sizing limit of the mobility spectrometer, such as an aerodynamic particle sizer (APS) or an optical particle counter (OPC). This feature can considerably improve the multiple-charge inversion result in the upper size range of the mobility spectrometer, for example, when substantial numbers of coarse particles are present. The program also yields a continuous size distribution from both sensors as an output. The algorithm is able to calculate the propagation of measurement errors, such as those based on counting statistics, into on the final particle number size distribution. As an additional aspect, the algorithm can perform all inversion steps under the assumption of non-spherical particle shape, including constant or size-dependent shape factors.

2012, Aas, W., Tsyro, S., Bieber, E., Bergström, R., Ceburnis, D., Ellermann, T., Fagerli, H., Frölich, M., Gehrig, R., Makkonen, U., Nemitz, E., Otjes, R., Perez, N., Perrino, C., Prévôt, A.S.H., Putaud, J.-P., Simpson, D., Spindler, G., Vana, M., Yttri, K.E.

The first EMEP intensive measurement periods were held in June 2006 and January 2007. The measurements aimed to characterize the aerosol chemical compositions, including the gas/aerosol partitioning of inorganic compounds. The measurement program during these periods included daily or hourly measurements of the secondary inorganic components, with additional measurements of elemental- and organic carbon (EC and OC) and mineral dust in PM1, PM2.5 and PM10. These measurements have provided extended knowledge regarding the composition of particulate matter and the temporal and spatial variability of PM, as well as an extended database for the assessment of chemical transport models. This paper summarise the first experiences of making use of measurements from the first EMEP intensive measurement periods along with EMEP model results from the updated model version to characterise aerosol composition. We investigated how the PM chemical composition varies between the summer and the winter month and geographically. The observation and model data are in general agreement regarding the main features of PM10 and PM2.5 composition and the relative contribution of different components, though the EMEP model tends to give slightly lower estimates of PM10 and PM2.5 compared to measurements. The intensive measurement data has identified areas where improvements are needed. Hourly concurrent measurements of gaseous and particulate components for the first time facilitated testing of modelled diurnal variability of the gas/aerosol partitioning of nitrogen species. In general, the modelled diurnal cycles of nitrate and ammonium aerosols are in fair agreement with the measurements, but the diurnal variability of ammonia is not well captured. The largest differences between model and observations of aerosol mass are seen in Italy during winter, which to a large extent may be explained by an underestimation of residential wood burning sources. It should be noted that both primary and secondary OC has been included in the calculations for the first time, showing promising results. Mineral dust is important, especially in southern Europe, and the model seems to capture the dust episodes well. The lack of measurements of mineral dust hampers the possibility for model evaluation for this highly uncertain PM component. There are also lessons learnt regarding improved measurements for future intensive periods. There is a need for increased comparability between the measurements at different sites. For the nitrogen compounds it is clear that more measurements using artefact free methods based on continuous measurement methods and/or denuders are needed. For EC/OC, a reference methodology (both in field and laboratory) was lacking during these periods giving problems with comparability, though measurement protocols have recently been established and these should be followed by the Parties to the EMEP Protocol. For measurements with no defined protocols, it might be a good solution to use centralised laboratories to ensure comparability across the network. To cope with the introduction of these new measurements, new reporting guidelines have been developed to ensure that all proper information about the methodologies and data quality is given.

2013, Wang, Z.B., Hu, M., Sun, J.Y., Wu, Z.J., Yue, D.L., Shen, X.J., Zhang, Y.M., Pei, X.Y., Cheng, Y.F., Wiedensohler, A.

Long-term measurements of particle number size distributions were carried out both at an urban background site (Peking University, PKU) and a regional Global Atmospheric Watch station (Shangdianzi, SDZ) from March to November in 2008. In total, 52 new particle formation (NPF) events were observed simultaneously at both sites, indicating that this is a regional phenomenon in the North China Plain. On average, the mean condensation sink value before the nucleation events started was 0.025 s−1 in the urban environment, which was 1.6 times higher than that at regional site. However, higher particle formation and growth rates were observed at PKU (10.8 cm−3 s−1 and 5.2 nm h−1) compared with those at SDZ (4.9 cm−3 s−1 and 4.0 nm h−1). These results implied that precursors were much more abundant in the polluted urban environment. Different from the observations in cleaner environments, the background conditions of the observed particle homogeneous nucleation events in the North China Plain could be characterized as the co-existing of a stronger source of precursor gases and a higher condensational sink of pre-existing aerosol particles. Secondary aerosol formation following nucleation events results in an increase of particle mass concentration, particle light scattering coefficient, and cloud condensation nuclei (CCN) number concentration, with consequences on visibility, radiative effects, and air quality. Typical regional NPF events with significant particle nucleation rates and subsequent particle growth over a sufficiently long time period at both sites were chosen to investigate the influence of NPF on the number concentration of "potential" CCN. As a result, the NPF and the subsequent condensable growth increased the CCN number concentration in the North China Plain by factors in the range from 5.6 to 8.7. Moreover, the potential contribution of anthropogenic emissions to the CCN number concentration was more than 50%, to which more attention should be drawn in regional and global climate modeling, especially in the polluted urban areas.