Introduction

Over the years, before the discovery of crude oil in Nigeria, Ekiti area was the least developed and neglected part of the former Western region of Nigeria despite the fact that the area was a major producer of cocoa that formed the backbone of the region’s economy. Most of the communities are predominantly rural with the exception of few township areas like Ado-Ekiti (the state capital) and others such as Ikere-Ekiti, Ise-Ekiti, Emure-Ekiti and Ikole-Ekiti. The communities are peasant farmers and depend mostly on shallow groundwater (hand dug wells) and surface waters (rivers, streams and springs) for sustainability of their farming activity and other domestic requirements. Surface waters in form of rivers and streams apart from being prone to contamination from anthropogenic activities are subject to seasonal variation as most of them dried off during the dry season. Hence, the most reliable source of water for the agricultural and domestic purposes in the study area is groundwater. However, the groundwater in the area were found within two major aquiferous units (the weathered and fractured layers) and tapped in form of shallow hand-dug wells and shallow hand-pump boreholes (Ademilua and Olorunfemi 2000).

Some of the existing wells were uncovered, unlined and in few cases even were covered with wooden plank or rusting iron sheet. Consequently, the incidences of waterborne diseases (guinea worm, typhoid fever, gastrointestinal diseases, poliomyelitis etc.,) were at times common in the area. Hence, there is need for hydrochemical assessment of the shallow groundwater system of the study area with particular reference to drinking and irrigation quality evaluation.

The groundwater was generally a product of infiltrated rainwater, channel precipitation, overland flow, subsurface baseflow, agricultural drains and even municipal effluents. In many cases, waters from the highlighted different sources usually have a distinct isotopic and geochemical composition. Increased knowledge about the geochemical processes that control groundwater chemical composition in weathered/fractured basement aquifer of Ekiti area could lead to improved understanding of hydrochemical systems in such area. Such improved understanding can contribute to effective management and utilization of the groundwater resources by clarifying relations among groundwater quality, aquifer lithology and recharge type. Furthermore, the spatial changes in the water types can equally be evaluated and used to trace the recharge source(s) to the groundwater system using stable isotopes of water (δ18O and δD). These stable isotopes (unlike other tracers of ‘water’ such as Cl which is still widely used) are ideal tracers of water sources and movement because they are constituents of the water molecules (Akiti 1986). The elements constituting the water molecule, hydrogen and oxygen, undergo changes in the various phases of hydrologic environments; atmosphere, hydrosphere, biosphere and the upper part of the earth’s crust. These changes results in isotopic fractionation to give isotopic “fingerprints” due to different physicochemical, biochemical, kinetic and thermodynamic effects (Das et al. 1998). Consequently this study also employed evaluations of stable isotopes to draw conclusions with respect to recharge source(s) and processes of the groundwater system in the study area.

Previously, there were a number of studies relating to chemical quality aspect of shallow groundwater systems in the Basement Complex area of southwestern Nigeria which focussed mostly on the township areas like Abeokuta, Ilesha, Ibadan and Ado-Ekiti (Loehnert 1981; Ako et al. 1990; Malomo et al. 1991; Sangodiyin 1993; Asubiojo et al. (1997); Olayinka et al. 1999; Abimbola et al. 1999; Tijani and Abimbola 2003; Olobaniyi et al. (2007); Omotoyinbo and Okafor 2008). Sangodiyin (1993) found poor bacteriological quality of groundwater close to sites of waste disposals in Abeokuta area, highlighting the potential for impacts from other pollutants. Asubiojo et al. (1997) carried out a survey of tap water and groundwater quality from various sources in Southern Nigeria and found that most elements investigated are present at concentrations below the respective WHO standard values with exception of Ba, Pb, Ni and Se.

Abimbola et al. (1999) identified the dominant ions in the groundwater of Abeokuta as Ca2+ and HCO3, and the water was characterized as Ca-HCO3 type having low total dissolved solids (TDS). Omotoyinbo and Okafor (2008) worked on the influence of rock mineralogy on subsurface water in Ado-Ekiti and concluded that dissolved minerals from rocks have great influence on the underground water property. The level of the dissolved minerals was within the acceptable limit for water quality standard values by WHO.

