SEARCH · USGS Science

Results for “Geologic Investigation”

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.

922 records · Page 34Linked to original sources

USGS2025

VIPER site analysis

We needed to evaluate available orbital data of NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission area in order to derive a variety of maps to help the science team identify scientifically interesting places for the rover to visit and to provide scientific context for our mission. Some of these maps also fulfilled engineering and mission design needs to enable safe and efficient landing and roving. We incorporated data from the Lunar Reconnaissance Orbiter Camera, the Lunar Orbital Laser Altimeter, the Mini-RF instrument, the Chandrayaan-2 Orbital High Resolution Camera, the Korean Pathfinder Lunar Orbiter’s Shadowcam, the Kaguya Spectral Profiler and Multiband Imager, and the Chandrayaan-1 Moon Mineralogy Mapper. We used a variety of techniques to build these maps, including stereogrammetry, shape-from-shading, ice stability depth and surface temperature calculations, and the horizon method for solar illumination and direct-to-Earth communications maps. Altogether, these maps allowed us to survey for boulders, evaluate features in permanently shadowed regions that VIPER might explore, provide mineralogic context for what VIPER’s instruments may learn, estimate the ages and radar properties of craters in the VIPER mission area, and evaluate the potential for gravity traverses with the rover. These data and techniques provided a rich set of information from which both the VIPER science team and engineering teams were able to draw in order to plan a safe landing and to plan a VIPER surface mission that will be both scientifically valuable and robust from an operational perspective.

The Planetary Science Journal↗

USGS2026

Geochemical, mineralogical, and isotopic evidence for multi-stage genesis of the Hicks Dome REE + Y-HFSE-fluorite deposit, Illinois, USA

Hicks Dome hosts breccias enriched in rare earth elements (REE), Y, Th, F, Ba, Ti, Nb, and Be, alongside spatially associated lamprophyre dikes (ca. 271 Ma). Hicks Dome is located within the Illinois–Kentucky Fluorspar District, which hosts fluorite, Pb–Zn, and barite resources. This study investigates the genetic relationships between Hicks Dome mineralization in breccias, alkaline magmatism, and Illinois–Kentucky Fluorspar District mineralization. Lamprophyre dikes are light REE–enriched with chondrite-normalized abundances decreasing from La to Lu. The Host Breccia exhibits middle and heavy REE–enriched patterns that mirror those of the principal REE–Th host minerals, including fluorapatite, xenotime, and thorite. Textural evidence suggests recrystallization of phosphates, sulfates, and Ti–Nb oxides in the Host Breccia. U–Pb geochronology constrains multiple mineralizing events, with ages of 277 ± 18 Ma from low-Th apatite interpreted as main-stage mineralization, and 121.6 ± 9.7 Ma from high-Th apatite indicating later overprinting. O–H–C stable isotope data provide evidence for multiple stages of fluid-rock interaction and fluid mixing: (1) early magmatic fluids dissolved limestone country rock, (2) mixing between magmatic fluids and basinal brines led to main-stage mineralization in the Host Breccia, and (3) late-stage mineralization occurred following mixing of meteoric water and basinal brine. These results indicate that heavy REEs, high field strength elements, and fluorine precipitated proximal to its alkaline magmatic source because of fluid–rock interactions and fluid mixing. Subsequent fluid mixing drove late-stage recrystallization and additional fluorite formation, a process that may be similar to mineralization in the Illinois-Kentucky Fluorspar District.

Illinois↗

USGS2025

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington↗

USGS2024

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington↗

USGS2025

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗

USGS2025

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts↗

USGS2024

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota↗

USGS1976

Water-quality study of a reach of the Merced River in Yosemite National Park and vicinity, California, April 1973 through September 1974

