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.     

Reducing Groundwater Pumping Cost Using Alternate Pulsed Pumping

This paper deals with water pumping cost minimization, in a confined infinite aquifer, proposing an alternate pulsed pumping schedule. The transient flow analysis is conducted for two wells with equal pumping rates. Specifically, two pumping schedules are analytically compared. In the first case, well users pump simultaneously, and in the second one they cooperate so that they pump alternately. This paper proves that the proposed alternate pumping schedule works as a stabilizer, reducing the high hydraulic drawdowns values, regardless of the aquifer characteristics. Moreover, pumping alternately is better in terms of pumping cost, compared to simultaneous pumping, though benefit become negligible as distance between wells becomes large. Two simplified equations are proposed, one to find the difference of the hydraulic drawdowns between the two pumping schedules and the other one to find the economic benefit of each well from cooperation. The equations are finally applied to a number of cases and their results are compared to the results obtained from an algorithm created to calculate the hydraulic drawdowns and the pumping cost, using the Theis equation. The results can be very useful in irrigation scheduling, as they can be applied to systems of well users/farmers, to reduce pumping cost.     

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.     

Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer,Cyprus

Abstract

We investigate the general methodology for an intensive development of coastal aquifers, described in a companion paper, through its application to the management of the Akrotiri aquifer, Cyprus. The Zakaki area of that aquifer, adjacent to Lemessos City, is managed such that it permits a fixed annual agricultural water demand to be met, as well as and a fraction of the water demand of Lemessos, which varies according to available surface water. Effluents of the Lemessos wastewater treatment plant are injected into the aquifer to counteract the seawater intrusion resulting from the increased pumping. The locations of pumping and injection wells are optimized based on least-cost, subject to meeting the demand. This strategy controls sea intrusion so effectively that desalting of only small volumes of slightly brackish groundwater is required over short times, while ～2.3 m³ of groundwater is produced for each 1 m³ of injected treated wastewater. The cost over the 20-year period 2000–2020 of operation is ～40 M€ and the unit production cost of potable water is under 0.2 €/m³. The comparison between the deterministic and stochastic analyses of the groundwater dynamics indicates the former as conservative, i.e. yielding higher groundwater salinity at the well. The Akrotiri case study shows that the proposed aquifer management scheme yields solutions that are preferable to the widely promoted seawater desalination, also considering the revenues from using the treated wastewater for irrigation.

Citation Koussis, A. D., Georgopoulou, E., Kotronarou, A., Mazi, K., Restrepo, P., Destouni, G., Prieto, C., Rodriguez, J. J., Rodriguez-Mirasol, J., Cordero, T., Ioannou, C., Georgiou, A., Schwartz, J. & Zacharias, I. (2010) Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer, Cyprus. Hydrol. Sci. J. 55(7), 1234–1245.     

Analysis of Nitrate-Nitrogen Movement Near High-Capacity Irrigation Wells

AGalerkin finite-element model coupled with a particle tracking routine was developed to analyze the flow and transport dynamics near a high-capacity irrigation well. The model was used to compute the head distribution around the pumping well, to determine the area of influence, and to define ground water flowlines during short-term pumping periods typical of those used to collect water quality samples from high-capacity wells. In addition to hypothetical example results, the model was used to qualitatively analyze data obtained from pump-and-sample experiments conducted in an unconfined alluvial aquifer within the Platte River valley of south-central Nebraska where nitrate-nitrogen (NO₃-N) contamination is prevalent.
Simulation results of both the hypothetical and field cases suggest that short-term pumping events, impact a limited volume of aquifer. The area of influence and flowlines are affected by aquifer anisotropy, pumping rate, and well construction characteristics). Ground water above or below the screened intervals does not enter a partially penetrating well in anisotropic aquifers. In aquifers where NO₃-N concentration varies vertically and horizontally, waler quality samples from an irrigation, or other high-capacity, well provide only limited information about ground water contamination. A numerical model is thus recommended for calculating the area of influence and determining flowlines around high-capacity wells so that information derived from water quality samples collected at the wellhead can be better interpreted.     

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.     

Estimating hydraulic properties using a moving-model approach and multiple aquifer tests

