Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geochemistry Geophysics Geosystems”

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.

174 records · Page 10Linked to original sources
USGS2025

The systematics of stable hydrogen (δ2H) and oxygen (δ18O) isotopes and tritium (3H) in the hydrothermal system of the Yellowstone Plateau volcanic field, USA

To improve our understanding of hydrothermal activity on the Yellowstone Plateau volcanic field, we collected and analyzed a large data set of δ 2 H, δ 18 O, and the 3 H concentrations of circum-neutral and alkaline waters. We find that (a) hot springs are fed by recharge throughout the volcanic plateau, likely focused through fractured, permeable tuff units. Previous work had stressed the need for light δ 2 H water recharge restricted to the northern part of the plateau or recharge during past cold periods. However, new data from the Y-7 drill hole suggests that recharge is not restricted to a certain area or a cold period. (b) δ 18 O values of thermal waters in the geyser basins are shifted from the global meteoric water line by temperature-dependent water-rock reactions with higher subsurface temperatures resulting in a greater shift. (c) Large temporal variations in the isotopic composition of meteoric water recharge and small temporal variability in the isotopic composition of hot spring discharge implies that the volume of groundwater in, and around the Yellowstone caldera is substantially larger than the volume of annual water recharge. (d) Hot springs discharged through different rhyolitic units correlate with identifiable differences in δ 2 H and δ 18 O compositions, 3 H concentrations, and water chemistry that imply equilibration at different temperatures and travel along different flow paths. (e) Based on measured 3 H concentrations, we calculate that hot spring waters in the central part of the geyser basins mostly contain <2% post-1950 meteoric water, whereas waters discharged at the basin margins contain larger fractions of post-1950s meteoric water.

Wyoming↗
USGS2025

The transition from melt accumulation to eruption initiation recorded by orthopyroxene Fe-Mg diffusion timescales in late Holocene rhyolites, South Sister volcano, Oregon Cascade Range

South Sister volcano, Oregon Cascade Range, USA, has repeatedly erupted rhyolite since ca. 40 ka. The youngest such eruptions are the ca. 2 ka Rock Mesa and Devils Chain rhyolites, erupted several hundred years apart from two multi-vent complexes separated by 3–6 km. Fe-Mg interdiffusion models of orthopyroxene rims from both rhyolites produce timescales up to several-thousand years, but dominantly decades-to-centuries. Notably, the timescales of step-normal zoned orthopyroxene rims (i.e., normally zoned with a steep chemical gradient) from the Rock Mesa rhyolite are longer than those of reversely zoned crystals, whereas the Devils Chain produced mostly decadal timescales for both zoning types. Despite the proximity and broadly similar products of these episodes, their respective timescales indicate distinct sequences of events leading up to each eruption. The Rock Mesa timescales record centuries of magma chamber growth followed by decades of predominantly magma rejuvenation, reorganization, and destabilization. In contrast, the Devils Chain episode was preceded by a single episode of coupled rhyolite extraction, rejuvenation, and hybridization. Rare, high-An plagioclase cores and evidence of reheating implicate cryptic emplacement of mafic magma at the base of the rhyolite reservoirs. However, the diffusion timescales do not unequivocally support a single magma recharge event that affected both. Fluid fluxing and the reorganization of melt into buoyant magma chambers likely provided the source of increasing pressurization that initiated each eruption after several decades. Geodetic models of ongoing deformation west of South Sister could consider these processes in addition to magma emplacement.

Oregon↗
USGS2026

Geochemical disequilibrium at the brittle-ductile transition

