Quantifying NOx emissions in Egypt using TROPOMI observations

Abstract

Urban areas and industrial facilities, which concentrate the majority of human activity and industrial production, are major sources of air pollutants, with serious implications for human health and global climate. For most of these pollutants, emission inventories are often highly uncertain, especially in developing countries. Spaceborne measurements from the TROPOMI instrument, on board the Sentinel-5 Precursor satellite, are used to retrieve nitrogen dioxide (${\mathrm{NO}}_{\mathrm{2}}$) column densities at high spatial resolution. Here, we use 2 years of TROPOMI retrievals to map nitrogen oxide (${\mathrm{NO}}_x$ $=$ NO $+$ ${\mathrm{NO}}_{\mathrm{2}}$) emissions in Egypt with a top-down approach using the continuity equation in steady state. Emissions are expressed as the sum of a transport term and a sink term representing the three-body reaction comprising ${\mathrm{NO}}_{\mathrm{2}}$ and hydroxyl radical (OH). This sink term requires information on the lifetime of ${\mathrm{NO}}_{\mathrm{2}}$, which is calculated with the use of the CAMS near-real-time temperature and OH concentration fields. We compare this derived lifetime with the lifetime inferred from the fitting of ${\mathrm{NO}}_{\mathrm{2}}$ line density profiles in large plumes with an exponentially modified Gaussian function. This comparison, which is conducted for different samples of ${\mathrm{NO}}_{\mathrm{2}}$ patterns above the city of Riyadh, provides information on the reliability of the CAMS near-real-time OH concentration fields; it also provides some hint on the vertical levels that best represent typical pollution sources in industrial areas and megacities in the Middle East region. In Egypt, total emissions of ${\mathrm{NO}}_x$ are dominated by the sink term, but they can be locally dominated by wind transport, especially along the Nile where human activities are concentrated. Megacities and industrial regions clearly appear as the largest sources of ${\mathrm{NO}}_x$ emissions in the country. Our top-down model infers emissions with a marked annual variability. By looking at the spatial distribution of emissions at the scale of different cities with different industrial characteristics, it appears that this variability is consistent with national electricity consumption. We detect lower emissions on Fridays, which are inherent to the social norm of the country, and quantify the drop in emissions in 2020 due to the COVID-19 pandemic. Overall, our estimations of ${\mathrm{NO}}_x$ emissions for Egypt are 7.0 % higher than the CAMS-GLOB-ANT_v4.2 inventory and significantly differ in terms of seasonality.

1 Introduction

Economic growth in developing countries has led to a strong increase in urban air pollution . Among the different pollutants, nitrogen oxides are key species. They are generally the products of fuel combustion, such as the burning of hydrocarbons in the air at high temperature. The main sources of these compounds are vehicle engines, but also heavy industrial facilities such as power plants, iron and steel mills , and cement kilns . Their accumulation in the lowest layers of the troposphere contributes to the formation of smog and acid rain . They also have a significant effect on human health, as they can cause various respiratory diseases . To deal with these phenomena, national and regional governments generally enact a series of air pollution control strategies, which typically take the form of bans on certain polluting technologies, with the aim of reducing the concentration of pollutants at the local level to targets that must be achieved within a given timeframe. These strategies, which also help drive technological innovation, have had a significant effect in Europe .

In Egypt, population growth, urbanisation, socio-economic development and the associated increase in the vehicle fleet led to a major degradation of air quality in the last decades, especially in highly populated areas such as Greater Cairo and the Nile Delta , which include the majority of the population. The Ministry of State for the Environment has thus initiated new policies that aim to reduce pollution levels throughout the country, through technical mitigation of emissions, emission standards for vehicles and intersectoral collaboration . However, Egypt, like most developing countries, lacks the local infrastructure to access detailed information on technical factors such as energy consumption or fuel type and technology, leading to discrepancies in inventories . As a consequence, the monitoring of emissions, which is important to evaluate the effects of air pollution control policies, is of limited reliability.

To overcome the uncertainties in the emission inventories, the use of independent observation systems is becoming increasingly prevalent. In this study, we investigate the use of satellite remote sensing of atmospheric concentrations to improve the quantification of ${\mathrm{NO}}_x$ emissions in Egypt. Spectrally resolved satellite measurements of solar backscattered radiation enable the quantification of ${\mathrm{NO}}_{\mathrm{2}}$ and other trace gases absorbing in the UV–visible spectral range based on their characteristic spectral absorption patterns. Tropospheric vertical column densities, i.e. vertically integrated ${\mathrm{NO}}_{\mathrm{2}}$ concentrations in the troposphere, have been providing information on the spatial distribution of tropospheric ${\mathrm{NO}}_{\mathrm{2}}$ at a global scale for nearly 30 years, allowing the identification of different sources of ${\mathrm{NO}}_x$ and the quantification of the associated emissions . In October 2017, the Sentinel-5 Precursor satellite was launched. Its main instrument is the TROPOspheric Monitoring Instrument (TROPOMI), which provides tropospheric ${\mathrm{NO}}_{\mathrm{2}}$ column densities at high spatial resolution with a large swath width and with a daily frequency . By applying the steady-state continuity equation , it is possible to build a top-down model that directly quantifies ${\mathrm{NO}}_x$ emissions from these ${\mathrm{NO}}_{\mathrm{2}}$ column densities, provided that some key parameters (wind, temperature, hydroxyl radical concentration and concentration ratio between ${\mathrm{NO}}_x$ and ${\mathrm{NO}}_{\mathrm{2}}$) are correctly estimated. This model is used to quantify the anthropogenic ${\mathrm{NO}}_x$ emissions in Egypt for a 2-year period, from November 2018 to November 2020.

This paper is organised as follows: Sect. 2 provides a description of the datasets used in this study. Section 3 explains the build-up and the limits of the top-down approach used to quantify ${\mathrm{NO}}_x$ emissions in Egypt. It also presents a method for validating the model parameters by using ${\mathrm{NO}}_{\mathrm{2}}$ line density profiles over Riyadh, Saudi Arabia. Section 4 presents the analysis of this validation method. It presents the location of the main ${\mathrm{NO}}_x$ sources in Egypt and evaluates the vertical sensitivity of the model. It also assesses the ability of the model to show less human activity on Fridays and during the lockdown that took place during the COVID-19 pandemic. It finally confronts the inferred emissions with different inventories in terms of amplitude and seasonality. Section 5 presents our conclusion and general remarks.

2 Instrumentation and data

2.1

TROPOMI ${\mathrm{NO}}_{\mathrm{2}}$ retrievals

The TROPOspheric Atmosphere Monitoring Instrument (TROPOMI), on board the European Space Agency's (ESA) Sentinel-5 Precursor (S-5P) satellite, provides measurements for atmospheric composition. TROPOMI is a spectrometer observing wavelengths in infrared, visible and ultraviolet light at around 13:30 LT (local time). The UV–visible spectral band at 405–465 nm is used to retrieve ${\mathrm{NO}}_{\mathrm{2}}$. Other spectral bands are used to observe methane, formaldehyde, sulfur dioxide, carbon monoxide and ozone, as well as aerosols and cloud physical properties. The very high spatial resolution offered by TROPOMI (originally $\mathrm{3.5}\times \mathrm{7}$ km${}^{\mathrm{2}}$ at nadir, improved to $\mathrm{3.5}\times \mathrm{5.5}$ km${}^{\mathrm{2}}$ since 6 August 2019) provides unprecedented information on tropospheric ${\mathrm{NO}}_{\mathrm{2}}$. Its large swath width ($\sim \mathrm{2600}$ km) makes it possible to construct ${\mathrm{NO}}_{\mathrm{2}}$ images on large spatial scales. Those images greatly improve the potential for detecting highly localised pollution plumes above the ground, identifying small-scale emission sources but also estimating emissions from megacities, industrial facilities and biomass burning. We use TROPOMI ${\mathrm{NO}}_{\mathrm{2}}$ retrievals (L2 data, OFFL stream, product version 1.0.0 and 1.1.0 successively) from November 2018 to November 2020 over Egypt. We also use them over Saudi Arabia, and more specifically over the city of Riyadh, to evaluate the reliability of other parameters. This will be explained in Sect. 3.3. Both countries have an arid climate, which offers a large number of clear-sky days throughout the year, enabling the calculation of monthly averages based on more than 20 observations. They are also the host to many large plumes of pollutants due to high human concentrations along rivers and around megacities, which allows us to observe high ${\mathrm{NO}}_{\mathrm{2}}$ concentration patterns with a high signal-to-noise ratio. TROPOMI products provide a quality assurance value $q_{\mathrm{a}}$, which ranges from 0 (no data) to 1 (high-quality data). For our analysis of concentrations, we selected ${\mathrm{NO}}_{\mathrm{2}}$ retrievals with $q_{\mathrm{a}}$ values greater than 0.75, which systematically correspond to clear-sky conditions . TROPOMI soundings are gridded at a spatial resolution of $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ with daily coverage. This resolution is lower than that of the instrument; the gridding thus provides a grid for which most ${\mathrm{NO}}_{\mathrm{2}}$ columns correspond to one or more measurements. The observed plumes remain correctly resolved. Cells without measurements are infrequent, which facilitates the calculation of derivatives.

2.2 Wind data

The horizontal wind $\mathit{w}=(u,v)$ is taken from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 data archive (fifth generation of atmospheric reanalyses) at a horizontal resolution of $\mathrm{0.25}{}^{\circ }\times \mathrm{0.25}{}^{\circ }$ on 37 pressure levels . The hourly values have been linearly interpolated to the TROPOMI orbit timestamp and re-gridded to a $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ resolution.

2.3 CAMS real-time fields

