Shape Of Hydrograph Research Articles

Calibrated Eckhardt’s filter versus alternative baseflow separation methods: A silica-based approach in a Brazilian catchment

The Eckhardt’s filter, a recursive digital filter widely used in hydrology to separate baseflow, relies on two parameters: the recession constant (a) and the maximum allowable long-term ratio of baseflow to streamflow (BFImax). While a can be determined through regression analyses, BFImax is often assigned an arbitrary value because of limited catchment information for a more objective determination. This study aimed to assess and calibrate BFImax for a 1.2 km2 rural catchment in southern Brazil using dissolved silica data. Calibration involved minimizing differences between filter-generated baseflow hydrographs and those derived from dissolved silica measurements and mass balance across 13 monitored rainfall-runoff events. The calibrated model was compared to three established models: a modified graphical method based on local minima, a backward filter model, and a one-parameter filter. The results yielded an optimal BFImax of 0.653 with high model performance (PBias ≈ 0, NSE = 0.85, KGE = 0.74, and NRMSD = 13 %) against silica-derived baseflow. Applying the calibrated filter to a six-year streamflow dataset yielded a long-term baseflow index (BFI) value of 54 %, aligned with regional and global averages. Comparisons with other models showed that the backward filter model produced unrealistic baseflow values nearly identical to streamflow values; the local minima method overestimated baseflow by +19 % and failed to accurately represent the smoothness and peaks of the baseflow hydrographs; and the one-parameter filter represented the shape and timing of the baseflow hydrographs well, but underestimated magnitude by −19 %. These findings validate the effectiveness of the Eckhardt’s filter in generating accurate baseflow hydrographs and confirm the importance of calibrating BFImax for optimal performance. The research also endorses the use of dissolved silica measurements and the mass balance method as a reliable approach for calibrating BFImax for long-term assessment of baseflows.

Innovative control strategy for flood mitigation through a combination of barrage management on impounded rivers and polders

The partially catastrophic effects of large floods, such as the events of 2002 and 2013 in Germany and Austria for example, clearly show the vulnerability of settlement structures along large rivers. In the case of impounded rivers with run-of-river hydropower plants, the main question is whether a reduction of floods can be achieved with flood-adapted, intelligent control of barrages - to date, the operational objective has been the safe discharge of floodwater; targeted reduction is currently an unoccupied field due to its complexity. Specifically, the interaction with the operation of a polder for flood retention is of special interest, as polders are a much-discussed and effective element of flood protection, especially at higher discharges, that directly influences the hydrograph shape. Extensive analyses and sensitivity investigations of a multitude of parameters and measured values, using a complex model system, are needed in order to investigate the reduction potential of the named technical flood protection measures and the development of control strategies relevant to practice. It is essential to base the developed controlling strategies on operationally available data, as this is the only way to make the determined reductions relevant for practical operation.

Technical note: Testing the connection between hillslope-scale runoff fluctuations and streamflow hydrographs at the outlet of large river basins

Abstract. A series of numerical experiments were conducted to test the connection between streamflow hydrographs at the outlet of large watersheds and the time series of hillslope-scale runoff yield. We used a distributed hydrological routing model that discretizes a large watershed (∼ 17 000 km2) into small hillslope units (∼ 0.1 km2) and applied distinct surface runoff time series to each unit that deliver the same volume of water into the river network. The numerical simulations show that distinct runoff delivery time series at the hillslope scale result in indistinguishable streamflow hydrographs at large scales. This limitation is imposed by space-time averaging of input flows into the river network that are draining the landscape. The results of the simulations presented in this paper show that, under very general conditions of streamflow routing (i.e., nonlinear variable velocities in space and time), the streamflow hydrographs at the outlet of basins with Horton–Strahler (H–S) order 5 or above (larger than 100 km2 in our setup) contain very little information about the temporal variability of runoff production at the hillslope scale and therefore the processes from which they originate. In addition, our results indicate that the rate of convergence to a common hydrograph shape at larger scales (above H–S order 5) is directly proportional to how different the input signals are to each other at the hillslope scale. We conclude that the ability of a hydrological model to replicate outlet hydrographs does not imply that a correct and meaningful description of small-scale rainfall–runoff processes has been provided. Furthermore, our results provide context for other studies that demonstrate how the physics of runoff generation cannot be inferred from output signals in commonly used hydrological models.

Hydrological dynamics of snowmelt induced streamflow in a high mountain catchment of the Pyrenees under contrasting snow accumulation and duration years

