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.     

Determination of aquifer parameters based on measurements of tidal effects on a coastal aquifer near Beihai,China

Measurements of the tide and groundwater levels in coastal zones are of importance in determining the properties of coastal aquifers. The solution to a one‐dimensional unsteady groundwater flow model in a coastal confined aquifer with sinusoidal fluctuation of the tide shows that the tidal efficiency decreases exponentially with distance and the time lag increases linearly with distance from the coast. The aquifer property described by the ratio of storage coefficient to transmissivity is determined if the damping constant of the tidal efficiency or the slope of the time lag with distance are obtained on the basis of tidal measurements. Hourly observations of the tide and groundwater levels at 10 wells on the northern coast near Beihai, China show that with distance from the coast, tidal efficiency decreases roughly exponentially and the time lag increases roughly linearly. The estimated ratio of storage coefficient to transmissivity of the confined aquifer ranges from 1·169 × 10^?6 d m^?2 to 1·83 × 10^?7 d m^?2. For a given transmissivity of 750 m² d^?1, the storage coefficient of the aquifer is 8·7675 × 10^?4 with the tidal efficiency method and 1·3725 × 10^?4 with the time lag method. The damping constant of the tidal efficiency with distance can be defined as the tidal propagation coefficient. The value of the confined aquifer is determined as 0·0018892 m^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Determination of longitudinal dispersivity in an unconfined sandy aquifer

Study of the mobility of contaminants in an aquifer is an important issue for the proper remediation of contaminated groundwater. Determination of associated solute transport parameters therefore is essential for investigation of the extent to which groundwater can be contaminated. This study aimed at determining solute transport parameters for an unconfined sandy aquifer at a laboratory scale through various tracer tests using a conservative solute as a tracer. Tracer tests consisted of both well‐tracer tests (single and double wells) and an aquifer tracer test using a plume‐capturing device such as time domain reflectometry (TDR). The results showed that longitudinal dispersivities estimated from the single and double well‐tracer tests were 2·2 cm and 13·5 cm for a travel distance of 9·3 cm and 13·5 cm from the injection point respectively. These results agreed reasonably well with the results of the aquifer tracer test. The solute transport parameters obtained at multiple points in the aquifer through the aquifer tracer test revealed that the dispersivity length was proportional to the travel distance by a factor of 0·3, which was moderately higher than the value of 0·1 given in the literature. Copyright © 2002 John Wiley & Sons, Ltd.     

Isotopes and geochemistry in a managed aquifer recharge scheme: a case study of fresh water injection at the Damascus University Campus,Syria

A combination of stable isotopes (¹⁸O and ²H) and hydrochemistry has been applied to investigate storage processes in relation to aquifer storage and recovery (ASR) of the shallow alluvial Quaternary aquifer in Damascus basin. The stored water, entirely taken from the Figeh springs during flood periods, was injected in a single well having a brackish groundwater. Water samples were collected from four observation wells drilled in the Damascus University Campus (DUC) site during a 3‐year period (2006–2008). The injectant water, which deviates in its chemical and isotopic signatures from that of the ambient groundwater, shows that the stored water plume remains within close proximity to the injection well (IW) (<≈ 100 m). Thus, only two wells (W13 and W14) located at a distance less than 80 m from the injection point were affected by this injection. The observation wells located at longer distances from the IW (≈145 m and ≈ 600 m for wells W15 and WHz, respectively) were completely unaffected by the injection. Although most of the chemical and isotopic parameters usefully reflected the mixing process that occurs between the injectant water and ambient groundwater, the stable isotope (¹⁸O) and chloride (Cl⁻) were the most sensitive parameters that quickly reflect this signature. Using a simple mass balance, the calculated proportion of injectant water reaching the well W13 was in the range of 50–90%. This proportion was even lower (30–55%) in the case of well W14. Although the drought event prevailing during this study did not much help to inject further amounts of water, higher than the injected volume (0·2416 M m³) and also not favourable to better evaluate the fate and subsurface hydrological processes, these findings offer encouragement to continue the ASR activities, as an alternative way for better management of water resources in this basin facing intensive problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Evaluation of an over-used coastal aquifer (Orissa,India) using statistical approaches

Abstract

