SEARCH · Geology Reports

Results for “New Zealand Journal of Geology and Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,690 records · Page 13Linked to original sources

USGS2026

From landslide susceptibility to risk assessment in the conterminous U.S.

Understanding the spatial distribution of landslide prone-areas and what consequences they may have is important for risk management and land-use planning. In the United States, although landslides occur in every state, a comprehensive landslide risk assessment is still missing. Existing efforts, such as the Federal Emergency Management Agency (FEMA)’s National Risk Index, rely on aggregated products and coarse cartographic units, limiting their geomorphological and practical accuracy. In this study, we present a methodological advance for landslide risk assessment across large areas with incomplete and sparse data. We apply our procedures to the conterminous United States by integrating geomorphologically meaningful partitions and spatial and temporal probability data-driven models. Landslide susceptibility is estimated using a Generalized Additive Mixed Model incorporating a bias capture/correction scheme to account for inventory inaccuracies (reference Area Under the Curve = 0.75). The exceedance probabilities of landslide occurrence are defined for three temporal scenarios (2, 5, and 10 year). Then, we explore the associated potential economic consequences for human settlements and agricultural areas. The findings indicate that the spatial variability of risk is primarily controlled by exposure rather than by susceptibility/hazard alone. The mean risk increases by ∼170% from the 2-year to the 10-year scenario. Beyond its quantitative outcomes, this study offers a blueprint for continental or sub-continental scale landslide risk assessments, demonstrating both the opportunities and current limitations.

Engineering Geology↗

USGS2025

Timing of hydrothermal alteration and Au-Sb-W mineralization, Stibnite-Yellow Pine district, Idaho

The Stibnite-Yellow Pine district of central Idaho was mined from the early 1900s until the 1990s, extracting gold, antimony, tungsten, and mercury from veins and disseminated and replacement ores in mountainous terrain along the headwaters of the Salmon River. Mining during the two World Wars supplied critical antimony and tungsten to the war efforts. Recent exploration has delineated mineral resources of over 187 metric tons Au, 274 metric tons Ag, and 93,000 metric tons Sb. Mineralization is hosted in Cretaceous Idaho batholith granitic rocks and a sequence of Neoproterozoic to Paleozoic metasedimentary strata of carbonate and siliciclastic compositions. Historical studies outlined some of the complex paragenesis but debated the absolute age of mineralization. New petrographic and geochronologic work documents a sequence of five hydrothermal events in the Stibnite-Yellow Pine district. Event 1 is related to Cretaceous magmatic and hydrothermal activity and includes events ranging in age from 86 to 75 Ma, including sparse quartz-molybdenite veins dated at 86 Ma. Disseminated gold mineralization of event 2 is associated with sericitic alteration and sulfidation of igneous biotite and replacement of plagioclase by potassium feldspar, largely in granodiorite. Gold is present in zoned arsenian pyrite in both disseminated ores and in crosscutting carbonate-quartz veins containing pyrite and arsenopyrite. The large Yellow Pine deposit, localized at a dilatant bend in the Meadow Creek fault, hosts such disseminated and vein gold. Event 2 is interpreted as the major gold-forming event; 40 Ar/ 39 Ar ages of sericite and potassium feldspar alteration range, respectively, from 70 to 59 and 66 to 56 Ma. The long span is interpreted to reflect the age of gold mineralization and local overprinting by event 3. A narrower range from 66 to 61 Ma is interpreted to date the peak of gold mineralization and alteration. Event 3, tungsten mineralization with scheelite, is texturally later than event 2 gold and localized along the Meadow Creek structure. Event 3 scheelite has been dated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb methods at 57 Ma. Event 4, best developed in the West End area, includes gold-silver–bearing quartz-carbonate-pyrite veins and breccias with epithermal textures and potassium feldspar alteration envelopes. Adularia from event 4 yields 40 Ar/ 39 Ar plateau ages of 52 to 51 Ma. Event 5 antimony and mercury mineralization consists of stibnite veins and breccia cements at the Yellow Pine and Hangar Flat deposits as well as cinnabar veins and replacements at the peripheral Fern and Hermes deposits; it is constrained by an LA-ICP-MS U-Pb date on scheelite (ca. 47 Ma) intergrown with stibnite. Minor propylitic and argillic alteration is evident in 47 Ma igneous dikes, which do not contain economic mineralization. The Au-Sb-W ores in the Stibnite-Yellow Pine mining district formed over an extended time period from about 70 to 45 Ma in multiple pulses that were localized along the Meadow Creek fault zone. Each event corresponds to episodes of magmatism and/or hydrothermal activity in the region. Insignificant event 1 skarn and molybdenum mineralization is similar in age to the Thompson Creek porphyry molybdenum deposit in central Idaho. Event 2 gold mineralization occurred during a magmatic gap in central Idaho but was synchronous with magmatism in the Bitterroot lobe further north; event 2 is similar in age to orogenic gold-arsenic mineralization at the Beartrack mine in eastern Idaho. Event 3 scheelite mineralization coincides with tungsten mineralization at the Quartz Creek deposit, late magmatism in the Bitterroot lobe, and rapid exhumation of the Atlanta lobe of the Idaho batholith. Event 4 gold mineralization is coincident with the onset of regional Challis magmatism and extension. Event 5 antimony and mercury mineralization is time-equivalent to epithermal gold mineralization in the nearby Thunder Mountain volcanic field and the peak of Challis magmatism.

