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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13Linked to original sources

USGS2005

Chemical openness and potential for misinterpretation of the solute environment of coastal sabkhat

Sabkha deposits in the geologic record are commonly used to interpret the environmental conditions of deposition. Implicit in this use is the assumption that the solute system is chemically closed, that is, the authigenic minerals represent the composition of the fluids in their environment of origin. Thermodynamic and mass-balance calculations based on measurements of water and solute flux of contemporary Abu Dhabi coastal sabkha system, however, demonstrate that the system is open for sodium and chloride, where nearly half of the input is lost, but closed for sulfur, where nearly 100% is retained. Sulfur and chloride isotopes were consistent with this observation. If these sabkha deposits were preserved in the geologic record, they would suggest a solute environment rich in sulfate and poor in chloride; yet the reverse is true. In most coastal-sabkha environments, capillary forces bring solutes and water to the surface, where the water evaporates and halite, carnallite, sylvite, and other soluble minerals are precipitated. Retrograde minerals, such as anhydrite, calcite, dolomite, and gypsum, however, precipitate and accumulate in the capillary zone beneath the surface of the coastal sabkha. Because they possess relatively low solubility and are below the surface, these retrograde minerals are protected from dissolution and physical erosion occurring from infrequent but intense rainfall events. Thus, they are more likely to be preserved in the geological record than highly soluble minerals formed on the surface. ?? 2004 Elsevier B.V. All rights reserved.

Chemical Geology↗

USGS2010

Parametric study of the physical properties of hydrate‐bearing sand, silt, and clay sediments: 2. Small‐strain mechanical properties

The small‐strain mechanical properties (e.g., seismic velocities) of hydrate‐bearing sediments measured under laboratory conditions provide reference values for calibration of logging and seismic exploration results acquired in hydrate‐bearing formations. Instrumented cells were designed for measuring the compressional (P) and shear (S) velocities of sand, silts, and clay with and without hydrate and subject to vertical effective stresses of 0.01 to 2 MPa. Tetrahydrofuran (THF), which is fully miscible in water, was used as the hydrate former to permit close control over the hydrate saturation S hyd and to produce hydrate from dissolved phase, as methane hydrate forms in most natural marine settings. The results demonstrate that laboratory hydrate formation technique controls the pattern of P and S velocity changes with increasing S hyd and that the small‐strain properties of hydrate‐bearing sediments are governed by effective stress, σ ′ v and sediment specific surface. The S velocity increases with hydrate saturation owing to an increase in skeletal shear stiffness, particularly when hydrate saturation exceeds S hyd ≈ 0.4. At very high hydrate saturations, the small strain shear stiffness is determined by the presence of hydrates and becomes insensitive to changes in effective stress. The P velocity increases with hydrate saturation due to the increases in both the shear modulus of the skeleton and the bulk modulus of pore‐filling phases during fluid‐to‐hydrate conversion. Small‐strain Poisson's ratio varies from 0.5 in soft sediments lacking hydrates to 0.25 in stiff sediments (i.e., subject to high vertical effective stress or having high S hyd ). At S hyd ≥ 0.5, hydrate hinders expansion and the loss of sediment stiffness during reduction of vertical effective stress, meaning that hydrate‐rich natural sediments obtained through pressure coring should retain their in situ fabric for some time after core retrieval if the cores are maintained within the hydrate stability field.

Journal of Geophysical Research↗

USGS2014

DigitalCrust – a 4D data system of material properties for transforming research on crustal fluid flow

