Search USGS⌕ Search

SEARCH · Search USGS

Search Search USGS

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.

At least 1,207 records · Page 67Linked to original sources
USGS2026

Luminescence dating of hydrothermal explosions in the Yellowstone Plateau volcanic field

Hydrothermal explosions are a significant geological hazard in some active volcanic systems; however, the timing and triggering mechanisms of these explosions are poorly constrained. This study applies luminescence dating techniques to hydrothermal explosion deposits in the Yellowstone Plateau volcanic field to constrain explosion chronologies and evaluate potential triggering mechanisms. We tested four luminescence dating techniques: K-feldspar post-infrared infrared stimulated luminescence (pIRIR225), quartz blue light optically stimulated luminescence (BLOSL), quartz blue thermoluminescence (BTL), and quartz red thermoluminescence (RTL). The pIRIR225 and RTL protocols produce consistent age estimates that agree with independent radiocarbon ages and with the timing of the Pinedale deglaciation. This study focuses on two craters, Mary Bay, along the northern shore of Yellowstone Lake, and Pocket Basin in Lower Geyser Basin. The mean pIRIR225 ages from Mary Bay deposits (11.99 ± 0.68 ka) agree with previous radiocarbon constraints. The mean pIRIR225 results from Pocket Basin deposits (13.44 ± 1.06 ka) suggest a history of explosion following Pinedale deglaciation, followed by recent hydrothermal alteration. Luminescence dating techniques are a promising tool for reconstructing the timing of hydrothermal explosions in the Late Pleistocene and Holocene, helping to constrain recurrence intervals of the largest hydrothermal systems, informing risk, and improving hazard assessments.

Idaho, Montana, Wyoming↗
USGS2026

Recent range expansion and documentation of a reproductive population of northern snakehead Channa argus (Cantor, 1842) in the Saint Francis River Drainage, Missouri

Northern snakehead Channa argus (Cantor, 1842) is an aquatic invasive fish species in the United States with first documented occurrence in the wild in the 2000s. Management efforts to control their populations in the eastern United States are ongoing. In the Mississippi River basin, limited resources have been allocated to control its distribution, after initial detection and rapid response in Arkansas were unsuccessful. Northern snakehead distribution in the Mississippi River basin was limited to Arkansas and Mississippi until 2019 when a single northern snakehead was detected on the southern border of Missouri in the Saint Francis River drainage, the furthest northern detection. Described here are additional northern snakehead detections following public reports and subsequent monitoring in the Mingo basin of the Saint Francis River drainage, a historical braided channel and floodplain habitat of the Mississippi River with intermittently flooded bottomland hardwood forests and wetlands, and other consistent aquatic habitats. These increasing captures document the recent range expansion of northern snakehead. Most of the 11,300 ha Mingo basin consists of Mingo National Wildlife Refuge and Duck Creek Conservation Area; these areas are protected aquatic ecosystems possessing sensitive species and serve as a potential example of prioritized areas for northern snakehead control efforts. Additionally, we highlight the significance of these detections in the Mingo basin which is connected via the Castor River water-control structure to the Upper Mississippi River and may facilitate further range expansion.

Missouri↗
USGS2026

Anuran occupancy of created wetlands within reforested legacy surface mines in Kentucky and West Virginia

Restoration of habitat is an important component of wildlife management. Surface coal mining has altered forest cover in the Central Appalachians, and most reclamation activities on previously mined lands result in non-native grasslands or shrub lands. The Forestry Reclamation Approach (FRA), which emphasizes non-native vegetation removal, soil decompaction, planting of native trees and shrubs, and, in some cases, wetland creation represents an alternative to grassland-based reclamation that may promote forest-associated wildlife on legacy surface mines. We used automated recording units to evaluate the response of anurans to created wetlands and the FRA in Kentucky and West Virginia, USA. We used a Bayesian community occupancy model to compare species occupancy and richness across a range of wetland sizes (6.6–252 m 2 ) and 3 site types: 1) younger FRA (1–6 yr), 2) older FRA (7–23 yr), and 3) unmined, mature forests (>100 yr). In addition, we evaluated factors influencing anuran detection probabilities. Mean occupancy and species richness in both states were highest in the younger FRA, suggesting anurans may rapidly colonize wetlands created on restored mines. In Kentucky, occupancy of several species (green frog [ Lithobates clamitans ], pickerel frog [ L. palustris ], and wood frog [ L. sylvaticus ]) was higher in larger created wetlands, but wetland size was not important in West Virginia. Daily minimum temperature positively influenced detection of anurans in both states, and time of day was also important for detection in West Virginia. Wetland creation and implementation of the FRA may be an important tool for wildlife managers aiming to improve anuran populations and diversity on legacy surface mines in the Central Appalachians.

