Backward Trajectory andMultifractal Analysis of Air Pollution in Zhengzhou Region of China

With the continuous promotion of industrialization and urbanization, China’s environmental pollution is becoming increasingly serious, which has caused considerable damage to the natural balance. Air pollution seriously harms people’s physical and mental health, the ecological environment, and the social sustainable development of society. In this study, the backward trajectory model and multifractal methods were adopted to analyze air pollution in Zhengzhou./e backward trajectory analysis showed that most clusters of air pollution were from southern Hebei, eastern Shandong, and mid-western Henan, which were then transported to Zhengzhou. For the PSCF and CWT analyses, we selected four representative cities to explore how close the air pollution of Zhengzhou is to other areas on the basis of air polluted concentration. /e results of several multifractal methods indicated that multifractality existed in the AQI time series of Zhengzhou and cross-correlations between Zhengzhou and each of the four cities. /e widths of multifractal spectra showed that the air pollution in Zhengzhou was closest to that in Jinan, followed by Shijiazhuang, Zibo, and Luoyang. /e CDFA analysis showed that carbon monoxide (CO), nitrogen dioxide (NO2), and inhalable particulate matter (PM10) had important influences on air pollution in Zhengzhou. /ese findings offer a useful reference for air pollution sources and their potential contributions in Zhengzhou, which can support policy makers in environmental governance and in achieving sustainable urban development.


Introduction
Aerosols are colloids of fine solid particles or liquid droplets suspended in the atmosphere that affect climate, weather, and air quality through a variety of mechanisms. Aerosols can significantly affect climate through the effects of aerosolradiation interactions (ARIs), aerosol-cloud interactions (ACIs), or both. Strong interaction between aerosols and meteorological variables occur in the planetary boundary layer (PBL). e PBL involves a number of chemical, physical, and dynamic processes that are closely related to the accumulation of aerosol pollution, which has a negative impact on the climate and environment and poses a threat to human health [1][2][3].For example, exposures to aeroallergens may cause allergy symptoms, especially for sensitized patients. Many people may be allergic to pollen, leading to seasonal rhinitis or skin diseases. Moreover, the transport of the pollen from distant sources may result in allergy symptoms in sensitized patients after the flowering periods [4]. e main factors that cause air pollution are natural sources and human activities. Natural sources mainly include plant pollen, soil dust, forest fire, and volcanic eruption, whereas human activities mainly include burning fuel, industrial emissions, and vehicle emissions. Dust storms (DSs) are natural phenomena that are prevalent in the areas with loose sands, especially in the Middle East, the Southwestern United States, Northern China, and the Saharan Desert. DSs cause not only environmental and human health problems but also decrease in visibility, which limits various activities and increases traffic accidents [5].
As for human activities, Bari and Kindzierski [6] investigated air quality characteristics and potential sources affecting PM 2.5 concentrations in the community of Hinton, Alberta, Canada. ey found that local industry emissions such as Hinton pulp and wood products had moderate influences from the southwest direction. Shipping emissions were regarded as an important factor associated with particulate matter, and the potential of offshore activities and shipping emissions may influence the air quality in urban areas located more than 1000 km away from the cities [7]. Anthropogenic activities combined with vehicular emissions influence the mixing ratios of CO and NO 2 due to office/ working hours and the shallowing and widening of boundary layer. NO 2 and CO pollutants from urban traffic help increase O 3 , which results in higher O 3 mixing ratios during the rush hours and daylight period. e planetary boundary layer (PBL) is crucial in determining O 3 mixing ratios. e PBL is shallow and wide during the winter and post-monsoon season, which leads to the stagnation of gaseous pollutants due to low atmospheric temperature. During summer, PBL is deep, with frequent mixing of the pollutants causing dilution of the pollutants due to higher solar radiation. ese mixing ratios indicate that the highest mean mixing ratios of O 3 during the day are greater than at night [8]. Outdoor barbecues may have a potential impact on the air quality in a city. Atmospheric aerosols produced by the burning of mixed hardwoods and agricultural wastes influence environmental conditions in another region [9]. With the continuous development of industrialization and urbanization in China, the environmental deterioration caused by aerosol pollution has been increasingly significant, seriously harming human health and the ecological environment. e huge economic advances achieved at the expense of the environment have undoubtedly increased the pressure on the environment, harmed the physical and mental health of people, and incurred large costs in terms of sustainable economic development [10].
In the study of air pollution and human health, scholars have found that air pollution damages the human immune system and respiratory systems [11][12][13]. For example, ozone is a poisonous gas and, if inhaled, can cause oxidative damages and inflammation, leading to neurodegenerative pathologies, especially Alzheimer's disease [12]. Exposure to NO 2 is also associated with leukocyte telomere length, a marker of molecular ageing, which leads to health problems [14]. SO 2 not only poses a serious threat to human health but also damages vegetation and ecosystems [15]. In 2020, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that caused the coronavirus disease (COVID-19) has spread rapidly across countries, leading to large numbers of human deaths and unpredictable economic loss. e severity of the symptoms of those with COVID-19 may be associated with PM exposure [16]: the potential mechanism is that inflammatory pathways in the small respiratory airways are activated when the COVID-19 patients are exposed to PM, which leads to inflammatory cell recruitment and affects the prognosis of COVID-19 patients [17]. Recent studies found that air quality has improved significantly during the lockdowns due to the COVID-19 pandemic. Li et al. [18] reported that changes in the Air Quality Index (AQI) corresponded to different response levels determined by the different levels of administrative restrictions. ey found that the concentrations of different pollutants were not statistically significantly different except for NO 2 during the second-and third-level pandemic responses. Xian et al. [19] found that the synergistic effect of meteorological conditions and heterogeneous chemistry on the particulate pollution formation was still uncertain. erefore, they could not confirm a certain relationship between the COVID-19 and PM exposure. e lockdowns may have had substantial short-term health benefits but did not improve the environment in the long term because local air quality is affected by local sources and transported pollutants [20]. Filonchyk and Peterson [21] confirmed that background sources significantly contributed to air pollution, even if human activities were severely reduced.
Air quality evolves with time and space evolution and so has different characteristics at different times and in different places. Huang et at. [22] found that PM 2.5 had different characteristics at different times and seasons. e daily and monthly PM 2.5 values showed U-shaped patterns, and the seasonal value of PM 2.5 was high in winter but low in summer, with spring and autumn values being in between. e spatial distribution of the PM 2.5 concentration mainly followed a northeast-southwest direction, and the air quality in North China was significantly worse than that in South China. e air quality of a region is affected not only by its own region's activities but also by those in other regions.
is is because the atmospheric circulation processes influence the extent and direction of air mass transfer; hence, air pollutants are transported from one place to another [4]. Numerous studies have been conducted to identify air pollutant transport and trajectories over different regions of the world [23][24][25][26][27][28][29]. Yu et al. [30] used the HYSPLIT model to study the pollution characteristics of 29 cities in the Central Plains Economic Zone (CPEZ). ey found that most serious air pollution may occur when air masses were transported from Shaanxi Province, Hubei Province, Anhui Province, and the northwest of the CPEZ. Hence, air pollution in Zhengzhou may originate from both long-distance transmission and local emissions. Filonchyk and Yan [31] studied the characteristics of six principal atmospheric pollutants in the city of Lanzhou. eir results from backward trajectory analysis showed that the air pollution may be transported from the north, northwest, and west directions. Li et al. [32] adopted the HYSPLIT model to explore air pollution in Qingdao, Shandong province. ey found that the main factors resulting in air pollution were coal combustion and vehicular, domestic, and industrial emissions. Vehicular emissions and liquefied petroleum gas (LPG) and solvent usage also led to air pollution in Xiamen city. e air pollutants from Southeast China were transported to Xiamen via Fuzhou and affected the local environment [33]. Other researchers also found that industrial emissions, vehicular emissions, and trans-boundary pollution were important factors influencing air quality [34,35]. Punsompong and Chantara [36] used the backward trajectory analysis to investigate the percent contribution of emissions from local and trans-boundary to air pollutant concentrations. e results showed that the major factor of increasing the PM 10 concentration in Northern ailand was open biomass burning in Myanmar. Domestic biomass-burning aerosols that were transported from Central ailand and Northeast ailand affected the ambient air in Bangkok, ailand, during the 2018-2019 haze episodes [37]. e long-range transport of the products of biomass burning affected ambient air in other regions [38,39].
Most studies used the potential source contribution function (PSCF) and concentration weighted trajectory (CWT) to identify the sources of air pollution in a region. However, these two methods do not consider the degree of influence of each source on the region. e multifractal method is non-linear and can be used to analyze the degree of influence of one time series on another by investigating the coupling correlations of two time series. Multifractal detrended fluctuation analysis (MF-DFA) proposed by Kantelhardt et al. [40] has been widely adopted to identify the source of time series [41,42]. ere are two main sources of multifractality: long-range correlations and fat-tailed probability distributions of time series. Zhou [43] firstly used the multifractal detrended cross-correlation analysis (MF-DCCA) to explore multifractal cross-correlations. e strength of multifractality can be regarded as an indicator to reflect the closeness of the connection between two time series, which could be used to obtain the degree of influence of one time series on another [44]. e coupling detrended fluctuation analysis (CDFA) proposed by Hedayatifar et al. [45] is an extension of MF-DCCA that can be used to reflect the degree of influence by investigating the coupled correlations among several time series. e multifractal cascade model was used to analyze the multifractal characteristics of ozone concentration time series in Taipei. e results showed that multifractal characteristics existed and other fractal methods should be used to further offer sound proof for the existence of a fractal nature [46]. Multifractal detrended fluctuation analysis (MF-DFA) was employed to analyze the multifractal nature of particulate matter (PM) concentrations from vehicle emissions in Hong Kong. e results indicated that particulate matters showed long-term correlation behaviors. e particulate matters in the range of d p < 5μm had strong multifractal characteristics, whereas the particulate matters in the range of d p ≥ 5μm had weak multifractal characteristics [47]. In addition, Liu et al. [48], Shen et al. [49], and Dong et al. [50] found multifractal characteristics based on the MF-DFA approach in mainland China. ere are a few studies based on the MF-DCCA method in the field of environmental pollution. Zhang et al. [51] used MF-DCCA and multifractal asymmetric detrended cross-correlation analysis (MF-ADCCA) to investigate the cross-correlations of PM 2.5 and meteorological factors. e results showed that the cross-correlations between PM 2.5 concentration and four meteorological factors were multifractal and anti-persistent, and the strength of multifractality in Beijing was stronger than that in Hong Kong. Stan et al. [52] employed the multifractal multiscale detrended cross-correlation analysis (MM-DCCA) to explore the cross-correlations among ozone and some selected precursors and between each precursor and temperature. e results showed that both individual data and correlated series showed multifractal properties and persistent long-range correlations. Hurst exponents, scaling exponents, and multifractal spectra confirmed that the two time series had multifractal characteristics in two time series, which implied the existence of non-linear cross-correlations between respiratory diseases and haze in South Korea. e main source of multifractality was fat-tailed distributions [53].
In this study, we adopted the backward trajectory model and multifractal methods to analyze the air pollution time series in Zhengzhou, Henan Province, China. We use the daily AQI index and six daily average concentrations of PM 2.5, PM 10 , SO 2 , CO, NO 2 , and O 3 in five cities from December 2019 to February 2020. e remainder of this paper is organized as follows. Section 2 introduces the methodology. Section 3 describes the data sources. Section 4 presents the empirical results, and Section 5 provides the conclusions.

HYSPLIT Model and Cluster Analysis.
In order to identify the air pollutant trajectory, we used the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model developed by the National Oceanic and Atmospheric Administration (NOAA), USA. According to the NOAA, the HYSPLIT model uses either puff or particle approaches to calculate trajectories, dispersion, and deposition. e model associates the mass of the pollutant species with the release of either puffs, particles, or a combination of both to calculate air concentrations. Air concentrations are calculated for each grid cell of advection and diffusion, and the advection and diffusion of particles are computed from their initial location. We collected data from the global data assimilation system (GDAS) from the NOAA National Centers for Environmental Prediction (NCEP) to calculate 72hour backward trajectories at a 500 m altitude. Cluster methods of air mass backward trajectories can be used to identify the major transport pathways. e three different approaches for time series clustering are raw-data-based approach, feature-based approach, and model-based approach. Trajectory clustering is a time series clustering problem. In this study, the raw-data-based method was used for the air pollutant trajectory clustering. We employed a K-means cluster analysis based on the Euclidean distance to divide the air pollutant trajectories into five categories.

PSCF and CWT Analyses.
In order to determine potential source regions, we used the potential source contribution function (PSCF) model and concentration weighted trajectory (CWT) model based on the HYSPLIT model. e PSCF is a conditional probability function; it can be calculated as the ratio of the number of trajectory endpoints that are higher than a given criterion value to the number of total endpoints in the grid cell. e method calculates the PSCF value in the grid cell (i, j) as follows:

Mathematical Problems in Engineering
where m ij is the number of trajectory endpoints greater than a given criterion value in the grid cell (i, j) and n ij is the total number of trajectory endpoints falling within the grid cell (i, j). We then used a weight function, W(n ij ), multiplied into the PSCF value to better reflect the uncertainty in the values for these cells. W(n ij ) is calculated as where N is the average number of trajectories points of all grid cells that contain at least one trajectory and n ij is the total number of trajectory endpoints in each grid cell. e PSCF method has been extensively used to identify the potential sources of air pollution [53,54]. However, the two grid cells may have the same PSCF value when sample concentrations are slightly or extremely higher than the criterion or extremely high. erefore, the PSCF method struggles to distinguish strong sources from moderate sources. e CWT method can overcome the limitation of the PSCF method because it can reflect the degree of influence of different potential sources by calculating the trajectory-weighted concentration.
e average weighted concentration in the grid cell (i, j) is calculated as where C ij is the average weighted concentration, l is the index of the trajectory, M is the total number of trajectories, and T ijl is the time spent in the grid cell (i, j) by trajectory l.

