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.     

Exploring the contribution of social capital in building resilience for climate change effects in peri-urban areas,Dar es Salaam,Tanzania

GeoJournal - This article aims to contribute to the literature on the quest for resilient cities by focusing on the climate change resilience building discourse in peri-urban areas, and...     

Similarity-based approaches in hydrogeology: proposal of a new concept for data-scarce groundwater resource characterization and prediction

Hydrogeology Journal - A new concept is proposed for describing, analysing and predicting the dynamic behaviour of groundwater resources based on classification and similarity. The concept makes...     

Modeling seawater intrusion in overexploited aquifers in the absence of sufficient data: application to the aquifer of Nea Moudania,northern Greece

In many coastal areas, overexploitation of groundwater resources has led both to the quantitative degradation of local aquifers and the deterioration of groundwater quality due to seawater intrusion. To investigate the behavior of coastal aquifers under these conditions, numerical modeling is usually implemented; however, the proper implementation of numerical models requires a large amount of data, which are often not available due to the time-consuming and costly process of obtaining them. In the present study, the investigation of the behavior of coastal aquifers under the lack of adequate data is attempted by developing a methodological framework consisting of a series of numerical simulations: a steady-state, a false-transient and a transient simulation. The sequence and the connection between these simulations constitute the backbone of the whole procedure aimed at adjusting the various model parameters, as well as obtaining the initial conditions for the transient simulation. The validity of the proposed methodology is tested through evaluation of the model calibration procedure and the estimation of the simulation errors (mean error, mean absolute error, root mean square error, mean relative error) using the case of Nea Moudania basin, northern Greece. Furthermore, a sensitivity analysis is performed in order to minimize the error estimates and thus to maximize the reliability of the models. The results of the whole procedure affirm the proper implementation of the developed methodology under specific conditions and assumptions due to the lack of sufficient data, while they give a clear picture of the aquifer’s quantitative and qualitative status.     

Estimating potential for exploitation of karst aquifer: case example on two Serbian karst aquifer

Rational use of karst water resources depends considerably on the existence of a successful approach to estimating karst groundwater potential for exploitation. The term groundwater potential for exploitation involves the application of regulatory measures affecting natural outflow regime. Generally, the concept of controlling the artificial discharge regime is based on two principles. The first one, the “borrowing” of water from deeper (storage) part of storage, is based mainly on the construction of groundwater wells. The second one, considers increasing dynamic reserves by providing additional artificial water storage thereby controlling the groundwater discharge regime, is based on the construction of underground dams. This paper examines how these two concepts can be used to increase potential groundwater usage from two karst springs in Serbia—St. Petka and Peru?ac. The ability to potentially increase the amount of karst ground water available for future exploitation may facilitate and give important directions to the future research.     

Understanding arsenic behavior in carbonate aquifers: Implications for aquifer storage and recovery (ASR)

Geochemical reactive transport modeling was coupled to bench-scale leaching experiments to investigate and verify the mobilization of geogenic arsenic (As) under a range of redox conditions from an arsenic-rich pyrite bearing limestone aquifer. Modeling and experimental observations showed similar results and confirmed the following: (1) native groundwater and aquifer matrix, including pyrite, were in chemical equilibrium, thus preventing the release of As due to pyrite dissolution under ambient conditions; (2) mixing of oxygen- and nitrate-rich surface water with oxygen-depleted native groundwater changed the redox conditions and promoted the dissolution of pyrite, and (3) the behavior of As along a flow path was controlled by a complex series of interconnected reactions. This included the oxidative dissolution of pyrite and simultaneous sorption of As onto neo-formed hydrous ferric oxides (HFO), followed by the reductive dissolution of HFO and secondary release of adsorbed As under reducing conditions. Arsenic contamination of drinking water in these systems is thus controlled by the re-equilibration of the system to more reducing conditions rather than a purely oxidative process.     

Methodology to evaluate the renewal period of carbonate aquifers: a key tool for their management in arid and semiarid regions,with the example of Becerrero aquifer,Spain