All these studies and many others have contributed immensely to groundwater development and management. However, most of the works were restricted to hydrogeochemistry involving major cations and anions with only few given consideration to heavy metals and other trace constituents. Substantial hydrogeological and hydrochemical research works have been carried out in the Basement Complex of southwestern Nigeria, but in the study area these research works were scanty and to a greater extent did not include isotope evaluation of the groundwater. In most cases, the understanding of the hydrochemical properties of the aquifers as a means of establishing the recharge regimes for the better management of the aquifer systems have not been fully exploited. This study integrates both hydrochemical and stable isotope studies to characterize shallow groundwater in a typical weathered basement aquifer terrain of Ekiti area, with the focus on assessment of water quality status, hydrochemical facies and associated geochemical processes as well as recharge source(s) and processes.

Location of the study area

The study area lies within latitudes 7°25′ and 7°52′N and longitudes 4°48′ and 5°30′E covering an approximate area of about 4,344 km2. Prominent towns and settlements in the area include Ado, Ikere, Ilawe, Ise, Emure, Ifaki, Aramoko, Ikogosi, Iyin, Igede, Ijan, Igbara-Odo, Ijero, Iworoko as well as other smaller settlements (Fig. 1).

Fig. 1
figure 1

Location Map of the study area showing sampling points

Geographically, the area falls within the rainforest in the southwestern part of Nigeria. The topography is generally undulating with most area lying above 250 m above sea level. The landscape is characterized by old plains, broken steep-sided outcrop dome rocks (Inselbergs) that may occur singularly or in crops of ridges. Such outcrops exist mainly in form of rugged hills at Ado-Ekiti (central part of study area) and Ikere-Ekiti in the southern part. The study area lies within the tropical climatic zone, characterized by alternating wet and dry season and drained by rivers such Elemi, Ureje, Ose, Ogbesse and their seasonal network of tributary and streams forming a somewhat dendritic drainage pattern (Fig. 1). The wet season was from April to October, dominated by heavy rainfall, while the dry season was from November to March. The mean annual rainfall was about 1,500 mm while the mean annual temperature ranged from 25 to 30 °C.

Geologic and hydrogeologic frameworks

Geologically, the study area lies within the Pre-Cambrian Basement Complex of SW-Nigeria which also constitutes part of the Pan-African mobile belt lying east of the West African Craton (Jones and Hockey 1964; Odeyemi 1976; Rahaman and Ocan 1988). The major rocks in the study area include five principal lithologic units; migmatite/migmatite gneiss, metasediment, Pan-African granite, charnockite and quartzite/quartz-schist (Fig. 2). The migmatite/migmatite gneiss rocks cover more than half of the study area, occurring as highly denuded hills of essentially fine-medium grained textures with closely spaced alternating bands of leucocratic and melanocratic minerals. They are characteristically low lying and conspicuously foliated with abundance of platy biotite minerals sandwiched into zones that are markedly distinguishable from the light-coloured quartzo-feldspartic portions. Quartzite occurs as ridge of steeply dipping, massively bedded fine-medium grained bodies or an environment characterized by the abundance of quartz rubbles. Granite and charnockite occur as minor intrusive bodies within the migmatite rock. The granite is essentially fine to medium grained and sometimes porphyritic texture acidic rocks. Pegmatite, aplite and quartz veins occur as minor intrusions into the granitic rocks. Field studies indicate that the charnockite occurs as dark greenish grey and greasy, medium to fine grained and weakly foliated with occasional gneissic banding. The rock outcrops are fresh, rough in surface appearance with little signs of weathering. The rock is in contact with migmatite and granite. With the exception of the small hills and inselbergs, the bedrocks are generally covered by weathered regolith composing mainly of clay and sandy materials which are capped by lateritic soil in places depending on the bedrock types.

Fig. 2
figure 2

Geological map of the study area

