Six-year source apportionment of submicron organic aerosols from near-1 continuous highly time-resolved measurements at SIRTA ( Paris area , France ) 2

17 Organic aerosol (OA) particles are recognized as key factors influencing air quality and climate 18 change. However, highly-time resolved long-term characterizations of their composition and 19 sources in ambient air are still very limited due to challenging continuous observations. Here, 20 we present an analysis of long-term variability of submicron OA using the combination of 21 aerosol chemical speciation monitor (ACSM) and multi-wavelength aethalometer from 22 November 2011 to March 2018 at a peri-urban background site of the Paris region (France). 23 Source apportionment of OA was achieved via partially constrained positive matrix 24 factorization (PMF) using the multilinear engine (ME-2). Two primary OA (POA) and two 25 oxygenated OA (OOA) factors were identified and quantified over the entire studied period. 26 POA factors were designated as hydrocarbon-like OA (HOA) and biomass burning OA (BBOA). 27 The latter factor presented a significant seasonality with higher concentrations in winter with 28

Abstract. Organic aerosol (OA) particles are recognized as key factors influencing air quality and climate change. However, highly time-resolved long-term characterizations of their composition and sources in ambient air are still very limited due to challenging continuous observations. Here, we present an analysis of long-term variability of submicron OA using the combination of an aerosol chemical speciation monitor (ACSM) and a multiwavelength Aethalometer from November 2011 to March 2018 at a peri-urban background site of the Paris region (France). Source apportionment of OA was achieved via partially constrained positive matrix factorization (PMF) using the multilinear engine (ME-2). Two primary OA (POA) and two oxygenated OA (OOA) factors were identified and quantified over the entire studied period. POA factors were designated as hydrocarbon-like OA (HOA) and biomass burning OA (BBOA). The latter factor presented a significant seasonality with higher concentrations in winter with significant monthly contributions to OA (18 %-33 %) due to enhanced residential wood burning emissions. HOA mainly originated from traffic emissions but was also influenced by biomass burning in cold periods. OOA factors were distinguished between their less-and more-oxidized fractions (LO-OOA and MO-OOA, respectively). These factors presented distinct seasonal patterns, associated with different atmospheric formation pathways. A pronounced increase in LO-OOA concentrations and contributions (50 %-66 %) was observed in summer, which may be mainly explained by secondary OA (SOA) formation processes involving biogenic gaseous precursors. Conversely, high concentrations and OA contributions (32 %-62 %) of MO-OOA during winter and spring seasons were partly associated with anthropogenic emissions and/or long-range transport from northeastern Europe. The contribution of the different OA factors as a function of OA mass loading highlighted the dominant roles of POA during pollution episodes in fall and winter and of SOA for highest springtime and summertime OA concentrations. Finally, long-term trend analyses indicated a decreasing feature (of about −175 ng m −3 yr −1 ) for MO-OOA, very limited or insignificant decreasing trends for primary anthropogenic carbonaceous aerosols (BBOA and HOA, along with the fossil-fuel and biomass-burning black carbon components) and no statistically significant trend for LO-OOA over the 6-year investigated period.

