Spatiotemporal variations of air pollutants and ozone prediction using machine learning algorithms in the Beijing-Tianjin-Hebei region from 2014 to 2021

Understanding the variability of ground-level ozone and fine particulate matter over the Tibetan plateau with data-driven approach

The Tibetan Plateau, known as the "Third Pole", is susceptible to ground-level ozone (O3) and fine particulate matter (PM2.5) pollution due to its unique high-altitude environment. This study constructed random forest regression models using multi-source data from ground measurements and meteorological satellites to predict variations in ground-level O3 and PM2.5 concentrations and their influencing factors across seven major cities in the Tibetan Plateau over two-year periods. The models successfully reproduced O3 and PM2.5 levels with satisfactory R-squared values of 0.71 and 0.73, respectively. Results reveal combustion-related carbon monoxide (CO) and nitrogen dioxide (NO2) as the most substantial influences on O3 and PM2.5 concentrations. Solar radiation, geographical factors, and meteorological variables also played crucial roles in driving pollutant variations. Conversely, transport-related and human activity factors exhibited relatively lower significance. High O3 and PM2.5 pollution occurred during pre-monsoon and post-monsoon/winter seasons, driven by solar radiation and emissions, respectively. While CO consistently contributed across cities and seasons, key influencing factors varied locally. This study unveils the key driving forces governing air pollutant variations across the Tibetan Plateau, shedding light on complex atmospheric processes in this unique high-altitude region.

Enhanced PM2.5 estimation across China: An AOD-independent two-stage approach incorporating improved spatiotemporal heterogeneity representations

In China, population growth and aging have partially negated the public health benefits of air pollution control measures, underscoring the ongoing need for precise PM2.5 monitoring and mapping. Despite its prevalence, the satellite-derived Aerosol Optical Depth (AOD) method for estimating PM2.5 concentrations often encounters significant spatial data gaps. Additionally, current research still needs better representation of PM2.5 spatiotemporal heterogeneity. Addressing these challenges, we developed a two-stage model employing the Extreme Gradient Boosting (XGBoost) algorithm. By incorporating improved spatiotemporal factors, we achieved high-precision and full-coverage daily 1-km PM2.5 mappings across China for the year 2020 without utilizing AOD products. Specifically, Model 1 develops improved temporal encodings and a terrain classification factor (DC), while Model 2 constructs an enhanced spatial autocorrelation term (Ps) by integrating observed and estimated values. Notably, Model 2 excelled in 10-fold sample-based cross-validation, achieving a coefficient of determination of 0.948, a mean absolute error of 3.792 μg/m³, a root mean square error of 7.144 μg/m³, and a mean relative error of 14.171%. Feature importance and Shapley Additive exPlanations (SHAP) analyses determined the relative importance of predictors in model training and outcome prediction, while correlation analysis identified strong links between improved temporal encodings, PM2.5 concentrations, and significant meteorological factors. Two-way Partial Dependence Plots (PDPs) further explored the interactions among these factors and their impact on PM2.5 levels. Compared to traditional methods, improved temporal encodings align more closely with seasonal variations and synergize more effectively with meteorological factors. Besides, the structured nature of DC aids in model training, while the improved Ps more effectively captures PM2.5's spatial autocorrelation, outperforming traditional Ps. Overall, this study effectively represents spatiotemporal information, thereby boosting model accuracy and enabling seamless large-scale PM2.5 estimations. It provides deep insights into variables and models, providing significant implications for future air pollution research.

Spatio-temporal evolution mechanism and dynamic simulation of nitrogen and phosphorus pollution of the Yangtze River economic Belt in China

