Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of the Mechanics and Physics of Solids”

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,261 records · Page 2Linked to original sources
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↗
USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗
USGS2026

Statewide surficial geologic map of Nebraska underscores Quaternary landscape evolution from the High Plains to the Central Lowland

Surficial geologic mapping in Nebraska has been conducted primarily at the 1:24,000 scale since the mid-1990s, although there have also been limited efforts to map generalized Quaternary and surficial geology within the state and the region. We compiled and evaluated disparate maps (1:24,000–1:1,000,000) and datasets—including geologic, soil and soil parent material, and geomorphic maps as well as LiDAR derivatives—to produce a single 1:500,000 scale surficial geologic map of Nebraska that is the first of its kind. This new map provides a coarse-scale surficial geologic map that will be incorporated into a nationwide U.S. Geological Survey Quaternary geologic map. It also reflects the variation and uniqueness of physical landscapes in the state, where the Great Plains and Central Lowland physiographic provinces meet, further developing a richer interdisciplinary understanding of regional geomorphology in the heart of North America.

Nebraska↗
USGS2026

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

California↗
USGS2025

Land application of drill waste: A scope analysis

Drilling fluid waste land application, a process where drilling wastes are spread and tilled into the land surface, has become common in some petroleum-producing states, however, the potential benefits and risks of this practice are not well studied. Drilling fluids can be water- or oil-based and can have high concentrations of total soluble salts and total petroleum hydrocarbons. Comprehensive chemical characterization of these fluids is not well documented in the literature, and the extent of land application is largely unknown. We hypothesized that the land application of drill waste would fluctuate over time due to economic factors. To begin to understand the extent of historical and potential future land application, we analyzed data from over 5,800 drilling fluid land application permits collected by the Oklahoma Corporation Commission for years 2000, 2005, 2010, and 2015–2020. During the years studied, drilling fluid wastes were applied to more than 250,000 acres in Oklahoma, with over 54,000 thousand barrels (Mbbl) of liquids and nearly 21,000 Mbbl of solids applied. Land application is widespread (occurring in 59/77 counties), however recent drilling activity, land availability, and the economics of transportation have created conditions favorable for land application specifically in the Anadarko Basin. Land application can co-occur with sensitive areas, such as important groundwater and surface-water drinking sources and agricultural fields used for subsistence or feed crop production. Our approach for quantifying the extent of land application, along with further chemical characterization studies, can aid operators and land managers who are considering this practice in assessing the associated benefits and risks.

Oklahoma↗
USGS2025

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal↗
USGS2026

Deaf, deafblind, and hard of hearing university student experiences with earthquake early warning in the United States: Evaluating language planning and technology access

The growing literature on deaf and hard of hearing (DHH+) populations and disasters demonstrates that emergency communication (including alerts) is not reaching global DHH + individuals with dangerous impacts for morbidity and mortality. This is the first research study in the U.S. to qualitatively explore the experiences of DHH + persons with earthquake early warning (EEW) through group-based dialogue sessions. The study investigates eight DHH + university students'past earthquake experiences, access to EEW alerts, and perceptions of ShakeAlert Ⓡ , an EEW system for detecting earthquakes and alerting residents of California, Oregon, and Washington. Findings highlight key gaps in disaster alert usability within four thematic areas: lack of messaging in participants' language(s), unclear alert messaging, deficient message delivery mechanisms for deafblind persons, and insufficient access to earthquake information and training that leads to dependence on informal information networks. Weaknesses identified in these four themes reduce DHH + trust in EEW systems and compromise the capacity of alert recipients to take swift protective action or to mentally prepare before shaking starts. The study also underscores structural factors such as insufficient linguistic representation in disaster language planning and technology design, which ignores the linguistic and sensory access needs of DHH + individuals. Building on disaster language planning frameworks, we recommend involving DHH + populations to co-develop EEW alerts. By centering DHH + perspectives, this research contributes to ongoing efforts to ensure that EEW systems reach everyone.

California↗
USGS2025

Expression of corticoid-regulatory genes in the gills of Atlantic salmon (Salmo salar) parr and smolt and during salinity acclimation

