USGS ScienceSearch

SEARCH · USGS Science

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

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

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

1,706 records · Page 24Linked to original sources

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra

Cancer risk and estimated lithium exposure in drinking groundwater in the US

Importance Lithium is a naturally occurring element in drinking water and is commonly used as a mood-stabilizing medication. Although clinical studies have reported associations between receiving lithium treatment and reduced cancer risk among patients with bipolar disorder, to our knowledge, the association between environmental lithium exposure and cancer risk has never been studied in the general population. Objectives To evaluate the association between exposure to lithium in drinking groundwater and cancer risk in the general population. Design, Setting, and Participants This cohort study included participants with electronic health record and residential address information but without cancer history at baseline from the All of Us Research Program between May 31, 2017, and June 30, 2022. Participants were followed up until February 15, 2023. Statistical analysis was performed from September 2023 through October 2024. Exposure Lithium concentration in groundwater, based on kriging interpolation of publicly available US Geological Survey data on lithium concentration for 4700 wells across the contiguous US between May 12, 1999, and November 6, 2018. Main Outcome and Measures The main outcome was cancer diagnosis or condition, obtained from electronic health records. Stratified Cox proportional hazards regression models were used to estimate the hazard ratios (HRs) and 95% CIs for risk of cancer overall and individual cancer types for increasing quintiles of the estimated lithium exposure in drinking groundwater, adjusting for socioeconomic, behavioral, and neighborhood-level variables. The analysis was further conducted in the western and eastern halves of the US and restricted to long-term residents living at their current address for at least 3 years. Results A total of 252 178 participants were included (median age, 52 years [IQR, 36-64 years]; 60.1% female). The median follow-up time was 3.6 years (IQR, 3.0-4.3 years), and 7573 incident cancer cases were identified. Higher estimated lithium exposure was consistently associated with reduced cancer risk. Compared with the first (lowest) quintile of lithium exposure, the HR for all cancers was 0.49 (95% CI, 0.31-0.78) for the fourth quintile and 0.29 (95% CI, 0.15-0.55) for the fifth quintile. These associations were found for all cancer types investigated in both females and males, among long-term residents, and in both western and eastern states. For example, for the fifth vs first quintile of lithium exposure for all cancers, the HR was 0.17 (95% CI, 0.07-0.42) in females and 0.13 (95% CI, 0.04-0.38) in males; for long-term residents, the HR was 0.32 (95% CI, 0.15-0.66) in females and 0.24 (95% CI, 0.11-0.52) in males; and the HR was 0.01 (95% CI, 0.00-0.09) in western states and 0.34 (95% CI, 0.21-0.57) in eastern states. Conclusions and Relevance In this cohort study of 252 178 participants, estimated lithium exposure in drinking groundwater was associated with reduced cancer risk. Given the sparse evidence and unknown mechanisms of this association, follow-up investigation is warranted.

contiguous United States

Introduction and methods of analysis for peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. Flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done by the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, to assess potential nonstationarity in peak flows in the north-central United States. This chapter summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flow. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The climate data consist of monthly time series estimates of temperature, precipitation, potential evapotranspiration, actual evapotranspiration, snowfall, soil moisture storage, snow water equivalent, and runoff on a 3.1-mile by 3.1-mile grid for the conterminous United States. Statistical and graphical analyses were used to investigate potential changes in hydrology and climate. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. To examine potential causal drivers of changes, the climate data were analyzed graphically and statistically. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Iowa, Illinois, Michigan, Missouri, Minnesota, Mon

Preserving and increasing water resources—Natural infrastructure in dryland streams in Baja California Sur, Mexico

