The TROPOMI surface UV algorithm

Abstract. The TROPOspheric Monitoring Instrument (TROPOMI), launched in October 2017 on board the Sentinel-5 Precursor (S5P) satellite, monitors the composition of the Earth's atmosphere at an unprecedented horizontal resolution as fine as 3.5 × 5.5 km2. This paper assesses the performances of the TROPOMI formaldehyde (HCHO) operational product compared to its predecessor, the OMI (Ozone Monitoring Instrument) HCHO QA4ECV product, at different spatial and temporal scales. The parallel development of the two algorithms favoured the consistency of the products, which facilitates the production of long-term combined time series. The main difference between the two satellite products is related to the use of different cloud algorithms, leading to a positive bias of OMI compared to TROPOMI of up to 30 % in tropical regions. We show that after switching off the explicit correction for cloud effects, the two datasets come into an excellent agreement. For medium to large HCHO vertical columns (larger than 5 × 1015 molec. cm−2) the median bias between OMI and TROPOMI HCHO columns is not larger than 10 % (< 0.4 × 1015 molec. cm−2). For lower columns, OMI observations present a remaining positive bias of about 20 % (< 0.8 × 1015 molec. cm−2) compared to TROPOMI in midlatitude regions. Here, we also use a global network of 18 MAX-DOAS (multi-axis differential optical absorption spectroscopy) instruments to validate both satellite sensors for a large range of HCHO columns. This work complements the study by Vigouroux et al. (2020), where a global FTIR (Fourier transform infrared) network is used to validate the TROPOMI HCHO operational product. Consistent with the FTIR validation study, we find that for elevated HCHO columns, TROPOMI data are systematically low (−25 % for HCHO columns larger than 8 × 1015 molec. cm−2), while no significant bias is found for medium-range column values. We further show that OMI and TROPOMI data present equivalent biases for large HCHO levels. However, TROPOMI significantly improves the precision of the HCHO observations at short temporal scales and for low HCHO columns. We show that compared to OMI, the precision of the TROPOMI HCHO columns is improved by 25 % for individual pixels and by up to a factor of 3 when considering daily averages in 20 km radius circles. The validation precision obtained with daily TROPOMI observations is comparable to the one obtained with monthly OMI observations. To illustrate the improved performances of TROPOMI in capturing weak HCHO signals, we present clear detection of HCHO column enhancements related to shipping emissions in the Indian Ocean. This is achieved by averaging data over a much shorter period (3 months) than required with previous sensors (5 years) and opens new perspectives to study shipping emissions of VOCs (volatile organic compounds) and related atmospheric chemical interactions.

TROPOspheric Monitoring Instrument observations of total column water vapour: Algorithm and validation

<p>We present the total column water vapor (TCWV) retrieval for the TROPOspheric Monitoring Instrument (TROPOMI) observations in the blue band. The retrieval was first developed to retrieve TCWV from Global Ozone Monitoring Experience 2 (GOME-2). We have modified the settings of the retrieval to adapt it for TROPOMI observations. The TROPOMI TCWV retrieval algorithm consists of two major steps. The first step is the retrieval of water vapor slant columns by applying the differential optical absorption spectroscopy (DOAS) technique to TROPOMI observations in the blue band. The retrieved water vapor slant columns are then converted to vertical columns using air mass factors (AMFs). An iterative optimization has been developed to dynamically find the optimal a priori water vapor profile for AMF calculation. The dynamic search algorithm makes use of the fact that the vertical distribution of water vapor is strongly correlated to the total column amount. This makes the algorithm better suited for climate studies compared to typical satellite retrievals with static a priori or vertical profile information from the chemistry transport model (CTM). Details of the TCWV retrieval are presented. The TCWV retrieval algorithm is applied to TROPOMI observations. The results are validated by comparing to Ozone Monitoring Instrument (OMI), GOME-2 and Special Sensor Microwave Imager Sounder (SSMIS) satellite observations. TCWV derived from TROPOMI observations agree well with the other data sets with Pearson correlation coefficient (R) ranging from 0.94 to 0.99. The correlation is slight better during winter time of the northern hemisphere. Small discrepancies are found among TROPOMI, OMI, GOME-2 and SSMIS observations. The discrepancies are mainly due to differences in measurement time and cloud filtering. More detailed validation against ground based sun-photometer observations are presented separately in this session*.</p><p> </p><p>*see the respective abstract by Katerina Garane.</p>