AbstractSnowmelt drives a large portion of streamflow in many mountain areas of the world. However, the water paths from snowmelt to the arrival of the water in the streams are still largely unknown. This work analyzes for first time the influence of snowmelt on spring streamflow with different snow accumulation and duration, in an alpine catchment of the central Spanish Pyrenees. This study presents the water balance of the main melting months (May and June). Piezometric values, water temperature, electrical conductivity and isotope data (δ18O) allow a better understanding of the hydrological functioning of the basin during these months. Results of the water balance calculations showed that snow represented on average 73% of the water available for streamflow in May and June while precipitation during these months accounted for only 27%. However, rainfall during the melting period was important to determine the shape of the spring hydrographs. On average, 78% of the sum of both the snow water equivalent (SWE) accumulated at the beginning of May and the precipitation in May and June converted into runoff during the May–June melting period. The average evaporation‐sublimation during the 2 months corresponded to 8.4% of the accumulated SWE and rainfall, so that only a small part of the water input was ultimately available for soil and groundwater storage. When snow cover disappeared from the catchment, soil water storage and streamflow showed a sharp decline. Consequently, streamflow electrical conductivity, temperature and δ18O showed a marked tipping point towards higher values. The fast hydrological response of the catchment to snow and meteorological fluctuations, as well as the marked diel fluctuations of streamflow δ18O during the melting period, strongly suggests short meltwater transit times. As a consequence of this hydrological behaviour, independently of the amount of snow accumulated and of melting date, summer streamflow remained always low, with only small runoff peaks driven by rainfall events.

The effects of vegetative feedbacks on flood shape, sediment transport, and geomorphic change in a dryland river: Moenkopi Wash, AZ

Since the 1950s, Moenkopi Wash, in Arizona, United States, has been transformed from a relatively wide river at low flow, to a narrow, heavily vegetated river that is less than half of its former width. We analyzed a ~95-years-long instantaneous-discharge record, an extensive sediment-transport record, oblique and aerial photographs, historical channel surveys, and historical stage-discharge rating relations to determine the primary mechanisms responsible for this transformation.Frequent, large floods dominated the early part of the discharge record between 1926 and 1940. A dramatic ~30 % decrease in annual mean discharge, and ~27 % decrease in the mean of the partial-duration flood series occurred in 1941. This decline did not result in widespread channel change but provided an opportunity for vegetation to establish along the channel margins. Widespread channel narrowing began after a second decline in mean and peak discharge in 1956 at which time the channel banks became heavily vegetated. Between 1952 and 2019, the river channel narrowed by 57–59 %. Approximately 2–4 m of bed aggradation occurred at most study sites. As the channel narrowed and became more vegetated, large floods maintained the narrower channel width but did not cause rewidening. Suspended-sediment transport data illustrate the sediment trapping effects of vegetation that lead to channel narrowing; suspended-sand concentrations decline, thereby indicating sediment deposition, as vegetation becomes progressively inundated during floods. Furthermore, dense channel-margin and floodplain vegetation provide increased roughness, resulting in flood-peak attenuation and alteration of hydrograph shape. Vegetation expansion causes positive feedbacks whereby sediment deposition during floods is exacerbated by the roughness and sediment-trapping effects of vegetation leading to further narrowing. These positive feedbacks have reduced sediment delivery to the Little Colorado River downstream. Channel widening is not likely to occur unless there are very large floods that exceed the erosional threshold of the channel-margin vegetation, or unless large-scale vegetation removal efforts are undertaken.

Flood Attenuation Potential of Italian Dams: Sensitivity on Geomorphic and Climatological Factors

In this work the attenuation potential of flood peaks of 265 large reservoirs all over Italy is analysed, considering a flood management that excludes gates opening and then configures strictly unsupervised attenuation effects. Key factors of dams and related basins are considered to develop a ranking method that can emphasize the interplay between dam geometry and the hydrological processes acting in the upstream watershed. To maintain a homogeneous approach in such a wide geographic area, the attenuation index is computed applying the numerical solution of the differential equation of lakes and only two different standardized hydrograph shapes have been used. An index design flood from the rational method is used as the incoming peak value for each dam, enhancing the use of the results of a recent analysis of all Italian rainfall extremes. Even with a very simple approach, twenty-four different design incoming floods are derived, by varying the shape of the incoming hydrograph and the parameters of the rational method. Exploring the ranking results in all the alternatives, the attenuation potential obtained for all dams demonstrates to be strongly sensitive to the assumptions on the time of concentration and to some rainfall features. On the other hand, the hydrograph shape seems to exert much less influence on the ranking outcome. Results obtained can be useful to studies of wide-area flood frequency analyses, as we highlighted the sensitivity of the rank of attenuation efficiency to hydrologic parameters widely used in the assessment of the design flood peaks in ungauged basins.

