Arabitol, mannitol and glucose as tracers of primary biogenic organic aerosol: influence of environmental factors on ambient air concentrations and spatial distribution over France

Abdoulaye Samaké, Jean-Luc Jaffrezo, Olivier Favez, Samuël Weber, Véronique Jacob, Trishalee Canete, Alexandre Albinet, Aurélie Charron, Véronique Riffault, Esperanza Perdrix, Antoine Waked, Benjamin Golly, Dalia Salameh, Florie Chevrier, Diogo Miguel Oliveira, Jean-Luc Besombes, Jean M.F. Martins, Nicolas Bonnaire, Sébastien Conil, Géraldine Guillaud, Boualem Mesbah, Benoit Rocq, Pierre-Yves Robic, Agnès Hulin, Sébastien Le Meur, Maxence Descheemaecker, Eve Chretien, Nicolas Marchand, and Gaëlle Uzu.

It is commonly acknowledged that SCs (particularly arabitol and mannitol) originate from primary biogenically derived sources such as bacterial, fungal spores, and plant materials Golly et al., 2018;Gosselin et al., 2016;Holden et al., 2011;Kang et al., 2018;Medeiros et al., 2006b;Wan et al., 2019;Yan et al., 2019;Yttri et al., 2007;Zhu et al., 2018a). Some studies have characterized the composition of SC in topsoil samples (for fractions larger than PM 10 ) from both natural (i.e., uncultivated) and agri-  (Medeiros et al., 2006a;Rogge et al., 2007;Simoneit et al., 2004;Wan and Yu, 2007). The authors suggested that the particulate arabitol, mannitol, and glucose are introduced into the atmosphere mainly through resuspended soils or dust particles and associated biota derived from natural soil erosion, unpaved road dust or agricultural practices. Conversely, Jia and Fraser (2011) reported higher concentrations of SC relative to PBOA in size-segregated aerosol samples collected at a suburban site (Higley, USA) compared to the local size-fractionated soils (equivalent to atmospheric PM 2.5 and PM 10 ). This suggested that direct emissions from biota (microbiota, vascular plant materials) could also be a significant atmospheric input process for SC at this suburban site.
A large database on SC concentrations was obtained over France in the last decade. It already allowed the investigation of the size distribution and seasonal variabilities of SC concentrations in aerosols at 28 French sites, notably showing that SCs are ubiquitous primary aerosols, accounting for a significant proportion of PM 10 organic matter (OM) mass (Samaké et al., 2019). Results confirmed that their ambient concentrations display a well-marked seasonality, with maximum concentrations from late spring to early autumn, followed by an abrupt decrease in late autumns and a minimum concentration during wintertime in France. This study also showed that the mean PBOA chemical profile is largely dominated by organic compounds, with only a minor contribution of dust particle fraction. The latter result indicated that ambient polyols could most likely be associated with direct biological particle emissions (active spore discharge, microbiota released from phylloplane or phyllosphere, etc.) rather than with the microorganism-containing soil resuspension. These observations call for more investigations of the predominant SC (and PBOA) emission sources.
Cellulose, a linear polymer composed of D-glucopyranose units linked by β-1,4 bonds, is the most frequent polysaccharide occurring in terrestrial environments (Ramoni and Seiboth, 2016). Plant materials contain cellulose, which has been reported as a suitable proxy to evaluate the vegetative debris contribution to OM mass (Bozzetti et al., 2016;Daellenbach et al., 2017;Glasius et al., 2018;Hiranuma et al., 2019;Puxbaum and Tenze-Kunit, 2003;Sánchez-Ochoa et al., 2007;Yttri et al., 2011b). The ambient PM 10 cellulose has been shown to be abundant in the European semirural or background environments (accounting for 2 % to 10 % of OM mass) (Glasius et al., 2018;Sánchez-Ochoa et al., 2007) and Nordic rural environments in Norway (contributing to 12 % to 18 % of total carbon mass) (Yttri et al., 2011b). Thus, simultaneous concentration measurements of cellulose and SC can provide essential information into their emission source dynamics.
As the continuation of our previous work (Samaké et al., 2019), the present paper aims to delineate the processes that drive the atmospheric concentrations of SC and then PBOA. This is achieved through (i) the analysis of simultaneous an-nual short-term time series of particulate SC concentrations over pairs of sites across multiple space ranges, including local, regional, and nationwide sites, and (ii) the investigation of links between concentrations and series key parameters such as meteorological and phenological ones. Simultaneous annual short-term concentration measurements of SC and cellulose were performed to better understand their sources correlations.
2 Material and methods

