ORIGINAL RESEARCH article

Front. Environ. Health, 11 March 2024

Sec. Air Quality and Health

Volume 3 - 2024 | https://doi.org/10.3389/fenvh.2024.1249457

Monitoring ammonia concentrations in more than 10 stations in the Po Valley for the period 2007–2022 in relation to the evolution of different sources

  • CC

    C. Colombi 1*

  • LD

    L. D’Angelo 1,2

  • BB

    B. Biffi 1

  • EC

    E. Cuccia 1

  • UD

    U. Dal Santo 1

  • GL

    G. Lanzani 1

  • 1. Air Quality Department, Environmental Protection Agency of Lombardy Region (ARPA Lombardia), Milan, Italy

  • 2. Institute for Atmospheric and Environmental Sciences, Goethe-University Frankfurt, Frankfurt am Main, Germany

Abstract

Regarding secondary aerosols, in addition to the significant contribution of anthropogenic gases such as NOx and SO2, atmospheric ammonia (NH3) plays a crucial role as the primary basic gaseous species capable of neutralizing acidic compounds. This acid–base reaction is responsible for the formation of ammonium nitrate (NH4NO3), constituting up to 60% of PM10 within the Po River basin in Italy. Ion chromatographic analyses performed on offline samples indicate that this secondary inorganic species exhibits minimal concentration variability over the Po Valley because of limited air circulation due to orography and mesoscale air circulation. Therefore, investigating gaseous precursors becomes crucial. From the northern to the southern part of Lombardy—the region at the center of the basin—NH3 emission amounts account for 2.5, 11.1, and 27.7 t/y/km2, mainly due to agriculture and livestock activities (∼97%). To study NH3 temporal and spatial variability, the Environmental Protection Agency of Lombardy Region has been monitoring NH3 concentrations across its territory since 2007, with 10 active monitoring sites. Annual and seasonal cycles are presented, along with a focus on different stations, including urban, low-mountain background, high-impact livestock, and rural background, highlighting the impact of various sources. Measurements indicate that within the Po basin, NH3 concentrations can reach up to 700 µg/m3 (as an hourly average) in proximity to the main gaseous NH3 source. Instrument intercomparisons among online monitors and passive vials, as well as different online monitors, are presented. Therefore, this paper provides crucial data to understand the formation of secondary inorganic aerosols in one of the most important hotspot sites for air pollution.

1 Introduction

Despite its substantial size, the Po Valley, located in northern Italy, is essentially a semi-closed basin surrounded by the Alps, the Apennines, and the Adriatic Sea. The orographic barriers, including the Alps and the Apennines, restrict the entry of external air and precipitation events, hindering the mixing of the planetary boundary layer. These conditions are especially severe during wintertime when mixing layer heights are extremely low (). Thermal inversion exacerbates the situation, trapping anthropogenic emissions close to the surface (). Consequently, air quality in the region is significantly impacted by stagnant low-atmospheric conditions, coupled with high anthropogenic gases and particle emissions from the most industrialized and cultivated area in Italy. As a result, the Po Valley is among the most polluted areas in Europe ().

Data from the Environmental Protection Agency of the Lombardy Region (L-EPA) highlight ongoing concerns regarding fine particles (particulate matter, PM10 and PM2.5), nitrogen dioxide (NO2), and O3. Italy has recently faced legal consequences, being sentenced by the European Court of Justice () for systematically and persistently exceeding the limit values for concentrations of PM10. Thanks to efforts by public authorities over the years, concentrations of these pollutants have gradually decreased (). However, further reduction efforts are necessary, as even a lockdown during the COVID-19 pandemic in the spring of 2020 did not result in significant improvements (, ).

Regarding atmospheric particles, various studies (–) emphasize the significance of secondary inorganic aerosols (SIA) in mass contribution, including SO42−, NO3−, and NH4+, which can constitute up to 60% of the PM10 mass during wintertime in the Po Valley. Specifically, on days when the regulatory limit of 50 µg/m3 is exceeded, ammonium nitrate (NH4NO3) is a fraction that significantly varies ().

In this context, the gaseous precursors are well-known. Sulfur dioxide (SO2) is no longer a critical issue, thanks to the ban on sulfur-containing gasoline, which has reduced the resulting SO42− concentrations. Furthermore, advancements in combustion technology, along with policy restrictions, have had a positive impact on NOx emissions from combustions. However, atmospheric ammonia (NH3) is the primary gaseous species capable of neutralizing inorganic acidic compounds resulting from the photochemistry of the aforementioned gases, with well-known reaction pathways (). In cold and humid episodes, which often occur in the Po Valley, the equilibrium of the following reaction tends to shift to the particulate phase, leading to the formation of NH4NO3 (, ).where NH3 and HNO3 are in the gas phase and NH4NO3 is in the aerosol phase. Nitric acid is a sticky gaseous compound whose formation during daylight hours is due to the reactivity between NO2 and OH radicals, and from N2O5 heterogeneous hydrolysis, and reaction of the nitrate radical during nighttime.

Various methods have been developed in an attempt to quantify NH3, including the filter-pack method, denuder, and tunable diode laser absorption (), among others (). However, these methods have exhibited limited agreement among themselves and a relative sensitivity to atmospheric concentrations that can range from a few parts per billion (ppbv) to a few parts per trillion (pptv), respectively, in polluted and clean air (). For several decades, highly performant instruments like the Chemical Ionization Mass Spectrometer (CIMS) (, ) have been available. For instance, Hanke et al. () conducted a measurement campaign on Monte Cimone, at the southern border of the Po Valley, with such instruments, showing detection limits between 20 and 50 pptv. The authors demonstrated that, over a limited period (3 June to 6 July), nitric acid concentrations could vary significantly depending on the wind direction, ranging from 400 pptv under conditions of wind from Africa and dust to 1.2 ppbv with air from the boundary layer of NW-Europe. Although CIMSs have proven suitable for such measurements, they currently pose challenges in terms of management and maintenance and are unsuitable for a branched and continuous network.

