USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal of Applied Earth Observations and Geoinformation”

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,248 records · Page 2Linked to original sources
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2025

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal↗
USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗
USGS2025

Relating systematic molecular and textural properties of graptolite pyrolyzed via gold tube hydrous pyrolysis: Implications for thermal proxies in lower Paleozoic marine shales

A series of gold tube pyrolysis experiments (72 h, 300–550 °C, 50 MPa) conducted on a graptolite-rich lower Paleozoic marine shale generated pyrolysis residues for a comprehensive evaluation of the molecular and structural variability of three types of graptolite periderm. Organic petrology, Raman spectroscopy, and field emission scanning electron microscopy (FE-SEM) with energy dispersive spectroscopy (EDS) were combined to evaluate the thermal evolution process. The three types of graptolite periderm, namely granular, non-granular, and nodular graptolite, were analyzed by Raman spectroscopy wherein point measurements were obtained after the maceral was identified and the location verified by organic petrology. Distinct thermal evolution pathways among non-granular, granular, and nodular graptolite periderms were recorded. The evolution patterns of the Raman parameters, particularly D1 and G bands, highlight the differences in geochemical composition of the graptolite periderm types and the alteration of molecular structure with increasing thermal maturity. Raman parameters D1 (position of the D1 peak), G-FWHM (full width at half maximum of the G peak), and ratios D1-FWHM/G-FWHM (full width at half maximum of the D1 peak ratioed to G-FWHM) and A D1 /A G (ratio of D1 and G peak intensities) showed effectiveness in assessing thermal maturity. Bireflectance with increasing gold tube pyrolysis temperature followed a hierarchy: non-granular > granular > nodular, reflecting different molecular alignment intensities. Qualitative FE-SEM evaluation showed that fine-grained mineral inclusions (primarily Fe-sulfide as determined via EDS) were associated with the graptolite populations, with granular graptolite containing greater amounts of coarser-grained (e.g., ∼300–1400 nm) mineral inclusions relative to non-granular and nodular graptolite, which contain finer-grained (e.g., ∼100–200 nm) inclusions difficult to resolve with optical microscopy. These findings are investigated to highlight the mechanisms that drive organic matter evolution within graptolite during thermal maturation, as well as to explore some of the limitations of using spectroscopic parameters as thermal maturity proxies.

International Journal of Coal Geology↗
USGS2025

Component identification of solid biomass fuels using reflected light microscopy: Interlaboratory study 2

As nations transition toward sustainable energy systems , biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

International Journal of Coal Geology↗
USGS2026

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin↗
USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗
USGS2025

International data gaps at the Center for Engineering Strong Motion Data

The Center for Engineering Strong Motion Data (CESMD) is utilized by seismologists, engineers, and disaster management professionals in the US and has historically achieved and distributed waveforms from across the globe for significant earthquakes. The increased access to the waveforms via Web API (Application Programming Interface) offers a unique opportunity to provide the community complete datasets, sampling a variety of tectonic environments and geologic conditions, increasing the number of available ground motion records for use in ground motion models (GMMs) and improving the accuracy of earthquake engineering evaluations. The objective of this study is to programmatically identify gaps in global event data from the past decade and backfill missing data gaps at CESMD. We first compare the CESMD catalog with the Advanced National Seismic System (ANSS) Comprehensive Earthquake Catalog identifying regions and time periods where strong-motion data is limited or inadequate. To backfill datasets at CESMD for significant events, we pinpoint regions and time intervals that lack information, creating a list of events for which we’d like to obtain data. An important facet of this work is identifying the source of data and metadata across earthquake repositories around the world and integrating these data repositories into our current strong-motion data processing workflow. In parallel with these newly processed datasets, we are developing a script to produce data origination citations to include provenance and attribution information to associate with respective datasets at CESMD. We showcase our methodology for identifying and filling data gaps at CESMD using three case studies (the 2018 Anchorage Alaska earthquake sequence, seismicity associated with the 2018 Hawaiian Kilauea volcano eruption, and several earthquakes in Turkey) and then outline our strategy to apply our data gap backfilling methods on an international scale.

Conference Paper↗
USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗
USGS2025

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth↗
USGS2024

Getting to the surface of the problem: A dynamic rupture benchmark for shallowly dipping faults near Earth’s surface

The Statewide California Earthquake Center (SCEC) Dynamic Rupture Group workshop “Getting to the Surface of the Problem: A Dynamic Rupture Benchmark for Shallowly-Dipping Faults Near Earth’s Surface” was convened virtually, in an online meeting on November 4, 2024. 52 people participated, including two who sent pre-recorded lightning talks earlier, but were unable to join us that day due to the international time-zone difference. One-half of our participants were students or postdocs, and the group included scientists from 11 countries (USA, Canada, China, Czech Republic, France, Germany, Iceland, Japan, New Zealand, Saudi Arabia, Switzerland).

