Low‐elevation coastal areas are increasingly vulnerable to seawater flooding as sea levels rise and the frequency and intensity of large storms increase with climate change. Seawater flooding can lead to the salinization of fresh coastal aquifers by vertical saltwater intrusion (SWI). Vertical SWI is often overlooked in coastal zone threat assessments despite the risk it poses to critical freshwater resources and salt‐intolerant ecosystems that sustain coastal populations. This review synthesizes field and modeling approaches for investigating vertical SWI and the practical and theoretical understanding of salinization and flushing processes obtained from prior studies. The synthesis explores complex vertical SWI dynamics that are influenced by density‐dependent flow and oceanic, hydrologic, geologic, climatic, and anthropogenic forcings acting on coastal aquifers across spatial and temporal scales. Key knowledge gaps, management challenges, and research opportunities are identified to help advance our understanding of the vulnerability of fresh coastal groundwater. Past modeling studies often focus on idealized aquifer systems, and thus future work could consider more diverse geologic, climatic, and topographic environments. Concurrent field and modeling programs should be sustained over time to capture interactions between physical processes, repeated salinization and flushing events, and delayed aquifer responses. Finally, this review highlights the need for improved coordination and knowledge translation across disciplines (e.g., coastal engineering, hydrogeology, oceanography, social science) to gain a more holistic understanding of vertical SWI. There also needs to be more education of communities, policy makers, and managers to motivate societal action to address coastal groundwater vulnerability in a changing climate.
Future increases in the frequency of tidal flooding due to sea level rise (SLR) are likely to affect pore water salinities in coastal aquifers. In this study, we investigate the impact of increased tidal flooding frequency on salinity and flow dynamics in coastal aquifers using numerical variable‐density variably‐saturated groundwater flow and salt transport models. Short (sub‐daily) and long (decadal) period tides are combined with SLR projections to drive continuous 80‐year models of flow and salt transport. Results show that encroaching intertidal zones lead to both periodic and long‐term vertical salinization of the upper aquifer. Salinization of the upper aquifer due to tidal flooding forces the lower interface seaward, even under SLR. System dynamics are controlled by the interplay between SLR and long period tidal forcing associated with perigean spring tides and the 18.6‐year lunar nodal cycle. Periodic tidal flooding substantially enhances intertidal saltwater‐freshwater mixing, resulting in a 6‐ to 10‐fold expansion of the intertidal saltwater‐freshwater mixing area across SLR scenarios. The onset of the expansion coincides with extreme high water levels resulting from lunar nodal cycling of tidal constituent amplitudes. The findings are the first to demonstrate the combined effects of gradual SLR and short and long period tides on aquifer salinity distributions, and reveal competing influences of SLR on saltwater intrusion. The results are likely to have important implications for coastal ocean chemical fluxes and groundwater resources as tidal flooding intensifies worldwide.
more » « less- Award ID(s):
- 1933058
- NSF-PAR ID:
- 10443976
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 12
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Knowledge of coastal groundwater flow is critical for managing coastal groundwater resources and quantifying submarine groundwater discharge (SGD), but this flow occurs over multiple scales that can be difficult to study in an integrated way. We designed a field and modeling study to investigate groundwater flow and the distribution of salinity during sea level rise in a domain that included beaches, salt marshes and the first major confined aquifer, which reached 10–15 km offshore. Numerical models were based on the flat‐lying, passive margin coastline of North Inlet, SC, and were constrained by field studies including subsurface resistivity surveys and hydraulic head observations. Simulations that included tidal fluctuations showed that the salt marsh generated more than three times as much SGD as the beach and inner shelf, per unit length of coastline. Groundwater exchange between scales was small, suggesting that physical fluxes of groundwater can be considered independently at different scales. However, salinization of the first major confined aquifer occurred by downward transport from overlying aquifers rather than intrusion from the seaward end, suggesting that studies of aquifer salinization should consider multiscale flow. During simulated sea level rise, fresh‐to‐brackish groundwater persisted in the first confined aquifer as far as the seaward end of the overlying confining unit, 10–20 km offshore. Total fluxes of SGD decreased significantly with future sea level rise, dominated by declining SGD in the salt marsh, and portending a marked decline in the flux of nutrients and carbon to estuaries and the coastal ocean.
-
more » « less
Abstract Intertidal aquifers are hotspots of biogeochemical cycling where nutrients and contaminants are processed prior to discharge to the ocean. The nature of the dynamic subsurface mixing zone is a critical control on mitigating reactions. Simulation of density‐dependent, variably saturated flow and salt transport incorporating realistic representations of aquifer heterogeneity was conducted within a Monte Carlo framework to investigate influence of nonuniform permeability on intertidal groundwater flow and salt transport dynamics. Results show that heterogeneity coupled with tides creates transient preferential flow paths within the intertidal zone, evolving multiple circulation cells and fingering‐type salinity distributions. Due to heterogeneity, strain‐dominated (intense mixing) and vorticity‐dominated (low mixing) flow regions coexist at small spatial scales, and their spatial extent reaches peaks at high tide and low tide. Such topological characteristics reveal complex tempo‐spatial mixing patterns for intertidal flow with localized areas of high and low mixing intensities, which have implications for intertidal biogeochemical processing.
-
more » « less
Abstract Preferential flow can result in rapid contamination of groundwater resources. This is particularly true in aquifers with connected, high permeability geologic structures and in coastal systems where the oceanic source of contamination is ubiquitous. We consider saltwater intrusion due to pumping in volcanic aquifers with lava tubes represented as connected high‐K structures and compare salinization responses to those of heterogeneous aquifers with different structure and equivalent homogeneous systems. Three‐dimensional simulations of variable‐density groundwater flow and salt transport show that conduits formed by lava flows create preferential groundwater flow in volcanic aquifers. These conduits allow fresh groundwater to extend further offshore than in other systems. However, onshore pumping causes saltwater to migrate landward quickly through the conduits relative to the other models, resulting in more severe saltwater intrusion, particularly at shallow depths. The geometry of geologic heterogeneity in volcanic aquifers leads to increased risk of salinization of fresh groundwater as well as substantial uncertainty due to significant spatial variation in saltwater intrusion. The findings illustrate the importance of considering geologic heterogeneity in assessing the vulnerability of coastal freshwater resources in volcanic and other aquifers with connected high‐permeability geologic structures.
-
more » « less
Abstract Coastal aquifers supply freshwater to nearly half the global population, yet they are threatened by salinization. Salinities are typically estimated assuming steady‐state, neglecting the effect of cyclical forcings on average salinity distributions. Here, numerical modeling is used to test this assumption. Multi‐scale fluctuations in sea level (SL) are simulated, from tides to glacial cycles. Results show that high‐frequency fluctuations alter average salinities compared with the steady‐state distribution produced by average SL. Low‐frequency forcing generates discrepancies between present‐day salinities estimated with and without considering the cyclical forcing due to overshoot effects. This implies that salinities in coastal aquifers may be erroneously estimated when assuming steady‐state conditions, since present distributions are likely part of a dynamic steady state that includes forcing on multiple timescales. Further, typically neglected aquifer storage characteristics can strongly control average salinity distributions. This has important implications for managing vulnerable coastal groundwater resources and for calibration of hydrogeological models.