This study evaluates an over-exploited aquifer (Balasor, India) and also explores the possibilities of sustainable management using several statistical approaches. First, bootstrap analysis indicates that groundwater pumping has resulted in the reduction of mean cultivated area as the average irrigation capacity per bore well dropped from 3.74 ha to 1.5 ha within a period of 10 years of operation. However, modelling the groundwater levels using the seasonal autoregressive integrated moving average (SARIMA) procedure showed no evidence of large-scale groundwater withdrawals. The derived models can be used for water table forecasting and also for infilling the missing observations. The empirical relationship among pumping test results indicated that well depth and aquifer thickness significantly influence the discharge from the aquifer. This discharge may have encouraged the resource-rich farmers to exploit the lower aquifer. Based on a contour plot, the zone of groundwater exploitation was estimated to vary from 25 to 60 m below the surface. Therefore, a uniform aquifer exploitation policy needs to be implemented to curb the vertical competition in exploiting the aquifer and to develop sustainable management of the groundwater.

Citation Panda, D. K. & Kumar, A. (2011) Evaluation of an over-used coastal aquifer (Orissa, India) using statistical approaches. Hydrol. Sci. J. 56(3), 486–497.     

Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal     

Isotopic tracing for conceptual models of groundwater hydrodynamics in multilayer aquifer systems of central and southern Tunisia

Abstract

New and previously published data sets including stable and radiogenic isotope measurements (¹⁸O, ²H, ³H, ¹³C and ¹⁴C) were used to investigate, conceptualize and compare groundwater hydrodynamics within three major multilayer aquifer systems located in central and southern Tunisia. It has been demonstrated that the investigated aquifer systems contain modern and palaeoclimatic waters. Modern groundwaters, which refer to contemporaneous and post-nuclear recharge waters, are characterized by enriched stable isotope contents, high carbon-14 activities and high to moderate tritium concentration. While palaeoclimatic groundwaters, which refer to Late Pleistocene and Early Holocene recharge waters, are distinguished by their depleted stable isotope contents, low carbon-14 activities and insignificant tritium concentrations. Established conceptual models have elucidated the groundwater hydrodynamics within the studied aquifer systems. They show that groundwater mixing occurs between end-members from the shallow and deep aquifers that migrate by downward and upward leakage towards the intermediate aquifer.

Editor D. Koutsoyiannis; Associate editor S. Faye

Citation Dassi, L. and Tarki, M., 2014. Isotopic tracing for conceptual models of groundwater hydrodynamics in multilayer aquifer systems of central and southern Tunisia. Hydrological Sciences Journal, 59 (6), 1240–1258. http://dx.doi.org/10.1080/02626667.2014.892206     

A stochastic optimization framework for the restoration of an over-exploited aquifer

ABSTRACT

This study investigates the impact of hydraulic conductivity uncertainty on the sustainable management of the aquifer of Lake Karla, Greece, using the stochastic optimization approach. The lack of surface water resources in combination with the sharp increase in irrigation needs in the basin over the last 30 years have led to an unprecedented degradation of the aquifer. In addition, the lack of data regarding hydraulic conductivity in a heterogeneous aquifer leads to hydrogeologic uncertainty. This uncertainty has to be taken into consideration when developing the optimization procedure in order to achieve the aquifer’s sustainable management. Multiple Monte Carlo realizations of this spatially-distributed parameter are generated and groundwater flow is simulated for each one of them. The main goal of the sustainable management of the ‘depleted’ aquifer of Lake Karla is two-fold: to determine the optimum volume of renewable groundwater that can be extracted, while, at the same time, restoring its water table to a historic high level. A stochastic optimization problem is therefore formulated, based on the application of the optimization method for each of the aquifer’s multiple stochastic realizations in a future period. In order to carry out this stochastic optimization procedure, a modelling system consisting of a series of interlinked models was developed. The results show that the proposed stochastic optimization framework can be a very useful tool for estimating the impact of hydraulic conductivity uncertainty on the management strategies of a depleted aquifer restoration. They also prove that the optimization process is affected more by hydraulic conductivity uncertainty than the simulation process.

Editor Z.W. Kundzewicz; Guest editor S. Weijs     

Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases,southern Tunisia: implications of return flow

