Next Article in Journal
Ensemble Classification for Anomalous Propagation Echo Detection with Clustering-Based Subset-Selection Method
Previous Article in Journal
Properties of the Long-Term Oscillations in the Middle Atmosphere Based on Observations from TIMED/SABER Instrument and FPI over Kelan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting PM2.5 Concentrations at a Regional Background Station Using Second Order Self-Organizing Fuzzy Neural Network

1
Faculty of Information Technology, Beijing University of Technology, Beijing 100124, China
2
Beijing Key Laboratory of Computational Intelligence and Intelligent System, Beijing University of Technology, Beijing 100124, China
*
Author to whom correspondence should be addressed.
Atmosphere 2017, 8(1), 10; https://doi.org/10.3390/atmos8010010
Submission received: 23 October 2016 / Revised: 3 January 2017 / Accepted: 5 January 2017 / Published: 12 January 2017

Abstract

:
This study aims to develop a second order self-organizing fuzzy neural network (SOFNN) to predict the hourly concentrations of fine particulate matter (PM2.5) for the next 24 h at a regional background station called Shangdianzi (SDZ) in China from 14 to 23 January 2010. The structure of the SOFNN was automatically adjusted according to the sensitivity analysis (SA) of model output and the parameter-learning phase was performed applying a second order gradient (SOG) algorithm. Principal component analysis (PCA) was employed to select the dominating factors for PM2.5 concentrations as the input variables for the SOFNN. It was found that the dominating variables (relative humidity (RH), pressure (Pre), aerosol optical depth (AOD), wind speed (WS) and wind direction (WD)) extracted by PCA agreed well with the characteristics of PM2.5 at SDZ where the PM2.5 concentrations were heavily affected by meteorological parameters and were closely related to AOD. The forecasting results showed that the proposed SOG-SASOFNN performed better than other models with higher coefficient of determination (R2) during both training phase and test phase (0.89 and 0.84, respectively) in predicting PM2.5 concentrations at SDZ. In conclusion, the developed SOG-SASOFNN provided satisfying results for modeling the hourly distribution of PM2.5 at SDZ during the studied period.

1. Introduction

China has suffered from conspicuously rising air pollution due to the rapid urbanization since 1980s. Gaseous and particulate emissions from various anthropogenic sources such as coal combustion, motor vehicles and industrial processes lead to poor air quality (AQ) [1]. One of the most harmful pollutants is known as fine particulate matter (PM2.5), a complex mixture of particles with aerodynamic diameters of 2.5 μm or less [2]. Chronic exposure to ambient PM2.5 causes the public to be faced with increased morbidity and mortality as it can lodge deeply into the lungs [3,4]. In Atlanta, a 10 μg/m3 increase of PM2.5 is associated with 3.3% increases in emergency department visits for cardiovascular disease [5] and 1.6% increases for respiratory disease [6]. In addition, aerosols affect climate change in a way that a dramatic emission reduction (35%–80%) in anthropogenic aerosols results in approximately 1 °C of extra warming and approximately 0.1 mm·day−1 of extra precipitation globally averaged by the end of the 21st century [7]. Moreover, visibility and PM2.5 are negatively correlated with Pearson correlations of −0.663 for a regional background station in the boundary of North China Plain (NCP) [8]. In order to take preventive and regulatory measures to manage AQ as well as to provide specific information about concentrations of PM2.5 with the public, it is necessary to develop methods that can predict PM2.5 concentrations precisely.
Modeling of the real world processes such as forecasting of AQ has proven to be tough due to their highly chaotic and nonlinear phenomenon [9,10]. Generally, models used to estimate air pollutants comprise process-based methods and data-driven statistical methods [11]. The first class of methods is based on chemical transport models (CTMs) which model the atmospheric and chemical processes involved in the production of an air pollutant [12]. Yu et al. [13] evaluated the performance of the Eta-Community Multiscale Air Quality (CMAQ) model in predicting PM2.5 and its chemical species over the eastern United States, where the model captured a majority (73%) of PM2.5 observations within a factor of 2. However, CTMs are complex to implement as they require detailed information on meteorology, emissions inventories and advanced understanding of chemical process. On the contrary, the latter class of methods is driven by satellite readings or ground monitoring measurements [14], thus there is no need for strong assumptions of the underlying air pollution processes during the modeling process. Among the statistical AQ forecasting models, we can mention those based on linear regressions (LR) and neural networks (NNs). Regression models show optimal results when relationship between the predictors (such as meteorological and AQ variables) and the predictand (pollutant concentrations of interest) are almost linear [15]. Regression models are also likely to under predict high concentrations and over predict low concentrations. To overcome some of these limitations, the artificial neural network (ANN), which can be used to derive non-linear functions relating to the predictand and the predictors, has been adopted and proven to be a reliable air pollution modeling tool in many studies [16,17,18]. Chaloulakou et al. [19] analyzed the performance of multiple linear regressions (MLR) and NNs on the forecasting of PM10 in Athens. Ordieres et al. [20] implemented three types of NNs as well as a LR model and a persistence model to forecast daily averages of PM2.5 in El Paso and Ciudad Juaŕez. The results of both studies clearly demonstrated that NNs were more accurate than LR models for forecasting of concentrations of air pollutants.
Unfortunately, the mapping rules in NNs are not visible and are difficult to understand. That is, one can only ascertain relationships among meteorology and pollutants concentrations data but never be sure about underlying causal mechanism included in the emission and dispersion pattern of air pollutants through NNs. In contrast, fuzzy system (FS) technologies often deal with issues such as reasoning on a high-level than NNs and can handle imprecise information through linguistic expressions [21,22]. Nevertheless, it is difficult to tune the fuzzy rules from the training data for a human operator since FS has limited learning capability. However, the combination of neural network (NN) and FS which leads to fuzzy neural network (FNN) can strength the prediction capabilities compared to using a single methodology [23]. Mishra et al. [24] developed an artificial intelligence based neuro-fuzzy (NF) model to forecast PM2.5 concentrations during haze hours in Delhi along with MLR and ANN, where NF model performed better in prediction capabilities compared to MLR and ANN.
For fuzzy neural networks (FNNs), the main training process consists of parameter-learning phase and structure-learning phase [25]. However, in most FNNs, only the learning phase of the parameters is determined using supervised or unsupervised learning algorithm and the structure identification is still difficult [26,27]. It is hard for designers to choose the appropriate number of fuzzy rules to forecast air pollutants, especially PM2.5, applying FNNs with empirically fixed structure. Several early research works have already dealt with self-organizing FNNs (SOFNNs), of which the structure can be adjusted automatically during training process [28]. Habbi et al. [29] used artificial bee colony optimization strategy to find the structure and the parameters of TS fuzzy systems simultaneously. Han et al. [30] developed a sensitivity analysis (SA) of model output based growing-and-pruning approach to generate fuzzy neural model with highly accurate and compact structure. However, as far as we know, few SOFNNs have been applied to estimate the distributions of air pollutants or to forecast the concentrations of PM2.5.
Therefore, in this paper, the SA based SOFNN (SASOFNN) was used for forecasting the hourly PM2.5 concentrations of the next 24 h at Shangdianzi (SDZ) station, which is a regional background station located in the boundary of the NCP, from 14 to 23 January 2010. The parameters of the SASOFNN were optimized by a second order gradient (SOG) algorithm with improved training speed and higher approximation accuracy [31,32], as the first order gradient (FOG) algorithm having limited search ability and the genetic algorithm (GA) being very time and computation expensive. Besides, in order to eliminate the irrelevant input parameters, the principal component analysis (PCA) was used to select dominating variables most correlated to PM2.5 from temperature (T), relative humidity (RH), wind speed (WS), wind direction (WD), pressure (Pre), visibility (Vis), aerosol optical depth (AOD), CO, NO2, O3 and SO2 as the input parameters for the SOG-SASOFNN [33]. It could be found that the dominating variables (RH, Pre, AOD, WS and WD) extracted by PCA coincided with the characteristics of PM2.5 at SDZ where there were no significant pollution sources within 30 km of the site and the PM2.5 concentrations were severely influenced by the meteorological parameters and were tightly associated with the AOD. In addition, the prediction performance of the SOG-SASOFNN was analyzed and compared to that of the SASOFNN with FOG algorithm (FOG-SASOFNN) and the echo state network (ESN) in the study of Xu et al. [11] as well as to the Eta-CMAQ modeling results for the hourly PM2.5 concentrations at the rural sites over the eastern United States from 15 July to 19 August 2004 [13]. The statistical parameters indexes calculated for the models above demonstrated that the SOG-SASOFNN performed better than the FOG-SASOFNN, the ESN and the Eta-CMAQ in predicting the hourly concentrations of PM2.5. It can be concluded that the SOG-SASOFNN presented here is valid for estimating the hourly distribution of PM2.5 24 h ahead at SDZ during the studied period.

