Impact of Sea-Level Rise on Sea Water Intrusion in Coastal Aquifers

Despite its purported importance, previous studies of the influence of sea-level rise on coastal aquifers have focused on specific sites, and a generalized systematic analysis of the general case of the sea water intrusion response to sea-level rise has not been reported. In this study, a simple conceptual framework is used to provide a first-order assessment of sea water intrusion changes in coastal unconfined aquifers in response to sea-level rise. Two conceptual models are tested: (1) flux-controlled systems, in which ground water discharge to the sea is persistent despite changes in sea level, and (2) head-controlled systems, whereby ground water abstractions or surface features maintain the head condition in the aquifer despite sea-level changes. The conceptualization assumes steady-state conditions, a sharp interface sea water-fresh water transition zone, homogeneous and isotropic aquifer properties, and constant recharge. In the case of constant flux conditions, the upper limit for sea water intrusion due to sea-level rise (up to 1.5 m is tested) is no greater than 50 m for typical values of recharge, hydraulic conductivity, and aquifer depth. This is in striking contrast to the constant head cases, in which the magnitude of salt water toe migration is on the order of hundreds of meters to several kilometers for the same sea-level rise. This study has highlighted the importance of inland boundary conditions on the sea-level rise impact. It identifies combinations of hydrogeologic parameters that control whether large or small salt water toe migration will occur for any given change in a hydrogeologic variable.     

Influence of Seasonal Variations in Sea Level on the Salinity Regime of a Coastal Groundwater–Fed Wetland

Seasonal variations in sea level are often neglected in studies of coastal aquifers; however, they may have important controls on processes such as submarine groundwater discharge, sea water intrusion, and groundwater discharge to coastal springs and wetlands. We investigated seasonal variations in salinity in a groundwater‐fed coastal wetland (the RAMSAR listed Piccaninnie Ponds in South Australia) and found that salinity peaked during winter, coincident with seasonal sea level peaks. Closer examination of salinity variations revealed a relationship between changes in sea level and changes in salinity, indicating that sea level–driven movement of the fresh water‐sea water interface influences the salinity of discharging groundwater in the wetland. Moreover, the seasonal control of sea level on wetland salinity seems to override the influence of seasonal recharge. A two‐dimensional variable density model helped validate this conceptual model of coastal groundwater discharge by showing that fluctuations in groundwater salinity in a coastal aquifer can be driven by a seasonal coastal boundary condition in spite of seasonal recharge/discharge dynamics. Because seasonal variations in sea level and coastal wetlands are ubiquitous throughout the world, these findings have important implications for monitoring and management of coastal groundwater–dependent ecosystems.     

Sea water intrusion by sea-level rise: scenarios for the 21st century

This study presents a method to assess the contributions of 21st-century sea-level rise and groundwater extraction to sea water intrusion in coastal aquifers. Sea water intrusion is represented by the landward advance of the 10,000 mg/L iso-salinity line, a concentration of dissolved salts that renders groundwater unsuitable for human use. A mathematical formulation of the resolution of sea water intrusion among its causes was quantified via numerical simulation under scenarios of change in groundwater extraction and sea-level rise in the 21st century. The developed method is illustrated with simulations of sea water intrusion in the Seaside Area sub-basin near the City of Monterey, California (USA), where predictions of mean sea-level rise through the early 21st century range from 0.10 to 0.90 m due to increasing global mean surface temperature. The modeling simulation was carried out with a state-of-the-art numerical model that accounts for the effects of salinity on groundwater density and can approximate hydrostratigraphic geometry closely. Simulations of sea water intrusion corresponding to various combinations of groundwater extraction and sea-level rise established that groundwater extraction is the predominant driver of sea water intrusion in the study aquifer. The method presented in this work is applicable to coastal aquifers under a variety of other scenarios of change not considered in this work. For example, one could resolve what changes in groundwater extraction and/or sea level would cause specified levels of groundwater salinization at strategic locations and times.     

A Geochemical Approach to Determine Sources and Movement of Saline Groundwater in a Coastal Aquifer