Fluid circulation in the Earth's crust plays an essential role in surface, near surface, and deep crustal processes. Flow pathways are driven by hydraulic gradients but controlled by material permeability, which varies over many orders of magnitude and changes over time. Although millions of measurements of crustal properties have been made, including geophysical imaging and borehole tests, this vast amount of data and information has not been integrated into a comprehensive knowledge system. A community data infrastructure is needed to improve data access, enable large-scale synthetic analyses, and support representations of the subsurface in Earth system models. Here, we describe the motivation, vision, challenges, and an action plan for a community-governed, four-dimensional data system of the Earth's crustal structure, composition, and material properties from the surface down to the brittle–ductile transition. Such a system must not only be sufficiently flexible to support inquiries in many different domains of Earth science, but it must also be focused on characterizing the physical crustal properties of permeability and porosity, which have not yet been synthesized at a large scale. The DigitalCrust is envisioned as an interactive virtual exploration laboratory where models can be calibrated with empirical data and alternative hypotheses can be tested at a range of spatial scales. It must also support a community process for compiling and harmonizing models into regional syntheses of crustal properties. Sustained peer review from multiple disciplines will allow constant refinement in the ability of the system to inform science questions and societal challenges and to function as a dynamic library of our knowledge of Earth's crust.

Geofluids↗

USGS1976

Submarine geothermal resources

Approximately 20% of the earth's heat loss (or 2 ?? 1012 cal/s) is released through 1% of the earth's surface area and takes the form of hydrothermal discharge from young (Pleistocene or younger) rocks adjacent to active seafloor-spreading centers and submarine volcanic areas. This amount is roughly equivalent to man's present gross energy consumption rate. A sub-seafloor geothermal reservoir, to be exploitable under future economic conditions, will have to be hot, porous, permeable, large, shallow, and near an energy-deficient, populated land mass. Furthermore, the energy must be recoverable using technology achievable at a competitive cost and numerous environmental, legal and institutional problems will have to be overcome. The highest-temperature reservoirs should be found adjacent to the zones of the seafloor extension or volcanism that are subject to high sedimentation rates. The relatively impermeable sediments reduce hydrothermal-discharge flow rates, forcing the heat to be either conducted away or released by high-temperature fluids, both of which lead to reservoir temperatures that can exceed 300??C. There is evidence that the oceanic crust is quite permeable and porous and that it was amenable to deep (3-5 km) penetration by seawater at least some time in the early stages of its evolution. Most of the heat escapes far from land, but there are notable exceptions. For example, in parts of the Gulf of California, thermal gradients in the bottom sediments exceed 1??C/m. In the coastal areas of the Gulf of California, where electricity and fresh water are at a premium, this potential resource lies in shallow water (< 200 m) and within sight of land. Other interesting areas include the Sea of Japan, the Sea of Okhotsk and the Andaman Sea along the margins of the western Pacific, the Tyrrhenian Sea west of Italy, and the southern California borderland and west flank of the Juan de Fuca Ridge off the west coast of the United States. Many questions remain to be answered about the physical and other characteristics of these systems before they can be considered a viable resource. Until several of the most promising areas are carefully defined and drilled, the problem will remain unresolved. ?? 1976.

Journal of Volcanology and Geothermal Research↗

USGS2006

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software↗

USGS2013

Compaction and gas loss in welded pyroclastic deposits as revealed by porosity, permeability, and electrical conductivity measurements of the Shevlin Park Tuff

Pyroclastic flows produced by large volcanic eruptions commonly densify after emplacement. Processes of gas escape, compaction, and welding in pyroclastic-flow deposits are controlled by the physical and thermal properties of constituent material. Through measurements of matrix porosity, permeability, and electrical conductivity, we provide a framework for understanding the evolution of pore structure during these processes. Using data from the Shevlin Park Tuff in central Oregon, United States, and from the literature, we find that over a porosity range of 0%–70%, matrix permeability varies by almost 10 orders of magnitude (from 10 –20 to 10 –11 m 2 ), with over three orders of magnitude variation at any given porosity. Part of the variation at a given porosity is due to permeability anisotropy, where oriented core samples indicate higher permeabilities parallel to foliation (horizontally) than perpendicular to foliation (vertically). This suggests that pore space is flattened during compaction, creating anisotropic crack-like networks, a geometry that is supported by electrical conductivity measurements. We find that the power law equation: k 1 = 1.3 × 10 –21 × ϕ 5.2 provides the best approximation of dominant horizontal gas loss, where k 1 = permeability, and ϕ = porosity. Application of Kozeny-Carman fluid-flow approximations suggests that permeability in the Shevlin Park Tuff is controlled by crack- or disk-like pore apertures with minimum widths of 0.3 and 7.5 μm. We find that matrix permeability limits compaction over short times, but deformation is then controlled by competition among cooling, compaction, water resorption, and permeable gas escape. These competing processes control the potential for development of overpressure (and secondary explosions) and the degree of welding in the deposit, processes that are applicable to viscous densification of volcanic deposits in general. Further, the general relationships among porosity, permeability, and pore geometry are relevant for flow of any fluid through an ignimbritic host.