The Los Planes watershed of Baja California Sur, Mexico, and its underlying aquifer are experiencing groundwater decline owing to low average annual rainfall (28.1 centimeters per year) and rising water demand from population growth and agricultural activities. This decline in water availability can lead to desertification—a process that changes arable land to desert by degrading soil and vegetation—and can pose serious challenges to livelihoods that depend on the land. To address these issues, a ranch in the Los Planes watershed has installed many natural infrastructures in dryland streams (NIDS) in channels for soil and water conservation. In 2022, the U.S. Geological Survey (USGS) began working with regional researchers and land managers to investigate the effects of NIDS on natural biological, geochemical, and physical processes and determine the efficacy of NIDS for water augmentation in the Los Planes watershed. The USGS also worked with local academic institutions and nonprofit organizations to create public educational opportunities focused on the area’s hydrogeology. These and other collaborative efforts with the U.S. Water Partnership and Innovaciones Alumbra aim at enhancing water resources in the Baja California Sur region and promoting water security and safeguarding community well-being. Resume La cuenca de Los Planes, ubicada en Baja California Sur, México, y su acuífero subyacente, están sufriendo una disminución de las aguas subterráneas debido a la baja precipitación media anual (28.1 centímetros por año) y la alta demanda de agua por parte de una población creciente y la actividad agrícola. Esta disminución de la disponibilidad de agua puede conducir a la desertificación—un proceso que por medio de la degradación del suelo y la vegetación convierte a la tierra cultivable en desierto—representando un serio desafío para los medios de vida de las personas. Para abordar estos problemas, un rancho en la cuenca de Los Planes ha instalado numerosas obras de Infraestructura Natural en Arroyos de Tierras Áridas (INATS) para conservación del suelo y del agua. En 2022, el Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) comenzó a trabajar con investigadores regionales y gestores de tierras para estudiar los efectos de INATS en los procesos biológicos, geoquímicos y físicos, y determinar su eficacia en el aumento de los recursos hídricos en la cuenca de Los Planes. El USGS se ha asociado con instituciones académicas y organizaciones locales sin fines de lucro para crear oportunidades educativas públicas centradas en la hidrogeología de la zona. Estos y otros esfuerzos colaborativos con la Asociación del Agua de Estados Unidos (U.S. Water Partnership) e Innovaciones Alumbra, tienen como objetivo mejorar el uso de los recursos hídricos en la región de Baja California Sur, promover la seguridad hídrica y proteger el bienestar de la comunidad.

Baja California Sur

ECCOE Landsat quarterly Calibration and Validation report—Quarter 1, 2024

Executive Summary The U.S. Geological Survey Earth Resources Observation and Science Calibration and Validation (Cal/Val) Center of Excellence (ECCOE) focuses on improving the accuracy, precision, calibration, and product quality of remote-sensing data, leveraging years of multiscale optical system geometric and radiometric calibration and characterization experience. The ECCOE Landsat Cal/Val Team continually monitors the geometric and radiometric performance of active Landsat missions and makes calibration adjustments, as needed, to maintain data quality at the highest level. This report provides observed geometric and radiometric analysis results for Landsats 8 and 9 for quarter 1 (January–March), 2024. All data used to compile the Cal/Val analysis results presented in this report are freely available from the U.S. Geological Survey EarthExplorer website: https://earthexplorer.usgs.gov . This quarterly report is the third to include analysis results for Landsat 9, which was launched in September 2021. The inclusion of Landsat 9 analysis results was dependent on two factors: a complete reprocessing of the Landsat 9 data archive and enough time elapsing to begin formulating lifetime trends. In April 2023, all Landsat 9 image data acquired since the satellite’s launch were reprocessed to take advantage of calibration updates identified by the ECCOE Landsat Cal/Val Team. Additional information about the Landsat 9 reprocessing effort is available at https://www.usgs.gov/landsat-missions/news/upcoming-reprocessing-all-landsat-9-data . Additional information about Landsat 9 prelaunch, commissioning, and early on-orbit imaging performance is available at https://www.mdpi.com/journal/remotesensing/special_issues/15B4V2K92K . This quarterly report is the first to not include analysis results for Landsat 7 because Enhanced Thematic Mapper Plus imaging was suspended on January 19, 2024, after the satellite transitioned into full sunlight. The satellite has been drifting since early 2022 after being lowered from the nominal orbit altitude, and the transition into full sunlight is a result of the satellite operating in its extended science mission. Additional information about the imaging suspension is available at https://www.usgs.gov/landsat-missions/news/landsat-7-imaging-suspended . Additional information about the Landsat 7 extended science mission is available at https://www.usgs.gov/landsat-missions/landsat-7-extended-science-mission .

Open-File Report