Ammonia is also a sticky gas, primarily originating from various sources, including wastewater treatment plants, coal combustion, solid waste incineration, vehicular exhaust, biomass burning, fertilizer production, and emissions from humans and pets (–). A geospatial visualization of NH3 concentration along the vertical column is proposed in the study by Clarisse et al. (), retrieved by the Infrared Atmospheric Sounding Interferometer sensor (installed on the meteorological platform MetOp-A) observation. According to the Regional Inventory Emission database (), the major contributor to NH3 in the Po Valley is agriculture and livestock activities (∼97%). When Lombardy is divided longitudinally into three large areas, NH3 emissions amount to 2.5, 11.1, and 27.7 t/y/km2, decreasing from the north to south within the region. This is attributed to the catabolism of physiologic proteins, primarily ending in urea or uric acid and ultimately producing NH3, as summarized by Bussink and Oenema (). Animal slurry, rich in nitrogen-containing and organic compounds, serves as a natural fertilizer for agricultural soil. However, NH3 volatilization from manure storage, as well as during and after spreading activities, occurs (–). Despite existing techniques aimed at reducing NH3 volatilization and losses in livestock farming and agriculture, many studies suggest that emissions are increasing due to rising demands for food production ().

In addition to the availability of its gaseous precursors, i.e., NH3 and HNO3 (), the formation of NH4NO3 is influenced by the thermodynamic conditions of the air, particularly high temperatures and low humidity ().

The analysis of variability in NH3 concentrations is crucial for estimating the impacts of the agricultural sector on ecosystem health (, ) and understanding the dynamics of atmospheric particulate matter formation. This understanding is essential for supporting policymakers in identifying effective measures to reduce atmospheric aerosol burst events and the resulting exceedances of legal limits for the protection of health, to which NH4NO3 strongly contributes. Several studies have demonstrated the link between exposure to high and prolonged concentrations of fine PM and damage to the cardiovascular system (). In addition, specific compounds within particulate matter have known toxic effects on human health ().

While Pietrogrande et al. () recently demonstrated that secondary inorganic compounds formed with the contribution of atmospheric NH3 (e.g., NH4NO3) do not play a role in the oxidative potential of atmospheric aerosols, gaseous NH3 is known to participate in reactivity with volatile organic compounds, as shown by Bones et al. (), Daellenbach et al. (), and Wang et al. (). Despite these findings, the role of NH3 in the toxicity of resulting compounds remains uncertain to date. Some authors [e.g., Babar et al. (), Laskin et al. (), Li et al. (), Updyke et al. (), and Smith et al. ()] have identified the role of NH3 in the organic aerosol browning process, but its specific contribution to the toxicity of the resulting compounds requires further investigation.

In this study, we analyze a multi-year dataset of 1-h time resolution measurements of ambient air NH3 concentrations in various sites, investigating the spatial–temporal variability of NH3, its sources, and formation/transport processes in the middle of the Po Valley. The study presents annual and seasonal cycles, with a specific focus on different stations, including urban, low-mountain background, high-impact livestock, and rural background sites. This approach highlights the impact of various sources, such as air transport from the Po plain and traffic in Milan, the influence of livestock and agricultural activities in the southern part of the region, and the dominance of mixing layer height in the pre-alpine site.

The paper also provides additional data on NH4NO3 and NH3 concentrations, although an unambiguous relation between the two is not observed. The process of NH4NO3 formation in the Po Valley remains uncertain. Nonetheless, this paper offers a crucial overview of NH3 concentrations, laying the groundwork for further studies aiming to understand its role in PM bursts in this region.

2 Materials and method

Over the years, L-EPA has adjusted its monitoring network in response to updated regulations and with the objective of optimizing monitoring efficacy. Despite these modifications, the sampling sites conform to the European Regulation (). The selection of sites with NH3 monitors was a compromise between source locations and existing regular monitoring sites, as detailed in Section 2.1. While no specific regulations pertain to gaseous NH3 measurements, L-EPA has implemented Quality Assurance/Quality Control (QA/QC) standards and a maintenance intervention schedule similar to other gas monitors, as outlined in Section 2.3.

To gain valuable insights into the connection between NH3 and secondary inorganic aerosol compounds, PM10 has been sampled and analyzed using ionic chromatographic techniques, as discussed in Section 2.4.

2.1 Monitoring sites

Since 2007, L-EPA has equipped a total of 14 stations with NH3 monitors, as detailed in Table 1. The classification of sites adheres to the European Directive, with six located in urban or suburban areas (BG, Co, CR-FbF, MI-PA, PV, and SnB), three in proximity to intensive livestock activities (CdC, Pi, and Be), three in rural agricultural areas (CR-GB, MV, and SKI), and one at an elevation of about 1,200 m in a grazing land (Mo). The site labeled MB-P is situated in the middle of one of the largest fenced parks in Europe, close to a small livestock farm within the park and police horse stables. The locations of these stations, which have been part of the NH3 monitoring network over the years, are depicted in Figure 1.

Table 1

