A new approach was presented to calculate the total drainage area (TDA) of many small hydraulic engineering projects on the upstream Wangkuai reservoir watershed considering flood scaling. By establishing the correlation of flood peak between Wangkuai and Fuping hydrological station, the event-based flood scaling in Wangkuai reservoir watershed was analyzed and the scaling exponent θ was 0.5688. Relations between flood peak and maximum 3 h rainfall amount were fitted for the undisturbed and disturbed periods to analyze the effect of small reservoirs on the flood peak. The average change ratio of flood peak was 15.05% and the TDA of small hydraulic structures was 939.88 km2.
In recent years, anthropogenic activity, like construction of hydraulic engineering projects and urban expansion, has changed land surface conditions, and therefore caused environmental change which has altered runoff generation mechanisms. It has also resulted in the nonstationarity in hydrological series which invalidate traditional flood frequency analysis based on stationary assumption (Milly et al. 2008). The environmental change which results from anthropogenic activity has destroyed the stationary hydrologic regularity, and its effect on the hydrologic cycle is necessary in research. Some research has studied the impact of environmental change on the generation of flood (Archer & Newson 2002; Nirupama & Simonovic 2007; Moel & Aerts 2011).
Land use change is one of most typical forms of environmental change that influence flood. Brath et al. (2006) used a distributed hydrological model to simulate flood in Samoggia River, Italy for different kinds of land cover, and the result showed that flood with a small return period was obviously impacted by land use change. Naef et al. (2002) found that flood was impacted by land use change in areas with fast runoff generation and the flood flows reduced in these places. Cammerer et al. (2013) analyzed flood risk under simulated future land use until 2030 and found settlement expansion inevitably increased risk; and the change rate depended on the land use conversion type. Some researchers have studied the effect of change of land use, such as afforestation and urbanization, on flood peak and duration (Robinson 1980, 1993; Archer & Newson 2002). It was found that urbanization and deforestation could increase flood peak in the present scenario, while afforestation demonstrated the opposite role (Liu et al. 2005). Saghafian et al. (2008) analyzed land use and flood peak of 1967 and 1996 of Golestan Province in Iran. They found that the agricultural area significantly increased and the peak flow under different return periods increased. Thus, the effect of environmental change is important and should not be neglected in flood risk and management analysis.
However, the effect of hydraulic engineering is an important contributing factor which was not considered in the above researches. Hydraulic engineering can change the mechanism of flood generation and runoff process directly by management measures. Many hydraulic engineering projects have been constructed all over the world for soil-water conservation and water resource management, but they have also influenced regional flood generation significantly and changed the natural runoff process downstream (Zahar et al. 2008). The effect of some hydraulic engineering projects was difficult to derive because some small structures, like check dams, lack design data and therefore their effect on the flood process is hard to quantify. In this case, the effect of hydraulic engineering needs to be estimated for further analysis. Li et al. (2013) used hourly rainstorm and flood data of Lengkou station in Luanhe River basin, China, and found flood volume and peak decreased by 5.1%–16.6% and 0–57.3%, respectively, due to the construction of soil and water conservation engineering upstream. Yu et al. (2014) further estimated water storage capacity as 7.8 mm and drainage area of the small hydraulic structures upstream of Lengkou were estimated as 73% of the whole catchment area.
Around 6,000 check dams have been built in the Wangkuai reservoir catchment. Due to there being no long-term program, observational data was difficult to obtain (Li et al. 2016). Also the design information was unknown. Li et al. (2016) presented an approach to calculate the total storage capacities (TSC) and drainage areas (TDA) of many check dams in Wangkuai watershed based on rainstorm and flood events. The results showed that the TSC was 14.89 million m3, and TDA was 1,304 km2 which accounts for 34.6% of the whole area of Wangkuai watershed.
However, Li et al. (2016) did not take flood scaling into consideration in calculating the drainage area of the check dams. Flood in similar regions was usually considered to exhibit a power law which was also the basis of regional flood frequency analysis (Gupta et al. 1996; Medhi & Tripathi 2015). By developing regression relationships between flood peak and basin attributes, flood scaling was often used to estimate flood frequency in ungaged basins where observed data are deficient (Haddad et al. 2011; Han et al. 2012; Furey et al. 2016). Due to good reflection of drainage characteristic, drainage area is widely used to analyze flood scaling (Jothityangkoon & Sivapalan 2001; Furey & Gupta 2007; Ishak et al. 2011). Flood scaling can be expressed by a power law as , where Q means the flood peak, A is the catchment area, logα means the intercept parameter, and θ means scaling exponent. Because t θ is generally less than 1, it should not be neglected in estimating the effect of many hydraulic structures.
The aim of this study is to (1) analyze event-based flood scaling in Wangkuai watershed and (2) estimate the total drainage area of small hydraulic structures upstream of Wangkuai reservoir considering flood scaling.
STUDY AREA AND DATA
Wangkuai (WK) reservoir, located in Daqinghe River basin, was chosen as the study region (Figure 1). Wangkuai reservoir watershed is located at 114°29′E longitude and 38°45′N latitude. Building of the reservoir started in the year 1958 and was completed in the year 1960, with a storage capacity and drainage area of 1.39 billion m3 and 3,770 km2, respectively. It was built for the multiple functions of flood controlling, irrigation, and electricity generation. The watershed has an average precipitation of 626.4 mm per year, most of which is in July to September (70–80%). The average annual temperature is around 12.5 °C and annual pan evaporation is around 1,265 mm. Grassland and forest land are the main land use in the study region, derived from the remote sensed land use data of the year 2014 (over 90%). Fuping hydrologic station is on the upstream of Wangkuai reservoir (Figure 1). Fuping (FP) station was built in June 1958, with a drainage area of 2,210 km2 and annual mean precipitation of 600 mm.
Due to the serious soil erosion in Wangkuai watershed before the 1980s, a great deal of construction and soil-water conservation have been conducted from the 1980s, which have significantly changed the land surface of the watershed. Taking the Fuping hydrologic station as an example, there are 11 small reservoirs located in the control area with a total storage capacity of 5.9 million m3. They have a total control area of about 124.95 km2. Moreover, more than 6,000 check dams have been built since 1980, and these small-scale hydraulic structures are used for irrigation and soil-water conservation. These measures have influenced the flood mechanism in Wangkuai watershed, but the accurate area of these soil-water conservation projects cannot be derived.
In this study, Wangkuai reservoir watershed and Fuping hydrological station (located on the upstream of Wangkuai watershed) are selected to estimate the area of many small hydraulic structures in the Wangkuai reservoir watershed. Data of the Wangkuai reservoir watershed are from 1955 to 2008 and Fuping from 1958 to 2002.
Flood scaling analysis
Calculating the catchment area of many check dams
Detection of nonstationarity in peak discharge series
In order to analyze the impact of many small check dams on the flood in Wangkuai watershed, the nonstationarity of the annual maximum flood peak discharge (AMFP) series (Figure 2) of Wangkuai and Fuping was identified. Non-parametric Mann–Kendall test (Mann 1945; Kendall 1975) and Pettitt test (Pettitt 1979) were applied for detecting the nonstationarity of the AMFP series. The statistical values of Mann–Kendall test were −3.31 and −2.82 for WK and FP, respectively, exceeding the critical value (−1.96) at significance level of 0.05. This illustrated that AMFP series in the WK watershed showed a significant downward trend. The results of the Pettitt test are shown in Figure 3. Both WK and FP exhibited significant change point in 1979. There were a mass of small reservoirs and check dams built around the 1980s in WK watershed, so the most probable change point was deemed to be 1979. Some other studies also agree with the result (Gong et al. 2012; Li & Tan 2015; Deng et al. 2016). The nonstationarity showed that the many check dams in WK have influenced the flood obviously and, indeed, caused nonstationarity in AMFP series.
Flood scaling in Wangkuai watershed
Estimation of catchment area of many check dams in the study region
|Maximum 3-hr rainfall amount (mm) .||Qpre (m3/s) .||Qpost (m3/s) .||Change ratio of flood peak (%) .||Average change ratio (%) .|
|Maximum 3-hr rainfall amount (mm) .||Qpre (m3/s) .||Qpost (m3/s) .||Change ratio of flood peak (%) .||Average change ratio (%) .|
Previous research considered that flood peak and drainage area showed a simple linear relation (Yu et al. 2014; Li et al. 2016). However, flood in watersheds exhibits a scaling effect in many studies (Goodrich et al. 1997; Ogden & Dawdy 2003; Ayalew et al. 2014). Hence, the former methods are unreasonable and need to be improved upon. In this paper, we analyzed the event-based flood peak scaling in Wangkuai reservoir watershed based on the correlation of event flood peak between Wangkuai reservoir and Fuping station according to the actual flood data. The flood peak of two watersheds showed significant linear correlation and the scaling exponent θ of WK watershed was calculated as 0.5688. While event-based flood scaling should be analyzed based on actual data of more sub-watersheds in Wangkuai reservoir watershed, due to a lack of data we only used flood data of two watersheds and the result may be inaccurate. Flood scaling is affected by factors like rainfall spatio-temporal distribution and land surface condition. In this paper, the coefficient of determination for fitted scaling curve is around 0.6 and rainfall spatial distribution should be considered in future research.
Based on historical hydrological data, Yu et al. (2014) and Li et al. (2016) estimated the water storage capacity of soil-water conservation projects upstream of Lengkou watershed and Wangkuai watershed, respectively. Due to scaling exponent θ generally less than 1, the result should be larger in the same study area if considering flood scaling according to Equation (3). However, Li et al. (2016), who did not consider flood scaling in Wangkuai reservoir watershed, found that the drainage area of the small hydraulic structures upstream of Wangkuai reservoir was 1,304 km2, that result being larger than the result in this paper (939.88 km2). The main reason is that more rainfall-flood events were selected in the undisturbed and disturbed period in this paper and, indeed, the influence of environmental change on flood peak was 15.05% which is larger than the research of Li et al. (2016) whose result was 34.6%. The larger sample size reduced the sampling error and the result was more credible in this paper.
Hydraulic engineering has a significant function of regulation to floods. The effect of large hydraulic engineering projects is usually easy to determine because the design data and scheduling rule are generally available. However, it is difficult to quantify the effect of small structures without sufficient data, such as small reservoirs and check dams. By using rainfall and flood data of Wangkuai reservoir and Fuping station upstream, we presented a new method to calculate the control area of many check dams considering flood scaling in Wangkuai watershed.
The nonstationarity was detected based on the AMFP series of WK and FP and a significant downward trend was diagnosed in both series. Moreover, 1979 was detected as the change point of AMFP series. Thus, in this paper, pre-1979 and post-1979 periods were considered as the undisturbed and disturbed period, respectively.
By establishing the correlation between WK and FP, the event-based flood scaling in Wangkuai watershed was analyzed. The flood peak showed significant linear correlation between WK and FP (R2 = 0.7436), and event-based flood scaling exponent θ was calculated as 0.5688. Then, rainfall-flood events in disturbed and undisturbed periods were selected to analyze the effect of the many check dams on flood. The exponential function between peak discharge and maximum 3-hr rainfall amount pre- and post-1979 were fitted and the average change ratio of flood peak was 15.05%. Therefore, the total drainage area of all the check dams upstream of Wangkuai reservoir was estimated as 939.88 km2.