2. Study Site and Data

2.1. Study Site

As shown in Figure 1, the SDZ (40°39’ N, 117°07’ E, 293.9 m a.s.l.) site selected for applying the proposed methodology is located in the northern part of NCP and is 100 km northeast of the urban area of Beijing. It is one of the regional Global Atmosphere Watch (GAW) stations in China. There are no densely populated and industrial areas within a distance of 30 km around the station, so the atmospheric pollution level at SDZ station represents the background concentration of atmospheric pollutants in the economically developed regions of north China. The prevailing winds at SDZ are influenced by the valley topography and are from east-northeast and west-southwest. When southwesterly winds arise, the polluted air masses from urban areas and satellite towns of Beijing can be transported to SDZ while relatively clean air masses arrive from other wind directions. There has been detailed information regarding SDZ site presented in previous studies [34,35].

2.2. Data Preparation

The data used in this paper have the time frame of 14 to 23 January 2010 during which a regional haze episode occurred in the Beijing, Tianjin and Hebei Provinces (BTH) area in the NCP from 16 to 19 January 2010. At SDZ station, concentrations of CO, NO2, O3, SO2 and PM2.5 were observed by a TE48C CO analyzer, TE 42CTL NOX analyzer, TE49C O3 analyzer, TE43C SO2 analyzer and a tapered element oscillating microbalance (TEOM 1400a), respectively. An automatic weather station installed at the SDZ meteorological station measured the hourly meteorological data including T, RH, WS, WD, Pre and Vis. The AOD retrieved with the lidar observation was also presented in the study of Zhao et al. [36]. The meteorological variables and air pollutants except PM2.5 as well as AOD constitute the predictors and the PM2.5 is the predictand. The entire data table of the predictors and predictand consists of 240 h of measurements and was uploaded as non-published material during the submission process. Considering that the values of AOD were mostly missed on 19 January due to the presence of cumulus clouds at the top of the planetary boundary layer (PBL), the rows from 13:00 18 January to 23:00 19 January were eliminated. The hourly measurements from 1:00 23 January to 23:00 23 January 2010 were also deleted owing to the missing corresponding predictand values. Finally, only 182 h of measurements were chosen for the present study. The units of most variables are summarized in Table 1 and the AOD is a non unit variable.

3. Theory and Methodology

PM2.5 concentrations are affected by emission sources, meteorological conditions and topographical characteristics of the area under study, which makes estimating the PM2.5 distribution difficult. In order to model the hourly concentrations of PM2.5 at SDZ precisely, the SOG-SASOFNN was developed to forecast the PM2.5 concentrations of the next 24 h employing the current hour values of the dominating parameters extracted by PCA.

3.1. Principal Component Analysis

PCA is a classical statistical technique that analyzes the covariance structure of multivariate variables and has been successfully applied for several tasks in AQ domain [37,38]. In this paper, PCA was used to provide the interdependencies of the data measured at SDZ station. PCA selected the dominating variables having the greatest impact on PM2.5 in the following steps:
  • Let 182 × 11 data matrix X = [ x 1 , x 2 , ... , x 11 ] denotes the 182 h of measurements of predictors and x 1 , x 2 , ... , x 11 are 182 × 1 data vectors of T, RH, WS, WD, Pre, Vis, AOD, CO, NO2, O3 and SO2, respectively. The data matrix should be transformed into a standardized form:
    Z = ( z i j ) 182 × 11 = ( x i j x ¯ j δ j ) 182 × 11 , i = 1 , 2 , ... , 182 ; j = 1 , 2 , ... , 11
    where Z is the standardized data matrix generated from X ; x i j and z i j are the values of predictor j in sample i before and after the standardization; and x ¯ j and δ j are the arithmetic mean value and the standard deviation for predictor j, respectively.
  • Calculate the correlation coefficient matrix using Equation (2):
    R = 1 182 Z T Z
  • Compute the eigenvalues λ 1 , λ 2 , ... , λ 11 and the corresponding eigenvectors γ 1 , γ 2 , ... , γ 11 of 11 × 11 correlation matrix R .
  • Reorder the eigenvalues in descending order to bring λ 1 ' > λ 2 ' > ... > λ 11 ' and readjust the eigenvectors as γ 1 ' , γ 2 ' , ... , γ 11 ' accordingly.
  • Obtain the unit orthogonal eigenvectors ρ 1 , ρ 2 , ... , ρ 11 using the Schmidt orthogonal method on γ 1 ' , γ 2 ' , ... , γ 11 ' .
  • Calculate the cumulative contribution rate θ 1 , θ 2 , ... , θ 11 of the eigenvalues λ 1 ' , λ 2 ' , ... , λ 11 ' and α variables will be extracted if θ α θ where θ is the preset extraction efficiency.
  • The data of the dominating variables is acquired by computing the projection of Z on the extracted unit orthogonal eigenvectors using Equation (3):
    Y = Z ρ
    where ρ = [ ρ 1 , ρ 2 , ... , ρ α ] .