Idaho↗

USGS2021

Geology, mineralogy, and cassiterite geochronology of the Ayawilca Zn-Pb-Ag-In-Sn-Cu deposit, Pasco, Peru

The Ayawilca deposit in Pasco, Peru, represents the most significant recent base-metal discovery in the central Andes and one of the largest undeveloped In resources globally. As of 2018, it hosts an 11.7 Mt indicated resource grading 6.9% Zn, 0.16% Pb, 15 g/t Ag, and 84 g/t In, an additional 45.0 Mt inferred resource grading 5.6% Zn, 0.23% Pb, 17 g/t Ag, and 67 g/t In, and a separate Sn-Cu-Ag inferred resource of 14.5 Mt grading 0.63% Sn, 0.21% Cu, and 18 g/t Ag. Newly obtained U–Pb dates for cassiterite by LA-ICP-MS (22.77 ± 0.41 and 23.05 ± 2.06 Ma) assign the Ayawilca deposit to the Miocene polymetallic belt of central Peru. The polymetallic mineralization occurs as up to 70-m-thick mantos hosted by carbonate rocks of the Late Triassic to Early Jurassic Pucará Group, and subordinately, as steeply dipping veins hosted by rocks of the Pucará Group and overlying Cretaceous sandstones-siltstones of the Goyllarisquizga Group. Relicts of a distal retrograde magnesian skarn and cassiterite (stage pre-A) were identified in the deepest mantos. The volumetrically most important mineralization at Ayawilca comprises a low-sulfidation assemblage (stage A) with quartz, pyrrhotite, arsenopyrite, chalcopyrite, Fe-rich sphalerite, and traces of stannite and herzenbergite. Stage A sphalerite records progressive Fe depletion, from 33 to 10 mol% FeS, which is compatible with the observed transition from low- to a subsequent intermediate-sulfidation stage (B) marked by the crystallization of abundant pyrite and marcasite. Finally, during a later intermediate-sulfidation stage (C) sphalerite (up to 11 mol% FeS), galena, native bismuth, Cu-Pb-Ag sulfosalts, siderite, Mn-Fe carbonates, kaolinite, dickite, and sericite were deposited. This paragenetic evolution shows striking similarities with that at the Cerro de Pasco Cordilleran-type polymetallic deposit, even if at Ayawilca stage C did not reach high-sulfidation conditions. The occurrence of an early retrograde skarn assemblage suggests that the manto bodies at Ayawilca formed at the transition between distal skarn and skarn-free (Cordilleran-type) carbonate-replacement mineralization. Mineral assemblages define a T- f S 2 evolutionary path close to the pyrrhotite-pyrite boundary. Buffering of hydrothermal fluids by underlying Devonian carbonaceous phyllites of the Excelsior Group imposed highly reduced conditions during stage A mineralization (log f O 2 < − 30 atm). The low f O 2 favored efficient Sn mobility during stages pre-A and A, in contrast to other known ore deposits in the polymetallic belt of central Peru, in which the occurrence of Sn minerals is minor. Subsequent cooling, progressive sealing of vein walls, and decreasing buffering potential of the host rocks promoted the shift from low- (stage A) to intermediate-sulfidation (stages B and C) states. LA-ICP-MS analyses reveal significant In contents in Fe-rich sphalerite (up to 1.7 wt%), stannite (up to 1908 ppm), and chalcopyrite (up to 1185 ppm). The highest In content was found in stage A sphalerite that precipitated along with chalcopyrite and stannite, thus pointing to the early, low-sulfidation assemblage as prospective for this high-tech metal in similar mineral systems. Indium was likely incorporated into the sphalerite crystal lattice via Cu + + In 3+ ↔ 2 Zn 2+ and (Sn, Ge) 4+ + (Ga, In) 3+ + (Cu + Ag) + ↔ 4 Zn 2+ coupled substitutions. Indium incorporation mechanisms into the stannite and chalcopyrite crystal lattices remain unclear.