A new method was developed for characterizing geohydrologic columns that extended >600 m deep at sites with as many as six discrete aquifers. This method was applied at 12 sites within the Southwest Florida Water Management District. Sites typically were equipped with multiple production wells, one for each aquifer and one or more observation wells per aquifer. The average hydraulic properties of the aquifers and confining units within radii of 30 to >300 m were characterized at each site. Aquifers were pumped individually and water levels were monitored in stressed and adjacent aquifers during each pumping event. Drawdowns at a site were interpreted using a radial numerical model that extended from land surface to the base of the geohydrologic column and simulated all pumping events. Conceptually, the radial model moves between stress periods and recenters on the production well during each test. Hydraulic conductivity was assumed homogeneous and isotropic within each aquifer and confining unit. Hydraulic property estimates for all of the aquifers and confining units were consistent and reasonable because results from multiple aquifers and pumping events were analyzed simultaneously.     

Interpretation of groundwater level monitoring results in karst aquifers: examples from the Dinaric karst

The paper presents an attempt to determine the characteristics of karst aquifers using information on groundwater level (GWL) in natural holes and boreholes with different data quantity and time resolution of GWL measurements. In this paper the particulars of karst aquifers were analysed for four examples from the Dinaric karst. In all four study areas, aquifers are formed in bare, deep and well‐developed Dinaric karst consisting of Cretaceous limestones. The first example represents a wide area of Imotsko polje in the karst. The aquifer was analysed on the basis of infrequent water level monitoring in natural karst water features (jamas, lakes, wells) and discharges of springs and rivers. The karst aquifer in this example is complex, non‐homogenous and variable in space and time, which is frequent in the Dinaric karst. Regardless of the aforementioned it was possible to determine its elementary characteristics. The second example represents 10 wells used for the water supply for the city of Pula. The GWL and salinity were measured once a week in the period between 1981 and 1996. Even though these measurements were relatively infrequent in space and time, they served as bases for assessment of average and maximum aquifer conditions as well as boundaries of saltwater intrusion. In the third example only a portion of aquifer of the karst spring Blaz, which is in the contact with the Adriatic Seas, has been analyzed. It is a spring with an intrusion of salt water. For purposes of study of saltwater intrusion, 26 piezometers were drilled in its vicinity in which GWL, salinity and temperature were measured once a day during 168 days, a period comprising one complete cycle of seawater intrusion and retreat. These measurements proved the existence of dispersed discharge from the aquifer into the sea and its non‐homogeneity in space. In the fourth example GWL was measured continuously in 10 deep (up to 300 m) piezometers in the hinterland of the Ombla Spring catchment. The measurement period lasted 2 years (January 1988 to December 1989). The analyses are made with hourly data. The results made it possible to determine numerous characteristics of the karst aquifer and a significant non‐homogeneity of groundwater distribution in karst aquifers, depending more on the underground karst phenomena than the surface karst forms. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

Geometry and dynamics of the freshwater—seawater interface in a coastal aquifer in southeastern Spain

Abstract

The contact between freshwater and seawater in coastal aquifers is studied using a relatively simple model for homogeneous aquifers. However, for real aquifers it is not so simple. The desalination plant built to supply water to the city of Almería is situated over the aquifer in the southern part of the River Andarax Delta. Its design capacity is 1100 L s^?1, and it is supplied from boreholes pumping water from beneath the freshwater—seawater contact in this aquifer. Well logs kept over a period of two years have allowed us to accurately define the interface geometry of the freshwater—seawater contact. Lithological data collected from 31 boreholes have also indicated the existence of strata with low hydraulic conductivity, within others of high conductivity. During a simultaneous pumping test of six wells with 690 L s^?1 total discharge, electrical conductivity measurements showed the influx of seawater 6–10 m below sea level and a drawdown of the interface in the piezometers closest to the pumping wells.     

The unusual and large drawdown response of buried-valley aquifers to pumping

The buried-valley aquifers that are common in the glacial deposits of the northern hemisphere are a typical case of the strip aquifers that occur in many parts of the world. Pumping from a narrow strip aquifer leads to much greater drawdown and much more distant drawdown effects then would occur in a sheet aquifer with a similar transmissivity and storage coefficient. Widely used theories for radial flow to wells, such as the Theis equation, are not appropriate for narrow strip aquifers. Previously published theory for flow to wells in semiconfined strip aquifers is reviewed and a practical format of the type curves for pumping-test analysis is described. The drawdown response of strip aquifers to pumping tests is distinctive, especially for observation wells near the pumped well. A case study is presented, based on extensive pumping test experience for the Estevan Valley Aquifer in southern Saskatchewan, Canada. Evaluation of groundwater resources in such buried-valley aquifers needs to take into account the unusually large drawdowns in response to pumping.     

