<figure><img src="image_1.png" /></figure>

HAL
open science

# Identifying the origin of springs in weathered-fractured crystalline aquifers using a hydrogeophysical approach

Kouassi Jean-Michel Kouassi, Patrick Lachassagne, Oi Mangoua Jules

Mangoua, Abé Parfait Sombo, Brou Dibi

► To cite this version:

Kouassi Jean-Michel Kouassi, Patrick Lachassagne, Oi Mangoua Jules Mangoua, Abé Parfait Sombo, Brou Dibi. Identifying the origin of springs in weathered-fractured crystalline aquifers using a hydrogeophysical approach. Scientific Reports, 2024, 14 (1), pp.12977. 10.1038/s41598-024-63748-8. hal-04870762

HAL Id: hal-04870762

https://hal.science/hal-04870762v1

Submitted on 8 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Identifying the origin of springs in weathered-fractured crystalline aquifers using a hydrogeophysical approach

Kouassi Jean-Michel Kouassi<sup>1,2,✉</sup>, Patrick Lachassagne<sup>2</sup>, Oi Mangoua Jules Mangoua<sup>1</sup>, Abé Parfait Sombo<sup>1</sup> & Brou Dibi<sup>1</sup>

Over the last few decades, important advances have been made in the development of relevant hydrogeological conceptual models for crystalline aquifers, and notably for weathered-fractured crystalline aquifers. Paradoxically and contrary to other types of aquifers, these researches never aimed at characterizing springs, the places were groundwater naturally outflows from such aquifers. With such an objective, our methodological approach consisted first of a lithological and hydrogeological description of the aquifer system based on borehole data and outcrops in a representative weathered-fractured crystalline aquifer (Daloa, Ivory Coast). Next, electrical resistivity tomography (ERT) has been used (after validating the appropriate inversion method) to provide the imagery of the weathering profile both below the plateaus and in the valleys where the springs outflow. Piezometric and river discharge data were also processed notably to determine the direction of groundwater flow. Results demonstrate unambiguously that the isalterites aquifer supplies the springs, and that the underlying fractured layer is not directly implied in this supply. ERT combined with borehole and field lithological data also shows that the lateritic formations (alloterites) present near surface below the plateaus, as well as the upper part of the isalterites, were eroded in the valleys, but not deep enough to let the fractured layer outcrop. This conceptual model for springs not only provides a basis for characterizing such complex aquifers, but also provides technical guidance for spring catchment and groundwater protection in these crystalline areas.

Water is essential for life, both within ecosystems and for humanity. Yet millions of people around the world struggle daily to find drinking water to meet their basic needs<fnref n="1" />,². Groundwater is increasingly used worldwide due to the impact of climate change and anthropogenic activities on surface water resources, both in terms of quantity and quality. Groundwater covers around 43% of irrigation water needs³, and also makes an important contribute to hydration, sanitation and hygiene in many parts of the world². Although groundwater is increas- ingly in demand throughout the world, its exploitation and management are often difficult and costly due to the complexity of the aquifer systems that contain it, notably those that are fractured⁴⁻⁸. The African continent is no exception. Moreover, over a large part of its surface, the only available resource of groundwater comes from aquifers from granitic and metamorphic rocks known as "crystalline"⁴,⁵. Worldwide, these crystalline aquifers make up more than 20% of the world's land surface⁵,⁹,¹⁰. They have a fracture permeability⁵ and, as a result, their structure and functioning are less well understood than those of other types of aquifers, notably porous aquifers, making it difficult and costly to exploit the groundwater they contain⁵.

In West Africa, groundwater accounts for a significant and steadily increasing proportion of the rural popula-
tion's water supply, particularly in the crystalline context. Indeed, the groundwater contained in these crystal-
line aquifers is geographically well distributed<sup>11-15</sup> and offers an alternative to surface water resources, which
are subject to contamination and are unavailable all or part of the year, particularly in semi-arid regions. These
groundwaters resources are generally accessed by drilling. However, the complexity of crystalline aquifers means
that between 25 and 60% of boreholes drilled are "dry"<sup>16-18</sup>. However, these West African crystalline rock regions
also contain numerous springs<sup>19,20</sup>, which could also be exploited to supply the population with water. Springs
are natural groundwater outflows that appear locally or diffusely on the surface of the ground. They play an

<footnote n="1">¹Laboratory of Environmental Sciences and Technologies, Univ. Jean Lorougnon Guédé, Daloa, Ivory Coast. ²HSM, Univ. Montpellier, CNRS, IRD, IMT Mines Alès, Montpellier, France. ✉email: kouassimichel48@gmail.com</footnote>
important role in socio-economic development in many parts of the world, where they are tapped for a variety of uses. Paradoxically, in crystalline aquifers, the processes and hydrogeological context that explains why groundwater outflows at springs is less well known than that of boreholes⁵. It is therefore important to improve the hydrogeological conceptual model of crystalline aquifers in areas where springs are present, in order to characterize the type of springs, their context of emergence, their vulnerability to contamination and to verify their seasonal sustainability²¹,²².

Over the last two decades⁵,²³, the hydrogeological conceptual model of crystalline aquifers has evolved considerably, from concepts invoking fractures of tectonic origin to a concept of fractures linked to weathering processes. Methods for characterizing the structure and functioning of these hydrosystems have evolved accordingly⁴,⁵,¹⁶,²¹,²³⁻²⁶. However, these methods do not provide relevant information on the origin and properties of springs in this type of geological environment. In this research, we therefore explore the use of electrical geophysical methods to develop the conceptual model of springs in a crystalline rock context, using an integrated approach in which geological, hydrogeological and geophysical data are jointly interpreted.

Over the last few decades, numerous advances have been made in the development of representative conceptual models of crystalline aquifers based on geophysical methods²⁷⁻³⁴. However, these researches never addressed the issue of spring characterization. Electrical resistivity tomography (ERT) methods are widely used in this type of context, as they provide 2D information, or even quasi 3D in the case of serial profiles, and as there is a good correlation between geophysical facies and those of the weathering profile³⁵⁻³⁸. The complexity of crystalline environments gives rise to biases that severely limit the validity of 1D methods²⁸,³⁷,³⁹,⁴⁰ which, for technical and low-cost reasons, are unfortunately still often used.

Building on feedback from previous research, this study aims to assess the effectiveness of ERT to characterize the context of groundwater springs in crystalline rock environments. Specifically, we aim to: (i) assess the impact of geophysical data inversion methods on the quality of hydrogeological models in a crystalline rock context, and (ii) leverage a large geological, hydrogeological and geophysical data set to improve the conceptualization of crystalline aquifers. To meet these objectives, we first described the lithology and hydrogeology of the aquifer system using borehole and outcrop data. The results of geophysical inversions were then calibrated on these geological data to provide a 2D conceptualization of the weathering profile. Piezometric data were also processed to determine the groundwater flow direction in the aquifer.

The research is being carried out on the aquifers from the catchment area of the Tétégbeu River (Fig. 1), a tributary of the Lobo River (Center-West, Ivory Coast). The Tétégbeu River starts flowing in the municipality of Daloa and drains an area of around 140 km². The climate belongs to the attenuated transitional equatorial regime, with average annual rainfall and temperature of 1200 mm and 26 °C respectively (see Appendix 2 in the Complementary Materials). The topography of the basin is relatively homogeneous (Fig. 1), with predominantly flat, gently sloping plateaus incised by the hydrographic network. The subsoil in the research area is composed of Middle Precambrian formations, dominated by granites and granodiorites⁴¹, with granodiorites accounting for 85% of the surface area (see Appendix 1 in the Complementary Materials). These plutonic rocks were deeply weathered during the Mesozoic to Quaternary eras, and are therefore covered by weathering profiles several tens of meters thick.