Pasco↗

USGS2025

Chemical and isotopic fractionation during melt inclusion formation

Melt inclusions are used to study the origin and evolution of magmas. The extent to which they represent equilibrium melt compositions, however, critically hinges on the ratio of crystal growth rate to diffusion rate in melts. If the rate of crystal growth is limited by the supply of nutrients and the dissipation of unwanted components, the trapped melt will be depleted in compatible elements and enriched in incompatible elements. Despite widespread recognition of the potential for melt inclusions to trap disequilibrium compositions, their identification remains a challenge. Here, we demonstrate how stable isotopes provide a solution to this problem. Melt inclusions in basaltic tephra from Kīlauea volcano, Hawai‘i, display Mg isotope fractionation up to 0.49‰ ( 26 Mg/ 24 Mg) relative to bulk rock. These observations indicate that kinetic processes, such as chemical diffusion, are at play, as equilibrium mineral-melt isotope fractionation at high temperatures is an order of magnitude smaller. We propose that the heavy Mg isotopic compositions of the melt inclusions were generated by the faster diffusion and thus preferential incorporation of 24 Mg by the growing olivine, leaving the slower diffusant ( 26 Mg) behind. The incompatible elements (e.g., Ca, Ti, and K) are predicted to display similarly large fractionations when boundary layer effects are significant. These findings show that stable isotopes can be a useful tool in identifying disequilibrium compositions in melt inclusions.

Geology↗

USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗

USGS2025

Simulation of the impacts of projected climate change on groundwater resources in the urban, semiarid Yucaipa Valley watershed, southern California using an integrated hydrologic model

Managing water resources in semiarid watersheds is challenging due to limited supply and uncertain future climate conditions. This paper examines the impact of future climate changes on an urban watershed in southern California using an integrated hydrologic model. GSFLOW modeling software is used to simulate the nonlinear relationships between climate trends and precipitation partitioning into ET, runoff, and subsurface storage. Four global circulation models (GCMs), each with two greenhouse-gas scenarios, RCP45 and RCP85 are used to project future climate conditions. GCMs include the CanESM2, CNRM-CM5, HadGEM2-ES, and MIROC5 models. The model's simulated hydrologic conditions are compared with historical data to assess changes in water budgets and groundwater supply. Results indicate decreased groundwater storage in most scenarios due to increased natural evapotranspiration, vegetation consumptive use, and streamflow out of the watershed. Only scenarios with substantially increased future precipitation show increased groundwater storage. The study also highlights increased future aridity despite the rise in precipitation and large precipitation events forecast by GCMs, which increase the risk of urban floods and decrease stream leakage and water available to vegetation.

