Articles | Volume 22, issue 4
Research article
01 Mar 2022
Research article |  | 01 Mar 2022

Variability of nitrogen oxide emission fluxes and lifetimes estimated from Sentinel-5P TROPOMI observations

Kezia Lange, Andreas Richter, and John P. Burrows

Satellite observations of the high-resolution TROPOspheric Monitoring Instrument (TROPOMI) on Sentinel-5 Precursor can be used to observe nitrogen dioxide (NO2) at city scales to quantify short time variability of nitrogen oxide (NOx) emissions and lifetimes on a daily and seasonal basis. In this study, 2 years of TROPOMI tropospheric NO2 columns, having a spatial resolution of up to 3.5 km× 5.5 km, have been analyzed together with wind and ozone data. NOx lifetimes and emission fluxes are estimated for 50 different NOx sources comprising cities, isolated power plants, industrial regions, oil fields, and regions with a mix of sources distributed around the world. The retrieved NOx emissions are in agreement with other TROPOMI-based estimates and reproduce the variability seen in power plant stack measurements but are in general lower than the analyzed stack measurements and emission inventory results. Separation into seasons shows a clear seasonal dependence of NOx emissions with in general the highest emissions during winter, except for isolated power plants and especially sources in hot desert climates, where the opposite is found. The NOx lifetime shows a systematic latitudinal dependence with an increase in lifetime from 2 to 8 h with latitude but only a weak seasonal dependence. For most of the 50 sources including the city of Wuhan in China, a clear weekly pattern of NOx emissions is found, with weekend-to-weekday ratios of up to 0.5 but with a high variability for the different locations. During the Covid-19 lockdown period in 2020, strong reductions in the NOx emissions were observed for New Delhi, Buenos Aires, and Madrid.

1 Introduction

Nitrogen oxides (NOx=NO+NO2) play a key role in atmospheric chemistry, air quality, and climate. They are emitted into the atmosphere by both natural processes and human activity. Natural sources include lightning, microbial processes in soils, and naturally occurring wildfires. The dominant source of NOx is fossil-fuel combustion from traffic, residential heating, cooking, and the industry and energy sectors. These sources are concentrated in cities and urban areas. In addition, biomass burning and the use of fertilizers are also significant sources of NOx (Jacob1999; Seinfeld and Pandis2006; Stocker2014).

These sources predominantly release nitrogen monoxide (NO). Nitrogen dioxide (NO2) is released in smaller amounts but is rapidly produced in the atmosphere, where NO reacts with ozone (O3). During the day NO2 is photolyzed, reforming NO and producing an oxygen atom O, which forms O3 in a termolecular reaction (Seinfeld and Pandis2006).


In the troposphere NOx are precursors of the health hazard and greenhouse gas O3, but NO2 is also an important health hazard in its own right; see Jacob (1999) and references therein. Consequently, monitoring and understanding its behavior is of particular importance in cities and urban agglomerations, where high NOx emissions from multiple sources are found in combination with high population density. As a result, a large part of the population is exposed to polluted air (Molina and Molina2004).

NOx are short lived in the atmosphere, with a lifetime of several hours in the boundary layer during daytime. This explains in part the high spatial and temporal variability observed for NO2. Other factors leading to large concentration gradients in NO2 are the heterogeneous distribution of sources and variations of meteorological parameters, such as wind speed, temperature, and illumination, which play key roles in determining the atmospheric lifetime and dilution of NO2 (Beirle et al.2003; Stavrakou et al.2008).

Measurements of NO2 with adequate spatial and temporal resolution are required to assess and compare the variability of NOx emissions and lifetimes around the world. Provided that tropospheric NO2 columns retrieved from satellite sensors have sufficient spatial resolution, the typically short daytime lifetime of NOx in urban agglomerations provides an opportunity to disentangle and quantify local sources of NOx and their variations over time. The TROPOspheric Monitoring Instrument (TROPOMI) on Sentinel-5 Precursor (S5P), which was launched in October 2017, provides higher spatial resolution than its predecessors, e.g., the Ozone Monitoring Instrument (OMI). Tropospheric columns of NO2 retrieved from TROPOMI thus offer the best opportunity so far to deconvolve urban daytime sources of NOx (Veefkind et al.2012; Griffin et al.2019). Making use of the NO2 tropospheric columns from TROPOMI but also from its predecessors, a number of studies on the variability of NOx have been carried out.

The high spatial resolution of TROPOMI makes it possible to identify much better NOx point sources and quantify their emissions. This provides the opportunity to control emission inventories and improve on them (Beirle et al.2019, 2021).

Investigation of seasonal emission estimates enables the disentanglement of the NOx sources in a region and the identification of their individual contributions (van der A et al.2008). The analysis of TROPOMI data for Paris showed that residential heating, rather than transport emissions, is the main source of NOx in winter. This is not well accounted for in current emission inventories, in which residential heating is underestimated during winter and overestimated during summer (Lorente et al.2019).

NOx emissions also change from day to day as a result of human behavior, especially between work and rest days. Such patterns are readily identified (Beirle et al.2003). A recent study by Stavrakou et al. (2020) investigated the weekly NO2 cycle and its trends using the long-term satellite observations of OMI and 1 year of TROPOMI NO2 column measurements. Significant weekly cycles were found. A weakening trend over Europe and the US could be observed. The opposite behavior was found for regions with increasing NOx emissions. The decreases or, respectively, increases in the contribution of anthropogenic NOx emissions to the observed NO2 levels are thus revealed. With TROPOMI data it is possible to use relatively short periods of data to analyze day-to-day variability. For Chicago a clear weekend effect with reduced NOx emissions of 30 % during weekends was found by analyzing one season of TROPOMI measurements (Goldberg et al.2019). Lorente et al. (2019) investigated the day-to-day variability of NOx emissions in Paris for individual days, with the highest emissions on cold weekdays and the lowest on warm weekend days.

Beginning in early 2020, first China and subsequently the majority of countries around the world took containment measures to limit the spread of Covid-19. These measures resulted in significant changes in human behavior, with reductions in traffic density and industrial activity and consequently anthropogenic emissions. Temporally fine-resolved emission estimates help to identify the different potential origins of such changes. This provides an approach and an opportunity to better quantify source contributions and to distinguish between different anthropogenic and natural sources of NOx. Several recent studies have analyzed satellite measurements of NO2 columns and report substantial decreases in the NO2 tropospheric column in early 2020 over China (Liu et al.2020; Bauwens et al.2020), northern Italy, South Korea, and the United States (Bauwens et al.2020). However, because of the high variability of NO2 columns these reductions cannot be simply directly attributed to a decrease in NOx emissions resulting from the Covid-19 containment measures. The tropospheric column of NO2 is influenced by behavioral patterns of anthropogenic activity, seasonality, and meteorology. To account for meteorology, which led to lower NO2 values in spring 2020 as compared to spring 2019 in North America, Goldberg et al. (2020) combined TROPOMI NO2 column data with meteorological data and a chemical transport model and calculated normalized NO2 changes that provide a better representation of Covid-19-related NOx emission changes.

In addition to NOx emissions, NOx lifetimes are determined from tropospheric NO2 columns (Leue et al.2001; Beirle et al.2003; Kunhikrishnan et al.2004). The NOx lifetime in daylight depends on the rate of loss of NOx which is attributed primarily to the reaction of the hydroxyl radical (OH) with NO2 to form HNO3. Consequently, there is a nonlinear relationship between the tropospheric concentrations of OH and NOx, which again depends in a complex fashion on the NOx concentration (Valin et al.2013; Stavrakou et al.2008). Laughner and Cohen (2019) showed that NOx lifetime changed significantly between 2005 and 2014 in North American cities, changes being of the same order of magnitude as those in NOx emissions. Another important factor influencing NOx lifetime is the actinic radiation photolyzing NO2, which varies diurnally and is modulated by the presence of clouds. Shorter lifetimes at smaller solar zenith angles in summer or at lower latitudes due to higher photolysis frequencies are expected (Stavrakou et al.2020). Typical lifetimes of NOx range from 2 to 8 h in polluted air masses and extend to about a day for cleaner, more rural background concentrations for which nighttime chemistry also has to be considered (Beirle et al.2011; Valin et al.2014; de Foy et al.2014; Seinfeld and Pandis2006). Studies for winter months and especially for winter months at higher latitudes are limited; analyses with the chemical transport model GEOS-Chem show a tendency towards longer lifetimes of about 1 d (Martin et al.2003; Shah et al.2020).

In this study, we use the method first developed by Beirle et al. (2011) and refined by later studies (Pommier et al.2013; Valin et al.2013) to estimate the NOx emissions and lifetimes from TROPOMI NO2 column amounts. Studies of NOx emissions and lifetimes have been reported using this method for OMI data, which has coarser spatial resolution, for a limited number of sources and using long periods of data (Ialongo et al.2014; de Foy et al.2015; Lu et al.2015; Liu et al.2016). Other studies which applied this method to TROPOMI data are reported by Goldberg et al. (2019) or Lorente et al. (2019). In this work, we provide the first study of NOx emissions and lifetimes for a large data set, comprising 50 NOx source regions. These are well distributed over the world, consisting of cities, isolated power plants, industrial regions, oil fields, and regions with a mix of sources. We show that the spatial resolution and good signal-to-noise ratio of TROPOMI NO2 columns allow us to investigate 2 years of data and additionally divide and study the following periods: seasons, weekend and working days, and times before and during Covid-19 restrictions to examine short-term variability. The focus of our investigation is to assess the variability of NOx emissions and lifetimes in space and time during the period of observation.

2 Data

2.1 TROPOMI NO2 tropospheric vertical column

In October 2017, the Copernicus S5P satellite with TROPOMI on board was launched into a Sun-synchronous orbit at 824 km altitude. TROPOMI comprises a hyperspectral spectrometer measuring radiation in the ultraviolet (270–320 nm), visible (310–500 nm), and infrared (675–775 nm, 2305–2385 nm) spectral regions (Veefkind et al.2012). The two-dimensional charge-coupled device (CCD) detectors measure spectra in 450 separate viewing directions across the 2600 km swath with an integration time of approximately 1 s. This results in a high spatial resolution of 3.5 km× 7 km in nadir with little variation across the swath. In August 2019, the pixel size was reduced further to 3.5 km× 5.5 km by reducing along-track averaging. One orbit around the Earth takes about 100 min, which, in combination with the wide swath, results in a daily global coverage. At higher latitudes, on some days up to three overpasses, with approximately 100 min in between the measurements, are available. In the rare case of three overpasses, two of them are on the outer edge of the swath. The Equator overpass time of S5P is 13:30 mean local solar time in the ascending node.

In this study, the operational NO2 product of TROPOMI is used (van Geffen et al.2018). The NO2 product from TROPOMI extends NO2 measurements from earlier instruments such as GOME (1995–2011, Burrows et al.1999), SCIAMACHY (2002–2012, Bovensmann et al.1999), GOME-2 (since 2007, Munro et al.2006), OMI (since 2004, Levelt et al.2006), and OMPS (since 2011, Dittman et al.2002). These instruments had increasing spatial coverage and resolution, OMI being the first instrument with daily global coverage. However, with a resolution of 13 km× 24 km at nadir, its spatial resolution is 1 order of magnitude poorer than that of TROPOMI in the center of the swath and much lower at the edges. Differences in resolution are larger at the edges, which is due to the fact that for TROPOMI averages of two pixels in the center of the scan, but not at the edges, are used to compensate for the geometric effect caused by the slanted view to the ground.