# Results

## Description of the weathering profile at borehole scale

Analysis of the lithological cross-sections of the boreholes (observations made during drilling, notably on cuttings) shows the succession of three very distinct horizons within the weathering profile (Fig. 2), which follow a classic vertical structuring⁵ with, from top to bottom: alloterites (10–40 m thick on complete, non-eroded profiles), isalterites (20–30 m thick on complete, non-eroded profiles), which lie on the fractured layer within which water inflows are observed during drilling. From bottom to top, below the soil horizon, the alloterites comprise the lateritic iron duricrust, locally dismantled into lateritic gravels, and lateritic spotted clays (Fig. 2). The alloterites, with their very low permeability, form an aquitard on this site and in general⁵. The isalterites, on the other hand, constitute a "medium" permeability aquifer (about $10^{-6}$ m/s)⁵, with an effective porosity that can exceed 5%⁵ (from our laboratory porosity analyses carried out using the imbibition method on intact soil samples). This aquifer is generally, and also in Daloa, tapped by traditional large-diameter shallow wells, usually located in valleys where the piezometric level is close to the ground surface. The fractured basement (the fractured layer) underlies the isalterites. In the Daloa region, statistical analysis of linear discharges and water inflow depths⁴,¹⁶ from 316 deep boreholes reveals that the first 35–40 m of the fractured layer constitute its most densely fractured, most productive part (See Appendix 7 in the Complementary Materials).

## Calibration and geoelectric signature of the formations constituting the weathering profile

The ERT cross sections were inverted with two different inversion methods (“Robust” and “Standard”) which are among the most widely used in geophysics (Fig. 3a, b). All ERT cross sections and inversions show a similar succession of stratiform layers of varying resistivity, with from top to bottom: a resistive layer ($R1 = 500–1600 Ω.m$) at the surface, up to 20 m thick, followed by a conductive layer ($C = 85–500 Ω.m$) 30–40 m thick, relatively continuous and homogeneous, and, below, very resistant formations ($R2 > 1600 Ω.m$).

Regardless of the inversion method used, comparison of geoelectric profile results with borehole lithology (e.g. borehole F4 in Fig. 4) shows that the resistant geoelectric layer observed at surface corresponds to the alloterites, characterized notably by lateritic formations. This layer is geologically and electrically heterogeneous, and predominantly resistant. Below this resistant layer, a more homogeneous conductive layer appears. This conductive layer corresponds to the lower, water-saturated, part of the saprolite, largely corresponding to the sandy saprolite (isalterite) and to the most superficial part (15–20 m) of the underlying fractured layer. Below this
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Location and topographic map of the study area. Location of the sites selected for ERT surveys. SP_F: site of ERT profiles SP_F4 and SP_F15 (the lithology of boreholes F4 and F15 was used to calibrate the SP_F4 and SP_F15 ERT profiles), STU: site of ERT profiles STU1 and STU2, SL: site of ERT profile SL, SA1: site of ERT profile SA1, SA2: site of ERT profiles SA2 and SZ: site of ERT profile SZ. The electrical cross sections for the SL, SA1, SA2 and SZ profiles are shown in Appendix 4 and 5 in the Complementary Materials. This graph was produced using ArcGIS 10.2.</figcaption></figure>

conductive layer, a highly resistant layer follows the conductive layer, via a resistivity gradient that is tighter with the Robust inversion method than with the Standard one. This latter layer corresponds to the middle and lower part of the fractured layer, as the depth of investigation of the ERT profiles does not allow to reach the base of the fractured layer and the underlying unweathered high resistivity rock. The saprolite-fractured layer interface lies at the heart of conductive layer C, frequently at the base of its most conductive part. This "isoline" (base of the most conductive part of conductive layer C) was used to extrapolate this interface; however, the precision of its positioning remains low (accuracy of the order of 5–10 m). The resistivity gradient observed within the fractured layer most likely corresponds to a fracture electrical connectivity threshold (in the sense of the percolation theory). Permeable fractures (identified as water occurrences during drilling) are not visible on the ERT.

Results from both inversion methods show the same sequence of resistive (R1)—conductive (C)—resistive (R2) layers (with R1-C-R2 from top to bottom) (Figs. 3 and 4). The depth and, above all, the progressiveness of the conductive-resistive boundaries (R1-C and C-R2) vary from one inversion method to the other. The upper boundary (R1-C), marked by a relatively abrupt transition, is fairly well defined on the ERT by both inversion methods, although the gradient is slightly stronger with the Standard inversion method than with the Robust one. On the basis of borehole data alone, it is difficult to determine whether this R1-C boundary corresponds to the contact between alloterite and isalterite or to the saturation limit within the weathering profile (piezometric level). However, as this boundary has abrupt variations and locally steep slopes, it cannot correspond to the piezometric level, whose spatial variations are much more regular. The piezometric level is thus not imaged by ERT. The lower limit (C-R2) is not defined in the same way by the two inversion methods. The resistivity gradient is stronger with the Standard method than with the Robust one. On this basis, the Robust inversion method was chosen for the interpretation of all the ERT profiles in this research. Nevertheless, inversions using the Standard method were also taken into consideration to avoid any artifacts associated with inversion.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Lithological profiles illustrating the weathering profile at the study site, based on borehole log observations. This graph was produced using Adobe Illustrator V25.1.</figcaption></figure>

## Geoelectrical signature of the weathering profile at springs

Given the limited geological information available at the springs, the interpretation of the ERT profiles was mostly based on the calibration of geophysical investigations carried out at the boreholes (see previous section) and on observations made at outcrops, notably at sand quarries located in the saprolite in valleys or close to the springs, or even at the springs themselves.

### Interpreting ERT profiles

The STU1 ERT is located 60 m upstream of the Tazibouo_U spring. It was sited perpendicular to the topographic slope and therefore perpendicular to the assumed groundwater flow axis (Fig. 5a). On this profile, the highest topographic heights are found at the highest abscissas (right-hand side of Fig. 5a). The second ERT, STU2, is located perpendicular to the STU1 profile and passes over the spring (Fig. 5b). At their point of intersection (abscissa 125 m on ERT STU1 and abscissa 120 m on ERT STU2), the two profiles show a very similar geoelectric signature and also similar depths for the contacts between the different geoelectric layers. This demonstrates the accuracy of the acquisitions and inversions.

The geoelectric layers observed on the STU1 and STU2 ERT are, from top to bottom, as follows: - a first conductive layer (resistivity 80-500 Ω.m) with a thickness fluctuating between 30 and 55 m is observed near the surface. Observations made during the building of the catchwork of this spring revealed the presence of granitic arenas (isalterites) at a depth of 3 m, topped by a tin sandy and sandy-clayey layer (soil), confirming the geological calibration of ERT carried out on the basis of drilling data (cf. previous section). The laterites and lateritic duricrust observed on the top of the plateau (particularly in borehole F4) would thus have been partially or completely eroded in this valley. Based on the calibration performed on the boreholes, the depth of the transition between the isalterites and the fractured layer would thus be located between 12 and 40 m below the topographic surface along the ERT, at elevations between 275 and 252 m asl. This interface is relatively flat on the STU2 profile. It shows variations of around 20 m in amplitude on the STU1 profile; - the underlying layer correspond to the resistivity gradient observed on the profiles described in the previous section, which images
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Interpretation of SP_F4 ERT profile with two inversion methods. (a): ERT profile inversed with the Robust inversion method; (b): location map of SP_F4, STU1 and STU ERT profiles (see the location of this area on Fig. 1); (c): ERT profile inversed with the Standard inversion method. Inversions were performed using Res2Div V-3.57.36.</figcaption></figure>

