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 181 records · Page 10Linked to original sources

USGS2011

Detecting aseismic strain transients from seismicity data

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate‐ and state‐dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space‐time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M ≥ 1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Science mission requirements for a globally ranging, riserless drilling vessel for U.S. Scientific Ocean Drilling

Through the collection and analysis of shallow and deep subseafloor sediments, rocks, fluids, and life, scientific ocean drilling has enriched our understanding of the complex Earth system. Among other achievements, scientific ocean drilling has documented the history of Earth’s climate, the waxing and waning of polar ice sheets, the past changes in ocean and atmospheric circulation, the existence and function of microbial life in the subseafloor, the compositional variations in Earth’s crust and underlying mantle, and the physical and chemical processes acting at subduction zones, including those associated with tsunamigenic earthquakes. Over the decades, more than 12,000 articles that depend on analyses of scientific ocean drilling samples and geophysical data have been published, many detailing breakthrough contributions to global knowledge about the Earth system. Approximately 45% of these publications were led by U.S.-affiliated authors (International Ocean Discovery Program Publication Services, 2021). Since the mid-1980s, the workhorse of this multidisciplinary, international research effort has been the riserless D/V JOIDES Resolution, operated by Texas A&M University with funding from the U.S. National Science Foundation (NSF). D/V JOIDES Resolution has conducted the vast majority of scientific ocean drilling expeditions and collected most of the scientific cores over that period, including 82% of the expeditions and 93% of the cores in the last decade alone, despite being one of three platforms that is operated within the International Ocean Discovery Program. However, D/V JOIDES Resolution is approaching the end of its useful life. With a strong commitment to continue scientific ocean drilling beyond the end of the current phase, the community developed a document outlining the research frontiers that should be pursued. Exploring Earth by Scientific Ocean Drilling: 2050 Science Framework (Koppers and Coggon, 2020) describes seven scientific strategic objectives that focus on understanding interconnections within the Earth system and five flagship initiatives that integrate these objectives into long-term research efforts that address issues facing society. Additional elements in the 2050 Science Framework, including STEM education, workforce development, technology development, and innovative applications of data analytics, will advance the goals of scientific ocean drilling. Addressing the 2050 Science Framework also requires building partnerships with allied U.S. and international science programs and strengthening existing ones. To implement a significant portion of the 2050 Science Framework, the U.S. scientific community seeks to lease or acquire a newly built, globally ranging, state-of-the art, riserless drilling vessel. The many and varied technical and human resources requirements for successful accomplishment of scientific and educational goals summarized in this document and described in detail in the 2050 Science Framework require broad community input and careful consideration. Following receipt of NSF’s formal Request for Assistance to the United States Science Support Program (USSSP), the U.S. scientific ocean drilling community conducted a one-year exercise to identify its national scientific needs and priorities in order to determine the Science Mission Requirements (SMRs) presented here. This community effort included: (1) a U.S. community-wide survey to identify the specific operational and technical capabilities critical to addressing science in the 2050 Science Framework; (2) a series of online workshops focusing on critical capabilities identified by the survey; and (3) a large in-person workshop to synthesize the results of the survey and the virtual workshops (Appendix 1). The approach was designed to reach as many participants as possible. Overall, 278 survey responses were received from U.S. community members, representing 104 unique institutions from 39 states and the District of Columbia, and 137 unique individuals participated in the workshops (Appendix 2). The results of this effort comprise two classes of SMRs: Foundational Science Mission Requirements and Primary Science Mission Requirements. Foundational SMRs define minimum criteria for a new riserless drilling vessel that can address significant portions of the 2050 Science Framework. Primary SMRs build upon the Foundational SMRs and will create more robust science opportunities and data collection capabilities, will increase progress in addressing the 2050 Science Framework objectives, and will provide more real-time ship-to-shore interaction to improve science productivity, engagement, and outreach. Modern safety and environmental standards, including meeting standards to access protected waters such as exclusive economic zones, extended continental shelves, or high latitudes, while being cognizant of the vessel’s environmental footprint. Safe and efficient operations in global locations and in water depths from 70 m to 6000+ m, with total drill string length of at least 7000 m. High-quality core and data collection from a range of key subseafloor environments. Advanced heave compensation, dynamic positioning, and drill pipe stability. Modern mud and cement/casing systems. Critical onboard measurements for safety, operational decision-making, documentation of ephemeral properties, mission-specific science, and long-term science goals that extend beyond a single expedition. Designated and appropriate space for sample and data preservation. Highly skilled onboard personnel, including technical staff for curation and core handling; support for safety, time-sensitive, and critical shipboard measurements; computer support; equipment and instrument repair; application support; and data assurance. Primary Science Mission Requirements include: Flexible shipboard space for laboratories and on deck to ensure safe, successful implementation of diverse science objectives and operations. Minimizing contamination of recovered samples. Over-the-side capabilities for science-supporting technology (e.g., remotely operated vehicles, water- column sampling, sediment-water interface sampling). Downhole logging and measurements. Expanded borehole observatory capabilities. Reliable and consistent ship-to-shore communications. NSF’s investment in a new globally ranging, riserless drilling vessel will have a powerful economic multiplier effect, including the infusion of additional science support funds in the United States for training and research, the development of new technologies and tools, and the associated scientific and technical workforce development. The skills and knowledge gained through scientific ocean drilling are translatable to careers in fields such as sustainable energy development (e.g., geothermal and offshore wind), carbon sequestration, data management and cyberinfrastructure, biotechnology, communications, science education, policy, hazard mitigation, and environmental management. The United States is a leader in a well-established and internationally collaborative scientific ocean drilling community. A modern, globally ranging, riserless drilling vessel will allow the United States to expand its leadership position, address broad scientific questions that current capabilities preclude, and cultivate equitable international, multidisciplinary collaborations that will ensure scientific ocean drilling’s future success.

