A multi-method approach for groundwater resource assessment in coastal carbonate (karst) aquifers: the case study of Sierra Almijara (southern Spain)

Understanding the transference of water resources within hydrogeological systems, particularly in coastal aquifers, in which groundwater discharge may occur through multiple pathways (through springs, into rivers and streams, towards the sea, etc.), is crucial for sustainable groundwater use. This research aims to demonstrate the usefulness of the application of conventional recharge assessment methods coupled to isotopic techniques for accurately quantifying the hydrogeological balance and submarine groundwater discharge (SGD) from coastal carbonate aquifers. Sierra Almijara (Southern Spain), a carbonate aquifer formed of Triassic marbles, is considered as representative of Mediterranean coastal karst formations. The use of a multi-method approach has permitted the computation of a wide range of groundwater infiltration rates (17–60%) by means of direct application of hydrometeorological methods (Thornthwaite and Kessler) and spatially distributed information (modified APLIS method). A spatially weighted recharge rate of 42% results from the most coherent information on physiographic and hydrogeological characteristics of the studied system. Natural aquifer discharge and groundwater abstraction have been volumetrically quantified, based on flow and water-level data, while the relevance of SGD was estimated from the spatial analysis of salinity, ²²²Rn and the short-lived radium isotope ²²⁴Ra in coastal seawater. The total mean aquifer discharge (44.9–45.9 hm³ year^?1) is in agreement with the average recharged groundwater (44.7 hm³ year^?1), given that the system is volumetrically equilibrated during the study period. Besides the groundwater resources assessment, the methodological aspects of this research may be interesting for groundwater management and protection strategies in coastal areas, particularly karst environments.     

Groundwater exploitation and recharge rate estimation of a quaternary sand aquifer in Dar-es-Salaam area,Tanzania

Dar-es-Salaam City gets water supply from surface water and groundwater. The groundwater is used to supplement surface water supply and has increasingly become a major source of water supply in the city. The study area comprises three major parts: the central coastal plain with quaternary fluvial–deltaic sediments, the deltaic Mio-Pliocene clay-bound sands and gravels in the northwest and southeast and the Lower Miocene fluviatile sandstones of Pugu Hills in the west of the study area. The main objective of this study was to quantify the integrated water balance. The major source of renewable groundwater in the aquifer is rainfall. Hence, the average recharge of 256.2 mm/year (for the year 2006) to the aquifer was estimated using the balance method of Thornthwaite and Mather, which is equal to 99.4 hm³/year for the whole alluvial aquifer. This value was balanced with total groundwater abstraction of 8.59 hm³/year, baseflow to rivers of 75.7 hm³/year and discharge into the sea (15.11 hm³/year).     

Saltwater intrusion and nitrate pollution in the coastal aquifer of Dar es Salaam, Tanzania

Dar es Salaam Quaternary coastal aquifer is a major source of water supply in Dar es Salaam City used for domestic, agricultural, and industrial uses. However, groundwater overdraft and contamination are the major problems affecting the aquifer system. This study aims to define the principal hydrogeochemical processes controlling groundwater quality in the coastal strip of Dar es Salaam and to investigate whether the threats of seawater intrusion and pollution are influencing groundwater quality. Major cations and anions analysed in 134 groundwater samples reveal that groundwater is mainly affected by four factors: dissolution of calcite and dolomite, weathering of silicate minerals, seawater intrusion due to aquifer overexploitation, and nitrate pollution mainly caused by the use of pit latrines and septic tanks. High enrichment of Na⁺ and Cl^? near the coast gives an indication of seawater intrusion into the aquifer as also supported from the Na–Cl signature on the Piper diagram. The boreholes close to the coast have much higher Na/Cl molar ratios than the boreholes located further inland. The dissolution of calcite and dolomite in recharge areas results in Ca–HCO₃ and Ca–Mg–HCO₃ groundwater types. Further along flow paths, Ca²⁺ and Na⁺ ion exchange causes groundwater evolution to Na–HCO₃ type. From the PHREEQC simulation model, it appears that groundwater is undersaturated to slightly oversaturated with respect to the calcite and dolomite minerals. The results of this study provide important information required for the protection of the aquifer system.     

Effect of diffuse recharge and wastewater on groundwater contamination in Douala,Cameroon

