Air quality data were obtained from the ChinaHighAirPollutants (CHAP) database, developed by our team and extensively applied in environmental health studies. The CHAP database integrates ground observations, satellite retrievals, and chemical transport model simulations, addressing gaps in satellite data while capturing nonlinear relationships between satellite column retrievals and surface pollutant concentrations through advanced spatiotemporal artificial intelligence (AI) models. This integration yields a long-term, gap-free dataset with complete spatial coverage and high spatiotemporal resolution. In this study, six conventional air pollutants from 2019 to 2022 were analyzed, including two particulate species (PM_2.5 and PM₁₀) and four gaseous pollutants (i.e., O₃, NO₂, SO₂, and CO), all uniformly resolved at 1 km spatial resolution on a daily basis.

These AQ products are derived from the latest version, featuring substantial enhancements to key input variables within a unified AI framework – the four-dimensional spatiotemporal deep forest (4D-STDF) model – which aims to minimize inconsistencies across datasets and to ensure higher accuracy and reliability. A key improvement over the original deep forest model is the incorporation of spatiotemporal information. Spatial variability is represented using three spherical coordinate features derived from longitude (Lon) and latitude (Lat), i.e., Ps∼[S1,S2,S3]=[sin⁡2πLon360, cos2πLon360sin2πLat180, cos2πLon360cos2πLat180], which capture spatial autocorrelation patterns. Temporal variability is represented using three helix-shaped trigonometric features based on the normalized day of the year (DOYN), i.e., Pt∼[T1,T2,T3]=[DOYN, cos2πDOYN, sin2πDOYN], which capture both linear temporal progression and periodic (seasonal) variations. Together, these features enable the model to better represent daily variations and seasonal cycles in air pollution.

Specifically, PM_2.5 and PM₁₀ were estimated using MODIS MAIAC 1 km AOD products as the primary predictors (Wei et al., 2021), whereas surface concentrations of O₃ (Yang et al., 2025), NO₂ (Wei et al., 2022b), SO₂ (Wei et al., 2023c), and CO (Wei et al., 2023c) were derived from downscaled 1 km TROPOMI vertical column density products utilizing a geostatistical technique (Cersosimo et al., 2020) implemented on the Google Earth Engine platform. Table S1 in the Supplement provides a summary of all data utilized for pollutant modeling in this study.