Oregon↗

USGS2025

Impact of gas/liquid phase change of CO2 during injection for sequestration

CO 2 sequestration in deep saline formations is an effective and important process to control the rapid rise in CO 2 emissions. The process of injecting CO 2 requires reliable predictions of the stress in the formation and the fluid pressure distributions – particularly since monitoring of the CO 2 migration is difficult – to mitigate leakage, prevent induced seismicity, and analyze wellbore stability. A key aspect of CO 2 is the gas–liquid phase transition at the temperatures and pressures of relevance to leakage and sequestration, which has been recognized as being critical for accurate predictions but has been challenging to model without ad hoc empiricisms. This paper presents a robust multiphase thermodynamics-based poromechanics model to capture the complex phase transition behavior of CO 2 and predict the stress and pressure distribution under super- and sub- critical conditions during the injection process. A finite element implementation of the model is applied to analyze the behavior of a multiphase porous system with CO 2 as it displaces the fluid brine phase. We find that if CO 2 undergoes a phase transition in the geologic reservoir, the spatial variation of the density is significantly affected, and the migration mobility of CO 2 decreases in the reservoir. A key feature of our approach is that we do not a priori assume the location of the CO 2 gas/liquid interface – or even if it occurs at all – but rather, this is a prediction of the model, along with the spatial variation of the phase of CO 2 and the change of the saturation profile due to the phase change.

Journal of the Mechanics and Physics of Solids↗

USGS2009

The length of channelized lava flows: Insight from the 1859 eruption of Mauna Loa Volcano, Hawai‘i

The 1859 eruption of Mauna Loa Volcano, Hawai'i, produced paired 'a'ā and pāhoehoe flows of exceptional length (51 km). The 'a'ā flow field is distinguished by a long (> 36 km) and well-defined pāhoehoe-lined channel, indicating that channelized lava remained fluid to great distances from the vent. The 1859 eruption was further unusual in initiating at a radial vent on the volcano's northwest flank, instead of along the well-defined rift zone that has been the source of most historic activity. As such, it presents an opportunity both to examine controls on the emplacement of long lava channels and to assess hazards posed by future flank eruptions of Mauna Loa. Here we combine evidence from historical chronicles with analysis of bulk compositions, glass geothermometry, and microlite textures of samples collected along the 1859 lava flows to constrain eruption and flow emplacement conditions. The bulk compositions of samples from the 'a'ā and pāhoehoe flow fields are bimodally distributed and indicate tapping of two discrete magma bodies during eruption. Samples from the pāhoehoe flow field have bulk compositions similar to those of historically-erupted lavas (< 8 wt.% MgO); lava that fed the 'a'ā channel is more primitive (> 8 wt.% MgO), nearly aphyric, and was erupted at high temperatures (1194–1216 °C). We suggest that the physical properties of proximal channel-fed lava (i.e., high-temperature, low crystallinity, and low bulk viscosity) promoted both rapid flow advance and development of long pāhoehoe-lined channels. Critical for the latter was the large temperature decrease (~ 50 °C) required to reach the point at which plagioclase and pyroxene started to crystallize; the importance of phase constraints are emphasized by our difficulty in replicating patterns of cooling and crystallization recorded by high-temperature field samples using common models of flow emplacement. Placement of the 1859 eruption within the context of historic activity at Mauna Loa suggests that the formation of radial vents and eruptions of high-temperature magma may not only be linked, but may also be a consequence of periods of high magma supply (e.g., 1843–1877). Flank eruptions could therefore warrant special consideration in models and hazard mitigation efforts.

Hawai'i↗

USGS2014

Modeling the influence of organic acids on soil weathering

