Research Article  Open Access
Jianhua Xu, Yiwen Xu, Chunan Song, "An Integrative Approach to Understand the ClimaticHydrological Process: A Case Study of Yarkand River, Northwest China", Advances in Meteorology, vol. 2013, Article ID 272715, 9 pages, 2013. https://doi.org/10.1155/2013/272715
An Integrative Approach to Understand the ClimaticHydrological Process: A Case Study of Yarkand River, Northwest China
Abstract
Taking the Yarkand River as an example, this paper conducted an integrative approach combining the DurbinWatson statistic test (DWST), multiple linear regression (MLR), wavelet analysis (WA), coefficient of determination (CD), and Akaike information criterion (AIC) to analyze the climatichydrological process of inland river, Northwest China from a multitime scale perspective. The main findings are as follows. (1) The hydrologic and climatic variables, that is, annual runoff (AR), annual average temperature, (AAT) and annual precipitation (AP), are stochastic and, no significant autocorrelation. (2) The variation patterns of runoff, temperature, and precipitation were scale dependent in time. AR, AAT, and AP basically present linear trends at 16year and 32year scales, but they show nonlinear fluctuations at 2year and 4year scales. (3) The relationship between AR with AAT and AP was simulated by the multiple linear regression equation (MLRE) based on wavelet analysis at each time scale. But the simulated effect at a larger time scale is better than that at a smaller time scale.
1. Introduction
The hydrological response to climate change is an important science issue. To well understand this issue, the coupled system of climatichydrological process should be thoroughly studied at different spatial and temporal scales.
Theoretically, a process can be evaluated to determine if they comprise an ordered, deterministic system, an unordered, random system, or a chaotic, dynamic system, and whether change patterns of periodicity or quasiperiodicity exist. However, it is difficult to achieve a thorough understanding of the mechanism of climatichydrological processes [1]. To date, these questions have not received satisfactory answers [2].
Case studies in different countries and regions have suggested that the climatichydrological process is a complex system [3–6]. Therefore, more studies are required to explore the mechanism of climatichydrological process from different perspectives and using different methods. As a result, the climatichydrological process has been explored using various analytical methods, including the fractal theory [7–9], selforganized criticality [10], wavelets analysis [11–13], and artificial neural networks [14, 15]. Although there were several effective methods available to reveal the variations in climatichydrological process [16–19], it has proven difficult to achieve a thorough understanding of the mechanism of climatichydrological process in inland river [2].
In the last 20 years, studies have been conducted to evaluate climate change and hydrological and ecological processes in the arid and semiarid regions in northwestern China [18–25]. Some studies have indicated that there was a visible transition in the hydroclimatic processes in the past halfcentury [24, 26–28]. This transition was characterized by a continual increase in temperature and precipitation, added river runoff volumes, increased lake water surface elevation and area, and elevated groundwater level. This transition may present a series of questions if these changes represent a localized transition to a warm and wet climate type in response to global warming, or merely reflect a centennial periodicity in hydrological dynamics. To date, these questions have not received satisfactory answers; therefore, more studies are required to explore the nonlinear characteristics of hydroclimatic process from different perspectives and using different methods [2, 15, 29].
Though some studies have shown that the inland river in northwest China (NW China), such as the Yarkand River, is mainly recharged by snowmelt, the main climatic factors affecting the streamflow are temperature and precipitation [20, 30, 31]. But due to the complexity of hydroclimatic system, it is difficult to understand the mechanism of climatichydrological process thoroughly [2]. For the above reasons, this paper did not involve the complex physical mechanisms but conducted an integrative approach combining statistics and wavelet analysis to understand the variation of annual runoff and its response to climatic factors at different time scales.
2. Materials and Methods
2.1. Study Area and Data
The Yarkand River is a typical representative of inland rivers, which is located in the Tarim River Basin of Xinjiang Uygur Autonomous Region, northwestern China (Figure 1), with a length of . The Yarkand River (, ) has a total basin area of , including as the mountain area, which accounts for 61.5%, and as the plain area, which takes up 38.5% [31]. The main stream of Yarkand River originates from Karakoram Pass in the north slope of Karakoram Mountain, which is full of towering peaks and glaciers, as well as the extremely rare precipitation in plain. Due to the special geographical conditions, the accumulation of ice and snow in high mountain is the only supply source for runoff. Therefore, the Yarkand River is a typical icesnow supply river, in which the multiyear average runoff in Kaqun hydrometric station consists of 64.0% from mean volume of glacial ablation, from rain and snow supply, and 22.6% from groundwater supply, respectively [32, 33].
For the Yarkand River is an inland river, no water recharges in the plain area, and its stream flow mainly comes from mountainous area, that is, the Pamir Mountains. In other words, the streamflow of the Yarkand River is mainly fed by glacier and snowmelt in the Pamir Mountains. Therefore, the climatic factors, especially temperature and precipitation, directly affect the annual changes in the runoff. So we use the runoff as well as temperature and precipitation data to analyze the climatichydrological process in Yarkand River. The runoff data were from the Kaqun hydrologic station, and temperature and precipitation data were from Tash Kurghan meteorological station. The two stations are located in the source areas of the river; the amount of water used by humans is minimal compared to the total discharge. Therefore, the observed hydrological and meteorological records reflect the natural conditions.
Longterm climate changes can alter the runoff production pattern, the timing of hydrological events, and the frequency and severity of floods, particularly in arid or semiarid regions. Therefore, a small change in precipitation and temperature may result in marked changes in runoff. To investigate the runoff and its related climatic effect, this study used the time series data of annual runoff (AR), annual average temperature (AAT), and annual precipitation (AP) from 1957 to 2008.
2.2. Methods
In order to study the variations of streamflow with regional climate change at different time scales, this paper conducted a comprehensive method including the DurbinWatson statistic test (DWST), multiple linear regression (MLR), wavelet analysis (WA), coefficient of determination (CD), and Akaike information criterion (AIC). Firstly, the DWST was used to explore the stochastic characteristic of hydrologic and climatic variables. Secondly, the WA was used to reveal the variation patterns of annual runoff (AR) and its related climatic factors at different time scales. Thirdly, the relationship between AR with AAT and AP was simulated by MLR based on WA at different time scales. Finally, the estimated effect of multiple linear regression equation (MLRE) at each time scale was tested by CD and AIC.
2.2.1. DurbinWatson Statistic Test
The DurbinWatson statistic is a test statistic used to detect the presence of autocorrelation (a relationship between values separated from each other by a given time lag) in the residuals (prediction errors) from a regression analysis [34, 35].
For a variable , the DurbinWatson statistic is where , and and are, respectively, the observed and predicted values of the response variable for individual ; is the number of observations.
To test for positive autocorrelation at significance α, the test statistic DW is compared to lower and upper critical values ( and ): if , there is statistical evidence that the error terms are positively autocorrelated; if , there is no statistical evidence that the error terms are positively autocorrelated; if , the test is inconclusive.
To test for negative autocorrelation at significance α, the test statistic is compared to lower and upper critical values ( and ): if , there is statistical evidence that the error terms are negatively autocorrelated; if , there is no statistical evidence that the error terms are negatively autocorrelated; if , the test is inconclusive.
Using the DurbinWatson statistic, we checked the autocorrelation of hydrological and climatic variables, such as temperature, precipitation, and runoff.
2.2.2. Wavelet Analysis
Wavelet transformation has been shown to be a powerful technique for characterization of the frequency, intensity, time position, and duration of variations in climate and hydrological time series [11, 12, 16, 36]. Wavelet analysis can also reveal the localized time and frequency information without requiring the time series to be stationary, as required by the Fourier transform and other spectral methods [37].
A continuous wavelet function that depends on a nondimensional time parameter can be written as [36] where denotes time, is the scale parameter, and is the translation parameter. must have a zero mean and be localized in both time and Fourier space [38]. The continuous wavelet transform (CWT) of a discrete signal, , such as the time series of runoff, temperature, or precipitation, is expressed by the convolution of with a scaled and translated , where * indicates the complex conjugate, and denotes the wavelet coefficient. Thus, the concept of frequency is replaced by that of scale, which can characterize the variation in the signal, , at a given time scale.
Selecting a proper wavelet function is a prerequisite for time series analysis [39, 40]. The actual criteria for wavelet selection include selfsimilarity, compactness, and smoothness [41]. Because the symlets are nearly symmetrical, orthogonal, and biorthogonal wavelets proposed by Daubechies as modifications to the db family [41], this study chose the symlets 8 to analyze the variation patterns of runoff and its related climatic factors in the computing environment of MATLAB.
For a time series, , it can be analyzed at multiple scales through wavelet decomposition on the basis of the discrete wavelet transform (DWT). The DWT is defined taking discrete values of and . The full DWT for signal, , can be represented as [42] where and are the flexing and parallel shift of the basic scaling function, , and the mother wavelet function, , and and are the scaling coefficients and the wavelet coefficients, respectively. Generally, scales and positions are based on powers of 2, which is the dyadic DWT.
Once a mother wavelet is selected, the wavelet transform can be used to decompose a signal according to scale, allowing separation of the finescale behavior (detail) from the largescale behavior (approximation) of the signal [43]. The relationship between scale and signal behavior is designated as follows: low scale corresponds to compressed wavelet as well as rapidly changing details, namely, high frequency; whereas high scale corresponds to stretched wavelet and slowly changing coarse features, namely, low frequency. Signal decomposition is typically conducted in an iterative fashion using a series of scales such as , with successive approximations being split in turn so that one signal is broken down into many lower resolution components.
The wavelet decomposition and reconstruction were used to approximate the variation patterns of AR and its related factors over the entire study period at the selected different time scales.
2.2.3. Multiple Linear Regression
For understanding the relationship between annual runoff with its related climatic factors at different time scales, we employed multiple linear regression (MLR) based on wavelet analysis. This method fits multiple linear regression equation (MLRE) between AR with AAT and AP by using multiple linear regression (MLR) based on the results of wavelet approximation [44].
The multiple linear regression model is where is dependent variable, is the independent variables; is the regression coefficient, which is generally calculated by method of least squares [45]. In this study, the dependent variable is the annual runoff (AR) and the independent variables are related climatic factors, such as the annual average temperature (AAT) and annual precipitation (AP).
2.2.4. Coefficient of Determination and Akaike Information Criterion
In order to identify the uncertainty of the estimated model for a given time scale, the coefficient of determination, also known as the goodness of fit, was calculated as follows: where is the coefficient of determination; and are the simulated value and actual data of runoff, respectively; is the mean of is the residual sum of squares; is the total sum of squares.
The coefficient of determination is a measure of how well the simulated results represent the actual data. A bigger indicates a higher certainty and lower uncertainty of the estimates [45].
To compare the relative goodness between the ANN and multiple linear regression (MLR) fit for a given timescale, we also used the measure of Akaike information criterion (AIC) [46]. The formula of AIC is as follows: where is the number of parameters estimated in the model; is the number of samples; RSS is the same as in formula (6). A smaller AIC indicates a better model.
For small sample sizes (i.e., ), the secondorder Akaike information criterion () should be used instead: where is the sample size. As the sample size increases, the last term of the approaches zero, and the tends to yield the same conclusions as the AIC [47].
3. Results
3.1. Check for Variable’s Autocorrelation
The premise of statistics indicates that models imply an assumption; that is, variables are stochastic and no significant autocorrelation is present. Is it really? This can be demonstrated by statistical check for variable’s autocorrelation (Table 1).
 
