A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

Proposed method for groundwater vulnerability mapping in carbonate (karstic) aquifers: the COP method

The ‘COP method’ has been developed for the assessment of intrinsic vulnerability of carbonate aquifers in the frame of the European COST Action 620. This method uses the properties of overlying layers above the water table (O factor), the concentration of flow (C factor) and precipitation (P factor) over the aquifer, as the parameters to assess the intrinsic vulnerability of groundwater. This method considers karst characteristics, such as the presence of swallow holes (C factor) and their catchment areas as well as karstic landforms, as factors which decrease the natural protection provided by overlying layers (O factor). The P factor allows for consideration of the spatial and temporal variability of precipitation, which is considered the transport agent of contamination. Two carbonate aquifers in the South of Spain, Sierra de Líbar (a conduit flow system) and Torremolinos (a diffuse flow system), have been selected for the application and validation of the method and the results have been compared with three methods widely applied in different aquifers around the world (AVI, GOD and DRASTIC). Comparisons with these methods and validation tools (hydrogeological data and tracer test) show the advantages of the COP method in the assessment of vulnerability of karstic groundwaters.     

Comparison of the AVI,modified SINTACS and GALDIT vulnerability methods under future climate-change scenarios for a shallow low-lying coastal aquifer in southern Finland

A shallow unconfined low-lying coastal aquifer in southern Finland surrounded by the Baltic Sea is vulnerable to changes in groundwater recharge, sea-level rise and human activities. Assessment of the intrinsic vulnerability of groundwater under climate scenarios was performed for the aquifer area by utilising the results of a published study on the impacts of climate change on groundwater recharge and sea-level rise on groundwater–seawater interaction. Three intrinsic vulnerability mapping methods, the aquifer vulnerability index (AVI), a modified SINTACS and GALDIT, were applied and compared. According to the results, the degree of groundwater vulnerability is greatly impacted by seasonal variations in groundwater recharge during the year, and also varies depending on the climate-change variability in the long term. The groundwater is potentially highly vulnerable to contamination from sources on the ground surface during high groundwater recharge rates after snowmelt, while a high vulnerability to seawater intrusion could exist when there is a low groundwater recharge rate in dry season. The AVI results suggest that a change in the sea level will have an insignificant impact on groundwater vulnerability compared with the results from the modified SINTACS and GALDIT. The modified SINTACS method could be used as a guideline for the groundwater vulnerability assessment of glacial and deglacial deposits in inland aquifers, and in combination with GALDIT, it could provide a useful tool for assessing groundwater vulnerability to both contamination from sources on the ground surface and to seawater intrusion for shallow unconfined low-lying coastal aquifers under future climate-change conditions.     

A combined methodology to assess the intrinsic vulnerability of aquifers to pollution from agrochemicals

The groundwater vulnerability indices are valuable tools for the development of agrochemicals management strategies based on environmental/agricultural policies. The groundwater vulnerability methods of LOS, SINTACS, DRASTIC, Pesticide DRASTIC, GOD and AVI were applied for the agricultural fields of Sarigkiol basin (Northern Greece). The results of the aforementioned methods were examined and discussed in order to show how the dissimilarities in the vulnerability assessment approaches may become an advantage. The results of the methods were used to propose a combined conceptual approach which adds another two dimensions (depth and time) in the current two-dimensional vulnerability mapping (longitude, latitude) procedures. The LOS method provided information about the intrinsic vulnerability of the topsoil (30 cm) to water (+conservative pollutants) and nitrogen losses, and the AVI method described the vulnerability of the unsaturated zone to allow pollutants to reach the aquifer while the aquifer vulnerability was analysed using SINTACS, DRASTIC, Pesticide DRASTIC and GOD. In this study, the results of the SINTACS method were found more accurate to describe the local aquifer conditions. The final conceptual approach provided a stratified vulnerability (dimension of depth) of the overall hydrogeologic system using LOS for the topsoil, AVI for unsaturated zone and SINTACS for the aquifer. The dimension of time was introduced by the LOS and AVI methods, which provide quantitative results in time. The use of LOS method also highlighted the basic limitation of the other methods to describe the potential contribution to pollution of areas (especially upland areas) which are out of the aquifer boundaries.     