Abstract

Major ions and stable isotopes in groundwaters of the Plio-Quaternary shallow aquifer of the Djerid oases, southern Tunisia, were investigated to elucidate the origin of groundwater recharge and the mineralization processes. It has been demonstrated that the groundwater composition is mainly controlled by the water–rock interaction, the encroachment of brines from the Chotts and the return flow of irrigation waters. The isotopically depleted groundwater samples suggest that the recharge waters derive from an old palaeoclimatic origin. However, the enriched groundwater samples reflect the presence of evaporated recharge water. Furthermore, the large negative deuterium-excess values indicate the effect of secondary evaporation processes, probably related to the return flow of irrigation waters pumped from the underlying aquifer.

Editor D. Koutsoyiannis; Associate editor E. Custodio

Citation Tarki, M., Dassi, L. and Jedoui, Y., 2012. Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases, southern Tunisia: implications of return flow. Hydrological Sciences Journal, 57 (4), 790–804.     

Modelling the water level of the alluvial aquifer of an ephemeral river in south-western Zimbabwe

ABSTRACT

Water from the alluvium of ephemeral rivers in Zimbabwe is increasingly being used. These alluvial aquifers are recharged annually from infiltrating floodwater. Nonetheless, the size of this water resource is not without limit and an understanding of the hydrological processes of an alluvial aquifer is required for its sustainable management. This paper presents the development of a water balance model, which estimates the water level in an alluvial aquifer recharged by surface flow and rainfall, while allowing for abstraction, evaporation and other losses. The model is coupled with a watershed model, which generates inflows from upland catchment areas and tributaries. Climate, hydrological, land cover and geomorphological data were collected as inputs to both models as well as observed flow and water levels for model calibration and validation. The sand river model was found to be good at simulating the observed water level and was most sensitive to porosity and seepage.     

Prediction of groundwater‐induced flooding in a chalk aquifer for future climate change scenarios

Current climate change models for the southeast UK predict changing rainfall patterns, with increased incidence of extreme events. The chalk aquifer in the UK and northern France is susceptible to groundwater‐induced flooding under such conditions. In this methodological study we apply a frequency domain analysis approach to the chalk aquifer to derive a transfer function between effective rainfall and groundwater level from 7 years of monitoring data from the North Heath Barn site, near Brighton. The derived transfer function was calibrated and validated against monitoring data and then used to predict groundwater level for rainfall models for high, medium and low emission scenarios from the UKCP09 database. The derived transfer function is most closely comparable to the linear aquifer model, despite evidence for both matrix and fracture or karst water flow in the chalk, with transmissivity and unconfined storativity at the catchment scale of 1548 m² day^?1 and 1.6 × 10^?2. The application of the transfer function to UKCP09 rainfall data suggests that groundwater‐induced flooding may be about four times more frequent by 2040–2069 compared with 1961–1990 and seven times more frequent by 2070–2099. The model data also suggest an increase in the duration of groundwater minima relative to the reference period. Compared to deterministic modelling which requires detailed knowledge of aquifer heterogeneity and processes, the transfer function approach, although with limitations, is simpler, incorporating these factors into the analysis through frequency and phase coefficients, and thus may have the potential for groundwater risk assessment in other areas. Copyright © 2015 John Wiley & Sons, Ltd.     

A double-porosity model for a fractured aquifer with non-Darcian flow in fractures

Abstract

Non-Darcian flow in a finite fractured confined aquifer is studied. A stream bounds the aquifer at one side and an impervious stratum at the other. The aquifer consists of fractures capable of transmitting water rapidly, and porous blocks which mainly store water. Unsteady flow in the aquifer due to a sudden rise in the stream level is analysed by the double-porosity conceptual model. Governing equations for the flow in fractures and blocks are developed using the continuity equation. The fluid velocity in fractures is often too high for the linear Darcian flow so that the governing equation for fracture flow is modified by Forcheimer's equation, which incorporates a nonlinear term. Governing equations are coupled by an interaction term that controls the quasi-steady-state fracture—block interflow. Governing equations are solved numerically by the Crank-Nicolson implicit scheme. The numerical results are compared to the analytical results for the same problem which assumes Darcian flow in both fractures and blocks. Numerical and analytical solutions give the same results when the Reynolds number is less than 0.1. The effect of nonlinearity on the flow appears when the Reynolds number is greater than 0.1. The higher the rate of flow from the stream to the aquifer, the higher the degree of nonlinearity. The effect of aquifer parameters on the flow is also investigated. The proposed model and its numerical solution provide a useful application of nonlinear flow models to fractured aquifers. It is possible to extend the model to different types of aquifer, as well as boundary conditions at the stream side. Time-dependent flow rates in the analysis of recession hydrographs could also be evaluated by this model.     