During the modeling process, ground-based measurements of each air pollutant served as the target variables, while predictors were specifically tailored for each pollutant. For PM_2.5, key input variables included MODIS MAIAC gap-filled AOD, GEOS-CF PM_2.5 simulations, and anthropogenic emissions (e.g., PM_2.5, NH₄, NO_x, SO₂, CO, and VOCs) from the daily 1 km anthropogenic emissions provided by Tsinghua's Air Benefit and Cost and Attainment Assessment System–Emission Inventory (ABaCAS-EI; Li et al., 2023). Additional auxiliary variables influencing PM_2.5 concentrations included ERA5 meteorological factors such as boundary-layer height (BLH), air temperature (TEM), relative humidity (RH), surface pressure (SP), total precipitation (TP), total evaporation (ET), wind speed (WS), and wind direction (WD). Further predictors comprised the normalized difference vegetation index (NDVI), digital elevation model (DEM), population density (POD), nighttime lights (NTL), and spatiotemporal terms (Ps and Pt): 1PM2.5∼f4D-STDF(MAIACAOD,GEOS-CFPM2.5,ABaCASPM2.5,NH4,NOx,SO2,CO,VOCs,ERA5BLH,TEM,SP,RH,TP,ET,WS,WD,NDVI,DEM,POD,NTL,PsPt) For PM₁₀, the input variables included MODIS MAIAC gap-filled AOD, CAMS PM₁₀ simulations, and ABaCAS-EI PM₁₀ anthropogenic emissions, along with the same meteorological, surface-related, population-related, and spatiotemporal variables used for PM_2.5, as expressed below: 2PM10∼f4D-STDF(MAIACAOD,CAMSPM10,ABaCASPM10,ERA5BLH,TEM,SP,RH,TP,ET,WS,WD,NDVI,DEM,POD,NTL,PsPt) For surface O₃, the main predictors included the TROPOMI gap-filled tropospheric NO₂ and HCHO columns, GEOS-CF surface O₃ simulations, and ABaCAS-EI anthropogenic emissions of NO_x, VOCs, and CO. Additionally, two key factors influencing photochemical reactions for ozone formation – satellite-derived surface downward shortwave radiation (DSR) and land surface temperature (LST) products – were integrated. The model also included the same auxiliary variables used for both PM_2.5 and PM₁₀: 3O3∼f4D-STDF(TROPOMINO2,HCHO,GEOS-CFO3,ABaCASNOx,VOCs,CO,DSR,LST,ERA5BLH,SP,RH,TP,ET,WS,WDNDVI,DEM,POD,NTL,PsPt) For other surface gaseous pollutants (i.e., NO₂, SO₂, and CO), the input variables included the TROPOMI gap-filled tropospheric or total column, GEOS-CF surface gaseous pollutant simulations, and ABaCAS-EI anthropogenic emissions, along with the same meteorological, surface-, population-related, and spatiotemporal variables used for PM_2.5 and O₃: 4NO2∼f4D-STDF(TROPOMINO2,GEOS-CFNO2,ABaCASNOx,ERA5BLH,TEM,SP,RH,TP,ET,WS,WD,NDVI,DEM,POD,NTL,PsPt)5SO2∼f4D-STDF(TROPOMISO2,GEOS-CFSO2,ABaCASSO2,ERA5BLH,TEM,SP,RH,TP,ET,WS,WD,NDVI,DEM,POD,NTL,PsPt)6CO∼f4D-STDF(TROPOMICO,GEOS-CFCO,ABaCASCO,ERA5BLH,TEM,SP,RH,TP,ET,WS,WD,NDVI,DEM,POD,NTL,PsPt) These estimated daily AQ data are comprehensively validated against surface measurements from over 1700 ground-based monitoring stations across China, as reported by the China National Environmental Monitoring Center (MEE, 2018). The evaluation utilizes three commonly used independent validation approaches for AI regression tasks: out-of-sample, out-of-station, and out-of-day ten-fold cross-validation (10-CV). These techniques are used to assess overall accuracy, spatial predictive accuracy, and temporal predictive accuracy, respectively. The latter two methods are specifically designed to evaluate the AI model's performance in predicting AQ data for areas and days without ground observations (Wei et al., 2023a, b, d).

Air Quality Index (AQI) is used to communicate daily air quality and is categorized into six levels, indicating different degrees of public exposure to ambient air pollution: Good (0≤AQI≤50), Moderate (50<AQI≤100), Unhealthy for Sensitive Groups (100<AQI≤150), Unhealthy (150<AQI≤200), Very Unhealthy (200<AQI≤300), and Hazardous (AQI >300). This index is established by calculating the Individual Air Quality Index (IAQI) for six major air pollutants: fine and coarse particulate matter (PM_2.5 and PM₁₀), and gaseous pollutants, including ozone (O₃), nitrogen dioxide (NO₂), sulfur dioxide (SO₂), and carbon monoxide (CO): 7IAQIp=Cp-Cp,l×IAQIp,j-IAQIp,j-1Cp,h-Cp,l+IAQIp,j-1 where IAQI_p represents the IAQI of pollutant p; Cp is the surface concentration of pollutant p. Specifically, we utilize raw floating-point pollutant concentrations (Cp) throughout the conversion process without any rounding, truncation, or exclusion. The thresholds Cp,1 and Cp,h are defined as the lower and upper surface concentration limits for pollutant p, respectively, selected based on proximity to Cp, while IAQI_p,j and IAQI_p,j−1 correspond to IAQI values for these respective limits (Table S2). All intermediate values are retained in floating-point precision to avoid discretization bias, and IAQI values are capped at 500.

