SEARCH · USGS Science

Results for “Regional Environmental Change”

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

SourceResearch connectionFrom yearThrough year

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

955 records · Page 46Linked to original sources

USGS2026

Subsurface water ice mapping on Mars: A probabilistic approach

Subsurface water ice deposits on Mars are an important resource for potential future human exploration. They are also an indicator of the planet’s past climate. However, the distribution of subsurface water ice in Mars’s midlatitudes is uncertain because spacecraft imagery cannot directly observe subsurface ice in most cases. Various spacecraft remote sensing instruments are sensitive to subsurface water ice, including thermal imaging spectrometers, radar sounders, and neutron spectrometers. Geomorphic analyses of images can also implicate subsurface ice. Building upon the data products from the Mars Subsurface Water Ice Mapping project, we provide a probabilistic framework to jointly interpret existing data and estimate the likelihood of subsurface water ice in the Martian midlatitudes between 60 ∘ S and 60 ∘ N with uncertainty. Broadly, we find that near-surface ice is likely present poleward of ∼45 ∘ in both the northern and southern hemispheres. However, closer to the equator, existing remote sensing data cannot uniquely constrain the presence of subsurface water ice. Our probabilistic results provide a framework for quantifying the abundance of ice on Mars, and our uncertainty estimates allow future analysis and exploration to target regions of high uncertainty.

Planetary Science Journal↗

USGS2023

Paleomagnetism of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia—Geologic unit correlations and geomagnetic cryptochron identifications

Paleomagnetic rock samples were collected from 173 drill sites in the Quaternary alkali basaltic volcanic field of northern Harrat Rahat, Kingdom of Saudi Arabia. Laboratory measurements on these samples established that lava flows and vent complexes—identified and mapped from field characteristics, rock types, and compositions as products of single or temporally close eruptions—typically record single, or very similar, directions of remanent magnetization. Correlations defined through geologic mapping, spatial association, geochemistry, geochronology, and identical mean remanent directions indicate at least 16 brief episodes of temporally clustered eruptions. These episodes had durations of a few centuries or less. Anomalous remanent magnetic directions were found for at least 13 mapped lavas of northern Harrat Rahat, which demonstrate that they were acquired during brief geomagnetic cryptochrons during the Brunhes Normal Polarity Chron. These uncharacteristic directions enhance the opportunity to identify common eruptive episodes, and to better understand and evaluate assessments of eruption ages based on 40 Ar/ 39 Ar geochronology. Combining paleomagnetic and regional archaeomagnetic results for the youngest eruptions allows us to evaluate their historical age assignments and, in one case, refute a previously assigned provisional age.

Harrat Rahat volcanic field↗

USGS2024

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon↗

USGS2025

The feasibility of using lidar-derived digital elevation models for gravity data reduction

Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Open-File Report↗

USGS2026

Use of precariously balanced rocks to constrain postglacial earthquake magnitudes in New England, United States

We use the age and fragility of 13 precariously balanced rocks (PBRs) in Massachusetts, New Hampshire, and southwestern Maine to evaluate the estimated probabilistic seismic hazard (PSH) from the U.S. Geological Survey’s National Seismic Hazard Model (NSHM), and also develop a map of maximum postglacial earthquake magnitudes in the region. The PBRs are erratic boulders left behind by the retreat of the Laurentide ice sheet between 13 and 18 ka, and published ice retreat maps are used to estimate their fragility ages. The fragilities and ages of the rocks appear to be consistent with the NSHM, in that the estimated earthquake magnitudes and resulting ground motions that would topple them are only expected on timescales longer than their fragility ages. This result, on the North American passive margin, contrasts with PBR studies at active plate boundaries, where the PBR constraints imply large reductions to estimated PSH. It may be that the PBRs are not old enough in the eastern United States to reveal the same PSH modeling issues as those revealed in the plate boundary settings. The study provides constraints on maximum magnitudes of postglacial earthquakes across much of Massachusetts, New Hampshire, and southwestern Maine, with the constraints limited to the ages of the PBRs.

Maine. Massachusetts, New Hampshire↗

USGS2025

Enhancing mineral systems exploration through geochronology, thermochronology, and isotope analysis: USGS Geochron and USGS Isotope databases

