Abstract

In this work, the air pollution index in three cities (Seoul, Busan, and Daegu) in South Korea was studied using multifractal detrended fluctuation analysis (MF-DFA). Hurst, Renyi, and Holder exponents were used to analyze the characteristics of the concentration time series of and . The results showed that multifractality exists in each season interval and the multifractal degree of is stronger than that of . To investigate the effects of the implementation of the “haze special law” on February 15, 2019, we analyzed the time series of during the time periods from February 15, 2018, to December 16, 2018, and February 15, 2019, to December 16, 2019. We found that the multifractal spectrum width after the implementation of the law was narrower than that before the law for all the cities, which shows that the enactment of the law has played a role in improving the efficiency of air pollution control in South Korea. We also conclude that the major effects of the law will be particularly visible in larger cities. To study the main causes of multifractality, the shuffled and phase-randomized series were analyzed using MF-DFA, and the results demonstrated that the fat-tailed distribution resulted in the multifractality of the time series before and after the implementation of the “haze special law” in Seoul and Daegu, whereas long-range correlation resulted in multifractality of the series before and after the implementation of the law in Busan.

1. Introduction

With the rapid development of a social economy and the acceleration of industrialization and urbanization, air pollution has gradually become an important problem affecting people’s lives and has therefore garnered considerable attention from the public and government. Air pollution is closely related to climate, ecology, and health. Studies have shown that sulfur oxides and nitrogen oxides can cause great harm to the respiratory system, even leading to respiratory failure in severe cases [1]. Therefore, many researchers have started to investigate the atmospheric particulate matters. Many studies have stated that and are the two major pollutants among various air pollutants [26]. Therefore, in this research, we focused on and and investigated the two pollutants based on the relationship between the air pollution index (API) and different seasons. We also checked whether the implementation of the “haze special law” on February 15, 2019, had an impact on improving the air quality in South Korea using multifractal detrended fluctuation analysis (MF-DFA).

We know that the MF-DFA model can unveil the multifractal properties hidden in nonstationary time series and that multifractality originates from temporal correlations [7, 8]. Therefore, the multifractal properties of stock markets [913], foreign exchange markets [1416], bitcoin markets [17], atmospheric sciences [18], and phase transitions [19] have been studied previously. Recently, multifractal analysis was used to study the daily air temperature time series [20]. Moreover, MF-DFA is also an efficient method in analyzing the human heart rate time series [21]. Based on the analysis of blood pressure and heart-rate complexity, Gender Castiglioni et al. [22] used the multifractal technique to investigate cardiovascular warning signals. Moreover, MF-DFA has been used in some studies to implement specific image analysis [23, 24].

In recent years, many studies on air pollutants have been conducted and the results suggested that the increased levels of are associated with a higher mortality and some negative effects on the lungs [25]. Although there are many studies on API, few studies have focused on the MF-DFA of the API time series. Through the hourly average concentration time series, Shi et al. [26] analyzed the multifractal nature at four air monitoring stations of Chengdu using MF-DFA. Zhang et al. [27] used the multifractal detrended cross-correlation analysis (MF-DCCA) to analyze the cross-correlations between and meteorological factors. Recently, Zhang et al. [28] analyzed the multifractal characteristics of the time series in Hong Kong using the empirical mode decomposition-based MF-DFA method. Multifractal property links between meteorological factors and pollutants in urban and rural areas have been verified by He [29].

Based on MF-DFA, the variations of multifractal characteristics of pollutants present in different seasons were investigated in this study. The effects of seasonal factors on and were analyzed, and the seasonal physical changes of and concentrations were confirmed using multifractality. Additionally, we conducted a multifractal comparison analysis of the effectiveness of air pollution renovation since the enactment of the “haze special law.” From this, we can provide a basis for formulating a scientific and effective comprehensive air pollution control policy.

The paper is organized as follows. We briefly outline the procedure of MF-DFA in Section 2. In the next section, we describe the data information. Section 4 illustrates the empirical results. Section 5 concludes the paper.

2. Methodology

In this section, we describe the utilization of the MF-DFA [30] to measure the multifractal behavior of haze in all seasons, and the multifractal characteristics of the API time series before and after the implementation of the “haze special law.” Kantelhardt summarized the technical details of MF-DFA as follows.

As a time series starts from to , where N is the length of the signal, the corresponding summation sequence is constructed by the following integration:where is the mean value of .

