Search USGSSearch

SEARCH · Search USGS

Results for “Physics of Fluids”

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 217 records · Page 12Linked to original sources

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Gas hydrate systems along the eastern continental margin of India

The primary objectives of the India National Gas Hydrate Program Expedition 02 (NGHP-02) were to obtain new data on the occurrence of gas hydrate systems and to advance the understanding of the controls on the formation of gas hydrate accumulations in the offshore of India. In accordance with the ultimate overall goal of the NGHP effort to assess the energy resource potential of marine gas hydrates in India, particular focus was placed on the exploration and evaluation of gas hydrate occurrences at high saturations in sand-rich systems. NGHP-02 operations were conducted from 3-March-2015 to 28-July-2015 off the eastern coast of India and included logging while drilling (LWD) operations at 25 locations, and coring and wireline logging operations at 10 locations, in the Krishna-Godavari and Mahanadi Basins. The formation of highly concentrated gas hydrate accumulations, which are more suitable for energy extraction, requires the presence of relatively coarse-grained sediments with porosity needed to support the migration and accumulation of gas, and the nucleation of gas hydrate. The results of downhole logging, coring and formation pressure testing operations during NGHP-02 have confirmed the presence of extensive sand-rich depositional systems throughout the deepwater portions of the Krishna-Godavari and Mahanadi Basins. Two areas of Krishna-Godavari Basin, referred to as Areas B and C, contain substantial gas hydrate accumulations in sand-rich systems and therefore represent ideal candidate sites for future gas hydrate production testing. This summary and technical report includes a comprehensive synthesis of the geologic, geophysical, geochemical, and physical property data acquired during NGHP-02 as it relates to the controls on gas hydrate occurrence, particularly with regards to sand-hosted accumulations. In the Mahanadi Basin, despite the confirmation of extensive reservoir capacity, gas supply at the NGHP-02 sites was insufficient to charge the reservoirs with gas hydrates. In the Krishna-Godavari Basin, extensive reservoir systems were confirmed with sediment grain-sizes ranging from coarse-silts to gravels. These reservoirs range from fully- to partially filled with gas hydrate. The gas is determined to be from only microbial sources, and in part migrated into the reservoirs from deeper systems. The controls on gas hydrate occurrence are complex and varied; and include substantial reservoir heterogeneity and sufficient permeability throughout the reservoirs and seals that allowed pervasive fluid flow into and through the hydrate-bearing systems. These discoveries are the most significant confirmation of the exploration approach that focuses on direct detection of hydrate reservoirs supported by comprehensive petroleum systems analyses.

Journal of Marine and Petroleum Geology

Trace element abundances in major minerals of Late Permian coals from southwestern Guizhou province, China

Fourteen samples of minerals were separated by handpicking from Late Permian coals in southwestern Guizhou province, China. These 14 minerals were nodular pyrite, massive recrystallized pyrite, pyrite deposited from low-temperature hydrothermal fluid and from ground water; clay minerals; and calcite deposited from low-temperature hydrothermal fluid and from ground water. The mineralogy, elemental composition, and distribution of 33 elements in these samples were studied by optical microscopy, scanning electron microscope equipped with energy-dispersive X-ray spectrometer (SEM-EDX), X-ray diffraction (XRD), cold-vapor atomic absorption spectrometry (CV-AAS), atomic fluorescence spectrometry (AFS), inductively coupled-plasma mass spectrometry (ICP-MS), and ion-selective electrode (ISE). The results show that various minerals in coal contain variable amounts of trace elements. Clay minerals have high concentrations of Ba, Be, Cs, F, Ga, Nb, Rb, Th, U, and Zr. Quartz has little contribution to the concentration of trace elements in bulk coal. Arsenic, Mn, and Sr are in high concentrations in calcite. Pyrite has high concentrations of As, Cd, Hg, Mo, Sb, Se, Tl, and Zn. Different genetic types of calcite in coal can accumulate different trace elements; for example Ba, Co, Cr, Hg, Ni, Rb, Sn, Sr, and Zn are in higher concentrations in calcite deposited from low-temperature hydrothermal fluid than in that deposited from ground water. Furthermore, the concentrations of some trace elements are quite variable in pyrite; different genetic types of pyrites (Py-A, B, C, D) have different concentrations of trace elements, and the concentrations of trace elements are also different in pyrite of low-temperature hydrothermal origin collected from different locations. The study shows that elemental concentration is rather uniform in a pyrite vein. There are many micron and submicron mosaic pyrites in a pyrite vein, which is enriched in some trace elements, such as As and Mo. The content of trace element in pyrite vein depends upon the content of mosaic pyrite and of trace elements in it. Many environmentally sensitive trace elements are mainly contained in the minerals in coal, and hence the physical coal cleaning techniques can remove minerals from coal and decrease the emissions of potentially hazardous trace elements. ?? 2002 Elsevier Science B.V. All rights reserved.