In teleost fishes, cortisol is the major corticoid and has both glucocorticoid and mineralocorticoid actions. However, how fish tissues discriminate between these distinct corticosteroid actions is unclear. In mammals, the major factors responsible for intracellular corticosteroid regulation are glucocorticoid receptors ( gr s) and the mineralocorticoid receptor ( mr ), but their role in osmoregulation of fish is unclear. 11β-hydroxysteroid dehydrogenases ( hsd11b s) control the levels of intracellular corticosteroids by converting from bioactive forms to inert forms. To investigate how Atlantic salmon ( Salmo salar ) respond to cortisol in different physiological or environmental conditions, we performed comparisons of parr and smolt, and osmotic challenge experiments to examine the physiological responses and gill transcript levels of genes underlying cortisol-signalling, including gr1 , gr2 , mr , hsd11b2 and hsd11b3 . Because cortisol may interact with growth hormone and prolactin during salinity changes, transcript levels encoding growth hormone receptors ( ghr1 , ghr2 ) and the prolactin receptor ( prlr ) were also examined. Hsd11b2 transcript levels in seawater-acclimated fish were consistently lower compared to fish acclimated to fresh water. After transfer to seawater, prlr transcript levels in fish significantly decreased and transcript levels of ghr1 , ghr2 and hsd11b3 showed no change or were slightly higher than those of freshwater control groups. Gr1 , gr2 and mr transcript levels were slightly but consistently higher in fish acclimated to fresh water relative to seawater. Our results indicate that changes in corticosteroid receptor and hsd11b2 transcript levels in the gills may be important mechanisms that regulate corticoid signals to achieve ion homeostasis in Atlantic salmon.

Journal of Fish Biology↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS1999

Water quality of Rob Roy Reservoir and Lake Owen, Albany County, and Granite Springs and Crystal Lake Reservoirs, Laramie County, Wyoming, 1997-98

The water quality of four reservoirs was assessed during 1997 and 1998 as a cooperative project between the Cheyenne Board of Public Utilities and the U. S. Geological Survey. The four reservoirs, Rob Roy, Lake Owen, Granite Springs, and Crystal Lake, provide approximately 75 percent of the public water supply for Cheyenne, Wyoming. Samples of water and bottom sediment were collected and analyzed for selected physical, chemical, and biological characteristics to provide data about the reservoirs. Water flows between the reservoirs through a series of pipelines and stream channels. The reservoirs differ in physical characteristics such as elevation, volume, and depth. Profiles of temperature, dissolved oxygen, specific conductance, and pH were examined. Three of the four reservoirs exhibited stratification during the summer. The profiles indicate that stratification develops in all reservoirs except Lake Owen. Stratification developed in Rob Roy, Granite Springs, and Crystal Lake Reservoirs by mid-July in 1998 and continued until September, with the thickness of the epilimnion increasing during that time. Secchi disk readings indicated Rob Roy Reservoir had the clearest water of the four reservoirs studied. The composition of the phytoplankton community was different in the upper two reservoirs from that in the lower two reservoirs. Many of the species found in Rob Roy Reservoir and Lake Owen are associated with oligotrophic, nutrient-poor conditions. In contrast, many of the species found in Granite Springs and Crystal Lake Reservoirs are associated with mesotrophic or eutrophic conditions. The total number of taxa identified also increased downstream. The chemical water type in the reservoirs was similar, but dissolved-solids concentrations were greater in the downstream reservoirs. Water in all four reservoirs was a calcium-bicarbonate type. In the fall of 1997, Rob Roy Reservoir had the lowest dissolved-solids concentration (19 milligrams per liter), whereas Crystal Lake Reservoir had the highest concentration (63 milligrams per liter). Relatively little differences in the concentrations of major-ion species were noted between samples collected near the surface and near the bottom of the same reservoir. In contrast, iron and manganese concentrations generally were higher in samples collected near the bottom of a reservoir than in near-surface samples collected from the same reservoir. Composite bottom-sediment samples from all four reservoirs contained similar concentrations of bulk constituents such as aluminum, iron, phosphorus and titanium, but varied in concentrations of trace elements. Trace-element concentrations in Rob Roy Reservoir and Lake Owen were similar to the crustal average, whereas in Granite Springs and Crystal Lake Reservoirs the concentrations were similar to granitic rocks.

Wyoming↗
USGS2023

3D Dynamic rupture modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 earthquake doublet using early observations

The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The M w 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral M w 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data‐driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

The Seismic Record↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2026

Bascom Laser Diffraction Sedimentology Laboratory, Reston, Virginia

