USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

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 52Linked to original sources
USGS2025

Characterization of stream water quality and groundwater levels in the Central Pine Barrens region, Suffolk County, New York, 2017–23

The area locally known as the “Central Pine Barrens” region, located in Suffolk County, New York, contains most of Long Island’s preserved and undeveloped land. This region overlays an aquifer system that provides potable groundwater for residents of Suffolk County. Between 2017 and 2023, the U.S. Geological Survey, in cooperation with the Central Pine Barrens Joint Planning & Policy Commission and the Town of Brookhaven, monitored groundwater levels and stream water quality in this region. Groundwater levels were measured monthly at five wells and continuously (15-minute intervals) at a sixth well. Water quality was monitored at five locations in the Carmans River and at two locations in the Peconic River, and samples were analyzed for major ions, trace elements, nutrients, pharmaceuticals, and pesticides. The major ion compositions at the sites were mainly sodium-chloride type waters, and compositions varied the most at the furthest upstream sites in both streams. Concentrations above aquatic-life criteria thresholds also occurred most frequently at the furthest upstream sites. The seasonal patterns of nutrient loads and concentrations varied between the Carmans and Peconic Rivers. Several organic compounds including pharmaceuticals, domestic use products, and pesticides were detected at low concentrations in both streams. Metformin was the most frequently detected pharmaceutical compound, and herbicides were the most frequently detected pesticide class. Water-quality conditions influenced by anthropogenic contributions are a result of current and historical land use, and these contributions include onsite wastewater disposal systems, commercial or domestic fertilizers and pesticides, and urban or industrial contaminants in road runoff. This study characterizes and improves understanding of the current hydrologic conditions in the Central Pine Barrens region and the study findings can help inform the development of plans to manage, protect, and restore water resources.

New York↗
USGS2024

Using geologic mapping to understand temporal and spatial relations of closely clustered to concurrent latest Holocene surface ruptures on two intersecting faults, south-central Mojave Desert, California

The Pinto Mountain Fault Zone (PMFZ) marks a major structural boundary between east-oriented sinistral faults of the eastern Transverse Ranges (to the south) and northwest-oriented dextral faults of the south-central Mojave Desert (to the north). These structural fault systems comprise sinistral and dextral deformational domains of the Eastern California Shear Zone (ECSZ) that intersect one another in the Copper Mountain and Twentynine Palms areas. The U.S. Geological Survey (USGS) is conducting detailed geologic mapping and geochronologic investigations designed to clarify geometric, kinematic, and temporal relations among the two domains, that are focused on the central portion of the left-lateral PMFZ near its intersection with the right-lateral Copper Mountain Fault (CMF) and Mesquite Lake Fault Zone (MLFZ).

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

Relationship between peak and cumulative ground motions from 49 Mw 3-6 earthquakes in the San Francisco Bay Area

We examine the relationship between peak ground velocity (PGV) and cumulative absolute displacement (CAD) for 49 M w 3 – 6 earthquakes in the San Francisco Bay Area (SFBA) and gain insight into the spatiotemporal partitioning of seismic energy in ground motion records with respect to source, path, and site effects. PGV and CAD are positively correlated, but there can be large deviations from the average trend. For example, ground motion records with either very long duration resonance or short pulse-like motions will have higher or lower CAD, respectively, but could have very similar PGV. We perform principal component analysis (PCA) on PGV-CAD for >7000 records in the SFBA with the goal of investigating what influences positive or negative anomalies in cumulative motions. PCA rotates the PGV-CAD datapoints into two principal components, where the one with the larger variance, which we call the “primary intensity” represents mostly the distance-dependence of ground motion amplitudes. The other principal component, which we call the “excess motion”, represents the deviation from cumulative motions that would be typical for a ground motion record with a given PGV. The excess motion will be positive in the case of records with long duration ringing and will be negative for short duration pulse-like ground motions. We find that excess motion is generally positive at sites in sedimentary basins and in soft sediments around the SF Bay. Excess motion is generally negative in the very near field, as well as at sites on hard bedrock. We discuss the findings here in terms of implications for seismic hazard applications and other wave propagation phenomena.

California↗
USGS2023

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗
USGS1977

Leachate plumes in ground water from Babylon and Islip landfills, Long Island, New York