Introduction
Organic aerosol (OA) particles account for a large mass fraction of submicron particulate matter (PM 1 ) in the atmosphere  and play a key role in regional air pollution and climate (Boucher et al., 2013). Primary OA (POA) originates from direct emissions of primary sources (e.g., fossil-fuel and biomass combustion). Secondary OA (SOA) is formed from atmospheric oxidation processes of gas precursors, i.e., volatile organic compounds (VOCs) (Kroll and Seinfeld, 2008;Hallquist et al., 2009;Nozière et al., 2015). Some typical SOA formation processes in the atmosphere, such as photochemistry (Xu et al., 2017), aqueous-phase oxidation (Gilardoni et al., 2016) and heterogeneous reactions (Xu et al., 2015), are observed. Due to their multiplicity and complexity, these various sources and physicochemical mechanisms remain poorly documented and understood.
Although numerous time-limited field campaigns allowed us to greatly improve our knowledge of OA properties in the last decade (e.g., Jimenez et al., 2009;Lanz et al., 2010;Zhang et al., 2011;Shrivastava et al., 2017;Li et al., 2017;Srivastava et al., 2018a, and references therein), similar studies performed on a long-term scale remain scarce and particularly challenging (Fröhlich et al., 2015a;Schlag et al., 2016;Bozzetti et al., 2017;Daellenbach et al., 2017;Sun et al., 2018). Long-term observations with high temporal resolution and source apportionment of OA are nevertheless necessary to better quantify the contribution of airborne OA particles to air quality and to set up scientifically sound emission control strategies. They can also contribute to a better understanding of the atmospheric fate of OA and reduce uncertainties associated with its (in)direct radiative forcing.
Online aerosol characterization techniques, such as aerosol mass spectrometry (AMS), have demonstrated their capacity to improve our knowledge of key aerosol chemical components -such as OA -by providing highly timeresolved mass spectral data for the nonrefractory PM 1 fraction (NR-PM 1 ) (Jayne et al., 2000;Canagaratna et al., 2007). Using receptor model approaches, especially positive matrix factorization (PMF) (Paatero and Tapper, 1994), OA measured by AMS techniques can be further portioned into various source factors using statistic models Zhang et al., 2011). For example, hydrocarbon-like OA (HOA) is frequently identified within urban environments and attributed to primary emissions from fuel consumption Jimenez et al., 2009), while biomass burning OA (BBOA) is often resolved specifically during cold seasons or within wild fire plumes Lanz et al., 2010;Zhou et al., 2017). Oxygenated OA (OOA), commonly considered a surrogate for SOA, is ubiquitously observed in urban, suburban and remote environments Srivastava et al., 2018a;Crippa et al., 2014). OOA can be further separated into different fractions, being for instance classified according to its atmospheric aging described as more oxidized (MO-OOA) or less oxidized (LO-OOA) compared to each other Ng et al., 2011a;Sun et al., 2018). Different OOA factors can also be identified as relevant to various sources of SOA precursors, such as anthropogenic activities (e.g., traffic and biomass burning emissions) (Gilardoni et al., 2016;Gentner et al., 2017) and biogenic emissions (e.g., isoprene and monoterpenes) (Xu et al., 2015;Zhang et al., 2018;Freney et al., 2018) in specific regions and/or seasons. Such source apportionment has the potential to assess the efficiency of air pollution mitigation by current emission control strategies.
Based on being better suited for long-term monitoring applications due to lower cost and easier maintenance than AMS, an aerosol chemical speciation monitor (ACSM) has been designed to provide continuous measurements of the main nonrefractory chemical species within submicron aerosols (Ng et al., 2011b). As for the AMS, OA mass spectra obtained by the ACSM can be used in PMF analysis for quantification of OA sources (e.g., Sun et al., 2012;Fröhlich et al., 2015b;Zhang et al., 2015). So far, several time-extended OA source apportionment studies have been reported based on ACSM measurements at various sites Fröhlich et al., 2015a;Schlag et al., 2016;Reyes-Villegas et al., 2016;Rattanavaraha et al., 2017;Sun et al., 2018). However, these studies have been limited to periods of up to 2 years durations.
The longest ACSM time series recorded so far (from end of 2011 onwards) is used here to investigate OA sources at a regional background site of the Paris region (France), which is one of the largest urbanized regions in Europe. It has already been demonstrated that OA plays a dominant role in controlling atmospheric pollution in this region (Bressi et al., 2013;Petit et al., 2015). Furthermore, time-limited (typically, 1-2 months) measurement campaigns demonstrated that primary fine aerosols are mainly influenced there by traffic emissions all over the year and residential wood burning during cold seasons, while secondary aerosols originate from both local production and regional transport (Sciare et al., 2011;Crippa et al., 2013a, b;Petit et al., 2014;Srivastava et al., 2018b). Furthermore, such a background site can be considered representative of air quality at a regional scale, including neighboring northwestern countries (Bressi et al., 2013(Bressi et al., , 2014. In the present study, main OA factors were identified and quantified from seasonal PMF analyses (total 25 seasons) on the 6-year ACSM datasets, with the objective of understanding sources and long-term temporal trends of these factors. In this respect, sporadic and/or minor OA sources were not accounted for in this study. The seasonal variations, weekly and diel cycles, and the long-term temporal trends of the major OA factors were investigated. The relative contributions of the various POA and SOA fractions were also plotted as a function of total submicron OA loadings with the objective to better identify the main OA sources responsible for regional pollution episodes. Finally, the geographical origins of high loadings of SOA factors were investigated using air mass back-trajectory analyses.

Sampling site and instrumentation
Long-term submicron aerosol online measurements used in this study were performed from 1 November 2011 to  Pandolfi et al., 2018). It is located 25 km southwest of Paris city center and is considered representative of the background air quality of the Paris region (Haeffelin et al., 2005;Petit et al., 2015). More information about the environmental features of the sampling site area is given in the introduction above. Major submicron aerosol chemical species, i.e., OA, nitrate, sulfate, ammonium, and chloride, were measured using a quadrupole ACSM. These measurements were performed continuously, always using the same instrument. Over the whole investigated period, the data capture is of about 87 %, and missing data correspond to two field campaigns performed elsewhere (in fall 2012 and March 2013) and to a few technical breakdown and maintenance periods. Briefly, fine aerosols are sampled into the ACSM system through a 100 µm diameter critical orifice mounted at the inlet of the PM 1 aerodynamic lens (Liu et al., 2007;Ng et al., 2011b). Then, submicron aerosol particles are impacted and vaporized at a temperature (T ) of about 600 • C and detected using electron impact (70 eV) ionization mass spectrometry. The ACSM was operated at a time resolution of about 30 min with a scan rate of 0.2 s amu −1 from m/z 12 to 150 amu (atomic mass unit). Coarse particles were removed upstream using an URG Corporation cyclone separator (with the size cutoff diameter of 2.5 µm). Calibrations of the detector response factor were performed regularly (typically every 6 months) using ammonium nitrate solutions (Ng et al., 2011b;Freney et al., 2019). The 1.4 default value was used for the OA relative ion efficiency for the whole dataset . The composition-dependent collection efficiency correction recommended by Middlebrook et al. (2012) was applied to the whole ACSM data used here. The accuracy of these ACSM measurements and the overall good working conditions of the instrument were verified through the participation in the ACTRIS ACSM intercomparison exercises that took place at SIRTA in November-December 2013 Fröhlich et al., 2015b) and March-April 2016 (Freney et al., 2019, and Fig. S1 in the Supplement).
Co-located multiwavelength Aethalometer (Magee Scientific) datasets were also available for the purpose of the study, providing complementary information on equivalent black carbon (eBC) concentrations and sources. Two Aethalometers were used successively: from November 2011 to February 2013 (AE31 model) and then from March 2013 to March 2018 (AE33 model). Both instruments measure aerosol light attenuation at seven wavelengths, i.e., 370, 470, 520, 590, 660, 880 and 950 nm. The detailed descriptions of the AE31 operation at SIRTA and Aethalometer data analysis can be found in Petit et al. (2015). The AE33 is an advanced Aethalometer version, which allows for better assessment and compensation of the filter-loading effect using two simultaneous light attenuation measurements performed at different rates of particle accumulation onto the filter tape (Drinovec et al., 2015(Drinovec et al., , 2017). An excellent agreement (R = 0.94, slope = 1.006 ± 0.006) between AE31 and AE33 for measuring eBC mass concentrations has been demonstrated by Drinovec et al. (2015), suggesting a negligible influence of measurement uncertainties between the two instruments on quantification of eBC concentrations. In this work, the mass concentration of eBC was estimated from an attenuation measurement performed at 880 nm as described by  and Zhang et al. (2018). A correction factor of 1.64 was applied to raw absorption data delivered by the instrument as recommended within the ACTRIS network (Zanatta et al., 2016). Furthermore, eBC could be discriminated between its two main combustion sources, i.e., fossil-fuel combustion (eBC ff ) and wood burning emissions (eBC wb ) using the Aethalometer model (Sandradewi et al., 2008;Sciare et al., 2011;Drinovec et al., 2015). For these calculations, absorption Angström exponent (AAE) valuesin the wavelength range 470-950 nm -of 0.9 and 1.7 for eBC ff and eBC wb , respectively, were estimated by the AAE distributions observed in this study (see Fig. S2). Meanwhile, those values are in a good agreement with the recommended values by Zotter et al. (2017).
In addition to ACSM and AE33 measurements, co-located offline analyses were performed from daily (24 h) PM 2.5 filter samples, collected and analyzed for their elemental and organic carbon contents (EC and OC, respectively) following the ACTRIS recommendations (Zanatta et al., 2017;Zhang et al., 2018). Briefly, filters were collected using a lowvolume sampler (Partisol model 2025; Thermo Scientific) equipped upstream with a VOC denuder system. Mass concentrations of OC and EC from August 2012 to March 2018 were then quantified using a Sunset Lab OC/EC analyzer implemented with the EUSAAR-2 thermal-optical protocol (Cavalli et al., 2010). As shown in Fig. S3, good agreements are obtained between eBC and EC measurements (r 2 = 0.79, slope = 0.94; N = 1185) and between OA and OC measurements (r 2 = 0.68). The slope of 2.14 obtained between submicron OA measured by the ACSM and PM 2.5 OC filterbased methods corresponded to the higher range of values generally observed at (sub)urban background sites -typically 1.6-2.2 (e.g., Bae et al., 2006;Aiken et al., 2008;Favez et al., 2010;Sun et al., 2011;Canagaratna et al., 2015, and references therein) -and may be partly explained by the fact that the filter sampling setup was designed to minimize positive sampling artifacts but does not prevent negative ones. Results obtained from these comparisons with filter-based measurements supported the validity of the datasets used in the present study.
Co-located measurements of nitric oxide (NO) and nitrogen dioxide (NO 2 ) were performed with a NO 2 /NO/NO x analyzer (model T200UP, Teledyne API, USA). Data mea-surements were used to further constrain traffic-related OA sources. The meteorological parameters, including temperature (T ), relative humidity (RH), wind speed (WS), wind direction (WD), boundary layer height (BLH) and precipitation, were obtained from the main SIRTA ground-based meteorological station (located at about 4 km northeast of the aerosol monitoring site). It should be noted that the BLH data were achieved in combining a diagnostic of the surface stability from high-frequency sonic anemometer measurements and light detection and ranging (lidar) attenuated backscatter gradients from aerosols and clouds (Pettie et al., 2015;Dupont et al., 2016).

PMF analysis
The OA source apportionment was performed using the PMF algorithm (Paatero and Tapper, 1994). Organic concentration and error matrices were exported from the ACSM Local software (v 1.5.11.2). Only m/z values ranging from 13 to 100 were applied in the PMF analysis due to larger uncertainties for larger m/z ions and large interferences of naphthalene (m/z 128) signals (Sun et al., 2012). Down-weighting of the m/z 44 group of ions for the PMF model analysis was performed following procedures implemented in the ACSM Local software and following the data treatment strategy proposed by Ulbrich et al. (2009).
When using PMF, it may be difficult to distinguish between factors with similar spectral profiles, especially for ACSM datasets, which are associated with larger uncertainties compared to AMS (Sun et al., 2012;Zhang et al., 2015;Fröhlich et al., 2015b). The source finder (SoFi) toolkit, implemented with the ME-2 solver , has been developed by  to better address this limitation. SoFi provides robust functions which allow us to constrain chosen factor profiles and/or time series. In particular, the a-value approach makes use of range-defining scalar values (with a values ranging from 0 to 1) in order to better elucidate specific PMF factor profile(s) with a chosen degree of freedom; the higher the a value, the less constrained the OA profile is . In the present work, this a-value approach has been used to constrain profiles of POA factors. Some previous studies have already been performed at SIRTA using high-resolution time-of-flight AMS (HR-ToF-AMS) along with PMF analysis during short-time campaigns (typically around 3-4 weeks), leading to the identification of HOA, BBOA, and a cooking OA (COA) factor (Crippa et al., 2013a, b;Fröhlich et al., 2015b). Mass spectra obtained from these studies were used here as references to constrain POA factors because of the prior knowledge of source information as constraint. Conversely, mass spectral profiles of possible OOA factors were left unconstrained. It should be noted that Crippa et al. (2013a) resolved up to three different types of OOA factors and/or a marine OA (MOA) factor when combining HR-ToF-AMS and proton-transferreaction mass spectrometer (PTR-MS) datasets obtained during a summer and a winter campaign at SIRTA.
OOA factor profiles may differ with time, notably due to seasonal variations of several parameters such as meteorological conditions, photochemistry, atmospheric lifetime, air mass origin and/or gaseous precursor origins. In order to better account for such variability, individual PMF analyses were performed on a 3-month basis, i.e., winter (December-January-February), spring (March-April-May), summer (June-July-August), and fall (September-October-November), with a total number of 25 different PMF runs (7 for winters and 6 for each of the other seasons). November 2011 and March 2018 data were included in the winter 2011-2012 and winter 2017-2018 analyses, respectively.
To evaluate the influence of the chosen temporal PMF window (i.e., time duration of data used in ME-2 runs) on the seasonal ME-2 model results, different time frames (i.e., 15, 30, 60, 90 and 120 d) were tested. As shown in Fig. S4 (with 2017 winter data as an example), the excellent consistency of those results from different scenarios suggests a very limited influence of the given PMF windows on determining the outputs of ME-2 analyses. To better assess the variations in primary and secondary OA in different seasons over the 6-year period and to allow for some degrees of freedom within the model runs, the main OA factors, including both POA factors (HOA and BBOA) and two SOA factors (a less oxidized OOA (LO-OOA) and a more oxidized OOA (MO-OOA)), were calculated as the average of 50 convergent ME-2 runs with the a values varying from 0 to 0.4. Moreover, results obtained with an a value of 0.2 were also compared to these results for sensitivity analyses (Fig. S5). The diagnostics of the final OA-factor solution are further discussed in Sect. 4.1.

Influence of biogenic SOA
Biogenic SOA (BSOA) might have a significant influence on OA loadings in mid-latitude regions during summertime and be further apportioned using AMS techniques (e.g., Leaitch et al., 2011;Canonaco et al., 2015). For that reason, the influence of this biogenic OA source was specifically investigated in the present study. To do so, BSOA derived from terpene emissions (BSOA t ) was taken as a surrogate for total BSOA, and the temperature dependence of the BSOA t formation process yield during summertime was simulated using a terpene emission model Schurgers et al., 2009;Leaitch et al., 2011, and references therein). This model is designed to quantify biogenic emissions over global and regional scales. The emission rate (γ ) is estimated by an exponential curve function (Eq. 1), which is describing the relation between terpene emission rate and leaf temperature. As we assumed changes in leaf temperature to be the same as ambient temperature, this could then result, in part, in uncertainties for the model calculation. In addi-tion, this temperature-dependent model reflects vapor pressure changes caused by temperature; however, changes in vapor pressure due to changes in the concentrations in the storage pool of terpene are not covered by the model (Schurgers et al., 2009). Therefore, this emission model is useful to simulate the short-term emissions because of temperature changes (Schurgers et al., 2009).
where γ 0 stands for the emission rate (µg g −1 h −1 ) at standard conditions, and β is an empirical constant chosen here to be equal to 0.09 K −1 (Schurgers et al., 2009;Leaitch et al., 2011). As reported by previous studies, biogenic terpene emissions could be a major source of such a PMF LO-OOA factor observed during summertime in western Europe (e.g., Canonaco et al., 2015;Daellenbach et al., 2017;Daellenbach et al., 2019). Given that BSOA t was assumed to be mainly included in the LO-OOA fraction in the present work, BSOA t estimated concentrations were compared to LO-OOA concentration data points corresponding to the daytime maximum temperature (at approximately 16:00-17:00 LT, local time) in summer. Assuming that LO-OOA could actually be mostly composed of BSOA t during this period of the day and following the procedure described by Leaitch et al. (2011), the daily mass concentrations of BSOA t were estimated as follows: where T min corresponds to the lowest daily maximum temperature observed across the investigated summer seasons (i.e., 12 • C ±1 • C) and LO-OOA (observed at T min) corresponds to the mean LO-OOA concentration obtained for these data points (0.7 ± 0.3 µg m −3 , N = 17).