California↗

USGS1999

Streamflow, base flow, and ground-water recharge in the Housatonic River basin, western Massachusetts and parts of eastern New York and northwestern Connecticut

Streamflows for selected flow durations from 1 to 99 percent and the August median streamflows were estimated for 11 long-term streamflow-gaging stations in and near the study area. Estimates of streamflow and associated standard errors were determined for selected flow durations from 50 to 99 percent and the August median streamflows for 21 low-flow partial-record stations and for selected flow durations from 1 to 99 percent and the August median streamflows for two partial-record stations and seven short-term discontinued streamflow-gaging stations. Median streamflows per square mile for the 10-, 50-, and 90-percent flow durations and the August median streamflows were 3.90, 1.01, 0.185, and 0.248 cubic feet per second per square mile. Streamflows per square mile at selected flow-duration discharges between 1 and 99 percent at the 41 stations were related to basin characteristics to explain differences in streamflow characteristics. Basin characteristics included basin elevations, extent of stratified-drift deposits, land use, aspect, and underlying bedrock geology types. Most streamflow differences were positively correlated to basin elevation differences, most likely because precipitation increases with elevation, and to stratified-drift deposits, which allow more precipitation to recharge the ground water and to discharge later than do till and bedrock deposits. Mean base flow was computed from continuous records of daily mean discharge at 11 long-term streamflow-gaging stations in and near the study area. Mean annual base flow ranged from 13.4 to 24.5 inches per year. Minimum annual base flow ranged from 45 to 72 percent of mean annual rates at the 11 long-term stations, and the ratio of base flow to streamflow (base-flow index) ranged from 0.55 to 0.80. Base-flow durations between 1 and 99 percent were calculated from streamflow records at the 11 long-term streamflow-gaging stations. Base flow accounted for 45.5 to 85.0 percent of total annual streamflow at the 1- and 99-percent flow durations. Ground-water-recharge rates were computed from continuous records of daily mean discharge at 11 long-term streamflow-gaging stations in and near the study area. Mean annual ground-water-recharge rates ranged from 17.5 to 22.4 inches per year at 10 of the 11 long-term stations. Mean annual ground-water-recharge rates ranged from 2 to 7 inches per year higher than base flow. Minimum annual ground-water-recharge rates ranged from 48 to 72 percent of mean annual ground-water-recharge rates. Mean annual potential ground-water recharge was estimated from monthly climatological data collected at six climatological stations in and near the study area. Mean potential ground-water recharge ranged from about 17.9 to 28.9 inches per year, with a median value of 22.6 inches per year. This median value compares well to that calculated by use of streamflow records at the 11 streamflow-gaging stations (20.0 inches per year). Streamflows per square mile for the 10-, 50-, and 90-percent flow durations at stations in and near the study area were similar to those computed for other unregulated long-term continuous streamflow-gaging stations in central and eastern Massachusetts. Base-flow and ground-water-recharge rates in the study area compared closely to results from other studies in southeastern Massachusetts and Rhode Island, which were based on the same computational methods.

Connecticut, Massachusetts, New York↗

USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗

USGS2025

Exploring the science and data foundation for Federal public lands decisions

Public lands provide diverse resources, values, and services worldwide. Laws and policies typically require consideration of science in public lands decisions, and resource managers are committed to science-informed decision-making. However, it can be challenging for managers to use, and document the use of, science and data in their decisions. To better understand science and data use in Federal public lands decisions in the United States, we assessed the number, type, and age of documents cited in 70 Environmental Assessments (EAs) completed by the Bureau of Land Management (BLM) in Colorado from 2015–2019. We focused on the BLM, as they manage the largest area of public lands in the United States. We selected Colorado as our study area, as actions proposed on BLM lands in Colorado are representative of those across the nation. Fifty percent of citations were categorized as science and 23% as data. EAs contained an average of 17 citations (range 0–111), with documents analyzing effects of oil and gas development and recreation actions including the highest and lowest mean number of citations (41 and 6, respectively). Of individual resource analysis sections within EAs, 24% contained ≥1 science citation and 21% contained ≥1 data citation. Journal articles were the most cited type of document (26% of citations) followed by non-BLM inventories (13%). Forty-seven percent of citations were relatively recent (2010 or later); the oldest citation was from 1927. Commonly analyzed resources with the highest mean number of citations were socioeconomics, mineral resources, and noise. Fourteen of 33 commonly analyzed resources included <1 citation on average. Actions and resources with no or few citations represent opportunities for strengthening the transparent use of science and data in public lands decision-making.

