Next Article in Journal
Impact of the Eclipsed Sun on Terrestrial Atmospheric Parameters in Desert Locations: A Comprehensive Overview and Two Events Case Study in Saudi Arabia
Previous Article in Journal
Multi-Scale Reconstruction of Turbulent Rotating Flows with Generative Diffusion Models
Previous Article in Special Issue
Comparison of 24 h Surface Ozone Forecast for Poland: CAMS Models vs. Simple Statistical Models with Limited Number of Input Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Hygroscopic Aerosol Growth Based on the OPC-N3 Counter

by
Katarzyna Nurowska
* and
Krzysztof M. Markowicz
Institute of Geophysics, Faculty of Physics, University of Warsaw, 00-927 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Atmosphere 2024, 15(1), 61; https://doi.org/10.3390/atmos15010061
Submission received: 14 November 2023 / Revised: 12 December 2023 / Accepted: 23 December 2023 / Published: 31 December 2023
(This article belongs to the Special Issue Aerosol Pollution in Central Europe)

Abstract

:
We present the first estimations of single hygroscopic parameter κ for Polish urban area. The results were obtained using an inexpensive OPC-N3 optical particulate matter counter and the ASC 1000 Acoem Aerosol Conditioning System. Our studies were carried out during the winter and spring seasons, between 2020 and 2022 in Warsaw, Poland. We study the difference of κ between these two seasons, as the aerosol has different origins. The mean value of the hygroscopicity parameter for spring (April–June) was 0.13 ± 0.13 and for winter (December-March) it was 0.04 ± 0.04 (obtained using PM 1 ). The mean value of the growth factor for spring (April–June) was 1.52 ± 0.23 and for winter (December–March) it was 1.16 ± 0.13 (obtained using PM 1 ). The values for winter period suggest that during the cold season the aerosol mixture in the atmosphere is dominated by non-hygroscopic particles. However events of higher κ > 0.5 where occurring (mostly for big aerosol particles 1–10 μ m in size), when the air mass was coming from North Atlantic carrying sea salt particles. Furthermore, based on κ we propose a method to remove the dependence of PM values on relative humidity on the OPC-N3 optical particulate counter.

1. Introduction

Aerosol particles directly and indirectly affect the Earth’s climate. The Intergovernmental Panel on Climate Change (IPCC) [1] assigned a great uncertainty to this influence. Aerosols scatters and absorb radiation, some act as the onset of cloud condensation nuclei furthermore they interacts with cloud droplets. The hygroscopicity is a key factor which influences those processes.
Particulate matter is one of the wildly monitored atmospheric compounds in the world. The high mass concentration of PM in the atmosphere can be associated with health problems [2,3]. For this reason several countries introduced limits of concentration for PM diameter lower than 2.5 μ m (PM 2.5 ) and 10 μ m (PM 10 ) [4,5]. The amount of PM changes rapidly in time and space. Polish cities are often exposed to values above the norm [6]. The concentration of PM can be different in distinct parts of the city (the reason for this may be the amount of traffic, the way of heating houses, and land use). For this reason, more dense measurement networks are created. To lower the cost of such systems, the low-cost apparatus is needed. The low-cost optical particle counters, like OPC-N3, can elevate the PM values in high RH. This can cause alerts of poor air conditions despite the fact that the actual level of PM is normal. For such reasons PM counters should operate with dehumidifiers or should be corrected on the RH of surrounding air. The OPC-N3 measures PM 1 , PM 2.5 and PM 10 . The measured PM values significantly relay on the air’s relative humidity.
When relative humidity (RH) is high, aerosols in the air have reservoir of water, therefore can grow larger and scatter more light. Parameters such as refractive index, angular scattering, size, or shape are sensitive to dousing aerosol particles with water. Typically, to remove the effects induced by hygroscopic growth the optical properties of the aerosol such as the scattering coefficient are obtained in low relative humidity (RH < 40%) [7]. The hygroscopicity of an aerosol can be characterized in several different ways, e.g., by the change of radius between wet and dry particles, by the mass ratio of wet and dry particles, or by the change of the light scattering coefficient between wet and dry particles. To assess the size of the hygroscopic effect one can calculate an enhancement factor, where the scattering coefficient in a given RH is divided by the scattering coefficient in dry conditions. Most scattering coefficient measurements are carried out in low RH conditions; however, for the proper operation of climate models, it is necessary to know the changes occurring when water begins to condense on aerosols. In fact, in the Arctic the enhancement effects is pronounced due to the presence of more sea salt in the air. Zieger et al. [8] demonstrate that the mean scattering enhancement factor in this region is 3.2.
Assuming incorrect values of growth factors in models can lead to biasing of the results. Zieger et al. [9] shows that the hygroscopicity of sea salt is different from the hygroscopicity of sodium chloride (many works assumed the hygroscopicity is the same). Furthermore correcting the κ parameter which describe the hygroscopicity in a global climate model provide differences in aerosol optical depth (AOD) in some areas by almost 20 % [9].
The Köhler theory describes how the particle will grow exposed to a certain RH level. To calculate the activity of cloud condensation nuclei (CCN) based on this theory, one needs to know the solute mass, dissociable ions, and other physicochemical properties. Petters and Kreidenweis [10] proposed a single dimensionless parameter κ to describe the hygroscopicity. Introduction of this parameter allowed to assign one parameter to a mixture of different aerosols or assigning a single value to molecules consisting of insoluble core immersed in aqueous solution.
Petters and Kreidenweis [10] calculated the values of κ for atmospheric aerosols. For typical aerosols found in the atmosphere values lie between 0.1 < κ < 0.9 . The κ = 0 denotes substances non-hygroscopic, while a highly higroscopic aerosol is often found in the marine atmosphere (sodium chloride reaches the value κ = 1.5 ). The article by Petters and Kreidenweis [10] neglected soot, which can be found in polluted urban areas. The values of κ for soot range between 0 and 10 3 [11]. The calculation of κ is possible by fitting nonlinear curve with one parameter to the relation between growth factor and air RH.
Aerosol present in the atmosphere is a mixture of particles with different chemical composition. In this case, one parameter κ is used to describe the resulting hygroscopicity of the mixture of different particles. The chemical compound of the aerosol is also crucial to describe its behavior under various RH conditions. Under the influence of high RH, some particles may be covered with water (adsorption) or dissolve in water (absorption), causing further changes in their properties [12].
Aerosols display various behaviors on the basis of their interaction with environmental conditions. Some aerosols remain insoluble in water and maintain a relatively stable size even when exposed to higher RH levels, such as soot. Others, such as sulfuric acid (H 2 SO 4 ), are soluble and undergo diameter changes when exposed to varying levels of RH. There are also aerosols that expand as humidity increases, yet they have specific markers in a humidogram known as the deliquescence point (DRH) and efflorescence point (ERH).
DRH and ERH represent critical junctures in the humidogram where the solid core of the aerosol dissolves or crystallizes, triggering sudden alterations in particle size [13]. It’s important to note that DRH and ERH are distinct points, emphasizing the significance of an aerosol’s history in determining its growth factor at a given RH [14].
The hygroscopic growth of the aerosol may have a significant impact on its undesirable effects on the human respiratory tract. Inhaling aerosols can deeply penetrate the respiratory system, impacting lung function and potentially aggravating respiratory conditions. Notably, in Poland, comprehensive studies regarding the hygroscopic properties of aerosols have not been conducted, creating a gap in understanding their behavior in the atmosphere. Aerosols’ hygroscopic nature significantly influences their behavior, altering their size, composition, and ability to act as cloud condensation nuclei. Moreover, the optical properties of aerosols undergo considerable modifications when exposed to varying RH levels. These changes in optical properties play a vital role in visibility, atmospheric radiative budget, and climate dynamics. Therefore, investigating aerosol hygroscopicity becomes critical for assessing health risks and understanding the broader atmospheric and climatic implications. Addressing this research gap within the context of Poland is crucial for accurately modeling climate changes induced by aerosol-humidity interactions.
The article has the following structure in Section 2 are presented instruments and in Section 3 digital tools used in data acquisition and analysis. Section 4 focus was on utilizing OPC-N3 alongside ACS 1000 (Aerosol Conditioning System) to gauge the hygroscopic nature of urban aerosols. The data were gathered from the κ and GF measurements taken at the Radiation Transfer Laboratory (RTL) site between 2020 and 2022. Given the absence of prior studies on aerosol hygroscopicity in Poland, our aim was to identify typical values. Additionally, in Section 4.3 we sought to uncover instances of heightened hygroscopicity and their underlying causes, which led us to present two case analyses. The Section 5 introduces a method for adjusting OPC-N3 PM values based on the surrounding RH. The latter part (Section 6) summarised the obtained results. Finally, Section 7 comprises the conclusions drawn from the results.

2. Instruments

2.1. OPC-N3

The OPC-N3 is a small (105 g) low-cost optical particle counter developed by Alphasense Ltd. It is composed of a diode laser of 658 nm wavelength and an elliptical mirror that reflects light to a photodiode detector. The flow is perpendicular to the laser beam and a fan is used to maintain a continuous flow. The device can measure the particle number concentration (PNC) in 24 bins, ranging from 0.35 to 40 μ m diameter. The PNC is then converted to PM 1 , PM 2.5 and PM 10 . As shown by Nurowska et al. [15], OPC-N3 is capable of registering water droplets, in addition to registering aerosols.
As stated in [16], previous studies report that Alphasense OPCs demonstrate strong correlations (R2 = 0.93–0.99) with PM10 in controlled laboratory investigations [17,18,19,20]. However, field-based studies have indicated relatively lower correlations (R2: 0.53–0.8) [21,22,23], primarily due to the variability in ambient meteorological conditions and fluctuations in PM compositions. More information on Alphasense OPC-N3 can be found in [24].

2.2. ACS 1000

