USGS Science⌕ Search

SEARCH · USGS Science

Results for “Climate of the Past”

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

Small-volume tephra deposits of the May 1924 explosions from Halemaʻumaʻu, Kīlauea volcano, and their origin

More than 50 explosive eruptions occurred from Halemaʻumaʻu at Kīlauea volcano over 17 days from May 11 to 27, 1924. Ballistics weighing as much as 14,000 kg were ejected and most landed within 2 km of the vent. Fine ash made up a major component of the tephra and was dispersed tens of kilometers downwind. Draining of the Halemaʻumaʻu lava lake occurred in late February 1924, with the crater floor eventually subsiding by a further ∼70 m (to ∼180 m below the crater rim) by the time the first explosions took place during the night of May 10–11. The largest explosions occurred on May 17–18 and smaller explosions continued until May 27, at which point Halemaʻumaʻu had more than doubled in width and depth. The explosions generated plumes reaching up to ∼10 km high with ballistics ejected up to 2 km from the crater. Almost 100 years later, we investigate and characterize the preserved tephra deposits within ∼3 km of the 1924 crater rim. Grain size and shape analyses were performed on 202 samples collected from 34 tephra profiles using dynamic image analysis, with a subset of layers from nine tephra profiles used for componentry (200 grains per layer in the 0.5–1 mm size fraction). Additionally, we characterize the average diameters (using the five largest clasts) at 216 locations and measure the average diameters of 2291 ballistics (largest per ∼100 m 2 area). Physical descriptions from fieldwork and grain size distributions were used to subdivide the tephra layers into five lithofacies: coarse homogeneous, fine homogenous, red ash, accretionary lapilli-bearing, and finely laminated. Grain size versus shape data show a range of values that demonstrate most grains are dense, smooth, and equant, in alignment with lithic clasts dominating the tephra componentry. The fine grained and accretionary lapilli-bearing nature of some of these lithofacies confirms that water influenced the style of the explosions. However, we also note juvenile clasts within many of the tephra layers, indicating that many of the layers were formed during phreatomagmatic explosions (sensu stricto), despite the eruptive mechanism being dominantly phreatic. Juvenile clasts are more abundant higher in the tephra profiles, suggesting that juvenile magma was more involved later in the explosive sequence. Thermal and hydrologic modeling indicate that groundwater inflow into a short-lived, small-diameter volcanic conduit (10-m to 120-m-diameter used for modeling) during the 78–85 days preceding the first explosion provides a physically plausible mechanism for this eruptive sequence.

Hawaii↗
USGS2026

Cambrian and Ordovician stratigraphy, conodont biostratigraphy, and microfacies analysis to support 1:24,000-scale geologic mapping of the southern Lake Champlain valley, New York and Vermont

Introduction Geologic mapping in the southern Lake Champlain valley of New York and Vermont (fig. 1) has required evaluation of stratigraphic nomenclature used since the early 1900s. The paleogeography of the Cambrian and Ordovician adjacent to the Adirondack Highlands, the high-relief mountains of Proterozoic igneous and metamorphic rocks of northern New York, has impacted the lithostratigraphic succession, facies changes, and several hiatuses that occur within the package of clastic and carbonate rocks that demonstrates overall deepening consistent from the southern to northern Appalachians. The stratigraphic nomenclature used for the geologic mapping of the Cambrian and Ordovician strata of the study area is derived from evaluation of names put forth by previous workers from areas of northern New York and western Vermont. Descriptions of type localities or type sections and other criteria set forth by the North American Stratigraphic Code (North American Commission on Stratigraphic Nomenclature [NACSN], 2021) were consulted. The North American Stratigraphic Code states that lithostratigraphic units, such as the fundamental unit, the formation, be defined by lithic character and be mappable (NACSN, 2021, Article 24). However, some early workers (Walcott, 1912; Kay, 1937; Fisher and Hanson, 1951; Oxley and Kay, 1959) have defined stratigraphic units based on fossil content, which may or may not be facies dependent, and hence may not be mappable as lithostratigraphic units. Although preservation of older names takes priority in stratigraphic nomenclature (NACSN, 2021, Article 7), some do not have adequately defined stratotypes or are located a distance away with different lithologies than that in the study area. Biostratigraphy can be a helpful tool for understanding correlations and facies changes. Herein, conodont biostratigraphy is utilized to date and correlate separate lithologic units. Further, microfacies analyses of units through inspection of petrographic thin sections provide paleogeographic information that helps in understanding the genesis of geologic units.

New York, Vermont↗
USGS2025

Relationship of basin structure and bedrock lithology to faulting in the 2019 Ridgecrest earthquake region, California, from gravity and aeromagnetic data

