Tropospheric NO 2 Pollution Monitoring with the GF-5 Satellite Environmental Trace Gases Monitoring Instrument over the North China Plain during Winter 2018–2019

: The Environmental Trace Gases Monitoring Instrument (EMI) is a high-spectral-resolution payload onboard the latest pathfinder mission GaoFen-5, designed specifically for the monitoring of global atmospheric trace gas compositions and trends. This study describes a comparative analysis of the tropospheric nitrogen dioxide (NO 2 ) columns over the North China Plain (NCP) from November 2018 to April 2019 based on EMI products. Validation of satellite products based on a cross-correlation analysis with data from four ground-based multi-axis differential optical absorption spectroscopy sites provided good correlation coefficients (r) ranging from 0.78 to 0.88. The distribution and monthly averaged tropospheric NO 2 columns revealed high pollution exposure levels during winter (November–January) and a decrease from February onward in the NCP. Moreover, a typical pollution event was analyzed in detail in combination with wind field statistics. The results indicated that variations of NO 2 concentrations in Beijing and Tianjin were highly correlated with the wind direction from 22.5–45.0 degrees west of south, especially during times of high NO 2 amounts. These findings highlight that the EMI payload on the GaoFen-5 (GF-5) satellite is useful for remote sensing of regional and global NO 2 detection.


Introduction
Nitrogen dioxide (NO 2 ) is a key trace gas in the atmosphere that has important impacts on atmospheric chemistry. It is not only an important precursor of ozone, photochemical smog, and acid rain [1][2][3], but also can exert negative effects on human health in its original chemical state [4][5][6]. Atmospheric NO 2 comes from both natural and anthropogenic sources. Natural emission sources include lightning, which promotes the transformation of nitrogen into nitrogen oxides, and forest fires, which lead to the combustion of biomass and the production of NO 2 [7]. Human activities, such as the burning of fossil fuels and industrial activities, are the main causes of increases in NO 2 concentrations in the troposphere [8]. Therefore, in areas such as eastern China, eastern US, eastern India, and the Persian Gulf, where the population is relatively large, transportation networks are well developed, and industrialization is relatively concentrated, high concentrations of NO 2 are typically found [9,10]. In the North China Plain (NCP), NO 2 concentrations can reach high levels during the winter in response to the burning of fossil fuels for heating [11,12].
Because trace gases such as NO 2 in the troposphere are closely related to the atmospheric environment and human activities, it is important to observe continuously and accurately the distribution and changes of these pollutants. Air pollution data from over 1500 monitoring stations show that significant pollution levels are widespread across northern and central China, and 92% of the population of China experienced unhealthy air conditions for a relatively long period from 5 April 2014, to 5 August 2014 [13]. Recently, there have been encouraging signs that NO 2 concentrations in megacities such as Beijing, Shanghai, and Guangzhou are showing downward trends, although NO 2 concentrations in the NCP have continued to increase [14,15].
Satellite remote sensing plays a crucial role in monitoring the atmospheric composition and its changes, by providing a matchless global perspective with consistency over long periods. Recently, instruments located on satellites in orbit, such as Global Ozone Monitoring Experiment (GOME2), Ozone Mapping and Profiler Suite (OMPS), TROPOspheric Monitoring Instrument (TROPOMI), and Ozone Monitoring Instrument (OMI), have been applied widely for NO 2 pollution detection. Nitrogen dioxide data derived from GOME2 and OMI show that the discharge of anthropogenic nitrogen dioxide in China is higher in winter and summer than in other seasons, and the emission ratios of weekdays to weekends have increased with the increase in nitrogen dioxide emissions [16]. During the 2020 coronavirus disease 2019 (COVID-19) pandemic, TROPOMI and OMI products were applied to observe air pollution in megacities worldwide. Due to the effects of the coronavirus lockdown adopted by the Chinese authorities, the NO 2 concentrations in different cities of China were significantly reduced [17,18]. However, the accuracy of pollution levels observed by spaceborne instruments can be strongly limited by the spatial resolution of the satellite payload in the target region, and the trace gas retrievals may not accurately reflect the pollution when the regional sampling range is large. Furthermore, the observation of diurnal cycles is limited by the small number of overpasses and temporal sampling for sun-synchronous orbit satellites. Ground-based differential optical absorption spectroscopy (DOAS) observations have the advantage of high temporal resolutions and the measurement of scattered sunlight from different, mostly slant elevation angles, which can be used to accurately retrieve the vertical profiles of aerosols and trace gases [19,20]. In addition, these data also provide direct verification and synergy analysis information for satellite observation results [21][22][23]. Therefore, a combination of spaceborne and groundbased observations can be very useful for deriving a complete picture of the spatiotemporal variability of particular pollutants.
As an important monitoring method, NO 2 monitoring using instruments on different platforms has attracted a great deal of scientific attention culminating in numerous studies, most of which have used single-type platforms for monitoring or comparative verification purposes [24][25][26]. Only a few studies have focused on the analysis of pollution trends through joint monitoring in the NCP [27,28]. Notably, numerical models and observational data can be combined to analyze meteorological conditions, transportation, and potential sources of heavy pollution processes [29][30][31][32]. As the most densely populated and industrially developed region, NO 2 concentrations have remained high in the NCP [33]. Dominated by the East Asian winter monsoon and the influence of the terrain, the pollution diffusion and transport in the NCP is more complicated [34]. Therefore, it is of great significance to identify the sources and sinks during heavy pollution events in the NCP.
We used China's first hyperspectral atmospheric trace gas detection payload, the Environmental Trace Gases Monitoring Instrument (EMI), to analyze the tropospheric NO 2 distribution and its changes in the NCP region from November 2018 to April 2019. We combined these data with ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) data for semi-annual comparative observations. In addition, with the aid of meteorological data and wind fields, the spatial variation of NO 2 , pollution trends and sources, and transportation contributions were analyzed by using statistics.

