Overview of the French operational network for in situ observation of PM chemical composition and sources in urban environments (CARA program)

The CARA program has been developed since 2008 by the French reference laboratory for air quality monitoring (LCSQA) and the regional monitoring networks to gain a better knowledge at the national level on the particulate matter (PM) chemistry and its diverse origins in urban environments. It results of strong collaborations with international-level academic partners, allowing to bring state-of-the-art, straightforward and robust results and methodologies within operational air quality stakeholders (and subsequently, decision makers). Here, we illustrate some of the main outputs obtained over the last decade thanks to this program, regarding methodological aspects (both in terms of measurement techniques and data treatment procedures) as well as acquired knowledge on the predominant PM sources. Offline and online methods are used following well-suited quality assurance and quality control procedures, notably including interlaboratory comparison exercises. Source apportionment studies are conducted using various Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 © 2021 by the author(s). Distributed under a Creative Commons CC BY license.

Abstract: The CARA program has been developed since 2008 by the French reference laboratory for air quality monitoring (LCSQA) and the regional monitoring networks to gain a better knowledge at the national level on the particulate matter (PM) chemistry and its diverse origins in urban environments. It results of strong collaborations with international-level academic partners, allowing to bring state-of-the-art, straightforward and robust results and methodologies within operational air quality stakeholders (and subsequently, decision makers). Here, we illustrate some of the main outputs obtained over the last decade thanks to this program, regarding methodological aspects (both in terms of measurement techniques and data treatment procedures) as well as acquired knowledge on the predominant PM sources. Offline and online methods are used following well-suited quality assurance and quality control procedures, notably including interlaboratory comparison exercises. Source apportionment studies are conducted using various air). These species account for a very small fraction of PM, while some others are found in much higher abundance. This is notably the case of the PM carbonaceous fractions -namely Elemental and Organic Carbon (EC and OC, respectively) -as well as some water-soluble inorganic species, notably including nitrate (NO3 -) , sulfate (SO4 2-), ammonium (NH4 + ), chloride (Cl -), sodium (Na + ) and calcium (Ca 2+ ) [8]. Actually, the contribution of these species to PM2.5 shall be monitored at a few remote stations of EU Member States -with no fixed limit nor target value yet -according to Directive 2008/50/CE (1 site per 100 000 km² in each country). This requirement has been elaborated to better assess transboundary air pollution phenomena, complementary to the co-operative Programme of the long-range transmission of air pollutants in Europe, also known as European Monitoring and Evaluation Programme (EMEP, www.emep.int, [9], [10]), under the UNECE Convention on Long-Range Transboundary Air Pollution (CLRTAP). This policy-driven program has been complemented with various research projects at the European scale over the last decades, notably leading to the current establishment of the Aerosol, Clouds and Trace gases Research Infrastructure (ACTRIS, www.actris.eu, [11], [12]). However, these initiatives are still mainly focused on background air pollution at a regional (or larger) scale, rather than on hotspots for population exposure to poor air quality (i.e., urban and/or heavily industrialized areas).
In this context, LCSQA has been implementing together with the MTE and AASQAs an ambitious operational network for in situ observation of urban PM chemical composition and sources (so-called CARA program, www.lcsqa.org/fr/le-dispositif-cara) since 2008 [13], with the following objectives:  Document the chemistry and main sources of PM pollution events;  Investigate peculiarities possibly leading to local PM10 threshold exceedances;  Achieve long-term (≥ 1 yr) source apportionment studies at various locations;  Provide datasets useful for the validation of chemical transport models (CTMs);  Trigger and/or support further research activities for a better knowledge on air quality fundamentals. These objectives, as well as related experimental strategies, are refined yearly to try and accommodate the needs of the MTE and AASQAs.
The present paper aims to describe the overall organization and activities of the CARA program, as well as to illustrate and/or summarize some of the main findings of the last decade. After a brief description of the methodologies in use, it presents major results obtained from some inter-laboratory comparison exercises, which are essential for data quality assurance and quality control in networking activities. A general picture of the PM chemical composition as obtained by the measurements of its main constituents is then proposed, after discussing in more details the prevalent primary PM sources that have been investigated so far in the frame of the CARA program, namely residential wood burning, road transport, primary biogenic particles, sea salts, and Saharan dust. Finally, we address the main origins of the secondary aerosols (formed in the atmosphere through oxidation processes and interactions with the gaseous phase) impacting French urban environments, i.e., ammonium nitrate, ammonium sulfate, and secondary organic aerosols (SOA).

PM10 sampling network
The list of past and current sampling stations included in the CARA program is given in Table  A1. These stations mainly correspond to (peri-)urban background stations, with however a few traffic or rural background sites. Except for specific measurement campaigns, most of these sites have been operated for quasi-continuous sampling over several years. Oldest ones -namely in Bordeaux, Lyon, Grenoble, Rouen, Strasbourg and Revin -have been collecting PM10 filter samples for more than 10 years. Here, results obtained from this dense observation network are supplemented by complementary results acquired during a specific pluriannual research program focusing on air pollution sources and mechanisms in the Arve valley of the French Alps (DECOMBIO, [14], [15]). 4 of 47 Data collected at the remote regional background site of the long-term Environmental Observatory (OPE, [16], [17]) -maintained by ANDRA -has also been added for purposes of comparison with results from urban environments.
Particles are collected onto quartz fiber filters (Whatman QMA in 2008-2010, and Pallflex TissuQuartz from 2010 onwards), pre-fired during 24 h at 500 °C before use, to lower filter blank values (especially for organic carbon concentrations). Filters are mainly collected by AASQAs following the requirements of the EN 12341 European standard, using sequential high-volume samplers (typically, DAH-80 models, Digitel, running at a flow rate of 30 m 3 /h), except for some intensive campaigns using low-volume samplers (e.g., Partisol 2025 model, Thermo, 1 m 3 /h, or SEQ47 model, Leckel, 2.3 m 3 /h). Quasi-continuous sampling is commonly performed in synergy with regulatory offline monitoring of PAHs, with continuous sampling of daily filters alternating between the latter measurements (every third or sixth day) and the CARA program filter collection (all the other days). Sampled filters are returned from station and then commonly sent to Ineris (or IGE, for Alpine and south-eastern sites) within cool boxes, and finally stored in freezers.
Only a small fraction (15-20%) of the collected filters is actually analyzed for its chemical composition, fitting the objectives described above. Examples of chemical analyses that have been performed in the frame of the CARA program can be found for instance in [18]. Each filter chemical analysis includes the measurement of the carbonaceous fractions -i.e., organic and elemental carbon (OC and EC) -as well as the main anions and cations (nitrate: NO3 -, sulfate: SO4 2-, chloride: Cl -, ammonium: NH4 + , sodium: Na + , potassium: K + , magnesium: Mg 2+ , and calcium: Ca 2+ ), as featuring major aerosol chemical species found in European ambient air, namely soot, organic matter, ammonium nitrate, ammonium sulfate, sea salt and mineral dust [8]. OC-EC analyses are performed following the EN 16909 European standard thermo-optical method (in transmission mode), using the so-called EUSAAR-2 protocol [19]. Anions and cations are measured by ion chromatography (IC), after extraction of filter sample punches in ultra-pure water and applying recommendations of the EN 16913 European standard. These IC analyses also allows for the quantification of some organic species, such as oxalate and methylsulfonic acid (MSA). Water-soluble sample extracts obtained from this analytical procedure are also commonly used for the measurement of some anhydrosugars considered as molecular markers of biomass burning emissions, including levoglucosan and its isomers, or sugars and sugar alcohol (glucose, mannitol, arabitol and sorbitol, the 3 last ones being often referred to as polyols) considered as indicators of primary biogenic emissions, using ion chromatography or high performance liquid chromatography followed by amperometric detection (IC-or HPLC-PAD). More recently, efforts have also been put in the measurements of cellulose, as a proxy of plant debris within the PM10 fraction [20]. For source apportionment purposes, the elemental composition (starting from sodium atomic mass) is additionally measured by inductively coupled plasma mass -or atomic emission -spectrometry (ICP-MS or ICP-AES) following sample preparation and analytical procedures defined by the EN 14902 European standard. Finally, various other analytical techniques, including liquid or gas chromatography coupled with mass spectrometry (LC or GC-MS), have also been used for a better apportionment of the organic aerosol fractions (e.g., [18], [20]- [23]).
Chemical analyses are conducted by Ineris and IMT Lille Douai as part of the LCSQA, as well as by academic partners (mostly IGE and LSCE). Each of these laboratories is regularly participating to inter-laboratory comparison (ILC) exercises at the international level, such as those performed in the frame of the CEN, EMEP, ACTRIS and/or GAW programs (e.g., [24]).
When using high-volume samplers, large circular filters (15 cm total diameter, corresponding to a sampled surface area of about 154 cm²) are collected, so that leftover filter surfaces generally remain available after sub-sampling for completion of the analyses mentioned above. These remaining filter parts are saved in freezers, with non-analyzed filters, contributing to enlarge the wide sets of PM10 samples that can be made available for possible further analyses (e.g., for oxidative potential measurements [25], or investigation of the links between PM chemistry and microbiology, [26]).