The level-1b spectra measured by TROPOMI are analyzed with the differential optical absorption spectroscopy (DOAS) technique in the fitting window 405–465 nm. The retrieved NO2 slant column densities are separated into a stratospheric part and a tropospheric part based on data assimilation by the TM5-MP global chemistry transport model. The resulting tropospheric slant columns are then converted to tropospheric vertical columns by tropospheric air mass factors (AMFs) using a look-up table of altitude-dependent AMFs, NO2 vertical profiles from the TM5-MP model, the OMI climatological surface albedo, and cloud characteristics derived using the FRESCO algorithm (van Geffen et al.2018). The final product is the tropospheric vertical column, defined as the vertically integrated number of molecules per unit area between the surface and the tropopause.

In this study, we use the operational level-2 TROPOMI tropospheric vertical column NO2 product from March 2018 to November 2020, the first 2 years of data for the general study, and additional months from 2020 for the Covid-19 impact study. This includes the reprocessed (RPRO) and offline-mode (OFFL) data of version V01.00.01 to version V01.03.02. Changes between the versions before V01.04 are minor, and the data products are comparable and can be mixed. V01.04 was activated on 2 December 2020 and implemented major changes, which led to a substantial increase in the tropospheric NO2 column over polluted areas for scenes with small cloud fractions. We only use the data up to the end of November 2020 for our analysis to ensure better comparability, since mixing data from before version V1.04 and after is not recommended. Further major updates in the TROPOMI NO2 retrieval algorithm were made in V02.02 released on 1 July 2021 (van Geffen et al.2021). A complete mission reprocessing will be performed to deliver a harmonized data set (Eskes and Eichmann2021). Data from 30 April 2018 onwards are freely available on (last access: 21 February 2022). Data before 30 April are only available on the S5P Copernicus Expert Hub and therefore not publicly accessible. The data with a spatial resolution of up to 3.5 km × 5.5 km are oversampled to a finer resolution of 0.01× 0.01. Each TROPOMI ground pixel is accompanied by a quality assurance value (qa_value). The qa_value ranges from zero (error, no output) to one (no errors or warnings) and indicates the quality of the processing and retrieval result. Based on the recommendation by van Geffen et al. (2018), measurements with a qa_value of lower than 0.75 are filtered and not used. A qa_value of 0.75 removes problematic retrievals, errors, partially snow-/ice-covered scenes, and measurements with cloud radiance fractions of more than 50 %.

2.2 Wind data

In addition to the TROPOMI NO2 data, wind speed and direction are required to determine the NOx emissions and lifetimes. The wind data used are provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis hourly data with a horizontal resolution of 0.25× 0.25 on model levels. To merge the wind data in space and time with the TROPOMI observations, they are interpolated to the overpass time and oversampled to the same 0.01× 0.01 resolution as the TROPOMI data. The chosen height of wind data can be critical for wind speed and direction, because of increasing wind speeds and turning with height, and therefore also for NOx emissions and lifetime estimates. Beirle et al. (2011) investigated the dependence of the estimated NOx emissions and lifetimes on the wind data level height by comparing the results calculated with wind fields averaged from ground up to 200, 500, and 1000 m and showed little dependence of the results on the wind-level height. The resulting emissions and lifetimes change by less than 2 % (5 %) on average when using 200 m (1000 m) instead of 500 m and by less than 15 % for individual sources. We use wind data averaged over the boundary layer, with the boundary layer height also provided by ERA5. For the early afternoon overpass of TROPOMI, we assume that the boundary layer is well mixed and that wind information averaged over this layer is representative of the NOx emissions and lifetimes investigated in this study. Using wind data averaged over the boundary layer instead of using a fixed height has the advantage that seasonality in wind data due to the seasonal variability of the boundary layer height is accounted for.

2.3 Ozone volume mixing ratios

Ozone volume mixing ratios, for converting the NO2 column measurements to NOx columns to estimate emissions in terms of NOx, were taken from the ECMWF ERA5 reanalysis. The ozone data have the same resolution as the used wind data, hourly data with a horizontal resolution of 0.25× 0.25 on model levels. They are averaged over the boundary layer, interpolated to TROPOMI overpass time, and oversampled to the same 0.01 resolution as the TROPOMI data in the same manner as for the wind data.

Figure 1(a) Plot of the mean TROPOMI NO2 tropospheric vertical column (VC) from 1 March 2018 to 29 February 2020 in the region of the power plants Medupi and Matimba in South Africa for days with wind speed >2m s−1. (b) In this plot each NO2 tropospheric VC is converted to an estimate of NOx VC and is rotated with its wind direction around the center of the source (cross) to an upwind–downwind pattern. The black lines indicate the sector of ± 50 km around the source relevant for the emission and lifetime estimates. (c) NOx VC integrated perpendicular to the wind direction (NOx line density) as a function of distance to the source calculated for the ± 50 km sector (gray) and fit results (black) with estimated NOx emissions and lifetime. Emission and lifetime uncertainties are 1 standard deviation retrieved by the fitting procedure.


2.4 Emission inventory

For comparison with the calculated TROPOMI NOx emissions, we use the Emissions Database for Global Atmospheric Research (EDGAR) v5.0 of 2015. The bottom-up emission data are available at 0.1× 0.1 gridded resolution. For calculation of NOx emissions in the specific source regions, the gridded emission inventory data have to be integrated over the area for which the top-down method based on TROPOMI data is sensitive. The EDGAR inventory is based on activity data (i.e., population, energy, fossil fuel production, industrial processes, agricultural statistics), mainly from the International Energy Agency (IEA), corresponding emission factors, and national and regional information on technology mix data and end-of-pipe measurements (Crippa et al.2018, 2019). Uncertainties in bottom-up inventories are inferred from the dependence of emission factors on the fuel type, technology and combustion condition, as well as the low-resolution activity data and emission factors. The uncertainty in NOx emissions from the latest EDGAR inventory version (v4.3.2) varies from 17 % to 69 % for different regions (Crippa et al.2018). The limited temporal coverage of bottom-up emissions results in additional uncertainties; 2015 is the most recent year available in EDGAR v5.0, which does not reflect the recent negative and positive trends, which were found from trend analysis of NO2 column satellite data (Georgoulias et al.2019). Since a large part of the regions analyzed in this study is located in industrialized and highly populated regions, where NO2 has generally decreased according to Georgoulias et al. (2019), we anticipate that the TROPOMI estimates for the majority of the analyzed regions will be lower than the EDGAR estimates for 2015.

3 Method

The method to estimate NOx emissions and lifetimes from satellite column data builds on the heritage of the method introduced by Beirle et al. (2011) and refined by later studies (Pommier et al.2013; Valin et al.2013). We use 2 years of TROPOMI NO2 tropospheric vertical column data from 1 March 2018 to 29 February 2020 for the general analysis, excluding data for which Covid-19 regulations influenced emissions. For the two Chinese cities included in this study, the analyzed period is limited to 22 January 2020 due to early Covid-19 regulations. The Covid-19 impact study in Sect. 4.5 is based on the months from January to November of 2019 and 2020. The three steps of the analysis are shown in Fig. 1 for the Medupi and Matimba power plants in South Africa as an example. First, the source region is selected. The choice of sources is described in more detail in Sect. 3.1. The next step is a rotation of the satellite measurements around the selected source with the corresponding ERA5 reanalysis wind data to a common wind direction resulting in an upwind–downwind pattern, which is described in more detail in Sect. 3.2. The NO2 column of each pixel is converted into an NOx column, and the mean NOx distribution is calculated. In the last step, we apply a fit with an exponentially modified Gaussian (EMG) function to the averaged NOx columns to calculate the NOx emissions and lifetime (see Sect. 3.3).

3.1 Selection of sources

To obtain a representative analysis of the variability of NOx emissions and lifetimes from urban and industrial areas around the globe, we have used the following guideline for the selection of targets. The targeted emission sources need to be well distributed around the world to provide information about the seasonal variations, climatic conditions, latitudinal dependence, and weekday behavior of the NOx lifetimes and emissions. By visually inspecting the global mean NO2 tropospheric vertical column distribution from March 2018 to February 2020 (shown in Fig. 2), we selected 50 target regions, which are marked with red circles, for the calculation of NOx emissions and lifetimes. The selected sources comprise a mix of cities with predominantly domestic and transport emissions, e.g., Paris (France), cities with more industrial emitters, e.g., Chelyabinsk (Russia), power plants like Medupi and Matimba (South Africa), or oil refinery regions like the Sarir Field (Libya). The method works best for isolated point sources with a high contrast between source and background. Regions with few clouds are preferred to maximize the number of satellite observations. In addition, a local meteorology for the target region, with rather homogeneous wind patterns, facilitates the detection of the outflow patterns, which is further enhanced by filtering out data with wind speeds of less than 2 m s−1 (Beirle et al.2011). Some initially promising sources like Santiago de Chile were omitted. This was a result of their location in coastal and mountainous regions with inhomogeneous terrain, resulting in inhomogeneous wind patterns, which are more difficult to interpret and lead to larger uncertainties in the NOx emission rates and lifetimes. For this study, many of the sources with high NO2 signals in China were not used, because the influence of large NOx-emitting sources nearby resulted in low contrast of the NO2 column amount between the target and the local background. For such conditions, other methods to determine NOx emission rates and lifetimes are more appropriate (Liu et al.2016). Despite efforts to select a broad range of regions in many respects, some selection bias could have been introduced, for example, by not analyzing regions with low wind speeds or mostly cloudy conditions or areas in regions with many emission sources and higher background levels.

Figure 2NO2 tropospheric vertical column of the level-2 S5P TROPOMI product from March 2018 to February 2020. Red circles mark NOx emission sources analyzed in this study. All analyzed source regions are listed, with more details in the Appendix, Table A1.

If the conditions are favorable (e.g., clear-sky, homogeneous, and strong wind conditions), even a single overpass can deliver valuable results for some regions (Lorente et al.2019; Goldberg et al.2019). To analyze the variability due to seasonal or weekday effects in a more general way, it is useful to have some statistics to prevent the introduction of a bias into the seasonal analysis by weekday variability. Due to the high spatial resolution and the good signal-to-noise ratio of the TROPOMI data, it is possible to use a data set of only 2 years' length to analyze successfully regions which are covered more frequently by clouds or which have low NO2 signals. Nevertheless, not every source region is used for every analysis in this study because of data availability after separating data into seasons or weekdays and weekend days. From the 50 targeted regions shown in Fig. 2, 8 are in the Southern Hemisphere. In order to integrate them into the analysis, they are mirrored in latitude and shifted by 6 months in season. All source regions are listed in the Appendix, Table A1. Surgut (Russia) is the NOx source with the highest latitude, at 61.25 N, and Singapore (Singapore) is the one closest to the Equator, at 1.3 N.

One example of the selected source regions is shown in Fig. 1. Figure 1a shows the average of NO2 tropospheric vertical columns for 2 years of data from 1 March 2018 to 29 February 2020 for days with wind speeds >2m s−1 in the region of the power plants Medupi and Matimba in South Africa. The NO2 distribution shows an isolated source and its plume, which has a high contrast to a low background concentration. In the southeastern part of the map, high tropospheric NO2 columns are visible. These originate from the South African Highveld conurbation near Johannesburg.

3.2 Rotation technique

To investigate the spatial pattern of NO2 column measurements, we combine the approach from Beirle et al. (2011) of a directional classification to determine the distribution, as a function of downwind distance, with a rotation of each TROPOMI measurement with its merged wind direction around the source to a common wind direction (Pommier et al.2013; Valin et al.2013).

Each TROPOMI observation is merged with the ERA5 wind data and rotated with its wind direction about a reference point (center of source, e.g., the city center or the power plant site) to a common wind direction, preventing a neutralization of outflow patterns of opposite wind directions. After the rotation, the TROPOMI measurement points are redistributed along the upwind–downwind direction, with the enhancement located in the downwind area of the source reference point and a clear outflow pattern. The measurement points maintain their upwind–downwind character and are analyzed simultaneously, independent of their wind direction. Figure 1b shows the mean of the rotated NOx tropospheric vertical column for the Medupi/Matimba example, where a pronounced plume with a clear upwind–downwind distribution is found. The black lines indicate a sector of ± 50 km around the source. Only data in this sector are used in the next step for the calculation of NOx emissions and lifetime. As an additional quality filter, only days where at least 50 % of the ground scenes in this area contain measurements and are not filtered because of clouds or wind speed are used in the analysis.

