Search USGS⌕ Search

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
USGS1967

Theory of friction based on brittle fracture

A theory of friction is presented that may be more applicable to geologic materials than the classic Bowden and Tabor theory. In the model, surfaces touch at the peaks of asperities and sliding occurs when the asperities fail by brittle fracture. The coefficient of friction, μ, was calculated from the strength of asperities of certain ideal shapes; for cone‐shaped asperities, μ is about 0.1 and for wedge‐shaped asperities, μ is about 0.15. For actual situations which seem close to the ideal model, observed μ was found to be very close to 0.1, even for materials such as quartz and calcite with widely differing strengths. If surface forces are present, the theory predicts that μ should decrease with load and that it should be higher in a vacuum than in air. In the presence of a fluid film between sliding surfaces, μ should depend on the area of the surfaces in contact. Both effects are observed. The character of wear particles produced during sliding and the way in which μ depends on normal load, roughness, and environment lend further support to the model of friction presented here.

Journal of Applied Physics↗
USGS2024

Laboratory hydrofractures as analogs to tectonic tremors

The fracture of Earth materials occurs over a wide range of time and length scales. Physical conditions, particularly the stress field and Earth material properties, may condition rupture in a specific fracture regime. In nature, fast and slow fractures occur concurrently: tectonic tremor events are fast enough to emit seismic waves and frequently accompany slow earthquakes, which are too slow to emit seismic waves and are referred to as aseismic slip events. In this study, we generate simultaneous seismic and aseismic processes in a laboratory setting by driving a penny-shaped crack in a transparent sample with pressurized fluid. We leverage synchronized high-speed imaging and high-frequency acoustic emission (AE) sensing to visualize and listen to the various sequences of propagation (breaks) and arrest (sticks) of a fracture undergoing stick-break instabilities. Slow radial crack propagation is facilitated by fast tangential fractures. Fluid viscosity and pressure regulate the fracture dynamics of slow and fast events, and control the inter-event time and the energy released during individual fast events. These AE signals share behaviors with observations of episodic tremors in Cascadia, United States; these include: (a) bursty or intermittent slow propagation, and (b) nearly linear scaling of radiated energy with area. Our laboratory experiments provide a plausible model of tectonic tremor as an indicative of hydraulic fracturing facilitating shear slip during slow earthquakes.

AGU Advances↗
USGS2000

Characterization methods for fractured glacial tills

This paper provides a literature review of methods successfully employed to characterize finegrained and fractured or unfractured glacial deposits. Descriptions and examples are given for four major categories of characterization methods: physical, hydraulic, chemical, and indirect. Characterization methods have evolved significantly within the past ten years; however, there still exists uncertainty about the reliability of individual characterization methods applied to till deposits. Therefore, a combination of methods is best, the choice of which depends on the objectives of the work. Sampling methods, sampling scales, and reporting methods are extremely important and should be considered when interpreting and comparing results between sites. Recognition of these issues is necessary to ensure that decisions regarding the transport of fluids in fractured tills are not based on the assumption that poorly permeable tills are always an inhibitor of subsurface flow.

Ohio Journal of Science↗
USGS2024

Triggering the 2022 eruption of Mauna Loa

Distinguishing periods of intermittent unrest from the run-up to eruption is a major challenge at volcanoes around the globe. Comparing multidisciplinary monitoring data with mineral chemistry that records the physical and spatio-temporal evolution of magmas fundamentally advances our ability to forecast eruptions. The recent eruption of Mauna Loa, Earth’s largest active volcano, provides a unique opportunity to differentiate unrest from run-up and improve forecasting of future eruptions. After decades of intermittent seismic and geodetic activity over 38 years of repose, Mauna Loa began erupting on 27 November 2022. Here we present a multidisciplinary synthesis that tracks the spatio-temporal evolution of precursory activity by integrating mineral and melt chemistry, fluid inclusion barometry, numerical modeling of mineral zoning, syn-eruptive gas plume measurements, the distribution and frequency of earthquake hypocenters, seismic velocity changes, and ground deformation. These diverse data indicate that the eruption occurred following a 2-month period of sustained magma intrusion from depths of 3–5 km up to 1–2 km beneath the summit caldera, providing a new model of the plumbing system at this very high threat volcano. Careful correlation of both the geochemistry and instrumental monitoring data improves our ability to distinguish unrest from the run-up to eruption by providing deeper understanding of the both the monitoring data and the magmatic system—an approach that could be applied at other volcanic systems worldwide.