PREDICTION OF RUNOFF BY SYNTHETIC UNIT HYDROGRAPH METHODS FOR THE DESIGN STORMS IN WARANA RIVER BASIN, MAHARASHTRA, INDIA

The authors' goal was to calculate average representative unit hydrographs (UHs) from a few isolated storms for the upper Warana River basin in Maharashtra, India. Based on the approaches used by Snyder, the Soil Conservation Service (SCS), the Central Water Commission (CWC), and Common, synthetic unit hydrographs were created. This study was utilised to choose appropriate correlation coefficients for Synthetic Unit Hydrograph (SUH) evaluation for the ungauged basins with identical hydrological conditions. The outcomes were then analyzed with those attained using other techniques and with the actual UH. In order to calibrate and validate the effectiveness of Snyder's model, the authors employed it on the sub-basins of the Warana river basin. The hydrograph's shape resembled the observed UH produced by the other methods, according to the results, and the peak discharge calculated using Snyder's technique had the lowest percentage error (0.26%)

Reservoir Release Allocation and Weir Design MATLAB® toolbox

ABSTRACT Designing reservoir outlets is a challenging subject in flood management among water resource managers. However, although several kinds of research are dedicated to flood mitigation strategies, studies considering flood uncertainties in reservoir outlet designs are rare in the literature. This study fills this knowledge gap by developing a MATLAB® Toolbox application, entitled Release Allocation and Weir Design (RAWD), based on MATLAB® App Designer, to obtain optimum reservoir outlet designs for effective flood management considering four outlets types, including an orifice as well as a proportional, a triangular, and a rectangular weir. The RAWD application handles uncertainties using non-linear interval number programming (NINP) and conditional value-at-risk (CVaR) approaches within an optimization model. Results show that the RAWD Toolbox application’s user can effectively and accurately meet downstream water demands and minimize water allocation shortages during the flood for a reservoir subjected to different shapes of inflow hydrographs.

Constraining flow and sediment transport intermittency in the geological past

Quantitative investigations of ancient rivers usually provide insights into either instantaneous or mean flow conditions. There is a critical gap between these time scales of investigation, which reflects the intermittency of flow and sediment transport, and closing this gap is crucial to fully explore the dynamics and evolution of ancient fluvial landscapes. Here, we combined fluvial stratigraphic data sets, flow and sediment transport models, and paleoclimate general circulation model (GCM) results to develop new methods to estimate intermittency in the geological past, specifically flow intermittency factors (Iw) and sediment transport intermittency factors (Is), and we show how they can be used to explore past hydrograph shapes. We illustrated these methods using the Upper Cretaceous Last Chance Ferron Sandstone in Utah, western United States. For sand-transporting flow conditions in Last Chance Ferron rivers, we estimated Iw values of 0.54−0.90, which imply that channel-forming conditions were sustained for the majority of the year, consistent with perennial systems in which relatively large discharges were sustained. In contrast, for gravel-transporting flow conditions, Iw values of 0.28−0.38 suggest that the largest formative flows may have occupied Last Chance Ferron rivers for nearly a third of the year, which could be explained by a monsoonal system in which high-magnitude discharge events were sustained, or a subtropical system in which high-magnitude discharge events had short durations but high frequencies. Meanwhile, Is values of 0.075−0.15 suggest that annual sediment budgets could have been transported in as little as 10 days and up to 2 months, if channel-forming conditions were sustained, and these values highlight that small changes in the duration of channel-forming conditions could significantly impact sediment budgets. Our results are consistent with independent facies- and proxy-based insights into Last Chance Ferron rivers, which point to a perennially wet system characterized by a monsoonal or subtropical discharge regime. Further, our results highlight new opportunities to use paleoclimate GCMs to constrain intermittency factors in the geological past. Going forward, paleoclimate GCMs will be particularly useful where the rock record is incomplete or inaccessible, or where stratigraphic approaches are limited, and they will enable us to tackle pertinent research questions pertaining to past surface processes on both Earth and other planets.

Stochastic Flood Simulation Method Combining Flood Intensity and Morphological Indicators