International Journal of Coal Geology

Thermal Methods for Investigating Ground-Water Recharge

Recharge of aquifers within arid and semiarid environments is defined as the downward flux of water across the regional water table. The introduction of recharging water at the land surface can occur at discreet locations, such as in stream channels, or be distributed over the landscape, such as across broad interarroyo areas within an alluvial ground-water basin. The occurrence of recharge at discreet locations is referred to as focused recharge, whereas the occurrence of recharge over broad regions is referred to as diffuse recharge. The primary interest of this appendix is focused recharge, but regardless of the type of recharge, estimation of downward fluxes is essential to its quantification. Like chemical tracers, heat can come from natural sources or be intentionally introduced to infer transport properties and aquifer recharge. The admission and redistribution of heat from natural processes such as insolation, infiltration, and geothermal activity can be used to quantify subsurface flow regimes. Heat is well suited as a ground-water tracer because it provides a naturally present dynamic signal and is relatively harmless over a useful range of induced perturbations. Thermal methods have proven valuable for recharge investigations for several reasons. First, theoretical descriptions of coupled water-and-heat transport are available for the hydrologic processes most often encountered in practice. These include land-surface mechanisms such as radiant heating from the sun, radiant cooling into space, and evapotranspiration, in addition to the advective and conductive mechanisms that usually dominate at depth. Second, temperature is theoretically well defined and readily measured. Third, thermal methods for depths ranging from the land surface to depths of hundreds of meters are based on similar physical principles. Fourth, numerical codes for simulating heat and water transport have become increasingly reliable and widely available. Direct measurement of water flux in the subsurface is difficult, prompting investigators to pursue indirect methods. Geophysical approaches that exploit the coupled relation between heat and water transport provide an attractive class of methods that have become widely used in investigations of recharge. This appendix reviews the application of heat to the problem of recharge estimation. Its objective is to provide a fairly complete account of the theoretical underpinnings together with a comprehensive review of thermal methods in practice. Investigators began using subsurface temperatures to delineate recharge areas and infer directions of ground-water flow around the turn of the 20th century. During the 1960s, analytical and numerical solutions for simplified heat- and fluid-flow problems became available. These early solutions, though one-dimensional and otherwise restricted, provided a strong impetus for applying thermal methods to problems of liquid and vapor movement in systems ranging from soils to geothermal reservoirs. Today?s combination of fast processors, massive data-storage units, and efficient matrix techniques provide numerical solutions to complex, three-dimensional transport problems. These approaches allow researchers to take advantage of the considerable information content routinely achievable in high-accuracy temperature work.

Professional Paper

Analysis of dispersion and attenuation of surface waves in poroelastic media in the exploration-seismic frequency band