Groundwater temperature and electrical conductivity as tools to characterize flow patterns in carbonate aquifers: The Sierra de las Nieves karst aquifer, southern Spain

In carbonate massifs, flow patterns are conditioned by karstification processes which develop a conduit network and preserve low permeability microfractured blocks. The Sierra de las Nieves karst massif (southern Spain) is subjected to a given climatic and geological context, and thus it is possible to analyse the spatial and temporal variability of the water temperature and electrical conductivity at its main karst outlets, which display different responses to rainfall episodes. In this experimental field area, conduit flow and diffuse flow drainage patterns have been distinguished by combining groundwater temperature and electrical conductivity data. Both parameters show large variations in water coming from conduit flow systems and low variations in water drained by springs draining diffuse flow systems. However, groundwater temperature displays the smallest variations, which seems to indicate that this parameter is less sensitive as regards characterising the degree of karstification, which is a key question in characterising the aquifer functioning.     

Comparative application of two methods (COP and PaPRIKa) for groundwater vulnerability mapping in Mediterranean karst aquifers (France and Spain)

A comparative test of two vulnerability mapping methods (COP and PaPRIKa) specifically dedicated to for karst aquifers was carried out on two Mediterranean carbonate aquifers. The vulnerability maps obtained for each aquifer present important differences. To identify and determine the origin of these differences, the results were statistically analyzed using sensitivity analysis, coefficients of determination and scatter graphs. In addition, the global vulnerability (Gv) parameter was used to measure the general vulnerability of the aquifer and to compare the results obtained. This statistical analysis led us to conclude that the main cause of differences between these two methods used to assess aquifer vulnerability lie in the relative importance of the parameters employed in calculating the vulnerability index. For the PaPRIKa method, the variable related to infiltration (slope and karst features) has the most influence, with less weight being assigned to the protective capacity of layers overlying the aquifer. For the COP method, the most influent variable is defined by the layers overlying the aquifer, together with infiltration characteristics, determined by the relative importance of different forms of infiltration in each aquifer. The vulnerability mappings performed using the COP method present greater coherence with the known hydrogeological behavior of the study areas, especially the Spanish aquifers. Nevertheless, further hydrogeological investigations are needed, such as ones to validate the obtained vulnerability maps.     

Groundwater recharge in semi-arid carbonate aquifers under intensive use: the Estepa Range aquifers (Seville, southern Spain)

Quantifying groundwater recharge in carbonate aquifers located in semi-arid regions and subjected to intensive groundwater use is no easy task. One reason is that there are very few available methods suitable for application under such climatic conditions, and moreover, some of the methods that might be applied were originally designed with reference to non-carbonate aquifers. In addition, it is necessary to take into account the fact that, in any given aquifer, groundwater recharge is modified by the groundwater exploitation. Here we focus on four methods selected to assess their suitability for estimating groundwater recharge in carbonate aquifers affected by intensive exploitation. The methods were applied to the Estepa Range aquifers of Seville, southern Spain, which are subjected to different degrees of exploitation. Two conventional methods were used: chloride mass balance and daily soil–water balance. These results were compared with the results obtained by means of two non-conventional methods, designed for application to the carbonate aquifers of southern Spain: the APLIS and ERAS methods. The results of the different methods are analogous, comparable to those obtained in nearby non-exploited carbonate aquifers, confirming their suitability for use with carbonate aquifers in either natural or exploited regimes in a semi-arid climate.     

Using non-conservative tracers to characterise karstification processes in the Merinos-Colorado-Carrasco carbonate aquifer system (southern Spain)