Sulfur dioxide (SO2) is released into the Earth’s atmosphere through natural and anthropogenic processes and the latter category amounts to the majority of the global SO2 emissions. Satellite SO2 observations have been used to monitor SO2 emissions at different regional and global scales and to detect large-point sources of SO2 emissions of diverse origins. In this study, we conducted atmospheric inversions at the global scale to estimate daily SO2 emissions at 1.26o×2.5o (latitude×longitude) spatial resolution over polluted regions for two years 2020 and 2021 using satellite SO2 total vertical column densities (TVCDs) obtained by the Sentinel 5p TROPOspheric Monitoring Instrument (TROPOMI) and AURA Ozone Monitoring Instrument (OMI). We used the global chemistry coupled transport model LMDz-INCA with 1.26o×2.5o (latitude×longitude) horizontal resolution and 79 hybrid σ-p vertical levels extending to the stratosphere to simulate the model SO2 TVCDs. The model uses a priori monthly global anthropogenic emission inventories from the open-source Community Emissions Data System (CDES). As the TROPOMI operational offline L2 SO2 data product has high noise levels, we used the TROPOMI COBRA SO2 data product in this study which has comparatively smaller noise. First, we evaluated the SO2 TVCDs from the LMDz-INCA model simulations for a reference year 2019 with the observed SO2 TVCDs from TROPOMI and OMI. The daily average (10-day running average) of the model simulated SO2 TVCDs over the major polluted regions like India, China, and the Middle East, and over less polluted regions like South Africa, and South America mostly follow the trend of the observed SO2 TVCDs from both TROPOMI and OMI. The model overestimates SO2 TVCDs over India and the Middle East and underestimates them over China, South Africa, and South America. For Europe and North America, the noise levels in the daily averaged TVCDs from both TROPOMI and OMI are too high for a meaningful comparison. In order to estimate anthropogenic SO2 emissions, we used a recently developed inversion approach (Zheng et al., 2020), which was previously used to estimate anthropogenic NOx emissions over China. We performed the model simulation for 2019 with 40% reduced anthropogenic SO2 emission to calculate the gridded local sensitivities of the TVCDs to the change in the anthropogenic SO2 emissions. The inversion approach combines these gridded local sensitivities and the relative change of the observed satellites and the modelled TVCDs to derive the relative change of anthropogenic SO2 emissions from the reference year 2019 to the inversion years 2020 and 2021. The estimated total SO2 emissions from TROPOMI and OMI observations for 2020 and 2021 are mostly higher compared to the reference year total emissions over the world and over the selected regions. The total SO2 emissions from TROPOMI and OMI observations at the common model grids for both inversion years are consistent with each other.

This study gives a systematic comparison of the Tropospheric Monitoring Instrument (TROPOMI) version 1.2 and Ozone Monitoring Instrument (OMI) QA4ECV tropospheric NO2 column through global chemical data assimilation (DA) integration for the period April−May 2018. DA performance is controlled by measurement sensitivities, retrieval errors, and coverage. The smaller mean relative observation errors by 16 % in TROPOMI than OMI over 60° N−60° S during April−May 2018 led to larger reductions in the global root mean square error (RMSE) against the assimilated NO2 measurements in TROPOMI DA (by 54 %) than in OMI DA (by 38 %). Agreements against the independent surface, aircraft-campaign, and ozonesonde observation data were also improved by TROPOMI DA compared to the control model simulation (by 12−84 % for NO2 and by 7−40 % for ozone), which were more obvious than those by OMI DA for many cases (by 2−70 % for NO2 and by 1−22 % for ozone). The estimated global total NOx emissions were 15 % lower in TROPOMI DA, with 2−23 % smaller regional total emissions, in line with the observed negative bias of the TROPOMI version 1.2 product compared to the OMI QA4ECV product. TROPOMI DA can provide city scale emission estimates, which were within 10 % differences with other high-resolution analyses for several limited areas, while providing a globally consistent analysis. These results demonstrate that TROPOMI DA improves global analyses of NO2 and ozone, which would also benefit studies on detailed spatial and temporal variations in ozone and nitrate aerosols and the evaluation of bottom-up NOx emission inventories.