In this study, we employed the refined Aggregated Air Quality Index (AAQI), which accounts for the combined effects of multiple criteria pollutants (Kyrkilis et al., 2007), offering a more comprehensive assessment of air quality by integrating various sources of pollutant exposure: 8AAQI=∑i=1nIAQIpρIρ where ρ is a constant that modulates the aggregation effect of individual IAQIs. Here, we set ρ to 2, allowing us to consider the combined effects of each pollutant while still emphasizing the impact of certain more severe pollutants in China (Hu et al., 2015). Similarly, the AAQI values are also capped at 500 to avoid skewing the assessment.

To comprehensively evaluate the impact of abrupt emission changes on air quality, the study period was divided into three distinct eras: the pre-reduction period (2019), the reduction period (2020), and the recovery period (2021–2022). First, air quality (AQ) data were adjusted to minimize the influence of short-term meteorological variations using our previously developed method. Next, Seasonal-Trend Decomposition (Cleveland et al., 1990) was applied to separate short-term fluctuations and seasonal effects, ensuring that observed trends reflect genuine changes in pollutant levels: 9yt=Tt+St+Rt where yt is the observed time series, Tt is the trend component, St is the seasonal component, and Rt is the residual (noise) component. The Seasonal and Trend components are extracted using the Locally estimated scatterplot smoothing (STL) model, a nonparametric regression method that fits a smooth curve to the data with local polynomial regression within a moving window and assigns weights based on the distance from the window's center: 10min⁡∑iωi(yi-(b0+b1xi-x+…))2 where ωi is the weight determined by a kernel function (e.g., tricubic kernel) for the smoothed value at a point x using the weighted least squares method.

The response duration – defined as the time required for pollutants to return to baseline levels after abrupt emission changes – was determined by identifying the start and end dates of emission perturbation phases based on intersections of pollutant levels between pre-reduction and reduction periods. Using this response duration, changes in air pollution conditions across different emission phases were quantified. To establish a common temporal baseline across years, the Chinese New Year was set as the start of the year, mitigating holiday-related impacts on pollutant levels.

Considering that certain factors (e.g., land use and surface elevation) remain relatively stable over short timescales, the drivers of air pollution conditions were categorized into meteorological and anthropogenic factors, with a focus on two major secondary pollutants. For PM_2.5, anthropogenic factors included surface concentrations of SO₂, NO₂, and CO from the CHAP daily AQ datasets, as well as surface VOCs from the CAMS reanalysis. For surface O₃, major precursors included CHAP daily surface NO₂ and CO concentrations together with CAMS surface VOC concentrations. Meteorological variables were sourced from the ERA5 reanalysis dataset (Hersbach et al., 2020; Muñoz-Sabater et al., 2021) and included boundary-layer height (BLH), relative humidity (RH), surface pressure (SP), total precipitation (TP), air temperature (TEM), and wind speed and direction for PM_2.5. For O₃, additional variables related to photochemical reactions were incorporated, including MODIS daily surface downward shortwave radiation (DSR) products (Yang et al., 2025) and ERA5 daily surface sensible heat flux and total cloud cover. Specifically, auxiliary datasets with coarser spatial resolutions were resampled to a common 1 km × 1 km grid using bilinear interpolation to match the CHAP pollutant data.

The population-weighted concentration (PWr,p) for pollutant p in region r is defined as: 11PWr,p=∑i∈cCi,p×Pi∑i∈cPi where Ci,p is the pollutant concentration Pi is the population in grid cell i. Population data are derived from the LandScan™ dataset and adjusted to match United Nations annual total population estimates.

To reveal the drivers of air pollution level changes under different emission phases, we employed SHapley Additive exPlanations (SHAP; Antwarg et al., 2021), an advanced XAI (eXplainable Artificial Intelligence) method. SHAP provides a flexible and interpretable framework for ML/DL modeling, especially when spatial and non-spatial effects interact (Cheng et al., 2025; Li, 2022; Wei et al., 2024). Key features of SHAP include consistent allocation of feature contributions and the ability to capture complex interactions while maintaining global and local interpretability. The contribution of each factor can be expressed as: 12ϕf=1N∑k=1N|∑S⊆NS!N-S-1!N![f(S∪i)-f(S))]| where ϕf represents the SHAP value for a given feature, quantifying its average contribution to the model's prediction. N is the total number of samples, and S is a subset of features with |S| as its size, and |N| as the total number of features. The weight factor S!N-S-1!N! measures the importance of the subset S within all possible feature permutations. f(S∪i) denotes the model output when feature i is included, while f(S) represents the output when it is excluded.