Kentucky, West Virginia↗
USGS2026

Genetic biocontrol strategy considerations for mosquito control in the Pacific Island Countries and Territories

Mosquito-borne diseases pose an existential threat to the health, economies, and unique ecosystems of Pacific Island Countries and Territories (PICTs). The remoteness of these islands, combined with the presence of highly competent mosquito vectors, complicates disease surveillance and vector control efforts. In response, the PICTs have become a focal point for the development and application of Genetic BioControl (GBC) technologies designed to break vector-borne disease transmission cycles. However, the application of GBC tools in this region warrants careful consideration of its unique history, including a legacy of colonialism, and challenges associated with logistic hurdles. Through meaningful community engagement and authentic collaborations, drawing from local knowledge, and building local capacity, the sustainable, efficient, and effective deployment of GBC tools may be achieved.

Current Opinion in Insect Science↗
USGS2026

Hydrologic dynamics of ephemerally flooded playas in a dryland environment

Ephemerally flooded playas are common in the southwestern United States and globally in drylands. Often formed in closed basins, playas are depressions which inundate infrequently from local precipitation and streamflow produced near the playa or from upland areas. Few studies have quantified the hydrologic connectivity between upland catchments and playas using observations. Here, we used rain gauge-corrected precipitation from weather radar and water level measurements in 18 playas of the Chihuahuan Desert to identify precipitation thresholds leading to playa inundation over a 6.4-year period. Geospatial data sets on topography, soil properties, and vegetation cover were employed to determine the controls on inundation. Only 9.4% of all precipitation events above 1 mm led to inundation, with 69.8% of all inundations occurring during the North American monsoon (NAM, July-September). Mean and standard deviations (Std) of runoff ratios at all playas were 2.74 ± 4.08% and 3.29 ± 5.19% for annual and NAM periods. At the annual scale, playa inundation occurred when mean precipitation thresholds of 18.3 ± 7.5 mm (event total) and 12.0 ± 4.5 mm/hr (60-min intensity) were exceeded. Across all playas, inundation occurrence and volume were related most strongly to precipitation metrics and catchment area, with secondary controls of soil and terrain properties. The explanatory power of the derived regressions describing the inundation response across the playas were significantly improved when considering their geological origin. As a result, the inundation response classification system could be applied to ephemeral playas in other arid and semiarid landscapes.

New Mexico↗
USGS2026

Coral reef protection may help avert risks to people, property, and economic activity caused by projected reef degradation

Degradation of coral reefs over the past several decades has caused regional-scale erosion of the shallow seafloor that serves as a protective barrier against coastal hazards along southeast Florida, USA. How future change in coral reefs may affect coastal flooding, however, has been less attended than other factors contributing to increasing risks such as sea-level rise and more intense storms. Here, the increased flooding hazard faced by Florida's coastal communities from the projected future degradation of its adjacent coral reefs is evaluated through oceanographic, coastal engineering, habitat, geospatial, and socioeconomic modeling. Risk-based valuation approaches were followed to map flood zones at 10-m 2 resolution along 430 km of Florida's reef-lined coast for the current and projected future coral reef conditions. The projected degradation of Florida's coral reefs can increase annual flooding to more than 8.77 km 2 of land and 4,980 km of roads, affecting more than 7,315 people, $412.5 million in damages to 1,400 buildings, and economic disruption of $438.1 million annually (2024 US dollars). The degradation of Florida's coral reefs would increase the annual risk to people and structures by more than 42% and 47%, respectively, but is spatially variable due to the heterogeneous alongshore nature and distribution of the reefs and communities: the increased risk exceeds $1 million/km annually to more than 17% of the coastline but also disproportionately would affect vulnerable populations. These results help identify areas where coral reef protection could help reduce the projected increased storm flooding risk to Florida's coastal communities.

Florida↗
USGS2026

Compounding of 100-year coastal floods by rainfall in an urban environment