3.2. Sensitivity Analysis Based Self-Organizing Fuzzy Neural Network with Second Order Gradient Algorithm

3.2.1. Architecture of the Proposed Model

The SOG-SASOFNN used to search a suitable nonlinear mapping between the dominating variables and the hourly concentrations of PM2.5 in the present study is based on radial basis function (RBF) neurons. As shown in Figure 2, the initial architecture of the SOG-SASOFNN has four layers, the input layer, the RBF layer, the normalized layer and the output layer.
Mathematically, each layer in the initial SOG-SASOFNN is described as follows:
  • Input layer: There are N neurons in this layer and the output value of the ith neuron can be expressed as follows:
    u i = r i , i = 1 , 2 , ... , N
    where r = [ r 1 , r 2 , ... , r N ] represents the dominating variables extracted from the predictors through PCA method.
  • RBF layer: The Gaussian membership functions (MFs) of every of RBF neurons in this layer is selected to deal with the input variables. Each RBF neuron represents an if-part of a fuzzy rule, and the outputs of RBF neurons are calculated in the following manner:
    ψ j = exp ( i = 1 N ( ( u i c i j ) 2 / 2 σ i j 2 ) ) , i = 1 , 2 , ... , N ; j = 1 , 2 , ... , M
    where ψ j is the output of the jth RBF neuron; c i j and σ i j are the center and width of the ith membership function (MF) in the jth neuron, respectively; and M is the total number of neurons in this layer.
  • Normalized layer: The number of the neurons in the normalized layer is the same as that in the RBF layer. The output values of nodes in this layer are given as follows:
    v l = exp ( i = 1 N ( ( u i c i l ) 2 / 2 σ i l 2 ) ) j = 1 M e x p ( i = 1 N ( ( u i c i j ) 2 / 2 σ i j 2 ) ) , i = 1 , 2 , ... , N ; j = 1 , 2 , ... , M ; l = 1 , 2 , ... , M
    where v l is the lth output value in the normalized layer.
  • Output layer: There is only one neuron in this layer, in which of the output represents the PM2.5 concentration that can be clarified through the gravity method given as follows:
    p = l = 1 M w l exp ( i = 1 N ( ( u i c i l ) 2 / 2 σ i l 2 ) ) j = 1 M exp ( i = 1 N ( ( u i c i j ) 2 / 2 σ i j 2 ) ) , i = 1 , 2 , ... , N ; j = 1 , 2 , ... , M ; l = 1 , 2 , ... , M
    where w l is the weight connecting the lth neuron in the normalized layer and the neuron in the output layer.
Moreover, 130 h of measurements of dominating variables and PM2.5 from the total 182 rows of data were randomly selected as training set of the SOG-SASOFNN, and the remaining data were regarded as test set. Thus, both the training set and the test set cover measurements from haze period (from 16 to 19 January 2010) and non-haze days (all other days from 14 to 23 January 2010). All data sets were normalized to the range of [0, 1] by linear scaling. Once the structure and the parameters of the SOG-SASOFNN were optimized through training, this optimized function was used to make prediction on the test data. The statistical parameters index of agreement (IA), coefficient of determination (R2), normalized mean bias (NMB), normalized mean gross error (NMGE), root mean square error (RMSE) and mean bias (MB) were used to assess the models performances between observed and predicted concentrations of PM2.5 [39,40]. The definitions of these statistical parameters are shown in Table 2, where p i refers to the ith predicted value and o i to the ith observed one for a total of n observations; and p ¯ and o ¯ are the averages of predicted values and observed values, respectively. Note that the normalized mean error (NME) used in the study of Yu et al. [13] and the NMGE described in our manuscript are equal in formula [39].

3.2.2. Sensitivity Analysis Method

In this present study, the SA was used to adjust the structure of the SOG-SASOFNN during training phase to understand the dynamic process that produces PM2.5 from pollution sources. SA is a method that evaluates the dependency of the system output on the input factors. The contribution of the output of each neuron in the normalized layer to the output of the SOG-SASOFNN was measured by global SA method in this paper. The measurement is based on the Fourier decomposition of variance [41,42]. The input parameters used for performing SA are the outputs of the normalized neurons expressed as v = [ v 1 , v 2 , ... , v M ] T . The SOG-SASOFNN output representing the PM2.5 concentration can be described as follows:
p = f ( v 1 , v 2 , ... , v M )
The first order sensitivity index of the output of the hth normalized neuron to the SOG-SASOFNN output is given by:
S h = V a r v h [ E ( p | v h = β h ) ] V a r ( p ) , h = 1 , 2 , ... , M
where E ( p | v h = β h ) denotes the expected value of p under the condition that the value of v h is equal to β h and the variance V a r v h takes over all the possible values of v h . V a r ( p ) represents the variance of p.
If the range of the input factor v h is [ a h , b h ] , it can be written as:
v h ( s ) = b h + a h 2 + b h a h π arcsin ( sin ( ω h s ) ) , h = 1 , 2 , ... , M
where ω h is the fundamental frequency of v h . The equation allows each factor to oscillate in a given range periodically as the scalar variable s varying in the range [−∞, ∞]. The M-factor model in Equation (8) can be described in the frequency domain using the following relationship:
f ( s ) = f ( v 1 ( s ) , v 2 ( s ) , ... , v M ( s ) )
The expanded Fourier series of f ( s ) is given by:
f ( s ) = ω = ( A ω cos ( ω s ) + B ω sin ( ω s ) )
where the Fourier coefficients at frequency ω are defined as:
A ω = 1 2 π π π f ( s ) cos ( ω s ) d s B ω = 1 2 π π π f ( s ) sin ( ω s ) d s
where s is on the range of [ π , π ] . The variance of p is calculated by the Fourier translation:
V a r ( p ) = 2 ω = 1 ( A ω 2 + B ω 2 )
The portion of the variance of p caused by v h alone is expressed as:
V a r v h [ E ( p | v h = β h ) ] = 2 m = 1 ( A m ω h 2 + B m ω h 2 ) , h = 1 , 2 , ... , M
where A m ω h and B m ω h denote the Fourier coefficients for the fundamental frequency of v h and its higher harmonics m ω h . Consequently, the expansion of the first order sensitivity index is given by:
S h = m = 1 ( A m ω h 2 + B m ω h 2 ) ω = 1 ( A ω 2 + B ω 2 ) , h = 1 , 2 , ... , M
Because the Fourier amplitude decreases as the frequency goes up, it is expected that the high order Fourier coefficients have negligible influence on the variation of the model output. Thus, the first order sensitivity index is approximated as:
S h m = 1 P ( A m ω h 2 + B m ω h 2 ) ω = 1 P ω max ( A ω 2 + B ω 2 ) , h = 1 , 2 , ... , M
where P is called the interference factor, which is usually set to 4 or 6 in the SA community, and ω m a x is the maximum value of the fundamental frequencies of all input factors.
In order to obtain the total sensitivity index, set ω h = 2 P max ( ω ~ h ) where max ( ω ~ h ) is the highest fundamental frequency of the remaining set of factors v ~ h (all the factors except the hth factor) to ensure that the frequencies generated by all interactions involving v h will not interfere with the frequencies induced by the nonlinear effect involving v ~ h alone. Then, the estimation of the total sensitivity index SW h is as follows:
SW h = ω = P max ( ω ~ h ) + 1 P ω h ( A ω 2 + B ω 2 ) ω = 1 P ω h ( A ω 2 + B ω 2 ) , h = 1 , 2 , ... , M
Because the outputs of neurons in the normalized layer are independent with each other, the total sensitivity index SW h can be acquired by the calculation of Fourier amplitude at fundamental frequency barely and is simplified as follows:
SW h = A ω h 2 + B ω h 2 j = 1 M ( A ω j 2 + B ω j 2 ) , h = 1 , 2 , ... , M
Based on the analysis above, the standardized total sensitivity index ST h of the output of the hth neuron in the normalized layer over the SOG-SASOFNN output is computed as follows:
ST h = SW h j = 1 M SW j , h = 1 , 2 , ... , M