The systematic sampling of the chemical composition of the groundwater from five karst springs (including an overflow spring) and one outflowing borehole have permitted to determine distinctive chemical changes in the waters that reflect the geochemical processes occurring in a carbonate aquifer system from southern Spain. The analysis of the dissolution parameters revealed that geochemical evolution of the karst waters basically depends on the availability of the minerals forming aquifer rocks and the residence time within the aquifers. In the three proposed scenarios in the aquifers, which include the preferential flow routines, the more important geochemical processes taking place during the groundwater flow from the recharge to the discharge zones are: CO₂ dissolution and exsolution (outgassing), calcite net dissolution, calcite and dolomite sequential dissolution, gypsum/anhydrite and halite dissolution, de-dolomitization and calcite precipitation. A detailed analysis of the hydrochemical data set, saturation indices of the minerals and partial pressure of CO₂ in the waters joined to the application of geochemical modelling methods allowed the elaboration of a hydrogeochemical model of the studied aquifers. The developed approach contributes to a better understanding of the karstification processes and the hydrogeological functioning of carbonate aquifers, the latter being a crucial aspect for the suitable management of the water resources.     

Hydrogeology of karst aquifer systems in SW Turkey and an assessment of water quality and contamination problems

The aim of this research is to determine the relationship between groundwater flow and water quality of different ground and surface water basins in the southwest Turkey. In addition, groundwater vulnerability is assessed taking into consideration groundwater flow and quality. The autochthonous Beydaglari limestone is the major karstic aquifer in the region. According to the groundwater level map of alluvium aquifers in the basins, groundwater discharge toward the carbonate aquifer is direct and indirect. The hydrogeological connection between ground and surface water basins occurs via the karstic aquifer located at the bottom of the alluvium bottom. In Egirdir lake, water also discharges in the karstic aquifer via karstic sinkholes at the western border of the lake. In the research area, general groundwater discharge is toward the Mediterranean Sea by means of autochthonous carbonate system, according to hydrogeological investigations, research of lineament and hydraulic conductivities. This result is supported by the locations of lineaments and shore springs discharging from the limestone. In addition, spreading of contaminants via karstic aquifer to great distance has been clearly identified.     

Methodology for groundwater recharge assessment in carbonate aquifers: application to pilot sites in southern Spain

Aquifer recharge can be determined by conventional methods such as hydrodynamic or hydrologic balance calculations, or numerical, hydrochemical or isotopic models. Such methods are usually developed with respect to detrital aquifers and are then used on carbonate aquifers without taking into consideration their hydrogeological particularities. Moreover, such methods are not always easy to apply, sometimes requiring data that are not available. Neither do they enable determination of the spatial distribution of the recharge. For eight regions in southern Spain, the APLIS method has been used to estimate the mean annual recharge in carbonate aquifers, expressed as a percentage of precipitation, based on the variables altitude, slope, lithology, infiltration landform, and soil type. The aquifers are representative of a broad range of climatic and geologic conditions. Maps of the above variables have been drawn for each aquifer, using a geographic information system; thus they can be superimposed to obtain the mean value and spatial distribution of the recharge. The recharge values for the eight aquifers are similar to those previously calculated by conventional methods and confirmed by discharge values, which corroborates the validity of the method.     

Coastal aquifers of Europe: an overview

Europe has a long coastline along which many human activities are concentrated. Coastal aquifer development is often intensive and subject to salinity problems as a result of seawater intrusion, upconing of deep saline water, and residual salinity in aquitards. Studies have been carried out since the end of the nineteenth century, and have set the basic principles of freshwater/marine-water relations. Some aquifers are carefully studied, monitored and managed; examples can be found in The Netherlands, Belgium, southern Italy and northeastern Spain. Often important coastal aquifers are of relatively small size but crucial to supply human needs, which vary from mostly urban and industrial in mid to high latitudes, to dominantly seasonal for irrigation in the Mediterranean area and the archipelagos. The different situations can be grouped into: recent deltas, other coastal detritic formations, carbonate formations, and the special situations in small islands. Comments on these situations are provided and some examples are introduced, just for illustration.     

Source vulnerability mapping in carbonate (karst) aquifers by extension of the COP method: application to pilot sites