We analyse dispersion and attenuation of surface waves at free surfaces of possible vacuum/poroelastic media: permeable-‘open pore’, impermeable-‘closed pore’ and partially permeable boundaries, which have not been previously reported in detail by researchers, under different surface-permeable, viscous-damping, elastic and fluid-flowing conditions. Our discussion is focused on their characteristics in the exploration-seismic frequency band (a few through 200 Hz) for near-surface applications. We find two surface-wave modes exist, R 1 waves for all conditions, and R 2 waves for closed-pore and partially permeable conditions. For R 1 waves, velocities disperse most under partially permeable conditions and least under the open-pore condition. High-coupling damping coefficients move the main dispersion frequency range to high frequencies. There is an f 1 frequency dependence as a constant- Q model for attenuation at high frequencies. R 1 waves for the open pore are most sensitive to elastic modulus variation, but least sensitive to tortuosities variation. R 1 waves for partially permeable surface radiate as non-physical waves (Im( k ) < 0) at low frequencies. For R 2 waves, velocities are slightly lower than the bulk slow P 2 waves. At low frequencies, both velocity and attenuation are diffusive of f 1/2 frequency dependence, as P 2 waves. It is found that for partially permeable surfaces, the attenuation displays - f 1 frequency dependence as frequency increasing. High surface permeability, low-coupling damping coefficients, low Poisson′s ratios, and low tortuosities increase the slope of the - f 1 dependence. When the attenuation coefficients reach 0, R 2 waves for partially permeable surface begin to radiate as non-physical waves.

Geophysical Journal International

SPEAR: The next generation GFDL modeling system for seasonal to multidecadal prediction and projection

We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR ( S eamless System for P rediction and EA rth System R esearch) is built from component models recently developed at GFDL—the AM4 atmosphere model, MOM6 ocean code, LM4 land model, and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0° (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1° to 0.25°. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM4 models but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM4 models.

Journal of Advances in Modeling Earth Systems

Review of magnetic field monitoring near active faults and volcanic calderas in California: 1974-1995

Differential magnetic fields have been monitored along the San Andreas fault and the Long Valley caldera since 1974. At each monitoring location, proton precession magnetometers sample total magnetic field intensity at a resolution of 0.1 nT or 0.25 nT. Every 10 min, data samples are transmitted via satellite telemetry to Menlo Park, CA for processing and analysis. The number of active magnetometer sites has varied during the past 21 years from 6 to 25, with 12 sites currently operational. We use this network to identify magnetic field changes generated by earthquake and volcanic processes. During the two decades of monitoring, five moderate earthquakes (M5.9 to M7.3) have occurred within 20 km of magnetometer sites located along the San Andreas fault and only one preseismic signal of 1.5 nT has been observed. During moderate earthquakes, coseismic magnetic signals, with amplitudes from 0.7 nT to 1.3 nT, have been identified for 3 of the 5 events. These observations are generally consistent with those calculated from simple seismomagnetic models of these earthquakes and near-fault coseismic magnetic field disturbances rarely exceed one nanotesla. These data are consistent with the concept of low shear stress and relatively uniform displacement of the San Andreas fault system as expected due to high pore fluid pressure on the fault. A systematic decrease of 0.8-1 nT/year in magnetic field has occurred in the Long Valley caldera since 1989. These magnetic field data are similar in form to observed geodetically measured displacements from inflation of the resurgent dome. A simple volcanomagnetic model involving pressure increase of 50 MPa/a at a depth of 7 km under the resurgent dome can replicate these magnetic field observations. This model is derived from the intrusion model that best fits the surface deformation data. ?? 1998 Elsevier Science B.V.

Physics of the Earth and Planetary Interiors

Site selection for DOE/JIP gas hydrate drilling in the northern Gulf of Mexico

In the late spring of 2008, the Chevron-led Gulf of Mexico Gas Hydrate Joint Industry Project (JIP) expects to conduct an exploratory drilling and logging campaign to better understand gas hydrate-bearing sands in the deepwater Gulf of Mexico. The JIP Site Selection team selected three areas to test alternative geological models and geophysical interpretations supporting the existence of potential high gas hydrate saturations in reservoir-quality sands. The three sites are near existing drill holes which provide geological and geophysical constraints in Alaminos Canyon (AC) lease block 818, Green Canyon (GC) 955, and Walker Ridge (WR) 313. At the AC818 site, gas hydrate is interpreted to occur within the Oligocene Frio volcaniclastic sand at the crest of a fold that is shallow enough to be in the hydrate stability zone. Drilling at GC955 will sample a faulted, buried Pleistocene channel-levee system in an area characterized by seafloor fluid expulsion features, structural closure associated with uplifted salt, and abundant seismic evidence for upward migration of fluids and gas into the sand-rich parts of the sedimentary section. Drilling at WR313 targets ponded sheet sands and associated channel/levee deposits within a minibasin, making this a non-structural play. The potential for gas hydrate occurrence at WR313 is supported by shingled phase reversals consistent with the transition from gas-charged sand to overlying gas-hydrate saturated sand. Drilling locations have been selected at each site to 1) test geological methods and models used to infer the occurrence of gas hydrate in sand reservoirs in different settings in the northern Gulf of Mexico; 2) calibrate geophysical models used to detect gas hydrate sands, map reservoir thicknesses, and estimate the degree of gas hydrate saturation; and 3) delineate potential locations for subsequent JIP drilling and coring operations that will collect samples for comprehensive physical property, geochemical and other analyses

