2.1 Simulations using the original emissions inventory
We evaluated the model accuracy achieved using the original extent of NOx emissions in the inventory. This involved comparing the tropospheric NO2 columns and surface NO2 concentrations observed during the months of March, April, and May (MAM) 2022 to those simulated by CMAQ using the a priori NOx emissions.
Figures 1a and 1b illustrate the monthly averages of hourly NO2 columns observed and modeled during daytime. Compared to GEMS NO2 columns, the model tended to overestimate the columns across most regions of South Korea (hereafter referred to as Korea), with the exception of the Seoul metropolitan area (SMA), and in densely populated urban areas in China, including Beijing and Shenyang, and those across the southeast region. The model also overestimated the columns in major urban centers in Southeast Asian countries (e.g., Ho Chi Minh City and Hanoi in Vietnam, Bangkok in Thailand, and Kuala Lumpur in Malaysia, and Singapore), in India (particularly in the northeastern region), and in Japan (except for Tokyo) (see Figure S1 for geographic labels). Conversely, the model generally underestimated the columns in the SMA, across the North China Plain (NCP) and the plains in the northeast region of China, and in the landlocked urban centers of northcentral and southwestern China. The model underestimated the columns in the Yangtze River Delta (YRD) in March, and then overestimated the columns in April and May. Similar patterns of monthly spatial discrepancies, but to a greater extent, were noted between the NO2 columns observed at 04:45 UTC (hereafter referred to as LEO proxy NO2 columns; see Methods) and the corresponding modeled NO2 columns (Figs. 1d and 1e).
Figure 2 shows a time series of the observed and modeled hourly surface NO2 concentrations at ground-based monitoring stations in China during MAM 2022. The model generally underestimated the concentrations during daylight hours, aligning with GEMS retrieval times, and began overestimating in the subsequent hours leading up to nighttime. A similar tendency was observed in Korea (Figure S1), indicating a regional modeling challenge in East Asia when using the a priori emissions. Overall, there was a fair agreement between the observed and modeled hourly NO2 concentrations, with correlation coefficients (R) ranging from 0.57 to 0.64 and indices of agreement (IOA) from 0.64 to 0.74 in China (Fig. 2) and those from 0.62 to 0.82 and 0.65 to 0.71, respectively, in Korea (Figure S1). The extent of overestimation was noticeable throughout the months, with normalized mean biases (NMB) ranging from 7.75–32.59%, mainly due to nighttime overestimation. However, during the daytime, the model severely underestimated NO2 concentrations, with NMB from − 13.17% to -29.31% in Korea and from − 5.08% to -23.58% in China. This suggests a potential underrepresentation of daytime NOx emissions in Asia, as discussed in several previous studies, which led to the overall underestimation of surface NO2 concentrations in East Asia27,34,36,52.
2.2. Top-down updates to the NOx emissions inventory
After the top-down updates to the NOx emissions inventory, we examined the subsequent changes in NOx emissions in Asia during MAM 2022. This involved comparing the a priori emissions with the a posteriori NOx emissions, which were obtained through the Bayesian inversions informed by the observation data from GEMS and LEO proxy (see Methods).
Figure 3 shows the monthly averages of hourly daytime NOx emissions, comparing the updated emissions to the a priori emissions. The GEMS-informed update led to spatial adjustments in emission quantities, aiming to offset the model’s prior biases discussed earlier. NOx emissions generally decreased in areas where the model previously underestimated the NO2 columns, including most regions across Korea (except for the SMA), Beijing, Shenyang, the YRD (in April and May), as well as in highly populated urban areas in southeastern China, Southeast Asian countries, India, Japan (except for Tokyo). Specifically, in Beijing and Shenyang, the decreases in NOx emissions at the pixels corresponding to the city centers were overwhelmed by the increased emissions in the surrounding greater metropolitan areas. Conversely, we observed substantial increases in NOx emissions in the SMA, the NCP, the plains in northeastern China, the YRD (in March), and urban centers in northcentral and southwestern China. In areas where the model previously underestimated and overestimated NO2 columns (Fig. 1c), we observed an average increase of 31.72% and a decrease of 22.78% in NOx emissions, respectively, during MAM (Fig. 3c). The spatial extent of areas with decreased emissions was slightly smaller than that with increased emissions by a factor of 0.998, which led to a slight increase in domain-averaged NOx emissions by 1.42%. This suggests that the localized increases in East Asia slightly outweighed the reductions elsewhere.
Upon observing similar spatial distributions, the LEO-informed update led to notably larger increases and decreases in NOx emissions. In areas where the model previously underestimated and overestimated NO2 columns (Fig. 1f), we noted an average increase of 35.39% and a decrease of 36.40% in NOx emissions, respectively (Fig. 3e). The total area with decreased emissions was slightly larger than that with increased emissions by a factor of 1.023, resulting in an average decrease of 17.10% in NOx emissions. This tendency towards more extensive adjustments, particularly towards decreasing the emissions, is attributed to the conventional inversion method practiced with LEO proxy NO2 columns, which relies on monthly averaged, time-specific discrepancies between observed and modeled NO2 columns (Fig. 1f) to adjust NOx emissions throughout the entire day.
Figure 4 shows a time series of the monthly averages of hourly daytime NOx emissions in Korea and China during MAM 2022 before and after the updates. We compared the absolute values of the a priori and a posteriori NOx emissions and the extents of adjustments at each hour. This enabled us to assess how effectively the GEMS- and LEO-informed inversions captured daytime diurnal variations in the emissions.
The GEMS-informed update generally increased NOx emissions in both Korea (Fig. 4a) and China (Fig. 4b), seemingly compensating for the previously underrepresented morning peak NO2 concentrations (Figs. 2 and S1). Such a tendency towards increasing the emissions was less pronounced in China during the first hour of the day (7 AM local time), which is attributable to GEMS’s partial coverage that reaches only the eastern regions of China at that time. After the day progressed into the afternoon hours, we observed smaller extents of adjustments towards increasing the emissions in general, and even decreases in the emissions in Korea. This suggests that the a priori NOx emissions for this midday time window, not associated with the morning or evening peaks, might have been overestimated. Such adjustments towards reducing the emissions were also noted after the LEO-informed update, but throughout the entire daylight hours. The LEO-informed update led to uniform decreases and increases in the emissions in Korea and China, respectively, the directions of which were determined based on the discrepancy between the observed and modeled NO2 columns at 04:45 UTC (Fig. 1f) during the inversion process. These distinct responses of NOx emissions to the GEMS- and LEO-informed inversions suggest the potential of geostationary observation data in refining the emissions in a temporally more nuanced manner, as further discussed in the following section.
2.3 Simulations using the updated NOx emissions inventory
To assess the effectiveness of our top-down approach, we compared the model performances in simulating tropospheric NO2 columns and surface NO2 concentrations in Asia, achieved before and after updating the NOx emissions inventory. This enabled us to determine the extent to which the model’s earlier biases could be reduced by exploiting more contemporary satellite observations.
Figures 5c and 5d illustrate the monthly averages of hourly daytime tropospheric NO2 columns simulated after the GEMS-informed and LEO-informed updates. The GEMS-informed update generally remedied the model’s earlier underestimation (Figs. 5a and 5b), particularly in capturing the high peaks in East Asia, resulting in a closer alignment of the modeled columns with GEMS tropospheric NO2 columns. NO2 columns decreased after the update in areas where overestimation was once prevalent. Also, we noticed some instances of overcompensation. While effective in addressing underestimations of NO2 columns, the GEMS-informed update introduced overestimations in some areas, such as the Yellow Sea, where the columns were already well-captured using the a priori emissions. Given that the inversion was not directly applied to NOx emissions over sea surfaces, the increased NO2 columns over the Yellow Sea are attributed to significant increases in nearby upwind inland areas, such as the NCP and the northeastern region of China.
The LEO-informed update also addressed the model’s prior biases, but with different spatial patterns across the domain. The update was generally more effective in constraining previously overestimated NO2 columns but was less effective in capturing their high peaks in areas where underestimation was once prevalent, such as the southern half of the NCP. This was considered to be caused by the inversion, which was more substantially directed towards reducing the extent of NOx emissions, as discussed earlier in Section 2.2. These spatial patterns closely reflected the discrepancies observed between LEO proxy at 04:45 UTC and the corresponding modeled NO2 columns (Fig. 1c and 1f), reaffirming the ongoing challenge of preventing over- and under-correction in emissions inventories, especially when relying solely on time-averaged observation data11,36,52–55.
Figures 6 and 7 compare the hourly daytime surface NO2 concentrations observed and modeled in Korea and China, respectively, during MAM 2022. The GEMS-informed update improved R from 0.70 to 0.71 in Korea and from 0.73 to 0.78 in China on average during MAM (Table 1). To a lesser extent, the LEO-informed update also resulted in a slight improvement in correlation, maintaining R at 0.70 in Korea and advancing R to 0.74 in China. We also observed improvements in IOA, where the extents of improvement were more pronounced after the GEMS-informed update. The GEMS-informed update improved IOA from 0.78 to 0.81 in Korea and from 0.79 to 0.83 in China, while the LEO-informed update improved IOA to 0.79 in Korea and maintained IOA at 0.79 in China.
The model previously underestimated NO2 concentrations by an average of 19.23% in Korea during MAM, and the GEMS-informed update moderated this negative bias to an average of 11.36% (Table 1). The smallest bias was achieved in April, with NMB of 0.02%, indicating a fairly close agreement between the modeled and observed concentrations. Similarly, in China, the extent of underestimation was substantially reduced from an average of 12.85–4.83% (Table 1) during MAM, with exceptional alignments in April and May, with NMB of 0.78% and 0.79%, respectively. The LEO-informed update was also generally effective in reducing the model’s prior biases but to a lesser extent. The extents of underestimation were slightly reduced to an average of 16.96% in Korea and 9.93% in China. Despite aiming to counterbalance the model’s prior biases, the LEO-informed update led to less pronounced improvements. This is attributable to the limited amount of information available from observation references, which were more beneficial towards constraining the overestimated extent of NOx emissions as discussed in Section 2.2.
Table 1
Descriptive statistics comparing observed and modeled hourly surface NO2 concentrations averaged at 459 AirKorea sites in Korea and 250 MEE sites in China during MAM 2022. Prior: modeled concentrations using the a priori emissions, LEO: modeled concentrations after the LEO-informed update, GEMS: modeled concentrations after the GEMS-informed update. R: Pearson’s correlation coefficient, IOA: index of agreement, NMB: normalized mean bias (%), and RMSE: root mean square error (ppb).
| | Korea | China |
| | March | April | May | Average | March | April | May | Average |
R | Prior | 0.56 | 0.75 | 0.79 | 0.70 | 0.75 | 0.71 | 0.72 | 0.73 |
LEO | 0.58 | 0.74 | 0.77 | 0.70 | 0.76 | 0.72 | 0.74 | 0.74 |
GEMS | 0.55 | 0.77 | 0.80 | 0.71 | 0.79 | 0.77 | 0.77 | 0.78 |
IOA | Prior | 0.66 | 0.84 | 0.85 | 0.78 | 0.78 | 0.77 | 0.81 | 0.79 |
LEO | 0.68 | 0.84 | 0.85 | 0.79 | 0.80 | 0.76 | 0.80 | 0.79 |
GEMS | 0.68 | 0.87 | 0.88 | 0.81 | 0.84 | 0.81 | 0.84 | 0.83 |
NMB (%) | Prior | -29.31 | -13.17 | -15.22 | -19.23 | -23.58 | -9.88 | -5.08 | -12.85 |
LEO | -27.75 | -9.74 | -13.38 | -16.96 | -19.28 | -6.98 | -3.52 | -9.93 |
GEMS | -24.84 | 0.02 | -9.27 | -11.36 | -14.82 | 0.78 | 0.79 | -4.42 |
MAE (ppb) | Prior | 5.90 | 2.80 | 2.19 | 3.63 | 4.63 | 4.09 | 3.01 | 3.91 |
LEO | 5.74 | 2.77 | 2.11 | 3.54 | 4.36 | 3.90 | 3.27 | 3.84 |
GEMS | 5.39 | 2.68 | 1.89 | 3.32 | 3.70 | 3.15 | 2.69 | 3.18 |
Figure 8 illustrates the daytime mean surface NO2 concentrations observed and modeled in Korea and China during MAM, aligned with the number of valid observations afforded by GEMS retrievals. While the GEMS-informed update was generally more effective in moderating the model’s prior biases, a noticeable aspect was the response of the a posteriori concentrations to the availability of valid observation references. On days with minimal valid data given (between 0 and 1 record on average), including March 13, 14, 17, 18, and April 13 in Korea, the GEMS-informed inversion led to minor adjustments in the corresponding NO2 concentrations due to limited top-down information (Fig. 8a). Meanwhile, the LEO-informed inversion allowed the modeled concentrations to be adjusted in a continuous manner, taking advantage of the monthly-averaged observation data applied during the inversion, but with less pronounced improvements. A similar pattern was observed sporadically between the modeled and observed concentrations in China, but to a less severe extent (Fig. 8b).