We investigate patterns of cumulative offsets on the faults that ruptured in 2019 and along the Garlock Fault in the Ridgecrest region, California using recently published gravity and aeromagnetic data. We also examine the relationship of basin structure and bedrock structure to the 2019 M7.1 Ridgecrest earthquake ruptures (Fig. 1A), which were primarily along a dextral northwest-striking fault system, and along a sinistral northeast-striking fault, which ruptured hours earlier with a M6.4 event.

California↗
USGS2026

Conventional hydrothermal power-producing systems of the Great Basin, USA

As part of the update to the electric-grade conventional hydrothermal assessment of the Great Basin, USA, Monte Carlo analyses of identified resources within explored regions will be performed to make estimates of discovered resources and associated uncertainty. Analyses use conditional statistics where estimates are conditioned upon a hydrothermal favorability map, allowing for the likelihood that more resources exist in regions of higher hydrothermal favorability. For these analyses, a dataset of identified hydrothermal systems is compiled, and the new compilation is described herein. Recognizing that a single hydrothermal system may be developed with multiple power plants, and that the hydrothermal upflow zone may be several kilometers across with many measurements characterizing a single hydrothermal system, a procedure was developed and employed to create clusters of points (power plants, measurements, etc.) that are associated with a single system, and a new central point was defined as the best estimator of the center of the hydrothermal system. Hydrothermal systems were uniquely identified by grouping electric-grade hydrothermal measurements and operating power plants within a distance of 10 km. Groups that are >10 km apart are assumed to be different electric-grade hydrothermal systems. While 10 km was used as the threshold, most systems were significantly further apart, and most points within groups were typically within 5 km of each other. A well measurement was considered an electric-grade measurement of a hydrothermal system if it had two properties: a measured temperature of >85 °C and evidence of hydrothermal convection. Other points that were added to the dataset are locations of operating powerplants or locations that have been classified as an electric-grade hydrothermal resource by either the U.S Geological Survey (USGS) or the Great Basin Center for Geothermal Energy. After all points are assigned to systems, new points were computed with the goal of identifying the center of the throat of the hydrothermal upflow zone. If operating powerplants exist for a system, then the arithmetic average of all power plant locations is used. Otherwise, if USGS made an estimate, that location is used. In the absence of both powerplants or USGS estimates, the arithmetic average of all electric-grade measurement locations is used. An example is shown of how these newly compiled locations might be ranked for uncertainty analyses, where higher confidence is assumed if measured temperature is higher and there are many supporting measurements indicating an electric-grade resource. In summary, 28 systems have operating power plants, an additional 78 systems are known identified electric-grade hydrothermal resources, and 99 new systems were identified as probable electric-grade systems with varying levels of confidence. These 205 locations are shown as a function of a recent hydrothermal favorability map, conceptually illustrating the conditional statistics that can be used to make estimates of the undiscovered resources of the Great Basin. An accompanying data release provides summaries of developed capacity by system and USGS estimates of likely total capacity and associated uncertainty.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2024

Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)

We developed Vp, Vs, Vp / Vs ratio, and Poisson’s ratio models of the uppermost crust (<4 km depth) from the eastern San Francisco (SF) Bay (California, USA) to near the Calaveras fault along a 15-km-long, linear profile. Upper crustal velocities are highly variable beneath, west, and well east of the Hayward fault. We observe eight notable features, from west to east: (1) Near San Francisco Bay, there is an ~2-km-wide structure with high Vp / Vs ratios (up to 5) and Poisson’s ratios (up to 0.48) extending from the surface to the base of our model, which we suggest the structure is a near-vertical fault that lies along a straight-line projection between the Silver Creek fault to the south and the Point Richmond fault to the north. The structure may be part of an ~90-km-long fault along the eastern SF Bay. (2) The western East Bay Plain, the lower lying area between the bay and the hills, includes up to 800 m of low-velocity sediments ( Vp ~1600–3000 m/s, Vs ~500 m/s to ~1000 m/s), underlain by higher velocity basement rocks ( Vp ~3000–5800 m/s; Vs ~1000–1500 m/s). (3) Between ~1 km and 3 km east of the Bay shoreline, sediments thin in a series of steps (likely faults) toward the Hayward fault. (4) Between ~3 km west and ~1 km east of the Hayward fault (at the East Chabot fault) at depths greater than 1 km, basement Vp (up to 6000 m/s) and Vs (up to 2800 m/s) are high, and Vp / Vs ratios (<2) and Poisson’s ratios (<0.3) are low, suggesting crystalline rocks. Furthermore, a near-vertical zone of low Vp / Vs ratios and Poisson’s ratios is between near-surface traces of the Hayward and East Chabot faults, likely corresponding to the San Leandro Gabbro of Ponce et al. (2003) . (5) Eastward of the East Chabot fault in the upper 1.5 km, basement Vp (~3000 m/s to ~4200 m/s) and Vs (~1200–2000 m/s) are lower than those west of the fault. (6) In the eastern Hayward/Oakland Hills, there are zones of laterally varying, high- and low-velocity ( Vp ~2500–3000 m/s) Jurassic–Cretaceous and Tertiary sediments in the shallow subsurface that likely extend much deeper than imaged. (7) Seismic energy that propagates westward from sources east of the Hayward fault (HF) appear weaker than energy that propagates eastward from sources west of the HF, suggesting that the HF acts as a partial barrier to shallow seismic energy propagation into the more populated eastern SF Bay area. (8) Unlike many fault zones, it appears that the active trace of the Hayward fault (in our study area) is not cored by a prominent, low-velocity zone relative to rocks to the east and west of the active trace. However, the active trace does mark a prominent change from relatively higher velocities to the west and lower velocities to the east.

California↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2026

Effects of a low-lipid diet on the gut microbiome and head kidney transcriptome of juvenile Chinook Salmon

Objective Pacific salmon Oncorhynchus spp. reared in production hatcheries are typically fed high-lipid, energy-dense diets to achieve large size and high body condition prior to release. In contrast, juveniles in natural environments tend to consume low-lipid, high-protein diets, and fish reared for research or conservation purposes are sometimes fed diets that are formulated to mimic natural diets and promote wild-like phenotypes. Understanding how these alternative diets affect fish health beyond growth and body condition could ultimately contribute to improving hatchery fish fitness. Methods In this work, we evaluated changes in the fecal microbiome and gene expression of juvenile Chinook Salmon O. tshawytscha on a standard high-lipid hatchery diet versus a low-lipid diet formulated to mimic the nutrition profile of natural-origin fish. To evaluate the time scale at which diet alters the fecal microbiome, we collected longitudinal samples over a 12-week period and switched the diets of a subset of fish twice during the experiment. We used 16S ribosomal RNA gene amplicon sequencing to characterize fecal microbiome differences between fish on the two diets as well as hatchery-reared fish at a production hatchery, hatchery fish that had been captured after release into a stream, and natural-origin, stream-reared fish of similar ages. Additionally, we conducted RNA sequencing on head kidney samples from laboratory-reared fish to evaluate changes in gene expression in this important immune organ. Results We found that the low-lipid diet and the hatchery diet resulted in microbiomes that differed from the microbiome of natural-origin fish and from each other and that diet-driven changes to the microbiome could occur in under 14 d. The low-lipid diet did not result in a microbiome that resembled the microbiome of naturally produced fish. Instead, the low-lipid diet resulted in a microbiome community that was distinct from those of fish reared on the hatchery diet and fish sampled from the wild. The RNA sequencing results indicated differential enrichment of pathways related to immunity, metabolism, and hormone synthesis between fish that were fed the two experimental diets. Conclusions The results suggest that additional environmental factors influence the microbiome more strongly than diet formulation or that the low-lipid diet has a smaller effect on the microbiome than a natural, ­invertebrate-based diet. Given that the gut microbiome and systemic immune function contribute significantly to disease resistance, our findings highlight the importance of understanding how diets fed to fish in captivity may affect fish health beyond growth and body condition metrics.

Journal of Aquatic Animal Health↗
USGS2026

Modeling future groundwater depletion to evaluate sustainability goals set under the Sustainable Groundwater Management Act in the critically overdrafted basins of the Central Valley, California, USA (2020–2070)

In 2014, California's Sustainable Groundwater Management Act (SGMA) mandated local agencies to devise and implement groundwater sustainability plans to address critically overdrafted conditions throughout the state's aquifers. However, the feasibility of these agencies' sustainability goals has not previously been assessed through a regional-scale, integrative lens. Here, we develop and analyze a novel, basin-wide database of 936 sustainability indicator wells located within Central Valley subbasins designated as critically overdrafted, most of which lie in the San Joaquin Valley. Our database shows 2040 groundwater elevation goals vary widely from 60 m above to 80 m below 2020 levels, with variability within and between adjacent subbasins. To evaluate the feasibility of achieving these goals, we coupled the database with a regional hydrologic model (Central Valley Hydrologic Model version 2) and simulated multiple future pumping scenarios. Results show that under increased groundwater demand, 60%–70% of indicator wells may fail to meet their 2040 goals. Even a 50% reduction from 2020 demand levels leaves nearly 40% of wells failing to meet their sustainability thresholds by 2040. Baseline models show that by 2070, up to 70% of wells could fail to meet their goals due to large-scale, spatially connected regions of groundwater depletion. This integrated framework, linking the first region-wide compilation of SGMA indicator wells with a regional groundwater model, demonstrates that many local sustainability goals may be unattainable with substantial (up to 50%) reductions in pumping. Additional management interventions, such as expanded recharge or coordinated demand reductions, may help achieve sustainability goals.

