A Framework to Regionalize Flow Information in a Catchment with Limited Hydrological Data ()
1. Introduction
The size and structure of watersheds as well as their topographic and geographic structures vary significantly due to geological, morphological, vegetational, soil, and climatic differences [1]. In response to human development, fast climate dynamics, and slow geologic processes, there has been a reciprocal evolution within catchments resulting in changes in the soils, vegetation, and topography, mediated by material and energy fluxes. This is exhibited in the catchment hydrologic partitioning at the landscape scale and the resulting runoff [2]. Anthropogenic activities by humans often (directly or indirectly) play an important role in altering these landscape characteristics [3] [4]. Quantitative or qualitative flow information is vital for many practical applications. These include water allocation, long-term planning, watershed management operations, flood forecasting, optimization of hydroelectric power production, and hydraulic structure design [5] [6].
Continuous observation of water discharge using the floats or the current meter is impractical, especially on small rivers due to time and resource constraints [7]. Subsequently, the determination of discharge by use of stage height and the stage-discharge rating curve is widely applicable to major rivers that are susceptible to floods or provide vital economic services to the larger ecosystem. In the absence of hydrological data, different approaches have been adopted to try and estimate hydrological aspects of the ungauged watershed using information from the gauged one(s) [8] [9] [10] [11]. One such method is the basic relationship regarding peak discharge at the ungauged site, whereby the gauged data are weighted as a ratio of drainage area with a slope exponent of the curve relating discharge to the watershed area for suitable gauges in the hydrological region [12]. Another technique is the use of a linear correction factor that accounts for the differences in the drainage areas between the gauged and ungauged sites [13]. The flow estimate for the ungauged site is determined by multiplying the correction factor for the ungauged site by the regional regression estimate for the gauged site. These could be applied to the drainage area of the ungauged site that has similar watershed characteristics and area ratios of 0.5 to 1.5. The use of the regional flow duration curve (FDC) has also been employed in the calibration of the ungauged sites [14]. However, the limitation of this method is that the absence of intangible flow timing information affects predictability in ungauged catchments [15].
Since the watershed differs appreciably in weather conditions, topography, vegetative cover, and geology, other researchers have utilized the spatial proximity and regression-based approaches to predict the flow characteristics of a river basin [8] [16] [17] [18]. The accuracy and validity of such flow estimates are directly hinged on the similarity of watershed characteristics (precipitation, drainage area and shape, orographic expression, aspect, vegetation, and lithology/geology) [19] [20]. The estimates at the ungauged sites can be traced back to the origin of the design flow at the ungauged site. Some researchers have determined discharge at the ungauged site by weighting the gauged data by the ratio of the drainage areas and multiplying by an exponent [19]. The exponent was the power function relating discharge to the watershed area (determined by regional regression of the flood region). Other researchers, [14] [21] recommended a different method using a linear correction factor for the difference in the drainage areas between the gauged and ungauged sites. In hydrologic regions [14] numerous gauges, regression equations have been developed and applied to extrapolate hydrologic information to the ungauged sites within such watersheds [22] [23].
The rainfall-runoff process is a complex, dynamic, and nonlinear process that is governed by multiple, often interrelated physicochemical factors [24]. The hydrologic model approach is widely used to quantify and qualify continuous flow of the rivers in response to human activities in the watershed [6] [25] [26]. Hydrologic models can provide information about rainfall-runoff processes in watersheds, although they remain abstractions of a real system. Any of these models cannot be assumed to exclusively provide accurate information about specific river basins and hydrologic conditions [27]. There are uncertainties in the direct implementation of hydrologic models at the ungauged catchments. One of the challenges is to determine whether the selected model structure reliably and adequately represents the hydrologic processes at the ungauged catchment of interest [28].
Rainfall-runoff models and hydrological model-independent methods have played a huge role in the regionalization of continuous stream flow [29]. Model parameters are used as a means of transferring hydrologic information from gauged to the ungauged basins, rather than a direct transfer of stream flow data through portioning methods [6]. In a process-based water quality model, such values may be numerous, and the problem of equifinality may arise, i.e., many model setups may yield similar model performance [30]. Due to heterogeneity, it is always very difficult to link all model parameters to the correct watershed characteristics [31]. The success of transferring hydrologic information between watersheds could be influenced by the spatial proximity and physical similarity between the donor (gauged) and recipient (ungauged) watersheds [32]. In other cases, regression equations that incorporate different gauges in a watershed with similar hydrologic characteristics are useful methods for extrapolating calibration parameters to ungauged sites [14] [29].
The SWAT (Soil and Water Assessment Tool) is an open-source distributed model [33] [34] that has a large and growing application in the hydrological characterizing of catchments [6]. Regarding previous studies, the tool has shown strength in discharge estimation in mountainous catchments [1] [35]. The model has successfully been used to estimate flow rates on a wide scale and in various hydrological regions. It was applied to evaluate the uncertainty and regionalization of hydrological information in numerous catchments [8] [9] [36] [37]. It is important to note that most of the above research involved the use of numerous existing gauges within a given watershed or pre-existing hydrologic regression equations.
The review report (the Decade of Prediction in the Ungauged Basins initiative (PUB 2003-2013)) indicates that the inability to make reliable predictions in ungauged watersheds will continue to impede sustainable water resources management and the development of effective flood and drought mitigation strategies, particularly in developing countries [31]. In the context of one or two gauges in a catchment with numerous river basins, the question is, “is calibrating a hydrological model using a single gauge’s data points an optimal strategy to estimate flows in other ungauged regions of the catchment?” Thus, this study endeavors to share a framework that employs a blend of spatial proximity (direct transfer of calibrated model parameter) and physical similarity (discharge ratio-scaling of various simulated scenarios) approaches to achieve better model predictability at the ungauged site.
2. Materials and Methods
2.1. The Study Area
Isahaya Reservoir (with an area of 26 km2 and 29 × 106 m3 holding capacity) was constructed at the innermost part of Isahaya Bay (Ariake Sea, Kyushu Island in southern Japan—Figure 1) to prevent flooding disasters and to protect reclaimed farmlands [38]. In April 1997, the inner part of the bay was separated from the main sea by a dike 7 km long and 200 m wide. Construction of the dike enclosed 35 km2 of Isahaya Bay’s tidal flats, which is 2% of the total area of the Ariake Sea [39].
Figure 1. The location of Isahaya Reservoir Catchment.
Table 1 summarizes the entire Isahaya Reservoir watershed (area 248.67 km2) with details of the ten river basins (211.34 km2). The entire catchment is ungauged apart from a single gauge upstream of the Honmyo River (R6G) basin, covering an area of 37.336 km2 with an average slope of 0.258 m/m. The 10 river basins are characterized by human interference with the course of the river for agricultural purposes, especially in the downstream sections. The climatic condition is temperate, with semi-humid, rainy summer seasons. Wet and high temperatures are experienced during the summer months, while dry and low temperatures prevail in the winter months. The average annual temperature in Isahaya is 16.3˚C. Over the course of the year, the temperature typically varies from 3˚C to 31˚C and is rarely below −1˚C or above 34˚C. On average, June is the wettest month with precipitation of about 350.0 mm and December is the driest month with total precipitation of about 65.0 mm. The region experiences annual total precipitation of 1800 - 2200 mm. Stream flow usually peaks in the middle of summer seasons because of heavy rainfall or typhoons. The information of interest for the behavior of the basin is therefore largely focused on the hydrographs of the summer season. Snow accumulation and ablation are rare occurrences due to low elevation in the Isahaya catchment. Thus, snow has a negligible effect on runoff variations in this watershed. It accounts for an insignificant portion of winter precipitation [40] [41].
2.2. Data Acquisition
Table 2 shows the available SWAT (Soil and Water Assessment Tool) model input datasets used in this research. 18 Intermittent flow rate data for each of the
Table 1. Details of Isahaya Reservoir Catchment.
R6G; Gauged River basin (upstream section of Honmyo River), AVG; Average slope, AGRR; Agricultural Land-Row Crops (mainly Vegetables-warm season), RICE; Rice farms, RNGB; Range-Brush, FRST; Forest-Mixed, FRSE; Forest-Evergreen, FRSD; Forest-Deciduous, URBN; Urban areas with high density.
NOTE: Crop management within Isahaya region was informed from Kyushu region and Nagasaki Prefectural reports and related research work. Such information encompassed, Land preparation period, types of crops, fertilization regimes, crop/land management operations, Harvesting and post-harvest operations (Fumikazu & Erdinc, 2013).
9 rivers were measured for 4 months. At the discharge point of each river, several depths and points across the river were recorded in meters to obtain the cross-sectional area and flow (velocity) captured in m/s. The measured cross-section area and the flow velocity were used to compute the daily flow rate at the discharge points in m3/s. Multiple measurements were made at each discharge point within the observational day and averaged to determine the daily mean flow rate (expounded in section 2.5).
2.3. Hydrological Characteristics of R6G Basin (Urayama Station Upstream)
The gauged part of the R6G basin was simulated using the above data set. The model SWAT was calibrated using gauge height data from the Urayama gauging station for the years (08/2016-07/2018). A sensitivity analysis was performed to determine the influence of each parameter included in the calibration on the flow rate. Model performance was then validated for 08/2018-07/2020 and 08/2019-07/2021. SWAT-CUP (a Calibration Uncertainty Program for SWAT based on the Sequential Uncertainty Conformity Algorithm-SUFI2) was used for calibration and validation and to determine the sensitivity of the calibration parameters. SUFI-2 is a semi-automatic optimization procedure that uses the Latin hypercube sampling method and is a highly efficient sampling method for obtaining optimal results. It can handle a large number of parameters, perform calibration on multiple measurement stations, and allows the use of different objective functions [42].
2.4. Extrapolation of Hydrological Characteristics from the Gauged R6G River Basin to the Ungauged River
Weather data from the three stations (Nagasaki, Omura, and Mt. Unzen) in the Nagasaki region were compared to determine seasonal variations in the region. Due to the high correlation results, the nearest weather information (Omura Meteorological Station) was adopted for the 10 streams. The 10 streams run downstream from the mountainous region east, south, and west of the reservoir with similar elevations at the discharge points (part of the Ariake sea). Based on the proportional similarity in soil characteristics, land uses, and geological history, the underground water flow characteristics, unless determined, were assumed to be like the gauged river throughout the study area.
The observed daily discharge data at the Urayama gauging-station (R6G) were denoted as QG (m3/s) and the resulting simulated flow rate of the gauged river was denoted as QSim (m3/s). Since there was no continuous data at the ungauged site, the parameters of the gauged site were integrated into the model building to obtain Qsim (m3/s). Due to the inexistence of a second gauge in the watershed, the use of methods such as global averaging or regression [43] to estimate discharge at the ungauged site was untenable. Considering the proximity and similarity between donor and recipient, the R6G parameters were globalized. Two methods were used to globalize (DG) hydrologic parameters within the watershed. First (DG), assuming similar groundwater characteristics and comparable soil hydrologic groups, the parameters were adopted from R6G (donor). For the second criterion (DG-1 & DG-2), an attempt was made to modify the groundwater lag value by considering the area of the watershed and/or the slope factor.
To determine the criterion for transplanting hydrologic flow signatures from the gauged site to the ungauged basin, different scenarios were tested separately with respect to land use/land cover and soil types to determine their influence on the flow rate. In the 9 streams, two or three soil types covering a reasonable area did not influence the changes in discharge rate. In most cases, altering the soil’s characteristics didn’t result in a significant change in flow signals. Changes in land use had a significant influence on the continuous discharge of the 10 river basins (including the gauged stream). Sensitivity analysis is critical because the new framework is based on the proportional influence of such factors on the transplantation of flow signals from the donor basin to the recipient basin.
The parameters from R6G were adopted by DG for model build-up at the 9-ungauged sites as baseline information for the watershed. Considering three land use scenarios and holding other factors constant, a continuous runoff was simulated independently for each of the following scenarios (the superscript indicates the gauged site (donor) and the lowercase indicates the (ungauged) river basin):
a) QF and Qf (m3/s)—The river basin was simulated with respective forests (deciduous, evergreen, and mixed) and high-density urban land use, taking into account the new hydrologic features. The agricultural land, i.e., rice, row crops (mainly warm-season vegetables), and Range bush (including citrus orchards), was classified as wasteland (barren land).
b) QH and Qh (m3/s)—The river basin was simulated with respect to human activities (agricultural land under crops in rows, rice, and citrus plantations) and Urban high-density land use. The land under the forests was set as barren.
c) QB and Qb (m3/s)—The river basin was simulated with barren land uses in a. and b., while the corresponding areas were maintained under high-density urban land uses.
Example R8: AGRR-24.4%, RICE 38.99% RNGB 2.92%, FRST 4.46%, FRSE 19.09% and URBN 9.00% of 24.95 km2.
a) Qf (m3/s): The FRST 4.46% and FRSE 19.09% of 24.95 km2 was maintained while AGRR-24.4%, Rice 38.99% and RNGB 2.92%, is set to barren. The URBN 9.00% of 24.95 km2, is also maintained.
b) Qh (m3/s): TheAGRR-24.4%, RICE 38.99% RNGB 2.92% of 24.95 km2, was maintained while the FRST 4.46% and FRSE 19.09%, is set to barren. The URBN 9.00% of 24.95 km2, is also maintained.
c) Qb (m3/s): The AGRR-24.4%, RICE 38.99% RNGB 2.92%, FRST 4.46%, FRSE 19.09% was set as barren. The URBN 9.00% of 24.95 km2, is also maintained.
The above scenarios were simulated independently for each river (R.No.), and the resulting continuous discharge was stored in a spreadsheet for the period 08/2016-07/2020 against the results of the calibrated river R6G. To map the dynamics of the discharge recorded by the gauging station to the ungauged basin and to use the influence of land use on the discharge characteristics (gauged and ungauged basins), the following correction factor was used. At time t, the relationship between the observed mean daily flow rate QG and the calculated mean daily discharge QT for the gauged basin was defined as:
(1)
Following that, the daily runoff at the discharge point is substantially influenced by the forest and by human activities. Taking account of the urban region as a fraction of the discharge obtained in the third scenario (QB or Qb), the total discharge at time T (Q(T)) for R6G and t (Q(t)) for the ungauged rivers was defined as:
(2a)
(2b)
where:
is a fraction of urban area with high density land use and/or fraction within each scenario if changed to barren.
Considering the differences in flow rate at each time t-step, a correction factor ratio
was determined as follows depending on the influence of the land uses within the donor and recipient river basins.
(3)
where:
is the correction factor ratio at time t and
and
are the fraction of the urban high density land use of the gauged and ungauged river basins, respectively.
The estimated daily flow at the ungauged discharge site was determined as follows to account for hydrologic differences in flow signatures within the donor and recipient when transferring flow characteristics from the observed data.
(4)
A comparison was made between the discharge predicted by the new framework method (N.F.), i.e., Qung (m3·s−1) and the 18 observed intermittent discharges (obtained between March 2021 and July 2021) to determine the performance of the applied method used at a given discharge point for each of the 9 river basins.
2.5. Stream Flow Measurement with the KENEK VPT-200-24
Discharge estimates were obtained by measuring river flow with the Kenek VPT-200-24 flow meter and cross-sectional area (using a depth staff and tape measure). The Kenek VPT-200-24 flow meter is an ideal instrument for measuring flow in rivers and lakes or reservoirs with a wide range of water quality (tap water quality to industrial wastewater).
The flow rate (Qm in m3/s) was calculated from the velocity profile using the mean-section method, which uses the depth dn (m) at a series of two verticals in a cross-section and the distance between the verticals bn (m) to determine the cross-sectional area for each subsection. The depth-averaged flow velocity
(m2·s−1) in each vertical was estimated by measuring the velocities at several points along the nth vertical. Cross-sectional discharges are computed as the product of the segment area and the depth-averaged velocity. Total discharge was computed as a summation of sectional discharges as defined in Equation (5) and Figure 2.
Figure 2. Sample of stream flow measurements (discharge point-R8) using VPT2-200-24.
(5)
where for example;
R8, 13 subsections were considered as represented in Figure 2.
Assuming the edges toward the river banks were approximately triangular, the edge coefficient of 0.3535 was adopted [44]. This applied to all other 8 river channels except for river R9 channel. Thus, at the edges, discharge was defined as:
(6)
and
(7)
The total discharge QT was then defined as:
(8)
At the discharge point of the river, R9 is a well-defined rectangular weir. The discharge was therefore determined to be equal to Qm with the numbering of verticals running from the end to the end of the riverbank (1, 2, 3 … nth).
2.6. Definition of the Objective Functions
To calibrate the model and determine how optimally the simulation reproduces the observed information, the objective functions: NSE, R2, and bR2 objective functions were considered.
1) Nash-Sutcliffe model efficiency coefficient (NSE)
The NSE was used to evaluate the predictive power of the model that replicates the hydrologic functions in the studied watersheds. It was defined as:
(9)
where:
is the mean of observed discharge,
is modeled discharge, and
is observed discharge at time t.
NSE range is −∞ ≥ 1; Values of NSE ≥ 0.5 were satisfactory.
2) The Pearson Coefficient of Correlation, R2
R2 is a measure, of how well-observed results are replicated by the model simulation based on the proportion of total variation in the results. The Coefficient of Correlation R2 was defined as:
(10)
where: Q is a variable (discharge), and m and s stand for measured and simulated, i is the ith measured or simulated data.
R2 range is 0 ≥ 1; Values of R2 ≥ 0.5 were satisfactory.
3) Modified coefficient of determination, bR2
The coefficient of determination R2 is multiplied by the slope of the regression line b between the simulation results and the measured data. This function takes into account both the discrepancy in the magnitude of two signals (depicted by b) and their dynamics (depicted by R2) [45].
(11)
bR2 range is 0 ≥ 1;Values of bR2 ≥ 0.5 were satisfactory.
Regarding to the merits and demerits associated with each one of the three objective functions [32] [45] [46], the three were used to complement each other’s shortcomings and evaluate the precision of the model prediction.
3. Results and Discussion
3.1. Calibration and Validation of the Gauged R6-G River Basin
SWAT-CUP SUF2 was used to calibrate and validate of the model output against the available gauged data. The performance of the model was recorded in Table 3 for the calibration periods (1 and 2). The resulting hydrograph (Figure 3) showed that the peak signatures, decline, and base flows were accurate compared to the observed data at the Urayama gauging station.
As recorded in Table 3, the model performance improved remarkably during validation period 1 & 2, depicting the flow rate and the trend at the gauge station (R6G). Although the number of daily flow rates covered within the 95PPU band decreased, the shrink in the band by half and the increase in NSE and R2 indicated the precision. Despite the difference in peak flow in the second validation, the model performed adequately. The high precision of these results could imply
Table 3. Model performance during calibration and validation for the gauged R6G basin.
Figure 3. The flow characteristics of R6G basin during calibration and 2nd validation period.
that: The model choice was appropriate for this mountainous river basin and that the input data, the crop schedule, and the weather data were true representatives of the region.
As discussed by various articles, this underscores the importance of model choice and correction of data from measuring instruments or other related accrued errors before use in modeling [47] [48]. Successful use of the distributed model to simulate river discharge depends on the characteristics of the input data, observational data for calibration, and the selection of an appropriate model for the study region. This ensures the robustness of the model results and adequate simulation of the flow responses of a watershed [46]. The choice of a model and the quality of the data used are reflected in a successful simulation, i.e., parameter values within a reasonable range and without overparameterization. This ensures proper calibration and validation of model results against a specific measured hydrologic feature [49]. The R6G river basin was adequately calibrated for 08/2016-07/2018 observed data and validated against a new set of data 08/2018-07/2020 and 08/2019-07/2020, as captured in Table 3, and partly in Figure 3, above.
During calibration, 9 parameters were used as recorded in Table 4. The river routing parameters were the most sensitive in the R6G basin. The sensitivity of the parameters during calibration follows the list in Table 4, column 1, denoted by a rank number. Groundwater delay (days) had the least influence, while the base flow factor for bank storage had the highest influence on the flow rate response during calibration of the SWAT model. The best simulation was realized when the model parameters were fitted at the values in Table 4, column 4. The mean daily flow for the basin post-calibration model was 2.48 m3/s with the maximum
Table 4. The parameters utilized in the calibration of the SWAT model for the R6G basin.
peak during the summer period (June-July) and lowest flow in the winter period (December-January) compared to 2.45 m3/s as observed.
3.2. Transplant of Hydrological Flow Characteristics from Gauged to Ungauged Basins
The parameters altered in calibration and validation of R6G provided basic hydrological characteristics of the Isahaya catchment (since exists only a single gauge in the whole catchment). The assumption of approximate homogeneity in the cases of DG makes the method very simplistic. Adopting all the parameters’ values (DG) as attained in R6G, the resultant hydrograph correlated reasonably with the observed intermitted dataset but overestimated or underestimate the mean daily flow rate. Modifying some of the globalized parameter(s) using an implicit equation could improve the accuracy of estimated discharge at an ungauged site. Using such implicit equations (Equation (12) or Equation (13)) as DG-1 and DG-2 below to alter the groundwater delay only improved precision in depicting the trends for R1, R2, R5, and R7 as compared to the observed dataset.
(12)
or
(13)
where; GW_D is ground water delay, R.No is river number (1, 2, 3…), A is the watershed area, SLP is the avarage slope.
Figure 4 and Table 5 are the results of DG, DG-1 & DG-2 attempts to characterize continuous flow at the ungauged site. The overestimation or underestimation of the flow rate and the absence of crests during large rainfall events
Figure 4. The flow characteristics at ungauged R1, R5, R7, and R8 using DG against observed data.
Table 5. Evaluation of flow characteristics at the ungauged site using direct globalization.
were some of the misrepresentations of the hydrography at the ungauged site using the methods DG, DG-1 & DG-2. Although the above methods represent the trend at ungauged sites, they are abstract because river basins have some heterogeneity despite their proximity to the donor. Changes in calibration parameters are relative to other parameters at a given site despite similarity and proximity. This was evident in the way different river basins responded to changes in groundwater delay value as per Equation (12) and Equation (13). In rare cases, a changing a single parameter can lead to an optimal simulation. In the case of a single gauging station, it was not possible to change parameters relative to other parameters due to scanty hydrological information about the watershed. Due to these uncertainties, the method described above is not reliable and cannot be applied to other sites with similar hydrological characteristics.
Using DG as baseline information that can be assigned to the watershed, flow characteristics were estimated at an ungauged station for each of the 9 river basins according to the new framework (N.F) described in Section 2.4. Precision and reproducibility were evaluated, and the baseline information from DG was incorporated into the modeling prior to applying the N.F. Although DG-1 and DG-2 showed better trend precision for R1, R2, R5, and R7, the resulting N.F hydrograph was close to the DG when used in model building for the N.F. To determine whether the N.F. estimates reflect the actual flow rate at the ungauged discharge point, the characteristics of the resulting hydrographs were statistically compared with the 18 intermittently observed flow records.
The continuous flow rate mapped from the gauged stream to the ungauged site by the N.F showed representative estimation as compared to the use of the DG method (Figures 5(a)-(d)). The N.F exhibited robustness in estimating base flow and peaks as a response to reasonable precipitation events. The problem of overestimating base flow or missing the peaks even with sufficient precipitation events was correctly represented by the N.F. method of transplanting signatures relative to the donor. This phenomenon was much pronounced in R1 & R8, R3 & R10, R5, and R9 as represented in Figures 5(a)-(d), respectively. DG is associated with certain levels of uncertainty, depending on the proximity and similarity between the recipient and the donor watershed [50]. Maintaining such rigidity of flows would adversely affect the physical and biochemical analysis of additional watersheds. Non-responsiveness of flow signals could be due to the assumption that the physical catchment attributes of the gauged basin adequately mirror those of the recipient, thus proportionately influencing changes in flow rate.
The N.F., as applied to the 9-ungauged sites, proved suitable for characterizing continuous flow rate in a data-scarce watershed (summarized in Table 6). Although the accuracy of the new framework was lowest for R1 and R2, the result was a reasonable estimate compared to the observed flux data (Figure 5(a) and Figure 5(b)). In the donor, R-6G, the watershed area is 71% forest and 21% farmlands, while the R1 basin is 90% forest and only 7% of the land is used for farming activities. The R1 basin, too, has the highest slope in the whole catchment. These factors and proximity to the donor basin may have influenced the results of DG and N.F. compared to estimates at discharge points on other rivers. However, the net effect was more pronounced with the DG method than with the N.F. method. The N.F. method’s prediction best depicted the trend and response of flow to precipitation events. The results are reliable for further analysis of the watershed. The R2 difference in land use coupled with a lower slope as compared to the donor might have influenced the transplant of flow signals.
The R3 river basin is relatively small (about the same size as R1 and R2), but the similarity in slope and proportion of land use to the donor catchment can be related to the high accuracy of the N.F. (Table 6 and Figure 5(b)). The R8 basin has the lowest gradient and is the least forested compared to the other 9 river basins draining into the Isahaya Reservoir. The land use characteristics are quite different from the gauged river basin (R6G), where more than 65% of the land is used for agriculture (for the other 7 rivers, more than 70% of the land is covered by forests). These differences didn’t influence the transplanting of flow signals from the donor. The sharp decline in the observed intermittent flow rate towards
the end of July was of concern and may indicate that much water was diverted to rice farms due to reduced rainfall. From the results of the 9 river basins shown in
Table 6. Summary of the N.F precision as applied at the ungauged basins to estimate flow.
W. Area: Watershed Area, D-flow: Daily Mean Flowrate, %FRST & %HUM.A: Percent Watershed area under Forests and Human activities (mainly farmlands).
Table 6, it is evident that the N.F provides reasonable and reliable results in transferring the hydrologic continuous discharge characteristics of the gauged donor basin to the ungauged recipient basin in the ungauged watershed.
Proportional scaling utilizing land use influence on discharge between the two basins showed adequate predictability and reduction in the loss of flow signatures in the event of direct globalization. Thus, the N.F. could be used in transferring runoff characteristics from gauged to ungauged river basins in a watershed with scarce flow data or without other hydrologic information. Despite the difference in size, slope, or land use, the method showed adequate flow response to weather and catchment characteristics at the ungauged site. The N.F. provides an alternative method for reasonably estimating discharge at an ungauged location when the number of gauges is limited, and other methods are inapplicable or ineffective.
The attributes used in the regionalization of flow characteristics in the pursuit of estimating the flow rate at the ungauged site include proximity, the catchment area, elevation, the slope of the basin or channels, land use, and weather conditions [29] [51]. In the Isahaya Reservoir watershed, area, slope, and land use are important aspects that affect the flow regime and influence the choice of method that can be used to transplant flow characteristics. The problem of data unavailability at the onset of this study raised the question; of what approaches could be adequate to estimate the flow rate at the ungauged basins with a single gauging station [52].
Using only a single gauge and several observation data, the N.F could be used to estimate flow characteristics of the ungauged basin. Where the direct parameter transfer method is impractical (e.g., lack of averaged watershed parameter values or regression formula relating to hydrological characteristics within or outside the watershed), the N.F would aid in refining the estimates by accounting for differences or similarities with the donor basin [29]. When a rainfall-runoff model is to be selected, the model efficiency and applicability are some of the vital factors that a researcher needs to consider in the optimization of flow estimates in such a catchment with limited data [16] [53] [54] [55]. Selecting parameters for calibration or optimization needs to be considered, as well as the inclusion of associated uncertainties [56] [57] [58]. Due to proximity and similarity within the river basins, DG yielded a fair trend of flow rate estimates at the ungauged site but was inconsistent in determining the actual flow rate with changes in weather. The N.F. was able to refine the results of DG and provide an adequate estimate of discharge at the ungauged site.
3.3. The Uncertainties
3.3.1. The Uncertainties Inherent in the Hydrological Modeling
In SUFI-2, the parameter uncertainty is determined by all the input and output sources of uncertainties. These include the uncertainty in the input rainfall data, the land use and soil type, parameters, and observed data. During simulation, uncertainty is quantified by the 95% prediction uncertainty (95PPU) which is referred to as the p-factor. The 95PPU was calculated at the 2.5% and 97.5% levels of the cumulative distribution function of the output variable obtained by Latin hypercube sampling. The best calibration and parameter uncertainty is measured based on its closeness to the p-factor of 100% (the observations bracketed by the prediction uncertainty) and the r-factor of 1 (size of the uncertainty band). When the two factors are within a satisfactory range of values, a uniform distribution in the parameter hypercube is explained by the following parameter distribution. The goodness fit in SUFI-2 was quantified by the coefficient of Nash-Sutcliff model Efficiency (NSE) and Pearson correlation coefficient (R2) between the observation data and the best simulation [36] [56].
Figure 6 shows a 95% prediction interval (dotted band) that indicates a 95% confidence interval for predicting a given day’s mean flow rate of a specific river using the N.F model. The prediction interval is wider than the confidence intervals (grey band) because predicting the flow of the ungauged river on a given day by transplanting the hydrologic characteristics of a gauged stream has greater uncertainty than predicted based on the mean flow, which is proportional to the mapped/transplanted annual flow characteristics.
3.3.2. The Uncertainty in Discharge Measurement
The measurement process and use of data reduction equations provide the basis needed for conducting the uncertainty analysis. The analysis is vital in deducing the confidence in discharge determined for other institutions/researchers that would base their analysis on the already conducted research. The following are the computed uncertainties involved in the flow measurement at each discharge point (Table 7).
Figure 6. Regression plots showing prediction uncertainty of the flow rate (Q m3/s) transplanted by the N.F in comparison to the observed intermitted datasets.
Table 7.1 Computed uncertainty for the discharge measurement at each site.
The low uncertainties at R9 (Table 7) are due to a triangular weir at the measurement site, implying a high probability of making similar measurements by different users. In all rivers, wading staff and ±2%/FS contributed significantly to the uncertainties. The uncertainties mentioned may have arisen during actual measurements at the discharge point. Other uncertainties, such as changes in flow outside the observation period, were not quantified. Such uncertainties should be considered when using the data or information from this work for further studies or analyses.
4. Conclusions
Natural systems in a watershed, such as rainfall-runoff processes, are complex, dynamic, nonlinear, and governed by multiple interrelated physicochemical components. Therefore, the successful use of simulations to describe the flow regime in a river basin requires the proper choice of a model, high-quality inputs, and observational data. The uncertainties associated with modeling such complex processes affect the reliability of the prediction of runoff at gauged and ungauged sites. Therefore, measured data are required for such models to force multiple variables to represent the process in the watershed and to ensure consistent and reasonable prediction of the stream flow. The unavailability of such data in many watersheds has been and continues to be an obstacle in evaluating river resources and or receiving reservoirs, dams, estuaries, and lakes.
In the case of limited gauges within a watershed, e.g., a single gauge in the Isahaya Reservoir watershed without available hydrologic information, direct regionalization of such discharge information to the ungauged site could be erratic and unreliable. This is due to differences in gradient, land use, soils, weather conditions, catchment area, and river basin shape. Therefore, the new framework presented in this study contributes to the hydrologic characterization of watersheds that have a limited number of gauges. Based on the differences in land use, slope, and size of the Isahaya Reservoir catchment, the new framework showed consistency and reliability in transplanting flow dynamics to the ungauged site.
In the Isahaya Reservoir watershed, differences in the river basins are not considered extreme. Therefore, further studies need to be conducted to determine the applicability of the new framework in cases where the donor and recipient river basins are vastly different. In addition, the framework should be further developed for watersheds with different weather conditions and subsurface hydrologic conditions. Depending on the assumptions made for this study, the applicability of the new framework was successful in the Isahaya catchment. This provides a flow regime for 10 river basins that would form the basis for further basin and/or water quality analyses. While this study does not provide an exclusive method for transplanting flow signatures from gauged river basin to ungauged river basin, it does provide a reliable method for estimating and characterizing continuous flow rate in a watershed with limited gauging stations and relatively similar weather conditions.
Acknowledgements
Support for this work was provided by the Nagasaki University, Engineering department. Authors are thankful to Riverine laboratory members; Mr. Nakamura, Mr. Mukai, and others for helping with field observations and emotional support. Authors also appreciate Dr. Kilonzi for proofreading our work.
Appendix
Appendix 1. Derivation of Equation (2)
Different land use influence characteristics of the daily flow rate at discharge point differently. The Characteristic of the hydrograph is the summation of partitioned individual land uses within the watershed. Within river basins in Isahaya catchment, Urban High-Density Land Use had no substantive influence on the characteristics of the mean daily flow rate. Agricultural lands (rice fields, and range bush (citrus fruits bush,)) which are major human activities in the catchment, and the forests (deciduous, Evergreen, and mixed) had a significant influence on the daily mean discharge for all the river basins. Arbitrary letting,
1) QF—Portion of discharge influenced by/from forest cover;
2) QH—Portion of discharge influenced by/from human activities related to the land use;
3) QU—Portion of discharge influenced by/from a fraction of watershed that is under Urban High Density;
4) QB—Portion of discharge influenced by/from a fraction of watershed that is set to barren during scenario simulations;
5)
—A fraction of watershed under Urban High Density/barren land use (only defined in scenario contest; was not among land use in the catchment).
Let the total discharge Q at the exit of the river be influenced by three land-uses as
(a)
In a scenario where a section of Forest or Human activities land use is changed to barren, then change in the total discharge was arbitrarily defined as
(b)
Given that the influence of QU is negligible within the study catchment, the urban area was assumed to have a similar influence on the hydrographs as barren land cover with respect to the definition and use of the formula. The Equation (b) was then rewritten as
(c)
where,
is a fraction of the watershed under urban high density and/or barren if so defined. The difference in discharge at time T (days) was defined by Equation (c) that could be re-written as
(d)
where
. At minimum before cosidering any of the three scernarios
Based on the scenarios in this study, the daily mean discharge is approximately equal to the QT as defined in Equation (d). The above equation Eqn.d. formed an integral part of the correction factor ratio Cf.
Appendix 2. Calculation of the Uncertainty in Flow Measurement
The following are some of the sources of the uncertainties involved in the estimation of the final discharge estimate.
1) Instrumental resolution
This entails a value corresponding to half of the last significant digit of the software output for velocity (0.001 m/s). The user has no capability to making any record beyond the last significant digit and primarily depends on digits outputted from the software. The same applies to the measuring tape and wading road; the last significant digit that can be read by the user.
2) Instrument accuracy
This was based on the manufacturers’ specifications for the instrument accuracy of ±2%/FS (±0.01 m/s). The accuracy is specified as % FS hence the error is applied as a fixed value no matter the output. During calibration of the machine and test in the laboratory, the precision was determined as ±1 cm/s. Using fine laboratory measurement and the use of the wading rod in the determination of depth, a discrepancy of ±0.001 and ±0.015 was noticed respectively.
3) Exposure/averaging time
These are fluctuations that appear as random noise in measurements taken by an instrument with high-frequency sampling capabilities such as the Kenek VPT-200-24. The ‘‘noise’’ is majorly due to flow turbulence. Other sources of such noise that could accumulate in the output signal due to instrument noise (electronic noise, Doppler noise), and environmental noise (wave effects, erratic operations by person). To make one observation at a given locale, the machine was set to 40 seconds for numerous samplings and then recorded output as averages. Records taken at a single point (continuous records, 40 seconds and reset for new reading) indicated average differences of 0.25% Table 6. (in the main text) summarizes the uncertainty sources associated with the stream discharge measurement using Kenek VPT-200-24 [59] [60].
4) Flow angle induced error
This is a flow-angle error due to incorrect inclination to flow and during calibration of the compass being used. The error was estimated as the manufacturer’s specification of ±3˚/360˚ for the calibrated flow meter. This is an important error to be considered since it causes overestimation in sub-section discharge, respectively as discussed by Huang, 2012 [61]. The inclination of the machine was adjusted in each point velocity before the determination of the discharge. However, the flow angle error arising from calibration was not considered as the instrument was well calibrated before use. It was assumed that the error was less than 1% thus negligible.
5) Operational conditions
This error is because of the non-uniformity of suspended scatters in the electromagnet beams, change in flow direction during measurement, machine deflections, disturbance of the water surface by waves, and operator-induced effects. This is influenced by procedures used in field measurements and/or natural conditions during measurement [62]. Ten repeated measurements at the same vertical were made in the field experiment to try and quantify this type A error for the flow meter. Using a wadding rod for measuring the depth, angle inclination, positioning of riverbed and miss-readings are some of the errors that could be incurred in field observation. Considering measurement at the same location by several users, the error was defined. The standard uncertainty associated with the main measured variables (mean velocity in the vertical, depth, and width) were then aggregated through the following relationships (scripts defined in Table 5):
(e)
(f)
(g)
The standard uncertainty for the measured discharges of R-8 (example) obtained with Equation (e) was expressed as combined standard uncertainty as:
(h)
The contribution of the correlated uncertainties in two width measurements in the sub-area was assumed to be zero for this research.
For a sub area discharge at the nth measurement, the functional relationship was defined as [62]:
(i)
Sensitivity coefficients, c, for each input were found by taking the partial derivatives for each independent variable in the above equation:
(j)
(k)
(l)
Combined standard uncertainty was then defined as
(m)
The expanded uncertainty was determined as follows.
(n)
The coverage factor of
was adopted corresponding to normal distribution giving the estimate a 95% confidence interval. Table 7 shows the uncertainty that was possibly propagated during flow measurement within the study period. The seasons determined how many sub-sections could be undertaken throughout the research time.
NOTES
1Appendix 2. Calculation of the uncertainty in flow measurement.