Geochemical evaluation of the sources and movement of saline groundwater in coastal aquifers can aid in the initial mapping of the subsurface when geological information is unavailable. Chloride concentrations of groundwater in a coastal aquifer near San Diego, California, range from about 57 to 39,400 mg/L. On the basis of relative proportions of major‐ions, the chemical composition is classified as Na‐Ca‐Cl‐SO₄, Na‐Cl, or Na‐Ca‐Cl type water. δ²H and δ¹⁸O values range from ?47.7‰ to ?12.8‰ and from ?7.0‰ to ?1.2‰, respectively. The isotopically depleted groundwater occurs in the deeper part of the coastal aquifer, and the isotopically enriched groundwater occurs in zones of sea water intrusion. ⁸⁷Sr/⁸⁶Sr ratios range from about 0.7050 to 0.7090, and differ between shallower and deeper flow paths in the coastal aquifer. ³H and ¹⁴C analyses indicate that most of the groundwater was recharged many thousands of years ago. The analysis of multiple chemical and isotopic tracers indicates that the sources and movement of saline groundwater in the San Diego coastal aquifer are dominated by: (1) recharge of local precipitation in relatively shallow parts of the flow system; (2) regional flow of recharge of higher‐elevation precipitation along deep flow paths that freshen a previously saline aquifer; and (3) intrusion of sea water that entered the aquifer primarily during premodern times. Two northwest‐to‐southeast trending sections show the spatial distribution of the different geochemical groups and suggest the subsurface in the coastal aquifer can be separated into two predominant hydrostratigraphic layers.     

Stack unit mapping of coastal aquifer to predict and control sea water intrusion using remote sensing and a geographical information system

Aquifers are inherently susceptible to contamination and coastal aquifers in specific are highly vulnerable to sea water intrusion. For efficient planning and management of coastal aquifers in Kayalpattu and Tiruchopuram villages, which extend over 4·05 km², it is essential to delineate and predict the extent of intrusion into the shallow aquifer. Management of ground water in coastal aquifers is composed of major elements that should be properly evaluated, and special attention is given to the sea water intrusion problem. Different data, like hydro‐geomorphological and depth‐wise iso‐apparent resistivity, are integrated spatially using a geographical information system. The stack‐unit mapping approach is used to delineate the zones with iso‐apparent resistivity of less than 10 Ω m have been found to be increasing in areal extent with reference to depth. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of Sea‐Level Rise on Salt Water Intrusion near a Coastal Well Field in Southeastern Florida

A variable‐density groundwater flow and dispersive solute transport model was developed for the shallow coastal aquifer system near a municipal supply well field in southeastern Florida. The model was calibrated for a 105‐year period (1900 to 2005). An analysis with the model suggests that well‐field withdrawals were the dominant cause of salt water intrusion near the well field, and that historical sea‐level rise, which is similar to lower‐bound projections of future sea‐level rise, exacerbated the extent of salt water intrusion. Average 2005 hydrologic conditions were used for 100‐year sensitivity simulations aimed at quantifying the effect of projected rises in sea level on fresh coastal groundwater resources near the well field. Use of average 2005 hydrologic conditions and a constant sea level result in total dissolved solids (TDS) concentration of the well field exceeding drinking water standards after 70 years. When sea‐level rise is included in the simulations, drinking water standards are exceeded 10 to 21 years earlier, depending on the specified rate of sea‐level rise.     

Remediation of Sea Water Intrusion: A Case Study

Sea water intrusion and remediation in the Upper Floridan Aquifer in South Carolina is simulated using the finite-element model SUTRA developed by the U.S. Geological Survey. A sensitivity analysis of the effect of the hydrogeologic parameters on the sea water recharge and seepage velocities is performed. An increase in confining unit and/or in aquifer conductivity results in an increase of the sea water recharge. An increase in aquifer porosity results in a decrease of the sea water recharge. Among the three remedial techniques simulated—reduced aquifer withdrawals, an injection well, and a combined injection and capture well—the reduced aquifer withdrawals and injection well are the best methods for preventing sea water intrusion.     

Modeling Coastal Salinity in Quasi 2D and 3D Using a DUALEM‐421 and Inversion Software