California↗
USGS2025

Sources of water and salts for the Zuni Salt Lake in west-central New Mexico

The Zuni Salt Lake is located in a maar in west-central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. There have been several investigations suggesting different sources for the water and salt to the lake. Springs, seeps, and ephemeral streamflow have all been observed to contribute freshwater to the lake, and brackish to hypersaline seeps have been documented along the banks of the lake. This report summarizes the findings of a study that characterizes the lake’s hydrology, its water and salinity sources, and the hydrogeologic conceptual model. Regional groundwater levels indicate that each of the aquifers in the area have the potential to discharge groundwater to the lake. There is also evidence of vertical groundwater flow pathways at the maar that were likely created by the igneous intrusion that fractured the intersecting aquifers. A detailed water budget was constructed from continuous lake stage, precipitation, and evaporation data to estimate the groundwater inflow to the Zuni Salt Lake. It was determined that groundwater inflow to the lake is 441 ±94 acre-feet per year, which composes as much as 77 percent of the total inflows. The high sodium and chloride concentrations measured in two hypersaline samples collected near the lake indicate that the majority of the dissolved solids entering the lake are from a hypersaline groundwater source. The geochemical and isotopic compositions measured in the lake and surrounding features support the interpretation that hypersaline groundwater is the primary source of salts to the lake, which is likely sourced from the older (and deeper) Permian units. The hypersaline groundwater samples collected during this investigation have a unique aqueous chemistry relative to each of the mapped aquifers, and variability in groundwater compositions is interpreted to result from differences in minerology and residence time.

New Mexico↗
USGS2025

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗
USGS1993

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

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

Tennessee↗
USGS2025

Geologic map of MTM −15032 and −20032 quadrangles, western Ladon basin, Mars

Our mapping efforts focused on the geomorphic terrains and geologic units contained within Mars Transverse Mercator (MTM) −15032 and −20032 quadrangles. These two quadrangles are located along the west side of Ladon basin and span lat −12.5° N. to −22.5° N. and long 325° E. to 330° E. The western part of Ladon basin and its bounding basin ring structures to the west preserved features that help to understand the long history of drainage across the Margaritifer Terra region of Mars. Our geologic map provides new insight into the extent of aqueous activity across western Ladon basin and the adjoining highlands, which includes the northern part of Ladon Valles, Arda Valles, numerous small valleys in the western highlands, and light-toned layered deposits associated with these fluvial features. The Mars Odyssey Thermal Emission Imaging System (THEMIS) infrared (IR) daytime mosaic (100 meters per pixel) was used as the primary base map. We constructed the geologic map of western Ladon basin at 1:1,000,000 scale. We identified 20 geologic units in the map area, which we divided into the following groups: crater units; volcanic units; chaotic units; basin fill units; crater, valley, and channel units; and plateau and highlands units. We observed the following sequence of events: (1) formation of Ladon and Holden impact basins in the middle Noachian, producing a mountainous unit; (2) landscape degradation and infilling during the Late Noachian to Late Hesperian forming terra units; (3) formation of Ladon Valles in the Late Noachian to Early Hesperian by catastrophic flooding, producing channel units; (4) accumulation of sediments in the Late Noachian to Early Amazonian, forming basin fill units; (5) formation of a smooth crater fill unit during the Early to Late Hesperian; (6) deposition during the Late Hesperian to Early Amazonian, creating light-toned layered units within Ladon Valles, Ladon basin, and other smaller valley networks along the western uplands, as well as formation of a light-toned unit along some crater floors; (7) formation a of chaotic unit in the Early to Middle Amazonian, formation of alluvial fans in the Early Amazonian, and eruption of a volcanic unit in the Middle Amazonian; and (8) formation of craters throughout the geologic history of the map region.

Scientific Investigations Map↗
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗
USGS2026

Legacy of the fumigant 1,2-dibromo-3-chloropropane (DBCP) in California groundwater