EMI Tropospheric NO 2 Data
The first full-spectrum hyperspectral sun-synchronous orbit satellite of China, GaoFen-5 (GF-5), was launched in May 2018, and it was placed into orbit to achieve comprehensive observations of the atmosphere and land. The EMI is an important sensor on the GF-5 satellite, which can scan the world in one day with high near-nadir resolution (nadir resolution Atmosphere 2021, 12, 398 3 of 17 is~13 km × 12 km) and a high spectral resolution (0.3 nm) [35,36]. The instrument is mainly used for the quantitative monitoring of the global atmospheric trace gas composition and trends, environmental quality supervision, climate change detection, and implementation of hyperspectral remote sensing applications such as those involving polluting gases and regional environmental air quality monitoring.
The EMI is a passive nadir-viewing imager with four individual spectral channels. It collects hyperspectral remote sensing data with a spectral range from ultraviolet to visible wavelengths. It uses the DOAS inversion algorithm to process the spectral data within the wavelength range 411-450 nm, as the efficiency of this algorithm has been confirmed for its global NO 2 retrievals [37,38].
Because near-real-time (NRT) data products may ignore the effects of wavelength shifts and instrument degradation, which may lead to a decrease in the accuracy of the results, this study used EMI offline (OFFL) data products. During the process of obtaining the OFFL products, spectral correction was performed on the collected spectra, and the optimal inversion band was subsequently determined by considering the fitting residual and the absorption cross-section (NO 2 , O 3 , H 2 O, etc.). The fitting wavelengths and absorption crosssections used in the DOAS fitting algorithm are described in Table 1. Finally, we selected the band in the range of 411-450 nm for the NO 2 retrievals. Combined with the solar irradiance data collected by the EMI in outer space, DOAS technology was used to obtain the total slant column density of NO 2 . All spectral data were provided by the National Environmental Protection Agency. Figure 1 shows an example of NO 2 inversion using DOAS technology with spectra collected by the EMI. The inversion error was significantly less than 10%. The tropospheric NO 2 vertical column density (VCD) has a direct relevance to the regional pollution transport and distribution patterns. Studies have shown that the NO 2 tropospheric columns in remote ocean areas are low [43]; therefore, the NO 2 columns in the stratosphere and troposphere can be separated by using the reference sector method (RSM) [44]; the tropospheric NO 2 slant column density (SCD) is then converted into the tropospheric NO 2 VCD by introducing the Air Mass Factor (AMF) using the radiance transfer model. This study focused on 33 • N-42 • N and 109 • E-125 • E, which covers the area of NCP. The tropospheric NO2 vertical column density (VCD) has a direct relevance to the regional pollution transport and distribution patterns. Studies have shown that the NO2 tropospheric columns in remote ocean areas are low [43]; therefore, the NO2 columns in the stratosphere and troposphere can be separated by using the reference sector method (RSM) [44]; the tropospheric NO2 slant column density (SCD) is then converted into the tropospheric NO2 VCD by introducing the Air Mass Factor (AMF) using the radiance transfer model. This study focused on 33°N-42°N and 109°E-125°E, which covers the area of NCP.

Ground-Based MAX-DOAS Data
The same MAX-DOAS measurement system deployed in the NCP has also been used to retrieve the tropospheric columns and vertical distributions of trace gases such as NO2 and formaldehyde, and to quantify the pollution at a certain point [45]; this system was manufactured by the Anhui Institute of Optics and Fine Mechanics, China. The system is mainly composed of a spectrum acquisition module and a data processing module ( Figure  2). It uses a stepping motor to control the steering of the lens collecting the spectrum at elevation angles of 5, 15, 30, 60, and 90, and measures the NO2 amounts on different optical paths. A full collection cycle requires approximately 5-15 min owing to the automatic integral time.

