Abstract

In the arid areas, the discharge in some reaches is unsurveyed, which makes related works trapped. It is important to estimate a long-term discharge in an ungauged river using available fragmentary data including those on human activities. In this work, the authors proposed a procedure of monthly discharge estimation for an ungauged inland river based on spatiotemporally discontinuous information, such as water diversion for irrigation, water allocation planning, and statistical relations of discharge at different locations. A typical ungauged river reach of Heihe River in northwestern China was selected as a case study. Monthly discharge from 1991 to 2007 at one gauged section was estimated using this procedure to validate the procedure, and the results show its potential in estimating discharge for an ungauged area. Monthly discharge from 1979 to 2014 at four unmonitored sections was estimated using the same procedure, and the trend analysis was conducted to check the facticity of the estimated discharge. The annual and monthly trends could well describe the effects of anthropogenic activities on the flow at different sections. In general, the method is feasible to provide support pertaining to water management or eco-hydrological researches and could be easily extended to other similar places.

1. Introduction

In most arid and semiarid regions, river laterally recharges the riparian groundwater to sustain the terrestrial ecosystem along the natural channel and impact interactions between food, water, and energy in oasis [13]. Discharge data especially in arid and semiarid regions is quite important for water resource management and related scientific researches. However, insufficient data or the lack of information renders many tasks difficult [4]. A recent study has focused on ecohydrological effects of stream-aquifer water interaction in an inland river basin [5]. But due to lacking of discharge data, the research of long-term riparian ecohydrological effects is trapped. Many hydrologists have tried to address this issue using simple regression methods to demonstrate the value of providing a rough control on water flow in regions characterized by scarce data [68]. A nonlinear regression model of a hydrological system could obtain rational outcome without requiring additional information, such as evaporation or soil moisture [9]. More complicated statistical models built on seasonal-based grouped catchments bring more robust results [1012]. Physical models of gauged and ungauged basins have been studied more thoroughly and include regionalizing model parameters, constraining hydrologic simulations by regionalized signatures [1316] and assigning parameters from similar hydrological catchments [17, 18]. These models are principally hampered by unsound calibration [19]. Although the transferability of parameters estimation methods from small subcatchments to a large area performs well in some case, the common contention of comparing physical and functional similarities may be invalid for many catchments [14, 20]. Diversity in climatic variability and geographical attributes causes different performances for each regionalization method in various regions [21]. To overcome this deficiency, more complicated models with adequate spatial and temporal information have been used, such as CLM, MIKE SHE, GSFLOW, and CATHY [2225]. Although these process-based models can better describe the exchange of surface water and groundwater, the high computational cost makes it difficult to apply to large basins [26]. Long-term simulations by these models depend on spatiotemporally continuous forcing, which is however difficult to meet, especially in ungauged areas.

The Heihe River Basin (HRB), the second largest inland river basin in China, is located in a transitional zone straddling semiarid and arid areas. There is an obvious interaction between human activity, such as water use, and basin ecosystem [27], leading to substantial research focused on hydrological processes [2832] and water resources management [3335]. The discharge is an important indicator of the availability of water for the ecosystem in the middle and lower reaches of HRB [36]. However, there is only one hydrologic station (i.e., the Gaoya hydrologic station) between the inlet (i.e., the Yingluoxia hydrologic station) and the outlet (i.e., the Zhengyixia hydrologic station) located in the midstream [37] generating a typical discharge-ungauged problem for this river reach [38]. It is imperative to estimate the discharge at ungauged sections for the past decades to strengthen the water management. Previous studies have produced several useful methods and data; few of them however have explicitly addressed the issue of the long-term discharge estimation at different ungauged locations in the reach.