Site nameSite short nameLongitude (WGS84)Latitude (WGS84)Altitude (m asl)Measuring periodNH3-monitor (technique)
Bergamo (UB)BG9.64366245.69103823219/02/2021—OngoingAPI-M201E (Chem.)
Bertonico (RB)Be9.6662575845.233496896504/03/2009—OngoingAPI-M201E (Chem.)
Colico (SB)Co9.39112346.12408722918/06/2014—OngoingAPI-M201E (Chem.)
Corte de Cortesi (RB)CdC10.0062052845.278489975701/01/2007—OngoingAPI-M201E (Chem.)
Cremona—Via Fatebenefratelli (UB)CR-FbF10.0438476745.142543584310/09/2011—OngoingTEI-17i (Chem.)
Cremona—via Gerre Borghi (RB)CR-GB10.0692418145.109542753611/10/2011—OngoingTEI-17i (Chem.)
Milano—Pascal (UB)MI-PA9.23166745.47834712210/07/2007—OngoingAP2E-ProCeas (CRD);
API-201 (Chem.)
Moggio (RB)Mo9.4975426145.912790971,19717/03/2007—03/03/2021(Chem.)
Monza Parco (SB)MP9.27615645.60307218131/03/2013—08/07/2019(Chem.)
Motta Visconti (SB)MV8.99126745.28077810013/07/2021—OngoingAP2E-ProCeas (CRD)
Pavia—via Folperti (UB)Pa9.1646491945.194682317712/12/2013—OngoingTEI-17i (Chem.)
Piadena (SB)Pi1037982404512980303209/07/2013—31/07/2018(Chem.)
Sannazzaro de’ Burgondi (UI)SdB8.9041874245.102774648518/10/2013—OngoingEnvea AC32e + CNH3 (Chem.)
Schivenoglia (RB)SKI11.0760972845.016880671616/02/2013—OngoingTEI-17i (Chem.)

Information on the L-EPA air quality network with regard to atmospheric NH3 monitoring.

U, urban; R, rural; I, industrial; T, traffic; B, background; R, remote; Chem., chemiluminescence; CRD, cavity ring down.

Figure 1

This paper focuses on four case studies involving the following sites: MI-PA[UB] (an urban background station in Milan, the largest city in the middle of the Po Valley), CdC[farmRB] (a site near a swine farm where commonly elevated concentrations are detected), SKI[RB] (under investigation in the LIFE Prepair project about aerosol chemical composition and source apportionment, along with MI-PA[UB] and other Po Valley sites), and Mo[Rem] (noteworthy for its position often above the mixing layer height). Acronyms in square brackets, denoting the site types, are provided alongside the site names. Refer to the caption in Table 1 for a complete explanation of the acronyms. For the Mo station, “Rem” is used to emphasize its more remote location compared to SKI. Similarly, CdC[farmRB] indicates the proximity of this site to a swine livestock facility.

2.2 Ammonia monitors

To detect ambient air concentrations of NH3, the L-EPA network is mainly equipped with chemiluminescence-based monitors (, ). Because chemiluminescence is identified as the reference method for NOx [UNI ()], this technique is widely used in air quality networks.

However, over the years, L-EPA has improved its instrumentation with additional equipment such as cavity ring-down spectroscopy [CRDS (, )]. Although every monitoring station listed in Table 1 is equipped with a chemiluminescence-based monitor, only MI-PA[UB] and MV[SB] allow a comparison between these two techniques (Section 2.3).

2.3 Quality assurance/quality control

Minimizing the error of measurements and correctly determining the uncertainty associated with the measurement means reducing both the stochastic and deterministic components; this is guaranteed by the application of QA/QC procedures, provided by 2008/50/CE and by the Italian ministerial decree (DM) of 30th March 2017 (). In particular, QC activities consist of a series of procedural actions implemented by the L-EPA technicians according to the indications of the aforementioned DM and help ensure the quality control of the measured data in compliance with the quality objectives of European legislation. QA procedures consist of a series of scheduled activities through blind testing and performance auditing.

L-EPA applies QC procedures also to NH3 analyzers through scheduled preventive maintenance activities, other than regular calibrations and zero and span checks. NH3 measurements are affected by several positive artifacts, and a major source of interference is the presence of high and variable water vapor in the ambient air (). As reported in the literature (, ) and, in particular, at low concentrations, other factors identified for affecting the measurements are the inlet design, the filter material and aging, and the quality of calibration standards. Based on own experience and knowledge of literature, L-EPA has developed further internal procedures. For that, in addition to the scheduled maintenance activities, various precautions have been taken during the installation of the instruments to minimize the interferers and maximize the efficiency of the measurement. Some of them are (i) all pneumatic connections must be made with chromatography-grade (passivated) stainless steel tubing, glass, Teflon®; (ii) it is preferred to avoid filters on the sample line or, at least, to use a PTFE filter; (iii) use of a siliconized glass fiber sleeve to minimize the condensation of water; (iv) short pneumatic connections are recommended; (v) diluters with capillaries (no Mass Flow Controller) for the production of the NH3 calibration sample are to be preferred; (vi) oscillations of the ambient working temperature have to stay within ±5°C; (vii) it is recommended to use three different lines for the zero sample and for the NO/GPT and NH3 calibration samples (manual switching); (viii) it is necessary to wait at least 12 h for the stabilization of the NH3 sample (rise time = 90% at 5 min). Finally, the calibration must be carried out in equipped laboratories; the CRDS also are checked in the laboratory with sample gas readings.

About QA procedures, several intercomparison campaigns have been conducted over the years to confirm the reliability of the measurements.

Using the reference method for analysis (–), radially symmetric diffusive samplers have been used compared to the online measurements (three in parallel for reproducibility). The use of diffusive or passive samplers for the study of air pollution in environments of work has been known since the 1970s (, ). Considering the high variability of the passive samplers used since 2007, the inherent variability of the method, and dependence on atmospheric conditions, the results shown in Supplementary Table S1 are to be considered acceptable. Within the same technology, i.e., chemiluminescence, tools from different brands have been also compared. In fact, despite being based on the same measurement principle, various manufacturers adopt small technological differences. These results are summarized in Supplementary Table S1.