A step from resource to source vulnerability mapping is presented, based on the European COST Action 620 approach for karst groundwater protection. Guidelines on vulnerability assessment of the horizontal groundwater flow path within the karst saturated zone (K factor) are proposed. By integrating this into the previously existing COP method for intrinsic resource vulnerability mapping, adequate source protection can be assessed. The proposed “Karst saturated zone (K) factor” assessment considers groundwater travel time (t subfactor), connection and contribution to the source (r subfactor) and active conduit or fissured network (n subfactor). The extended COP method was applied in two carbonate aquifers in southern Spain with different geological, hydrogeological and climate settings. The results are coherent with previous research results of the studied areas. Moreover, they are consistent with the occasional groundwater contamination detected in one of the springs. On the other hand, an absence of contamination, despite high risk, justifies the lower degrees of vulnerability assigned to the sources surveyed. The source vulnerability maps obtained can thus be used as a basis for the delineation of protection zones.     

Intrinsic groundwater vulnerability assessment: issues,comparison of different methodologies and correlation with nitrate concentrations in NW Italy

The concept of aquifer vulnerability is certainly useful in the field of groundwater protection. Nevertheless, within the scientific community, the definition of vulnerability is still under debate and lacks standardisation. As a consequence, the methods for evaluating the vulnerability degree are numerous and often lead to conflicting results. Thus, in this study, three methods that are commonly used in groundwater vulnerability assessments due to their easy application (namely DRASTIC, GOD and TOT) were utilised in four areas of the Piedmont region (NW Italy). The results obtained by the different methods were compared and correlated with the nitrate concentrations in the groundwater. The aims of the study were (i) to evaluate the effectiveness of the adopted methods and their comparability, (ii) to discuss the limits of the intrinsic vulnerability methods and (iii) to verify the applicability of nitrate as a tracer in the assessment of groundwater vulnerability or explain the reasons why it is not applicable. It was observed that the three intrinsic vulnerability methods are not able to uniquely identify the most or least vulnerable areas. Additionally, the comparison of the intrinsic vulnerability indexes only occasionally showed a reasonable correlation. Furthermore, there was no clear correlation between the vulnerability indexes and nitrate concentrations in the groundwater. These results could be explained by several reasons: (1) the methods are mostly based on the level of groundwater protection provided by the overlaying lithologies and do not consider the physical processes occurring in the aquifer; (2) the intrinsic vulnerability methods only consider vertical pathways for contaminants, but a pre-existing contaminant could be present in the aquifer; (3) groundwater nitrate concentrations are affected by the nitrate input and surplus; and (4) nitrates are subject to physical and biological attenuation in aquifers and cannot necessarily be considered stable tracers in the assessment of groundwater vulnerability.     

An integrative geological and geophysical approach to characterize a superficial deltaic aquifer in the Camargue plain,France

Deltaic aquifers are complex due to the important heterogeneity of their structure and their hydrogeological functioning. Auger drilling provides localized, but very robust geological and hydrogeological information, while geophysical surveys can provide integrated subsurface information. An effective, easy-to-use and low-cost methodology combining geological/hydrogeological information from Auger drillings and the results from three geophysical techniques (Electromagnetic mapping, Electrical Sounding and Electrical Resistivity Tomography) is being developed to characterize the structure of a typical Mediterranean deltaic aquifer. A first hypothesis about hydrodynamic properties of the aquifer is also obtained. The study area is located in the Rhone delta (Middle Camargue/southern France). Integration of geophysical and geological techniques allowed identifying the presence and lateral extension of the Saint-Ferréol paleochannel, the vertical contact between lagoonal–fluvial deposits and the marine clayed silt that separate the superficial aquifer from the deeper aquifer. Likewise, high north–south heterogeneity and east–west homogeneity were highlighted in the study area. Presence of clay in sandy deposits in the low areas implies changes in lateral hydraulic permeability. This fact, jointed to the low hydraulic gradient, suggests a slow groundwater flow in the local system. The Rhone delta presents a typical configuration of a Mediterranean deltaic aquifer, thus this methodological approach can be used for similar deltaic Mediterranean systems.     

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.     

