Geology Reports⌕ Search

SEARCH · Geology Reports

Search Geology Reports

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,748 records · Page 51Linked to original sources
USGS2025

Factors regulating the potential for freshwater mineral soil wetlands to function as natural climate solutions

There are increasing global efforts and initiatives aiming to tackle climate change and mitigate its impacts via natural climate solutions (NCS). Wetlands have been considered effective NCS given their capacity to sequester and retain atmospheric carbon dioxide (CO 2 ) while also providing a myriad of other ecosystem functions that can assist in mitigating the impacts of climate change. However, wetlands have a dual impact on climate, influencing the atmospheric concentrations of both CO 2 and methane (CH 4 ). The cooling effect associated with wetland CO 2 sequestration can be counterbalanced by the warming effect caused by CH 4 emissions from wetlands. The relative ability of wetlands to sequester CO 2 versus emit CH 4 is dependent on a suite of interacting physical, chemical, and biological factors, making it difficult to determine if/which wetlands are considered important NCS. The fact that wetlands are embedded in landscapes with surface and subsurface hydrological connections to other wetlands (i.e., wetlandscapes) that flow over and through geochemically active soils and sediments adds a new layer of complexity and poses further challenges to understanding wetland carbon sequestration and greenhouse gas fluxes at large spatial scales. Our review demonstrates how additional scientific advances are required to understand the driving mechanisms associated with wetland carbon cycling under different environmental conditions. It is vital to understand wetland functionality at both wetland and wetlandscape scales to effectively implement wetlands as NCS to maximize ecological, social, and economic benefits.

Wetlands↗
USGS2025

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska↗
USGS2025

Generalized Bancroft algorithm for locating earthquakes with P- and S-wave arrival times

Because of similarities between locating an earthquake with seismic stations and locating a Global Positioning System (GPS) receiver from satellites, the Bancroft algorithm developed for GPS processing can be used to locate earthquakes. Such an approach to earthquake location differs from the conventional method of choosing an initial or trial solution and then iteratively improving the solution until convergence. The Bancroft algorithm has the advantage of being a direct, noniterative solution but with the disadvantage of only being able to accommodate a homogeneous velocity model. An additional limitation of the standard Bancroft algorithm is that it considers arrival times in a medium with a single propagation velocity. This poses no problem for GPS processing because electromagnetic waves travel at the speed of light; however, for seismic waves it means the algorithm can be applied to collections of either P ‐ or S ‐wave arrival times. Here, I show how the Bancroft algorithm can be generalized to handle both P ‐ and S ‐wave arrival‐time measurements simultaneously. I also show how to accommodate depth‐varying P ‐ and S ‐wave velocity models. I apply the generalized Bancroft algorithm to microearthquakes beneath Tanaga Volcano in Alaska and compare standard locations from the widely used HYPOINVERSE location code to Bancroft locations and to the output of HYPOINVERSE when setting the trial location to the Bancroft location. I find the Bancroft locations outperform the results from the other methods for shallow earthquakes near sea level, where a quantity known as the geometric dilution of precision is large and linearized approaches such as HYPOINVERSE are expected to struggle.

Alaska↗
USGS2025

Map of topographic lineaments interpreted as recent surface ruptures along the Bennett Valley and Southern Maacama Fault Zones, Sonoma County, California

This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.

California↗
USGS2025

Groundwater-storage change in the north Phoenix aquifer, Arizona, 2020–23