the upper and middle parts of the fractured layer, whose resistivity is between 500 and 1600 Ω.m in its upper part and > 1600 Ω.m beyond that.

The spring thus emerges from the isalterites, which constitutes the aquifer that feeds it (Fig. 5d). The transition between the isalterites and the fractured layer lies several tens meters below the spring.

Electrical imaging of the SL, SA1, SA2 and SZ profiles (see Appendix 4 and 5 in the Complementary Materials), which are also located on or near springs, show the same geoelectrical succession as the STU ERT, with a conductive medium near the surface, with resistivities between 80 and 500 Ω.m, corresponding to isalterites. Auger sampling at several of these sites confirmed the presence of coarse sand at 1 m depth (below the pedological layer). Beneath this conductive layer, resistivity increases as on the STU profiles and on profiles located at a distance from the springs. The geoelectrical profiles observed on the plateaus and in the valleys, and therefore their different weathering profiles, are therefore similar. The only difference is that weathering profiles in the valleys were partly eroded. Everywhere on this site, springs emerge in the heart of the isalterites (Fig. 6). Unlike springs emerging from sandstones, where springs outflow from conduits located where fractures have been enlarged by weathering and/or erosion due to the groundwater flow⁶,⁷, at the outlet of the springs in the saprolite (and elsewhere in the saprolite), there are no fractures. Consequently, there are no conduits either.

In detail, we also note, on almost all of these ERT, the presence of intermediate resistivity structures (resistivity: 500–1600 Ω.m) within the subsurface conductive layer (C). These structures are associated with local lower depth or deepening of the base of the conductive layer (C) at the expense of the underlying resistant layer (R2), which corresponds to the median part of the fractured layer (Fig. 6). These vertical variations indicate the presence of heterogeneities within the weathering profile linked to differential weathering of the parent rock: quartz veins less weatherable than their host rock (some were observed in the field), pegmatite or green rocks veins also observed on some drill logs, boulders, etc.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Vertical distribution of resistivity as a function of the inversion method.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Interpretation of ERT profiles at the Tazibouo_U spring (Robust inversion). (a): ERT STU1; (b): Location of ERT STU1, STU2 and SP_F4 (see the location of this area on Fig. 1); (c): ERT profile STU2; (d): Conceptual model of the Tazibouo_U spring.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Lateral evolution of geoelectric signatures from the plateau to the valleys. The SL and SA1 ERT are presented in Appendix 4 of Complementary Materials.</figcaption></figure>

## Typology of springs in a weathered crystalline rock context

Figure 7 shows that the elevation of springs, of the piezometric level in boreholes and the elevation of the perennial parts of watercourses observed during low-water periods are perfectly consistent each with each other (the same applies to high-water periods, but only when there is no more runoff in the river due to a recent rainfall event). As groundwater supplies the surface water flow, by taking all these data into account to draw up the piezometric map, we were able to obtain a faithful representation of reality (Fig. 7). This approach proved much more accurate than the one using only the piezometric elevation of the 14 boreholes surveyed. The piezometric map (Fig. 7) shows that the groundwater flow follows the topography, as observed in most crystalline aquifers⁵. The perennial parts of the Tétégbeu River drain the aquifer, that is unconfined. The river therefore owes its low-water flow to the groundwater outflow. Springs also drain the aquifer. The piezometric gradient is of the order of 0.5%, which corresponds to a fairly high value, consistent with the modest permeability of the isalterite aquifer.

The combined analysis of information from borehole geological cross-sections, ERT, topographic and piezometric maps provides a coherent hydrogeological conceptual model (Fig. 8) of the study area, enabling us to understand the context of spring emergence in a weathered crystalline rock. With regard to the geological structure of the weathering profile, it is shown that the lateritic formations observed on the plateaus (relics of the weathering paleosurface) have been eroded into the valleys. In the valleys, only isalterites outcrop and, possibly, recent iron crusts of late secondary weathering of this eroded weathering profile. The springs emerge from and are fed by the isalterite aquifer. The fractured layer, usually tapped by boreholes, is located some tens of meters deeper. Springs are tapped in the upstream part of the area, where the piezometric level begins to tangent the topographic surface. Downstream wetlands are also zones of diffuse emergence of the isalterite groundwater.

## Discussion

The combination of results from ERT on the plateaus and in spring zones, with borehole lithological cross-sections and outcrop observations at a weathered-crystalline aquifer site has enabled us to characterize the architecture of the aquifer that supplies the springs. ERT profiles revealed a weathering profile several tens of meters thick, composed of three main geoelectric layers below the plateaus and two geoelectric layers in the spring emergence zones (valleys). These results are in line with those of previous studies carried out in a similar geological environment, but which were not aimed at understanding the context of spring emergence²⁶,³²,⁴². From top to bottom, we observe: a resistant layer (500–1600 Ω.m) including locally relatively more conductive zones, then a highly conductive layer (80–500 Ω.m) and a highly resistant layer (resistivity > 1600 Ω.m). The resistant layer observed near the surface on the plateaus doesn't exist in the valleys, notably near the springs,
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Low-water piezometric map of the Tétégbeu river watershed (Dec. 2021). The Grass Green dots and associated numerical values represent the location and elevation (in m amsl) of groundwater at the springs, respectively. The Ruby Red points and their associated numerical values represent the location and piezometric level of boreholes, respectively. Blue dots and associated numerical values represent the location and water level in leveled points in the river, respectively. Refer to the text for the methodology used for the contouring. This graph was produced using Surfer 19 and was manually improved.</figcaption></figure>

where isalterite outcrops are also frequently observed. This clearly shows that the laterites and iron crust observed on the plateaus, as well as a part or the totality of the underlying spotted clays have been eroded in the valleys. The slightly more resistant subsurface soils observed are relicts of these iron crust or, more likely, late secondary ones⁴³,⁴⁴. The zones of high resistivity (resistivity > 1600 Ω.m) that stand out within the resistant layer at the surface on the plateaus would correspond to such duricrust observed at the surface at many points, where it is not covered by a gravelly or lateritic layer (duricrust in the process of being dismantled). The wide variations in resistivity of this resistant layer reflect the highly heterogeneous nature of the lateritic formations. This shows that, in these areas of weathered crystalline rock, the alloterite layer comprising these lateritic formations is highly heterogeneous. Based on direct field observations, zones of resistivity above 1600 Ω.m can be attributed to areas where the duricrust is highly indurated and intact, or to the presence of quartz veins within these lateritic formations (See Appendix 6 in the Complementary Materials). Zones of lower resistivity characterize areas where the duricrust has been fractured and cracked by dismantling, forming laterites, or lateritic gravels. Such
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Generalized conceptual model of springs in a weathered-fractured crystalline rock context. This graph was produced using Adobe Illustrator V25.1.</figcaption></figure>