Adapting existing experience with aquifer vulnerability and groundwater protection for Africa

Whilst groundwater vulnerability mapping and the delineation of resource and source protection zones have become an appropriate set of management tools for Britain as incorporated in European policy, much of semi-arid Africa is still dealing with more pressing issues centred on water supply coverage. There are a number of fundamental differences with Britain which disallow conventional vulnerability mapping and land zonation in much of Africa. Firstly, the scale of groundwater occurrence in weathered basement aquifers does not encourage vulnerability mapping to be undertaken at a field scale, whereas the Karoo and some of the larger areas of unconsolidated sedimentary aquifers could more readily be zoned according to aquifer vulnerability. Secondly, analysis needs to disregard the productivity (or recharge potential) of the aquifer so that poorly productive but socially important aquifers can be assessed. Such practice also avoids the need to identify a value for effective rainfall, a problematical value in semi-arid climates given current uncertainties. Thirdly, it is difficult to protect the many small and dispersed groundwater supply sources typical of many African communities when compared with the fewer and larger sources characteristic of Britain and Europe. Some aspects of European groundwater management practice have been transferred to the African context, notably by South Africa, but there are other practices that should not be attempted. Lessons from experience in South Africa highlight capacity to implement as a key inhibiting factor. Examples of vulnerability assessment and land zonation in a variety of African settings indicate only limited success so far with standard vulnerability assessment and land zonation techniques. Alternative means of quantifying the problem of the optimum proximity of the pit latrine from the well are highlighted, with a minimum separation of 10 m suggested for typical weathered basement rocks.     

Estimation of groundwater recharge in a major sand and gravel aquifer in Ireland using multiple approaches

Groundwater recharge was investigated in the most extensive sand and gravel aquifer (area of approximately 200 km²) in the Republic of Ireland as part of a wider study seeking to derive recharge estimates using aquifer vulnerability mapping. The proportion of effective rainfall (total rainfall minus actual evapotranspiration) that leads to recharge is known as the recharge coefficient. The recharge investigation involved a variety of approaches, including soil moisture budgeting, well hydrograph analysis, numerical modelling and a catchment water balance. The adoption of multiple techniques provided insights on recharge and also on aquifer properties. Comparison of two soil moisture budgeting approaches (FAO Penman-Monteith with Penman-Grindley) showed how variations in the effective rainfall values from these methods influence groundwater levels simulated in a numerical groundwater model. The catchment water balance estimated the recharge coefficient to be between 81 and 85%, which is considered a reasonable range for this aquifer, where overland flow is rarely observed. The well hydrograph analysis, using a previous estimate of specific yield (0.13), gave recharge coefficients in the range of 40–80%, considered low for this aquifer: a revised specific yield of 0.19 resulted in a more reasonable range of recharge coefficients of between 70 and 100%.     

Correlation between nitrate concentration in groundwater and parameters affecting aquifer intrinsic vulnerability

This paper is the result of a study which was carried out in order to verify if the traditional methods to evaluate the intrinsic vulnerability or vulnerability related parameters, are able to clarify the problem of nitrate pollution in groundwater. In particular, the aim was to evaluate limitations and problems connected to aquifer vulnerability methods applied to nitrate contamination prevision in groundwater. The investigation was carried out by comparing NO₃ ⁻ concentrations, measured in March and November 2004 in the shallow aquifer, and the vulnerability classes, obtained by using GOD and TOT methods. Moreover, it deals with a comparison between NO₃ ⁻ concentrations and single parameters (depth to water table, land use and nitrogen input). The study area is the plain sector of Piemonte (Northern Italy), where an unconfined aquifer nitrate contamination exists. In this area the anthropogenic presence is remarkable and the input of N-fertilizers and zootechnical effluents to the soil cause a growing amount of nitrates in groundwater. This approach, used in a large area (about 10,000 km²) and in several monitoring wells (about 500), allowed to compare the efficiency of different vulnerability methods and to verify the importance of every parameter on the nitrate concentrations in the aquifer. Furthermore it allowed to obtain interesting correlations in different hydrogeological situations. Correlations between depth to water table, land use and nitrogen input to the soil with nitrate concentrations in groundwater show unclear situations: in fact these comparisons describe the phenomenon trend and highlight the maximum nitrate concentrations for each circumstance but often show wide ranges of possible nitrate concentrations. The same situation could be observed by comparing vulnerability indexes and nitrate concentrations in groundwater. These results suggest that neither single parameters nor vulnerability methods (GOD and TOT) are able to describe individually the complex phenomena affecting nitrate concentrations in soil, subsoil and groundwater. In particular, the traditional methods for vulnerability analysis do not analyze physical processes in aquifers, such as denitrification and nitrate dilution. According to a recent study in the shallow unconfined aquifer of the Piemonte plain, dilution can be considered as the main cause for nitrate attenuation in groundwater.     

