SEARCH · USGS Science

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2Linked to original sources

USGS2023

Two-dimensional inverse energy cascade in a laboratory surf zone for varying wave directional spread

Surfzone eddies enhance the dispersion and transport of contaminants, bacteria, and larvae across the nearshore, altering coastal water quality and ecosystem health. During directionally spread wave conditions, vertical vortices (horizontal eddies) are injected near the ends of breaking crests. Energy associated with these eddies may be transferred to larger-scale, low-frequency rotational motions through an inverse energy cascade, consistent with two-dimensional turbulence. However, our understanding of the relationships between the wave conditions and the dynamics and energetics of low-frequency surfzone eddies are largely based on numerical modeling. Here, we test these relationships with remotely sensed and in situ observations from large-scale directional wave basin experiments with varying wave conditions over alongshore-uniform barred bathymetry. Surface velocities derived with particle image velocimetry were employed to assess the spatial scales of low-frequency surfzone eddies and compute structure functions with alongshore velocities. Second-order structure functions for directionally spread waves ( ⁠

Physics of Fluids↗

USGS2012

Migration rates and formation injectivity to determine containment time scales of sequestered carbon dioxide

Supercritical carbon dioxide exhibits highly variable behavior over a range of reservoir pressure and temperature conditions. Because geologic sequestration of supercritical carbon dioxide is targeted for subsurface injection and containment at depths ranging from approximately 3,000 to 13,000 feet, the investigation into the physical properties of this fluid can be restricted to the pressure and temperature conditions likely encountered in the sedimentary strata within this depth interval. A petrophysical based approach was developed to study the widest range of formation properties potentially encountered in sedimentary strata. Fractional porosities were varied from 5 to 95 percent, in 5-percent increments, and permeability values were varied over thirteen orders of magnitude, from 10.0 darcys down to 1.0 picodarcy. Fluid-flow modeling incorporated two constitutive equations from fluid dynamics: hydraulic diffusivity for near-surface applications, and Darcy's Law for deeper formations exhibiting higher pressure gradients. Based on the flow modeling results, first-order approximations of carbon dioxide lateral migration rates were determined. These first-order approximations enable the establishment of a permeability classification system for dividing the subsurface into flow units that provide short, moderate, and long-term containment of carbon dioxide. These results enable a probabilistic determination of how fluids will enter and be contained in a subsurface storage formation, which is a vital step in the calculation of the carbon dioxide storage capacity of a reservoir. Additionally, this research establishes a methodology to calculate the injectivity of a target formation. Because injectivity describes the pressure increase due to the introduction of fluids into a formation, the relevant application of injectivity is to determine the pressure increase, due to an injection volume and flow rate, that will induce fractures in the reservoir rocks. This quantity is defined mathematically as the maximum pressure differential between the hydrostatic gradient and the fracture gradient of the target formation. Injectivity is mathematically related to the maximum pressure differential of the formation, and can be used to determine the upper limit for the pressure increase that an injection target can withstand before fracturing.

Open-File Report↗

USGS2008

Containment of fluid samples in the hydrothermal diamond-anvil cell without the use of metal gaskets: Performance and advantages for in situ analysis

Metal gaskets (Re, Ir, Inconel, or stainless steel) normally used to contain fluid samples in the hydrothermal diamond-anvil cell (HDAC) are sometimes undesirable due to possible contamination and to gasket deformation at high pressures and temperatures resulting in nonisochoric behavior. Furthermore, in x-ray spectroscopic experiments, metal gaskets may attenuate the incident x-ray beam and emitted fluorescence x-rays, and the interaction of scattered radiation with the gasket may produce fluorescence that interferes with the x-ray spectrum of the sample. New arrangements and procedures were tested for the operation of the HDAC without using the metal gaskets. Distilled, de-ionized water was loaded into the sample chamber, a laser-milled recess 300 ??m in diameter and ???50 ??m deep centered in the 1.0 mm face of the lower diamond anvil, and sealed by pressing the top diamond anvil face directly against the lower one without a metal gasket in between. A maximum sample pressure of 202 MPa at 617 ??C was maintained for a duration of 10 min without evidence of leakage. A small change in fluid density was observed in one experiment where the sample was held at 266 MPa at 708 ??C for 10 min. The gasketless HDAC was also employed in x-ray absorption spectroscopy experiments, where, in addition to the sample chamber in the lower diamond, two grooves were milled at a 90?? angle to each other around the sample chamber to minimize the attenuation of incident and fluorescent x rays. With a minimum distance between the sample chamber and the grooves of 80 ??m, a pressure of 76 MPa at 500 ??C was maintained for 2 h with no change in the original fluid density. ?? 2008 American Institute of Physics.

Review of Scientific Instruments↗

USGS2014

Debris flows: Behavior and hazard assessment

Debris flows are water-laden masses of soil and fragmented rock that rush down mountainsides, funnel into stream channels, entrain objects in their paths, and form lobate deposits when they spill onto valley floors. Because they have volumetric sediment concentrations that exceed 40 percent, maximum speeds that surpass 10 m/s, and sizes that can range up to ~109 m3, debris flows can denude slopes, bury floodplains, and devastate people and property. Computational models can accurately represent the physics of debris-flow initiation, motion and deposition by simulating evolution of flow mass and momentum while accounting for interactions of debris' solid and fluid constituents. The use of physically based models for hazard forecasting can be limited by imprecise knowledge of initial and boundary conditions and material properties, however. Therefore, empirical methods continue to play an important role in debris-flow hazard assessment.

Geology Today↗

USGS1978

Pressure gradients and boiling as mechanisms for localizing ore in porphyry systems

Fluid inclusions in ore zones of porphyry systems indicate that extensive boiling of hydrothermal fluids accompanies deposition of ore and gangue minerals. The boiling commonly accompanied a change from a lithostatic to a hydrostatic environment during evolution of an epizonal stock. Pressure gradients near the margin of the stock can determine whether ore or only a diffuse zone of mineralization is formed. A sharp drop in pressure in an epizonal environment is more likely to cause extensive boiling than a comparable change in a deeper environment, as the slope of the boiling curve steepens with an increase 'in pressure. The drop in pressure causes the hydrothermal fluids to boil and creates a crackle (stockwork) breccia, which hosts the veinlets of gangue quartz and ore minerals. The boiling selectively partitions CO 2 , H 2 S, and HCl into the vapor phase, changing the pH, composition, ionic strength, and thus the solubility product of metal complexes in the remaining liquid and causing the ore and gangue to come out of solution. Fluid inclusions trapped from boiling solutions can exhibit several forms, depending on the physical and chemical conditions of the hydrothermal fluid from which they were trapped. In one case, inclusions when heated can homogenize to either liquid or vapor at the same temperature, which is the true boiling temperature. In another case, homogenization of various inclusions can occur through a range of temperatures. The latter case results from the trapping of mixture of liquid and vapor. Variations in salinity can result from boiling of the hydrothermal fluid, or intermittent incorporation of high-salinity fluids from the magma, or trapping of fluids of varying densities at pressure-temperature conditions above the critical point of the fluid. In places, paleopressure-temperature transition zones can be recognized by fluid-inclusion homogenization temperatures and phase relationships and by the presence of anhydrite daughter minerals. Boiling of a hydrothermal fluid in the porphyry environment affects light stable isotopes. Hydrogen is preferentially fractionated into the vapor phase from water boiling below 223°C; above this temperature deuterium is selectively enriched in the vapor phase. In certain environments boiling creates a vapor-dominated system in which the condensate is swept away by meteoric waters and the H/D in the residual fluids is progressively increased through time.

Utah↗

USGS2021

Forecasting induced earthquake hazard using a hydromechanical earthquake nucleation model

In response to the dramatic increase in earthquake rates in the central United States, the U.S Geological Survey began releasing 1 yr earthquake hazard models for induced earthquakes in 2016. Although these models have been shown to accurately forecast earthquake hazard, they rely purely on earthquake statistics because there was no precedent for forecasting induced earthquakes based upon wastewater injection data. Since the publication of these hazard models, multiple physics‐based methods have been proposed to forecast earthquake rates using injection data. Here, we use one of these methods to generate earthquake hazard forecasts. Our earthquake hazard forecasts are more accurate than statistics‐based hazard forecasts. These results imply that fluid injection data, where and when available, and the physical implications of fluid injection should be included in future induced earthquake hazard forecasts.

Kansas, Oklahoma↗

USGS2025

Machine learning provides reconnaissance-type estimates of carbon dioxide storage resources in oil and gas reservoirs

Oil and gas reservoirs represent suitable containers to sequester carbon dioxide (CO 2 ) in a supercritical state because they are accessible, reservoir properties are known, and they previously contained stored buoyant fluids. However, planners must quantify the relative magnitude of the CO 2 storage resource in these reservoirs to formulate a comprehensive strategy for CO 2 mitigation. Even reconnaissance-type estimates of CO 2 storage resources of known oil and gas reservoirs may require complicated calculations involving 1) estimates of recoverable oil and gas, 2) reservoir properties (depth, temperature, pressure, etc.), and 3) the physical qualities of the retained fluids. We demonstrate the application of machine learning (ML) algorithms to bypass these computations to yield more rapid estimates of CO 2 storage resources in reservoirs capable of hosting CO 2 in a supercritical state. ML algorithms are computationally efficient because they do not impose the strong assumptions on the data-generating process that standard statistical or engineering procedures require. Further, ML algorithms can capture highly complex, particularly nonlinear, relationships among predictor variables. We demonstrate the application of four different ML algorithms using data from onshore and offshore oil and gas reservoirs in Europe, and show they perform well when predictions are compared to engineering estimates. The proposed methods and models provide an effective and novel way to more rapidly and directly determine the subsurface CO 2 storage capacity of oil and gas reservoirs around the world, information that operators, researchers, and policymakers alike require to meet energy transition and decarbonization goals.

