Spatial Characteristics Analysis for Coupling Strength among Air Pollutants during a Severe Haze Period in Zhengzhou, China

This paper investigates the multifractal characteristics of six air pollutants using the coupling detrended fluctuation analysis method. The results show that coupling correlations exist among the air pollutants and have multifractal characteristics. The sources of multifractality are identified using the chi square test. The coupling strengths between different pollutants are quantified. In addition, the coupling contribution of a series in the haze system is calculated, and SO2, as the main pollutant, plays a key role in the pollution system. Moreover, the Kriging interpolation method is used to analyze the spatial characteristic on coupling contribution of SO2. The spatial analysis of coupling strength for air pollutants will provide an effective approach for pollution control.


Introduction
With the rapid development of China's economy, industrialization and urbanization are progressing rapidly, and environmental problems have become increasingly serious in China since the 21st century. Especially, the persistent haze appears in a large scale and high frequency [1,2] in recent years, which has become an important environmental problem. The haze pollution has aroused widespread attention in the scientific community [3].
Some studies explored the cross-correlations between non-stationary financial time series using the detrended cross-correlation method (DCCA) proposed by Podobnik et al. (2008) [18]. Subsequently, Zhou et al. (2008) [19] combined the multifractal detrended fluctuation analysis and DCCA method (the combination of which is known as the MF-DCCA method) to investigate multifractal cross-correlations. The MF-DCCA method, based on q-order detrended covariance [20][21][22][23], is a combination of multifractal analysis and detrended cross-correlation analysis. The DCCA and MF-DCCA methods have been widely applied in many fields, such as financial data [24,25], traffic flow [26,27], river 1. In this paper, the coupling relationship of six pollutants in a severe haze in Zhengzhou City from 17 December 2016 to 26 December 2016 was studied by using CDFA method. The coupling strength between different pollutants are quantified in the haze system, and the contribution of one of the series in the coupling of the others, are quantified and analyzed.