Rising sea levels, owing to climate change, are a threat to fresh water coastal aquifers. This is because saline intrusions are caused by increases and intensification of medium‐large scale influences including sea level rise, wave climate, tidal cycles, and shifts in beach morphology. Methods are therefore required to understand the dynamics of these interactions. While traditional borehole and galvanic contact resistivity (GCR) techniques have been successful they are time‐consuming. Alternatively, frequency‐domain electromagnetic (FEM) induction is potentially useful as physical contact with the ground is not required. A DUALEM‐421 and EM4Soil inversion software package are used to develop a quasi two‐ (2D) and quasi three‐dimensional (3D) electromagnetic conductivity images (EMCI) across Long Reef Beach located north of Sydney Harbour, New South Wales, Australia. The quasi 2D models discern: the dry sand (<10 mS/m) associated with the incipient dune; sand with fresh water (10 to 20 mS/m); mixing of fresh and saline water (20 to 500 mS/m), and; saline sand of varying moisture (more than 500 mS/m). The quasi 3D EMCIs generated for low and high tides suggest that daily tidal cycles do not have a significant effect on local groundwater salinity. Instead, the saline intrusion is most likely influenced by medium‐large scale drivers including local wave climate and morphology along this wave‐dominated beach. Further research is required to elucidate the influence of spring‐neap tidal cycles, contrasting beach morphological states and sea level rise.     

Coastal Resiliency Groundwater Considerations

The potential for rising groundwater is an important consideration in any coastal resiliency assessment. Unlike other groundwater modeling that focuses mostly on contaminant tracking, coastal groundwater resiliency assessments are primarily concerned with the potential for groundwater emergence induced by sea level rise. This provides more options for modelers that range from simplified water table elevation models to fully integrated groundwater and storm water models. The selection is dependent on available data and project needs. However, despite the relative simplicity of some of the techniques, all the methods benefit from a professional with hydrogeological training.     

Vulnerability indicators of sea water intrusion

In this paper, simple indicators of the propensity for sea water intrusion (SWI) to occur (referred to as "SWI vulnerability indicators") are devised. The analysis is based on an existing analytical solution for the steady-state position of a sharp fresh water-salt water interface. Interface characteristics, that is, the wedge toe location and sea water volume, are used in quantifying SWI in both confined and unconfined aquifers. Rates-of-change (partial derivatives of the analytical solution) in the wedge toe or sea water volume are used to quantify the aquifer vulnerability to various stress situations, including (1) sea-level rise; (2) change in recharge (e.g., due to climate change); and (3) change in seaward discharge. A selection of coastal aquifer cases is used to apply the SWI vulnerability indicators, and the proposed methodology produces interpretations of SWI vulnerability that are broadly consistent with more comprehensive investigations. Several inferences regarding SWI vulnerability arise from the analysis, including: (1) sea-level rise impacts are more extensive in aquifers with head-controlled rather than flux-controlled inland boundaries, whereas the opposite is true for recharge change impacts; (2) sea-level rise does not induce SWI in constant-discharge confined aquifers; (3) SWI vulnerability varies depending on the causal factor, and therefore vulnerability composites are needed that differentiate vulnerability to such threats as sea-level rise, climate change, and changes in seaward groundwater discharge. We contend that the approach is an improvement over existing methods for characterizing SWI vulnerability, because the method has theoretical underpinnings and yet calculations are simple, although the coastal aquifer conceptualization is highly idealized.     

Application of an Analytical Solution as a Screening Tool for Sea Water Intrusion

Sea water intrusion into aquifers is problematic in many coastal areas. The physics and chemistry of this issue are complex, and sea water intrusion remains challenging to quantify. Simple assessment tools like analytical models offer advantages of rapid application, but their applicability to field situations is unclear. This study examines the reliability of a popular sharp‐interface analytical approach for estimating the extent of sea water in a homogeneous coastal aquifer subjected to pumping and regional flow effects and under steady‐state conditions. The analytical model is tested against observations from Canada, the United States, and Australia to assess its utility as an initial approximation of sea water extent for the purposes of rapid groundwater management decision making. The occurrence of sea water intrusion resulting in increased salinity at pumping wells was correctly predicted in approximately 60% of cases. Application of a correction to account for dispersion did not markedly improve the results. Failure of the analytical model to provide correct predictions can be attributed to mismatches between its simplifying assumptions and more complex field settings. The best results occurred where the toe of the salt water wedge is expected to be the closest to the coast under predevelopment conditions. Predictions were the poorest for aquifers where the salt water wedge was expected to extend further inland under predevelopment conditions and was therefore more dispersive prior to pumping. Sharp‐interface solutions remain useful tools to screen for the vulnerability of coastal aquifers to sea water intrusion, although the significant sources of uncertainty identified in this study require careful consideration to avoid misinterpreting sharp‐interface results.     