Sampling sites
Daily PM 10 concentrations reported in the present work were obtained from different research and monitoring programs conducted over the last 6 years in France. Within the framework of the present study, we carefully selected sites sharing at least 1 complete year of concurrent monitoring with another one, to be representative of the annual variation cycles. The final dataset includes data from 16 sites, which are distributed in different regions of France ( Fig. 1) and cover several main types of environmental conditions in terms of site topography, local vegetation, and climate. The characteristics and data available at each sampling site are listed in Table S1 of the Supplement, together with the information on the annual average concentrations of aerosol chemical composition (Table S2). Detailed information on the sampling conditions can be found in Samaké et al. (2019), such as the campaign periods, number of collected PM samples, sampling flow rates, sample storage and handling, etc. Note that, the previous database (Samaké et al., 2019) has been updated here with arabitol and mannitol in PM 10 collected at the suburban site of Nogent-sur-Oise for a series covering the years 2013 to 2017.

Chemical analyses
Daily (24 h) PM 10 samples were collected onto prebaked quartz fiber filters (Tissuquartz PALL QAT-UP 2500 150 mm diameter) every third or sixth day but not concurrently at all sites. They were then analyzed for various chemical species using subsampled fractions of the collection filters and a large array of analytical methods. Details of all the chemical analysis procedures are reported elsewhere (Golly et al., 2018;Samaké et al., 2019;Waked et al., 2014;Weber et al., 2018). Briefly, primary sugar compounds were extracted from filter aliquots (punches typically about 10 cm 2 ) into ultrapure water. The extracts are then filtered using a 0.22 µm Acrodisc filter. Depending on the site, analyses were conducted either by the IGE (Institut des Géosciences de l'Environnement) or by the LSCE (Laboratoire des Sciences du Climat et de l'Environnement) (Samaké et al., 2019). At the IGE, extraction was performed during 20 min in a vortex shaker, and analyses were achieved using highperformance liquid chromatography with pulsed amperomet- ric detection (HPLC-PAD). A first set of equipment was used until March 2016, consisting of a Dionex DX500 equipped with a three-column Metrosep (Carb 1-Guard + A Supp 15-150 + Carb 1-150); the analytical program was isocratic with 70 mM sodium hydroxide (NaOH) as an eluent for 11 min, followed by a gradient cleaning step with a 120 mM NaOH as an eluent for 9 min. This procedure allows the analysis of arabitol, mannitol, and glucose (Waked et al., 2014). A second set of equipment was used after March 2016, with a Thermo-Fisher ICS 5000+ HPLC equipped with a 4 mm diameter Metrosep Carb 2 × 150 mm column and a 50 mm pre-column. The analytical run was isocratic with 15 % of an eluent of sodium hydroxide (200 mM) and sodium acetate (4 mM) and 85 % water, at 1 mL min −1 . At the LSCE, extraction was performed for 45 min by sonication and analyses were achieved using an ion chromatography instrument (IC, DX600, Dionex) with pulsed amperometric detection (ICS3000, Thermo-Fisher). In addition, a CarboPAC MA1 column was used (4 × 250 mm, Dionex) along with an isocratic analytical run with 480 mM sodium hydroxide eluent. This analytical technique allows us to quantify arabitol, mannitol, and glucose (Srivastava et al., 2018). Examples of standard solution and sample raw HPLC-PAD chromatograms are presented in Fig. S1.
For cellulose quantification, we used an optimized protocol based on that described by Kunit and Puxbaum (1996) and Puxbaum and Tenze-Kunit (2003), in which the cellulose contained in the lignocellulosic material is enzymatically hydrolyzed into glucose units before analysis. Since the alkaline peroxide pretreatment step used to remove lignin in the original protocol results in a loss of sample material, it has been avoided in this study. Therefore, only the "free cellulose" is reported in our samples. Note that Sánchez-Ochoa et al. (2007) consider that this free cellulose could represent only about 70 % of the total cellulose in air samples and that the total cellulose could represent only about 50 % of the "plant debris" content of atmospheric PM. Very few other results are available on this topic (Bozzetti et al., 2016;Glasius et al., 2018;Vlachou et al., 2018;Yttri et al., 2011b). The protocol has been improved to increase sensitivity and accuracy by reducing the contribution of glucose in the blanks and by using an HPLC-PAD as the analytical method for the determination of glucose concentrations. Trichoderma reesei cellulase (> 700 U g −1 , Sigma Aldrich) and Aspergillus niger glucosidase (> 750 U g −1 , Sigma Aldrich) were used as saccharification enzymes. The protocol is detailed in Sect. S2 in the Supplement.
Field blank filters (about 10 % of samples) were handled as real samples for quality assurance. The present data have been corrected from field blanks. The reproducibility of the analysis of primary sugar compounds (polyols, glucose) and cellulose, estimated from the analysis of sample extracts from 10 punches of the same filters, were in the range of 10 %-15 %. The quantification limits of primary sugar compounds and cellulose ranged from 0.63 to 0.89 ng m −3 . About 2800 samples are considered in this work for the polyols and glucose series, while 290 samples (from the sites of Grenoble_LF and OPE-ANDRA) are considered for the cel-Atmos. Chem. Phys., 19, 11013-11030, 2019 www.atmos-chem-phys.net/19/11013/2019/ lulose series. Hereafter, the term "polyols" is used to refer uniquely to the sum of arabitol and mannitol concentrations.