Ground-Based MAX-DOAS Data
The same MAX-DOAS measurement system deployed in the NCP has also been used to retrieve the tropospheric columns and vertical distributions of trace gases such as NO 2 and formaldehyde, and to quantify the pollution at a certain point [45]; this system was manufactured by the Anhui Institute of Optics and Fine Mechanics, China. The system is mainly composed of a spectrum acquisition module and a data processing module ( Figure 2). It uses a stepping motor to control the steering of the lens collecting the spectrum at elevation angles of 5, 15, 30, 60, and 90, and measures the NO 2 amounts on different optical paths. A full collection cycle requires approximately 5-15 min owing to the automatic integral time. We used the zenith measurement spectrum recorded at noon locally as the reference spectrum to obtain the NO2 differential slant column density (DSCD), which significantly minimized the interference from the stratospheric atmosphere on the tropospheric measurements. The best fitting wavelength intervals and the absorption cross-sections considered are described in Table 1. Then, by using the DSCD of different angles to obtain the local NO2 profile, we combined the differential air mass factor (dAMF) hold by the radiation transfer model to obtain the tropospheric NO2 VCDs. We used the zenith measurement spectrum recorded at noon locally as the reference spectrum to obtain the NO 2 differential slant column density (DSCD), which significantly Atmosphere 2021, 12, 398 5 of 17 minimized the interference from the stratospheric atmosphere on the tropospheric measurements. The best fitting wavelength intervals and the absorption cross-sections considered are described in Table 1. Then, by using the DSCD of different angles to obtain the local NO 2 profile, we combined the differential air mass factor (dAMF) hold by the radiation transfer model to obtain the tropospheric NO 2 VCDs. Table 2 shows the details of the four ground-based MAX-DOAS sites in the BTH region. Due to the need for daily maintenance (calibration) of the spectrometer, data on certain days were missing. In addition, the observation azimuth of the telescope was taken into account to prevent direct sunlight and building blockages. As the DOAS observation system uses scattered sunlight as the light source for spectral detection, cloudy or foggy weather significantly decreases the accuracy of the NO 2 inversion [46]. Therefore, all data mentioned in this study from the EMI and ground-based sensors were selected only on days when the satellite passed the ground station and the relative cloud cover was less than 0.5. The relative cloud cover was characterized by the cloud reflectivity obtained from EMI radiance products [47].

