Geology Reports⌕ Search

SEARCH · Geology Reports

Search Geology Reports

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,748 records · Page 22Linked to original sources
USGS2026

Prioritizing resource protection and understanding potential susceptibility of springs to surficial changes in a low-temperature geothermal system

Geothermal systems are vulnerable to changes in water budget and composition, requiring science-based management. This study uses a dataset of spring water temperatures, time series of groundwater residence time tracers (tritium and carbon-14), and stable isotopes of water to understand geothermal flow in a low-temperature geothermal system in north west Colorado, United States (Steamboat Springs). The geothermal system is bisected by the Yampa River, necessitating a stream mass balance approach to quantify total discharge. Time series analysis of water temperature data provides a ranked list of features more susceptible to surficial changes, which is corroborated using time series of tritium which indicate spatially distinct patterns of mixing between modern and pre-modern groundwater. All springs contain a portion of pre-modern groundwater that is thousands to tens of thousands of years old, a period coinciding with melting of extensive Pleistocene glaciers that was likely one of the recharge sources to the geothermal system. Stream mass balance indicates that greater than 80% of the total geothermal discharge is derived from diffuse or small springs, highlighting the extensive nature of the geothermal outflow zone and the association with local geologic structures. This study provides baseline data to support management of the Steamboat Springs geothermal system and indicates the utility of these approaches in developing science-based geothermal management.

Colorado↗
USGS2026

Surface variable‐based machine learning for scalable arsenic prediction in undersampled areas

In the United States, private wells are not federally regulated, and many households do not test for Arsenic (As). Chronic exposure is linked with multiple health outcomes, and risk can change sharply over short distances and with well depth. Coarse maps or sparse sampling often miss exceedances. Most existing models operate at ∼1 km resolution and use groundwater chemistry or detailed geologic logs, which limits their use in undersampled areas where improved guidance is most needed. We overcome these limitations by developing a machine learning model for Minnesota, USA, that predicts As exposure risk using only surficial variables from remote sensing and global data sets. Variables related to surface water hydrology and geomorphology are selected based on mechanistic links that control redox conditions and As mobilization. Local training was essential, and surficial geology variables that are more sensitive to local conditions were needed to maximize model accuracy. The resulting complete model was sufficiently sensitive to generate accurate and detailed risk maps and depth profiles of As concentrations above the 10 μg/L maximum contaminant level. Accuracy depended on local training data density. We identified a training data density of 0.07 wells/km 2 as a practical target for stable county-level performance. Maps of exceedance probabilities highlight priority areas for testing that are particularly important in rural communities that have received less sampling. These results support public health action by guiding where to install wells and where to test them, how much new sampling is needed, and where treatment outreach is most urgent.

Minnesota↗
USGS2026

Earthquake-hazard exposure of residents with potential access and functional needs in the United States

Earthquake response plans and earthquake early warning (EEW) systems designed for general populations may not consider potential access and functional needs (AFN) of individuals with physical, sensory, cognitive, or social limitations. Previous efforts to map the distribution of these populations have focused on social-vulnerability indices that ignore or oversimply these limitations. The descriptive and exploratory analysis summarized in this United States (U.S.) case study addresses this gap by identifying and integrating spatially explicit data for AFN-related residential populations, earthquake hazards, and county and county equivalents for the conterminous U.S., Alaska, Puerto Rico, and Hawaii. We focus on 13 AFN-related attributes that relate to an individual's ability to access information contained in an EEW alert, to understand and process earthquake information or observed ground shaking, and to take self-protective actions based on this information and physical cues of an earthquake. Depending on the demographic attribute, there are millions to tens of millions of U.S. residents with AFN-related attributes in areas considered to have varying likelihoods (2%, 10%, and 50%) of exceedance of a damaging earthquake in the next 50 years. Although these amounts represent low percentages at the national level, the percentage of individuals with AFN-related attributes in many counties and county equivalents substantially exceeds national percentages. No one county, county equivalent, U.S. state, or U.S. territory has the highest percentage of individuals in all AFN-related attributes; therefore, future efforts to increase individual resilience to earthquakes may benefit from understanding the local context of individuals with potential access and functional needs.