3.3 Line density calculation

The average outflow pattern of the NO2 tropospheric vertical column with a decay of the signal with distance from the reference source point reflects transport and nonlinear effects of atmospheric chemistry. Beirle et al. (2011) proposed a model to estimate the NOx emissions and lifetimes by integrating the mean NO2 columns in the across-wind direction and thereby reducing the two-dimensional maps to one-dimensional so-called line densities with unit molec. cm−1. In this study, we first converted the NO2 columns for each pixel into NOx columns to obtain NOx line densities from which NOx emissions are calculated directly (comparable to Beirle et al.2021) instead of applying the commonly used fixed [NOx] / [NO2] ratio of 1.32. Assuming that the Leighton photostationary state applies for the polluted air masses investigated, the NO2 is considered a surrogate for NOx, and concentrations of NO and NO2 are coupled by

(1) NO x NO 2 = 1 + NO NO 2 = 1 + J NO 2 k NO + O 3 n O 3 ,

with JNO2 the photolysis frequency of NO2 and kNO+O3 the rate constant for the reaction of NO with O3 (Seinfeld and Pandis2006).

Ozone data are taken from ERA5 and interpolated to the S5P overpass as described in Sect. 2.3.

For clear-sky conditions the photolysis frequency in the boundary layer is parameterized as a function of solar zenith angle (SZA) as proposed by Dickerson et al. (1982):

(2) J NO 2 = 0.0167 exp - 0.575 cos ( SZA ) ( s - 1 ) .

The rate constant kNO+O3 can in general be well represented by the Arrhenius expression, following the recommendation by Atkinson et al. (2004):

(3) k NO + O 3 ( T ) = 2.07 × 10 - 12 exp - 1400 T ( cm - 3 molec . - 1 s - 1 ) ,

with temperature T (K) taken from the ERA5 reanalysis hourly data with a horizontal resolution of 0.25× 0.25, averaged about the boundary layer, interpolated to TROPOMI overpass time, and oversampled to the same spatial resolution.

From the averaged NOx columns NOx line densities are calculated by integrating perpendicularly to the wind direction. To reduce the influence of possible surrounding sources, the line density is only calculated in a sector around the source (black vertical lines in Fig. 1b), which does not cut the plume and thereby misses emissions. Typical values used for the sector size are between ± 15 and ± 70 km across the plume, depending on its width, and up to 200 km upwind and up to 400 km downwind of the source, depending on plume length and the presence of other influencing sources. Figure 1c shows the calculated line density as a function of distance to the source. The calculated line density is shown in gray and the fit in black. Due to mixing and transport processes, the maximum is shifted in the wind direction, and the line density curve is steep upwind and less steep downwind with an exponential decay. From this line density curve, NOx lifetime and emissions can be estimated. The fitting model M as a function of distance to the source x is described by (similarly to Beirle et al.2011, Supplement)

(4) M ( x ) = E e G ( x ) + B ,

with a convolution of the exponential e and the Gaussian G function scaled by a multiplicative emission factor E and shifted by a background concentration offset B. The exponential function describes transport and chemical decay:

(5) e ( x ) = exp - ( x - X ) x 0 ,

with x>X (downwind) and else zero, where X is the location of the apparent source relative to the source reference point and x0 the distance over which the line density decreases by a factor of e (e-folding distance). The Gaussian function represents the broadening of the source by spatial smoothing with the Gaussian function width σ, which accounts for spatial smoothing caused by the extent of the spatial source, the TROPOMI pixel size, and wind variations:

(6) G ( x ) = 1 2 π σ exp - x 2 2 σ 2 .

This results in

(7) M ( x ) = E ( exp - ( x - X ) x 0 1 2 π σ exp - x 2 2 σ 2 ) + B = E 2 exp σ 2 2 x 0 2 - x - X x 0 erfc σ 2 - x 0 ( x - X ) 2 σ x 0 + B .

The fitted e-folding distance x0 and the mean wind speed w from the line density sector are then used to calculate the mean lifetime:

(8) τ = x 0 w .

The calculated lifetime includes effects of deposition, chemical conversion, and wind advection but must be considered an effective mean dispersion lifetime. This is because downwind changes, for example, due to a changing [NO2] / [NOx] ratio or the effects of nonlinearities on the NOx lifetime in the plume, are not considered in the method used.

The multiplicative emission factor E characterizes the total amount of NOx near the source and is used together with the mean wind speed w from the line density sector to derive the NOx flux in mol s−1:

(9) E NO x = E N A w ,

with the Avogadro constant NA=6.02214076×1023mol−1.

The described method to calculate NOx emissions and lifetimes is the EMG method. Uncertainties and error bars for NOx emission and lifetime estimates presented in the following results are based on 1 standard deviation derived by the EMG fitting procedure and are calculated with error propagation. The total uncertainty of NOx emissions and lifetimes is influenced additionally by many different error sources, especially the TROPOMI NO2 tropospheric vertical column itself, which are mainly systematic and lead to an overall low bias of the derived NOx emissions which is discussed in more detail in Sect. 4.6.

4 Results and discussion

The EMG method was applied to the mean TROPOMI NO2 column data of the selected regions, and NOx emissions and lifetimes were calculated for all available TROPOMI measurements of the 2-year period from 1 March 2018 to 29 February 2020 and also separated into seasons, working days and weekends, and pre-Covid-19 times and the Covid-19 pandemic. Not all source regions are used for all analyses due to sometimes poor statistics when separating the 2 years of data into specific periods.

4.1 Comparison to similar studies and emission inventories

The NOx emission estimates from this study are first compared to results from other recent studies (Beirle et al.2019; Goldberg et al.2019; Lorente et al.2019), which estimated NOx emissions with TROPOMI data. The studies used for comparison focused only on specific regions and used different periods, and the methods differ slightly but have in common that all used TROPOMI data for their emission calculations. Table 1 compares the NOx emission estimates for all source regions used in the comparative studies to the emission estimates retrieved in this study. For the study of Goldberg et al. (2019), which uses a time window from May to September 2018, it is possible to estimate and compare the results for the same period; values are given in parentheses. Due to different periods and slightly different methods and input parameters, differences are expected, but in general the emissions are in reasonable agreement.

Beirle et al. (2019)Goldberg et al. (2019)Lorente et al. (2019)

Table 1NOx emissions estimated for seven source regions are compared to the emissions derived from other studies using TROPOMI data. Values in parentheses are derived for the same data period as in Goldberg et al. (2019). Errors are 1 standard deviation derived by the EMG fitting procedure.

Download Print Version | Download XLSX

Beirle et al. (2019) used modified operational TROPOMI NO2 data for the period from December 2017 to October 2018, and NOx emission estimates are based on the continuity equation. For Riyadh NOx emissions are estimated for the urban area and a greater area of 250 × 250 km2 around Riyadh. Emissions retrieved in this study are slightly higher than those for Riyadh and significantly higher than those for Medupi/Matimba given in Beirle et al. (2019). Possible reasons for this are different wind fields and the conversion from NO2 to NOx. Beirle et al. (2019) used wind information from a fixed vertical level at about 450 m and a fixed [NOx] / [NO2] ratio of 1.32 compared to the wind data averaged over the boundary layer and the conversion based on the assumed photostationary state with ERA5 ozone data used in this study. For Riyadh the estimated mean [NOx] / [NO2] value is 1.41 and therefore slightly higher than the 1.32 used by Beirle et al. (2019). For Medupi/Matimba it is at 1.54 significantly higher and can explain in large part the higher NOx emissions estimated in this study compared to those in Beirle et al. (2019).

Goldberg et al. (2019) used operational and modified NO2 data from May to September 2018, from which here only the emissions based on the operational data product are used. They found that plumes can be misaligned by up to 30 after rotation based on the ERA5 wind data. Consequently the authors made a visual inspection of every day and manually adjusted the data. Due to the focus on summer months in Goldberg et al. (2019) and generally lower NOx emissions in summer than in winter (see also Sect. 4.2), lower values are expected in comparison to the 2-year average used in this study. For a better comparison we added estimates for the same period from May to September 2018, and, as expected, the agreement improves. Nevertheless, differences remain and are especially significant for New York; possibly the different wind information has a greater influence here.

Lorente et al. (2019) used 36 orbits of the operational NO2 product obtained on 34 different days between February and June 2018 to investigate the NOx emissions from Paris. Their NOx emissions of 53 mol s−1 agree well with the estimated emissions for Paris of 56.2 mol s−1 from this study.

The comparison indicates that the NOx emissions from the source regions are in agreement with previous estimates taking into account differences in the analysis techniques.

The retrieved NOx emissions for the 2 years of data have also been compared to those reported in the EDGAR emission database. In Fig. 3 the estimated NOx emissions for the 50 source regions are plotted versus the respective emissions received from the EDGAR database for the year 2015. Most NOx source regions have higher emissions in the EDGAR database than estimated by the EMG method in this study. One possible explanation for differences is the different periods the two methods are based on, with 2015 for the EDGAR database and March 2018 to February 2020 for this study. In addition, the EMG method is only sensitive to daytime emissions on nearly cloudless days close to the times of measurements, whereas the EDGAR database gives 24 h annual averages. Another possible explanation is the known low bias of current TROPOMI tropospheric NO2 columns as compared to ground-based or aircraft measurements of the tropospheric NO2 column. This underestimation is more pronounced for regions with larger NO2 columns (Verhoelst et al.2021), which is in agreement with our finding that differences from EDGAR are largest for the source regions with high NOx emissions such as Tokyo, Wuhan, Moscow, and New York. This underestimation of TROPOMI will probably be reduced using the reprocessed data set with V02.02 and is discussed in more detail in Sect. 4.6. The estimated NOx emissions for the Medupi and Matimba power plants are an exception as they show significantly higher emissions estimated by the EMG method than reported in the EDGAR database. One likely explanation for this observation is that the Medupi power plant was put in operation only in 2015, the year to which the EDGAR database refers. Consequently, it was not or only partially considered in EDGAR. This also shows the potential of the EMG method to identify NOx point sources that are not included in emission inventories.

Figure 3NOx emissions derived from 2 years of TROPOMI data (1 March 2018 to 29 February 2020) for 50 sources calculated with the EMG method compared to the NOx emissions derived from the EDGAR emission database (v.5.0, 2015). The dashed line shows the 1:1 ratio. Error bars are the standard deviation derived by the EMG fitting procedure for emission estimates.


Emissions from power plants in the United States are measured continuously at the stacks and are made available as hourly mean values by the Environmental Protection Agency's (EPA's) Continuous Emission Monitoring System (CEMS). The NOx emissions for the Colstrip power plant (Montana, USA) are reported in the CEMS database and can also be analyzed with the EMG method based on TROPOMI data. Figure 4a shows a time series of the NOx emissions reported by the CEMS for the Colstrip power plant from February 2018 to March 2020 (blue) and reported data filtered for days with TROPOMI measurements that were used for the EMG method (green). Assuming that a typical NOx lifetime is around 3 h on average for the Colstrip power plant, reported NOx emissions are averaged between 11:00 and 14:00 local time, approximately 3 h before a typical TROPOMI overpass. The reported NOx emissions for the Colstrip power plant show a high variability over the course of the year.

