Technique for concentrating nannoplankton from the Tertiary rocks of the California Coast and Peninsular Ranges
No abstract available.
Journal of Research of the U.S. Geological Survey↗
SEARCH · Search USGS
Results for “Journal of Geotechnical and Geoenvironmental Engineering”
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.
USGS 2025 critical minerals review
The United States Geological Survey (USGS) provides scientific information for the Department of Interior and the nation, consistent with its original mission expressed in the Organic Act of 1879 (43 U.S.C. 31): “the classification of the public lands and examination of the geological structure, mineral resources, and products within and outside the national domain.” Legislation such as the Energy Act of 2020 and the 2022 Infrastructure Investment and Jobs Act (43 USC 31l) and recent executive actions (Executive Orders 14154 , 14153, 14241 Secretary’s Orders 3417, 3418, 3422, 3436) underscore the importance of mineral resources and focus USGS activities on mapping and assessing mineral resources, with a particular focus on those presently identified as critical, both in ground and above ground in mine wastes. This article reviews selected activities and accomplishments by the USGS Mineral Resources Program related to critical minerals in 2025. Highlights include a new List of Critical Minerals, a first-ever national mine waste inventory, international minerals partnerships, mineral resource assessment advancements, and national data collection activities and outcomes of the Earth Mapping Resources Initiative (Earth MRI). The selected contributions are not comprehensive but are intended to demonstrate USGS leadership in critical mineral mapping and assessment, the importance of domestic and global partnerships, and the breadth of research activities that are responsive to national needs and priorities.
Mining Engineering↗
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↗
HyFlood: A surrogate-model-based framework for compound coastal flooding
Compound coastal flooding is a major threat to low-lying coastal regions and is expected to intensify under future climate change projections. However, modeling the joint interaction of waves, storm surge, tides, and rainfall remains computationally demanding, limiting the development of fast and reliable forecast tools. Here we present HyFlood, a hybrid statistical-numerical downscaling framework capable of computing and mapping high-resolution compound flood hazards while substantially reducing the computational cost compared with fully process-based hydrodynamic modeling. HyFlood combines statistical sampling and selection algorithms with a cascade of reduced-complexity surrogate models that emulate nearshore wave transformation, surf-zone hydrodynamics, and coastal, fluvial, and pluvial flooding. The surrogate models employ machine-learning and regression algorithms applied to a low-dimensional representation of the flooding outputs, obtained through statistical dimensionality reduction. The framework is demonstrated in southern O'ahu, Hawai'i, a region exposed to elevated sea levels driven by tides, waves, and storm surge along with frequent precipitation-driven flash flooding. Validation of the surrogates against the physics-based model outputs demonstrates that HyFlood accurately reproduces daily maxima of spatially distributed flooding depths. This hybrid approach offers a scalable and efficient tool to better quantify how changes in flooding drivers translate into hazard and impact assessments, and to support compound-flood risk assessments and climate-change adaptation planning.
Hawaii↗
From landslide susceptibility to risk assessment in the conterminous U.S.
Understanding the spatial distribution of landslide prone-areas and what consequences they may have is important for risk management and land-use planning. In the United States, although landslides occur in every state, a comprehensive landslide risk assessment is still missing. Existing efforts, such as the Federal Emergency Management Agency (FEMA)’s National Risk Index, rely on aggregated products and coarse cartographic units, limiting their geomorphological and practical accuracy. In this study, we present a methodological advance for landslide risk assessment across large areas with incomplete and sparse data. We apply our procedures to the conterminous United States by integrating geomorphologically meaningful partitions and spatial and temporal probability data-driven models. Landslide susceptibility is estimated using a Generalized Additive Mixed Model incorporating a bias capture/correction scheme to account for inventory inaccuracies (reference Area Under the Curve = 0.75). The exceedance probabilities of landslide occurrence are defined for three temporal scenarios (2, 5, and 10 year). Then, we explore the associated potential economic consequences for human settlements and agricultural areas. The findings indicate that the spatial variability of risk is primarily controlled by exposure rather than by susceptibility/hazard alone. The mean risk increases by ∼170% from the 2-year to the 10-year scenario. Beyond its quantitative outcomes, this study offers a blueprint for continental or sub-continental scale landslide risk assessments, demonstrating both the opportunities and current limitations.
Engineering Geology↗
Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends
Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.
Engineering Geology↗
Development of liquefaction-targeted design basis in U.S. seismic provisions
Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.
conterminous United States↗
Spatio-temporal modeling for assessing geoenergy resources: A workflow applied to gas in place variation in coal beds
The ability to estimate spatio-temporal changes in hydrocarbon reservoir properties and energy resources within pore volumes is essential for optimizing production, reservoir management, geologic energy storage, and safety in underground mining operations. In coal seams, predicting remaining methane gas-in-place (GIP) is critical for quantifying producible gas and improving mine safety and productivity through effective ventilation planning. Although such changes are commonly evaluated using physics-based numerical simulation models, these approaches often require extensive data, calibration effort, and time. This study presents a spatio-temporal geostatistical modeling approach that bridges the gap between purely spatial models and full numerical simulations. The method is applied to a case study of coal seam degasification in the Mary Lee coal group, Black Warrior Basin, Alabama, USA, to estimate GIP evolution over time within a selected mining district. The analysis uses published data from prior natural gas production history-matching of degasification using vertical wells. Empirical spatial and temporal statistics were calculated for reservoir pressure and water saturation, and spatio-temporal variogram models were fitted to experimental variograms. These models provided the structural basis for spatio-temporal kriging, integrated with spatial estimates of time-invariant parameters (porosity, density, and thickness) to estimate GIP. This approach enabled estimation of GIP changes over time, including periods without data. Boxplots of GIP estimates indicated systematic depletion and decreasing spatial variability, reflecting the impacts of degasification. Comparison with cumulative gas production from empirical well records showed approximately 85% agreement based on a relative similarity metric. Spatio-temporal GIP estimates were also used to estimate methane emissions to longwall ventilation systems and compared with reported emissions from the U.S. EPA Greenhouse Gas Reporting Program, showing similar distributions (≈80%) given data limitations. Overall, this integrated modeling approach provides time-dependent GIP estimates with broader implications for resource assessment applications.
Alabama↗
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↗
Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings
Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.
California↗
Multi-scale predictors of Northern Long-eared Bat (Myotis septentrionalis) occupancy in the United States
Historically, Myotis septentrionalis (Northern Long eared Bat) was among the most common forest-interior species in North America. Largely due to high mortality from white-nose syndrome, this species has experienced severe population declines across its range. To create an updated species distribution map representing summer occupancy probabilities from 2017 to 2022, we integrated stationary acoustic data with live-capture data from the database of the North American Bat Monitoring Program into a multi-scale, multi-method occupancy modeling framework. Our results provide data-driven predictions with quantified uncertainty for summer occupancy probabilities for Northern Long-eared Bats at 2 spatial scales across the range of the species, while also accounting for inherent observation biases (e.g., imperfect detection).
Journal of North American Bat Research↗
Imaging of seismic discontinuities using an adjoint method
For imaging of seismic discontinuities at depth, reverse time migration (RTM) is a powerful method to apply to recordings of seismic events. It is especially powerful when an extensive receiver array, numerous seismic sources, or both, permit adequate reconstruction of incident and scattered wavefields at depth. Reconstructing either the incident or scattered wavefield at depth becomes less accurate when relatively few recordings of seismic events are available. Here we explore an inverse scattering approach to imaging discontinuities based on an adjoint method, employing sensitivity kernels (Frechet derivatives) that represent jumps in material properties across seismic-discontinuity surfaces. When combined with ray-based requirements on scattering geometry, it constitutes a powerful approach to determining the locations and amplitudes of the discontinuities, recovering only those properties that can be resolved by a spatially limited source and/or receiver distribution. This is illustrated by synthetic examples with local sources followed by a field example in a subduction zone setting
Washington↗
Deep syntectonic burial of the Anthracite belt, Eastern Pennsylvania
Fluid inclusion microthermometry and Raman spectroscopy of fluid inclusions in quartz veins from the Pennsylvanian rocks of the Anthracite belt, eastern Pennsylvania support a deep burial model of coalification in favor of focused orogenic hot fluid flow. High-temperature (250 to 255 °C) trapping of CH 4 ± CO 2 saturated aqueous fluids and CH 4 ± CO 2 inclusions indicate fluid trapping at depths of 11.5 to 13.4 km under a cover of Pennsylvanian to Permian(?) syntectonic load. In the folded rocks to the south of the Anthracite belt, CH 4 ± CO 2 fluid inclusions indicate a sediment load that was up to 16.3 km thick. Re -equilibrated aqueous fluid inclusions from veins in Silurian through Devonian rocks give the same range of trapping conditions but a wide range of fluid salinities suggesting that folding, fracturing, and meteoric recharge resulted in the intermixing of fluids from throughout the stratigraphic succession.
Pennsylvania↗
Solute transport and modeling of water quality in a small stream
An injection of chloride, sodium, and stable strontium was made at a constant rate for 3 hours into Uvas Creek, Santa Clara County, Calif., to determine the mass transport processes in a small stream. Five observation points were selected within a 610-metre reach of the stream below the injection site. Water samples were collected at the observation points during and immediately after the injection. A mathematical model of the stream was obtained by solving analytically and optimally the one-dimensional mass transport equation of the solutes in the stream. Comparison of field results with a simplified mathematical model indicates the dominance of convection in the behavior of sodium and chloride. The concentration of chloride and sodium can be closely simulated by the model. However, strontium cannot be well represented by the simplified model, which contains a first-order decay-type sink.
California↗
Remote compositional analyses of space-weathered lunar maria
Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.
Planetary Science Journal↗
Inference of pattern-based geological CO2 sequestration and oil recovery potential in a commingled main pay and residual oil zone CO2-EOR flood
Several detailed studies have shown that residual oil zones (ROZs) can present significant resources for additional hydrocarbon recovery as well as subsurface carbon dioxide (CO 2 ) sequestration via enhanced oil recovery by injecting CO 2 (CO 2 -EOR). Field development strategies included new wells drilled dedicated to main pay zones (MPZ) and ROZs, or existing wells in MPZs deepened to ROZs for commingled injection-production using different well patterns. The latter presented a challenge when discerning the injection and production from each of the zones, and for subsequent quantification of CO 2 sequestration and EOR potential from different patterns and from the field. In this paper, an innovative method for analyzing commingled injections and productions from MPZs and ROZs, with application to pattern-based data from four staggered line drive patterns in Wasson Field's Denver Unit, Texas, USA, was developed. Decline curve and ratio-trend methods were used as means of history-matching and forecasting. Cumulative production-time and cumulative production-rate data for oil, gas, and water, as well as water-oil ratio (WOR) and gas-oil ratio (GOR), were analyzed along with injection data for time intervals covering major injection events in MPZ, or MPZ and ROZ combined. A combined analysis enabled inference of allocation of fluids into different zones during WAG (water alternating gas) injection and thereby estimation of CO 2 storage, utilization, and retention in different zones as a function of total injection. Results show that ROZs generally present higher CO 2 sequestration potential compared to MPZs, and a comparable incremental oil recovery factor of ∼20%, on average. Results based on ratio analysis further show that while the WOR trend of the pattern production is mostly dominated and controlled by ROZ, GOR is controlled by both intervals. Although the method relying on decline curves and the approach used in zonal fluid allocations are subject to their limitations, this study presents a practical and innovative well-pattern-based method to infer and forecast CO 2 sequestration and oil recovery quantities and fluid ratios from MPZs and ROZs in commingled operations and highlight the added potential offered by ROZs.
Texas↗
Statewide surficial geologic map of Nebraska underscores Quaternary landscape evolution from the High Plains to the Central Lowland
Surficial geologic mapping in Nebraska has been conducted primarily at the 1:24,000 scale since the mid-1990s, although there have also been limited efforts to map generalized Quaternary and surficial geology within the state and the region. We compiled and evaluated disparate maps (1:24,000–1:1,000,000) and datasets—including geologic, soil and soil parent material, and geomorphic maps as well as LiDAR derivatives—to produce a single 1:500,000 scale surficial geologic map of Nebraska that is the first of its kind. This new map provides a coarse-scale surficial geologic map that will be incorporated into a nationwide U.S. Geological Survey Quaternary geologic map. It also reflects the variation and uniqueness of physical landscapes in the state, where the Great Plains and Central Lowland physiographic provinces meet, further developing a richer interdisciplinary understanding of regional geomorphology in the heart of North America.
Nebraska↗
Preface to the focus section on earthquake-induced ground deformation
Ground failure due to earthquake-induced liquefaction, surface fault rupture, and landsliding has caused substantial damage in the past and continues to pose risks to critical infrastructure in the future. The consideration of ground failure is central to engineering analysis and design, ranging from site-specific engineering assessments of ground failure potential to broader regional-scale ground failure assessments of near-real-time earthquake impacts or longer-term portfolio risk. Uncertainties in all aspects of the problem complicate the prediction of earthquake-induced ground failure, including the identification of hazard-susceptible geomaterials, characterizing their spatial extent and physical properties, spatio-temporal variability in groundwater conditions, characterization of earthquake ground motions, assessing ground failure severity, and linking ground failure severity to infrastructure damage and loss.
Seismological Research Letters↗