Landfills operated by the Towns of Babylon and Islip in southwest and central Suffolk County contain urban refuse, incinerated garbage, and scavenger (cesspool) waste; some industrial refuse is deposited at the Babylon site. The Islip landfill was started in 1933, the Babylon landfill in 1947. The landfills are in contact with and discharge leachate into the highly permeable upper glacial aquifer (hydraulic conducitvity 190 to 500 ft/d [feet per day]). The aquifer is 74 feet thick at the Babylon landfill and 170 feet thick at the Islip landfill. The leachate-enriched water occupies the entire thickness of the aquifer beneath both landfills, but hydrologic boundaries retard downward migration of the plumes to deeper aquifers. The Babylon plume is 1,900 feet wide at the landfill and narrows to about 700 feet near its terminus 10,000 feet from the landfill. The Islip plume is 1,400 feet wide at the landfill and narrows to 500 feet near its terminus 5,000 feet from the landfill. Hydrochemical maps and sections show the distribution of the major chemical constituents of the plumes. The most highly leachate-enriched ground water obtained was from the Babylon site; it contained 860 mg/L (milligrams per liter) sodium, 110 mg/L potassium, 565 mg/L calcium, 100 mg/L magnesium, 2,700 mg/L bicarbonate, and 1,300 mg/L chloride. Sulfate was notably absent or in low concentration in most parts of both plumes. Nitrogen in plume water was mostly in the form NH4+, and concentrations as high as 90 mg/L were found; concentrations of nitrogen as N in the plume were less than 10 mg/L. As much as 440 mg/L iron and 190 mg/L manganese were found in the leachate-enriched water. Samples were also tested for arsenic, boron, cadmium, cobalt, chromium, copper, mercury, nickel, lead, selenium, strontium, and zinc. Boron was more or less ubiquitous and was found in concentrations as high as 2 mg/L. Organic carbon was found in concentrations as high as 2,250 mg/L in the most highly leachate-enriched water but attenuated rapidly to less than 20 mg/L. Dissolved-solids concentrations near the landfills were between 400 and 3,000 mg/L at Babylon and between 500 and 1,500 mg/L at Islip. Ground-water temperatures near the landfills exceed those in ambient water by as much as 7°C at Babylon and 16°C at Islip. Heat contributed by the landfills was mostly dissipated within 0.4 mi of the landfill, butat Islip, the warm leachate-enriched water extended 0.5 mi downgradient. The entrance of leachate into the less dense ground water as pulsations after rainfall may explain the presence of high leachate enrichment near the bottom of the aquifer. A comparison of the physical characteristics of leachate-enriched ground water with those of ambient water suggests that the downward movement of leachate results from its greater density. Simulation of the movement and dispersion of the Babylon plume with a mathematical dispersion model indicated the coefficient of longitudinal dispersion to be about 60 ft 2 /d (feet squared per day) and the ground-water velocity to be 1 ft/d. However, the velocity determined from the hydraulic gradient and public-supply wells in the area was 4 ft/d, which would cause a plume four times as long as that predicted by the mathmatical dispersion model. At the Islip site, the plume was one-third the length calculated on the basis of the age of the landfill. The shortness of the plumes has not been explained; it may be a result of the leachate's having been too dilute to form a plume during the early years of the landfills.

New York↗
USGS2025

U.S. Geological Survey Integrated Water Availability Assessment—2010–20

This professional paper is a multichapter report that assesses water availability in the United States for water years 2010–20. This work was conducted as part of the fulfillment of the mandates of Subtitle F of the Omnibus Public Land Management Act of 2009 (Public Law 111-11), also known as the SECURE Water Act. As such, this work examines the spatial and temporal distribution of water quantity and quality in surface water and groundwater, as related to human and ecosystem needs and as affected by human and natural influences. Chapter A introduces the National Integrated Water Availability Assessment and provides important background and definitions for how the report characterizes water availability and its components. Chapter A also presents the key findings of Chapters B–F and thus acts as a summary of the entire report. Chapter B is a national assessment of water supply, which is the quantity of water supplied through climatic inputs. Chapter C is a national assessment of water quality, which is the chemical and physical characteristics of water. Chapter D assesses water use including withdrawals and consumptive use in the conterminous United States. Chapter E presents an analysis of factors affecting future water availability under changing climate conditions. The National Integrated Water Availability Assessment culminates with Chapter F, which is an integrated assessment of water availability that considers the amount and quality of water coupled with the suitability of that water for specific uses. Together, these six chapters constitute the National Integrated Water Availability Assessment for water years 2010–20.

Professional Paper↗
USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS2026

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater↗
USGS2026

Simulation of groundwater flow in the Silurian aquifer, eastern Iowa, 2020–45

