Seasonal variation and origins of volatile organic compounds observed during 2 years at a western Mediterranean remote background site (Ersa, Cape Corsica)

Abstract. An original time series of about 300 atmospheric measurements of a wide range of volatile organic compounds (VOCs) was obtained at a
remote Mediterranean station on the northern tip of Corsica (Ersa, France) over 25 months from June 2012 to June 2014. This study presents the
seasonal variabilities of 35 selected VOCs and their various associated sources. The VOC abundance was largely dominated by oxygenated VOCs (OVOCs)
along with primary anthropogenic VOCs with a long lifetime in the atmosphere. VOC temporal variations were then examined. Primarily of local
origin, biogenic VOCs exhibited notable seasonal and interannual variations,
related to temperature and solar radiation. Anthropogenic compounds showed
increased concentrations in winter (JFM months) followed by a decrease in
spring/summer (AMJ/JAS months) and higher winter concentration levels in 2013 than in 2014 by up to 0.3 µg m−3 in the cases of propane,
acetylene and benzene. OVOC concentrations were generally high in
summertime, mainly due to secondary anthropogenic/biogenic and primary
biogenic sources, whereas their lower concentrations during autumn and winter were potentially more influenced by primary/secondary anthropogenic sources.
Moreover, an apportionment factorial analysis was applied to a database
comprising a selection of 14 individual or grouped VOCs by means of the
positive matrix factorization (PMF) technique. A PMF five-factor solution was taken on. It includes an anthropogenic factor (which contributed 39 % to
the total concentration of the VOCs selected in the PMF analysis) connected
to the regional background pollution, three other anthropogenic factors
(namely short-lived anthropogenic sources, evaporative sources, and
long-lived combustion sources, which together accounted for 57 %) originating from either nearby or more distant emission areas (such as Italy
and south of France), and a local biogenic source (4 %). Variations in
these main sources impacting VOC concentrations observed at the Ersa station
were also investigated at seasonal and interannual scales. In spring and
summer, VOC concentrations observed at Ersa were the lowest in the 2-year
period, despite higher biogenic source contributions. During these seasons,
anthropogenic sources advected to Ersa were largely influenced by chemical
transformations and vertical dispersion phenomena and were mainly of
regional origins. During autumn and winter, anthropogenic sources showed higher contributions when European air masses were advected to Ersa and
could be associated with potential emission areas located in Italy and
possibly more distant ones in central Europe. Higher VOC winter
concentrations in 2013 than in 2014 could be related to contribution
variations in anthropogenic sources probably governed by their emission
strength with external parameters, i.e. weaker dispersion phenomena and the pollutant depletion. High-frequency observations collected during several intensive field campaigns conducted at Ersa during the three summers
2012–2014 confirmed findings drawn from bi-weekly samples of the 2-year period
in terms of summer concentration levels and source apportionment. However,
they also suggested that higher sampling frequency and temporal resolution,
in particular to observe VOC concentration variations during the daily
cycle, would have been necessary to confirm the deconvolution of the
different anthropogenic sources identified following the PMF approach.
Finally, comparisons of the 25 months of Ersa observations with VOC
measurements conducted at 17 other European monitoring stations highlighted
the representativeness of the Ersa station for monitoring seasonal
variations in VOC regional pollution impacting continental Europe.
Nevertheless, VOC winter concentration levels can significantly vary between
sites, pointing out spatial variations in anthropogenic source
contributions. As a result, Ersa concentration variations in winter were
more representative of VOC regional pollution impacting central Europe.
Moreover, interannual and spatial variations in VOC winter concentration
levels were significantly impacted by synoptic phenomena influencing
meteorological conditions observed in continental Europe, suggesting that
short observation periods may reflect the variability of the identified
parameters under the specific meteorological conditions of the study period.


Abstract. An original time series of about 300 atmospheric measurements of a wide range of volatile organic compounds (VOCs) was obtained at a remote Mediterranean station on the northern tip of Corsica (Ersa, France) over 25 months from June 2012 to June 2014. This study presents the seasonal variabilities of 35 selected VOCs and their various associated sources. The VOC abundance was largely dominated by oxygenated VOCs (OVOCs) along with primary anthropogenic VOCs with a long lifetime in the atmosphere. VOC temporal variations were then examined. Primarily of local origin, biogenic VOCs exhibited notable seasonal and interannual variations, related to temperature and solar radiation. Anthropogenic compounds showed increased concentrations in winter (JFM months) followed by a decrease in spring/summer (AMJ/JAS months) and higher winter concentration levels in 2013 than in 2014 by up to 0.3 µg m −3 in the cases of propane, acetylene and benzene. OVOC concentrations were generally high in summertime, mainly due to secondary anthropogenic/biogenic and primary biogenic sources, whereas their lower concentrations during autumn and winter were potentially more influenced by primary/secondary anthropogenic sources. Moreover, an apportionment factorial analysis was applied to a database comprising a selection of 14 individual or grouped VOCs by means of the positive matrix factorization (PMF) technique. A PMF five-factor solution was taken on. It includes an anthropogenic factor (which contributed 39 % to the total concentration of the VOCs selected in the PMF analysis) connected to the regional background pollution, three other anthropogenic factors (namely short-lived anthropogenic sources, evaporative sources, and long-lived combustion sources, which together accounted for 57 %) originating from either nearby or more distant emission areas (such as Italy and south of France), and a local biogenic source (4 %). Variations in these main sources impacting VOC concentrations observed at the Ersa station were also investigated at seasonal and interannual scales. In spring and summer, VOC concentrations observed at Ersa were the lowest in the 2-year period, despite higher biogenic source contributions. During these seasons, anthropogenic sources advected to Ersa were largely influenced by chemical transformations and vertical dispersion phenomena and were mainly of regional origins. During autumn and winter, anthropogenic sources showed higher contributions when European air masses were advected to Ersa and could be associated with potential emission areas located in Italy and possibly more distant ones in central Europe. Higher VOC winter concentrations in 2013 than in 2014 could be related to contribution variations in anthropogenic sources probably governed by their emission strength with external parame-Published by Copernicus Publications on behalf of the European Geosciences Union. ters, i.e. weaker dispersion phenomena and the pollutant depletion. High-frequency observations collected during several intensive field campaigns conducted at Ersa during the three summers 2012-2014 confirmed findings drawn from bi-weekly samples of the 2-year period in terms of summer concentration levels and source apportionment. However, they also suggested that higher sampling frequency and temporal resolution, in particular to observe VOC concentration variations during the daily cycle, would have been necessary to confirm the deconvolution of the different anthropogenic sources identified following the PMF approach. Finally, comparisons of the 25 months of Ersa observations with VOC measurements conducted at 17 other European monitoring stations highlighted the representativeness of the Ersa station for monitoring seasonal variations in VOC regional pollution impacting continental Europe. Nevertheless, VOC winter concentration levels can significantly vary between sites, pointing out spatial variations in anthropogenic source contributions. As a result, Ersa concentration variations in winter were more representative of VOC regional pollution impacting central Europe. Moreover, interannual and spatial variations in VOC winter concentration levels were significantly impacted by synoptic phenomena influencing meteorological conditions observed in continental Europe, suggesting that short observation periods may reflect the variability of the identified parameters under the specific meteorological conditions of the study period.

