USGS Science⌕ Search

Geology topics

B.A. Bekins

Publications and source records attributed to B.A. Bekins.

At least 19 recordsLinked to original sources
USGS2010

Reanalysis of in situ permeability measurements in the Barbados décollement

A cased and sealed borehole in the Northern Barbados accretionary complex was the site of the first attempts to measure permeability in situ along a plate boundary décollement. Three separate efforts at Hole 949C yielded permeability estimates for the décollement spanning four orders of magnitude. An analysis of problems encountered during installation of the casing and seals provides insights into how the borehole conditions may have led to the wide range of results. During the installation, sediments from the surrounding formation repeatedly intruded into the borehole and casing. Stress analysis shows that the weak sediments were deforming plastically and the radial and tangential stresses around the borehole were significantly lower than lithostatic. This perturbed stress state may explain why the test pressure records showed indications of hydrofracture at pressures below lithostatic, and permeabilities rose rapidly as the estimated effective stress dropped below 0.8 MPa. Even after the borehole was sealed, the plastic deformation of the formation and relatively large gap of the wire wrapped screen allowed sediment to flow into the casing. Force equilibrium calculations predict sediment would have filled the borehole to 10 cm above the top of the screen by the time slug tests were conducted 1.5 years after the borehole was sealed. Reanalysis of the slug test results with these conditions yields several orders of magnitude higher permeability estimates than the original analysis which assumed an open casing. Overall the results based on only the tests with no sign of hydrofracture yield a permeability range of 10 −14 –10 −15 m 2 and a rate of increase in permeability with decreasing effective stress consistent with laboratory tests on samples from the décollement zone.

Geofluids↗
USGS2010

Milankovitch-scale correlations between deeply buried microbial populations and biogenic ooze lithology

The recent discoveries of large, active populations of microbes in the subseafloor of the world's oceans supports the impact of the deep biosphere biota on global biogeochemical cycles and raises important questions concerning the functioning of these extreme environments for life. These investigations demonstrated that subseafloor microbes are unevenly distributed and that cell abundances and metabolic activities are often independent from sediment depths, with increased prokaryotic activity at geochemical and/or sedimentary interfaces. In this study we demonstrate that microbial populations vary at the scale of individual beds in the biogenic oozes of a drill site in the eastern equatorial Pacific (Ocean Drilling Program Leg 201, Site 1226). We relate bedding-scale changes in biogenic ooze sediment composition to organic carbon (OC) and microbial cell concentrations using high-resolution color reflectance data as proxy for lithology. Our analyses demonstrate that microbial concentrations are an order of magnitude higher in the more organic-rich diatom oozes than in the nannofossil oozes. The variations mimic small-scale variations in diatom abundance and OC, indicating that the modern distribution of microbial biomass is ultimately controlled by Milankovitch-frequency variations in past oceanographic conditions. ?? 2010 Geological Society of America.

Geology↗
USGS2009

A critical evaluation of crustal dehydration as the cause of an overpressured and weak San Andreas Fault

Many plate boundary faults, including the San Andreas Fault, appear to slip at unexpectedly low shear stress. One long-standing explanation for a "weak" San Andreas Fault is that fluid release by dehydration reactions during regional metamorphism generates elevated fluid pressures that are localized within the fault, reducing the effective normal stress. We evaluate this hypothesis by calculating realistic fluid production rates for the San Andreas Fault system, and incorporating them into 2-D fluid flow models. Our results show that for a wide range of permeability distributions, fluid sources from crustal dehydration are too small and short-lived to generate, sustain, or localize fluid pressures in the fault sufficient to explain its apparent mechanical weakness. This suggests that alternative mechanisms, possibly acting locally within the fault zone, such as shear compaction or thermal pressurization, may be necessary to explain a weak San Andreas Fault. More generally, our results demonstrate the difficulty of localizing large fluid pressures generated by regional processes within near-vertical fault zones. ?? 2009 Elsevier B.V.

Earth and Planetary Science Letters↗
USGS2008

Limited occurrence of denitrification in four shallow aquifers in agricultural areas of the United States

