Meteorology impact on PM 2 . 5 change over a receptor region in the regional transport of air pollutants : observational study of recent emission reductions in central China

As an important issue in atmospheric environment, the contributions of anthropogenic emissions and meteorological conditions to air pollution have been little assessed over receptor regions in regional transport of air pollutants. In the present study of 5-year observations and modeling, we targeted the Twain-Hu Basin (THB), a large region of heavy PM2.5 pollution in central China, to assess the effects of meteorology on PM2.5 change over a receptor region in the regional transport of air pollutants. Based on observations of environment and meteorology over 2015–2019, the Kolmogorov–Zurbenko (KZ) filter was performed to decompose the PM2.5 variations into multiple timescale components over the THB, where the short-term, seasonal and long-term components accounted for, respectively, 47.5 %, 41.4 % and 3.7 % of daily PM2.5 changes. The short-term and seasonal components dominated the day-to-day PM2.5 variations with long-term component determining the change trend of PM2.5 concentrations over recent years. As the emissionand meteorology-related long-term PM2.5 components over the THB were identified, the meteorological contribution to the declining PM2.5 trend presented a distinct spatial pattern over the THB with northern positive rates up to 61.92 % and southern negative rates down to −24.93 %. The opposite effects of meteorology on PM2.5 pollution could accelerate and offset the effects of emission reductions in the northern and southern THB, which is attributed to the upwind diffusion and downward accumulation of air pollutants over receptor regions in regional PM2.5 transport. It is noteworthy that the increasing conversion efficiencies of SO2 and NO2 to sulfate and nitrate for secondary PM2.5 could have offset the effect of PM2.5 emission reduction on air pollution in the THB during recent years, revealing the enhancing contribution of gaseous precursor emissions to PM2.5 concentrations under control of anthropogenic emissions of PM2.5 and the gaseous precursors over receptor regions in the regional transport of air pollutants. Our results highlight the effects of emission mitigation and meteorological changes on the source–receptor relationship of the regional transport of air pollutants with the implication of long-range transport of air pollutants for regional and global environment changes. Published by Copernicus Publications on behalf of the European Geosciences Union. 3580 X. Sun et al.: Meteorology impact on PM2.5 change over a receptor region


Introduction
Haze pollution with high levels of PM 2.5 (fine particulate matters with aerodynamic diameters equal to or less than 2.5 µm) has been a serious problem in the atmospheric environment (Peng et al., 2016;Wang et al., 2016) with adverse influences on air quality and human health (Cao et al., 2012;Crouse et al., 2012). In recent years, large areas in central and eastern China (CEC) have experienced haze pollution with unprecedentedly high PM 2.5 levels in the regions covering the North China Plain (NCP), Yangtze River Delta (YRD), Pearl River Delta (PRD) and Sichuan Basin (SB) (Zhang et al., 2012;Guo et al., 2017). In order to improve air quality by reducing air pollutant emissions, the Chinese government has implemented an action plan for controlling anthropogenic emissions since September 2013 (http://www.gov. cn/xinwen/2018-02/01/content_5262720.htm, last access: 17 January 2022). Surface PM 2.5 concentrations have exhibited 30 %-40 % decreases in CEC over recent years (Xue et al., 2019;Zhang et al., 2019). However, the changes of air pollution are generally co-determined by air pollutant emissions and meteorological conditions. The contributions of changes in meteorology and anthropogenic emissions to the improvement of air quality need to be comprehensively investigated. PM 2.5 includes primary particles emitted directly from various sources and secondary particles generated by homogeneous and heterogeneous chemical reactions of gaseous precursor in the atmosphere, depending on the emissions of the primary PM 2.5 and its gaseous precursors (Y. Du et al., 2020). In addition to the emissions of air pollutants, meteorological conditions can alter the local accumulation, regional transport, chemical conversion and wet and dry depositions of air pollutants (Lu et al., 2017;. Severe haze pollution always occurs in the wintertime under stagnant meteorological conditions with weak near-surface winds, strong temperature inversion and high relative humidity in the atmospheric boundary layer, which are favorable for the accumulation of air pollutants to form air pollution Miao et al., 2015;Tang et al., 2016). Meteorological conditions are closely governed by synoptic circulation by modulating the atmospheric physical and chemical processes, including regional transport of air pollutants (Miao et al., 2017;Ning et al., 2019). By regulating synoptic circulation patterns, the climate changes of East Asian monsoons greatly influence the seasonal and interannual variations of aerosol concentrations for air pollution over China (Zhu et al., 2012;Jeong and Park, 2017).
Assessments of contributions of anthropogenic emissions and meteorological changes to air quality improvement are an important issue in environmental changes (Pearce et al., 2011;Zhang et al., 2018;Chen et al., 2019). Chemical transport models have been widely used to quantify meteorology's effect on PM 2.5 variations by using a linear additive relationship between sensitivity and base simulations (Mueller and Mallard, 2011;Zhang et al., 2020). The contribution of meteorological changes to PM 2.5 decreases was estimated at the averages of 10 %-20 % with interannual fluctuations of about 5 % in CEC from 2015 to 2019 through a model-based environmental meteorology index (Gong et al., 2021). The accuracy of modeling assessments can be influenced by the uncertainties in emission inventories and the incomplete chemical and physical mechanisms in air pollution simulation (Li et al., 2011). Based on the statistical analysis of the long-term observational data, it was quantified that the emission controls could explain more of the variances in PM 2.5 than meteorology , and 12 % of the observed PM 2.5 decrease was attributed to meteorological drivers in China since 2013 (Zhai et al., 2019). However, modeling and observational studies have mostly assessed the contribution of emissions and meteorology to regional PM 2.5 variations in the emission source regions with high anthropogenic emissions of air pollutants, and there have been few assessments of multi-scale changes of atmospheric environment over receptor regions in the regional transport of air pollutants.
The Twain-Hu Basin (THB), featuring the lowlands (mainly less than 200 m a.s.l.) of two central Chinese provinces, Hubei and Hunan (Fig. 1), is surrounded by the high air pollution regions of NCP, YRD, PRD and SB. As such, it is a key receptor region in the regional transport of air pollutants from the upstream region driven by East Asian monsoonal winds over CEC (Shen et al., 2020). Heavy air pollution in the THB with a unique "non-stagnation" atmospheric boundary layer is aggravated by the regional PM 2.5 transport over CEC Yu et al., 2020). In combination with the heavy pollution region of NCP through distinct transport channels, the regional transport from northern China to central China contributed 70.5 % PM 2.5 concentrations to a wintertime heavy pollution episode in the THB (Hu et al., 2021). Thus, the contributions of air pollutant emissions and meteorological conditions to air quality change over this region in central China need to be specifically assessed with the long-term observations over recent years.
In this study, we investigated the multi-scale changes of PM 2.5 concentrations over the THB, a key receptor region of regional PM 2.5 transport in China, from 2015 to 2019 by establishing a statistic model with a Kolmogorov-Zurbenko (KZ) filter. We then evaluated the contributions of anthropogenic emissions and meteorological changes to the declining trends in PM 2.5 concentrations over this receptor region in regional PM 2.5 transport over CEC during the past 5 years of emission controls. The analysis of THB's multi-scale air quality changes can improve the understanding of the effects of emission mitigation and meteorological changes on environmental change with regional transport of air pollutants.