High water demand for domestic use in Douala with over 3 million inhabitants is met mainly by shallow groundwater. Field measurements and water sampling in January 2015 were carried out to examine the major controls on the groundwater composition and spatial view of ions in the water, timing of recharge and link between the recharge process and quality of the water. Fifty-two water samples were analysed for major ions and stable hydrogen and oxygen isotopes. Low pH values (3.61–6.92) in the groundwater indicated an acidic aquifer; thus, prone to acidification. The dominant water type was Na–Cl. Nitrate, which exceeded the WHO guide value of 50 mg/l in 22% of the groundwater, poses a health problem. Mass ratios of Cl^?/Br^? in the water ranged from 54 to 3249 and scattered mostly along the mixing lines between dilute waters, septic-tank effluent and domestic sewage. A majority of the samples, especially the high NO₃ ^? shallow wells, clustered around the septic-tank effluent-end-member indicating high contamination by seepage from pit latrines; hence, vulnerable to pollution. Stable isotopes in the groundwater indicated its meteoric origin and rapid infiltration after rainfall. The δ¹⁸O values showed narrow ranges and overlaps in rivers, springs, open wells and boreholes. These observations depict hydraulic connectivity, good water mixing and a homogeneous aquifer system mainly receiving local direct uniform areal recharge from rainfall. The rapid and diffused recharge favours the leaching of effluent from the pit toilets into the aquifer; hence, the high NO₃ ^? and Cl^? in shallow wells. Silicate weathering, ion exchange and leaching of waste from pit toilets are the dominant controls on the groundwater chemistry. Drilling of deep boreholes is highly recommended for good-quality water supply. However, due the hydraulic connection to the shallow aquifer, geochemical modelling of future effects of such an exploitation of the deeper aquifer should support groundwater management and be ahead of the field actions.     

Hydrochemistry of groundwater in the Savelugu–Nanton District, Northern Ghana

Calcite (CaCO₃), dolomite [CaMg(CO₃)₂], silicate dissolution, ion-exchange and reverse ion-exchange reactions are the predominant processes influencing groundwater quality in the Savelugu?CNanton District. The main objective of this study is to characterize groundwater and delineate water?Crock interactions responsible for the chemical evolution of groundwater in the District. Eighty-one (81) boreholes were sampled for quality assessment. Results showed that, the pH of the boreholes are slightly acidic to basic ranging from 6.1?C8.3?pH units. Conductivity values are low to high ranging from 147?C23,200???S/cm with, 23.5?% of groundwater within the study area being either brackish (1,500?C5,000???S/cm) or saline (>5,000???S/cm) and therefore, unsuitable for potable purposes. TDS values in groundwater varied widely, ranging from 62?C11,900?mg/L. 61.7?% of groundwater within the study area are fresh (TDS?<?500?mg/L). The chemical constituents generally, have low concentrations and are within the WHO (Guidelines for drinking water quality. Revision of the 1993 Guidelines. Final Task Group Meeting. Geneva, 2004) Guideline values. The relative abundance of cations and anions in the groundwater are in the order: Ca²⁺?>?Na⁺?>?Mg²⁺?>?K⁺ and HCO₃ ^??>?Cl^??>?SO₄ ^2??>?NO₃ ^? respectively. Multivariate statistical analysis showed expected process-based relationship derived mainly from the geochemical and biochemical processes within the aquifer. Hydrochemical facies using piper plot of major ions showed one major hydrochemical water type. The Ca?CMg?CHCO₃ water type. Due to the high cost of drilling of boreholes coupled with the high percentage (×1?%) of people without access to potable water in the Northern Region, it is recommended that the Government of Ghana and other stakeholders within the Water Sector take immediate measures to reduce (to recommended limits for potable uses) the levels of dissolved solids either by installing Reverse Osmosis equipments on such boreholes or employ other relatively known cheaper methods to reduce the dissolved solids to recommended limits. High yielding boreholes with hydrochemical parameters within WHO guideline limits in the Savelugu?CNanton District could also be mechanized to serve a wider area.     

Large-scale hydraulic conductivity distribution in an unconfined carbonate aquifer using the ocean tidal propagation

