USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal for Disaster Risk Reduction”

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.

1,443 records · Page 45Linked to original sources
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2026

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

Reconciling bias in moderate magnitude earthquake ground motions predicted by numerical simulations

Recent studies found a significant underprediction in ground motion intensity measures for finite-fault simulations of moderate magnitude events in southern California relative to established ground motion models. This study aims to understand the source(s) of this bias by evaluating ground motion residuals. For this, simulations have been performed for a total of 27 well-recorded earthquakes in southern California. Systematic efforts have been employed to identify the source(s) of bias by ruling out factors that are insignificant. Preliminary findings indicate that the magnitude-area scaling used in the simulations is the likely major cause of the observed bias. Adjustment in the source attributes on event-by-event basis is underway to study if the observed bias can be reconciled.

Conference Paper↗
USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2024

Bioconcentration of per- and polyfluoroalkyl substances and precursors in fathead minnow tissues environmentally exposed to aqueous film-forming foam-contaminated waters

Exposure to per- and polyfluoroalkyl substances (PFAS) has been associated with toxicity in wildlife and negative health effects in humans. Decades of fire training activity at Joint Base Cape Cod (MA, USA) incorporated the use of aqueous film-forming foam (AFFF), which resulted in long-term PFAS contamination of sediments, groundwater, and hydrologically connected surface waters. To explore the bioconcentration potential of PFAS in complex environmental mixtures, a mobile laboratory was established to evaluate the bioconcentration of PFAS from AFFF-impacted groundwater by flow-through design. Fathead minnows ( n = 24) were exposed to PFAS in groundwater over a 21-day period and tissue-specific PFAS burdens in liver, kidney, and gonad were derived at three different time points. The ∑PFAS concentrations in groundwater increased from approximately 10,000 ng/L at day 1 to 36,000 ng/L at day 21. The relative abundance of PFAS in liver, kidney, and gonad shifted temporally from majority perfluoroalkyl sulfonamides (FASAs) to perfluoroalkyl sulfonates (PFSAs). By day 21, mean ∑PFAS concentrations in tissues displayed a predominance in the order of liver > kidney > gonad. Generally, bioconcentration factors (BCFs) for FASAs, perfluoroalkyl carboxylates (PFCAs), and fluorotelomer sulfonates (FTS) increased with degree of fluorinated carbon chain length, but this was not evident for PFSAs. Perfluorooctane sulfonamide (FOSA) displayed the highest mean BCF (8700 L/kg) in day 21 kidney. Suspect screening results revealed the presence of several perfluoroalkyl sulfinate and FASA compounds present in groundwater and in liver for which pseudo-bioconcentration factors are also reported. The bioconcentration observed for precursor compounds and PFSA derivatives detected suggests alternative pathways for terminal PFAS exposure in aquatic wildlife and humans.

Environmental Toxicology and Chemistry↗
USGS2026

Rapid earthquake magnitude classification via P-wave strains from borehole strainmeters and Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) offers a promising approach for earthquake early warning (EEW) in settings where seismic networks are costly to maintain. By repurposing fiber-optic cables as dense strainmeter arrays, DAS enables real-time earthquake detection wherever those fibers are accessible. However, poor azimuthal coverage and challenges in estimating magnitude from strain measurements remain key hurdles in applying for earthquake monitoring. Here, we develop a machine learning method to distinguish large (M≥5.4) earthquakes from smaller ones within the first 4 seconds of a strain waveform after a P-wave arrival without determining location. Using ensemble decision tree models trained on borehole strainmeter data (3.5≤M≤7.1) and tested on onshore DAS waveforms (including the 2024 M7 Offshore Cape Mendocino earthquake), we find that low-frequency (0.2–0.5 Hz) continuous wavelet transform coefficients are the strongest predictors of magnitude, in addition to strain amplitude. Both DAS and borehole strainmeters effectively capture long-period strain signals, making these findings valuable for EEW systems. Our method shows high precision compared to the real-time EEW system, ShakeAlert®, supporting the position that DAS is a viable technology for earthquake monitoring and magnitude classification.

California↗
USGS2026

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri↗
USGS2026

Simulation of groundwater flow in the Silurian aquifer, eastern Iowa, 2020–45

