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 703 records · Page 39Linked to original sources

Integration of stable carbon isotope, microbial community, dissolved hydrogen gas, and 2 H H 2 O tracer data to assess bioaugmentation for chlorinated ethene degradation in fractured rocks

An in situ bioaugmentation (BA) experiment was conducted to understand processes controlling microbial dechlorination of trichloroethene (TCE) in groundwater at the Naval Air Warfare Center (NAWC), West Trenton, NJ. In the BA experiment, an electron donor (emulsified vegetable oil and sodium lactate) and a chloro-respiring microbial consortium were injected into a well in fractured mudstone of Triassic age. Water enriched in 2 H was also injected as a tracer of the BA solution, to monitor advective transport processes. The changes in concentration and the δ 13 C of TCE, cis-dichloroethene (cis-DCE), and vinyl chloride (VC); the δ 2 H of water; changes in the abundance of the microbial communities; and the concentration of dissolved H2 gas compared to pre- test conditions, provided multiple lines of evidence that enhanced biodegradation occurred in the injection well and in two downgradient wells. For those wells where the biodegradation was stimulated intensively, the sum of the molar chlorinated ethene (CE) concentrations in post-BA water was higher than that of the sum of the pre-BA background molar CE concentrations. The concentration ratios of TCE/(cis-DCE + VC) indicated that the increase in molar CE concentration may result from additional TCE mobilized from the rock matrix in response to the oil injection or due to desorption/diffusion. The stable carbon isotope mass-balance calculations show that the weighted average 13 C isotope of the CEs was enriched for around a year compared to the background value in a two year monitoring period, an effective indication that dechlorination of VC was occurring. Insights gained from this study can be applied to efforts to use BA in other fractured rock systems. The study demonstrates that a BA approach can substantially enhance in situ bioremediation not only in fractures connected to the injection well, but also in the rock matrix around the well due to processes such as diffusion and desorption. Because the effect of the BA was intensive only in wells where an amendment was distributed during injection, it is necessary to adequately distribute the amendments throughout the fractured rock to achieve substantial bioremediation. The slowdown in BA effect after a year is due to some extend to the decrease abundant of appropriate microbes, but more likely the decreased concentration of electron donor.

New Jersey↗

A mass balance approach to investigating geochemical controls on secondary water quality impacts at a crude oil spill site near Bemidji, MN

Secondary water quality impacts can result from a broad range of coupled reactions triggered by primary groundwater contaminants. Data from a crude-oil spill research site near Bemidji, MN provide an ideal test case for investigating the complex interactions controlling secondary impacts, including depleted dissolved oxygen and elevated organic carbon, inorganic carbon, CH 4 , Mn, Fe, and other dissolved ions. To better understand these secondary impacts, this study began with an extensive data compilation of various data types, comprising aqueous, sediment, gas, and oil phases, covering a 260 m cross-sectional domain over 30 years. Mass balance calculations are used to quantify pathways that control secondary components, by using the data to constrain the sources and sinks for the important redox processes. The results show that oil constituents other than BTEX (benzene, toluene, ethylbenzene, o -, m - and p -xylenes), including n -alkanes and other aromatic compounds, play significant roles in plume evolution and secondary water quality impacts. The analysis underscores previous results on the importance of non-aqueous phases. Over 99.9% of the Fe 2+ plume is attenuated by immobilization on sediments as Fe(II) and 85–95% of the carbon biodegradation products are outgassed. Gaps identified in carbon and Fe mass balances and in pH buffering mechanisms are used to formulate a new conceptual model. This new model includes direct out-gassing of CH 4 and CO 2 from organic carbon biodegradation, dissolution of directly produced CO 2 , and sorption with H + exchange to improve pH buffering. The identification of these mechanisms extends understanding of natural attenuation of potential secondary impacts at enhanced reductive dechlorination sites, particularly for reduced Fe plumes, produced CH 4 , and pH perturbations.

Minnesota↗

Comparison of surficial CO2 efflux to other measures of subsurface crude oil degradation