In Supplementary Figure S1, a comparison is shown between the data collected from a CRDS monitor (Proceas AP2E) and a chemiluminescence monitor (TEI-17i) near an agricultural field. The time series, with a 10-min resolution in the graph, highlights the good agreement between the two instruments. The comparison campaign between the two technologies consists of approximately 7,900 h of measurement. The results of the linear regression are indicated in the text box of the figure.

2.4 PM analysis

Starting from 2013 for MI-PA[UB] and 2018 for SKI[RB], a time series for PM10 chemical composition is available. Aerosol samples are collected daily on quartz fiber filters (Pall TissuQuartz®, Ø = 47 mm) and Teflon filters (Pall, Ø = 47 mm, PMP ring, 2 µm porosity) by means of gravimetric samplers [Skypost PM, TCR-TECORA, Cogliate (MB), Italy, or Lifetek PMS, Megasystem, Bareggio (MI)] equipped with a PM10 sampling head (1 m3/h, 24 h). After collection, filters are stored in the darkness and at low temperature to prevent photochemical reactions and compound volatilizations. A punch of 1.5 cm2 is removed from each filter, and the water-soluble compounds are extracted through ultra-pure water (Sartorius® Arium Mini, resistivity 18.2 MΩ) in a sonic bath (20 min). The resultant solution is then filtrated (Nylon or PTFE Syringe Filter—pore size 0.45 µm) and injected into an ion chromatographer (Metrohm 930 and 881) for the determination of anions (Cl−, NO2−, Br−, NO3−, PO43−, SO42−) and cations (Na+, NH4+, K+, Mg2+, Ca2+). Another punch of 1.5 cm2 is removed from each filter for the determination of the carbonaceous fraction by means of thermal–optical analysis (Sunset Laboratory Inc., Tigard, OR, USA), according to the NIOSH-like and EUSAAR-2 protocols. Teflon filters are used for elemental composition determination by dispersive x-Ray fluorescence (ED-XRF, Epsilon 4, Malvern Panalytical) for elements with atomic number Z higher than 11 (Mg, Al, Si, P, S, Cl, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Cu, Zn, As, Br, Rb, Cd, Pb, Sr, Sn, Sb, Ba). Starting from 2013, in MI-PA[UB], hourly Black Carbon measurements with a Multi-Angle Absorption Photometer (MAAP, Thermo Scientific Model 5012, PM2.5 cutoff, no dryer) that measures the aerosol absorption coefficient at a wavelength of 637 nm are also available (), and from it, the BC is computed considering the deposit area and sampling air flow and using a mass-specific absorption coefficient of 6.6 m2/g. In addition, during the latter part of winter 2022 (from 22 February to 15 March), a higher time resolution sampling campaign was performed to observe the variability of secondary inorganic compounds and the role of gaseous NH3 at various sites in Lombardy. For this reason, PM10 samples were collected (Pall TissuQuartz®, Ø = 47 mm, 1 m3/h, 6 h) at six sites of the L-EPA air quality network and equipped with NH3 monitors, i.e., MI-PA[UB], SKI[RB], CdC[farmRB], MV[RB], Be[RB], and SdB[UI].

3 Results

In this section, data concerning NH3 concentrations are presented. An overview of the data collected at all the stations is provided in Table 2.

Table 2

Statistics (µg/m3)BGBeCoCdCCR-FBFCR-GBMI-PAMoMBMVPaPiSdBSKI
Average
(±CI)		7.9 (±0.0)31.6 (±0.2)4.3 (±0.0)53.5 (±1.1)8.2 (±0.1)14.7 (±0.2)9.5 (±0.1)2.8 (±0.0)9.1 (±0.1)13.9 (±0)7.8 (±0.0)22.7 (±0.1)8.3
(±0.0)		16.1
(±0.2)
Geometric average5.925.42.937.45.310.57.51.16.612.85.319.85.912.8
Range0.1–46.10.1–432.70.1–58.60.1–7080.1–138.90.1–6430.1–182.90.1–36.80.1–238.22.0–152.80.1–850.1–463.10.1–660.1–739.7
10th–90th
percentile		1.9–13.811.6–59.80.7–8.313.2–109.81.4–16.94.2–26.92.9–16.30.1–6.82.5–17.87.7–21.31.7–15.710.4–38.22.2–15.36.0–27.5
Average February to April10.735.35.555.29.315.410.12.210.516.77.230.38.715.9
Average
May to July		8.224.13.952.179.910.13.75.812.67.613.36.632.1
Average
September to November		5.7363.566.39.523.69.33.110.611.88.22710.115.7
Average
Other periods		6.125.54.932.56.49.88.10.98.88.96.1227.514.2

Statistics concerning ammonia measurements based on hourly average native data since 2007.

CI, confidence interval at 95% significance level.

The annual cycle of average NH3 concentrations is depicted in Figure 2, along with the maximum and minimum values, the 75th and 25th percentiles, and the number of stations measuring NH3 for each year. It is evident that the average concentration remains almost constant over the years. The overall average is impacted not only by rural stations but also by the one in a remote area, which exhibits very low concentrations (Supplementary Figure S2). The maximum values are always attributed to the rural station close to husbandry activity, described later (Section 3.3). We speculate that the decrease in maximum values over the years is due to an improvement in technologies used in zootechnical and agricultural activities or a decrease in the surrounding agricultural activities. This decrease is coherent with INEMAR's databases, which report a decline in NH3 emissions in this area from 101,779 t/y in 2014 to 90,727 t/y in 2019.

Figure 2