Excess nitrogen and phosphorus inputs are the main causes of aquatic environmental deterioration. Accurately quantifying and dynamically assessing the regional nitrogen and phosphorus pollution emission (NPPE) loads and influencing factors is crucial for local authorities to implement and formulate refined pollution reduction management strategies. In this study, we constructed a methodological framework for evaluating the spatio-temporal evolution mechanism and dynamic simulation of NPPE. We investigated the spatio-temporal evolution mechanism and influencing factors of NPPE in the Yangtze River Economic Belt (YREB) of China through the pollution load accounting model, spatial correlation analysis model, geographical detector model, back propagation neural network model, and trend analysis model. The results show that the NPPE inputs in the YREB exhibit a general trend of first rising and then falling, with uneven development among various cities in each province. Nonpoint sources are the largest source of land-based NPPE. Overall, positive spatial clustering of NPPE is observed in the cities of the YREB, and there is a certain enhancement in clustering. The GDP of the primary industry and cultivated area are important human activity factors affecting the spatial distribution of NPPE, with economic factors exerting the greatest influence on the NPPE. In the future, the change in NPPE in the YREB at the provincial level is slight, while the nitrogen pollution emissions at the municipal level will develop towards a polarization trend. Most cities in the middle and lower reaches of the YREB in 2035 will exhibit medium to high emissions. This study provides a scientific basis for the control of regional NPPE, and it is necessary to strengthen cooperation and coordination among cities in the future, jointly improve the nitrogen and phosphorus pollution tracing and control management system, and achieve regional sustainable development.

Elucidating contributions of volatile organic compounds to ozone formation using random forest during COVID-19 pandemic: A case study in China

Ozone has been reported to increase despite nitrogen oxides reductions during the COVID-19 pandemic, and ozone formation needs to be revisited using volatile organic compounds (VOCs), which are rarely measured during the pandemic. Here, a total of 98 VOCs species were monitored in an economy-active city in China from January 2021 to August 2022 to assess contributions to ozone formation during the pandemic. Total VOCs concentrations were 35.55 ± 21.47 ppb during the entire period, among which alkanes account for the largest fraction (13.78 ppb, 38.0%), followed by aromatics (6.16 ppb, 16.8%) and oxygenated VOCs (OVOCs, 5.69 ppb, 15.7%). Most VOCs groups (e.g., alkenes, OVOCs) and individual species (e.g., isoprene, methyl vinyl ketone) display obvious seasonal and diurnal variations, which are related to their sources and reactivities. No weekend effects of VOCs suggest limited influences from traffic emissions during pandemic. Aromatics and alkenes are the major contributors (39% and 33%) to ozone formation potential, largely driven by o/m/p-xylene (21%), ethylene (15%), toluene (9%). Secondary organic aerosol formation potential is dominated by toluene (>50%) despite its low proportion (5%). Further inclusion of VOCs and meteorology in the Random Forest model shows good ozone prediction performance (R2 = 0.77-0.86, RMSE = 11.95-19.91 μg/m3, MAE = 8.89-14.58 μg/m3). VOCs and NO2 contribute >50% of total importance with the largest difference in importance ratio of VOCs/NO2 in the summer and winter, implying ozone formation regime may vary. No seasonal variations in importance of meteorology are observed, while importance of other variables (e.g., PM2.5) is highest in the summer. This work identifies critical VOCs groups and species for ozone formation during the pandemic, and demonstrates the feasibility of machine learning algorithms in elucidation of ozone formation mechanisms.

A machine learning-based ensemble model for estimating diurnal variations of nitrogen oxide concentrations in Taiwan

Air pollution is inextricable from human activity patterns. This is especially true for nitrogen oxide (NOx), a pollutant that exists naturally and also as a result of anthropogenic factors. Assessing exposure by considering diurnal variation is a challenge that has not been widely studied. Incorporating 27 years of data, we attempted to estimate diurnal variations in NOx across Taiwan. We developed a machine learning-based ensemble model that integrated hybrid kriging-LUR, machine-learning, and an ensemble learning approach. Hybrid kriging-LUR was performed to select the most influential predictors, and machine-learning algorithms were applied to improve model performance. The three best machine-learning algorithms were suited and reassessed to develop ensemble learning that was designed to improve model performance. Our ensemble model resulted in estimates of daytime, nighttime, and daily NOx with high explanatory powers (Adj-R2) of 0.93, 0.98, and 0.94, respectively. These explanatory powers increased from the initial model that used only hybrid kriging-LUR. Additionally, the results depicted the temporal variation of NOx, with concentrations higher during the daytime than the nighttime. Regarding spatial variation, the highest NOx concentrations were identified in northern and western Taiwan. Model evaluations confirmed the reliability of the models. This study could serve as a reference for regional planning supporting emission control for environmental and human health.