Notes. AR: annual runoff, AAT: annual average temperature, and AP: annual precipitation; the subscripts, and , represent time. 
For (1*), (2*), and (3*) in Table 1, their degree of freedom is, respectively, equals 2 and (i.e., ) equals 51. Upper and lower critical values of the DurbinWatson Statistic (DW) are equals 1.509 and equals 1.58 when significance level (α) equals 0.01. Because the values of DW are between and , it is obvious that annual runoff (AR), annual average temperature (ATT), and annual precipitation (AP) indicate no autocorrelation.
For (4*) in Table 1, its degree of freedom is, respectively, equals 4 and equals 51. Upper and lower critical values of DW are equals 1.25350 and equals 1.49384 when significance level equals 0.01. For (5*), its degree of freedom is, respectively, equals 6 and equals 51. Upper and lower critical values of DW are equals 1.17372 and equals 1.58811 when significance level equals 0.01. Thereby they indicate no autocorrelation either.
In fact, it can be determined that variables and model reveal nonautocorrelation because the value of DW is close to 2 for each regression equation shown in Table 1. Therefore, the assumption of our model is logical.
3.2. Variation Patterns of ClimaticHydrological Process at Different Time Scales
Our previous study indicated that [44] the annual average temperature and annual precipitation are the most important factors that related with the annual runoff. The result was also supported by the other studies for the headwaters of the Tarim River Basin [20, 30–33].
The raw data of AR, AAT, and AP showed fluctuation. It is difficult to identify any patterns simply based on the raw data. In order to show the scaledependent with time for the climatichydrological process of the Yarkand River, the wavelet analysis was used. The nonlinear variation for the annual runoff process and the related climate factors were analyzed at multipleyear scales through wavelet decomposition on the basis of the discrete wavelet transform (DWT).
The wavelet decomposition for the time series of annual runoff at five time scales resulted in five variation patterns (Figure 2(a)). The S1 curve retains a large amount of residual from the raw data, and drastic fluctuations exist in the period from 1957 to 2008. These characteristics indicate that although the runoff varied greatly throughout the study period, there was a hidden increasing trend. The S2 curve still retains a considerable amount of residual, as indicated by the presence of 4 peaks and 4 valleys. However, the S2 curve is much smoother than the S1 curve, which allows the hidden increasing trend to be more apparent. The S3 curve retains much less residual, as indicated by the presence of 2 peaks and 2 valleys. Compared to S2, the increase in runoff over time is more apparent in S3. Finally, the S5 curve presents an ascending tendency, whereas the increasing trend is obvious in the S4 curve.
(a)
(b)
(c)
Accordingly, Figures 2(b) and 2(c) provide us with a method for comparing the variation patterns of annual average temperature and annual precipitation at different time scales. The wavelet decomposition for the time series at five time scales resulted in five variation patterns, respectively. These five time scales are also designated as S1 to S5. The curves present an ascending tendency despite drastic fluctuations in S1 and S2. Then, the curves are getting much smoother and the increasing trend becomes even more obvious as the scale level increases.
The upper analysis showed that the nonlinear variations of runoff, temperature, and precipitation of the Yarkand River basin were dependent on time scales. The annual runoff, annual average temperature and annual precipitation at five time scales resulted in five patterns of nonlinear variations, respectively.
3.3. Relationship between Streamflow and Climate Factors
Based on the raw data of AAT, AP, and AR, multiple linear regression equation (MLRE) was developed as follows:
Equation (9) reveals a positive correlation between the annual runoff and the annual average temperature. These results are readily supported by the fact that the majority of streamflow comes from glacial melt and snowmelt, which have been occurring at increased rates as the temperature increases. These results have been confirmed by other studies [48]. However, (9) also indicates the existence of a weak, negative correlation between the annual runoff and the annual precipitation, which does not seem reasonable. Indeed, this finding conflicts with the results of other studies, which have suggested that both the temperature and precipitation series in the Tarim Basin have been increasing in a pattern similar to that of annual runoff over the past 50 years [20, 30].
Additionally, the coefficient of determination of (9) is as low as 0.093. Furthermore, the average absolute error and average relative error for predicted results is as high as and 14.11%, respectively. All this means that the regression (9) is not authentic. What is the reason for this? It is possible that this inconsistency is caused by randomicity in the raw timeseries data, which should be filtered out via wavelet decomposition based on the discrete wavelet transform [17, 18].
To understand the response of the runoff to regional climatic change at different time scales, based on the results of wavelet decomposition (Figure 2), multiple linear regression equation (MLRE) at each time scale was fitted for describing the relationship among annual runoff, annual average temperature, and annual precipitation (Table 2).
 