Meteorological data and leaf area index (LAI) measurements
Ambient weather data were not available at all monitoring sites (see Table S1). In this study, data including daily relative humidity (%), nighttime temperature ( • C), average and maximum temperatures ( • C), wind speed (m s −1 ), solar radiation (W m −2 ), and rainfall level (mm) for the sites of Marnaz and OPE-ANDRA ( Fig. 1), representing different climatic regions and environmental conditions, were obtained from the French meteorological data sharing service system (Météo-France) and ANDRA (French national radioprotective agency, in charge of the OPE-ANDRA site). The LAI, which is defined as the projected area of leaves over a unit of land, is an important measure of the local vegetation density variation (Heald and Spracklen, 2009;Yan et al., 2016a, b). For this study, we used the MODIS Collection 6 LAI product because it is considered to have the highest quality among all the MODIS LAI products (Yan et al., 2016a, b). The MCD15A3H product uses both Terra and Aqua reflectance observations as inputs to estimate daily LAI at 500 m spatial resolution, and a 4 d composite is calculated to reduce the noise from abiotic factors. Using a 2 × 2 km grid box around the monitoring site, the local vegetation density variation was retrieved from LP DAAC (https: //lpdaac.usgs.gov/, last access: 15 March 2019) for the sites of Marnaz, OPE-ANDRA, and Grenoble_LF.

Data analyses
All the statistical analyses were carried out using the opensource R software (R studio interface, version 3.4.1). Several statistical analyses were performed on the concentrations to identify the spatial patterns of emission sources and the potential parameters of influence as explained below.
The normalized cross-correlation (NCC) test was chosen to examine the potential similarities among the monitoring sites for particulate SC concentrations, in terms of short-term temporal trends (e.g., synchronized periods of increase or decrease, simultaneous fluctuations during specific episodes). The main advantage of NCC over the traditional correlation tests is that it is less sensitive to linear changes in the amplitudes of the two time series compared. Therefore, to reduce the possibility of spurious "anti-correlation" due to highly variable concentration ranges, data were amplitudenormalized prior to correlation analysis. A thorough discussion of the normalized cross-correlation method can be found elsewhere (Bardal and Saetran, 2016;Dai and Zhou, 2017;Eisner et al., 2009;Kaso, 2018;Lainer et al., 2016;Le Pichon et al., 2019). To achieve pair-wise correlation analysis between the sampling sites collected during the same periods, the original raw daily measurements were processed as follows: starting on identical days for each pairs of sites, the arrangement of the original daily data into consecutive 3 d intervals (or 6 d intervals in the case of OPE-ANDRA) and the calculation of the average concentration values for the middle day were performed. The resultant data were used for correlation analysis between site pairs (Table S3).
Multiple linear regression (MLR) was used to assess the strength of the relationships between atmospheric concentrations of particulate SC and local environmental factors including the daily mean relative humidity, nighttime temperature, average and maximum temperature, wind speed, solar radiation, rain levels, and LAI. Because the LAI is a 4 d composite, daily values of the other variables were rescaled to consecutive 4 d averaged values. The linear regression (linear model or lm) package in R was employed for multiple regression analyses. The concentration data were log-transformed to obtain regression residual distributions as close as possible to the normal Gaussian one (Fig. S2). Stepwise forward selection was used to select the predictors that explain well the temporal variation in SC concentrations at the site of Marnaz.
It should be noted that due to the limited availability of external parameters, the environmental factors driving SC atmospheric levels have been extensively investigated for only two monitoring sites with contrasted characteristics: the urban background site of Marnaz located in an Alpine valley, and the rural OPE-ANDRA site surrounded by field crop areas spreading over several tens of kilometers.
3 Results and discussion 3.1 Example of spatial coherence of the concentrations at different scales Our previous work (Samaké et al., 2019) showed that particulate polyols and glucose are ubiquitous primary compounds with nonrandom spatial and seasonal variation patterns over France. Here, an inter-site comparison of their short-term concentration evolutions has been carried out at different space scales (from local to national) for the pairs that can be investigated in our database. Figure 2 presents some of these comparisons for three spatial scales (15, 120, and 205 km). The daily average concentrations of polyols (defined as the sum of arabitol and mannitol) and glucose display highly synchronous evolutional trends (i.e., homogeneity in the concentrations, the timing of concentration peaks, simultaneity of the daily specific episodes of increase or decrease in concentrations) over three neighboring monitoring sites located 15 km apart in the Grenoble area ( Fig. 2a and b). Interestingly, remarkable synchronous patterns both for the short term (near-daily) and longer term (seasonal) still occur for sites located 120 km apart, as exemplified for two sites in Alpine environments (Grenoble and Marnaz) ( Fig. 2c  and d). However, as shown in Fig. 2e and f, the evolutions of concentrations become quite dissimilar and asynchronous in terms of seasonal and daily fluctuations for more distant sites (Grenoble and Nice, 205 km apart) that are located in different climatic regions (Alpine for Grenoble, Mediterranean for Nice). This contrasts with results from the rural background site of OPE-ANDRA and the suburban site of Nogent-sur-Oise, both located in a large field crop region of extensive agriculture and about 230 km apart from each other (Fig. 2g). Indeed, they present very similar variations in daily concentrations for multiyear series, despite their distance apart, with concentration peaks generally more pronounced at the rural site of OPE-ANDRA.
The following sections are dedicated to the investigation of the processes that can lead to these similarities and differences according to these spatial scales. 3.2 Inter-site correlations and spatial-scale variability Figure 3a and b provide an overview of the cross-correlation coefficients for the daily evolution of concentrations (for polyols and glucose (SC)) between pairs of sites located at multiple increasing space scales across France (Table S3). Time series of concentrations for both SCs show a clear distance-dependent correlation. The strength of the correlations is highly significant for distances up to 150-190 km (R > 0.72, p < 0.01) and gradually decreases with increasing inter-site distances. One exception is the pair OPE-ANDRA and Nogent-sur-Oise (high correlation for a distance above 230 km), both sites being located in highly impacted agricultural areas. This overall pattern suggests that the processes responsible for the atmospheric concentrations of SC present a spatial homogeneity over typical areas of at least several tens of kilometers.
Unlike SC, ambient air concentrations of sulfate (Fig. 3c) and ammonium (Fig. 3d), associated with long-range aerosol transport (Abdalmogith and Harrison, 2005;Amato et al., 2016;Coulibaly et al., 2015;Pindado and Perez, 2011;Waked et al., 2014), and levoglucosan (Fig. 3e), associated with biomass burning Xiao et al., 2018), display stronger positive correlations (R > 0.72-0.98, p < 0.01) at all pairs of sites considered in the present work. The concentrations of levoglucosan and those of SC clearly display very different annual atmospheric evolution cycles: i.e., higher concentrations of levoglucosan in France are observed in the coldest season (winter) due to the increased biomass burning, while those of SC are observed in warm seasons and coinciding with negligible ambient concentrations of levoglucosan (Fig. S3). Moreover, ambient concentrations of calcium (Fig. 3f), associated with local fugitive dust sources or/and long-range aerosol transport (Ram et al., 2010;Wan et al., 2019) display random correlation patterns. These results are in agreement with Zhu et al. (2018b), who also reported nonsignificant correlations between SC and sulfate in PM 2.5 aerosols measured at Shanghai, China. The distinct spatial behaviors between sulfate (or Ca 2+ ) and SC in the present work further suggest a dominant regional influence for atmospheric SC, as opposed to processes associated with either local sources for calcium or long-range transport for sulfate.
Mannitol and arabitol are well-known materials of fungal spores, serving as osmoregulatory solutes (Medeiros et al., 2006b;Simoneit et al., 2004;Verma et al., 2018;Xiao et al., 2018;Zhang et al., 2015). Based on parallel measurements of spore counts and PM 10 polyol concentrations at three sites within the area of Vienna (Austria), Bauer et al. (2008) found an average arabitol and mannitol content per fungal spores of, respectively, 1.2 pg per spore (range 0.8-1.8 pg per spore) and 1.7 pg per spore (range 1.2-2.4 pg per spore). Mannitol and arabitol have also often been identified in the green algae and lower plants Di Filippo et al., 2013;Gosselin et al., 2016;Vélëz et al., 2007;Xu et al., 2018;Zhang et al., 2010). Gosselin et al. (2016) observed a relatively low (R 2 = 0.31) to high (R 2 = 0.84) coefficient of determination between mannitol and arabitol for total suspended particles (TSPs) collected at a pine-forested area during dry and rainy periods, respectively. High correlation in rainy periods possibly suggested that both chemical species in the TSP fraction in this pineforested area could have been derived mainly from the same sources, i.e., actively wet-discharged ascospores and basidiospores, while the relatively poor correlation in dry periods could likely have been due to more complex sources (dry discharged spores, plants, algae, etc.). Being important chemical species for the metabolism of the microorganisms (Shcherbakova, 2007), it may well be that the concentration ratio of mannitol to arabitol could deliver some information on the spatial or temporal evolution of their emission processes (Gosselin et al., 2016). The annual average mannitolto-arabitol ratio at all sites is about 1.15 ± 0.59, with ratios for the warm period (June-September) being 1 to 2 times higher than those in the cold period (December-May) (Table S1). These ratios are within the range of those previously reported for PM 10 aerosols collected at various urban and rural background sites in Europe (Bauer et al., 2008;Yttri et al., 2011b). Similarly, Burshtein et al. (2011) also reported comparable ratios for PM 10 aerosols collected during autumn and winter from a Mediterranean region in Israel.
Similarly, the annual average glucose-to-polyols ratio at all sites is about 0.79 ± 0.77. No literature data are currently available for comparison. Further work is needed to relate these variations with microorganism communities and plant growing stages.
However, as evidenced in Fig. 4, both mannitol-toarabitol and glucose-to-polyols ratios show a clear distancedependent correlation, with higher correlations (R = 0.64 to 0.98, p < 0.01) observed for pairs of sites within 150-190 km distance. This spatial consistency highlights once again that the dominant emission processes should be effective regionally, rather than being specific local input processes, and that atmospheric dynamics of the concentration levels (i.e., driven by the interplay of emission and removal processes) are determined by quite similar environmental factors (e.g., meteorological conditions, vegetation, land use) at such a regional scale. This implies that local events and phenomena, such as the mechanical resuspension of topsoil and associated biota (like bacteria, fungi, plant materials), might not be their major atmospheric input processes, particularly in urban background areas typically characterized by less bare soil and with a variable nature of the unpaved topsoil at the regional scale (Karimi et al., 2018). Furthermore, Karimi et al. (2018) also recently reported heterogeneous topsoil microbial structure within patches of 43 to 260 km across different regions of France. It follows that the hypotheses of emissions related to the mechanical resuspension of topsoil particles and associated biota or microbiota emitted actively from surface soil into the air generally assumed in most pioneering reports (Medeiros et al., 2006b;Rogge et al., 2007;Simoneit et al., 2004;Wan and Yu, 2007)   Alternatively, the vegetation leaves have also been suggested as sources of atmospheric SCs (Bozzetti et al., 2016;Golly et al., 2018;Jia et al., 2010;Myriokefalitakis et al., 2017;Pashynska et al., 2002;Sullivan et al., 2011;Verma et al., 2018;Wan et al., 2019). In fact, vascular plant leaf surfaces is an important habitat for endophytic and epiphytic microbial communities (Kembel and Mueller, 2014;Lindow and Brandl, 2003;Lymperopoulou et al., 2016;Mhuireach et al., 2016;Whipps et al., 2008). Our results are more in agreement with a dominant atmosphere entrance process closely linked to vegetation, which is more homogeneous than topsoil at the climatic regional scale. Consistent with this, Sullivan et al. (2011) also observed evident distinct regional patterns for daily PM 2.5 polyol and glucose concentrations at 10 urban and rural sites located in the upper Midwest (USA). The authors attributed such a spatial pattern to the differences in vegetation types and microbial diversity over distinct geographical regions. Accordingly, the vegetation struc-ture and composition have previously been shown to play essential roles on airborne microbial variabilities in nearby areas (Bowers et al., 2011;Laforest-Lapointe et al., 2017;Lymperopoulou et al., 2016;Mhuireach et al., 2016).