The variability of the four case studies is described in further detail, as they exhibit characteristics representative of typical monitoring site classes. As mentioned earlier, MI-PA[UB] is a background station placed in an urban background area of Milan, the biggest city in the Po Valley. The nearest agricultural and livestock activities are approximately 5 km from the station. Therefore, MI-PA[UB] can be considered a background site also concerning the main NH3 sources. Conversely, the CdC[farmRB] site is a few tens of meters away from a swine farm and within an area primarily intended for agriculture. For this reason, this station could be considered an “industrial site” in relation to NH3 sources, especially when considering the presence of husbandry activity as an anthropic (industrial) influence. In a completely different context, the SKI[RB] station is placed in an agricultural area but with few animal husbandries in the surroundings. Therefore, it can be defined as an agricultural background site in the same manner that MI-PA[UB] is an urban background site. Finally, Mo[Rem] is located in an isolated small pre-alpine valley at 1,200 m a.s.l. Some nearby areas are sporadically and occasionally used as grazing land for small herds of bovines. For this reason, Mo[Rem] is rarely affected by NH3 emission, and the station may be regarded as a rural background and remotely located.

Ammonia concentrations measured at the four stations described above are presented in detail below, focusing on the cycle over the years.

3.1 Annual cycle

A valuable overview of NH3 concentrations is obtained by calculating the average daily concentrations (as the geometric mean) across the years for each day of the year. Figure 3 illustrates the annual cycle for the four selected monitoring sites as case studies, with dark lines representing the geometric mean. Generally, the arithmetic mean is more susceptible to spikes associated with local sources, whereas the geometric mean is more effective at representing background conditions ().

Figure 3

Ammonia concentrations measured at the MI-PA[UB] station are depicted in Figure 3A, where values of the geometric mean comparable to the literature for measurements in urban areas (, ) are presented. Although the average annual cycle does not show significant variability throughout the year (the 10th–90th percentiles range is between 3 and 40 µg/m3, with a geometric mean value of 7.5 µg/m3), the measurements suggest three periods with an increase in NH3 concentrations: late winter/beginning of spring, late springtime, and autumn. At the SKI[RB] station (Figure 3B), as for MI-PA[UB], the time series shows higher values in three periods, albeit slightly different compared to the previously described site: winter/early springtime, from June to August, and from mid-October to early November. The growth in July of the annual cycle is dominated by high peaks: this is due to very high hourly average concentrations (up to 740 µg/m3) observed every second year since 2018. L-EPA technicians verified that fertilization operations of nearby agricultural fields were in progress during those events. Nevertheless, the geometric annual mean for the NH3 level amounts to 13 µg/m3.

The highest concentrations are recorded at the CdC[farmRB] station, as expected. At this site, the yearly geometric mean level reaches 37 µg/m3. Although the maximum concentration measured on an hourly basis is 708 µg/m3, comparable to that measured in SKI[RB], the CdC[farmRB] site significantly differs from SKI[RB]. This difference is evident when considering the 90th percentile value of the two sites (28 µg/m3 at SKI[RB] compared to 110 µg/m3 at CdC[farmRB]), as reported in Table 2. In addition, the annual cycle at CdC[farmRB] station shows two main periods of rising concentrations: from February to the beginning of April and from July to November, the latter one is preceded by 2 months of a modest increase in concentrations. At the Mo[Rem] site, the annual cycle shows a distinctly different behavior. Broadly speaking, NH3 concentrations are higher during summertime and very low (sometimes even below the detection limit) during winter, and the overall shape is an upside-down “U.” However, data collected at the Mo[Rem] station from 2007 to 2020 suggest even here the influence of NH3 in three distinct periods. Rises in NH3 concentration pattern are observed from February to April, during June and July, and during September and October. Nevertheless, the annual mean concentration is about 1 µg/m3.

Averaged concentrations during three representative periods, i.e., from mid-February to mid-April, from mid-May to mid-June, and from mid-September to mid-November, are calculated and shown for each site. Further statistics on an hourly basis are summarized in Table 2.

In Figure 4, the atmospheric NH3 concentrations are presented as weekly arithmetic means for the four sites previously described. The weekly time resolution allows appreciation of the variability of the shown data, avoiding excessive scatter. On the right axes, the maximum hourly data in the specific average week are also shown. Although two high values were detected in MI-PA[UB] in 2009 and 2012 (up to 180 µg/m3 as maximum hourly), NH3 weekly average concentrations typically remain within a range between 3 and 30 µg/m3. As shown with the dashed red line in the figure, the time series indicates a negligible cycle over the years for MI-PA[UB]. As previously mentioned, SKI[RB] confirms to be characterized by low variability as well, besides rare events connected to soil fertilization (the direct link between fertilization and high levels is shown in Section 3.3). For this site, a negligible cycle is observed too, even though a slight increase of the interannual trend is shown in this case. A more appreciable reduction of NH3 interannual trend levels is suggested by the time series of CdC[farmRB] site, also displayed by the mean annual cycle (Figure 5). The Mo[Rem] site, instead, clearly shows that values above the detection limit are mainly detected when the warmer season starts and again decrease after it.

Figure 4

Figure 5

In general, a clear interannual trend is not evident; on the whole, the NH3 values measured by the L-EPA Air Quality Network seem quite stable, as shown in Figures 4, 5. The data, therefore, suggest not great changes in the activities emitting NH3 into the atmosphere, over the years, with a few exceptions. As already mentioned, the reduction in annual average concentrations measured at CdC[farmRB], which affects the whole maximum hourly dataset depicted in Figure 2, can be attributed to local improvements in waste management practices and in a change in the number of pig livestock. On the other hand, the monthly average cycle (Figure 6) confirms the two main periods of rising concentrations, with the highest values in March and August to September. At SKI[RB], concentrations exhibit a slight positive pattern, although strongly influenced by high concentrations observed in July 2018, 2020, and 2022, as also highlighted in Figure 4. MI-PA[UB] is influenced by transport events that are unlikely to result in a similar increase in values as observed in locations near direct emissions. Mo[Rem] shows a completely different graph than other sites, showing a bell-trend: the highest concentrations are measured during warmer periods reaching the maximum value in June.