As part of a major initiative to provide a complete scientifically based planning for sustainable water management in arid and semiarid ecosystems, a monthly discharge estimation procedure based on spatiotemporally discontinuous data, was proposed using diversion for irrigation water, water allocation planning, water budget of the inland river, and statistical relations of discharge at different locations. In this work, monthly discharge from 1979 to 2014 at four unmonitored river sections and that from 1991 to 2007 at one gauged section were estimated using this procedure. The HRB is taken as a case study to propose an appropriate procedure for discharge estimation in ungauged rivers and provide data for long-term discharge to support water resource management and further ecohydrological studies in HRB.

2. Methodology

For ungauged river sections in arid and semiarid basins, we established a procedure of monthly discharge estimation with spatiotemporally discontinuous data information including those on human activities. Under natural circumstance, the main equilibrium of the river water volume is determined by the inflow, outflow, precipitation, evaporation, river leakage, and recharge from groundwater. For a river flowing in an arid area, the inflow from the upper mainstream principally dictates the total river water quantity of the middle and lower basin reaches. Thus, knowing the accurate inflow will provide the upper limit of water flow. Because of complicated hydrogeological distribution in arid and semiarid zones, factors such as river leakage and recharge from groundwater become more important than those in humid regions [39]. The lack of reliable data at the river scale produces uncertainties in determining different components. Because of limited data on precipitation and evaporation, these are given as a fixed seasonal change. Moreover, a large number of the population relies on the HRB for their subsistence so the effect of human disturbance on streamflow cannot be ignored [40]. Different demands for water diversion should be considered for irrigation and nonirrigation periods. The water diversion activity is intensive on some river reaches whereas there is no obvious human influence on others. Intense human activity in oasis segments occurring in the middle reach can severely restrict sustainable eco-development in the lower reaches of the HRB [41]. To ensure a sufficient water delivery to downstream ecosystems, the government instances would usually take action to limit the use of water [42]. This operation can affect the natural discharge along the river. However, it is difficult to obtain the detailed water allocation data. Hence, we need to incorporate a suitable parameter into the water balance equation to reflect the effect of governmental management on the quantity of river drainage. A water balance framework for river flowing in arid and semiarid regions is given in Figure 1. Because the monthly discharge estimation procedure is based on spatiotemporally discontinuous data, any term of the water balance could be calculated from any possible data. For example, if there are surveyed data or known statistical relations, they can be directly used for calculation. Otherwise, the indirect relations between any terms of the water balance at different or same locations should be established on estimated data at previous step. Main equation is shown below: means the monthly discharge at the section of outlet along the river reach; denotes recharge from precipitation; is the discharge at the section of inlet along the river reach; is the water extracted from river for irrigation by human; is the evaporation from the upper section to the lower section; means the distance between the two locations; is the water leakage rates from riverbed or recharge rate from groundwater per distance, which is calculated by statistical formulas; is the change of river water storage in each month. Water balance law is suitable for all segments of study area although and could be ignored in some river reaches without significant human activities and exchange of surface and subsurface water. Discharge estimation in this paper is primarily based on the equation and other statistical relations of discharge between different river sections. It means some terms, such as , would be estimated by statistical relations investigated from previous studies by others. Moreover, if the direct relation of discharge between two sites exists, it would also be used for calculation. Then the results calculated, respectively, from water balance equations and statistical relations should be averaged to be a robust result.

3. Numerical Simulation

3.1. Study Area

The HRB covers a total area of 143,000 km2 and is the second largest inland river basin in China, located at 98°E–102°E and 37°N–42°N. Water flowing from the headwater to downstream controls the landscape. Glaciers and permafrost are present in the high-mountain areas while forests and alpine meadows cover the middle-mountain regions. A desert-oasis-river landscape occurs in the middle and lower basin reaches, with well-developed agriculture in the midstream and a riparian ecosystem distributed downstream along the river. In the upstream basin, the mean annual precipitation is approximately 300–700 mm with distinct vertical zonal divisions. The upper reaches of the HRB contain nearly the entire water resource of the whole basin [43]. In the middle reaches, the elevation decreases from 2000 m to 1000 m and the precipitation correspondingly decreases from south to north from 200 mm to 100 mm [44]. The lower reaches having an average altitude of 1000 m constitute an arid region with a mean annual precipitation of 42 mm according to statistical data gathered from meteorological stations and simulations from regional climate models [45]. Conflicts between districts related to specific types of water demands have long prevailed in the region. The national water allocation plan effectively solved the conflicts [46] by closing all channels so that enough water can arrive at the lower reaches.

