USGS Science⌕ Search

SEARCH · USGS Science

Results for “One Ecosystem”

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,216 records · Page 41Linked to original sources
USGS2024

The U.S. Geological Survey National Water Quality Network—Groundwater—2023

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, eight of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,114 wells in the network.

conterminous United States↗
USGS2024

Debris-flow entrainment modelling under climate change: Considering antecedent moisture conditions along the flow path

Debris-flow volumes can increase along their flow path by entraining sediment stored in the channel bed and banks, thus also increasing hazard potential. Theoretical considerations, laboratory experiments and field investigations all indicate that the saturation conditions of the sediment along the flow path can greatly influence the amount of sediment entrained. However, this process is usually not considered for practical applications. This study aims to close this gap by combining runout and hydrological models into a predictive framework that is calibrated and tested using unique observations of sediment erosion and debris-flow properties available at a Swiss debris-flow observation station (Illgraben). To this end, hourly water input to the erodible channel is predicted using a simple, process-based hydrological model, and the resulting water saturation level in the upper sediment layer of the channel is modelled based on a Hortonian infiltration concept. Debris-flow entrainment is then predicted using the RAMMS debris-flow runout model. We find a strong correlation between the modelled saturation level of the sediment on the flow path and the channel-bed erodibility for single-surge debris-flow events with distinct fronts, indicating that the modelled water content is a good predictor for erosion simulated in RAMMS. Debris-flow properties with more complex flow behaviour (e.g., multiple surges or roll waves) are not as well predicted using this procedure, indicating that more physically complete models are necessary. Finally, we demonstrate how this modelling framework can be used for climate change impact assessment and show that earlier snowmelt may shift the peak of the debris-flow season to earlier in the year. Our novel modelling framework provides a plausible approach to reproduce saturation-dependent entrainment and thus better constrain event volumes for current and future hazard assessment.

Illgraben basin↗
USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗
USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2026

Scenarios to assess the future water availability in the Mississippi River Valley Alluvial Aquifer for the Cache River and Grand Prairie Regions of Arkansas

The U.S. Geological Survey, as part of the Arkansas Groundwater Initiative, developed forecast scenarios using previously calibrated MODFLOW 6 groundwater models that focused on the Cache and Grand Prairie Critical Groundwater Areas to assess the impact of future climate and water management strategies on the Mississippi River Valley alluvial aquifer. A Soil Water Balance model was used to forecast recharge and irrigation water use. The forecast scenario period was from January 1, 2019, through December 31, 2055, with monthly stress periods. Twenty scenarios were simulated and included seven alternate climate forecasts, five 13 general groundwater pumping reduction scenarios (round 1), and groundwater pumping reduction scenarios by crop type and for the Bayou Meto Water Management Project and Grand Prairie Area Demonstration Project (round 2). Declines in saturated thickness within the Cache Critical Groundwater Area were larger for 18 of the 20 scenarios as compared to outside of the Critical Groundwater Area. The largest average increase in saturated thickness inside the Critical Groundwater Area was 6.4 m which occurred for the round 1, 50 percent reduction scenario. Automatic reductions in groundwater pumping by MODFLOW 6 in the Cache simulation ranged from 0.02 to 13.1 percent of total groundwater pumping. For the Grand Prairie model domain, the average change in saturated thickness of the Mississippi River Valley alluvial aquifer inside the Critical Groundwater Area for the forecast period ranged between -6.6 to 1.7 m. The average saturated thickness of the Mississippi River Valley alluvial aquifer inside the Grand Prairie Critical Groundwater Area declined for 16 of the 20 scenarios. The average reduction in requested groundwater pumping for all scenarios inside the Grand Prairie Critical Groundwater Area was 25.1 percent, and the largest reduction was 46.5 percent.

ESS Open Archive↗
USGS2025

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts↗
USGS2025

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California↗
USGS2026

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗
USGS1998

Technique for simulating peak-flow hydrographs in Maryland

