A new methodology for performing long-term source apportionment (SA) using positive matrix factorization (PMF) is presented. The method is implemented within the SoFi Pro software package and uses the multilinear engine (ME-2) as a PMF solver. The technique is applied to a 1-year aerosol chemical speciation monitor (ACSM) dataset from downtown Zurich, Switzerland. The measured organic aerosol mass spectra were analyzed by PMF using a small (14 d) and rolling PMF window to account for the temporal evolution of the sources. The rotational ambiguity is explored and the uncertainties of the PMF solutions were estimated. Factor–tracer correlations for averaged seasonal results from the rolling window analysis are higher than those retrieved from conventional PMF analyses of individual seasons, highlighting the improved performance of the rolling window algorithm for long-term data. In this study four to five factors were tested for every PMF window. Factor profiles for primary organic aerosol from traffic (HOA), cooking (COA) and biomass burning (BBOA) were constrained. Secondary organic aerosol was represented by either the combination of semi-volatile and low-volatility organic aerosol (SV-OOA and LV-OOA, respectively) or by a single OOA when this separation was not robust. This scheme led to roughly 40 000 PMF runs. Full visual inspection of all these PMF runs is unrealistic and is replaced by predefined user-selected criteria, which allow factor sorting and PMF run acceptance/rejection. The selected criteria for traffic (HOA) and BBOA were the correlation with equivalent black carbon from traffic (eBC_tr) and the explained variation of $m/z$ 60, respectively. COA was assessed by the prominence of a lunchtime concentration peak within the diurnal cycle. SV-OOA and LV-OOA were evaluated based on the fractions of $m/z$ 43 and 44 in their respective factor profiles. Seasonal pre-tests revealed a non-continuous separation of OOA into SV-OOA and LV-OOA, in particular during the warm seasons. Therefore, a differentiation between four-factor solutions (HOA, COA, BBOA and OOA) and five-factor solutions (HOA, COA, BBOA, SV-OOA and LV-OOA) was also conducted based on the criterion for SV-OOA. HOA and COA contribute between 0.4–0.7 µg m⁻³ (7.8 %–9.0 %) and 0.7–1.2 µg m⁻³ (12.2 %–15.7 %) on average throughout the year, respectively. BBOA shows a strong yearly cycle with the lowest mean concentrations in summer (0.6 µg m⁻³, 12.0 %), slightly higher mean concentrations during spring and fall (1.0 and 1.5 µg m⁻³, or 15.6 % and 18.6 %, respectively), and the highest mean concentrations during winter (1.9 µg m⁻³, 25.0 %). In summer, OOA is separated into SV-OOA and LV-OOA, with mean concentrations of 1.4 µg m⁻³ (26.5 %) and 2.2 µg m⁻³ (40.3 %), respectively. For the remaining seasons the seasonal concentrations of SV-OOA, LV-OOA and OOA range from 0.3 to 1.1 µg m⁻³ (3.4 %–15.9 %), from 0.6 to 2.2 µg m⁻³ (7.7 %–33.7 %) and from 0.9 to 3.1 µg m⁻³ (13.7 %–39.9 %), respectively. The relative PMF errors modeled for this study for HOA, COA, BBOA, LV-OOA, SV-OOA and OOA are on average ±34 %, ±27 %, ±30 %, ±11 %, ±25 % and ±12 %, respectively.

