Assessing “Urban Karst” Effects from Groundwater–Storm Sewer System Interaction in a Till Aquitard

Storm sewer systems and their associated utility trenches may strongly influence the effects of urbanization on a groundwater system. This study was undertaken to identify the causes of district-wide basement infiltration in an aquitard system. It comprised widespread continuous monitoring of utility trench wells and dye tracing from storm sewer system exfiltration tests. The results indicate that a major effect of urbanization on shallow groundwater is related to storm sewer system exfiltration, which is marked by a characteristic pattern of head variations in the aquitard unrelated to distributed surface infiltration. The aquitard constrains flow from storm sewer system exfiltration to the utility trench, creating an urban flow path for groundwater discharge. Temporary buildup of water levels in the utility trench drives relatively high-velocity flow through the permeable sewer bedding material of the utility trench to a separate foundation drainage collector system, ultimately causing a severe “urban karst” effect that produces system surcharging and widespread basement water infiltration. The main conditions causing the “urban karst” are the large hydraulic conductivity ratio between the utility trench material and the aquitard, and the shallow depth and low gradient of the storm sewer system imposed by a very flat drainage basin.     

Uncertainty in sewer sediment deposit modelling: Detailed vs simplified modelling approaches

The paper presents the results of a study in which the uncertainty levels associated with a detailed and a simplified/parsimonious sewer sediment modelling approach have been compared. The detailed approach used an Infoworks CS sewer network model combined with a user developed sediment transport code and the simplified approach used a conceptual sewer flow and quality model. The two approaches have been applied to a single case study sewer network and the simulation results compared. The case study was selected as moderate storm events had occurred during a 2 year rainfall and sewer flow monitoring period. Flooding had been observed and this was thought to be caused by significant solids accumulation in the sewer network. As a result sediment deposit measurements were carried out over a 6 month period. Model simulations were made of this period and predictions obtained of sediment deposit location and depth. The uncertainty analysis of both modelling approaches was carried out using Monte Carlo based computational methods. This was a limitation for the detailed approach with regards to computational time. Use of the simplified model was not constrained by this issue and so a more conventional assessment of the uncertainty was possible. The simplified approach, due to its structure, only provided a temporal estimate of uncertainty at the final section of the catchment. The detailed approach enabled an assessment of uncertainty at an individual pipe scale but only at the end of the simulation period. A comparison of the uncertainty estimations from both methods at the final section of the catchment and the end of the simulation period indicated comparable values of predicted uncertainty. Therefore a complementary use of both approaches would allow reasonably comparable estimations of levels of uncertainty at both a spatial and temporal scale. The use of such modelling approaches may provide a useful decision-making tool for sewer system management.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

The design storm concept in flood control design and planning

The design storm approach, where the subject criterion variable is evaluated by using a synthetic storm pattern composed of identical return frequencies of storm pattern input, is shown to be an effective approximation to a considerably more complex probabilistic model. The single area unit hydrograph technique is shown to be an accurate mathematical model of a highly discretized catchment with linear routing for channel flow approximation, and effective rainfalls in subareas which are linear with respect to effective rainfall output for a selected “loss” function. The use of a simple “loss” function which directly equates to the distribution of rainfall depth-duration statistics (such as a constant fraction of rainfall, or a ?-index model) is shown to allow the pooling of data and thereby provide a higher level of statistical significance (in estimating T-year outputs for a hydrologic criterion variable) than use of an arbitrary “loss” function. The above design storm unit hydrograph approach is shown to provide the T-year estimate of a criterion variable when using rainfall data to estimate runoff.     

Evaluation of the parameters affecting the roughness coefficient of sewer pipes with rigid and loose boundary conditions via kernel based approaches