lower resistivity zones also correspond to lateritic clays²⁹,⁴⁵. This resistant layer and its heterogeneities observed in the subsurface on ERT sited on plateaus do not appear on the ERT sited in the valleys, particularly at or near spring emergence zones, where the underlying conductive layer characterizes isalterites near the topographic surface. This unambiguously shows that the springs emerge from the isalterites, at 20–30 m above its contact with the underlying fractured layer. The saprolite thus constitutes the aquifer supplying the springs, with the upstream emergence zone corresponding to the intersection of the piezometry with the topographic surface. The statistical analysis of linear flow rates and the vertical distribution of water inflows within the fractured layer revealed that the most densely fractured and therefore productive part of the fractured horizon is located in the first 35–40 m below the base of the saprolite. These are emergence springs⁴⁶, the discharge of which depends on the thickness of the aquifer slice tapped. The information gathered during our various field missions reveals that most of these springs are developed by local populations, and generally capture locally a thickness of 2–3 m of the aquifer. This is consistent with the low discharge (0.3–2 L·s⁻¹) measured at these springs during our field-work (during low water stage). These springs, rarely isolated, are well distributed along the outcropping area of the saprolite aquifer (Fig. 1), at medium elevation in the valleys. It is noticeable that no spring was observed near the main watercourses in the downstream part of the watershed; this is due either to the fact that there are diffuse groundwater outflows along the watercourse and in the wetlands, and/or that the river eroded there the weathering profile below the base of the saprolite. Spatially distributed river discharge measurements and groundwater modeling should answer this question in a further stage of the research. Detailed discharge measurements performed during the 2022 low water stage in the Tazibouo and Balouzon area (Fig. 1) show that the 4 springs of this area have a total discharge of 2 L·s⁻¹ whereas the stream, immediately downstream them, has a 4.6 L·s⁻¹ discharge (the total discharge of this area being about 6.6 L·s⁻¹, which means about 1.6 L·s⁻¹·km⁻²; or about 50 mm/year, which gives by the way a minorant of the annual aquifer recharge). These data confirm that the springs only partially drain the aquifer, about two thirds of the groundwater “escaping” them and diffusely emerging in the wetlands and along the watercourse. A more suitable development of the springs, by increasing for instance the thickness of the slice of aquifer tapped, would therefore increase their discharge, but transform them into kind of dugwels, and reduce the stream discharge. This increase of the springs’ discharge would be a maximum of about 4 L·s⁻¹. It means that in these areas the potential for spring development is about 1.5 L·s⁻¹·km⁻² (or 50 mm/year). Compared with boreholes, spring exploitation is subject to the disadvantages of reduced discharge during the dry season; moreover, springs buffer inter-annual recharge variations to a lesser extent and
does not, as borehole does, allow for peaks in water uses. From a quality point of view, the isalterite aquifer is exposed at the surface; its groundwater is vulnerable. Other parts of this research project will deal with groundwater quality and clarify this statement.

The rare resistant anomalies visible on the ERT in the isalterite are interpreted either as boulders or to the presence of veins. Indeed, numerous quartz blocks, possibly the result of the dismantling of veins, are visible within the lateritic layers (See Appendix 6 in the Complementary Materials). On this site, both outcrop observations and the geoelectric signature rule out the presence of a laminated layer at the base of the saprolite³⁷.

## Conclusion

The development of a conceptual hydrogeological model of springs in this research responds to the urgent need to mobilize and effectively manage groundwater resources to support and maintain the supply of drinking water to populations in the crystalline regions of Côte d'Ivoire. This study explored the effectiveness of ERT in characterizing the context of groundwater emergence in weathered crystalline environments. To this end, ERT were sited on the plateaus and in the valleys where springs appear. The geoelectric layers revealed were compared with outcrops and observations performed during boreholes' drilling. A statistical analysis of the linear discharge of 316 deep boreholes and of the depths of water inflows was also carried out across the Daloa region to determine the thickness of the most productive part of the fractured layer. The main conclusions are as follows:

* ERT is effective in describing the stratigraphic structures and hydrogeological characteristics of the weathering profile. However, ERT cannot image the individual fractures of the fractured layer that were observed during drilling.

* The clayey sand, sand and granitic arenas that constitute the isalterites aquifer supplies the springs. The underlying fractured layer, which top is 20-30 m deeper than the springs level, is not directly implied in this supply.

* Statistical analysis of linear discharge and of the vertical distribution of water inflows beneath the saprolite shows that the most permeable part of the fractured layer is located in the first 35-40 m below the base of the saprolite. This layer is tapped by most borewells.

Although 2D electrical resistivity techniques are effective in identifying weathering profile structures, including zones of intense fracturing (the fractured horizon), they are limited in imaging individual fracture planes and quantifying the volume of mobilizable water in the aquifer. We propose further geophysical investigations using 3D electromagnetic imaging and Proton Magnetic Resonance (PMR) to quantify the volume of water stored in the aquifer, and hydrogeological modeling to quantify the volume of mobilizable groundwater.

## Methods

The methodological approach adopted consisted in describing the hydrogeological structure (weathering profile) of the aquifer system using electrical resistivity tomography (ERT) calibrated on borehole and outcrop data, then characterizing the piezometry.

### Description of weathering profile at borehole scale

The description of the vertical structuring of the weathering profile is based on the interpretation of 316 borehole data sheets with an average depth of 68 m distributed over the Lobo River watershed, including 16 in the study area (some of which are shown in Fig. 2). Drilling data sheets were collected from the Human Hydraulics Department (DHH), based in Daloa. These drilling data sheets are similar to those described in details by Aoulou et al.⁴. Interpretation of the lithologies was enriched by the observation of fractured and weathered crystalline rock outcrops and sand quarries.

### Fractured layer thickness

The method implemented is based on recent concepts of permeability linked to weathering in crystalline aquifers and the method proposed by Courtois et al¹⁶ and adopted by Aoulou et al.⁴ in a geological environment similar to the study site. One of the main results of the method was the characterization (thickness, productivity) of the fractured layer in crystalline aquifers. Productivity is expressed in terms of linear discharge (m³/h/m) and its variation was observed as a function of drilling depth beneath the base of the saprolite. On the basis of the conceptual model of hard rock aquifers, Courtois et al.¹⁶ defined a new parameter that they named "linear discharge" (Eq. 1):

$$q_i(l_i) = \frac{Q_i}{l_i}, i = 1, n \quad (1)$$

where $Q_i$ is the instantaneous discharge (m³/h) of borehole number i; $l_i$ is the length of the borehole i below the base of the saprolite (m); $q_i(l_i)$ is the resulting linear discharge of the borehole i (m³/h/m) and n is the number of borehole considered.

The cumulative percentage of linear discharge under the saprolite base was determined and plotted as a cumulative curve from Eq. (2):

$$p_q(L) = \frac{\sum_{l=L_{\min}}^{l=L} q_i(l)}{\sum_{l=L_{\min}}^{l=L_{\max}} q_i(l)} \quad (2)$$

where $q_i(l)$ is the linear discharge of borehole i calculated with borehole length l under the saprolite base (m³/h/m) and $p_q(L)$ is the cumulative percentage of linear discharge (%) calculated for the sample composed of boreholes whose length under the saprolite base is less than L (m). The cumulative percentage $p_q(L)$ is calculated
for all possible values of L present in the borehole dataset, ranging from Lmin to Lmax. By definition, the cumulative percentage ranges from p<sub>q</sub> (L<sub>min</sub>) = 0% to p<sub>q</sub> (L<sub>max</sub>) = 100%.

A graphical representation of cumulative linear discharge as a function of borehole length l below the base of the saprolite has been produced as proposed by Courtois et al.<sup>16</sup>. Within the cumulative linear discharge curve, a median part of the curve with a linear trend is generally distinguished, with a relatively gentle slope, corresponding to the highest linear discharge, which characterizes the thickness and productivity of the most productive part of the fractured layer. Beyond this point, linear discharges decrease and become zero as the base of the fractured layer is passed. The latter part of the curve characterizes the deepest part of the fractured layer and any underlying permeable discontinuities. The first part of the curve is sometimes steeper at the beginning. This is due to uncertainty as to the location of the saprolite base. The average productivity of the most transmissive part of the fractured layer can be quantified by Eqs. (3) and (4):

$$q_M(L_u) = \sum_{l=\text{Min}(l_i, i=1, n)}^{l=\text{Max}(l_i, i=1, n)} q_i(l)/j(L_u), \quad i=1, n \quad (3)$$