The Copernicus Atmospheric Monitoring Service (CAMS) global near-real-time service provides analyses and forecasts for reactive gases, greenhouse gases and aerosols on 25 pressure levels with a horizontal resolution of $\mathrm{0.4}{}^{\circ }\times \mathrm{0.4}{}^{\circ }$ and a temporal resolution of 3 h . For the calculation of ${\mathrm{NO}}_x$ emissions from TROPOMI observations, we use CAMS concentration fields of nitrogen oxides (NO and ${\mathrm{NO}}_{\mathrm{2}}$) and the hydroxyl radical (OH). We also use the CAMS temperature field $T$. NO and ${\mathrm{NO}}_{\mathrm{2}}$ concentrations are used to account for chemical processes that take place in polluted air. Anthropogenic activities produce mainly NO, which is transformed into ${\mathrm{NO}}_{\mathrm{2}}$ by reaction with ozone ${\mathrm{O}}_{\mathrm{3}}$. ${\mathrm{NO}}_{\mathrm{2}}$ is then photolysed during the day, reforming NO . This photochemical equilibrium between NO and ${\mathrm{NO}}_{\mathrm{2}}$ can be highlighted with the ${\mathrm{NO}}_x:{\mathrm{NO}}_{\mathrm{2}}$ concentration ratio, whose value depends on local conditions, allowing us to perform a conversion from ${\mathrm{NO}}_{\mathrm{2}}$ production to ${\mathrm{NO}}_x$ emissions. The reason for the use of OH is different. OH is the main oxidant that controls the ability of the atmosphere to remove pollutants such as ${\mathrm{NO}}_{\mathrm{2}}$ . It is mainly produced during daylight hours by interaction between water and atomic oxygen produced by ozone dissociation . In air that is directly influenced by pollution, another source of OH is due to a reaction between NO and ${\mathrm{HO}}_{\mathrm{2}}$. This reaction, referred to as the ${\mathrm{NO}}_x$ recycling mechanism, illustrates the nonlinear dependence of the OH concentration on ${\mathrm{NO}}_{\mathrm{2}}$ . Since the OH lifetime is typically less than a second, its concentration in the troposphere is very low and difficult to measure. As a consequence, global analyses, which estimate OH concentrations from other variable species , provide a representation for OH concentrations with high associated uncertainties. Therefore, the CAMS OH concentrations are used here to account for the ${\mathrm{NO}}_{\mathrm{2}}$ oxidation to form nitric acid (${\mathrm{HNO}}_{\mathrm{3}}$), whose representation is explained in Sect. 3.1. Finally, the temperature field is used to control variations in the kinetic parameters . The hourly values are also linearly interpolated to the TROPOMI orbit timestamp and re-gridded to a $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ resolution.

2.4 Background removal

Detecting traces of anthropogenic emissions in TROPOMI ${\mathrm{NO}}_{\mathrm{2}}$ images is not a straightforward process. The ${\mathrm{NO}}_{\mathrm{2}}$ signal from a sparsely populated area or a small industrial facility may be covered by numerical noise or by the signal generated by natural ${\mathrm{NO}}_x$ emissions. In the absence of anthropogenic sources, TROPOMI observes ${\mathrm{NO}}_{\mathrm{2}}$ concentrations, which constitute a tropospheric background of $\sim \mathrm{0.5}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$. At the global scale, this background is the result of different sources. In the lower troposphere, natural ${\mathrm{NO}}_x$ emissions are mostly due to fires and soil emissions . In the upper troposphere however, sources include lightning, convective injection and downwelling from the stratosphere , but the factors controlling the resulting concentrations are poorly understood. According to state-of-the-art estimates, anthropogenic ${\mathrm{NO}}_x$ accounts for most of the emissions at the global scale, whereas natural emissions from fires, soils and lightning are less significant at the global scale and do not exceed a share of 35 % combined , although associated errors can be very high. In eastern China, the non-anthropogenic share of total ${\mathrm{NO}}_x$ emissions is variable but does not exceed 20 % . With Egypt being a desertic region and not very conductive to lightning, we expect the share of those non-anthropogenic emissions to be smaller. To estimate anthropogenic ${\mathrm{NO}}_x$ emissions, it is therefore necessary to remove this share.

Figure 1

(a) Part of Egypt centred on Nile River. Within this domain, pink cells represent locations with an average human density above 100 inhabitants per square kilometre, brown cells represent locations with industrial facilities outside cities and yellow cells represent locations in their vicinity. These cells constitute the mask used to calculate anthropogenic emissions. Outside this mask, green cells represent areas which do not correspond to any of the three criteria, considered to be void of human activity. Five large cities in the country and the industrial area of Ain Sokhna are denoted with stars. Two smaller domains centred around the cities of Asyut and Aswan are represented with dotted lines; their use is presented in Sect. 4.6. (b) TROPOMI observation of ${\mathrm{NO}}_{\mathrm{2}}$ slant column densities above Nile valley on 3 January 2019. White pixels correspond to areas with low-quality data ($q_{\mathrm{a}}<\mathrm{0.75}$) or no data.

[Figure omitted. See PDF]

With an atmospheric lifetime of about a few hours, the presence of ${\mathrm{NO}}_{\mathrm{2}}$ is relatively short. Consequently, the majority of ${\mathrm{NO}}_{\mathrm{2}}$ is not transported far downwind from its sources. Thus, near-surface ${\mathrm{NO}}_{\mathrm{2}}$ concentrations are generally high over industrial facilities and densely populated areas that need to be identified. Because Egypt's population is almost entirely located along the Nile River and its delta, the study of ${\mathrm{NO}}_x$ emissions in this country cannot therefore be reduced to the study of a small number of point sources, as would be the case for several other parts of the Middle East region, and must be carried out in the form of a mapping of the country. Further discussion is provided in Sect. 3.3. To identify urban areas in Egypt, we use the Socioeconomic Data and Applications Center (SEDAC) GRUMP (Global Rural-Urban Mapping Project) data archive, which comprises eight global datasets, including a population density grid provided at a resolution of 30 arcsec, with population estimates normalised for the year 2000 . We combine this database with field data giving the location of industrial facilities from energy-intensive industries in the region. Data have been retrieved from the Global Energy Observatory for five oil- and gas-fired power plants, from the work of for flaring sites, and from the work of Steven J. Davis and Dan Tong (unpublished data) for cement plants; links are at the end of this article. The locations of other industrial facilities, such as aluminium and iron smelters, were obtained from various sources.

These datasets are used to remove the non-anthropogenic part of the ${\mathrm{NO}}_x$ emissions signal. We conduct this removal by subtracting the mean emissions over areas without human activity from the mean emissions over industrial and densely populated areas. In order to perform this distinction between these two types of areas, our study is carried out using a mask within a $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ grid. A grid cell is considered to be part of the mask if it has a population density higher than a threshold of 100 inhabitants per square kilometre, or if its centre is close to an industrial facility. Otherwise, the cell is considered to be part of the “background”, i.e. outside the mask. In order to avoid any smearing that would correspond to abnormally high emissions outside urban and industrial centres (which can happen if the wind is poorly estimated), transition cells (in the immediate vicinity of the mentioned mask cells) are also considered to be mask cells. Figure  shows the distinction between mask cells and background cells on our domain in Egypt that lies between parallels 23 and 32${}^{\circ }$ N and meridians 29 and 34${}^{\circ }$ E. Most of the mask cells are located in the Nile area. Some mask cells are also found on the coast or in isolated parts in the desert. They correspond to remote industrial facilities, including major flaring sites, or sparsely populated industrial centres such as Ain Sokhna's industrial area. The domain comprises $n_{\mathrm{m}}=\mathrm{949}$ mask cells and $n_{\mathrm{b}}=\mathrm{3692}$ background cells. The mathematical description of the background removal is outlined in Sect. 3.4.

2.5 Emission inventories

We compare TROPOMI-derived ${\mathrm{NO}}_x$ emissions to the Emissions Database for Global Atmospheric Research (EDGARv5.0) for 2020 and the CAMS global anthropogenic emissions (CAMS-GLOB-ANT_v4.2) inventory released in 2020. Both datasets provide $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ gridded emissions for different sectors on a monthly basis. EDGARv5.0 emissions are based on activity data (population, energy production, fossil fuel extraction, industrial processes, agricultural statistics, etc.) derived from the International Energy Agency (IEA) and the Food and Agriculture Organization (FAO), corresponding emission factors, national and regional information on technology mix data, and end-of-pipe measurements. The inventory covers the years 1970–2015 and differs from the previous version EDGARv4.3.2, which does not use splitting factors derived from the Energy Information Administration (EIA) data on fuel consumption of coal, oil and natural gas for specific countries . CAMS-GLOB-ANT_v4.2 is developed within the framework of the Copernicus Atmospheric Monitoring Service . For this inventory, ${\mathrm{NO}}_x$ emissions are based on various existing sectors in the EDGARv4.3.2 emissions from 2000–2012, which are used as a basis for 2010 emissions and are extrapolated to the current year using 2011–2014 sector-based trends from the Community Emissions Data System (CEDS) inventory . From one inventory to another, the names and definitions of the sectors may vary. In EDGARv5.0 and CAMS-GLOB-ANT_v4.2, the emissions for a given country are derived from the type of technologies used, the dependence of emission factors on fuel type, combustion conditions, and activity data and low-resolution emission factors .

3 Method

3.1

Calculation of ${\mathrm{NO}}_{\mathrm{2}}$ production from TROPOMI observations

As a first step, we use tropospheric ${\mathrm{NO}}_{\mathrm{2}}$ vertical column densities ${\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}$ to derive top-down ${\mathrm{NO}}_{\mathrm{2}}$ production maps. Vertical column densities ${\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}$ are obtained from TROPOMI slant column densities using an air mass factor (AMF) which is also provided by TROPOMI. Previous studies have shown that the use of the AMF is a source of structural uncertainty in ${\mathrm{NO}}_{\mathrm{2}}$ retrievals . In polluted environments, this source of uncertainty becomes non-negligible. Here, the AMF does not vary much temporally throughout the studied period and is around 1.6 for mask cells and around 1.9 for background cells. The difference between the two types of cells is probably due to a different albedo between the urban environment and desert areas. Using the horizontal wind $\mathit{w}$, the ${\mathrm{NO}}_{\mathrm{2}}$ flux is given as ${\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}\mathit{w}$. The divergence of this flux can be added to the local time derivative $\frac{\partial {\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}}{\partial t}$ to balance ${\mathrm{NO}}_{\mathrm{2}}$ sources $e_{{\mathrm{NO}}_{\mathrm{2}}}$ and sinks $s_{{\mathrm{NO}}_{\mathrm{2}}}$ according to the continuity equation 1$${\displaystyle \frac{\partial {\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}}{\partial t}}+\mathrm{div}\left({\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}\mathit{w}\right)=e_{{\mathrm{NO}}_{\mathrm{2}}}-s_{{\mathrm{NO}}_{\mathrm{2}}}.$$ In steady state, the time derivative disappears and the mass balance is reduced to three terms. The ${\mathrm{NO}}_{\mathrm{2}}$ production can thus be estimated by taking the combined effect of atmospheric transport losses and the different sinks. For the transport term, we calculate numerical derivatives with a fourth-order central-finite difference scheme for each cell of the domain. Moreover, since the local overpass time of TROPOMI occurs in the middle of the day, ${\mathrm{NO}}_{\mathrm{2}}$ losses are largely dominated by the three-body reaction involving ${\mathrm{NO}}_{\mathrm{2}}$ and OH . Two channels have been identified for this reaction , leading to the production of nitric acid ${\mathrm{HNO}}_{\mathrm{3}}$ and pernitrous acid HOONO: $$\begin{array}{cc}{\displaystyle }& {\displaystyle }{\mathrm{NO}}_{\mathrm{2}}+\mathrm{OH}+M\to {\mathrm{HNO}}_{\mathrm{3}}+M,\\ {\displaystyle }& {\displaystyle }{\mathrm{NO}}_{\mathrm{2}}+\mathrm{OH}+M\to \mathrm{HOONO}+M.\end{array}$$ For the OH concentrations that are considered in this region (1–$\mathrm{20}\times {\mathrm{10}}^{\mathrm{6}}$ molec. cm${}^{-\mathrm{3}}$), the reactions above follow first-order kinetics. The total sink term can therefore be calculated as $s_{{\mathrm{NO}}_{\mathrm{2}}}={\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}/\mathit{\tau }$ with 2$$\mathit{\tau }={\displaystyle \frac{\mathrm{1}}{k_{\mathrm{mean}}(T,[M\left]\right)\cdot \left[\mathrm{OH}\right]}}.$$ $\mathit{\tau }$ appears here as the characteristic mixed lifetime of ${\mathrm{NO}}_{\mathrm{2}}$ in the atmosphere. The reaction rate $k_{\mathrm{mean}}$ characterises the reactions between ${\mathrm{NO}}_{\mathrm{2}}$ and OH and depends on atmospheric conditions. provide a general expression of this rate as a function of both temperature $T$ and total air concentration [$M$]. Note that HOONO can be rapidly decomposed back to ${\mathrm{NO}}_{\mathrm{2}}$ and OH in the lower troposphere. We assume here that this decomposition is slow and does not affect the ${\mathrm{NO}}_{\mathrm{2}}$ horizontal gradients. Both pathways are therefore taken into account, and the value of $k_{\mathrm{mean}}$ represents the total loss of ${\mathrm{NO}}_{\mathrm{2}}$ due to OH, with a contribution of the HOONO-forming reaction between 5 % and 15 % under atmospheric conditions . Thus, the ${\mathrm{NO}}_{\mathrm{2}}$ production can be calculated as the sum of a transport term and a sink term: 3$$e_{{\mathrm{NO}}_{\mathrm{2}}}=\mathrm{div}\left({\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}\mathit{w}\right)+{\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}/\mathit{\tau }.$$

