Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

An analysis of 30 years of surface ozone concentrations in Austria: temporal evolution, changes in precursor emissions and chemical regimes, temperature dependence, and lessons for the future

Monika Mayer *a, Stefan F. Schreier a, Wolfgang Spangl b, Christoph Staehle a, Heidelinde Trimmel a and Harald E. Rieder a
aInstitute of Meteorology and Climatology, Department of Water, Atmosphere, and Environment (WAU), University of Natural Resources and Life Sciences (BOKU), Vienna, Austria. E-mail: monika.mayer@boku.ac.at
bEnvironment Agency Austria, Vienna, Austria

Received 11th January 2022 , Accepted 12th April 2022

First published on 19th April 2022


Abstract

Despite substantial reductions in anthropogenic emissions of nitrogen oxides (NOx) and non-methane volatile organic compounds (NMVOCs) in Austria over the 30 year time period 1990–2019, summertime surface ozone (O3) concentrations still exceed frequently and over wide areas the ozone maximum 8 hour mean target value for the protection of human health. We present a detailed analysis of in situ observations of O3 and NOx to (1) disentangle the main processes propelling O3 formation such as precursor emissions and meteorology and (2) quantify the impact of NOx reductions and (3) estimate the effect of climate warming. The temperature sensitivity of surface O3 production is assessed separately for NOx and VOC limited regimes. The temperature sensitivity of ozone increases with temperature in spring and summer. On average, the evaluated absolute values of the sensitivities are a factor of 2.5 larger in summer than in spring. The analysis of ambient O3 burdens during hot summers indicates that rising temperatures in a warming climate might largely offset the benefit of future emission reductions. MAX-DOAS formaldehyde (HCHO) measurements in Vienna from 2017 to 2019 are used as a proxy for VOC emissions. The seasonal and the temperature dependence of the observed HCHO mixing ratios indicate that biogenic VOCs (BVOCs) are the dominant source of hydrocarbons in the urban setting during the ozone season. The result agrees well with VOC emission estimates that show BVOCs to be the dominant VOC source in Austria since the early 2000s. Accordingly, anthropogenic NOx emission reductions remain, outside of urban cores, the most effective instrument for policymakers to lower surface ozone concentrations in the short term.



Environmental significance

Elevated surface ozone concentrations continue to be of concern for public health in many countries, despite substantial reductions of anthropogenic precursor emissions during recent years. Progressive climate warming accompanied by increased biogenic volatile organic compound emissions propels ozone production. On a regional basis the contribution of changes in individual factors affecting surface ozone concentrations such as anthropogenic and natural emissions, solar radiation or temperature is not well constrained. In our work, we quantify the effect of NOx emission reductions, and the temperature sensitivity of surface ozone based on 30 years of observational data for NOx and VOC limited chemical regimes in Austria. Overall, we find that the climate penalty outweighs the gains achieved by emission reductions. Applying our analysis scheme would allow other research teams to investigate the climate penalty on ozone formation in other countries/world regions. Furthermore, our results provide a first-order estimate regarding the magnitude of changes in the surface ozone burden that can be expected in light of climate warming and emission changes over the next few decades, thereby also providing a context for results obtained with chemistry–climate and chemistry–transport models.

1 Introduction

Surface ozone is harmful at elevated levels in surface air to human health1–4 and vegetation, e.g. ref. 5–7, and also impacts surfaces of built structures.8 Furthermore, O3 is an important tropospheric greenhouse gas belonging to the group of short-lived climate forcers.9–12 The global average radiative forcing due to increases in tropospheric O3 since pre-industrial times is estimated to be +0.35 ± 0.15 W m−2.13 While the average atmospheric lifetime of free tropospheric O3 is in the order of weeks,14 it reduces to a few hours for surface O3 in areas with elevated precursor emissions, e.g. ref. 15, leading to a highly variable spatial and temporal ozone burden.

In surface air, O3 is produced through a series of chemical reactions involving sunlight and O3 precursors such as nitrogen oxides (NOx: NO and NO2), non-methane volatile organic compounds (NMVOCs), carbon monoxide (CO) and methane (CH4).15 These precursors are emitted from anthropogenic (e.g., traffic or industry) and natural (e.g., soil, trees, lightning) sources, e.g. ref. 4. Ozone production rates depend strongly on the chemical environment (i.e., NOx or VOC limitation) and ambient meteorological conditions (i.e., solar radiation, temperature, and humidity).15,16 The O3 loss occurs through deposition at the surface and uptake by plants, heterogeneous reactions involving aerosols, or O3 titration. The latter occurs close to combustion sources, e.g. in urban environments, during the night when excess amounts of NO react with O3 to form NO2, e.g. ref. 17. The highest O3 concentrations are typically found in or downwind of major municipalities, e.g. ref. 16, and not in the urban core (due to the aforementioned titration effects). Mountain sites also generally exhibit elevated O3 concentrations because of enhanced short-wave radiation, intrusions of stratospheric O3, e.g. ref. 18 and 19, and an increased chemical lifetime of O3 in the free atmosphere.

As surface O3 is harmful to human health, national governments have implemented the European ozone air quality legislation in recent decades.20 The compliance with air quality standards is evaluated with monitoring networks applying state-of-the-art precision measurement techniques. The European Union (EU) defines an O3 target value and a long-term objective for protecting human health in the directive 2008/50/EC.21 The ozone target value for the protection of human health is defined as the maximum daily eight-hour average value (MDA8), not to exceed 120 μg m−3 on more than 25 days per calendar year averaged over three years. Although the EU MDA8 target value lies well above the World Health Organization's (WHO's) air quality guideline (i.e., MDA8 not to exceed 100 μg m−3), it is frequently not attained in Austria despite substantial efforts to curb precursor emissions, and the subsequent reductions achieved.22 During the 3 year evaluation period, 2017–2019, the EU target value was exceeded at 57% of the Austrian monitoring sites on more than 25 days per year.23 Furthermore, in 2019, the 120 μg m−3 threshold value was violated at least once at every station in Austria.

To understand changes in surface O3, it is essential to disentangle the factors influencing O3 concentrations. Important factors are anthropogenic and biogenic precursor emissions, the meteorological conditions that affect precursor dispersion, and photochemistry. Further, the temperature dependence of surface ozone production is relevant in the context of global warming. Today it is well understood that ozone chemistry depends nonlinearly on the above-mentioned drivers.15

In this paper, we focus on the effect of reduced precursor emissions (particularly NOx emissions) and increasing temperatures (in spring and summer) on Austrian surface O3 concentrations during 1990–2019. We exclude data for the year 2020, given the widespread impact of COVID-19 related restrictions on anthropogenic emissions and thus ambient air quality in Austria. Our findings provide a comprehensive characterisation of the evolution of the past O3 burden in Austria. They also enable us to estimate potential future developments in surface O3 burdens in the light of projections of further declining NOx emissions22 and concurrent rising temperatures due to climate change.24,25

2 Data and methods

2.1. Air quality data and metrics

Nationwide air pollution monitoring, including O3 and NOx, started around 1990 in Austria. The current monitoring network is operated by the Environment Agency Austria (Umweltbundesamt) and the governments of the nine Austrian provinces.

For surface O3 the monitoring network comprised in 2019 a total of 99 background sites in rural (50), suburban (40) and urban (9) regions, as well as additional monitors at traffic (6) and industrial sites (4) (note that monitoring sites are classified according to the guidelines of the European Environment Agency, EEA).26 The highest short-term O3 values (1 hour mean values) are observed in the densely populated region around the capital Vienna; the highest long-term levels (8 hour-mean values, annual mean values) are recorded in medium and high altitude alpine regions.23

For the present analysis, we consider only background sites located at altitudes below 1500 m a.s.l. and stations that were continuously operated throughout 1990–2019. This selection yields a total of 54 background sites (6 urban, 20 suburban, and 28 rural), for which we calculate the MDA8 O3.

To address the issue of meteorological variability and its impact on surface ozone formation, we aggregate data over 10 year periods,27 resulting in three time periods 1990–1999, 2000–2009, and 2010–2019. This enables us to study changes (in the order of several μg m−3) that are smaller than the inter-annual fluctuations caused by the variability of ambient meteorology. We analyse for each monitor category (1) the number of days on which the 120 μg m−3 EU threshold for MDA8 O3 is exceeded on an annual basis and for the warm seasons (JJA and MAM) of the year, which show peak ozone abundance, and (2) the spring and summertime probabilistic 1-, 3-, and 5-year return levels (RLs) of MDA8 O3. To this end we fit a peak-over-threshold model,28 with the threshold u = 120 μg m−3, to seasonal MDA8 O3 data for the decadal time slices defined above. These n-year RLs represent the MDA8 O3 values, exceeded in probabilistic terms at least once in n-years over the decadal time slice considered.29 Furthermore, we investigate the impact of temperature (a climate penalty) on MDA8 O3 by contrasting results for full 10 year time periods and the same time periods excluding data of hot years.

