SEARCH · USGS Science

Results for “International Journal for Disaster Risk Reduction”

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,443 records · Page 44Linked to original sources

USGS2025

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2022

Potentiometric-surface and depth-to-water maps for spring 2022 were created for the Mississippi River Valley alluvial aquifer (MRVA) using groundwater-altitude data from 1,136 wells completed in the MRVA and from the altitude of the top of the water surface in area rivers from 160 streamgages. The potentiometric-surface and depth-to-water maps for 2022 were created to support investigations to characterize the MRVA as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the MRVA for spring 2022 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 340 feet (ft) above the North American Vertical Datum of 1988. The regional direction of groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient was into the depression, and near rivers, where the groundwater gradient can be from aquifer to the river or from the river into the aquifer. There are large depressions in the potentiometric-surface map in the lower one-half of the Cache region and in much of the Grand Prairie and Delta regions. Depth to water in the MRVA, spring 2022, by well ranged from 5.00 ft above land surface to 145.66 ft below land surface.

Mississippi River Valley alluvial aquifer↗

USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗

USGS2025

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗

USGS2026

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗

USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗

USGS1995

The distribution of bromide in water in the Floridan aquifer system, Duval County, northeastern Florida

Although Duval County, Florida, has ample ground-water resources for public supply, the potential exists for a problem with excessive disinfectant by-products. These disinfectant by-products result from the treatment of raw water containing low concentrations of bromide and naturally occurring organic compounds. Because of this potential problem, the relation of bromide concentrations to aquifer tapped, well location and depth, and chemical characteristics of water in the Floridan aquifer system underlying Duval County were studied to determine if these relations could be applied to delineate water with low-bromide concentrations for future supplies. In 1992, water samples from 106 wells that tap the Floridan aquifer system were analyzed for bromide and major dissolved constituents. A comparison of bromide concentrations from the 1992 sampling with data from earlier studies (1979-80) indicates that higher bromide concentrations were detected during the earlier studies. The difference between the old and new data is probably because of a change in analytical methodology in the analysis of samples. Bromide concentrations exceeded the detection limit (0.10 milligrams per liter) in water from 28 of the 106 wells (26 percent) sampled in 1992. The maximum concentration was 0.56 milligrams per liter. There were no relations between bromide and major dissolved constituents, well depth, or aquifer tapped that would be useful for determining bromide concentrations. Areal patterns of bromide occurrence are not clearly defined, but areas with relatively high bromide concentrations tend to be located in a triangular area near the community of Sunbeam, Florida, and along the St. Johns River throughout Duval County.

Florida↗

USGS2016

Sedimentary exhalative (sedex) zinc-lead-silver deposit model

This report draws on previous syntheses and basic research studies of sedimentary exhalative (sedex) deposits to arrive at the defining criteria, both descriptive and genetic, for sedex-type deposits. Studies of the tectonic, sedimentary, and fluid evolution of modern and ancient sedimentary basins have also been used to select defining criteria. The focus here is on the geologic characteristics of sedex deposit-hosting basins that contain greater than 10 million metric tons of zinc and lead. The enormous size of sedex deposits strongly suggests that basin-scale geologic processes are involved in their formation. It follows that mass balance constraints of basinal processes can provide a conceptual underpinning for the evaluation of potential ore-forming mechanisms and the identification of geologic indicators for ore potential in specific sedimentary basins. Empirical data and a genetic understanding of the physicochemical, geologic, and mass balance conditions required for each of these elements are used to establish a hierarchy of quantifiable geologic criteria that can be used in U.S. Geological Survey national assessments. In addition, this report also provides a comprehensive evaluation of environmental considerations associated with the mining of sedex deposits.

Scientific Investigations Report↗

USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗

USGS2024

Parameter ESTimation with the Gauss–Levenberg–Marquardt algorithm: An intuitive guide

In this paper, we review the derivation of the Gauss–Levenberg–Marquardt (GLM) algorithm and its extension to ensemble parameter estimation. We explore the use of graphical methods to provide insights into how the algorithm works in practice and discuss the implications of both algorithm tuning parameters and objective function construction in performance. Some insights include understanding the control of both parameter trajectory and step size for GLM as a function of tuning parameters. Furthermore, for the iterative Ensemble Smoother (iES), we discuss the importance of noise on observations and show how iES can cope with non-unique outcomes based on objective function construction. These insights are valuable for modelers using PEST, PEST++, or similar parameter estimation tools.

Groundwater↗

USGS2024

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper↗

USGS1990

Floods of February 1989 in Tennessee