Hawaii↗
USGS2014

Assessment of the geoavailability of trace elements from selected zinc minerals

This assessment focused on five zinc-bearing minerals. The minerals were subjected to a number of analyses including quantitative X-ray diffraction, optical microscopy, leaching tests, and bioaccessibility and toxicity studies. Like a previous comprehensive assessment of five copper-bearing minerals, the purpose of this assessment was to obtain structural and chemical information and to characterize the reactivity of each mineral to various simulated environmental and biological conditions. As in the copper minerals study, analyses were conducted consistent with widely accepted methods. Unless otherwise noted, analytical methods used for this study were identical to those described in the investigation of copper-bearing minerals. Two sphalerite specimens were included in the zinc-minerals set. One sphalerite was recovered from a mine in Balmat, New York; the second came from a mine in Creede, Colorado. The location and conditions of origin are significant because, as analyses confirmed, the two sphalerite specimens are quite different. For example, data acquired from a simulated gastric fluid (SGF) study indicate that the hydrothermally formed Creede sphalerite contains orders of magnitude higher arsenic, cadmium, manganese, and lead than the much older metamorphic Balmat sphalerite. The SGF and other experimental results contained in this report suggest that crystallizing conditions such as temperature, pressure, fluidization, or alteration processes significantly affect mineral properties—properties that, in turn, influence reactivity, solubility, and toxicity. The three remaining minerals analyzed for this report—smithsonite, hemimorphite, and hydrozincite—are all secondary minerals or alteration products of zinc-ore deposits. In addition, all share physical characteristics such as tenacity, density, streak, and cleavage. Similarities end there. The chemical composition, unit-cell parameters, acid-neutralizing potential, and other observable and quantifiable properties indicate very different minerals. Only one of each of these minerals was studied. Had this assessment included multiples of these minerals, geochemical and mineralogical distinctions would have emerged, similar to the results for the two sphalerite specimens.

Arizona;Chihuahua;Colorado;New York↗
USGS2024

Real-time lava flow forecasting during the 2022 Mauna Loa eruption response

On November 27, 2022, Mauna Loa (Hawai‘i) erupted for the first time in 38 years, initially producing lava flows that covered the floor of its summit caldera, Moku‘āweoweo. Over the first 12 h following the summit eruption, four main fissures opened on Mauna Loa’s Northeast Rift Zone, with “fissure 3” quickly becoming the dominant source of lava flows. For the next 12 days, fissure 3 produced a 19-km-long lava flow to the north, crossing the Mauna Loa Weather Observatory access road and coming within 2.8 km of inundating the Daniel K. Inouye Highway (Saddle Road). Within 40 min of fissure 3 opening, inundation modeling efforts had begun. For the duration of the eruption, the computational fluid dynamics model Lava2d was run in real time, using flow front locations and other field observations to make sequential forecast improvements, eventually producing a set of models which accurately predicted the routing and arrival times of lava from fissure 3. These models were used to inform timing estimates of possible future inundation of the Saddle Road. As the eruption progressed, almost 4000 Lava2d models were made using high-performance computing resources, providing critical information on uncertainty in multi-week forecasts. To our knowledge, this was the first ever real-time physics-based ensemble lava flow modeling and forecasting effort. Here, we present a chronology of these real-time efforts, focusing on the successes and limitations of this approach.