Biological inputs and organic matter cycling have long been regarded as important factors in the physical and chemical development of soils. In particular, the extent to which low molecular weight organic acids, such as oxalate, influence geochemical reactions has been widely studied. Although the effects of organic acids are diverse, there is strong evidence that organic acids accelerate the dissolution of some minerals. However, the influence of organic acids at the field-scale and over the timescales of soil development has not been evaluated in detail. In this study, a reactive-transport model of soil chemical weathering and pedogenic development was used to quantify the extent to which organic acid cycling controls mineral dissolution rates and long-term patterns of chemical weathering. Specifically, oxalic acid was added to simulations of soil development to investigate a well-studied chronosequence of soils near Santa Cruz, CA. The model formulation includes organic acid input, transport, decomposition, organic-metal aqueous complexation and mineral surface complexation in various combinations. Results suggest that although organic acid reactions accelerate mineral dissolution rates near the soil surface, the net response is an overall decrease in chemical weathering. Model results demonstrate the importance of organic acid input concentrations, fluid flow, decomposition and secondary mineral precipitation rates on the evolution of mineral weathering fronts. In particular, model soil profile evolution is sensitive to kaolinite precipitation and oxalate decomposition rates. The soil profile-scale modeling presented here provides insights into the influence of organic carbon cycling on soil weathering and pedogenesis and supports the need for further field-scale measurements of the flux and speciation of reactive organic compounds.

California↗

USGS2014

Integration of seismic-reflection and well data to assess the potential impact of stratigraphic and structural features on sustainable water supply from the Floridan aquifer system, Broward County, Florida

The U.S. Geological Survey and Broward County water managers commenced a 3.5-year cooperative study in July 2012 to refine the geologic and hydrogeologic framework of the Floridan aquifer system (FAS) in Broward County. A lack of advanced stratigraphic knowledge of the physical system and structural geologic anomalies (faults and fractures originating from tectonics and karst-collapse structures) within the FAS pose a risk to the sustainable management of the resource. The principal objective of the study is to better define the regional stratigraphic and structural setting of the FAS in Broward County. The objective will be achieved through the acquisition, processing, and interpretation of new seismic-reflection data along several canals in Broward County. The interpretation includes integration of the new seismic-reflection data with existing seismic-reflection profiles along Hillsboro Canal in Broward County and within northeast Miami-Dade County, as well as with data from nearby FAS wellbores. The scope of the study includes mapping the geologic, hydrogeologic, and seismic-reflection framework of the FAS, and identifying stratigraphic and structural characteristics that could either facilitate or preclude the sustainable use of the FAS as an alternate water supply or a treated effluent repository. In addition, the investigation offers an opportunity to: (1) improve existing groundwater flow models, (2) enhance the understanding of the sensitivity of the groundwater system to well-field development and upconing of saline fluids, and (3) support site selection for future FAS projects, such as Class I wells that would inject treated effluent into the deep Boulder Zone.

Florida↗

USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗

USGS2019

Physical property characteristics of gas hydrate-bearing reservoir and associated seal sediments collected during NGHP-02 in the Krishna-Godavari Basin, in the offshore of India

India’s National Gas Hydrate Program Expedition 02 (NGHP-02), was conducted to better understand geologic controls on gas hydrate occurrence and morphology, targeting potentially coarse-grained sediments near the base of the continental slope offshore eastern India. This study combines seismic, logging-while-drilling data, and a petroleum systems approach to provide a regional geologic context for the core- and grain-scale analyses. This multi-scale approach provides insight on the gas hydrate distribution, morphology and anticipated system response to depressurization-induced methane extraction. The study area, NGHP-02 Area B in the Krishna-Godavari Basin, contains a buried anticline/syncline structure that hosts fracture-filling gas hydrate in fine-grained sediment overlying coarser sediments with pore-occupying gas hydrate. Core- and grain-scale measurements show fine-grained sediment exerts a primary control on the distribution and morphology of gas hydrate in Area B. Diatoms in the fine-grained overburden cause porosity to increase with depth, reaching ~70% at the underlying reservoir contact. High porosity, combined with near-vertical faults, suggests the overlying sediment is an imperfect seal. This allows methane to escape the gas hydrate reservoir sediments and form primarily grain-displacing gas hydrate veins in the fine-grained overburden. Within the reservoir, fine-grained layers are interbedded with coarser-grained gas hydrate reservoir sands. Even in the reservoir sands, however, a soil classification study shows the fines content is high enough to control hydraulic and mechanical properties, such as permeability, compressibility and shear strength. Fluid motion during methane extraction from gas hydrates can mobilize those fines, which can then clog pore throats, limiting production rates. Pore-water freshening during gas hydrate dissociation can increase fines mobilization, particularly given the smectite identified in the fine-grained interbeds. Accounting for fines content and specific fines mineralogy throughout the gas hydrate petroleum system is important for predicting production efficiency from gas hydrate occurrences along the crest of the anticline in NGHP-02 Area B.