The ability of natural attenuation to mitigate agricultural nitrate contamination in recharging aquifers was investigated in four important agricultural settings in the United States. The study used laboratory analyses, field measurements, and flow and transport modeling for monitoring well transects (0.5 to 2.5 km in length) in the San Joaquin watershed, California, the Elkhorn watershed, Nebraska, the Yakima watershed, Washington, and the Chester watershed, Maryland. Ground water analyses included major ion chemistry, dissolved gases, nitrogen and oxygen stable isotopes, and estimates of recharge date. Sediment analyses included potential electron donors and stable nitrogen and carbon isotopes. Within each site and among aquifer-based medians, dissolved oxygen decreases with ground water age, and excess N 2 from denitrification increases with age. Stable isotopes and excess N 2 imply minimal denitrifying activity at the Maryland and Washington sites, partial denitrification at the California site, and total denitrification across portions of the Nebraska site. At all sites, recharging electron donor concentrations are not sufficient to account for the losses of dissolved oxygen and nitrate, implying that relict, solid phase electron donors drive redox reactions. Zero-order rates of denitrification range from 0 to 0.14 μmol N L −1 d −1 , comparable to observations of other studies using the same methods. Many values reported in the literature are, however, orders of magnitude higher, which is attributed to a combination of method limitations and bias for selection of sites with rapid denitrification. In the shallow aquifers below these agricultural fields, denitrification is limited in extent and will require residence times of decades or longer to mitigate modern nitrate contamination.

Journal of Environmental Quality↗
USGS2008

Nitrogen fluxes through unsaturated zones in five agricultural settings across the United States

The main physical and chemical controls on nitrogen (N) fluxes between the root zone and the water table were determined for agricultural sites in California, Indiana, Maryland, Nebraska, and Washington from 2004 to 2005. Sites included irrigated and nonirrigated fields; soil textures ranging from clay to sand; crops including corn, soybeans, almonds, and pasture; and unsaturated zone thicknesses ranging from 1 to 22 m. Chemical analyses of water from lysimeters and shallow wells indicate that advective transport of nitrate is the dominant process affecting the flux of N below the root zone. Vertical profiles of (i) nitrogen species, (ii) stable isotopes of nitrogen and oxygen, and (iii) oxygen, N, and argon in unsaturated zone air and correlations between N and other agricultural chemicals indicate that reactions do not greatly affect N concentrations between the root zone and the capillary fringe. As a result, physical factors, such as N application rate, water inputs, and evapotranspiration, control the differences in concentrations among the sites. Concentrations of N in shallow lysimeters exhibit seasonal variation, whereas concentrations in lysimeters deeper than a few meters are relatively stable. Based on concentration and recharge estimates, fluxes of N through the deep unsaturated zone range from 7 to 99 kg ha −1 yr −1 Vertical fluxes of N in ground water are lower due to spatial and historical changes in N inputs. High N fluxes are associated with coarse sediments and high N application rates.

Journal of Environmental Quality↗
USGS2007

Dissolution of biogenic ooze over basement edifices in the equatorial Pacific with implications for hydrothermal ventilation of the oceanic crust

Recent observations indicate that curious closed depressions in carbonate sediments overlying basement edifices are widespread in the equatorial Pacific. A possible mechanism for their creation is dissolution by fluids exiting basement vents from off-axis hydrothermal flow. Quantitative analysis based on the retrograde solubility of calcium carbonate and cooling of basement fluids during ascent provides an estimate for the dissolution capacity of the venting fluids. Comparison of the dissolution capacity and fluid flux with typical equatorial Pacific carbonate mass accumulation rates shows that this mechanism is feasible. By maintaining sediment-free basement outcrops, the process may promote widespread circulation of relatively unaltered seawater in the basement in an area where average sediment thicknesses are 300–500 m. The enhanced ventilation can explain several previously puzzling observations in this region, including anomalously low heat flux, relatively unaltered seawater in the basement, and aerobic and nitrate-reducing microbial activity at the base of the sediments.