The city of Phoenix, Arizona, relies primarily on surface water for municipal water supply. The city also maintains wells to withdraw groundwater, particularly in times of drought and reduced surface-water supply, and to recharge groundwater when excess surface water is available. As of 2023, withdrawals from the aquifer in the northeastern part of the city are a small volume of water, less than 3,000 acre-feet in most years. Each year a similar volume of water is recharged through injection wells. The withdrawal wells are permitted to produce in total more than 20,000 acre-feet per year; increased pumping could lead to future groundwater-storage declines. To better understand groundwater-storage change in the north Phoenix aquifer, a repeat microgravity monitoring network was established in 2020. Measurements of changes in Earth’s gravity provide a direct, non-invasive measurement of subsurface mass change. Groundwater-storage changes were small during the 2020–23 study period, consistent with the relatively small volumes of pumping and recharge. Groundwater levels measured in monitoring wells were stable or increased slightly during this period, although the number of monitoring wells within the area of the gravity monitoring network is sparse. In total, about 15,000 acre-feet of water were pumped and 31,000 acre-feet recharged through injection in the north Phoenix aquifer during the 2020–22 period. Within the monitored area, groundwater storage increased by a small amount, about 1,000 acre-feet, in 2020, and decreased by a larger amount, about 6,000 acre-feet, each year in 2021 and 2022. Storage decreased at 89 of 102 stations from 2020 to 2023. Groundwater-storage decreases greater than the volume of net pumping indicate down-gradient subflow from the aquifer is greater than recharge plus incoming subflow, drying of the unsaturated zone resulting from decreased land-surface recharge, or both. At present (2023), the aquifer appears able to store and supply the relatively small amounts of water needed without excessive drawdown or harmful effects, such as land subsidence and surface inundation. If pumping or recharge increases in the future, the established repeat microgravity network is well suited to capture the spatial extent and magnitude of groundwater-storage changes in the aquifer.

Arizona↗
USGS2025

Metal-rich lacustrine sediments from legacy mining perpetuate copper exposure to aquatic-riparian food webs

Historic copper mining left a legacy of metal-rich tailings resulting in ecological impacts along and within Torch Lake, an area of concern in the Keweenaw Peninsula, Michigan, USA. Given the toxicity of copper to invertebrates, this study assessed the influence of this legacy on present day nearshore aquatic and terrestrial ecosystems. We measured the metal (Co, Cu, Ni, Zn, Cd) and metalloid (As) concentrations in sediment, pore water, surface water, larval and adult insects, and two riparian spider taxa collected from Torch Lake and a nearby reference lake. Overall, elevated metal and metalloid concentrations, particularly Cu, were measured in all sediment samples and some surface and pore water samples collected from Torch Lake. For instance, Cu concentrations in the Torch Lake sediment were ∼200% higher than the reference lake and all measured concentrations exceeded predicted effects concentrations by at least ninefold. Within larval insect tissues, we observed 160% higher Cu concentrations than measured in the reference lake, and Cu was the only measured element above predicted effects concentrations in Torch Lake. Adult insects collected at both lakes had similar metal concentrations irrespective of exposure levels. Yet we found 100% higher copper concentrations in Torch Lake riparian spiders, demonstrating elevated exposure risk to insectivores across the aquatic-terrestrial boundary. Our results highlight that other metals in the mixture may not be as concerning to adjacent riparian ecosystems, but copper remains a contaminant of concern in Torch Lake 60 years after mining ceased.

Michigan↗
USGS2025

Woods Hole Coastal and Marine Science Center—2023 annual report

The 2023 annual report of the U.S. Geological Survey Woods Hole Coastal and Marine Science Center highlights accomplishments of 2023, includes a list of 2023 publications, and summarizes the work of the center, as well as the work of each of its science groups. This product allows readers to gain a general understanding of the focus areas of the center’s scientific research and learn more about specific projects and progress made throughout 2023, all while enjoying photographs taken in various environments and laboratories, and applicable maps and figures.

Circular↗
USGS2025

Evaluation of solid bitumen created from marine oil shale bituminite under hydrous and anhydrous pyrolysis conditions