The hydraulic conductivity of an unconfined carbonate aquifer at the uplifted atoll of Minami-Daito, Japan, was evaluated by a combination of cross-spectral analysis, analytical solution, and density-dependent groundwater modeling based on observed groundwater levels in 15 wells and at sea level. The island area was divided into 10 subregions based on island morphology and on inland propagation of ocean tides. The hydraulic conductivity was obtained for each subregion using analytical solutions based on phase lags of M₂ constituents of ocean tides at each well by assuming two aquifer thicknesses (300 and 1,800 m) and two effective porosities (0.1 and 0.3). The density-dependent groundwater model evaluated the hydraulic conductivity of the subregions by reproducing observed groundwater levels. The hydraulic conductivity in the subregions was estimated as 3.46?×?10^?3 to 6.35?×?10^?2 m/s for aquifer thickness of 300 m and effective porosity of 0.1, and as 1.73?×?10^?3 to 3.17?×?10^?2 m/s for aquifer thickness of 1,800 m and the effective porosity of 0.3. It was higher in southern and northern areas, and higher in interior lowland than in the western and eastern areas. Fissures and dolomite distributions on the island control differences of the omnidirectional ocean tidal propagation and cause these differences in hydraulic conductivity. The method used for this study may also be applicable to other small islands that have few or no data for hydraulic conductivity.

     

Groundwater flow modeling of Quaternary aquifer Ras Sudr,Egypt

Recently, Ras Sudr (the delta of Wadi Sudr) area received a great amount of attention due to different development expansion activities directed towards this area. Although Quaternary aquifer is the most prospective aquifer in Ras Sudr area, it has not yet completely evaluated. The present work deals with the simulation of the Quaternary groundwater system using a three-dimensional groundwater flow model. MODFLOW code was applied for designing the model of the Ras Sudr area. This is to recognize the groundwater potential as well as exploitation plan of the most prospective aquifer in the area. The objectives were to determine the hydraulic parameters of the Quaternary aquifer, to estimate the recharge amount to the aquifer, and to determine the hydrochemistry of groundwater in the aquifer. During this work, available data has been collected and some field investigation has been carried out. Groundwater flow model has been simulated using pilot points’ method. SEAWAT has been also applied to simulate the variable-density flow and sea water intrusion from the west. It can be concluded that: (1) the direction of groundwater flow is from the east to the west, (2) the aquifer system attains a small range of log-transformed hydraulic conductivity. It ranges between 3.05 and 3.35 m/day, (3) groundwater would be exploited by about 6.4 × 10⁶ m³/year, (4) the estimated recharge accounts for 3 × 10⁶ m³/year, (5) an estimated subsurface flow from the east accounts for 2.7 × 10⁶ m³/year, (6) the increase of total dissolved solids (TDS) most likely due to dilution of salts along the movement way of groundwater from recharge area to discharge area in addition to a contribution of sea water intrusion from the west. Moreover, it is worth to note that a part of TDS increase might be through up coning from underlying more saline Miocene sediments. It is recommended that: (1) any plan for increasing groundwater abstraction is unaffordable, (2) reliable estimates of groundwater abstraction should be done and (3) automatic well control system should be made.     

Using tritium to document the mean transit time and sources of water contributing to a chain-of-ponds river system: Implications for resource protection

Documenting the interaction between groundwater and rivers is fundamental to understanding hydrological systems. While many studies have examined the location and magnitude of groundwater inflows to rivers, much less is known about the transit times of water in catchments and from where in the aquifer the groundwater originates. Resolving those questions is vital for protecting riverine ecosystems, assessing the impact of contamination, and understanding the potential consequences of groundwater pumping. This study uses tritium (³H) to evaluate the mean transit times of water contributing to Deep Creek (southeast Australia), which is a chain-of-ponds river system. ³H activities of river water vary between 1.47 and 2.91 TU with lower ³H activities recorded during cease-to-flow periods when the river comprises isolated groundwater-fed pools. Regional groundwater 1–2.5 km away from Deep Creek at depths of 7.5–46.5 m has ³H activities of between <0.02 and 0.84 TU. The variation in ³H activities suggest that the water that inflows into Deep Creek is dominated by near-river shallow groundwater with the deeper groundwater only providing significant inflows during drier periods. If the water in the catchment can be represented by a single store with a continuum of ages, mean transit times of the river water range between <1 and 31 years whereas those of the groundwater are at least 75 years and mainly >100 years. Alternatively the variation in ³H activities can be explained by mixing of a young near-river water component with up to 50% older groundwater. The results of this study reinforce the need to protect shallow near-river groundwater from contamination in order to safeguard riverine ecosystems and also illustrate the potential pitfalls in using regional bores to characterise the geochemistry of near-river groundwater.     

