USGS Science⌕ Search

SEARCH · USGS Science

Results for “Nature Reviews Earth & Environment”

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.

635 records · Page 28Linked to original sources
USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗
USGS2023

WellSTIC: A cost-effective sensor for performing point dilution tests to measure groundwater velocity in shallow aquifers

Many individual measurement points are required to characterize groundwater velocity within an aquifer. Groundwater velocity is most commonly measured using a network of >5 cm diameter monitoring wells, which, if not already present at a site, are expensive and labor-intensive to install. Drive-point piezometers—simple, cost-effective wells that can be installed by hand—are a common tool for sampling groundwater in shallow, alluvial aquifers, but most groundwater velocity measurement techniques require equipment that is too large for these narrow (usually <2 cm inside diameter) piezometers. In this technical note, we introduce a low-cost sensor and well packer system (<$90 USD) for performing point dilution tests in narrow piezometers. Field data show that the magnitude of groundwater velocity measured with this technique agrees with velocities computed from natural gradient tracer tests. Additionally, with proper calibration, these sensors can be used to continuously monitor in-well specific conductance, either during inter-well tracer tests with saline tracers or for water quality monitoring. This system is a viable tool for rapid assessment of the magnitude of groundwater velocity in shallow aquifers.

Water Resources Research↗
USGS2025

Reconstruction of Holocene and Last Interglacial vegetation dynamics and wildfire activity in Southern Siberia

Wildfires are a rapidly increasing threat to boreal forests. While our understanding of the drivers behind wildfires and their environmental impact is growing, it is mostly limited to the observational period. Here we focus on the boreal forests of southern Siberia and exploit a U–Th-dated stalagmite from Botovskaya Cave, located in the upper Lena region of southern Siberia, to document wildfire activity and vegetation dynamics during parts of two warm periods: the Last Interglacial (LIG; specifically part of the Last Interglacial maximum between 124.1 and 118.8 ka) and the Holocene (10–0 ka). Our record is based on levoglucosan (Lev), a biomarker sensitive to biomass burning, and on lignin oxidation products (LOPs) that discriminate between open and closed forest and hard- or softwood vegetation. In addition, we used carbonate carbon stable isotope ratios ( δ 13 C ), which reflect a dominant control of the host rock, to evaluate soil respiration and local infiltration changes. Our LOP data suggest that, during the Last Interglacial, the region around Botovskaya Cave was characterised by open forest, which by ca. 121.5 ka underwent a transition from fire-resistant hardwood to fire-prone softwood. The Lev record indicates that fire activity was high and increased towards the end of Last Interglacial just before 119 ka. In contrast, the Holocene was characterised by a closed-forest environment with mixed hard- and softwood vegetation. Holocene fire activity varied but at a much lower level than during the Last Interglacial. We attribute the changes in wildfire activity during the intervals of interest to the interplay between vegetation and climate. The open forests of the Last Interglacial were more likely to ignite than their closed Holocene equivalents, and their flammability was aided by warmer and drier summers and a stronger seasonal temperature contrast due to the increase in seasonal insolation difference compared to the Holocene. Our comparison of the last two interglacial intervals suggests that, with increasing global temperatures, the boreal forest of southern Siberia may become progressively more vulnerable to higher wildfire activity.

Botovskaya Cave, Siberia↗
USGS2026

Multiple-well monitoring site adjacent to the Midway- Sunset and Buena Vista Oil Fields, Kern County, California