2.
The coupling contribution of the main pollutants in the haze system is studied via spatial interpolation, and the spatial characteristics of the pollutants are analyzed according to the spatial distribution. The spatial analysis of coupling intensity of air pollutants provides an effective method for exploring the causes, so as to provide some scientific basis for pollution control. As the provincial capital of Henan Province, Zhengzhou is an important city on the new Eurasian Continental Bridge, a national central city and a historical and cultural city. The atmospheric environment in Zhengzhou presents a significant compound pollution situation in recent years. Under the combined effects of industrial pollution, motor vehicle emission, building activities and life area source, etc., fine particle pollution in the atmosphere in Zhengzhou is increasing. In 2016, severe pollution weathers frequently happened in Zhengzhou, the scope of regional pollution has expanded, and the characteristics of compound pollution in large cities have been becoming increasingly prominent. What's more, air pollution is more severe during the winter heating season, from late 2016 to early 2017, the severe regional haze pollution for many days caused by high concentrations of pollutants happened again and again. The haze pollution can cause pulmonary sclerosis, asthma and bronchitis, even cardiovascular disease. The general public have attached great importance to haze pollution. How to propose effective control measures on haze weather is a difficult problem to be solved at present. More profound study is needed on the prevention and control of air pollution in Zhengzhou.
There are nine automatic monitoring stations (as seen in Figure 1), including the Gang Li Reservoir (GLR) located in Huiji district, Water Supply Company (WSC) located in Gaoxin district, First Affiliated Hospital of Zhengzhou University (FAHZU) located in Erqi district, Administrative Committee of Economic Development Zone (ACEDZ) located in Jingkai district, Monitoring Station of Zhengzhou (MSZZ) located in Zhongyuan As the provincial capital of Henan Province, Zhengzhou is an important city on the new Eurasian Continental Bridge, a national central city and a historical and cultural city. The atmospheric environment in Zhengzhou presents a significant compound pollution situation in recent years. Under the combined effects of industrial pollution, motor vehicle emission, building activities and life area source, etc., fine particle pollution in the atmosphere in Zhengzhou is increasing. In 2016, severe pollution weathers frequently happened in Zhengzhou, the scope of regional pollution has expanded, and the characteristics of compound pollution in large cities have been becoming increasingly prominent. What's more, air pollution is more severe during the winter heating season, from late 2016 to early 2017, the severe regional haze pollution for many days caused by high concentrations of pollutants happened again and again. The haze pollution can cause pulmonary sclerosis, asthma and bronchitis, even cardiovascular disease. The general public have attached great importance to haze pollution. How to propose effective control measures on haze weather is a difficult problem to be solved at present. More profound study is needed on the prevention and control of air pollution in Zhengzhou.
There are nine automatic monitoring stations (as seen in Figure 1), including the Gang Li Reservoir (GLR) located in Huiji district, Water Supply Company (WSC) located in Gaoxin district, First Affiliated Hospital of Zhengzhou University (FAHZU) located in Erqi district, Administrative Committee of Economic Development Zone (ACEDZ) located in Jingkai district, Monitoring Station of Zhengzhou (MSZZ) located in Zhongyuan district, The 47th Middle School (MS47) located in Zhengdong New district, Tobacco Factory (TF) located in Guancheng district, Banking Institute (BI) and Zhengzhou Textile Machinery Limited Company (ZTMLC) located in Jinshui district.
Near the MSZZ and TF stations are Zhengzhou railway station, Zhengzhou residential and traffic intensive areas, where there is a thermal power plant next to the MSZZ station. The TF, BI and MS47 stations are mainly surrounded by schools, parks and residential areas, and they are close to the heat source plant. Near the ZTMLC station, there are some enterprises with serious pollution, such as a porcelain factory, Hengtai aluminum alloy factory, Zhengzhou first steel plant and Yuxiang electric power engineering company, etc. Around the FAHZU and BI stations, there are mainly residential areas, hospitals and schools, and they are far away from the industrial pollution sources. The GLR station is mainly surrounded by reservoirs and Yellow River scenic spots. The WSC station is close to the industrial zone and heat source plant in the western of Zhengzhou. Around the ACEDZ station, there are mainly residential areas and some industrial areas, and it is close to the famous Zhongyuan futa. These urban monitoring stations are evenly distributed in the important areas of Zhengzhou, and the pollution data from that can well reveal the whole pollution situation about Zhengzhou.
The hourly data of NO 2 , CO, O 3 , PM 10 , PM 2.5 and SO 2 in Zhengzhou at nine automatic monitoring stations from 24:00 17 December 2016 to 23:00 26 December 2016, are supplied by the China National Environmental Monitoring Centre [43] in this paper. The missing data have been replaced by the mean between two of the nearest data. In order to learn the general situation during the severe haze period in Zhengzhou, the average value of each air pollutant time series at nine stations is calculated, as shown in Figure 2. As seen from Figure 2, the haze continues for about 9 days (about 240 h), and it is the severest haze pollution in Zhengzhou in 2016, which the average maximum of PM 2.5 is close to 800 µg/m 3 .
district, The 47th Middle School (MS47) located in Zhengdong New district, Tobacco Factory (TF) located in Guancheng district, Banking Institute (BI) and Zhengzhou Textile Machinery Limited Company (ZTMLC) located in Jinshui district.
Near the MSZZ and TF stations are Zhengzhou railway station, Zhengzhou residential and traffic intensive areas, where there is a thermal power plant next to the MSZZ station. The TF, BI and MS47 stations are mainly surrounded by schools, parks and residential areas, and they are close to the heat source plant. Near the ZTMLC station, there are some enterprises with serious pollution, such as a porcelain factory, Hengtai aluminum alloy factory, Zhengzhou first steel plant and Yuxiang electric power engineering company, etc. Around the FAHZU and BI stations, there are mainly residential areas, hospitals and schools, and they are far away from the industrial pollution sources. The GLR station is mainly surrounded by reservoirs and Yellow River scenic spots. The WSC station is close to the industrial zone and heat source plant in the western of Zhengzhou. Around the ACEDZ station, there are mainly residential areas and some industrial areas, and it is close to the famous Zhongyuan futa. These urban monitoring stations are evenly distributed in the important areas of Zhengzhou, and the pollution data from that can well reveal the whole pollution situation about Zhengzhou.
The hourly data of NO2, CO, O3, PM10, PM2.5 and SO2 in Zhengzhou at nine automatic monitoring stations from 24:00 17 December 2016 to 23:00 26 December 2016, are supplied by the China National Environmental Monitoring Centre [43] in this paper. The missing data have been replaced by the mean between two of the nearest data. In order to learn the general situation during the severe haze period in Zhengzhou, the average value of each air pollutant time series at nine stations is calculated, as shown in Figure 2. As seen from Figure

