A study on the complex processes of hydrogeological transformation of recharge into discharge within a karst system is described. Given the absence of boreholes or points of direct observation, the ...applied methodology is based on study of the natural response of the main outlet of the karst aquifer, considered here as a black-box. Time series analysis and recession curve analysis are applied to daily time series of rainfall as input, and spring discharge as the output of the aquifer. Daily variations of temperature and electrical conductivity of spring water are considered as well. The present study, an integrated research endeavor, is the first quantitative approach toward understanding the hydrogeological functioning of a large and not-well-known karst aquifer and improving existing knowledge about its internal physical characteristics. The case study is the Degracias-Sicó karst aquifer, in Portugal’s central-western region, which is relatively unknown from the hydrogeological perspective. The response of the karst aquifer’s main spring shows an important component of quick flow after heavy rainfall events and then a very slow recession. This bimodal behavior is proof of the significant spatial complexity and heterogeneity of the karst aquifer, marked by the presence of a conduit-dominated flow throughout a connected conduit network system and the dominance of base flow that confirms the aquifer’s large storage capacity. The results also demonstrate a significant component of delayed flow.
•The StreamGEM model was calibrated using monthly nitrate data from four catchments.•Posterior parameter distributions were derived by Bayesian analysis with DREAMZS.•Predicted flow path water and ...nitrate fluxes were consistent with catchment geology.•Use of daily nitrate data in one catchment did not improve model predictions.•Four years of monthly samples was sufficient to characterise catchment dynamics.
The common practice of infrequent (e.g., monthly) stream water quality sampling for state of the environment monitoring may, when combined with high resolution stream flow data, provide sufficient information to accurately characterise the dominant nutrient transfer pathways and predict annual catchment yields. In the proposed approach, we use the spatially lumped catchment model StreamGEM to predict daily stream flow and nitrate concentration (mg L−1 NO3-N) in four contrasting mesoscale headwater catchments based on four years of daily rainfall, potential evapotranspiration, and stream flow measurements, and monthly or daily nitrate concentrations. Posterior model parameter distributions were estimated using the Markov Chain Monte Carlo sampling code DREAMZS and a log-likelihood function assuming heteroscedastic, t-distributed residuals. Despite high uncertainty in some model parameters, the flow and nitrate calibration data was well reproduced across all catchments (Nash-Sutcliffe efficiency against Log transformed data, NSL, in the range 0.62–0.83 for daily flow and 0.17–0.88 for nitrate concentration). The slight increase in the size of the residuals for a separate validation period was considered acceptable (NSL in the range 0.60–0.89 for daily flow and 0.10–0.74 for nitrate concentration, excluding one data set with limited validation data). Proportions of flow and nitrate discharge attributed to near-surface, fast seasonal groundwater and slow deeper groundwater were consistent with expectations based on catchment geology. The results for the Weida Stream in Thuringia, Germany, using monthly as opposed to daily nitrate data were, for all intents and purposes, identical, suggesting that four years of monthly nitrate sampling provides sufficient information for calibration of the StreamGEM model and prediction of catchment dynamics. This study highlights the remarkable effectiveness of process based, spatially lumped modelling with commonly available monthly stream sample data, to elucidate high resolution catchment function, when appropriate calibration methods are used that correctly handle the inherent uncertainties.
Abstract
The objective of this study was to experimentally evaluate the difference in the transport of uniform (5.17, 10.35, 14, 20.7 mm) and graded sediment (mixture of these rounded particles with ...equal weight proportions) under different unsteady flow hydrographs in a 12 m long, 0.5 m wide and deep glass-walled flume. There was a lag time between fractions and uniform particles, such that peaks of coarser and finer fraction particles occurred before and after the peak of uniform sediment with the same size, respectively. Comparison between uniform particles and fractions in graded sediment showed that the sediment transport rate of fine and coarse fractions was lower and higher than their counterpart uniform particles, respectively. Overall, the uniform particles demonstrated a clockwise hysteresis loop and graded sediment had a counterclockwise hysteresis loop. The mobility of coarser fractions increased during the rising limb of hydrograph, whereas the mobility of finer fractions increased during the falling limb. In general, the mobility of coarse fractions increased and that of fine fractions reduced. Result of transported sediment showed that average particle size collected in traps (Db50) was coarser than bed material (Ds50) on both limbs. The relative transport ratio for uniform and graded sediment is higher and lower than 1, respectively.
Abstract
Boyong River is one of many rivers originating from Mt. Merapi, flows across three autonomy administrative of Sleman Regency, Yogyakarta City, and Bantul Regency. The river experiences flood ...in the form of lava flow several times, and the 1994 and 2010 occurrences were considered the biggest ones along with the river history. In line with the rapid development of information and communication technology, efforts to develop the early warning system due to the Mt. Merapi disaster have been implemented by the Hydraulic Laboratory of Gadjah Mada University 2006. This paper presents the study results of Boyong River flow behavior by analyzing the data obtained from the monitoring system. The Gemawang Weir at Boyong River was selected as the river control point understudy; those include the catchment boundary, the catchment characteristics, and the hydraulic features. Monitoring equipment consists of an automatic water level recorder (AWLR), the flow visualization using a Brinno camera, and the hydrophone monitoring system. The flow hydrograph characteristics and its corresponding sediment transport rate are considered two parameters for identifying the flow behavior. The results show that the precursive and recession times of the flood hydrograph are about 1-3 hours and 3.5-5 hours, respectively.
Certain relevant variables for dam safety and downstream safety assessments are analyzed using a stochastic approach. In particular, a method to estimate quantiles of maximum outflow in a dam ...spillway and maximum water level reached in the reservoir during a flood event is presented. The hydrological system analyzed herein is a small mountain catchment in north Spain, whose main river is a tributary of Ebro river. The ancient Foradada dam is located in this catchment. This dam has no gates, so that flood routing operation results from simple consideration of fixed crest spillway hydraulics. In such case, both mentioned variables (maximum outflow and maximum reservoir water level) are basically derived variables that depend on flood hydrograph characteristics and the reservoir’s initial water level. A Monte Carlo approach is performed to generate very large samples of synthetic hydrographs and previous reservoir levels. The use of extreme value copulas allows the ensembles to preserve statistical properties of historical samples and the observed empirical correlations. Apart from the classical approach based on annual periods, the modelling strategy is also applied differentiating two subperiods or seasons (i.e., summer and winter). This allows to quantify the return period distortion introduced when seasonality is ignored in the statistical analysis of the two relevant variables selected for hydrological risk assessment. Results indicate significant deviations for return periods over 125 years. For the analyzed case study, ignoring seasonal statistics and trends, yields to maximum outflows underestimation of 18% for T = 500 years and 29% for T = 1000 years were obtained.
In the wake of increasing flood disasters, there is an increasing use of flood inundation models to assess risks and impacts at different temporal and spatial scales. Assessing the impacts of extreme ...climatic rainfall events will require developing design rainfall profiles to represent rainfall under different conditions. Rainfall profiles of different return periods were developed using the Flood Estimation Handbook (FEH) methodology for a small rural catchment of Scotland, to assess flood risks at a catchment scale. Rainfall induced runoff flows were estimated based on a set of catchment characteristics. The channel and floodplain flows were modelled using a two-dimensional hydrodynamic model-TUFLOW. The main channel was represented by a one-dimensional linear channel based on surveyed data and the floodplain topography, was represented by a digital terrain model based on Light Detection and Ranging (LiDAR). A range of hydrological events with different return periods are simulated. Results show that many residential houses and an extensive area of agricultural land are at risk of flooding from extreme events such as a 1 in 100 year flood.
Digital elevation model (DEM) resolution and the assigned threshold for river network delineation affect the results of rainfall-runoff models. In this study, the effects of these 2 issues on the ...extracted geomorphologic parameters of watersheds and the performance of a kinematic wave based model, called KW-GIUH, are investigated. The results show that by decreasing the DEM resolution at fixed thresholds, parameters such as subbasin mean slope and the number of streams decrease and the area of the ith order subbasins and the mean length of the overland flow increase. Moreover, the results indicate that the reduction of the DEM resolution at a fixed threshold causes the peak flow and hydrograph time base to decrease up to the cell size of 100 m and then, after experiencing a jump, again decrease with the increase of the cell size. According to the achieved results, above the threshold of 2%, the difference between the peak flows of different hydrographs at different resolutions is meaningful. The KW-GIUH sensitivity to DEM resolutions and thresholds is sharper in peak flow and then in hydrograph time base and time to peak. At a fixed threshold, the value of time to peak is independent of DEM resolution.
Physically based modeling approach has been widely developed in recent years for simulation of dam failure process resulting from overtopping flow. Due to the lack of field data, there exist few ...applications to natural quake dams with complex erosion mechanisms. This paper presents a physically based simulation of the failure process of the Tangjiashan Quake Dam formed as a result of the “May 12, 2008” Wenchuan earthquake in China. The one-dimensional model adopted features as cost saving but enables capturing the main characteristics of the failure process, where selective sediment transport and gravitational slope collapse are accounted for. The simulated flow hydrograph and breach progression process are generally in good agreement with the observed data. Unsteadiness and non-uniformity are found to be substantial characteristics of breach progression during the failure process of natural quake dams. Sensitivity analysis showed that the Manning resistance coefficient and the erodibility coefficient in Osman and Thorne’s (J Hydraul Eng 114(2):134–150,
1988
) model significantly influences the flow peak discharge but has less influence on its occurrence time, while the velocity lag coefficient associated with bed-load transport may affect the two breaching parameters substantially.
This study developed a genetic algorithm model to predict flow rates at sites receiving significant lateral inflow. It predicts flow rate at a downstream station from flow stage measured at upstream ...and downstream stations. For this purpose, it constructed two different models: First is analogous to the rating curve model (RCM) of Moramarco et al. Moramarco, M., Barbetta, S., Melone, F., Singh, V.P., 2005. Relating local stage and remote discharge with significant lateral inflow. J. Hydrologic Eng., ASCE, 10(1) and the second is based on summation of contributions from upstream station and lateral inflows using kinematic wave approximation. The model was applied to predict flow rates at three different gauging stations located on Tiber River, Upper Tiber River Basin, Italy. The model used average wave travel time for each river reach and obtained average set of parameter values for all the events observed in the same river reach. The GA model was calibrated, for each river reach and for each formulation, by three events and tested against three other events. The results showed that the GA model produced satisfactory results and it was superior over the most recently developed rating curve method. This study further analyzed the case where only water surface elevation data were used in the input vector to predict flow rates. The results showed that using elevation data produces satisfactory results. This has an implication for predicting flow rates at ungauged river sites since the surface elevation data can be obtained without needing the detailed geometry of river section which could change significantly during a flood.