Prediction of air pollutant concentrations based on the long short-term memory neural network

In recent years, environmental problems caused by air pollutants have received increasing attention. Effective prediction of air pollutant concentrations is an important way to protect the public from harm. Recently, due to extreme climate and social development, the forest fire frequency has increased. During the biomass combustion process caused by forest fires, the content of particulate matter (PM) in the atmosphere increases significantly. However, most existing air pollutant concentration prediction methods do not consider the considerable impact of forest fires, and effective long-term prediction models have not been established to provide early warnings for harmful gases. Therefore, in this paper, we collected a daily air quality data set (aerodynamic diameter smaller than 2.5 µm, PM2.5) for Heilongjiang Province, China, from 2017 to 2023 and A novel Long Short-Term Memory (LSTM) model was proposed to effectively predict the situation of air pollutants. The model could automatically extract information of the effective time step from the historical data set and combine forest fire disturbance and climate data as auxiliary data to improve the model prediction ability. Moreover, we created artificial neural network (ANN) and permissive regression (support vector machine, SVR) models for comparative experiments. The results showed that the precision accuracy of the developed LSTM model is higher. Unlike the other models, the LSTM neural network model could effectively predict the concentration of air pollutants in long-term series. Regarding long-term observation missions (7 days), the proposed model performed well and stably, with R2 reaching over 88%.

Long-term spatiotemporal variations in surface NO2 for Beijing reconstructed from surface data and satellite retrievals

Remote sensing data from the Ozone Monitoring Instrument (OMI) and the TROPOspheric Monitoring Instrument (TROPOMI) play important roles in estimating surface nitrogen dioxide (NO2), but few studies have compared their differences for application in surface NO2 reconstruction. This study aims to explore the effectiveness of incorporating the tropospheric NO2 vertical column density (VCD) from OMI and TROPOMI (hereafter referred to as OMI and TROPOMI, respectively, for conciseness) for deriving surface NO2 and to apply the resulting data to revisit the spatiotemporal variations in surface NO2 for Beijing over the 2005-2020 period during which there were significant reductions in nitrogen oxide emissions. In the OMI versus TROPOMI performance comparison, the cross-validation R2 values were 0.73 and 0.72, respectively, at 1 km resolution and 0.69 for both at 100 m resolution. The comparisons between satellite data sources indicate that even though TROPOMI has a finer resolution it does not improve upon OMI for deriving surface NO2 at 1 km resolution, especially for analyzing long-term trends. In light of the comparison results, we used a hybrid approach based on machine learning to derive the spatiotemporal distribution of surface NO2 during 2005-2020 based on OMI. We had novel, independent passive sampling data collected weekly from July to September of 2008 for hindcasting validation and found a spatiotemporal R2 of 0.46 (RMSE = 7.0 ppb). Regarding the long-term trend of surface NO2, the level in 2008 was obviously lower than that in 2007 and 2009, as expected, which was attributed to pollution restrictions during the Olympic Games. The NO2 level started to steadily decline from 2015 and fell below 2008's level after 2017. Based on OMI, a long-term and fine-resolution surface NO2 dataset was developed for Beijing to support future environmental management questions and epidemiological research.

Contrasting effects of clean air actions on surface ozone concentrations in different regions over Beijing from May to September 2013-2020