In addition to SHAP, the traditional multiple linear regression (MLR) method was performed for comparison, using standardized regression coefficients to evaluate the average contribution of each factor. The same set of input variables used in the XAI analysis was applied, along with additional feature testing to address multicollinearity issues (Tables S3 and S4), evaluated using the Variance Inflation Factor (VIF) method (Wei et al., 2022a) across national, regional, and urban scales. A threshold of VIF <10 was adopted. Although some variables showed relatively high VIF values (e.g., 9.86 for air temperature and 8.79 for surface pressure in the national O₃ model), all remained below this threshold. Therefore, no variables were excluded, and multicollinearity is unlikely to substantially affect the results.

Our developed uniform DL model can effectively estimate daily ambient particulate and gaseous pollutants, demonstrating high sample-based cross-validation coefficients of determination (CV-R2) [low root mean-square errors (RMSE)] values – 0.92 (16.87 µgm-3) for PM₁₀, 0.93 (7.58 µgm-3) for PM_2.5, 0.94 (11.48 µgm-3) for O₃, 0.87 (6.34 µgm-3) for NO₂, 0.79 (3.93 µgm-3) for SO₂, and 0.80 (0.18 mgm-3) for CO – validated against over 2.2 million ground-based observations collected from 2019 to 2022 at the national scale (Fig. 1). Strong model performance is also observed at individual ground-based monitoring stations across the domain (e.g., 52 %–99 % exhibiting moderate to strong CV-R2 values >0.5), particularly in the eastern regions (Fig. S1 in the Supplement). Additionally, the model maintains consistent reliability in both spatial and temporal predictions of various air pollutant species across China, with minor declines in station-based and day-based cross-validated R² values (ranging from 0.63 to 0.92, and 0.64 to 0.8, respectively) when randomly dropping ground-based monitoring stations and observation days (Figs. S2–S3). These results underscore the strong potential of deep learning in monitoring ambient air pollution from space, particularly in areas and on days lacking direct ground measurements. However, the performance of specific trace gas retrievals (e.g., SO₂ and CO) is relatively poorer than that of particulate matter pollutants, primarily due to the weaker signals and reduced sensitivity of satellite observations.

Furthermore, year-specific validation demonstrates stable model performance over time, showing consistent sample-based CV-R2 (RMSE) values during 2019–2022, ranging from 0.89–0.93 (12.51–20.21 µgm-3) for PM₁₀, 0.93–0.94 (6.40–8.79 µgm-3) for PM_2.5, 0.93–0.95 (10.11–13.59 µgm-3) for O₃, 0.86–0.87 (5.80–7.19 µgm-3) for NO₂, 0.76–0.81 (3.13–4.87 µgm-3) for SO₂, and 0.75–0.82 (0.16–0.20 mgm-3) for CO, indicating consistent reliability despite interannual variations in emissions and meteorology. In addition, the number of monitoring stations shows a moderate increase over time, ranging from 1474 to 1896 (Fig. S4), while the spatial distribution remains broadly stable, ensuring that interannual model stability is not confounded by changes in sampling density or spatial coverage.

At the regional scale (Fig. S5), both pollutants show strong performance, with sample-based CV-R2 (slope) exceeding 0.93 (0.90) in both Northern and Southern China, and RMSE (MAE) ranges from 5.45–9.42 (3.68–5.26) µgm-3 for PM_2.5 and 10.84–11.78 (7.33–8.18) µgm-3 for O₃. At the city scale, Wuhan also demonstrates high accuracy, with CV-R2 of 0.91 (slope = 0.91) for PM_2.5 and 0.94 (0.91) for O₃, along with average RMSE (MAE) values of 7.12 (5.05) and 12.71 (8.95) µgm-3, respectively. These results confirm the model's robustness across spatial scales, including its applicability to fine-scale urban analyses.