Analytical Design Curves to Maximize Pumping or Minimize Injection in Coastal Aquifers

Explicit algebraic equations are derived to determine approximate maximum pumping rates or minimum injection rates to limit sea water intrusion to a prespecified distance from the coastline. The equations are based on Strack's (1976) single-potential solution. The maximum pumping rates and minimum injection rates applied at wells with uniform spacing to control the inland movement of the fresh water-salt water interface in a coastal aquifer could be calculated from Strack's (1976) solution without the need of a numerical optimization algorithm. When wells are distributed in a simple fashion, the maximum intrusion location can be identified precisely for pumping cases and approximately for injection cases. For pumping cases, critical points are the limit of allowable salt water intrusion, whereas no such limit exists for injection cases. Once an application site is identified, a series of design curves for pumping and injection rates can be developed for arbitrary intrusion limits. When a user is interested only in the largest pumping rates associated with critical points, one design curve can yield complete information.     

Closed-form analytical solutions incorporating pumping and tidal effects in various coastal aquifer systems

Pumping wells are common in coastal aquifers affected by tides. Here we present analytical solutions of groundwater table or head variations during a constant rate pumping from a single, fully-penetrating well in coastal aquifer systems comprising an unconfined aquifer, a confined aquifer and semi-permeable layer between them. The unconfined aquifer terminates at the coastline (or river bank) and the other two layers extend under tidal water (sea or tidal river) for a certain distance L. Analytical solutions are derived for 11 reasonable combinations of different situations of the L-value (zero, finite, and infinite), of the middle layer’s permeability (semi-permeable and impermeable), of the boundary condition at the aquifer’s submarine terminal (Dirichlet describing direct connection with seawater and no-flow describing the existence of an impermeable capping), and of the tidal water body (sea and tidal river). Solutions are discussed with application examples in fitting field observations and parameter estimations.     

Maximizing Net Extraction Using an Injection‐Extraction Well Pair in a Coastal Aquifer

In this study, we examine the maximum net extraction rate from the novel arrangement of an injection‐extraction well pair in a coastal aquifer, where fresh groundwater is reinjected through the injection well located between the interface toe and extraction well. Complex potential theory is employed to derive a new analytical solution for the maximum net extraction rate and corresponding stagnation‐point locations and recirculation ratio, assuming steady‐state, sharp‐interface conditions. The injection‐extraction well‐pair system outperforms a traditional single extraction well in terms of net extraction rate for a broad range of well placement and pumping rates, which is up to 50% higher for an aquifer with a thickness of 20 m, hydraulic conductivity of 10 m/d, and fresh water influx of 0.24 m²/d. Sensitivity analyses show that for a given fresh water discharge from an inland aquifer, a larger maximum net extraction is expected in cases with a smaller hydraulic conductivity or a smaller aquifer thickness, notwithstanding physical limits to drawdown at the pumping well that are not considered here. For an extraction well with a fixed location, the optimal net extraction rate linearly increases with the distance between the injection well and the sea, and the corresponding injection rate and recirculation ratio also increase. The analytical analysis in this study provides initial guidance for the design of well‐pair systems in coastal aquifers, and is therefore an extension beyond previous applications of analytical solutions of coastal pumping that apply only to extraction or injection wells.     

An Application of HFE-D for Evaluating Sea Water Intrusion in Coastal Aquifers of Southern Vietnam

The coastal aquifers and inland waters of the Long Xuyen Quadrangle and Ca Mau Peninsula of southern Vietnam have been significantly impacted by sea water intrusion (SI) as a result of recent anthropogenic activities. This study identified the evolution and spatial distribution of hydrochemical conditions in coastal aquifers at this region using Hydrochemical Facies Evolution Diagram (HFE-D) and Geographical Information System mapping. Hydraulic heads and water chemistry were measured at 31 observation wells in four layered aquifers during dry and rainy seasons in early (2005), and more recent (2016), stages of agricultural development. Hydrochemical facies associated with intrusion or freshening stages were mapped in each aquifer after assigning mixing index values to each facies. The position of groundwater freshening and SI phases differed in Holocene, Upper Pleistocene, Middle Pleistocene, and Lower Pleistocene aquifers. The geographic position of freshening and intrusion fronts differ in dry and rainy seasons, and shifted after 11 years of groundwater abstraction in all four aquifers. The spatial and temporal differences in hydrochemical facies distributions according to HFE-D reflect the relative impact of SI in the four aquifers. The study results provide a better understanding of the evolution of groundwater quality associated with SI in a peninsular coastal aquifer system, and highlight the need for improving groundwater quality and management in similar coastal regions.     

