USGS ScienceSearch

SEARCH · USGS Science

Results for “The Leading Edge”

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,663 records · Page 29Linked to original sources

Evaluating the potential of co-located non-petroleum energy resources in the Gulf Coast using play fairwayaAnalysis

Geological resources critical to the energy transition, such as sedimentary geothermal, carbon storage potential, and lithium in brines, commonly struggle for economic feasibility as stand-alone developments but can have greater viability when the potential for more than one of these resources exist within the same reservoir or stacked in different stratigraphic intervals. There may also be instances where development of one resource inhibits development of others and decisions need to be made on how to best prioritize the use of those resources in the subsurface. Subsurface data sets were analyzed and integrated to evaluate the distribution of non-petroleum energy and related resources in the onshore and nearshore U.S. Gulf Coast. Temperature, pressure, brine composition (lithium content), and reservoir quality data for thirty-four depositional units have been compiled and visualized to high-grade areas where multiple resource opportunities likely coexist. For sedimentary geothermal, possible resource areas are defined as low potential (<90°C), moderate potential (90–150°C), and high potential (>150°C). For CO 2 storage, high potential areas exhibit supercritical CO 2 conditions less than 80% of the fracture gradient. Areas with pore pressure between 80% and 100% of the fracture gradient are considered to exhibit moderate potential and areas where the fracture gradient is equal to or greater than 100% are assigned low potential. Lithium resources in brines were defined by lithium concentrations as low potential (<100ppm), moderate potential (100-200ppm), or high potential (>200ppm). Reservoir quality affects the viability of all three of these resources and is evaluated using depositional environment maps of each unit. The resulting play fairway maps can be used for regional scale screening evaluations of these resources and to identify areas of interest where more detailed, prospect-scale studies can be undertaken.

Gulf Coast

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

Regional variations in sea ice and primary productivity in the Bering Sea during Marine Isotope Stage 11

Marine Isotope Stage (MIS) 11 (424-374 ka) has long been an analog for Holocene climate, because it is the most recent interglacial period with similar orbital conditions. However, there is significant global and regional variability in the climate response to MIS 11 warmth. Here, we review sediment core records from across the Bering Sea to investigate changes in paleoceanographic conditions during Marine Isotope Stages 12-10. Sea ice was present over much of the Bering Sea during MIS 11, but today, none of the sites investigated are ever ice-covered. This suggests that sea ice regimes in the Bering Sea during MIS 11 were different to those of the Holocene. There are also regional differences in the response of sea ice to MIS 11 warming. At the Umnak Plateau, Southeastern Bering Sea, sea ice concentrations decline during deglaciation, but they remain high at the slope sites until Peak MIS 11. Sea ice re-advances over the Umnak Plateau during peak interglacial warmth, at the same time that it declines over the slope sites. Late MIS 11 is characterized by high concentrations of seasonal sea ice at the Umnak Plateau, whilst sea ice at the slope sites fluctuates between consolidated and unconsolidated ice cover. This east-west dichotomy may be explained by changes in the behavior of the Aleutian Low. Productivity increases dramatically during deglaciation due to increased upwelling and sea level rise bringing fresh nutrients into the oceans. This is characterized by increased diatom productivity, increased terrestrial carbon deposition, and laminations at all sites.

Alaska

California State Waters Map Series—Benthic habitat characterization in the region offshore Humboldt Bay, California

Coastal and Marine Ecological Classification Standard (CMECS) geoform, substrate, and biotic component geographic information system (GIS) products were developed for the California State Waters of northern California in the region offshore of Humboldt Bay. The study was motivated by interest in development of offshore wind-energy capacity and infrastructure in Federal waters offshore. This project, carried out by the U.S. Geological Survey (USGS), resulted in four data releases for individual map blocks that are part of the “California State Waters Map Series”: (1) Offshore of Arcata, (2) Offshore of Eureka, (3) Offshore of the Eel River, and (4) Offshore of Cape Mendocino. The study area consists of 436 square kilometers of multibeam echo sounder (MBES) data acquired by Fugro Pelagos, Inc., in 2007. Towed camera-sled video was acquired in 2009 and 2010 to supervise the classification of the MBES data into habitats, and single channel sparker data were collected to calculate sediment thickness above the transgressive unconformity. Using video observations of habitat as ground truth, derivatives of the MBES data were classified into 3 seafloor character types (hard-rugose, hard-flat, and soft-flat), 26 induration-slope-depth groups, and 15 geoforms. The study area substrate is predominantly soft-flat sediment (mud and fine sand) covering 73.6 percent of the area. Hard-flat substrate areas, predominantly coarse sediment in scour depressions, cover 5.4 percent of the study area. The hard-rugose substrate areas are primarily outcrops of layered sedimentary bedrock and constitute 20.9 percent of the study area. Fifteen geoforms were identified in the analysis. The predominant geoforms mirror the seafloor character results, shelf geoforms, rock outcrop geoforms, and scour depression geoforms. Rock and scour areas are restricted to the southern portion of the study area off Cape Mendocino where uplift has exposed bedrock. On the flat shelf area post-transgressive sediment varies in thickness from 1.7 meters (m) nearshore to 28.1 m offshore.