The existing flood stochastic simulation methods are mostly applied to the stochastic simulation of flood intensity characteristics, with less consideration for the randomness of the flood hydrograph shape and its correlation with intensity characteristics. In view of this, this paper proposes a flood stochastic simulation method that combines intensity and morphological indicators. Using the Foziling and Xianghongdian reservoirs in the Pi River basin in China as examples, this method utilizes a three-dimensional asymmetric Archimedean M6 Copula to construct stochastic simulation models for peak flow, flood volume, and flood duration. Based on K-means clustering, a multivariate Gaussian Copula is employed to construct a dimensionless flood hydrograph stochastic simulation model. Furthermore, separate two-dimensional symmetric Copula stochastic simulation models are established to capture the correlations between flood intensity characteristics and shape variables such as peak shape coefficient, peak occurrence time, rising inflection point angle, and coefficient of variation. By evaluating the fit between the simulated flood characteristics and the dimensionless flood hydrograph, a complete flood hydrograph is synthesized, which can be applied in flood control dispatch simulations and other related fields. The feasibility and practicality of the proposed model are analyzed and demonstrated. The results indicate that the simulated floods closely resemble natural floods, making the simulation outcomes crucial for reservoir scheduling, risk assessment, and decision-making processes.

Detection of decadal time‐series changes in flow hydrology in eastern Australia: Considerations for river recovery and flood management

AbstractSince European colonisation, many coastal rivers of New South Wales (NSW) have been highly modified by channelisation, flow regulation, sediment mining, intensive grazing and agriculture, deforestation, and riparian vegetation and wood removal. Since the late 1980s, nearly 55% of these rivers have undergone a significant ‘re‐greening’ and geomorphic recovery with changes to instream and riparian roughness. However, little research has been undertaken to quantify if there have been any coeval changes in flow hydrology. Approximately 7000 flow hydrographs with one‐hour time‐steps from 117 gauges on 45 study rivers (17 of 20 coastal catchments of NSW) were used to assess changes to in‐channel and overbank hydrology over decadal timeframes. Hydrograph shape and attenuation characteristics for three morphologically defined flow stages were analysed: in‐channel fresh, high flow and overbank flood. Time‐series analysis was used to quantify changes in flow hydrology from the early‐20th century to present and identify trends occurring between known degradation phases (1910s–1940s and 1940s–1980s) and a recovery phase (1990s to present) for these rivers. Our findings indicate that, on average, and for some key examples, changes in flow hydrology are occurring. Flows are slowing and attenuating as indicated by decreases in flood wave celerity and hydrograph kurtosis, skewness and rate of rise and increases in peak‐to‐peak travel time, flood peak attenuation and flood wave attenuation index. The most significant changes have occurred since the 1980s with the most noticeable effects on the behaviour of high flows (around the bankfull stage). However, there is not a consistent direction change in these indicators in all places. We use these findings to categorise the flow mitigation signal of the study regions and catchments and discuss the implication for river and flood management.

The persistence of snow on the ground affects the shape of streamflow hydrographs over space and time: a continental-scale analysis

Snow persistence (SP) is a widely available remotely-sensed measure of snowpack accumulation and ablation, reflecting the duration of snow presence on the ground in a given year. Available local-scale studies showed that SP is associated with the average magnitude of streamflow. However, despite the intuitive relationship between SP and catchment storage/release functioning, the spatial and temporal links between the persistence of snow on the ground and the shape and functionality of streamflow hydrographs were not studied empirically and were not generalized to diverse climatic settings. This study empirically explores the spatial and temporal links that SP has with measures of hydrograph shape and variability during low-flow and high-flow conditions across continent-wide gradients of aridity and seasonality. In arid in-phase and wet out-of-phase climates, higher SP is spatially associated with a damper (i.e., less flashy) streamflow hydrograph during low-flow and high-flow conditions. This is shown by a larger ratio of baseflow to average flow, a larger ratio of extreme low-flow to average flow, lower low-flow variability, and lower high-flow variability. While SP is spatially associated with a damped hydrograph in both arid/in-phase and wet/out-of-phase climates, this effect is stronger in the former region. For example, the size of the nonlinear impact of SP on reducing low-flow and high-flow variabilities is larger in arid in-phase climates (−7.64, −3.44, respectively) than in wet out-of-phase climates (−4.34, −2.02, respectively). Temporal analyses for “typical snow-rich” catchments show that years with relatively higher SP may lead to relatively flashier streamflow hydrographs, with lower baseflow indices, lower ratios of extreme low-flow to average flow, higher ratios of extreme high-flow to average flow and higher high-flow variability. Such results 1) demonstrate the utility of SP as a globally available descriptor of streamflow hydrograph shape and variability in a wide diversity of climatic conditions, 2) highlight that climate-driven snow loss may lead to substantial changes to hydrograph form and functionality, and 3) indicate that space-time symmetry may not be a valid assumption in hydrology.