We collected 1 year of aerosol chemical speciation monitor (ACSM) data in Magadino, a village located in the south of the Swiss Alpine region, one of Switzerland's most polluted areas. We analysed the mass spectra of organic aerosol (OA) by positive matrix factorisation (PMF) using Source Finder Professional (SoFi Pro) to retrieve the origins of OA. Therein, we deployed a rolling algorithm, which is closer to the measurement, to account for the temporal changes in the source profiles. As the first-ever application of rolling PMF with multilinear engine (ME-2) analysis on a yearlong dataset that was collected from a rural site, we resolved two primary OA factors (traffic-related hydrocarbon-like OA (HOA) and biomass burning OA (BBOA)), one mass-to-charge ratio ($m/z$ ) 58-related OA (58-OA) factor, a less oxidised oxygenated OA (LO-OOA) factor, and a more oxidised oxygenated OA (MO-OOA) factor. HOA showed stable contributions to the total OA through the whole year ranging from 8.1 % to 10.1 %, while the contribution of BBOA showed an apparent seasonal variation with a range of 8.3 %–27.4 % (highest during winter, lowest during summer) and a yearly average of 17.1 %. OOA (sum of LO-OOA and MO-OOA) contributed 71.6 % of the OA mass, varying from 62.5 % (in winter) to 78 % (in spring and summer). The 58-OA factor mainly contained nitrogen-related variables which appeared to be pronounced only after the filament switched. However, since the contribution of this factor was insignificant (2.1 %), we did not attempt to interpolate its potential source in this work. The uncertainties (σ) for the modelled OA factors (i.e. rotational uncertainty and statistical variability in the sources) varied from ±4 % (58-OA) to a maximum of ±40 % (LO-OOA). Considering that BBOA and LO-OOA (showing influences of biomass burning in winter) had significant contributions to the total OA mass, we suggest reducing and controlling biomass-burning-related residential heating as a mitigation strategy for better air quality and lower PM levels in this region or similar locations. In Appendix A, we conduct a head-to-head comparison between the conventional seasonal PMF analysis and the rolling mechanism. We find similar or slightly improved results in terms of mass concentrations, correlations with external tracers, and factor profiles of the constrained POA factors. The rolling results show smaller scaled residuals and enhanced correlations between OOA factors and corresponding inorganic salts compared to those of the seasonal solutions, which was most likely because the rolling PMF analysis can capture the temporal variations in the oxidation processes for OOA components. Specifically, the time-dependent factor profiles of MO-OOA and LO-OOA can well explain the temporal viabilities of two main ions for OOA factors, $m/z$ 44 (CO${}_{\mathrm{2}}^{+}$ ) and $m/z$ 43 (mostly C₂H₃O⁺). Therefore, this rolling PMF analysis provides a more realistic source apportionment (SA) solution with time-dependent OA sources. The rolling results also show good agreement with offline Aerodyne aerosol mass spectrometer (AMS) SA results from filter samples, except for in winter. The latter discrepancy is likely because the online measurement can capture the fast oxidation processes of biomass burning sources, in contrast to the 24 h filter samples. This study demonstrates the strengths of the rolling mechanism, provides a comprehensive criterion list for ACSM users to obtain reproducible SA results, and is a role model for similar analyses of such worldwide available data.

Kraków is routinely affected by very high air pollution levels, especially during the winter months. Although a lot of effort has been made to characterize ambient aerosol, there is a lack of online and long-term measurements of non-refractory aerosol. Our measurements at the AGH University of Science and Technology provide the online long-term chemical composition of ambient submicron particulate matter (PM₁) between January 2018 and April 2019. Here we report the chemical characterization of non-refractory submicron aerosol and source apportionment of the organic fraction by positive matrix factorization (PMF). In contrast to other long-term source apportionment studies, we let a small PMF window roll over the dataset instead of performing PMF over the full dataset or on separate seasons. In this way, the seasonal variation in the source profiles can be captured. The uncertainties in the PMF solutions are addressed by the bootstrap resampling strategy and the random a-value approach for constrained factors. We observe clear seasonal patterns in the concentration and composition of PM₁, with high concentrations during the winter months and lower concentrations during the summer months. Organics are the dominant species throughout the campaign. Five organic aerosol (OA) factors are resolved, of which three are of a primary nature (hydrocarbon-like OA (HOA), biomass burning OA (BBOA) and coal combustion OA (CCOA)) and two are of a secondary nature (more oxidized oxygenated OA (MO-OOA) and less oxidized oxygenated OA (LO-OOA)). While HOA contributes on average 8.6 % ± 2.3 % throughout the campaign, the solid-fuel-combustion-related BBOA and CCOA show a clear seasonal trend with average contributions of 10.4 % ± 2.7 % and 14.1 %, ±2.1 %, respectively. Not only BBOA but also CCOA is associated with residential heating because of the pronounced yearly cycle where the highest contributions are observed during wintertime. Throughout the campaign, the OOA can be separated into MO-OOA and LO-OOA with average contributions of 38.4 % ± 8.4 % and 28.5 % ± 11.2 %, respectively.