Colorado↗

USGS2026

Keeping pace: Migrating sea lampreys adjust swim speed in response to shifting water velocity to maintain a constant ground speed

Upstream migrating fish must balance energetic costs with environmental and ecological constraints. Theoretical models predict that fish can conserve energy by adjusting swim speed in response to changing water velocity, but detailed field tests of these predictions remain rare. We used fine-scale acoustic telemetry and hydrodynamic modeling to track adult sea lamprey (Petromyzon marinus) migrating through a river and evaluated how movement speeds varied with environmental conditions and individual traits. Swim speed increased with faster water velocities, while ground speed remained constant, supporting the hypothesis that lampreys regulate swim effort to maintain a fixed migration pace. The average ground speed aligned with a predicted optimum of ∼1 body length per second. Both swim and ground speeds declined slightly in deeper water, suggesting a trade-off between movement efficiency and predator avoidance. These findings provide field-based support for a theoretical model of energy optimization and underscore how migratory fish adaptively adjust behavior to local hydrodynamic conditions. Our results highlight the value of using high-resolution tracking and environmental modeling to reveal behavioral tactics that underlie efficient migration strategies.

Michigan↗

USGS2025

Geologic map of the Greater Antilles and Virgin Islands

Introduction The geologic map of the Greater Antilles and Virgin Islands is a compilation of information from the literature, integrated to provide a seamless geologic map of the region. This map was prepared to serve as a base map for a mineral resource assessment of the region. Several small-scale regional geologic maps of the region have been prepared in the past. This report supersedes an earlier version of the geologic map of the Greater Antilles and the Virgin Islands that was released as U.S. Geological Survey (USGS) Open-File Report 2019–1036. For this report, the regional geologic overview shown on sheet 1 covers Cuba; the island of Hispaniola, which includes Haiti and Dominican Republic; Jamaica; the Cayman Islands; Puerto Rico; and the U.S. and British Virgin Islands. Sheet 2 shows the geology of Cuba and Cayman Islands, sheet 3 the geology of Jamaica, sheet 4 the geology of Hispaniola, and sheet 5 the geology of Puerto Rico and Virgin Islands. Accompanying the maps are three tables: table 1 lists the map units (in alphabetical order by map-unit label), their ages, and the countries or territories in which they can be found. These are listed in alphabetical order for ease in locating unit names from the symbols on the maps. Table 2 is a list of sources for the geologic map compilation of Puerto Rico. Table 3 shows formations in the Devil’s Race Course Group of Jamaica.

Greater Antilles, Virgin Islands↗

USGS2026

Identifying overlap between native fish movements and a possible seasonal Grass Carp deterrent in the Sandusky River, Ohio, USA

Riverine fishes are vulnerable to habitat fragmentation caused by interrupted access to vital habitats. Fragmentation using artificial structures can be used to limit the spread of invasive fishes by inhibiting movement between critical habitats. Thus, tradeoffs exist between maintaining connectivity for native species and restricting movement of invasive fishes. A nonphysical deterrent has been proposed to prevent invasive Grass Carp ( Ctenopharyngodon idella ) from reaching spawning habitats in the Sandusky River, a tributary to Lake Erie, during late spring–summer. However, this timing overlaps with spawning periods of a suite of native fishes. We used acoustic telemetry and ichthyoplankton surveys to examine movement and spawning activity of multiple native species and to assess potential impacts of deterrent operation on reproduction. All species used the proposed deterrent location during the Grass Carp spawning period, and their movements corresponded with larval fish collections, indicating upstream spawning. These findings suggest the deterrent could unintentionally restrict native species’ access to essential habitats. Balancing habitat connectivity and invasive species control will likely depend on species-specific responses and careful timing of deterrent operation.