Marine and Petroleum Geology↗

USGS2019

Multi-scale preferential flow processes in an urban streambed under variable hydraulic conditions

Spatially preferential flow processes occur at nested scales at the sediment-water interface (SWI), due in part to sediment heterogeneities, which may be enhanced in flashy urban streams with heavy road sand influence. However, several factors, including the flow-rate dependence of preferential hyporheic flow and discrete groundwater discharge zones are commonly overlooked in reach-scale models of groundwater/surface water exchange. Using a series of controlled-head tracer-injection experiments coupled with cm-scale geophysics within the highly reactive upper 30 cm of the hyporheic zone of an urban stream, we quantified the flow dependence of local less-mobile porosity volume, mass-transfer rate coefficient, and the resulting local residence time in the less-mobile pore space at three controlled downward fluid fluxes (0.8, 2, and 3 m/d). Experiments were performed in two adjacent streambed locations, representing different sediment bulk vertical permeability. Less-mobile porosity parameters were generally substantial and similar between the two streambed locations; though a more competent, thin, organic layer at ∼15 cm depth in one location strongly impacted tracer loading, flushing dynamics, and local residence times. Increased downward flux led to (1) a decrease in less-mobile porosity residence time in all experiments, and (2) an increase in less-mobile porosity fraction for most experiments. Additionally, at the larger stream reach-scale, surface electrodes for electrical resistivity measurement were installed along 22 m of the wetted stream channel. These surface electrode measurements were collected during a natural storm flow event, which revealed widespread, short-term, flushing (e.g. <3 h) of the hyporheic zone with stream water, followed by longer-term (e.g. >60 h) flushing of the SWI with riparian zone groundwater. Flow dependence of preferential hyporheic zone flowpaths, like in the controlled tracer experiments, was also observed in these reach-scale electrical resistivity tomography measurements. Our findings reveal that the spatial and temporal dependence of preferential flow processes create highly dynamic SWI conditions that will affect the physical and coupled biogeochemical functions of the SWI in urbanized, sand-impacted streams.

Journal of Hydrology↗

USGS2024

Energy-related rare earth element sources