$$Q_M(L_u) = q_M(L_u) \times L_u \quad (4)$$

where L<sub>u</sub> (m) is the "useful thickness or thickness of the most productive part" of the fractured layer (defined either by the slope method or the percentile method); q<sub>M</sub> (L<sub>u</sub>) is the average linear discharge (m<sup>3</sup>/h/m) evaluated for the "useful thickness (L<sub>u</sub>)"; j (L<sub>u</sub>) is the corresponding number of boreholes; and Q<sub>M</sub> (L<sub>u</sub>) is the resulting average discharge (m<sup>3</sup>/h) for "theoretical" boreholes that would intersect the entire "useful thickness" (length L<sub>u</sub>) of the fractured layer. In the case of this study, L<sub>u</sub> was determined using the slope method as recommended by Aoulou et al.<sup>4</sup>, who have worked in geological environments similar to the study site.

The database contains information on the depth of water occurrences observed during drilling. Water inflows are observed exclusively within the fractured layer (below the base of the saprolite). They are linked to the existence of permeable fractures within the fractured layer. The database contains 316 boreholes well distributed over the Lobo River watershed, with a total of 593 water inflows observed. Firstly, the depth of each water inflow below the base of the saprolite was calculated as the difference between the depth of the water inflow in question and the depth of the base of the saprolite in the borehole in question. Next, the density of water inflows below the saprolite base was estimated using the R-4.2.1 software by the kernel method<sup>47</sup>, which is a generalization of the histogram estimation method. Kernel estimation is a non-parametric method for estimating the probability density of a random variable. The estimated water density follows a Gaussian distribution with parameters (mean = 23.64 m; standard deviation = 17.22 m). Finally, the joint analysis of the distribution of water inflows and the linear discharge curve enabled us to determine the most productive part of the fractured layer. The results obtained are presented in Appendix 7 in Complementary Materials.

## Description of weathering profile with ERT

The methodological approach consisted in creating electrical panels at ridge level (on boreholes F4 and F15, known as "test boreholes") to describe the weathering profile of plateau zones. Next, electrical panels were installed at the edge of the valleys (preferential groundwater zones) to describe the weathering profile of the spring zones. Finally, a comparative analysis of the geoelectric signatures was carried out to deduce the geological formations at the origin of the springs. 2D electrical resistivity tomography (ERT) was used to describe the weathering profiles. The implementation of this technique is divided into two main stages: field investigations and data processing (inversion, calibration and classification).

### Field investigations

The most commonly used arrays for 2D resistivity measurements are the Wenner, Wenner-Schlumberger, pole-pole, dipole-dipole and pole-dipole arrays<sup>48</sup>. In this research, the electrode configuration chosen was the pole-dipole, with 48 electrodes spaced 5 m apart. This choice stemmed from the strong urbanization of the study area, limiting the space available to extend electrical cables due to roads, buildings and buried electrical cables. For this reason, the choice fell on the pole-dipole array, known for its ability to differentiate complex geological structures at great depth with short ERT cables, as highlighted notably by Kumar et al.<sup>49</sup>. In addition, the pole-dipole array has a significantly higher signal strength than other arrays to obtain high-resolution 2D resistivity data, as well as higher vertical sensitivity and greater depth of investigation on short linear arrays, as confirmed notably by Kumar et al.<sup>49</sup> and Olivier et al.<sup>32</sup> who worked in similar geological environments. At site level, two ERT profiles (SP_F4 and SP_F15) were carried out at the plateau level to characterize the geoelectrical succession in the vicinity of boreholes F4 and F15. F4 and F15 are two boreholes in the University of Daloa's water catchment area, designed to provide the campus with a self-sufficient supply of drinking water. Six other ERT, namely SA1, SA2, SZ, STU1 and STU2, passing close to the Labia, Abattoir 2, Abattoir 2 SUD B, Zaragoua and Tazibouo-U springs respectively, were carried out in the dry season, during the low-water period, to characterize the weathering profile in the areas where the springs emerge. The aim of the ERT carried out at plateau level is to ensure that all stratiform geological formations in the weathering profile encountered in the area are intersected. This makes it possible to assess the geophysical signature describing the different lithological formations encountered in these regions, and to make a comparison between the geophysical signature of the ERT made near the test boreholes and that of the ERT made in the areas where the springs emerge. This technique made it possible to calibrate the model using a set of inversion parameter adjustments, while searching for the model that best describes the lithology of the test boreholes.

ERT were acquired using a Syscal-Pro Switch48 resistivity meter (Iris Instruments, France) with 48 metal electrodes regularly spaced 5 m apart over 235 m. The 5-m spacing between electrodes was chosen to ensure
both adequate resolution of subsurface resistivity at shallow depths and a maximum depth of investigation of around 90 m. Electrode contact resistance is generally less than 3 kΩ for all profiles. Acquisition parameters for all profiles and different configurations were as follows: acquisition time 30 min (1181 measurements in direct mode and 1181 measurements in reverse mode); Q < 1% (where Q is the standard deviation on resistance measurements, a repeatability characteristic) and injection voltage of 400 Volt. A high Q factor is generally an indication of a low signal-to-noise ratio. The pole-dipole protocol, with a remote injection electrode, maximizes the measurements signal-to-noise ratio and increases the depth of investigation. The position of the remote electrode was determined in accordance with the recommendations proposed by Razafindratsima and Lataste⁵⁰ to minimize its influence on resistivity measurements, with a minimum distance from the remote electrode of 1000 m (perpendicular to the ERT profile).

### Data processing: inversion, calibration and classification

After acquisition, measurements were filtered to eliminate anomalous values. Raw data from electrical measurements (apparent resistivity) were filtered and processed (correction of inconsistencies,) with Prosys II software (Iris Instruments) on the following criteria: V (signal) > 0.3 mV and Q < 1%.

Data processing includes data inversion, calibration and classification. Data inversion consists in obtaining iteratively interpreted resistivities from the apparent resistivities measured in the field. The inversion routine generally used is based on smoothing using the least squares constraint⁵¹. The objective function below (Eq. (5)) constrains the inversion in such a way that spatial variations in the calculated resistivities are smoothed.

$$f = (\rho_{\text{cal}} - \rho_{\text{mes}})^2 + \lambda \left( \frac{d\rho}{dx} + Z \frac{d\rho}{dy} \right) \quad (5)$$

With $\rho_{\text{cal}}$: calculated resistivity, $\rho_{\text{mes}}$: measured resistivity; $\rho$: resistivity of the reference model. The parameters $\lambda$ and $Z$ are, respectively, the damping factor used to define the intensity of the "smoothing" stress in the inversion and the parameter controlling the anisotropy.

The quality of the inversion is monitored by the RMS (Root Mean Squared)⁴⁸,⁵¹. RMS measures the sum of the square of the differences between measured and calculated apparent electrical resistivities. The Root Mean Square (RMS) percentage is used to evaluate inversion performance³¹. RMS is defined by the following Eq. (6):

$$RMS = \sqrt[100]{\frac{1}{N} \sum_{1}^{N} \left[ \frac{(h_0)_i - (h_m)_i}{(h_0)_i} \right]^2} \quad (6)$$

where $h_0$ is the observed depth (i.e. the true model) of the soil layer at location i, $h_m$ is the model (i.e. the inverted model) and $N$ is the number of observations (the total number of measured data sets).

If RMS is equal to 0, it is clear that the results simulated by the model correspond perfectly to the observations; if RMS is greater than 100%, then the model is no better than a predictor using zero³¹,⁵².