One of the important issues in water transport and sewer systems is determining the flow resistance and roughness coefficient.An accurate estimation of the roughness coefficient is a substantial issue in the design and operation of hydraulic structures such as sewer pipes,the calculation of water depth and flow velocity,and the accurate characterization of energy losses.The current study,applies two kernel based approaches [Support Vector Machine(SVM) and Gaussian Process Regression(GPR)] to develop roughness coefficient models for sewer pipes.In the modeling process,two types of sewer bed conditions were considered:loose bed and rigid bed.In order to develop the models,different input combinations were considered under three scenarios(Scenario 1:based on hydraulic characteristics,Scenarios2 and 3:based on hydraulic and sediment characteristics with and without considering sediment concentration as input).The results proved the capability of the kernel based approaches in prediction of the roughness coefficient and it was found that for prediction of this parameter in sewer pipes Scenario 3 performed better than Scenarios 1 and 2.Also,the sensitivity analysis results showed that Dgr(Dimensionless particle number) for a rigid bed and w_b/y(ratio of deposited bed width,w_b,to flow depth,y) for a loose bed had the most significant impact on the modeling process.     

Simulation and design of pumped storm water systems

Flow and damage in pumped storm water systems are simulated, taking into consideration surcharged flow, basement and street storage. The simulation model is used in a design model to obtain practical global optimum values for the sewer depth, slope and diameter as well as sump and pump sizes.     

From marginal to conditional probability functions of parameters in a conceptual rainfall–runoff model: an event-based approach

ABSTRACT

A parameter estimation strategy for a conceptual rainfall–runoff (CRR) model applied to a storm sewer system in an urban catchment (Chassieu, Lyon, France) is proposed on the basis of event-by-event Bayesian local calibrations. The marginal distribution formed by locally-estimated parameters is divided into conditional functions, clustering the event-based parameters based on their transferability to similar rainfall events. The conditional functions showed to be consistent with an observed bimodality in the marginal representation, reflecting two different hydrological conditions mainly related to the magnitude of the rainfall intensities (high or low). The improvements achieved by expressing the parameter probability functions into a conditional form are shown in terms of accuracy (Nash-Sutcliffe efficiency criterion), precision (average relative interval length) and reliability (percentage of coverage) for simulating flow rate in 255 and 110 calibration/verification events.     

Investigating hydrologic responses to spatio‐temporal characteristics of storms using a Dynamic Moving Storm generator

A synthetic storm generator—Dynamic Moving Storm (DMS)—is developed in this study to represent spatio‐temporal variabilities of rainfall and storm movement in synthetic storms. Using an urban watershed as the testbed, the authors investigate the hydrologic responses to the DMS parameters and their interactions. In order to reveal the complex nature of rainfall–run‐off processes, previously simplified assumptions are relaxed in this study regarding (a) temporal variability of rainfall intensity and (b) time‐invariant flow velocity in channel routing. The results of this study demonstrate the significant contribution of storm moving velocity to the variation of peak discharge based on a global sensitivity analysis. Furthermore, a pairwise sensitivity analysis is conducted to elucidate not only the patterns in individual contributions from parameters to hydrologic responses but also their interactions with storm moving velocity. The intricacies of peak discharges resulting from sensitivity analyses are then dissected into independent hydrologic metrics, that is, run‐off volume and standard deviation of run‐off timings, for deeper insights. It is confirmed that peak discharge is increased when storms travel downstream along the main channel at the speed that corresponds to a temporal superposition of run‐off. Spatial concentration of catchment rainfall is found to be a critical linkage through which characteristics of moving storms affect peak discharges. In addition, altering peak timing of rainfall intensity in conjunction with storm movement results in varied storm core locations in the channel network, which further changes the flow attenuation effects from channel routing. For future directions, the DMS generator will be embedded in a stochastic modelling framework and applied in rainfall/flow frequency analysis.     

Integration of transfer function model and back propagation neural network for forecasting storm sewer flow in Taipei metropolis

Typhoons and storms have often brought heavy rainfalls and induced floods that have frequently caused severe damage and loss of life in Taiwan. Our ability to predict sewer discharge and forecast floods in advance during storm seasons plays an important role in flood warning and flood hazard mitigation. In this paper, we develop an integrated model (TFMBPN) for forecasting sewer discharge that combines two traditional models: a transfer function model and a back propagation neural network. We evaluated the integrated model and the two traditional models by applying them to a sewer system of Taipei metropolis during three past typhoon events (NARI, SINLAKU, and NAKR). The performances of the models were evaluated by using predictions of a total of 6 h of sewer flow stages, and six different evaluation indices of the predictions. Finally, an overall performance index was determined to assess the overall performance of each model. Based on these evaluation indices, our analysis shows that TFMBNP yields accurate results that surpass the two traditional models. Thus, TFMBNP appears to be a promising tool for flood forecasting for the Taipei metropolis sewer system. For publication in Stochastic Environmental Research and Risk Analysis.     

