Impact of Meteorological Conditions on PM 2.5 Pollution in Changchun and Associated Health Risks Analysis

: The escalating concern regarding increasing air pollution and its impact on the health risks associated with PM 2.5 in developing countries necessitates attention. Thus, this study utilizes the WRF-CMAQ model to simulate the effects of meteorological conditions on PM 2.5 levels in Changchun, a typical city in China, during January 2017 and January 2020. Additionally, it introduces a novel health risk-based air quality index (NHAQI) to assess the influence of meteorological parameters and associated health risks. The findings indicate that in January 2020, the 2-m temperature (T2), 10-m wind speed (WS10), and planetary boundary layer height (PBLH) were lower compared to those in 2017, while air pressure exhibited a slight increase. These meteorological parameters, characterized by reduced wind speed, heightened air pressure, and lower boundary layer height—factors unfavorable for pollutant dispersion—collectively contribute to the accumulation of PM 2.5 in the atmosphere. Moreover, the NHAQI proves to be more effective in evaluating health risks compared to the air quality index (AQI). The annual average decrease in NHAQI across six municipal districts from 2017 to 2020 amounts to 18.05%. Notably, the highest health risks are observed during the winter among the four seasons, particularly in densely populated areas. The pollutants contributing the most to the total excess risk (ERtotal) are PM 2.5 (45.46%), PM10 (33.30%), and O3 (13.57%) in 2017, and PM 2.5 (67.41%), PM10 (22.32%), and O3 (8.41%) in 2020. These results underscore the ongoing necessity for PM 2.5 emission control measures while emphasizing the importance of considering meteorological parameters in the development of PM 2.5 reduction strategies.


Introduction
In recent years, rapid urbanization and industrialization have propelled the country's economy forward and enhanced the convenience of people's lives.However, this progress has also led to the emergence of air pollution issues in many Chinese cities, particularly concerning fine particulate matter (PM 2.5 ) pollution [1,2].PM 2.5 refers to fine particulate matter with a diameter less than 2.5 µm, originating primarily from human activities' primary PM 2.5 emissions and natural sources, as well as secondary PM 2.5 generated through chemical reactions with other substances upon entering the atmospheric environment [3,4].Numerous studies have demonstrated that meteorological conditions play a significant role in influencing the mass concentration of PM 2.5 and, consequently, urban air quality.For instance, a study conducted in Northern China revealed a positive correlation between relative humidity and PM 2.5 levels, whereas wind speed exhibited a negative correlation with PM 2.5 [5].In addition to relative humidity and wind speed, PM 2.5 levels were found to be negatively correlated with surface temperature and planetary boundary layer height (PBLH) in previous studies [5,6].The study by Zhai et al. confirmed that meteorological conditions significantly contribute to the variability of air quality in many regions of China [4].Air quality models are frequently employed to gain a comprehensive understanding of and forecast air pollution issues.For instance, the Weather Research and Forecasting/Community Multiscale Air Quality (WRF-CMAQ) model enables the quantification of the impacts of anthropogenic emissions and meteorological conditions on PM 2.5 [7,8].Related studies have demonstrated that the interannual variation in atmospheric pollutant concentrations (e.g., PM 2.5 , NO 2 , and O 3 ) and inter-city pollution transport processes in the Beijing-Tianjin-Hebei and Eastern provinces regions are influenced not only by anthropogenic emissions but also by meteorological conditions that cannot be overlooked [9][10][11].
Furthermore, excessive air pollutants not only degrade the quality of the air environment but also exert a significant impact on human health and climate change [12,13].Past epidemiological research has indicated that atmospheric pollutants contribute to a range of human diseases, including lung diseases, respiratory ailments, cardiovascular disorders, and even malignancies [14][15][16].Moreover, they lead to elevated morbidity and premature mortality among exposed populations.A study conducted in the Beijing-Tianjin-Hebei region revealed that exposure to PM 2.5 resulted in 120,000 premature deaths (95% CI: 80,000, 150,000) in the region in 2013 [17].In a nationwide study conducted across China, it was determined that 1.8 million premature deaths were attributable to PM 2.5 exposure in 2017 (95% CI: 1.6 million, 2 million), marking a 31% increase compared to 2005 [18].Atmospheric pollutants also influence the spread of epidemics, as evidenced by several studies correlating the number of novel coronavirus infections (COVID-19) with air pollutants [19][20][21][22].
In 2012, China introduced the latest Ambient Air Quality Standards (AAQS), which provide detailed guidelines on the utilization of the Air Quality Index (AQI).AQI serves as a crucial indicator for assessing air quality, and numerous studies have explored the relationship between AQI changes and pollutants [23,24].However, due to AQI's reliance on the maximum calculated value among six pollutants, it may not adequately characterize overall air quality or assess the adverse effects of individual pollutants on public health.To address this limitation, Kyrkilis et al. and Hu et al. proposed the aggregate air quality index (AAQI) and health risk-based air quality index (HAQI) as alternatives to AQI [25,26].Subsequently, Ma et al. introduced an enhanced version of HAQI to improve its accuracy.Several subsequent studies utilizing HAQI have demonstrated its superior effectiveness over AQI in assessing population health [27][28][29].
The rapid onset and global dissemination of COVID-19, which emerged around the Chinese New Year in 2020, prompted stringent quarantine measures by the Chinese government.These measures included the implementation of town and city lockdowns and restrictions on transportation.Consequently, emissions of most air pollutants decreased, leading to significant changes in air quality.This provided a unique opportunity to investigate the impacts of anthropogenic emissions and meteorological conditions on air quality [30,31].Some studies have revealed that in certain areas, air pollution levels were more pronounced in January 2020 compared to previous years, indicating a correlation between reduced social activity and heightened pollution levels [32,33].This observation warrants further attention and investigation.
As the capital city of Jilin Province, Changchun is situated in the heart of Northeast China.Its primary industries revolve around automobile manufacturing and agricultural product processing.However, the city also grapples with environmental pollution and associated public health issues due to frequent biomass burning and high coal consumption, factors that cannot be overlooked [32].However, most recent studies have predominantly focused on regions such as Beijing-Tianjin-Hebei, the Yangtze River Delta, the Pearl River Delta, and Southwest China, leaving fewer studies conducted in other areas.In light of this gap, this paper selects Changchun as the study area and utilizes the WRF-CMAQ model to examine the influence of meteorological conditions on PM 2.5 concentrations in Changchun during January 2017 and January 2020.Additionally, it assesses the health risks associated with air pollutants in these two time periods based on the enhanced HAQI.The findings of this study aim to provide valuable insights for policymakers in formulating more effective emission reduction strategies.