International Journal of Disaster Risk Reduction↗
USGS2026

Miocene evolution of the Humboldt Current

Diatom records from the East Pisco Basin (EPB) of southern Peru and offshore Ocean Drilling Program (ODP) Hole 682 A reveal stepwise increases in the primary productivity of the Humboldt Current during the middle and late Miocene. Although diatoms are present back through the late middle Eocene, successively enhanced diatom production occurs during the Miocene in four steps. The first step between 14.2 and 13.8 Ma marks the onset of diatom deposition in the Pisco-0 sequence. A second step at ∼12.9 Ma coincides with a major drop in global sea level during which diatom deposition ceased in the EPB but continued in ODP 682 A as well as in Chile, Ecuador, and Colombia. Beginning at 10.4 to 10.0 Ma, a major increase in diatom sedimentation rates in both the EPB and in ODP 682 A signals an abrupt intensification of the Humboldt Current productivity, which may be related to the closure of the Central American Seaway to deeper water circulation. A fourth step coincides with the transgressive deposition of the diatom-rich Pisco-2 sequence at 8.4 Ma and is enhanced by the onset of the global Late Miocene Biogenic Bloom (8.0–4.5 Ma). Between 7 and 6 Ma, common subtropical diatoms in ODP 682 A contrast with high diatom deposition rates in the EPB, suggesting variable El Niño-La Niña conditions. During this time, fossil vegetation evidence from southern Peru supports the presence of seasonal periods of enhanced rainfall typical of El Niño conditions. The Miocene trend of increasing sea mammal diversity in the EPB parallels these steps, further supporting stepwise enhancement of primary productivity in the Humboldt Current during the Miocene.

Humboldt Current↗
USGS2026

Mountain goat declines in a protected, interior, native population

A shifting climate poses threats to alpine-adapted species including mountain goats. We used long-term (12 years) citizen science monitoring data and Bayesian N-mixture modeling to estimate population trends and drivers of population metrics among mountain goats in Glacier National Park (GNP). Median goats per site ( n = 37 sites) declined by 45% (95% credible interval [CRI] = 32%, 57%) from 77.8 (95% CRI = 64.4, 95.1) in 2008 to 42.3 (95% CRI = 34.3, 52.2) in 2019, with consistent declines from 2008 until 2015, when the number of estimated goats stabilized. The decline exceeds IUCN criteria for classifying a population as vulnerable, >30% declines over only two generations. Across years, relatively few goats occupied northwestern GNP. Goat numbers declined the most at northeastern sites, trended toward decline in most southern sites, and increased at only two west-central sites. The proportion of permanent snow and glaciers, the presence of natural mineral licks, and habituation strongly increased the initial abundance of goats in the area. Weather variables had the greatest influence on population growth rates, particularly precipitation between May 15 and June 15 of the previous summer, the neonatal period. Lower growth occurred with less snow water equivalent and lower mean winter temperature, early summer temperature, and early summer precipitation. Projected reductions of permanent snow, increasing spring and summer temperatures, and insufficient and variable spring precipitation raise concerns for the future of native goats in this region. Our analyses reveal ways to improve detection rates of goats during surveys, which is important for optimizing the precision of estimates and the power to detect future trends. Detection increased with goat habituation, retention of observers with experience, use of binoculars, and conducting surveys at lower temperatures and earlier dates. Improving detection will be particularly important given the lower number of goats currently observed in the park. Research to estimate park-wide population size, evaluate genetic structure and diversity, assess changing habitat, human recreation levels and forage, and forward-project climate effects on persistence will be crucial to understanding the context of these results and conserving this iconic, metapopulation at the southern edge of the distribution of native mountain goats.

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

Microtextural characteristics of adularia in banded quartz veins from the Midas low-sulfidation epithermal deposit, Nevada