Rainfall amounts of over 5 inches the night of February 13 and the morning of February 14, 1989, caused flooding in areas of Middle and West Tennessee. The towns of Lebanon in Middle Tennessee and Obion in West Tennessee were most severely affected. Most of the business district in Lebanon and many residential areas in Obion were flooded. Recurrence intervals for 24-hour rainfall totals were as high as 25 years at some sites but most peak discharges had recurrence intervals of less than 10 years. Rainfall amounts for the period February 13-20, 1989, peak stages and discharges for this flood, the peak of record, and a list of discharge measurements made during the flood are documented in this report.

Tennessee↗

USGS1976

A comparison of electronic particle counting and pipette techniques in routine mud analysis

Electronic particle counting and pipette size-analysis techniques are compared for their relative effectiveness in large-scale routine mud analysis. Optimum laboratory procedures for conducting electronic analyses with a Model TA Coulter Counter are presented. Comparative electronic and pipette analyses conducted on the same samples resulted in equivalent precision for both techniques. Results of electronic analyses consistently showed size distributions significantly coarser than those obtained by pipette analyses. The laboratory efficiency of the relatively rapid electronic technique is far superior to that of the pipette technique, although more extensive precautionary measures are required to maintain a high level of analytical accuracy and precision. The quality of electronic analyses largely depends upon the effectiveness of standardized operational procedures established within the laboratory.

Open-File Report↗

USGS2026

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗

USGS2023

Stratigraphic notes—Volume 1, 2022

This is the first volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 . This first volume ("Stratigraphic notes—Volume 1, 2022") includes articles that provide guidance for those who wish to submit papers to “Stratigraphic Notes,” as well as information on how to make your manuscripts compliant for geologic names reviews and how to organize your paper’s content to facilitate inclusion of new or revised names in Geolex. This volume also includes some specific guidance on conducting geologic names reviews of geologic and hydrogeologic reports.

Professional Paper↗

USGS2025

The impact of burial diagenesis on soil-formed minerals in paleosols using stable isotopes of phyllosilicates and carbonate clumped isotopes

To understand the effects of burial diagenesis on the stable isotope geochemistry of soil-formed clay and carbonate minerals in paleosols, samples were collected from seven cores, spanning middle- to upper-Pennsylvanian strata of the Illinois Basin, with varied maximum burial depths of 1–3 km. Mixed-layer illite-smectite and kaolinite mixtures give δ 2 H and δ 18 O values of −83 ‰ to −36 ‰ and 11.9 ‰ to 21.1 ‰ (VSMOW), respectively. After carbonates were screened petrographically for diagenetic textures using transmitted light and cathodoluminescence, measured clumped isotope Δ 47 values range from 0.504 to 0.563 ‰ (I-CDES). Resulting mineral formation temperatures for phyllosilicate mineral mixtures are 28 to 66 °C (mean = 47 °C), whereas T(Δ 47 ) estimates for calcites are 36 to 61 °C (mean = 45 °C). Calculated δ 18 O water values from which phyllosilicate minerals and calcites precipitated under isotopic equilibrium ranges from −7.1 to −1.2 ‰ and − 1.4 to +4.9 ‰, respectively. Closed and open-system phyllosilicate-fluid exchange modeling indicates that phyllosilicate alteration occurred in the presence of a low temperature brine or meteoric water and is interpreted to occur in a layer-by-layer illitization transformation. Due to the lack of diagenetic textures and positively correlated T(Δ 47 ) and δ 18 O water , calcites are interpreted to have undergone solid-state bond reordering. Despite low to moderate temperatures (<125 °C) and varying depths of shallow burial (1–3 km), solid-state transformation of phyllosilicates and calcites indicates paleosols had prolonged exposure to burial conditions which has implications for the use of paleosol minerals for paleoenvironmental reconstructions.

Illinois, Indiana, Kentucky↗

USGS2025

Comparisons of shoreline positions from satellite-derived and traditional field- and remote-sensing techniques

Satellite-derived shorelines (SDS) have the potential to help researchers answer critical coastal science questions and support work to predict coastal change by filling in the spatial and temporal gaps present in current field-based and remote-sensing data collection methods. The U.S. Geological Survey conducted comparison analyses of traditionally sourced shorelines and SDS in diverse coastal landscapes to determine how SDS could be used in ongoing and future work across varied coastal environments and provided some initial findings that could be used for implementation. Using CoastSeg, a browser-based program for SDS detection and mapping, SDS for the period 1984–2023 for multiple locations across the United States were compared to shoreline positions from traditionally sourced shoreline data. In this report, the authors present these comparisons alongside lessons learned and challenges encountered when building SDS workflows in different coastal locations. Results show that individual SDS have larger uncertainty and yet produced similar linear trends to sparser, traditionally sourced shoreline data; because SDS methods provide orders of magnitude more data than traditional shoreline-detection methods, they can be used to evaluate shoreline behaviors. Refining average scalar slopes used in tidal corrections did not result in substantial decreases in uncertainty. Using lessons from this work to outline needs for regional implementation, initial setup time would be considerable, being on the order of weeks. However, once complete, shoreline detections and analyses are fast (on the order of minutes to hours) and achievable using a desktop computer.

Alaska, Florida, Massachusetts, Washington↗

USGS2024

Geophysical identification of potential groundwater – surface water interactions in the Bonita Peak Mining District, San Juan County, Colorado

This report includes the results and interpretation of geophysical investigations to aid in the location, identification, and conceptual site model (CSM) development of the near surface groundwater and surface water hydrology of the Upper Cement Creek and California Gulch of the Animas Rivers within the Bonita Peak Mining Districts (BPMD). Mine influenced water (MIW) from numerous mine adits within the BPMD has altered the geochemistry and ecology of these ground and surface waters. Many studies of varying scales have investigated the geology, hydrogeology, geochemistry, and aqueous chemistry of these adits and impacted areas. Field investigations using non-invasive geophysical methods were conducted in the summer of 2019 to aid in identifying potential groundwater inputs to these streams and further the CSM of the near surface within and along these streams. Two field deployments were completed in August and September of 2019. In August 2019, Phase 1 was completed using multifrequency electromagnetic induction (EMI) and thermal infrared imaging (TIR) to measure the bulk electrical conductivity, magnetic susceptibility, and the temperature response at known and suspected groundwater seeps into the Upper Cement Creek and California Gulch stretches of these alpine streams. Results and interpretation from Phase 1 were used to co-locate stretches of these streams for further detailed fiber optic distributed temperature system (FODTS) monitoring of the temperature dynamics for submeter identification of potential groundwater seeps during Phase 2. Phase 2 included additional EMI measurements, FODTS measurements, and repeated some Phase 1 measurements for improved data quality. All these data were used for the combined qualitative interpretation identifying seep locations. The data from both Phase 1 and Phase 2 are available in the USGS data release at: https://www.sciencebase.gov/catalog/item/5f18a30082cef313ed84879b. These data include detailed near-surface electrical conductivity and magnetic susceptibility mapping, discrete TIR imaging, and several kilometers of FODTS data. These data were georeferenced and used to identify potential groundwater seeps into these streams. The results indicate a complex near surface geophysical and hydrogeologic regime. Seeps are classified and spatially correlated to the EMI, TIR, and geochemical sampling results. This report is a thorough compilation of the geophysical field efforts in the summer of 2019, the data processing, and interpretations of these data. These data alone should not be used for final interpretations of the complex hydrogeology within these alpine systems. Rather, these data should be combined with other known data from various investigations to further refine the CSM and guide decisions for future investigations, watershed protective measures, and other site decisions. Geographical information system (GIS) layers of these results are available to project managers for such integration and combined interpretations.

Colorado↗

USGS2026

Increasing artificial light at night enhances salmonid predator foraging efficiency in an urbanized lake

Artificial light at night (ALAN) poses a threat to ecosystems globally. It includes both direct and indirect light, or skyglow, which occurs when ALAN scatters in the atmosphere, extending beyond its original source. We analyzed ALAN trends in Lake Washington, WA, from 2014 to 2023 using Visible Infrared Imaging Radiometer Suite (VIIRS) nighttime light measurements, evaluated the relationship between in situ and satellite measurements, and modeled predator search volumes for a juvenile salmon predator, cutthroat trout ( Oncorhynchus clarki ), as a function of ambient light and turbidity conditions. Open water regions experienced significant increases in ALAN, while nearshore areas primarily showed no or negative trends, revealing the role of skyglow in shaping open water light environments. Using a visual foraging model for cutthroat trout, we found that juvenile salmon at the shallow southern pelagic site experienced light 28 times brighter, resulting in a 168% greater predation vulnerability than those at the northern site. In the Ship Canal, a narrow corridor for outmigrating salmon, predator search volumes were 249% higher than at the southern site. These contrasts in predation vulnerability demonstrate how local conditions influence predator–prey dynamics and provide critical insight for targeting mitigation of both nearshore and distant light sources.

Washington↗