Report↗

USGS1959

Lithologic log and drilling information for the Marble exploration hole 3, U15 area, Nevada Test Site, Nye County, Nevada

The Marble exploration hole 3 is in the northern part of the Tippipah Spring NE 7-1/2 minute quadrangle at the north end of Yucca Valley (fig.1). The collar elevation of the hole is about 5,316 feet and has the approximate Nevada State coordinates of N. 903,093 and E. 674,870. The hole was drilled to determine the physical and chemical properties of the marble within a radius of 200 feet of a point 950 feet below the surface. The exploration hole is 978 feet deep and was drilled from June 12 to July 8, 1959, by Minerals Engineering Co. of Grand Junction, Colorado, using conventional NX coring equiipment and a Portadrill 2,500 drill rig. Water-base bentonitic mud was used as the drilling fluid in the upper part of the hole. Clear water was used after circulation became unrestorable below 200 feet. The hole was cased from surface to 20 feet. The drilling rate during coring ranged from 2.0 to 19.3 minutes per foot and averaged 4.6 minutes per foot (table 1). The core recovered was excellent in quality and quantity. The average core recovery is 92.7 percent (table 2).

Nevada↗

USGS2024

Relatively stable pressure effects and time-increasing thermal contraction control Heber geothermal field deformation

Due to geological complexities and observational gaps, it is challenging to identify the governing physical processes of geothermal field deformation including ground subsidence and earthquakes. In the west and east regions of the Heber Geothermal Field (HGF), decade-long subsidence was occurring despite injection of heat-depleted brines, along with transient reversals between uplift and subsidence. These observed phenomena contradict current knowledge that injection leads to surface uplift. Here we show that high-yield production wells at the HGF center siphon fluid from surrounding regions, which can cause subsidence at low-rate injection locations. Moreover, the thermal contraction effect by cooling increases with time and eventually overwhelms the pressure effects of pressure fluctuation and poroelastic responses, which keep relatively stable during geothermal operations. The observed subsidence anomalies result from the siphoning effect and thermal contraction. We further demonstrate that thermal contraction dominates long-term trends of surface displacement and seismicity growth, while pressure effects drive near-instantaneous changes.

