USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37Linked to original sources

Distribution and transport of selected anthropogenic organic compounds on Mississippi River suspended sediment (U.S.A.), May/June 1988

The distribution and transport of selected hydrophobic halogenated organic compounds on suspended sediment from the lower Mississippi River were determined using discharge-weighted sampling with concurrent discharge measurements. Trends in compound concentration from upstream to downstream and the effects of selected tributaries were determined. The compounds identified on the suspended sediment include pentachlorobenzene, hexachlorobenzene, pentachloroanisole, dacthal, chlordane (trans-), nonachlor (trans-), chlorthalonil, and penta-, hexa-, hepta- and octachlorobiphenyls. Loads of most of the compounds increased from upstream to downstream on the main stem of the Mississippi River. Of the tributaries studied, the Ohio River had the most significant effect on the loads.The distribution and transport of selected hydrophobic halogenated organic compounds on suspended sediment from the lower Mississippi River were determined using discharge-weighted sampling with concurrent discharge measurements. Trends in compound concentration from upstream to downstream and the effects of selected tributaries were determined. The compounds identified on the suspended sediment include pentachlorobenzene, hexachlorobenzene, pentachloroanisole, dacthal, chlordane (trans-), nonachlor (trans-), chlorthalonil, and penta-, hexa-, hepta- and octachlorobiphenyls. Loads of most of the compounds increased from upstream to downstream on the main stem of the Mississippi River. Of the tributaries studied, the Ohio River had the most significant effect on the loads.

Journal of Contaminant Hydrology↗

A top specified boundary layer (TSBL) approximation approach for the simulation of groundwater contamination processes

This paper presents improvements in the 'classical boundary layer' (CBL) approximation method to obtain simple but robust initial characterization of aquifer contamination processes. Contaminants are considered to penetrate into the groundwater through the free surface of the aquifer. The improved method developed in this study is termed the 'top specified boundary layer' (TSBL) approach. It involves the specification of the contaminant concentration at the top of the contaminated 'region of interest' (ROI), which is simulated as a boundary layer. the TSBL modification significantly improves the ability of the boundary layer method to predict the development of concentration profiles over both space and time. The TSBL method can be useful for the simulation of cases in which the contaminant concentration is prescribed at the aquifer's free surface as well as for cases in which the contaminant mass flux is prescribed at the surface.

Journal of Contaminant Hydrology↗

Characterization of metal adsorption variability in a sand and gravel aquifer, Cape Cod, Massachusetts, U.S.A

Several geochemical properties of an aquifer sediment that control metal-ion adsorption were investigated to determine their potential use as indicators of the spatial variability of metal adsorption. Over the length of a 4.5-m-long core from a sand and gravel aquifer, lead (Pb2+) and zinc (Zn2+) adsorption at constant chemical conditions (pH 5.3) varied by a factor of 2 and 4, respectively. Pb2+ and Zn2+ were adsorbed primarily by Fe- and Al-oxide coatings on quartz-grain surfaces. Per unit surface area, both Pb2+ and Zn2+ adsorption were significantly correlated with the amount of Fe and Al that dissolved from the aquifer material in a partial chemical extraction. The variability in conditional binding constants for Pb2+ and Zn2+ adsorption (log KADS) derived from a simple non-electrostatic surface complexation model were also predicted by extracted Fe and Al normalized to surface area. Because the abundance of Fe- and Al-oxide coatings that dominate adsorption does not vary inversely with grain size by a simple linear relationship, only a weak, negative correlation was found between the spatial variability of Pb2+ adsorption and grain size in this aquifer. The correlation between Zn2+ adsorption and grain size was not significant. Partial chemical extractions combined with surface-area measurements have potential use for estimating metal adsorption variability in other sand and gravel aquifers of negligible carbonate and organic carbon content.

Massachusetts↗

Effect of cell physicochemical characteristics and motility on bacterial transport in groundwater