Acoem (previously known as Ecotech) Aerosol Conditioning System ACS 1000 (Aerosol Conditioning System) is an indoor instrument capable of measuring the hygroscopic properties of aerosol particles. The ariflow is from the top to the bottom of the aparatus; at the top, air enters a tube where it is dehumidified (RH < 40 % ) and then divided into two separate flows. At the end of two tubes, in which the air can be subjected to various conditions, there are outlets for connecting an external apparatus to carry out air measurements. One flow passes through a humidifier tube in which RH can be changed in a range from RH 10 % above ambient relative humidity (minimum RH 40 % ) up to 90 % (we refer to this tube as wet channel). The second flow, without any additional conditions, is directed to the connected measuring instruments. The flow rate can be adjusted in the range from 1 to 10 L·min 1 . The ACS 1000 can change RH constantly or pass through some points in cycles. One cycle is set as a continuous change in relative humidity from the minimum level RH m i n to the maximum level RH m a x , then returns to the minimum relative humidity level RH m i n . To perform faster cycles the return from maximum to minimum RH was abrupt, therefore the data collected only during the changing of RH from low to high were used for later analysis.

2.3. LAS

The Laser Aerosol Spectrometer (LAS) model 3340, manufactured by TSI, is a device that is used to measure the size distribution of atmospheric aerosol particles in the range from 90 to 7500 nm. It is capable of measuring 100 particle size classes. The measurements were conducted in a 5-min time resolution. This instrument is connected to the air intake system in the Radiation Transfer Laboratory (22 m a.g.l.).

2.4. GIOŚ

Polish Chief Inspectorate of Environmental Protection (GIOŚ) has air quality measurement station at Ave. Niepodległości, Warsaw (52.22 ° N 21.00 ° E) located 1.74 km from RTL. The air inlet of the station is approximately 5 m a.g.l. This is an automatic station measuring: benzene, carbon monoxide, nitrogen monoxide, nitrogen dioxide, nitrogen oxides, PM 10 , PM 2.5 with time resolution one hour (data is available to download [25]).

2.5. MERRA-2

The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) [26] is a long-term global reanalysis from NASA Global Modelling and Assimilation Office, which assimilates space-based observations of aerosols. It has a resolution of 0.5° lat × 0.625° lon × 72 pressure levels. In this paper, we used the 3D aer_NV collection. The data were interpolated for the point of location of the measurements. We extracted the data about the aerosol mass mixing ratios from the lowest model layer (72: The vertical layers of MERRA-2 are arranged top-down). Concentrations of PM 10 and PM 2.5 have been calculated using Equations (1) and (2) [27]:
PM 1 [ μ gm 3 ] = ( 0.7 [ DU 1 ] + [ SS 1 ] + [ SS 2 ] + [ OC p h i l i c ] + [ OC p h o b i c ] + [ BC p h i l i c ] + [ BC p h o b i c ] + 1.375 [ SU ] ) ρ a i r 1 × 10 9
PM 10 [ μ gm 3 ] = ( [ DU 1 ] + [ DU 2 ] + [ DU 3 ] + 0.74 [ DU 4 ] + [ SS 1 ] + [ SS 2 ] + [ SS 3 ] + [ SS 4 ] + [ OC p h i l i c ] + [ OC p h o b i c ] + [ BC p h i l i c ] + [ BC p h o b i c ] + 1.375 [ SU ] ) ρ a i r 1 × 10 9
Brackets denote mass concentration of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SU—sulfate aerosol. To calculate the GF from MERRA-2 data we used the hygroscopic growth factors as in [28]. Table 1 show the GF at different RH for OC, BC, SU, SS [28].

2.6. HYSPLIT

HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) is a model developed in 1998 by the National Oceanic and Atmospheric Administration (NOAA) and the Australian Bureau of Meteorology Research Centere. HYSPLIT is capable of computation of air parcel trajectories. In this article it was used to calculate the back trajectories of air parcels that flowed over Warsaw. We calculated the origin of air parcels 24 and 48 h before arriving at the Radiation Transfer Laboratory, Warsaw. The calculations were done for 00:00UTC and 12:00UTC for altitude 500 m using meteo data GDAS.

3. Methods

3.1. Data Acquisition

The experiment was carried out at the Institute of Geophysics of the University of Warsaw. The building is located in an urban area, location 52°12′40.7″ N 20°58′58.8″ E. The RTL is situated on the roof of the Faculty of Physics, where measurements of the optical and microphysical characteristics of atmospheric aerosols and clouds are constantly taken (e.g., with LAS 3340), as well as components of radiation fluxes and sensible and latent heat fluxes at the Earth’s surface. The basic atmospheric parameters are collected by the Vaisala WXT 520 weather station. The platform is located 20 m above the ground.
The inlet for the sampled air was located on the roof of the building (22 m a.g.l) and connected to the ACS 1000. At the inlet, air was drained (RH < 40 % ) and then divided into two routes. In the wet pipe, the RH was changed in cycles (the RH was changed from around 55 to 90% and back to 55%—we refer to such a run as one cycle). In the other tube, the airflow was not modified (dry channel). At the outlets of the wet and dry channels, two identical devices were connected-OPC-N3.
The data analyzed come from the period between 2020 and 2022. The data were not assembled continuously (during autumn measurements at another location were performed). Data were collected in months between December and June; Table 2 presents the respective measurement times.

3.2. Data Calibration

By manual ACS 1000 should be able to operate in the range of RH between RH 40 % ( 10 % above ambient RH) up to 90 % . During winter, it was possible to obtain an average minimum RH below 55 % and for spring 60 % (The obtained minimum RH for each cycle was with standard deviation ± 4%). We were always trying to reach the minimum possible RH in the cycle.
The OPC-N3s were calibrated to each other seven times (listed in Table 2), approximately every new month of measurements. The calibration of both OPC-N3 was performed for RH below the threshold listed in Table 2.
Data from OPC-N3 were reported every 10 s and ACS 1000 every 6 s. In the first step, the PM 1 , PM 2.5 and PM 10 from OPC-N3 were interpolated with the data from ACS 1000. Second, both OPC-N3 sensors were calibrated to each other by the least squares fit, the fitting coefficients are listed in Table A1. To denoise the data, a running average of 12 points was used.

3.3. κ -Köhler Theory

To obtain information about hygroscopicity of atmosphere aerosol the κ parameter proposed by Crilley et al. [21] was calculated. This parameter was obtained for each cycle performed by ACS100. As stated above (Section 2) the data collected during changing of RH from low to high were used for analysis. Figure 1d show one cycle of ACS100 with the RH from dry/wet brunch, Figure 1e,f show the PM 1 , PM 2.5 and PM 10 respectively from both brunches.
The growth factor (GF) is defined by the equation:
GF = PM w e t PM d r y ,
where PM denotes particulate matter and w e t denotes conditions in which RH is high and d r y indicates some reference level in which RH is low (RH < 40 % ).
Crilley et al. [21] proposed a formula to determine the dimensionless hygroscopicity coefficient κ .
GF = 1 + ρ w ρ o p c κ 1 + 100 RH ,
where GF is an aerosol growth factor, RH denotes relative humidity, ρ w —water density 1.0 [g·cm 3 ] and ρ o p c —internal aerosol density parameter set in OPC-N3 sensors at 1.65 [g·cm 3 ]. The developed curve for κ describe the relation of hygroscopic growth of particle and RH above the deliquescent point. Extending the fit to lower RH values overestimates the hygroscopic growth of the aerosol (see Figure 1), the same effect was observed by Crilley et al. [21]. In Figure 1a,b, the deliquescence point (DRH) is visible, when the core of the aerosol particle dissolves in the water that pours over the body of the particle.
The κ was calculated for each cycle. Figure 1a–c shows obtained κ for one cycle from PM 1 , PM 2.5 and PM 10 data respectively. The non-linear fit of Formula (4) was performed on the GF data as a function of RH. The Levenberg-Marquardt non-linear least squares algorithm was used [29]. Sometimes the fitted κ was below 0. For cases that κ was between −0.05 and 0, the value of κ was changed to 0 (This refers to a case when particles were not hygroscopic and did not grow in high RH). If fitted κ was lower than 0.05 such cases were excluded from further analysis.
To calculate the GF, the following method was used for each cycle (graphically illustrated in Figure 2). The data from the wet brunch were used to calculate the GF. For each cycle, the average value of PM was calculated at time ( t 0 ), defined as a moment when the RH in the wet chamber was up to 5 percentiles.
As during the cycle, the external conditions could change, data from dry brunch were used to correct data at the wet brunch. It was determined what is the change in PM d r y between a given moment and t 0 . The difference thus determined was added to the PM w e t at a given time point. The data in the next step was divided by the mean PM w e t value at the time t 0 .
GF = PM w e t ( t ) PM w e t ( t 0 ) + Δ PM d r y ( t )
Δ PM d r y ( t ) = PM d r y ( t ) PM d r y ( t 0 )

4. Results

4.1. κ Half Year Variability