Groundwater evolution and mean water age inferred from hydrochemical and isotopic tracers in a tropical confined aquifer

The coastal confined aquifer in the Gulf of Urabá (Colombia) is an important water source for the banana agro‐industry as well as for urban and rural communities. However, the main processes controlling recharge and mixing in the aquifer are still poorly understood. Hydrochemical analyses and stable isotope monitoring were conducted to (a) determine groundwater recharge origin, mean groundwater age, and the main processes governing groundwater chemistry and the potential mixing of marine water and the influence of diffusive processes from the two surrounding aquitard layers. Hydrochemical data indicate that the main processes affecting the dissolved chemical composition include cation exchange, dissolution of carbonated and CO₂, and silicate weathering. δ¹⁸O and δ²H compositions combined with ¹⁴C data highlight the differences in climatic conditions between the recharge zone and the confined section of the aquifer, which is close to the Atlantic Ocean. Groundwater samples with ¹⁴C ages from recent to 28,300 years BP show a depleted isotopic trend ranging from ?6.43‰ to ?9.14‰ in δ¹⁸O and from ?43.2‰ to ?65.7‰ in δ²H. The most depleted δ¹⁸O and δ²H compositions suggest a cooler recharge climate than the current conditions (corresponding to the last glacial period of the late Pleistocene). Depleted δ¹³C values in the total dissolved inorganic carbon indicate the existence of organic material oxidation processes within the geologic formation. These results can be used or transferred to enhance groundwater modelling efforts in other confined coastal aquifers of South America where scarcity of long‐term monitoring data limits water resources planification under a changing climate.     

Natural 222Rn as a tracer of mixing and volatilization in a shallow aquifer during a CO2 injection experiment

This study aims to evaluate the application of ²²²Rn in groundwater as a tracer for monitoring CO₂ plume migration in a shallow groundwater system, which is important to detect potential CO₂ leakage in the carbon capture and storage (CCS) project. For this research, an artificial CO₂-infused water injection experiment was performed in a shallow aquifer by monitoring hydrogeochemical parameters, including ²²²Rn. Radon in groundwater can be a useful tracer because of its sensitivity to sudden changes in subsurface environment. To monitor the CO₂ plume migration, the data were analysed based on (a) the influence of mixing processes on the distribution of ²²²Rn induced by the artificial injection experiment and (b) the influence of a carrier gas role by CO₂ on the variation of ²²²Rn. The spatio-temporal distributions of radon concentrations were successfully explained in association with horizontal and vertical mixing processes by the CO₂-infused water injection. Additionally, the mixing ratios of each monitoring well were calculated, quantitatively confirming the influence of these mixing processes on the distribution of radon concentrations. Moreover, one monitoring well showed a high positive relationship between ²²²Rn and Total dissolved inorganic carbon (TIC) by the carrier gas effect of CO₂ through volatilization from the CO₂ plume. It indicated the applicability of ²²²Rn as a sensitive tracer to directly monitor CO₂ leakage. When with a little effect of carrier gas, natural ²²²Rn in groundwater can be used to compute mixing ratio of CO₂-infused water indicative of CO₂ migration pathways. CO₂ carrier gas effect can possibly increase ²²²Rn concentration in groundwater and, if fully verified with more field tests, will pose a great potential to be used as a natural tracer for CO₂.     

Basin structure and water quality of the Oligocene aquifer in the Sebkhet El Behira basin (central Tunisia)

ABSTRACT

