USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Modelling & Software”

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.

834 records · Page 22Linked to original sources
USGS2025

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗
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

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2025

Central Valley Hydrologic Model version 2 (CVHM2): Decision support tool for groundwater and land subsidence management

The San Joaquin Valley (SJV) of California is one of the world’s most productive agricultural regions. Reliance on groundwater has led to some of the greatest rates of human-induced land subsidence in the world in the 20th century, as well as more recently. The United States Geological Survey (USGS) has recently developed an integrated surface–subsurface hydrologic model, the Central Valley Hydrologic Model 2 (CVHM2), that represents the major components of the hydrologic system of California’s Central Valley. In this study, CVHM2 was applied as a decision support tool while simulating various management strategies to mitigate the land subsidence caused by the extraction of groundwater. CVHM2 was extended through to 2073 and applied to simulate management scenarios in terms of three primary drivers and their impact on subsidence along the Delta–Mendota Canal (DMC), a critical piece of infrastructure in the western SJV. The drivers considered were agricultural water demands, managed aquifer recharge (MAR), and changes in future climate. The results show that future subsidence is most sensitive to water demands, second most sensitive to future changes in climate, and relatively insensitive to MAR when it is applied as a surface application in the western SJV. However, we demonstrate via proof-of-concept scenarios that the MAR is capable of arresting subsidence when implemented via injection below the Corcoran Clay Member of the Tulare Formation instead of as a surface application. We also examine the uncertainty that is the result of climate variability and how to use the tool to identify the most appropriate strategies to constrain future subsidence to acceptable levels.

California↗
USGS2025

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗
USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗
USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗
USGS2025

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↗
USGS2025

Predictive regressive models of recent marsh sediment thickness improve the quantification of coastal marsh sediment budgets

Coastal marsh wetlands experience variations in vertical gains and losses through time, which have allowed them to infill relict topography and record variations in drivers. The stratigraphic unit associated with the development of the marsh also reflects the long-term importance of key ecosystem services supplied by the marsh environment, including carbon storage and storm mitigation. Mapping these coastal wetland sediments and the marsh unit thickness is challenging as traditional coastal geophysical tools are not easily deployable (acoustic methods) or are unreliable in saline-soil environments (e.g., ground-penetrating radar), leaving core-based methods the most viable mapping method. In the present study, we utilized prior information on the geologic architecture of the region to select spatial and physical metrics that likely persisted throughout evolution of the marsh during the late Holocene. We then assessed the individual and collective power of these metrics to predict marsh thickness observed from cores. Employing regressive predictive models powered by these data, we improve the quantification of marsh thickness for a coastal fringing marsh within the Grand Bay estuary in Mississippi and Alabama (USA). The information gained from this approach yields improved estimates of the carbon stocks in this environment. Additionally, the stored sediment masses reflect the past, and potential future, persistence of the Grand Bay marsh under historical and present marsh-estuarine sediment exchange fluxes. Such improvements to both the sediment budget of recent marsh stratigraphic units and the spatial extent provide new resources for comparison with large-scale landscape models, the latter of which may be used, when validated, to predict future change and ecosystem transformations.

Alabama. Mississippi↗
USGS2026

Groundwater and surface-water interactions in the Waihe‘e-Kahalu‘u watershed, O‘ahu, Hawai‘i—Analysis of historical data and numerical groundwater-model simulations

