SEARCH · USGS Science

Results for “Applied and Environmental Microbiology”

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.

1,657 records · Page 58Linked to original sources

USGS2025

Implications of physics-based M9 ground motions on liquefaction-induced damage in the Cascadia Subduction Zone: Looking forward and backward

Given the likelihood of future M9 Cascadia Subduction Zone (CSZ) earthquakes, various estimates of the resulting, regional ground motions have been made, including a suite of 30 physics-based simulations that reflect key modeling uncertainties. However, because the last CSZ interface rupture occurred in 1700 CE, the shaking expected in such an event is especially uncertain, as are the impacts to the built and living environments. Like other coseismic impacts, soil liquefaction poses a significant threat and must be considered by any scenario study used to inform planning and response, or to focus mitigation resources. Liquefaction is also notable for its potential to “ground truth” ground-motion estimates, given that its presence or absence in the geologic record can provide constraint on the intensities of shaking in past events. It is thus an important phenomenon looking both forward and backward. Accordingly, using recent physics-based simulations, this study (1) predicts liquefaction in M9 CSZ ruptures at 400 locations in Oregon, Washington, and British Columbia (BC) using an array of cone-penetration-test based models and (2) uses paleoliquefaction evidence at ten sites spanning from Southern Oregon to Vancouver, BC to constrain possible ground-motion intensities experienced in the 1700 CE earthquake. The forward predictions indicate that liquefaction in M9 events could be pervasive in the region and affect numerous population hubs, with the potential for damage across hundreds of square kilometers. The backward analyses suggest that 1700 CE ground-motion intensities may have been less than expected from M9 simulations in some northern portions of the CSZ (e.g. Seattle), given the paucity of 1700 CE liquefaction evidence in these areas. Ultimately, further discovery and analysis of CSZ paleoliquefaction, or lack thereof, will confirm or modify this possibility and the conclusions drawn herein.

Oregon, Washington↗

USGS2026

Albuquerque Seismological Laboratory strategic vision

This circular presents a strategic outlook for the U.S. Geological Survey (USGS) Albuquerque Seismological Laboratory (ASL) for the next 10 years (2026–36). The ASL is a USGS field office in the Geological Hazards Science Center that operates portions of the Advanced National Seismic System and the Global Seismographic Network and focuses on fundamental research for instrumentation testing and data quality. The strategic outlook is categorized into two types of tasks: “Foundational Tasks” and “Aspirational Tasks.” Foundational Tasks are those that maintain the laboratory’s basic operations and services, including regional and global seismic monitoring, improving data quality, and providing instrument testing and support. A suite of Aspirational Tasks is also articulated; these can be considered priority targets of ASL that could improve ASL’s seismic monitoring capabilities and mission. Such tasks include improvements to remote stations, testing capabilities of nonseismic geophysical instruments, detection threshold monitoring, rapid aftershock deployments, and expanding seismic monitoring networks. This report was written with input from the USGS Geological Hazards Science Center, the USGS Earthquake Hazards Program (EHP), and colleagues with an interest in the work done by the ASL. Although the details of these tasks may change, this document can provide guidance on the overarching tasks at the ASL from 2026 to 2036 and an overview of the various components of the ASL and how they fit into the EHP and the Global Seismographic Network Program.

Circular↗

USGS2024

Metal mobilization from thawing permafrost to aquatic ecosystems is driving rusting of Arctic streams

Climate change in the Arctic is altering watershed hydrologic processes and biogeochemistry. Here, we present an emergent threat to Arctic watersheds based on observations from 75 streams in Alaska’s Brooks Range that recently turned orange, reflecting increased loading of iron and toxic metals. Using remote sensing, we constrain the timing of stream discoloration to the last 10 years, a period of rapid warming and snowfall, suggesting impairment is likely due to permafrost thaw. Thawing permafrost can foster chemical weathering of minerals, microbial reduction of soil iron, and groundwater transport of metals to streams. Compared to clear reference streams, orange streams have lower pH, higher turbidity, and higher sulfate, iron, and trace metal concentrations, supporting sulfide mineral weathering as a primary mobilization process. Stream discoloration was associated with dramatic declines in macroinvertebrate diversity and fish abundance. These findings have considerable implications for drinking water supplies and subsistence fisheries in rural Alaska.

Alaska↗

USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗

USGS2024

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

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

Open-File Report↗

USGS2023

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia↗

USGS2024

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

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

Open-File Report↗

USGS2025

Airborne radiometric data map alteration of porphyry copper systems in the Elkhorn district, MT

