Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “ACS ES&T Water”

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.

808 records · Page 11Linked to original sources
USGS2024

Cross section N–N' through the Valley and Ridge province of the southern Appalachian basin, from Greene County, west-central Alabama, to Bibb County, central Alabama

Introduction Geologic cross section N–N′ is the sixth in a series of geologic cross sections constructed by the U.S. Geological Survey to document and improve understanding of the geologic framework and petroleum systems of the Appalachian basin. Cross section N–N′ provides a regional view of the structural and stratigraphic framework of the Appalachian basin in the Valley and Ridge province in western and central Alabama; it spans approximately 69 miles (mi) (111 kilometers [km]). This geologic cross section is a companion to geologic cross sections E–E′ , D–D′ , C–C′ , I–I′ , and A–A′ that are located approximately 350 to 550 mi (563 to 885 km) to the northeast. Cross section N–N' complements earlier geologic cross sections through the Alabama part of the Appalachian basin. Although some of the other cross sections show more structural and stratigraphic detail, they are of more limited extent geographically and stratigraphically. Cross section N–N′ contains information that is useful for evaluating energy resources in the Appalachian basin. Although the Appalachian basin petroleum systems are not shown on the cross section, many of their key elements (such as source rocks, reservoir rocks, seals, and traps) can be inferred from lithologic units, unconformities, and geologic structures shown on the cross section. Other aspects of petroleum systems (such as the timing of petroleum generation and petroleum migration pathways) may be evaluated by burial history, thermal history, and fluid flow models based on what is shown on the cross section. In addition, cross section N–N′ may be used as a reconnaissance tool to identify plausible geologic structures and strata for the subsurface storage of liquid waste or for the sequestration of carbon dioxide.

Alabama↗
USGS2025

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗
USGS2026

Indirect mineral import reliance and provenance

Mineral commodity supply chain analyses rely on international trade data reported by individual countries as quantities of a mineral commodity form imported from (or exported to) a partner. However, export quantities frequently exceed a country’s domestic production, or occur when no production data are reported, suggesting that the trade partner is merely an intermediary in a transshipment. These discrepancies can result in misleading conclusions regarding supply chain vulnerabilities and dependencies. We present a two-stage methodology to reconcile gaps between reported material sources and actual producers. First, we construct trade networks for specific mineral forms, treating production as a type of import to distinguish producing nations from entrepôts. By tracing flows through these networks, we attribute a target country’s imports to original producers via both direct (in a single trade link) and indirect (transferring through intermediaries) pathways. Second, these production-attributed flows are incorporated into multi-stage supply chains to determine the upstream provenance of feedstock for domestic refining and processing. This approach provides a more representative picture of trade reliance. For example, while the United States (U.S.) Geological Survey reports no imports of unwrought antimony metal from Russia in 2022 (U.S. Geological Survey (2025). Mineral Commodity Summaries 2025. 10.3133/mcs2025), our analysis reveals that over 16% of U.S. imports can be traced back to Russian mining through intermediate processing in countries such as China, India, and Vietnam. Additionally, our analysis of the aluminum supply chain shows that while the U.S. is reported as 52% net import reliant on aluminum materials in 2022, it is 100% reliant on foreign bauxite, 7% of which arrived indirectly. This unreported reliance, which is predominantly tied to bauxite mined in Brazil (43%) and Jamaica (28%), highlights our methods ability to capture the supply chain’s dependence on foreign feedstock that may be missing in single-stage trade data.

Mineral Economics↗
USGS1998

Design, revisions, and considerations for continued use of a ground-water-flow model of the coastal plain aquifer system in Virginia

