SEARCH · Geology Reports

Results for “Frontiers for Young Minds”

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.

1,406 records · Page 8Linked to original sources

USGS2025

The feasibility of using lidar-derived digital elevation models for gravity data reduction

Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Open-File Report↗

USGS2026

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↗

USGS2025

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan↗

USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗

USGS1995

Profiles of gamma-ray and magnetic data for aerial surveys over parts of the Western United States from longitude 108 to 126 degrees W. and from latitude 34 to 49 degrees N.

This CD-ROM contains images generated from geophysical data, software for displaying and analyzing the images and software for displaying and examining profile data from aerial surveys flown as part of the National Uranium Resource Evaluation (NURE) Program of the U.S. Department of Energy. The images included are of gamma-ray data (uranium, thorium, and potassium channels), Bouguer gravity data, isostatic residual gravity data, aeromagnetic anomalies, topography, and topography with bathymetry. This publication contains image data for the conterminous United States and profile data for the conterminous United States within the area longitude 108 to 126 degrees W. and latitude 34 to 49 degrees N. The profile data include apparent surface concentrations of potassium, uranium, and thorium, the residual magnetic field, and the height above the ground. The images on this CD-ROM include graytone and color images of each data set, color shaded-relief images of the potential-field and topographic data, and color composite images of the gamma-ray data. The image display and analysis software can register images with geographic and geologic overlays. The profile display software permits the user to view the profiles as well as obtain data listings and export ASCII versions of data for selected flight lines.

Data Series↗

USGS1976

Digital telemetry system network installation for the National Earthquake Information Service

A Digital Telemetry System Network, developed by the Albuquerque Seismological Laboratory, was installed for the USGS National Earthquake Information Service (NEIS) during January 1975. NEIS is located in Golden, Colorado. The initial Digital Telemetry System Network consists of four independent digital telemetry systems that provide wide area coverage over most of the United States. The four systems were installed between: Albuquerque, New Mexico to NEIS; St. Louis, Missouri to NEIS; Salt Lake City, Utah to NEIS; and Menlo Park, California to NEIS. Each digital telemetry system will transmit fourteen channels of high accuracy seismic data in a digital format over a standard C-1 grade telephone line. Digitized data words will have a resolution of eleven digital bits plus one sign bit. Each digital telemetry receiver, located at NEIS, will have the capacity to provide fourteen channels of seismic data in a digital format suitable for either on-line digital computer analysis or magnetic tape recording. In addition the data can be presented in an analog format suitable for recording on helicorders or develocorders. This network is the first large scale digital telemetry system used by the U.S. Geological Survey for the transmission of real time seismic data.

Open-File Report↗

USGS2026

Methodology for quantitative spatial sensitivity analysis of volcanic geodetic networks

Introduction This report introduces a methodology for assessing the state of the U.S. Geological Survey Volcano Observatories’ geodetic monitoring networks that measure how volcanoes deform or change shape. This new method uses a model-based approach that considers the uniqueness of the instrument environments at each volcano. This report focuses on simplified volcanic sources, is independent of the shape or size of the volcano, or the network geometry, and thus highlights the strengths and potential vulnerabilities of each volcano’s geodetic network in an actionable visual format. This analysis can help observatories to make informed decisions about whether volcanoes have an adequate level of geodetic monitoring and indicate where improvements are needed.

Lassen Peak, Mount Shasta↗

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

USGS2024

Evaluating the impact of uncertainty in ground motion forecasts for post-earthquake impact modeling applications

The US Geological Survey’s (USGS) ShakeMap system provides a rapid characterization of strong ground shaking in areas directly affected by an earthquake. This study focuses on studying the aggregate effects of macroseismic shaking estimates from ShakeMap, expressed in terms of modified Mercalli intensity (MMI), when accounting for the uncertainty in forecasted ground motions. We use a Monte Carlo approach to generate numerous spatially correlated realizations of ground motions by utilizing a combination of circulant embedding and kriging techniques for efficiently handling the correlations. We then assessed the aggregate effects of shaking by looking at bin counts across these realizations. We demonstrate that the aggregate shaking regarding the mean macroseismic intensity estimates (from the ShakeMap output) is a biased representation of the aggregate shaking when shaking uncertainty is included. Incorporating shaking uncertainty can help to improve various downstream earthquake impact applications, such as the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) overall earthquake fatality distribution or estimates of shaking-induced ground failure impacts from consequential earthquakes.