The EMG method is a widely used method but is mostly used for data around summer. In this study we will also extensively use winter data which bring greater uncertainties to the analysis, for example, due to longer lifetimes in winter, and additionally data availability is often worse during winter than summer. Therefore, we wanted to verify whether the EMG method can reproduce the variability seen in the reported emissions independently of season. To analyze this variability, we defined seven periods of a range between 2 and 5 months during February 2018 to March 2020 (dashed black lines). CEMS data reported on TROPOMI measurement days are averaged for these periods and are shown as horizontal green lines with standard deviation. Estimates from the EMG method for these periods are shown as horizontal red lines with standard deviation derived by the EMG fitting. The averaged CEMS data and EMG estimates show a similar pattern of NOx emissions over the year with an underestimation of the CEMS data by the TROPOMI EMG results for all periods. This underestimation could be explained by a loss of NOx after its measurement in the power plant stack or in the early stage of the NO2 plume formation and to a large extent by the discussed known low bias of the TROPOMI NO2 data. The low bias is expected to be particularly pronounced for power plants due to wrong AMFs from the TM5 model vertical profiles (Beirle et al.2019, Supplement). Figure 4b shows a comparison of the two data sets for the seven periods in a scatter plot. The 1:1 line is indicated with the gray dashed line. The solid black line indicates the linear regression, with a slope of 0.4 ± 0.09 and a correlation of 0.9, reflecting the underestimation of the CEMS data by the TROPOMI EMG results but also the well-matching emission patterns.

Figure 4(a) CEMS NOx emissions for the Colstrip (Montana, USA) power plant compared to the NOx emissions estimated by the EMG method, using TROPOMI NO2 data. CEMS data are daily mean values, averaged between 11:00 and 14:00 local time, shown in blue and filtered for days with TROPOMI measurements shown in green. Means over seven selected periods, marked by the dashed black lines, are shown by the horizontal green lines for the CEMS data and red lines for the TROPOMI EMG estimates. (b) Scatter plot of CEMS-measured NOx emissions versus TROPOMI EMG NOx emission estimates for the seven periods. The 1:1 line is indicated with the gray dashed line. The solid black line indicates the linear regression. Error bars are the standard deviation derived by the EMG fitting procedure for the EMG emission estimates and the monthly standard deviation for the CEMS mean values.


From this comparison we can conclude that power plant emissions can vary strongly over time and that the TROPOMI-based NOx emissions can be lower by about a factor of 2 for power plants as already discussed in previous publications (Beirle et al.2019). However, from this comparison, we can conclude that the EMG method applied to TROPOMI data, at least for this power plant, reproduces the temporal variability seen in the CEMS data reasonably well and does not show a clear seasonal bias. Therefore, it can be assumed that the method also gives reasonable results when analyzing winter data.

4.2 Seasonality of NOx emissions

To investigate the seasonal variation of NOx emissions, the 2 years of TROPOMI data are separated into four seasons, and emissions are calculated separately for each season. Winter months are defined as December to February (DJF, Southern Hemisphere JJA), spring as March to May (MAM, Southern Hemisphere SON), summer as June to August (JJA, Southern Hemisphere DJF), and fall as September to November (SON, Southern Hemisphere MAM). The 2-year data set provides for each season a maximum of 2 times 3 months. Due to cloud cover, inhomogeneous wind patterns, and resulting partly poor statistics after separating into seasons, only 34 of the 50 source regions are analyzed. For Wuhan, which is already affected by Covid-19 regulations when analyzing data from March 2018 to February 2020, the analysis is limited to data from March 2018 to 22 January 2020.

In Fig. 5 the summer-to-winter ratios of the NOx emissions are plotted for these 34 source regions. For most of the source regions, it was found that the NOx emissions are higher during winter months than during summer. Naples (Italy), Wuhan (China), and Novosibirsk and Krasnoyarsk (Russia) are among the cities with the lowest summer-to-winter ratio of around 0.3 and accordingly clearly higher NOx emissions in winter than during summer. This is expected in places where domestic heating in winter contributes significantly to NOx emissions. For some source regions the found ratios are unexpected, for example, for Casablanca (Morocco), with a summer-to-winter ratio of 0.32 ± 0.17. Due to the location on the Atlantic coast, which mitigates temperature fluctuations, summers are hot but typically not very hot and winters mild, and smaller seasonal variation in emissions is expected. No explanation for this observation could be found so far. For some regions, the NOx emissions are actually higher in summer than in winter, for example, Medupi/Matimba in South Africa. Interestingly, the regions with a summer-to-winter ratio larger than 1 or close to 1 in Saudi Arabia (Riyadh, Rabigh, Tabuk, Buraidah), Libya (Sarir Field), and Sudan (Khartoum) are all located in hot desert climate regions which are dry and warm all year round but where, especially in summer, long heat periods are common. A likely explanation is therefore that the NOx emissions are higher in summer than in winter due to higher power consumption caused by the use of air conditioning in summer.

Figure 5Summer-to-winter ratio of retrieved NOx emission data (1 March 2018 to 29 February 2020). Data are separated into winter (Northern Hemisphere: DJF, Southern Hemisphere: JJA) and summer (Northern Hemisphere: JJA, Southern Hemisphere: DJF) for 34 source regions from left to right in increasing ratio. The red line indicates the 1:1 line where summer and winter NOx emissions are equal; below the line winter emissions are larger, and over the line the summer emissions are predominant compared to the winter emissions. Error bars are the standard deviation derived by the EMG fitting procedure for emission estimates.


Beirle et al. (2011) reported that the derived NOx emissions are stable throughout the year for Singapore and Madrid and that variability found for other cities, which could reflect a change in emissions, was not significant. Most other emission studies focused on summer months only, and therefore an investigation of the seasonality of NOx emissions was not possible (Beirle et al.2019; Goldberg et al.2019; Ialongo et al.2014). Lorente et al. (2019) found the highest NOx emissions to occur on cold weekdays in February and the lowest emissions on warm weekend days in spring 2018. Combined with comparisons to emission inventories, the authors concluded that this indicated the importance of the contributions from the residential heating sector to NOx emissions in winter. Using the full 2-year data set, we found a summer-to-winter ratio of 0.51 ± 0.2 for Paris, which supports the results of Lorente et al. (2019).

Overall, the investigation of seasonality shows higher NOx emissions in winter than in summer for the majority of source regions, which is expected for the location of the majority of the selected regions at mid latitudes and in temperate climate. Not all ratios are significantly different from unity, and some results are unexpected, but the trend from the source regions at high latitudes with higher NOx emissions in winter to the regions in desert climate with higher emissions in summer is plausible. In some part our results differ from those previously reported in the literature (Crippa et al.2020; Zheng et al.2018), where generally weaker seasonality is found. This may reflect the use of emission inventories and the analysis of average values for countries in the literature that probably does not reflect the specific situation in cities. The composition of sources for a specific city is expected to be different from that for a country. Further detailed investigations and comparisons are needed to understand the discrepancies. Nevertheless, it should be considered that the uncertainties of the EMG NOx emission estimates are relatively large, and especially the analysis of winter data has not been performed often and is influenced even more by uncertainties (see Sect. 4.1 and 4.6).

Figure 6(a) The seasonal NOx line densities for the city of Madrid. The line densities are in light colors and the fit in dark colors. (b) Estimated lifetimes as a function of season for 10 source regions. (c) Estimated NOx lifetimes for the selected 50 sources as a function of latitude. (d) Estimated lifetimes in dependence on latitude (black) separated into winter (blue), spring (green), summer (red), and fall (yellow). All plots are based on TROPOMI NO2 data from 1 March 2018 to 29 February 2020. Sources from the Southern Hemisphere are mirrored in latitude and season. Given uncertainties and error bars are the standard deviation derived by the EMG fitting procedure.


4.3 Latitudinal and seasonal dependence of lifetimes

Another parameter retrieved by the EMG method is the mean effective lifetime of NOx. Thus, in addition to the seasonality of NOx emissions, the seasonality and latitudinal dependence of lifetimes have been investigated. When data averaged over the 2 years are considered, lifetimes can be calculated for all 50 source regions. Separating the data set into seasons reduces the number of available data, and lifetimes can be estimated for 34 of the 50 source regions in total.

Figure 6a shows the line density in dependence on distance to the source, separated into four seasons for the city of Madrid. The light blue line is the calculated line density and the dark blue one the fitted curve for the winter months (DJF). For the two winters (6 months) between 1 March 2018 and 29 February 2020, 25 d could be used for the analysis with a mean wind speed on these days of 5.5 m s−1. For winter, this results in NOx emissions of 84.2 ± 8.8 mol s−1 and a lifetime of 2.41 ± 0.27 h. The results for spring (MAM) are shown in green, the summer results (JJA) in red, and the fall months (SON) in yellow. The calculated lifetimes vary between 1.59 ± 0.16 h in summer and 2.41 ± 0.27 h in winter. The average lifetime over the full 2-year period is 2.03 ± 0.15 h.

In Fig. 6b estimated NOx lifetimes for 10 source regions are shown as a function of season. Only a selection of the sources is shown here, which are a mix of regions from southern and northern latitudes, distributed over many latitudes and with different source types. In general these source regions show similar behavior, with the longest lifetimes in winter, the shortest ones in summer, and quite similar lifetimes in spring and fall that are between the summer and winter estimates. Also, deviations are visible, as could already be seen for the case of Madrid, where the fall lifetime is much shorter than in spring. Seasonal differences are strongest for the higher latitudes and lowest for the sources closest to the Equator. In general, the results are in agreement with values shown in Beirle et al. (2011), where eight different source regions were analyzed, yielding lifetimes within a range of 2 to 6 h and a maximum of 8.5 h during wintertime for Moscow.

The averaged lifetimes over the full 2-year period for these 50 source regions as a function of latitude are shown in Fig. 6c. As a result of the lower Sun and thus reduced photolysis and likely lower OH at higher latitudes, one would expect an increase in NOx lifetime with latitude (Martin et al.2003; Stavrakou et al.2013). This is broadly confirmed by the results, averaged lifetimes increasing from approximately 2 h for source regions at low latitudes near the Equator to about 6 to 8 h for source regions at high latitudes of around 60.

The retrieved lifetimes can be used to estimate atmospheric OH concentrations. Assuming that the decay of NOx is determined by the termolecular reaction of OH with NO2 with the rate constant kOH+NO2+M for 298 K and 1 atm reported by Burkholder et al. (2020), the mean OH concentration ranges from 0.3×107 to 1× for a range of NOx lifetimes of 2 to 8 h. This is in a reasonable range for OH concentrations at midday (Holland et al.2003; Smith et al.2006; Lu et al.2013).

Figure 6d also shows, in addition to the 2-year mean, the lifetimes separated by seasons. As a result of the reduced number of data points, not all source regions (34 out of 50) could be included. All seasons show similar latitudinal dependence as described for the mean lifetime, but in addition a seasonal dependence is visible, with the shortest lifetimes in summer and longer lifetimes in winter.

We can only detect a weak seasonality of the NOx lifetime, which is in agreement with Beirle et al. (2011) but significantly lower than would be expected from the analyses with the chemical transport model GEOS-Chem by Martin et al. (2003) and Shah et al. (2020), which showed lifetimes of about 1 d in winter compared to around 6 h in summer. Possible explanations could be as yet insufficient statistics, a clear-sky bias, and also the midday observation time of TROPOMI, which all lead to more balanced lifetimes in summer and winter. Furthermore, it should be considered that the total uncertainty for lifetime estimates is larger than just the standard deviation derived by the EMG fit, especially in winter due to longer lifetimes and often reduced data availability (Sect. 4.6). The large data set used in this study reveals a clear but complex latitudinal dependence of NOx lifetimes. This can be used for studies similar to those of Beirle et al. (2019, 2021), where an assumption about the lifetime is necessary to calculate NOx emissions and can also provide relevant observational constraints for model simulations of NOx lifetimes.

4.4 Weekend effect