ECCOE Landsat quarterly Calibration and Validation report—Quarter 2, 2024

Executive Summary The U.S. Geological Survey Earth Resources Observation and Science Calibration and Validation (Cal/Val) Center of Excellence (ECCOE) focuses on improving the accuracy, precision, calibration, and product quality of remote-sensing data, leveraging years of multiscale optical system geometric and radiometric calibration and characterization experience. The ECCOE Landsat Cal/Val Team continually monitors the geometric and radiometric performance of active Landsat missions and makes calibration adjustments, as needed, to maintain data quality at the highest level. This report provides observed geometric and radiometric analysis results for Landsats 8 and 9 for quarter 2 (April–June) of 2024. All data used to compile the Cal/Val analysis results presented in this report are freely available from the U.S. Geological Survey EarthExplorer website: https://earthexplorer.usgs.gov . This is the fourth quarterly report to include analysis results for Landsat 9, which was launched in September 2021. The inclusion of Landsat 9 analysis results was dependent on two factors: a complete reprocessing of the Landsat 9 data archive and enough time elapsing to begin formulating lifetime trends. In April 2023, all Landsat 9 image data acquired since the satellite’s launch were reprocessed to take advantage of calibration updates identified by the ECCOE Landsat Cal/Val Team. Additional information about the Landsat 9 reprocessing effort is available at https://www.usgs.gov/landsat-missions/news/upcoming-reprocessing-all-landsat-9-data . Additional information about Landsat 9 prelaunch, commissioning, and early on-orbit imaging performance is available at https://www.mdpi.com/journal/remotesensing/special_issues/15B4V2K92K . This is the second quarterly report that does not include analysis results for Landsat 7 because Enhanced Thematic Mapper Plus imaging was suspended on January 19, 2024, after the satellite transitioned into full sunlight. The satellite has been drifting since early 2022 when it was lowered from the nominal orbit altitude, and the transition into full sunlight is a result of the satellite operating in its extended science mission. Additional information about the imaging suspension is available at https://www.usgs.gov/landsat-missions/news/landsat-7-imaging-suspended . Additional information about the Landsat 7 extended science mission is available at https://www.usgs.gov/landsat-missions/landsat-7-extended-science-mission .

Open-File Report

Investigation of drivers and sources of potential harmful algal blooms in Lake Anna, Virginia, 2023-24

Algal blooms in freshwater reservoirs are an increasing concern because of the risk to aquatic ecosystem health, recreational use, and water quality. Among the taxa that comprise an algal bloom, cyanobacteria are of particular concern due to the potential to produce toxic compounds which can cause acute and chronic illness in humans if ingested through contaminated water or shellfish. When toxins are above a toxin-specific threshold in the water column, the Commonwealth of Virginia categorizes the bloom as a harmful algal bloom (HAB), distinguishing it from an algal bloom without toxins present. Toxins are of particular concern in bodies of water which serve as major recreational areas, such as Lake Anna in central Virginia. Lake Anna is a 27 km-long and 47 km2 reservoir in central Virginia. Since the Commonwealth of Virginia revised the HAB monitoring framework in 2018, algal bloom advisories for potentially toxigenic blooms have been issued in the western, riverine zone of the reservoir. There were detections of microcystin and anatoxin production but no detections of cylindrospermopsin and saxitoxin production in Lake Anna; however, toxin concentrations were low, at values less than 0.65 µg/L. In 2023, the U.S. Geological Survey, in cooperation with Virginia Department of Environmental Quality, initiated a 19-month intensive and extensive study of western Lake Anna and its tributaries (North Anna River and Pamunkey Creek) to (1) characterize the algal community, biomass, and toxin production; (2) identify potential drivers leading to algal bloom initiation, persistence, and decline; and (3) assess the primary sources of the drivers contributing to algal blooms. Continuous and discrete monitoring of water quality within the lake, watershed inputs, meteorological conditions, and algal communities were monitored. Multiple modeling scenarios identified that the formation and persistence of algal blooms in Lake Anna are driven by interactions among chemical, macro- and micronutrient, and physical factors. Key drivers include total nitrogen, total phosphorus, water temperature, suspended sediment concentration, wind speed, copper, iron, molybdenum, nickel, sodium, and alkalinity. During this study, winter stormflow events delivered the highest loading of macro- and micronutrients from the tributaries to the lake. The lake is a sink for these nutrients in the water column or in the lake-bed sediments. During high-occupancy and recreation periods on the lake in the summer months, boating activity can promote internal loading by the resuspension of the lake-bed sediments, thereby increasing the nutrient availability for algal uptake. Thermal stratification and hypoxic conditions were observed, which can contribute to the passive release of nutrients into the water column, further facilitating algal bloom proliferation. These processes may create a positive feedback loop that supports the formation, persistence, and annual recurrence of algal blooms in the western region of Lake Anna. Continued water-quality monitoring of early warning indicators may improve the detection and forecasting of algal growth in the lake and help managers proactively manage water quality and algal growth.