Frontiers in Enviornmental Science↗

USGS1973

Hydraulic sand-model studies of miscible-fluid flow

Hydraulic sand models are useful physical tools in the investigation of the transition zone that occurs between salt and fresh ground water in coastal aquifers. Such models are used to demonstrate the effects of transport mechanisms that influence the dynamic behavior and the shape of the transition zone. The techniques employed in obtaining in-place measurements of solute concentrations are generally the stumbling block in generating data for two-dimensional dispersion systems. Two in-place measurement techniques were used in the studies described: (1) conductivity probes when salt was used as a tracer; and (2) photoelectric cells when organic dye was used as a tracer. Results indicate that conductivity methods are more reliable; however, care must be exercised inasmuch as the probes tend to disturb the fluid flow.

Journal of Research of the U.S. Geological Survey↗

USGS2004

Shear wave velocity, seismic attenuation, and thermal structure of the continental upper mantle

Seismic velocity and attenuation anomalies in the mantle are commonly interpreted in terms of temperature variations on the basis of laboratory studies of elastic and anelastic properties of rocks. In order to evaluate the relative contributions of thermal and non-thermal effects on anomalies of attenuation of seismic shear waves, Q −1 s , and seismic velocity, V s , we compare global maps of the thermal structure of the continental upper mantle with global Q −1 s and V s maps as determined from Rayleigh waves at periods between 40 and 150 s. We limit the comparison to three continental mantle depths (50, 100 and 150 km), where model resolution is relatively high. The available data set does not indicate that, at a global scale, seismic anomalies in the upper mantle are controlled solely by temperature variations. Continental maps have correlation coefficients of <0.56 between V s and T and of <0.47 between Q s and T at any depth. Such low correlation coefficients can partially be attributed to modelling artefacts; however, they also suggest that not all of the V s and Q s anomalies in the continental upper mantle can be explained by T variations. Global maps show that, by the sign of the anomaly, V s and Q s usually inversely correlate with lithospheric temperatures: most cratonic regions show high V s and Q s and low T , while most active regions have seismic and thermal anomalies of the opposite sign. The strongest inverse correlation is found at a depth of 100 km, where the attenuation model is best resolved. Significantly, at this depth, the contours of near-zero Q s anomalies approximately correspond to the 1000 °C isotherm, in agreement with laboratory measurements that show a pronounced increase in seismic attenuation in upper mantle rocks at 1000–1100 °C. East–west profiles of V s , Q s and T where continental data coverage is best (50°N latitude for North America and 60°N latitude for Eurasia) further demonstrate that temperature plays a dominant, but non-unique, role in determining the value of lithospheric V s and Q s . At 100 km depth, where the resolution of seismic models is the highest, we compare observed seismic V s and Q s with theoretical V T s and Q T s values, respectively, that are calculated solely from temperature anomalies and constrained by experimental data on temperature dependencies of velocity and attenuation. This comparison shows that temperature variations alone are sufficient to explain seismic V s and Q s in ca 50 per cent of continental regions. We hypothesize that compositional anomalies resulting from Fe depletion can explain the misfit between seismic and theoretical V s in cratonic lithosphere. In regions of active tectonics, temperature effects alone cannot explain seismic V s and Q s in the lithosphere. It is likely that partial melts and/or fluids may affect seismic parameters in these regions. This study demonstrates that lithospheric temperature plays the dominant role in controlling V s and Q s anomalies, but other physical parameters, such as compositional variations, fluids, partial melting and scattering, may also play a significant role in determining V s and Q s variations in the continental mantle.