Earthquake Spectra Journal↗

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

Flood-inundation maps for the Cuyahoga River at Jaite, Ohio, 2024

Digital flood-inundation maps for a nearly 6-mile reach of the Cuyahoga River at Jaite, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. The maps depict estimates of the extent and depth of flooding corresponding to selected water levels (stages) at USGS streamgage 04206425 on the Cuyahoga River at Jaite, Ohio. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The hydraulic model was calibrated to the current USGS streamgage data and then used to compute 15 water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from 6 to 20 feet, which correspond to below “action stage” to “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used with a geographic information system digital elevation model derived from light detection and ranging data to delineate the areas flooded at each stage. These maps, along with current stage data from the USGS streamgage and forecasted high-flow stages from the National Weather Service, can provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Ohio↗

USGS1990

Preconditioned Conjugate-Gradient 2 (PCG2), a computer program for solving ground-water flow equations

This report documents PCG2: a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model. PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options: modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers. Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Water-Resources Investigations Report↗

USGS2025

Understanding market sensitivity: Estimation of supply and demand elasticities for non-fuel minerals

In today’s rapidly changing economic landscape, understanding market responsiveness to price changes and the factors influencing commodity prices has become increasingly relevant. Price elasticities serve as indicators of how variations in market conditions affect supply and demand, providing insights into the sensitivity of commodity markets to price fluctuations. This paper presents a comprehensive analysis of price elasticities of supply and demand for 74 non-fuel mineral commodities including precious metals, base metals, minor metals, and industrial minerals that are utilized across various industries. We employ various econometric techniques, including fixed effects models for panel data and two-stage dynamic ordinary least squares (2S-DOLS) alongside autoregressive distributed lag (ARDL) models for time series analysis, to derive robust estimates of price elasticities. Our findings reveal variability in elasticities among different commodities and indicate that all studied mineral commodities exhibit price inelastic supply and demand in the short run, which we define as one year for the purposes of our analysis, given that the data is all annual. This research provides original estimates of price responsiveness for a wide range of commodities that have not been previously addressed in the literature, thereby enhancing the understanding of market dynamics in the mineral sector. Given that price elasticities can be influenced by factors such as market structure, technological advancements, mining costs, and industry-specific demand drivers, we use variables that serve as proxies for these factors.

Mineral Economics↗

USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗

USGS2026

Where will the cat cross the road? Comparing camera and GPS-based models for identifying wildlife corridors

Designing effective wildlife corridors is a critical conservation challenge in fragmented landscapes. GPS-based step selection functions strongly predict dispersal corridors and connectivity, but GPS collaring can be expensive and invasive. Camera-based occupancy models are widely used for connectivity analyses but may involve trade-offs in data resolution. Despite widespread use of both approaches, few studies have directly compared them using concurrent datasets. We developed a stacked single-species, single-season occupancy model and a Circuitscape connectivity surface for mountain lions (Puma concolor) on Washington’s Olympic Peninsula, USA, and compared them with a connectivity surface from an existing integrated step selection function. Both models predicted mountain lion GPS locations well, with binned Spearman rank correlations of 1 for Circuitscape and 0.96 for the step selection function, though step selection better identified habitat use by dispersers. Connectivity predictions were moderately correlated across the landscape ( r = 0.26), but agreement was strongest in human-dominated areas most critical for corridor planning. We conclude that GPS-based approaches are advantageous when data collection is feasible and the focus is on dispersal or fine-scale movement. However, camera-based approaches may be preferable for multi-species monitoring, large spatial and temporal scales, noninvasive sampling, when resources are limited, or when fine-scale or dispersal-specific inference is not required.

Washington↗

USGS2025

Prospectivity mapping for geologic hydrogen