The presence of a weekend effect in NO2 on a global scale was first shown by Beirle et al. (2003) with GOME measurements. Anthropogenic activities have their maximum during the week and are reduced during the weekend. Thus, we expect fewer NOx emissions in cities on weekends. This behavior should be reflected in a comparison of NOx emissions on weekdays and weekends. How large the difference between weekends and weekdays is depends on the types of NOx sources and the different patterns of anthropogenic activity in the source region. To investigate the weekly cycle, the TROPOMI data are separated into weekdays and weekend days, and NOx emissions and lifetimes were calculated separately. Weekend days can be one or two days, and those days can also differ according to religious tradition. For source regions, for example, in Europe and the United States, weekend days were set to be Saturday and Sunday, and for Saudi Arabia weekend days are Friday and Saturday (see also Table A1). Often one weekend day can be different from the other; for example, emissions on Saturdays are often not as low as those on Sundays, even if they are already lower than those from Monday to Friday (Crippa et al.2020; Goldberg et al.2021). If the weekend would be limited to one day in the analysis, this would lead to a strengthening of the weekend-to-week ratio for some cities.

Figure 7 shows the weekend-to-weekday ratio, with higher NOx emissions during weekdays than on the weekend for most of the source regions albeit with rather high variability. For Paris (France), the NOx emissions are reduced by 40 % on the weekend, which agrees with Lorente et al. (2019). Chicago only shows a reduction of 16 % on the weekend versus weekdays, which is less than the weekend reductions of 30 % found in Goldberg et al. (2019, 2021). This could be attributed to the fact that in Goldberg et al. (2019, 2021) Sunday was assumed to be the only weekend day. This results in a more pronounced weekend reduction than with Saturday and Sunday as weekend days, as done in this study.

Source regions not showing any reductions on the weekend are mostly dominated by power plants. Rabigh in Saudi Arabia is a small city with a large gas-fired power plant. Medupi and Matimba are two large coal-fired power plants in South Africa. Sarir Field is a large oil field in Libya. Hwange in Zimbabwe is a small town with Zimbabwe's biggest power plant. Chelyabinsk (Russia) is a city with nearly 1.2 million inhabitants but also a center of heavy industry, and Hassi Messaud in Algeria is an oil refinery town. Source regions with less industry which are dominated by domestic and transport emissions, like Tokyo, Brisbane, or Madrid, on the other hand, show large emission reductions on the weekend of 50 % to 60 %.

Despite showing large NO2 columns, Chinese cities did not show a weekend effect in previous studies (Beirle et al.2003; Stavrakou et al.2020). Stavrakou et al. (2020) showed an average weekend-to-weekday ratio in the NO2 column over all large Chinese cities of 0.97 ± 0.02 using 2005–2017 OMI NO2 columns and 1 year of TROPOMI NO2 columns (May 2018–April 2019). Ratios of NO2 columns are related but not identical to ratios in NOx emissions as discussed here, as effects from meteorology and lifetime are not accounted for, and tropospheric background is contained, which is removed in the EMG method by fitting the background. Using the EMG method to investigate the weekly cycle, we calculated a weekend-to-weekday ratio of 0.79 ± 0.3 in NOx emissions for Wuhan (China). The ratio includes 27 weekend days and 67 weekdays during the period March 2018–22 January 2020. Due to the small number of isolated point sources and because of the limited statistics from using 2 years of data, it was not yet possible to determine weekly cycles for other Chinese cities. The ratio for Wuhan indicates a reduction of NOx emissions on rest days also in China in recent years, in contrast to earlier studies showing no such effect. This can be analyzed in more detail as more TROPOMI data become available.

Figure 7Weekend-to-weekday ratio of NOx emission data (1 March 2018 to 29 February 2020) from left to right in increasing ratio. Data are separated into working and weekend days for 44 source regions. Weekend can be one or two days, and those days can also differ according to religious tradition, which is considered. The red line indicates the 1:1 line where weekend and week NOx emissions are equal; below the line emissions during the week are predominant, and above the line weekend emissions are higher. Error bars are the standard deviation derived by the EMG fitting procedure for emission estimates.


A side-effect of the reduced NOx emissions on weekends might be lower OH levels due to the slower production rate of OH from the NO+HO2 reaction and therefore longer NOx lifetimes (Stavrakou et al.2008). However, in our data set the retrieved lifetimes for weekdays and weekend days do not show a clear enhancement of lifetimes on weekends (see Fig. A1). This may possibly be due to an insufficient number of data points required for a statistically significant result, and this question should be revisited once a larger TROPOMI data set becomes available.

4.5 Covid-19 effect

In early 2020, several countries took containment measures to prevent the spread of the coronavirus outbreak (Covid-19), which caused reductions in industrial activities and traffic volume. Due to the link of NOx emissions to human activities, it is possible to investigate the impact of the Covid-19-induced activity reductions with the TROPOMI NO2 column data (Bauwens et al.2020; Liu et al.2020; Goldberg et al.2020). Part of the observed decreases in the NO2 columns may result from effects other than the measures designed to limit transmission of Covid-19, e.g., meteorological variability, seasonal variability, or environmental policy regulations. Consequently, it is problematic to identify a clear Covid-19 effect using only NO2 column amounts. Besides the possibility of using models to separate the Covid-19 effect from other factors (Goldberg et al.2020), it is also possible to use the EMG method. This approach accounts for wind conditions and NOx lifetime, which influence the NO2 columns observed for similar NOx emissions. Since the EMG method can only be used to investigate point sources, the range of possible study areas is limited. In addition, the analysis of the Covid-19 effect is also much more limited compared to the rest of the study, because we decided to compare only monthly means from 2 different years to each other to minimize influence by meteorology. We focused on the cities of Buenos Aires (Argentina), New Delhi (India), and Madrid (Spain), which are considered to be point sources and have an appropriate number of days with satellite data available during the comparative periods. The EMG method was used, and monthly means of NOx emissions from 2019 and 2020 were calculated and compared. If only a few days per month are available, fit results can be of low quality if conditions on these days are not ideal and statistics are missing, so that, for example, the variation of NOx emissions on weekdays can also have a significant influence. Consequently, even for the selected cities not all months can be considered.

Figure 8 shows the monthly means from January to November of the calculated NOx emissions of TROPOMI data for 2019 and 2020 for (a) Buenos Aires, (b) New Delhi, and (c) Madrid. The same months in 2019 and 2020 are compared and thus pre-Covid-19 periods with periods of containment measures. For analyses where NOx emissions from different months are compared, seasonality of NOx emissions must be considered. The NOx emissions retrieved for Buenos Aires are low during summer from January to March (months 1–3), increase in the winter months with a maximum in July, and then decrease again towards summer. The emissions for New Delhi do not show such a strong seasonality as those for Buenos Aires but in general are also higher during the winter months January and February (months 1–2). The seasonality is also clearly visible for Madrid, which must be considered in possible comparisons.

On 20 March 2020, a strong nationwide lockdown began in Argentina and lasted in Buenos Aires for more than 7 months, ending on 8 November. From January to March, NOx emissions from 2020 are comparable to those from the same months in 2019. In April 2020, the first complete month in lockdown, the emissions are 57 % lower than in April 2019, and in May 2020, the emissions are also 36 % reduced compared to 2019. In June, however, emissions are higher in 2020 than in 2019, although there has been no major change in the containment measures by the government. A possible explanation is that the June 2019 emissions are lower than expected from the seasonal cycle compared to May and July 2019. It is also possible that June 2020 emissions are unexpectedly high due to a cold winter month and additional emissions from heating. This is also indicated by ERA5 reanalysis temperature data. They were averaged over the boundary layer in the Buenos Aires target area for the measurement days used in our study in May, June, and July 2019 and 2020. They show that June 2020 was on average 3 C colder than 2019. In comparison, temperatures for May and July 2019 and 2020 are very comparable. The NOx emissions in July behave in a manner similar to those in March and April but have a somewhat smaller decrease of 26 % compared to 2019. These reductions are similar to the reductions of 31 % and 44 % in August and September. In October the NOx emissions are almost equal for both years, and in November 2020 they are higher than in 2019, possibly due to the end of lockdown on 8 November.

Figure 8Monthly NOx emissions calculated with the EMG method using TROPOMI NO2 tropospheric column data for 2019 (blue) and 2020 (orange) from January to November for (a) Buenos Aires (Argentina), (b) New Delhi (India), and (c) Madrid (Spain). The numbers in the bars represent the number of available days for the monthly mean. Due to an insufficient data availability of less than 3 d, comparisons are not possible for some months. Error bars are the standard deviation derived by the EMG fitting procedure for emission estimates.


In India, a nationwide strict lockdown started on 24 March 2020. In January the estimated NOx emissions are almost equal for both years. The significantly higher emissions in February 2020 compared to 2019 probably cannot be explained by a Covid-19 effect. They are more likely a result of the limited number of clear-sky days for TROPOMI observations, of on average only 5 d during February 2019, which will reveal more natural variability. In April 2020, the first complete month in lockdown, the emissions are 87 % lower than in April 2019, and also in May and June the emissions are 54 % and 31 % lower than in 2019. For July to September, comparisons are not possible, as less than 3 d of measurements per month are available due to cloud coverage, and statistical and fit results are not good. In October and November 2020, emissions are almost equal to 2019 levels, suggesting that NOx emissions are similar to those before the Covid-19 containment measures.

In Europe, Madrid was one of the cities strongly impacted by measures to prevent the spread of Covid-19. A strict lockdown was enacted in the middle of March 2020, some restrictions were lifted from the middle of April, and in May the government followed a plan for easing lockdown restrictions slowly. Due to again a rising number of cases, new restrictions started in early October 2020. Emission comparisons are not possible for January, April, and November due to a lack of data caused by cloud coverage. The significantly higher emissions in February 2019 compared to 2020 are probably mainly caused by two factors, first a strong synoptic meteorological variability in Europe and, second, February is also a month with typically persistent cloud cover, resulting in a reduced number of TROPOMI observations which will reveal more natural variability. Similar findings are also described in Bauwens et al. (2020) and Levelt et al. (2021). For the months May and October 2020, when Covid-19 restrictions were in place, the calculated NOx emissions in Madrid are 43 % and 63 % lower than in 2019. This is comparable to results from NO2 column comparisons which showed a reduction of around 30 % from the middle of March to early April 2020 relative to 2019 (Bauwens et al., 2020). The larger reductions in October may be due to stricter restrictions resulting in larger emission reductions. However, even with monthly averages, the high variability of NOx emissions is still not negligible and could also have an influence here.

Despite the shortness of the periods available for analysis, it is possible to investigate short-term variability of NOx emissions induced by Covid-19 with TROPOMI NO2 data and the EMG method. Strong decreases due to lockdown measures of 57 % in Buenos Aires and 87 % in New Delhi are shown for April 2020 compared to April 2019 as well as a general tendency towards lower emissions in 2020 after the start of the Covid-19 pandemic than in 2019. Nevertheless, even with monthly averages, the high variability of NOx must still be considered for particular months. These emission estimates account for wind conditions and NOx lifetimes and can therefore give a better estimate of the Covid-19 measures on NOx emissions than just comparing NO2 column measurements. For some cities and months, the number of days in monthly means is limited due to cloud cover. This is also a problem when comparing monthly NO2 column levels.

4.6 Uncertainties

The uncertainties and error bars for emission and lifetime estimates given in this work are based on the standard error (1 standard deviation) derived by the fitting procedure and are calculated with error propagation. For emissions this results in

(10) σ E = w σ E' N A + E σ w N A

and for lifetimes in

(11) σ τ = σ x 0 w + x 0 σ w w 2 ,

with the emission factor E, the e-folding distance x0, the mean wind speed w, Avogadro's constant NA (see Sect. 1), and the standard deviations σE of E and σx0 of x0 derived from the fit and σw derived from the wind field in the line density sector. These estimates are based on the fitting uncertainties.