Using machine learning, we derived daily, gapless 1 km maps of six ambient air pollutants across China from satellite observations, providing spatially continuous estimates at each location and for each day, particularly in regions lacking ground-based measurements. Figure 2 presents a representative example on 1 January 2020, illustrating pronounced pollutant-specific spatial heterogeneity. Particulate matter exhibits strong regional clustering, with PM₁₀ showing elevated concentrations in arid and semi-arid regions of northwestern China, largely influenced by natural dust emissions and surface wind erosion, while PM_2.5 displays prominent hotspots over major urbanized and industrialized regions such as the North China Plain, the Fenwei Plain, and the Sichuan Basin. These PM_2.5 enhancements are closely associated with intensive anthropogenic emissions and unfavorable meteorological conditions during winter, including shallow boundary layers and weak ventilation, which promote pollutant accumulation.

In contrast, gaseous pollutants exhibit distinct spatial patterns that reflect differences in emission sources and chemical lifetimes. NO₂ exhibits strong urban-scale gradients and is tightly linked to traffic and industrial combustion sources, whereas SO₂ concentrations are generally low nationwide due to effective emission control policies, with residual enhancements confined to localized industrial regions. CO presents a relatively homogeneous spatial distribution owing to its longer atmospheric lifetime and regional transport. O₃ shows higher concentrations in western and high-altitude regions, while remaining lower in eastern urban areas during winter, reflecting weaker photochemical production and enhanced NO titration. Together, these results highlight the ability of high-resolution, gapless datasets to capture multi-scale characteristics of air pollution and their underlying drivers across China.

These high-resolution AQ datasets enable a detailed examination of short-term pollutant responses to abrupt emission changes (Fig. 3). During the period of stringent emission reductions, marked by elevated values of the Oxford Coronavirus Government Response Tracker (OxCGRT) index in 2020 (Hale et al., 2021), most pollutants exhibited substantial declines, although the magnitude and temporal response varied markedly among species. The Individual Air Quality Index (IAQI) for PM_2.5 showed a sharp decline of 30 % over 22 d, while PM₁₀ decreased by 15 %, taking twice as long – 47 d – to return to normal levels. NO₂ experienced the largest drop, dropping by 32 % over a 52 d response period, while SO₂ and CO demonstrated moderate reductions of 21 % and 19 % over 45 and 25 d, respectively, driven by substantial cuts in anthropogenic emissions from transportation and industrial sources (Wei et al., 2023a, 2022b) In contrast, O₃ concentrations increased by 18 % and exhibited the fastest rebound time (17 d), largely attributable to the substantial reduction in nitrogen oxides (NO_x) emissions under a VOC-limited (NO_x-saturated) chemical regime across much of eastern China (Liu et al., 2021; Wei et al., 2022a). Overall, the Average Air Quality Index (AAQI) declined by 17 % during the most stringent reduction period, indicating widespread improvements in air pollution conditions across China.

During the subsequent emission recovery phase, pollutant-specific responses varied considerably. Approximately half of the pollutant species showed rebounds consistent with a partial resumption of anthropogenic emissions. NO₂ increased most notably, rising by 25 % relative to the reduction phase, reflecting the recovery of traffic and industrial activities. PM₁₀ increased by 16 %, influenced by spring dust events and renewed human activities. O₃ exhibited a slight increase of 1 %, highlighting the continued complexity of ozone formation driven by changes in precursor emissions, particularly NO_x. In contrast, PM_2.5 decreased marginally by 1 %, suggesting more persistent improvements, while SO₂ and CO decreased by 3 % and 8 %, respectively, indicating sustained control of primary emissions. Overall, the AAQI increased by 10 % during the recovery phase, reflecting a partial decline in air pollution conditions; however, conditions remained better than pre-reduction levels, with the aggregated indicator still 9 % lower on average, particularly for PM_2.5 (31 % reduction) and its major precursors (16 %–26 % reductions). These results suggest that although some benefits diminished as emissions resumed, structural reductions and ongoing policy measures contributed to lasting improvements in air quality.