California

Fault displacement model for surface principal rupture of strike-slip faults

The probability distribution model for principal displacement accommodated on the surface main trace is a critical input to the fault displacement hazard analysis. This article presents a new model for strike-slip ruptures in the moment magnitude ( M ) range of 6 to 8.3. The new model is the outcome of a multi-year research effort to update the widely used model developed by Petersen and others in 2011. Updates include the adoption of the Fault Displacement Hazard Initiative database and enhancements to rupture and displacement data preparation. Statistical formulation and estimation have also been updated substantially. A three-parameter modified normal distribution that we refer to as the negative Exponentially Modified Gaussian distribution is adopted to model the probability distribution of the natural logarithm of principal displacement. Formulation for the mean parameter of the modified normal includes a random earthquake term, a nonlinear scaling relation with M , and an ellipse function for along-main-trace variation. The aleatory variability of the updated model now depends on M as well as site’s along-main-trace position. These updates not only significantly improve the fit to the distribution of the observed displacements but also yield reasonable 95th percentile predictions for M > 7.5 events. Alternative models representing the estimation uncertainty of the M -scaling relation are also developed. These new models are compared to the previous model in terms of percentile predictions and the calculated hazard curves. The steeper hazard curves from the new models yield a lower exceedance rate than the normal-distribution based model developed previously by Petersen and others.

Earthquake Spectra

Magnitude and frequency of floods in the Suwannee River Water Management District, Florida

Flood-frequency statistics for 2-, 5-, 10-, 25-, 50-. 100-, 200-, and 500-year recurrence intervals, based on three methods of analysis, are presented for 25 continuous-record and seven peak flow partial-record gaging stations in the Suwannee River Water Management District. The first method, for gaged stations, utilizes station records; the second method, for ungaged sites, utilizes regional regression analysis; and the third method uses a weighted combination of the station and regional values. Because the weighted values utilize two more or less independent estimates of the peak flow statistic, they are considered more accurate than the station estimates or the regression estimates alone. Also, the use of another weighting scheme to improve estimates of flood frequency statistics at ungaged sites is demonstrated. The karstic nature of much of the Suwannee River Water Management District significantly attenuates flood peaks in some streams by providing substantial subsurface storage when river stages are high. At such times, springs discharging into rivers may reverse flow temporarily and become sinks.

Florida

An improved empirical model for predicting postfire debris-flow volume in the western United States

Reliable estimates of debris-flow volume can be used to help predict the magnitude of debris-flow hazards following wildfire in the western United States. In this study, we compiled and used a database of 227 postfire debris-flow volumes that were collected across the western United States to develop a multiple linear regression model for predicting postfire debris-flow volume. We explored 36 predictor variables related to rainfall, terrain, and fire characteristics, and selected the model with the combination of variables that yielded the most accurate predictions of debris-flow volume. We evaluated model performance against the entire volume database, as well as against four subsets of volume data from southern California, the Intermountain West, the Southwest, and regions with limited volume data, such as northern California and Washington. We also compared model performance against 3 existing postfire debris-flow volume models that were developed for use in southern California, the Intermountain West, and the Southwest. We demonstrate that the new volume model performs as well as the regional models in the regions for which they were developed and outperforms existing models when applied to volumes from data-limited regions in the western United States. These results indicate that the debris-flow volume model introduced in this study can be used to improve postfire hazard assessments across the western United States, especially outside of southern California.

Arizona, California, Colorado, New Mexico, Utah, W

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

A decade of shaking in the Garden City: The dynamics of preparedness, perceptions, and beliefs in Canterbury, New Zealand, and implications for earthquake information