Influence of the vegetation on polyol and glucose concentrations
The relationships between SC PM 10 concentrations and vegetation (plant materials) can be examined at the site of Grenoble Les Frênes (Grenoble_LF) by comparing the annual evolutions of SC and the free atmospheric cellulose concentrations, together with LAI ones. The daily ambient concentration levels of SC and cellulose range, respectively, from 5.0 to 301.9 ng m −3 (with an average of 41.2±39.9 ng m −3 ) and 0.7 to 207.2 ng m −3 (with an average of 52.9 ± 44.2 ng m −3 ), which corresponds to, respectively, to 0.1 % to 6.6 % and 0.01 % to 5.3 % of total OM mass in PM 10 . These values are comparable to those previ-ously reported for various sites in Europe (Daellenbach et al., 2017;Sánchez-Ochoa et al., 2007;Vlachou et al., 2018;Yttri et al., 2011b). Thus, a major part of PBOA could possibly be ascribed to cellulose and SC-derived sources.
As evidenced in Fig. 5a, ambient free cellulose concentrations vary seasonally, with maximum seasonal average values observed in summer (81.4 ± 47.6 ng m −3 ) and autumn (64.2 ± 49.2 ng m −3 ), followed by spring (52.6 ± 37.8 ng m −3 ) and lower levels in winter (23.0±19.9 ng m −3 ). This is the same global pattern for polyols, which are also more abundant in summer (82.4 ± 47.4 ng m −3 ) and autumn (48.7 ± 41.6 ng m −3 ), followed by spring (24.9 ± 16.3 ng m −3 ) and winter (10.2 ± 9.6 ng m −3 ) in the Grenoble area. On a daily scale, the episodic increases or decreases in polyols in PM 10 are very often well synchronized with that in cellulose (Fig. 5a). Moreover, the maximum atmospheric concentrations of polyols also mainly occur when the vegetation density (LAI) is at its highest in late summer (Fig. 5b). Similar global behaviors are also observed for atmospheric particulate glucose and LAI ( Fig. 5a and b). To further assess the relationships between SC PM 10 concentrations and vegetation in a rural area, a 2-year measurement of cellulose concentrations at the highly impacted agricultural rural site of OPE-ANDRA has been conducted. The average concentration of cellulose at OPE-ANDRA (197.9 ± 217.8 ng m −3 ) is 3.5 times higher than that measured in the urban area of Grenoble. In terms of temporal dynamics, the evolution cycles (i.e., peaks and decreases) of both polyols and glucose are also very often well synchronized with that of cellulose at OPE-ANDRA (Fig. 5c).
Altogether, these findings highlight that SC in PM 10 and cellulose in both urban background and rural agricultural areas most probably share a common source related to the vegetation. This is additional evidence in support of the hypothesis suggested in previous studies (Bozzetti et al., 2016;Burshtein et al., 2011;Daellenbach et al., 2017;Pashynska et al., 2002;Verma et al., 2018;Vlachou et al., 2018;Wan and Yu, 2007;Yttri et al., 2007). It is also in line with studies indicating that the PBOA source profile identified using offline aerosol mass spectrometry (offline-AMS) correlates very well with coarse cellulose concentrations (Bozzetti et al., 2016;Vlachou et al., 2018). A noticeable contribution of cellulose to PBOA mass (26 %) at the rural background site of Payerne (Switzerland), during summer 2012 and winter 2013, was reported by Bozzetti et al. (2016).
As also evidenced in Fig. 5, the cellulose concentration peaks are not systematically correlated to those of polyols. The development stage of the plants (developing or mature leaves, flowering plants) in addition to the metabolic activities of endophytic and epiphytic biota (growth, sporulation), all closely related to meteorological conditions (Bodenhausen et al., 2014;Bringel and Couée, 2015;Lindow and Brandl, 2003;Pirttilä and Frank, 2011;Reddy et al., 2017), could explain such observations. The influence of local meteorological conditions for an urban Alpine valley site is discussed in Sect. 3.4. Consistent with our observations, previous studies conducted at various urban background sites in Europe have suggested that particulate polyols are associated to mature plant leaves and microorganisms (bacterial and fungal spores), while glucose, which is a monomer of cellulose, would most likely be linked to the developing leaves (Bozzetti et al., 2016;Burshtein et al., 2011;Pashynska et al., 2002;Yttri et al., 2007;Zhu et al., 2015).

