Observing lingering hyporheic storage using electrical resistivity: variations around stream restoration structures,Crabby Creek,PA

Time‐lapse geophysical surveys can map lingering hyporheic storage by detecting changes in response to saline tracer. Tracer tests were conducted in Crabby Creek, an urban stream outside Philadelphia, to examine the influence of stream restoration structures and variable sediment thickness. We compared electrical resistivity surveys with extensive well sampling (57 wells) in two 13.5‐m‐long reaches, each with a step drop created by a J‐hook. The two step drops varied in tracer behaviour, based on both the well data and the geophysical data. The well data showed more variation in arrival time where the streambed sediment was thick and was more uniform where sediment was thin. The resistivity in the reach with thin sediment showed lingering tracer in the hyporheic zone both upstream and downstream from the J‐hook. In the second reach where the sediment was thicker, the lingering tracer in the hyporheic zone was more extensive downstream from the J‐hook. The contrasting results between the two reaches from both methods suggested that sediments influenced hyporheic exchange more than the step at this location. Resistivity inversion differed from well data in both reaches in that it showed evidence for tracer after well samples had returned to background, mapping lingering tracer either upstream or downstream of a step. We conclude that resistivity surveys may become an important tool for hyporheic zone characterization because they provide information on the extent of slow moving fluids in the hyporheic zone, which have the potential to enhance chemical reactions. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of theory and experiment for solute transport in weakly heterogeneous bimodal porous media

In this work, the influence of non-equilibrium effects on solute transport in a weakly heterogeneous medium is discussed. Three macro-scale models (upscaled via the volume averaging technique) are investigated: (i) the two-equation non-equilibrium model, (ii) the one-equation asymptotic model and (iii) the one-equation local equilibrium model. The relevance of each of these models to the experimental system conditions (duration of the pulse injection, dispersivity values…) is analyzed. The numerical results predicted by these macroscale models are compared directly with the experimental data (breakthrough curves). Our results suggest that the preasymptotic zone (for which a non-Fickian model is required) increases as the solute input pulse time decreases. Beyond this limit, the asymptotic regime is recovered. A comparison with the results issued from the stochastic theory for this regime is performed. Results predicted by both approaches (volume averaging method and stochastic analysis) are found to be consistent.     

Laboratory and numerical investigation of flow and transport near a seepage-face boundary

Laboratory and numerical modeling investigations were completed to study the unconfined ground water flow and transport processes near a seepage-face boundary. The laboratory observations were made in a radial sand tank and included measurements of the height of the seepage face, flow velocity near the seepage face, travel time distribution of multiple tracer slugs, and streamlines. All the observations were reliably reproduced with a three-dimensional, axi-symmetric, variably saturated ground water flow model. Physical data presented in this work demonstrate and quantify the importance of three-dimensional transport patterns within a seepage-face zone. The results imply that vertically averaged flow models that employ Dupuit approximations might introduce error in the analysis of localized solute transport near a seepage-face boundary. The experimental dataset reported in this work will also be of interest for those who are attempting to validate a numerical algorithm for solving ground water and contaminant discharge patterns near a surface-water boundary.     

A numerical modelling and experimental study of flow width dynamics on alluvial fans

Alluvial fans are dynamic landforms, the evolution of which is controlled by both external environmental forcing (climate, tectonics and base level change) and internal process‐form feedbacks. The latter include changes in flow configuration (between sheetflow and channelized flow states), driven by aggradation and degradation, which may in turn promote changes in sediment transport capacity. Recent numerical modelling indicates that such feedbacks may lead to dramatic and persistent fan entrenchment in the absence of external forcing. However, the parameterization of flow width within such models is untested to date and is subject to considerable uncertainty. This paper presents results from an experimental study of flow width dynamics on an aggrading fan in which spatial and temporal patterns of fan inundation are monitored continuously using analysis of digital vertical photography. Observed flow widths are compared with results from a simple theoretical model developed for non‐equilibrium (aggradational) conditions. Results demonstrate that the theoretical model is capable of capturing the first‐order characteristics of width adjustment over the course of the experiment, and indicate that flow width is a function of fan aggradation rate. This illustrates that models of alluvial flow width derived for equilibrium conditions may have limited utility in non‐equilibrium situations, despite their widespread use to date. Copyright © 2009 John Wiley & Sons, Ltd.     

Transport of three herbicides in ground water at Twin Lake test site,Chalk River,Ontario, Canada