Along with ozone, we analyse co-located measurements of NO and NO2 available from the Austrian air quality monitoring network. Continuous observations are available at 39 background sites during 1990–2019 (5 urban, 18 suburban, and 16 rural sites). For these data, we focus hereinafter on the daily average (DA) NOx mixing ratio (expressed in ppb).

A table of all O3 and NOx stations used in this study is provided in the ESI in Table S1.Fig. 1a shows the location of the selected ozone monitoring sites; here, circles with a white dot indicate sites with co-located NOx measurements. For illustrative purposes, we select another six representative stations (highlighted with red circles in Fig. 1a and marked with an asterisk in Table S1: one urban, two suburban and three rural sites, located in or in the vicinity of Austria's largest municipalities Vienna and Graz). We use these sites to illustrate the chemical processes driving O3 formation and decay in more detail.


image file: d2ea00004k-f1.tif
Fig. 1 (a) Overview of the location of background O3 monitoring sites (light grey refers to rural, dark grey to suburban and black to urban sites). Numbered sites marked with red circles are selected as representative examples for site categories. (b) Estimates of Austrian NOx, NMVOC and biogenic VOC emissions in 1990–2019. Blue numbers at the bottom of the panel indicate the ratio of NOx emissions, and yellow numbers indicate that of NMVOC emissions, in a particular 10 year time period relative to the base period 1990–1999.

By analysing the co-evolution of O3 (MDA8) and NOx (DA) under clear sky conditions we determine whether O3 chemistry is NOx or VOC limited. If the observed MDA8 O3 and DA NOx data are directly proportional, we deduce a NOx limited regime. If O3 reveals no correlation with NOx, we find O3 chemistry to be VOC limited.30 The distinction between the chemical regimes can only be achieved if the total amount of global radiation is large enough not to limit ozone production. For both chemical regimes, peak ozone is known to be temperature sensitive, but to different extents.31 Therefore, we include mean temperature (Tmean) as a covariate to obtain the temperature sensitivity of MDA8 O3 separately for spring and summer (see Section 4.3).

Formaldehyde (HCHO) is a short-lived oxidation product of VOCs (including methane – CH4) and can therefore be used as a proxy for VOC emissions, e.g. ref. 32 and 33. The HCHO background concentration is determined by the oxidation of CH4,34,35 a globally well-mixed trace gas, and other long-lived VOCs.33 Observed variations of HCHO concentrations originate either from VOC oxidation processes (secondary source) or primary HCHO emissions (mainly combustion processes). BVOCs are highly reactive, resulting in a short atmospheric lifetime in the order of hours.36,37 As a result, ambient HCHO concentrations are substantially influenced on short time and small spatial scales by BVOC emissions. Today it is well understood that BVOCs are an important HCHO source, even in urban areas.38 The partitioning between secondary and primary sources in urban areas is about 95% to 5%, e.g. ref. 39 and 40.

Since June 2017, HCHO has been measured continuously in Vienna with a series of MAX-DOAS instruments that are part of the VINDOBONA network (https://www.doas-vindobona.at/).41 Here we use data retrieved from the instrument, stationed on the roof of the University of Natural Resources and Life Sciences, Vienna (BOKU). We analyse MAX-DOAS retrievals in the UV wavelength range (338–370 nm) for a light path that stretches horizontally across Vienna's city center (median path length is 9 km) at an altitude of approximately 150 m above ground for solar zenith angles smaller than 75°. We then calculate from the available daytime retrievals the DA HCHO mixing ratio.

HCHO MAX-DOAS observations are further combined with DA NOx, MDA8 O3 and Tmean data for Vienna's city center (Stephansplatz) to investigate O3 chemistry at the urban background site. This enables us to derive O3 isopleths directly from observations, using HCHO as a proxy for VOCs. Additionally, we examine the relationship between Tmean and HCHO mixing ratios. HCHO observations in Vienna allow us to focus on three objectives, (1) to confirm our approach to determining the O3 chemical regime, (2) to derive an ozone isopleth plot purely from observations, and (3) to estimate the impact of BVOCs in the largest metropolitan region in Austria.

2.2. Evolution of precursor emissions

In tandem with emissions in most European regions, Austrian anthropogenic NOx and VOC emissions dropped significantly between 1990 and 2019.4,42Fig. 1b provides an overview of the Austrian national annual total anthropogenic emissions of NOx and NMVOCs over the last few decades.22 The gridded data we analyse here are available at the EMEP Centre on Emission Inventories and Projections (https://www.ceip.at/the-emep-grid/gridded-emissions). BVOC emissions are derived from various models such as Bauwens et al. (2018)43 for isoprene (https://emissions.aeronomie.be; BIRA IASB, 2018) and Simpson and Winiwarter (1998)44 for other BVOCs. Two data sets for NOx and NMVOCs are given: one provides emission estimates based on fuel sold (dashed line) and one based on fuel used (solid line). The latter is used as a reference in our analysis as it represents the better proxy for actual emissions in Austria. We note that the difference between the two scenarios is minimal for NMVOCs since traffic is not a dominant source. Austrian NOx emissions decreased between 1990 and 2019 by about 35%, while NMVOC emissions decreased by about 66%. If we analyse average NOx emissions over the 10 year periods, we find in 2000–2009 a decline of only 2% and in 2010–2019 a decline of 18% compared to the reference time-period 1990–1999 (for NMVOCs: −33% in 2000–2009 and −53% in 2010–2019, respectively).

NMVOC emissions reached in 2002 approximately the same level as estimated BVOC emissions (see Fig. 1b). The error bar given in Fig. 1b for the year 2006 is derived from the work of Skjoth et al. (2008)45 and Oderbolz et al. (2013).46 The isoprene share of BVOC emissions was determined by a modelling effort driven by reanalysis data and constrained by HCHO satellite observations.32,43 Monoterpenes and oxygenated volatile organic compound emissions are assumed to be constant over time.44 Since 2014 BVOC emissions have been thought to be twice as large as national anthropogenic NMVOC emissions.22 Concerning surface O3, VOC emissions during spring and summer are most relevant. It was estimated that about 80% of the annual emissions occur during summer.33,47–49 Future Austrian anthropogenic NMVOC emissions are projected to reach levels slightly lower than today,22,50 while BVOC emissions are likely to increase due to their temperature dependence, e.g. ref. 48, 51 and 52. As Austrian NOx emissions are expected to decrease substantially in the coming decades,53,54 it is crucial to understand the response of surface ozone concentrations to these reductions while anthropogenic NMVOC emissions remain essentially unchanged. Hence, we focus explicitly on the time-period 2010–2019, for which we assess ozone's temperature sensitivity at rural, suburban, and urban sites.

To investigate the correlation between annual emission estimates and observed DA NOx mixing ratios during the prime ozone season (March to August) we analyse the ECDFs (empirical cumulative distribution functions) of measured DA NOx. The median of observed DA NOx mixing ratios decreased across all monitor categories by 30–37% during the spring and summer season in 2010–2019 compared to 1990–1999. This decline in the 10 year average observed DA NOx mixing ratio is larger than the corresponding emission reduction (−18%). In the ESI in Fig. S1, we show the ECDFs of observed DA NOx mixing ratios for spring (MAM) and summer (JJA) for 1990–1999, 2000–2009 and 2010–2019. For completeness, Table S2 in the ESI lists the observed changes for the 50th and 90th quantiles for each site category for MAM and JJA, respectively.

2.3. Meteorological data

As meteorological covariates, we consider air temperature (2 m above ground) and total global radiation (i.e., the daily sum of global radiation), available from gridded data sets (SPARTACUS (1 km × 1 km grid)55,56 and APOLIS (100 × 100 m grid)57) provided by the Austrian meteorological service Zentralanstalt für Meteorologie und Geodynamik (ZAMG). Given the importance of photochemistry for O3 formation, we specify radiation levels by selecting clear sky days when investigating the impact of NOx and temperature on O3 concentrations. To quantify radiation, we calculate the daily total global radiation for each monitoring site. The annual cycle of the daily sums of the total global radiation is plotted, and the envelope of the 30 years of observational data is fitted with a polynomial of the 5th order. This envelope-function approximates our maximum feasible total global radiation for each day of the year. “Clear sky” conditions are assumed if the observed total daily global radiation reaches 90% of the polynomial fit on a specific day. The selection of clear sky days indicates that mean and maximum temperatures are correlated. We choose the daily mean temperature (Tmean) at the site level from the gridded data set as variable because we anticipate a significant variance in the ozone productivity solely due to the difference in prevailing temperatures in spring and summer.

3 Results

3.1. Evolution of MDA8 O3 between 1990 and 2019

While hourly peak O3 values have substantially decreased over the last three decades,23 exceedances of the MDA8 O3 threshold of 120 μg m−3 are still regularly observed nationwide, particularly during the summer season. Fig. 2a illustrates the time evolution of the annual average number of exceedance days of the MDA8 O3 target value for the protection of human health for urban, suburban, and rural monitoring sites. The fewest exceedances are observed at urban sites, the most at rural stations. Overall exceedances at all site categories show a negative trend between 1990 and 2019. Despite some difference in the precise time periods considered, our results show a strong agreement with those provided by Chang et al. (2017)58 for European sites contained in the tropospheric ozone assessment report (TOAR) database. In Fig. 2b, we plot the percentage of sites exceeding the MDA8 O3 in each category. The apparent large year to year variability can be attributed to different meteorological conditions. The high number of exceedance days in 2003, 2015, 2017, 2018 and 2019 is related to the five hottest summers (JJA) in Austria since the beginning of national records in 1767 (deviation in summer temperature relative to the average for 1981–2010 is in 2003: +2.8 °C, 2019: +2.7 °C, 2015: +2.4 °C, 2017: +2.1 °C, and 2018: +2.0 °C (ref. 59)). As high temperatures propel ozone production, their occurrence is directly reflected by an increased number of days above the MDA8 O3 threshold value. To study the influence of hot seasons such as the aforementioned hot summers, we exclude these years in an additional assessment, i.e. for 2000–2009 we omit 2003, and for 2010–2019 we omit 2015, 2017, 2018, and 2019, respectively.
image file: d2ea00004k-f2.tif
Fig. 2 (a) Annual average number of exceedance days of the MDA8 O3 threshold value of 120 μg m−3 for each category of monitoring sites in Austria. (b) The percentage of sites for each category that exceeded. (c–e) Boxplots of the annual number of exceedance days of the MDA8 O3 threshold value for urban, suburban, and rural background sites for 1990–1999, 2000–2009 and 2010–2019. Boxplots in orange colour in (d) and (e) exclude measurements for 2003 and 2015, 2017, 2018, and 2019, respectively. The colour-coded solid lines in panel (a) provide linear trend estimates per site category. The dotted line in panels (a) and (c)–(e) indicates the target threshold of 25 days.

In Fig. 2c–e, we show the spread in the annual average number of exceedances per site location category (urban, suburban, and rural background) for all three selected time periods. The evaluations without data for the hot years are plotted in orange, respectively. Most outliers (at all site location categories) in the MDA8 O3 exceedances during 2000–2009 and 2010–2019 occurred during hot summers (Fig. 2d and e).

Overall, we find an increase in exceedance frequencies at all monitoring stations during 2000–2009 compared to the 1990–1999 period. At rural monitors, we observe a decrease in 2010–2019 (on average minus 7 exceedance days) relative to the base period, whereas at suburban and urban sites, the average number of exceedance days converged to the limit value of 25 days. If we omit years with hot summers from our analysis (orange boxplots in Fig. 2d and e), the number of exceedance days in the period 2000–2009 is slightly higher compared to 1990–1999 but substantially lower in the period 2010–2019 (median: minus 5 days at suburban and urban, minus 13 days at rural sites). The actual magnitude of the change in surface O3 is defined by the underlying ozone production mechanisms. These are mainly VOC limited in central urban areas, and NOx limited in rural areas during summer. The NOx to VOC ratio is predominantly altered during heat waves by increased BVOC emissions caused by the heat stressed plants (e.g. observed isoprene mixing ratios in the UK were found to be a factor of 4 higher during heat spells60). Heatwaves are commonly accompanied by stagnant, clear sky conditions, leading to elevated radiation levels under which photochemical reactions are accelerated. Furthermore, under such conditions, O3 dry deposition typically diminishes because plants close their stomata and reduce therefore the O3 uptake.61

At rural monitoring sites, i.e., in NOx limited areas, we expect NOx emission reductions to translate directly into a reduction in the number of exceedance days. The observed q50 of DA NOx mixing ratios at rural sites was reduced by 36% between 1990–1999 and 2010–2019 (see Table S2) during the ozone season. The number of days with MDA8 O3 exceedances decreased for the final analysis period by 18% (minus 7 days) at rural sites reflecting the adverse impact of a warming climate. If we exclude years with hot summers, the median of exceedance days reduces in 2010–2019 by an additional 6 days, yielding a total of minus 13 days (−33%). This reduction is very close to the 36% decrease in DA NOx mixing ratios that we derive from measurements at rural sites.

Precursor concentrations at suburban sites represent the transition zone between NOx and VOC limited regimes, depending on pollutant transport (enhanced NOx concentrations downwind of urban cores) and VOC levels. The similarity of the evolution of exceedance days at suburban and urban monitors indicates that peak MDA8 O3 concentrations occur primarily in the VOC limited regime.

To study the time evolution of the MDA8 O3 and not only exceedances of 120 μg m−3, we investigate spring (March, April, May – MAM) and summer (June, July, August – JJA) separately. O3 chemistry may differ seasonally at a particular site due to the impact of BVOCs with their seasonal cycle. Changes in the seasonal MDA8 O3 burden mimic the temporal evolution of exceedance days discussed above. On average, surface ozone concentrations are elevated in 2000–2009 for both seasons compared to 1990–1999 throughout all categories of background sites. Increasing O3 abundances in 2000–2009 have been followed by an overall decline of MDA8 O3 concentrations in 2010–2019 in MAM as well as JJA. The median MDA8 O3 levels at rural sites in the last decadal time slice match concentrations observed in the base period (in spring and summer). For suburban and urban monitors, we find slight increases in the median values in the order of 2–5 μg m−3 (4 μg m−3 in JJA for suburban and urban, 5 μg m−3 in MAM for urban, and 2 μg m−3 in MAM for suburban sites) between the periods 1990–1999 and 2010–2019. This rise in the MDA8 O3 median values is primarily driven by increases in the lower tail of the MDA8 O3 distribution. In contrast, changes in the upper tail are less pronounced (see Fig. S2 in the ESI for the evolution of the MDA8 O3 concentrations for MAM and JJA).

At all site categories, the range of concentrations decreased in the middle and last 10 year time period with respect to the first. This can be attributed to the reduction of precursor emissions resulting in narrower O3 distribution functions (see Fig. S2, where the length of whiskers in the boxplots becomes shorter with time). The most pronounced peak MDA8 O3 reductions are found at rural sites where ozone production is NOx limited (maximum MDA8 O3 ∼minus 15 μg m−3 in both seasons in the last period compared to the first period). The analysis without hot years allows us to attribute changes linked to higher temperatures. We find an increase of the median MDA8 O3 of about 2 μg m−3 in 2000–2009 and about 6 μg m−3 in 2010–2019 because of the enhanced ozone productivity at elevated summertime temperatures. In spring the temperature effect is much smaller, 0–1 μg m−3 in 2000–2009 and 1–3 μg m−3 in 2010–2019 despite spring temperatures in 2003, 2015, 2017, and 2018 being well above the climatological mean of 1981–2010 (2003: +1.5 °C, 2015: +1.5 °C, 2017: +2.4 °C, 2018: +2.7 °C, 2019: +0.4 °C (ref. 59)). Elevated temperatures affect peak ozone during spring to a much smaller extent than during the summer season because the whole temperature range is too low for efficient ozone production. A detailed temperature sensitivity analysis of MDA8 O3 is presented in Section 4.3.

To assess changes in peak MDA8 O3 in more detail, we apply methods from extreme value theory (EVT) in the next section.

3.2. Extreme value analysis – return levels of MDA8 O3

Since summer months are most prone to high surface O3 levels, we focus in the extreme value assessment on MDA8 O3 during JJA. Applying a peak-over-threshold (POT) model with a threshold of 120 μg m−3, we derive the probabilistic MDA8 O3 return levels at intervals of 1, 3, and 5 years.

The RLs corresponding to the MDA8 O3 value exceeded at least once in 1, 3 or 5 years, respectively, thus illustrating the changes in peak ozone abundances. In Fig. 3, boxplots of the estimated summertime MDA8 O3 RLs for urban, suburban, and rural sites are shown for 2000–2009 and 2010–2019. The median RLs are generally highest in rural areas, while suburban sites show the largest spread in site-level RLs.


image file: d2ea00004k-f3.tif
Fig. 3 Boxplots of 1-, 3-, and 5-year return levels for JJA MDA8 O3 in 2000–2009 for (a) urban, (b) suburban and (c) rural sites. (d)–(f) as (a)–(c) but for 2010–2019. Boxplots in orange colour refer to the analysis without the hot years: 2003 (a–c) and 2015, 2017, 2018, and 2019 (d–f).

To estimate the impact of heatwaves on MDA8 O3 RLs, we analyse and contrast MDA8 O3 RLs for 2000–2009 and 2010–2019, including/excluding the years with hot summers (2003, 2015, 2017, 2018, and 2019). In Fig. 3, the orange boxplots indicate the analysis without the extreme years with generally smaller MDA8 O3 RLs. The 1 year RLs are lowered on average by ∼4 μg m−3 in 2000–2009 and by ∼6 μg m−3 in 2010–2019 when excluding hot years. This illustrates the temperature penalty on O3 during hot summers. The median 1 year RLs of the decadal time slices decreased by ∼4 μg m−3 at urban, by ∼2 μg m−3 at suburban and by ∼11 μg m−3 at rural sites in 2010–2019 compared to 2000–2009. For higher-order return levels, decreases particularly at rural and suburban sites are even more pronounced, yielding up to ∼6 μg m−3 for 3 year RLs and ∼14 μg m−3 for 5 year RLs.

4 Relationship between MDA8 O3, DA NOx, total global radiation, and daily mean surface temperature

As we focus on peak O3 concentrations, we study ozone chemistry only under clear sky conditions (see the definition in Section 2.3). Observations of NOx mixing ratios, total global radiation, and mean temperature allow us to investigate ozone chemistry at the site level and its temperature dependence in greater detail. The subsequent analysis is performed on a seasonal basis for MAM and JJA for the period from 2010 to 2019, as this time period is most informative for potential future changes in surface ozone following emission controls. During this period, annual anthropogenic NMVOC emissions decreased slightly and approached a value comparable to (near-)future estimates.22 In a first approximation, we can assume anthropogenic NMVOC emission to be constant on an annual basis. Hence the analysis of this period allows us to study the impact of NOx reductions on ozone concentrations which will be relevant for future scenarios. First, we select 6 representative sites (see Fig. 1a, red circled stations): one urban, two suburban, and three rural sites in close proximity in two Austrian regions, to demonstrate the approach. Subsequently, we use the results of all monitors to generate a contour plot of MDA8 O3 temperature sensitivity with respect to observed NOx mixing ratios and Tmean.

4.1. MDA8 O3vs. DA NOx

Scatterplots of (all year) MDA8 O3 and DA NOx data for 2010–2019 are shown in Fig. 4 for selected sites, colour coded according to radiation abundances. Note the wide range in DA NOx mixing ratios within each site category, specifically for suburban and rural sites. Generally varying mixing layer heights determine NOx mixing ratios at each site throughout the year. The red dotted line in all panels indicates the EU threshold value of 120 μg m−3 for MDA8 O3.
image file: d2ea00004k-f4.tif
Fig. 4 Scatterplots of MDA8 O3versus DA NOx for (a) the urban station Vienna Stephansplatz, (b and c) suburban sites Hermannskogel and Graz Nord, and (d–f) rural sites Illmitz, Klöch, and Masenberg for 2010–2019. The red dotted line in all panels indicates the MDA8 O3 threshold value of 120 μg m−3. Colour coding is according to daily total global radiation abundances.

Across sites, we find enhanced ozone production if total global radiation exceeds 4 kW per day per m2 (3.46 MJ m−2, yellow to red data points). Conversely, for low radiation values (blue to green data points), the observed MDA8 O3 is determined by the ambient O3 background concentration, or O3 is even found to be depleted if the monitor is close to NOx emission sources (titration). O3 titration is seen at the urban site Vienna (Fig. 4a) and the suburban site Graz (Fig. 4d) for days with total global radiation below 2 kW per day per m2 (1.73 MJ m−2).

4.2. Clear sky MDA8 O3 chemistry

We filter measurements with respect to clear sky conditions to assess if ozone chemistry is NOx or VOC limited. In the case of NOx limited chemistry, we expect ozone concentrations to follow NOx proportionally. Conversely, for VOC limited conditions, we anticipate no relationship between O3 and NOx.30 For both regimes, ambient temperature is an important covariate. Fig. 5 shows the seasonal (MAM and JJA) scatterplot of MDA8 O3versus DA NOx, colour-coded by Tmean, for an urban (Vienna Stephansplatz), a suburban (Hermannskogel), and a rural (Illmitz) site for 2010–2019. The corresponding illustrations for the other 3 exemplary stations – Graz Nord, Klöch, and Masenberg – are given in ESI Fig. S3.
image file: d2ea00004k-f5.tif
Fig. 5 MDA8 O3versus daily average (DA) NOx during summer (JJA) and spring (MAM) under clear sky conditions for (a) the urban station Vienna Stephansplatz, (b) the suburban site Hermannskogel, and (c) the rural site Illmitz. The red dotted line in all panels indicates the threshold of 120 μg m−3. MDA8 O3 is coloured according to Tmean. The grey lines are polynomial fits of the 4th order to the binned data envelopes.

To indicate the chemical regimes, splines are calculated as polynomial fits of the 4th order to the envelopes of the scatterplots of MDA8 O3versus DA NOx. For the envelope fitting, we bin the observations in 20 equally sized bins using data between the 10% and 90% quantiles (to avoid introducing arbitrary effects by outliers) and fit to bin maxima. The splines obtained are plotted in Fig. 5 (and in the ESI in Fig. S3) as grey lines. For the urban station in Vienna Stephansplatz (Fig. 5a), three domains based on NOx levels can be identified during spring: low-range NOx (DA NOx < 8 ppb) which is NOx limited, medium-range NOx (8 ppb < DA NOx < 40 ppb) which is VOC limited, and high-range NOx (DA NOx > 40 ppb) where ozone concentrations decline with increasing NOx (ozone titration because of relatively low absolute radiation in early spring). During summertime, ozone chemistry is predominantly VOC limited. In this regime, the highest ozone concentrations are observed. Only for DA NOx below 8 ppb NOx limitation occurs (∼20% of the observations). Generally, rising temperatures facilitate enhanced O3 production, i.e., for a given DA NOx value the MDA8 O3 yield increases with rising temperature.

For suburban sites, we expect a mix of VOC and NOx limited conditions from theory. Fig. 5b provides seasonal scatterplots of MDA8 O3versus DA NOx for the monitoring site Hermannskogel. Ozone production during MAM appears mainly VOC limited (4 ppb < DA NOx < 18 ppb). However, at very low NOx concentrations (DA NOx < 4 ppb) the chemical regime shifts towards NOx-limitation. In summer, NOx concentrations are found to be significantly lower than in MAM. As a result, NOx limited ozone chemistry becomes more important. About 65% of the observations reside in the NOx limited regime in JJA. In Fig. 5c, seasonal scatterplots for the rural station Illmitz are shown. Here a NOx limitation of ozone chemistry is observed during summertime, while NOx limitation emerges in spring only for very low NOx mixing ratios (<4 ppb). For DA NOx greater than 5 ppb, ozone chemistry is VOC limited, indicating the limited availability of BVOCs at that time of the year.

For the metropolitan area of Vienna, MAX-DOAS HCHO observations (starting 06/2017) are available. This data allows us to analyse the correlation of HCHO mixing ratios and Tmean. If BVOC emissions are the dominant source, observed HCHO mixing ratios are expected to rise with temperature during the growing season. In Fig. 6a, we plot DA HCHO mixing ratios versus Tmean for clear sky days (time-period 06/2017–12/2019). Data are fit with a linear model. From April to October, observed MDA8 O3 concentrations increase with temperature. From November to March, ozone and HCHO are negatively correlated. Colder days and lower mixing layer heights result in higher HCHO mixing ratios. This illustrates that BVOCs control hydrocarbon chemistry in Vienna's city centre during the growing season, which overlaps with the ozone season. The HCHO measurements allow us also to investigate ozone chemistry directly by assessing ozone isopleth plots. Spring and summertime measurements of HCHO (MAX-DOAS), NOx, and O3 (in situ, urban site Stephansplatz) are combined to obtain the O3 isopleth plot displayed in Fig. 6b. This illustration clearly reveals that peak MDA8 O3 concentrations are in the VOC limited regime. Most of the data employed for the isopleth plot stem from hot springs and summers (with average temperature lying more than 2 °C above the climatological mean for 1981–2010, except spring 2019, which was only +0.4 °C warmer). Hence the isopleths may provide some guidance regarding estimates for future ozone chemistry in Vienna in a warming climate.


image file: d2ea00004k-f6.tif
Fig. 6 (a) Temperature (Tmean) dependence of observed HCHO mixing ratios in Vienna (city centre) colour-coded according to MDA8 O3 for the period from 06/2017 to 12/2019. (b) MDA8 O3 isopleth plot for MAM and JJA derived from path averaged DA HCHO mixing ratios (approx. 9 km across Vienna's city centre) and in situ O3 and NOx measurements.

4.3. MDA8 O3 temperature sensitivity

The temperature sensitivity of the observed surface ozone concentrations is predominantly determined by the temperature sensitivity of chemical kinetics and the temperature dependence of BVOC emissions. Increased BVOC concentrations due to higher temperatures will lead to enhanced ozone concentrations in VOC limited areas, e.g. ref. 62, such as the Vienna metropolitan region. Based on our results for Vienna, we assume BVOCs to be the dominant source in the other urban environments across Austria.

The temperature dependence of the ozone kinetics is determined by the absolute NOx concentration.31 It is weak for the low NOx range (NOx limited) and pronounced for the VOC limited range. For both regimes, a non-linear relation between ozone production and temperature is observed. To determine ozone's temperature sensitivity, we plot MDA8 O3versus Tmean for each site for clear sky days and fit the data with a linear model (LM). The linear fit works well because we distinguish between seasons, MAM and JJA, respectively. Therefore, the whole temperature range is partitioned into segments, and within each segment, Tmean and MDA8 O3 show a linear correlation. For convenient reference, we provide scatterplots with the linear fits for the 6 sample stations in ESI Fig. S6 for MAM and JJA. LM slopes are generally larger in JJA than in MAM, the broadest range of observed temperature sensitivities (2.5 to 5 μg per m3 per °C Tmean) is found at suburban sites in JJA, and the highest values occur at urban sites in JJA (6 μg per m3 per °C Tmean).

In the next step, the temperature sensitivity of MDA8 O3 defined by the LM slopes is related to the median of the corresponding DA NOx mixing ratios. The contour plots in Fig. 7 summarise the slopes of the linear fits – changes of the MDA8 O3 in μg °C−1 m−3 – for both seasons for rural, suburban, and urban sites for 2010–2019. In spring, the temperature sensitivity is much smaller than in summer (∼2 μg °C−1 m−3). During summer, absolute values at all site categories are found in the range of 4.4–5.6 μg °C−1 m−3 but reveal a very different correlation with respect to DA NOx mixing ratios.


image file: d2ea00004k-f7.tif
Fig. 7 MDA8 O3 temperature sensitivities for rural, suburban, and urban monitors for various DA NOx mixing ratios and mean temperatures. The ozone temperature sensitivities for 2010–2019 under clear sky conditions are shown for spring (MAM) in the top panels (a–c) and for summer (JJA) in the bottom panels (d–f). The contours of the temperature sensitivity of ozone are given in μg °C−1 m−3 with respect to changes of Tmean. Plots for 1990–1999 and 2000–2009 can be found in ESI Fig. S8 and S9.

At rural sites under low NOx conditions (DA NOx < 5 ppb), the MDA8 O3 sensitivity is constant for all observed Tmean values, whereas for DA NOx > 5 ppb, the temperature sensitivity increases with Tmean. In the low NOx range, as defined by Vogel et al. (1999),31 the dominant temperature-sensitive reaction is the release of NO2 from peroxyacetyl nitrate (PAN). In our study, rural monitoring sites can be predominately ascribed to the low NOx range. Hence, the observed temperature sensitivity at rural sites can be understood by considering the decomposition of PAN, and the VOC to NOx ratio changes due to increasing BVOC emissions with rising temperatures.

The altered NOx/VOC ratio with the additional BVOCs leads to a higher ozone production rate and therefore enforces the effect of PAN decomposition. For suburban monitors, the temperature sensitivity generally increases with temperature, reflecting the temperature dependence of VOC dominated chemistry31,63 and the temperature dependence of BVOC emissions. The observed DA NOx mixing ratios for suburban monitors (>5 ppb) also indicate mainly VOC limited ozone chemistry, as illustrated in Fig. 5b for Hermannskogel, and Fig. S3a for Graz Nord. The strongest increase in the temperature sensitivity (25%) is observed at rural sites with DA NOx mixing ratios ranging from 5 ppb to 20 ppb. For Tmean ∼10 °C, the sensitivity is found to be 4.4 μg °C−1 m−3 and rises for Tmean ∼20 °C to 5.4 μg °C−1 m−3 in summer.

A one °C increase in Tmean results on average in a 5 μg m−3 increment in MDA8 O3. This is a substantial increase, particularly considering the decrease of MDA8 O3 median concentrations achieved due to precursor emission reductions since 1990, which are of the same order of magnitude (2–5 μg m−3 when hot years are excluded from the analysis, see Fig. S2d–f). Within the last 30 years, Austrian NOx concentrations declined by ∼30–40%. Simultaneously, NMVOC emissions are estimated to have decreased by ∼60%. As illustrated in Fig. 2e, the series of hot summers in the period from 2010 to 2019 offset almost completely the MDA8 O3 reductions obtained by decreasing precursor emissions. Our results indicate that the (MDA8) O3 temperature sensitivity will pose a continuing challenge for air pollution regulators in suburban and urban environments, where the ozone production regime is VOC limited and (temperature-dependent) BVOCs are nationwide the dominant VOC source.

In the ESI in Fig. S8 and S9, the temperature sensitivity of ozone is provided for the time periods 1990–1999 and 2000–2009. While we do not find a pronounced change in ozone-temperature sensitivity patterns for the individual site categories, we find a decrease in the overall magnitude of the ozone-temperature sensitivity in spring and an increase in summer. This is particularly pronounced in spring, where absolute values decreased from ranging between 1.9 and 2.8 μg °C−1 m−3 in 1990–1999 to 0.8 and 1.4 μg °C−1 m−3 in 2010–2019. Ozone production is predominantly VOC limited in spring. Thus, the reduction of anthropogenic NMVOC emissions propels less efficient ozone production regimes, which reflects also in the decrease of the overall temperature sensitivity of chemical ozone production.

As in spring, temperature sensitivity patterns do not change during summer over time except for urban sites in 2000–2009, which can be related to altered NOx to VOC ratios during that time (i.e., more significant reduction in NMVOCs compared to NOx). In absolute terms, the temperature sensitivity increased over time during summer by about 10% at urban sites.

5 Discussion and conclusions

In this study, we present an analysis of 30 years of ozone observations available from rural, suburban, and urban background sites in Austria. Our analysis focuses on days with MDA8 O3 exceeding the national and the European Union's target threshold for the protection of human health. To reduce the noise stemming from meteorological variability in the analysis, we assess changes over three 10 year time periods covering the years 1990–1999, 2000–2009 and 2010–2019. MDA8 O3 concentrations larger than 120 μg m−3 are still frequently observed at all site categories (rural, suburban, and urban) during spring and summer despite substantial precursor emission declines: NOx emissions were cut by 18% (observed concentrations declined by ∼30–40%) and anthropogenic NMVOC emissions by 53% in 2010–2019 compared to 1990–1999.

We investigated several factors affecting MDA8 O3 to understand the changes occurring in the past and assess possible future developments. We studied the modified VOC to NOx ratios emerging from diverse reductions of anthropogenic precursors over the 30 year study time-period. Altered precursor ratios explain the different impacts of emission abatements at various site categories. Our results show that during summer (JJA), ozone production is predominantly NOx limited at rural sites and VOC limited at urban sites. Furthermore, for suburban environments, we identify over time a shift in the ozone chemical regime towards NOx limitation. Focusing on the MDA8 O3 distribution, the most significant reduction of the upper tail (75% quantile, q75) MDA8 O3 is found comparing 2010–2019 to 1990–1999 at (NOx limited) rural sites (q75 minus 8 μg m−3vs. suburban (urban) sites: q75 minus 4 (4) μg m−3). Similarly, the largest reduction of exceedance days is found at rural sites (minus 7 exceedance days), while slight increases are found at suburban (and urban) sites (plus 1 (2) exceedance days). The narrowing of the MDA8 O3 probability distribution function is caused by lower anthropogenic NOx and VOC emissions. As a result, MDA8 O3 maxima decreased while minima increased simultaneously. The extent of the site-level changes depends on the prevailing ozone chemistry. In the VOC limited regime, the MDA8 O3 median increased between the initial and the final time period by 4 μg m−3 (initially at 94 (96) μg m−3 for suburban (urban) sites). At NOx limited rural sites, the median remained unchanged at 104 μg m−3, and overall rural sites continued to show the highest ozone concentrations across all site categories.

Concerning the temperature sensitivity of ozone production, our study shows that the temperature sensitivity of MDA8 O3 is one of the main drivers of elevated ozone levels. Increasing temperatures have been boosting ozone production to such an extent that in recent hot summers (2015, 2018, and 2019) the number of annual exceedance days (MDA8 O3 > 120 μg m−3) was well above the target of 25 days at all site categories. While we do not find a pronounced change in ozone-temperature sensitivity patterns for the individual site categories, we see a decrease in the overall magnitude of the ozone-temperature sensitivity in spring and an increase in summer. The analysis shows further that the temperature increase annihilated the benefits from emission reductions. The latter is clearly seen when we repeat our analysis and exclude hot summers. For MDA8 O3, we find a decrease in the median, peak values (q75), and exceedance days at all site categories (rural: median −5 μg m−3, q75 −16 μg m−3, exceedance days −13; suburban/urban: median −2 μg m−3, q75 −14 μg m−3, exceedance days −5) in 2010–2019 compared to 1990–1999.

Given the strong temperature sensitivity of MDA8 O3, we analyse MDA8 O3 under clear sky conditions for NOx and VOC limited regimes separately. We do so to quantify the effect concerning changes in DA NOx mixing ratio and Tmean. Our results show that the temperature sensitivity of ozone formation increases wtih rising Tmean. This pattern is found in spring and in summer. Still, the absolute temperature sensitivity is, on average, a factor of 2.5 larger in the summertime than in spring at all site categories (4.4–5.6 μg per m3 per °C Tmean in JJA). Besides enhanced chemical kinetics at elevated temperatures, increased BVOC emissions, reduced dry deposition rates, and faster PAN decomposition (NO2 source, specifically relevant in the NOx limited regime) result in increased ozone concentrations. The enhancement of MDA8 O3 due to higher ambient temperatures is in the same order as the ozone reduction achieved due to cuts in anthropogenic precursor emissions over the last 30 years.

For the city of Vienna, we employ MAX-DOAS HCHO measurements (06/2017–12/2019) to evaluate whether BVOCs are the dominant source in the metropolitan area in spring and summer. Throughout the growing season (April–October), we find a distinct positive relationship between HCHO and temperature, indicative of the hydrocarbons' biogenic source. Additionally, we derive ozone isopleths by combining the MAX-DOAS HCHO measurements and in situ O3 and NOx observations (Vienna, Stephansplatz). The isopleths confirm that peak MDA8 O3 is VOC limited.

In conclusion, our analysis shows that during the last three decades, rising air temperatures have been counteracting the effects of ozone precursor reductions and altering the NOx to VOC ratios by modifying the contribution of BVOCs to the total VOC burden. In particular, the latter poses a challenge for air quality managers, given the temperature dependence of both BVOC emission and ozone production, in the light of progressive climate change. The assessment of potential changes in the future Austrian ozone burden – and the relative contribution of changes in temperature, and anthropogenic and biogenic precursor emissions to it – is left for future modelling work, following the representative concentration pathways.53

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

M. Mayer, C. Staehle and H. E. Rieder acknowledge support by the Austrian Climate and Energy Fund via project ACRP11-KR18AC0K14686. The authors thank the VINDOBONA team for helping to establish a MAX-DOAS measurement network in Vienna. The BOKU MAX-DOAS instrument as part of the VINDOBONA project was funded by the German Research Foundation (grant no. Ri 1800/6-1) and Stefan F. Schreier received funding from the Austrian Science Fund (grant no. I 2296-N29).

References

  1. M. L. Bell, F. Dominici and J. M. Samet, A meta-analysis of time-series studies of ozone and mortality with comparison to the national morbidity, mortality, and air pollution study, Epidemiology, 2005, 16, 436–445 CrossRef PubMed.
  2. M. L. Bell, A. McDermott, S. L. Zeger, J. M. Samet and F. Dominici, Ozone and short-term mortality in 95 US urban communities, 1987-2000, JAMA, J. Am. Med. Assoc., 2004, 292, 2372–2378 CrossRef CAS PubMed.
  3. U. Confalonieri, B. Menne, R. Akhtar, K. L. Ebi, M. Hauengue, R. S. Kovats, B. Revich and A. Woodward, Human Health. Climate Change 2007: Impacts, Adaptation and Vulnerability, Cambridge, 2007 Search PubMed.
  4. EEA, Air Quality in Europe — 2020 Report, EEA Report No 09/2020, 2020,  DOI:10.2800/786656.
  5. A. Arneth, S. P. Harrison, S. Zaehle, K. Tsigaridis, S. Menon, P. J. Bartlein, J. Feichter, A. Korhola, M. Kulmala, D. O'Donnell, G. Schurgers, S. Sorvari and T. Vesala, Terrestrial biogeochemical feedbacks in the climate system, Nat. Geosci., 2010, 3, 525–532 CrossRef CAS.
  6. IPCC, in Fourth Assessment Report of the Intergovernmental Panel on Climate Change 2007, pp. 391–431 Search PubMed.
  7. D. L. Mauzerall and X. P. Wang, Protecting agricultural crops from the effects of tropospheric ozone exposure: Reconciling science and standard setting in the United States, Europe, and Asia, Annual Review of Energy and the Environment, 2001, 26, 237–268 CrossRef.
  8. P. Kumar and B. Imam, Footprints of air pollution and changing environment on the sustainability of built infrastructure, Sci. Total Environ., 2013, 444, 85–101 CrossRef CAS PubMed.
  9. IPCC, in Intergovernmental Panel on Climate Change Climate Change: the Physical Science Basis,Contribution of Working Group I to the Fourth Assessment Report of the IPCC, Cambridge University Press, 2007, pp. 79–131 Search PubMed.
  10. D. Shindell, J. C. I. Kuylenstierna, E. Vignati, R. van Dingenen, M. Amann, Z. Klimont, S. C. Anenberg, N. Muller, G. Janssens-Maenhout, F. Raes, J. Schwartz, G. Faluvegi, L. Pozzoli, K. Kupiainen, L. Hoglund-Isaksson, L. Emberson, D. Streets, V. Ramanathan, K. Hicks, N. T. K. Oanh, G. Milly, M. Williams, V. Demkine and D. Fowler, Simultaneously Mitigating Near-Term Climate Change and Improving Human Health and Food Security, Science, 2012, 335, 183–189 CrossRef CAS.
  11. D. S. Stevenson, P. J. Young, V. Naik, J. F. Lamarque, D. T. Shindell, A. Voulgarakis, R. B. Skeie, S. B. Dalsoren, G. Myhre, T. K. Berntsen, G. A. Folberth, S. T. Rumbold, W. J. Collins, I. A. MacKenzie, R. M. Doherty, G. Zeng, T. P. C. van Noije, A. Strunk, D. Bergmann, P. Cameron-Smith, D. A. Plummer, S. A. Strode, L. Horowitz, Y. H. Lee, S. Szopa, K. Sudo, T. Nagashima, B. Josse, I. Cionni, M. Righi, V. Eyring, A. Conley, K. W. Bowman, O. Wild and A. Archibald, Tropospheric ozone changes, radiative forcing and attribution to emissions in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 2013, 13, 3063–3085 CrossRef.
  12. V. Naik., S. Szopa, B. Adhikary, P. Artaxo, T. Berntsen, W. D. Collins, S. Fuzzi, L. Gallardo, A. Kiendler Scharr, Z. Klimont, H. Liao, N. Unger and P. Zanis, Short-lived climate forcers, in Climate Change 2021: the Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, 2021 Search PubMed.
  13. R. B. Skeie, G. Myhre, O. Hodnebrog, P. J. Cameron-Smith, M. Deushi, M. I. Hegglin, L. W. Horowitz, R. J. Kramer, M. Michou, M. J. Mills, D. J. L. Olivie, F. M. O. Connor, D. Paynter, B. H. Samset, A. Sellar, D. Shindell, T. Takemura, S. Tilmes and T. W. Wu, Historical total ozone radiative forcing derived from CMIP6 simulations, npj Clim. Atmos. Sci., 2020, 3, 1–10 CrossRef.
  14. P. J. Young, A. T. Archibald, K. W. Bowman, J. F. Lamarque, V. Naik, D. S. Stevenson, S. Tilmes, A. Voulgarakis, O. Wild, D. Bergmann, P. Cameron-Smith, I. Cionni, W. J. Collins, S. B. Dalsoren, R. M. Doherty, V. Eyring, G. Faluvegi, L. W. Horowitz, B. Josse, Y. H. Lee, I. A. MacKenzie, T. Nagashima, D. A. Plummer, M. Righi, S. T. Rumbold, R. B. Skeie, D. T. Shindell, S. A. Strode, K. Sudo, S. Szopa and G. Zeng, Pre-industrial to end 21st century projections of tropospheric ozone from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 2013, 13, 2063–2090 CrossRef.
  15. P. S. Monks, A. T. Archibald, A. Colette, O. Cooper, M. Coyle, R. Derwent, D. Fowler, C. Granier, K. S. Law, G. E. Mills, D. S. Stevenson, O. Tarasova, V. Thouret, E. von Schneidemesser, R. Sommariva, O. Wild and M. L. Williams, Tropospheric ozone and its precursors from the urban to the global scale from air quality to short-lived climate forcer, Atmos. Chem. Phys., 2015, 15, 8889–8973 CrossRef CAS.
  16. X. Lu, L. Zhang and L. Shen, Meteorology and Climate Influences on Tropospheric Ozone: a Review of Natural Sources, Chemistry, and Transport Patterns, Curr. Pollut. Rep., 2019, 238–260 CrossRef.
  17. X. Wang, Z. X. Shen, J. J. Cao, L. M. Zhang, L. Liu, J. J. Li, S. X. Liu and Y. F. Sun, Characteristics of surface ozone at an urban site of Xi'an in Northwest China, J. Environ. Monit., 2012, 14, 116–126 RSC.
  18. P. Cristofanelli, P. Bonasoni, L. Tositti, U. Bonafe, F. Calzolari, F. Evangelisti, S. Sandrini and A. Stohl, A 6-year analysis of stratospheric intrusions and their influence on ozone at Mt. Cimone (2165 m above sea level), J. Geophys. Res.: Atmos., 2006, 111, 1–11 CrossRef.
  19. M. Lin, A. M. Fiore, L. W. Horowitz, A. O. Langford, S. J. Oltmans, D. Tarasick and H. E. Rieder, Climate variability modulates western US ozone air quality in spring via deep stratospheric intrusions, Nat. Commun., 2015, 6, 7105 CrossRef CAS PubMed.
  20. Z. L. Fleming, R. M. Doherty, E. von Schneidemesser, C. S. Malley, O. R. Cooper, J. P. Pinto, A. Colette, X. B. Xu, D. Simpson, M. G. Schultz, A. S. Lefohn, S. Hamad, R. Moolla, S. Solberg and Z. Z. Feng, Tropospheric Ozone Assessment Report: Present-day ozone distribution and trends relevant to human health, Elementa: Science of the Anthropocene, 2018, 6, 12 Search PubMed.
  21. EUR-Lex, EU: Directive 2008/50/EC of the European Parliament and of the Council of 21 May 2008 on ambient air quality and cleaner air for Europe, Official Journal of the European Union, L 152, 11 June 2008, 2008, 51, 1–44, https://data.europa.eu/eli/dir/2008/50/oj Search PubMed.
  22. Umweltbundesamt, Austria's Annual Air Emissions Inventory 1990–2019. Emissions of SO2, NOx, NMVOC, NH3 and PM2.5, with assistance of M. Anderl, S. Haider, C. Lampert, D. Perl, M. Pinterits, S. Poupa, M. Purzner, W. Schieder, G. Schmidt, B. Schodl, M. Titz and M. Wieser, Reports, Bd. REP-0760, 2021 Search PubMed.
  23. Umweltbundesamt, Jahresbericht der Luftgütemessungen in Österreich 2019, W. Spangl and C. Nagl, 2020 Search PubMed.
  24. B. Chimani, G. Heinrich, M. Hofstätter, M. Kerschbaumer, S. Kienberger, A. Leuprecht, A. Lexer, S. Peßenteiner, M. S. Poetsch, M. Salzmann, R. Spiekermann, M. Switanek and H. Truhetz, ÖKS15 – Klimaszenarien für Österreich. Daten, Methoden und Klimaanalyse, Vienna, Austria, 2016 Search PubMed.
  25. IPCC, Climate Change 2014: Impacts, Adaptation, and Vulnerability – Part B: Regional Aspects. Contribution of Working Group II to the Fifth AssessmentReport of the Intergovernmental Panel on Climate Change, Cambridge, United Kingdom and New York, NY, USA, 2014 Search PubMed.
  26. EUR-Lex, 2011/850/EU: Commission Implementing Decision of 12 December 2011 laying down rules for Directives 2004/107/EC and 2008/50/EC of the European Parliament and of the Council as regards the reciprocal exchange of information and reporting on ambient air quality (notified under document C(2011) 9068), Official Journal of the European Union, 2011, L335 (OJ L 335), 54, 86–106, https://eur-lex.europa.eu/eli/dec_impl/2011/850/oj Search PubMed.
  27. B. Brown-Steiner, N. E. Selin, R. G. Prinn, E. Monier, S. Tilmes, L. Emmons and F. Garcia-Menendez, Maximizing ozone signals among chemical, meteorological, and climatological variability, Atmos. Chem. Phys., 2018, 18, 8373–8388 CrossRef CAS.
  28. A. C. Davison, Encyclopedia of Biostatistics: 4 Extreme Values, John Wiley & Sons, Ltd., USA, 2005 Search PubMed.
  29. H. E. Rieder, A. M. Fiore, L. M. Polvani, J. F. Lamarque and Y. Fang, Changes in the frequency and return level of high ozone pollution events over the eastern United States following emission controls, Environ. Res. Lett., 2013, 8, 014012 CrossRef CAS.
  30. S. Sillman and D. Y. He, Some theoretical results concerning O-3-NOx-VOC chemistry and NOx-VOC indicators, J. Geophys. Res.: Atmos., 2002, 107, ACH 26-1–ACH 26-15 CrossRef.
  31. B. Vogel, N. Riemer, H. Vogel and F. Fiedler, Findings on NOy as an indicator for ozone sensitivity based on different numerical simulations, J. Geophys. Res.: Atmos., 1999, 104, 3605–3620 CrossRef CAS.
  32. M. Bauwens, T. Stavrakou, J. F. Müller, I. De Smedt, M. Van Roozendael, G. R. van der Werf, C. Wiedinmyer, J. W. Kaiser, K. Sindelarova and A. Guenther, Nine years of global hydrocarbon emissions based on source inversion of OMI formaldehyde observations, Atmos. Chem. Phys., 2016, 16, 10133–10158 CrossRef CAS.
  33. G. Curci, P. I. Palmer, T. P. Kurosu, K. Chance and G. Visconti, Estimating European volatile organic compound emissions using satellite observations of formaldehyde from the Ozone Monitoring Instrument, Atmos. Chem. Phys., 2010, 10, 11501–11517 CrossRef CAS.
  34. H. Singh, Y. Chen, A. Staudt, D. Jacob, D. Blake, B. Heikes and J. Snow, Evidence from the Pacific troposphere for large global sources of oxygenated organic compounds, Nature, 2001, 410, 1078–1081 CrossRef CAS.
  35. B. Franco, E. A. Marais, B. Bovy, W. Bader, B. Lejeune, G. Roland, C. Servais and E. Mahieu, Diurnal cycle and multi-decadal trend of formaldehyde in the remote atmosphere near 46A degrees aEuro N, Atmos. Chem. Phys., 2016, 16, 4171–4189 CrossRef CAS.
  36. R. Atkinson and J. Arey, Gas-phase tropospheric chemistry of biogenic volatile organic compounds: a review, Atmos. Environ., 2003, 37, S197–S219 CrossRef CAS.
  37. A. Calogirou, B. R. Larsen and D. Kotzias, Gas-phase terpene oxidation products: a review, Atmos. Environ., 1999, 33, 1423–1439 CrossRef CAS.
  38. Y. N. Lee, X. Zhou, L. I. Kleinman, L. J. Nunnermacker, S. R. Springston, P. H. Daum, L. Newman, W. G. Keigley, M. W. Holdren, C. W. Spicer, V. Young, B. Fu, D. D. Parrish, J. Holloway, J. Williams, J. M. Roberts, T. B. Ryerson and F. C. Fehsenfeld, Atmospheric chemistry and distribution of formaldehyde and several multioxygenated carbonyl compounds during the 1995 Nashville Middle Tennessee Ozone Study, J. Geophys. Res.: Atmos., 1998, 103, 22449–22462 CrossRef CAS.
  39. D. D. Parrish, T. B. Ryerson, J. Mellqvist, J. Johansson, A. Fried, D. Richter, J. G. Walega, R. A. Washenfelder, J. A. de Gouw, J. Peischl, K. C. Aikin, S. A. McKeen, G. J. Frost, F. C. Fehsenfeld and S. C. Herndon, Primary and secondary sources of formaldehyde in urban atmospheres: Houston Texas region, Atmos. Chem. Phys., 2012, 12, 3273–3288 CrossRef CAS.
  40. W. Su, C. Liu, Q. Hu, S. Zhao, Y. Sun, W. Wang, Y. Zhu, J. Liu and J. Kim, Primary and secondary sources of ambient formaldehyde in the Yangtze River Delta based on Ozone Mapping and Profiler Suite (OMPS) observations, Atmos. Chem. Phys., 2019, 19, 6717–6736 CrossRef CAS.
  41. S. F. Schreier, A. Richter, E. Peters, M. Ostendorf, A. W. Schmalwieser, P. Weihs and J. P. Burrows, Dual ground-based MAX-DOAS observations in Vienna, Austria: Evaluation of horizontal and temporal NO2, HCHO, and CHOCHO distributions and comparison with independent data sets, Atmos. Environ.: X, 2020, 5, 100059 CAS.
  42. EEA, Emissions of the main air pollutants in Europe, https://www.eea.europa.eu/data-and-maps/indicators/main-anthropogenic-air-pollutant-emissions/assessment-6 Search PubMed.
  43. M. Bauwens, T. Stavrakou, J. F. Muller, B. Van Schaeybroeck, L. De Cruz, R. De Troch, O. Giot, R. Hamdi, P. Termonia, Q. Laffineur, C. Amelynck, N. Schoon, B. Heinesch, T. Holst, A. Arneth, R. Ceulemans, A. Sanchez-Lorenzo and A. Guenther, Recent past (1979-2014) and future (2070-2099) isoprene fluxes over Europe simulated with the MEGAN-MOHYCAN model, Biogeosciences, 2018, 15, 3673–3690 CrossRef CAS.
  44. D. Simpson and W. Winiwarter, Emissions from Natural Sources, Environment Agency Austria, Vienna, Austria, 1998 Search PubMed.
  45. C. A. Skjoth, C. Geels, M. Hvidberg, O. Hertel, J. Brandt, L. M. Frohn, K. M. Hansen, G. B. Hedegard, J. H. Christensen and L. Moseholm, An inventory of tree species in Europe - An essential data input for air pollution modelling, Ecol. Modell., 2008, 217, 292–304 CrossRef CAS.
  46. D. C. Oderbolz, S. Aksoyoglu, J. Keller, I. Barmpadimos, R. Steinbrecher, C. A. Skjoth, C. Plass-Dulmer and A. S. H. Prevot, A comprehensive emission inventory of biogenic volatile organic compounds in Europe: improved seasonality and land-cover, Atmos. Chem. Phys., 2013, 13, 1689–1712 CrossRef.
  47. G. Curci, M. Beekmann, R. Vautard, G. Smiatek, R. Steinbrecher, J. Theloke and R. Friedrich, Modelling study of the impact of isoprene and terpene biogenic emissions on European ozone levels, Atmos. Environ., 2009, 43, 1444–1455 CrossRef CAS.
  48. M. Karl, A. Guenther, R. Koble, A. Leip and G. Seufert, A new European plant-specific emission inventory of biogenic volatile organic compounds for use in atmospheric transport models, Biogeosciences, 2009, 6, 1059–1087 CrossRef CAS.
  49. R. Steinbrecher, G. Smiatek, R. Koble, G. Seufert, J. Theloke, K. Hauff, P. Ciccioli, R. Vautard and G. Curci, Intra- and inter-annual variability of VOC emissions from natural and semi-natural vegetation in Europe and neighbouring countries, Atmos. Environ., 2009, 43, 1380–1391 CrossRef CAS.
  50. K. Riahi, S. Rao, V. Krey, C. Cho, V. Chirkov, G. Fischer, G. Kindermann, N. Nakicenovic and P. Rafaj, RCP 8.5—A scenario of comparatively high greenhouse gas emissions, Clim. Change, 2011, 109, 33 CrossRef CAS.
  51. A. Guenther, T. Karl, P. Harley, C. Wiedinmyer, P. I. Palmer and C. Geron, Estimates of global terrestrial isoprene emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature), Atmos. Chem. Phys., 2006, 6, 3181–3210 CrossRef CAS.
  52. J. Peñuelas and M. Staudt, BVOCs and global change, Trends Plant Sci., 2010, 15, 133–144 CrossRef.
  53. D. P. van Vuuren, J. Edmonds, M. Kainuma, K. Riahi, A. Thomson, K. Hibbard, G. C. Hurtt, T. Kram, V. Krey, J.-F. Lamarque, T. Masui, M. Meinshausen, N. Nakicenovic, S. J. Smith and S. K. Rose, The representative concentration pathways: an overview, Clim. Change, 2011, 109, 5 CrossRef.
  54. Umweltbundesamt, Austria's National Air EmissionPollutants: NOx, SO2, NMVOC, NH3 and PM2.5 : Projections 2021 for 2020, 2025 and 2030, 2021 Search PubMed.
  55. J. Hiebl and C. Frei, Daily temperature grids for Austria since 1961-concept, creation and applicability, Theor. Appl. Climatol., 2016, 124, 161–178 CrossRef.
  56. J. Hiebl and C. Frei, Daily precipitation grids for Austria since 1961-development and evaluation of a spatial dataset for hydroclimatic monitoring and modelling, Theor. Appl. Climatol., 2018, 132, 327–345 CrossRef.
  57. M. Olefs and W. Schöner, presented in part at the EGU General Assembly Vienna, Austria, 22.04–27.04.2012, 2012 Search PubMed.
  58. K. L. Chang, I. Petropavlovskikh, O. R. Cooper, M. G. Schultz and T. Wang, Regional trend analysis of surface ozone observations from monitoring networks in eastern North America, Europe and East Asia, Elementa: Science of the Anthropocene, 2017, 5, 50 Search PubMed.
  59. M. Stangl, H. Formayer, A. Höfler, K. Andre, M. Kalcher, J. Hiebl, M. Hofstätter, A. Orlik and C. Michl, Klimastatusbericht Österreich 2019, 2020 Search PubMed.
  60. V. Ferracci, C. G. Bolas, R. A. Freshwater, Z. Staniaszek, T. King, K. Jaars, F. Otu-Larbi, J. Beale, Y. Malhi, T. W. Waine, R. L. Jones, K. Ashworth and N. R. P. Harris, Continuous Isoprene Measurements in a UK Temperate Forest for a Whole Growing Season: Effects of Drought Stress During the 2018 Heatwave, Geophys. Res. Lett., 2020, 47, e2020GL088885 CrossRef CAS.
  61. M. Y. Lin, L. W. Horowitz, Y. Y. Xie, F. Paulot, S. Malyshev, E. Shevliakova, A. Finco, G. Gerosa, D. Kubistin and K. Pilegaard, Vegetation feedbacks during drought exacerbate ozone air pollution extremes in Europe, Nat. Clim. Change, 2020, 10, 444–451 CrossRef CAS.
  62. T. Zenone, C. Hendriks, F. Brilli, E. Fransen, B. Gioli, M. Portillo-Estrada, M. Schaap and R. Ceulemans, Interaction between isoprene and ozone fluxes in a poplar plantation and its impact on air quality at the European level, Sci. Rep., 2016, 6 CAS.
  63. J. Coates, K. A. Mar, N. Ojha and T. M. Butler, The influence of temperature on ozone production under varying NOx conditions - a modelling study, Atmos. Chem. Phys., 2016, 16, 11601–11615 CrossRef CAS.

Footnote

Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d2ea00004k

This journal is © The Royal Society of Chemistry 2022