Characterization of the subsurface urban heat island and its sources in the Milan city area,Italy

Urban areas are major contributors to the alteration of the local atmospheric and groundwater environment. The impact of such changes on the groundwater thermal regime is documented worldwide by elevated groundwater temperature in city centers with respect to the surrounding rural areas. This study investigates the subsurface urban heat island (SUHI) in the aquifers beneath the Milan city area in northern Italy, and assesses the natural and anthropogenic controls on groundwater temperatures within the urban area by analyzing groundwater head and temperature records acquired in the 2016–2020 period. This analysis demonstrates the occurrence of a SUHI with up to 3 °C intensity and reveals a correlation between the density of building/subsurface infrastructures and the mean annual groundwater temperature. Vertical heat fluxes to the aquifer are strongly related to the depth of the groundwater and the density of surface structures and infrastructures. The heat accumulation in the subsurface is reflected by a constant groundwater warming trend between +0.1 and?+?0.4 °C/year that leads to a gain of 25 MJ/m² of thermal energy per year in the shallow aquifer inside the SUHI area. Future monitoring of groundwater temperatures, combined with numerical modeling of coupled groundwater flow and heat transport, will be essential to reveal what this trend is controlled by and to make predictions on the lateral and vertical extent of the groundwater SUHI in the study area.

     

Assessing seasonal variations and aquifer vulnerability in coastal aquifers of semi-arid regions using a multi-tracer isotopic approach: the case of Grombalia (Tunisia)

The Grombalia aquifer (NE Tunisia) is an example of an important source of water supply for regional and national development, where the weak controls over abstraction, fertilizer application and waste disposal, coupled with limited knowledge of aquifer dynamics, is causing aquifer over-exploitation and water quality degradation. Assessing the key role of groundwater in water-resources security is therefore of paramount importance to support new actions to preserve water quality and quantity in the long-run. This study presents one of the first investigations targeted at a complete assessment of aquifer dynamics in the Grombalia aquifer. A multi-tracer hydrogeochemical and isotopic (δ²H, δ¹⁸O and ³H) approach was used to study the influence of seasonal variation on piezometric levels, chemical and isotopic compositions, and groundwater recharge. A total of 116 samples were collected from private wells and boreholes during three periods in a 1 year monitoring campaign (February–March 2014, September 2014 and February 2015). Results revealed the overall unsuitability of groundwater for drinking and irrigation purposes (NO₃?>?50 mg/L in 51% of the wells; EC >1,000 μS/cm in 99% of the wells). Isotopic balance coupled to piezometric investigation indicated the contribution of the shallow aquifer to deep groundwater recharge. The study also revealed the weakness of ‘business as usual’ management practices, highlighting possible solutions to tackle water-related challenges in the Grombalia region, where climate change, population growth and intensive agricultural activities have generated a large gap between demand and available water reserves, hence becoming a possible driver for social insecurity.     

Naturally occurring arsenic: Mobilization at a landfill in Maine and implications for remediation

Elevated levels of dissolved arsenic (∼300 μg L⁻¹) have been detected beneath and in groundwater plumes extending away from a closed landfill in southern Maine. This study sought to determine the source of arsenic to the aquifer, the processes responsible for arsenic mobilization, and to evaluate the effectiveness of remediation efforts that have occurred at this site. The As appears to originate in the natural (glacial) aquifer solids, which contain ∼5 mg kg⁻¹ As on a dry weight basis. This conclusion is supported by the relatively uniform distribution of As in sediment samples, results of laboratory batch incubation experiments, and comparisons with groundwaters in nearby wetlands, which also have high levels of dissolved As that do not appear to originate within the landfill. The As is mobilized in the subsurface by strongly reducing conditions beneath the landfill and in nearby wetlands. In the aquifer beneath the landfill, the average oxidation–reduction potential (ORP) is −95 mV (Eh + 105 mV), and these reducing conditions were primarily induced by landfill leachate. Remediation efforts at this site have included installation of a low permeability clay cap; groundwater extraction, oxidation, and re-injection; and subsurface oxidation by injection of magnesium peroxide. The natural source of arsenic within the aquifer solids, coupled with widespread reducing conditions, has severely limited the effectiveness of these interventions on groundwater arsenic concentrations.     