Introduction At the Bascom Laser Diffraction Sedimentology Laboratory, which is located in the Florence Bascom Geoscience Center at U.S. Geological Survey (USGS) headquarters in Reston, Virginia, scientists use physical sedimentology and particle characterization techniques to conduct detailed sediment characterization. Scientists address research problems in collaboration with other USGS science centers, State geological surveys, commercial industry, universities, and other partners. Laboratory capabilities include laser diffractometry for quantitative particle-size analysis, portable x-ray fluorescence (XRF) analysis for determining elemental abundances in rock or sediment samples, petrographic analysis of geologic media, and mechanical sieve analysis. These methods are used to analyze soil and sediment core material from terrestrial, marine, and lacustrine environments, surface sediments from coastal regions, and calcareous materials. Work done by the laboratory supports geologic mapping, resource assessments, land change studies, and geohazard analyses.

Fact Sheet↗
USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗
USGS2026

Deep subsurface organic-rich shale supports abundant, diverse, and novel fungi

As Earth’s principal reservoir of organic carbon and microbial biomass, the deep subsurface hosts microorganisms capable of mobilizing this once-sequestered carbon. Contrary to standard assumptions of eukaryotic scarcity, this study documents abundant fungal communities, ranging from 4.2 × 10 3 to 6.8 × 10 3 fungal cells mL −1 , across a methane-producing organic-rich shale 247–556 meters below the surface. Although fungal:bacterial cell ratios ranged from 1:7028 to 1:713, application of biomass conversion factors developed for oceanic systems yielded a median fungal:bacterial biomass ratio of 1:4.7. 16S rRNA gene amplicons revealed bacterial and archaeal communities mirroring those found in well-characterized extremophilic, carbon-degrading environments, while sequencing of 18S rRNA gene and ITS rRNA spacer amplicons collectively identified a eukaryotic hotspot with 689 fungal OTUs across six phyla. The dominant fungal classes, Agaricomycetes and Dothideomycetes, are well-established degraders of recalcitrant carbon compounds at the surface, suggesting they may similarly contribute to organic matter degradation and ecosystem maintenance in the subsurface. Cultivation and isolation efforts yielded 205 fungal strains, including 13 candidate novel taxa, underscoring the deep subsurface as an underexplored eukaryotic habitat. Stable carbon isotopes indicate methane is predominantly generated via microbial conversion of the fossil carbon, while water isotopes suggest in situ geochemical conditions have been relatively stable since the Late Pleistocene, with subglacial recharge as a plausible mechanism for microbial introduction. Collectively, these findings suggest that fungi are underrecognized contributors to organic matter transformation and functional diversity in the deep biosphere, revealing a critical gap in our understanding of deep subsurface ecosystem processes.

Indiana, Michigan, Ohio↗
USGS2026

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2026

The petrogenesis of Þingmúli volcano, East Fjords, Iceland

In this work we revisit Þingmúli volcano (Þ = Th), a classic locality known as an example of a complete tholeiitic differentiation. Þingmúli is a ~ 9.5 Ma extinct central volcano located in the East Fjords of Iceland, in which the whole compositional spectrum from basalt to rhyolites have erupted. These volcanic products have been previously considered as petrogenetically related by an ideal fractionation trend, regardless any temporal relationship or volumetric considerations. Here we report new whole-rock geochemistry, mineral chemistry, isotope analyses, estimation of residence times of the different eruptive deposits, and an update of the original petrogenetic model. Our results highlight that an enriched source, likely spinel lherzolites, generated transitional-alkaline basaltic melts after 15–20% of partial melting at depths of 40–45 km. Many of these basaltic melts erupted at various stages of the volcano's history, while others remained longer in the volcanic plumbing system. These evolved by fractional crystallisation into basaltic andesite magmas with a residence time of ~5 years based on the crystal size distribution of the plagioclase population. Isotopic differences between the basalts/basaltic andesites ( 87 Sr/ 86 Sr ~ 0.7034; 143 Nd/ 144 Nd ~ 0.51315) and the erupted rhyolites ( 87 Sr/ 86 Sr ~ 0.7037; 143 Nd/ 144 Nd ~ 0.51304) indicate that the latter are not petrogenetically related to the former. Therefore, instead of a fractional crystallisation mechanism to generate the rhyolites, we propose the partial melting of ignimbrite layers located beneath the volcano. The broad range of trace element concentrations in andesites and dacites and their different isotopic values compared to the basalts strongly suggest that these magmas have been generated by magma mixing between basaltic and rhyolitic melts, similar to modern day Icelandic volcanoes such as Hekla. These results highlight the need to revisit previously studied Icelandic classic localities and reassess their traditionally proposed petrogenetic models.

Þingmúli volcano, East Fjords↗
USGS2025

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