Articles | Volume 21, issue 4
Research article
02 Mar 2021
Research article |  | 02 Mar 2021

Decoupling of urban CO2 and air pollutant emission reductions during the European SARS-CoV-2 lockdown

Christian Lamprecht, Martin Graus, Marcus Striednig, Michael Stichaner, and Thomas Karl

Lockdown and the associated massive reduction in people's mobility imposed by SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) mitigation measures across the globe provide a unique sensitivity experiment to investigate impacts on carbon and air pollution emissions. We present an integrated observational analysis based on long-term in situ multispecies eddy flux measurements, allowing for quantifying near-real-time changes of urban surface emissions for key air quality and climate tracers. During the first European SARS-CoV-2 wave we find that the emission reduction of classic air pollutants decoupled from CO2 and was significantly larger. These differences can only be rationalized by the different nature of urban combustion sources and point towards a systematic bias of extrapolated urban NOx emissions in state-of-the-art emission models. The analysis suggests that European policies, shifting residential, public, and commercial energy demand towards cleaner combustion, have helped to improve air quality more than expected and that the urban NOx flux remains to be dominated (e.g., >90 %) by traffic.

1 Introduction

Managing air pollution and climate change are among the most important environmental challenges of modern society. As urban population continues to grow, emissions from metropolitan areas play an increasingly important role. For example, European cities already host about 74 % of the population (UN, 2019) and are a major contributor to air pollutant and greenhouse gas emissions. Urban growth, along with socioeconomic development and without mitigation, can lead to substantial increases in anthropogenic emissions. Many cities are committing to sustainable development goals, and improvement of air pollution and mitigation of climate change are emerging as key sustainability priorities across the globe. Quantifying the diversity of urban emissions is often one of the most uncertain components of complex atmospheric models, and development of a robust predictive capability requires accurate data and careful evaluation of bottom-up emissions (Blain et al., 2019; NAS, 2016).

During the last 2 decades Europe's policy to reduce mid-term carbon emissions has fostered the proliferation of diesel-driven vehicles (EU-EUR-Lex, 2008). While soot emissions can be successfully removed with a diesel particulate filter, the reduction of NOx from diesel exhaust has been more challenging and was at the center of “Dieselgate” (Franco et al., 2014). As a consequence, European NOx concentrations have declined less rapidly than elsewhere (Carslaw and Rhys-Tyler, 2013; Im et al., 2015; Karl et al., 2017) and put the EU-28 (28 member states of the European Union) emission target for NOx reductions (2005–2030: −63 %) in jeopardy (EU-EUR-Lex, 2008). Nitrogen oxides have therefore emerged as a primary public health concern (Anenberg et al., 2017). European suppression measures due to the SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) outbreak provide a unique opportunity to track drastic changes in urban mobility during the lockdown phase and combined with eddy flux methods allow for investigating the sensitivity towards emission changes directly.

After the initial SARS-CoV-2 emergence in China in late 2019, the World Health Organization declared the outbreak of a global pandemic on 11 March 2020. Worldwide measures to mitigate or suppress exponential growth of SARS-CoV-2 have resulted in an unprecedented global intervention on mobility and industrial activity (WHO, 2020), allowing for studying a number of environmental aspects (e.g., F. Liu et al., 2020; Schiermeier, 2020; Le Quéré et al., 2020). A growing number of studies document changes of regional and global air composition (e.g., Menut et al., 2020; Bao and Zhang, 2020) with respect to lockdown measures, including remote sensing observations and aspects of adequate data-processing strategies (e.g., F. Liu et al., 2020; Sussmann and Rettinger, 2020).

In Europe, most countries have implemented suppression strategies involving a more or less extensive lockdown of public life. At the beginning of the pandemic, the level of suppression varied among different countries, with some imposing very early (“China-like”) lockdown measures (e.g., Austria), while others shifted from gradual social-distancing measures to a lockdown after reconsideration of alternative strategies (e.g., the UK). Depending on the magnitude of the outbreak, European countries put increasingly stringent measures in place. The extent of different lockdown measures has been assessed early on via cell phone activity tracking. For example, Google mobility reports published in March 2020 suggested an 80 % reduction of retail and recreational activities across Europe. Traffic count data show a 60 % reduction of urban mobility due to a state-wide quarantine in the state of Tyrol early during the pandemic. Such a drastic mobility reduction during the suppression period allows for performing a granular assessment of processes impacting emissions and the distribution of air pollutant and climate gases.

A direct and quantitative way to assess air pollutant and climate gas emission changes can be based on the eddy covariance method (Aubinet et al., 2012; Dabberdt et al., 1993). Briefly, in its simplest form for stationary conditions and neglecting horizontal advection, the turbulent surface–atmosphere flux (measured at height h) can represent the diffusive flux at the surface:

(1) w c h = - D c z 0 ,

where w represents the vertical fluctuation of wind speed, D is the molecular diffusion coefficient, and c is the concentration fluctuation. The turbulent flux at measurement height h (left side) equals the diffusive surface flux (right side), which we are usually interested in. Brackets denote the averaging interval. The ensemble average is typically 30 min. Eddy covariance measurements have been extensively used in atmospheric sciences (Foken and Wichura, 1996; Oncley et al., 2007; Patton et al., 2011) and biogeochemistry (Aubinet et al., 2012; Baldocchi et al., 1988; Fowler et al., 2009; Rannik et al., 2012) (e.g., AmeriFlux:, last access: 18 December 2020; Euroflux:, last access: 18 December 2020). The method has also become more tractable for reactive trace gases such as NMVOCs (non-methane volatile organic compounds; Karl et al., 2001; Rinne et al., 2001; Spirig et al., 2005) or NOx (Lee et al., 2015) and has been used at urban sites (Christen, 2014; Langford et al., 2009; Velasco et al., 2005; Squires et al., 2020). Urban eddy covariance methods can monitor aggregated emission changes in real time. Here we build on a set of long-term multispecies flux and concentration datasets for NOx, O3, aromatic NMVOCs, and CO2 (Karl et al., 2020). Being inspired by early empirical persistence models used in atmospheric chemistry and ecology, we propose a new quantitative way for the analysis of urban fluxes during an intervention experiment by combining eddy covariance data with a boosted regression tree model (Duffy and Helmbold, 2002). This method allows for directly assessing changes of surface fluxes for different trace gases in response to the SARS-CoV-2 lockdown and rebounding effects.

2 Methods

2.1 Eddy covariance

