2016 International Conference on Environment, Climate Change and Sustainable Development (ECCSD 2016) ISBN: 978-1-60595-358-8

Multiscale Correlation Analysis between Runoff and Sediment and Possible Underlying Causes in in a Changing Environment Jun SHAO*, Ming-long DAI, Shuo OUYANG and Ke-xu FAN Bureau of Hydrology, Changjiang Water Resources Commission, Wuhan 430010 *Corresponding author

Keywords: Sediment, runoff, Correlation analysis, Cross wavelet transform, Jinsha River.

Abstract. Great changes have taken place in the sediment processes in Jinsha River because of human activities. The wavelet analysis and the cross wavelet transform was applied to analyze the variation of streamflow and sediment in Shigu, Panzhihua, and Pinshan hydrometric station. The multiscale correlation between Jinsha River sediment and runoff data was discussed. The result showed that the power spectrum of streamflow had the same variation in Shigu, Panzhihua, and Pinshan stations. The trends exhibited a great difference for sediment changes. For Shigu and Panzhihua stations, visible regions of six-month and below six-month bands in the entire time coordinate, which is a great significance period for sediment observation, were found. The significant variation was especially powerful from 1987 to 2009. These findings indicated that soil erosion in the upper and middle streams of Jinsha River became more serious starting in the 1980s. The sediment load in Pinshan station exhibited another trend over the significant periods. The six-month fluctuations disappeared from 2000 to 2009. This result showed that sediment content decreased from the Panzhihua to the Pinshan range. The sediment content of Jinsha River will continue to decrease in the future with ongoing water and soil conservation projects and hydropower development.

Introduction Jinsha River spans for more than ten latitudes and its altitude differs by more than 4000 meters. This river is one of the major sediment-formation areas in the upper reaches of River. Great changes have taken place in the streamflow and sediment characteristics of Jinsha River in the past few decades because of natural factors and human activities. These changes directly affected reservoir operation, water resource development, and aquatic ecosystem protection. Variations in runoff and sediment concentration in Jinsha River affected the incoming runoff and sediment condition in the Three Gorges reservoir (Chen et al., 2008; Feng et al., 2008). Many studies have been conducted on the characteristics of sediments in Jinsha River. The sediment loads in Jinsha River originated from different source areas. The effects of precipitation, soil and water conservation, and reservoir on sediment production at various source areas were also different (Chen et al., 2008). Fluvial sediments were chiefly produced in the lower reaches of the river, especially in the main river valley of the lower reaches (Liu et al., 1999). The annual runoff 8 and sediment transport in Pingshan station of Jinsha River increased by about 0.35×10 t compared with their value before 1990. The percentage also increased from 51.8% to 71.9% compared with that of Yichang station (Xu et al., 2004). Suspended sediment load and concentration in Jinsha River showed no trending variation from 1954 to 1975. However, suspended sediment load and concentration increased after this period (Huang et al., 2002; Xu et al., 2007, 2008). The annual suspended load in Jinsha River from 1976 to 1996 was higher than that in 1954 to 1975 on the condition that annual river flows for both periods were equal. This condition was due to the increase in soil loss by large-scale mining and the construction of the iron-steel industry and communication lines (Xu et al., 2007, 2008; Zhang et al., 2007). After 2000, the annual sediment load in Jinsha River decreased (Fu et al., 2006; Xu, 2009). Spectral analysis, particularly wavelet transform, has been successfully applied in sediment research in Yangtze River (Zhang et al., 2006, 2008). However, it’s not easy for wavelet transform to describe the relationship between two time series. We need to determine the bivariate extension of a wavelet transform when analyzing the relationship between two time series or comparing two different variables such as streamflow and sediment load. A cross wavelet transform could show common power and relative phase in both time and frequency space (Torrence et al., 1998; Grinsted et al., 2004). The objectives of this study are to determine the multiple scale characteristics of runoff and sediment and to elucidate the relationship between sediment concentration and human activities. We achieved these objectives using wavelet transform and cross wavelet transform approaches.