Abstract. Nitrogen dioxide (NO2) is one of the main data products measured by the Tropospheric Monitoring Instrument (TROPOMI) on the Sentinel-5 Precursor (S5P) satellite, which combines a high signal-to-noise ratio with daily global coverage and high spatial resolution. TROPOMI provides a valuable source of information to monitor emissions from local sources such as power plants, industry, cities, traffic and ships, and variability of these sources in time. Validation exercises of NO2 v1.2–v1.3 data, however, have revealed that TROPOMI's tropospheric vertical column densities (VCDs) are too low by up to 50 % over highly polluted areas. These findings are mainly attributed to biases in the cloud pressure retrieval, the surface albedo climatology and the low resolution of the a priori profiles derived from global simulations of the TM5-MP chemistry model. This study describes improvements in the TROPOMI NO2 retrieval leading to version v2.2, operational since 1 July 2021. Compared to v1.x, the main changes are the following. (1) The NO2-v2.2 data are based on version-2 level-1b (ir)radiance spectra with improved calibration, which results in a small and fairly homogeneous increase in the NO2 slant columns of 3 % to 4 %, most of which ends up as a small increase in the stratospheric columns. (2) The cloud pressures are derived with a new version of the FRESCO cloud retrieval already introduced in NO2-v1.4, which led to a lowering of the cloud pressure, resulting in larger tropospheric NO2 columns over polluted scenes with a small but non-zero cloud coverage. (3) For cloud-free scenes a surface albedo correction is introduced based on the observed reflectance, which also leads to a general increase in the tropospheric NO2 columns over polluted scenes of order 15 %. (4) An outlier removal was implemented in the spectral fit, which increases the number of good-quality retrievals over the South Atlantic Anomaly region and over bright clouds where saturation may occur. (5) Snow/ice information is now obtained from ECMWF weather data, increasing the number of valid retrievals at high latitudes. On average the NO2-v2.2 data have tropospheric VCDs that are between 10 % and 40 % larger than the v1.x data, depending on the level of pollution and season; the largest impact is found at mid and high latitudes in wintertime. This has brought these tropospheric NO2 closer to Ozone Monitoring Instrument (OMI) observations. Ground-based validation shows on average an improvement of the negative bias of the stratospheric (from −6 % to −3 %), tropospheric (from −32 % to −23 %) and total (from −12 % to −5 %) columns. For individual measurement stations, however, the picture is more complex, in particular for the tropospheric and total columns.

Nitrogen dioxide (NO2) is an important trace-gas pollutant and climate agent whose presence also leads to spectral interference in ocean color retrievals. NO2 column densities have been retrieved with satellite UV-Vis spectrometers such as the Ozone Monitoring Instrument (OMI) and Tropospheric Monitoring Instrument (TROPOMI) that typically have spectral resolutions of the order of 0.5 nm or better and spatial footprints as small as 3.5 km × 5 km. These NO2 observations are used to estimate emissions, monitor pollution trends, and study effects on human health. Here, we investigate whether it is possible to retrieve NO2 amounts with lower spectral resolution hyper-spectral imagers such as the Ocean Color Instrument (OCI) that will fly on the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) satellite set for launch in early 2024. OCI will have a spectral resolution of 5 nm and a spatial resolution of ∼1 km with global coverage in 1–2 days. At this spectral resolution, small scale spectral structure from NO2 absorption is still present. We use real spectra from the OMI to simulate OCI spectra that are in turn used to estimate NO2 slant column densities (SCDs) with an artificial neural network trained on target OMI retrievals. While we obtain good results with no noise added to the OCI simulated spectra, we find that the expected instrumental noise substantially degrades the OCI NO2 retrievals. Nevertheless, the NO2 information from OCI may be of value for ocean color retrievals, as our simulations suggest that it will be of similar or slightly better quality as compared with TROPOMI NO2 data at TROPOMI spatial resolution on a daily basis and will be available simultaneously. OCI retrievals can also be temporally averaged over time-scales of the order months to reduce noise and provide higher spatial resolution maps that may be useful for downscaling information provided by lower spatial resolution instruments such as OMI and TROPOMI, for high resolution emissions estimates, and other applications. In addition, we explore the possibility of using an extended fitting window for NO2 retrievals as compared with traditional approaches. We demonstrate that the use of an extended spectral fitting window can reduce random errors in a current state-of-the-art OMI NO2 SCD product. Machine learning approaches with extended fitting windows, once trained, can also substantially speed up NO2 spectral fitting algorithms as applied to OMI, TROPOMI, and similar instruments that are flying or will soon fly in geostationary orbit.