In this work, four ungauged river cross-sections were selected, namely, the Caotanzhuang site, the S213 Bridge, the G312 Bridge, and the Tielu Bridge. The Gaoya station and Zhengyixia station, both of which are the lower sites compared to Tielu Bridge in the middle reach, were selected to validate the procedure. All are located in the midstream area and shown in Figure 2. The middle reach experiences different hydrogeological conditions. From the Yingluoxia Gorge to the location near the G312 Bridge, the river level is higher than groundwater level so that groundwater recharge occurs through leakage of the riverbed [35]. Along this river reach, there is intensive human activity leading to water diversion from the Yingluoxia Gorge to the Caotanzhuang, whereas there is little water diversion in other parts. The river is again affected by water diversion from the G312 Bridge to the Tielu Bridge. The groundwater recharges the river from the G312 Bridge to the outlet (Zhengyixia Gorge) in the midstream [35, 47]. Then diversion for irrigation water reappears in some river reaches.

3.2. Data

Some important data, such as measured discharge at Yingluoxia, Gaoya, and Zhengyixia, were collected from the Hydrological Bureau of the Gansu Province, the Heihe River Bureau [48, 49]. Other data, such as estimated evaporation and monthly water diversion from river, were extracted mainly from the literature [29, 50]. Note that replenishment by precipitation is sometimes ignored because of the scant rainfall in the entire middle reaches [45]. Data would be introduced in detail in the following.

3.3. Estimation of Monthly Discharge at Ungauged River Sections

Downstream calculations, shown in Figure 3, were used to estimate the monthly discharge. The mainstream was divided into four parts. Part 1 runs from the Yingluoxia Gorge to the Caotanzhuang section (~10 km), Part 2 spans the Caotanzhuang site to the S213 Bridge (~12 km), Part 3 starts from the S213 Bridge and ends at the G312 Bridge (~12 km), and Part 4 goes from the G312 Bridge to the Tielu Bridge (~3 km). We initially used the monitored 1979 to 2010 monthly discharge, the 2011 to 2014 simulated monthly discharge, and the observed annual discharge [48, 49] at the Yingluoxia Gorge to generate a prolonged and precise series of monthly discharges. We then estimated monthly discharges at the Caotanzhuang, the S213 Bridge, the G312 Bridge, and the Tielu Bridge based on statistical and water balance methods that included water diversion for irrigation, river leakage, evaporation, and the national water allocation plan.