Energy-related materials such as coal, coal-bearing wastes, and coal combustion products are traditionally thought of as sources or by-products of electric power generation. Increasingly, these materials are considered resources for their content of rare earth elements (REEs) and other useful constituents. In this chapter, we examine the distribution, modes of occurrence, and relative extractability of REEs from coal-derived materials. We also consider economic factors associated with recovery of REEs from these sources. While several coal-derived sources show promise for REE recovery at the pilot scale, in all cases, REE contents are much below those of primary ores, such that extraction and concentrating the REEs require new and innovative approaches that are largely developmental. Among coal-related sources, fly ash is the most REE-enriched, as REEs from coal are strongly retained in these refractory solids remaining after coal combustion. Partitioning of coal-derived elements into fly ash has been known for decades but this has yet to be commercially exploited. A key drawback shown in this chapter is that a significant fraction of REEs in fly ash is contained in highly insoluble aluminosilicate glasses that make up the largest portion of this material. In addition to testing chemical or physical pretreatment approaches to help improve the extractability of REEs from fly ash, current research is applying modern analytical approaches to better understand the distribution of REEs on increasingly smaller scales, in the interest of targeting their recovery. Next-most REE-enriched among coal-related materials are solid waste products of coal mining and wastes from coal preparation, both of which are REE-enriched relative to coal itself. These waste coals concentrate mineralogical constituents that are excluded during mining or removed during coal preparation because they do not contribute to the heating value of coal for power generation. Recovery of REEs from coal waste has shown promise at the pilot scale and has the added benefit of converting a waste into useful constituents. Total REE contents of commercial coals are, on average, much below the 300 parts per million interest level for REE recovery set by the U.S. Department of Energy (DOE). However, as reviewed in this chapter, certain horizons within coal beds show preferential REE enrichment and could be targeted by selective mining. Beyond this, certain coals are REE-enriched overall due to their unique geologic histories involving derivation from REE-enriched sediment sources, deposition of volcanic ash during coal formation, or interaction of coal with REE-bearing fluids. Acidic drainage from abandoned coal mines is produced by the breakdown of pyrite (FeS 2 ), which is unstable in oxygenated conditions. While these acidic fluids have lower REE contents than any of the coal-based solids described above, they are proportionally enriched in certain heavy rare earths, especially yttrium (Y). Precipitates from coal-based acid-mine drainage concentrate REEs to levels that are of interest for recovery, and these are also promising sources for extraction at the pilot scale.

Book chapter↗

USGS2011

Physical properties of sediment from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

This study characterizes cored and logged sedimentary strata from the February 2007 BP Exploration Alaska, Department of Energy, U.S. Geological Survey (BPXA-DOE-USGS) Mount Elbert Gas Hydrate Stratigraphic Test Well on the Alaska North Slope (ANS). The physical-properties program analyzed core samples recovered from the well, and in conjunction with downhole geophysical logs, produced an extensive dataset including grain size, water content, porosity, grain density, bulk density, permeability, X-ray diffraction (XRD) mineralogy, nuclear magnetic resonance (NMR), and petrography. This study documents the physical property interrelationships in the well and demonstrates their correlation with the occurrence of gas hydrate. Gas hydrate (GH) occurs in three unconsolidated, coarse silt to fine sand intervals within the Paleocene and Eocene beds of the Sagavanirktok Formation: Unit D-GH (614.4 m-627.9 m); unit C-GH1 (649.8 m-660.8 m); and unit C-GH2 (663.2 m-666.3 m). These intervals are overlain by fine to coarse silt intervals with greater clay content. A deeper interval (unit B) is similar lithologically to the gas-hydrate-bearing strata; however, it is water-saturated and contains no hydrate. In this system it appears that high sediment permeability (k) is critical to the formation of concentrated hydrate deposits. Intervals D-GH and C-GH1 have average "plug" intrinsic permeability to nitrogen values of 1700 mD and 675 mD, respectively. These values are in strong contrast with those of the overlying, gas-hydrate-free sediments, which have k values of 5.7 mD and 49 mD, respectively, and thus would have provided effective seals to trap free gas. The relation between permeability and porosity critically influences the occurrence of GH. For example, an average increase of 4% in porosity increases permeability by an order of magnitude, but the presence of a second fluid (e.g., methane from dissociating gas hydrate) in the reservoir reduces permeability by more than an order of magnitude.

Alaska↗

USGS2021

Velocity-porosity relations in carbonate and siliciclastic subduction zone input materials