CDFA (Coupling Detrended Fluctuation Analysis) Method
To reveal the multi fractal characteristics of CO, NO2, O3, PM10, PM2.5 and SO2 in Zhengzhou, we use the CDFA approach proposed by Hedayatifar et al. [37]. The method is well suited for testing the mutual coupling of over two series, and it can be described as follows: where i = 1, 2, N, j = 1, 2, n, , as N is often not be an integral multiple of s, cumulative column tail may have partial residual data not for calculation. in order to take into account such remaining data, the partitioning process can be repeated at the tail of the cumulative series, therefore, a total of 2Ns small boxes have been obtained.

CDFA (Coupling Detrended Fluctuation Analysis) Method
To reveal the multi fractal characteristics of CO, NO 2 , O 3 , PM 10 , PM 2.5 and SO 2 in Zhengzhou, we use the CDFA approach proposed by Hedayatifar et al. [37]. The method is well suited for testing the mutual coupling of over two series, and it can be described as follows:

1.
Series containing n time series {x 1 m , . . . , x j m , . . . , x n m } in N is considered, and, m is the member of each time series. Constructing cumulative series as:

2.
Cumulative series X j (i) is divided into N boxes not overlapping each other in s length, and N s = int(N/s), as N is often not be an integral multiple of s, cumulative column tail may have partial residual data not for calculation. in order to take into account such remaining data, the partitioning process can be repeated at the tail of the cumulative series, therefore, a total of 2N s small boxes have been obtained.

3.
Then, the least square method is used to fit the cumulative series in each small box to gain the local trend x j λ (i) of 2Ns boxes. Then the detrended multicovariance is calculated as The qth order fluctuation function is calculated as follows: 6 of 19

5.
If time series x 1 m , . . . , x j m , . . . , x n m is involved in long-range power law correlation, the F x 1 m ,...,x n m (q, s) will increase as a power law of s.
where h x 1 m ,...,x n m (q) is the generalized CDFA scaling exponent for all the time series, which can be obtained by calculating the slope of the log-log plots of F x 1 m ,...,x n m (q, s) versus s.
When n = 1, q = 2, it is simplified to the DFA method, which can be used to study the long-term persistence characteristics of only one time series. When n = 2, q = 2, it is known as the DCCA method, and the method can be used to analyze the correlation between two time series. If h x 1 m ,...,x n m (q) = c, c is a constant, then the coupling relation is monofractal, while its value changing along with q indicates the presence of multifractality. More specifically, if h x 1 m ,...,x n m (q) > 0.5, the coupled correlation is positive, and it means that a large (or small) increment in one variable is more likely to be followed by a large (or small) increment in other variables. If h x 1 m ,...,x n m (q) < 0.5, the coupling correlation is negative, i.e., anti-correlated, which means that a large (or small) increment in one variable is more likely to be followed by a large (or small) decrement in other variables. If h x 1 m ,...,x n m (q) = 0.5, there is no coupling correlation among time series, and the changes in one variable cannot affect the behavior of other variables.
In addition, for q > 0, h(q) describes the scaling behaviour of the segments with large fluctuations. For q < 0, h(q) describes the scaling behaviour of segments with small fluctuations. Therefor the multifractality gives information about the relative importance of various fractal exponents present in the time series. The strength of the multifractality of a time series can be characterized by ∆h = h max (q) − h min (q), where ∆h is the range of generalized Hurst exponents h(q). The larger the ∆h is, the stronger the multifractality degree become, and vice versa.

Multifractal Contribution Sources
Researchers have shown interest in the nature of multifractal behavior in time series. There are two multifractal contribution sources which can be found in the time series [44,45]. One is the cross-correlation of different fluctuations, and the other is the variable fat-tailed probability distribution [46]. The two sources can be realized using a shuffling procedure (shuffled data) and phase randomization procedure (surrogate data), so as to study the multifractal characteristics of time series. The shuffling procedure can eliminate crosscorrelation of original data, and keep the original data distribution unchanged. Thus, the shuffled data can be used to study multifractal contribution from cross-correlation in time series. On the contrary, the phase randomization procedure can eliminate non-Gaussian distribution and only retain correlation of original data. Thus, the surrogate data can be used to study the multifractal contribution of variable fat-tailed probability distribution in time series.
According to the chi square statistical test method proposed by Hedayatifar et al. [37], We can quantify the coupling strength of long-term memory or fat-tail probability as follows: The " " symbol represents "surrogate" or "shuffled", and the Y means that the series has been phase randomized or shuffled. We can use the following local Relation (7) to quantify the contribution of sequence y in the other coupling. where The "all" specifies all the series are shuffled, and the σ(q) is the generalized standard deviation for q.