Here we analyze air quality data based on the eddy covariance method (e.g., Aubinet et al., 2012; Dabberdt et al., 1993), which represents the most direct meteorological method to determine surface emissions (Baldocchi et al., 1988; Fowler et al., 2009). The method is widely established in biogeosciences (e.g., Euroflux,, last access: 30 January 2021; AmeriFlux,, last access: 30 January 2021). A number of studies investigated eddy fluxes of chemical species and aerosols in urban settings (Nemitz et al., 2008; Velasco et al., 2009; Rantala et al., 2016; Lee et al., 2015; Karl et al., 2017; Vaughan et al., 2017; Striednig et al., 2020). Briefly, the method relies on the conservation equation of a scalar, which under homogenous conditions and can be simplified to

(2) C t + F z = S ,

where dC/dt represents the storage term, dF/dz is the measured vertical turbulent flux as in Eq. (1), and S is sources and sinks between the surface and height z.

Integration of Eq. (2) yields

(3) 0 h C t d z + w c h = w c 0 := F s ,

where h is the measurement height (39 m above street level) and Fs represents the surface flux. In this context the turbulent flux term usually dominates the left-hand side. The storage term typically accounted for 5 %–7 % of the fluxes on average, and we consider it a minor component in our analysis. Similarly we neglect advection fluxes. We take advantage of the fact that our analysis is based on relative changes of different air pollutant fluxes normalized by a boosted regression tree model. Any systematic bias therefore cancels out under the assumption of comparable source distributions verified by the emission inventory. The datasets are analyzed with the innFLUX code (Striednig et al., 2020), which outputs standard parameters used for filtering flux data as described by Foken and Wichura (1996). In addition to raw data filters (spike removals and weather flags), we applied standard criteria using u* and the stationarity criterion for all species (Foken and Wichura, 1996). We specifically do not apply tests on integral turbulence, as parameterizations for urban areas are not available or accurate. After applying the abovementioned filters, 73 % of the flux data were used for the training dataset, and 82 % of the flux data were used for the intervention period. Systematic errors due to attenuation caused by slow sensor response was assessed previously (Karl et al., 2017). It is considered minor due to the large eddy size above the urban roughness layer, and for the trace gases considered here it is on the order of 2 %–5 %. A detailed description of errors and data treatment for this site was published by Striednig et al. (2020).

2.2 Flux footprint and IAO observations

A site description of the Innsbruck Atmospheric Observatory (IAO), instrumentation, and site validation were previously extensively described (Karl et al., 2020). The flux footprint (Fig. 1) was calculated according to Kljun et al. (2015). For the measurement–inventory comparison we mapped the two-dimensional climatological footprint (March–May) onto the spatially disaggregated Austrian EMIKAT emission inventory (, last access: 30 January 2021). The relative seasonal variability was accounted for by scaling total yearly traffic emissions to measured seasonal traffic activity (state of Tyrol, Austria) and total yearly residential, commercial, and public (RCP) emissions to measured seasonal natural gas (NG) consumption (TIGAS, Austria,, last access: 30 January 2021). The land surface distribution within the flux footprint (Fig. S1 in the Supplement) is dominated by roads and building surfaces with a fraction between 70 % and 88 % depending on the wind sector. For comparisons the district level emission distribution from the inventory was mapped onto the land surface distribution and then weighted according to the footprint function. Traffic counts used in the data comparisons were based on a conductive loop measuring directional traffic flows along Innrain provided by the state of Tyrol, a main street dissecting the flux tower footprint, and are considered representative of traffic activity surrounding the flux tower. The inductive loop provides rudimentary information on light- vs. heavy-duty vehicles and suggests that 95 % of traffic is caused by vehicles <3.5 t. We assume that all fuel types used for heating appliances and warm-water consumption track relative changes of NG consumption, which is largely a function of base load and heating degree days (Fig. S2). Since many commercial buildings (e.g., shops, restaurants, and retail) are not clearly separable from residential buildings in European cities (e.g., upper floors are used for housing, and the ground floor houses shops or restaurants using shared heating), it is hard to clearly separate the RCP sector into individual components in the urban core. Overall, heating energy supply in the RCP sector is comprised of district heating (8.7 %), oil (34.5 %), natural gas (34 %), biomass (16 %), electricity (6.1 %), and alternative sources (0.7 %).

Figure 1Flux footprint surrounding the IAO tower plotted on top of a gridded land use map derived from OpenStreetMap (© OpenStreetMap contributors 2020. Distributed under a Creative Commons BY-SA License).

NOx measurements were based on a dual-channel chemiluminescence instrument (CLD 899 Y, ECO PHYSICS). The instrument was operated in flux mode, acquiring data at about 5 Hz. A NO standard was periodically introduced for calibration. Zeroing was performed once a day close to midnight. The chemiluminescent instrument is equipped with a metal oxide (i.e., molybdenum) converter. It has been shown (Steinbacher et al., 2007) that this can result in an overestimation of NO2 due to decomposition of NOy species. For Innsbruck we evaluated the accuracy using side-by-side measurements with a cavity ring-down spectrometer in 2015 (Karl et al., 2017). Both independent techniques agreed to within 6 %, confirming that this problem plays a minor role for polluted sites. CO2 and H2O were measured with a closed-path eddy covariance system (CPEC200, short inlet, enclosed IRGA – infrared gas analyzer – design, Campbell Scientific) along with three-dimensional winds. Calibration for CO2 was performed once a day. Aromatic NMVOCs (i.e., benzene, toluene, xylenes and ethylbenzene, and C9 benzenes) were measured with a PTR-TOF 6000 X2 mass spectrometer (IONICON, Austria), operated in hydronium mode at standard conditions in the drift tube of about 110 Td (Townsend). The instrument was set up to sample ambient air from a turbulently purged 3/8 in. Teflon line. Zero calibrations were performed by providing NMVOC-free air from a continuously purged catalytical converter though a setup of software-controlled solenoid valves. In addition, daily calibrations were performed using known quantities of a suite of NMVOCs from a 1 ppm calibration gas standard (Apel Riemer, USA) that were added to the NMVOC-free air and dynamically diluted into low-ppbv mixing ratios. Errors arising from analytical uncertainty mainly stem from calibration procedures. For NMVOCs these are estimated as 10 % for aromatic-NMVOC compounds based on a calibration standard; similarly the uncertainty of NOx is 2 %, and that of CO2 is 5 %, respectively.

This study builds on long-term NOx and CO2 flux measurements that have run operationally since 1 June 2018. NMVOC fluxes were measured during a field campaign from 11 March to 9 April 2019 and during the SARS-CoV-2 lockdown, when measurements started on 16 March 2020. The NMVOC analysis presented in this paper spans from 16 March to 1 May 2020.

2.3 Boosted regression tree model