A methodological procedure is proposed for determining the renewal period (RP), which expresses the ratio of total storage to recharge of carbonate aquifers, and it was applied to the overexploited moderate-size Becerrero carbonate aquifer (southern Spain). To this end, geological and subsurface data—time domain electromagnetic (TEM) soundings and borehole logs—were integrated to construct a three-dimensional (3D) geological model of the aquifer. The interconnected porosity was estimated by analyzing 73 rock samples. The resulting 3D geometrical model makes it possible to quantify the fractions of the aquifer having a confined or unconfined behaviour. Based on the total storage capacity (179?·?10⁶–514?·?10⁶ m³) and available aquifer recharge estimation (4.8?·?10⁶–6.4?·?10⁶ m³/year), an RP between 37 and 106 years is obtained. In view of the RP, an exploitation rate slightly lower than the average recharge of the system is recommended, so that the piezometric level will be stable but below the discharge head that is produced through the springs in natural conditions. The proposed methodology to obtain an aquifer RP and the management strategies designed accordingly are of broad interest, especially for carbonate aquifers, which are abundant in arid and semiarid regions.     

Development of a hydrogeological conceptual wetland model in the data-scarce north-eastern region of Kilombero Valley,Tanzania

Understanding groundwater/surface-water interactions in wetlands is crucial because wetlands provide not only a high potential for agricultural production, but also sensitive and valuable ecosystems. This is especially true for the Kilombero floodplain wetland in Tanzania, which represents a data-scarce region in terms of hydrological and hydrogeological data. A comprehensive approach combining hydrogeological with tracer-based assessments was conducted, in order to develop a conceptual hydrogeological wetland model of the area around the city of Ifakara in the north-eastern region of Kilombero catchment. Within the study site, a heterogeneous porous aquifer, with a range of hydraulic conductivities, is underlain by a fractured-rock aquifer. Groundwater chemistry is mainly influenced by silicate weathering and depends on groundwater residence times related to the hydraulic conductivities of the porous aquifer. Groundwater flows from the hillside to the river during most of the year. While floodwater close to the river is mainly derived from overbank flow of the river, floodwater at a greater distance from the river mainly originates from precipitation and groundwater discharge. Evaporation effects in floodwater increase with increasing distance from the river. In general, the contribution of flood and stream water to groundwater recharge is negligible. In terms of an intensification of agricultural activities in the wetland, several conclusions can be drawn from the conceptual model. Results of this study are valuable as a base for further research related to groundwater/surface-water interactions and the conceptual model can be used in the future to set up numerical flow and transport models.     

Recovery scenarios for deep over-exploited aquifers with limited recharge: methodology and application to an aquifer in Belgium

Recharge of deep-seated aquifers must be provided through leakage from overlying geological formations and can be very limited. Although hydraulic characteristics of these aquifers may be favorable for extensive exploitation, the renewability of the groundwater resource may be very restricted and only significant on long time scales. Over-exploitation of such aquifers leads to steadily declining piezometric levels and water balancing and steady state conditions are not to be expected on the short term. Recovery of such systems is very difficult, also because of the socio-economic dependencies of the water resource and the long time it requires to replenish the system. Management plans for restoration should be based on the transient hydrodynamical behaviour. A recovery plan should be developed based on results of a groundwater flow model. In a first step must be simulated how long it will take for the aquifer system to fully recover by natural recharge. As this will most often be a very long time it can be cancelled out as a realistic scenario. An alternative objective is to be formulated which could be achieved in the near future on a time scale of years or decades. Using model simulations, reduction schemes of exploitation rates shall be quantified that will allow to reach these defined goals, such as raising piezometric levels above the top of the aquifer itself or above the top of the overlying geological unit. The methodology is illustrated with an example of a deep hard rock aquifer in Belgium where piezometric levels have dropped below the top. The objective here was to raise the levels again above the top within the next 50 years. However, this is not accomplished by balancing aquifer recharge and discharge, but is reached by a lateral redistribution of the water over the aquifer extent. A full recovery of this aquifer would require at least a century with only a very limited amount of exploitation.     

Radial Dupuit interface flow to assess the aquifer storage and recovery potential of saltwater aquifers

A new accurate numerical solution is presented for aquifer storage and recovery (ASR) systems in coastal aquifers; flow is approximated as radial Dupuit interface flow. The radial velocities of points on the interface are a function of time, the vertical coordinate, and the dimensionless parameter D (the discharge of the well divided by the product of the hydraulic conductivity, the square of the aquifer thickness, and the dimensionless density difference). The recovery efficiency of an ASR system (the ratio of the recovered volume of water divided by the injected volume of water) is determined by D and by the relative lengths of the injection, storage and recovery periods. Graphs are produced for the recovery efficiency as a function of parameter D for ASR operations with and without storage periods and for multiple cycles. The presented solutions and graphs are to be used as screening tools to assess the feasibility of specific injection, storage and recovery scenarios of planned ASR systems in saltwater aquifers without having to run complicated flow and transport models. When the screening tool indicates that recovery efficiencies are acceptable, the consideration of other features such as mixing and chemistry is warranted.     

Use of molecular approaches in hydrogeological studies: the case of carbonate aquifers in southern Italy

Waterborne pathogens represent a significant health risk in both developed and developing countries with sensitive sub-populations including children, the elderly, neonates, and immune-compromised people, who are particularly susceptible to enteric infections. Annually, approximately 1.8 billion people utilize a faecally contaminated water source, and waterborne diseases are resulting in up to 2.1 million human mortalities globally. Although groundwater has traditionally been considered less susceptible to contamination by enteric pathogens than surface water due to natural attenuation by overlying strata, the degree of microbial removal attributable to soils and aquifers can vary significantly depending on several factors. Thus, accurate assessment of the variable presence and concentration of microbial contaminants, and the relative importance of potentially causative factors affecting contaminant ingress, is critical in order to develop effective source (well) and resource (aquifer) protection strategies. “Traditional” and molecular microbiological study designs, when coupled with hydrogeological, hydrochemical, isotopic, and geophysical methods, have proven useful for analysis of numerous aspects of subsurface microbial dynamics. Accordingly, this overview paper presents the principal microbial techniques currently being employed (1) to predict and identify sources of faecal contamination in groundwater, (2) to elucidate the dynamics of contaminant migration, and (3) to refine knowledge about the hydrogeological characteristics and behaviours of aquifer systems affected by microbial contamination with an emphasis on carbonate aquifers, which represent an important global water supply. Previous investigations carried out in carbonate aquifers in southern Italy are discussed.     

Experimental observations of aquifer storage and recovery in brackish aquifers using multiple partially penetrating wells

Hydrogeology Journal - Aquifer storage and recovery systems using multiple partially penetrating wells (MPPW-ASR) can form a viable solution to the problem of freshwater buoyancy when using...     

Development of an efficient surrogate model based on aquifer dimensions to prevent seawater intrusion in anisotropic coastal aquifers,case study: the Qom aquifer in Iran

Coastal aquifers are usually exposed to saltwater intrusion. Therefore, groundwater extracted from these aquifers should be regulated considering their dimensions and effective parameters. In this paper, optimum discharge from a large number of exploitation wells is evaluated according to variations of width, length, and anisotropy coefficient in the Qom aquifer near the salt lake in central Iran. First, the wells are divided into clusters to decrease the number of decision variables. Then, the location and discharge from each cluster is obtained using SEAWAT and charged system search (CSS) simulation–optimization model with the assumption of three-dimensional variable density flow. The maximum discharge considering various anisotropy rates is computed based on different values of lengths and widths of the aquifer. Finally, an M5-tree model is trained using the obtained samples to derive a linear relationship between input and output data. Based on the results, for various ranges of width and length of an aquifer with impermeable boundaries, different linear equations for optimum discharge are obtained. Also, it was found that for an aquifer with a small width, the critical discharge is a function of the length while the effect of the boundaries is negligible. Sensitivity analysis of the anisotropy coefficient reveals that with increasing the anisotropy rate, thickness and slope of the transition zone decrease and as the maximum discharge increases consequently. However, the sensitivity of the discharge to anisotropy rate is not remarkable. A comparison between the results of this study with those of the analytical method based on sharp interface assumption is carried out. For the critical condition, the best agreement between analytical equation (\(\overline {L} =0.87\overline {W} +0.62\)) and proposed method (\(\overline {L} =0.83\overline {W} - 1.41\)) is achieved for the anisotropic aquifer when the 50% isochlor is assumed as the measure of salt water intrusion.     

Threshold behavior of karst aquifers: the example of the Lurbach karst system (Austria)

Threshold behavior in hydrological systems generally involves a qualitative change of a single process, the system response or the functioning of the system. Different types of thresholds and their underlying controls are examined using the example of the Lurbach karst system (Austria). This karst system receives allogenic recharge from the sinking stream Lurbach, which under low-flow conditions only resurges at the Hammerbach spring. Yet, under medium- to high-flow conditions an overflow toward another spring, the Schmelzbach outlet occurs. Thresholds in physicochemical spring responses and their underlying controls are identified from the analysis of heat and solute transport processes in karst conduits. Applying this concept to the Hammerbach spring suggests that the threshold controlling the response of the spring water temperature was changed in the time period from 2006 to 2009 relative to the years before. At the same time, changes are observed in the behavior of the spring hydrograph and the discharge threshold at which the overflow to the Schmelzbach system is activated. All of these observations can be consistently explained by a decreased diameter of the conduit pathways within the indicated time period, presumably caused by the redistribution of sediments due to a flood event in 2005. Thus, thresholds in the physicochemical spring response were successfully employed to support the identification of a change in the functioning of the Lurbach karst system, which occurred possibly because a threshold related to the sediment transport within the karst conduits was crossed.     