3.2.3. Second Order Gradient Algorithm

To further improve the accuracy of prediction of hourly PM2.5 distribution, the parameter-learning phase of the SOG-SASOFNN was performed through the SOG algorithm proposed in the study of Xie et al. [32]. The SOG algorithm does not suffer from enormous Jacobian matrix and its side effects when training data are huge. All parameters such as centers, widths and weights of the SOG-SASOFNN adjusted by the SOG algorithm are expressed as follows:
Φ ( t ) = [ c 11 ( t ) , ... , c N M ( t ) , σ 11 ( t ) , ... , σ N M ( t ) , w 1 ( t ) , ... , w M ( t ) ]
where t is the current time (or training step). Φ ( t ) is the parameter vector of the SOG-SASOFNN at time t.
Following the Levenberg–Marquardt (LM) algorithm [31], the update rule is given by:
Φ ( t + 1 ) = Φ ( t ) + ( Q ( t ) + μ ( t ) I ) 1 g ( t )
where Q ( t ) is the quasi Hessian matrix computed as the sum of the sub matrix q k ( t ) for the kth training pattern:
Q ( t ) = k = 1 K q k ( t ) , q k ( t ) = j k T ( t ) j k ( t )
and the gradient vector g ( t ) is calculated as the sum of the sub vector η k ( t ) for the kth training pattern:
g ( t ) = k = 1 K η k ( t ) , η k ( t ) = j k T ( t ) e k ( t )
where K is the amount of rows of training pattern with value of 130. j k ( t ) is one of the rows of Jacobian matrix for pattern k described as follows:
j k ( t ) = [ e k ( t ) c 11 ( t ) , ... , e k ( t ) c N M ( t ) , e k ( t ) σ 11 ( t ) , ... , e k ( t ) σ N M ( t ) , e k ( t ) w 1 ( t ) , ... , e k ( t ) w M ( t ) ]
where e k ( t ) is the error calculated by:
e k ( t ) = o k ( t ) p k ( t )
where o k ( t ) and p k ( t ) are the observed output and predicted output when the kth training pattern is presented at time t, respectively.
Be aware that the combination coefficient μ ( t ) = g ( t ) and I is the identity matrix.

3.2.4. Design of the Second Order Sensitivity Analysis Based Self-Organizing Fuzzy Neural Network

The SOG-SASOFNN applied to forecast the hourly concentrations of PM2.5 for the next 24 h at SDZ from 14 to 23 January was trained with SA method and SOG algorithm. The computation procedure of the SOG-SASOFNN is described as follows:
  • Initialization of the SOG-SASOFNN: The initial SOG-SASOFNN is with random number of neurons in the normalized layer and the inputs of the SOG-SASOFNN are dominating variables selected through PCA method. There is one neuron in the output layer in which of the output represents the PM2.5 concentration. The parameters such as centers, widths, and weights of the SOG-SASOFNN are initially distributed on the random range of [0, 1].
  • Parameter learning: Adjust the parameters of the SOG-SASOFNN using Equation (22) with all training set for several training steps.
  • Growing phase: After some time ( Θ ) steps, calculate the standardized total sensitivity index of output of each normalized neuron to the network output using Equation (20). The hth normalized neuron is overactive and will be spilt into two new normalized neurons if ST h is larger than ε 1 . In order to guarantee the convergence, the outputs of the SOG-SASOFNN before and after the structure has been adjusted must be identical and the initial parameters of the two new normalized neurons are set as follows:
    c n e w 1 = c n e w 2 = c h ( t ) σ n e w 1 = σ n e w 2 = σ h ( t ) w n e w 1 = τ w h ( t ) , w n e w 2 = ( 1 τ ) w h ( t )
    where new1 and new2 denote the two new normalized neurons. c n e w 1 , σ n e w 1 and w n e w 1 are the center vector, width vector and weight of neuron new1, respectively. c n e w 2 , σ n e w 2 and w n e w 2 are the center vector, width vector and weight of neuron new2, respectively. c h ( t ) , σ h ( t ) and w h ( t ) are the center vector, width vector and weight of the hth normalized neuron before the structure has been adjusted at step t, respectively. τ is a random number which is distributed in the range of [0, 1].
  • Pruning phase: The hth normalized neuron is useless and will be pruned if ST h is less than ε 2 . To reduce the fluctuation of the output of the network, the parameters of the nearest neuron are compensated as:
    c n e a = c n e a ( t ) σ n e a = σ n e a ( t ) w n e a = w n e a ( t ) + w h ( t ) v h ( t ) / v n e a ( t )
    where nea is the normalized neuron with the minimum Euclidean distance to the hth normalized neuron and ST n e a ε 2 . c n e a , σ n e a and w n e a are the center vector, width vector and weight of neuron nea after pruning at step t, respectively. c n e a ( t ) , σ n e a ( t ) and w n e a ( t ) are the center vector, width vector and weight of neuron nea before pruning at step t, respectively. w h ( t ) is the weight of the hth normalized neuron before pruning when training to step t. v h ( t ) and v n e a ( t ) are the output of the hth normalized neuron and neuron nea before pruning when training to step t, respectively.
  • Relearning of parameters: Turn the algorithm to procedure 2 to make the parameters under relearning applying Equation (22). The training process terminates when the process achieves the expected training RMSE Ed or reaches the pre-set running step Rmax.
  • Test stage: Once the SOG-SASOFNN is optimized by the training set, this optimized nonlinear function is used to make prediction on the test data.

4. Results and Discussion