A mineral systems approach to mineral exploration provides a comprehensive framework for understanding ore deposit formation by examining the geodynamic, magmatic, hydrothermal, and sedimentary processes responsible for mineralization, alteration, and remobilization of economic mineral deposits. Temporal and thermal constraints on ore genesis are crucial for refining mineral system models and guiding predictive exploration strategies. Geochronology and thermochronology offer invaluable insights into the timing and thermal evolution of ore-forming processes, whereas isotopic analyses provide critical information on the source and geochemical history of ore-forming fluids. Combining these methodologies have proven highly effective for mineral exploration in regions like Australia, however, their combined application has been limited in the United States. To apply these tools to mineral systems-based exploration, the U.S. Geological Survey (USGS) has developed two products: (1) The USGS Geochron Database, and (2) the USGS Isotope Database. These databases provide centralized repositories of geo/thermochronological dates and data (Geochron Database) and both radiogenic and stable isotope data (Isotope Database) generated by the USGS and partners over the past decades. Integrating these datasets together and with traditional exploration approaches provides the mineral exploration community with powerful tools for determining the temporal and thermal histories of ore systems and identifying metallogenic source provinces.

Continental United States↗

USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗

USGS2024

Timing and geometry of the Chemehuevi Formation reveal a late Pleistocene sediment pulse into the Lower Colorado River

The Chemehuevi Formation is a distinctive 50−150-m-thick wedge-shaped Pleistocene sedimentary unit deposited by the Colorado River. It lines the perimeters of the river’s floodplains and bedrock canyons for more than 600 km between the mouth of the Grand Canyon and the delta region in the Gulf of California. The formation is composed of a basal tan to light-yellowish-brown and pale-orange mud-dominated facies overlain and interbedded by a light-yellow-brown sand-dominated facies. The unit is one of two extensively exposed aggradational packages in the Lower Colorado River corridor, in addition to a series of other smaller alluvial terrace deposits. The Chemehuevi Formation appears to represent the response of a fully integrated Colorado River system to a significant perturbation, in contrast to the Bullhead Alluvium, which is likely a unique result of Pliocene river integration. The aggradation of the Chemehuevi Formation in the Lower Colorado River corridor may be similarly due to a unique event in the Colorado River system, or it may instead be a well-preserved sedimentary sequence recording typical behavior of the Colorado River below the Grand Canyon in the late Pleistocene. As such, multiple causal mechanisms have been proposed, but no study to date has conclusively explained the Chemehuevi Formation. To help resolve its timing, duration, and origin, we applied post-infrared infrared stimulated luminescence, carbonate U-Th series, and zircon sensitive high-resolution ion microprobe U-Th series geochronology to determine the ages of key exposures of the unit over a wide spatial area. These new data demonstrate that the Chemehuevi Formation was deposited ca. 110−90 ka. The depositional ages collectively overlap, suggesting that deposition occurred rapidly relative to the resolution of the geochronometers. The new depositional timing coincides with a shift from glacial to interglacial conditions after the marine isotope stage 5-6 transition. This observation is consistent with a climate-induced sediment pulse as a causal mechanism, yet correlations with similar deposits in the Colorado River headwaters or in neighboring catchments appear elusive. Potentially, climate transitions between glacial and interglacial periods induced a sediment pulse from hillslopes of the Colorado River system that resulted in the Chemehuevi Formation. An alternative or additional explanation is that the Chemehuevi Formation represents release of lava dam−impounded sediment in the Grand Canyon. The surface geometry of the Chemehuevi Formation projects upstream to the approximate location of lava dams, and the largest possible lava dam impoundment (the Upper Prospect dam) is comparable in volume to the formation. The lava dam hypothesis appears to be a possible explanation for the Chemehuevi Formation. However, tying deposition to a specific lava dam or series of lava dams remains challenging due to discrepancies in timing and volume. The combined effects of a series of lava dams may have led to the Chemehuevi Formation, as the last Pleistocene lava dam eruption coincides with the onset of deposition. Alternatively, the formation may result from the combined effects of both regional climate transitions and the lava dams that created a transient reservoir to compound a climate transition−driven sediment pulse. The geochronologic data presented here do not allow us to distinguish between the lava dam or climate transition hypotheses but will need to be reconciled with any future proposed depositional model.

Arizona, California, Nevada↗

USGS2024

Flood-inundation maps for the Cuyahoga River in and near Independence, Ohio, 2024

Digital flood-inundation maps for a 9.9-mile reach of the Cuyahoga River in and near Independence, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The model was calibrated to the current (2024) stage-streamflow relation (rating curve 43.0) for the USGS streamgage 04208000, Cuyahoga River at Independence, Ohio. The resulting hydraulic model was then used to compute 13 water-surface profiles for water levels (flood stages) ranging from 14.00 to 26.00 feet. The flood stages range from “action stage” to above “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used in combination with a digital elevation model derived from light detection and ranging data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood mitigation strategies and to define flood hazard areas to protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid postflood recovery efforts.