This study explored earthquake preparedness over time - before, during, and 10 years after the Canterbury Earthquake Sequence (CES) in Aotearoa New Zealand (NZ). Surveys of Canterbury residents were conducted in 2009, 2013 and again in 2021, using variables derived from Community Engagement Theory (CET). The surveys measured earthquake perceptions and beliefs, participation and engagement, and preparedness actions. Results were compared across the three samples. Findings indicate that perceptions and beliefs (e.g. risk perception, outcome expectancy beliefs), and types of preparedness actions taken (e.g. collection of survival items, structural preparedness, community and agency relationships), differed over time, depending on people’s experiences before, during, and after the CES. For example, during and after the CES people were more likely to believe that preparing provided a benefit to daily life, but less likely to think it could reduce property damage, perhaps due to people’s experiences of disruption and damage during the earthquakes. An understanding of such dynamics can assist with the provision and timing of risk and preparedness information. This study highlights the importance of providing applicable and actionable preparedness information, that is relevant to people’s experiences, throughout an earthquake sequence. Such information might evolve and change in focus over time depending on risks and needs. Focus could also be given to information that builds peoples beliefs and capacities to undertake preparedness in evolving situations. Understanding preparedness in the context of different experiences and timeframes is useful in helping update models such as the CET, where the dynamics of time might be better incorporated.

Canterbury

Evaluating the U.S. Geological Survey’s earthquake shaking hazard forecasts and their implications on seismic risk

We analyze the last six update cycles of the long-term probabilistic earthquake shaking hazard forecast from the U.S. Geological Survey (USGS) and discuss the changes in hazard estimates from the 1996 to the latest 2023 update of the National Seismic Hazard Model (NSHM) for the conterminous U.S. We summarize how our understanding of earthquake hazards has evolved over the last two decades and quantify the implications of changing hazards on people, buildings, and lifeline infrastructure. The net effect of changes in hazard estimates, with reference to the mean hazard from the 2023 NSHM update that has 2% probability of exceedance in 50 years, suggests an overall increase in total geographic areas facing very strong shaking levels (modified Mercalli Intensity (MMI) of VII or more) when compared with any of the previous cycles. This increase puts ~159.2 million people (an increase from 111.8 million when using the 2018 NSHM), ~47.3 million residential buildings (an increase from 32.5 million with the 2018 NSHM), and ~$25.6 trillion of economic exposure of all buildings (an increase from $17.4 trillion when using the 2018 NSHM) at risk from earthquakes. We demonstrate that small changes in hazard estimates do not necessarily imply small changes in risk estimates. Therefore, changes in risk estimates can be used to highlight key USGS NSHM updates for improving risk mitigation efforts.

Earthquake Spectra Journal

Aeromagnetic and magnetotelluric imaging of west-central Idaho and the Stibnite-Yellow Pine mining district: A regional to district perspective

Aeromagnetic and magnetotelluric (MT) data are used to better understand the geology and mineral resources near the Stibnite-Yellow Pine mining district in central Idaho. The reduced-to-pole (RTP) transformation of regional-scale aeromagnetic data shows that allochthonous island-arc rocks west of the Salmon River suture are significantly more magnetic than the Laurentian continental rocks east of the suture and that the granitoids of the Idaho batholith have moderate to low magnetization in both early, metaluminous, and late, peraluminous phases. Application of tilt derivative to aeromagnetic data highlights major crustal-scale structures. The 5-km upward continued magnetic data indicate island-arc rocks have deep magnetic sources. The 110-km-long MT profile images resistivity structure to depths around 30 km. At shallow depths, resistivity corresponds to mapped geologic units, with moderate resistivities underlying volcanic and roof-pendant metasedimentary rocks and moderate to high resistivities occurring beneath the Idaho batholith. Crustal-scale moderate resistivities beneath the suture image the results of tectonomagmatic processes that accompanied suturing and translating allochthonous terranes. Low resistivity values beneath and fringing the batholith are derived from metasedimentary rocks that may have served as a melt source and reductant during melt generation and provided metals during later ore formation. In the Stibnite-Yellow Pine mining district, a high-resolution aeromagnetic compilation is shown to correlate with mapped lithologies and mineral deposit-related structures. The RTP transform distinguishes magnetic and nonmagnetic granitoid phases of the Idaho batholith. The tilt derivative highlights metasedimentary rocks, some of which are favorable ore hosts. The Meadow Creek fault hosts the Stibnite and Hangar Flats deposits and is imaged as a magnetic low due to hydrothermal alteration. Reconstructions of magnetic anomaly offsets and orebodies indicate around 3 km of post-95 Ma dextral separation, with some or all of the offset inferred to postdate the main Au mineralization episode (61–66 Ma).

Idaho

Parameter ESTimation with the Gauss–Levenberg–Marquardt algorithm: An intuitive guide

In this paper, we review the derivation of the Gauss–Levenberg–Marquardt (GLM) algorithm and its extension to ensemble parameter estimation. We explore the use of graphical methods to provide insights into how the algorithm works in practice and discuss the implications of both algorithm tuning parameters and objective function construction in performance. Some insights include understanding the control of both parameter trajectory and step size for GLM as a function of tuning parameters. Furthermore, for the iterative Ensemble Smoother (iES), we discuss the importance of noise on observations and show how iES can cope with non-unique outcomes based on objective function construction. These insights are valuable for modelers using PEST, PEST++, or similar parameter estimation tools.