Online monitoring network
In such a sampling network, filter collection and transportation followed by subsequent offline chemical speciation and data treatment can be achieved within a relatively short timeframe (e.g., down to 3 days, [27]), but do not allow for the assessment of the nature of aerosol pollution episodes in nearreal time. Moreover, albeit allowing for extensive characterization of the PM chemistry at the molecular level, documentation of seasonal and inter-annual variabilities at multiple sites using such a methodology is very much labor-intensive and costly on the long-term basis. To tackle these issues, online and continuous measurements shall be envisaged, using simple but robust technologies. In this context, the LCSQA/Ineris started to evaluate and implement instrumental candidates for in situ automatic measurements of the aerosol chemical composition as soon as in the early 2010's, in collaboration with LSCE. In particular, two instruments appeared to be well-suited for these purposes [28], [29]: the multi-wavelength aethalometer (AE33 model, Magee Scientific) and the aerosol chemical speciation monitor (ACSM, Aerodyne Res. Inc.), which are briefly described hereafter.
A full description of the AE33 measurement principle is available elsewhere [30], [31]. Briefly, this instrument allows for the determination of the aerosol light absorption coefficients at seven different wavelengths (370, 470, 520, 590, 660, 880 and 950 nm). Ambient concentrations of equivalent black carbon (BC) are then estimated using data obtained at 880 nm and assumptions on BC mass absorption efficiency (MAE). Multi-wavelength measurements are also allowing for the discrimination between the two main BC sources, i.e., liquid fuel burning and solid fuel burning, using the so-called Aethalometer model [32]. Equivalent BC concentrations (eBC) -as well as its source apportioned fractions (named here by convention eBCff and eBCwb for liquid and solid fuel emissions, respectively) -which are presented and discussed in the present paper have been obtained in the PM2.5 fraction and using the default calculation parameters proposed by the manufacturer.
The ACSM has been designed to provide continuous measurements of the main non-refractory chemical species within submicron aerosols [33]. It is based on a less advanced but more robust technology than Aerosol Mass Spectrometry (AMS) nowadays widely used for state-of-the-art research experiments in atmospheric science. Briefly, fine aerosols are sampled into the ACSM system through a critical orifice and an aerodynamic lens, allowing to focus them -under high vacuumonto a vaporizer usually heated at a temperature of about 600 ˚C. Vaporized molecular fragments can then be detected using electron impact (70 eV) ionization mass spectrometry. Each ACSM operated in the frame of the CARA program has been equipped with a PM1 aerodynamic lens. Measured species notably include organic aerosols (OA), nitrate, sulfate, ammonium, and chloride. The determination of their ambient mass concentrations requires cautious calibrations of the instrument response factor and relative ion efficiency [34], which are predominantly performed together with the LCSQA. The consistency of these ACSM measurements at the national level is further assessed by regular participation to ILC exercises organized by the aerosol chemical monitor calibration center (ACMCC), as a unit of the ACTRIS topical centre for aerosol in situ measurements (see section 2.4).
The combination of AE33 and ACSM measurements is then allowing to determine both the refractory and non-refractory major chemical species within the fine aerosol mode. Following successful evaluation of the operability of this monitoring strategy at the ACTRIS-SIRTA [28], it could be applied to selected urban background stations operated by AASQAs. This deployment has been achieved gradually since then, leading to a current number of 12 stations being equipped with AE33 and ACSM in the frame of the CARA program (Table A1 and Figure 1). Besides SIRTA, 10 of these sites are part of 7 different AASQAs and the twelfth one corresponds to a new ACTRIS site operated by IMT Lille Douai in Lille. The current distribution of sites equipped with both ACSM and AE33 devices is then offering a relatively dense covering of the French metropolitan area. Considering the investment costs associated with ACSM equipment, and that it measures species displaying rather homogeneous background concentration at regional scale, this ACSM network nowadays offers very close-to-optimum number of sites. In contrast to ACSM, AE33 is measuring primary species which can have a larger spatial heterogeneity, so that additional urban areas have been equipped with AE33 only (Table A1 and Figure  1).
The LCSQA is coordinating the quality assurance and quality control (QA/QC) strategy related to this operational online monitoring program. This notably includes calibration of ACSM devices (as mentioned above), production and dissemination of guidelines and standard operating procedures (for both ACSM and AE33), training sessions, and ILC organization (see next sub-section). In return, feedbacks and expertise brought back on a regular basis by the AASQAs, by the manufacturers and the French ACSM seller (ADDAIR) are ensuring lively dissemination of the information between the CARA program partners and efficient improvement of the QA/QC protocols and procedures.