Due to the nonlinear impacts of meteorology and precursors, the response of ozone (O3) trends to emission changes is very complex over different regions in megacity Beijing. Based on long-term in-situ observations at 35 air quality sites (four categories, i.e., urban, traffic, northern suburban and southern suburban sites) and satellite data, spatiotemporal variability of O3, gaseous precursors, and O3-VOCs-NOx sensitivity were explored through multiple metrics during the warm season from 2013 to 2020. Additionally, the contribution of meteorology and emissions to O3 was separated by a machine-learning-based de-weathered method. The annual averaged MDA8 O3 and O3 increased by 3.7 and 2.9 μg/m3/yr, respectively, with the highest at traffic sites and the lowest in northern suburb, and the rate of Ox (O3 + NO2) was 0.2 μg/m3/yr with the highest in southern suburb, although NO2 declined strongly and HCHO decreased slightly. However, the increment of O3 and Ox in the daytime exhibited decreasing trends to some extent. Additionally, NOx abatements weakened O3 loss through less NO titration, which drove narrowing differences in urban-suburban O3 and Ox. Due to larger decrease of NO2 in urban region and HCHO in northern suburb, the extent of VOCs-limited regime fluctuated over Beijing and northern suburb gradually shifted to transition or NOx-limited regime. Compared with the directly observed trends, the increasing rate of de-weathered O3 was lower, which was attributed to favorable meteorological conditions for O3 generation after 2017, especially in June (the most polluted month); whereas the de-weathered Ox declined except in southern suburb. Overall, clean air actions were effective in reducing the atmospheric oxidation capacity in urban and northern suburban regions, weakening local photochemical production over Beijing and suppressing O3 deterioration in northern suburb. Strengthening VOCs control and keeping NOx abatement, especially in June, will be vital to reverse O3 increase trend in Beijing.

Analysis of severe ozone-related human health and weather influence over China in 2019 based on a high-resolution dataset

Ozone pollution in 2019 in China is particularly severe posing a tremendous threat to the health of Chinese inhabitants. In this study, we constructed a more reliable and accurate 1-km gridded dataset for 2019 with as many sites as possible using the inverse distance weight interpolation method to analyze spatiotemporal ozone pollution characteristics and health burden attributed to ozone exposure from the perspective of different diseases and weather influence. The accuracy of this new dataset is higher than other public datasets, with the coefficient of determination of 0.84 and root-mean-square error of 8.77 ppb through the validation of 300 external sites which were never used for establishing retrieval methods by the datasets mentioned-above. The averaged MDA8 (the daily maximum 8 h average) ozone concentrations over China was 43.5 ppb, and during April-July, 83.9% of total grids occurred peak-month ozone concentrations. Overall, the highest averaged exceedance days (60 days) and population-weighted ozone concentrations (55.0 ppb) both concentrated in central-eastern China including 9 provinces (only 11.4% of the national territory); meanwhile, all-cause premature deaths attributable to ozone exposure reached up to 142,000 (54.9% of national total deaths) with higher deaths for cardiovascular and respiratory, and the provincial per capita premature mortality was 0.27~0.44‰. The six most polluted weather types in the central-eastern China are in order as follows: westerly (SW and W), cyclonic, northerly, and southerly (NW, N, and S) types, which accounts for approximately 73.2% of health burden attributed to daily ozone exposure and poses the greatest public health risk with mean daily premature deaths ranging from 466 to 610. Our findings could provide an effective support for regional ozone pollution control and public health management in China.

Assessing uncertainty and heterogeneity in machine learning-based spatiotemporal ozone prediction in Beijing-Tianjin- Hebei region in China

