Next Article in Journal
Impact of COVID-19 Pandemic on Daily Life, Physical Exercise, and General Health among Older People with Type 2 Diabetes: A Qualitative Interview Study
Previous Article in Journal
Educational, Emotional, and Social Impact of the Emergency State of COVID-19 on Romanian University Students
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regional Prediction of Ozone and Fine Particulate Matter Using Diffusion Convolutional Recurrent Neural Network

1
Center for Intelligent Transportation Systems and Unmanned Aerial Systems Applications Research, State Key Laboratory of Ocean Engineering, School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
2
International Center for Adaptation Planning and Design, College of Design, Construction and Planning, University of Florida, P.O. Box 115706, Gainesville, FL 32611, USA
3
Department of Computer Science, Viterbi School of Engineering, University of Southern California, Los Angeles, CA 90089, USA
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(7), 3988; https://doi.org/10.3390/ijerph19073988
Submission received: 13 January 2022 / Revised: 13 March 2022 / Accepted: 25 March 2022 / Published: 27 March 2022
(This article belongs to the Section Environmental Science and Engineering)

Abstract

:
Accurate air quality forecasts can provide data-driven supports for governmental departments to control air pollution and further protect the health of residents. However, existing air quality forecasting models mainly focus on site-specific time series forecasts at a local level, and rarely consider the spatiotemporal relationships among regional monitoring stations. As a novelty, we construct a diffusion convolutional recurrent neural network (DCRNN) model that fully considers the influence of geographic distance and dominant wind direction on the regional variations in air quality through different combinations of directed and undirected graphs. The hourly fine particulate matter (PM2.5) and ozone data from 123 air quality monitoring stations in the Yangtze River Delta, China are used to evaluate the performance of the DCRNN model in the regional prediction of PM2.5 and ozone concentrations. Results show that the proposed DCRNN model outperforms the baseline models in prediction accuracy. Compared with the undirected graph model, the directed graph model considering the effects of wind direction performs better in 24 h predictions of pollutant concentrations. In addition, more accurate forecasts of both PM2.5 and ozone are found at a regional level where monitoring stations are distributed densely rather than sparsely. Therefore, the proposed model can assist environmental researchers to further improve the technologies of air quality forecasts and could also serve as tools for environmental policymakers to implement pollution control measures.

1. Introduction

In recent decades, developing countries such as China have experienced rapid economic growth and urbanization, and the substantial problems of urban air pollution have also emerged [1,2]. For example, frequent occurrences of haze weather have attracted worldwide concerns due to the deterioration of urban particulate pollution closely related to intensive emissions of fine particles (PM2.5) and coarse particles (PM10). According to epidemiological studies, long-term exposure to higher concentration levels of particulate matter (PM) can cause serious health risks, such as cardiovascular diseases, respiratory diseases, and even deaths [3,4]. Additionally, PM-related visibility reduction also brings negative effects on human production and daily life [5]. Tropospheric ozone is another air pollutant and one of the most important greenhouse gases. Ozone can participate in various atmospheric photochemical processes and contribute to the indirect formation of secondary particulate matter, which also exhibits very detrimental effects on human health [6,7]. Therefore, the accurate prediction of particulate matter and ozone in a high spatial and temporal resolution is essential in assisting governmental departments to design and implement the related emission and pollution control policies, and the refined regional forecasts of PM and ozone concentrations could strongly enhance the protection of public health.
Related research on air quality forecasts has mainly been conducted with deterministic models and empirical models. The deterministic models are developed based on atmospheric physics and mechanisms, such as the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) [8,9] and Community Multi-scale Air Quality (CMAQ) model [10]. The models can help researchers understand the physical and chemical formation mechanisms of urban air pollution but require lots of parameters including meteorological conditions and emission inventories as inputs, thus bringing high uncertainty in air quality forecasts. Furthermore, emission inventories may change over time and need to be updated from time to time, which are usually accompanied by larger amounts of costs and difficulties. By contrast, the empirical models are built based on statistics and machine learning approaches. Although the models could not deeply explore the meteorological and chemical coupling patterns behind air pollution, they often only require a few extra parameters as model inputs in addition to historical data of air pollutants. Hence, statistical methods and machine learning models are widely used to achieve the time series predictions of urban air pollution, such as Auto-regressive Integrated Moving Average (ARIMA) [11], Support Vector Machine (SVM) [12], Classification and Regression Tree (CART) [13], etc. These models could fully cover different types of inputs including air pollutants, meteorology, land use, etc., but the spatial characteristics of air pollution data cannot be sufficiently incorporated in the process of model construction.
Deep learning models have presented excellent performances in the fields of air quality forecasts and environmental assessments, such as the Recurrent Neural Network (RNN) [14], Long Short-Term Memory (LSTM) [15,16,17,18,19], and Gated Recurrent Unit (GRU) [20,21]. However, these RNN-based models cannot fully consider the physical characteristics of topology networks among air quality monitoring stations to characterize the spatial correlations [22]. Motivated by the potentials of the Convolutional Neural Network (CNN) in capturing the spatial relationships, CNN-LSTM [22] and Convolutional LSTM [23] are employed to conduct air quality predictions by extracting both spatial and temporal features in grid-structured data (e.g., images). However, the CNN-based models only extract the spatial features of the research target areas from the input grid-structured data and are incapable of modeling complex topological relationships among large-scale air quality monitoring network. The graph neural network (GNN) is an emerging deep learning model that could map the complex topological relationships of certain spatial areas into a low-dimensional matrix. With this potential, the GNN have been widely applied in air quality forecast tasks [24,25], and outperforms the common deep learning models mentioned above. Generally, these GNN models [26] use undirected graphs to capture the topological relationships among air quality monitoring network, but the undirected edges between nodes cannot consider the effects of wind direction in the graph. Considering that air pollutants can be transported following the direction of the wind, an integration of the wind factors into the GNN model may improve the prediction performance. Thus, an exploration of the GNN model driven by a directed graph is necessary to model the influence of wind factors and further improve regional air quality prediction.
To fill the above research gaps, we develop a novel diffusion convolutional recurrent neural network (DCRNN) model for the regional prediction of PM2.5 and ozone concentrations. Specifically, different graph construction methods including undirected and directed graphs are separately integrated into the proposed DCRNN model to fully consider the spatial relationships among air quality monitoring stations. The main difference between undirected and directed graphs lies in that the directed graph considers the network-level dominant wind direction as an extra and important spatial dependency among monitoring stations while the undirected graph does not. Then, we further evaluate the performances of different DCRNN models in forecasting the PM2.5 and ozone concentrations and compare them with baseline models for various prediction lengths of time. Finally, we discuss the influence of spatial characteristics within the proposed DCRNN model on the accuracy of air quality forecasts.