Statistical persistence and regression models have a long history in atmospheric chemistry (Robeson and Steyn, 1990) for predicting empirical trends of pollutants (e.g., ozone) that factor in meteorological and chemical processes. These approaches have been used to forecast local surface ozone (Cobourn, 2007; Prybutok et al., 2000) and more recently trends of other atmospheric pollutants (Grange and Carslaw, 2019). Here we developed a boosted regression tree model using machine learning that is widely used in ecological modeling (Elith et al., 2008): for each variable we based the model on the following key astronomical and environmental driving variables: day of year (DOY), time of day (TOD), weekday or holiday (WDY), cartesian wind vectors (north–south and west–east directions), temperature (T), relative humidity (RH), global radiation (GR), and pressure (P). The model is set up using the machine learning toolbox in MATLAB (MathWorks Inc, USA) and trained for individual datasets until 29 February 2020 or during key reference periods (Table S1 in the Supplement). The model performance was assessed by comparing predicted and observed quantities using reference datasets (Table S2). To obtain a quantitative measure of emission changes, the differences between observed and predicted fluxes are integrated from the beginning of the lockdown period. As the predicted and observed quantities diverge, the integrated relative difference serves as a quantitative measure of emission (or activity) alteration (e.g., reduction).

2.4 Multispecies pollutant model

Based on two major and distinct urban pollution sources (i.e., road traffic and energy production in the residential, public, and commercial sectors) proportional contributions to the observed flux changes can be attributed based on a two-end-member mixing model: traffic emissions are primarily related to exhaust from internal combustion engines. The Austrian passenger car fleet is comprised of 43 % gasoline- and 55 % diesel-driven cars (Statistik Austria, 2019,, last access: 30 January 2021), with the latter being a key player for urban NOx emissions. The second significant emission source stems from fossil energy production in the residential, public, and commercial sectors with a significant contribution of natural gas combustion. In its simplest form we can therefore aggregate the observed flux changes into two main emission source categories using a two-end-member mixing model:

(4) δ F ( s ) F ( s ) = a s δ T T + b s δ R R + ε ,

where δF(s)F(s) is the measured relative flux difference between the output of the boosted regression tree model and the actual flux observations of species s (e.g., NOx, CO2, and aromatic NMVOCs), δTT is the traffic load difference determined from traffic count data, δRR is the residential energy consumption change, as and bs are proportionality terms, and ε is an error term. The proportionality terms (as and bs) represent the area weighted emission factors of the fleet average traffic (as) and RCP sector (bs). By definition as+bs:=1, if only two sources are considered.

3 Results

The urban NO–NO2–O3 triad. Due to the short atmospheric lifetime (e.g., up to 7 h, Laughner and Cohen, 2019) nitrogen oxides can serve as a gauge to assess air pollution changes as their atmospheric concentrations rapidly respond to shifting surface fluxes. The quantitative assessment of NOx emissions based on ambient-air concentrations however remains challenging due to non-linearities within the NONO2O3 triad in polluted regions (Lenschow et al., 2016). Under sunlight conditions and high NOx pollution the cycling between the NONO2O3 triad is described by the following reaction sequence:


The chemical timescale of the NOx triad (Eqs. 5 to 7) can be derived (Lenschow and Delany, 1987) as

(8) τ = 2 j 2 + k 3 2 O 3 - NO 2 + 2 j k 3 O 3 + NO + 2 NO 2 .

For typical conditions encountered during this study, this equates to timescales of about 100 s, comparable to the vertical turbulent exchange time in cities. Due to the rapid interconversion, the partitioning between NO and NO2 is typically dominated by fast chemistry, and the bulk of NO2 in the urban atmosphere is produced secondarily via the reaction of NO and O3. In the urban atmosphere this leads to a non-linear relationship between NO2 and NOx concentrations as depicted in Fig. 2. A repartitioning can be observed during the suppression phase for example, when the NO2-to-NOx trajectory moves from an urban NOx-saturated regime to a more NOx-limited regime. During the SARS-CoV-2 lockdown this shift was more pronounced than for typical weekend–weekday variations (Fig. 2b and c). As a consequence the relationship between changes of NOx fluxes and NO2 concentrations becomes a non-linear function of NOx concentrations when moving from NOx-saturated to NOx-limited conditions. Data from a nearby air quality station support these conclusions, showing significantly different NOx concentrations during the 2020 lockdown compared to the previous 5 years (i.e., a 50 % reduction of NOx) but no significant change for Ox (:=NO2+O3) based on the z hypothesis test. This chemical repartitioning and vertical redistribution in the surface layer needs to be accounted for when quantifying changing NOx emissions from concentrations. A more quantitative picture of changing NOx emissions can be obtained from direct flux measurements that are intrinsically linked to surface emissions (Vaughan et al., 2016).

Figure 2(a) Time series of mixing ratios of ambient NO2 and NOx before and during the lockdown. Gray-shaded vertical bars indicate weekends. The gray vertical solid line depicts the start of lockdown measures on 16 March 2020. (b) NO2 vs. NOx during weekdays (Tuesday to Thursday). (c) NO2 vs. NOx on Sundays.


Figure 3Diurnal variations of CO2 (a) and NOx (b) fluxes since 2018. Please note that the date format in this figure is the two-digit year followed by the month.


Figure 4Observed changes of air pollutant fluxes, CO2 flux, and traffic during the course of the first SARS-CoV-2 wave. (a) Normalized traffic counts, daily infection rate, and daily average flux reduction. (b) Cumulative reduction of NOx and CO2 fluxes and traffic activity. (c) Cumulative reduction of aromatic VOCs (AVOCs; volatile organic compounds), toluene (TOL), and benzene (BEN) fluxes. Red vertical lines indicate the start of university closure, Austrian lockdown (LO), school closure (SCL), and quarantine (QA) in the state of Tyrol. Green vertical lines show the lifting of mobility restrictions. Light-shaded areas represent the uncertainty of the boosted regression tree model analysis (Supplement). Please note that the date format in this figure is month day (mm/dd).


