Abstract

Characteristics of annual and seasonal precipitation variation are explored in the upstream of Minjiang River (UMR), Southwestern China, spanning from 1960 to 2015. The moment of method (MOM), linear regression method, Mann–Kendall test, sequential cluster analysis, and Morlet wavelet analysis were utilized. The results clearly show the following: (1) Distribution of precipitation is uneven in space, with more in the south and less in the southeast. Decade average of annual precipitation reached the lowest in the 2000s and increased during 2010–2015 at all gauging stations and UMR. (2) Areal annual precipitation exhibited an insignificant decreasing trend with a rate of 4.47 mm/10a, which was mainly attributable to decreased summer precipitation. Spring precipitation exhibited an insignificant increasing trend and winter precipitation remained unchanged. (3) The change points mainly appeared in the 1980s and 1990s. And the almost periods of study area were generally 2–5 years, 7–11 years, and 15–20 years. (4) The increasing trend of annual precipitation is relatively obvious at higher altitudes, while the decreasing trend is more significant at low altitude stations.

1. Introduction

The issue of global climate change has become increasingly prominent, which is one of the world-wide problems that should be faced during economic and social development [13]. As is known to all, climate change has a relatively large influence on changes in meteorological factors such as temperature and precipitation. In recent years, the precipitation of China has changed abnormally under climate change condition [46]. Changes in the precipitation pattern may lead to droughts, floods, soil erosion, and loss of biodiversity and agricultural productivity. Hence, analyzing statistical properties of precipitation plays a crucial role not only in drought and flood risk assessments but also in water resources planning and management.

Trend analysis and abrupt change detection in hydrometeorological time series have been investigated using statistical methods by many researchers [715]. Sayemuzzaman et al. [16] utilized the Mann–Kendall test and the Theil–Sen approach to detect the trends of annual and seasonal precipitation sequences in time and space across the state of North Carolina, United States. Shirvani [17] detected the change point in the annual Persian Gulf sea surface temperature anomalies time series during 1951–2013 based on Student’s t parametric and Mann–Whitney nonparametric change point models. Zarenistanak et al. [18] analyzed trends and change points of annual and seasonal precipitation and temperature series over southwest Iran with Pettitt’s test, Mann–Kendall test spanning from 1950 to 2007.

The wavelet method has good time and frequency multiresolution, which can effectively diagnose the main frequency component and abstract local information of the time series. Some researchers have applied wavelet analysis to investigate the periods of time series such as runoff, precipitation, and temperature [1921]. Jemai et al. [22] revealed the periods of rain in Gabes Basin based on wavelet and coherence wavelet analyses. Sang et al. [23] used the improved continuous wavelet transform method to study the periodic characteristics of several typical hydrological series, including runoff and precipitation data.

Due to the impact of climate change, the frequency and intensity of droughts and floods in the southwest of China have been increasing [24, 25]. As one of the core areas of southwest China, the UMR, located in the transitional area between the Qinghai–Tibet Plateau and the Sichuan Basin, is the main source of water for Chengdu plain, and its runoff is mainly recharged from precipitation. However, frequent droughts and floods in this area have hindered its economic development and resulted in severe environmental problems [26]. Thus, the study of hydrological and meteorological variation characteristics in the UMR is increasingly being paid more and more attention. Huang et al. [27] detected the characteristics of hydroclimatic change in the UMR based on seven hydrological stations and six meteorological stations over 1956–2007, including precipitation, temperature, and runoff. The results indicated that there was a significant decreasing trend in precipitation, a significant increasing trend in annual average temperature, and a significant decreasing trend in total runoff. Zhou et al. [28] explored the spatial and temporal distribution of precipitation in Minjiang River Basin from 1956 to 2001 and revealed that precipitation in flood season presented a decreasing trend, and its periodic variations are 2 years, 5–7 years, and 14–15 years. Ma et al. [29] analyzed trends and abrupt changes of monthly and annual average temperature and precipitation of Songpan County. To the best of our knowledge, previous studies are almost before 2010 and do not cover the entirety of UMR as a study area to systematically and deeply detect precipitation variation.