Hawaii↗
USGS2008

Scientific results from Gulf of Mexico Gas Hydrates Joint Industry Project Leg 1 drilling: Introduction and overview

The Gulf of Mexico Gas Hydrates Joint Industry Project (JIP) is a consortium of production and service companies and some government agencies formed to address the challenges that gas hydrates pose for deepwater exploration and production. In partnership with the U.S. Department of Energy and with scientific assistance from the U.S. Geological Survey and academic partners, the JIP has focused on studies to assess hazards associated with drilling the fine-grained, hydrate-bearing sediments that dominate much of the shallow subseafloor in the deepwater (>500 m) Gulf of Mexico. In preparation for an initial drilling, logging, and coring program, the JIP sponsored a multi-year research effort that included: (a) the development of borehole stability models for hydrate-bearing sediments; (b) exhaustive laboratory measurements of the physical properties of hydrate-bearing sediments; (c) refinement of new techniques for processing industry-standard 3-D seismic data to constrain gas hydrate saturations; and (d) construction of instrumentation to measure the physical properties of sediment cores that had never been removed from in situ hydrostatic pressure conditions. Following review of potential drilling sites, the JIP launched a 35-day expedition in Spring 2005 to acquire well logs and sediment cores at sites in Atwater Valley lease blocks 13/14 and Keathley Canyon lease block 151 in the northern Gulf of Mexico minibasin province. The Keathley Canyon site has a bottom simulating reflection at ???392 m below the seafloor, while the Atwater Valley location is characterized by seafloor mounds with an underlying upwarped seismic reflection consistent with upward fluid migration and possible shoaling of the base of the gas hydrate stability (BGHS). No gas hydrate was recovered at the drill sites, but logging data, and to some extent cores, suggest the occurrence of gas hydrate in inferred coarser-grained beds and fractures, particularly between 220 and 330 m below the seafloor at the Keathley Canyon site. This paper provides an overview of the results of the initial phases of the JIP work and introduces the 15 papers that make up this special volume on the scientific results related to the 2005 logging and drilling expedition.

Marine and Petroleum Geology↗
USGS2005

Performance of four turbulence closure models implemented using a generic length scale method

A two-equation turbulence model (one equation for turbulence kinetic energy and a second for a generic turbulence length-scale quantity) proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235] is implemented in a three-dimensional oceanographic model (Regional Oceanographic Modeling System; ROMS v2.0). These two equations, along with several stability functions, can represent many popular turbulence closures, including the k-kl (Mellor-Yamada Level 2.5), k-??, and k-?? schemes. The implementation adds flexibility to the model by providing an unprecedented range of turbulence closure selections in a single 3D oceanographic model and allows comparison and evaluation of turbulence models in an otherwise identical numerical environment. This also allows evaluation of the effect of turbulence models on other processes such as suspended-sediment distribution or ecological processes. Performance of the turbulence models and sediment-transport schemes is investigated with three test cases for (1) steady barotropic flow in a rectangular channel, (2) wind-induced surface mixed-layer deepening in a stratified fluid, and (3) oscillatory stratified pressure-gradient driven flow (estuarine circulation) in a rectangular channel. Results from k-??, k-??, and gen (a new closure proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235]) are very similar for these cases, but the k-kl closure results depend on a wall-proximity function that must be chosen to suit the flow. Greater variations appear in simulations of suspended-sediment concentrations than in salinity simulations because the transport of suspended-sediment amplifies minor variations in the methods. The amplification is caused by the added physics of a vertical settling rate, bottom stress dependent resuspension, and diffusive transport of sediment in regions of well mixed salt and temperature. Despite the amplified sensitivity of sediment to turbulence models in the estuary test case, the four closures investigated here all generated estuarine turbidity maxima that were similar in their shape, location, and concentrations.

Ocean Modelling↗
USGS2008

Comment on the treatment of residual water content in “A consistent set of parametric models for the two-phase flow of immiscible fluids in the subsurface” by L. Luckner et al.