Geophysical Journal International↗

USGS1987

Evidence for Late-Paleozoic brine migration in Cambrian carbonate rocks of the central and southern Appalachians: Implications for Mississippi Valley-type sulfide mineralization

Many Lower Paleozoic limestones and dolostones in the Valley and Ridge province of the central and southern Appalachians contain 10 to 25 weight percent authigenic potassium feldspar. This was considered to be a product of early diagenesis, however, 40 Ar 39 Ar "> 40 Ar 39 Ar analyses of overgrowths on detrital K-feldspar in Cambrian carbonate rocks from Pennsylvania, Maryland, Virginia, and Tennessee yield Late Carboniferous-Early Permian ages (278–322 Ma). Simple mass balance calculations suggest that the feldspar could not have formed isochemically, but required the flux of multiple pore volumes of fluid through the rocks, reflecting regional fluid migration events during the Late-Paleozoic Alleghanian orogeny. Microthermometric measurements of fluid inclusions in overgrowths on detrital K-feldspar and quartz grains from unmineralized rocks throughout the study area indicate homogenization temperatures from 100° to 200°C and freezing point depressions of −14° to −18.5°C (18–21 wt.% NaCl equiv). The apparent similarity of these fluids to fluid inclusions in ore and gangue minerals of nearby Mississippi Valley-type (MVT) deposits suggests that the regional occurrences of authigenic K-feldspar and MVT mineralization may be genetically related. This hypothesis is supported by the discovery of authigenic K-feldspar intergrown with sphalerite in several mines of the Mascot-Jefferson City District, E. Tennessee. Regional potassic alteration in unmineralized carbonate rocks and localized occurrences of MVT mineralization are both explainable by a gravity-driven flow model, in which deep brines migrate towards the basin margin under a hydraulic gradient established during the Alleghanian orogeny. The authigenic K-feldspar may reflect the loss of K during disequilibrium cooling of the ascending brines. MVT deposits are probably localized manifestations of the same migrating fluids, occurring where the necessary physical and chemical traps are present.