A 25-mile (40-kilometre) reach of the Merced River was studied to provide the National Park Service with water-quality data prior to the operation of a new sewage-treatment plant. Data were collected periodically at four Merced River and two inflow stations. Discharge varied seasonally with highest flow occurring at the time of spring snowmelt runoff. Lowest flow occurred during late summer, coinciding with high park visitor use and elevated water temperature. Water-quality variables investigated were discharge, temperature, major chemical constituents, selected trace metals, specific conductance, suspended sediment, turbidity, dissolved oxygen, pH, alkalinity, carbon dioxide, total organic carbon, nitrogen, phosphorus, algal growth potential, coliform bacteria, periphyton, and benthic invertebrates. Diel measurements were made at some stations to delineate daily fluctuations in selected water-quality variables. Results indicate that water in the reach of the Merced River studied was of good quality. The river had low concentrations of dissolved solids (8 to 41 milligrams per litre) and low alkalinity values (2 to 24 milligrams per litre of alkalinity as CaCO 3 ). Mean concentrations of nitrogen and phosphorus were highest in the Merced River at Rancheria Flat. Algal growth potential tests showed that phosphorus was limiting algal growth at Happy Isles Bridge, El Capitan Bridge, and South Fork Merced River, whereas nitrogen was limiting growth at Rancheria Flat. Input of excessive nitrogen and phosphorus concentrations in the Merced River during the summer low-flow period, and mixing with warm tributary inflow could cause excessive algal growth and production. The algal types found were diatoms, green algae, and blue-green algae. Biomass/chlorophyll a ratios in the Merced River were low (0.3 to 7.7). Benthic invertebrates consisted chiefly of Diptera, Ephemeroptera, and Plecoptera.

California↗

USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗

USGS2026

Geochemical geodatabase of sedimentary strata (coal, coal-adjacent rocks, tuffaceous oil shale, phosphate-rich rocks) and produced water in the Uinta region, Utah and Colorado

The Geochemical Geodatabase of Sedimentary Strata (Coal, Coal-adjacent Rocks, Tuffaceous Oil Shale, Phosphate-rich Rocks) and Produced Water in the Uinta Region, Utah and Colorado, consists of compiled datasets acquired as part of the Carbon Ore, Rare Earth, and Critical Mineral (CORE-CM) Uinta Region assessment funded by the U.S. Department of Energy (DEFE0032046, 2021–2024; Birgenheier et al., 2024). The CORE-CM assessment focused on providing comprehensive geological and geochemical characterization of current and prospective sedimentary-hosted resources including coal, oil shale, phosphatic limestone, and produced water from oil and gas targets present in eastern Utah and northwestern Colorado (Figure 1). This Data Series includes a geodatabase that consists of analytical geochemical data collected September 2021 through December 2024 via portable X-ray fluorescence (pXRF), and laboratory measured analyses produced by inductively coupled plasma mass spectrometry (ICP-MS) and inductively coupled plasma optical emission spectroscopy (ICP-OES). The coal-related geochemical data are derived primarily from the Cretaceous Blackhawk Formation and Ferron Sandstone of Utah, and the Mesaverde Group of Colorado. Additional non-coal resources assessed include oil shale-bearing strata of the Eocene upper Green River Formation (Utah and Colorado), phosphate-rich limestone of the Permian Park City Formation (Utah) and produced water from oil and gas-bearing strata of the Eocene Green River and Wasatch Formations (Uinta Basin) and the Pennsylvanian Paradox Formation (Paradox Basin) (Table 1). The CORE-CM assessment included a wide range of lithologies present in the coal, oil shale, and phosphate geologic resource systems whether or not the specific lithology has current economic value. Geochemical analyses of produced water from oil and gas wells focused on current and emerging hydrocarbon targets in the central Uinta Basin and northern Paradox Basin. A total of 13,092 geochemical analyses from these geologic systems is provided in the included geodatabase. A series of coal quality data (e.g., composition and maceral analyses) is also included in the database and was digitized from archived coal samples from the Utah Geological Survey (Appendix A).

Colorado, Utah↗

USGS2026

Assessing environmental drivers of denitrification in restored riverine floodplains

