2.1 Study Area
The study area is located west of Lusaka City as shown in Fig. 1. The area comprises private farm lands and dotted densely-populated unplanned residential area lacking amenities and a reticulation system with a population of 36,532 (Ministry of Health, 2019). The area is characterized by a few boreholes with numerous springs mainly utilized for gardening and domestic use. The community heavily depends on groundwater sources for drinking, accessed through springs, shallow wells and boreholes.
The climate is mild and sub-tropical, (James et al., 2009) with three seasons: a hot dry season from August to October with temperatures between 26 and 38 ºC; a rainy season from November to April with temperatures between 27 and 34 ºC; and a cool dry season from May to August with temperatures from 13 to 20 ºC. The area lies in region II of the ecological regions of Zambia with average rainfall of between 800 and 1000 mm.
2.2 Geology and Hydrogeology of the Study Area
Lusaka lies topographically on a plateau of about 1,280 m above sea level (Bäumle & Kang'omba, 2012). It is underlain by a thick sequence of metasedimentary rocks of the Katanga Supergroup (Nkhuwa, 1996). Geologically, the lowest geological unit is the Chunga Formation (Fig. 2) that comprises crystalline rock of low permeability. The Cheta Formation then overlies the Chunga comprises mainly of schist and quartzite but with some more permeable dolomitic limestone, exploited to limited extent as a water source (Nyambe & Maseka, 2000). The Lusaka Formation is the uppermost with a dolomitic limestone and marble unit. With its complex interconnective of dissolution conduits that follow joint set planes and sheet fissures, it is the primary aquifer for Lusaka. It has a cutter and pinnacle epikarst surface topography extending to a depth of 25 m (Nkhuwa, 2003). Groundwater abstraction from the Lusaka Aquifer has led to a progressive decline in the water table elevation and is suggested to be overexploited (Mpamba et al., 2008). It has generally a shallow water table depths from 6 to 15 m.
Direct recharge is predominantly from the Lusaka Forest Reserve, southwest of Lusaka, an area undergoing rapid deforestation and settlement development. Substantial localized recharge has also be reported from sources such as overflowing sewers, adding to the public health risk (Karen et al., 2019). Recharge to the groundwater during the wet season from November to April averages 186 mm/yr, which is approximately 27% of the annual rainfall (De Waele & Follesa, 2003; Nyambe & Maseka, 2000). Groundwater flow direction is generally southeast to northwest through the Lusaka area but tends to be in northeast and southwest orientation in some parts (Museteka & Bäumle, 2009; Mpamba et al., 2008; Von Hoyer et al., 1978). These general, large-scale flow directions, however, are not definitive in a Karstic system and thus require detailed assessment as they are based on limited groundwater observations. Groundwater discharge is on the schist-dolomite contact through intermittent springs and in low-lying swampy areas (Museteka & Bäumle, 2009). Using the underlying geology, soil type and soil thickness, Kang'omba and Bäumle (2013) classified groundwater contamination vulnerability classes for the Lusaka area and its environs. Most of the area associated with the Lusaka Formation (Fig. 2) was designated as extremely vulnerable to groundwater pollution. Of great concern is that pit latrines and drinking water sources may be directly connected via highly conductive Karst conduits, a situation which can lead to bacterial and nitrate contamination and the spread of waterborne diseases (Reaver et al., 2021).
2.3 Dye tracing
2.3.1 Field setup
A reconnaissance survey was undertaken to identify potential dye-injection points (pit latrines) and observation points within the study area. The selection of the dosing points was done assuming that the springs (Fig. 3) represent focused discharge points for the groundwater therefore dosing points should be placed up stream (Benischke, 2021). The general methodological approach undertaken was qualitative over a 10 week maximum period. The observation in the dry season was done in between third week of September and second week of November 2018, while the observations for wet season was between second week of March and second week of April 2019. Dye tracing was done using Fluorescein and Rhodamine. Fluorescein is a reddish-brown powder that turns vivid yellow-green in water. It is photochemically unstable and its concentration decreases at a pH of less than 5.5 (Mull et al., 1988). However, Fluorescein has moderate to high detection, generally has a low sorptive tendency and diffusity (Benischke, 2021). Rhodamine on the other hand, is highly detectable with a low sorptive and diffusivity capacity. It turns vivid orange-pink when in water and is photochemically stable. Dye-recovery sites (springs) were selected based on the assumption that they will intersect groundwater flow from the injection points (pit latrines). Further, we also considered if the spring was perennial and accessibility throughout the year.
Dye receptors were placed in all of the area springs. They were first placed in the springs and sampled before any dye injection to establish background dye signals The receptors were made of wire mesh envelopes measuring about 8 x 10 cm containing two tablespoons of activated charcoal (Fig. 4). Sampling was done in both the wet and dry season. The quantity of dye used in the injection into pit latrines (Fig. 5) was estimated based on the estimated flow conditions, dilution and distance from the injection point to the remotest recovery point. A study conducted by Field (2003), has recommended several formulae for estimating quantities of injection. However, these formulae are based on a number of parameters such as average discharge, distance between injection and sampling place, expected time of the first detection and the peak of a dye breakthrough, desired peak concentration, specific coefficients for the type of tracer, and others parameters. Furthermore, these formulas were developed by the authors from experience and no one can cover all the various actual conditions during a test.
In this study, we used the approach recommended by Benischke, 2021 summarized as: ‘as much as necessary, as little as possible’. It was judiciously decided that 300 mg Fluorescein and 958 mg Rhodamine would be dosed in the pit latrines.
The initial charcoal receptors were deployed to measure background data and presence of compounds that fluoresce at similar wavelengths as rhodamine and fluorescein. These background receptors were left in place for four days prior to any dye injection. This was done to detect any presence of dyes and their concentrations. New receptors were replaced at the recovery points just after dye injection to detect the dye breakthroughs and taken to the laboratory for measurement and replaced with new receptors every four days
2.3.2 Laboratory Analysis
Once collected, the receptors were kept in clean Zip-loc ® plastic bags which were labeled with the time collected, site number, date of placement and collection. The receptors were kept in a cool place out of the light. An AquaFluor® handheld fluorimeter, was used to analyze the samples. This instrument has the capacity to measure both Rhodamine and Fluorescein by switching channels. The minimum detection ranges for Fluorescein and Rhodamine was 0.4 ppb with linear range of 1–400 ppb. Before analysis, the Aqua Fluor® fluorimeter was calibrated using blank samples and known standards of the dyes. An eluent composed of propanol (5 parts), distilled water (3 parts) and ammonia hydroxide (2 parts) was prepared as described by Mull et al, (1988). On collection of receptors, 5 g of the charcoal was removed from the mesh and placed in a beaker. 10 mL of the eluent was poured into the beaker to immerse the charcoal. The elution was then allowed to stand for one hour after which the eluent was poured into a cuvette for analysis. This was repeated three time to ensure replication of results for quality assurance. A detailed description of the method can be found in Currens, 2013. The tracer concentration results were then used to plot breakthrough curves to evaluate tracer transport in the karst aquifer system.
To calculate the average travel time using the measured concentrations per recovery site, and can be determined using two methods namely firstly by calculating the concentration against the distance from the dosing point. Secondly a formula is employed to calculate time of travel using concentrations Eq. 1 (Benischke, 2021, Chapra, 2008) as below:
\(t=\frac{\sum _{i=0}^{n-1 }{\left({C}_{i}{t}_{i}+{C}_{i+1}{t}_{i+1}\right)}_{{(t}_{i+1}-{ t}_{i}) }}{{\sum }_{i=0}^{n-1}{\left({C}_{i}+{C}_{i+1}\right)}_{{(t}_{i+1}-{t}_{i})}}\) Equation 1
Where: \(t\) = Average travel time
Ci = Concentration (mg/L)
ti = Discrete time when sampled (min)
n = Number of weeks
I = Batch numbers
Average groundwater velocity was then calculated by dividing the distance from the dosing pit latrine to the recovery points with the average time.