Ohio↗

USGS2026

Subsurface water ice mapping on Mars: A probabilistic approach

Subsurface water ice deposits on Mars are an important resource for potential future human exploration. They are also an indicator of the planet’s past climate. However, the distribution of subsurface water ice in Mars’s midlatitudes is uncertain because spacecraft imagery cannot directly observe subsurface ice in most cases. Various spacecraft remote sensing instruments are sensitive to subsurface water ice, including thermal imaging spectrometers, radar sounders, and neutron spectrometers. Geomorphic analyses of images can also implicate subsurface ice. Building upon the data products from the Mars Subsurface Water Ice Mapping project, we provide a probabilistic framework to jointly interpret existing data and estimate the likelihood of subsurface water ice in the Martian midlatitudes between 60 ∘ S and 60 ∘ N with uncertainty. Broadly, we find that near-surface ice is likely present poleward of ∼45 ∘ in both the northern and southern hemispheres. However, closer to the equator, existing remote sensing data cannot uniquely constrain the presence of subsurface water ice. Our probabilistic results provide a framework for quantifying the abundance of ice on Mars, and our uncertainty estimates allow future analysis and exploration to target regions of high uncertainty.

Planetary Science Journal↗

USGS2026

Opportunities for the U.S. Geological Survey’s National Seismic Hazard Model to improve seismic risk assessment of critical infrastructure.

As fragility and risk modeling techniques and computational capabilities evolve, complemented by moving toward more routine and systematic seismic risk assessment of all buildings and critical infrastructure, the authors pose a few critical questions to investigate how the U.S. Geological Survey (USGS) National Seismic Hazard Models (NSHMs) can be used and enhanced further to serve such issues. In this paper, we use three examples from multiple sectors to (1) identify the role of USGS NSHMs in evaluating seismic risks to critical infrastructure, (2) quantify potential impacts from NSHM enhancements (i.e., [i] hazard curves for the vertical component of ground motion, [ii] stochastic event sets, and [iii] maps of probabilistic ground failure hazards), and (3) clarify the feasibility of relevant NSHM improvements. We illuminate that NSHMs are commonly used in location-specific performance assessments, whereas earthquake effects on critical infrastructure can be widespread across large geospatial regions. Further, we found that without the NSHM extensions considered here, risk can be severely underestimated, e.g., neglecting ground failure hazards can underestimate regional loss by a factor of two or more. Although many challenges remain, we developed example prototypes to clarify the feasibility of the NSHM extensions, which can facilitate improved management of risks to critical infrastructure.

Earthquake Spectra Journal↗

USGS2026

Variability and consistency in wildfire susceptibility: Insights from a national compilation

Background Wildfire risk in the United States is rising and remains a land management priority. The quantitative wildfire risk assessment (QWRA) framework integrates fuels, topography, weather and values at risk to estimate the potential change in value from wildfire. Within this, response functions (RFs) represent how values respond to fire intensity. These are often based on expert judgment, but variation across assessments is unclear. Aims This study uses data from the US Geological Survey (USGS) Wildfire Hazard and Risk Assessment Clearinghouse to characterize consistency and variation across categories and contexts. Methods We applied descriptive statistics to summarize RFs, using tables, box-and-whisker plots and heat maps stratified by highly valued resource or asset (HVRA) category and spatial scale. Key results RFs and value definitions vary, especially for ecosystem-related resources. Some functions, such as for buildings in the wildland–urban interface (WUI), translate well across contexts, while others require more input. Conclusions Some functions are broadly transferable, while others need customization. This analysis provides references and starting points for improvement to RFs in QWRAs. Interpretations Expanding the clearinghouse and dataset and building more transparency in expert elicitation can build trust among communities, agencies and end-users, and can support efficient use of limited resources to mitigate wildfire risk.

