Modeling of Vertical Solute Dispersion in a Sediment Bed Enhanced by Wave‐Induced Interstitial Flow1

Abstract: Mass (solute) transport in a stream or lake sediment bed has a significant effect on chemical mass balances and microbial activities in the sediment. A “1D vertical dispersion model” is a useful tool to analyze or model solute transfer between river or lake water and a sediment bed. Under a motionless water column, solute transfer into and within the sediment bed is by molecular diffusion. However, surface waves or bed forms create periodic pressure waves along the sediment/water interface, which in turn induce flows in the pores of the sediment bed. The enhancement of solute transport by these interstitial periodic flows in the pores has been incorporated in a 1D depth‐dependent “enhanced dispersion coefficient (D_E).” Typically, D_E diminishes exponentially with depth in the sediment bed. Relationships have been developed to estimate D_E as a function of the characteristics of sediment (particle size, hydraulic conductivity, and porosity) and pressure waves (wave length and height). In this paper, we outline and illustrate the calculation of D_E as well as the penetration depth (d_p) of the flow effect. Sample applications to illustrate the computational procedure are provided for dissolved oxygen transfer into a stream gravel bed and release of phosphorus from a lake bed. The sensitivity of the results to input parameter values is illustrated, and compared with the errors obtained when interstitial flow is ignored. Maximum values of D_E near the sediment surface can be on the order of 1 cm²/s in a stream gravel bed with standing waves, and 0.001 cm²/s in a fine sand lake bed under progressive surface waves, much larger than molecular diffusion coefficients.     

GEOSTATISTICAL ESTIMATION OF HORIZONTAL HYDRAULIC CONDUCTIVITY FOR THE KIRKWOOD‐COHANSEY AQUIFER1

ABSTRACT: The Kirkwood‐Cohansey aquifer has been identified as a critical source for meeting existing and expected water supply needs for southern New Jersey. Several contaminated sites exist in the region; their impact on the aquifer has to be evaluated using ground water flow and transport models. Ground water modeling depends on availability of measured hydrogeologic data (e.g., hydraulic conductivity, for parameterization of the modeling runs). However, field measurements of such critical data have inadequate spatial density, and their locations are often clustered. The goal of this study was to research, compile, and geocode existing data, then use geostatistics and advanced mapping methods to develop a map of horizontal hydraulic conductivity for the Kirkwood‐Cohansey aquifer. Spatial interpolation of horizontal hydraulic conductivity measurements was performed using the Bayesian Maximum Entropy (BME) Method implemented in the BMELib code library. This involved the integration of actual measurements with soft information on likely ranges of hydraulic conductivity at a given location to obtain estimate maps. The estimation error variance maps provide an insight into the uncertainty associated with the estimates, and indicate areas where more information on hydraulic conductivity is required.     

SPATIAL VARIABILITY OF CONDUCTWITY AS APPLIED TO DISTRIBUTED PARAMETER INFILTRATION MODELS1

ABSTRACT: Two infiltration models, called VVSIM (variable variance stochastic infiltration model) and EVVSIM (enhanced variable variance stochastic infiltration model), are developed in this study. A distributed parameter infiltration model can estimate the amount of infiltration over a field area by computing the infiltration over zones of the field area. Hydraulic conductivity is the most important parameter determining infiltration in simulations by infiltration models. The performance of an infiltration model depends on how well the model accommodates the complicated spatial distribution of hydraulic conductivity. The two proposed models include the effects of spatial correlation of the conductivity distribution. Virtual conductivity fields are generated using the turning bands method. Monte Carlo simulations show that the proposed models give infiltration estimates more accurate than those obtained by the other models employed.     

ESTIMATING STREAMFLOW AND ASSOCIATED HYDRAULIC GEOMETRY,THE MID‐ATLANTIC REGION,USA1

