Seismic anomalies expand Alaska’s Umiat Anticline potential
Recent seismic advancements, including AVO analysis, have redefined Alaska's Umiat field potential, revealing deeper reservoirs and increasing estimated oil in place.
Alaska↗
SEARCH · USGS Science
Results for “Frontiers in Marine Science”
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.
3-D mapping of the conterminous U.S. within the USGS National Cooperative Geologic Mapping Program: Progress and future prospects
The U.S. Geological Survey (USGS) National Cooperative Mapping Program (NCGMP) is bringing together subsurface and three-dimensional information at multiple scales for the conterminous United States from data produced throughout the USGS and by federal and state partners. Components of this work include data inventory and catalog development, data integration and database development, and construction of subregional and basin-scale three-dimensional (3-D) geologic models. At small scales, the NCGMP is compiling subsurface data to create a limited number of subsurface layers that map the majority of the lower 48 States. At larger scales, 3-D models of stratigraphic units are intended to be the subsurface analog and extension of the surface geologic maps made within the Program. The Mapping Program is developing methods of storing, visualizing, and distributing 3-D models and subsurface data from multiple sources in easily shareable, queryable, non-proprietary format(s). We report on progress in 3-D mapping of the conterminous U.S. and on future prospects.
conterminous United States↗
Hypothetical CO2 leakage into, and hydrological plume management within, an underground source of drinking water at a proposed CO2 storage facility, Kemper County, Mississippi, USA
A large Geologic Carbon Sequestration (GCS) hub has been proposed in Kemper County, Mississippi. The target injection interval consists of numerous Cretaceous-aged deep saline aquifers overlain by a competent and extensive regional sealing layer. Above the seal, the deepest Underground Source of Drinking Water (USDW) at the site is the Eutaw aquifer of the Eutaw Group and McShan Formation, undifferentiated. To assess potential risks of leakage from the deep sequestration reservoir, a model of a portion of the Cretaceous Eutaw Group was constructed in this study. Simulations tested various permeabilities, hypothetical leakage rates, and plume mitigation strategies utilizing existing wells. Results suggest that, under the influence of regional groundwater flow fields, leaking CO 2 would effectively bypass the existing wells, and to influence this migration would require very large water extraction rates. Therefore, to ensure plume detection, monitoring for leakage at the injection wells themselves is very important.
Mississippi↗
Temperate wetlands lose climate-cooling capacity under warming
Temperate inland wetlands have long been viewed as climate regulators, with sustained CO₂ sequestration offsetting continuous CH₄ emissions over centennial–millennial timescales. However, the net climatic effect of these opposing fluxes remains uncertain under anthropogenic warming as they differ in both magnitude and atmospheric lifetime. We combine 19,000 CH₄ chamber measurements from 202 temperate inland wetlands with carbon burial rates from 16 dated sediment cores in North America’s Prairie Pothole Region to estimate greenhouse gas dynamics over two centuries. Here we show that CH₄ emissions rise by ~1.39 kg CH₄ ha⁻¹ yr⁻¹ over the historical period and are projected to approach 950 kg CH₄ ha⁻¹ yr⁻¹ under high-forcing scenarios by 2100, while CO₂ sequestration is expected to decline by ~0.3 kg CO₂ ha⁻¹ yr⁻¹, driving a shift toward persistent net CO₂-equivalent source behaviour and a weakening of the long-term cooling effect that wetlands have provided under natural conditions.
Prairie Pothole Region↗
Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake
A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.
British Columbia, California, Oregon, Washington↗
Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States
The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.
Arizona, California, Colorado, Idaho, Montana, Nev↗
Automated, near real-time ground-motion processing at the U.S. Geological Survey
We describe automated ground‐motion processing software named gmprocess that has been developed at the U.S. Geological Survey (USGS) in support of near‐real‐time earthquake hazard products. Because of the open‐source development process, this software has benefitted from the involvement and contributions of a broad community and has been used for a wider range of applications than was initially envisioned. Here, we give an overview and introduction to the software, including how it has leveraged other open‐source libraries. We highlight some key features that gmprocess provides, compare response spectra calculated with the automated processing approach of gmprocess to the response spectra provided by the Next Generation Attenuation projects, and summarize projects that have utilized gmprocess. These use‐cases demonstrate that this software development effort has been successfully leveraged in earthquake research activities both within and outside the USGS.
Seismological Research Letters↗
Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023
The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.
Alaska↗
Zircon as a pathfinder to REE mineralization
Carbonatites and alkaline silicate rocks are major primary sources of the rare earth elements (REE) and other critical metals, such as Nb. Despite the economic significance of these rocks, their formation and the processes of REE enrichment are poorly understood. Here, statistical analysis of a global dataset demonstrates that zircon geochemistry is a powerful recorder of REE metallogenesis and a potential pathfinder for REE deposits. Zircons from REE and Nb fertile intrusions lack Eu anomalies and have elevated Gd/Yb and Th/Yb, indicating they crystallised from magmas that originated from deep, oxidised and enriched mantle sources. Complexes with Nb enrichment have low U/Nb, reflecting an enriched mantle source, whereas high U/Nb in REE-only fertile intrusions suggest a subduction-metasomatised mantle source. Machine learning models demonstrate high accuracy in classifying zircon from barren and fertile deposits. Classification of detrital zircons shows that REE-enriched deposits correlate with supercontinent assembly, whereas Nb fertile complexes are associated with supercontinent breakup. This approach offers a new, mineral to global scale, petrologic and exploration tool that enhances understanding of REE metallogenesis.
Geochemical Perspectives Letters↗
A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes
Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.
Seismological Research Letters↗
Dating a medieval tsunami with uranium-series techniques on Caribbean corals
Uranium-series dates from coral boulders constrain the timing of a medieval tsunami from the Puerto Rico Trench. Previously reported evidence for this tsunami includes hundreds of coral boulders that came to rest hundreds of meters inland on Anegada, British Virgin Islands. New U-series dates on these coral boulders provide limiting dates for the tsunami. The narrowest limits were by dating interior bands of a coral that retains the hemispherical form of a living coral colony, and which include adjustments for the number of annual density band couplets between the dated samples and the boulder exteriors. By those limits, the tsunami dates between 1381 and 1391 CE, and likely occurred during summer or fall. The tsunami is important as the only reported sign that the eastern Puerto Rico Trench has produced a great earthquake. The dating may aid in defining the earthquake source and in communicating tsunami hazards.
British Virgin Islands, Puerto Rico↗
Ground magnetic and VLF studies at Midnite uranium mine, Stevens County, Washington
Ground magnetic and E-mode VLF measurements were made at more than 200 stations at the Midnite uranium Mine, Stevens County, Washington. The equipment was found to be light and foolproof, and each set of data was collected in one day, exclusive of surveying and marking the station points. The maps of contoured magnetic and VLF data resembled geologic maps of the area, showing many common features, as well as a few unique features on each map. Each showed some evidence of a small normal fault cutting the region; taken together, the evidence for this fault was judged to be much stronger, so that further structural details were inferred. Because of soil cover and inadequate outcrops, the presence of this fault had hitherto been unsuspected. A target area for further uranium mineralization north of the fault was indicated by this quick, inexpensive geophysical survey. The drilling program which Midnite Mine management had underway at that time had already outlined a mineralized zone in the indicated area.
Washington↗
Field evidence and indicators of rockfall fragmentation and implications for mobility
Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.
Albacete province, Lombardy and Aosta Valley, Yose↗
The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes
The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.
Bulletin of the Seismological Society of America↗
Groundwater budget for the Mountain Home area, southern Idaho, 2022–23
The U.S. Geological Survey, with funding from the Idaho Department of Water Resources, developed a groundwater budget for the Mountain Home area in southern Idaho for irrigation year 2023 (November 1, 2022–October 31, 2023). This study focused on the water balance across the Cinder Cone Butte Critical Groundwater Area (CGWA), Mountain Home Groundwater Management Area (GWMA), and the rest of the study area (RoSA), compiling data from various sources, including precipitation records, groundwater level measurements, metered groundwater pumpage data, surface water diversions and evapotranspiration (ET) estimates derived from remote sensing satellite imagery, and ground-based reference data. Key inflow components included recharge from applied surface water irrigation (which incorporates incidental recharge from irrigation practices and conveyance losses), estimated tributary streamflow, and estimated mountain block recharge. The key outflow components were groundwater pumpage for irrigation, municipal, industrial, and domestic uses, and ET. Recharge from applied irrigation and mountain block recharge were the largest inflows, and groundwater pumpage for irrigation was the largest outflow. The CGWA had a positive groundwater budget residual of 2,170 acre-feet (acre-ft), which contrasts with observed long-term groundwater level declines and historical trends of storage depletion. This positive residual is likely associated with unquantified outflows, including lateral groundwater flow out of the subregion, or other complexities, such as overestimated tributary contributions relative to the actual recharge for the 2023 water budget. The GWMA exhibited a positive residual of 56,563 acre-ft, primarily owing to recharge from applied surface water irrigation and areal recharge during a wetter-than-average year, which allowed irrigation entities to deliver more water from in-basin and out-of-basin reservoirs. The RoSA showed a large positive residual of 124,933 acre-ft. The interpretation of these positive residuals must account for significant uncertainties, including estimations of areal recharge, tributary streamflow (particularly losses and diversions), ET, the volume of surface water loss to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and the unquantified groundwater discharge to the Snake River. These uncertainties, in combination with the complex hydrogeologic controls on water movement and limitations of remotely sensed data, directly affect the accuracy of water availability assessments. Future data collection efforts would help reduce these uncertainties and support water resource management decisions in the Mountain Home area. Key efforts could include installing additional streamflow gaging stations (particularly to quantify tributary losses and gains and surface water losses to the Snake River), improving groundwater pumpage metering, and validating remotely sensed ET data with ground-based measurements. Furthermore, to better quantify unrepresented or highly uncertain fluxes, focused investigations on groundwater discharge to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and a more robust determination of the actual influence and volume of mountain block recharge would help refine future water availability assessments for the Mountain Home area.
Idaho↗
Mapping a Carrington storm
A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.
contiguous United States↗
The saline groundwater legacy of a large buried coastal paleo-estuary
Elevated groundwater salinity in coastal regions threatens the beneficial use of fresh groundwater. Coastal groundwater management typically focuses on preventing intrusion from modern sources of seawater; however, past geological processes can also leave a legacy of saline groundwater now hidden in the subsurface. Here, multiple extensive airborne electromagnetic surveys provide detailed evidence of residual salinity from a paleo-estuary filling a late Pleistocene incised valley impacting more than 10,000 km 2 that is now hidden beneath coastal Louisiana’s deltaic plain. Our results show that the three-dimensional pattern of saline groundwater beneath Louisiana mimics that of near-surface aquifers surrounding the modern Delaware Bay estuary, fingerprinting the signature of the past drowning of a large, incised valley of the Mississippi River following post-glacial sea-level rise. These findings demonstrate a new framework for understanding legacy sources of saltwater critical for managing stressed water resources along global coastlines.
Louisiana↗
Waters divided: A history of alluvial fan research and a view of its future
Flows exiting confined valleys tend to deposit sediment in fan-shaped landforms. Where deposition is wholly or largely by the tractive forces of flowing water, these landforms are called alluvial fans. They are the product of the progressive division of water and sediment downfan, from slopes that may exceed 0.10 to distal slopes that may be below 0.01. Channel depths also tend to decline, from values that approach one to several meters at steep fanheads, to a few decimeters at distal fan margins. The result is a radiating, depositional ramp where confined or unconfined flows transport sediment from source basins to bounding streams, subsiding basins, or stable platforms. Where streams or subsiding basins consume the sediment supply from the source basin, the fan may approach a steady form whose extent and distal slope are set by stream location or subsidence rate. Where boundary conditions do not remove sediment, the fan may prograde out to long distances and low slopes (<0.01). Theoretical and experimental work over the past several decades support the notion that alluvial fan long-profiles become steeper as sediment supply increases or transport capacity decreases, and increasingly concave upward as the rate of bed material deposition decreases downfan. Grainsize distributions of alluvial fans seem to span the range observed in alluvial rivers, with no processes that uniquely identify them, apart from the distributary pattern of deposition. Bed sand cover tends to increase downfan in arid-region fans, with an absence of systematic downfan fining of coarser grain sizes. Surficial mapping and geochronology have demonstrated that fan deposition varies greatly through time, arguably from climate variations that alter hillslope sediment supply. The combination of surficial mapping and hydraulic modeling with high-resolution topography can now produce detailed flood susceptibility maps. The effective use of these maps to protect lives and property, however, depends on answering many of the enduring questions about the mechanics of how water and sediment divide down alluvial fans.
Book chapter↗