EarthArXiv

Identifying transportation data and system needs for a Federal lands transportation data platform

Executive Summary Modern transportation and land-use planning efforts include information from many sources to address topics such as safety, efficiency, commercial, and social needs. This wide breadth of topics provides opportunities for collaboration and development of common tools for diverse users. In many cases, different information systems provide the spatial data and geographic content necessary for transportation and land-use planners to consider multiple lines of evidence. The Federal Highway Administration Office of Federal Lands Highway (FLH) and Federal Land Management Agency partners use detailed spatial and quantitative data to inform transportation decisions. However, logistic challenges to data sharing exist because data are often managed by separate agencies; data-exchange frameworks and interagency data agreements are insufficient; and consistency from aggregated data requires maintenance, coordination, and supporting infrastructure. The FLH and U.S. Geological Survey collaboratively examined (1) use and availability of spatial data for transportation planning and (2) a possible mechanism to use more shared and consistent data in a common planning environment. The goals of this collaborative effort were to describe data needs from the perspective of planners and to identify opportunities for shared data resources. Results presented here focus on two workshops and a subsequent investigation of data and tools available from partner agencies. The objectives of this report are to (1) describe information used in transportation planning with geographic data; (2) identify spatially explicit data that inform transportation plans and could be shared among all partners; and (3) describe current platforms, planning and administrative opportunities, and potential barriers to developing an integrated planning tool. Key information and data needs were identified in three major classes: system, user, and influential factors. System data are parts of the transportation network and information about the condition of individual segments and the network. User data provide details about the function of the system and insights into potential needs; for example, user trips between source and destinations inform road and network demands that can lead to congestion and safety issues (in the future, user data might also include scenarios and projections based on land-use plans). Influential data represent social and environmental factors that influence transit demands and network conditions. These factors could be popular locations or seasonal events that influence demand and congestion; wildlife habitat or migration intersections that affect safety and management priorities; or geologic features that influence hazards, maintenance, and safety. Responses described here provide specific information for web-tool design and give a framework for interagency communication and cooperation to address specific information needs for integrated planning. Existing web-mapping and web-services, and the data that inform them, are also described. Commonly, these data are created and published by one agency, and the core users are outside of that agency; for example, threatened species distributions are published by the U.S. Fish and Wildlife Service for consideration by planners in advance of National Environmental Policy Act (42 U.S.C. 4321 et seq.) evaluation. This report is provided to inform FLH leaders and Federal Land Management Agency partners by articulating user needs and requirements for integrated planning tool(s). Programmers creating a secure web-based data-sharing platform (with data-viewing, -analysis and -download functions) can use the information presented here to organize data and user interfaces. This integrated perspective can help FLH and Federal Land Management Agency partners develop transportation networks that better serve the needs of people in local communities and across States and the Nation.

Open-File Report

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater. This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks. Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

California

Geology, mineralogy, and cassiterite geochronology of the Ayawilca Zn-Pb-Ag-In-Sn-Cu deposit, Pasco, Peru