The Waihe‘e-Kahaluʻu watershed and surrounding areas lie in a dike-impounded groundwater setting of the Ko‘olau Range, O‘ahu, Hawai‘i, where groundwater withdrawal from wells and tunnels have resulted in reductions in streamflow. Resource managers and surface-water users seek information that can be used to balance the needs of competing uses of groundwater and streamflow in the watershed. This study assesses the effects of four groundwater developments in the Waihe‘e-Kahalu‘u watershed—Kahaluu Tunnel, Waihee Tunnel, Kahaluu Well, and Waihee Incline Wells—on Kahalu‘u and Waihe‘e Streams. Analyses of historical streamflow, withdrawal, and rainfall data indicate that base flow (the component of streamflow that is supplied from groundwater discharge) in Kahalu‘u and Waihe‘e Streams decreased in response to withdrawals from wells and tunnels within and outside the watershed. Groundwater-model simulations indicate that shutdown of selected tunnels and wells in the Waihe‘e-Kahaluʻu watershed will result in increases of groundwater discharge, but not all of the increase will occur in streams within the watershed; increases will also occur in streams outside the basin, the ocean, and tunnels other than those that were shut down. The analysis of historical data and model simulations is consistent in showing that the effects of withdrawal changes on streams are not constrained by the topographic boundaries of watersheds. Stream base flows in the Waihe‘e-Kahalu‘u watershed are subject to a complex interaction of multiple withdrawals from wells and tunnels installed at various times within and outside the watershed. The correspondence between withdrawals and base-flow changes is consistent with the conceptual model and understanding of the connection between groundwater and streamflow in the dike-impounded groundwater setting of the Waihe‘e-Kahalu‘u watershed and the Ko‘olau Range.

Hawaii↗
USGS1999

Estimation of magnitude and frequency of floods for streams in Puerto Rico: New empirical models

Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.

Puerto Rico↗
USGS2024

Geophysical modeling of a possible blind geothermal system near Battle Mountain, NV

The northeastern portion of the Reese River basin in north-central Nevada is the focus of detailed geophysical and geological studies as part of the INGENIOUS project, which aims to identify new, commercially viable hidden geothermal systems in the Great Basin region of the western U.S. This location, herein referred to as Argenta Rise, occupies a broad (~15km wide) left-step between major range-front fault systems along the northwestern edge of the Shoshone Range and Argenta Rim, with numerous ENE-striking intra-basin faults presumably accommodating sinistral-normal oblique slip across the step-over. Four discrete regions have been identified within the study area that have favorable structural settings for hosting a blind hydrothermal system. However, with no definitive or extensive surface manifestations of an active hydrothermal system (e.g., geysers, steam vents, sinter, etc.), detailed geophysical studies are necessary to resolve subsurface geology and structure, and identify zones of enhanced structural complexity that may promote hydrothermal fluid flow. Hence, we collected high-resolution gravity, MT, and rock property data (density, magnetic susceptibility), and analyzed the recently acquired GeoDAWN aeromagnetic data to characterize potential geothermal resources in this region. Using the new geophysical datasets, we jointly modeled gravity and magnetic data along a series of intersecting 2D profiles that integrated information from recent, local-scale fault mapping. Rock property measurements performed on outcrops and hand samples throughout the study area constrained the models. The MT data were used to construct a 3D resistivity model that highlights the location of inferred alteration and fluids in the subsurface. Combined MT and potential field results reveal which structures may be most important for controlling hydrothermal fluid migration, as well as which geologic units may host hydrothermal fluids. Our gravity derived depth to basement surface coincides well with the base of shallow conductive anomalies, suggesting hydrothermal fluids may be confined to basin fill sediments and volcanics. This work supports our development of 3D geophysical and geologic models that are focused along the western flank of the northern Shoshone Range and aids the process of selecting sites for temperature gradient drilling.

Nevada↗
USGS2025

Utilizing downhole datasets for modelling the aeromagnetic signature of the Iron Creek Co-Cu deposit in the Idaho Cobalt Belt