The efficient design and management of many bridges, culverts, embankments, and flood-protection structures may require the estimation of time-of-inundation and (or) storage of floodwater relating to such structures. These estimates can be made on the basis of information derived from the peak-flow hydrograph. Average peak-flow hydrographs corresponding to a peak discharge of specific recurrence interval can be simulated for drainage basins having drainage areas less than 500 square miles in Maryland, using a direct technique of known accuracy. The technique uses dimensionless hydrographs in conjunction with estimates of basin lagtime and instantaneous peak flow. Ordinary least-squares regression analysis was used to develop an equation for estimating basin lagtime in Maryland. Drainage area, main channel slope, forest cover, and impervious area were determined to be the significant explanatory variables necessary to estimate average basin lagtime at the 95-percent confidence interval. Qualitative variables included in the equation adequately correct for geographic bias across the State. The average standard error of prediction associated with the equation is approximated as plus or minus (+/-) 37.6 percent. Volume correction factors may be applied to the basin lagtime on the basis of a comparison between actual and estimated hydrograph volumes prior to hydrograph simulation. Three dimensionless hydrographs were developed and tested using data collected during 278 significant rainfall-runoff events at 81 stream-gaging stations distributed throughout Maryland and Delaware. The data represent a range of drainage area sizes and basin conditions. The technique was verified by applying it to the simulation of 20 peak-flow events and comparing actual and simulated hydrograph widths at 50 and 75 percent of the observed peak-flow levels. The events chosen are considered extreme in that the average recurrence interval of the selected peak flows is 130 years. The average standard errors of prediction were +/- 61 and +/- 56 percent at the 50 and 75 percent of peak-flow hydrograph widths, respectively.

Maryland↗
USGS2024

Status and understanding of groundwater quality in the Mojave Basin Domestic-Supply Aquifer study unit, 2018—California GAMA Priority Basin Project

Groundwater quality in the western part of the Mojave Desert in San Bernardino County, California, was investigated in 2018 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The Mojave Basin Domestic-Supply Aquifer study unit (MOBS) region was divided into two study areas—floodplain and regional—to assess differences between the two major aquifers used for drinking water supply in the area. This assessment characterized the quality of ambient groundwater and not the quality of treated drinking water. The study included three components: (1) a status assessment, which characterized the quality of groundwater resources used for domestic drinking-water supply in the floodplain and regional study areas; (2) a brief understanding assessment, which evaluated factors that could potentially affect the quality of groundwater used by domestic wells in the region; and (3) a comparative assessment between the groundwater resources used by domestic wells and public-supply wells in the two study areas. The domestic-well assessment was based on data collected by the U.S. Geological Survey from 48 domestic wells in January–May 2018. The public-supply assessment was based on data for samples from 322 public-supply wells in 2008–18, either collected by the U.S. Geological Survey or compiled from the California State Water Resources Control Boards Division of Drinking Water publicly available database. Concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. Relative concentrations, defined as the measured concentration divided by the benchmark concentration, were classified as high (greater than 1.0), moderate (greater than 0.5 for inorganic constituents or 0.1 for organic and special-interest constituents, and not high), or low (concentrations lower than moderate). The floodplain and regional study areas were divided into 15 and 35 grid cells, respectively, and grid-based methods were used to compute the areal proportions of the two study areas with high, moderate, or low relative concentrations of individual constituents and classes of constituents. For the domestic-supply assessment, one or more inorganic constituents with health-based benchmarks were detected at high relative concentrations in 58 percent of the regional study area and 13 percent of the floodplain study area. The inorganic constituents with health-based benchmarks detected at high relative concentrations in the regional study area were arsenic, chromium and hexavalent chromium, fluoride, adjusted gross alpha particle activity, uranium, molybdenum, strontium, and nitrate; only arsenic was detected at high relative concentrations in the floodplain study area. One or more inorganic constituents with secondary maximum contaminant level benchmarks were detected at high concentrations in 15 and 6.7 percent of the regional and floodplain study areas, respectively. The constituents detected at high relative concentrations in the regional study area were total dissolved solids, chloride, sulfate, and iron; only total dissolved solids and sulfate were detected at high relative concentrations in the floodplain study area. Organic constituents were not detected at moderate or high relative concentrations in either the regional or floodplain study areas. Volatile organic compounds were detected at low relative concentrations in 21 and 27 percent of the regional and floodplain study areas, respectively, and pesticides were detected at low relative concentrations in 9.1 and 20 percent of the regional and floodplain study areas, respectively. The only individual organic constituent detected in more than 10 percent of either study area was the trihalomethane trichloromethane. Total coliform bacteria were detected in 15 and 27 percent of the grid wells in the regional and floodplain study areas, respectively. The greater prevalence of high relative concentrations of many inorganic constituents in the regional study area compared to the floodplain area likely indicates the greater diversity of geologic material at depth in aquifer material and generally finer-grained alluvium compared to the floodplain study area combined with generally older groundwater that has had more contact time with aquifer materials. In general, trace element concentrations (1) increased with increasing groundwater age, (2) increased with distance from recharge sources in the mountains, and (3) increased with closer proximity to some types of geological units. In general, groundwater from domestic wells in the floodplain study area is young, with most samples containing a component of modern groundwater based on tritium and unadjusted carbon-14 activities, whereas groundwater from domestic wells in the regional study area generally is old, with most samples having unadjusted carbon-14 ages of 5,000–40,000 years. Public-supply wells in MOBS generally were deeper than domestic wells and presumably are in contact with older, more weathered alluvium that may have more mobile trace elements, such as arsenic or uranium. However, only 26 percent of the public-supply regional study area had high relative concentrations of inorganic constituents, compared to 58 percent for the domestic regional study area. The percentages of the public-supply and domestic floodplain study areas with high relative concentrations of inorganic constituents were 11 and 13 percent, respectively. The ages of groundwater used by public-supply and domestic wells in each study area were similar, which was not expected given the greater depth of the public-supply wells. Three potential factors may contribute to these results: (1) greater spatial footprint of domestic well network, which may result in domestic wells pumping groundwater from fractured bedrock or mineralized areas not used by public-supply wells; (2) greater pumping rates in public-supply wells, resulting in more water being withdrawn from coarse-grained, heterogeneous alluvium than finer-grained layers, which may have higher concentrations of (or more mobile) inorganic constituents; and (3) a greater degree of well management with public-supply wells, which may include pausing use of or decommissioning wells if treating or blending water is not feasible to lower constituent concentrations.