Data processing
At the European level, PM10 is known to be generally and predominantly composed of carbonaceous matter (i.e., elemental carbon, EC, and organic matter, OM), inorganic matter notably including nitrate, sulfate, and ammonium, sea salts and mineral dust [8]. Out of extremely rare exceptions, each site of the CARA program is equipped with PM10 (and also often PM2.5) automated measurement systems, so that quality control of the chemical speciation datasets includes a systematic chemical mass closure exercises [35], where the sum of major components (typically, PMchem = EC + OM + nitrate + sulfate + ammonium + sea salts + mineral dust) is validated by comparison with regulatory PM10 measurements performed by the AASQAs (e.g., [36]).
For a further investigation of the PM sources, receptor models are nowadays widely used in atmospheric sciences. In particular, a growing number of studies have been applying Positive Matrix Factorization (PMF, [37]) to chemically speciated datasets over the last decade (e.g., [38]- [40]). PMF is a multivariate tool which makes use of measured concentrations and their uncertainties to solve the mass balance equation X = G×F+E. where X corresponds to the input data matrix, F represents a matrix whose vectors stand for the profiles of a p number of factors (that might eventually be Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 attributed to specific emission sources or secondary aerosol types), G corresponds to the contributions of the p sources to total PM, and E represents the residual matrix. Detailed and practical descriptions of the way PMF has been used in the frame of the CARA program, and this program participate in the development of receptor models, can be found elsewhere [18], [36], [41]- [43]. All the PMF results presented in this paper were obtained using the most recent US EPA PMF v5.0 software and following recommendations of the EU guidance document [44]. Source apportionment of PM in near-real-time is a tremendous challenge. ACSM (and AMS) measurements commonly allow for discriminating between different types of submicron organic aerosols, based on PMF approaches applied to OA mass spectra ( [45], [46]). Despite recent advances and current developments towards automatic source apportionment analysis [47], applying receptor models on ACSM datasets is still highly labor-intensive. Moreover, while outputs are quite interesting for atmospheric sciences, they might be misleading for general public and policy makers when trying to link OA components and actual emission sources (e.g., [48], [49]). For these reasons, the application of near-real time source apportionment of ACSM datasets within such a network as the CARA program is still requiring further extensive research and harmonization efforts (see also sections 5 and 12).
By contrast, AE33 can provide direct information on solid fuel and liquid fuel combustions. Indeed, for French urban environments, eBCff has been shown to be predominantly related to road transport emissions, while eBCwb is mainly originating from wood burning for residential heating in cold seasons [50], [51]. These fractions can then be used to estimate PM loadings that may be attributable to both sources (i.e., PMff and PMwb, respectively), using conversion factors such as: PMff = a x eBCff and PMwb = b x eBCwb, where a and b can be retrieved from the literature and/or further site-specific a priori knowledge. In practice, a is considered here as constant over the metropolitan territory (with a = 2, roughly corresponding to OC/EC ratio of 0.5-0.6 and hydrocarbon-like organic aerosol O-to-C ratio of 1.2-1.4 within vehicular exhaust, [52], [53]) whereas fit-for-purpose b values have been determined at each site through comparisons with independent offline filter-based source apportionment analyses [41], [54]. For stations where no independent source apportionment study allowed to determine a site-specific b value, the LCSQA currently recommends to AASQAs estimating wintertime PMwb from AE33-based brown carbon measurement and using a constant conversion factor, as explained in [55] and [56]. We do not pretend in this paper that these simple calculations could be accurately applied in other countries, and overall uncertainties on the order of 30-40% can be considered here for PMff and PMwb loading estimates. Yet, they proved to compare well with independent methods (examples given for PMwb in section 5) and offers very useful information to AASQAs in near-real time, especially during pollution episodes (section 4.3). Finally, it shall also be underlined that PMff do not include non-exhaust vehicular emissions nor the influence of tailpipe emissions on the formation of secondary aerosols from gaseous precursor emissions (e.g., volatile organic compounds (VOCs), NOx and NH3), whose influence on PM cannot be assessed in near-real time.
A next data treatment step, currently under investigation at some of the CARA supersites, is the combination between PMff and PMwb with ACSM chemical species concentrations to refine the identification of the fine aerosol mode main components within PM10. To do so, ammonium nitrate and ammonium sulfate concentrations are estimated based on ACSM measurements, assuming that molecular stoichiometry between nitrate and ammonium and between sulfate and ammonium is respected. Then, considering that most of the ammonium sulfate, ammonium nitrate and secondary organic aerosols are mainly present in the PM2.5 fraction, it is possible to extrapolate the ACSM data (in PM1) to deduce a concentration of these latter species in the PM10 fraction. This is done by multiplying the ACSM concentrations by the PM2.5/PM1 ratio measured online by a FIDAS instrument (based on optical particle spectrometry, PALAS). The concentrations of ammonium nitrate and ammonium sulfate in PM10 are then assumed to be equal to their loadings estimated in this way. In parallel, secondary organic aerosols (SOA) present in PM2.5 are obtained from the difference between total OA concentrations calculated in PM2.5 and primary organic aerosol contributions assessed to be contained within PMff and PMwb. Using such a methodology, chemically-speciated PM2.5 components as well as coarse mode aerosols (corresponding to the difference between PM10 and PM2.5) can be plotted in a single graph (see section 4.3) in real-time, which again is providing, as a first approach, useful information to AASQAs and decision makers during PM pollution episodes.

Inter-Laboratory Comparisons (ILC) exercises
For QA/QC purposes, the laboratories involved in the CARA program are not only participating in ILCs but also co-organizing some of them at the national and international levels (e.g., starting with EC-OC thermo-optical methods, [57]). During these ILCs, a statistical analysis is generally performed, in accordance with the EN ISO/CEI 17043 standard. For this analysis, the analytical or measurement performances are evaluated using the Z-score calculated as follows: where Zi is the Z score obtained by instrument i, xi is the mean of the measurements performed by instrument i, x* is the robust mean of all the instruments participating in the ILC, and σ^ is the standard deviation of fitness estimated from the entire dataset. It is generally accepted that a Z-score between -2 and +2 is indicative of good performance while a Z-score below -2 or above +2 raises questions about accuracy of the analytical procedures or the proper functioning of the online instrument.
Hereafter, we provide examples of results obtained during three ILCs, focusing on different types of measurements, namely levoglucosan (and isomers) offline analyses, as well as AE33 and ACSM online measurements.

Levoglucosan
Levoglucosan is known as a tracer for biomass burning aerosols for several decades [58], and has been widely used to estimate the concentration of organic matter related to this emission source by applying a multiplicative conversion factor [59]. Its accurate quantification is then critical to avoid significant bias when using such a mono-tracer method. In this context, LCSQA and IGE organized an ILC dedicated to the measurement of this compound and its isomers (mannosan and galactosan), including results obtained by 13 European research laboratories [60]. The main purpose of this ILC was to evaluate the analytical repeatability and reproducibility standard deviations obtained by the participants using their own analytical method. Four test materials were sent to the participants including two ambient air filter samples (corresponding to high winter range concentrations), one laboratory blank filter and one certified reference material (SRM 1649b, urban dust, NIST). Most of participants obtained satisfactory Z-scores (-2 < Z-score < +2). Furthermore, the analytical procedure used by participant appeared to have no significant influence on the quantification of levoglucosan and its isomers, and relatively low expanded uncertainties could then be calculated (e.g., 35 -52 % for levoglucosan). Results from this ILC also provided a robust concentration value of levoglucosan (177.0 µg g -1 with 52 % expanded uncertainty) and its isomers for the SRM 1649b ( Figure 2), consistent with previous values available in the literature (e.g., [61]). Finally, the use of cost-effective amperometry-based methodologies (PAD), which were developing in early 2010's, could be validated for application within the CARA program.

AE33
Deployment of AE33 devices within this observational network from 2013 onwards has also been accompanied with ILC exercises at the national scale. First campaign has been performed at the SIRTA facility in Fall 2014, gathering 19 instruments and 10 AASQAs. Instruments were installed and individually verified by the participants -including leak tests, calibration of sample flow rates and mitigation measures -before being operated in ambient air for a period of 2 weeks (25 Sept.-6 Oct. 2014). Weather conditions were relatively mild during this intercomparison period, with average quarter-hourly temperatures ranging from 5°C to 30°C. A wide range of eBC concentrations ranging from nearly 0 to about 8 µg/m 3 (calculated here as the average of the data provided by all AE33) could be sampled.  Figure 3 presents the Z-scores obtained for each instrument for the following parameters: eBC, eBCff, and eBCwb concentrations, as well as the measured Angström exponent (α). Only two instruments (among the first ones to be produced by the manufacturer) presented a Z-score requiring corrective action (i.e., shipment to the manufacturer for maintenance and upgrade), confirming the good general behavior of the vast majority of AE33 devices used within the CARA program. A very good agreement between other instruments could indeed be observed for each investigated parameter. In particular, and despite the light influence of biomass combustion on eBC levels (less than 10%, on average over the campaign), a satisfactory overall agreement was obtained for the eBCwb concentrations related to this source.
It may be noted that a similar ILC exercise was performed in Spring 2019, leading to similar results as those obtained during this first campaign (unpublished results).

