TRANSPORT CHARACTERISTICS OF PHOSPHORUS IN CHANNELIZED AND MEANDERING STREAMS1

ABSTRACT: Comparisons were made between rates of movement of orthophosphate in a canal and a meandering stream. The meander system had greater algal and macrophyte phosphate uptake rates, and lower plankton and sediment release rates compared to the canal. Chemical precipitation and direct rainfall influences on orthophosphate movement were insignificant relative to other terms. The major source of phosphorus to both systems was from upland runoff. The impact of this source was greater on the meandering system due to the smaller channel volume. When secondary effects of meandering were considered such as marsh inundation, the net orthophosphate movement within the meandering channel was less than that for the canal; due to the lower concentrations of phosphorus in marsh effluent waters. Field experiments were conducted to compare the longitudinal dispersion coefficient between a canal and meandering river system; the meandering stream had a dispersion coefficient over 17 times that measured for the canal. Rates of orthophosphate movement were combined into a single mass transport equation, and a numerical solution was obtained. Internal river and canal channel processes were overshadowed by external point source loadings.     

LINE SOURCE DISPERSION WITH APPLICATION TO MIXING IN RIVER CHANNELS1

ABSTRACT: Two-dimensional solutions for transient dispersion of nonconservative dispersants in uniform flow resulting from a transverse line source of variable concentration are obtained using multiple integral transformations. In general, the solutions are in integral forms, which can be efficiently evaluated using Filon's quadratures. Examples are presented for cases of practical interest. Applicability of the solution for modeling dispersion in natural river channel where the distribution of flows across the channel are nonuniform is discussed.     

Chaotic Behavior and Pollution Dispersion Characteristics in Engineered Tidal Embayments: A Numerical Investigation1

Abstract: This paper presents a numerical investigation of approaches to enhance the mixing and dispersion processes in tidal areas by effecting changes in the natural estuary system. It compares the impact of various estuary modifications stemming from human intervention to pollutant dispersion and chaotic flow within the estuary including the implications of alteration of the original channel shape, change of the channel bathymetry, and modification of the tidal signal. Our findings indicate that chaotic flow analysis is similar in many regards, but not all, to conventional dispersion analysis. Specifically, we conclude that (1) simplification of the flow regime reduces chaotic flow patterns and tracer particle dispersion, (2) creation of extensively protruding barriers and/or installation of barriers on opposite sides of the main stem of the estuary enhances particle dispersion and chaotic mixing, (3) installation of underwater berms has relatively minor beneficial, but highly localized, effects on chaotic regime and particle dispersion, and (4) increasing the tidal signal amplitude was shown to increase chaotic and dispersion properties of the estuarine system. A parametric study investigating the effect of several geometrical configurations and tidal signals on characteristics of chaotic flows concludes the paper.     

Chromatographic partitioning of impurities (H2S) contained in a CO2 stream injected into a deep saline aquifer: Part 2. Effects of flow conditions

Laboratory experiments reported in a companion paper were carried out to examine the chromatographic partitioning of impurities contained in a stream of CO₂ injected into a deep saline aquifer. The solubility of the impurity gas in the CO₂ stream compared to that of CO₂, the in situ conditions of pressure, temperature and water salinity, and the concentration of the impurity gas affect the partitioning of the two gases. For CO₂ streams containing H₂S, numerical simulations reported here have successfully replicated the laboratory results including the breakthrough of CO₂ ahead of H₂S. Sensitivity analysis performed with the numerical model has shown that flow conditions, controlled by such parameters as medium permeability, pressure gradient, dispersion, gas mobility and flow direction, affect the breakthrough time and separation of the two gases, leading to delayed or earlier breakthrough, and increasing or decreasing the time lag between the breakthrough of the two gases. Vertical bottom-up flow leads to earlier breakthrough, while top-down flow leads to delayed breakthrough. These results are important in establishing monitoring strategies at CO₂ storage sites and in evaluating the risks associated with the possible leakage of injected CO₂ that contains impurities.     