Accurate prediction of spatiotemporal ozone concentration is of great significance to effectively establish advanced early warning systems and regulate air pollution control. However, the comprehensive assessment of uncertainty and heterogeneity in spatiotemporal ozone prediction remains unknown. Here, we systematically analyze the hourly and daily spatiotemporal predictive performances using convolutional long short term memory (ConvLSTM) and deep convolutional generative adversarial network (DCGAN) models over the Beijing-Tianjin-Hebei region in China from 2013 to 2018. In extensive scenarios, our results show that the machine learning-based (ML-based) models achieve better spatiotemporal ozone concentration prediction performance with multiple meteorological conditions. A further comparison to the air pollution model-Nested Air Quality Prediction Modelling System (NAQPMS) and monitoring observations, the ConvLSTM model demonstrates the practical feasibility of identifying high ozone concentration distribution and capturing spatiotemporal ozone variation patterns at a high spatial resolution (here 15 km × 15 km).

VAR-tree model based spatio-temporal characterization and prediction of O3 concentration in China

Ozone (O₃) pollution in the atmosphere is getting worse in many cities. In order to improve the accuracy of O₃ prediction and obtain the spatial distribution of O₃ concentration over a continuous period of time, this paper proposes a VAR-XGBoost model based on Vector autoregression (VAR), Kriging method and XGBoost (Extreme Gradient Boosting). China is used as an example and its spatial distribution of O₃ is simulated. In this paper, the O₃ concentration data of the monitoring sites in China are obtained, and then a spatial prediction method of O₃ mass concentration based on the VAR-XGBoost model is established, and finnally its influencing factors are analyzed. This paper concludes that O₃ features the highest correlation with PM_2.5 and the lowest correlation with SO₂. Among the measurement factors, wind speed and temperature are the most important factors affecting O₃ pollution, which are positively correlated to O₃ pollution. In addition, precipitation is negatively correlated with 8-hour ozone concentration. In this paper, the performance of the VAR-XGBoost model is evaluated based on the ten-fold cross-validation method of sample, site and time, and a comparison with the results of XGBoost, CatBoost (categorical boosting), ExtraTrees, GBDT (gradient boosting decision tree), AdaBoost (adaptive boosting), RF (random forest), Decision tree, and LightGBM (light gradient boosting machine) models is conducted. The result shows that the prediction accuracy of the VAR-XGBoost model is better than other models. The seasonal and annual average R² reaches 0.94 (spring), 0.93 (summer), 0.92 (autumn), 0.93 (winter), and 0.95 (average from 2016 to 2021). The data show that the applicability of the VAR-XGBoost model in simulating the spatial distribution of O₃ concentrations in China performs well. The spatial distribution of O₃ concentrations in the Chinese region shows an obvious feature of high in the east and low in the west, and the spatial distribution is strongly influenced by topographical factors. The mean concentration is clearly low in winter and high in summer within a season. The results of this study can provide a scientific basis for the prevention and control of regional O₃ pollution in China, and can also provide new ideas for the acquisition of data on the spatial distribution of O₃ concentrations within cities.

Combined exposure of heat stress and ozone enhanced cognitive impairment via neuroinflammation and blood brain barrier disruption in male rats

BACKGROUND

Heat stress (HS) exposure has been linked to cognitive dysfunction. In reality, high temperature does not occur alone in environment, and ozone (O₃) and heatwaves usually co-exist in atmospheric environment. However, whether O₃ exposure exacerbates HS-induced cognitive impairment and the potential underlying mechanisms have not been explored experimentally. The aim of this study was to determine the co-effects and mechanisms of HS and O₃ on the cognitive dysfunction.

METHODS

48 Sprague Dawley male rats were randomly divided into 4 groups: control, HS, O₃ and HS plus O₃ (HO₃) groups. Rats in HS and HO₃ group were exposed to 40 °C every morning from 9:00 to 12:00 for 15 consecutive days. While rats in O₃ and HO₃ groups were exposed to 0.7 ppm O₃ the same day from 14:00 to 17:00 for 15 days. Cognitive performance was examined with Morris water maze test. Neurodegeneration, glial activation, neuroinflammation, blood brain barrier (BBB) disruption and apoptosis were evaluated by Western blot, Elisa, immunohistochemistry and immunofluorescence staining.

RESULTS