To test the influence of environmental conditions on aromaticity of solid bitumen generated during petroleum generation, four organic-rich (26–36 wt% total organic carbon) oil shale samples collected from the Neoproterozoic–Lower Cambrian restricted marine Salt Range Formation in the upper Indus Basin, Pakistan, were pyrolyzed under hydrous and anhydrous conditions. Experiments used closed system batch reactors at subcritical water temperatures between 300 and 370°C for 72 h. Thermal conversion of bituminite in the Salt Range oil shales created a newly formed solid bitumen, similar to previous observations from experiments on the Eocene lacustrine Green River Formation Mahogany zone oil shale. Newly formed solid bitumen in the Salt Range Formation oil shales generally has higher reflectance (R o ) in hydrous residues compared to anhydrous experiments subjected to the same time-temperature conditions, also similar to prior observations. This finding supports the hypothesis that radical disproportionation is favored in the presence of hydrogen contributed by water, promoting aromatization in the solid bitumen residue with concomitant increased expulsion of saturated hydrocarbons. Indigenous solid bitumen (and vitrinite in a comparison sample) also shows higher reflectance values in hydrous versus anhydrous residues, indicating that additional aromatization in the presence of exogenous hydrogen occurs both in newly formed organic matter and in organic matter that is present throughout the experiment. Despite similarities in their bulk rock geochemical screening parameters, R o evolution shows different trajectories amongst the four Salt Range oil shales, suggesting as-yet undetermined differences in kinetic properties which are probably related to differences in a priori chemical composition. These results have implications for the use of solid bitumen reflectance (BR o ) as a thermal proxy, suggesting BR o values and appearance could vary as a function of the concentration of water. Variation in water concentration may be present at the reservoir or formation scale, but may also be present at a much finer scale in tight oil shales where permeability is several orders of magnitude lower than conventional reservoirs. Therefore, local variations in the presence of water potentially could explain substantial variation in BR o values and appearance in closely spaced source rock samples and even within an individual microscope field.

Journal of Analytical and Applied Pyrolysis↗
USGS2025

Smectite-rich horizons in Inceptisols trigger shallow landslides in tropical granitic terranes

Puerto Rico was affected by >70,000 landslides in the wake of 2017 Hurricane Maria, and landslide prevalence was especially high in the Utuado region in the Cordillera Central. Landslide density was highest where soil parent material is granodiorite; landslide slip surfaces tended to be shallow (<60 cm), and often were mobilized rapidly and with long runout distances. This study combines field observations with soil mineralogy (bulk and clay fractions), soil geochemistry (bulk fraction), and soil strength as determined by field cone penetrometer testing (CPT) to test the hypothesis that clay-rich subsoil horizons function as slip planes when water-logged. Soil pits were excavated to depths of ∼200 cm in Ultisols on an undulating plateau and to ∼100 cm in Inceptisols on steep slopes (36-43 o ) that flank the plateau and cone penetrometer tests (CPT) were done within 2 m of the soil pit. Six pits were located adjacent to scarps from previous landslides, enabling analysis of soil profiles downward through extrapolated slip surfaces. Results from X-ray diffraction (XRD), X-ray fluorescence (XRF) and thermogravimetric analysis (TGA) indicate that soils are heterogeneous, often with subsoil horizons enriched in clay minerals and immobile elements (Al, Fe, Ti). Inceptisols on steep slopes often contain smectite-rich horizons at 30–60 cm depth that appear to function as slip surfaces; in other Inceptisols, such horizons are not present and landslide susceptibility is potentially lower. In Ultisols, soil mineralogy is dominated by kaolinite with minor halloysite, and increased kaolinite content at soil depths ≥80 cm at some sites suggests potential slip surfaces enhancing probability of landslides. The origin of clay-rich horizons appears to be (1) fractures in granodiorite that facilitate water flow and leaching, accelerating mineral dissolution during early weathering stages, and (2) smectite-rich buried soils under permeable colluvium likely deposited by a prior mass wasting event. Where clay-rich layers occur beneath more-permeable horizons, rapid infiltration then absorption of water in clay-rich subsoil horizons causes decreased shear strength and increased landslide susceptibility.

Puerto Rico↗
USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗
USGS2025

U.S. Geological Survey research and assessments supporting carbon dioxide removal