Geochimica et Cosmochimica Acta↗

USGS2011

Imaging the San Andreas fault with explosion and earthquake sources

Mounting evidence suggests that fault zone heterogeneity may play a crucial role in the localization of rupture in earthquakes [ Aki , 1995]. The heterogeneity can take several forms: spatial variations in physical properties (elastic properties, pore fluid pressure, etc.) or complexity in the fault surface (bends, offsets, etc.). High-resolution, three-dimensional models of the P and S wave velocity ( V p and V s ) structure and accurate hypocenters provide tools for studying the effects of fault zone heterogeneity on rupture localization. A dense, passive array of 48 seismic instruments was deployed south of Hollister, California, from mid-November 1994 to late May 1995. The array obtained local-earthquake seismograms for high-resolution, three-dimensional imaging of the ( V p and V s ) structure in the San Andreas fault zone. In mid-May, anactive seismic experiment was also carried out. Figure 1 shows a 7-km wide, 3- to 4-km deep zone of very low V p imaged during the experiment. The zone is bounded on the southwest by the San Andreas fault and the adjacent Northern Gabilan Range. The San Andreas fault appears to be approximately vertical to at least 5 km depth (the limit of good model resolution), contrary to an existing model based on gravity data.

California↗

USGS2006

Phase equilibria constraints on the chemical and physical evolution of the campanian ignimbrite

