A new TROPOMI product for tropospheric NO2 columns over East Asia with explicit aerosol corrections

We present a new product with explicit aerosol corrections, POMINO-TROPOMI, for tropospheric nitrogen dioxide (NO2) vertical column densities (VCDs) over East Asia, based on the newly launched TROPOspheric Monitoring Instrument with an unprecedented high horizontal resolution. Compared to the official TM5-MP-DOMINO (OFFLINE) product, POMINO-TROPOMI shows stronger concentration gradients near emission source locations and better agrees with MAX-DOAS measurements (R2 = 0.75; NMB= 0.8 % versus R2 = 0.68, NMB=−41.9 %). Sensitivity tests suggest that implicit aerosol corrections, as in TM5-MP-DOMINO, lead to underestimations of NO2 columns by about 25 % over the polluted northern East China region. Reducing the horizontal resolution of a priori NO2 profiles would underestimate the retrieved NO2 columns over isolated city clusters in western China by 35 % but with overestimates of more than 50 % over many offshore coastal areas. The effect of a priori NO2 profiles is more important under calm conditions.

POspheric Monitoring Instrument (TROPOMI), which was jointly developed by the Netherlands and Europe Space Agency (ESA) (Veefkind et al., 2012) and was launched on 13 October 2017, is a UV-visible (Vis)-near-infraredshortwave infrared backscattering sensor on board the sunsynchronous Sentinel-5 Precursor (S5P) satellite with an overpass time of 13:30 local solar time. With a wide swath of 2600 km and an unprecedentedly high horizontal resolution of 3.5 km×7 km, (3.5 km×5.5 km since 6 August 2019) TROPOMI achieves daily global coverage. This high horizontal resolution and good spatial coverage, combined with the very high signal-to-noise ratio, enables the instrument to resolve NO 2 pollution from point sources, medium-size urban centers, highways, or rivers, tasks that were difficult to achieve before.
Retrievals of tropospheric NO 2 vertical column densities (VCDs) in the UV-Vis spectral range from satellite instruments consist of three steps: (1) using differential optical absorption spectroscopy (DOAS) to fit the slant column density (SCD) of NO 2 along the light path, (2) subtracting the stratospheric contribution from the SCD to obtain the tropospheric SCD, and (3) converting the tropospheric SCD to the tropospheric VCD by using the calculated air mass factor (AMF). For TROPOMI, the random uncertainty in the total SCDs is ∼ 0.6 × 10 15 molec. cm −2 , considerably smaller than for the Ozone Monitoring Instrument (OMI; ∼ 0.8 × 10 15 molec. cm −2 ; Zara et al., 2018). The (total or stratospheric) bias is generally between 0 % and −10 % according to SAOZ observations , which meets the error requirements as defined in the S5P Calibration and Validation Plan (Goryl et al., 2017). The calculation of the AMF introduces the dominant source of error in the retrieved tropospheric NO 2 VCDs over polluted areas (Boersma et al., 2004(Boersma et al., , 2011Lorente et al., 2017;Lin et al., 2014). The median negative biases of the daily comparisons between tropospheric VCDs from the Dutch official TM5-MP-DOMINO (OFFLINE) product and MAX-DOAS measurements are generally less than 50 % (the allowable bias is 25 %-50 %; Goryl et al., 2017) but quite variable with the stations and NO 2 levels, especially at polluted locations van Geffen et al., 2019).
TM5-MP-DOMINO uses implicit aerosol corrections by assuming aerosols to be "effective clouds", as assumed in most satellite NO 2 products except POMINO (Liu et al., 2019;Lin et al., 2014Lin et al., , 2015. In addition, TM5-MP-DOMINO employs a priori NO 2 profiles from the TM5 model at a relatively coarse horizontal resolution (1 • × 1 • ; Williams et al., 2017). The implicit aerosol corrections (Lin et al., 2014;Lorente et al., 2017;Liu et al., 2019) and the coarse horizontal resolution of a priori NO 2 profile data (Laughner et al., 2016;Russell et al., 2011;McLinden et al., 2014) may be the largest sources of the large biases observed between TM5-MP-DOMINO and MAX-DOAS. Based on previous studies for OMI, implicit aerosol corrections can lead to uncertainties of more than 50 % over polluted ar-eas with high aerosol loadings like China (Lin et al., 2014;Lorente et al., 2017;Liu et al., 2019). Eskes et al. (2018) showed that using a priori NO 2 profiles from the regional CAMS model at a 12 km × 12 km resolution to replace the TM5 NO 2 profiles increases the retrieved NO 2 VCDs by ∼ 0 % to 50 % over western Europe depending on the location.
Here we present a new TROPOMI tropospheric NO 2 VCD product over East Asia, namely POMINO-TROPOMI. This product is based on our POMINO algorithm (Liu et al., 2019;Lin et al., 2014Lin et al., , 2015 previously applied to OMI. POMINO-TROPOMI improves upon TM5-MP-DOMINO by employing explicit aerosol corrections and using high-resolution (∼ 25 km) a priori NO 2 profiles, among other improvements. POMINO-TROPOMI NO 2 VCD data over July-October 2018 are presented and validated by MAX-DOAS measurements, along with additional sensitivity tests of the effects of aerosol representations and a priori NO 2 profiles.
2 Method and data 2.1 POMINO-TROPOMI retrieval algorithm and product As one of the UV-Vis backscatter instruments to observe NO 2 , TROPOMI inherits much of the design of OMI (Veefkind et al., 2012). Thus the POMINO-TROPOMI algorithm here largely follows our previous POMINO algorithm (Liu et al., 2019), with some modifications to adapt to its high horizontal resolution and different cloud retrieval procedure. The POMINO-TROPOMI algorithm focuses on improving the calculation of tropospheric AMF. It thus takes the tropospheric SCD data from TM5-MP-DOMINO (OFFLINE), which are obtained by fitting the 405-465 nm wavelength range with the DOAS method. Our tropospheric AMF calculation is done for 437.5 nm, following TM5-MP-DOMINO.
We use the parallelized linearized discrete ordinate radiative transfer (LIDORT)-driven AMFv6 package to derive tropospheric AMFs via online pixel-specific radiative transfer calculations, with no use of lookup tables. Our algorithm explicitly accounts for aerosol optical effects and anisotropic properties of surface reflectance, uses daily a priori aerosol and NO 2 profiles from the simulation of nested GEOS-Chem model (0.25 • lat × 0.3125 • long; Zhang et al., 2016), and further uses aerosol optical depth (AOD) data from Moderate Resolution Imaging Spectroradiometer (MODIS/Aqua) to correct GEOS-Chem simulated aerosols on a monthly basis. Figure 1 shows the procedure of using the AMFv6 package to derive the tropospheric NO 2 VCDs of POMINO-TROPOMI. Table S1 in the Supplement shows the retrieval parameters in POMINO-TROPOMI, in comparison with those in TM5-MP-DOMINO and POMINO v2.
The independent pixel approximation (IPA) is used to calculate AMF as a linear combination of a cloudy AMF (M cld ) (1) w is the cloud radiation fraction (CRF) calculated by where I cld denotes the radiance from the cloudy part of the pixel, I clr the radiance from the clear-sky part of the pixel, f eff the cloud fraction (CF), and R the total scene radiance. Retrieval of cloud properties is a prerequisite for NO 2 retrieval. We take the effective cloud pressure (CP) from the FRESCO-S algorithm (van Geffen et al., 2019), which uses the O 2 A band (around 758 nm) for TROPOMI trace gas retrievals. We recalculate w and f eff at the NO 2 fitting wavelength (437.5 nm). The online CF calculation is similar to that for TM5-MP-DOMINO (Arnoud et al., 2017;van Geffen et al., 2019) but with an explicit aerosol correction to be consistent with the following NO 2 retrieval and to remove the aerosol signal from the retrieved CF data. For explicit aerosol corrections in this study, we take daily aerosol simulation results from the GEOS-Chem v9-02 nested model over East Asia, followed by a monthly AOD correction using MODIS/Aqua C6.1 AOD data. Our future study will use the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aerosol extinction vertical profiles to further improve the modeled aerosol profiles. Figure 2b shows the spatial distribution of AOD in July 2018 used in clouds and NO 2 retrievals. The AOD distribution is consistent (R = 0.42; N = 1447) with that of near-surface PM 2.5 mass concentration measurements (Fig. 2a) taken from the Ministry of Ecology and Environment of China (MEE); the difference between AOD and near-surface PM 2.5 is expected because they represent different parameters of aerosols. The criteria to select valid pixels in this study are as follows. We exclude pixels with viewing zenith angles (VZAs) greater than 80 • , with high albedos caused by ice or snow on the ground, or with a quality flag (from TM5-MP-DOMINO) of less than 0.5. To screen out cloudy scenes, we discard pixels with CRFs greater than 50 % in the POMINO-TROPOMI product.
In addition to our formal POMINO-TROPOMI product (referred to as Case REF), we use sensitivity cases to evaluate the impacts of aerosol corrections (explicit versus implicit) and the horizontal resolution of a priori NO 2 profiles (Cases 1 and 2 in Table 1). Two additional cases (Cases 3 and 4) concern the treatment of CP in combination with the choice of aerosols and surface reflectance. Specifically, using the CP data directly from FRESCO-S means that our retrieval algorithm does not perfectly account for the effect of aerosols on clouds. Our retrievals consider the anisotropy of the surface via (bidirectional reflectance distribution function) BRDF effects, while a Lambertian surface is used in  deriving the FRESCO-S CP. To ensure sampling consistency, the pixels used in all cases are selected based on the CRF values in Case REF.

Ground-based MAX-DOAS measurements
We use ground-based MAX-DOAS measurements to validate the POMINO-TROPOMI NO 2 product. The MAX-DOAS measurements are from two suburban stations (Xuzhou and Nanjing) and one remote station (Fukue; Kanaya et al., 2014). Table 2 shows the geographical and time information of each MAX-DOAS site, and Sect. S1 in the Supplement describes each instrument in detail. Although Xuzhou and Nanjing are both classified as suburban sites located at university campuses, the NO 2 spatial distributions around the two sites are very different. The spatial distribution of NO 2 VCDs is relatively smooth around Xuzhou, whereas the VCDs exhibit a strong spatial gradient around Nanjing (Fig. 3).
A consistent spatiotemporal sampling is crucial in comparing satellite measurements and MAX-DOAS data Lin et al., 2014;Liu et al., 2019;Wang et al., 2017). We average all valid MAX-DOAS data within ±1 h of the TROPOMI overpass time to obtain daily values for comparison. To reduce the influence of local events, we exclude all MAX-DOAS data whose standard deviations within the 2 h exceed 20 % of their mean values. We average all valid pixels within 5 km of each MAX-DOAS site to represent the respective daily satellite data. Section S2 shows how the validation results are affected by the sampling choice.  Fig. 4a visualize the near-surface NO 2 concentrations observed at the MEE sites at the overpass time of TROPOMI. In both the VCD and the near-surface concentration maps ( Fig. 4a and b), hotspots like urban centers and isolated sources can be seen clearly, due to the short lifetimes of NO x in summer. The spatial correlation is about 0.55 (N = 1458) between the VCD and the near-surface concentration distributions. Figure 4c shows the spatial distribution of TM5-MP-DOMINO (OFFLINE) NO 2 VCDs for comparison. The general distribution of TM5-MP-DOMINO is consistent with that of POMINO-TROPOMI with a correlation coefficient of 0.97 (N = 1 091 154). However, TM5-MP-DOMINO NO 2 VCD values are lower than POMINO-TROPOMI by about 35 % averaged over the whole domain (Fig. 4d), by −37 %-68 % over cleaner areas (POMINO-TROPOMI < 5 × 10 15 molec. cm −2 ), and by 0 %-66 % over more polluted areas (POMINO-TROPOMI ≥ 5 × 10 15 molec. cm −2 ). TM5-MP-DOMINO does not show strong local signals at pollution hotspots like the urban center of Beijing (Fig. 4c). Over these hotspot locations, TM5-MP-DOMINO is lower than POMINO-TROPOMI by up to 5 × 10 15 molec. cm −2 (Fig. 4e). Figure 5a and b present the two products over Beijing and surrounding areas, showing a much weaker spatial gradient of NO 2 VCDs from Beijing urban center to its outskirts in  TM5-MP-DOMINO than in POMINO-TROPOMI. The corresponding histograms and Gaussian fittings in Fig. 5c also show a lower mean value and a smaller standard deviation of TM5-MP-DOMINO than POMINO-TROPOMI. These results highlight the important differences between the two products at fine scales.

POMINO-TROPOMI NO 2 VCDs over East Asia
We further compare the two satellite products with ground-based MAX-DOAS measurements at three sites. The scatterplots in Fig. 6a and b compare the NO 2 VCDs on 63 d (from 109 pixels) over July-October 2018 with their MAX-DOAS counterparts. Different colors differentiate the sites. POMINO-TROPOMI captures the day-to-day variability in MAX-DOAS data (R 2 = 0.75) and shows a small normalized mean bias (NMB = 0.8 %). The reduced major axis (RMA) regression shows a slope of 0.70, mainly because of the underestimate on high-NO 2 days. TM5-MP-DOMINO is also correlated with MAX-DOAS (R 2 = 0.68), although the correlation is weaker than our retrieval. The NMB of TM5-MP-DOMINO is much more significant (−41.9 %), and the RMA regression slope is much smaller (0.42). Similar underestimates of TM5-MP-DOMINO have been discussed in their "Readme" document  and the algorithm theoretical basis document (ATBD) file (van Geffen et al., 2019) in general, in Griffin et al. (2019) for Canada. Major plausible causes of the underestimate of TM5-MP-DOMINO include coarse-resolution climatological surface albedo data, coarse-resolution (1 • ×1 • ) a priori NO 2 profiles, implicit aerosol corrections, and uncertainties in CP from FRESCO-S.
Cloud pressures from FRESCO-S are found to be too high -i.e., the cloud top is too close to the ground, especially over China (van Geffen et al., 2019). We examine this effect by excluding the pixels with CP > 850 hPa when comparing with MAX-DOAS data. With this additional criterion, the number of valid days drops to 20. Figure 6c and d show the scatterplots and corresponding RMA results. The NMB of TM5-MP-DOMINO is reduced slightly to −38.2 %, and its R 2 for day-to-day variation is increased from 0.68 to 0.85. POMINO-TROPOMI still outperforms TM5-MP-DOMINO: 0.85 versus 0.85 for R 2 , 13.8 % versus −38.2 % for NMB, and 0.82 versus 0.56 for RMA regression slope. The improvement from excluding CP > 850 hPa scenes is larger in TM5-MP-DOMINO (with implicit aerosol corrections) than in POMINO-TROPOMI (with explicit aerosol corrections). The averaged CF of data excluding CP > 850 hPa is 0.13 (AOD = 0.63), which is much larger than the averaged value (CF = 0.06) in Fig. 6a and b (AOD = 0.57). This appears to imply that the overestimated CP may be partly because the FRESCO-S cloud algorithm might misinterpret heavy aerosol loadings near the ground as clouds, a common issue in satellite data (Lin and Li, 2016).

Influences of aerosol correction approaches and horizontal resolutions of a priori NO 2 profiles
To investigate the causes of difference between POMINO-TROPOMI and TM5-MP-DOMINO (OFFLINE), we conduct two sensitivity retrievals based on the POMINO-TROPOMI algorithm (Cases 1 and 2 in Table 1). Figure 7 shows the relative (panels a-c) and absolute (panels d-f) dif-ferences between retrieval cases (REF, Case 1, and Case 2). Case 1 adopts the implicit aerosol correction for both clouds and NO 2 retrievals, as in TM5-MP-DOMINO, so the difference between REF and Case 1 indicates the effect of aerosol representation (explicit versus implicit) (Fig. 7a, d).
In Case 2, the high-resolution (0.25 • lat × 0.3125 • long) NO 2 profiles are replaced by low-resolution (2 • lat × 2.5 • long) profiles simulated by the GEOS-Chem global model; aerosols are represented implicitly as in Case 1. Case 2 thus mimics TM5-MP-DOMINO, which uses an implicit aerosol correction and coarse-resolution NO 2 profiles. Thus, the difference between Case 1 and Case 2 arises from the a priori NO 2 profiles (Fig. 7b, e). The difference between Case 2 and REF further indicates the joint effect of using coarseresolution a priori NO 2 profiles and an implicit aerosol correction (Fig. 7c, f). Figure 7 shows that the individual influences of aerosol representations (explicit versus implicit) and a priori NO 2 profiles (fine versus coarse resolution) vary substantially from one location to another. The impacts of aerosol corrections are most evident over the areas of heavy aerosol loadings including East China, India, and parts of Southeast Asia. Figure 6. Scatterplot for daily NO 2 VCDs (10 15 molec. cm −2 ) between MAX-DOAS and two TROPOMI NO 2 VCD products. Each colored dot represents a day, and each color denotes a station. For each day, the satellite data are averaged over all pixels. Panels (c) and (d) are results of two TROPOMI NO 2 VCD products with effective cloud pressures ≤ 850 hPa.
The implicit aerosol correction (Case 1) tends to result in lower NO 2 VCDs by 0 %-50 % over urban areas compared to the explicit aerosol correction (Case REF) (Fig. 7a, d). By comparison, the impacts of NO 2 profiles are spatially more heterogeneous ( Fig. 7b and e). Over the offshore coastal areas, using coarse-resolution NO 2 profiles (Case 2) tends to overestimate the NO 2 VCDs by 30 %-100 % relative to when high-resolution profiles are used (Case 1), due to (horizontal) over-smoothing of NO 2 .
Below, we discuss these differences over two key areas including northern East China and Xinjiang. Northern Eastern China (bounded by the black rectangle in Fig. 7a) is a region with heavy aerosol loading region (Fig. 2). Over this region, an implicit rather than explicit aerosol representation results in lower NO 2 VCDs by ∼ 25 % (Fig. 7a, d), while the effect of NO 2 profiles is weaker (Fig. 7b, e). The joint effect is dominated by the effect of aerosol representation (Fig. 7c, f).
That the impact of a priori NO 2 profiles is relatively small over northern Eastern China is partly because of the smearing effect by wind. Figure 8 differentiates the effects of NO 2 profiles between 23 windy (daily average wind speed under 500 m > 2 m s −1 ) days and 7 calm (wind speed < 2 m s −1 ) days. The dataset of wind is taken from the National Aeronautics and Space Administration-Global Modeling and Assimilation Office's (NASA-GMAO's) "forward-processing" (GEOS-FP) data product with the horizontal resolution at 0.25 • lat × 0.3125 • long. In the windy cases, the NO 2 VCDs are much more smoothed even at 0.3125 • resolution, and thus the difference of NO 2 resolutions is very small. In contrast, NO 2 VCDs exhibit much stronger horizontal gradient in calm situations, which are better retrieved when highresolution NO 2 profiles are used. As calm situation is more favorable for pollutant accumulation, while windy days help to dilute concentration of NO 2 , so much higher NO 2 VCDs are found in Fig. 8a and b. Xinjiang (bounded by the red rectangle in Fig. 7a) is a deserted region in West China. Over Xinjiang, the resolution of a priori NO 2 profiles affects the retrieved NO 2 VCDs much more than the aerosol representation does (Fig. 7b and e). Compared to Case 1 (with high-resolution NO 2 profiles), Case 2 (with coarse-resolution profiles) leads to much lower NO 2 VCDs over the isolated urban areas which are not resolved by the coarse model.

