Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “AGU Advances”

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.

101 records · Page 4Linked to original sources
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↗
USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗
USGS2025

A review of post-wildfire adaptations of surface-water-quality models: Synthesis, gaps, and opportunities

As wildfires increasingly affect water-supply watersheds, the demand for models to predict water-quality responses is increasing. This work reviews and synthesizes existing post-wildfire applications of water-quality models in the context of geographic and ecohydrological distribution, hydrologic and water-quality response process representation, model parameterization, model and input data scales, model calibration data availability, as well as calibration and performance evaluation approaches. Emphasis is placed on models that simulate water-quality output, rather than sediment and erosional response as the primary focus. Here, identified gaps and opportunities to advance the post-wildfire application of water-quality models include: 1. applying models in under-represented geographic and ecohydrologic regions, 2. simulating multiple streamflow generation mechanisms, including groundwater, with an emphasis on shifting dominant flow pathways as the landscape recovers following wildfire, 3. adding studies that include the simulation of metals, 4. incorporating more biogeochemical and in-stream processes to model applications, 5. applying finer spatial and temporal resolution of precipitation data input as well as finer spatial resolution hydrologic response units, 6. implementing fully distributed grid or element models or finer resolution response units to capture burn severity heterogeneity, 7. collecting enhanced water-quality data for model calibration and validation, 8. conducting model-intercomparison studies, and 9. developing model parameter value guidance in post-wildfire applications. These identified gaps and opportunities may assist users in deciding on key processes and approaches to consider in modeling post-wildfire water-quality conditions.

Science of the Total Environment↗
USGS2026

Regional variations in sea ice and primary productivity in the Bering Sea during Marine Isotope Stage 11

Marine Isotope Stage (MIS) 11 (424-374 ka) has long been an analog for Holocene climate, because it is the most recent interglacial period with similar orbital conditions. However, there is significant global and regional variability in the climate response to MIS 11 warmth. Here, we review sediment core records from across the Bering Sea to investigate changes in paleoceanographic conditions during Marine Isotope Stages 12-10. Sea ice was present over much of the Bering Sea during MIS 11, but today, none of the sites investigated are ever ice-covered. This suggests that sea ice regimes in the Bering Sea during MIS 11 were different to those of the Holocene. There are also regional differences in the response of sea ice to MIS 11 warming. At the Umnak Plateau, Southeastern Bering Sea, sea ice concentrations decline during deglaciation, but they remain high at the slope sites until Peak MIS 11. Sea ice re-advances over the Umnak Plateau during peak interglacial warmth, at the same time that it declines over the slope sites. Late MIS 11 is characterized by high concentrations of seasonal sea ice at the Umnak Plateau, whilst sea ice at the slope sites fluctuates between consolidated and unconsolidated ice cover. This east-west dichotomy may be explained by changes in the behavior of the Aleutian Low. Productivity increases dramatically during deglaciation due to increased upwelling and sea level rise bringing fresh nutrients into the oceans. This is characterized by increased diatom productivity, increased terrestrial carbon deposition, and laminations at all sites.

Alaska↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2025

Effect of land cover type on 3D deformation recovery from synthetically deformed high-resolution satellite optical imagery

The limits of detection for earthquake surface deformation in the spatial domain have improved with advances in remote sensing imagery data availability, resolution, and analysis. Sub-pixel correlation and digital elevation model (DEM) differencing from sub-meter, earthquake-spanning satellite optical imagery has enhanced surface rupture mapping and deformation measurements. However, knowledge of measurement accuracy and uncertainty is limited. To address this, we construct orthophotos and digital elevation models (DEMs) from repeat high resolution (∼0.5 m) satellite optical imagery along two sections of the Garlock fault, California with clear fault geomorphology and differing land cover. We deform later sets of DEMs and images with synthetic earthquakes containing both diffuse and discrete horizontal and vertical displacements. Sub-pixel image correlation and DEM differencing demonstrate how vegetation degrades recovered displacement accuracy. In barren land cover, horizontal displacements are detectable to an expected ∼1/10th-pixel size. With shrubs, trees, and grass, detectable displacements increase to >1/2-pixel size, and filtering results by correlation score and using elevation values as input rather than image values improves accuracy. Vertical displacement detection thresholds remain lower in vegetation, at >1-pixel size. Higher slope angles degrade displacement recovery, worsened by vegetation. Diminishing seasonal separation improves accuracy over vegetated regions, though not to the level achieved in barren environments. These results will inform research and operational efforts on the utility of high resolution satellite optical imagery for detecting deformation in varied land cover. Furthermore, they reveal where alternative measurements, such as from LiDAR or radar interferometry, are required to mitigate the effects of vegetation and capture fine-scale crustal deformation.