California↗
USGS2026

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California↗
USGS1998

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada↗
USGS1998

Water resources of Mellette and Todd counties, South Dakota

Mellette and Todd Counties are located in south-central South Dakota and have a combined area of 2,694 square miles. The White River and its tributaries, which include the Little White River, drain Mellette County and about one-half of Todd County. Tributaries to the Niobrara River, which include the Keya Paha River, drain the other one-half of Todd County. The average discharge of the Little White River is about 56 cubic feet per second as the river enters Todd County and is about 131 cubic feet per second as it discharges to the White River in northern Mellette County. The average discharge of the Keya Paha River just outside Todd County is about 39 cubic feet per second. The average annual runoff for Mellette and Todd Counties ranges from 0.94 to 2.36 inches based on records from nine streamflow-gaging stations in and near the counties. The average annual runoff is 1.62 inches, which compares with the average annual precipitation of about 19 inches. In Todd County, shallow wells completed in the alluvial, Ogallala, Arikaree, and White River aquifers generally can supply water that has low concentrations of dissolved solids, is fresh, and is soft to moderately hard. Ground water from shallow aquifers is limited in Mellette County; therefore, deep wells, often greater than 1,000 feet, are sometimes installed. The Pierre Shale often is used to supply rural domestic and stock wells in Mellette County even though well yields are low and the water has high dissolved solids, is moderately saline, and is very hard. Alluvial aquifers are present in both counties and store an estimated 1.6 million acre-feet of water. The water quality of the alluvial aquifers is dependent on the underlying deposits, and generally the water has low concentrations of dissolved solids, is fresh, and is soft to moderately hard where underlain by the Ogallala and Arikaree Formations; has moderate concentrations of dissolved solids, is slightly saline, and is hard where underlain by the White River Group; and has high concentrations of dissolved solids, is saline, and is very hard where underlain by the Pierre Shale. Also, yields often are lower where the alluvial aquifers are underlain by the Pierre Shale. The Ogallala aquifer is present in only Todd County, and the Arikaree aquifer is present throughout most of Todd County and southwestern and south-central Mellette County. The Ogallala aquifer contains an estimated 17 million acre-feet of water in storage, and the Arikaree aquifer contains an estimated 50 million acre-feet of water in storage. Both aquifers generally are suitable for irrigation, and yields from these aquifers are sometimes greater than 1,000 gallons per minute. Nitrate concentrations in 13 out of 92 water samples collected from the Ogallala aquifer exceeded the Primary Drinking Water Maximum Contaminant Level (MCL) of 10 milligrams per liter. In 11 out of 46 samples collected from the Arikaree aquifer, arsenic concentrations exceeded the MCL of 50 micrograms per liter. The White River aquifer, where present, is usually the shallowest source of ground water in Mellette County. The White River aquifer also is used in northern Todd County where the Ogallala and Arikaree aquifers are not present. The White River aquifer contains an estimated 50 million acre-feet of water in storage. Reported yields from the aquifer range from 1 to 30 gallons per minute, which generally is insufficient to support irrigation in most areas. However, yields are sufficient for livestock-watering and rural-domestic purposes. In both counties, the Pierre Shale is the shallowest bedrock aquifer and is exposed at the land surface throughout most of Mellette County. This aquifer is used primarily in Mellette County. Although the aquifer contains an estimated maximum of 1.5 million acre-feet of water in storage, it is not a viable source of ground water because the aquifer is relatively impermeable, yields are low, and water usually can be obtained from shallower sources, especially in Todd County. Reported yields from the Pierre Shale aquifer range from 1 to 8 gallons per minute. Because few test holes and wells penetrate below the Pierre Shale, little is known about the extent of the deeper bedrock aquifers. All wells completed in the Dakota Sandstone, Inyan Kara, and Minnelusa and Madison aquifers in the counties are used for stock-watering purposes. High concentrations of dissolved solids and hard water are characteristic of the water quality in the bedrock aquifers. Depths to the top of the deeper bedrock aquifers range from 1,270 feet to greater than 2,000 feet below land surface.