A field tracer test performed under natural flow conditions at the Twin Lake test site, Chalk River Laboratories of the Atomic Energy of Canada Ltd. in Chalk River, Ontario, Canada, using tritium and three herbicides (Chlortoluron, Terbuthylazine, and Pendimethalin) was interpreted using the dispersion equation with a combined reaction model. The reaction model couples an instantaneous equilibrium reaction governed by a linear adsorption isotherm with a reversible or irreversible kinetic reaction of the first order, and decay. An improved interpretation method consists of a simultaneous fitting of theoretical concentration and mass-recovery curves to the experimental data, which leads to a more reliable determining of reaction models and improves the accuracy of fitting. Tritium served as the reference tracer to determine the flow velocity, dispersivity, and the recovery of the herbicides. Chlortoluron was slightly delayed by equilibrium exchange with strongly reduced concentration due to an irreversible kinetic reaction and/or decay. Terbuthilazine was slightly delayed by equilibrium exchange, with strongly reduced concentration due to a reversible kinetic reaction with some influence of decay. A strong equilibrium reaction and a strong reversible kinetic reaction without degradation governed the transport of Pendimethalin, reducing considerably its concentration. The results obtained show that simulations based only on Kd and decay constant, especially if these parameters are found in the laboratory, may considerably differ from those performed with reaction parameters determined in properly performed field tests. The dominant reaction types, and the values of parameters found in the study, supply useful information on the transport of the investigated herbicides in sandy aquifers under natural flow conditions.     

Variability of unsaturated Bromide fluxes as measured through a layered volcanic vadose zone in New Zealand

The measured drainage fluxes through a layered volcanic vadose zone exhibited high spatial variability as a consequence of heterogeneous flow conditions. The drainage flux variability was quantified using automated equilibrium tension lysimeters, installed in close‐proximity and resulted in high variability in the Br masses recovered from a conservative tracer experiment. The primary cause of the heterogeneous flow was attributed to textural changes occurring at the interface between volcanic layers, resulting in development of funnel‐flow patterns, and further enhanced by the existence of hydrophobic conditions. The Br recoveries in individual automated equilibrium tension lysimeters were used to determine the corresponding variable sizes of the surface areas contributing drainage to the lysimeters. The tracer experiment confirmed the existence of unsaturated lateral transport occurring at the sloping interface of the coarse Taupo Ignimbrite material with the silty Palaeosol layer at approximately 4.2 m depth. This study demonstrates that measurements of both flux and solute concentrations at multiple locations are essential when heterogeneous flow is suspected to be present, to be able to determine reliable estimates of contaminant leaching through the vadose zone at the plot scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of moving storms on attainment of equilibrium discharge

Although rainfall is assumed spatially uniform in conventional hydrological modelling for rainfall–runoff simulations, moving storms have been shown to have substantial influence on flow hydrographs. In this study, criteria for attainment of the equilibrium discharge from watersheds subjected to moving storms were examined. Non-linear numerical kinematic-wave models were developed to simulate runoff from an overland plane and from a V-shaped catchment. Dimensional analysis was applied to obtain the independent variables to be used as control factors in performing a series of numerical tests. The results indicate that, for storms moving downstream, runoff can attain equilibrium discharge even though the storm length is shorter than the watershed length and the rainfall duration is less than the time to equilibrium of the watershed for stationary uniform storms. The phenomenon of attainment of equilibrium discharge from watersheds subjected to moving storms is contradictory to conventional hydrologic design, which assumes the storm duration must equal the time to equilibrium to attain the maximum discharge. Copyright © 2007 John Wiley & Sons, Ltd.     

Imaging hyporheic zone solute transport using electrical resistivity

Traditional characterization of hyporheic processes relies upon modelling observed in‐stream and subsurface breakthrough curves to estimate hyporheic zone size and infer exchange rates. Solute data integrate upstream behaviour and lack spatial coverage, limiting our ability to accurately quantify spatially heterogeneous exchange dynamics. Here, we demonstrate the application of near‐surface electrical resistivity imaging (ERI) methods, coupled with experiments using an electrically conductive stream tracer (dissolved NaCl), to provide in situ imaging of spatial and temporal dynamics of hyporheic exchange. Tracer‐labelled water in the stream enters the hyporheic zone, reducing electrical resistivity in the subsurface (to which subsurface ERI is sensitive). Comparison of background measurements with those recording tracer presence provides distributed characterization of hyporheic area (in this application, ∼0·5 m²). Results demonstrate the first application of ERI for two‐dimensional imaging of stream‐aquifer exchange and hyporheic extent. Future application of this technique will greatly enhance our ability to quantify processes controlling solute transport and fate in hyporheic zones, and provide data necessary to inform more complete numerical models. Copyright © 2010 John Wiley & Sons, Ltd.     

