Surface Water Groundwater Research Articles

Storm surge, seawater flooding, and sea-level rise paradoxically drive fresh surface water expansion

Abstract Coastal storms and sea-level rise are expected to increase seawater flooding in low-elevation coastal zones. High sea levels and seawater flooding can drive groundwater table rise via ocean-aquifer connections. These dynamics are often overlooked but can cause groundwater flooding and salinization hazards, increasing freshwater security challenges for coastal communities and driving ecosystem transgressions. Field data and numerical modeling were used to evaluate how heavy rainfall, storm surge, and seawater flooding and infiltration during Hurricane Fiona (September 2022) and projected sea-level rise impact groundwater levels, inland surface waters, and saltwater intrusion on Sable Island National Park Reserve, Canada. During the passage of Hurricane Fiona, precipitation increased groundwater and pond levels before seawater flooded the beach. Seawater flooding and infiltration caused a sharp rise in beach groundwater levels, which in turn caused inland pond levels to rise without coincident direct inputs from precipitation or seawater. Model simulations reveal that seawater infiltration on beaches flooded the subsurface and drove the observed inland groundwater rise and freshwater pond expansion. Simulations of projected sea-level rise show that seawater flooding will only inundate a small area of land along the coast; however, inland groundwater rise and flooding, which is less well-studied, may inundate up to 30 times more land area. Further, groundwater flooding driven by rising sea levels decreases hydraulic gradients and increases saltwater intrusion via freshwater lens contraction. Findings demonstrate that seawater flooding from coastal storms and sea-level rise paradoxically cause concurrent fresh surface water expansion but freshwater lens contraction. This study provides new insights into the spatiotemporal dynamics of island freshwater resources and highlights that unseen and often overlooked groundwater-surface water exchanges are critical to consider when evaluating coastal flooding and groundwater salinization hazards and management strategies for low-elevation coastlines.

Groundwater-surface water exchanges in an alluvial plain in southern France subjected to pumping: A coupled multitracer and modeling approach

Study regionThe study was conducted in an alluvial plain between the Rhône and the Ouvèze Rivers (in the southeast of France) extensively exploited for drinking water. The research area is characterized by significant groundwater-surface interactions influenced by groundwater pumping activities. Study focusThe aim of this study is to enhance the understanding of interactions between rivers and alluvial aquifers by combined multi-tracer and numerical modeling approaches. Over an 18-month period, groundwater temperature, piezometric levels, and river surface water levels were continuously monitored. Field campaigns focused on conductivity, stable isotopes of water, and radon-222 activity concentration in both groundwater and surface water. Radon-222 was used to quantify water exchanges between the river and the aquifer. A MODFLOW model, calibrated using piezometric data and PEST, was employed to simulate groundwater flow and reactive transport of radon-222 using MT3DMS. New hydrological insights for the regionThe study reveals that river water recharges the aquifer, with radon-222 data delineating this recharge zone. The methodology extended the interpretation of periodic groundwater temperature signals to isotopic signals, allowing the identification of dispersivity and Darcy's velocity. The Ouvèze River was found to contribute approximately 55 % of the pumping water supply, alongside the Rhône. These findings provide valuable insights for sustainable water resource management, demonstrating the relevance of using natural tracers in scenarios where artificial tracers are impractical.

Methods for Quantifying Interactions Between Groundwater and Surface Water

Driven by the need for integrated management of groundwater (GW) and surface water (SW), quantification of GW–SW interactions and associated contaminant transport has become increasingly important. This is due to their substantial impact on water quantity and quality. In this review, we provide an overview of the methods developed over the past several decades to investigate GW–SW interactions. These methods include geophysical, hydrometric, and tracer techniques, as well as various modeling approaches. Different methods reveal valuable information on GW–SW interactions at different scales with their respective advantages and limitations. Interpreting data from these techniques can be challenging due to factors like scale effects, heterogeneous hydrogeological conditions, sediment variability, and complex spatiotemporal connections between GW and SW. To facilitate the selection of appropriate methods for specific sites, we discuss the strengths, weaknesses, and challenges of each technique, and we offer perspectives on knowledge gaps in the current science.

Spatio-temporal and depth-oriented evaluation of hyporheic zone processes in headwater catchments

The hyporheic zone (HZ) is vital for nutrient cycling and overall stream ecosystem health across groundwater-surface water interfaces, with biogeochemical and physical variations influencing exchange processes. To gain a depth-oriented insight into the various hyporheic functioning, isotopic (18O and 2H) and chemical analyses, along with physical parameters, were performed in two streams named Ahna and Losse in North Hesse, Germany. Multi-level interstitial probes were installed to extract subsurface pore water samples up to 0.45 m depth. We identified three distinct HZ in Ahna upstream, downstream and Losse. The Ahna upstream HZ was a less permeable barrier with low water fluxes but reactive, while the downstream HZ was less reactive with higher permeability and vertical fluxes. The Losse HZ, remained non-reactive but well-mixed with good hydraulic circulation due to fast downward flows. Isotopic composition revealed distinctive features, indicating increased evaporation and intricate mixing with other water sources and intensified subsurface flow in Ahna downstream compared to upstream. Hydro-chemical dynamics in Ahna exposed increased ion concentrations downstream, driven by anthropogenic factors along the course. Upstream denitrification was dominant under low oxygen conditions due to the slow infiltration and low mixing rates. The primary water source for both streams and HZs was groundwater, with the potential for rainwater to serve as a secondary source. This study shows the spatial heterogeneity of HZ processes and highlights that the specific HZ conditions need to be examined locally.

A multi-objective optimization and evaluation framework for LID facilities considering urban surface runoff and shallow groundwater regulation

Currently, the goals of sponge city construction mostly focus on urban surface water, with few considerations of groundwater indicators, and there is a lack of simulation methods that consider the interaction between surface water and groundwater under low impact development (LID) scenarios. Taking a sponge campus in Xi'an as the research object, this paper proposes a method to realize the interaction between surface water and groundwater at different spatial and temporal scales, so as to construct a SWMM-MODFLOW coupled numerical model to carry out the analysis of hydrological effect of sponge city comprehensively taking into account both surface water and groundwater. Simultaneously, the SWMM is coupled with the NSGA-III algorithm for multi-objective optimization of LID facility configurations, aiming to identify a globally optimal solution set for LID facility configurations under both long and short-duration rainfall scenarios. A novel framework for a comprehensive evaluation system of surface water-groundwater benefits is also developed to comprehensively assess LID facility optimization schemes that meet the goals of stormwater runoff infiltration and reduction, pollution interception and purification, and groundwater recharge conservation. The simulation results show that after optimization, the LID scheme has good surface runoff regulation effects and pollution reduction capabilities, reducing the runoff peak from 2.297 m3/s to 2.128 m3/s, significantly lowering the risk of inundation and effectively reducing the total suspended solids (SS) load by 42.38 kg at the outfall. Meanwhile, the groundwater recharge effect is particularly significant in September. Compared to the baseline scenario, the average groundwater levels from June to September were raised by 3.966 cm, 5.120 cm, 6.743 cm, and 7.639 cm, respectively. The maximum daily groundwater level rise reached 7.82 cm, while the average groundwater level for the whole year of 2023 increased by 2.61 cm. The maximum transport area of pollutants decreased from 2359.44 m2 to 1796.54 m2, and the maximum transport distance reduced from 79.985 m to 69.977 m. Additionally, the central concentration of pollutants also decreased. This optimization evaluation framework will aid decision-makers in selecting LID facility management strategies that balance the dual benefits of surface water runoff and groundwater regulation.

Historical patterns of well drilling and groundwater depth in Arizona considering groundwater regulation and surface water access

Historical patterns of well drilling and groundwater depth in Arizona considering groundwater regulation and surface water access

Comprehensive, Multi‐Scale Evaluation of Field Methods for Assessing Stream–Aquifer Interactions Along Channelised Lowland Streams

ABSTRACTStream–aquifer interactions (SAIs) play a critical role in effective groundwater management, yet their complex dynamics remain poorly understood in channelized lowland perennial streams. This study presents a multi‐scale, multi‐technique investigation of SAIs along two long‐stream reaches in Tennessee, United States. The goal is to define a suitable methodology for characterising SAIs in this specific hydrological setting, serving as a starting point for developing a more standardised and replicable approach. The methodology includes an initial evaluation of various field techniques, followed by extensive surveys using potentiomanometers, electromagnetic induction (EMI), vertical temperature profilers (VTPs) and complementary methods such as seepage meters, bank tests and well‐data analyses. Results reveal distinct hydrogeomorphic behaviours across and along the streams, challenging the SAI‐homogeneity notions typically assumed in groundwater models. Nonconnah Creek exhibited streambed colmation and negligible hydraulic gradients, resulting in disconnection from the aquifer during low flows, except for a 300‐m losing reach with high downward gradients. In contrast, the Loosahatchie River displayed relatively homogeneous streambed properties and small, upward hydraulic gradients, suggesting uniform SAIs along the surveyed reaches. EMI proved highly effective for mapping streambed sediments quickly, while potentiomanometers accurately measured small head differences critical for understanding SAI dynamics. VTPs were less practical due to the extended data‐collection times required and their vulnerability to flooding. This study emphasises the importance of multi‐scale investigations using diverse techniques to accurately characterise SAIs in lowland streams, highlighting the confounding influences of geological formations, anthropogenic alterations and depositional processes on groundwater–surface water interactions. The findings contribute to refining local water balances, informing groundwater management strategies and underscoring the need for incorporating local‐scale field data into regional groundwater models. The proposed methodology serves as a foundation for developing a standardised approach for characterising SAIs in lowland channelized perennial streams, adaptable for similar stream systems worldwide.

Tracing Groundwater-Surface Water Sources and Transformation Processes in the Ba River Basin through Dual Isotopes and Water Chemistry

Tracing Groundwater-Surface Water Sources and Transformation Processes in the Ba River Basin through Dual Isotopes and Water Chemistry

Interaction of pans and groundwater: A case study in the Khakhea- Bray Transboundary aquifer portion in South Africa

Understanding the groundwater-surface water interaction is important for the sustainable management of water resources and ecology. Most of the studies on groundwater-surface interactions are confined to aquifer-river and aquifer-lake interactions. The interactions of water pans and groundwater have not received much attention. The study aims to provide insights into the interaction of the pans and groundwater at a case study in the Khakhea-Bray dolomite Transboundary Aquifer (TBA) shared by Botswana and South Africa. The study uses field observations, geology, hydrogeophysics, and stable isotopes complimentary tools. Complimentary tools provide several evidence to corroborate the findings. The hydrogeophysics taps into the recently developed Telluric Electric Frequency Selection Method (TEFSM). Groundwater discharge zones/springs were delineated at the geological contact of the boundary of the Khakhea-Bray dolomite transboundary aquifer and the granite of the Kalahari group. 4 of the 33 research pans are in the groundwater discharge zone and are classified as groundwater dependent. The findings highlight the benefits of complimentary approaches to corroborate the findings on groundwater-surface water interactions. The TEFSM hydrogeophysics can be able to delineate groundwater discharge zones, aquifer boundaries, and the mapping of TBAs. More studies are recommended to test the TEFSM hydrogeophysical approach to study groundwater-surface water in different hydrological environments.

Constraining seasonal and spatial ambient and urban groundwater contributions to streamflow across a mixed land-use regional-scale Precambrian Shield watershed

Groundwater-surface water (GW-SW) interactions are complex phenomena that vary naturally over space and time and are being influenced by environmental change. The Whitson River watershed is a mixed land-use Precambrian shield watershed located in Northeastern Ontario, Canada in which groundwater is a source of municipal water supply and in which groundwater and surface waters are at risk of potential impacts from urban runoff, ongoing municipal drainage projects, periodic flooding, and climate change. This study used watershed-scale synoptic surveys of stable isotopes (δ18O and δ2H) and geochemistry, and corresponding mixing-model approaches (two-component and three-component) to quantify the seasonal and spatial variation in surface water, ambient, and urban groundwater contributions to streamflow. Multiple mixing models were generated based on isotope and geochemical source water identification, and the performance of these models was compared. Mixing-model analyses showed that groundwater contributes a critical, sustaining source to streamflow (> 50 % or more) during summer baseflow conditions, dropping to ∼ 15 % in fall when connection to shield lakes and wetlands dominate. Three-component mixing models showed similar results to the two-component models, refining groundwater contributions into ambient and urban groundwater sources. Estimated urban groundwater contributed 30 % or more to summer baseflow, dropping to ∼ 4 % in the fall. Baseflow estimates derived from graphical hydrograph separation were similar to mixing model groundwater estimates in the summer but suggest this approach overestimates groundwater contributions during the fall. These results demonstrate the complexity of source water contributions to streamflow and the challenges in their determination in mixed-land-use Precambrian Shield landscapes where communities are largely located. Assessment of source water contributions to streamflow across the Whitson River sub-watershed has generated a region-specific conceptualization (e.g. urban versus ambient groundwater) of a basin that is experiencing increased urban development, providing information that will be of value for long term management.

Delineating groundwater-surface water interaction in a lake watershed (Sapanca Lake, NW Turkey)

Delineating groundwater-surface water interaction in a lake watershed (Sapanca Lake, NW Turkey)

Groundwater-Surface water interactions research: Past trends and future directions

Interactions between groundwater and surface water sustain groundwater-dependent ecosystems and regulate river temperature and biogeochemical cycles, amongst many other processes. These interactions occur in freshwater environments including rivers, springs, lakes, and wetlands, and in coastal environments via tidal pumping, submarine groundwater discharge, and seawater intrusion. Here, we explore groundwater-surface water interactions research using bibliometric analyses of titles, abstracts, and keywords from 20,275 journal papers published between 1970 and 2023 extracted from Scopus. Analyses show that research into groundwater-surface water interactions is highly multi-disciplinary, with growing contributions from the social and biological sciences. The number of groundwater-surface water interactions papers is rapidly increasing with over 1200 papers published per year since 2020. Drawing on our data-driven approach and expert knowledge, we synthesise current research trends and identify critical future research directions. Despite the thousands of papers on groundwater-surface water interactions, important processes are still difficult to quantify or predict at meaningful spatial scales to inform water-resources management. We see benefits in future groundwater-surface water interactions research focusing on: (1) using new technologies including internet-of-things-based sensors, uncrewed vehicles, and remote-sensing approaches for data collection to inform groundwater-surface water interactions at large scales, (2) seeking approaches to upscale site-specific findings to better inform management, and (3) continuing the movement towards multi-disciplinary investigations to better inform the understanding of groundwater-surface water interactions and processes that will enable better management outcomes.

Model simplification to simulate groundwater recharge from a perched gravel-bed river

Gravel-bed rivers are an important source of groundwater recharge in some regions of the world. Their interactions with groundwater are complex and highly variable in space and time, with considerable water storage in the riverbed sediments. In losing river sections, where most of the groundwater recharge occurs, the river can be separated from the regional groundwater system by an unsaturated zone (i.e., perched). The complexity of groundwater–surface water interactions in these environments calls for the use of 3D fully integrated hydrological models to represent them, but their computational intensity limits their practicality for parameter inference, uncertainty quantification and regional scale problems. On the other hand, the simple groundwater–surface water exchange functions currently implemented in regional scale groundwater models are not suited to represent complex gravel-bed river systems such as braided rivers. There is therefore a need for developing groundwater–surface water exchange functions tailored to gravel-bed rivers that can be used in regional scale models.To address this issue, we developed a model simplification framework that combines a 3D integrated surface and subsurface hydrological model, a 2D cross-sectional river-aquifer model and a 1D conductance-based analytical model. We aim at broadly simplifying the 3D model while ensuring the appropriate simulation of groundwater recharge. We demonstrate our modelling approach on the Selwyn River (New Zealand) using piezometric data and groundwater recharge estimates derived from field observations and satellite imagery.Our results indicates that groundwater recharge from this river can be simulated using a simple 1D analytical model, which can easily be implemented in regional groundwater models (e.g., MODFLOW models). However, to represent properly the time variability of groundwater recharge, it is essential to use the groundwater level in the shallow aquifer associated with the river as input to the regional groundwater model. Our approach is generally transferable to other gravel-bed rivers but requires some observations of river losses for proper calibration.

Influence of Dissolved Oxygen and Temperature on Nitrogen Transport and Reaction in Point Bars of River

Point bars are crucial elements of river systems, significantly enhancing the nitrogen cycle in riparian zones by facilitating hyporheic exchange between surface water and riparian zones. This study investigated the impact of dissolved oxygen (DO) concentration and temperature on nitrogen transport and reactions in river point bars. A two-dimensional coupled surface water–groundwater model was developed to analyze nitrogen distribution, variations, and reaction rates in rivers with point bars. The model considered three chemical reactions controlling nitrogen transformation: aerobic respiration, nitrification, and denitrification, with DO and temperature as independent variables. The results indicated that DO variations have a limited effect on solute migration depth, whereas increased temperature reduces solute migration depth. At surface water DO concentrations of 0.1, 0.2, and 0.4 mol/m3, nitrate removal in the riparian zone was 0.022, 0.0064, and 0.0019 mol/m, respectively. At riparian temperatures of 5 °C, 15 °C, and 25 °C, nitrate removal was 0.012, 0.041, and 0.16 mol/m, respectively. Nitrogen removal is more sensitive to temperature variations than to changes in DO concentration. In this research, the decrease in DO concentrations and the temperature increase greatly enhanced the riparian zone’s denitrification effect. This study improves our understanding of how riparian zones impact nitrogen cycling under various environmental conditions.

A coupled regional-scale numerical model for hydrological processes and interactions between groundwater and surface water in a controlled drainage district

Controlled drainage (CD) has emerged as an effective strategy to prevent field waterlogging and mitigate drought during crop growth by altering the hydrological process, while few studies have quantified its effect on groundwater-surface water interactions and groundwater recharge at a regional scale. In this study, a new model is developed for the coupled numerical simulation of water flow in ditches, soil, and underground aquifers under CD conditions. The domain is partitioned into horizontal sub-areas and two vertical layers, with ditches discretized into segments based on groundwater flow grids. The Saint-Venant equations, Richards equation, and MODFLOW are applied to describe water movement processes in ditches, soil, and underground aquifers, respectively. The proposed model is calibrated and validated using five years of observations of ditch water levels (DWL) and groundwater levels, demonstrating excellent agreement with observed data. Subsequently, water balance components of the shallow aquifer below 1 m depth (GW), 1 m of root layer soil profile (SW), and ditch water (DW) under seven scenarios with different control schemes during flood and dry seasons are analyzed to quantify variations in hydrological processes caused by CD. The results show that CD significantly impacts water within SW, GW, and DW, soil evaporation, infiltration, and net recharge from GW to SW (GtoS) by altering regional groundwater table depth (GWT) through the exchange between DW and GW (DtoG and GtoD). GWT and DW show moderate correlation with GtoD, while their correlation with DtoG varies significantly under different rainfall conditions. Smaller precipitation results that precipitation and GWT have a higher linear correlation with water balance components during the dry season compared to the flood season. The correlation values (r) between transpiration (T) and GWT range from −0.99 to 0.96 and −1 to −0.85 during the flood and dry seasons, respectively, indicating a pronounced influence of CD on crop growth. On average, a 1 m increase in DWL control scheme leads to a 0.51 m increase in DWL, and regional GWT decreases by 0.35 m and 0.24 m during flood and dry seasons, respectively. Furthermore, for every 1 m decrease in GWT, net GtoS increases by 36.1 mm and 28.1 mm, respectively, resulting in an increase in SW by 13.9 mm and 11.8 mm, respectively. Considering the varying main crop water stress and the impact of CD under different rainfall conditions, an appropriate CD scheme needs to be adjusted accordingly. These findings provide a quantitative reference for the effect of CD on regional hydrological process and serve as a typical case for similar areas aiming to implement CD strategies for waterlogging and drought control.

Nitrogen fate in riparian zones: Insights from experiments and analysis of sediment porosity and surface water-groundwater exchange

Riparian zones play a vital role in the river ecosystem. Solutes in vertical riparian zones are transported being by alternating hydraulic gradients between river water and groundwater, due to natural or human activities. This study investigates the impacts of porous sediments and alternating rate of surface water-groundwater on nitrogen removal in the riparian zone through experiments based on the field sampled. The experimental results, combined with dimensionless numbers (Péclet and Damköhler) and Partial Least Squares-Path Modeling, analyze the nitrogen fate responding to hydrodynamics changes. The results show that increased sediment porosity contributes to the ammonium removal, particularly when the oxygen content of river water is low, with the removal rate up to 72.57%. High ammonium content and dissolved organic carbon (DOC) in rural rivers lead to a constant low-oxygen condition (4 mg/L) during surface water-groundwater alternation, and promote denitrification. This threatens groundwater with ammonium pollution and causes accumulation at the top of vertical riparian zones during upwelling, potentially causing secondary river pollution. However, increasing the alternating rate hinders the nitrate denitrification and drastically changes in the redox environment of the riparian zone, despite contributing to ammonium removal. Rapid oxygen consumption during aerobic metabolism and nitrification in groundwater-surface water exchange created favorable conditions for denitrification. Floodplains sediment porosity is unfavorable for nitrification. This study improves understanding of coupled hydrologic and solute processes in vertical riparian zones, informing strategies for optimizing nitrogen attenuation and riparian zone construction.

Uncovering the hidden world of riverbed sediments: The role of sediment heterogeneity and cross-bar channel fills in the hydrogeochemical dynamics of the hyporheic zone

Groundwater-surface water interaction (hyporheic exchange) is critical in numerous hydrogeochemical processes; however, hyporheic exchange is difficult to characterize due to the various spatial (e.g., sedimentary architecture) and temporal (e.g., stage fluctuations) variables that influence it. This interdisciplinary study brings forth novel insights by integrating various methodologies including geophysical surveys, physical and chemical sediment characterization, and water chemistry analysis to explore the interplay of the numerous facets governing hyporheic zone processes within a compound bar deposit. The findings reveal distinct sedimentary facies and geochemical zones within the compound bar, driven by the sedimentary architecture. Cross-bar channel fills are identified as critical structures influencing hydrogeochemical dynamics, acting as baffles to groundwater flow and modulating nutrient transformations. Geophysical imaging and hydrogeochemical analyses highlight the complex interplay between sediment characteristics and subsurface hydraulic connectivity, emphasizing the role of sediment heterogeneity in controlling hyporheic exchange and solute mixing. The study concludes that sediment heterogeneity, particularly the presence of cross-bar channel fills, plays a pivotal role in the hydrogeochemical dynamics of the hyporheic zone. These structures significantly influence hyporheic flow paths, solute residence times, and nutrient cycling, underscoring the necessity to consider the fine-scale sedimentary architecture in models of hyporheic exchange. The findings contribute to a deeper understanding of riverine ecosystem processes, offering insights that can inform management strategies for water quality and ecological integrity.

Iraqi Meteoric Water Line and its Application in Groundwater-Surface Water Connection

The LMWL was re-plotted from the isotopic analysis data (δ2H and δ18O) for rainwater samples collected from three stations in Baghdad city in the current study and data from previous studies. Four samples of shallow groundwater and five samples from the Tigris River were isotopically analyzed for 2H and 18O using a Liquid-Water Stable Isotope Analyzer (LWSIA). The samples were collected from the Al-Jaderia area within the University of Baghdad. The analysis was done against Vienna Standard Main Oceanic Water (VSMOW). The isotopic composition of groundwater ranges between -43.91‰ and -37.57 ‰ for 2H with an average of -40.06±2.7 ‰, while it ranges between -6.28‰ and -6.01‰ for 18O with an average of -6.14±0.1‰. These values are slightly more depleted compared to the isotopic composition of the Tigris River water (-38.79 ‰ to -37.87 ‰) for 2H with an average of -38.05±0.9‰ and between -6.23‰ to -5.11 ‰ for 18O with an average of -5.74±0.4 ‰. This is a logical result since the river is usually fed with fresh new rainwater, which isotopically is relatively more enriched. The current study confirms that the Local Meteoric Water Line (LMWL) has a formula of δ 2H‰ = 7.5 δ 18O‰ + 13. The LMWL is located above the Global MWL by +13.9‰ toward the Eastern Mediterranean (EMMWL δ2H‰ = 8 δ 18O‰ + 22), indicating the origin rainwater from the Mediterranean Sea and the Atlantic Ocean. The groundwater and surface water isotopic compositions were below the LMWL due to the evaporation of shallow groundwater. There is interconnection between river water and groundwater, directly or by recharge from irrigation water. Iraq may have two different MWLs, one for the northeastern region and the other for the middle, western and southern Iraq. The MWL for the northern eastern region of Iraq is close to that of the EMMWL, indicating Mediterranean climate and vapor source, while the rest of the areas of Iraq have MWL with various isotopic compositions affected by more hot climate in these regions and by the continental tropical air masses.