Multifractal Method.
We used coupling detrended fluctuation analysis (CDFA) to investigate the coupled correlations between several non-stationary time series. According to the CDFA method, we can perform the multifractal detrended fluctuation analysis (MF-DFA) with one time series of length N and perform the multifractal detrended cross-correlation analysis (MF-DCCA) with two time series of length N. e CDFA method can be summarized as follows. ere are several time series where n is the length of the series and t is the tth member of each time series. e CDFA method can be summarized as follows.
First, calculate the profile , where x j t represents the mean value of the kth time series.
Second, divide the profile X j (i) into N s � int(N/s) nonoverlapping segments of equal length s. Since the length N of each series is often not a multiple of the considered time scale s, a short part at the end of the profile may remain. In order not to disregard this part of the series, the same procedure is repeated starting from the opposite end. us, 2N s segments are obtained altogether. e CDFA method calculates the local trend x j v (i) for each of the 2N s segments by a least-square fit of each series. en, the detrended covariance is determined by for v � 1, 2, . . . , N s and for ird, average over all segments and obtain the qth order fluctuation function as for q ≠ 0 and for q � 0. Fourth, analyze the scaling behavior of the fluctuation function by observing the log-log plots of F x 1 ,...,x n (q, s) versus s. If the time series . , x n t is long-range power law correlated, F x 1 ,...,x n (q, s) will increase for large values of s, as a power law of where h x 1 ,...,x n (q) is the generalized scaling exponent of CDFA which can be obtained by computing the slope of the log-log plots of F x 1 ,...,x n (q, s) versus s. When h x 1 ,...,x n (q) is constant for all q, coupled correlations will be monofractal. On the contrary, when h x 1 ,...,x n (q) varies with different q, coupled correlations indicate the presence of multifractality. Specifically, if h x 1 ,...,x n (q) > 0.5, the coupled correlation is persistent (positive), which indicates that an increase (decrease) is always followed by another increase (decrease). If h x 1 ,...,x n (q) < 0.5, the coupled correlation is anti-correlated (negative), which means that an increase (decrease) is always followed by another decrease (increase). If h x 1 ,...,x n (q) � 0.5, the time series displays random walk behavior. Fifth, calculate the singularity spectrum to qualify the multifractality by using the relationship between h(q) and the Renyi exponent τ(q).
rough a Legendre transform, the multifractal spectrum f(α) is defined as where α is the Holder exponent index and f(α) is the singularity spectrum. e width of the singularity spectrum can be used to estimate the strength of multifractality as 4 Mathematical Problems in Engineering where Δα is the width of singularity spectrum. e greater the Δα value, the stronger the multifractality. In other words, a broad spectrum provides evidence for a strong multifractality nature, whereas a narrow one evidences for a weak multifractality nature.

