USGS Science⌕ Search

SEARCH · USGS Science

Results for “Annual Review of Environment and Resources”

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

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

1,696 records · Page 69Linked to original sources
USGS2025

Integrating marine historical ecology into management of Alaska’s Pacific cod fishery for climate readiness

The Pacific cod ( Gadus macrocephalus ) fishery was closed in 2020 after a rapid decline in biomass caused by the marine heat waves of 2014–2019. Pacific cod are exceptionally thermally sensitive and management of this fishery is now challenged by increasingly unpredictable climate conditions. Fisheries monitoring is critical for climate readiness, but short-term monitoring data may be inadequate for recognizing and anticipating change under rapid climate changes. We propose an interdisciplinary, marine historical ecology framework that looks to long-term records (local and traditional knowledge, history, archaeology, and paleoclimatology) to capture a long range of ecological variability and provide historical context for management. In order to connect to contemporary fisheries management, this framework must be built on a common vocabulary and an understanding of the key metrics used in fisheries stock assessments. Here, we propose metrics derived from Pacific cod stock assessment and synthesize information relevant to understanding the effects of past warming periods on cod populations across the Gulf of Alaska and Bering Sea. This case study provides a framework for thinking about how to use these historical records in the context of fisheries management under rapidly changing climate conditions.

Alaska↗
USGS2026

Chronic exposure to waterborne nickel significantly reduced growth of juvenile crayfish (Faxonius virilis)

Crayfish are critical functional components of aquatic ecosystems. Previous research has documented adverse effects of mineral extraction on crayfish. Here, we characterize potential risks of mining-derived waterborne nickel (Ni) to crayfish by documenting the effects of dissolved Ni on growth and food consumption of juvenile virile crayfish ( Faxonius virilis) in a 28-day chronic laboratory exposure. Nominal Ni concentrations ranged from 31.25 to 500 micrograms per liter (µg/L; pH = 7.96 ± 0.20, hardness = 150 ± 1 milligrams per liter as calcium carbonate). Crayfish survival, carapace length, and wet weight were measured. After 28 days of exposure, a 24-h feeding trial was performed to determine differences in food consumption. During the growth trial, 99% of crayfish survived. Change in wet weight and final wet weight were the most sensitive endpoints, with 20% effect concentrations of 24.8 and 22.6 µg/L Ni, respectively. Crayfish exposed to an average of 438 µg/L Ni consumed 41% less, and weighed 65.1% less, than control crayfish. These results suggest chronic, sublethal exposure to waterborne Ni may have negative effects on crayfish growth. Reduced growth and consumption rates in crayfish could have wide-ranging consequences throughout aquatic ecosystems since crayfish are consumers, prey, keystone trophic regulators, and ecosystem engineers. Finally, these results could inform bioenergetics and may be coupled with population models to predict potential changes in population sizes of native and invasive crayfishes.

Ecotoxicology↗
USGS2024

California State Waters Map Series—Benthic habitat characterization in the region offshore Humboldt Bay, California

Coastal and Marine Ecological Classification Standard (CMECS) geoform, substrate, and biotic component geographic information system (GIS) products were developed for the California State Waters of northern California in the region offshore of Humboldt Bay. The study was motivated by interest in development of offshore wind-energy capacity and infrastructure in Federal waters offshore. This project, carried out by the U.S. Geological Survey (USGS), resulted in four data releases for individual map blocks that are part of the “California State Waters Map Series”: (1) Offshore of Arcata, (2) Offshore of Eureka, (3) Offshore of the Eel River, and (4) Offshore of Cape Mendocino. The study area consists of 436 square kilometers of multibeam echo sounder (MBES) data acquired by Fugro Pelagos, Inc., in 2007. Towed camera-sled video was acquired in 2009 and 2010 to supervise the classification of the MBES data into habitats, and single channel sparker data were collected to calculate sediment thickness above the transgressive unconformity. Using video observations of habitat as ground truth, derivatives of the MBES data were classified into 3 seafloor character types (hard-rugose, hard-flat, and soft-flat), 26 induration-slope-depth groups, and 15 geoforms. The study area substrate is predominantly soft-flat sediment (mud and fine sand) covering 73.6 percent of the area. Hard-flat substrate areas, predominantly coarse sediment in scour depressions, cover 5.4 percent of the study area. The hard-rugose substrate areas are primarily outcrops of layered sedimentary bedrock and constitute 20.9 percent of the study area. Fifteen geoforms were identified in the analysis. The predominant geoforms mirror the seafloor character results, shelf geoforms, rock outcrop geoforms, and scour depression geoforms. Rock and scour areas are restricted to the southern portion of the study area off Cape Mendocino where uplift has exposed bedrock. On the flat shelf area post-transgressive sediment varies in thickness from 1.7 meters (m) nearshore to 28.1 m offshore.