Subsequently, the profile Y is further divided into nonoverlapping windows of equal length s. In most cases, the scale s is not a multiple of the time series and a short part at the end of series Y exists. To overcome the problem of information being lost in the division process, the same process is repeated starting from the other side of the series. Thus, windows are acquired.

Next, the least squares method is used to fit the data for evaluating the local trend of each window , and the fitting polynomial in the window is denoted by . The variance is determined by the detrended time series, which is calculated as the difference between Y and y, and the resultant equation is described asif , andif .

Finally, computing the mean of windows, the q order wave function is obtained asif ; according to L’Hôpital’s rule,

Through the analysis of double log plots of versus s and a varied q, the scaling behavior of the fluctuation is determined by the power-law , and from this, a family of scaling exponents , which are generalized Hurst components, can be obtained. is the standard DFA, if . The Hurst exponents provide information on the time series, such as power-law correlated behavior, and when , it indicates that the time series has a negative or antipersistence property. When , then the time series has a positive persistence, and indicates that the time series has an uncorrelated Brownian process.

The Renyi exponent , which is related to the general Hurst exponent, can be expressed byIn addition,where α represents the Holder exponent and characterizes the singularity strength and is a fractal dimension of the set of points with particular α. In the plotted curve between α and , the shape resembles an inverted parabola and the degree of their complexity is denoted by the width of their fractal strength [31].

3. Data Collection

We used the API time series of Seoul to study the multifractal characteristics of API sequences in different seasons. To gauge the impact of the implementation of the “haze special law,” our data set covered 3 cities in South Korea. The 3 cities are located in different parts of South Korea, and there is a certain distance between any two cities. The results calculated from the samples selected in this manner fully reflect the impact and effect of the implementation of the “haze special law” over the entirety of South Korea.

The experimental data sources for the and concentrations were the “Korean Ministry Of Environment” and “Korea Environment Corporation,” and more specific data were provided by the monitoring stations in South Korea. For more detailed information, please refer to https://www.airkorea.or.kr/web and https://aqicn.org/city/seoul/kr. Table 1 shows information of the target cities and the corresponding monitoring stations. For the multifractal analysis of API in different seasons, the samples were recorded from March 3, 2018, to February 25, 2019. We considered March to May, June to August, September to November, and December to February as the four seasons: spring, summer, autumn, and winter, respectively. Each interval was assigned a time series of 90 daily data points. For analyzing the efficiency of the implementation of the “haze special law” on February 15, 2019, we selected a time series of over the same period. The data were recorded from February 15, 2018, to December 16, 2018, and from February 15, 2019, to December 16, 2019, and each interval was assigned a time series of 305 daily data points. The descriptive statistics for each time series are shown in Table 2.

4. Experiment Results

We first analyze the multifractal properties for two time series of each season and compare the multifractal degree of with that of . Second, we conduct multifractal analysis for two time series in the same season and confirm whether the implementation of the “haze special law” has any implications for the air quality.

In this study, we choose the minimum segment scale to be , and the maximum segment scale , and q is taken from to 10. The range of q was similarly chosen for the multifractal analysis of stock prices and returns in Latin-American stock markets [32]. Besides, the range of s was selected according to suggestions in [33, 34].

4.1. Multifractal Analysis of and

From the daily API data selected in Section 3, we draw the time series of and in Figure 1. From Figure 1, we can see that a pattern of high concentrations on both sides and low concentrations in the middle emerges, which indicates that due to seasonal factors, both air pollutants have the characteristics of low concentrations in summer and high concentrations in winter. In summer (June to August), the concentration of air pollutants tends to rise, while in winter (December to February of the following year), the concentration of air pollutants tends to decrease.

To check the existence of a scaling range, long-range correlations, and multifractality of both indices in four seasons, we first show the log-log plots of versus time scale s of and . In Figure 2, the continuous line represents the double plots between the calculated and time scale s, and the dashed line denotes the corresponding fitting line. As the dashed line shows, the decreasing gradient with q increases from to 10, which indicates that and time series of all seasons have multifractal properties.

To measure the multifractality quantitatively, we apply MF-DFA to estimate properties for two pollutants of indices and . First, we calculate the generalized Hurst exponent and then obtain the Renyi exponent and the Holder exponent α for both indices.