Trend analysis
The multiyear trends of monthly mean OA factors and total OA, as well as other chemical components (including eBC wb , eBC ff , nitrate, sulfate and total PM 1 ), were analyzed using the Mann-Kendall (MK) trend test (Mann, 1945). The trend slope was calculated using Theil-Sen estimator (Sen, 1968). Before performing the MK trend test, the normality and seasonality of the OA factors were examined, respectively. The normality of the mass concentrations of the OA factors was examined by the Shapiro-Wilk normality test (Shapiro and Wilk, 1965). As a result of the Shapiro-Wilk normality test, all datasets of the mass concentrations of the four OA factors were not normally distributed cases. The MK test associated with Sen's estimator of slope is insensitive to outliers, while it is not appropriate for the chosen dataset with significant seasonality. The Kruskal-Wallis test (Kruskal and Wallis, 1952) was performed to evaluate the seasonality of monthly average datasets at the 5 % significance level. If the seasonality of the data is insignificant, the MK test was used for the trend analysis, while the seasonal MK test was then applied for the data with significant seasonality. In addition, to further compare the differences between the MK test and the seasonal MK test in our trend analysis, both methods have been applied for all datasets (see Table S1 in the Supplement). The trend computation was performed here using a R trend package (Pohlert, 2018).

Air mass back-trajectory analysis
The HYbrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) (Draxler and Rolph, 2003;Stein et al., 2015) was applied to calculate 72 h back trajectories hourly arriving at SIRTA at a height of 100 m above ground level, based on the Global Data Assimilation System (GDAS) meteorological data. The potential source contribution function model (PSCF) (Polissar et al., 1999) was used in this study to investigate the potential source origins that may contribute to high concentrations of OA factors at SIRTA. This analysis was achieved with a resolution of 0.2 • × 0.2 • for each grid cell, using the ZeFir toolkit . The probability function for a given grid cell (i, j ), where i stands for the latitude and j for the longitude, is related to observed concentrations that are higher than a threshold value, which is defined by Eq. (3): where m ij is the total number of selected trajectory endpoints (i, j ) associated with receptor concentrations of PMF factors higher than the threshold value, and n ij is the total number of back-trajectory endpoints at each grid cell (i, j ). The 75th percentile of each OA factor during the entire study was used as the threshold value to calculate m ij . To reduce uncertainty caused by small n ij values for the PSCF modeling, an arbitrary weighting function (w ij ) was applied (Waked et al., 2014). To minimize the influence of some trajectories on the possible pathways of air mass transport, observed data points associated with low-wind-speed conditions (WS < 4 m s −1 ) were filtered out. In addition, observed data points at SIRTA during the period with any hourly precipitation events (precipitation > 0 mm) were removed to reduce influence of wet deposition on ambient aerosol concentrations.  (Crippa et al., 2013a, b;Freutel et al., 2013;Petit et al., 2014;Fröhlich et al., 2015b). Another POA source, i.e., COA, has also been identified using HR-ToF-AMS during previous campaigns in the Paris region (Crippa et al., 2013a, b;Fröhlich et al., 2015b). However, the distinction between COA and HOA factors based solely on ACSM measurements remains challenging due to highly similar mass spectra and uncertainties associated with the ACSM low mass spectral resolution Fröhlich et al., 2015b).
To better assess a potential role of COA in our source apportionment study, several ME-2 runs were conducted constraining either three POA factors (HOA, BBOA, COA) or two (HOA, BBOA). In these tests, POA reference mass spectra determined by Results obtained from these preliminary individual PMF runs showed very good consistency between them with two unconstrained OOA factors -MO-OOA and LO-OOA -always appearing in the 4-factor (with constrained HOA and BBOA factor) and 5-factor (with constrained HOA, BBOA and COA factor) solutions. Conversely, 3-and 6-factor PMF analyses generally led to unsatisfactory solutions.
Figures 1 and S6 present results obtained for the 4-and 5factor solutions, respectively, for the winter 2017-2018 period, taken here as an example. In both cases, mass spectra are in good agreement with those reported in the literature. However, the COA and BBOA factors are displaying very similar diel patterns, leading to surprisingly good correlations between the two factors (Fig. S7). In order to further evaluate possible COA contribution at SIRTA, we applied an m/z-tracer algorithm (Mohr et al., 2012), trying to identify pure cooking aerosol signals. The distribution of the estimated COA signals is centered at about 0, as illustrated by Fig. S8. This could be probably explained by the very little pure cooking influence that could not be quantified by the lower resolution quadrupole ACSM than AMS, which is logically in agreement with negligible cooking source at the sampling site area nearby. Altogether, it could then be concluded that the constrained COA-like aerosols at SIRTA were primarily linked with wood burning emissions, while pure cooking aerosols were probably present loadings that were too low to be properly quantified within the present study. This assumption is consistent with conclusions drawn by other studies performed at SIRTA, e.g., based on an online (ACSM) dataset ) and a combining PMF method using online (ACSM) and offline (4 h filter sampling) Stick lines indicate average values over all selected ME-2 runs. Averaged a values for the constrained factors during the ME-2 runs are also shown. In (b), the upper and lower boundaries of boxes indicate the 75th and 25th percentiles, respectively; the vertical lines within the box correspond to median values; the whiskers above and below boxes refer to 95th and 10th percentiles, respectively; and solid colored lines represent mean values.
Therefore, the 4-factor solution, including two constrained POA factors (BBOA and HOA) and two unconstrained factors, was chosen here as the "best estimate" for the PMF runs performed over the long-term dataset. A total of 25 seasonal and individual PMF analyses were then conducted using a similar procedure. The seasonal OOA factor mass spectra are presented in Fig. S9, showing high seasonal consistency for each OA factor. Moreover, as shown in Fig. S10, the distribution of residuals derived from these 4-factor solution ME-2 runs was sharply centered at around 0, suggesting an insignificance of possible unresolved OA factor(s). The correlations of OA factors with their tracers were examined to globally evaluate the 4-factor PMF solution (see Figs. S11 and S12). As shown in Fig. S11a, HOA is correlated well (r 2 = 0.54) with NO x , a common tracer of primary combustion sources (e.g., traffic emissions). While HOA shows a relatively weaker correlation (r 2 = 0.33) with eBC ff (Fig. S11b), this could be explained by two possible reasons: (i) uncertainties of the Aethalometer model that could not be evaluated by the present study and (ii) the HOA factor here could not be reprehensive for pure fossil-fuel combustion POA. BBOA presents an overall good correlation (r 2 = 0.50) with eBC wb (Fig. S11c), suggesting important influence of wood burning emissions on this factor production. Based on the filter-based dataset, primary OC (POC) and secondary OC (SOC) were calculated using a method of OC-to-EC mass ratio (see Fig. S12). Overall, POA (sum of HOA and BBOA) versus POC (r 2 = 0.47) and SOA (sum of LO-OOA and MO-OOA) versus SOC (r 2 = 0.38) have acceptable correlations during the entire filter measurement period. Thus, all of these comparison results could additionally support our "best estimation" for selecting such 4-factor PMF solution across the entire period.