California↗
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↗
USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗
USGS2024

How does the onset of offset influence geologic slip rates?

Geologic slip rates are typically based on the displacement accrued by a geomorphic or stratigraphic feature and the age of the offset feature. Because slip rates are commonly calculated by dividing the displacement of a faulted marker by its age, they contain two open time intervals: the elapsed time between the age of an offset feature and the age of the earthquake that displaced the feature, and the time between the present‐day and the most recent earthquake. Here, we explore the influence of including unconstrained open intervals in geologic slip rate calculations. We test the degree to which these open intervals affect geologic slip rates and their uncertainties, and we find that their influence depends primarily on mean earthquake recurrence intervals (RIs). Slip rates on faults with longer RIs, such as the Wasatch fault, can be greatly influenced by an increase of up to 20% when accounting for open intervals. In contrast, slip rates on faults with shorter RIs, such as the San Andreas fault, are only slightly influenced by the assumption that slip rates calculated over open intervals approximate those calculated over closed intervals. Our analyses indicate that faults with moderate slip rates (∼0.2–5 mm/yr) are sensitive to both open interval effects themselves, as well as methods to quantify and account for these effects. We re‐evaluate how slip rates are calculated and defined in displacement–time space using published deformation records. We explore the utility of assigning a probability distribution to the initiation of offset of the oldest faulted feature and the timing of the most recent earthquake (MRE). We find that calculating geologic slip rates without using probability distributions that capture the timing of the MRE and the onset of offset of the oldest faulted feature, especially on slow‐to‐moderate slip rate faults, can lead to systematic underestimation of average geologic slip rates.

Seismological Research Letters↗
USGS2025

Cave records reveal recent origin of North America’s deepest canyon

We explore how and when Hells Canyon, North America’s deepest river gorge (~2,400 m deep), formed, addressing these fundamental questions first posed by W. Lindgren [ The Gold Belt of the Blue Mountains of Oregon (1901)]. Existing hypotheses about the canyon’s formation and timing of incision remain speculative due to a lack of direct constraints and geomorphic analysis in the canyon. Herein, we combine cosmogenic nuclide dating of cave-bound river deposits, river profile analysis, and numerical modeling to provide the first direct age constraints and systematic analysis of incision processes at work in Hells Canyon. Our study reveals a significant drainage capture triggered rapid incision at ~2.1 ± 1.0 Ma, establishing the Snake River’s modern route into the Columbia River system. The increased drainage area and subsequent increase in stream power resulted in the rapid incision of Hells Canyon and the formation of tributary knickpoints (KPs) that decrease in elevation away from the capture location. Cosmogenic dating of cave deposits indicates incision rates increased from ~0.01 to ~0.16 mm y −1 . Numerical modeling of the stream capture supports these observations, demonstrating how abrupt drainage area increase drives rapid river incision. Our findings from Hells Canyon provide a well-constrained example of how drainage capture can dramatically shape the evolution of a major river gorge.

Idaho, Oregon, Washington↗
USGS2025

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land↗
USGS2025

World minerals outlook—Cobalt, gallium, helium, lithium, magnesium, palladium, platinum, and titanium through 2029