Nitrogen dioxide (NO2) is an important trace-gas pollutant and climate agent whose presence also leads to spectral interference in ocean color retrievals. NO2 column densities have been retrieved with satellite UV-Vis spectrometers such as the Ozone Monitoring Instrument (OMI) and Tropospheric Monitoring Instrument (TROPOMI) that typically have spectral resolutions of the order of 0.5 nm or better and spatial footprints as small as 3.5 km × 5 km. These NO2 observations are used to estimate emissions, monitor pollution trends, and study effects on human health. Here, we investigate whether it is possible to retrieve NO2 amounts with lower spectral resolution hyper-spectral imagers such as the Ocean Color Instrument (OCI) that will fly on the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) satellite set for launch in early 2024. OCI will have a spectral resolution of 5 nm and a spatial resolution of ∼1 km with global coverage in 1–2 days. At this spectral resolution, small scale spectral structure from NO2 absorption is still present. We use real spectra from the OMI to simulate OCI spectra that are in turn used to estimate NO2 slant column densities (SCDs) with an artificial neural network trained on target OMI retrievals. While we obtain good results with no noise added to the OCI simulated spectra, we find that the expected instrumental noise substantially degrades the OCI NO2 retrievals. Nevertheless, the NO2 information from OCI may be of value for ocean color retrievals, as our simulations suggest that it will be of similar or slightly better quality as compared with TROPOMI NO2 data at TROPOMI spatial resolution on a daily basis and will be available simultaneously. OCI retrievals can also be temporally averaged over time-scales of the order months to reduce noise and provide higher spatial resolution maps that may be useful for downscaling information provided by lower spatial resolution instruments such as OMI and TROPOMI, for high resolution emissions estimates, and other applications. In addition, we explore the possibility of using an extended fitting window for NO2 retrievals as compared with traditional approaches. We demonstrate that the use of an extended spectral fitting window can reduce random errors in a current state-of-the-art OMI NO2 SCD product. Machine learning approaches with extended fitting windows, once trained, can also substantially speed up NO2 spectral fitting algorithms as applied to OMI, TROPOMI, and similar instruments that are flying or will soon fly in geostationary orbit.

Spaceborne NO2 column observations from two high‐resolution instruments, Tropospheric Monitoring Instrument (TROPOMI) on board Sentinel‐5 Precursor and Ozone Monitoring Instrument (OMI) on Aura, reveal unprecedented NO2 decreases over China, South Korea, western Europe, and the United States as a result of public health measures enforced to contain the coronavirus disease outbreak (Covid‐19) in January–April 2020. The average NO2 column drop over all Chinese cities amounts to −40% relative to the same period in 2019 and reaches up to a factor of ~2 at heavily hit cities, for example, Wuhan, Jinan, while the decreases in western Europe and the United States are also significant (−20% to −38%). In contrast with this, although Iran is also strongly affected by the disease, the observations do not show evidence of lower emissions, reflecting more limited health measures.

Solar-Induced chlorophyll Fluorescence (SIF) retrieved from spaceborne hyperspectral measurements can provide a synoptic perspective on the functioning of photosynthetic active vegetation and algae. The TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 precursor mission provides data for the retrieval of SIF with an unprecedented resolution in the spatial and temporal dimensions. In this study, we propose an approach for the retrieval of the red SIF peak value (at around 685 nm) over both land and ocean from TROPOMI measurements. This approach exploits the solar Fraunhofer lines located in the spectral window of 669–686 nm to disentangle SIF from the solar radiation reflected by the surface–atmosphere system. First, through simulated TROPOMI-like datasets, we determine the proper parameter settings and demonstrate the feasibility of the approach for the red SIF retrieval from TROPOMI data. Then, the approach is applied to real TROPOMI measurements. The red SIF retrieved with the proposed approach for two months’ TROPOMI data is compared with another TROPOMI red SIF dataset and the marine normalized fluorescence line height (nFLH) product of the Aqua Moderate Resolution Imaging Spectroradiometer (MODIS). The comparisons display similar spatial distributions over both land and ocean for the weekly global composites. Especially, the comparison of the two red SIF datasets demonstrates good agreement overall, indicating consistency of the two retrieval approaches. The retrieval uncertainty for the weekly global composite is about 13% and 10% of the peak SIF value over land/ocean, respectively, which can be considered as satisfactory error thresholds for global SIF observations.

