Search USGS⌕ Search

Geology topics

H.I. Essaid

Publications and source records attributed to H.I. Essaid.

14 recordsLinked to original sources
USGS2010

Simulation of fluid distributions observed at a crude oil spill site incorporating hysteresis, oil entrapment, and spatial variability of hydraulic properties

Subsurface oil, water, and air saturation distributions were determined using 146 samples collected from seven boreholes along a 120-m transect at a crude oil spill site near Bemidji, Minnesota. The field data, collected 10 years after the spill, show a clearly defined oil body that has an oil saturation distribution that appears to be influenced by sediment heterogeneities and water table fluctuations. The center of the oil body has depressed the water-saturated zone boundary and the oil appears to have migrated laterally within the capillary fringe. A multiphase cross-sectional flow model was developed and used to simulate the movement of oil and water at the spill site. Comparisons between observed and simulated oil saturation distributions serve as an indicator of the appropriateness of using such models to predict the actual spread of organic immiscible liquids at spill sites. Sediment hydraulic properties used in the model were estimated from particle size data. The general large-scale features of the observed oil body were reproduced only when hysteresis with oil entrapment and representations of observed spatial variability of hydraulic properties were incorporated into the model. The small-scale details of the observed subsurface oil distribution were not reproduced in the simulations. The discrepancy between observed and simulated oil distributions reflects the considerable uncertainty in model parameter estimates and boundary conditions, three-phase capillary pressure-saturation-relative permeability functions, representations of spatial variability of hydraulic properties, and hydrodynamics of the groundwater flow system at the study site.

Water Resources Research↗
USGS2010

In situ measurements of volatile aromatic hydrocarbon biodegradation rates in groundwater

Benzene and alkylbenzene biodegradation rates and patterns were measured using an in situ microcosm in a crude-oil contaminated aquifer near Bemidji, Minnesota. Benzene-D6, toluene, ethylbenzene, o-, m- and p-xylenes and four pairs of C3- and C4-benzenes were added to an in situ microcosm and studied over a 3-year period. The microcosm allowed for a mass-balance approach and quantification of hydrocarbon biodegradation rates within a well-defined iron-reducing zone of the anoxic plume. Among the BTEX compounds, the apparent order of persistence is ethylbenzene > benzene > m,p-xylenes > o-xylene ≥ toluene. Threshold concentrations were observed for several compounds in the in situ microcosm, below which degradation was not observed, even after hundreds of days. In addition, long lag times were observed before the onset of degradation of benzene or ethylbenzene. The isomer-specific degradation patterns were compared to observations from a multi-year study conducted using data collected from monitoring wells along a flowpath in the contaminant plume. The data were fit with both first-order and Michaelis-Menten models. First-order kinetics provided a good fit for hydrocarbons with starting concentrations below 1 mg/L and Michaelis-Menten kinetics were a better fit when starting concentrations were above 1 mg/L, as was the case for benzene. The biodegradation rate data from this study were also compared to rates from other investigations reported in the literature.

Journal of Contaminant Hydrology↗
USGS2009

Benthic nutrient sources to hypereutrophic Upper Klamath Lake, Oregon, USA

Three collecting trips were coordinated in April, May, and August 2006 to sample the water column and benthos of hypereutrophic Upper Klamath Lake (OR, USA) through the annual cyanophyte bloom of Aphanizomenon flos‐aquae . A pore‐water profiler was designed and fabricated to obtain the first high‐resolution (centimeter‐scale) estimates of the vertical concentration gradients of macro‐ and micronutrients for diffusive‐flux determinations. A consistently positive benthic flux for soluble reactive phosphorus (SRP) was observed with solute release from the sediment, ranging between 0.4 and 6.1 mg/m 2 /d. The mass flux over an approximate 200‐km 2 lake area was comparable in magnitude to riverine inputs. An additional concern related to fish toxicity was identified when dissolved ammonium also displayed consistently positive benthic fluxes of 4 to 134 mg/m 2 /d, again comparable to riverine inputs. Although phosphorus was a logical initial choice by water quality managers for the limiting nutrient when nitrogen‐fixing cyanophytes dominate, initial trace‐element results from the lake and major inflowing tributaries suggested that the role of iron limitation on primary productivity should be investigated. Dissolved iron became depleted in the lake water column during the course of the algal bloom, while dissolved ammonium and SRP increased. Elevated macroinvertebrate densities, at least of the order of 10 4 individuals/m 2 , suggested that the diffusive‐flux estimates may be significantly enhanced by bioturbation. In addition, heat‐flux modeling indicated that groundwater advection of nutrients could also significantly contribute to internal nutrient loading. Accurate environmental assessments of lentic systems and reasonable expectations for point‐source management require quantitative consideration of internal solute sources.