A digital numerical model of the ground-water-flow system in the Virginia Coastal Plain was constructed as part of the Regional Aquifer System Analysis program of the U.S. Geological Survey and subsequently revised based on reinterpretations of the hydrogeologic framework. The revised model has been incorporated by the Virginia Department of Environmental Quality as a means to evaluate the effects of existing and proposed ground-water withdrawals. The revised model, which is documented here, includes many design features of the original model. Principal modifications are (1) reassignment of some model layers to represent different aquifers, (2) replacement of applied recharge and the uppermost surface-water layer by constant-head cells to represent the unconfined aquifer, and (3) revision of the lateral extents and transmissivity and vertical leakance values assigned to model layers. Ground-water levels simulated by the revised model under steady-state conditions, using withdrawal rates representing the period 1978-80, generally are similar to the water levels simulated by the original model under transient conditions using the same withdrawal rates. Differences probably result largely from changes in transmissivity and vertical leakance values. T he revised model enables the continued evaluation of development of the large and complex aquifer system. Improved information on parts of some aquifers has been incorporated. Only large, regional-scale trends in ground-water water levels, flow directions, and flow rates are simulated, however, and not local trends and (or) short-term changes in water level and flow. Changing hydraulic stresses, knowledge of hydrogeologic conditions, and resource-management needs require periodic revision of the model to maintain its usefulness for the management of ground-water resources in the Virginia Coastal Plain.

Virginia↗
USGS2024

A partially nonergodic ground-motion model for Fourier amplitude spectra for the San Francisco Bay area, California, USA

We develop a partially nonergodic ground-motion model (GMM) for Fourier amplitude spectra for the San Francisco Bay Area, California, USA, using the Bayless and Abrahamson (2019) GMM as a reference ergodic GMM and developing location-dependent adjustments to the predicted median and variance. We compile regional ground-motion data from moment magnitude (𝑀 w ) >3 earthquakes occurring during 2000–2022 for which magnitude information is available in the U.S. Geological Survey Comprehensive Catalog (Guy et al., 2015). The data set predominantly consists of records from 𝑀 w 3.5–4.5 earthquakes but includes three well-recorded 𝑀 w > 5 events. Ground-motion residuals are evaluated using the time-averaged shear-wave velocity in the top 30 m (𝑉 S30 ) from the California-specific map of Thompson et al. (2018) and basin-depth site parameters from the seismic velocity model of Aagaard and Hirakawa (2021). The 𝑉 S30 dependence and basin-depth scaling of the reference ergodic GMM of Bayless and Abrahamson (2019) are evaluated and modified with the updated data set. We compute maps of site adjustments using a varying-coefficient model that considers the spatial correlation structure and uncertainties at each observation location. The spatial covariance model is developed using ground-motion residuals that are standardized by the uncertainty model, which allows for consideration of the aleatory variability in developing the site adjustments. The covariance model is fit considering the means and standard deviations of the site terms at all locations. The use of partially nonergodic median adjustments results in modified variance components of the within-event variability. Due to the low number of large-magnitude earthquakes that control seismic hazard in the data set, we do not modify between-event variance; however, we present adjustments to site-to-site variability for use in partially nonergodic hazard assessments.

California↗
USGS2025

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama↗
USGS2025

Beyond optimality: Dryland ecosystems infrequently use water efficiently for carbon gain

Optimality theory assumes plants maximize carbon gain per unit water lost and is often implemented to scale leaf-level carbon gain and water use to regional and global scales. Optimality theory is often mathematically represented by assuming plant water-use efficiency (WUE) scales with VPD k , where k = ½ represents expected optimal behavior. It is unclear, however, if this relationship holds in arid and semi-arid ecosystems that are strongly impacted by soil and atmospheric moisture status. We used data from seven flux tower sites along an aridity gradient in New Mexico to answer: how does the relationship between WUE and VPD compare to expectations based on optimality theory? To address this question, we integrated the Dynamic Evapotranspiration Partitioning Approach for Rapid Timescales with a stochastic antecedent model to estimate ecosystem-level WUE (GPP/T) and the net sensitivity of WUE to VPD, or k Dynamic , which we compare to the theoretical optimal sensitivity of k = ½. Our results show that optimality theory is not always appropriate, and k Dynamic often deviates from ½, especially at some of the more arid sites or during periods of low soil moisture. At less arid, higher elevation sites, k Dynamic is most consistent with optimality theory at moderate VPD levels, but not at high VPD. In general, the sensitivity of WUE to VPD is highly variable such that k Dynamic exhibits notable daily and seasonal variability, suggesting highly dynamic stomatal behavior. These results emphasize that representing plant water-use strategies as dynamic in time and space is critical to improving large-scale estimates of plant water use.