ACSM
Expertise developed on ACSM measurements at the SIRTA facility and then within the CARA program has led to the establishment of the ACMCC, which has been selected as a unit of the ACTRIS topical center for aerosol in situ measurements. ACMCC is notably organizing ILCs within the ACTRIS [34], [62], as well as for ACSM operated by AASQAs. Here, we briefly present results obtained during the last ILC campaign performed at the national scale in June 2019. This campaign included 8 ACSM of the CARA program and has been performed following 3 successive steps: 1) Instruments were installed and operated by AASQAs in ambient air for 3 days.
2) Response factor and relative ion efficiency calibrations as well as general tuning operations were carried out by ACMCC personnel. 3) Instruments were operated in ambient air for another period of 5 days, to check calibration and tuning accuracies as well as to compare instrument performances before and after servicing.  Figure 4 reports, for each measured chemical species, the evolution of the standard deviations obtained from the 8 instruments as a function of the mean concentration, before and after calibrations. First, we note that these standard deviations are exponentially decreasing with concentrations, clearly illustrating higher uncertainties at lower loadings, for each chemical species. Nevertheless, these standard deviations are reaching an asymptotic plateau relatively rapidly towards lowest standard deviation values, pointing to reproducibility measurement uncertainties below 20% for nitrate, sulfate and ammonium loadings typically higher than 0.5 µg/m 3 , and for OA concentrations typically higher than 1.5 µg/m 3 . Similar decreasing features are also observed for chloride, with however much higher dispersion of the data points, due to extremely low concentrations (below 0.1 µg/m 3 ). Finally, these ACSM measurement standard deviations are only slightly improved after calibrations and tunings, illustrating the overall good performance of the ACSM at their respective measurement stations before the ILC exercise.
4. Aerosol chemical composition: overall picture 4.1. Seasonality and spatial distribution of the main PM10 components Figure 5 summarizes the longest (semi-)continuous datasets that could be obtained so far from offline chemical analyses of the main PM10 components within the program. Based on monthly-mean concentrations, it exemplifies long-term variations of these components at four French urban background environments, for which every day (Rouen), every second day (Bordeaux), or every third day (Grenoble and Creil) filter samples were analyzed routinely during the last decade.
In addition, the spatial distribution and the seasonality of some key chemical species are provided for 18 French sites in Appendix B.
Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 OM appears as the predominant component at all sites, with typical PM10 relative contributions of 40-45% on a yearly basis. It is also the main driver of PM temporal variations, displaying a clear seasonality with maxima during the winter season. These observations reflect the diversity and the significance of the various OM origins, including a large zoo of primary (natural and anthropogenic) emission sources but also multiple secondary formation processes in the atmosphere (again from both biogenic and anthropogenic gaseous precursors) [39]. Winter maxima can be partly explained by generally lower planetary boundary layer height at low ambient air temperatures, concentrating all pollutants in a narrow mixing layer. This feature is also related to additional primary sources compared to the rest of the year, especially residential wood burning for heating purposes (section 5). A first insight into the strong influence of this specific source is given by the highest wintertime OM relative contribution in Grenoble compared to the three other sites presented in Figure 5, as well as the highest wintertime OC and levoglucosan absolute concentrations recorded for the different Alpine sites, compared to other French sites (Appendix B). Besides, Bordeaux surprisingly reveals biomass burning influences almost as important as those observed for Alpine sites, e.g., with similar levoglucosan loadings, which has been confirmed by comprehensive studies [42], [54], [63]. Other major contributors to OM include vehicular exhaust emissions (section 6), primary biogenic particles (section 7), and SOA (section 12).
Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 Mineral dust happens to be the second largest contributor to PM10 in Creil, Bordeaux and Grenoble. As for OM, this component is assumed to represent a wide diversity of natural and anthropogenic emission sources, including construction works, non-exhaust vehicular emissions (section 6), transport of soil material from arable lands to urban environments, as well as transport of Saharan dust over the metropolitan territory but also to the French West Indies (section 10) [9], [64]- [66]. Amongst other possible emission sources, one should also mention rare volcanic ash emission events, such as the one which led to the interruption of airplane traffic but also influenced air quality in north-eastern France in Spring 2010 [27]. Overall, maximum mineral dust contributions to PM10 are assessed here from spring to autumn, notably due to lower concentrations of some other PM10 components (starting with OM) but also due to more intense emissions and drier ambient air conditions during this period of the year.
Nitrate can be considered as the main second driver of PM10 seasonal cycles, with highest concentrations and relative contributions in late-winter/early-spring, especially in the northern part of France ( Figure 5 and Appendix B). This is related to the formation mechanisms and the semivolatile nature of ammonium nitrate. Indeed, this secondary species results from the oxidation of nitrogen oxides (NOx) and byproduct combination and/or condensation with ammonia (NH3) in the particulate phase [67]. Its concentration in ambient air is therefore notably depending on the availability of its gaseous precursors and on the strength of (photo-)chemical atmospheric processes. In late-winter/early-spring, intense agricultural activities, including the spreading of fertilizers, under anticyclonic conditions induce a substantial increase of nitrate and ammonium concentrations within the PM2.5 fraction (section 10). Such a phenomenon furthermore relies on the meteorological conditions as ammonium nitrate is submitted to partitioning between the gaseous and the particulate phases as a function of ambient temperature, hygrometry and internal mixing with other compounds [67], which also partly explains the spatial and seasonal nitrate distributions illustrated in this paper.
As nitrate, sulfate originates almost exclusively from secondary formation mechanisms involving gaseous precursors (SO2) emitted by natural sources and anthropogenic activities. Note that primary sulfate contained in sea salts has been removed to total sulfate concentrations presented in Figure 5, therefore illustrating the behavior of non-sea salt sulfate (nss-sulfate) at the four sites exemplified here. Contrary to nitrate, ammonium sulfate is a low volatile species, so that it can stay in the particulate phase at elevated ambient air temperature, i.e., in summer, when enhanced photochemical processes are favoring its formation with subsequent increase of sulfate loadings in PM (section 11).
It should be mentioned that sulfate (as well as nitrate and SOA) can also be present in the coarse aerosol mode (i.e., 2.5-10 µm), notably resulting from the adsorption of gaseous species via acid-base surface reaction on crustal particles or sea salts (whose spatial distribution is discussed further in section 8).
Finally, if discarding concentration peaks observed in December 2016, Figure 5 seems to picture the decreasing long-term trends that have been described elsewhere for regulated PM10 and/or PM2.5 total mass concentrations over France and Europe [7], [68]. In agreement with results obtained from a preliminary analysis performed using the SIRTA online measurement datasets for the 2011-2018 period [46], these PM decreasing features could be mainly driven by the reduction of the ambient concentrations of both nitrate and OM (especially its hydrocarbon-like and its highly oxidized fractions). However, sound statistical analyses require longer datasets (e.g., > 10 years) that the ones discussed here to accurately evaluate long-term trends for each of the main chemically-speciated PM10 component.

Diel variations of the main components of the fine aerosol mode
As mentioned above, high-time resolution AE33 and ACSM measurements benefits the CARA program to continuously document the chemical composition of the fine aerosol mode (PM1/PM2.5). They also allow refining our understanding of the main anthropogenic PM sources, as well as of the influences of atmospheric dynamics and meteorological conditions on urban air quality.  Figure 6 summarizes the yearly-average diel variations of OA, nitrate, sulfate, ammonium and eBC measured by ACSM and AE33 at seven sites of the CARA program over the years 2017-2018. As for the PM10 fraction, OA is the main contributor to the fine aerosol fraction. ACSM data exhibits an evening peak at all sites, more prominent in winter due to both residential heating emissions of wood burning and a lower boundary layer height. This also strongly reflects in the eBC diel profiles ( Figure  7). A smaller peak of OA in the morning at some sites, concomitant to an eBC peak, points out primary emissions due to traffic exhaust. Being representative of background air quality of the Paris region, the SIRTA is less influenced by these local emissions and displays smoother variations. For inorganics, very similar profiles are obtained for both ammonium and nitrate, illustrating the overall predominance of ammonium nitrate over ammonium sulfate in French urban atmospheres. Due to its semi-volatile nature, ammonium nitrate is very much dependent on the meteorological conditions, with a peak in the late-night/early-morning before the daytime temperature and boundary layer heights increase, leading to reduced particulate phase Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 concentrations. Interestingly, the profiles exhibit a decreasing gradient from North to South and East to West, with Lille showing the highest concentrations and Marseille the lowest ones ( Figure 6). As mentioned above, sulfate concentrations in the atmosphere mostly originate from the conversion of SO2, whose sources can be natural (marine emissions) or anthropogenic (industrial, shipping). Far from emission sources, the diel profile reflects rather long-range transport of nonvolatile ammonium sulfate (Paris region, Lille, Creil, Poitiers), whereas more local/regional sources are expected at the sites in Lyon, Bordeaux and Marseille, where enhanced photochemistry in summer also likely plays a role.