Development and evaluation of the channel routing model and parameters within the National Water Model

AbstractThe National Water Model (NWM) was deployed by the National Oceanic and Atmospheric Administration to simulate operational forecasts of hydrologic states across the continental United States. This paper describes the geospatial river network (“hydro‐fabric”), physics, and parameters of the NWM, elucidating the challenges of extrapolating parameters a large scale with limited observations. A set of regression‐based channel geometry parameters are evaluated for a subset of the 2.7 million NWM reaches, and the riverine compound channel scheme is described. Based on the results from regional streamflow experiments within the broader NWM context, the compound channel reduced the root mean squared error by 2% and improved median Nash–Sutcliffe efficiency by 16% compared with a non‐compound formulation. Peak event analysis from 910 peak flow events across 26 basins matched from the US Flash Flood Observation Database revealed that the mean timing error is 3 h lagged behind the observations. The routing time step was also tested, for 5‐min (default, operational setting) and 1‐h increments. The model was computationally stable and able to convey the flood peaks, although the hydrograph shape and peak timing were altered.

Generation of type-specific synthetic design flood hydrographs

ABSTRACT Design flood hydrographs are essential for hydraulic engineering and flood risk assessment. Typically, they are derived from quantiles estimated from flood frequency analyses of flood peaks. This usually does not take into account that one and the same flood peak can occur with different volumes and various hydrograph shapes due to different meteorological and catchment-specific processes. In view of this issue, we generate design flood scenarios, i.e. type-specific design hydrographs are created for a flood quantile of a certain return period based on type-specific flood volumes and hydrograph shapes. Additionally, the probability of each flood type is derived. Different probability density functions are used as templates for generating different hydrograph shapes by varying their parameters. An ensemble of synthetic design flood hydrographs is obtained, which allows the consideration of best- and worst-case scenarios. A case study for 162 catchments in Germany demonstrates the applicability of the approach.

River flood routing using artificial neural networks

Flood routing through rivers is very important for the design of protection measures. Sensitivity analysis was performed using HEC-RAS 1D to show the effect of Manning’s coefficient, reach length, peak inflow hydrograph, base time of the inflow hydrograph, channel bed slope, channel side slope, and channel bed width on the peak discharge and the shape of the downstream hydrograph. About 1,600 synthetic realizations were generated using HEC-RAS 1D to train, validate, and test the artificial neural networks (ANNs). The network was trained to estimate the peak discharge from river reach and inflow hydrograph parameters. Testing this network shows a very good agreement between the predicted and actual peak discharges. Another network is trained to predict the whole outflow hydrograph from inflow hydrograph for a certain reach. Hence, flood routing through rivers can be performed using ANNs easily, accurately and quickly.

POLA DISTRIBUSI HUJAN DAN PENGARUHNYA TERHADAP PENYIMPANGAN HIDROGRAF BANJIR DAS JANGKOK

Determining the rain pattern distribution is believed to affect the design flood. Providing hourly rainfall observations can be obtained, and the distribution patterns can be determined. However, in areas that do not have one, flood discharge calculations are carried out by distributing daily rainfall using empirical methods. This study distributes daily rainfall into hourly rainfall in the Jangkok watershed using ABM and Mononobe methods and calculates the flood hydrograph using the Nakayasu model. The synthetic unit hydrograph obtained from the actual rain distribution is then compared with the unit hydrograph generated from the empirical distribution, then the deviation is measured. The results showed that the Mononobe and ABM rain distribution gave the same hydrograph shape as the observation one, only for 2 hours of rain duration. As long as the rain lasts 3-6 hours, the peak flood discharge (Qp) tends to be lower, 7-20% for Mononobe and 2-7% for ABM. Then the Qp becomes overestimated (5-12%) for a rain duration of >6 hours. Mononobe gives the same flood peak time (Tp) as observation rain for 1-5 hours rain duration and becomes 1 hour earlier for >6 hours. In contrast, ABM produces peak times 1 hour longer than hydrographs with observation rain. Generally, the two methods provide a more significant hydrograph deviation for a longer duration of rain.

Effects of Within-Storm Variability on Allochthonous Flash Flooding: A Synthetic Study