2. Data and Methods

2.1. Data Description

The Yangtze River Delta region is selected as the study area covering three provincial units (i.e., Zhejiang, Jiangsu, and Shanghai), and its latitude and longitude range from 27° N to 35° N and 116° E to 123° E, respectively. Figure 1 shows the geographical location of the study area and the distributions of 123 air quality monitoring stations. In this study, hourly data of six air pollutants (e.g., PM2.5, O3, PM10, SO2, NO2, and CO) from monitoring stations between January 2015 and December 2018 are used to feed the proposed deep learning model. Hourly grid-level weather data (e.g., temperature, humidity, air pressure, precipitation, and wind speed at both X- and Y-axis) during the same period are generated by the Weather Research and Forecasting (WRF) model with a grid resolution of 5 km × 5 km, as extra model inputs. For convenience of graph construction, the air pollutant and meteorology datasets are divided into two groups according to seasonal discrepancies: one group contains the summer and autumn data from April to September (higher ozone and lower PM2.5 concentrations) and the other group includes the winter and spring data from October to March (lower ozone and higher PM2.5 concentrations).

2.2. Diffusion Convolution

The Convolutional Neural Network (CNN) is a widely-used network structure that uses the filters containing convolutional kernels to extract spatial features from grid-structured data such as images. With this motivation, a similar idea can be extended into the graph-structured data to extract spatial features from the data and build a model, which is also the essential idea of diffusion convolution.
Diffusion convolution [24] is defined as a combination of the diffusion processes with different steps over the graph. Specifically, the K diffusion steps represent the distance of each node in the graph from the current forecasting position, i.e., how many edges are passed to reach the center point, as shown in Figure 2. For each node, the model calculates the neighbors from 0 to k steps away from itself separately and computes the corresponding transition matrix for each step. The probability θ is a learnable parameter to combine all transition matrices into a diffusion convolution filter when training the model. Here, the diffusion convolution operator ★ B over a graph signal X I × N and the filter f θ are defined as:
X : ,   n B f θ = k = 0 K 1 ( θ k , 1 ( D u p 1 A ) k + θ k , 2 ( D d o w n 1 A T ) k ) X : ,   n     for   n { 0 , , N }
where θ K × 2 is the probability parameter for the filter, D u p and D d o w n represent the in-degree and out-degree diagonal matrix of a graph, and D u p 1 A and D d o w n 1 A T represent the transition matrices of the forward and backward diffusion processes, respectively. Particularly, for undirected graphs, D u p is equal to D d o w n , and D u p 1 A is also equal to D d o w n 1 A T . With the diffusion convolution filter and activation function, the diffusion convolution layer in a neural network can map N -dimensional features to M -dimensional outputs.

2.3. Diffusion Convolutional Recurrent Neural Network (DCRNN)

The DCRNN model captures the spatial dependencies extracted by diffusion convolution and integrates them into the Recurrent Neural Network (RNN) model to handle time series data. The basic principle of RNN is to consider the current inputs as hidden states and process the information from the previous inputs with a multi-gate mechanism [27]. Specifically, the initial GRU sets two internal gated recurrent units to capture the long-term dependencies from time series data. The gate signals of GRU are first computed as follows:
r ( t ) = δ ( W x r X ( t ) + W h r H ( t 1 ) + b r )
z ( t ) = δ ( W x z X ( t ) + W h z H ( t 1 ) + b z )
where r ( t ) is the reset gate and z ( t ) is the update gate at time t ; W x r , W h r , W x z , W h z represent different weight parameters; b r and b z are the biases; the δ denotes logistics sigmoid function. Then, the hidden state H ( t ) at time t is computed as follows:
C ( t ) = t a n h ( W x c X ( t ) + W h c ( r ( t ) H ( t 1 ) ) + b c )
H ( t ) = z ( t ) H ( t 1 ) + ( 1 z ( t ) ) C ( t )
where   C ( t ) represents reset hidden states at time t ; W x c , W h c represent weight parameters; b c is the bias; the operator refers to the Hadamard product of two matrices; t a n h denotes hyperbolic tangent function. Then, the matrix multiplication in the GRU is replaced with the diffusion convolution to build the DCRNN model as follows:
r ( t ) = δ ( Θ r B [ X ( t ) ,   H ( t 1 ) ] + b r )
z ( t ) = δ ( Θ z B [ X ( t ) , H ( t 1 ) ] + b z )
C ( t ) = t a n h ( Θ C B [ X ( t ) , ( r ( t ) H ( t 1 ) ) ] + b c )
H ( t ) = z ( t ) H ( t 1 ) + ( 1 z ( t ) ) C ( t )
where ★ B represents the diffusion convolution as defined in Equation (1) and Θ r , Θ z , Θ C are learnable parameters of the diffusion convolutional filters.
To perform multi-step air quality prediction, the DCRNN model utilizes the Sequence to Sequence (Seq2Seq) architecture, which is a typical Encoder-Decoder architecture based on RNN units [28,29]. In training the DCRNN model, we feed the input sequences (i.e., all the historical features X I × N ) into the encoder and initialize the decoder using the final state of the encoder. Then, the decoder emits the predictions based on the observations. When testing the DCRNN model, we input the test set into the encoder and compare the corresponding prediction results generated by the decoder with the measured data to evaluate the proposed model. Generally, the DCRNN model can achieve accurate air quality predictions by simultaneously capturing the spatial dependencies of topological features among air quality monitoring network and temporal dependencies of multi-source inputs of air quality time series data.

