A DFA-based bivariate regression model for estimating the dependence of PM2.5 among neighbouring cities

by   Fang Wang, et al.

On the basis of detrended fluctuation analysis (DFA), we propose a new bivariate linear regression model. This new model provides estimators of multi-scale regression coefficients to measure the dependence between variables and corresponding variables of interest with multi-scales. Numerical tests are performed to illustrate that the proposed DFA-based regression estimators are capable of accurately depicting the dependence between the variables of interest and can be used to identify different dependence at different time scales. We apply this model to analyze the PM2.5 series of three adjacent cities (Beijing, Tianjin, and Baoding) in Northern China. The estimated regression coefficients confirmed the dependence of PM2.5 among the three cities and illustrated that each city has different influence on the others at different seasons and at different time scales. Two statistics based on the scale-dependent t-statistic and the partial detrended cross-correlation coefficient are used to demonstrate the significance of the dependence. Three new scale-dependent evaluation indices show that the new DFA-based bivariate regression model can provide rich information on studied variables.




Coefficient Decomposition of Spatial Regressive Models Based on Standardized Variables

Spatial autocorrelation analysis is the basis for spatial autoregressive...

A Stratification Approach to Partial Dependence for Codependent Variables

Model interpretability is important to machine learning practitioners, a...

Max-sum tests for cross-sectional dependence of high-demensional panel data

We consider a testing problem for cross-sectional dependence for high-di...

Explanation of multicollinearity using the decomposition theorem of ordinary linear regression models

In a multiple linear regression model, the algebraic formula of the deco...

Are Instrumental Variables Really That Instrumental? Endogeneity Resolution in Regression Models for Comparative Studies

We provide a justification for why, and when, endogeneity will not cause...

Estimating The Proportion of Signal Variables Under Arbitrary Covariance Dependence

Estimating the proportion of signals hidden in a large amount of noise v...

Exploring and measuring non-linear correlations: Copulas, Lightspeed Transportation and Clustering

We propose a methodology to explore and measure the pairwise correlation...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

In recent years, air pollution has become a more and more serious problem around the world. The new air quality model presented by the World Health Organization in 2016 confirmed that 92% of the world’s population lives in areas where air quality levels exceed their limits [1]. Fortunately, more and more governments have realized the importance of managing air pollution and some actions have been placed. Nowadays, a common topic around the world is the governance of the air pollution source such as smog (the main ingredient is fine particulate matter). Many researchers have been involved in the study on the cause and propagation of smog [2, 3, 4, 5, 6, 7, 8, 9]. Modern statistic methods provide some new perspectives to assess smog trends and propagation characteristics. Among them, most studies have focused on studying the correlations among various air pollution indicators including air pollution index (API), air quality index (AQI), fine particulate matter of PM2.5 (diameter ) concentrations, and PM10 (diameter ) concentrations, and very limited studies considered the correlations among neighboring areas. A common sense is that smog produced at one source place can spread to surrounding areas [6, 7, 8, 9, 10]. Therefore, it is more practical to explore the dependence of air pollution indicators among adjacent cities as it helps assess the causes of local smog and its spread behavior. It has been found by a newly proposed time-lagged cross-correlation coefficient in Ref. [10] that there are different degrees of correlation for PM2.5 series between four neighboring cities in Northern China. However, what has not been investigated is how the PM2.5 series of one city depends on those of the neighbouring cities. In this work, we will develop a detrended fluctuation analysis (DFA)-based bivariate regression model to investigate this dependence.

The simplest and maturest method to describe the dependence of variables is the linear regression. However, the information gained from the traditional linear regression cannot fully meet our need of investigation on the dependence among different variables at different time periods. On the other hand, note that the DFA proposed in 1990s [11, 12] performs excellently in analyzing the long-range correlations [13] of a nonstationary series with fractality and multifractality [14, 15] at different time-scales. To obtain the cross-correlation between two nonstationary series, DFA was extended to the detrended cross-correlation analysis (DCCA) [16]. By defining scale-dependent detrended fluctuation functions, the methods of DFA and DCCA together with their extensions have been applied in a wide range of disciplines [17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]

. Since the ordinary least squares (OLS) method expresses the estimated parameters of standard regression framework as a form of variances and covariances, it builds a bridge between the regression framework and the family of DFA/DCCA as the latter can also produce variances and covariances. Then, the idea of estimating multiple time scale regression coefficients can be achieved by the DFA/DCCA. Recently, Kristoufek 

[31] constructed a simple DFA-based regression framework exactly by this bridge. The selected examples show the relationship between the pair of variables varies strongly across scales.

In this work, we focus on the interaction of PM2.5 series of three adjacent cities in Northern China, namely, Beijing, Tianjin, and Baoding. The three cities form a triangle shape in the map. The distances between Beijing and Tianjin, Beijing and Baoding, and Tianjin and Baoding are about 115km, 140km, and 150km, respectively. All three cities have a population of more than 10 million and have been greatly affected by heavy smog in recent years. The real-time data of PM 2.5 series of these three cities from December 1, 2013 to November 30, 2016 are chosen for our consideration, which are taken from the Ministry of Environmental Protection of the People’s Republic of China (http://datacenter.mep.gov.cn

). The original data show an obvious periodic characteristic and roughly similar trends among the three cities, which imply that there is a possible relevance between per two cities of them. To verify that, the partial correlation technique is employed to get the intrinsic relations between two cities by deleting the interference from the third variable. Four seasons, classified as winter (December, January, and February), spring (March, April, and May), summer (June, July, and August), and fall (September, October, and November), are considered. The results are listed in Table 


winter spring summer fall
Beijing vs. Tianjin 0.3048 0.3072 0.2625 0.1660
Beijing vs. Baoding 0.2745 0.4545 0.4815 0.4468
Tianjin vs. Baoding 0.5570 0.4517 0.3461 0.5992

Note: * indicates statistical significance with 0.01 significance level.

Table 1: Partial correlation coefficients and -statistics between per two cities of Beijing, Tianjin, and Baoding in four seasons.

In Table 1, we also list the -statistics (, where denotes the partial correlation coefficient between the first and second variables eliminating the effects of the third one,

is the degree of freedom) of the partial correlation coefficients to access the statistical significance at the given significance level. Unsurprisingly, Table 

1 shows that the correlations of PM2.5 between both per two cities are of statistical significance. It explains that the air quality in one city of Northern China cannot be irrelevant to that of its neighbouring cities, which implies potential dependence among the three cities. However, we also note in Table 1 that the degree of relevance is different among different cities and in different seasons though all of them are significant.

To fully detect and quantify the dependence among the PM2.5 series of the above-mentioned three cities, in this work, we construct a new bivariate regression framework which prevails the DFA method and allows us to investigate the dependence of three nonstationary series with multiple time scales. With the DFA-based variance instead of the standard variance, this new DFA bivariate regression model provides more information on the dependence among variables at different time scales. We organize the rest of this paper as follows. The performance of the proposed DFA regression model and the results on the application to PM2.5 series analysis are reported and discussed in the following section, which is followed by our conclusions. The methodologies including the standard regression method, the DFA method, and the DFA-based regression method are introduced at the end of this paper.

2 Results and Discussions

Performance of DFA estimators

The bivariate DFA-based regression model produces two time scale-based regression coefficients. This allows us to detect the dependence of a response variable and two independent variables at different time scales. In order to examine the validity of the model and show its advantages, in this section, we perform two numerical tests on the non-stationary bivariate regression frameworks

In the first test, we investigate the performance of the DFA estimators under different levels of long-term dependence in , , and . According to  [31], the setting I is given as below: two artificial series and with length are generated by ARFIMA() process with identical fractional integration parameter () and independent Gaussian noises ( and 2 ) as . The quantity is defined by , where is the Gamma function. The error-term is set as a standard Gaussian noise so that the response variable has the same parameter as the two independent variables. The regression coefficients are set as and

. Fig. 1 shows mean values and standard deviation of the two DFA estimators

(1 and 2) for the generated series with ranging from to (at the step size of ). The estimators are averaged over scales between and with a logarithmic isometric step. Each case is run times to eliminate the noise interference. It is clear that the two estimators locate the two given regression coefficients of (Fig. 1a) and 2 (Fig. 1b) unbiasedly, and are independent of the value of . In addition, the standard deviations of both estimators decrease with the increasing memory. The good performance shows that the method is feasible. On the other hand, to investigate the performance of the DFA estimators faced with a long-range dependent error-term , we use setting II given as: the memory parameter is fixed at 0.4 for both and , and the is produced by an ARFIMA process with varying from to . Other settings are as those in setting I. Fig. 2 records similar information as that in Fig. 1. Although the fluctuation of DFA estimators increases with , which is expected due to an increasing weight of the error-term in the dynamics of with the increasing memory of the error-term, we are satisfied to find that the two estimators are still unbiased pointing to the given values with a narrow range for each level of memory of the error-terms.

Our second numerical test aims to show that the DFA estimators are able to identify the dependence of studied variables at different time scales whereas the classical method cannot. To this end, a binomial multifractal series (BMFs) is employed to be regarded as the independent variable , which is constructed as , , , , , where the parameter (We take in our test), denotes the number of digit in the binary representation of the index . The variable is a Gauss variable with mean and standard deviation. Both and are of length . The bivariate regression framework is set with the same coefficients as the first test ( and ). The error-term is the Gauss noise of the same strength as . For the BMFs , we remove all values smaller than

so that only a few of the largest elements are left. In their places, we substitute Gaussian distributed random numbers with

mean and standard deviation. Then we obtain a binomial cascade series embedded in random noise. We analyze the dependence between the response variable and two independent variables and find that the estimated is unbiased at with a desirable error bar for every time scale, as shown in Fig. 3. However, the performance of has changed a lot. The dependence between and is obviously less than the given value at the smaller scales contrary to the larger ones. This is because in the smaller scales, the dependency has been destroyed by the random noise. Our DFA estimators have the capability to recognize this effect while the classical estimators fail to do so (see the errorbar with circle symbol in Fig. 3).

Performance of the three models’ regression coefficients

As mentioned above, air pollution in Northern China is very serious in recent years. Fine particulate matter from industrial exhaust and smoke dust forms smog to fill in the air. We now apply our DFA regression model to investigate the dependence of PM2.5 series in these three cities. We build three bivariate models for Beijing, Tianjin, and Baoding, respectively. In Model I, the dependent variable () is the PM2.5 of Beijing while the two independent variables are the PM2.5 of Tianjin () and Baoding (); in Model II, is the PM2.5 of Tianjin, is the PM2.5 series of Beijing and is the PM2.5 series of Baoding; in Model III, is the PM2.5 of Baoding, and stand for the PM2.5 series in Beijing and Tianjin, respectively. In this section, we first show the performance of the regression coefficients at different scales in the three models and then make two statistical tests for the two regression coefficients in each model. Some evaluations for the DFA-based regression and the standard regression are conducted at the end of this section.

The two regression coefficient estimators together with their standard deviations of the three models are sketched in Figs. 4–6, respectively. As expected, the effect is obviously positive. However, a strong variation across scales is found in different seasons. More specifically,

  • In the Beijing’s model, Tianjin () has strongly positive effect in every season, especially for the larger time scales. On the contrary, Baoding () has different effects on Beijing. Compared to spring and summer, the effect is quite weak in the other two seasons, especially in winter, is nearly when the scale is more than hours.

  • In the Tianjin’s model, Baoding () presents more unstable effect at different scales. Particularly in summer, is close to from the smaller scale to the larger scale at about days ( hours), which implies that the positive correlation between Tianjin and Baoding can last less than days. In addition, the two coefficients are less than in most days, which indicates that Beijing and Baoding have little impact on the PM2.5 in Tianjin.

  • For the model of Baoding, the effect of Tianjin () to Baoding is similar to that of Baoding to Tianjin in model II. However, the fact that after approximately days ( hours) the effect reaches the value greater than indicates that an increase of unit PM2.5 concentration of Tianjin will lead to the increase of more than

    unit PM2.5 concentration in Baoding. In this regard, Tianjin has more impact on Baoding. In addition, the narrow confidence intervals and low standard deviations (less than

    ) shown in all sub-plots suggest satisfied reliability of the estimates.

Statistic significance tests of regression coefficients

As mentioned above, the estimated is able to theoretically describe the dependence between the impulse variables and the response variables at different time scales. In theory, as long as is not equal to zero, the independent variable will affect . However, for finite time series, is not always equal to even in the absence of relationship between and due to the size limitation. Therefore, we perform a hypothesis test for the estimated

to ensure the significance. The standard regression analysis provides a so-called

statistic defined as (, 2) for this purpose. We have for the bivariate regression model as . In general, if with a given

, we should reject the null hypothesis of

and the dependence between and is considered to be statistically significant. However, since lots of time scales are taken accounted in the DFA regression model, using a single critical value of is inappropriate. A correct way is to generate a critical value for each time scale. To this end, inspired by the idea proposed by Podobnik et al. [32], we shuffle the considered PM2.5 series and repeat the DFA regression calculations for

times. Then let the integral of probability distribution function (PDF) from

to be equal to (here, we take ). As an example, we show the PDF of with five given ’s produced by the shuffled PM2.5 series of fall in Fig. 7.

As expected, the symmetrical PDF of

converges to a Gaussian distribution according to the central limit theorem. In addition, the critical value increases as

increases. This implies that large time scale may strengthen dependence between two variables. By using , we can determine whether the dependence between the impulse variable and the response variable is significant or not. In practice, the dependence between and is present when is larger than . For the four seasons, the scale-dependent -statistics of regression coefficient together with the scale-dependent critical value are presented in Fig. 8.

Note that in Model I (for Beijing), the -statistics of (Tianjin’s coefficient) is equal to that of (Beijing’s coefficient) in Model II (for Tianjin), the -statistics of (Baoding’s coefficient) is equal to that of (Beijing’s coefficient) in Model III (for Baoding), and in Model II, the -statistics of Baoding’s coefficient is equal to that of Tianjin’s coefficient in Model III (for Baoding). Here the three colored lines with different symbols represent the -statistics between each per two cities while the black dashed line stands for . The partial DCCA coefficient is recently developed to uncover the intrinsic relation for two nonstationary series at different time scales. We also calculate the partial DCCA coefficients of Beijing and Tianjin, Beijing and Baoding, and Tianjin and Baoding, respectively, and present the results in Fig. 9. For the same purpose of testing the statistical significance, we also produce a critical value for the four seasons. Similarly, the PM2.5 data are shuffled times in the PDCCA calculations repeatedly, and thus for 99% confidence level is obtained, which is also shown in Fig. 9.

Comparing results in Fig. 8 and Fig. 9 gives amazing similarities, which are also in agreement with the results shown in Figs. 4–6. Based on the results, we can draw the following three main points.

  • The dependence between Beijing and Tianjin (the blue square line) gradually increases with the increasing time scales in all seasons. However, the dependence between the two cities is lower than other cities. This finding uncovers that the reason for the serious air pollution in these two cities are mainly due to their own heavy smog or are impacted by other cities.

  • The dependence between Beijing and Baoding (the green triangle line) is significant in spring, summer, and fall. In winter, however, the dependence disappears at long time scale, which implies that the two cities can only affect each other at relatively short term. Moreover, compared to winter and fall, the dependence is much stronger in spring and summer, especially at long time scales, which indicates that they affect much longer in warm weather.

  • In spring and summer, the -statistics and of Tianjin vs. Baoding (the red circle line) go down through the critical lines of and , respectively at about hours. This suggests that the dependence between Tianjin and Baoding will disappear when it’s more than one month. However, the exact opposite occurs in winter and fall. In these two seasons, both -statistic and increase with the increasing time scales, which demonstrates that the interaction of bad air quality between the two cities will last longer in cold days.

Evaluations of DFA-based regression model

To evaluate our estimated DFA-based bivariate regression model, we plot the scale-dependent determination coefficient , and the beta coefficient and the average elasticity coefficient in Fig. 10, and Figs. 11–13, respectively.

To show the new model provides more information than the standard regression model does, we also include the three corresponding coefficients of standard bivariate regression model in these figures. As seen from Fig. 10 that is superior to the standard at most time scales. The good performance illustrates that one will gain richer information in explaining the response variable when using our DFA-based regression model. On the other hand, we can conclude from Figs. 11–13 that (1) Baoding has more influence than Tianjin on Beijing in all seasons except for winter. (2) Tianjin is more sensitive to Baoding’s changes in air quality than Beijing’s in winter and fall. (3) Tianjin affects Baoding more than Beijing does in winter and fall, but less in the other two seasons. In addition, Figs. 11–13 illustrate that the standard and can be seen as the mean values of the DFA-based and , respectively. This means that and are able to measure the dependence degree of the studied independent variable on the dependent variable in all directions. Thus one can access the measurement according to his/her needs. For example, in winter of Model I, we find that the and are larger than and , respectively, at smaller scales but much smaller at larger scales, which shows that the sensitivity of to (Baoding) is greater than that of to (Tianjin) for short term ( hours) but Tianjin is more sensitive to Beijing at the long term. This can help air quality inspectors make the correct analysis for Beijing’s PM2.5 at different periods.

3 Conclusions

The study of dependence between variables helps expose the causal relationship and correlation of the variables of interest in the real world. The linear regression model is undoubtedly one of the simplest methods among many approaches. However, single variety of regression coefficient and evaluation index cannot show all aspects of the dependence between independent variables and dependent variable. As a meaningful extension, we design a new framework for bivariate regression model using the prevailing DFA method. The proposed bivariate DFA regression model allows us to estimate multi-scale regression coefficients and other corresponding scale-dependent evaluation indicators. It has been shown via two artificial tests that these DFA-based regression coefficients are able to describe the dependence between the response variable and two independent variables exactly; and can capture different dependence at different time scales.

An application of the new model to the study of dependence of PM2.5 series among three heavily air polluted cities in Northern China unveils that huge difference of the dependence exists in per two cities in different seasons and at different periods. Three new indicators of the scale-dependent determination coefficient, the scale-dependent beta coefficient, and the scale-dependent elasticity coefficient are proposed, which turned out to be more practical than those in standard regression models. Three main points can be concluded as (1) Beijing and Baoding have little impact on the PM2.5 in Tianjin while Tianjin takes more impact on Baoding and the air quality of Beijing is more sensitive to the changes in Baoding. (2) In contrast, the air quality in Beijing and Tianjin is not significantly relevant, while the air quality in Tianjin and Baoding has a very significant impact on each other especially in the cold weather. (3) In comparison, the fluctuation of PM2.5 in Baoding has the greatest impact on the other two cities in most days. While Baoding’s air quality is more sensitive to Beijing’s changes in spring and summer, and is more sensitive to Tianjin’s changes in winter and fall. These findings may provide some useful insights on understanding air pollution sources among cities in Northern China.

4 Methods

The standard bivariate regression model

To study the dependence of air quality among three neighboring cities, we consider a bivariate linear regression model as


where is a dependent variable, and are two independent variables, is a Gaussian error term with zero mean value, and (, ) is the partial regression coefficient characterizing the dependence on . The most critical work in empirical studies is to estimate and . The OLS method gives


where denotes the mean value of the whole time period, , , and . Then the estimator of residuals can be determined by . With it one can obtain the estimators of variance of the two regression coefficients as

The variance illustrates the accuracy of the estimated parameters. The estimated regression coefficients together with their variances can be further employed for hypothesis test and model evaluation. As an important indicator to evaluate the regression model, the determination coefficient is defined by


with the range of . measures a proportion of variance of explained by and and higher value of implies better model interpretation ability. Moreover, to quantify sensitivity of explained variable to each explaining variable, two quantities, namely, the beta coefficient (denoted as ) and the average elasticity coefficient (denoted as ), are defined




which can explain the relative importance of variables and to . According to [31], the advantage of translating the standard notation into variance and covariance shown on the right-hand side of Eqs. (4.2)-(4.5) is available to use the DFA/DCCA methods based on the same idea.

The DFA-based variance and DCCA-based covariance functions

DFA and DCCA methods are described as follows. For a time series , , 2, , , we split its profile into nonoverlapping segments of equal length , denoted as , , 2, , . The same procedure is repeated starting from the opposite end to avoid disregarding a short part of the series in the end and thus segments are obtained altogether. In the segment, we have for , 2, , and for , , , , where , 2, , . In each segment, the local linear (or other) trend  [33, 34] can be fitted as (in our work, we use order polynomial to fit the trend in each segment). Fluctuation function is then defined for each segment as


Averaging the fluctuation over all segments yields


which is the so-called DFA-based scale-dependent variance function. To obtain the scale-dependent covariance of two equal length series and , , 2, , , we only need to translate the univariate fluctuation function in each segment and average fluctuation into the bivariate case directly,


The scale-characteristic fluctuation is the so-called DCCA-based covariance, which expresses the cross-correlation fluctuations between the series of and . Thus we have obtained all objects to create the DFA-based regression model. But for purpose of testing, we need some accessories of the DFA process. The DCCA cross-correlation coefficient , proposed by Zebende [35], can measure the cross-correlation between two nonstationary series at multiple time scales, which is defined as


To access intrinsic relations between the two time series on time scales of , Yuan et al. [36] and Qian et al. [37] developed a so-called partial DCCA coefficient independently, which applies partial correlation technique to delete the impact of other variables on the two currently studied variables. This coefficient is defined as


where is the inverse matrix of the cross-correlation matrix produced by of , , , and subscripts and stand respectively for the row and column of the location of .

The DFA-based bivariate regression model

We now translate the standard bivariate regression process described above into the DFA-based bivariate regression model. The two estimators in Eq. (4.2) can be extended to the scale-dependent estimators in the following way using the scale-dependent variance and covariance defined in Eqs. (4.7) and (4.9),


Similarly, the scale-dependent residuals are

with zero mean value. Inserting the calculated into the DFA process, we obtain the fluctuation to estimate the variances of and via Eq. (4.12) as


Then Eqs. (4.3)-(4.5) can be translated into the DFA regression form as




Comparing to the standard , , and , the scale-dependent , , and express more abundant information on model interpretation from multiple time scales.


  • [1] WHO releases country estimates on air pollution exposure and health impact. http://www.who.int/mediacentre/news/releases/2016/air-pollution-estimates/en/. Sep. 27, 2016.
  • [2] Wang, S. & Hao, J. Air quality management in china: Issues, challenges, and options. Journal of Environmental Sciences 24, 2–13 (2012).
  • [3] Han, L. J., Zhou, W. Q., Li, W. F. & Li, L. Impact of urbanization level on urban air quality: A case of fine particles (pm2.5) in chinese cities. Environmental Pollution 194, 163–170 (2013).
  • [4] Han, L. J., Zhou, W. Q. & Li, W. F. City as a major source area of fine particulate (pm 2.5) in china. Environmental Pollution 206, 183–187 (2015).
  • [5] Han, L. J., Zhou, W. Q. & Li, W. F. Increasing impact of urban fine particles (pm2.5) on areas surrounding chinese cities. Scientific Reports 5, 12467 (2015).
  • [6] Shen, C. H. & Li, C. An analysis of the intrinsic cross-correlations between api and meteorological elements using dpcca. Physica A 446, 100–109 (2016).
  • [7] Shi, K. Detrended cross-correlation analysis of temperature, rainfall, pm 10 and ambient dioxins in hong kong. Atmospheric Environment 97, 130–135 (2014).
  • [8] Shen, C. H. A new detrended semipartial cross-correlation analysis: Assessing the important meteorological factors affecting api. Physics Letters A 379, 2962–2969 (2015).
  • [9] Zeng, M., Zhang, X. N. & Li, J. H. Dcca cross-correlation analysis of 3d wind field signals in indoor and outdoor environments. Intelligent Control and Automation (WCICA) 2791–2796 (2016).
  • [10] Wang, F., Wang, L. & Chen, Y. M. Detecting PM2.5’s Correlations between Neighboring Cities Using a Time-Lagged Cross-Correlation Coefficient. Scientific Reports 7, 10109 (2017).
  • [11] Peng, C. K., Buldyrev, S. V., Goldberger, A. L., Havlin, S., Simon, M. & Stanley, H. E. Finite-size effects on long-range correlations: Implications for analyzing DNA sequences. Physical Review E 47, 3730 (1993).
  • [12] Peng, C. K., Buldyrev, S. V. , Havlin, S., Simons, M., Stanley, H. E. & Goldberger, A. L. Mosaic organization of DNA nucleotides. Phys. Rev. E 49, 1685-1689 (1994).
  • [13] Kantelhardt, J. W., Koscielny-Bunde, E., Rego, H. H. A., Havlin, S.& Bunde, A. Detecting long-range correlations with detrended fluctuation analysis. Physica A 295, 441–454 (2001).
  • [14] Kantelhardt J. W., Fractal and multifractal time series[M]//Mathematics of complexity and dynamical systems, Springer, New York (2012).
  • [15] Kantelhardt, J. W. et al. Multifractal detrended fluctuation analysis of nonstationary time series. Physica A 316, 87-114 (2002).
  • [16] Podobnik, B. & Stanley, H. E. Detrended cross-correlation analysis: a new method for analyzing two nonstationary time series. Phys. Rev. Lett 100, 084102 (2008).
  • [17] Podobnik, B., Horvatic, D., Petersen, A. M. & Stanley, H. E. Cross-Correlations between Volume Change and Price Change. Proc. Natl. Acad. Sci. USA 106, 22079-22084 (2009).
  • [18] Zhou, W. X. Multifractal detrended cross-correlation analysis for two nonstationary signals. Phys. Rev. E 77, 066211 (2008).
  • [19] Wang, F., Liao, G. P., Zhou, X.  Y. & Shi, W. Multifractal detrended cross-correlation analysis for power markets. Nonlinear Dynamics 72, 353-363 (2013).
  • [20] Wang, F., Liao, G. P., Li, J. H., Li X. C. & Zhou T. J. Multifractal detrended fluctuation analysis for clustering structures of electricity price periods. Physica A 392, 5723-5734 (2013).
  • [21] Wang, F., Fan, Q. J. & Stanley H. E. Multiscale multifractal detrended-fluctuation analysis of two- dimensional surfaces. Physical Review E 93, 042213 (2016).
  • [22] Wang, F. A novel coefficient for detecting and quantifying asymmetry of California electricity market based on asymmetric detrended cross-correlation analysis. Chaos 26, 063109 (2016).
  • [23] Wei, Y. L, Yu, Z. G., Zou, H. L. & Anh V. Multifractal temporally weighted detrended cross-correlation analysis to quantify power-law cross-correlation and its application to stock markets. Chaos 27, 063111 (2017).
  • [24] Oświȩcimka, P., Drożdż, S., Forczek, M., Jadach, S. & Kwapień, J. Detrended cross-correlation analysis consistently extended to multifractality. Phys. Rev. E 89, 023305 (2014).
  • [25] Jiang, Z. Q. & Zhou, W. X. Multifractal detrending moving-average cross-correlation analysis. Phys. Rev. E 84, 016106 (2011).
  • [26] Lin, A. J., Shang, P. J. & Zhao, X. J. The cross-correlations of stock markets based on dcca and time-delay dcca. Nonlinear Dyn 67, 425-435 (2012).
  • [27] Kristoufek, L. Measuring correlations between non-stationary series with dcca coefficient. Physica A 402, 291–298 (2014).
  • [28] Kristoufek, L. Finite sample properties of power-law cross-correlations estimators. Physica A 419, 513-525 (2015).
  • [29] Yu, Z. G, Leung, Y., Chen, Y. Q.,Zhang, Q., Anh, V. & Zhou, Y. Multifractal analyses of daily rainfall in the Pearl River basin of China. Physica A 405, 193-202 (2014).
  • [30] Li, Z. W & Zhang, Y. K, Quantifying fractal dynamics of groundwater systems with detrended fluctuation analysis. Journal of hydrology 336.1-2, 139-146 (2007).
  • [31] Kristoufek, L. Detrended fluctuation analysis as a regression framework: Estimating dependence at different scales. Phys. Rev. E 91, 022802 (2015).
  • [32] Podobnik, B., Jiang, Z. Q., Zhou, W. X. & Stanley, H. E. Statistical tests for power-law cross-correlated processes. Phys. Rev. E 84, 066118 (2011).
  • [33] Oświȩcimka, P., Drożdż, S., Kwapień, J. & Górski, A. Z. Effect of detrending on multifractal characteristics. Acta Physica Polonica A 123, 597-603 (2012).
  • [34] Ludescher, J., Bogachev, M. I., Kantelhardt, J. W., Schumann, A. Y. & Bunde, A. On spurious and corrupted multifractality: The effects of additive noise, short-term memory and periodic trends. Physica A 390(13), 2480-2490 (2011).
  • [35] Zebende, G. Dcca cross-correlation coefficient, quantifying level of cross-correlation. Physica A 390, 614-618 (2011).
  • [36] Yuan, N. M., Fu, Z. T., Zhang, H.,Piao, L.,Xoplaki, E.,Luterbacher, J. Detrended partial-cross-correlation analysis: a new method for analyzing correlations in complex system. Scientific Reports 5, 8143 (2015).
  • [37] Qian, X. Y., Liu, Y. M., Jiang, Z. Q., Podobnik, B., Zhou, W. X. & Stanley, H. E. Detrended partial cross-correlation analysis of two nonstationary time series influenced by common external forces. Phys. Rev. E 91, 062816 (2015).