We investigate the microtextural, microchemical, and isotopic effects of late-stage ductile deformation in quartzite mylonites and kyanite–muscovite–quartz veins from the Raft River shear zone (Utah). Quartz microstructures record pervasive disequilibrium, expressed by unannealed features including undulatory extinction, deformation lamellae, and poorly defined fabrics, typical of waning deformation in shear zones. Microchemical disequilibrium is best preserved in kyanite quartzites, where CL-zoned kyanite records repeated fractures, overgrowth, and mineral precipitation, and in muscovite from muscovite-poor quartzite mylonites that shows minor-element zoning consistent with syn-deformational overgrowth on detrital cores. In contrast, muscovite from muscovite-rich kyanite quartzites exhibits minimal chemical zoning. These microchemical variations correlate with 40 Ar/ 39 Ar age systematics. Chemically zoned muscovite preserves variable single-step ages, including ∼150 Ma ages, reflecting retention of detrital cores. In contrast, syndeformational muscovite consistently yields Miocene ages, indicating recrystallization and new growth below argon closure temperatures that reset inherited isotopic signatures. Similar trends are observed in quartzite mylonites, where increasing quartz recrystallization and stronger crystallographic preferred orientations occur toward deeper structural levels. Together, these observations indicate increasing retrograde deformation and recrystallization with depth in the Raft River shear zone and demonstrate that strain-driven recrystallization exerts a first-order control on muscovite ages. We suggest that apparent thermochronologic gradients in retrograde shear zones may reflect recrystallization gradients rather than temperature gradients. Where cooling limits the thermal driving force for recrystallization, isotopic relics are preserved, and local deformation and fluid availability control re-equilibration. Consequently, isotopic disequilibrium—particularly in the 40 Ar/ 39 Ar system—may be the rule rather than the exception of retrograde tectonic environments.

Idaho, Utah↗
USGS2026

Springtime formation of laminated soil carbonate rinds and changes in fluvial terrace soils on orbital timescales at Rio Mesa, Utah, USA

Laminated soil carbonate rinds are a Quaternary paleoclimate archive whose isotope composition is linked to soil formation conditions. At Rio Mesa, Utah (USA), we investigated the fidelity of rind records in a river terrace setting by determining the seasonal timing of rind formation and testing for inter-record replication. We infer soil carbonate formed in the spring season, contrasting with our prior inference of summer formation at Teasdale, Utah, ≈200 km distant. This apparent discrepancy occurs because of differences in the timing of the largest annual infiltration (spring vs. summer). At Rio Mesa, modern soil data show that soil carbonate δ 13 C would have high values (−2 to 2‰ VPDB) regardless of seasonal activity of C 3 versus C 4 plants because respiration rate is a strong control. We accordingly suggest reassessment of published records interpreting soil carbonate δ 13 C only via C 3 versus C 4 plant abundance. Three rind δ 13 C and δ 18 O records generally replicated. Intriguingly, rind δ 13 C may inversely correlate with summer insolation, evidence for global-scale influence on soils. Rind δ 18 O is not as clearly correlated with published western USA paleoclimate records, potentially due to regional differences in climate and because rinds record soil-specific processes. Our results support the fidelity of the soil carbonate rind paleoarchive and suggest that because rind formation seasonality is intimately tied to infiltration seasonality, spatial transects of rind records might be used to delineate boundaries between areas dominated by spring and summer infiltration, permitting reconstruction of the geographic extent of large-scale hydrologic phenomena such as the North American Monsoon.

Utah↗
USGS2026

Regional variations in sea ice and primary productivity in the Bering Sea during Marine Isotope Stage 11

Marine Isotope Stage (MIS) 11 (424-374 ka) has long been an analog for Holocene climate, because it is the most recent interglacial period with similar orbital conditions. However, there is significant global and regional variability in the climate response to MIS 11 warmth. Here, we review sediment core records from across the Bering Sea to investigate changes in paleoceanographic conditions during Marine Isotope Stages 12-10. Sea ice was present over much of the Bering Sea during MIS 11, but today, none of the sites investigated are ever ice-covered. This suggests that sea ice regimes in the Bering Sea during MIS 11 were different to those of the Holocene. There are also regional differences in the response of sea ice to MIS 11 warming. At the Umnak Plateau, Southeastern Bering Sea, sea ice concentrations decline during deglaciation, but they remain high at the slope sites until Peak MIS 11. Sea ice re-advances over the Umnak Plateau during peak interglacial warmth, at the same time that it declines over the slope sites. Late MIS 11 is characterized by high concentrations of seasonal sea ice at the Umnak Plateau, whilst sea ice at the slope sites fluctuates between consolidated and unconsolidated ice cover. This east-west dichotomy may be explained by changes in the behavior of the Aleutian Low. Productivity increases dramatically during deglaciation due to increased upwelling and sea level rise bringing fresh nutrients into the oceans. This is characterized by increased diatom productivity, increased terrestrial carbon deposition, and laminations at all sites.