2.4. Graph Construction

Another important step of the DCRNN model is graph construction, which usually reflects the spatial relationships among geospatial data. In this paper, we map the air quality monitoring network with node and edge properties into one graph and calculate the weight matrix among edges over the graph. Generally, one element w i , j of the weight matrix is a reflection of the spatial correlation between nodes v i and v j . Here, we use two types of graph construction methods: undirected and directed graphs. The construction of the undirected graph only considers the geographic distance between two monitoring stations as the spatial relationship:
d i j = d g e o ( ( x i , y i ) , ( x j , y j ) )
W i , j = { exp ( d i j 2 σ 2 ) , d i j < K ( t h r e s h o l d ) 0 ,   o t h e r w i s e
where ( x i , y i ) represents the latitude and longitude coordinates of the node v i , σ and K ( t h r e s h o l d ) are two user-defined hyperparameters.
In terms of the directed graph, we consider the effects of wind direction because wind direction is an important factor that greatly affects the dispersion of air pollutants, as revealed by previous studies [12]. In the undirected graph, W i , j is equal to W j , i , but generally they are not equal in the directed graph due to the consideration of wind factors. Here, we use multiple transformations of the wind direction and geographic distance between nodes to construct different directed graphs and explore their influences on the prediction performances of the DCRNN model. Figure 3 illustrates the relationship between wind direction and directed graph construction, and Equations (12)–(16) clearly present five calculation methods of the edge weight matrix for the construction of different directed graphs as follows:
Directed graph 1:
W i , j = { exp ( d i j 2 sec 2 θ i j σ 2 ) , 90 ° θ i j 90 ° 0 ,   o t h e r w i s e
Directed graph 2:
W i , j = { exp ( d i j 2 sin 2 θ i j σ 2 ) , 90 ° θ i j 90 ° 0 ,   o t h e r w i s e
Directed graph 3:
W i , j = { exp ( d i j 2 ( 1 + sin 2 θ i j ) σ 2 ) , 90 ° θ i j 90 ° 0 ,   o t h e r w i s e
Directed graph 4:
W i , j = exp ( d i j 2 ( 2 cos θ i j ) σ 2 )
Directed graph 5:
W i , j = { exp ( d i j 2 sin θ i j σ 2 ) , 90 ° θ i j 90 ° 0 ,   o t h e r w i s e

2.5. Experimental Design

The experiment is conducted on a server with Ubuntu 16.04 Linux system, 128 GB memory, and NVIDIA Titan RTX (24GB GDDR5 VRAM) graphics card. Python 3.6, Pandas, NumPy, TensorFlow, and Keras are used for data processing and model configuration. In the experiment, the dataset is first divided into two groups according to the seasonal discrepancies for model construction and verification. Roughly 70%, 10%, and 20% of the dataset in each group are separately used for training, validating, and testing the proposed DCRNN model. Specifically, the training set and validation set are used to train the model and evaluate the model during the training process, respectively. The test set is only applied to provide an unbiased evaluation of a final model fit on the training set.
The loss function uses Mean Absolute Error (MAE) and the optimizer adopts adaptive moment estimation (Adam) to minimize the absolute error between the predicted and measured data. Hyperparameters are determined according to the model performances on the validation set. The early stopping technique is used for model training to improve training efficiency and avoid overfitting. Specifically, when the validation error cannot be further improved within the pre-specified number of cycles, the algorithm will terminate early, which can help reduce the computational costs.
In this paper, different statistical indicators including MAE, Root Mean Squared Error (RMSE), and Pearson Correlation Coefficient (r) are used to evaluate the prediction performances of the proposed model, as computed below:
MAE = 1 n i = 1 n | O i P i |
RMSE = 1 n i = 1 n ( O i P i ) 2
r = i = 1 n ( O i O ¯ i ) ( P i P ¯ i ) i = 1 n ( O i O ¯ i ) 2 i = 1 n ( P i P ¯ i ) 2
where O i and O ¯ i , respectively, refer to the observed values and their mean value, and P i and P ¯ i , separately, refer to the prediction value and their mean value.

3. Results and Discussion

3.1. Prediction Performances of DCRNN Using Different Graph Construction Methods

Table 1 and Table 2 present the prediction performances of the DCRNN models using different graph construction methods on the PM2.5 and ozone datasets, respectively. Figure 4 further illustrates the results of Table 1 and Table 2. Overall, the DCRNN models using any graph construction method all exhibit smaller error and higher precision for regional prediction of PM2.5 and ozone concentrations than other baseline models.
As seen in Figure 4a, the PM2.5 prediction results of the DCRNN models using directed graphs and undirected graphs both show smaller MAE than GRU, LSTM, bidirectional LSTM, and Seq2Seq models. Figure 4b indicates that the evaluation metric RMSE presents similar results with MAE in terms of model comparison. Furthermore, the MAE and RMSE from the winter and spring data group are both smaller for the directed graph model compared with the undirected graph model, particularly for the DCRNN models using directed graphs 3, 4, and 5 as shown in Equations (14)–(16). The result is partly attributed to the fact that the construction of directed graphs 3, 4, and 5 all consider the influences of geographic distance and wind direction on the variations in pollutant concentrations. By contrast, the prediction errors of the undirected graph model based on the summer and autumn data group are low enough so that the prediction performance of the directed graph model hardly shows extra superiority to the undirected graph model.
Similarly, the DCRNN models show less errors for regional prediction of ozone than the GRU, LSTM, bidirectional LSTM, and Seq2Seq models, as shown in Figure 4d,e. This result indicates that the DCRNN model can be well applied to regional prediction over a wide range of air pollutants, whether particulate matter or gaseous pollutants. The main difference between PM2.5 and ozone forecasts lies in that the performances of the directed graph model and the undirected graph model depend more on the dataset itself for ozone than PM2.5. Specifically, the directed graph models 3, 4, and 5 present smaller errors for the ozone forecasts on the winter and spring data group than the undirected graph model. However, the directed graph model does not perform significantly better than the undirected graph model on the summer and autumn data group. Overall, Figure 4f shows that directed graph model 5 exhibits the best agreement with the measured values among the DCRNN models with different graph construction methods.
In summary, the DCRNN model using directed graphs outperforms that using undirected graphs in the PM2.5 forecasts based on the winter and spring data group. There exist slight differences for the PM2.5 forecasts on the summer and autumn data group, as well as the ozone forecasts on the two data groups. However, the DCRNN model using directed graph 5 (Gauss Vector Weight), widely used in recent studies [12], generally brings the lowest prediction errors in both PM2.5 and ozone forecasts on all data groups. Therefore, we select the DCRNN model using directed graph 5 as the optimal directed graph model to conduct subsequent model comparison.

3.2. Multi Time-Step Prediction

Table 3 and Table 4 show the performances of the DCRNN model in different time-step predictions separately for PM2.5 and ozone. The prediction errors calculated from the directed graph and undirected graph models both demonstrate an increasing trend with a rise of the prediction time steps. The DCRNN model exhibits the smallest errors between the predicted and measured PM2.5 and ozone data at the 1st hour and shows the largest prediction errors at the 24th hour.
Figure 5 provides a more intuitive visualization of the differences in model performances at different time steps. The DCRNN models using the undirected and directed graphs both show significant differences for the PM2.5 forecasts at different time steps. As shown in Figure 5a,b, the directed graph model on the winter and spring data group has smaller errors at long time steps above 12 h, while the undirected graph model has smaller errors at short time steps below 12 h. However, only when the time step is increased to 24 h, the directed graph model on the summer and autumn data group has similar performances to the undirected graph model.
Figure 5d,e illustrate that the ozone forecasts at different time steps present a similar pattern to the PM2.5 predictions. The directed graph model on the winter and spring data group performs better at long time steps beyond 12 h, while the undirected graph model performs better at short time steps below 12 h. In addition, the undirected graph model on the summer and autumn data group performs better at most time steps except for the 24 h time step. The above results suggest that the performance difference of the DCRNN models using undirected and directed graphs at different time steps could not be solely limited to the PM2.5 and ozone forecasts.
Overall, for short-term forecasts for the 1st to 8th hour, the DCRNN model using the undirected graph exhibits significantly smaller prediction errors than that using directed graphs. However, the discrepancy in prediction errors between the undirected graph and directed graph models gradually decreases with an increase in the prediction time steps. In contrast, for long-term forecasts in the next 12–24 h, the directed graph model has less prediction errors than the undirected graph model for the regional prediction of PM2.5 and ozone on the winter and spring data group. However, the directed graph model does not outperform the undirected graph model in the PM2.5 and ozone forecasts on the summer and autumn data group until the time step was increased beyond 24 h.

3.3. Spatial Distributions of PM2.5 and Ozone Forecasts based on the DCRNN Model

To further evaluate the spatial differences in the model performance, we divide the whole dataset into three groups according to the administrative provinces which the monitoring stations belong to, i.e., Shanghai, Zhejiang, and Jiangsu provinces. Table 5 and Table 6 respectively show the comparison of PM2.5 and ozone concentration prediction of the DCRNN model in the three provinces, and the above results are computed based on 24 h time steps. Figure 6 shows the visualization of the comparison of the above results. Figure 7 presents the results of the spatial distributions of MAE within the study area, using kriging interpolation to calculate the average MAE between predictions and observations in different geographical locations.
In terms of PM2.5 prediction, as shown in Table 5 and Figure 6, the DCRNN models based on directed and undirected graphs show a similar trend, with the lowest prediction error in Zhejiang Province, followed by Shanghai, and the largest prediction error in Jiangsu Province. In Figure 7a,b,e,f, the PM2.5 prediction errors of the undirected and directed graphs show the similar spatial distributions in the same seasons, but the northern province notably presents more prediction errors than the south. The larger prediction error in the northern province, i.e., Jiangsu province, could be related to two factors: (1) The relatively sparse distribution of air quality monitoring stations in the northern part of Jiangsu province increases the difficulty of providing refined data support for accurate air quality prediction. (2) Northern regions are usually accompanied by heavier particulate matter pollution and more variable PM2.5 concentrations, thus leading to larger PM2.5 forecast errors especially on heavily polluted days.
In terms of ozone prediction, as shown in Table 6 and Figure 6, the ozone prediction errors of the directed graph model are smaller than that of the undirected graph model in Jiangsu province, which indicates that the ozone concentration variation in Jiangsu province could be closely related to wind factors and transport. In Shanghai, an opposite trend of ozone prediction error happens. The undirected graph model performs significantly better than the directed graph model for summer ozone concentration prediction. The results suggest that summer ozone pollution in Shanghai could be more related to local source emissions. In Figure 7c,d,g,h, the spatial difference of ozone prediction errors varies less than PM2.5 at the regional level, with smaller differences (MAE and RMSE) in Zhejiang, Shanghai, and Jiangsu provinces. Meanwhile, compared with the PM2.5 prediction results, the directed and undirected graph models both show stronger spatial variability in the ozone prediction errors.

4. Conclusions

In this study, we employ different construction methods of directed and undirected graphs to establish a novel diffusion convolutional recurrent neural network (DCRNN) model for the regional prediction of PM2.5 and ozone concentrations. The model can fully consider the spatial relationships between nodes within air quality monitoring network by integrating the combined effects of station-level geographic distance and dominant wind direction in the study area. Then, hourly PM2.5 and ozone data collected from 123 air quality monitoring stations in the Yangtze River Delta region are used to train, validate, and test the proposed DCRNN model. Several meaningful findings are summarized as follows:
(1)
The DCRNN model outperforms the baseline models (e.g., GRU and LSTM) in PM2.5 and ozone forecasts.
(2)
The DCRNN model using directed graphs with an integration of wind factors outperforms the undirected graph model in the long-term prediction of PM2.5 and ozone.
(3)
The undirected graph model could achieve better performance in the short-term forecasts, particularly for the next 1st hour prediction.
(4)
The prediction errors of the DCRNN model using undirected and directed graphs both suggest an upward trend with an increase in the prediction time steps, particularly for the undirected graph model.
(5)
The monitoring stations that are sparsely distributed or located in heavily polluted areas could both cause lower prediction accuracy.
In terms of applications, the proposed model could assist environmental researchers in further improving the technologies of air quality prediction and serve as tools for environmental policymakers to implement related pollution control policies. The comparison results between the directed and undirected graph-based models for specific regions (e.g., provinces) and the inferences about the effects of wind factors on pollution in different regions could provide decision support for accurate pollution control.
One limitation of this study is that the DCRNN model restricts the dynamic characterization of the wind direction factor and just uses the weighted-average vector representing wind direction. In future studies, we could consider more advanced spatiotemporal prediction methods based on dynamic graph structures to model the dynamic effects of wind direction and further strengthen air quality prediction.

Author Contributions

Conceptualization, D.W., H.-W.W., K.-F.L., and Z.-R.P.; methodology, D.W. and H.-W.W.; formal analysis, D.W.; software, D.W., H.-W.W. and J.Z.; data curation, D.W. and J.Z.; writing—original draft preparation, D.W.; writing—review and editing, D.W., H.-W.W., K.-F.L. and Z.-R.P.; visualization, D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Planning Office of Philosophy and Social Science (No. 16ZDA048).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available from the corresponding author on reasonable request.

Acknowledgments

The authors thank Hong-Di He and Wan-Jin Cai for their valuable advice on this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chan, C.K.; Yao, X. Air Pollution in Mega Cities in China. Atmos. Environ. 2008, 42, 1–42. [Google Scholar] [CrossRef]
  2. Gao, J.; Woodward, A.; Vardoulakis, S.; Kovats, S.; Wilkinson, P.; Li, L.; Xu, L.; Li, J.; Yang, J.; Li, J.; et al. Haze, Public Health and Mitigation Measures in China: A Review of the Current Evidence for Further Policy Response. Sci. Total Environ. 2017, 578, 148–157. [Google Scholar] [CrossRef] [PubMed]
  3. Laden, F.; Schwartz, J.; Speizer, F.E.; Dockery, D.W. Reduction in Fine Particulate Air Pollution and Mortality: Extended Follow-up of the Harvard Six Cities Study. Am. J. Respir. Crit. Care Med. 2006, 173, 667–672. [Google Scholar] [CrossRef] [PubMed]
  4. Song, Y.; Wang, X.; Maher, B.A.; Li, F.; Xu, C.; Liu, X.; Sun, X.; Zhang, Z. The Spatial-Temporal Characteristics and Health Impacts of Ambient Fine Particulate Matter in China. J. Clean. Prod. 2016, 112, 1312–1318. [Google Scholar] [CrossRef]
  5. Pui, D.Y.H.; Chen, S.C.; Zuo, Z. PM2.5 in China: Measurements, Sources, Visibility and Health Effects, and Mitigation. Particuology 2014, 13, 1–26. [Google Scholar] [CrossRef]
  6. Brauer, M.; Freedman, G.; Frostad, J.; Van Donkelaar, A.; Martin, R.V.; Dentener, F.; Dingenen, R.V.; Estep, K.; Amini, H.; Apte, J.S.; et al. Ambient Air Pollution Exposure Estimation for the Global Burden of Disease 2013. Environ. Sci. Technol. 2016, 50, 79–88. [Google Scholar] [CrossRef]
  7. Wang, T.; Xue, L.; Brimblecombe, P.; Lam, Y.F.; Li, L.; Zhang, L. Ozone Pollution in China: A Review of Concentrations, Meteorological Influences, Chemical Precursors, and Effects. Sci. Total Environ. 2017, 575, 1582–1596. [Google Scholar] [CrossRef] [PubMed]
  8. Skamarock, W.C.; Klemp, J.; Dudhia, J.; Gill, D.O.; Barker, D.; Wang, W.; Powers, J.G. A Description of the Advanced Research WRF Version 3. Univ. Corp. Atmos. Res. 2008, 27, 3–27. [Google Scholar]
  9. Grell, G.A.; Peckham, S.E.; Schmitz, R.; McKeen, S.A.; Frost, G.; Skamarock, W.C.; Eder, B. Fully Coupled “Online” Chemistry within the WRF Model. Atmos. Environ. 2005, 39, 6957–6975. [Google Scholar] [CrossRef]
  10. Byun, D.W.; Schere, K.L. Review of the Governing Equations, Computational Algorithms, and Other Components of the Models-3 Community Multiscale Air Quality (CMAQ) Modeling System. Appl. Mech. Rev. 2006, 59, 51–77. [Google Scholar] [CrossRef]
  11. Ni, X.Y.; Huang, H.; Du, W.P. Relevance Analysis and Short-Term Prediction of PM2.5concentrations in Beijing Based on Multi-Source Data. Atmos. Environ. 2017, 150, 146–161. [Google Scholar] [CrossRef]
  12. Yang, W.; Deng, M.; Xu, F.; Wang, H. Prediction of Hourly PM2.5 Using a Space-Time Support Vector Regression Model. Atmos. Environ. 2018, 181, 12–19. [Google Scholar] [CrossRef]
  13. Shang, Z.; Deng, T.; He, J.; Duan, X. A Novel Model for Hourly PM2.5 Concentration Prediction Based on CART and EELM. Sci. Total Environ. 2019, 651, 3043–3052. [Google Scholar] [CrossRef] [PubMed]
  14. Feng, R.; Zheng, H.J.; Gao, H.; Zhang, A.; Huang, C.; Zhang, J.; Luo, K.; Fan, J. Recurrent Neural Network and Random Forest for Analysis and Accurate Forecast of Atmospheric Pollutants: A Case Study in Hangzhou, China. J. Clean. Prod. 2019, 231, 1005–1015. [Google Scholar] [CrossRef]
  15. Li, X.; Peng, L.; Yao, X.; Cui, S.; Hu, Y.; You, C.; Chi, T. Long Short-Term Memory Neural Network for Air Pollutant Concentration Predictions: Method Development and Evaluation. Environ. Pollut. 2017, 231, 997–1004. [Google Scholar] [CrossRef] [PubMed]
  16. Zhou, Y.; Chang, F.; Chang, L.; Kao, I.; Wang, Y. Explore a Deep Learning Multi-Output Neural Network for Regional Multi-Step-Ahead Air Quality Forecasts. J. Clean. Prod. 2019, 209, 134–145. [Google Scholar] [CrossRef]
  17. Huang, Y.; Shen, L.; Liu, H. Grey Relational Analysis, Principal Component Analysis and Forecasting of Carbon Emissions Based on Long Short-Term Memory in China. J. Clean. Prod. 2019, 209, 415–423. [Google Scholar] [CrossRef]
  18. Li, L.L.; Wen, S.Y.; Tseng, M.L.; Wang, C.S. Renewable Energy Prediction: A Novel Short-Term Prediction Model of Photovoltaic Output Power. J. Clean. Prod. 2019, 228, 359–375. [Google Scholar] [CrossRef]
  19. Chang, Y.S.; Chiao, H.T.; Abimannan, S.; Huang, Y.P.; Tsai, Y.T.; Lin, K.M. An LSTM-Based Aggregated Model for Air Pollution Forecasting. Atmos. Pollut. Res. 2020, 11, 1451–1463. [Google Scholar] [CrossRef]
  20. Athira, V.; Geetha, P.; Vinayakumar, R.; Soman, K.P. DeepAirNet: Applying Recurrent Networks for Air Quality Prediction. Procedia Comput. Sci. 2018, 132, 1394–1403. [Google Scholar] [CrossRef]
  21. Zhang, K.; Thé, J.; Xie, G.; Yu, H. Multi-Step Ahead Forecasting of Regional Air Quality Using Spatial-Temporal Deep Neural Networks: A Case Study of Huaihai Economic Zone. J. Clean. Prod. 2020, 277, 123231. [Google Scholar] [CrossRef]
  22. Wen, C.; Liu, S.; Yao, X.; Peng, L.; Li, X.; Hu, Y.; Chi, T. A Novel Spatiotemporal Convolutional Long Short-Term Neural Network for Air Pollution Prediction. Sci. Total Environ. 2019, 654, 1091–1099. [Google Scholar] [CrossRef] [PubMed]
  23. Shi, X.; Chen, Z.; Wang, H.; Yeung, D.-Y.; Wong, W.; Woo, W. Convolutional LSTM Network: A Machine Learning Approach for Precipitation Nowcasting. Adv. Neural Inf. Processing Syst. 2015, 1, 802–810. [Google Scholar]
  24. Lin, Y.; Mago, N.; Gao, Y.; Li, Y.; Chiang, Y.Y.; Shahabi, C.; Ambite, J.L. Exploiting Spatiotemporal Patterns for Accurate Air Quality Forecasting Using Deep Learning. In Proceedings of the 26th ACM SIGSPATIAL International Conference on Advances in Geographic Information Systems, Seattle, WA, USA, 6–9 November 2018; pp. 359–368. [Google Scholar] [CrossRef]
  25. Ouyang, X.; Yang, Y.; Zhang, Y.; Zhou, W. Spatial-Temporal Dynamic Graph Convolution Neural Network for Air Quality Prediction. In Proceedings of the 2021 International Joint Conference on Neural Networks, Shenzhen, China, 18–22 July 2021; pp. 1–8. [Google Scholar] [CrossRef]
  26. Qi, Y.; Li, Q.; Karimian, H.; Liu, D. A Hybrid Model for Spatiotemporal Forecasting of PM2.5 Based on Graph Convolutional Neural Network and Long Short-Term Memory. Sci. Total Environ. 2019, 664, 1–10. [Google Scholar] [CrossRef] [PubMed]
  27. Li, Y.; Yu, R.; Shahabi, C.; Liu, Y. Diffusion Convolutional Recurrent Neural Network: Data-Driven Traffic Forecasting. In 6th International Conference on Learning Representations, Proceedings of the ICLR 2018—Conference Track Proceedings, Vancouver, BC, Canada, 30 April–3 May 2018; pp. 1–16.
  28. Cho, K.; van Merrienboer, B.; Gulcehre, C.; Bahdanau, D.; Bougares, F.; Schwenk, H.; Bengio, Y. Learning Phrase Representatxions Using RNN Encoder–Decoder for Statistical Machine Translation. In Proceedings of the 2014 Conference on Empirical Methods in Natural Language Processing (EMNLP), Doha, Qatar, 25–29 October 2014; pp. 1724–1734. [Google Scholar] [CrossRef]
  29. Sutskever, I.; Vinyals, O.; Le, Q.V. Sequence to Sequence Learning with Neural Networks. Adv. Neural Inf. Processing Syst. 2014, 2, 3104–3112. [Google Scholar]
Figure 1. Geographical location of the study area and spatial distributions of 123 air quality monitoring stations.
Figure 1. Geographical location of the study area and spatial distributions of 123 air quality monitoring stations.
Ijerph 19 03988 g001
Figure 2. Illustration of the diffusion convolution process with K diffusion steps.
Figure 2. Illustration of the diffusion convolution process with K diffusion steps.
Ijerph 19 03988 g002
Figure 3. Spatial relationship between wind direction and directed graph construction. Wind vector shows the vector weighted average wind direction of each data group (i.e., Winter and Spring, Summer and Autumn).
Figure 3. Spatial relationship between wind direction and directed graph construction. Wind vector shows the vector weighted average wind direction of each data group (i.e., Winter and Spring, Summer and Autumn).
Ijerph 19 03988 g003
Figure 4. Model comparison between the DCRNN model and baseline models. (a) PM2.5, MAE; (b) PM2.5, RMSE; (c) PM2.5, r; (d) ozone, MAE; (e) ozone, RMSE; (f) ozone, r.
Figure 4. Model comparison between the DCRNN model and baseline models. (a) PM2.5, MAE; (b) PM2.5, RMSE; (c) PM2.5, r; (d) ozone, MAE; (e) ozone, RMSE; (f) ozone, r.
Ijerph 19 03988 g004
Figure 5. Regional prediction of PM2.5 and ozone based on the DCRNN model at different time steps. (a) PM2.5, MAE; (b) PM2.5, RMSE; (c) PM2.5, r; (d) ozone, MAE; (e) ozone, RMSE; (f) ozone, r.
Figure 5. Regional prediction of PM2.5 and ozone based on the DCRNN model at different time steps. (a) PM2.5, MAE; (b) PM2.5, RMSE; (c) PM2.5, r; (d) ozone, MAE; (e) ozone, RMSE; (f) ozone, r.
Ijerph 19 03988 g005
Figure 6. Model comparison between the DCRNN model and baseline models. (a) PM2.5, MAE; (b) PM2.5, RMSE; (c) ozone, MAE; (d) ozone, RMSE.
Figure 6. Model comparison between the DCRNN model and baseline models. (a) PM2.5, MAE; (b) PM2.5, RMSE; (c) ozone, MAE; (d) ozone, RMSE.
Ijerph 19 03988 g006
Figure 7. Spatial distributions of the mean MAE of 24 h air quality predictions based on the DCRNN model in the Yangtze River Delta region: (a) undirected graph, PM2.5 in winter and spring; (b) undirected graph, PM2.5 in summer and autumn; (c) undirected graph, ozone in winter and spring; (d) undirected graph, ozone in summer and autumn; (e) directed graph, PM2.5 in winter and spring; (f) directed graph, PM2.5 in summer and autumn; (g) directed graph, ozone in winter and spring; (h) directed graph, ozone in summer and autumn.
Figure 7. Spatial distributions of the mean MAE of 24 h air quality predictions based on the DCRNN model in the Yangtze River Delta region: (a) undirected graph, PM2.5 in winter and spring; (b) undirected graph, PM2.5 in summer and autumn; (c) undirected graph, ozone in winter and spring; (d) undirected graph, ozone in summer and autumn; (e) directed graph, PM2.5 in winter and spring; (f) directed graph, PM2.5 in summer and autumn; (g) directed graph, ozone in winter and spring; (h) directed graph, ozone in summer and autumn.
Ijerph 19 03988 g007
Table 1. Model comparison between the DCRNN model and baseline models in the PM2.5 forecasts.
Table 1. Model comparison between the DCRNN model and baseline models in the PM2.5 forecasts.
ModelWinter and SpringSummer and Autumn
MAE (μg/m3)RMSE (μg/m3)rMAE (μg/m3)RMSE (μg/m3)r
GRU23.0733.350.6211.6916.360.53
LSTM22.8232.930.6312.1516.600.53
Bidirectional LSTM19.7928.980.7310.3114.710.64
Seq2seq18.5028.710.759.8414.440.65
DCRNN (undirected graph)18.0530.110.798.7612.920.73
DCRNN (directed graph 1)17.0126.230.798.9513.240.72
DCRNN (directed graph 2)17.7327.230.799.0813.330.72
DCRNN (directed graph 3)16.8225.730.808.9213.100.72
DCRNN (directed graph 4)16.3725.130.818.9213.090.73
DCRNN (directed graph 5)17.2025.740.808.8513.110.73
Table 2. Model comparison between the DCRNN model and baseline models in the ozone forecasts.
Table 2. Model comparison between the DCRNN model and baseline models in the ozone forecasts.
ModelWinter and SpringSummer and Autumn
MAE (μg/m3)RMSE (μg/m3)rMAE (μg/m3)RMSE (μg/m3)r
GRU21.6028.120.7028.1837.400.72
LSTM21.8428.490.6828.4437.890.71
Bidirectional LSTM20.0326.590.7226.8336.020.74
Seq2seq19.8726.820.7025.3534.590.75
DCRNN (undirected graph)18.3025.110.7622.9532.440.78
DCRNN (directed graph 1)18.8526.080.7522.7132.070.80
DCRNN (directed graph 2)18.4525.210.7623.9433.720.78
DCRNN (directed graph 3)17.7424.340.7723.3433.060.79
DCRNN (directed graph 4)17.9924.610.7723.4032.970.79
DCRNN (directed graph 5)17.9224.530.7723.0032.240.80
Table 3. The performance of the DCRNN model in forecasting PM2.5 concentrations at different time steps.
Table 3. The performance of the DCRNN model in forecasting PM2.5 concentrations at different time steps.
Data GroupTime-StepUndirected GraphDirected Graph
MAE (μg/m3)RMSE (μg/m3)rMAE (μg/m3)RMSE (μg/m3)r
Winter
and Spring
1 h5.178.000.986.7510.160.97
2 h6.4010.160.977.8011.850.96
4 h8.2613.200.959.4814.530.94
8 h11.0317.400.9211.8018.020.91
12 h13.3321.210.8813.5720.580.88
24 h18.0530.110.7817.2025.740.80
Summer
and Autumn
1 h3.515.540.964.226.370.95
2 h4.236.770.944.867.470.92
4 h5.228.260.915.808.920.89
8 h6.4310.170.876.9210.740.85
12 h7.2911.250.847.6411.660.82
24 h8.7712.920.738.8513.110.73
Table 4. The performance of the DCRNN model in forecasting ozone concentrations at different time steps.
Table 4. The performance of the DCRNN model in forecasting ozone concentrations at different time steps.
Data GroupTime-StepUndirected GraphDirected Graph
MAE (μg/m3)RMSE (μg/m3)rMAE (μg/m3)RMSE (μg/m3)r
Winter
and Spring
1 h6.259.350.957.2010.650.93
2 h7.6411.410.928.4112.340.90
4 h9.4713.770.8810.0414.420.86
8 h11.3315.930.8311.6416.280.82
12 h12.5117.350.8012.6717.420.80
24 h18.3025.110.7517.9224.530.77
Summer
and Autumn
1 h7.2911.030.958.6812.740.93
2 h8.7913.120.9210.0114.530.90
4 h10.6115.420.8811.6616.610.86
8 h12.3817.780.8213.2518.790.80
12 h14.7121.020.8215.3121.580.81
24 h22.9532.440.7823.0032.240.80
Table 5. The performance of the DCRNN model in forecasting PM2.5 concentrations in different geographical areas.
Table 5. The performance of the DCRNN model in forecasting PM2.5 concentrations in different geographical areas.
Data GroupProvinceUndirected GraphDirected Graph
MAE (μg/m3)RMSE (μg/m3)rMAE (μg/m3)RMSE (μg/m3)r
Winter
and Spring
Shanghai16.7928.630.7615.4423.100.77
Jiangsu20.6233.070.7919.7029.510.75
Zhejiang15.7827.220.7515.0621.970.79
Summer
and Autumn
Shanghai8.5813.430.789.0113.170.79
Jiangsu9.9114.130.729.9414.330.71
Zhejiang7.6811.670.727.8011.830.71
Table 6. The performance of the DCRNN model in forecasting ozone concentrations in different geographical areas.
Table 6. The performance of the DCRNN model in forecasting ozone concentrations in different geographical areas.
Data GroupProvinceUndirected GraphDirected Graph
MAE (μg/m3)RMSE (μg/m3)rMAE (μg/m3)RMSE (μg/m3)r
Winter
and Spring
Shanghai19.3325.810.6919.0025.090.72
Jiangsu18.3625.100.7517.2523.510.79
Zhejiang18.1125.030.7618.4125.380.76
Summer
and Autumn
Shanghai22.1632.000.7724.2533.880.75
Jiangsu23.8233.460.7823.6632.930.80
Zhejiang22.2331.500.7922.2031.340.79
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, D.; Wang, H.-W.; Lu, K.-F.; Peng, Z.-R.; Zhao, J. Regional Prediction of Ozone and Fine Particulate Matter Using Diffusion Convolutional Recurrent Neural Network. Int. J. Environ. Res. Public Health 2022, 19, 3988. https://doi.org/10.3390/ijerph19073988

AMA Style

Wang D, Wang H-W, Lu K-F, Peng Z-R, Zhao J. Regional Prediction of Ozone and Fine Particulate Matter Using Diffusion Convolutional Recurrent Neural Network. International Journal of Environmental Research and Public Health. 2022; 19(7):3988. https://doi.org/10.3390/ijerph19073988

Chicago/Turabian Style

Wang, Dongsheng, Hong-Wei Wang, Kai-Fa Lu, Zhong-Ren Peng, and Juanhao Zhao. 2022. "Regional Prediction of Ozone and Fine Particulate Matter Using Diffusion Convolutional Recurrent Neural Network" International Journal of Environmental Research and Public Health 19, no. 7: 3988. https://doi.org/10.3390/ijerph19073988

APA Style

Wang, D., Wang, H. -W., Lu, K. -F., Peng, Z. -R., & Zhao, J. (2022). Regional Prediction of Ozone and Fine Particulate Matter Using Diffusion Convolutional Recurrent Neural Network. International Journal of Environmental Research and Public Health, 19(7), 3988. https://doi.org/10.3390/ijerph19073988

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop