SEARCH · USGS Science

Results for “Journal of Waterway, Port, Coastal, and Ocean Engineering”

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

USGS2024

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗

USGS2025

The role of alkali bicarbonate-sulfate brines in the genesis of carbonatite REE resources at the Bear Lodge Alkaline Complex, Wyoming

Rare-earth element (REE) resources in the Bear Lodge Alkaline Complex, Wyoming, are hosted in a variably leached carbonatite dike swarm spatially related to bodies of diatreme breccia. This study examines fluid inclusions in carbonatite dikes, peripheral fluorite breccias, and smoky quartz veins to reconstruct the physiochemical conditions of REE mineralization. Results reveal a multi-stage evolution of carbonatitic fluids: 1) high-temperature (330–432°C), CO 2 -rich magmatic fluids, 2) REE-enriched alkali bicarbonate-sulfate brines, and 3) low-temperature meteoric (117–182°C) or diluted magmatic fluids. Multiphase inclusions (stage-2) contain burbankite, nahcolite, and alkali sulfate daughter crystals, linking the brines to early burbankite [(Na,Ca) 3 (Sr,Ba,REE) 3 (CO 3 ) 5 ] mineralization in carbonatite. REE mobility has been linked to the stability of aqueous complexes with ligands such as Cl-, SO 4 2- , and HCO 3 - . However, light REEs were not initially transported significant distances, rather they crystallized as burbankite from brines or brine-melts associated with the carbonatite magma. Crystallization of burbankite removed Na from the system, leaving residual fluids with high K/Na ratios and promoting potassic fenitization throughout the complex. Fractionation of light and heavy REEs and involvement of low temperature (<300°C), high K/Na bicarbonate brines may have facilitated transport and formation of peripheral occurrences that are more enriched in heavy REEs at Bear Lodge.

Wyoming↗

USGS2026

Complex carbonate ore mineralogy in the Mountain Pass carbonatite rare earth element deposit, USA

Economic concentrations of rare earth element (REE) minerals are uncommon in the Earth’s crust, with most occurring in carbonatites. Unlike most igneous rocks composed of silicate minerals, carbonatites are dominated by carbonate minerals, some of which can incorporate significant light REEs (LREEs; La, Ce, Pr, Nd). Technological applications of REEs are numerous and they have been identified as some of the most critical mineral commodities to the global economy. The Mountain Pass carbonatite stock in the Mojave Desert of California is the most economically significant REE deposit in the USA and contains a few to tens of percent (by volume) of the carbonate REE ore mineral bastnäsite. Despite the economic significance of the Mountain Pass deposit, studies of its ore mineralogy are limited. Here we present new carbonate ore mineralogy data for a compositionally diverse suite of carbonatitic rocks from the Mountain Pass stock and related dikes. Whole-rock geochemical data are integrated with mineral-scale textural and chemical data obtained by scanning electron microscopy (SEM), electron probe microanalysis (EPMA), and microRaman spectroscopy. Our results document a complex spectrum of REE-bearing carbonate minerals and intermediate mixed-layer structures. Mineral species include bastnäsite [REE(CO 3 )F], hydroxylbastnäsite [REE(CO 3 )OH], parisite [Ca(REE) 2 (CO 3 ) 3 F 2 ], synchysite [Ca(REE)(CO 3 ) 2 F], röntgenite [Ca 2 (Ce,La) 3 (CO 3 ) 5 F 3 ], and sahamalite [(Mg,Fe 2+ )(REE) 2 (CO 3 ) 4 ]. Carbonate ore mineralogy is heterogeneous within and between samples, including at the intracrystal scale. Complexly zoned crystals exhibit as many as five to six different compositional domains and syntaxial intergrowths, commonly with the more Ca-rich varieties (parisite, synchysite) forming crystal rims that surround relict bastnäsite cores. We attribute the phenocryst variability to changes in the chemistry and temperature of primary carbonatite magmas and evolved/exsolved fluids. Cross-cutting vein textures of calcite, celestine and various REE carbonate minerals, interstitial bastnäsite crystallization, breccia blocks lined by fine-grained bastnäsite, and the presence of hydroxylbastnäsite and partially hydroxylated bastnäsite point to the role of secondary hydrothermal processes in REE mineralization. Fluorcarbonate mineral compositions demonstrate that La and Ce are more structurally abundant in bastnäsite, whereas the more Ca-rich species (parisite, synchysite) contain a greater proportion of REE heavier than Pr (Nd, Sm, Eu, Gd) and Y. Atomic ratios of Pr/(Nd + Pr) are likewise variable, with the highest average value for bastnäsite (0.25) compared to parisite (0.22) and sychysite (0.21). This finding has geometallurgical implications, given that current mining operations are focused on recovery of Nd and Pr for high field strength permanent magnets and the Nd/Pr ratios are a critical factor in ore processing and magnet manufacture.