In this paper, the dominating variables most correlated to PM2.5 were firstly selected by the PCA method using the data measured at SDZ from 14 to 23 January 2010. Then, the SOG-SASOFNN was developed to forecast the PM2.5 concentrations of the next 24 h with the current hour values of dominating variables at SDZ during the studied period. Finally, the prediction performance of the SOG-SASOFNN was compared with that of the FOG-SASOFNN and the ESN as well as with the performance of the Eta-CMAQ in modeling the hourly PM2.5 concentrations at the rural sites over the eastern United States from 15 July to 19 August 2004 to verify the effectiveness of using SOG-SASOFNN to estimate the hourly PM2.5 distribution at SDZ.

4.1. Variation of PM2.5 Concentrations with Meteorological Conditions and Aerosol Optical Depth at SDZ

A regional haze episode occurred in the BTH region in the NCP from 16 to 19 January 2010 and was caused by a surface high-pressure system from 16 to 18 January 2010 and a low-pressure system on 20 January, which was unfavorable for the dispersion of pollutants. Figure 3 shows the measured temporal variation of PM2.5 and the modeled results for the SOG-SASOFNN, the FOG-SASOFNN and the ESN. According to Figure 3, the PM2.5 concentrations at SDZ during the haze period are much higher than that in non-haze days although there were no significant pollution sources within 30 km of the site, mainly owing to the regional transportation from the southern urban area. The surface wind at SDZ turning from east-northeast to southwest in the afternoon of 18 January, which carried pollutants to this area, led to the highest aerosol loading on 19 January. During the haze period, the T increased continuously day by day and the RH remained high at SDZ. PM2.5 concentrations in the haze episode were also high enough to extinct the light and the low Vis recorded in daytime indicated the haze phenomenon. Similar to the variation of PM2.5 concentration, the AOD began to increase on 16 January and increased significantly in the afternoon of 18 January to reach its maxima on 19 January. Under the Mongolia anticyclone with strong northerly wind, the haze episode was finally terminated on 20 January. Evidently, the variation of the PM2.5 concentrations at SDZ station relies seriously on the meteorological variables and is closely related to the AOD. Besides, we could find that the modeled result of the SOG-SASOFNN is more reasonable than that of the FOG-SASOFNN and the ESN.

4.2. Dominating Variables Selected by Principal Component Analysis

PCA was performed to extract the dominating factors for PM2.5 among the predictors at SDZ station. The percentages of process variance explained by principal components (PCs) defined as PCs = [ PCs 1 , PCs 2 , ... , PCs 11 ] are shown in Figure 4 and Table 3.
As shown in Table 3, over 90% of the variation within the data can be explained by the former five PCs, and the most important variables (dominating variables) for the former five PCs are RH, Pre, AOD, WS and WD according to the coefficients of the selected PCs given in Equation (29):
PCs 1 = [ 0.06 0.07 0.24 0.06 0.04 0.27 0.42 0.22 0.31 0.20 0.11 ] PCs 2 = [ 0.33 0.10 0.34 0.05 0.49 0.04 0.34 0.38 0.25 0.40 0.16 ] PCs 3 = [ 0.28 0.39 0.16 0.15 0.49 0.18 0.61 0.25 0.04 0.09 0.04 ] PCs 4 = [ 0.12 0.13 0.80 0.36 0.42 0.07 0.08 0.01 0.06 0.03 0.01 ] PCs 5 = [ 0.20 0.38 0.27 0.77 0.15 0.25 0.11 0.22 0.06 0.03 0.01 ]
It could be noted that the dominating variables (RH, Pre, AOD, WS and WD) extracted through PCA are consistent with the characteristics of PM2.5 at SDZ where the PM2.5 concentrations are severely correlated to the meteorological parameters and the AOD.

4.3. Modeling: Training and Validation

The SOG-SASOFNN was developed to predict the hourly PM2.5 concentrations 24 h ahead at SDZ along with the FOG-SASOFNN and the ESN in the present study. For SOG-SASOFNN, the first 130 h of 182 rows of randomized and normalized measurements of dominating variables and PM2.5 were used as training set and the rest 52 observations were taken as test set. The input layer of the SOG-SASOFNN consisted of five parameters containing RH, Pre, AOD, WS and WD selected by PCA. There were three neurons in the normalized layer before the input data have been loaded. The parameters (centers, widths and weights) of the SOG-SASOFNN were initialized on the random range of [0, 1]. The structure of the SOG-SASOFNN was adjusted based on the contribution of the output of each normalized neuron to the output of network measured by global SA method. The parameter-learning phase was performed using Equation (22). The network was trained for 100 steps with Θ = 10 unless the training RMSE was less than 0.01. The input parameters of the FOG-SASOFNN and the ESN were the same as that of the SOG-SASOFNN. The FOG-SASOFNN was optimized after it was trained for 500 steps or the training RMSE was less than 0.01. The weight-learning process for the ESN was implemented applying linear regression equation [11]. The description and the evaluation of the PM2.5 forecasts over the eastern United States using the Eta-CMAQ model was presented in detail in the study of Yu et al. [13].
The change of the training RMSE of the SOG-SASOFNN and the FOG-SASOFNN during learning phase is shown in Figure 5. The growing and pruning process of the neurons in the normalized layer of the SOG-SASOFNN and the FOG-SASOFNN for the training period has been shown in Figure 6. In Figure 5 and Figure 6, we could find that the SOG-SASOFNN converges with lower training RMSE and more compact structure under much less training steps compared to the FOG-SASOFNN.
Figure 7 shows the scatter plots between the observed and predicted concentrations of PM2.5 during training phase and test phase for the SOG-SASOFNN, the FOG-SASOFNN and the ESN. As can be seen in Figure 7, the SOG-SASOFNN possesses stronger learning ability and better generalization performance compared to the FOG-SASOFNN and the ESN.
The models’ performances were evaluated by the calculation of several statistical parameters indexes described in Table 2. The statistical performances for the SOG-SASOFNN, the FOG-SASOFNN, the ESN and the Eta-CMAQ are given in Table 4. It can be found that the values of R2, NMB, NMGE and MB for the SOG-SASOFNN in both training phase and test phase are closer to their corresponding ideal values contrasted to that for the FOG-SASOFNN, the ESN as well as the Eta-CMAQ. The training RMSE and test RMSE of the SOG-SASOFNN are larger than the Eta-CMAQ RMSE, which may be due to the data used in our paper covers wider PM2.5 concentrations scope compared with the data for the Eta-CMAQ. Through the analysis above, we can conclude that the SOG-SASOFNN performs better than the FOG-SASOFNN, the ESN and the Eta-CMAQ in predicting the hourly concentrations of PM2.5. Overall, the performance of the proposed SOG-SASOFNN model is satisfactory.

5. Conclusions