California↗

USGS1988

Facies composition calculated from the sonic, neutron, and density log suite, upper part of the Minnelusa Formation, Powder River basin, Wyoming

Sandstones and dolomites of the Permian upper part of the Minnelusa Formation are treated here as four-component systems consisting of fluid-filled pore space, quartz, dolomite, and anhydrite. Response equations of sonic, neutron, and density logs form a system of four simultaneous equations. With four equations and four unknowns, the composition of upper Minnelusa facies is defined by the three-log suite and can be calculated by solving a 4 ?? 4 matrix. Such calculations of facies composition help in establishing subsurface correlations and yield information on the diagenesis and physical character of upper Minnelusa sandstones and dolomites. Applications of composition calculations are illustrated by examples drawn from the area of the West Mellott field (T52N, R68W), where the upper Minnelusa is at depths of about 7000 ft (2100m). -from Authors

Mountain Geologist↗

USGS2023

Simulation of regional groundwater flow and advective transport of per- and polyfluoroalkyl substances, Joint Base McGuire-Dix-Lakehurst and vicinity, New Jersey, 2018

A three-dimensional numerical model of groundwater flow was developed and calibrated for the unconsolidated New Jersey Coastal Plain aquifers underlying Joint Base McGuire-Dix-Lakehurst (JBMDL) and vicinity, New Jersey, to evaluate groundwater flow pathways of per- and polyfluoroalkyl substances (PFAS) contamination associated with use of aqueous film forming foam (AFFF) at the base. The regional subsurface flow model spans an area of approximately 518 square miles around JBMDL and is based on a previously developed hydrogeologic framework of the area. Steady-state flow in the unconsolidated aquifers was simulated using the MODFLOW 6 groundwater flow model, which is able to account for hydrostratigraphic pinchouts and discontinuities in the Coastal Plain aquifers underlying JBMDL. To account for local patterns of fluid flow driving advective subsurface migration of PFAS, the grid was refined using quadtree meshes spanning 21 areas where historical AFFF use was identified, five off-site reconnaissance areas identified by AFCEC as areas in which the occurrence of PFAS is most likely to pose a potential danger to local drinking water supplies, and along streams that behave as drains in the base-flow-dominated Coastal Plain. Following grid refinement, four physical processes known to govern subsurface flow were introduced to the model. These included effective precipitation recharge, discharge to streams and stream-connected wetlands, regional inflows and outflows along the model bottom, and withdrawals from wells, each of which were incorporated into the model as either external or internal boundary conditions. To account for effective precipitation recharge, a specified-flow boundary was assigned along the top of the model. Similarly, regional flows predicted using the modified U.S Geological Survey’s New Jersey Coastal Plain Regional Aquifer System Analysis model were treated as specified-flow boundary conditions along the bottom of the model. Base-flow losses were treated as drains along streams delineated using a 10-foot LiDAR dataset. Drains were also assigned to cells falling within stream-connected National Hydrologic Database wetlands. Finally, well-pumpage data mined from the New Jersey Water Transfer database were added to the model to account for extraction of groundwater through pumping from industrial-supply and drinking-water-supply wells. Along model edges established at groundwater divides, where the net flux of water across the boundary is equal to zero, natural no-flow boundary conditions were imposed. The refined flow model was calibrated using the parameter-estimation (PEST) program, which adjusts model parameters by performing a gradient search over the sum-of-squared-error objective function until the parameter set that produces simulated water levels and base flows most closely matches 544 water levels and 20 estimated base flows and closely adheres to initial parameter estimates. Based on the analysis of calibration residuals, the model did not appear to be affected by significant model structural error. The MODPATH particle-tracking algorithm was used to estimate advective transport paths of PFAS in the vicinity of JBMDL. Forward tracking was used to determine paths of PFAS away from AFFF source areas to streams, wetlands, pumping wells, and geographic areas that PFAS may contaminate. Additionally, reverse tracking was used to determine particle pathlines away from off-site PFAS reconnaissance areas, or areas within which all sources of PFAS might be advectively transported into subsurface drinking-water supplies, to locations at land surface that may indicate a source of PFAS. The coupled and calibrated groundwater flow and particle-tracking transport model provide valuable tools for predicting the relative extent of PFAS contamination from onsite legacy source areas. The calibrated model also provides measures of water-level and base-flow observation influence that can help guide future data-collection efforts related to groundwater and surface water sampling for PFAS.