The influence of physicochemical characteristics and motility on bacterial transport in groundwater were examined in flow-through columns. Four strains of bacteria isolated from a crystalline rock groundwater system were investigated, with carboxylate-modified and amidine-modified latex microspheres and bromide as reference tracers. The bacterial isolates included a gram-positive rod (ML1), a gram-negative motile rod (ML2), a nonmotile mutant of ML2 (ML2m), and a gram-positive coccoid (ML3). Experiments were repeated at two flow velocities, in a glass column packed with glass beads, and in another packed with iron-oxyhydroxide coated glass beads. Bacteria breakthrough curves were interpreted using a transport equation that incorporates a sorption model from microscopic observation of bacterial deposition in flow-cell experiments. The model predicts that bacterial desorption rate will decrease exponentially with the amount of time the cell is attached to the solid surface. Desorption kinetics appeared to influence transport at the lower flow rate, but were not discernable at the higher flow rate. Iron-oxyhydroxide coatings had a lower-than-expected effect on bacterial breakthrough and no effect on the microsphere recovery in the column experiments. Cell wall type and shape also had minor effects on breakthrough. Motility tended to increase the adsorption rate, and decrease the desorption rate. The transport model predicts that at field scale, desorption rate kinetics may be important to the prediction of bacteria transport rates. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

Direct-push geochemical profiling for assessment of inorganic chemical heterogeneity in aquifers

Discrete-depth sampling of inorganic groundwater chemistry is essential for a variety of site characterization activities. Although the mobility and rapid sampling capabilities of direct-push techniques have led to their widespread use for evaluating the distribution of organic contaminants, complementary methods for the characterization of spatial variations in geochemical conditions have not been developed. In this study, a direct-push-based approach for high-resolution inorganic chemical profiling was developed at a site where sharp chemical contrasts and iron-reducing conditions had previously been observed. Existing multilevel samplers (MLSs) that span a fining-upward alluvial sequence were used for comparison with the direct-push profiling. Chemical profiles obtained with a conventional direct-push exposed-screen sampler differed from those obtained with an adjacent MLS because of sampler reactivity and mixing with water from previous sampling levels. The sampler was modified by replacing steel sampling components with stainless-steel and heat-treated parts, and adding an adapter that prevents mixing. Profiles obtained with the modified approach were in excellent agreement with those obtained from an adjacent MLS for all constituents and parameters monitored (Cl, NO3, Fe, Mn, DO, ORP, specific conductance and pH). Interpretations of site redox conditions based on field-measured parameters were supported by laboratory analysis of dissolved Fe. The discrete-depth capability of this approach allows inorganic chemical variations to be described at a level of detail that has rarely been possible. When combined with the mobility afforded by direct-push rigs and on-site methods of chemical analysis, the new approach is well suited for a variety of interactive site-characterization endeavors. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

Degradation of 1,1,2,2-tetrachloroethane and accumulation of vinyl chloride in wetland sediment microcosms and in situ porewater: Biogeochemical controls and associations with microbial communities

The biodegradation pathways of 1,1,2,2-tetrachloroethane (TeCA) and 1,1,2-trichloroethane (112TCA) and the associated microbial communities in anaerobic wetland sediments were evaluated using concurrent geochemical and genetic analyses over time in laboratory microcosm experiments. Experimental results were compared to in situ porewater data in the wetland to better understand the factors controlling daughter product distributions in a chlorinated solvent plume discharging to a freshwater tidal wetland at Aberdeen Proving Ground, Maryland. Microcosms constructed with wetland sediment from two sites showed little difference in the initial degradation steps of TeCA, which included simultaneous hydrogenolysis to 112TCA and dichloroelimination to 1,2-dichloroethene (12DCE). The microcosms from the two sites showed a substantial difference, however, in the relative dominance of subsequent dichloroelimination of 112TCA. A greater dominance of 112TCA dichloroelimination in microcosms constructed with sediment that was initially iron-reducing and subsequently simultaneously iron-reducing and methanogenic caused approximately twice as much vinyl chloride (VC) production as microcosms constructed with sediment that was methanogenic only throughout the incubation. The microcosms with higher VC production also showed substantially more rapid VC degradation. Field measurements of redox-sensitive constituents, TeCA, and its anaerobic degradation products along flowpaths in the wetland porewater also showed greater production and degradation of VC with concurrent methanogenesis and iron reduction. Molecular fingerprinting indicated that bacterial species [represented by a peak at a fragment size of 198 base pairs (bp) by MnlI digest] are associated with VC production from 112TCA dichloroelimination, whereas methanogens (190 and 307 bp) from the Methanococcales or Methanobacteriales family are associated with VC production from 12DCE hydrogenolysis. Acetate-utilizing methanogens (acetotrophs) appear to be involved in the biodegradation of VC. The relative abundance of Methanosarcinaceae, the only methanogen group with acetotrophic members, doubled in microcosms in which degradation of VC was observed. In addition, molecular analyses using primers specific for known dehalorespiring bacteria in the Dehalococcoides and Desulfuromonas groups showed the presence of these bacteria in microcosm slurry from the site that showed the highest VC production and degradation. Determination of biogeochemical controls and microbial consortia involved in TeCA degradation is leading to a better understanding of the heterogeneity in biodegradation rates and daughter product distribution in the wetland, improving capabilities for developing remediation and monitoring plans.

