Abstract
Surface ozone is an important air pollutant detrimental to human health and vegetation productivity, particularly in China. However, high resolution surface ozone concentration data is still lacking, largely hindering accurate assessment of associated environmental impacts. Here, we collected hourly ground ozone observations (over 6 million records), remote sensing products, meteorological data, and social-economic information, and applied recurrent neural networks to map hourly surface ozone data (HrSOD) at a 0.1° × 0.1° resolution across China during 2015–2020. The coefficient of determination (R2) values in sample-based, site-based, and by-year cross-validations were 0.72, 0.65 and 0.71, respectively, with the root mean square error (RMSE) values being 11.71 ppb (mean = 30.89 ppb), 12.81 ppb (mean = 30.96 ppb) and 11.14 ppb (mean = 31.26 ppb). Moreover, it exhibits high spatiotemporal consistency with ground-level observations at different time scales (diurnal, seasonal, annual), and at various spatial levels (individual sites and regional scales). Meanwhile, the HrSOD provides critical information for fine-resolution assessment of surface ozone impacts on environmental and human benefits.
Similar content being viewed by others
Background & Summary
Ozone (O3) is an important constituent of the atmosphere and is ubiquitously present in both the troposphere and the stratosphere. Stratospheric ozone protects life on Earth by absorbing harmful solar ultraviolet rays1,2,3. Tropospheric ozone is a major gaseous pollutant produced in a series of complex reactions between volatile organic compounds (VOCs) and nitrogen oxides (NOx) in the presence of sunlight4. Exposure to high-concentration surface ozone can cause severe impacts on human health, inducing high morbidity in respiratory, cardiopulmonary, and cardiovascular diseases5,6,7. Moreover, surface ozone of high concentrations could damage the leaf cell structure of plants and thus decrease natural vegetation productivity, crop yield and quality8,9,10,11.
In the past decades, the number of ozone pollution events has increased significantly, particularly in highly populated and developed regions12,13,14,15. Real-time surface ozone monitoring networks have been established on a regional basis around the world. But their coverage is still insufficient in both space and time, due to uneven distribution of monitoring sites and lack of mid- to long-term continuous records in the majority of the world10,16. In contrast, satellite remote sensing can monitor the spatial and temporal variability of ozone at regional to global scales. For instance, the Ozone Monitoring Instrument (OMI) on the Aura satellite, launched in 2004, provides global daily total column ozone retrievals. Nonetheless, satellite-based estimates of surface ozone concentrations are not available at high spatial and temporal resolutions17,18. Hence, various models have been developed to extrapolate site observations, refine satellite retrievals, or fuse them to generate long-term, high-quality surface ozone datasets19,20.
These models, according to their underlying principles, can be generally grouped into chemical transport models (CTMs), geostatistical models, and machine learning models. CTMs are physics-based, accounting for atmospheric chemical reactions, emission inventories, meteorological conditions and transport of atmospheric pollutants, but usually are prone to high uncertainties in emission inventories and model assumptions21,22,23. Geostatistical models, such as Kriging interpolation24, land-use regression (LUR), Bayesian maximum entropy25 (BME), and geographically weighted regression26 (GWR), estimate surface ozone by fitting its relationships with the influential factors. However, collinearity (the non-independence of predictor variables) in these geostatistical models usually makes them difficult to estimate accurately19,27. Machine learning models, such as neural network, random forest (RF) and extreme gradient boosting (XGBoost), are widely used due to their strong data-mining capabilities. Among them, deep learning algorithms utilize more precise hidden layer structures for data-driven prediction, resulting in higher prediction accuracy than traditional regression and neural network models28, and have been developing rapidly and show great potential for predicting atmospheric pollutions including surface ozone concentrations. For instance, Eslami et al.29 utilized a deep convolutional neural network (CNN) to predict hourly ozone concentrations in Seoul, South Korea in 2017. Cheng et al.30 used a hybrid deep learning model to explore the complex nonlinear relationships between meteorological factors and ozone concentrations and applied it to hourly and daily forecasts of ozone concentrations in China.
In recent years, surface ozone pollution in China has become increasingly serious, with frequent large-scale high ozone pollution events31,32,33. Since 2013, China has established a national ozone observation network10, utilizing which several gridded surface ozone products were generated34,35. Liu et al.19 utilized the XGBoost algorithm in combination with monitoring station data, concurrent ozone retrievals, aerosol reanalysis, meteorological parameters, and land use data to predict maximum daily average 8-hour ozone (MDA8) concentration across China from 2015 to 2020. At the daily level, the coefficient of determination (R2) values for cross validation (CV) were 0.61–0.78. Wang et al.33 used a space-time extremely randomized trees (STET) model, with solar radiation intensity and air temperature as the main predicting factors, combined with ground observation data, meteorological data, and emission inventory data, to simulate MDA8 data across China from 2013 to 2020, with R2 of 0.87 and the root mean square error (RMSE) of 17.10 µg m−3. However, some input variables, particularly those related to ozone precursor emission inventories, were found to contribute less significantly than originally anticipated20. Moreover, the predictions were mostly focused on daily ozone concentrations, such as MDA8. Although there have been some exceptional datasets of hourly surface ozone concentrations36, long-term gridded hourly products of high accuracy are still lacking in China. Such a data gap impedes accurate assessment of environmental and human health impacts of surface ozone. For example, in estimating ozone damage to crops, hourly ozone data is usually required for stomatal ozone flux models37 or generating ozone exposure index38,39. Moreover, hourly ozone data is advantageous over that at coarser temporal resolution in determining ozone exposure of humans40.
To address the issue, here we developed a deep learning model based on the Long Short-Term Memory (LSTM) recurrent neural networks to generate hourly surface ozone data (HrSOD) at a spatial resolution of 0.1° × 0.1° from 2015 to 2020 over China. The model utilized a parsimonious set of predictor variables (excluding co-linear variables and ozone precursor emission inventories), including meteorological factors, remote sensing data, socio-economic and land use data, and more than six million ground station monitoring records as references.
Methods
Data
Surface ozone observation data
Over six million records of hourly surface ozone concentration measurements during June 2014 to February 2021 were obtained from the real-time air quality monitoring platform of the China National Environmental Monitoring Centre (CNEMC; https://air.cnemc.cn:18007/) and the archived data was uploaded to the Zenodo repository41 (https://doi.org/10.5281/zenodo.10911197). The monitoring network was expanded to more than 1500 monitoring sites from 2013 to 2020, covering 31 provinces and 368 cities across mainland China. However, these monitoring sites are mainly located in the eastern region of China, with a much lower site distribution density in the northwest and the Qinghai-Tibet Plateau (Fig. 1).
Hourly ozone concentrations are measured at all monitoring sites by continuous monitoring instruments, and the unit of ozone reported by CNEMC is µg m−3 (standard atmospheric conditions at a temperature of 273.0 K and a pressure of 1013.25 hPa; 1 µg m−3 = 0.467 ppb). According to the Ambient Air Quality Standard42 and the Technical Specification for Ambient Air Quality Assessment43 set by the Ministry of Ecology and Environment of China (MEE) for ozone concentration data norms and standards, the ozone data was screened by removing outliers and null values. The multi-year mean hourly ozone concentrations ranged from 14–48 ppb during 2015–2020 in China, with areas of high ozone concentrations mainly in eastern China, especially in four densely populated megacity clusters of China, including the Beijing-Tianjin-Hebei (BTH) region, the Pearl River Delta (PRD), the Sichuan Basin (SCB) and the Yangtze River Delta (YRD).
Predictor variables
The predictor variables include satellite retrieved ozone products, meteorological factors, land use, population, and gross domestic product (see Supplementary Table S1).
-
Remote sensing data
The OMI carried by the Earth Observing System (EOS) Aura satellite was launched by the United States in 2004. Its primary mission is to monitor trace gases in the atmosphere, such as ozone, sulfur dioxide and nitrogen dioxide, while also collecting information on aerosols, clouds, ozone profiles, etc. The OMI sensor operates in a wavelength range of 270 to 500 nm with a spectral resolution of 0.5 nm. It has a swath width of 2600 km and provides a spatial resolution of 13 km × 24 km. OMI can complete a global scan in just one day, measuring column concentrations and profiles of O3, NO2, SO2, as well as data on aerosols, clouds, surface ultraviolet radiation, and various other parameters18,44. Previous studies have shown that the OMI ozone column concentrations and profile data exhibit a reasonable consistency with lower- to mid-troposphere ozone across the world17,45. Similarly, the OMI ozone data for different cities in China also manifests a high consistency with ground measurements46,47, facilitating a wide range of applications in atmospheric ozone research48,49.
We collected remote sensing data including OMI Level 3 global daily total ozone grid product50 (OMTO3G; https://disc.gsfc.nasa.gov) and ozone profile products (PROFOZ; v0.9.3, level 2), which is derived using backscattered radiation within the sensitive ultraviolet spectral range for various atmospheric constituents51. The OMI provides daily ozone column concentration (0.25° × 0.25°) data, and the ozone profile product contains 18 vertical layers52, of which the first layer (air pressure of 1000 hPa) was selected to represent surface ozone in this study. We also calculated the average percentage of days with valid OMI data for each grid cell from 2015 to 2020 (Supplementary Figure S1). Most grid cells had a relative high percentage of qualified OMI retrievals, ranging from 64% to 83%. Specifically, the central and eastern regions had an average percentage between 70% and 75%, while the northeastern region had a lower percentage. The percentage in the southern region showed a larger spatial variability.
-
Climate data
Meteorological conditions are key driving factors in the formation process of surface ozone at short time-scales53,54,55. A total of seven climatic variables (solar radiation, air temperature, relative humidity, surface pressure, horizontal wind velocity, vertical wind velocity, and precipitation) were obtained from the ERA5-Land reanalysis data (Supplementary Table S1). Air temperature and solar radiation, which contribute to photochemical reactions, have strong positive correlations with ozone concentration56,57,58, whereas there is a significant negative correlation between ozone concentrations and atmospheric pressure. When the near surface is controlled by low pressure, pollutants from surrounding areas converge towards the center, driven by high-pressure air masses, resulting in a sharp increase in ozone concentrations in the center of the low pressure59. The relative humidity is negatively correlated with ozone concentrations because high relative humidity generally corresponds to precipitation, fog, and other weathers that do not have strong UV radiation, which is not conducive to the occurrence of photochemical reactions and the further development of ozone pollution60. The impact of wind speed on surface ozone concentration is complex. High wind speeds can lead to the dilution of local ozone concentrations, resulting in a negative correlation with the concentrations. However, high wind speeds can also enhance the transport of pollutants downwind, resulting in a positive correlation with downwind ozone concentrations61. As for precipitation, it facilitates the removal of pollutants such as ozone58,62.
The ERA5-Land reanalysis dataset has a spatial resolution of 0.1° × 0.1° (about 9 km) and an hourly time-step, produced by the European Centre for Medium-Range Weather Forecasts (ECMWF; https://www.ecmwf.int/en/forecasts). The ERA5 reanalysis data combines land surface model simulations with ground and satellite observations63,64, and has been widely used across the world65. It has also been validated in China, showing good performance in predicting air temperature66, solar radiation67, and precipitation68.
-
Auxiliary data
Socio-economic data reflects human living and production activities, which is the major source of ozone precursors4 (VOCs and NOx). Existing emission inventories have significant uncertainties and low temporal resolutions, largely restricting their use in predicting hourly surface ozone concentrations69. we used socio-economic data and land use data as an input. We obtained population distribution data and Gross Domestic Product (GDP) data with a 1 km spatial resolution from the Resource and Environmental Science and Data Center, Chinese Academy of Sciences70 (https://www.resdc.cn/DOI). The data has a time interval of five years and is available for two years (2015 and 2019) during the study period. The nationwide land use data was derived from the Moderate Resolution Imaging Spectroradiometer71 (MODIS; https://lpdaac.usgs.gov/products/mcd12c1v006/) product at a resolution of 0.05°. Additionally, we also included latitude and longitude as predictor variables.
-
Data processing
We constructed a 0.1° × 0.1°grid over China and averaged all the concurrent surface ozone measurements of monitoring sites within each grid cell to obtain grid-level surface ozone concentrations. Finally, we obtained a total of 643 grid cells with surface ozone observations across China (Fig. 2). The Thin Plate Spline (TPS) method was used to fill the missing values in OMI total column ozone data (Supplementary Figure S2). TPS has been proven to be effective in interpolating meteorological data72 and handling missing values in OMI remote sensing data, including total column ozone73 and aerosol optical depth74. Correspondingly, all predictor variables (including satellite retrievals, climate, land use, population distribution, and GDP data) were aggregated or resampled to the targeted grid resolution of 0.1° × 0.1° using the nearest neighbor interpolation or the bilinear interpolation approach. To avoid high collinearity among predictor variables, we conducted variance inflation factor (VIF) tests to all the predictor variables and only those with a VIF value less than 8.0 were retained30 (Supplementary Figure S3).
Model development
The long short-term memory network model
The long short-term memory network is a special type of recurrent neural networks (RNN) that differs from traditional ones. The traditional artificial neural network (ANN) is fully connected between layers and has no connection within a specific layer, whereas the hidden layers of RNN are connected75. The output of RNN is not only affected by the current input features but also influenced by the output of the previous moment, and thus RNN generally has a better performance in estimating time-series and has been widely used to proceed sequence data76.
The LSTM can further overcome the limitations of conventional RNNs that they could be trapped by a vanishing gradient or exploding gradient during training75. It excels through integrating input gates, forgetting gates, and output gates into the cell structure. The input gates control whether a cell value can be added to a memory cell, the forgetting gates determine the weight of the value, and the output gates determine which information eventually is output from the cell. The LSTM has a long-term memory capability, which is ideal for predicting long time-series of historical ozone concentrations.
Specifically, based on LSTM, we built a five-layer neural network model for surface ozone concentration prediction (Supplementary Figure S4). It consists of an input layer, two LSTM layers, one Dense layer (also called fully connected layers), and an output layer. The data specification for the model’s input layer is in a 3-dimensional format (n_samples, n_time_steps, n_features), where n_samples represents the batch size for training, n_time_steps is the time window of 24 hours (to determine the optimal time window for training, we conducted several experiments with eight different lookback windows the detailed experiment results are shown in the Supplementary Table S2), representing the first 24 hours’ ozone sequence to predict the ozone at the 24th hour, and n_features is the number of 12 variables in the training set.
To determine the optimal hyperparameters (including epoch, batch size, number of neurons and optimizer), we first conducted a sensitivity analysis to identify the importance of each hyperparameter. Specifically, each hyperparameter was assigned a prior range, and the whole dataset was partitioned into the training data and the validation data using a ratio of 9:1. We adopted a one-at-a-time (OAT) strategy, i.e., changing one parameter at a fixed interval while keeping others unchanged, to avoid consuming too many high-performance computer resources. The results showed that changes in hyperparameters had minor effects on model performance (the R2 and RMSE values were nearly stable, being around 0.7 and 10.00 ppb, respectively). Thus, the mean value of the specific range for each of the hyperparameters was used. Supplementary Figure S5 shows the convergence of the loss function using the final hyperparameters. The number of neurons in each hidden layer is 50, and we used mean absolute error (MAE) as the loss function and the Adaptive moment estimation (Adam) as the optimization algorithm. The model was trained for 50 epochs with a batch size of 3000. The CNEMC ground measurements were used as the target for the model training and validation (Table 1).
Model training, validation and test
All selected predictor variables that passed the VIF tests were used as inputs to train the LSTM model. The importance of different variables in the model was determined using the permutation importance method77, which measures the degree of decline in the model’s performance score after the random rearrangement of different features, and also represents the importance of each variable in estimating the concentration of ozone in the model. The feature importance scores of all selected variables in the pre-trained model are shown in Supplementary Figure S6. Air temperature, surface pressure and relative humidity were the top three factors affecting the spatiotemporal variability of surface ozone concentrations in China. In addition, longitude, day of year (DOY), latitude, downwelling surface radiation, wind speed, land use data, socio-economic data, and OMI’s SFO3 product also have significant impacts on ozone estimation. Finally, total column concentration ozone data and precipitation have relatively weaker influence in the model.
We divided the original data into a training set (more than 600 grid cells during 2015–2020) and a testing set (for the periods of June to December in 2014 and January to February in 2021). To determine the best model and its corresponding hyperparameters, the 10-fold CV approach was utilized to evaluate the performance of the LSTM model on the training set (from 2015 to 2020), with three sampling strategies, namely sample-based CV, site-based CV and by-year CV, corresponding to the model’s performances on capturing overall, spatial, and temporal patterns, respectively. We performed a by-year CV by dividing the dataset into six folds, each representing one year from 2015 to 2020. In each iteration, five folds were used as the training set and the remaining fold was used as the validation set. The training and validation processes were repeated six times. In the other two sampling strategies, 90% of the total surface ozone observations were randomly sampled for training, and the rest 10% was used for validation, the process of which was repeated 10 times. In addition, ozone monitoring station data obtained in 2014 and 2021 was used as a test data set to evaluate the generic capability of the optimal LSTM model. The specific process is shown in Fig. 3.
The R2, RMSE, MAE, linear regression slope, and intercept were calculated to evaluate the performance of the model.
Where the subscript i represents the pairing of n observed ozone concentrations pi and their corresponding predictions oi, and \(\bar{p}\) represents the arithmetic mean of the observed ozone concentrations.
Data Records
The HrSOD dataset78 is available on the Zenodo repository at https://doi.org/10.5281/zenodo.7415326. The gridded ozone concentration data are provided in NetCDF format at 0.1° spatial resolution and hourly temporal resolution during 2015–2020 in ppb. The file size is 40 GB. The hourly data is a NetCDF file and the file is named “YYYYMMDD.nc”, where “YYYY”, “MM” and “DD” refer to the year, month, and data of the file. We have uploaded all the ozone site measurements data41 to the product repository. And this data can be accessed via the link: https://doi.org/10.5281/zenodo.10911197.
This study’s external data include OMI satellite remote sensing data for total column ozone50 and surface ozone concentrations51, available at 10.5067/Aura/OMI/DATA2025 and 10.5067/Aura/OMI/DATA2026 respectively. Climate data65 were sourced from ERA5-land at 10.24381/cds.e2161bac. Socio-economic information, such as population distribution and GDP data70, is accessible through http://www.resdc.cn/DOI. Nationwide land use data71 was derived from MODIS, available at https://lpdaac.usgs.gov/products/mcd12c1v006/.
Technical Validation
Model evaluation
At the hourly time-scale, the LSTM model obtained R2 values of 0.72, 0.65 and 0.71 using three CV sampling methods (sample-based, site-based and by-year), respectively. The corresponding RMSE values were 11.71 ppb, 12.81 ppb, 11.14 ppb, and MAE values were 8.80 ppb, 9.64 ppb and 8.44 ppb (Fig. 4a–c). At the daily time-step, the model’s performance improved with R2 values of 0.71, 0.63, and 0.71 (sample-based, site-based, and by-year), RMSE values of 8.53 ppb, 9.61 ppb, and 7.97 ppb, and MAE values of 6.42 ppb, 7.24 ppb, and 6.09 ppb (Fig. 4d–f). The predictive ability of the model further improved at the monthly time-step, with higher R2 values of 0.82, 0.72, and 0.84 (sample-based, site-based, and by-year), smaller RMSE values of 5.14 ppb, 6.54 ppb, and 4.39 ppb and MAE values of 3.69 ppb, 4.69 ppb, and 3.35 ppb (sample-based, site-based, and by-year) (Fig. 4g–i).
Among the three CV sampling strategies, the site-based CV (Fig. 4b,e,h) had slightly lower R2 values than the sample-based CV (Fig. 4a,d,g) and by-year CV (Fig. 4c,f,i) R2 values, while the RMSE and MAE values were slightly higher than the sample-based CV RMSE and MAE values and by-year CV RMSE and MAE values. It is worth noting that the model tended to underestimate surface ozone when it was at high concentrations, but this bias was largely reduced at the monthly time-step (Fig. 4g–i). In addition, we compared the performance of LSTM with two other commonly used machine learning methods (RF and XGboost) using the same input data at an hourly time-step. The results show that the LSTM model performed better, particularly in terms of R2, RMSE, and slope values (Supplementary Figure S7).
The spatial prediction accuracy of the LSTM model was evaluated based on the values of CV R2 (Fig. 5a), MAE (Fig. 5b), and RMSE (Fig. 5c), which were 0.66, 8.45 ppb, and 11.03 ppb, respectively. The R2 values at around 75% of the monitoring sites ranged from 0.61 to 0.87, and 75% of the monitoring sites had MAE values less than 9.25 ppb and RMSE values less than 12.00 ppb. The model showed a better hourly ozone prediction ability in the North China Plain and the Southwest region of China, with R2 values generally higher than 0.70 (Fig. 5a). Furthermore, the MAE and RMSE values in the southwest region are lower than those in other regions (Fig. 5b,c). However, the model’s uncertainty was higher in the central and eastern regions of China, with MAE values ranging from 8.00 to 11.00 ppb and RMSE values ranging from 11.00 to 14.00 ppb.
The independent test set mainly comprised hourly ozone concentration records obtained in 2014 (June-December; 946 sites) and 2021 (January-February; 1720 sites). Supplementary Figure S8 shows the HrSOD performance (R2 = 0.64, RMSE = 15.44 ppb, MAE = 10.66 ppb) across China at the hourly time-scale. Despite the differences in data distribution and sample size between the test set and the validation set, the model continued to perform well on the test, suggesting that the LSTM model could accurately capture the spatiotemporal patterns of surface ozone concentrations.
In light of lacking direct observations in some regions, we figured out a workaround to validate the reliability of LSTM model. Specifically, the OMI remotely sensed surface ozone concentrations were taken as a benchmark for the whole country. Despite criticism for its low accuracy18, the OMI surface ozone concentration product has a consistent performance both spatially and temporally in areas with and without monitoring stations. Therefore, the OMI surface ozone concentration product is an appropriate choice for evaluating the consistency of the HrSOD product. Figure 6 shows that the HrSOD product has a highly congruous performance against the OMI product in regions with (R² = 0.25, RMSE = 8.18 ppb; nationwide) and without site measurements (R² = 0.23, RMSE = 7.74 ppb; nationwide). Hence, we can conclude that the HrSOD product demonstrates consistent performance across China. We also compared the spatial ozone patterns from HrSOD and OMI data in 2015 (Supplementary Figure S9). The result shows that except in northeastern China, HrSOD and OMI generally show a consistent pattern with higher ozone concentration in the south and lower ozone in the north. Shen et al.18 also observed this pattern in their comparison of surface ozone observations with OMI enhancements, and found that OMI data exhibits relatively weak retrieval sensitivity in the north due to greater upper tropospheric ozone variability there than in the south.
Spatiotemporal variations of surface ozone across China
The diurnal, monthly mean and monthly surface ozone concentrations predicted by the LSTM model were consistent with those observed across China from 2015 to 2020 (Fig. 7). The diurnal variations of mean hourly surface ozone concentrations across China exhibited a unimodal curve. Specifically, the national average for hourly ozone concentrations gradually increased from around 9:00–10:00 (UTC + 8) and peaked at approximately 15:00 (UTC + 8) with a value of about 48.42 ppb based on either site measurements or the HrSOD product (Fig. 7a). Then the mean hourly ozone concentrations gradually decreased to about 20.00–25.00 ppb. Similarly, the mean monthly ozone concentration in China also displayed a unimodal pattern from 2015 to 2020 (Fig. 7b). The ozone concentration gradually increased and peaked in May at 41.73 ppb. Subsequently, the concentration gradually decreased until December, when the surface ozone concentration reached its minimum at about 17.28 ppb. The surface ozone concentrations (Fig. 7c) across China showed regular seasonal changes from 2015 to 2020, and the concentrations gradually increased over this period. It is noteworthy that the ozone concentrations in May 2017 (46.30 ppb) and June 2018 (47.28 ppb) were higher than those in other months.
Figure 8 shows that the spatial distribution of surface ozone concentrations observed cross China from 2015 to 2020 is generally consistent with HrSOD at different time scales. The highest ozone concentration was observed at 15:00 (Fig. 8c,g), while the southwestern region had lower ozone concentrations compared to the North China Plain region at all four times. The multi-year mean seasonal ozone concentrations were predicted to be 37.64 ± 3.35, 39.16 ± 2.37, 28.40 ± 3.17, and 25.07 ± 2.60 ppb in spring (March–May), summer (June–August), autumn (September–November), and winter (December, January, and February), respectively. In springs (Fig. 8i,m), northern and eastern China had higher ozone concentrations. In summer (Fig. 8j,n), the areas with high ozone concentrations (>45.00 ppb) were north China, northwestern China, and southern Inner Mongolia. The hotspot areas with high ozone concentrations in autumns (Fig. 8k,o) decreased and spread to the southeast coast. During winters (Fig. 8l,p), the areas with high ozone concentrations (>30.00 ppb) almost disappeared in southeastern China. The strong spatial and temporal variations in surface O3 concentrations could be attributed to multiple drivers. In densely populated regions, particularly in northern China, industrial air pollutions contribute more to O3 production than in other regions79. In contrast, in northwest China, despite relatively low population, the special topography and strong radiation affect atmospheric diffusion conditions and strengthen photochemical reactions, resulting in generally higher background surface O3 concentrations. In summer time, high temperature and dry air accelerates photochemical reactions in northern China80. Correspondingly, in southern China, due to the influences of monsoon climate in summer, frequent cloud cover and relatively high-water vapor benefit removal of surface ozone. However, surface O3 concentrations become higher in southern China in autumns because of intensified solar radiation during the season. In springs, stratospheric ozone intrusions and photochemical reactions of winter-accumulated precursor compounds contribute to high surface O3 concentrations in some regions, notably in southwest and northeast China81.
Upon comparing predicted and observed mean hourly surface ozone concentrations, The discrepancies mostly fell within a range of −5 to 5 ppb at the hourly scale. The HrSOD estimates tended to underestimate hourly surface ozone concentrations in the majority of China, except some overestimation in the southeast part. Such overestimation was particularly manifest during summers and autumns, with the bias reaching up to approximately 5–10 ppb (Fig. 9f,g).
Surface ozone changes in key regions
Among the four megacity clusters, mean annual surface ozone concentrations in BTH (mean = 32.35 ppb), YRD (mean = 32.78 ppb), and PRD (mean = 27.59 ppb) regions were higher than in the SCB (mean = 25.62 ppb) region during 2015–2020. In the BTH region (Fig. 10a), surface ozone concentrations showed a continuous increase from 28.45 ppb in 2015 to 34.92 ppb in 2018, before decreasing to 33.80 ppb in 2020. In the PRD and YRD regions (Fig. 10c,d), the annual mean surface ozone concentrations showed an obvious increase from 25.32 and 29.54 ppb in 2015 to 34.58 ppb in 2017 and 29.45 ppb in 2019, respectively, and then declined to 33.03 ppb and 28.00 ppb in 2020. Similar to BTH, both regions experienced a growth in surface ozone concentrations before 2017. In contrast, annual mean surface ozone concentrations in the SCB region were relatively low, which showed an increase from 23.35 ppb in 2015 to 27.95 ppb in 2018, a decrease in 2019, followed by a slight rebound in 2020.
The seasonal patterns of surface ozone concentrations varied across the four key regions (Supplementary Figure S10a). From April to July, the monthly mean ozone concentrations were higher than 38.00 ppb across China, while they were less than 24.00 ppb in January, November, and December. In BTH, the ozone concentrations followed a unimodal distribution and gradually increased over time, peaking in June (57.30 ppb) before declining to their lowest value in December (11.81 ppb). Conversely, the other three regions (YRD, SCB, and PRD) showed a bimodal pattern, with the first peak occurring in May (44.22 ppb in YRD, 29.64 ppb in PRD, and 37.80 ppb in SCB), and the second peak occurring in September (40.09 ppb in YRD), October (36.02 ppb in PRD), and August (38.95 ppb in SCB), respectively. The lowest surface ozone concentrations were found to be 16.31 ppb (YRD in December), 22.25 ppb (PRD in January), and 11.45 ppb (SCB in December).
We conducted additional partial correlation analysis to investigate the relationships between regional surface ozone concentrations and meteorological factors at hour scales (Supplementary Figure S10b). The results indicate that temperature and relative humidity are the primary controlling factors of regional surface ozone concentration at the hourly scale. Besides BTH, radiation is relatively important. Relative humidity dominates in the YRD region, horizontal wind speed, temperature, and relative humidity co-regulated surface ozone concentration in the PRD region, and solar radiation dominated in the BTH region.
Comparison with previous studies
We conducted a comparison between HrSOD and two other datasets: the long-term hourly surface ozone mixing ratios dataset (1.25° × 1.875°) estimated by the UK Earth System Model 1-0-LL (UKESM1-0-LL) under the Coupled Model Intercomparison Project Phase 6 (CMIP6), and the ERA5 reanalysis hourly surface ozone dataset simulated by the atmospheric model at a resolution of 0.25° × 0.25°. The validation results are presented in Supplementary Figure S11. It is obvious that there exist large uncertainties in ozone estimates from ERA5 and CMIP6. The CMIP6 datasets (R2 = 0.01, RMSE = 32.06 ppb and MAE = 25.77 ppb) exhibited an overestimation of surface ozone concentrations in western China and underestimated ozone concentrations in the North China Plain (Supplementary Figure S12a), and mainly due to uncertainty in emission inventories, deposition processes or vertical mixing82. Similarly, the ERA5 simulated surface ozone dataset (R2 = 0.01, RMSE = 54.19 ppb and MAE = 49.97 ppb) showed significant deviations from the observed values and displayed a decreasing trend from north to south (Supplementary Figure S12b). This discrepancy primarily arises from the simplified representation of ozone chemistry mechanisms in the ERA5 simulations83. In contrast, our HrSOD product (Supplementary Figure S12c) demonstrated a significant improvement compared to the aforementioned products, exhibiting high consistency with sites measurements (e.g., R2 = 0.71, RMSE = 11.14 ppb and MSE = 8.44 ppb). A notable phenomenon is the comparison of surface ozone concentration in Tibet against Qinghai or Xinjiang regions. While Wei et al.20 reported that ozone concentrations over the Tibet were comparable to those in Xinjiang and Qinghai, many studies19,35,84 align with this study. Even seen from the limited number of ground observations in Tibet, the central and western parts of Tibet had much lower surface ozone concentrations compared to Xinjiang and northeast Qinghai. The OMI remote sensing data can further support this phenomenon. After all, both natural and anthropogenic conditions are so different in Tibet from those in Xinjiang and northeast Qinghai. Such a difference is more likely to be caused by the different algorithms and more ground observations are needed in areas of data scarcity to constrain estimation uncertainties in surface ozone concentrations.
In traditional approaches for predicting surface ozone concentrations, meteorological variables and ozone precursors such as NO2 and VOCs emissions inventory data have been commonly included19,20. While these data have played a significant role in predicting ozone concentrations, it is worth noting that in our study, we did not utilize these specific data as predictive variables. Instead, we relied mainly on hourly and daily satellite retrievals and some auxiliary information (population and GDP data), and the results proved that this strategy of parsimonious inputs could obtain satisfying outcomes. This result may be explained by the estimation of VOCs and NOX emissions inventory data, which is mainly obtained by emission inventory models using emission factors for different emission sources, including power plants, industrial plants, as well as residential, transportation, and agricultural sectors85. Therefore, there may be a strong correlation between ozone precursor emission data and population and GDP distribution data, which may cause data redundancy if they are used as predictor variables simultaneously.
Moreover, there exist scale effects for different environmental variables influencing surface ozone concentrations. At hourly time-scale (diurnal), meteorological conditions, such as radiation, air temperature and air humidity, are critical in ozone formation and destruction22,32. In contrast, precursor emissions19, warming86 and climate patterns87 play more important roles in regulating surface ozone concentrations at daily, seasonal and yearly time-scales. Thus, in the case of adopting instantaneous and particularly the daily remotely sensed ozone products, which have implicitly reflected the long-term changes in emission levels of ozone precursors, there is no necessity to include ozone precursor emissions in the predictive model. In addition, ozone precursor emission data are mostly monthly scale data with low spatial resolution, and due to the diversity of emission sources and the lack of reliable measurement methods, the data uncertainty is high, and it contributes relatively low in the previous ozone prediction models19. In future research, obtaining high-resolution and high-precision data on ozone precursors may become a key focus and direction for ozone estimation. As for whether this finding is limited to ozone or also occurs to other air pollutants, e.g., PM2.5, it is a topic worth discussing. At the hourly to daily time-scales, wind speed and local emissions are the major players in determining PM2.5 concentrations88. Compared to gaseous atmospheric pollutants (such as SO2, NO2, O3 and NH3, which are PM2.5 precursors), impacts of meteorological conditions (except wind) on PM2.5 are relatively minor in China89, which is exactly opposite to those found for ozone20,85. Thus, it is reasonable to conclude that influencing factors critical in predicting different air pollutants at an hourly time-scale should be different.
The driving factors for surface ozone concentrations in key regions
Short-term changes in ozone concentration are particularly sensitive to meteorological conditions54,90. For example, at hourly scale, surface ozone concentrations are significantly influenced by different meteorological conditions and synoptic type, especially strong solar radiation, high temperature and high humidity, which play a key role in photochemical reactions involved in ozone formation32,91. Here relative humidity is identified as the most influential variable in YRD and PRD, air temperature is the most influential variable in BTH, while the two factors show comparable effects in SCB. Notably, the positive surface ozone concentration-air temperature relationship may become increasingly important in the context of global warming, which could lead to an increase in surface ozone concentrations, i.e., “ozone climate penalty”91. For example, Wu et al.86 estimated that ozone levels will increase by 2–5 ppb in a warmer environment by 2050. The mechanism lies in that warming can lead to increased emissions of biogenic VOCs from natural sources, accelerate peroxyacetyl nitrate dissociation, enhance natural soil NOx emissions, influence the efficiency of ozone dry deposition, and thus finally affect atmospheric ozone dynamics87.
In 2020 when COVID-19 lockdown/control occurred, surface ozone dynamics differed in the four major megacity regions. Specifically, surface ozone concentrations in the PRD region decreased, primarily due to a reduction in anthropogenic emissions over a short period and the influence of decreased solar radiation intensity92. This result is consistent with previous research, as reduction in NOx emissions from road traffic could promote surface ozone concentration (NO titration effect; NO + O3 = NO2 + O2). Moreover, lower emissions of inhalable particulate matter, along with higher solar radiation, favor ozone formation93,94. However, this mechanism cannot explain surface ozone concentration changes in all regions. In YRD and BTH regions, surface ozone concentration remained stable, while it even experienced an increase in the PRD region in 2020. This phenomenon suggests that the interactive effects of various natural and anthropogenic factors affecting surface ozone concentrations are complex and efforts for reducing surface ozone pollution should account for this.
Uncertainties and limitations
In this research, uncertainties exist in several aspects. Firstly, the monitoring stations were mainly concentrated in the central-eastern region of China, which may limit the model’s ability to fully capture the relationship between surface ozone concentrations and environmental factors in western China. Additionally, the majority of monitoring stations were located in urban areas, which may restrict the model’s accuracy in estimating surface ozone concentrations in natural and agricultural ecosystems. A typical example is the Taklimakan Desert in northwestern China, which is enclosed by high mountains. Due to the accumulation of ozone precursors emitted from surrounding oasis areas95, surface ozone concentration in this region was relatively high (Supplementary Figure S9) but with its temporal variations more related to natural factors, including solar radiation96 and air temperature97. In comparison with site observations from the center of the desert (Tazhong station, 38°58′N, 83°39′E), the HrSOD values (Fig. 8f,g) were consistent with the observed maximum and minimum hourly surface ozone concentrations (69.2 ppb and less than 20 ppb, respectively, during July 2010-Dec 201797), and the mean daily surface ozone concentration (49.0 ppb during June 2010-March 201296 vs. 51.5 ppb by HrSOD). Nevertheless, the HrSOD could not fully capture the temporal variations at the station during 2010–2017. For instance, in 2015, the mean annual surface ozone concentration estimated by HrSOD (54.6 ppb; Supplementary Figure S9) was rather higher than the observations97. Thus, more extensive and continuing observations are required in the future to improve the accuracy of surface ozone predictions, particularly in regions of data scarcity. Secondly, uncertainties can arise from the input data. For example, ERA5 reanalysis data underestimates surface temperatures in the coastal urban agglomerations of southeast China and the Tibetan Plateau66,98, which may lead the model to underestimate ozone concentrations. Enhancing the accuracy of meteorological data, land use maps, and socio-economic data is necessary to further improve ozone estimation accuracy. Furthermore, the mismatch in temporal resolution between OMI remote sensing data and ozone measurements may also affect the final estimation accuracy.
While the LSTM networks effectively capture the temporal variations of surface ozone concentrations, spatial information such as changes in pollutant concentrations due to the emission and transport of surrounding pollutants is not fully considered. The underestimation of the LSTM model in southeast China but overestimation in other parts (Figs. 8, 9) also underlines there may exist some deficit in the trained model in capturing the spatial heterogeneity of surface ozone concentrations across different environmental conditions. Therefore, to enhance the current deep learning model, combining it with other algorithms that could effectively extract spatial dependencies within data may be beneficial. The most widely used framework is integrating CNN with LSTM to leverage the strengths of temporal memory by LSTM and feature representation by CNN for improved prediction accuracy. To validate the plausibility of this methodology, we conducted extra simulations using the Convolutional Long Short-Term Memory (ConvLSTM) algorithm. However, ConvLSTM performed only slightly better than LSTM, at the cost of much more model parameters and computation resource consumption (Supplementary Table S3) and with the spatial biases not resolved (Supplementary Figure S13). This unexpected result suggests that more efforts are warranted in developing novel algorithms to address the fundamental challenge in considering both spatial and temporal information inherently embedded in environmental datasets.
Potential applications of HrSOD
Compared to the currently available surface ozone products in China, HrSOD offers several advantages. It covers a longer time range and has a higher temporal resolution, enabling more robust historical environmental impact and human health risk assessments. HrSOD can be used to derive various ozone exposure indicators (Supplementary Figure S14), such as seasonal 7-hour mean ozone concentrations (M7), seasonal 12-hour mean ozone concentrations99 (M12), sum of all hourly average concentrations >60 μg kg−1 (SUM06)100, cumulative ozone exposure index based on sigmoid-weighted daytime ozone concentrations101 (W126), and accumulated hourly ozone concentration over a threshold of X μg kg−1 during daylight hours102 (AOTX). Therefore, HrSOD can cater to the requirements of ozone impact models and provide flexibility for assessing ozone effects on ecosystems38 and epidemiological studies103.
Code availability
The code is available on GitHub (https://github.com/Wenxiu0902/Ozone_prediction) primarily using Python and R languages. It includes data preprocessing, model training, testing, prediction, and visualization sections. Additionally, sample model input data is also provided.
References
Norval, M. et al. The effects on human health from stratospheric ozone depletion and its interactions with climate change. Photochemical & Photobiological Sciences 6, 232–251, https://doi.org/10.1039/B700018A (2007).
Slaper, H., Velders, G. J., Daniel, J. S., de Gruijl, F. R. & van der Leun, J. C. Estimates of ozone depletion and skin cancer incidence to examine the Vienna Convention achievements. Nature 384, 256–258, https://doi.org/10.1038/384256a0 (1996).
van der Leun, J., Tang, X. & Tevini, M. Environmental effects of ozone depletion and its interactions with climate change: 2002 assessment. Photochemical & Photobiological Sciences 2, vii–vii, https://doi.org/10.1039/b211913g (2003).
Wang, T. et al. Ozone pollution in China: A review of concentrations, meteorological influences, chemical precursors, and effects. Science of The Total Environment 575, 1582–1596, https://doi.org/10.1016/j.scitotenv.2016.10.081 (2017).
Berman, J. D. et al. Health benefits from large-scale ozone reduction in the United States. Environmental Health Perspectives 120, 1404–1410, https://doi.org/10.1289/ehp.1104851 (2012).
Li, H. et al. Short-term effects of various ozone metrics on cardiopulmonary function in chronic obstructive pulmonary disease patients: Results from a panel study in Beijing, China. Environmental Pollution 232, 358–366, https://doi.org/10.1016/j.envpol.2017.09.030 (2018).
Magzamen, S., Moore, B. F., Yost, M. G., Fenske, R. A. & Karr, C. J. Ozone-related respiratory morbidity in a low-pollution region. Journal of Occupational and Environmental Medicine 59, 624–630, https://doi.org/10.1097/jom.0000000000001042 (2017).
Cooper, O. et al. Global distribution and trends of tropospheric ozone: An observation-based review. Elementa: Science of the Anthropocene 2, 000029, https://doi.org/10.12952/journal.elementa.000029 (2014).
Giles, J. Hikes in surface ozone could suffocate crops. Nature 435, 7–7, https://doi.org/10.1038/435007a (2005).
Lu, X. et al. Severe surface ozone pollution in China: a global perspective. Environmental Science & Technology Letters 5, 487–494, https://doi.org/10.1021/acs.estlett.8b00366 (2018).
Tian, H. et al. Climate extremes and ozone pollution: a growing threat to china’s food security. Ecosystem Health and Sustainability 2, e01203, https://doi.org/10.1002/ehs2.1203 (2016).
Huang, J., Pan, X., Guo, X. & Li, G. Health impact of China’s air pollution prevention and control action plan: an analysis of national air quality monitoring and mortality data. The Lancet. Planetary health 2, e313–e323, https://doi.org/10.1016/s2542-5196(18)30141-4 (2018).
Ma, Z. et al. Significant increase of surface ozone at a rural site, north of eastern China. Atmos. Chem. Phys. 16, 3969–3977, https://doi.org/10.5194/acp-16-3969-2016 (2016).
Maji, K. J., Ye, W.-F., Arora, M. & Nagendra, S. M. S. Ozone pollution in Chinese cities: Assessment of seasonal variation, health effects and economic burden. Environmental Pollution 247, 792–801, https://doi.org/10.1016/j.envpol.2019.01.049 (2019).
Sahu, S. K., Liu, S., Liu, S., Ding, D. & Xing, J. Ozone pollution in China: Background and transboundary contributions to ozone concentration & related health effects across the country. Science of The Total Environment 761, 144131, https://doi.org/10.1016/j.scitotenv.2020.144131 (2021).
Chang, K.-L., Petropavlovskikh, I., Cooper, O. R., Schultz, M. G. & Wang, T. Regional trend analysis of surface ozone observations from monitoring networks in eastern North America, Europe and East Asia. Elementa: Science of the Anthropocene 5, https://doi.org/10.1525/elementa.243 (2017).
Liu, X., Bhartia, P. K., Chance, K., Spurr, R. J. D. & Kurosu, T. P. Ozone profile retrievals from the Ozone Monitoring Instrument. Atmos. Chem. Phys. 10, 2521–2537, https://doi.org/10.5194/acp-10-2521-2010 (2010).
Shen, L. et al. An evaluation of the ability of the Ozone Monitoring Instrument (OMI) to observe boundary layer ozone pollution across China: application to 2005–2017 ozone trends. Atmos. Chem. Phys. 19, 6551–6560, https://doi.org/10.5194/acp-19-6551-2019 (2019).
Liu, R. et al. Spatiotemporal distributions of surface ozone levels in China from 2005 to 2017: A machine learning approach. Environment International 142, 105823, https://doi.org/10.1016/j.envint.2020.105823 (2020).
Wei, J. et al. Full-coverage mapping and spatiotemporal variations of ground-level ozone (O3) pollution from 2013 to 2020 across China. Remote Sensing of Environment 270, 112775, https://doi.org/10.1016/j.rse.2021.112775 (2022).
Liu, H. et al. Ground-level ozone pollution and its health impacts in China. Atmospheric Environment 173, 223–230, https://doi.org/10.1016/j.atmosenv.2017.11.014 (2018).
Sun, L. et al. Impacts of meteorology and emissions on summertime surface ozone increases over central eastern China between 2003 and 2015. Atmos. Chem. Phys. 19, 1455–1469, https://doi.org/10.5194/acp-19-1455-2019 (2019).
Travis, K. R. et al. Why do models overestimate surface ozone in the Southeast United States? Atmos. Chem. Phys. 16, 13561–13577, https://doi.org/10.5194/acp-16-13561-2016 (2016).
Adam-Poupart, A., Brand, A., Fournier, M., Jerrett, M. & Smargiassi, A. Spatiotemporal modeling of ozone levels in Quebec (Canada): A comparison of Kriging, Land-Use Regression (LUR), and Combined Bayesian Maximum Entropy–LUR approaches. Environmental Health Perspectives 122, 970–976, https://doi.org/10.1289/ehp.1306566 (2014).
Chen, L. et al. A hybrid approach to estimating long-term and short-term exposure levels of ozone at the national scale in China using land use regression and Bayesian maximum entropy. Science of The Total Environment 752, 141780, https://doi.org/10.1016/j.scitotenv.2020.141780 (2021).
Zhang, X. Y., Zhao, L. M., Cheng, M. M. & Chen, D. M. Estimating ground-level ozone concentrations in eastern China using satellite-based precursors. IEEE Transactions on Geoscience and Remote Sensing 58, 4754–4763, https://doi.org/10.1109/TGRS.2020.2966780 (2020).
Jumin, E. et al. Machine learning versus linear regression modelling approach for accurate ozone concentrations prediction. Engineering Applications of Computational Fluid Mechanics 14, 713–725, https://doi.org/10.1080/19942060.2020.1758792 (2020).
Pak, U., Kim, C., Ryu, U., Sok, K. & Pak, S. A hybrid model based on convolutional neural networks and long short-term memory for ozone concentration prediction. Air Quality, Atmosphere & Health 11, 883–895, https://doi.org/10.1007/s11869-018-0585-1 (2018).
Eslami, E., Choi, Y., Lops, Y. & Sayeed, A. A real-time hourly ozone prediction system using deep convolutional neural network. Neural Computing and Applications 32, 8783–8797, https://doi.org/10.1007/s00521-019-04282-x (2020).
Cheng, M. et al. Spatio-temporal hourly and daily ozone forecasting in China using a Hybrid machine learning model: Autoencoder and generative adversarial networks. Journal of Advances in Modeling Earth Systems 14, e2021MS002806, https://doi.org/10.1029/2021MS002806 (2022).
Li, G. et al. Widespread and persistent ozone pollution in eastern China during the non-winter season of 2015: observations and source attributions. Atmos. Chem. Phys. 17, 2759–2774, https://doi.org/10.5194/acp-17-2759-2017 (2017).
Mousavinezhad, S., Choi, Y., Pouyaei, A., Ghahremanloo, M. & Nelson, D. L. A comprehensive investigation of surface ozone pollution in China, 2015–2019: Separating the contributions from meteorology and precursor emissions. Atmospheric Research 257, 105599, https://doi.org/10.1016/j.atmosres.2021.105599 (2021).
Wang, W.-N. et al. Assessing spatial and temporal patterns of observed ground-level ozone in China. Scientific Reports 7, 3651, https://doi.org/10.1038/s41598-017-03929-w (2017).
Li, M. et al. Large scale control of surface ozone by relative humidity observed during warm seasons in China. Environmental Chemistry Letters 19, 3981–3989, https://doi.org/10.1007/s10311-021-01265-0 (2021).
Xue, T. et al. Estimating spatiotemporal variation in ambient ozone exposure during 2013–2017 using a data-fusion model. Environmental Science & Technology 54, 14877–14888, https://doi.org/10.1021/acs.est.0c03098 (2020).
Kong, L. et al. A 6-year-long (2013–2018) high-resolution air quality reanalysis dataset in China based on the assimilation of surface observations from CNEMC. Earth Syst. Sci. Data 13, 529–570, https://doi.org/10.5194/essd-13-529-2021 (2021).
Feng, Z. et al. A stomatal ozone flux–response relationship to assess ozone-induced yield loss of winter wheat in subtropical China. Environmental Pollution 164, 16–23, https://doi.org/10.1016/j.envpol.2012.01.014 (2012).
Ren, W. et al. Effects of tropospheric ozone pollution on net primary productivity and carbon storage in terrestrial ecosystems of China. 112, https://doi.org/10.1029/2007JD008521 (2007).
MILLS, G. et al. Evidence of widespread effects of ozone on crops and (semi-)natural vegetation in Europe (1990–2006) in relation to AOT40- and flux-based risk maps. Global Change Biology 17, 592–613, https://doi.org/10.1111/j.1365-2486.2010.02217.x (2011).
Niu, Y. et al. Long-term ozone exposure and small airway dysfunction: The China Pulmonary Health (CPH) study. American journal of respiratory and critical care medicine 205, 450–458, https://doi.org/10.1164/rccm.202107-1599OC (2022).
CNEMC (China National Environmental Monitoring Centre). Hourly surface ozone observations across China. Zenodo. https://doi.org/10.5281/zenodo.10911197 (2024).
Ambient Air Quality Standard. GB3095-2012. (Ministry of Ecology and Environment, 2012).
Technical Specification for Ambient Air Quality Assessment. HJ663-2013. (Ministry of Ecology and Environment, 2013).
Bak, J. et al. Temporal variability of tropospheric ozone and ozone profiles in the Korean Peninsula during the East Asian summer monsoon: insights from multiple measurements and reanalysis datasets. Atmos. Chem. Phys. 22, 14177–14187, https://doi.org/10.5194/acp-22-14177-2022 (2022).
Huang, G. et al. Validation of 10-year SAO OMI Ozone Profile (PROFOZ) product using ozonesonde observations. Atmos. Meas. Tech. 10, 2455–2475, https://doi.org/10.5194/amt-10-2455-2017 (2017).
Antón, M. et al. Validation of OMI-TOMS and OMI-DOAS total ozone column using five Brewer spectroradiometers at the Iberian peninsula. Journal of Geophysical Research: Atmospheres 114, https://doi.org/10.1029/2009JD012003 (2009).
Hu, Y. et al. Study on calculation and validation of tropospheric ozone by ozone monitoring instrument – microwave limb sounder over China. International Journal of Remote Sensing 41, 9101–9120, https://doi.org/10.1080/01431161.2020.1800124 (2020).
Chen, Y. et al. Research on the ozone formation sensitivity indicator of four urban agglomerations of China using Ozone Monitoring Instrument (OMI) satellite data and ground-based measurements. Science of The Total Environment 869, 161679, https://doi.org/10.1016/j.scitotenv.2023.161679 (2023).
Zhao, F. et al. Ozone profile retrievals from TROPOMI: Implication for the variation of tropospheric ozone during the outbreak of COVID-19 in China. Science of The Total Environment 764, 142886, https://doi.org/10.1016/j.scitotenv.2020.142886 (2021).
Pawan, K. OMI/Aura Ozone (O3) total column daily L2 global ridded 0.25 degree × 0.25 degree. GESDISC https://doi.org/10.5067/Aura/OMI/DATA2025 (2012).
Johan, D. & Pepijn, V. OMI/Aura Ozone (O3) Profile 1-Orbit L2 Swath 13x48km V003. GESDISC https://doi.org/10.5067/Aura/OMI/DATA2026 (2009).
McPeters, R. et al. Validation of the Aura Ozone Monitoring Instrument total column ozone product. J. Geophys. Res 113, D15S14, https://doi.org/10.1029/2007JD008802 (2008).
Cheng, N. et al. Ground ozone variations at an urban and a rural station in Beijing from 2006 to 2017: Trend, meteorological influences and formation regimes. Journal of Cleaner Production 235, 11–20, https://doi.org/10.1016/j.jclepro.2019.06.204 (2019).
Shan, W., Yin, Y., Zhang, J., Ji, X. & Deng, X. Surface ozone and meteorological condition in a single year at an urban site in central–eastern China. Environmental Monitoring and Assessment 151, 127–141, https://doi.org/10.1007/s10661-008-0255-0 (2009).
Tarasova, O. A. & Karpetchko, A. Y. Accounting for local meteorological effects in the ozone time-series of Lovozero (Kola Peninsula). Atmos. Chem. Phys. 3, 941–949, https://doi.org/10.5194/acp-3-941-2003 (2003).
Lee, Y. C. et al. Increase of ozone concentrations, its temperature sensitivity and the precursor factor in South China. Tellus B: Chemical and Physical Meteorology 66, 23455, https://doi.org/10.3402/tellusb.v66.23455 (2014).
Xu, W. Y. et al. Characteristics of pollutants and their correlation to meteorological conditions at a suburban site in the North China Plain. Atmos. Chem. Phys. 11, 4353–4369, https://doi.org/10.5194/acp-11-4353-2011 (2011).
Chen, Z. et al. Understanding the causal influence of major meteorological factors on ground ozone concentrations across China. Journal of Cleaner Production 242, 118498, https://doi.org/10.1016/j.jclepro.2019.118498 (2020).
Kovač-Andrić, E., Brana, J. & Gvozdić, V. Impact of meteorological factors on ozone concentrations modelled by time series analysis and multivariate statistical methods. Ecological Informatics 4, 117–122, https://doi.org/10.1016/j.ecoinf.2009.01.002 (2009).
Belan, B. D. & Savkin, D. E. The role of air humidity in variations in near-surface ozone concentration. Atmospheric and Oceanic Optics 32, 586–589, https://doi.org/10.1134/S1024856019050038 (2019).
Topçu, S. & Incecik, S. Surface ozone measurements and meteorological influences in the urban atmosphere of Istanbul. International Journal of Environment and Pollution 17, 390–404, https://doi.org/10.1504/IJEP.2002.000680 (2002).
Li, K. et al. Meteorological and chemical impacts on ozone formation: A case study in Hangzhou, China. Atmospheric Research 196, 40–52, https://doi.org/10.1016/j.atmosres.2017.06.003 (2017).
Albergel, C. et al. ERA-5 and ERA-Interim driven ISBA land surface model simulations: which one performs better? Hydrol. Earth Syst. Sci. 22, 3515–3532, https://doi.org/10.5194/hess-22-3515-2018 (2018).
Hersbach, H. et al. The ERA5 global reanalysis. Q J R Meteorol Soc. 146, 1999–2049, https://doi.org/10.1002/qj.3803 (2020).
Muñoz-Sabater, J. et al. ERA5-Land hourly data from 1950 to present. CDS https://doi.org/10.24381/cds.e2161bac (2019).
Zou, J. et al. Performance of air temperature from ERA5-Land reanalysis in coastal urban agglomeration of Southeast China. Science of The Total Environment 828, 154459, https://doi.org/10.1016/j.scitotenv.2022.154459 (2022).
Jiang, H., Yang, Y., Bai, Y. & Wang, H. Evaluation of the total, direct, and diffuse solar radiations from the ERA5 reanalysis data in China. IEEE Geoscience and Remote Sensing Letters 17, 47–51, https://doi.org/10.1109/LGRS.2019.2916410 (2020).
Jiang, Q. et al. Evaluation of the ERA5 reanalysis precipitation dataset over Chinese Mainland. Journal of Hydrology 595, 125660, https://doi.org/10.1016/j.jhydrol.2020.125660 (2021).
Li, S. et al. Emission trends of air pollutants and CO2 in China from 2005 to 2021. Earth Syst. Sci. Data 15, 2279–2294, https://doi.org/10.5194/essd-15-2279-2023 (2023).
Xu, X. Spatial distribution of GDP in China with km grid dataset. RESDC http://www.resdc.cn/DOI (2017).
Friedl, M. & Sulla-Menashe, D. MCD12C1 MODIS/Terra + Aqua Land Cover Type Yearly L3 Global 0.05Deg CMG V006. NASA EOSDIS Land Processes DAAC https://doi.org/10.5067/MODIS/MCD12C1.006 (2015).
Zhang, X. et al. Decadal trends in wet sulfur deposition in China estimated from OMI SO2 columns. 123, 10,796–710,811, https://doi.org/10.1029/2018JD028770 (2018).
Tcherkezova, E., Kaleyna, P. & Mukhtarov, P. J. B. G. J. Modelling spatial distribution of global total column ozone in QGIS and GRASS GIS environment. 39, 26-37, http://www.niggg.bas.bg/en/about-us/periodicals/bulgarian-geophysical-journal/2013-vol-39/ (2013).
Halos, S. H., Al-Jiboori, M. H., Al-Taai, O. T. & Halos, S. J. I. J. O. N. S. Aerosol optical properties estimation over Iraq and surrounding regions using best GIS spatial interpolation method. Indian Journal Of Natural Sciences 7, 11648–11654, https://www.researchgate.net/publication/344387764_Aerosol_Optical_Properties_Estimation_over_Iraq_and_Surrounding_Regions_using_Best_GIS_Spatial_Interpolation_Method (2016).
Hochreiter, S. & Schmidhuber, J. Long Short-Term Memory. Neural Computation 9, 1735–1780, https://doi.org/10.1162/neco.1997.9.8.1735 (1997).
Goodfellow, I., Bengio, Y. & Courville, A. in Deep Learning (MIT Press, 2016).
François, D., Wertz, V., & Verleysen, M., The permutation test for feature selection by mutual information, in Proceedings of European Symposium on Artificial Neural Networks, Bruges, Belgium, 26–28 April, 239–244, 2006.89 (2006).
Zhang, W. X., Liu, D. & Shi, H. Hourly Surface Ozone data (HrSOD) across China during 2005–2020. Zenodo https://doi.org/10.5281/zenodo.7415326 (2022).
Yin, C. Q. et al. Geographical distribution of ozone seasonality over China. Science of The Total Environment 689, 625–633, https://doi.org/10.1016/j.scitotenv.2019.06.460 (2019).
Wang, W. et al. Long-term trend of ozone pollution in China during 2014–2020: distinct seasonal and spatial characteristics and ozone sensitivity. Atmos. Chem. Phys. 22, 8935–8949, https://doi.org/10.5194/acp-22-8935-2022 (2022).
Stohl, A. et al. Stratosphere-troposphere exchange: A review, and what we have learned from STACCATO. Journal of Geophysical Research: Atmospheres 108, https://doi.org/10.1029/2002JD002490 (2003).
Turnock, S. T. et al. Historical and future changes in air pollutants from CMIP6 models. Atmos. Chem. Phys. 20, 14547–14579, https://doi.org/10.5194/acp-20-14547-2020 (2020).
Knowland, K. E., Ott, L. E., Duncan, B. N. & Wargan, K. Stratospheric intrusion-influenced ozone air quality exceedances investigated in the NASA MERRA-2 reanalysis. Geophysical Research Letters 44, 10,691–610,701, https://doi.org/10.1002/2017GL074532 (2017).
Li, K. et al. Increases in surface ozone pollution in China from 2013 to 2019: anthropogenic and meteorological influences. Atmos. Chem. Phys. 20, 11423–11433, https://doi.org/10.5194/acp-20-11423-2020 (2020).
Wang, Q. G., Han, Z., Wang, T. & Zhang, R. Impacts of biogenic emissions of VOC and NOx on tropospheric ozone during summertime in eastern China. Science of The Total Environment 395, 41–49, https://doi.org/10.1016/j.scitotenv.2008.01.059 (2008).
Wu, S. et al. Effects of 2000–2050 global change on ozone air quality in the United States. Journal of Geophysical Research: Atmospheres 113, https://doi.org/10.1029/2007JD008917 (2008).
Fu, H. et al. Investigating PM2.5 responses to other air pollutants and meteorological factors across multiple temporal scales. Scientific Reports 10, 15639, https://doi.org/10.1038/s41598-020-72722-z (2020).
Feng, R., Gao, H., Luo, K. & Fan, J.-R. Analysis and accurate prediction of ambient PM2.5 in China using Multi-layer Perceptron. Atmospheric Environment 232, 117534, https://doi.org/10.1016/j.atmosenv.2020.117534 (2020).
Pu, X. et al. Enhanced surface ozone during the heat wave of 2013 in Yangtze River Delta region, China. Science of The Total Environment 603-604, 807–816, https://doi.org/10.1016/j.scitotenv.2017.03.056 (2017).
Shu, L. et al. Integrated studies of a regional ozone pollution synthetically affected by subtropical high and typhoon system in the Yangtze River Delta region, China. Atmospheric Chemistry and Physics 16, 15801–15819, https://doi.org/10.5194/acp-16-15801-2016 (2016).
Rasmussen, D. J., Hu, J., Mahmud, A. & Kleeman, M. J. The ozone–climate penalty: past, present, and future. Environmental Science & Technology 47, 14258–14266, https://doi.org/10.1021/es403446m (2013).
Liu, S. et al. Distinct regimes of O3 response to COVID-19 lockdown in China. Atmosphere 12, 184, https://doi.org/10.3390/atmos12020184 (2021).
Sicard, P. et al. Amplified ozone pollution in cities during the COVID-19 lockdown. The Science of the total environment 735, 139542, https://doi.org/10.1016/j.scitotenv.2020.139542 (2020).
Wang, N. et al. Air Quality During COVID-19 Lockdown in the Yangtze River Delta and the Pearl River Delta: two different responsive mechanisms to emission reductions in China. Environmental Science & Technology 55, 5721–5730, https://doi.org/10.1021/acs.est.0c08383 (2021).
Liu, X. et al. The variation characteristics and effect factors of surface ozone concentration in the Taklimakan Desert hinterland. Sci. Cold Arid. Reg. 6, 81–88 (2014).
Wang, H., Ma, J., Shen, Y. & Wang, Y. Assessment of ozone variations and meteorological influences at a rural site in northern Xinjiang. Bulletin of Environmental Contamination and Toxicology 94, 240–246, https://doi.org/10.1007/s00128-014-1451-y (2015).
Liu, X. et al. Observational study of ground-level ozone in the desert atmosphere. Bulletin of Environmental Contamination and Toxicology 108, 219–224, https://doi.org/10.1007/s00128-021-03444-9 (2022).
Li, Y. Z. et al. Evaluation of long-term and high-resolution gridded precipitation and temperature products in the Qilian mountains, Qinghai-Tibet Plateau. Frontiers in environmental science 10, https://doi.org/10.3389/fenvs.2022.906821 (2022).
Legge, A. H. et al. Ambient ozone and adverse crop response: An evaluation of North American and European data as they relate to exposure indices and critical levels. J Appl Bot Food Qual 69, 192–205, https://www.cabidigitallibrary.org/doi/full/10.5555/19960704827 (1995).
Lefohn, A. S. & Foley, J. K. NCLAN results and their application to the standard-setting process: Protecting vegetation from surface ozone exposures. Journal of the Air & Waste Management Association 42, 1046–1052, https://doi.org/10.1080/10473289.1992.10467049 (1992).
Fuhrer, J., Skärby, L. & Ashmore, M. R. Critical levels for ozone effects on vegetation in Europe. Environmental Pollution 97, 91-106, https://doi.org/10.1016/S0269-7491(97)00067-5 (1997).
LRTAR Convention, Draft Chapter III: Mapping critical levels for vegetation, of the manual on methodologies and criteria for modelling and mapping critical loads and levels and air pollution effects, risks and trends, http://icpmapping.org/Mapping_Manual (2015).
Huangfu, P. & Atkinson, R. Long-term exposure to NO2 and O3 and all-cause and respiratory mortality: A systematic review and meta-analysis. Environment International 144, 105998, https://doi.org/10.1016/j.envint.2020.105998 (2020).
Acknowledgements
This work was supported by the National Key Research and Development Program of China (Grant No. 2023YFF1303700 and Grant No. 2018YFA0606001). H.T. acknowledges funding support from the National Science Foundation of the United States (Grant No. 1903722) and Andrew Carnegie fellowship Program (Grant No. G-F-19-56910).
Author information
Authors and Affiliations
Contributions
W.Z. and D.L. contributed equally, with W.Z. performing the data curation, modelling, validation, and writing the original draft of the paper. D.L. was responsible for conceptualization, providing data, and reviewed the manuscript. H.S. was the lead and corresponding author of this work, supported and supervised the study and reviewed the paper. H.T. and S.W. also conceptualized the project and supervised the simulations and analyses. R.Y., W.T. and H.M. provided coding and developed the model. H.T., N.P., J.Y., F.L., B.D. and S.W. contributed to the writing of the manuscript. Correspondence to Hao Shi (haoshi@recces.ac.cn).
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, W., Liu, D., Tian, H. et al. Parsimonious estimation of hourly surface ozone concentration across China during 2015–2020. Sci Data 11, 492 (2024). https://doi.org/10.1038/s41597-024-03302-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41597-024-03302-3