The Ayawilca deposit in Pasco, Peru, represents the most significant recent base-metal discovery in the central Andes and one of the largest undeveloped In resources globally. As of 2018, it hosts an 11.7 Mt indicated resource grading 6.9% Zn, 0.16% Pb, 15 g/t Ag, and 84 g/t In, an additional 45.0 Mt inferred resource grading 5.6% Zn, 0.23% Pb, 17 g/t Ag, and 67 g/t In, and a separate Sn-Cu-Ag inferred resource of 14.5 Mt grading 0.63% Sn, 0.21% Cu, and 18 g/t Ag. Newly obtained U–Pb dates for cassiterite by LA-ICP-MS (22.77 ± 0.41 and 23.05 ± 2.06 Ma) assign the Ayawilca deposit to the Miocene polymetallic belt of central Peru. The polymetallic mineralization occurs as up to 70-m-thick mantos hosted by carbonate rocks of the Late Triassic to Early Jurassic Pucará Group, and subordinately, as steeply dipping veins hosted by rocks of the Pucará Group and overlying Cretaceous sandstones-siltstones of the Goyllarisquizga Group. Relicts of a distal retrograde magnesian skarn and cassiterite (stage pre-A) were identified in the deepest mantos. The volumetrically most important mineralization at Ayawilca comprises a low-sulfidation assemblage (stage A) with quartz, pyrrhotite, arsenopyrite, chalcopyrite, Fe-rich sphalerite, and traces of stannite and herzenbergite. Stage A sphalerite records progressive Fe depletion, from 33 to 10 mol% FeS, which is compatible with the observed transition from low- to a subsequent intermediate-sulfidation stage (B) marked by the crystallization of abundant pyrite and marcasite. Finally, during a later intermediate-sulfidation stage (C) sphalerite (up to 11 mol% FeS), galena, native bismuth, Cu-Pb-Ag sulfosalts, siderite, Mn-Fe carbonates, kaolinite, dickite, and sericite were deposited. This paragenetic evolution shows striking similarities with that at the Cerro de Pasco Cordilleran-type polymetallic deposit, even if at Ayawilca stage C did not reach high-sulfidation conditions. The occurrence of an early retrograde skarn assemblage suggests that the manto bodies at Ayawilca formed at the transition between distal skarn and skarn-free (Cordilleran-type) carbonate-replacement mineralization. Mineral assemblages define a T- f S 2 evolutionary path close to the pyrrhotite-pyrite boundary. Buffering of hydrothermal fluids by underlying Devonian carbonaceous phyllites of the Excelsior Group imposed highly reduced conditions during stage A mineralization (log f O 2 < − 30 atm). The low f O 2 favored efficient Sn mobility during stages pre-A and A, in contrast to other known ore deposits in the polymetallic belt of central Peru, in which the occurrence of Sn minerals is minor. Subsequent cooling, progressive sealing of vein walls, and decreasing buffering potential of the host rocks promoted the shift from low- (stage A) to intermediate-sulfidation (stages B and C) states. LA-ICP-MS analyses reveal significant In contents in Fe-rich sphalerite (up to 1.7 wt%), stannite (up to 1908 ppm), and chalcopyrite (up to 1185 ppm). The highest In content was found in stage A sphalerite that precipitated along with chalcopyrite and stannite, thus pointing to the early, low-sulfidation assemblage as prospective for this high-tech metal in similar mineral systems. Indium was likely incorporated into the sphalerite crystal lattice via Cu + + In 3+ ↔ 2 Zn 2+ and (Sn, Ge) 4+ + (Ga, In) 3+ + (Cu + Ag) + ↔ 4 Zn 2+ coupled substitutions. Indium incorporation mechanisms into the stannite and chalcopyrite crystal lattices remain unclear.

Pasco

Spatial distribution of API gravity and gas/oil ratios for petroleum accumulations in Upper Cretaceous strata of the San Miguel, Olmos, and Escondido Formations of the south Texas Maverick Basin—Implications for petroleum migration and charge history