Luckner et al. [1989] (hereinafter LVN) present a clear summary and generalization of popular formulations used for convenient representation of porous media fluid flow characteristics, including water content ( θ ) related to suction ( h ) and hydraulic conductivity ( K ) related to θ or h. One essential but problematic element in the LVN models is the concept of residual water content ( θ r ; in LVN, θ w , r ). Most studies using θ r determine its value as a fitted parameter and make the assumption that liquid flow processes are negligible at θ values less than θ r . While the LVN paper contributes a valuable discussion of the nature of θ r , it leaves several problems unresolved, including fundamental difficulties in associating a definite physical condition with θ r , practical inadequacies of the models at low θ values, and difficulties in designating a main wetting curve.

Water Resources Research↗
USGS2003

Pallid sturgeon in the Lower Mississippi Region: Hematology and genome information

This project (Project 1448-43270-2M-002) has been coordinated through the Natchitoches National Fish Hatchery (NNFH) and the U.S. Geological Survey’s National Wetlands Research Center (NWRC). From November 2001 to April 2002, over 280 sturgeon of the genus Scaphirhynchus (including pallid sturgeon, shovelnose, and their hybrids) were sampled from the outflow channel of the Old River Control Structure Complex (ORCC) in Concordia Parish, La. In the overall project, several datasets were collected (see Appendix), including species identification by using microsatellites and morphometric characters, food habits, physical anomalies, information on blood cells, and pathologic evidence of iridovirus – the first indication in the lower Mississippi population of pallid sturgeon. In this study, data on blood cells were obtained from the sturgeon collected monthly from approximately 20 different animals at each sampling time. This report presents preliminary information on differential blood cell identifications in sturgeon, data on comparative genomic DNA content and DNA degradation, and summaries and interpretations of data collected in light of available scientific literature addressing blood parameters of fish and sturgeon, in particular. Results obtained from collection and examination of blood and body fluids are often essential in establishing the health of fish (Blaxhall, 1972; Fange, 1992). Blood cells and sperm cells can be obtained nondestructively from fishes, even from small specimens that weigh less than 100 g (Stoskopf, 1992a). For flow cytometry assays, whereby cells are analyzed individually in a fluid stream, less than 1 :L of blood is needed. Examinations of blood by microscopy and flow cytometry were performed at NWRC in assisting in the efforts directed at recovery of the pallid sturgeon population in the Lower Mississippi River Basin.

Open-File Report↗
USGS2018

Induced earthquake families reveal distinctive evolutionary patterns near disposal wells

The timing of events in seismic sequences can provide insights into the physical processes controlling fault slip. In southern Kansas, the rate of earthquakes rose rapidly starting in 2013 following expansion of energy production into the area, demanding the disposal of large volumes of wastewater into deep wells. Seismicity catalogs that are complete to low magnitudes can provide insights into the physical processes that induce seismicity near wastewater disposal. We develop a catalog of over 130,000 earthquakes recorded in southern Kansas from mid‐March 2014 through December 2017 by applying a matched filter algorithm to an original catalog of 5,831 template earthquakes. Detections have nearly identical waveforms to their associated template event and represent slip on nearly co‐located sections of a fault. We select template events with at least 100 associated detections and examine the characteristics of these prolific families of earthquakes. We find that families located close (<10 km) to areas with significant volumes of injected fluids have near‐Poissonian interevent times and the families remain active over longer durations. Families farther from high‐volume injection wells show strong clustering of interevent times and shorter sequence durations. We conclude that increasing pore fluid pressures from nearby disposal of large volumes of wastewater is the primary driver of these long duration episodes, with earthquake‐earthquake interactions driving sequences at greater distance from the wells.

Kansas↗
USGS2005

Impacts of petroleum production on ground and surface waters: Results from the Osage-Skiatook Petroleum Environmental Research A site, Osage County Oklahoma