A STUDY ON THE EQUILIBRIUM PROFILE FOR THE LUOSHAN-HANKOU REACH IN THE MIDDLE YANGTZE RIVER

1 INTRODUCTION Large-scale flood disasters have frequently occurred in the middle Yangtze River since the 1990抯. The Jingjiang River and Dongting Lake (Fig. 1) comprise the most serious area of flood disasters. The main characteristic of recent disasters is low discharge and high water stage. Recent research has begun to pay more attention to the important role of sediment deposition (Li and Ni, 1998). Though the Yangtze River is not an overloaded river, the amount of sediment trans…     

Density contrast effects on tracer dispersion in variable aperture fractures

We investigate the development of preferential flow paths and anomalous dispersion resulting from weak density contrasts in the course of tracer experiments in a tortuous natural fracture. The processes are first documented by the non-invasive measurement of the fracture aperture and of the time-resolved distribution of the tracer using Positron Emission Projection Imaging. Then, numerical simulations of the three-dimensional tracer transport in the fracture are performed to explore the parameters that control the development and the persistence of the tracer localization, as a function of the density contrast between the tracer and the resident solution. Results reveal that density contrasts representative of what could be expected in borehole and laboratory tracer tests can induce irreversible localization along preferential channels. As density contrast increases, the correlation between velocity and aperture distributions vanishes whereas (i) velocity field increasingly correlates with the fracture median plan elevation and (ii) the longitudinal dispersion coefficient increases. The anomalous velocity distribution may persist well after the injection stops due to the occurrence of tracer trapped zone.     

Drainage in finite-sized unsaturated zones

After the initiation of gravity drainage, water is often assumed to be either (a) draining under unit gradient, or (b) at capillary/gravity equilibrium. Both of these simplifications can be useful, but the regimes of validity of each assumption must be delineated. Water pressures are measured versus time and distance as water drains out of a 1.6 m long sand column to determine the relative effects of capillary and gravitational forces during drainage. For medium sized sands (0.15–0.3 mm in diameter), the capillary pressure is constant in space in a large region of the column for over 12 days, and the water continues to flow under unit gradient for relatively long time scales. Similar results are seen for finer sands, but with a much faster approach to equilibrium. Numerical simulations and analytical estimates are presented and compare favorably to the measurements. Together, the experimental, theoretical and analytical results are used to calculate when capillary/gravity equilibrium is reached as a function of porous media properties and length of the unsaturated zone. The ratio of the length of the unsaturated zone to the bubbling pressure is a key parameter in determining the drainage regime, and that even for relatively short unsaturated zones the equilibrium time scale can be on the order of years.     

Progress in numerical modeling of subducting plate dynamics

The core concerns of plate tectonics theory are the dynamics of subducting plates, which can be studied by integrating multidisciplinary fields such as seismology, mineral physics, rock geochemistry, geological formation studies, sedimentology, and numerical simulations. By establishing a theoretical model and solving it with numerical methods, one can replicate the dynamic effects of a subducting plate, quantifying its evolution and the surface response. Simulations can also explain the observations and experimental results of other disciplines. Therefore, numerical models are among the most important tools for studying the dynamics of subducting plates. This paper provides a review on recent advances in the numerical modeling of subducting plate dynamics. It covers various aspects, namely, the origin of plate tectonics, the initiation process and thermal structure of subducting slab, and the main subduction slab dynamics in the upper mantle, mantle transition zone, and lower mantle. The results of numerical models are based on the theoretical equations of mass, momentum, and energy conservation. To better understand the dynamic progress of subducting plates, the simulation results must be verified in comparisons with the results from natural observations by geology, geophysics and geochemistry. With the substantial increase in computing power and continuous improvement of simulation methods, numerical models will become a more accurate and efficient means of studying the frontier issues of Earth sciences, including subducting plate dynamics.     

The effect of transient storage on nitrate uptake lengths in streams: an inter‐site comparison

Surface water–groundwater interaction in the hyporheic zone may enhance biogeochemical cycling in streams, and it has been hypothesized that streams exchanging more water with the hyporheic zone should have more rapid nitrate utilization. We used simultaneous conservative solute and nitrate addition tracer tests to measure transient storage (which includes hyporheic exchange and in‐stream storage) and the rate of nitrate uptake along three reaches within the Red Canyon Creek watershed, Wyoming. We calibrated a one‐dimensional transport model, incorporating transient storage (OTIS‐P), to the conservative solute breakthrough curves and used the results to determine the degree of transient storage in each reach. The nitrate uptake length was quantified from the exponential decrease in nitrate concentration with distance during the tracer tests. Nitrate uptake along the most downstream reach of Red Canyon Creek was rapid (turnover time K^?1_c = 32 min), compared with nitrate uptake reported in other studies (K^?1_c = 12 to 551 min), but other sites within the watershed showed little nitrate retention or loss. The uptake length S_w‐NO^?₃ for the most downstream reach was 500 m and the mass transfer coefficient V_f‐NO^?₃ was 6·3 m min^?1. Results from 15 other nitrate‐addition tracer tests were used to create a regression model relating transient storage and measures of stream flow to nitrate uptake length. The model, which includes specific discharge and transient storage area, explains almost half the variability in nitrate uptake length (adjusted R² = 0·44) and is most effective for comparing sites with very different stream characteristics. Although large differences in specific discharge and storage zone area explain inter‐site differences in nitrate uptake, other unmeasured variables, such as available organic carbon and microbial community composition, are likely important for predicting differences in nitrate uptake between sites with similar specific discharge rates and storage zone areas, such as when making intra‐site comparisons. Copyright © 2007 John Wiley & Sons, Ltd.     

Quantitative analysis of numerically induced mixing in a coastal model application

In recent years, various attempts have been made to estimate the amount of numerical mixing in numerical ocean models due to discretisation errors of advection schemes. In this study, a high-resolution coastal model using the ocean circulationmodel GETM is applied to the Western Baltic Sea, which is characterised by energetic and episodic inflows of dense bottom waters originating from the Kattegat. The model is equipped with an easy-to-implement diagnostic method for obtaining the numerical mixing which has recently been suggested. In this diagnostic method, the physical mixing is defined as the mean tracer variance decay rate due to turbulent mixing. The numerical mixing due to discretisation errors of tracer advection schemes is defined as the decay rate between the advected square of the tracer variance and the square of the advected tracer, which can be directly compared to the physical variance decay. The source and location of numerical mixing is further investigated by comparing different advection schemes and analysing the amount of numerical mixing in each spatial dimension during the advection time step. The results show that, for the setup used, the numerically and physically induced mixing have the same orders of magnitude but with different vertical and horizontal distributions. As the main mechanism for high numerical mixing, vertical advection of tracers with strong vertical gradients has been identified. The main reason for high numerical mixing is due to bottom-following coordinates when density gradients, especially for regions of steep slopes, are advected normal to isobaths. With the bottom-following coordinates used here, the horizontal gradients are reproduced by a spurious sawtooth-type profile where strong advection through, but not along, the vertical coordinate levels occurs. Additionally, the well known relation between strong tracer gradients and high velocities on the one and high numerical mixing on the other hand is approved quantitatively within this work.     

STUDY ON THE MATHEMATICAL MODEL OF HEAVY－METAL POLLUTANT TRANSPOPT－TRANSFORMATION IN RIVERS

ＳＴＵＤＹＯＮＴＨＥＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＯＦＨＥＡＶＹ－ＭＥＴＡＬＰＯＬＬＵＴＡＮＴＴＲＡＮＳＰＯＰＴ－ＴＲＡＮＳＦＯＲＭＡＴＩＯＮＩＮＲＩＶＥＲＳ￥Ｓｕｉ－ｌｉａｎｇＨｕａｎｇａｎｄＺｈａｏ－ｈｕｉＷａｎ（ＰｏｓｔｄｏｃｔｏｒａｌＦｅ...     

Geometric characteristics and evolution of a tidal channel network in experimental setting

This paper reports on a laboratory study that aims to reproduce a tidal channel network, in order to enhance the understanding of the morphodynamic evolution of the channel characteristics as the network expands and finally reaches equilibrium. A high‐resolution laser system scanned the bed topography at different time steps creating multiple digital elevation models of the channel network. Two hundred and seventy individual channel segments are analyzed and cross‐correlated in terms of their width, depth and length. The laboratory results show positive linear correlations between depth and width as well as between length and width of channel segments of the network configuration at final equilibrium. In a downstream direction, channels appear to widen more than they deepen, indirectly a sign that discharge has a stronger control on channel width than on depth. In contrast to fluvial drainage networks that commonly display fractal and scale‐invariant behavior, the geometric properties of the experimental tidal creek network shows scale dependence. Channel attributes exhibit consistent patterns of exponentially decreasing abundance, with increasing creek length, depth and width. The nature of the observed exponential distributions within creek attributes (width, depth, length) allows for statistical predictability of relative creek attribute dimensions downstream and through time. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluating nitrate uptake in a Rocky Mountain stream using labelled 15N and ambient nitrate chemistry