Figure 3 shows the relationship between and q. The generalized Hurst exponents for both time series in any season are not constants, which indicate that the series are not monofractal. We note that decreases monotonously with an increase in q, which shows that the time series of and are not single fractals. Figures 3(a) and 3(b) represent and , respectively.

As illustrated in Table 3, the values of for all time series are larger than 0.5, implying that the fluctuations have significant positive persistence. Both indices of the time series in some seasons are even larger than 1, which is consistent with the conclusions drawn in some previous studies [26, 35]. For the time series, the persistence in spring is stronger than for other seasons, and for the time series, autumn has the strongest persistence. We compute for each time series using . From Figure 3 and Table 3, we observe that in winter is the highest for the time series, which implies its’ multifractal characteristic is stronger in winter, and the multifractalities of the other three seasons are relatively stable. However, for the time series, the multifractality in spring is higher than for all other seasons.

Furthermore, the multifractality strength of a time series can also be measured with a nonlinear Renyi exponent curve [36]. As shown in Figure 4, all curves of the Renyi exponent are nonlinear, which reflects the evidence of multifractality of the and series. Moreover, the different nonlinear behaviors of the curves of these series show that the multifractal strengths have different degrees.

Then, we calculate the Holder exponent α and the fractal dimension using equation (8). The width of α also stands for the degree of the multifractality. A larger indicates a stronger multifractal nature. In Figure 5 and Table 3, the conclusion we described above about the Hurst exponent and the Renyi exponent is discussed. Moreover, as shown in Figure 5, the multifractal spectra are strongly asymmetrical, and as an important indicator of fractal organization, the asymmetry coefficient, which is used to quantify the asymmetry of the multifractal spectrum, has been estimated in [37].

Based on Table 3, the multifractal features can be clearly confirmed for all time series, and it can be seen that except for spring, the Hurst exponents and multifractal spectrum widths of are obviously larger than those of and the multifractal degrees of are stronger than those of .

4.2. Implications of the “Haze Special Law”

In this section, we focus on the enactment of the “haze special law” and study whether the law is effective in the governance of atmospheric pollution. We select daily data from the same period described in Section 3. As we studied above, the multifractal degree of is stronger than that of . Thus, in this section, we focus on the comparison of concentrations. We show the time series before and after the law was implemented in Seoul, Busan, and Daegu in Figure 6.

As the fluctuation statistics are of asymptotic feature, a sufficiently long-length sequence guarantees the accurate estimation of the Hurst exponent, and for a small-length sequence, the asymptotic value in the slope of versus will be lacked [38]. As the total number of observations is selected to be 305, to eliminate the doubt of whether the selected time series can give representative results for a multifractal analysis in the investigation, we dealt with this concern by computing confidence intervals of Hurst exponent to check whether there is a difference caused by the law for all the three cities. The number of samples increased from 50 to 300, and every additional interval is 10. The confidence level is selected as . For each set of sample lengths, we generate the shuffled time series 1000 times, and then, the confidence interval is plotted in Figure 7.

As shown in Figure 7 and as expected, the domain size decreases with an increase in the length of the time series. However, when the length of the time series is larger than 200, the decrease in the domain size of the confidence intervals is rather slow, and of a distinct asymptotic behavior. Besides, we see that the confidence intervals do not overlap for Seoul, while for Busan ([0.34, 0.39] before and [0.32, 0.37] after) and Daegu ([0.55, 0.64] before and [0.58, 0.69] after), the confidence intervals overlap a lot. Therefore, we can conclude that there is a difference caused by the law for Seoul; however, the same conclusion is inconclusive for Busan and Daegu.

Now, we first depict the log-log plots of the fluctuation function versus q and s for the time series. As illustrated in Figure 8, for all series, the fluctuation value increases linearly with s, indicating that a power-law behavior and long-range correlations exist in each series pair. The decreasing gradient shows that all time series have multifractal characteristics.

The generalized Hurst exponent , Renyi exponent , and Holder exponent α of the time series before and after the enactment of the “haze special law” for all the three cities are shown in Figures 911, respectively.

In Figure 9, it is shown that the generalized Hurst exponents of of all three cities decrease with varied q, thereby showing the multifractal behaviors of these two series. We present for these series in Table 4. All Hurst exponents are confirmed to be larger than 0.5, which indicates that all the time series exhibit positive persistence. Besides, we observe that before the implementation of the law is larger than after the law, implying that the fluctuations before the law have more significant persistence. Subsequently, we calculate the range of for each time series and list them in Table 4. We find that for all three cities, after the law is larger, which indicates that the multifractal characteristics are stronger, and that the multifractalities of the time series before the law are more stable.