Electrical resistivity inversions were performed using Res2Dirv V-3.57.36 software. The aim of the inversion process is to minimize the objective function, which is the root-mean-square (RMS) error between the modeled and measured resistivities. Minimization is achieved by an iterative process in which the objective function must decrease. Convergence of the minimization is a mandatory criterion. This is because there is an infinite number of models available corresponding to measurements with the same RMS error⁴⁸,⁴⁹,⁵³. It is therefore important to reduce the number of mathematical solutions to models that are more geologically realistic, using a so-called inversion procedure and incorporating all available prior information into the inversion process. Caterina et al.⁵⁴ recommend adding prior information such as structural constraints, which seems to be particularly appropriate when the boundary is located in areas that are not very sensitive to ERT. However, such prior information is generally absent, and it is obvious that there is no specific set of parameters giving the "real" 2D resistivity distribution. It is therefore proposed to perform various inversions with different methods and specific constraints, so that each inversion imagery is particularly suited to determining the lithological and structural characteristics of the hydrogeological model³⁷,⁴⁸,⁵⁵⁻⁵⁷. Following this recommendation, we used two inversion methods for comparison: robust inversion, or L1-norm, and standard least-squares constraint inversion, or L2-norm⁵⁸. In this study, we refer to robust inversion (L1-norm) as "robust inversion methods" and standard least-squares constraint inversion (L2-norm) as "standard inversion methods". The residual values of the models for L1-norm and L2-norm inversion are expressed as absolute errors (Abs) and root mean square errors (RMS) respectively⁵⁸. According to Marescot⁴⁸, the robust inversion method (ROB) is well suited to environments with strong resistivity contrasts. It minimizes absolute error and provides clearer boundaries between varied layers. ROB inversion provides more accurate results for determining the boundary between the weathering profile and the unweathered hard rock. The standard inversion method (STD) minimizes squared error and is a smoothed inversion. Combined with an anisotropy factor (flatness ratios) of 0.3, this inversion reinforces the horizontal continuity of the geoelectric layers. In this study, each selected inversion method (Robust or Standard) is used to identify and characterize both vertical resistivity anomalies and recognize subhorizontal resistivity contrasts related to the stratiform weathering profile, in both cases using an anisotropy factor of 1. The structures thus revealed with each inversion method are compared with the lithology of the test boreholes. The resistivity values and geometries are then compared with the geological facies of the test boreholes to confirm the reality of the structures, particularly those that best describe the lithology of these boreholes. The resistivity patterns obtained after inversion were classified according to the resistivity ranges corresponding to each geological formation observed on the geological model of the test boreholes. For each spring ERT, resistivity ranges are established and compared with the resistivities and lithology of the test boreholes, before proceeding with the geological and hydrogeological interpretation.
## Piezometric mapping

Groundwater piezometric levels were measured using a probe whose sensor is reactive to the electrical conductivity of the water. When the probe reaches the level of the air/water interface, electrical contact is established between two metal rods, triggering a sound and light signal. The measurements were carried out in December 2021 (low-water period) on 14 boreholes. The boreholes monitored are little disturbed by domestic water abstraction, so the piezometric levels are not significantly influenced. The first step after the measurement campaign was to determine the piezometric values for each borehole surveyed. The piezometric values were determined using Eq. (7). The piezometric map was produced using the kriging method in Surfer 19, integrating the elevations of perennial springs (28 springs), permanent watercourses (210 points digitalized along the perennial part of watercourses) and the piezometric level measured on the 14 boreholes. Elevations were obtained from the Digital Elevation Model (DEM) extracted from STRM (Shuttle Radar Topographic Mission), with a XY resolution of 30 m. It is important to note that the survey of perennial springs and permanent rivers were carried out between December 2021 and February 2022, during the low-flow period. In this study, the type of kriging used is “point kriging” using the “ordinary kriging method”, as it is the most frequently used by many authors due to its ability to give accurate results⁵⁹,⁶⁰. Interpolation was carried out using the “standard ordinary kriging” option in Surfer. Prediction efficiency was calculated from cross-validation for the entire dataset. Analysis of the cross-validation parameters indicates a Root Mean Square (RMS) of 1.13 m, a variance of 1.28 m, a standard deviation of 1.13 m and a mean of 0.17 m. These accuracy characteristics prove to be satisfactory for good interpolation⁶¹. Although the cross-validation parameters are satisfactory, the interpolation is an estimate of point values, not exact values at all points. As a result, the value extrapolated from one point to another is not necessarily identical to the value that would be found if a well were drilled at that point. Moreover, interpolations at the edges of the map have a greater uncertainty than those at the center, as there are fewer field data around them. On the other hand, the accuracy is better in the center of the map. Then, on that basis, we manually adjusted the piezometric contours, using the Spline Polyline tool in Surfer, to account for the relationships between the piezometry and the topography that was not considered during kriging.

$$H = Z - (P - H_m) \qquad (7)$$

With: H: Piezometric elevation (m); Z: Ground elevation or elevation of natural terrain (m); P: Measured depth (m); Hm: well curb height.

## Data availability

The data used in the study are available on request from the first author.

Received: 3 July 2023; Accepted: 31 May 2024

Published online: 05 June 2024

## References

1. WWAP (United Nations World Water Assessment Programme). The United Nations World Water Development Report 2015: Water for a Sustainable World. Paris, UNESCO; Preprint at https://unesdoc.unesco.org/ark:/48223/pf0000231823 (2015).

2. Giao, N. T. et al. Groundwater quality assessment for drinking purposes: A case study in the Mekong Delta, Vietnam. Sci. Rep. 13, 4380. https://doi.org/10.1038/s41598-023-31621-9 (2023).

3. Siebert, S. et al. Groundwater use for irrigation - a global inventory. Hydrol. Earth Syst. Sci. 14, 1863-1880. https://doi.org/10.5194/hess-14-1863-2010 (2010).

4. Aoulou, K. A., Pistre, S., Oga, Y. M. S., Dewandel, B. & Lachassagne, P. Improving the methods for processing hard rock aquifers boreholes' databases. Application to the hydrodynamic characterization of metamorphic aquifers from Western Côte d'Ivoire. Water 13, 3219. https://doi.org/10.3390/w13223219 (2021).

5. Lachassagne, P., Dewandel, B. & Wyns, R. Review: Hydrogeology of weathered crystalline/hard-rock aquifers-guidelines for the operational survey and management of their groundwater resources. Hydrogeol. J. 29, 2561-2594. https://doi.org/10.1007/s10040-021-02339-7 (2021).

6. Tellam, J. H. & Barker, R. D. Towards prediction of saturated-zone pollutant movement in groundwaters in fractured permeable-matrix aquifers: The case of the UK Permo-triassic sandstones. Geol. Soc. London, Special Publications 263(1), 1-48. https://doi.org/10.1144/GSL.SP.2006.263.01.01 (2006).

7. Medici, G. & West, L. J. Review of groundwater flow and contaminant transport modelling approaches for the Sherwood sandstone Aquifer, UK; insights from analogous successions worldwide. Q. J. Eng. Geol. Hydrogeol. 55(4), qjegh2021-176. https://doi.org/10.1144/qjegh2021-176 (2022).

8. Medici, G. & West, L. J. Reply to discussion on 'review of groundwater flow and contaminant transport modelling approaches for the sherwood sandstone Aquifer, UK; insights from analogous successions worldwide' by Medici and West (QJEGH, 55, Qjegh 2021-176). Q. J. Eng. Geol. Hydrogeol. 56(1), 022-097. https://doi.org/10.1144/qjegh2022-097 (2023).