Figure 3 gives an overview of NOx and CO2 fluxes which have been continuously measured at the study site in central Europe since 2018. In addition, we have performed regular field campaigns augmenting these long-term datasets with NMVOC flux measurements (Karl et al., 2020). While atmospheric concentrations of primary air pollutants often exhibit strong surface maxima due to inversion layers during winter and spring, the corresponding emission fluxes typically track urban emission source activity and reflect changes in emission strengths and flux footprint. Turbulent fluxes typically exhibit midday maxima, reflecting increases in urban emission sources, which in the case of nitrogen oxides closely follow traffic load patterns (Karl et al., 2017). Urban CO2 fluxes follow these general trends but are to some extent less pronounced (e.g., weekend–weekday effect). During the vegetation period, CO2 emission fluxes can be suppressed (Ward et al., 2015) due to photosynthetic uptake by urban plants. For Innsbruck, we have assessed this effect previously and find that within the flux footprint the contribution of vegetation is relatively small (i.e., only about 10 % of the urban surface within the flux footprint is covered by plants). Urban CO2 fluxes are therefore primarily controlled by combustion processes. The flux site is situated in a valley with two dominant wind sectors, which cover a typical inner-city residential and business district (Fig. 1) with no significant industrial activities. In order to quantitatively investigate emission flux changes in response to SARS-CoV-2 intervention measures, we implemented a boosted regression tree model to define a business-as-usual scenario of the observed fluxes (Duffy and Helmbold, 2002). The model allows for factoring in differences in weather patterns (e.g., meteorological variations such as temperature, wind direction, and flux footprint) and describes changes that can be primarily attributed to the intervention itself. Accounting for seasonal differences is key to an accurate analysis of emission alterations due to lockdown measures. The essential time period of pre- and post-lockdown measures in Europe spans from March to about May 2020. Weather patterns in Europe can be particularly variable during this period, as the continent transitions from winter to summer. The climate of Tyrol is fairly representative of central Europe, where the transitional period between March to May can exhibit significant synoptic variability. For example, average monthly temperatures in March 2020 were about 0.9 K colder than in 2019. April and May 2020 tended to be 1.8 and 3.2 K warmer than 2019. Warmer temperatures in spring 2020 resulted in 24 % fewer heating degree days (HDDs) than in the year 2019 (Supplement). Consistent with these observations, natural gas consumption in Tyrol (Supplement) was reported to be 25 % lower during this period than in 2019. We can quantify changes of the observed fluxes due to the lockdown intervention in spring 2020 by referencing actual flux measurements to results from a trained boosted regression tree model (Fig. 4).

Shortly after the European SARS-CoV-2 outbreak first sparked in Italy, which was among the first European countries, the greater part of the Alps was under lockdown by mid-March to inhibit cross-border transmission. Tyrol implemented extensive measures of shelter in place and a state-wide quarantine (QA) on top of the Austrian lockdown (LO) on 16 March, 1 week after all universities closed. At the same time, Europe-wide measures of border control impacted all major north–south transport corridors to Italy. These measures resulted in massively reduced local mobility in combination with significant disruptions of one of the major transport routes across the Alps. As a consequence, average traffic loads in Innsbruck decreased by ∼64 %. The traffic data allow for partitioning traffic into “all vehicles”, “truck-similar vehicles”, “HDVs” (heavy-duty vehicles), and “semi-trailer trucks”. The reductions were 64 % (all vehicles), 40 % (truck-similar vehicles), 35 % (HDVs), and 21 % (semi-trailer trucks). Since it is an inner-city location, the fraction of passenger cars dominate. In absolute numbers, the distribution is dominated by passenger cars (<3.5 t), amounting to 95 % of all traffic, with the remainder attributed to the truck categories. The Austrian rate of infections reached a peak of 900 newly confirmed cases per day in mid-March and started to decline at the end of March. Along with efforts to reduce SARS-CoV-2 transmission, the shelter-in-place legislation resulted in a rapid decline of NOx, CO2, and aromatic NMVOCs (benzene, toluene, xylenes and ethylbenzene, and C9 benzenes) fluxes (Fig. 4a) reaching significantly lower emission fluxes relative to the business-as-usual scenario. The cumulative reduction of surface emissions of air pollutants (NOx and aromatic NMVOCs) closely follows traffic (Fig. 4b and c), declining by about 64 % during the lockdown period. At the end of the Austrian lockdown, traffic counts and integrated emissions of NOx and aromatic NMVOCs were −64 %, −59 %, and −56 % lower compared to the business-as-usual scenario. This reduction is significantly lower than observed for CO2 fluxes, leveling out at about −38 %. Notably benzene emission declines were also less pronounced than toluene and higher aromatic NMVOCs, which track NOx and traffic loads more closely. These different sensitivities indicate a non-linear relationship between the reduction of carbon dioxide and air pollution gases due to different urban combustion sources. Particularly reductions of NOx and CO2 exhibit quite different emission trajectories during the lockdown phase (Fig. 5). The observed reduction of air pollution gases, such as NOx, is significantly larger than estimated by early bottom-up model predictions for expected NOx-to-CO2 emission changes (Le Quéré et al., 2020). Can these observations be reconciled with bottom-up emission projections?

Figure 5Daily change of CO2 and NOx fluxes during the lockdown period. Flux observations are depicted by the blue dots. Emission model projections are represented by the solid orange line (Austrian emission inventory) and the dashed red line (Le Quéré et al., 2020).


4 Discussion

Our analysis indicates that the reduction of classic air pollutant emissions during the SARS-CoV-2 lockdown was more significant than that of CO2, which comes as a surprise. Comparable to most European countries, Austria-specific bottom-up emission models typically attribute 40 % of CO2 emissions to traffic and 19 % to the residential, commercial, and public (RCP) sector (UBA, 2019). For NOx, Austrian and European bottom-up emission projections predict similar contributions (i.e., 58 % from traffic and 12 % from the RCP sector). In its simplest form, by using a two-end-member pollutant model, we can test these assumptions in more detail and compare our observations with an Austrian state-of-the-art emission model (, last access: 30 January 2021) used for national emission reporting. For the analysis we take advantage of the fact that the seasonal influence on pollutant fluxes is factored in by referencing the flux analysis to the trained boosted regression tree model. Further, measured relative reductions of vehicle counts are assumed to represent the decrease of traffic activity reasonably well. We are then left with constraining the intervention-specific changes in the RCP sector. We argue that these must not have changed much because (a) heating appliances are primarily driven by temperature (Z. Liu et al., 2020) (accounted for by our analysis), (b) changes in electricity needs do not enter the local pollutant budget, and (c) less time spent in commercial and public buildings was compensated by more time in residential buildings. Google mobility reports (Alphabet Inc., 2020) based on cell phone tracking suggest a 20 % increase in time spent in the residential sector and a 30 % decrease in the commercial and public sector for Tyrol during the lockdown period. The energy mix in Innsbruck for heating demand is partitioned into residential and “other” (state of Tyrol). The relative contributions to the energy mix for heating in these two broad categories are comparable. In the residential sector it is comprised of 9 % district heating, 34 % oil, 34 % natural gas, 16 % biomass, and 6 % electricity, with the remainder (1 %) being attributed to alternative energy. The category “other” (i.e., everything else) is comprised of 4 % district heating, 37 % oil, 42 % gas, 11 % biomass, and 4 % electricity, with the remainder (2 %) being attributed to alternative energy. Z. Liu et al. (2020) estimated a decline of commercial and residential emissions by 3.6 %; Le Quéré et al. (2020) assumed an increase of residential emissions by 4 % and a decrease in the public sector by 33 % for Europe. As a conservative estimate we bracket changes in the RCP sector activity between 0 % and −20 %, with a best estimate based on the local Google mobility index (−10 %). The observed flux changes can then be partitioned into NOx emissions from vehicular traffic (94-11+2%) and the RCP sector (6-2+11%) accordingly. For CO2, benzene, toluene, and the sum of aromatic NMVOCs we calculate 59-10+7%, 70-7+5%, 94-11+2%, and 87-11+2% arising from vehicular traffic emissions and 41-11+7%, 30-8+5%, 6-1+11%, and 13-3+2%, respectively, coming from the RCP sector. These results suggest that NOx is dominated by vehicular traffic emissions and that CO2 is partitioned more equally between the traffic and RCP sectors. In contrast, urban NMVOC emissions are generally more diverse (Karl et al., 2018). Here we investigate aromatic NMVOCs, which are closely linked to combustion processes and fossil fuel use (EPA, 1998). We observe that toluene and higher aromatic NMVOCs closely track reductions of NOx emissions and vehicular traffic activity. Benzene declined less readily, suggesting that benzene emissions could be more prevalent from the RCP sector. Speciated NMVOC emission factors from residential gas and oil combustion are still quite uncertain, but recent reports from shale gas operations in the US for example indicate a higher contribution of benzene than toluene emissions from natural gas combustion when compared to traffic sources (Gilman et al., 2013; Halliday et al., 2016; Helmig et al., 2014).