At a spill site near Bemidji, Minnesota, crude oil at the water table has been undergoing anaerobic biodegradation for over 30 years. Previous work at this site has shown that methane produced from biodegradation of the oil migrates upward and is oxidized in a methanotrophic zone midway between the water table and the surface. To compare microbial activity measurement methods from multiple locations in the oil body, surficial carbon dioxide efflux, methanogen and methanotroph concentrations, and oil degradation state were collected. Carbon dioxide effluxes over the oil body averaged more than four times those at the background site. Methanotrophic bacteria concentrations measured using pmoA were over 10 5 times higher above the oil-contaminated sediments compared with the background site. Methanogenic archaea measured using mcrA ranged from 10 5 to over 10 7 in the oil and were below detection in the background. Methanogens correlated very well with methanotroph concentrations (r = 0.99), n-alkylcyclohexane losses as a proxy for degradation state (r = − 0.96), and somewhat less well with carbon dioxide efflux (r = 0.92). Carbon dioxide efflux similarly correlated to methanotroph concentrations (r = 0.90) and n-alkylcyclohexane losses (r = − 0.91).

Minnesota↗

High-resolution delineation of chlorinated volatile organic compounds in a dipping, fractured mudstone: depth- and strata-dependent spatial variability from rock-core sampling

Synthesis of rock-core sampling and chlorinated volatile organic compound (CVOC) analysis at five coreholes, with hydraulic and water-quality monitoring and a detailed hydrogeologic framework, was used to characterize the fine-scale distribution of CVOCs in dipping, fractured mudstones of the Lockatong Formation of Triassic age, of the Newark Basin in West Trenton, New Jersey. From these results, a refined conceptual model for more than 55 years of migration of CVOCs and depth- and strata-dependent rock-matrix contamination was developed. Industrial use of trichloroethene (TCE) at the former Naval Air Warfare Center (NAWC) from 1953 to 1995 resulted in dense non-aqueous phase liquid (DNAPL) TCE and dissolved TCE and related breakdown products, including other CVOCs, in underlying mudstones. Shallow highly weathered and fractured strata overlie unweathered, gently dipping, fractured strata that become progressively less fractured with depth. The unweathered lithology includes black highly fractured (fissile) carbon-rich strata, gray mildly fractured thinly layered (laminated) strata, and light-gray weakly fractured massive strata. CVOC concentrations in water samples pumped from the shallow weathered and highly fractured strata remain elevated near residual DNAPL TCE, but dilution by uncontaminated recharge, and other natural and engineered attenuation processes, have substantially reduced concentrations along flow paths removed from sources and residual DNAPL. CVOCs also were detected in most rock-core samples in source areas in shallow wells. In many locations, lower aqueous concentrations, compared to rock core concentrations, suggest that CVOCs are presently back-diffusing from the rock matrix. Below the weathered and highly fractured strata, and to depths of at least 50 meters (m), groundwater flow and contaminant transport is primarily in bedding-plane-oriented fractures in thin fissile high-carbon strata, and in fractured, laminated strata of the gently dipping mudstones. Despite more than 18 years of pump and treat (P&T) remediation, and natural attenuation processes, CVOC concentrations in aqueous samples pumped from these deeper strata remain elevated in isolated intervals. DNAPL was detected in one borehole during coring at a depth of 27 m. In contrast to core samples from the weathered zone, concentrations in core samples from deeper unweathered and unfractured strata are typically below detection. However, high CVOC concentrations were found in isolated samples from fissile black carbon-rich strata and fractured gray laminated strata. Aqueous-phase concentrations were correspondingly high in samples pumped from these strata via short-interval wells or packer-isolated zones in long boreholes. A refined conceptual site model considers that prior to P&T remediation groundwater flow was primarily subhorizontal in the higher-permeability near surface strata, and the bulk of contaminant mass was shallow. CVOCs diffused into these fractured and weathered mudstones. DNAPL and high concentrations of CVOCs migrated slowly down in deeper unweathered strata, primarily along isolated dipping bedding-plane fractures. After P&T began in 1995, using wells open to both shallow and deep strata, downward transport of dissolved CVOCs accelerated. Diffusion of TCE and other CVOCs from deeper fractures penetrated only a few centimeters into the unweathered rock matrix, likely due to sorption of CVOCs on rock organic carbon. Remediation in the deep, unweathered strata may benefit from the relatively limited migration of CVOCs into the rock matrix. Synthesis of rock core sampling from closely spaced boreholes with geophysical logging and hydraulic testing improves understanding of the controls on CVOC delineation and informs remediation design and monitoring.

