A Straightforward Random Walk Model for Fast Push‐Pull Tracer Test Evaluation

In this article, we present a straightforward random walk model for fast evaluation of push‐pull tracer tests. By developing an adaptive algorithm, we overcome the problem of manually defining how many particles have to be used to simulate the transport problem. Beside this, we validate the random walk model by evaluating a push‐pull tracer test with drift phase and confirm the results with MT3DMS. The random walk model took less than 1% of computational time of MT3DMS, thus allowing a remarkable faster evaluation of push‐pull tracer tests.     

Improved Method for Estimating Reaction Rates During Push-Pull Tests

The breakthrough curve obtained from a single-well push-pull test can be adjusted to account for dilution of the injection fluid in the aquifer fluid. The dilution-adjusted breakthrough curve can be analyzed to estimate the reaction rate of a solute. The conventional dilution-adjusted method assumes that the ratios of the concentrations of the nonreactive and reactive solutes in the injection fluid vs. the aquifer fluid are equal. If this assumption is invalid, the conventional method will generate inaccurate breakthrough curves and may lead to erroneous conclusions regarding the reactivity of a solute. In this study, a new method that generates a dilution-adjusted breakthrough curve was theoretically developed to account for any possible combination of nonreactive and reactive solute concentrations in the injection and aquifer fluids. The newly developed method was applied to a field-based data set and was shown to generate more accurate dilution-adjusted breakthrough curves. The improved dilution-adjusted method presented here is simple, makes no assumptions regarding the concentrations of the nonreactive and reactive solutes in the injection and aquifer fluids, and easily allows for estimating reaction rates during push-pull tests.     

Evaluating the Use of In‐Well Heat Tracer Tests to Measure Borehole Flow Rates

Recent research has demonstrated the use of in‐well heat tracer tests monitored by a fiber optic distributed temperature sensing (DTS) system to characterize borehole flow conditions in open bedrock boreholes. However, the accuracy of borehole flow rates determined from in‐well heat tracer tests has not been evaluated. The purpose of the research presented here is to determine whether borehole flow rates obtained using DTS‐monitored in‐well heat tracer tests are reasonable, and to evaluate the range of flow rates measureable with this method. To accomplish this, borehole flow rates measured using in‐well heat tracer tests are compared to borehole flow rates measured in the same boreholes using an impeller or heat pulse flowmeter. A comparison of flow rates measured using in‐well heat tracer tests to flow rates measured with an impeller flowmeter under the same conditions showed good agreement. A comparison of in‐well heat tracer test flow rate measurements to previously‐collected heat pulse flowmeter measurements indicates that the heat tracer test results produced borehole flow rates and flow profiles similar to those measured with the heat pulse flowmeter. The results of this study indicate that borehole flow rates determined from DTS‐monitored in‐well heat tracer tests are reasonable estimates of actual borehole flow rates. In addition, the range of borehole flow rates measurable by in‐well heat tracer tests spans from less than 10^?1 m/min to approximately 10¹ m/min, overlapping the ranges typically measurable with an impeller flowmeter or a heat pulse flowmeter, making in‐well heat tracer testing a versatile borehole flow logging tool.     

Estimating Groundwater Mounding in Sloping Aquifers for Managed Aquifer Recharge

Design of managed aquifer recharge (MAR) for augmentation of groundwater resources often lacks detailed data, and simple diagnostic tools for evaluation of the water table in a broad range of parameters are needed. In many large‐scale MAR projects, the effect of a regional aquifer base dip cannot be ignored due to the scale of recharge sources (e.g., wadis, streams, reservoirs). However, Hantush's (1967) solution for a horizontal aquifer base is commonly used. To address sloping aquifers, a new closed‐form analytical solution for water table mound accounts for the geometry and orientation of recharge sources at the land surface with respect to the aquifer base dip. The solution, based on the Dupiuit‐Forchheimer approximation, Green's function method, and coordinate transformations is convenient for computing. This solution reveals important MAR traits in variance with Hantush's solution: mounding is limited in time and space; elevation of the mound is strongly affected by the dip angle; and the peak of the mound moves over time. These findings have important practical implications for assessment of various MAR scenarios, including waterlogging potential and determining proper rates of recharge. Computations are illustrated for several characteristic MAR settings.     

Quasi-Saturated Layer: Implications for Estimating Recharge and Groundwater Modeling