Coupling Detrended Fluctuation Analysis (CDFA) of Air Pollutants Series
In order to study the multifractal characteristics and quantify the coupling correlations of CO, NO 2 , O 3 , PM 10 , PM 2.5 and SO 2 , we first need to process the original data of six pollutants by shuffling and surrogate procedures, respectively. Figure 3 shows the result surfaces obtained from Formula (4)  it proves the conclusion mentioned in Section 3 that the shuffling procedure eliminates cross-correlation of original data and keeps original data distribution unchanged. In addition, the surface in Figure 3c becomes more smoother than that in Figure 3a,b, and it means that the surrogate procedure eliminate non-Gaussian distribution and retain correlation of original data. Based on the above conclusions and analysis, the shuffled and surrogate processing for the six pollutants original data in this section are very effective. Additionally, the shuffled and surrogate data can meet the requirements of CDFA multifractal analysis in this paper well.

Analysis on Multifractal Contribution Sources of Pollution Series
According to Section 3.2, there are two sources, i.e., cross-correlation and variable fattailed probability distribution, which can affect the multifractal characteristics of time series. However, for a specific system, the two sources may not all have an important influence on the multifractal. So, we first analyze which source has a greater impact on the It can be seen that the surface of Figure 3b is more wrinkled and less smooth than that of Figure 3a. And as seen from Figure 3b, the value of h x 1 m ,...,x n m (q), which can be obtained by calculating the slope of ln(F x 1 m ,..., n m (q, s)) versus ln(s) plots, is near 0.5, and it proves the conclusion mentioned in Section 3 that the shuffling procedure eliminates cross-correlation of original data and keeps original data distribution unchanged. In addition, the surface in Figure 3c becomes more smoother than that in Figure 3a,b, and it means that the surrogate procedure eliminate non-Gaussian distribution and retain correlation of original data.
Based on the above conclusions and analysis, the shuffled and surrogate processing for the six pollutants original data in this section are very effective. Additionally, the shuffled and surrogate data can meet the requirements of CDFA multifractal analysis in this paper well.

Analysis on Multifractal Contribution Sources of Pollution Series
According to Section 3.2, there are two sources, i.e., cross-correlation and variable fat-tailed probability distribution, which can affect the multifractal characteristics of time series. However, for a specific system, the two sources may not all have an important influence on the multifractal. So, we first analyze which source has a greater impact on the haze pollution system studied in this paper, and we also prepare for the analysis on coupling correlation and contribution of six pollutants in next section.
The steps of calculation and analysis in this section are: First, the time series of CO, NO 2 , O 3 , PM 10 , PM 2.5 and SO 2 at nine stations are processed by shuffling and phase randomization procedures, respectively, then we can obtain the shuffled and surrogate series of them. Second, Using CDFA method to calculate h x 1 m ,...,x n m (q) (simplify for h(q)) values versus q of the six pollutants series based on the Formula (4) in Section 3.1, and plot the q~h(q) curves as shown in Figure 4. Finally, by comparing the changes of shuffled and surrogate series curves relative to the original series curve, it can be concluded that the curve with greater changes is the source with greater impact on multifractal. As shown in Figure 4, the CDFA analysis of the shuffled series is significantly different from that of the original series at nine stations. It shows that the shuffling procedure is effective in cleaning data correlation, and the h(q) value of shuffled series is close to 0.5. This implies that almost all the series after being shuffled have no coupling correlation, and shows again that it is very effective to study the coupling strength and contribution of each pollutant by shuffling procedure.  In order to more accurately study the influence of different sources on multifractal, the behaviors of q~h(q) of original series, shuffled series and surrogate series of six pollutants at nine stations are given using CDFA method in Figure 4, and the q~h(q) curves are calculated by the Formula (4). The black line with " " is the CDFA analysis on the original series of six pollutants. The red line with " NO2, O3, PM10, PM2.5 and SO2 at nine stations are processed by shuffling and phase randomization procedures, respectively, then we can obtain the shuffled and surrogate series of them. Second, Using CDFA method to calculate 1 ,..., ( ) n m m x x h q (simplify for h(q)) values versus q of the six pollutants series based on the Formula (4) in Section 3.1, and plot the q~h(q) curves as shown in Figure 4. Finally, by comparing the changes of shuffled and surrogate series curves relative to the original series curve, it can be concluded that the curve with greater changes is the source with greater impact on multifractal.
In order to more accurately study the influence of different sources on multifractal, the behaviors of q~h(q) of original series, shuffled series and surrogate series of six pollutants at nine stations are given using CDFA method in Figure 4, and the q~h(q) curves are calculated by the Formula (4). The black line with "★" is the CDFA analysis on the original series of six pollutants. The red line with "◆" is the CDFA analysis on the shuffled series of six pollutants, and the blue line with "•"is the CDFA analysis on the surrogate series of six pollutants.
" is the CDFA analysis on the shuffled series of six pollutants, and the blue line with "•" is the CDFA analysis on the surrogate series of six pollutants.
As shown in Figure 4, the CDFA analysis of the shuffled series is significantly different from that of the original series at nine stations. It shows that the shuffling procedure is effective in cleaning data correlation, and the h(q) value of shuffled series is close to 0.5. This implies that almost all the series after being shuffled have no coupling correlation, and shows again that it is very effective to study the coupling strength and contribution of each pollutant by shuffling procedure.
As shown in Figure 4, the behaviors of q~h(q) for shuffled and surrogate series are different from the original series. All the q~h(q) curves of surrogate series are close to the original series; however, the q~h(q) curves of shuffled series are far away from the original series. It is obvious to see from Figure 4 that there is h shu f f (q) < h surr (q) relationship for each monitoring station. This means that during the period of severe haze, the mutual coupling of pollutants time series is not only affected by cross-correlation, but also affected by the leptokurtic and fat tailed feature. The cross-correlation play a dominant role in the coupling relationship among the pollutants time series. The coupling correlation between pollutants is mainly caused by the large and small fluctuations, while the leptokurtic and fat tailed effect for the peak value are not large. Therefore, the influence of cross-correlation for coupling strength of pollutants is analyzed in the next section.