The Elkhorn district in southwestern Montana is known for past production of high-grade base and precious metals from skarn and carbonate replacement deposits. The district also hosts a porphyry Cu-Mo deposit. A radiometric survey flown over the district and surrounds provides maps of surficial concentrations of potassium, thorium, and uranium. Handheld gamma ray spectrometer measurements on altered and unaltered rocks in the district and surrounding region provide a means to ground truth the airborne data. We demonstrate that the airborne and ground-based measurements show similar map patterns. We calculate ratios of the radioelement concentrations and present a potassium enhancement map that combines potassium and ratios of potassium to thorium and uranium. The results highlight the rocks previously mapped as having K-feldspar-biotite and quartz-sericite alteration assemblages in the Elkhorn district and map additional hydrothermal systems in the region, including known alteration in the Radersburg district. The data enhancement techniques can be used as a screening tool for mapping additional porphyry copper systems.

Montana↗

USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗

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↗

USGS2024

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

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

Massachusetts↗

USGS2025

National population exposure and evacuation potential in the United States to earthquake-generated tsunami threats

Previous efforts to characterize tsunami threats to people have focused primarily on individual scenarios in specific areas but have not recognized multiple scenarios across an entire country. This study addresses this gap by quantifying population exposure and evacuation potential in the United States to 102 earthquake-related, tsunami-hazard zones, including 92 local scenarios, 8 distant scenarios, and 2 probabilistic products. Geospatial path-distance modeling quantified evacuation potential and the influence of departure delays. We focused on residents to support other national, multi-hazard risk analyses. Millions of residents are in distant-tsunami zones, and hundreds of thousands of residents are in local-tsunami zones. In 41 scenarios, there is at least one resident that may have insufficient time to evacuate before wave arrival. Tens of thousands of residents may have insufficient time to evacuate from local tsunamis that impact the U.S. Pacific Northwest or Puerto Rican coastlines. The largest improvements in evacuation potential may come from reducing departure delays in some areas but may involve vertical-evacuation structures or changing land use in other areas.

International Journal of Disaster Risk Reduction↗

USGS2025

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin↗

USGS2025

Road salt collection and redistribution at an urban rain garden on sandy soil, Gary, Indiana

Rain gardens installed as green infrastructure to divert storm runoff from entering combined sewers also collect dissolved constituents and particulates. An urban rain garden in northwestern Indiana, USA, was continuously monitored from November 2019 to May 2021 to evaluate the fate of dissolved constituents entering the rain garden in runoff. Physical and chemical properties of soils in the rain garden were also monitored, along with underlying groundwater. Linear regression models relating specific conductance to chloride concentration indicated that the 0.0371-ha (3998 square feet) rain garden collected approximately 1490 kg (3285 pounds) of road salt from the surrounding 0.2228 ha (24,500 square feet) of impervious surfaces. Soils and groundwater were seasonally affected by road salt application but carryover from year to year was not indicated. Rain garden soil permeability (5.20 × 10 −5 to 9.72 × 10 −5 m/s) remained unchanged during the study period and soil organic carbon generally increased under native vegetation. The results suggest that a rain garden built on sandy soil can divert substantial quantities of runoff and dissolved constituents from combined sewers; however, chloride is transported to sub-infrastructure groundwater that eventually discharges to adjacent waterways with concentrations lower than those observed in runoff.

Indiana↗