The mechanical, physical, and frictional properties of incoming materials play an important role in subduction zone structure and slip behavior because these properties influence the strength of the accretionary wedge and megathrust plate boundary faults. Incoming sediment sections often show an increase in compressional wave speed (Vp) and a decrease in porosity with depth due to consolidation. These relations allow seismic-velocity models to be used to elucidate properties and conditions at depth. However, variations in these properties are controlled by lithology and composition as well as cementation and diagenesis. We present an analysis of shipboard measurements of Vp and porosity on incoming sediment cores from International Ocean Discovery Program (IODP) expeditions at the Hikurangi Margin, Nankai Trough, Aleutian Trench, Middle America Trench, and Sunda Trench. Porosity for these samples ranges from 5% to 85% and Vp ranges from 1.5 to 6 km/s. Vp-porosity relations developed by Erikson & Jarrad (1998), https://doi.org/10.1029/98JB02128 and Hoffman & Tobin (2004) https://10.2973/odp.proc.sr.190196.355.2004, with a critical porosity of ∼30%, can represent carbonate-poor (<50 wt% CaCO3), mainly hemipelagic, incoming sediment regardless of the margin. But these relations tend to underestimate porosity in incoming sediments with carbonate content greater than 50 wt%, which appear to have a critical porosity of between 45% and 50%. This discrepancy will lead to inaccuracy in estimates of fluid budget and overpressure in subduction zones. The velocity-porosity relation in carbonate sediments is non-unique due to the complexity that results from the greater susceptibility of carbonate rocks to diagenetic processes.

Geochemistry, Geophysics, Geosystems↗

USGS2000

Late-kinematic timing of orogenic gold deposits and significance for computer-based exploration techniques with emphasis on the Yilgarn Block, Western Australia

Orogenic gold deposits are a widespread coherent group of epigenetic ore deposits that are sited in accretionary or collisional orogens. They formed over a large crustal-depth range from deep-seated low-salinity H 2 O–CO 2 ±CH 4 ±N 2 ore fluids and with Au transported as thio-complexes. Regional structures provide the main control on deposit distribution. In many terranes, first-order faults or shear zones appear to have controlled regional fluid flow, with greatest ore-fluid fluxes in, and adjacent to, lower-order faults, shear zones and/or large folds. Highly competent and/or chemically reactive rocks are the most common hosts to the larger deposits. Focusing of supralithostatic ore fluids into dilatant zones appears to occur late during the evolutionary history of the host terranes, normally within D 3 or D 4 in a D 1 –D 4 deformation sequence. Reactivation of suitably oriented pre-existing structures during a change in far-field stress orientation is a factor common to many deposits, and repeated reactivation may account for multiple mineralization episodes in some larger deposits. Absolute robust ages of mineralization support their late-kinematic timing, and, in general, suggest that deposits formed diachronously towards the end of the 100 to 200 m.y. long evolutionary history of hosting orogens. For example, in the Yilgarn Block, a region specifically emphasised in this study, orogenic gold deposits formed in the time interval between 40 and 90 m.y., with most about 60 to 70 m.y., after the youngest widespread basic-ultrabasic volcanism and towards the end of felsic magmatism. The late timing of orogenic gold deposits is pivotal to geologically-based exploration methodologies. This is because the present structural geometries of: (i) the deposits, (ii) the hosting goldfields, and (iii) the enclosing terranes are all essentially similar to those during gold mineralization, at least in their relative position to each other. Thus, interpretation of geological maps and cross-sections and three-dimensional models can be used to accurately simulate the physical conditions that existed at the time of ore deposition. It is particularly significant that the deposits are commonly related to repetitive and predictable geometries, such as structural heterogeneities within or adjacent to first-order structures, around rigid granitoid bodies, or in specific “locked-up” fold-thrust structures. Importantly, the two giant greenstone-hosted goldfields, Kalgoorlie and Timmins, show a remarkably similar geometry at the regional scale. Computer-based stress mapping and GIS-based prospectivity mapping are two computer-based quantitative methodologies that can utilize and take advantage of the late timing aspect of this deposit type to provide important geological aids in exploration, both in broad regions and more localized goldfields. Both require an accurate and consistent solid geology map, stress mapping requires knowledge of the far-field stresses during mineralization, and the empirical prospectivity mapping requires data from a significant number of known deposits in the terrane. The Kalgoorlie Terrane, in the Yilgarn Block, meets these criteria, and illustrates the potential of these methodologies in the exploration for orogenic gold deposits. Low minimum stress anomalies, interpreted to represent dilational zones during gold-related deformation, coincide well with the positions of known goldfields rather than individual gold deposits in the terrane, and there are additional as-yet unexplained anomalies. The prospectivity analysis confirms that predictable and repetitive factors controlling the siting of deposits are: (i) proximity to, and orientation and curvature of, granitoid-greenstone contacts, (ii) proximity to segments of crustal faults which strike in a preferred direction, (iii) proximity to specific lithological contacts which have similar preferred strike, (iv) proximity to anticlinal structures, and (v) the presence of preferred reactive host rocks (e.g., dolerite). The prospectivity map defines a series of anomalous areas, which broadly conform to those of the stress map (>78% correspondence). The most prospective category on this map covers less than 0.3% of the greenstone belts and yet hosts 16% of the known deposits, which have produced>80% of known gold. Thus, it discriminates in favour of the larger economically more-attractive deposits in the terrane. The successful application of stress mapping and prospectivity mapping to geology-based exploration for orogenic gold deposits indicates that more quantitative analysis of geological map data is a profitable line of research. The computer-based nature of these methodologies is ideal for the production of an ultimate, integrated, deposit target map, which can be compared to other, more conventional, targeting parameters such as geophysical and geochemical anomalies. Such an integrated strategy appears the way forward in the increasingly difficult task of cost-effective global exploration for orogenic gold deposits in poorly exposed terranes.