The treatment for ${\mathrm{NO}}_x$ removal is simplified here. ${\mathrm{NO}}_x$ concentrations are influenced by other sinks. showed that the reaction between ${\mathrm{NO}}_{\mathrm{2}}$ and OH forming ${\mathrm{HNO}}_{\mathrm{3}}$ accounted for most of total ${\mathrm{NO}}_x$ loss at the global scale, but with high uncertainties associated with other sinks. Here, the features of the climate in Egypt during daytime hinder many processes to have a significant effect. The following ${\mathrm{NO}}_x$ sinks, which can be of notable importance at the global scale, are not taken into account here.

• NO${}_{\mathrm{2}}$ deposition through the leaf stomata of vegetation. This sink can be significant in forested areas. In Egypt, the leaf area index is very low, except in the croplands of the Nile Delta where it is comparable to that of southern Europe or the western United States , for which the corresponding lifetime was about 10–100 h , i.e. about an order of magnitude larger than the lifetimes calculated here. To our knowledge, there are no studies focusing on the corresponding lifetimes for croplands, and we therefore do not take this sink into account.

• NO${}_x$ oxidation by organic radicals to produce alkyl and multifunctional nitrates . This sink increases with the concentration of volatile organic compounds (VOCs), whose presence cannot be excluded in Egypt. Different models have estimated low biogenic isoprene emissions in the region . These emissions are concentrated around the Nile River and its delta and do not exceed 15 mg m${}^{-\mathrm{2}}$ d${}^{-\mathrm{1}}$. They are certainly noticeable and higher in summer than in winter, and contrast with the rest of the country, but they remain low compared to other regions in the world. They are, for instance, about an order of magnitude lower than in the forested areas of the eastern US, where the corresponding sink accounts for between 30 % and 60 % of the total ${\mathrm{NO}}_x$ sink . Furthermore, at large ${\mathrm{NO}}_{\mathrm{2}}$ concentrations (compared to VOC concentrations), the share of this sink in the total ${\mathrm{NO}}_x$ loss is weakened compared to that of ${\mathrm{HNO}}_{\mathrm{3}}$ . The effect of biogenic emissions of VOC can therefore be considered minor. However, VOC emissions can also be of anthropogenic origin, especially in urban areas, where they are difficult to estimate. To our knowledge, there is no study evaluating the competition of the two sinks in Egypt or in a region with similar features, and we therefore do not account for this reaction in our calculations.

• NO reaction with HO${}_{\mathrm{2}}$ to produce nitric acid , whose yield is strongly enhanced in presence of water vapour . Here, we neglect this reaction as the corresponding reaction rate is lower by a factor of 3 to 8 in dry conditions .

• NO conversion to NO${}_{\mathrm{3}}$, the latter being in thermal equilibrium with NO${}_{\mathrm{2}}$ and N${}_{\mathrm{2}}{O}_{\mathrm{5}}$. This sink, which takes place via heterogeneous processes, has a significant contribution during nighttime in the Mediterranean region and is neglected at 13:30 LT when OH is close to its daily maximum.

• NO${}_{\mathrm{2}}$ reversible reaction with peroxyacetyl radical to produce peroxyacetylnitrate . In the Nile Delta region, PAN concentrations in the lower troposphere are significantly below the global average , possibly due to high temperatures favouring short PAN lifetimes. Moreover, its production peaks in the late afternoon and early evening . We therefore do not consider this sink in the representation of ${\mathrm{NO}}_x$ emissions at 13:30 LT.

• NO${}_{\mathrm{2}}$ uptake onto black carbon particles . This uptake is of a limited amount in the Mediterranean region .

With all these processes not being accounted for, the reaction between ${\mathrm{NO}}_{\mathrm{2}}$ and OH is the only sink that is considered in our calculations to provide an indication of ${\mathrm{NO}}_x$ emissions. Section 4.7 details the consequences of not considering these various minor sinks on the results.

3.2 Interpolation to daily average emissions