The Sebkhet El Behira–Garaat El Majdoul multilayer aquifer system is composed of Mio-Plio-Quaternary and Oligocene groundwater exploited in Sidi Bouzid and Kairouan governorates. Annual withdrawal volumes from El Behira Oligocene groundwater were about 0.37 hm³ in 2005 and 0.36 hm³ in 2008. The present study of the Oligocene reservoir in the Sebkhet El Behira–Garaat El Majdoul basin is based on various data including 2D seismic sections, petroleum wells, field geological cross-sections and geochemical analysis. The gathered data allowed us to establish the tectonic framework and to define its influence on the structure of the aquifers seated in the deep Oligocene sandstone reservoirs. Three Oligocene sub-basin aquifers are defined showing different depth, thickness and petrophysical characteristics: lower salinity (<1.5 g/L) and higher porosity (30%). Two field sections in the northern and in the southern parts of Cherahil anticline exhibit the presence of four to five sandstone levels. The anticlines limit the Sebkhet El Behira–Garaat El Majdoul basin. In addition, the morphostructural configuration controls the piezometry of underground flows in the Oligocene and Mio-Plio-Quaternary unconfined aquifers.

Editor D. Koutsoyiannis; Associate editor M. Besbes     

Analysis of benzene transport in a two‐dimensional aquifer model

In this study, the fate and transport of aqueous benzene was investigated in a laboratory‐scale homogeneous aquifer by conducting a two‐dimensional plume test. Benzene solution was introduced as a pulse type along the width of the aquifer model through a recharge zone situated at the upper‐left part of the model and followed by a steady state flow. Solution samples were collected at various locations on the front side of the model to capture two‐dimensional plumes at discrete time intervals. The benzene plumes showed a moderate retardation relative to chloride plumes observed from the previous study conducted for the same aquifer model. The retardation factor was obtained from the ratio of travel distances of benzene peaks to chloride peaks from the injection point, computed using a line integral method. Mass recovery of aqueous benzene revealed that there was a significant reduction of benzene mass, indicating the occurrence of volatilization and/or irreversible sorption during transport. Thus, retardation along with volatilization and/or irreversible sorption may be important processes affecting the fate and transport of aqueous benzene in the aquifer model. Copyright © 2005 John Wiley & Sons, Ltd.     

A combined methodology for estimating the potential natural aquifer recharge in an arid environment

ABSTRACT

An innovative methodology that combines an indirect physiography-based method for determining the runoff coefficient at a sub-basin scale and a water balance model applied on a daily time scale was developed to calculate the natural groundwater recharge in three watersheds within the Oum Zessar arid area, Tunisia. The effective infiltration was calculated as part of the water surplus by considering the average available water content (AWC) of soil and an average runoff coefficient for each sub-basin. The model indicates that the sub-basins covered mainly by the “artificial” soils of tabias and jessour, characterized by average AWC values greater than 150 mm, did not contribute to natural groundwater recharge over the 10-year period (2003–2012) considered. The estimated volume for the Triassic aquifer amounted to about 4.5 hm³ year^?1, which is consistent with previous studies. For the Jurassic and Cretaceous aquifers, the estimated volumes amounted to about 200 dm³ year^?1.     

A new method for aquifer system identification and parameter estimation

The standard practice for assessing aquifer parameters is to match groundwater drawdown data obtained during pumping tests against theoretical well function curves specific to the aquifer system being tested. The shape of the curve derived from the logarithmic time derivative of the drawdown data is also very frequently used as a diagnostic tool to identify the aquifer system in which the pumping test is being conducted. The present study investigates the incremental area method (IAM) to serve as an alternative diagnostic tool for the aquifer system identification as well as a supplement to the aquifer parameter estimation procedure. The IAM based diagnostic curves for ideal confined, leaky, bounded and unconfined aquifers have been derived as part of this study, and individual features of the plots have been identified. These features were noted to be unique to each aquifer setting, which could be used for rapid evaluation of the aquifer system. The effectiveness of the IAM methodology was investigated by analyzing field data for various aquifer settings including leaky, unconfined, bounded and heterogeneous conditions. The results showed that the proposed approach is a viable method for use as a diagnostic tool to identify the aquifer system characteristics as well as to support the estimation of the hydraulic parameters obtained from standard curve matching procedures. Copyright © 2012 John Wiley & Sons, Ltd.     

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.