Source attribution
BBOA mass spectra are quite constant throughout the seasons and present characteristic peaks at m/z 29, 60 and 73, indicative of biomass burning combustion (Fig. S9). As shown in Fig. 2a, BBOA diel cycles displayed wellmarked patterns with strong nighttime maxima. This confirms the predominance of residential wood burning activities on BBOA concentrations at SIRTA and in the Paris region, as already shown previously (e.g., Favez et al., 2009;Sciare et al., 2011;Crippa et al., 2013b;Petit et al., 2014). As expected, BBOA diel cycles are similar to the ones obtained for eBC wb , except for small eBC wb morning peaks that were not observed for BBOA (possibly due to uncertainties of the Aethalometer model). Interestingly, the higher concentrations of eBC wb and BBOA were observed on Saturday and Sunday, reflecting the weekend effect likely due to enhanced residential wood burning emissions.
Compared to BBOA, HOA shows a more complex weekly diel pattern (Fig. 2b). Its pattern is generally similar to eBC ff and NO x (both being considered here as markers for traffic emissions). HOA presents two peaks during the working day, one in the morning and another in the evening. Morning peaks, occurring during traffic rush hour, are clearly indicative of road transport contributions, confirming HOA as a proxy for traffic emissions. However, HOA evening peaks occur globally later than eBC ff and NO x ones (21:00-22:00 LT vs. 19:00 LT , respectively) and much lower ratios are observed between HOA and eBC ff in the morning than in the evening. This might be partly explained by (i) higher eBC traffic emission factor in the morning and/or (ii) impacts of residential heating sources, e.g., wood and/or heating oil burning (Lin et al., 2018), on the HOA concentrations in the late evening. Moreover, eBC ff shows a clear weekend effect, with a less-pronounced pattern on Saturday and Sunday due to road transport reduction, while HOA displays intense nighttime peaks during the weekend. This HOA mean pattern was substantially influenced by winter data, whereas sum- mertime patterns display better consistency between HOA, eBC ff and NO x (Fig. S13). Altogether, these results suggest that this HOA is considered a mixed factor partly composed of both traffic and residential heating aerosols. This statement is in good agreement with conclusions from complementary studies showing wood burning contribution to HOA at the same site Srivastava et al., 2019).
As presented in Figs. 1a and S9, MO-OOA mass spectra present a strong peak at m/z 44. In fact, this spectrum has been widely reported as low volatility OOA (LV-OOA) and considered composed of highly oxidized and aged SOA Ulbrich et al., 2009;Zhang et al., 2011;Ng et al., 2011a). Overall, MO-OOA had a weak correlation (r 2 = 0.23) with sulfate during the entire period, supporting their different source origins to some extent. Compared to the poorly pronounced diel variability of sulfate, this MO-OOA factor exhibits a slight enhancement at nighttime (Fig. 2c), suggesting a possible local formation mechanism involving nighttime chemistry, on top of its overall regional feature.