California↗
USGS2026

A review and synthesis of post-wildfire shifts in hydrologic processes and streamflow generation mechanisms

Critical water supply watersheds in the western United States (WUS) are impacted by wildfires, with potential negative effects on water quality and quantity. Scientific understanding is currently insufficient to deliver estimates of wildfire consequences for water quantity that are regionally accurate. Regional variability in the directionality and magnitude of post-wildfire shifts in streamflow generation fuels uncertainty in estimates of wildfire effects on water supply. In this work we provide a narrative review of wildfire effects on hydrologic processes and the resulting changes in streamflow generation mechanisms with a focus on the WUS, incorporating other global regions when pertinent. A conceptual model summary of wildfire effects on streamflow generation emphasizes: (1) precipitation seasonality, (2) synchrony of precipitation and potential evapotranspiration, (3) net shifts in interception, evaporation, and transpiration relative to total annual precipitation, (4) vegetation changes, including compensatory uptake and type conversion, (5) degree of overlap in rainfall rates and infiltration, (6) fire extent and severity, (7) burn scar positioning (e.g. in headwaters or proximal to watershed outlet), (8) scale-dependent groundwater leakage, (9) near-surface water storage reduction, and (10) soil to groundwater connectivity. Ongoing gaps and challenges include separating the influences of precipitation variability, water withdrawals, and post-fire land management; compound and overlapping disturbances; and lack of pre-fire data. Notable future opportunities include: harnessing ever-improving gridded and remotely sensed precipitation and fire-effects data; linking geophysical, isotopic tracer, and geochemical signatures to diagnose hydrologic changes; leveraging physically based and data-driven model advancements; and analyzing streamflow generation recovery trajectories across diverse watersheds.

western United States↗
USGS2026

Quantifying methane emissions from a rich fen with uncrewed aircraft systems in boreal Alaska

Thawing of permafrost in northern latitudes is accelerating, potentially releasing substantial amounts of methane (CH 4 ) as forested permafrost plateaus transition into wetlands. This ecosystem shift alters the carbon exchange between the soil and atmosphere, influencing the permafrost-carbon feedback. Monitoring these changes may require measurement platforms operating across varied spatial and temporal scales. Recent advancements in small uncrewed aircraft systems (sUAS) enable high resolution CH 4 flux quantification in remote, complex terrains; however, comparisons with established methods such as eddy covariance flux towers remain limited. We used a hexacopter sUAS to quantify CH 4 emissions from the Alaska Peatland Experiment, a wetland within the Bonanza Creek Experimental Forest. Using an ensemble of methods to define the background CH 4 concentration, along with near surface emissions from soil chambers, helped constrain our flux estimates. The sUAS method yielded an average flux of 0.0077 ± 0.0019 mol s −1 CH 4 , within a factor of two concurrent tower-derived total source flux estimates (0.0036 ± 0.00042 mol s −1 CH 4 ). To assess spatial drivers of observed fluxes, we conducted a 2D footprint analysis and overlaid the results with high-resolution hyperspectral land cover classification, quantifying vegetative contributions within each footprint. This revealed higher fen representation in sUAS measurements (73.8%) than in tower footprints (58.8%), and lower tussock meadow representation (15.6% and 30.3%, respectively). These differences were consistent with known variation in vegetation-specific CH 4 emissions. Our results highlight that combining footprint modeling with land cover characterization can enhance interpretations of CH 4 fluxes and guide cross-platform comparisons.