The Silurian aquifer is an important water source for municipalities, industry, and rural households and communities in eastern Iowa, including Johnson County. Increasing demand for groundwater from the aquifer indicated that analytical tools may be beneficial to quantify groundwater resources and inform water-management decisions for the aquifer. The U.S. Geological Survey, in cooperation with the Johnson County Board of Supervisors, developed conceptual and numerical groundwater models to simulate water levels in the Silurian aquifer, determine groundwater budgets, and forecast changes in groundwater levels through 2045. A MODFLOW numerical model was constructed to match water levels for 2020 through 2022. Model performance was assessed by comparing simulated and observed water-level hydrographs at several pumping and monitoring wells. Hydrograph comparison indicates that the model simulates water levels accurately at monitoring wells that are not near areas of pumping and generally overpredicts water levels near pumping wells, while accurately simulating the drawdown trend over time at pumping wells. This transient model framework was then used to predict water levels through 2045 in a scenario of moderate drought and increased groundwater withdrawals. Results from the predictive model indicate as much as 13 meters of additional drawdown in 2045 from 2020 water levels in the area of pumping wells.

Iowa↗
USGS2025

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2024

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2024, shaded depictions of estimated water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 301 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 184 ft above NAVD 88. The largest decline in water-level altitudes depicted by the 1977–2024 long-term water-level-change map was in south-central Montgomery County. In comparison, the 1990–2024 long-term water-level-change map depicts the largest declines in water-level altitudes in an area northwest of The Woodlands and in an area of northern Waller County. The largest rise in water-level altitudes for 1977–2024 is depicted in an area of east-central Harris County, whereas the largest rise in water-level altitudes for 1990–2024 is depicted in an area of central Harris County. The 5-year short-term water-level-change map depicts the largest declines in several parts of the study area, but these declines are concentrated primarily in northern Fort Bend County, southwestern Harris County, and south-central Montgomery County. The largest rise for 2019–24 is depicted at a well in northern Fort Bend County. The 1-year short-term water-level-change map depicts the largest declines at a well in northern Fort Bend County and a well in west-central Harris County. The largest rise for 2023–24 is depicted at a well in east-central Fort Bend County. In 2024, shaded depictions of estimated water-level altitudes for the Jasper aquifer ranged from about 255 ft below NAVD 88 to about 321 ft above NAVD 88. The 2000–24 long-term water-level-change map depicts the largest water-level decline in an area of central San Jacinto County; the largest rise is depicted in an area of central Grimes County. The 5-year short-term water-level-change map depicts the largest declines across parts of central and southern Montgomery County and at one well in north-central Harris County. The largest rise for 2019–24 is depicted at a well centered on the Montgomery-Grimes County line. The 1-year short-term water-level-change map depicts the largest declines at two wells in south-central Montgomery County and one well in northwestern Montgomery County on the west side of Lake Conroe. The largest rises during 2023–24 are depicted at one well in northwestern Montgomery County and one well in south-central Montgomery County.

Texas↗
USGS2025

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

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

Seismological Research Letters↗
USGS2026

Water use permits as of July 2024 and reported water use near the North Unit of Theodore Roosevelt National Park, North Dakota, 1980–2023

Starting in the early 2000s, increasing oil and gas development in western North Dakota created a need for additional water resources from surface-water and groundwater sources near the North Unit of Theodore Roosevelt National Park. To summarize the use of water in that area, the U.S. Geological Survey, in cooperation with the National Park Service, developed a map of surface-water and groundwater resources, aquifers, and water-use diversions, and plotted water-use trends from 1980 to 2023. Reported water used from permits in the map area has more than doubled since 2020, increasing from about 750 acre-feet in 2020 to about 2,300 acre-feet in 2022 and 2,000 acre-feet in 2023. Surface water provided the primary source of reported water used for the study period with an average of about 410 acre-feet per year from 1980 through 2017 and about 1,330 acre-feet per year from 2018 through 2023. After 2011, groundwater sourced from the Little Missouri River, Tobacco Garden Creek, Fox Hills, Fort Union, and Dakota aquifers became a larger portion of total annual reported water use from permits in the map area. From 1980 through 2015, water use for irrigation averaged 86 percent of the total annual reported surface-water and groundwater use in the map area. Starting in 2011, however, industrial uses became a proportionally larger total use of water, and in 2015, became the highest reported volume of water use in the map area. From 2011 to 2023, industrial use designated for water depots increased from 50 acre-feet to about 1,370 acre-feet, accounting for about 70 percent of total reported water use in the map area in 2023.

North Dakota↗
USGS2026

Production of mineral commodities and geospatial map of the mineral industries and related infrastructure of China