Y. Zhang et al.: Six-year source apportionment of submicron organic aerosols
The geographic origins of the MO-OOA factor are further discussed in Sect. 4.2.1 for each season.
The mass spectra of LO-OOA in this study present a higher m/z 43 and a lower m/z 44 (Figs. 1a and S9) compared to MO-OOA, which is consistent with the mass spectral pattern of previously reported freshly formed semivolatile OOA (SV-OOA) Ng et al., 2010). The diel variations in LO-OOA display higher concentrations during nighttime than daytime (Fig. 2d), with relative variations much more pronounced than for the MO-OOA diel pattern, highlighting important roles of nighttime chemistry and/or gas-particle partitioning in the LO-OOA formation. These results support different formation pathways of the two OOA fractions. In addition, LO-OOA and nitrate present different diel cycles, suggestive of different formation processes and sources between each other. Different diel cycles of LO-OOA in different seasons have been also observed, which are further discussed in Sect. 4.2.1. In winter, LO-OOA mass spectra have higher contributions of m/z 29 and elevated m/z, i.e., starting from m/z 60, than during other seasons (Fig. S9). Such characteristics suggest a major influence of biomass burning emissions onto the LO-OOA factor during wintertime, as previously proposed from measurements in the Paris region (e.g., Crippa et al., 2013b). Conversely, in summer, this factor may be significantly influenced by BSOA formation Daellenbach et al., 2017). To investigate this possible origin, we checked if the summertime LO-OOA concentrations at higher daily temperature were following the temperature dependence similar to the one expected for the formation of terpene SOA, as explained in Sect. 3.2. Results of these calculations are presented in Fig. 3. The LO-OOA concentrations substantially increase with temperature, showing a good agreement with the estimated BSOA t formation exponential profiles. However, when comparing with estimation derived from Eq. (2) (referred to Fig. 3), the observed LO-OOA displays substantially higher loadings than estimated BSOA t at highest concentration range. This could be partly due to the influence of regional transport and atmospheric dilution on aerosol loadings and some possible uncertainties (such as unclear formation mechanism of biogenic SOA at SIRTA) which were not considered in the BSOA t estimation. These comparison results between observation and estimation indicates that the LO-OOA factor observed in summer might be mainly associated with biogenic sources. This is aligned with the VOCs seasonal patterns observed in the Paris region (Baudic et al., 2016), although the underlying SOA formation mechanism is still unclear and needs to be further investigated (Beekmann et al., 2015). Further discussion about seasonality of the LO-OOA factor is given in Sect. 4.2.1.  (Schurgers et al., 2009;Leaitch et al., 2011).  Figure 5 illustrates monthly average concentrations obtained for each OA factor over the studied period. HOA monthly concentrations vary from 0.4 to 1.3 µg m −3 and display a statistically insignificant seasonal trend (p > 0.05, Fig. 5a). Nevertheless, the mass concentration of HOA is nearly twice higher during cold months (in the range of 0.9-1.3 µg m −3 , from November to March) than in other months (in the range of 0.4-0.5 µg m −3 from April to October). This monthly cycle of HOA could be partially explained by lower BLH conditions and the influence of more intense emissions of biomass burning in cold seasons than in warm seasons (Figs. 4 and S14). As illustrated by Fig. S15, HOA clearly presents two peaks (in the morning and late evening) for each season. The evening HOA peak is about 1.3-1.5 times higher than the morning peak in winter and fall seasons when high loadings of BBOA are observed as well. Although dy- namic processes (establishment of a stable nighttime boundary layer) cannot be excluded, these results point to a possible contribution of biomass burning emissions to the HOA factor in the evening during cold months, as discussed before from the diel cycles of OA factors.

