1. Introduction
Traditional design flood for hydraulic structures such as reservoirs is based on a stationary hypothesis, meaning that the driving factors (e.g., climate change, urbanization and reservoir flood regulation) act to generate the hydrological variables in the same way as in the past, present and likely the future [
1,
2,
3,
4,
5]. However, the statistical characteristics of flood series might alter due to a changing environment [
6,
7]. If hydrological engineers do not fully consider the nonstationarity of the hydrological series, the results of the conventional stationary flood frequency analysis would be inaccurate in practice [
8]. López and Francés [
9] used climate and reservoir indices as external covariates in nonstationary flood frequency analysis. Yan et al. [
10] considered climate change and population growth to explain the nonstationary properties of hydrological time series. Global warming, the primary factor that drives climate change, has altered the meteorological regimes in some regions [
11,
12].
Population growth will not only lead urbanization but also lead to increased water consumption. Rapid urbanization over recent decades has significantly changed the regional hydrological characteristics of catchments [
13,
14,
15]. The river network systems have been obviously influenced by the process of urbanization, aggravating the hazards of floods and water degradation [
16,
17]. Water consumption is broadly reflected in daily life and productions, which can also affect the hydrological regime in the catchment. Reservoirs, one of the effective facilities for flood control, hydropower generation and other social functions, have gradually formed cascade reservoir systems [
18,
19]. The great impact brought by reservoir flood regulation has tremendously altered the hydrological characteristics in rivers [
20]. To sum up, it is difficult to find a river basin that is not impacted by global warming and anthropogenic activities, particularly in rapidly developing China.
Milly et al. have elucidated the challenges about how to deal with design floods and water resources management under nonstationary conditions [
1]. The commonly used technique to gain the changing knowledge of flood regimes is nonstationary flood frequency analysis [
9,
21]. Strupczewski et al. [
22,
23,
24] presented a nonstationary approach to at-site flood frequency analysis, in which the distribution parameters are represented as the functions of explanatory variables to explain the nonstationary hydrological series [
25,
26]. Stasinopoulos and Rigby proposed “Generalized Additive Models for Location, Scale and Shape” (GAMLSS), which is a powerful implement for nonstationary frequency analysis of time series [
27]. The meteorological factors (such as precipitation), population growth and reservoir index (RI) are widely used as explanatory variables incorporated in GAMLSS as covariates [
9,
28,
29]. Taking RI for an example, López and Francés [
9] considered the effects of RI from two aspects: the data used for flood frequency analysis are observed daily flow series’ that have been affected by upstream reservoirs; some of the reservoir characteristic parameters, such as the catchment area and the reservoir total storage capacity, are integrated to form RI. However, the RI series are piecewise constants in spite of reservoir operation rules [
30]. Although RI has been improved (by a hand of studies) for greater performance in nonstationary model fitting, the strategies of reservoir operation are hard to be quantified with the GAMLSS-based nonstationary flood frequency analysis framework [
8,
31]. Therefore, the covariates grounded in (or modified by) RI are unable to accurately consider the impact of cascade reservoir regulation.
The flood regional composition (FRC) combined with reservoir operation rules [
32,
33,
34] can overcome this drawback. The aim of the FRC is to study the flood generation mechanism at the investigated downstream site. The inflow of the target reservoir consists of the first upstream reservoir inflow and all interval inflows between adjacent reservoirs. Then, the outflow of the downstream reservoir can be obtained through the reservoir operation rules. Among all possible compositions based on the water balance equation, an appropriate FRC needs to be selected. Guo et al. [
33] proposed the most likely regional composition (MLRC) method and derived theoretical formula for triple cascade reservoirs. The MLRC method presumes that FRCs are diverse with their occurrence probabilities, which can be quantified by the multivariate probability density function (PDF) of flood events occurring at all sub-basins, and the FRC with the largest occurrence probability should be chosen for representing the actual spatial correlation of flood events. With a rigorous statistical basis, the MLRC method has been successfully applied in the cascade reservoirs in the upper Yangtze River [
35]. In contrast to the nonstationary flood frequency analysis applied with RI, the natural flow data restored by the observed data is employed in the FRC framework. Although the MLRC-based approach takes into full consideration both FRC and reservoir operation rules, it is based upon a stationary assumption condition so that the other variables such as climate change and population growth are ignored.
Under nonstationary conditions, the copulas with time-varying dependence structures have been primarily applied for modeling coincidence probabilities such as the joint return periods [
8,
36,
37,
38]. In fact, the changing environments might also alter the statistical correlation of FRC, and the dependence structure of the MLRC-based method is unable to catch the nonstationary spatial correlation. Therefore, we propose a nonstationary MLRC (NS-MLRC) method and then compare it with the MLRC method under stationary distribution conditions. Equivalent Reliability (ER) [
39] is employed in this study to concatenate the stationary and nonstationary design criteria grounded in a given return period [
40].
In a word, this study focuses on three objectives: (1) to propose a nonstationary design flood estimation framework; (2) to develop and verify the NS-MLRC method; and (3) to estimate design flood hydrographs at downstream sites. The rest of the paper is organized as follows:
Section 2 describes the methodology used in this study.
Section 3 briefly introduces the study area and data acquisition.
Section 4 analyzes the nonstationary design flood estimation results.
Section 5 discusses the nonstationary characteristics and the worst regional flood composition. Finally,
Section 6 ends with conclusions.
3. Study Area and Data
The Han River basin in China is located between 106–115° E and 30–35° N (see
Figure 3), and has a total length of 1530 km. As one of the most important tributaries of the Yangtze River, The Han River rises in the southern of Qinling mountains, and flows from the northwest to the southeast. This mountainous region lies in the humid zone with a subtropical monsoonal climate. The annual average temperature is between 14 and 16 °C. The annual precipitation varies from 700 to 1100 mm, and about 70 to 80% of the annual precipitation occurs during the flood season (May to October), in which heavy rains in early summer and continuous rainfall in autumn often cause major floods.
3.1. Cascade Reservoirs
The Ankang (AK) and Danjiangkou (DJK) cascade reservoirs are located at the upper and middle reach of the Han River basin, respectively. The AK reservoir was built in 1992 and provides hydropower generation and flood control, while the DJK reservoir was built in 1973 and its primary functions are flood control, water supply, hydropower generation, and irrigation. These two reservoirs are selected for the case study because their inflow data is lengthy enough (over 60 years) and they have large storage capacities.
The basic information of AK and DJK cascade reservoirs is listed in
Table 3. The characteristic parameter values and current flood control operation rules of the AK and DJK reservoirs are provided by the Changjiang (Yangtze River) Water Resources Commission (CWRC), Ministry of Water Resource [
70]. More information about AK and DJK reservoirs can be found in the references [
71,
72].
3.2. Dataset
Four categories of data series were collected, including restored streamflow data, observed hydro-climate data, population growth data and GCMs outputs from CMIP5. The observed data series provide information up to 2020, while a projected dataset of the future period from 2021 to 2095 is also used in this study.
(1) Restored mean daily streamflow data from both the inflow of AK reservoir and the inflow of DJK reservoir were provided by the Hydrology Bureau of the Changjiang (Yangtze River) Water Resources Commission during 1954–2020. The restoration of streamflow data can be taken as a natural flow series that eliminates the regulation impact of reservoirs.
(2) Observed daily precipitation series from 27 stations during 1951–2020 were obtained from the National Climate Center of the China Meteorological Administration (source:
http://data.cma.cn/ (accessed on 16 April 2021)).
(3) Given the unavailability of population data at the basin scale, the total registered population of all the prefecture-level cities amidst in the Han River basin was collected. These cities include Wuhan, Shiyan, Jingmen, Xiangyang in Hubei province, Hanzhong, Ankang and Shangluo in Shanxi province, and Nanyang in Henan province. The annual registered population data in Han River basin were obtained from the China Compendium of Statistics 1949–2008 [
73], the website of the National Bureau of Statistics of China (source:
http://www.stats.gov.cn/tjsj/ndsj/ (accessed on 1 July 2021)) and the websites of statistical bureaus of the provinces and cities mentioned above. For future projection, the logistic growth model that adapts the growth restriction resulting from limited natural resources is applied to predict the growth of population [
74]. Based on the logistic growth model, the evolution of the population for the period 1950–2100 is illustrated in
Figure 4.
(4) Future daily precipitation is simulated by GCM using climate change scenarios [
75]. Our lab research team, Tian et al. (2021) used 10 different GCMs (see
Table 4) and two representative concentration pathways (RCPs) of 4.5 and 8.5 from the IPCC Fifth Assessment to project climate change for the Han River basin, which are employed in this study [
76,
77,
78]. The RCP 4.5 scenario represents a future of medium emission where climate policies limit [
79]. Without climate change policies, RCP 8.5 scenario presumes that high emissions of greenhouse gases continue in the future.
The outputs of the GCMs not only involve the historical period before 2006 as a reference, but also cover 2021–2095 for future projection. A daily bias correction method is applied in this study for statistical downscaling [
80]. Six statistics containing mean, standard deviation, 85th, 90th, 95th and 99th percentiles of future precipitation series are used to test the performance of the daily bias correction method [
80]. Taking the BCC-CSM1-1 model as an example, the bias of raw and the corrected model for daily precipitation series during 1991–2005 are illustrated in
Figure 5, where the horizontal and vertical coordinates follow the 27 meteorological stations and the bias about the six statistics, respectively. These findings indicate the great performance of the statistical downscaling method so that it can be employed for future projection.
Figure 6 shows the projected future precipitation under RCP 4.5 and RCP 8.5 scenarios. The arithmetic mean values of the 10 GCMs are employed in this study for nonstationary flood frequency analysis.
6. Conclusions
There is an increasing need to develop an effective design flood estimation framework to deal with nonstationary data series caused by climate change and anthropogenic activities. In this study, the univariate flood frequency analysis, nonstationary hydrologic design criteria and the NS-MLRC method were adopted to derive the nonstationary design flood volumes in the Han River. The design flood hydrographs at the downstream site were estimated after cascade reservoir regulation. The main conclusions are summarized as follows:
(1) The proposed NS-MLRC method can be effectively implemented as an extension of the MLRC method for explaining the nonstationary spatial correlation of the flood events. The multiple nonstationary driving forces, i.e., climate change and population growth, can be captured and precisely quantified by the proposed design flood estimation framework.
(2) The slow-to-change impacts of climate change and population growth are presented in design flood volumes according to the nonstationary flood frequency analysis method and ER criteria. The long-lasting driving factors imply the larger risks of the flood hazard. The 1000-year design W15 of the DJK reservoir under the stationary distribution scenario (S1), RCP 4.5-based (S2), and the RCP 8.5-based (S3) nonstationary scenario are 20.041, 21.680 and 22.352 billion m3, respectively.
(3) The swift effects of cascade reservoirs are reflected in design flood hydrographs with lower peaks and less volumes based on the NS-MLRC method and flood control operation rules. For instance, the 1000-year design Qmax of the AK (DJK) downstream site under the stationary distribution scenario (S1), RCP 4.5-based (S2), and RCP 8.5-based (S3) nonstationary scenario decrease by 10.97% (49.43%), 5.06% (53.25%) and 5.31% (54.66%), respectively due to cascade reservoir regulation.
As the dominant element affecting flood design, the current cascade reservoir operation strategies can be improved to accommodate the slow-to-change nonstationarity in further research.