USGS2024

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2019–2021

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Water availability is an important issue in the Black Mesa area because of the arid climate, past industrial water use, and continued water requirements for municipal use by a growing population. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year, depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents the results of data collected as part of the monitoring program in the Black Mesa area from calendar years 2020–2021 and, additionally, uses streamflow statistics from November and December 2019. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2020, total groundwater withdrawals were estimated to be 2,680 acre-feet (acre-ft), and, in 2021, total withdrawals were estimated to be 2,570 acre-ft. Total withdrawals during 2021 were about 65 percent less than total withdrawals in 2005 because the Peabody Western Coal Company discontinued its use of water to transport coal in a coal slurry pipeline after 2005 and ceased mining operations in 2019. Owing to Navajo Nation and Hopi Reservation access restrictions during the Coronavirus pandemic, water levels were not collected from municipal wells in 2020 or 2021. Water levels measured in 2021 from wells completed in the unconfined areas of the N aquifer within the Black Mesa area showed a decline in 7 of 13 wells when compared with water levels from the prestress period (prior to 1965). The changes in water levels across all 13 wells ranged from +8.4 feet (ft) to −42.4 ft, and the median change was −0.4 ft. Water levels also showed decline in 11 of 12 wells measured in the confined area of the aquifer when compared to the prestress period. The median change for the confined area of the aquifer was −25.9 ft, with changes across all 12 wells ranging from +17.3 ft to −133.7 ft. Spring flow was measured at four springs between 2020 and 2021. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso, Arizona. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s, and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976–2021), Dinnebito Wash near Sand Springs 09401110 (1993–2020), Polacca Wash near Second Mesa 09400568 (1994–2020), and Pasture Canyon Springs 09401265 (2004–2021). Median winter flows (November through February) of each winter were used as an estimate of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Polacca Wash and Pasture Canyon Springs, whereas a decreasing trend was observed at Moenkopi Wash and Dinnebito Wash. In 2020 and 2021, water samples were collected from a total of four springs in the Black Mesa area and analyzed for selected chemical constituents. Results from the four springs were compared with previous analyses from the same springs. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. However, concentrations of chloride from Pasture Canyon Spring show a diminishing trend. Concentrations of dissolved solids, chloride, and sulfate at Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Concentrations of dissolved solids at Burro Spring have varied for the period of record, but there is no statistical trend in the data. However, concentrations of chloride and sulfate from Burro Spring show a trend towards lower concentrations.

Arizona↗

USGS2025

A spatiotemporal deep learning approach for predicting daily air-water temperature signal coupling and identification of key watershed physical parameters in a montane watershed

Seasonal shifts from runoff to groundwater dominance influence daily headwater stream temperatures, especially where local groundwater input is strong. This input buffers temperature during hot periods, supporting cold-water habitats. Recent studies use air–water temperature signal metrics to identify zones of strong stream–groundwater connectivity. While Previous studies used air–water signal ratios as proxies for groundwater influence but were limited to specific sites and periods, without dynamic forecasting. This study is the first to forecast daily A r as a spatiotemporal signal using a Graph Convolutional Network–Long Short-Term Memory (GCN-LSTM) model. The model was trained using hydroclimate data (air temperature, precipitation, shortwave radiation, streamflow) and watershed physical features (e.g., sand content, slope). Results showed high predictive skill, achieving R 2 (NSE, RMSE) of 0.86 (0.73, 0.0004) for one-day-ahead to 0.52 (0.50, 0.0009) for seven-days ahead forecasts. Prior studies often have not explicitly incorporated spatial hydrogeologic drivers, but this model explicitly incorporates them to assess their impact on A r forecasting and stream-groundwater connectivity. Feature analysis identified mean sand, elevation, slope, clay, and TWI as key predictors of A r . Stronger groundwater signals appeared in hillslopes, elevations, and tributaries, highlighting watershed influence on streamflow. However, limitations include reliance on historical air–water temperature patterns for training and limited representation of extreme climate conditions. Despite these limitations, unlike previous studies relying on measured in-situ stream and air temperature, this study forecasts A r directly from climate and physiographic features after training, avoiding in-situ data requirements. Findings aiding predictions of stream ecosystem resilience.

New York↗

USGS2025

Insights into widespread landsliding in southern Appalachia from Hurricane Helene

Between 23 and 27 September 2024, antecedent rain followed by Hurricane Helene produced one of the most damaging weather events in southern Appalachia history. The back-to-back storm events resulted in a maximum cumulative rainfall of 848 mm and hurricane-force wind gusts over 170 km/h in western North Carolina, eastern Tennessee, and southwestern Virginia. The resulting regional flooding, landslides, and tree blowdown caused over 100 fatalities, damaged or destroyed critical infrastructure and thousands of structures, and severed connectivity across the region. Over the next several weeks, a multi-agency landslide response produced a rapid hazard assessment and mapped 2217 landslides, 55% of which damaged infrastructure or property. Orographic uplift enhanced rainfall, resulting in concentrated landsliding along the ~250 km swath of the Blue Ridge escarpment in western North Carolina. Landslides initiated predominantly on windward-facing (southeast-facing) slopes, and localized clustering of initiation points indicated a strong influence of hillslope-scale meteorological and geomorphic factors. Many shallow landslides mobilized into larger, highly mobile, and damaging debris flows that graded into floods. Here, we put our preliminary observations in the context of historical storm-driven landslide events and open new avenues for investigating the nature and extent of landslides and their effects in southern Appalachia and similar environments.

Georgia, North Carolina, South Carolina, Tennessee↗

USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