Restoration of impaired floodplains is an increasingly prevalent strategy for alleviating water quality concerns and reducing downstream flooding at watershed scales. Floodplains temporarily store water and slow flow velocity to promote sedimentation during overbank flooding and remove inorganic nitrogen from floodwater and groundwater via denitrification. Evaluating the impacts of different restoration strategies on denitrification can inform more strategic investments into floodplain modifications that improve water quality outcomes. Our research investigates how denitrification rates in floodplains respond to environmental factors that are actionable from an engineering perspective through design and water resources management. We seasonally measured soil denitrification enzyme activity and various environmental characteristics in 4 floodplains with different restoration design and management approaches at the confluence of the Wabash and Tippecanoe Rivers in Indiana, United States. Our results showed that denitrification rates in an agricultural floodplain were significantly lower than in restored floodplains with native vegetation. Certain soil conditions characteristic of floodplain wetlands were associated with higher denitrification, particularly elevated total nitrogen, moisture, silt, and organic matter contents. Vegetation species composition was correlated with denitrification rates. This link may reflect the direct effects of vegetation on soil conditions, such as supplying labile organic carbon, or indirect effects, such as vegetation acting as an indicator of hydrologic regime and land use. Denitrification seasonally varied, peaking in winter when nitrate supply from rivers draining agricultural watersheds in the region is also high. Substrate limitation of soil denitrification enzyme activity was most significant during the summer when overbank flooding, which replenishes soil nitrogen stocks, rarely occurs. Our findings indicate that denitrification capacity will likely be maximized in riverine floodplains that are restored as wetlands with diverse native vegetation and enhanced hydrologic connectivity. Such restoration activities promote higher denitrification rates via elevated moisture, fine sediment deposition, and soil organic matter.

Tippecanoe River, Wabash River↗

USGS2024

A seismic nodal deployment to understand magmatic structure in the vicinity of the Pahala earthquake swarm

In summer-fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawaiʻi, with the goal of improving seismicity catalogs, and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at depths of ~ 25 – 40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the dataset for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short period instruments, we find that the nodes produce high quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data is likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawaiʻi. While this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawaiʻi continue to evolve in the future.

Hawaii↗

USGS2024

Preliminary ground and airborne-based geophysical mapping and modelling of an active hydrothermal system at Mammoth Lakes, California

Mammoth Lakes, California hosts a productive hydrothermal system within the seismically active south moat of Long Valley Caldera. Surficial evidence of the shallow hydrothermal system includes discrete zones of tree-kill dispersed between Shady Rest Park and the Casa Diablo Geothermal Power Plant (40 MW), as well as east of the power plant. The tree-kill areas are associated with elevated diffuse CO₂ emissions, heated ground, hydrothermal alteration, diffuse soil H₂S emissions, and gas vents. Previous mapping delineates prominent north and northwest trending structures within the south moat along the southwestern edge of the resurgent dome that may accommodate gas and fluid flow at the Shady Rest Park and Basalt Canyon Tree Kill Areas (SRTKA and BCTKA, respectively). Both tree-kill areas are also located along contacts between resurgent rhyolite, mafic lavas, and surficial deposits which may provide additional pathways for gas and fluid migration in the shallow subsurface. Characterizing structure and lithology using geophysical anomalies is critical to determining primary structural controls on the hydrothermal system and the extent of subsurface alteration at these sites. We conducted ground and airborne-based potential field geophysical surveys to map gravity and magnetic anomalies. These anomalies are then used to model subsurface geology, structure, and hydrothermal alteration. Here we present our preliminary geophysical mapping and modelling results at both tree-kill locations. Gravity and magnetic data suggest complex structural intersections are coincident with heated ground and gas emissions at the SRTKA and BCTKA. Hydrothermal systems are often observed or interpreted to exploit fault intersections which can serve as highly permeable pathways for hydrothermal fluid and gas discharge, enabling economic geothermal energy production. Geophysical mapping and modelling are an effective means of investigating such structural complexity at Mammoth Lakes due to the presence of unidentified and concealed structures.

California↗

USGS2026

Creation of aerial camera system boresight installations in Bangor, Maine; Monroe, Louisiana; and Bassville Park, Florida