Coupling Correlation and Contribution of Pollution Series in the Haze System
From the analysis in the above section, we have learned that the multifractal source of this severe haze pollution in Zhengzhou is cross-correlation. Let's take SO 2 as an example to explain the analysis process of the CDFA method. If we want to study the coupling contribution of SO 2 in the haze system using the CDFA method, we can take the following steps: (1) Processing original time series of SO 2 (named as Origin_SO 2 ) by shuffling procedure, and we obtain the shuffled time series of SO 2 (named as Shuff_SO 2 ). (2) Substituting the Shuff_SO 2 , Origin_CO, Origin_NO 2 , Origin_O 3 , Origin_PM 10 , Origin_PM 2. 5 into Formula (4) of Section 3.1, and we calculate the value of h(q) versus q. (3) Plotting the q~h(q) curve, we name the curve as CDFA_SO 2 , which represents the generalized scaling exponents of CDFA when only SO 2 is shuffled. (CDFA_CO, CDFA_NO 2 , CDFA_O 3 , CDFA_PM 10 , and CDFA_ PM 2.5 can be calculated in the same way). (4) Replace the Shuff_SO 2 with Origin_SO 2 in step (2) above, and we calculate the q~h(q) again. Plotting the new q~h(q) curve which named CDFA_Origin, it denotes the generalized scaling exponents of CDFA for all the original time series. (5) The coupling contribution of SO 2 (or CO, NO 2 , O 3 , PM 10 , PM 2.5 ) can be analyzed by using the relevant parameters in Section 3 finally.
We calculate and plot the q~h(q) curves of CDFA_Origin, CDFA_CO, CDFA_NO 2 , CDFA_O 3 , CDFA_PM 10 , and CDFA_PM 2.5 at nine stations, respectively, according to the above steps as seen in Figure 5. It shows the behaviors of h(q) when applying the CDFA for CO, NO 2 , O 3 , PM 10 , PM 2.5 and SO 2 time series, only the mentioned time series is shuffled and the other original time series are kept unchanged. By using the shuffling procedure, only the correlation effect is washed out and the probability density function effect exists. From Figure 5, we can see that the h(q) value changing along with q indicates the presence of multifractality. Meanwhile, the shuffled series shows distinct differences for 0 q < and 0 q > , that is, the coupling strength among the pollution series shows distinct differences for the small fluctuations ( 0 q < ) and the large fluctuations ( 0 q > ) in the original series. For example, in the GLR and ACEDZ stations, when 0 q < , the O3 is the farthest pollutant from the original series after being shuffled, which illustrates that it has the largest influence in small fluctuations of the haze system; when 0 q > , it is the nearest pollutant from the original series after being shuffled. This illustrates that O3 has the least influence in large fluctuations of the haze system. For SO2, the influence of SO2 in 0 q > is obvious greater than that in 0 q < , and this indicates that the SO2 has a greater influence in the large fluctuations, and has a less influence in the small fluctuations of the haze system. The similar conclusions can be obtained for other pollutants as well.
We can quantify the coupling strength and contribution of cross-correlation by calculating the values of 2 shuff χ and 2 shuff χ make use of Formulas (5) and (7). Table 1  respectively, that is, the CO and SO2 have the strongest coupling correlation and largest coupling contribution, as observed from  From Figure 5, we can see that the h(q) value changing along with q indicates the presence of multifractality. Meanwhile, the shuffled series shows distinct differences for q < 0 and q > 0, that is, the coupling strength among the pollution series shows distinct differences for the small fluctuations (q < 0) and the large fluctuations (q > 0) in the original series. For example, in the GLR and ACEDZ stations, when q < 0, the O 3 is the farthest pollutant from the original series after being shuffled, which illustrates that it has the largest influence in small fluctuations of the haze system; when q > 0, it is the nearest pollutant from the original series after being shuffled. This illustrates that O 3 has the least influence in large fluctuations of the haze system. For SO 2 , the influence of SO 2 in q > 0 is obvious greater than that in q < 0, and this indicates that the SO 2 has a greater influence in the large fluctuations, and has a less influence in the small fluctuations of the haze system. The similar conclusions can be obtained for other pollutants as well.
We can quantify the coupling strength and contribution of cross-correlation by calculating the values of χ 2 shu f f and χ 2 shu f f make use of Formulas (5) and (7). Table 1 gives the average of χ 2 shu f f and χ 2 shu f f values for six pollutants in nine stations. To be convenient for analysis, the values are expressed by the high-dimension graph for 0 ≤ q ≤ 20, as shown in Figure 6. during the periods from 2012 to 2015 in Zhengzhou, and Chemical Mass Balance results show that the contributions of major sources (i.e., nitrate, sulfate, biomass, carbon and refractory material, coal combustion, soil dust, vehicle, and industry) of PM2.5 are 13%, 16%, 12%, 2%, 14%, 8%, 7%, and 8% in heavily polluted days, with sulfate and nitrate as major ions. In 2016, the sum of the concentrations of SO , NO 3− , and NH 4+ increase with the aggravation of pollution level, and the secondary aerosols is the main source of PM2.5, which accounts for 38.4% on heavy pollution days in Zhengzhou [49].