However, the NOx emissions and lifetimes derived from TROPOMI NO2 data are influenced by additional error sources. The most important contribution directly influencing our estimates is the accuracy of the TROPOMI NO2 tropospheric vertical column itself. This uncertainty is dominated by the accuracy of the tropospheric AMF and is estimated to be on the order of 30 % (Boersma et al., 2004; Bucsela et al., 2013). Recent studies comparing the TROPOMI NO2 column with co-located ground-based or aircraft measurements reported a low bias for TROPOMI NO2 columns, which is most likely caused by a priori information such as the surface albedo, cloud-top height, cloud fraction, and the NO2 vertical profile used for tropospheric AMF calculations. This bias differs for different regions and is more pronounced for regions with larger NO2 columns (Griffin et al.2019; Ialongo et al.2020; Judd et al.2020; Dimitropoulou et al.2020; Verhoelst et al.2021). Some studies scaled up the measured NO2 columns with a factor of 1.33 for Paris (Lorente et al.2019) up to a factor of 1.98 for Germany (Beirle et al.2019). As this suspected underestimation is not yet fully characterized, and as it is not clear without independent measurements how much the various regions used in our study are affected, we decided not to correct the operational product. Therefore, NOx emissions derived in this study are systematically biased low. A new operational NO2 product version, V01.04, activated after the analysis period of this study on 2 December 2020, implemented the first major changes, which led to a substantial increase in the tropospheric NO2 column over polluted areas for scenes with small cloud fractions. Further major updates in the TROPOMI NO2 retrieval algorithm are done in V02.02 released on 1 July 2021. Tropospheric vertical columns are between 10 % and 40 % larger than the v1.x data depending on the level of pollution and season; the largest impact is found in wintertime at mid and high latitudes (van Geffen et al.2021). A complete mission reprocessing will be performed to get a harmonized data set (Eskes and Eichmann2021). The use of the reprocessed data set for this study will increase NOx emissions and probably affect the seasonality analysis but not so much the weekend-to-week comparison results.

To calculate NOx emissions, we applied a conversion of each TROPOMI pixel from NO2 column to NOx column by assuming that the Leighton photostationary state applies for the polluted air masses. This is more accurate than using a fixed conversion factor as in many earlier studies, especially for our analysis over a large latitudinal range and for different seasons. Nevertheless, the photolysis frequencies have to be parameterized and the temperatures for the rate constant and the ozone concentrations taken from the monthly ERA5 data set interpolated to the S5P overpass time. Thus, the conversion from NO2 to NOx adds systematic errors in the emission estimates.

In addition, using satellite data introduces a systematic clear-sky bias, because only measurements from nearly cloudless days are used, which favor specific emission patterns. These may differ from those of cloudier and therefore often cooler days. The limitation to nearly cloudless measurements also influences lifetime estimates, which are systematically lower due to higher photolysis rates on cloudless days. Furthermore, the retrieved NOx emissions and lifetimes are based on measurements in the early afternoon and are therefore systematically biased due to the measurement time. The variability in time can be further analyzed using follow-up sensors on geostationary satellites, for example, GEMS, Tempo, or Sentinel-4.

The NO2 tropospheric columns are strongly affected by the wind fields. This affects the calculation of NOx emissions and lifetimes. We filtered the NO2 measurements depending on the corresponding wind speed, and only data with wind speeds >2m s−1 are included in the analysis. As a result of the short lifetime of NO2, the observed NO2 distribution should, in general, be dominated by the wind conditions around the satellite overpass. On days with rapidly changing wind directions around the time of measurement, the spatial patterns and thus the estimates of emissions and lifetimes may be affected. An effect observed for some locations on some days is curved plumes, which for the case of strong curvature leads to a part of the plume being outside of the line density calculation sector and an underestimation of both NOx emissions and lifetime. This has a large effect when analyzing estimates for single days and is generally larger for longer lifetimes, i.e., especially on winter days. It has a smaller influence on the overall result analyzing a larger average if not too many days are affected by rapidly changing wind directions. Nevertheless, the wind field is the largest uncertainty (random and systematic) influencing our estimates after the NO2 column itself. Lorente et al. (2019) modified wind speeds by 20 % and found that emissions changed by 20 %, which demonstrates the strong influence of wind speed on NOx emissions. However, reliable global wind information is hard to obtain. Beirle et al. (2011) estimated an uncertainty of 30 % for the wind data. The uncertainties due to the chosen wind fields will vary for different source regions. Overall, we assume an uncertainty of 30 %.

To avoid interference from sources of NOx surrounding the target region, only rather isolated source regions are chosen for the analysis. Since almost no site is perfectly isolated, sectors were defined in which the line density was calculated by integration and fitted with the EMG method to minimize interference between different sources. Due to the rotation of measurements with their corresponding wind direction around the source, the NO2 signal from sources in the surroundings is smeared around the source region in their distance to the source location. To exclude these contaminants, the sector size used for the EMG method has to be chosen carefully. In order to have an adequate amount of data for a robust EMG fit, the sector must first be large enough in both the downwind and across-wind directions, but the size is also influenced by other factors. The sector length in the wind direction is mainly determined by the influence of other sources but also the spatial extent of the source region itself and wind speeds. A typical total size is 300 km with 100 km upwind and 200 km downwind of the source and is adjusted visually, if necessary, by inspecting the NO2 distribution and line density. If the influence of surrounding sources is negligible or becomes negligible by adjusting the sector size, the EMG method is robust in variation of the sector size in the wind direction. The sector width in the across-wind direction is mainly influenced by the geographical extent of the source region. If the sector width is chosen too small, some of the NOx emissions are outside of the sector due to dilution by wind or due to curved plumes, as described above. This obviously leads to an underestimation of the calculated emissions which acts as an additional apparent loss, leading to the e-folding distance x0 being systematically biased low, and lifetimes, defined as τ=x0w-1, with mean wind speed w, are also underestimated. Typical sector widths vary between 30 and 140 km and are determined by visually inspecting the NO2 distribution after rotation. Beirle et al. (2019) estimated the uncertainty of NOx emissions and lifetimes due to sector size as 10 %.

The EMG method is well suited to investigating point sources. In reality, most of the analyzed sources deviate from the assumption of a point source. Isolated sources, such as the Colstrip power plant in Montana, USA, or the Sarir Oil Field in Libya are close to being point sources. In previous studies, the Four Corners and San Juan power plants in New Mexico, USA, which are located 13 km apart, were investigated as one source using the EMG method (Beirle et al.2011; Goldberg et al.2019). With the higher resolution of TROPOMI compared to OMI, it becomes clear that the situation with the two plumes is more complex. This is evidenced by strong irregularities in the line density and fit, indicating that it should not be treated as one point source. In this particular case, re-gridding the TROPOMI data to a coarser resolution would potentially allow an analysis. Other sources such as Tokyo, Moscow, or Chicago, all cities with emissions originating from a larger area, are treated as extended point sources, but additional uncertainties must be considered. For example, it is not possible to consider a change in the instantaneous NOx lifetime downwind of the source with the EMG method, and estimated lifetimes should be interpreted as an effective mean lifetime (Beirle et al.2011). This effect is particularly pronounced in spatially extended source areas and can lead to low-biased lifetimes.

In summary, the total uncertainty of NOx emissions and lifetimes derived from TROPOMI NO2 data is influenced, in addition to the 1-sigma uncertainties, derived by the fitting procedure, by many different error sources, which are mainly systematic and lead to an overall low bias of the derived emissions dominated by the systematic errors in the TROPOMI NO2 tropospheric column itself (30 %–50 %) and in the wind field (30 %). The sector size can lead to low-biased emissions and lifetimes and is estimated to introduce an uncertainty of 10 %. Further uncertainties due to the measurement time, the clear-sky bias, and assumptions about point sources tend to have low bias in the lifetimes and emissions. More analysis is needed to resolve these issues. The error contributions summed in quadrature result in an overall uncertainty estimate for lifetimes and emissions in the range 43 %–62 %.

5 Conclusions

In this study, we present investigations of the variability of NOx emissions and lifetimes estimated from S5P TROPOMI observations of selected NOx source areas around the world. Similarly to earlier studies, we combine TROPOMI NO2 tropospheric vertical column data with wind information from ECMWF ERA5 reanalysis for the EMG method. TROPOMI measurements with their high spatial resolution and high signal-to-noise ratio enable the detailed analysis of NOx emissions and lifetimes using only 2 years of data. Investigating small emission sources not analyzed before and also monitoring the variability on a short-term temporal basis has been possible. Here, a total of 50 NOx sources from different regions, located between the Equator and 61 latitude, have been investigated. Despite efforts to select a broad range of regions well distributed around the globe and a mix of different sources, predominantly cities, isolated power plants, industrial regions, oil fields, and regions with a mix of sources, some selection bias could have been introduced, for example, by not analyzing regions with low wind speeds or mostly cloudy conditions or areas in regions with many emission sources and higher background NO2 levels.

Comparisons to other studies using TROPOMI data and similar methods for emission estimates show small differences and in general good agreement taking into account the differences of the analyses. The emission estimates are compared to the EDGAR (v.5.0 2015) emission inventory. For most of the source regions the emissions are higher than those estimated with the EMG method. These differences are particularly strong for the source regions with the highest emissions. As the operational TROPOMI tropospheric NO2 product has been reported to be low in comparison to independent measurements, part of the apparent overestimation by EDGAR could be related to a TROPOMI low bias. On the other hand, NOx emission trends over the last 5 years are likely to also have an impact. EPA CEMS data for the Colstrip power plant show that power plant emissions can vary strongly over time. The EMG method applied to TROPOMI data, at least for this power plant, reproduces the temporal variability seen in the CEMS data reasonably well and does not show a significant seasonal bias. Consequently, we consider that the method also gives reasonable results analyzing winter data. However, the TROPOMI-based NOx emissions are lower by about a factor of 2 for this power plant as compared to the CEMS data, in agreement with earlier findings for other power plant sources.

The seasonal separation of the emission estimates indicates that the highest emissions occur during winter, especially for source regions at higher latitudes, whereas in hot desert climates the highest emissions occur in summer. This is best explained by the different contributions to NOx emissions depending on source region, which are typically dominated by domestic heating in winter or air conditioning in hot summer months, depending on the climatic conditions of the source region. In some respects our results differ from those in previous studies using emission inventories. These studies have analyzed emissions for countries, where generally weaker seasonality is found. However, this probably does not reflect the specific situation in cities, where the composition of sources is expected to be different. Further detailed investigations and comparisons would be helpful to understand the discrepancies between emission estimates using emission inventories and the EMG method with satellite data. Generally, it should be considered that the uncertainties for emission estimates from the EMG method can be large, and especially the analysis of winter data has not been performed often and is influenced even more by uncertainties.

The investigation of the seasonal and latitudinal dependence of the NOx lifetime shows an increase in lifetime from 2 to 6 h for increasing latitudes but only a weak seasonal dependence with longer lifetimes in winter than in summer. Larger differences are expected based on modeling studies. The weak seasonal dependence found in this study could be explained by as yet insufficient statistics, a clear-sky bias, and also the midday observation time of TROPOMI, which all lead to more balanced lifetimes in summer and winter. Uncertainties which tend to have a low bias in lifetime estimates, especially in winter due to longer lifetimes and the aforementioned points, should be considered.

Separating NOx emission estimates into working and weekend days indicates that, for most NOx source regions, emissions are higher during weekdays than on the weekend, albeit with rather high variability in the weekend-to-weekday ratios. Only source regions, which are dominated by power plants or industry, do not show any reductions in NOx emissions on the weekend. The largest emission reductions on weekends are found for source regions with little industry, which are dominated by city emissions mostly from traffic. Our results indicate reductions of NOx emissions for Wuhan (China) during weekends, in contrast to earlier studies showing no such effect in China. This study can be extended to other Chinese cities as more high-resolution tropospheric NO2 columns become available from TROPOMI and GEMS in the next years. Separately calculated lifetimes for weekend and working days do not show longer lifetimes during weekends. This result also requires further investigation.