New Jersey↗

USGS2003

Modeling radium and radon transport through soil and vegetation

A one-dimensional flow and transport model was developed to describe the movement of two fluid phases, gas and water, within a porous medium and the transport of 226Ra and 222Rn within and between these two phases. Included in this model is the vegetative uptake of water and aqueous 226Ra and 222Rn that can be extracted from the soil via the transpiration stream. The mathematical model is formulated through a set of phase balance equations and a set of species balance equations. Mass exchange, sink terms and the dependence of physical properties upon phase composition couple the two sets of equations. Numerical solution of each set, with iteration between the sets, is carried out leading to a set-iterative compositional model. The Petrov-Galerkin finite element approach is used to allow for upstream weighting if required for a given simulation. Mass lumping improves solution convergence and stability behavior. The resulting numerical model was applied to four problems and was found to produce accurate, mass conservative solutions when compared to published experimental and numerical results and theoretical column experiments. Preliminary results suggest that the model can be used as an investigative tool to determine the feasibility of phytoremediating radium and radon-contaminated soil. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Contaminant Hydrology↗

USGS2012

Isotropic, anisotropic, and borehole washout analyses in Gulf of Mexico Gas Hydrate Joint Industry Project Leg II, Alaminos Canyon well 21-A

Through the use of three-dimensional seismic amplitude mapping, several gas hydrate prospects were identified in the Alaminos Canyon area of the Gulf of Mexico. Two of the prospects were drilled as part of the Gulf of Mexico Gas Hydrate Joint Industry Program Leg II in May 2009, and a suite of logging-while-drilling logs was acquired at each well site. Logging-while-drilling logs at the Alaminos Canyon 21–A site indicate that resistivities of approximately 2 ohm-meter and P-wave velocities of approximately 1.9 kilometers per second were measured in a possible gas-hydrate-bearing target sand interval between 540 and 632 feet below the sea floor. These values are slightly elevated relative to those measured in the hydrate-free sediment surrounding the sands. The initial well log analysis is inconclusive in determining the presence of gas hydrate in the logged sand interval, mainly because large washouts in the target interval degraded well log measurements. To assess gas-hydrate saturations, a method of compensating for the effect of washouts on the resistivity and acoustic velocities is required. To meet this need, a method is presented that models the washed-out portion of the borehole as a vertical layer filled with seawater (drilling fluid). Owing to the anisotropic nature of this geometry, the apparent anisotropic resistivities and velocities caused by the vertical layer are used to correct measured log values. By incorporating the conventional marine seismic data into the well log analysis of the washout-corrected well logs, the gas-hydrate saturation at well site AC21–A was estimated to be in the range of 13 percent. Because gas hydrates in the vertical fractures were observed, anisotropic rock physics models were also applied to estimate gas-hydrate saturations.

Scientific Investigations Report↗

USGS2021

Correlation of porosity variations and rheological transitions on the southern Cascadia megathrust

The unknown onshore extent of megathrust earthquake rupture in the Cascadia subduction zone represents a key uncertainty in earthquake hazard for the Pacific Northwest that is governed by the physical state and mechanical properties of the plate interface. The Cascadia plate interface is segmented into an interseismically locked zone located primarily offshore that is expected to rupture in large earthquakes, a region of aseismic slow slip at greater depth, and an intervening transition zone of uncertain rupture potential. Here we image the evolution of the ratio of seismic compressional to shear wave velocities from the locked zone to the transition zone, which is related to changes in fluid content of the plate boundary zone, using a dense onshore–offshore seismic dataset from southernmost Cascadia. The locked zone shows evidence of high fluid content implying a high porosity, yet the downdip transition zone shows an order of magnitude lower porosity. This strong variation is consistent with models that contain a ductile region between the earthquake rupture and slow slip zones that would inhibit onshore propagation of future large earthquake ruptures and hence reduce seismic hazard.

