SEARCH · Search USGS

Results for “Earth Surface Processes and Landforms”

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,665 records · Page 4Linked to original sources

USGS2024

Ambient flow and transport in long-screened, sand-packed wells: Insights into cross contamination and wellbore flow

The presence of long-screened wells with a surrounding sand pack can have a major effect on the redistribution of contaminants in groundwater, particularly when the wells are set in low-hydraulic conductivity aquifers. Such redistribution, or cross contamination, can occur through vertical flow and advective transport or by in-well mixing via multiple non-advective transport processes. A multi-method approach, including the use of single borehole dilution tracer (SBDT) logging, was undertaken to estimate vertical transport of trichloroethylene (TCE) in 8 discontinued remedial extraction wells, all constructed with long screens (100 ft, or 30.6 m long) and surrounding sand packs, at Site 25, Edwards Air Force Base, California. The site is within an enclosed drainage basin that is underlain primarily by quartz-monzonite-granitic rocks in various states of weathering. Prior to this study, little information was available on the depths of fracture zones intersecting the wellbores. Results indicate that because of in-well mixing processes, a potential redistribution of TCE of up to 9 g/d per well occurs as a consequence of leaving the wells inactive (unpumped) and unsealed, as measured by SBDT logging. Simulations of flow made with a generic model of the site show that if the wells were to be sealed with well liners, with the intent of reducing vertical TCE transport but the sand pack left intact, TCE transport decreases by 53% overall compared to leaving the wells unlined.

Environmental Earth Sciences↗

USGS2026

Distributed volcanic fields—Patterns, products, processes, and plumbing

Distributed volcanic fields—also referred to as “monogenetic fields” or simply “volcanic fields”—are found on every continent, on the seafloors, and in every tectonic setting on Earth. Yet, our understanding of what constitutes a distributed volcanic field, and our knowledge of the common products, spatiotemporal patterns in activity, and plumbing systems of these fields, are commonly poor and (or) debated. Here we provide an overview of the current understanding of distributed volcanic fields (DVFs) through review and compilation of published data, with an emphasis on understanding the characteristics of fields that are essential to prepare for future hazardous volcanic activity. Distributed fields can erupt restricted to variable magma compositions. Activity in these fields produces nearly all types of volcanoes, including scoria and spatter cones, maars, low shields, lava fields, stratovolcanoes, and calderas, but monogenetic volcanoes are the defining characteristic of DVFs. The number of vents and lifespans of DVFs vary greatly; vent numbers range from 2 to greater than 1,000, and fields can remain active for 150,000 to 10 million years. Some characteristics of DVFs appear field dependent. Although the volumes of single monogenetic eruptions within DVFs globally are variable (from less than 0.1 to a few cubic kilometers), each field may have a characteristic range in eruption volumes. Eruptions within distributed volcanic fields can have a random distribution and timing, but spatial and (or) temporal vent clustering and alignment are common, and eruptions may be time or volume predictable. Average eruption recurrence intervals in DVFs span over an order of magnitude, from less than 3 to greater than 100 thousand years. Exposures of the shallow portions of older systems reveal that many magma batches end their ascent before reaching the surface, with some dikes stopping within only meters of the surface. Petrologic insights into the deeper plumbing systems of DVF volcanoes indicate variable depths of magma crystallization and (or) stalling, and that magma transport and storage systems may be unique to a particular DVF and its underlying crustal structure. Both seismic and petrologic data for mafic DVF eruptions corroborate intrusions of magma into storage systems in the months to years before eruption, as well as variable ascent rates through the crust (days to weeks). Few distributed fields are closely monitored, as they erupt infrequently, but these observations imply that the amount of time available to install monitoring instrumentation after the first felt seismicity may be short, as precursory activity is frequently on the order of days to months for modern DVF eruptions. Satellite-based geodesy may provide useful warning in such cases. This chapter concludes with outstanding questions and areas for future research that will aid in our understanding and preparedness for future DVF eruptions.