Alaska↗
USGS2026

Hydrothermal manganese-oxide mineralization of a carbonate ooze, Samoan hotspot region, South Pacific Ocean

Low-temperature hydrothermal manganese oxides occur throughout the global oceans. However, the hydrothermal replacement of a carbonate ooze by manganese oxides is described here for the first time. The 24 samples dredged from three locations around the Territory of American Samoa in the South Pacific Ocean include six samples with remnant carbonate and volcaniclastic sediments, which we refer to as “low Mn,” and 18 samples in which the mineralization is pervasive and has replaced most original sediment, termed “high Mn.” The 18 high-Mn samples exhibit a mean Mn content of 51 wt.%. Higher Li and Mn contents and lower Fe contents in the high-Mn samples indicate a hydrothermal origin and distinguish these samples from hydrogenetic ferromanganese crusts. Mn-oxide layers are up to 90 mm thick, some with columns to 44 mm long and 10 mm wide, the magnitude of which has not been described previously. Samples are composed of birnessite and 10 Å phyllomanganate minerals. Textures of the thickest Mn layers indicate mineralization below the seabed from ascending fluids during multiple hydrothermal pulses. Mineralization took place by complete to partial replacement and cementation of foraminiferal sediments intermixed with volcaniclastic sediments in varying amounts. Our results highlight the production of carbonate sediment-hosted hydrothermal Mn oxides from multiple hydrothermal sources within the Samoan volcanic chain. The potential extensive distribution on the regional scale of this newly described mineralization process and unique element enrichments raise questions about its broader distribution globally and potential importance to hydrothermal processes, element mass balance, and seabed mineral resources.

Geochemistry, Geophysics, Geosystems↗
USGS2026

Constraining sources of mid-Pleistocene to present explosive volcanism in the Gulf of Alaska using machine learning and compositional data analysis

The long-term eruptive record of a region helps better elucidate magmatic processes at depth as well as volcanic hazards at the surface. Typically, reconstructing such records is done using proximal tephrostratigraphy and linking individual tephras to source volcanoes. These records, however, can only be accurately constructed if they are both well preserved and correctly linked to source volcanoes—not a trivial task, especially in areas such as Alaska that have experienced numerous glaciation events since the Pleistocene. This ultimately necessitates another way of assessing the long-term volcanic record such that these histories may be better discerned. Here we present data from 70 marine core tephras from the Gulf of Alaska, which have a virtually uninterrupted depositional record going back through the mid-Pleistocene. We utilize compositional data analysis techniques to quantify 37 eruptions over the span of eight cores, machine learning classification and conformal prediction algorithms to assign the most probable volcanic source(s) to each eruption, and multivariate distance-based metrics when machine learning classification algorithms are inappropriate. We find that the Mount Katmai magmatic system is the most probable volcanic source for analyzed tephras and that large volcanic centers such as Mount Katmai, Fisher Caldera, and Emmons Lake volcanic center have produced nearly invariant incompatible trace element ratio magmas, allowing for them to be confidently identified in long-term tephra records. This highlights the utility of trace elements for tephra studies, especially when paired with compositional data analysis, multivariate statistical tests, and petrologically informed discriminants.

Alaska↗
USGS2005

Structure and mechanics of the San Andreas–San Gregorio fault junction, San Francisco, California