Parameterization and multi‐criteria calibration of a distributed storm flow model applied to a Mediterranean agricultural catchment

The principal challenge in the parameterization of storm flow models for agricultural catchments with an artificial drainage network and fields with different degrees of tillage lies in the parsimonious definition of distributed model parameters in a way that reduces the number of calibration parameters to a justifiable minimum. This paper presents a comprehensive case study for the parameter estimation of a distributed storm flow model applied to an agricultural catchment (0.91 km²) in the Mediterranean region. Model parameterization was combined with procedures for multi‐criteria, multi‐storm calibration, where we automatically calibrated three parameters related to flow velocity and infiltration, and compared single and multi‐storm criteria that are based on discharge volume, peak flow, and the Nash–Sutcliffe coefficient. Multi‐storm calibration yielded a set of parameter values for the simulation batch with best multi‐storm overall performance, which are close to the median values in the pre‐calibration of individual storms. Our results suggest that flow velocities and proportionality of the channel infiltration rate do not vary significantly over the course of 11 years. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of storm direction and duration on infiltrating planar flow with partial coverage

Using kinematic wave equations analytical solutions are derived for flow resulting from a storm moving either up or down an infiltrating plane but not fully covering it. By comparing the flow resulting from this storm with that from a stationary storm of the same duration the influence of storm duration, direction and velocity is investigated. It is found that the direction of storm movement, duration and velocity of storms, as well as basin infiltration, have a pronounced effect on the discharge hydrograph. Copyright © 2004 John Wiley & Sons, Ltd.     

The effect of dispersion on the stability of density-driven flows in saturated homogeneous porous media

This work is a continuation of Musuuza et al. [37] in which a stability criterion for density-driven flow in a saturated homogeneous medium was derived. The criterion predicted the stability of a system as a function of the density and viscosity contrasts, the magnitude of the flow velocity and the concentration gradients for flow aligned orthogonal to gravity. It could not accurately predict stability transition with increasing velocity, a failure we attributed to dispersion effects that were not included. Small-scale dispersion and molecular diffusion, the main stabilising mechanisms in homogeneous media can stabilise flow if the instability wavelengths are below a certain cutoff. The width of the mixing zone is also central in controlling the range of wavelengths that persist into fingers. We propose a method of quantifying the cutoff wavelength and the width of the mixing zone, which are incorporated into the earlier criterion as constituents of the dispersive part. The earlier criterion is reformulated in terms of the Rayleigh number and with the dispersive part added, we attempt to predict the number of fingers formed which is directly related to the physical stability of the system. The inclusion of the Rayleigh number and dispersion into a single stability criterion provides new insight in the way dispersion affects vertical flow systems. Stability numbers computed with the new criterion are in reasonable agreement with numerical simulations for a range of physical variables. The numerical computations are performed with the software package d³f, which uses the cell-centred finite volume and the implicit Euler methods for the spatial and temporal discretisations, respectively. The admission of the density and dispersivities as inputs into the criterion makes it usable in practical problems.     

Continuous fluorescence assessment of organic matter variability on the Bournbrook River,Birmingham, UK