As described above, there is intense human activity related to water diversion from Yingluoxia to Caotanzhuang while there is nearly no water diversion in other parts of the river reach from Yingluoxia to Tielu Bridge. Accurate discharge estimation from the Yingluoxia Gorge to the Caotanzhuang hydrojunction will dictate the reliability of all discharge estimates downstream due to the water diversion caused by human activity. The water diversion related to irrigation and the national water allocation plan will also substantially modify the total discharge. Therefore, we adopted a statistical regression method and a water balance equation to calculate the discharge and then averaged the two results to reduce the uncertainty related to the methodologies. In the statistical regression method, we applied different calculations during the irrigation periods and the nonirrigation periods of January, February, and December. For the nonirrigation periods, a linear regression, shown as (2), was defined by the measured monthly discharge from the Yingluoxia and Caotanzhuang sites during the closed channels periods of 2011 and 2012 [51]. The unit is m3/s. The coefficient of determination is larger than 0.90 and it could be applied in January, February, and December. Another empirical equation, shown as (3), was found by the measured discharge from the Yingluoxia and Caotanzhuang sites during 1998 to 2001 [29] for the irrigation period. The unit is 108 m3/month. The coefficient of determination is about 0.90 and it could be used during other months. Due to lacking of monitoring, these equations from literatures could be deemed to be relatively believable in this paper. denotes the monthly discharge at the Yingluoxia section; means the monthly discharge at the Caotanzhuang site. Because the statistical computation was performed for the transitional period tentatively implementing the national water allocation plan, it is necessary to tune the results. To satisfy this constraint, we recalculated the discharge at the Zhengyixia site when it reached 15.8 × 108 m3 at the Yingluoxia location. We then calculated the ratios of the computed annual discharge of the 1979 to 1997 and 2002 to 2014 periods to the average annual discharge from 1998 to 2001 from the Zhengyixia site. The specific equations are shown in (4). Monthly discharge before 1998 was then corrected by the obtained ratios. The data after 2001 period were corrected using these determined ratios only if current monthly discharge given at Yingluoxia location grew more than 5% relative to the sum of the multiyear averaged monthly discharge and the standard deviation both of which were calculated by data from 1990 to 1999. The value of 5% is an assumed threshold by experience. indicates the revised monthly discharge at the Caotanzhuan location; stands for the correction ratio; signifies the number of months when the current monthly discharge at the Yingluoxia is 5% higher than the sum of the monthly mean and the standard deviation. To ensure the reliability of the results, we need check the computations in a different way. To get a more credible discharge at the Caotanzhuang site, we substituted the river seepage equations, shown as (5) and (6) with the coefficient of determination of 0.89 and 0.98, respectively, and provided by literatures [28, 29] and the 1990 to 2007 monthly extracted water data [50] of five irrigated areas (Daman, Yingke, Xigan, Gansuo, and Shangsan) into the river water balance shown as (7). The final river leakage is obtained by averaging the calculated results by two river leakage descriptions.The left variables of (5) and (6) are the leakage rates per distance. The units of and are (%  ×  km−1) and (km−1), respectively. Equation (7) is the river water balance model in study area. The term in (1) is ignored because it is small in river reach with a very short distance for one month. is set to a constant extracted from the literature. denotes the irrigation water diversion. Using these equations, we obtain the monthly discharge from 1990 to 2007 at the Caotanzhuang site based on the water balance method. Estimated discharges from 1990 to 2007 based on statistical regression and the water balance method applied to the Caotanzhuang site are shown in Figure 4. The estimated results from the statistical regression and water balance methods are consistent. The calculated Nash-Sutcliffe coefficient (NSE) [52] is 0.86, the correlation coefficient is 0.95, and the Root Mean Square Error (RMSE) is 11.38 m3/s, which shows the results provided by the statistical regression, and water balance methods are reasonable and comparable. The regression equations were applied using the monthly discharge from 1990 to 1999 and from 2000 to 2007 to extend the results to the range from 1979 to 2014. As shown in Figures 4(b) and 4(c), the linear regression equations are listed with coefficient of determination of 0.93 and 0.90, respectively, for the 1990–1999 and 2000–2007 periods. Both regression equations passed the significance test. Estimated monthly discharge for the 1979–1989 and 2008–2014 periods by statistical modeling was applied to complement the discharge obtained by the water balance method and linear regression equations. The two calculated results were then averaged to get a more reliable data. Once the monthly discharge at the Caotanzhuang site was determined, the discharge at the S213 Bridge and G312 Bridge river sections was calculated in the similar way. The main equations used for these river reaches are listed as follows:The in and out mark the location of the inlet and outlet of the river reach from Caotanzhuang to S213 Bridge or that from S213 Bridge to G312 Bridge. Because there is hardly water diversion from the Caotanzhuang site to the G312 Bridge, (10) omits . The river reach from the G312 Bridge to the Zhengyixia location is a water gaining segment, but there is no detailed information on the segment from the G312 Bridge to the Tielu Bridge. The distance between the G312 Bridge and the Tielu Bridge is only 3 km; we ignored both the river leakage and the recharge from groundwater and integrated the water extracted for Wujiang irrigation district. The estimated monthly discharge at four ungauged sections is shown in Figure 5.