A short-term reduction in emissions attributable to the measures introduced to limit the spread of Covid-19 infections was found by comparing NOx emissions estimated with the EMG method. Since the method accounts for wind conditions and NOx lifetimes, it gives a better estimate of the impact of the Covid-19 measures on NOx emissions than comparisons of tropospheric NO2 column measurements as done in other studies. Strong NOx emission reductions during the first lockdown phase and a general tendency to lower emissions in 2020 than in 2019 are shown for Buenos Aires, New Delhi, and Madrid.

We conclude that the EMG method in combination with the high-spatial-resolution TROPOMI NO2 tropospheric column measurements allows us to investigate the high variability of NOx emissions and lifetimes on a global scale and on short time frames. Depending on the goal of the analysis, already a few days of measurements can be sufficient; for seasonal studies, depending on the local meteorological conditions, one to two seasons are sufficient, and it is no longer necessary to average over several years. The ability to estimate emissions over short periods will also allow policy makers to evaluate NOx emission regulations better and more quickly. The presented high variability should be further investigated using follow-up sensors on geostationary satellites, for example GEMS, Tempo, or Sentinel-4, which have the potential to also investigate the diurnal variability.

Appendix A

Table A1NOx source regions sorted with increasing latitude with mean NOx lifetime and emissions, wind speed, available days, season flag and weekly cycle flag (X: possible, –: not possible). Errors are the standard deviation derived by the EMG fitting procedure.

Download Print Version | Download XLSX

Figure A1Weekend-to-weekday ratio of NOx lifetimes (1 March 2018 to 29 February 2020) from left to right in increasing weekend-to-weekday NOx emission ratio (see Fig. 7). Data are separated into working and weekend days for 44 source regions. Weekend can be one or two days, and those days can also differ according to religious tradition, which is considered. The red line indicates the 1:1 line where weekend and weekday lifetimes are equal; below the line weekday lifetimes are longer, and over the line weekend lifetimes are longer. Error bars are the standard deviation derived by the EMG fitting procedure for lifetime estimates.


Data availability

TROPOMI data from July 2018 onwards are freely available via (S5P Data Hub2022). The wind, ozone, temperature, and boundary layer height data from the ERA5 reanalysis are freely available from the Copernicus Climate Change (C3S) climate data store (CDS) (, Copernicus2022). EDGAR v5.0 Global Air Pollutant Emissions are available at (EDGAR2022). CEMS data from power plants can be downloaded here: (EPA CEMS2022).

Author contributions

KL, AR, and JPB designed the study. KL performed the analysis and wrote the paper with contributions from JPB and AR.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Copernicus Sentinel-5P level-2 NO2 data are used in this study. Sentinel-5 Precursor is a European Space Agency (ESA) mission on behalf of the European Commission (EC). The TROPOMI payload is a joint development by ESA and the Netherlands Space Office (NSO). The Sentinel-5 Precursor ground-segment development has been funded by the ESA and with national contributions from the Netherlands, Germany, Belgium, and UK. Meteorological data are provided by ECMWF ERA5. EDGAR v5.0 Global Air Pollutant Emissions are provided by (last access: 12 February 2022), (last access: 12 February 2022).

Financial support

This research has been supported by the Deutsches Zentrum für Luft- und Raumfahrt (grant no. 50 EE 1709A).

The article processing charges for this open-access publication were covered by the University of Bremen.

Review statement

This paper was edited by Tanja Schuck and reviewed by three anonymous referees.


Atkinson, R., Baulch, D. L., Cox, R. A., Crowley, J. N., Hampson, R. F., Hynes, R. G., Jenkin, M. E., Rossi, M. J., and Troe, J.: Evaluated kinetic and photochemical data for atmospheric chemistry: Volume I – gas phase reactions of Ox, HOx, NOx and SOx species, Atmos. Chem. Phys., 4, 1461–1738,, 2004. a

Bauwens, M., Compernolle, S., Stavrakou, T., Müller, J.-F., Van Gent, J., Eskes, H., Levelt, P., van der A, R., Veefkind, J., Vlietinck, J., and Zehner, C.: Impact of coronavirus outbreak on NO2 pollution assessed using TROPOMI and OMI observations, Geophys. Res. Lett., 47, e2020GL087978,, 2020. a, b, c, d

Beirle, S., Platt, U., Wenig, M., and Wagner, T.: Weekly cycle of NO2 by GOME measurements: a signature of anthropogenic sources, Atmos. Chem. Phys., 3, 2225–2232,, 2003. a, b, c, d, e

Beirle, S., Boersma, K. F., Platt, U., Lawrence, M. G., and Wagner, T.: Megacity emissions and lifetimes of nitrogen oxides probed from space, Science, 333, 1737–1739, 2011. a, b, c, d, e, f, g, h, i, j, k, l, m, n

Beirle, S., Borger, C., Dörner, S., Li, A., Hu, Z., Liu, F., Wang, Y., and Wagner, T.: Pinpointing nitrogen oxide emissions from space, Sci. Adv., 5, eaax9800,, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n

Beirle, S., Borger, C., Dörner, S., Eskes, H., Kumar, V., de Laat, A., and Wagner, T.: Catalog of NOx emissions from point sources as derived from the divergence of the NO2 flux for TROPOMI, Earth Syst. Sci. Data, 13, 2995–3012,, 2021. a, b, c

Bovensmann, H., Burrows, J., Buchwitz, M., Frerick, J., Noël, S., Rozanov, V., Chance, K., and Goede, A.: SCIAMACHY: Mission objectives and measurement modes, J. Atmos. Sci., 56, 127–150, 1999. a

Burkholder, J., Sander, S., Abbatt, J., Barker, J., Cappa, C., Crounse, J., Dibble, T., Huie, R., Kolb, C., Kurylo, M. J., Orkin, V. L., Percival, C., Wilmouth, D. M., and Wine, P. H. : Chemical kinetics and photochemical data for use in atmospheric studies; evaluation number 19, Tech. Rep., Pasadena, CA, Jet Propulsion Laboratory, National Aeronautics and Space Administration, (last access: 24 February 2022), 2020. a

Burrows, J. P., Weber, M., Buchwitz, M., Rozanov, V., Ladstätter-Weißenmayer, A., Richter, A., DeBeek, R., Hoogen, R., Bramstedt, K., Eichmann, K.-U., Eisinger, M., and Perner, D.: The global ozone monitoring experiment (GOME): Mission concept and first scientific results, Atmos. Sci., 56, 151–175, 1999. a

Copernicus: Copernicus Climate Data Store,, 2022. a

Crippa, M., Guizzardi, D., Muntean, M., Schaaf, E., Dentener, F., van Aardenne, J. A., Monni, S., Doering, U., Olivier, J. G. J., Pagliari, V., and Janssens-Maenhout, G.: Gridded emissions of air pollutants for the period 1970–2012 within EDGAR v4.3.2, Earth Syst. Sci. Data, 10, 1987–2013,, 2018. a, b

Crippa, M., Oreggioni, G., Guizzardi, D., Muntean, M., Schaaf, E., Lo Vullo, E., Solazzo, E., Monforti-Ferrario, F., Olivier, J. G., and Vignati, E.: Fossil CO2 and GHG emissions of all world countries, Luxemburg, Publication Office of the European Union,, 2019. a

Crippa, M., Solazzo, E., Huang, G., Guizzardi, D., Koffi, E., Muntean, M., Schieberle, C., Friedrich, R., and Janssens-Maenhout, G.: High resolution temporal profiles in the Emissions Database for Global Atmospheric Research, Sci. Data, 7, 1–17,, 2020. a, b

de Foy, B., Wilkins, J. L., Lu, Z., Streets, D. G., and Duncan, B. N.: Model evaluation of methods for estimating surface emissions and chemical lifetimes from satellite data, Atmos. Environ., 98, 66–77, 2014. a

de Foy, B., Lu, Z., Streets, D. G., Lamsal, L. N., and Duncan, B. N.: Estimates of power plant NOx emissions and lifetimes from OMI NO2 satellite retrievals, Atmos. Environ., 116, 1–11, 2015. a

Dickerson, R. R., Stedman, D. H., and Delany, A. C.: Direct measurements of ozone and nitrogen dioxide photolysis rates in the troposphere, J. Geophys. Res.-Ocean., 87, 4933–4946, 1982. a

Dimitropoulou, E., Hendrick, F., Pinardi, G., Friedrich, M. M., Merlaud, A., Tack, F., De Longueville, H., Fayt, C., Hermans, C., Laffineur, Q., Fierens, F., and Van Roozendael, M.: Validation of TROPOMI tropospheric NO2 columns using dual-scan multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements in Uccle, Brussels, Atmos. Meas. Tech., 13, 5165–5191,, 2020. a

Dittman, M. G., Ramberg, E., Chrisp, M., Rodriguez, J. V., Sparks, A. L., Zaun, N. H., Hendershot, P., Dixon, T., Philbrick, R. H., and Wasinger, D.: Nadir ultraviolet imaging spectrometer for the NPOESS Ozone Mapping and Profiler Suite (OMPS), in: Earth Observing Systems VII, Vol. 4814, 111–119, International Society for Optics and Photonics,, 2002. a

EDGAR: EDGAR v5.0 GlobalAir Pollutant Emissions,, last access: 21 February 2022. a

EPA CEMS: Environmental Protection Agency’s Continuous Emission Monitoring System,, last access: 21 February 2022. a

Eskes, H. and Eichmann, K.: S5P MPC Product Readme Nitrogen Dioxide, Tech. rep., Report S5P-MPC-KNMI-PRF-NO2, version 2.1, 17 November 2021, ESA, (last access: 21 February 2022), 2021. a, b

Georgoulias, A. K., van der A, R. J., Stammes, P., Boersma, K. F., and Eskes, H. J.: Trends and trend reversal detection in 2 decades of tropospheric NO2 satellite observations, Atmos. Chem. Phys., 19, 6269–6294,, 2019. a, b

Goldberg, D. L., Lu, Z., Streets, D. G., de Foy, B., Griffin, D., McLinden, C. A., Lamsal, L. N., Krotkov, N. A., and Eskes, H.: Enhanced Capabilities of TROPOMI NO2: Estimating NOX from North American Cities and Power Plants, Environ. Sci. Technol., 53, 12594–12601, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m

Goldberg, D. L., Anenberg, S. C., Griffin, D., McLinden, C. A., Lu, Z., and Streets, D. G.: Disentangling the impact of the COVID-19 lockdowns on urban NO2 from natural variability, Geophys. Res. Lett., 47, e2020GL089269, 2020. a, b, c

Goldberg, D. L., Anenberg, S. C., Kerr, G. H., Mohegh, A., Lu, Z., and Streets, D. G.: TROPOMI NO2 in the United States: A detailed look at the annual averages, weekly cycles, effects of temperature, and correlation with surface NO2 concentrations, Earth's Future, 9, e2020EF001665,, 2021. a, b, c

Griffin, D., Zhao, X., McLinden, C. A., Boersma, F., Bourassa, A., Dammers, E., Degenstein, D., Eskes, H., Fehr, L., Fioletov, V., Hayden, K., Kharol, S. K., Li, S.-M., Makar, P., Martin, R. V., Mihele, C., Mittermeier, R. L., Krotkov, N., Sneep, M., Lamsal, L. N., ter Linden, M., van Geffen, Jos., Veefkind, P., and Wolde, M.: High-resolution mapping of nitrogen dioxide with TROPOMI: First results and validation over the Canadian oil sands, Geophys. Res. Lett., 46, 1049–1060, 2019. a, b