Pollution episodes commonly observed within the CARA program
As discussed in section 3.1, highest monthly-mean PM concentrations are quasi-systematically observed during a period covering Winter and early Spring. Logically, this period of the year also corresponds to the predominant occurrences of persistent and large-scale pollution episodes (i.e., with PM10 daily concentrations above 50 µg/m 3 for at least three consecutive days, affecting several administrative regions at the same time). The chemistry of such PM pollution episodes is highly dependent on meteorological conditions and on the type of emission sources which are prevailing while they occur. Figure 8 presents two contrasted examples of such episodes, clearly illustrating their diversity.  Figure 8a pictures a typical springtime PM pollution episode dominated by secondary aerosols and more specifically ammonium nitrate. The time-evolution of this episode, which occurred in March 2014, is extensively detailed elsewhere for the Paris region and using a whole set of in situ and remote sensing measurements at the SIRTA facility [69]. The chemical composition of a similar episode occurring the following year (March 2015) could be further documented using multi-sites ACSM and AE33 measurements [29]. Briefly, this type of episode starts with extremely stable weather conditions associated with very low wind speed, favoring the increase of PM10 in northern France and in the Rhône-Alpes area. Then, anticyclonic conditions associated with (North-)East wind sectors transports polluted continental air masses from neighboring countries, contributing to a sharp deterioration of the situation in France. High loadings of (ammonium) nitrate can be observed at all sites in the northern and eastern parts of the metropolitan territory, which can be related to enhanced ammonia emissions from agricultural activities (manure spreading) in this period of the year. This observation generally applies to all types of stations (traffic, urban and rural background), leading to much more homogeneous pollution distribution than during early-wintertime episodes (Figure 8b). This demonstrates the regional character of this type of PM pollution episode, with large-scale production of secondary aerosols being its main determining factor. Since ammonium nitrate formation needs NH3 emissions but also NOx emissions, the latter ones primarily originating from Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 road transport and other combustion activities (e.g., industries), agricultural activities cannot be considered as the only ones responsible for the occurrence of this type of episode. Such a latewinter/early-spring episode commonly ends with the breakdown of anticyclonic conditions, bringing oceanic air masses from the western wind sector over the North of France whereas the pollution plumes evacuates passing through South-West French regions. This type of persistent and large-scale episode occurred almost every year around the 2010's [70], which, in fact, initially motivated the setting-up of the CARA program to better understand the origins of semi-volatile particles that could be recorded when equipping the conventional monitoring TEOM (tapered element oscillating microbalance) instruments with FDMS (filter dynamics measurement system) devices in AASQAs.
Interestingly, it has been barely observed during the 2016-2020 period, or at much lower amplitude, possibly reflecting the favorable impacts of nitrogen species emission reductions at the European scale [68], [71]. Nevertheless, intense but short-term (1 or 2 day-long) ammonium nitrate-dominated episodes are still recorded in the late-winter/early-spring period (Figure 9) under atmospheric chemistry and dynamic processes that need to be further investigated, but seems to be different from the previous ones. Figure 8b exemplifies a fully different type of PM pollution episodes, which can develop in wintertime under anticyclonic regime associated with stable and cold weather conditions. Such meteorological conditions favor the accumulation of pollutants at a local scale, trapping them in very thin inversion layers. This results in pollution maps presenting a heterogeneous spatial distribution of PM concentration levels, with highest loadings at the vicinity of prevalent emission zones, i.e., urban areas. Since happening during cold periods, such pollution episodes are concomitant with a greater use of domestic heating, especially residential wood burning. During these events, the PM chemical composition is then dominated by carbonaceous matter, with OM-to-PM10 contributions as high as 70-80% in some cases (e.g., Bordeaux in Figure 8b). This impressive burden of (primary) carbonaceous aerosols is further assessed to be mainly influenced by domestic human activities. Figure 10 illustrates monitoring data that could be obtained in near-real time by AASQAs using PM10 and multi-wavelength absorption measurements at different urban background stations during such an early-winter pollution episode which occured in December 2016. Due to the temporality of residential heating emissions, highest PMwb relative contributions, as well as highest PM10 loadings, are generally recorded at nighttime.
Clearly, local-scale emission abatement strategies might help reducing the impact of such wintertime episodes on urban air quality. Associated action plans which have been developed in recent years all over France, e.g., incentive measures to replace open fireplaces (wood stoves) with more efficient (pellet-)burners, should hopefully allow to substantially downscale the magnitude of this type of PM pollution episode in the near-future. Further insights into the overall influence of wood burning emissions on air quality in French urban environments, and the way to evaluate it, is proposed in the following section.

Residential wood burning
Wood burning has been described as a major source of wintertime PM pollution in some of the largest French cities for more than ten years [50], [51], [72], [73]. Yet, it is still a subject of great scientific and political concern, notably due to possible conflicting interest in policy strategies promoting the use of renewable energies.
The high contribution of this primary source to ambient PM loadings is mainly linked with poorly efficient combustion processes, leading to the emission of high amounts of organic species (including PAHs, semi-volatile or intermediate-volatile organic compounds and highly oxidized organic aerosols). Figure 11 summarizes some follow-up results obtained from the comprehensive work presented in [42], which gathers and synthetizes some of the main outcomes of the so-called SOURCES project ( [70], [74]) to analyze pre-existing datasets relevant for urban air quality PM10 source apportionment using PMF in a harmonized way. This figure ranks, in order of importance, the contribution of each identified sources to measured OM concentrations at a yearly-average basis. Although this figure is not focusing on the winter season only, biomass burning still appears clearly as the first source of organic aerosols influencing overall urban air quality in France, with yearlyaverage concentrations of roughly 2.5 ± 1.2 µg/m 3 (representing about 30 ± 10% of total OM loadings). The overwhelming influence of residential wood-burning combustion for heating purposes on this biomass burning fraction -compared to other types of biomass burning -is notably reflected in eBCwb diel variations obtained from AE33 measurements. Indeed, as exemplified for Poitiers in Figure 7, these patterns display well-marked diel cycles with eBCwb concentration maxima at late evening, just following eBCff (traffic emissions) peaks, i.e., when individuals are coming back home. A week-end effect might even be observed at stations located in residential suburban areas, with highest mean concentration peaks observed for Sunday night in winter [28], [46]. Nevertheless, some events of open green waste burning may also temporarily influence urban air quality, especially in fall [75]. Right: average relative contribution (%). Adapted from work described in Weber et al. [42] As mentioned above, levoglucosan is commonly used as a specific tracer of biomass burning aerosols. This organic compound appears to represent a fairly stable fraction (i.e., 8 ± 2 %) of the PM10 biomass burning factor obtained from harmonized PMF analyses described in [42]. These results can be used to check the consistency of PMwb estimates obtained from AE33 measurements. This is exemplified here in Figure 12  Despite rather homogeneous contributions of levoglucosan (and OC) within the biomass burning factor from one urban site to another, other species (such as EC, potassium, chloride, and nitrate) display higher variabilities [42], which can be partially explained by fluctuations in the type of wood which is used and/or in burning conditions. Furthermore, biomass burning emissions in ambient air may favor i) nighttime mechanisms leading to the condensation of ammonium nitrate in the particulate phase [76], ii) the build-up of SOA onto wet aerosols [77], and iii) the formation of organo-nitrate (or nitrated organic compounds) and/or organo-sulfate [41], [78], probably increasing the impacts of wood burning emissions onto wintertime urban air quality up to levels of complexity and of magnitude which are still to be better understood and quantified. As some of the secondary compounds formed through these processes are absorbing light at near-ultraviolet wavelengths [77], [79], it may be assumed that they are partially detected and somehow accounted along with primary biomass burning aerosols when using the Aethalometer model [51], [79]. For ACSM measurements, the diversity of emitted species and secondary processes potentially involving them, makes it difficult to acknowledge that they can always be fully captured in a single (BBOA) factor obtained from classical PMF analysis applied to organic mass spectra [48], [80]- [82]. In this context, results potentially obtained for this BBOA fraction from a classical and automated ACSM organic mass spectra PMF data treatment procedure within an operational network such as the CARA program may happen to be misleading for the general public and decision makers. Further investigations and methodological developments, e.g., using constraints from independent and complementary methodologies, are still to be tested before implementation in a routine way (see also section 12).