Coastal and pluvial flooding are both becoming more prevalent and severe due to climate change and urbanization in floodplains. The co-occurrence of these flood drivers is generally assumed to exacerbate the resulting flood impacts, a result referred to as compound flooding. However, few observational or modeling studies have investigated the circumstances under which this occurs. Here, we study the impacts of these combined flood drivers and evaluate the implicit hypothesis of official flood maps, which is that rainfall has a negligible impact on the flood depth and flooded area due to a 100 year coastal flood. A coastal system model, configured to capture coastal and pluvial flood drivers, is used. We evaluate the flooding for different urban landform types, including coastal landfill (human-made land), convergent areas (topographic depressions) and other urban terrain, within a model domain covering the Jamaica Bay watershed of New York City. A scenario-based strategy is adopted with a 100 year coastal flood as a control simulation, to which we add a set of realistic scenarios of rainfall data from historical tropical cyclones. We also apply a joint probability analysis framework with historical data to evaluate the probability of these compound coastal-pluvial scenarios. Results reveal cases where the pluvial driver compounds the coastal flood through expansion of the flood zone, with a 17% chance of rainfall increasing the flood area by 6%–38%, and a 5% chance of an increase of 61%–73%. It is rare that floods are significantly deepened but when deepening occurs, it is more common for the convergent zone than for the coastal landfill. These findings quantitatively assess the potential of the pluvial driver to exacerbate flooding, which may influence emergency management strategies such as evacuation plans, shelter arrangements, and related preparedness measures.

New York↗
USGS2026

Gerrard strain Rainbow Trout are not piscivorous in a shallow, polymictic reservoir

Objective Georgetown Lake (Montana, USA) is renowned for its multispecies salmonid fishery. Although many anglers target kokanee Oncorhynchus nerka in Georgetown Lake, the body length of kokanee has typically been considered unsatisfactory. To reduce the density of kokanee and increase the average size, Montana Fish, Wildlife & Parks began stocking the piscivorous Gerrard strain of Rainbow Trout O. mykiss (hereafter “Gerrard”) in 2015 to consume kokanee. In 2020, preliminary evidence suggested that the average length of kokanee greater than 200 mm increased by almost 40 mm. Methods To assess the efficacy of biocontrol through the introduction of a piscivore to predate on kokanee, we examined the diet composition of all potential predators (three strains of Rainbow Trout and Brook Trout Salvelinus fontinalis ) to determine the amount of predation on kokanee and to understand the feeding ecology of all potential predators. Results There was little evidence of piscivory in any predator and no evidence of Gerrards consuming kokanee. Gerrards exhibited a generalist feeding strategy, and there was dietary overlap and no difference in trophic position among Gerrards and other predators. Conclusions From a biological-control perspective, there is no benefit in continuing to stock Gerrards to reduce kokanee densities in Georgetown Lake. These findings underscore the challenge of predicting predator–prey relationships, as diet plasticity can lead to unexpected dynamics based on resource availability and habitat structure within aquatic ecosystems, highlighting the importance of evaluating management interventions, such as biocontrols, in different ecosystems.

Montana↗
USGS2026

An integrated mudstone facies classification scheme and revised interpretation of the sedimentologic processes driving carbon burial in the Cenomanian–Turonian Greenhorn Formation, Colorado, U.S.A.

Standardizing facies descriptions has proven key to integrating interpretations of depositional processes and environments from sedimentologic observations with geochemistry data for mudstone lithologies. Because of their fine-grained nature, high degree of compaction, and heterogeneous composition, standardizing methods for mudstone descriptions has proven difficult, but it is critical to formulating meaningful interpretations of the processes that govern the accumulation of organic-rich lithologies and their role in both petroleum systems and the global carbon cycle. In this study, we have developed a modified facies classification scheme for mudstone lithologies that incorporates sedimentologic and compositional observation at the hand-sample and thin-section scales with geochemical measurements, including bulk organic and inorganic geochemistry, to characterize these rocks and their variability more completely for improved interpretations of depositional environments during a low-order sea-level transgression. The facies described in this study are of the Cenomanian–Turonian Greenhorn Formation in the USGS #1 Portland Core drilled in Fremont County, Colorado. Strata of the Greenhorn Formation span Oceanic Anoxic Event 2 (OAE-2) and the preceding interval. Lithologies range from organic-rich argillaceous mudstones with varied sedimentary structures to organic-lean, highly bioturbated limestones. Six facies were identified, each differentiated by varied sedimentary structures and geochemical composition. These facies occur in a predictable stratigraphic stacking pattern that represents a low-order sea-level transgression with interpreted depositional environments ranging from terrigenous-dominated pro-delta and muddy continental shelf at the base of the interval to pelagic offshore marine at the top of the Greenhorn Formation. Though the facies are consistent with previous interpretations of depositional environments at this locale in the Cretaceous Western Interior Seaway during the Greenhorn cyclothem, the sedimentary processes governing the accumulation of organic-rich strata that have defined this interval are significantly revised. Variability in the proximity and intensity of bottom currents driven by storms and geostrophic flows were key to the accumulation of each facies, with significant sediment transport occurring even through deposition in the most oxygen-depleted bottom waters. The methodology and interpretations provided here are now being employed to basin-scale predictions of organic enrichment utilizing calibrated petrophysical methods. The approach and results from this study improve understanding of how organic and inorganic carbon was sequestered during perturbations to the global carbon cycle associated with events such as OAE-2.