Continuous monitoring of dissolved organic matter (DOM) character and concentration at hourly resolution is rare, despite the importance of analysing organic matter variability at high‐temporal resolution to evaluate river carbon budgeting, river water health by detecting episodic pollution and to determine short‐term variations in chemical and ecological function. The authors report a 2‐week experiment performed on DOM sampled from Bournbrook, Birmingham, UK, an urban river for which spectrophotometric (fluorescence, absorbance), physiochemical (dissolved organic carbon [DOC], electrical conductivity, pH) and isotopic (D/H) parameters have been measured at hourly frequency. Our results show that the river had sub‐daily variations in both organic matter concentration and characteristics. In particular, after relatively high‐magnitude precipitation events, organic carbon concentration increased, with an associated increase in intensity of both humic‐like and tryptophan‐like fluorescence. D/H isotopic ratio demonstrates different hydrological responses to different rainfall events, and organic matter character reflects this difference. Events with precipitation < 2 mm typically yielded isotopically heavy water with relatively hydrophilic DOM and relatively low specific absorbance. Events with precipitation > 2 mm had isotopically lighter water with higher specific absorbance and a decrease in the proportion of microbially derived to humic‐like fluorescence. In our heavily urbanized catchment, we interpret these signals as one where riverine DOM is dominated by storm sewer‐derived ‘old’ organic matter at low‐rainfall amounts and a mixed signal at high‐precipitation amounts where ‘event’ surface runoff‐derived organic matter dominate during storm sewer and combined sewer overflow routed DOM. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of the duration and direction of storm movement on infiltrating planar flow with full areal coverage

Using kinematic wave equations, analytical solutions are derived for flow due to a storm moving up or down an infiltrating plane and covering it completely. The storm duration is assumed in two ways. First, the plane is covered everywhere for the same duration by the storm. Second, the plane is covered in a linearly decreasing manner from the beginning of its coverage of the plane to the other end of the plane. By comparing the flow due to this storm with that due to a stationary storm of the same duration, the influence of storm duration, direction and velocity on flow hydrograph is investigated. It is found that storm movement has a pronounced effect on runoff hydrograph. Copyright © 2002 John Wiley & Sons, Ltd.     

Comparison of catchment and network delineation approaches in complex suburban environments: application to the Chaudanne catchment,France

Suburban areas are subject to strong anthropogenic modifications, which can influence hydrological processes. Sewer systems, ditches, sewer overflow devices and retention basins are introduced and large surface areas are sealed off. The knowledge of accurate flow paths and watershed boundaries in these suburban areas is important for storm water management, hydrological modelling and hydrological data analysis. This study proposes a new method for the determination of the drainage network based on time efficient field investigations and integration of sewer system maps into the drainage network for small catchments of up to 10 km². A new method is also proposed for the delineation of subcatchments and thus the catchment area. The subcatchments are delineated using a combination of an object‐oriented approach in the urban zone and geographical information system–based terrain analysis with flow direction forcing in the rural zone. The method is applied to the Chaudanne catchment, which belongs to the Yzeron river network and is located in the suburban area of Lyon, France. The resulting subcatchment map gives information about subcatchment response and contribution. The method is compared with six other automatic catchment delineation methods based on stream burning, flow direction forcing and calculation of subcatchments for inlet points. None of the automatic methods could correctly represent the catchment area and flow paths observed in the field. The watershed area calculated with these methods differs by as much as 25% from the area computed with the new method. Copyright © 2012 John Wiley & Sons, Ltd.     

The use of agent based modelling techniques in hydrology: determining the spatial and temporal origin of channel flow in semi‐arid catchments

Information on the spatial and temporal origin of runoff entering the channel during a storm event would be valuable in understanding the physical dynamics of catchment hydrology; this knowledge could be used to help design flood defences and diffuse pollution mitigation strategies. The majority of distributed hydrological models give information on the amount of flow leaving a catchment and the pattern of fluxes within the catchment. However, these models do not give any precise information on the origin of runoff within the catchment. The spatial and temporal distribution of runoff sources is particularly intricate in semi‐arid catchments, where there are complex interactions between runoff generation, transmission and re‐infiltration over short temporal scales. Agents are software components that are capable of moving through and responding to their local environment. In this application, the agents trace the path taken by water through the catchment. They have information on their local environment and on the basis of this information make decisions on where to move. Within a given model iteration, the agents are able to stay in the current cell, infiltrate into the soil or flow into a neighbouring cell. The information on the current state of the hydrological environment is provided by the environment generator. In this application, the Connectivity of Runoff Model (CRUM) has been used to generate the environment. CRUM is a physically based, distributed, dynamic hydrology model, which considers the hydrological processes relevant for a semi‐arid environment at the temporal scale of a single storm event. During the storm event, agents are introduced into the model across the catchment; they trace the flows of water and store information on the flow pathways. Therefore, this modelling approach is capable of giving a novel picture of the temporal and spatial dynamics of flow generation and transmission during a storm event. This is possible by extracting the pathways taken by the agents at different time slices during the storm. The agent based modelling approach has been applied to two small catchments in South East Spain. The modelling approach showed that the two catchments responded differently to the same rainfall event due to the differences in the runoff generation and overland flow connectivity between the two catchments. The model also showed that the time of travel to the nearest flow concentration is extremely important for determining the connectivity of a point in the landscape with the catchment outflow. Copyright © 2007 John Wiley & Sons, Ltd.     

Sewers as a Source and Sink of Chlorinated‐Solvent Groundwater Contamination,Marine Corps Recruit Depot,Parris Island,South Carolina

Groundwater contamination by tetrachloroethene and its dechlorination products is present in two partially intermingled plumes in the surficial aquifer near a former dry‐cleaning facility at Site 45, Marine Corps Recruit Depot, Parris Island, South Carolina. The northern plume originates from the vicinity of former above‐ground storage tanks. Free‐phase tetrachloroethene from activities in this area entered the groundwater. The southern plume originates at a nearby, new dry‐cleaning facility, but probably was the result of contamination released to the aquifer from a leaking sanitary sewer line from the former dry‐cleaning facility. Discharge of dissolved groundwater contamination is primarily to leaking storm sewers below the water table. The strong influence of sanitary sewers on source distribution and of storm sewers on plume orientation and discharge at this site indicates that groundwater‐contamination investigators should consider the potential influence of sewer systems at their sites.     

Inference model derivation with a pattern analysis for predicting the risk of microbial pollution in a sewer system

Developing a mathematical model for predicting fecal coliform bacteria concentration is very important because it can provide a basis for water quality management decisions that can minimize microbial pollution risk to the public. This paper introduces a hybrid modeling methodology which is a combined use of a neural network-based pattern analysis and an evolutionary process model induction system. The neural network-based pattern analysis technique is applied to extract knowledge on inter-relationships between fecal coliform concentrations and other measurable variables in a sewer system. Based on the result of neural network-based pattern analysis, an evolutionary process model induction system is used to derive mathematical inference models that can predict fecal coliform bacteria concentration from easily measurable variables instead of directly measuring fecal coliform bacteria concentration in a sewer system. The neural network-based pattern analysis extracts that temperature and ammonia concentration are the most important driving forces leading to an increase in fecal coliform bacteria concentration in the sewer system at Paraparaumu City, New Zealand. Fecal coliform bacteria concentration is also positively correlated with dissolved phosphorus and inversely with flow rate. The multivariate inference models that are able to predict fecal coliform bacteria concentration are successfully derived as functions of flow rate, temperature, ammonia, and dissolved phosphorus in the form of understandable mathematical formulae using the evolutionary process model induction system, even if a priori mathematical knowledge of the dynamic nature of fecal coliform bacteria is poor. The multivariate inference models evolved by the evolutionary process model induction system produce a slightly better performance than the multi-layer perceptron neural network model.     

Aggregate settling velocity of combined sewage overflow

In order to better model the fate and transport of material in a combined sewer overflow (CSO), hydrologic and geochemical measurements were made during a storm event in Flushing Bay, NY, USA. Particle size, total suspended solids concentrations, and CSO solid density are used to calculate the settling velocity of aggregates of the CSO material. Roughly half of the material that is discharged from the outfall sinks rapidly to the bottom, while the other half forms a turbid surface plume in fresher water. The dry solid density of the CSO material was 1589 kg m(-3) and the volume fraction of solid material to total aggregate volume ranged from 0.1 to 1.0. The settling velocity estimates of aggregated material in the surface plume ranged from about 0.5 to 9.0 mm s(-1) with a median settling velocity of 2.9 mm s(-1). A strong inverse relationship between the size of the aggregates and the aggregate density resulted in a fractal dimension of the aggregates of 2.3.     

Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal