RECHARGE TO ALLUVIAL VALLEY AQUIFERS FROM OVERBANK FLOW AND EXCESS INFILTRATION1

ABSTRACT: Recharge is an important parameter for models that simulate water and contaminant transport in unconfined aquifers. Unfortunately, measurements of actual recharge are not usually available causing recharge to be estimated or possibly added to the calibration procedure. In this study, differences between observed water-table elevations and water-table elevations simulated with a model based on the one-dimensional Boussinesq equation were used to identify both the timing and quantity of recharge to an alluvial valley aquifer. Observed water table elevations and river stage data were recorded during a five-year period from 1991 to 1995 at the Ohio Management Systems Evaluation Area located in south-central Ohio. Direct recharge attributed to overbank flow during and shortly after flood conditions accounted for 65 percent of the total recharge computed during the five-year study period. Recharge of excess infiltration to the aquifer was intermittent and occurred soon after large rainfall events and high river stage. Specification of constant recharge with time values in ground-water simulation models seems inappropriate for stream-aquifer systems given the strong influence of the river on water table elevations in these systems.     

An innovative modeling approach using Qual2K and HEC-RAS integration to assess the impact of tidal effect on River Water quality simulation

Water quality modeling has been shown to be a useful tool in strategic water quality management. The present study combines the Qual2K model with the HEC-RAS model to assess the water quality of a tidal river in northern Taiwan. The contaminant loadings of biochemical oxygen demand (BOD), ammonia nitrogen (NH₃-N), total phosphorus (TP), and sediment oxygen demand (SOD) are utilized in the Qual2K simulation. The HEC-RAS model is used to: (i) estimate the hydraulic constants for atmospheric re-aeration constant calculation; and (ii) calculate the water level profile variation to account for concentration changes as a result of tidal effect. The results show that HEC-RAS-assisted Qual2K simulations taking tidal effect into consideration produce water quality indices that, in general, agree with the monitoring data of the river. Comparisons of simulations with different combinations of contaminant loadings demonstrate that BOD is the most import contaminant. Streeter-Phelps simulation (in combination with HEC-RAS) is also performed for comparison, and the results show excellent agreement with the observed data. This paper is the first report of the innovative use of a combination of the HEC-RAS model and the Qual2K model (or Streeter-Phelps equation) to simulate water quality in a tidal river. The combination is shown to provide an alternative for water quality simulation of a tidal river when available dynamic-monitoring data are insufficient to assess the tidal effect of the river.     

Modeling Flow and Sediment Transport in a River System Using an Artificial Neural Network

A river system is a network of intertwining channels and tributaries, where interacting flow and sediment transport processes are complex and floods may frequently occur. In water resources management of a complex system of rivers, it is important that instream discharges and sediments being carried by streamflow are correctly predicted. In this study, a model for predicting flow and sediment transport in a river system is developed by incorporating flow and sediment mass conservation equations into an artificial neural network (ANN), using actual river network to design the ANN architecture, and expanding hydrological applications of the ANN modeling technique to sediment yield predictions. The ANN river system model is applied to modeling daily discharges and annual sediment discharges in the Jingjiang reach of the Yangtze River and Dongting Lake, China. By the comparison of calculated and observed data, it is demonstrated that the ANN technique is a powerful tool for real-time prediction of flow and sediment transport in a complex network of rivers. A significant advantage of applying the ANN technique to model flow and sediment phenomena is the minimum data requirements for topographical and morphometric information without significant loss of model accuracy. The methodology and results presented show that it is possible to integrate fundamental physical principles into a data-driven modeling technique and to use a natural system for ANN construction. This approach may increase model performance and interpretability while at the same time making the model more understandable to the engineering community.     

Hydrodynamic Modeling Improves Green River Reconnection with Floodplain Wetlands for Endangered Fish Species Recovery