Locating the zone of saline intrusion in a coastal karst aquifer using springflow data

Coastal fresh water aquifers are an increasingly desirable resource. In a karstic aquifer, sea water intrusion occurs as a salt water wedge, like in porous media. However, preferential flow conduits may alter the spatial and temporal distribution of the salt water. This is typically the case when the outlet of the aquifer is a brackish spring. This paper shows that salinity and flow rate variations at a spring, where salinity is inversely proportional to discharge, can help to understand the hydrodynamic functioning of the aquifer and to locate the fresh water-sea water mixing zone deep inside the aquifer. The volume of water-filled conduit between the sea water intrusion zone and the spring outlet is calculated by the integral over time of the flow rate during the time lag between the flow rate increase and the salinity decrease as measured at the spring. In the example of the spring at Almyros of Heraklio (Crete, Greece), this time lag is variable, depending on the discharge, but the volume of water-filled conduit appears to be constant, which shows that the processes of salt water intrusion and mixing in the conduit are constant throughout the year. The distance between the spring and the zone where sea water enters the conduit is estimated and provides an indication of the position where only fresh water is present in the conduit.     

Poroelastic Effects on Fracture Characterization

Reverse water‐level fluctuations have been widely observed in aquitards or aquifers separated from a pumped confined aquifer (Noordbergum effect) immediately after the initiation of pumping. This same reverse fluctuation has been observed in a fractured crystalline‐rock aquifer at the Coles Hill uranium site in Virginia in which the reverse water‐level response occurs within a pumped fracture and results from an instantaneous strain response to pumping that precedes the pore‐pressure response in observation wells of sufficient distance from the pumped well. This response is referred to as the Mandel‐Cryer effect. The unique aspect of this water level rise during a controlled 24 h pumping test was that the reverse water levels lasted for approximately 100 min and reached a magnitude of nearly 1 cm prior to a typical drawdown response. The duration and magnitude of the response reflects the poromechanical properties of the fractured host rock and hydraulic properties of the pumped fracture. An axisymmetric flow and deformation model were developed using Biot2 in an effort to simulate the observed water‐level response along an assumed 0.5 to 1.0 cm aperture horizontal fracture 176 m from the pumping well and to identify the importance of the poroelastic effect. Results indicate that traditional aquifer‐testing methods that ignore the poromechanical response are not significantly different than results that include the response. However, the poroelastic effect allows for more accurate and efficient parameter calibration.     

Efficiency of joint use of MRS and VES to characterize coastal aquifer in Myanmar

The productivity and the water quality of coastal aquifers can be highly heterogeneous in a complex environment. The characterization of these aquifers can be improved by hydrogeological and complementary geophysical surveys. Such an integrated approach is developed in a non-consolidated coastal aquifer in Myanmar (previously named Burma).A preliminary hydrogeological survey is conducted to know better the targeted aquifers. Then, 25 sites are selected to characterize aquifers through borehole drillings and pumping tests implementation. In the same sites, magnetic resonance soundings (MRS) and vertical electrical soundings (VES) are carried out. Geophysical results are compared to hydrogeological data, and geophysical parameters are used to characterize aquifers using conversion equations. Finally, combining the analysis of technical and economical impacts of geophysics, a methodology is proposed to characterize non-consolidated coastal aquifers.Depth and thickness of saturated zone is determined by means of MRS in 68% of the sites (evaluated with 34 soundings). The average accuracy of confined storativity estimated with MRS is ± 6% (evaluated over 7 pumping tests) whereas the average accuracy of transmissivity estimation with MRS is ± 45% (evaluated using 15 pumping tests). To reduce uncertainty in VES interpretation, the aquifer geometry estimated with MRS is used as a fixed parameter in VES inversion. The accuracy of groundwater electrical conductivity evaluation from 15 VES is enough to estimate the risk of water salinity. In addition, the maximum depth of penetration of the MRS depends on the rocks' electrical resistivity and is between 20 and 80 m at the study area.     

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.