[1] The right‐lateral San Gregorio and San Andreas faults meet west of the Golden Gate near San Francisco. Coincident seismic reflection and refraction profiling across the San Gregorio and San Andreas faults south of their junction shows the crust between them to have formed shallow extensional basins that are dissected by parallel strike‐slip faults. We employ a regional finite element model to investigate the long‐term consequences of the fault geometry. Over the course of 2–3 m.y. of slip on the San Andreas‐San Gregorio fault system, elongated extensional basins are predicted to form between the two faults. An additional consequence of the fault geometry is that the San Andreas fault is expected to have migrated eastward relative to the San Gregorio fault. We thus propose a model of eastward stepping right‐lateral fault formation to explain the observed multiple fault strands and depositional basins. The current manifestation of this process might be the observed transfer of slip from the San Andreas fault east to the Golden Gate fault.

California↗
USGS2006

Temperature variations at diffuse and focused flow hydrothermal vent sites along the northern East Pacific Rise

In the decade following documented volcanic activity on the East Pacific Rise near 9°50′N, we monitored hydrothermal vent fluid temperature variations in conjunction with approximately yearly vent fluid sampling to better understand the processes and physical conditions that govern the evolution of seafloor hydrothermal systems. The temperature of both diffuse flow (low-temperature) and focused flow (high-temperature) vent fluids decreased significantly within several years of eruptions in 1991 and 1992. After mid-1994, focused flow vents generally exhibited periods of relatively stable, slowly varying temperatures, with occasional high- and low-temperature excursions lasting days to weeks. One such positive temperature excursion was associated with a crustal cracking event. Another with both positive and negative excursions demonstrated a subsurface connection between adjacent focused flow and diffuse flow vents. Diffuse flow vents exhibit much greater temperature variability than adjacent higher-temperature vents. On timescales of a week or less, temperatures at a given position within a diffuse flow field often varied by 5°–10°C, synchronous with near-bottom currents dominated by tidal and inertial forcing. On timescales of a week and longer, diffuse flow temperatures varied slowly and incoherently among different vent fields. At diffuse flow vent sites, the conceptual model of a thermal boundary layer immediately above the seafloor explains many of the temporal and spatial temperature variations observed within a single vent field. The thermal boundary layer is a lens of warm water injected from beneath the seafloor that is mixed and distended by lateral near-bottom currents. The volume of the boundary layer is delineated by the position of mature communities of sessile (e.g., tubeworms) and relatively slow-moving organisms (e.g., mussels). Vertical flow rates of hydrothermal fluids exiting the seafloor at diffuse vents are less than lateral flow rates of near-bottom currents (5–10 cm/s). The presence of a subsurface, shallow reservoir of warm hydrothermal fluids can explain differing temperature behaviors of adjacent diffuse flow and focused flow vents at 9°50′N. Different average temperatures and daily temperature ranges are explained by variable amounts of mixing of hydrothermal fluids with ambient seawater through subsurface conduits that have varying lateral permeability.

Geochemistry, Geophysics, and Geosystems↗
USGS2021

The Alaska convergent margin backstop splay fault zone, a potential large tsunami generator between the frontal prism and continental framework

The giant tsunami that swept the Pacific from Alaska to Antarctica in 1946 was generated along one of three Alaska Trench instrumentally recorded aftershock areas following great and giant earthquakes. Aftershock areas were investigated during the past decade with multibeam bathymetry, ocean bottom seismograph wide‐angle seismic, reprocessed legacy, and new seismic reflection images. Summarized and updated here are previous papers and additional data. Tectonic structures collocated with aftershock area boundaries indicate possible lengths of rupture in future great earthquakes. NE aftershock area boundaries relate to subducted lower plate structures whereas the SW zone upper plate retains Beringian structural relicts. The lower to middle slope transition separating a stronger continental framework rock from a weaker accreted prism occurs along splay fault zones previously interpreted as backstops in seismic images. Damage zones along splay faults are generally 1‐km‐wide dipping typically 21°. Splays form slip paths from the plate interface to the seafloor much shorter than the 3–4° dipping plate interface beneath the frontal prism. Associated seafloor vent structures indicate overpressured fluids at depth. Splay fault dip and its rigid hanging wall impart greater seafloor uplift than the accreted prism per unit of slip making them effective tsunami generators. Backstop splay fault zones (BSFZs) run along the entire Alaska Trench. Beneath the frontal prism, active bend faults add rugosity to the plate interface and km high relief is commonly imaged in reprocessed legacy and new seismic data. The 1946 Unimak great (M8.6) earthquake epicenter is located near the BSFZ.