New Mexico↗
USGS2024

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https:​//ca.water​.usgs.gov/​bay-​delta/​2018-​delta-​wide-​mapping-​surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California↗
USGS2025

Geologic models underpinning the 2024 U.S. Geological Survey assessment of undiscovered oil and gas resources in the Hosston and Travis Peak Formations of the onshore Gulf Coast region, U.S.A.

The Early Cretaceous (Berriasian–Hauterivian) Hosston Formation in Louisiana and eastward is time correlative to the Travis Peak Formation of Texas and southern Arkansas. The formation is a first-order clastic sequence with a regional carbonate transgressive surface as an upper contact. The Hosston and Travis Peak formations contain conventional natural gas and oil accumulations that have been produced for nearly a century. These mature reservoirs contain terrigenous fluvial-deltaic, shore-zone, and paralic deposits across the productive trend; organic-lean mudstone and siltstone lithologies are found outboard of the Lower Cretaceous shelf margin. Producing reservoirs exhibit normal pressure gradients from 0.43 to 0.55 psi/ft (9.7 to 12.4 kpa/m), depths from 4000 to over 20,000 ft (1220 to 6100 m), and temperatures from 150 to 385°F (65 to 196°C). Wells are primarily vertical completions. The number of new field wildcats has been declining since the late 1990s. This paper presents comprehensive geologic models, which include lithofacies maps, structure and isopach maps, burial history models, regional seismic interpretations, and events charts that underpin the recently completed U.S. Geological Survey assessment of undiscovered, technically recoverable hydrocarbons within the Hosston and Travis Peak formations. This study also provides geographic and stratigraphic distributions of Hosston–Travis Peak reservoir properties, including geopressure, reservoir temperature, porosity, permeability, API gravity, and gas-oil ratios. Results indicate estimated undiscovered, technically recoverable mean resources of 28 million barrels of oil and 35.8 trillion cubic ft of gas in conventional and continuous accumulations within the Lower Cretaceous Hosston and Travis Peak formations of the onshore U.S. Gulf Coast region. Quantitative assessment results are detailed in U.S. Geological Survey Fact Sheet 2025–3021 and associated Data Release.

GCAGS Journal↗
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Integration of rupture directivity models for the US National Seismic Hazard Model

Several rupture directivity models (DMs) have been developed in recent years to describe the near-source spatial variations in ground motion amplitudes related to propagation of rupture along the fault. We recently organized an effort towards incorporating these directivity effects into the USGS National Seismic Hazard Model (NSHM), by first evaluating the community's work and potential methods to implement directivity adjustments into probabilistic seismic hazard analysis (PSHA). Guided by this evaluation and comparison among the considered DMs, we selected an approach that can be readily implemented into the USGS hazard software, that provides an azimuthally varying adjustment to the median ground motion and its aleatory variability. This method allows assessment of the impact on hazard levels and provides a platform to test the DM amplification predictions using a generalized coordinate system, necessary for consistent calculation of source-to-site distance terms for complex ruptures. We give examples of the directivity-related impact on hazard, progressing from a simple, hypothetical rupture, to more complex fault systems, composed of multiple rupture segments and sources. The directivity adjustments were constrained to strike-slip faulting, where DMs have good agreement. We find that rupture directivity adjustments using a simple median and aleatory adjustment approach can impact hazard both from a site perspective and on a regional scale, increasing shaking off the end of the fault trace up to 30--40\% and potentially reducing it for sites along strike. Statewide hazard maps of California show that the change in shaking along major faults can be a factor to consider for assessing long-period (>ls) near-source effects within the USGS NSHM going forward, reaching up to 10--20\%. Finally, we suggest consideration of minimum parameter ranges and baseline requirements as future DMs are developed to minimize single approach adaptations, to enable more consistent application within both ground motion and hazard studies.

Earthquake Spectra↗
USGS2026

Tringa flavipes (Lesser Yellowlegs) from separate breeding sites subdivides the Prairie Pothole Region in space and time during southbound migration

Some staging regions support multiple groups of the same migratory species, each of which may use the region differently. Characterizing the ways, in which separate groups use such regions can therefore help to identify vulnerabilities during this sensitive period of the annual cycle. The Prairie Pothole Region (PPR) is a massive wetland complex in the northern Great Plains of North America used by ∼11 million shorebirds during migration. The PPR has been heavily modified by agriculture and is experiencing varied effects of global climate change, threatening the health of the shorebirds that rely on it. Here, we used 6 seasons of southbound tracking data of Tringa flavipes (Lesser Yellowlegs)—a long-distance migratory shorebird species with an estimated population decline of 63% over the last 4 decades—from 9 sites across their breeding range to explore differences in migratory behavior within this important staging region. We found that 75% of tracked individuals used the region during southbound migration, and T. flavipes from different breeding sites detoured 110–875 km from their most direct migratory route to access the PPR. Individuals that arrived later stayed longer and made more stops within the region than those that arrived early. Individuals originating from different breeding sites also displayed spatial and temporal segregation within the region: T. flavipes from southwest and central Alaska relied heavily on the northwestern PPR, while those from Canada used the central and southeastern portions of the PPR. Finally, timing of use varied among groups, but the southeastern PPR became increasingly important over the course of the southbound migratory window, as other wetlands likely dried out. Our study highlights the portions of the PPR of critical importance to migrating T. flavipes and the diversity of ways, in which different groups from within the same species can use a single staging region.

Alaska, Manitoba, Northwest Territories, Ontario, ↗
USGS2026

The 2025 Puerto Rico and Virgin Islands U.S. National Seismic Hazard Model Update: Ground motion model selection and comparison

We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.

Puerto Rico, Virgin Islands↗
USGS2026

Cambrian and Ordovician stratigraphy, conodont biostratigraphy, and microfacies analysis to support 1:24,000-scale geologic mapping of the southern Lake Champlain valley, New York and Vermont

Introduction Geologic mapping in the southern Lake Champlain valley of New York and Vermont (fig. 1) has required evaluation of stratigraphic nomenclature used since the early 1900s. The paleogeography of the Cambrian and Ordovician adjacent to the Adirondack Highlands, the high-relief mountains of Proterozoic igneous and metamorphic rocks of northern New York, has impacted the lithostratigraphic succession, facies changes, and several hiatuses that occur within the package of clastic and carbonate rocks that demonstrates overall deepening consistent from the southern to northern Appalachians. The stratigraphic nomenclature used for the geologic mapping of the Cambrian and Ordovician strata of the study area is derived from evaluation of names put forth by previous workers from areas of northern New York and western Vermont. Descriptions of type localities or type sections and other criteria set forth by the North American Stratigraphic Code (North American Commission on Stratigraphic Nomenclature [NACSN], 2021) were consulted. The North American Stratigraphic Code states that lithostratigraphic units, such as the fundamental unit, the formation, be defined by lithic character and be mappable (NACSN, 2021, Article 24). However, some early workers (Walcott, 1912; Kay, 1937; Fisher and Hanson, 1951; Oxley and Kay, 1959) have defined stratigraphic units based on fossil content, which may or may not be facies dependent, and hence may not be mappable as lithostratigraphic units. Although preservation of older names takes priority in stratigraphic nomenclature (NACSN, 2021, Article 7), some do not have adequately defined stratotypes or are located a distance away with different lithologies than that in the study area. Biostratigraphy can be a helpful tool for understanding correlations and facies changes. Herein, conodont biostratigraphy is utilized to date and correlate separate lithologic units. Further, microfacies analyses of units through inspection of petrographic thin sections provide paleogeographic information that helps in understanding the genesis of geologic units.

New York, Vermont↗
USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗
USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗
USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