Given the rapid expansion in the demand for mineral commodities that underpin worldwide economic growth and technological advancement, information regarding expected country-level mine production and production capacity is becoming increasingly important to industry stakeholders, end users, and policymakers. Production capacity can limit future supply, depending on how rapidly that capacity is able to expand. Current capacity can be evaluated on the basis of past production. Decreases to future capacity can be taken into account from announcements of planned shutdowns of mines or processing facilities, which are frequently publicized well in advance of such closures. Likewise, capacity expansions, which usually involve multiple stages—such as permitting, financing, and construction (all of which take time)—can also be estimated. As such, it is possible to evaluate midterm future capacity based on estimates of today’s capacities along with consideration of future investment plans. This World Minerals Outlook provides estimated capacities for cobalt, gallium, helium, lithium, magnesium, palladium, platinum, and titanium for 2025 through 2029. The results of this analysis indicate that two mineral commodities important to the manufacture of lithium-ion batteries—cobalt and lithium—are expected to have large capacity growth in the next few years, likely owing to expectations for increased demand for these batteries. For gallium, helium, palladium, and platinum, capacity is expected to remain stable or exhibit moderate growth. Still, these expected capacity levels are higher than current production, allowing for future production growth. The production capacity outlook is opaque for magnesium and titanium metal, which have a significant fraction of current production in nonmarket economies, such as China and Russia. Ultimately, though, where free market conditions prevail, full utilization of capacity potential for those commodities is likely to depend on supply deficits and prices that are above production costs.

Scientific Investigations Report↗
USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2024

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W–trending Central China Orogenic Belt has been overprinted by the N-S–trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau's weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W–trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene–Oligocene) and reworking of the deep lithosphere (Stage II, middle–late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

GSA Bulletin↗
USGS2025

Ultramafic float rocks at Jezero crater (Mars): Excavation of lower crustal rocks or mantle peridotites by impact cratering?

Based on observation and data from meteorites and in situ scientific missions, experiments as well as models, the Martian mantle is assumed to share some compositional and mineralogical affinity with the terrestrial mantle. However, there might be subtle differences like the Martian mantle being more ferroan. Yet, we do not have any direct analysis of a Martian mantle rock to confirm this assumption. NASA’s Perseverance rover found olivine-rich boulder-sized float rocks on the upper Jezero fan (Mars). These boulders have an ultramafic composition and their mineralogy is dominantly composed of Fo 73±3 olivine with high-Mg orthopyroxene, Cr-rich Ti-Fe oxides and minor plagioclase and high-Ca pyroxene. Microtextural and petrological analysis reveals that these minerals crystallized at equilibrium. In addition, these boulders are different from all the bedrocks analyzed by Perseverance along its traverse which are crustal igneous rocks and sediments. Comparing our data to Martian meteorites and available Mars bulk silicate models (BSM), we discuss that these boulders could represent primitive melts and/or lower crustal material, and we specifically hypothesize that they could be mantle peridotites. We propose that these putative mantle rocks could have been excavated by the succession of impacts from the shallow mantle or lower crust in the Isidis region where Jezero crater is located. These olivine-rich boulders could thereby constitute the first direct analysis of a Martian mantle rock.

Earth and Planetary Science Letters↗
USGS2026

Modeled groundwater and surface-water interactions surrounding Mobile Bay, Alabama, 2008–15

The U.S. Geological Survey, in cooperation with the Gulf Coast Ecosystem Restoration Council, has used MODFLOW 6 to develop a groundwater-flow model to simulate groundwater and surface-water interactions in the Mobile Bay, Alabama, area. The model results indicated that, on average, groundwater discharge near the coastline is equal to 2.5 percent of the surface water that flows into the bay. The model was also used to determine how changes in recharge, sea level, and groundwater pumping affect groundwater levels and discharge rates. The results indicate that more groundwater discharge occurred in the winter and spring when recharge was higher, sea level was lower, and groundwater pumping was lower than during the summer and fall. Additionally, the amount of emergent groundwater was closely related to sea level; when sea level was higher, there was more area with emergent groundwater. Furthermore, the depth of nonemergent groundwater was related to trends in recharge and pumping. During periods of increased recharge and reduced pumping, a significant portion of the model area exhibited a depth to the water table of less than 1 meter below the land surface. Conversely, during periods of decreased recharge and heightened pumping, much of the model area showed a depth to the water table ranging from 1 to 5 meters.

Alabama, Mississippi↗
USGS2024

Current status of the community sensor model standard for the generation of planetary digital terrain models

The creation of accurate elevation models (topography) from stereo images are critical for a large variety of geospatial activities, including the production of digital orthomosaics, change detection, landing site analysis, geologic mapping, rover traverse planning, and spectral analysis. The United Stated Geological Survey, Astrogeology Science Center, continues to transition the supported planetary sensor models to the Community Sensor Model (CSM) standard. This paper describes the current state of use for this photogrammetric standard, supported sensor model types, and qualitatively compares derived topography between SOCET SET and SOCET GXP ( ® BAE Systems) using HiRISE stereo images of Mars. Our transition to the CSM standard will ensure an uninterrupted capability to make these valuable products for Mars and many other extraterrestrial planets and moons.

Remote Sensing↗
USGS2024

Groundwater and surface-water interactions in the Lower Duwamish Waterway, Seattle, Washington

The U.S. Geological Survey (USGS), in cooperation with the Washington State Department of Ecology (Ecology), conducted a study to describe the current understanding of the regional groundwater system of the lower Duwamish River valley and groundwater and surface-water interactions in the lower Duwamish Waterway. The lower Duwamish Waterway is the final 5-mile (mi) reach of the Duwamish River before it empties into Elliott Bay in Puget Sound near Seattle, Washington. A nearshore site (hereinafter referred to as “Nearshore Site” to distinguish the particular site from general discussions of nearshore areas) along the western shoreline of the Duwamish River, about 1.5 mi upstream from the river mouth, was selected for focused groundwater data collection by USGS. Data loggers were deployed in seven groundwater wells and one stilling well in the Duwamish River to measure specific conductance, temperature, and depth at 15-minute intervals for a period of about 2 years. At the Nearshore Site during 2020–22, water levels in the shallow wells were 3–8 feet (ft) higher than water levels in the deep wells, providing evidence for a low-permeability layer between the shallow and deep aquifers in this area. The shallow wells had a pronounced seasonal variability, with high water levels in winter and low water levels in summer. Data from the deep wells showed far less seasonal variability, with slight increases in winter and a near-constant water level from spring to autumn. The deep wells had a strong hydraulic connection to the Duwamish River, as evidenced by the synchronous water-level variability during the tidal cycle, whereas the shallow wells had minimal to no tidal response. The potentiometric maps developed for the Nearshore Site and surrounding areas indicate large differences in groundwater-flow directions for the shallow and deep aquifers at low and high tides. For the shallow aquifer, flow is toward the lower Duwamish Waterway near the Nearshore Site, regardless of the tidal condition. For the deep aquifer, a potentiometric trough forms parallel to the shoreline during high tide, indicating that groundwater flow converges from the uplands to the west and the Duwamish River to the east. The geometry of the potentiometric surfaces between the nearshore-most well and the shoreline is complex and is further confounded by intermittent shoreline armoring and other buried infrastructure, which could serve as either a barrier or a conduit to flow. Groundwater and surface-water interactions in the lower Duwamish Waterway are inherently complex as a result of three overarching factors. First, water levels in the lower reaches of the Duwamish River vary daily by 11–16 ft because of tides from Puget Sound, which create large swings in the hydraulic gradient in the nearshore groundwater system. Second, the density and chemical composition of water in the Duwamish River change daily with the tides and seasonally, which constrains how river water entering the nearshore sediments interacts with discharging groundwater. Third, the nearshore subsurface and shoreline conditions are heterogenous because of extensive shoreline armoring over the past century, which governs the flow of groundwater and infiltrating river water. These unique features of groundwater and surface-water interactions in the lower Duwamish Waterway thus govern the transport of terrestrial contaminants to the lower Duwamish Waterway. Furthermore, the heterogenous aquifer properties in the lower Duwamish Waterway contribute to spatially and temporally dynamic contaminant-transport processes.

Washington↗