Environmental Toxicology and Chemistry↗
USGS2008

Using heat to characterize streambed water flux variability in four stream reaches

Estimates of streambed water flux are needed for the interpretation of streambed chemistry and reactions. Continuous temperature and head monitoring in stream reaches within four agricultural watersheds (Leary Weber Ditch, IN; Maple Creek, NE; DR2 Drain, WA; and Merced River, CA) allowed heat to be used as a tracer to study the temporal and spatial variability of fluxes through the streambed. Synoptic methods (seepage meter and differential discharge measurements) were compared with estimates obtained by using heat as a tracer. Water flux was estimated by modeling one-dimensional vertical flow of water and heat using the model VS2DH. Flux was influenced by physical heterogeneity of the stream channel and temporal variability in stream and ground-water levels. During most of the study period (April–December 2004), flux was upward through the streambeds. At the IN, NE, and CA sites, high-stage events resulted in rapid reversal of flow direction inducing short-term surface-water flow into the streambed. During late summer at the IN site, regional ground-water levels dropped, leading to surface-water loss to ground water that resulted in drying of the ditch. Synoptic measurements of flux generally supported the model flux estimates. Water flow through the streambed was roughly an order of magnitude larger in the humid basins (IN and NE) than in the arid basins (WA and CA). Downward flux, in response to sudden high streamflows, and seasonal variability in flux was most pronounced in the humid basins and in high conductivity zones in the streambed.

Journal of Environmental Quality↗
USGS2005

Progression of methanogenic degradation of crude oil in the subsurface

Our results show that subsurface crude-oil degradation rates at a long-term research site were strongly influenced by small-scale variations in hydrologic conditions. The site is a shallow glacial outwash aquifer located near Bemidji in northern Minnesota that became contaminated when oil spilled from a broken pipeline in August 1979. In the study area, separate-phase oil forms a subsurface oil body extending from land surface to about 1 m (3.3 ft) below the 6–8-m (20–26 ft)-deep water table. Oil saturation in the sediments ranges from 10–20% in the vadose zone to 30–70% near the water table. At depths below 2 m (6.6 ft), degradation of the separate-phase crude oil occurs under methanogenic conditions. The sequence of methanogenic alkane degradation depletes the longer chain n-alkanes before the shorter chain n-alkanes, which is opposite to the better known aerobic sequence. The rates of degradation vary significantly with location in the subsurface. Oil-coated soils within 1.5 m (5 ft) of land surface have experienced little degradation where soil water saturation is less than 20%. Oil located 2–8 m (6.6–26 ft) below land surface in areas of higher recharge has been substantially degraded. The best explanation for the association between recharge and enhanced degradation seems to be increased downward transport of microbial growth nutrients to the oil body. This is supported by observations of greater microbial numbers at higher elevations in the oil body and significant decreases with depth in nutrient concentrations, especially phosphorus. Our results suggest that environmental effects may cause widely diverging degradation rates in the same spill, calling into question dating methods based on degradation state.

Environmental Geosciences↗
USGS2003

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m , p -xylene, and o -xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o -xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m , p -xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o -xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.

Minnesota↗
USGS2001

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↗
USGS2001

Progression of natural attenuation processes at a crude oil spill site: II. Controls on spatial distribution of microbial populations