Monthly and seasonal variations in OA factors
As shown in Fig. 5b, BBOA displays a statistically significant seasonal pattern trend (p < 0.0001) with higher monthly mean concentrations (1.1-1.9 µg m −3 ) during cold months (November-March) than during the April-September period (0.3-0.5 µg m −3 ). This seasonal dependence of wood burning emissions is associated with the residential heating activities in the Paris region. BBOA presents a seasonal dependence of its diel cycle, as presented in Fig. S15. In particular, BBOA shows an evident peak at evening and nighttime in winter, spring and fall, while it presents a stable diel cycle during summertime. The highest seasonally averaged nighttime peak (up to 2.4 µg m −3 ) is observed in winter, highlighting a significant enhancement of wood burning emissions and influence of meteorological conditions (such as low BLH) during this season (Fig. S14).
Monthly average mass concentrations of MO-OOA present a significant seasonal trend (p < 0.05), varying from 1.0 in September to 3.5 µg m −3 in March (Fig. 5c), in agreement with previous studies performed in Europe (Schlag et al., 2016;Daellenbach et al., 2017;Bozzetti et al., 2017). The highest MO-OOA mass concentrations observed in the cold months are somehow similar to the seasonal variation in BBOA. MO-OOA diel cycles also present a seasonal variation, with significant increase during evening and nighttime in winter, spring and fall (Fig. S15). In order to minimize the effect of atmospheric dilution and regional transport, the mass concentration of MO-OOA was normalized to sulfate, the latter one being considered a regional secondary production marker (Petit et al., 2015 and Fig. S16). As shown in Fig. S16, the correlations between MO-OOA and sulfate are found to be strongly BBOA and wind speed dependent. For high wind speed and low BBOA concentrations, the mean MO-OOA-to-sulfate ratio is close to 1, while it reaches up to 8 under high BBOA and low-to-medium wind speed. This is consistent with the assumption of an enhancement of MO-OOA formation in the presence of substantial biomass burning emissions, which have been reported as a major anthropogenic SOA source (Heringa et al., 2011;Tiitta et al., 2016;Bertrand et al., 2017;Stavroulas et al., 2019;Daellenbach et al., 2019). Actually, both MO-OOA and LO-OOA factors may be significantly influenced by wood burning emissions as they are displaying similar correlations with eBC wb for highest MO-OOA-to-sulfate ratios during wintertime (Fig. 6).
As shown in Fig. S17, high concentrations of MO-OOA are generally observed at high RH (> 80 %) and low temperature (< 0 • C) conditions during wintertime. This low air temperature condition could be associated with a possible scenario for increase in the MO-OOA precursor emissions from biomass burning by residential heating activities during wintertime. As Fig. S18 shows, the MO-OOA-to-sulfate ratio shows a significant enhancement as a function of RH in winter, suggesting that the aqueous-phase heterogeneous processes may represent an important pathway for the local MO-OOA formation in winter as proposed by Gilardoni et al. (2016). Conversely, there are no obvious RHtemperature dependent patterns for the MO-OOA in spring (Fig. S17), indicative of more complex formation processes during this season. In summer, MO-OOA displays an evident increase from early afternoon to evening (Fig. S15), suggesting significant local photochemical production of SOA particles in summer with higher temperature and increased solar radiation . As a matter of fact, MO-OOA presents high concentrations under high temperature (> 25 • C) and low RH (< 65 %) summertime conditions (Fig. S17). In conclusion, and despite relatively constant mass spectra all over the year, MO-OOA appears to originate from various season-dependent formation pathways and sources (such as biomass burning and biogenic sources) that should still be investigated in more detail.
The LO-OOA mass spectra with high f 43 /f 44 ratios are frequently observed in spring, summer and fall, whereas a lower ratio is obtained for winter (Fig. S9). These different mass spectra of LO-OOA could be partially explained by season-dependent formation mechanisms and sources. The monthly mean mass concentrations of LO-OOA vary from 0.8 to 3.6 µg m −3 (Fig. 5d) and show a statistically significant seasonality (p < 0.001) with higher concentrations during warm months and lower during cold months. As discussed above, the high summertime LO-OOA concentrations are assessed to be mainly linked with BSOA formation. As presented in Fig. S17, the temperature-RH dependence of the LO-OOA factor is very different according to the season. In particular, the highest wintertime LO-OOA concentrations are mainly observed at low temperature and high RH conditions, suggesting that gas-particle partitioning may play an important role in LO-OOA formation during this season. As shown in Fig. 6b, LO-OOA shows a good correlation (r 2 = 0.55) with eBC wb at high MO-OOA-to-sulfate ratio conditions during wintertime, in agreement with the case of MO-OOA, which suggests that wood burning emissions could also contribute to the LO-OOA formation during winter. This is consistent with freshly formed SOA originating from primary biomass burning emissions, as reported by previous studies (Crippa et al., 2013b;Zhang et al., 2015;Sun et al., 2018;Stavroulas et al., 2019). In summer, the LO-OOA concentrations present a strong positive temperature dependence, while the RH dependence is not clear, indicating that photochemical production of LO-OOA became more important in summer than in winter. Moreover, high concentrations of the summertime LO-OOA are observed during the two distinct time periods in a day, i.e., early afternoon (around 12:00-15:00 LT) and nighttime (around 22:00-05:00 LT), which is different from the diel variations in other seasons with high concentrations only during nighttime (Fig. S15). These LO-OOA diel variations may reflect different formation pathways across a day in summer. Photochemical processes might dominate the LO-OOA production at daytime, while nighttime chemistry and/or gas-particle partitioning might promote its formation at low temperature conditions at night. Figure 7 presents the results obtained from the trend analysis of the 6-year monthly time series of the four OA factors, total OA, the two eBC components, secondary inorganic species (sulfate and nitrate) and total PM 1 (the sum of NR-PM 1 and eBC). The significance and magnitude of these trends were examined using the MK p value and Sen's slope, respectively. Annual mean mass concentrations of these OA factors and aerosol chemical species are given in Table 1.  BBOA presents a statistically significant decreasing trend (p < 0.05) with a Sen's slope of about −65 ng m −3 per year in the Paris region. On the other hand, eBC wb concentration trends appear quite stable over the investigated period. Two possible reasons may explain the discrepancy in trend results between BBOA and eBC wb . It may be hypothesized that a limited overall improvement of woodstove performances in the Paris region could have influenced BBOA emission factors more than eBC wb ones (Saleh et al., 2014) but no evidence has been found to support this assumption. Similarly, but in the opposite way, eBC ff was found to have a significant decreasing trend, while the HOA trend was found to be statistically insignificant (p value > 0.05). However, if removing the high concentration peaks observed in December (for which an important contribution of wood burning HOA can be expected), the MK p value is reduced to be 0.03, which would be indicative of a significant decreasing trend (with a related Sen's slope of −59 ng m −3 yr −1 ). These results would be in line with a reduction of PM 1 traffic emissions (−37 %) over the past years (2012)(2013)(2014)(2015)(2016)(2017) in France, as estimated by the French emission inventory state opera-tor (CITEPA, 2018). However, such trend analyses should be performed on longer datasets for a much better evaluation of the pollution control strategies (both on road transport and residential heating emissions) in the Paris region.

Long-term temporal trends
MO-OOA shows a significant decreasing trend (p < 0.05) with a Sen's slope of −175 ng m −3 per year. Considering the overwhelming secondary origin of this factor, this significant decreasing trend may be partially explained by an overall reduction of anthropogenic VOCs emissions (−13 %) in France (CITEPA, 2018) and even in a larger spatial scale, e.g., the western European region, during 2012-2017. LO-OOA presents no significant trend (with p > 0.05). As discussed above, higher LO-OOA loadings may be linked to BSOA formation, especially at summertime. The stability of LO-OOA concentrations over time may be linked to limited changes in biogenic VOC emissions and/or in relevant oxidant concentrations that control the SOA burden in the atmosphere. The effect of anthropogenic-biogenic interaction mechanisms on biogenic SOA formation -e.g., involving NO x , as reported by previous studies in urban regions (Budisulistiorini et al., 2015;Zhang et al., 2017) -could Figure 7. Temporal trends of monthly mass concentrations of different chemical speciation, including eBC wb , eBC ff , SO 4 (sulfate), NO 3 (nitrate), four OA factors, total OA and total PM 1 (the sum of NR-PM 1 and eBC). The MK or seasonal MK tests associated with estimated Sen's slope (µg m −3 yr −1 ) were used for the trend analysis. also partially explain the limited changes for the long-term trend of LO-OOA at SIRTA. Detailed LO-OOA formation processes involved here still need to be further investigated. Nevertheless, it may be assumed that reductions of anthropogenic VOC emissions only cannot be sufficient to weaken the total SOA background concentrations in the Paris area.
The trends are more significant for total OA (p < 0.005, about −382 ng m −3 yr −1 ), as well as for nitrate (p of about 0.01 and approximately −145 ng m −3 yr −1 ) and total PM 1 (p of about 0.002, and approximately −644 ng m −3 yr −1 ), than for sulfate (with a p value around 0.5117). In addition to SO 2 emitted from anthropogenic sources (e.g., industrial and shipping emissions) (Hoesly et al., 2018), natural sources (e.g., volcanic emissions) (Boichu et al., 2019) could also in-fluence the sulfate budget in western Europe. This suggests that regional aerosol chemistry modeling simulations by using different SO 2 emission sectors may help to further explain the temporal trends of sulfate. Meanwhile, it should be noted that the sulfate trend here could be probably influenced by the ACSM measurement uncertainties Freney et al., 2019), which could not be fully quantified here. Overall, these decreasing trends could reflect the response of the PM concentrations to the decrease in anthropogenic source emissions during these last years in Europe. Reduction in NO x (−19 %) came with a negligible change in NH 3 emissions (+2 %) over the French region during recent years (2012)(2013)(2014)(2015)(2016)(2017), which may support the fact that the decreasing trend in particulate nitrate was likely driven by the NO x emissions control in the Paris region (CITEPA, 2018). A continuous effort to reduce POA emissions and SOA precursors (VOCs) may lead to the decrease in both the total OA and fine PM budgets (EMEP, 2016;CITEPA, 2018). Figure 8 presents the contribution of the four OA sources as a function of total submicron OA mass loadings or each season along with percent changes in meteorological conditions. In winter and fall ( Fig. 8e and h), all meteorological parameters -except limited changes in RH -show negative relationships as a function of the OA mass concentrations, confirming the coincidence of low temperature, low WS and/or low BLH in the formation of pollution episodes (Dupont et al., 2016). POA contributions gradually increase from 35 % (27 %) up to 64 % (70 %) as a function of OA mass concentrations in winter (fall) ( Fig. 8a and d). These results illustrate the major role of primary sources during periods with high OA concentrations during the cold seasons. In particular, BBOA contribution gradually increase from 21 % (15 %) to 41 % (40 %) in winter (fall) along with OA mass loading increase. In spring, OA composition is radically changed and is dominated by the two OOA fractions, with almost constant average contributions (68 %-77 %), regardless of OA concentration levels (Fig. 8b), indicating the major role in SOA during this season. MO-OOA presents higher contributions to OA (45 %-53 %) than LO-OOA (15 %-31 %), suggesting that the formation of aged SOA plays a key role on the buildup of episodes with high OA concentrations during springtime. As shown in Fig. 8f, the percent changes in temperature, WS and BLH gradually decrease with increasing OA concentrations. By contrast, RH shows a positive relationship with OA mass concentrations, with the largest RH enhancement (16 %) at the highest OA-loading bin (> 25 µg m −3 ). This may suggest that high RH is the most favorable environmental condition for SOA formation during springtime OA pollution episodes, as supported by a high contribution of OOA factors at the highest OA concentration level (Fig. 8b). In addition, although BBOA contributions remained relatively limited, they increase from 11 % to 17 % when OA increased from less than 10 to > 25 µg m −3 . This may reveal a non-negligible influence of wood burning emissions during early spring pollution episodes.