Influence of meteorological parameters on ambient concentrations of polyols and glucose
We used here a multiple linear regression analysis (MLR) approach to gain further insight into the environmental factors influencing the annual and short time variation cycles of atmospheric SC concentrations. This tentative MLR analysis is focused on the urban background site of Marnaz only since meteorological and other data are readily available for this site and are not influenced too much by some large city effects. Several variables were tested, which are already mentioned in the literature as drivers of SC concentrations. It includes the ambient relative humidity, rainfall level, wind speed, solar radiation, nighttime temperature, average (or maximum) temperature, and LAI. Nighttime temperature was selected since the time series in Marnaz and Grenoble indicate that the major drop in concentrations in late fall (Fig. 2c) is related to the first night of the season with nighttime temperature below 5 • C. The use of the night temperature is also consistent with the bimodal distribution of polyols during nighttime and daytime found in previous studies (Claeys et al., 2004;Graham et al., 2003;Yan et al., 2019;Yttri et al., 2011a). Overall, the environmental factors including the mean nighttime temperature, relative humidity, wind speed, and the leaf area index explain up to 82 % (adjusted R 2 = 0.82; see Table 1) of the annual temporal variation cycles of SC concentrations. The mean nighttime temperature and LAI contribute, respectively, 54 % and 37 % of the observed annual variabilities of SC concentrations. The atmospheric humidity is also a driver for these chemical species (3 % of the explained variation). These results are consistent with previous studies showing that concentrations of mannitol (in both PM 10 and PM 2.5 size fractions) linearly correlate best with the LAI, atmospheric water vapor, and temperature (Heald and Spracklen, 2009;Hummel et al., 2015;Myriokefalitakis et al., 2017). All of these drivers have been previously shown to induce the initial release and influence the long-term airborne microbial (i.e., bacteria, fungi) concentrations (China et al., 2016;Elbert et al., 2007;Grinn-Gofroń et al., 2019;Jones and Harrison, 2004;Rathnayake et al., 2017;Zhang et al., 2015).
Besides, the wind speed (range of 0.2 to 5.6 m s −1 ) seems an additional effective driver affecting the contribution of the local vegetation to SC concentrations in the atmosphere. Although enough air movement is required to passively re-Atmos. Chem. Phys., 19, 11013-11030, 2019 www.atmos-chem-phys.net/19/11013/2019/ Contributions of predictor are normalized to sum 1. The "Relaimpo package under R" was used to compute bootstrap confidence intervals for the importance of effective predictors (n = 1000) (Grömping, 2006 lease microorganisms along with plant debris into the atmosphere, strong air motions induce higher dispersion. These observations are in good agreement with those previously reported (Jones and Harrison, 2004;Liang et al., 2013;Zhang et al., 2010Zhang et al., , 2015Zhu et al., 2018b). For instance Liang et al. (2013) have found a negative correlation between wind speed and polyol concentrations, and the highest atmospheric fungal spores concentrations were observed for a wind speed range of 0.6 to 1.0 m s −1 .
One of the limitations of this study is that 4 d averaged observations do not allow us to evaluate the driver contributions that might explain some short-term events for which the influence of meteorological parameters such as rainfall or solar radiation could also be significant (Grinn-Gofroń et al., 2019;Heald and Spracklen, 2009;Jones and Harrison, 2004). However, such simple parameterizations could be a first step in considering SC concentrations in chemical-transport models (CTMs), and further work is required in this direction in order to generate a robust parameterization of the emissions.

Specific case of a highly impacted agricultural area
This section focuses on evidencing the environmental drivers of PM 10 SC concentrations specific to agricultural areas. To achieve this objective, the site of OPE-ANDRA was selected because it is extensively impacted by agricultural activities, without being too prone to influences by other sources. OPE-ANDRA is a specific rural background observatory located about 230 km east of Paris at an altitude of 392 m. It is characterized by a low population density (< 22 inhabitants km −2 within an area of 900 km 2 ), with no surrounding major transport road or industrial activities. The air monitoring site itself lies in a "reference sector" of 240 km 2 , in the middle of a field crop area (tens of kilometers in all directions). The daily agricultural practices within this reference sector are recorded and made available by ANDRA. The parcels within the agricultural area are submitted to a 3-year crop-rotation system. The major crops are wheat, barley, rape, pea, and sunflower. Additionally, OPE-ANDRA is also characterized by a homogeneous type of soil, with a predominance of superficial clay limestone. Figure 6 shows the daily evolution of polyol concentrations in the PM 10 fraction at OPE-ANDRA from 2012 to 2018, together with the agricultural activities recorded daily and averaged over 12 d.
Although the concentrations of polyols fluctuate from one year to another, they display clear annual variation cycles, with higher values in the warm periods (June-November) and lower concentration values in the cold periods (October-May). Interestingly, the annual concentrations of polyols in 2015 (4.2-111.7 ng m −3 ; annual average: 37.0 ± 29.1 ng m −3 ) are significantly lower than those observed for the other years (0.6-1084.6 ng m −3 ; annual average: 62.9 ± 96.8 ng m −3 ). Similar inter-annual evolution trends, but with variable intensities, are also observed for glucose concentrations (Fig. 6). The year 2015 was found to be particularly hot and dry at OPE-ANDRA (Fig. 7), whereas the locally averaged wind conditions are quite stable over the years within the period of study, suggesting that the wind conditions are not the main driver of the observed inter-annual variability. These results highlight that ambient air temperature and humidity are key meteorological drivers of the annual variation cycles of polyol and glucose concentrations. Hot and dry ambient air conditions may decrease the metabolic activity of the microorganisms (e.g., microbial growth and sporulation) (Fang et al., 2018;Liang et al., 2013;Meisner et al., 2018).  Finally, maximum ambient concentration levels for both SC and cellulose are observed in excellent temporal agreement with the harvest periods (late summer) at the OPE-ANDRA site (Fig. 6). Harvesting activities have been previously reported as the major sources for particulate polyols and glucose to the atmosphere in agricultural and nearby urbanized areas (Golly et al., 2018;Rogge et al., 2007;Simoneit et al., 2004). Hence, the resuspension of plant materials (crop detritus, leaves debris) and associated microbiota (e.g., bacteria, fungi) originating from cultivated lands is the most likely major input process of PM 10 polyols and glucose at field crop sites.

Conclusions
The short-term temporal (daily) and spatial (local to nationwide) evolutions of particulate polyol (defined here as the sum of arabitol and mannitol) and glucose concentrations are rarely discussed in the current literature. The present work aimed at investigating the spatial behavior of these chemicals and evidencing their major effective environmental drivers. The major results mainly showed the following.
The short-term evolution of ambient polyol and glucose concentrations is highly synchronous across an urban city scale and remains very well correlated throughout the same geographic areas of France, even if the monitoring sites are situated in different cities about 150-190 km apart. However, Atmos. Chem. Phys., 19, 11013-11030, 2019 www.atmos-chem-phys.net/19/11013/2019/ sampling sites located in two distinct geographic areas are poorly correlated. This indicates that emission sources of these chemicals are uniformly distributed, and their accumulation and removal processes are driven by quite similar environmental parameters at the regional scale. Therefore, local phenomena, such as the atmospheric resuspension of topsoil particles and associated microbiota and microbial direct emissions (e.g., sporulation), cannot be the main emission processes of particulate polyols and glucose in urban areas not directly influenced by agricultural activities. The atmospheric concentrations of polyols (or glucose) and cellulose display remarkably synchronous temporal evolution cycles at the background urban site of Grenoble, indicating a common source related to plant debris.
Higher ambient concentrations of polyols and glucose at the rural site of OPE-ANDRA occur during each harvest period, pointing out resuspension processes of plant materials (crop detritus, leaves debris) and associated microbiota for agricultural and nearby urbanized areas. This is associated with higher PM 10 cellulose concentration levels, as high as 0.4 to 2.0 µg m −3 on a daily basis (accounting up to 7.5 % to 32.4 % of the OM mass).
Multiple linear regression analysis of the yearly series from the site of Marnaz gave insightful information on which parameter controls the ambient concentrations of polyols and glucose. Ambient air nighttime temperature, relative humidity, and vegetation density are the most important drivers, whilst wind speed conditions tend to affect the contribution of local vegetation.
Altogether, these results improve our understanding of the spatial behavior tracers of PM 10 PBOA emission sources in France and in general, which is imperative for further implementation of this important mass fraction of OM into chemical transport models. Further investigations of the airborne microbial fingerprint (bacteria and fungi) are ongoing, which may deepen our understanding of the PBOA source profile.
Data availability. All relevant data for this paper are archived at the IGE (Institut des Géosciences de l'Environnement), and are available upon request from the corresponding author (Jean-Luc Jaffrezo).
Author contributions. JLJ was the (co-)supervisor for the PhDs for AS, FC, and SW and for the post-docs of DS, BG, and AW. He directed all the personnel who performed the analysis at IGE. He is the coordinator for the CNRS LEFE-EC2CO CAREMBIOS program that is funding the work of AS. GU and JMFM were the cosupervisor for the PhDs of AS and SW. EP, OF, and VR supervised the PhD of DMO, who investigated the sites in northern France. OF, JLJ, JLB, AA, and NM were the coordinators or partners of the different initial programs for the collection and chemical analysis of the samples. VJ developed the analytical techniques for polyol and cellulose measurements. TC performed the cellulose measurements. Sample analyses at LSCE were performed by NB. AC gave advice on the statistical aspects of the data processing. AS and JLJ processed the data and wrote up the paper. SW participated in the visualization of the results. SC supervises the OPE station and collected the agricultural activities records. All authors from AASQA (author affiliation nos. 9 to 16) are representatives of each network that conducted the sample collection and the general supervision of the sampling sites. All authors reviewed and commented on the paper.
Competing interests. The authors declare that they have no conflict of interest.
Acknowledgements. We would like to express special acknowledgement of Pierre Taberlet (LECA, Grenoble, France) for fruitful discussions about the importance of endophytic and epiphytic biota for aerobiology. The PhDs of Abdoulaye Samaké and Samuël Weber are funded by the Government of Mali and ENS Paris, respectively. We gratefully acknowledge the LEFE-CHAT and EC2CO programs of the CNRS for financial support of the CAREMBIOS multidisciplinary project and the LEFE-CHAT program for the MECEA project for the development of the atmospheric cellulose measurements. Samples were collected and analyzed within the framework of many different programs funded by ADEME, Primequal, the French Ministry of Environment, the CARA program led by the French Reference Laboratory for Air Quality Monitoring (LCSQA), ANDRA, and actions funded by many AASQAs and IMT Lille Douai (especially Labex CaPPA ANR-11-LABX-0005-01 and CPER CLIMIBIO projects). Analytical aspects were supported at IGE by the Air-O-Sol platform within Labex OSUG@2020 (ANR10 LABX56). We acknowledge the work of many engineers in the lab at IGE for the analyses (Aude Wack, Céline Charlet, Fany Donaz, Fany Masson, Sylvie Ngo, Vincent Lucaire, Claire Vérin, and Anthony Vella). Finally, the authors would like to kindly thank the dedicated efforts of many other people at the sampling sites and in the laboratories for collecting and analyzing the samples.
Financial support. This research has been supported by the Centre National de la Recherche Scientifique (LEFE-CHAT and EC2CO programs: CAREMBIOS multidisciplinary project).
Review statement. This paper was edited by Alex Huffman and reviewed by three anonymous referees.