During one day several cycles were performed, therefore to show half year variability the data were averaged for each day. The Figure 3 show the calculation of κ parameter for three consecutive years (2020–2022) for months December–June, each panel corresponds to calculation of κ from a different particle matter size: panel (a) PM 1 , panel (b) PM 2.5 , panel (c) PM 10 .
In the figure it is possible to observe that during the months December-March (later refereed as winter) κ PM 1 ( κ obtained from PM 1 data) is mostly below 0.05. However, more variation is visible for κ PM 2.5 , κ PM 10 ( κ obtained, respectively, from PM 10 , PM 2.5 data). For winter period only mixture of large particles (above 1 μ m) is subject to hygroscopic growth, and mixture of particles which diameter is smaller than 1 μ m is consisting mostly of non hygroscopic aerosols ( κ below 0.1).
The same behavior can be seen if instead of κ parameter is calculated GF (see Figure 4). There is a very high Pearson correlation coefficient between GF and κ calculated from the same PM sizes. The correlation coefficient between κ PM 1 vs. GF PM 1 is 0.98, κ PM 2.5 vs. GF PM 2.5 is 0.97 and κ PM 10 vs. GF PM 10 is 0.89. All possible correlation coefficients between GF and kappa for different particle sizes are presented in Table 3. There is a high correlation coefficient between GF PM 1 vs. GF PM 2.5 (0.78) and κ PM 1 vs. κ PM 2.5 (0.77) and medium correlation coefficient between GF PM 10 vs. GF PM 2.5 (0.63), κ PM 1 vs. κ PM 2.5 (0.62). Correlations of GF and κ obtained from PM 10 shows a weak correlation with respect to GF and κ obtained from PM 1 . Due to the low concentration of particles of several micrometers in size, the PM 10 data show significant Poisson noise. The lower correlation may also be due to the different chemical composition of small and large particles, which is why they are subject to different hygroscopic growth processes.
During winter κ PM 1 varied from 0 to 0.11 (mean value 0.04 ± 0.04 ), κ PM 2.5 0 to 0.24 ( 0.06 ± 0.07 ) and κ PM 10 0 to 0.60 (mean value 0.12 ± 0.17 ). During the spring κ PM 1 varied from 0.05 to 0.23 (mean value 0.13 ± 0.13 ), κ PM 2.5 0.06 to 0.23 (mean value 0.13 ± 0.13 ) and κ PM 10 0.11 to 0.48 (mean value 0.13 ± 0.13 ). During the winter GF PM 1 varied from 1 to 1.48 (mean value 1.16 ± 0.13 ), GF PM 2.5 1.02 to 2.00 ( 1.25 ± 0.27 ) and GF PM 10 1.02 to 4.58 (mean value 1.51 ± 0.76 ). During the spring GF PM 1 varied from 1.24 to 1.97 (mean value 1.52 ± 0.23 ), GF PM 2.5 1.28 to 2.00 (mean value 1.53 ± 0.23 ) and GF PM 10 1.32 to 2.80 (mean value 1.53 ± 0.23 ). The κ /GF PM 1 and κ /GF PM 2.5 for winter were mostly below 0.05/1.3 and for spring mostly above 0.05/1.3.

4.2. κ Diurnal Variation

We were curious if the higroscopicity depends on diurnal variation. The Figure 5 presents the diurnal variation of κ parameter. Periods of winter and spring were depicted separately as the air aerosol composition differ between those seasons. For winter we found out (using one-way analysis of variance ANOVA) a statistically-significant difference in κ PM 1 , κ PM 2.5 and κ PM 10 according to time of day (F PM 1 (statistical F-test)(7) = 5.1, F PM 2.5 (7) = 3.15, F PM 10 (7) = 3.02, p < 0.05). A Tukey’s honestly significant difference criterion (TukeyHSD) revealed significant pairwise differences during winter season Table 4:
  • for κ PM 1 between time slot 18–21 and 3–6, 6–9, 9–12, with a difference of κ above 0.01 (p < 0.05),
  • for κ PM 1 between time slot 3–6 and 15–18, 18–21, 21–00, 00–03, with a difference of κ above 0.01 (p < 0.05),
  • for κ PM 2.5 between time slot 18–21 and 3–6, 9–12, with an average difference of κ 0.03 (p < 0.05),
  • for κ PM 10 between time slot 9–12 and 15–18, 18–21, with an average difference of κ 0.07 (p < 0.05).
For spring we found out a statistically-significant difference in κ PM 1 , κ PM 2.5 according to time of day (F PM 1 (7) = 3.51, F PM 2.5 (7) = 2.91, p < 0.05). A TukeyHSD revealed significant pairwise differences:
  • for κ PM 1 between time slot 12–15 and 18–21, 21–00, with an average difference of κ 0.05 (p < 0.05),
  • for κ PM 2.5 between time slot 18–21 and 9–12, 12–15 with an average difference of κ 0.05 (p < 0.05),
For κ PM 10 during spring period no statistically-significant difference was found.
For winter, there is a significant difference in mean between the 9–12 and 18–21 time slot; and during spring, between the 12–15 and 18–00 time slots. Elevated values of κ during sun operating hours may be related to the production of photocemical particulate matter. As shown for the Los Angeles area by Dale A. Lundgren [30], during dense photochemically produced smog, particulate matter was mainly a water-soluble nitrate compound, which was very hygroscopic. In the Figure A1 we present the diurnal variation of NO x and NO 2 with regard to time of a day. The correlation between mean daily concentration of NO x /NO 2 and κ is not statistically significant, although it is slightly positive for spring and negative for winter. Based on this, it’s difficult to determine a direct link between aerosol hygroscopicity and public transportation emission.

4.3. Case Study

4.3.1. Winter 2021/2022 (15 December 2021–31 January 2022)

To show the picture of winter period regarding hygroscopicity of aerosol we have chosen a period from 16 December 2021 till 31 January 2022. During winter time values of κ obtained in this study are low, close to zero. However, events occur where these values arise. As illustrated in Figure 6 and Figure 7 during episodes the values of κ PM 1 are slightly increasing, κ PM 2.5 more, and the biggest change is in κ PM 10 . The additional figures presenting meteorological variables are shown in Appendix C Figure A2 and Figure A3.
First event of higher κ occurred between 16 December and 21 December. During 16–20 December the air pollution does not exceed norms except PM 10 on December 17 (24-h average norm for PM 10 concentration recommended by the World Health Organization (WHO) is 50 μ g·m 3 and for PM 2.5 25 μ g·m 3 [5]). Between 16–19 December 2021, the weather over Poland was shaped by a large high located over the British Isles (pressure in the high center 1040 hPa). Together with the wide low initially located over northern Finland and moving over western Russia, they caused an increase in the pressure gradient, and thus an increase in wind speed and gusts bringing moist polar sea air over Poland. The wind was blowing from the north-west direction. On the night of 19–20 December, a cold front passed through Poland from north to south, which caused the exchange of air masses from polar moderately warm, stable air to polar, colder, thermodynamically unstable, contributing to convective phenomena. The cold front cleaned the air values of PM 10 and PM 2.5 were below 15 μ g·m 3 till midday of 20 December, when PM 10 values increased to almost 100 μ g·m 3 and next day till 200 μ g·m 3 . From 21 December, the expanding high pressure from Great Britain began to move towards the east, over the North Sea and western Germany, and later over Poland, Romania. Polar Sea air was replaced by the influx of arctic air–much colder air masses [31,32] (The synoptic overview was made on the basis of synoptic maps created by Polish Institute of Meteorology and Water Management–National Research Institute [32] and commentary of the forecaster of the numerical weather forecast model created by the Interdisciplinary Center for Mathematical and Computational Modeling of the University of Warsaw [31]).
During 16–20 December, the registered values of κ for PM 10 ranged mostly from 0.1 to 0.6. Using HYSPLIT model was calculated where the air coming to RTL was 48 h prior. For 16–20 December it was over the region of Northern Sea and Norwegian Sea. HYSPLIT back trajectories show that mostly 50 % time the air trajectories passed over the sea (Figure 6d) (To calculate the path over the sea the coastline map of 10 m resolution was used [33]). For the winter period the Pearson correlation coefficient between κ or GF was calculated and the percentage of the distance traveled over the sea 24/48 h earlier (for κ PM 10 it is 0.32/0.39 and for GF 0.34/0/39 respectively). The correlation is better for PM 10 than for PM 1 or PM 2.5 probably due to the different composition of the particulate. The Pearson correlation coefficient for κ PM 10 and GF PM 10 and the percentage of sea salt in the mixture of particles up to 10 μ m id 0.41 and 0.39 respectively. All Pearson correlation coefficients for κ /GF and effective radius, GF from MERRA-2 meteorogical variables are listed in Table A2.
During elevated values of κ PM 10 the κ PM 1 does not show a similar behavior. We retrieved from MERRA-2 the reanalysis of the composition of PM and calculated the GF for RH = 85 % (Figure 6d). The chemical composition of MERRA-2 for PM 1 and PM 10 is also different. During elevated values of κ PM 10 in the reanalysis is visible a greater contribution to PM 10 mass of sea salt (Figure 6g) and the contribution of sea salt to PM 1 mass is small (Figure 6f). The calculated GFs from MERRA-2 show elevated values of GF PM 10 and no elevated values of GF PM 1 , which is in line with our observations.
We were expecting to see some variability of r e f f in time. However the range of variability of r e f f is relatively small.
From 21 December around 10 a.m., the values of κ PM 10 dropped to values close to zero (0.02). After days with temperatures above zero, the arctic air brought negative temperatures. The arrival of a high-pressure over Poland resulted in weak baric gradient, slower flow of air masses and slower wind gusts. With the strengthening of the high pressure over Poland on 22 December, the influx of air from the north weakened. From 22 December the PM 10 and PM 2.5 values till 25th exceeded the standards set by WHO, PM 10 reached maximum of 140 μ g·m 3 on 11 a.m. 23 December. The advance of the low pressure over Finland from the north of Norway caused a reversal of the wind direction from NW to west on 22 to 23 December, and then to the south-west (24 December). On 24 December, warm air flowed over Poland, causing snow, sleet, rain and freezing rain. During those days the κ value for PM 1 , PM 2.5 and PM 10 was mostly around 0.02.
On 25 December, over northern Finland, a low tide developed with a gulf stretching over Belarus. It brought frosty Arctic sea air over northern Poland. This caused another short event of κ PM 10 above 0.2, while the air was clear, the values of PM 1 , PM 2.5 dropped below 25 μ g·m 3 .
From 26 December, there was a high over Poland, providing sunny and windless weather on 27 December. On 28 December, the high was displaced by the Atlantic low approaching from Great Britain. Together with the Mediterranean low, they caused an increase in wind speed on 28 December, incoming air from south-eastern directions [31,32]. The weather from 28 December (till 6 January) was mainly controlled by the appearance of successive Atlantic lows that favored the influx of warm air, mainly from the south-west and west directions. These lows moved along the northern side of the quasi-stationary high of the Azores. From 30 December, the temperatures were above zero degrees. When the coming air to Poland was from south-west directions the κ values where increasing from 0 at 27 December to 0.13 on 31 December. During days 26 and 31 December the air pollution was high, the norms of PM 1 , PM 2.5 where exceeded.
The last performed cycle in 2021 was during the New Year celebration. Before New Years firework shoot the κ PM 1 , κ PM 2.5 , κ PM 10 where respectively 0.05, 0.08, 0.09. For the last cycle of κ PM 1 , κ PM 2.5 and κ PM 10 were slightly elevated (0.09, 0.12 and 0.14 respectively). After midnight, the values returned to values similar to before midnight (0.06, 0.07 and 0.09 κ PM 1 , κ PM 2.5 , κ PM 10 respectively). The Figure A4 show the humidograms for growth factor of PM 1 , PM 2.5 and PM 10 as well as registered values of PM and RH in ACS 1000. The atmosphere before New Year was clear (PM 1 around 1.5 and PM 2.5 2 [ μ g·m 3 ]), we observed a jump of PM 1 to 6.5 and PM 2.5 to 8 [ μ g·m 3 ] 6 min after midnight. At that time, the RH in ACS 1000 was 81.5%.
The particulate matter emitted did not show an abnormal change in GF for this cycle. This is in contrast to previously reported works in which fireworks were identified to consist of ion salts, mostly sulfate, potassium and chloride, sodium, and magnesium, which show high hygroscopicity [34,35].
In the subsequent month, January 2022, similar behavior was observed. Whenever the MERRA-2 reanalysis shows more than 50% of sea salt contribution into PM 10 mass the κ PM 10 increases from around 0.05 to values greater than 1.0 and the GF 85 % for PM 10 from MERRA-2 also increases to values greater than 2.2. When periods such as 8–13 January appear, when the air masses did not pass over sea (the sea salt contribution is low), then the values of κ are below 0.05. During this period, the concentration of PM 2.5 were above 25 μ g·m 3 and PM 10 above 50 μ g·m 3 .
Figure 8 presents the air flow path with the assigned value κ PM 10 (denoted by colored circle). The grater the length of the path is the quicker the air parcel was moving towards RTL at Warsaw. Values of κ PM 10 above 0.2 are mostly associated with the fast flow of air from the North Atlantic, Northern Sea, Norwegian Sea, or Barents Sea.

