Geology ReportsSearch

SEARCH · Geology Reports

Results for “Physical Review Letters”

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.

303 records · Page 17Linked to original sources

Technique for estimating magnitude and frequency of peak flows in Maryland

A convenient and reliable technique for estimating flood magnitudes is required for effective flood-plain management and for the efficient design of bridges, culverts, embankments, and flood-protection structures. Methods are presented for estimating peak-flow magnitudes of selected frequencies, ranging from 2 to 500 years, for all nontidal drainage basins in Maryland. The methods were developed by generalized least-squares regression techniques using data from 219 gaged basins in and near Maryland. The State is divided into five hydrologic regions: the Appalachian Plateaus and Allegheny Ridges region, the Blue Ridge and Great Valley region, the Piedmont region, the Western Coastal Plain region, and the Eastern Coastal Plain region. These regions correspond to the physiographic provinces of the State, with the exceptions that (1) the Coastal Plain Province is divided into two hydrologic regions, and (2) there is no distinct hydrologic region corresponding to the Valley and Ridge Province as it is divided into its constituent Allegheny Ridges and Great Valley subdivisions. Sets of equations for calculating peak discharges based on physical basin characteristics are provided for each of the regions. Based on the peak-flow equations, methods for estimating peak flows are presented for ungaged and gaged streams in Maryland. The methods and equations are supported by generalized least-squares analysis of basin and flood-frequency characteristics data from 219 drainage basins in and near Maryland. estimates for each of the five regions are calculated using combinations of the fol-lowing basin characteristics: drainage area, forest cover, basin relief, carbonate rock coverage, storage, and runoff-curve number. Drainage area contributes to the estimate in all five study regions. Carbonate rock coverage is used only in the Blue Ridge and Great Valley region. Storage and runoff-curve number are used solely in the Eastern Coastal Plain region. All other basin characteristics are used in two or more regions. Standard errors of estimate for the regression equations range from 19 to 31 percent in the Appalachian Plateaus and Allegheny Ridges region, 34 to 47 percent in the Blue Ridge and Great Valley region, 33 to 48 percent in the Piedmont region, 45 to 64 percent in the Western Coastal Plain region, and 36 to 42 percent in the Eastern Coastal Plain region.

Maryland

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California

Properties of solid bitumen formed during hydrous, anhydrous, and brine pyrolysis of oil shale: Implications for solid bitumen reflectance in source-rock reservoirs

Solid bitumen reflectance (BR o in %) is widely used as a thermal maturity proxy in source-rock reservoirs, yet solid bitumen texture and physical attributes may be affected by varying environmental constraints during its formation, e.g., water saturation, mineral catalysis, or salinity. Herein we investigated the development of solid bitumen properties during artificial maturation using three diverse (lacustrine to marine) oil shale samples containing abundant amorphous H-rich sedimentary organic matter (bituminite). The oil shales were treated via pyrolysis (320 °C, 72 h) using hydrous, anhydrous, and brine conditions, causing the development of a newly formed solid bitumen in the experiment residues from the thermal conversion of bituminite. The properties of the newly formed solid bitumen then were evaluated via geochemical screening tests, optical and electron microscopy, and infrared spectroscopy. Experimental residues also were treated via solvent extraction, allowing characterization of the effects of extraction on solid bitumen. Results showed that bituminite with higher resin and asphaltene components evolved to a solid bitumen with higher reflectance (as a percentage of the original value) when exposed to the same heat treatment, relative to bituminite with higher saturate and aromatic components. Aromatization of solid bitumen also was generally more pronounced in the presence of deionized water relative to anhydrous conditions, supporting prior observations. These results suggest the compositions of primary sedimentary organic matter, and the local concentration of water, affect the development of solid bitumen during thermal advance, potentially explaining the origin of multiple solid bitumen populations with varying reflectance in source-rock reservoir samples or in a single microscope field. Experiments using brine were inconclusive regarding enhanced/decreased aromatization, which could be attributed to the salinities of the brines used in the tests. Extraction of residues caused a consistent increase in BR o values, suggesting that migrating or expelled oils could cause an increase in BR o via natural solvation and absent an increase in temperature. This work provides insights into the development and evolution of BR o in source-rock reservoirs as a function of the composition of the original bituminite and changing environmental conditions, with potentially broad application for petroleum prospecting and resource estimation.

Colorado

Tephra from Kīlauea’s 2008–2018 lava lake eruption—Proximal deposits and dispersal characteristics