Geology↗
USGS2005

A simple method for calculating growth rates of petroleum hydrocarbon plumes

Consumption of aquifer Fe(III) during biodegradation of ground water contaminants may result in expansion of a contaminant plume, changing the outlook for monitored natural attenuation. Data from two research sites contaminated with petroleum hydrocarbons show that toluene and xylenes degrade under methanogenic conditions, but the benzene and ethylbenzene plumes grow as aquifer Fe(III) supplies are depleted. By considering a one‐dimensional reaction front in a constant unidirectional flow field, it is possible to derive a simple expression for the growth rate of a benzene plume. The method balances the mass flux of benzene with the Fe(III) content of the aquifer, assuming that the biodegradation reaction is instantaneous. The resulting expression shows that the benzene front migration is retarded relative to the ground water velocity by a factor that depends on the concentrations of hydrocarbon and bioavailable Fe(III). The method provides good agreement with benzene plumes at a crude oil study site in Minnesota and a gasoline site in South Carolina. Compared to the South Carolina site, the Minnesota site has 25% higher benzene flux but eight times the Fe(III), leading to about one‐sixth the expansion rate. Although it was developed for benzene, toluene, ethylbenzene, and xylenes, the growth‐rate estimation method may have applications to contaminant plumes from other persistent contaminant sources.

Ground Water↗
USGS2005

A comparison of zero-order, first-order, and monod biotransformation models

Under some conditions, a first-order kinetic model is a poor representation of biodegradation in contaminated aquifers. Although it is well known that the assumption of first-order kinetics is valid only when substrate concentration, S, is much less than the half-saturation constant, K(s), this assumption is often made without verification of this condition. We present a formal error analysis showing that the relative error in the first-order approximation is S/K(S) and in the zero-order approximation the error is K(s)/S. We then examine the problems that arise when the first-order approximation is used outside the range for which it is valid. A series of numerical simulations comparing results of first- and zero-order rate approximations to Monod kinetics for a real data set illustrates that if concentrations observed in the field are higher than K(s), it may better to model degradation using a zero-order rate expression. Compared with Monod kinetics, extrapolation of a first-order rate to lower concentrations under-predicts the biotransformation potential, while extrapolation to higher concentrations may grossly over-predict the transformation rate. A summary of solubilities and Monod parameters for aerobic benzene, toluene, and xylene (BTX) degradation shows that the a priori assumption of first-order degradation kinetics at sites contaminated with these compounds is not valid. In particular, out of six published values of KS for toluene, only one is greater than 2 mg/L, indicating that when toluene is present in concentrations greater than about a part per million, the assumption of first-order kinetics may be invalid. Finally, we apply an existing analytical solution for steady-state one-dimensional advective transport with Monod degradation kinetics to a field data set.A formal error analysis is presented showing that the relative error in the first-order approximation is S/KS and in the zero-order approximation the error is KS/S where S is the substrate concentration and KS is the half-saturation constant. The problems that arise when the first-order approximation is used outside the range for which it is valid are examined. A series of numerical simulations comparing results of first- and zero-order rate approximations to Monod kinetics for a real data set illustrates that if concentrations observed in the field are higher than KS, it may be better to model degradation using a zero-order rate expression.

Ground Water↗
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↗
USGS2003

Topographically driven groundwater flow and the San Andreas heat flow paradox revisited

