Source apportionment of atmospheric PM10 oxidative potential: synthesis of 15 year-round urban datasets in France

Abstract

Reactive oxygen species (ROS) carried or induced by particulate matter (PM) are suspected of inducing oxidative stress in vivo, leading to adverse health impacts such as respiratory or cardiovascular diseases. The oxidative potential (OP) of PM, displaying the ability of PM to oxidize the lung environment, is gaining strong interest in examining health risks associated with PM exposure. In this study, OP was measured by two different acellular assays (dithiothreitol, DTT, and ascorbic acid, AA) on ${\mathrm{PM}}_{\mathrm{10}}$ filter samples from 15 yearly time series of filters collected at 14 different locations in France between 2013 and 2018, including urban, traffic and Alpine valley site typologies. A detailed chemical speciation was also performed on the same samples, allowing the source apportionment of PM using positive matrix factorization (PMF) for each series, for a total number of more than 1700 samples. This study then provides a large-scale synthesis of the source apportionment of OP using coupled PMF and multiple linear regression (MLR) models. The primary road traffic, biomass burning, dust, MSA-rich, and primary biogenic sources had distinct positive redox activity towards the OP${}^{\text{DTT}}$ assay, whereas biomass burning and road traffic sources only display significant activity for the OP${}^{\text{AA}}$ assay. The daily median source contribution to the total OP${}^{\text{DTT}}$ highlighted the dominant influence of the primary road traffic source. Both the biomass burning and the road traffic sources contributed evenly to the observed OP${}^{\text{AA}}$. Therefore, it appears clear that residential wood burning and road traffic are the two main target sources to be prioritized in order to decrease significantly the OP in western Europe and, if the OP is a good proxy of human health impact, to lower the health risks from PM exposure.

1 Introduction

Air quality has become a major public health issue, being considered the fourth largest global cause of mortality, with 7 million premature deaths worldwide per year due to both indoor and outdoor exposure . Possibly driving 90 $\mathrm{\%}$ of this health impact , particulate matter (PM) is one of the key pollutants in the air linked to health outcomes, although the exact mechanism leading to toxicity is not yet fully understood . Many urbanized areas, mainly located in low- or middle-income countries, are exposed to PM concentrations far higher than the recommendation guideline of the WHO.

Although PM is now monitored in most countries and large efforts are observed to document ambient concentrations, the underlying processes leading to the observed concentrations in the atmosphere, and particularly the understanding of emission sources, are still active fields of research . In recent years, strong focus has been put worldwide on source-apportionment methods in order to better understand the processes leading to the airborne concentrations and the accumulation of PM in the atmosphere. This includes direct modeling approaches such as chemistry transport models (CTMs) using tagged species or field studies coupled with receptor models (RMs) , notably positive matrix factorization (PMF). PMF can be based either on an aerosol mass spectrometer (AMS) time-resolved spectrum or on filter analysis or a mix of these different measurement techniques . Results of these studies indicate that PM originates from a wide variety of sources, not only natural (volcano, sea spray, soil dust, vegetation, bacteria, pollen, etc.) or anthropogenic (road traffic, residential heating, industry, etc.) sources, but it is also formed as a secondary product and condensed from the gaseous phase (ammonium nitrate and ammonium sulfate, etc.). As a result, the chemistry, size distribution or reactivity of PM widely vary from location to location and season to season, which induces large changes in the health impacts depending on all of these parameters .

Furthermore, the mass of PM is not the most relevant metric when dealing with health impacts of airborne particles since major properties (chemistry, shape, size distribution, solubility, speciation) driving PM toxicity are not taken into account within this single mass metric. It is now believed that the measurement of the reactive oxygen species (ROS) issued from PM may be more closely linked to the potential adverse health effects of atmospheric PM, since oxidative stress is a key factor in the inflammatory response of the organism, leading for instance to respiratory diseases or when exposed for a long period of time, cardiovascular diseases or even cancer . Therefore, the oxidative potential (OP) of PM being an indirect measure of the ability of the particles to induce ROS in a biological medium has been proposed as a potential proxy of the health impacts of atmospheric PM exposure. Indeed, even if the clear demonstration of the OP as a good proxy of health impact is still needed, some recent studies already established associations between OPs and different possible health outcomes or cellular stress in vitro .