The key objective of this study is to investigate the space-time distribution, trends, abrupt changes, and periods of precipitation in the UMR by using historical observed daily precipitation from 1960 to 2015 at annual and seasonal time scale. At the same time, the cause of precipitation change is also discussed. Such a study is beneficial for reducing the risk of disaster, reasonably allocating water resources, and supporting the decisions of water resource management under climate change condition. It also plays a crucial role in providing a reliable basis for government decision-making and helping the regional social and economic sustainable development. After Introduction, the study area and the selected data are given in Section 2. Section 3 lists the appropriate methodology. In Section 4, the obtained results are presented and discussed. Finally, conclusions are summarized in Section 5.

2. Study Area and Data

2.1. Study Area

As one of the major tributaries of the Yangtze River in China, Minjiang River originates in the southern part of Minshan Mountain on the border between Sichuan and Gansu provinces, crosses Chengdu Plain, and runs into the Yangtze River from Yibin city. The upstream of Minjiang River (UMR) refers to the river above Zipingpu Reservoir, which is a plateau climate and subtropical climate region. It stretches from longitude 102°45′E to 103°96′E and from latitude 30°80′N to 33°20′N. It covers an area of approximately 23000 km2 (16.9% of the entire Minjiang Basin), and the river course is 341 km long. The location of the study area is shown in Figure 1.

2.2. Data

There are five meteorological stations (Maoxian, Songpan, Wenchuan, Lixian, and Heishui) and one hydrological station (Zipingpu) available within the study basin, and they are distributed over the study area as shown in Figure 1. Detailed information regarding the stations is summarized in Table 1. Daily precipitation data of above 6 gauging stations were collected from the China Meteorological Data Sharing Service System and the Bureau of Hydrology and Water Resources of Sichuan Province from the period 1960 to 2015. The little missing data were interpolated by linear regression of the adjacent station data. Daily values were summed up to obtain seasonal and annual totals of precipitation. Areal precipitation was obtained by the Thiessen polygon method. The season in a year was defined as follows: winter (December, January, and February), spring (March, April, and May), summer (June, July, and August), and autumn (September, October, and November).

3. Methodology

3.1. Statistical Parameters Estimated Based on MOM

The statistical parameters, mean and deviation coefficient () of seasonal and annual precipitation, have been estimated with the moment of method (MOM) [30]. The formulas of MOM are not listed here.

3.2. Trend Detection Method
3.2.1. Linear Regression Method

The linear regression method is extensively used to determine the trend of climatic series [31]. In this study, unitary linear equation about precipitation to time is expressed aswhere is the time variable, is the observed precipitation value at time , is a constant, and is the slope which denotes the change rate. In detail, a positive value of indicates an increasing trend, and a negative value of indicates a decreasing trend. We used Student's t-test to examine the significance level of parameter , and the details are not listed here.

3.2.2. Mann–Kendall Test of Trend

The nonparametric Mann–Kendall (MK) test [32, 33] has been widely used to diagnose trends in time series. The statistic of MK test is defined aswhere and are the observed precipitation values at time i and j, respectively, is the length of the data, and sgn( ) is the sign function.

The variance of S and the test statistic Z are calculated as follows:

A positive value of indicates an increasing trend; a negative value of indicates a decreasing trend. If |Z| < Z1−α/2, the null hypothesis is accepted, and the trend is not significant. The null hypothesis is rejected at the α-significance level if |Z| ≥ Z1−α/2. When the α-significance level is 0.05, the corresponding value of Z1−α/2 is 1.96.

In the MK test of trend, another useful index is the Kendall slope, β, which is an unbiased estimator of monotonic trend magnitude:where 1 < < < . A positive slope value (β > 0) denotes an increasing trend; a negative slope value (β < 0) denotes a decreasing trend.