California↗

USGS2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater-dependent ecosystems in areas with industrial land use are at risk of exposure to a PFAS chemicals. We investigated one such system with several known PFAS source areas, where high and low permeability sediments (glacial) coupled with groundwater-lake and groundwater/surface-water interactions created complex ‘source to seep’ dynamics. Using heat-tracing and chemical methods, numerous preferential groundwater discharge zones were identified and sampled across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, downgradient of Joint Base Cape Cod (JBCC). Surface-water and groundwater samples were analyzed for 40 PFAS compounds between March and October 2022. Samples were collected from groundwater seeps identified as preferential discharge points (PDPs), wells upgradient of the stream-wetland system, contributing flow-through kettle lakes, and along Quashnet River surface-waters. PFAS from sampled waters contained perfluorinated carboxylates (PFPea, PFHxA, PFNA), perfluorinated sulfonates (PFBS, PFPeS, PFHxS, PFOS), fluorotelomer sulfonates (6:2, 8:2 FtS), and perfluoroalkyl sulfonamides (PFOSA). Samples from PDPs and wells had measured PFAS concentrations ranging from non-detect to 4677 ng/L ng/L (mean = 418 ng/L, std. = 709 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to groundwater flowpath emergence at PDPs. Correlations ( p < 0.01) between deuterium excess, %PFAS precursors, and terminal PFAS compounds highlighted potential precursor transformations associated with lake-groundwater exchange along flowpaths sourcing PDPs. However, some seepages had higher total PFAS concentrations (>1000 ng/L) than upgradient kettle lakes despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. Results from these synoptic surveys address gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations and mass loading in preferential groundwater discharge zones.

Massachusetts↗

USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗

USGS2025

Tracing metal sources and groundwater flow paths in the Upper Animas River watershed using rare earth elements and stable isotopes

Groundwater flow paths and processes that govern metal mobility and transport are difficult to characterize in mountainous bedrock watersheds. Despite the difficulty in holistic characterization, conceptual understanding of subsurface hydrologic and geochemical processes is key to developing remediation plans for locations affected by acid mine drainage, such as the Upper Animas River watershed in southwestern Colorado, USA. Stable isotopes of water and rare earth elements were utilized to evaluate groundwater flow and metal sources within this complex catchment. Stable isotope samples collected from draining mine adits and springs display systematic spatial variation wherein sample sites at higher elevations have greater seasonal variability than sites at lower elevations. The Upper Cement Creek watershed, where multiple draining mines are present, displays the lowest seasonal variation in stable isotopic signatures, potentially indicating the presence of a large, well-mixed volume of groundwater storage or interbasin groundwater flow. Rare earth elements display statistically significant variation between different alteration styles in the catchment. Overprinting of regional propylitic alteration is evident based on enrichment of middle rare earth elements in acidic springs and mines that are not spatially associated with surficial exposures of acid generating alteration styles. Europium anomaly and middle rare earth enrichment signatures from two flooded mine tunnels on opposite sides of a watershed divide indicate connections to the same subsurface flooded mine workings.

Colorado↗

USGS1976

Corwin Springs Known Geothermal Resources Area, Park County, Montana

The Corwin Springs Known Geothermal Resources Area (KGRA) is contiguous to Yellowstone National Park along a part of the northern boundary of the park near Gardiner, Park County, Mont. The area contains two known sites of hot-spring activity--LaDuke Spring, 2.8 km southeast of the small resort community of Corwin Springs, and Bear Creek Spring, 2.6 km east of Gardiner. LaDuke Spring issues from brecciated quartzite and has a flow rate of 380 1/min, a surface water temperature of 65°C, a silica-geothermometer temperature of 68.7°C, and a Na-K-Ca geothermometer temperature of 76.8°C. Bear Creek Spring issues from limestone and has a flow rate of 4 1/min, a surface water temperature of 32°C, a silica temperature of 46.7°C, and a Na-K-Ca temperature of 87.2°C. The springs, which are actively depositing travertine, are on or near the trace of the Gardiner fault-a high-angle reverse fault which forms the southwestern boundary of the Beartooth uplift. Features that make the Corwin Springs KGRA potentially significant include: (1) Proximity of the Corwin Springs area to significant Pliocene and Pleistocene volcanism in Yellowstone National Park and related geothermal activity, (2) existence of a large potential geothermal reservoir in the northern extension of the Sepulcher Mountain graben, (3) localization of thermal activity along the Gardiner fault in the Corwin Springs area, (4) occurrence of negative gravity anomaly centered over the northern end of the Sepulcher Mountain graben, (5) existence of negative magnetic anomaly in the Corwin Springs area, (6) location of the area within the Intermountain seismic belt, (7) evidence of recent tectonism in the area, as demonstrated by Pleistocene and Holocene faulting, and (8) observed surface temperatures and estimated geochemical temperatures of the two known hot springs (features which do not indicate high subsurface temperatures).

Montana↗

USGS2025

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California↗

USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗

USGS2026

Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows

Groundwater pumping-induced reductions in streamflow (known as ‘streamflow depletion’) have been documented worldwide, but potential impacts of streamflow depletion on stream temperature are not well understood. Here, we use two types of models to identify potential impacts of pumping on stream temperature across the conterminous United States (CONUS) to determine which aspects of a stream's annual thermograph (thermal signatures) can be used to monitor and manage streamflow depletion impacts on stream temperature. We used long-term streamflow and stream temperature data from 30 streamgages across CONUS and surrogate models of streamflow depletion to analyse potential stream temperature impacts at each site. We compared two different stream temperature modelling approaches: (i) a process-based energy balance model and (ii) statistical regression models based on air temperature and stream discharge. We calculated a suite of thermal signatures under depleted and non-depleted conditions for each stream and found that maximum annual 7-day temperature and annual temperature range are potentially the most sensitive to streamflow depletion, with potential changes of at least 2°C at > 70% of the sites when using the process-based model. We also found that the regression-based models predicted much less sensitivity of stream temperature to streamflow depletion than the process-based model. This work provides an initial evaluation and sensitivity analysis of the potential impacts of streamflow depletion on stream temperature. We demonstrate that stream temperature may be most sensitive to pumping in streams with a high proportion of flow sourced from relatively cold groundwater inputs, and that regression-based stream temperature models may underpredict stream temperature changes caused by streamflow depletion.

conterminous United States↗

USGS2025

lasertram: A Python library for time resolved analysis of laser ablation inductively coupled plasma mass spectrometry data

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) data has a wide variety of uses in the geosciences for in-situ chemical analysis of complex natural materials. Improvements to instrument capabilities and operating software have drastically reduced the time required to generate large volumes of data relative to previous methodologies. Raw data from LA-ICP-MS, however, is in counts per unit time (typically counts per second), not elemental concentrations and converting these count ratesto concentrations requires additional processing. For complex materials where the ablated volume may contain a range of material compositions, a moderate amount of user input is also required if appropriate concentrations are to be accurately calculated. In geologic materials such as glasses and minerals that potentially have numerous heterogeneities (e.g., microlites or other inclusions) within them, this is typically determiningwhether the total ablation signal should be filtered to remove these heterogeneities. This necessitates that the LA-ICP-MS data processing pipeline is one that is not automated, but is also designed to enable rapid and efficient processing of large volumes of data. Here we introduce , a Python library for the time resolved analysis of LA-ICP-MS data. We outline its mathematical theory, code structure, and provide an example of how it can be used to provide the time resolved analysis necessitated by LA-ICP-MS data of complex geologic materials. Throughout the pipeline we show how metadata and data are incrementally added to the objects created such that virtually any aspect of an experiment may be interrogated and its quality assessed. We also show, that when combined with other Python libraries for building graphical user interfaces, it can be utilized outside of a pure scripting environment. can be found at https://doi.org/10.5066/P1DZUR3Z