JGR Atmospheres↗
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗
USGS2025

Landsliding follows signatures of wildfire history and vegetative regrowth in a steep coastal shrubland

Five years after the deadly and destructive 9 January 2018 Montecito debris flows (Santa Barbara County, California, USA), an atmospheric river storm on 9 January 2023 triggered widespread landsliding that affected many of the same drainages in the Santa Ynez Mountains. Using high-resolution aerial imagery, we identified >10,000 landslides over an ∼160 km 2 area. Most of the landslides were shallow (<1 m in depth) translational debris slides that initiated on steep (∼40°), south-facing hillslopes, with the highest incidence of landsliding in a sandstone-dominated bedrock unit. The landslides mobilized into debris flows and delivered substantial quantities of sediment downstream, which contributed to costly infrastructure impairments. We detected order-of-magnitude differences in landslide density across the study area that could not be attributed to variations in geomorphology (topographic aspect and slope), geology (bedrock type), or hydrology (seasonal antecedent rainfall, peak hourly storm rainfall intensity, total storm rainfall), which are usually considered relevant factors for shallow landsliding. Rather, we found that vegetation regrowth following wildfire was likely a relevant factor associated with the highly variable landslide densities. Hillslopes with less and different types of vegetation regrowth after fire appear to have been more susceptible to shallow landslides. We identify a possible vegetation control on postfire landsliding, which highlights an opportunity for hypothesis testing using more advanced techniques to track the evolution of vegetation cover and vegetation type in steep shrubland environments following wildfire.

California↗
USGS2025

Nanometer-scale relationships between sedimentary organic matter molecular composition, fluorescence, cathodoluminescence, and reflectance: The importance of oxygen content at low thermal maturities

Molecular characterization of sedimentary organic matter (SOM), termed macerals, is a common goal when seeking to understand petroleum generation as well as other geologic processes in deep time. However, unambiguous measurement of discrete macerals is challenging due to the small size of organic particles in sedimentary rocks, the proximity of different organic matter types to one another, mineral-organic matter interactions, and maceral mixing that occurs during SOM isolation prior to ex situ analysis. The recent advent of infrared spectrometers capable of nanometer-scale resolution and the application of these technologies to geologic samples has enabled advances in rapid, in situ molecular characterization of SOM allowing for insights into paleoenvironmental processes, such as organic matter productivity and preservation, among others. Here we employ one such technology, optical photothermal infrared (OPTIR) spectroscopy, to map SOM functional group distributions at 500-nm resolution in a sample from the Lower Cretaceous Sunniland Limestone of the South Florida Basin. Examined fields of view include occurrences of amorphous organic matter (AOM), inertinite, micrinite, solid bitumen, telalginite, and vitrinite. OPTIR data from these macerals are compared against traditional organic petrographic data from the same organic grains including fluorescence intensity and white light reflectance as well as against cathodoluminescence response, an emerging organic petrographic approach. Maceral oxygen content (using carbonyl functional group abundance as a proxy) is observed to vary widely between maceral types but correlates strongly with fluorescence and cathodoluminescence intensity as well as against reflectance. These findings highlight the important role that oxygen content plays in determining the optical properties of SOM and further demonstrate the ability of OPTIR to discriminate subtle molecular differences between SOM types.

Organic Geochemistry↗
USGS2026

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗
USGS2025

Reflections on a trio of North American earthquakes in 1925