Quantification of groundwater–seawater interaction in a coastal sandy aquifer system: a study from Panama,Sri Lanka

The Panama coastal aquifer system is an important water resource in the southeast coast of Sri Lanka that provides adequate supplies of water for agriculture and domestic uses. One of the biggest threats to these fragile aquifers is the sea water intrusion. In this study, recharging mechanism and geochemical evaluation of groundwater in the coastal sandy aquifer of Panama were evaluated using chemical and stable isotope techniques. Thirty groundwater samples were collected and analyzed for their major ion concentrations and stable isotope ratios of oxygen (¹⁸O/¹⁶O) and hydrogen (D/H). All studied samples showed a ranking of major anions in the order Cl^? > HCO ₃ ^?  > SO ₄ ^2?  > N-NO₃ ^? while cations showed a decreasing order of abundance with Na⁺ > Ca²⁺ > Mg²⁺ > K⁺. Dominant groundwater hydrogeochemical types were Na–Cl and mixed Ca–Mg–Cl. Results of saturation index calculations indicate that the investigated groundwater body was mostly saturated with respect to calcite, dolomite and gypsum. In addition, stable isotope and geochemical data suggest that fresh groundwater in the aquifer is recharged mainly by local precipitation with slight modification from evaporation and saline water intrusions. Isotope data suggest that mixing of salt water with freshwater occurs in aquifers which are located towards the lagoon. Since the communities in the study area depend entirely on groundwater, an understanding of the hydrogeochemical characteristics of the aquifer system is extremely important for the better water resource management in the region.     

Hydraulic response of a tidally forced coastal aquifer, Pontal do Paraná, Brazil

The unconfined aquifer supported by the extensive Holocene coastal plain of southern Brazil is found dissected at several places by tidally forced streams. One such unit of coastal plain, falling within the watershed limits of Perequê Stream in Paraná state, was analyzed using MODFLOW to assess the surface and subsurface water components under different stress conditions. The steady-state simulation suggests that ~23,000 m³/day of water forms the total budget of the designated aquifer zone and this zone discharges ~16,000 m³/day across the interface between the freshwater aquifer zone and the sea as direct submarine discharge. The maximum evaporation loss was assumed to be 2,700 m³/day. Any change in effective aquifer storage was insignificant, as there is no large-scale exploitation of local groundwater reserves. A calibration exercise, based on field measurements lasting 1 year, supports validity of the model considered. Simulations representing two different transient conditions point to the dynamic hydraulic response exhibited by the aquifer. A schematic map, showing directions and velocity of subsurface flows, and a three-dimensional model of the groundwater reserve helped to visualize the hydrogeologic changes and to formulate management plans.     

Behaviour of boron and strontium isotopes in groundwater–aquifer interactions in the Cornia Plain (Tuscany,Italy)

The Cornia Plain alluvial aquifer, in Tuscany, is exploited intensely to meet the demand for domestic, irrigation and industrial water supplies. The B concentration of groundwater, however, is often above the European limit of 1 mg L⁻¹, with the result that exploitation of these water resources requires careful management. Boron and Sr isotopes have been used as part of a study on the origin and distribution of B dissolved in groundwater, and indirectly as a contribution to the development of appropriate water management strategies.The geochemistry of the Cornia Plain groundwater changes from a HCO₃ facies in the inland areas to a Cl facies along the coastal belt, where seawater intrusion takes place. The B concentration of groundwater increases towards the coastal areas, while the ¹¹B/¹⁰B ratio decreases. This indicates that there is an increasing interaction between dissolved B and the sediments forming the aquifer matrix, whose B content is in the order of 100 mg kg⁻¹. Adsorption–desorption exchanges take place between water and the sediment fine fraction rich in clay minerals, with a net release of B from the matrix into the groundwater, and a consequent δ¹¹B shift from positive to negative values. The aquifer matrix sediments therefore seem to be the major source of B dissolved in the groundwater.The groundwater–matrix interactions triggered by the ionic strength increase caused by seawater intrusion can also be detected in the Ca–Na ion exchanges. Dissolved Sr follows a trend similar to that of Ca, while the ⁸⁷Sr/⁸⁶Sr ratio is equal to that of the exchangeable Sr of the aquifer matrix and therefore does not change significantly.These results have helped to define a new strategy for groundwater exploitation, with the final objective of reducing B concentration in the water extracted from the aquifer.     