South Dakota↗
USGS2025

Geochemical and hydrological investigations of historical data collected at the Lee Acres Landfill and Giant Bloomfield Refinery, New Mexico, 1985–2020

The Lee Acres Landfill and Giant Bloomfield Refinery are adjacent properties near the City of Farmington, New Mexico, each having undergone monitoring and remediation related to historical site activities. At the landfill, site cleanup has included the installation of a capillary barrier over former liquid waste lagoons and periodic monitoring of groundwater elevations and groundwater quality. At the refinery, remediation has focused on several petrochemical and crude oil release areas and included soil excavation, groundwater treatment, and regular monitoring of groundwater elevations and quality. Groundwater at both sites has higher concentrations of volatile organic compounds and trace metals than background aquifer concentrations. In 2022, the U.S. Geological Survey compiled the Lee Acres-Giant Bloomfield Refinery Database (LAGBRD), which contains publicly available groundwater-elevation data and organic and inorganic groundwater-quality data from both sites, spanning from 1985 to 2020. Data from the LAGBRD and precipitation data from other sources were used to better understand the cause of relatively high manganese concentrations observed in some groundwater wells at the site through comparison of groundwater chemistry to chemical end members, interpretation of spatial and temporal patterns in the groundwater chemistry, and interpretation of groundwater flow properties. In this study, elevated chloride concentrations in groundwater downgradient from the landfill have been attributed to landfill leachate based on the temporal and spatial variability of chloride concentrations and chloride-to-bromide ratios. Installation of a capillary barrier and surface-water runoff controls at the landfill in 2005 appears to have altered infiltration patterns at that site, resulting in a decrease in chloride at some wells but an increase in chloride and dissolved manganese at others. The timing and relation among groundwater elevation, chloride concentration, and manganese concentration suggest that leachate stored in the vadose zone provides a continued source of contamination to groundwater.