All parameters are evaluated at 13:30 LT, which means that the ${\mathrm{NO}}_{\mathrm{2}}$ production is calculated at the same moment. In Egypt, the maximum and minimum electricity consumption is reached around 20:00 and 06:00 LT respectively, and inter-daily consumption differences have been weakened by the increasing sales of air conditioning and ventilation systems in the past decades . The daily load profiles provided by the National Egyptian Electricity Holding Company show that the mean daily electricity consumption corresponds approximately to the consumption at 13:30 LT in the country . The difference between the two quantities being small both in summer (about $+\mathrm{2}$ % to $-\mathrm{3}$ %) and winter (about $-\mathrm{2}$ % to $-\mathrm{6}$ %), we consider our inferred emissions to be representative of the average activity in Egypt at any time of the year. This assumes that electricity consumption dominates the emissions of the country, or that the other emitting sectors have a similar daily profile. This can be justified. According to CAMS-GLOB-ANT_v4.2, the power sector accounts for 50 % to 60 % of total ${\mathrm{NO}}_x$ emissions in Egypt. EDGARv5.0 presents a lower share (40 % to 45 % of total emissions). Moreover, for both inventories, the transport sector accounts for the majority of the remaining emissions. According to the traffic congestion index in Cairo (https://www.tomtom.com/en_gb/traffic-index/cairo-traffic/, last access: 21 April 2021), the congestion level around 13:30 LT seems to be slightly higher than during the morning peak, but lower than during the night peak. Traffic emissions at this moment of the day could therefore be representative of the average traffic emissions as well.

3.3 Validation of CAMS OH concentration using line density calculations for Riyadh

When the transport term is integrated over large spatial scales, it cancels out due to the mass balance in the continuity equation between ${\mathrm{NO}}_{\mathrm{2}}$ sources and ${\mathrm{NO}}_{\mathrm{2}}$ sinks. At first order, the integration of the inferred emissions over the whole domain (of about 490 000 km${}^{\mathrm{2}}$) thus reflects chemical losses of the sink term. In this term, the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime calculation involves the reaction rate $k_{\mathrm{mean}}$, whose annual variability is low due to small changes in Egyptian midday temperatures throughout the year, and OH concentration, whose annual variability is highly marked. In Egypt, tropospheric OH concentrations, which are strongly correlated with solar ultraviolet radiation and ${\mathrm{NO}}_x$ emissions, are higher in summer than in winter. To ensure an appropriate representation of the OH field by CAMS data, we select a large number of TROPOMI images characterised by a homogeneous wind field, in which we calculate the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime according to Eq. (), where [OH] corresponds to the near-real-time CAMS data and $k_{\mathrm{mean}}$ is calculated with the formula from . We compare this value with the lifetime determined by a method initially developed by , and expanded by by introducing a rotation of the image to refine the chemical lifetime. This method consists in fitting an exponentially modified Gaussian function (EMG) to ${\mathrm{NO}}_{\mathrm{2}}$ line density profiles. These profiles correspond to the integrated ${\mathrm{NO}}_{\mathrm{2}}$ columns along the mean wind direction in the pollution pattern and centred around the source. They are obtained by rotating TROPOMI images in the mean wind direction and using the values of the nearest columns in a 100 km${}^{\mathrm{2}}$ area. Line density profiles are generated on a span of 300 km. An example is given in Fig. . Within the average profile, the ${\mathrm{NO}}_{\mathrm{2}}$ burden and lifetime can be derived from the parameters that describe the best statistical fit. The EMG model is expressed as follows :

4$$\begin{array}{rl}⟨{\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}⟩& \left(x|B,A,x_{\mathrm{0}},\mathit{\mu },\mathit{\sigma }\right)=B+{\displaystyle \frac{A}{\mathrm{2}{x}_{\mathrm{0}}}}\exp \left({\displaystyle \frac{\mathit{\mu }-x}{x_{\mathrm{0}}}}+{\displaystyle \frac{{\mathit{\sigma }}^{\mathrm{2}}}{\mathrm{2}{x}_{\mathrm{0}}^{\mathrm{2}}}}\right)\\ & \mathrm{erfc}\left(-{\displaystyle \frac{\mathrm{1}}{\sqrt{\mathrm{2}}}}\left({\displaystyle \frac{x-\mathit{\mu }}{\mathit{\sigma }}}-{\displaystyle \frac{\mathit{\sigma }}{x_{\mathrm{0}}}}\right)\right).\end{array}$$ Here, $x$ is the distance in the downwind–upwind direction; $B$ is the ${\mathrm{NO}}_{\mathrm{2}}$ background; $A$ is the total number of ${\mathrm{NO}}_{\mathrm{2}}$ molecules observed in the vicinity of the point source; $x_{\mathrm{0}}$ is the $e$-folding distance downwind, representing the exponential length scale of ${\mathrm{NO}}_{\mathrm{2}}$ decay; $\mathit{\mu }$ is the location of the apparent source relative to the centre of the point source; and $\mathit{\sigma }$ is the standard deviation of the Gaussian function, representing the length scale of Gaussian smoothing. Using a non-linear least-squares fit, we estimate the five unknown parameters: $A$, $B$, $x_{\mathrm{0}}$, $\mathit{\mu }$ and $\mathit{\sigma }$. From the mean wind module $w_{\mathrm{mean}}$ in the domain, the mean effective ${\mathrm{NO}}_{\mathrm{2}}$ lifetime ${\mathit{\tau }}_{\mathrm{fit}}$ can be estimated using the fitted parameters: 5$${\mathit{\tau }}_{\mathrm{fit}}={\displaystyle \frac{x_{\mathrm{0}}}{w_{\mathrm{mean}}}}.$$ The geography of Egypt does not suit the method described here. The Egyptian population is contiguously concentrated along the Nile, which makes it difficult to define point sources isolated from human activity. Furthermore, large isolated cities such as Alexandria or Suez are too close to the coast for the wind to be considered homogeneous. We therefore use the city of Riyadh, Saudi Arabia (24.684${}^{\circ }$ N, 46.720${}^{\circ }$ E), to perform the comparison between the CAMS-induced lifetime and the fit-obtained lifetime. Riyadh has been the focus of anterior studies and is particularly suitable for several reasons. Firstly, Riyadh is a city within the latitudinal extent of Egypt (1600 km from Cairo) with a climate that is similar to the typical Egyptian climate. Secondly, ${\mathrm{NO}}_{\mathrm{2}}$ tropospheric columns over Riyadh are high ($\sim \mathrm{9}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$), leading to retrievals with a high signal-to-noise ratio. Thirdly, Riyadh is far from the coast, and its flat terrain makes the surrounding wind fields rather homogeneous during most of the year.

As the fitting algorithm is very sensitive to any disturbance that might be induced by ${\mathrm{NO}}_{\mathrm{2}}$ production from other point sources, it is necessary to identify heavy industrial facilities in the area. Riyadh is also an industrial area, with several power plants located close to the city centre. Figure  shows the location of the most important emitters in the region, which include five gas-fired power plants (PP7, PP8, PP9, PP10 and PP14), one oil-fired power plant (PP4) and one cement plant (CP). The five gas power plants, with a total capacity of more than 10 GW, are located in the periphery of the city. These power plants are sufficiently far away from the city centre for TROPOMI to distinguish their own emissions from those of Riyadh’s centre with a resolution of $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$, which is not the case for CP and PP4, which are located in the city centre. It is therefore appropriate to restrict the study of ${\mathrm{NO}}_{\mathrm{2}}$ patterns over Riyadh to days for which the emissions from the city centre and those from the gas power plants do not overlap. This is the case when the wind blows steadily and homogeneously in a north–south direction. Within about 150 km around the city centre, we thus calculate the average wind given by ERA5 and consider the observation to be reliable if the mean angle $⟨\mathit{\theta }⟩$ of the observations deviates by less than 40${}^{\circ }$ from the north or the south, with a standard deviation ${\mathit{\sigma }}_{\mathit{\theta }}$ of less than 36${}^{\circ }$. This condition generally leads to a selection of observations with large wind speeds, low winds speeds being often associated with more variable directions. This ensures the ${\mathrm{NO}}_{\mathrm{2}}$ decay to be dominated by chemical removal and not by the variability of the wind . Finally, we select observations with clear-sky conditions downstream of the flow (with at least 80 % downstream cells with $q_{\mathrm{a}}>\mathrm{0.75}$).

Figure 2

Map of Riyadh's city centre with the surrounding power plants (PP4, 7, 8, 9, 10 and 14) and cement plant (CP). The map has been extracted from © Google Maps.

[Figure omitted. See PDF]

Our $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ gridding ensures that retrieved lifetimes are governed by physical decay of ${\mathrm{NO}}_{\mathrm{2}}$ and not an artefact of the spatial resolution . The fitting procedure is very sensitive to the wind direction. Instead of manually correcting the ERA5 wind field for individual ${\mathrm{NO}}_{\mathrm{2}}$ patterns, the curve fitting is performed for every sample with three different rotation angles, corresponding to the wind direction with a correction of $-\mathrm{10}$, 0 or 10${}^{\circ }$. A record is kept if one of these three fits leads to a correlation with the corresponding ${\mathrm{NO}}_{\mathrm{2}}$ line density whose coefficient is greater than 0.97. Among the remaining samples, we keep those with a value of ${\mathit{\tau }}_{\mathrm{fit}}$ greater than 1.0 h (considered sufficiently high to be relevant). An example of curve fitting is given in Fig. .

Figure 3

Estimation of the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime from a pattern above Riyadh on 11 March 2020: (a) ${\mathrm{NO}}_{\mathrm{2}}$ plume rotated with its wind direction around the source (star) to an upwind–downwind pattern. Grey boxes centred around black points indicate the extent of the spatial integration of ${\mathrm{NO}}_{\mathrm{2}}$ columns to obtain the ${\mathrm{NO}}_{\mathrm{2}}$ line density. Values of cardinal points are noted in black. (b) Corresponding line densities (points) representing the downwind evolution of ${\mathrm{NO}}_{\mathrm{2}}$ as a function of the distance to Riyadh's city centre and the corresponding fit according to the exponentially modified Gaussian function (dashed line).

[Figure omitted. See PDF]

The phenomena under study here take place in the planetary boundary layer (PBL), which in this region has a midday height of about 2 km . TROPOMI observations only provide information on the total ${\mathrm{NO}}_{\mathrm{2}}$ content of the tropospheric column, without providing information on the vertical distribution of concentrations. Extracting emissions from concentrations therefore requires a selection on the height at which wind, temperature and OH data are taken. and conducted similar studies using the boundary layer average wind, while chose a vertical level of about 450 m above ground. Because vertical transport of ${\mathrm{NO}}_x$, which is emitted mainly from combustion engines and industrial stacks, is generally minor compared to horizontal transport, ${\mathrm{NO}}_x$ is confined to the first few hundred metres above ground level. Using PBL-averaged data poses a problem of consistency as wind, temperature and OH concentration values significantly vary within the PBL. As a consequence, we compare the CAMS-induced lifetime $\mathit{\tau }$ and the fit-induced lifetime ${\mathit{\tau }}_{\mathrm{fit}}$ using the parameters ($\mathit{w}$, [OH] and $T$) at two different vertical levels: a medium level $\mathcal{A}$ at 925 hPa (about 770 m a.g.l. – above ground level) and a bottom level $\mathcal{B}$ at 987.5 hPa (about 210 m). These levels are interpolated from four and two ECMWF or CAMS consecutive pressure levels respectively (1000–850 hPa for level $\mathcal{A}$ and 1000–975 hPa for level $\mathcal{B}$). Since most mask cells have an altitude between 0 and 150 m, the corresponding pressure variations are small (up to $\sim \mathrm{16}$ hPa), which allows us to neglect the effects of topography on the position of pressure levels. Figure  sums up the selection method for the comparison of lifetimes.

Figure 4

Selection method for ${\mathrm{NO}}_{\mathrm{2}}$ patterns over Riyadh. Datasets (yellow-orange) are used to calculate the quantities (light green) that are submitted to different tests (grey). A total of 731 patterns are progressively conserved (green arrows) or rejected (red arrows). At each stage, the number of conserved or rejected patterns are noted within brackets (the value is only given for calculations performed at level $\mathcal{B}$). This selection process compares the lifetimes estimated by the EMG function fitting with TROPOMI line density profiles to the lifetimes calculated according to Eq. () with CAMS data.

[Figure omitted. See PDF]

3.4

Calculation of anthropogenic ${\mathrm{NO}}_x$ emissions and comparison with inventories

We calculate ${\mathrm{NO}}_x$ emissions on the entire domain from ${\mathrm{NO}}_{\mathrm{2}}$ production by using CAMS NO and ${\mathrm{NO}}_{\mathrm{2}}$ concentrations. These are not intended to replace TROPOMI observations; they are used to apply the concentration ratio $\left[{\mathrm{NO}}_x\right]/\left[{\mathrm{NO}}_{\mathrm{2}}\right]=\left(\right[\mathrm{NO}]+[{\mathrm{NO}}_{\mathrm{2}}\left]\right)/\left[{\mathrm{NO}}_{\mathrm{2}}\right]$ to account for the conversion of ${\mathrm{NO}}_{\mathrm{2}}$ to NO and vice versa. As diurnal NO concentrations in urban areas generally range from 10 to 150 ppb , the characteristic stabilisation time of this ratio never exceeds a few minutes . Since this time is lower than the order of magnitude of the inter-mesh transport time (about 30 min considering the resolution used and the mean wind module in the region), we can reasonably neglect the effect of the stabilisation time of the conversion factor on the total composition of the emissions and treat each cell of the grid independently from its neighbours. found an annual average of 1.32 for this conversion factor, but CAMS data show small deviations from this value over Egyptian urban areas. We therefore calculate ${\mathrm{NO}}_x$ emissions for each cell of the domain as follows: 6$$e_{{\mathrm{NO}}_x}={\displaystyle \frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}}{e}_{{\mathrm{NO}}_{\mathrm{2}}}.$$

For convenience, quantities $\frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}\mathrm{div}\left({\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}\mathit{w}\right)$ and $\frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}{\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}/\mathit{\tau }$ represent the respective contributions of the transport and the sink terms to total ${\mathrm{NO}}_x$ emissions. We finally obtain the emissions related to human activity $E_{{\mathrm{NO}}_x}$ by removing the arithmetic mean value of ${\mathrm{NO}}_x$ emissions above background cells from total emissions: 7$$E_{{\mathrm{NO}}_x}=e_{{\mathrm{NO}}_x}-{\displaystyle \frac{\mathrm{1}}{n_{\mathrm{b}}}}\sum _{i=\mathrm{1}}^{n_{\mathrm{b}}}{e}_{{\mathrm{NO}}_x,i}.$$

