# Non-linear response of PM 2.5 to changes in NO x and NH 3 emissions in the Po basin (Italy): consequences for air quality plans

Philippe Thunis, Alain Clappier, Matthias Beekmann, Jean Philippe Putaud, Cornelis Cuvelier, Jessie Madrazo, Alexander de Meij

► To cite this version:

Philippe Thunis, Alain Clappier, Matthias Beekmann, Jean Philippe Putaud, Cornelis Cuvelier, et al. Non-linear response of PM 2.5 to changes in NO x and NH 3 emissions in the Po basin (Italy): consequences for air quality plans. Atmospheric Chemistry and Physics, 2021, 21 (12), pp.9309-9327. 10.5194/acp-21-9309-2021. hal-03442074

HAL Id: hal-03442074

https://hal.science/hal-03442074

Submitted on 22 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Non-linear response of PM<sub>2.5</sub> to changes in NO<sub>x</sub> and NH<sub>3</sub> emissions in the Po basin (Italy): consequences for air quality plans

Philippe Thunis<sup>1</sup>, Alain Clappier<sup>2</sup>, Matthias Beekmann<sup>3</sup>, Jean Philippe Putaud<sup>1</sup>, Cornelis Cuvelier<sup>4</sup>, Jessie Madrazo<sup>5</sup>, Alexander de Meij<sup>6</sup>

Corresponding author: Philippe Thunis (philippe.thunis@ec.europa.eu)

**Abstract:** Air pollution is one of the main causes of damages to human health in Europe with an estimate of about 380 000 premature deaths per year in the EU28, as the result of exposure to fine particulate matter (PM<sub>2.5</sub>) only. In this work, we focus on one specific region in Europe, the Po basin, a region where chemical regimes are the most complex, showing important non-linear processes, especially those related to interactions between NO<sub>x</sub> and NH<sub>3</sub>. We analyse the sensitivity of PM<sub>2.5</sub> concentration to NO<sub>x</sub> and NH<sub>3</sub> emissions by means of a set of EMEP model simulations performed with different levels of emission reductions, from 25% up to a total switch-off of those emissions. Both single and combined precursor reduction scenarios are applied to determine the most efficient emission reduction strategies and quantify the interactions between NO<sub>x</sub> and NH<sub>3</sub> emission reductions. The results confirmed the peculiarity of secondary PM<sub>2.5</sub> formation in the Po basin, characterised by contrasting chemical regimes within distances of few (hundreds of) kilometres, as well as non-linear responses to emission reductions during wintertime. One of the striking results is the slight increase of the PM<sub>2.5</sub> concentration levels when NO<sub>x</sub> emission reductions are applied in NO<sub>x</sub>-rich areas, such as the surroundings of Bergamo. The increased oxidative capacity of the atmosphere is the cause of the increase of PM<sub>2.5</sub> induced by a reduction in NO<sub>x</sub> emission. This process can have contributed to the absence of significant PM<sub>2.5</sub> concentration decrease during the COVID-19 lockdowns in many European cities. It is important to account for this process when designing air quality plans, since it could well lead to transitionary increases in PM<sub>2.5</sub> at some locations in winter as NO<sub>x</sub> emission reduction measures are gradually implemented. While PM<sub>2.5</sub> chemical regimes, determined by the relative importance of the NO<sub>x</sub> vs. NH<sub>3</sub> responses to emission reductions, show large variations seasonally and spatially, they are not very sensitive to moderate (up to 50-60%) emission reductions. Beyond 25% emission reduction strength, responses of PM2.5 concentrations to NOx emission reductions become non-linear in certain areas of the Po basin mainly during wintertime.

**Keywords:** urban air pollution, air quality planning, non-linearity, chemical regimes

## 1. Introduction

Air pollution is one of the main causes of damages to human health in Europe with an estimate of about 380 000 premature deaths per year in the EU28, as the result of exposure to fine particulate matter (PM<sub>2.5</sub>) only (EEA, 2020). Many of the exceedances to the EU limit values occur in urban areas where most of the population is exposed.

<footnote n="1"><sup>1</sup> European Commission, Joint Research Centre, Ispra, Italy</footnote>

<footnote n="2"><sup>2</sup> Université de Strasbourg, Laboratoire Image Ville Environnement, Strasbourg, France</footnote>

<footnote n="3"><sup>3</sup> Laboratoire Interuniversitaire des Systèmes Atmosphériques, UMR CNRS 7583, Université Paris Est Créteil et Université de Paris, Institut Pierre Simon Laplace, Créteil, France</footnote>

<footnote n="4"><sup>4</sup> Ex European Commission, Joint Research Centre, Ispra, Italy</footnote>

<footnote n="5"><sup>5</sup> Signa Terre SA, Geneva, Switzerland</footnote>

<footnote n="6"><sup>6</sup> MetClim, Varese, Italy</footnote>
PM<sub>2.5</sub> is partly emitted directly (primary particles) and partly formed through photo-chemical reactions that involve gaseous precursors like SO<sub>x</sub>, NO<sub>x</sub>, NH<sub>3</sub> and non-methane volatile organic compounds (NMVOC) to form secondary inorganic and organic aerosol (SIA and SOA). The secondary fraction is often dominating the total concentration of particulate matter in urban areas as shown by e.g., Beekmann et al. (2015) for the Greater Paris region, De Meij et al. (2006, 2009) or Larsen et al. (2012) in northern Italy, hence the importance to understand the complex chemical processes that lead to its formation. In particular, it is key to identify the precursors involved in these reactions in order to target the right sectors of activity in air quality plans to effectively reduce pollution levels. According to the EDGAR estimates for 2015 for Italy, about 90% of the NH<sub>3</sub> is directly emitted in the atmosphere by the agriculture sector while SO<sub>x</sub> precursors are predominantly released by the energy production and use (industrial) sectors (EDGAR, 2020). For NO<sub>x</sub>, emissions are spread among various sectors, with transport (50%), industry (40%) and agriculture (4%) being the main ones. The gaseous precursors of secondary organic aerosols (SOA) include a vast range of NMVOCs among which biogenic terpenes and anthropogenic aromatics. The main sources of aromatics in Italy were in 2012 transport (58%), use of fuels and solvents (32%), and domestic heating (15%) (EDGAR, 2020).

Regarding SIA, early works used box models with thermodynamic schemes to address the sensitivity of ammonium nitrate and sulfate concentrations to gaseous NH<sub>3</sub>, NO<sub>x</sub>, and SO<sub>2</sub> emissions (Watson et al., 1994; Blanchard and Hidy, 2003; Pozzer et al. 2017; Guo et al. 2018; Nenes et al. 2020). These models were later on integrated into chemical transport models (CTM), in particular to address the benefit of additional NH<sub>3</sub> emissions reductions in addition to already ongoing SO<sub>2</sub> and NO<sub>x</sub> emission reductions. For North America, Makar et al. (2009) simulated with a regional CTM that a 30% reduction of ammonia emissions would lead to about 1 µg m<sup>-3</sup> reduction in PM<sub>2.5</sub>. For Europe, Bessagnet et al. (2014) simulated the effects of a 30% NH<sub>3</sub> emission reduction in addition to those foreseen by the Gothenburg Protocol for 2030, and found that the G-ratio defined as the ratio between free ammonia and total nitrate (Ansari and Pandis, 1998) was a good predictor for the efficiency of NH<sub>3</sub> reductions on SIA concentrations. These sensitivities to emission reductions are often governed by complex chemical mechanisms. A well-known phenomenon is the release of free ammonia as a result of decreased SO<sub>2</sub> emissions and sulfate formation, which allows for the formation of additional particulate nitrate, as described for example by Blanchard and Hidy (2003) and Shah et al. (2018) for wintertime PM<sub>2.5</sub> over the Eastern US. For eastern China, Fu et al. (2017) and Lachatre et al. (2019) showed both from modelling and satellite observations that such processes lead to strongly increased ammonia tropospheric columns. Finally, several works compared CTM simulations to specific observations. For instance, Pay et al. (2012) showed that the G-ratio was generally underestimated over Europe, inducing that the Caliope model they used could probably overestimate the efficiency of NH<sub>3</sub> emission reductions. Petetin et al. (2016) came to a similar conclusion comparing CHIMERE CTM simulations to observations in the Paris region. They ascribed this underestimation to missing NH<sub>3</sub> emissions especially during warmer periods. The formation of SOA results from even more complex reactions involving photo-chemical oxidation (as for SIA), nitration, fragmentation, and oligomerisation of gaseous precursors or secondary products (Kroll and Seinfeld, 2008; Shrivastava et al., 2017). Models generally use simplified parameterizations to calculate the SOA formation yield from various classes of parent VOCs (Tsigaridis et al., 2014), and comparison with measurements often show that SOA sources are still missing in models (Huang et al., 2020; Tsimpidi, 2016).

In this work, we focus on one specific region in Europe, the Po basin. In a companion paper (Clappier et al., 2021) that analyses PM secondary formation chemical regimes across Europe, the Po basin is clearly identified as a peculiar area where the chemical regime distributions are the most complex, showing non-linear processes (Thunis et al., 2013 and 2015; Carnevale, 2020; Bessagnet, 2014), especially those related to interactions between NO<sub>x</sub> and NH<sub>3</sub>. The Po basin is also one of the pollution hot spots in Europe where the number of days above the limit values prescribed by the European Ambient Air Quality Directives (AAQD) for PM<sub>10</sub> is yet largely exceeded (EEA, 2020). This situation results from the high emission density in this region and also from the geographical setting of the area, in border of the Alps and Apennines mountain ranges that lead to very weak winds in the area, favouring the accumulation of atmospheric pollutants.