However, there is also no consensus on a standardized method to measure the OP of PM, and many assays and protocols co-exist (DTT, GSH, AA, ESR, ${}^{\circ }$OH or ${\mathrm{H}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$, among others), with samples extracted with different methods (water, methanol, simulated lung fluid (SLF), etc.) and not always with a constant mass of PM. The dithiothreitol (DTT) and ascorbic-acid (AA) assays are widely used in association with health endpoints , even if the exact methodologies differ from one study to the other. Results can also differ for the seasonality of OP based on these two assays, and some studies report strong seasonality of OP, whereas others do not . Finally, several studies have already shown that different sources of PM have different reactivities to OP tests . In particular, sources with high concentrations of transition metals, such as road traffic, appear to have a higher intrinsic oxidative potential (i.e., OP per microgram of PM) than other sources of PM. However, the number of these studies is still limited, and they do not always take into account complete seasonal cycles and therefore may not encompass the variety of sources for a given site, possibly omitting some important sources. Also, spatial variability at a country scale is currently unknown and requires homogeneous sampling and analysis methodologies for all filters and time series.

Figure 1

Locations of the 14 sampling sites. Color codes denote the typologies of the sites: red, urban; orange, urban valley; magenta, industrial; blue, traffic. Background photography from NASA's Earth Observatory.

[Figure omitted. See PDF]

For a comprehensive investigation of the intrinsic OPs of various PM sources, we built up an extensive dataset of about 1700 samples from 14 sites consisting of 15 year-round time series of observations over continental France, collected during research programs conducted between 2013 and 2018. In each of these samples, we concurrently measured the OP with the DTT and AA assays together with an extensive chemical characterization allowing PM source apportionment using a harmonized PMF approach . Then, we apportioned the OP measured by the DTT and AA assays to the emission sources using a multilinear regression approach following . In this way, we can estimate the oxidizing capacity of each microgram of PM from the different identified emission sources but also the relative contribution of the different sources to the OP${}^{\text{DTT}}$ and OP${}^{\text{AA}}$ on seasonal and daily bases. These results are presented in this paper.

2 Materials and methods

2.1 Site descriptions

The selected sites had to fulfill three conditions: (1) a yearly sampling period, (2) the required chemical analysis to perform a harmonized PMF analysis and (3) enough filter surface left to assess the OP measurements. A total of 14 sites were included in this study (one being sampled twice at a 5-year interval) taken from different research programs. These sites reflect the diversity of typology we could encounter in western Europe: urban (NGT, TAL, AIX, MRS-5av, NIC), urban Alpine valley (GRE-cb, GRE-fr, VIF, CHAM, MNZ, PAS), industrial (PdB) and traffic (RBX, STG-cle) (see Table S1 in the Supplement) covering different areas of France (Fig. ). We can note, however, the absence of remote or rural sites in our current dataset. The air quality at all of these sites is operated by the local air quality agencies (Atmo Sud, Atmo Auvergne Rhône-Alpes, Atmo Nouvelle Aquitaine, and Atmo Hauts de France), and most of them are part of the so-called CARA program .

2.2 Sample analysis

Even if it has been shown that mainly ${\mathrm{PM}}_{\mathrm{2.5}}$ is deposited in lung alveoli , ${\mathrm{PM}}_{\mathrm{10}}$ is still a public health concern and under regulation in the EU and France (Directive ). Moreover, recent studies also highlight the role of the coarse fraction of PM in health impacts . ${\mathrm{PM}}_{\mathrm{10}}$ has the advantage of encompassing all parts of PM potentially reaching the lower respiratory tract.

2.2.1 Chemical speciation

The ${\mathrm{PM}}_{\mathrm{10}}$ concentrations were measured at each site by means of an automatic analyzer, according to EN 16450:2017 , and daily (24 $\mathrm{h}$) filter samples were collected every third day. Samplings were achieved on pre-heated quartz fiber filters using a high-volume sampler (DA80, Digitel), following EN 12341:2014 procedures . Offline chemical analysis performed on these filters has been fully described previously . Briefly, the elemental and organic carbon fractions (EC and OC) were measured via thermo-optical analysis (Sunset Lab. Analyzer ) using the EUSAAR-2 protocol . Major water-soluble inorganic contents (Cl${}^{-}$, ${\mathrm{NO}}_{\mathrm{3}}^{-}$, ${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$, ${\mathrm{NH}}_{\mathrm{4}}^{+}$, Na${}^{+}$, K${}^{+}$, Mg${}^{\mathrm{2}+}$, and Ca${}^{\mathrm{2}+}$) and methanesulfonic acid (MSA) were determined using ion chromatography . Metals or trace elements (e.g., Al, Ca, Fe, K, As, Ba, Cd, Co, Cu, La, Mn, Mo, Ni, Pb, Rb, Sb, Sr, V, and Zn) were measured by inductively coupled plasma atomic emission spectroscopy or mass spectrometry (ICP-AES or ICP-MS) . Finally, anhydrosugars and sugar alcohols (including levoglucosan, mannosan, arabitol, sorbitol, and mannitol) were analyzed using liquid chromatography followed by pulsed amperometric detection (LC-PAD) .

2.2.2 OP assays

Identical methodologies were applied for all the OP measurements of the collected filters . The extraction of PM was performed using a simulated lung fluid (SLF: Gamble $+$ DPPC) to simulate the bio-accessibility of PM and to closely simulate exposure conditions. In order to take into account the nonlinearity of the OP with PM mass and to have comparable results between sites, the extraction has been carried out at the iso-mass concentration (10 or 25 $\mathrm{\mu }\mathrm{g}{\mathrm{mL}}^{-\mathrm{1}}$ of PM, depending on the site, both values being in the low mass range of linear dose response), by adjusting the surface of the filter extracted. The filter extraction method allowed us to include both soluble and insoluble particles in the extracts. After the SLF extraction, particles removed from the filter were not filtrated, and the whole extract was injected in a 96-well plate for analysis. Samples were processed using the AA and DTT assays, as described below.

DTT depletion when in contact with PM extracts was determined by dosing the remaining amount of DTT with DTNB (dithionitrobenzoic acid) at different reaction times (0, 15 and 30 $\min$), and absorbency was measured at 412 $\mathrm{nm}$ using a plate spectrophotometer (Tecan, M200 Infinite). An initial concentration of 12.5 $\mathrm{nmol}$ of DTT is used (50 $\mathrm{\mu }\mathrm{L}$ of 0.25 mM DTT solution in phosphate buffer) to react with 205 $\mathrm{\mu }\mathrm{L}$ of phosphate buffer and 40 $\mathrm{\mu }\mathrm{L}$ of PM suspension.

The AA assay is a simplified version of the synthetic respiratory tract lining fluid (RTFL) assay , where only AA is used. A mix of 80 $\mathrm{\mu }\mathrm{L}$ of PM suspension with 24 $\mathrm{nmol}$ of AA (100 $\mathrm{\mu }\mathrm{L}$ of 0.24 mM AA solution in Milli-Q water) is used, and AA depletion was read continuously for 30 $\min$ by absorbency at 265 $\mathrm{nm}$ (TECAN, M1000 Infinite). The depletion rate of AA was determined by linear regression of the linear section data. For both assays, the 96-well plate was auto-shaken for 3 $\mathrm{s}$ before each measurement and kept at 37 ${}^{\circ }\mathrm{C}$.

Three filter blanks (laboratory blanks) and three positive controls (1,4-naphthoquinone, 24.7 $\mathrm{\mu }\mathrm{mol}{\mathrm{L}}^{-\mathrm{1}}$) were included in each plate (OP${}^{\text{AA}}$ and OP${}^{\text{DTT}}$) of the protocol. The average values of these blanks were then subtracted from the sample measurements of the given plate. The detection limit (DL) value was defined as 3 times the standard deviation of laboratory blank measurements (laboratory blank filters in Gamble $+$ DPPC solution).

Hereafter, the OP${}^{\text{DTT}}$ and OP${}^{\text{AA}}$ normalized by air volume are denoted ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, respectively, with unit $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$.

2.3 Source apportionment

The source apportionment of the OP can be performed in two main ways: (1) by including the OP as an input variable for the RM or (2) by conducting source attribution to the PM mass and then, using a multilinear regression (MLR) model, assigning OP to each of the sources from the source–receptor model . We decided to use the second approach since adding the OP variable to the PMF may change the source-apportionment solution. Further, the first method would imply that the intrinsic OP would be positive by construction and not due to relevant physical properties (see below). Moreover, the OP apportionment in a two-step process (PM mass source apportionment and then OP apportionment) allowed us to potentially use different model types for the OP apportionment and re-use results of the PMF for other OP assays later on.

2.3.1 PM mass apportionment: positive matrix factorization

Methodological background

The PM mass source apportionment for the 15 yearly series was conducted using the U.S. Environmental Protection Agency (US-EPA) EPA PMF 5.0 software with the ME-2 solver from . Briefly, the PMF was introduced by and is now one of the most common approaches used for PM source-apportionment studies . It aims at solving the receptor model equation

1$${\displaystyle }{\displaystyle }\mathbf{X}=\mathbf{G}\cdot \mathbf{F},$$ where $\mathbf{X}$ is the $n\times m$ observation matrix, $\mathbf{G}$ is the $n\times p$ contribution matrix and $\mathbf{F}$ is the $p\times m$ factor profile matrix (or source, despite some factors not being proper emission sources but possibly reflecting secondary processes), with $n$ the number of samples, $m$ the number of measured chemical species and $p$ the number of profiles. Hereafter, the $\mathbf{G}$ matrix is expressed in $\mathrm{\mu }\mathrm{g}{\mathrm{m}}^{-\mathrm{3}}$ and the $\mathbf{F}$ matrix in $\mathrm{\mu }\mathrm{g}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ of PM.

PMF setup

Some of the PMF analyses included in this paper have been performed during previous programs, namely, SOURCES (http://pmsources.u-ga.fr, last access: 27 July 2021, ), DECOMBIO , or MobilAir (https://mobilair.univ-grenoble-alpes.fr/, last access: 27 July 2021, ). In order to get comparable PM source profiles from a common set of input species and constraints in the model, PMF analyses have been run with a harmonized methodology, as previously reported . Then, the PMF of MNZ and PAS were run specifically for this study, and the ones of MobilAir (GRE-fr_2017, GRE-cb and VIF) are the “classic” runs in . The other ones (AIX, CHAM, GRE-fr_2013, MRS-5av, NGT, NIC, PdB, RBX, STG-cle and TAL) directly come from the SOURCES program .

The input species were slightly different from one study to another but always included carbonaceous compounds (OC and EC), ions (${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$, ${\mathrm{NO}}_{\mathrm{3}}^{-}$, Cl${}^{-}$, ${\mathrm{NH}}_{\mathrm{4}}^{+}$, K${}^{+}$, Mg${}^{\mathrm{2}+}$, and Ca${}^{\mathrm{2}+}$), organic compounds (levoglucosan, mannosan, arabitol and mannitol – the latter two summed and referred to as polyols) and MSA – and a set of trace metals for a total of about 30 species. The list of metals used within the PMF analysis was not the same for each of the sites due to overly low concentrations (lower than the quantification limit) on some filters leading to a very low signal-to-noise ratio (see Table S2). The uncertainties were estimated following the method proposed by and were tripled if the signal-to-noise ratio was below 2 (classified as “weak” in the PMF software). Between 8 and 10 factors were identified at the different sites and are summarized in Table S3. For each of the PMF analyses, the possibility of using constraints to refine some of the chemical profiles of factors was considered in order to better disentangle possible mixing between factors and reduce the rotational ambiguity, based on a priori expert knowledge of the geochemistry of the sources. A PMF solution was considered valid if it followed the recommendation of the “European guide on air pollution source apportionment with receptor models” , also requiring a proper geochemical identification of the various factors. Estimation of the uncertainties of the PMF was obtained on both the base and constrained runs using the bootstrap (BS) and displacement (DISP) functions of the EPA PMF5.0 .

2.3.2 Similarity assessment of the PMF factors

Since PMF resolved site-specific PM factors, we checked whether a given factor had a consistent chemical profile over the different sites. For this purpose, as presented in , a similarity assessment of all PMF factor profiles was run following the DeltaTool approach . Using this tool, we compared pairs of factor profiles based on their mass-normalized chemical compounds using two different metrics, namely, the Pearson distance (PD) and the standardized identity distance (SID) . The first one is defined as

2$${\displaystyle }{\displaystyle }\text{PD}=\mathrm{1}-r^{\mathrm{2}},$$ with $r^{\mathrm{2}}$ the Pearson correlation coefficient, and might be strongly influenced by individual extreme points. The second one, SID, can expressed as follows: 3$${\displaystyle }{\displaystyle }\text{SID}={\displaystyle \frac{\sqrt{\mathrm{2}}}{m}}\sum _{j=\mathrm{1}}^m{\displaystyle \frac{\left|x_j-y_j\right|}{x_j+y_j}},$$ where $x$ and $y$ are two different factor profiles expressed in relative mass, and $m$, the number of common species in $x$ and $y$, is evenly sensitive to every species since it includes a normalization term.

2.4 OP apportionment

The computation was done thanks to the statsmodels 0.12 python package , and the graphics were produced with matplotlib 3.3.1 and seaborn 0.11 .

2.4.1 Apportionment using multilinear regression (MLR)

MLR was conducted independently at each site, with results from the two (DTT and AA) OP assays being the dependent variables and the source contributions obtained from the PMF being the explanatory variables, following the equation (Eq. ), similarly to :

4$${\displaystyle }{\displaystyle }{\text{OP}}_{\text{obs}}=\mathbf{G}\times \mathit{\beta }+\mathit{\epsilon },$$ where OP${}_{\text{obs}}$ is a vector of size $n\times \mathrm{1}$ of the observed ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ or ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ in $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$, $\mathbf{G}$ is the matrix ($n\times (p+\mathrm{1})$) of the mass contribution of PM sources obtained from the PMF in $\mathrm{\mu }\mathrm{g}{\mathrm{m}}^{-\mathrm{3}}$ and a constant unit term for the intercept (no unit), $\mathit{\beta }$ are the coefficients (i.e., intrinsic OP of the source and the intercept) of size ($(p+\mathrm{1})\times \mathrm{1}$) in $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ for the intrinsic OP and in $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$ for the intercept. The residual term $\mathit{\epsilon }$ ($n\times \mathrm{1}$) accounts for the misfit between the observations and the model.

The intercept was not forced to zero on purpose. Indeed, if the system is well constrained, the intercept should spontaneously be close to zero, and conversely a non-zero intercept would point out missing explanatory variables.

A weighted least-square regression (WLS) was finally used to consider the uncertainties of the OP measurements. The uncertainties of the coefficients $\mathit{\beta }$ given by the MLR were estimated by bootstrapping the solutions 500 times, randomly selecting 70 $\mathrm{\%}$ of the samples each time to account for possible remaining extreme events or seasonal variations of the intrinsic OP per source. The uncertainty of the PMF result $\mathbf{G}$ was however not considered because the EPA PMF software only returns to the user the uncertainties associated with the profile matrix $\mathbf{F}$ (see , for a first-order estimation of the $\mathbf{G}$ uncertainties).

2.4.2 Contribution of the sources to the OP

The contribution ${\mathbf{G}}^{\text{OP}}$ in $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$ of the sources to the OP was computed at each site independently and was calculated following Eq. ():

5$${\displaystyle }{\displaystyle }{\mathbf{G}}_k^{\text{OP}}={\mathbf{G}}_k\times {\mathit{\beta }}_k,$$ where $k$ is the source considered, $\mathbf{G}$ the PMF sources' contribution in mass concentration ($\mathrm{\mu }\mathrm{g}{\mathrm{m}}^{-\mathrm{3}}$) and $\mathit{\beta }$ the intrinsic OP of the sources in $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$. The uncertainties of ${\mathbf{G}}^{\text{OP}}$ were computed using the uncertainties of $\mathit{\beta }$ estimated from the 500 bootstraps.

2.5 Focus on the main PMF factors

This study focuses on the main drivers of OP at the regional scale. For this reason, we decided to include in the main discussion only the PMF factors identified in at least two-thirds of the series (i.e., 10 out of 15 series), namely, the aged salt, biomass burning, dust, MSA-rich, nitrate-rich, primary biogenic, primary road traffic and sulfate-rich factors. However, the remaining sources, often local, barely contributed to the total PM mass, and important uncertainties were often attached to them. The only notable exception is the HFO (heavy fuel oil) profile identified at some coastal sites, discussed hereafter in its own section.

3 Results and discussion

As a large set of results has been obtained in the present study and cannot be exhaustively presented here, an interactive visualization tool providing details on PM and OP source time series and apportionment outputs is available online at http://getopstandop.u-ga.fr/ (last access: 27 July 2021, and also http://pmsources.u-ga.fr, last access: 27 July 2021, for the chemical part of the SOURCES ones) and is proposed as the Supplement of this paper including all factors for all series.

Note that, since there are more samples where PMF has been run compared to available OP measurements (around 1700 concomitant OP and PMF samples compared to 2048 samples with a PMF solution), the discussion hereafter on the source contributions to the OP takes into account the whole PMF analysis, including days when models were not trained but predicted by the above-mentioned approach.

3.1 PMF results

In this section, we summarize the main results acquired from the harmonized PMF approach conducted for the present paper, but we invite the reader to refer to the previous study and to the website to have a more complete view of the results.

3.1.1 PMF source apportionment

The list of the identified factors at each site is given in Table S3, and individual profiles and time series together with uncertainties can be found at http://getopstandop.u-ga.fr/results?component=pmf_profile_and_contribution (last access: 27 July 2021). Table  summarizes the main PMF factors found at least at 10 out of the 15 series. Shortly, we obtained PMF factors corresponding to biomass burning (mainly from residential heating), primary road traffic, mineral dust, secondary inorganic nitrate-rich and sulfate-rich, salt (fresh and aged) as well as primary biogenic and MSA-rich. Some other local sources were also identified at some sites, targeting some local heavy-loaded metal sources with a very low contribution to the total PM mass – supposedly linked to the industrial process – which contained a wide variety of chemicals but shared a common set of metals (Al, As, Cd, Mn, Mo, Pb, Rb, Zn). Finally, a factor related to shipping emission (namely, HFO) was obtained at some coastal sites.

Table 1

Main PMF factors identified (at least at two-thirds of the series, i.e., 10 out of 15) and species used as proxies for the determination.

Factor name	Number of	Main species used as tracers	General remarks
	sites identified
Biomass burning	15	Levoglucosan, OC, EC, K${}^{+}$, Rb	High in Alpine valley, strong seasonality
Nitrate-rich	15	${\mathrm{NO}}_{\mathrm{3}}^{-}$, ${\mathrm{NH}}_{\mathrm{4}}^{+}$	Mostly in spring
Primary biogenic	15	Polyols (arabitol, mannitol), OC	Strong seasonality
Road traffic	14	OC, EC, Cu, Fe, Sb, Sn	Mixed exhaust and non-exhaust emission
MSA-rich	14	MSA	Strong seasonality
Dust	13	Al, Ti, Ca${}^{\mathrm{2}+}$	Episodic, some OC is present
Sulfate-rich	13	${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$, ${\mathrm{NH}}_{\mathrm{4}}^{+}$, Se	Some OC is present
Aged salt	12	Na${}^{+}$, Mg${}^{\mathrm{2}+}$, ${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$, ${\mathrm{NO}}_{\mathrm{3}}^{-}$	Some OC is present

3.1.2 PMF similarities between sites

The similarities between chemical profile composition estimated by the PD and SID metrics are presented in Fig. S2. We observed a strong similarity for the main sources of PM, namely, biomass burning, nitrate-rich, primary biogenic, sulfate-rich and to a lesser extent road traffic. The dust, aged salt and MSA-rich were often identified and presented acceptable SID but also showed large values for the PD metric. As the PD is sensitive to “extreme points”, this translates in our case to different contributions for the chemical species contributing most to the PM mass (mainly OC and EC). MSA-rich is the most variable factor, and a detailed analysis of its chemistry profile indicated many differences from site to site for the concentrations of EC but also ${\mathrm{NO}}_{\mathrm{3}}^{-}$ and ${\mathrm{NH}}_{\mathrm{4}}^{+}$. This factor being essentially a secondary organic factor, this variability may be explained by different formation or evolution pathways or different levels of aging. We can also point out that the industrial source had a very diverse chemical composition since it is related to different local industrial processes.

Nevertheless, the geochemical stability of the majority of PMF factors on a regional scale allows us to consider that these emission sources are rather homogeneous over France.

Figure 2

Box plot of ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ seasonal values. The numbers in the $x$ axis indicate the number of observations. Each box represents one decile, and the black horizontal lines indicate the median of the distributions. Some values greater than $\mathrm{15.5}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$ are not displayed for graphical purposes.

[Figure omitted. See PDF]

3.2 OP results and seasonality

The 15 time series for both ${\text{OP}}^{\text{DTT}}$ and ${\text{OP}}^{\text{AA}}$ at each site are presented on the website (http://getopstandop.u-ga.fr/results?component=rd_ts, last access: 27 July 2021). A monthly aggregated view is given in Fig.  for ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, respectively. As the dataset covered complete full years, including the influence of different PM sources with different seasonal activities, the results obtained are representative of spatio-temporal patterns of the OPs at least over France and probably over large parts of western Europe.

As reported previously by , we observed a seasonality of both ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, with higher OP values during the colder months (October–March) compared to the warmer months (April–September). We also noted that during the winter period, the statistical distribution of OP values did not follow a normal distribution, and a significant variability was observed. This was especially the case for the sites located in the Alpine area (GRE-fr, GRE-cb, VIF, CHAM, MNZ, PAS) showing stronger seasonality compared to the other locations. Such specificity was already reported previously by together with some rapid variation of the OP${}^{\text{DTT}}$ and OP${}^{\text{AA}}$, with drastic increase or decrease within the time frame of a few days, similarly to the ${\text{PM}}_{\mathrm{10}}$ mass concentration. This behavior may be related to the formation of thermal inversion layers in such valleys, leading to the accumulation of pollutants and the promotion of the secondary processes inducing increased formation of secondary organic aerosol (SOA) and of key organic species, like polyaromatic quinones or humic-like substances (HULIS) with a significant impact on the OP.

Several sites exhibit much lower seasonality in the OP values, especially traffic sites (RBX and STG-cle), the urban traffic site (NIC), or the industrial one (PdB). The lack of seasonality for some sites excludes the hypothesis of the OP being driven only by synoptic meteorological parameters such as sunshine or temperature, as it would impact all the sites similarly. It is clearly the difference in PM chemical compounds and reactivity, together with the timing of emission, that induces the seasonality of OP values when it is observed.

For comparison with previous studies, the Spearman correlation between chemical species, source contributions and OPs is also reported in Appendix  of this paper.

3.3 Results of OP models for the main PMF sources

3.3.1 Model accuracy and linear limitation

The MLR statistical validation was carried out by a residual analysis between the OP observed and the OP reconstructed by the model. For this evaluation, the intrinsic OP of the sources was set to the mean of the 500 bootstrap values. All but two sites present a very good correlation between observed and reconstructed OP ($r^{\mathrm{2}}>\mathrm{0.7}$) and a regression line close to unity (see Table S4; details and individual scatter plots are given at http://getopstandop.u-ga.fr/results?component=op_obsvsmodel, last access: 27 July 2021). We therefore consider our models valid, and each intrinsic OP (i.e., coefficient of the regression) may be explored individually to geochemically explain the observed OP.

However, despite our models being able to reproduce most of the observations with normally distributed residuals, it also tends to underestimate the highest values, and the residuals are often heteroscedastics (i.e., the higher the values, the higher the uncertainties). Then, the underlying hypothesis of linearity between endogenous variables (PM concentration of the sources) and exogenous variables (OPs) may be deemed invalid. It is also important to note that nonlinear processes are strongly suspected of the source apportionment of OP, as already noted by or and . As a result, future development on OP apportionment models should focus on this suspected nonlinearity, either by introducing covariation terms or using nonlinear models such as neural networks, for instance .

3.3.2 Intrinsic OP of the main PMF sources

Even if the models reproduce the observations correctly, this does not guarantee that the geochemical meaning extracted is the same for each of the models; i.e., the intrinsic OPs of the sources may completely differ from site to site. The question is then to identify whether a given source contributes similarly to the OP at all the sites. In other words, do all the models extract any general geochemical information relative to the OP?

Figure  presents the intrinsic OP${}^{\text{DTT}}$ and OP${}^{\text{AA}}$ for the selected subset of sources in terms of median and interquartile range (Q1 and Q3) to investigate the general information given by the different model. The values of mean and standard deviation and details per station for all sources are given in Tables S5–S7.

Figure 3

Median and interquartile range (solid bars and values in the figure) and mean (cross) of intrinsic OP${}^{\text{DTT}}$ and OP${}^{\text{AA}}$ values for the sources identified for at least two-thirds of the site (i.e., 10 series). The number of data considered are $n\times N$, with $n$ the number of series where the source was identified and $N=\mathrm{500}$ bootstraps: Road traffic: 7000, biomass burning: 7500, dust: 6500, primary biogenic: 7500, nitrate-rich: 7500, sulfate-rich: 7500, MSA-rich: 7000, aged salt: 6000. The values of means and standard deviations for each site are given in the Supplement.

[Figure omitted. See PDF]

For most of the sources identified, positive intrinsic OP values were observed considering the whole dataset. Small negative median values were only observed for the MSA-rich factor for the AA assay associated with large variance. Such results highlighted again that airborne particles, whatever their sources of origin, have significant oxidative properties.

We also observed a clear distinction between the intrinsic OP values for the different PM sources, with the median ranging from 0.04 to $\mathrm{0.21}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ for the OP${}^{\text{DTT}}$ and from $-\mathrm{0.00}$ to $\mathrm{0.17}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ for the OP${}^{\text{AA}}$. Such results agree with previous studies reporting different reactivities (or intrinsic OPs) for different sources based on receptor-model techniques .

Concerning the difference between the two tests, a road traffic source is the most reactive source towards OP${}^{\text{DTT}}$, with a value of about $\mathrm{0.21}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$. Such a value is almost twice as high as the ones observed for other significantly reactive sources, namely, primary biogenic sources (0.13 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$), biomass burning (0.13 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$), MSA-rich (0.12 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$), mineral dust (0.11 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$) and sulfate-rich (0.09 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$). Based on the intrinsic OP${}^{\text{AA}}$ results, a small number of PM sources shows significant redox activity, as already pointed out in previous studies and references therein. Only the biomass burning and road traffic sources show intrinsic OP values significantly higher than 0 (0.17 and 0.14 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$, respectively). Overall, the OP${}^{\text{DTT}}$ is sensitive to more sources than the OP${}_{\text{AA}}$, as already pointed out by and , and seems to target all the sources containing either metals or organic species. However, it is not sensitive to the nitrate-rich source.

We then confirm what previous studies found for these two assays, either by direct OP measurements at the emission source or by source apportionment. It is however hard to directly compare the absolute values from our results to the literature since the measurement protocols used are highly variable from one study to another.

3.3.3 Profile of intrinsic OP sources

The coefficients of variation (CVs, standard deviation over the mean) of the intrinsic OPs are lowest for the biomass burning and primary road traffic for the DTT assay, with values of 0.50 and 0.38, respectively, as well as for the AA assay with values of 0.52 and 0.67, respectively. However, since several intrinsic OPs are close to 0, the interquartile range (IQR) is also used to quantify the variability of intrinsic OP.

Biomass burning

The variability of the biomass burning intrinsic OP is somewhat site-dependent, with a low uncertainty within a given site but with slightly different intrinsic OPs between sites. This suggests that the variability is not linked to uncertainties of the MLR but may be due to actual local variations of the chemical composition of this profile.

This result contrasts with the fact that the biomass burning was identified as displaying relatively homogeneous chemical profiles between site profiles, with PD $<$ 0.1 and SID $<$ 0.7 (Fig. S2). Actually, the variability of OP-intrinsic values may be related to species that are not directly measured in our dataset. That is, polycyclic aromatic hydrocarbons (PAHs), oxy-PAHs, OH-PAHs, nitro-PAHs and especially polyaromatic quinones could not be measured at each of the sites used in the present study. These compounds are known to generate ROS and to contribute to the OP , have short lifetimes and are heavily influenced by the climatic condition . Similarly, the importance of the HULIS fraction within different biomass burning emissions might play a role in OP variability.

Road traffic

For the DTT assay, the uncertainty of the road-traffic-intrinsic OP at each site lies in the uncertainties of the other sites (but one, CHAM). Then, even if known with some uncertainties, these sources behave similarly at a large scale with regard to the OP${}^{\text{DTT}}$ and represent the most redox-active source.

For the OP${}^{\text{AA}}$, even if the CV is one of the lowest due to the high intrinsic OP of this source on average, the IQR is high. The different MLR models do not agree on a common intrinsic OP, with values ranging from 0.02 to 0.35 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$, without clear distinction by typology or groups of sites.

The chemical profile of the road traffic presents noticeable variability between sites, and no distinction between exhaust and non-exhaust emission is available with the current PMF analysis. Although all the profiles agree with a high loading of metals (notably Cd, Cu, Fe, Mo, Sb, and Sn) and organic carbon, the exact amount per microgram of PM differs from site to site (http://getopstandop.u-ga.fr/results?component=pmf_profiles, last access: 27 July 2021).

Also, the higher variability of the AA test outputs, compared to the DTT, suggests that some chemical species that are not measured here may influence more the OP${}^{\text{AA}}$ than the OP${}^{\text{DTT}}$ of vehicular emissions, as for biomass burning aerosols. These unmeasured species might vary in terms of loading between sites due to the different traffic-related sources, notably the extent of SOA or that of exhaust and non-exhaust emissions, which could have a high impact on OP, as shown by .

Nitrate-rich

The nitrate-rich factor, one of the most stable factors in terms of chemical fingerprint, is mainly composed of nitrate and ammonium, and very low amounts of OC or metals are present. At some sites, sulfate is also partly apportioned in this factor. Since nothing of this species is known to contribute to the OP, this factor presents logically almost null intrinsic OP for both assays.

Interestingly, this factor correlates fairly well with OPDTTv ($r=\mathrm{0.43}$) but is associated with the lowest intrinsic OP${}^{\text{DTT}}$ value (0.04 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$). This emphasizes the need to apply a relevant technique for OP source apportionment, beyond simple correlation.

We then confirm previous analysis of the low impact of secondary inorganic aerosols (SIA), at least for the nitrate part, on OP measurement in favor of the subsequent moderate role of SIA in human toxicity .

Sulfate-rich

Similarly to the nitrate-rich factor, this secondary inorganic factor displays rather stable chemical composition at all the sites. It is mainly composed of sulfate and ammonium but also organic carbon (between 5 $\mathrm{\%}$ and 33 $\mathrm{\%}$ of total OC${}^{*}$) and Se. The latter two species present some site-to-site discrepancies, and it should be mentioned that, when adding complementary organic markers such as 3-MBTCA to some sets of PMF input data, the sulfate-rich factor has been shown to split into two factors, revealing a biogenic SOA factor . Without these molecular markers, such mixing in the present study might explain a fraction of the variability of the OC apportioned and then of the intrinsic OP.

Interestingly, this factor does present an intrinsic OP for the DTT assay but almost none for the AA assay. Then, most probably the OC${}^{*}$ apportioned in this factor affects only the DTT assay and not the AA one, which contrasts with the biomass burning and road traffic sources (see above).

Mineral dust

The mineral dust displays a rather low variability when considering all datasets, but a deeper analysis showed two groups of sites: AIX-RBX-VIF vs. all other sites. The first group presents a high intrinsic OP for both assays (around 0.25 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ OP${}^{\text{DTT}}$ and 0.16 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ from OP${}^{\text{AA}}$), whereas sites from the second group display halved (around 0.1 $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$ for OP${}^{\text{DTT}}$) or almost null (OP${}^{\text{AA}}$) intrinsic OPs.

A first conclusion is that 80 $\mathrm{\%}$ of the sites present a common intrinsic OP for the dust source, hence the low IQR observed for this source. This factor apportions some OC${}^{*}$ but most of the ${\mathrm{Ca}}^{\mathrm{2}+}$, Al, Co, Cs, Mn and Sr. Then, none of these species appear to present intrinsic redox activity for the AA assay. It is not clear however whether the DTT is reacting with the organic species, with the metals, or both.

The differences between the two groups of sites may be further explained by difference in the chemical profile of the PMF factors. Indeed, since the dust factor at VIF highly differs from the other dust factors with PD $>$ 0.75 (mostly due to a very low amount of OC${}^{*}$) when compared to the other sites, we expect to have different intrinsic OPs. The biggest difference between VIF and other sites for the dust factor relies on the very low amount of PM apportioned in this factor, linked with important uncertainties. It is then expected that this factor will present an uncertain intrinsic OP for this source at this site.

For the sites of VIF, AIX and RBX, it appears that they also have higher loading in metals compared to other sites. Also, all three sites have important uncertainties associated with different metals (notably Cu, Sn or V; see http://getopstandop.u-ga.fr/results?component=pmf_unc, last access: 27 July 2021). Finally, it is also hypothesized that unmeasured trace species might be responsible for these differences, possibly coming from road dust resuspension and/or secondary processes leading to oxy-PAHs or HULIS .

MSA-rich (secondary organic aerosols)

It has been shown that biogenic SOA species generally contribute to the OP${}^{\text{DTT}}$ or to the generation of ROS . In our study, the MSA-rich factor is the only one strictly included in this category. However, the intrinsic OP of the MSA-rich source presents high variability between sites, with CVs of 3.1 and 7.8 for the DTT and AA assays, respectively, and the highest IQR, all sites combined, with no consistent similarity by typology or proximity. For instance, GRE-cb and VIF are within 15 $\mathrm{km}$ with similar sampling periods but present, respectively, the highest and lower intrinsic OP${}^{\text{DTT}}$ for this factor. We also observe important variability within each MLR model result at every site. This secondary organic source appears to be the most variable source in terms of intrinsic OPs, notably for the DTT assay.

The MSA-rich factor is identified thanks to a single marker (MSA), and very few studies only have reported it so far . As a result, we do not know for instance the loading of HULIS, quinone or isoprene-derived compounds contained in this factor nor the amount of aging it presents at each site. The exact primary sources or processes leading to this factor are still under discussion for non-Arctic regions . Hence, these uncertainties in the additional chemical compounds included in this factor, despite the excellent tracer capability of the MSA itself, might explain the diversity of its observed intrinsic OPs.

Moreover, this factor represents a low amount of PM mass. Its important uncertainties apportioned by the different PMF for this factor (between 0.7 $\mathrm{\%}$ and 5.5 $\mathrm{\%}$ of the total ${\mathrm{PM}}_{\mathrm{10}}$ mass) make it unclear whether the variability is mainly due to the sources' chemical component variability or to the inversion method used. Hence, we prefer not to conclude a clear effect of this source on the OP.

Primary biogenic organic aerosols

The primary biogenic source, mainly traced by polyols, presents some variability for the OP${}^{\text{DTT}}$. highlighted that spores or bacteria do contribute to the OP${}^{\text{DTT}}$ and OP${}^{\text{AA}}$ activities, even when the microbial cells are dead. However, the authors also present the inhibition of the DTT loss rate in the presence of 1,4-naphthoquinone or Cu. The presence of both synergistic and antagonistic effects between species and microbiota might explain the variability of intrinsic OP${}^{\text{DTT}}$ observed in Fig. , reflecting the different local microbiology carried by the PM, or covariations of the primary biogenic source with other metals or quinone-rich sources for instance.

Another hypothesis to explain the variability of the intrinsic OPs might be the “aging” of this factor, since pointed out that some secondary species may be incorporated in this factor at some sites, possibly making it a mix of primary biogenic and SOA (see also ). SOA mixed in the primary biogenic factor may then influence the intrinsic OP of this factor in different ways, similarly to our hypothesis for the MSA-rich factor.

Aged salt

Being mostly composed of Cl${}^{-}$, ${\mathrm{NO}}_{\mathrm{3}}^{-}$, ${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}-}$, Na${}^{+}$ and Mg${}^{\mathrm{2}+}$, this factor barely reacts with anti-oxidants (in both the DTT and AA assays). However, the variability observed for this factor is suspected of coming from the aging process, leading to the aggregation of organic compounds (between 0 $\mathrm{\%}$ and 9 $\mathrm{\%}$ of the total OC${}^{*}$ apportioned in this factor).

Local PMF factors

Some other factors were obtained in a limited set of PMF outputs, and their chemical profiles are somewhat variable. It is however interesting to discuss their local impact on OP.

The HFO source is identified at MRS-5av and PdB, both sites being located close to harboring activities in large ports on the Mediterranean coast. It presents an intrinsic OP${}^{\text{DTT}}$ of $\mathrm{0.51}\pm \mathrm{0.14}$ and $\mathrm{0.21}\pm \mathrm{0.04}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$, respectively, and an intrinsic OP${}^{\text{AA}}$ of $\mathrm{0.04}\pm \mathrm{0.02}$ and $\mathrm{0.11}\pm \mathrm{0.03}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$, respectively. The intrinsic OP${}^{\text{DTT}}$ is then on average higher than the road traffic one, making HFO the second contributor of the daily mean and median source contribution at MRS-5av for the OP${}^{\text{DTT}}$ contribution and the fourth one for the OP${}^{\text{AA}}$ contribution (see the website). For the PdB site, the contributions are a bit lower. Although only two sites presented an HFO factor, and similarly to previous studies , it suggests that the PM that originated from this source may significantly contribute to the total OP in and around harbor cities.

For the industrial factor, and as already stated, the chemical composition profile highly varies from site to site for the six sites where it is determined. As a matter of fact, the intrinsic OPs of this profile are also highly variable. The intrinsic OP${}^{\text{DTT}}$ are high for GRE-cb and GRE-fr_2017 ($\mathrm{0.52}\pm \mathrm{0.30}$ and $\mathrm{0.37}\pm \mathrm{0.27}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$, respectively), as are the intrinsic OP${}^{\text{AA}}$ ($\mathrm{0.82}\pm \mathrm{0.29}$ and $\mathrm{0.61}\pm \mathrm{0.17}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}\mathrm{\mu }{\mathrm{g}}^{-\mathrm{1}}$, respectively). However, both are close to 0 for the other sites where this factor is found (PdB, AIX, TAL and VIF). The high intrinsic OP again suggests the role of metals in the OP of PM; however, since this factor has strong uncertainties associated with the PMF results and then with the intrinsic OP, further work on the source profiles is needed to draw firmer conclusions.

Overall geochemical agreement

Finally, despite the differences within PMF solutions from one site to another, including the slightly different number of source factors and the different OP signals, the rather low variability of the intrinsic OP determined for a given source suggests that most of the sources of PM behave coherently with regards to the OP over a large geographical area in France. It then supports the idea that, at the national scale, the sources described above have a rather stable intrinsic OP. We also stressed the importance of the underlying source apportionment method (here the PMF) and the investigation of its uncertainties in order to interpret the variabilities of the intrinsic OP.

3.4 Contribution of the sources to the OPs

The relative importance of the contributions of the sources to the total ${\mathrm{PM}}_{\mathrm{10}}$ OPs is weighted by their different intrinsic OPs. The question, therefore, is to what extent the contributions of the sources to the OPs differ from their contributions to the PM mass concentration. In this part, we present an aggregated view of the seasonal contribution of the sources to the OP in Fig.  and the daily contribution in Figs.  and , considering all sites. Details per sites are presented in the website provided herewith as the Supplement.

Figure 4

Mean monthly contribution of the eight main sources to the (a) ${\text{PM}}_{\mathrm{10}}$ mass, (b) ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and (c) ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, taking into account each source contribution of every site; their respective normalized contributions are presented in (d) ${\text{PM}}_{\mathrm{10}}$ mass, (e) ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and (f) ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$.

[Figure omitted. See PDF]

Figure 5

Averaged daily contribution of the sources to (a) the ${\text{PM}}_{\mathrm{10}}$ mass, (b) the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and (c) the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$. The bars represent the mean, and the error bars represent the 95 $\mathrm{\%}$ confidence interval of the mean.

[Figure omitted. See PDF]

Figure 6

Median daily contribution of the sources to (a) the PM mass, (b) the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and (c) the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$. The bars represent the mean and the error bars represent the 95 $\mathrm{\%}$ confidence interval of the median.

[Figure omitted. See PDF]

We would like to stress here that our dataset included an important proportion of Alpine sites as well as urban sites. Then, the extrapolation to the whole of France or to other regions of western Europe should be done cautiously.

3.4.1 Seasonality of the contribution by mass or OPs

As already shown by some previous studies in France , the seasonal mean contributions to the ${\mathrm{PM}}_{\mathrm{10}}$ mass show the importance of the biomass burning source, followed by secondary inorganic (sulfate-rich and nitrate-rich), dust and road traffic. As a direct consequence of the different intrinsic OPs for these sources, we do observe a redistribution of their relative importance for the total OP. That is, the nitrate-rich source that may contribute a significant amount to the ${\mathrm{PM}}_{\mathrm{10}}$ mass, notably in spring, barely contributes to the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ or to the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, in agreement with . Conversely, the road traffic contributes about 15 $\mathrm{\%}$ during summer to the mean ${\mathrm{PM}}_{\mathrm{10}}$ mass but represents more than 50 $\mathrm{\%}$ of the mean ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ in the same period (Fig. d–f). However, the biomass burning sources are still major contributors to both the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ during the winter months. We note that the primary biogenic source also contributes to the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ but to a lesser extent to the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$. Finally, the dust source is an important contributor to the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ but not to the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$. These results confirm and extend what previous studies have already found . Overall, the main contributors to OPs are the three factors suspected of including anthropogenic SOA (biomass burning, road traffic and dust, including possible resuspension of road dust for the latter one). It follows that, considering the seasonality of OPs, regulations should target the biomass burning emission in order to decrease the ${\mathrm{PM}}_{\mathrm{10}}$ OPs during winter by a large amount but also the road traffic that contributes homogeneously to both OPs around the whole year.

3.4.2 Daily mean and median contribution: insights for exposure

A concise view of the same results, this time on the daily-aggregated basis, is provided in Figs.  and , presenting the contributors to ${\mathrm{PM}}_{\mathrm{10}}$ mass and OPs, ranked in decreasing order. Figure  reports the typical “mean” daily value, a parameter generally used in the atmospheric community, while Fig.  reports the “median” daily value, often used in epidemiological studies in order to discuss the chronic exposure of the population and avoid the high impact of infrequent low or high events that highly influence the mean value.

Due to the non-normality of the contribution and the high contributions of some sources at some sites (for instance, the biomass burning source in Alpine valleys), the results highly differ when considering the mean or the median contributions, and the two statistical indicators may not address the same questions (the mean is more related to the identification of the major source contributors and the median to the exposure of the population). The skewness of the distribution is not surprising as some high-${\mathrm{PM}}_{\mathrm{10}}$ events (i.e., a short time event in the dataset) were present in our measurements. This is also specifically anticipated in Alpine areas (CHAM, PAS, MNZ, VIF, GRE-cb, GRE-fr_2013 and GRE-fr_2017), where the frequent development of atmospheric thermal inversion layers in winter is causing increased pollutant concentrations.

We observed (Fig. ) a redistribution of the daily mean contribution sources' rank between the ${\mathrm{PM}}_{\mathrm{10}}$ mass, ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ similarly to the monthly mean contribution discussed above. The biomass burning source, being an important contributor to the ${\mathrm{PM}}_{\mathrm{10}}$ mass, contributes also significantly to both OP and is ranked as the first contributor to both ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ mean daily contribution (mean $\mathrm{0.51}$ and $\mathrm{0.72}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$, respectively). The road traffic source contribution, due to its highest intrinsic OP in both assays, presents almost the same daily mean contribution as the biomass burning for the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and is the second contributor to the daily mean ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, with half the contribution of the biomass burning (mean $\mathrm{0.50}$ and $\mathrm{0.34}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$, respectively). The other sources barely contribute to the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ ($<\mathrm{0.1}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$). For the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$, the dust is the third contributor (mean $\mathrm{0.34}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$), followed by sulfate-rich and primary biogenic ($\mathrm{0.23}$ and $\mathrm{0.16}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$, respectively). The nitrate-rich, aged sea salt and MSA-rich present a low contribution (mean $<\mathrm{0.1}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$) due either to their low intrinsic OP or to their low contribution to the PM mass.

However, for the daily median contribution, due to the high seasonality of the biomass burning source and the consistent contribution throughout the year of the primary road traffic, sulfate-rich and dust sources, the ranks of the sources are drastically redistributed between the three metrics considered (Fig. ). Moreover, the absolute values of the contributions are also lowered compared to the mean daily contribution due to low frequencies of highly loaded PM events. The major source contributing to the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ is now the primary road traffic (median $\mathrm{0.36}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$), contributing more than twice as much as the second source, namely, the dust one (median $\mathrm{0.16}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$), followed by sulfate-rich (median $\mathrm{0.13}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$) and then biomass burning (median $\mathrm{0.11}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$). For the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$, the two dominant sources are the primary road traffic (median $\mathrm{0.29}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$) and the biomass burning (median $\mathrm{0.19}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$), all other sources being now negligible (aged salt $\mathrm{0.016}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$, primary biogenic $\mathrm{0.015}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$ and the others contributing less than $\mathrm{0.01}$ $\mathrm{nmol}{\min }^{-\mathrm{1}}{\mathrm{m}}^{-\mathrm{3}}$).

The high differences between the mean and median contributions could have strong implication for air quality policies. Indeed, as previously shown, the biomass burning may contribute more than 50 $\mathrm{\%}$ of the high OPs during winter, and even more for some days. However, such events do not represent well the daily exposure of the population over the full year. Even if the regulations should target those events to prevent acute exposure, they should also strongly take into account the long-term exposure to a lower but constant level of pollutant, since there is no threshold below which PM is no longer harmful . In this respect, the emissions from the road traffic become a major concern as well, showing that this source might actually be the most important one to be targeted in order to decrease adverse health effects due to chronic exposure to PM pollutants.

4 Limitations of the study

In this study, we focused on major sources and trends and hence limited our study to some aspects. Notably, the PMF standardized approach allows common source identification at the national scale but may also damper some site specificities. Also, the choice to focus on the main sources of PM and to discuss the aggregated results shortens the discussion on some local specificity, notably potential local sources that contribute to the OP (for instance HFO or industry) and that may be relevant for some sites but that are not applicable to a wider area.

One main limitation is also the use of linear regression tools, whereas it has been shown that OP is not fully proportional to the mass of compounds. The residual analysis seems to agree with this experimental finding since the highest OP samples are underestimated by the MLR model. The addition of a covariation term or even the use of nonlinear regression may be the next step to better explain the OP of the sources .

Moreover, if the intrinsic OP results from the MLR can be extrapolated to any given site with a similar regional background to the urbanized area used in this study, the source contribution extrapolation should be taken cautiously since our dataset displays an over-representation of the Alpine sites with regard to the whole of France.

5 Conclusions

To our knowledge, this study gathers the most important database of OP samples, with concomitant observations of chemistry analysis, source apportionment through PMF, and the measure of two OP assays (DTT and AA) for 15 yearly time series over France spanning between 2013 and 2018 for a total of $>\mathrm{1700}$ samples.

We demonstrated that source apportionment of OP through a “simple” multilinear regression without any constraint on the coefficient provides good statistical results and can explain the observed ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$.

• The intrinsic OPs of the main regional sources present values in the same range at each site, especially for the primary road traffic, biomass burning, nitrate-rich, dust, and sulfate-rich PMF factors. Biogenic and MSA-rich factors present higher discrepancy according to the site together with the highest uncertainties at each site.

• Some site-specific sources might have an important intrinsic OP and can account for a non-negligible part of the observed OP (notably in harbor cities or near industrial sites).

• Different sensitivities for the two OP assays to a given source are highlighted. The DTT appears to be sensitive to a wide range of sources, whereas the AA targets mainly the biomass burning and primary road traffic factors.

• With consistency at the regional scale, the primary road traffic and biomass burning factor are the main absolute OP contributors, together with dust and secondary sulfate for OP${}^{\text{DTT}}$ to a lesser extent. Conversely, the nitrate-rich factor barely contributes to OPs.

• In order to assess the chronic population exposure, the median daily contributions of sources to the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ are also reported and present important differences in ranking compared to the mean contributions. The importance of the primary road traffic source drastically increases, notably for the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$, whereas the biomass burning contribution is lowered. However, only the road traffic and biomass burning sources contribute to the daily median of the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$.

Finally, the relatively stable intrinsic OP at a large geographical scale for the main PM sources allows future work on the implementation of the OP in the regional chemistry transport model. This step would allow a quantitative estimation of the population exposure OP, expanding potential cross-over studies with epidemiology.

Appendix A Correlation between OP, chemical species, and sources

The Spearman correlations between the chemical species and ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ are presented in Fig. , while the Spearman correlations between the source mass apportionment from the PMF and the measured ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ are presented in Fig. . All samples from all sites were considered in the results presented here.

First, we can see that both OP assays correlated well with each other ($r_{{\text{OP}}^{\text{DTT}}-{\text{OP}}^{\text{AA}}}=\mathrm{0.61}$) but do not present an exact similitude. Notably, the OP${}_{\mathrm{v}}^{\text{AA}}$ has stronger seasonality than OP${}_{\mathrm{v}}^{\text{DTT}}$, and higher correlations are found during winter than summer. Details of the individual time series are given on the website.

Second, the only source that strongly correlates with one OP ($r>\mathrm{0.6}$) is the biomass burning with the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$. Some low to mild correlations ($\mathrm{0.3}<r<\mathrm{0.6}$) are found for the ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ vs. road traffic, biomass burning, nitrate-rich and dust and for the ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ vs. nitrate-rich and road traffic. The copper, mostly apportioned by the road traffic, is the most correlated metal with the OP${}_{\mathrm{v}}^{\text{AA}}$ and the second for OP${}_{\mathrm{v}}^{\text{DTT}}$. For the levoglucosan and mannosan, strong correlations are found with the OP${}^{\text{AA}}$. These results are in agreement with previous studies, either with the source correlation or with the proxies of sources (namely, levoglucosan for biomass burning and EC, and iron, copper or PAH for road traffic) (Fig.  and ).

Figures  and  also indicate that the nitrate and nitrate-rich source concentrations mildly correlate with both OPs. However, neither nitrate nor the chemical species included in the chemistry profile of the nitrate-rich source present redox-active capabilities. Conversely, crustal elements (mainly Ti and Ca${}^{\mathrm{2}+}$) present no to low correlation with OPs. Although mineral dust has been reported to contribute to the OP, the episodic event of Saharan wind may be insufficient to be reflected in the simple $r^{\mathrm{2}}$ value. Also, as already stated, the strong seasonal cycle leads to negative correlation between some sources or species and the OP (MSA or polyols species and MSA-rich and primary biogenic factors). When considering only the warm period, they appear positively correlated with both OPs due to the exclusion of the strong impact of the biomass burning source (not shown here).

These examples indicate that it is hard to assess robustly the links between OPs and chemistry when using only correlations. Moreover, even a good correlation may not reflect any causality, and the multilinear regression should disentangle possible covariation due to meteorological effects (accumulation or long-range transport).

Figure A1

Spearman correlation coefficients between ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ and ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ and the chemical species observed at each site. The numbers of samples are specified in parentheses.

[Figure omitted. See PDF]

Figure A2

Spearman correlation coefficient between ${\text{OP}}_{\mathrm{v}}^{\text{DTT}}$ or ${\text{OP}}_{\mathrm{v}}^{\text{AA}}$ and the different PM sources identified for at least two-thirds of the sites. All the sites are merged, and whole time series of measurements were considered. The number of observations considered for the different sources are as follows: aged salt: 1430; biomass burning: 1700, dust: 1489, MSA-rich: 1595, nitrate-rich: 1700, primary biogenic: 1700, road traffic: 1587, sulfate-rich: 1524.

[Figure omitted. See PDF]

Code and data availability

Data available on request at [email protected]. Website code available at https://gricad-gitlab.univ-grenoble-alpes.fr/pmall/app_op (last access: 27 July 2021) (commit fb99ed42e3f704464e88deca9548b74128fd7b24) .

The supplement related to this article is available online at: https://doi.org/10.5194/acp-21-11353-2021-supplement.

Author contributions

OF, JLJ, GU and JLB were in charge of the coordination of different research programs and funding acquisitions. DS did the data curation and ran the PMF for the SOURCES program, FC and JA did the data curation and ran the PMF for the DECOMBIO program, and SW and LJSB did the data curation and ran the PMF for the MobilAir program. AC and GU set up the two OP assay methodologies. SW designed the methodology, did the formal analysis and prepared the present manuscript and figures. GU and JLJ designed, reviewed and edited the first draft of the manuscript. All the co-authors read and edited the manuscript.

Competing interests

The authors declare that there is no conflict of interest.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

This work was partially funded by ANSES for OP measurements (ExPOSURE program, grant 2016-CRD-31), an IDEX UGA grant for innovation 2017 ROS-ONLINE and CDP IDEX UGA MOBILAIR (ANR-15-IDEX-02). It was also supported by the French Ministry of Environment as part of the national reference laboratory for air quality monitoring (LCSQA, program CARA), for some of the chemical analysis related to the GRE-fr, TAL, RBX, STG-cle, MRS-5av and NIC stations. The study in CHAM, MNZ and PAS was funded by ADEME and PRIMEQUAL within the DECOMBIO program (1362C0028). The studies in GRE-cb and VIF were funded by the UGA IDEX MobilAir program (ANR-19-CE34-0002-01) and the ADEME QAMECS program (1262c0011). Regional monitoring networks (namely, Atmo AuRA, Atmo Sud, Atmo HdF, Atmo NA and Atmo GE) financially contributed to sampling and/or chemical analyses for samples from their respective sites. The PhD of Samuël Weber was funded by a grant from ENS Paris. This study was also supported by direct funding by IGE (technician salary), the LEFE CHAT (program 863353: “Le PO comme proxy de l'impact sanitaire”), and LABEX OSUG@2020 (ANR-10-LABX-56) (for funding of analytical instruments).

The authors wish to thank all the numerous people (who could not be listed exhaustively here) from the different laboratories (IGE and Air-O-Sol analytical platform, Ineris and PTAL analytical platform of EDYTEM) and from the regional air quality monitoring networks listed as co-author affiliations, who actively participated in filter sampling and/or analysis. The Institute Mines Telecom (IMT-Lille Douai, PI: Laurent Y. Alleman) and the Laboratoire des Sciences du Climat et de l'Environnement (LSCE, PI: Nicolas Bonnaire) are also warmly thanked for their contribution to chemical analysis of samples collected at RBX, NGT TAL and STG-Cle.

Financial support

This research has been supported by the Agence Nationale de Sécurité Sanitaire de l'Alimentation, de l'Environnement et du Travail (grant no. 2016-CRD-31), the Université Grenoble Alpes (grant nos. ROS-ONLINE and ANR-19CE34-0002-01), the Agence Nationale de la Recherche (grant nos. ANR-15-IDEX-02 and ANR-10-LABX-56), and the Agence de l'Environnement et de la Maîtrise de l'Energie (grant nos. 1362C0028 and 1262C0011).

Review statement

This paper was edited by Manabu Shiraiwa and reviewed by two anonymous referees.