New Jersey, New York, Pennsylvania↗

Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer-aquitard complexes

This study evaluates the role of the Peclet number as affected by molecular diffusion in transient anomalous transport, which is one of the major knowledge gaps in anomalous transport, by combining Monte Carlo simulations and stochastic model analysis. Two alluvial settings containing either short- or long-connected hydrofacies are generated and used as media for flow and transport modeling. Numerical experiments show that 1) the Peclet number affects both the duration of the power-law segment of tracer breakthrough curves (BTCs) and the transition rate from anomalous to Fickian transport by determining the solute residence time for a given low-permeability layer, 2) mechanical dispersion has a limited contribution to the anomalous characteristics of late-time transport as compared to molecular diffusion due to an almost negligible velocity in floodplain deposits, and 3) the initial source dimensions only enhance the power-law tail of the BTCs at short travel distances. A tempered stable stochastic (TSS) model is then applied to analyze the modeled transport. Applications show that the time-nonlocal parameters in the TSS model relate to the Peclet number, P e . In particular, the truncation parameter in the TSS model increases nonlinearly with a decrease in P e due to the decrease of the mean residence time, and the capacity coefficient increases with an increase in molecular diffusion which is probably due to the increase in the number of immobile particles. The above numerical experiments and stochastic analysis therefore reveal that the Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer–aquitard complexes.

Journal of Contaminant Hydrology↗

Relating subsurface temperature changes to microbial activity at a crude oil-contaminated site

Crude oil at a spill site near Bemidji, Minnesota has been undergoing aerobic and anaerobic biodegradation for over 30 years, creating a 150–200 m plume of primary and secondary contaminants. Microbial degradation generates heat that should be measurable under the right conditions. To measure this heat, thermistors were installed in wells in the saturated zone and in water-filled monitoring tubes in the unsaturated zone. In the saturated zone, a thermal groundwater plume originates near the residual oil body with temperatures ranging from 2.9 °C above background near the oil to 1.2 °C down gradient. Temperatures in the unsaturated zone above the oil body were up to 2.7 °C more than background temperatures. Previous work at this site has shown that methane produced from biodegradation of the oil migrates upward and is oxidized in a methanotrophic zone midway between the water table and the surface. Enthalpy calculations and observations demonstrate that the temperature increases primarily result from aerobic methane oxidation in the unsaturated zone above the oil. Methane oxidation rates at the site independently estimated from surface CO 2 efflux data are comparable to rates estimated from the observed temperature increases. The results indicate that temperature may be useful as a low-cost measure of activity but care is required to account for the correct heat-generating reactions, other heat sources and the effects of focused recharge.

Minnesota↗

Non-invasive flow path characterization in a mining-impacted wetland

Time-lapse electrical resistivity (ER) was used to capture the dilution of a seasonal pulse of acid mine drainage (AMD) contamination in the subsurface of a wetland downgradient of the abandoned Pennsylvania mine workings in central Colorado. Data were collected monthly from mid-July to late October of 2013, with an additional dataset collected in June of 2014. Inversion of the ER data shows the development through time of multiple resistive anomalies in the subsurface, which corroborating data suggest are driven by changes in total dissolved solids (TDS) localized in preferential flow pathways. Sensitivity analyses on a synthetic model of the site suggest that the anomalies would need to be at least several meters in diameter to be adequately resolved by the inversions. The existence of preferential flow paths would have a critical impact on the extent of attenuation mechanisms at the site, and their further characterization could be used to parameterize reactive transport models in developing quantitative predictions of remediation strategies.

Colorado↗