Holland, F., Hofzumahaus, A., Schäfer, J., Kraus, A., and Pätz, H.-W.: Measurements of OH and HO2 radical concentrations and photolysis frequencies during BERLIOZ, J. Geophys. Res.- Atmos., 108, PHO 2-1–PHO 2-23,, 2003. a

Ialongo, I., Hakkarainen, J., Hyttinen, N., Jalkanen, J.-P., Johansson, L., Boersma, K. F., Krotkov, N., and Tamminen, J.: Characterization of OMI tropospheric NO2 over the Baltic Sea region, Atmos. Chem. Phys., 14, 7795–7805,, 2014. a, b

Ialongo, I., Virta, H., Eskes, H., Hovila, J., and Douros, J.: Comparison of TROPOMI/Sentinel-5 Precursor NO2 observations with ground-based measurements in Helsinki, Atmos. Meas. Tech., 13, 205–218,, 2020. a

Jacob, D. J.: Introduction to atmospheric chemistry, Princeton University Press, ISBN 9781400841547, 1999. a, b

Judd, L. M., Al-Saadi, J. A., Szykman, J. J., Valin, L. C., Janz, S. J., Kowalewski, M. G., Eskes, H. J., Veefkind, J. P., Cede, A., Mueller, M., Gebetsberger, M., Swap, R., Pierce, R. B., Nowlan, C. R., Abad, G. G., Nehrir, A., and Williams, D.: Evaluating Sentinel-5P TROPOMI tropospheric NO2 column densities with airborne and Pandora spectrometers near New York City and Long Island Sound, Atmos. Meas. Tech., 13, 6113–6140,, 2020. a

Kunhikrishnan, T., Lawrence, M. G., von Kuhlmann, R., Richter, A., Ladstätter-Weißenmayer, A., and Burrows, J. P.: Analysis of tropospheric NOx over Asia using the model of atmospheric transport and chemistry (MATCH-MPIC) and GOME-satellite observations, Atmos. Environ., 38, 581–596, 2004. a

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

Leue, C., Wenig, M., Wagner, T., Klimm, O., Platt, U., and Jähne, B.: Quantitative analysis of NOx emissions from Global Ozone Monitoring Experiment satellite image sequences, J. Geophys. Res.-Atmos., 106, 5493–5505, 2001. a

Levelt, P. F., van den Oord, G. H., Dobber, M. R., Malkki, A., Visser, H., de Vries, J., Stammes, P., Lundell, J. O., and Saari, H.: The ozone monitoring instrument, IEEE Trans. Geosci. Remote Sens., 44, 1093–1101, 2006. a

Levelt, P. F., Stein Zweers, D. C., Aben, I., Bauwens, M., Borsdorff, T., De Smedt, I., Eskes, H. J., Lerot, C., Loyola, D. G., Romahn, F., Stavrakou, T., Theys, N., Van Roozendael, M., Veefkind, J. P., and Verhoelst, T.: Air quality impacts of COVID-19 lockdown measures detected from space using high spatial resolution observations of multiple trace gases from Sentinel-5P/TROPOMI, Atmos. Chem. Phys. Discuss. [preprint],, in review, 2021. a

Liu, F., Beirle, S., Zhang, Q., Dörner, S., He, K., and Wagner, T.: NOx lifetimes and emissions of cities and power plants in polluted background estimated by satellite observations, Atmos. Chem. Phys., 16, 5283–5298,, 2016. a, b

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, eabc2992,, 2020. a, b

Lorente, A., Boersma, K., Eskes, H., Veefkind, J., Van Geffen, J., de Zeeuw, M., van der Gon, H. D., Beirle, S., and Krol, M.: Quantification of nitrogen oxides emissions from build-up of pollution over Paris with TROPOMI, Sci. Rep., 9, 1–10,, 2019. a, b, c, d, e, f, g, h, i, j, k, l

Lu, K. D., Hofzumahaus, A., Holland, F., Bohn, B., Brauers, T., Fuchs, H., Hu, M., Häseler, R., Kita, K., Kondo, Y., Li, X., Lou, S. R., Oebel, A., Shao, M., Zeng, L. M., Wahner, A., Zhu, T., Zhang, Y. H., and Rohrer, F.: Missing OH source in a suburban environment near Beijing: observed and modelled OH and HO2 concentrations in summer 2006, Atmos. Chem. Phys., 13, 1057–1080,, 2013. a

Lu, Z., Streets, D. G., de Foy, B., Lamsal, L. N., Duncan, B. N., and Xing, J.: Emissions of nitrogen oxides from US urban areas: estimation from Ozone Monitoring Instrument retrievals for 2005–2014, Atmos. Chem. Phys., 15, 10367–10383,, 2015. a

Martin, R. V., Jacob, D. J., Chance, K., Kurosu, T. P., Palmer, P. I., and Evans, M. J.: Global inventory of nitrogen oxide emissions constrained by space-based observations of NO2 columns, J. Geophys. Res.-Atmos., 108, D17,, 2003. a, b, c

Molina, M. J. and Molina, L. T.: Megacities and atmospheric pollution, J. Air Waste Manage., 54, 644–680, 2004. a

Munro, R., Eisinger, M., Anderson, C., Callies, J., Corpaccioli, E., Lang, R., Lefebvre, A., Livschitz, Y., and Albinana, A. P.: GOME-2 on MetOp, in: Proc. of The 2006 EUMETSAT Meteorological Satellite Conference, Helsinki, Finland, Vol. 1216, 48–53, 2006. a

Pommier, M., McLinden, C. A., and Deeter, M.: Relative changes in CO emissions over megacities based on observations from space, Geophys. Res. Lett., 40, 3766–3771, 2013. a, b, c

Seinfeld, J. H. and Pandis, S. N.: Atmospheric Chemistry and Physics, A Wiley-Inter Science Publication, John Wiley & Sons Inc., Hoboken, New Jersey, ISBN 978-0-471-72018-8, 2006. a, b, c, d

S5P Data Hub: Sentinel-5P Pre-Operations Data Hub,, last access: 21 February 2022. a

Shah, V., Jacob, D. J., Li, K., Silvern, R. F., Zhai, S., Liu, M., Lin, J., and Zhang, Q.: Effect of changing NOx lifetime on the seasonality and long-term trends of satellite-observed tropospheric NO2 columns over China, Atmos. Chem. Phys., 20, 1483–1495,, 2020. a, b

Smith, S. C., Lee, J. D., Bloss, W. J., Johnson, G. P., Ingham, T., and Heard, D. E.: Concentrations of OH and HO2 radicals during NAMBLEX: measurements and steady state analysis, Atmos. Chem. Phys., 6, 1435–1453,, 2006. a

Stavrakou, T., Müller, J.-F., Boersma, K. F., De Smedt, I., and Van Der A, R.: Assessing the distribution and growth rates of NOx emission sources by inverting a 10-year record of NO2 satellite columns, Geophys. Res. Lett., 35, 2–6,, 2008. a, b, c

Stavrakou, T., Müller, J.-F., Boersma, K. F., van der A, R. J., Kurokawa, J., Ohara, T., and Zhang, Q.: Key chemical NOx sink uncertainties and how they influence top-down emissions of nitrogen oxides, Atmos. Chem. Phys., 13, 9057–9082,, 2013. a

Stavrakou, T., Müller, J.-F., Bauwens, M., Boersma, K., and van Geffen, J.: Satellite evidence for changes in the NO2 weekly cycle over large cities, Sci. Rep., 10, 1–9, 2020. a, b, c, d

Stocker, T.: Climate change 2013: the physical science basis: Working Group I contribution to the Fifth assessment report of the Intergovernmental Panel on Climate Change, Cambridge University Press, (last access: 24 February 2022), 2014. a

Valin, L., Russell, A., and Cohen, R.: Variations of OH radical in an urban plume inferred from NO2 column measurements, Geophys. Res. Lett., 40, 1856–1860, 2013. a, b, c, d

Valin, L. C., Russell, A. R., and Cohen, R. C.: Chemical feedback effects on the spatial patterns of the NOx weekend effect: a sensitivity analysis, Atmos. Chem. Phys., 14, 1–9,, 2014. a

van der A, R. J., Eskes, H. J., Boersma, K. F., van Noije, T. P. C., Van Roozendael, M., De Smedt, I., Peters, D. H. M. U., and Meijer, E. W.: Trends, seasonal variability and dominant NOx source derived from a ten year record of NO2 measured from space, J. Geophys. Res.-Atmos., 113, D04302,, 2008. a

van Geffen, J., Eskes, H., Boersma, K., Maasakkers, J., and Veefkind, J.: TROPOMI ATBD of the total and tropospheric NO2 data products (issue 1.2. 0), Royal Netherlands Meteorological Institute (KNMI), De Bilt, the Netherlands, s5P-KNMI-L2-0005-RP, (last access: 24 February 2022), 2018.  a, b, c

van Geffen, J., Eskes, H., Compernolle, S., Pinardi, G., Verhoelst, T., Lambert, J.-C., Sneep, M., ter Linden, M., Ludewig, A., Boersma, K. F., and Veefkind, J. P.: Sentinel-5P TROPOMI NO2 retrieval: impact of version v2.2 improvements and comparisons with OMI and ground-based data, Atmos. Meas. Tech. Discuss. [preprint],, in review, 2021. a, b

Veefkind, J., Aben, I., McMullan, K., Förster, H., De Vries, J., Otter, G., Claas, J., Eskes, H., De Haan, J., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P. F.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83, 2012. a, b

Verhoelst, T., Compernolle, S., Pinardi, G., Lambert, J.-C., Eskes, H. J., Eichmann, K.-U., Fjæraa, A. M., Granville, J., Niemeijer, S., Cede, A., Tiefengraber, M., Hendrick, F., Pazmiño, A., Bais, A., Bazureau, A., Boersma, K. F., Bognar, K., Dehn, A., Donner, S., Elokhov, A., Gebetsberger, M., Goutail, F., Grutter de la Mora, M., Gruzdev, A., Gratsea, M., Hansen, G. H., Irie, H., Jepsen, N., Kanaya, Y., Karagkiozidis, D., Kivi, R., Kreher, K., Levelt, P. F., Liu, C., Müller, M., Navarro Comas, M., Piters, A. J. M., Pommereau, J.-P., Portafaix, T., Prados-Roman, C., Puentedura, O., Querel, R., Remmers, J., Richter, A., Rimmer, J., Rivera Cárdenas, C., Saavedra de Miguel, L., Sinyakov, V. P., Stremme, W., Strong, K., Van Roozendael, M., Veefkind, J. P., Wagner, T., Wittrock, F., Yela González, M., and Zehner, C.: Ground-based validation of the Copernicus Sentinel-5p TROPOMI NO 2 measurements with the NDACC ZSL-DOAS, MAX-DOAS and Pandonia global networks, Atmos. Meas. Tech., 14, 481–510,, 2021. a, b

Zheng, B., Tong, D., Li, M., Liu, F., Hong, C., Geng, G., Li, H., Li, X., Peng, L., Qi, J., Yan, L., Zhang, Y., Zhao, H., Zheng, Y., He, K., and Zhang, Q.: Trends in China's anthropogenic emissions since 2010 as the consequence of clean air actions, Atmos. Chem. Phys., 18, 14095–14111,, 2018. a

Short summary
In this study, we investigated short time variability of NOx emissions and lifetimes on a global scale. We combined 2 years of satellite Sentinel-5P TROPOMI tropospheric NO2 column data with wind data. Fifty NOx sources distributed around the world are analyzed. The retrieved emissions show a clear seasonal dependence. NOx lifetime shows a latitudinal dependence but only a week seasonal dependence. NOx emissions show a clear weekly pattern which in contrast is not visible for NOx lifetimes.
Final-revised paper