Estimating the distribution of PM2.5 has proven to be tough due to the highly chaotic and nonlinear phenomena existing in the atmospheric and chemical processes that result in air pollution. In order to understand the complex time series of PM2.5 precisely, it is necessary to use the SOFNN to analyze the underlying dynamic process that produces PM2.5. However, as far as we know, few studies have applied the SOFNN to forecast the PM2.5 concentrations. In this paper, the SOG-SASOFNN with SA for structure-learning phase and SOG algorithm for parameter-learning phase was developed to model the hourly concentrations of PM2.5 of the next 24 h from 14 to 23 January 2010 at SDZ. The input parameters for the SOG-SASOFNN were selected through PCA method to eliminate the irrelevant variables. It is worth mentioning that the dominating variables (RH, Pre, AOD, WS and WD) extracted by PCA were consistent with the characteristics of PM2.5 at SDZ where there were no significant pollution sources within 30 km of the site and the PM2.5 concentrations showed strong correlation with the meteorological parameters and the AOD. The prediction results showed that the SOG-SASOFNN performed better than the FOG-SASOFNN, the ESN and the Eta-CMAQ in estimating the hourly distribution of PM2.5. Obviously, the developed SOG-SASOFNN model gives satisfying results for prediction of hourly concentrations of PM2.5 at SDZ during the studied period. Moreover, the SOG-SASOFNN performance could be further improved by considering the effect of the local anthropogenic emission activities on PM2.5 concentrations and utilizing more sufficient measurements for simulation, which can be guidance for our further work.
Ultimately, it is expected that the SOG-SASOFNN proposed here can be applied to interpret the complex time series of any pollutant affected by emission sources, meteorological conditions and topographical characteristics of the area under study to help in forecasting of the air pollution over other regions of the world.

Acknowledgments

The authors would like to acknowledge the support from the National Natural Science Foundation of China (61533002 and 61622301), the Beijing Municipal Education Commission Foundation (km201410005001 and KZ201410005002) and the China Postdoctoral Science Foundation (2014M550017 and 2015M570911). We also want to acknowledge the editing and reviewing for this paper provided by Wenjing Li from Beijing University of Technology. The authors wish to thank the anonymous reviewers for their constructive remarks that helped improve the manuscript.

Author Contributions