The Silurian aquifer is an important water source for municipalities, industry, and rural households and communities in eastern Iowa, including Johnson County. Increasing demand for groundwater from the aquifer indicated that analytical tools may be beneficial to quantify groundwater resources and inform water-management decisions for the aquifer. The U.S. Geological Survey, in cooperation with the Johnson County Board of Supervisors, developed conceptual and numerical groundwater models to simulate water levels in the Silurian aquifer, determine groundwater budgets, and forecast changes in groundwater levels through 2045. A MODFLOW numerical model was constructed to match water levels for 2020 through 2022. Model performance was assessed by comparing simulated and observed water-level hydrographs at several pumping and monitoring wells. Hydrograph comparison indicates that the model simulates water levels accurately at monitoring wells that are not near areas of pumping and generally overpredicts water levels near pumping wells, while accurately simulating the drawdown trend over time at pumping wells. This transient model framework was then used to predict water levels through 2045 in a scenario of moderate drought and increased groundwater withdrawals. Results from the predictive model indicate as much as 13 meters of additional drawdown in 2045 from 2020 water levels in the area of pumping wells.

Iowa↗
USGS2025

The impact of burial diagenesis on soil-formed minerals in paleosols using stable isotopes of phyllosilicates and carbonate clumped isotopes

To understand the effects of burial diagenesis on the stable isotope geochemistry of soil-formed clay and carbonate minerals in paleosols, samples were collected from seven cores, spanning middle- to upper-Pennsylvanian strata of the Illinois Basin, with varied maximum burial depths of 1–3 km. Mixed-layer illite-smectite and kaolinite mixtures give δ 2 H and δ 18 O values of −83 ‰ to −36 ‰ and 11.9 ‰ to 21.1 ‰ (VSMOW), respectively. After carbonates were screened petrographically for diagenetic textures using transmitted light and cathodoluminescence, measured clumped isotope Δ 47 values range from 0.504 to 0.563 ‰ (I-CDES). Resulting mineral formation temperatures for phyllosilicate mineral mixtures are 28 to 66 °C (mean = 47 °C), whereas T(Δ 47 ) estimates for calcites are 36 to 61 °C (mean = 45 °C). Calculated δ 18 O water values from which phyllosilicate minerals and calcites precipitated under isotopic equilibrium ranges from −7.1 to −1.2 ‰ and − 1.4 to +4.9 ‰, respectively. Closed and open-system phyllosilicate-fluid exchange modeling indicates that phyllosilicate alteration occurred in the presence of a low temperature brine or meteoric water and is interpreted to occur in a layer-by-layer illitization transformation. Due to the lack of diagenetic textures and positively correlated T(Δ 47 ) and δ 18 O water , calcites are interpreted to have undergone solid-state bond reordering. Despite low to moderate temperatures (<125 °C) and varying depths of shallow burial (1–3 km), solid-state transformation of phyllosilicates and calcites indicates paleosols had prolonged exposure to burial conditions which has implications for the use of paleosol minerals for paleoenvironmental reconstructions.

Illinois, Indiana, Kentucky↗
USGS1998

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida↗
USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗
USGS2025

Quantitative mineral resource assessment of lithium pegmatite deposits in the Appalachian Orogen, USA

Lithium is classified as a U.S. critical mineral commodity, and its demand is projected to drastically increase through 2040, driven by electric vehicle production and energy storage applications (IEA 2021).Most global lithium production is not in the United States increasing vulnerability to a supply disruption. The U.S. Geological Survey is actively assessing domestic lithium deposits including lithium-bearing pegmatites in the Appalachian orogen. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical, and mineral occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. Estimates were then integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was used to estimate the amount of potentially recoverable undiscovered resources. Preliminary computations for the northern Appalachians, including application of the economic filter to the median recoverable contained resource, yields 900,000 metric tons of Li 2 O that correspond to enough Li 2 O to replace 127 years of import reliance at the current rate (7,100 t Li 2 O/yr; USGS, 2025). For the southern Appalachians, preliminary computations yielded 1,430,000 metric tons of Li 2 O, which corresponds to 201 years of import reliance.

Alabama, Connecticut, Delaware, Georgia, Maine, Ma↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2024

Perchlorate, metals, organic compounds, and lead isotopes in groundwater, surface water, shallow groundwater, and soil within and near the Middleton Municipal Airport–Morey Field (C29), Middleton, Wisconsin, 2022