After mapping NOx and CO2 emissions from a spatially disaggregated emission model on the seasonal flux footprint (Supplement), the observationally inferred results from above can be compared to the relative attribution of inventory-based emission projections. As for NOx and CO2, the official local bottom-up emission inventory apportions 78 % of NOx fluxes coming from vehicular traffic and 21 % from the RCP sector. For CO2 these relative contributions are 54 % (traffic sector) and 46 % (RCP sector), respectively. These inventory-based results are roughly in line with a recently published bottom-up assessment for CO2 emissions (Le Quéré et al., 2020). We also find that CO2 fluxes are consistent with the relative emission attribution in the inventory but that NOx emissions are significantly overestimated from the RCP sector (e.g., 21 % vs. 6 %) in favor of traffic (Fig. 4). This suggests cleaner NOx combustion sources in the RCP sector and higher NOx emissions from the traffic sector.

The European gas demand has increased significantly over the past decades (European Commission, 2020). As an example, consumption of natural gas increased by about a factor of 4 in Austria (Statistik Austria, 2019) since 1965 and has expanded to 9×109m3. Across Europe growing demand has increased dependence on gas imports, triggering fierce competition between major gas-producing nations (European Commission, 2020). Apart from the power sector, residential demand has contributed significantly to the overall consumption growth across Europe (European Commission, 2020). While residential gas consumption per inhabitant varies quite drastically across European countries, many countries have invested in developing the residential sector towards a higher fraction of natural gas by fuel subsidy policies. Particularly urban areas, where gas infrastructure is in place, have seen significant growth. As an example, the residential energy sector has seen a doubling of the natural gas share for space-heating appliances in western Austria over the last 9 years (Statistik Austria, 2019). In parallel, oil consumption and solid fuel consumption have decreased by about 40 % in the residential sector over the same period. On average, gas represents about a third of the final energy consumption in the residential sector in Austria and across Europe (European Commission, 2020). One of the reasons for promoting natural gas through subsidies in the past was that gas combustion releases about 25 % less CO2 than oil and 40 % less than solid fuels (IEA, 2020). In addition to more efficient energy production, natural gas combustion releases fewer air toxins, such as NOx, CO, NMVOCs, and SO2, when compared to biomass and other solid fuels (EEA, 2019). However, emissions from the RCP sector are quite uncertain and often rely on Tier 1 upscaling methodology (Blain et al., 2019; EEA, 2019). As the European community is committed to transitioning to a carbon-neutral economy (OECD/IEA/NEA/ITF, 2015), the air quality benefit from natural gas in the residential sector needs to be considered, particularly when introducing renewable alternatives such as wood and pellet combustion on a large scale. Our data suggest that the air quality benefit for the release of reactive nitrogen in the RCP sector might have been underestimated in bottom-up emission inventories used for policy making. Official inventory data show that the increase of natural gas combustion in the RCP sector played a significant role in Europe's energy policy. Wood combustion on the other hand would release significant amounts of reactive nitrogen in the gas and aerosol phase depending on fuel N content (Roberts et al., 2020). While pellet combustion is considered cleaner than wood combustion, Tier 1 emission factors for NOx are still about twice as high compared to natural gas combustion, and the release of aerosols is of particular concern (EEA, 2019). When transitioning to a climate-neutral economy, the air quality penalty arising from some renewables needs to be sustainable. From the present analysis we find that the biggest gain for the reduction of urban NOx in Europe remains in the mobility sector and that NOx emissions from the RCP sector are significantly lower than expected. Europe's push towards a diesel-driven car fleet has helped to curb CO2 emissions in the mobility sector but created excess emissions of nitrogen oxides. While the extent of cheating devices used in cars to simulate lower than actual NOx emissions is still unraveling, aggressive reductions of nitrogen oxides are needed to meet Europe's air quality goals (EU-EUR-Lex, 2008). A significant NOx emission reduction in the mobility sector could help counteract potential increases of air pollutants from promoted renewables such as biomass combustion in the future. Urban eddy flux methods present a top-down methodology allowing for quantifying and testing urban sustainability goals of air pollution and climate gas emissions. In combination with an intervention experiment as shown here they can provide a unique and independent verification method of anticipated air quality and climate policy targets.

Code availability

Codes used to analyze eddy covariance data are open source and can also be requested from the corresponding author.


The supplement related to this article is available online at:

Author contributions

TK conceived the overall analysis. TK, CL, and MG designed and performed the field experiments and interpreted the data. MaS conducted the NMVOC flux analysis. MiS assisted with the field experiments. All authors contributed to writing the paper.

Competing interests

The authors declare that they have no conflict of interest.


We thank Christoph Haun of the state of Tyrol, Johann Haun of TIGAS, and Florian Haidacher of the state of Tyrol for providing official emissions, gas consumption, and traffic data referenced in this study.

Financial support

This work was supported by the Austrian Science Fund (FWF, grant no. P30600).

Review statement

This paper was edited by Ralf Sussmann and reviewed by two anonymous referees.


Alphabet Inc.: Google LLC Community Mobility Reports, 2020, available at: (last access: 30 January 2021), 2020. 