Ohio↗

USGS1999

Concentrations and distribution of manmade organic compounds in the Lake Tahoe Basin, Nevada and California, 1997-99

The U.S. Geological Survey, in cooperation with the Tahoe Regional Planning Agency and the Lahontan Regional Water-Quality Control Board, sampled Lake Tahoe, major tributary streams to Lake Tahoe, and several other lakes in the Lake Tahoe Basin for manmade organic compounds during 1997-99. Gasoline components were found in all samples collected from Lake Tahoe during the summer boating season. Methyl tert -butyl ether (MTBE), benzene, toluene, ethylbenzene, and xylenes (BTEX) were the commonly detected compounds in these samples. Most samples from tributary streams and lakes with no motorized boating had no detectable concentrations of gasoline components. Motorized boating activity appears to be directly linked in space and time to the occurrence of these gasoline components. Other sources of gasoline components to Lake Tahoe, such as the atmosphere, surface runoff, and subsurface flow, are minor compared to the input by motorized boating. Water sampled from Lake Tahoe during mid-winter, when motorized boating activity is low, had no MTBE and only one sample had any detectable BTEX compounds. Soluble pesticides rarely were detected in water samples from the Lake Tahoe Basin. The only detectable concentrations of these compounds were in samples from Blackwood and Taylor Creeks collected during spring runoff. Concentrations found in these samples were low, in the 1 to 4 nanograms per liter range. Organochlorine compounds were detected in samples collected from semipermeable membrane devices (SPMD's) collected from Lake Tahoe, tributary streams, and Upper Angora Lake. In Lake Tahoe, SPMD samples collected offshore from urbanized areas contained the largest number and highest concentrations of organochlorine compounds. The most commonly detected organochlorine compounds were cis- and trans-chlordane, p , p' -DDE, and hexachlorobenzene. In tributary streams, SPMD samples collected during spring runoff generally had higher combined concentrations of organochlorine compounds than those collected during baseflow conditions. Upper Angora Lake had the fewest number of organochlorine compounds detected of all lake samples. Dioxins and furans were not detected in SPMD samples from two sites in Lake Tahoe or from two tributary streams. The number of polycyclic aromatic hydrocarbon (PAH) compounds and their combined concentrations generally were higher in samples from Lake Tahoe than those from tributary streams. Areas of high-motorized boating activity at Lake Tahoe had the largest number and highest concentrations of PAH's. PAH compounds were detected in samples from SPMD's in four of six tributary streams during spring runoff, all tributary streams during baseflow conditions, and at all lake sites. The most commonly detected PAH's in tributary streams during spring runoff were phenanthrene, fluoranthene, pyrene, and chrysene, and during baseflow conditions were phenanthrene, 1-methylphenanthrene, diethylnaphthalene, and pyrene. Upper Truckee River, which has an urban area in its drainage basin, had the largest number and highest combined concentration of PAH's of all stream samples. Bottom-sediment from Lake Tahoe had detectable concentrations of p-cresol, a phenol, in all but one sample. A sample collected near Chambers Lodge contained phenol at an estimated concentration of 4 micrograms per kilogram (µg/kg). Bottom-sediment samples from tributary streams had no detectable concentrations of organochlorine or PAH compounds. Several compounds were detected in bottom sediment from Upper Angora Lake at high concentrations. These compounds and their concentrations were p , p' -DDD (10 µg/kg), p , p' -DDE (7.4 µg/kg), 2,6-dimethylnaphthalene (estimated at 190 µg/kg), pentachlorophenol (3,000 µg/kg), and p-cresol (4,400 µg/kg).

California, Nevada↗

USGS2024

Compositional and structural mapping of Northwest Africa 15507 angrite

Angrite meteorites represent interesting sampling of planetary crustal environments. Quench-textured angrites with strong crystal zoning originated from the shallow surface region, with evidence of reducing conditions during solidification. Plutonic angrites have more coarse-grained igneous and metamorphic textures with comparatively less zoning and are interpreted as having equilibrated at greater depth. Plutonic angrites contain the minerals magnetite Fe 3 O 4 comprised of Fe 2+ and Fe 3+ , and rhönite, where Fe 3+ is required by inspection of mineral stoichiometry. NWA 15507 is a plutonic angrite with a microgabbroic texture (mean grainsize ∼1.4 mm) composed of zoned Al-Ti-augite, Ca-bearing olivine, anorthite, with accessory kirschsteinite, rhönite, hercynite, low-Ni kamacite, merrillite, Ti-free magnetite and troilite [ 1 ]. Rhönite analyzed by electron-probe microanalysis (EPMA) has the formula Ca 2.04 (Mg 0.32 Fe 2+ 4.25 Fe 3+ 0.47 Ti 0.33 Al 0.61 )(Si 3.74 Al 2.26 )O 20 , where Fe 3+ was estimated by stoichiometric analysis. During preliminary cathodoluminescence (CL) analysis, complex subgrain and oscillatory zoning was observed in the anorthite. In this study we use a combined approach of electron probe microanalysis (EPMA), cathodoluminescence (CL), electron-backscatter diffraction (EBSD), and laser ablation inductively-coupled plasma mass spectrometry (LA-ICP-MS) to further investigate the anorthite and distribution of Fe 3+ in NWA 15507.

Microscopy and Microanalysis↗

USGS2024

Flood-inundation maps for the Cuyahoga River at Jaite, Ohio, 2024

Digital flood-inundation maps for a nearly 6-mile reach of the Cuyahoga River at Jaite, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. The maps depict estimates of the extent and depth of flooding corresponding to selected water levels (stages) at USGS streamgage 04206425 on the Cuyahoga River at Jaite, Ohio. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The hydraulic model was calibrated to the current USGS streamgage data and then used to compute 15 water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from 6 to 20 feet, which correspond to below “action stage” to “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used with a geographic information system digital elevation model derived from light detection and ranging data to delineate the areas flooded at each stage. These maps, along with current stage data from the USGS streamgage and forecasted high-flow stages from the National Weather Service, can provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Ohio↗

USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗

USGS2023

The Saudi Geological Survey-U.S. Geological Survey northern Harrat Rahat project—Styles, rates, causes, and hazards of volcanism near Al Madīnah al Munawwarah, Kingdom of Saudi Arabia

Active volcanic systems pose serious hazards to people and property including inundation and incineration by lava, blanketing by tephra (volcanic ash), exposure to noxious volcanic gases, and damage from shallow earthquakes triggered by ascending molten material (magma). To improve understanding of volcanism and associated seismicity on the western Arabia Plate, the Saudi Geological Survey and the U.S. Geological Survey conducted a multi-year investigation of the northern Harrat Rahat volcanic field adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia. Project components included creation of a high-resolution digital topographic base; interpretation of eruptive history supported by detailed geologic mapping, paleomagnetism, and abundant high-precision geochronology of volcanic deposits; assessments of eruptive styles and volcanic hazards by physical volcanology; investigation of the origins of magmas in the mantle and of their differentiation in the crust revealed by chemical and isotopic petrology; gravity and magnetotelluric surveys to reveal crustal structures and to search for magma reservoirs; and regional and local seismic tomography and analyses of seismic hazards. Project results are presented in this Professional Paper as chapters written for technical scientific audiences. This initial chapter introduces the project and briefly summarizes results in plain language for readers who have more general backgrounds.

Al Madīnah al Munawwarah↗

USGS2026

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America↗

USGS2024

Deep lithospheric controls on surface deformation and seismicity around the East Anatolian Fault Zone and A3 Triple Junction

The East Anatolian Fault Zone (EAFZ) is a plate-bounding strike-slip fault capable of hosting large earthquakes, as demonstrated by the extremely damaging February 2023 M w 7.8 and M w 7.7 mainshocks of the Kahramanmaraş earthquake sequence. Deformation related to this boundary, part of the Anatolia-Arabia-Africa (A 3 ) Triple Junction, is diffuse, as was shown by part of this earthquake sequence occurring on a northern splay of the EAFZ (the Sürgü-Çardak Fault Zone; SCFZ). Controls on surface deformation are commonly linked to stress in the brittle upper crust, but the complex deformation and seismicity patterns in this region may also reflect deeper processes, such as variations in the location and extent of the strong Arabian Plate lithospheric mantle. Seismic tomography indicates that the Arabian Plate underthrusts Anatolia as far north as the SCFZ and extends as far west as the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced M w >6 earthquakes. By investigating the relationship between deformation since the inception of the EAFZ (ca. 5 Ma), seismic structure, and seismicity, we infer that the SCFZ will become the future SE boundary of the Anatolian Plate as part of the evolving A 3 Triple Junction.

GSA Today↗

USGS2025

Controls on water quality below a reclaimed surface coal mine, southeastern Montana

Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Montana↗

USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