Spatial Characteristics Analysis on Coupling Strength of Air Pollutants
Through the above analysis on coupling parameter values in each station, it can be found that they have obvious spatial characteristics. Therefore, it is necessary for us to make spatial analysis on pollutant coupling strength parameters to mine more useful conclusions. The coupling correlations and contributions of SO2 is the greatest for this haze among the four gaseous pollutants except PM10 and PM2.5. The PM10 and PM2.5 are fine particles, and their formation are diverse and complex. However, SO2 mainly comes from the direct emission of pollution sources, and its emission is also a cause of the formation As seen from Figure 6a,b, in the ACEDZ, MSZZ and MS47 stations, the χ 2 shu f f and χ 2 shu f f values of PM 10 are the largest, which indicates that the coupling correlation and contribution of PM 10 with other pollutants in large fluctuation is the greatest. In the GLR and BI stations, the pollutants with the largest χ 2 shu f f and χ 2 shu f f values are CO and SO 2 respectively, that is, the CO and SO 2 have the strongest coupling correlation and largest coupling contribution, as observed from Table 1 Based on the above analysis, SO 2 has the greatest impact in the haze system other than PM 10 and PM 2.5 . SO 2 not only pollutes the atmosphere in gaseous form, but more than 30% of it can be converted into sulfate aerosols to cause an increase of atmospheric particulate concentration. It is also a critical precursor reactant for secondary aerosol. In the winter, the increase of coal emissions results the increase of SO 2 emissions, while the SO 2 is the key gaseous precursor of secondary inorganic aerosol, so it directly affects the large fluctuation part in the haze system, such as the violent fluctuations of PM 2.5 . Therefore, the SO 2 has the largest impact and plays a dominant role in the large fluctuation part of the haze coupling system except PM 10 and PM 2.5 . The finding is consistent with the source apportionment experiment of PM 2.5 in Zhengzhou [47][48][49]. In 2010, the dominant components of PM 2.5 in Zhengzhou are secondary ions (sulfate and nitrate) and carbon fractions, which soluble ions and total carbon contribute 41% and 13% of PM 2.5 mass, respectively [47]. Jiang et al. (2017) [48] conduct the sample source apportionment experiment of PM 2.5 during the periods from 2012 to 2015 in Zhengzhou, and Chemical Mass Balance results show that the contributions of major sources (i.e., nitrate, sulfate, biomass, carbon and refractory material, coal combustion, soil dust, vehicle, and industry) of PM 2.5 are 13%, 16%, 12%, 2%, 14%, 8%, 7%, and 8% in heavily polluted days, with sulfate and nitrate as major ions. In 2016, the sum of the concentrations of SO 2− 4 , NO 3− , and NH 4+ increase with the aggravation of pollution level, and the secondary aerosols is the main source of PM 2.5 , which accounts for 38.4% on heavy pollution days in Zhengzhou [49].