3.3. Test Method of Abrupt Change
3.3.1. Mann–Kendall Test of Abrupt Changes

A rank series is constructed for a time series x with n samples [34, 35]:where

can be considered as normal distribution under the null hypothesis that a trend does not occur, and the mean and variance of can be calculated as follows:

The statistic is computed by

In order to search an abrupt change, it is necessary to perform a similar analysis on the reverse time series. A retrograde is expressed as

The intersection point of and located between the confidence lines is the time when an abrupt change takes place. A typical confidence level of 95% is utilized in the detection of the precipitation series.

3.3.2. Sequential Cluster Analysis

Sequential cluster analysis is an effective method for estimating the abrupt change point τ in hydrological time series. The essence of sequential cluster analysis is to find the optimal split point [36, 37]. The objective function is calculated asin which, , .

When is the minimum value, its corresponding τ is the optimal split point, which is the most likely change point.

3.4. Morlet Wavelet Analysis

Wavelet analysis is employed to detect the periods of time series [37]. For time series f(t), the continuous wavelet transform is calculated aswhere is the wavelet transform coefficient, is the mother wavelet, is the complex conjugate functions of , a is the scale factor, and b is a time parameter. In this study, we have chosen Morlet wavelet function to be the mother wavelet. The Morlet wavelet function is expressed as

4. Results and Discussion

4.1. Variation of Basic Statistical Parameters

Basic statistical parameters (mean and Cv) of annual and seasonal (spring, summer, autumn, and winter) precipitation are obtained for 6 gauging stations and the UMR during the period of 56 years (Table 2).

On the annual basis, precipitation in the UMR is uneven spatially and increases from the northeast to southwest for both upstream and downstream (Figure 2(e)); the maximum value of mean precipitation is obtained at Zipingpu (1179 mm) in the south, and minimum value of mean precipitation is obtained at Maoxian (475 mm) in the southeast. The varies from 0.11 (Songpan) to 0.20 (Zipingpu) with average value 0.10 of the entire basin, which indicates that interannual variation is stable in the study area (Figure 3). Figure 4 presents the decadal variation of annual precipitation. The mean annual precipitation of UMR was higher in the 1960s, 1980s, 1990s, and 2010–2015 and lower in the 1970s and 2000s. Especially, in the 2000s, mean annual precipitation reached the lowest value, nearly 47.8 mm lower than 680 mm of the interannual level. During 2010–2015, it reached the highest value, 22.8 mm higher than the average value in the studied period, which is consistent with the corresponding studies [38]. Moreover, the decade average precipitation reached the lowest in the 2000s and increased during 2010–2015 at all gauging stations. A preliminary conclusion was drawn that the precipitation in the study area was relatively abnormal after the 21st century.

Space distribution of seasonal precipitation is uneven, with more in the south and less in the southeast (Figure 2(2)2(d)). Mean values of precipitation vary from 132 mm (Maoxian) to 224 m (Heishui) in spring, from 230 mm (Maoxian) to 647 mm (Zipingpu) in summer, from 103 mm (Maoxian) to 279 mm (Zipingpu) in autumn, and from 10 mm (Maoxian ands Wenchuan) to 49 mm (Zipingpu) in winter. We can conclude that the intraannual variation of annual precipitation is inhomogeneous with lager difference in four seasons; summer precipitation is the most and winter precipitation is the least. Cv ranges from 0.19 (Wenchuan, Lixian, and Heishui) to 0.25 (Zipingpu) in spring, 0.17 (Songpan) to 0.34 (Zipingpu) in summer, 0.22(Heishui) to 0.35 (Zipingpu) in autumn, and 0.33 (Zipingpu) to 0.52 (Wenchuan) in winter. It is obvious that the interannual variations of seasonal precipitation are different; the change is the largest in winter, following autumn, summer and spring (Figure 3).

4.2. Trend Analysis
4.2.1. Trends Based on Linear Regression Method