Aerial camera systems rely on boresights for defining internal and external metrics. The U.S. Fish and Wildlife Service’s Migratory Bird Surveys Branch required boresights to calibrate their camera systems for creating georeferenced imagery that can then be used in geographic information systems. The U.S. Geological Survey created high-accuracy boresight installations in Bangor, Maine, in July 2024; Monroe, Louisiana, in September 2024; and Bassville Park, Florida, in January 2025 to provide this metric-defining capacity to the U.S. Fish and Wildlife Service. This U.S. Geological Survey Data Report details the methods and results of the boresight installation process and provides camera operators with the data needed to use each of the installations.

Florida, Louisiana, Maine↗

USGS2024

Current status of the community sensor model standard for the generation of planetary digital terrain models

The creation of accurate elevation models (topography) from stereo images are critical for a large variety of geospatial activities, including the production of digital orthomosaics, change detection, landing site analysis, geologic mapping, rover traverse planning, and spectral analysis. The United Stated Geological Survey, Astrogeology Science Center, continues to transition the supported planetary sensor models to the Community Sensor Model (CSM) standard. This paper describes the current state of use for this photogrammetric standard, supported sensor model types, and qualitatively compares derived topography between SOCET SET and SOCET GXP ( ® BAE Systems) using HiRISE stereo images of Mars. Our transition to the CSM standard will ensure an uninterrupted capability to make these valuable products for Mars and many other extraterrestrial planets and moons.

Remote Sensing↗

USGS2024

Debris-flow entrainment modelling under climate change: Considering antecedent moisture conditions along the flow path

Debris-flow volumes can increase along their flow path by entraining sediment stored in the channel bed and banks, thus also increasing hazard potential. Theoretical considerations, laboratory experiments and field investigations all indicate that the saturation conditions of the sediment along the flow path can greatly influence the amount of sediment entrained. However, this process is usually not considered for practical applications. This study aims to close this gap by combining runout and hydrological models into a predictive framework that is calibrated and tested using unique observations of sediment erosion and debris-flow properties available at a Swiss debris-flow observation station (Illgraben). To this end, hourly water input to the erodible channel is predicted using a simple, process-based hydrological model, and the resulting water saturation level in the upper sediment layer of the channel is modelled based on a Hortonian infiltration concept. Debris-flow entrainment is then predicted using the RAMMS debris-flow runout model. We find a strong correlation between the modelled saturation level of the sediment on the flow path and the channel-bed erodibility for single-surge debris-flow events with distinct fronts, indicating that the modelled water content is a good predictor for erosion simulated in RAMMS. Debris-flow properties with more complex flow behaviour (e.g., multiple surges or roll waves) are not as well predicted using this procedure, indicating that more physically complete models are necessary. Finally, we demonstrate how this modelling framework can be used for climate change impact assessment and show that earlier snowmelt may shift the peak of the debris-flow season to earlier in the year. Our novel modelling framework provides a plausible approach to reproduce saturation-dependent entrainment and thus better constrain event volumes for current and future hazard assessment.

Illgraben basin↗

USGS2024

Extending the Boore and Abrahamson (2023) modified square-root-impedance method for the development of site amplifications consistent with the full-resonance approach to a range of VS30 values

The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity ( V S ) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to 1.3 across a wide frequency range (Boore, 2013). Recently, Boore and Abrahamson (2023; hereafter, BA23) related this difference in the SRI and FR methods to differences in the exponent η of the ratio of seismic impedances between the two methods. They proposed the implementation of a modified frequency-dependent η in the SRI method to improve its match to the FR site amplifications. This modified η was derived using only five V S profiles. We investigate the performance of the BA23 η for a wide range of realistic gradient V S profiles with V S30 ranging from 180 to 1500 m/s. These gradient V S profiles are constructed using two power-law functions of depth and are constrained by the assigned VS30 value, the depth and velocity of the half-space, and depths to shear-wave velocity horizons of 1.0 and 2.5 km/s ( Z 1.0 and Z 2.5 ) based on western United States sites. Despite observing a V S30 dependence of η, we find that the BA23 η generally works reasonably well for the range of V S profiles analyzed. Using the VS30 -dependent η derived in this study results in improvements in matching the FR site amplification compared to using the BA23 η. These improvements are more pronounced for the soft-site conditions and become modest to negligible for the stiff site conditions

Bulletin of the Seismological Society of America↗

USGS2024

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan↗