These removed emissions are linked to the ${\mathrm{NO}}_{\mathrm{2}}$ background estimated by TROPOMI. This background, which is mostly located in the upper troposphere, is inconsistent with the use of other parameters which are calculated in the lower troposphere. As such, these emissions do not correspond to anthropogenic emissions, but they provide the value of what must be subtracted from the estimates to obtain emissions related to human activity. Such a removal assumes that other processes involved in ${\mathrm{NO}}_x$ budgets lead to similar emissions inside and outside the mask, which is not evident, as the majority of background cells are located in the desert or the ocean while the majority of mask cells are located near the Nile River. However, as the processes involved in natural ${\mathrm{NO}}_x$ sources lead to emissions much smaller than anthropogenic emissions in polluted areas, we neglect this difference in the following calculations. An alternative would be to calculate an average concentration for the background cells and subtract the corresponding value from the column densities before calculating the emissions. This would pose further reliability problems: for instance, very high ${\mathrm{NO}}_{\mathrm{2}}$ concentrations could appear outside the mask due to wind transport (an example is shown in Fig. ). They would lead to an overestimation of the ${\mathrm{NO}}_{\mathrm{2}}$ background and thus to an underestimation of the anthropogenic emissions.

Neglecting the part of the country that lies outside the domain, total emissions from the anthropogenic activity of Egypt can then be obtained by integrating $E_{{\mathrm{NO}}_x}$ on the whole domain. For robust statistics, these derived emissions can be averaged monthly, enabling a month-by-month comparison with bottom-up inventories. The linearity of the averaging processes ensures the interchangeability of temporal and spatial averages. A monthly average is relevant because it aggregates enough data to limit the impact of the outliers due to uncertainties in wind and OH representation. In addition, it enables the study of monthly ${\mathrm{NO}}_x$ emission profiles, which reflect changes in human activities throughout the year due to temperature changes, economic constraints and cultural norms.

4 Results and discussion

4.1

Line densities and ${\mathrm{NO}}_{\mathrm{2}}$ lifetime

We compare the results of the TROPOMI line density fits for Riyadh to the lifetime calculated by Eq. () using CAMS OH data. The 2 years of TROPOMI observations (from November 2018 to November 2020) provide a wide variety of profiles. For level $\mathcal{B}$, Fig.  also provides the number of samples that are being kept at each stage of the process. Of the 731 observations available, 203 have a wind direction in the cone with a north–south orientation with an aperture of 40${}^{\circ }$ (i.e. between 340 and 20${}^{\circ }$ or between 160 and 200${}^{\circ }$). Of the remaining observations, 166 occurred with a sufficiently clear sky to be retained. The criterion of weak variability for the wind direction brings to 91 the number of observations that are kept by the method. With these 91 observations, the line density profiles are calculated and the fits applied. According to Eq. (), the lifetime is calculated using the mean wind module around the point source. The two lifetimes are calculated with the parameters taken at the medium level $\mathcal{A}$ or at the top level $\mathcal{B}$. Of the 91 fits obtained, 51 are of high quality (correlation coefficient between fit function and line density profile greater than 0.97) for level $\mathcal{A}$ and 52 for level $\mathcal{B}$. A total of 39 of these fits lead to a lifetime ${\mathit{\tau }}_{\mathrm{fit}}$ greater than 1.0 h for level $\mathcal{A}$ and 41 for level $\mathcal{B}$. All remaining samples correspond to atmospheric conditions with moderate to fast winds, with a module ranging between 2 and 11 m s${}^{-\mathrm{1}}$ (with an average of 5.9 m s${}^{-\mathrm{1}}$) for level $\mathcal{A}$ and between 3 and 8 m s${}^{-\mathrm{1}}$ (with an average of 5.4 m s${}^{-\mathrm{1}}$) for level $\mathcal{B}$. These lifetimes are compared to the corresponding lifetimes obtained from CAMS data in Fig. , which is divided into seasons for a more convenient comparison. The use of level $\mathcal{A}$ leads to notable underestimations of the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime in autumn compared to the lifetime calculated with the fitting method. This same level leads to lifetime overestimations in winter. This trend is not found with the use of level $\mathcal{B}$, which leads to a better reproduction of the lifetimes calculated with the fitting method for the available seasons. Figure  shows a linear regression between the two calculated lifetimes. The results are scattered, with a correlation coefficient higher for level $\mathcal{B}$ (0.408) than for level $\mathcal{A}$ (0.220). When the intercept of the regression line is forced to zero, the resulting slope is closer to 1 for level $\mathcal{B}$ (0.998) than for level $\mathcal{A}$ (1.071).

Figure 5

(a, c) Comparison between CAMS-derived ${\mathrm{NO}}_{\mathrm{2}}$ lifetimes and lifetimes from ${\mathrm{NO}}_{\mathrm{2}}$ line density fittings with EMG function above Riyadh city centre, for levels $\mathcal{A}$ (a) and $\mathcal{B}$ (c). The samples presented correspond to patterns in clear-sky conditions for which the mean wind is in the north–south direction with a low variance, and for which the correlation between line density profile and fit gives a correlation coefficient of more than 0.97 and a lifetime of more than 1.0 h. ${\mathrm{NO}}_{\mathrm{2}}$ patterns do not have these conditions during the summers of 2019 and 2020. Dashed lines separate the groups of observations by season. (b, d) Comparison between the two calculated lifetimes for levels $\mathcal{A}$ (b) and $\mathcal{B}$ (d). A linear regression with an intercept forced to be zero is displayed with a green dashed line.

[Figure omitted. See PDF]

Although both correlations are weak, level $\mathcal{B}$ leads to a better match between the lifetime calculated with Eq. () and the lifetime calculated from line densities. The results that are presented in the following sections (except for Sect. 4.3) are therefore results of calculations performed with parameters ($\mathit{w}$, [OH], $T$ and $\left[{\mathrm{NO}}_x\right]/{\mathrm{NO}}_{\mathrm{2}}]$) estimated at level $\mathcal{B}$. Nevertheless, it should be noted that no summer observations were included in the comparison. The main reason for this is the wind direction: of the 188 summer days observed, 178 have a mean wind direction outside the north–south cone over central Riyadh. On the remaining 10 d (one for summer 2019 and nine for summer 2020), the ERA5 wind direction is too variable for the fit to be considered relevant, or the fit results in a correlation coefficient below 0.97. It is not clear how correctly the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime would be calculated during both summer periods by Eq. (). With OH concentrations being the main driver of this lifetime, we cannot assess the relevance of the representation of OH concentrations by CAMS data during summer days in the study.

4.2

Mapping of Egypt's ${\mathrm{NO}}_x$ emissions

As a first step, we try to map ${\mathrm{NO}}_x$ emissions in Riyadh using parameters estimated at level $\mathcal{B}$. For the period from December 2017 to October 2018 and using a constant lifetime of 4 h, estimated at 6.66 kg s${}^{-\mathrm{1}}$ the emissions of the corresponding urban area and a mean rate density of about 3.7 $\mathrm{\mu }$g m${}^{-\mathrm{2}}$ s${}^{-\mathrm{1}}$ for power plants PP9 and PP10–14, with the transport term accounting for about 80 % to 90 % of this budget. Using the same domain for December 2018 to October 2019 with our method, we found a mean lifetime of 2.94 h and mean emissions of 5.92 kg s${}^{-\mathrm{1}}$ for the urban area. We also found smaller rate densities for the power plants (about 3.4 $\mathrm{\mu }$g m${}^{-\mathrm{2}}$ s${}^{-\mathrm{1}}$ for PP9 and 3.0 $\mathrm{\mu }$g m${}^{-\mathrm{2}}$ s${}^{-\mathrm{1}}$ for PP10–14), with a smaller contribution of the transport term (about 70 %). Despite differences in resolution, AMF calculation, lifetime variability and background removal, the two methods give similar results.

Figure 6

${\mathrm{NO}}_x$ emissions above most of Egypt's territory: (a–d) transport term (a), sink term (b), resulting emissions (c) and the corresponding anthropogenic emissions after non-anthropogenic background removal (d) for January 2019. (e–h) Transport term (e), sink term (f), resulting emissions (g) and the corresponding anthropogenic emissions after background removal (h) for July 2019.

[Figure omitted. See PDF]

The top-down emission model is then applied to the Egyptian domain with CAMS OH concentration and temperature fields for lifetime calculations. For each cell, ${\mathrm{NO}}_x$ emissions are calculated according to Eq. (), resulting in a mapping of Egypt's emissions. The obtained values are averaged monthly from November 2018 to November 2020. Figure  shows a composition of the emissions map with the transport term and the sink term for the months of January and July 2019. The corresponding anthropogenic emissions, calculated according to Eq. (), are added. The Nile appears on transport term maps: the divergence calculation complies with what is expected from a line of emitters, i.e. a clear separation of zones of positive divergence from zones of negative divergence with a separation line corresponding to the course of the river. The fact that areas of negative and positive divergence are respectively located to the east and the west of the river is the result of the zonal component of the wind being positive most of the time. Some point sources like Cairo, Alexandria, Asyut or Aswan are easily identifiable. The sink term, directly proportional to the TROPOMI column densities, also highlights these cities. However, unlike the transport term, which has a similar spatial pattern from month to month, the sink term is clearly stronger in summer than in winter. This is mainly due to a higher lifetime in winter than in summer (4.94 h on average in January 2019 and 2.62 h in July 2019) while the average TROPOMI ${\mathrm{NO}}_{\mathrm{2}}$ concentrations are slightly higher during winter ($\mathrm{1.071}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ for January 2019 and $\mathrm{0.899}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ for July 2019). Over the whole domain, the mean transport term varies throughout the studied period between $-\mathrm{0.014}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$ (December 2019) and $\mathrm{0.015}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$ (May 2020). Thus, it hardly contributes to the ${\mathrm{NO}}_x$ emission budget, with the mean chemical sink term alone varying between $\mathrm{0.223}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$ (December 2019) and $\mathrm{0.534}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$ (September 2020).

Table 1

Comparison between the transport term and the sink term above different cities in Egypt, as well as the industrial region of Ain Sokhna located 45 km southwest of Suez for January and July 2019. TROPOMI vertical ${\mathrm{NO}}_{\mathrm{2}}$ columns, ${\mathrm{NO}}_x$ emissions and population densities correspond to average values within 18 km from the city centre. Unit $\mathcal{M}$ stands for a quantity of 10${}^{\mathrm{15}}$ molecules (${\mathrm{NO}}_{\mathrm{2}}$ or ${\mathrm{NO}}_x$).

City	Population density	Jan 2019				Jul 2019
	(inhabitants per	${\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}$	Transport	Sink		${\mathrm{\Omega }}_{{\mathrm{NO}}_{\mathrm{2}}}$	Transport	Sink
	square kilometre)	(${\mathcal{M}}_{{\mathrm{NO}}_{\mathrm{2}}}$ cm${}^{-\mathrm{2}}$)	(${\mathcal{M}}_{{\mathrm{NO}}_x}$ cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$)			(${\mathcal{M}}_{{\mathrm{NO}}_{\mathrm{2}}}$ cm${}^{-\mathrm{2}}$)	(${\mathcal{M}}_{{\mathrm{NO}}_x}$ cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$)