Groundwater investigation in Kuala Selangor using vertical electrical sounding (VES) surveys

Integrated geoelectrical and hydrochemical surveys were used to investigate and delineate different types of groundwater in the Kuala Selangor alluvial aquifer. Previous hydrogeological borehole investigation showed that this aquifer contains several types of groundwater in relation to its salinity. The high salinity of the groundwater in some areas is believed to be due to either saltwater intrusion from the nearby sea or river infiltration during high tide season. The vertical electrical sounding (VES) method was employed to study and map the subsurface variation of resistivity in the area. For each sounding measurement, a total spread length of 300 m was obtained with a vertical depth penetration of about 60 to 75 m. Chemical analysis of the groundwater samples taken from both shallow and deep boreholes was carried out for the water quality determination. A total of 45 VES stations were succesfully established along three parallel roads with a direction almost perpendicular to the coastal line. The distance between stations varies from 1 to 2 km with a maximum length of about 60 km surveyed line. Results of the vertical electrical soundings as well as the hydrochemistry of the groundwater samples show that the soil and groundwater in the study area can be grouped into fresh and brackish water zones. The subsurface resisitivity sections derived from the VES study suggest that the area is dominated by brackish soil and groundwater zones, especially in the area towards the coast. This result appears to agree well with the groundwater pumped from boreholes scattered around in the area. Water drawn from boreholes near the coast showed higher salinity compared to the water pumped from inland boreholes. Chloride values greater than 250 mg/L are considered to represent the brackish zones whilst values less than 250 mg/L represents zones of fresh soil and groundwater.     

Origin of groundwater salinity and hydrogeochemical processes in a confined coastal aquifer: Case of the Rhône delta (Southern France)

The Rhône delta, South of France (Camargue, 750 km²) is a coastal saline wetland located along the Mediterranean Sea. The confined aquifer of this delta shows high values of electrical conductivity rising from the north (4 mS/cm) to the shoreline (58 mS/cm). This work aims to identify the origin of groundwater salinity and the geochemical processes occurring in this coastal confined aquifer according to the degree of salinity. A natural tracing approach is considered using monthly sampling in 8 piezometers for chemical and isotopic analyses (¹⁸O, ²H, ¹³C_TDIC). Ionic and isotopic ratios demonstrate that strong salinities are due to a simple mixing between Mediterranean seawater and freshwater; seawater contribution reaches up to 98% at 8 km from the shoreline. Seawater intrusion induces a particular groundwater chemistry which varies with the degree of seawater contribution: (1) In the less saline part of the aquifer (seawater contribution <20%), the intrusion induces an increase of Na⁺ in groundwater leading to Ca²⁺/Na⁺ exchange processes. The δ¹³C_TDIC analyses show that matrix exchange processes most likely occur for the less saline samples. (2) In the saline part of the aquifer (seawater contribution >20%), the intrusion induces SO₄ reduction which is confirmed by highly depleted δ¹³C_TDIC values (up to −19‰). The δ¹³C_TDIC also reveals that methanogenesis processes may occur in the most reductive part of the aquifer. Due to SO₄ reduction, the intrusion induces a shift in carbonate equilibrium leading to supersaturation with respect to dolomite and/or magnesian calcite. Thus carbonate precipitation may occur in the area strongly influenced by seawater.     

Modelling surface–groundwater interaction in the Ruataniwha basin, Hawke’s Bay, New Zealand

Modelling groundwater and surface water is important for integrated water resources management, especially when interaction between the river and the aquifer is high. A transient groundwater and surface water flow model was built for Ruataniwha basin, New Zealand. The model covers a long-time period; starting in 1990, when water resources development in the area started, to present date. For a better resolution, the simulation period was divided into 59 stress periods, and each stress period was divided to 10 time steps. The model uses data obtained from surface water, and groundwater collected over the last 20 years. Rivers and streams were divided into 28 segments and flow and streambed data at the beginning and end of each segment was used. Parameter estimation and optimisation ‘PEST’ was used for automatic calibration of hydraulic conductivity, groundwater recharge and storativity; whereas riverbed conductance was manually calibrated. Model results show that the rivers gain from the aquifer considerably more than the river losses. The cumulative groundwater abstraction over the last 20 years is approximately 210 million m³. This amount is very low compared to other water budget components; however, the effect of groundwater abstraction on storage is significant. Based on the results of this study, it was found that the loss of storage over the last 20 years is more than 66 million m³. Results also reveal that the effect of groundwater abstraction on rivers and springs flow is significant. The rivers gain from the groundwater system, and the springs flow have been decreasing.     