New Mexico↗
USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗
USGS2025

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

Restoring dryland water cycles for precipitation feedback and climate stability; a review

Drylands across the globe are experiencing intensifying water scarcity, land degradation, and hydroclimatic extremes. This review integrates evidence from multidecadal field studies, hydrologic monitoring, geomorphic and ecological assessments, remote sensing, and land–atmosphere science to evaluate how restoration influences key components of the terrestrial water cycle. Low-tech natural infrastructure in dryland streams (NIDS)—including check dams, leaky weirs, one-rock dams, and gabions—has emerged as a promising but under-synthesized nature-based solution for restoring hydrologic function in these environments. We describe the mechanisms through which these interventions modify runoff detention, infiltration, sediment and alluvial storage, shallow-groundwater recharge, vegetation recovery, and surface-energy partitioning, and we summarize outcomes across diverse dryland settings. Findings consistently show increased water residence time, enhanced soil-moisture storage, expanded riparian vegetation, extended flow duration, and shifts toward greater latent-heat flux—producing localized cooling and strengthened ecohydrological feedbacks. Building on these localized effects, we articulate a hypothesis that links the spatial extent of restoration, the density of NIDS per unit drainage area, and the magnitude of the latent-to-sensible-heat contrast generated by wetter post-rainfall conditions. Specifically, we hypothesize that when NIDS are implemented at densities permitted by topography and across areas large enough to maintain elevated soil moisture after storm events, the resulting increases in latent heat flux, surface cooling, and boundary-layer moistening may enhance moisture convergence and boundary-layer development, potentially increasing the likelihood or stability of convective precipitation, analogous to how reductions in these processes have contributed to regional drought intensification. These land–atmosphere feedbacks remain untested at scale but represent an important research Frontier. By integrating hydrologic, geomorphic, ecological, and atmospheric perspectives, this review provides a comprehensive framework for considering how low-tech, landscape-scale interventions can strengthen watershed resilience and contribute to climate-relevant nature-based solutions.

Arizona↗
USGS2026

Climate-adaptive urban planning: Quantitative assessment of drought impacts and practical strategies for climate-resilient urban green spaces

Urban green spaces (UGSs) are vital for enhancing a city’s resilience and livability; however, their functionality is increasingly jeopardized by drought, particularly in water-scarce regions. This study evaluates drought impact on UGSs in Metropolitan Adelaide, Australia, a representative semi-arid urban system, using satellite-derived Normalized Difference Vegetation Index (NDVI) time-series data spanning 2000–2020. Vegetation dynamics were analyzed through Seasonal-Trend decomposition using Loess (STL), standardized anomaly assessment, lagged Pearson correlation, Ordinary Least Squares (OLS) regression, and Mann–Kendall trend analysis. To isolate climatically sensitive signals, 29 urban lawn patches were examined separately from mixed urban canopy, given their shallow root systems and direct dependence on surface moisture. NDVI declined by approximately 0.09 units during the Millennium Drought (2001–2009), with summer greenness deficits reaching 24% below the 20-year benchmark. Temperature was the dominant driver of lawn NDVI variability (r = −0.863, R 2 = 74.5%), substantially exceeding the effect of rainfall (r = 0.156, R 2 = 2.4%). El Niño–Southern Oscillation (ENSO) cycles modulated vegetation responses, with La Niña years supporting recovery and El Niño years amplifying decline. Post-drought recovery remained incomplete, with NDVI deficits of 8–20% persisting through 2020; full recovery was observed only in 2017, coinciding with the highest recorded summer rainfall. No significant directional trend was detected over the full study period (Mann–Kendall τ = 0.005, p = 0.908). These findings demonstrate that heat, rather than water limitation alone, is the primary driver of vegetation stress in urban systems, highlighting the benefits of integrated management strategies that address both warming and moisture deficits to sustain urban green infrastructure under future climate conditions. We introduce the concept of “urban greenery drought,” referring to a form of vegetation stress in managed urban landscapes where greenness is reduced primarily by elevated temperature and atmospheric demand despite water availability.

Adelaide↗