<p>Along with monitoring of seismic activity and ground deformation, the measurement of volcanic gas emissions and composition plays a key role in the surveillance of active volcanoes and the mitigation of volcanic hazards. Volcanic gas emissions also potentially impact the environment, human health and climate, providing further motivation for study. Currently, volcano observatories typically employ ground-based or airborne techniques to monitor volcanic gas emissions, mainly sulfur dioxide (SO<sub>2</sub>) fluxes and its ratios over other species (e.g., CO<sub>2</sub>, H<sub>2</sub>S). However, in recent years there have been significant breakthroughs in satellite observations of passive volcanic SO<sub>2</sub> emissions, including high-resolution ultraviolet (UV) measurements from the Tropospheric Monitoring Instrument (TROPOMI) on the Sentinel-5 Precursor (S5P) satellite, and the development of long-term records of volcanic SO<sub>2</sub> degassing from the Ozone Monitoring Instrument (OMI) on NASA’s Aura satellite. Satellite measurements offer some advantages over traditional gas monitoring techniques, e.g., synoptic coverage of large regions, relative immunity to variations in wind direction, and ability to map the spatial extent and dispersion of volcanic SO<sub>2</sub> plumes with applications for health hazard mitigation. Although these satellite datasets are potentially valuable for active volcano monitoring and as a supplement to other gas monitoring techniques, significant barriers remain to their use at many volcano observatories, particularly in low-income countries. Notably, the increasing volume of satellite datasets (NASA’s database is bigger than 3 petabytes) and the demands of data processing represent challenges to their operational use at observatories with limited internet connectivity or computational capacity. Here, we present an ongoing effort to develop open-source Python software to access and process SO<sub>2</sub> data directly through NASA’s Earthdata portal Application Processing Interface (API), in order to streamline the satellite SO<sub>2</sub> data processing workflow for a volcano observatory. By allowing server-side satellite data subsetting around the volcano of interest, this API greatly reduces the processing burden and only requires an internet connection to the NASA server hosting the required datasets (including S5P/TROPOMI, Aura/OMI and many others). We present some examples of software output and potential applications. Our current goal is to deploy and test the software for operational use in a volcano observatory.  </p>

A new method is presented to estimate urban hydroxyl radical (OH) concentrations using the downwind decay of the Tropospheric Monitoring Instrument (TROPOMI) derived nitrogen dioxide (NO2) / carbon monoxide (CO) ratio combined with Weather Research Forecast (WRF) simulations. Seasonal OH concentrations, nitrogen oxides (NOx) and CO emissions for summer (June to October, 2018) and winter (November, 2018 to March, 2019) are derived for Riyadh. WRF is able to spatially simulate NO2 and CO urban plumes over Riyadh as observed by TROPOMI. However, WRF-simulated NO2 plumes close to center of the city are overestimated by 25 % in summer and 40 to 50 % in winter compared to TROPOMI observations. WRF simulated CO plumes differ by 10 % with TROPOMI in both seasons. The differences between model and TROPOMI are used to optimize the OH concentration, NOx and CO emissions iteratively using a least squares method. For summer, both the NO2 / CO ratio optimization and the XNO2 optimization imply that the OH prior from the Copernicus Atmospheric Monitoring Service (CAMS) has to be increased by 32.03±4.0 % . The OH estimations from the NO2 / CO ratio and the XNO2 optimization differ by 10 %. Summer Emission Database for Global Atmospheric Research v4.3.2 (EDGAR) NOx and CO emissions over Riyadh need to be increased by 42.1±8.7 % and 100.8±9.5 %. For winter, the optimization method increases OH by ~52.0±5.3 %, while reducing NOx emission by 15.45± 3.4 % and doubling the CO emission. TROPOMI derived OH concentrations and pre-existing Exponentially Modified Gaussian function fit (EMG) method differ by 18 % in summer and 7.5 % in winter, confirming that urban OH concentrations can be reliably estimated using the TROPOMI-observed NO2 / CO ratio.