Background aqueous chemistry and ¹⁵N_nitrate tracer injection methods were used to calculate in‐stream nitrate uptake metrics at Red Canyon Creek, a third‐order stream in the Rocky Mountains in the state of Wyoming, United States. ‘Net’ nitrate uptake lengths, which reflect both nitrate uptake and regeneration, and ‘gross’ nitrate uptake lengths, which exclude re‐mineralization, were quantified separately from background nitrate chemistry and ¹⁵N labelling tracer data, respectively. Gross nitrate uptake lengths, from tracer injections of ¹⁵N labelled nitrate, ranged from 502 to 3140 m. Net nitrate uptake lengths, from background nitrate chemistry downstream of a point source, ranged from 1170 to 4330 m. Diurnal changes in uptake lengths suggest the importance of nitrate utilization by autotrophs in the stream and benthic zone. The differences between net and gross nitrate uptake lengths along lower reaches of Red Canyon Creek allowed us to estimate the nitrate regeneration rate, which was 0·056–0·080 µmol m^?2 s^?1 during the day and 0·0062–0·0083 µmol m^?2 s^?1 at night. Spatial patterns of streambed pore water chemistry indicate those areas of the hyporheic zone where denitrification was likely occurring. Permanent log dams generated stronger redox gradients in the hyporheic zone than areas with transient beaver dams. By combining isotopically labelled nitrate additions, estimates of uptake from background aqueous nitrate chemistry and characterization of redox conditions in the hyporheic zone, we were able to determine the nitrate regeneration rate and the redox processes responsible for nitrogen cycling in the hyporheic zone. Copyright © 2010 John Wiley & Sons, Ltd.     

Influence of Outcrop Scale Fractures on the Effective Stiffness of Fault Damage Zone Rocks

We combine detailed mapping and microstructural analyses of small fault zones in granodiorite with numerical mechanical models to estimate the effect of mesoscopic (outcrop-scale) damage zone fractures on the effective stiffness of the fault zone rocks. The Bear Creek fault zones were active at depths between 4 and 15 km and localize mesoscopic off-fault damage into tabular zones between two subparallel boundary faults, producing a fracture-induced material contrast across the boundary faults with softer rocks between the boundary faults and intact granodiorite outside the boundary faults. Using digitized fault zone fracture maps as the modeled fault geometries, we conduct nonlinear uniaxial compression tests using a novel finite-element method code as the experimental “laboratory” apparatus. Map measurements show that the fault zones have high nondimensional facture densities (>1), and damage zone fractures anastamose and intersect, making existing analytical effective medium models inadequate for estimation of the effective elastic properties. Numerical experiments show that the damage zone is strongly anisotropic and the bulk response of the fault zone is strain-weakening. Normal strains as small as 2% can induce a reduction of the overall stiffness of up to 75%. Fracture-induced effective stiffness changes are large enough to locally be greater than intact modulus changes across the fault due to juxtaposition of rocks of different lithologies; therefore mesoscopic fracturing is as important as rock type when considering material or bimaterial effects on earthquake mechanics. These results have important implications for earthquake rupture mechanics models, because mesoscopic damage zone fractures can cause a material contrast across the faults as large as any lithology-based material contrast at seismogenic depths, and the effective moduli can be highly variable during a single rupture event.     

Linkages between bedload displacements and topographic change

Changes in bed topography that build and maintain channel morphology are driven by the displacements of individual particles, either though their entrainment or deposition. However, the linkages between these topographic changes and individual grain displacements have not been comprehensively addressed, as many historical tracer studies have not included coincident topographic data. In this study, we compare the movements of bedload tracers to the differences in repeat topographic surveys across four gravel-bed river reaches. To do this, we apply a 1-D Bayesian survival process model to the starting and ending locations of tracers. This model estimates downstream trapping probabilities, which represent the likelihood that a given segment of channel will “trap” an entrained particle. We then adapt this model to estimate downstream trapping probabilities using digital elevation models of difference and compare the results. The estimates from the tracer and topographic trapping models showed general alignment, meaning that tracers were preferentially trapped in segments that experienced deposition along the channel. Thus, tracers in this study were able to identify downstream differences in bedload transport. The comparison also highlighted that tracer-estimated trapping probabilities were larger than topographically estimated ones. This supports previous observations that sediment travel distances estimated using tracers are shorter than those estimated using morphological methods. We find that the differences between these two estimates vary systematically across study environments. These variations are attributable to either study design (i.e., tracers being larger than the median size of the sediment that deforms the bed) or differences in compensating scour and fill. We explore potential causes for differences in compensating scour and fill, including hydrograph shape, sediment delivery regime, channel deformation style, and channel width, highlighting that morphodynamics needs to be considered in designing bedload tracer studies.