International Journal of Wildland Fire↗

USGS2024

Taking heat (downstream): Simulating groundwater and thermal equilibrium controls on annual paired air–water temperature signal transport in headwater streams

Headwater stream temperature often exhibits spatial variation at the kilometer-scale, but the relative importance of the underlying hydrogeological processes and riverine perturbations remains poorly understood. In this study, we investigated the relative importance of groundwater (GW) and other processes on downstream annual stream temperature signal characteristics using deterministic heat budget model (HFLUX) scenarios within an idealized stream reach representative of mountainous forested conditions. We summarized annual stream thermal regimes from the relationship of paired sinusoidal air and water temperature signals (amplitude ratio, phase lag, and mean ratio). Results showed that downstream changes in annual temperature depended on the thermal gradient between water and the hypothetical equilibrium temperature (where all heat fluxes sum to zero). GW inflow, riparian shading, and the boundary input signal were the most significant factors affecting downstream annual water temperature signals, while flow volume and channel dimensions impacted how quickly annual temperature signals changed. Effects of GW were dominated by advective rather than conductive heat exchange processes, but conduction played a larger role when GW input was more spatially diffuse. Our results indicated several mechanisms by which local processes may affect stream thermal resilience to disturbances and can help guide management of wildfire and climate change.

Journal of Hydrology↗

USGS2026

Post-fire soil hydrologic response and recovery in northern California (USA)

Background Wildfires abruptly change landscapes by altering soil properties and vegetation cover. These changes are thought to reduce soil infiltration capacity, making landscapes susceptible to runoff and erosion. However, post-fire soil response is complex and likely varies across locations and time. Aims Here, we aim to understand regional post-fire soil response and recovery by tracking changes across different northern California (USA) lithology and vegetation types. Methods We conducted repeat in situ soil infiltration tests for 3 years post-fire at 31 burned and 10 unburned sites spanning the 2021 Dixie, 2020 LNU Lightning Complex, 2020 Walbridge and 2020 Glass fires. Key results Our two main findings are: (1) burned chaparral soils have increased hydraulic conductivity compared with unburned sites, and (2) infiltration rates return to pre-fire conditions within 3 years across most lithologies and vegetations. Conclusions Recovery might be generalizable by vegetation and lithology but differ regionally, making it important to identify meaningful hydrologic response units (HRUs). Multi-year studies with paired burned and unburned measurements can constrain the recovery timeline and provide information missed by observations solely of burned soils. Implications Understanding where, and for how long, soil remains susceptible to runoff and erosion can help prioritize areas and time periods most in need of mitigation.

California↗

USGS2026

A process-based model for forecasting wave runup along the coast of Georgia

Wave runup is an important nearshore process that impacts total water level, sediment transport, and coastal design. Current methods for forecasting wave runup implement an empirical model that considers offshore wave height, wave period, and generalized beach slope. In this study, the authors generated wave runup forecasts from offshore wave conditions and a system of polynomial equations derived from numerical simulations at three different still water datums for each beach profile. They developed a process-based methodology that incorporated site-specific cross-shore topobathy into the phase-resolving numerical model. A comparison between the system of equations, deterministic hydrodynamic simulations, and observed high-water marks was made using Hurricanes Matthew (2016) and Irma (2017) for 12 cases, and it showed that the polynomials were capable of being consistent with the results from full simulation runs, while not requiring hours of runtime when a forecast was needed—the differences between the polynomial and the observed high water marks ranged from 3 to 32 cm for the Irma hindcast and 9–70 cm for Matthew. Then, using forcings from Hurricanes Ian and Nicole (2022), the model predicted the occurrence of dune collision, overwash, and inundation for the coast of Georgia and suggested that wave runup was impacted by the still water level and local topobathy.

Georgia↗