Cairo	18 064	9.415	2.903	3.684		5.618	2.022	4.879
Alexandria	9133	3.034	1.179	0.975		1.674	0.410	1.421
Asyut	1644	1.708	0.679	0.718		2.137	1.236	1.520
Aswan	319	0.976	0.182	0.473		0.871	0.308	0.523
Beni Suef	2056	2.950	0.548	1.080		2.321	0.428	1.591
Ain Sokhna	5	3.133	1.256	1.115		2.561	1.346	1.757

Figure 7

Share of the absolute value of the transport term in the sum of the sink term and the absolute value of the transport term above most of Egypt's territory, indicating the local importance of the transport term in ${\mathrm{NO}}_x$ emissions above mask cells. The average of this ratio is shown for January 2019 (a) and July 2019 (b).

[Figure omitted. See PDF]

Several cities in the country appear as the main emitters of the country, such as Cairo, Alexandria, Beni Suef, Asyut or Aswan. The industrial area of Ain Sokhna, located southwest of Suez, also appears as a main emitter. Table  compares the monthly values for the sink term and the absolute value of the transport term above five major cities of the country, with populations ranging from 200 000 to 20 million inhabitants, as well as Ain Sokhna's area. The mean values for TROPOMI column densities are also provided. According to the results, the capital city of Cairo is by far the largest emitter in the country, largely due to its large population, resulting in high traffic emissions, but also to its intensive industrial activity. Alexandria, the country's second largest city, is not necessarily the second largest emitter, as its emissions are comparable to those of smaller cities such as Beni Suef or Asyut. However, the three cities concentrate a large amount of industrial activity: Alexandria hosts several oil and gas power plants and a small number of cement factories, while Beni Suef is close to several oil and gas power plants and hosts several flaring sites. Similarly, the city centre of Asyut is close to three oil and gas-fired power plants and a cement factory. This seems to indicate that industrial activity might be the main cause of ${\mathrm{NO}}_x$ emissions differences between these cities, before population size. This explains why ${\mathrm{NO}}_x$ emissions from these three cities are comparable to those of the industrial area of Ain Sokhna, which includes several cement facilities, iron smelters and oil and gas plants. It might also explain why Aswan, which has a population that is comparable to Beni Suef or Asyut, but which does not have any major industrial site, has slightly lower emissions than the two other cities. An additional analysis of the differences between Asyut and Aswan is provided in Sect. 4.6. Finally, the Gulf of Suez displays relatively large emissions, which might be attributed to the shipping sector, the region being a major gateway for international trade. Because it also hosts several flaring sites, these emissions might also be due to the oil and gas extraction activity.

Table 2

Analysis of the effect of a vertical change in the parameters used to estimate the mean sink term in ${\mathrm{NO}}_x$ emissions: temperature, hydroxyl radical concentration and ${\mathrm{NO}}_x:{\mathrm{NO}}_{\mathrm{2}}$ concentration ratio. The comparison is conducted between the estimated quantities at level $\mathcal{B}$ and level $\mathcal{A}$ for mask cells (MASK) and background cells (BKGD) for January, April, October and July 2019. Values within brackets represent the variation from the base case for which all quantities are estimated at level $\mathcal{B}$.

		Sink term (10${}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$)
Level $\mathcal{B}$	Level $\mathcal{A}$	19 Jan	19 Jan	19 Apr	19 Apr	19 Jul	19 Jul	19 Oct	19 Oct
(987.5 hPa)	(925 hPa)	(MASK)	(BKGD)	(MASK)	(BKGD)	(MASK)	(BKGD)	(MASK)	(BKGD)
$T$, $[$OH$]$, $\frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}$	–	0.859	0.253	1.072	0.345	1.125	0.376	0.932	0.277
$[$OH$]$, $\frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}$	$T$	0.899	0.264	1.127	0.361	1.185	0.394	0.887	0.264
		($+\mathrm{4.7}$ %)	($+\mathrm{4.2}$ %)	($+\mathrm{5.2}$ %)	($+\mathrm{4.6}$ %)	($+\mathrm{5.3}$ %)	($+\mathrm{4.9}$ %)	($+\mathrm{4.8}$ %)	($+\mathrm{4.5}$ %)
$T$, $\frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}$	$[$OH$]$	0.769	0.219	1.013	0.324	1.129	0.375	0.853	0.251
		($-\mathrm{10.5}$ %)	($-\mathrm{13.6}$ %)	($-\mathrm{5.5}$ %)	($-\mathrm{6.0}$ %)	($+\mathrm{0.4}$ %)	($-\mathrm{0.3}$ %)	($-\mathrm{8.5}$ %)	($-\mathrm{9.5}$ %)
$T$, $[$OH$]$	$\frac{\left[{\mathrm{NO}}_x\right]}{\left[{\mathrm{NO}}_{\mathrm{2}}\right]}$	0.872	0.257	1.094	0.352	1.143	0.383	0.904	0.271
		($+\mathrm{1.6}$ %)	($+\mathrm{1.4}$ %)	($+\mathrm{2.1}$ %)	($+\mathrm{2.0}$ %)	($+\mathrm{1.6}$ %)	($+\mathrm{1.9}$ %)	($+\mathrm{3.1}$ %)	($+\mathrm{2.2}$ %)

Although these cities and areas can be described as high-emission sites, the terms responsible for these emissions differ from one site to the other. Figure  shows the contribution of the transport term (taken in absolute value) to total emissions for January and July 2019. Because wind fields are relatively homogeneous along the Nile on spatial scales of less than 100 km, ${\mathrm{NO}}_{\mathrm{2}}$ concentration gradients perceived by TROPOMI in the region mainly contribute to the increase in the transport term, which can reach similar values as the sink term. However, it is never significantly higher than the sink term: due to a spread of the emissions over large urban areas, the behaviour of these cities is therefore different from that of a point source for which the transport term would be very high .

Desert areas such as the Libyan Desert, the Eastern Desert and the Sinai region (located respectively to the west, east and northeast of the Nile) show a very low value for the transport term compared to the sink term, due to the homogeneity of both the wind field and the detected ${\mathrm{NO}}_{\mathrm{2}}$ concentrations in these areas. In the case of the Gulf of Suez, the transport term can be 1 to 2 times higher than the sink term, which varies between 0.4 and $\mathrm{1.2}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$. Those values are slightly higher than the average emissions above background cells areas due to the sink term (about $\mathrm{0.2}-\mathrm{0.6}\times {\mathrm{10}}^{\mathrm{15}}$ molec. cm${}^{-\mathrm{2}}$ h${}^{-\mathrm{1}}$), but remain quite low compared to the emissions in large cities. This relative predominance of the transport term is explained by a visible gradient of the TROPOMI ${\mathrm{NO}}_{\mathrm{2}}$ column densities. The region thus acts as a very thin line of emitters. Nevertheless, this predominance might also be partly due to a poor representation of the wind field. The low resolution of ERA5 (about 26 km in this region, which is the same order of magnitude as the width of the channel) may misrepresent the wind near the coast, creating artificial gradients.

4.3 Vertical analysis

Here we investigate the influence of the choice of the vertical level in the representation of the different model parameters. This influence is of considerable importance, as ${\mathrm{NO}}_x$ sources in urban areas can be located at different altitudes. For instance, emissions from the road sector from tailpipes are located at ground level, whereas ${\mathrm{NO}}_x$ from power plants and industrial facilities can be emitted from stacks, which are usually located between 50 and 300 m a.g.l. (above ground level). Section 4.1 results showed that level $\mathcal{B}$ was more appropriate for the representation of the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime. This level is therefore chosen as a reference for the comparison. We study the effect of a transition from level $\mathcal{B}$ to level $\mathcal{A}$ for each of the three parameters involved in the representation of the sink term, i.e. temperature $T$, hydroxyl radical concentration [OH] and concentration ratio $\left[{\mathrm{NO}}_x\right]/\left[{\mathrm{NO}}_{\mathrm{2}}\right]$. The results for the averages over mask cells and background cells are given for the months of January, April, July and October 2019 in Table . As the wind field is only involved in the transport term whose spatial integration nearly leads to zero, the influence of this parameter is not studied.

The transition to level $\mathcal{A}$ generally results in a decrease in temperature, leading to an increase in the reaction rate $k_{\mathrm{mean}}$ and thus an increase in the emissions from the sink term. This transition has only a small influence on the total ${\mathrm{NO}}_x$ emission estimates, with mask and background cell emissions increasing by 4 % to 6 %. The influence of OH goes in the opposite direction: its concentration decreases with altitude, weakening the sink term. This weakening is particularly visible during winter months, for which the emissions are lower by up to 14 %. In summer however, the effect is hardly noticeable. Finally, the influence of the ${\mathrm{NO}}_x:{\mathrm{NO}}_{\mathrm{2}}$ ratio is negligible on the ${\mathrm{NO}}_x$ emission estimates. Thus, the transition to level $\mathcal{A}$ results in an increase in the sink term of 1 % to 4 %, due to a decrease in concentrations of both NO and ${\mathrm{NO}}_{\mathrm{2}}$ with respect to the vertical but with a greater decrease for ${\mathrm{NO}}_{\mathrm{2}}$. This vertical study confirms the crucial importance of the OH concentration for the accurate representation of ${\mathrm{NO}}_x$ emissions. It appears here as an important driver of the sink term, which is much more sensitive to vertical differences than temperature or the ${\mathrm{NO}}_x:{\mathrm{NO}}_{\mathrm{2}}$ concentration ratio.

4.4 Weekly cycle

In Egypt, the official rest day is Friday, and the economic activity of the country is a priori lower during this day than during the other days of the week. We therefore try to characterise this feature, by evaluating the weekly cycle of ${\mathrm{NO}}_x$ emissions. We use the TROPOMI-inferred emissions to obtain averages per day of the week. We use the quality assurance $q_{\mathrm{a}}$ of TROPOMI retrievals to ignore the days for which more than 20% of the domain has low-quality data (this happens 43 times in 2018/19 and 28 times in 2019/20). Such a filtering avoids accounting for the days when a large part of the urban and industrial areas is covered by clouds. However, it misses situations where small clouds are localised over large emitters, in which case the corresponding emissions are under-estimated. Figure  shows the resulting daily emissions for the period November 2018–November 2019 and November 2019–November 2020. ${\mathrm{NO}}_x$ emissions are expressed as NO and in kilotons per day. A Friday minimum is observed, defining a weekly cycle. This trend is also observed for mean ${\mathrm{NO}}_{\mathrm{2}}$ column densities, for which no intra-weekly variation is observed. Over the 2018–2019 period, Fridays have average emissions of $\mathrm{0.811}\pm \mathrm{0.408}$ kt, which is lower than average emissions for the rest of the week, which reach $\mathrm{0.997}\pm \mathrm{0.533}$ kt. A similar trend is observed in 2019–2020, for which the average for Fridays is $\mathrm{0.704}\pm \mathrm{0.357}$ kt and the average for other days is $\mathrm{0.921}\pm \mathrm{0.449}$ kt. The difference in emissions between the two periods is due to smaller emissions in December 2019, January 2020 and February 2020 that are discussed in Sect. 4.5. On average, Friday emissions correspond to a ratio of $\mathrm{0.83}:\mathrm{7}$ (i.e. a value of 0.83 after normalisation on the seven days of the week) for the entire domain. This result is consistent with the values obtained by , who used TROPOMI data and another emission model to calculate a ratio of $\mathrm{0.71}:\mathrm{7}$ for Cairo and $\mathrm{0.89}:\mathrm{7}$ for Alexandria in 2017.