California↗

USGS1971

Sulfur isotope study of the "B" limestone and Galena fissure ore deposits of the U.S. mine, Bingham mining district, Utah

The delta S 34 values for 261 sulfide concentrates from the "B" limestone replacement and Galena fissure deposits of the U. S. mine range from +3.8 to --8.0 permil; the mean is +0.56 permil. The data indicate that (1) sulfides of the "B" limestone ore deposits are depleted in S 34 by 0.4 to 0.7 permil relative to those of the Galena fissure ore deposits; (2) sulfides of both deposits exhibit a primary fractionation trend characterized by increasing S 34 depletion in the sequence pyrite, sphalerite, tetrahedrite, and galena; (3) delta values for coexisting sulfide pairs of the "B" limestone deposit increase with increasing altitude; (4) sulfides of both deposits are increasingly depleted in S 34 with increasing altitude; (5) the paragenetic order of sulfide deposition is accompanied by an apparent depletion of S 34 that is independent of mineral species; and (6) pyrite disseminated in wall rocks is depleted in S 34 relative to pyrite from adjacent ore deposits. The similarity of delta S 34 values and isotopic trends to those of the central porphyry-type deposit at Bingham (Field, 1966) suggests that the hydrothermal deposits of the district were formed by fluids that originated from a common source during one episode of sequential mineralization. Gross isotopic differences between the deposits, consistent differences between coexisting sulfides, and systematic variations of delta values are attributed to temperature effects imposed upon sulfide exchange equilibria; isotopic temperatures suggest a lateral gradient between the "B" limestone and adjacent Galena fissure deposits and a possible vertical gradient within the "B" limestone deposit. The vertical zonation of delta S 34 values and trends related to pyrite physical occurrence and sulfide paragenesis are provisionally attributed to a pH-Eh control. Although the apparent trends are imperfectly resolved because of textural complexities, the marked depletion of S 34 in pyrite disseminated in country rocks warrants additional study as a potential exploration method.

Utah↗

USGS2016

Dynamics of permeability evolution in stimulated geothermal reservoirs

Spatially and temporally evolving permeability fields are fundamentally associated with the operation of enhanced geothermal systems. This chapter explores the resulting magnitude and patterns of permeability alteration as well as the coupled physical processes that control the evolution of permeability during shear stimulation and long-term evolution of a geothermal reservoir. It also explores mechanisms involved in the enhancement of fracture permeability, targeting injection pressures that encourage self-propping shear stimulation alongside elastic pressure and thermal dilatation of fracture sets. Simulation of a short-term, high-injection-rate enhanced geothermal system stimulation explored the importance of hydraulic versus thermal effects and the impact of a temperature differential between fluid in the fractures and the adjoining rock matrix. The effects of variations in fracture spacing on the evolution of permeability is at least as important as the magnitude of initial permeability, and quantifying the driving mechanisms is critical to understanding reservoir response.

Book chapter↗

USGS1975

Flows of impact melt at lunar crater

Lavalike materials that were emplaced in a fluid state occur in and around lunar impact craters whose diameters range from 3 km to more than 200 km and whose ages span a time interval of at least 3.5 b.y. Evidence of fluid emplacement includes flow lobes and leveed channels, a veneer mantling rough crater topography, level-surfaced pools, and complex contraction fissuring. The distribution of the lavalike deposits conforms to asymmetries of other ejecta from the same craters, and the material is concentrated downrange to distances as great as a crater radius. The character and distribution of the lavalike materials support the idea that they formed by impact melting rather than by volcanism. Returned samples indicate that materials with appropriate physical characteristics are generated by partial melting of feldspathic rocks by impact. The geologic evidence at lunar craters suggests that there is more melt rock in and near the craters than is predicted by experiment and theory.

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

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↗