Journal of Contaminant Hydrology↗

Mechanisms of electron acceptor utilization: Implications for simulating anaerobic biodegradation

Simulation of biodegradation reactions within a reactive transport framework requires information on mechanisms of terminal electron acceptor processes (TEAPs). In initial modeling efforts, TEAPs were approximated as occurring sequentially, with the highest energy-yielding electron acceptors (e.g. oxygen) consumed before those that yield less energy (e.g., sulfate). Within this framework in a steady state plume, sequential electron acceptor utilization would theoretically produce methane at an organic-rich source and Fe(II) further downgradient, resulting in a limited zone of Fe(II) and methane overlap. However, contaminant plumes often display much more extensive zones of overlapping Fe(II) and methane. The extensive overlap could be caused by several abiotic and biotic processes including vertical mixing of byproducts in long-screened monitoring wells, adsorption of Fe(II) onto aquifer solids, or microscale heterogeneity in Fe(III) concentrations. Alternatively, the overlap could be due to simultaneous utilization of terminal electron acceptors. Because biodegradation rates are controlled by TEAPs, evaluating the mechanisms of electron acceptor utilization is critical for improving prediction of contaminant mass losses due to biodegradation. Using BioRedox-MT3DMS, a three-dimensional, multi-species reactive transport code, we simulated the current configurations of a BTEX plume and TEAP zones at a petroleum- contaminated field site in Wisconsin. Simulation results suggest that BTEX mass loss due to biodegradation is greatest under oxygen-reducing conditions, with smaller but similar contributions to mass loss from biodegradation under Fe(III)-reducing, sulfate-reducing, and methanogenic conditions. Results of sensitivity calculations document that BTEX losses due to biodegradation are most sensitive to the age of the plume, while the shape of the BTEX plume is most sensitive to effective porosity and rate constants for biodegradation under Fe(III)-reducing and methanogenic conditions. Using this transport model, we had limited success in simulating overlap of redox products using reasonable ranges of parameters within a strictly sequential electron acceptor utilization framework. Simulation results indicate that overlap of redox products cannot be accurately simulated using the constructed model, suggesting either that Fe(III) reduction and methanogenesis are occurring simultaneously in the source area, or that heterogeneities in Fe(III) concentration and/or mineral type cause the observed overlap. Additional field, experimental, and modeling studies will be needed to address these questions. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

Field investigation into unsaturated flow and transport in a fault: Model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model. ?? 2004 Published by Elsevier B.V.

Journal of Contaminant Hydrology↗

Numerical methods for improving sensitivity analysis and parameter estimation of virus transport simulated using sorptive-reactive processes

Using one- and two-dimensional homogeneous simulations, this paper addresses challenges associated with sensitivity analysis and parameter estimation for virus transport simulated using sorptive-reactive processes. Head, flow, and conservative- and virus-transport observations are considered. The paper examines the use of (1) observed-value weighting, (2) breakthrough-curve temporal moment observations, and (3) the significance of changes in the transport time-step size. The results suggest that (1) sensitivities using observed-value weighting are more susceptible to numerical solution variability, (2) temporal moments of the breakthrough curve are a more robust measure of sensitivity than individual conservative-transport observations, and (3) the transport-simulation time step size is more important than the inactivation rate in solution and about as important as at least two other parameters, reflecting the ease with which results can be influenced by numerical issues. The approach presented allows more accurate evaluation of the information provided by observations for estimation of parameters and generally improves the potential for reasonable parameter-estimation results. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

Reply to “Commentary: Assessment of past infiltration fluxes through Yucca Mountain on the basis of the secondary mineral record—Is it a viable methodology?” by Y.V. Dublyansky and S.Z. Smirnov