Evidence for a weak San Andreas Fault includes (1) borehole heat flow measurements that show no evidence for a frictionally generated heat flow anomaly and (2) the inferred orientation of ??1 nearly perpendicular to the fault trace. Interpretations of the stress orientation data remain controversial, at least in close proximity to the fault, leading some researchers to hypothesize that the San Andreas Fault is, in fact, strong and that its thermal signature may be removed or redistributed by topographically driven groundwater flow in areas of rugged topography, such as typify the San Andreas Fault system. To evaluate this scenario, we use a steady state, two-dimensional model of coupled heat and fluid flow within cross sections oriented perpendicular to the fault and to the primary regional topography. Our results show that existing heat flow data near Parkfield, California, do not readily discriminate between the expected thermal signature of a strong fault and that of a weak fault. In contrast, for a wide range of groundwater flow scenarios in the Mojave Desert, models that include frictional heat generation along a strong fault are inconsistent with existing heat flow data, suggesting that the San Andreas Fault at this location is indeed weak. In both areas, comparison of modeling results and heat flow data suggest that advective redistribution of heat is minimal. The robust results for the Mojave region demonstrate that topographically driven groundwater flow, at least in two dimensions, is inadequate to obscure the frictionally generated heat flow anomaly from a strong fault. However, our results do not preclude the possibility of transient advective heat transport associated with earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Microbial ecology of a crude oil contaminated aquifer

Detailed microbial analyses of a glacial outwash aquifer contaminated by crude oil provide insights into the pattern of microbial succession from iron reducing to methanogenic in the anaerobic portion of the contaminant plume. We analysed sediments from this area for populations of aerobes, iron reducers, fermenters and methanogens, using the most probable number method. On the basis of the microbial data the anaerobic area can be divided into distinct physiological zones dominated by either iron-reducers or a consortium of fermenters and methanogens. Chemistry and permeability data show that methanogenic conditions develop first in areas of high hydrocarbon flux. Thus, we find methanogens both in high permeability horizons and also where separate-phase crude oil is present in either the saturated or unsaturated zone. Microbial numbers peak at the top of the separate-phase oil suggesting that growth is most rapid in locations with access to both hydrocarbons and nutrients infiltrating from the surface.

IAHS-AISH Publication↗
USGS2001

Modeling consolidation and dewatering near the toe of the northern Barbados accretionary complex

At the toe of the northern Barbados accretionary complex, temperature and pore water chemistry data indicate that fluid flow is channeled along the décollement and other shallow thrust faults. We examine mechanisms that may prevent consolidation and maintain high permeability over large sections of the décollement. High-resolution bulk density data from five boreholes show that the décollement is well consolidated at some sites while other sites remain underconsolidated. Underconsolidated décollement behavior is associated with kilometer-scale negative-polarity seismic reflections from the décollement plane that have been interpreted to be fluid conduits. We use a coupled fluid flow/consolidation model to simulate the loading response of a 10-km-long by 680-m-thick slice of sediment as it enters the accretionary complex. The simulations capture 185 ka (5 km) of subduction, with a load function representing the estimated effective stress of the overriding accretionary prism (3.8° taper angle). Simulation results of bulk density in the décollement 3.2 km arcward of the deformation front are compared with observations. The results show that persistent high pore pressures at the arcward edge of the simulation domain can explain underconsolidated behavior. The scenario is consistent with previous modeling results showing that high pore pressures can propagate intermittently along the décollement from deeper in the complex. Simulated seaward fluxes in the décollement (1–14 cm yr −1 ) lie between previous estimates from modeling studies of steady state (<1 cm yr −1 ) and transient (>1 m yr −1 ) flow. Maximum simulated instantaneous fluid sources (2.5×10 −13 s −1 ) are comparable to previous estimates. The simulations show minor swelling of incoming sediments (fluid sources ∼−3×10 −15 s −1 ) up to 3 km before subduction that may help to explain small-scale shearing and normal faulting proximal to the protodécollement.

Atlantic Ocean, Northern Barbados Accretionary Com↗
USGS2001

Smectite diagenesis, pore-water freshening, and fluid flow at the toe of the Nankai wedge

The presence of low-chloride fluids in the lowermost sediments drilled at Ocean Drilling Program Site 808, at the Nankai accretionary wedge, has been considered as prime evidence for long-distance, lateral fluid flow from depth. Here, we re-evaluate the potential role of in situ reaction of smectite (S) to illite (I) in the genesis of this low chloride anomaly. This reaction is known to be occurring at Site 808, with both the S content and S to I ratio in the mixed layer clays decreasing substantially with depth. We show that the bulk of the chloride anomaly can generate by in situ clay dehydration, particularly if pre-reaction smectite abundances (Ai) approach ?? 10-15% of the bulk sediment. The Ai values, however, are not well constrained. At Ai values < 10-15%, an additional source of low-Cl fluid centered close to the de??collement could be required. Thus, there remains the important possibility that the observed low-Cl anomaly is a compound effect of both lateral flow and in situ smectite dehydration. ?? 2001 Elsevier Science B.V. All rights reserved.