Gulf of Mexico

Inferring shallow groundwater flow in saprolite and fractured rock using environmental tracers

The Ridge and Valley Province of eastern Tennessee is characterized by (1) substantial topographic relief, (2) folded and highly fractured rocks of various lithologies that have low primary permeability and porosity, and (3) a shallow residuum of medium permeability and high total porosity. Conceptual models of shallow groundwater flow and solute transport in this system have been developed but are difficult to evaluate using physical characterization or short‐term tracer methods due to extreme spatial variability in hydraulic properties. In this paper we describe how chlorofluorocarbon 12, 3 H, and 3 He were used to infer groundwater flow and solute transport in saprolite and fractured rock near Oak Ridge, Tennessee. In the shallow residuum, fracture spacings are <0.05 m, suggesting that concentrations of these tracers in fractures and in the matrix have time to diffusionally equilibrate. The relatively smooth nature of tracer concentrations with depth in the residuum is consistent with this model and quantitatively suggests recharge fluxes of 0.2 to 0.4 m yr −1 . In contrast, groundwater flow within the unweathered rock appears to be controlled by fractures with spacings of the order of 2 to 5 m, and diffusional equilibration of fractures and matrix has not occurred. For this reason, vertical fluid fluxes in the unweathered rock cannot be estimated from the tracer data.

Water Resources Research

Assessing factors affecting the thermal properties of a passive thermal refuge using three-dimensional hydrodynamic flow and transport modeling

Everglades restoration activities may cause changes to temperature and salinity stratification at the Port of the Islands (POI) marina, which could affect its suitability as a cold weather refuge for manatees. To better understand how the Picayune Strand Restoration Project (PSRP) may alter this important resource in Collier County in southwestern Florida, the USGS has developed a three-dimensional hydrodynamic model for the marina and canal system at POI. Empirical data suggest that manatees aggregate at the site during winter because of thermal inversions that provide warmer water near the bottom that appears to only occur in the presence of salinity stratification. To study these phenomena, the environmental fluid dynamics code simulator was used to represent temperature and salinity transport within POI. Boundary inputs were generated using a larger two-dimensional model constructed with the flow and transport in a linked overland-aquifer density-dependent system simulator. Model results for a representative winter period match observed trends in salinity and temperature fluctuations and produce temperature inversions similar to observed values. Modified boundary conditions, representing proposed PSRP alterations, were also tested to examine the possible effect on the salinity stratification and temperature inversion within POI. Results show that during some periods, salinity stratification is reduced resulting in a subsequent reduction in temperature inversion compared with the existing conditions simulation. This may have an effect on POI’s suitability as a passive thermal refuge for manatees and other temperature-sensitive species. Additional testing was completed to determine the important physical relationships affecting POI’s suitability as a refuge.

Florida

Parametric study of the physical properties of hydrate-bearing sand, silt, and clay sediments: 1. Electromagnetic properties