Figure 6

3.2 Impact of the agricultural sector on NH3 ambient air

As demonstrated in the previous section, specific periods of the year are typically marked by higher NH3 concentrations, particularly in agricultural areas. As suggested by emission inventory databases, these variations could indicate that some specific agricultural and zootechnical activities contribute significantly to NH3 emissions. For this reason, starting from 2017, an examination of the impact of the agricultural sector on NH3 concentrations has been conducted. In this context, four monitoring campaigns at different animal husbandry activities along the Lombardy region have been carried out to compare different fertilization techniques on agricultural fields. The considered techniques were (1) surficial spreading, (2) direct injection within the soil, (3) nebulization by means of a pivot 20 cm above the ground level of N-abate and micro-filtrated slurry, and (4) fertigation of N-abate and micro-filtrated slurry. Among these, the first two are the most used techniques in the Po Valley. In this section, we briefly describe the results of the first technique. The monitoring campaign was conducted in an agricultural area between Milan and Bergamo cities in September 2018. The comparison among the different techniques will be discussed in a work that is in progress.

Ammonia concentrations were measured by a chemiluminescence monitor installed on a mobile laboratory positioned near the border of the agricultural field (Figure 7A). Data were acquired with a 1-min time resolution starting from 3 days before the fertilization event, dated 26 September 2018. During this time range, observations displayed NH3 concentrations lower than 20 µg/m3. The superficial spreading fertilization used in this area consists of a rotating disk, which spreads the fertilizer on the soil from a height approximately 3 m above ground level. Due to the emission and strong volatilization caused by the spreading technique, atmospheric NH3 concentrations increased from the above-mentioned value up to 1,300 µg/m3. Conversely, when the monitor was upwind of the fertilized field, concentrations were lower than 50 µg/m3. In this regard, the polar plot in Figure 7B reports the probability that the measured concentrations above the 50th percentile come from the target fertilized area. Ammonia concentrations rapidly decreased during the following days, thanks to the regulations that order such fertilizers to be buried within 48 h after spreading. This led the concentrations to decrease below 100 µg/m3 on 29 September.

Figure 7

3.3 Secondary inorganic compounds in the Po Valley

As already mentioned, the Po Valley is a hotspot regarding PM concentrations. Many studies focused their attention on the chemical characterization of atmospheric particles in this basin and provided averaged PM composition both on an annual (–) and seasonal (–) basis. Following the results shown in the Life PrepAIR—Interim Report (), in Milan, the averaged composition of PM10 in the last 10 years is due to SIA (34%), organic carbon (OC) (20%), crustal matter (12%), elemental carbon (EC) (4%), and other trace elements (2%). Carbon compounds remain quite constant in percentage from summer to winter, while SIA increases up to 38% in the cold season. By focusing solely on the results of chemical analyses associated with days surpassing the 50 µg/m3 limit, the influence of NH4NO3 becomes distinctly evident. In instances where the limit is exceeded, NH4NO3 constitutes 28 ± 10% of the PM10, contrasting with the 17 ± 13% observed in cases below the limit. Beyond a mere distinction between below and above the limit, it is evident that NH4NO3 plays a crucial role in determining PM10 levels in the Po Valley. Thus, the study of NH3 concentrations in relation to its sources and meteorological phenomena is important to understand its role in secondary inorganic compounds formation. In this section, data regarding secondary inorganic compounds are presented. Figure 8 shows the results of the chromatographic analysis performed on about 4,000 PM10 samples. Daily concentrations were aggregated in monthly time resolution and then used to highlight the annual cycle. The ammonium sulfate compound shows very low variability, ranging from 2.0 to 3.5 µg/m3 and no seasonality (Figure 8C). It is worth noticing that MI-PA[UB] and SKI[RB] show the same annual cycle and the same absolute value of this compound.

Figure 8

On the contrary, NH4NO3 is found to contribute up to 40% of PM10 mass on a monthly basis in colder periods, which can reach up to 60% of PM10 on a daily average basis (Figures 8D,E).

To identify a clearer relationship within the non-linear (, ) NH3–NOx–NH4NO3 system, further investigation was carried out to observe the cycle of NH4NO3 during one of the periods with the highest probability of elevated NH3 concentrations, based on the levels of the previous years. The investigation has been conducted through an intensive monitoring campaign between February and March 2022 at six different sites, as explained in Section 2.4. Figure 9 displays the variability of two parameters, gaseous NH3 and NH4NO3 on PM10, over a 6-h sampling period. Gaseous NH3 concentrations exhibit high variability among the monitoring sites (Figure 9A), suggesting that the emissive source and its distance from the measurement site have the greatest impact on the determination of atmospheric concentrations. Within the considered period, the CdC[farmRB] site shows the highest values (57 ± 35 µg/m3) followed by Be[RB] (47 ± 23 µg/m3), whereas the lowest ones are at MI-PA[UB] (9 ± 2 µg/m3) and SdB[UI] (6 ± 3 µg/m3) sites. On the contrary, NH4NO3 (Figure 9B) shows an extremely low variability. The mean concentrations are in a range from 7 ± 6 µg/m3 [at SKI(RB)] to 11 ± 7 µg/m3 [CdC(farmRB)]. Figure 10 allows us to observe the NH3 and NH4NO3 cycles during the intensive campaign in CdC[farmRB]. Although the determination coefficient is very low (R2 = 0.2), the data point out a common pattern in specific events. This pattern is observed across all individual measurement sites, indicating a potential direct relation between these two parameters.

Figure 9

Figure 10

4 Discussion

In this section, we delve into the results of NH3 monitoring at the four selected sites, examining their cycle based on their sources. Subsequently, we compare the NH3 cycle with NH4NO3 concentrations and critically analyses them based on the results of chemical analyses of PM10 samples.