We focus the present analysis on the NH<sub>3</sub>-NO<sub>x</sub> chemical processes and describe their spatial and seasonal variability, which could help to design more effective mitigation strategies. We start by describing the modelling set-up and detail the series of simulations required to perform our analysis. Section 3 provides a brief overview of the modelled base-case concentrations. In Section 4, we analyse the sensitivity of PM<sub>2.5</sub> concentrations to NH<sub>3</sub> and NO<sub>x</sub> emissions. Section 5 provides an analysis of the non-linearity in PM<sub>2.5</sub> response to these emissions while in
95

Section 6, we discuss the implications of these results in terms of mitigation measures and design of air quality plans. Conclusions are finally proposed.

## 2. Methodology

### 2.1 Modelling set-up

100

The modelling study is performed with the EMEP/MSC-W air quality model, version rv4_17 (Simpson et al., 2012). The emission input consists of gridded annual national emissions ($SO_2$, NO, $NO_2$, $NH_3$, NMVOC, CO and primary PM$_{2.5}$) at 0.1 × 0.1 degrees resolution, based on data reported every year by parties to the Convention on Long Range Transboundary Air Pollution (CLRTAP). These emissions are provided for 10 anthropogenic source-sectors classified by SNAP (Selected Nomenclature for Air Pollution) codes (EMEP, 2003). Meteorological input data are based on forecasts from the Integrated Forecast System (IFS), a global operational forecasting model from the European Centre for Medium-Range Weather Forecasts (ECMWF). Meteorological fields are retrieved at a 0.1 × 0.1 degree longitude latitude resolution and are interpolated to the 50 × 50 km$^2$ polar-stereographic grid projection (EMEP, 2011).

105

The gas-phase chemistry is based on the evolution of the so-called “EMEP scheme” as described in Simpson et al. (2012) and references therein. The chemical scheme couples the sulfur and nitrogen chemistry to the photochemistry using about 140 reactions between 70 species (Andersson-Sköld and Simpson, 1999; Simpson et al. 2012). In the EMEP Status Report 1/2004 (Fagerli et al., 2004) the reactions are described that cover acidification, eutrophication and ammonium chemistry. The aqueous phase chemistry describes the formation of sulfate in clouds via $SO_2$ oxidation by ozone and $H_2O_2$ and catalysed by metal ions. An important pathway of particulate nitrate formation is through the hydrolysis of $N_2O_5$ on wet aerosol surfaces that converts $NO_x$ into $HNO_3$. More information on the chemical equations is given in Simpson et al. (2012), section 7.

110

115

The EMEP model has two size fractions for aerosols, fine aerosol (PM$_{2.5}$) and coarse aerosol (PM$_{10-2.5}$). The aerosol components presently accounted for are sulfate ($SO_4^{2-}$), nitrate ($NO_3^-$), ammonium ($NH_4^+$), anthropogenic primary PM and sea salt.

For inorganic aerosols, EMEP uses the MARS equilibrium module to calculate the partitioning between gas and fine-mode aerosol phase in the system of $SO_4^{2-}$, $HNO_3$, $NO_3^-$, $NH_3$ and $NH_4^+$ (Binkowski and Shankar, 1995). Aerosol water is calculated to account for particle water within the PM$_{2.5}$ mass, which depends on the mass of soluble PM fraction and on the type of salt mixture in particles. Sea salt (sodium chloride) and dust components are not accounted for by MARS, which might lead to PM underestimations close to coastal sites and where the dust contribution is important. More information on the gas and aerosol partitioning is given in Simpson et al. (2012), section 7.6.

120

125

Regarding secondary organic aerosols (SOA), the EmChem09soa scheme is used, which is a simplified version of the so-called volatility basis set (VBS) approach (Robinson et al., 2007; Donahue et al., 2009). The VBS mechanism is discussed in detail in Bergström et al. (2012). The main differences between the VBS schemes and EmChem09soa is that all primary organic aerosol (POA) emissions are treated as non-volatile in EmChem09soa. This is done to keep the emission totals of both PM and VOC components the same as in the official emission inventories. The semi-volatile biogenic and anthropogenic SOA species are assumed to further oxidise (also known as ageing process) in the atmosphere by OH-reactions. This will lead to a reduction in volatility for the SOA, and thereby increased partitioning to the particle phase. More information on SOA is given in Simpson et al. (2012), section 7.7.

130