4.3.2. Spring 2020 (28 April–5 May, 16 May–24 May)

To show the picture of the spring measurements regarding the aerosol hygroscopicity a period was chosen from 28 April to 24 May 2020 (with a break between 5 and 16 May). Figure 9 presents the values κ for spring 2020 and Figure A5 the meteorological conditions at the RTL station during that time. Obtained values of κ were oscillating near value 0.2. Most of the time κ PM 1 , κ PM 2.5 where having similar values. The events when κ was greater than 0.4 where occurring only for κ PM 10 this was for days 20–24 May.
The 24-h mean of the concentrations of PM 2.5 and PM 10 from the GIOś station did not exceed WHO limits. In the winter case (apart from smog events), the concentrations of PM 2.5 and PM 10 were similar; in the spring case the lines in Figure 9a are separated, suggesting that there are more large particles (above 2.5 μ m) in the air aerosol compared to the winter period. The composition of the PM mass is different than in winter. For PM 1 contribution of sulphates is more than 40 % , with visible diurnal cycle, black carbon provide less to the total mass than in winter season. However, dust gives a higher contribution to total mas of PM than in winter. Nevertheless MERRA-2 reanalysis show events of sea salt particles in the air there is no visible correlation of it with κ . Using the HYSPLIT model, the percentage of incoming air over the sea was calculated; the correlation between kappa and distance over the sea is also low. Figure 10 presents what was the path of incoming air to RTL in previous 48 h. There is no specific region for which would be seen that the κ is grater.

5. Calibrating OPC

As we have shown, the PM values are biased by high RH. To calibrate the OPC we transformed Equation (4) to obtain PM as it would be registered if RH was zero (refereed as PM RH = 0 ). The OPC-N3 built-in hygrometer shows lower RH values than those registered by ASC1000. Our guess is that the OPC-N3 hygrometer is near to the control board, which causes heating, leading to underestimation of RH. Therefore first we calibrated the RH from OPC-N3 versus ACS 1000 and then calculated the PM RH = 0 .
Figure 11b presents relative humidity registered by OPC-N3 versus RH form ACS 1000. For Period nr 1 (Spring 2020), we fit a line to the RH data from OPC-N3 versus ACS 1000 (orange line Figure 11b). We obtained the following fitting:
RH A C S 100 ± Δ = a · RH O P C + b ± Δ = 1.17 · RH O P C + 18.33 ± 1.40 ,
By transforming Equation (4), we get the formula for PM d r y :
PM RH = 0 = ( 1 a · RH O P C + b ) · PM w e t 1 + ( ρ w ρ o p c κ 1 ) · ( a · RH O P C + b ) ,
where PM w e t is the PM measured by OPC-N3, RH O P C is the relative humidity from OPC-N3, a and b are the estimated parameters from Equation (7). For calibration we used the mean κ value for the period 29–30 April 2020 ( κ = 0.16 ).
Figure 11 shows the comparison of PM w e t and PM d r y over a period of two days. The upper panel shows the data from OPC-N3 PM w e t (blue color) compared to the PM data from OPC-N3 PM d r y (orange color) and the calibrated PM RH = 0 (green color). Lower panel presents the cycles of RH in ACS100: in pink the RH from dry brunch and in black from the wet brunch.
Figure 11 displays the PM 1 data obtained from both the wet brunch (blue line in panel a) and the dry channel (orange line in panel a). In Figure 11c, the graph illustrates the temporal changes in relative humidity within the wet brunch (depicted by the black line), while the pink line represents the RH in the dry channel.
Analysis of Figure 11 indicates that the PM values measured at high RH, following the proposed correction (green line in Figure 11a), closely resemble the PM levels measured at RH below 30% (orange line in Figure 11a). In particular, on 29 April, PM 1 appeared more hygroscopic compared to 30 April, suggesting that applying the same correction factor ( κ ) for both days works better for the latter day.
Despite the presence of residual variations with respect to RH after the correction, these fluctuations are considerably smaller.
In our study, the hygroscopicity of the atmospheric aerosol mixture obtained from PM 1 varied from almost 0 to κ = 0.2. We wanted to know how much differ the values of PM 1 , R H obtained at high ambient RH (where RH denotes = 60 / 65 / 70 / 75 / 80 / 85 / 90 % ) with the PM 1 , R H < 40 % obtained in RH below 40 % . We chose three different days 27 December 2020, 29 March 2020 and 4 May 2020 where κ was respectively 0.01, 0.13 and 0.20. Table A3 shows for those days the PM 1 registered at RH < 40 % and at RH = 60 / 65 / 70 / 75 / 80 / 85 / 90 % , apart from values in [ μ g· m 3 ] the PM 1 , R H is presented as a percentage of value PM 1 , R H < 40 % . For a day when κ was 0.01 even for the RH of 85 % the error made is within 10%, however, if κ was 0.13 the error is greater. For κ = 0.13 the error is bigger than 20% for RH exceeding 75 % , in case of RH = 90% the error is 61.33%. For κ = 0.20 the the error bigger than 20% is already for RH exceeding 70 % . When the RH = 90% the obtained values of PM 1 are more than twice bigger than registered in dry conditions.
In addition to checking the error generated by recording PM values without prior drying the air, we wanted to check how correcting the values of PM registered in high RH using κ helps to improve the results. For each day 27 December 2020, 29 March 2020 and 4 May 2020 we corrected the PM 1 , R H using κ = 0.05 / 0.10 / 0.15 / 0 / 20 / 0.25 , the results are shown in the same Table A3. For each day when using κ closest to the estimated value, the percentage error between the corrected value and the measured PM 1 , R H < 40 % value was below 10 % . In the case of estimated κ = 0.01 using for correction larger values of κ caused an underestimation of PM 1 by up to 20%. The biggest errors occur when selecting the wrong kappa for correction when the actual kappa value is large.

6. Summary of the Results

In this work we show that a set of two OPC-N3 and an ACS 1000 can be used to study the hygroscopicity of an aerosol. During three years, we measured κ (a dimensionless parameter specifying the hygroscopicity of the aerosol [21]) during the months December–June. This were the first measurements of outdoor aerosol hygroscopicity conducted in Poland. Overall registered values of κ in the spring season (April–July) were higher— 0.13 ± 0.13 , 0.13 ± 0.13 , 0.13 ± 0.13 ( κ PM 1 , κ PM 2.5 , κ PM 10 respectively) than in winter (January–March)— 0.04 ± 0.04 , 0.06 ± 0.07 , 0.12 ± 0.17 ( κ PM 1 , κ PM 2.5 , κ PM 10 respectively). There is a high correlation between κ and GF calculated from the same PM size: 0.98, 0.97 and 0.89 (PM 1 , PM 2.5 , PM 10 respectively). Using the GF parameter has better variability, so the error values are not the same magnitude as the given value. For spring, the GF was for PM 1 , PM 2.5 , PM 10   1.52 ± 0.23 , 1.53 ± 0.23 , 1.53 ± 0.23 respectively; and for winter it was 1.16 ± 0.13 , 1.25 ± 0.27 , 1.63 ± 0.26 , respectively. We have shown that
  • the hygroscopicity parameter ( κ ) for PM 1 during winter (December–March) is below 0.05. These periods were when pollution was high (PM 2.5 and PM 10 WHO norms were exceeded). This is likely that the air pollution consists mostly of soot, from home heating systems, which is non-hygroscopic.
  • During winter there are events when the hygroscopicity of PM 1 is low and the hygroscopicity of PM 10 is high up to 1. These events are probably associated with fast flow of air masses from the North Atlantic Ocean during low pressure circumstances. The PM 10 mass in more than 50 % consists of sea salt particles.
  • There is a dependence of κ on the time of day. During the sun’s operating hours, the values of κ are statistically higher than in the evening/night.
  • During the New Year’s Eve midnight (when the fireworks are shoot) the κ parameter for PM 1 , PM 2.5 and PM 10 raised from values 0.05, 0.08, 0.09 to 0.09, 0.12 and 0.14 respectively. We did not observe higher values of κ during fireworks shooting.