This study presents an extension of the concept of “quasi-saturation” to a quasi-saturated layer, defined as the uppermost dynamic portion of the saturated zone subject to water table fluctuations. Entrapped air here may cause substantial reductions in the hydraulic conductivity (K) and fillable pore water. Air entrapment is caused by a rising water table, usually as a result of groundwater recharge. The most significant effects of entrapped air are recharge overestimation based on methods that use specific yield (S_y), such as the water table fluctuation method (WTF), and reductions in K values. These effects impact estimation of fluid flow velocities and contaminant migration rates in groundwater. In order to quantify actual groundwater recharge rates and the effects of entrapped air, numerical simulations with the FEFLOW (Version 7.0) groundwater flow model were carried out using a quasi-saturated layer for a pilot area in Rio Claro, Brazil. The calculated recharge rate represented 16% of the average precipitation over an 8-year period, approximately half of estimates using the WTF method. Air entrapment amounted to a fillable porosity of 0.07, significant lower that the value of 0.17 obtained experimentally for S_y. Numerical results showed that the entrapped air volume in the quasi-saturated layer can be very significant (0.58 of the air fraction) and hence can significantly affect estimates of groundwater recharge and groundwater flow rates near the water table.     

Comparison of Phytoscreening and Direct‐Push‐Based Site Investigation at a Rural Megasite Contaminated with Chlorinated Ethenes

The reliable characterization of subsurface contamination of spatially extended contaminated sites is a challenging task, especially with an unknown history of land use. Conventional technologies often fail due to temporal and financial constraints and thus hinder the redevelopment of abandoned areas in particular. Here we compare two site screening techniques that can be applied quickly at relatively low cost, namely Direct Push (DP)‐based groundwater sampling and tree core sampling. The effectiveness of both methods is compared for a rural megasite contaminated with chlorinated hydrocarbons. Unexpected pollution hot spots could be identified using both of these methods, while tree coring even enabled the delineation of the contaminant plume flowing into an adjacent wetland inaccessible for DP units. Both methods showed a good agreement in revealing the spatial pattern of the contamination. The correlation between groundwater concentrations and equivalent concentrations in wood was linear and highly significant for trichloroethene. Correlation was less obvious for its metabolite cis‐dichloroethene, but still significant. As outcome of our study we recommend tree coring and for initial screening in combination with a DP sampling to retrieve quantitative data on groundwater pollutants in order to assess the contamination situation of a non‐ or only partly investigated site. The subsequent placement of monitoring wells for long‐term monitoring of contamination levels is recommended. A combination of methods would achieve more relevant information at comparable or possibly even lower efforts in comparison to a conventional site investigation.     

Conducting Slug Tests in Mini‐Piezometers

Slug tests performed using mini‐piezometers with internal diameters as small as 0.43 cm can provide a cost effective tool for hydraulic characterization. We evaluated the hydraulic properties of the apparatus in a laboratory environment and compared those results with field tests of mini‐piezometers installed into locations with varying hydraulic properties. Based on our evaluation, slug tests conducted in mini‐piezometers using the fabrication and installation approach described here are effective within formations where the hydraulic conductivity is less than 1 × 10^?3 cm/s. While these constraints limit the potential application of this method, the benefits to this approach are that the installation, measurement, and analysis is cost effective, and the installation can be completed in areas where other (larger diameter) methods might not be possible. Additionally, this methodology could be applied to existing mini‐piezometers previously installed for other purposes. Such analysis of existing installations could be beneficial in interpreting previously collected data (e.g., water‐quality data or hydraulic head data).     

Homotopy Perturbation Method‐Based Analytical Solution for Tide‐Induced Groundwater Fluctuations

The groundwater variations in unconfined aquifers are governed by the nonlinear Boussinesq's equation. Analytical solution for groundwater fluctuations in coastal aquifers under tidal forcing can be solved using perturbation methods. However, the perturbation parameters should be properly selected and predefined for traditional perturbation methods. In this study, a new dimensional, higher‐order analytical solution for groundwater fluctuations is proposed by using the homotopy perturbation method with a virtual perturbation parameter. Parameter‐expansion method is used to remove the secular terms generated during the solution process. The solution does not require any predefined perturbation parameter and valid for higher values of amplitude parameter A/D, where A is the amplitude of the tide and D is the aquifer thickness.     

Efficacy of Hollow‐Fiber Ultrafiltration for Microbial Sampling in Groundwater