Kharaka has been a research hydrogeochemist with the U.S. Geological Survey since1975. His current research covers the fields of water-rock-gas interactions, CO 2 sequestration, contamination from agricultural drainage and petroleum produced water, and fluid-fault interactions. Kharaka has authored more than 100 scientific papers and book chapters and has delivered about 200 presentations. He received his Ph.D. from the University of California, Berkeley, in 1971. Thordsen is a geologist, employed with the U.S. Geological Survey since 1981 and presently working in the project team Hydrogeochemistry of Fine Grained Sediments, in the Water Resources Division, National Research Program. His recent research covers water-rock interactions and geochemical modeling in geothermal systems, fault-fluid systems, and petroleum produced water brine impacts. He received his M.S. in geology from the Ohio State University in 1988. Kakouros has been a hydrogeologist in the U.S. Geological Survey Water Resources Division in Menlo Park, California, since 1999. His current research covers the fields of Se contamination from agricultural drainage, CO 2 sequestration, and contamination from petroleum and petroleum-produced water. Kakouros received his M.Sc. from the San Jose State University in 2001. Herkelrath is a physicist in the U.S. Geological Survey's Water Resources Division National Research Program in Menlo Park, California. Herkelrath received his Ph.D. in physics from the University of Wisconsin at Madison. His recent research interests include multiphase flow and solute transport modeling in shallow ground-water systems where oil or brine is present.

Oklahoma↗
USGS2004

Lattice-Boltzmann simulation of coalescence-driven island coarsening

The first-order phase separation in a thin fluid film was simulated using a two-dimensional lattice-Boltzman model (LBM) with fluid-fluid interactions. The effects of the domain size on the intermediate asymptotic island size distribution were also discussed. It was observed that the overall process is dominated by coalescence which is independent of island mass. The results show that the combined effects of growth, coalescence, and Ostwald ripening control the phase transition process in the LBM simulations.

Journal of Chemical Physics↗
USGS2006

An effective medium inversion algorithm for gas hydrate quantification and its application to laboratory and borehole measurements of gas hydrate-bearing sediments

The presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments.

Geophysical Journal International↗
USGS1997

Practical use of video imagery in nearshore oceanographic field studies

An approach was developed for using video imagery to quantify, in terms of both spatial and temporal dimensions, a number of naturally occurring (nearshore) physical processes. The complete method is presented, including the derivation of the geometrical relationships relating image and ground coordinates, principles to be considered when working with video imagery and the two-step strategy for calibration of the camera model. The techniques are founded on the principles of photogrammetry, account for difficulties inherent in the use of video signals, and have been adapted to allow for flexibility of use in field studies. Examples from field experiments indicate that this approach is both accurate and applicable under the conditions typically experienced when sampling in coastal regions. Several applications of the camera model are discussed, including the measurement of nearshore fluid processes, sand bar length scales, foreshore topography, and drifter motions. Although we have applied this method to the measurement of nearshore processes and morphologic features, these same techniques are transferable to studies in other geophysical settings.

IEEE Journal of Oceanic Engineering↗
USGS2000

Debris flow monitoring in the Acquabona watershed on the Dolomites (Italian Alps)