Earth and Planetary Science Letters↗
USGS2001

Progression of natural attenuation processes at a crude-oil spill site . I. Geochemical evolution of the plume

A 16-year study of a hydrocarbon plume shows that the extent of contaminant migration and compound-specific behavior have changed as redox reactions, most notably iron reduction, have progressed over time. Concentration changes at a small scale, determined from analysis of pore-water samples drained from aquifer cores, are compared with concentration changes at the plume scale, determined from analysis of water samples from an observation well network. The small-scale data show clearly that the hydrocarbon plume is growing slowly as sediment iron oxides are depleted. Contaminants, such as ortho -xylene that appeared not to be moving downgradient from the oil on the basis of observation well data, are migrating in thin layers as the aquifer evolves to methanogenic conditions. However, the plume-scale observation well data show that the downgradient extent of the Fe 2+ and BTEX plume did not change between 1992 and 1995. Instead, depletion of the unstable Fe (III) oxides near the subsurface crude-oil source has caused the maximum dissolved iron concentration zone within the plume to spread at a rate of approximately 3 m/year. The zone of maximum concentrations of benzene, toluene, ethylbenzene and xylene (BTEX) has also spread within the anoxic plume. In monitoring the remediation of hydrocarbon-contaminated ground water by natural attenuation, subtle concentration changes in observation well data from the anoxic zone may be diagnostic of depletion of the intrinsic electron-accepting capacity of the aquifer. Recognition of these subtle patterns may allow early prediction of growth of the hydrocarbon plume.

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

Microbial populations in contaminant plumes

Efficient biodegradation of subsurface contaminants requires two elements: (1) microbial populations with the necessary degradative capabilities, and (2) favorable subsurface geochemical and hydrological conditions. Practical constraints on experimental design and interpretation in both the hydrogeological and microbiological sciences have resulted in limited knowledge of the interaction between hydrogeological and microbiological features of subsurface environments. These practical constraints include: (1) inconsistencies between the scales of investigation in the hydrogeological and microbiological sciences, and (2) practical limitations on the ability to accurately define microbial populations in environmental samples. However, advances in application of small-scale sampling methods and interdisciplinary approaches to site investigations are beginning to significantly improve understanding of hydrogeological and microbiological interactions. Likewise, culture-based and molecular analyses of microbial populations in subsurface contaminant plumes have revealed significant adaptation of microbial populations to plume environmental conditions. Results of recent studies suggest that variability in subsurface geochemical and hydrological conditions significantly influences subsurface microbial-community structure. Combined investigations of site conditions and microbial-community structure provide the knowledge needed to understand interactions between subsurface microbial populations, plume geochemistry, and contaminant biodegradation.

Hydrogeology Journal↗
USGS1999

Fluid budgets at convergent plate margins: Implications for the extent and duration of fault-zone dilation

Faults at convergent plate boundaries are important conduits for fluid escape, and recent evidence suggests that fluid expulsion along them is both transient and heterogeneous. For the Nankai and Barbados convergent margins, we have used numerical models to investigate the long-term partitioning of expelled fluids between diffuse flow and flow along connected high-permeability fault conduits. For a simple case of spatial heterogeneity, we estimated the extent of high-permeability conduits necessary to maintain a balance between incoming and expelled fluids. For the case of transient expulsion, we constrained the duration of elevated permeability required to balance the fluid budgets. Comparison of modeled and observed geochemical profiles suggests that the initiation of connected flow conduits is delayed with respect to the time of accretion into each accretionary complex and may be related to burial below a critical depth, either where the overlying wedge is sufficiently thick to prevent fluid escape to the sea floor or where sediments behave brittlely.

Geology↗