Colorado↗
USGS2026

Monitoring recreation on federally managed lands and waters—Aspects of visitor use

Federally managed public lands and waters receive about 1 billion recreational visits each year. Data on these visitors can aid in guiding policy decisions, managing resources effectively, and communicating the economic contributions of lands and waters. This report explores the methods used by agencies to collect data on aspects of recreational visitor use to Federal lands and waters (apart from visitation numbers, which are the focus of a companion publication). Aspects of recreational visitor use include visitor demographics, recreational activity participation, visitor satisfaction, visitor attitudes and experiences, trip characteristics, and economic contributions. We review practices used to understand aspects of visitor use across seven Federal agencies, revealing similarities such as the use of visitor intercept surveys and coverage of similar topic area, and differences in how survey programs are operationalized and how specific questions on visitor surveys are worded. We also evaluate emerging technologies, such as geolocated social media and mobile device location data, for their potential to aid in understanding aspects of visitor use. This report concludes with potential opportunities to enhance data collection and coordination, ensuring cost-effective data collection and informed decision-making.

SocArXiv↗
USGS2026

Conducting feasibility assessments of potential conservation reintroductions: A case study with the imperiled foothill yellow-legged frog, Rana boylii

Conservation translocations are an increasingly common and often necessary component of recovering species that have become extirpated from portions of their range. Understanding and ameliorating potential threats that reduce the likelihood of successful population establishment at recipient sites is a key component of successful translocation planning. We examined multiple potential threats, including pathogens, contaminants, and invasive species, as well as habitat suitability and food resources, to assess the feasibility of reintroducing threatened, stream-obligate foothill yellow-legged frogs, Rana boylii , to Pinnacles National Park. Foothill yellow-legged frogs were extirpated from this protected area more than half a century ago. Although invasive species, disease, contaminants, food resources, and water temperatures are unlikely to inhibit foothill yellow-legged frog population establishment, potential recipient streams at Pinnacles National Park had shorter hydroperiods and much higher canopy cover than reference streams with extant foothill yellow-legged frog populations. Although the exact cause of extirpation of foothill yellow-legged frogs from Pinnacles National Park is unknown, translocations of foothill yellow-legged frogs to the park are more likely to succeed if riparian canopy cover is reduced and stream hydroperiods increased to better match those at nearby populations. Thoroughly understanding the threats to and characteristics of potential recipient sites could improve the likelihood of success of translocation outcomes in natural areas.

California↗
USGS2026

A review and synthesis of post-wildfire shifts in hydrologic processes and streamflow generation mechanisms

Critical water supply watersheds in the western United States (WUS) are impacted by wildfires, with potential negative effects on water quality and quantity. Scientific understanding is currently insufficient to deliver estimates of wildfire consequences for water quantity that are regionally accurate. Regional variability in the directionality and magnitude of post-wildfire shifts in streamflow generation fuels uncertainty in estimates of wildfire effects on water supply. In this work we provide a narrative review of wildfire effects on hydrologic processes and the resulting changes in streamflow generation mechanisms with a focus on the WUS, incorporating other global regions when pertinent. A conceptual model summary of wildfire effects on streamflow generation emphasizes: (1) precipitation seasonality, (2) synchrony of precipitation and potential evapotranspiration, (3) net shifts in interception, evaporation, and transpiration relative to total annual precipitation, (4) vegetation changes, including compensatory uptake and type conversion, (5) degree of overlap in rainfall rates and infiltration, (6) fire extent and severity, (7) burn scar positioning (e.g. in headwaters or proximal to watershed outlet), (8) scale-dependent groundwater leakage, (9) near-surface water storage reduction, and (10) soil to groundwater connectivity. Ongoing gaps and challenges include separating the influences of precipitation variability, water withdrawals, and post-fire land management; compound and overlapping disturbances; and lack of pre-fire data. Notable future opportunities include: harnessing ever-improving gridded and remotely sensed precipitation and fire-effects data; linking geophysical, isotopic tracer, and geochemical signatures to diagnose hydrologic changes; leveraging physically based and data-driven model advancements; and analyzing streamflow generation recovery trajectories across diverse watersheds.