OA source contribution as a function of OA concentrations
In summer, OA was also dominated by the two OOA fractions (around 80 %-85 % at all OA-loading bins). The LO-OOA contribution gradually increased from 51 % to 69 % as a function of OA mass loading associated with a significant increase in temperature ( Fig. 8c and g). Other meteorological variables (i.e., RH, WS and BLH) showed relatively stable changes across different OA mass loadings (Fig. 8g). These results confirm that high OA concentrations during summer are strongly determined by temperature-driven biogenic SOA formation processes. Figure 9 presents seasonal wind-dependent variations of the POA factors. Wind rose plots (Fig. 9a) present the greatest frequency (up to approximate 32 %) of winds associated with southwestern wind direction in all seasons. In spring, the occurrences of air masses from the west, northwest, north and northeast sectors are comparable (about 15 %). As presented in Fig. 9b, BBOA shows high concentrations associated with the eastern wind sectors, except in summer, which is in agreement with the location of the residential area on the eastern and western parts of SIRTA. Compared to other seasons, high concentrations of BBOA are also observed linking to western wind sectors in winter, which may imply more intense biomass burning from larger scales during colder months. As discussed above, HOA is a mixed factor with biomass burning aerosols during wintertime, which therefore presents a similar wind-dependent pattern as BBOA ( Fig. 9b and c). In spring, summer and fall, HOA presents a distinct pattern with high concentrations associated with northeastern wind sectors from urban area of Paris, suggesting that the short-range transport from the urban Paris area may strongly impact the HOA concentrations at SIRTA. Figure 10 shows maps of the most probable geographic origins of the two OOA factors for each season based on PSCF analysis. In winter, MO-OOA presents high PSCF values over the Benelux, Germany and Poland, showing a major influence of long-range transport of OA from northeastern sectors. Those results may suggest more intense SOA production and aging processes at regional scale for continental air masses. As a matter of fact, MO-OOA shows wider potential source regions than LO-OOA, which is assessed as fresh SOA and could be mainly formed at more local scale in winter. Moreover, the impact of transport from northeastern regions -hosting intense anthropogenic activities (e.g., industries) -onto MO-OOA concentrations may also support a significant anthropogenic origin for this aged SOA factor.