The Maverick Basin of south Texas is currently undergoing active exploration and production of gas and oil from tight sandstone reservoirs. The most productive tight sandstones in the basin are in the Upper Cretaceous San Miguel, Olmos, and Escondido Formations. These units are second only to the Eagle Ford Shale in terms of cumulative production volumes. The structural history of the Maverick Basin, from rifting to subsidence to exhumation, has had a profound effect on the characteristics of these reservoirs and the petroleum resources contained therein. This U.S. Geological Survey review of the production history of these strata reflects a recent shift from conventional production to horizontal drilling (unconventional) that exploits low permeability reservoirs in previously overlooked areas of existing oil and gas fields in southern Texas, typically outside of established field boundaries. To investigate the physical properties of the Maverick Basin hydrocarbon accumulations, this case study compiled American Petroleum Institute (API) gravity measurements and calculated cumulative gas/oil ratios (GOR) for thousands of producing wells from the San Miguel, Olmos, and Escondido Formations. Maps were generated from the compiled well production data to show the spatial heterogeneity of API gravity and GOR values for the three formations within the Maverick Basin and immediately outside the basin to the northeast. Within the Maverick Basin, the spatial patterns of API gravity values indicate lighter oils downdip towards the southern basin edge. GOR values indicative of wet and dry gases within the basin are seen interspersed, with values that correspond to black and heavy oils. Differences in the spatial patterns of the petroleum properties within the Maverick Basin are interpreted as effects of Eocene basin inversion caused by Laramide orogenic deformation, and the resulting reservoir exhumation of basin strata. East of the Maverick Basin, spatial distributions of API gravity and GOR values show progressively heavier oils updip to the northwest, grading to dry gases downdip to the southeast, which correlates to the oil and gas windows of the underlying Eagle Ford Shale. Correlation of API gravity and GOR values from the San Miguel, Olmos, and Escondido Formations with thermal maturity data from the Eagle Ford Shale suggests that the Eagle Ford Shale is the petroleum source, and that petroleum migration was approximately vertical for areas to the east of the Maverick Basin. The discontinuity of API gravity and GOR properties within the Maverick Basin implies a complex petroleum charge history, possibly involving the remigration of petroleum and the addition of petroleum from other source intervals in Mexico, to the southwest. Depressurization of exhumed, overpressured reservoirs of the San Miguel, Olmos, and Escondido Formations can explain the intermittent occurrence of gas production throughout the southern Maverick Basin by exsolution of gas from formation brines and the resulting dry gas flushing of hydrocarbon-charged reservoirs. The introduction of dry gas through flushing can, in turn, explain why the patterns of API gravity and GOR values are so dissimilar in the Maverick Basin. This process has implications for possible future production of unconventional resources from undiscovered tight-gas reservoirs in strata of the San Miguel, Olmos, and Escondido Formations, and a different approach to petroleum exploration may be needed in the Maverick Basin relative to exploration techniques applied in other basins within the northern Gulf of Mexico.

Texas

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming

Evidence of nitrate attenuation in intertidal and subtidal groundwater in a subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2015–16