The goal of this study was to test hollow‐fiber ultrafiltration as a method for concentrating in situ bacteria and viruses in groundwater samples. Water samples from nine wells tapping a shallow sandy aquifer in a densely populated village in Bangladesh were reduced in volume approximately 400‐fold using ultrafiltration. Culture‐based assays for total coliforms and Escherichia coli, as well as molecular‐based assays for E. coli, Bacteroides, and adenovirus, were used as microbial markers before and after ultrafiltration to evaluate performance. Ultrafiltration increased the concentration of the microbial markers in 99% of cases. However, concentration factors (CF = post‐filtration concentration/pre‐filtration concentration) for each marker calculated from geometric means ranged from 52 to 1018 compared to the expected value of 400. The efficiency was difficult to quantify because concentrations of some of the markers, especially E. coli and total coliforms, in the well water (WW) collected before ultrafiltration varied by several orders of magnitude during the period of sampling. The potential influence of colloidal iron oxide precipitates in the groundwater was tested by adding EDTA to the pre‐filtration water in half of the samples to prevent the formation of precipitates. The use of EDTA had no significant effect on the measurement of culturable or molecular markers across the 0.5 to 10 mg/L range of dissolved Fe²⁺ concentrations observed in the groundwater, indicating that colloidal iron did not hinder or enhance recovery or detection of the microbial markers. Ultrafiltration appears to be effective for concentrating microorganisms in environmental water samples, but additional research is needed to quantify losses during filtration.     

A Borehole Test for Chlorinated Solvent Diffusion and Degradation Rates in Sedimentary Rock

We present a new field measurement and numerical interpretation method (combined termed “test”) to parameterize the diffusion of trichloroethene (TCE) and its biodegradation products (DPs) from the matrix of sedimentary rock. The method uses a dual-packer system to interrogate a low-permeability section of the rock matrix adjacent to a previously contaminated borehole and uses the borehole monitoring history to establish the pretest condition. TCE and its DPs are removed from the groundwater between the packers at the onset of the testing. The parameters estimated by fitting a radial diffusion model to the concentration history and borehole concentration data, also termed back diffusion, are the tortuosity factor and sorption coefficients of TCE and DPs in the rock matrix and the TCE and DP biodegradation rate coefficients in the borehole. We demonstrate the equipment design and the interpretive method using a borehole accessing the gray mudstone at a TCE contaminated site in the Newark Basin. In this test, both nonreactive (bromide) and reactive (trichlorofluoroethene) tracers are used to constrain the estimated parameters; however, the bromide tracer was not needed to estimate the parameters in this test. The parameters estimated from the field test are consistent with values measured independently in laboratory experiments using field samples of similar lithology. From the interpretation, we compute the TCE and DP concentration distributions in the rock matrix prior to the test to illustrate how the results can be used to enhance understanding of contaminant distribution in the rock matrix.     

Accumulation of Pesticides in Anaerobic Clayey Till‐Controls and Implications for Groundwater

Pesticide residuals after point‐source pesticide spills in clay‐rich aquitards may potentially affect underlying groundwater for many decades due to slow release of accumulated pollution in the clayey matrix material of the aquitard. In this study, we evaluated factors behind different degrees of accumulation of phenoxy acids (MCPP, dichlorprop, 2,4‐D MCPA) and triazines (simazine and terbutylazine) observed in an old pesticide pollution described by Jørgensen et al. (2016a, this issue). By using leaching experiments, it was shown that a zone of maximum concentrations of MCPP and dichlorprop (mg/L) observed by Jørgensen et al. (2016a, this issue) represented accumulated potentially mobile pollution in anaerobic, however largely immobile pore water of the clayey matrix in the upper few meters of the unoxidized aquitard zone. By using sorption experiments, we determined 9 to16 times higher mobility by diffusion and flow for the phenoxy acids (R = 1 to 2) than for the triazines (R = 9 to 16) in the clayey matrix material of the aquitard. This indicated that more rapid and greater accumulation could occur for the phenoxy acids in the clayey matrix than for the triazines. In contrast, the relative mobility of the phenoxy acids and triazines was much closer in sand‐filled fractures and thin sand layers/lamina in the clay, suggesting that the migration of the same compounds along these textural preferential flow paths into the underlying aquifer was less different. Despite that a greater mass had originally been spilled of 2,4‐D and MCPA having similar mobility as the accumulated MCPP and dichlorprop, these compounds were not accumulated in the zone of maximum concentrations. It is suggested that the controls, which initially allowed for the observed separate accumulation of MCPP and dichlorprop as a zone of maximum pollution, were due to the combination of high persistence and high mobility for these specific pesticides in the clayey till matrix material of the aquitard. The investigation showed that over time the initial high concentrations of the accumulated phenoxy acids (MCPP, dichlorprop) transformed into high concentrations of related mobile degradation products (e.g., 4‐CPP and 2‐CPP), which extended the total time of groundwater pollution beyond the disappearance of the original phenoxyacids.     