Hydrogeologically, like in other typical basement complex terrain, unweathered or fresh crystalline rocks have, on their own, hardly any potential in terms of groundwater occurrence and flow. However, the mode of groundwater occurrence in the study area is through the development of secondary porosities or permeabilities as a result of fracturing, jointing, shearing and deep weathering (Ademilua and Olorunfemi 2000; Oyinloye and Ademilua 2005). The fractures might have probably developed from transpressive forces during tectonic movement, pressure relief due to erosion of overburden rock, shrinking during cooling of the rock mass and the compression and tensional forces caused by regional tectonic stresses. Most of the groundwater in the study area occur in the weathered regolith/fractured parent rock zone (Fig. 3). Investigation during field work activity in the course of this study revealed that the varied distribution of the rock units in terms of mineralogy and the resulting differential weathering are responsible for the variable thickness of the in situ weathered regolith, which is also responsible for the differences in the groundwater potential at different locations. In SW-Nigeria, such variations in thickness usually in the range of 2–25 m are controlled by local and regional factors such as bedrock type, the degree of resistance to weathering, intensity and fracture density among others (Ako et al. 1990; Malomo et al. 1991; Tijani and Abimbola 2003).

Fig. 3
figure 3

Typical weathered soil profile developed on charnockitic bedrock from the study area

The overall hydrogeologic setting of the study area, as revealed by the hand-dug wells, is characterised primarily by the weathered regolith or saprolite units with varied thickness ranging from 4.3 to 15 m, located entirely within the saprolite unit with characteristically low yield of 6–246 l/min. Hence, the groundwater in the study area could be said to occur at shallow depths in localised and disconnected weathered regolith aquifers existing mainly as unconfined to semi-confined groundwater table conditions. The variable aquifers are tapped through shallow hand-dug wells that terminate mostly in the saprolite unit. This reported hydrogeologic profile is consistent with other research work reported elsewhere within the Basement Complex of SW-Nigeria (Tijani 1994; Tijani and Abimbola 2003).

Methodology

A total number of 73 water samples were collected mainly from domestic hand-dug wells within the study area for hydrochemical analysis, out of which 27 samples were selected for stable isotope analysis. The sampled locations are shown in Fig. 1. For the hydrochemical analysis, 73 samples were categorized into five classes based on rock formations, namely migmatite (30), porphyritic granite (11), fine to medium grained granite (17), quartzite (4) and charnockitic bedrock (11) settings while the isotope samples in the same order of bedrocks were eleven (11), four (4), seven (7), two (2) and three (3). All water samples were collected in 50-ml preconditioned polyethylene bottles. They were conditioned by rinsing initially with 10 % nitric acid and followed by several rinsing with distilled water to ensure that the sample bottles were free from any impurities.

At each sampling location, samples were collected into polyethylene bottles in triplicate. Two sets of samples were used for cation and anion analyses and 27 selected samples, the third set, for stable isotope analyses. Samples for cations analysis were acidified with concentrated nitric acid. Prior to this, field in situ parameters (pH, temperature and electrical conductivity) were measured using Multi-parameter TestrTm 35 series Meter. The samples were then preserved in a refrigerator prior to the laboratory analyses. The cation metals were analysed using the 36-element Inductively Coupling Plasma Oxygen Emission Spectrometer (ICP-OES) analytical method at ACTLABS Laboratory, Canada while the major anions were analysed at the Federal Ministry of Water Resources Laboratory, Alagbaka Akure, Nigeria using Ion Chromatography method. The analyses of the stable isotopes (18O and D) of the 27 water samples were carried out at Hydroisotop Laboratory, Schweitenkirchen in Germany using mass spectrometry method. The δ18O of groundwater samples from the study area were determined by equilibrating CO2(g) with water at 25 °C (Epstein and Mayeda 1953). The 18O/16O ratio of the CO2(g) was then measured on a Micromass MS903 mass spectrometer and the δ18O reported as per mil (‰) deviation from the Vienna Standard Mean Ocean Water (VSMOW) using the delta (δ)-scale with measurement accuracy of ±0.15. The mass ratio of deuterium in the groundwater samples of the study area was determined by reduction of the water molecules to gas following the Zn-reduction procedures outlined in Coleman et al. (1982). Water reacts with zinc at high temperature (400–450 °C) which proceeds by the following reactions:

$$ {\text{H}}_{ 2} {\text{O }} + {\text{ Zn}} \to {\text{H}}_{ 2} + {\text{ ZnO}} $$
(1)

Having gasified the hydrogen elements in the samples, consequently, the masses of hydrogen isotopes were determined by supplying the obtained hydrogen gas into a Micromass 602 mass spectrometer with the δD reported as  ‰ relative to V-SMOW.

Result of the hydrochemical data were subjected to graphical evaluation using Piper (1944) and Gibbs (1970) to define the various chemical evolution of groundwater in the study area. Furthermore, scatter plot of δ18O versus δD with respect to Global Meteoric Water Line (GMWL) was employed to determine the source of groundwater in the study area.

Results and discussion

Hydrochemical profiles

The results of the physico-chemical field measurements are presented in Table 1 while those of the hydrochemical analyses are shown in Table 2. The overall summary is presented in Table 3 along with WHO (2004) standard for drinking water quality. The pH of groundwater varied from 7.0 to 8.9 with a mean of 8.0. The temperature varied from 25.4 to 30.2 °C with a mean of 28.2 °C, which reflects the average air temperature during the sampling period (dry season) while the pH values revealed that the groundwater in the study area is neutral to slightly alkaline. Temperature and pH influence dissolution and precipitation of minerals in groundwater. Thus, the average pH value (av 8.0) is indicative of the presence of free CO2 and that the dissolved anions exist almost entirely in HCO3 ion form (Freeze and Cherry 1979). The electrical conductivity (EC) ranges between 109 and 939 μS/cm with a mean of 431.96 μS/cm while the total dissolved solid (TDS) values were generally less than 1,000 mg/L, ranging from 32.25 to 746.25 mg/L with an average value of 323.97 mg/L. The low TDSs are indications of low mineralized water with limited migratory history and water–rock interactions.

Table 1 Physico-chemical parameters analysed in water samples
Table 2 Hydrochemical parameters analyzed in water samples
Table 3 Descriptive summary of hydrochemical data in the study area

The total hardness (TH) varied from 31.28 to 295.10 mg/L with an average value of 128.49 mg/L. The guidelines for Canadian drinking water quality hardness (Health Canada 2012) divided hardness in equivalent concentration of CaCO3 into the following categories: Soft, <60 mg/L; Medium hard, 60 mg/L to <120 mg/L; Hard, 120 to <180 mg/L and Very hard, 180 mg/L or greater. Based on this classification, the groundwater system in the study area can be classified as soft to very hard water with about 70 % of the sample locations having TH > 60 mg/L which is the limit of soft water. The probable source of hardness in the groundwater of the study area could be from the weathering of calcium bearing minerals in most bedrocks of the area. However, hard water is mainly of aesthetic concern because of the unpleasant taste that a high concentration of calcium and other ions give to water. It also reduces the ability of soap to produce lather, and causes scale formation in pipes and on plumbing fixtures. The pH and TH have positive correlations with HCO3 with correlation coefficients of 0.50 and 0.94 respectively, while pH and TH are equally positively correlated with correlation coefficient of 0.56. These positive correlation values are clear indication of the interdependency and control of pH on the CO2 dissolution process in the groundwater system of the study area.

Calcium is the dominant cation with its values ranging from 1.5 to 106 mg/L (av 33.20 mg/L), while magnesium ranges from 0.4 to 39.2 mg/L (av 8.77 mg/L) in the study area. With respect to the alkali metals, the concentrations of Na ranged from 1.1 to 91.5 mg/L (av 25.21 mg/L) while that of K varied from 1.6 to 133 mg/L with an average of 25.13 mg/L. Bicarbonate is the dominant anion with concentration of 17.69–456.28 mg/L and average value of 164.85 mg/L. This is followed by Cl concentrations ranging from 2.8 to 61.5 mg/L with an average value of 32.9 mg/L. Sulphate has concentrations in the range of 1.5–76.10 mg/L with an average value of 21.06 mg/L, while NO3 concentrations ranged between 1.4 and 19 mg/L with an average value of 5.6 mg/L.