Nitrogen dynamics in intertidal and nearshore subtidal groundwater (subterranean estuary) adjacent to the Seacoast Shores peninsula, Falmouth, Massachusetts, were investigated during 2015–16 by the U.S. Geological Survey. The peninsula is a densely populated residential area with septic systems and cesspools that are substantial sources of nitrogen to groundwater. The study area is in the Eel River, an estuarine saltwater embayment connected to the ocean adjacent to the western shore of the peninsula, that was the subject of an earlier study by Colman and others (2018, https://doi.org/10.3133/sir20185095 ) on nitrogen transport and transformations in groundwater between onshore and offshore locations. The previous study documented the distribution of nitrate concentrations and nitrate attenuation reactions in fresh groundwater beneath the peninsula and the estuary. The current study extended those observations with more detailed sampling and analysis of shallow groundwater from wells near discharge sites beneath the estuary. The current field investigation included sampling of existing wells and installation and sampling of clusters of wells and temporary sampling points in the subterranean estuary, including (1) shallow transects 0.3 to 1.2 meters (m) deep extending from 1 to 13.5 m offshore and (2) deeper wells (from 1.83 to 4.88 m deep) extending from 4.3 to 14.3 m offshore. Measurements of hydraulic-head gradients 2–5 m below the sediment/water interface in the intertidal and nearshore subtidal zones indicated that groundwater flow generally was upwards (towards the estuary) under all tide conditions in October 2016. The magnitude of the gradient was greatest during low tide conditions, indicating that groundwater discharge likely decreased during high tides. Measurements of specific conductance in shallow groundwater in the subterranean estuary in three transects perpendicular to shore were consistent with the existence of saltwater flow cells (infiltration of overlying saline water, mixing with fresh groundwater, and discharge to the overlying saline water) in the intertidal and nearshore subtidal regions. The size of these flow cells was variable in space and time and dependent on the elevation of the tide (spring or neap). At this location in the Eel River subterranean estuary, and offshore to at least 13.5 m, offshore flow of fresh groundwater apparently prevented a deeper saltwater wedge from discharging to the surface. Nitrate concentrations in shallow groundwater (30 to 122 centimeters [cm] depth) were variable in space and time, ranging from not detectable to 600 micromoles per liter (μmol/L) (8.4 milligrams per liter as N), and were highest in June 2016 at depths from 61 to 122 cm below the sediment/water interface and from 4 to 9 m offshore. Nitrate generally was not detectable in saline shallow groundwater at 30-cm depth or at any depth from 30 to 122 cm from 10 to 13.5 m offshore. Dissolved oxygen concentrations were suboxic (less than 16 μmol/L) in 60 percent of the sampled subterranean groundwater beneath the intertidal and subtidal zones. In the remaining sites, the range of dissolved oxygen concentrations was from 18 to 272 μmol/L and the median concentration was 43 μmol/L. Evidence for microbial nitrate reduction (denitrification and possibly anammox) was provided by the distribution of the reaction product nitrogen gas (excess N 2 , or N 2MIC ), as determined from analysis of the dissolved nitrogen gas and argon gas (Ar) concentrations in groundwater samples. Excess nitrogen gas provided evidence for nitrate reduction in shallow groundwater below the subtidal and, to a lesser extent, intertidal zones adjacent to the Seacoast Shores peninsula. These zones, where evidence for nitrate reduction was detected, were in fresh and brackish groundwater near subtidal or intertidal saltwater cells where discharging fresh groundwater mixed with infiltrating saline water. Infiltrating seawater may have supplied organic carbon, one of several potential electron donors that are required for denitrification. Other potential electron donors, such as organic carbon, iron, manganese, hydrogen, methane, ammonium, elemental sulfur, or sulfide phases, may have been supplied by the estuarine sediments. Drainage from surface runoff near the shore also may have supplied organic carbon to fresh groundwater near the intertidal saltwater cell. The highest amounts of nitrate converted to excess nitrogen gas were estimated to be in the range of 230 to 430 μmol/L in nearly fresh groundwater near the subtidal saltwater cell at depths of 61 to 122 cm below the sediment/water interface and from 10 to 13.5 m offshore. Evidence of denitrification within 10 m of the shore was sparse (generally limited to less than 50 μmol/L of N 2 -N) despite the presence of high nitrate concentrations. The spatial distribution of estimated nitrate reduction in the intertidal and nearshore subtidal fresh and brackish groundwater may be related to local variability in the distribution of reactive electron donors in those zones. Variations in the amount of nitrate reduction to nitrogen gas were not clearly related to potential aqueous electron donors such as dissolved organic carbon, nor to potential reaction products such as alkalinity, but may have been controlled by combinations of aqueous and solid-phase reactants. The distribution of relatively shallow fresh groundwater containing nitrate could indicate potential nitrate discharge areas in the lower intertidal zone and uncertain locations farther offshore; however, the data did not extend all the way to the sediment/water interface or to the offshore freshwater limit. This study confirmed substantial loss of nitrate from some of the fresh and brackish groundwater in shallow subestuarine sediments prior to discharge but did not quantify how much nitrate eventually discharged to the estuary.

Massachusetts

Insights and strategic opportunities from the USGS 2024 Per- and Polyfluoroalkyl Substances (PFAS) Interagency Workshop

Introduction In 2021, the U.S. Geological Survey (USGS) published Circular 1490 titled, “Integrated Science for the Study of Perfluoroalkyl and Polyfluoroalkyl Substances (PFAS) in the Environment: A Strategic Science Vision for the U.S. Geological Survey” (Tokranov and others, 2021). Circular 1490 was created to be a resource for USGS scientists prioritizing and planning research related to per- and polyfluoroalkyl substances (PFAS) and to be a guide for developing partnerships with other scientists, State and Federal agencies, and stakeholders engaged in PFAS research and management and mitigation of the environmental and human-health effects of PFAS. This USGS PFAS Strategic Science Vision document was intended to be the foundation for a “living strategic vision,” periodically providing updates on the state of USGS PFAS research, emerging PFAS data gaps and needs, and progress on interagency and stakeholder PFAS partnerships and priorities. To meet this objective, the USGS planned to host an Interagency and Stakeholder PFAS Workshop every 2–3 years. During September 10–12, 2024, the USGS hosted the first Interagency and Stakeholder PFAS Workshop in Reston, Virginia. The Workshop brought together experts from other Federal agencies (U.S. Environmental Protection Agency, National Institute of Environmental Health Sciences, Food and Drug Administration, Department of Defense [Air Force, Army]), State agencies (Washington Fish and Wildlife, Virginia Department of Transportation), and academia (Harvard University, University of Maryland) to address key challenges relating to the measurement and modeling of PFAS and the implications for environmental health. Participants engaged in in-depth discussions centered around six pivotal topics related to PFAS: (1) sampling protocols, methods and interpretation; (2) environmental sources, source apportionment, and occurrence; (3) environmental fate and transport; (4) human and wildlife exposure routes and risk; (5) bioconcentration, bioaccumulation, and biomagnification; and (6) ecotoxicology and effects. Each topic had three breakout sessions. A recurrent theme of workshop discussions was how data on a nationwide scale for PFAS occurrence in various environmental matrices, including air, water, food crops, biota, soil, and streambed sediment could help to advance scientific understanding. Participants noted significant geospatial data gaps, particularly in the midwestern and southern United States and the Pacific Northwest. PFAS data collection tends to be more robust along the eastern seaboard and in California. Participants stressed how enhancing the integration of large and small datasets across various agencies could help to support national scale understanding of PFAS. To address these gaps, attendees suggested leveraging datasets from Federal entities like the USGS and the U.S. Department of Defense, State agencies, and municipal utility services to develop predictive contaminant detection and transport models. Improved coordination between water quality programs and USGS research could help to facilitate access to valuable data, leading to comprehensive databases that inform PFAS point (wastewater treatment plants and landfills) and nonpoint (runoff from land, atmospheric deposition, food packaging) sources, environmental transport mechanisms, environmental detection and concentrations, potential exposure routes, and health effects on different biota, including humans. A specific request was made to develop a map demarking the depth of modern (1953 or later) groundwater, which is susceptible to surface-derived anthropogenic (that is, human-made) contamination, based on tritium-age dating. Emphasis was placed on incorporation of hydrology, groundwater flow paths, groundwater–surface water interactions, and landscape factors in predictive statistical models as a step to improve contaminant source identification and tracking. Molecular fingerprinting approaches garnered attention as techniques to link specific PFAS mixtures detected in a sample to environmental sources and levels in biota (Dávila-Santiago and others, 2022). Integrating data from abiotic (that is, water, soil, and air) and biotic (that is, living organisms) systems identified as a research opportunity. For example, understanding the composition of soils and sediments, which include a mixture of mineral, plant, and animal components, could advance understanding of exposure pathways. The discussions highlighted opportunities to explore and understand the potential redistribution and biotic exposures of PFAS from biosolid and wastewater treatment plant effluent land application practices, in addition to atmospheric releases and discharges from landfill and wastewater treatment plants. Participants identified research gaps surrounding how these sources may contribute to contamination and may affect surrounding ecosystems, including a better definition of anthropogenic background concentrations. Moving forward, the collection of co-occurrence data was noted as a means to improve understanding of complex mixtures and to leverage companion modeling efforts focused on areas with high and low contamination levels to identify areas of concern and unaffected resources. Participants emphasized how centralized USGS databases and the establishment of sample-metadata archives can help to ensure that samples are preserved and accessible for future research. In conclusion, the workshop participants identified opportunities to bridge data gaps and improve measurement techniques, modeling frameworks, databases, and communication, to enhance the understanding of PFAS and their effects on environmental and human health. Upon completion of the workshop, participants indicated an interest in developing strategic data collection, modeling, and analytical approaches to address these challenges.

Open-File Report

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York