As part of the U.S. Geological Survey’s (USGS) mission to distribute global mineral information and analyze supply chains, this study provides a comprehensive review of the global significance of China’s mineral production and capacity in 2023. Of 77 mineral commodities in the USGS dataset, China produced 74 and was the world’s first-ranked producer for 39 of the 74. Compared to the high share of global mineral production, including up to 98 percent of global gallium production, the country’s share of global mineral reserves was relatively small, ranging from 20 percent (zinc ore) to 52 percent (tungsten ore). China’s imports of metal ores, slag, and ash accounted for 64 percent of global imports of such commodities by value. The country’s exports of base metals and articles of base metal accounted for 17 percent of the global exports. To help nongeographic information system users assess the spatial distribution of mineral mines, processing facilities, and ports for trades in China, this study created a geospatial (also called “georeferenced”) portable document format (GeoPDF) map. In addition, the GeoPDF contains mineral resource tracts (such as antimony, copper, potash, coal, and oil and gas), exploration sites, and energy infrastructure based on the preexisting USGS data.

Open-File Report↗
USGS2026

Age and colony variation in Adélie penguin metapopulation vital rates: Insights from a 25-year mark–recapture study

Understanding how vital rates vary with age, life-history stage, and among populations is fundamental for predicting the demographic consequences of environmental change, especially in longer-lived species with complex life histories. These species often exhibit delayed maturity and iteroparity, making nuanced demographic insights critical for assessing their long-term viability. This study investigated age- and colony-related variation in Adélie penguin vital rates including survival, recruitment, and breeding propensity. We used mark–recapture data collected over 25 years (1996–2020) from three Adélie penguin breeding colonies that differed in population sizes and trends but comprised a metapopulation located at capes Royds, Bird, and Crozier on Ross Island, Antarctica. We used multi-state models to estimate survival and detection rates relative to reproductive state and breeding colony and estimated transition probabilities reflecting movements between reproductive states and colonies. Apparent survival varied by reproductive state, colony, and age and averaged 0.80 (SD = 0.02) at Bird, and 0.72 (SD = 0.03) and 0.73 (SD = 0.03) at Crozier and Royds, respectively, for pre-breeders age 2–7 years with strong declines in pre-breeder survival after age 8. We observed less age-related variation in survival of breeders and non-breeders, but we observed differences between colonies with lower survival for breeders (0.72 to 0.80) compared to non-breeders (0.75 to 0.82). The average probability of surviving the first 2 years after fledging ranged from 0.43 (SD = 0.14) at Royds and Crozier (0.43, SD = 0.10) to 0.55 (SD = 0.15) at Bird. Movement between colonies was highest for pre-breeders (0.00%–12.00% depending on age and colony) and lowest for breeders (<0.20%). We observed the lowest age-related recruitment rates at Royds, with recruitment at Crozier almost twice as high, and intermediate at Bird. Breeding propensity was highest at Crozier and lowest at Bird. Colony-specific variation in vital rates likely contributed differently to population trajectories, suggesting that care must be taken to extrapolate vital rate estimates across colonies even within a metapopulation. These findings also highlight the importance of considering age, life-history stage, and geographic variation when assessing population vital rates in species with complex life histories.

Frontiers in Ecology and Environment↗
USGS2024

Debris-flow entrainment modelling under climate change: Considering antecedent moisture conditions along the flow path

Debris-flow volumes can increase along their flow path by entraining sediment stored in the channel bed and banks, thus also increasing hazard potential. Theoretical considerations, laboratory experiments and field investigations all indicate that the saturation conditions of the sediment along the flow path can greatly influence the amount of sediment entrained. However, this process is usually not considered for practical applications. This study aims to close this gap by combining runout and hydrological models into a predictive framework that is calibrated and tested using unique observations of sediment erosion and debris-flow properties available at a Swiss debris-flow observation station (Illgraben). To this end, hourly water input to the erodible channel is predicted using a simple, process-based hydrological model, and the resulting water saturation level in the upper sediment layer of the channel is modelled based on a Hortonian infiltration concept. Debris-flow entrainment is then predicted using the RAMMS debris-flow runout model. We find a strong correlation between the modelled saturation level of the sediment on the flow path and the channel-bed erodibility for single-surge debris-flow events with distinct fronts, indicating that the modelled water content is a good predictor for erosion simulated in RAMMS. Debris-flow properties with more complex flow behaviour (e.g., multiple surges or roll waves) are not as well predicted using this procedure, indicating that more physically complete models are necessary. Finally, we demonstrate how this modelling framework can be used for climate change impact assessment and show that earlier snowmelt may shift the peak of the debris-flow season to earlier in the year. Our novel modelling framework provides a plausible approach to reproduce saturation-dependent entrainment and thus better constrain event volumes for current and future hazard assessment.