We performed two‐dimensional (2D) hydrodynamic modeling to aid recovery of the endangered razorback sucker (Xyrauchen texanus) by reconnecting the Green River with its historic bottomland floodplain wetlands at Ouray National Wildlife Refuge, Utah. Reconnection allows spring flood flows to overtop the river levee every two to three years, and passively transport razorback sucker larvae to the wetlands to grow in critical habitat. This study includes (1) river hydrologic analysis, (2) simulation of a levee breach/weir, overtopping of river flood flows, and 2D flow through the wetlands using Hydrologic Engineering Center River Analysis System 2D, and (3) modeling flow and restoration scenarios. Indicators of hydrologic alteration were used to evaluate river flow metrics, in particular flood magnitudes, frequency, and duration. Results showed a target spring flow of 16,000 cfs (453 m³/s) and a levee breach elevation of 4,663 ft (1,421 m) amsl would result in a median flow >6,000 acre‐feet (7.4 million m³) over five days into the wetlands, which is adequate for razorback sucker larvae transport and rearing. Modeling of flow/restoration scenarios showed using gated water control structures and passive low‐water crossings between wetland units can provide adequate control of flow movement into and storage in multiple units. Levee breaching can be a relatively simple, cost‐effective method to reconnect rivers and historic floodplains, and hydrodynamic modeling is an important tool for analyzing and designing wetland reconnection.     

Conceptual Design of a System for Selecting Appropriate Groundwater Models in Groundwater Protection Programs

/ An effective groundwater protection program requires understanding of water flow and contaminant transport processes in the subsurface. Although many mathematical models have been developed to simulate the processes, few actually are used in groundwater protection programs due to the difficulties in data collection, model selection, and model implementation. This study presents a conceptual design of a GIS-supported model selection system that evaluates available data and mathematical models to facilitate groundwater protection programs. Steady-state groundwater and contaminant transport models applied in isotropic aquifers are placed into four classes to simulate conservative or nonconservative contaminant transports in simple or complex geohydrological conditions. After analyzing specific study objectives, available data, and model requirements, the proposed system selects a class of models that can be used in simulation and recommends any need for additional data collection. This study initiates an effort to integrate GIS, mathematical models, and expert knowledge in one system to promote the application of appropriate groundwater models. The new technology of GIS and digital data-base management makes it possible to develop such a system in practice.KEY WORDS: Groundwater models; Geographic information systems     

PARAMETER ESTIMATES FOR A GROUND WATER MODEL1

ABSTRACT: This paper presents a parameter sensitivity study of a two-dimensional flow and transport model of a contaminated site. Hydrogeological and site data from previous investigations were used for calibration. The USGS contaminant transport model (MOC) was used. After flow calibration to establish a reference model, parameters were varied to examine the effect each had on predictions of a contaminant plume. Hydrogeological parameters and a step size parameter were incrementally varied individually. Each result was compared to the reference model output to evaluate changes in concentration values and contaminant plume configuration. The study indicated that a generally predictable trend can be established for some parameters not affected by pumping or similar high stresses. Ranges were identified to relate concentration error or plume change to the amount of parameter error. Some parameter perturbations produced distorted model responses at high stress locations. Porosity and anisotropy were found to be the most influential of the model parameters studied on the plume predictions. (KEY TERMS: ground water hydrology; hydrogeology; pollution modeling; water quality; model calibration.)     

HYDRODYNAMIC MODELING OF THE CARSON RIVER AND LAHONTAN RESERVOIR,NEVADA1

ABSTRACT: The RIVMOD hydrodynamic model was used to route upstream flows through a 115 km section of the Carson River and Lahontan Reservoir, Nevada. RIVMOD results will later be used to predict sediment movement and ultimately to determine mercury transport within the river/reservoir system. Significant modifications to the model computer code were necessary to represent the narrow, steeply sloping rectangular channel and relatively shallow sloping floodplain of the Carson River and its confluence with the Lahontan Reservoir. These changes include expansion of the continuity and momentum equations to account for rapidly changing channel widths along with the characterization of a complex cross-sectional shape. This modified version of the RIVMOD model can handle shallower side slopes and much more severe flood flow simulations than the original version. A 0.25 km spatial increment was required in the zone of confluence between the river and reservoir. Model predictions show excellent agreement with observed downstream flow and reservoir stage for the entire 1986 water year, which includes one of the most severe flood events of recent record. (KEY TERMS: hydraulics; modeling; simulation; surface water hydrology.)     