Anenberg, S. C., Miller, J., Minjares, R., Du, L., Henze, D. K., Lacey, F., Malley, C. S., Emberson, L., Franco, V., Klimont, Z., and Heyes, C.: Impacts and mitigation of excess diesel-related NOx emissions in 11 major vehicle markets, Nature, 545, 467–471,, 2017. 

Aubinet, M., Vesala, T., and Papale, D.: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, Springer,, 2012. 

Baldocchi, D. D., Hincks, B. B., and Meyers, T. P.: Measuring Biosphere-Atmosphere Exchanges of Biologically Related Gases with Micrometeorological Methods, Ecology, 69, 1331–1340,, 1988. 

Bao, R. and Zhang, A.: Does lockdown reduce air pollution? Evidence from 44 cities in northern China, Sci. Total Environ., 731, 139052,, 2020. 

Blain, W., Buendia, C., Fuglestvedt, E., and Al., E.: Short lived climate forcers, edited by: Blain, W. D., Calvo Buendia, E., Fuglestvedt, J. S., Gómez, D., Masson-Delmotte, V., Tanabe, K., Yassaa, N., Zhai, P., Kranjc, A., Jamsranjav, B., Ngarize., S., Pyrozhenko, Y., Shermanau, P., Connors, S., and Moufouma-Okia, W., Institute for Global Environmental Strategies, 1–66, available at: (last access: 30 January 2021), 2019. 

Carslaw, D. C. and Rhys-Tyler, G.: New insights from comprehensive on-road measurements of NOx, NO2 and NH3 from vehicle emission remote sensing in London, UK, Atmos. Environ., 81, 339–347,, 2013. 

Christen, A.: Atmospheric measurement techniques to quantify greenhouse gas emissions from cities, Urban Clim., 10, 241–260,, 2014. 

Cobourn, W. G.: Accuracy and reliability of an automated air quality forecast system for ozone in seven Kentucky metropolitan areas, Atmos. Environ., 41, 5863–5875,, 2007. 

Dabberdt, W. F., Lenschow, D. H., Horst, T. W., Zimmerman, P. R., Oncley, S. P., and Delany, A. C.: Atmosphere-Surface Exchange Measurements, Science, 260, 1472–1481,, 1993. 

Duffy, N. and Helmbold, D.: Boosting Methods for Regression, Mach. Learn., 47, 153–200,, 2002. 

EEA: EMEP/EEA air pollutant emission inventory guidebook, available at: (last access: 30 January 2021), 2019. 

Elith, J., Leathwick, J. R., and Hastie, T.: A working guide to boosted regression trees, J. Anim. Ecol., 77, 802–813,, 2008. 

EPA: Locating and estimating air emissions from sources of benzene, EPA report: EPA-454/R-98-011, 1998. 

EU-EUR-Lex: Directives 1999/94/EC and 2008/50/EC, available at: (last access: 30 January 2021), 2008. 

European Commission: Gas and electricity market reports, available at: (last access: 30 January 2021), 2020. 

Foken, T. and Wichura, B.: Tools for quality assessment of surface-based flux measurements, Agr. Forest. Meteorol., 78, 83–105,, 1996. 

Fowler, D., Pilegaard, K., Sutton, M. A., Ambus, P., Raivonen, M., Duyzer, J., Simpson, D., Fagerli, H., Fuzzi, S., Schjoerring, J. K., Granier, C., Neftel, A., Isaksen, I. S. A., Laj, P., Maione, M., Monks, P. S., Burkhardt, J., Daemmgen, U., Neirynck, J., Personne, E., Wichink-Kruit, R., Butterbach-Bahl, K., Flechard, C., Tuovinen, J. P., Coyle, M., Gerosa, G., Loubet, B., Altimir, N., Gruenhage, L., Ammann, C., Cieslik, S., Paoletti, E., Mikkelsen, T. N., Ro-Poulsen, H., Cellier, P., Cape, J. N., Horváth, L., Loreto, F., Niinemets, U., Palmer, P. I., Rinne, J., Misztal, P., Nemitz, E., Nilsson, D., Pryor, S., Gallagher, M. W., Vesala, T., Skiba, U., Br ggemann, N., Zechmeister-Boltenstern, S., Williams, J., ODowd, C., Facchini, M. C., de Leeuw, G., Flossman, A., Chaumerliac, N., and Erisman, J. W.: Atmospheric composition change: Ecosystems Atmosphere interactions, Atmos. Environ., 43, 5193–5267,, 2009. 

Franco, V., Posada Sanches, F., German, J., and Mock, P.: Real-world exhaust emissions from modern diesel cars, ICCT, 1–52, available at: (last access: 26 February 2021), 2014. 

Gilman, J. B., Lerner, B. M., Kuster, W. C., and de Gouw, J. A.: Source Signature of Volatile Organic Compounds from Oil and Natural Gas Operations in Northeastern Colorado, Environ. Sci. Technol., 47, 1297–1305,, 2013. 

Grange, S. K. and Carslaw, D. C.: Using meteorological normalisation to detect interventions in air quality time series, Sci. Total Environ., 653, 578–588,, 2019. 

Halliday, H. S., Thompson, A. M., Wisthaler, A., Blake, D. R., Hornbrook, R. S., Mikoviny, T., Mueller, M., Eichler, P., Apel, E. C., and Hills, A. J.: Atmospheric benzene observations from oil and gas production in the Denver-Julesburg Basin in July and August 2014, J. Geophys. Res.-Atmos., 121, 1111–5574,, 2016. 

Helmig, D., Thompson, C. R., Evans, J., Boylan, P., Hueber, J., and Park, J.-H.: Highly Elevated Atmospheric Levels of Volatile Organic Compounds in the Uintah Basin, Utah, Environ. Sci. Technol., 48, 4707–4715,, 2014. 

IEA: Fuels and technologies, available at: (last access: 30 January 2021), 2020. 

Im, U., Bianconi, R., Solazzo, E., Kioutsioukis, I., Badia, A., Balzarini, A., Baro, R., Bellasio, R., Brunner, D., Chemel, C., Curci, G., Flemming, J., Forkel, R., Giordano, L., Jimenez-Guerrero, P., Hirtl, M., Hodzic, A., Honzak, L., Jorba, O., Knote, C., Kuenen, J. J. P., Makar, P. A., Manders-Groot, A., Neal, L., Perez, J. L., Pirovano, G., Pouliot, G., San Jose, R., Savage, N., Schroder, W., Sokhi, R. S., Syrakov, D., Torian, A., Tuccella, P., Werhahn, J., Wolke, R., Yahya, K., Zabkar, R., Zhang, Y., Zhang, J., Hogrefe, C., and Galmarini, S.: Evaluation of operational on-line-coupled regional air quality models over Europe and North America in the context of AQMEII phase 2. Part I: Ozone, Atmos. Environ., 115, 404–420,, 2015. 

Karl, T., Guenther, A., Jordan, A., Fall, R., and Lindinger, W.: Eddy covariance measurement of biogenic oxygenated VOC emissions from hay harvesting, Atmos. Environ., 35, 491–495,, 2001. 

Karl, T., Graus, M., Striednig, M., Lamprecht, C., Hammerle, A., Wohlfahrt, G., Held, A., Von Der Heyden, L., Deventer, M. J., Krismer, A., Haun, C., Feichter, R., and Lee, J.: Urban eddy covariance measurements reveal significant missing NOx emissions in Central Europe, Sci. Rep., 7, 2536,, 2017. 

Karl, T., Striednig, M., Graus, M., Hammerle, A., and Wohlfahrt, G.: Urban flux measurements reveal a large pool of oxygenated volatile organic compound emissions, P. Natl. Acad. Sci. USA, 115, 1186–1191,, 2018. 

Karl, T., Gohm, A., Rotach, M., Ward, H., Graus, M., Cede, A., Wohlfahrt, G., Hammerle, A., Haid, M., Tiefengraber, M., Lamprecht, C., Vergeiner, J., Kreuter, A., Wagner, J., and Staudinger, M.: Studying Urban Climate and Air quality in the Alps – The Innsbruck Atmospheric Observatory, B. Am. Meteorol. Soc., 101, E488–E507,, 2020. 

Kljun, N., Calanca, P., Rotach, M. W., and Schmid, H. P.: A simple two-dimensional parameterisation for Flux Footprint Prediction (FFP), Geosci. Model Dev., 8, 3695–3713,, 2015. 

Langford, B., Davison, B., Nemitz, E., and Hewitt, C. N.: Mixing ratios and eddy covariance flux measurements of volatile organic compounds from an urban canopy (Manchester, UK), Atmos. Chem. Phys., 9, 1971–1987,, 2009. 

Laughner, J. L. and Cohen, R. C.: Direct observation of changing NOx lifetime in North American cities, Science, 366, 723–727,, 2019. 

Lee, J. D., Helfter, C., Purvis, R. M., Beevers, S. D., Carslaw, D. C., Lewis, A. C., Moller, S. J., Tremper, A., Vaughan, A., and Nemitz, E. G.: Measurement of NOx Fluxes from a Tall Tower in Central London, UK and Comparison with Emissions Inventories, Environ. Sci. Technol., 49, 1025–1034,, 2015. 

Lenschow, D. H. and Delany, A. C.: An analytic formulation for NO and NO2 flux profiles in the atmospheric surface layer, J. Atmos. Chem., 5, 301–309,, 1987. 

Lenschow, D. H., Gurarie, D., and Patton, E. G.: Modeling the diurnal cycle of conserved and reactive species in the convective boundary layer using SOMCRUS, Geosci. Model Dev., 9, 979–996,, 2016. 

Le Quéré, C. , Jackson, R. B., Jones, M. W., Smith, A. J. P., Abernethy, S., Andrew, R. M., De-Gol, A. J., Willis, D. R., Shan, Y., Canadell, J. G., Friedlingstein, P., Creutzig, F., and Peters, G. P.: Temporary reduction in daily global CO2 emissions during the COVID-19 forced confinement, Nat. Clim. Change, 10, 647–653,, 2020. 

Liu, F., Page, A., Strode, S. A., Yoshida, Y., Choi, S., Zheng, B., Lamsal, L. N., Li, C., Krotkov, N. A., Eskes, H., van der A, R., Veefkind, P., Levelt, P. F., Hauser, O. P., and Joiner, J.: Abrupt decline in tropospheric nitrogen dioxide over China after the outbreak of COVID-19, Sci. Adv., 6, 1–5,, 2020. 

Liu, Z., Ciais, P., Deng, Z., Lei, R., Davis, S. J., Feng, S., Zheng, B., Cui, D., Dou, X., Zhu, B., Guo, R., Ke, P., Sun, T., Lu, C., He, P., Wang, Y., Yue, X., Wang, Y., Lei, Y., Zhou, H., Cai, Z., Wu, Y., Guo, R., Han, T., Xue, J., Boucher, O., Boucher, E., Chevallier, F., Tanaka, K., Wei, Y., Zhong, H., Kang, C., Zhang, N., Chen, B., Xi, F., Liu, M., Breon, F., Lu, Y., Zhang, Q., Guan, D., Gong, P., Kammen, D., He, K., and Schellnhuber, H.: Near-real-time monitoring of global CO2 emissions reveals the effects of the COVID-19 pandemic, Nat. Commun., 11, 5172,, 2020. 

Menut, L., Bessagnet, B., Siour, G., Mailler, S., Pennel, R., and Cholakian, A.: Impact of lockdown measures to combat Covid-19 on air quality over western Europe, Sci. Total Environ., 741, 140426,, 2020. 

NAS: The Future of Atmospheric Chemistry Research: Remembering Yesterday, Understanding Today, Anticipating Tomorrow, The National Academies Press, Washington, D.C., 2016. 

Nemitz, E., Jimenez, J. L., Huffman, J. A., Ulbrich, I. M., Canagaratna, M. R., Worsnop, D. R., and Guenther, A. B.: An Eddy-Covariance System for the Measurement of Surface/Atmosphere Exchange Fluxes of Submicron Aerosol Chemical Species – First Application Above an Urban Area, Aerosol Sci. Tech., 42, 636–657,, 2008. 

OECD/IEA/NEA/ITF: Aligning Policies for a Low-carbon Economy, OECD Publishing, Paris, ISBN 978-92-64-23326-3,, 2015. 

Oncley, S. P., Foken, T., Vogt, R., Kohsiek, W., DeBruin, H. A. R., Bernhofer, C., Christen, A., van Gorsel, E., Grantz, D., Feigenwinter, C., Lehner, I., Liebethal, C., Liu, H., Mauder, M., Pitacco, A., Ribeiro, L., and Weidinger, T.: The Energy Balance Experiment EBEX-2000. Part I: overview and energy balance, Bound.-Lay. Meteorol., 123, 1–28,, 2007. 

Patton, E. G., Horst, T. W., Sullivan, P. P., Lenschow, D. H., Oncley, S. P., Brown, W. O. J., Burns, S. P. P., Guenther, A. B., Held, A., Karl, T., Mayor, S. D., Rizzo, L. V., Spuler, S. M., Sun, J., Turnipseed, A. A., Awine, E. J., Edburg, S. L., Lamb, B. K., Avissar, R., Calhoun, R. J., Kleissl, J., Massman, W. J., Paw U, K. T., and Weil, J. C.: The canopy horizontal array turbulence study, B. Am. Meteorol. Soc., 92, 593–611, 2011. 