Apart from hygroscopicity measurements we have tested how to correct low-cost particulate matter counters such as OPC-N3 for RH. Our recommendations for performing air quality measurements with the OPC-N3 device without using a dehumidifier:
  • The data can be corrected using κ , however appointing the appropriate kappa for corrections is difficult.
  • If one suspect that the aerosol may be hygroscopic, it is worth taking into account the kappa correction, because PM values registered at high RH with respect to those measured in RH < 40 % may differ more than twice.
  • The RH for correction should be taken from a device different from OPC-N3. The RH values from OPC-N3 are significantly lower than ambient and cannot be easily corrected, as the black OPC-N3 body absorbs sun light, changing the amount of bias.
  • For cases when κ is small (for example for winter) registering PM at ambient RH causes a small error up to 10%.
  • We do not recommend using PM values from OPC-N3 for ambient RH above 80 % without any correction.

7. Conclusions

This study underscores the feasibility of employing OPC-N3 devices alongside an ACS 1000 to evaluate aerosol hygroscopicity, especially in different seasons. Our findings highlight distinct seasonal variations in κ values, with elevated hygroscopicity observed during spring compared to winter, indicative of varying aerosol compositions and sources.
Crucially, our observations regarding the values of κ during high pollution periods in winter suggest a prevalence of hydrophobic particles (such as soot), possibly from domestic heating systems. Additionally, meteorological events, such as rapid air mass movements from the North Atlantic Ocean, significantly influenced aerosol hygroscopicity.
The temporal dependency of κ values throughout the day and the impact of specific events, further emphasize the dynamic nature of aerosol hygroscopicity.
Furthermore, our study provided practical recommendations for correcting low-cost particulate matter counters like OPC-N3 for RH, emphasizing the significance of considering κ values in correcting aerosol measurements, especially in conditions of varying RH.
In conclusion, this research contributes valuable insights into the understanding of aerosol hygroscopicity patterns in the Polish urban area, emphasizing the importance of seasonal variations, meteorological influences, and the potential for improved correction methods for RH-dependent aerosol measurements.

Author Contributions

Conceptualization K.M.M. and K.N.; methodology, K.N.; formal analysis, K.N.; investigation, K.N.; writing—original draft preparation, K.N.; writing—review and editing, K.N.; supervision, K.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research is part of the OPUS project Impact of aerosol on the microphysical, optical and radiation properties of fog, which was funded by the National Science Center grant number UMO-2017/27/B/ST10/00549.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ongoing transformation within our Faculty to consolidate information into a single database.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACS 1000Acoem Aerosol Conditioning System
AODaerosol optical depth
BCBlack carbon
DRHDeliquescent point
DUDust
ERHEfflorescence point
FF-test (statistics)
GFGrowth factor
HYSPLITHybrid Single-Particle Lagrangian Integrated Trajectory model
κ Dimensionless parameter proposed by Crilley et al. [21]
IPCCIntergovernmental Panel on Climate Change
LASLaser Aerosol Spectrometer
MERRA-2Modern-Era Retrospective analysis for Research and Applications, Version 2
OCOrganic carbon
OPC-N3Low-cost optical particle counter developed by Alphasense Ltd.
pp-value (statistics) or pressure [hPa]
PM 1 , PM 2.5 , PM 10 Particulate matter of diameter less than or equal to
1 μ m, 2.5 μ m and 10 μ m, respectively
PNCparticle number concentration
r e f f Effective radius
RHRelative humidity [%]
RH m i n Minimum relative humidity possible to achive with ACS 1000 system
RH m a x Maximum relative humidity possible to achive with ACS 1000 system
ρ w water density 1.0 [g ·cm 3 ]
ρ o p c internal aerosol density parameter set in OPC-N3 sensors at 1.65 [g ·cm 3 ]
RTLRadiative Transfer Labolatory, Geophysic Insitute, Faculty of Physics,
University of Warsaw
SEStandard error
SUSulfate aerosol
SSSea salt
TTemperature [°C]
UAVUnmanned aerial vehicle
VWind [km·h 1 ]

Appendix A. GIOŚ Station Ave. Niepodległości, Warsaw, NOx

Figure A1. Boxplot of diurnal variation of NO x (panel a,c) and of NO 2 (panel b,d) obtained for GIOŚ station Ave. Niepodległości, Warsaw based on data from 2020–2022. Upper panels presents data collected for winter (December–March) lower panels presents data for spring (April–June).
Figure A1. Boxplot of diurnal variation of NO x (panel a,c) and of NO 2 (panel b,d) obtained for GIOŚ station Ave. Niepodległości, Warsaw based on data from 2020–2022. Upper panels presents data collected for winter (December–March) lower panels presents data for spring (April–June).
Atmosphere 15 00061 g0a1

Appendix B. OPC-N3 Calibration

Table A1. Linear fit factors y = ax + b used to calibrate OPC-N3 data for PM 1 , PM 2.5 , PM 10 . Only values of SE needs to be multiplied by 10 2 .
Table A1. Linear fit factors y = ax + b used to calibrate OPC-N3 data for PM 1 , PM 2.5 , PM 10 . Only values of SE needs to be multiplied by 10 2 .
PeriodCalibrationa ± SE ∗ 10−2b ± SE ∗ 10−2
nrfor
1PM 1 0.64 ± 0.11−0.12 ± 1.03
PM 2.5 0.67 ± 0.170.03 ± 2.00
PM 10 0.57 ± 0.440.89 ± 6.88
2PM 1 0.85 ± 0.100.20 ± 0.50
PM 2.5 0.82 ± 0.130.29 ± 0.81
PM 10 0.79 ± 0.270.47 ± 1.80
3PM 1 0.88 ± 0.050.4 ± 0.57
PM 2.5 0.86 ± 0.060.6 ± 0.79
PM 10 0.84 ± 0.100.82 ± 1.39
4PM 1 0.40 ± 0.120.03 ± 0.18
PM 2.5 0.43 ± 0.310.18 ± 0.78
PM 10 0.24 ± 0.711.52 ± 3.31
5PM 1 1.22 ± 0.080.09 ± 0.36
PM 2.5 1.09 ± 0.090.17 ± 0.51
PM 10 1.00 ± 0.150.45 ± 0.89
6PM 1 1.24 ± 0.05−0.06 ± 0.12
PM 2.5 1.11 ± 0.09−0.09 ± 0.27
PM 10 0.81 ± 0.190.61 ± 0.74
7PM 1 1.23 ± 0.120.01 ± 0.39
PM 2.5 1.13 ± 0.16−0.01 ± 0.63
PM 10 0.75 ± 0.421.18 ± 2.12

Appendix C. Winter 2021–2022

Figure A2. December 2021: Panels represents changes of variables during winter 2021/2022 case. (a) Kappa obtained from data from ACS 1000 and OPC-N3 in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 . (b) wind velocity (blue line) and wind gusts (orange line), (c) wind direction, (d) pressure, (e) temperature, (f) relative humidity. Panels (b,df) presents one hour average from Vaisala WXT 520 instrument, Panels (c) presents 12 h average from Vaisala WXT 520 instrument.
Figure A2. December 2021: Panels represents changes of variables during winter 2021/2022 case. (a) Kappa obtained from data from ACS 1000 and OPC-N3 in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 . (b) wind velocity (blue line) and wind gusts (orange line), (c) wind direction, (d) pressure, (e) temperature, (f) relative humidity. Panels (b,df) presents one hour average from Vaisala WXT 520 instrument, Panels (c) presents 12 h average from Vaisala WXT 520 instrument.
Atmosphere 15 00061 g0a2
Figure A3. January 2022: Panels represents changes of variables during winter 2021/2022 case. (a) Kappa obtained from data from ACS 1000 and OPC-N3 in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 . (b) wind velocity (blue line) and wind gusts (orange line), (c) wind direction, (d) pressure, (e) temperature, (f) relative humidity. Panels (b,df) presents one hour average from Vaisala WXT 520 instrument, Panels (c) presents 12 h average from Vaisala WXT 520 instrument.
Figure A3. January 2022: Panels represents changes of variables during winter 2021/2022 case. (a) Kappa obtained from data from ACS 1000 and OPC-N3 in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 . (b) wind velocity (blue line) and wind gusts (orange line), (c) wind direction, (d) pressure, (e) temperature, (f) relative humidity. Panels (b,df) presents one hour average from Vaisala WXT 520 instrument, Panels (c) presents 12 h average from Vaisala WXT 520 instrument.
Atmosphere 15 00061 g0a3
Figure A4. Figure shows last cycle of the ACS 1000 and registered PM by OPC-N3 in 2021. The change of RH in the ACS 1000 is shown on panel (a), the registered PM 1 , PM 2.5 are shown in panel (b) and in panel (c) PM 10 (after calibration and applying running average of 12 points). Plots (df) presents GF versus RH for one cycle obtained from PM 1 , PM 2.5 , PM 10 respectively, the black line shows fitted curve defined by Equation (4). Time is in UTC, 23:00 UTC is 00:00 at Warsaw.
Figure A4. Figure shows last cycle of the ACS 1000 and registered PM by OPC-N3 in 2021. The change of RH in the ACS 1000 is shown on panel (a), the registered PM 1 , PM 2.5 are shown in panel (b) and in panel (c) PM 10 (after calibration and applying running average of 12 points). Plots (df) presents GF versus RH for one cycle obtained from PM 1 , PM 2.5 , PM 10 respectively, the black line shows fitted curve defined by Equation (4). Time is in UTC, 23:00 UTC is 00:00 at Warsaw.
Atmosphere 15 00061 g0a4
Table A2. Pearson correlation coefficient between variables listed in first column and κ , GF. Colors divide values by strength of correlation.
Table A2. Pearson correlation coefficient between variables listed in first column and κ , GF. Colors divide values by strength of correlation.
r e f f 0.190.310.250.250.400.28
GF PM 1 MERRA-20.140.230.340.160.240.28
GF PM 10 MERRA-2−0.010.270.420.050.300.40
dist sea 24 h0.180.240.320.180.240.34
dist sea 48 h0.150.280.390.160.310.39
% SS−0.010.260.410.040.290.39
T0.560.350.490.480.390.34
p0.060.160.250.040.210.33
RH0.27−0.06−0.010.18−0.03−0.14
V0.350.160.090.340.15−0.03
V g u s t 0.530.330.360.510.350.17
κ PM 1 κ PM 2.5 κ PM 10 GF PM 1 GF PM 2.5 GF PM 10
statistically insignificant
0.0 ≤ |r| ≤ 0.2—no correlation
0.2 < |r| ≤ 0.4—weak correlation
0.4 < |r| ≤ 0.7—medium correclation
0.7 < |r| ≤ 0.9—high correlation
0.9 < |r| ≤ 1.0—very high correlation