Numerical Simulations of Seasonally Oscillated Groundwater Dynamics in Coastal Confined Aquifers

Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10⁻⁴ m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10⁻⁴ m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10⁻⁴ m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10⁻⁴ m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.     

Boundary condition effects on maximum groundwater withdrawal in coastal aquifers

Prevention of sea water intrusion in coastal aquifers subject to groundwater withdrawal requires optimization of well pumping rates to maximize the water supply while avoiding sea water intrusion. Boundary conditions and the aquifer domain size have significant influences on simulating flow and concentration fields and estimating maximum pumping rates. In this study, an analytical solution is derived based on the potential-flow theory for evaluating maximum groundwater pumping rates in a domain with a constant hydraulic head landward boundary. An empirical correction factor, which was introduced by Pool and Carrera (2011) to account for mixing in the case with a constant recharge rate boundary condition, is found also applicable for the case with a constant hydraulic head boundary condition, and therefore greatly improves the usefulness of the sharp-interface analytical solution. Comparing with the solution for a constant recharge rate boundary, we find that a constant hydraulic head boundary often yields larger estimations of the maximum pumping rate and when the domain size is five times greater than the distance between the well and the coastline, the effect of setting different landward boundary conditions becomes insignificant with a relative difference between two solutions less than 2.5%. These findings can serve as a preliminary guidance for conducting numerical simulations and designing tank-scale laboratory experiments for studying groundwater withdrawal problems in coastal aquifers with minimized boundary condition effects.     

Effect of climate change on sea water intrusion in coastal aquifers

There is increasing debate these days on climate change and its possible consequences. Much of this debate has focused in the context of surface water systems. In many arid areas of the world, rainfall is scarce and so is surface runoff. These areas rely heavily on groundwater. The consequences of climate change on groundwater are long term and can be far reaching. One of the more apparent consequences is the increased migration of salt water inland in coastal aquifers. Using two coastal aquifers, one in Egypt and the other in India, this study investigates the effect of likely climate change on sea water intrusion. Three realistic scenarios mimicking climate change are considered. Under these scenarios, the Nile Delta aquifer is found to be more vulnerable to climate change and sea level rise. Copyright © 1999 John Wiley & Sons, Ltd.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Does sea-level rise have an impact on saltwater intrusion?

Climate change effects are expected to substantially raise the average sea level. It is widely assumed that this raise will have a severe adverse impact on saltwater intrusion processes in coastal aquifers. In this study we hypothesize that a natural mechanism, identified here as the “lifting process,” has the potential to mitigate, or in some cases completely reverse, the adverse intrusion effects induced by sea-level rise. A detailed numerical study using the MODFLOW-family computer code SEAWAT was completed to test this hypothesis and to understand the effects of this lifting process in both confined and unconfined systems. Our conceptual simulation results show that if the ambient recharge remains constant, the sea-level rise will have no long-term impact (i.e., it will not affect the steady-state salt wedge) on confined aquifers. Our transient confined-flow simulations show a self-reversal mechanism where the wedge which will initially intrude into the formation due to the sea-level rise would be naturally driven back to the original position. In unconfined systems, the lifting process would have a lesser influence due to changes in the value of effective transmissivity. A detailed sensitivity analysis was also completed to understand the sensitivity of this self-reversal effect to various aquifer parameters.     

The Interface Configuration of the Fresh‐/Dead Sea Water – Theory and Measurements

The high‐density Dead Sea water (1.235 g/cm³) forms a special interface configuration with the fresh groundwater resources of its surrounding aquifers. The fresh groundwater column beneath its surroundings is around one tenth of its length compared to oceanic water. This fact alone indicates the vulnerability of the fresh groundwater resources to the impacts of changes in the Dead Sea level and to saltwater migration. Ghyben‐Herzberg and Glover equations were used to calculate the volumes of water in coastal aquifers which were replaced by freshwater due to the interface seaward migration as a result of the drop in the level of the Dead Sea. For that purpose, the dynamic equation of Glover approach has been integrated to accommodate that type of interface readjustment. The calculated amounts of freshwater which substituted salt Dead Sea water due to the migration of interface are 3.21 · 10¹¹ m³, from a Dead Sea level of –392 m to τ411 m below sea level. The average porosity of coastal aquifers was calculated to range from 2.8 to 2.94%. Geoelectric sounding measurements showed that areas underlying the coastal aquifers formerly occupied by the Dead Sea water are gradually becoming flushed and occupied by freshwater. The latter is becoming salinized due to the residuals of Dead Sea water in the aquifer matrix, the present salinity of which is lower than that of the Dead Sea water. At the same time salt dissolution from the Lisan Marl formation is causing collapses along the shorelines in the form of sinkholes, tens of meters in diameter and depth.     

Relative role and extent of marine and groundwater inundation on a dune‐dominated barrier island under sea‐level rise scenarios

Climate change and sea‐level rise will have severe impacts on coastal water resources around the world. However, whereas the influence of marine inundation is well documented in the literature, the impact of groundwater inundation on coastal communities is not well known. Here, core analysis, groundwater monitoring, and ground penetrating radar are utilized to assess the groundwater regime of the surficial aquifer on Bogue Banks Barrier Island (USA). Then, geospatial techniques are used to assess the relative roles and extents of groundwater and marine inundation on the dune‐dominated barrier island under sea‐level rise scenarios of 0.2, 0.5, and 1.0 m above current conditions by 2100. Additionally, the effects of rising water tables on onsite wastewater treatment systems (OWTS) are modelled using the projected sea‐level rise scenarios. The results indicate that the surficial aquifer comprising fine to medium sands responds quickly to precipitation. Water‐level measurements reveal varying thicknesses of the vadose zone (>3 to 0 m) and several groundwater mounds with radial flow patterns. Results from projected sea‐level rise scenarios suggest that owing to aquifer properties and morphology of the island, groundwater inundation may occur at the same rate as marine inundation. Furthermore, the area inundated by groundwater may be as significant as that affected by marine inundation. The results also show that the proportion of land in the study area where OWTS may be perpetually compromised by rising water tables under worst case scenarios may range from ~43 to ~54% over an 86‐year‐period. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of Recharge Wells and Flow Barriers on Seawater Intrusion

The installation of recharge wells and subsurface flow barriers are among several strategies proposed to control seawater intrusion on coastal groundwater systems. In this study, we performed laboratory‐scale experiments and numerical simulations to determine the effects of the location and application of recharge wells, and of the location and penetration depth of flow barriers, on controlling seawater intrusion in unconfined coastal aquifers. We also compared the experimental results with existing analytical solutions. Our results showed that more effective saltwater repulsion is achieved when the recharge water is injected at the toe of the saltwater wedge. Point injection yields about the same repulsion compared with line injection from a screened well for the same recharge rate. Results for flow barriers showed that more effective saltwater repulsion is achieved with deeper barrier penetration and with barriers located closer to the coast. When the flow barrier is installed inland from the original toe position however, saltwater intrusion increases with deeper barrier penetration. Saltwater repulsion due to flow barrier installation was found to be linearly related to horizontal barrier location and a polynomial function of the barrier penetration depth.     

An Excel Macro to Plot the HFE‐Diagram to Identify Sea Water Intrusion Phases

A hydrochemical facies evolution diagram (HFE‐D) is a multirectangular diagram, which is a useful tool in the interpretation of sea water intrusion processes. This method note describes a simple method for generating an HFE‐D plot using the spreadsheet software package, Microsoft Excel. The code was applied to groundwater from the alluvial coastal plain of Grosseto (Tuscany, Italy), which is characterized by a complex salinization process in which sea water mixes with sulfate or bicarbonate recharge water.