Professional Paper↗

USGS2024

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface↗

USGS2025

The feasibility of using lidar-derived digital elevation models for gravity data reduction

Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Open-File Report↗

USGS2025

Relationship of atmospheric nitrogen deposition to soil nitrogen cycling along an elevation gradient in the Colorado Front Range

Microbial processing of atmospheric nitrogen (N) deposition regulates the retention and mobilization of N in soils, with important implications for water quality. Understanding the links between N deposition, microbial communities, N transformations, and water quality is critical as N deposition shifts toward reduced N and remains persistently high in many regions. Here, we investigated these connections along an elevation transect in the Colorado Front Range. Although rates of N deposition and pools of extractable N increased down the elevation transect, soil microbial communities and N transformation rates did not follow clear elevational patterns. The subalpine microbial community was distinct, corresponding to a high C:N ratio and low pH, while the microbial communities at the lower elevation sites were all very similar. Net nitrification, mineralization, and nitrification potential rates were highest at the Plains (1,700 m) and Montane (2,527 m) sites, suggesting that these ecosystems mobilize N. In contrast, the net immobilization of N observed at the Foothills (1,978 m) and Subalpine (3,015 m) sites suggests that these ecosystems retain N deposition. The contrast in N transformation rates between the plains and foothills, both of which receive elevated N deposition, may be due to spatial heterogeneity not captured in this study and warrants further investigation. Stream N concentrations from the subalpine to the foothills were consistently low, indicating that these soils are currently able to process and retain N deposition, but this may be disrupted if drought, wildfire, or land-use change alter the ability of the soils to retain N.

Colorado↗

USGS2024

The addition of 144Nd atomic mass to routine ICP-MS analysis as a Quick Screening Tool for Approximating Rare Earth Elements (Q-STAR) in natural waters

Rare earth elements (REEs) are a class of critical minerals, all of which can have supply chain vulnerability that impacts economic security. These elements are widely measured in environmental matrices via inductively coupled plasma mass spectrometry (ICP-MS); however, successful quantification can require time-consuming, sample-specific optimization. While a sample-by-sample approach is appropriate for targeted quantification studies, this approach is not suitable for mineral exploration efforts where rapidly screening thousands of samples for the presence of REEs is desired. Here, we demonstrated the use of a Quick Screening Tool for Approximating REEs (Q-STAR) to detect REEs in surface water and groundwater matrices, collected as part of existing environmental studies. A mass-to-charge ratio of 144 ( m / z = 144) was added to an ICP-MS method to screen for REEs in filtered water samples submitted for metals analyses to the U.S. Geological Survey (USGS) National Water Quality Laboratory. We detected the presence of REEs above a reference threshold of 1200 counts per second in 18 % of pre-selected 6626 samples. Using this screened dataset, we mapped estimated dissolved REE concentrations across the United States in relation to ecoregions and underlying geology. Data are constrained to where sample collection took place but nevertheless show estimated aqueous dissolved REE concentrations on a geographic scale that has not yet been studied. To validate Q-STAR, REEs were measured in a USGS standard reference sample, a subset of 88 archived filtered water samples, and in fresh filtered surface water samples. Our targeted analyses demonstrated a strong linear relationship between Q-STAR predicted and measured values in all archived samples for Nd (r 2 = 0.94), and light REEs (LREEs) such as lanthanum (La) (r 2 = 0.93), praseodymium (Pr) (r 2 = 0.94) and samarium (Sm) (r 2 = 0.94). Using Q-STAR screen values, nine field sites were identified and surface water samples recollected to confirm the continued presence of Nd and LREEs. Q-STAR can be used to screen an unlimited number of water samples for the presence of REEs prior to time-intensive and costly quantitative analyses and to generate large REE datasets for further investigation.

Journal of Geochemical Exploration↗

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↗

USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗

USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗

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↗

USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗

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↗

USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗

USGS2025

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska↗

USGS2025

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

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

Washington↗

USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗

USGS2024

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science↗

USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