<strong class="journal-contentHeaderColor">Abstract.</strong> A new method is presented to estimate urban hydroxyl radical (OH) concentrations using the downwind decay of the Tropospheric Monitoring Instrument (TROPOMI) derived nitrogen dioxide (NO₂) / carbon monoxide (CO) ratio combined with Weather Research Forecast (WRF) simulations. Seasonal OH concentrations, nitrogen oxides (NOx) and CO emissions for summer (June to October, 2018) and winter (November, 2018 to March, 2019) are derived for Riyadh. WRF is able to spatially simulate NO₂ and CO urban plumes over Riyadh as observed by TROPOMI. However, WRF-simulated NO₂ plumes close to center of the city are overestimated by 25 % in summer and 40 to 50 % in winter compared to TROPOMI observations. WRF simulated CO plumes differ by 10 % with TROPOMI in both seasons. The differences between model and TROPOMI are used to optimize the OH concentration, NOx and CO emissions iteratively using a least squares method. For summer, both the NO₂/ CO ratio optimization and the XNO₂ optimization imply that the OH prior from the Copernicus Atmospheric Monitoring Service (CAMS) has to be increased by 32.03&plusmn;4.0 % . The OH estimations from the NO₂/ CO ratio and the XNO₂ optimization differ by 10 %. Summer Emission Database for Global Atmospheric Research v4.3.2 (EDGAR) NOx and CO emissions over Riyadh need to be increased by 42.1&plusmn;8.7 % and 100.8&plusmn;9.5 %. For winter, the optimization method increases OH by ~52.0&plusmn;5.3 %, while reducing NOx emission by 15.45&plusmn; 3.4 % and doubling the CO emission. TROPOMI derived OH concentrations and pre-existing Exponentially Modified Gaussian function fit (EMG) method differ by 18 % in summer and 7.5 % in winter, confirming that urban OH concentrations can be reliably estimated using the TROPOMI-observed NO₂/ CO ratio.

Nitrogen oxides (NOx) are air pollutants critical to ozone and fine particle production in the troposphere. Here, we present fuel‐based emission inventories updated to 2018, including for mobile source engines using the Fuel‐based Inventory of Vehicle Emissions (FIVEs) and oil and gas production using the Fuel‐based Oil and Gas (FOG) inventory. The updated FIVE emissions are now consistent with the NEI17 estimates differing within 2% across the contiguous US (CONUS). Tropospheric NO2 columns modeled by the Weather Research and Forecasting with Chemistry model (WRF‐Chem) are compared with those observed by TROPOspheric Monitoring Instrument (TROPOMI) and Ozone Monitoring Instrument (OMI) during the summer of 2018. Modeled NO2 columns show strong temporal and spatial correlations with TROPOMI (OMI), identified with biases of −3% (−21%) over CONUS, and +8% (−6%) over point sources plus urban regions. Taking account of the negative bias (∼20%) in early version of TROPOMI over polluted regions, WRF‐Chem shows good performance with updated FIVE and FOG emissions. Our model tends to under‐predict the tropospheric NO2 columns over background and rural regions (bias of −21% to −3%). Through model sensitivity analyses, we demonstrate the important roles of emissions from soils (11.7% average over CONUS), oil and gas production (4.1%), wildfires (10.6%), and lightning (2.3%) with greater contributions at regional scales. This work provides a roadmap for satellite‐based evaluations for emission updates from various sources.

The variation of tropospheric nitrogen dioxide (NO2) vertical column densities (VCDs) indirectly reflects the difference in pollution emissions from industrial production and transportation. Accurately analyzing its pollution sources and driving factors plays an important role in energy conservation, emission reduction, and air pollution reduction. NO2 concentration products of Sentinel-5P (Sentinel-5 Precursor) TROPOMI (TROPOspheric Monitoring Instrument) from 2019 to 2021 and Aura OMI (Ozone Monitoring Instrument) from 2009 to 2021, combined with China’s main energy consumption, the growth value of the industry, Gross Domestic Product (GDP), and other data were used to analyze the influencing factors of NO2 variations. Firstly, NO2 tropospheric vertical column densities (NO2 TVCDs) of China increased by 14.72% and 3.26% in 2021 and 2020 compared with the 2019. The secondary and tertiary industry and the national energy consumption increased synchronously, which was highly related to the increase in NO2 TVCDs. Secondly, the impact of COVID-19 (coronavirus disease 2019) on China’s industrial production and residents was mainly concentrated in the first quarter of 2020, which leading to a decline in the annual average NO2 concentration in densely populated areas in 2020 compared to the same period in 2019. The industrial production scale and production capacity has gradually recovered since April 2020, and the NO2 concentration has gradually reached or exceeded the level of the same period of 2019. Finally, atmospheric pollution prevention and control measures played a positive role in the decline of NO2 of China.