The b values of linear regression method are listed in Table 3. On the annual basis, the precipitation series of the entire basin showed an insignificant decreasing trend with a rate of −5.61  mm/10a, and large fluctuations occurred in the 1970s and 2000s (Figure 5(e)). For seasonal scale, an insignificant increasing trend was found in spring and winter precipitation series. In summer and autumn, the areal precipitation series had an insignificant decreasing trend.

4.2.2. Trends Based on MK Test

The MK test was utilized to determine the trend using two-tailed hypothesis at 5% significant levels. The MK test results (Z) along with the magnitude of trend (β) are presented in Table 3. Comparing with the linear regression method, their results are almost consistent. Box diagram of Z at 6 stations is displayed in Figure 6.

For annual scale, areal precipitation showed an insignificant decreasing trend, and the Kendall slope was −4.47 mm/10a, which was consistent with the previous studies [28, 38, 39]. A significant decreasing trend was only found at Zipingpu station. Analyzing the precipitation series on seasonal scale, areal precipitation displayed an insignificant increasing trend in spring (β = 3.45 mm/10a); an insignificant decreasing trend in summer (β = −8.00 mm/10a) and autumn (β = −3.60 mm/10a) [40]. In comparison with other seasons, winter precipitation was relatively stable with little annual variability at a rate of 0.35 mm/10a in the past 56 years. Comprehensive analysis leads to the conclusion that the decreased annual precipitation was mainly attributable to decreased summer precipitation in the upstream of Minjiang River [27]. There was a significant increasing trend of spring precipitation at Songpan station and a significant decreasing trend of summer precipitation at Maoxian, Wenchuan, and Zipingpu stations.

4.3. Abrupt Changes Based on Sequential Cluster Analysis

Detecting the precipitation change point is an essential step that can enable a better interpretation and more accurate forecasts of hydrological data. The change points of the annual and seasonal precipitation time series are presented in Table 4. Parts of diagnosed graphs are given in Figure 7.

For annual precipitation of UMR, there were several intersections between the UF and UB series in the 1980s and 1990s (Figure 7(e)); and the optimal split point was in 1995 (Figure 7(j)). Although the diagnosis results are not same, it is more appropriate to divide the year of abrupt change into around 1995. The mean precipitation was 693 mm during 1960–1995 and decreased to 657 mm during 1996–2015. Therefore, the results of MK test of abrupt change and sequential cluster analysis indicated that there was a sudden decrease in annual precipitation around 1995 [41]. As for spring precipitation of UMR, the intersection between the UF and UB series was in 2001, but the optimal split point was in 1995. There was a large difference of time-interval, so we believed that no abrupt change occurred. The abrupt change result of winter precipitation was the same as spring. By selecting the same point of two methods, the change points took place in 1994 during summer, in 1976 during autumn. In sum, the change points mainly appeared in the 1980s and 1990s, which was similar to the past studies [35].

There exists certain degree of abrupt changes of the annual and seasonal precipitation series; and change points show different characteristics among different gauging stations and UMR. Such spatial variability might be derived from the different natural topography, geographical location, and anthropic activities in the study area. And the specific reasons need to be further analyzed.

4.4. Periods Based on Morlet Wavelet Analysis

The periodic fluctuation of precipitation results in the up and down of runoff, so it is important to correctly understand the periodic variation of precipitation. Morlet wavelet analysis was used to identify periods of precipitation series. And the results are listed in Table 5. Time frequency graphs of wavelet transform in UMR are presented in Figure 8.

As shown in Figure 8(e), the time parameter b is the abscissa, and the frequency parameter a is the ordinate. The upper part is low frequency, and the isolines are relatively sparse, which corresponds to a long period of oscillation. The lower part is high frequency, and the isolines are relatively dense, which corresponds to a shorter scale cycle. It can be seen that areal annual precipitation had periods of 2, 5, and 17 years. The almost periods of seasonal precipitation were 2, 4, and 19 years in spring; 5 years in summer; 2, 5, 10, and 17 years in autumn; 3 and 15 years in winter. The detailed information of periods at gauging stations is no longer elaborated. In conclusion, the almost periods of study area were generally 2–5 years, 7–11 years, and 15–20 years [28].