Furthermore, the cation concentrations was in the order of Ca2+ > K> Na> Mg2+ while that of the anions was HCO3 > Cl > SO42− > NO3. The higher concentration of Ca2+ and (K++Na+) in the groundwater of the study area may probably be due to the dissolution of plagioclase feldspars. HCO3 might have been generated in the soil zone en route to the groundwater zone as a result of decomposition of organic matter, which releases carbon dioxide that reacts with water in the soil zone. The reaction generates weak carbonic acid (H2CO3) that aids the breakdown of minerals in the rocks resulting in dissolution and the release of the ions into the groundwater which was responsible for its hydrochemical characteristics.

Chemical character and groundwater quality

All over the world, water chemistry is controlled by some few salient factors which include precipitation, mineral weathering/dissolution as well as evaporation–recrystallisation processes (Gibbs 1970). In addition, the concentration of dissolved ions in groundwater samples are generally governed by lithology, nature of geochemical reactions and solubility of interacting rocks (Nosrat and Asghar 2010). In many cases, water within the various phases of the hydrological cycle is affected differently by these controlling processes. Apart from this, other processes such as cation–anion exchange and mixing processes may also exert some considerable influence on the groundwater chemistry (Tijani and Abimbola 2003). Thus, water–rock interaction contributes significantly to evolved chemical composition of the groundwater system in the study area in addition to contribution from precipitation. Results of hydrochemical analysis in the study area (Tables 1, 2) clearly reflect varied chemical composition with respect to the measured parameters which is a reflection of the varied lithologic units in the area as well as differential weathering/dissolution in such bedrocks.

Strong correlation exists among the major cations and TDS with correlation coefficients >0.60 (Table 4) which is a clear indication of the contribution of these ionic components to the overall mineralisation (TDS). Similarly, with the exception of NO3 versus TDS (r = −0.17), all the major anions also have positive correlation coefficients (r > 0.70) (Table 4). This observed chemical character suggested that NO3 concentration in the groundwater system was from different source possibly from anthropogenic contamination (household contaminants and waste pits). Hence, most ionic input into the groundwater system of the study area was mostly from geogenic source. However, few cross plots of the groundwater ionic parameters in the study area as indicated in Fig. 4(a–f) signify distinctly two types of water. The first type of water has low TDS while the second type has fairly high TDS values. Hudson and and Golding (1997) pointed out that the concentrations of sodium, magnesium, calcium and silica in groundwater system of any area are controlled by mineral weathering of the bedrocks. This was clearly demonstrated in the study that the percolated CO2-charged rain water turned to weak carbonic acid initiating water–rock interactions during which weathered minerals released cations (Ca, Mg, Na and K) that reacted with the weak carbonic acid to produce bicarbonate waters. However, the low TDS in the groundwater system of the study area was a reflection of the general characteristics of groundwater in basement terrain. The groundwater in such basement terrain has limited migratory history and its composition was controlled by recharge water from rainfall and mineral dissolution from slow weathering of the minerals constituting the bedrocks (Tijani 1994).

Table 4 Correlation matrix of major ions in groundwater of the study area
Fig. 4
figure 4

Cross plots of major ionic parameters in the analysed groundwater samples from the study area

Furthermore, a plot of TDS against Na+/(Na++Ca2+) ratio referred to as Gibbs plot (Fig. 5) assessed the functional source of dissolved ions in waters as evaporation dominance, rock dominance (dilution and weathering) and precipitation dominance which controls the water chemistry (Gibbs 1970). The chemical data of the groundwaters from the study area plotted mainly around the weathering zone and therefore indicate that the chemical compositions of these waters were controlled primarily by rock dominance (weathering reactions). Therefore, the waters of the study area have their chemistry controlled by the weathered materials derived from the underlying bedrocks.

Fig. 5
figure 5

Controlling mechanisms for groundwater quality in the study area (after Gibbs 1970)

Regarding water quality for potability and usability, all the various chemical components of the groundwater system in the study area fell within the WHO (2004) acceptable limits for drinking water standard. The irrigation quality assessment was carried out using EC (μS/cm), SAR (Richards 1954), Kelly ratio (Kelly 1963) and magnesium adsorption ratio (MAR) (Paliwal 1972). Detailed formula employed in estimating irrigation quality parameters is highlighted below with ionic concentrations in meq/L:

$$ {\text{SAR}}\, = \,\frac{{{\text{Na}}^{ + } }}{{\sqrt {\frac{{({\text{Ca}}^{2 + } ) + ({\text{Mg}}^{2 + } )}}{2}} }} $$
(2)
$$ {\text{KR}}\,{ = }\,\frac{\text{Na}}{{{\text{Ca}} + {\text{Mg}}}} $$
(3)
$$ {\text{MAR}}\,{ = }\,\frac{{{\text{Mg}}^{2 + } }}{{({\text{Ca}}^{2 + } + {\text{Mg}}^{2 + } )}} $$
(4)

Fifty-nine (59) samples out of total seventy-three (73) samples, representing about (81 %) of the water samples, have EC < 700 μS/cm except for one sample with EC > 1,000 μS/cm viz., PG11, 1,041 μS/cm. Normally, irrigation water with an EC < 700 μS/cm causes little or no threat to most crops while EC > 3,000 μS/cm may limit their growth (Tijani 1994). Esimated Kelly ratio <1 for all water samples implies groundwater suitability for irrigation (Raihan and Alam 2008), MAR for the water samples ranges from 3.41 to 46.17 which is within permissible limit of 50 (Paliwal 1972). Mandel and Shiftan (1991) stated that water with SAR < 10 can be applied on all soil for irrigation usability and the groundwater system in the study area fell into this category. Thus, based on these evaluations, groundwater in the study area was suitable for agricultural purposes.

Hydrochemical facies

The Piper trilinear diagram (1944) is used to infer hydrogeochemical facies. These plots include two triangles, one for plotting cations and the other for plotting anions. The cation and anion fields are combined to show a single point in a diamond-shaped field, from which inference is drawn on the basis of hydrogeochemical facies concept (Back and Hanshaw 1965). The trilinear diagram is useful in bringing out chemical relationships among groundwater samples in more definite terms compared to other convectional plotting methods. The concept of hydrochemical facies was developed to understand and identify the water composition in different classes. Facies are recognizable parts of different characters belonging to any genetically related system. Hydrochemical facies are distinct zones that possess cation and anion concentration categories (Rajendra Prasad et al. 2009).

Plots of the hydrochemical parameters of the water samples on the Piper (1944) trilinear diagram (Fig. 6) classified the groundwater system in the study area into Ca(Mg)HCO3 (dominant) and Na(K)HCO3 (minor) waters. Ca(Mg)HCO3 water facies represent about 89 % of total water samples from the study area. The two main water facies cut across the entire bedrocks of the area except in the Na(K)HCO3 water where only water samples from migmatite, quartzite and fine-medium grained granite bedrocks were represented. The mineralogical composition of the various rocks especially (migmatite and granite) in the study area could easily weather to produce CaMgHCO3 water type. However, since the study area is located in the rain forest zone where total annual rainfall is about 1,500 mm, the granites weathered to form thick layers of clay deposit and as such cation–anion exchange reactions were therefore suspected to contribute to the groundwater evolution of the study area.

Fig. 6
figure 6

Water type classification using the Piper trilinear diagram (Back and Hanshaw 1965). A Calcium type,B sodium or potassium type,C magnesium type,D no dominant type, E bicarbonate type,F chloride type,G sulfate type,H no dominant I = Ca + Mg type, J = HCO3 + CO3 type, K = Na + K type, L = Cl + SO4 type M = Ca(Mg) HCO3 type, N = Ca(Mg) Cl SO4 type, O = Na(K) HCO3 type, P = Na(K) Cl(SO4) type

Stable isotopes profiles (δ18O and δD ‰)

Isotope techniques can be used to solve problems such as identification of the origin of groundwater, determination of its age, flow velocity and direction, interrelations between surface waters and groundwaters, possible connections between different aquifers, local porosity, transmissivity and dispersivity of an aquifer (Faure 1998). Stable isotopes of 18O and Deuterium are ideal tracers of water sources and movement because they are constituents of the water molecules (Akiti 1986). In addition, 18O and Deuterium undergo changes in the various phases of hydrologic environments, resulting in isotopic fractionation that gives isotopic “fingerprints” (Das et al. 1998).