Notes. AR: annual runoff, AAT: annual average temperature, and AP: annual precipitation. 
3.4. Comparison of the Estimated Results at Different Time Scales
Though all MLREs at each time scale in Table 2 got across the statistical test at the significant level of 0.01 or 0.001, the predicted error of MLRE at the chosen time scales is different. Figure 3 shows the comparison for the simulated value by MLREs and original data of AR at different time scales. The predicted error of MLRE at the time scale of S1 and S2 (i.e., 2year and 4year scales) is large, that at the time scale of S3 (i.e., 8year scale) is also fairish, and wee predicted error of MLRE at the time scale of S4 and S5 (i.e., 4year and 5year scales) appears. These results show that MLRE only can well fit the relationship between runoff and climate factors at large time scale such as at 16year and 32year scales.
(a) S1 (the time scale of 2year)
(b) S2 (the time scale of 4year)
(c) S3 (the time scale of 8year)
(d) S4 (the time scale of 16year)
(e) S5 (the time scale of 32year)
By comparing the and AIC value in Table 3, we can know the estimated effect (good or bad) of models at different time scales.

Table 3 tells us that the for MLRE at the time scale of S1 and S2 (i.e., 2year and 4year scales) is lower (only 0.361 and 0.416, resp.) and that at the time scale of S3 (i.e., 8year scales) is higher, reaching 0.894. Only the MLRE at the time scale of S4 and S5 (i.e., 16year and 32year scales) has the high coefficient of determination, which is as high as 0.975 and 0.996, respectively.
The lower AIC value indicates better model, which tells us that the MLRE at time scale of S5 is the best, that at time scale of S4 is better, that at time scale of S3 is moderate, that at time scale of S2 is the penult, and that at time scale of S1 is the worst.
Overall, the relationship between AR with AAT and AP was simulated by the multiple linear regression (MLR) based on wavelet analysis at different time scales, but the simulated effect at a larger time scale is better than that at a smaller time scale.
4. Discussion and Conclusions
Many studies indicated that the climatichydrological process is a complex system with nonlinearity, but it is difficult to understand the mechanism of climatichydrological process thoroughly [2]. Our results showed the following fact: the simulated effect at large time scale is better than that at small time scale, and the estimated precision at large time scale is higher than that at small time scale. What are the causes for this? It is difficult to thoroughly answer the question because of the nonlinear complicated climatichydrological process, which is essentially difficult to precisely predict [15].
Our study revealed that the climatichydrological process at larger time scale (e.g., 16year or 32year scales) basically presented a linear process, but that at smaller time scale (e.g., 2year or 4year scales) is essentially nonlinear process with complicated causations. Because the time series of runoff are essentially monotonic trends related to longterm climatic changes at large time scale (e.g., 16year and 32year scales), the estimated precision is much higher. Otherwise, due to the difficulty to accurately predict nonlinear climatichydrological process at small time scales (e.g., 2year or 4year scale), the estimated precision and simulated effect are not satisfactory.
Nevertheless, the comprehensive method conducted by this paper provides a method to understand the climatichydrological process in the Yarkand River from the perspective of multiscale, which may be used to explore the climatichydrological process in other inland rivers of northwest China.
The main conclusions of this work can be summarized as follows.(1)The hydrologic and climatic variables, that is, annual runoff (AR), annual average temperature (AAT), and annual precipitation (AP) are stochastic and show no significant autocorrelation.(2)The variation pattern of runoff, temperature, and precipitation was scale dependent with time. The annual runoff (AR), annual average temperature (AAT), and annual precipitation (AP) basically present linear trends at 16year and 32year scales, but they show nonlinear fluctuations at 2year, 4year, and 8year scales. (3)The relationship between AR with AAT and AP was simulated by the multiple linear regression equation (MLRE) based on wavelet analysis at each time scale. The results showed that the AR is basically monotonic trend related to longterm climatic changes at a larger time scale (e.g., 16year or 32year scales), and the estimated precision is much higher. But due to an essentially nonlinear climatichydrological process at a smaller time scale (e.g., 2year or 4year scales), the estimated precision is lower than that at a larger time scale.
Acknowledgments
This work was supported by the Director Fund of the Key Lab of GIScience of the Education Ministry PRC. The authors are grateful to the editor and referees whose comments helped us in improving the article’s quality.
References
 A. J. Cannon and I. G. McKendry, “A graphical sensitivity analysis for statistical climate models: application to Indian monsoon rainfall prediction by artificial neural networks and multiple linear regression models,” International Journal of Climatology, vol. 22, no. 13, pp. 1687–1708, 2002. View at: Publisher Site  Google Scholar
 J. H. Xu, Y. N. Chen, W. H. Li, Q. Nie, Y. L. Hong, and Y. Yang, “The nonlinear hydroclimatic process in the Yarkand River, northwestern China,” Stochastic Environmental Research and Risk Assessment, vol. 27, no. 2, pp. 389–399, 2013. View at: Google Scholar
 R. Ibbitt and R. Woods, “Rescaling the topographic index to improve the representation of physical processes in catchment models,” Journal of Hydrology, vol. 293, no. 1–4, pp. 205–218, 2004. View at: Publisher Site  Google Scholar
 W. G. Strupczewski, V. P. Singh, S. Weglarczyk, K. Kochanek, and H. T. Mitosek, “Complementary aspects of linear flood routing modelling and flood frequency analysis,” Hydrological Processes, vol. 20, no. 16, pp. 3535–3554, 2006. View at: Publisher Site  Google Scholar
 B. Sivakumar, “Nonlinear determinism in river flow: prediction as a possible indicator,” Earth Surface Processes and Landforms, vol. 32, no. 7, pp. 969–979, 2007. View at: Google Scholar
 J. H. Xu, Y. N. Chen, W. H. Li, and S. Dong, “Longterm trend and fractal of annual runoff process in mainstream of Tarim River,” Chinese Geographical Science, vol. 18, no. 1, pp. 77–84, 2008. View at: Publisher Site  Google Scholar
 B. P. Wilcox, M. S. Seyfried, and T. H. Matison, “Searching for chaotic dynamics in snowmelt runoff,” Water Resources Research, vol. 27, no. 6, pp. 1005–1010, 1991. View at: Publisher Site  Google Scholar
 J. W. Kantelhardt, D. Rybski, S. A. Zschiegner et al., “Multifractality of river runoff and precipitation: comparison of fluctuation analysis and wavelet methods,” Physica A, vol. 330, no. 12, pp. 240–245, 2003. View at: Publisher Site  Google Scholar
 J. H. Xu, Y. N. Chen, W. H. Li, M. H. Ji, and S. Dong, “The complex nonlinear systems with fractal as well as chaotic dynamics of annual runoff processes in the three headwaters of the Tarim River,” Journal of Geographical Sciences, vol. 19, no. 1, pp. 25–35, 2009. View at: Publisher Site  Google Scholar
 Z. L. Wang and C. Y. Huang, “Selforganized criticality of rainfall in central China,” Advances in Meteorology, vol. 2012, Article ID 203682, 8 pages, 2012. View at: Publisher Site  Google Scholar
 L. C. Smith, D. L. Turcotte, and B. L. Isacks, “Streamflow characterization and feature detection using a discrete wavelet transform,” Hydrological Processes, vol. 12, no. 2, pp. 233–249, 1998. View at: Google Scholar
 C. M. Chou, “Efficient nonlinear modeling of rainfallrunoff process using wavelet compression,” Journal of Hydrology, vol. 332, no. 34, pp. 442–455, 2007. View at: Google Scholar
 J. H. Xu, Y. N. Chen, M. H. Ji, and F. Lu, “Climate change and its effects on runoff of Kaidu River, Xinjiang, China: a multiple timescale analysis,” Chinese Geographical Science, vol. 18, no. 4, pp. 331–339, 2008. View at: Publisher Site  Google Scholar
 C. H. Hu, Y. H. Hao, T. C. J. Yeh, B. Pang, and Z. N. Wu, “Simulation of spring flows from a karst aquifer with an artificial neural network,” Hydrological Processes, vol. 22, no. 5, pp. 596–604, 2008. View at: Publisher Site  Google Scholar
 J. H. Xu, Y. N. Chen, W. H. Li et al., “Combining BPANN and wavelet analysis to simulate hydroclimatic process—a case study of the Kaidu River, NW China,” Frontiers of Earth Science, 2013. View at: Publisher Site  Google Scholar
 J. H. Xu, Y. N. Chen, W. H. Li, M. H. Ji, S. Dong, and Y. L. Hong, “Wavelet analysis and nonparametric test for climate change in Tarim River Basin of Xinjiang during 1959–2006,” Chinese Geographical Science, vol. 19, no. 4, pp. 306–313, 2009. View at: Publisher Site  Google Scholar
 J. H. Xu, W. H. Li, M. H. Ji, F. Lu, and S. Dong, “A comprehensive approach to characterization of the nonlinearity of runoff in the headwaters of the Tarim River, western China,” Hydrological Processes, vol. 24, no. 2, pp. 136–146, 2010. View at: Publisher Site  Google Scholar
 J. H. Xu, Y. N. Chen, F. Lu, W. H. Li, L. J. Zhang, and Y. L. Hong, “The nonlinear trend of runoff and its response to climate change in the Aksu River, western China,” International Journal of Climatology, vol. 31, no. 5, pp. 687–695, 2011. View at: Google Scholar
 J. H. Xu, Y. N. Chen, W. H. Li, Y. Yang, and Y. L. Hong, “An integrated statistical approach to identify the nonlinear trend of runoff in the Hotan River and its relation with climatic factors,” Stochastic Environmental Research and Risk Assessment, vol. 25, no. 2, pp. 223–233, 2011. View at: Publisher Site  Google Scholar
 X. M. Hao, Y. N. Chen, C. C. Xu, and W. H. Li, “Impacts of climate change and human activities on the surface runoff in the Tarim River Basin over the last fifty years,” Water Resources Management, vol. 22, no. 9, pp. 1159–1171, 2008. View at: Publisher Site  Google Scholar
 Y. N. Chen and Z. X. Xu, “Plausible impact of global climate change on water resources in the Tarim River Basin,” Science in China D, vol. 48, no. 1, pp. 65–73, 2005. View at: Google Scholar
 Y. N. Chen, K. Takeuchi, C. C. Xu, Y. P. Chen, and Z. X. Xu, “Regional climate change and its effects on river runoff in the Tarim Basin, China,” Hydrological Processes, vol. 20, no. 10, pp. 2207–2216, 2006. View at: Publisher Site  Google Scholar
 Y. N. Chen, C. C. Xu, X. M. Hao et al., “Fiftyyear climate change and its effect on annual runoff in the Tarim River Basin, China,” Quaternary International, vol. 208, no. 12, pp. 53–61, 2009. View at: Google Scholar
 J. Wang, H. Li, and X. Hao, “Responses of snowmelt runoff to climatic change in an inland river basin, Northwestern China, over the past 50 years,” Hydrology and Earth System Sciences, vol. 14, no. 10, pp. 1979–1987, 2010. View at: Publisher Site  Google Scholar
 Q. Zhang, C. Y. Xu, H. Tao, T. Jiang, and Y. D. Chen, “Climate changes and their impacts on water resources in the arid regions: a case study of the Tarim River basin, China,” Stochastic Environmental Research and Risk Assessment, vol. 24, no. 3, pp. 349–358, 2010. View at: Publisher Site  Google Scholar
 Y. F. Shi, Y. P. Shen, E. S. Kang et al., “Recent and future climate change in northwest china,” Climatic Change, vol. 80, no. 34, pp. 379–393, 2007. View at: Google Scholar
 B. F. Li, Y. N. Chen, Z. S. Chen, and W. H. Li, “Trends in runoff versus climate change in typical rivers in the arid region of northwest China,” Quaternary International, vol. 282, pp. 87–95, 2012. View at: Google Scholar
 B. F. Li, Y. N. Chen, and X. Shi, “Why does the temperature rise faster in the arid region of northwest China?” Journal of Geophysical Research, vol. 117, no. 16, Article ID 16115, 2012. View at: Publisher Site  Google Scholar
 Y. Yang, J. H. Xu, Y. L. Hong, and G. H. Lv, “The dynamic of vegetation coverage and its response to climate factors in Inner Mongolia, China,” Stochastic Environmental Research and Risk Assessment, vol. 26, no. 3, pp. 357–373, 2012. View at: Publisher Site  Google Scholar
 Y. N. Chen, K. Takeuchi, C. C. Xu, Y. P. Chen, and Z. X. Xu, “Regional climate change and its effects on river runoff in the Tarim Basin, China,” Hydrological Processes, vol. 20, no. 10, pp. 2207–2216, 2006. View at: Publisher Site  Google Scholar
 B. G. Sun, W. Y. Mao, Y. R. Feng et al., “study on the change of air temperature, precipitation and runoff volume in the Yarkant River basin,” Arid Zone Research, vol. 23, no. 2, pp. 203–209, 2006 (Chinese). View at: Google Scholar
 M. Sabit and A. Tohti, “An analysis of water resources and it's hydrological characteristic of Yarkend River Valley,” Journal of Xinjiang Normal University, vol. 24, no. 1, pp. 74–78, 2005 (Chinese). View at: Google Scholar
 T. L. Liu, Q. Yang, R. Qin, Y. P. He, and R. Liu, “Climate change towards warmingwetting trend and its effects on runoff at the headwater region of the Yarkand River in Xinjiang,” Journal of Arid Land Resources and Environment, vol. 22, no. 9, pp. 49–53, 2008 (Chinese). View at: Google Scholar
 J. Durbin and G. S. Watson, “Testing for serial correlation in least squares regression. I,” Biometrika, vol. 37, no. 34, pp. 409–428, 1950. View at: Google Scholar
 J. Durbin and G. S. Watson, “Testing for serial correlation in least squares regression. II,” Biometrika, vol. 38, no. 12, pp. 159–178, 1951. View at: Google Scholar
 D. Labat, “Recent advances in wavelet analyses: part 1. A review of concepts,” Journal of Hydrology, vol. 314, no. 1–4, pp. 275–288, 2005. View at: Google Scholar
 C. Torrence and G. P. Compo, “A practical guide to wavelet analysis,” Bulletin of the American Meteorological Society, vol. 79, no. 1, pp. 61–78, 1998. View at: Google Scholar
 M. Farge, “Wavelet transforms and their applications to turbulence,” Annual Review of Fluid Mechanics, vol. 24, no. 1, pp. 395–457, 1992. View at: Google Scholar
 J. B. Ramsey, “Regression over timescale decompositions: a sampling analysis of distributional properties,” Economic Systems Research, vol. 11, no. 2, pp. 163–183, 1999. View at: Google Scholar
 J. H. Xu, Y. Lu, F. L. Su, and N. S. Ai, “R/S and wavelet analysis on the evolutionary process of regional economic disparity in China during the past 50 years,” Chinese Geographical Science, vol. 14, no. 3, pp. 193–201, 2004. View at: Google Scholar
 I. Daubechies, “Orthonormal bases of compactly supported wavelets,” Communications on Pure and Applied Mathematics, vol. 41, pp. 909–996, 1988. View at: Google Scholar
 S. G. Mallat, “Theory for multiresolution signal decomposition: the wavelet representation,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 11, no. 7, pp. 674–693, 1989. View at: Publisher Site  Google Scholar
 L. M. Bruce, C. H. Koger, and J. Li, “Dimensionality reduction of hyperspectral data using discrete wavelet transform feature extraction,” IEEE Transactions on Geoscience and Remote Sensing, vol. 40, no. 10, pp. 2331–2338, 2002. View at: Publisher Site  Google Scholar
 J. H. Xu, Y. N. Chen, M. H. Ji, and F. Lu, “Climate change and its effects on runoff of Kaidu River, Xinjiang, China: a multiple timescale analysis,” Chinese Geographical Science, vol. 18, no. 4, pp. 331–339, 2008. View at: Publisher Site  Google Scholar
 J. H. Xu, Mathematical Methods in Contemporary Geography, Higher Education Press, Beijing, China, 2002, (in Chinese).
 D. R. Anderson, K. P. Burnham, and W. L. Thompson, “Null hypothesis testing: problems, prevalence, and an alternative,” Journal of Wildlife Management, vol. 64, no. 4, pp. 912–923, 2000. View at: Google Scholar
 K. P. Burnham and D. R. Anderson, Model Selection and Multimodel Inference: A Practical InformationTheoretic Method, Springer, New York, NY, USA, 2nd edition, 2002.
 X. Wang, Z. C. Xie, S. Y. Liu, D. H. Shangguan, J. J. Tao, and Y. L. Yang, “Prediction on the variation trend of glacier system in the source region of Tarim River responding to climate change,” Journal of Mountain Research, vol. 24, no. 6, pp. 641–646, 2006 (Chinese). View at: Google Scholar
Copyright
Copyright © 2013 Jianhua Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.