The Renyi exponents plotted in Figure 10 are nonlinear along q, which provides further evidence of the existence of multifractality. We notice that the curvature of the Renyi exponents before the law are higher, which implies that the time series before the implementation of the law have stronger multifractal characteristics. All of these are consistent with the results calculated from the generalized Hurst exponents.

At last, the multifractal spectrums of these time series are examined. We depict the curves between α and in Figure 11. We see that the multifractal spectra are not shown as points, which indicates all the time series have multifractality. Then, are calculated and presented in Table 4. From Table 4, it can be seen that and of the time series after the enactment of law are lower than those before the implementation of the law for all three cities, implying that the multifractality of the time series becomes more stable with the measures taken by the government and that the efficiency of the “haze special law” in three cities is demonstrable; this conclusion can be referred to [16, 39], in which it was pointed out that the decrease in multifractality means the efficiency is improved when some measures are taken. Moreover, we also notice that the efficiency of the “haze special law” in Seoul is most obvious, when compared with Busan, and in last place is Daegu. Interestingly, the order is the same as the scale and prosperity of the cities. We can therefore conclude that the law has a higher effect on a larger city.

4.3. Sources of Multifractal Features

The two universally acknowledged sources of multifractality are long-range correlations and a fat-tailed distribution [7, 40]. Now, we investigate the major source of multifractality of the concentration series. To check the contribution of long-range correlations and fat-tailed distribution quantitatively, we shuffle and phase-randomize the time series.

Subsequently, we calculate the Hurst exponent , the Renyi exponent , and the Holder exponent of the original, shuffled, and phase-randomized time series of the three cities. As shown in Figures 1214, all the time series are strongly multifractal. From Table 5, all of the original series are larger than those of the shuffled and phase-randomized time series, which indicates that the persistence has been moved after shuffling and phase-randomization. Moreover, the original series of all three cities has the greatest spectrum width , while the multifractality decreased after shuffling and phase-randomizing the series. The findings can also be confirmed in Figures 1214.

To examine the main source of multifractality, the first and third rows of Figure 14 show that the multifractal spectra of phase-randomized series are the narrowest ones, which means that the main source of multifractality of the series before and after the law in Seoul and Daegu are the fat-tailed distributions. The second row in Figure 14 shows that the multifractal spectrum of the shuffled series is the narrowest, which implies that the long-range correlation results in multifractality of the series before and after implementation of the law in Busan.

5. Conclusions

In this study, we examined the multifractal characteristics of and time series of all seasons. We determined that multifractality existed in each season, and the multifractal property of was stronger than that of in each season. We also validated the effectiveness of the “haze special law,” which was implemented to improve the governance of air pollution in South Korea. We checked whether there is a difference caused by the law for these cities and calculated confidence intervals of Hurst exponent . The results showed that the confidence intervals do not overlap for Seoul, while the confidence intervals overlap a lot for Busan and Daegu, indicating that there is a difference caused by the law for Seoul; however, the same conclusion is inconclusive for Busan and Daegu. By comparing the Hurst exponent, Renyi exponent, and Holder exponent, the time series of the concentration before the implementation of the law was found to have a higher multifractal degree, which decreased after the enactment of the law. This phenomenon reflected the fact that the law has played a role in improving the efficiency of air pollution control in South Korea. We also concluded that the effect of the law will be more significant in a larger city. To explore the major causes of multifractality, we shuffled and phase-randomized the original series of . By analyzing the width of the multifractal spectrum , the results showed that the fat-tailed distribution contributed to the multifractality of both times series before and after the implementation of the “haze special law” in Seoul and Daegu, whereas long-range correlations resulted in the multifractality of the series before and after the implementation of the law in Busan. We concluded that the implementation of the law was very successful and effective for the improvement of levels. Therefore, we believe that the “haze special law” is a possible cause for the change in API in South Korea.

Data Availability

The data will be available upon request with the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The first author (Jian Wang) was supported by the China Scholarship Council (201808260026). J. S. Kim expresses gratitude for the support from the BK21 PLUS program. The corresponding author Wei Shao thanks Professor Darae Jeong for helpful discussion.