During the period of stringent emission reductions, Wuhan, a major urban center with intensive anthropogenic activities, experienced more pronounced and rapid air quality changes than the national average (Fig. 4). PM_2.5 and PM₁₀ levels declined by 37 % and 34 %, respectively, with shorter response times of 35 and 38 d, indicating a stronger and faster particulate matter response to abrupt emission decreases. In contrast, O₃ levels increased sharply by 73 %, about four times the national average, reflecting an intensified ozone formation response under substantially reduced NO_x emissions. NO₂ exhibited the largest decrease (47 %) and the longest rebound period (65 d), consistent with the strong suppression of traffic and industrial emissions in the city. SO₂ and CO experienced moderate declines of 28 % and 12 %, respectively, with response times 2.3–2.4 times shorter than national averages. Overall, the AAQI in Wuhan decreased by 26 % over an average recovery period of 61 d. Notably, the end of this recovery period, identified from satellite observations, closely matched the timing of the city's emission control easing, differing by only ∼5 d, demonstrating the capability of space-based measurements to capture city-scale responses to emission interventions.

During the subsequent emission recovery phase, Wuhan's ambient air pollution levels partially rebounded, with NO₂ and PM₁₀ increasing by 58 % and 35 %, respectively, approximately 2.2–2.3 times the corresponding national averages. PM_2.5 levels also rose by 24 %, contrasting with the slight national decline. In comparison, O₃ exhibited a modest increase of 6 %, while SO₂ and CO continued to decrease at rates similar to national trends. Consequently, the composite index increased by 21 %, reflecting the partial resumption of urban activities. Compared with national patterns, Wuhan experienced more pronounced and faster reductions in particulate matter and NO₂ during the emission reduction phase, a substantially larger O₃ enhancement, and a stronger rebound in NO₂ during the recovery phase, highlighting the effects of stricter local emission controls followed by a rapid return of anthropogenic emissions. Despite this rebound, overall air pollution conditions in Wuhan remained improved relative to pre-intervention levels: except for O₃, which increased by 69 %, all other pollutants decreased by 10 %–21 %, resulting in an overall AAQI reduction of 11 %. These results indicate that, although emissions partially recovered, lasting air quality improvements persisted following the targeted emission control measures.

During the most stringent emission reduction period, only 13 % and 9 % of inhabited areas (population density >1 person km⁻²) and counties experienced an increase in PM₁₀ compared to pre-reduction levels (Figs. 5a and S6a), primarily in the desert-dominated northwestern regions affected by enhanced dust events. In contrast, PM_2.5 decreased in 74 % (76 %) of inhabited areas (counties), particularly across the densely populated Northern China Plain, with reductions exceeding 50 %. Pronounced declines in NO₂, SO₂, and CO were observed in 92 % (98 %), 85 % (90 %), and 72 % (87 %) of inhabited areas (counties), respectively, with the largest reductions occurring in urbanized regions of eastern China, especially Wuhan city, reflecting substantial suppression of direct emissions from transportation and industry. O₃ levels, conversely, increased in 76 % (84 %) of inhabited areas (counties), consistent with the spatial reduction of NO₂, which promotes ozone formation under a VOC-limited chemical regime. Overall, atmospheric conditions improved markedly, with 73 % of inhabited areas and 87 % of counties showing decreases in the Aggregated Air Quality Index (AAQI).

During the recovery period, fine-scale air quality patterns exhibited pronounced heterogeneity compared to the reduction period (Figs. 5b and S6b). PM₁₀ increased in 87 % (90 %) of inhabited areas (counties), particularly in Northern and Southwestern China, while PM_2.5 worsened in approximately half of the domain (49 % and 52 % of inhabited areas and counties), mainly in urbanized southeastern regions, including the Yangtze and Pearl River Deltas. NO₂ experienced the largest increase, affecting 90 % of inhabited areas and 96 % of counties, including rural and less industrialized regions. This trend was particularly evident in Hubei province and surrounding areas, where increases exceeded 50 %, reflecting the resumption of construction, transportation, and industrial activity. Moderate increases in SO₂ and CO (47 % and 18 % of inhabited areas, and 48 % and 23 % of counties) were primarily observed in specific industrial regions, suggesting a rebound in coal-dependent industries and partial relaxation of emission controls. O₃ rose in 55 % (48 %) of inhabited areas (counties), particularly in southern regions, likely due to high solar radiation and elevated precursor concentrations. While western regions generally experienced continued air quality improvements, central and eastern China saw noticeable declines. Overall, 72 % of inhabited areas and 70 % of counties showed composite index increases, highlighting the impact of emission recovery on ambient air pollution levels.