Rainfall is a spatiotemporally variated process and one of the key elements to accurately capture both catchment runoff response and floodplain extents. Flash floods are the result of intense rainfall, typically associated to highly variable rain in both space and time, such as convective storms. In this work, the extent within-storm variability affects runoff and flooding is explored. The Kan catchment (Tehran, Iran) is used as base topography for the simulations. The allochthonous nature of floods in the catchment and how they interact with the effects of storm variability are further investigated. For this, 300 synthetic rainfall signals with different hyetograph variabilities are generated and imposed on a 1D/2D hydrodynamic model. Additionally, a set of simulations with different levels of spatial variability are performed. The results suggest that temporal and spatial variability affect the runoff response in different degrees. Peak discharge and hydrograph shapes, as well as flooded areas, are affected. The effect of storm temporal variability is shown to be significantly higher than storm spatial variability and storm properties such as return period, duration, and volume. Further on the influence of storm spatiotemporal variability on stream discharge and flood response is seen to be strongly dependent on the location within the drainage network at which it is assessed.

Rainfall–streamflow response times for diverse upland UK micro-basins: quantifying hydrographs to identify the nonlinearity of storm response

AbstractWhile it is known that antecedent conditions and rainfall profiles contribute to the nonlinearity of streamflow response and that hydrograph shape can be dependent on the nature of rainfall inputs, how antecedent conditions (with similar rainfall inputs) impact hydrograph shape is less known. Here, a data-based mechanistic (DBM) approach is applied to quantify hydrograph shape, in terms of timing and volume, for the purposes of comparing hydrographs across 17 micro-basins at selected localities in upland UK over a 4-year period. The analysis demonstrates the nonlinearity of storm response for small catchments and revealed that with low antecedent conditions and/or small rainfall inputs there was a high variance in hydrograph shape quantifiers and that these variances decrease (at rates micro-basin dependent) as the micro-basins became wetter or as the storms increased in size, potentially converging to a more stable response.

Symptoms of Performance Degradation During Multi‐Annual Drought: A Large‐Sample, Multi‐Model Study

AbstractHydrologic models are essential tools to understand and plan for the effect of changing climates; however, they underperform in transitory climate conditions. Existing research identifies models' inadequacy to perform during prolonged drought, but falls short on pinpointing which specific aspects of model performance are affected. We study five conceptual rainfall‐runoff models and their performance in 155 Australian catchments which recently experienced a 13‐year long drought. We use a wide range of performance metrics and a methodology based on ranked differences to a benchmark to fairly compare levels of degradation across metrics and periods. We show model performance degrading extensively during and after the drought, largely driven by overestimation of flow. Representation of shape and variability of hydrograph and flow‐duration curve are more resilient to the prolonged dry climate and rarely more degraded during the multi‐annual drought that on isolated dry years in the pre‐drought record. Conversely, volumetric error suffers from significant exacerbation over the multiple subsequent dry years. This indicates that catchment retention times and rates of storage depletion storage are significantly less affected by the drought than amounts of streamflow produced, pointing to a mismatch between reduction of influxes and out‐fluxes during the drought. We also identify a deficiency of models to delay and remove flow before it reaches the stream and keep track of moisture deficits over multiple dry seasons. By promoting rigorous investigation of models' shortcomings, we hope to foster the development of more robust model structures and/or calibration frameworks to improve applicability within climate change scenarios.

Artificial viscosity technique for direct runoff calculation

In this study, an Artificial Viscosity (AV) technique was developed for direct runoff calculation, instead of using the conventional Synthetic Unit Hydrograph (SUH) methods. We solved the Shallow Water Equations (SWE) with second-order accurate Godunov finite volume model and fourth-order Runge-Kutta explicit scheme. The AV technique was devised with a Laplacian and a biharmonic operator, and employed to solve the convective terms of the SWE. This technique was applied to rainfall-runoff laboratory cases with the measured rainfall and observed direct runoff values. For comparison purpose, several SUH methods commonly used such as SCS, Snyder, GAMA-1, ITB-1, ITB-2, and Nakayasu, were also used to compute the direct runoff values. The results showed that the AV technique could predict three parameters (i.e., peak discharge, time-to-peak, and shape of the direct runoff hydrograph) accurately. Meanwhile, significant discrepancies were shown by the SUH methods in estimating such parameters, thus indicating that the SWE modeling with an AV technique is significantly more accurate than the SUH methods in predicting the direct runoff hydrograph. This study shows an interesting example of how modern numerical computations can improve the runoff prediction and may be included as a standard technique for runoff computation in the future.