The Station for Measuring Ecosystem–Atmosphere Relations (SMEAR) II, located within the boreal forest of Finland, is a unique station in the world due to the wide range of long-term measurements tracking the Earth–atmosphere interface. In this study, we characterize the composition of organic aerosol (OA) at SMEAR II by quantifying its driving constituents. We utilize a multi-year data set of OA mass spectra measured in situ with an Aerosol Chemical Speciation Monitor (ACSM) at the station. To our knowledge, this mass spectral time series is the longest of its kind published to date. Similarly to other previously reported efforts in OA source apportionment from multi-seasonal or multi-annual data sets, we approached the OA characterization challenge through positive matrix factorization (PMF) using a rolling window approach. However, the existing methods for extracting minor OA components were found to be insufficient for our rather remote site. To overcome this issue, we tested a new statistical analysis framework. This included unsupervised feature extraction and classification stages to explore a large number of unconstrained PMF runs conducted on the measured OA mass spectra. Anchored by these results, we finally constructed a relaxed chemical mass balance (CMB) run that resolved different OA components from our observations. The presented combination of statistical tools provided a data-driven analysis methodology, which in our case achieved robust solutions with minimal subjectivity. Following the extensive statistical analyses, we were able to divide the 2012–2019 SMEAR II OA data (mass concentration interquartile range (IQR): 0.7, 1.3, and 2.6 µg m⁻³) into three sub-categories – low-volatility oxygenated OA (LV-OOA), semi-volatile oxygenated OA (SV-OOA), and primary OA (POA) – proving that the tested methodology was able to provide results consistent with literature. LV-OOA was the most dominant OA type (organic mass fraction IQR: 49 %, 62 %, and 73 %). The seasonal cycle of LV-OOA was bimodal, with peaks both in summer and in February. We associated the wintertime LV-OOA with anthropogenic sources and assumed biogenic influence in LV-OOA formation in summer. Through a brief trajectory analysis, we estimated summertime natural LV-OOA formation of tens of ng m⁻³ h⁻¹ over the boreal forest. SV-OOA was the second highest contributor to OA mass (organic mass fraction IQR: 19 %, 31 %, and 43 %). Due to SV-OOA's clear peak in summer, we estimate biogenic processes as the main drivers in its formation. Unlike for LV-OOA, the highest SV-OOA concentrations were detected in stable summertime nocturnal surface layers. Two nearby sawmills also played a significant role in SV-OOA production as also exemplified by previous studies at SMEAR II. POA, taken as a mix of two different OA types reported previously, hydrocarbon-like OA (HOA) and biomass burning OA (BBOA), made up a minimal OA mass fraction (IQR: 2 %, 6 %, and 13 %). Notably, the quantification of POA at SMEAR II using ACSM data was not possible following existing rolling PMF methodologies. Both POA organic mass fraction and mass concentration peaked in winter. Its appearance at SMEAR II was linked to strong southerly winds. Similar wind direction and speed dependence was not observed among other OA types. The high wind speeds probably enabled the POA transport to SMEAR II from faraway sources in a relatively fresh state. In the event of slower wind speeds, POA likely evaporated and/or aged into oxidized organic aerosol before detection. The POA organic mass fraction was significantly lower than reported by aerosol mass spectrometer (AMS) measurements 2 to 4 years prior to the ACSM measurements. While the co-located long-term measurements of black carbon supported the hypothesis of higher POA loadings prior to year 2012, it is also possible that short-term (POA) pollution plumes were averaged out due to the slow time resolution of the ACSM combined with the further 3 h data averaging needed to ensure good signal-to-noise ratios (SNRs). Despite the length of the ACSM data set, we did not focus on quantifying long-term trends of POA (nor other components) due to the high sensitivity of OA composition to meteorological anomalies, the occurrence of which is likely not normally distributed over the 8-year measurement period. Due to the unique and realistic seasonal cycles and meteorology dependences of the independent OA subtypes complemented by the reasonably low degree of unexplained OA variability, we believe that the presented data analysis approach performs well. Therefore, we hope that these results encourage also other researchers possessing several-year-long time series of similar data to tackle the data analysis via similar semi- or unsupervised machine-learning approaches. This way the presented method could be further optimized and its usability explored and evaluated also in other environments.

In this study, we investigate the occurrence of primary biological aerosol particles (PBAP) over all sectors of the Southern Ocean (SO) based on a 90‐day data set collected during the Antarctic Circumnavigation Expedition (ACE) in austral summer 2016–2017. Super‐micrometer PBAP (1–16 µm diameter) were measured by a wide band integrated bioaerosol sensor (WIBS‐4). Low (3σ) and high (9σ) fluorescence thresholds are used to obtain statistics on fluorescent and hyper‐fluorescent PBAP, respectively. Our focus is on data obtained over the pristine ocean, that is, more than 200 km away from land. The results indicate that (hyper‐)fluorescent PBAP are correlated to atmospheric variables associated with sea spray aerosol (SSA) particles (wind speed, total super‐micrometer aerosol number concentration, chloride and sodium concentrations). This suggests that a main source of PBAP over the SO is SSA. The median percentage contribution of fluorescent and hyper‐fluorescent PBAP to super‐micrometer SSA was 1.6% and 0.13%, respectively. We demonstrate that the fraction of (hyper‐)fluorescent PBAP to total super‐micrometer particles positively correlates with concentrations of bacteria and several taxa of pythoplankton measured in seawater, indicating that marine biota concentrations modulate the PBAP source flux. We investigate the fluorescent properties of (hyper‐)fluorescent PBAP for several events that occurred near land masses. We find that the fluorescence signal characteristics of particles near land is much more variable than over the pristine ocean. We conclude that the source and concentration of fluorescent PBAP over the open ocean is similar across all sampled sectors of the SO.