Data Processing.
In this study, we used the daily Air Quality Index (AQI) from 1 December 2019 to 29 February 2020. Air quality was estimated through the calculation of the daily Air Quality Index which was constituted by six main pollutants (PM 2.5 , PM 10 , O 3 , SO 2 , NO 2 , and CO). e AQI and index of the six main pollutants index in Zhengzhou, Shijiazhuang, Jinan, Zibo, and Luoyang were considered in our study. ey were all obtained from the air quality publishing platform of China (https://www.resdc.cn/). Table 1 lists the total observation stations in the five cities. According to Table 1, there were 9 observation stations in Zhengzhou, 8 in Shijiazhuang, 11 in Jinan, 6 in Zibo, and 7 in Luoyang. We indicate the locations of the observation stations mentioned in a geographic map in Figure 1. Figure 2 shows the time series variation in the daily AQI in the five cities. e AQI value in Shijiazhuang is generally higher than that in the other cities. A similarity in the trends exists in the five AQI time series. e descriptive statistics of the daily AQI are listed in Table 2. e standard deviation of the AQI time series in

Mathematical Problems in Engineering
Shijiazhuang was the largest, indicating that the daily air quality had highest volatility. e skewness and kurtosis coefficient of each time series showed that each of the five AQI time series had non-normal distribution.

HYSPLIT Analysis.
We calculated the three-day (72hour) backward trajectories for every hour at an arrival altitude of 500 m in Zhengzhou city in Henan Province for December 2019, January 2020, and February 2020 using the HYSPLIT model. Zhengzhou, located in the hinterland center of the central plains in China, is heavily polluted and has high levels of particulate air pollution given the expanding population and rapidly developing economy.
We then divided the trajectories into five groups, Cluster1, Cluster2, Cluster3, Cluster4, and Cluster5, according to their origin and length using a K-means cluster analysis based on the Euclidean distance. Figures 3-5 show the five groups of trajectory clusters, respectively.
In December 2019, the clusters came from four different directions. Cluster 1 accounted for 22.31% of air masses mainly coming from the middle area of Shandong Province. Cluster 2 (19.76%) had short transport patterns, mainly coming from the middle of Henan Province. Cluster 3 (12.63%) was from northern Shanxi and passed through southern Hebei. Trajectories in Cluster 4 (28.23%) were typically long trajectories, which originated in the mid-west area of Inner Mongolia and were transported across northern Ningxia, Northern Shaanxi, and southern Shanxi. Clusters 3 and 4 are accounted for 40.86% of the air masses, indicating that air masses mostly moved from northwestern direction. Cluster 5 (17.07%) was from the middle area of Shaanxi Province and is across southern Shanxi.
We identified five air mass clusters coming from three directions in January, 2020. Cluster 1 accounted for 8.33%, which originated from Northern Shaanxi and passed through the mid-west area of Shanxi Province. Cluster 2 (16.13%) was mainly from the middle-south area of Hebei Province and transported across northern Henan. Cluster 3 (4.44%) was characterized by the longest trajectories, being transported from the middle-west area of Inner Mongolia to Zhengzhou via northern Ningxia, Northern Shaanxi, and southern Shanxi. Cluster 4 (47.58%) refers to the main pollution masses starting from northern Henan. Cluster 5 (23.52%) was mainly from northwestern part of Henan. Clusters 3 and 4 were from different directions in Henan, showing that the air pollution was mainly from the local area.
In February 2020, we identified five air mass clusters coming from four directions. Cluster 1 (17.67%) was mainly from the Bohai Rim region and eastern Shandong and passed through eastern Henan. Cluster 2 (12.64%), was from the middle area of Inner Mongolia, traveled through Northern Shaanxi and southern Shanxi, and then arrived at Zhengzhou. Cluster 3 (8.33%) was mainly from southwestern part of Hebei. Cluster 4 (17.82%) was from Northern Shaanxi and passed through southern Shanxi prior to arriving in Zhengzhou. Cluster 4 and Cluster 2 were from the same direction. Cluster 5, accounting for 43.53% of the air masses and mainly starting from the southwestern part of Henan, refers to the main pollution masses. e PSCF and CWT models were then applied to backward trajectories in order to identify potential source regions.
e results of PSCF and CWT analyses for air pollution of Zhengzhou are shown in Figure 6. e highest PSCF values were obtained from the middle-western Henan, southern Hebei, and eastern Shandong. e plausible explanation is that the air pollutants originating from other regions may have accumulated and remained in this inland area for a while in winter due to the thermal stability and slow air movement. On the other hand, China is one of the countries that use coal as the main energy, and coal plays an important role in economic and social development. e enhanced coal combustion from heating in winter increases emission rate of air pollutants during the heating period [55][56][57]. As a result, the anthropogenic emission from coal burning has a major impact on pollutant concentrations. We further used the CWT model to reveal the potential source areas. e results also showed that the major potential source areas were the middle-western Henan, southern Hebei, and eastern Shandong areas.
Since Shijiazhuang, Jinan, Zibo, and Luoyang showed higher concentrations of the air pollution than other areas in the middle-western Henan, southern Hebei, and eastern Shandong, we choose these four representative cities to explore how close the air pollution of Zhengzhou is to other areas.

Fuzzy Multifractal Detrended Fluctuation Analysis.
After identifying the source regions, we used the multifractal methods to investigate how close the air pollution of Zhengzhou was to that of each of other four cities. We first is implies that a power-law function relation exists between F(q, s) and q. Hence, the AQI time series has a multifractality nature. e generalized scaling exponent, h(q), can be obtained by calculating the slopes of ln(F(q, s)) against time scale ln(s) based on equations (4)- (8). e right part of Figure 7 shows the h(q) value when q varies from −10 to 10 for the AQI time series. All h(q) values decrease as q varies from −10 to 10, which also confirms that the AQI time series has multifractal characteristics.
e Renyi exponent τ(q) and singularity spectrum f(α) were calculated based on equations (9) and (10). e plot of τ(q) versus fluctuation order q and f(α) versus α is displayed in Figure 8. e plot of τ(q) is convex rather than a straight line. e non-linear nature proves the existence of   multifractality. e width of the multifractal spectrum, f(α), is significantly non-zero, indicating that the AQI time series are multifractal in nature.

Multifractal Detrended Cross-Correlation Analysis.
In order to investigate how close the air pollution of Zhengzhou is to that of each of other four cities, we used the MF-DCCA method in this section. We firstly investigated if the cross-correlations existed between the air pollution in Zhengzhou and that in the other four cities using a cross-correlation statistic proposed by Podobnik et al. [58] e cross-correlation statistic between two series x(i) and y(i) of equal length N was calculated as follows: en, the cross-correlation test statistic was calculated as where Q cc (m) is approximately x 2 (m) distributed with m degrees of freedom. e null hypothesis is that none of the first m cross-correlation coefficients is different from zero. If the cross-correlations do not exist between two time series, the cross-correlation test agrees well with the x 2 (m) distribution. If the cross-correlations test exceeds the critical value of the x 2 (m) distribution, the cross-correlations are significant at a certain significance level. Figure 9 shows the critical value of the x 2 (m) distribution at 5% level of significance for the degrees of freedom varying from 1 to 1000. e results showed that the cross-correlation statistic Q cc (m) of the four pairs was larger than (or close to) the critical value of the x 2 (m) at 5% level of significance. is indicated that cross-correlations existed between the AQI time series in Zhengzhou and that of each of the other four cities.
Secondly, we used the MF-DCCA method to quantitatively present the cross-correlations quantitatively between the air pollution in Zhengzhou and that in the other four cities. e range of time scale s was also set ranging from 10 events to N/4 events, where N is the length of each AQI time series in the five cities. h xy (q) values with q varying from −10 to 10 are shown in Table 3. All the cross-correlation exponents h xy (q) decrease with q varying from −10 to 10, indicating that the air pollution of Zhengzhou and that of four cities have multifractal characteristics. Moreover, the four h xy (q) values are greater than 0.5 when q is equal to two, which means that the two correlated series are persistent or show long-range dependence.
irdly, we calculated the width of the multifractal spectrum based on equation (11). e results are shown in Table 4. e multifractal spectrum provides evidence of the existence of multifractality. Moreover, its width can be used to estimate the strength of multifractality. e greater the Δα value, the closer the connection of the air pollution in Zhengzhou to that of a specific city. e Δα value was the largest in Jinan, followed by Shijiazhuang, Zibo, and Luoyang. Hence, the air pollution in Zhengzhou is closest to that in Jinan, followed by Shijiazhuang, Zibo, and Luoyang. According to the air quality publishing platform in China, Shijiazhuang and Jinan are two of the top ten most polluted cities in China, and Zibo and Luoyang are also heavy polluted cities. e multifractal results agree with the statistical results reported by the air quality publishing platform in China. We therefore concluded that Shijiazhuang and Jinan play an important role in affecting the air pollution of Zhengzhou, and Zibo and Luoyang are less important.

Coupling Detrended Fluctuation Analysis.
e MF-DCCA method was used to qualify the cross-correlations of the air pollution in Zhengzhou and that in the other four cities after identifying the potential source regions using the HYSPLIT analysis. In this part, we used the coupling detrended fluctuation analysis (CDFA) method to further  Table 5, which shows that the h x 1 ,...,x n (q) values in the four cities decreased when the q value increased. is means that the six air pollutant time series showed multifractality in each of the four cities. e two main sources of multifractality are long-range temporal correlations for small and large fluctuations and fat-tailed probability distributions in fluctuations. We shuffled and phase-randomized the original series to identify the contribution of each source. e shuffling procedure can preserve the probability distribution of the original series but destroy any temporal correlations. e shuffled series showed single fractal characteristics, that is, when the sources of multifractality were only from long-range temporal correlations, the generalized scaling exponent index h(q) was 0.5. e phase-randomization procedure can weaken the non-Gaussianity of the time series while leaving the correlation of the time series unchanged. e generalized scaling exponent index h(q) processed by phase randomization was independent of q.
In order to identify the sources of multifractalities in the four cities, we shuffled and phase-randomized the original series and obtained a shuffled series and the surrogated series, respectively. e width of the multifractal spectrum of each time series is reported in Table 6. e spectral widths of both the shuffled and surrogated series are smaller than those of the original series, indicating that long-range temporal correlations and fat-tailed probability distributions are the common causes of multifractality in the four cities. Moreover, the spectra widths of the shuffled series are smaller than those of the corresponding surrogated series in each city, which implies that long-range temporal correlations play an important role in multifractality, while fattailed probability distributions are less important.        According to the above principle, we used either the shuffling or phase-randomization procedure to modify one of the air pollutants time series to identify the sources of the multifractality in the four cities. e difference in the generalized scaling exponent value between the original time series and the modified time series is due to either the longrange temporal correlations or fat-tailed probability distributions. e larger the difference is, the larger the contribution of one time series to the coupled correlations of six air pollutant time series will be. Hence, we investigated the role of the given air pollutant in air pollution of the four cities. Figures 10-13 show generalized scaling exponent values based on the shuffling procedure.
For Shijiazhuang and Jinan, carbon monoxide (CO) had the strongest impact on the air pollution, while sulfur dioxide (SO 2 ) had the smallest impact on the air pollution. Nitrogen dioxide (NO 2 ) played an important role in the air pollution of Jinan, but less so in Shijiazhuang. e impact of inhalable particulate matter (PM 10 ) on the air pollution in Zibo was the largest, followed by fine particulate matter (PM 2.5 ) and carbon monoxide (CO). Ozone (O 3 ) had the smallest impact on the air pollution in Zibo; SO 2 played a major role in affecting the air pollution in Luoyang.
Although fat-tailed probability distributions were less important than long-range temporal correlations in multifractality, Figures 14-17 show the h x 1 ,...,x n (q) values based on the phase-randomization procedure. Carbon monoxide (CO) had the greatest impact on the air pollution, and sulfur dioxide (SO 2 ) had the smallest impact on the air pollution in Shijiazhuang and Jinan. Nitrogen dioxide (NO 2 ) more strongly affected the air pollution in Jinan than in Shijiazhuang. e result is inconsistent with the findings based on the shuffled procedure. e impact of inhalable particulate matter (PM 10 ) on the air pollution was the largest in Zibo. We could not differentiate the contributions of SO 2, PM 2.5 , and CO according to the figure in the h x 1 ,...,x n (q) value between the original time series, and each of these modified time series was approximately same. As for Luoyang city, SO 2 had the largest impact on the air pollution. We still could not conclude which air pollutants contributed more to the air pollution. e average concentrations of the air pollutants in Shijiazhuang, Jinan, Zibo, and Luoyang were higher than the national average. e average concentrations of carbon monoxide (CO) were the highest in Shijiazhuang, followed by Zibo and Jinan. e average concentrations of nitrogen dioxide (NO 2 ), inhalable particulate matter (PM 10 ), and fine particulate matter (PM 2.5 ) were generally high in the four cities. e average concentration of sulfur dioxide (SO 2 ) was the highest in Luoyang. e air pollution in these four cities could be characterized by obvious seasonal variation, that is, the air pollution is serious in winter, and less so in spring and autumn. e main reason for this finding is that these four inland cities are located in the north of China, where people burn coal for heating in winter. In winter, the weather is dry with little precipitation, light wind, and frequent inversion weather, which results in the minimal diffusion of the emissions of motor vehicles and coal difficult to diffuse. In addition to industrial emissions, vehicle emissions and adverse meteorological conditions also exacerbate air pollution.
In summary, Jinan and Shijiazhuang have important impacts on the air pollution in Zhengzhou, followed by Zibo and Luoyang. Among these source regions, carbon monoxide (CO), nitrogen dioxide (NO 2 ), and inhalable particulate matter (PM 10 ) had important influences on the air pollution. Sulfur dioxide (SO 2 ) had a small impact on the air pollution in Shijiazhuang and Jinan but a strong impact on the air pollution in Luoyang. Ozone (O 3 ) had the smallest impact on the air pollution in Zibo. e air pollutants may be transported to Zhengzhou and affect local air quality. Carbon monoxide (CO), nitrogen dioxide (NO 2 ), and inhalable particulate matter (PM 10 ) played a major role in

Discussion
In this study, we used multifractal methods combined with the backward trajectory model to study the air pollution in Zhengzhou. Compared to previous studies, we adopted several methods and investigated several air pollutants.
Using the HYSPLIT analysis, we calculated three-day backward trajectories for every hour at an arrival altitude of 500 m in Zhengzhou; subsequently, we classified the trajectories into five clusters. e results showed that the clusters of air masses mainly came from Shandong Province, Shaanxi Province, Shanxi Province, Inner Mongolia, Hebei Province, and Henan Province. Similar to the results of our present study, a previous study found that some trajectories arriving at Hefei, an inland city in Central China, were mainly from Shandong Province, in addition to Liaoning Province, Jiangsu Province, and Hefei city [59].
Comparing the results in three different months, the directions of the five air mass clusters did not change significantly. In December 2019, the clusters were from four different directions: northwest, northeast, north, and south.
ere were five air mass clusters from three directions in January 2020, which are northwest, northeast, and west. In February 2020, there were five air mass clusters from four directions of northwest, northeast, and south.
Previous studies have reported that the air pollutant concentration may increase because of anthropogenic combustion for heating purposes. e wind direction is an important factor affecting the aerosol transported from areas  where biomass is openly burned area, leading to the transboundary and local air pollution transport [26,37,60]. For this reason, the large amount of burning coal to create heat in winter causes environmental pollution. In our study, the results of PSCF and CWT analyses showed that the potential source areas are mainly middle-western Henan, southern Hebei, and eastern Shandong. e sample period was December, January, and February, which was winter in north Yangtze River area in China. e emission rate of air pollutants may be increased by the increased amount of coal combusted for heating, and air pollutants may accumulate and remain in this area due to meteorological conditions [55][56][57]. Additionally, most backward trajectories indicated that the air pollutants were from neighboring provinces. is is consistent with the results of a previous study [54]. e log-log plots of F(q, s) versus s and the generalized Hurst exponent index were calculated with the MF-DFA method.
e results showed that the AQI time series of Zhengzhou has multifractal characteristics. e Renyi exponent index and singularity spectrum confirmed that the results are stable. In the previous studies, persistent powerlaw correlations in PM 10 and O 3 fluctuations were detected in Baltimore, USA, which reflected short-and long-term meteorological features related to the self-organized critical behavior of the atmosphere [61]. Diosdado et al. [62] found that air pollutants' time series, inhalable particulate matter (PM 10 ), sulfur dioxide (SO 2 ), nitrogen dioxide (NO 2 ), ozone (O 3 ), and carbon monoxide (CO) were multifractal in Mexico. Dong et al. [50] verified the existence of a significant power-law behavior and long-range correlations for the PM 2.5 and PM 10 series in Shanghai. Shen et al. [49] found that the daily Air Pollution Index (API) in Nanjing was multifractal. is indicates that the findings of our study are consistent with those of previous studies, proving important complementary information to previous studies.
We used the multifractal methods and the backward trajectory model to explore the cross-correlations between the air pollution in Zhengzhou and that in four other cities. Our multifractal study considered meteorological conditions, similar to a previous study [51]. ey also investigated the multifractal characteristics of PM 2.5 concentration combined with meteorological conditions in Beijing and Hong Kong. In contrast, they used specific meteorological factors such as temperature, air pressure, relative humidity, and wind speed to explore the cross-correlations between PM 2.5 concentration and each of meteorological factors. Stan et al. [52] also presented the results of the cross-correlations between temperature and each of the several gases. ey found that all correlated series had multifractal properties and persistent long-range correlations. Additionally, they conducted a fuzzy multifractal analysis of ozone and some precursors (nitrogen oxides (NOXs), non-methane hydrocarbons (NMHCs), total hydrocarbons (THCs), and methane (CH 4 )) during the cold season using a generalization of multifractal detrended crosscorrelation analysis (MF-DCCA). e results showed that the cross-correlations for different pairs of ozone precursors had multifractal characteristics and the generalized Hurst exponents for different pairs differed. e MF-DCCA is a powerful method that can be used to investigate the cross-correlations between two time series. However, it cannot be used to investigate simultaneously the coupled correlations of several time series simultaneously. For example, the MF-DCCA can be used to analyze the cross-correlation of all air pollutant time series with PM 2.5 and PM 10 . However, it cannot investigate the coupled correlations of nine air pollution time series such as fine particulate matter (PM 2.5 ), inhalable particulate matter (PM 10 ), sulfur dioxide (SO 2 ), nitrogen dioxide (NO 2 ), ozone (O 3 ), carbon monoxide (CO), nitric oxide (NO), nitrogen oxide (NOx), and ammonia(NH 3 ) at the same time [63]. e coupling detrended fluctuation analysis (CDFA) is an extended MF-DXA that is used to investigate the coupled correlations of several time series simultaneously [45].
In this study, we used this method to explore what is the role of a given air pollutant in the air pollution in four cities. After the shuffling procedure and phase-randomization, we found that carbon monoxide (CO), nitrogen dioxide (NO 2 ), inhalable particulate matter (PM 10 ), and sulfur dioxide (SO 2 ) were important air pollutants in affecting air quality, whereas fine particulate matter (PM 2.5 ) and ozone (O 3 ) were less important. Our findings are different from those reported by Hedayatifar et al. [45]. Using the shuffling procedure, they found that ozone (O 3 ) played an important role in affecting air quality, whereas nitrogen dioxide (NO 2 ) was less important. e importance of nitrogen dioxide (NO 2 ) in affecting air quality increases when the original series is phase-randomized.

Conclusions
In this study, we used the backward trajectory model and multifractal methods to analyze the air pollution in Zhengzhou.
e HYSPLIT model was used to divide the trajectories of air masses into five groups: Cluster 1, Cluster 2, Cluster 3, Cluster 4, and Cluster 5. Most clusters in each month originated from southern Hebei, eastern Shandong, and mid-western Henan and were then transported to Zhengzhou. We further adopted the PSCF and CWT models to identify potential source regions. e results are consistent with the HYSPLIT analysis. Four representative cities were selected to explore how close the air pollution of Zhengzhou is to other areas on the basis of air polluted concentration.
We used the MF-DFA method to explore if the air polluted time series in Zhengzhou had multifractal characteristics. Log-log plots of fluctuation functions, generalized Hurst exponent, Renyi exponent, and singularity spectrum all provided evidence of the multifractality of the AQI time series in Zhengzhou. Hence, we found that the multifractal method was feasible for our purpose. We then used the MF-DCCA method to investigate how close the AQI series of Zhengzhou is to that of each of other four cities. All of the cross-correlation exponents h xy (q) decreased with q varying from −10 to 10, which implied that the AQI series of Zhengzhou and those of the other four cities had multifractal nature. e widths of the multifractal spectra showed that the air pollution in Zhengzhou was closest to that in Jinan, followed by Shijiazhuang, Zibo, and Luoyang. Finally, we used the CDFA method to further explore what is the role of the given air pollutants, PM 2.5 ,PM 10 ,SO 2 ,NO 2 ,O 3 , and CO, in the air pollution of the four cities. e results showed that carbon monoxide (CO) had the strongest impact on the air pollution, while sulfur dioxide (SO 2 ) had the least impact on the air pollution in Shijiazhuang and Jinan. Nitrogen dioxide (NO 2 ) was important to the air pollution in Jinan but less so in Shijiazhuang. As for Zibo, inhalable particulate matter (PM 10 ) had the strongest impact on the air pollution, while ozone (O 3 ) had the weakest impact on the air pollution. SO 2 affected the air pollution in Luoyang the most. e findings of this research provide useful references for air pollution sources and their potential contributions in Zhengzhou. Each city should take measures to control industrial enterprises' pollution emissions, to manage motor vehicle pollution emissions, and to establish no-coal-fired areas and basically no-coal-fired areas. Enterprises should improve their technological level, reform the production process to realize ecological benefits, and optimize their energy structure. Importantly, air pollution control and management is not one city's independent responsibility; instead, many cities must coordinate to achieve this goal. Hence, cities should collaborate to improve the ecological environment and achieve sustainable development.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.