The isotopic composition of oxygen and hydrogen is reported in terms of differences of 18O/16O and D/H (2H/1H) ratios relative to a standard called Standard Mean Ocean Water (SMOW) (Craig 1961).

$$ \delta^{ 1 8} {\text{O}}\, = \,\frac{{\left( {^{18} {\text{O}}/^{16} {\text{O}} -^{18} {\text{O}}/^{16} {\text{O}}_{\text{SMOW}} } \right)\, \times \,10^{3} }}{{^{18} {\text{O/}}^{16} {\text{O}}_{\text{SMOW}} }} $$
(5)
$$ \delta \,{\text{D}}\,{ = }\,\frac{{\left( {^{2} {\text{H/}}^{1} {\text{H}} -^{2} {\text{H}}/^{1} {\text{H}}_{\text{SMOW}} } \right)\, \times \,10^{3} }}{{^{2} {\text{H}}/^{1} {\text{H}}_{\text{SMOW}} }} $$
(6)

Consequently, positive values of δ18O and δD indicate enrichment of the sample relative to δ18O and δD as compared to SMOW, whereas negative values imply depletion of these samples relative to the standard. On the basis of large numbers of meteoritic water collected at different latitudes, it has been shown that δ18O and δD values of meteoritic water are linearly related and can be represented by the equation δD = 8δ18O + 10 (Craig 1961). Data plotted on this suggest that such water is derived from the meteoritic origin. Large deviations also imply either enrichment or depletion.

For this study, 27 locations were sampled and analysed for the stable isotopes (18O and Deuterium). The results and statistical summary of the stable isotopes compositions in δ-notation are presented in Table 4 with all values reported in per mil (‰) in reference to V-SMOW and measurement accuracy of ±0.15 for both δD and δ18O. The relationship between δD and δ18O in the classical δD versus δ18O plot (Fig. 7) alongside with the global meteoric water line (GMWL) defined by Craig (1961) usually reveals indications about origin and recharge source(s) of groundwater system as well as possible evaporation effect (Murherjee et al. 2007).

Fig. 7
figure 7

Stable isotopes of Ekiti groundwaters

The measured δD values of the analysed water samples range from −22.3 to −10.2 ‰ with an average value of −15 while that of δ18O range from −4.2 to −2.1 ‰ with an average value of −3.1 (Table 5). With the exception of location C2 on charnockite bedrock that deviate slightly from GMWL, all sampled water do not only show depletion compared to that of the sea water (δD = 0  ‰; δ18O = 0  ‰) but also plotted along the GMWL as defined by the equation δD = 8δ18O + 10 (Fig. 7). The general implication of this is that the origin of groundwater in the study area is related to meteoric source, i.e. recent precipitation water recharging the associated shallow basement aquifer, with little or no imprint of kinetic evaporation. This result is further supported by the fact that deuterium excess (d-excess) values are generally in the range of 10 ± 2, i.e. close to 10 which support the observation of recent precipitation with little or no imprint of kinetic evaporation (Table 5). Similar isotopic signature had been reported for shallow groundwater in other parts of the Basement Complex of southwestern Nigeria (Loehnert 1988; Tijani and Abimbola 2003).

Table 5 Result of stable isotopes in the study area

The observed isotopic signature was in agreement with the early conclusion from hydrochemical data, i.e. evolution of the groundwater was related to infiltration of CO2-charged precipitation with limited migratory history characterised by low TDS. The shift of the sample on location C2 slightly away from GMWL with d-excess of 6.8 can be the result of limited degree of kinetic evaporation before or during infiltration process. Furthermore, as shown in Fig. 7, the groundwater samples appeared to group in a narrow range signifying a well mixed system with relatively constant isotopic composition. This observation indicates that the meteoric water recharging the groundwater system in the study area is relatively homogenous with evaporation playing insignificant role on the infiltrated water.