Figure 8

Weekly profiles of anthropogenic ${\mathrm{NO}}_x$ emissions for Egypt using TROPOMI observations in 2018–2019 (purple line) and 2019–2020 (green line). Thin grey lines represent individual weeks. Days for which less than 80 % of the domain counts low-quality observations ($q_{\mathrm{a}}<\mathrm{0.75}$) are not represented.

[Figure omitted. See PDF]

4.5 Impacts of lockdown during COVID-19

The ongoing global outbreak of COVID-19 forced many countries around the world to implement unprecedented public health responses, including travel restrictions, quarantines, curfews and lockdowns. Such measures have helped to counter the spread of the virus and have, meanwhile, caused high reductions in global demand for fossil fuels . They also led to a fall in the levels of ${\mathrm{NO}}_{\mathrm{2}}$ and other air pollutants across the globe . To prevent the spread of COVID-19, Egyptian authorities ordered a partial lockdown from 15 March till 30 June 2020, closing all public areas (e.g. sport centres, nightclubs, restaurants and cafes) and suspending religious activities in all mosques and churches throughout the country. They also implemented more drastic measures such as a full lockdown during Easter (20 April) and Eid (23 to 25 May), before lifting some restrictions on 1 June . In addition to the effect of containment on the activity of the country, the global decline in consumption led to a drop in the production of certain industrial products.

Several studies have estimated the impact of these events on the air pollution levels in the urban centres of the country: from in situ measurements, estimated that ${\mathrm{NO}}_{\mathrm{2}}$ concentrations had dropped by 25.9 % in Alexandria’s city centre after the start of the lockdown on 13 March, while used OMI retrievals to estimate a 45.5 % reduction of ${\mathrm{NO}}_{\mathrm{2}}$ concentrations for the entire country during the spring compared to 2018 and 2019 average values. However, due to a changing lifetime of ${\mathrm{NO}}_{\mathrm{2}}$, reductions in the concentrations of ${\mathrm{NO}}_{\mathrm{2}}$ might not be entirely due to a decrease in ${\mathrm{NO}}_x$ emissions, which leads us to focus on the variation in ${\mathrm{NO}}_x$ emissions during this singular period. Using our top-down emission model, reductions in total ${\mathrm{NO}}_x$ emissions of 20.1 %, 11.8 % and 13.5 % are observed for the months of March, April and May 2020 compared to the equivalent months in 2019. This drop of emissions in 2020 compared to 2019 calculated by the model also correspond to a decrease in observed ${\mathrm{NO}}_{\mathrm{2}}$ columns: TROPOMI retrievals above mask cells show a decrease in ${\mathrm{NO}}_{\mathrm{2}}$ column densities of 21.6 % over the same period. However, these effects observed for the months of March, April and May 2020 are not repeated in June 2020, for which emissions show an increase of 15.8 % compared to June 2019. This rise is largely the result of an increase in the difference between average estimates inside and outside the mask. Indeed, emissions within the mask in June 2020 are higher than those of June 2019, due to an increase in TROPOMI concentrations above mask cells ($+\mathrm{7.7}$ %) while the ${\mathrm{NO}}_{\mathrm{2}}$ lifetime is almost unchanged ($+\mathrm{3.3}$ %). Emissions outside the mask vary in the opposite direction: a decrease in TROPOMI background concentrations ($-\mathrm{5.4}$ %) is observed while ${\mathrm{NO}}_{\mathrm{2}}$ lifetime increases strongly ($+\mathrm{16.0}$ %). This increase in June emissions seems to indicate that the lift on restrictions allowed a catch-up of the economic activity which was sufficiently strong to generate higher emissions in 2020 than in 2019. Note that CAMS OH concentrations during the lockdown periods do not show significant variations from previous years, although concentration values are slightly lower in 2020 than in 2019 (about 5.5 % lower over the mask cells for the period of March, April and May). The near-real-time CAMS system did not take into account the decrease in anthropogenic emissions in the representation of its OH concentrations. However, the satellite constraints inherent in the system may have modulated the lockdown effects locally or globally. Given the non-linearity of the chemistry but also given the large reactivity of OH with other species whose concentrations have varied differently during the lockdown, it is difficult to determine how these observations have impacted OH concentrations.

4.6 Annual cycle and comparison to inventories

Here, we attempt to compare our TROPOMI-derived ${\mathrm{NO}}_x$ emissions in Egypt to emissions from CAMS-GLOB-ANT_v4.2 and EDGARv5.0 inventories. Figure  shows the total anthropogenic ${\mathrm{NO}}_x$ emissions over the mask cells from November 2018 to November 2020, calculated according to Eq. (). As indicated in Sect. 3.2, the emissions, calculated at 13:30 LT, are representative of the average daily consumption in Egypt. The total calculated for each month therefore corresponds to the ${\mathrm{NO}}_x$ production by human activities in the country. After aggregating the different sectors of activity, CAMS and EDGAR inventories directly provide the anthropogenic ${\mathrm{NO}}_x$ emissions over the same domain. All ${\mathrm{NO}}_x$ emissions are expressed in mass terms as NO. We note that the EDGAR inventory does not cover the period 2018–2020 (the last available year of the inventory is 2015). In Fig. , EDGAR emissions corresponding to the period between November 2013 and November 2015 are displayed.

Figure 9

(a) Comparison of TROPOMI-derived anthropogenic ${\mathrm{NO}}_x$ emissions in Egypt (light blue), with the corresponding emissions from EDGAR (green with stripes) and CAMS (yellow) inventories. EDGAR data are shown for comparison purposes and cover the years 2013–2015. (b) TROPOMI-derived anthropogenic ${\mathrm{NO}}_x$ emissions for the cities of Asyut (dark red) and Aswan (light blue). The corresponding domains are displayed in Fig. .

[Figure omitted. See PDF]

TROPOMI-derived emissions are higher than the CAMS inventory estimates. The top-down model estimates total ${\mathrm{NO}}_x$ emissions of 697.6 kt over the 24 months, which is 45.9 kt higher than CAMS for the same period (651.6 kt). This difference is primarily localised in the first 12 months, for which TROPOMI-inferred emissions are always higher than the inventories and show higher values in summer than during the rest of the year. The next 12 months show similar emissions in summer but much lower values in winter. In particular, the difference is significant in December 2019 and January 2020 (respectively 56.5 % and 66.5 % of CAMS levels). These emissions also contrast with other winter emissions, with a total of 31.7 kt for December 2019–January 2020 against 53.3 kt for December 2018–January 2019 and 57.7 kt for December 2020–January 2021. In the computations, this drop for winter 2019/20 is mainly due to a relatively low value of the OH concentration, which reaches $\mathrm{4.61}\times {\mathrm{10}}^{\mathrm{6}}$ molec. cm${}^{-\mathrm{3}}$ on average for December 2019 and January 2020, with $\mathrm{4.29}\times {\mathrm{10}}^{\mathrm{6}}$ molec. cm${}^{-\mathrm{3}}$ above mask cells and $\mathrm{4.69}\times {\mathrm{10}}^{\mathrm{6}}$ molec. cm${}^{-\mathrm{3}}$ over background cells. They were respectively 5.29, 5.74 and $\mathrm{5.18}\times {\mathrm{10}}^{\mathrm{6}}$ molec. cm${}^{-\mathrm{3}}$ for the previous year (December 2018–January 2019) and 5.11, 4.90 and $\mathrm{5.16}\times {\mathrm{10}}^{\mathrm{6}}$ molec. cm${}^{-\mathrm{3}}$ for the subsequent year (December 2020–January 2021). A decrease in tropospheric columns ($-\mathrm{18.5}$ % for mask cells and $-\mathrm{7.6}$ % for background cells compared to winter 2018/19) also contributes to this drop. The accuracy of the inferred emissions for winter 2019/20 can therefore be questioned.

At first sight, the annual variability of TROPOMI-inferred emissions, which describes a 1-year cycle with higher emissions in summer, seems to be correlated with power emissions which dominate the use of fossil fuels in Egypt . These power emissions are due to the country's residential electricity consumption . They also meet the needs of industry. Summer peaks in electricity consumption are mostly driven by temperature, as illustrated by the increasing sales of air conditioning and ventilation systems for several decades . The use of air conditioning in cars, which requires an additional amount of fuel, could also contribute to the increase in ${\mathrm{NO}}_x$ emissions in summer. To support this hypothesis, we use our model on two smaller domains centred around the two cities of Asyut and Aswan. The corresponding domains are displayed in Fig. . Both cities have similar demographic features, with populations of about 467 000 and 315 000 inhabitants in 2021 and human densities of about 3000 and 1600 inhabitants per square kilometre respectively. However, their industrial features largely differ. There is no large fossil-fuel-fired power plant in Aswan, where most of the electricity is produced by a hydroelectric dam, whereas Asyut counts three oil and gas power plants of various capacities (90, 650 and 1500 MW) in its urban area. Both cities have a cement plant, but the one in Asyut has a larger production capacity (5.7 Mt yr${}^{-\mathrm{1}}$ in Asyut, 0.8 Mt yr${}^{-\mathrm{1}}$ in Aswan). Our model is used following the same procedure as for the main domain. The background removal is done at the scale of the country. A seasonal cycle appears for Asyut, with a minimum for winter months and a maximum for summer months. This cycle seems slightly shifted from the one observed for the entire country, for which May emissions are as important as those of summer months. We also note that the decrease in emissions for winter 2019/20 is less marked than for the emissions of the whole country, and of a similar value to the previous winter. This suggests that national ${\mathrm{NO}}_x$ emissions are indeed lower during winter, but that the values obtained for winter 2019/20 are particularly low. We also find that the seasonality of the emissions is more pronounced for the Asyut domain than for the country as a whole. The case of Aswan is different. Emissions within the corresponding domain are significantly lower than for Asyut. The signal-to-noise ratio being higher, it is difficult to characterise an annual cycle, but the results do not seem to indicate low emissions in winter and high emissions in summer. This identification of a seasonal cycle identical to that of the entire country for a city with several power plants, and the absence of such a cycle in a city without any, strengthens the hypothesis that the power sector plays a major role in Egyptian ${\mathrm{NO}}_x$ emissions.