Assessing aquifer vulnerability to pollutants by electrical resistivity tomography (ERT) at a nitrate vulnerable zone in NE Spain

There is an increasing demand for groundwater vulnerability maps which illustrate the exposure of aquifers against pollution. These maps show areas of greatest potential for groundwater contamination on the basis of local subsurface conditions. Parameters affecting vulnerability are mainly permeability and thickness of each protective layer. For unconsolidated sediments, the permeability is strongly related to the clay content, which can be deduced from indirect resistivity methods, like electrical-imaging. Such geophysical methods can be of great help in groundwater vulnerability studies because they disturb neither the structure nor the dynamics of the soil. Sensibility analysis was performed of the electrical resistivity tomography method for accurately mapping soil media. Managers and public administrators may effectively use this method for assessing the potential risk of groundwater contamination. In the studied zone, electrical resistivity exhibits a wide range of variability that can be easily correlated to soil parameters, such as clay content and hydraulic conductivity. A numerical index of protection has been assessed from the geophysical information derived from 2D electrical resistivity tomography. This work represents a preliminary approach on the natural vulnerability evaluation of shallow aquifers at the Empordà basin (NE Spain) that is highly affected by diffuse pollution by nitrates.     

Groundwaters of the Central Ethiopian Rift: diagnostic trends in trace elements, δ<Superscript>18</Superscript>O and major elements

The Ethiopian Rift (a major portion of the Great East African Rift) is characterized by a narrow elongated depression bounded by highlands from both sides. This topographic configuration leads to a monsoon redistribution which resulted in an arid rift floor and humid high rainfall highlands. The rifting and associated volcanism also caused a thinning of the crust and facilitates influx of CO₂ and other mantle gases as diffuse sources or along faults from deeper sources. Groundwaters in the rift floor are usually of high mineral content (high F, U, As and salinity) while those on the plateau are of low mineral content. Among many factors, groundwater availability and quality in the rift floor aquifers is the function of their connection to the aquifers in the high rainfall plateau and the residence time of groundwater prior to reaching the rift floor. This entails the need for addressing one basic hydrologic question in such a setting: at what depth and rate does recharge from the high rainfall highland reach the lowland rift aquifers? This study uses spatial variations in trace elements and relates them to ¹⁴C variations, thereby investigating the suitability of using trace elements as proxies for residence time estimation of groundwaters of relatively short (1,000–2,000 years) residence time. This work also investigates the behavior of trace element trends along the groundwater flow path in a rifted setting and compares them with such trends in sedimentary aquifers elsewhere. The comparison shows a clear difference in behavior of trace elements along the groundwater flow path when compared with such variations in big sedimentary basins with no prominent rifting and volcanism, suggesting the need of calibrating the relation between trace elements and any direct residence time indicators. An integrated use of major elements, trace elements, and environmental isotopes reveals that the main recharge of the aquifers originates from mountain blocks and that recharge takes place via fractures with no evidence of evaporation prior to recharge. Redox processes appear to play a limited role in trace element geochemistry of groundwaters in the region. Progressive trends in trace element composition along the groundwater flow path suggest continuous groundwater flow from the plateau.