In 1925, three moderately large damaging earthquakes occurred in North America over four months: the 28 February (local time; LT) M 6.2 Charlevoix, 27 June (LT) M 6.6 Montana, and 29 June M 6.5 Santa Barbara earthquakes. The centennial anniversaries of these events motivated this retrospective consideration focused on the ground motions generated by the three events, including a reconsideration of early intensity assignments for the Montana earthquake. At the time, these three earthquakes appeared to support the arguments of some geologists who downplayed the severity of seismic hazard in southern California relative to other parts of the country. Some of the arguments advanced at that time, for example that Los Angeles “has the least to fear from ‘Acts of God’ of any city under the American flag,” ( Hill, 1928 ) sound naïve if not laughable now, but a comparison of well‐constrained shaking distributions for the three earthquakes reveals the dramatic difference in wave propagation efficiency in western versus eastern North America (ENAM), which leads to moderate ENAM events being felt to much larger distances. At M 6.2, the 1925 Charlevoix earthquake was a notably large event in ENAM. This earthquake was the largest event in eastern Canada since 1870 and caused damage in the epicentral region in addition to towns as far away as 200 km, with felt shaking extending over 1000 km. In contrast, felt shaking from the Santa Barbara earthquake barely extended beyond ∼200 km. Compiling published intensity distributions for larger ENAM earthquakes, we show that perceptible earthquake shaking is not uncommon in ENAM over century time scales, but experience with weakly felt shaking may incline people to downplay potential earthquake risk.

California, Montana, Quebec↗
USGS2026

Lithium brine potential of North Dakota: Results from the 2025 sampling program

Lithium (Li), classified as a critical mineral and key in energy storage applications because of its high energy density, faces unprecedented demand. This has driven interest in unconventional sources of Li, especially brines from sedimentary basins. Accordingly, the spatial distribution of basinal brinehosted Li resources across North America is becoming better defined, including the Williston Basin underlying Saskatchewan, Canada. This has led to growing investment and exploration by industry throughout the 2020s; however, despite this interest, the Li potential of basinal brines of the Williston Basin in North Dakota remain underexplored. This work presents the results from the 2025 joint brine sampling campaign between the North Dakota Geological Survey, the United States Geological Survey (USGS), and The University of Texas at Austin. Twenty-seven fluid samples were collected from producing oil and water source wells from eight stratigraphic intervals across western North Dakota and analyzed for major cation and anion concentrations. Similar to southeastern Saskatchewan, the highest Li concentrations were identified in the Devonian Duperow Formation, with concentrations up to 179 mg/L being observed in a well commingled between the Duperow and Red River Formations. Some samples from the overlying Birdbear Formation exceeded 75 mg/L, whereas Li concentrations were lower in the Madison Group (33-44 mg/L) and the Three Forks (50 mg/L), Dawson Bay (43 mg/L), Winnipegosis (23 mg/L), and Red River (40.0 mg/L) Formations. This work indicates Li resource potential exists across the Williston Basin and advances our understanding of the Li distribution in basinal brines which could provide a new domestic source of this critical mineral.

North Dakota↗
USGS2026

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗
USGS2025

The 3D National Topography Model Call for Action—Part 2: The Next Generation 3D Elevation Program

The three-dimensional (3D) National Topography Model initiative to integrate elevation and hydrography data includes the next generation of hydrography data from the 3D Hydrography Program and the next generation of elevation data from the 3D Elevation Program (3DEP). The first-ever collection of light detection and ranging (lidar) data for the nation (IfSAR for Alaska) provides a critical baseline reference, and the addition of multiple repeat elevation mapping projects as part of the next generation of 3DEP would substantially expand analysis capabilities. As the U.S. Geological Survey (USGS) is closing in on our goal of complete coverage with 3DEP data available or in progress for 98.3 percent of the Nation at the end of fiscal year 2024, the USGS is already transitioning to the next generation of 3DEP. Based on the 3D Nation Study results and input from a broad range of stakeholders, the USGS National Geospatial Program has finalized a new design for 3DEP that provides increased lidar quality levels and refresh rates. The new program is designed with more flexibility to meet changing user needs and take advantage of improvements in mapping technologies. The program will aim to expand the level of interagency coordination for topobathymetric lidar acquisition for inland rivers. The next generation of 3DEP will also aim to emphasize research, including advancing program design, products, and services and engaging and leveraging the evolving 3D industry. Research goals also include becoming more flexible in meeting user needs and taking advantage of evolving remote-sensing technologies. The program also plans to move from focusing on producing standard products to producing a concept of a 3D Nation Ecosystem with a variety of inputs, products, and services.

Circular↗