The Idaho Cobalt Belt in east-central Idaho is host to some of the largest domestic Co resources, including the informal Iron Creek deposit. The two main ore zones of this deposit, the Iron Creek and the Ruby, are hosted in greenschist-grade interbedded argillite/siltstone and quartz-rich units of the Mesoproterozoic Apple Creek Formation of Lemhi Group. The primary ore mineral is cobalt-bearing pyrite, which occurs with pyrrhotite, chalcopyrite, and magnetite. This study integrates mineralogical and core-scale geophysical measurements with regional-scale aeromagnetic surveys. The high magnetite content within the Ruby zone produces elevated magnetic susceptibility, but the relatively limited spatial footprint of the ore zone would produce a small-scale anomaly that may be overlooked in regional surveys. The low magnetite content in the Iron Creek zone results in low magnetic susceptibility, creating a relatively low amplitude geophysical response. By characterizing the magnetic properties and mineralogy of these ore zones, this study enhances the interpretation of aeromagnetic data, enabling the identification of small or faint anomalies as potential Co targets. These findings improve can improve exploration strategies, both within the Idaho Cobalt Belt as well as for similar deposit types globally.

Idaho↗
USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗
USGS2025

Comparing magmatism and hydrothermal alteration using magnetic modelling and stable isotopes at the Clementine porphyry copper prospect, Montana, USA

Recent mapping of the underexplored Clementine prospect in southwestern Montana has revealed evidence of hydrothermal alteration and mineralized breccia vein gossan interpreted to represent the upper expression of a Cretaceous, sediment-hosted copper porphyry system. The prospect is at the nexus of several Cretaceous granites, including the pre-mineralization Butte Granite and Big Hole Canyon plutons and local pre- to post-mineralization granite intrusions. Here, 3D magnetic inversions and stable isotope data are used to evaluate the spatial and genetic relationship of mineralization and Cretaceous magmatism. Magnetic inversions reveal a zone of high magnetic susceptibility beneath the prospect that is related to an exposure of unaltered, likely post-mineralization granite. The granite appears to connect to the Big Hole Canyon pluton at depth, suggesting it is a late-stage differentiate of the cooling, pre-mineralization pluton that may have been one of several potential heat sources for the mineral system at Clementine. This is supported by carbon isotope values of graphite that indicate organic material in shale of the Cretaceous Kootenai Formation was subject to metamorphism. Preliminary light hydrogen isotope values of actinolite hydrothermal orbs suggest formation following magmatic outgassing during thermal advection from an underlying intrusion.

Montana↗
USGS2026

Hydrogeologic framework and conceptual groundwater-flow model of the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, 1998–2022

This study was conducted by the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, to update the hydrogeologic framework and conceptual flow model for the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, which together compose the Ogallala aquifer focus area. The study included the construction of a potentiometric surface, and available geologic and hydrologic data were used to evaluate saturated thickness of the aquifer. The water budget for the updated conceptual groundwater-​flow model was based on estimated inflows and outflows for the 1998–2022 study period. Saturated thickness of the Ogallala aquifer averaged 127 and 116 feet for the panhandle and northwest parts, respectively. Groundwater withdrawals from the Ogallala aquifer for 1998–2022 averaged 422,054 and 39,645 acre-feet per year (acre-ft/yr) for the panhandle and northwest parts, respectively. Recharge, the primary inflow, was estimated at 0.63 inch per year for the 1998–2022 study period, with the panhandle part of the Ogallala aquifer receiving 175,068 acre-ft/yr and the northwest part of the Ogallala aquifer receiving 49,376 acre-ft/yr. Additional inflows included irrigation return flows, estimated at 8,111 and 642 acre-ft/yr for the panhandle and northwest parts, respectively, of the Ogallala aquifer. Net lateral groundwater flows, considered to be aquifer outflows, were estimated to account for 31,908 acre-ft/yr for the Ogallala aquifer focus area. Streambed seepage, which was an outflow of 5,535 acre-ft/yr, was only present in the northwest part of the Ogallala aquifer. Vertical leakage and saturated-zone evapotranspiration were considered negligible outflows. These findings provide a revised conceptual groundwater-flow model water budget for the Ogallala aquifer focus area in Oklahoma.

Colorado, Kansas, Oklahoma, Texas↗
USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗
USGS2026

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li 2 O, respectively, for global Li pegmatite deposits ( n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs 2 O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb 2 O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Ore Geology Reviews↗