High-grade ores at the Miocene Midas low-sulfidation epithermal deposit in northern Nevada are confined to crustiform quartz veins containing abundant adularia. Micro-X-ray fluorescence elemental mapping reveals that adularia is a common gangue mineral occurring in colloform bands, bands showing bladed textures, and bands with dendritic terminations. The adularia aggregates have delicate shapes and are comprised of stacked, submillimeter crystals hosted by fine-grained quartz. The textural evidence suggests that the adularia aggregates originally formed within a gel-like, noncrystalline silica matrix, which subsequently transformed into quartz. This indicates that the adularia did not precipitate in open space along the vein walls. Correlative microscopy, involving scanning electron microscopy-based automated mineralogy and optical petrography, demonstrates that bands containing abundant adularia are not the primary host to ore minerals. The ore minerals occur in different bands within the crustiform veins, implying that adularia and ore mineral precipitation did not always occur simultaneously. It is hypothesized here that fluid flow at Midas involved intermittent short-lived events of fluid flashing, causing rapid solute supersaturation in the liquid. During each flashing event, different amounts of vapor were produced along a given vein. Compositional differences between adjacent bands in the crustiform quartz veins may, therefore, be linked to variations in the amount of vapor formed during each flash event.

Nevada↗
USGS2026

Bedrock geologic map of the Eagle Lake quadrangle, Essex County, New York

The bedrock geology of the 7.5-minute Eagle Lake quadrangle, Essex County, New York, consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) predated AMCG magmatism and is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the paragneiss host rocks. A sample of hornblende quartz syenite from the AMCG suite, collected from an abandoned railroad cut on Old Furnace Road, yielded a U-Pb zircon age of 1,149±10 million years before present. D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into the previously deformed lithologies. The LMG crystallized at approximately 1.15 to 1.14 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation of the Marcy anorthosite began during D3 along a mylonitic extensional detachment, as a type of core complex. Protracted D3 produced F3 folds exhibited in regional domes and basins, such as the Hammondville antiform, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. Lower Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. The Potsdam Sandstone preserves the Cambrian stratigraphy in outliers above the Great Unconformity. The Paleozoic rocks are weakly folded and block faulted. Parts of the quadrangle are covered by undifferentiated glacial deposits, but much of the quadrangle contains only a variably thick, veneer of unmapped glacial till over significant areas of exposed bedrock. The map also shows waste rock piles and locations of historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondack Mountains, and provide a modern context for historical mines. This Scientific Investigations Map of the Eagle Lake 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, and two cross sections. The explanatory pamphlet includes a discussion of the geology.

New York↗
USGS2026

Widespread terrestrial ecosystem disruption at the onset of the Paleocene–Eocene Thermal Maximum

The Paleocene–Eocene Thermal Maximum (PETM, ~56 Mya) interval was marked by massive 13 C-depleted carbon emissions into the ocean/atmosphere system, manifested as a negative carbon isotope excursion (CIE) in sedimentary components, and ~5 °C global average warming. Episodes of hydrological perturbations and soil-erosion have been widely documented for the PETM but their link with vegetation- and carbon cycle changes remain poorly constrained. Here, we present organic microfossil evidence showing a strong increase in fern-dominated pioneer vegetation that replaced coniferous forests on the margin of the Norwegian Sea during the first millennia of the CIE. With the present stratigraphic constraints, the “fern spike” occurred simultaneously in terrestrial settings along the North Sea, Arctic Ocean, the US east coast and in southern Australia, indicating that pioneer vegetation persisted for several millennia following a partial collapse of previously stable terrestrial ecosystems. Both the ferns and influx of microcharcoal imply recurrent physical disturbance, including soil destabilization and erosion, potentially linked to droughts, wildfires, and strong hydrological forcing resulting from extreme climate change. Together with evidence for reworked clay minerals and ancient organic matter (kerogen), these findings show that highly disturbed terrestrial ecosystems were widespread across mid- and high-latitude regions globally. Carbon cycle model simulations suggest that a substantial loss of standing and buried biomass, along with oxidation of soil organic matter, acted as important positive feedbacks during the onset of the CIE. Additionally, enhanced kerogen weathering likely contributed as another major positive feedback throughout both the onset and main phase of the CIE.