Study Area and Data Jinsha River, including Tuotuo River and Tongtian River, lie between Qinghai-Tibet Plateau, Hengduan Mountains, Yunnan-Guizhou Plateau, and Basin. Jinsha River is an important part of Yangtze River. The former has a total length of 3500 km and a drainage basin area of 500000 km². Jinsha River extends across portions of Qinghai, Tibet, Sichuan, Yunnan, and Guizhou. The upper reaches of Jinsha River extend from Zhimenda to Shigu, with a length of 984 km and a drainage area covering 76500 km². The middle Jinsha River basin that is upstream of Panzhihua is about 564 km long with a drainage area of 45000 km². The downstream of Panzhihua includes the lower reaches of Jinsha River, which extends for 768 km with a drainage area of 214000 km². Annual and monthly suspended sediment concentration (SSC) and streamflow data from the three main gauging stations of Jinsha River were analyzed in this paper. Shigu, Panzhihua, and Pinshan stations represented the upper, middle, and lower reaches of Jinsha River, respectively.

Methodology A wavelet transform is a useful time-frequency analysis method that is widely used to analyze nonstationary hydrological time series at different frequencies (Daubechies, 1990; Kumar et al., 1997; Torrence et al., 1998). A continuous wavelet transform Wn is the convolution of a vector x (with time dimension n) with a wavelet function Ψ