ABSTRACT: Methods to estimate streamflow and channel hydraulic geometry were developed for unpaged streams in the Mid‐Atlantic Region. Observed mean annual streamflow and associated hydraulic geometry data from 75 gaging stations in the Appalachian Plateau, the Ridge and Valley, and the Piedmont Physiographic Provinces of the Mid‐Atlantic Region were used to develop a set of power functions that relate streamflow to drainage area and hydraulic geometry to streamflow. For all three physiographic provinces, drainage area explained 95 to 98 percent of the variance in mean annual streamflow. Relationships between mean annual streamflow and water surface width and mean flow depth had coefficients of determination that ranged from R²= 0.55 to R²= 0.91, but the coefficient of determination between mean flow velocity and mean annual streamflow was lower (R²= 0.44 to R²= 0.54). The advantages of using the regional regression models to estimate streamflow over a conceptual model or a water balance model are its ease of application and reduced input data needs. The prediction of the regression equations were tested with data collected as part of the U.S. Environmental Protection Agency (USEPA) Environmental Monitoring and Assessment Program (EMAP). In addition, equations to transfer streamflow from gaged to ungaged streams are presented.     

Model Development and Calibration of a Saltwater Intrusion Model in Southern California1

Abstract: A nine‐layered confined‐unconfined flow and transport model is developed for the Alamitos saltwater intrusion barrier in Southern California. The conceptual model is based on the geological structure of the coastal aquifer system. The key parameters in the flow and transport models are calibrated using a two‐phase procedure which matches the types of data available for calibration. Because of the abundance of point measurements of hydraulic conductivity, the heterogeneous and random hydraulic conductivity field for each of the five aquifers is estimated by the geostatiscal method of natural‐neighbor‐kriging in Phase 1. In Phase 2, the longitudinal and transverse dispersivities in the transport model are estimated by a traditional inverse procedure that minimizes the least‐squares error for concentration (LSE‐CON). The minimum LSE‐CON is achieved near 15.2 and 1.52 m for the longitudinal and transverse dispersivities, respectively. Additional simulations with increasing transport parameter complexity did not yield significant improvements in LSE‐CON. Also, tracking least‐squares error for head while parametrically varying the transport parameters revealed there is a negligible interaction between predicted head and transport parameters.     

RESPONSE OF GREEN‐AMPT MEIN‐LARSEN SIMULATED RUNOFF VOLUMES TO TEMPORALLY AGGREGATED PRECIPITATION1

ABSTRACT: Data collection frequency in automated systems is user determined and can range from seconds to hours or days. Currently, there is no standard or recommended frequency interval for collecting precipitation data from automated systems for input to event‐based models such as Green‐Ampt Mein‐Larsen (GAML). Data from 47 storm events at seven locations were used to simulate the response of GAML excess rainfall to temporally aggregated precipitation data. No difference in model efficiency was recognized when comparing one‐minute interval data (R²= 1.00) to five‐minute data (R²= 1.00). Very little model efficiency was lost at a 10‐minute (R²= 0.96) interval. After 10‐minutes, decline in efficiency became more rapid with R²= 0.16 at one hour. The combined effect of time interval with respect to drainage area, hydraulic conductivity, maximum 30‐minute intensity, and total precipitation also revealed similar results.     

TRANSIENT STORAGE ASSESSMENTS OF DYE‐TRACER INJECTIONS IN RIVERS OF THE WILLAMETTE BASIN,OREGON1

ABSTRACT: Rhodamine WT dye‐tracer injections in rivers of the Willamette Basin yield concentration‐time curves with characteristically long recession times suggestive of active transient storage processes. The scale of drainage areas contributing to the stream reaches studied in the Willamette Basin ranges from 10 to 12,000 km². A transient storage assessment of the tracer studies has been completed using the U.S. Geological Survey's One‐dimensional Transport with Inflow and Storage (OTIS) model, which incorporates storage exchange and decay functions along with the traditional dispersion and advection transport equation. The analysis estimates solute transport of the dye. It identifies first‐order decay coefficients to be on the order of 10^?5/sec for the nonconservative Rhodamine W.T. On an individual subreach basis, the first‐order decay is slower (typically by an order of magnitude) than the transient storage process, indicating that nonconservative tracers may be used to evaluate transient storage in rivers. In the transient storage analysis, a dimensionless parameter (A_s/A) expresses the spatial extent of storage zone area relative to stream cross section. In certain reaches of Willamette Basin pool‐and‐riffle, gravel‐bed rivers, this parameter was as large as 0.5. A measure of the storage exchange flux was calculated for each stream subreach in the simulation analysis. This storage exchange is shown subjectively to be higher at higher stream discharges. Hyporheic linkage between streams and subsurface flows is the probable physical mechanism contributing to a significant part of this inferred active transient storage. Hyporheic linkages are further suggested by detailed measurements of river discharge with an Acoustic Doppler Current Profiler system delineating zones in two large rivers where water alternately enters and leaves the surface channels through gravel‐and‐cobble riverbeds. Measurements show patterns of hyporheic exchange that are highly variable in time and space.     