Appendix D. Spring 2020

Figure A5. April–May 2020: Panels represents changes of variables during spring 2020 case. (a) Kappa obtained from data from ACS 1000 and OPC-N3 in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 . (b) wind velocity (blue line) and wind gusts (orange line), (c) wind direction, (d) pressure, (e) temperature, (f) relative humidity. Panels (b,df) presents one hour average from Vaisala WXT 520 instrument, Panels (c) presents 12 h average from Vaisala WXT 520 instrument.
Figure A5. April–May 2020: Panels represents changes of variables during spring 2020 case. (a) Kappa obtained from data from ACS 1000 and OPC-N3 in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 . (b) wind velocity (blue line) and wind gusts (orange line), (c) wind direction, (d) pressure, (e) temperature, (f) relative humidity. Panels (b,df) presents one hour average from Vaisala WXT 520 instrument, Panels (c) presents 12 h average from Vaisala WXT 520 instrument.
Atmosphere 15 00061 g0a5

Appendix E. OPC-N3 Correction

Table A3. In the Table are presented values of PM 1 registered during specific date by ACS 1000 and OPC-N3. PM 1 , R H < 40 % refers to data from dry brunch of ACS 1000 where RH was not exceeding 40 % and PM 1 , R H refers to PM 1 measured at specific RH listed in each row. Values are given in [ μ g·m 3 ] or as a percentage of PM 1 , R H < 40 % . Columns from 5 to 9 were obtained by applying Equation (8) to data PM 1 , R H with assumption of one of the κ = 0.05 / 0.10 / 10.15 / 0.20 / 0.25 . At the top of each subtable is presented for which day was created and what was the value of κ obtained for that day.
Table A3. In the Table are presented values of PM 1 registered during specific date by ACS 1000 and OPC-N3. PM 1 , R H < 40 % refers to data from dry brunch of ACS 1000 where RH was not exceeding 40 % and PM 1 , R H refers to PM 1 measured at specific RH listed in each row. Values are given in [ μ g·m 3 ] or as a percentage of PM 1 , R H < 40 % . Columns from 5 to 9 were obtained by applying Equation (8) to data PM 1 , R H with assumption of one of the κ = 0.05 / 0.10 / 10.15 / 0.20 / 0.25 . At the top of each subtable is presented for which day was created and what was the value of κ obtained for that day.
Particulate Matter κ = 0.01Day: 27 December 2020
RH PM 1 , R H PM 1 , R H < 40 % ESTIMATED Value Using   κ
Measured ValueMeasured Value0.050.100.150.200.25
60 %[ μ g·m 3 ]4.204.114.073.963.853.753.65
[%]2 −1−4−6−9−11
65 %[ μ g·m 3 ]4.294.114.154.013.893.773.66
[%]4 1−2−6−8−11
70 %[ μ g·m 3 ]4.414.174.234.073.913.773.64
[%]6 2−2−6−9−13
75 %[ μ g·m 3 ]4.444.134.224.023.843.673.52
[%]8 2−3−7−11−15
80 %[ μ g·m 3 ]4.514.244.274.063.873.693.53
[%]6 1−4−9−13−17
85 %[ μ g·m 3 ]4.534.244.284.053.843.653.48
[%]7 1−5−9−14−18
90 %[ μ g·m 3 ]-------
[%]-------
Particulate Matter κ = 0.13Day: 29 March 2020
RH PM 1 , R H PM 1 , R H < 40 % ESTIMATED value using κ
Measured valueMeasured value0.050.100.150.200.25
60 %[ μ g·m 3 ]15.6014.8814.9114.2713.6913.1512.65
[%]5 0−4−8−12−15
65 %[ μ g·m 3 ]16.0014.9615.1214.3413.6312.9912.41
[%]7 1−4−9−13−17
70 %[ μ g·m 3 ]13.9612.5713.0212.2011.4810.8310.26
[%]11 4−3−9−14−18
75 %[ μ g·m 3 ]15.2212.7613.9012.7911.8411.0310.32
[%]19 90−7−14−19
80 %[ μ g·m 3 ]16.4813.1814.7213.3112.1411.1710.34
[%]25 121−8−15−22
85 %[ μ g·m 3 ]19.0113.0616.0313.8612.2110.919.86
[%]46 236−7−16−25
90 %[ μ g·m 3 ]22.6214.0218.3115.3713.2511.6410.38
[%]61 3110−6−17−26
Particulate Matter κ = 0.20Day: 4 May 2020
RH PM 1 , R H PM 1 , R H < 40 % ESTIMATED value using κ
Measured valueMeasured value0.050.100.150.200.25
60 %[ μ g·m 3 ]2.051.901.961.871.791.721.66
[%]8 3−2−6−9−13
65 %[ μ g·m 3 ]2.872.532.722.592.462.352.25
[%]14 82−3−7−11
70 %[ μ g·m 3 ]2.872.322.682.522.372.242.13
[%]24 1692−3−8
75 %[ μ g·m 3 ]3.812.833.483.202.962.762.58
[%]35 23135−2−9
80 %[ μ g·m 3 ]4.152.763.693.313.012.762.54
[%]51 342090−8
85 %[ μ g·m 3 ]4.922.604.133.563.132.792.52
[%]89 5937207−3
90 %[ μ g·m 3 ]5.732.624.593.833.282.882.56
[%]119 75462510−2

