Search USGS⌕ Search

SEARCH · Search USGS

Results for “Chemical Geology”

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.

238 records · Page 5Linked to original sources
USGS2026

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li 2 O, respectively, for global Li pegmatite deposits ( n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs 2 O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb 2 O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Ore Geology Reviews↗
USGS1986

Regionalization of low-flow characteristics of Tennessee streams

Procedures for estimating 3-day 2-year, 3-day 10-year, 3-day 20-year, and 7-day 10-year low flows at ungaged stream sites in Tennessee are based on surface geology and drainage area size. One set of equations applies to west Tennessee streams, and another set applies to central and east Tennessee streams. The equations do not apply to streams where flow is significantly altered by activities of man. Standard errors of estimate of equations for west Tennessee are 24 to 32% and for central and east Tennessee 31 to 35%. Streamflow recession indexes, in days/log cycle, are used to account for effects of geology of the drainage basin on low flow of streams. The indexes in Tennessee range from 32 days/log cycle for clay and shale to 350 days/log cycle for gravel and sand, indicating different aquifer characteristics of the geologic units that sustain streamflows during periods of no surface runoff. Streamflow recession rate depends primarily on transmissivity and storage characteristics of the aquifers, and the average distance from stream channels to basin divides. Geology and drainage basin size are the most significant variables affecting low flow in Tennessee streams according to regression analyses. (Author 's abstract)

Tennessee↗
USGS2024

Field evidence and indicators of rockfall fragmentation and implications for mobility

Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.

Albacete province, Lombardy and Aosta Valley, Yose↗
USGS2026

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana↗
USGS2024

Insights on gas hydrate formation and growth within an interbedded sand reservoir from well logging at the Qiongdongnan Basin, South China Sea

Although variable well log resolution and its control on saturation estimation has been studied, it has not been directly applied to a specific location to explore the nature of gas hydrate within a sand reservoir. We applied in-situ measurements of resistivities, neutron porosity, and gamma ray at two sites in the Qiongdongnan Basin, South China Sea (QDN-W05–2021 and QDN-W08–2021) to investigate the reservoir parameters of a hydrate-bearing sand reservoir. Our results show that gas hydrate is distributed in 5 zones with a total thickness of 10.7 m and an average saturation of 69% at the QDN-W05–2021 site, while they are distributed in 2 zones with a total thickness of 4.3 m and an average saturation of 49% at the QDN-W08–2021 site. We found that variances in saturations estimated from lateral-extra deep button (RX), phase shift (P40H-P40L), and attenuation (A40H-A40L) resistivities within the laterally mapped continuous sand body were affected by the nature of gas hydrate occurrences. Results indicate gas hydrate forms and accumulates at the center of the sand layer and tends to be less or not present toward the top and base. Integrated with seismic data, the in-situ measurements provide insights in the evolution of a mushroom-shaped, hydrate-gas reservoir system. In the system, free gas is likely horizontally transported from the top-center of the gas chimney to the surrounding areas in the early stage dominated by a warm-gas environment, whereas hydrate forms in the opposite pathway starting from the surrounding areas in the following stage with temperature reducing. Our study suggests that high-resolution in-situ measurements not only are a tool to identify the physical properties, but also can be used to help explain the physical process of hydrate growth and accumulation.

Marine Geology↗
USGS2025

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

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

New York↗
USGS2025

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

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

International Journal of Coal Geology↗
USGS2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater-dependent ecosystems in areas with industrial land use are at risk of exposure to a PFAS chemicals. We investigated one such system with several known PFAS source areas, where high and low permeability sediments (glacial) coupled with groundwater-lake and groundwater/surface-water interactions created complex ‘source to seep’ dynamics. Using heat-tracing and chemical methods, numerous preferential groundwater discharge zones were identified and sampled across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, downgradient of Joint Base Cape Cod (JBCC). Surface-water and groundwater samples were analyzed for 40 PFAS compounds between March and October 2022. Samples were collected from groundwater seeps identified as preferential discharge points (PDPs), wells upgradient of the stream-wetland system, contributing flow-through kettle lakes, and along Quashnet River surface-waters. PFAS from sampled waters contained perfluorinated carboxylates (PFPea, PFHxA, PFNA), perfluorinated sulfonates (PFBS, PFPeS, PFHxS, PFOS), fluorotelomer sulfonates (6:2, 8:2 FtS), and perfluoroalkyl sulfonamides (PFOSA). Samples from PDPs and wells had measured PFAS concentrations ranging from non-detect to 4677 ng/L ng/L (mean = 418 ng/L, std. = 709 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to groundwater flowpath emergence at PDPs. Correlations ( p < 0.01) between deuterium excess, %PFAS precursors, and terminal PFAS compounds highlighted potential precursor transformations associated with lake-groundwater exchange along flowpaths sourcing PDPs. However, some seepages had higher total PFAS concentrations (>1000 ng/L) than upgradient kettle lakes despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. Results from these synoptic surveys address gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations and mass loading in preferential groundwater discharge zones.

Massachusetts↗
USGS2025

Solution-collapse breccia pipe uranium deposits of the southern Colorado Plateau, northwestern Arizona, USA

Some of the highest-grade uranium deposits in the United States occur in breccia pipes that formed by solution and collapse of sedimentary strata, which occur in the southern portion of the Colorado Plateau in northwestern Arizona. The host breccia pipes are up to 1200 m in vertical extent, average about 90 m in diameter, and can cross-cut strata from their base in the Mississippian Redwall Limestone to as stratigraphically high on some plateaus as the Triassic Chinle Formation. These uranium-base metal deposits are up to 600 m thick and formed within the breccia pipes where they transect the Permian Coconino Sandstone, Hermit Formation, and the Esplanade Sandstone. Of the hundreds of breccia pipes identified across this region, only a small percentage are known to contain mineralization. The main uranium ore mineral is uraninite that is intergrown with at least 20 base-metal sulfide minerals, which contribute Fe, Cu, Co, As, Pb, Zn, Ni, and Ag to the deposits. This study considered regional stratigraphy, sulfur isotope systematics, mineralogy, in situ dating, and compilation and analysis of previous work on the deposits. A comprehensive deposit model has not been published for these deposits. This analysis identified new additions to update the deposit model for these unusual, possibly unique deposits. Proposed modifications to the model include: (1) the source, mechanisms, timing of the base-metal sulfide mineral assemblages, and (2) the source, mechanism, and timing of the uranium mineralization. Sulfide and uranium deposition are shown to be separate mineralization events. The study proposes the possible role of gypsum as a source of sulfur for the sulfide minerals in the deposits. Groundwaters carrying uranium encountered the preexisting sulfides in breccia pipes, reducing the uranyl ions, and precipitating U oxide (as uraninite). Analysis of the regional stratigraphy recognized that numerous beds of gypsum are in the strata that lie only tens of meters above the breccia pipe deposits. In the breccia pipe region, if these stratigraphic units (Toroweap and Kaibab Formations) do not contain gypsum layers then the underlying pipes are not mineralized; where these Permian gypsum layers do occur, breccia pipes can host mineralization. This new understanding should be useful in identifying the prospective region for mineralized pipes.

Arizona↗
USGS2025

Confocal laser-scanning microscopy (CLSM)-based thermal maturity of Tasmanites and progress in standardization of fluorescence microspectrometry

Evaluation of thermal maturity in vitrinite-free or vitrinite-deficient sediments via fluorescence microspectrometry can provide relevant information related to petroleum exploration and thermal history assessment. However, variation in spectral fluorescence properties of alginite macerals with increasing thermal maturity is largely underexplored. Here, authors of this study have applied confocal laser-scanning microscopy (CLSM) in conjunction with fluorescence microspectrometry to a maturity series of marine Upper Devonian Tasmanites algae from the Ohio Shale (Huron Member) and a single sample from the Marcellus Formation of the Appalachian Basin. Spectral fluorescence properties of Tasmanites were evaluated in relation to orientation, measurement location, and the number of measurements per sample, and were compared to published literature. Emission spectra of Tasmanites from continuous wave laser excitation (405 nm) were acquired from sections perpendicular and parallel to bedding and at different positions within individual Tasmanites bodies. The results showed a progressive red-shift in emission maxima (λ max ) in a large sample sized maturity series ( N = 19), e.g., 493 to 578 nm for the perpendicular section at middle position. Further, blue-shifted apex and mineral-adjacent positions within sections perpendicular to bedding were observed, with the latter being reported here for the first time. While blue-shift at apex positions can be attributed to mechanical deformation-induced reorientation of photoselected fluorophores, the blue-shifted mineral-adjacent positions could result from strain loading and development of a plastic deformation region at the mineral contact zone with Tasmanites . A decrease in standard deviation with increasing number of measured emission maxima is well-observed, and 15 to 20 individual measurements per sample appears sufficient for low standard deviation and coefficient of variance. CLSM-derived thermal maturity parameters indicated that a moderate positive correlation of red/green quotient (Q 650/500 ; R 2 = 0.67) with solid bitumen reflectance (BR o in %) exists. For reproducible results, the determination of λ max and Q 650/500 should be conducted exclusively in the middle position at perpendicular and parallel sections of the polished whole-rock pellets, where the lowest standard deviation in measurement was observed. These results strengthen the suitability and relevance of the CLSM technique in thermal maturity studies of dispersed organic matter (DOM) and contribute to the standardization of fluorescence microspectrometry methods in organic petrology investigation.

Kentucky, Ohio, Virginia, West Virginia↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗
USGS2025

Assimilation of reduced carbon triggers platinum alloy saturation in mafic and ultramafic magmas

It is generally observed that magmatic sulfide ores have higher ratios of Pd/Pt than the mantle-like values of their parental magmas. This discrepancy has defied simple explanation because the partitioning behavior of both elements between sulfide and silicate liquids is very similar. Assimilation of sulfur- and carbon-rich country rocks by mafic and ultramafic magmas is considered a critical, if not essential, step in the formation of magmatic base metal sulfide deposits. Although there is general consensus that the assimilation of external sulfur and carbon promotes sulfide saturation, the effect of carbon assimilation on the solubilities of platinum-group elements in natural S-bearing silicate melt has been overlooked. In this study, we investigate the variations of platinum and palladium solubilities during assimilation of graphite and methane through thermodynamic modeling, in comparison with data from an array of highly distinctive magmatic sulfide ore systems representing ages from Archean to Paleozoic, melt compositions from komatiite to basalt, and magmatic settings including lavas, hypabyssal intrusions, plutonic continental arc roots, and plutonic layered intrusions, namely: Raglan, Norilsk-Talnakh, Lac des Iles, and the J-M Reef of the Stillwater Complex. We model assimilation-fractional crystallization processes to estimate the reduction of oxygen fugacity ( ⁠fO2⁠ ) of the melt due to incorporation of graphite and methane. The simulations show that although Pd remains highly soluble during the progressive assimilation of reduced carbon, Pt solubility decreases significantly as the silicate melt becomes increasingly reduced. With less than 8 % of sediment assimilation, Pt alloy may saturate and then deviate from sulfide-undersaturated silicate melts, concomitantly increasing the Pd/Pt value of the remaining melts of the Raglan and Norilsk-Talnakh systems. For the Lac des Iles and Stillwater systems, a higher extent of assimilation is needed to reach Pt saturation because of the relatively carbon-poor nature of the lower crustal rocks. The assimilation of methane volatiles is shown to be more effective than graphite assimilation, and it provides a pathway to Pt alloy fractionation in the absence of detectable amounts of bulk host-rock assimilation. High Pd/Pt values have been documented in many world-class magmatic sulfide deposits whose parental magmas have demonstrably experienced crustal contamination. Our model suggests that although anomalous Pd/Pt values may be explained by other mechanisms such as incongruent melting of preexisting sulfide or differences in the diffusivities of the metals within achieving equilibration, the assimilation of graphite or methane may play an important role in the global occurrence of magmatic sulfide ores with elevated Pd/Pt values.

Economic Geology↗
USGS2024

Complex landslide patterns explained by local intra-unit variability of stratigraphy and structure: Case study in the Tyee Formation, Oregon, USA

Lithology and geologic structure are important controls on landslide susceptibility and are incorporated into many regional landslide hazard models. Typically, metrics for mapped geologic units are used as model input variables and a single set of values for material strength are assumed, regardless of spatial heterogeneities that may exist within a map unit. Here we describe how differences in bedding thickness, grain size, inferred uniaxial compressive strength, and bedding dip control the inherent susceptibility of slopes to deep-seated failure within a single mapped geologic unit - the Tyee Formation of Oregon, USA. The Tyee, which covers over 15,000 km2 and underlies much of the Oregon Coast Range, comprises gently folded alternating beds of sandstone and siltstone deposited as turbidites, forming a 2-km thick Eocene submarine fan which has been uplifted and exhumed through the Cenozoic. Deep-seated landslides are widespread in the Tyee, but form a complex spatial pattern such that landslide density ranges from 0 to 24% of the total landscape area. These slides are often extensive and sufficiently deep to reduce local hillslope gradients, resulting in a strong negative correlation between landslide density and mean local slope. Mean annual precipitation and predicted strong ground motions from Cascadia earthquake scenarios also fail to explain the spatial distribution of deep-seated landslides. Consequently, landslide stability models, which are strongly influenced by landscape slope, pore-water pressure, and seismic acceleration, yield landslide susceptibility maps which are broadly anti-correlated with mapped deep-seated landslide density. Through a multivariable linear regression model, we show that much of the variance in deep-seated landslide density can be explained by variability of intra-unit stratigraphic and structural characteristics, which we measure at 128 sites across two study areas totaling ∼3000 km2. Our results suggest bedding dip is only weakly correlated to landslide density, but strongly influences landslide failure style. Subtle increases in bedding dip, even in the gently folded Tyee Formation, result in a substantially higher likelihood of a landslide being cataclinal, or parallel to bedding. Overall, we find a slight majority of landslides fail within these cataclinal slopes, and that these landslides tend to be larger than non-cataclinal landslides. We also show that the lithological and structural properties that influence landslide susceptibility are distinct for these two populations of landslides. Our results demonstrate how localized, intra-unit, geologic variability can exert strong control on landslide susceptibility and failure style. This suggests that in some locations, landslide hazard models could be significantly improved by incorporating detailed, spatially variable, geologic properties rather than relying solely on generalized geologic map units.

Oregon↗
USGS2024

Geologic map and structure sections along the southern part of the Bartlett Springs Fault Zone and adjacent areas from Cache Creek to Lake Berryessa, northern Coast Ranges, California

Introduction Located in the Coast Ranges of northern California, the Bartlett Springs Fault Zone is the easternmost fault in the San Andreas Fault system in northern California. The fault is a right-lateral, strike-slip structure considered capable of producing an earthquake of moment magnitude 7. The purpose of this mapping is to better characterize the geology and earthquake hazards associated with the southern part of the Bartlett Springs Fault Zone and to help identify any evidence of active uplift on the faults bounding the Coast Ranges. Although the area immediately surrounding the Bartlett Springs Fault Zone is sparsely populated, its southern segment presents a potential seismic hazard to northern California communities as far away as the San Francisco Bay region and Sacramento. There are also nearby water resources, mineral resources, and public lands used for public recreation. The Coast Ranges of northern California are a series of northwest-southeast-oriented mountain ranges and valleys located north of the San Francisco Bay region, between the Pacific Ocean to the west and the Sacramento Valley to the east. The region has rugged terrain, high mountain peaks that reach more than 2,400 meters above sea level, isolated and narrow valley bottoms on which most human settlements are located, and large drainage systems that tend to follow the northwest-southeast-oriented topographic grain. The physiographic character of the region is shaped by its bedrock geology, deformational history, and active faulting. The basement rocks of the northern Coast Ranges consist of the Franciscan Complex and the Great Valley complex, the latter of which consists of two informal units, the Coast Range ophiolite and the Great Valley sequence. The Franciscan Complex and the Great Valley complex are in structural contact along the Coast Range Fault, a regional-scale structure and fundamental crustal boundary. The Franciscan Complex and the Great Valley complex are superposed by active, northwest-southeast-striking strike-slip faults that are associated with seismicity swarms. These active strike-slip faults can produce moderate to large earthquakes that have moment magnitudes of 7–8. In places, these active structures bound large ranges and valleys, suggesting that much of the modern topographic expression is the result of active deformation processes. This report contains new 1:24,000-scale geologic mapping along the southern part of the Bartlett Springs Fault Zone between Clear Lake and Lake Berryessa. The map area spans 738 square kilometers in northern Napa County, southern Lake County, and parts of Yolo and Colusa Counties. The south and east borders of the map are 90 kilometers north of San Francisco and 70 kilometers west of Sacramento, respectively. The map area is within the Knoxville mining district, which has a history of mercury and gold mining dating back to the mid-19th century. The two main towns in the region, Lower Lake and Clearlake, California, are west-northwest of the map area. Approximately 71,000 people live in the cities and rural communities located within a 40-kilometer radius of the center of the map area. The bedrock geology, cross sections, and structural data presented herein are critical for evaluating the long-term evolution of the Bartlett Springs Fault Zone. This work will supplement studies on local seismic hazards, liquefaction potential, landslide hazards, earthquake geology, natural resources, groundwater resources, engineering geology, and tectonic history by providing the background information for site-specific investigations on these subjects.

California↗
USGS2025

Variable partitioning of lithium in rhyolitic melt during decompression and ascent

The partitioning behavior of Li in magmatic systems is increasingly being investigated due to the economic importance of Li in the transition to sustainable energy resources (e.g., batteries). However, at upper crustal pressures, it remains uncertain whether Li preferentially partitions into the vapor or liquid (brine) phase or remains in the silicate melt. This complicates our ability to determine where Li resides—silicate melt, minerals, or fluid phase—upon eruption, a crucial factor for understanding its postdepositional movement and concentration into a brine or volcano-sedimentary deposit. Here, we present a novel investigation into the behavior of Li within natural evolved melts during continuous magma decompression and ascent using melt embayments (open melt inclusions). Mineral-hosted melt embayments preserve records of the evolving composition of the exterior melt, including degassing pathways and ascent timescales, when paired with appropriate diffusion coefficients. Lithium concentration profiles were measured in quartz-hosted melt embayments from the rapidly quenched eruptive phases of five rhyolitic, caldera-forming eruptions to investigate the behavior of Li during magma decompression and ascent, where vapor partitioning and ascent dynamics were previously established by investigating H 2 O and CO 2 profiles. We find that in four systems, embayments contain lower interior Li concentrations than the coerupted melt inclusions; the fifth system contains the same Li concentrations in embayments and melt inclusions. However, many of these embayments contain gradients, with 84 % preserving Li enrichment near the melt-bubble interface, as compared to their interior concentration. We interpret these characteristics to represent two distinct stages of Li partitioning during magma decompression and ascent, in contrast to existing literature that proposes only one type of partitioning behavior. The first stage is interpreted as melt depletion of Li, likely driven by partitioning into an exsolved supercritical fluid phase, supported by the strong correlation between the extent of Li depletion and Cl concentration in the melt, as well as the decompression rate. This behavior then fundamentally shifts, where Li reenriches in the melt, postulated to be driven by the unmixing of the supercritical fluid phase at shallow pressures. For the one system that did not develop Li gradients through decompression, we attribute this to the lower values of Na and Cl in the melt, potentially inhibiting the partitioning of Li into a fluid phase. Importantly, the behavior of Li during decompression is not consistent within or between volcanic centers, highlighting the need for systematic experimental investigation in variable composition melts at pressures relevant to conduit dynamics. This knowledge would improve our ability to model Li profiles to understand magma decompression, and predict where Li resides (e.g., stored in volcanic glass, gas, or crystals) upon eruption prior to any later extraction.

Economic Geology↗
USGS2026

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10 6 m 3 of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H 2 O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Hawaii↗
USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