Spatial Characteristics Analysis on Coupling Strength of Air Pollutants
Through the above analysis on coupling parameter values in each station, it can be found that they have obvious spatial characteristics. Therefore, it is necessary for us to make spatial analysis on pollutant coupling strength parameters to mine more useful conclusions. The coupling correlations and contributions of SO 2 is the greatest for this haze among the four gaseous pollutants except PM 10 and PM 2.5 . The PM 10 and PM 2.5 are fine particles, and their formation are diverse and complex. However, SO 2 mainly comes from the direct emission of pollution sources, and its emission is also a cause of the formation for PM 10 and PM 2.5 . SO we think it is more meaningful to study the spatial characteristics of SO 2 coupling strength.
In the severe haze of Zhengzhou, the cross-correlation has stronger impact than the fat tailed distribution on multifractal characteristics for pollutants. In addition, the coupling contribution is good enough to characterize the coupling strength of pollutants in the haze power system, so we use ArcGIS software to analyze the spatial distribution of the average concentration and coupling contribution χ 2 shu f f value of SO 2 in the nine monitoring stations by using Kriging interpolation. The spatial distribution of the average concentration and χ 2 shu f f value of SO 2 in Zhengzhou city by contour interpolation method when 0 ≤ q ≤ 20 is given in the Figure 7a  for PM10 and PM2.5. SO we think it is more meaningful to study the spatial characteristics of SO2 coupling strength.
In the severe haze of Zhengzhou, the cross-correlation has stronger impact than the fat tailed distribution on multifractal characteristics for pollutants. In addition, the coupling contribution is good enough to characterize the coupling strength of pollutants in the haze power system, so we use ArcGIS software to analyze the spatial distribution of the average concentration and coupling contribution 2  As seen in Figure 7a, the spatial distribution of SO2 concentration in Zhengzhou presents decreasing trends from the central zone to north and south ends. Firstly, the west and east of the city are factory concentrated areas, and factory exhaust emissions are a major source of SO2 pollutant; secondly, the area is also a densely populated area, and the winter is a heating period. Therefore, the pollution gases emitted by traffic, coal burning -are also an important source of SO2; Finally, from Figure 7a, it can be seen that the production of SO2 should be closely related to the distribution of thermal power plants, heat source plants and waste disposal plants. Throughout the western region, especially in the ZTMLC, there exists maximum SO2 concentration. On the one hand, the western region is an industrial concentration area, on the other hand, the heavily polluting enterprises, such as City porcelain plant, City Hengtai aluminum alloy factory, Zhengzhou No.1 steel mill, As seen in Figure 7a, the spatial distribution of SO 2 concentration in Zhengzhou presents decreasing trends from the central zone to north and south ends. Firstly, the west and east of the city are factory concentrated areas, and factory exhaust emissions are a major source of SO 2 pollutant; secondly, the area is also a densely populated area, and the winter is a heating period. Therefore, the pollution gases emitted by traffic, coal burning -are also an important source of SO 2 ; Finally, from Figure 7a, it can be seen that the production of SO 2 should be closely related to the distribution of thermal power plants, heat source plants and waste disposal plants. Throughout the western region, especially in the ZTMLC, there exists maximum SO 2 concentration. On the one hand, the western region is an industrial concentration area, on the other hand, the heavily polluting enterprises, such as City porcelain plant, City Hengtai aluminum alloy factory, Zhengzhou No.1 steel mill, Yuxiang power engineering company are distributed near the ZTMLC, and this is also the direct cause of the serious SO 2 pollution. Although the actual concentration of SO 2 in the eastern and northern parts of the city is not small as well, this may be related to winter heating by coal combustion and motor vehicle exhaust emissions, and the emission sources are relatively less than those in the western region, so the SO 2 concentration is relatively low here.
As seen from Figure 7b, the spatial distribution of SO 2 coupling contribution in Zhengzhou presents a decreasing trend from central and southeast regions (FAHZU, BI and ZTMLC) to northwest region. The spatial distribution is different from that of SO 2 concentration, and this indicates that the concentration of SO 2 is high, and its coupling contribution may not be necessarily large, which may also be affected by the weather and chemical reaction between pollutants. That is, there may be many factors influence the coupling chemical reaction of SO 2 . As seen from the Figure 7b, in the vicinity of ZTMLC, the concentration of SO 2 is very large where coupled chemical reactions are more likely to happen, thus resulting in the relatively large coupling contribution in this area. In the vicinity of FAHZU and BI, the actual concentration value is relatively small, because it is the main residential areas, hospitals and schools, away from industrial pollution sources. However, nearby the FAHZU, as seen from Figure 6b, the coupling contribution of NO 2 is the least, the regional haze is mainly caused by SO 2 , and the SO 2 plays a main role in the haze system. As seen also from Figure 6b, the maximum value of coupling contributions for NO 2 , SO 2 , PM 2.5 and PM 10 appear in the BI station, which indicates that in this area these pollutants are involved in a series of coupling chemical reactions. The pollutants of NO 2 , SO 2 and O 2 produce nitric acid, sulfuric acid and O 3 under the light and chemical actions, so the coupling contribution of O 3 is maximum. What's more, the nitric acid and sulfuric acid produced by SO 2 form aerosols, then convert to PM 2.5 and PM 10 , therefore, the coupling contribution of both PM 2.5 and PM 10 have a local maximum. During this haze, the ground wind in winter mainly tends north and northwest in Zhengzhou, and it has a positive effect on the coupling of SO 2 with other pollutants in the Southeast. It makes the SO 2 have a great coupling contribution in the south central Zhengzhou. In the GLR station, the coupling contribution of SO 2 is minimum. This is mainly because the area around the GLR station are reservoirs, as well as the Yellow River scenic spots, and there are no sources of SO 2 emissions. As seen from Figure 7a, the SO 2 concentration is also the lowest due to the small coupling contribution of SO 2 in the station.