Illgraben basin↗
USGS2026

Characterizing changes in postfire debris-flow hazard as burned areas recover

Emergency assessments of postfire debris-flow hazards that are performed by the U.S. Geological Survey (USGS) provide estimates of debris-flow likelihood and rainfall triggering conditions that are used for evaluating and managing runoff-generated debris-flow hazards in recently burned areas throughout the western United States. Although the immediate postfire period, within roughly one year after fire, is typically the most susceptible to runoff-generated debris flows, the hazard evolves in time and space as the burned area recovers. The recovery trajectory a given burned area will take depends on local climate and weather and can be difficult to predict. Some burned areas recover quickly, whereas others experience debris flows for multiple years after fire. As a result, extending our ability to update debris-flow likelihood estimates and rainfall thresholds based on observed recovery of the burned area would be beneficial. We present a method for multi-year runoff-generated debris-flow hazard assessment that leverages the USGS “M1” debris-flow likelihood model and integrates updated, satellite-derived, normalized burn ratio data to estimate vegetation recovery. We predict recovery-aware rainfall thresholds and validate them against a multi-year debris-flow hazard prediction and could be adapted for use with other debris-flow models that incorporate burn severity data.

Arizona, California, Colorado, New Mexico, Washing↗
USGS2026

In-situ geochronology using LA-ICP-MS/MS: Application of the Lu-Hf system in carbonate, apatite and fluorite

Geochronology is a fundamental tool for reconstructing earth history and constraining the timing of ore deposit formation. Traditionally, U-rich minerals like zircon, titanite, monazite and apatite have been analysed to determine the timing of mineralization and alteration events. However, not all deposits or alteration assemblages contain these minerals, and the U-Pb system can be disturbed due to open system behaviour. We present a methodology for obtaining in-situ Lu-Hf dates from various minerals to demonstrate the potential for dating ore deposits that lack U-rich minerals. Carbonate, apatite and fluorite minerals from different tectonic settings and geological environments were analysed. We acquired Lu-Hf and REE (rare earth element) data from carbonate minerals from a carbonatite orebody and dike samples from the Mountain Pass REE carbonatite deposit in California, USA. We also tested a ‘double-dating’ method, where U-Th-Pb dates and Lu-Hf dates were obtained from the same analysis of apatite from rocks in the Michigan’s Great Lakes Tectonic Zone. Preliminary results for Lu-Hf dating of fluorite from carbonatite rocks in the Lemitar Mountains of New Mexico are also included. The analytical methodology shows promising results in each of these cases where there is sufficient parent Lu and U and measurable daughter Hf and Pb respectively.

Conference Paper↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗
USGS2025

A review of abrupt permafrost thaw: Definitions, usage, and a proposed conceptual framework

Purpose of Review We review how ‘abrupt thaw’ has been used in published studies, compare these definitions to abrupt processes in other Earth science disciplines, and provide a definitive framework for how abrupt thaw should be used in the context of permafrost science. Recent Findings We address several aspects of permafrost systems necessary for abrupt thaw to occur and propose a framework for classifying permafrost processes as abrupt thaw in the future. Based on a literature review and our collective expertise, we propose that abrupt thaw refers to thaw processes that lead to a substantial persistent environmental change within a few decades. Abrupt thaw typically occurs in ice-rich permafrost but may be initiated in ice-poor permafrost by external factors such as hydrologic change (i.e., increased streamflow, soil moisture fluctuations, altered groundwater recharge) or wildfire. Summary Permafrost thaw alters greenhouse gas emissions, soil and vegetation properties, and hydrologic flow, threatening infrastructure and the cultures and livelihoods of northern communities. The term ‘abrupt thaw’ has emerged in scientific discourse over the past two decades to differentiate processes that rapidly impact large depths of permafrost, such as thermokarst, from more gradual, top-down thaw processes that impact centimeters of near-surface permafrost over years to decades. However, there has been no formal definition for abrupt thaw and its use in the scientific literature has varied considerably. Our standardized definition of abrupt thaw offers a path forward to better understand drivers and patterns of abrupt thaw and its consequences for global greenhouse gas budgets, impacts to infrastructure and land-use, and Arctic policy- and decision-making.

Current Climate Change Reports↗