We note that some features of the industrial activities in the country might be counteracting this trend. For some sectors such as cement or steel, production is lower in summer, due to the physical wear experienced by workers due to heat, but also due to certain periods of leave. Given the importance of industrial activities in the production of ${\mathrm{NO}}_x$ shown in Sect. 4.2, this aspect cannot be neglected. The transport sector could also counteract the observed trend: although the use of air conditioning in cars increases ${\mathrm{NO}}_x$ emissions of the sector, the observed mean traffic in the country is higher between November and February and lower between June and August, especially in Cairo, which includes most of the population. In the absence of additional data, it therefore seems difficult to conclude on the amplitude of the seasonal cycle produced by our top-down model. This caution is all the more necessary as CAMS and EDGAR show seasonal cycles for ${\mathrm{NO}}_x$ emissions, with different dynamics than those displayed by TROPOMI emissions: while the EDGAR inventory predicts a maximum of emissions in December or January and a minimum in April, the CAMS inventory shows two local maxima each year in May and November and two local minima in February and September. The amplitude of the corresponding cycles is much lower in those inventories, representing 14.2 % of the average value for emissions estimates for EDGAR and 12.4% for CAMS. Those values must be compared to the amplitude displayed by TROPOMI-inferred emissions, for which the maximum-to-minimum ratio is about 1.8 if winter 2019/20 is excluded and 2.7 if it is included.

4.7 Uncertainties and assessments of results

The estimation of ${\mathrm{NO}}_x$ emissions is based on the use of several quantities with varying uncertainties. The error bars shown in Figs.  and  are thus calculated from uncertainty statistics whose references are presented in this section. Since these references do not specify the exact nature of these statistics, we assume they correspond to standard deviations. The uncertainty of tropospheric ${\mathrm{NO}}_{\mathrm{2}}$ columns under polluted conditions is dominated by the sensitivity of satellite observations to lower tropospheric air masses, expressed by the tropospheric air-mass factor (AMF). The column relative uncertainty due to the AMF is of the order of 30 % . S-5P validation activities indicate that TROPOMI tropospheric ${\mathrm{NO}}_{\mathrm{2}}$ columns are systematically biased low by about 30 %–50 % over cities , which is most likely related to the a priori profiles used within the operational retrieval that do not reflect the ${\mathrm{NO}}_{\mathrm{2}}$ peak close to ground well. For the Middle East region, the impact of the a priori profile is less critical, as surface albedo is generally high and cloud fractions are generally low. Thus, we expect no such bias and consider a relative uncertainty of 30 % for the tropospheric column. Other uncertainties must be taken into account: the transition from ${\mathrm{NO}}_{\mathrm{2}}$ TROPOMI columns to ${\mathrm{NO}}_x$ emissions requires parameters which appear in Eqs. () and (). For the wind module, uncertainties are generally of about 1 m s${}^{-\mathrm{1}}$ for components taken at precise altitudes . Here, we assume an uncertainty of 3 m s${}^{-\mathrm{1}}$ for both zonal and meridional wind components. For [OH], the analysis of different methods conducted by showed smaller differences for low latitudes than for extratropics, but still significant. We thus take a relative uncertainty of 30 % for OH concentration. For the reaction rate $k_{\mathrm{mean}}$, the value of the corresponding relative uncertainty has been estimated by . Finally, we use the sensitivity tests performed in Sect. 4.3 to assess the uncertainty associated with the choice of the vertical level. The cumulative effects on the final emissions of the three parameters studied, in particular the OH concentration, lead to a relative uncertainty that varies from month to month between 7 % and 18 %. The propagation of these different uncertainties on the monthly estimates of ${\mathrm{NO}}_x$ emissions in Egypt leads to an expanded uncertainty between 47 % and 51 %. For lifetimes calculated with the EMG function fitting, the corresponding expanded uncertainty ranges between 18 % and 79 %.

We acknowledge the fact that our treatment of ${\mathrm{NO}}_x$ is simplified. Many minor sinks highlighted in Sect. 3.1 are not taken into account. In particular, anthropogenic VOC emissions, which remove ${\mathrm{NO}}_x$ from the atmosphere, compete with the oxidation by OH for the representation of ${\mathrm{NO}}_x$ loss. These emissions are difficult to estimate, and the corresponding sink is complex to model. Taking this reaction into account would a priori lead to a strengthening of the sink term and thus to an increase in the ${\mathrm{NO}}_x$ emissions estimates. Moreover, due to the coarse resolution of CAMS data, OH gradients might also be underestimated, especially in the southern part of the domain, leading to a local under-estimation of the sink term and the corresponding emissions. Other assumptions in the model are also simplifications. For instance, obtaining anthropogenic emissions by subtracting the average emissions over background cells assumes that the non-anthropogenic sources of ${\mathrm{NO}}_{\mathrm{2}}$ are similar inside and outside the mask, which is not true, since a large part of the mask cells correspond to croplands. For these cells, soil emissions may play a non-negligible role in the natural ${\mathrm{NO}}_{\mathrm{2}}$ budget. As a consequence, mean background emissions that are removed from ${\mathrm{NO}}_x$ emissions estimates above mask cells might be under-estimated. Finally, the reliability of the data used can be questioned. The representation of the wind is crucial to avoid creating artificial patterns in the transport term. The OH concentration, which is proportional to the intensity of the sink term, is also important. We have shown that OH concentrations are partially responsible for an important drop in ${\mathrm{NO}}_x$ emissions in the winter of 2019/20 that may be unrealistic. Because this decrease is largely due to variations in OH concentrations provided by CAMS, whose reliability has been evaluated for Riyadh, the transposability hypothesis between Riyadh and Egypt may be subject to further discussion.

5 Conclusions

In this study, we investigated the potential of a top-down model of ${\mathrm{NO}}_x$ emissions based on TROPOMI retrievals at high resolution over Egypt. The model is based on the study of a transport term and a sink term that require different parameters to be calculated. Among those parameters, the concentration of OH, involved in the calculation of the ${\mathrm{NO}}_{\mathrm{2}}$ mixed lifetime, is of fundamental importance. The comparison between ${\mathrm{NO}}_{\mathrm{2}}$ lifetimes derived from OH concentrations and ${\mathrm{NO}}_{\mathrm{2}}$ lifetimes derived from EMG function fittings of line density profiles shows that the OH concentration provided by CAMS is reasonably reliable for the country. Parameters are taken in the first 200 m of the planetary boundary layer, because it is where OH shows the best consistency. However, the vertical sensitivity linked to this parameter remains high. Results illustrate the importance of the transport term at a local scale, which is of the same order of magnitude as the sink term above large cities and industrial facilities; it ceases to be relevant only at the scale of the whole country. The top-down model is able to characterise declines in human activities due to restrictions during the COVID-19 pandemic or to Friday rest. It also estimates higher emissions during summer. These high emissions may be interpreted by a higher consumption of electricity driven by air-conditioning during hot days, but it remains unclear whether this pattern clearly reproduces changes in human activity, in particular because the emission inventories show different seasonalities. These inventories also differ in the amount of total emissions: the average value for TROPOMI-derived ${\mathrm{NO}}_x$ emissions is 7.0 % higher than CAMS-GLOB-ANT_v4.2 estimates. This discrepancy could be resolved by comparing the results of the model and inventory estimates to industrial production or electricity consumption data at the scale of countries or regions.

Here, our estimation of ${\mathrm{NO}}_x$ emissions benefited from favourable conditions. Egypt has a desertic climate, allowing us to neglect many ${\mathrm{NO}}_x$ loss mechanisms for the sink term calculation, a flat terrain in most of its territory, limiting wind field errors for the transport term calculation and a large population concentrated in a small number of cities, providing ${\mathrm{NO}}_{\mathrm{2}}$ maps with large signal-to-noise ratios above urban and industrial areas. For other regions of the world that do not have such features, the method presented here must be modified accordingly. However, we expect this method to be applicable to most countries similar to Egypt without substantial changes. For Middle East countries, this study thus demonstrates the potential of TROPOMI data for evaluating ${\mathrm{NO}}_x$ emissions. More generally, it demonstrates the importance of the contribution of independent observation systems to overcome the weaknesses of emission inventories, provided that the local chemistry is well understood and modelled. The development of similar applications for different species is likely to allow a better monitoring of global anthropogenic emissions, therefore helping companies and countries to report their emissions of air pollutants and greenhouse gases as part of their strategies and obligations to tackle air pollution issues and climate change.

Code availability

The software code for this research is not in a stand-alone user-friendly form but can be shared on request to the authors.

Data availability

The TROPOMI NO${}_{\mathrm{2}}$ product is publicly available on the TROPOMI Open hub (http://www.tropomi.eu/data-products/data-access, ). The offline mode (OFFL) data of version 1.0.0 to version 1.1.0 were used. CAMS data can be downloaded from the Copernicus Climate Data Store (https://ads.atmosphere.copernicus.eu/cdsapp#!/dataset/cams-global-atmospheric-composition-forecasts, ). The European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis be downloaded from the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels-monthly-means, ). The Global Rural-Urban Mapping Project (GRUMP) dataset is available from the NASA Socioeconomic Data and Applications Center at https://sedac.ciesin.columbia.edu/data/collection/grump-v1 . Ground data for Egypt can be taken from http://globalenergyobservatory.org/ for oil and gas power plants and from https://eogdata.mines.edu/download_global_flare.html for flaring sites. CAMS-GLOB-ANT_v4.2 data are available at https://eccad3.sedoo.fr/#CAMS-GLOB-ANT . EDGAR v5.0 Global Air Pollutant Emissions are provided by https://edgar.jrc.ec.europa.eu/overview.php?v=50_AP .

Author contributions

AR analysed the data, prepared the main software code and wrote the paper. FC provided the TROPOMI NO${}_{\mathrm{2}}$ data product and corresponding gridded maps. PC contributed to building the database of emitters for Egypt. DH contributed to the interpretation of the use of OH for NO${}_x$ loss. GB, TC, JK and JS contributed to the improvement of the method and the interpretation of the results. All authors read and agreed on the published version of the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

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

Acknowledgements

The authors would like to thank Steven J. Davis (University of California, Irvine) and Dan Tong (Tsinghua University) for their contribution to the construction of our emitter database.

Financial support

This research has been supported by the Horizon 2020 (EMME-CARE (grant no. 856612)).

Review statement

This paper was edited by Bryan N. Duncan and reviewed by two anonymous referees.