The role of alluvial aquifer sediments in attenuating a dissolved arsenic plume

In a crude-oil-contaminated sandy aquifer at the Bemidji site in northern Minnesota, biodegradation of petroleum hydrocarbons has resulted in release of naturally occurring As to groundwater under Fe-reducing conditions. This study used chemical extractions of aquifer sediments collected in 1993 and 2011–2014 to evaluate the relationship between Fe and As in different redox zones (oxic, methanogenic, Fe-reducing, anoxic-suboxic transition) of the contaminated aquifer over a twenty-year period. Results show that 1) the aquifer has the capacity to naturally attenuate the plume of dissolved As, primarily through sorption; 2) Fe and As are linearly correlated in sediment across all redox zones, and a regression analysis between Fe and As reasonably predicted As concentrations in sediment from 1993 using only Fe concentrations; 3) an As-rich “iron curtain,” associated with the anoxic-suboxic transition zone, migrated 30 m downgradient between 1993 and 2013 as a result of the hydrocarbon plume evolution; and 4) silt lenses in the aquifer preferentially sequester dissolved As, though As is remobilized into groundwater from sediment after reducing conditions are established. Using results of this study coupled with historical data, we develop a conceptual model which summarizes the natural attenuation of As and Fe over time and space that can be applied to other sites that experience As mobilization due to an influx of bioavailable organic matter.

Journal of Contaminant Hydrology↗

Porosity and pore size distribution in a sedimentary rock: Implications for the distribution of chlorinated solvents

Characterizing properties of the rock matrix that control retention and release of chlorinated solvents is essential in evaluating the extent of contamination and the application of remediation technologies in fractured rock. Core samples from seven closely spaced boreholes in a mudstone subject to trichloroethene (TCE) contamination were analyzed using Mercury Intrusion Porosimetry to investigate porosity and pore size distribution as a function of mudstone characteristics, and depth and lateral extent in the aquifer; organic carbon content was also evaluated to identify the potential for adsorption. Porosity and retardation factor varied over two orders of magnitude, with the largest porosities and largest retardation factors associated with carbon-rich mudstone layers. Larger porosities were also measured in the shallow rock that has been subject to enhanced groundwater flow. Porosity also varied over more than an order of magnitude in spatially continuous mudstone layers. The analyses of the rock cores indicated that the largest pore diameters may be accessible to entry of the nonaqueous form of TCE. Although the porosity associated with the largest pore diameters is small (~ 0.1%), that volume of TCE can significantly affect the total TCE that is retained in the rock matrix. The dimensions of the largest pore diameters may also be accessible to microbes responsible for reductive dechlorination; however, the small percentage of the pore space that can accommodate microbes may limit the extent of reductive dechlorination in the rock matrix.

Journal of Contaminant Hydrology↗

Relative contributions of microbial and infrastructure heat at a crude oil-contaminated site

Biodegradation of contaminants can increase the temperature in the subsurface due to heat generated from exothermic reactions, making temperature observations a potentially low-cost approach for determining microbial activity. For this technique to gain more widespread acceptance, it is necessary to better understand all the factors affecting the measured temperatures. Biodegradation has been occurring at a crude oil-contaminated site near Bemidji, Minnesota for 39 years, creating a quasi-steady-state plume of contaminants and degradation products. A model of subsurface heat generation and transport helps elucidate the contribution of microbial and infrastructure heating to observed temperature increases at this site. We created a steady-state, two-dimensional, heat transport model using previous-published parameter values for physical, chemical and biodegradation properties. Simulated temperature distributions closely match the observed average annual temperatures measured in the contaminated area at the site within less than 0.2 °C in the unsaturated zone and 0.4 °C in the saturated zone. The model results confirm that the observed subsurface heat from microbial activity is due primarily to methane oxidation in the unsaturated zone resulting in a 3.6 °C increase in average annual temperature. Another important source of subsurface heat is from the active, crude-oil pipelines crossing the site. The pipelines impact temperatures for a distance of 200 m and contribute half the heat. Model results show that not accounting for the heat from the pipelines leads to overestimating the degradation rates by a factor of 1.7, demonstrating the importance of identifying and quantifying all heat sources. The model results also highlighted a zone where previously unknown microbial activity is occurring at the site.