The Middleton Municipal Airport–Morey Field (C29) is in the City of Middleton and adjacent to the towns of Middleton and Springfield, Wisconsin. Nearby homes in the towns rely on private drinking water wells, and residents are concerned about the potential contamination of groundwater and surface water by airport activities, including flights by small aircraft that use leaded aviation fuel and a fireworks display in July 2021. The U.S. Geological Survey, in cooperation with the Town of Middleton, completed a study in 2022 to characterize the occurrence and sources of perchlorate, metals (including lead), and organic compounds in samples of groundwater, surface water, shallow groundwater, and soils within and near the airport. Lead isotopes were also measured to determine sources of lead by comparing samples to environmental references. Magnitudes of concentrations from samples of water and soil collected in 2022, and their spatial patterns across site locations, indicate the fireworks display in July 2021 was a likely source of perchlorate and metals in the airport study area. The highest perchlorate concentration was measured in surface water at the southeastern corner of the airport near the fireworks launch site; the highest concentrations of fireworks-associated metals were measured in shallow groundwater near the same location. Fireworks were not the only possible source of perchlorate and metals in the airport study area because both were also detected upgradient and away from the fireworks launch site. Ratios of lead isotopes indicate that lead measured in water and soil within the airport study area was primarily sourced from background atmospheric lead deposition or Wisconsin galena lead ore. However, two groundwater samples (one upgradient and one downgradient from the airport; both with concentrations less than 1 microgram per liter) had isotopic signatures matching leaded aviation fuel sold at the airport.

Wisconsin↗
USGS2025

Central Valley Hydrologic Model version 2 (CVHM2): Decision support tool for groundwater and land subsidence management

The San Joaquin Valley (SJV) of California is one of the world’s most productive agricultural regions. Reliance on groundwater has led to some of the greatest rates of human-induced land subsidence in the world in the 20th century, as well as more recently. The United States Geological Survey (USGS) has recently developed an integrated surface–subsurface hydrologic model, the Central Valley Hydrologic Model 2 (CVHM2), that represents the major components of the hydrologic system of California’s Central Valley. In this study, CVHM2 was applied as a decision support tool while simulating various management strategies to mitigate the land subsidence caused by the extraction of groundwater. CVHM2 was extended through to 2073 and applied to simulate management scenarios in terms of three primary drivers and their impact on subsidence along the Delta–Mendota Canal (DMC), a critical piece of infrastructure in the western SJV. The drivers considered were agricultural water demands, managed aquifer recharge (MAR), and changes in future climate. The results show that future subsidence is most sensitive to water demands, second most sensitive to future changes in climate, and relatively insensitive to MAR when it is applied as a surface application in the western SJV. However, we demonstrate via proof-of-concept scenarios that the MAR is capable of arresting subsidence when implemented via injection below the Corcoran Clay Member of the Tulare Formation instead of as a surface application. We also examine the uncertainty that is the result of climate variability and how to use the tool to identify the most appropriate strategies to constrain future subsidence to acceptable levels.

California↗
USGS2026

Distributed faulting of the northern West Napa Fault Zone in Napa Valley, California

Mapped surface ruptures from the 24 August 2014 M w 6.0 South Napa earthquake in the Napa Valley, California, show a 2‐km‐wide zone of distributed faulting in the southern and central West Napa fault zone (WNFZ). In the northern WNFZ at Hendry Winery (HW), however, the mapped 2014 surface ruptures encompass an ∼100‐m‐wide zone, implying significant narrowing of the near‐surface fault zone to the north. We present a tomographic shear‐wave velocity ( ⁠⁠ V S ) model and guided‐wave data that indicate the northern WNFZ is at least 400‐m wide, with multiple near‐surface fault traces. Our V S model shows that the 2014 surface ruptures are underlain by discrete low‐velocity zones (LVZs), and coincident guided‐wave data show that the LVZs carry fault‐zone guided waves. If nearby (<500 m) mapped faults to the east of HW are part of the WNFZ, the entire WNFZ is more than 1 km wide in the northern Napa Valley. WNFZ guided waves travel up to 38% slower than S body waves, and low‐strain guided‐wave shaking is up to five times stronger than the associated body‐wave shaking. Our data suggest that guided waves, traveling along distributed faults, may result in an increased shaking hazard over a 1‐km‐wide area of the northern Napa Valley during future significant earthquakes. In places, the 2014 surface ruptures were difficult to find one year after the earthquake, and paleoseismic trenching showed only weak evidence for faulting, which may not have been identified in trenches if the locations of the 2014 surface ruptures had not been previously mapped ( Prentice et al. , 2015 ). Guided‐wave and V S tomography data, however, show strong evidence for faulting beneath the 2014 surface ruptures and at locations to the east. Although paleoseismic trenching is the gold standard for identifying near‐surface faulting, methods such as peak ground velocities of guided waves may better identify immature near‐surface fault traces.

California↗