The Campanian Ignimbrite is a > 200 km3 trachyte-phonolite pyroclastic deposit that erupted at 39.3 ?? 0.1 ka within the Campi Flegrei west of Naples, Italy. Here we test the hypothesis that Campanian Ignimbrite magma was derived by isobaric crystal fractionation of a parental basaltic trachyandesitic melt that reacted and came into local equilibrium with small amounts (5-10 wt%) of crustal rock (skarns and foid-syenites) during crystallization. Comparison of observed crystal and magma compositions with results of phase equilibria assimilation-fractionation simulations (MELTS) is generally very good. Oxygen fugacity was approximately buffered along QFM+1 (where QFM is the quartz-fayalite-magnetite buffer) during isobaric fractionation at 0.15 GPa (???6 km depth). The parental melt, reconstructed from melt inclusion and host clinopyroxene compositions, is found to be basaltic trachyandesite liquid (51.1 wt% SiO2, 9.3 wt% MgO, 3 wt% H2O). A significant feature of phase equilibria simulations is the existence of a pseudo-invariant temperature, ???883??C, at which the fraction of melt remaining in the system decreases abruptly from ???0.5 to < 0.1. Crystallization at the pseudo-invariant point leads to abrupt changes in the composition, properties (density, dissolved water content), and physical state (viscosity, volume fraction fluid) of melt and magma. A dramatic decrease in melt viscosity (from 1700 Pa s to ???200 Pa s), coupled with a change in the volume fraction of water in magma (from ??? 0.1 to 0.8) and a dramatic decrease in melt and magma density acted as a destabilizing eruption trigger. Thermal models suggest a timescale of ??? 200 kyr from the beginning of fractionation until eruption, leading to an apparent rate of evolved magma generation of about 10-3 km3/year. In situ crystallization and crystal settling in density-stratified regions, as well as in convectively mixed, less evolved subjacent magma, operate rapidly enough to match this apparent volumetric rate of evolved magma production. ?? Copyright 2007 Oxford University Press.

Journal of Petrology↗

USGS2009

Geophysical log analysis of selected test holes and wells in the High Plains Aquifer, Central Platte River Basin, Nebraska

The U.S. Geological Survey in cooperation with the Central Platte Natural Resources District is investigating the hydrostratigraphic framework of the High Plains aquifer in the Central Platte River basin. As part of this investigation, a comprehensive set of geophysical logs was collected from six test holes at three sites and analyzed to delineate the penetrated stratigraphic units and characterize their lithology and physical properties. Flow and fluid-property logs were collected from two wells at one of the sites and analyzed along with the other geophysical logs to determine the relative transmissivity of the High Plains aquifer units. The integrated log analysis indicated that the coarse-grained deposits of the alluvium and the upper part of the Ogallala Formation contributed more than 70 percent of the total transmissivity at this site. The lower part of the Ogallala with its moderately permeable sands and silts contributed some measureable transmissivity, as did the fine-grained sandstone of the underlying Arikaree Group, likely as a result of fractures and bedding-plane partings. Neither the lower nor the upper part of the siltstone- and claystone-dominated White River Group exhibited measurable transmissivity. The integrated analysis of the geophysical logs illustrated the utility of these methods in the detailed characterization of the hydrostratigraphy of the High Plains aquifer.