3.4. Validation of the Procedure and the Estimated Discharge

The proposed procedure and the estimated discharge were validated by two steps. Firstly, the observed monthly data at Zhengyixia station, which belongs to the segment lower than that from Yingluoxia to Tielu Bridge, was used to validate the procedure. Then the trend analysis of estimated annual and monthly discharge was implemented to indirectly check the facticity of the estimated data.

The proposed procedure was used to estimate the monthly discharge at the outlet (Zhengyixia station) of the middle reach based on the limited data of the segment from Gaoya to Zhengyixia. Difference from the river reaches mentioned in Section 3.3 is that the groundwater chiefly recharges the river in this part. Recharge from groundwater was estimated based on the literature [29]. The total water extracted from river was computed based on shared data [50]. Then the input data from 1991 to 2000, such as precipitation, evaporation, observed discharge at Gaoya station, recharge from groundwater, and water diversion for irrigation, were collected and processed for water balance computing. Similar to the method used for the segment from Yingluoxia to Caotanzhuang, a linear regression relation was built on the observed data from 1998 to 2001 at Gaoya and Zhengyixia stations. The average of the data calculated by the two ways was regarded as the final result. The final estimated data from 1991 to 1997 and from 2002 to 2007 were compared with the observation at Zhengyixia station to validate the proposed procedure. The simulated result is shown in Figure 6. The results of NSE in two periods for validation are, respectively, about 0.70, which is satisfactory. Although the simulation does not match well with the observation during some seasons, the main change of monthly discharge could be captured. The gap between simulation and observation may be mainly attributed to three aspects. The first is the uncertainty of data related to human activities. The data just describe the total water extracted from river but ignore the distribution at cultivated land so that the lag of backing to river by overland flow may be difficult to be considered. The second aspect is the recharge from groundwater. Because the distance from Gaoya to Zhengyixia is about four times as long as that from Yingluoxia to Tielu Bridge, there may be more uncertainties of calculating recharge rate along the river. The third aspect is the regression relation based on linear assumption, which is however not in accord with the true world. In general, according to this validation, for the ungauged inland river in an arid and semiarid area, the design idea of the discharge estimation is reliable.