A network of ten buckets was established early in the 2008–2018 summit eruption at Kīlauea to collect proximal tephra ejected from the new, informally named the “Overlook crater”; the buckets were emptied on most days of the eruption thereafter. This report summarizes the results of more than 2,400 different sampling intervals (most 1–3 days long) during the eruption, focusing on the physical and dispersal characteristics of the tephra deposits. The network was within about 300 meters south of the vent to capture tephra dispersed by the dominant northeast trade wind. The juvenile tephra mainly reflected spattering at the southeast (SE) sink, a downwelling area in the southeastern part of the lava lake in the Overlook crater that remained in the same area throughout the eruption, with admixtures of solid rock and secondary minerals derived from the wall of the crater. The proportion of juvenile material to lithic material ranged widely early in the eruption but was generally greater than 90 percent for the last 6 years of the eruption as lake level rose and the crater walls decreased in height and became more stable. The accumulation rate of tephra at each bucket was strongly dependent on the location of the bucket and reflects the interplay between lava lake level and wind direction and speed. The mass per unit area (m/a) of collected tephra was a maximum of about 97 kilograms per square meter for the entire eruption, equivalent to a thickness of about 75 millimeters. Thirty-two explosive events with a volcanic explosivity index of −2 to −4 deposited much of the tephra in the network. Fifteen of these rock-fall-induced events occurred within 6 days of one another. Pele’s hair, a volcanic glass predominantly associated with quieter activity at the vent in the final half of the eruption, was dispersed more than 60 kilometers downwind from the lava lake and formed a nearly continuous deposit near the Overlook crater. This eruption was probably the most frequently sampled long-lasting eruption in history, but most of the deposits are ephemeral. The collected samples, although generally of small mass, are retained by the Hawaiian Volcano Observatory and are available for detailed study.

Hawaii

Introduction to recommended capabilities and instrumentation for volcano monitoring in the United States

Introduction The National Volcano Early Warning System (NVEWS) was authorized and partially funded by the U.S. Government in 2019. In response, the U.S. Geological Survey (USGS) Volcano Hazards Program asked its scientists to reflect on and summarize their views of best practices for volcano monitoring. The goal was to review and update the recommendations of a previous report (Moran and others, 2008) and to provide a more detailed analysis of capabilities and instrumentation for monitoring networks for U.S. volcanoes. This Scientific Investigations Report and its chapters reflect those USGS scientists’ views and summaries and will serve as a guide for future network upgrades funded through NVEWS. Given the well-documented hazards posed by volcanoes to population centers and aviation (for example, Blong, 1984; Scott, 1989; Neal and others, 1997, 2019; Guffanti and others, 2010; Shroder and Papale, 2014; Prata and Rose, 2015; Palmer, 2020), volcano monitoring is critical for ensuring public safety and for mitigating the impacts of volcanic activity. Accurate and timely forecasts are facilitated by well-designed monitoring networks that are in place long enough to allow for background behavior to be recognized and understood. Because precursory signals may be limited and unrest may progress rapidly to an eruption, our goal is to deploy monitoring systems that enable detection of the reactivation of dormant volcanoes as early as possible, allowing for public safety and risk mitigation. NVEWS planning is also informed by the results of Ewert and others (2005, 2018), whereby 161 U.S. volcanoes are currently categorized and ranked commensurate with their relative threat. In each chapter, author(s) considered the need for some redundancy of instrumentation and telemetry, given the likelihood of occasional equipment failure, particularly in extreme and remote environments. Establishing digital telemetry networks requires advanced planning, sighting, radio-shot testing, and, inevitably, troubleshooting in the field. This is harder to achieve rapidly during a crisis; thus, an important goal for monitoring U.S. volcanoes is to establish digital telemetry backbones with redundancy and extra capacity to absorb additional instruments should a volcano begin to exhibit signs of unrest (fig. A1). The National Telecommunications and Information Administration (NTIA) imposed new regulations in the United States, eliminating the use of older analog radios for many purposes, which had been one previous means for redundant data delivery. However, the resulting conversion from analog to digital systems usefully enables stations to accommodate new and multivariate real-time data streams (for example, Global Navigation Satellite System [GNSS] receivers, infrasound arrays, gas spectrometers, visible and infrared cameras, and broadband seismometers). We note that other USGS and broader national and international hazard programs can leverage NVEWS instrumentation plans. Examples of this include the following: Improved seismic coverage of volcanoes will increase the capability of the USGS Earthquake Hazards Program to detect and locate earthquakes, estimate ground shaking, and provide timely early warnings through the ShakeAlert Earthquake Early Warning System (Given and others, 2018). The National Oceanic and Atmospheric Administration’s Tsunami Program will benefit from additional seismic stations, particularly within the sparsely instrumented Aleutian Islands, Northern Mariana Islands, and American Samoa. Infrasound stations can detect signals from landslides, debris flows and lahars, floods, and weather events, providing benefits to the National Weather Service and the USGS Landslide Hazards Program.