Alaska↗
USGS2020

Structural controls on slope failure within the western Santa Barbara Channel based on 2D and 3D seismic imaging

The Santa Barbara Channel, offshore California, contains several submarine landslides and ample evidence for incipient failure. This region hosts active thrust and reverse faults that accommodate several mm/yr of convergence, yet the relationships between tectonic deformation and slope failure remain unclear. We present 3‐D and 2‐D multichannel seismic reflection (MCS) data sets, multibeam bathymetry, and chronostratigraphic constraints to investigate the controls on slope failure. Splay faulting along the North Channel Deformation Trend (NCDT) coincides with a distinct zone of compressional uplift and onlapping of steeply dipping Quaternary strata. The NCDT is spatially correlated with seafloor fissures, and 3‐D seismic analyses reveal an intricate system of en echelon reverse faults that offset sediments younger than ~25 ka. Localized uplift zones are located between faults, one of which underlies the Gaviota landslide headscarp. We observe a direct relationship between slope failure and along‐strike variations in the tectonostratigraphic framework. Based on geophysical properties at Ocean Drilling Program (ODP) Site 893, we predict a trend in compaction and porosity reduction in the basin that drives pore fluids up‐dip, toward the zone of onlap above the NCDT, thus reducing slope stability. This interplay between tectonic, sedimentary, and fluid‐flow processes along the NCDT has created a confluence of preconditioning factors, with Gaviota and Goleta landslides being distinguished from the surrounding slopes by their position above the NCDT. The distribution of seafloor fissures suggests sections of the slope remain unstable and are prone to future landsliding. These results provide insights into the processes and 3‐D feedbacks that lead to slope instability along other convergent margins.

California↗
USGS2023

The Mojave section of the San Andreas fault (California), 1: Shaping the terrace stratigraphy of Littlerock Creek through the competition between rapid strike-slip faulting and lateral stream erosion over the last 40ka.

To determine the post-40 ka slip-rate along the Mojave section of the San Andreas Fault (MSAF) we re-analyze the sedimentary record preserved where Little Rock (LR) Creek flows across the fault. At this location, interaction between the northeast-flowing stream and right-lateral fault has resulted in the abandonment and preservation of 11 strath terraces and one paleo-floodplain in the downstream trailing corner of the river, two of which are also preserved upstream to provide cross-fault matches. A new model of fault-induced river deflection, together with standard terrace riser restoration, yields strike-slip displacements of 1,140 ± 160 m for the older terrace and 360 ± 70 m for the younger one. When combined with new 10 Be dating and reinterpretation of prior measurements the displaced terraces yield right-lateral slip-rates of 27.7 +6.9/−3.5 and 26.8 +3.4/−3.0 mm/yr over the last 23 k.y. and last 40 k.y., where uncertainties are at 95% credible intervals. These new rate determinations are consistent with independent late Holocene estimates, indicating that the long-term rate of strain accumulation along the MSAF is relatively fast and does not vary significantly when averaged over timescales of 15–20 k.y. Using our new model of stream deflection, we find that the fluvial sequence was emplaced in two distinct periods, each characterized by a temporally stable but markedly different deflected river geometry. Each period coincides with a distinct stage of erosive power along LR Creek determined from independent paleoclimate proxies. Importantly, application of the new river-deflection model allows strike-slip displacements to be determined in the absence of upstream piercing points.

California↗