Final Report↗
USGS2025

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2026

Hyperspectral narrowband imaging spectroscopy: A new paradigm for Earth observation

This editorial introduces the Special Issue entitled “Hyperspectral Narrowband Imaging Spectroscopy: A New Paradigm for Earth Observation” in the August 2026 issue of Photogrammetric Engineering & Remote Sensing (PE&RS), the flagship journal of the American Society for Photogrammetry and Remote Sensing (ASPRS). This volume represents the fourth dedicated hyperspectral special issue published in PE&RS, following earlier contributions by Thenkabail et al. (2025, 2024a, 2024b), and continues ASPRS’s commitment to advancing cutting‑edge imaging spectroscopy research and its applications across Earth system science. Remote sensing is undergoing one of the most profound transformations in its history. The emergence of hyperspectral narrowband (HNB) imaging spectroscopy data, capable of acquiring hundreds of contiguous, narrow spectral bands, has shifted the discipline from observing Earth in a handful of broad spectral windows to capturing continuous spectral signatures of the Planet. This transition marks a decisive break from the multispectral paradigm that has dominated satellite remote sensing for nearly five decades, driven by the advent of new orbital imaging spectrometers such as EnMAP, PRISMA, and NASA’s EMIT, and by the forthcoming Surface Biology and Geology (SBG) mission ((Pires Silva et al., 2026; Bourriz et al., 2025; Thenkabail et al., 2025; Chabrillat et al., 2024; Aneece et al., 2024; Dave et al., 2024; Thenkabail et al., 2024a; Thenkabail et al., 2024b; Thenkabail, 2024a; Thenkabail, 2024b; Thompson et al., 2022; Kokaly et al., 2022; Aneece & Thenkabail, 2022; Cawse Nicholson et al., 2021; Guanter et al., 2021; Vangi et al., 2021; Thenkabail et al., 2021). These missions (e.g., Table 1) deliver unprecedented spectral fidelity, improved signal to noise ratios, and global coverage capabilities, enabling a new era of quantitative, spectroscopy based Earth observation. Where multispectral broadbands (MBBs) provide only a few discrete measurements along the electromagnetic spectrum, HNB systems deliver rich, diagnostic information that enables scientists to characterize Earth’s surface with unprecedented biochemical, biophysical, and structural detail (Figure 1a, 1b). The implications for environmental monitoring, agriculture, water resources, and mineral exploration are profound. Several overarching themes emerge: • Spectral fidelity matters. The ability to preserve subtle absorption features is essential for mineral mapping, vegetation trait retrieval, and biochemical modeling. • AI and deep learning are indispensable. From destriping to classification, modern analytics must be scalable, label‑efficient, and capable of exploiting the full spectral–spatial richness of HNB data. • Physics‑based and data‑driven approaches must converge. Radiative transfer models such as PROSAIL, enhanced with localized soil parameterizations, remain foundational for biophysical retrievals and model‑based inference. • Dimensionality reduction and feature extraction are critical. Techniques such as L1‑ISOMAP demonstrate that intelligent manifold learning can unlock the structure of fused, high‑dimensional datasets. • Next‑generation architectures must be interactive and multimodal. ICTNet exemplifies the future of hyperspectral classification: hybrid, synergistic, and capable of modeling both local textures and global spectral dependencies.

Photogrammetric Engineering and Remote Sensing (PE↗
USGS2026

Remote compositional analyses of space-weathered lunar maria

Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.

Planetary Science Journal↗
USGS2026

The damability function: A probabilistic approach to regional landslide dam susceptibility analysis applied to the Oregon Coast Range, USA

Landslides can dam rivers and require rapid response to mitigate catastrophic outburst floods. Here, we present a workflow to map landslide dam formation susceptibility at a regional scale. We define a probabilistic function that combines river valley width and landslide volume to efficiently determine the likelihood of a landslide dam or “damability”. We combine damability values with landslide susceptibility to estimate landslide dam susceptibility. The valley width measurements are automated using a new elevation threshold-based algorithm. Landslide volume is represented as a statistical distribution from mapped landslides. We validate and apply our approach to the Oregon Coast Range, USA and find that 36 % of river stretches exceed a dam potential threshold; these are in river headwaters and steeper terrain, which in this case correlate with more resistant lithologies. We also estimate volumes of the potential dammed lakes and find that most rivers with high dam susceptibility are less likely to impound large lakes because they have low drainage areas. However, widespread susceptibility, and the potential impacts from exceptionally large landslides, suggest that this hazard should be considered in the Pacific Northwest. The damability function workflow can ingest new data and be applied more broadly to assess future landslide dam hazards.

Oregon↗
USGS2025

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