Groundwater

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

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

Marine Geology

The Long Island Sound and Watershed Metadata map application

The Long Island Sound and its watershed encompass an area of about 17,000 square miles and include the Connecticut, Housatonic, and Thames Rivers, which all drain to the sound. Dozens of organizations from government agencies, nonprofits, and Tribal Nations have developed projects and monitoring programs to analyze and protect the water resources of the watershed and sound. The abundance of data and lack of an existing searchable index require a centralized metadata repository to allow users to find water resources data more efficiently. The U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency and the Long Island Sound Study, has created an interactive map application to visualize and search for metadata information across organizations working to monitor and protect the Long Island Sound.

Connecticut, Massachusetts, New Hampshire, New Yor

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California

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

Quality of groundwater used for domestic supply in the Gilroy-Hollister basin and surrounding areas, California, 2022

More than 2 million Californians rely on groundwater from domestic wells for drinking-water supply. This report summarizes a 2022 California Groundwater Ambient Monitoring and Assessment Priority Basin Project (GAMA-PBP) water-quality survey of 33 domestic and small-system drinking-water supply wells in the Gilroy-Hollister Valley groundwater basin and the surrounding areas, where more than 20,000 residents are estimated to utilize privately owned domestic wells. The study area includes the Llagas subbasin in the north, the North San Benito subbasin in the south, and the surrounding uplands. The study was focused on groundwater resources used for domestic drinking-water supply, which are mostly drawn from shallower parts of aquifer systems rather than those of groundwater resources used for public drinking-water supply in the same area. This assessment characterized the quality of ambient groundwater in the aquifer before filtration or treatment, rather than the quality of drinking water delivered to the tap. To provide context, the measured concentrations of constituents in groundwater were compared to Federal and California State regulatory and non-regulatory benchmarks for drinking-water quality. A grid-based method was used to estimate the areal proportions of groundwater resources used for domestic drinking wells that have water-quality constituents present at high concentrations (above the benchmark), moderate concentrations (between one-half of the benchmark and the benchmark for inorganic constituents, or between one-tenth of the benchmark and the benchmark for organic constituents), and low concentrations (less than one-half or one-tenth the benchmark for inorganic and organic constituents, respectively). This method provides statistically representative results at the study-area scale and permits comparisons to other GAMA-PBP study areas. In the study area, inorganic constituents in groundwater were greater than regulatory benchmarks (U.S. Environmental Protection Agency [EPA] or State of California maximum contaminant levels [MCLs]) for public drinking-water quality in 24 percent of domestic groundwater resources. The inorganic constituents present at concentrations greater than MCLs for drinking water were nitrate (as nitrogen), barium, chromium, and selenium. Total dissolved solids (TDS) or manganese were present at concentrations greater than the secondary maximum contaminant levels (SMCLs) that the State of California uses as aesthetic-based benchmarks in 48 percent of domestic groundwater resources. No volatile organic compounds or pesticide constituents were present at concentrations greater than regulatory benchmarks. Total coliform bacteria and enterococci were detected in 4 percent of domestic groundwater resources. Per- and polyfluoroalkyl substances (PFAS) were detected in 19 percent of domestic groundwater resources, and 10 percent had concentrations greater than recently enacted (April 2024) EPA MCLs. Physical and chemical factors from natural and anthropogenic sources that could affect the groundwater quality were evaluated using results from statistical testing of associations between constituent concentrations and potential explanatory variables. In this study, relevant physical factors include well construction characteristics, groundwater age, site proximity to groundwater recharge or discharge zones, and potential sources of contamination. Relevant chemical factors include the initial chemistry of the recharge water, the mineralogy of the aquifer sediments, and the subsequent shifts in chemistry as biologic and geologic reactions alter groundwater in the subsurface. Nitrate concentrations were correlated to agricultural land use, distance from the boundary of the Gilroy-Hollister Valley groundwater basin, and the proportion of modern (post-1950s) water captured by the well. Denitrification under anoxic redox conditions can mitigate some nitrate derived from fertilizer application. Total dissolved solids primarily were derived from water-rock interactions with soils and aquifer materials in the study area, but there were high concentrations where agricultural practices contributed additional TDS. Mineralogy of aquifer sediments and rocks also affect barium, selenium, boron, and chromium concentrations in the Gilroy-Hollister Valley groundwater basin. PFAS were positively correlated with urban land use and the proportion of modern water captured by the well.

California