Using MODFLOW code to approach transient hydraulic head with a sharp‐interface solution

Saltwater intrusion problems have been usually tackled through analytical models because of its simplicity, easy implementation and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, which neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. This paper provides insight into the validity of a sharp‐interface approximation defined from a steady state solution when applied to transient seawater intrusion problems. The validation tests have been performed on a 3D unconfined synthetic aquifer, which include spatial and temporal distribution of recharge and pumping wells. Using a change of variable, the governing equation of the steady state sharp‐interface problem can be written with the same structure of the steady confined groundwater flow equation as a function of a single potential variable (?). We propose to approach also the transient problem solving a single potential equation (using also the ? variable) with the same structure of the confined groundwater flow equation. It will allow solving the problem by using the classical MODFLOW code. We have used the parameter estimation model PEST to calibrate the parameters of the transient sharp‐interface equation. We show how after the calibration process, the sharp‐interface approach may provide accurate enough results when applied to transient problems and improve the steady state results, thus avoiding the need of implementing a density‐dependent model and reducing the computational cost. This has been proved by comparing results with those obtained using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. The comparison was performed in terms of piezometric heads, seawater penetration, transition zone width and critical pumping rates. Copyright © 2014 John Wiley & Sons, Ltd.     

Optimal and sustainable extraction of groundwater in coastal aquifers

Four examples are investigated for the optimal and sustainable extraction of groundwater from a coastal aquifer under the threat of seawater intrusion. The objectives and constraints of these management scenarios include maximizing the total volume of water pumped, maximizing the profit of selling water, minimizing the operational and water treatment costs, minimizing the salt concentration of the pumped water, and controlling the drawdown limits. The physical model is based on the density-dependent advective-dispersive solute transport model. Genetic algorithm is used as the optimization tool. The models are tested on a hypothetical confined aquifer with four pumping wells located at various depths. These solutions establish the feasibility of simulating various management scenarios under complex three-dimensional flow and transport processes in coastal aquifers for the optimal and sustainable use of groundwater.     

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.     

Modeling Steady Sea Water Intrusion with Single‐Density Groundwater Codes

Steady interface flow in heterogeneous aquifer systems is simulated with single‐density groundwater codes by using transformed values for the hydraulic conductivity and thickness of the aquifers and aquitards. For example, unconfined interface flow may be simulated with a transformed model by setting the base of the aquifer to sea level and by multiplying the hydraulic conductivity with 41 (for sea water density of 1025 kg/m³). Similar transformations are derived for unconfined interface flow with a finite aquifer base and for confined multi‐aquifer interface flow. The head and flow distribution are identical in the transformed and original model domains. The location of the interface is obtained through application of the Ghyben‐Herzberg formula. The transformed problem may be solved with a single‐density code that is able to simulate unconfined flow where the saturated thickness is a linear function of the head and, depending on the boundary conditions, the code needs to be able to simulate dry cells where the saturated thickness is zero. For multi‐aquifer interface flow, an additional requirement is that the code must be able to handle vertical leakage in situations where flow in an aquifer is unconfined while there is also flow in the aquifer directly above it. Specific examples and limitations are discussed for the application of the approach with MODFLOW. Comparisons between exact interface flow solutions and MODFLOW solutions of the transformed model domain show good agreement. The presented approach is an efficient alternative to running transient sea water intrusion models until steady state is reached.     

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.     

Groundwater modelling and marine intrusion in the semi‐arid systems (Cap‐Bon,Tunisia)

This article studies the effect of drought and pumping discharge on groundwater supplies and marine intrusion in the Korba aquifer (Cap‐Bon peninsula, Tunisia). The Groundwater Modelling System has been used to model the groundwater flow and to simulate the seawater intrusion. The calibration is based on the groundwater levels in the steady state from 1963, and in the transient state from the groundwater levels from 1963 to 2005. The main objective is to quantify the components of the groundwater mass balance and to estimate the hydraulic conductivity distribution. The impact of pumping discharge on the groundwater level evolution has been examined by two pumping scenarios P1 (no. 8420) and P2 (no. 8862) wells. The hydrodynamic modelling shows the increasing drawdowns after 14 years of pumping: 4 m in P1 well and about 5 m in P2 well below sea level. The drawdowns are accompanied by the inverse hydraulic gradient. The numerical model was used to discuss the management of the groundwater resources of Cap‐Bon. As the population continues to grow and the demand for groundwater pumping intensifies beyond the 1963 level, it can be expected that the actual extent of seawater intrusion in the future would be more severe than the model prediction. Better strategies for groundwater development and management will be necessary to protect the freshwater aquifers to the marine intrusion. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling study on the impact of deep building foundations on the groundwater system

Coastal areas are usually the preferred place of habitation for human beings. Anthropogenic activities such as the construction of high‐rise buildings and underground transport systems usually require extensive deep foundations and ground engineering works, which may unintentionally modify the coastal groundwater system because the construction materials of foundations are usually of low hydraulic conductivity. In this paper, the impact of these building foundations on the groundwater regime is studied using hypothetical flow and transport models. Various possible realizations of foundation distributions are generated using stochastic parameters derived from a topographical map of an actual coastal area in Hong Kong. The effective hydraulic conductivity is first calculated for different realizations and the results show that the effective hydraulic conductivity can be reduced significantly. Then a hypothetical numerical model based on FEFLOW is set up to study the change of hydraulic head, groundwater discharge, and saltwater‐fresh water interface. The groundwater level and flow are modified to various degrees, depending on the foundations percentage and the distribution pattern of the buildings. When the foundations percentage is high and the building foundations are aggregated, the hydraulic head is raised significantly and the originally one‐dimensional groundwater flow field becomes complicated. Seaward groundwater discharge will be reduced and some groundwater may become seepage through the ground surface. The transport model shows that, after foundations are added, overall the seawater and fresh groundwater interface moves landward, so extensive foundations may induce seawater intrusion. It is believed that the modification of the coastal groundwater system by building foundations may have engineering and environmental implications, such as submarine groundwater discharge, foundation corrosion, and slope stability. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling of colloidal particle and heavy metal transfer behaviours during seawater intrusion and refreshing processes

The proper management of coastal aquifers commonly requires an understanding of regional mass flow and complete seawater–freshwater circulation. In this study, time series observations of seawater intrusion and refreshing were conducted using a column experiment based on natural flow conditions in coastal groundwater and a sampled medium from a coastal sandy aquifer without chemical treatment. Ranges of hydrodynamic and hydrochemical variables were tested and analysed. The results showed that the zeta potential of suspended colloids in aqueous solution in an aquifer polluted with 0.5 g/kg of heavy metals exhibited an isoelectric point for pH values ranging from 5.70 to 6.07 when freshwater or seawater completely occupied the aquifer pores, which is representative of natural hydrochemical conditions. In this scenario, a high background concentration of heavy metals induced colloidal immobilization. Otherwise, seawater–freshwater circulation enabled colloid mobilization due to ionic strength and pH fluctuations. The migration of multiple heavy metals occurred at a characteristic time of approximately 1 pore volume after each intrusion stage began and when the peak rate of colloid release was reached. At these times, the colloid behaviour determined the quantity and pathway of heavy metal transport. On the basis of the influences of seawater and freshwater interactions, the quantity of mobilized particles generally decreased and was uniformly distributed in each fraction due to particle loss and decreased porous connectivity. We speculate that the decrease in the total surface area of the migratory colloids may cause colloid‐associated heavy metal transport to decrease. The experimental results provide a useful basis for testing coastal groundwater flow and mass transport models because these phenomena require full characterization to precisely evaluate the associated fluxes from the field scale to the microscopic dimension.     

Estimation of the inland extending length of the freshwater–saltwater interface in coastal unconfined aquifers

ABSTRACT

The inland extending length of the freshwater–saltwater interface toe is useful in studies of seawater intrusion in coastal areas. The submarine fresh groundwater discharge in coastal zones is affected not only by hydraulic conductivity and hydraulic gradient of the aquifer, but also by the position of the interface. Two observation wells at different distances from the coast are required to calculate the fresh groundwater flow rate in coastal unconfined aquifers. By considering that the submarine groundwater discharge is equal to the groundwater flow rate, the length of the interface toe extending inland can be estimated when the groundwater flow is at a steady-flow state. Aquifers with horizontal and sloping confined beds and without/with unique surface vertical infiltration are considered. Examples used to illustrate the application of these methods indicate that the inland extending lengths of the interface toe in aquifers with vertical surface infiltration are much shorter than those in aquifers without vertical surface infiltration, and the length of the interface in aquifers with a horizontal confining lower bed are smaller than those in aquifers with a confining lower bed sloping towards the sea. The extent of the interface on the northwestern coast near the city of Beihai in southern Guangxi, China, on 18 January 2013 was estimated as 471–478 m.

Editor M.C. Acreman Associate editor not assigned     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

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.     

Understanding and modelling spatial drain–aquifer interactions in a low‐lying coastal aquifer—the Thurne catchment,Norfolk, UK

Seawater intrusion into fresh groundwater formations generally results inadvertently from human activities, such as over‐abstraction from coastal aquifers. This article describes the data analysis to quantify drain–aquifer interactions in a low‐lying pump‐drained coastal aquifer, which is subject to saline intrusion due to widespread land drainage, and the resulting development and application of a numerical groundwater model to understand the spatial groundwater system behaviour (including groundwater salinity fluxes). Without measured flow data in this pump‐drained catchment, a novel groundwater head‐dependent approach to hydrograph separation is described. Time‐variant and time‐invariant MODFLOW analyses are utilised to examine the flow processes. A new approach to calculate drain coefficients, which represent the extensive network of drainage ditches in the regional model, using field information, is described; the sum of the drainage coefficients are close to the values independently estimated from the head‐dependent hydrograph separation. Results show that (1) the groundwater flows into the drainage systems are well reproduced using the new drain coefficients, (2) particle tracking of fresh and saline water can explain observed spatial salinity distribution within drainage networks and (3) the modelled flow of seawater across the coast is approximately 25% greater than that discharged by the pumps, demonstrating the need for drainage management to be aware of the slow response of groundwater systems to past drainage system changes. The article demonstrates that numerical groundwater modelling can produce the improved understanding needed to inform management decisions in such complex environments. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

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.     