The use of vertical electrical sounding resistivity method for the location of low salinity groundwater for irrigation in Chaj and Rachna Doabs

A geoelectrical resistivity survey using vertical electrical sounding (VES) was conducted at Chaj Doab (land between rivers Jhelum and Chenab, Pakistan) and Rachna Doab (land between rivers Chenab and Ravi, Pakistan), with the objective of investigating groundwater conditions. A total of 90 sites were selected with 43 sites in Chaj and 47 sites in Rachna Doabs. The resistivity meter (ABEM Terrameter SAS 4000, Sweden) was used to collect the VES data by employing a Schlumberger electrode configuration, with half current electrode spacings (AB/2) ranging from 2 to 180 m and the potential electrode (MN) from 1 to 40 m. The field data were interpreted using the Interpex IX1D computer software and the resistivity versus depth models for each location was estimated. The outputs of subsurface layers with resistivities and thickness presented in contour maps and 3-D views by using SURFER software were created. A total of 102 groundwater samples from nearby hydrowells at different depths were collected to develop a correlation between the aquifer resistivity of VES and the electrical conductivity (EC) of the groundwater and to confirm the resulted geophysical resistivity models. From the correlation developed, it was observed that the groundwater salinity in the aquifer may be considered low and so safe for irrigation if resistivity >45 Ω m, and marginally fit for irrigation having resistivity between 25 and 45 Ω m. The study area has resistivities from 3.9 to 2,222 Ω m at the top of the unsaturated layer, between 1.21 and 171 Ω m, in the shallow aquifers, and 0.14–152 Ω m in the deep aquifers of the study area. The results indicate that the quality of groundwater is better near the rivers and in the shallow layers compared to the deep layers.     

Salinity mapping of coastal groundwater aquifers using hydrogeochemical and geophysical methods: a case study from north Kelantan,Malaysia

Integrated hydrogeochemical and geophysical methods were used to study the salinity of groundwater aquifers along the coastal area of north Kelantan. For the hydrogeochemical investigation, analysis of major ion contents of the groundwater was conducted, and other chemical parameters such as pH and total dissolved solids were also determined. For the geophysical study, both geoelectrical resistivity soundings and reflection seismic surveys were conducted to determine the characteristics of the subsurface and groundwater contained within the aquifers. The pH values range from 6.2 to 6.8, indicating that the groundwater in the study area is slightly acidic. Low content of chloride suggests that the groundwater in the first aquifer is fresh, with an average concentration of about 15.8 mg/l and high geoelectrical resistivity (>45 ohm m). On the other hand, the groundwater in the second aquifer is brackish, with chloride concentration ranging from 500 mg/l to 3,600 mg/l and very low geoelectrical resistivity (<45 ohm m) as well as high concentration of total dissolved solids (>1,000 mg/l). The groundwater in the third aquifer is fresh, with chloride concentrations generally ranging from 2 mg/l to 210 mg/l and geoelectrical resistivity of greater than 45 ohm m. Fresh and saltwater interface in the first aquifer is generally located directly in the area of the coast, but, for the second aquifer, both hydrogeochemical and geoelectrical resistivity results indicate that the fresh water and saltwater interface is located as far as 6 km from the beach. The considerable chloride ion content initially suggests that the salinity of the groundwater in the second aquifer is probably caused by the intrusion of seawater. However, continuous monitoring of the chloride content of the second aquifer indicated no significant changes with time, from which it can be inferred that the salinity of the groundwater is not affected by seasonal seawater intrusion. Schoeller diagrams illustrate that sulphate concentrations of the groundwater of the second aquifer are relatively low compared to those of the recent seawater. Therefore, this result suggests that the brackish water in the second aquifer is probably from ancient seawater that was trapped within the sediments for a long period of time, rather than due to direct seawater intrusion.