In 1997 a field monitoring system was installed in Acquabona Creek in the Dolomites (Eastern Italian Alps) to observe the hydrologic conditions for debris flow occurrence and some dynamic properties of debris flow. The monitoring system consists of three remote stations: an upper one located at the head of a deeply-incised channel and two others located downstream. The system is equipped with sensors for measuring rainfall, pore pressures in the mobile channel bottom, ground vibrations, debris flow depth, total normal stress and fluid pore-pressure at the base of the flow. Two video cameras record events at the upper channel station and one video is installed at the lowermost station. During summer 1998, three debris flows (volumes from less than 1000 m 3 up to 9000 m 3 ) occurred at Acquabona. The following results were obtained from a preliminary analysis of the data: 1) All of the flows were triggered by rainfalls of less than 1 hour duration, with peak rainfall intensities ranging from 4.8 to 14.7 mm / 10 minute. 2) Debris flows initiated in several reaches of the channel, including the head of the talus slope. 3) The initial surges of the mature flows had a higher solid concentration and a lower velocity (up to 4 m/s) than succeeding, more dilute surges (more than 7 m/s). 4) Total normal stress and pore fluid pressures measured at the base of the flow. (mean depth about 1.1 m) were similar (about 15 kPa), indicating a completely liquefied flow. 5) Peak flows entrained debris at a rate of about 6 m 3/m of channel length and channel bed scouring was proportional to the local slope gradient and was still evident in the lower channel where the slope was 7°.

Acquabona Creek, Alps, Boite River, Dolomites↗
USGS2004

Methane hydrate formation in partially water-saturated Ottawa sand

Bulk properties of gas hydrate-bearing sediment strongly depend on whether hydrate forms primarily in the pore fluid, becomes a load-bearing member of the sediment matrix, or cements sediment grains. Our compressional wave speed measurements through partially water-saturated, methane hydrate-bearing Ottawa sands suggest hydrate surrounds and cements sediment grains. The three Ottawa sand packs tested in the Gas Hydrate And Sediment Test Laboratory Instrument (GHASTLI) contain 38(1)% porosity, initially with distilled water saturating 58, 31, and 16% of that pore space, respectively. From the volume of methane gas produced during hydrate dissociation, we calculated the hydrate concentration in the pore space to be 70, 37, and 20% respectively. Based on these hydrate concentrations and our measured compressional wave speeds, we used a rock physics model to differentiate between potential pore-space hydrate distributions. Model results suggest methane hydrate cements unconsolidated sediment when forming in systems containing an abundant gas phase.

American Mineralogist↗
USGS1998

Seasonal uranium distributions in the coastal waters off the Amazon and Mississippi Rivers

The chemical reactivity of uranium was investigated across estuarine gradients from two of the world’s largest river systems: the Amazon and Mississippi. Concentrations of dissolved (<0.45 μm) uranium (U) were measured in surface waters of the Amazon shelf during rising (March 1990), flood (June 1990) and low (November 1991) discharge regimes. The dissolved U content was also examined in surface waters collected across estuarine gradients of the Mississippi outflow region during April 1992, August 1993, and November (1993). All water samples were analyzed for U by isotope dilution inductively coupled plasma mass spectrometry (ICP-MS). In Amazon shelf surface waters uranium increased nonconservatively from about 0.01 μg I −1 at the river’s mouth to over 3 μg I −1 at the distal site, irrespective of river discharge stage. Observed large-scale U removal at salinities generally less than 15 implies a) that riverine dissolved U was extensively adsorbed by freshly-precipitated hydrous metal oxides (e.g., FeOOH, MnO 2 ) as a result of flocculation and aggregation, and b) that energetic resuspension and reworking of shelf sediments and fluid muds on the Amazon shelf released a chemically reactive particle/colloid to the water column which can further scavenge dissolved U across much of the estuarine gradient. In contrast, the estuarine chemistry of U is inconclusive within surface waters of the Mississippi shelf-break region. U behavior is most likely controlled less by traditional sorption and/or desorption reactions involving metal oxides or colloids than by the river’s variable discharge regime (e.g., water parcel residence time during estuarine mixing, nature of particulates, sediment storage and resuspension in, the confined lower river), and plume dispersal. Mixing of the thin freshwater lens into ambient seawater is largely defined by wind-driven rather than physical processes. As a consequence, in the Mississippi outflow region uranium predominantly displays conservative behavior; removal is evident only during anomalous river discharge regimes. ‘Products-approach’ mixing experiments conducted during the Flood of 1993 suggest the importance of small particles and/or colloids in defining a depleted U versus salinity distribution.

Louisiana↗