4.1 Concentrations of NH3 and their sources

The selected locations enable us to investigate how concentrations of this gaseous compound may vary depending on the proximity of the primary source identified by the emission inventory. Figure 11’s polar plots illustrate the relationship between NH3 levels, wind speed, and direction.

Figure 11

SKI[RB] and CdC[farmRB] sites are surrounded by agricultural fields. Source apportionment analysis at SKI[RB] suggests that traffic and/or industrial activities collectively contribute about 10% to the air quality impact. Monitoring at CdC[farmRB] in 2014 revealed low concentrations of SO2 and NOx, proxies for industrial and combustion sources. The regional emission inventory (INEMAR) emphasizes that these two sources contribute less than 5% of total NH3 emissions in the Schivenoglia and Corte de’ Cortesi administrative areas. Conversely, NH3 emissions are primarily attributed to the agricultural sector, particularly swine husbandries, accounting for over 75%. Thus, SKI[RB] and CdC[farmRB] sites reinforce the idea that NH3 concentration levels are strongly influenced by the proximity of the emissive source.

At CdC[farmRB], the polar plot indicates that the highest values are measured under no-wind or low-wind conditions, supporting a local source, such as animal husbandries. As wind speed increases, concentrations rapidly decline, though remaining high due to a persistent high background condition. Supplementary Figure S3 further supports this evidence.

The conditional functional plot, generated by eliminating extreme cases, enables the determination of the source of concentrations below the 80th percentile (Supplementary Figure S3A). This allows observation of the homogeneity of the probability of such concentrations across all wind directions and speeds, indicating that the calculated values are consistent with the local agricultural background in the Schivenoglia area. On the other hand, a local source is also suggested to affect the SKI[RB] site with high NH3 levels: the polar plot (Supplementary Figure S3B) displays the highest concentrations in low-wind conditions but also highlights that transport events can increase the background values. Limiting the observation only to the data over the 90th percentile, or 27 µg/m3, confirms that the source is indeed local, but also the concentrations in these cases are indicative of a transport process with wind speeds below 4 m/s. The little variability in the direction of origin of such concentrations suggests that the monitoring station detects activities carried out at west–northwest directions.

Figures 3, 4 demonstrate that hourly NH3 concentrations for SKI[RB] and CdC[farmRB] can reach 700 g/m3. This aligns with the findings of the campaign described in Section 3.2: soil fertilization by means of animal manure or slurry strongly affects the detected levels of NH3. Table 2 reports the arithmetical averages calculated on an hourly time resolution data in four periods, i.e., mid-February to mid-April, mid-May to mid-July, mid-September to mid-November, and the remaining periods of the year. Based on the information about agricultural practices, the observations found an explanation referring to the detected values. These are the periods in which soil fertilization occurs, although each agricultural area has its own particularity. In addition, meteorological parameters (such as temperature) affect resulting NH3 concentrations due to volatilization variability. Both these facts explain the lowest values detected for SKI[RB] and CdC[farmRB] during December and January (and for the entire L-EPA network): in these 2 months, spreading procedures are prohibited or heavily limited within the entire Po basin, and NH3 volatilization is prevented by the lower temperature that, at the same time, enhances gas-to-particles partitioning in the aerosol phase.

The Mo[Rem] site, as depicted in Figure 6D, shows distinct variability compared to other sites. In particular, the highest NH3 concentrations are measured during warmer periods and the lowest during colder ones. This observed pattern is attributed to the mixing layer height (Supplementary Figure S4): considering its location, i.e., about 1,200 m a.s.l, the site is always above the mixing layer from October to April, with the exception of some isolated events where meteorological conditions favor convective motions vertically, raising the mixing height to levels compatible with that of the site. Otherwise, during this period, the vertical stability of the atmosphere hinders compounds emitted by the Po basin from reaching altitudes above the mixing layer height. As a result, NH3 concentrations are frequently below the instrument's quantification limit. On the contrary, the maximum monthly average is in June (Figure 6D) when the site is within the mixing layer, allowing for the measurement of concentrations of air masses transported from the Po Basin (, –, , ).

MI-PA[UB] (Figure 11A) demonstrates that the main source of atmospheric NH3 comes from the plain located ESE of the site, whereas the lowest values are detected in cases of wind reinforcement that cleans the atmosphere. The graph also indicates a slight increase in concentrations under light or gentle breezes (2–4 m/s) from the west, suggesting that these increases may originate from livestock activities on the opposite side of the city of Milan or from vehicular traffic.

It is known that in urban sites, a significant source of NH3 is vehicular traffic. NH3 emissions from gasoline vehicles equipped with a three-way catalyst (TWC) are an important source of NH3 in areas with heavy traffic (, , ), since it is generated as a side product in the NOx reduction process (). Furthermore, the recent introduction of the selective catalytic reduction (SCR) system with the addition of urea or NH3 in heavy-duty vehicles (HDV), and mandatory since 2016 for Euro 5 and Euro 6 vehicles, resulted in increased NH3 emissions from traffic (), which needs further investigation.

To verify that a fraction of the NH3 at the MI-PA[UB] site may be influenced by this source, the concentrations measured in December and January were compared with Black Carbon values (Section 2.4; Supplementary Figure S5), a well-known marker for vehicular combustion (). In this period, regional limitations ban spreading activities according to the Nitrates Directive [and its Italian regulatory transposition DM 5046/2016 ()]. The resulting correlation demonstrates a good level of agreement between these two parameters (R2adj = 0.661), confirming a relative contribution from this source. With the same hypotheses, an attempt was made to verify whether the addition of urea had repercussions on NH3 concentrations in Milan. The historical pattern (Figure 4) does not show an increase. It is necessary to remember that the restrictions introduced during the COVID-19 pandemic have certainly had a positive effect in reducing the impact of vehicular traffic. However, similar to before, taking into consideration the months of December and January when there is a good correlation between atmospheric NH3 and traffic, and excluding only December 2020 and January 2021, which suffered the restrictions of the second wave of COVID-19, the NH3 concentrations were averaged over the two periods: 2007–2016 and 2017–2022 with 6.4 and 11.8 µg/m3, respectively. This result cannot be used as verification of an increase in NH3 concentrations due to the addition of urea but suggests the possibility of further investigations.