The fumigant pesticide 1,2-dibromo-3-chloropropane (DBCP) was widely used in California agriculture during the 1960s and 1970s before being banned in 1979. Despite this ban, DBCP continues to contaminate groundwater due to its persistence and mobility. This study evaluates the distribution, historical trends, and projected persistence of DBCP in California using data from over 13,000 public supply wells and additional domestic, irrigation, and observation wells (1980-2022). Since 2010, DBCP has been detected in 9% of public supply wells statewide, with higher frequencies in the San Joaquin Valley (21%) and upper Santa Ana River watershed (13%), where DBCP use was most prevalent. Approximately 70% of wells had decreasing concentration trends, whereas increases were more common in deeper wells, indicating downward vertical migration of the DBCP front. Groundwater age estimates show that recharge timing aligns with the 1960s–1970s loading period, enabling reconstruction of peak inputs and providing a basis for age based modeling. To estimate future persistence, we applied a one dimensional advection–dispersion model that simulates long term declines in peak concentrations based on groundwater age, historical loading, and a 38 year degradation half life. Model projections suggest that concentrations above the maximum contaminant level may persist in a declining number of wells until approximately 2080 (range: 2048–2109), with longer persistence in the San Joaquin Valley. The simplified modeling framework, based on age distributions typical of wells capturing peak concentrations, can provide practical regional scale assessment of non-point source contaminants where long-term monitoring exists. This study highlights how the legacy of DBCP contamination will likely affect California's groundwater resources throughout the 21st century.

California↗
USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗
USGS2026

A three-dimensional geologic framework model of the northern Great Plains region of Montana, North Dakota, South Dakota, and Wyoming, USA

This report presents a new three-dimensional geologic framework model (GFM) of the northern Great Plains region, encompassing parts of Montana, North Dakota, South Dakota, and Wyoming. The model provides a regionally consistent, geographic information system (GIS)-ready representation of Phanerozoic sedimentary strata, major fault systems, and Precambrian basement geometry across two sedimentary basins and adjacent uplifts. More than 300,000 geologic and geophysical data inputs were synthesized to model 41 stratigraphic horizons and 47 faults, yielding an internally coherent, sealed-volume interpretation of the subsurface. The modeling workflow developed for this study demonstrates an efficient and scalable approach for constructing basin-to regional-scale GFMs in geologically complex and data-variable settings. Although model fidelity varies with data density and quality, the resulting geometry is broadly consistent with 1:500,000-scale geologic mapping and highlights areas where additional geologic study is most needed. The three-dimensional GFM provides a foundational framework to support groundwater, energy, and mineral resource assessments, and offers a transferable methodology for potential future U.S. Geological Survey efforts to build large-area subsurface models in underexplored regions of the United States.

Montana, North Dakota, South Dakota, Wyoming↗
USGS2024

Geomorphic map of the Umatilla River corridor, Oregon

This map portrays the distribution of landforms along the Umatilla River in northeastern Oregon and covers a corridor 127 kilometers long from the confluence of the Umatilla River with the Columbia River upstream to Meacham Creek. The map encompasses the valley bottom and extends about 1 kilometer up the adjoining hillslopes. Map data are intended to support water quality and fisheries enhancement efforts pursuant to the First Foods, a resource-management approach that focuses on traditionally gathered foods including water, fish, big game, roots, and berries and calls attention to the reciprocity between people and the foods upon which humans depend. The Umatilla River drains about 6,300 square kilometers on the northwest slope of the Blue Mountains in northeast Oregon. Most of the drainage basin is underlain by Miocene basalt flows of the Columbia River Basalt Group. Younger, weakly lithified, late Miocene and early Pliocene gravel deposits of local origin (for example, McKay Formation) are mapped in a few places. Upland surfaces are mantled with windborne silt (loess) correlative with deposits elsewhere known as the Palouse Formation. Surfaces below an elevation of about 340 meters were inundated repeatedly by large Pleistocene glacial outburst floods, most emanating from glacial Lake Missoula in western Montana. In backflooded areas such as the lower Umatilla River valley, Missoula floods deposited extensive slack-water silt. Areas mapped as open water, active channel and tie channel, flood basin, valley bottom, and modified land constitute the geomorphic floodplain: the area subject to occasional inundation by the Umatilla River. Deposits and landforms within the floodplain are inset into Missoula flood deposits and hence postdate the 20–15-kilo-annum Missoula floods. Some floodplain deposits are no more than a few centuries old, as indicated by substantial erosion and deposition during the Umatilla River flood of February 2020, the largest since systematic measurements began in October 1903. Deposits and landforms of the floodplain are transient features within the longer-term incision of the Umatilla River into mid-Miocene flood basalts and younger gravel of the McKay Formation.

Oregon↗
USGS2025

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