9. Dewandel, B. et al. Upscaling and regionalizing hydraulic conductivity and effective porosity at watershed scale in deeply weathered crystalline aquifers. J. Hydrol. 416(471), 83-97. https://doi.org/10.1016/j.jhydrol.2011.11.038 (2012).

10. Gleeson, T., Novakowski, K. & Kyser, T. K. Extremely rapid and localized recharge to a fractured rock aquifer. J. Hydrol. 376, 3-4. https://doi.org/10.1016/j.jhydrol.2009.07.056 (2009).

11. Lods, G. & Gouze, Ph. WTFM, software for Well Test analysis in Fractured Media combining fractional flow with double porosity and leakance approaches. Comput. Geosci. J. 30, 937-947. https://doi.org/10.1016/j.cageo.2004.06.003 (2004).

12. Lachassagne, P. et al. Recent improvements in the conceptual model of hard rock aquifers and its application to the survey, management, modelling and protection of groundwater. Groundwater and Climate in Africa (Proceedings of the Kampala Conference, June 2008) IAHS Publ. 334; https://www.researchgate.net/publication/228817904 (2009).

13. Dewandel, B. et al. Development of a tool for managing groundwater resources in semi-arid hard rock regions: Application to a rural watershed in South India. Hydrol. Process. https://doi.org/10.1002/hyp.7696; 10.1002/hyp.7696 (2010).

14. Hu, K., Chen, H., Nie, Y. & Wang, K. Seasonal recharge and mean residence times of soil and epikarst water in a small karst catchment of southwest China. Sci. Rep. 5, 10215. https://doi.org/10.1038/srep10215 (2015).

15. Dewandel, B. et al. A methodology for regionalizing 3D effective porosity at watershed scale in crystalline aquifers. Hydrol. Process. 31, 2277-2295. https://doi.org/10.1002/hyp.11187 (2017).
16. Courtois, N. et al. Large-scale mapping of hard-rock aquifer properties applied to Burkina Faso. *Ground Water* **48**, 269–283. https://doi.org/10.1111/j.1745-6584.2009.00620.x (2010).

17. Vouillamoz, J. M., Lawson, F. M. A., Yalo, N. & Desclaitres, M. The use of magnetic resonance sounding for quantifying specific yield and transmissivity in hard rock aquifers: The example of Benin. *J. Appl. Geophys.* **107**, 16–24. https://doi.org/10.1016/j.jappgeo.2014.05.012 (2014).

18. Maurice, L. et al. Characteristics of high-intensity groundwater abstractions from weathered crystalline bedrock aquifers in East Africa. *Hydrogeol. J.* **27**, 459. https://doi.org/10.1007/s10040-018-1836-9 (2018).

19. Ligban, R., Gone, D. L., Kamagate, B., Saley, M. B. & Biemi, J. Processus hydrogéochimiques et origine des sources naturelles dans le degré carré de Daloa. *Int. J. Biol. Chem. Sci.* **3**(1), 38–47. https://doi.org/10.4314/ijbcs.v3i1.42733 (2009).

20. Adjiri, O. A., Koné, B., Aka, N., Djabakate, I. & Dibi, B. Caractérisation physico-chimique et source de la minéralisation des eaux souterraines des départements de Daloa et Zoukougbeu, Côte d'Ivoire. *Int. J. Biol. Chem. Sci.* **13**(4), 2388–2401. https://doi.org/10.4314/ijbcs.v13i4.40 (2019).

21. Maréchal, J. C., Dewandel, B. & Subrahmanyam, K. Use of hydraulic tests at different scales to characterize fracture network properties in the weathered-fractured layer of a hard rock aquifer. *Water Resour. Res.* **40**(11), 17. https://doi.org/10.1029/2004W R003137 (2004).

22. Dewandel, B., Lachassagne, P., Zaidi, F. K. & Chandra, S. A. Conceptual hydrodynamic model of a geological discontinuity in hard rock aquifers: Example of a quartz reef in granitic terrain in South India. *J. Hydrol.* **405**(2011), 474–487. https://doi.org/10.1016/j.jhydrol.2011.05.050 (2011).

23. Lachassagne, P., Wyns, R. & Dewandel, B. The fracture permeability of hard rock aquifers is due neither to tectonics, nor to unload- ing, but to weathering processes. *Terra Nova.* **23**, 145–161. https://doi.org/10.1111/j.1365-3121.2011.00998.x (2011).

24. Wyns, R. et al. Application of proton magnetic resonance soundings to groundwater reserve mapping in weathered basement rocks (Brittany, France). *Bulletin de la Société Géologique de France* **175**(1), 21–34. https://doi.org/10.2113/175.1.21 (2004).

25. Dewandel, B., Lachassagne, P., Wyns, R., Maréchal, J. C. & Krishnamurthy, N. S. A generalized hydrogeological conceptual model of granite aquifers controlled by single or multiphase weathering. *J. Hydrol.* **330**, 260–284. https://doi.org/10.1016/j.jhydrol.2006.03.026 (2006).

26. Worthington, S. R. H., Davies, G. J. & Alexander, E. C. J. Enhancement of bedrock permeability by weathering. *Earth-Sci. Rev.* **160**, 188–202. https://doi.org/10.1016/j.earscirev.2016.07.002 (2016).

27. St, J. et al. Geophysical imaging reveals topographic stress control of bedrock weathering. *Science* **350**(6260), 534–538. https://doi.org/10.1126/science.aab2210 (2015).

28. Allé, C. et al. Why 1D electrical resistivity techniques can result in inaccurate siting of boreholes in hard rock aquifers and why electrical resistivity tomography must be preferred: The example of Benin, West Africa. *J. African Earth Sci.* **139**, 341–353. https://doi.org/10.1016/j.jafrearsci.2017.12.007 (2018).

29. Soro, D. D. et al. Geophysical demonstration of the absence of correlation between lineaments and hydrogeologically usefull fractures: Case study of the Sanon hard rock aquifer (central northern Burkina Faso). *J. Afr. Earth Sci.* **129**, 842–852. https://doi.org/10.1016/j.jafrearsci.2017.02.025 (2017).

30. Alfy, M. E. et al. Quantitative hydro-geophysical analysis of a complex structural karst aquifer in Eastern Saudi Arabia. *Sci. Rep.* **9**, 2825. https://doi.org/10.1038/s41598-019-39192-4 (2019).

31. Cheng, Q. et al. Evaluation of electrical resistivity tomography (ERT) for mapping the soil-rock interface in karstic environments. *Environ. Earth Sci.* **78**, 439. https://doi.org/10.1007/s12665-019-8440-8 (2019).