The marked decrease in bulk electrical conductivity of sediments in the presence of gas hydrates has been used to interpret borehole electrical resistivity logs and, to a lesser extent, the results of controlled source electromagnetic surveys to constrain the spatial distribution and predicted concentration of gas hydrate in natural settings. Until now, an exhaustive laboratory data set that could be used to assess the impact of gas hydrate on the electromagnetic properties of different soils (sand, silt, and clay) at different effective stress and with different saturations of hydrate has been lacking. The laboratory results reported here are obtained using a standard geotechnical cell and the hydrate-formed tetrahydrofuran (THF), a liquid that is fully miscible in water and able to produce closely controlled saturations of hydrate from dissolved phase. Both permittivity and electrical conductivity are good indicators of the volume fraction of free water in the sediment, which is in turn dependent on hydrate saturation. Permittivity in the microwave frequency range is particularly predictive of free water content since it is barely affected by ionic concentration, pore structure, and surface conduction. Electrical conductivity (or resistivity) is less reliable for constraining water content or hydrate saturation: In addition to fluid-filled porosity, other factors, such as the ionic concentration of the pore fluid and possibly other conduction effects (e.g., surface conduction in high specific surface soils having low conductivity pore fluid), also influence electrical conductivity.

Journal of Geophysical Research

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface

A reassessment of the role of tidal dispersion in estuaries and bays

The role of tidal dispersion is reassessed, based on a consideration of the relevant physical mechanisms, particularly those elucidated by numerical simulations of tide-induced dispersion. It appears that the principal influence of tidal currents on dispersion occurs at length scales of the tidal excursion and smaller; thus the effectiveness of tidal dispersion depends on the relative scale of the tidal excursion to the spacing between major bathymetric and shoreline features. In estuaries where the typical spacing of topographic features is less than the tidal excursion, tidal dispersion may contribute significantly to the overall flushing. In estuaries and embayments in which the typical spacing between major features is larger than the tidal excursion, the influence of tidal dispersion will be localized, and it will not markedly contribute to overall flushing. Tidal dispersion is most pronounced in regions of abrupt topographic changes such as headlands and inlets, where flow separation occurs. The strong strain rate in the region of flow separation tends to stretch patches of fluid into long filaments, which are subsequently rolled up and distorted by the transient eddy field. The dispersion process accomplished by the tides varies strongly as a function of position and tidal phase and thus does not lend itself to parameterization by an eddy diffusion coefficient.

Estuaries

Chemical hydrogeology in natural and contaminated environments

Chemical hydrogeology, including organic and inorganic aspects, has contributed to an increased understanding of groundwater flow systems, geologic processes, and stressed environments. Most of the basic principles of inorganic-chemical hydrogeology were first established by investigations of organic-free, regional-scale systems for which simplifying assumptions could be made. The problems of groundwater contamination are causing a shift of emphasis to microscale systems that are dominated by organic-chemical reactions and that are providing an impetus for the study of naturally occurring and manmade organic material. Along with the decrease in scale, physical and chemical heterogeneity become major controls. Current investigations and those selected from the literature demonstrate that heterogeneity increases in importance as the study site decreases from regional-scale to macroscale to microscale. Increased understanding of regional-scale flow systems is demonstrated by selection of investigations of carbonate and volcanic aquifers to show how applications of present-day concepts and techniques can identify controlling chemical reactions and determine their rates; identify groundwater flow paths and determine flow velocity; and determine aquifer characteristics. The role of chemical hydrogeology in understanding geologic processes of macroscale systems is exemplified by selection of investigations in coastal aquifers. Phenomena associated with the mixing zone generated by encroaching sea water include an increase in heterogeneity of permeability, diagenesis of minerals, and formation of geomorphic features, such as caves, lagoons, and bays. Ore deposits of manganese and uranium, along with a simulation model of ore-forming fluids, demonstrate the influence of heterogeneity and of organic compounds on geochemical reactions associated with genesis of mineral deposits. In microscale environments, importance of heterogeneity and consequences of organic reactions in determining the distribution and concentrations cf. constituents are provided by several studies, including infiltration of sewage effluent and migration of creosote in coastal plain aquifers. These studies show that heterogeneity and the dominance of organically controlled reactions greatly increase the complexity of investigations.Current investigations and those selected from the literature demonstrate that heterogeneity increases in importance as the study site decreases from regional-scale to macroscale to microscale. Increased understanding of regional-scale flow systems is demonstrated by selection of investigations of carbonate and volcanic aquifers to show how application of present-day concepts and techniques can identify controlling chemical reactions and determine their rates; identify groundwater flow paths and determine flow velocity; and determine aquifer characteristics. The role of chemical hydrogeology in understanding geologic processes of macroscale systems is exemplified by selection of investigations in coastal aquifers. Ore deposits of manganese and uranium, along with a simulation model of ore-forming fluids, demonstrate the influence of heterogeneity and of organic compounds on geochemical reactions associated with genesis of mineral deposits. In microscale environments, importance of heterogeneity is illustrated by studies of infiltration of sewage effluent and migration of creosote in coastal plain aquifers.