Many of the comments by Dublyansky and Smirnov (2005) on Marshall et al. (2003) reflect a longstanding debate over the origin of secondary calcite and opal deposits found in cavities and on fracture surfaces at Yucca Mountain, Nevada, site of a proposed high-level nuclear waste repository (US Department of Energy, 2001). These comments require consideration of data and interpretations beyond the scope of Marshall et al. (2003). Dublyansky et al. (2004) and Dublyansky et al. (2005) also have commented on papers published by Whelan et al. (2002) and Wilson et al. (2003), and we will refer to the replies to those comments (Whelan et al., 2004, Wilson and Cline, 2005) in addressing the comments that go beyond the scope of Marshall et al. (2003).

Nevada↗

Parameter and observation importance in modelling virus transport in saturated porous media - Investigations in a homogenous system

This paper evaluates the importance of seven types of parameters to virus transport: hydraulic conductivity, porosity, dispersivity, sorption rate and distribution coefficient (representing physical-chemical filtration), and in-solution and adsorbed inactivation (representing virus inactivation). The first three parameters relate to subsurface transport in general while the last four, the sorption rate, distribution coefficient, and in-solution and adsorbed inactivation rates, represent the interaction of viruses with the porous medium and their ability to persist. The importance of four types of observations to estimate the virus-transport parameters are evaluated: hydraulic heads, flow, temporal moments of conservative-transport concentrations, and virus concentrations. The evaluations are conducted using one- and two-dimensional homogeneous simulations, designed from published field experiments, and recently developed sensitivity-analysis methods. Sensitivity to the transport-simulation time-step size is used to evaluate the importance of numerical solution difficulties. Results suggest that hydraulic conductivity, porosity, and sorption are most important to virus-transport predictions. Most observation types provide substantial information about hydraulic conductivity and porosity; only virus-concentration observations provide information about sorption and inactivation. The observations are not sufficient to estimate these important parameters uniquely. Even with all observation types, there is extreme parameter correlation between porosity and hydraulic conductivity and between the sorption rate and in-solution inactivation. Parameter estimation was accomplished by fixing values of porosity and in-solution inactivation.

Journal of Contaminant Hydrology↗

Simulation of solute transport across low-permeability barrier walls

Low-permeability, non-reactive barrier walls are often used to contain contaminants in an aquifer. Rates of solute transport through such barriers are typically many orders of magnitude slower than rates through the aquifer. Nevertheless, the success of remedial actions may be sensitive to these low rates of transport. Two numerical simulation methods for representing low-permeability barriers in a finite-difference groundwater-flow and transport model were tested. In the first method, the hydraulic properties of the barrier were represented directly on grid cells and in the second method, the intercell hydraulic-conductance values were adjusted to approximate the reduction in horizontal flow, allowing use of a coarser and computationally efficient grid. The alternative methods were tested and evaluated on the basis of hypothetical test problems and a field case involving tetrachloroethylene (PCE) contamination at a Superfund site in New Hampshire. For all cases, advective transport across the barrier was negligible, but preexisting numerical approaches to calculate dispersion yielded dispersive fluxes that were greater than expected. A transport model (MODFLOW-GWT) was modified to (1) allow different dispersive and diffusive properties to be assigned to the barrier than the adjacent aquifer and (2) more accurately calculate dispersion from concentration gradients and solute fluxes near barriers. The new approach yields reasonable and accurate concentrations for the test cases.

Journal of Contaminant Hydrology↗

Recharge processes drive sulfate reduction in an alluvial aquifer contaminated with landfill leachate