Road transport (primary emissions)
On annual average and overall, road transport isthe second contributor to OM loadings affecting air quality in French urban environments ( Figure 11). Related organic aerosols are assessed to originate from tailpipe emissions (along with roughly equivalent EC loadings), but also from nonexhaust emissions, including resuspension of various types of particles which may deposit on the road [66], [83]. Substantial amounts of metallic species are also expected to be related to various traffic-related abrasion phenomena (from road surfaces, tires, brake pads, etc.). In practice, the discrimination between vehicular exhaust and non-exhaust particles appears to be hardly possible when applying classical PMF approach on chemical speciation data obtained at urban background conditions. Similarly, gasoline and diesel tailpipe emissions are barely separated by receptor models at background sites. The mixing of these different emission types in a single PMF factor is leading to rather heterogeneous primary traffic chemical profiles obtained from a given site to another [42].
However, source apportionment studies conducted at roadside monitoring stations might allow to delineate them between two (or more) specific fingerprints, as exemplified hereafter with results recently obtained for a traffic-oriented station (so-called Renéville) of Fort-de-France, as part of a recent study conducted in Martinique Island (see section 9 and [84]).
Briefly, 89 PM10 daily filter samples collected by Madininair AASQA at Renéville over 2018 were analyzed offline by the LCSQA. A relatively small set of chemical variables (31) were selected and used within the input data matrix to fulfill equilibrium requirements between the number of samples and the number of chemical species [44]. Six factors could be obtained from this PMF analysis. Figure  13 presents the relative contributions of each of these six factors to PM10 concentrations as well as to each of the individual chemical species considered here. Besides Saharan dust, sea salts, biomass burning and industrial activities, two factors are assumed to be related to road transport. The first one -'Road transport 1' -exhibits a relatively low EC/OC ratio (0.5) and includes almost all antimony and tin measured at this station, whereas the second one -'Road transport 2' -displays a much higher EC/OC ratio (about 2.6) and a high chromium contribution. These profiles suggest a possible differentiation of i) light-duty vehicles equipped with gasoline engine and metalloid-lined brakes, and ii) other types of (heavy-duty) vehicles equipped with diesel engine and/or chrome plated cast iron brakes [85]. Alternatively, it might also be hypothesized that the first factor mainly corresponds to non-exhaust vehicular emissions while the second one is primarily composed of exhaust emissions (with chromium, but also other elements such as molybdenum or copper, possibly originating from catalyst converters). Further works, including detailed investigations of the elemental and isotopic composition of non-exhaust road transport emissions [86], might greatly help to better evaluate the impact of these sources on urban air quality. Overall, primary particles related to road transport emissions are assessed to represent about 15% of annual mean PM10 concentrations recorded in French cities [42]. It is reminded here that gaseous tailpipe emissions contribute also substantially to ammonium nitrate (section 10) and SOA (section 12) formation in urban environments. Furthermore, road transport is a significant source of ultrafine particles (UFP) whose monitoring in French urban environment will be developed at the operational level in coming years [87].

Primary biogenic organic aerosols (PBOA)
Besides anthropogenic activities, vegetation may also be a substantial contributor to primary OM emissions [88], [89]. The large number of CARA stations, combined with many other series of samples from specific complementary programs, has made possible to produce a national synthesis on the importance of such emissions in the PM10 fraction. Based on the use of polyols as PBOA markers [22], [90], it has been shown that primary OM emissions associated with fungal species might represent up to about 40% of the total summertime-averaged OM concentration in PM10 at some urban background stations ( Figure 14). With much lower influence at some other stations, an overall summertime contribution of about 25% can be considered on average at the national level during this season. On a yearly basis, OM comprised in fungal spores are assessed to represent an annual mean concentration of about 1 µg/m 3 in French urban background environments. Size-segregated offline chemical speciation measurements also indicate that these particles are essentially present in the coarse fraction of PM10, with relatively similar temporality and high spatial homogeneity of concentrations at sub-regional scale (at least 100 km). Furthermore, these chemical species have been demonstrated to be strongly associated with the abundance of specific fungal genera, which are distinct in different climates and are carried by vegetation rather than from the soil. This is particularly true in rural areas of field crops, where agricultural activities are highly emitting (especially during harvest) [26], [91].
It should also be underlined that fungal spores are not the only type of PBOA. In particular, it is likely that another substantial part of the bioaerosol loadings is made up of plant debris that can be measured directly via the free cellulose in PM10. Some of the first published annual cycle measurements show that concentrations of the order of several hundred ng/m 3 may be present, again in the coarse fraction of PM10, throughout the year [20]. The spatial variability of these concentrations at the urban scale is quite high, indicating nearby sources for this fraction.

Sea salts
At the global scale, sea salts represent the second largest contributor to the total aerosol burden in the atmosphere (after mineral dust), playing a key role in the overall water cycle and possibly moderating climate change [6]. Logically, they can then make up another large fraction of natural aerosols observed in urban environments, mostly in the coarse PM mode. When located close to the shore, monitoring stations might even be over-influenced by these particles, which can also cause practical measurement issues related to their high hygroscopicity and corrosivity [92].
Sea salts are primarily composed of Cl -, Na + , SO4 2-, Mg 2+ , and Ca 2+ . However, they can rapidly age during transport, notably due to the mixing with anthropogenic emissions over lands or before reaching them. This phenomenon commonly induces chloride replacement by sulfate and nitrate via acid-base reactions. The mixing between sea spray and heavy fuel oil combustion residues -e.g., emitted from shipping or industrial activities -also leads to the build-up of an aged sea-salt component detected in source apportionment studies through the presence of a few chemical species generally attributed to anthropogenic emissions, such as EC and Ni (see Figure 13) [36].
As illustrated in Appendix B, seasonal-average Na + concentrations up to 1 µg/m 3 are observed along all coastlines of the French metropolitan territory (i.e., the Atlantic Ocean, the English Channel, and the Mediterranean Sea), with large influences inland up to distances of several hundreds of km. Lowest concentrations are detected in the most eastern part of France and in Alpine valleys. A wellmarked seasonal cycle, with winter maxima, is observed at all sites located in northern France, whereas a reverse pattern tends to be seen for southern regions. Finally, depending on air mass origins (and on the strength of the weather front in the westerlies over the Atlantic), high loading events can be observed, where daily average fresh + aged sea salt concentrations typically of 10-15 µg/m 3 are not uncommon up to a couple of hundreds km from the western coasts [36], [42].