Potential geographic origins
As shown in Fig. 10c and d, both MO-OOA and LO-OOA present high springtime PSCF values originating from the northeastern regions too, which can participate in pollution episodes frequently observed during this season Srivastava et al., 2018b). Therefore, mitigation of VOC emissions at the regional scale could help to reduce the substantial influence of OA on PM limit value exceedances during this season. A narrower distribution of potential source regions was observed in summer and fall, compared to winter and spring. MO-OOA presents potential source regions mainly from the northeast in summer, while it has a high potential source region originating from the south in fall (Fig. 10g).
All these results indicated that significant reduction of the SOA burden in the Paris region does not only require the limitation of local source emissions but also needs a synergistic control strategy for the regional sources, especially from northeastern European regions. In this respect, they confirmed conclusions reached by previous short-term campaigns (e.g., Sciare et al., 2010;Crippa et al., 2013b;Freutel et al., 2013;Beekman et al., 2015).

Conclusions
A comprehensive OA source apportionment has been achieved over the region of Paris from November 2011 to March 2018. Four factors, comprising HOA, BBOA, MO-OOA and LO-OOA, have been identified and selected to ensure consistency of PMF factor solution over 6 years in this study. Mean annual contributions of these factors to OA were of 11 %-16 % (HOA), 14 %-19 % (BBOA), 25 %-42 % (LO-OOA), and 30 %-45 % (MO-OOA), respectively. BBOA presented a statistically significant seasonal pattern with highest concentrations during cold months, due to residential wood burning emissions. The contribution of BBOA increased with increasing concentration of OA mass in winter and fall -along with decreasing boundary layer height and wind speed -highlighting the importance of biomass burning emissions for OA pollution under stagnant meteorological conditions. HOA presented temporal variations similar to BBOA in cold months, which was partly related to the fact that wood burning emissions also contributed to HOA burden. BBOA and HOA exhibited very limited (< 70 ng m −3 yr −1 ) and not significant (at the 5 % significance level) decreasing trends. These results imply that specific mitigation strategies (e.g., emission control), especially for residential wood burning, are still necessary for substantial improvement of air quality in the cold season in the Paris region. Moreover, residential biomass burning emissions could be assumed as an important air pollution source over western Europe, while such a source remains largely unregulated.
LO-OOA and MO-OOA presented different seasonal variations, reflecting different formation mechanisms and/or precursor sources. LO-OOA displayed a pronounced seasonal cycle, with highest contributions to total OA in summer  . Maps for potential source origins of regional transport that may contribute to SOA (including MO-OOA and LO-OOA) burdens at SIRTA. Observed data points with wind speed (less than 4 m s −1 ) and in the presence of precipitation events are filtered for the PSCF calculation. The solid black point in each plot represents the location of the sampling site.
(50 %-66 %) and lowest ones in winter (12 %-19 %). Enhanced LO-OOA production during the warm season was assessed to be mainly driven by biogenic SOA formation. This factor showed no significant long-term trend (p > 0.05) for the studied period. MO-OOA presented a higher contribution to OA in wintertime (35 %-51 %) and springtime (32 %-62 %) than during the rest of the year. PSCF analyses suggested a high probability of MO-OOA long-range transport from northeastern Europe towards the Paris region. MO-OOA displayed a significant (p < 0.05) decreasing trend (of about −175 ng m −3 yr −1 ), which might reflect the effect of the emission control strategy of anthropogenic SOA precursors at the regional scale over the last decade. However, future work is needed to fully understand chemical properties of these SOA factors corresponding to different origins over different seasons in the Paris region and to quantify the impact of emission control on ambient SOA burden.
Data availability. The data have been presented in the text, table and figures as well as the Supplement. Additional related data can be made available upon request (olivier.favez@ineris.fr).
Author contributions. OF, AA and VG designed and led the study. YZ conducted the data analyses. J-EP, FT, NB, VC, TA and JS provided the field observations. FC and AP supported the source apportionment analyses. YZ and OF interpreted the data and wrote the article, with inputs from all coauthors.
Competing interests. The authors declare that they have no conflict of interest.