4.2 Ammonia vs. ammonium nitrate in the Po Valley

The representative cycle of the six chosen monitoring sites, as detailed in Section 3.3, appears to confirm a qualitative correlation between gaseous NH3 and NH4NO3 in the aerosol phase. However, the low coefficient of determination equally indicates that the relationship between the two compounds is not linear (Section 3.3). The intensive period spans 21 days (from late February to mid-March) and is characterized by meteorological stability, with no rainfall, an average wind speed of 2 m/s, and atmospheric pressure of 1,014 hPa. In addition, thermodynamic air parameters, including humidity and temperature, are at levels that could favor the partitioning of NH4NO3 into the particulate phase.

The variability of NH4NO3 during this period (Figure 9B) highlights three events with a significant increase in concentrations occurring on (I) 25 February, (II) 4 March, and (III) 10–11 March. Compared to an average over the period of 22% of NH4NO3 in PM10, the contributions were 37%, 42%, and 33%, respectively. These three periods were investigated by analyzing the PM10 and NOx patterns, the mixing layer height, and relative humidity (Figure 12). The graphs present cycles as moving averages, helping to smooth out the time series curve by computing the average of all data points in a fixed-length window. It can be observed that in the first and third episodes, the increase in concentrations was followed by a significant decrease in the mixing layer height, reaching the lowest value in the campaign between 10 and 11 March. The second episode of increased NH4NO3 occurs instead in conjunction with an increase in the mixing layer height, which remained quite high even in the earlier time slots. The relative humidity also increases together with NH4NO3, while precursors decrease. These factors suggest the possible occurrence of an NH4NO3 formation event.

Figure 12

On the other hand, Figure 13 shows the average daily concentration cycle for both compounds for 1 year (2019 was chosen as an example). During the cold seasons, particularly for many peak episodes in spring and fall, it can be noticed that when NH3 increases, NH4NO3 also increases. However, several events contradict a direct cause–effect connection: in warmer periods, the condensation into NH4NO3 is inhibited by the temperature, which favors the evaporation of the nitrates. Moreover, the formation of NH4NO3 is also caused by accumulation phenomena, which occur frequently in the Po Valley, and by combustion sources. Ammonium nitrate formation, and secondary inorganic aerosol formation in general, is a complex process influenced not only by the concentration of its precursors but also by thermodynamics and meteorological conditions (, ), as already discussed. Nenes et al. () and Thunis et al. () published two different works about the PM2.5 response to HNO3 or NOx and NH3 emissions changes. Their modeling approaches converge to similar results, i.e., a variation in one of the two gaseous precursors does not lead to a linear change in PM2.5 concentrations and, in some cases, it could lead to the opposite result.

Figure 13

Nevertheless, it is worth noting that the findings presented in this study are based on offline analyses conducted on PM10 samples. These methodologies can be prone to negative artifacts, which may significantly underestimate the semi-volatile component of particulate matter. Past studies by Minguillón et al. () and Poulain et al. () have effectively compared offline and online analyses, the latter utilizing data from the Aerosol Chemical Speciation Monitor (ACSM, Aerodyne Research Inc.). Importantly, within the scope of this study, it is essential to highlight that offline and online nitrate correlations exhibit slopes greater than 6 during the summer months. This suggests that offline results obtained during sampling under high-temperature conditions could consistently underestimate NH4NO3 concentrations, thereby reducing the observable relationship between concentrations of gaseous precursors and the resulting particulate phase.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in https://www.arpalombardia.it/temi-ambientali/aria/form-richiesta-dati-stazioni-fisse.

Author contributions

BB, LD, CC, EC and UD conceived and designed the study, acquired data and interpreted the results, and wrote the paper. CC, LD and BB analyzed the data. BB and UD supplied the meteorological observations. LD, UD, and CC produced the figures. LD and BB produced the tables. EC carried out the IC analyses and helped with the interpretation of the chemical speciation. GL and CC interpreted the data and reviewed and edited paper. All authors contributed to the article and approved the submitted version.

Funding

The intensive campaign reported in Section 3.3 was partially financed by the General Directorate of Agriculture of Region Lombardia.

Acknowledgments

The authors would like to thank all the technician colleagues who with their professionalism and effort make the measurement network efficient.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvh.2024.1249457/full#supplementary-material

References

Summary

Keywords

ammonia, animal husbandry, ammonium nitrate, Po Valley, PM, secondary aerosol, gaseous precursors

Citation

Colombi C, D’Angelo L, Biffi B, Cuccia E, Dal Santo U and Lanzani G (2024) Monitoring ammonia concentrations in more than 10 stations in the Po Valley for the period 2007–2022 in relation to the evolution of different sources. Front. Environ. Health 3:1249457. doi: 10.3389/fenvh.2024.1249457

Received

28 June 2023

Accepted

21 February 2024

Published

11 March 2024

Volume

3 - 2024

Edited by

Xavier Querol, Spanish National Research Council (CSIC), Spain

Reviewed by

Cristina Mangia, National Research Council (CNR), Italy

Sailesh Behera, Shiv Nadar University, India

Updates

Copyright

© 2024 Colombi, D'Angelo, Biffi, Cuccia, Dal Santo and Lanzani.

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: C. Colombi

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Share on WeChat

Scan with WeChat to share this article

Article metrics