Applying Zonation Methods and Tabu Search to Improve the Ground‐Water Modeling1

Abstract: An optimization procedure combining zonation methods with the Tabu Search method is proposed to identify the spatial distribution of hydraulic conductivity field. Three zonation methods, Voronoi diagram (VD), multiplicatively weighted Voronoi diagram (MWVD), and pattern zonation are adopted for the parameterization purposes. With the consideration of the homogeneity and the heterogeneity, there are four spatial distributions of hydraulic conductivity designed to test whether the parameter structure can be successfully identified. The fitting residual error is first considered to determine an adequate number of zones without over parameterization. Then, the parameter uncertainty is evaluated the decision of the number of zones. The results indicate that the MWVD performs better than other two methods because the MWVD has better flexibility in describing the zonal boundaries with small number of decision variables.     

MODELING INFILTRATION WITH VARYING HYDRAULIC CONDUCTWITY UNDER SIMULATED RAINFALL CONDITIONS1

ABSTRACT: Physically-based models are extensively used to simulate the infiltration process under varied field conditions. Most models are based on the deterministic nature of input parameters related to the flow process (such as hydraulic conductivity). These models yield poor predictions of infiltration rates because they do not include the field-scale variations of flow parameters. The paper presents an approach for integrating the field-scale variability of hydraulic conductivity with an infiltration model to simulate infiltration under the rainfall conditions. A model describing the spatial structure of hydraulic conductivity has been developed using stochastic techniques. The stochastic structure of hydraulic conductivity was then incorporated in the Green-Ampt and Mein-Larson infiltration model. The model outputs on the instantaneous infiltration rates and cumulative infiltration were evaluated using the field infiltration data measured under simulated rainfall conditions. The results show that the combined model is capable of rep. resenting the instantaneous infiltration rates and cumulative infiltration of the study soils. The model may, therefore, be used to simulate the rainfall infiltration process for spatially-variable soils under the field conditions.     

Estimating Hydraulic Properties of Volcanic Aquifers Using Constant‐Rate and Variable‐Rate Aquifer Tests1

Abstract: In recent years the ground‐water demand of the population of the island of Maui, Hawaii, has significantly increased. To ensure prudent management of the ground‐water resources, an improved understanding of ground‐water flow systems is needed. At present, large‐scale estimations of aquifer properties are lacking for Maui. Seven analytical methods using constant‐rate and variable‐rate withdrawals for single wells provide an estimate of hydraulic conductivity and transmissivity for 103 wells in central Maui. Methods based on constant‐rate tests, although not widely used on Maui, offer reasonable estimates. Step‐drawdown tests, which are more abundantly used than other tests, provide similar estimates as constant‐rate tests. A numerical model validates the suitability of analytical solutions for step‐drawdown tests and additionally provides an estimate of storage parameters. The results show that hydraulic conductivity is log‐normally distributed and that for dike‐free volcanic rocks it ranges over several orders of magnitude from 1 to 2,500 m/d. The arithmetic mean, geometric mean, and median values of hydraulic conductivity are respectively 520, 280, and 370 m/d for basalt and 80, 50, and 30 m/d for sediment. A geostatistical approach using ordinary kriging yields a prediction of hydraulic conductivity on a larger scale. Overall, the results are in agreement with values published for other Hawaiian islands.     

A FINITE ELEMENT MODEL OF CONTAMINANT MOVEMENT IN GROUNDWATER1

ABSTRACT. Transient, two-dimensional solutions are developed which describe the movement and distribution of a conservative substance in a stream-aquifer system. The solutions are obtained by solving sequentially the groundwater flow and mass transport equations. A variational approach in conjunction with the finite element method is used to solve the groundwater flow equation. Galerkin's approach coupled with the finite element method is used to solve the mass transport equation. Linear approximated triangular elements and a centered scheme of numerical integration are employed to calculate the hydraulic head distribution and the concentration of solute in the flow region. The linear approximation used to define the concentration function within each element is not appropriate for cases involving steep concentration gradients. For such cases, higher order approximations are necessary to assure the continuity of gradients across interelemental boundaries. Numerical examples that illustrate the applicability of the model are presented.     

IDENTIFICATION OF AQUIFER DISPERSIVITIES: METHODS OF ANALYSIS AND PARAMETER UNCERTAINTY1

ABSTRACT: This paper presents a method for estimating aquifer dispersivities in solute transport models. Sensitivity equations are derived for the calculation of sensitivity coefficients. A modified Gauss-Newton algorithm is used to perform the least-squares minimization. A statistical procedure is outlined to assess reliability of the estimated parameters. The solute transport model is solved by the upstream weighted, multiple cell balance method which combines the concepts of local mass balance and finite element approximations. A one-dimensional solute transport problem in a vertical column system is first used to illustrate the inverse technique. A second example considers the parameter identification problem for three-dimensional solute transport with a unidirectional steady and uniform flow field. The third example solves the parameter identification problem in a three-dimensional, stream-aquifer, solute transport system with steady state flow. Numerical experiments are conducted to study data requirements for parameter identification.     

Contribution of In‐Channel Processes to Sediment Yield of an Urbanizing Watershed1

Abstract: A study was conducted between September 2003 and September 2006 to obtain baseline sediment inventories and monitor sediment transport and storage along a 3.7 km length of the channel of Valley Creek within Valley Forge National Historical Park, Pennsylvania. Valley Creek is a tributary of the Schuylkill River and drains an urbanizing 60.6 km² watershed that currently has 18% impervious land cover. Numerous field methods were employed to measure the suspended sediment yield, longitudinal profile, cross‐sections, banklines, and particle size distribution of the streambed. Suspended sediment yield for the watershed was measured at a USGS gage located just upstream of the park boundary between July 2004 and July 2005, the period corresponding to field surveys of bank erosion and channel change. The estimated suspended sediment yield of 95.7 t/km²/year is representative of a year with unusually high discharge, including a storm event that produced a peak of 78 m³/s, the second highest discharge on record for the USGS gage. Based on the median annual streamflow for the 24 years of record at the USGS gage from 1983 to 2006, the median annual sediment yield is estimated to be closer to 34 t/km²/year, considerably lower than median and mean values for other sites within the region. The mass of silt, clay, and fine sand derived from bank erosion along the 3.7 km study reach during the field survey period accounts for an estimated 2,340 t, equivalent to about 43% of the suspended sediment load. The mass of fine sediment stored in the bed along the study reach was estimated at 1,500 t, with about 330 t of net erosion during the study period. Although bank erosion appears to be a potentially dominant source of sediment by comparison with annual suspended sediment load, bed sediment storage and potential for remobilization is of the same order of magnitude as the mass of sediment derived from bank erosion.     

Nitrate Removal in a Restored Spring‐Fed Wetland,Pennsylvania, USA1

ABSTRACT:  In 2001, the 1.04‐ha Hornbaker wetland in south‐central Pennsylvania was restored by blocking an artificial drainage ditch to increase water storage and hydraulic retention time (HRT). A primary goal was to diminish downstream delivery of nitrate that enters the wetland from a limestone spring, its main source of inflow. Wetland inflow and outflow were monitored weekly for two years to assess nitrate flux, water temperature, pH, and specific conductivity. In Year 2, spring discharge was measured weekly to allow calculation of nitrate loads and hydraulic retention time. Surface runoff was confirmed to be a small fraction of wetland inflows via rainfall‐runoff modeling with TR‐55. The full dataset (n = 102) was screened to remove 13 weeks in which spring discharge constituted < 85% of total inflows because of high precipitation and surface runoff. Over two years (n = 89), mean nitrate‐nitrogen concentrations were 7.89 mg/l in inflow and 3.68 mg/l in outflow, with a mean nitrate removal of 4.19 mg/l. During Year 2 (n = 47), for which nitrate load data were available, the wetland removed an average of 2.32 kg N/day, 65% of the load. Nitrate removal was significantly correlated with HRT, water temperature, and the concentration of nitrate in inflow and was significantly greater during the growing season (5.36 mg/l, 64%) than during the non‐growing season (3.23 mg/l, 43%). This study indicates that hydrologic restoration of formerly drained wetlands can provide substantial water quality benefits and that the hydrologic characteristics of spring‐fed wetlands, in particular, support effective nitrogen removal.     

基于分形理论的皖江城市带城镇体系结构研究

运用分形理论对皖江城市带城镇体系的规模分布和城镇空间结构进行了具体分析,得出如下结论:在城镇规模分布方面,城市规模结构分布较为集中,中等位序的城市较多,首位城市作用不显著,人口分布低水平均衡;在城镇空间结构方面,城镇之间的空间关联程度、空间相互作用一般,空间布局较分散。在此基础上,提出了未来皖江城市带城镇体系发展的相应对策。     

Spatial Convergence in Major Dissolved Ion Concentrations and Implications of Headwater Mining for Downstream Water Quality

Spatial patterns in major dissolved solute concentrations were examined to better understand impact of surface coal mining in headwaters on downstream water chemistry. Sixty sites were sampled seasonally from 2012 to 2014 in an eastern Kentucky watershed. Watershed areas (WA) ranged from 1.6 to 400.5 km² and were mostly forested (58%–95%), but some drained as much as 31% surface mining. Measures of total dissolved solutes and most component ions were positively correlated with mining. Analytes showed strong convergent spatial patterns with high variability in headwaters (<15 km² WA) that stabilized downstream (WA > 75 km²), indicating hydrologic mixing primarily controls downstream values. Mean headwater solute concentrations were a good predictor of downstream values, with % differences ranging from 0.55% (Na⁺) to 28.78% (Mg²⁺). In a mined scenario where all headwaters had impacts, downstream solute concentrations roughly doubled. Alternatively, if mining impacts to headwaters were minimized, downstream solute concentrations better approximated the 300 μS/cm conductivity criterion deemed protective of aquatic life. Temporal variability also had convergent spatial patterns and mined streams were less variable due to unnaturally stable hydrology. The highly conserved nature of dissolved solutes from mining activities and lack of viable treatment options suggest forested, unmined watersheds would provide dilution that would be protective of downstream aquatic life.     

Hydraulics Near Unscreened Diversion Pipes in Open Channels: Large Flume Experiments

Most of the water diversions on the Sacramento and San Joaquin Rivers (California, United States) and their tributaries are currently unscreened. These unscreened diversions are commonly used for irrigation and are potentially harmful to migrating and resident fishes. A large flume (test section: 18.29 m long, 3.05 m wide and 3.20 m high) was used to investigate the hydraulic fields near an unscreened water diversion under ecologically and hydraulically relevant diversion rates and channel flow characteristics. We investigated all combinations of three diversion rates (0.28, 0.42, and 0.57 m³/s) and three sweeping velocities (0.15, 0.38, and 0.61 m/s), with one additional test at 0.71 m³/s and 0.15 m/s. We measured the three‐dimensional velocity field at seven cross sections near a diversion pipe and constructed regression equations of the observed maximum velocities near the pipe. Because the velocity components in three directions (longitudinal, transverse, and vertical) were significantly greater near the diversion pipe inlet compared with those farther from it, they cannot be neglected in the modeling and design of fish guidance and protection devices for diversion pipes. Our results should be of great value in quantifying the hydraulic fields that are formed around fish guidance devices to design more effective protection for fishes from entrainment into unscreened water‐diversion pipes.     

Empirical Estimation of Stream Discharge Using Channel Geometry in Low‐Gradient,Sand‐Bed Streams of the Southeastern Plains

Manning's equation is used widely to predict stream discharge (Q) from hydraulic variables when logistics constrain empirical measurements of in‐bank flow events. Uncertainty in Manning's roughness (n_M) is the major source of error in natural channels, and sand‐bed streams pose difficulties because flow resistance is affected by flow‐dependent bed configuration. Our study was designed to develop and validate models for estimating Q from channel geometry easily derived from cross‐sectional surveys and available GIS data. A database was compiled consisting of 484 Q measurements from 75 sand‐bed streams in Alabama, Georgia, South Carolina, North Carolina (Southeastern Plains), and Florida (Southern Coastal Plain), with six New Zealand streams included to develop statistical models to predict Q from hydraulic variables. Model error characteristics were estimated with leave‐one‐site‐out jackknifing. Independent data of 317 Q measurements from 55 Southeastern Plains streams indicated the model (Q = A_cR_H^0.6906S^0.1216; where A_c is the channel area, R_H is the hydraulic radius, and S is the bed slope) best predicted Q, based on Akaike's information criterion and root mean square error. Models also were developed from smaller Q range subsets to explore if subsets increased predictive ability, but error fit statistics suggested that these were not reasonable alternatives to the above equation. Thus, we recommend the above equation for predicting in‐bank Q of unbraided, sandy streams of the Southeastern Plains.     

Solute movement through an allophanic soil

Allophanic soils are widespread around the world, but little research has been done on their transport properties. This study reveals the effect of two soil water potential heads and two water-flow regimes of continuous and intermittent flow on solute transport through undisturbed soil columns of Horotiu silt loam (Typic Hapludand), an allophanic soil. Two different methods--breakthrough curves (BTCs) and time domain reflectometry (TDR)--were employed to determine the extent of preferential solute transport in the topsoil. The TDR data were also used to look at the depth dependence of the transport properties. The convection-dispersion equation (CDE) with the appropriate boundary conditions adequately described the movement of both Br and Cl under the various flow conditions. Although no preferential flow was found under the imposed unsaturated flow conditions, the flow of water and transport of solute became more uniform with depth. The results show that both Br and Cl are retarded in this allophanic soil. Retardation values range from 1.5 to 1.9, and, as the TDR data showed, increase from the depth of 5.0 to 10.0 cm. Intermittent leaching results showed that there was no effect on solute concentrations in the leachate following no-flow periods. This suggests that water and solute transport in this soil were either relatively uniform or that transverse mixing during flow was already fast enough to eliminate concentration gradients between regions of different "mobility."     

HYSTAR Sediment Model: Distributed Two‐Dimensional Simulation of Watershed Erosion and Sediment Transport Using Time‐Area Routing

An erosion and sediment transport component incorporated in the HYdrology Simulation using Time‐ARea method (HYSTAR) upland watershed model provides grid‐based prediction of erosion, transport and deposition of sediment in a dynamic, continuous, and fully distributed framework. The model represents the spatiotemporally varied flow in sediment transport simulation by coupling the time‐area routing method and sediment transport capacity approach within a grid‐based spatial data model. This avoids the common, and simplistic, approach of using the Universal Soil Loss Equation (USLE) to estimate erosion rates with a delivery ratio to relate gross soil erosion to sediment yield of a watershed, while enabling us to simulate two‐dimensional sediment transport processes without the complexity of numerical solution of the partial differential governing equations. In using the time‐area method for routing sediment, the model offers a novel alternative to watershed‐scale sediment transport simulation that provides detailed spatial representation. In predicting four‐year sediment hydrographs of a watershed in Virginia, the model provided good performance with R² of 0.82 and 0.78 and relative error of ?35% and 11% using the Yalin and Yang's sediment transport capacity equations, respectively. Prediction of spatiotemporal variation in sediment transport processes was evaluated using maps of sediment transport rates, concentrations, and erosion and deposition mass, which compare well with expected behavior of flow hydraulics and sediment transport processes.