Minnesota↗

Variability of organic carbon content and the retention and release of trichloroethene in the rock matrix of a mudstone aquifer

Contaminants diffusing from fractures into the immobile porosity of the rock matrix are subject to prolonged residence times. Organic contaminants can adsorb onto organic carbonaceous materials in the matrix extending contaminant retention. An investigation of spatial variability of the fraction of organic carbon ( f oc ) is conducted on samples of rock core from seven closely spaced boreholes in a mudstone aquifer contaminated with trichloroethene (TCE). A total of 378 samples were analyzed at depths between 14 and 36 m below land surface. Mudstone units associated with deep water deposition have the largest f oc , with a maximum value of 0.0396, and units associated with shallow water deposition have the smallest f oc . Even though f oc correlates with depositional conditions, f oc still varies over more than an order of magnitude in continuous mudstone layers between boreholes, and there is large variability in f oc over short distances perpendicular to bedding. Simulations of diffusion and linear equilibrium adsorption of TCE using spatially variable f oc in the rock matrix show order of magnitude variability in the adsorbed TCE over short distances in the matrix and residence times extending to hundreds of years following remediation in adjacent fractures. Simulations using average values of f oc do not capture the range of TCE mass that can be retained in a rock matrix characterized by spatially variable f oc . Bounds on TCE mass within the rock matrix can be obtained by simulations with spatially uniform values of f oc equal to the maximum and minimum values of f oc for a given mudstone unit.

New Jersey↗

The complex spatial distribution of trichloroethene and the probability of NAPL occurrence in the rock matrix of a mudstone aquifer

Methanol extractions for chloroethene analyses are conducted on rock samples from seven closely spaced coreholes in a mudstone aquifer that was subject to releases of the nonaqueous phase liquid (NAPL) form of trichloroethene (TCE) between the 1950's and 1990's. Although TCE concentration in the rock matrix over the length of coreholes is dictated by proximity to subhorizontal bedding planefractures, elevated TCE concentrations in the rock matrix are not continuous along the most permeable bedding plane fractures. A complex configuration of subvertical and subhorizontal fractures appears to be responsible for the TCE distribution from prior TCE releases at land surface. Phase partitioning calculations of TCE in the rock matrix show that most TCE is adsorbed to solid surfaces because of the large fraction of organic carbon ( f oc ) in the mudstone. Large TCE content in some cores indicate the likely presence of the NAPL form of TCE in the rock matrix. Using average values of porosity ( n ) and f oc in phase partitioning calculations identifies a number of locations of possible NAPL occurrence in the rock matrix. Samples of mudstone analyzed for n and f oc show variability in these properties over several orders of magnitude. Accounting for this variability in phase partitioning calculations identifies a probability of NAPL occurrence, P NAPL . The spatial variability of P NAPL along coreholes identifies a configuration that may be attributed to a TCE source zone that has evolved after emplacement due to NAPL dissolution, adsorption, and matrix diffusion.

New Jersey↗

Sequential biodegradation of 1,2,4-trichlorobenzene at oxic-anoxic groundwater interfaces in model laboratory columns