Data
In order to analyze air quality changes in the THB, the observational data of hourly NO 2 , SO 2 and PM 2.5 concentrations from 2015 to 2019 were collected from the national air quality monitoring network (http://www.mee.gov.cn/, last access: 17 January 2022). The air quality observation data are under quality control, based on China's national standard of air quality observation.
The data of meteorological observations in the THB were sourced from the weather monitoring network of China Meteorological Administration (http://data.cma.cn/, last access: 17 January 2022), including air temperature (T ), relative humidity (RH), sea level pressure (SLP), wind speed (WS) and precipitation (Pre), with temporal resolutions of 3 h.

KZ filter
To better understand the multiple timescale variations of PM 2.5 and the relations to air pollutant emissions and meteorological drivers, a KZ filter (Rao and Zurbenko, 1994;Seo et al., 2018) is used to separate the daily data into multiscale components, based on an iterative moving average that removes high-frequency variations in the data with applications in the study of air pollutants, especially O 3 and PM 2.5 variations Ma et al., 2016;Seo et al., 2014;Zheng et al., 2020).
The KZ filter KZ m,p with a length of moving average window m and number of iterations p can remove the highfrequency component of periods smaller than the effective filter width N (≥ m×p 1/2 ). The KZ filter is applicable to the time series with missing data owing to the iterative moving average process, which provides high accuracy when compared with the wavelet transform method (Eskridge et al., 1997). By comparing different sets of moving average m and number of iterations p, it was found that the decomposed time series using the KZ 15,5 filter exhibited no white noise (short-term component), and the trend of the long-term component derived with the KZ 365,3 filter corresponded approximately to the interannual trend of the original data (Rao and Zurbenko, 1994;Eskridge et al., 1997). Based on the spectral decompositions of the daily observational data and three components, the power spectra of daily observational data in periods of less than 33 d and longer than 632 d (1.7 years) have been well reproduced by short-term and long-term components, and the seasonal component well represents the seasonal variations, i.e., periods between 33 d and 1.7 years . Thus we applied KZ 15,5 and KZ 365,3 filters to remove the variations with periods shorter than 33 d and 1.7 years in this study.
A meteorological or environmental variable X (t) observed in time series t can be decomposed into the short-term component X ST (t) and the baseline component X BL (t) presenting as (1) The baseline component X BL (t) is obtained by applying the KZ (15,5) filter to X (t), removing the short-term component X ST (t) with the temporal period shorter than 33 d from the observed data, expressed as X BL (t) = KZ (15,5) The baseline component X BL (t) also can be separated into the daily climatic averages X clm BL over the study period, occupying most of the seasonality in X BL (t) and the residual ε(t): To obtain the long-term component X LT (t) by removing the variations with the temporal period shorter than 1.7 years, the KZ 365,3 filter is applied to ε(t) expressed as follows: with the short-term component and the seasonal component The KZ filter was used to separate the daily surface PM 2.5 , NO 2 and SO 2 concentrations into short-term, seasonal and long-term components in this study. The short-term component presents a synoptic-scale variation of meteorological influences, which could control local accumulation and 3582 X. Sun et al.: Meteorology impact on PM 2.5 change over a receptor region regional transport of air pollutants (Seo et al., 2017), partly associated with short-term fluctuations in air pollutant emissions (Russell et al., 2010). The seasonal and long-term components are attributable to the variations in air pollutant emissions related to human activities as well as the seasonal and interannual changes in meteorological conditions .

Multiple linear regression of air pollutant changes with meteorological variables
By altering the local accumulation, regional transport, chemical conversion and wet and dry depositions of air pollutants, the meteorological factors such as wind, RH, T , air pressure and Pre could exert significant impacts on PM 2.5 changes (Sun et al., 2013;Z. Y. Chen et al., 2020). Therefore, with the multiple factors of the baseline components of 10 m WS, 2 m RH, 2 m T , SLP and Pre calculated by Eq.
(2), a multiple linear regression equation was stepwise established for the baseline component of PM 2.5 as follows: where MET BL i (t) (i [1,5]) is the baseline component of the meteorological variable i with i = 1, 2, 3, 4, 5, respectively, for WS BL (t), RH BL (t), T BL (t), SLP BL (t), Pre BL (t). We fit the regression coefficient a i for each meteorological variable and the intercept a 0 . The residual ε PM 2.5 between PM 2.5BL and PM 2.5BL MLR regressed with the multiple linear Eq. (7) is given as ε PM 2.5 contains not only the variability of PM 2.5 related to long-term changes in air pollutant emissions but also the minor seasonal change of PM 2.5 attributable to unconsidered meteorological influences in the multiple linear regression. By removing the minor seasonal change from ε PM 2.5 with the KZ 365,3 filter, the emissions-related long-term component PM 2.5 emiss LT (t) can be isolated as follows: Here the long-term component of surface PM 2.5 concentrations can be further separated into the emission-and meteorology-related long-term components with Eqs. (9) and (4) . Similarly, the multiple timescale variations in SO 2 and NO 2 with long-term variations related to changes in air pollutant emissions and meteorological drivers are decomposed by the KZ filter with multiple linear regression. Seo et al. (2018) described the details of this method. and (c) long-term components, and (d) total contribution to the total variances of daily PM 2.5 changes observed at 14 sites in the THB with regional averages of 47.5 %, 41.4 %, 3.7 % and 92.7 %, respectively.

Verification of PM 2.5 decompositions in multi-scale variations
The daily PM 2.5 concentrations observed at 14 sites over the THB (Fig. 1) were decomposed into short-term, seasonal and long-term components with Eqs. (4), (5) and (6) of the KZ filter. To verify the decomposition results, the spatial distribution of the total contributions of the shortterm, seasonal and long-term PM 2.5 components to the total variances of observed daily changes in PM 2.5 concentrations over 2015-2019 are shown in Fig. 2a. The larger the total variance, the more independent the three components are of each other . The sum of the long-term, seasonal and short-term components contributed 91.4 %-94.4 % of the total variance, with regional averages of 92.7 % (Fig. 2), reflecting a satisfactory verification of the KZ filtering results. Based on the PM 2.5 decomposition results of the KZ filter, the short-term, seasonal and long-term components, respectively, accounted for 34.8 %-53.8 %, 29.2 %-56.3 % and 0.2 %-9.8 % of the total variances of daily PM 2.5 changes in the THB over recent years (Fig. 2b, c and d), reflecting the different patterns of multiple timescale variations of PM 2.5 over this region in central China with diverse effects of emissions and meteorology. The regional contributions of the averaged short-term, seasonal and long-term compo-nents were, respectively, with 47.5 %, 41.4 % and 3.7 % of the daily PM 2.5 changes over the THB (Fig. 2), and thus it could be reasonably verified that the daily variation in atmospheric pollutant was generally dominated by the shortterm and seasonal components, with the long-term component determining the change trend (Ma et al., 2016;Yin et al., 2019a).
The short-term, seasonal and long-term PM 2.5 components were averaged for 14 sites of the THB to characterize the temporal variations of three components in the THB for 2015-2019 (Fig. 3). The correlation coefficients of 0.05, 0.01 and 0.04 among the decomposed short-term, seasonal and long-term components were near zero, indicating the orthogonal decomposition of multiple timescale components (Eskridge et al., 1997). According to the decomposed long-term, seasonal and short-term components demonstrated in Fig. 3, the notable peaks of decomposed seasonal and short-term components were highly consistent with the peaks of PM 2.5 concentrations in the original observed data, which further proved the reasonable decomposition of the multi-scale components of PM 2.5 change over 2015-2019.
The observed daily PM 2.5 exhibited a distinct daily variation, with an overlapping of high-frequency variations, which could be caused by mesoscale and synoptic-scale meteorological processes (Ma et al., 2016). The short-term component of PM 2.5 fluctuated frequently with a significantly positive correlation to the daily change of PM 2.5 (r = 0.68, p<0.05), indicating the important role of the short-term component with the temporal period <33 d in the day-to-day variations of PM 2.5 concentrations in the THB (Fig. 3a).
The notable peaks of PM 2.5 seasonal component that emerged in winters were highly consistent with the peaks of observed daily PM 2.5 concentrations (Fig. 3b). A close linkage with the significant correlation coefficient of 0.75 (p<0.05) was found between the changes of PM 2.5 seasonal component and daily PM 2.5 concentrations, which could reflect a significant modulation of the PM 2.5 seasonal oscillations in the day-to-day variations of PM 2.5 , driven by the seasonal shift of East Asian summer and winter monsoons, as well as the seasonal change of anthropogenic emissions (Zhu et al., 2012;Jeong and Park, 2017). The change of the long-term component of PM 2.5 exhibited a steadily declining trend over 2015-2019 (Fig. 3c), which was consistent with the interannual trend of observed PM 2.5 concentrations under the sustained impact of emission controls Xu et al., 2020). The correlation coefficient (r = 0.24, p<0.05) of the long-term PM 2.5 component with the observed daily PM 2.5 change was much smaller than those of the short-term and seasonal PM 2.5 components, implying less influence of emission reduction on the daily PM 2.5 change and air pollution frequency, although the declining trend in PM 2.5 was determined by anthropogenic emission reduction.
In previous studies, chemical transport models and statistical methods were both used to assess the changes in air pol-lution attributable to emissions and meteorology (Xiao et al., 2021). Significant declines in emissions-related PM 2.5 concentrations occurred in central China , and the meteorology offset the impact of emission reduction in typical years of unfavorable meteorological conditions (Xu et al., 2020;Gong et al., 2021). The regional averaged emission-and meteorologyrelated long-term components and the long-term component over the THB are displayed in Fig. S1a in the Supplement, and suggest the steadily declining trend of PM 2.5 and the dominant impact of emission reduction on long-term PM 2.5 changes, which is consistent with previous studies using multiple linear regression models for central China (Fig. S1b). The meteorology-related long-term component has a positive value in certain periods, implying the significant modulating effect of meteorology on PM 2.5 decline in the THB.
3.2 Multiple linear regressions of PM 2.5 , SO 2 and NO 2 with meteorological drivers Since the short-term variations in meteorological variables were excluded, the correlations between baseline components of PM 2.5 and meteorological variables were only related to their seasonal and long-term components, affected by the regional climate of East Asian monsoons rather than synoptic-scale meteorological processes. Based on our understanding of chemical and physical processes of diffusive transport, chemical transformation, emissions and depositions of PM 2.5 in the atmosphere, the dominant meteorological factors for changing PM 2.5 concentrations over china are wind speed, relative humidity, air temperature, atmospheric pressure and precipitation (Z. Y. . We examined the significant correlations between baseline components of air pollutant concentrations and selected a set of meteorological factors, including air temperature, wind speed, precipitation, relative humidity and air pressure (Tables S1-S3 in the Supplement). The meteorological parameters selected in this study are consistent with previous studies (Z. Y. . Generally, the baseline components of air pollutants were negatively correlated with the baseline components of wind speed (WS BL ) and positively correlated with the baseline components of sea level pressure (SLP BL ; Tables S1-S3), which could be attributed to the ventilation effect of wind and the stagnant conditions of meteorology in high-pressure systems, restraining the horizontal and vertical dispersions of air pollutants (Hsu and Cheng, 2016;Wang et al., 2016;Miao et al., 2017). Although wind speed exerts a negative influence of on PM 2.5 concentrations over the emission source region, increasing wind speed might cause the accumulation of PM 2.5 concentrations over the downwind region of emission sources (Z. Y. , which led to the inconsistent influence of WS BL in the region of central China (Tables S1-S3). Under high surface air temperature conditions, there are strong thermal activities such as turbulence, mak- ing an accelerated dispersion of air pollutants (Y. Yang et al., 2016). The negative influence of RH BL and T BL on PM 2.5BL , SO 2BL and NO 2BL mainly reflected the effect of the seasonal cycle in East Asian winter and summer monsoons, whereas the influence of precipitation on air pollutants was more straightforward than other meteorological parameters, negatively influencing surface pollutant concentrations through the precipitation washout of air pollutants (Tables S1-S3).
To isolate emissions-related long-term components from long-term components of PM 2.5 , NO 2 and SO 2 , the stepwise multiple linear regressions of PM 2.5BL , SO 2BL and NO 2BL , respectively, with the baseline components of the meteorological parameters (T BL , WS BL , RH BL , SLP BL and Pre BL ) were conducted with Eq. (7) for 14 sites, by adding and deleting meteorological variables based on their independent statistical significance to obtain the best model fit (Draper, 1998). We evaluated PM 2.5BL , SO 2BL and NO 2BL fitted by multiple linear regression models with KZ decomposition (Table 1). The multiple linear regressions explained PM 2.5BL , SO 2BL and NO 2BL with the adjusted determination coefficients (Adj. R 2 ) of 0.5695-0.8093, 0.0630-0.4592 and 0.6304-0.8669, passing the confidence level of 99 % in all the THB sites, confirming the reasonable construction of the multiple linear regressions. The Adj. R 2 of multiple linear regression for SO 2BL were lower than those of PM 2.5BL and NO 2BL , which might be attributed to the larger impact of SO 2 emission controls on the seasonal and long-term SO 2 variations. In general, the variations of meteorological drivers can well reproduce the meteorology-related seasonal and longterm variations of PM 2.5 , SO 2 and NO 2 in the THB (Table 1).

Interannual variations in air pollutants observed over
the THB PM 2.5 consists of chemical components generated in the complex physical and chemical processes (S. . Primary particles are emitted directly from anthropogenic (e.g., industry, power plants and vehicles) and natural (e.g., outdoor biomass burning and dust storms) sources. Secondary particles (e.g., sulfate and nitrate) are converted in the chemical reactions of the precursor gases (e.g., SO 2 and NO x ), which are mainly produced by human activities (S. . Therefore, in addition to the reductions in primary particulate emissions, control of the secondary aerosol precursor emissions is of great importance in mitigating air pollution. The interannual variations of the ratios in annual mean PM 2.5 , SO 2 and NO 2 concentrations relative to the annual averages in 2015 over the THB are displayed in Fig. 4. The declines of PM 2.5 and SO 2 in 2019 averaged over the THB were −26 % and −68 % relative to 2015, while the decrease ratio in NO 2 was only −8 % over this region. The observed SO 2 concentrations had a steeper decrease than PM 2.5 and NO 2 , possibly because the dominant source sectors (i.e., power and industry) of SO 2 significantly reduced their emissions . The power sector was the major contributor to emission reduction but only accounted for onethird of NO x emissions and the contribution of transportation to NO x emissions was estimated to have increased over recent years . The interannual variations in emissions for China were calculated from MEIC , the annual total emissions of SO 2 and NO x and the PM in the THB region reported by the National Bureau of Statistic of China (http://www.stats.gov.cn/tjsj/ndsj/, last access: 17 January 2022), presenting a more rapid decline of SO 2 emissions in the THB than the changes of PM 2.5 and NO x emissions (Fig. S2 in the Supplement). The declining trend of anthropogenic emissions estimated from emission inventories can support the explanation of the changes in air pollutant concentrations. Figure 5 shows the spatial distributions of 5-year averaged concentrations, the linear trends and the change rates in interannual variations of PM 2.5 , SO 2 and NO 2 observed in the THB over 2015-2019. The change rates (% yr −1 ) were calculated with the linear trends by dividing with temporalmean concentrations of air pollutants at the observation sites for the analysis period in Fig. 5. The 5-year averaged PM 2.5 concentrations over the THB exceeded the Chinese National secondary air quality standard of 35 µg m −3 for annual mean PM 2.5 concentrations (Fig. 5a), while SO 2 and NO 2 concentrations reached the secondary standards of 60 and 40 µg m −3 in annual mean SO 2 and NO 2 concentrations at most sites over the THB ( Fig. 5d and g). Specifically, the 5-year averaged NO 2 concentrations exceeded 40 µg m −3 in WH (Wuhan), a megacity in central China, which might be attributable to the large amounts of traffic transportation. From 2015 to 2019, both PM 2.5 and SO 2 decreased at all sites over the THB (Fig. 5b and e), whereas NO 2 trends changed from mostly negative to positive in some sites (Fig. 5h), possibly due to the spatial disparity of NO x emissions in traffic sectors . The comparison among the change rates of PM 2.5 , SO 2 and NO 2 in the THB presented the largest decreases of SO 2 with −40 %-−20 % yr −1 over the 5 years (Fig. 5c, f and i), reflecting the effective control of SO 2 emissions in terms of primary gaseous pollutants.
There were obvious decreases in regional mean PM 2.5 , SO 2 and NO 2 concentrations over the THB (Fig. 4), while the declining degree of PM 2.5 and SO 2 varied from site to site over the THB and the change trends in NO 2 were weak negative and even positive in certain sites (Fig. 5c, f and  i). These interannual changes of air pollutants in the THB over recent years were investigated using the emission-and meteorology-related long-term components of the air pollutants in the next sections.

Effects of NO 2 and SO 2 emission reductions on PM 2.5 change trends
The declining trend of PM 2.5 in China could be partly attributed to the reduced NO x and SO 2 concentrations for producing secondary aerosols . The reduction rates of anthropogenic emissions markedly accelerated after 2013, decreasing by 59 % for SO 2 , 21 % for NO x and 33 % for PM 2.5 during 2013-2017 over China . In order to assess the effect of changing precursor pollutant emissions on PM 2.5 reductions, we compared the linear trends of emissions-related long-term components of PM 2.5 , NO 2 and SO 2 decomposed based on Eq. (9) over the THB for 2015-2019 (Fig. 6). The distinct declining trends of emissions-related long-term PM 2.5 and SO 2 components, as well as the variable trends of emissions-related long-term NO 2 components were distributed basically consistently with the positive and negative trends in the interannual variations of air pollutant concentrations in the THB (middle column of Fig. 5; Fig. 6), demonstrating that the local emissions of air pollutants could spatially dominate the long-term variations of air pollutants in central China, especially the increasing trends in NO 2 at some THB sites. PM 2.5 concentrations are changed by emissions of both primary PM 2.5 and its gaseous precursors. As major gaseous precursors, SO 2 and NO 2 can be oxidized to convert nitrate and sulfate for secondary PM 2.5 (S. . To investigate the effects of emission reductions on the interannual variations of PM 2.5 , NO 2 and SO 2 over recent years, the ratios of change trends in long-term (k LT ) and emissionsrelated long-term (k emiss ) components of PM 2.5 , SO 2 and NO 2 , in the THB over 2015-2019 were demonstrated in Fig. 7, where the long-term and emissions-related long-term components of PM 2.5 , SO 2 and NO 2 were calculated with Eqs. (4) and (9). The trend ratios k LT /k emiss <1 indicated the more obvious downward trend of emissions-related longterm variations than the long-term trend of air pollutant concentrations, which might be attributed to the offsetting effect of meteorological conditions on emission reduction in the air quality change, whereas the long-term trend of air pollutant concentrations was more significant than the emissionsrelated long-term trend as k LT /k emiss >1, reflecting the synchronous impacts of anthropogenic emissions and meteorology on the long-term trend in air pollutant change. In addition, the trend ratios k LT /k emiss >1 and k LT /k emiss <1 of the gaseous precursors of PM 2.5 , SO 2 and NO 2 , could reflect the high and weak efficiencies of SO 2 and NO 2 , converting to sulfate and nitrate in the production of secondary PM 2.5 during air pollutant emission reduction. The notable differences in Fig. 7 were spatially distributed with the trend ratios k LT /k emiss >1 and k LT /k emiss <1 in PM 2.5 , SO 2 and NO 2 concentrations under the same meteorological conditions, indi- cating the different influences of emissions on the long-term variations of PM 2.5 , SO 2 and NO 2 in the THB during recent years. The reduction in PM 2.5 emissions was a primary cause of the long-term declines in PM 2.5 concentrations in the THB, even though the meteorological changes might offset the effects of emission reduction on air quality improvement over the southern THB (Figs. 6 and 7). It is noteworthy that the trend ratios k LT /k emiss <1 of PM 2.5 were accompanied with k LT /k emiss >1 of SO 2 and NO 2 at the downwind southern THB sites with both negative k LT and k emiss (Fig. 7, Table S4 in the Supplement), which could imply the increasing conversion efficiency of SO 2 and NO 2 to sulfate and nitrate for secondary PM 2.5 during the reductions of air pollutant emissions over recent years. In the upwind northern THB sites, the k LT /k emiss >1 of PM 2.5 were accompanied with k LT /k emiss >1 of SO 2 and NO 2 with an obvious facilitating effect of meteorology on the PM 2.5 decline (Fig. 7,  Table S4), revealing the underlying effect of regional transport of air pollutants on the spatial distribution of conversion efficiency of gaseous precursor to secondary PM 2.5 .
In order to further assess the effect of gaseous precursor emissions on PM 2.5 declines during recent 5-year air pollution mitigation, we selected seven and nine sites in the THB with decreasing trends of emissions-related long-term SO 2 and NO 2 components below −0.5 and 0.0 µg m −3 100 d −1 , respectively, (Table S4) to compare the trend ratios k LT /k emiss of PM 2.5 , NO 2 and SO 2 for 2015-2019 (Fig. 8). The significantly negative linear correlations between the changes in k LT /k emiss of gaseous precursors (SO 2 , NO 2 ) and PM 2.5 could present the connection of k LT /k emiss >1 of NO 2 and SO 2 with k LT /k emiss <1 of PM 2.5 , which confirms the fact that the high conversion efficiency of SO 2 and NO 2 to sulfate and nitrate could have have offset the role of PM 2.5 emission reduction in controlling PM 2.5 pollution. The decreasing emissions of gaseous precursors drove faster oxidation of NO 2 and SO 2 to nitrate and sulfate components of PM 2.5 in the source regions of air pollution in China (Zhai et al., 2021;Huang et al., 2021). This study identified the enhancing contribution of gaseous precursors to PM 2.5 concentrations with reduced anthropogenic emis-sions of air pollutants over receptor regions in regional PM 2.5 transport.
There are a few sources of uncertainty in the discussion of chemical transformation, for example in the separation of emission-and meteorology-related long-term components and in the selection of observational sites. Our results point to a better understanding of the offsetting effect of SO 2 and NO 2 oxidized to secondary particles in PM 2.5 mitigation during the emission reduction in the THB. The long-term changes in PM 2.5 are also caused by the emission variations of primary components like black and organic carbon, in addition to the chemical transformation of gaseous precursors. The difference in the emission of different primary pollutants may also lead to modifications in k LT k emiss of PM 2.5 . However, due to the current lack of long-term observations of PM 2.5 components in the THB, the influence of emission variations of primary components on long-term changes in PM 2.5 concentrations is not assessed in our study. Further work with long-term observational data of PM 2.5 components like black and organic carbon could be conducted to quantify the influence of emissions of primary components and chemical transformation of gaseous precursors on PM 2.5 changes.

Meteorological contribution to PM 2.5 change trends
As the air pollutant change trend is assumed to generally consist of emission-and meteorology-related changes Yin et al., 2019b), the meteorological contribution rate Con met to the long-term PM 2.5 change trend is calculated with the following equation: Here, Con met (in %) is estimated with the linear trends k LT of the long-term component PM 2.5LT (t) and k emiss of the emissions-related long-term component PM 2.5 emiss LT (t). PM 2.5LT (t) and PM 2.5 emiss LT (t) are, respectively, calculated with Eqs. (4) and (9).  To quantitatively assess the meteorological contribution to the declining PM 2.5 trends, the linear trends k LT and k emiss with the meteorological contribution rate Con met in Eq. (10) were presented in Table S5 in the Supplement for 14 sites over the THB during 2015-2019. All the trends k LT and k emiss , respectively, in PM 2.5LT (t) and PM 2.5 LT emiss (t) were negative over the THB (Table S5), indicating the significant effect of emission reductions on declining PM 2.5 trends for improving regional air quality in central China. By comparing the declining PM 2.5 trends k emiss and k LT from site to site (Table S5), the positive and negative contributions of meteorological variations to PM 2.5 change trends over recent years were determined with the positive and negative differences between k emiss and k LT with the distinct meteorological influences on the change of the THB regional environment.
The spatial distribution of meteorological contribution rates Con met regarding the long-term declining PM 2.5 trend presented a unique pattern of northern positive and southern negative values over the THB (Fig. 9), with high positive contributions in the northern sites XY (61.92 %) and EZ (37.31 %), and as low negative contributions in the southern sites CD (−24.93 %) and CS (−23.03 %). It is worth mentioning that the contribution rates of meteorological variations show great spatial disparities at a small scale, i.e., EZ, HG and HS, which seems not be induced by the variation in synoptic weather or meteorological conditions. The underlying surface conditions dominate the near-surface meteorological conditions in the atmospheric boundary layer at a small scale . The topography and land use of HG, HS, EZ and the surrounding regions vary distinctly with the underlying surface conditions of plains, lakes and hilly areas. The underlying surface of observational sites with different near-surface meteorology effectively influences the local accumulation, chemical transformation, and dry and wet depositions of air pollutants (Bai et al., 2022). Therefore, the heterogeneity of meteorological contributions to PM 2.5 at such a small spatial scale might be attributed to the local meteorological conditions in the atmospheric boundary layer, which are largely affected by the underlying surface changes. Compared with the statistical studies using synthetic data of meteorological influence on regional PM 2.5 changes in central China and the meteorological contribution from −45.5 % to 29.0 % over recent years (Gong et al., 2021;, the PM 2.5 pollution over the THB was affected contrarily by meteorological drivers with the northern positive and southern negative contributions from 2015 to 2019 (Fig. 9). The meteorological change could accelerate and offset the effects of emission reductions on declining PM 2.5 trends in the northern and southern THB, which might be attributed to the regional transport of air pollutants conducive to the upwind diffusion and downward accumulation of air pollutants, respectively, over the northern and southern THB under the declining wind of East Asian monsoons over recent years (Hu et al., 2020;Zhong et al., 2019).
3.6 Meteorological contribution to PM 2.5 changes validated with WRF-Chem modeling The above observational study investigated the meteorological influence on the changes in PM 2.5 concentrations in the THB using the KZ filter, concluding on the large impact of meteorology on the PM 2.5 changes over 2015-2019. To validate this conclusion of analyses with the KZ filter, we designed three sets of modeling experiments CTRL, SENS-MET and SENS-EMI (Table S6 in  Air pollutant emission inventories, modeling configuration, experiment design and modeling verification are described in the Supplement. The modeling verification of experiments CTRL indicated that PM 2.5 and meteorology were reasonably reproduced by the WRF-Chem simulation (Figs. S4-S5, Table S7 in the Supplement), and the three designed sets of modeling experiments CTRL, SENS-MET and SENS-EMI can be used in the further analyses of emission and meteorological impact on PM 2.5 change over 2015-2019 to confirm the results of the KZ filter.
We derived the impact of meteorology by comparing the simulated PM 2.5 concentrations in three sets of experiments CTRL, SENS-MET and SENS-EMI (Table S6). The relative contribution of meteorology to the interannual changes of PM 2.5 concentrations was calculated with a linear additive relationship of contributions of meteorology and emission in the following equations: RCon MET = Con MET Con MET + Con EMI × 100 %.
k CTRL , k MET and k EMI represent the trends in interannual changes of PM 2.5 concentrations simulated by the experiments CTRL, SENS-MET and SENS-EMI, respectively. Con MET and Con EMI are the contributions of meteorology and emission, and RCon MET is the contribution rate (%) of meteorology to interannual changes of PM 2.5 concentrations . Based on WRF-Chem modeling experiments, we assessed the impact of meteorological changes on interannual PM 2.5 variations from 2015 to 2019 with Eqs. (11)-(13). The relative contribution of meteorology to interannual PM 2.5 variations displayed a regional pattern of northern positive and southern negative values over the THB (Fig. 10), confirming the impact of meteorological changes by accelerating and offsetting the effects of emission reductions on declining PM 2.5 trends in the northern and southern THB, respectively. The general spatial distribution of meteorological contribution rates regarding declining PM 2.5 trends from the WRF-Chem simulation was consistent with the results using the KZ filter (Figs. 9 and 10), validating these results and confirming that meteorological drivers exerted a contrary impact between the northern positive and southern negative contributions to long-term changes of PM 2.5 concentrations in the THB.

Conclusions
The meteorological effect on multi-scale changes of the atmospheric environment has been rarely assessed for the re- ceptor region in the regional transport of air pollutants. In this study of observations and modeling, we targeted the THB, a large region of heavy PM 2.5 pollution in central China, to assess the effects of meteorology on PM 2.5 changes over a receptor region in the regional transport of air pollutants during recent years. The study results provide insights into the effects of emission mitigation and meteorological changes on the source-receptor relationship of the long-range transport of air pollutants for regional and global environment changes.
This study decomposed the observed PM 2.5 concentrations into multiple timescale components with a modified KZ filter to better understand the PM 2.5 variations, with the short-term, seasonal and long-term components accounting for, respectively, 47.5 %, 41.4 % and 3.7 % of observed PM 2.5 changes. The short-term and seasonal PM 2.5 components dominated the daily PM 2.5 changes and longterm component determined the trend of PM 2.5 change over recent years. The long-term components of PM 2.5 , SO 2 and NO 2 were further isolated into emission-and meteorologyrelated long-term components with multiple linear regressions in order to determine the contributions of emission and meteorology to PM 2.5 decline in the THB over 2015-2019. The reduction in anthropogenic emissions was the primary cause of the long-term decline in PM 2.5 concentrations, and the meteorological changes moderated the PM 2.5 variations in the THB. As the receptor region of regional PM 2.5 transport, the impact of diverse meteorological conditions on long-term trend of PM 2.5 changes displayed a unique regional pattern of northern positive rates up to 61.92 % and southern negative rates down to −24.93 %. The change of meteorological conditions could accelerate and offset the effects of the emission reductions on declining PM 2.5 trends in the northern and southern THB, which can be attributed to the upwind diffusing and downward accumulating roles of regional transport pathways in air pollutants in the THB. In terms of gaseous precursor emissions, the increasing conversion efficiency of SO 2 and NO 2 to sulfate and nitrate for secondary PM 2.5 could have offset the role of emission reduction in controlling air pollution, and the contribution of gaseous precursors to secondary PM 2.5 was enhanced with the reducing anthropogenic emissions of air pollutants over this receptor region.
This study revealed the impact of anthropogenic emissions and meteorological conditions on PM 2.5 decline over a receptor region in the regional transport of air pollutants in central China. The effect of regional transport on PM 2.5 pollution over receptor regions is found to differ from that over source regions with high anthropogenic emissions. We considered the uncertainties induced by statistical methods as systematic biases and explained the offsetting effect of enhancing oxidation of gaseous precursors to secondary particles on PM 2.5 decline during stringent emission controls. The changes in data coverage and meteorological parameter selection largely influence the quantitative estimation of the contributions of meteorology and emissions. Further work should be undertaken regarding climate analysis of the longterm data of fine meteorological and environmental observations and more comprehensive modeling of chemical and physical processes in the atmosphere to generalize the assessment of the effects of emission mitigation and meteorological changes on the source-receptor relationship of the regional transport of air pollutants.