Closed‐Form Approximations for Two‐Dimensional Groundwater Age Patterns in a Fresh Water Lens

Simple closed‐form approximations are presented for calculating the steady‐state groundwater age distribution in two‐dimensional vertical cross sections of idealized fresh water lenses overlying salt water, for aquifers that are vertically semi‐infinite and of finite thickness. The approximations are developed on the basis of existing one‐dimensional analytical solutions for travel‐time calculation in fresh water lenses and approximate streamline formulations. The two‐dimensional age distributions based on the closed‐form solutions match convincingly with numerical simulations. As expected, notable deviations from the numerical solution are encountered at the groundwater flow divide and when submarine groundwater discharge occurs. Ratios of recharge over hydraulic conductivities are varied to explore how the magnitude of the deviations changes, and it is found that the approximate closed‐form solutions perform well over a range of conditions found in natural systems.     

Total Dissolved Gas Pressure Measurements for Assessing Contaminant Degradation Processes in Groundwater

The total dissolved gas pressure (P_TDG ) probe has been used in groundwater studies for over a decade, but rarely in assessing contaminant degradation, despite the many degradation reactions that produce or consume dissolved gases. Here we present three studies to demonstrate the application of P_TDG measurements to groundwater experiencing contaminant degradation, with discussion of its benefits and limitations. The first study is a pilot‐scale laboratory experiment simulating dissolved ethanol contamination of an anaerobic sand aquifer. Continuous monitoring of P_TDG showed the rapid onset of microbial hydrocarbon degradation via denitrification and fermentation. The subsequent formation of a gas phase was revealed when P_TDG began mimicking the bubbling pressure (P_G *; sum of hydrostatic and atmospheric pressure), fluctuating with atmospheric pressure. Some deviations of P_TDG above P_G * occurred also, which may hold promise for signalling substantial changes in the rate or type of degradation process (here, the onset of methanogenesis). In the second study, synoptic field measurements at a petroleum plume site demonstrated how elevated P_TDG could identify wells with evidence of hydrocarbon degradation (denitrification and/or methanogenesis). And finally, combined field measurements of dissolved oxygen (DO) and P_TDG in monitoring wells of a nitrate‐contaminated aquifer (Abbottsford‐Sumas) revealed areas where denitrification was likely occurring. Limitations to P_TDG use identified in these studies included the masking of degradation processes by the presence of a gas phase, as when trapped following water table fluctuations or formed from rigorous degradation reactions, and confounded assessment of P_TDG patterns from other natural or anthropogenic processes that can also influence groundwater P_TDG .     

Estimating depth‐averaged velocities in rough ­channels

Profiles of streamwise velocity obtained from North Boulder Creek, Colorado, typically are non‐logarithmic in form and exhibit the strong influence of form drag associated with coarse bed roughness. The spatially averaged profile is consistent with recent theoretical profile forms suggested for rough channels that are based on a partitioning of the total stress between a fluid part and a part associated with form drag on bed particles. Estimates of local depth‐averaged velocity using algorithms that are based on several measurements in the flow column improve with explicit Riemann averaging, versus simple averaging, of the measurements. Estimates based on a single‐point measurement at 0·6 of the flow depth, assuming a logarithmic or approximately logarithmic velocity profile, are the least reliable. Copyright © 2000 John Wiley & Sons, Ltd.     

Recognition of Regional Water Table Patterns for Estimating Recharge Rates in Shallow Aquifers

We propose a new method for groundwater recharge rate estimation in regions with stream-aquifer interactions, at a linear scale on the order of 10 km and more. The method is based on visual identification and quantification of classically recognized water table contour patterns. Simple quantitative analysis of these patterns can be done manually from measurements on a map, or from more complex GIS data extraction and curve fitting. Recharge rate is then estimated from the groundwater table contour parameters, streambed gradients, and aquifer transmissivity using an analytical model for groundwater flow between parallel perennial streams. Recharge estimates were obtained in three regions (areas of 1500, 2200, and 3300 km²) using available water table maps produced by different methods at different times in the area of High Plains Aquifer in Nebraska. One region is located in the largely undeveloped Nebraska Sand Hills area, while the other two regions are located at a transition zone from Sand Hills to loess-covered area and include areas where groundwater is used for irrigation. Obtained recharge rates are consistent with other independent estimates. The approach is useful and robust diagnostic tool for preliminary estimates of recharge rates, evaluation of the quality of groundwater table maps, identification of priority areas for further aquifer characterization and expansion of groundwater monitoring networks prior to using more detailed methods.     