HS induced cognitive decline and neuronal damage in rats. Further studies showed that exposure of rats to HS could also induce glial activation, neuroinflammation and neuronal apoptosis in hippocampus, and decrease in the expressions of ZO-1, claudin-5 and occluding, indicative of BBB disruption. Impressively, the neuronal effects induced by HS, as depicted above, could be worsened by co-exposure to O₃ in rats.

CONCLUSIONS

Co-exposure to O₃ promotes HS-induced cognitive impairment in rats possibly through glial-mediated neuroinflammation and BBB disruption.

Tracking long-term population exposure risks to PM2.5 and ozone in urban agglomerations of China 2015-2021

China has experienced severe air pollution in the past decade, especially PM2.5 and emerging ozone pollution recently. In this study, we comprehensively analyzed long-term population exposure risks to PM2.5 and ozone in urban agglomerations of China during 2015-2021 regarding two-stage clean-air actions based on the Ministry of Ecology and the Environment (MEE) air monitoring network. Overall, the ratio of the population living in the regions exceeding the Chinese National Ambient Air Quality Standard (35 μg/m3) decreases by 29.9 % for PM2.5 from 2015 to 2021, driven by high proportions in the Middle Plain (MP, 42.3 %) and Lan-Xi (35.0 %) regions. However, this ratio almost remains unchanged for ozone and even increases by 1.5 % in the MP region. As expected, the improved air quality leads to 234.7 × 103 avoided premature mortality (ΔMort), mainly ascribed to the reduction in PM2.5 concentration. COVID-19 pandemic may influence the annual variation of PM2.5-related ΔMort as it affects the shape of the population exposure curve to become much steeper. Although all eleven urban agglomerations share stroke (43.6 %) and ischaemic heart disease (IHD, 30.1 %) as the two largest contributors to total ΔMort, cause-specific ΔMort is highly regional heterogeneous, in which ozone-related ΔMort is significantly higher (21 %) in the Tibet region than other urban agglomeration. Despite ozone-related ΔMort being one order of magnitude lower than PM2.5-related ΔMort from 2015 to 2021, ozone-related ΔMort is predicted to increase in major urban agglomerations initially along with a continuous decline for PM2.5-related ΔMort from 2020 to 2060, highlighting the importance of ozone control. Coordinated controls of PM2.5 and O3 are warranted for reducing health burdens in China during achieving carbon neutrality.

Changes in Air Quality and Drivers for the Heavy PM2.5 Pollution on the North China Plain Pre- to Post-COVID-19

Under the clean air action plans and the lockdown to constrain the coronavirus disease 2019 (COVID-19), the air quality improved significantly. However, fine particulate matter (PM_2.5) pollution still occurred on the North China Plain (NCP). This study analyzed the variations of PM_2.5, nitrogen dioxide (NO₂), sulfur dioxide (SO₂), carbon monoxide (CO), and ozone (O₃) during 2017-2021 on the northern (Beijing) and southern (Henan) edges of the NCP. Furthermore, the drivers for the PM_2.5 pollution episodes pre- to post-COVID-19 in Beijing and Henan were explored by combining air pollutant and meteorological datasets and the weighted potential source contribution function. Results showed air quality generally improved during 2017-2021, except for a slight rebound (3.6%) in NO₂ concentration in 2021 in Beijing. Notably, the O₃ concentration began to decrease significantly in 2020. The COVID-19 lockdown resulted in a sharp drop in the concentrations of PM_2.5, NO₂, SO₂, and CO in February of 2020, but PM_2.5 and CO in Beijing exhibited a delayed decrease in March. For Beijing, the PM_2.5 pollution was driven by the initial regional transport and later secondary formation under adverse meteorology. For Henan, the PM_2.5 pollution was driven by the primary emissions under the persistent high humidity and stable atmospheric conditions, superimposing small-scale regional transport. Low wind speed, shallow boundary layer, and high humidity are major drivers of heavy PM_2.5 pollution. These results provide an important reference for setting mitigation measures not only for the NCP but for the entire world.