Per- and polyfluoroalkyl substances (PFAS) fate and transport across a groundwater-surface water interface

Per- and polyfluoroalkyl substances (PFAS) are emerging contaminants of concern whose fate and transport in environmental media are incompletely understood. In the 1960s, PFAS were dumped in the House Street Disposal Site, an unlined landfill on the crest of a glacial end moraine near Rockford, Michigan, USA. In 2017, PFAS were discovered in groundwater and subsequently, a network of monitoring wells delineated a 2 mi (3 km) PFAS plume migrating downgradient toward the Rogue River. Today, the Michigan Department of Natural Resources (MDNR) operates fish-rearing ponds in the area where the plume intersects the groundwater-surface water interface (GSI). Each year, the MDNR fills these man-made ponds using water from a nearby creek. Springs in the ponds prevent them from draining completely at the end of fish-rearing each fall.We sampled surface water and modeled groundwater flow to investigate PFAS transport across the GSI. Numerical models constructed with and without the fishponds did not substantially change MODFLOW model calibration curves or predicted MODPATH flow lines, indicating that PFAS transport is dominated by the regional flow system with limited influence from semiannual changes to boundary conditions at the GSI. Surface water samples collected from five locations within and adjacent to the fishponds were analyzed using EPA Draft Method 1633. PFAS were detected at all locations with the highest total PFAS >60 ng/L in the fishponds. Mixing models based on total PFAS indicate that approximately 10 % of the fishpond water is sourced by groundwater. However, similar analyses with perfluoroalkyl carboxylic acids (PFCA) and perfluoroalkyl sulfonic acids (PFSA) imply that groundwater comprises as much as 30 % of water in the ponds, suggesting differential movement of individual PFAS across the groundwater-surface water interface. Additional investigation of PFAS within the pond sediments is needed to better understand partitioning and differential transport behavior across the GSI.

Strontium isotopes in the atmosphere, geosphere and hydrosphere: Developing a systematic “fingerprinting” framework of rocks and water in sedimentary basins in eastern Australia

Understanding the connection between aquifers, aquitards, and groundwater-dependant ecosystems remains a key challenge when developing a conceptual hydrogeological model. The aim of this study was to develop a systematic strontium isotope (87Sr/86Sr) fingerprinting framework of rocks and water within the sedimentary Surat and Clarence-Moreton basins (SCM basins) in eastern Australia – an area of extensive coal seam gas development and high potential for aquifer and groundwater-surface water connectivity. To do this, new groundwater samples (n = 298) were collected, analyzed and integrated with published data (n = 154) from the basins' major sedimentary, volcanic and alluvial aquifers, including the major coal seam gas target, the Walloon Coal Measures. Samples were also analyzed from rainfall (n = 2) and surface water (n = 40). In addition, rock core samples (n = 39) from exploration and stratigraphic wells were analyzed to determine the range of Sr isotope composition from host rocks. The analyses of cores demonstrate a distinct and systematic contrast in 87Sr/86Sr between different hydrogeological units. This confirms that all major hydrogeological units have a narrow range with unique 87Sr/86Sr population characteristics that are useful for guiding conceptual model development. Comparison with selected hydrochemical and groundwater age tracers (14C and 36Cl) suggests only limited changes of 87Sr/86Sr from recharge beds to the deeper parts of the basins or with a decrease in natural 14C and 36Cl tracer content along flow paths. Stream sampling during baseflow conditions confirms that 87Sr/86Sr in surface waters are similar to those of the underlying bedrock formations. We demonstrated that 87Sr/86Sr analyses of rocks and water provide a powerful hydrostratigraphic and chemostratigraphic fingerprinting framework in the SCM basins, enabling reliable assessments of plausible aquifer and groundwater-surface water interconnectivity pathways. Applied in other complex multi-aquifer sedimentary basins in Australia, and globally, a similar approach can help to constrain conceptual hydrogeological models and facilitate improved water resource management.