Auxiliary Data
When analyzing the pollution distribution, it is necessary to consider the influence of atmospheric transport. In this study, wind speed and wind direction data were extracted from the Copernicus Atmosphere Monitoring Service (CAMS) reanalysis data [48]. The smog levels were provided by the State Environmental Protection Administration (EPA) to screen NO 2 data from the EMI and ground-based MAX-DOAS observations in the NCP (http://www.pm25.in/, accessed on 19 March 2021).

Methodology
A direct comparison of the trend in NO 2 observed by the EMI and MAX-DOAS will cause the results with large differences to not be directly reflected in the correlation. Therefore, this study analyzed the results of the mean difference (Equation (1)) and mean relative difference (Equation (2)) of NO 2 VCD: where MD is the mean absolute difference, and MRD is the mean relative difference. A positive (negative) MD is thus an indication of EMI overestimation (underestimation), and a small (large) MRD value reflects high (low) agreement between the EMI and the groundbased observation results. Additionally, the Pearson correlation coefficient and fitting slope were used to validate the results and to explain the correspondence between the EMI NO 2 vertical column densities and the corresponding averages from the ground-based MAX-DOAS measurements.
In the analysis of the impact of the wind field on pollution change in the NCP region using EMI data products, a combination of cluster analysis and two-way analysis of Atmosphere 2021, 12, 398 6 of 17 variance (ANOVA) was performed [49]. Because different daily wind speeds and directions might easily cause large degrees of freedom as independent variables, the wind speed, as a meteorological factor, was rounded according to the actual wind speed, and the wind direction data were divided into 16 clusters based on the direction of 22.5 • . Finally, taking the NO 2 column as the dependent variable, a two-way analysis of variance was used to quantitatively analyze the NO 2 changes.

Spatial Distribution of Tropospheric NO 2 in the NCP
The average distribution of the NO 2 VCD in the troposphere of mainland China from November 2018 to April 2019 was measured via spectrum inversion using EMI data ( Figure 3). The NO 2 distribution in China was uneven. For example, some economically developed regions such as Urumqi (Figure 3a), Beijing-Tianjin-Hebei (BTH) (Figure 3b ground-based observation results. Additionally, the Pearson correlation coefficient and fitting slope were used to validate the results and to explain the correspondence between the EMI NO2 vertical column densities and the corresponding averages from the groundbased MAX-DOAS measurements. In the analysis of the impact of the wind field on pollution change in the NCP region using EMI data products, a combination of cluster analysis and two-way analysis of variance (ANOVA) was performed [49]. Because different daily wind speeds and directions might easily cause large degrees of freedom as independent variables, the wind speed, as a meteorological factor, was rounded according to the actual wind speed, and the wind direction data were divided into 16 clusters based on the direction of 22.5°. Finally, taking the NO2 column as the dependent variable, a two-way analysis of variance was used to quantitatively analyze the NO2 changes.

Spatial Distribution of Tropospheric NO2 in the NCP
The average distribution of the NO2 VCD in the troposphere of mainland China from November 2018 to April 2019 was measured via spectrum inversion using EMI data (Figure 3). The NO2 distribution in China was uneven. For example, some economically developed regions such as Urumqi (Figure 3a), Beijing-Tianjin-Hebei (BTH) (Figure 3b   A detailed assessment of NO 2 pollution also is useful for local-scale analyses to better evaluate the role of pollution sources and pollution changes. For these purposes, NO 2 pollution in the NCP (32 • N-42 • N, 110 • E-120 • E), where the capital city of China is located, was analyzed in detail ( Figure 4).
The NO 2 distribution maps from November 2018 to January 2019 (Figure 4a-c) revealed that the level of pollution could persist for a couple of months over large areas across the NCP. This period corresponded to the normal conditions of economic activities and winter heating in northern China, which are characteristic of a region with a dense population.
As a result of increases in temperature, NO 2 levels decreased considerably from February to April (Figure 4d-f), which resulted in an improvement of the air quality. In addition, the high NO 2 distribution in the NCP showed considerable reductions in terms of the spatial extent, and only some megalopolis hotspots and isolated hotspots with heavy industries remained.
The spatial coverage areas and changes in NO 2 clearly indicated that most parts of the NCP are subjected to air pollution in the winter. The highest NO 2 mainly appeared in the southern part of Beijing, in most of Tianjin, in the central and south parts of Hebei, and in northeast Shandong, where the tropospheric NO 2 column density maintained high values of over 3 × 10 16 molecules cm −2 . These results indicate that when a large area of pollution recedes, some hot spots of pollution will remain, which can be attributed to vehicle traffic and intensive industrial activities. However, the highest monthly mean NO 2 VCD was found to be less than 2.5 × 10 16 molecules cm −2 .
Atmosphere 2021, 12, 398 7 of 17 A detailed assessment of NO2 pollution also is useful for local-scale analyses to better evaluate the role of pollution sources and pollution changes. For these purposes, NO2 pollution in the NCP (32°N-42°N, 110°E-120°E), where the capital city of China is located, was analyzed in detail ( Figure 4). As a result of increases in temperature, NO2 levels decreased considerably from February to April (Figure 4d-f), which resulted in an improvement of the air quality. In addition, the high NO2 distribution in the NCP showed considerable reductions in terms of the spatial extent, and only some megalopolis hotspots and isolated hotspots with heavy industries remained.
The spatial coverage areas and changes in NO2 clearly indicated that most parts of the NCP are subjected to air pollution in the winter. The highest NO2 mainly appeared in the southern part of Beijing, in most of Tianjin, in the central and south parts of Hebei, and in northeast Shandong, where the tropospheric NO2 column density maintained high values of over 3 × 10 16 molecules cm −2 . These results indicate that when a large area of pollution recedes, some hot spots of pollution will remain, which can be attributed to vehicle traffic and intensive industrial activities. However, the highest monthly mean NO2 VCD was found to be less than 2.5 × 10 16 molecules cm −2 .

Monthly Averaged NO2 Columns of the Four Main Cities
The monthly average columns of NO2 monitored by the EMI for the region within a 0.15 latitude and longitude (~16.7 km × 13 km in the NCP) distance of the four ground stations (Beijing, Baoding, Cangzhou, and Tianjin city) are illustrated in Figure 5.

Monthly Averaged NO 2 Columns of the Four Main Cities
The monthly average columns of NO 2 monitored by the EMI for the region within a 0.15 latitude and longitude (~16.7 km × 13 km in the NCP) distance of the four ground stations (Beijing, Baoding, Cangzhou, and Tianjin city) are illustrated in Figure 5. Temporal variations in data were observed in the EMI NO2 datasets. High NO2 pollution levels were found during winter in the four urban regions. However, the months in which the mean NO2 columns were highest in different regions varied. This variation was likely related to the meteorological conditions (e.g., monsoons) and economic activity. The anomalous high points in Figure 5 represent high values of the tropospheric NO2  Temporal variations in data were observed in the EMI NO 2 datasets. High NO 2 pollution levels were found during winter in the four urban regions. However, the months in which the mean NO 2 columns were highest in different regions varied. This variation was likely related to the meteorological conditions (e.g., monsoons) and economic activity. The anomalous high points in Figure 5 represent high values of the tropospheric NO 2 VCD, which correspond to the most polluted days for the analyzed period. In contrast, the outlier points for the lower points may correspond to weather parameters that favored pollution reductions, such as a particular wind regime or rain wash-out.

Error Analysis
The uncertainty of tropospheric NO 2 inversion mainly includes two parts: (a) the uncertainty of inversion of SCD, (b) the uncertainty of AMF in the VCD calculation. For EMI NO 2 observations, as it involves the separation of NO 2 columns in the stratosphere and troposphere, the NO 2 VCD inversion in the troposphere of the spaceborne platform also includes stratosphere-troposphere separation error. Assuming that these uncertainties are independent, the total uncertainty of the tropospheric NO 2 VCD for MAX-DOAS and EMI can be quantified as follows: With σ DSCD , σ SCDtotal , and σ SCDstrat the slant column error σ AMF , σ dAMF are the error of tropospheric AMF based on EMI and the difference AMF estimation error of ground-based MAX-DOAS, respectively.
For σ DSCD , σ SCDtotal , originate from the DOAS fitting residuals, which can be obtained directly in the spectral fitting. This uncertainty is mainly dominated by the instrument noise, and other systematic errors, such as the uncertainty of fitting the absorption cross-sections and the uncertainty of spectral correction. In this study, the mean vertical errors of fitting MAX-DOAS for Beijing, Tianjin, Baoding, and Cangzhou foundations during the observation period were 1.1 × 10 15 molecules cm −2 , 9.2 × 10 14 molecules cm −2 , 1.2 × 10 15 molecules cm −2 , and 8.6 × 10 14 molecules cm −2 , respectively. σ SCDstrat is stratospheric slant column error. Because this study uses RSM to separate the NO 2 columns in the stratosphere and the troposphere, this part of the error is related to the RSM 7.6 × 10 14 molecules cm −2 [50].
The tropospheric AMF and dAMF were calculated by the radiation transfer model, and depend on a priori gas profile and model parameters (cloud cover, surface reflectivity, etc.). In addition, AMF is also related to the geometric parameters such as the solar zenith angle, the relative azimuth angle and the viewing zenith angle, but the measurement geometry is known with high accuracy and therefore does not contribute significantly to the AMF errors. Therefore, σ AMF , σ dAMF mainly depends on the surface albedo, NO 2 profile, aerosol and cloud cover.
In the comparative study, the NO 2 data of MAX-DOAS and EMI were selected according to cloud cover, and the contribution of cloud cover can be ignored. In addition, most of the input parameters in the AMF simulation come from the chemical transfer model and the weather website, errors will be dampened to some extent. For the impact of the NO 2 profile on AMF, the NO 2 profile in the 2 × 2.5 • geographic grid simulated by the GEOS-Chem model was used as the input parameter of the radiation transfer model, then, simulate the maximum standard deviation of the AMF as the error introduced by the difference in NO 2 profile. In the same way, the perturbation error of the surface albedo to the AMF is also based on the same calculation method. Different types of aerosols and the altitude of the Atmosphere 2021, 12, 398 9 of 17 aerosol layer have different assessments of AMF. In the calculation of AMF, the aerosol optical thickness was assumed to be 0.06. To evaluate the influence of different aerosol scenarios on the AMF under spaceborne and ground-based observation, the actual daily average aerosol data product (https://aeronet.gsfc.nasa.gov/, accessed on 19 March 2021) was used as the input parameter of the radiation transfer model to simulate AMF. The results show that, during the observation period, the influence of aerosols on the AMF of the spaceborne and ground-based observation methods is 0.02% and 0.2%, respectively.
The results in Table 3 provide a general estimate of overall retrieval uncertainties that may be expected for tropospheric NO 2 data under different conditions.

Comparison of EMI and Ground-Based MAX-DOAS Observations
The comparison of the tropospheric NO 2 between the EMI and MAX-DOAS data provided useful information about the changes in NO 2 within the NCP. Figure 6 shows the time series of the tropospheric NO 2 VCD variations derived from the EMI and MAX-DOAS in four cities, and these data cover the period from November 2018 to April 2019. The MAX-DOAS NO 2 vertical columns are shown in their original time resolution (gray dots) as well as averaged 60 min around the EMI overpass (red star points) for all cloud-free days (EMI cloud fraction ≤ 0.5). The latter were used in the comparison for the EMI NO 2 tropospheric vertical columns (green star points). Only the NO 2 results of the cloud-screened qualitycontrolled data meeting the conditions were considered for the comparison.
Through the comparison of the EMI and MAX-DOAS measurements, it was observed that the collocated vertical columns mostly followed the same day-to-day variability. We noted that slight differences occurred between the EMI and ground-based vertical columns, in which the MAX-DOAS results were larger than the EMI results in some cases but not in others, and these findings correspond to the fact that the results did not match perfectly for the EMI and MAX-DOAS measurements. This was expected because of the influence of cloud cover and smog, as well as the differences in spatial representation; in addition, the greater period averaging of ground-based data within 60 min compared to the EMI overpass sampling time could have contributed to the differences.
The EMI tropospheric NO 2 column was statistically compared with datasets derived from MAX-DOAS NO 2 columns from selected potential sources of pollution or contaminated areas. The validation results showed that there was a high correlation between ground-based and satellite-based data for the selected sites during half of a year. Specifically, the correlation coefficients were 0.88, 0.83, 0.78, and 0.85 for Beijing, Tianjin, Baoding, and Cangzhou, respectively, as shown in Figure 7. Moreover, the MD values were positive in Tianjin and Cangzhou, which means that, on average, the EMI slightly overestimated the NO 2 columns in comparison with the ground-based data in Tianjin and Cangzhou. Correspondingly, the MD values for Beijing and Baoding were negative, as the EMI products underestimated the data in Beijing and Baoding in comparison with the ground-based products. When considering the MRD for the four datasets, the biases were all positive and the values were not large, but variations remained within the uncertainties. dots) as well as averaged 60 min around the EMI overpass (red star points) for all cloudfree days (EMI cloud fraction ≤ 0.5). The latter were used in the comparison for the EMI NO2 tropospheric vertical columns (green star points). Only the NO2 results of the cloudscreened quality-controlled data meeting the conditions were considered for the comparison. Through the comparison of the EMI and MAX-DOAS measurements, it was observed that the collocated vertical columns mostly followed the same day-to-day variability. We noted that slight differences occurred between the EMI and ground-based vertical columns, in which the MAX-DOAS results were larger than the EMI results in some cases but not in others, and these findings correspond to the fact that the results did not match perfectly for the EMI and MAX-DOAS measurements. This was expected because of the influence of cloud cover and smog, as well as the differences in spatial representation; in addition, the greater period averaging of ground-based data within 60 min compared to the EMI overpass sampling time could have contributed to the differences.
The EMI tropospheric NO2 column was statistically compared with datasets derived from MAX-DOAS NO2 columns from selected potential sources of pollution or contaminated areas. The validation results showed that there was a high correlation between ground-based and satellite-based data for the selected sites during half of a year. Specifi- cally, the correlation coefficients were 0.88, 0.83, 0.78, and 0.85 for Beijing, Tianjin, Baoding, and Cangzhou, respectively, as shown in Figure 7. Moreover, the MD values were positive in Tianjin and Cangzhou, which means that, on average, the EMI slightly overestimated the NO2 columns in comparison with the ground-based data in Tianjin and Cangzhou. Correspondingly, the MD values for Beijing and Baoding were negative, as the EMI products underestimated the data in Beijing and Baoding in comparison with the ground-based products. When considering the MRD for the four datasets, the biases were all positive and the values were not large, but variations remained within the uncertainties. In fact, according to some other satellite-based and ground-based NO2 product comparisons, the correlation values can be high between the EMI and ground-based MAX-DOAS NO2 observations. For instance, the NO2 columns derived by MAX-DOAS installed at the Royal Meteorological Institute, Brussels, were used to analyze the correlation with In fact, according to some other satellite-based and ground-based NO 2 product comparisons, the correlation values can be high between the EMI and ground-based MAX-DOAS NO 2 observations. For instance, the NO 2 columns derived by MAX-DOAS installed at the Royal Meteorological Institute, Brussels, were used to analyze the correlation with Tropomi products, classified according to seasons, and the correlation between the two datasets reached up to 0.85, while the average value was around 0.7 [51]. Meanwhile, the correlation between the OMI and ground-based MAX-DOAS data was found to be relatively low, with the correlation coefficients in the Mohali area between 0.38 and 0.56 [52].
Overall, the EMI, a state-of-the-art hyperspectral remote sensing instrument of China, can grasp the distribution of NO 2 , and EMI products can play an important role in the analyses of trace gases and the daily variation trends of pollution on a large scale. Meanwhile, ground-based MAX-DOAS measurements have high temporal resolutions, which is of great significance for analyzing the daily and monthly changes of local pollutants.

Application of EMI to NO 2 Pollution Transfer Monitoring
As mentioned in Section 3.1, during winter, the NO 2 amounts in North China are significantly elevated. For the analysis of the tropospheric NO 2 pollution process in the NCP, we considered the following three main factors: the emission of nitrogen oxides, the lifetime of NOx, and the meteorological conditions. First, NOx is produced through coal burning. During winter, heating in the north corresponds to increased coal consumption. Second, due to the low temperatures in winter, the lifetime of NOx is longer than that in other seasons, which delays the progression of atmospheric chemical reactions, slows down the degradation rate of NO 2 , and causes NO 2 build-up in the atmosphere. Finally, worsening meteorological conditions in winter accelerate pollution. Driven by the effects of lower temperatures, the atmospheric boundary layer is lower in winter, and thus, most of the NO 2 remains in the lower troposphere.
In this study, the correlation between the wind field and NO 2 changes was analyzed. The wind field is an important factor leading to NO 2 changes. Under the influence of the winter monsoon, NO 2 can be transported over long distances, which can lead to severe NO 2 pollution events even in areas with few local industrial activities in the NCP. The impact of the wind pattern on NO 2 pollution in the NCP area was analyzed in detail. Statistical methods were used to reveal the sources of pollutant transport at each site. Figure 8 shows the statistical wind rose diagrams between the wind fields and NO 2 columns in Beijing, Tianjin, Baoding, and Cangzhou from November 2018 to April 2019. The direction of the wind field in the rose diagram denotes the real wind direction. Different colored bars represent the daily maximum column of NO 2 measured by MAX-DOAS.
The statistical results illustrate the correlations between the NO 2 transport and the wind field. According to wind field statistics, these four regions are less exposed to easterly weather in winter. When the wind direction was southwest, which is an important transport route, the NO 2 column amount showed a higher value than that for the other wind directions in Beijing and Tianjin. However, the results indicated that the correlation between NO 2 and the wind field was weak in Cangzhou.
To analyze the influence of wind directions on NO 2 pollution, a two-way analysis of variance was adopted. The variation of NO 2 only under the effects of different wind speeds and wind directions was explained. The wind direction was treated as the independent variable and wind speed as the covariate. Table 4 shows the results for the two-way analysis of variance of wind fields in Beijing, Tianjin, Baoding, and Cangzhou based on the wind direction and wind speed. By comparing the p value of the respective variables at a level of 0.05, the statistical analysis revealed that the wind direction had a significant impact on the NO 2 columns in Beijing and Tianjin. tistical methods were used to reveal the sources of pollutant transport at each site. Figure  8 shows the statistical wind rose diagrams between the wind fields and NO2 columns in Beijing, Tianjin, Baoding, and Cangzhou from November 2018 to April 2019. The direction of the wind field in the rose diagram denotes the real wind direction. Different colored bars represent the daily maximum column of NO2 measured by MAX-DOAS. The statistical results illustrate the correlations between the NO2 transport and the wind field. According to wind field statistics, these four regions are less exposed to easterly weather in winter. When the wind direction was southwest, which is an important transport route, the NO2 column amount showed a higher value than that for the other wind directions in Beijing and Tianjin. However, the results indicated that the correlation between NO2 and the wind field was weak in Cangzhou.
To analyze the influence of wind directions on NO2 pollution, a two-way analysis of variance was adopted. The variation of NO2 only under the effects of different wind speeds and wind directions was explained. The wind direction was treated as the independent variable and wind speed as the covariate. Table 4 shows the results for the twoway analysis of variance of wind fields in Beijing, Tianjin, Baoding, and Cangzhou based  The marginal mean of the main effect based on the analysis of variance was computed. The results showed that when the wind direction was approximately 22.5-45 degrees south by west, the maximum marginal mean columns of NO 2 was 8.647 × 10 16 molecules cm −2 . That is, the wind field with this wind direction had the most apparent influence on the NO 2 changes in Beijing. Likewise, the correlation between the change in NO 2 and wind direction was also obvious in Tianjin. However, no apparent evidence that the NO 2 changes have a high correlation with the wind field was found in Baoding and Cangzhou, which indicates that local emissions rather than transport are the most important sources of NO 2 pollution in these two regions.

Local Scale NO 2 Pollution Mapping and Pollution Transport: A Case Study of Beijing and Tianjin
For Beijing and Tianjin, the influence of wind fields on the NO 2 changes in December 2018 was determined by analyzing the NO 2 column amounts and meteorological conditions, including the temperature, humidity, wind speed, and direction at 850 hPa atmospheric pressure.
As mentioned before, high-pollution episodes occur frequently in Beijing and Tianjin in December. Figure 9 shows a time series diagram of the NO 2 changes retrieved by MAX-DOAS and the meteorological conditions at the Beijing and Tianjin sites during 22-25 December 2018. The NO 2 columns exhibited an obvious variation trend. The change process was strongly associated with the wind direction, that is, when the wind direction was a southeast wind, NO 2 showed an increasing trend. The EMI NO2 products also showed good capture of the pollution process. Figure 10 shows the NO2 distribution map and the wind field of the NCP during 22-25 December 2018. The aggregation of NO2 was obvious in Beijing when the wind field was a southwest wind, and this area is surrounded by the Taihang Mountains in the north and west. However, results showed that when the wind field was northwest, there was no obvious pollution transportation from the northeast direction of Beijing. Meanwhile, the NO2 columns measured by the EMI in Beijing and Tianjin were low on those days. Wind fields and NO2 variation reveal a potential pollution transport channel in the southwest direction of Beijing and Tianjin. The EMI NO 2 products also showed good capture of the pollution process. Figure 10 shows the NO 2 distribution map and the wind field of the NCP during 22-25 December 2018. The aggregation of NO 2 was obvious in Beijing when the wind field was a southwest wind, and this area is surrounded by the Taihang Mountains in the north and west. However, results showed that when the wind field was northwest, there was no obvious pollution transportation from the northeast direction of Beijing. Meanwhile, the NO 2 columns measured by the EMI in Beijing and Tianjin were low on those days. Wind fields and NO 2 variation reveal a potential pollution transport channel in the southwest direction of Beijing and Tianjin.

Conclusions
In recent years, the Chinese government has put forth significant efforts toward energy conservation and emission reductions. In this respect, the 12th and 13th Five-Year Plans have been formulated to reduce the impacts of air pollution such as that caused by nitrogen oxides and other greenhouse gases on ecosystems and the environment. However, during particular months, heavy NO2 pollution can still be observed in some areas. Therefore, NO2 pollution monitoring is an important task for identifying emissions and developing effective reduction strategies.
The recently launched Environmental Trace Gases Monitoring Instrument (EMI) payload on the GF-5 satellite allows for the monitoring of trace gases, and these data can be used to analyze atmospheric pollution on different scales from large cities to global coverage. In this study, we used the new satellite EMI payload to observe the NO2 amounts in North China from November 2018 to April 2019, in combination with groundbased multi-axis differential optical absorption spectroscopy (MAX-DOAS) observation data and numerical models. The main conclusions can be summarized as follows: (a) Using six months' worth of data over the North China Plain (NCP), the correlation coefficient of tropospheric NO2 between the EMI and MAX-DOAS located at four sites reached 0.88. The spatial distribution and changes in tropospheric NO2 of EMI agreed well with ground-based observations. The good agreement of NO2 products shows that the EMI can be used to accurately monitor tropospheric NO2 pollution.
(b) The distribution of NO2 in China showed obvious temporal and spatial variation, in which higher amounts of pollution were detected in winter. The NO2 in the NCP from November to January increased owing to the large quantities of fossil fuel consumption in North China for heating.

Conclusions
In recent years, the Chinese government has put forth significant efforts toward energy conservation and emission reductions. In this respect, the 12th and 13th Five-Year Plans have been formulated to reduce the impacts of air pollution such as that caused by nitrogen oxides and other greenhouse gases on ecosystems and the environment. However, during particular months, heavy NO 2 pollution can still be observed in some areas. Therefore, NO 2 pollution monitoring is an important task for identifying emissions and developing effective reduction strategies.
The recently launched Environmental Trace Gases Monitoring Instrument (EMI) payload on the GF-5 satellite allows for the monitoring of trace gases, and these data can be used to analyze atmospheric pollution on different scales from large cities to global coverage. In this study, we used the new satellite EMI payload to observe the NO 2 amounts in North China from November 2018 to April 2019, in combination with ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) observation data and numerical models. The main conclusions can be summarized as follows: (a) Using six months' worth of data over the North China Plain (NCP), the correlation coefficient of tropospheric NO 2 between the EMI and MAX-DOAS located at four sites reached 0.88. The spatial distribution and changes in tropospheric NO 2 of EMI agreed well with ground-based observations. The good agreement of NO 2 products shows that the EMI can be used to accurately monitor tropospheric NO 2 pollution.
(b) The distribution of NO 2 in China showed obvious temporal and spatial variation, in which higher amounts of pollution were detected in winter. The NO 2 in the NCP from November to January increased owing to the large quantities of fossil fuel consumption in North China for heating.
(c) Based on observations and statistical analysis, the wind field is one of the important factors affecting the distribution of NO 2 pollution in the NCP region. The southwest wind in Beijing and Tianjin has a correlation with the increase of NO 2 to a certain degree. This research is a good demonstration of how a comprehensive understanding of regional NO 2 pollution can be obtained through joint observations of the EMI and ground-based MAX-DOAS, and this technique will be applied to more research projects in the future.

Data Availability Statement:
The EMI and MAX-DOAS data used in this study can be made available upon request (sifuqi@aiofm.ac.cn).