Data Collection
The hourly mass concentration data of PM 2.5 pollutants utilized in this study were sourced from the publication platform of the China National Environmental Monitoring Center (CNEMC, http://106.37.208.233:20035/(accessed on 7 April 2023)).Specifically, the data were collected from 10 state-controlled automatic atmospheric environmental monitoring stations and one meteorological monitoring station located in Changchun (Table 1).All data underwent pre-processing in accordance with the Technical Regulation for Ambient Air Quality Assessment.This involved excluding spatial and temporal anomalies in the monitoring results to ensure data quality.The number of valid data obtained adhered to the relevant regulations outlined in the Ambient Air Quality Standards (AAQS).Ground-based meteorological monitoring data in Changchun were obtained from the China Meteorological Data Service Center (CMDSC, http://data.cma.cn(accessed on 10 April 2023)), and rigorous quality checks were conducted on the data to evaluate the WRF model.

Model Configuration
In this study, the WRF-CMAQ model was employed to simulate the spatial and temporal distribution of PM 2.5 concentration (versions WRF v3.9.1 and CMAQ v5.3.2). Figure 1 illustrates the three-level nested domains of the study area.The outermost domain encompassed the three northeastern provinces with a grid resolution of 27 km × 27 km; the middle domain covered Jilin province with a grid resolution of 9 km × 9 km; and the innermost domain focused on Changchun with a grid resolution of 3 km × 3 km.The initial and boundary conditions for the meteorological inputs in the WRF model were obtained from the National Centers for Environmental Prediction (NCEP) final analysis data (FNL), featuring a temporal resolution of 6 h and a spatial resolution of 1.0 • × 1.0 • .The gridded emission inventory required for the CMAQ model was derived from the Multiscale Emission Inventory Model for China (MEIC, http://meicmodel.org/(accessed on 16 March 2023)), developed by Tsinghua University.This study utilized the 2017 regional anthropogenic source pollutant and greenhouse gas (GHG) emission data from Tsinghua University's emission source inventory for China, encompassing SO 2 , NOx, CO, NH 3 , VOCs, PM 2.5 , PMcoarse, BC, OC, CO 2 , and other pollutants.VOCs were allocated based on the CB06 atmospheric chemistry mechanism.The inventory provided gridded emission data for various pollutants at a horizontal resolution of 0.25 • , categorized into five sectors: power, industry, residential, transport, and agriculture [34,35].To conform the emission inventory data to the input format required by the CMAQ model, the Inventory Spatial Allocation Tools (ISAT) were utilized for further spatial and temporal allocation of the emission inventory [36].Spatial allocation coefficients were established based on town grid data for electricity sources, road network grid data for transportation sources, farmland grid data for agricultural sources, industrial grid data for industrial and mining sources, and GDP grid data for residential sources.
Atmosphere 2024, 15, x FOR PEER REVIEW 4 of 20 anthropogenic source pollutant and greenhouse gas (GHG) emission data from Tsinghua University's emission source inventory for China, encompassing SO2, NOx, CO, NH3, VOCs, PM2.5, PMcoarse, BC, OC, CO2, and other pollutants.VOCs were allocated based on the CB06 atmospheric chemistry mechanism.The inventory provided gridded emission data for various pollutants at a horizontal resolution of 0.25°, categorized into five sectors: power, industry, residential, transport, and agriculture [32,33].To conform the emission inventory data to the input format required by the CMAQ model, the Inventory Spatial Allocation Tools (ISAT) were utilized for further spatial and temporal allocation of the emission inventory [34].Spatial allocation coefficients were established based on town grid data for electricity sources, road network grid data for transportation sources, farmland grid data for agricultural sources, industrial grid data for industrial and mining sources, and GDP grid data for residential sources.

Evaluation Indicators
In order to evaluate the effectiveness of WRF-CMAQ simulations, six statistical indicators, namely correlation coefficient I, mean fraction bias (MFB), mean fraction error (MFE), normalized mean bias (NMB), normalized mean error (NME), and Root Mean Square Error (RMSE), are used to validate the system simulation results.R indicates the degree of agreement between the simulated value and the trend of the monitored value, with a value closer to 1 indicating better simulation quality.MFB reflects the fractional bias of the simulated value from the mean of the monitored value, while MFE represents the mean absolute error of the simulated value from the mean of the monitored value.NMB reflects the average direction and degree of bias from the monitored value, while NME indicates the average absolute degree of bias from the monitored value for each simulated value.RMSE is a commonly used measure of the accuracy of a model's predictions and is used to measure the degree of discrepancy between observations and model predictions.The meteorological parameters evaluated included temperature at 2 m (T2), wind speed at 10 m (WS10), and wind direction at 10 m (WD10), and the pollutant evaluated was PM2.5.Emery et al. [35] and Boylan and Russell et al. [36] proposed that the model performance criteria are considered satisfactory when −60% ≤ MFB ≤ 60% and MFE ≤ 75%,

Evaluation Indicators
In order to evaluate the effectiveness of WRF-CMAQ simulations, six statistical indicators, namely correlation coefficient I, mean fraction bias (MFB), mean fraction error (MFE), normalized mean bias (NMB), normalized mean error (NME), and Root Mean Square Error (RMSE), are used to validate the system simulation results.R indicates the degree of agreement between the simulated value and the trend of the monitored value, with a value closer to 1 indicating better simulation quality.MFB reflects the fractional bias of the simulated value from the mean of the monitored value, while MFE represents the mean absolute error of the simulated value from the mean of the monitored value.NMB reflects the average direction and degree of bias from the monitored value, while NME indicates the average absolute degree of bias from the monitored value for each simulated value.RMSE is a commonly used measure of the accuracy of a model's predictions and is used to measure the degree of discrepancy between observations and model predictions.The meteorological parameters evaluated included temperature at 2 m (T2), wind speed at 10 m (WS10), and wind direction at 10 m (WD10), and the pollutant evaluated was PM 2.5 .Emery et al. [37] and Boylan and Russell et al. [38] proposed that the model performance criteria are considered satisfactory when −60% ≤ MFB ≤ 60% and MFE ≤ 75%, or when NMB ≤ ±30% and NME ≤ 50%, and R ≥ 0.40.The specific formulas for the six indicators are as follows [39]: where N is the number of samples, m i is the simulation result at moment i, O i is the monitoring result at moment i, m is the mean value of the simulation result, and O is the mean value of the monitoring result.

Health Assessment 2.4.1. AQI
Referring to the Ambient Air Quality Standards (AAQS) and the Technical Regulation on Ambient Air Quality Index, the air quality sub-index (AQI i ) for individual pollutant i is calculated using Equation (7), and the total AQI for six pollutants is calculated using Equation (8), as follows: where i is the pollutant category, j is the health risk category, C i,m is the observed concentration of pollutant i, AQI i,j and AQI i,j−1 are the maximum limits of AQI corresponding to the jth and j−1th risk category nearest to the observed concentration of pollutant i, and C i,j and C i,j−1 are the maximum limits of concentration corresponding to the jth and j−1th risk category nearest to the observed concentration of pollutant i, respectively.The range of AQI values and the corresponding pollutant concentration limits and health categories obtained from the Ministry of Ecology and Environment of China are shown in Table 2.The excess risk (ER) of the pollutant was introduced to represent HAQI [26].In this case, the relative risk of pollutant i (RRi) is calculated as follows: where β i is the exposure-response relationship coefficient that can represent the excess health risk associated with each unit increase in concentration of pollutant i.According to Shang et al. [40], each 1 µg/m 3 increase in the concentration of six pollutants: PM 2.5 , PM 10 , SO 2 , NO 2 , and O 3 , corresponds to 0.038%, 0.032%, 0.081%, 0.13%, and 0.048%, of β i ; each 1 mg/m 3 increase in the concentration of CO corresponds to β i is 3.7%.C i,m is the observed concentration of pollutant i and C i,0 is the baseline concentration of pollutant i, which is the maximum health concentration limit, using the upper limit of pollutant concentration of the AAQS 24-h secondary standard (Table 2).In particular, when C i,m is less than or equal to C i,0 , RR i = 1, which is considered a health risk at this time.The ER i for pollutant i is calculated as follows: The total ER for all six pollutants is calculated as follows: where ER total was classified as an arbitrary index between 0 and 10 in the study by Cairncross et al. [41] to represent the excess health risk from air pollution.to represent the equivalent concentration of pollutant i when ER i equals ER total [26], and if we use Equation ( 13) to calculate C * i,m directly, it may lead to high results in subsequent calculations of HAQI, so this study used the novel HAQI (NHAQI) instead of HAQI [29].Using the segmentation function, when C i,m ≤ C i,0 , i.e., the observed concentration of pollutant i does not cause additional health risk, the observed concentration Equation ( 14) is still used; when C i,m > C i,0 , i.e., the observed concentration of pollutant i causes additional health risk, Equation ( 13) is used to calculate.In addition, the equivalent relative risk (RR * i ) of pollutant i is calculated as follows: The equivalent concentration of pollutant i, C * i,m is used instead of C i,m to calculate NHAQI i , which is calculated as follows:

Scenario Settings
To evaluate the influence of meteorological conditions on the change of PM 2.5 in Changchun City in January 2017 and January 2020, this study utilized the emission inventory of 2017 and meteorological data of 2017 for WRF-CMAQ model simulation in January 2017.Similarly, for January 2020, the emission inventory of 2017 and meteorological data of 2020 were used for the WRF-CMAQ model simulation.This approach allows for the investigation of the impact of meteorological conditions on air quality in 2020 through simulations.The details are summarized in Table 3.In this study, the accuracy of the WRF model simulation is initially verified, focusing on meteorological parameters including T2, WS10, and WD10.The comparison of the WRF model simulation results with the monitoring results is presented in Table 4 and Figure 2. It is evident that the WRF model effectively replicates the peaks and trends of meteorological parameters.Regarding meteorological parameters in 2017, the simulated value of T2 in January was recorded as −16.72 • C, which is 3.57 • C lower than the observed value.Conversely, the simulated mean values of wind speed and direction in January are higher than the monitored mean by 0.62 m/s and 17.7 • , respectively, indicating minor differences.drastic meteorological changes in July 2020.Notably, the NMB for T2 in 2020 reaches 74.22%, which may be related to limitations in the planetary boundary layer scheme and microphysical scheme used in the present model for temperature simulation [43].Overall, the WRF model effectively captures meteorological parameters such as T2 and WS10, providing relatively accurate input data for the CMAQ model.

CMAQ Model
The CMAQ model was employed to simulate the same time period as the WRF model.Figure 3 presents box plots of simulated and monitored PM2.5 values in Changchun for January 2017 and 2020.It is evident that the simulated PM2.5 values for both 2017 and 2020 are slightly lower than the observed values.Specifically, the simulated average PM2.5 value in January 2017 was 71.95 µg/m 3 , whereas the observed average value was 89.23 µg/m 3 .In January 2020, the simulated average PM2.5 value was 91.14 µg/m 3 , while the observed average value was 115.10 µg/m 3 .Notably, the simulated average PM2.5 value in January 2020 was 26.7% higher than that in January 2017, while the observed average value was 28.9% higher.The difference between the two is relatively small.
To evaluate the simulation results of the CMAQ model, PM2.5 hourly monitoring data from five national air quality monitoring stations located in the built-up area of Changchun were selected for validation, as shown in Figure 4. Table 5 presents the statistical results of PM2.5 concentration simulation compared to monitoring values.It is evident that the simulation performance for January 2017 is significantly better than that for January 2020, with the former outperforming the latter in all evaluation indicators.Overall, the NMB values ranged from 0.042 to 0.41, and the NME values ranged from 0.48 to 0.61 for all sites in both years, with an average NMB value of 0.26 and an average NME value of 0.59.The higher NME may be attributed to uncertainties in the parameter program and In 2020, T2 followed a similar trend as observed in 2017, displaying lower values than the observed ones.This underestimation of T2 by the model is consistent with findings from other studies, suggesting a possible correlation with the physicochemical scheme employed in the model [42,43].Overall, the R of T2, WS10, and WD10 are all above 0.5, indicating a strong correlation.The R of WS10 and WD10 are lower than those of T2, which is consistent with the results of other scholars, suggesting relatively poorer simulation accuracy for wind speed and direction compared to T2 due to the subsurface effects on the surface wind field.In terms of NMB, NME, MFB, MFE, and RMSE, January 2017 is simulated better than January 2020 for all three meteorological parameters.This difference in simulation quality between the two time periods may be attributed to frequent and drastic meteorological changes in July 2020.Notably, the NMB for T2 in 2020 reaches 74.22%, which may be related to limitations in the planetary boundary layer scheme and microphysical scheme used in the present model for temperature simulation [44].Overall, the WRF model effectively captures meteorological parameters such as T2 and WS10, providing relatively accurate input data for the CMAQ model.

CMAQ Model
The CMAQ model was employed to simulate the same time period as the WRF model.Figure 3 presents box plots of simulated and monitored PM 2.5 values in Changchun for January 2017 and 2020.It is evident that the simulated PM 2.5 values for both 2017 and 2020 are slightly lower than the observed values.Specifically, the simulated average PM 2.5 value in January 2017 was 71.95 µg/m 3 , whereas the observed average value was 89.23 µg/m 3 .In January 2020, the simulated average PM 2.5 value was 91.14 µg/m 3 , while the observed average value was 115.10 µg/m 3 .Notably, the simulated average PM 2.5 value in January 2020 was 26.7% higher than that in January 2017, while the observed average value was 28.9% higher.The difference between the two is relatively small.
To evaluate the simulation results of the CMAQ model, PM 2.5 hourly monitoring data from five national air quality monitoring stations located in the built-up area of Changchun were selected for validation, as shown in Figure 4. Table 5 presents the statistical results of PM 2.5 concentration simulation compared to monitoring values.It is evident that the simulation performance for January 2017 is significantly better than that for January 2020, with the former outperforming the latter in all evaluation indicators.Overall, the NMB values ranged from 0.042 to 0.41, and the NME values ranged from 0.48 to 0.61 for all sites in both years, with an average NMB value of 0.26 and an average NME value of 0.59.The higher NME may be attributed to uncertainties in the parameter program and the emission inventories [45,46].Furthermore, the MFB and MFE of all sites satisfy the criteria of −60% ≤ MFB ≤ 60%, and MFE ≤ 75%.In summary, the validation results of the CMAQ model are acceptable for subsequent analysis.
the emission inventories [44,45].Furthermore, the MFB and MFE of all sites satisfy the criteria of −60%≤ MFB ≤ 60%, and MFE ≤ 75%.In summary, the validation results of the CMAQ model are acceptable for subsequent analysis.Observing Figure 4, it is evident that the CMAQ model effectively captures the peak variations of pollutants, and the observed and simulated values exhibit similar trends.Notably, the model performs well in simulating heavy pollution episodes from 15-17 January 2017.However, it is observed that the IPT, LP, and GMO sites tend to overestimate PM2.5 pollution levels.Overall, the simulated values for each site tend to be lower, result-  Observing Figure 4, it is evident that the CMAQ model effectively captures the peak variations of pollutants, and the observed and simulated values exhibit similar trends.Notably, the model performs well in simulating heavy pollution episodes from 15-17 January 2017.However, it is observed that the IPT, LP, and GMO sites tend to overestimate PM 2.5 pollution levels.Overall, the simulated values for each site tend to be lower, resulting in an underestimation of PM 2.5 pollution from 6-8 January 2017 and 24-28 January 2020.This discrepancy may stem from deviations between the inventory building results and actual source emissions, particularly in terms of spatial and temporal distribution.Additionally, some bias may be present in the topographic resolution and the WRF-CMAQ model's representation of complex topographic meteorological fields and pollutant concentration distribution fields.Therefore, in alignment with previous findings, it can be concluded that meteorological conditions in January 2020 were more unfavorable for pollutant dispersion compared to January 2017, resulting in higher PM 2.5 concentrations [32,33].

The Influence of Meteorological Parameters on PM 2.5
To investigate the influence of meteorological parameters on PM 2.5 concentrations, this paper compares the differences between 2017 and 2020 based on meteorological parameters such as T2, WS10, air pressure, and PBLH. Figure 5a depicts the temperature trends for January 2017 and 2020.The average temperature in January 2017 was −16.7 • C, whereas the average temperature in January 2020 was −17.8 • C, indicating a 1.1 • C lower temperature in 2020 compared to 2017. Figure 5b,c illustrate that the overall wind speed is lower in 2020 than in 2017, with the average wind speed in January 2017 being 4.04 m/s and the average wind speed in January 2020 being 3.47 m/s.Notably, during the period of severe pollution from 24-28 January 2020, wind speeds were significantly lower compared to the same period in 2017.Lower wind speeds are more conducive to the accumulation of pollutants, which raises concerns regarding air quality during this period.
can be concluded that meteorological conditions in January 2020 were more unfavorable for pollutant dispersion compared to January 2017, resulting in higher PM2.5 concentrations [30,31].

The Influence of Meteorological Parameters on PM2.5
To investigate the influence of meteorological parameters on PM2.5 concentrations, this paper compares the differences between 2017 and 2020 based on meteorological parameters such as T2, WS10, air pressure, and PBLH. Figure 5a depicts the temperature trends for January 2017 and 2020.The average temperature in January 2017 was −16.7 °C, whereas the average temperature in January 2020 was −17.8 °C, indicating a 1.1 °C lower temperature in 2020 compared to 2017. Figure 5b,c illustrate that the overall wind speed is lower in 2020 than in 2017, with the average wind speed in January 2017 being 4.04 m/s and the average wind speed in January 2020 being 3.47 m/s.Notably, during the period of severe pollution from 24-28 January 2020, wind speeds were significantly lower compared to the same period in 2017.Lower wind speeds are more conducive to the accumulation of pollutants, which raises concerns regarding air quality during this period.The barometric pressure values depicted in Figure 6a reveal that the difference between 2017 and 2020 is not significant, but there are higher barometric pressure values during periods of heavy PM 2.5 pollution, such as 6-8 January and 15-18 January 2017, and 24-28 January 2020.During periods of high pressure, low wind speeds accompany the sinking movement of air, which tends to create a sinking inversion and hinder the upward diffusion of pollutants.Consequently, PM 2.5 pollution is more likely to occur.In Figure 6b, the PBLH for January 2017 and 2020 are compared.The monthly average value of PBLH in January 2017 was 259 m, whereas the monthly average value of PBLH in January 2020 was 138 m, representing an 87.7% decrease compared to 2017.Notably, during the pollution period of 24-28 January 2020, the PBLH was significantly lower than in 2017.There is a significant negative correlation between PBLH and PM 2.5 pollution.From the above analysis, it is evident that the meteorological parameters contributing to higher PM 2.5 pollution in January 2020 compared to January 2017 include lower wind speeds, higher air pressure, and lower boundary layer heights.These conditions favor the accumulation of PM 2.5 in the atmosphere and hinder its dispersion, thus exacerbating pollution levels.Therefore, it is imperative to consider the influence of meteorological parameters when developing strategies to mitigate PM 2.5 pollution.
The barometric pressure values depicted in Figure 6a reveal that the difference between 2017 and 2020 is not significant, but there are higher barometric pressure values during periods of heavy PM2.5 pollution, such as 6-8 January and 15-18 January 2017, and 24-28 January 2020.During periods of high pressure, low wind speeds accompany the sinking movement of air, which tends to create a sinking inversion and hinder the upward diffusion of pollutants.Consequently, PM2.5 pollution is more likely to occur.In Figure 6b, the PBLH for January 2017 and 2020 are compared.The monthly average value of PBLH in January 2017 was 259 m, whereas the monthly average value of PBLH in January 2020 was 138 m, representing an 87.7% decrease compared to 2017.Notably, during the pollution period of 24-28 January 2020, the PBLH was significantly lower than in 2017.There is a significant negative correlation between PBLH and PM2.5 pollution.From the above analysis, it is evident that the meteorological parameters contributing to higher PM2.5 pollution in January 2020 compared to January 2017 include lower wind speeds, higher air pressure, and lower boundary layer heights.These conditions favor the accumulation of PM2.5 in the atmosphere and hinder its dispersion, thus exacerbating pollution levels.Therefore, it is imperative to consider the influence of meteorological parameters when developing strategies to mitigate PM2.5 pollution.
significant reduction in all six pollutants (PM2.5, PM10, SO2, NO2, O3, CO) by 8.39% (3.87 µg/m 3 ), 24.68% (19.94 µg/m 3 ), 62.05% (15.57µg/m 3 ), 19.62% (7.45 µg/m 3 ), 10.38% (9.43µg/m 3 ), and 36.65%(0.41 mg/m 3 ).Among them, the reduction in SO2 is the most significant, while the reduction in PM2.5 is not satisfactory, which may be related to the severe meteorological conditions in 2020.The AQI quantifies air quality, and the level of pollution is positively correlated with the magnitude of the Air Quality Index (AQI) values.Figure 8 shows the percentage of days in 2017 and 2020 for each of the six AQI health categories.Compared to 2017, air quality at all ten sites in 2020 will improve significantly, with an average increase of 84.2 days for excellent (AQI < 50) and 29.9 days for excellent and good (AQI < 100).While the number of days for excellent pollution increases significantly, the number of days for serious and very severe pollution (AQI > 200) increases by an average of 4.2 days from 2017 to 2020.We then analyze the average AQI hourly concentrations in Changchun for 2017 and 2020, as shown in Figure 9.The four seasons in Changchun are divided into spring (4.20-6.30),summer (7.01-8.10),autumn (8.11-10.10),and winter (1.1-4.19 and 10.11-12.31)according to the climatic seasonal division method.It can be seen that the overall pattern of air quality in Changchun is more polluted in the winter and less polluted in the summer.In 2017, the monthly average AQI value was highest in January at 118.87, followed by October at 99.2.In addition, it is worth noting that the highest hourly AQI value of the year occurred in May, mainly because of a strong dust storm weather process in Changchun starting on 6 May 2017.The dusty weather came from the Inner Mongolia Plateau and affected Changchun from 4 May.AQI from 6-8 was in the range of 116 to 256, making Changchun have the highest AQI hourly value for the year in May.The hourly AQI values in 2020 remain low in summer, with the most severe pollution in January, when the monthly average AQI value is 152.77,28.5% higher year-on-year compared to 2017.The number of days in January 2017 when the daily average AQI reached moderately polluted and above (AQI > 150) was 9, and the number of days in January 2020 when the daily average AQI reached moderately polluted and above (AQI > 150) was 14, which suggests that pollution is significantly worse in January 2020 than in 2017.The monthly average AQI value in April 2020 is only lower than that in January, at 132.56, which is 48.2% higher than the same period in 2017, due to the impact of the COVID-19 epidemic, which led to an increase in emission sources as work resumed and production resumed The AQI quantifies air quality, and the level of pollution is positively correlated with the magnitude of the Air Quality Index (AQI) values.Figure 8 shows the percentage of days in 2017 and 2020 for each of the six AQI health categories.Compared to 2017, air quality at all ten sites in 2020 will improve significantly, with an average increase of 84.2 days for excellent (AQI < 50) and 29.9 days for excellent and good (AQI < 100).While the number of days for excellent pollution increases significantly, the number of days for serious and very severe pollution (AQI > 200) increases by an average of 4.2 days from 2017 to 2020.We then analyze the average AQI hourly concentrations in Changchun for 2017 and 2020, as shown in Figure 9.The four seasons in Changchun are divided into spring (4.20-6.30),summer (7.01-8.10),autumn (8.11-10.10),and winter (1.1-4.19 and 10.11-12.31)according to the climatic seasonal division method.It can be seen that the overall pattern of air quality in Changchun is more polluted in the winter and less polluted in the summer.In 2017, the monthly average AQI value was highest in January at 118.87, followed by October at 99.2.In addition, it is worth noting that the highest hourly AQI value of the year occurred in May, mainly because of a strong dust storm weather process in Changchun starting on 6 May 2017.The dusty weather came from the Inner Mongolia Plateau and affected Changchun from 4 May.AQI from 6-8 was in the range of 116 to 256, making Changchun have the highest AQI hourly value for the year in May.The hourly AQI values in 2020 remain low in summer, with the most severe pollution in January, when the monthly average AQI value is 152.77,28.5% higher year-on-year compared to 2017.The number of days in January 2017 when the daily average AQI reached moderately polluted and above (AQI > 150) was 9, and the number of days in January 2020 when the daily average AQI reached moderately polluted and above (AQI > 150) was 14, which suggests that pollution is significantly worse in January 2020 than in 2017.The monthly average AQI value in April 2020 is only lower than that in January, at 132.56, which is 48.2% higher than the same period in 2017, due to the impact of the COVID-19 epidemic, which led to an increase in emission sources as work resumed and production resumed in April in various places.And the spring of 2020 saw a large, high-intensity open straw burning in the northeast, which was the main cause of heavy pollution.This was accompanied by static meteorological conditions and the presence of inversions, which aggravated the accumulation of pollutants.Changchun Municipal Government issued an emergency yellow alert for heavy air pollution on 15 April, and the air quality improved on the 19th.burning in the northeast, which was the main cause of heavy pollution.This was accom panied by static meteorological conditions and the presence of inversions, which aggra vated the accumulation of pollutants.Changchun Municipal Government issued an emer gency yellow alert for heavy air pollution on 15 April, and the air quality improved on th 19th.The NHAQI used in this paper is more effective for representing health risks com pared to the AQI, and Figure 10 shows the distribution of days in each category for th NHAQI based on the health classification categories of the AQI, including 2017 and 2020 in April in various places.And the spring of 2020 saw a large, high-intensity open straw burning in the northeast, which was the main cause of heavy pollution.This was accompanied by static meteorological conditions and the presence of inversions, which aggravated the accumulation of pollutants.Changchun Municipal Government issued an emergency yellow alert for heavy air pollution on 15 April, and the air quality improved on the 19th.The NHAQI used in this paper is more effective for representing health risks compared to the AQI, and Figure 10 shows the distribution of days in each category for the NHAQI based on the health classification categories of the AQI, including 2017 and 2020.The NHAQI used in this paper is more effective for representing health risks compared to the AQI, and Figure 10 shows the distribution of days in each category for the NHAQI based on the health classification categories of the AQI, including 2017 and 2020.When the AQI is excellent or good (AQI < 100), the NHAQI is also excellent or good because there is no excess risk ER at this time.When the AQI is light pollution (100 < AQI < 150), the percentage of days when the NHAQI is light or moderate pollution is 92% and 8% in 2017 and 100% and 0% in 2020, respectively.When the AQI is moderate pollution (150 < AQI < 200), the percentage of days with moderate and serious pollution in the NHAQI is 59% and 41% in 2017 and 79% and 21% in 2020, respectively.When AQI is serious pollution (200 < AQI < 300), the percentage of days with NHAQI of serious and very severe pollution is 78% and 22% in 2017 and 58% and 41% in 2020, respectively.When AQI is very severe pollution (AQI > 300), NHAQI is also very severe pollution.Compared with 2017, the number of days with NHAQI of light and moderate pollution in 2020 has decreased and the air quality has improved, but the number of days with serious pollution and above has increased, which can bring an increased health risk to the exposed population.Overall, AQI is inadequate in characterizing the health risk from air pollution, and there is an underestimation.
When the AQI is excellent or good (AQI < 100), the NHAQI is also excellent or good be cause there is no excess risk ER at this time.When the AQI is light pollution (100 < AQI < 150), the percentage of days when the NHAQI is light or moderate pollution is 92% and 8% in 2017 and 100% and 0% in 2020, respectively.When the AQI is moderate pollution (150 < AQI < 200), the percentage of days with moderate and serious pollution in the NHAQI is 59% and 41% in 2017 and 79% and 21% in 2020, respectively.When AQI i serious pollution (200 < AQI < 300), the percentage of days with NHAQI of serious and very severe pollution is 78% and 22% in 2017 and 58% and 41% in 2020, respectively.When AQI is very severe pollution (AQI > 300), NHAQI is also very severe pollution.Compared with 2017, the number of days with NHAQI of light and moderate pollution in 2020 ha decreased and the air quality has improved, but the number of days with serious pollution and above has increased, which can bring an increased health risk to the exposed popula tion.Overall, AQI is inadequate in characterizing the health risk from air pollution, and there is an underestimation.The NHAQI in the spring of both years is good (50 < NHAQI < 100), with a 29.43% decrease in 2020 compared to 2017, the largest decrease among the four seasons, but high NHAQI is observed in all areas on 6 May 2017 due to dust storms, and in April 2020 due to straw burning in the suburbs of Changchun, from 5 April lasted until 18 April, during which the NHAQI is serious pollution and above.The NHAQI values in summer and autumn are excellent or good (NHAQI < 100), decreasing by 18.20% and 18.14%, respec tively, with LY, NG, and ED having relatively higher NHAQI in summer and KC, LY, CY and ED having relatively higher NHAQI in autumn, which is mainly due to the short and relatively stable meteorology in summer and autumn in Changchun, where the air quality The NHAQI in the spring of both years is good (50 < NHAQI < 100), with a 29.43% decrease in 2020 compared to 2017, the largest decrease among the four seasons, but high NHAQI is observed in all areas on 6 May 2017 due to dust storms, and in April 2020 due to straw burning in the suburbs of Changchun, from 5 April lasted until 18 April, during which the NHAQI is serious pollution and above.The NHAQI values in summer and autumn are excellent or good (NHAQI < 100), decreasing by 18.20% and 18.14%, respectively, with LY, NG, and ED having relatively higher NHAQI in summer and KC, LY, CY, and ED having relatively higher NHAQI in autumn, which is mainly due to the short and relatively stable meteorology in summer and autumn in Changchun, where the air quality index is mainly influenced by anthropogenic emissions.And these five areas are all located in the main urban area, where anthropogenic emissions are more intensive.In the long winter of 2017, NHAQI values were above light pollution (NHAQI > 100) except for NG, which can pose a risk to population health, probably due to the need for coal burning for thermal heating in winter and frequent atmospheric stagnation in meteorological conditions.In 2020, although it decreased by 8.241% compared to 2017, the NHAQI in KC still showed light pollution, which may be related to the prevailing northwesterly winds in winter in Changchun.In general, NG shows low NHAQI in any season due to the presence of National Forests in the area and the large area of green space, in addition to the fact that pollution is more severe in spring and winter than in summer and fall, especially in densely populated areas.
index is mainly influenced by anthropogenic emissions.And these five areas are all located in the main urban area, where anthropogenic emissions are more intensive.In the long winter of 2017, NHAQI values were above light pollution (NHAQI > 100) except for NG, which can pose a risk to population health, probably due to the need for coal burning for thermal heating in winter and frequent atmospheric stagnation in meteorological conditions.In 2020, although it decreased by 8.241% compared to 2017, the NHAQI in KC still showed light pollution, which may be related to the prevailing northwesterly winds in winter in Changchun.In general, NG shows low NHAQI in any season due to the presence of National Forests in the area and the large area of green space, in addition to the fact that pollution is more severe in spring and winter than in summer and fall, especially in densely populated areas.Figure 12 shows the ERtotal of six administrative regions in Changchun for 2017 and 2020.the average ERtotal for 2017 is 0.65%, and among the six administrative regions, NG has the smallest ER value of 0.40%, and KC, CY, and LY have higher than average ER values of 0.84%, 0.72%, and 0.67%, respectively.The pollutants that contribute most to the ERtotal are PM2.5, PM10, and O3, with 45.46%, 33.30%, and 13.57%, respectively, and CO does not cause excess health risk.The average ERtotal in 2020 decreased by 0.11% compared to 2017, and all six regions decreased significantly except NG, which remained unchanged, and KC and CY, which had above-average ER values of 0.63% and 0.58%, respectively.The pollutants contributing most to ERtotal are still PM2.5, PM10, and O3, with 64.71%, 22.32%, and 8.41%, respectively, which shows that the control of PM10 and O3 is very effective, but the control of PM2.5 has weakened.Figure 12 shows the ER total of six administrative regions in Changchun for 2017 and 2020.the average ER total for 2017 is 0.65%, and among the six administrative regions, NG has the smallest ER value of 0.40%, and KC, CY, and LY have higher than average ER values of 0.84%, 0.72%, and 0.67%, respectively.The pollutants that contribute most to the ER total are PM 2.5 , PM 10, and O 3, with 45.46%, 33.30%, and 13.57%, respectively, and CO does not cause excess health risk.The average ER total in 2020 decreased by 0.11% compared to 2017, and all six regions decreased significantly except NG, which remained unchanged, and KC and CY, which had above-average ER values of 0.63% and 0.58%, respectively.The pollutants contributing most to ER total are still PM 2.5 , PM 10, and O 3, with 64.71%, 22.32%, and 8.41%, respectively, which shows that the control of PM 10 and O 3 is very effective, but the control of PM 2.5 has weakened.

Conclusions
This paper investigates the effect of meteorological conditions on PM2.5 through the simulation results of the WRF-CMAQ model and assesses the exposure risk of pollutants to the population using the NHAQI index and ER.The results show that the higher PM2.5 in January 2020 is mainly due to the influence of meteorological conditions.By studying the effect of meteorological parameters on PM2.5, we found that T2, WS10, and PBLH were all lower in January 2020 than in January 2017, while the air pressure was slightly higher than in January 2017.During the heavy PM2.5 pollution period from 24-28 January 2020, T2, WS10, and PBLH were negatively correlated with PM2.5 concentrations, while air pressure was positively correlated.The lower wind speed and larger air pressure, as well as the lower boundary layer height, are factors that favor the accumulation of PM2.5 in the atmosphere and are not easily diffused, causing pollution.The average number of days with excellent and good air quality (AQI < 100) in 2020 is 29.9 days more than the annual average in 2017, but AQI is lacking in evaluating health risks, so NHAQI is introduced.The NHAQI of all six municipal districts in Changchun decreased significantly from 2017 to 2020, but the four districts with a relatively high population (i.e., KC, LY, CY, and ED) had relatively high NHAQI values and the highest NHAQI in winter among the four seasons.The pollutants that contribute most to the ERtotal shift are PM2.5, PM10, and NO2 in 2017 and PM2.5 and PM10 in 2020.This paper will contribute to the future pollutant control strategy of Changchun City.
However, there are some issues in this study that still lead to the creation of uncertainty.First, due to the limitations of the inventory data, the 2017 MEIC inventory data is used as the inventory input for case 2, which leads to errors in the subsequent simulations, and this study will continue to research updated inventory files to improve them.Secondly, this study do not take into account the effect of meteorological conditions on PM2.5 without considering the effect of long-distance transmission on the local area, which may lead to uncertainty in the results and still requires further research.Thirdly, the national exposure-response relationship coefficients are used in this paper to represent Changchun City, and we will further refine the coefficients in the future by combining epidemiologic studies and clinical data.Finally, the calculation of NHAQI and ER may generate statistical errors, and further research is needed for NHAQI to account for the interactions between pollutants.

Conclusions
This paper investigates the effect of meteorological conditions on PM 2.5 through the simulation results of the WRF-CMAQ model and assesses the exposure risk of pollutants to the population using the NHAQI index and ER.The results show that the higher PM 2.5 in January 2020 is mainly due to the influence of meteorological conditions.By studying the effect of meteorological parameters on PM 2.5 , we found that T2, WS10, and PBLH were all lower in January 2020 than in January 2017, while the air pressure was slightly higher than in January 2017.During the heavy PM 2.5 pollution period from 24-28 January 2020, T2, WS10, and PBLH were negatively correlated with PM 2.5 concentrations, while air pressure was positively correlated.The lower wind speed and larger air pressure, as well as the lower boundary layer height, are factors that favor the accumulation of PM 2.5 in the atmosphere and are not easily diffused, causing pollution.The average number of days with excellent and good air quality (AQI < 100) in 2020 is 29.9 days more than the annual average in 2017, but AQI is lacking in evaluating health risks, so NHAQI is introduced.The NHAQI of all six municipal districts in Changchun decreased significantly from 2017 to 2020, but the four districts with a relatively high population (i.e., KC, LY, CY, and ED) had relatively high NHAQI values and the highest NHAQI in winter among the four seasons.The pollutants that contribute most to the ER total shift are PM 2.5 , PM 10, and NO 2 in 2017 and PM 2.5 and PM 10 in 2020.This paper will contribute to the future pollutant control strategy of Changchun City.
However, there are some issues in this study that still lead to the creation of uncertainty.First, due to the limitations of the inventory data, the 2017 MEIC inventory data is used as the inventory input for case 2, which leads to errors in the subsequent simulations, and this study will continue to research updated inventory files to improve them.Secondly, this study do not take into account the effect of meteorological conditions on PM 2.5 without considering the effect of long-distance transmission on the local area, which may lead to uncertainty in the results and still requires further research.Thirdly, the national exposureresponse relationship coefficients are used in this paper to represent Changchun City, and we will further refine the coefficients in the future by combining epidemiologic studies and clinical data.Finally, the calculation of NHAQI and ER may generate statistical errors, and further research is needed for NHAQI to account for the interactions between pollutants.

Figure 2 .
Figure 2. Simulation results of WRF meteorological parameters in Changchun in 2017 and 2020 (2017 on the left, 2020 on the right).

Figure 2 .
Figure 2. Simulation results of WRF meteorological parameters in Changchun in 2017 and 2020 (2017 on the left, 2020 on the right).

Figure 4 .
Figure 4. CMAQ simulation results for each monitoring station in Changchun City in January 2017 and 2020 (2017 on the left, 2020 on the right).

Figure 4 .
Figure 4. CMAQ simulation results for each monitoring station in Changchun City in January 2017 and 2020 (2017 on the left, 2020 on the right).

Figure 6 .
Figure 6.AP (a) and PBLH (b) in January from 2017 to 2020 in Changchun.

Figure 7
Figure 7 shows the comparison of the annual average mass concentrations of six pollutants in Changchun in 2017 and 2020 (Data from 10 state-controlled atmospheric environment automatic monitoring stations in Changchun).From 2017 to 2020, there was a

Figure 6 .
Figure 6.AP (a) and PBLH (b) in January from 2017 to 2020 in Changchun.

Figure 7 .
Figure 7. Annual average mass concentrations of six pollutants in Changchun in 2017 and 2020 (The left y-axis is applicable to PM2.5, PM10, SO2, NO2, and O3, and the right y-axis is applicable to CO).

Figure 7 .
Figure 7. Annual average mass concentrations of six pollutants in Changchun in 2017 and 2020 (The left y-axis is applicable to PM 2.5 , PM 10 , SO 2 , NO 2 , and O 3 , and the right y-axis is applicable to CO).

Figure 8 .
Figure 8. Percentage of days in six AQI health categories in Changchun in 2017 and 2020.

Figure 9 .
Figure 9. Hourly change in air quality index in Changchun in 2017 and 2020.

Figure 8 .
Figure 8. Percentage of days in six AQI health categories in Changchun in 2017 and 2020.

Figure 8 .
Figure 8. Percentage of days in six AQI health categories in Changchun in 2017 and 2020.

Figure 9 .
Figure 9. Hourly change in air quality index in Changchun in 2017 and 2020.

Figure 9 .
Figure 9. Hourly change in air quality index in Changchun in 2017 and 2020.

Figure 10 .
Figure 10.Distribution of NHAQI days in each category for 2017 and 2020 based on AQI health category determination.

Figure 10 .
Figure 10.Distribution of NHAQI days in each category for 2017 and 2020 based on AQI health category determination.The spatial distribution of the seasonal average NHAQI is studied based on six administrative districts of Changchun city (KC: Kuancheng district, LY: Lvyuan district, CY: Chaoyang district, NG: Nanguan district, ED: Erdao district, and SY: Shuangyang district), as shown in Figure 11.All four seasons show a significant improvement in air quality from 2017 to 2020, with the annual average NHAQI decreasing from 80.87 in 2017 to 66.27 in 2020.The NHAQI in the spring of both years is good (50 < NHAQI < 100), with a 29.43% decrease in 2020 compared to 2017, the largest decrease among the four seasons, but high NHAQI is observed in all areas on 6 May 2017 due to dust storms, and in April 2020 due to straw burning in the suburbs of Changchun, from 5 April lasted until 18 April, during which the NHAQI is serious pollution and above.The NHAQI values in summer and autumn are excellent or good (NHAQI < 100), decreasing by 18.20% and 18.14%, respectively, with LY, NG, and ED having relatively higher NHAQI in summer and KC, LY, CY, and ED having relatively higher NHAQI in autumn, which is mainly due to the short and relatively stable meteorology in summer and autumn in Changchun, where the air quality index is mainly influenced by anthropogenic emissions.And these five areas are all located in the main urban area, where anthropogenic emissions are more intensive.In the long winter of 2017, NHAQI values were above light pollution (NHAQI > 100) except for NG, which can pose a risk to population health, probably due to the need for coal burning for thermal heating in

Figure 12 .
Figure 12.Average Total Excess Health Risk Values for Changchun City in 2017 and 2020.

Figure 12 .
Figure 12.Average Total Excess Health Risk Values for Changchun City in 2017 and 2020.

Table 1 .
List of the 10 individual air quality monitoring sites and the meteorological station.

Table 2 .
Ranges of AQI values and the corresponding pollutant concentration limits and health categories.

Table 3 .
Different scenarios were considered in this study.

Table 4 .
Statistical metrics for WRF model evaluation over Changchun in 2017 and 2020.

Table 4 .
Statistical metrics for WRF model evaluation over Changchun in 2017 and 2020.

Table 5 .
Error statistics of simulated PM2.5 concentrations and monitoring values in Changchun.

Table 5 .
Error statistics of simulated PM 2.5 concentrations and monitoring values in Changchun.