4.5. Cause of Precipitation Change

Global warming is the background factor that influences the precipitation variation in the UMR. Over the twentieth century, especially since the 1990s, continuously increasing concentration of atmospheric carbon dioxide is resulting in global warming. Higher average air temperatures result in higher evaporation rates, higher water vapor contents, and consequently an accelerated hydrologic cycle. An important consequence of global warming is increased atmospheric moisture levels, thunderstorm activity, and/or large-scale storm activity. In the context of global warming, the precipitation pattern has changed in the UMR, which is influenced by atmospheric circulations, altitude, and anthropic activities.

The large-scale atmospheric circulation leads to the variation of precipitation over the UMR [40,42,43]. The advance and retreat of East Asian summer monsoon determine to a large extent the timing of the precipitation season and the amount of precipitation [44]. Other special monsoon-like circulations, such as the western Pacific subtropical high, Somali jet, and tropical cyclones, together with a blocking high and low trough in the mid-high latitudes, also play active roles in precipitation variation over the study area [45].

Altitude is another cause of precipitation variation [46]. The statistic Z value of annual precipitation has a good relationship with altitude, and the correlation coefficients are 0.99 (Figure 9). That is, the increasing trend of precipitation is relatively obvious at higher altitude, while the decreasing trend is more significant at low-altitude stations. On the whole, the Z values in UMR increased with altitude, which is in accordance with the previous studies [47].

Some considerations show that the influence of anthropic activities on precipitation is robust [48]. Anthropic activities not only modify the exchange of heat, water, and momentum between the land surface and overlying atmosphere but also change the composition of the atmosphere over urban areas. Land use change driven by anthropic activities influences the local hydrometeorological processes and sometimes affects the precipitation significantly.

The physics mechanisms of precipitation variation are a challenge to deeply understand. Further investigation of the multiscale interaction among the atmospheric circulation, altitude, and anthropic activities may also be important and essential to understand the precipitation variation.

5. Conclusion

In this study, space-time distribution characteristics, trends, abrupt changes, and periods of annual and seasonal precipitation series in the UMR are analyzed during 1960–2015. According to the above analysis, the main conclusions are summarized as follows.

Distribution of precipitation is uneven in space, with more in the south and less in the southeast. The interannual variation of precipitation showed a clear regional and seasonal difference in the UMR. The decade average precipitation reached the lowest in the 2000s at all gauging stations and UMR and increased during 2010–2015. A preliminary conclusion was drawn that the precipitation in the study area was relatively abnormal after the 21st century.

Areal annual precipitation exhibited an insignificant decreasing trend from 1960 to 2015, which was mainly attributable to decreased summer precipitation. Spring precipitation exhibited an insignificant increasing trend, and winter precipitation remained unchanged.

The change points mainly appeared in 1980s and 1990s. And the almost periods of study area were generally 2–5 years, 7–11 years, and 15–20 years. Due to the different natural topography, geographical location, and anthropic activities at the gauging stations, the different stations show different results of abrupt changes and periods.

Under the background of global warming, the precipitation variation of the UMR is influenced by atmospheric circulations, altitude, and anthropic activities. The increasing trend of annual precipitation is relatively obvious at higher altitudes, while the decreasing trend is more significant at low altitude stations. Further analysis is needed to quantitatively and systematically investigate the interactions among atmospheric circulations, altitude, and anthropic activities. Related research will provide scientific theoretical basis for water issues under climate change and help us to improve areal hazard mitigation and water resources management.

Data Availability

Most of the data series used to support the findings of this study is collected from the China Meteorological Data Sharing Service System, which is available on the web at http://data.cma.cn/. Only the data series at Zipingpu Station is supplied from the Bureau of Hydrology and Water Resources of Sichuan Province and cannot be made freely available.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Science Foundation Committee of China (No. 51679155).