Halogenated organic solvents such as chlorobenzenes (CBs) are frequent groundwater contaminants due to legacy spills. When contaminated anaerobic groundwater discharges into surface water through wetlands and other transition zones, aeration can occur from various physical and biological processes at shallow depths, resulting in oxic-anoxic interfaces (OAIs). This study investigated the potential for 1,2,4-trichlorobenzene (1,2,4-TCB) biodegradation at OAIs. A novel upflow column system was developed to create stable anaerobic and aerobic zones, simulating a natural groundwater OAI. Two columns containing (1) sand and (2) a mixture of wetland sediment and sand were operated continuously for 295 days with varied doses of 0.14-1.4 mM sodium lactate (NaLac) as a model electron donor. Both column matrices supported anaerobic reductive dechlorination and aerobic degradation of 1,2,4-TCB spatially separated between anaerobic and aerobic zones. Reductive dechlorination produced a mixture of di- and monochlorobenzene daughter products, with estimated zero-order dechlorination rates up to 31.3 µM/hr. Aerobic CB degradation, limited by available dissolved oxygen, occurred for 1,2,4-TCB and all dechlorinated daughter products. Initial reductive dechlorination did not enhance the overall observed extent or rate of subsequent aerobic CB degradation. Increasing NaLac dose increased the extent of reductive dechlorination, but suppressed aerobic CB degradation at 1.4 mM NaLac due to increased oxygen demand. 16S-rRNA sequencing of biofilm microbial communities revealed strong stratification of functional anaerobic and aerobic organisms between redox zones including the sole putative reductive dechlorinator detected in the columns, Dehalobacter. The sediment mixture column supported enhanced reductive dechlorination compared to the sand column at all tested NaLac doses and growth of Dehalobacter populations up to 4.1×108 copies/g (51% relative abundance), highlighting the potential benefit of sediments in reductive dechlorination processes. Results from these model systems suggest both substantial anaerobic and aerobic CB degradation can co-occur along the OAI at contaminated sites where bioavailable electron donors and oxygen are both present.

Journal of Contaminant Hydrology↗

Per- and polyfluoroalkyl substances (PFAS) mass flux and mass balance at an aqueous film-forming foam release site in semiarid eastern New Mexico, USA

Passive flux meters (PFMs) directly measure groundwater chemistry mass flux and Darcy flux, providing insight into contaminant source-zone architecture and transport properties. This study uses PFMs to characterize PFAS flux in groundwater at a semiarid site with a thick (greater than 90-m) unsaturated zone where groundwater has been contaminated with per - and polyfluoroalkyl substances (PFAS) related to the use of aqueous film-forming foam (AFFF) for fire training and fire suppression. PFAS mass discharge (PFAS mass flux integrated over a control plane) in groundwater downgradient from several PFAS release areas is calculated using PFM results. In groundwater downgradient from fire-training areas, total PFAS mass discharge (summed across 14 compounds) was estimated to be between 6.0 and 31 g per day in 2020 and between 5.9 and 23 g per day in 2021. Site-specific documentation, generic information on AFFF properties, and literature values of PFAS concentration in AFFF are used to estimate site-specific PFAS-application rates at fire-training areas. These PFAS-application rates are compared to groundwater PFAS-discharge rates. Results suggest that transformation processes (exact pathways unknown) have led to increased discharge of measured PFAS in groundwater relative to initial AFFF formulations. The mass balance approach has broad applicability as a high-level approach that can provide insight into PFAS transport at AFFF sites.

New Mexico↗

Contaminated stormwater sediment source tracking for polychlorinated biphenyls in an urban watershed of the Chesapeake Bay, United States

Fine-grained sediment in stormwater acts as a vector for persistent organic pollutants, like polychlorinated biphenyls (PCBs), through mobilization from sources within drainage areas of impacted urban watersheds. This study implemented a novel approach to identify the relative contributions of various landscape and stream sources of sediment from the Back River watershed in eastern Baltimore, Maryland, and investigated the applicability of using trace PCBs found in an urban environment as discriminants between each source type. Trace PCBs were found to be poor discriminants when identifying the relative sediment contributions of watershed-scale land use categories. When excluding PCBs in the development of a sediment fingerprinting model and instead utilizing trace elements and carbon only, sediment fingerprint modeling successfully differentiated green spaces and eroding streambanks as the most significant contributors to stormwaters sediment (37.1 % and 44.0 %, respectively) of the total sediment contributions of all considered source categories. In all samples collected from various landscape sources, storms, and cores detectable concentrations of PCBs were measured. The results of this study indicate that sediment fingerprinting may not be an effective method in predicting where PCBs may be found within an impacted watershed.

Maryland↗

Distinguishing natural from mining-related metal sources by including streambank groundwater data in a stream mass loading study