A multidisciplinary study of a crude-oil contaminated aquifer shows that the distribution of microbial physiologic types is strongly controlled by the aquifer properties and crude oil location. The microbial populations of four physiologic types were analyzed together with permeability, pore-water chemistry, nonaqueous oil content, and extractable sediment iron. Microbial data from three vertical profiles through the anaerobic portion of the contaminated aquifer clearly show areas that have progressed from iron-reduction to methanogenesis. These locations contain lower numbers of iron reducers, and increased numbers of fermenters with detectable methanogens. Methanogenic conditions exist both in the area contaminated by nonaqueous oil and also below the oil where high hydrocarbon concentrations correspond to local increases in aquifer permeability. The results indicate that high contaminant flux either from local dissolution or by advective transport plays a key role in determining which areas first become methanogenic. Other factors besides flux that are important include the sediment Fe(II) content and proximity to the water table. In locations near a seasonally oscillating water table, methanogenic conditions exist only below the lowest typical water table elevation. During 20 years since the oil spill occurred, a laterally continuous methanogenic zone has developed along a narrow horizon extending from the source area to 50–60 m downgradient. A companion paper [J. Contam. Hydrol. 53, 369–386] documents how the growth of the methanogenic zone results in expansion of the aquifer volume contaminated with the highest concentrations of benzene, toluene, ethylbenzene, and xylenes.

Journal of Contaminant Hydrology↗
USGS1998

Ground water contamination by crude oil near Bemidji, Minnesota

Ground-water contamination by crude oil, and other petroleum-based liquids, is a widespread problem. An average of 83 crude-oil spills occurred per year during 1994-96 in the United States, each spilling about 50,000 barrels of crude oil (U.S. Office of Pipeline Safety, electronic commun., 1997). An understanding of the fate of organic contaminants (such as oil and gasoline) in the subsurface is needed to design innovative and cost-effective remedial solutions at contaminated sites.

Minnesota↗
USGS1992

Determination of subsurface fluid contents at a crude-oil spill site

Measurement of the fluid-content distribution at sites contaminated by immiscible fluids, including crude oil, is needed to better understand the movement of these fluids in the subsurface and to provide data to calibrate and verify numerical models and geophysical methods. A laboratory method was used to quantify the fluid contents of 146 core sections retrieved from boreholes aligned along a 120-m longitudinal transect at a crude-oil spill site near Bemidji, Minnesota, U.S.A. The 47-mm-diameter, minimally disturbed cores spanned a 4-m vertical interval contaminated by oil. Cores were frozen on site in a dry ice-alcohol bath to prevent redistribution and loss of fluids while sectioning the cores. We gravimetrically determined oil and water contents using a two-step method: (1) samples were slurried and the oil was removed by absorption onto strips of hydrophobic porous polyethylene (PPE); and (2) the samples were oven-dried to remove the water. The resulting data show sharp vertical gradients in the water and oil contents and a clearly defined oil body. The subsurface distribution is complex and appears to be influenced by sediment heterogeneities and water-table fluctuations. The center of the oil body has depressed the water-saturated zone boundary, and the oil is migrating laterally within the capillary fringe. The oil contents are as high as 0.3cm 3 cm −3 , which indicates that oil is probably still mobile 10 years after the spill occurred. The thickness of oil measured in wells suggests that accumulated thickness in wells is a poor indicator of the actual distribution of oil in the subsurface. Several possible sources of error are identified with the field and laboratory methods. An error analysis indicates that adsorption of water and sediment into the PPE adds as much as 4% to the measured oil masses and that uncertainties in the calculated sample volume and the assumed oil density introduce an additional ±3% error when the masses are converted to fluid contents.

Journal of Contaminant Hydrology↗
USGS1986

A comparison of the coupled fresh water-salt water flow and the Ghyben-Herzberg sharp interface approaches to modeling of transient behavior in coastal aquifer systems

A quasi-three dimensional finite difference model which simulates coupled, fresh water and salt water flow, separated by a sharp interface, is used to investigate the effects of storage characteristics, transmissivity, boundary conditions and anisotropy on the transient responses of such flow systems. The magnitude and duration of the departure of aquifer response from the behavior predicted using the Ghyben-Herzberg, one-fluid approach is a function of the ease with which flow can be induced in the salt water region. In many common hydrogeologic settings short-term fresh water head responses, and transitional responses between short-term and long-term, can only be realistically reproduced by including the effects of salt water flow on the dynamics of coastal flow systems. The coupled fresh water-salt water flow modeling approach is able to reproduce the observed annual fresh water head response of the Waialae aquifer of southeastern Oahu, Hawaii.

Hawaii↗