Conclusions
In this paper, the coupling relationship among six pollutants (CO, NO 2 , O 3 , PM 10 , PM 2.5 and SO 2 ) in a severe haze weather from 17 December 2016 to 26 December 2016 in Zhengzhou is analyzed by the CDFA method, as well as the effects of various factors on city haze pollution. The CDFA method can be used to quantify the strength of the coupling correlations and their coupling contributions of different pollutants in haze weather, thus a detailed understanding on the role and status of different pollutants in the haze power system can be gained.
From the experimental results, we can see that the SO 2 has the greatest influence in the large fluctuation part of haze coupling system. It is concluded that the primary pollutant of SO 2 plays a very important role in the typical haze in Zhengzhou, and it is the main contribution to the haze system. This is consistent with that Zhengzhou City belongs to SO 2 Pollution Control Zone. Thus, the coupling contribution of SO 2 in the haze system is studied by spatial distribution using Kriging interpolation, and the spatial characteristics of the pollutants are analyzed according to the spatial distribution.
The spatial analysis on coupling intensity of air pollutants will provide a good approach for analyzing pollution sources. It is noteworthy that there is an increasing trend for the NO 2 pollution in 2016, and the contribution of NO 2 pollutants to the haze cannot be ignored. The spatial interpolation on coupling contribution of pollutants in each station is studied. It is a new attempt and innovation to study the distribution of coupling contribution of pollutants in each station by spatial interpolation method. Through spatial analysis, we can obtain the spatial distribution of pollutants concentration or relevant parameter values, and further analyze the causes of air pollution in different regions in combination with other influencing factors, such as urban planning, geographical conditions or population density, so as to provide some scientific basis for prevention and control of air pollution. It provides an effective method for exploring the causes of air pollution.
Due to the recurrence of the COVID-19 in the past two years, people have greatly reduced their production and living activities, and affected by which, the air quality is significantly improved, but the economic development is hit hard. Therefore, how to reduce air pollution as much as possible while ensuring economic development is the focus of our attention. The economic development mode of China will not change greatly in a short time, so we are committed to finding a better way to study the mechanism of air pollution, and make our own efforts for improving the atmospheric environment. In the future, whether the air pollution has an impact on extreme weather deserves our attention. And whether the continuous air pollution in Zhengzhou has an impact on the 7.20 severe rainstorm last year, is a subject worthy of study also.