Geologic, or naturally occurring, hydrogen has the potential to become a new, low-carbon, primary energy resource. Often referred to as “white” or “gold” hydrogen, this gas occurs naturally in the Earth’s subsurface, similar to petroleum resources. However, unlike petroleum, which releases carbon dioxide when burned, burning hydrogen only produces water as a byproduct. Exploration for geologic hydrogen remains in an early stage and discoveries of high concentrations of subsurface hydrogen are still relatively rare. To facilitate research and exploration for this potential resource, this report presents the first publicly available prospectivity map of geologic hydrogen accumulations in the conterminous United States. Prospective regions are those regions in which all major components necessary for a hydrogen accumulation likely are present—a source of sufficient hydrogen generation, porous reservoirs for storage, and seals to prevent leakage. The midcontinent region of the United States and the central California coast are revealed as having high prospectivity. This analysis also identifies previously unrecognized prospective regions that may be favorable due to long distance lateral migration of subsurface hydrogen, such as the offshore eastern seaboard of the United States, and can provide a linkage between surface observations of hydrogen degassing and far-field source regions. The methodology developed to create this map is expandable and flexible and may be adapted to incorporate new concepts in the hydrogen system and for application to other regions of the world.

conterminous United States↗

USGS2025

Meet the people where they are: Assessing user needs for aftershock forecast products in El Salvador, Mexico and the United States

Aftershock forecasts can help communities reduce their seismic risk by conveying how many aftershocks can be expected following a large earthquake, and how the expected number of aftershocks and their corresponding ground shaking evolves over time and space. Prior work finds that graphical forecast products may communicate such information better than only text or numbers. To identify which visual products can serve multiple user groups, we held workshops with members of several professions, including emergency managers, engineers, critical infrastructure operators, public health specialists, science communicators, and more. We conducted these workshops in El Salvador, Mexico and the United States to understand which forecast products may be effective across different countries. In these workshops, users performed small-group activities to elicit the types of aftershock forecast information that would support decisions in their respective roles and how this information would optimally be displayed. Maps of shaking hazards were frequently requested across all professions and countries, even for dissimilar forecast uses. The design of these maps, and other product needs, were differentiated by profession, country and other factors. Other forecast products, including those showing forecasts about the magnitudes and time periods of future aftershocks, also served a variety of users, but for different types of decisions. We found a greater variation in user needs by profession than by country, and that user needs also vary with time, communication channel and other contextual factors. We discuss practical implications for user-centered visual communication of operational aftershock forecasts.

International Journal for Disaster Risk Reduction↗

USGS2026

Re-examining the growing degree day minimum for grass carp spawning initiation in Lake Erie: Have we missed the mark?

Understanding the timing of life history processes of invasive fish species, such as adult aggregations during spawning, allows natural resource managers to optimize targeted population control efforts. Temperature-based estimates of the timing of spawning initiation for grass carp ( Ctenopharyngodon idella ) exist from its native range, but have not previously been developed for grass carp reproduction in its invaded range in the Laurentian Great Lakes of North America. We accounted for differences among methods across studies and calculated cumulative growing degree days (GDD) for the dates of grass carp egg capture in Lake Erie tributaries (Sandusky, Maumee, and Huron rivers) during 2015–2024 and compared these values to the existing GDD base 15°C requirement estimate from the native range (633 GDD–15°C). The earliest egg detection in Lake Erie (Sandusky River) was 19 GDD–15°C or 361 GDD–5°C (latest egg detection was 2285 GDD–5°C from the Sandusky River), indicating that spawning was initiated up to 25 days earlier than would be predicted by the native range GDD. Moreover, spawning in Lake Erie tributaries occurred earlier than predicted for 58% of spawning events examined based on the native range GDD. The GDD values for spawning initiation in Lake Erie tributaries can be used to refine the timing of control efforts targeting pre-spawn and spawning adult fish, inform the timing of egg sampling, and provide linkages to existing early life stage models to infer the length of the growing season available for age-0 fish to reassess tributary suitability.

Lake Erie, Maumee, River, Sandusky river↗