Pan-European rural monitoring network shows dominance of NH 3 gas and NH 4 NO 3 aerosol in inorganic atmospheric pollution load

. A comprehensive European dataset on monthly atmospheric NH 3 , acid gases (HNO 3 , SO 2 , HCl), and aerosols (NH + 4 , NO − 3 , SO 2 − 4 , Cl − , Na + , Ca 2 + , Mg 2 + ) is presented and analysed. Speciated measurements were made with a low-volume denuder and ﬁlter pack method (DEnuder for Long-Term Atmospheric sampling, DELTA ® ) as part of the EU NitroEurope (NEU) integrated project. Altogether, there were 64 sites in 20 countries (2006–2010), coordinated between seven European laboratories. Bulk wet-deposition measurements were carried out at 16 co-located sites (2008– 2010). Inter-comparisons of chemical analysis and DELTA ® measurements allowed an assessment of comparability between laboratories. The form northern across 4


Introduction
Air quality policies and research on atmospheric sulfur (S) and nitrogen (N) pollutant impacts on ecosystems and human health have focused on the emissions, concentrations, and depositions of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ), ammonia (NH 3 ), and their secondary inorganic aerosols (SIAs; ammonium sulfate, (NH 4 ) 2 SO 4 ; ammonium nitrate, NH 4 NO 3 ) (ROTAP, 2012;EMEP, 2019). The aerosols, formed through neutralization reactions between the alkaline NH 3 gas and acids generated in the atmosphere by the oxidation of SO 2 and NO x (Huntzicker et al., 1980;AQEG, 2012), are a major component of fine particulate matter (PM 2.5 ) (AQEG, 2012;Vieno et al., 2016a) and precipitation (ROTAP, 2012;EMEP, 2019). The negative effects of these pollutants on sensitive ecosystems are mainly through acidification (excess acidity) and eutrophication (excess nutrient N) processes that can lead to a loss of key species and decline in biodiversity (e.g. Hallsworth et al., 2010;Stevens et al., 2010). They are also implicated in radiative forcing and influence climate change through inputs of nitrogen that can alter the carbon cycle (Reis et al., 2012;Sutton et al., 2013;Zaehle and Dalmonech, 2011).
A number of EU policy measures (e.g. 2008/50/EC Ambient Air Quality Directive: EU, 2008;/2284/EU National Emissions Ceilings Directive, NECD: EU, 2016 and wider international agreements (e.g. Gothenburg Protocol; UNECE, 2012) are targeted at abating the emissions and environmental impacts of SO 2 , NO x , and NH 3 . The largest emissions reductions have been achieved for SO 2 , which decreased by 82 % across the 33 member countries of the European Environment Agency (EEA-33) since 1990, to 4743 kt SO 2 in 2017 (EEA, 2019). Reductions in NO x emissions have been more modest, at 45 % over the same period, with emissions in 2017 of 8563 kt NO x exceeding those of SO 2 . By contrast, the reductions in NH 3 emissions (of which over 90 % come from agriculture) have been more modest, decreasing by only 18 %. Here, the decrease was largely driven by reductions in fertilizer use and livestock numbers, in particular from eastern European countries, rather than through implementation of any abatement or mitigation measures. More worryingly, the decreasing trend has reversed in recent years, with emissions increasing by 5 % since 2010, to 4788 kt NH 3 in 2017 (EEA, 2019).
In recent assessments, critical loads of acidity were exceeded in about 5 % of the ecosystem area across Europe in 2017 (EMEP, 2018). While the substantial decline in SO 2 emissions has allowed the recovery of ecosystems from acid rain, NH 3 from agriculture and NO x from transport are increasingly contributing to a larger fraction of the acidity load. Although NH 3 is not an acid gas, nitrification of NH 3 and ammonium (NH + 4 ) releases hydrogen ions (H + ) that acidify soils and fresh water. The deposition of reactive N (N r ; including oxidized N: NO x , HNO 3 , NO − 3 , and reduced N: NH 3 , NH + 4 ) and its contribution to eutrophication effects have also been identified by the European Environment Agency (EEA) as the most important impact of air pollutants on ecosystems and biodiversity (EEA, 2019). The deposition of N r throughout Europe remains substantially larger than the level needed to protect ecosystems, with critical-load thresholds for eutrophication from N exceeded in around 62 % of 877 the EU-28 ecosystem area and in almost all countries in Europe in 2017 (EMEP, 2018).
Following emission, atmospheric transport and fate of the gases are controlled by the following processes: short-range dispersion and deposition, chemical reaction and formation of NH + 4 aerosols, and the long-range transport and deposition of the aerosols ROTAP, 2012). Atmospheric S and N r inputs from the atmosphere to the biosphere occur though (i) dry deposition of gases and aerosols, (ii) wet deposition in rain, and (iii) occult deposition in fog and cloud (Smith et al., 2000;ROTAP, 2012). The deposition processes contribute very different fractions of the total S or N r input and different chemical forms of the pollutants at different spatial scales. NH 3 is a highly reactive, water-soluble gas and deposits much faster than NO x (which is not very water-soluble and has low deposition velocity). Dry N deposition by NH 3 therefore contributes a significant fraction of the total N deposition to receptors close to source areas and will often exert the larger ecological impacts compared with other N pollutants (Cape et al., 2004;Sutton et al., 1998Sutton et al., , 2007. Numerous studies have shown that N r deposition in the vicinity of NH 3 sources is dominated by dry NH 3 -N deposition (e.g. Pitcairn et al., 1998;Sheppard et al., 2011), with removal of NH 3 close to a source controlled by physical, chemical, and ecophysiological processes (Flechard et al., 2011;Sutton et al., 2007Sutton et al., , 2013. Unlike NO x , HNO 3 (from oxidation of NO x ) is very water-soluble, while NO − 3 particles can act as cloud condensation nuclei (CCN) so that they are both scavenged quickly and removed efficiently by precipitation. Since NO x is inefficiently removed by precipitation, wet deposition of NO x near a source is small and only becomes important after NO x has been converted to HNO 3 and NO − 3 . Because of the large numbers of atmospheric N species and their complex atmospheric chemistry, quantifying the deposition of N r is hugely complex and is a key source of uncertainty for ecosystem effect assessment (Bobbink et al., 2010;Fowler and Reis, 2007;Schrader et al., 2018;Sutton et al., 2007). Input by dry deposition can be estimated using a combination of measured and/or modelled concentration fields with high-resolution inferential models (e.g. Smith et al., 2000;Flechard et al., 2011) or by making direct flux measurements (e.g. Fowler et al., 2001;Nemitz et al., 2008). Although it is possible to measure N r deposition directly (e.g. Skiba et al., 2009), the flux measurement techniques are complex and resource-intensive, unsuited to routine measurements at a large number of sites. The "inferential" modelling approach provides a direct estimation of deposition from N r measurements by applying a land-use-dependent deposition velocity (V d ) to measured concentrations (Dore et al., 2015;Flechard et al., 2011;Simpson et al., 2006;Smith et al., 2000).
At present, there are limited atmospheric measurements that speciate the gas-and aerosol-phase components at multiple sites over several years. On a European scale, atmo-spheric measurements of sulfur (SO 2 , particulate SO 2− 4 ) and nitrogen (NH 3 , HNO 3 , particulate NH + 4 , NO − 3 ) have been made by a daily filter pack method across the European Monitoring and Evaluation Program (EMEP) networks since 1985, providing data for evaluating wet-and dry-deposition models (EMEP, 2016;Tørseth et al., 2012). The method, however, does not distinguish between the gas-and aerosolphase N species. Consequently, these data are reported as total inorganic ammonium (TIA = sum of NH 3 and NH + 4 ) and total inorganic nitrate (TIN = sum HNO 3 and NO − 3 ), limiting the usefulness of the data. Speciated measurements by an expensive and labour-intensive daily annular denuder method are also made (Tørseth et al., 2012) but are necessarily restricted to a small number of sites due to the high costs associated with this type of measurement. There are also networks with a focus on specific N components, for example, the national NH 3 monitoring networks in the Netherlands (LML; van Zanten et al., 2017) and in the UK (National Ammonia Monitoring Network, NAMN; Tang et al., 2018a) or compliance monitoring across Europe in the case of SO 2 and NO x . The UK is unique in having an extensive set of speciated gas and aerosol monitoring data from the Acid Gas and Aerosol Network (AGANet), with measurements from 1999 to the present (Tang et al., 2018b).
In this context, there is an ongoing need for cost-effective, easy-to-operate, time-integrated atmospheric measurement for the respective gas and aerosol phases at sufficient spatial scales. Such data would help to (1) improve estimates of N deposition; (2) contribute to development and validation of long-range transport models, e.g. EMEP (Simpson et al., 2006) and EMEP4UK (Vieno et al., , 2016b; (3) interpret interactions between the gas and aerosol phases; and (4) interpret ecological responses to nitrogen (e.g. ecosystem biodiversity or net carbon exchange). To contribute to this goal, a "three-level" measurement strategy in the EU Framework Programme 6 Integrated Project "NitroEurope" (NEU; http://www.nitroeurope.ceh.ac.uk/, last access: 29 July 2020) between 2006 and 2010 delivered a comprehensive integrated assessment of the nitrogen cycle, budgets, and fluxes for a range of European terrestrial ecosystems Skiba et al., 2009). At the most intensive level (Level 3), state-of-the-art instrumentation for high-resolution, continuous measurements at just 13 "flux super sites" provided detailed understanding on atmospheric and chemical processes . By contrast, manual methods with a low temporal frequency (monthly) at the basic level (Level 1) provided measurements of N r components at a large number of sites (> 50 sites) in a cost-efficient way in a pan-European network . Key species of interest included NH 3 , HNO 3 , and ammonium aerosols ((NH 4 ) 2 SO 4 , NH 4 NO 3 ).
In this paper, we present and discuss 4 years of monthly reactive gas (NH 3 , HNO 3 , HCl) and aerosol (NH + 4 , NO − 3 , SO 2− 878 Y. S. Tang et al.: Pan-European rural monitoring network shows dominance of NH 3 gas Level 1 network set up under the NEU integrated project (Fig. 2). A harmonized measurement approach with a simple, cost-efficient time-integrated method, applied with high spatial coverage, allowed a comprehensive assessment across Europe. The gas and aerosol network was complemented by 2 years of wet-deposition data made at a subset of the sites (Fig. 3). The intention of the smaller bulk wet-deposition network was two-fold: (i) to provide wet-deposition estimates at DELTA ® (DEnuder for Long-Term Atmospheric sampling) sites that do not already have such measurements on site and (ii) to compare the relative importance of reduced and oxidized N versus sulfur in the atmospheric pollution load. Measurements across the network were coordinated between multiple European laboratories. The measurement approach and the operations of the networks, including the implementation of annual inter-comparisons to assess comparability between the laboratories, are described. The data are discussed in terms of spatial and temporal variation in concentrations, relative contribution of the inorganic nitrogen and sulfur components to the inorganic pollution load, and changes in atmospheric concentrations of acid gases and their interactions with NH 3 gas and NH + 4 aerosol.

NEU Level 1 DELTA ® network
The NitroEurope (NEU) Level 1 network was operated between November 2006 and December 2010 to deliver the core measurements of reactive nitrogen gases (NH 3 , HNO 3 ) and aerosols (NH + 4 , NO − 3 ) for the project (Fig. 1). A lowvolume denuder filter pack method, the "DEnuder for Long-Term Atmospheric sampling" system (DELTA ® ; Sutton et al., 2001a;Tang et al., 2009Tang et al., , 2018b, with time-integrated monthly sampling was used, which made implementation at a large number of sites possible. Other acid gases (SO 2 , HCl) and aerosols (SO 2− 4 , Cl − , Na + , Ca 2+ , Mg 2+ ) were also collected at the same time and measured by the DELTA ® method. DELTA ® measurements were co-located with all NEU Level 3 sites with advanced flux measurements  and with the network of main CarboEurope-IP CO 2 flux monitoring sites (http://www.carboeurope.org/ ceip/, last access: 5 January 2020) (Flechard et al., 2011. Two of the UK sites in the NEU DELTA ® network are existing UK NAMN (Tang et al., 2018a) and AGANet sites (Tang et al., 2018b). These are Auchencorth Moss (UK-Amo) and Bush (UK-EBu), located in southern Scotland. Monthly gas and aerosol data at the two sites, made as part of the UK national networks, were included in the NEU network. NEU network N r data were used, together with a range of dry-deposition models, to model dry-deposition fluxes (Flechard et al., 2011) and to assess the influence of N r on the C cycle, potential C sequestration, and the greenhouse gas balance of ecosystems using CO 2 exchange data from the Reaction scheme for the formation of ammonium aerosols from interaction of NH 3 with acid gases HNO 3 , SO 2 , and HCl, showing the components (green) that were measured in the Ni-troEurope (NEU) DELTA ® network. Dry deposition of the gas and aerosol components was estimated by inferential modelling (Flechard et al., 2011), while wet deposition (blue) was measured in the NEU bulk wet-deposition network at a subset of the DELTA ® sites.
co-located CarboEurope sites . Other measurements made at the Level 1 sites included estimation of wet-deposition fluxes (Sect. 2.3) and also soil and plant bioassays (Schaufler et al., 2010).
Altogether, the DELTA ® network covered a wide distribution of sites across 20 countries and four major ecosystem types: crops, grassland, semi-natural, and forests. These sites can be described as "rural" and were chosen to provide a regionally representative estimate of air composition. The network site map is shown in Fig. 2, with site details given in Table S1 in the Supplement. Further information on the sites are also provided in Flechard et al. (2011). Network establishment started in November 2006, with 57 sites operational from March 2007 onwards. Over the course of the network, some sites closed or were relocated due to infrastructure changes, and new sites were also added. A total of 64 sites provided measurements at the end of the project, with 45 of the sites operational the entire time. In addition, replicated DELTA ® measurements were made at four sites: Ecosystem types are C (crops), G (grassland), F (forests), and SN (short semi-natural; includes moorland, peatland, shrubland, and unimproved and upland grassland). Replicated (P: parallel) DELTA measurements are made at four sites: SK04/SK04P; UK-AMo/UK-AMoP (NH 3 and NH + 4 only), UK-Bu/UK-BuP, and FR-Fgs/FR-FgsP (NaCl-coated denuders instead of K 2 CO 3 -glycerol in sample train).

Coordinating laboratories
A team of seven European laboratories shared responsibility for running the network. Measurement was on a monthly timescale, with each laboratory preparing and analysing the monthly samples with documented analytical methods (see Table S3 for information on analytical methods and limit of detection (LOD)) for between 5 and 16 DELTA sites (Fig. 2). The use of a harmonized DELTA ® methodology, coupled to defined quality protocols , ensured comparability of data between the laboratories (see later in Sect. 3.1 and 3.2). A network of local site operators representing the science teams of each site performed the monthly sample changes and posted the exposed samples back to their designated laboratories for analysis. Air concentration data were submitted by the laboratories for their respective sites in a standard reporting template to UKCEH. Following data checks against defined quality protocols , the finalized dataset was uploaded to the NEU database (http://www.nitroeurope.ceh.ac.uk/, last access: 29 July 2020). Establishment of the network, including the first year of measurement results on N r com-880 Y. S. Tang et al.: Pan-European rural monitoring network shows dominance of NH 3 gas ponents, is reported in Tang et al. (2009). Information on co-located measurements and agricultural activities at each of the sites was also collected and is accessible from the NEU website (http://www.nitroeurope.ceh.ac.uk/, last access: 29 July 2020).

DELTA ® methodology
The DELTA ® method used in the NEU Level 1 network is based on the system developed for the UK Acid Gas and Aerosol Monitoring Network (AGANet, Tang et al., 2018b). Full details of the DELTA ® method and air concentration calculations in the NEU network are provided by Tang et al. (2009Tang et al. ( , 2018b. The method uses a small 6 V air pump to deliver low air-sampling rates of between 0.2 and 0.4 L min −1 , a high-sensitivity gas meter to record the typical monthly volume of air collected, and a DELTA ® denuder filter pack sample train to collect separately the gasand aerosol-phase components. The sample train is made up of two pairs of base-and acid-impregnated denuders (15 and 10 cm long) to collect acid gases and NH 3 , respectively, under laminar conditions. A two-stage filter pack with baseand acid-coated cellulose filters collects the aerosol components downstream of the denuders. The base coating used was K 2 CO 3 -glycerol, which is effective for the simultaneous collection of HNO 3 , SO 2 , and HCl (Ferm, 1986), while the acid coating was either citric acid for temperate climates or phosphorous acid for Mediterranean climates (Allegrini et al., 1987;Ferm, 1979;Perrino et al., 1990;Fitz, 2002). In this way, artefacts between gas-and aerosol-phase concentrations are minimized (Ferm et al., 1979;Sutton et al., 2001a). The DELTA ® air inlet has a particle cut-off of ∼ 4.5 µm, which means fine-mode aerosols in the PM 2.5 fraction and some of the coarse-mode aerosols < PM 4.5 will be collected (Tang et al., 2015).
A low-voltage version of the AGANet DELTA ® system was built centrally by UKCEH and sent to each of the European sites, where they were installed by local site contacts. These systems operated on either 6 V (off mains power with a transformer) or 12 V from batteries (wind-and solarpowered). Air sampling was direct from the atmosphere without any inlet lines or filters to avoid potential loss of components -in particular HNO 3 , which is very "sticky" -to surfaces. Sampling height was 1.5 m above ground or vegetation in open areas. In forested areas, the DELTA ® equipment was set up either in large clearings or on towers 2-3 m above the canopy (see Flechard et al., 2011).

Calculation of gas and aerosol concentrations
Atmospheric gas and aerosol concentrations in the DELTA ® method are calculated from the number of inorganic ions (NH + 4 , NO − 3 , SO 2− 4 , Cl − , and base cations) in the denuder and aerosol aqueous extracts and the volume of air sampled (from gas meter readings), which is typically 15 m 3 for a monthly sample. The volumes of deionized water used to extract acid-coated denuders and aerosols filters are 3 and 4 mL, respectively. For the base-coated denuders and aerosol filters, the extract volume in both cases is 5 mL An example is shown here for calculating the atmospheric concentrations of NH 3 (gas) (Eq. 1) and NH + 4 (aerosol) (Eq. 2) from the aqueous extracts, based on an air volume of 15 m 3 collected in a typical month.
Pairs of base-and acid-coated denuders are used to collect the acid gases and alkaline NH 3 gas, respectively. This allows denuder collection efficiency of, for example, NH 3 (Eq. 3) to be assessed as part of the data quality assessment process. An imperfect acid coating on the denuders for example can lead to lower capture efficiencies (Sutton et al., 2001a;Tang and Sutton, 2003).
Denuder collection efficiency, NH 3 (%) A correction, based on the collection efficiency, is applied to provide a corrected air concentration (χ a (corrected); Eq. 4) (Sutton et al., 2001a;Tang et al., 2018aTang et al., , 2018b. With a collection efficiency of 95 %, the correction amounts to 0.3 % of the corrected air concentration. For an efficiency below 60 %, the correction amounts to more than 50 % and is not applied. The air concentration (χ a ) of NH 3 is then determined as the sum of NH 3 in denuders 1 and 2 (Tang et al., 2018a). By applying the infinite series correction, the assumption is that any NH 3 (and other gases) that is not captured by the denuders will be collected on the downstream aerosol filter. To avoid double-counting, the estimated amount of "NH 3 breakthrough" is subtracted from the NH + 4 concentrations on the aerosol filter. Sea salt SO 2− 4 (ss-SO 2− 4 ) in aerosol was estimated according to Eq. (5), based on the ratio of the mass concentrations of

Artefact in HNO 3 determination
Results from the first DELTA ® inter-comparison in the NEU network ) (see also Sect. 2.5) and further work by Tang et al. (2015Tang et al. ( , 2018b have shown that HNO 3 concentrations may be overestimated on the carbonatecoated denuders used due to co-collection of other oxidized nitrogen components, most likely from nitrous acid (HONO). In the UK AGANet, HNO 3 data are corrected with an empirical factor of 0.45 derived by Tang et al. (2015). Since the correction factor for HNO 3 is uncertain (estimated to be ±30 %) and derived for UK conditions, no attempt has been made to correct the HNO 3 data from the NEU network. The DELTA ® method remained unchanged throughout the entire network operation and provided a consistent set of measurements by the same protocol. The caveat is that the HNO 3 data presented in this paper also include an unknown fraction of oxidized N, most probably HONO, and therefore represent an upper limit in the determination of HNO 3 . Contribution from NO 2 is likely to be small since this is collected with a low efficiency on carbonate-coated denuders (Bai et al., 2003;Tang et al., 2015), and the network sites are rural, where NO x concentrations are expected to be in the low parts per billion. At the French Fougéres parallel site (FR-FgsP), NaCl-coated denuders were used to measure HNO 3 to compare with results from K 2 CO 3 -glycerol-coated denuders at the main site (FR-Fgs) (see Sect. 2.1 for methodology and Sect. 3.3.1 for data inter-comparison).

NEU bulk wet-deposition network
The NEU bulk wet-deposition network (Fig. 3, Table S2) was established to provide wet-deposition data on NH + 4 and NO − 3 . It was set up 2 years after the establishment of the NEU DELTA ® network, with sites located at a subset of DELTA ® sites that did not already have on-site wetdeposition measurements. Sampling commenced at some sites in January 2008, with 14 sites operational from March 2008. Site changes also occurred during the operation of this network, again with some site closures and new site additions over time. In total, 12 sites provided 2 years of monthly data, with a further 6 sites providing 1 year of monthly data between 2008 and October 2010, when measurements ended.
The type of bulk precipitation collector used was a Rotenkamp sampler (Dämmgen et al., 2005), mounted 1.5 m above ground, or in the case of forest sites, either in clearings or above the canopy. Each unit has two collectors providing replicated samples, comprising a pyrex glass funnel (aperture area = 84.9 cm 2 ) with vertical sides, connected directly to a 3 L collection bottle (material: low-density polyethylene), which was changed monthly. Thymol (5-methyl-2-(1methylethyl)phenol) (150 mg) was added as a biocide  to a clean, dry pre-weighed bottle at the start of each collection period. This provided a minimum thymol concentration of 50 mg L −1 for a full bottle to preserve the sample against biological degradation of labile nitrogen compounds during the month-long sampling.
Three European laboratories shared management and chemical analysis for the network (Fig. 3). The laboratories were CEAM (all three Spanish sites); INRAE (French Renon site); and SHMU, designated the main laboratory responsible for all other sites. A full suite of precipitation chemistry analyses were carried out that included pH, conductivity, NH + 4 , NO − 3 , SO 2− 4 , PO 3− 4 , Cl − , Na + , K + , Ca 2+ , and Mg 2+ . Rain volumes and precipitation chemistry data were submitted in a standard template to UKCEH for checking and then uploaded to the NEU database (http://www.nitroeurope.ceh.ac.uk/, last access: 29 July 2020). Samples with high P (> 1 µg L −1 PO 3− 4 ), high K + , and/or NH + 4 values that are indicative of bird contamination were rejected. Annual wet deposition (e.g. kg N ha −1 yr −1 ) was estimated from the product of the 882 Y. S. Tang et al.: Pan-European rural monitoring network shows dominance of NH 3 gas species concentrations and rain volume. Determinations of organic N were also carried out on some of the rain samples in a separate investigation reported by Cape et al. (2012).

Laboratory inter-comparisons: chemical analysis
All laboratories in the DELTA ® and bulk wet-deposition networks participated in water chemistry proficiency testing (PT) schemes in their own countries as well as the EMEP (once annually; http://www.emep.int, last access: 8 February 2010) and/or WMO-GAW (twice annually; http://www. qasac-americas.org, last access: 12 January 2020) laboratory inter-comparison schemes. PT samples for analysis are synthetic precipitation samples for determination of pH, conductivity, and all the major inorganic ions at trace levels. In addition, UKCEH also organized an annual PT scheme for the duration of the project (NEU-PT) to compare laboratory performance in the analysis of inorganic ions at higher concentrations relevant for DELTA ® measurements. This comprised the distribution of reference solutions containing known concentrations of ions that were analysed by the laboratories as part of their routine analytical procedures.

Field inter-comparisons: DELTA measurements
Prior to the NEU DELTA ® network establishment, a workshop was held to provide training to participating laboratories on sample preparation and analysis. This was followed by a 4-month inter-comparison exercise (July to October 2006) between six laboratories at four test sites (Montelibretti, Italy; Braunschweig, Germany; Paterna, Spain, and Auchencorth, UK). Results of the inter-comparison on N r components were reported by Tang et al. (2009), which demonstrated good agreement under contrasting climatic conditions and atmospheric concentrations of the N r gases and aerosols. The first DELTA ® inter-comparison allowed the new laboratories to gain experience in making measurements and was an extremely useful exercise to check how the whole system works, starting with coating of denuders and filters and DELTA ® train preparation, sample exchange via post, sample handling, and inter-comparing laboratory analytical performance. Further DELTA ® inter-comparisons between laboratories were conducted each year for the duration of the project, details of which are summarized in Table 1. At each test site, DELTA ® systems were randomly assigned to each of the participating laboratories. All laboratories provided DELTA ® sampling trains for each of the inter-comparison sites and carried out chemical analysis on the returned exposed samples. Measurement results were returned in a standard template to UKCEH, the central coordinating laboratory for collation and analysis.

European emissions data
With the exception of Russia and Ukraine, official reported national emissions data on SO 2 , NO x , and NH 3 are avail-able for all other 18 countries in the NEU network from the European Environment Agency (EEA) website (EEA, 2020). Emissions data for the period 2007 to 2010 were extracted, and the emission densities of each gas (t km −2 yr −1 ) in each country were derived by dividing the 4-year-averaged total emissions by the land area (km 2 ). Gridded emissions data (at 0.1 • × 0.1 • resolution) for SO 2 , NO x , and NH 3 are available from the EMEP emissions database (EMEP, 2020). The 0.1 • × 0.1 • gridded data for the period 2007 to 2010 were downloaded and were used to estimate national total emissions (sum of all grid squares in each country) and 4-yearaveraged emission densities (t km −2 yr −1 ) for Russia and Ukraine. As a check, total emissions for the other 18 countries were also calculated by this method and were the same as the national emission totals reported by the EEA (EEA, 2019).  , 2015). The focus of the MAN network is to provide site-based NH 3 concentrations for the nature conservation sites rather than a representative spatial-concentration field for the country. Hourly NH 3 and SO 2 data which were also available from the eight sites in the LML network were downloaded from the Rijksinstituut voor Volksgezondheid en Milieu (RIVM, the Dutch National Institute for Public Health and the Environment) website (http://www.lml.rivm.nl/gevalideerd/index.php, last access: 6 November 2018). The 4-year-averaged NH 3 and SO 2 concentrations for the period 2007 to 2010 were calculated and used to complement measurement data from the four Dutch sites in the NEU DELTA ® network.

UK NAMN and AGANet network data
Atmospheric NH 3 , acid gases, and aerosols are measured in the UK NAMN (since 1996) and AGANet (since 1999) (Tang et al., 2018a, b). The UK approach is a high-density network with low-time-resolution (monthly) measurements, combining an implementation of the DELTA ® method used in the present NEU DELTA ® network and a passive ALPHA ® method  to increase network coverage in NH 3 measurements (Sutton et al., 2001b;Tang et al., 2018a). Monthly and annual data for the overlapping period of the project were extracted from the UK-AIR website (https://uk-air.defra.gov.uk/, last access: 25 November 2019) and nested with the NEU network data for analysis in this paper.
3 Results and discussion 3.1 Laboratory inter-comparison results: chemical analysis Figure 4 compares the percentage deviation of results from reference solution concentrations ("true value") reported by the laboratories for different chemical components in the EMEP, WMO-GAW, and NEU proficiency testing (PT) schemes, combined from 2006 to 2010. Each data point is colour-coded in the graphs according to the laboratory providing the measurements. Altogether, results from the combined PT schemes produced > 100 observations for each reported chemical component over the 4-year period. The performances of laboratories in Fig. 4 can be summarized in terms of the percentage of reported results agreeing within 10 % of the true values (see summary table below Fig. 4), where the true values represent the nominal concentrations in the aqueous test solutions. The best agreements were for SO 2− 4 and NO − 3 , with an average of 92 % and 87 % of all reported results agreeing within 10 % of the true value across the concentration range covered in the PT schemes. In the case of NH + 4 , while an average of 90 % of reported results were within 10 % of the reference at 1 mg L −1 NH + 4 , laboratory performance was poorer (68 % agreeing within 10 %) at lower concentrations (0.1-0.9 mg L −1 ). Poorer performance at the low concentrations was largely due to two laboratories (CEAM and SHMU), with > 50 % of their results reading high. For Na + and Cl − , the percentages of results agreeing within 10 % of the refer-ence were 81 % and 86 %, respectively, across the full range of PT concentrations. At concentrations above 1 mg L −1 , the agreement improved and increased to 89 % for Na + and 96 % for Cl − . A larger spread around the reference values was provided for the base cations Ca 2+ and Mg 2+ at low concentrations (< 1 mg L −1 ). The percentage of results passing at low concentrations below 1 mg L −1 was 36 % (Ca 2+ ) and 59 % (Mg 2+ ), increasing to 80 % (Ca 2+ ) and 90 % (Mg 2+ ) above 1 mg L −1 . The larger scatter at low concentrations is likely due to uncertainty in the chemical analysis at or close to the method limit of detection and reflects challenges of measuring base cations, in particular Ca 2+ as this is very "sticky" and adsorbs-desorbs from surfaces, leading to analytical artefacts.
To show what the PT reference solution concentrations would correspond to if they were a denuder and/or aerosol extract, equivalent gas (Eq. 1) and/or aerosol concentrations (Eq. 2) (Sect. 2.2.1) are calculated for each of the ions and provided in the summary table in Fig. 4. A 0.5 mg L −1 NH + 4 solution, for example, is equivalent to an atmospheric concentration of 0.09 µg NH 3 m −3 (gas) or 0.13 µg NH + 4 m −3 (aerosol) for a monthly sample. In Fig. 5, scatterplots are shown comparing all NEU laboratory-reported results with PT reference, where all ion concentrations (mg L −1 ) from Fig. 4 have been converted to equivalent gas and aerosol concentrations (µg m −3 ), based on a typical volume of 15 m 3 over a month. With the exception of a small number of outliers, most data points are close to the 1 : 1 line, with laboratory results agreeing within ±0.05 µg m −3 in equivalent gas and/or aerosol concentrations. These are low ambient concentrations and show that the measurement uncertainty in the analysis of very low concentrations in the PT schemes will be small for the majority of sites in the network, where concentrations were found to be much higher (see Fig. 6).

Field inter-comparison results: DELTA ® measurements
Results from 4 years of annual DELTA ® field intercomparisons (2006)(2007)(2008)(2009), for all field sites, are combined and summarized in Fig. 6. The gas and aerosol concentrations measured and reported by each of the laboratories are compared with the median estimate of all laboratories in each of the scatterplots, with the colour of the symbols identifying the laboratory providing the measurements. Regression results (slope and R 2 ) in the table below the plots provide the main features of the inter-comparison. The slope is equivalent to the mean ratio of each laboratory against the median value, where values close to unity indicate closer agreement to the median value. Overall, the scatterplots show good agreement between the laboratories, with some laboratories showing very close agreement to the median estimates and more scatter observed from the others. The occurrence of outliers in some of the individual monthly values indicates that caution needs to be exercised in the interpretation of these data points in the intercomparison. To average out the influence of a few individ-ual outliers, the mean concentrations from each of the seven laboratories for each of the four field sites were calculated and compared with averaged median estimates of all laboratories for each site. A summary of the mean concentrations and the percentage difference from the median is presented in Table 2. Since the INRAE laboratory did not join the NEU network until 2008, averaged median values from the 2008 and 2009 inter-comparisons are used to compare with the INRAE results, included in the table for clarity. The mean concentrations between laboratories are broadly comparable. Each of the laboratories were also able to resolve the main differences in mean concentrations at the four field sites, ranging from the lowest concentrations at Auchencorth (e.g. median = 1.4 µg NH 3 m −3 ) to higher concentrations, representing a more polluted site at Paterna (e.g. median = 5.2 µg NH 3 m −3 ) for the test periods (Table 2). Larger differences for HCl, Ca 2+ , and Mg 2+ are due to clear outliers from one or two laboratories at the very low concentrations of these species encountered and may be related to measurement uncertainties at the low air concentrations. The compa- rability between laboratories for each of the components is next considered in turn.

Inter-comparisons: NH
The best agreement between laboratories was for the N r gases (NH 3 , HNO 3 ) and aerosol species (NH + 4 , NO − 3 ), with slopes within ±10 % of the median values and R 2 > 0.9 in the regression analysis from five of the laboratories (Fig. 6, Table 2). This is important since N r species were the primary focus for the NEU DELTA ® network. Slightly poorer agreement for NH 3 and NH + 4 was provided by CEAM and MHSC laboratories, with data points both above and below the 1 : 1 line (Fig. 6). The outliers above the 1 : 1 line from MHSC were from the 2006 inter-comparison exercise. Removal of these 2006 outliers improved the MHSC regression slope for NH 3 from 1.21 (R 2 = 0.87, n = 41) to 0.99 (R 2 = 0.99, n = 10) (Fig. S1). While this seems to suggest that the performance of MHSC for NH 3 improved following the first intercomparison exercise, the regression slope for aerosol NH + 4 increased instead from a slope of 1.26 (R 2 = 0.83, n = 41) to 1.48 (R 2 = 0.93, n = 10), suggesting an overestimation of NH + 4 concentrations (Fig. S1). A possible cause may be the quality and/or variability in the aerosol filter blank values for NH + 4 as laboratory blanks are subtracted from exposed samples to estimate aerosol NH + 4 concentrations. While the laboratory blanks reported by MHSC for aerosol NH + 4 were low (mean = 0.48 µg NH + 4 ) and smaller than other laboratories (mean = 0.64-1.20 µg NH + 4 ) ( Fig. S2), their field blanks in the 2006 DELTA inter-comparison exercise were on average 5.5 times larger than the laboratory blanks. This is likely due to extensive delays in getting samples released from customs in Slovakia at the start of the network. Another possibility is a breakthrough of NH 3 from the acid-coated denuders onto the aerosol filters. The denuder collection efficiency of NH 3 gas (Eq. 3, Sect. 2.2.1) reported by MHSC was on average 88 % for all years and 91 % where 2006 data have been excluded (Table S3). This is comparable with the mean collection efficiencies of all laboratories (91 % and 90 %) (Table S4), which makes NH 3 breakthrough an unlikely explanation for the higher readings. The assessment of NH + 4 is however more uncertain from the reduced number of data points (n = 10).
For the CEAM laboratory, reported NH 3 concentrations were on average 16 % lower (n = 41) than the median, with a slope of 0.89 (R 2 = 0.87), and particulate NH + 4 was on average 13 % lower (n = 41) than the median, with a slope of 0.42 (R 2 = 0.22) (Fig. 6). A need to improve the NH + 4 analysis (indophenol colorimetric assay) in the acid-coated denuders and aerosol filters by the CEAM laboratory was identified from the 2006 inter-comparison ). The Table 2. Inter-comparison of results from seven European laboratories at four different field test sites for all years (2006)(2007)(2008)(2009)(2010). The results shown are the mean concentrations from each laboratory for each site and the averaged median estimates derived from all laboratories for each site. indophenol method for aqueous NH + 4 determination is pHsensitive. Calibration solutions and quality control checks for the colorimetric assays are made up in deionized water (pH 7), whereas the aqueous extracts from the DELTA ® acidcoated denuders and cellulose filters are acidic (pH ∼ 3). Determination of NH + 4 in the denuder extracts may therefore be underestimated if the pH of the indophenol reaction has not been adjusted for the increased acidity in the sample extracts. When the 2006 data are excluded from the regression analysis, the slopes for NH 3 and NH + 4 increased to 1.02 (R 2 = 0.94, n = 12) and 0.98 (R 2 = 0.51, n = 12), respectively (Fig. S1). The improved agreement with other lab-oratories after the 2006 inter-comparison suggests that the method under-read was largely resolved, reflected in an improvement in the slope. Despite some uncertainties in the NH 3 and NH + 4 measurements, the laboratories were able to clearly resolve the main differences in mean concentrations at the four different field sites in all years ( Table 2). The results presented here for CEAM and MHSC highlight the importance of the initial inter-comparison exercise in identifying and resolving sampling and analytical issues at the start of the project.
3.2.2 Inter-comparisons: SO 2 , SO 2− 4 Six laboratories provided slopes within 12 % of the median values in the regression analysis for SO 2 (Fig. 6). The smaller R 2 values were from two laboratories (CEAM and SHMU; R 2 < 0.7), with data points both above and below the 1 : 1 line. For INRAE, the larger slope of 1.6 (R 2 = 9) was due to a single high SO 2 reading reported for Auchencorth of 2.0 µg SO 2 m −3 , compared with the median of 1.4 µg SO 2 m −3 . When the mean SO 2 concentrations measured by INRAE are compared with the median, the difference was on average 13 %, providing acceptable agreement, which suggests that the high reading may just be an outlier. There was more scatter in the inter-comparison for SO 2− 4 , although the majority of points are still close to the 1 : 1 line (Fig. 6). Six laboratories provided slopes within 12 % of the median values in the regression analysis also for SO 2− 4 . The regression slope from CEAM for SO 2− 4 was 1.2 (R 2 = 0.9), which is still within 20 % of the median. The SO 2 and SO 2− 4 measurements were broadly comparable between the laboratories, with mean concentrations agreeing on average within 6 % of the median (Table 2).

Inter-comparisons: HCl, Cl −
The HCl inter-comparison shows clear outliers from the CEAM laboratory, with concentrations that were on average up to 2 times higher than other laboratories (slope = 1.8). For example, a mean concentration of 1.8 µg HCl m −3 was reported by CEAM for Paterna, compared with a median of 0.7 µg HCl m −3 . Apart from CEAM, the mean concentrations of HCl reported by the other laboratories were generally comparable ( Table 2). The larger percentage differences between the measured mean and median at each site reflect the challenges of measuring the very low concentrations of HCl at these sites of < 0.5 µg HCl m −3 (slightly higher at Paterna). HCl results were reported by NILU for the 2008 inter-comparison exercise only, limiting the number of measurements (n = 4) available for comparison.
The comparison for Cl − showed better agreement of the CEAM laboratory results with other laboratories in both the inter-comparison of individual monthly values (Fig. 6) and the mean concentrations (Table 2). Like HCl, larger percentage differences between the measured concentrations and median at each site may be attributed to higher measurement uncertainties at the low concentrations of Cl − . For NILU, there were only two data points for Cl − from the Auchencorth site in the 2008 inter-comparison. Overall, the intercomparison for HCl and Cl − showed that the laboratories were able to resolve the main differences in mean concentrations at the different sites even at the low concentrations encountered.

Inter-comparisons: base cations (Na
Measurements of Ca 2+ and Mg 2+ were the most uncertain, with the largest scatter in the inter-comparisons (Fig. 6). Despite the trace levels of these base cations at all field sites, four laboratories (INRAE, UKCEH, SHMU, vTI) provided data close to the 1 : 1 line, demonstrating close agreement between these laboratories. The clear outliers above the 1 : 1 line are from CEAM, MHSC, and NILU, with slopes > 2. While MHSC over-read Ca 2+ and Mg 2+ , its results for Na + were in better agreement with other laboratories, with a slope of 0.9 (R 2 = 0.5) (Fig. 6). There was a lot of scatter in the data however, with outlier points both above and below the 1 : 1 line, suggesting measurement uncertainties in their base cation measurements. For NILU, the only base cation results reported by the laboratory were for the 2008 DELTA ® inter-comparisons at Auchencorth and Braunschweig. This accounts for the low number of data points (n = 4) from the NILU laboratory. The median concentrations of Ca 2+ and Mg 2+ at both field sites were very low (< 0.1 µg m −3 ), which makes comparison with the few data reported from NILU highly uncertain. Like NILU, CEAM also did not report base cation results for all of the DELTA ® inter-comparison. Base cation results provided by CEAM were for 2007-2009 only.

Comparisons according to ecosystem types
Annual averaged concentrations of gases and aerosols measured in the NEU DELTA ® network are presented in Fig. 7, with sites grouped according to each of the four major ecosystem types: crops, grassland, forests, and semi-natural. These are the classifications used in dry-deposition models, where ecosystem-specific deposition velocities (V d ) are combined with measurement data to produce estimates of N r dry deposition (Flechard et al., 2011). A total of 64 sites from 20 different countries, including replicated measurements at four of the sites, are compared in Fig. 7. Not all of the sites however were operational all of the time or at the same time. Changes in the numbers and locations of sites occurred over the duration of the network, for example, due to site closures, relocations, and/or new site additions. The annual averaged concentrations plotted for each  (2007)(2008)(2009)(2010) of sites in the NEU DELTA ® network, grouped according to ecosystem type: crops (n = 10), grassland (n = 9 + 1 parallel), semi-natural (n = 11 + 1 parallel), and forests (n = 34 + 2 parallel). (b) Percentage composition of gas and aerosol components measured at NEU DELTA ® network sites (n = 64 + 4 parallel sites) (mean of all annual mean concentrations from 2007 to 2010). Years with < 7 months of data, including 2006, are excluded. Where the number of years contributing to the annual average is < 4, the number is shown in brackets beside the site data. Ca 2+ and Mg 2+ data are not included as these were mostly at or below the limit of detection. Replicated DELTA measurements are made at four sites: FR-Fgs/FR-FgsP (NaCl instead of K 2 CO 3 -glycerol-coated denuders; HCl not measured), SK04/SK04P, UK-Ebu/UK-EBuP, and UK-AMo/UK-AMoP (NH 3 and NH + 4 only).
site are the mean of all available annual means. Where the annual averaged concentration is derived from less than 4 full years of data, the number of years providing the mean is shown, in brackets, next to the site data in the graph. To avoid bias in the calculation of annual means, due to seasonality in the data (see later in Sect. 3.5), years with incomplete data coverage (< 7 months of data in any year) were excluded. Applying these data exclusions, the number of sites that provided annual data was 55 for 2007, 57 for 2008, 54 for 2009, and 55 for 2010. The number of sites that provided annual data for each year over the entire period was 45 sites. Sites with parallel (P) DELTA ® measurements were Auchencorth Moss (UK-AMoP), Easter Bush (UK-EBuP), Fougéres (FR-FgsP), and SK04P (EMEP site in Slovakia) (Fig. 7). Overall, good reproducibility in DELTA ® measurements was demonstrated by the parallel measurements (Figs. S3-S6). At the Auchencorth Moss parallel site (UK-AMoP), NH 3 and NH + for these two components was on average within 4 % at the low concentrations measured at this site (annual mean: 0.5-0.9 µg NH 3 m −3 and 0.3-0.5 µg NH + 4 m −3 ) (Table S5). Parallel measurements at Easter Bush (UK-EBuP) stopped in March 2010. With the exception of Ca 2+ and Mg 2+ , the comparison of annual mean data from the replicated measurements for 2006 to 2009 provided excellent agreement of 4 % (NO − 3 ) to 12 % (NH − 3 ) at Easter Bush (Table S6). At Fougéres (Table S7), HNO 3 concentration measured on K 2 CO 3 -glycerol-coated denuders (FR-Fgs) was about 2-fold higher than on NaCl-coated denuders in the parallel DELTA ® system (FR-FgsP), consistent with overestimation of HNO 3 (on average 45 %) on carbonate-coated denuders (see Sect. 2.2.3). The disadvantage of a NaCl coating, however, is that it can only collect HNO 3 and not the other acid gases. A third carbonate denuder is necessary in the sample train to collect and measure SO 2 since SO 2 is only partially captured, and HCl cannot be measured on NaCl denuders (Tang et al., , 2018b. This explains the smaller SO 2 concentrations reported by the FR-FgsP site, with breakthrough of SO 2 (inefficiently captured by NaCl denuders) onto the aerosol filters resulting in larger particulate SO 2− 4 concentrations than the Fr-Fgs site. For the SK04 site, measurement reproducibility for the 4 years of parallel data for the N and S components was good, with agreement ranging from 1.2 % (NH + 4 ) to 9 % (SO 2− 4 ) ( Table S8). HCl and Na + and determinations were however more uncertain, with differences of 67 % and 43 %, respectively (Table S8). It has to be noted, however, that the concentrations of the two components were very low, at < 0.2 µg HCl m −3 and < 0.4 µg Na + m −3 . The differences in concentrations are therefore actually within ±0.1 µg m −3 for HCl and within ±0.2 µg m −3 for Na + .
A key feature in Fig. 7 is the dominance of N over S species at most sites, when expressed as micrograms per cubic metre of the element. The mean percentage contribution of sum N r (NH 3 -N, HNO 3 -N, NH + 4 -N, NO − 3 -N) concentrations to the total mass of gas and aerosol species measured is 52 % (range = 24 %-80 %), twice as much as from sum S (SO 2 -S and SO 2− 4 -S; mean = 23 %, range = 7 %-53 %) (Fig. 8). This is consistent with more substantial reductions in SO 2 emissions (−72 %) than achieved with NO x (−43 %) or NH 3 (−18 %) in Europe between 1991-2010 (EEA, 2019). The differences in atmospheric composition of S and N species in the present assessment therefore reflected changes in emissions of the precursor gases and are also in agreement with a recent assessment of air quality trends showing important changes in S and N composition in air and rain across the EMEP networks (EMEP, 2016).
Secondary NH + 4 particles are mainly in the "fine" mode with diameters of less than 2.5 µm (PM 2.5 ) and estimated to contribute between 10 % and 50 % of ambient PM 2.5 mass concentration in some parts of Europe (Putaud et al., 2010, Schwartz et al., 2016. An assessment by Hendriks et al. (2013) found that secondary NH + 4 contributed 10 %-20 % of the PM 2.5 mass in densely populated areas in Europe and even higher contributions in areas with intensive livestock farming. Concentrations of PM 2.5 continue to exceed the EU limit values of 25 µg m −3 annual mean in large parts of Europe in 2017 (EEA, 2019). Particulate NH + 4 data presented from the DELTA ® network therefore highlight the potential contribution of NH 3 of agricultural origin to fine NH + 4 aerosols in PM 2.5 . The formation and transport of these secondary aerosols pose a serious risk to human health since PM 2.5 is linked to increased mortality from respiratory and cardiopulmonary diseases (AQEG, 2012).
A considerable fraction of the aerosol components measured was made up of sea salt (Na + and Cl − ), with contributions from the sum of Na + and Cl − ranging from 4 % of the total aerosol loading at the inland Höglwald site in Germany (DE-Hog) to 43 % at Dripsey (IE-Dri), a coastal site in Ireland (Fig. 7). With the reduction in European emissions and concentrations of the gases SO 2 , NO x , and NH 3 for formation of NH + 4 -containing aerosols, sea salt is therefore assuming a proportionate increase in the aerosol composition, consistent with observations from a recent European assessment of composition and trends in long-term EMEP measurements (EMEP, 2016). The concentrations of Ca 2+ and Mg 2+ were very low across the network, with values (mean of all sites < 0.1 µg m −3 ) that were at or below method limit of detection (LOD = ∼ 0.1 µg m −3 ) (Table S3). These data are also considered to be underestimated due to the DELTA particle sampling cut-off (∼PM 4.5 ), and they were excluded from further assessment in this paper.

Comparisons with national gas emissions
In Fig. 9, the annual averaged gas and aerosol concentrations of grouped sites from each country are plotted with the corresponding national emission densities derived for NH 3 , NO x , and SO 2 . The emissions data in the graphs are the 4-year averages for the period 2007 to 2010, expressed as emissions per unit area of the country per year (t km −2 yr −1 ) (see Sect. 2.6) and ranked in order of increasing emission densities. From the visual comparisons, national mean measured concentrations in each country appear to scale reasonably well with the ranked emission densities. This is supported by further regression analyses which showed significant cor-relation between annual averaged concentrations of NH 3 , NO x , and SO 2 with emission densities of NH 3 (R 2 = 0.49, p< 0.001; Fig. 10a), NO x (R 2 = 0.20, p< 0.05; Fig. 10b), and SO 2 (R 2 = 0.65, p< 0.001; Fig. 10c), respectively (Table 3). The particulate components NH + 4 and NO − 3 were also correlated with emission densities of NH 3 and HNO 3 (Table 3). By contrast, there was no relationship between SO 2− 4 with emission densities of any of the three gases, possibly because of contributions to SO 2− 4 from long-range transport. All regression plots of concentrations against emission densities, including summary statistics, are provided in Fig. S7.  (Fig. S10). Since results from this analysis were similar to the comparisons with national emission densities, they are not included for further discussions in this paper. The purpose of the ranked emission densities is to compare the pollution climate in terms of primary gas emissions (SO 2 , NO 2 , NH 3 ) across the 20 European countries and to see if this is matched by the  DELTA ® measurements. Despite the complex relationship between emissions and concentrations, the pollution gradient in Europe is clearly captured by the present data. At the same time, it also demonstrated the potential application of the DELTA ® approach in providing national concentration fields as evidence to compare against spatial and long-term trends in the national emissions data.

Spatial variability across geographical regions
The form and concentrations of the different gas and aerosol components measured also varied according to geographic region across Europe (Fig. 11). The smallest concentrations (with the exception of SO 2− 4 and Na + ) were in northern Europe (Scandinavia), with broad elevations across other regions. Gas-phase NH 3 and particulate NH + 4 were the dominant species in all regions (Fig. 11). NH 3 showed the widest range of concentrations, with the largest concentrations in western Europe (mean = 2.4 NH 3 m −3 , range = 0.2-7.1 µg NH 3 m −3 , n = 26 in four countries). By contrast, HNO 3 and SO 2 concentrations were largest in high-NO xand high-SO 2 -emitting countries in central and eastern Europe (Sect. 3.3.2). Particulate SO 2− 4 concentrations were however more homogeneous between regions, which may be attributed to atmospheric dispersion and long-range transboundary transport of this stable aerosol between countries in Europe (Szigeti et al., 2015;Schwarz et al., 2016). In the aerosol components, the spatial correlations between NO − 3 , NH + 4 , and NH 3 illustrate the potential for NH 3 emissions to drive the formation and thus regional variations in NH + 4 and NO − 3 aerosol. Particulate SO 2− 4 concentrations in northern Europe (Scandinavia) were similar to other countries, despite having the smallest SO 2 and NH 3 emissions and concentrations (Fig. 9). By comparison, the smaller particulate NH + 4 and NO − 3 concentrations in northern Europe are consistent with the smallest emissions (NH 3 and NO x ) and concentrations of NH 3 and HNO 3 (Fig. 9). As discussed later in Sect. 3.4, the larger SO 2− 4 concentrations reported in northern Europe were flagged up as anomalous from ion balance checks (ratio of NH + 4 : sum anions).

Comparisons by grouped components
In the following sections, variations in concentrations of the different gas and aerosol components according to ecosystem type (crops, grassland, forests, and semi-natural) and in relation to emissions (NH 3 , NO x , and SO 2 ) are further discussed. For ease of interpretation, components are grouped as follows: reduced N (NH 3 , NH + 4 ), oxidized N (HNO 3 , NO − 3 ), S (SO 2 , SO 2− 4 ), HCl, Na + , and Cl − .

Reduced N (NH 3 and NH + 4 )
Broad differences in NH 3 concentrations are observed between the grouped sites, with the largest concentrations at cropland sites, as expected, as these are intensively managed agricultural areas dominated by NH 3 emissions (Fig. 7a). Borgo Cioffi (IT-BCi) is an ecosystem station located in a 15 ha field (arable crops) on the Sele Plain, an agricultural area with intensive buffalo farming in southern Italy, and this provided the highest 4-year average of 8.1 µg NH 3 -N m −3 (cf. group mean = 3.8 µg NH 3 -N m −3 , n = 10) ( Table 4, Table S9). Next highest in this group are the German Gebesee (DE-Geb) and the Belgian Lonzée (BE-Lon) sites, with 4year-average concentrations of 4.9 and 4.8 µg NH 3 -N m −3 , respectively (Table S9). At Gebesee, a decrease in NH 3 concentrations was observed over the 4-year period, falling almost 2-fold from an annual mean of 8.8 µg NH 3 -N in 2007 to 4.8 µg NH 3 -N in 2010 (Table S9). Annual mean concentrations in 2008 (2.9 µg NH 3 -N m −3 ) and 2009 (3.2 µg NH 3 -N m −3 ) were similar but smaller than in 2010. This illustrates the large inter-annual variability in concentrations that can occur even over a short time period. Variability between years may reflect the impact of changes in meteorological conditions on emissions from potential sources, with for example warmer, drier years increasing emissions and concentrations, contrasting with lower emissions and concentrations from the same source in a colder and wetter year. Episodic pollution events can also have a large influence on the annual mean concentration rather than the direct effects of changes in anthropogenic emissions over this short timescale. This suggests that for compliance assessment, an average over several years would provide a more robust basis than individual years. The assessment of trends also needs a longer time series of at least 10 years (Tang et al., 2018a, b;Tørseth et al., 2012;van Zanten et al., 2017). Grassland sites, with NH 3 emissions from grazing and fertilizers, provided the next highest concentrations, with annual averaged concentrations of 2.2 µg NH 3 -N m −3 from the 10 sites in this group (Table 4). Cabauw in the Netherlands (NL-Cab) in this group was the second-highest NH 3 concentration site in the DELTA ® network, after Borgo Cioffi (IT-BCi), with a 4-year-averaged annual concentration of 5.9 µg NH 3 -N m −3 (Table S9). Unlike the Gebesee site (DE-Geb), annual NH 3 concentrations were consistent between years at Cabauw, ranging from annual mean of 6.3 µg NH 3 -N m −3 in 2017 to 5.8 µg NH 3 -N m −3 in 2010 (Table S9).
With the exception of the Czech Republic, the annual averaged NH 3 concentrations scaled reasonably well with the 4-year-averaged mean NH 3 emission density in each country (Figs. 9, 10a1, 10b1) (see also Sect. 3.3.2). In the Czech Republic, measurement was made at a single site, BKFores (CZ-BK1), located at a remote forest location. The 4-yearaveraged emissions in the EMEP grid (0.1 • × 0.1 • ) containing the site are very small, at 2 t NH 3 -N yr −1 , compared with an average of 68 t NH 3 -N yr −1 (range = < 0.01 to 567 t NH 3 -N yr −1 ) across the Czech Republic (Fig. S9). The low emissions, combined with the small concentrations measured at BKFores (0.5 µg NH 3 -N m −3 ), suggests it is highly likely to represent concentrations at the low end of the range of NH 3 concentrations that might be expected to be encountered in the Czech Republic. By comparison, Belgium has a similar emission density to the Czech Republic, but the mean concentrations from three sites (2.6 µg NH 3 -N m −3 ) encompassed sites located in cropland areas (Lonzée BE-Lon, 4.7 µg NH 3 -N m −3 ) and forest sites (Braschaat BE-Bra, 2.8 µg NH 3 -N m −3 , and Vielsalm BE-Vie, 0.4 µg NH 3 -N m −3 ) ( Table S9). The markedly high concentrations of NH 3 across the NEU network indicates that contributions by emission and deposition of NH 3 would be a major contributor to the effects of N r on sensitive habitats. In comparing the annual averaged NH 3 concentration with the revised UNECE "critical levels" of NH 3 concentrations , the lower limit of 1 µg NH 3 m −3 annual mean for the protection of lichens and bryophytes was exceeded at 63 % of sites (40 sites in 15 countries) (Table S10). Even the higher 3 µg NH 3 m −3 annual mean for the protection of vegetation was still exceeded at 27 % of sites (17 sites in 10 countries) (Table S10). Most notably, all four sites from the Netherlands were in exceedance of both the 1 and the 3 µg NH 3 m −3 thresholds. The large concentrations in the Netherlands highlight the high levels of NH 3 that semi-natural and forest areas are exposed to within an intensive agricultural landscape, where 117 out of the 166 Natura 2000 areas were reported to be sensitive to nitrogen input (Lolkema et al., 2015). A recent assessment estimated that critical loads for eutrophication were exceeded in virtually all European countries and over about 62 % of the European ecosystem area in 2016 (EMEP, 2018). In particular, the highest exceedances occurred in the Po Valley (Italy), the Dutch-German-Danish border areas, and north-western Spain, where the highest NH 3 concentrations have been measured in this network. Since NH 3 is preferentially deposited to semi-natural and forests (high V d to these ecosystem types; Sutton et al., 1995), then NH 3 will dominate dry N deposition and exert the larger ecological impact. In Flechard et al. (2011), dry NH 3 -N deposition from the first 2 years of NH 3 measurement in the NEU DELTA ® network was esti-mated to contribute between 25 % and 50 % of total dry N deposition in forests, according to models. The fraction is larger in short semi-natural vegetation since V d for NH + 4 and NO − 3 is smaller in short vegetation than in forests (Flechard et al., 2011).

Comparison with NH 3 data from the Dutch LML network
The 4-year-averaged NH 3 concentrations from the Dutch LML air quality network (see Sect. 2.7.1) for the period 2007 to 2010 are plotted alongside the NH 3 measurements made at the four Dutch sites in the DELTA ® network (Fig. 9a). The 4year-averaged concentrations from the eight LML sites were between 1.5 and 15 µg NH 3 -N m −3 , highlighting the high concentrations and spatial variability in concentrations in the Netherlands. The mean NH 3 concentrations measured at the four Dutch sites in the DELTA ® network of 2.9 µg NH 3 -N m −3 (Horstermeer, NL-Hors; semi-natural) to 5.9 µg NH 3 -N m −3 (Cabauw, NL-Cab; grassland) were within the range of concentrations measured in the Dutch LML network.

Comparison with NH 3 data from the UK NAMN network
The 4-year-averaged NH 3 concentrations calculated from the 72 sites in the NAMN (see Sect. 2.7.2) for the period 2007 to 2010 were smaller than the Dutch LML network, ranging from 0.05 to 6.7 µg NH 3 -N m −3 , consistent with smaller NH 3 emissions from the UK (Fig. 9a). In a joint collaboration between the UK and Dutch networks, inter-comparison of NH 3 measurements by the DELTA ® method (monthly) with the AMOR wet chemistry system (hourly; van Zanten et al., 2017) was carried out at the Zegweld site (ID 633) in the Dutch LML network (van Zanten et al., 2017) between 2003 and 2015. Good agreement was provided, lending support for comparability between the independent measurements reported in Tang et al. (2018a).
Particulate NH + 4 concentrations across the 64 sites were more homogeneous than NH 3 , varying over a narrower range between 0.13 µg NH + 4 -N m −3 at Sodankylä (Finland, FI-Sod) and 2.1 µg NH + 4 -N m −3 at Borgo Cioffi (Italy, IT-BCi) (Fig. 7, Table S11). By comparison, the difference in NH 3 between the smallest (0.07 µg NH 3 -N m −3 at Lompolojänkkä, Finland, FI-Lom) and largest (8.1 µg NH 3 -N m −3 at Borgo Cioffi, Italy, IT-BCi) concentrations varied by a factor of 110 (Fig. 7, Table S10). Secondary aerosols have longer atmospheric lifetimes and will therefore vary spatially much less than their precursor gas concentrations. While the concentrations of NH 3 vary at a local to regional level owing to large numbers of sources at ground level and high deposition in the landscape, NH + 4 is less influenced by proximity to NH 3 emission sources and varies in concentration at regional scales Tang et al., 2018a).
The UK and Irish sites have the next smallest NH + 4 concentrations of 0.4 and 0.5 µg NH + 4 -N m −3 (cf. mean of all countries = 0.74 µg NH + 4 -N m −3 ). Particulate NH + 4 data from the UK NAMN (Tang et al., 2018a) are also included for comparison. The 4-year-average concentrations from the 30 sites (0.5 µg NH + 4 -N m −3 , range = 0.14 to 1.0 µg NH + 4 -N m −3 ) are comparable with the mean of 0.40 µg NH + 4 -N m −3 (range = 0.2 to 0.9 µg NH + 4 -N m −3 ) from just 4 sites in the NEU network. A combination of lower emissions of precursor gases (Fig. 9) and being farther away from the influence of long-range transport of NH + 4 aerosols from the higher-emission countries on mainland Europe may be contributing factors to the small NH + 4 concentrations measured in the UK and Ireland.

Oxidized N (HNO 3 and NO − 3 )
The percentage mass contribution of oxidized N (sum of HNO 3 and NO − 3 ; µg N m −3 ) to the total gas and aerosol species measured was on average 13 % (range = 2 %-24 %) (Fig. 8). This compares with 41 % (range = 17 %-70 %) from reduced N (sum NH 3 and NH + 4 ; µg N m −3 ) and 23 % (range = 7 %-53 %) from sulfur (sum of SO 2 and SO 2− 4 ; µg S m −3 ) (Fig. 8). DELTA ® measurements of HNO 3 also include contributions from co-collected oxidized N species such as HONO (see Sect. 2.2.3) and are therefore an upper estimate that may in some cases be twice as large as the actual HNO 3 concentration, based on observations in the UK (Tang et al., 2018b; correction factor of 0.45) and from the parallel DELTA ® measurements made at Fougéres (FR-FgsP) (Fig. S5). At this site, HNO 3 measurement with NaClcoated denuders provided an annual mean concentration of 0.08 µg HNO 3 -N m −3 , compared with 0.19 µg HNO 3 -N m −3 measured on carbonate-coated denuders from the main site (FR-Fgs) (Table S7). With this caveat in mind, uncorrected annual mean HNO 3 concentrations were in the range of 0.03 µg HNO 3 -N at Kaamenan (Finland; FI-Kaa) to 0.47 µg HNO 3 -N at Braschaat (Belgium; BE-Bra) (Table S7). In Fig. 9b, HNO 3 concentrations are compared with NO x emissions, the precursor gas for secondary formation of HNO 3 . Overall, a weak but significant correlation was observed between concentrations of HNO 3 and NO x emission densities across the 20 countries (R 2 = 0.2, p< 0.05) (Fig. 10b, Table 3, Fig. S2). Russia has the lowest NO x emission densities (0.04 t NO x -N yr −1 ), but HNO 3 from the single site (0.15 µg HNO 3 -N m −3 ) is larger than the smallest concentrations measured in Finland, Norway, and Sweden (annual average < 0.1 µg HNO 3 -N m −3 ). HNO 3 formation by photochemical processes may be enhanced in hotter, sunnier summer weather in Russia. Since SO 2 concentrations (mean = 0.49 µg SO 2 -S) at the Russian site (RU-Fyo) are in molar excess over the low levels of NH 3 (mean = 0.32 µg NH 3 -N m −3 ), removal of HNO 3 by reaction with NH 3 will also be limited. HNO 3 concentrations in the UK and Ireland are marginally higher than the Scandinavian countries. Here, the annual averaged concentrations of HNO 3 are similar (0.10 vs. 0.09 µg m −3 ) (Table S12), despite NO x emissions density (t km −2 yr −1 ) in the UK being 3 times larger than in Ireland (Fig. 9b). HNO 3 concentrations on the European continent were generally higher (0.2-0.4 µg HNO 3 -N m −3 ).
In the UK, HNO 3 data are also available on a wider spatial scale from the AGANet (Sect. 2.7.2; Tang et al., 2018b). The 4-year-average concentrations of HNO 3 from 30 sites in the AGANet are plotted alongside the NEU HNO 3 data from the four UK sites in its network in Fig. 9b. The UK HNO 3 data on the UK-AIR database (https://uk-air.defra.gov.uk/, last ac-cess: 25 November 2019) have been corrected for HONO interference with a 0.45 correction factor (see Tang et al., 2018b). For consistency in Fig. 9b, the UK raw uncorrected HNO 3 data are used for the present comparison. The 30-site mean (0.17 µg HNO 3 -N m −3 ) was higher than from just 4 UK sites in the NEU network (0.10 µg HNO 3 -N m −3 ). The range of concentrations was also wider, from 0.03 µg HNO 3 -N m −3 at a remote background site in Northern Ireland to 0.77 µg HNO 3 -N m −3 at a central London urban site, where interference from HONO and NO x in HNO 3 determination is likely to be larger 2018b).
Like particulate NH + 4 , NO − 3 concentrations are also correlated with emission densities of NH 3 (R 2 = 0.57, p< 0.001, n = 20) and NO x (slope = 0.15, R 2 = 0.44, p< 0.01, n = 20) but not with SO 2 (Table 3, Fig. S7). The smallest NO − 3 concentrations were again in Sweden, Norway, and Finland with low NH 3 and NO x emissions and also smallest concentrations of HNO 3 , SO 2 , and NH + 4 in the network (Fig. 9). The largest NO − 3 concentrations were measured in the Netherlands, with a mean of 0.92 µg NO − 3 -N m −3 , compared with a network average of 0.39 µg NO − 3 -N m −3 (Fig. 9e, Table S13). The higher NO − 3 concentrations correlated well with the high NH 3 , HNO 3 , and NH + 4 concentrations in the Netherlands (Fig. 9). This suggests that concentrations of NO − 3 are linked to local formation of NH 4 NO 3 , which is dependent on concentrations of NH 3 and HNO 3 , and also to the influence of meteorology on transport of NH 4 NO 3 between countries on mainland Europe and export out of Europe. Countries in Scandinavia such as Sweden, Norway, and Finland and in the British Isles are farthest from the influence of long-range transboundary transport from Europe, with concentrations of NH 4 NO 3 that are smaller than on the continent.

Sulfur (SO 2 and SO 2− 4 )
Annual averaged SO 2 concentrations measured across the network were between 0.9 and 2.3 µg SO 2 -S m −3 (Fig. 9c, Table S14). By comparison, the EMEP network of 58 urbanbackground sites reported annual mean concentrations of 0.03 and 5.5 µg SO 2 -S m −3 over the same period, with the largest SO 2 concentrations from North Macedonia and Serbia. Since these high-emitting countries were not included in the DELTA ® network, the range of SO 2 concentrations is smaller. Together, the small SO 2 concentrations reflect the substantial reductions in SO 2 emissions across Europe (−74 % between 1990 and 2010) and large reductions in ambient concentrations and deposition of sulfur species across Europe during the last decades (EMEP, 2016). SO 2 concentrations were also correlated with SO 2 emission density (R 2 = 0.65, p< 0.001, n = 20) in each country (Fig. 10c, Table 3). The smallest and largest SO 2 annual average concentrations corresponded with the lowest emissions in Norway and highest in the Czech Republic (Fig. 9c). By contrast, SO 2 concentrations from the only measurement site, Bugac, in Hungary (HU-Bug) are much higher than expected on the basis of SO 2 emission density estimated for the country. This suggests that Bugac is likely to be affected by proximity to sources. This contrasts with the BKFores site in the Czech Republic (CZ-BK1), which had smaller NH 3 concentrations due to its location away from sources.
Following emission, SO 2 disperses and undergoes chemical oxidation in the atmosphere to form SO 2− 4 both in the gas phase and in cloud and rain droplets (Baek et al., 2004;Jones and Harrison, 2011). Particulate SO 2− 4 produced is generally associated with NH + 4 and NO − 3 (see Sect. 3.4). The regional pattern of SO 2− 4 was similar to, and correlated well with, particulate NH + 4 and NO − 3 (Fig. 9g), suggesting well-mixed air on the continent since (NH 4 ) 2 SO 4 is stable and long-lived. Countries in the British Isles (UK and Ireland) and in Scandinavia (Sweden, Norway, Finland) have smaller concentrations of SO 2− 4 (Table S15). They are located far enough away from sources and activities on continental Europe such that they are less influenced by the emissions from central Europe.
As discussed earlier, particulate NH + 4 and NO − 3 concentrations were smallest in the Scandinavian countries, which corresponded with low emission densities of the precursor gases NH 3 and NO x . By analogy, since these countries also have the lowest emission densities of SO 2 (Fig. 9c), particulate SO 2− 4 concentrations would be expected to be similarly low. Particulate SO 2− 4 in Finland and Norway (mean = 0.34 µg SO 2− 4 -S m −3 ) and Sweden (mean = 0.37 µg SO 2− 4 -S m −3 ) were however comparable with concentrations on mainland Europe (range = 0.33 to 1.0 µg SO 2− 4 -S m −3 ) and larger than the UK (0.18 µg SO 2− 4 -S m −3 ) and Ireland (0.24 µg SO 2− 4 -S m −3 ) (Fig. 9f). An ion balance check on the ratio of equivalent concentrations of NH + 4 to the sum of NO − 3 and SO 2− 4 (see Sect. 3.4) was less than 0.5. Since NH + 4 is a counter-ion to NO − 3 and SO 2− 4 formation, the imbalance suggests that SO 2− 4 concentrations may be overestimated at the sites in Sweden, Norway, and Finland.
HCl, Cl − , and Na + The average concentrations of HCl across the network were of low magnitude, with limited variability, ranging from 0.07 in Russia to 0.36 µg HCl-Cl − m −3 in Portugal (Fig. 9d). At a site level, HCl concentrations varied between 0.06 at Renon (Italy; IT-Ren -inland location) to 0.48 µg HCl-Cl − m −3 at Espirra (Portugal, PT-Esp -coastal location) (Table S16). In the UK AGANet network, the highest concentrations of HCl were found in the source areas in the south-east and south-west of England and also in central England, north of a large coal-fired power station (Tang et al., 2018b). HCl emissions and concentrations in the atmosphere are mostly derived from combustion of fossil fuels (coal and oil), biomass burning, and the burning of municipal and domestic waste in municipal incinerators (Roth and Okada 1998;McCulloch et al., 1999;Ianniello et al., 2011). Several manufacturing processes, including cement production, also emit HCl (Mc-Culloch et al., 1999). At coastal sites, HCl released from the reaction of sea salt with HNO 3 and H 2 SO 4 can be a significant source (Roth and Okada 1998;Keene et al., 1999;McCulloch et al., 1999;Ianniello et al., 2011). UK is the only country with available HCl emission estimates (https: //naei.beis.gov.uk/data/, last access: 3 January 2020). Emissions of HCl in the UK (mainly from coal burning in power stations) have declined to very low levels, from 74 kt in 1999 to 5.7 kt in 2015. The 4-year-averaged emission density for HCl for the period 2007 to 2010 was just 0.05 tonnes HCl-Cl − km −2 yr −1 , although HCl emissions could still pose a threat to sensitive habitats close to sources (Evans et al., 2011). The low HCl concentrations measured in the network would suggest that the shift in Europe's energy system from coal to other sources has contributed to low HCl emissions (UK) and concentrations (observed across the network).
Particulate Cl − on the other hand is predominantly marine in origin, with sea salt (NaCl) as the most significant source (Keene et al., 1999). Molar concentrations of Cl − and Na + are seen to be similar in most countries, demonstrating close coupling between the two components (Fig. 9h). The largest concentrations of Na + and Cl − occurred in coastal countries such as the UK, Ireland, Netherlands, and Portugal, with the highest of country-averaged annual concentrations of 1.6 µg Cl − m −3 and 0.9 µg Na + m −3 from Ireland (Tables S16 and S17). Data from the 30 sites in the UK AGANet network showed a wider range of Cl − and Na + concentrations (Fig. 9h). The AGANet site with the largest 4-year-averaged annual concentrations of 3.8 µg Cl m −3 and 2.0 µg Na + m −3 is Lerwick from the east coast of the Shetland Islands, exposed to the North Atlantic.
Farther away from the coastal influence of marine aerosol, the smallest concentrations of Cl − and Na + were measured in landlocked countries such as Germany (mean of all sites = 0.27 µg Cl − m −3 and 0.15 µg Na + m −3 ). Concentrations in Hungary, Poland, the Czech Republic, and Russia were also low, but inferences about these countries are necessarily limited by measurements at a single site in each of these countries. At coastal sites in Norway (NO-Bir) and Sweden (SE-Nor and SE-Sk2), the very low particulate Cl − concentrations (< 0.1-0.3 µg m −3 ) and high Na : Cl molar ratios (3-5) are anomalous. It is possible for sea salt to be depleted in Cl − (through the loss of HCl gas) by the reaction of NaCl particles with atmospheric acids (Finalyson-Pitts and Pitts, 1999;Keene et al., 1999), leading to high Na : Cl ratios for sea salts transported over long distances. The coastal locations of these sites (Fig. 2) suggest that they are more likely to be influenced by freshly generated marine aerosols (cf. coastal sites in UK and Ireland), and larger concentrations of sea salt (Na + and Cl − ) and a 1 : 1 relationship between Na + and Cl − are expected. The Cl − concentrations are likely to be underestimated at these sites (see Sect. 3.2.3) and are further discussed in the next section (Sect. 3.4).

Correlations between gas and aerosol components
Regression analyses were carried out between the mean molar-equivalent concentrations of all inorganic gas and aerosol components measured at each site (n = 66; Fr-FgsP and UK-AMoP excluded) in the NEU network, with summary statistics provided in Table 5. With the exception of SO 2 vs. HCl (R 2 = 0.05, p> 0.05), the gases were positively correlated with each other, possibly due to similarities in the regional distribution of their emissions and concentrations. Comparing the mean molar concentrations of NH 3 with SO 2 and HNO 3 showed that NH 3 was on average 6-fold and 7fold higher, respectively, whereas molar concentrations of SO 2 and HNO 3 were similar (Table 6, Fig. 12). The molar ratio of NH 3 to the sum of all acid gases (SO 2 , HNO 3 , and HCl) was on average 3 (Table 6, Fig. 12), confirming that there is a surplus of the alkaline NH 3 gas to neutralize the atmospheric acids in the atmosphere, similar to that observed in the UK (Tang et al., 2018b). With the more substantial decline in emissions of SO 2 , compared with a more modest reduction in NO x , the concentrations of SO 2 are at a level where it is no longer the dominant acid gas, such that HNO 3 and HCl are together contributing a larger fraction of the total acidity in the atmosphere in the present assessment.
Non-sea-salt SO 2− 4 (nss-SO 2− 4 ) was also estimated from the SO 2− 4 and Na + data (see Sect. 2.2.1). The nss-SO 2− 4 is estimated to comprise on average 25 % (range = 3 %-83 %, n = 187) of the measured total SO 2− 4 aerosol (Table 8). This demonstrates that sea salt SO 2− 4 (ss-SO 2− 4 ) aerosol makes up a large and variable fraction of the total SO 2− 4 measured, consistent with observations of the contribution by ss-SO 2− 4 to the total SO 2− 4 in precipitation observed in the wet-deposition measurements in this study (Fig. 11) and across Europe (ROTAP, 2012). Regression of nss-SO 2− 4 vs.
NH + 4 (slope = 0.27, R 2 = 0.30) was not significantly different from the regression of SO 2− 4 vs. NH + 4 (slope = 0.27, R 2 = 0.28) ( Table 5). This suggests that NH + 4 is mainly associated with the nss-SO 2− 4 . Correlation between NH + 4 and the sum of anions (NO − 3 + SO 2− 4 ) is an important point of discussion (Table 7) as the ion balance serves as a quality check for the aerosol measurement. Due to some outliers in the comparison, the correlation between NH + 4 and SO 2− 4 (R 2 = 0.28; Fig. 13b) is weaker than between NH + 4 and NO − 3 (R 2 = 0.75; Fig. 13c, Table 7). The outliers were measurements made by NILU and CEAM, although these vary according to monitor- Significance level: * p < 0.05, * * p < 0.01, * * * p < 0.001; ns: non-significant (p > 0.05). 3 n cording to measurements made by each of the seven laboratories (Fig. S11). Regressions for CEAM and NILU show the vast majority of data points below the 1 : 1 line, indicating a systematic underestimation of particulate Cl − concentrations. The other five laboratories (INRAE, MHSC, SHMU, UKCEH, and vTI) all have data points close to the 1 : 1 line, with larger scatter both above and below the 1 : 1 line at lower concentrations. In Fig. 13f, a new regression line has therefore also been fitted, where outlier data with Na : Cl ratios > 2 from NILU (13 out of 16 sites) and CEAM (all 3 sites) have been removed. Exclusion of the outlier data points provided a regression line that is not significantly different from unity (slope = 1.02), with an R 2 value of 0.95 (p< 0.001). The near-1 : 1 relationship between particulate Na + and Cl − is consistent with their origin from sea salt (NaCl).
The ion balance checks, together with the regular PT exercises and field inter-comparisons, therefore provided the platform against which to assess data quality and comparability of measurements between laboratories. This shows that overall, with the exception of a few identified outlier measurements, the laboratories are performing well and providing good agreement.

Seasonal variability in gases and aerosol
The time series of monthly averaged concentrations for the period 2006 to 2010 have been plotted to examine seasonality in the different gas and aerosol components according to ecosystem types (crops, grassland, semi-natural, and forests) (Fig. 14) and geographical regions (Fig. 15). Distinct seasonality was observed in the data, influenced by seasonal changes in emissions, chemical interactions, and the influence of meteorology on partitioning between the main inorganic gases and aerosol species.

NH 3
Distinctive and contrasting features in the seasonal cycle are observed, with the largest concentrations at cropland sites and smallest at semi-natural and forest sites (Fig. 14a). Similar to those observed in the annual mean concentrations (Fig. 9, 11), the monthly concentrations are also smallest in northern Europe and largest in western Europe (Fig. 15a).

Semi-natural sites
There are two distinct peaks in the seasonal cycle of grouped semi-natural sites: in April (mean = 2.2 µg NH 3 m −3 , n = 12) and in July (mean = 1.9 µg NH 3 m −3 , n = 12) (Fig. 14a). Since these sites are located away from agricultural sources, the seasonality in NH 3 concentrations is mostly governed by changes in environmental conditions and regional changes in NH 3 emissions. The differences in concentrations between the summer and winter at these sites were by a factor of 3, with the smallest concentrations in wintertime (December and January), when low temperatures and wetter conditions decrease NH 3 emissions from regional agricultural sources while favouring a thermodynamic shift from gaseous NH 3 to the aerosol NH + 4 phase. Conversely, warm, dry conditions in summer increases surface volatilization of NH 3 from lowdensity grazing livestock and wild animals and favour a thermodynamic shift to the gaseous (NH 3 ) phase, producing the summer peak. Vegetation is another potential source at these background sites under the right conditions Massad et al., 2010). A complex interaction between atmospheric NH 3 concentrations and vegetation can lead to both emission and deposition fluxes, known as "bidirectional exchange", dependent on relative differences in concentrations. This process is controlled by the so-called "compensation point", defined as the concentration below which growing plants start to emit NH 3 into the atmosphere (Flechard et al., 1999;Massad et al., 2010;Sutton et al., 1995). At sites distant from intensive farming and emissions, the bidirectional exchange with vegetation will partly control NH 3 concentrations. Inclusion of bidirectional exchange in dispersion modelling of NH 3 by incorporating a "canopy compensation point" is shown to improve model results for NH 3 concentrations in remote areas (e.g. Smith et al., 2000;Flechard et al., 1999Flechard et al., , 2011Massad et al., 2010). The larger peak in April at these sites on the other hand suggests the influence of emissions from agricultural sources, e.g. from land spreading of manures.

Forest sites
The average seasonal cycle from the forest sites is similar to that of the semi-natural sites but diverged over the summer months (Fig. 14a). Here, the seasonal profile is characterized by the absence of any peaks in summer, with concentrations plateauing between May and August. Studies have shown that atmospherically deposited N is taken up by forest canopies since growth in forest ecosystems is commonly limited by the availability of N (Sievering et al., 2007), and tree canopies are a potential sink for atmospheric NH 3 (Fowler et al., 1989;Theobald et al., 2001). The capture and uptake of NH 3 during the growing seasons over the summer period could therefore account for the absence of a summer peak in NH 3 concentrations at forest monitoring sites, although a similar effect would also be expected for semi-natural sites.

Cropland sites
Fertilizers and arable crops are significant sources of NH 3 emissions and concentrations in an intensive agricultural landscape. Sites in this group showed considerably higher monthly mean monitored NH 3 concentrations than the other groups (Fig. 14a). A more complex seasonal pattern can be seen, with three peaks in NH 3 concentrations. Concentrations here are also lowest in the winter, although the wintertime concentrations are 3 times larger than semi-natural and forest sites, reflecting the elevated regional background in NH 3 concentrations located within agricultural landscapes. This rises rapidly with improving weather conditions and peaks in the spring to coincide with the main period for manure spreading and fertilizer application before the sowing of arable crops (Hellsten et al., 2007). The distinct spring-time maxima in NH 3 also reflects implementation of the Nitrates Directive (EC, 1991), which prohibits manure spreading in winter. In summer, the second peak in NH 3 concentrations may be associated with increased land surface emissions promoted by warm, dry conditions and possibly with the application of fertilizers. The smaller autumn peak is also expected to be related to seasonal farming activities (e.g. manure spreading). The key drivers for seasonal variability in NH 3 concentrations at crops sites are therefore a combination of seasonal changes in agricultural practices (e.g. timing of fertilizer/manure applications) and climate that will affect emissions, concentrations, transport, and deposition of NH 3 .

Grassland sites
An additional major source of NH 3 in this group of sites is expected to come from grazing emissions and housed livestock (e.g. cattle). Concentrations in this group of sites were generally 2-3 times larger than semi-natural sites (Fig. 14a), attributed to the increased emissions and concentrations from livestock (Hellsten et al., 2007). The spring peak is related to the practice of fertilizer and manure being spread on grazing fields to aid spring grass growth, which will be cut for hay and silage later in the year. NH 3 concentrations in June and July are smaller than in spring or late summer, possibly because grass will be actively growing with possible uptake and removal of NH 3 from the atmosphere . The concentrations are also larger in summer than winter, with warmer conditions promoting NH 3 volatilization and thermodynamic shift in NH 4 NO 3 to the gas phase.

European regions
The seasonal profiles of NH 3 for central and western European regions were similar, characterized by a large peak in spring that is likely to be agriculture-related (Fig. 15a), as observed at cropland sites (Fig. 14a). While the peak concentrations in both regions are of comparable magnitude (central = 2.6 µg NH 3 m −3 , western = 2.8 µg NH 3 m −3 ), winter concentrations in central Europe (0.6 µg NH 3 m −3 ) were 3 times smaller than the west (1.5 µg NH 3 m −3 ). This may be related to lower regional background in NH 3 concentrations and/or suppressed emissions in colder temperatures of central Europe in winter. By contrast, eastern and southern European regions have a broad peak in summer, although the eastern region also has a second peak in October (likely agriculture-related). The smallest concentrations were found in northern Europe, with the lowest NH 3 emissions (Fig. 9). The three peaks in the profile show elevated concentrations in summer driven by warming temperatures, with the spring and autumn peaks attributed to influence from NH 3 emissions from agricultural sources.

HNO 3
The seasonal distribution in HNO 3 is similar between the different ecosystem groups, varying only in magnitude of concentrations (Fig. 14c), and reflects the secondary nature of this component that is formed from oxidation of NO x (Fahey et al., 1986;ROTAP, 2012). Since the HNO 3 data are actually the sum of HNO 3 and HONO, with a small contribution from NO 2 (see Sect. 2.2.3), the temporal patterns seen are likely to be the superimposed profiles of both HNO 3 and HONO. In the studied region, NO 2 is predominantly from vehicular emissions, which are not expected to show large seasonal variations and should therefore exert a negligible effect on the temporal patterns in HNO 3 . With this caveat in mind, HNO 3 concentrations in the crops group are up to 2 times larger than the grassland group, while the smallest concentrations are in the semi-natural group. This is likely related to the proximity of sites in the different groups to combustion sources. A weak seasonal cycle is seen in the secondary HNO 3 air pollutant in all cases, with slightly higher concentrations in late winter, spring, and summer and lowest in March and November. The reaction of NO 2 with the OH radical is an important source of HNO 3 during daytime, whereas N 2 O 5 hydrolysis is considered an important source of HNO 3 at nighttime (Chang et al., 2011). Larger HNO 3 concentrations in summer are therefore from increased OH radicals for reaction with NO 2 to form HNO 3 . Similarly, higher concentrations of ozone in spring in Europe (EMEP, 2016) can potentially increase HNO 3 concentrations in springtime. Conversely, HNO 3 concentrations are lower in winter, when oxidative capacity is less. Seasonal variability in HNO 3 will also be influenced by gas-aerosol-phase equilibrium. In the atmosphere, HNO 3 reacts reversibly with NH 3 , forming the semi-volatile NH 4 NO 3 aerosol if the necessary concentration product [HNO 3 ] · [NH 3 ] is exceeded (Baek et al., 2004;Jones and Harrison et al., 2011). Because of this process, the prime influences upon HNO 3 concentrations at sites where NH 4 NO 3 is formed are expected to be ambient temperature, relative humidity, and NH 3 concentrations that affect the partitioning between the gas and aerosol phase (Allen et al., 1989;Stelson and Seinfeld, 1982). The availability of surplus NH 3 in spring (Sect. 3.5.1) would tend to reduce HNO 3 and increase NH 4 NO 3 formation, which is reflected in the reduced HNO 3 concentrations observed in March, when NH 3 is at a maximum. In summer, warmer, drier conditions promote volatilization of the NH 4 NO 3 aerosol, increasing the gas-phase concentrations of HNO 3 and NH 3 relative to the aerosol phase. Seasonality in HNO 3 is therefore complex, related to traffic and industrial emissions, photochemistry, and HNO 3 : NH 4 NO 3 partitioning.
An analysis of the same data grouped according to geographical region revealed distinctive cycles in HNO 3 in eastern and southern Europe (Fig. 15c). These two regions showed the highest concentrations in summer and lowest in winter, consistent with enhanced photochemistry in warmer, sunnier climates, and thermodynamic equilibrium favouring gas-phase HNO 3 (Fig. 15c). Summertime peak concentrations in NH 3 were also observed in these two regions (Fig. 15a). In comparison, the seasonal profiles of HNO 3 in other regions were similar to those described for different ecosystem types (Fig. 14c). and Cl − in the measurements are mainly sea salt (see also spatial correlation in Sect. 3.4). The highest concentrations of Na + and Cl − during winter months would be consistent with increased generation and transport of sea salt generated by more stormy weather from marine sources during those periods (O'Dowd and de Leeuw, 2007).

Bulk wet-deposition measurements
Annual mean wet deposition of chemical species measured at the NEU bulk sampling sites was estimated by combining measured concentrations with annual precipitation. Site changes also occurred during the operation of the bulk wetdeposition network, with some sites closed and new sites added. At Mitra (PT-Mi3), contamination of the rain samples from bird strikes resulted in the rejection of a large proportion of the monthly data, and this site was excluded from the data analysis. In total, 12 sites provided 2 years of monthly data, with a further 5 sites providing 1 year of monthly data over the period 2008 to 2010. Due to differences in start and end dates for bulk measurements between the sites, the annual mean data derived are for 12-month periods or 2 × 12month periods and not from calendar years.
Annual mean wet-deposition data for the 17 sites from seven countries (Belgium, France, Germany, Italy, Poland, Spain, and Switzerland) are summarized in Fig. 16. Using Na + as a tracer for sea salt (Keene et al., 1986), nss-SO 2− 4 concentrations were also estimated from the total SO 2− 4 (see Sect. 2.2.2) and are included for comparison. Since the measurements were made at a limited number of sites across Europe, there is insufficient information to make inferences about spatial differences in concentrations. Detailed assess-ments of extensive precipitation chemistry across Europe are made elsewhere, for example from the EMEP wet-deposition networks (EMEP, 2016;Trseth et al., 2012). What the NEU bulk network data clearly shows is that N r components in rain also exceed those of S (Fig. 16), as was observed in the atmospheric data. The mean proportional contribution of total N (NH + 4 and NO − 3 ) to the sum total of all wet-deposited species measured (by mass) was 19 % (range = 3 %-39 %), compared with a smaller 9 % (range = 1 %-19 %) contribution from nss-SO 2− 4 (Table S14). Wet-deposited N (NH + 4 and NO − 3 ) was on average 2 times higher than nss-SO 2− 4 , similar to that seen in the relative proportion of total N r (sum of NH 3 , NH + 4 , HNO 3 , NO − 3 ) to total S (sum of SO 2 , SO 2− 4 ) in the atmospheric data (Sect. 3.3.4). Similar to the atmospheric data (Sect. 3.3.4), a considerable fraction of the wet-deposited components was made up of sea salt (Na + and Cl − ), with the sum of Na + and Cl − contributing on average 50 % of the total wet-deposited components (range = 20 %-84 %, n = 17). Contributions by the other base cations Ca 2+ and Mg 2+ gave a further 20 % (range = 8 %-41 %, n = 17) ( Table S14).
The wet-deposition data on NH + 4 and NO − 3 , combined with a wider precipitation chemistry dataset (e.g. from EMEP and other national precipitation networks), were used to estimate total N r deposition to a site (Flechard et al., 2011. Together, the dry (DELTA ® network) and wet N r estimates (NEU bulk network, combined with data from other national precipitation chemistry networks) are used to compare with EMEP models and to examine the interactions between N r supply and greenhouse gas exchange at the NEU DELTA ® sites, presented in a separate paper by Flechard et al. (2020). The wet-deposition measurements in this paper highlight where DELTA ® and bulk wet-deposition data are co-located and provide parallel information on gas and aerosol concentrations (for dry-deposition modelling) and wet deposition at those sites. The co-located data are important for deriving N budgets and linking to ecosystem response (e.g. Flechard et al., 2020) and invaluable for modellers.

Conclusions
The NitroEurope DELTA ® network has provided for the first time a comprehensive quality-assured multi-annual dataset on reactive gases (NH 3 , HNO 3 , SO 2 , HCl) and aerosols (NH + 4 , NO − 3 , SO 2− 4 , Cl − ) across the major gradients of emission densities, ecosystem type, and climatic zones of Europe. By sharing the method and protocol with several European laboratories and developing synergies with established infrastructure (e.g. CarboEurope network and EMEP field sites), it has proven possible to establish a large-scale network within a relatively short timescale and with low costs. Key elements were a harmonized methodology and the implementation of quality protocols that included regular lab-oratory and field inter-comparisons to monitor and improve performance.
At the same time, the concurrent measurement of the gas and aerosol components permitted an assessment of the atmospheric composition and the spatial and seasonal characteristics in the gas and aerosol phase of these components. The dataset has also been used to develop estimates of sitebased N r dry-deposition fluxes across Europe, including supporting the development and validation of long-range transport models. Combined with estimates of wet deposition (from NEU bulk wet-deposition network and other networks such as EMEP), an assessment of the interactions between N supply and greenhouse gas exchange was addressed in a separate paper by Flechard et al. (2020) using N r and CO 2 flux data from the co-location of the NEU DELTA ® with Car-boEurope Integrated Project sites.
Two key features have emerged in the data. The first is the dominance of NH 3 as the largest single component at the majority of sites, with molar concentrations exceeding those of HNO 3 and SO 2 combined. As expected, the largest NH 3 concentrations were measured at cropland sites, in intensively managed agricultural areas dominated by NH 3 emissions. The smallest concentrations were at remote semi-natural and forest sites, although concentrations in the Netherlands, Italy, and Germany were up to 45 times larger than similarly classed sites in Finland, Norway, and Sweden (< 0.6 µg NH 3 -N m −3 ), illustrating the high NH 3 concentrations that sensitive habitats are exposed to in intensive agricultural landscapes in Europe. The second key feature is the dominance of NH 4 NO 3 over (NH 4 ) 2 SO 4 , with on average twice as much NO − 3 as SO 2− 4 (on a molar basis). A change to an atmosphere that is more abundant in NH 4 NO 3 will likely increase the atmospheric lifetimes and extend the footprint of the NH 3 and HNO 3 gases by the re-volatilization of NH 4 NO 3 in warm weather.
Temporally, peak concentrations in NH 3 for crops and grassland sites occurred in spring, reflecting the implementation of the EU Nitrates Directive that prohibits winter manure spreading. The spring agriculture-related peak was seen even at semi-natural and forest sites, highlighting the influence of NH 3 emissions at sites that are more distant from sources. Summer peaks, promoted by increased volatilization of NH 3 but also by gas-aerosol-phase thermodynamics under warmer, drier conditions, were seen in all ecosystem groups, except at forest sites. The seasonality in the NH 3 concentrations thus provided important insights into both the relationship to occurrence of emissions and possible abatement measures to target peak emission periods.
Seasonality in the other gas and aerosol components is also driven by changes in emission sources, chemical interactions, and changes in environmental conditions influencing partitioning between the precursor gases (SO 2 , HNO 3 , NH 3 ) and secondary aerosols (SO 2− 4 , NO − 3 , NH + 4 ). Seasonal cycles in SO 2 were mainly driven by emissions (combustion), with concentrations peaking in winter, except in southern Europe, where the peak occurred in summer. HNO 3 concentrations were more complex, as seasonal variations were affected by photochemistry, meteorology, and gas-aerosol-phase equilibrium. Southern and eastern European regions provided the clearest seasonal cycle for HNO 3 , with the highest concentrations in summer and the lowest in winter, attributed to increased photochemistry in the summer months in hotter climates. In comparison, a weaker seasonal cycle is seen in other regions, with marginally elevated concentrations in late winter, spring, and summer and lowest in March and November. Increased ozone in spring is likely to enhance oxidation of NO x to HNO 3 for forming the semi-volatile NH 4 NO 3 by reaction with a surplus of NH 3 . Cooler, wetter conditions in spring also favour the formation of NH 4 NO 3 , and more of the NH 4 NO 3 remains in the aerosol or condensed phase. This accounts for the higher concentrations of NH + 4 and NO − 3 in spring and the absence of a HNO 3 peak at this time of year. Conversely, increased partitioning to the gas phase in summer decreases NH 4 NO 3 concentrations relative to gas-phase NH 3 and HNO 3 . Particulate SO 2− 4 showed large peaks in concentrations in summer in southern and also eastern Europe, contrasting with much smaller peaks occurring in early spring in other regions. The peaks in particulate SO 2− 4 coincided with peaks in NH 3 concentrations, illustrating the importance of NH 3 in driving the formation of (NH 4 ) 2 SO 4 . Since NH 4 NO 3 is more abundant than (NH 4 ) 2 SO 4 , the seasonality of NH + 4 is likely to be influenced more by the temperature and humidity dependence of the semi-volatile NH 4 NO 3 than by the stable (NH 4 ) 2 SO 4 . This is supported by similarity in the seasonal profiles of NH + 4 and NO − 3 at all sites, demonstrating temporal as well as regional correlation between these two components.
Data from the network showed that critical levels of 1 and 3 µg NH 3 m −3 for the protection of lichens and bryophytes and vegetation were exceeded at 62 % and 27 % of the sites, respectively. At the same time, NH 3 dry deposition will also contribute to a significant fraction of deposited acidity and total N deposition to sensitive habitats, along with NH + 4 and HNO 3 dry deposition and wet-deposited NH + 4 and NO − 3 . Although the concentrations of SO 2 have fallen to very low levels at all sites (< 1 µg SO 2 -S m −3 ), SO 2 will continue to be important in contributing to the exceedance of acidification in European ecosystems (EEA, 2019) since SO 2 has a higher acidification potential than NO x (0.70 kg SO 2 = 1 kg equivalent NO 2 in acidity) (see Hauschild and Wenzel, 1998). Changes in the relative concentrations of the pollutant gases captured in the data suggest that the deposition rates of SO 2 and NH 3 will increasingly be controlled by the molar ratio of NH 3 to combined acidity (sum of SO 2 , HNO 3 , and HCl), and deposition models should take these changes into account. Indications from the current and projected trends in emissions of SO 2 , NO x , and NH 3 are that NH 3 and NH 4 NO 3 will continue to dominate the inorganic pollution load over the next decades, contributing to ecosystem effects through acid and N deposition. The growing relative importance of NH 3 and NH + 4 to total acidic and total N deposition indicates that strategies to tackle acidification and eutrophication need to include measures to abate emissions of NH 3 (Sutton and Howard, 2018).
There is still a lack of NH 3 and speciated monitoring of the inorganic gas and aerosol composition across the EU. An implementation of the DELTA ® approach across Europe would provide cost-efficient monitoring of the gas-and aerosol-phase pollutants for which reduction commitments are set out in Annex II to the NECD. Monitoring of NH 3 and the interacting acid gases and aerosols is needed to assess contributions of NH 3 to PM 2.5 , known to be harmful to human health. In addition, such monitoring will also provide the baseline and evidence against which any changes and potential recovery in ecosystem response to changes in emissions of the pollutant gases can be assessed, as required under Article 9 of the NECD. Issues such as human health impacts from fine ammonium aerosols will also drive policy decisions since controlling NH 3 should also reduce PM concentrations (Backes et al., 2016).
Data availability. Summary data are provided in the Supplement. The full dataset is available from http://www.nitroeurope.ceh.ac.uk/ .
Author contributions. YST coordinated the establishment of the networks, measurement, and collection of data with the support of several European laboratories. A large number of research institutes provided monitoring sites and local support for installation of equipment and carrying out the monthly exchange of samples. MAS conceived the NEU project and the DELTA ® network. IS helped with designing and building the low-voltage DELTA ® equipment. EN helped with network logistics and provided science advice. NvD helped with running proficiency testing schemes and inter-comparisons. UlrikD provided GIS support and science advice. CP, MJS, and UlricD facilitated and hosted the DELTA ® intercomparisons at their field sites. UlricD also sourced Rotenkamp bulk collectors for the project. JNC provided advice on bulk wet-deposition measurements and calculations. CRF, UlricD, SV, MM, HTU, and MJS led the chemical laboratories that shared the DELTA ® and wet-deposition measurements. DL developed the NEU database and provided support in data submission. Several of the authors contributed to measurements, network operations, and equipment and site maintenance. YST performed all the data collection and data analysis (including statistics) and wrote the manuscript, with input from all co-authors. MAS, MRH, CRF, Ul-ricD, MF, and JNC provided valuable advice on the interpretation of results and feedback on the manuscript.
Competing interests. The authors declare that they have no conflict of interest.