Prybutok, V. R., Yi, J., and Mitchell, D.: Comparison of neural network models with ARIMA and regression models for prediction of Houstons daily maximum ozone concentrations, Eur. J. Oper. Res., 122, 31–40,, 2000. 

Rannik, Ü., Altimir, N., Mammarella, I., Bäck, J., Rinne, J., Ruuskanen, T. M., Hari, P., Vesala, T., and Kulmala, M.: Ozone deposition into a boreal forest over a decade of observations: evaluating deposition partitioning and driving variables, Atmos. Chem. Phys., 12, 12165–12182,, 2012. 

Rantala, P., Järvi, L., Taipale, R., Laurila, T. K., Patokoski, J., Kajos, M. K., Kurppa, M., Haapanala, S., Siivola, E., Petäjä, T., Ruuskanen, T. M., and Rinne, J.: Anthropogenic and biogenic influence on VOC fluxes at an urban background site in Helsinki, Finland, Atmos. Chem. Phys., 16, 7981–8007,, 2016. 

Rinne, H. J. I., Guenther, A. B., Warneke, C., de Gouw, J. A., and Luxembourg, S. L.: Disjunct eddy covariance technique for trace gas flux measurements, Geophys. Res. Lett., 28, 3139–3142,, 2001. 

Roberts, J. M., Stockwell, C. E., Yokelson, R. J., de Gouw, J., Liu, Y., Selimovic, V., Koss, A. R., Sekimoto, K., Coggon, M. M., Yuan, B., Zarzana, K. J., Brown, S. S., Santin, C., Doerr, S. H., and Warneke, C.: The nitrogen budget of laboratory-simulated western US wildfires during the FIREX 2016 Fire Lab study, Atmos. Chem. Phys., 20, 8807–8826,, 2020. 

Robeson, S. M. and Steyn, D. G.: Evaluation and comparison of statistical forecast models for daily maximum ozone concentrations, Atmos. Environ. B, 24, 303–312,, 1990. 

Schiermeier, Q.: Why pollution is plummeting in some cities, but not others, Nature, 580, 313–314,, 2020. 

Spirig, C., Neftel, A., Ammann, C., Dommen, J., Grabmer, W., Thielmann, A., Schaub, A., Beauchamp, J., Wisthaler, A., and Hansel, A.: Eddy covariance flux measurements of biogenic VOCs during ECHO 2003 using proton transfer reaction mass spectrometry, Atmos. Chem. Phys., 5, 465–481,, 2005. 

Squires, F. A., Nemitz, E., Langford, B., Wild, O., Drysdale, W. S., Acton, W. J. F., Fu, P., Grimmond, C. S. B., Hamilton, J. F., Hewitt, C. N., Hollaway, M., Kotthaus, S., Lee, J., Metzger, S., Pingintha-Durden, N., Shaw, M., Vaughan, A. R., Wang, X., Wu, R., Zhang, Q., and Zhang, Y.: Measurements of traffic-dominated pollutant emissions in a Chinese megacity, Atmos. Chem. Phys., 20, 8737–8761,, 2020. 

Statistik Austria: Energiebilanzen, available at: (last access: 30 January 2021), 2019. 

Steinbacher, M., Zellweger, C., Schwarzenbach, B., Bugmann, S., Buchmann, B., Ordóñez, C., Prevot, A. S. H., and Hueglin, C.: Nitrogen oxide measurements at rural sites in Switzerland: Bias of conventional measurement techniques, J. Geophys. Res., 112, D11307,, 2007.  

Striednig, M., Graus, M., Märk, T. D., and Karl, T. G.: InnFLUX – an open-source code for conventional and disjunct eddy covariance analysis of trace gas measurements: an urban test case, Atmos. Meas. Tech., 13, 1447–1465,, 2020. 

Sussmann, R. and Rettinger, M.: Can We Measure a COVID-19-Related Slowdown in Atmospheric CO2 Growth? Sensitivity of Total Carbon Column Observations, Remote Sens., 12, 1–22,, 2020. 

UBA: Bundeslaender Luftschadstoffinventur 1990–2017, available at: (last access: 30 January 2021), 2019. 

UN: World Population Prospects 2019, Volume II: Demographic Profiles, available at: (last access: 30 January 2021), 2019. 

Vaughan, A. R., Lee, J. D., Misztal, P. K., Metzger, S., Shaw, M. D., Lewis, A. C., Purvis, R. M., Carslaw, D. C., Goldstein, A. H., Hewitt, C. N., Davison, B., Beevers, S. D., and Karl, T. G.: Spatially resolved flux measurements of NOx from London suggest significantly higher emissions than predicted by inventories, Faraday Discuss., 189, 455–472, 2016. 

Vaughan, A. R., Lee, J. D., Shaw, M. D., Misztal, P. K., Metzger, S., Vieno, M., Davison, B., Karl, T. G., Carpenter, L. J., Lewis, A. C., Purvis, R. M., Goldstein, A. H., and Hewitt, C. N.: VOC emission rates over London and South East England obtained by airborne eddy covariance, Faraday Discuss., 200, 599–620,, 2017. 

Velasco, E., Pressley, S., Allwine, E., Westberg, H., and LAmb, B.: Measurements of CO fluxes from the Mexico City urban landscape, Atmos. Environ., 39, 7433–7446,, 2005. 

Velasco, E., Pressley, S., Grivicke, R., Allwine, E., Coons, T., Foster, W., Jobson, B. T., Westberg, H., Ramos, R., Hernández, F., Molina, L. T., and Lamb, B.: Eddy covariance flux measurements of pollutant gases in urban Mexico City, Atmos. Chem. Phys., 9, 7325–7342,, 2009. 

Ward, H. C., Kotthaus, S., Grimmond, C. S. B., Bjorkegren, A., Wilkinson, M., Morrison, W. T. J., Evans, J. G., Morison, J. I. L., and Iamarino, M.: Effects of urban density on carbon dioxide exchanges: Observations of dense urban, suburban and woodland areas of southern England, Environ. Pollut., 198, 186–200,, 2015. 

WHO: Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19), available at: (last access: 30 January 2021), 2020. 

Short summary
The first European SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) wave and associated lockdown provided a unique sensitivity experiment to study air pollution. We find significantly different emission trajectories between classical air pollution and climate gases (e.g., carbon dioxide). The analysis suggests that European policies, shifting residential, public, and commercial energy demand towards cleaner combustion, have helped to improve air quality more than expected.
Final-revised paper