Distinguishing stream metal loading caused by mine features from that caused by natural background sources remains challenging, yet this distinction is essential for making effective remedial decisions at many legacy mine sites. We combine a stream tracer injection and synoptic sampling study with data from shallow near-stream groundwater wells to estimate left-bank versus right-bank metal loading contributions at the 100-m spatial scale. The study was performed in the East Mancos River, a mountain headwater stream in Colorado, USA. The dominant source of elevated stream metal concentrations could be either groundwater infiltration through right-bank Doyle Mine waste piles or natural acid rock drainage from hydrothermally altered bedrock located mainly on the left bank. For the five metals of concern (Cu, Al, Zn, Cd, and Mn), we find that 15 % of the load contributed by diffuse groundwater inputs in the section potentially influenced by Doyle mine waste originates from the right bank. This right-bank potential mine contribution equates to only 3 % of the total watershed instream load for these metals. Furthermore, apparent 3 H/ 3 He groundwater ages in segments contributing most of the right-bank metal loading are sufficiently old (9–12 yr) to suggest that infiltration through the waste piles, located only 140–180 m from the stream, is unlikely. Estimated potential Doyle mine loading contributions can therefore be considered maximum values. Study results thus indicate that Doyle mine waste piles are a minor source of metal loading under low-flow conditions, and streambank groundwater data can provide valuable additional information in stream mass loading studies.

Colorado↗

A functional relation for field-scale nonaqueous phase liquid dissolution developed using a pore network model

A pore network model with cubic chambers and rectangular tubes was used to estimate the nonaqueous phase liquid (NAPL) dissolution rate coefficient, K diss a i , and NAPL/water total specific interfacial area, a i . K diss a i was computed as a function of modified Peclet number( Pe ′) for various NAPL saturations ( S N ) and a i during drainage and imbibition and during dissolution without displacement. The largest contributor to a i was the interfacial area in the water-filled corners of chambers and tubes containing NAPL. When K diss a i was divided by a i , the resulting curves of dissolution coefficient, K diss versus Pe ′ suggested that an approximate value of K diss could be obtained as a weak function of hysteresis or S N . Spatially and temporally variable maps of K diss a i calculated using the network model were used in field-scale simulations of NAPL dissolution. These simulations were compared to simulations using a constant value of K diss a i and the empirical correlation of Powers et al. [Water Resour. Res. 30(2) (1994b) 321]. Overall, a methodology was developed for incorporating pore-scale processes into field-scale prediction of NAPL dissolution.

Journal of Contaminant Hydrology↗

Model coupling intraparticle diffusion/sorption, nonlinear sorption, and biodegradation processes

Diffusion, sorption and biodegradation are key processes impacting the efficiency of natural attenuation. While each process has been studied individually, limited information exists on the kinetic coupling of these processes. In this paper, a model is presented that couples nonlinear and nonequilibrium sorption (intraparticle diffusion) with biodegradation kinetics. Initially, these processes are studied independently (i.e., intraparticle diffusion, nonlinear sorption and biodegradation), with appropriate parameters determined from these independent studies. Then, the coupled processes are studied, with an initial data set used to determine biodegradation constants that were subsequently used to successfully predict the behavior of a second data set. The validated model is then used to conduct a sensitivity analysis, which reveals conditions where biodegradation becomes desorption rate-limited. If the chemical is not pre-equilibrated with the soil prior to the onset of biodegradation, then fast sorption will reduce aqueous concentrations and thus biodegradation rates. Another sensitivity analysis demonstrates the importance of including nonlinear sorption in a coupled diffusion/sorption and biodegradation model. While predictions based on linear sorption isotherms agree well with solution concentrations, for the conditions evaluated this approach overestimates the percentage of contaminant biodegraded by as much as 50%. This research demonstrates that nonlinear sorption should be coupled with diffusion/sorption and biodegradation models in order to accurately predict bioremediation and natural attenuation processes. To our knowledge this study is unique in studying nonlinear sorption coupled with intraparticle diffusion and biodegradation kinetics with natural media.

Journal of Contaminant Hydrology↗