Influences of directly using the CP data from FRESCO-S
As we take the CP data directly from the FRESCO-S retrieval rather than re-retrieving CP (as done for CF), two main issues arise. First, the FRESCO-S retrieved CP may be affected by aerosols, and thus using such CP data in our explicit aerosol corrections (Case REF) may lead to overcorrection of aerosol effects. To estimate the effect of such overcorrection on retrieved NO 2 VCDs, we employ in an additional sensitivity case (Case 3 in Table 1) a "semi-explicit" aerosol correction approach. This approach explicitly includes aerosols in the calculation of AMF for the clear-sky portion (M clr ) of a pixel (as in Case REF) but excludes aerosols for the cloudy-sky portion (M cld ) of that pixel. Correspondingly, CF is recalcu-lated on the basis that the radiance at 437.5 nm received by TROPOMI is contributed from the aerosol-contained clearsky portion and the no-aerosol, cloudy-sky portion. Table 3 shows that in July 2018, on a pixel basis, the derived NO 2 profiles in Case 3 are larger than those in Case REF, with an average difference increasing from 3.1 % at relatively clean situations (NO 2 VCDs in Case REF < 5×10 15 molec. cm −2 ) to 11.2 % for polluted situations (NO 2 VCDs in Case REF ≥ 15 × 10 15 molec. cm −2 ). The spatial distributions in Fig. S1a and S1b also show higher NO 2 VCDs in Case 3 relative to Case REF. The corresponding increases in CF (Fig. S1c versus Fig. S1d) are because in Case 3 the scattering contributions to the radiance from aerosols in the cloudy-sky portion (that would have occurred) are accounted for with higher CFs. The enhanced "shielding effect" of clouds (due to higher CFs) result in lower NO 2 AMFs and higher VCDs.  For surface reflectance, Case REF considers the BRDF effect instead of Lambertian equivalent reflectivity (LER) which is used in deriving FRESCO-S clouds. This leads to inconsistency between when the CP is derived and when it is used. The LER data used by FRESCO-S are gener-ated at 758 and 772 nm (based on the Global Ozone Monitoring Experiment-2), rather than at the 437.5 nm for the NO 2 retrieval. Thus Case 4 adopts the OMI LER data from TM5-MP-DOMINO -a 5-year, monthly-based climatology at 440 nm -and recalculates CFs and NO 2 AMFs with ex-  Lin et al. (2010), Lin et al. (2014), Boersma et al. (2004Boersma et al. ( , 2011Boersma et al. ( , 2018 and this study. The error accounts for the effect of horizontal resolutions and the vertical process in GEOS-Chem. b This accounts for the expectation that our explicit aerosol correction leads to more reasonable CFs. c Based on van Geffen et al. (2019). Our estimated NO 2 errors related to the use of FRESCO CP data (instead of recalculating it) are within this error range. d The estimated parameter error is taken from Zhou et al. (2010), and the corresponding AMF error is based on Lin et al. (2014) and Case 4. e Based on Tian et al. (2018), which compared the MODIS Merged AOD C6.1 product with AERONET in China. f Based on Lin et al. (2015), which compared GEOS-Chem simulated SSA with Lee et al. (2007). g Based on Liu et al. (2019), which compared GEOS-Chem simulated aerosol extinction profiles with CALIOP. h This is a tentative error estimate based on Lin et al. (2015), Lorente et al. (2017), Liu et al. (2019), and this study. With explicit aerosol corrections, the errors in heavily polluted situations are expected to be smaller than when assuming implicit aerosol corrections.
plicit aerosol corrections and FRESCO-S CP (Table 1). Here, the ice/snow flag in the TM5-MP-DOMINO product is used to exclude the possible ice/snow contamination, and only the pixels with blue-sky albedos (derived from the BRDF data in Case REF) less than 0.3 are taken into consideration. The resulting NO 2 VCDs in Case 4 are lower than Case REF by 3.7 % on a pixel-based average for relatively clean situations and by 8.3 % for polluted situations (Table 3). Figure S2a and S2b further show the spatial distributions of the relative and absolute differences in derived monthly mean NO 2 VCDs between Case 4 and Case REF.  S2d).

Error estimates for POMINO-TROPOMI NO 2 AMFs
It is difficult to derive an overall AMF error for each pixel with our algorithm, particularly because of the threedimensional aerosol parameters used, the interlinkage between clouds and other parameters (aerosols and surface reflectance), and the exclusion of lookup tables (LUTs; which leads to too computationally expensive error estimates). Table 4 provides a preliminary estimate of the NO 2 AMF errors with respect to uncertainties in individual parameters. We follow the ATBD of TM5-MP-DOMINO ( van Geffen et al., 2019) and make use of error estimates in previous studies. Individual errors with respect to CF, CP, and BRDF coefficients are within 10 %. Errors with respect to aerosols are consid-ered to be larger, due to uncertainty in AOD, single-scattering albedo (SSA), and aerosol vertical profiles (Liu et al., 2019), as well as the fact that using the CP from FRESCO-S rather than deriving it here may lead to an additional error in the NO 2 AMFs. Note that these individual errors are not fully independent due to the aforementioned linkage between parameters. The magnitude of potential systematic bias in NO 2 may be smaller than the quadrature sum of errors in individual parameters, as suggested by the slight mean bias relative to MAX-DOAS data (Fig. 6a).

Conclusion and discussion
The POMINO algorithm to retrieve tropospheric NO 2 VCDs has been successfully applied to TROPOMI over East Asia. The resulting POMINO-TROPOMI product shows higher tropospheric NO 2 VCDs (by about 35 % averaged over East Asia) and much clearer urban and other hotspot signals, compared to the TM5-MP-DOMINO (OFFLINE) product. Further evaluation using independent MAX-DOAS measurements indicates very good performance of POMINO-TROPOMI in capturing the day-to-day variation (R 2 = 0.75; N = 63) and mean value (NMB = 0.8 %) of NO 2 , better than TM5-MP-DOMINO (0.68 and −41.9 %, respectively).
Over regions with heavy aerosol loading, the accuracy of retrieved NO 2 VCDs is affected substantially by how aerosols are represented in the retrieval process (implicit or explicit). The implicit aerosol representation tends to underestimate the NO 2 VCDs by 0 %-50 % over most urban areas in East Asia and by about 25 % over northern East China. Using a priori NO 2 profile data at a horizontal resolution of ∼ 25 km, POMINO-TROPOMI captures the cityscale NO 2 hotspots. Reducing the horizontal resolution of a priori profiles to what is typically set up by global models (100-200 km) underestimates the NO 2 hotspots and the spatial gradient surrounding them, and the effects are more pronounced in calm than in windy situations. Overall, our results provide useful information to improve TROPOMI retrieval algorithms and offer insight for applications to the upcoming geostationary satellite instruments including GEMS, TEMPO, and Sentinel-4 Precursor.
Further work can be done to improve the retrieval algorithm for TROPOMI. First, the hybrid cloud retrieval method used in POMINO-TROPOMI is not optimal, because only cloud fraction but not cloud pressure is recalculated with explicit aerosol corrections and BRDF effects. The uncertainty caused by inconsistent assumptions of aerosols and albedos in cloud pressure and NO 2 VCD retrievals are initially estimated in this study. Recalculation of cloud pressure will be done in the future using the O 2 -O 2 method when the O 2 -O 2 SCD data are available. Second, correcting the simulated aerosol extinction vertical profiles will further improve the clouds and NO 2 retrievals (Liu et al., 2019). Third, the horizontal resolution of a priori NO 2 profiles (∼ 25 km) does not match the fine footprint of TROPOMI, and further increasing the resolution will help achieve better accuracy for analyses of very fine scale pollution characteristics such as along the highways and rivers and within urban centers.
Data availability. The POMINO-TROPOMI NO 2 data are available at the ACM group product website (http: //www.pku-atmos-acm.org/acmProduct.php/, Lin et al., 2020). The TM5-MP-DOMINO product can be download via the TEMIS website at http://www.temis.nl/airpollution/no2.html Royal Netherlands Meteorological Institute, 2020). The near-surface data of NO 2 and PM 2.5 can be downloaded from http://www.cnemc.cn/sssj/cskqzl/ (Ministry of Ecology and Environment, 2019). The ground-based MAX-DOAS observations would be provided after the applications of users are approved by corresponding owners.
Author contributions. JL conceived the research. ML and JL designed the experiment. ML performed all calculations with some code support from HK and HW. ML and JL wrote the paper with inputs from KFB and HE. RS provided LIDORT. YK, QH, XT, KQ, and PX provided the MAX-DOAS observations. RN helped to process surface observations. YY and JW helped to analyze the relationship between meteorological and NO 2 VCD variations. All authors commented on the paper.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "TROPOMI on Sentinel-5 Precursor: first year in operation (AMT/ACP inter-journal SI)". It is not associated with a conference.
Financial support. This research has been supported by the National Natural Science Foundation of China (grant no. 41775115) and the second Tibetan Plateau Scientific Expedition and Research Program (grant no. 2019QZKK0604).
Review statement. This paper was edited by Ben Veihelmann and reviewed by two anonymous referees.