Junfei Qiao and Honggui Han conceived and supervised the topic of this paper. Honggui Han and Jie Cai proposed the methods and the data. Jie Cai and Jianxian Cai analyzed the results and wrote the paper. All authors read and approved the submitted manuscript and accepted the version for publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Song, Y.; Zhang, Y.H.; Xie, S.D.; Zeng, L.M.; Zheng, M.; Salmon, L.G.; Shao, M.; Slanina, S. Source apportionment of PM2.5 in Beijing by positive matrix factorization. Atmos. Environ. 2006, 40, 1526–1537. [Google Scholar] [CrossRef]
  2. Lv, B.L.; Hu, Y.T.; Chang, H.H.; Russell, A.G.; Bai, Y.Q. Improving the accuracy of daily PM2.5 distributions derived from the fusion of ground-level measurements with aerosol optical depth observations, a case study in north China. Environ. Sci. Technol. 2016, 50, 4752–4759. [Google Scholar] [CrossRef] [PubMed]
  3. Pope, C.A., III; Dockery, D.W. Health effects of fine particulate air pollution: Lines that connect. J. Air Waste Manag. 2006, 56, 709–742. [Google Scholar] [CrossRef]
  4. Lepeule, J.; Laden, F.; Dockery, D.; Schwartz, J. Chronic exposure to fine particles and mortality: An extended follow-up of the Harvard Six Cities study from 1974 to 2009. Environ. Health Perspect. 2012, 120, 965–970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Metzger, K.B.; Tolbert, P.E.; Klein, M.; Peel, J.L.; Flanders, W.D.; Todd, K.; Mulholland, J.A.; Ryan, P.B.; Frumkin, H. Ambient air pollution and cardiovascular emergency department visits. Epidemiology 2004, 15, 46–56. [Google Scholar] [CrossRef] [PubMed]
  6. Peel, J.L.; Tolbert, P.E.; Klein, M.; Metzger, K.B.; Flanders, W.D.; Todd, K.; Mulholland, J.A.; Ryan, P.B.; Frumkin, H. Ambient air pollution and respiratory emergency department visits. Epidemiology 2005, 16, 164–174. [Google Scholar] [CrossRef] [PubMed]
  7. Levy, H., II; Horowitz, L.W.; Schwarzkopf, M.D.; Ming, Y.; Golaz, J.C.; Naik, V.; Ramaswarmy, V. The roles of aerosol direct and indirect effects in past and future climate change. J. Geophys. Res. 2013, 118, 4521–4532. [Google Scholar]
  8. Zhao, P.S.; Zhang, X.L.; Xu, X.F.; Zhao, X.J. Long-term visibility trends and characteristics in the region of Beijing, Tianjin, and Hebei, China. Atmos. Res. 2011, 101, 711–718. [Google Scholar] [CrossRef]
  9. Valverde, V.; Pay, M.T.; Baldasano, J.M. Circulation-type classification derived on a climatic basis to study air quality dynamics over the Iberian Peninsula. Int. J. Climatol. 2015, 35, 2877–2897. [Google Scholar] [CrossRef] [Green Version]
  10. Raga, G.B.; Moyne, L.L. On the nature of air pollution dynamics in Mexico City—I. Nonlinear analysis. Atmos. Environ. 1996, 30, 3987–3993. [Google Scholar] [CrossRef]
  11. Xu, Z.; Xia, X.P.; Liu, X.N.; Qian, Z.G. Combining DMSP/OLS nighttime light with echo state network for prediction of daily PM2.5 average concentrations in Shanghai, China. Atmosphere 2015, 6, 1507–1520. [Google Scholar] [CrossRef]
  12. Chemel, C.; Fisher, B.E.A.; Kong, X.; Francis, X.V.; Sokhi, R.S.; Good, N.; Collins, W.J.; Folberth, G.A. Application of chemical transport model CMAQ to policy decisions regarding PM2.5 in the UK. Atmos. Environ. 2014, 82, 410–417. [Google Scholar] [CrossRef] [Green Version]
  13. Yu, S.C.; Mathur, R.; Schere, K.; Kang, D.W.; Pleim, J.; Young, J.; Tong, D.; Pouliot, G.; McKeen, S.A.; Rao, S.T. Evaluation of real-time PM2.5 forecasts and process analysis for PM2.5 formation over the eastern United States using the Eta-CMAQ forecast model during the 2004 ICARTT study. J. Geophys. Res. 2008, 113, D06204. [Google Scholar] [CrossRef]
  14. Fernando, H.J.S.; Mammarella, M.C.; Grandoni, C.; Fedele, P.; Marco, R.D.; Dimitrova, R.; Hyde, P. Forecasting PM10 in metropolitan areas: Efficacy of neural networks. Environ. Pollut. 2012, 163, 62–67. [Google Scholar] [CrossRef] [PubMed]
  15. Elbayoumi, M.; Ramli, N.A.; Yusof, N.F.F.M.; Yahaya, A.S.B.; Madhoun, W.A.; UI-Saufie, A.Z. Multivariate methods for indoor PM10 and PM2.5 modeling in naturally ventilated schools buildings. Atmos. Environ. 2014, 94, 11–21. [Google Scholar] [CrossRef]
  16. Niska, H.; Hiltunen, T.; Karppinen, A.; Ruuskanen, J.; Kolehmainen, M. Evolving the neural network model for forecasting air pollution time series. Eng. Appl. Artif. Intell. 2004, 17, 159–167. [Google Scholar] [CrossRef]
  17. Nagendra, S.M.S.; Khare, M. Artificial neural network approach for modelling nitrogen dioxide dispersion from vehicular exhaust emissions. Ecol. Model. 2006, 190, 99–115. [Google Scholar] [CrossRef]
  18. Corani, G. Air quality prediction in Milan: Feed-forward neural networks, pruned neural networks and lazy learning. Ecol. Model. 2005, 185, 513–529. [Google Scholar] [CrossRef]
  19. Chaloulakou, A.; Grivas, G.; Spyrellis, N. Neural network and multiple regression models for PM10 prediction in Athens: A comparative assessment. J. Air Waste Manage. 2003, 53, 1183–1190. [Google Scholar] [CrossRef]
  20. Ordieres, J.B.; Vergara, E.P.; Capuz, R.S.; Salazar, R.E. Neural network prediction model for fine particulate matter (PM2.5) on the US-Mexico border in El Paso (Texas) and Ciudad Juaŕez (Chihuahua). Environ. Model. Softw. 2005, 20, 547–559. [Google Scholar] [CrossRef]
  21. Lin, C.J. An efficient immune-based symbiotic particle swarm optimization learning algorithm for TSK-type neuro-fuzzy networks design. Fuzzy Sets Syst. 2008, 159, 2890–2909. [Google Scholar] [CrossRef]
  22. Assimakopoulos, M.N.; Dounis, A.; Spanou, A.; Santamouris, M. Indoor air quality in a metropolitan area metro using fuzzy logic assessment system. Sci. Total Environ. 2013, 449, 461–469. [Google Scholar] [CrossRef] [PubMed]
  23. Mishra, D.; Goyal, P. Neuro-fuzzy approach to forecast NO2 pollutants addressed to air quality dispersion model over Delhi, India. Aerosol Air Qual. Res. 2016, 16, 166–174. [Google Scholar] [CrossRef]
  24. Mishra, D.; Goyal, P.; Upadhyay, A. Artificial intelligence based approach to forecast PM2.5 during haze episodes: A case study of Delhi, India. Atmos. Environ. 2015, 102, 239–248. [Google Scholar] [CrossRef]
  25. Wu, S.Q.; Er, M.J.; Gao, Y. A fast approach for automatic generation of fuzzy rules by generalized dynamic fuzzy neural networks. IEEE Trans. Fuzzy Syst. 2001, 9, 578–594. [Google Scholar] [CrossRef]
  26. Wu, S.Q.; Er, M.J. Dynamic fuzzy neural networks-a novel approach to function approximation. IEEE Trans. Syst. Man. Cybern. B. Cybern. 2000, 30, 358–364. [Google Scholar] [PubMed]
  27. Kupongsak, S.; Tan, J. Application of fuzzy set and neural network techniques in determining food process control set points. Fuzzy Sets Syst. 2006, 157, 1169–1178. [Google Scholar] [CrossRef]
  28. Leng, G.; McGinnity, T.M.; Prasad, G. Design for self-organizing fuzzy neural networks based on genetic algorithms. IEEE Trans. Fuzzy Syst. 2006, 14, 755–766. [Google Scholar] [CrossRef]
  29. Habbi, H.; Boudouaoui, Y.; Karaboga, D.; Ozturk, C. Self-generated fuzzy systems design using artificial bee colony optimization. Inf. Sci. 2015, 295, 145–159. [Google Scholar] [CrossRef]
  30. Han, H.G.; Qiao, J.F. A self-organizing fuzzy neural network based on a growing-and-pruning algorithm. IEEE Trans. Fuzzy Syst. 2010, 18, 1129–1143. [Google Scholar] [CrossRef]
  31. Wilamowski, B.M.; Yu, H. Improved computation for Levenberg-Marquardt training. IEEE Trans. Neural Netw. 2010, 21, 930–937. [Google Scholar] [CrossRef] [PubMed]
  32. Xie, T.T.; Yu, H.; Hewlett, J.; Rozycki, P.; Wilamowski, B. Fast and efficient second-order method for training radial basis function networks. IEEE Trans. Neural Netw. Learn. Syst. 2012, 23, 609–619. [Google Scholar] [PubMed]
  33. Voukantsis, D.; Karatzas, K.; Kukkonen, J.; Räsänen, T.; Karppinen, A.; Kolehmainen, M. Intercomparison of air quality data using principal component analysis, and forecasting of PM10 and PM2.5 concentrations using artificial neural networks, in Thessaloniki and Helsinki. Sci. Total Environ. 2011, 409, 1266–1276. [Google Scholar] [CrossRef] [PubMed]
  34. Lin, W.; Xu, X.; Zhang, X.; Tang, J. Contributions of pollutants from North China Plain to surface ozone at the Shangdianzi GAW Station. Atmos. Chem. Phys. 2008, 8, 5889–5898. [Google Scholar] [CrossRef]
  35. Zhao, X.J.; Zhang, X.L.; Xu, X.F.; Xu, J.; Meng, W.; Pu, W.W. Seasonal and diurnal variations of ambient PM2.5 concentration in urban and rural environments in Beijing. Atmos. Environ. 2009, 43, 2893–2900. [Google Scholar] [CrossRef]
  36. Zhao, X.J.; Zhao, P.S.; Xu, J.; Meng, W.; Pu, W.W.; Dong, F.; He, D.; Shi, Q.F. Analysis of a winter regional haze event and its formation mechanism in the North China Plain. Atmos. Chem. Phys. 2013, 13, 5685–5696. [Google Scholar] [CrossRef]
  37. Chavent, M.; Guegan, H.; Kuentz, V.; Patouille, B.; Saracco, J. PCA- and PMF-based methodology for air pollution sources identification and apportionment. Environmetrics 2009, 20, 928–942. [Google Scholar] [CrossRef]
  38. Viana, M.; Querol, X.; Alastuey, A.; Gil, J.I.; Menendez, M. Identification of PM sources by principal component analysis (PCA) coupled with wind direction data. Chemosphere 2006, 65, 2411–2418. [Google Scholar] [CrossRef] [PubMed]
  39. Carslaw, D.C.; Ropkins, K. Openair-an R package for air quality data analysis. Environ. Model. Softw. 2012, 27, 52–61. [Google Scholar] [CrossRef]
  40. Elangasinghe, M.A.; Singhal, N.; Dirks, K.N.; Salmond, J.A.; Samarasinghe, S. Complex time series analysis of PM10 and PM2.5 for a coastal site using artificial neural network modeling and k-means clustering. Atmos. Environ. 2014, 94, 106–116. [Google Scholar] [CrossRef]
  41. Lauret, P.; Fock, E.; Mara, T.A. A node pruning algorithm based on a Fourier amplitude sensitivity test method. IEEE Trans. Neural Netw. 2006, 17, 273–293. [Google Scholar] [CrossRef] [PubMed]
  42. Han, H.G.; Qiao, J.F. A structure optimisation algorithm for feedforward neural network construction. Neurocomput 2013, 99, 347–357. [Google Scholar] [CrossRef]