Comparing the recovery period to pre-reduction levels (Figs. 5c and S6c), PM₁₀ worsened in 49 % (43 %) of the domain (counties), particularly in central-northern and southern regions, while PM_2.5 IAQI increased in 26 % (24 %) of areas (counties), mainly in densely populated and industrialized southern regions, reflecting the influence of resumed human activities. NO₂ increased in 32 % (17 %) of inhabited areas (counties), while SO₂ and CO showed smaller increases (16 % (8 %) and 13 % (5 %), respectively), indicating that emission controls on vehicles and certain industrial sectors remained effective in many areas. In contrast, O₃ rose sharply in 87 % and 92 % of inhabited areas and counties, particularly in southern and southeastern China, suggesting that elevated precursor concentrations (NO_x and VOCs) and favorable meteorological conditions enhanced ozone formation. Overall, air quality improved in 53 % of inhabited areas and 66 % of counties relative to pre-reduction levels, especially in densely populated central and eastern regions, although other regions experienced localized declines in air pollutant concentrations.

Leveraging the daily, gapless air pollutant dataset, we assessed changes in the number of days at different AAQI levels during the first quarter (January–March), which typically experiences severe pollution in China. For example, 40 % (55 %), 34 % (44 %), and 33 % (48 %) of inhabited areas (counties) experienced unhealthy air quality (AAQI >100) for half of the quarter during the pre-reduction, reduction, and recovery periods, respectively (Figs. 6 and S7). Central and eastern China were particularly affected, with pollution exposure frequency exceeding three-fourths of the quarter. Compared to the pre-reduction period (Fig. S8), “Good” days (AAQI <50) increased in 12 % and 38 % of inhabited areas and counties, primarily in southern China, while “Moderate” days (AAQI = 51–100) increased in 72 % and 86 % of areas, mostly in central and eastern regions. More notably, “Unhealthy for Sensitive Groups” (AAQI = 101–150) days decreased in 66 % (62 %) of inhabited areas (counties), especially in central and southeastern China, while reductions in “Unhealthy” (AAQI = 151–200) and “Very Unhealthy” days reached 85 % (86 %) and 93 % (91 %), with hotspots in the North China Plain. Overall, healthy (AAQI <100) days increased in 76 % and 87 % of inhabited areas and counties (Figs. 6d and S7d), reflecting widespread improvement in ambient air pollution conditions during the emission reduction period.

During the recovery period, “Moderate” days showed a mixed pattern, with increases in northern China but decreases in southern and eastern regions (amplitude >20 %), resulting in an overall increase in 65 % of inhabited areas and 58 % of counties (Fig. S9). In contrast, higher pollution levels – “Unhealthy for Sensitive Groups” and “Unhealthy” days – substantially increased in 59 % (68 %) and 50 % (72 %) of areas, mainly in central and southern China, consistent with the partial resumption of anthropogenic emissions. “Good” and “Very Unhealthy” categories showed minimal changes, whereas “Hazardous” days (AAQI = 301–500) increased in 57 % (50 %) of inhabited areas (counties). These patterns indicate an overall deterioration in air pollution conditions during emission recovery (Figs. 6e, S7e), with unhealthy days (AAQI >100) rising in 63 % and 71 % of inhabited areas and counties.