N −1 − δ = Ψ (n ' n ) t  Wn( s ) ∑ x n '   n'= 0 s  where s is the scale or dilation, n’-n shows the number of points from the time series origin, δt is the time interval, N is the number of points, and the overbar indicates the complex conjugate. We chose the complex Morlet wavelet function.

Results and Discussion Runoff and Sediment Load Changes in the Three Stations The continuous wavelet power spectrum of the monthly runoff and monthly SSC represent the average variance over a range of scales. This representation provides an efficient approach in examining river streamflow and sediment load variations. Fig. 1 shows the wavelet power spectrum for Shigu station. Fig. 1a represents runoff and Fig. 1b represents sediment load. The thick contours depict the 95% confidence level of local power relative to red noise. Figure 1a reveals that the runoff had many kinds of quasi-periodic oscillations, which include twelve-month, six-month, and below six-month fluctuations. The twelve-month band penetrated the entire time period from 1958 to 2009, which was the main oscillation of Shigu station. The six-month and below six-month fluctuations were sporadically distributed along the runoff series. Fig. 1b typically shows the temporal fluctuations of SSC. The figure gives a measure that is similar to temporal fluctuations for the runoff, but the fluctuation was not continuous over the entire time series unlike that of the runoff sediment load. The 5% significant region indicates that the time

series 1961 to 1967, 1971 to 1976, and 1980 to 2009 exhibited intervals of higher sediment load variance, whereas the time series 1958 to 1960, 1968 to 1970, and 1977 to 1979 demonstrated lower SSC variance.

(a) (b) Figure 1. Continuous wavelet power spectrum of Shigu station (a) monthly runoff and (b) monthly SSC standardized time series using the Morlet wavelet.

(a) (b) Figure 2. Continuous wavelet power spectrum of Panzhihua station (a) monthly runoff and (b) monthly SSC standardized time series using the Morlet wavelet.

(a) (b) Figure 3. Continuous wavelet power spectrum of Pinshan station (a) monthly runoff and (b) monthly SSC standardized time series using the Morlet wavelet. The power spectrum for the runoff and sediment load of Panzhihua station (Figures. 2a, 2b) and Pinshan station (Figures. 3a, 3b) was the same as that of Shigu station, with almost in the twelve-month band, and some sparsely distributed below six-month fluctuations. For the three hydrological stations, the runoff exhibited similar variance law as the trends showed a great amount of difference in sediment changes. In Panzhihua station (Figure. 3b), high significant regions were found in the 12-month band, which ran through the entire study period. High power was centered from 1983 to 2005, at individual six-month periodicity. Stronger 12-month oscillations in runoff and sediment variation were recorded in Pinshan station (Figure. 3b) from 1954 to 2003. The 12-month periodicity from 2003 to 2009 was not significant enough at the 95% confidence level. This result revealed that sediment load started to decrease after 2003 in Pinshan station. High-power fluctuations always corresponded with flood sessions whether in streamflow or in sediment concentration wavelet transform. Sediment concentration wavelet transform is more complicated than streamflow, which explains why watershed sediment load dynamics is a complex system. Relationship between Runoff and Sediment Load Applying cross wavelet transform between streamflow and SSC is useful in determining the relationship between the two variables. Fig. 4 shows the cross wavelet spectrum of Shigu, Panzhihua, and Pinshan stations. The thick gray areas stand out as significant, and phase relationship is depicted as arrows in these areas. Figure 4(a) indicates that a significant covariance of the 12-month band was found between monthly runoff and SSC from 1958 to 2009. The phase relationship between the variables was quite easy to ascertain within the significant areas. Visible regions with statistical significance of six-month and below six-month bands were found in the entire time coordinate. The significance from 1987 to 2009 was more powerful because soil erosion in upstream Jinsha River became more serious starting in the 1980s. Fig. 4(b) shows the cross wavelet transform spectrum for Panzhihua station. Power was broadly distributed from the six-month to the twelve-month band in Panzhihua station. A growing trend in sediment concentration was also recorded in the late 1980s. Sediment concentration in Pinshan station presented another trend over the significant periods. Aside from the same regularity exhibited by the two previous stations, the six-month fluctuations also disappeared from 2000 to 2009 in Pinshan station. This finding showed that sediment content decreased from the Panzhihua to Pinshan range.

(a) (b)

(c) Figure 4. Cross wavelet transform between monthly runoff and SSC: (a) Shigu station, (b) Panzhihua station, and (c) Pinshan station. The effect of watershed sediment production was mainly caused by natural factors and human activities (Xu et al., 2004). Natural factors include precipitation, geomorphology, climate, and soil and vegetation conditions. Human activities include reservoir construction, soil and water conservation engineering, and river sand mining. Studies have shown that landslides and mudslides brought a huge amount of sediments into the river, which was the main reason for the high rate of sediment production in Jinsha River. The upper reaches of Jinsha River were less affected by human activities. Thus, the sediment concentration in that area was low. The middle and lower reaches of Jinsha River, which are located in Hengduan Mountains, were the transition areas for both Qinghai-Tibet Plateau and Yunnan-Guizhou Plateau to the Sichuan Basin. The sediments were mainly from the Panzhihua to Pingshan range. This area accounted for only 7.8% of Yichang station but produced 35.5% of the sediments. More seasonal and torrential rains led to severe ground erosion. Mudslides, landslides, and other geological disasters resulted in heavy sediment concentration, aside from surface erosion

and gully erosion. Human activities further increased sediment concentration especially in the 1960s, such as mining, mineral processing, and iron-to-steel manufacturing. While a large-scale mining and steel enterprise was constructed in the region. This is the now famous Panzhihua Iron and Steel Base, which is a medium-sized industrial city in Jinsha River valley. The aforementioned process was also accompanied by large-scale railway and highway construction. Transportation, mining, and city construction not only discharged a huge amount of rock and soil debris into Jinsha River and its tributaries, but also generated a number of mudslides and landslides. These factors greatly increased sediment concentration in Jinsha River. However, the amount of sediments coming to Jinsha River was reduced since 2000. The key soil erosion control project in upstream Yangtze River basin and the construction of cascade reservoirs were mainly responsible for the reduction. Terracing, tree planting, protecting existent vegetation, and returning farmland to forest also promoted soil and water conservation. Thus, the amount of sediments deposited into the river was reduced with surface runoff. By the end of 2008, more than 5000 small watersheds and 96000 km 2 of soil erosion area had been completely governed. In the past 20 years, soil erosion area in Jinsha River has decreased by 15%. Meanwhile, Ertan Reservoir, with a capacity of 5.79 billion cubic meters, began to be filled with water in 1998. In addition, many medium-sized reservoirs were built on the tributaries of Jinsha River. Thus, the amount of sediments in Pinshan station gradually reduced because of reservoir storage effect. In the future, the state plans to construct more cascade hydropower stations on Jinsha River. Several of these stations will be operating soon. The sediment concentration in Jinsha River will continue to decrease with the ongoing water and soil conservation projects and hydropower development.

Conclusions Runoff and sediment load are key elements in river dynamics system. In Jinsha River, great changes in redistribution have taken place in sediment processes because of human activities, especially the construction of a large number of cascade hydropower stations. These new features have brought many changes in reservoir operation, water resource development, and environmental protection. This paper addressed an important scientific problem about the multiscale relationship between sediment load and streamflow. With the development of water resources, hydropower stations, and soil-erosion control projects, human influence became the most important factor that affects the relationship between runoff and sediment load. The interaction and coupling mechanism among precipitation, sediment generation, sediment transportation, reservoir storage effect, and soil erosion control measures need further study.

References [1] Chen. S. S., Zhang. O. Y., Chen. Z. F., Peng. W. B. 2008. Variations of runoff and sediment load of the Jinsha river. Advances in water science , 19(4), 475-482. (in Chinese) [2] Daubechies. I. 1990. The wavelet transform time-frequency localization and signal analysis. IEEE Transactions on Information Theory , 36, 961-1004. [3] Feng. X. F., Yang. Q. Y., Zhang. O. Y., Wang. X. K. 2008. Sediment trap of Ertan reservoir and its effect on sediment budget of Jinsha River. Journal of Sichuan University (Engineering Science Edition) , 40(6), 37-42. (in Chinese) [4] Fu. R. S., Qi M. L., Fang. H. W., Chen. X. Q. 2006. Incoming runoff and sediment runoff to Yichang hydrometric station affected by the projects on the upstream of Yangtze River, Journal of Hydroelectric Engineering , 25(6),103-110, 118. (in Chinese) [5] Huang. C., Lou. X. P., Liu. Y. Y. 2002. Analysis of the sediment development and trend in Jinshajiang River drainage area. Journal of Chongqing University(Natural Science Edition) , 25(1), 21-23. (in Chinese) [6] Kumar. P., Foufoula. G. E. 1997. Wavelet analysis for geophysical applications. Reviews of Geophysics , 35(4), 385-412. [7] Labata. D., Ababoua. R., Manginb. A. 2000. Rainfall-runoff relations for karstic springs. Part II: continuous wavelet and discrete orthogonal multiresolution analyses. Journal of Hydrology , 238, 149-178. [8] Liu. S. Q., Chen. Z. J., Chen. G. J., Yang. D. G. 1999. Surface erosion and fluvial silt in the Jinsha River watershed. Resources and Environment in the Yangtze Basin , 8(4), 423-428. (in Chinese) [9] Torrence. C., Compo. G. P. 1998. A practical guide to wavelet analysis. Bulletin of the American Meteorological Society , 79(1), 61-78. [10] Xu, Q. X., Shi, G. Y., Chen, Z. F. 2004. Analysis of recent changing characteristics and tendency runoff and sediment transport in the upper reach of Yangtze River. Advances in water science , 15(4), 420-426.(in Chinese) [11] Xu. J. X., Sun. J. 2007. Sediment yield in major sediment source areas of the upper Changjiang River basin in response to human activities. Scientia Geographica Sinica , 27(2), 211-218. (in Chinese) [12] Xu. J. X., Sun. J. 2008. Temporal variation in suspended sediment concentration of the upper Changjiang River and its tributaries. Geographical Research , 27(2), 332-342. (in Chinese) [13] Xu. J. X. 2009. Recent variations in water and sediment in relation with reservoir construction in the upper Changjiang River Basin. Journal of Mountain Science , 27(4), 385-393. (in Chinese) [14] Zhang. L. L., Chen. J. 2007. Analysis on change of water flow and sediment in upstream of Yangtze River. Journal of Yangtze River Scientific Research Institute , 24(6), 34-37. (in Chinese) [15] Zhang. Q., Xu. C. Y., Becker. S., Jiang. T. 2006. Sediment and runoff changes in the Yangtze River basin during past 50 years. Journal of Hydrology , 331, 511-523. [16] Zhang. Q., Chen. G. Y., Su. B. D., Disse. M., Jiang. T,. Xu. C. Y. 2008. Periodicity of sediment load and runoff in the Yangtze River basin and possible impacts of climatic changes and human activities. Hydrological Sciences Journal , 53(2), 457-465.