Support of sustainable management of nitrogen contamination due to septic systems using numerical modeling methods

Nitrogen contamination is a serious concern to sustainable environmental management, and one important source of nitrogen contaminant is due to wastewater treatment using onsite sewage treatment and disposal systems (OSTDS, a.k.a., septic systems). This paper presents a study in which numerical modeling is used to support sustainable decision-making and management of nitrogen contamination by utilizing a recently developed GIS-based software, VZMOD, a Vadose Zone MODel for simulating nitrogen transformation and transport in vadose zone between drainfield of septic systems and water table. VZMOD is based on a physical model of unsaturated flow and nitrogen transformation and transport, and the model is solved numerically using the finite element methods. This is the major difference between VZMOD and other GIS-based software of nitrogen modeling. Using GIS techniques, VZMOD considers spatial variability of a number of hydrogeologic parameters such as hydraulic conductivity and porosity. A unique feature of VZMOD is that VZMOD addresses spatial variability of water table by using VZMOD together with ArcNLET, an ArcGIS-based software developed to simulate groundwater flow and nitrate load from septic systems to surface water bodies. VZMOD is designed to execute in different modes to be compatible with different levels of data availability in various management projects of nitrogen contamination. This paper presents an application of VZMOD at a neighborhood with hundreds of septic systems and heterogeneous hydraulic conductivity, porosity, and water table depth. The modeling results indicate that using septic systems at the considered neighborhood is unsustainable and more management means are necessary.     

Dynamic floodplain vegetation model development for the Kootenai River, USA

The Kootenai River floodplain in Idaho, USA, is nearly disconnected from its main channel due to levee construction and the operation of Libby Dam since 1972. The decreases in flood frequency and magnitude combined with the river modification have changed the physical processes and the dynamics of floodplain vegetation. This research describes the concept, methodologies and simulated results of the rule-based dynamic floodplain vegetation model "CASiMiR-vegetation" that is used to simulate the effect of hydrological alteration on vegetation dynamics. The vegetation dynamics are simulated based on existing theory but adapted to observed field data on the Kootenai River. The model simulates the changing vegetation patterns on an annual basis from an initial condition based on spatially distributed physical parameters such as shear stress, flood duration and height-over-base flow level. The model was calibrated and the robustness of the model was analyzed. The hydrodynamic (HD) models were used to simulate relevant physical processes representing historic, pre-dam, and post-dam conditions from different representative hydrographs. The general concept of the vegetation model is that a vegetation community will be recycled if the magnitude of a relevant physical parameter is greater than the threshold value for specific vegetation; otherwise, succession will take place toward maturation stage. The overall accuracy and agreement Kappa between simulated and field observed maps were low considering individual vegetation types in both calibration and validation areas. Overall accuracy (42% and 58%) and agreement between maps (0.18 and 0.27) increased notably when individual vegetation types were merged into vegetation phases in both calibration and validation areas, respectively. The area balance approach was used to analyze the proportion of area occupied by different vegetation phases in the simulated and observed map. The result showed the impact of the river modification and hydrological alteration on the floodplain vegetation. The spatially distributed vegetation model developed in this study is a step forward in modeling riparian vegetation succession and can be used for operational loss assessment, and river and floodplain restoration projects.     

OPTIMIZATION OF INTERMITTENT PUMPING SCHEDULES FOR AQUIFER REMEDIATION USING A GENETIC ALGORITHM1

ABSTRACT: Using a genetic algorithm (GA), optimal intermittent pumping schedules were established to simulate pump‐and‐treat remediation of a contaminated aquifer with known hydraulic limitations and a water miscible contaminant, located within the Duke Forest in Durham, North Carolina. The objectives of the optimization model were to minimize total costs, minimize health risks, and maximize the amount of contaminant removed from the aquifer. Stochastic ground water and contaminant transport models were required to provide estimates of contaminant concentrations at pumping wells. Optimization model simulations defined a tradeoff curve between the pumping cost and the amount of contaminant extracted from the aquifer. For this specific aquifer/miscible contaminant combination, the model simulations indicated that pump‐and‐treat remediation using intermittent pumping schedules for each pumping well produced significant reductions in predicted contaminant concentrations and associated health risks at a reasonable cost, after a remediation time of two years.     

Modeling Long‐Term Trends of Chlorinated Ethene Contamination at a Public Supply Well

A mass‐balance solute‐transport modeling approach was used to investigate the effects of dense nonaqueous phase liquid (DNAPL) volume, composition, and generation of daughter products on simulated and measured long‐term trends of chlorinated ethene (CE) concentrations at a public supply well. The model was built by telescoping a calibrated regional three‐dimensional MODFLOW model to the capture zone of a public supply well that has a history of CE contamination. The local model was then used to simulate the interactions between naturally occurring organic carbon that acts as an electron donor, and dissolved oxygen (DO), CEs, ferric iron, and sulfate that act as electron acceptors using the Sequential Electron Acceptor Model in three dimensions (SEAM3D) code. The modeling results indicate that asymmetry between rapidly rising and more gradual falling concentration trends over time suggests a DNAPL rather than a dissolved source of CEs. Peak concentrations of CEs are proportional to the volume and composition of the DNAPL source. The persistence of contamination, which can vary from a few years to centuries, is proportional to DNAPL volume, but is unaffected by DNAPL composition. These results show that monitoring CE concentrations in raw water produced by impacted public supply wells over time can provide useful information concerning the nature of contaminant sources and the likely future persistence of contamination.     

Investigating the Fate and Transport of Escherichia coli in the Charles River,Boston, Using High‐Resolution Observation and Modeling1

Abstract: The processes affecting the fate and transport of Escherichia coli in surface waters were investigated using high‐resolution observation and modeling. The concentration patterns in Boston’s Charles River were observed during four sampling events with a total of 757 samples, including two spatial surveys with two along‐river (1,500 m length) and three across‐river (600 m length) transects at approximately 25‐m intervals, and two temporal surveys at a fixed location (Community Boating) over seven days at hourly intervals. The data reveal significant spatial and temporal structure at scales not resolved by typical monitoring programs. A mechanistic, time‐variable, three‐dimensional coupled hydrodynamic and water quality model was developed using the ECOMSED and RCA modeling frameworks. The computational grid consists of 3,066 grid cells with average length dimension of 25 m. Forcing functions include upstream and downstream boundary conditions, Stony Brook, and Muddy River (major tributaries) combined sewer overflow (CSO) and non‐CSO discharge and wind. The model generally reproduces the observed spatial and temporal patterns. This includes the presence and absence of a plume in the study area under similar loading, but different hydrodynamic conditions caused by operation of the New Charles River Dam (downstream) and wind. The model also correctly predicts an episode of high concentrations at the time‐series station following seven days of no rainfall. The model has an overall root mean square error (RMSE) of 250 CFU/100 ml and an error rate (above or below the USEPA‐recommended single sample criteria value of 235 CFU/100 ml) of 9.4%. At the time series station, the model has an RMSE of 370 CFU/100 ml and an error rate of 15%.     

A STATISTICAL APPROACH TO ASSESS FACTORS AFFECTING WATER CHEMISTRY USING MONITORING DATA1

ABSTRACT: A method to partition the variation in concentrations of water chemistry parameters in a river is described. The approach consists of fitting a family of curves for each chemical parameter. Each curve indicates the response of the parameter to river flow for a particular time period or location. An analysis of covariance is then used to identify statistically significant differences between curves. Such differences result largely from two factors: (1) the discharge of effluents and (2) river flow-concentration relationships. The deviations from the fitted curves indicate month-to-month variations unrelated to river flow that are controlled by factors such as temperature-related seasonal patterns. Underlying statistical assumptions are discussed with respect to water chemistry data. The technique is applied to a data set consisting of monthly samples of 22 water chemistry parameters from the Sulphur River of Texas and Arkansas. Several patterns of response to river flow and to two effluent discharges were revealed.     

MANAGEMENT OF FLOOD CONTROL SUMPS AND POLLUTANT TRANSPORT1

ABSTRACT: Levee sump systems are used by many riverine communities for temporary storage of urban wet weather flows. The hydrologic performance and transport of stormwater pollutants in sump systems, however, have not been systematically studied. The objective of this paper is to present a case study to demonstrate development and application of a procedure for assessing the hydraulic performance of flood control sumps in an urban watershed. Two sumps of highly variable physical and hydraulic characteristics were selected for analysis. A hydrologic modeling package was used to estimate the flow hydrograph for each outfall as part of the flow balance for the sump. To validate these results, a water balance was used to estimate the total runoff using sump operational data. The hydrologic model calculations provide a satisfactory estimate of the total runoff and its time‐distribution to the sump. The model was then used to estimate pollutant loads to the sump and to the river. Although flow of stormwater through a sump system is regulated solely by flood‐control requirements, these sumps may function as sedimentation basins that provide purification of stormwater. A sample calculation of removals of several conventional pollutants in the target sumps using a mass balance approach is presented.     

Constructing an Interdisciplinary Flow Regime Recommendation1

Bartholow, John M., 2010. Constructing an Interdisciplinary Flow Regime Recommendation. Journal of the American Water Resources Association (JAWRA) 1-15. DOI: 10.1111/j.1752-1688.2010.00461.x Abstract: It is generally agreed that river rehabilitation most often relies on restoring a more natural flow regime, but credibly defining the desired regime can be problematic. I combined four distinct methods to develop and refine month-by-month and event-based flow recommendations to protect and partially restore the ecological integrity of the Cache la Poudre River through Fort Collins, Colorado. A statistical hydrologic approach was used to summarize the river’s natural flow regime and set provisional monthly flow targets at levels that were historically exceeded 75% of the time. These preliminary monthly targets were supplemented using results from three Poudre-specific disciplinary studies. A substrate maintenance flow model was used to better define the high flows needed to flush accumulated sediment from the river’s channel and help sustain the riparian zone in this snowmelt-dominated river. A hydraulic/habitat model and a water temperature model were both used to better define the minimum flows necessary to maintain a thriving cool water fishery. The result is a range of recommended monthly flows and daily flow guidance illustrating the advantage of combining a wide range of available disciplinary information, supplemented by judgment based on ecological principles and a general understanding of river ecosystems, in a highly altered, working river.     

Hydrogeological study of an anti-tank range

The Arnhem Anti-Tank Range (Canadian Forces Base [CFB] Valcartier, Canada, in operation since the 1970s) has been characterized, including the drilling, installation, and characterization of 25 wells and a ground-penetrating radar survey. The observed particular features of this site include highly variable flow velocities (from < 3 to 1200 m/yr) and transient flow regime in the regional aquifer below the contaminant source zone of the impact area, sharp flow direction shifts, discontinuous stratigraphy and a local perched aquifer. A transient ground water flow model permitted us to understand how the complex hydrogeological setting shapes contaminant transport in the regional aquifer. The model explains the highly variable energetic material (EM) concentrations measured in the plume with peaks associated to spring and to a lesser extent to fall recharge events. As a conclusion from this work, the authors suggest that the characterization of contaminant sources on slopes should extend over all seasons to be sure to detect potential transient flow conditions and variable contaminant concentrations.     

Aging effects on cadmium transport in undisturbed contaminated sandy soil columns

Limited information is available on the effects of contaminant aging (i.e., the contact time of Cd with the soil) on Cd transport in soils. We conducted displacement experiments in which indigenous Cd and freshly applied Cd were leached simultaneously from undisturbed samples of three Spodosol horizons. Sorption of Cd was described using Freundlich isotherms, whereas transport was described as a convection-dispersion process. Parameter optimization analysis using a mobile-immobile transport model applied to nonsorbing tracer displacement data showed that 16 to 22% of the water in the columns was immobile. The low dimensionless mass transfer coefficients in the mobile-immobile model were indicative of diffusion-limited transfer between mobile and immobile water, and hence physical nonequilibrium. A two-site kinetic sorption model could be fitted closely to breakthrough curves of the non-aged Cd for three soil horizons. No conclusive evidence was found that contaminant aging in soil affects cadmium transport. On the one hand, predictions of aged Cd leaching, using parameters estimated from displacement experiments with nonaged Cd, differed from those for the aged Cd in the E horizon. On the other hand, no meaningful differences in transport behavior between aged and non-aged Cd were found for the humus Bh and Bh/C horizons. The two-site kinetic rate coefficient alphac was found to depend on water flux, further indicating that mass transfer between sorption sites and the liquid is limited by diffusion rather than by kinetic sorption.     

SPATIALLY DISTRIBUTED MODELING OF STREAM FLOW DURING STORM EVENTS1

ABSTRACT: The purpose of this paper is to present a new approach for the spatially distributed modeling of water flow during storm events. Distributed modeling of flow during storm events is an important basis for any environmental modeling, including turbidity or sediment transport. During the initial phase of a rainstorm, surface runoff is the main contributor of flow. To provide the spatial components for distributed hydrological modeling a Geographic Information System (GIS) was used to map and visualize contributing areas around a stream channel. Stream segments were defined using the hydrologic response unit (HRU) concept. Lateral flows were derived from GIS output for each segment of the stream and at each time interval of the rain storm and were routed using the kinematic routing equation. This approach is new in hydrological modeling and can be used to enhance many existing simulations. The model is also unique in the fine time scale (i.e., intervals are on the order of minutes). Model results showed good correlation with measured discharge values; however, further studies of contributing area behavior, its relationship with soil types and slope categories, and the influence of watershed size are needed to improve model performance. This model will be used in the future as the basis to model turbidity in streams.     

Water quality modeling to determine minimum instream flow for fish survival in tidal rivers

The Hsintien Stream is one of the major branches of the Danshuei River system, which runs through the metropolitan capital city of Taipei, Taiwan and receives a large amount of wastewater. The dissolved oxygen concentration is generally low in the tidal portion of the Hsintien Stream. Hypoxia/anoxia occurs often, particularly during the low-flow period when the Feitsui Reservoir, Chingtan Dam and Chihtan Dam impound the freshwater for municipal water supply. Fish kills happen from time to time. This paper describes the application of a numerical hydrodynamic and water quality model to the Danshuei River system, with special attention to the tidal portion of the Hsintien Stream. The model is recalibrated with the prototype conditions of the year 2000. The hydrodynamic portion of the model is recalibrated with measured surface elevation and velocity at various stations in the river system. The water quality portion of the model is recalibrated with respect to the field data provided by Taiwan EPA. The input data of point and nonpoint sources are also estimated. The model simulates the concentrations of various forms of nutrients, CBOD and dissolved oxygen. A series of sensitivity runs was conducted to investigate the effects of point source loadings and river flow on the DO level in the river. It is demonstrated that the augmentation of river flow has as much effect on raising DO level as the reduction of point source loadings. The completion of the Taipei sewer project is expected to reduce the point source loadings by at least 75%. Under these reduced loadings, if the daily instream flow is maintained above the monthly Q75 flow throughout the year, the minimum DO concentration in the river would not fall below 1mg/L, which is the suffocation level for most fish species in the Hsintien Stream. (Q75 is the flow which is equaled or exceeded 75% of the days in the month.) The Feitsui Reservoir, Chingtan Dam and Chihtan Dam may impound water during the high flow periods and release freshwater to maintain the flow at the Q75 value in the Hsintien Stream during the drought periods.