Saltwater intrusion estimation in a karstified coastal system using density-dependent modelling and comparison with the sharp-interface approach

Abstract

Saltwater intrusion is a naturally occurring phenomenon that is exacerbated significantly by excessive groundwater exploitation in coastal aquifers. In order to determine the extent of saltwater intrusion in a karstified coastal aquifer in Crete, Greece, a three-dimensional, density-dependent groundwater flow and transport model was developed and compared to the more traditional sharp-interface approach. The karstified medium was modelled using a combination of the equivalent porous medium approach (for lower-order fractures) and a discrete fracture approach (for the main fractures/faults). The model takes into consideration the geomorphologic characteristics of the karstic system, such as the depth and orientation of the fault network, and the diffusion phenomena associated with the variable densities of freshwater and saltwater—parameters that create a complex system, inducing uncertainty in the model. The model results showed that the orientation of the fractures, the pumping activity and the fluid density effects drive the seawater intrusion front asymmetrically inland.

Editor Z.W. Kundzewicz

Citation Dokou, Z. and Karatzas, G.P., 2012. Saltwater intrusion estimation in a karstified coastal system using density-dependent modelling and comparison with the sharp-interface approach. Hydrological Sciences Journal, 57 (5), 985–999.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

Brief Review on Methods of Estimation of the Location of a Fresh Water–Salt Water Interface with Hydraulic Heads or Pressures in Coastal Zones

Methods of estimation of the location of a sharp fresh water–salt water interface with hydraulic heads or pressures are relatively simple and are widely used. Progress has been made in the recent decade toward the mathematical relations describing the position of the sharp interface using hydraulic heads or pressures in coastal zones. This paper reviews several methods for estimation of the location of fresh water–salt water interface in coastal aquifers, including the classical Ghyben–Herzberg relation. The location of the fresh water–salt water interface in a coastal homogeneous, isotropic unconfined aquifer can be estimated based on piezometric heads at two points in the same vertical line tapping, respectively, the salt water zone (including the interface) and the fresh water zone (from the water table to the interface) when the groundwater flow system is in a steady state and satisfies the Dupuit assumption. If pressures are measured at two points in the fresh water and salt water zones in the same vertical line in the coastal aquifer under the same assumption, then the position of the interface can still be estimated with the pressure data. If the Dupuit assumption is not met in coastal aquifers and the vertical fresh water head gradients can be approximated with a straight line, the position of the interface can roughly be estimated by using the water level data in a partially penetrating well during drilling of the well.     

Estimation of seawater retreat timescales in homogeneous and confined coastal aquifers based on dimensional analysis

Reliable estimation of the timescales of seawater retreat (SWR) in coastal aquifers has significant implications for sustainable groundwater management in coastal areas. Due to the complexity of coastal dynamics, analytical estimates of SWR timescales are limited in number, and numerical models are mainly utilized. Although numerical models allow for accurate analysis of the coastal systems, they concern particular aquifer cases, thus their results are not generally applicable. Here, we perform a dimensional analysis for the depth-integrated freshwater/seawater continuity equations and identify the dimensionless parameters that affect SWR in confined, homogeneous and flux-controlled aquifers. Based on numerical simulations, we gain insight into how these parameters affect SWR, and produce type curves which can be used as a simple tool for estimating SWR timescales in idealized coastal systems, independent of the hydraulic and geometric characteristics of the aquifer. The reliability of our estimates in real-world applications is also discussed.     

Numerical study of salinity variation in a coastal aquifer: a case study of the Motooka region in western Japan

The Motooka region in the Fukuoka prefecture in western Japan is a coastal area, where groundwater is utilized as the main water resource for greenhouse agriculture and domestic use. Over-exploitation of groundwater has resulted in seawater intrusion and thus in the contamination of the freshwater aquifer. Fluctuations in electric conductivities caused by such intrusion are a crucial problem, since even slight changes in electric conductivities of the water used for agricultural purposes significantly affect the crops’ growth and yield. However, no study has thus far been conducted on the electric conductivity fluctuations caused by groundwater pumping and seasonal recharge of groundwater in the Motooka region. To this end, an attempt is made in the present study to develop a numerical variable-density solute transport model and then apply it to simulate the electric conductivity fluctuations with groundwater pumping and rainwater recharge. This model is developed under a finite difference scheme, and the method of characteristics is used as the numerical technique for solving the advection term of the advection–dispersion solute transport equation. The results from this numerical model are compared with the field measurements.     

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.