References

  1. IPCC. Summary for Policymakers. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T., Qin, D., Plattner, G.K., Tignor, M., Allen, S., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; Book Section SPM; pp. 1–30. [Google Scholar] [CrossRef]
  2. Laden, F.; Neas, L.; Dockery, D.; Schwartz, J. Association of Fine Particulate Matter from Different Sources with Daily Mortality in Six U.S. Cities. Environ. Health Perspect. 2000, 108, 941–947. [Google Scholar] [CrossRef] [PubMed]
  3. Maciejewska, K. Short-term impact of PM2.5, PM10, and PMc on mortality and morbidity in the agglomeration of Warsaw, Poland. Air Qual. Atmos. Health 2020, 13, 659–672. [Google Scholar] [CrossRef]
  4. European Union. Directive 2008/50/EC of the European Parliament and of the Council of 21 May 2008 on ambient air quality and cleaner air for Europe. Off. J. Eur. Union 2008, L152, 1–44. [Google Scholar]
  5. World Health Organization. WHO Air Quality Guidelines for Particulate Matter, Ozone, Nitrogen Dioxide and Sulfur Dioxide. 2006. Available online: https://apps.who.int/iris/handle/10665/69477 (accessed on 13 November 2023).
  6. European Environment Agency. Air Quality in Europe — 2020 Report No 9/2020; Technical Report; European Environment Agency: Copenhagen, Denmark, 2020. [Google Scholar]
  7. WMO/GAW. Aerosol Measurement Procedures Guidelines and Recommendations; World Meteorological Organization: Geneva, Switzerland, 2003. [Google Scholar]
  8. Zieger, P.; Fierz-Schmidhauser, R.; Gysel, M.; Ström, J.; Henne, S.; Yttri, K.E.; Baltensperger, U.; Weingartner, E. Effects of relative humidity on aerosol light scattering in the Arctic. Atmos. Chem. Phys. 2010, 10, 3875–3890. [Google Scholar] [CrossRef]
  9. Zieger, P.; Väisänen, O.; Corbin, J.; Partridge, D.; Bastelberger, S.; Mousavi-Fard, M.; Rosati, B.; Gysel, M.; Krieger, U.; Leck, C.; et al. Revising the hygroscopicity of inorganic sea salt particles. Nat. Commun. 2017, 8, 15883. [Google Scholar] [CrossRef] [PubMed]
  10. Petters, M.D.; Kreidenweis, S.M. A single parameter representation of hygroscopic growth and cloud condensation nucleus activity. Atmos. Chem. Phys. 2007, 7, 1961–1971. [Google Scholar] [CrossRef]
  11. Moore, R.H.; Ziemba, L.D.; Dutcher, D.; Beyersdorf, A.J.; Chan, K.; Crumeyrolle, S.; Raymond, T.M.; Thornhill, K.L.; Winstead, E.L.; Anderson, B.E. Mapping the Operation of the Miniature Combustion Aerosol Standard (Mini-CAST) Soot Generator. Aerosol Sci. Technol. 2014, 48, 467–479. [Google Scholar] [CrossRef]
  12. Titos, G.; Cazorla, A.; Zieger, P.; Andrews, E.; Lyamani, H.; Granados-Muñoz, M.i.; Olmo, F.; Arboledas, L. Effect of hygroscopic growth on the aerosol light-scattering coefficient: A review of measurements, techniques and error sources. Atmos. Environ. 2016, 141, 494–507. [Google Scholar] [CrossRef]
  13. Vlasenko, S.S.; Su, H.; Pöschl, U.; Andreae, M.O.; Mikhailov, E.F. Tandem configuration of differential mobility and centrifugal particle mass analysers for investigating aerosol hygroscopic properties. Atmos. Meas. Tech. 2017, 10, 1269–1280. [Google Scholar] [CrossRef]
  14. Wise, M.E.; Martin, S.T.; Russell, L.M.; Buseck, P.R. Water Uptake by NaCl Particles Prior to Deliquescence and the Phase Rule. Aerosol Sci. Technol. 2008, 42, 281–294. [Google Scholar] [CrossRef]
  15. Nurowska, K.; Mohammadi, M.; Malinowski, S.; Markowicz, K. Applicability of the low-cost OPC-N3 optical particle counter for microphysical measurements of fog. Atmos. Meas. Tech. 2023, 16, 2415–2430. [Google Scholar] [CrossRef]
  16. Kaur, K.; Kelly, K.E. Performance evaluation of the Alphasense OPC-N3 and Plantower PMS5003 sensor in measuring dust events in the Salt Lake Valley, Utah. Atmos. Meas. Tech. 2023, 16, 2455–2470. [Google Scholar] [CrossRef]
  17. Sousan, S.; Koehler, K.; Hallett, L.; Peters, T.M. Evaluation of the Alphasense optical particle counter (OPC-N2) and the Grimm portable aerosol spectrometer (PAS-1.108). Aerosol Sci. Technol. 2016, 50, 1352–1365. [Google Scholar] [CrossRef] [PubMed]
  18. Samad, A.; Melchor Mimiaga, F.E.; Laquai, B.; Vogt, U. Investigating a Low-Cost Dryer Designed for Low-Cost PM Sensors Measuring Ambient Air Quality. Sensors 2021, 21, 804. [Google Scholar] [CrossRef]
  19. Sousan, S.; Regmi, S.; Park, Y.M. Laboratory Evaluation of Low-Cost Optical Particle Counters for Environmental and Occupational Exposures. Sensors 2021, 21, 4146. [Google Scholar] [CrossRef]
  20. Dubey, R.; Patra, A.K.; Joshi, J.; Blankenberg, D.; Kolluru, S.S.R.; Madhu, B.; Raval, S. Evaluation of low-cost particulate matter sensors OPC N2 and PM Nova for aerosol monitoring. Atmos. Pollut. Res. 2022, 13, 101335. [Google Scholar] [CrossRef]
  21. Crilley, L.R.; Shaw, M.; Pound, R.; Kramer, L.J.; Price, R.; Young, S.; Lewis, A.C.; Pope, F.D. Evaluation of a low-cost optical particle counter (Alphasense OPC-N2) for ambient air monitoring. Atmos. Meas. Tech. 2018, 11, 709–720. [Google Scholar] [CrossRef]
  22. Bílek, J.; Bílek, O.; Maršolek, P.; Buček, P. Ambient Air Quality Measurement with Low-Cost Optical and Electrochemical Sensors: An Evaluation of Continuous Year-Long Operation. Environments 2021, 8, 114. [Google Scholar] [CrossRef]
  23. Dubey, R.; Patra, A.K.; Joshi, J.; Blankenberg, D.; Nazneen. Evaluation of vertical and horizontal distribution of particulate matter near an urban roadway using an unmanned aerial vehicle. Sci. Total Environ. 2022, 836, 155600. [Google Scholar] [CrossRef]
  24. Hagan, D.H.; Kroll, J.H. Assessing the accuracy of low-cost optical particle sensors using a physics-based approach. Atmos. Meas. Tech. 2020, 13, 6343–6355. [Google Scholar] [CrossRef]
  25. Air Quality Portal of the Chief Inspectorate of Environmental Protection— Measurement Data Archives. Available online: https://powietrze.gios.gov.pl/pjp/current/station_details/archive/530?lang=en (accessed on 5 September 2023).
  26. Global Modeling and Assimilation Office. MERRA2 tavg1_2d_aer_Nx: 2d,1-Hourly, Time-Averaged, Single-Level, Assimilation, Aerosol Diagnostics, Version 5.12.4. 2015. Available online: https://disc.gsfc.nasa.gov/datasets/M2T1NXAER_5.12.4/summary (accessed on 31 August 2023).
  27. MERRA-2 Frequently Asked Questions. Available online: https://gmao.gsfc.nasa.gov/reanalysis/MERRA-2/FAQ/ (accessed on 5 September 2023).
  28. Chin, M.; Ginoux, P.; Kinne, S.; Torres, O.; Holben, B.N.; Duncan, B.N.; Martin, R.V.; Logan, J.A.; Higurashi, A.; Nakajima, T. Tropospheric Aerosol Optical Thickness from the GOCART Model and Comparisons with Satellite and Sun Photometer Measurements. J. Atmos. Sci. 2002, 59, 461–483. [Google Scholar] [CrossRef]
  29. Seber, G.; Wild, C. Nonlinear Regression; Wiley Series in Probability and Statistics; Wiley: Hoboken, NJ, USA, 2003. [Google Scholar]
  30. Lundgren, D.A. Atmospheric Aerosol Composition and Concentration as a Function of Particle Size and of Time. J. Air Pollut. Control. Assoc. 1970, 20, 603–608. [Google Scholar] [CrossRef] [PubMed]
  31. Klejnowski, R.; Olędzki, R.; Ostrowski, M. Weatherman Commentary. Available online: https://meteo.icm.edu.pl/komentarze (accessed on 5 September 2023).
  32. Weather Service of the Polish Institute of Meteorology and Water Management—National Research Institute. Available online: http://meteo.imgw.pl/ (accessed on 5 September 2023).
  33. Natural Earth-Free Vector and Raster Map Data. Available online: https://www.naturalearthdata.com/downloads/10m-physical-vectors/10m-ocean/ (accessed on 5 September 2023).
  34. ten Brink, H.; Henzing, B.; Otjes, R.; Weijers, E. Visibility in the Netherlands during New Year’s fireworks: The role of soot and salty aerosol products. Atmos. Environ. 2018, 173, 289–294. [Google Scholar] [CrossRef]
  35. Drewnick, F.; Hings, S.S.; Curtius, J.; Eerdekens, G.; Williams, J. Measurement of fine particulate and gas-phase species during the New Year’s fireworks 2005 in Mainz, Germany. Atmos. Environ. 2006, 40, 4316–4327. [Google Scholar] [CrossRef]