Journal of Hydrology

Reconnaissance geotechnical study of the Lucerne Granite, Maine

Physical properties, including density and strength, of the Lucerne Granite (Wing, 1958) vary because of inhomogeneities caused by fracturing and differences in texture. Two granitic facies were found: a coarse-grained variety in the southern part of the body and near the edges of the pluton and a medium-grained porphyritic variety in the central and northern sections. Point-load strength tests indicate that the medium-grained porphyritic granite is stronger than the coarse-grained porphyritic or non-porphyritic granite. Results of bulk and grain density tests suggest that rock strength in the weaker granite samples is controlled by microfractures. In the pluton, joint trends changed from a scattered pattern in the north to a preferred orientation in the south. Because it has a simple joint pattern and is homogeneous in texture, the southern part of the Lucerne Granite is a more suitable area for excavation and fluid storage than the northern section.

Maine

Making fired bricks with spent equilibrium catalyst-a technical feasibility study

Fluid catalytic cracking in an oil refinery uses a catalyst, such as an alumino-silicate zeolite, in the conversion of heavy hydrocarbons to light hydrocarbons. A small fraction of the catalyst is continually replaced with fresh catalyst to maintain activity. In North America, more than 400 tons of spent alumino-silicate equilibrium catalyst (spent e-cat), and worldwide, more than 1,100 tons, are generated daily, most of which is disposed of in landfills (municipal and on-site facilities). In this study, three spent e-cat samples were tested in a value-added application that would utilize this waste in the manufacturing of fired bricks. The results of this study indicate that spent e-cat is a technically feasible raw material substitute for the clay and shale commonly used in fired brick production. Fired bricks produced with up to 30 wt% of spent e-cat showed good physical appearance and their water absorption properties met the ASTM C 62 specifications for building bricks of either the moderate- or severe-weathering grade.

Journal of Solid Waste Technology and Management

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S 2- ) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

California

Earthquake rupture process recreated from a natural fault surface

What exactly happens on the rupture surface as an earthquake nucleates, spreads, and stops? We cannot observe this directly, and models depend on assumptions about physical conditions and geometry at depth. We thus measure a natural fault surface and use its 3D coordinates to construct a replica at 0.1 m resolution to obviate geometry uncertainty. We can recreate stick-slip behavior on the resulting finite element model that depends solely on observed fault geometry. We clamp the fault together and apply steady state tectonic stress until seismic slip initiates and terminates. Our recreated M~1 earthquake initiates at contact points where there are steep surface gradients because infinitesimal lateral displacements reduce clamping stress most efficiently there. Unclamping enables accelerating slip to spread across the surface, but the fault soon jams up because its uneven, anisotropic shape begins to juxtapose new high-relief sticking points. These contacts would ultimately need to be sheared off or strongly deformed before another similar earthquake could occur. Our model shows that an important role is played by fault-wall geometry, though we do not include effects of varying fluid pressure or exotic rheologies on the fault surfaces. We extrapolate our results to large fault systems using observed self-similarity properties, and suggest that larger ruptures might begin and end in a similar way, though the scale of geometrical variation in fault shape that can arrest a rupture necessarily scales with magnitude. In other words, fault segmentation may be a magnitude dependent phenomenon and could vary with each subsequent rupture.

Journal of Geophysical Research B: Solid Earth

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California