Saharan dust
The French West Indies are regularly impacted by Saharan dust events from West Africa (Sahel) or North Africa (Sahara) between March and October [93]. These desert dust particles are transported by the trade winds across the Atlantic Ocean, reaching the Caribbean Sea in one to two weeks [94], resulting in ambient particles concentrations above the PM10 EU air quality standards. This is the case, for instance, in Martinique Island where the 50 µg/m 3 daily mean concentration threshold has been exceeded every year since 2005, except for 2008. In addition to their global climate and ecosystem impact [95], they have been associated with adverse health effects, including cardiorespiratory mortality and morbidity [96]. Recently, Viel et al. [97] also evidenced an increasing risk of preterm births with the levels of PM10 and the frequency of intense Saharan dust events over Guadeloupe Island.
In this context, filter samples collected in the French West Indies have been occasionally analyzed since 2010 in the frame of the CARA program, systematically indicating a predominant influence of Saharan dust during exceedances of the PM10 daily threshold value [98], [99]. A more comprehensive study has been conducted for the year 2018, with daily PM10 filter sampled semicontinuously at 5 monitoring stations in Martinique Island [84]. More than 300 filters were then selected for extensive offline chemical characterization, during and out of PM10 daily threshold exceedances. Outcomes of this study notably indicate a homogeneous spatial distribution of the Saharan dust loadings during PM10 concentration peaks and confirm the major role played by natural particles (Saharan dust but also relatively high sea salt background levels) onto air quality over the island. Figure 15 summarized some the results obtained from the PMF analysis conducted for the Renéville roadside monitoring station located in Fort-de-France (see section 6). On average for daily PM10 concentrations above 50 µg/m 3 , Saharan dust represents about half the PM10. It should be noted that this relative contribution raises up to more than 70% during the highest PM10 daily peaks (above 100 µg/m 3 ). Here, measures taken towards the abatement of anthropogenic emissions will not lead to significant reduction of the number of days exceeding the PM10 threshold value. Nevertheless, this may help lowering adverse health effects induced by the overall ambient aerosol burden. Finally, it should be noticed that Saharan dust outbreaks also influence the French metropolitan territory, especially its southeastern part (e.g., [64], [100]), although less impressively than in the French West Indies and more often evidenced in the southern European countries [8], [9].

Ammonium nitrate
In order to delineate the sources of nitrogen species during springtime PM pollution episodes dominated by ammonium nitrate (see section 4.3), we tried to perform isotopic measurements of nitrogen, conducted separately on nitrate and ammonium [101], [102]. Indeed, it is known that N produced during biological processes is depleted in heavy isotopes compared to atmospheric N [103], [104]. This may be a way to identify the source of N coming from enhanced ammonia emissions in late-winter/early-spring due to agricultural manure spreading [105]. This has been notably performed within the INACS research project, notably taking advantage of yearly semi-continuous series of filter samples collected in the frame of the CARA program [106]. Obtained results are exemplified here with the case of Lyon, including a persistent large-scale pollution episode in spring 2012 ( Figure 16). It clearly shows a strong depletion in δ 15 N of ammonium concomitant to a very large increase of its concentrations during this overall episode, underlining the strong impact of the agricultural source with biological ammonia fertilizers during this type of event. More surprisingly, relatively low ammonium concentrations recorded during the rest of the year reveal to be mainly influenced by other types of ammonia sources (with much higher δ 15 N values). This might be related to the importance of still poorly documented specific urban emitters, such as Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 vehicular exhaust and garbage piling, and subsequent rapid formation of ammonium nitrate at the local scale [67], [107]. Nevertheless, long-range transport of ammonium nitrate from the continental sector has often been evidenced, notably using backward air mass trajectory analyses or CTMs [29], [108]. To illustrate such a phenomenon, a comprehensive trajectory analysis has been carried out here using ACSM datasets available over the years covering the 2016-2017 period (corresponding to data collected at SIRTA, Creil, Metz, Poitiers, Bordeaux, Lyon and Marseille stations). The analysis consisted in the combination of single site concentration-weighted trajectory analyses using the ZeFir software package [109], as already performed and explained in [110] to investigate particulate sulfate formation from SO2 volcanic emissions. Briefly, the HYSPLIT model [111] was used to calculate 5-day back-trajectories every 6h, from 1°x1° Global Data Assimilation System (GDAS) files. Trajectory endpoints were discarded if significant precipitation occurred (> 1 mm/h). Finally, a sigmoidal weighing function was applied on the multi-site results in order to downweight grid cells with low trajectory density. For nitrate, this analysis reveals hotspots mainly located in central Europe, northern Italy and the Benelux regions (Figure 17a), in good agreement with previous studies focusing on shorter time periods and/or for a given receptor urban area (e.g., [29], [36], [41]).

Ammonium sulfate
The same analysis has been performed for sulfate concentrations, also known to be highly influenced by long-range advection from continental Europe, at least for northern France [28], [36]. As for nitrate, central Europe and the Benelux regions, corresponding to abundant SO2 emissions [112], are also clearly captured (Figure 17b). Additional regions are also highlighted in the Mediterranean and -less importantly -the Channel areas, suggesting a substantial role of shipping emissions. Marseille represents here a peculiar case with a clear hotspot, probably also related to intense industrial and petrochemical activities [113]. Indeed, Marseille is directly downwind of the large Fos-Berre petrochemical area (located 40 km away from city center, at its north-western side) and the city harbor basins. Depending on the wind direction, the urban background site is frequently impacted by industrial and shipping traffic plumes [114]. A typical summertime case study is provided below to exemplify this phenomenon. Figure 18 presents the temporal variation of the wind direction, of key chemical species concentrations (OA, SO4 2-, NH4 + , NO3 -, Cl -, eBC, V, Ni and Pb in the particulate phase, as well as SO2) and UFP (considered here in the 3-158 nm size range) concentration and size distribution. In the early morning (07:00 UTC), very intense plumes of SO2 associated with UFP, OA and metals are observed when sea breeze sets in (190°-240°), advecting air masses impacted by industrial and shipping emissions over the city. Several smallest plumes advection occurs through the day under this sea breeze system. These events are followed by enhanced ambient sulfate concentrations in the city center, assessed to originate from SO2 photochemical oxidation. Finally, in the evening (19:00/20:00 UTC), the land breeze brings back the processed air masses to the city, with another increase of sulfate loadings.
This example shows a large increase on short time scale of ammonium sulfate concentrations, in relation with emissions from heavy fuel oil combustion emissions (which can be traced by vanadium). Elevated concentrations of both sulfate and vanadium can also be seen on large temporal (yearly) and spatial (Mediterranean coastline from Nice to Port-de-Bouc) on the maps in Appendix B, with concentrations larger than anywhere else in France. Constraining the contribution of refineries and shipping emissions remains a tremendous challenge in this whole region.