Tests on low‐ductility RC frames under high‐ and low‐frequency excitations

The response of low‐ductility reinforced concrete (RC) frames, designed typically for a non‐seismic region, subjected to two frequencies of base excitations is studied. Five half‐scaled, two‐bay, two‐storey, RC frames, each approximately 5 m wide by 3.3 m high, were subjected to both horizontal and/or vertical base excitations with a frequency of 40 Hz as well as a lower frequency of about 4 Hz (close to the fundamental frequency) using a shake table. The imposed acceleration amplitude ranged from 0.2 to 1.2g. The test results showed that the response characteristics of the structures differed under high‐ and low‐frequency excitations. The frames were able to sustain high‐frequency excitations without damage but were inadequate for low‐frequency excitations, even though the frames exhibited some ductility. Linear‐elastic time‐history analysis can predict reasonably well the structural response under high‐frequency excitations. As the frames were not designed for seismic loads, the reinforcement detailing may not have been adequate, based on the crack pattern observed. The effect of vertical excitation can cause significant additional forces in the columns and moment reversals in the beams. The ‘strong‐column, weak‐beam’ approach for lateral load RC frame design is supported by experimental observations. Copyright © 2001 John Wiley & Sons, Ltd.     

Using DNA‐Stable Isotope Probing to Identify MTBE‐ and TBA‐Degrading Microorganisms in Contaminated Groundwater

Although the anaerobic biodegradation of methyl tert‐butyl ether (MTBE) and tert‐butyl alcohol (TBA) has been documented in the laboratory and the field, knowledge of the microorganisms and mechanisms involved is still lacking. In this study, DNA‐stable isotope probing (SIP) was used to identify microorganisms involved in anaerobic fuel oxygenate biodegradation in a sulfate‐reducing MTBE and TBA plume. Microorganisms were collected in the field using Bio‐Sep® beads amended with ¹³C₅‐MTBE, ¹³C₁‐MTBE (only methoxy carbon labeled), or ¹³C₄‐TBA. ¹³C‐DNA and ¹²C‐DNA extracted from the Bio‐Sep beads were cloned and 16S rRNA gene sequences were used to identify the indigenous microorganisms involved in degrading the methoxy group of MTBE and the tert‐butyl group of MTBE and TBA. Results indicated that microorganisms were actively degrading ¹³C‐labeled MTBE and TBA in situ and the ¹³C was incorporated into their DNA. Several sequences related to known MTBE‐ and TBA‐degraders in the Burkholderiales and the Sphingomonadales orders were detected in all three ¹³C clone libraries and were likely to be primary degraders at the site. Sequences related to sulfate‐reducing bacteria and iron‐reducers, such as Geobacter and Geothrix, were only detected in the clone libraries where MTBE and TBA were fully labeled with ¹³C, suggesting that they were involved in processing carbon from the tert‐butyl group. Sequences similar to the Pseudomonas genus predominated in the clone library where only the methoxy carbon of MTBE was labeled with ¹³C. It is likely that members of this genus were secondary degraders cross‐feeding on ¹³C‐labeled metabolites such as acetate.     

Dual‐Screened Vertical Circulation Wells for Groundwater Lowering in Unconfined Aquifers

A new type of vertical circulation well (VCW) is used for groundwater dewatering at construction sites. This type of VCW consists of an abstraction screen in the upper part and an injection screen in the lower part of a borehole, whereby drawdown is achieved without net withdrawal of groundwater from the aquifer. The objective of this study is to evaluate the operation of such wells including the identification of relevant factors and parameters based on field data of a test site and comprehensive numerical simulations. The numerical model is able to delineate the drawdown of groundwater table, defined as free‐surface, by coupling the arbitrary Lagrangian–Eulerian algorithm with the groundwater flow equation. Model validation is achieved by comparing the field observations with the model results. Eventually, the influences of selected well operation and aquifer parameters on drawdown and on the groundwater flow field are investigated by means of parameter sensitivity analysis. The results show that the drawdown is proportional to the flow rate, inversely proportional to the aquifer conductivity, and almost independent of the aquifer anisotropy in the direct vicinity of the well. The position of the abstraction screen has a stronger effect on drawdown than the position of the injection screen. The streamline pattern depends strongly on the separation length of the screens and on the aquifer anisotropy, but not on the flow rate and the horizontal hydraulic conductivity.