Comparing the recovery to the pre-reduction period (Fig. S10), “Good” days remained similar to the reduction period, with minimal changes in southern regions. “Moderate” days increased in 63 % and 73 % of inhabited areas and counties, with widespread growth across central and southern China, especially in Hubei and surrounding areas, while declines were observed in southern provinces such as Guangdong and Guangxi. Continuous reductions were observed in “Unhealthy for Sensitive Groups,” “Unhealthy,” and “Very Unhealthy” days across 52 % (43 %), 73 % (68 %), and 82 % (79 %) of areas (counties), particularly in the Shandong Peninsula and nearby regions. Meanwhile, “Hazardous” days showed minor spatial changes (<5 %). These results suggest that, although emission levels partially recovered, overall air quality has generally improved relative to pre-reduction levels (Figs. 6f, S7f), with a clear north–south contrast – increasing healthy days in northern China but decreasing in southern regions (65 % vs. 35 % of inhabited areas, 73 % vs. 27 % of counties).

Using eXplainable Artificial Intelligence (XAI), we examined the drivers of changes in PM_2.5 and O₃ concentrations across the pre-reduction (2019), reduction (2020), and recovery (2021–2022) periods. For PM_2.5, population-weighted concentrations declined by 32 % during the reduction period, primarily driven by a 22 % decrease in anthropogenic emissions (Fig. 7a). Sharp emission reductions from industrial activities and transportation led to marked decreases in particulate matter and its precursors. Meteorological conditions contributed a secondary yet important role, accounting for an additional 10 % reduction through enhanced dispersion and increased precipitation. During the recovery period, PM_2.5 stabilized at 40.5 µgm-3, with minor changes attributed to meteorology (1 % reduction) and a modest rebound in anthropogenic emissions (4 % increase). These results highlight the dominant role of human-driven emissions in controlling PM_2.5 levels, even as air quality improvements persisted post-reduction.

For O₃, the response differed: during the reduction period, population-weighted concentrations increased by 18 % relative to pre-reduction levels, largely due to meteorological factors (15 %) and partially to emission changes (3 %) (Fig. 7b). This increase reflects the complex photochemical response to reduced NO_x emissions during lockdowns: Lower NO_x emissions during the reduction period decreased O₃ titration in urban areas, resulting in higher ambient O₃ through nonlinear photochemical processes. During the recovery period, O₃ remained elevated at 84.6 µgm-3, with negligible additional contributions from meteorology or anthropogenic emissions. These findings underscore the dominant influence of meteorology on O₃ variations, modulated by chemical responses to emission reductions.

Considering regional differences north and south of the Qinling-Huaihe Line, we conducted separate XAI analyses (Fig. 7c–f). During the reduction period, PM_2.5 changes mirrored national trends, dominated by anthropogenic emission decreases (23 % in northern China and 11 % in southern China), while O₃ variations were primarily meteorology-driven (12 % and 13 %, respectively). In the recovery period, PM_2.5 and O₃ continued to decline in northern China by 5 % and 6 %, with anthropogenic factors dominating, whereas both pollutants increased in southern China. In this region, PM_2.5 increases were driven by both human activities and meteorology, while O₃ rises were largely meteorology-driven, influenced by high temperature and humidity.

At the city scale, Wuhan exhibited distinctive trends (Fig. 7g–h). During the reduction period, PM_2.5 decreased by 37 %, mainly due to human emission reductions (22 %), whereas O₃ increased by 40 %, dominated by meteorological factors. During the recovery period, PM_2.5 rebounded by 26 %, largely driven by meteorology (∼16 %), while O₃ decreased by 5 %, with meteorology again as the primary driver. These XAI results are consistent with traditional Multiple Linear Regression (MLR) analyses (Fig. S11) but provide a more nuanced understanding of nonlinear interactions between pollutants and their drivers at national, regional, and urban scales.

Overall, the results highlight the contrasting responses of O₃ and PM_2.5 to changes in emissions: reductions in primary pollutants yield clear benefits for particulate matter but can lead to increased O₃ through complex photochemistry. This underscores the need for balanced emission control strategies that simultaneously address primary and secondary pollutants, especially in densely populated and industrialized regions. The analysis further shows that while emission reductions can rapidly improve air quality, adaptive regulatory frameworks are essential for sustainable management of both PM_2.5 and O₃.