Natural attenuation of contaminants in groundwater depends on an adequate supply of electron acceptors to stimulate biodegradation. In an alluvial aquifer contaminated with leachate from an unlined municipal landfill, the mechanism of recharge infiltration was investigated as a source of electron acceptors. Water samples were collected monthly at closely spaced intervals in the top 2 m of the saturated zone from a leachate-contaminated well and an uncontaminated well, and analyzed for δ 18 O, δ 2 H, non-volatile dissolved organic carbon (NVDOC), SO 4 2− , NO 3 − and Cl − . Monthly recharge amounts were quantified using the offset of the δ 18 O or δ 2 H from the local meteoric water line as a parameter to distinguish water types, as evaporation and methanogenesis caused isotopic enrichment in waters from different sources. Presence of dissolved SO 4 2− in the top 1 to 2 m of the saturated zone was associated with recharge; SO 4 2− averaged 2.2 mM, with maximum concentrations of 15 mM. Nitrate was observed near the water table at the contaminated site at concentrations up to 4.6 mM. Temporal monitoring of δ 2 H and SO 4 2− showed that vertical transport of recharge carried SO 4 2− to depths up to 1.75 m below the water table, supplying an additional electron acceptor to the predominantly methanogenic leachate plume. Measurements of δ 34 S in SO 4 2− indicated both SO 4 2− reduction and sulfide oxidation were occurring in the aquifer. Depth-integrated net SO 4 2− reduction rates, calculated using the natural Cl − gradient as a conservative tracer, ranged from 7.5 × 10 − 3 to 0.61 mM·d − 1 (over various depth intervals from 0.45 to 1.75 m). Sulfate reduction occurred at both the contaminated and uncontaminated sites; however, median SO 4 2− reduction rates were higher at the contaminated site. Although estimated SO 4 2− reduction rates are relatively high, significant decreases in NVDOC were not observed at the contaminated site. Organic compounds more labile than the leachate NVDOC may be present in the root zone, and SO 4 2− reduction may be coupled to methane oxidation. The results show that sulfur (and possibly nitrogen) redox processes within the top 2 m of the aquifer are directly related to recharge timing and seasonal water level changes in the aquifer. The results suggest that SO 4 2− reduction associated with the infiltration of recharge may be a significant factor affecting natural attenuation of contaminants in alluvial aquifers.

Journal of Contaminant Hydrology↗

Investigating the role of gas bubble formation and entrapment in contaminated aquifers: Reactive transport modelling

In many natural and contaminated aquifers, geochemical processes result in the production or consumption of dissolved gases. In cases where methanogenesis or denitrification occurs, the production of gases may result in the formation and growth of gas bubbles below the water table. Near the water table, entrapment of atmospheric gases during water table rise may provide a significant source of O 2 to waters otherwise depleted in O 2 . Furthermore, the presence of bubbles will affect the hydraulic conductivity of an aquifer, resulting in changes to the groundwater flow regime. The interactions between physical transport, biogeochemical processes, and gas bubble formation, entrapment and release is complex and requires suitable analysis tools. The objective of the present work is the development of a numerical model capable of quantitatively assessing these processes. The multicomponent reactive transport code MIN3P has been enhanced to simulate bubble growth and contraction due to in-situ gas production or consumption, bubble entrapment due to water table rise and subsequent re-equilibration of the bubble with ambient groundwater, and permeability changes due to trapped gas phase saturation. The resulting formulation allows for the investigation of complex geochemical systems where microbially mediated redox reactions both produce and consume gases as well as affect solution chemistry, alkalinity, and pH. The enhanced model has been used to simulate processes in a petroleum hydrocarbon contaminated aquifer where methanogenesis is an important redox process. The simulations are constrained by data from a crude oil spill site near Bemidji, MN. Our results suggest that permeability reduction in the methanogenic zone due to in-situ formation of gas bubbles, and dissolution of entrapped atmospheric bubbles near the water table, both work to attenuate the dissolved gas plume emanating from the source zone. Furthermore, the simulations demonstrate that under the given conditions more than 50% of all produced CH 4 partitions to the gas phase or is aerobically oxidised near the water table, suggesting that these processes should be accounted for when assessing the rate and extent of methanogenic degradation of hydrocarbons.

Journal of Contaminant Hydrology↗

Analysis of a mesoscale infiltration and water seepage test in unsaturated fractured rock: Spatial variabilities and discrete fracture patterns