Groundwater undergoing evaporation will have a positive correlation between δ18O and electrical conductivity (Gibrilla et al. 2010). In Fig. 8a showing the plot of δ18O versus EC, it was observed that a positive correlation of r = 0.54 supports earlier observations that few of the groundwater samples plotted slightly away from the GMWL, which are indicative of possible mild evaporation prior to recharge or kinetic interaction with the subsurface medium. With EC < 300 μS/cm, δ18O (‰) changed minimally as shown by samples in circle B. Samples in circle A illustrate evaporative effects with mineralization accompany enrichment of δ18O (‰) (Fig. 8a). Nitrate content in the groundwater at any point or area in the saturated aquifer is controlled by the water recharging through the unsaturated zone and the lateral flow from surrounding areas. Therefore, depending on the recharge/infiltration characteristics and the period of stay of water on the land surface (and thereby subjected to evaporation) at different locations, some kinds of correlation are expected between δ18O and NO3 contents of groundwater (Akiti 1986). It is observed from Fig. 8b that correlation coefficient of 0.25 is obtained from the plot of δ18O versus NO3. In circle I of Fig. 8b, enriched δ18O signifies evaporation while that of circle II indicates low nitrate groundwater associated with depleted δ18O water implying anaerobic conditions at deeper levels in the soil profile (Datta et al. 1997). Some of the wells were recharged by significant quantities of evaporated (18O enriched) rainwater infiltrating along with NO3 salts from the soil to the groundwater system. In circle III, fairly constant δ18O with relatively high NO3 values support argument of surface source of NO3. Finally, assessment of the deuterium excess (d-excess) reflects the conditions that lead to kinetic isotope fractionation between water and vapour during primary evaporation in the oceans (Dansgaard 1964). It shows the extent of deviation of a given sample from the Global meteoric water line. The calculated deuterium excess of the groundwater was found ranging from 6.80 to 11.00 ‰. As the δ18O increases (more enriched), the deuterium excess in most of the samples decreases gradually (Fig. 8c and Table 5) signifying the slight evaporative effects on the groundwater system. However, the fact that the groundwater samples plot in close range to GMWL imply a modern recharge source to the groundwater system of the study area.

Fig. 8
figure 8

Cross-plots of δ18O (‰) versus EC(μS/cm), NO3(mg/L) and d-excess

Conclusion

Assessment and evaluation of the hydrochemical data from the study area have thrown light on the hydrochemistry, quality and stable isotope characterisation of shallow groundwater system from the weathered basement aquifer of Ekiti area, southwestern Nigeria. The physical and chemical parameters of all sampled water fell within the WHO (2004) guideline values for drinking water and revealed that the groundwater was potable. Further, evaluation of the hydrochemical results showed that the groundwater system in the study area was fresh (32.25 ≤ TDS ≤ 746.25 mg/L), low mineralised with limited water–rock interactions. The groundwater system was soft to moderately hard (31.28 ≤ TH ≤ 295.10 mg/L) and slightly alkaline (7.0 ≤ pH ≤ 8.9). Quality evaluation in terms of agricultural uses revealed that the estimated irrigation indices using SAR, MAR and Kelly’s ratio indicated that the groundwater was suitable for irrigation purposes.

The principal source of the dissolved chemical species in the groundwater system of the study area was from mineral weathering/dissolution process facilitated by the action of infiltrating rain water rich in atmospheric and/or biogenic CO2. This mineral weathering/dissolution was the main process responsible for the hydrochemical evolution of the groundwater system in the study area. This was further supported with the fact that most water samples plotted in the rock dominance part of the Gibb’s diagram. The principal hydrochemical facies in the study area were Ca-(Mg)-HCO3 and Na-(K)-HCO3 waters. Ca-(Mg)-HCO3 water represents geogenic characters of the shallow groundwater system in the study area which was the result of weathering/dissolution of aluminosilicate minerals. This is also consistent with the results of stable isotope analysis which revealed meteoric source of groundwater recharge (i.e. recent precipitation with insignificant degree of kinetic evaporation recharging the associated shallow basement aquifer).

Finally, the fact that the aquifers in the study area are shallow and localised and that recharge is primarily from infiltrating precipitation, warrants the need for proper design and construction of sewage pits to forestall anthropogenic pollution due to seepage/leakage that may arise from possible hydraulic connection with the shallow groundwater system in the study area.