Nebraska↗

USGS2005

Evaluation of ground-penetrating radar to detect free-phase hydrocarbons in fractured rocks: Results of numerical modeling and physical experiments

The suitability of common-offset ground-penetrating radar (GPR) to detect free-phase hydrocarbons in bedrock fractures was evaluated using numerical modeling and physical experiments. The results of one- and two-dimensional numerical modeling at 100 megahertz indicate that GPR reflection amplitudes are relatively insensitive to fracture apertures ranging from 1 to 4 mm. The numerical modeling and physical experiments indicate that differences in the fluids that fill fractures significantly affect the amplitude and the polarity of electromagnetic waves reflected by subhorizontal fractures. Air-filled and hydrocarbon-filled fractures generate low-amplitude reflections that are in-phase with the transmitted pulse. Water-filled fractures create reflections with greater amplitude and opposite polarity than those reflections created by air-filled or hydrocarbon-filled fractures. The results from the numerical modeling and physical experiments demonstrate it is possible to distinguish water-filled fracture reflections from air- or hydrocarbon-filled fracture reflections, nevertheless subsurface heterogeneity, antenna coupling changes, and other sources of noise will likely make it difficult to observe these changes in GPR field data. This indicates that the routine application of common-offset GPR reflection methods for detection of hydrocarbon-filled fractures will be problematic. Ideal cases will require appropriately processed, high-quality GPR data, ground-truth information, and detailed knowledge of subsurface physical properties. Conversely, the sensitivity of GPR methods to changes in subsurface physical properties as demonstrated by the numerical and experimental results suggests the potential of using GPR methods as a monitoring tool. GPR methods may be suited for monitoring pumping and tracer tests, changes in site hydrologic conditions, and remediation activities.The suitability of common-offset ground-penetrating radar (GPR) to detect free-phase hydrocarbons in bedrock fractures was evaluated using numerical modeling and physical experiments. The results of one- and two-dimensional numerical modeling at 100 megahertz indicate that GPR reflection amplitudes are relatively insensitive to fracture apertures ranging from 1 to 4 mm. The numerical modeling and physical experiments indicate that differences in the fluids that fill fractures significantly affect the amplitude and the polarity of electromagnetic waves reflected by subhorizontal fractures. Air-filled and hydrocarbon-filled fractures generate low-amplitude reflections that are in-phase with the transmitted pulse. Water-filled fractures create reflections with greater amplitude and opposite polarity than those reflections created by air-filled or hydrocarbon-filled fractures. The results from the numerical modeling and physical experiments demonstrate it is possible to distinguish water-filled fracture reflections from air- or hydrocarbon-filled fracture reflections, nevertheless subsurface heterogeneity, antenna coupling changes, and other sources of noise will likely make it difficult to observe these changes in GPR field data. This indicates that the routine application of common-offset GPR reflection methods for detection of hydrocarbon-filled fractures will be problematic. Ideal cases will require appropriately processed, high-quality GPR data, ground-truth information, and detailed knowledge of subsurface physical properties. Conversely, the sensitivity of GPR methods to changes in subsurface physical properties as demonstrated by the numerical and experimental results suggests the potential of using GPR methods as a monitoring tool. GPR methods may be suited for monitoring pumping and tracer tests, changes in site hydrologic conditions, and remediation activities.

Ground Water↗

USGS1986

A numerical investigation of head waves and leaky modes in fluid- filled boreholes