western United States↗
USGS2026

Low-frequency earthquakes track the motion of a captured slab fragment

Accurate tectonic models are essential for assessing seismic hazard and fault interactions. However, the plate configuration at the complex Mendocino triple junction, where the San Andreas Fault and the Cascadia subduction zone meet, remains uncertain. We analyzed fault slip associated with a recently identified zone of tectonic tremor and low-frequency earthquakes (LFEs) near the southern edge of the subducting Gorda slab. Based on tidal sensitivity and P-wave first motions, we show that the LFEs are generated by dipping, strike-slip motion. This suggests that a former Farallon slab fragment, now captured by the Pacific plate, is translating northward beneath westernmost North America. This geometry effectively extends the slab interface fault, challenging prevailing interpretations of slab window formation and creating a potential unaccounted earthquake hazard in this region.

Science↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗
USGS2026

Groundwater tracing used to delineate recharge areas and map karst groundwater pathways for subterranean streams at Oregon Caves National Monument and Preserve

Oregon Caves National Monument and Preserve in southwestern Oregon is a 4,554-​acre area managed by the National Park Service that is home to several cave systems, including Oregon Caves, which is the longest cave in Oregon, with 3.03 miles of mapped passages. Because of the interconnected nature of karst hydrologic systems, it is critical to understand the areas that can influence water quality and quantity in karst environments. Toward this goal, dye tracing was conducted by the U.S. Geological Survey from 2021 to 2024 to better understand the pathways that karst groundwater follows at Oregon Caves National Monument and Preserve and to delineate recharge areas for two caves, Oregon Caves and Cave Next Door. During the project, eight dye injections were conducted, delineating a 0.51-​square-​mile recharge area for Oregon Caves and a 0.69-​square-​mile recharge area for Cave Next Door. Additionally, the study helped to identify three resurgences associated with Oregon Caves that were previously unknown and showed that the recharge areas for the two caves were distinct from one another. The dye traces also illuminated some unique recharge characteristics of the karst at Oregon Caves, including a high variance in karst groundwater velocities, retention within the karst aquifers, and a significant diffuse-​flow component.

Oregon↗
USGS2026

Unveiling a legacy of fish introductions to mountain lakes using historical records and eDNA surveys in a National Park

Across the western United States, introductions of non-native fish into historically fishless mountain lakes have impacted native biota. Understanding the impacts of fish introductions is essential for conservation in Olympic National Park, a Biosphere Reserve. We reconstructed fish plantings using records dating back to 1930, followed by environmental DNA (eDNA) surveys to estimate the current distribution of fish and amphibians in 117 remote mountain lakes. We used Bayesian multiscale occupancy models to determine how lake attributes and planting history related to fish and amphibian occupancy. The most frequently detected species were Brook Trout, Rainbow Trout, Cascades Frog, and Northwestern Salamander. eDNA sampling revealed 52 lakes with amphibians only, 45 with fish and amphibians, 14 with fish only, and 6 unoccupied. Of the 53 lakes with planting records, 38 had fish eDNA detected. Fish eDNA was also detected in 21 lakes lacking planting records, which could reflect incomplete records, unauthorized plantings, and false positive detections. Of the three species planted, Cutthroat Trout had the highest failure rate and did not become established in 23 of 28 historically planted lakes. In a subset of 9 lakes sampled for up to 7 years, those with known fish and amphibian presence showed consistent eDNA detections over time. The number of times a lake was stocked was the best predictor of occupancy for Brook and Rainbow trout, while higher occupancy for Brook Trout was also associated with lower elevations, lower solar radiation, and larger lake area. We did not observe widespread negative associations between amphibian occupancy and fish presence, although there was a negative relationship between fish presence and Rough-skinned Newt and Long-toed Salamander occupancy. Cascades Frog occupancy showed no relationship to fish presence or lake traits. Our results suggest mechanisms of fish persistence over time and highlight areas where native amphibians are impacted by introduced fish. These results can guide management options like targeted fish removals that benefit native fauna while still supporting recreational fishing. More broadly, our work demonstrates the value of combining historical records with contemporary surveys and the utility of eDNA for broad-scale surveys of species distribution in remote wilderness areas.

Washington↗