Proceedings of the National Academy of Sciences↗
USGS2026

Monitoring recreation on federally managed lands and waters—Aspects of visitor use

Federally managed public lands and waters receive about 1 billion recreational visits each year. Data on these visitors can aid in guiding policy decisions, managing resources effectively, and communicating the economic contributions of lands and waters. This report explores the methods used by agencies to collect data on aspects of recreational visitor use to Federal lands and waters (apart from visitation numbers, which are the focus of a companion publication). Aspects of recreational visitor use include visitor demographics, recreational activity participation, visitor satisfaction, visitor attitudes and experiences, trip characteristics, and economic contributions. We review practices used to understand aspects of visitor use across seven Federal agencies, revealing similarities such as the use of visitor intercept surveys and coverage of similar topic area, and differences in how survey programs are operationalized and how specific questions on visitor surveys are worded. We also evaluate emerging technologies, such as geolocated social media and mobile device location data, for their potential to aid in understanding aspects of visitor use. This report concludes with potential opportunities to enhance data collection and coordination, ensuring cost-effective data collection and informed decision-making.

SocArXiv↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗
USGS2026

Low-frequency earthquakes track the motion of a captured slab fragment

Accurate tectonic models are essential for assessing seismic hazard and fault interactions. However, the plate configuration at the complex Mendocino triple junction, where the San Andreas Fault and the Cascadia subduction zone meet, remains uncertain. We analyzed fault slip associated with a recently identified zone of tectonic tremor and low-frequency earthquakes (LFEs) near the southern edge of the subducting Gorda slab. Based on tidal sensitivity and P-wave first motions, we show that the LFEs are generated by dipping, strike-slip motion. This suggests that a former Farallon slab fragment, now captured by the Pacific plate, is translating northward beneath westernmost North America. This geometry effectively extends the slab interface fault, challenging prevailing interpretations of slab window formation and creating a potential unaccounted earthquake hazard in this region.

Science↗
USGS2026

A review and synthesis of post-wildfire shifts in hydrologic processes and streamflow generation mechanisms

Critical water supply watersheds in the western United States (WUS) are impacted by wildfires, with potential negative effects on water quality and quantity. Scientific understanding is currently insufficient to deliver estimates of wildfire consequences for water quantity that are regionally accurate. Regional variability in the directionality and magnitude of post-wildfire shifts in streamflow generation fuels uncertainty in estimates of wildfire effects on water supply. In this work we provide a narrative review of wildfire effects on hydrologic processes and the resulting changes in streamflow generation mechanisms with a focus on the WUS, incorporating other global regions when pertinent. A conceptual model summary of wildfire effects on streamflow generation emphasizes: (1) precipitation seasonality, (2) synchrony of precipitation and potential evapotranspiration, (3) net shifts in interception, evaporation, and transpiration relative to total annual precipitation, (4) vegetation changes, including compensatory uptake and type conversion, (5) degree of overlap in rainfall rates and infiltration, (6) fire extent and severity, (7) burn scar positioning (e.g. in headwaters or proximal to watershed outlet), (8) scale-dependent groundwater leakage, (9) near-surface water storage reduction, and (10) soil to groundwater connectivity. Ongoing gaps and challenges include separating the influences of precipitation variability, water withdrawals, and post-fire land management; compound and overlapping disturbances; and lack of pre-fire data. Notable future opportunities include: harnessing ever-improving gridded and remotely sensed precipitation and fire-effects data; linking geophysical, isotopic tracer, and geochemical signatures to diagnose hydrologic changes; leveraging physically based and data-driven model advancements; and analyzing streamflow generation recovery trajectories across diverse watersheds.

western United States↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2026

Groundwater tracing used to delineate recharge areas and map karst groundwater pathways for subterranean streams at Oregon Caves National Monument and Preserve