The modelling domain covers the entire Po basin (Figure 1) with a 0.1 by 0.1 degree resolution (polar stereographic projection centred at 60 °N) and includes 20 vertical levels. The initial and background concentrations for ozone are based on Logan (1998) climatology, as described in Simpson et al. (2003). For the other species, background/initial conditions are set within the model using functions based on observations (Simpson et al., 2003; Fagerli et al., 2004). The simulations cover the entire meteorological year 2015. We will not discuss the validation of the base-case simulation, as this is available in other publications (e.g. Simpson et al., 2012) and regular status reports by EMEP (https://emep.int/mscw/mscw_publications.html).

135

140
In this work, we simulated a series of 24 scenarios in which NOₓ and NH₃ emissions were reduced independently or simultaneously by 25, 50, 75 and 100% from the base case reference levels. Emission reductions were applied over the entire Po-basin domain for a complete meteorological year (2015).

## 2.2 Spatial and temporal focus

Results are generally presented in terms of maps but three locations within the domain were selected for a more detailed analysis. The locations are Bergamo (Be) in the northern part of the domain, Mantova (Ma) in the central eastern part of the Po basin and Bologna (Bo) in its southern part. As described in the following sections, we will see that these locations show very different behaviours in terms of response to emission changes.

We also aggregate results into 2 seasons: winter and summer that cover the period from November to February and from April to September, respectively. These two seasons are characteristic of different chemical regimes as illustrated in the following sections. The process to define the temporal bounds of these two seasons is discussed in Annex 1. The two remaining months (March and October) represent transition periods and are not considered in our analysis.

As we only analyse processes involving inorganic gas-phase precursors, our focus is on secondary inorganic PM₂.₅ although most of the results are expressed in terms of total PM₂.₅ concentrations. The impact of NOₓ emission reductions on SOA concentration is only briefly discussed in section 5.

## 2.3 Indicators

To describe the interactions between NH₃ and NOₓ emissions, we use the relationship proposed by Stein and Alpert (1993) and Thunis and Clappier (2014). This relation expresses the change of concentration resulting from a reduction of both precursors NOₓ and NH₃ simultaneously, as the sum of two single concentration changes and an interaction term, as follows:

$$ \Delta C_{NO_xNH_3}^{\alpha} = \Delta C_{NO_x}^{\alpha} + \Delta C_{NH_3}^{\alpha} + \hat{C}_{NO_xNH_3}^{\alpha} \quad (1) $$

Where ΔC stands for the PM₂.₅ concentration change (reference minus scenario) for a percentage α emission reduction (thus the term ΔC/α is defined positive for a concentration reduction, consecutive to an emission reduction) and ˆC for the interaction term. We then scale each of these terms by the emission reduction (α) to generate potential impacts (P) (Thunis and Clappier, 2014).

$$ \frac{\Delta C_{NO_xNH_3}^{\alpha}}{\alpha} = P_{NO_xNH_3}^{\alpha} = P_{NO_x}^{\alpha} + P_{NH_3}^{\alpha} + \hat{P}_{NO_xNH_3}^{\alpha} = \frac{\Delta C_{NO_x}^{\alpha}}{\alpha} + \frac{\Delta C_{NH_3}^{\alpha}}{\alpha} + \frac{\hat{C}_{NO_xNH_3}^{\alpha}}{\alpha} \quad (2) $$

This division by the factor α is a mean to extrapolate virtually the impact resulting from any percentage emission reduction to 100%. Potential impacts facilitate the comparison of concentration changes obtained for different emission reduction levels. Indeed, equal potential impacts imply a linear relationship between emission reductions and concentration changes. For example, $P^{\alpha} = P^{\beta} \Rightarrow \Delta C^{\alpha} = \frac{\alpha}{\beta} \Delta C^{\beta}$, for α and β, two emission reduction levels.

The overall potential impact is therefore the sum of two single potential impacts and one interaction term.

A relation between the potential impacts of combined emission reductions at two levels of intensity is obtained by writing (2) for two reduction levels α and β and subtract the two equations. This leads to the following relation:

$$ P_{NO_xNH_3}^{\beta} = P_{NO_xNH_3}^{\alpha} + \frac{P_{NO_x}^{\beta} - P_{NO_x}^{\alpha}}{\hat{P}_{NO_xNH_3}^{\beta-\alpha}} + \frac{P_{NH_3}^{\beta} - P_{NH_3}^{\alpha}}{\hat{P}_{NH_3}^{\beta-\alpha}} + \frac{\hat{P}_{NO_xNH_3}^{\beta} - \hat{P}_{NO_xNH_3}^{\alpha}}{\hat{P}_{NO_xNH_3}^{\beta-\alpha}} \quad (3) $$

or equivalently via equation 2:
$$P_{NO_xNH_3}^{\beta} = P_{NO_x}^{\alpha} + P_{NH_3}^{\alpha} + \frac{\hat{P}_{NO_xNH_3}^{\alpha} + \hat{P}_{NO_x}^{\beta-\alpha} + \hat{P}_{NH_3}^{\beta-\alpha} + \hat{P}_{NO_xNH_3}^{\beta-\alpha}}{\text{non-linear terms}} \quad (4)$$

While the two first terms on the right hand side of Eq. 4 represent the single potential impacts, the remaining right hand side terms quantify the magnitude of non-linearities. $\hat{P}_{NO_xNH_3}^{\alpha}$ quantifies the NOx-NH3 interaction at level $\alpha$, $\hat{P}_{NO_x}^{\beta-\alpha}$ and $\hat{P}_{NH_3}^{\beta-\alpha}$ are the single non-linearities associated to NOx and NH3 emissions, respectively, between levels $\alpha$ and $\beta$, whereas $\hat{P}_{NO_xNH_3}^{\beta-\alpha}$ represents the incremental change of the NOx-NH3 interaction between levels $\alpha$ and $\beta$.

Information about non-linearity is important to design air quality plans as it informs on the robustness of a given response, i.e. whether or not this response remains valid over a certain range and type of emission reductions.

Because air quality models often provide responses for a limited set of scenarios that are then used as a basis to interpolate/extrapolate the responses to other emission reduction levels, robustness shall always be carefully assessed.

In the next section, we present the baseline results in terms of spatial and temporal variations.

### 3. Baseline concentrations of PM<sub>2.5</sub> and gaseous inorganic precursors

Before analysing the impact of emission changes on concentrations, it is worth having a look at the baseline concentration fields. In Figure 1, the yearly averaged PM<sub>2.5</sub> concentration fields show a widespread pollution plume covering most of the area, with peak values extending in its central part. The maximum modelled yearly values reach 29 µg m<sup>-3</sup>, that represent an average between maximum winter values (maximum of 59 µg m<sup>-3</sup>) and minimum summer values (17 µg m<sup>-3</sup>).

The seasonal fields of PM<sub>2.5</sub> clearly show that high yearly average values mostly result from the winter season contributions when more stable atmospheric conditions lead to stagnant conditions, favouring the accumulation of particulate matter in the area (Pernigotti et al., 2014; Raffaelli et al., 2020). The increased emissions from the residential sector (heating, especially wood burning) foster this process (Ricciardelli et al., 2017; Hakimzadeh et al., 2020). Wintertime low temperatures also favour the partitioning of semi-volatile components (e.g. ammonium nitrate) towards the particulate phase. Overall, the relative contribution of secondary inorganic particles (SIA) ranges between 40 and 50%, regardless of the season and is quite homogeneously distributed spatially over the entire area. Strategies targeting SIA have therefore the potential to abate about half of the total PM<sub>2.5</sub> concentration.

As mentioned earlier, the secondary inorganic fraction of PM<sub>2.5</sub> results from complex atmospheric processes that involve gaseous precursors (mainly SO<sub>2</sub>, NO<sub>x</sub> and NH<sub>3</sub>), that can be summarized by the two following chemical pathways:

$$ \begin{aligned} \mathrm{NO_2(g)} + \mathrm{OH} + \mathrm{M} &\rightarrow \mathrm{HNO_3(g)} + \mathrm{M} \\ \mathrm{HNO_3(g)} + \mathrm{NH_3(g)} &\leftrightarrow \mathrm{NH_4^+(p)} + \mathrm{NO_3^-(p)} \end{aligned} \quad (5) $$

$$ \begin{aligned} \mathrm{SO_2(g)} + \mathrm{OH} &\rightarrow \mathrm{H_2SO_4(g)} \\ \mathrm{SO_2(aq)} + \mathrm{O_3}, \mathrm{H_2O_2(aq)} &\rightarrow \mathrm{H_2SO_4(aq)} \\ \mathrm{H_2SO_4(g)} &\rightarrow \mathrm{H_2SO_4(p)} \rightarrow 2\mathrm{H^+(p)} + \mathrm{SO_4^{2-}(p)} \\ \mathrm{H_2SO_4(p)} + 2\mathrm{NH_3(g)} &\rightarrow 2\mathrm{NH_4^+(p)} + \mathrm{SO_4^{2-}(p)} + \mathrm{H_2O} \end{aligned} \quad (6) $$

where (g) means “gas phase”, (aq) aqueous phase, while (p) means “particulate matter”, and the character → symbolizes a chemical pathway that summarizes a set of underlying reactions.

The second of these pathways is generally slower than the first one, the NO<sub>2</sub> oxidation specific time constant being typically some hours to a day, and that of SO<sub>2</sub> one to several days (Seinfeld and Pandis, 2006).
The spatial fields for the seasonal average concentrations of these precursors (Figure 2) reflects their emissions spatial patterns resulting in a NO₂-rich area that comprises Milan plus its northern districts (and to a lesser extent, Turin), while NH₃ is more abundant in the central part of the Po basin, east of Milan, where intensive agriculture practices take place. Finally, high SO₂ concentrations are collocated with the NO₂-rich areas nearby Milan but with an additional zone around the harbour city of Genova (along the south coast), reflecting the more important SO₂ emissions from the shipping sector there. However, SO₂ concentrations are about one order of magnitude below those of NO₂. Seasonal variations are well marked for NO₂ and SO₂ concentrations, not so much in terms of minimum and maximum values but rather in terms of spread with an extended high concentration zone during wintertime. In contrast, NH₃ concentrations remain very similar in summer and winter, both in terms of values and spatial distribution.

In next sections, we simulate a series of emission reduction scenarios to analyse the response of PM₂.₅ concentration to single and combined reductions of NH₃ and NOₓ.

# 4. Analysis of the SIA formation chemical regimes

## 4.1 Seasonal trends

From a strategic point of view, it is important to know whether NOₓ (mostly emitted by the transport, industry and residential sectors) or NH₃ (mostly emitted by agriculture) need to be reduced in priority in order to reach effective results on particulate pollution mitigation. In Clappier et al. (2021), we have observed a great heterogeneity in SIA formation chemical regimes across the Po-basin, different regimes being present in limited geographical areas. Here we intend to look at these regimes in more detail.

To analyse in details the chemical regimes in the Po basin, we compare the two single potential impacts P<sub>NOx</sub><sup>25%</sup>, P<sub>NH3</sub><sup>25%</sup> obtained for moderate emission reductions of 25%. Figure 3 provides a spatial overview of the difference between these two potential impacts (P<sub>NOx</sub><sup>25%</sup> - P<sub>NH3</sub><sup>25%</sup>). This indicator tells whether reductions of NOₓ or NH₃ will lead to the greatest PM₂.₅ concentration abatement, i.e. if the regime is rather NOₓ- or NH₃-sensitive, with positive and negative values, respectively.

* During summer time (Figure 3 – left), the entire area is under weak NOₓ-sensitive conditions with a maximum intensity in its central part, between Bergamo and Mantova.

* During wintertime (Figure 3 – right), the situation is contrasted with a wide and intense NH₃-sensitive area that appears around and south-eastwards of Bergamo. This area includes big cities like Milan. Other (not as marked) NH₃ sensitive regime zones appear nearby coastal areas. Most strongly NOₓ-sensitive areas are located in the eastern parts of the domain, north of Bologna and Venice. NH₃-sensitive regimes are generally collocated with the NO₂- and SO₂- rich areas (Figure 2), whereas NOₓ-sensitive regimes coincide with NH₃-rich areas. The cases of the three selected cities (Bergamo, Mantova and Bologna) representative of the NH₃-sensitive, NOₓ-sensitive and neutral regimes, respectively, are further analysed below.

The chemical regimes deduced from the results of emission reduction scenarios can be compared with the maps of the G-ratio (Figure 4), defined by Ansari and Pandis (1998) as the ratio between free ammonia (NH₃ and NH₄⁺) and total nitrate (HNO₃ + NO₃⁻) after neutralization of H₂SO₄. Values of the G-ratio below 1 indicate a NH₃-limited chemical regime, while values above 1 characterize a HNO₃-limited chemical regime.

$$G = \frac{NH_3(g) + NH_4^+(p) - 2SO_4^{2-}(p)}{HNO_3(g) + NO_3^-(p)} \quad (7)$$

During summer, the G-ratio values well above unity indicate a HNO₃ limited chemical regime across the Po basin. This corresponds to a NOₓ sensitive chemical regime in this region (Figure 3). Moreover, the location of the G-ratio maximum between Bergamo and Mantova spatially coincides with the most pronounced NOₓ sensitive regime, and to a maximum of NH₃ concentrations of about 20 µg m⁻³ (Figure 2). Indeed, NOₓ emission reductions lead to HNO₃ concentration reductions that is the limiting factor in NH₄NO₃ formation according to the G-ratio. During winter, the
G-ratio still shows large values in the region south-east of Bergamo, but contrary to summer, the chemical regime is clearly NH₃ sensitive (Figure 3). More generally, G-ratio values remain above unity over the whole Po basin while both NOₓ and NH₃ sensitive chemical regimes prevail in different areas. Therefore, the G-ratio, related to the abundance of total free ammonia and total nitrate, provides information which differs from that obtained by determining the distribution of the NH₃ and NOₓ emission sensitive chemical regimes. These differences illustrate the impossibility to directly use the G-ratio for air quality management, an interesting result in itself. We will further discuss this interesting behaviour later when addressing non-linearity in section 4.3.

## 4.2 Impact of the emission reduction strength

In this section, we repeat the analysis of Section 4.1 for yearly average concentrations but looking at the step changes of regimes as we progressively reduce emissions from the base case situation. (Figure 5). Chemical regimes are well in place for a 25% level reduction (top left) and are only slightly perturbed from 25 to 50% with a reinforcement of the NOₓ limited regime (top right). Despite this slight change in intensity, the regimes keep therefore the same spatial patterns. From 50% onward, chemical regimes tend to attenuate and reverse themselves from 75% to 100% (bottom right).

In other words, locations that are NH₃-sensitive for the first steps emission reductions will become NOₓ-sensitive for the last steps emission reductions, and vice versa.

## 4.3 A summarized overview: the PM₂.₅ isopleths

Like isopleth plots that show the variations in the O₃ concentrations as a function of NOₓ and VOC concentrations (Dodge, 1977), similar plots can be created for PM₂.₅ concentrations as a function of NOₓ and NH₃ emissions. Simulation results have indeed often been presented as 2D isopleths of PM₂.₅ or nitrate as a function of precursor emissions, which allows showing in a comprehensive manner their sensitivity, and also, in a qualitative manner, non-linear effects (for example Watson et al., 1994 over USA; Xing et al., 2018 over the Beijing-Tianjin-Hebei region in China). Figure 6 shows the PM₂.₅ isopleths obtained through an interpolation among the 25 simulation concentration values (these 25 simulations correspond to the white square symbols in each isopleth diagram) at the three locations Be, Ma, and Bo previously defined. The X- and Y-axes represent the strengths of the NH₃ and NOₓ emission sources, respectively. With this type of graphical representation, it is possible to visualize the response of PM₂.₅ to a NOₓ emission change by moving vertically, the reaction to a NH₃ emission change by moving horizontally or the reaction to a combined NOₓ-NH₃ emission change by moving diagonally. The larger the number of isopleths we cross on the path (high gradient), the larger the expected impact from an emission reduction will be. A simple theoretical model to generate and interpret these isopleths is proposed in Annex 2.

From the analysis of the isopleths, we note the following points:

* In general, the isopleths show a regular pattern with a progressive decrease of PM₂.₅ concentration when either the NOₓ or NH₃ emissions are reduced, with the only exception of Bergamo during wintertime, where NOₓ reductions up to 70% lead to a small increase of PM₂.₅ (Figure 6 - top left), whatever the reduction in NH₃ emission is. We discuss this particular feature later in this section.

* The diagram areas can be divided into two zones separated by a ridge (dashed line in Figure 6). Above the ridge line, PM₂.₅ is more sensitive to NH₃ while below the ridge line, it is more sensitive to NOₓ emission reductions. The orientation of the ridge (tending to vertical or horizontal) informs on the type of chemical regime (NOₓ or NH₃ sensitive, respectively).

* The efficiency of a NOₓ vs NH₃ emission reduction varies across locations. We can compare the efficiency of a given reduction by looking at the horizontal (for NOₓ) and vertical (for NH₃) gradients. To support this comparison, we included in each diagram dashed oblique lines that connect similar PM₂.₅ concentration values for single NOₓ and NH₃ emission reductions. The more vertical are these lines, the larger is the NH₃ abatement impact compared to the NOₓ abatement impact. Conversely, the more horizontal they are, the larger is the NOₓ abatement impact compared to the NH₃ abatement impact. For moderate emission reductions (up to 50%, top right corner), different behaviours are observed: while in Bergamo PM₂.₅ is more sensitive to NH₃ reductions, it is more sensitive to NOₓ reductions in Mantova, and equally sensitive to both precursors in Bologna. This corresponds to the spatial patterns of NOₓ- and NH₃-sensitivive regimes depicted in Fig. 3.
* Winter- and summertime isopleths show completely different patterns in Bergamo whereas at the two other locations, they remain similar.

* At Mantova where moderate NOx reductions (e.g. 50%) are the most efficient among the 3 sites, NH3 emission reductions are more efficient than NOx emission reductions for larger additional reductions (going for example from 75% to 100%). At Bergamo, NH3 reductions are the most effective for moderate reductions whereas NOx reductions become more effective for larger reductions, as seen by the isopleths spacing. This confirms the finding of reversed chemical regimes for larger additional emission reductions detailed in the previous section and illustrated in Figure 5.

The special pattern of Bergamo's PM2.5 isopleths during wintertime needs some additional discussion. The increase of the inorganic fraction of PM2.5 as a response to NOx reductions during wintertime has already been noted by several authors (e.g. Le et al., 2020; Sheng et al., 2018). It has been related to an increase in the oxidizing capacity of the atmosphere and in particular to increased ozone levels. This is due to the prevailing titration of O3 by NO in wintertime high NOx conditions and in the absence of photochemical ozone production due to reduced solar radiation (Kleinman et al., 1991).

$$NO(g) + O_3(g) \rightarrow NO_2(g) + O_2(g) \quad (8)$$

The impact of NOx emission reductions on the concentration of various pollutants in Bergamo during wintertime is illustrated in Figure 7. As expected, a 50% reduction in NOx emissions leads to a decrease in NO2 concentration (from 47 to 28 µg m⁻³, i.e. a factor of 1.7). In contrast, O3 concentration increases from 8 to 16 µg m⁻³, roughly a factor of two. These compensating changes result in a small increase in NO3 radical production (Eq. 9), the initial step of the major pathway of wintertime HNO3 and particulate nitrate formation (Kenagy et al., 2018).

$$NO_2(g) + O_3(g) \rightarrow NO_3(g) + O_2(g) \quad (9)$$

In this pathway, the NO3 radical formation is followed by combination with NO2 to form N2O5, a reversible process, and heterogeneous HNO3 formation on wet particle surfaces.

$$NO_3(g) + NO_2(g) \leftrightarrow N_2O_5(g) \quad (10)$$

$$N_2O_5(g) + H_2O(aerosol surface) \rightarrow 2HNO_3(g) \quad (11)$$

The NO3 radical has three major rapid sinks, reaction with NO, photolysis, and reaction with NMVOCs especially terpenes.

$$NO_3(g) + NO(g) \rightarrow 2NO_2(g) \quad (12)$$

$$NO_3(g) + hv \rightarrow NO_2(g) + O(g) \quad (13)$$

$$NO_3(g) + hv \rightarrow NO_2(g) + O_2(g) \quad (14)$$

$$NO_3(g) + NMVOC(g) \rightarrow \text{products} \quad (15)$$

Reactions 8 to 15 induce additional dependence of HNO3 formation on NOx species on top of (5), but which partly cancel out, as they are both involved in formation and sink processes.

SOA is formed through a series of chemical reactions of gaseous precursors (mainly volatile, intermediate volatile or semi-volatile organic compounds (VOCs) with the oxidants O3, OH and nitrate radical (NO3) (Li et al., 2011).
$$ (S, I)\text{VOC(g)} + \text{oxidants} \rightarrow (S, I)\text{OVOC(g)} \leftrightarrow \text{SOA(s)} \quad (16) $$

Putting all the arguments together, it follows that wintertime ammonium nitrate formation over Bergamo is most probably controlled by $NO_3$ radical formation (9). The fact that this behaviour is observed in Bergamo and not in Mantova or Bologna is due to the much larger $NO_2$ levels in the Bergamo - Milano area (above 50 µg m⁻³ during winter, Fig.2). Such large $NO_2$ levels are also simulated locally over the Turin area, and also lead to a slightly $NH_3$ sensitive regime there despite a G-ratio well above unity. Beyond 50% $NO_x$ reduction, $NH_4NO_3$ formation decreases because $NO_2$ decreases more rapidly than ozone increases up to its maximum (at 75% $NO_x$ emission reduction, see Figure 7).

The negative response of $NH_4NO_3$ to $NO_x$ emission reductions in Be during wintertime explains the apparent discrepancies with the analysis of G-ratio, which indicates that $NH_4NO_3$ is strongly $HNO_3$ limited. Simply, the $HNO_3$-limited chemical regime cannot be extrapolated to sensitivity to $NO_x$ emissions in case of the above shown negative response. Total nitrate is less abundant than free ammonia (defined as $NH_3(g) + NH_4^+(p) - 2SO_4^{2-}(p)$), but $NO_x$ emission reductions up to about 50% do not reduce its concentration, and $NH_3$ emission reductions are thus more efficient. In this respect, the G-ratio cannot provide information about negative responses.

At Bergamo during winter, the increase in $PM_{2.5}$ (+ 1.8 µg m⁻³) arising from a 50% reduction in $NO_x$ emission also results from an increase in sulfate (+ 0.3 µg m⁻³) and in SOA (+ 0.6 µg m⁻³) concentrations. Both sulfate and SOA concentrations are closely related to $O_3$ concentrations (Figure 7). The sulfate increase is comparable in magnitude to the nitrate increase, even if sulfate levels are much smaller than nitrate ones. Figure 7 actually shows a strikingly similar response of sulfate, SOA and ozone to $NO_x$ emission reductions (given that $SO_2$ and NMVOC emissions are held constant). Indeed the prevailing wintertime aqueous production of $H_2SO_4$ requires oxidants and in particular ozone (Le et al., 2020; Sheng et al., 2018). In addition, the formation of SOA in both the gas and particulate phases also requires oxidants (Vahedpour et al., 2011; Huang et al., 2020; Feng et al., 2016; Li et al., 2011; Tsimpidi et al., 2010).

Pinder et al. (2008) also note an oxidant limitation for SIA formation over Eastern United states for the 2000 to 2020 period, but in their simulations, it mainly affects sulfate that increases as a result of $NO_x$ emission reductions while nitrate decreases. This is due to a more important sulfate to nitrate ratio in eastern US than over the Po basin. Fu et al. (2020) derive from combined measurements and modelling that wintertime nitrate during haze events in the North China Plain (NCP) are nearly insensitive to 30% $NO_x$ emission reductions, because increased ozone levels increase the $NO_x$ to $HNO_3$ conversion efficiency. Following these authors, this conversion also involves the homogeneous $HNO_3$ formation via the $NO_2 + OH$. This reaction also could play a role in the Po basin, in addition to the heterogeneous pathway. Also Leung et al. (2020) simulate that wintertime nitrate abatement in the NCP is buffered with respect to emission reductions by increased oxidant build-up, but also by sulfate to nitrate conversion by liberating free $NH_3$ through sulfate concentration reduction, which can then enhance nitrate formation. Womack et al. (2019) find an oxidant limitation of nitrate formation over wintertime Utah (USA), and show that nitrate concentration diminishes when reducing VOC emissions.

## 5. Analysis of non-linearities

Clappier et al. (2021) highlighted the specificities of the Po basin area within Europe. They showed that non-linearities are present in this region. One peculiarity of the Po basin is a marked difference between the chemical regimes encountered within a confined area, which have implications on the linearity of $PM_{2.5}$ responses to emission changes. In this section, we analyse in more details these non-linearities.

Information on the $NO_x-NH_3$ interaction term at the reduction level $\alpha=25\%$: $P_{NO_xNH_3}^{25\%} = [P_{NO_xNH_3}^{25\%} - (P_{NO_x}^{25\%} + P_{NH_3}^{25\%})]$ (first non-linear term in Eq. 4) is provided in Figure 8. At $\alpha=25\%$, the interaction term is negative (or null) across the entire modelling domain (most data points are below the 1:1 line) regardless of the chemical regimes, and averages to approximately -10%, as indicated by the linear fit slope (= 0.9). In relative terms, this interaction term is also almost constant regardless of the season (not shown).
This negativity can be explained by the fact that a reduction of only NOₓ implies a reduction of both NO₃⁻ and NH₄⁺ and the same happens when reducing only NH₃, therefore a simultaneous reduction of both precursors is lower than the sum of the two. Single impacts would therefore lead to an overestimation (of about 10%) in PM₂.₅ reduction if added up to extrapolate linearly the impact of combined 25% NOₓ and NH₃ emission reductions on yearly averaged PM₂.₅ concentrations. This result is expected for what concerns particulate NH₄NO₃, as a consequence of the gas/particle equilibrium described in Eq. 5, although non-linear relationships between NOₓ emissions and HNO₃ concentrations also play a role. Qualitatively, this negative interaction is also highlighted by the hyperbolic shapes of the PM₂.₅ isopleths determined for 3 different sites of the domain (Figure 6). As discussed in Annex 2, linearity would result in isopleths parallel to the descending diagonal lines.

When emission reductions increase from 25% to 50%, three additional non-linear terms are generated (three last terms in equation 4). Figure 9 and 10 provide an overview of these non-linear terms during wintertime and summertime, respectively. The top left panel of each figure represents their sum, i.e. the total non-linearity generated between 25 and 50% emission reduction. The right column shows the non-linearities associated to NOₓ and NH₃ while the bottom left panel reports the non-linear interaction between the two precursors.

Overall, non-linearities are more important during wintertime than during summertime. This is true both in absolute and relative (not shown) terms. Non-linearities tend to be the largest in between areas that are characterised by well-marked NH₃- or NOₓ-sensitive regimes (indicated by the blue and red drawn contours). This can be explained by the fact that when one of the two components (NOₓ or NH₃) is in large excess (compared to the other one), reductions of this compound have then generally only little impact, implying that both single and combined reductions only involve one compound and are therefore similar.

During wintertime, the overall non-linearity (Figure 9 top left) is largely dominated by the single NOₓ-related non-linearity ($\hat{P}_{NO_x}^{50-25\%}$, Figure 9 top right), a singularity in Europe as the Po basin is the only area where this occurs to this extent (Clappier et al., 2021). In the region of Bergamo, the NOₓ non-linearity remains weak despite the peculiar PM₂.₅ responses to NOₓ emission reductions (i.e. an increase of PM₂.₅ concentrations for NOₓ emission reduction up to 50%, Figure 7). In this NH₃ sensitive region, this behaviour can be explained by the strong oxidant limitation of HNO₃ formation outlined above (Figure 7). It is worth mentioning that although atypical, this behaviour is quasi-linear with PM₂.₅ responses that remain proportional to the emission reduction strength (up to 50%) but with a negative slope. Similarly to NOₓ, NH₃ single non-linearities (Figure 9, bottom right) are positive but weaker, indicating slightly larger potential impacts for emission reductions in the range 25-50% than in the 0-25% range. Finally, the NOₓ-NH₃ non-linear interaction terms (Figure 9, bottom left) are mostly negative, indicating that the interaction term for 50% emission reductions is more negative than the corresponding term for 25%, pointing out to a strengthening of the NOₓ-NH₃ non-linearity when more intense emission reductions are considered.

In relative terms, the overall wintertime non-linearity terms increase by about 30% when emission reductions increase from 25 and 50% (Figure 11 top, blue points). Note that these non-linearity terms reach larger values in some places of the modelling domain as highlighted by the data point dispersion. In a previous work, Thunis et al. (2015) quantified the non-linearity of model responses to emission reductions in three areas in Europe, among which the Po Basin. One of their conclusions was that non-linearities remain relatively low for yearly averaged responses. Although the results presented here show important nonlinearities, these occur mainly during wintertime and are limited to specific areas. It is also worth noting that these non-linear responses (for moderate emission reductions up to 50%) only occur in the Po Basin (Clappier et al., 2021).

During summertime, the magnitude of non-linear terms is smaller than during wintertime. The overall non-linearity term (top left) is dominated by the NH₃-NOₓ non-linear interaction term ($\hat{P}_{NO_xNH_3}^{50-25\%}$, Figure 10, bottom left). In summer, PM₂.₅ concentration is NOₓ sensitive in almost all the domain (Figure 3) and emission reductions do not lead to shifts in the chemical regimes. In this case, the interaction terms become more important. The non-linear interaction terms are negative everywhere in the Po basin, implying again that the interaction terms at 50% is more negative than at 25% emission reduction (reinforcement of the non-linearities when more intense emission reductions are considered). Most negative values appear in the western and northern part of domain.
Figure 11 shows the increase of the overall non-linear terms for emission reduction steps starting from different points, from 25 to 50%, from 50 to 75%, and finally from 75 to 100%. Regardless of the emission reduction step, summer time non-linearities remain small all over the domain with regression slope close to 0.9 and very limited data point dispersions (i.e. low RMSE). Wintertime non-linearities further increase significantly from 50% to 75% reduction levels (regression fit parameter close to 1.21) but tend to stabilize between 75 and 100% reduction level (regression fit parameter close to 1.05). It is interesting to note that potential impacts in winter increase for all segments (all winter points are above the 1:1 line) indicating that the same percentage reduction (25%) gains progressively more impact when more intense reductions are considered.

## 6. Discussion

When designing air quality plans, it is important to identify the key precursors on which to act in priority to hit a specific air quality target, but also to understand the consequences of these choices for various seasons (temporal variations), locations (spatial variability), emission reduction levels (strength) and strategies (combined or single emission reductions). As the information to make this decision is generally incomplete, assessing the robustness of the available model responses is essential. From the results presented here, a few key points appear.

* The seasonal and spatial variabilities in the response of PM<sub>2.5</sub> to the reduction of NO<sub>x</sub> and NH<sub>3</sub> emissions are extremely large, with different and sometimes opposite responses to emission changes. Yearly averages do not represent the appropriate time window to evaluate the impact of such emission reductions and a focus on wintertime (November to February) seems to be the right option, especially because concentrations are larger during this period of the year.

* The responses of PM<sub>2.5</sub> to emission reduction plans that cover the whole area (i.e. uniform emission reductions are applied everywhere in the domain) vary from location to location: opposite responses occur within a few hundreds of kilometres for some reduction levels. In the region of Bergamo, PM<sub>2.5</sub> response to NO<sub>x</sub> emission reductions can be negative, meaning an increase of PM<sub>2.5</sub> when reducing the NO<sub>x</sub> emissions. It is important to combine NO<sub>x</sub> and NH<sub>3</sub> emission reductions in winter, or to go for stronger emission reductions to make sure these unwanted effects are limited.

* Despite quite important non-linearities, PM<sub>2.5</sub> responses to emission reductions are not chaotic. Indeed, regardless of the emission reduction level, the non-linear terms related to NH<sub>3</sub> emission reduction and to NO<sub>x</sub>-NH<sub>3</sub> interactions are quite uniform spatially. This is not the case of NO<sub>x</sub> emission reduction, for which care must be taken to ensure that the detailed response of PM<sub>2.5</sub> is captured.

* Although they are location specific, PM<sub>2.5</sub> isopleths represent an interesting tool to assess the impact of different NO<sub>x</sub> and NH<sub>3</sub> emission reductions on PM<sub>2.5</sub> concentration. They indeed allow visualising in one single diagram the impact of any type of reductions on concentrations in a single grid cell (or set of grid cells). We must however remember that these isopleths derive from uniform emission reduction over the whole domain. Comparing sites where PM<sub>2.5</sub> responses to the same "domain-wide" policy are different, it appears challenging to define a single domain-wide policy efficiently reducing PM at all locations.

## 7. Conclusions

In this work, we analysed the sensitivity of PM<sub>2.5</sub> to NO<sub>x</sub> and NH<sub>3</sub> emissions by means of a set of EMEP simulations performed with different levels of emission reductions, from 25% up to a total switch-off of those emissions. Both single and combined precursor reduction scenarios were applied to determine the most efficient emission reduction strategies and quantify the interactions between NO<sub>x</sub> and NH<sub>3</sub> emission reductions. Our results confirm the peculiarity of secondary inorganic PM<sub>2.5</sub> formation in the Po basin suggested by Clappier et al. (2021), characterised by contrasting chemical regimes within distances of few (hundreds of) kilometres, as well as non-linear responses to emission reductions during wintertime. One of the striking results is the increase of the PM<sub>2.5</sub> concentration levels when NO<sub>x</sub> emission reductions are applied in NO<sub>x</sub>-rich areas, such as the surroundings of Bergamo. The isopleths in the graphs showing PM<sub>2.5</sub>, nitrate, sulfate, SOA and O<sub>3</sub> concentrations as a function of NH<sub>3</sub> and NO<sub>x</sub> emissions (Figure 7) indicate that the increased oxidative capacity of the atmosphere is the cause of the increase of PM<sub>2.5</sub> induced by a reduction in NO<sub>x</sub> emission of up to 50%. This process can have contributed to the absence of
significant PM<sub>2.5</sub> concentration decrease during the COVID-19 lockdowns in many European cities (EEA, 2020; Putaud et al., 2020). It is important to account for this process when designing air quality plans, since it could well lead to transitionary increases in PM<sub>2.5</sub> at some locations in winter as NO<sub>x</sub> emission reduction measures are gradually implemented. At this type of location, mitigation measures that would be optimal in the long-term might not be efficient in the short-term.

Joint analyses of PM sensitivity to emissions and the G-ratio can give a clue if a NH₃ sensitive chemical regime is due to either a lack of NH₃ or to a non-linear and negative response of HNO₃ concentration to NOₓ emissions. In this latter case, the chemical regime is NH₃ sensitive in terms of NH₃ and NOₓ emission reductions, but it can be HNO₃ limited in terms of G-ratio, as observed for the Bergamo – Milano region. Inversely, a G-ratio greater than 1, indicating HNO₃ limitation to particulate nitrate formation, does not necessarily indicate a larger sensitivity to NOₓ than to NH₃ emissions. Thus, our results show the impossibility to directly use the G-ratio for air quality management, an interesting result in itself. While PM<sub>2.5</sub> chemical regimes (determined by the relative importance of the NOₓ vs. NH₃ responses to emission reductions) show large variations seasonally and spatially, they are not very sensitive to moderate (up to 50%) emission reductions. Beyond 25% emission reduction strength, responses of PM<sub>2.5</sub> concentrations to NOₓ emission reductions become non-linear in certain areas of the Po basin mainly during wintertime.

Since sulfate concentrations are low in the Po basin, the impact of SO₂ emission reductions was not evaluated here. However, the simulations performed in Clappier et al. (2021) indicate that air quality improvement plans addressing SO₂ emissions may still lead to additional PM<sub>2.5</sub> decreases. Further works should also test if NMVOC emission would further affect the concentration of oxidants, and subsequently of nitrate (and sulfate) during winter. This depends on the fraction of ozone formed photochemically in the Po basin, compared to the one transported from outside by advection or entrainment.

Finally, it would be important to compare the results obtained in this work from the EMEPrv4_17 model with similar results obtained from other models. With its complex setting, the Po basin represents a good test case for such inter-comparisons.

## References

* Andersson-Sköld, Y., and Simpson, D.: Comparison of the chemical schemes of the EMEP MSC-W and the IVL photochemical trajectory models, Atmos. Environ., 33, 1111–1129, 1999.

* Ansari, A.S., and Pandis, S., 1998: Response of inorganic PM to precursor concentrations. Environ. Sci. Technol. 32, 2706–2714.

* Beekmann, M., Prévôt, A. S. H., Drewnick, F., Sciare, J., Pandis, S. N., Denier van der Gon, H. A. C., Crippa, M., Freutel, F., Poulain, L., Ghersi, V., Rodriguez, E., Beirle, S., Zotter, P., von der Weiden-Reinmüller, S.-L., Bressi, M., Fountoukis, C., Petetin, H., Szidat, S., Schneider, J., Rosso, A., El Haddad, I., Megaritis, A., Zhang, Q. J., Michoud, V., Slowik, J. G., Moukhtar, S., Kolmonen, P., Stohl, A., Eckhardt, S., Borbon, A., Gros, V., Marchand, N., Jaffrezo, J. L., Schwarzenboeck, A., Colomb, A., Wiedensohler, A., Borrmann, S., Lawrence, M., Baklanov, A., and Baltensperger, U.: In situ, satellite measurement and model evidence on the dominant regional contribution to fine particulate matter levels in the Paris megacity, Atmos. Chem. Phys., 15, 9577–9591, https://doi.org/10.5194/acp-15-9577-2015, 2015.

* Bergström, R., Denier van der Gon, H. A. C., Prévôt, A. S. H., Yttri, K. E., and Simpson, D.: Modelling of organic aerosols over Europe (2002–2007) using a volatility basis set (VBS) framework: application of different assumptions regarding the formation of secondary organic aerosol, Atmos. Chem. Phys., 12, 8499–8527, https://doi.org/10.5194/acp-12-8499-2012, 2012.

* Bessagnet, B., Beauchamp, M., Guerreiro, C., et al: Can further mitigation of ammonia emissions reduce exceedances of particulate matter air quality standards? Environmental Science & Policy, 44, 149–163. DOI:10.1016/j.envsci.2014.07.011, 2014.
Blanchard, C.L., Hidy, G.M.: Effects of Changes in Sulfate, Ammonia, and Nitric Acid on Particulate Nitrate Concentrations in the South-eastern United States; J. Air & Waste Manage. Assoc., 53, 283-290, 2003.

Carnevale, C., Pisoni, E., and Volta, M.: A non-linear analysis to detect the origin of PM10 concentrations in Northern Italy, Science of the Total Environment 409 (2010) 182-191, 2020.

Clappier, A., Thunis, P., Putaud, J. P., Beekmann, M., and De Meij, A.: Analysis of the chemical regimes and non-linearity associated to PM2.5 in Europe. Submitted in Environment International, 2021.

De Meij, A., Krol, M., Dentener, F., Vignati, E., and Cuvelier, C.: The sensitivity of aerosol in Europe to two different emission inventories and temporal distribution of emissions, Atmos. Chem. Phys., 6, 4287-4309, 2006.

De Meij, A., Thunis, P., Bessagnet, B., and Cuvelier, C.: The sensitivity of the CHIMERE model to emissions reduction scenarios on air quality in Northern Italy., Atmos. Env. Volume 43, Issue 11, Pages 1897-1907, 2009.

Dodge, M. C.: Effect of selected parameters on predictions of a photochemical model, US Environmental Protection Agency, EPA-600/3-77/048, Research Triangle Park, NC, USA, 1977.

EDGAR2020, Emission Database for Global Atmospheric research, https://edgar.jrc.ec.europa.eu/

EEA2020, Air quality in Europe — 2020 report, EEA Report No 11/2020

EMEP 2003, Transboundary acidification and eutrophication and ground level ozone in Europe: Unified EMEP model description" EMEP/MSC-W Report 1/2003.
https://emep.int/publ/reports/2003/emep_report_1_part1_2003.pdf

EMEP 2011, EMEP/MSC-W Model - User's Guide July 2011.

https://wiki.met.no/_media/emep/page1/userguide_062011.pdf

Fagerli, H., Simpson, D., and Tsyro, S.: Unified EMEP model: Updates. In EMEP Status Report 1/2004, Transboundary acidification, eutrophication and ground level ozone in Europe. Status Report 1/2004, pages 11–18. The Norwegian Meteorological Institute, Oslo, Norway, 2004.

Feng, T., Li, G., Cao, J., Bei, N., Shen, Z., Zhou, W., Liu, S., Zhang, T., Wang, Y., Huang, R.-J., Tie, X. and Molina, L. T.: Simulations of organic aerosol concentrations during springtime in the Guanzhong Basin, China, Atmos. Chem. Phys., 16(15), 10045–10061, doi:10.5194/acp-16-10045-2016, 2016.

Feng, T., Zhao, S., Bei, N., Wu, J., Liu, S., Li, X., Liu, L., Qian, Y., Yang, Q., Wang, Y., Zhou, W., Cao, J., and Li, G.: Secondary organic aerosol enhanced by increasing atmospheric oxidizing capacity in Beijing-Tianjin-Hebei (BTH), China, Atmos. Chem. Phys., 19, 7429–7443, https://doi.org/10.5194/acp-19-7429-2019, 2019.

Fu, X., Wang, S., Xing, J., Zhang, X., Wang, T., and Hao, J.: Increasing Ammonia Concentrations Reduce the Effectiveness of Particle Pollution Control Achieved via SO2 and NOX Emissions Reduction in East China, Environ. Sci. Tech. Let., 4, 221–227, https://doi.org/10.1021/acs.estlett.7b00143, 2017.

Fu X, Wang T, Gao J, Wang P, Liu Y, Wang S, Zhao B, Xue L. Persistent Heavy Winter Nitrate Pollution Driven by Increased Photochemical Oxidants in Northern China. Environ Sci Technol. 2020 Apr 7;54(7):3881-3889. doi: 10.1021/acs.est.9b07248. Epub 2020 Mar 16. PMID: 32126767.

Guo, H., Otjes, R., Schlag, P., Kiendler-Scharr, A., Nenes, A., and Weber, R. J.: Effectiveness of ammonia reduction on control of fine particle nitrate, Atmos. Chem. Phys., 18, 12241-12256, 10.5194/acp-18-12241-2018, 2018.

Hakimzadeh, M., E., Soleimanian, A., Mousavi, A., Borgini, C., De Marco, A., Ruprecht, C., Sioutas. The impact of biomass burning on the oxidative potential of PM2.5 in the metropolitan area of Milan, Atmospheric Environment 224 (2020) 117328, 2020.

Huang, L., Q., Wanga, Y., Wang, C., Emeryc, A., Zhua, Y., Zhua, S., Yina, G., Yarwood, K., Zhang, L., Li: Simulation of secondary organic aerosol over the Yangtze River Delta region: The impacts from the emissions of
intermediate volatility organic compounds and the SOA modeling framework, Atmospheric Environment, https://doi.org/10.1016/j.atmosenv.2020.118079, 2020.

Kenagy, H. S., Sparks, T. L., Ebben, C. J., Wooldrige, P. J., Lopez-Hilfiker, F. D., Lee, B. H., et al.: NOx lifetime and NOy partitioning during WINTER. Journal of Geophysical Research: Atmospheres, 123, 9813–9827. https://doi.org/10.1029/2018JD028736, 2018.

Kleinman, L.I.: Seasonal dependence of boundary layer peroxide concentration: The low and high NOx regimes, J. Geophys. Res, vol.96, pp.20721-20733, 1991.

Köble, R., and Seufert, G.: Novel maps for forest tree species in Europe. In: A changing atmosphere, 8th European Symposium on the Physicochemical Behaviour of Atmospheric Pollutants, 17–20 September, Torino, 2001.

Kroll, J.H., and Seinfeld, J. H.: Chemistry of secondary organic aerosol: Formation and evolution of low-volatility organics in the atmosphere, Atmospheric Environment 42 (2008) 3593–3624, 2008.

Lachatre, M., Fortems-Cheiney, A., Foret, G., Siour, G., Dufour, G., Clarisse, L., Clerbaux, C., Coheur, P.-F., Van Damme, M., and Beekmann, M.: The unintended consequence of SO2 and NO2 regulations over China: increase of ammonia levels and impact on PM2.5 concentrations, Atmos. Chem. Phys., 19, 6701–6716, https://doi.org/10.5194/acp-19-6701-2019, 2019

Larsen, B., Gilardoni, S., Stenström, K., Niedzialek, J., Jimenez, J., and Belis, C.: Sources for PM air pollution in the Po Plain, Italy: II. Probabilistic uncertainty characterization and sensitivity analysis of secondary and primary sources, Atmos. Environ., 50, 203–213, https://doi.org/10.1016/j.atmosenv.2011.12.038, 2012.

Le, T., Wang, Y., Liu, L., Yang, J., Yung, Y., Li, G., and Seinfeld, J.: Unexpected air pollution with marked emission reductions during the COVID-19 outbreak in China, Science, 369, eabb7431, 10.1126/science.abb7431, 2020.

Leung, D., Shi, H., Zhao, B., Wang, J., Ding, E.M., Gu, Y., Zheng, H., Chen, G., Liou, K. N., Wang, S., Fast, J. D., Zheng, G., Jiang, J., Li, X., and Jiang, J. H.: Wintertime Particulate Matter Decrease Buffered by Unfavorable Chemical Processes Despite Emissions Reductions in China, GRL, Volume47, Issue14, https://doi.org/10.1029/2020GL087721, 2020.

Li, G., Zavala, M., Lei, W., Tsimpidi, A. P., Karydis, V. A., Pandis, S. N., Canagaratna, M. R., and Molina, L. T.: Simulations of organic aerosol concentrations in Mexico City using the WRF-CHEM model during the MCMA-2006/MILAGRO campaign, Atmos. Chem. Phys., 11, 3789-3809, https://doi.org/10.5194/acp-11-3789-2011, 2011.

Makar, P. A., Moran, M. D., Zheng, Q., Cousineau, S., Sassi, M., Duhamel, A., Besner, M., Davignon, D., Crevier, L.-P., and Bouchet, V. S.: Modelling the impacts of ammonia emissions reductions on North American air quality, Atmos. Chem. Phys., 9, 7183–7212, https://doi.org/10.5194/acp-9-7183-2009, 2009.

Nenes, A., Pandis, S. N., Weber, R. J., and Russell, A.: Aerosol pH and liquid water content de-termine when particulate matter is sensitive to ammonia and nitrate availability, Atmos. Chem. Phys., 20, 3249-3258, 10.5194/acp-20-3249-2020, 2020.

Pay, M.T., Jimenez-Guerrero, P., Baldasano, J.M.: Assessing sensitivity regimes of secondary inorganic aerosol formation in Europe with the CALIOPE-EU modeling system. Atmos. Environ. 51, 146–164, 2012.

Pernigotti, D., Georgieva, E., Thunis, P., Bessagnet, B: Impact of meteorological modelling on air quality: Summer and winter episodes in the Po valley (Northern Italy). Int. J. Envir. Poll., 50, 111–119, 2014.

Petetin, H., Sciare, J., Bressi, M., Gros, V., Rosso, A., Sanchez, O., Sarda-Estève, R., Petit, J.-E., and Beekmann, M.: Assessing the ammonium nitrate formation regime in the Paris megacity and its representation in the CHIMERE model, Atmos. Chem. Phys., 16, 10419–10440, https://doi.org/10.5194/acp-16-10419-2016, 2016.
610

Pinder, R. W., Gilliland, A. B., and Dennis, R. L.: Environmental impact of atmospheric NH3 emissions under present and future conditions in the eastern United States, Geophysical Research Letters, vol. 35, L12808, 2008

Pozzer, A., Tsimpidi, A. P., Karydis, V. A., de Meij, A., and Lelieveld, J.: Impact of agricultural emission reductions on fine-particulate matter and public health, Atmospheric Chemistry and Physics, 17, 12813-12826, 10.5194/acp-17-12813-2017, 2017.

615

Putaud, J.P., Pozzoli, L., Pisoni, E., Martins Dos Santos, S., Lagler, F., Lanzani, G., Dal Santo, U., Colette, A., Impacts of the COVID-19 lockdown on air pollution at regional and urban background sites in northern Italy, ACPD, https://doi.org/10.5194/acp-2020-755,2020.

Raffaelli, K.; Deserti, M.; Stortini, M.; Amorati, R.; Vasconi, M.; Giovannini, G. Improving Air Quality in the Po Valley, Italy: Some Results by the LIFE-IP-PREPAIR Project. Atmosphere, 11, 429, 2020.

620

Ricciardelli, I., D., Bacco, M., Rinaldi, G., Bonafe, F., Scotto, A., Trentini, G., Bertacci, P., Ugolini, C., Zigola, F., Rovere, C., Maccone, C., Pironi, V., Poluzzi: A three-year investigation of daily PM2.5 main chemical components in four sites: the routine measurement program of the Supersito Project (Po Valley, Italy), Atmospheric Environment 152 (2017) 418e430, 2017.

Seinfeld, J. H., and Pandis, S.N: Atmospheric Chemistry and Physics: From Air Pollution to Global Change, 2nd edn., J. Wiley and Sons, New York, 2006.

625

Shah, V., Jaeglé, L., Thornton, J. A., et al: Chemical feedbacks weaken the wintertime response of particulate sulfate and nitrate to emissions reductions over the eastern United States. Proceedings of the National Academy of Sciences of the United States of America. Aug;115(32):8110-8115. DOI: 10.1073/pnas.1803295115, 2018.

Sheng, F., Jingjing, L., Yu., C., Fu-Ming, T., Xuemei, D., Jing-yao, L.: Theoretical study of the oxidation reactions of sulfurous acid/sulfite with ozone to produce sulfuric acid/sulfate with atmospheric implications. RSC Advances. 8. 7988-7996. 10.1039/C8RA00411K, 2018.

630

Shrivastava, M., et al.: Recent advances in understanding secondary organic aerosol: Implications for global climate forcing, Rev. Geophys., 55, 509–559, doi:10.1002/2016RG000540., 2017.

635

Simpson, D., Fagerli, H., Jonson, J. E., Tsyro, S., Wind, P., and Tuovinen, J. P.: The EMEP Unified Eulerian Model. Model Description. EMEP MSC-W Report 1/2003. The Norwegian Meteorological Institute, Oslo, Norway, 2003.

640

Simpson, D., Benedictow, A., Berge, H., Bergström, R., Emberson, L. D., Fagerli, H., Flechard, C. R., Hayman, G. D., Gauss, M., Jonson, J. E., Jenkin, M. E., Nyíri, A., Richter, C., Semeena, V. S., Tsyro, S., Tuovinen, J.-P., Valdebenito, Á., and Wind, P.: The EMEP MSC-W chemical transport model – technical description, Atmos. Chem. Phys., 12, 7825–7865, https://doi.org/10.5194/acp-12-7825-2012, 2012.

Stein, U., and Alpert, P.: Factor separation in numerical simulations, J. Atmos. Sci., 50, 2107-2115, 1993.

645

Thunis, P., and Clappier, A.: Indicators to support the dynamic evaluation of air quality models, Atmospheric Environment, 98, 402-409, 2014.

Thunis, P., Clappier, A., Pisoni, E., and Degraeuwe, B.: Quantification of non-linearities as a function of time averaging in regional air quality modeling applications, Atmospheric Environment, 103, 263-275, 2015.

650

Thunis, P., Pernigotti, D., Cuvelier, C., Georgieva, E., Gsella, A., De Meij, A., Pirovano, G., Balzarini, A., Riva, G. M., Carnevale, C., Pisoni, E., Volta, M., Bessagnet, B., Kerschbaumer, A.,Viaene, P., De Ridder, K., Nyiri, A., and Wind, P.: POMI: A Model Intercomparison exercise over the Po valley, Air Quality, Atmosphere and health, DOI: 10.1007/s11869-013-0211-1, 2013.
Tsimpidi, A. P., Karydis, V. A., Zavala, M., Lei, W., Molina, L., Ulbrich, I. M., Jimenez, J. L., and Pandis, S. N.: Evaluation of the volatility basis-set approach for the simulation of organic aerosol formation in the Mexico City metropolitan area, Atmos. Chem. Phys., 10, 525–546, https://doi.org/10.5194/acp-10-525-2010, 2010.

Tsimpidi, A. P., Karydis, V. A., Pandis, S. N., and Lelieveld, J.: Global combustion sources of organic aerosols: model comparison with 84 AMS factor-analysis data sets, Atmos. Chem. Phys., 16, 8939–8962, 10.5194/acp-16-8939-2016, 2016.

Tsigaridis, K., Daskalakis, N., Kanakidou, M., Adams, P.J., Artaxo, P., Bahadur, R., et al. The AeroCom evaluation and intercomparison of organic aerosol in global models. Atmos Chem Phys. 14(19), 10845–95. doi.org/10.5194/acp-14-10845-2014, 2014.

Vahedpour, M., M., Goodarzi, N., Hajari, and F., Nazari: Struct. Chem., 22, 817–822, 2011.

Venkatram, A., and Pleim, J., The electrical analogy does not apply to modelling dry deposition of particles, Atmos. Environ., 33, 3075–3076. doi: 10.1016/S1352-2310(99)00094-1, 1999.

Watson, J.G., Chow, J., C., Lurmann, F. W., and Musarra, S. P.: Ammonium Nitrate, Nitric Acid, and Ammonia Equilibrium in Wintertime Phoenix, Arizona, Air &\#x26; Waste, 44:4, 405-412, DOI: 10.1080/1073161X.1994.10467262, 1994.

Womack, C. C., McDuffie, E. E., Edwards, P. M., Bares, R., de Gouw, J. A., Docherty, K. S., Dube, W. P., Fibiger, D. L., Franchin, A., Gilman, J. B., Goldberger, L., Lee, B. H., Lin, J. C., Long, R., Middlebrook, A. M., Millet, D. B., Moravek, A., Murphy, J. G., Quinn, P. K., Riedel, T. P., Roberts, J. M., Thornton, J. A., Valin, L. C., Veres, P. R., Whitehill, A. R., Wild, R. J., Warneke, C., Yuan, B., Baasandorj, M., and Brown, S. S.: An odd oxygen framework for wintertime ammonium nitrate aerosol pollution in urban areas: NOx and VOC control as mitigation strategies, Geophys. Res. Lett., 46, 4971–4979, https://doi.org/10.1029/2019GL082028, 2019.

Xing, J., Ding, D., Wang, S., Zhao, B., Jang, C., Wu, W., Zhang, F., Zhu, Y., and Hao, J.: Quantification of the enhanced effectiveness of NOx control from simultaneous reductions of VOC and NH3 for reducing air pollution in the Beijing–Tianjin–Hebei region, China, Atmos. Chem. Phys., 18, 7799–7814, https://doi.org/10.5194/acp-18-7799-2018, 2018.
## Annex 1: Selection of seasons

Chemical regimes show a great variability with time. To select the extension of the seasons, we analyse the monthly behaviour of the NH₃ and NOₓ responses as well as the interaction terms (Figure 12). The analysis is performed for two locations: Bergamo which lies in a NH₃ sensitive zone and Mantova which lies within a NOₓ sensitive zone. On these plots, we identify two major seasons with a consistent behaviour: winter from November to February and summer, from April to September. These two seasons are similar at both locations. The two remaining months, March and October are transition months and are not considered in the analysis. It is interesting to note that these two transition periods correspond more or less to the switching time between the NOₓ and NH₃ concentration time profiles (Figure 12 - right panel). On the latter figures, the SO₂ and NO₂ temporal evolutions are almost identical, in contrast to NH₃.

## Annex 2: A theoretical example for the isopleths

To facilitate the interpretation of the isopleths diagrams, we use a simple theoretical example that mimics the complex reactions process schematized by equations (5) and (6) above. Our simplified process is described by the following relation: $C_c(x, y) = \min(E_A[x], E_B[y])$ where $C_c$ is the concentration of a compound “c” that is given by the minimum between two emitted species A and B. The concentration depends on the strengths of these emissions, specified by the parameters x and y. For each emission strength (x or y), the two emission species are proportional to their full-scale value (100%): $E_A(x) = xE_A[100]$ and $E_B(y) = yE_B[100]$, respectively. If we choose $E_A(100) >> E_B(100)$, we create a B-sensitive environment (Figure 13 - left column) and inversely (Figure 13 - right column). If we select mixed situations, representing for example an average of days, during which we alternate between A- and B-sensitive regimes, we obtain the two bottom isopleth diagrams that represent cases where a larger number of A-sensitive (right) or B-sensitive (left) events are recorded. Although extremely simple, these diagrams illustrate properties that are observed on the real test cases.

Let's take the example of a A-sensitive regime. Similar observations can be made in the case of a B-sensitive regime. We note the following points:

* The diagram area can roughly be divided into two zones separated by a ridge: a top-left triangle where the sensitivity to emission reductions of species “B” dominate and a bottom-right triangle where the sensitivity to A dominates. The slope of the ridge (larger or less than one) informs on the type of regime.

* In the case of a single A-sensitive day (top right) with $E_B(100)=2E_A(100)$, the concentration of compound C remains unchanged for emission reductions of B up to 50% while its concentration react in a linear way to emission changes of A from 0 to 100%. Between the base case and a reduction level of 50%, the A-gradient is therefore larger than the B-gradient. This implies that the B-gradient is larger than the A-gradient between the 50 and 100% reduction levels because we know that for a 100% reduction of A or B, the concentration must be zero. In our simple example, the gradient of B is zero from 0 to 50% but is twice as large as A between 50 to 100%.

* While the combination of several events (e.g. days) characterised by different regimes leads to smoother isopleths (bottom), the same characteristics can be noted. In particular, the inclination (tending to the horizontal or vertical) provides information on the type of chemical regime.
<figure><img src="image_2.png" /><figcaption>Figure 1: Seasonal variations of the base case PM<sub>2.5</sub> (in µg m<sup>-3</sup>, top row) and SIA (in % of the PM<sub>2.5</sub> concentration, bottom row). The symbols: Be, Ma and Bo represent the locations selected for more detailed analysis: Bergamo (Be), Mantova (Ma) and Bologna (Bo). Other cities are indicated by their first letter for convenience: Venice, Milan, Turin and Genova.</figcaption></figure>

725

<figure><img src="image_3.png" /><figcaption>Figure 2: Summer (top) and winter (bottom) concentration base-case fields for SO<sub>2</sub> (left column), NO<sub>2</sub> (central column) and NH<sub>3</sub> (right column) expressed in µg m<sup>-3</sup>. The symbols: Be, Ma and Bo represent the locations selected for more detailed analysis: Bergamo (Be), Mantova (Ma) and Bologna (Bo). Other cities are indicated by their first letter for convenience: Venice, Milan, Turin and Genova.</figcaption></figure>

730
<figure><img src="image_4.png" /></figure>

Figure 3: Winter (right) and Summer (left) chemical regimes obtained from an emission reduction level of α=25%. The maps represent the $P_{NO_x}^{25\%} - P_{NH_3}^{25\%}$ (unit: µg m⁻³) indicator that shows the NOₓ- and NH₃-sensitive areas in yellow/red and blue, respectively. The symbols Be, Ma and Bo indicate the location of Bergamo, Mantova and Bologna, respectively. Other cities are indicated by their first letter for convenience: Venice, Milan, Turin and Genova.

735

<figure><img src="image_5.png" /></figure>

Figure 4: G-ratio for winter (right) and summer (left) times. The symbols Be, Ma and Bo indicate the location of Bergamo, Mantova and Bologna, respectively. Other cities are indicated by their first letter for convenience: Venice, Milan, Turin and Genova.

740
<figure><img src="image_6.png" /><figcaption>Figure 5: Yearly averaged chemical regimes obtained from a 25% emission reduction starting at different levels of emissions corresponding to α=0, 25, 50 and 75. The maps represent the ($P_{NO_x} - P_{NH_3}$) between the starting and ending levels (unit: µg m⁻³) showing the NOx-and NH₃-sensitive areas in yellow/red and blue, respectively. The symbols Be, Ma and Bo indicate the location of Bergamo, Mantova and Bologna, respectively. Other cities are indicated by their first letter for convenience: Venice, Milan, Turin and Genova.</figcaption></figure>

745

750
<figure><img src="image_7.png" /><figcaption>Figure 6: PM<sub>2.5</sub> isopleths during Winter (top row) and Summer (bottom row) at the three locations of interest (see maps). PM<sub>2.5</sub> concentrations are expressed in µg m<sup>-3</sup> as a function of the intensity of the NO<sub>x</sub> (Y-axis) and NH<sub>3</sub> (X-axis) emissions, respectively. The overlaid dashed oblique lines on each diagram connect similar PM<sub>2.5</sub> concentration values for single NO<sub>x</sub> and NH<sub>3</sub> reductions. The more vertical are these lines, the larger is the NH<sub>3</sub> abatement impact compared to the NO<sub>x</sub> abatement impact; The more horizontal they are, the larger is the NO<sub>x</sub> abatement impact compared to the NH<sub>3</sub> abatement impact. The dashed line delineates the ridge between the NH<sub>3</sub> and NO<sub>x</sub> sensitive areas.</figcaption></figure>

755
<figure><img src="image_8.png" /><figcaption>Figure 7: Wintertime isopleths in Bergamo for the species: dry PM<sub>2.5</sub>, O<sub>3</sub>, NO<sub>2</sub>, sulfate (SO<sub>4</sub>), nitrate (NO<sub>3</sub>), ammonium (NH<sub>4</sub>), organic matter (OM), anthropogenic and biogenic secondary aerosols (ASOA and BSOA, respectively). The two numbers on the vertical axis indicate concentration values for the base case and at the 50% NOx emission level.</figcaption></figure>

760
<figure><img src="image_9.png" /></figure>

Figure 8: Yearly mean non-linear interaction term for a 25% reduction in NOₓ and NH₃. The scatter plot compares the sum of the potential impacts of the two-single reductions (X-axis) with the potential impact of the combined emission reduction (unit: µg m⁻³) (Y-axis). Departure from the 1:1 line quantifies the overall non-linearity. Each data point represents the yearly average values for a grid cell in the modelling domain.

765

<figure><img src="image_10.png" /></figure>

Figure 9: Wintertime maps of the overall non-linearity term (top left) and of its components expressed as potential impacts between 25 and 50%: the single NOₓ non-linearity term (top right), the single NH₃ non-linearity term (bottom

770
right) and the NOₓ-NH₃ interaction term (bottom left). The three locations of interest: Bergamo, Mantova and Bologna are indicated by their first two letters while other cities are indicated with their first letter for convenience (Venice, Milan, Turin and Genova). The hand drawn contours roughly indicate the central position of the NH₃ (blue) and NOₓ (red) sensitive regime areas (see Figure 3 - right).

<figure><img src="image_11.png" /><figcaption>Figure 10: same as Figure 9 but for summertime.</figcaption></figure>
<figure><img src="image_12.png" /></figure>

Figure 11: Changes of the overall non-linearity terms from 25 to 50% (top), from 50 to 75% (middle) and from 75 to 100% (bottom) in NOₓ and NH₃. The overall non-linearity is visualised as the distance from the 1:1 diagonal, i.e. the difference between the overall potential impact at two levels of emission reduction, X- and Y-axis. Each point represents one land grid cell within the domain for wintertime (blue) and summer time (red). The “fit” parameter indicates the slope of the regression line while R² and RMS provide information on the coefficient of determination and the root mean square error, respectively.

780

785
<figure><img src="image_13.png" /><figcaption>Figure 12: Monthly averaged responses to NH₃ (blue), NOₓ (red) reductions (25%) and interaction terms (grey) at two locations: Bergamo (top) and Mantova (bottom). The right panel shows the monthly evolution of the concentrations of NO₂, NH₃ and SO₂ at those two locations. Note that concentrations are normalised by their average values.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Figure 13: Isopleths for a simple theoretical system consisting of two emission precursors (A and B) competing through non-linear reactions to the concentration of a pollutant. See details in the text.</figcaption></figure>