Although synthetic borehole seismograms can be computed for a wide range of borehole conditions, the physical nature of shear and compressional head waves in fluid-filled boreholes is poorly understood. This paper presents a series of numerical experiments designed to explain the physical mechanisms controlling head-wave propagation in boreholes. These calculations demonstrate the existence of compressional normal modes equivalent to shear normal modes, or pseudo-Rayleigh waves, with sequential cutoff frequencies spaced between the cutoff frequencies for the shear normal modes. Major contributions to head-wave spectra occur in discrete peaks at frequencies just below mode cutoff for both compressional and shear modes. This result is confirmed by calculations with synthetic waveforms at frequencies corresponding to mode cutoff, and by branch-cut integrals designed to yield independent spectra for the compressional mode. For soft formations where shear velocity falls below acoustic velocity in the borehole fluid, leaky compressional normal modes attain properties similar to those of shear normal modes in the case of hard rock. In the limit of vanishing S-wave velocity, this result is formally related to a fluid-fluid waveguide with undamped compressional normal modes. Synthetic waveforms demonstrate that high-amplitude arrivals, traveling at velocities less than the acoustic velocity of the borehole fluid and at frequencies above a few kilohertz, represent the Airy phase of the compressional mode and are not a tube wave. Comparison of synthetic waveforms with waveforms obtained in soft sea sediments indicates that the predicted Airy phase arrivals are present in the experimental data.

Geophysics↗

USGS1984

Systematics of stretching of fluid inclusions; I. Fluorite and sphalerite at 1 atmosphere confining pressure

More than 1,300 measurements on fluid inclusions in fluorite and sphalerite indicate that stretching proceeds systematically and predictably. In order to generate internal pressures that are sufficiently high to cause stretching, most inclusions must be heated beyond their initial homogenization temperatures, i.e., overheated. The amount of overheating necessary to initiate stretching depends on the P-V-T-X properties of the inclusion fluid, the inclusion size and shape, physical properties of the host mineral, and the confining pressure. In the range of homogenization temperatures and salinities examined in this study, the amount of overheating necessary to initiate stretching of inclusions in fluorite is inversely related to inclusion volume. The results of this study show that the measured homogenization temperature of a fluid inclusion may be considerably higher than the true homogenization temperature if the internal pressure reached sufficiently high values during previous testing in the laboratory or, less likely, during postentrapment thermal events in nature. The systematic relationship between the internal pressure necessary to initiate stretching and the inclusion volume provides a means of recognizing previously stretched inclusions and estimating the magnitude of postentrapment thermal events. Furthermore, reproducibility of homogenization temperature measurements, lack of microscopically observable fractures in the inclusion walls, and the lack of a noticeable increase in the vapor bubble size after heating are insufficient to prove that an inclusion has not stretched.

Economic Geology↗

USGS2024

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science↗

USGS2012

The physical hydrogeology of ore deposits

Hydrothermal ore deposits represent a convergence of fluid flow, thermal energy, and solute flux that is hydrogeologically unusual. From the hydrogeologic perspective, hydrothermal ore deposition represents a complex coupled-flow problem—sufficiently complex that physically rigorous description of the coupled thermal (T), hydraulic (H), mechanical (M), and chemical (C) processes (THMC modeling) continues to challenge our computational ability. Though research into these coupled behaviors has found only a limited subset to be quantitatively tractable, it has yielded valuable insights into the workings of hydrothermal systems in a wide range of geologic environments including sedimentary, metamorphic, and magmatic. Examples of these insights include the quantification of likely driving mechanisms, rates and paths of fluid flow, ore-mineral precipitation mechanisms, longevity of hydrothermal systems, mechanisms by which hydrothermal fluids acquire their temperature and composition, and the controlling influence of permeability and other rock properties on hydrothermal fluid behavior. In this communication we review some of the fundamental theory needed to characterize the physical hydrogeology of hydrothermal systems and discuss how this theory has been applied in studies of Mississippi Valley-type, tabular uranium, porphyry, epithermal, and mid-ocean ridge ore-forming systems. A key limitation in the computational state-of-the-art is the inability to describe fluid flow and transport fully in the many ore systems that show evidence of repeated shear or tensional failure with associated dynamic variations in permeability. However, we discuss global-scale compilations that suggest some numerical constraints on both mean and dynamically enhanced crustal permeability. Principles of physical hydrogeology can be powerful tools for investigating hydrothermal ore formation and are becoming increasingly accessible with ongoing advances in modeling software.

Economic Geology↗