Applied Computing and Geosciences↗

USGS2024

Vulnerability assessment of groundwater influenced ecosystems in the Northeastern United States

Groundwater-influenced ecosystems (GIEs) are increasingly vulnerable due to groundwater extraction, land-use practices, and climate change. These ecosystems receive groundwater inflow as a portion of their baseflow or water budget, which can maintain water levels, water temperature, and chemistry necessary to sustain the biodiversity that they support. In some systems (e.g., springs, seeps, fens), this connection with groundwater is central to the system’s integrity and persistence. Groundwater management decisions for human use often do not consider the ecological effects of those actions on GIEs. This disparity can be attributed, in part, to a lack of information regarding the physical relationships these systems have with the surrounding landscape and climate, which may influence the environmental conditions and associated biodiversity. We estimate the vulnerability of areas predicted to be highly suitable for the presence of GIEs based on watershed (U.S. Geological Survey Hydrologic Unit Code 12 watersheds: 24–100 km 2 ) and pixel (30 m × 30 m pixels) resolution in the Atlantic Highlands and Mixed Wood Plains EPA Level II Ecoregions in the northeastern United States. We represent vulnerability with variables describing adaptive capacity (topographic wetness index, hydric soil, physiographic diversity), exposure (climatic niche), and sensitivity (aquatic barriers, proportion urbanized or agriculture). Vulnerability scores indicate that ~26% of GIEs were within 30 m of areas with moderate vulnerability. Within these GIEs, climate exposure is an important contributor to vulnerability of 40% of the areas, followed by land use (19%, agriculture or urbanized). There are few areas predicted to be suitable for GIEs that are also predicted to be highly vulnerable, and of those, climate exposure is the most important contributor to their vulnerability. Persistence of GIEs in the northeastern United States may be challenged as changes in the amount and timing of precipitation and increasing air temperatures attributed to climate change affect the groundwater that sustains these systems.

Connecticut, Delaware, District of Columbia, Maine↗

USGS2024

Potential corrosivity of untreated groundwater in Louisiana

Corrosive groundwater can cause lead, copper, and other metals to leach from pipes and plumbing fixtures in water distribution systems. Metals, if ingested, could lead to serious health implications to the nearly 2.9 million people in Louisiana who obtain their drinking water from groundwater sources. Four indices—the Langelier Saturation Index (LSI), Ryznar Stability Index (RSI), Puckorius Scaling Index (PSI), and the Potential to Promote Galvanic Corrosion (PPGC)—in addition to an analysis which normalized the results from the existing indices, the Combined Index (CI), were used to assess the corrosivity of groundwater in Louisiana and identify areas within eight major aquifers and aquifer systems with moderate to high corrosivity potential. The purpose of this study is to provide State and local governments, public water system managers, and the nearly 500,000 private well owners in Louisiana with information needed to manage drinking-water supplies and mitigate potential health risks related to leaching of metals from water pipes and fixtures. The average scores of untreated groundwater samples from approximately 375 wells by index are as follows: LSI, −1.28; RSI, 9.78; PSI, 9.34; and CI, 4.14. The PPGC does not produce a numerical score, but the total percentage of class counts can be used to assign a classification; overall, samples in Louisiana were classified as significant concern. The percentages of groundwater samples from wells classified as potentially corrosive, by index, are as follows: LSI, 53 percent; RSI, 94 percent; PSI, 81 percent; PPGC, 98 percent; and CI, 81 percent. The percentages of samples classified as indeterminate, by index, are as follows: LSI, 46 percent; RSI, 5 percent; PSI, 12 percent; PPGC, 0 percent; and CI, 18 percent.

Louisiana↗

USGS2024

Satellite interferometry landslide detection and preliminary tsunamigenic plausibility assessment in Prince William Sound, southcentral Alaska

Regional mapping of actively deforming landslides, including measurements of landslide velocity, is integral for hazard assessments in paraglacial environments. These inventories are also critical for describing the potential impacts that the warming effects of climate change have on slope instability in mountainous and cryospheric terrain. The objective of this study is to identify slow-moving landslides in the Prince William Sound region, southcentral Alaska, United States, which has had rapid deglaciation since the mid-1800s, and assess their tsunamigenic plausibility. We use an automated time series persistent scatterer interferometric synthetic aperture radar processing method with 7 years of Sentinel-1 data (2016–22) to identify 43 slow-moving slopes with average velocities ranging from approximately 0.2 to 21 millimeters per year. Landslide presence is confirmed using aerial imagery and previous landslide inventory records. We assess the tsunamigenic plausibility of the landslides using empirically derived estimates of landslide mobility based on modeled landslide volumes. Of the identified landslides, our preliminary analysis suggests that 11 have tsunamigenic potential if they were to fail rapidly and catastrophically. Although our estimate of tsunamigenic plausibility is preliminary and can be refined with additional observations and analyses, it can be used to prioritize ongoing and future hazard assessment, surveillance, and research efforts.

Alaska↗

USGS1983

Base flow and ground water in upper Sweetwater Valley, Tennessee

The upper Sweetwater Valley area has a flow system with complex interaction between surface and ground water. A water budget study indicated that during dry years approximately three-fourths of the annual flow to Sweetwater Creek may be derived from ground-water sources. Hydrograph analysis showed seasonal variation of recharge to the ground-water flow system. Streamflow records were analyzed to estimate the frequency of low flow of Sweetwater Creek at river mile 16.7, and indicated the lowest average flow for 1 day in 20 years to be about 5.1 cubic feet per second. Two periods of base-flow measurements of Sweetwater Creek identified channel reaches with significant gains and losses of streamflow. Base flow measurements also showed interbasin transfer of water among sub-basins of the valley. Major flow surpluses were associated with areas in which the majority of flow originated at a spring. Topographically low areas adjacent to the main stem of Sweetwater Creek generally have surplus flow. Topographically higher areas generally have deficient surface outflow unless significant spring flow occurs in the basin. Ground-water recharge occurs by water draining into sinkholes, faults, and fractures. Ground-water flow is regionally diffused across formation strikes from the topographically low areas unless the water is exposed to highly permeable formations or impervious formations. Ground water infiltrates the highly permeable formations and flows along strike. Ground water encountering impervious formations may discharge at small springs at the contact, or may reroute along the contact if the upgradient rock is sufficiently permeable or has well developed secondary porosity. Groundwater discharges to streams at innumerable seeps and at a few large springs. Areas of ground-water flow up-gradient of large springs are hypothesized as likely areas of significant ground-water reservoirs.

Tennessee↗

USGS2026

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting. The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore. We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors. The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2024

Spatial differences in soil nutrients along a hydrographic gradient on floodplains in Dongting Lake

The spatial heterogeneity of soil nutrients is crucial for the water bird and whole floodplain wetland ecosystem in large lakes, and it is influenced by the dramatic water level changes and sedimentation progress in West Dongting Lake (WDL). Soil samples were collected at various soil depths along the Yuan River and Li River that feed into WDL. The concentrations of soil total organic carbon (TOC), total nitrogen (TN), total phosphorus (TP), and soil grain size were tested. The stoichiometric ratios of C, N, P, and the mean value of soil grain size ( Mz ) were calculated. The differences of soil TOC, TN, TP and the stoichiometric ratio at different sites and soil depths were compared. Linear regression was used to explore the relationships of Mz and nutrient concentrations, and relationships between TOC, TN, and TP. Redundancy analysis was used to explore the relationship between soil nutrients, heavy metal concentrations, and plant community diversity. The results showed that the distributions of soil TOC, TN, and TP concentrations differed across regions in west Dongting Lake along the Yuan and Li Rivers. Total organic carbon concentration differed at different sedimentation depths. Soil grain size showed negative effect with soil TOC, TN, and TP concentrations in this region. Plant community diversity correlated positively with soil TOC and negatively with Hg. West Dongting Lake was N limited despite the high wet deposition of N. It could potentially be attributed to the insufficient presence of aerobic environments for microbes during intermittent flooding of the floodplain, coupled with feeble mineralization. This study can provide valuable insights for the conservation of water bird habitats and wetland ecosystems.

West Dongting Lake↗

USGS2026

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