Both carbon capture and storage (CCS) and carbon dioxide removal (CDR) are methods to limit future global temperature rise and ocean acidification. CCS sequesters (stores) carbon dioxide (CO2) captured from industrial sources thereby preventing the CO2 from reaching the atmosphere. The captured CO2 is injected into underground geologic reservoirs or is converted into stable industrial products. CDR enhances natural processes or uses technical approaches to remove excess CO2 from the atmosphere. Carbon dioxide is sequestered in natural sinks such as vegetation and soil, or injected into underground reservoirs or converted into durable products. The U.S. Geological Survey (USGS) conducts research and resource assessments that support biologic and geologic CDR approaches. This report will review some of the ongoing USGS CDR research and assessment efforts.

Conference Paper↗
USGS2025

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska↗
USGS2025

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗
USGS2025

Submarine avalanche deposits hold clues to past earthquakes

Earthquakes and other natural events sometimes shake the seafloor near coastlines severely enough to cause underwater avalanches that rush down steep slopes, scouring the seabed and carrying sediment to greater depths. These fast-moving sediment-laden flows, called turbidity currents , have at times damaged underwater infrastructure like pipelines and communications cables, as they did, for example, in snapping transatlantic cables off the coast of Newfoundland after the 1929 Grand Banks earthquake.

EOS Transactions↗
USGS2025

Using hydrogeologic context and water budgets to evaluate the potential for groundwater contributions to contaminants in Lake Superior

This study presents a synthesis of the hydrogeology in the U.S. Lake Superior watershed and the contribution of groundwater to the water budget of the U.S. Lake Superior basin. Much of the shoreline of Lake Superior in Minnesota and Michigan is composed of hydrogeologic units contributing very little direct groundwater discharge to the lake. Groundwater in watersheds adjacent to the lake typically flows in short, local flow systems characterized by thin glacial sediments with active groundwater flow in fractured bedrock within the top 60–90 m below land surface. The complex groundwater system in Wisconsin’s Bayfield Peninsula has the largest groundwater reservoir near the lake, characterized by thick sand and gravel glacial deposits and underlying sandstone aquifer. Although these thick sandy glacial deposits are not in direct contact with the lake at the shoreline, groundwater discharge may still be significant via subsurface exposures beyond the shoreline or flow through the underlying sandstone aquifer. Overall, most groundwater in the watershed is contributed as indirect base flow in streams around the lakeshore and comprises about 60 percent of the upland water budget. Direct groundwater flow to the shoreline contributes 2 to 9 percent of the inflow. Identifying possible contamination sources through direct sampling of groundwater would be an inefficient way to detect problems if sources are unknown, particularly for some chemicals of concern such as PFAS, pesticides, PCBs, chloride, and nutrients. Evaluating the chemical characteristics of contamination is also important to consider in evaluating how groundwater may contribute to pollution in Lake Superior.

Lake Superior↗
USGS2025

Quaternary-active faults and the role of inherited structures in the Sacramento-San Joaquin Delta, western Central Valley, northern California

Seismic sources and their associated hazards within the Sacramento-San Joaquin Delta region of north-central California are relatively poorly characterized as compared to other, more heavily studied regions of northern California, such as the San Francisco Bay Area. Here we present a synthesis of subsurface, bedrock geology, and geodetic datasets from the Delta and from the Coast Ranges and Diablo Range to the northwest and southwest, respectively. We integrate these data and our own surface geologic and geomorphic observations to present a comprehensive review of faults in the Delta that exhibit Quaternary activity. Structural geologic data from the surrounding region highlight the significant influence that Late Cretaceous-to-Paleogene forearc structures exert on the geometry and kinematics of major Quaternary-active structures within the Delta. These inherited structures — including the Pittsburg-Kirby Hills Fault, Midland Fault, and Great Valley Fault System — exhibit a range of geometries and kinematics. Analysis of geomorphology along these structures suggests that these structures combine to accommodate Quaternary strain across the Delta region. A clearer understanding of subsurface geometries and structural relationships, built upon the regional tectonic history, provides insight into modern deformation accommodated on older structures and helps inform interpretations of seismic hazard within the Delta.

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

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