Introduction
The main trace pollutants in the atmosphere encompass a multitude of volatile organic compounds (VOCs), with lifetimes varying from minutes to months (e.g. Atkinson, 2000). Their distribution is principally due to (i) multiple natural and anthropogenic sources, which release VOCs directly into the atmosphere. At a global scale, natural emissions are quantitatively larger than anthropogenic ones (Guenther et al., 2000), and the largest natural source is considered to be the vegetation (Finlayson-Pitts and Pitts, 2000;Guenther et al., 2000Guenther et al., , 2006. In urban areas, numerous anthropogenic sources can abundantly emit various VOCs (Friedrich and Obermeier, 1999). Once in the atmosphere, VOC temporal and spatial variabilities are notably influenced by (ii) mixing processes along with (iii) removal processes or chemical transformations (Atkinson, 2000;Atkinson and Arey, 2003). Accordingly, with a view to thoroughly characterizing VOC sources, it is meaningful to examine their chemical composition and identify the factors controlling their variations at different timescales.
VOC regional distributions change considerably as a result of various confounding factors, namely the emission strength of numerous potential sources, diverse atmospheric lifetimes and removal mechanisms, transport processes and fluctua-tions in meteorological conditions. Therefore, these elements underline the necessity of carrying out long-term VOC measurements. Growing efforts are currently being made to conduct European background measurements over several seasons (e.g. Seco et al., 2011), 1 year (such as Helmig et al., 2008;Legreid et al., 2008) and even several years (Solberg et al., 1996, and Tørseth et al., 2012, at several European sites; Hakola et al., 2006, andHellén et al., 2015, in Scandinavia;Dollard et al., 2007;Grant et al., 2011, andMalley et al., 2015, in the United Kingdom;Borbon et al., 2004;Sauvage et al., 2009, andWaked et al., 2016, in France;Plass-Dülmer et al., 2002, in Germany;Navazo et al., 2008, in the Iberian Peninsula;Lo Vullo et al., 2016 in Italy). These research studies principally explored the effectiveness of emission regulations and links between tropospheric ozone production and VOC concentration levels. They also assessed seasonal variations and regional distributions of VOC concentrations. Nonetheless, investigations of principal factors governing temporal and spatial variations in VOC concentration levels in the European background atmosphere remain scarce. However, the consideration of the influence of (i) source emission strength variations (built upon a factorial analysis -e.g. Lanz et al., 2009;Lo Vullo et al., 2016), (ii) long-range transport of pollution (e.g. by the examination of air-mass trajectories combined with concentrations measured at a study site; Sauvage et al., 2009) and (iii) fluctuations in meteorological conditions (which are prone to dispersing the pollutants over long distances by convective and advective transport) can supply relevant information to deal more in depth with the evaluation of seasonal variations and regional distribution of VOC concentrations in the European background atmosphere.
Particulate and gaseous pollutants detrimentally affect the Mediterranean atmosphere. Accordingly, they are prone to increasing aerosol and/or ozone concentration levels in the Mediterranean, regularly higher compared to most regions of continental Europe, and primarily during summer (Doche et al., 2014;Nabat et al., 2013;Safieddine et al., 2014). The Mediterranean region is known to be a noteworthy climate change "hotspot" which is expected to go through severe warming and drying in the 21st century (Giorgi, 2006;Kopf, 2010;Lelieveld et al., 2014). As a consequence, this can have serious consequences for the release of VOCs from biogenic and anthropogenic sources along with their fate in the atmosphere, with uncertain predicted impacts (Colette et al., 2012(Colette et al., , 2013Jaidan et al., 2018). Actually, the examination of air composition, concentration levels and trends in the Mediterranean region continues to be challenging, primarily due to the lack of extensive in situ observations. Given this context, as part of the multidisciplinary regional research programme MISTRALS (Mediterranean Integrated Studies at Regional and Local Scales; http://mistrals-home. org/, last access: 11 October 2020), the project ChArMEx (the Chemistry-Aerosol Mediterranean Experiment; Dulac, 2014) aims at assessing the current and future state of the at-mospheric environment in the Mediterranean along with examining its repercussions for the regional climate, air quality and marine biogeochemistry. Within the framework of ChArMEx, several observation periods were conducted at the Ersa station, a remote site considered to be representative of the north-western Mediterranean basin, in order to explain variations in VOC concentrations affecting the western Mediterranean atmosphere. Michoud et al. (2017) characterized the variations in VOC concentrations observed at Ersa in summer 2013 (from 15 July to 5 August 2013) by identifying and examining their sources.
In this article, we have presented and discussed the factors controlling seasonal and interannual variations of a selection of VOCs observed at the Ersa station over more than 2 years (from early June 2012 to late June 2014). To this end, this study describes (i) the concentration levels of the targeted VOCs, (ii) their temporal variations at seasonal and interannual scales, (iii) the identification and characteristics of their main sources by statistical modelling, (iv) the evaluation of their source contributions on seasonal bases, together with (v) the representativeness of the Ersa station in terms of seasonal variations in VOC concentrations impacting continental Europe.

Study site
Located in the north-western part of the Mediterranean Sea, Corsica is a French territory situated 11 km north of the Sardinian coasts, 90 km east of Tuscany (Italy) and 170 km south of the French Riviera (France). The fourth largest Mediterranean island, its land corresponds to an area of 8681 km 2 encompassed by around 1000 km of coastline (Encyclopaedia Britannica, 2018). Corsica contrasts to other Mediterranean islands due to the importance of its forest cover (about a fifth of the island).
Within the framework of the ChArMEx project, an enhanced observation period was set up at a ground-based station in the north of Corsica (Ersa; 42.969 • N, 9.380 • E) over 25 months, from early June 2012 to late June 2014. The aim was to provide a high-quality controlled climatically relevant gas/aerosol database following the recommendations and criteria of international atmospheric chemistry networks, i.e. the Aerosol, Clouds and Trace gases Research Infrastructure (ACTRIS -https://www.actris.eu/, last access: 11 October 2020), the European Monitoring and Evaluation Program (EMEP -http://www.emep.int/, last access: 11 October 2020; Tørseth et al., 2012), and the Global Atmosphere Watch of the World Meteorological Organization (WMO-GAW -http://www.wmo.int/pages/prog/arep/ gaw/gaw_home_en.html, last access: 11 October 2020). The Ersa remote site is part of the Corsican Observatory for Research and Studies on Climate and Atmosphere-ocean en-vironment (CORSiCA -https://corsica.obs-mip.fr/, last access: 11 October 2020; Lambert et al., 2011) and is located on the highest point of a ridge equipped with windmills (see the orographic description of the surroundings in Cholakian et al., 2018), at an altitude of 533 m above sea level (a.s.l.). Given its position in the north of the 40 km long Cape Corsica peninsula (Fig. 1), the Mediterranean Sea is clearly visible from the sampling site on the western, northern, and eastern sides (2.5-6 km from the sea; see also the figure presented in Michoud et al., 2017). The station was initially set up in order to monitor and examine pollution advected to Ersa by air masses transported over the Mediterranean and originating from the Marseille-Fos-Berre region (France; Cachier et al., 2005), the Rhone Valley (France), and the Po Valley (Italy; Royer et al., 2010), namely largely industrialized regions. The Ersa station is about 30 km north of Bastia (Fig. 1), the second largest Corsican city (44121 inhabitants; census 2012) and the main harbour. An international airport (Bastia-Poretta) is located 16 km further south of Bastia city centre. More than 2 million passengers transited in Corsica via Bastia during the tourist season (May-September) in 2013 (ORT Corse, 2013). However, as the Cape Corsican peninsula benefits in the south from a mountain range (peaking between 1000 and 1500 m a.s.l.) acting as a natural barrier, the sampling site is therefore not affected by transported pollution originating from Bastia. Only small rural villages and a small local fishing harbour (Centuri) are found within 5 km of the measurement site. Additionally, the Ersa station is accessible by a dead-end road serving only the windmill site, surrounded by vegetation made up of Mediterranean maquis, a shrubland biome characteristically consisting of densely growing evergreen shrubs, and also roamed by a herd of goats from a nearby farm. Some forests (78 % of holm oaks, with some cork oaks and chestnuts) are also located nearby, thus ensuring that local anthropogenic pollution does not contaminate in situ observations. As a result, the Ersa station can be characterized as a remote background Mediterranean site.

VOC measurements
During a study period of 2 years, non-methane hydrocarbons (NMHCs) and OVOCs (carbonyl compounds) were measured, routinely employing complementary offline methods. Four-hour-integrated (09:00-13:00 or 12:00-16:00 UTC) ambient air samples were collected bi-weekly (every Monday and Thursday) into steel canisters and on sorbent cartridges. The inlets were roughly 1.5 m above the roof of a container mainly housing trace gas analysers. Table 1 describes VOC measurements set up throughout the observation period and Fig. S1 specifies their collection periods.
As generally realized in the EMEP network, 24 C 2 -C 9 NMHCs were collected into Silcosteel canisters of a volume of 6 L, conforming to the TO-14 technique, which is consid- ered adequate for the measurement of many non-polar VOCs (US-EPA, 1997); 152 air samples were realized with a homemade device (PRECOV) for sampling air at a steady flow rate regulated to 24 mL min −1 by canisters previously placed under vacuum. NMHC analysis was performed by a gas chromatograph coupled with a flame ionization detector (GC-FID) within 3 weeks following the sampling. Separation was performed by a system of dual-capillary columns supplied with a switching device: the first column was a CP Sil5CB (50 m × 0.25 mm × 1 µm), suitable for the elution of VOCs from six to nine carbon atoms, and the other one was a plot Al 2 O 3 / Na 2 SO 4 (50 m × 0.32 mm × 5 µm), in order to effectively elute VOCs from two to five carbon atoms. Four main steps constituted the quality assurance/quality control programme: (i) the implementation of standard operating procedures, (ii) canister cleaning and certification (blank levels < 0.02 ppb), (iii) regular intercomparison exercises and (iv) sampling tests carried out in field conditions and concomitant with in situ measurements (Sauvage et al., 2009).
About 150 air samples were gathered using sorbent cartridges (63 air samples on multi-sorbent cartridges and 89 additional ones on 2,4-dinitrophenylhydrazine -DNPHcartridges), by means of an automatic clean room sampling system (ACROSS, TERA Environment, Crolles, France). C 1 -C 16 VOCs were collected via a 0.635 cm diameter 3 m long PFA line. They were then trapped into one of the two cartridge types: a multi-sorbent one consisting of carbopack C (200 mg) and carbopack B (200 mg; marketed under the name of carbotrap 202 by Perkin-Elmer, Wellesley, Massachusetts, USA) or a Sep-Pak DNPH-Silica one (proposed by Waters Corporation, Milford, Massachusetts, USA). These offline techniques are further characterized in Detournay et al. (2011), and their satisfying use in situ has already been discussed by Detournay et al. (2013) and Ait-Helal et al. (2014). Succinctly here, the sampling  Detournay et al. (2013) Detournay (2011), Detournay et al. (2011) of 44 C 5 -C 16 NMHCs, comprising alkanes, alkenes, aromatic compounds and six monoterpenes, as well as six C 6 -C 11 n-aldehydes, was conducted at a flow rate fixed at 200 mL min −1 and using the multi-sorbent cartridges. These latter ones were preliminary prepared by means of a RTA oven (French abbreviation for "régénérateur d'adsorbant thermique" -manufactured by TERA Environment, Crolles, France) in order to condition them during 24 h with purified air heated to 250 • C and at a flow rate regulated at 10 mL min −1 . Concomitantly, 15 additional C 1 -C 8 carbonyl compounds were collected at a flow rate fixed at 1.5 L min −1 using the DNPH cartridges. During the field campaign, several ozone scrubbers were successively inserted into the sampling lines in order to limit any eventual ozonolysis of the measured VOCs: a MnO 2 ozone scrubber was retained for the multi-sorbent cartridges, while KI ozone scrubber was placed upstream of the DNPH cartridges. Moreover, stainless-steel particle filters of 2 µm diameter porosity (Swagelok) were installed in order to prevent particle sampling. Then, VOC samples were transferred to the laboratory to be analysed within 6 weeks of their collection using a GC-FID (for the multi-adsorbent cartridges) or by a highperformance liquid chromatograph connected to an ultraviolet detector (HPLC-UV; for the DNPH cartridges). The reproducibility of each analytical instrument was frequently checked by analysing a standard and examining results by plotting them on a control chart realized for each compound. The VOC detection limit was determined as 3 times the standard deviation of the blank variation. Detection limits in this study were all below 0.05 µg m −3 for the steel canisters and the DNPH cartridges and 0.01 µg m −3 for the multi-sorbent cartridges. The uncertainties for each species were evaluated respecting the ACTRIS-2 guidelines for the uncertainty evaluation (Reimann et al., 2018), considering precision, detection limit and systematic errors in the measurements. Relative uncertainties assessed in this study ranged from 7 % to 43 % for the steel canisters, from 7 % to 73 % for the multi-sorbent cartridges and from 6 % to 41 % for the DNPH cartridges. Finally, the VOC dataset was validated following the ACTRIS protocol (Reimann et al., 2018).
Among the 71 different VOCs monitored at Ersa during the observation period, 35 VOCs were finally selected in this study following the methodology described in Sect. S1 of the Supplement.

Ancillary measurements
Other trace gases (CO and O 3 ) and meteorological parameters were ancillary monitored at the Ersa site during the observation period. CO was measured from 22 November 2012 to 16 December 2013 by a cavity ring-down spectroscopy analyser (G2401; Picarro, Santa Clara, California, USA) at a time resolution of 5 min. O 3 was measured from 31 May 2012 to 26 December 2013 using an ultraviolet absorption analyser (TEI 49i manufactured by Thermo Environmental Instruments Inc., Waltham, Massachusetts, USA) at a time resolution of 5 min. Meteorological parameters (temperature, pressure, relative humidity, wind speed, wind direction and total -direct and diffuse -solar radiation) were measured every minute from 8 June to 14 August 2012 and every 5 min from 15 August 2012 to 11 July 2014, with a weather station (CR1000 manufactured by Campbell Scientific Europe, Antony, France) placed at approximately 1.5 m above an adjacent container roof to the one which housed trace gas instruments. Note that ancillary trace gas and meteorological parameter results presented in this study are 4 h averages concurrent with VOC sampling periods (see Fig. S1).

Identification and contribution of major sources of VOCs
In order to characterize NMHC concentrations measured at Ersa with steel canisters (the reasons for this VOC selection are presented in Sect. S2 of the Supplement), we apportioned them within their sources in this study using the positive matrix factorization approach (PMF; Paatero, 1997;Paatero and Tapper, 1994). The PMF mathematical theory has already been presented in Debevec et al. (2017) and is therefore mentioned in Sect. S2. We used the PMF version 5.0, an enhanced tool developed by the Environmental Protection Agency (EPA) and including the multilinear engine program version 2 (ME-2; Paatero, 1999), and followed the guidance on its use (Norris et al., 2014). Using NMHC inputs composed of 152 atmospheric data points of 14 variables (13 single primary NMHCs and another 1 resulting from the grouping of C 8 aromatic compounds) and following the methodology presented in Sect. S2, a five-factor PMF solution has been selected in this study.

Classification of air-mass origins
In order to identify and classify air-mass origins, we analysed back-trajectories calculated by the online version of the HYSPLIT Lagrangian model (the Hybrid Single Particle Lagrangian Integrated Trajectory Model developed by the National Oceanic and Atmospheric Administration -NOAA -Air Resources Laboratory; Draxler and Hess, 1998;Stein et al., 2015) using Ersa as the receptor site (arrival altitude: 600 m a.s.l.). For each 4 h-atmospheric data point of the field campaign used for the factorial analysis, five back-trajectories of 48 h were computed using GDAS onedegree resolution meteorological data, in order to follow the same methodology as Michoud et al. (2017). The first backtrajectory of a set corresponds to the hour when the air sampling was initiated (i.e. 09:00 or 12:00 UTC -see Table 1) and the four other ones were calculated every following hour. The time step between each point along the back-trajectories was fixed at 1 h. Then, back-trajectories were visually classified. Having several back-trajectories per sample allowed us to check whether air masses transported to the station over 4 h were of the same origin. Samples associated with air masses showing contrasted trajectories (e.g. due to a transitory state between two different origins) were classified as being of mixed origins (9 % of the air masses) and discarded from this study. Remaining air masses were then manually classified into five trajectory clusters (marine, Corsica-Sardinia, Europe, France and Spain - Fig. 2 and Table 2) depending on their pathway when they reached the Ersa station, their residence time over each potential source region and the length of their trajectories. Additionally, air masses of each cluster were subdivided depending on their distance travelled during their 48 h course in order to highlight potential more distant sources from local ones. This subdivision is also given in Table 2 in order to pinpoint differences in transport times.

Identification of potential emission areas
PMF source contributions were coupled with backtrajectories in order to investigate potential emission regions contributing to long-distance pollution transport to the Ersa site. To achieve this, the concentration field (CF) statistical method established by Seibert et al. (1994) was chosen in the present study. The CF principle has already been presented in Debevec et al. (2017) and is therefore mentioned in Sect. S3 of the Supplement. Seventy-two-hour back-trajectories together with meteorological parameters of interest (i.e. precipitation) were retrieved from the GDAS meteorological fields with a PCbased version of the HYSPLIT Lagrangian model (version 4 revised in April 2018), following the same methodology as that used for the 48 h back-trajectories previously presented. The arrival time of trajectories at the Ersa station corresponds to the hour when half of the sampling was carried out (i.e. 11:00 or 14:00 UTC -see Table 1). Note that longer backtrajectories were considered for CF analyses than those for air-mass origin classification (Sect. 2.4.1), in order to be in the same conditions as Michoud et al. (2017) and hence to have comparative results between the two Ersa VOC studies.
CF analyses applied to VOC source contributions were carried out by means of the ZeFir tool (version 3.70; Petit et al., 2017). Back-trajectories were shortened (i.e. the Ze-Fir tool considered shorter back-trajectories than 72 h) when precipitation higher than 0.1 mm was encountered along the trajectory (Bressi et al., 2014). As done by Michoud et al. (2017), back-trajectories have also been shortened when air-mass altitudes exceeded 1500 m a.s.l. in order to discard biases related to the significant dilution impacting air masses reaching the free troposphere. A better statistical significance of the CF results is commonly considered for grid cells with a high number of crossing trajectory points. As a result, some studies applied an empirical weighing function so as to limit the possible influence of elevated concentrations, which may be observed during occasional episodes (e.g. Bressi et al., 2014;Waked et al., 2014Waked et al., , 2018, on cells with a low number of trajectory points. We initially tried to apply this weighing function in this study. Exploratory tests revealed that CF results with the empirical weighing function only highlighted local contributions, given the number of air masses of this study. The farther a cell was from the Ersa station, the lower was its corresponding n ij value (number of points of the total number of back-trajectories contained in the ij th grid cell, Sect. S3 of the Supplement) and the more the weighing function tended toward downweighting the results related to this cell. Therefore, CF results discussed in this study were realized without weighing, and these limitations should be taken into account when examining CF analyses, which are hence considered to be indicative information.
Finally, the spatial coverage of grid cells was set from (9 • W, 32 • N) to (27 • E, 54 • N), with a grid resolution of 0.3 • × 0.3 • . Allocated contributions were smoothed following a factor (corresponding to the strength of a Gaussian filter) set to 5 to take into account the uncertainties in the backtrajectory path (Charron et al., 2000).

Meteorological conditions
Monthly variations in meteorological parameters are depicted in Fig. 3. As the field measurement period covered 2 years (i.e. from June 2012 to June 2014), their interannual variations are also shown in Fig. 3b.
Air temperatures showed typical seasonal variations, i.e. the highest recorded in summer (i.e. the months of July to September) and the lowest in winter (i.e. the months of January to March). They were globally in the range of normal values over the period 1981-2010 determined by Météo-France (the French national meteorological service; normal values correspond to minimal and maximal mean values determined for Bastia and are available for consultation at http://www.meteofrance.fr/climat-passe-et-futur/ climathd, last access: 11 October 2020). Moreover, June temperatures were colder in 2013 than in 2012 and 2014 (mean temperatures of 24.7 ± 5.8, 19.4 ± 4.1 and 22.5 ± 5.4 • C in 2012, 2013 and 2014, respectively), which could have influenced biogenic emissions. Winter temperatures were colder in 2013 than in 2014 (mean temperatures of 7.0 ± 4.1 and 9.7 ± 1.5 • C in 2013 and 2014, respectively). This finding could be explained by different winter climatic events which affected a large part of continental Europe in 2013 and 2014. On the one hand, the European winter was particularly harsh in 2013, caused by changes in air-flux orientation originally due to the sudden stratospheric warming of the stratospheric polar vortex (Coy and Pawson, 2015). On the other hand, most of the western European countries experienced a mild winter in 2014 characterized by its lack of cold outbreaks and nights and caused by an anomalous atmospheric circulation (Rasmijn et al., 2016;Van Oldenborgh et al., 2015;Watson et al., 2016). Solar radiation followed typical seasonal variations, with the highest recorded from May to August and the lowest in December and January. Spring (i.e. the months of April to June) solar radiations were higher in 2014 than in 2013 (average values of 371 ± 157 and 478 ± 153 W m −2 in 2013 and 2014, respectively). Summer solar radiations were higher in 2013 than in 2012 (average values of 332 ± 164 and 395 ± 128 W m −2 in 2012 and 2013, respectively). These solar radiation variations may have affected biogenic VOC (BVOC) emissions and photochemical reactions.
Relative humidity followed opposite seasonal variations in temperature and solar radiation. Air in June was drier in 2012 than in 2013 and 2014 (mean relative humidities of 57 ± 15 %, 77 ± 16 % and 67 ± 33 % in 2012, 2013 and 2014, respectively). The wind speed did not show a clear seasonal variation over the 2 years studied. Wind speeds were slightly higher in April and May, which could induce higher dispersion of air pollutants and favour their advection to the Ersa station by the most distant sources.

Air-mass origins
Occurrences of air-mass origins which influenced Ersa throughout the observation period are indicated in Table 2. The Ersa station was predominantly influenced by continental air masses coming from Europe (corresponding to cluster 3, 31 %), France (cluster 4, 26 %), Corsica-Sardinia (cluster 2, 14 %) and Spain (cluster 5, 5 %) and to a lesser extent by air masses predominantly of marine origin (cluster 1, 15 %). Each of these five clusters is mostly associated with a particular trajectory sector (e.g. south for air masses originating from Corsica and/or Sardinia; see Fig. 2) and is defined by a different transit time from continental coasts. According to Michoud et al. (2017), transit time can be viewed as an indicator of the last time when an air mass could have been enriched in anthropogenic sources (Table 2). Continental air masses spent less time over the sea than marine ones. Transit times of continental air masses over the sea differed depending on how they were categorized. Air masses originating from Corsica-Sardinia, France and Europe spent 0-8 h (median values - Table 2) above the sea before reaching the Ersa station, while the air masses originating from Spain spent about 36 h. These contrasting transit times may denote both distinctive atmospheric processing times for air masses and different oceanic source influences on VOC concentrations observed at the Ersa station.
European and French air masses showed lower transit times over the sea (median values of 6 and 8 h, respectively; Table 2) when their trajectories were categorized as long, compared to short ones (23 and 19 h, respectively). These findings are based on the fact that an air-mass trajectory classified as short has a closer distance between two of its successive trajectory points compared to another one classified as long. Due to the Ersa location in the Mediterranean Sea, air masses with trajectories categorized as long spent longer periods above the sea before reaching the Ersa site. Note that European and French air masses were more frequently characterized by long trajectories (accounting for 20 % of the air masses observed at Ersa during the studied period, for each) than short ones (11 % and 6 %, respectively). Moreover, marine air masses with short and long trajectories both showed long transit times (median values of 40-48 h -Table 2). Corsican-Sardinian air masses were only characterized by long trajectories. radiation in W m −2 , relative humidity in % and wind speed in m s −1 ) represented by box plots; the blue solid line, the red marker, and the box represent the median, the mean, and the interquartile range of the values, respectively. The bottom and top of the box depict the first and third quartiles and the ends of the whiskers correspond to the first and ninth deciles. (b) Their monthly average concentrations as a function of the year. Note that meteorological parameter data used in this study were restricted to periods when VOC measurements were realized.

VOC mixing ratios
Statistical results of concentrations of the 35 VOCs selected in this study (see Sect. S1 in the Supplement) are summarized in Table 3. Their average concentration levels as a function of the measurement sampling time (09:00-13:00 or 12:00-16:00 UTC) are also indicated in Table S1. These VOCs were organized into three principal categories: biogenic, anthropogenic, and oxygenated VOCs (5, 16 and 14 targeted species, respectively; Table 3). Isoprene and four monoterpenes were classified into BVOCs, while the other primary hydrocarbons (alkanes, alkenes, alkynes and aromatic compounds) were included in anthropogenic NMHCs, since their emissions are especially in connection with human activities. OVOCs have been presented separately, as these compounds come from both biogenic and anthropogenic (primary and secondary) sources. OVOCs were the most abundant, accounting for 65 % of the total concentration of the 35 VOCs selected in this study. They were mainly composed of acetone (contribution of 51 % to the OVOC cumulative concentration). Anthropogenic NMHCs also contributed significantly (26 %) to the total VOC concentration and principally consisted of ethane and propane (which represented 34 % and 17 % of the anthropogenic NMHC cumulative concentration, respectively) as well as n-butane (7 %). The high contribution of species which generally have the longest lifetimes in the atmosphere (see Sect. 3.5) is consistent with the remote location of the Ersa site and is in agreement with Michoud et al. (2017). BVOCs contributed little to the total VOC concentration (9 % on annual average, 13 % in summer). They were mainly composed of isoprene and α-pinene (contributions of 44 % and 32 % to the BVOC mass, respectively). These compounds are among the major BVOCs in terms of emission intensity for the Mediterranean vegetation (Owen et al., 2001) and accounted for half of isoprenoid concentrations recorded during the intensive field campaign conducted at Ersa in summer 2013 (Debevec et al., 2018;Kalogridis, 2014). By contrast, a larger α-terpinene contribution was noticed during the summer 2013 intensive field campaign than the 2-year observation period. Note that speciated monoterpenes were measured differently during the summer 2013 field campaign by means of an automatic analyser (see Sect. S4 in the Supplement).

VOC variability
Monthly and interannual variations of primary (anthropogenic and biogenic) NMHCs along with OVOCs selected in this study (Sect. S1) are discussed in this section. Seasonal VOC concentration levels are indicated in Table 4. Note that the results of comparison between the VOC monitoring measurements investigated in this study and concurrent field campaign measurements performed during the summers 2012-2014 (presented in Sect. S4 of the Supplement) sup-ported the representativeness of the 2-year observation period with regard to its summer concentration levels.

Biogenic VOCs
Concentration variations of three selected BVOCs, isoprene, α-pinene and camphene, were analysed at different timescales (monthly/interannual variations; Fig. 4). These BVOCs exhibited high concentrations from June to August consistently with temperature and solar radiation variations (see Sect. 3.1). Indeed, throughout the summer 2013 observation period, Michoud et al. (2017) and Kalogridis (2014) observed that emissions of isoprene and the sum of monoterpenes were mainly governed by temperature and solar radiation, supported both by the diurnal variations in their concentrations (Geron et al., 2000a, b;Guenther et al., 2000) and their correlations with environmental parameters. Additionally, significant concentrations of α-pinene were noticed from September to November (Fig. 4), while isoprene concentrations were close to the detection limit and temperature and solar radiation were decreasing. However, solar radiation decreased much faster than temperature during these months (Fig. 3), which could suggest that additional emissions (Laothawornkitkul et al., 2009), dependent only on temperature contrary to those prevailing in summer, have influenced α-pinene concentrations during these months. Ozone concentrations, lower in autumn (i.e. the months of October to December) than in summer (O 3 concentration variations are depicted in Fig. S2 of the Supplement), also pointed out a weaker degradation of α-pinene in autumn.
Biogenic compounds showed significant interannual variations over the study period, linked to temperature and solar radiation variations. Higher June concentrations of isoprene and α-pinene were noticed in 2012 (average concentrations of 1.0 ± 1.1 and 2.6 ± 1.4 µg m −3 for isoprene and α-pinene, respectively) and 2014 (0.7 ± 0.5 and 0.2 µg m −3 ) than in 2013 (0.2 ± 0.2 and < 0.1 µg m −3 ). Higher June concentrations of camphene (and α-terpinene; not shown) were also noticed in 2014 than in 2013 (Fig. 4). These concentration levels may be related to the fact that temperature and solar radiation were more favourable to enhancing June biogenic emissions in 2012 and 2014 than in 2013 (Sect 3.1). Due to air relative humidity values observed in June (Sect. 3.1), we cannot rule out that an increase in BVOC concentrations may be linked to a transient modification of BVOC emissions induced by drought stress (Ferracci et al., 2020;Loreto and Schnitzler, 2010;Niinemets et al., 2004). Moreover, isoprene and α-pinene concentrations in July and August were higher in 2013 (average concentrations of 0.5 ± 0.3 and 1.1 ± 0.4 µg m −3 for isoprene and α-pinene, respectively) than in 2012 (0.3 ± 0.2 and 0.6 ± 0.3 µg m −3 ). High concentrations of camphene and α-terpinene were also noticed in August 2013 (0.2 ± 0.1 and 0.3 ± 0.3 µg m −3 , respectively; Fig. 4). Solar radiations were lower in July and August 2012, temperatures were slightly lower in July 2012 and mean wind speed was slightly higher in July 2012 (Fig. 3). These meteorological conditions may have affected biogenic emissions and have favoured their dispersion and their dilution by marine air masses, owing to the position of the Ersa station (Sect. 2.1). Note that the interpretation of interannual variations in BVOC measurements is based on a limited number of sampling days during the study period and different collection times (Table 1 and Sect. 2.2.1). It should therefore be considered cautiously, given variable day-to-day and strong diurnal BVOC concentration variations which were observed during the summer 2013 observation period (Kalogridis, 2014).

Anthropogenic VOCs
Variations of a selection of NMHCs are depicted in Fig. 5. These compounds illustrate contrasted reactivity, based on their atmospheric lifetimes (estimated according to their photochemical reaction rates with OH radicals defined in Atkinson, 1990, andArey, 2003). Despite lifetimes in the atmosphere ranging from a few hours to some days, all selected NMHCs were characterized by similar seasonal variations, with an increasing winter trend followed by a decrease in spring/summer ( Fig. 5 and Table 4), with the exception of n-hexane, propene and C 8 aromatics (which were the most reactive species of the NMHCs selected in this study and had the lowest concentrations; Tables 3 and 4). NMHC concentrations were higher in winter than in summer, up to

Species
Winter Spring Summer Autumn Anthropogenic Ethane 2.9 ± 0.5 1.8 ± 0.6 1.0 ± 0.2 1.9 ± 0.5 NMHCs Propane NMHCs exhibited different concentration levels during the two studied winter periods (Fig. 5). Mean winter NMHC concentrations were higher in 2013 than in 2014 by up to 0.3 µg m −3 in the cases of propane, acetylene and benzene (relative differences of 15 %, 42 % and 42 %, respectively). These compounds and ethane had the longest lifetimes among those selected in this study. However, ethane concentrations recorded at Ersa did not show any interannual variation over the study period (Fig. 5).

Oxygenated VOCs
Variations of selected OVOCs, illustrating contrasted reactivity (following the same methodology as that applied to NMHCs; Sect. 3.4.2), are presented in Fig. 6. Formaldehyde, acetaldehyde, glyoxal, methylglyoxal and C 6 -C 11 aldehydes have relatively short lifetimes into the atmosphere, and hence they can result from relatively close sources. On the other hand, acetone and methyl ethyl ketone (MEK) have the longest atmospheric lifetime among the OVOCs selected in this study, and hence they can also result from distant sources and/or be formed within emission-enriched air masses advected to the Ersa station.
Firstly, formaldehyde, methylglyoxal and n-hexanal have shown similar seasonal variations (Fig. 6), with high summer and spring concentrations (Table 4), suggesting an important contribution of primary/secondary biogenic sources to their concentrations. Fu et al. (2008) found that the largest global sources for methylglyoxal were isoprene and to a lesser extent acetone. Besides photochemical production, nhexanal and formaldehyde can be notably emitted by many plant species (Guenther et al., 2000;Kesselmeier and Staudt, 1999;Wildt et al., 2003). Interannual variations in formaldehyde, methylglyoxal and n-hexanal summer concentrations confirmed their links with biogenic sources. For instance, the highest concentrations of methylglyoxal were observed in June 2012 (average concentration of 0.7 µg m −3 ), simi- Acetaldehyde and acetone have shown an increase in their concentrations more marked in summer than in winter (Fig. 6), suggesting they were probably mainly of both secondary (anthropogenic/biogenic) and primary biogenic origins. Acetaldehyde is known to be mainly produced through the chemical transformation of anthropogenic and biogenic VOCs (Rottenberger et al., 2004;Schade and Goldstein, 2001;Seco et al., 2007;Wolfe et al., 2016), particularly in clean and remote areas. Acetaldehyde can also be released by plants (Jardine et al., 2008;Rottenberger et al., 2008;Winters et al., 2009). Acetone emissions are thought to be globally of biogenic rather than anthropogenic origin (Goldstein and Schade, 2000;Schade and Goldstein, 2006). Acetone can also result from the oxidation of various VOCs (Goldstein and Schade, 2000;Jacob et al., 2002;Singh et al., 2004), and roughly half of its concentrations measured at diverse urban or rural sites have been assigned to regional background pollution by several studies (e.g. Debevec et al., 2017;de Gouw et al., 2005;Legreid et al., 2007; with regional contributions at a scale of hundreds of kilometres). Additionally, acetaldehyde and acetone concentration variations in winter (e.g. mean February concentration in 2013 was 0.5 and 2.4 µg m −3 higher than in 2014, respectively) also pinpointed primary/secondary anthropogenic origins (Sect. 3.4.2).
Glyoxal and MEK showed an increase in their concentrations in both summer and winter ( Fig. 6 and Table 4), suggesting they were probably produced by several biogenic and anthropogenic sources. Glyoxal increases were in similar proportions (Fig. 6 and Table 4  indicate that primary/secondary anthropogenic sources primarily contributed to MEK concentrations. Fu et al. (2008) found that the largest global sources for glyoxal were isoprene and to a lesser extent acetylene and direct emissions. MEK can be emitted from both vegetation and numerous anthropogenic sources and can result from the oxidation of various VOCs (Yáñez-Serrano et al., 2016, and references therein). Glyoxal and MEK exhibited different concentration levels during the two studied winter periods since their mean February concentration in 2013 was 65 %-75 % higher than in 2014, confirming their links with anthropogenic sources. These two OVOCs also showed different interannual variations in late spring and summer (Fig. 6). The highest concentrations of glyoxal were measured in June 2012 (average concentration of 0.2 µg m −3 ), similarly to isoprene and α-pinene Finally, n-nonanal did not show clear seasonal variations. Its March and April concentrations were higher in 2014 than in 2013. An increase in MEK concentrations was also noticed in March and April 2014 in lower proportion. These findings suggest contributions from different sources in winter/early spring, compared to the rest of the year. Moreover, a slight increase in n-nonanal concentrations was noted in August 2013, which is consistent with a biogenic source contribution (Matsunaga et al., 2003;Wildt et al., 2003).

Major NMHC sources
In the coming section, major sources which have impacted primary NMHC concentrations measured at Ersa were iden-tified using a PMF five-factor solution (from simulations presented in Sect. S2). Figure 7 depicts factor contributions to the species chosen as inputs for the PMF tool along with NMHC contributions to the five factors defined by the factorial analysis. Figure 8 and Table 5 show factor contribution time series and their relative contributions to the total concentration of the NMHCs selected in this factorial analysis, respectively. Average factor relative contributions as a function of the measurement sampling time are also indicated in Table 5. In the present section, lifetimes were assessed from kinetic rate constants of the reactions of the selected NMHCs with OH (Atkinson and Arey, 2003) given an average OH concentration of 0.5 × 10 6 and 2.5 × 10 6 molecules cm −3 in winter and summer, respectively (Spivakovsky et al., 2000).
As VOC concentrations resulted from direct emissions, chemistry, transport and mixing, each individual computed factor cannot be attributed solely to one source category, especially for a remote receptor site like Ersa. A part of them may not be specifically associated with emission profiles but should rather be considered aged profiles originating from several source regions comprising several source categories (Sauvage et al., 2009). The PMF analysis was hence performed in this study to define co-variation factors of NMHCs that were characteristic of aged or local primary emissions along with secondary photochemical transformations taking place during the transport of air masses observed at this remote site (Michoud et al., 2017).
Since the low number of NMHCs considered in the factorial analysis of this study, PMF result relevance was checked, benefiting from a previous PMF analysis performed with the Ersa VOC summer 2013 dataset (Michoud et al., 2017; see Sect. 3.5.6). Experimental strategies to characterize VOC concentration variations at a remote background such as Ersa are also discussed in Sect. 3.5.6.

Local biogenic source (factor 1)
Factor 1 contributions to measured concentrations of the NMHCs selected in the factorial analysis were 0.2 ± 0.4 µg m −3 on average (corresponding to 4 % of the sum of these NMHCs - Table 5) and peaked at 3.1 µg m −3 on 20 June 2012. In late spring/summer, it was one of the main factors observed (16 % on average and up to 53 %; Fig. 8). The chemical profile of factor 1 depicts an elevated contribution of isoprene, known as a chemical marker for biogenic emissions, with its variability fully related to this factor. The isoprene relative contribution to factor 1 is 70 %. The estimated tropospheric lifetime of isoprene was quite short (winter: 6 h and summer: ∼ 1 h), indicating that this compound was emitted mostly by local vegetation. Consequently, factor 1 is labelled "local biogenic source". Note that the factor 1 contribution to selected NMHC concentrations may be slightly influenced by the two different sampling times used during the 2-year observation period (Table 5) as well as the number of VOCs and data points considered in the PMF analysis (see Sects. 3.4.1 and S2). However, Michoud et al. (2017) provided additional information on this local primary biogenic source in summer, like the contributions of additional primary BVOCs (the sum of monoterpenes) and some OVOCs (carboxylic acids, methanol and acetone), and presented its clear diurnal variations.

Short-lived anthropogenic sources (factor 2)
Factor 2 contributed 19 % to the sum of the selected NMHCs (Table 5; average contribution of 1.0 ± 0.5 µg m −3 ). This factor mainly consists of toluene (73 % of its variability attributed to this factor; Fig. 7), C 8 aromatic compounds (93 %), ethylene (48 %) and propene (83 %), typically emitted by combustion processes and with short-to-medium lifetimes (winter: 1-4 d; summer: 5-20 h). The cumulative relative contribution of these VOCs to factor 2 is 66 %. Besides road traffic, toluene is also a good marker for solvents generated by industrial sources (Buzcu and Fraser, 2006;Leuchner et al., 2015;Zhang et al., 2014), suggesting that this profile could also be imputed to industrial sources. Additionally, a significant proportion of C 5 -C 6 alkanes, i.e. i-pentane (32 % explained), n-pentane (37 %) and n-hexane (51 %), typically emitted by gasoline evaporation and with medium lifetimes (winter: 4-6 d; summer: 21-31 h), also contributes 19 % to this factor. Factor 2 did not show a good correlation with CO (Pearson correlation coefficient only of 0.2). As a conclusion, factor 2 is hence attributed to the grouping of several short-lived anthropogenic sources, partly related to gasoline combustion and/or evaporation and solvent use. Note that the PMF model did not reconstruct well 5 of the species comprising this factor (propene, n-pentane, n-hexane, toluene, and C 8 aromatics; Sect. S2 in the Supplement), especially their concentration variations over short periods. As a result, factor 2 contributions over short periods may be underestimated. However, given the objectives of this study, the examination of factor 2 contribution variations will be limited to seasonal and interannual ones (Sect. 4). Factor 2 contribution variations over a short period were previously investigated in Michoud et al. (2017).
Regarding factor 2 contributions coupled with air-mass origins (Fig. 9), more elevated contributions were noticed at Ersa under the influence of continental air masses coming from France and Europe. The distinction of short trajectories from long ones (see Sect. 3.2) highlighted that factor 2 was potentially influenced by relatively close sources when Ersa received air masses from continental France, whereas other continental European sources were probably more distant (Fig. 9). Furthermore, CF analysis applied to factor 2 contributions (see Fig. S3 of the Supplement) confirmed that this factor was influenced by various potential emission areas, either located in Italy (the Po Valley and central Italy), France (south-eastern region) or possibly in central Europe (western Hungary, Croatia and Slovenia). Given the CF results, ship emission contribution cannot be discarded as well, as already . Factor contributions to each species (µg m −3 ) and the percent of each species apportioned to the factor are displayed as a grey bar and a coloured circle, respectively. Factor 1 -local biogenic source; factor 2 -short-lived anthropogenic sources; factor 3 -evaporative sources; factor 4 -long-lived combustion sources; factor 5 -regional background. Table 5. Average relative factor contributions ± standard deviations (%) for the whole period and as a function of the measurement sampling time (see Table 1).
(with lifetimes of 10-11 d in winter and ∼ 2 d in summer), with more than 69 % of their variabilities explained by factor 3 (Fig. 7), along with i-/n-pentanes, n-hexane and propane (42 %-59 %; lifetimes of 4-21 d in winter and 21 h-4 d in summer). The cumulative relative contribution of these alkanes to factor 3 is up to 88 %. C 3 -C 6 alkanes are identified in the gasoline composition and evaporation sources (storage, extraction and distribution of gasoline or liquid petroleum gas; Sauvage et al., 2009, and references therein). Additionally, propane can be viewed as a relevant profile signature of natural gas transport, storage and use (Leuchner et al., 2015). Consequently, this factor can be viewed as "evaporative sources". Regarding factor contributions as a function of air-mass origins (Fig. 9), more elevated contributions were noticed at the Ersa station when aged air masses originated from France and Europe and were probably transported toward Ersa by relatively distant sources. According to the CF analysis ap-plied to factor 3 contributions (Fig. S3 of the Supplement), the Po Valley (especially Emilia-Romagna, an Italian region centre for food and automobile production), central Italy, the south-east of France and the Sardinian region seemed to be identified as the main potential emission areas for factor 3, as well as possibly more distant emission areas in central Europe such as western Hungary (i.e. the western Transdanubian region specializing in the automotive and machinery industries).

Long-lived combustion sources (factor 4)
The average contribution of factor 4 to the selected NMHC concentrations is evaluated at 0.9 ± 0.7 µg m −3 (corresponding to 16 % of their sum; Table 5). Its profile displays an important contribution from acetylene (100 % explained; Fig. 7), benzene (49 %) and propane (37 %), with lifetimes of 19-26 d in winter and 4-5 d in summer. The average cumulative relative contribution of these VOCs to factor 4 is 80 %. Figure 9. NMHC factor contributions (µg m −3 ) as a function of air-mass origin represented by box plots; the blue solid line, the red marker, and the box represent the median, the mean, and the interquartile range of the values, respectively. The bottom and top of the box depict the first and third quartiles and the ends of the whiskers correspond to the first and ninth deciles. Factor 2 -short-lived anthropogenic sources; factor 3 -evaporative sources; factor 4 -long-lived combustion sources; factor 5 -regional background. Air masses originating from France and Europe are subdivided into short and long trajectories to highlight local and more distant contributions (see Sect. 3.2). Note that the NMHC dataset used for the PMF analysis included different sampling time hours (09:00-13:00 or 12:00-16:00 UTC) following shifts during the 2-year period (see Table 1).
Aromatic compounds and acetylene are generally associated with combustion sources, such as vehicle exhaust (e.g. Badol et al., 2008;Pang et al., 2014). However, factor 4 is characterized by a loading of benzene far superior to toluene's (49 % and 3 %, respectively), suggesting it is more related to a residential heating source than a traffic one (Elbir et al., 2007;Leuchner et al., 2015;Sauvage et al., 2009). The factor 4 pro-file, mainly composed of long-lived species together with a low contribution of shorter-lived species, may indicate that it is related to partly aged air masses which were advected towards the Ersa station. These suggestions are consistent with the fact that this factor correlated particularly well with CO (Pearson correlation coefficient of 0.8). As a result, this fac-Atmos. Chem. Phys., 21, 1449Phys., 21, -1484Phys., 21, , 2021 https://doi.org/10.5194/acp-21-1449-2021 tor can be viewed as "long-lived combustion sources", including residential heating. Factor 4 showed higher contributions when the Ersa station received European air masses (see Fig. 9), particularly those with a long trajectory, which is consistent with factor 4 profile examination. The CF analysis depicted in Fig. S3 of the Supplement only pointed out western Hungary and to a lesser extent the Po Valley as the main potential emission areas for factor 4.

Regional background (factor 5)
Contributing 39 % to the selected NMHC concentrations (average contribution of 2.2 ± 1.0 µg m −3 ), factor 5 corresponds to the dominant source detected at Ersa during the study period. Its profile is principally composed of ethane (with its variability fully explained by factor 5; Fig. 7) and propane (18 % explained; lifetimes of 21-93 d in winter and 4-19 d in summer). Their relative contribution to factor 5 is up to 96 % and they typically result from the transport, storage and use of natural gas. Additional anthropogenic NMHCs with shorter lifetimes are attributed to this factor, including ethylene and propene (12 %-16 % explained; lifetimes of 21 h-3 d in winter and 4-13 h in summer), despite a low contribution to factor 5 mass (∼ 3 %). Given the high abundance of long-lived species in its profile, factor 5 can be viewed as a regional contribution of diverse remote sources of the Mediterranean region, thus indicating the continental regional background (Hellén et al., 2003;Leuchner et al., 2015;Sauvage et al., 2009; regional contributions at a scale of hundreds of kilometres). As a conclusion, we associated factor 5 with the "regional background".
Factor 5 showed slightly higher contributions when the Ersa station was under the influence of European air masses, particularly those with a long trajectory (potentially connected to distant emission areas; Fig. 9). As expected, the CF analysis applied to factor 5 contributions did not clearly pinpoint a specific potential emission area (Fig. S3 of the Supplement), apart from maybe western Hungary and to a lesser extent the Po Valley.

Towards the best experimental strategy to characterize variations in VOC concentrations observed at a remote background site
In this study, we decided to restrain the number of VOCs in the factorial analysis in order to have a better temporal representativeness of factor contributions (see Sect. S2). This choice has undoubtedly had some consequences for the PMF solution results of this study (in terms of factor number, chemical profiles and variations). Consequently, to support our results, we compared our five-factor PMF solution modelled with the 2-year NMHC dataset (14 variables, 152 4 h resolution observations), with the six-factor PMF solution modelled with the summer 2013 VOC dataset (42 variables, 329 1.5 h resolution observations; Michoud et al., 2017). The comparative study (see Sect. S5 in the Supplement) pointed out a good representativeness of the primary sources identified in this study. Overall, sources identified in both factorial analyses as influencing NMHC concentrations at Ersa had similar chemical compositions, regarding NMHCs in common. This finding supports the relevance of the 17 NMHCs selected in the factorial analysis of this study to identify major sources which have influenced Ersa NMHC concentrations over 2 years and suggests that they include key source tracers.
In order to support these statements, we computed a summer 2013 PMF solution only with the NMHCs selected in the factorial analysis of this study, and we compared it with the summer 2013 PMF solution modelled with a higher VOC number by Michoud et al. (2017). Results (see Sect. S6 in the Supplement) revealed a global consistency of the two summer 2013 PMF solutions in terms of factor chemical composition, regarding NMHCs in common between the two factorial analyses, and contribution variations, with the exception of those of short-lived anthropogenic sources (which may also depend on additional VOCs; Sect. S6). PMF factorial analyses computed with a limited VOC number therefore show limitations in the investigation of short-lived anthropogenic sources influencing Ersa VOC concentrations given the remote location of the Ersa station.
The two studies on Ersa VOC measurements also helped to discuss here the best experimental strategy to characterize variations in VOC concentrations measured at a remote background site surrounded by vegetation such as Ersa (Sect. 2.1). Firstly, the consideration of monoterpenes in the 2-year PMF analysis would significantly increase the contribution of primary biogenic sources (see Sect. 4.1). A better temporal resolution and coverage of the VOC measurements would have been beneficial in this study to further assess the contribution of local biogenic sources, considering BVOC day-today variations (as assessed with the summer 2013 dataset; Sect. 3.4.1). Additional measurements of speciated monoterpenes, sesquiterpenes and BVOC oxidation products would even help in identifying different local sources associated with various emitting Mediterranean plants (Owen et al., 2002). This statement is supported by the fact that different interannual variations in BVOC concentrations were noticed between isoprene and selected speciated monoterpenes in this study (Sect. 3.4.1). Contrarily, concentrations of nonspeciated monoterpenes and isoprene measured in summer 2013 were mainly attributed to the same primary biogenic source in Michoud et al. (2017). However, concentration variations of the sum of monoterpenes may have mirrored those related to abundant biogenic species such as α-pinene and thus masked those characterizing more reactive species such as α-terpinene (Kalogridis, 2014;Debevec et al., 2018).
The limited number of VOCs included in our PMF study had the advantage of a better deconvolution by the PMF model of primary sources and a better reconstruction of mea-sured VOCs with low concentrations (Sects. S2 and S6), without the need to minimize the weight of species with elevated background levels (unlike ethane, methanol and acetone in Michoud et al., 2017). These low concentrations thus represented a higher proportion of the total VOC concentrations that the model was supposed to explain. The longer timescale of VOC measurements (i.e. the 2-year period) even helped to deconvolve long-lived combustion sources from regional background (Sect. S5). However, PMF comparison results also raised the importance of a finer time resolution and a higher temporal coverage of VOC measurements conducted over a long period at a remote background site such as Ersa to support source apportionment results, in terms of decorrelation of anthropogenic sources. Indeed, anthropogenic sources identified in this study with bi-weekly observations have shown some significant correlations between them (Sect. S5) as a result of similar seasonal variations in concentrations of NMHCs composing them (Sects. 3.4.2 and 4). The consideration of diurnal variations in the 2-year factorial analysis would help in limiting this potential statistical bias.
The incorporation of 19 OVOCs in the source apportionment of Michoud et al. (2017) had little impact on the deconvolution of the main sources influencing NMHC concentrations observed at the Ersa station (Sects. S5 and S6). However, it can modify their relative contributions to the total VOC concentration, emphasizing those of local biogenic/anthropogenic sources and decreasing those of regional anthropogenic sources (Michoud et al., 2017). VOCs were largely dominated at Ersa by OVOCs (Sect. 3.3), especially in summer (Michoud et al., 2017). OVOC abundance observed at Ersa was further increased by the measurement of alcohols and carboxylic acids, carried out only in summer 2013 with automatic analysers (Sect. S4). The incorporation of OVOCs in the 2-year PMF analysis would have helped to better identify co-variations between these species (Sect. 3.4.3) and their proportion resulting from primary sources, but the limited number of samples realized with DNPH cartridges did not favour it ( Fig. S1 and Sect. S2). Formaldehyde has shown in this study high concentration levels (Table 3)

Discussions on the seasonal variability of VOC concentrations
In this section, seasonal accumulated concentrations of the 35 VOCs selected in this study (21 NMHCs and 14 OVOCs; Sect. S1) are examined (Fig. 10) so as to identify the prevailing drivers of their variations. Measured concentrations of the 17 NMHCs selected in the factorial analysis (Sect. 3.5) were distributed among their major sources in Fig. 10. Figure S4 presents separately variations in NMHC factor con-tributions. Based on the work of Michoud et al. (2017) and Kalogridis (2014), it has been assumed in this study that the concentrations of the four monoterpenes (Table 1), not included in the PMF analysis, were fully emitted by the local vegetation. Biogenic sources refer here to the grouping of the PMF biogenic source and speciated monoterpenes. The coming discussion has been structured in two steps: on the one hand the factors explaining the VOC concentration levels and their variations in spring and summer and on the other hand the factors controlling those in autumn and winter.

VOC concentration variations in spring and summer
Firstly, the 21 NMHCs selected in this study have shown low concentrations during spring and summer periods (average seasonal accumulated concentration of 4.6 ± 0.1 µg m −3 ; Fig. 10), while the 14 OVOCs exhibited high concentrations (8.0 ± 1.8 µg m −3 ). PBL height can be high in these seasons (von Engeln and Teixeira, 2013), favouring vertical dispersion. OVOC concentration levels in spring and summer are in agreement with a high photochemical decay of VOCs, which were rapidly depleted. Consequently, the most reactive VOCs composing emission profiles of the anthropogenic sources can have reacted and/or been dispersed before reaching the Ersa station. The clear seasonal cycles of these sources, with a decline in their contributions in spring and summer (Fig. S4), confirmed this finding. Regional background contributed in spring and summer 24 % to 53 % to total NMHC concentrations. As natural gas sources were attributed to the regional background (Sect. 3.5.5), a decrease in their emissions can presumably occur in the hot season, enhancing the decline in regional background contributions (Fig. S4). These regional background contributions also suggest that aged emissions advected by air masses to the Ersa station significantly influenced VOC concentrations observed during these seasons. This statement is also in agreement with the high acetone concentration levels in spring and summer (4.8 ± 1.5 µg m −3 ; Fig. 10). As a reminder, this carbonyl compound can result from distant sources and/or be formed within polluted air masses before reaching the Ersa station (Sect. 3.4.3). As shown in Fig. 10, biogenic sources contributed to total NMHC concentrations especially in summer (up to 41 %), in connection with temperature and solar radiation variations (Sect. 3.4.1). This finding is also supported by the high summer concentration levels of formaldehyde, methylglyoxal, acetaldehyde and acetone ( Fig. 10; Sect. 3.4.3). Short-lived anthropogenic sources explained in spring and summer 16 % to 24 % of the NMHC concentrations. Despite the high reactivity of the VOCs comprising these anthropogenic sources (Sect. 3.5.2), their contributions were not reduced as sharply as those of other anthropogenic sources (Fig. S4), which may be related to high influences of ship transport and relatively close potential emission areas (Sect. 3.5.2). Evapo- NHMCs selected in the factorial analysis were apportioned to the five modelled NMHC sources. NMHC seasonal measured concentrations which were not modelled by the PMF tool were lower than 0.09 µg m −3 and are not reported here. Winter: 1 January-31 March periodsspring: 1 April-30 June periods -summer: 1 July-30 September periods -autumn: 1 October-31 December periods. Note that the VOC dataset included different sampling time hours (09:00-13:00 or 12:00-16:00 UTC) following shifts during the 2-year period (see Table 1). rative and long-lived combustion sources only contributed 6 % to 16 % in spring and summer, as they likely originated from emission areas relatively distant from the Ersa station (Sects. 3.5.3 and 3.5.4). Low contributions of longlived combustion sources can also be partly explained by a lower source strength (especially for a residential heating source; Sect. 3.5.4).
Looking now at the interannual variations during the spring and summer periods, Fig. 10 highlights that total NMHC concentrations were in the same range during the two summer periods as well as the two spring ones (absolute difference was 0.2 µg m −3 between summers 2012 and 2013 and below 0.1 µg m −3 between springs 2013 and 2014). As depicted in Fig. S4, spring and summer contributions of evaporative sources, long-lived combustion sources, and regional background were in the same range over the study period, as well as MEK concentrations (0.3-0.4 µg m −3 ; Fig. 10), which have been attributed mostly to anthropogenic (primary/secondary) sources (Sect. 3.4.3). These findings can suggest that these anthropogenic sources, originating from distant emission areas, were largely influenced during these seasons by pollutant depletion and vertical/horizontal dispersion. Spring and summer contributions of short-lived anthropogenic sources appear to be more variable over the study period (0.7-1.1 µg m −3 ; Figs. 10 and S4). This finding suggests that these sources were largely influenced by origins of air masses, which hence advected to Ersa numerous emissions, potentially of variable strength and from various locations relatively close to Ersa. Furthermore, biogenic sources showed slightly higher summer contributions in 2013 than in 2012 (1.9 and 1.4 µg m −3 , respectively), especially in connection with concentrations of monoterpenes (Fig. 10) and consistent with temperature and solar radiation variations (Sect. 3.4.1). Summer concentrations of formaldehyde and acetone were also higher in 2013 than in 2012 (differences of 1.5-1.6 µg m −3 ), while acetaldehyde and methylglyoxal ones were in the same range (Fig. 10), in agreement with isoprene. Considering regional background contributions and MEK concentrations in summer periods, interannual variations in acetone summer concentrations seem to be related to a change in its biogenic primary/secondary contributions rather than a change in its background levels. These findings suggest enhanced emissions from the local vegetation in summer 2013, partly related to the release of monoterpenes, and/or higher biogenic secondary contributions. For similar reasons, spring contributions of biogenic sources were slightly higher in 2014 than in 2013 (Fig. 10). Transient drought stress of the vegetation could also have influenced BVOC emissions in spring 2014 (Sect. 3.4.1).

VOC concentration variations in autumn and winter
During autumn and winter periods, total NMHC concentrations increased (7.8 ± 1.2 µg m −3 ), while OVOC concentrations declined (7.1 ± 0.8 µg m −3 ). The decline in OVOC concentration levels mainly results from the decrease in acetone and formaldehyde concentrations (3.3 ± 0.8 µg m −3 and 1.0 ± 0.3 µg m −3 , contributing 34 %-52 % and 10 %-16 % to total OVOC concentrations, respectively), partly related to a weaker biogenic activity. Chemical reaction rates also dropped in autumn and winter as a consequence of the decrease in OH concentrations related to the decline in available UV light and temperatures (Fig. 3). As a result, VOCs were not removed from the atmosphere as quickly as in spring and summer, in agreement with OVOC concentration levels. PBL height also decreased during these seasons, impacting pollutant vertical dispersion. These conditions, favouring pollutant advection and accumulation, explained the higher contributions of anthropogenic sources perceived at Ersa in autumn and winter, especially those from sources mainly of regional origins (Fig. S4).
The regional background and evaporative sources were identified as the dominant contributors to total NMHC concentrations in autumn and winter (contributions of 31 %-45 % and 22 %-29 %, respectively). Long-lived combustion sources also contributed, particularly in winter (explaining 18 %-24 % of total NMHC concentrations), since the low ambient temperatures typical of winter (Fig. 3) may lead to an increased use of residential heating (Sect. 3.5.4). In order to better identify regional influences, average contributions of the anthropogenic sources in autumn and winter were investigated as a function of air-mass origin in Fig. 11. During these seasons, the Ersa station was mostly influenced by continental air masses coming from Europe and France (their cumulative occurrences were 36 % to 70 %). NMHC anthropogenic sources showed higher accumulated contributions when European air masses were advected to Ersa (Fig. 11), and they could be attributed to potential emission areas located in Italy and possibly more distant ones in central Europe (Sect. 3.5). To a lesser extent, high accumulated contributions of anthropogenic sources were also noticed when Ersa received air masses originating from continental France (potential emission areas located in the south-east of France) and Corsica-Sardinia. Figure 10 highlights notable interannual variations in VOC concentrations in winter (absolute difference of 3.5 µg m −3 ), while they were in the same range in autumn (difference below 0.1 µg m −3 ). Indeed, NMHC winter concentrations were higher in 2013 (9.4 µg m −3 ) than in 2014 by up to 1.9 µg m −3 , similarly to OVOC winter concentrations (7.5 and 6.0 µg m −3 in 2013 and 2014, respectively). Even though winter contributions of long-lived combustion sources, shortlived anthropogenic sources and evaporative sources were significantly reduced in 2014 compared to in 2013 (absolute differences from 0.3 to 0.9 µg m −3 ; Fig. 10), the seasonal patterns of their variations were similar in 2013 and 2014, as depicted in Fig. S4. These findings could be evidence of homogenous regional background pollution distribution at a synoptic scale. Regional background winter contributions observed in 2013 and 2014 (absolute difference below 0.1 µg m −3 ; Fig. 10) are also in agreement with this suggestion. Acetone and MEK winter concentrations were lower in 2014 than in 2013 (Fig. 10), which could be related to variations in their anthropogenic primary/secondary sources rather than in their background levels.
The interannual variations in VOC winter concentrations could be partly explained by occurrences and origins of air masses advected to the Ersa station, their enrichment in different anthropogenic sources (based on potential emission areas observed in winters 2013 and 2014; Fig. 12) and their respective contributions (Fig. 11). On the one hand, Ersa was more under the influence of European air masses in winter 2013 than in winter 2014 (occurrences of 37 % and 18 %, respectively). When continental European and French air masses were advected to Ersa, anthropogenic sources showed higher accumulated winter contributions in 2013 (10.8 and 9.6 µg m −3 , respectively; Fig. 11) than in 2014 (8.2 and 7.8 µg m −3 ). In winter 2013, the main potential emission areas for the anthropogenic NMHC sources were located in Italy (Tuscan coasts and the Po Valley), central Europe (Slovenia and western Hungary) and to a lesser extent the south of France. Contrariwise, in winter 2014, Ersa did not seem to be influenced at all by air masses originating from central Europe (especially from Slovenia and Hungary; Fig. 12), and potential emission areas located in the Po Valley and the south-east of France did not contribute significantly to VOC concentrations. On the other hand, the Ersa station received air masses originating from Corsica-Sardinia, i.e. relatively close potential emission areas, more frequently in winter 2014 than in winter 2013 (occurrences of 24 % and 0 %, respectively). The cumulative contribution of anthropogenic sources in winter 2014 under Corsican-Sardinian air-mass influence was only 7.0 µg m −3 (Fig. 12). These findings can suggest less aged (but rather low) emissions have influenced Ersa concentrations in winter 2014, partly explaining why OVOC winter concentrations observed at Ersa were lower in 2014 than in 2013.
Moreover, most countries of western Europe experienced different winters in 2013 and 2014 (see Sect. 3.1 and Fig. S5 of the Supplement). As a reminder, winter 2013 was considered rather rigorous, since e.g. last access: 11 October 2020). On the other hand, winter 2014 was rather mild and temperatures were the hottest of the 1951-2014 period (the mean temperature for Europe was ∼ 2.6 • C higher than the pre-industrial period; EEA, 2015). As a consequence, the interannual variations in VOC winter concentrations may also result from different influences of meteorological conditions. Indeed, meteorological winter conditions were probably more favourable to inducing higher OH concentrations in 2014 than in 2013, leading to higher photochemical decay and thus lower NMHC concentrations observed at the Ersa station. Nevertheless, OVOC winter concentrations were higher in 2013 than in 2014, which can therefore highlight that they were mainly due to increased contributions from anthropogenic sources of regional origins, thus offsetting

VOC concentration variations in continental Europe
From June 2012 to June 2014, NMHC measurements were concurrently conducted at 17 other European background monitoring stations (described in Sect. S7), allowing us (i) to examine the representativeness of the Ersa station in terms of seasonal variations in NMHC concentrations impacting continental Europe and (ii) to provide some insights into dominant drivers for VOC concentration variations in Europe built on what we have learned from Ersa's VOC observations. Figure 13 depicts monthly concentration time series of a selection of NMHCs measured at the 18 European monitoring stations (including Ersa). Concentrations of NMHCs with typically medium-to-long lifetimes in the atmosphere, Figure 12. Potential source area contributions to NMHC anthropogenic sources during winters 2013 and 2014. Contributions are expressed in µg m −3 . Factor 2 -short-lived anthropogenic sources; factor 3 -evaporative sources; factor 4 -long-lived combustion sources; factor 5regional background. Winter: 1 January-31 March periods.
i.e. ethane, propane, n-butane, acetylene and benzene, were examined here since they can be significantly influenced by regional contributions, as observed in this study (Sect. 3.5) and elsewhere (Debevec et al., 2017;Michoud et al., 2017;Sauvage et al., 2009). As a reminder, concentrations observed at Ersa were mainly explained in this study by regional background for ethane, long-lived combustion sources for propane, acetylene and benzene, and evaporative sources for n-butane (see Sect. 3.5 and Fig. 7). The study of concentration variations of these source tracers can thus help to highlight temporal and spatial variations in source contributions to NMHC concentrations observed in most of continen-tal Europe. In addition, despite its shorter lifetime compared to other NMHCs of the selection, Fig. 13 also presents ethylene concentration variations in order to investigate shortlived anthropogenic source importance and variability in continental Europe.
Globally, concentrations of these selected NMHCs measured at the 17 European stations showed the same seasonal variations as observed at Ersa (Fig. 13  ever the location and the typology of the station. The highest absolute differences between NMHC summer concentrations measured at two stations was 0.4-0.7 µg m −3 for ethane, 0.1-0.2 µg m −3 for acetylene, 0.1-0.7 µg m −3 for propane and benzene, 0.2-0.6 µg m −3 for n-butane and 0.3-1.2 µg m −3 for ethylene. They suggest that the temperature was the main driver in regulating summer concentration levels, linked to photochemistry processes and the vertical dispersion (Sect. 4.1). Ethane concentration levels were still relatively important during summer (above 1.0 µg m −3 ), suggesting that long-range transport (up to intercontinental pollution transport) was among the main parameters governing VOC summer concentrations in continental Europe. On the other hand, anthropogenic NMHC monthly concentration levels in winter appear to be more spatially variable in continental Europe. The highest absolute differences between NMHC winter concentrations measured at two stations were 1.3-2.6 µg m −3 for ethane, 0.6-1.6 µg m −3 for propane, acetylene and benzene, 0.6-1.4 µg m −3 for nbutane, and 1.2-4.1 µg m −3 for ethylene. They may probably highlight spatial variations in anthropogenic source contributions to NMHC winter concentrations observed in continen- tal Europe. The lowest NMHC winter concentrations were observed at stations located in southern and western Europe, including Ersa, other French sites and high-altitude ones (see Fig. 13   However, NMHC winter concentration levels observed at the Ersa and north-western French stations, along with their interannual variations, were more similar to those in central Europe than those in south-western France and southern Europe (Fig. 14). Indeed, NMHC winter concentrations measured in central Europe were generally significantly higher in 2013 than in 2014 (total differences of 2.6-3.2 µg m −3 ), with the exception of stations located in north-eastern Germany (i.e. WAL, NGL and ZGT stations -Sect. S7; total differences of −2.8-0.7 µg m −3 ). These interannual variations seem to be particularly linked to contribution variations of long-lived combustion sources (reduction of 0.7-1.1 µg m −3 , i.e. 21 %-44 %) and short-lived anthropogenic sources (reduction of 0.6-1.3 µg m −3 , i.e. of 21 %-43 %). NMHC winter concentrations related to evaporative sources also decreased in 2014, but to a lesser extent (reduction of 0.3-0.5 µg m −3 , i.e. of 13 %-24 %). These findings are consistent with winter variations in anthropogenic source contributions impacting VOC concentrations at Ersa (Sect. 4.2), especially when the stations received air masses originating from central Europe (Fig. 11). They also highlight that interannual variations in NMHC winter concentrations observed in central Europe were influenced by both regional and local contributions. Furthermore, due to synoptic phenomena that occurred in winters 2013 and 2014 (Sect. 3.1), temperatures in central Europe were, respectively, colder and warmer compared to average values for the 1981-2010 period (Fig. S5 of the Supplement). Interannual variations in winter temperatures were more marked in central Europe than in southern France and southern Europe (Fig. S5 of the Supplement), which could partly explain their respective interannual variations. Moreover, NMHC winter concentrations monitored in Scandinavia (represented by PAL station results in Fig. 14) were higher in 2014 than in 2013 (total difference of 0.7 µg m −3 ) as well as those measured at stations located in north-eastern Germany (especially NGL and ZGT stations; differences of 1.7-2.8 µg m −3 ). These increases in winter 2014 in these regions concerned in particular NMHCs related to long-lived combustion (0.2-0.5 µg m −3 , i.e. 9 %-25 %) and evaporative sources (0.3-0.4 µg m −3 , i.e. 11 %-18 %). Even though these regions experienced a cold winter in 2013 (Fig. S5 of the Supplement), early winter 2014 in northern Europe was also colder than average values for the 1981-2010 period, due to an intense cold wave resulting from a strong anticyclone centred on western Russia and extending from Finland to Crimea. Additionally, winter precipitation in these regions was less frequent and/or intense than average values for the 1981-2010 period (Fig. S5 of the Supplement) in 2014 than in 2013, which may have favoured accumulation and transport and hence induced higher NMHC concentrations measured at nearby monitoring stations.

Conclusions
Within the framework of the ChArMEx project, a monitoring station has been set up and operated from June 2012 to June 2014 at a remote background site of Corsica (Ersa) in the north-western part of the Mediterranean in order to examine seasonal variabilities of 35 selected VOCs (21 NMHCs and 14 carbonyl compounds).
Biogenic VOCs were principally imputed to the local vegetation and exhibited high concentrations and interannual variations in late spring and summer, related to temperature and solar radiation variations. Anthropogenic NMHCs, with atmospheric lifetimes ranging from a few hours to some days, have shown low concentrations in spring and summer since the Ersa station was mostly influenced by aged advected air masses. Anthropogenic sources identified in this study were hence largely influenced by the pollutant depletion and vertical dispersion before reaching the Ersa station, especially those originating from distant emission areas. As a result, the regional background and short-lived anthropogenic sources (from various locations relatively close to Ersa) mainly contributed to NMHC concentrations in spring and summer. Anthropogenic NMHCs were also characterized by high concentration levels in autumn and winter and winter concentrations higher in 2013 than in 2014. The weaker pollutant depletion and vertical dispersion in autumn and winter have favoured the advection to Ersa of air masses enriched in primary anthropogenic sources originat- ing from more distant emission areas (located in France, Italy and to a lesser extent central Europe) than in spring and summer. As a result, the regional background, evaporative sources and long-lived combustion sources were identified as the dominant contributors to NMHC autumn and winter concentrations. Interannual variations in their winter contributions were mainly governed by occurrences and origins of air masses advected to the Ersa station as well as by meteorological conditions (temperature and precipitation) occurring before they reached the Ersa station. OVOCs largely prevailed in the VOC abundance during the 2-year monitoring period. They have also shown different seasonal and interannual variations, suggesting several source contributions. OVOC concentrations measured at Ersa were generally higher in summer, which could be the result of a high contribution of secondary anthropogenic and primary/secondary biogenic sources, whereas their concentrations during autumn and winter declined and were potentially more influenced by primary/secondary anthropogenic sources, more specifically for MEK and glyoxal. OVOCs showed higher winter concentrations in 2013 than in 2014 as a result of the enrichment of air masses advected to Ersa in primary anthropogenic sources of regional origins.
Selected NMHC variabilities in continental Europe supported Ersa being considered a good regional representative station. However, NMHC winter concentration levels can significantly vary temporally and spatially, pointing out spatial variations in anthropogenic source contributions to NMHC concentrations observed in continental Europe in relation to meteorological conditions. Ersa winter concentration variations were found to be more representative of central Europe than southern/western or northern Europe. These findings point out the interest in conducting multi-site and multi-year measurements that are sufficiently representative of interannual and spatial variations in regional pollution impacting continental Europe in winter. They also pointed out that the PMF-derived factors controlling VOC concentration variations at remote sites may be mainly controlled by the meteorological conditions that occurred during the study period when issued from short observation periods (i.e. up to 2 months).
Concurrent datasets of VOC concentrations from three summer campaigns performed at Ersa helped to support the representativeness of the 2-year monitoring period in terms of summer concentration levels, variations, and source apportionment. They also highlight limitations of this study in the investigation of seasonal variations in Ersa VOC concentrations, in connection with the VOC range and the temporal resolution of their measurements considered in the PMF analysis. After this work, some questions remain regarding the contribution of local BVOC and OVOC sources at seasonal and interannual scales. It would be interesting to conduct additional long-term VOC measurements at Ersa, including OVOCs and tracers of various primary sources, at a finer time resolution and a higher temporal coverage. This would help to complete the understanding of determinants governing VOC concentration variations initiated by both Michoud et al. (2017) and this study.
Author contributions. SS, NL, JS and FD designed the research and were involved in the logistics and the collection of VOC samples in the field. TS calculated uncertainties related to VOC measurements conducted with DNPH cartridges and steel canisters and validated them following the ACTRIS protocol. CD, SS, VG and NL analysed VOC data, conducted and interpreted the VOC PMF analysis and examined dominant factors controlling VOC concentrations. CD wrote the manuscript. All the co-authors were involved in data discussion and edited the paper. son, Vincent Gaudion, Laurence Depelchin and Isabelle Fronval for their contributions to obtaining VOC datasets, from the organization and the setting of instruments at Ersa to the analysis of VOC samples at the laboratory, but also Thierry Bourianne and the Qualitair Corse team for their involvement on site concerning VOC sample routines as well as the associated logistic procedure. The authors finally thank the ChArMEx project manager Eric Hamonou for his logistical help.
Financial support. This study received financial support from the MISTRALS/ChArMEx programme funded by CNRS/INSU, CEA, and Météo-France, from ADEME (grant nos. 1262c0009, 1362c0031, and 1462c0031), the French Environmental Ministry, the CaPPA projects, and the Communauté Territoriale de Corse (CORSiCA project). This research has also been supported by the CaPPA project (Chemical and Physical Properties of the Atmosphere) funded by the French National Research Agency (ANR) through the PIA (Programme d'Investissement d'Avenir, grant no. ANR-11-LABX-0005-01) and by the Regional Council Nord-Pas de Calais and the European Funds for Regional Economic Development (FEDER). This research was also funded by the European Union Seventh Framework Programme under the DEFI-VOC project (grant agreement no. 293897) and by CARBO-SOR/Primequal (grant no. 10-MRES-PRIMEQUAL-5-CVS-007no. 2100082986) and SAF-MED (grant no. ANR-12-BS06-0013-02).
Review statement. This paper was edited by Xavier Querol and reviewed by two anonymous referees.