Groundwater quality in and around oil fields in the Southern San Joaquin Valley is of interest to many California residents that rely heavily on groundwater for domestic, commercial, and agricultural use. To help assess the effects of historical oil-field activities and natural geologic sources on groundwater near the southwest margins of the Kern County Groundwater Subbasin, a multiple-well monitoring site was installed near the administrative boundary between the Midway-Sunset and Buena Vista Oil Fields in Kern County, California. The installation of the Midway-Sunset Buena Vista multiple-well monitoring site (MSBV) supports regional analysis of the relations of oil and gas sources to groundwater quality by providing information about the geology, hydrology, geophysical properties, and water quality of the alluvial and upper Tulare aquifers in areas where groundwater data were limited. Data collected from the site included drill cuttings, whole core samples, sidewall core samples, mud-gas analysis, borehole geophysical logs, depth to water measurements, and water quality samples. Whole cores were scanned using dual energy computed tomography. Subsamples of selected cores were analyzed for density, porosity, specific retention, and bulk minerology. Thin sections of the subsamples were prepared, photographed, and examined. Two samples were analyzed using scanning electron microscope technology to examine the microporosity of diatomite laden sediment. Instrumentation installed in the wells collect hourly depth to water measurements. Analysis of the data show there is 355 feet of alluvium overlying the Tulare Formation at the well site. The contact between the two formations is an aquitard resulting in a perched aquifer in the alluvium and unconfined aquifer in the Tulare Formation. The alluvium is more heterogenous and finer grained than the Tulare Formation resulting in markedly higher porosity in the alluvium compared to the Tulare Formation. Higher specific retention observed in the alluvium is attributed to the finer grained sediment and greater abundance of reworked diatomite (as represented by opal-CT [cristobalite-tridymite]) compared to the Tulare Formation. Total dissolved solids (TDS) approached or exceeded 10,000 milligrams per liter (mg/L) in the alluvium from approximately 176 to 242 feet below land surface and at the top of the Amnicola clay at approximately 670 feet below land surface within the Tulare Formation. Elevated TDS, chloride, and boron concentrations in the alluvium and on top of the Amnicola clay likely reflect groundwater that is mixed with oil-field water. Water chemistry and modern-aged groundwater in the alluvial monitoring well (MSBV #3) are consistent with the oil-field water in the alluvium being derived from documented historical surface disposal of oil-field water upslope (northwest) of the site. Water chemistry and pre-modern groundwater age in the deeper Tulare monitoring well (MSBV #1) on top of the Amnicola clay are consistent with oil-field fluids derived from upslope natural geologic sources or old oil wells that leak in the subsurface. Shallow groundwater in the Tulare (MSBV #2) is not affected by mixing with oil-field sources.

EarthArXiv↗
USGS2024

Recent advances in characterizing the crustal stress field and future applications of stress data: Perspectives from North America

The stress field controls patterns of crustal deformation, including which faults are likeliest to cause earthquakes or transmit fluids. Since the 1950s, maps of maximum horizontal stress ( S Hmax ) orientations have advanced dramatically, and the style of faulting (relative principal stress magnitudes) has recently been mapped in some regions as well. This perspectives paper summarizes developments in characterizing stress orientations and (relative) magnitudes, including new seismic and borehole methods, as well as progress in identifying the causes of stress variations. Despite these advances, adding far more spatiotemporal detail would allow geoscientists to address many of today's key challenges regarding natural hazards, energy development, and geodynamics. In particular, it is critically important to characterize stress heterogeneity at multiple scales while also recognizing the coherent variability of the stress field. The second part of the paper considers how more detailed stress datasets could prove essential to addressing some of the grand questions in geoscience, including deciphering the poorly understood feedbacks between crustal dynamics and surface processes, improving earthquake and eruption forecasts, and determining the origins and shared properties of plate boundaries.

Geological Society, London, Special Publications↗
USGS1993

Geochemistry of and radioactivity in ground water of the Highland Rim and Central Basin aquifer systems, Hickman and Maury counties, Tennessee

A reconnaissance of the geochemistry of and radioactivity in ground water from the Highland Rim and Central Basin aquifer systems in Hickman and Maury Counties, Tennessee, was conducted in 1989. Water in both aquifer systems typically is of the calcium or calcium magnesium bicarbonate type, but concentrations of calcium, magnesium, sodium, potassium, chloride, and sulfate are greater in water of the Central Basin system; differences in the concentrations are statistically significant. Dissolution of calcite, magnesium-calcite, dolomite, and gypsum are the primary geochemical processes controlling ground-water chemistry in both aquifer systems. Saturation-state calculations using the computer code WATEQF indicated that ground water from the Central Basin system is more saturated with respect to calcite, dolomite, and gypsum than water from the Highland Rim system. Geochemical environments within each aquifer system are somewhat different with respect to dissolution of magnesium-bearing minerals. Water samples from the Highland Rim system had a fairly constant calcium to magnesium molar ratio, implying congruent dissolution of magnesium-bearing minerals, whereas water samples from the Central Basin system had highly variable ratios, implying either incongruent dissolution or heterogeneity in soluble constituents of the aquifer matrix. Concentrations of radionuclides in water were low and not greatly different between aquifer systems. Median gross alpha activities were 0.54 picocuries per liter in water from each system; median gross beta activities were 1.1 and 2.3 picocuries per liter in water from the Highland Rim and Central Basin systems, respectively. Radon-222 concentrations were 559 and 422 picocuries per liter, respectively. Concentrations of gross alpha and radium in all samples were substantially less than Tennessee’s maximum permissible levels for community water-supply systems. The data indicated no relations between concentrations of dissolved radionuclides (uranium, radium-226, radium-228, radon-222, gross alpha, and gross beta) and any key indicators of water chemistry, except in water from the Highland Rim system, in which radon-222 was moderately related to pH and weakly related to dissolved magnesium. The only relation among radiochemical constituents indicated by the data was between radium-226 and gross alpha activity; this relation was indicated for water from both aquifer systems.

Tennessee↗
USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗
USGS2024

Framework for implementing damping scaling factors in U.S. Geological Survey National Seismic Hazard Models

Traditionally, probabilistic seismic hazard analysis (PSHA) has focused on calculating ground motion hazard curves for elastic, 5%-damped pseudo spectral accelerations, Sa(T,5%), which are used as the basis for engineering design parameters and targets for ground motion selection and modification. However, structures and geotechnical systems can exhibit a wide range of damping ratios both above and below the 5% level, depending on the construction material, structural system, nonstructural elements, or subsurface soil properties. When spectral parameters at such damping levels are required for certain applications, 5%-damped accelerations have traditionally been extracted from PSHA-based hazard curves and adjusted outside of the hazard integral using damping scaling factors (DSF) such as those from Newmark & Hall (1982). Recent advances in the development of more rigorous and comprehensive damping scaling models (e.g., Rezaeian et al., 2014; Rezaeian et al., 2021) have allowed for the modeling of means and standard deviations of DSFs as functions of earthquake source and path properties for crustal, intraslab, and subduction interface tectonic environments. These DSF models can be applied to ground motion model (GMM) estimates of Sa(T,5%) for a given earthquake rupture scenario to produce a corresponding mean and standard deviation Sa at a specified damping ratio β, Sa(T,β). In this study, the DSF models of Rezaeian et al. (2014) and Rezaeian et al. (2021) are implemented within the U.S. Geological Survey National Seismic Hazard Model (NSHM) PSHA framework to calculate probabilistic hazard curves for spectral accelerations at damping ratios from 0.5% to 30%. The DSF models are applied directly to the mean and standard deviation of Sa(T,5%) predictions from each GMM in the NSHM logic tree. Resulting hazard curves and uniform hazard and risk spectra for Sa(T,β) are presented for several geographic locations and compared with corresponding spectra estimated using current design practices by applying the same DSFs outside of the PSHA calculation. Key differences between the two methods for estimating Sa(T,β) are discussed, and potential strategies are presented for the implementation and usage of the hazard-consistent Sa(T,β) in building codes. Comparing the results to those from DSFs used in current design practices that are mainly based on Newmark & Hall (1982) is not explored in this study.

Conference Paper↗
USGS2024

Peri-Gondwanan sediment in the Arkoma Basin derived from the north: The detrital zircon record of a uniquely concentrated non-Laurentian source signal in the late Paleozoic

During the assembly of Pangea, peri-Gondwanan terranes collided with the eastern and southern margins of Laurentia and brought with them unique detrital zircon U-Pb signatures. Discriminating between individual peri-Gondwanan terranes in the detrital record is difficult due to their similar geologic histories. However, characterization of this provenance is critical for understanding late Paleozoic sediment routing during development of Pangea. Along southeastern Laurentia, in the Arkoma Basin (present-day Arkansas and eastern Oklahoma, southeastern United States), we identified Middle Pennsylvanian (Desmoinesian) strata that exhibit a concentrated peri-Gondwanan detrital zircon signature (e.g., ca. 800–550 Ma). Although several southern peri-Gondwanan terranes (e.g., Maya, Suwannee) are closer to the Arkoma Basin, geologic data, such as predominantly north-to-south paleocurrents and proximal-to-distal facies relationships in these Desmoinesian strata, support a northern source (e.g., Ganderia, Avalonia, Meguma). Further evidence of a northern source comes from detrital zircon source mapping, which reveals the persistence of this peri-Gondwanan signal in depocenters to the north of the basin after the signal had diminished in the Arkoma Basin. To this end, bottom-up detrital zircon source modeling, source mapping, regional stratigraphy, paleocurrent data, and sandstone petrography allow us to reconstruct the evolution of this Middle Pennsylvanian (Desmoinesian) sediment pathway in the context of intraplate and plate-margin tectonic activity. This reconstruction documents processes affecting Earth’s surface (e.g., tectonics, climate) during the assembly of Pangea and describes in detail part of a dynamic continental-scale drainage system.

Geosphere↗
USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗
USGS2026

Preliminary geologic map of the Sparta East, Sparta West, and parts of the Glade Valley and Whitehead 7.5-minute quadrangles, North Carolina and Virginia, and the epicentral area of the August 9, 2020, Mw 5.1 earthquake near Sparta, North Carolina

Introduction New bedrock and surficial geologic mapping in the Sparta East, Sparta West, and parts of the Glade Valley and Whitehead 7.5-minute quadrangles, North Carolina and Virginia, investigates the geologic framework and causative mechanisms of the August 9, 2020, Mw 5.1 earthquake near Sparta, North Carolina. The mapping documents (1) the coseismic surface rupture from the 2020 earthquake and related brittle structures in the bedrock; (2) the fault contact between the western Blue Ridge and eastern Blue Ridge; (3) lithostratigraphy in the Lynchburg Group, Ashe Metamorphic Suite, and Alligator Back Metamorphic Suite; (4) the nature of the contact between the Lynchburg Group, Ashe Metamorphic Suite, and Alligator Back Metamorphic Suite; and (5) surficial deposits.

North Carolina, Virginia↗
USGS2024

Deep-ocean macrofaunal assemblages on ferromanganese and phosphorite-rich substrates in the Southern California Borderland

Mineral-rich hardgrounds, such as ferromanganese (FeMn) crusts and phosphorites, occur on seamounts and continental margins, gaining attention for their resource potential due to their enrichment in valuable metals in some regions. This study focuses on the Southern California Borderland (SCB), an area characterized by uneven and heterogeneous topography featuring FeMn crusts, phosphorites, basalt, and sedimentary rocks that occur at varying depths and are exposed to a range of oxygen concentrations. Due to its heterogeneity, this region serves as an optimal setting for investigating the relationship between mineral-rich hardgrounds and benthic fauna. This study characterizes the density, diversity, and community composition of macrofauna (>300 μm) on hardgrounds as a function of substrate type and environment (depth and oxygen ranges). Rocks and their macrofauna were sampled quantitatively using remotely operated vehicles (ROVs) during expeditions in 2020 and 2021 at depths above, within, and below the oxygen minimum zone (OMZ). A total of 3,555 macrofauna individuals were counted and 416 different morphospecies (excluding encrusting bryozoans and hydrozoans) were identified from 82 rocks at depths between 231 and 2,688 m. Average density for SCB macrofauna was 11.08 ± 0.87 ind. 200 cm −2 and mean Shannon-Wiener diversity per rock (H′ [loge] ) was 2.22 ± 0.07. A relationship was found between substrate type and macrofaunal communities. Phosphorite rocks had the highest H′ of the four substrates compared on a per-rock basis. However, when samples were pooled by substrate, FeMn crusts had the highest H′ and rarefaction diversity. Of all the environmental variables examined, water depth explained the largest variance in macrofaunal community composition. Macrofaunal density and diversity values were similar at sites within and outside the OMZ. This study is the first to analyze the macrofaunal communities of mineral-rich hardgrounds in the SCB, which support deep-ocean biodiversity by acting as specialized substrates for macrofaunal communities. Understanding the intricate relationships between macrofaunal assemblages and mineral-rich substrates may inform effects from environmental disruptions associated with deep-seabed mining or climate change. The findings contribute baseline information useful for effective conservation and management of the SCB and will support scientists in monitoring changes in these communities due to environmental disturbance or human impact in the future.

California↗
USGS2024

Seabed maps showing topography, ruggedness, backscatter intensity, sediment mobility, and the distribution of geologic substrates in quadrangle 2 of the Stellwagen Bank National Marine Sanctuary region offshore of Boston, Massachusetts

The U.S. Geological Survey, in cooperation with the National Marine Sanctuary Program of the National Oceanic and Atmospheric Administration, has conducted seabed mapping and related research in the Stellwagen Bank National Marine Sanctuary (SBNMS) region since 1993. The area being mapped using geophysical and geological data includes the SBNMS and the surrounding region, which totals approximately 3,700 square kilometers (km 2 ) and is subdivided into 18 quadrangles. The seabed is a glaciated terrain that is topographically and texturally diverse. Quadrangle 2, the subject of this scientific investigations map, has a mapped area of 209 km 2 and has water depths that range from about 19 meters (m) on the Stellwagen Bank crest to about 68 m in the Stellwagen Basin. Seven map types, each at a scale of 1:25,000, depict seabed topography, ruggedness, backscatter intensity, distribution of geologic substrates, sediment mobility, distribution of fine- and coarse-grained sand, and substrate mud content. These maps show the distribution of geologic substrates across the southwestern part of Stellwagen Bank, in Stellwagen Basin to the west and southwest of the bank, and in Little Stellwagen Basin and the western part of Race Point Channel to the south of the bank. Interpretations of multibeam sonar bathymetric and seabed backscatter imagery, photographs, video imagery, and grain-size analyses were used to create the geology-based maps. Data from 733 stations were analyzed, including 656 sediment samples. The geologic substrate maps of quadrangle 2 show the distribution of 19 geologic substrates that represent a wide range of textures, such as rippled and immobile sand, immobile muddy sand and sandy mud, sand that partially veneers gravel, and a boulder ridge. Mapped substrates are characterized by sediment grain-size composition, surface morphology, substrate layering, the mobility or immobility of substrate surfaces, and water depth range. This scientific investigations map portrays the major geological elements (substrates, topographic features, and processes) of environments in quadrangle 2. It is intended to provide a foundation for research into present and past sediment transport processes in a complex terrain, provide insights into the ecological requirements of invertebrate and vertebrate species that utilize the various substrates, and support seabed management in the region.

Massachusetts↗
USGS2026

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii↗
USGS2025

VIPER site analysis

We needed to evaluate available orbital data of NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission area in order to derive a variety of maps to help the science team identify scientifically interesting places for the rover to visit and to provide scientific context for our mission. Some of these maps also fulfilled engineering and mission design needs to enable safe and efficient landing and roving. We incorporated data from the Lunar Reconnaissance Orbiter Camera, the Lunar Orbital Laser Altimeter, the Mini-RF instrument, the Chandrayaan-2 Orbital High Resolution Camera, the Korean Pathfinder Lunar Orbiter’s Shadowcam, the Kaguya Spectral Profiler and Multiband Imager, and the Chandrayaan-1 Moon Mineralogy Mapper. We used a variety of techniques to build these maps, including stereogrammetry, shape-from-shading, ice stability depth and surface temperature calculations, and the horizon method for solar illumination and direct-to-Earth communications maps. Altogether, these maps allowed us to survey for boulders, evaluate features in permanently shadowed regions that VIPER might explore, provide mineralogic context for what VIPER’s instruments may learn, estimate the ages and radar properties of craters in the VIPER mission area, and evaluate the potential for gravity traverses with the rover. These data and techniques provided a rich set of information from which both the VIPER science team and engineering teams were able to draw in order to plan a safe landing and to plan a VIPER surface mission that will be both scientifically valuable and robust from an operational perspective.

The Planetary Science Journal↗
USGS2026

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗
USGS2017

Field-trip guide to Mount Hood, Oregon, highlighting eruptive history and hazards

This guidebook describes stops of interest for a geological field trip around Mount Hood volcano. It was developed for the 2017 International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Scientific Assembly in Portland, Oregon. The intent of this guidebook and accompanying contributions is to provide an overview of Mount Hood, including its chief geologic processes, magmatic system, eruptive history, local tectonics, and hazards, by visiting a variety of readily accessible localities. We also describe coeval, largely monogenetic, volcanoes in the region. Accompanying the field-trip guidebook are separately authored contributions that discuss in detail the Mount Hood magmatic system and its products and behavior (Kent and Koleszar, this volume); Mount Hood earthquakes and their relation to regional tectonics and the volcanic system (Thelen and Moran, this volume); and young surface faults cutting the broader Mount Hood area whose extent has come to light after acquisition of regional light detection and ranging coverage (Madin and others, this volume). The trip makes an approximately 175-mile (280-kilometer) clockwise loop around Mount Hood, starting and ending in Portland. The route heads east on Interstate 84 through the Columbia River Gorge National Scenic Area. The guidebook points out only a few conspicuous features of note in the gorge, but many other guides to the gorge are available. The route continues south on the Mount Hood National Scenic Byway on Oregon Route 35 following Hood River, and returns to Portland on U.S. Highway 26 following Sandy River. The route traverses rocks as old as the early Miocene Eagle Creek Formation and overlying Columbia River Basalt Group of middle Miocene age, but chiefly lava flows and clastic products of arc volcanism of late Miocene to Holocene age.

Oregon↗
USGS2025

Preventing overfitting when using tree-based methods for mapping hydrothermal favorability

Ensemble tree-based algorithms are robust tools for estimating sparsely distributed resources with non-linear dependencies (e.g., hydrothermal systems). These algorithms naturally accommodate the threshold conditions necessary to enable and support hydrothermal systems (e.g., having sufficient heat and permeability) and are simpler than many other non-linear machine learning strategies (e.g., artificial neural networks), which is an advantage when working with few labeled examples from which to learn. In previous work, we used eXtreme Gradient Boosting (XGBoost) to produce regional prediction and uncertainty maps of hydrothermal favorability; however, recent studies suggest that, even when properly applied, XGBoost has some risk of overfitting when there are few labeled examples from which to learn. To evaluate overfitting when constructing hydrothermal favorability maps with tree-based methods, we compare XGBoost with Extremely Randomized Trees (ExtraTrees), another ensemble tree-based algorithm that has the potential to underfit when using few labeled examples. We hold all other modeling parameters constant, resulting in two contrasting favorability maps of conventional geothermal resources for the Great Basin. Our results indicate that ExtraTrees demonstrably reduces overfitting compared with XGBoost. After considering overall performance, we conclude that ExtraTrees provides a more suitable modeling approach than XGBoost for the purposes of conventional hydrothermal resource assessments.

Conference Paper↗