Scientific Investigations Report

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S 2- ) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

California

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

Hyperspectral narrowband imaging spectroscopy: A new paradigm for Earth observation

This editorial introduces the Special Issue entitled “Hyperspectral Narrowband Imaging Spectroscopy: A New Paradigm for Earth Observation” in the August 2026 issue of Photogrammetric Engineering & Remote Sensing (PE&RS), the flagship journal of the American Society for Photogrammetry and Remote Sensing (ASPRS). This volume represents the fourth dedicated hyperspectral special issue published in PE&RS, following earlier contributions by Thenkabail et al. (2025, 2024a, 2024b), and continues ASPRS’s commitment to advancing cutting‑edge imaging spectroscopy research and its applications across Earth system science. Remote sensing is undergoing one of the most profound transformations in its history. The emergence of hyperspectral narrowband (HNB) imaging spectroscopy data, capable of acquiring hundreds of contiguous, narrow spectral bands, has shifted the discipline from observing Earth in a handful of broad spectral windows to capturing continuous spectral signatures of the Planet. This transition marks a decisive break from the multispectral paradigm that has dominated satellite remote sensing for nearly five decades, driven by the advent of new orbital imaging spectrometers such as EnMAP, PRISMA, and NASA’s EMIT, and by the forthcoming Surface Biology and Geology (SBG) mission ((Pires Silva et al., 2026; Bourriz et al., 2025; Thenkabail et al., 2025; Chabrillat et al., 2024; Aneece et al., 2024; Dave et al., 2024; Thenkabail et al., 2024a; Thenkabail et al., 2024b; Thenkabail, 2024a; Thenkabail, 2024b; Thompson et al., 2022; Kokaly et al., 2022; Aneece & Thenkabail, 2022; Cawse Nicholson et al., 2021; Guanter et al., 2021; Vangi et al., 2021; Thenkabail et al., 2021). These missions (e.g., Table 1) deliver unprecedented spectral fidelity, improved signal to noise ratios, and global coverage capabilities, enabling a new era of quantitative, spectroscopy based Earth observation. Where multispectral broadbands (MBBs) provide only a few discrete measurements along the electromagnetic spectrum, HNB systems deliver rich, diagnostic information that enables scientists to characterize Earth’s surface with unprecedented biochemical, biophysical, and structural detail (Figure 1a, 1b). The implications for environmental monitoring, agriculture, water resources, and mineral exploration are profound. Several overarching themes emerge: • Spectral fidelity matters. The ability to preserve subtle absorption features is essential for mineral mapping, vegetation trait retrieval, and biochemical modeling. • AI and deep learning are indispensable. From destriping to classification, modern analytics must be scalable, label‑efficient, and capable of exploiting the full spectral–spatial richness of HNB data. • Physics‑based and data‑driven approaches must converge. Radiative transfer models such as PROSAIL, enhanced with localized soil parameterizations, remain foundational for biophysical retrievals and model‑based inference. • Dimensionality reduction and feature extraction are critical. Techniques such as L1‑ISOMAP demonstrate that intelligent manifold learning can unlock the structure of fused, high‑dimensional datasets. • Next‑generation architectures must be interactive and multimodal. ICTNet exemplifies the future of hyperspectral classification: hybrid, synergistic, and capable of modeling both local textures and global spectral dependencies.

Photogrammetric Engineering and Remote Sensing (PE

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

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

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

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

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

New York

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

Wind River subbasin restoration: Annual report of U.S. Geological Survey activities January 2024 through December 2024

We sampled juvenile wild Oncorhynchus mykiss (Steelhead Trout) in headwater streams of the Wind River, Washington, to characterize population attributes and investigate life-history metrics, particularly migratory patterns, and early life-stage survival. Additionally, this sampling provides an index of abundance of non-native Brook Trout Salvelinus confluentus . We used passive integrated transponder (PIT) tagging and a series of instream PIT-tag interrogation systems (PTISs) to track juveniles and adults. The Wind River subbasin is considered a wild Steelhead refuge by Washington Department of Fish and Wildlife (WDFW). No hatchery Steelhead Trout have been released in the Wind River subbasin since 1997, and hatchery adults are estimated at less than one percent of spawners in most years. Over 25 years of Steelhead Trout status and trend monitoring and research in the subbasin is contributing to understanding of population response to numerous restoration actions in the subbasin, including removal of Hemlock Dam from Trout Creek in 2009, which had an inadequate adult ladder and contributed to increased water temperatures, to near lethal levels for salmonids. Data from our study, and companion work by WDFW, are contributing to the Columbia Basin Fish and Wildlife Program (2009) Research, Monitoring, and Evaluation (RM&E) Strategy of Fish Population Status Monitoring. Specifically, this work addresses the sub-strategies of 1) Assessing the Status and Trends of Diversity of Natural Origin Fish Populations and Uncertainties Research regarding differing life histories of a wild Steelhead Trout population, 2) Assessing the Status and Trend of Adult Natural Origin Fish Populations, and 3) Monitoring and Evaluating the Effectiveness of Tributary Habitat Actions Relative to Environmental, Physical, or Biological Performance Objectives. During summer and fall 2024, we PIT-tagged 2,075 Steelhead Trout parr (age-0 and age-1), in the Trout Creek and upper Wind River watersheds. Age-0 and age-1 parr were at high abundance due to strong spawner numbers in 2023 and 2024 (Estimate of 932 in 2023; and a moderate number of spawners in 2024 though an estimate is not available at this time; 24-year median = 619; Charlie Cochran, WDFW Fish Biologist, personal commun., 2025). An additional 100 parr (age-0 and age-2 or older) were tagged to provide fish for estimating detection efficiencies at PTISs. Steelhead Trout parr were recaptured and detected through repeat headwater sampling, smolt trapping, instream PTISs and Columbia River mainstem PIT-tag detection. We maintained and upgraded six instream PTISs to detect PIT-tagged Steelhead Trout parr, smolts, and adults, providing data for population assessments, and life-cycle research. Detection data from PIT-tagged adult Steelhead Trout at our four primary PTISs allow assessment of adult escapement to tributary watersheds within the Wind River subbasin. Most adult Steelhead detected were tagged by WDFW at an adult trap in the lower Wind River and some as smolts at screw traps. Adult Steelhead Trout detection efficiency estimates at our primary PTIS in Trout Creek have been greater than 99 percent during eight of the past ten years and have exceeded 97% at our primary PTIS in the Wind River during nine of the past ten years. Adult escapement estimates to tributary watersheds are helping evaluate the efficacy of the 2009 removal of Hemlock Dam from rkm 2.0 of Trout Creek. Detections at the instream PTISs have shown trends of age-0 and age-1 Steelhead Trout parr emigration from natal areas during summer and fall, in addition to the expected movement of parr and smolts in spring. Our data suggest that most fish from a cohort that migrate downstream will do so at age-1 for additional rearing in the subbasin. We have estimated that from 15% to 56% of parr tagged as age-0 fish in headwater areas make downstream migrations at age-1 for additional rearing (median = 29%). We have estimated that from 1% to 27% of Steelhead Trout parr, tagged as age-1 fish, make downstream migrations during fall (median = 3.5%): this is especially pronounced in the upper Wind River portion of the watershed (median 10%). These findings raise questions about parr rearing habitat and whether migrations are density- or habitat-quality driven, and answers to such questions are long-term goals of this study as active and passive habitat restoration occurs. Repeat sampling at sites in the subbasin within and between years has enabled assessment of juvenile Steelhead Trout growth patterns. Growth rates (relative change in weight) of age-0 PIT-tagged parr during summer have been similar across the subbasin, though slightly lower in the Trout Creek watershed (Trout Creek sites, 0.57 - 1.25%/day, Wind River sites = 0.94 – 2.48%/day). Summer growth rates have been lower for age-1 parr in the Trout Creek watershed than the upper Wind River watershed (Trout Creek sites, = 0.03 - 0.14%/day, Wind River sites = 0.05 – 0.19%/day). Yearly relative growth was similar across the subbasin for both age-0 and age-1 tagged parr. During 2024 range of mean daily relative change in weight from the previous year for age-0 tagged fish was: Trout Creek sites = 0.96 – 1.41%/day, Wind River sites = 0.80 – 1.30 %/day; and for age1 tagged fish was: Trout Creek sites = 0.37 - 0.50 %/day, Wind River sites = no age-1 tagged recaptures. Non-native Brook Trout are present in the subbasin, chiefly the Trout Creek watershed. Mean percent-of-catch that is Brook Trout, at four sample sites in Trout Creek, has declined from the period 1998 – 2003 to the period 2011 – 2024 (range of decrease = 19 - 82%). Percent-of-catch (of 64 samples, 52 < 10%; 3 > 20%) and number of Brook Trout at the Trout Creek sites from 2011 through 2024, though variable, have generally declined. Habitat restoration evaluation efforts are critical to ensure efficient use of money and resources but require long-term monitoring of population metrics. Detailed viable salmonid population and life-history data, such as that provided by PIT-tagging and instream PTIS networks inform fisheries policy and management and enable assessment of long-term effects of habitat restoration actions such as the removal of Hemlock Dam on Trout Creek and proposed major instream habitat restoration in the upper Wind River.

Washington