Western Australia↗

USGS1996

Simulation of subsurface storage and recovery of treated effluent injected in a saline aquifer, St. Petersburg, Florida

The potential for subsurface storage and recovery of treated effluent into the uppermost producing zone (zone A) of the Upper Floridan aquifer in St. Petersburg, Florida, is being studied by the U.S. Geological Survey, in cooperation with the city of St. Petersburg and the Southwest Florida Water Management District. A measure of the success of this practice is the recovery efficiency, or the quantity of water relative to the quantity injected, that can be recovered before the water that is withdrawn fails to meet water-quality standards. The feasibility of this practice will depend upon the ability of the injected zone to receive, store, and discharge the injected fluid. A cylindrical model of ground-water flow and solute transport, incorporating available data on aquifer properties and water quality, was developed to determine the relation of recovery efficiency to various aquifer and fluid properties that could prevail in the study area. The reference case for testing was a base model considered representative of the saline aquifer underlying St. Petersburg. Parameter variations in the tests represent possible variations in aquifer conditions in the area. The model also was used to study the effect of various cyclic injection and withdrawal schemes on the recovery efficiency of the well and aquifer system. A base simulation assuming 15 days of injection of effluent at a rate of 1.0 million gallons per day and 15 days of withdrawal at a rate of 1.0 million gallons per day was used as reference to compare changes in various hydraulic and chemical parameters on recovery efficiency. A recovery efficiency of 20 percent was estimated for the base simulation. For practical ranges of hydraulic and fluid properties that could prevail in the study area, the model analysis indicates that (1) the greater the density contrast between injected and resident formation water, the lower the recovery efficiency, (2) recovery efficiency decreases significantly as dispersion increases, (3) high formation permeability favors low recovery efficiencies, and (4) porosity and anisotropy have little effect on recovery efficiencies. In several hypothetical tests, the recovery efficiency fluctuated between about 4 and 76 percent. The sensitivity of recovery efficiency to variations in the rate and duration of injection (0.25, 0.50, 1.0, and 2.0 million gallons per day) and withdrawal cycles (60, 180, and 365 days) was determined. For a given operational scheme, recovery efficiency increased as the injection and withdrawal rate is increased. Model results indicate that recovery efficiencies of between about 23 and 37 percent can be obtained for different subsurface storage and recovery schemes. Five successive injection, storage, and recovery cycles can increase the recovery efficiency to about 46 to 62 percent. There is a larger rate of increase at smaller rates than at larger rates. Over the range of variables studied, recovery efficiency improved with successive cycles, increasing rapidly during initial cycles tyhen more slowly at later cycles. The operation of a single well used for subsurface storage and recovery appears to be technically feasible under moderately favorable conditions; however, the recovery efficiency is higly dependent upon local physical and operational parameters. A combination of hydraulic, chemical, and operational parameters that minimize dispersion and buoyancy flow, maximizes recovery efficiency. Recovery efficiency was optimal where resident formation water density and permeabilities were relatively similar and low.

Water-Resources Investigations Report↗