32. Olivier, K. J. J. et al. Elaboration of a hydrogeological conceptual model by application of electrical resistivity tomography: Case of the Lobo catchment (Centre-Western Côte d'Ivoire). *Sci. African* **16**(7), 1–13. https://doi.org/10.1016/j.sciaf.2022.e01234 (2022).

33. Pradhan, R. M. et al. Structural controls on bedrock weathering in crystalline basement terranes and its implications on ground-water resources. *Sci. Rep.* **12**, 11815. https://doi.org/10.1038/s41598-022-15889-x (2022).

34. Acharya, B. S. et al. Vegetation controls on the spatio-temporal heterogeneity of deep moisture in the unsaturated zone: A hydro- geophysical evaluation. *Sci. Rep.* **7**, 1499. https://doi.org/10.1038/s41598-017-01662-y (2017).

35. Desclaitres, M. et al. Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding. *Hydrol. Process.* **22**(3), 384–394. https://doi.org/10.1002/hyp.6608 (2008).

36. Braun, J. J. et al. Regolith mass balance inferred from combined mineralogical, geochemical and geophysical studies: Mule Hole gneissic watershed South India. *Geochimica et Cosmochimica Acta* **73**, 935–961. https://doi.org/10.1016/j.gca.2008.11.013 (2009).

37. Belle, P. et al. Characterization and location of the laminated layer within hard rock weathering profiles from electrical resistivity tomography: Implications for water well siting. *Geol. Soc. Lond.* **479**, 187–205. https://doi.org/10.1144/SP479.7 (2019).

38. Ugbor, C. C., Ikwuagwu, I. E. & Ogboke, O. J. 2D inversion of electrical resistivity investigation of contaminant plume around a dumpsite near Onitsha expressway in southeastern Nigeria. *Sci. Rep.* **11**, 11854 (2021).

39. Zhou, B. & Dahlin, T. Properties and effects of measurement errors on 2D resistivity imaging surveying. *Near Surface Geophys.* **1**, 105–117. https://doi.org/10.3997/1873-0604.2003001 (2003).

40. Doro, K. O., Adegboyega, C. O., Aizebeokhai, A. P. & Oladunjoye, M. A. The Ibadan hydrogeophysics research site (IHRS)-An observatory for studying hydrological heterogeneities in A crystalline basement aquifer in southwestern Nigeria. *Water* **15**, 433. https://doi.org/10.3390/w15030433 (2023).

41. Delor, C., Siméon Y., Vidal M., Zeade, Z., Koné, Y. & Adou M. Carte géologique de la Côte d'Ivoire à 1/200 000, feuille Séguéla, Mémoire n°9 de la Direction des Mines et de la Géologie, (Abidjan, Côte d'Ivoire), p. 19 (1995).

42. Nagaiah, E., Sonkamble, S. & Chandra, S. Electrical geophysical techniques pin-pointing the bedrock fractures for groundwater exploration in granitic hard rocks of Southern India. *J. Appl. Geophys.* **199**, 10460. https://doi.org/10.1016/j.jappgeo.2022.104610 (2022).

43. Koita, M., Jourde, H., Koffi, K. J. P., Da Silvera, K. S. & Biaou, A. Characterization of weathering profile in granites and volcano- sedimentary rocks in West Africa under humid tropical climate conditions. Case of the Dimbokro Catchment (Ivory Coast). *J. Earth Syst. Sci.* **122**, 841–854. https://doi.org/10.1007/s12040-013-0290-2 (2013).

44. Lachassagne, P., Aunay, B., Frissant, N., Guilbert, M. & Malard, A. High-resolution conceptual hydrogeological model of complex basaltic volcanic islands: A Mayotte, Comoros, case study. *Terra Nova* **26**, 307–321. https://doi.org/10.1111/ter.12102 (2014).

45. Blot, A. Caractérisation des chapeaux de fer en milieu latéritique cuirassé. *C. R. Geosci.* **336**(16), 1473–1480. https://doi.org/10.1016/j.crte.2004.07.008 (2004).

46. Archambault, J. Le captage des sources. Preprint at https://www.pseau.org/outils/ouvrages/gret_afvp_le_captage_des.sources_1987.pdf (1987).

47. Simon, J. S. Density estimation. *Statist. Sci.* **19**(4), 588–597. https://doi.org/10.1214/088342304000000297 (2004).

48. Marescot, L. Imagerie électrique pour géologues: Acquisition, traitement, interprétation. Cours. Preprint at http://www.tomoquest.com/attachments/File/Geol_Intro_ERT.pdf (2008).

49. Kumar, D., Rao, V. A. & Sarma, V. S. Hydrogeological and geophysical study for deeper groundwater resource in quartzitic hard rock ridge region from 2D resistivity data. *J. Earth Syst. Sci.* **123**, 531–543. https://doi.org/10.1007/s12040-014-0408-1 (2014).

50. Razafindratsima, S. & Lataste, J. F. Estimation of the error made in pole-dipole electrical resistivity tomography depending on the location of the remote electrode: Modeling and field study. *J. Appl. Geophys.* **100**, 44–57. https://doi.org/10.1016/j.jappgeo.2013.10.008 (2014).
51. Günther, T. Inversion methods and resolution analysis for the 2D/3D reconstruction of resistivity structures from DC measure-
ments. PhD Thesis, University of Mining and Technology of Freiberg, Germany, 160 p; https://www.researchgate.net/publication/
44231115 (2004).

52. Cheng, Q. B. et al. Using maximum likelihood to derive various distance-based goodness-of-fit indicators for hydrologic modeling assessment. *Stoch. Environ. Res. Risk Assess* **32**, 949–966. https://doi.org/10.1007/s00477-017-1507-8 (2018).

53. Zhou, J. et al. Image-guided inversion of electrical resistivity data. *Geophys. J. Int.* **197**(1), 292–309. https://doi.org/10.1093/gji/ggu001 (2014).

54. Caterina, D., Hermans, T. & Nguyen, F. Case studies of incorporation of prior information in electrical resistivity tomography; comparison of different approaches. *Near Surface Geophys.* **12**, 451–465. https://doi.org/10.3997/1873-0604.2013070 (2014).

55. Dahlin, T. & Zhou, B. A numerical comparison of 2D resistivity imaging with ten electrode arrays. *Geophys. Prospect.* **52**, 379–398.
https://doi.org/10.1111/j.1365-2478.2004.00423.x (2004).

56. Jia, Z. et al. The deep origin of ground fissures in the Kenya Rift Valley. *Sci. Rep.* **13**, 3672. https://doi.org/10.1038/s41598-023-30918-z (2023).

57. Priou, J., Chevreuil, M. & Lecieux, Y. Utilisation d'un méta-modèle pour l'identification paramétrique: Application à la mesure de résistivité électrique. Preprint at https://www.researchgate.net/publication/336718766 (2019).

58. Guedes, V., Lima, V., Borges, W. & Cunha, L. Comparison of the geoelectric signature with different electrode arrays at the jockey club landfill of Brasília, DF, Brazil. *Brazil. J. Geophys.* **38**, 41. https://doi.org/10.22564/rbgf.v38i1.2034 (2020).

59. Raid, I., Kusnezowa, T. & Seewig, J. Application of ordinary kriging for interpolation of micro-structured technical surfaces. *Measur. Sci. Technol.* **24**(9), 095201. https://doi.org/10.1088/0957-0233/24/9/095201 (2013).

60. Mueller, T. G. et al. Map quality for ordinary kriging and inverse distance weighted interpolation. *Soil Sci. Soc. Am. J.* **68**, 2042–2047.
https://doi.org/10.2136/sssaj2004.2042 (2004).

61. Fencik, R., & Vajsáblová, M. Parameters of interpolation methods of creation of digital model of landscape. *Ninth AGILE Conference on Geographic Information Science* **20**, https://www.researchgate.net/publication/228438864 (2006).

Acknowledgements

Funding for this work was provided by PASRES (PASRES is the Strategic Support Program for Scientific Research in Ivory Coast). This study was also supported by UEMOA's scholarship program for excellence in training and research (UEMOA is West African Economic and Monetary Union). The authors would like to thank the funders for their help and support. We would also like to thank the Scholarship Department of the Ministry of Higher Education and Scientific Research of Ivory Coast for granting us a scholarship to study in France as part of this research project. The funders played no role in the design of the study, the collection and processing of data, the decision to publish or the preparation of the manuscript.

Author contributions

K.K.J.M., D.B., M.O.M.J. and L.P designed this research project. K.K.J.M., M.O.M.J. and S.A.P have collected the field data. K.K.J.M. has analyzed and processed the data. K.K.J.M., D.B. and L.P have interpreted the results and designed the visual elements of the document. All authors critically edited the manuscript.

Competing interests

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-63748-8.

Correspondence and requests for materials should be addressed to K.J.-M.K.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024