This work was partially supported by National Natural Science Foundation of China (No.31501227, No. 11401577), and NSERC of Canada. The authors would like to thank two anonymous reviewers and the handling editor for their comments and suggestions, which led to a great improvement to the presentation of this work.

Author Contributions

F.W. designed the framework and performed the statistical analysis. L.W. wrote the manuscript and checked the analysis. Y.C. reviewed the analysis.

Additional Information

Competing Interests: The authors declare no competing interests.

Figure 1: Bivariate DFA regression of Beijing. Main planes of subplots (a), (b), (c), and (d) show estimated DFA regression coefficients and of winter, spring, summer, and fall, respectively. Gray zones denote 95% confidence intervals. Inserts are standard deviations of and . Subscripts 1 and 2 denote Tianjin and Baoding, respectively.
Figure 2: Bivariate DFA regression of Tianjin with the same legend as in Fig. 1. Here, subscripts 1 and 2 denote Beijing and Baoding, respectively.
Figure 3: Bivariate DFA regression of Baoding with the same legend as in Fig. 1. Here, subscripts 1 and 2 denote Beijing and Tianjin, respectively.
Figure 4: PDF of critical points - statistics critical values at different scales for the statistical test with times of the shuffled PM series of fall.
Figure 5: -statistical test of the estimated DFA-based bivariate regression coefficients. (a)-(d) are for winter, spring, summer and fall, respectively. The dashed line represents the with significant levels. Above this line means decline of the null hypothesis .
Figure 6: Statistical test of DPCCA coefficients among the three cities. (a)-(d) are for winter, spring, summer, and fall, respectively. The dashed line represents the critical value of which is obtained from times Monte-Carlo simulations with 99% confidence level. Below this line suggests no cross-correlated significance.
Figure 7: Determination coefficients of bivariate DFA and standard regression model. (a)-(d) are for models of Beijing, (e)-(h) are for models of Tianjin, and (i)-(l) are for models of Baoding. The solid line denotes and the dashed line denotes .
Figure 8: Beta coefficients and elasticity coefficients of bivariate DFA and standard regression model of Beijing. The four columns from left to right are for winter, spring, summer, and fall, respectively. The subscript 1 of and denotes Tianjin and 2 denotes Baoding.
Figure 9: Beta coefficients and elasticity coefficients of bivariate DFA and standard regression model of Tianjin with the same legend as in Fig. 8. Here the subscripts 1 and 2 denote Beijing and Baoding, respectively.
Figure 10: Beta coefficients and elasticity coefficients of bivariate DFA and standard regression model of Baoding with the same legend as in Fig. 8. Here the subscripts 1 and 2 denote Beijing and Tianjin, respectively.
Figure 11: Beta coefficients and elasticity coefficients of bivariate DFA and standard regression model of Beijing. The four columns from left to right are for winter, spring, summer, and fall, respectively. The subscript 1 of and denotes Tianjin and 2 denotes Baoding.
Figure 12: Beta coefficients and elasticity coefficients of bivariate DFA and standard regression model of Tianjin with the same legend as in Fig. 11. Here the subscripts 1 and 2 denote Beijing and Baoding, respectively.
Figure 13: Beta coefficients and elasticity coefficients of bivariate DFA and standard regression model of Baoding with the same legend as in Fig. 11. Here the subscripts 1 and 2 denote Beijing and Tianjin, respectively.