Figure 1. Location of Shangdianzi (SDZ) site.
Figure 1. Location of Shangdianzi (SDZ) site.
Atmosphere 08 00010 g001
Figure 2. Initial architecture of sensitivity analysis based self-organizing fuzzy neural network with second order gradient algorithm (SOG-SASOFNN).
Figure 2. Initial architecture of sensitivity analysis based self-organizing fuzzy neural network with second order gradient algorithm (SOG-SASOFNN).
Atmosphere 08 00010 g002
Figure 3. Time series variation of measured PM2.5 concentrations and the modeled results of SOG-SASOFNN, SASOFNN with first order gradient algorithm (FOG-SASOFNN) and echo state network (ESN) at SDZ from 14 to 23 January 2010.
Figure 3. Time series variation of measured PM2.5 concentrations and the modeled results of SOG-SASOFNN, SASOFNN with first order gradient algorithm (FOG-SASOFNN) and echo state network (ESN) at SDZ from 14 to 23 January 2010.
Atmosphere 08 00010 g003
Figure 4. Percentages of process variance explained by principal components (PCs).
Figure 4. Percentages of process variance explained by principal components (PCs).
Atmosphere 08 00010 g004
Figure 5. Training root mean square error (RMSE) values of SOG-SASOFNN and FOG-SASOFNN.
Figure 5. Training root mean square error (RMSE) values of SOG-SASOFNN and FOG-SASOFNN.
Atmosphere 08 00010 g005
Figure 6. Dynamic structure process of SOG-SASOFNN and FOG-SASOFNN.
Figure 6. Dynamic structure process of SOG-SASOFNN and FOG-SASOFNN.
Atmosphere 08 00010 g006
Figure 7. Scatter plots between the observed and the predicted PM2.5 concentrations for models: (a) SOG-SASOFNN for training set; (b) SOG-SASOFNN for test set; (c) FOG-SASOFNN for training set; (d) FOG-SASOFNN for test set; (e) ESN for training set; and (f) ESN for test set.
Figure 7. Scatter plots between the observed and the predicted PM2.5 concentrations for models: (a) SOG-SASOFNN for training set; (b) SOG-SASOFNN for test set; (c) FOG-SASOFNN for training set; (d) FOG-SASOFNN for test set; (e) ESN for training set; and (f) ESN for test set.
Atmosphere 08 00010 g007aAtmosphere 08 00010 g007b
Table 1. Units of variables measured at SDZ.
Table 1. Units of variables measured at SDZ.
VariablesTRHWSWDPreVisCONO2O3SO2PM2.5
Units°C%m/s°hPakmppmppbppbppbμg/m3
Table 2. Definitions of statistical parameters used for assessing the models performances.
Table 2. Definitions of statistical parameters used for assessing the models performances.
Statistical ParameterDescriptionMathematical Function
IAExpresses the difference between predicted and observed values I A = 1 i = 1 n | p i o i | 2 i = 1 n ( | p i o ¯ | + | o i o ¯ | ) 2
R2A measure of linear relationship between predicted and observed values R 2 = ( i = 1 n ( p i p ¯ ) ( o i o ¯ ) ) 2 i = 1 n ( p i p ¯ ) 2 i = 1 n ( o i o ¯ ) 2
NMBIndicates over or under estimation of the model N M B = i = 1 n ( p i o i ) i = 1 n o i
NMGEIndicates mean error regardless of it is over or under estimation N M G E = i = 1 n | p i o i | i = 1 n o i
RMSEProvides an overall measure of how close predicted values and observed values are R M S E = i = 1 n ( p i o i ) 2 n
MBMeasure of model bias M B = i = 1 n ( p i o i ) n
Table 3. Percent and additive percent of process variance caused by PCs.
Table 3. Percent and additive percent of process variance caused by PCs.
PCs1234567891011
Percent (%)56.5515.7110.005.984.012.341.711.501.130.860.21
Additive percent (%)56.5572.2682.2688.2492.2594.5996.3097.8098.9399.79100
Table 4. Statistical performances of SOG-SASOFNN, FOG-SASOFNN, ESN and Eta-Community Multiscale Air Quality (Eta-CMAQ) model.
Table 4. Statistical performances of SOG-SASOFNN, FOG-SASOFNN, ESN and Eta-Community Multiscale Air Quality (Eta-CMAQ) model.
Statistical ParameterIdeal ValueSOG-SASOFNNFOG-SASOFNNESNEta-CMAQ
TrainingTestTrainingTestTrainingTest
IA10.970.950.910.860.800.70~
R210.890.840.720.700.510.300.22 *
NMB0−0.01−0.05−0.17−0.23−0.260.31−0.32 *
NMGE00.250.370.420.430.520.570.51*
RMSE (μg/m3)013.5617.9026,8429.3135.3539.7911.6 *
MB (μg/m3)0−0.16−1.86−2.01−3.6−5.616.33−5.2 *
~ The value does not shown in original paper; * The results are the same as the original paper.

Share and Cite

MDPI and ACS Style

Qiao, J.; Cai, J.; Han, H.; Cai, J. Predicting PM2.5 Concentrations at a Regional Background Station Using Second Order Self-Organizing Fuzzy Neural Network. Atmosphere 2017, 8, 10. https://doi.org/10.3390/atmos8010010

AMA Style

Qiao J, Cai J, Han H, Cai J. Predicting PM2.5 Concentrations at a Regional Background Station Using Second Order Self-Organizing Fuzzy Neural Network. Atmosphere. 2017; 8(1):10. https://doi.org/10.3390/atmos8010010

Chicago/Turabian Style

Qiao, Junfei, Jie Cai, Honggui Han, and Jianxian Cai. 2017. "Predicting PM2.5 Concentrations at a Regional Background Station Using Second Order Self-Organizing Fuzzy Neural Network" Atmosphere 8, no. 1: 10. https://doi.org/10.3390/atmos8010010

APA Style

Qiao, J., Cai, J., Han, H., & Cai, J. (2017). Predicting PM2.5 Concentrations at a Regional Background Station Using Second Order Self-Organizing Fuzzy Neural Network. Atmosphere, 8(1), 10. https://doi.org/10.3390/atmos8010010

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