MORPHOLOGIC RESPONSE OF SUBALPINE STREAMS TO TRANSBASIN FLOW DIVERSION1

ABSTRACT: Flow has been diverted from Rocky Mountain streams for many years with little documentation of the impacts on physical form and biological function of the losing stream system. This study addresses whether differences in form can be detected in subalpine step-pool, plane bed, or pool-riffle channels and linked to changes in flow regime from diversion. Total annual discharge was reduced between 20 and 60 percent and average annual peak flow was reduced up to 45 percent in the subalpine systems assessed; channels were diverted between 20 and 100 years. Expected impacts include reduced conveyance and changes in vegetation growth patterns because formerly active surfaces are colonized by riparian species, effectuating shrinking channel capacity. In this study, reduced channel width is used as an indicator of morphologic response. Observed reductions in width, ranging from 35 to 50 percent at some sites, resulted not only from vegetation invasion of stable surfaces but also from the development of an inset beside former cut banks. This observation, however, was restricted to wider pool-riffle channels with gravel bars. Outside of these areas, morphologic changes were either subtle or absent. The absence of widespread response is attributed primarily to periodic “flooding” of the diverted channel. During wetter-than-average years when excess water is available, minimal flow is diverted and the hydrograph resembles a free-flowing regime. The release of high flow to the natural channel potentially offsets changes in form incurred during “dry” periods. The stable nature and structure of subalpine channels also contributes to the absence of reduced capacity.     

MODELING THE LONG TERM IMPACTS OF USING RIGID STRUCTURES IN STREAM CHANNEL RESTORATION1

Abstract: Natural channel designs often incorporate rigid instream structures to protect channel banks, provide grade control, promote flow deflection, or otherwise improve channel stability. The long term impact of rigid structures on natural stream processes is relatively unknown. The objective of this study was to use long term alluvial channel modeling to evaluate the effect of rigid structures on channel processes and assess current and future stream channel stability. The study was conducted on Oliver Run, a small stream in Pennsylvania relocated due to highway construction. Field data were collected for one year along the 107 m reach to characterize the stream and provide model input, calibration, and verification data. FLUVIAL-12 was used to evaluate the long term impacts of rigid structures on natural channel adjustment, overall channel stability, and changing form and processes. Based on a consideration of model limitations and results, it was concluded that the presence of rigid structures reduced channel width-to-depth ratios, minimized bed elevation changes due to long term aggradation and degradation, limited lateral channel migration, and increased the mean bed material particle size throughout the reach. Results also showed how alluvial channel modeling can be used to improve the stream restoration design effort.     

EFFECT OF FLOW VELOCITY ON SEDIMENT OXYGEN DEMAND: COMPARISON OF THEORY AND EXPERIMENTS1

ABSTRACT: Theoretical equations that establish the relationship between sediment oxygen demand (SOD) in a lake and the flow velocity and dissolved oxygen concentration in the bulk water already exist. These theoretical equations for oxygen consumption in the sediment express biological consumption with Michaelis-Menten kinetics, and chemical consumption by a first order reaction. Data from laboratory experiments that were conducted to validate the theoretical equations also exist. These experiments were performed in a laboratory channel with well defined flow characteristics for three types of sediments. Herein, the theoretical equations are used to model the experimental data for the three types of sediments. The values used for the parameters in the theoretical equations are determined by iteration until a best fit is obtained for the relationship of SOD to flow velocity from both the theoretical model and experimental data. The goodness of fit is measured by the standard error of prediction and the regression coefficient.     

ANALYSIS OF ROCK STRESS AROUND A LIQUID WASTE INJECTION SYSTEM IN ILLINOIS1

: The danger to the environment associated with the injection of liquid industrial wastes into a deep, confined, subsurface rock formation may arise from the transport of the waste laterally or vertically in the formation. The pattern of lateral transport, which can take place as a result of convection as well as dispersion and diffusion, can be determined by an approximate analytical solution to the mass transport equation. Vertical transport may take place through both natural fractures and fractures created by hydrostatic stresses generated around the well during injection. To determine the stresses, we used the finite element method to get a numerical solution of the flow equation. We applied a solution of the flow equation to calculate the stress buildup and decay for the Jones & Laughlin Steel Corporation's injection well near Hennepin in Putnam County, Illinois. According to our computations, the stress buildup due to injection is about 0.16 pounds per square inch per foot - psi - (0.362 Newton per square centimeter per meter), which, added to normal pressure, makes an estimated total stress of 0.60 psi/ft (1.36 Newton/cm²/m). That pressure is insufficient to cause fracturing of the Cambrian Eau Claire aquitard, the confining bed for the disposal zone.     

VELOCITY AND CONCENTRATION DISTRIBUTIONS OF SEDIMENT‐LADEN OPEN CHANNEL FLOW1

ABSTRACT: A mathematical model to predict both velocity and concentration distributions for sediment‐laden open channel flow is developed. Velocity profiles are derived by theoretical analysis and numerical method. Logarithmic law and semi‐empirical wake function concept are not adopted. An empirical equation for the ratio of sediment exchange and fluid diffusion coefficients is considered to solve the diffusion equation for suspended‐sediment concentration profiles. Four sets of experimental data from previous researchers are compared to numerical calculation. In the engineering applications, velocity and concentration profiles of sediment‐laden flow can be predicted simultaneously by the present model with the measured velocity and sediment‐concentration at reference level.     

Modeling Hydrologic Benefits of Low Impact Development: A Distributed Hydrologic Model of The Woodlands,Texas

Low Impact Development (LID) is alternative design approach to land development that conserves and utilizes natural resources to minimize the potential negative environmental impacts of development, such as flooding. The Woodlands near Houston, Texas is one of the premier master‐planned communities in the United States. Unlike in a typical urban development where riparian corridors are often replaced with concrete channels, pervious surfaces, vegetation, and natural drainage pathways were preserved as much as possible during development. In addition, a number of detention ponds were strategically located to manage runoff on site. This article uses a unique distributed hydrologic model, Vflo?, combined with historical (1974) and recent (2008 and 2009) rainfall events to evaluate the long‐term effectiveness of The Woodlands natural drainage design as a stormwater management technique. This study analyzed the influence of LID in The Woodlands by comparing the hydrologic response of the watershed under undeveloped, developed, and highly urbanized conditions. The results show that The Woodlands drainage design successfully reflects predeveloped hydrologic conditions and produces peak flows two to three times lower than highly urbanized development. Furthermore, results indicate that the LID practices employed in The Woodlands successfully attenuate the peak flow from a 100‐year design event, resulting in flows comparable to undeveloped hydrologic conditions.     

一个研究街道峡谷流场及浓度场特征的三维数值模式

目前，研究街道峡谷内流场及机动车排放污染物的扩散行为特征所采用的主要方法为：采用野外测试法和物理模拟法，而采用三维数值模拟方法研究此问题的工作很少。本文创建了一个研究微尺度街道峡谷内流场及机动车排放污染物扩散特征的三维数值模式，即首次采用伪不定常方法，利用K——E闭合方案，建立了一个模拟城市街道峡谷内流场及污染物扩散特征与街道峡谷风场、街道几何结构及两侧建筑物高度对称性之间的复杂关系的三维数值模式。经过与实际监测资料及风洞实验对比，结果表明此三维模式具有较好的模拟精度，能够很好模拟峡谷内的风场及街谷几何结构对街道峡谷内流场及浓度场特征的影响，有很强的实用性。     

BEDROCK CONTROLS ON STREAM CHANNEL ENLARGEMENT WITH URBANIZATION,NORTH CENTRAL TEXAS1

Loss due to channel erosion in the Dallas, Texas, area is estimated to approach one-half million dollars in the last several years. Hydrogeomorphic analysis of natural and urban chalk and shale watersheds was performed in the central Texas area on watersheds ranging in size from 0.5 to 10 square miles in an effort to more adequately predict channel enlargement due to urbanization. Chalk watersheds were found to have greater drainage density, greater channel slope, lower sinuosity, and greater discharge per unit area than similar sized shale watersheds under natural conditions. With subsequent urbanization of the watersheds, chalk channel enlargement was from 12 to 67 percent greater than shale channel enlargement for similar sized watersheds. Greater enlargement in chalk channels is attributed to greater channel velocities and unit tractive force. Vegetation seems to play a significant role in influencing channel adjustments to the new flow regimes brought on by urbanization. Channel response to urbanization is documented and specific nonstructural guidelines are proposed which could reduce structural loss along urban stream channels.     

Sorption, mobility, and transformation of estrogenic hormones in natural soil

Potent estrogenic hormones are consistently detected in the environment at low concentration, yet these chemicals are strongly sorbed to soil and are labile. The objective of this research was to improve the understanding of the processes of sorption, mobility, and transformation for estrogens in natural soils, and their interaction. Equilibrium and kinetic batch sorption experiments, and a long-term column study were used to study the fate and transport of 17beta-estradiol and its primary metabolite, estrone, in natural soil. Kinetic and equilibrium batch experiments were done using radiolabeled 17beta-estradiol and estrone. At the concentrations used, it appeared that equilibrium sorption for both estrogens was achieved between 5 and 24 h, and that the equilibrium sorption isotherms were linear. The log K(oc) values for 17beta-estradiol (2.94) and estrone (2.99) were consistent with previously reported values. Additionally, it was found that there was rate-limited sorption for both 17beta-estradiol (0.178 h(-1)) and estrone (0.210 h(-1)). An approximately 42 h long, steady-flow, saturated column experiment was used to study the transport of radiolabeled 17beta-estradiol, which was applied in a 5.00 mg L(-1) solution pulse for 44 pore volumes. 17beta-estradiol and estrone were the predominant compounds detected in the effluent. The effluent breakthrough curves were asymmetric and the transport modeling indicated that sorption was rate-limited. Sorption rates and distributions of the estrogens were in agreement between column and batch experiments. This research can provide a better link between the laboratory results and observations in the natural environment.     

Impact of Temporal Resolution of Flow‐Duration Curve on Sediment Load Estimation1

Chen, Li, Rina Schumer, Anna Knust, and William Forsee, 2011. Impact of Temporal Resolution of Flow‐Duration Curve on Sediment Load Estimation. Journal of the American Water Resources Association (JAWRA) 48(1): 145‐155. DOI: 10.1111/j.1752‐1688.2011.00602.x Abstract: Estimates of a channel’s annual sediment transport capacity typically incorporate annualized flow‐duration curves. Average daily flow data, commonly used to develop flow‐duration curves, may not adequately describe sediment‐transporting flows in arid and semiarid ephemeral streams. In this study, we examined impacts of varied temporal resolution flow data on annual sediment load estimation. We derived flow‐duration curves for eight sites in the Southwestern United States based on both 15‐min and daily‐averaged flow data. We then estimated sediment loads for both flow‐duration curves using the Sediment Impact Analysis Method, implemented in HEC‐RAS. When average daily flow is used to generate flow‐duration curves, sediment load estimation is lower by up to an order of magnitude. This trend is generally unaffected by uncertainty associated with sediment particle size or hydraulic roughness. The ratio of sediment loads estimated by 15‐min versus daily‐averaged flow‐duration curves is strongly correlated with channel slope, being greater on steep‐slope channels. Sediment loads estimated by the two types of flow‐duration curves are closely correlated, suggesting possible relationships for improving predictions when high‐temporal resolution data are unavailable. Results also suggest that the largest flow contributes significantly to total sediment load, and thus will greatly impact ephemeral stream geomorphology in arid and semiarid regions.     

EFFECTS OF TIME STEP SIZE IN IMPLICIT DYNAMIC ROUTING1

ABSTRACT The effects of the size of the Δt time step used in the integration of the implicit difference equations of unsteady open-channel flow are determined for numerous typical hydrographs with durations in the order of days or even weeks. Truncation errors related to the size of the Δt time step cause a numerical distortion (dispersion and attenuation) of the computed transient. The magnitude of the distortion is related directly to the size of the time step, the length of channel reach, and the channel resistance and inversely to the time of rise of the hydrograph. The type of finite difference expression which replaces spatial derivatives and non-derivative terms in the partial differential equations of unsteady flow has an important influence on the magnitude of the numerical distortion, as well as the numerical stability of the implicit difference equations. Time step sizes in the range of 3 to 6 hrs generally tend to minimize the combination of required computation time and numerical distortion of transients having a time of rise of the order of several days.     

Transport of sulfadiazine in undisturbed soil columns: effects of flow rate, input concentration and pulse duration

Antibiotics reach soils via spreading of manure or sewage sludge. Knowledge on the transport behavior of antibiotics in soils is needed to assess their environmental fate. The effect of flow rate and applied mass, i.e., input concentration and pulse duration, on the transport of 14C-sulfadiazine (SDZ; 4-aminoN-pyrimidin-2-yl-benzenesulfonamide) was investigated with soil column experiments and numerical studies. Sulfadiazine was applied in pulses (6.8, 68 or 306 h) under steady-state (0.051 and 0.21 cm h(-1)) and intermittent flow conditions and at two input concentrations (0.57 and 5.7 mg L(-1)). Breakthrough curves (BTCs) of 14C were measured and for one experiment concentrations of SDZ, and its transformation products 4-(2-iminopyrimidin-1(2H)-yl)aniline (An-SDZ) and N(1)-2-(4-hydroxypyrimidinyl)benzenesulfanilamide (4-OH-SDZ) were determined. After finalizing the leaching experiments, 14C was quantified in different slices of the columns. A lower flow rate led to remarkably lower eluted masses compared with the higher flow rates. All BTCs could be described well using a three-site attachment-detachment model for which a common set of parameters was determined. However, the BTC obtained with the high input concentration was slightly better described with a two-site isotherm-based model. The prediction of the concentration profiles was good with both model concepts. The fitted sorption capacities decreased in the order SDZ > 4-OH-SDZ > An-SDZ. Overall, the experiments reveal the presence of similar mechanisms characterizing SDZ transport. The dependence of model performance on concentration implies that although the three-site attachment-detachment model is appropriate to predict the transport of SDZ in soil columns, not all relevant processes are adequately captured.     

EMERGENCY WATER SUPPLIES FROM GROUND WATER IN HUMID REGIONS1

ABSTRACT: This paper focuses on the development and testing of a mathematical model of an emergency ground water supply operated principally during periods of low streamflow. The process of ground water withdrawal and recharge is simulated taking account of streamflow, water demand, evapotranspiration, natural and artificial recharge and increased evapotranspiration due to artificial recharge, ground water pumpage, and streamflow contribution to pumped water. The model determines whether natural recharge is possible in less time than the return period of drought and also whether artificial recharge is needed. By simulating operation over a long period of time, the model can examine different droughts of short and long duration and can test the operating rules for ground water storage development in an area. Submodels analyze the components of the operating process including ground water flow into the stream, seepage losses, stream portion of well discharge due to induced infiltration and recharge from rainfall or water spreading. The model has been tested for areas in the humid northeastern United States.     

RIVER DYE GAGING1

Procedures have been developed (a) to inject a tracer at a constant rate below the water surface at selected points across a stream and (b) to deal with suspended sediment. Mixing remained far from complete in relatively long channels, owing to channel and flow divergence with uncertainty where to sample downstream and which marginal sample values to include for flow calculation. These problems are encountered when mixing is largely dependent on transverse diffusion. Accurate and replicable results were obtained where dye was injected upstream and detected downstream from riffles that induced thorough turbulent mixing. Dye gaging should be practical in gorges or wherever flow is turbulent across the whole width of a channel.