Integration of geoelectric and hydrochemical approaches for delineation of groundwater potential zones in alluvial aquifer

Geoelectric and hydrochemical approaches are employed to delineate the ground-water potential zones in District Okara, a part of Bari Doab, Punjab, Pakistan. Sixty-seven VES surveys are conducted with the Electrical Resistivity Meter. The resultant resistivity verses depth model for each site is estimated using computer-based software IX1D. Aquifer thickness maps and interpreted resistivity maps were generated from interpreted VES results. Dar-Zarrouk parameters, transverse resistance (TR), longitudinal conductance (SL) and anisotropy (λ) were also calculated from resistivity data to delineate the potential zones of aquifer. 70% of SL value is ≤3S, 30% of SL value is > 3S. According to SL and TR values, the whole area is divided into three potential zones, high, medium and low potential zones. The spatial distribution maps show that north, south and central parts of study area are marked as good potential aquifer zones. Longitudinal conductance values are further utilized to determine aquifer protective capacity of area. The whole area is characterized by moderate to good and up to some extent very good aquifer protective area on the basis of SL values. The groundwater samples from sixty-seven installed tube wells are collected for hydro-chemical analysis. The electrical conductivity values are determined. Correlation is then developed between the EC (μS/cm) of groundwater samples vs. interpreted aquifer resistivity showing R2 value 0.90.     

Strontium,SO42−/Cl− and Mg2+/Ca2+ ratios as tracers for the evolution of seawater into coastal aquifers: the example of Castell de Ferro aquifer (SE Spain)

The strontium content and the SO₄^2?/Cl^? and Mg²⁺/Ca²⁺ ratios were used as natural tracers of the residence time of seawater intrusion into the Castell de Ferro aquifer. Analysis of these parameters indicated the existence of two principal flowpaths in the aquifer. The first flows through the eastern part of the aquifer, through the karstified Castell de Ferro massif; it accommodates a larger and more rapid flow, so that the residence time is shorter, leading to lower SO₄²⁺/Cl^? ratios, lower Sr²⁺ content and higher Mg²⁺/Ca²⁺ ratios. The second flowpath is in the western sector, and flows exclusively through alluvial deposits; the flow here is slower, particularly that flowing towards the sea. Thus the residence time of the water here will be longer and there is scant flushing of the intruded seawater; this is manifested in the high Sr²⁺ content, high SO₄²⁺/Cl^? and low Mg²⁺/Ca²⁺ ratios. To cite this article: P. Pulido-Leboeuf et al., C. R. Geoscience 335 (2003).     

Relative importance and chemical effects of diffuse and focused recharge in an eogenetic karst aquifer: an example from the unconfined upper Floridan aquifer, USA

Karst aquifer studies often focus on allogenic water inputs and large conduit flow. However, diffuse recharge can be significant, particularly in unconfined eogenetic karst aquifers that retain high matrix permeability. This study examines an unconfined region of the upper Floridan aquifer (USA) that hosts a sinking stream, its resurgence, and a large conduit system. Daily diffuse recharge was approximated using a water-budget method and ranged from 17% of precipitation during a low precipitation year to >53% during the highest precipitation year, illustrating the highly variable nature of diffuse recharge in this region. The total allogenic input via the sinking stream over the 5 years of the study was significantly larger than the volume of diffuse recharge. However, only about 2% of the allogenic recharge flows from the conduit into the surrounding aquifer. That flow is restricted to storm events when hydraulic heads in the conduits exceed those in the surrounding aquifer. The estimated volume of dissolution is similar for allogenic recharge and diffuse recharge to the unconfined region surrounding the conduits, but dissolution from the diffuse recharge is distributed over a larger area than dissolution from allogenic recharge. These results exemplify how recharge type impacts flow and water–rock interactions in eogenetic karst aquifers.     

Management of saltwater intrusion in coastal aquifers using different wells systems: a case study of the Nile Delta aquifer in Egypt

Hydrogeology Journal - Saltwater intrusion (SWI) is a type of pollution that adversely affects the quality of groundwater in coastal aquifers. The Nile Delta aquifer (NDA) in Egypt contains a large...