Figure 1. Figure shows one cycle of the ACS 1000 and registered PM by OPC-N3 from which was calculated κ . The change of RH in the ACS 1000 is shown on panel (a), the registered PM 1 , PM 2.5 are shown in panel (b) and in panel (c) PM 10 (after calibration and applying running average of 12 points). Plots (df) presents GF versus RH for one cycle obtained from PM 1 , PM 2.5 , PM 10 respectively, the black line shows fitted curve defined by Equation (4). Pink oval presents a jump in GF indicating DRH point which occured at RH = 84%.
Figure 1. Figure shows one cycle of the ACS 1000 and registered PM by OPC-N3 from which was calculated κ . The change of RH in the ACS 1000 is shown on panel (a), the registered PM 1 , PM 2.5 are shown in panel (b) and in panel (c) PM 10 (after calibration and applying running average of 12 points). Plots (df) presents GF versus RH for one cycle obtained from PM 1 , PM 2.5 , PM 10 respectively, the black line shows fitted curve defined by Equation (4). Pink oval presents a jump in GF indicating DRH point which occured at RH = 84%.
Atmosphere 15 00061 g001
Figure 2. Figure represents steps performed to calculate GF for each cycle. (a) First, the region where RH was up to 5 percentile was calculated (blue region). (b) The PM w e t (t 0 ) was calculated. (c) The Δ PM d r y (t) was calculated as a difference between PM d r y (t) and PM d r y (t 0 ).
Figure 2. Figure represents steps performed to calculate GF for each cycle. (a) First, the region where RH was up to 5 percentile was calculated (blue region). (b) The PM w e t (t 0 ) was calculated. (c) The Δ PM d r y (t) was calculated as a difference between PM d r y (t) and PM d r y (t 0 ).
Atmosphere 15 00061 g002
Figure 3. Daily average κ obtained for months between December and June, for years 2020–2022. Panel (a) calculation done using PM 1 , (b) PM 2.5 and (c) PM 10 . The dashed line presents level κ = 0.05.
Figure 3. Daily average κ obtained for months between December and June, for years 2020–2022. Panel (a) calculation done using PM 1 , (b) PM 2.5 and (c) PM 10 . The dashed line presents level κ = 0.05.
Atmosphere 15 00061 g003
Figure 4. Daily average GF obtained for months between December and June, for years 2020–2022. Panel (a) calculation done using PM 1 , (b) PM 2.5 and (c) PM 10 . The dashed line presents the level of GF = 1.3 .
Figure 4. Daily average GF obtained for months between December and June, for years 2020–2022. Panel (a) calculation done using PM 1 , (b) PM 2.5 and (c) PM 10 . The dashed line presents the level of GF = 1.3 .
Atmosphere 15 00061 g004
Figure 5. Boxplot presenting 3 h averaged mean of κ obtained from PM 1 (a,d), PM 2.5 (b,e), PM 10 (c,f) for winter months upper panels, for spring bottom panels. The central line indicates the median, the ‘x’ symbol indicates the mean.
Figure 5. Boxplot presenting 3 h averaged mean of κ obtained from PM 1 (a,d), PM 2.5 (b,e), PM 10 (c,f) for winter months upper panels, for spring bottom panels. The central line indicates the median, the ‘x’ symbol indicates the mean.
Atmosphere 15 00061 g005
Figure 6. December 2021: The panels represent changes of: (a) PM 2.5 (orange line) and PM 10 (black line) from GIOŚ station located 1.74 km from RTL; (b) κ obtained from ACS 1000 and OPC-N3 data in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 ; (c) effective radius of particles measured by the LAS instrument; (d) GF for RH 85% obtained from MERRA-2; (e) percentage over the sea of air back trajectory obtained from HYSPLIT model, in the panel (f,g) is presented the mass contribution of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SO—sulfate aerosol to the PM 1 and PM 10 respectively obtained from MERRA-2.
Figure 6. December 2021: The panels represent changes of: (a) PM 2.5 (orange line) and PM 10 (black line) from GIOŚ station located 1.74 km from RTL; (b) κ obtained from ACS 1000 and OPC-N3 data in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 ; (c) effective radius of particles measured by the LAS instrument; (d) GF for RH 85% obtained from MERRA-2; (e) percentage over the sea of air back trajectory obtained from HYSPLIT model, in the panel (f,g) is presented the mass contribution of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SO—sulfate aerosol to the PM 1 and PM 10 respectively obtained from MERRA-2.
Atmosphere 15 00061 g006
Figure 7. January 2022: The panels represent changes of: (a) PM 2.5 (orange line) and PM 10 (black line) from GIOŚ station located 1.74 km from RTL; (b) κ obtained from ACS 1000 and OPC-N3 data in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 ; (c) effective radius of particles measured by the LAS instrument; (d) GF for RH 85% obtained from MERRA-2; (e) percentage over the sea of air back trajectory obtained from HYSPLIT model, in the panel (f,g) is presented the mass contribution of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SO—sulfate aerosol to the PM 1 and PM 10 respectively obtained from MERRA-2.
Figure 7. January 2022: The panels represent changes of: (a) PM 2.5 (orange line) and PM 10 (black line) from GIOŚ station located 1.74 km from RTL; (b) κ obtained from ACS 1000 and OPC-N3 data in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 ; (c) effective radius of particles measured by the LAS instrument; (d) GF for RH 85% obtained from MERRA-2; (e) percentage over the sea of air back trajectory obtained from HYSPLIT model, in the panel (f,g) is presented the mass contribution of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SO—sulfate aerosol to the PM 1 and PM 10 respectively obtained from MERRA-2.
Atmosphere 15 00061 g007
Figure 8. Gray lines represents the back trajectories calculated, for 0000UTC and 1200UTC, using HYSPLIT ending at 500 m a.s.l. at RTL. The origin of air (48 h prior) observed over the RTL at Warsaw is marked by the lower dot in the stuck. Each stuck point represents the measurements κ PM 10 for the period closest to the simulation of air flow from this region. Two plots represent, respectively, (a) December 2021 and (b) January 2022.
Figure 8. Gray lines represents the back trajectories calculated, for 0000UTC and 1200UTC, using HYSPLIT ending at 500 m a.s.l. at RTL. The origin of air (48 h prior) observed over the RTL at Warsaw is marked by the lower dot in the stuck. Each stuck point represents the measurements κ PM 10 for the period closest to the simulation of air flow from this region. Two plots represent, respectively, (a) December 2021 and (b) January 2022.
Atmosphere 15 00061 g008
Figure 9. April–May 2020: The panels represent changes of: (a) PM 2.5 (orange line) and PM 10 (black line) of the GIOŚ station located 1.74 km from RTL; (b) κ obtained from the ACS 1000 and OPC-N3 data in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 ; (c) effective radius of particles measured by the LAS instrument; (d) GF for RH 85% obtained from MERRA-2; (e) percentage over the sea of air back trajectory obtained from HYSPLIT model, in the panel (f,g) is presented the mass contribution of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SO—sulfate aerosol to the PM 1 and PM 10 respectively obtained from MERRA-2.
Figure 9. April–May 2020: The panels represent changes of: (a) PM 2.5 (orange line) and PM 10 (black line) of the GIOŚ station located 1.74 km from RTL; (b) κ obtained from the ACS 1000 and OPC-N3 data in RTL for each cycle. Blue line PM 1 , orange line PM 2.5 , black line PM 10 ; (c) effective radius of particles measured by the LAS instrument; (d) GF for RH 85% obtained from MERRA-2; (e) percentage over the sea of air back trajectory obtained from HYSPLIT model, in the panel (f,g) is presented the mass contribution of DU—dust, SS—sea salt, OC—organic carbon, BC—black carbon, SO—sulfate aerosol to the PM 1 and PM 10 respectively obtained from MERRA-2.
Atmosphere 15 00061 g009
Figure 10. Gray lines represents the back trajectories calculated, for 0000UTC and 1200UTC, using HYSPLIT ending at 500 m a.s.l. at RTL. The origin of air (48 h prior) observed over the RTL at Warsaw is marked by the lower dot in the stuck. Each dot on the stuck represents the κ PM 10 measurements for the period closest to the simulation of air flow from this region.
Figure 10. Gray lines represents the back trajectories calculated, for 0000UTC and 1200UTC, using HYSPLIT ending at 500 m a.s.l. at RTL. The origin of air (48 h prior) observed over the RTL at Warsaw is marked by the lower dot in the stuck. Each dot on the stuck represents the κ PM 10 measurements for the period closest to the simulation of air flow from this region.
Atmosphere 15 00061 g010
Figure 11. (a) PM 1 , w e t , PM 1 , d r y and PM 1 , R H = 0 obtained from PM 1 w e t for 29 April–1 May 2020; (b) RH from OPC-N3 versus RH from ACS 1000 for Period nr 1 (spring 2020); (c) RH from ACS 1000 for 29 April–1 May 2020.
Figure 11. (a) PM 1 , w e t , PM 1 , d r y and PM 1 , R H = 0 obtained from PM 1 w e t for 29 April–1 May 2020; (b) RH from OPC-N3 versus RH from ACS 1000 for Period nr 1 (spring 2020); (c) RH from ACS 1000 for 29 April–1 May 2020.
Atmosphere 15 00061 g011
Table 1. GF at different RH for OC, BC, SU, SS [28].
Table 1. GF at different RH for OC, BC, SU, SS [28].
RH [%]0507080909599
GF S U 1.01.41.51.61.81.92.2
GF O C 1.01.21.41.51.61.82.2
GF B C 1.01.01.01.21.31.51.8
GF S S 1.01.61.82.02.42.94.8
Table 2. Periods in which the data about hygroscopicity were collected by ACS 1000 and OPC-N3. When the setting of ACS 1000 was changed, a new calibration of OPC-N3 was performed.
Table 2. Periods in which the data about hygroscopicity were collected by ACS 1000 and OPC-N3. When the setting of ACS 1000 was changed, a new calibration of OPC-N3 was performed.
PeriodCalibration PerformedBeginingEnding
nrfor RH [%] Below:[YYYY.MM.DD][YYYY.MM.DD]
1602020.03.262020.05.23
2552020.12.172021.01.05
3502021.01.062021.01.17
4602021.06.042021.06.13
5502021.12.162022.01.20
6502022.01.202022.02.28
7502022.03.012022.03.15
Table 3. Pearson correlation coefficient between variables listed in first column and row ( κ , GF calculated from PM 1 , PM 2.5 , PM 10 ). Colors divide values by strength of correlation. All correlations have a p-value lower than the assumed statistical significance 5 % .
Table 3. Pearson correlation coefficient between variables listed in first column and row ( κ , GF calculated from PM 1 , PM 2.5 , PM 10 ). Colors divide values by strength of correlation. All correlations have a p-value lower than the assumed statistical significance 5 % .
GF PM 1 1.00
GF PM 2.5 0.781.00
GF PM 10 0.380.631.00
κ PM 1 0.980.760.361.00
κ PM 2.5 0.770.970.610.771.00
κ PM 10 0.390.620.890.380.621.00
GF PM 1 GF PM 2.5 GF PM 10 κ PM 1 κ PM 2.5 κ PM 10
0.9 < |r| ≤ 1.0—very high correlation
 0.7 < |r| ≤ 0.9—high correlation
 0.4 < |r| ≤ 0.7—medium correlation
 0.2 < |r| ≤ 0.4—weak correlation
 0.0 ≤ |r| ≤ 0.2—no correlation
Table 4. Mean values of GF and κ with their standard deviations calculated for PM 1 , PM 2.5 , PM 10 for winter and spring season for years 2020–2022.
Table 4. Mean values of GF and κ with their standard deviations calculated for PM 1 , PM 2.5 , PM 10 for winter and spring season for years 2020–2022.
SeasonYear κ PM 1 κ PM 2.5 κ PM 10
winter (January–March)20200.04 ± 0.040.05 ± 0.050.06 ± 0.07
20210.04 ± 0.030.05 ± 0.070.08 ± 0.13
20220.04 ± 0.050.07 ± 0.080.16 ± 0.20
Over all years0.04 ± 0.040.06 ± 0.070.12 ± 0.17
spring (April–July)20200.14 ± 0.150.15 ± 0.150.15 ± 0.15
20210.11 ± 0.090.09 ± 0.090.09 ± 0.09
2022
Over all years0.13 ± 0.130.13 ± 0.130.13 ± 0.13
SeasonYearGF PM 1 GF PM 2.5 GF PM 10
winter (January–March)20201.16 ± 0.151.19 ± 0.161.22 ± 0.29
20211.14 ± 0.121.21 ± 0.261.35 ± 0.63
20221.16 ± 0.131.30 ± 0.301.71 ± 0.87
Over all years1.16 ± 0.131.25 ± 0.271.51 ± 0.76
spring (April–July)20201.58 ± 0.261.63 ± 0.261.63 ± 0.26
20211.43 ± 0.141.38 ± 0.141.38 ± 0.14
2022
Over all years1.52 ± 0.231.53 ± 0.231.53 ± 0.23
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nurowska, K.; Markowicz, K.M. Determination of Hygroscopic Aerosol Growth Based on the OPC-N3 Counter. Atmosphere 2024, 15, 61. https://doi.org/10.3390/atmos15010061

AMA Style

Nurowska K, Markowicz KM. Determination of Hygroscopic Aerosol Growth Based on the OPC-N3 Counter. Atmosphere. 2024; 15(1):61. https://doi.org/10.3390/atmos15010061

Chicago/Turabian Style

Nurowska, Katarzyna, and Krzysztof M. Markowicz. 2024. "Determination of Hygroscopic Aerosol Growth Based on the OPC-N3 Counter" Atmosphere 15, no. 1: 61. https://doi.org/10.3390/atmos15010061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop