Geochemical imaging of flow near an artificial recharge facility, Orange County, California

Critical for the management of artificial recharge operations is detailed knowledge of ground water dynamics near spreading areas. Geochemical tracer techniques including stable isotopes of water, tritium/helium-3 (T/3He) dating, and deliberate gas tracer experiments are ideally suited for these investigations. These tracers were used to evaluate flow near an artificial recharge site in northern Orange County, California, where approximately 2.5 x 10(8) m3 (200,000 acre-feet) of water are recharged annually. T/3He ages show that most of the relatively shallow ground water within 3 km of the recharge facilities have apparent ages < 2 years; further downgradient apparent ages increase, reaching > 20 years at approximately 6 km. Gas tracer experiments using sulfur hexafluoride and xenon isotopes were conducted from the Santa Ana River and two spreading basins. These tracers were followed in the ground water for more than two years, allowing subsurface flow patterns and flow times to be quantified. Results demonstrate that mean horizontal ground water velocities range from < 1 to > 4 km/year. The leading edges of the tracer patch moved at velocities about twice as fast as the center of mass. Leading edge velocities are important when considering the potential transport of microbes and other "time sensitive" contaminants and cannot be determined easily with other methods. T/3He apparent ages and tracer travel times agreed within the analytical uncertainty at 16 of 19 narrow screened monitoring wells. By combining these techniques, ground water flow was imaged with time scales on the order of weeks to decades.     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

EFFECT OF SHOCK WAVES ON RILL FORMATION

Hydraulic mechanism of fill formation was studied theoretically and experimentally. It was assumed that the impact of varied boundary on overland flow results in fluctuating of water surface, and shock waves that may contribute to the formation of fills. Both theoretical derivation and laboratory experiments were used to compare the hydraulic characteristics of flows with and without shock waves. Results showed that shock waves can lead to an increase in flow depth, flow velocity, and turbulence intensity. Consequently, flow shear stress or stream energy increase dramatically and fill headcuts may occur where shock waves converge.     

Capabilities of Large-scale Particle Image Velocimetry to characterize shallow free-surface flows

Irrespective of their spatial extent, free-surface shallow flows are challenging measurement environments for most instruments due to the relatively small depths and velocities typically associated with these flows. A promising candidate for enabling measurements in such conditions is Large-scale Particle Image Velocimetry (LSPIV). This technique uses a non-intrusive approach to measure two-dimensional surface velocity fields with high spatial and temporal resolutions. Although there are many publications documenting the successful use of LSPIV in various laboratory and field open-channel flow situations, its performance has not been equally substantiated for measurement in shallow flows. This paper aims at filling in this gap by demonstrating the capabilities of LSPIV to: (a) accurately evaluate complex flow patterns in shallow channel flows; and (b) estimate depth in shallow flows using exclusively LSPIV measurements. The demonstration is provided by LSPIV measurements in three shallow flow laboratory situations with flow depths ranging from 0.05 to 0.31 m. The obtained measurements illustrate the LSPIV flexibility and reliability in measuring velocities in shallow and low-velocity (near-zero) flows. Moreover, the technique is capable to evaluate and map velocity-derived quantities that are difficult to document with alternative measurement techniques (e.g. vorticity and shear stress distributions and mapping of large-scale structure in the body of water).     

Modeling changes in rill erodibility and critical shear stress on native surface roads

This study investigated the effect of cumulative overland flow on rill erodibility and critical shear stress on native surface roads in central Idaho. Rill erodibility decreased exponentially with increasing cumulative overland flow depth; however, critical shear stress did not change. The study demonstrated that road erodibility on the studied road changes over the course of one or more consecutive overland flow events. Therefore, model simulations that fail to take into consideration this change will probably over-estimate sediment yields. An exponential function describing the relationship between rill erodibility and cumulative overland flow depth is presented as a basis for future model development for simulating erosion on native surface roads. Copyright © 2008 John Wiley & Sons, Ltd.     

Size characteristics of sediment in interrill overland flow on a semiarid hillslope,Southern Arizona

This study examines the size characteristics of sediment removed from a semiarid hillslope by interrill overland flow. Rainfall simulation experiments were conducted on a runoff plot 18 m wide and 35 m long established on a piedmont hillslope in southern Arizona. The top of the plot coincided with the hillslope divide, and its outlet was located within a shallow rill. Samples of runoff were obtained from two cross-sections located in the interrill portion of the plot upslope of the rill and from a calibrated flume through which was directed interrill overland flow reaching the bottom of the plot. Analyses of sediment contained in these samples showed that sediment in interrill flow is finer than the matrix soil. The fineness of the interrill sediment compared to the matrix soil appears to be due to the inability of interrill overland flow to transport the coarser fraction of the sediment supplied to it by raindrop detachment. This finding implies that the rate of soil erosion in interrill areas is not. as is commonly supposed, limited by the rate at which raindrops can detach sediment but by the rate at which they detach sediment of a size that the overland flow is competent to transport. The relative fineness of sediment eroded from this hillslope is consistent with other evidence for the recent evolution of shrub-covered hillslopes in southern Arizona.     

The European Soil Erosion Model (EUROSEM): a dynamic approach for predicting sediment transport from fields and small catchments

The European Soil Erosion Model (EUROSEM) is a dynamic distributed model, able to simulate sediment transport, erosion and deposition over the land surface by rill and interill processes in single storms for both individual fields and small catchments. Model output includes total runoff, total soil loss, the storm hydrograph and storm sediment graph. Compared with other erosion models, EUROSEM has explicit simulation of interill and rill flow; plant cover effects on interception and rainfall energy; rock fragment (stoniness) effects on infiltration, flow velocity and splash erosion; and changes in the shape and size of rill channels as a result of erosion and deposition. The transport capacity of runoff is modelled using relationships based on over 500 experimental observations of shallow surface flows. EUROSEM can be applied to smooth slope planes without rills, rilled surfaces and surfaces with furrows. Examples are given of model output and of the unique capabilities of dynamic erosion modelling in general. © 1998 John Wiley & Sons, Ltd.     

Seasonal vegetation and management influence overland flow velocity and roughness in upland grasslands

There is considerable interest in how headwater management may influence downstream flood peaks in temperate humid regions. However, there is a dearth of data on flow velocities across headwater hillslopes and limited understanding of whether surface flow velocity is influenced by seasonal changes in roughness through vegetation cycles or management. A portable hillslope flume was used to investigate overland flow velocities for four common headwater grassland habitats in northern England: Low-density Grazing, Hay Meadow, Rank Grassland and Juncus effusus Rush pasture. Overland flow velocity was measured in replicate plots for each habitat, in response to three applied flow rates, with the experiments repeated during five different periods of the annual grassland cycle. Mean annual overland flow velocity was significantly lower for the Rank Grassland habitat (0.026 m/s) followed by Low-density Grazing and Rushes (0.032 and 0.029 m/s), then Hay Meadows (0.041 m/s), which had the greatest mean annual velocity (examples from 12 L/min flow rate). Applying our mean overland flow velocities to a theoretical 100 m hillslope suggests overland flow is delayed by >1 hr on Rank Grassland when compared to Hay Meadows in an 18 mm storm. Thus grassland management is important for slowing overland flow and delaying peak flows across upland headwaters. Surface roughness was also strongly controlled by annual cycles of vegetation growth, decay, grazing and cutting. Winter overland flow velocities were significantly higher than in summer, varying between 0.004 m/s (Rushes, November) and 0.034 m/s (Rushes, June); and velocities significantly increased after cutting varying between 0.006 m/s (Hay meadows, July) and 0.054 m/s (Hay meadows, September). These results show that seasonal vegetation change should be incorporated into flood modelling, as cycles of surface roughness in grasslands strongly modify overland flow, potentially having a large impact on downstream flood peak and timing. Our data also showed that Darcy-Weisbach roughness approximations greatly over-estimated measured flow velocities.     

DYNAMIC AND KINEMATIC WAVE MODELS FOR 1-D RILL FLOW: A COMPARISON

1INTRODUCTIONInprocess-basedoverlandflowandsoilerosionmodels,surfacerunoffonahillslopeisoftenrepresentedaseitherbroadsheetfloworflowinrillswithassumedrectangularchannelcrosssections(e.g.,Bairdelal1992,NSERL1995).Inmostcasesthehydraulicsofoverlandflowiscalculatedbyusingthekinematicwavemodel,whichisasirnplificationofthedynamicwavemodel(theequationsystemofSaint-Venantequationandequationofcontinuity).Forabroadsheetoverlandflowonhillslopeduetorainfallexcess,theequationofcontinuityiswherehis…     

INITIATION OF SEDIMENT MOTION IN LAMINAR OVERLAND FLOW

IINTRODUCTIONWhileriverflowsareusuallydeepandturbulent,overlandflowisextremelyshallowandcanbelaminar,transitionalandturbulent.Becauseoftheshallownessoftheflolw,overlandflowhydraulicsisgreatlyaffectedbysurfaceroughness,raindropimpact,andinthecaseoflaminarflow,flui(Iviscosity.Theinitiationofsedimentmovementinoverlandflowisthereforeexpectedtodifferfromthatinriverflows.InriverstUdies,bedshearStressgbhastraditionallybeenusedtocharacterizethecriticalflowconditionatwhichsedimentbeginstomove.At…     

Interstitial flow through preferential flow paths in the hyporheic zone of the Oberer Seebach,Austria

We investigated interstitial flow velocities in the Oberer Seebach, Austria, with NaCl tracer injections at a sediment depth of 30 cm to estimate the hydraulic conditions experienced by invertebrates inhabiting the hyporheic zone. Flow velocity measured with tracers is taken as travel time of the water along a straight line between injection and sampling points, although the water flows around sediment particles, and thus travels a somewhat longer distance. From sections of stream sediment in which the interstitial spaces were replaced by concrete, we estimated that this difference amounts, on average, to 27% and used this factor to correct the results of our velocity measurements. Corrected interstitial water velocities ranged from 0.01 to 1.32 cm s^-1 and were independent of surface discharge. We also studied spatial flow patterns in the bed sediments with long-term tracer injections. The three-dimensional distribution of tracer concentrations 24 hours after the start of the injection indicated that interstitial water preferentially flows in a complex network of areas of high hydraulic connectivity. Reynolds numbers for flow in the hyporheic pore space ranged from 0.1 to 489, implying that the flow environment varies from laminar up to the zone of transition to turbulent flow. Therefore, invertebrates may have a size-related active choice of areas where either friction drag or pressure drag predominates. The consequence of flow patterns, such as those observed in our study, is that small-scale variability of hydraulic conditions may be an important determinant of the patchy invertebrate distribution in bed sediments.     

Tracing variability of run‐off generation in mountainous permafrost of semi‐arid north‐eastern Mongolia

The headwaters of mountainous, discontinuous permafrost regions in north‐eastern Mongolia are important water resources for the semi‐arid country, but little is known about hydrological processes there. Run‐off generation on south‐facing slopes, which are devoid of permafrost, has so far been neglected and is totally unknown for areas that have been affected by recent forest fires. To fill this knowledge gap, the present study applied artificial tracers on a steppe‐vegetated south‐facing and on two north‐facing slopes, burned and unburned. Combined sprinkling and dye tracer experiments were used to visualize processes of infiltration and water fluxes in the unsaturated zone. On the unburned north‐facing slope, rapid and widespread infiltration through a wet organic layer was observed down to the permafrost. On the burned profile, rapid infiltration occurred through a combusted organic and underlying mineral layer. Stained water seeped out at the bottom of both profiles suggesting a general tendency to subsurface stormflow (SSF). Ongoing SSF could directly be studied 24 h after a high‐intensity rainfall event on a 55‐m hillslope section in the burned forest. Measurements of water temperature proved the role of the permafrost layer as a base horizon for SSF. Repeated tracer injections allowed direct insights into SSF dynamics: A first injection suggested rather slow dispersive subsurface flow paths; whereas 18 h later, a second injection traced a more preferential flow system with 20 times quicker flow velocities. We speculate that these pronounced SSF dynamics are limited to burned slopes where a thermally insulating organic layer is absent. On three south‐facing soil profiles, the applied tracer remained in the uppermost 5 cm of a silt‐rich mineral soil horizon. No signs of preferential infiltration could be found, which suggested reduced biological activity under a harsh, dry and cold climate. Instead, direct observations, distributed tracers and charcoal samples provided evidence for the occurrence of overland flow. Copyright © 2014 John Wiley & Sons, Ltd.     

Determining friction coefficients for interrill flows: the significance of flow filaments and backwater effects

Friction coefficients in overland flows are customarily estimated from mean flow properties (depth, velocity, slope) that subsume spatial variations in flow arising from two major causes: microtopography and obstacles. This paper uses laboratory experiments in shallow flumes to examine the extent of non‐uniformity in flow conditions associated with each cause. Randomly placed emergent obstacles in a flume with a shallow axial channel generally yielded higher hydraulic roughness than the same pattern of obstacles on a planar flume, as well as greater variation in roughness as the obstacle locations were altered. In both flumes, hydraulic roughness fell with increasing Reynolds number for 10% obstacle cover, showed a flattening trend at 20% cover, and exhibited a convex‐downward trend at 30% obstacle cover. These results indicate the progressive onset of flow controls at narrow gaps in the obstacle field. In such flows, the use of mean flow properties conceals the existence of two main subdivisions of flow: flow filaments and backwater flows. In the experiments, flow filaments involved velocities more than twice the overall mean, whereas backwater flows were much slower than the mean. The existence of fast‐moving flow filaments may be significant in understanding soil transport in surface runoff, and backwater depths may modify splash detachment. Similarly, friction coefficients that fail to reflect these important non‐uniform flow components may not be optimal for hydraulic calculations or in erosion models. It is concluded that new approaches to observing and processing flow data may be required, in order to avoid the loss of important flow detail that is entailed in assuming uniform flow conditions. Copyright © 2003 John Wiley & Sons, Ltd.     

Unsteady sedimentation in nonuniform rills

Runoff and sediment were measured from agricultural land exposed to controlled simulated rainfall. We extended the kinematic unsteady overland sedimentation theory of prismatic channels to the experiments by considering both hydraulics and sediment dynamics of rill flow for changing flow geometries of nonuniform and unsteady rill development. The characteristic unimodal concentration peak observed in the experiments and the changing channel geometry were interpreted in theoretical terms. Overland sedimentation in unsteady nonprismatic rills under uniform rainfall can be described with kinematic models of flow, entrainment and deposition applied to developing flow geometries; this is not possible with sheet-flow models. Other interpretations are considered and experimental needs are identified.     

Stochastic overland flows Part 1: Physics-based evolutionary probability distributions

A theoretical solution framework to the nonlinear stochastic partial differential equations (SPDE) of the kinematic wave and diffusion wave models of overland flows under stochastic inflows/outflows, stochastic surface roughness field and stochastic state of flows was obtained. This development was realized by means of an eigenfunction representation of the time-space overland flow depths, and by transforming the problem into the phase space. By using Van Kampen's lemma and the cumulant expansion theory of Kubo-Van Kampen-Fox, the deterministic partial differential equation (PDE) for the evolutionary probability density function (pdf) of overland flow depths was finally obtained. Once this deterministic PDE is solved for the time-varying pdf of overland flow depths, then the time-space varying pdf of overland flow depths can be obtained by a transformation given in the text. In this solution framework it is possible to incorporate the stochastic dynamic behavior of the parameters and of the forcing functions of the overland flow process. For example, not only the individual rainfall duration and fluctuating rain intensity characteristics but also the sequential behavior of rainfall patterns is incorporated into the evolutionary probability density function of overland flow depths.     

Rill length and plot‐scale effects on the hydrogeomorphologic response of gravelly roadbeds

Although unpaved roads are well‐recognized as important sources of Hortonian overland flow and sediment in forested areas, their role in agriculturally‐active rural settings still lacks adequate documentation. In this study, we assessed the effect of micro‐catchment size, slope, and ground cover on runoff and sediment generation from graveled roadbeds servicing a rural area in southern Brazil. Fifteen replications based on 30‐min‐long simulated rainfall experiments were performed at constant rainfall intensities of 22–58 mm h^?1 on roadbeds with varying characteristics including ~3–7 m² micro‐catchment areas, 2–11° slopes, 2–9.7‐m‐long shallow rill features, and 30–100% gravel cover. The contributions of micro‐catchment size and rill length were the most important physical characteristics affecting runoff response and sediment production; both the size of the micro‐catchment and the length of the rills were inversely related to sediment loss and this contradicts most of the rill erosion literature. The effect of micro‐catchment size on runoff and sediment response suggests a potentially problematic spatial‐scale subjectivity of experimental plot results. The inverse relationship between rill length and sediment generation is interpreted here as related to the predominance of coarse fragments within rills, the inability of the shallow flows generated during the simulations to erode this sediment, and their role as zones of net sediment storage. Copyright © 2015 John Wiley & Sons, Ltd.     

A velocity projection framework for inferring shallow water currents from surface tracer fields

To the extent that sea surface temperature and colors can be considered passive tracers, their motions can be tracked to estimate the current velocities, or a conservation equation can be invoked to relate their temporal variations to the velocities. We investigate the latter, the so-called tracer inversion problem, with a particular focus on (1) the conditions under which the problem can be rendered over-determined for least squares solutions, (2) the possibility of using the tracer conservation equation within the “velocity projection” framework to estimate subsurface current profiles in shallow coastal waters, and (3) the accuracy of the tracer inversion calculation in terms of the data resolution and noise. The velocity projection framework refers to relating surface motion, either measured directly or made visible by tracers, to the subsurface current motion through the equations of motion. The accuracy of the tracer inversion calculation is quantified in terms of the spatial and temporal resolution of the tracer distribution. In the presence of irreducible tracer noise, the accuracy of the inversion rapidly degrades, and it is shown that the inversion with velocity projection can help improve accuracy. The tracer inversion method developed in this study is applied to the satellite sea surface temperature data, and the velocity result is compared to the velocity measurements made with the shore-based HF Coastal Current Radar. The potential of improving the velocity estimation with the present approach is indicated.     

Effect of saltating sediment on flow resistance and bed roughness in overland flow

The acceleration of saltating grains by overland flow causes momentum to be transferred from the flow to the grains, thereby increasing flow resistance and bed roughness. To assess the impact of saltating sediment on overland flow hydraulics, velocity profiles in transitional and turbulent flows on a fixed sand-covered bed were measured using hot-film anemometry. Five discharges were studied. At each discharge, three flows were measured: one free of sediment, one with a relatively low sediment load, and one with a relatively high sediment load. In these flows from 83 to 90 per cent of the sediment was travelling by saltation. As a result, in the sediment-laden flows the near-bed velocities were smaller and the velocity profiles steeper than those in the equivalent sediment-free flows. Sediment loads ranged up to 87·0 per cent of transport capacity and accounted for as much as 20·8 per cent of flow resistance (measured by the friction factor) and 89·7 per cent of bed roughness (measured by the ratio of the roughness length to median grain diameter). It is concluded that saltating sediment has a considerable impact on overland flow hydraulics, at least on fixed granular beds. Saltation is likely to have a relatively smaller effect on overland flow on natural hillslopes and agricultural fields where form and wave resistance dominate. Still, saltation is generally of greater significance in overland flow than in river flow, and for this reason its effect on overland flow hydraulics is deserving of further study. © 1998 John Wiley & Sons, Ltd.     

Line-source multi-tracer test for assessing high groundwater velocity

Segmented line-source multi-tracer injection is suggested as an effective method for assessing groundwater velocities and flow directions in subsurfaces characterized by high water flux. Modifying the common techniques of injecting a tracer into a well became necessary after point-source natural and forced gradient tracer tests ended with no reliable information on the local groundwater flow. The tracer's line-source increases the likelihood of success of the test and could provide additional information regarding the lateral heterogeneity of the aquifer. In a field experiment conducted in the northwestern part on the Dead Sea coast, tracers were injected into an 8-m-long line injection system perpendicular to the assumed flow direction. The injection system was divided into four separate segments with four different tracers. An array of five boreholes located within a 10 × 10 m area downstream was used for monitoring the tracers' transport. Two dye tracers (uranine and Na naphthionate) were injected in a long pulse of several hours into two of the injection pipe segments. Two other tracers (Rhenium oxide and Gd-DTPA) were instantaneously injected into the other two segments. The tracers were detected 0.7 to 2.3 h after injection in four of the five observation wells, located 2.3 to 10 m away from the injection system. The groundwater velocity was determined to be ～80 to 170 m/d, based on the recoveries of the tracers. The groundwater flow direction was derived based on the arrival of the tracers and was found to be quite consistent with the apparent direction of the hydraulic gradient.