Oregon Caves National Monument and Preserve in southwestern Oregon is a 4,554-​acre area managed by the National Park Service that is home to several cave systems, including Oregon Caves, which is the longest cave in Oregon, with 3.03 miles of mapped passages. Because of the interconnected nature of karst hydrologic systems, it is critical to understand the areas that can influence water quality and quantity in karst environments. Toward this goal, dye tracing was conducted by the U.S. Geological Survey from 2021 to 2024 to better understand the pathways that karst groundwater follows at Oregon Caves National Monument and Preserve and to delineate recharge areas for two caves, Oregon Caves and Cave Next Door. During the project, eight dye injections were conducted, delineating a 0.51-​square-​mile recharge area for Oregon Caves and a 0.69-​square-​mile recharge area for Cave Next Door. Additionally, the study helped to identify three resurgences associated with Oregon Caves that were previously unknown and showed that the recharge areas for the two caves were distinct from one another. The dye traces also illuminated some unique recharge characteristics of the karst at Oregon Caves, including a high variance in karst groundwater velocities, retention within the karst aquifers, and a significant diffuse-​flow component.

Oregon↗
USGS2026

Widespread abyssal turbidites record megathrust earthquake-triggered landslides and coseismic deformation in the Cascadia subduction zone

Abyssal marine turbidites provide some of the longest and most spatially extensive records of subduction zone earthquake recurrence globally; however, correlation of these deposits over long distances and interpretation of synchronous emplacement requires both an understanding of the turbidite generating systems and precise dating. Here, we present an integrated suite of high-resolution bathymetry, subbottom profiles, and sediment cores from combined autonomous underwater vehicle, remotely operated vehicle, and ship-based studies at a key paleoseismic site in the southern Cascadia subduction zone. We demonstrate how widespread, earthquake-triggered landslides on the lower slope deposit discrete, proximal mass transport deposits (MTDs) that grade offshore into complex, interfingered abyssal turbidites, which correspond to records of megathrust earthquake history. We propose accretion and oversteepening of thrust folds on the lower slope both preconditions the slope to fail and provides a perpetual source of unstable material to fail during every earthquake cycle. Furthermore, we suggest the periodic and pervasive landsliding indicates coseismic deformation of the outer accretionary wedge during megathrust rupture.

California, Oregon↗
USGS2026

Distributed faulting of the northern West Napa Fault Zone in Napa Valley, California

Mapped surface ruptures from the 24 August 2014 M w 6.0 South Napa earthquake in the Napa Valley, California, show a 2‐km‐wide zone of distributed faulting in the southern and central West Napa fault zone (WNFZ). In the northern WNFZ at Hendry Winery (HW), however, the mapped 2014 surface ruptures encompass an ∼100‐m‐wide zone, implying significant narrowing of the near‐surface fault zone to the north. We present a tomographic shear‐wave velocity ( ⁠⁠ V S ) model and guided‐wave data that indicate the northern WNFZ is at least 400‐m wide, with multiple near‐surface fault traces. Our V S model shows that the 2014 surface ruptures are underlain by discrete low‐velocity zones (LVZs), and coincident guided‐wave data show that the LVZs carry fault‐zone guided waves. If nearby (<500 m) mapped faults to the east of HW are part of the WNFZ, the entire WNFZ is more than 1 km wide in the northern Napa Valley. WNFZ guided waves travel up to 38% slower than S body waves, and low‐strain guided‐wave shaking is up to five times stronger than the associated body‐wave shaking. Our data suggest that guided waves, traveling along distributed faults, may result in an increased shaking hazard over a 1‐km‐wide area of the northern Napa Valley during future significant earthquakes. In places, the 2014 surface ruptures were difficult to find one year after the earthquake, and paleoseismic trenching showed only weak evidence for faulting, which may not have been identified in trenches if the locations of the 2014 surface ruptures had not been previously mapped ( Prentice et al. , 2015 ). Guided‐wave and V S tomography data, however, show strong evidence for faulting beneath the 2014 surface ruptures and at locations to the east. Although paleoseismic trenching is the gold standard for identifying near‐surface faulting, methods such as peak ground velocities of guided waves may better identify immature near‐surface fault traces.

California↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