A mesoscale (21??m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3??m ?? 4??m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten ?? parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum-based modeling of unsaturated flow in fractured rock at mesoscale or a larger scale is not necessarily conditional explicitly on discrete fracture patterns. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

Tracer test with As(V) under variable redox conditions controlling arsenic transport in the presence of elevated ferrous iron concentrations

To study transport and reactions of arsenic under field conditions, a small-scale tracer test was performed in an anoxic, iron-reducing zone of a sandy aquifer at the USGS research site on Cape Cod, Massachusetts, USA. For four weeks, a stream of groundwater with added As(V) (6.7 μM) and bromide (1.6 mM), was injected in order to observe the reduction of As(V) to As(III). Breakthrough of bromide (Br − ), As(V), and As(III) as well as additional parameters characterizing the geochemical conditions was observed at various locations downstream of the injection well over a period of 104 days. After a short lag period, nitrate and dissolved oxygen from the injectate oxidized ferrous iron and As(V) became bound to the freshly formed hydrous iron oxides. Approximately one week after terminating the injection, anoxic conditions had been reestablished and increases in As(III) concentrations were observed within 1 m of the injection. During the observation period, As(III) and As(V) were transported to a distance of 4.5 m downgradient indicating significant retardation by sorption processes for both species. Sediment assays as well as elevated concentrations of hydrogen reflected the presence of As(V) reducing microorganisms. Thus, microbial As(V) reduction was thought to be one major process driving the release of As(III) during the tracer test in the Cape Cod aquifer.

Journal of Contaminant Hydrology↗

Insights into the use of time-lapse GPR data as observations for inverse multiphase flow simulations of DNAPL migration

Perchloroethylene (PCE) saturations determined from GPR surveys were used as observations for inversion of multiphase flow simulations of a PCE injection experiment (Borden 9??m cell), allowing for the estimation of optimal bulk intrinsic permeability values. The resulting fit statistics and analysis of residuals (observed minus simulated PCE saturations) were used to improve the conceptual model. These improvements included adjustment of the elevation of a permeability contrast, use of the van Genuchten versus Brooks-Corey capillary pressure-saturation curve, and a weighting scheme to account for greater measurement error with larger saturation values. A limitation in determining PCE saturations through one-dimensional GPR modeling is non-uniqueness when multiple GPR parameters are unknown (i.e., permittivity, depth, and gain function). Site knowledge, fixing the gain function, and multiphase flow simulations assisted in evaluating non-unique conceptual models of PCE saturation, where depth and layering were reinterpreted to provide alternate conceptual models. Remaining bias in the residuals is attributed to the violation of assumptions in the one-dimensional GPR interpretation (which assumes flat, infinite, horizontal layering) resulting from multidimensional influences that were not included in the conceptual model. While the limitations and errors in using GPR data as observations for inverse multiphase flow simulations are frustrating and difficult to quantify, simulation results indicate that the error and bias in the PCE saturation values are small enough to still provide reasonable optimal permeability values. The effort to improve model fit and reduce residual bias decreases simulation error even for an inversion based on biased observations and provides insight into alternate GPR data interpretations. Thus, this effort is warranted and provides information on bias in the observation data when this bias is otherwise difficult to assess. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

The chemical response of particle-associated contaminants in aquatic sediments to urbanization in New England, U.S.A.

Relations between urbanization and particle-associated contaminants in New England were evaluated using a combination of samples from sediment cores, streambed sediments, and suspended stream sediments. Concentrations of PAHs, PCBs, DDT, and seven trace metals (Cd, Cr, Cu, Hg, Ni, Pb, Zn) were correlated strongly with urbanization, with the strongest relations to percentage commercial, industrial, and transportation (CIT) land use. Average PAH and metal concentrations in the most urbanized watersheds were approximately 30 and 6 times the reference concentrations, respectively, in remote, undeveloped watersheds. One-quarter to one-half of sampling sites had concentrations of PAHs, Cu, Pb, or Zn above the probable effects concentration, a set of sediment quality guidelines for adverse effects to aquatic biota, and sediments were predicted to be toxic, on average, when CIT land use exceeded about 10%. Trends in metals in cores from urban watersheds were dominantly downward, whereas trends in PAHs in a suburban watershed were upward. A regional atmospheric-fallout gradient was indicated by as much as order-of-magnitude-greater concentrations and accumulation rates of contaminants in cores from an undeveloped reference lake in Boston compared to those from remote reference watersheds. Contaminant accumulation rates in the lakes with urbanization in their watersheds, however, were 1-3 orders of magnitude greater than those of reference lakes, which indicate the dominance of local sources and fluvial transport of contaminants to urban lakes. These analyses demonstrate the magnitude of urban contamination of aquatic systems and air sheds, and suggest that, despite reductions in contaminant emissions in urban settings, streams and lakes will decline in quality as urbanization of their watersheds takes place. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