Secondary organic aerosols (SOA)
Secondary organic aerosols are known to make up a very significant amount of OM, and then also of PM, especially in the fine aerosol mode. [39]. Yet, they are barely represented within the outputs of classical filter-based source apportionment studies. This is reflected in Figure 11 (for instance), where only one PMF factor (namely, marine SOA) is directly linked to its precursor and secondary origins (using MSA, which is already not so commonly used in such type of analysis). To bridge the gap between real-world atmospheric chemistry and receptor modelling results, increasing efforts have been made to determine molecular organic compounds, that can be identified as typical oxidation by-products of biogenic or anthropogenic SOA precursors and then used as specific markers in source apportionment studies [20], [115], [116].
Recently, some of these SOA markers could be monitored offline at a couple of sites within the CARA program (i.e., Grenoble and SIRTA). Their use in PMF analysis allowed the deconvolution and apportionment of different SOA sources not commonly resolved in the literature, highlighting substantial contributions of specific biogenic and anthropogenic SOA factors [18], [20], [21], [117]. More precisely, a peculiar secondary anthropogenic fraction could be evidenced in Grenoble, possibly related with Fenton-like chemistry and self-amplification cycle of SOA formation during Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 wintertime PM pollution episodes [18]. At SIRTA, various anthropogenic SOA fractions could be linked to different secondary oxidation processes involving the following precursors: toluene and naphthalene, phenolic compounds, and PAHs (notably leading to the formation of DHOPA and phthalic acid, methylnitrocatechols and oxy/nitro-PAHs, respectively). These precursors are assumed to be mainly emitted by biomass burning and vehicular exhaust, and the sum of corresponding SOA loadings might represent up to 40% of daily PM10 concentrations in Greater Paris during the winter season [21], [117]. Regarding biogenic SOA, and besides the MSA-rich factor, one additional fraction could be attributed to the oxidation of alkenes (i.e., isoprene and terpenes) emitted by terrestrial vegetation at both sites, more importantly at summertime. Nevertheless, these PMF analyses did not allow to discriminate between both types of alkene byproducts, and then between VOC emissions from deciduous and resinous trees. This could be done for the SIRTA dataset using another type of source apportionment method, i.e. the so-called SOA-tracer method [115], which, in return, still needs fine optimizations to better account for anthropogenic SOA [117], [118]. Further efforts are then required to refine the source apportionment of the various SOA fractions, along with other PM contributors, using receptor models. It should however be noted that the search for proper biogenic tracers implies very sensitive and labor-intensive (and then costly) offline chemical analyses [116], that may not be applied easily and routinely within operational networks. Figure 19: Comparison of PMF outputs obtained using a conventional approach applied to ACSM OA mass spectra and a novel methodology combining these mass spectra with key molecular markers in the PMF input data matrix. Adapted from Srivastava et al. [82] On the other hand, classical source apportionment methods applied to AMS (or ACSM) organic mass spectra commonly reveals overwhelming contributions of different types of oxygenated organic aerosols (OOA), that are usually attributed to SOA [119], [120]. While allowing for a specific discrimination of the full submicron OA burden, results obtained using such methodologies rarely consider interactions with inorganics and do not trivially link the different OOA fractions to specific SOA precursors and/or formation mechanisms. Advanced data treatment strategies combining OA mass spectra with datasets obtained from complementary co-located measurements could provide further insight into these PM origins (e.g., [120], [121]). In particular, adding specific SOA markers to ACSM organic mass spectra proved to allow a better understanding of OA sources during a PM pollution events documented at SIRTA in March 2015. As detailed in [82], this novel methodology lead to the deconvolution of ten different PMF factors, among which two could be clearly attributed to biogenic SOA and four could be attributed to relatively fresh anthropogenic SOA fractions mainly linked to traffic and/or biomass burning emissions at the urban scale (as revealed by diel variations of these factors). These factors are gathered in Figure 19 according to the main origin of their precursors (i.e., biogenic vs. anthropogenic SOA), for a simpler comparison with the classical ACSM approach. Since it requires the offline chemical analysis of a limited number of key molecular markers only, such a source apportionment strategy might be applied regularly for the CARA program supersites to refine our knowledge on the origins of both the primary and secondary aerosol fractions.

Conclusions and perspectives
This paper presents some of the main activities and results of the French CARA program, which has been developed for more than ten years. Being operated in the field by regional air quality monitoring networks in close collaboration with the LCSQA (and other academic partners) for so long, this program is unique in Europe. It is designed on simple but robust methodologies to document the PM chemical composition and sources in urban environment. It also includes -or facilitates -state-of-the-art research activities, which, in return, allow stepwise upgrades and optimization of conventional ambient air monitoring systems in France. In this regard, the CARA program may be considered as a pioneer operational tool for the implementation of European policydriven air quality strategies.
This program is notably based on the application of sound QA/QC procedures, including ILC exercises for key measurement techniques. Initially developed for offline chemical analyses of daily PM10 samples collected quasi-continuously or during specific campaigns, it now also benefits from dense networks of AE33 and ACSM monitoring sites. These highly time-resolved and online measurements allow to extensively document the major chemical species that compose the fine aerosol mode. Moreover, they greatly help for a better understanding, in near-real time, of the main anthropogenic emission sources, of the potential influences of secondary processes and/or of geographical origins of PM pollution episodes. The fine aerosol mode being mainly influenced by anthropogenic pollutant emissions, this online monitoring set-up provide direct information to AASQAs (and to the LCSQA), which might then be considered by decision makers when it comes to implement and prioritize selected short-term legal decisions and long-term action plans.
Scientific results presented in this paper notably illustrate the predominant influences of i) wintertime residential wood burning, and ii) year-long (exhaust and non-exhaust) vehicular emissions, on the overall urban air quality in France. They also exemplify the importance of additional sources at local scale (such as shipping and industrial emissions in Marseille) or during critical periods of the year (e.g., ammonia emission by agricultural activities under stable lowpressure meteorological systems in late-winter/early-spring). Moreover, recent works on natural particles evidenced the substantial contribution of PBOA in the coarse aerosol mode, especially in summer, as well as the overwhelming impacts of Saharan dust outbreaks on the PM10 limit value exceedances in the French West Indies. Furthermore, some efforts have been dedicated to the refinement of our knowledge on the origins of secondary species, notably using new organic molecular markers in receptor models.
Development perspectives also comprise near-real time source apportionment methodologies, along with a wider characterization of the PM chemical and physical properties, as well as of inorganic gaseous precursors playing a prominent role in the formation of secondary aerosols. UFP (or, at least, total particle number concentration) and ammonia measurement systems are currently deployed by the AASQAs in a growing number of urban sites. Considering the oxidative potential, and probable associated health outcomes, of non-exhaust PM vehicular emissions [3], high timeresolved measurements of a wide range of heavy metals [122] should be envisaged in near-future. This would also greatly help for a better source apportionment of other types of mineral particles, which substantially contribute to the PM10 burden in urban areas.
These near-future milestones shall be elaborated and reinforced in synergetic approaches with European programs and research infrastructures. In particular, the ACTRIS has a vocation to provide guidance for the monitoring of novel air quality parameters in urban environments, including advanced in situ observations but also key remote sensing variables describing atmospheric dynamics. Closer interactions with the EMEP (and its French diversification, i.e., MERA) are also underway. In this regard, the transition of some offline monitoring site towards selected remote stations (which could also progressively host online instrumentation) has been decided in 2020. This Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 January 2021 doi:10.20944/preprints202101.0182.v1 will allow for a better characterization of (sub)regional background air quality and to set-up so-called 'paired sites' strategies to refine the determination of long-range transport vs. local air pollution [10]. Finally, extended datasets and filter collections that have been gathered in the frame of the CARA program are already made available for complementary studies, such as comparison exercises between in situ observations and CTM outputs or remote sensing measurements (e.g., [105]), as well as cross-over studies on adverse health effects [3], [25]. Such initiatives shall be further developed in a broader landscape (e.g., including near-real time data provision to the European Copernicus Programme) while continuous online measurements at multi-sites will also set the basis for achieving long-term investigations, like epidemiological studies and air quality trend analyses. Funding: This work has been mainly funded by the French Ministry of Environment via its annual financial supports to the LCSQA, to the AASQAs, and to INERIS research programs. It has also been partially funded by ADEME (notably for the SOURCES project, under grant agreement 1462C0064, and for the DECOMBIO (1362C0028), PM-DRIVE (1162C0002), CAMERA (1062c0008), INACS (1262c0011), and QAMECS (1262c0011) programs), as well as by targeted initiatives led by the AASQAs, and by ANDRA. Offline chemical analyses are also supported at IGE by the Air-O-Sol platform with Labex OSUG@2020 (ANR10 LABX56) fundings and at IMT LD by the Labex CaPPA (ANR-11-LABX-0005-01) and CPER CLIMIBIO projects (also supporting ACSM and AE33 measurements at Lille supersite). Measurements at SIRTA have been partially supported by CNRS, CEA and ACTRIS-France (notably through the 'CLAP national observation service'). They were also part of the EUFP7 and H2020 ACTRIS projects (grant agreements Nos. 262254 and 654109). Finally, the PhD of Samuël Weber was funded by ENS Paris.
Acknowledgments: Authors acknowledge the precious works in the field and/or in the lab. over the last decade of so many technicians, engineers and students from AASQAs, INERIS, IMT LD, IGE, and LSCE, who we prefer not to try and mention individually here, fearing we might unfortunately miss some of them. We strongly hope they may all recognize themselves and accept our warm acknowledgements. Greetings and personal supports from various individuals at funding and leading institutes are also highly appreciated.

Conflicts of interest:
The authors declare no conflict of interest.