Based on regression analysis, presented in Figure 7, an increasing trend of discharge from 1979 to 2014 is observed. The amount of water at the Yingluoxia Gorge controls the upper limit of the downstream sectional flow. The annual rates of increase at the Yingluoxia, the Caotanzhuang hydrojunction, the S213 Bridge, the G312 Bridge, and the Tielu Bridge are 0.0695 × 108 m3/a, 0.0603 × 108 m3/a, 0.0549 × 108 m3/a, 0.0465 × 108 m3/a, and 0.0436 × 108 m3/a, respectively. Declining rates are in turn observed especially in the segments from the Yingluoxia to Caotanzhuang locations and from the S213 Bridge to the G312 Bridge. This may reflect an important water diversion and river leakage. This result is in agreement with the mainstream exchange flux variation determined by Tian et al. (2015) [35]. The result of the sequential Mann-Kendall test described in the literature [53] is shown in Figure 8. Both UF and UB in the figure are sequential statistics. If the value of UF is greater than 0, there is a rising trend and vice versa. If it is greater than 1.96 or less than −1.96, the change trend is statistically significant. The intersection of UF and UB denotes a change point over the time series. From the figure, it indicates a rising trend of annual runoff at the Yingluoxia location from 1979 to 1990, a downward trend from 1991 to 2006 followed by an increasing trend from 2007 to 2014. Annual runoff at other sections consistently displays a downward trend from 1985 to 2007 followed by an increase after 2007. The difference in runoff among ungauged sections is not apparent. A noticeable variation occurs in 2006 at the Yingluoxia site while those at other sites lag behind. This may be explained by climatic factors controlling the discharge at the Yingluoxia location whereas the discharge at midstream is affected by human activity and geological features. We further analyzed the variation in multiyear averaged monthly discharge by dividing the whole period into three periods. The river was less affected by human activity during the period of 1979–1989. Intensive effect of human activity existed in the period of 1990–2000. The national water allocation plan was carried out during the 2001–2014 period. Figure 9 shows the discharge at all sections exhibiting significant seasonal changes. The discharge at the Yingluoxia site varies slightly for the three periods. There is a peak discharge of 140 m3/s consistently occurring in July. During the high flow periods of June, July, August, and September, the ungauged sections discharge varies between 60 and 110 m3/s in the first period, declines from 40 to 80 m3/s in the second, and rises to a higher range in the last period. At Gaoya station, the monthly discharge is larger than that at Yingluoxia station in winter and spring. This is perhaps due to the fact that the groundwater recharges the streamflow around that place. Since 2000, the national water allocation plan was gradually implemented to generate an amount of ecologically satisfactory water level in the lower reaches of the HRB. On the whole, although the proposed monthly discharge estimation method could not fairly accurately calculate the discharge, the estimated streamflow at ungauged river could be appropriate and satisfactory.

4. Conclusions and Discussion

In the arid and semiarid areas, streamflow plays an important role in riparian water-energy-food nexus under frequent human activities. However, the unsurveyed discharge data may make related works trapped. It is important to take full advantage of fragmentary data to estimate a long-term discharge in an ungauged river reach. In this work, a procedure of monthly discharge estimation for an ungauged inland river was proposed based on spatiotemporally discontinuous data information including those on human activities. Monthly discharge from 1979 to 2014 at four unmonitored river sections (the Caotanzhuang, the S213 Bridge, the G312 Bridge, and the Tielu Bridge) and that from 1991 to 2007 at one gauged section (Zhengyixia) were estimated using this procedure. The observed data at Zhengyixia was used to validate the procedure. The Nash-Sutcliffe coefficient and Root Mean Square Error are, respectively, about 0.70 and 13.00 m3/s, which shows its potential in estimating discharge for an ungauged area. The trend analysis for discharge at four ungauged river sections was conducted to indirectly check the facticity of the estimated data. From annual trends, the calculated data probably reflects surface and subsurface water interaction. The monthly trends show that the national water allocation plan improved the ecological water supply to the lower reaches of the Heihe River Basin. Annual and monthly trend could well describe the effects of anthropogenic activities on the flow at different sections. In summary, the proposed discharge estimation procedure for ungauged rivers is reliable and this work also provided the long-term discharge data for different ungauged river sections to support water resource management and scientific researches of riparian water-energy-food nexus in HRB. But even more importantly, this procedure could be easily extended to other similar ungauged inland rivers.

Most of available data were employed to estimate the discharge in the study area. However, those data did not include factors such as Land Use/Cover Change [54] and the economic policy [55]. Moreover, a long-term continuous flow monitoring at ungauged sections [36] should be set up to calibrate and validate our method or other hydrological models. A recent airborne thermal infrared remote sensing study indicates the possibility of quantifying the groundwater discharge of a river [56], which may be useful for discharge estimation.

Additional Points

The estimated data can be obtained upon request from the authors and would be permanently archived at the Cold and Arid Regions Science Data Center at Lanzhou, http://westdc.westgis.ac.cn/data/9fb70559-fe79-4f99-ba11-3f02cdd184b0 [57].

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grants 41575096 and 91125016) and by the Key Research Program of Frontier Sciences, CAS (Grants QYZDY-SSW-DQC012).