C. R.

Geoscience 339 (2007) 418–429

http://france.elsevier.com/direct/CRAS2A/

Géosciences de surface (Hydrologie–Hydrogéologie)

Processus et bilan des flux hydriques d’un bassin versant de
milieu tropical de socle au Bénin (Donga, haut Ouémé)
Bamory Kamagaté a, Luc Séguis a, Guillaume Favreau a,*, Jean-Luc Seidel a,
Marc Descloitres b, Pascal Affaton c
a
UMR HydroSciences (IRD, CNRS, UMI, UMII), université Montpellier-2, place Eugène-Bataillon, CC MSE,
34095 Montpellier cedex 5, France
b
IRD, UMR LTHE, BP 53, 38041 Grenoble cedex 09, France
c
Université Aix-Marseille, UMR CEREGE, Europôle de l’Arbois, BP 80, 13545, Aix-en-Provence cedex 4, France
Reçu le 26 juillet 2006 ; accepté après révision le 17 avril 2007
Disponible sur Internet le 30 mai 2007

Présenté par Ghislain de Marsily

Résumé
L’étude d’un bassin versant de 586 km2 au Bénin (Afrique de l’Ouest) par des suivis hydrodynamiques et hydrochimiques, ainsi
que par reconnaissance géophysique de subsurface, a permis de déterminer les processus et le bilan des flux hydriques sur deux
années à pluviométrie contrastée. La recharge de la nappe phréatique s’effectue par infiltration directe des précipitations et
représente de 5 à 24 % de la pluviométrie. L’écoulement à l’exutoire, limité à la saison des pluies ainsi que la minéralisation faible
des eaux indiquent des flux nappe–rivière négligeables. Cette faible contribution de la nappe implique que les variations
interannuelles du coefficient d’écoulement (de 14 à 28 %) sont principalement gouvernées par des flux de surface et de subsurface.
Pour citer cet article : B. Kamagaté et al., C. R. Geoscience 339 (2007).
# 2007 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés.

Abstract
Hydrological processes and water balance of a tropical crystalline bedrock catchment in Benin (Donga, upper Ouémé
River). Hydrodynamic, geochemical, and subsurface geophysical investigations, for two consecutive years with contrasting rainfall
conditions, were used to characterize the hydrological processes occurring, and the water balance of a 586-km2 watershed in Benin
(Africa). The water table’s monitoring shows that recharge occurs by direct infiltration of rainfall, and represents between 5 to 24%
of the annual rainfall. Both surface water outflow, limited to the rainy season, and water chemistry indicate a weak groundwater
contribution to river discharge. This implies that the calculated variations in annual runoff coefficients (of 14 and 28%) are mainly
governed by surface and subsurface flows. To cite this article: B. Kamagaté et al., C. R. Geoscience 339 (2007).
# 2007 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés.

Mots clés : Bénin ; Aquifère libre ; Hydrochimie ; Hydrodynamique ; Rivière tropicale

Keywords: Benin; Groundwater; Hydrochemistry; Hydrodynamics; Tropical river

* Auteur correspondant.
Adresse e-mail : [email protected] (G. Favreau).

1631-0713/$ – see front matter # 2007 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés.
doi:10.1016/j.crte.2007.04.003
 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429 419

Abridged English version basin. Our surveys indicate a water use of 21 L

inhab day 1, a value close to the 17 L inhab day 1
Introduction estimated at a regional scale [15]. Considering the
population density, this implies a human-induced
In West Africa, the rainfall deficit of about 20 to 30% discharge from the aquifer close to 0.2 mm yr 1.
observed since the early 1970s [18,21] has led to a
higher decrease in river discharge [23]. This implicitly Data
non-linear relationship was mainly studied for large
river basins, by indirect or conceptual methods [23,24]; Electrical resistivity data were obtained in May 2003,
the rationale is that lower rainfall resulted in lower with a view to characterize the weathered substratum;
groundwater discharge into the rivers. However, studies changes were expected for differences in rock type, clay,
performed at a finer scale (e.g., [6]) suggested more and/or water content [8]. The electrical tomography
complex surface–subsurface flow processes. method [22] was implemented at a scale of several
The AMMA project (African Monsoon Multidisci- hundred meters on few representative slopes; a geolo-
plinary Analysis, http://amma-international.org) was gical survey helped in constraining the interpretation.
designed to better understand the African monsoon, and Hydrodynamics data records were obtained on 16
its hydrological impacts. For this goal, surveys have rain gauges (each 0.5 mm) and 26 wells (each 30 min
been performed since 2003 in a 586-km2 catchment in or three times a day). Additional data from an upstream
central Benin (Fig. 1), a key-region for understanding site were obtained in 2004 from nine piezometers
the impact of climate [19], and land-use changes [14] on screened at different depths (2, 10, and 20 m; Fig. 1). For
the regional hydrological cycle. surface waters, three limnimetric records (time steps of 5
This paper presents the first 2-year hydrodynamic, or 30 min) represent catchments of respectively 13, 105,
geochemical, and subsurface geophysical surveys and 586 km2 (Fig. 1). A differential GPS campaign was
performed within this basin. A first estimate of the also performed to characterize the hydraulic potentials
water balance is obtained for these two years, for which [16].
contrasting rainfall amounts were recorded. Hydrochemical data (temperature, pH, EC 25 8C)
were obtained monthly or quarterly on 20 wells; some
Study area measurements were also performed on bulk rainfall
events. Measurements on surface waters were made for
The Donga River (Fig. 1) is an affluent of the Ouémé springs, tributaries and on the Donga River. The
River (523 km), in central Benin [17]. The tropical limnimetric station K was equipped in June 2003 with
climate is characterized by a dry and a rainy season, an EC recorder (time step of 30 min). Other EC recorders
with a mean annual rainfall measured at Djougou of with the same time step were settled in April 2004 on the
1195 mm (1950–2004). The potential evapotranspira- nested piezometers on the Ara sub-catchment. Chemical
tion (ETP, Penman–Monteith [10]) is of 1390 mm yr 1 analysis (Ca2+, Mg2+, Na+, K+ for cations, HCO3 ,
(2002–2005). SO42 , NO3 , Cl for anions) were performed on
Altitude ranges from 520 to 340 m for a river slope of samples from rainwater (18 analyses), groundwater (23
1.7 m km 1. Outcropping siliceous rocks are of wells and 15 boreholes), seasonal springs (8 analyses),
metamorphic origin [1,5]. The weathered bedrock has and water from the Donga River (a dozen of analyses for
a thickness of 10 to 25 m, and represents the matrix of each station); in addition, three main flood events were
the phreatic aquifer. In depth, some faults also contain surveyed by sequential sampling in 2003. All the samples
deeper groundwater, in hydraulic continuity with the were filtered in the field and subsequently analysed by
overlaying reservoir [16]. The soils are classified as capillary electrophoresis (cations and SO42 , NO3 , Cl )
polycyclic ferruginous soils [12], but some details may or acidic titration (HCO3 ), with a precision of about 5%.
change depending on the topography, rock type, or land
use [11]. The natural wooded savannah has been largely Results
cleared, and most of the land is now covered by rain-fed
crop fields; a forest gallery of a few tens of metres runs Subsurface geophysics
along most of the streams. A typical example of the resistivity profiles performed
The population density is of 30 inhab km 2, and has within the limits of the basin [8,16] is displayed in Fig. 2;
increased recently by 4.9% yr 1 (1992–2002). Most of it reveals an important heterogeneity of the resistivity at
the domestic water comes from wells scattered in the the hectometric scale. The bedrock lies from 15 to 25 m
420 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429

below the soil surface. Clayey materials (electrically Discussion

conductive) and/or indurated covers, or quartzite out-
crops (both resistive) occur at surface and may locally Hydrological processes
outcrop; the clayey layers act as preferential places for Hydrodynamics, hydrochemical, and geophysical
the formation of perched groundwaters. data were crossed to decipher the natural processes in
play in the basin (Fig. 6).
Hydrodynamics For groundwater flow, the water table is always in a
The total rainfall amount was of 25% lower in 2004 higher position than the drainage network. The
than in 2003, but the difference is even higher for total weathered bedrock acts as a storage reservoir, with a
runoff, with a decrease of 64% (Table 1). The river flow weak transmissivity (estimated of the order of
shows rapid changes, with daily discharge reaching a 10 4 m2 s 1, [16]). During the rainy season, a water
maximum of 140 m3 s 1 at the basin outlet; at the table rise is observed in each point, thus indicating that
seasonal scale, the discharge period was shorter in 2004 recharge occurs mainly by vertical infiltration of
than in 2003, under the influence of the seasonal rainfall rainfall.
distribution (Fig. 3). The seasonal discharge of poorly mineralized river
The depth-to-water table ranges from less than 1 m waters suggests that most of the river water originates
(rainy season) up to 15 m (dry season), without any from the subsurface. Observations of springs on
spatial relationship within the basin. The potentiometric hillslopes and the occurrence of perched groundwaters
levels are always higher than the riverbed; this both suggest important subsurface flow. The river
observation is also valid for the more densely monitored discharge may therefore be considered as a mixing of
Ara subcatchment (Fig. 1), where seasonal perched surface runoff (Hortonian runoff and flow over saturated
groundwater was recognized. The water table peak is soils) with subsurface flow from the upper metres of the
reached in August, and the drop begins simultaneously soil (Fig. 6). The aquifer discharge may occur mainly
with the rainfall (Fig. 3C). The water table fluctuations through evapotranspiration processes, although locally,
range between 3.3 and 12.5 m in 2003 and between 2.4 groundwater discharge could also occur through some
and 10.0 m in 2004, and appear synchronous at the deep fractures of higher porosity [1].
catchment scale.
Water balance
Hydrochemistry Two methods were applied to constrain the hydro-
Sampled waters are slightly acidic (pH 5–7), present a logical balance. The first one considers natural water
low mineralization (<400 mS cm 1; Fig. 4), and are of table fluctuation induced from rainfall events over a
HCO3–Ca, Mg type (some groundwaters enriched in Cl short period [7]: R = h  n, where R represents recharge
and NO3 near villages were considered polluted). The (L T 1), h the water table rise (L T 1), and n the specific
lowest EC is observed in rainfall waters (10  2 mS yield (%). For our study case, h is the average seasonal
cm 1). Except for the very first flood events that may rise of the water table (Fig. 3C). The calculated recharge
flush the salt accumulated for the dry season, EC for an n value therefore represents a minimum estimate,
measured in river waters is low (<70 mS cm 1). The EC considering that discharge may occur during the rainy
displays typical changes during flood events, with low season. A range of values of 1 to 5% [9,13] was
values at the peak flow, followed by a recovery of the considered to quantify the influence of n on R (Table 1).
initial value (Fig. 5A). The phreatic groundwater EC is of The estimated recharge ranges between 5 to 24% of the
155  70 mS cm 1, but EC is higher for the deeper total rainfall amount, in agreement with other regional
fractured bedrock (288  80 mS cm 1). Perched estimates [13]. These recharge rates are compatible with
groundwaters show lower EC, within the range of 30 discharge by evapotranspiration during the dry season.
to 130 mS cm 1 (Fig. 4). At slope scale, an estimate of groundwater discharge
Ionic ratios show little changes through time and, for using Darcy’s law (T of 10 4 m2 s 1, piezometric
river waters, elements related to chemical weathering gradients of 2–3%) shows that ETP rates of about 15–
(Si, Ca, Mg, Na, and HCO3) vary proportionally to EC 20 mm day 1 for forest galleries would be necessary to
(Fig. 5A); potassium is an exception, probably because prevent any inflow to the river (Fig. 6). Although this
of its sensitivity to biochemical cycles [4]. Changes in range appears to be higher than the computed ETP (5–
ionic fluxes appear to mimic the peak flow with a time 6 mm/day 1), clothesline and oasis effects are known to
lag of a few hours (Fig. 5B), a possible consequence of increase ETP rates in forest gallery, possibly at a factor
subcatchments with different weathering rates. higher than 2 [10]; further research would however be
 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429 421

needed to better constrain surface water–groundwater ont jusqu’à présent été menées à partir de données
relationships. hydrodynamiques acquises sur de grands bassins versants
The second approach, based on hydrograph decom- (de l’ordre de 10 000 à 100 000 km2), par des approches
position [2], was used to discriminate direct runoff and indirecte ou conceptuelle (tests statistiques ; calculs de
subsurface flow in river discharge. A two-component, coefficients de tarissement [23,24]). L’explication
one-tracer (EC) hydrograph separation model was proposée est que la diminution relativement plus
considered (EC being representative of the ionic load: importante des débits trouve son origine dans une
e.g., R2 of 0.91 between EC and Ca, Fig. 5A). The réduction durable des apports en eau souterraine. Si la
different signatures were assumed roughly constant in validation de cette hypothèse nécessite d’être effectuée
time and space; computations were done for the K station par une analyse des processus physiques, à une échelle
(Fig. 1). The end-members values (direct runoff and sub- comparable à celle des observations de surface (données
surface flux) were estimated to be of 10 and 55 mS cm 1, piézométriques disponibles sur plusieurs milliers de
with a standard-deviation of 2 and 6, respectively. A kilomètres carrés, sur la période pluridécennale des
Monte Carlo approach was used to estimate the observations de débits), des éléments de réponse peuvent
uncertainty; sub-surface flux makes up to 68  13 and cependant être apportés par la caractérisation des
83  13% of the stream flow in 2003 and 2004, processus physiques affectant les différents réservoirs
respectively. This suggests that the river flow deficit hydriques d’un bassin versant représentatif, étudié sur
has not a single origin, but results in a decrease in both des années à pluviométrie et écoulements contrastés. En
surface and subsurface flows. This conclusion differs Afrique subsaharienne, les rares études existant à fine
from those of other studies performed in the Sahel, where échelle (bassins versants expérimentaux de l’ordre du
Hortonian runoff was shown to have increased following kilomètre carré (e.g. [6]), montrent que la discrimination
soil crusting (e.g., [26]). des écoulements est améliorée par une approche
pluridisciplinaire des différents termes du bilan.
Conclusion Le programme international AMMA (Analyse
multidisciplinaire de la mousson africaine ; http://
In central tropical Benin, at the 586-km2 Donga amma-international.org) a été initié en vue d’améliorer
basin’s scale, runoff and subsurface flow are the main les connaissances sur la mousson ouest-africaine et de
contributors to stream flow, with little or no input from son impact hydrologique. Différents sites instrumentés
the aquifer. This result supplements previous studies font l’objet, depuis 2003, de suivis renforcés, dont celui
performed in smaller (30 km2) neighbouring catch- du bassin versant de la Donga (Ouémé supérieur) au
ments [14], where land clearance was shown to have Bénin (Fig. 1). Sa superficie (586 km2) en fait un objet
increased surface runoff, in accordance with the global
trend [25]. These results will first help in improving the
relevancy of the regional hydrological models
[3,20,27]; on a broader perspective, they will also
contribute to a better understanding of the origin of the
decrease in river discharge observed for the past
decades in West Africa [18,23].

1. Introduction

En Afrique de l’Ouest, sous régime de mousson, la

sécheresse en cours depuis le début des années 1970,
caractérisée par une diminution de l’ordre de 20 à 30 % de
la pluviométrie [18,21], a engendré une diminution
encore plus importante des débits [23]. Au centre du
Bénin, sur la zone d’étude considérée, le déficit des
Fig. 1. Localisation du bassin versant de la Donga au Bénin (encart,
écoulements, de l’ordre de 40 %, correspond à une symbole étoilé) et détails du dispositif de suivi hydrologique (2003–
diminution de la pluviométrie de 15 à 20 % [18]. Cette 2004).
relation implicitement non linéaire entre pluie et débit, Fig. 1. Location of the Donga watershed in Benin (black star symbol,
ainsi que sa variabilité spatiotemporelle, sont encore mal small inset) and details of the hydrological survey network (2003–
comprises à l’échelle régionale ; les principales études 2004).
422 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429

hydrologique intermédiaire entre les grands bassins supérieur et constituées de granites, de granito-gneiss et
régionaux et les bassins expérimentaux ; cette échelle de schistes [1,5]. Les altérites, principalement silteuses
rend possible une caractérisation des processus hydro- et argileuses sur une épaisseur de 10 à 25 m, constituent
logiques tout en permettant, en terme de bilan, une l’aquifère de la nappe libre. Localement et plus en
intégration suffisamment représentative des variations profondeur, des fractures du socle peuvent constituer un
spatiotemporelles naturelles. À l’échelle régionale, de réservoir pérenne, en continuité hydraulique avec la
même que l’analyse de la variation spatiotemporelle de couverture altérée du socle [16].
la pluviométrie [19] ou de la caractérisation spatiale des Les sols régionaux, de type « ferrugineux tropicaux
sols [14], la connaissance des processus surface– lessivés » [12] sont de nature variable en fonction de leur
souterrain est un élément fondamental pour appréhen- position topographique, de la nature des roches sous-
der la variabilité temporelle des ressources hydriques. jacentes ou de l’anthropisation du site (sols localement
Dans cette note, nous présentons les premiers hydromorphes, indurés, peu lessivés ou lessivés [11]). La
résultats des suivis hydrodynamiques (piézométrie et végétation climacique est constituée d’une savane à
débits) et hydrochimiques (physicochimie et ions Isoberlinia sp. à couvert herbacé continu (Andropogon
majeurs) réalisés sur le bassin versant de la Donga au sp., Hypparhenia sp.). Cette végétation naturelle est
cours de deux années à pluviométrie contrastée aujourd’hui largement dégradée en une mosaı̈que de
(1514 mm en 2003, 1140 mm en 2004). L’analyse champs (igname, manioc, sorgho, maı̈s, coton). À l’aval
des processus en jeu sur le bassin, appuyée par des des versants, les principaux axes de drainage sont bordés
mesures de géophysique de subsurface (résistivité), sur 10 à 20 m par une dense forêt-galerie.
permet d’obtenir une première estimation du bilan L’activité socio-économique, hormis la ville de
hydrologique. Le schéma des flux hydriques proposé Djougou (20 000 habitants) est essentiellement rurale
relativise, à l’échelle de l’étude, l’impact de l’aquifère et basée sur des cultures pluviales (absence d’irrigation).
sur les flux de surface et suggère une importance La densité de la population, proche de 30 hab km 2,
prépondérante des écoulements de surface et de connaı̂t une croissance annuelle de 4,9 % par an (1992–
subsurface comme facteurs de la variabilité interan- 2002) ; l’essentiel des ressources en eau à usage
nuelle des débits observés à l’exutoire. domestique provient d’une multitude de puits villageois,
disséminés sur l’ensemble du bassin. Des suivis de terrain
2. Zone d’étude des puisages quotidiens, effectués en 2003 dans six
villages représentatifs, indiquent une consommation
Le bassin de la rivière Donga est un affluent de domestique de l’ordre de 21 L hab 1 j 1, proche de celle
l’Ouémé, le plus grand fleuve du Bénin (523 km), dont la de 17 L hab 1 j 1 estimée pour l’ensemble du haut
haute vallée couvre une superficie de 14 300 km2 [17]. Il Ouémé [15]. Compte tenu de la population du bassin de la
se situe au Centre-Ouest du Bénin (Fig. 1), sous climat Donga, cela implique des prélèvements anthropiques
tropical de type soudano-guinéen, à saisons contrastées faibles, de l’ordre de 0,2 mm an 1.
(alternance d’une saison sèche et d’une saison plu-
vieuse). Sur ce bassin versant, la pluviométrie moyenne 3. Données
annuelle observée à la station météorologique de
Djougou, sur la période 1950–2004, est de 1195 mm Des données préliminaires sur l’organisation de la
(1280 mm de 1950 à 1969, 1150 mm de 1970 à 2004) ; couverture d’altération ont été obtenues en mai 2003 par
l’essentiel des pluies (60 %) se concentre entre juillet et des prospections géophysiques basées sur la résistivité
septembre. L’évapotranspiration potentielle (ETP) cal- électrique [8]. Dans ce type de milieu, ce paramètre
culée par la méthode de Penman–Monteith [10] est de présente des contrastes importants en fonction de la
1390 mm an 1 (2002–2005) ; à l’échelle mensuelle, elle nature des roches, des teneurs en eau ou en argile. La
montre une faible variabilité saisonnière, avec un méthode de la tomographie électrique [22] a été mise en
maximum de 5 à 6 mm j 1 en saison sèche et un œuvre sur plusieurs versants du bassin, à une échelle
minimum de 3 à 4 mm j 1 en saison des pluies. pluri-hectométrique. En parallèle aux mesures, des
La topographie du bassin est vallonnée, avec des relevés géologiques ont été effectués pour contraindre
altitudes comprises entre 520 m à l’amont (ouest) et l’interprétation.
340 m à l’exutoire (est), avec une pente moyenne de la Les données hydrodynamiques sont issues de 16
rivière de 1,7 m km 1. Les roches sont essentiellement pluviomètres, de 26 puits suivis en piézométrie et de
des formations métamorphiques de la bordure est trois stations limnimétriques ; ces stations représentent,
du craton Ouest-Africain, datées du Protérozoı̈que d’amont en aval, des bassins emboı̂tés, de superficies
 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429 423

croissantes (Ara-pont (A), 13 km2 ; Donga-Kolokondé été réalisées sur des eaux de pluie (18 analyses), de la
(K), 105 km2 ; et Donga-pont BV, 586 km2 ; Fig. 1). En nappe (23 puits et 15 forages), de sources temporaires (8
complément, à plus fine échelle, des données issues analyses) et du cours principal de la Donga (une douzaine
d’un site à l’amont du bassin où neuf piézomètres ont d’analyses par station limnimétrique). Pour les eaux de la
été crépinés à différentes profondeurs (2, 10 et 20 m) Donga, un échantillonnage à haute fréquence a égale-
sont disponibles depuis 2004 (Fig. 1). Les pluviomètres ment pu être effectué sur trois crues majeures en 2003.
enregistreurs permettent de discriminer les événements L’ensemble des échantillons, filtrés in situ, a été analysé
par cumuls de 0,5 mm. Les niveaux de la nappe libre ont avec une précision de 5 % par électrophorèse capillaire
été enregistrés à des pas de temps de 30 min (13 sites) (cations et SO42 , NO3 , Cl ) ou titration acide
et/ou par lectures tri-journalières (23 sites). Les trois (HCO3 ) au laboratoire HydroSciences de Montpellier
stations limnimétriques considérées, étalonnées, ont (France).
bénéficié d’enregistrements au pas de temps de 5 min
pour les deux stations à l’amont (A, K) et de 30 min 4. Résultats
pour la station à l’aval (exutoire). L’ensemble des
piézomètres et des échelles limnimétriques ont été 4.1. Géophysique de subsurface
nivelés par GPS différentiel pour caractériser la
distribution spatiale des potentiels hydrauliques [16]. Plusieurs profils de résistivité électrique transversaux
Les données hydrochimiques sont constituées de aux talwegs ont été réalisés par tomographie en différents
suivis mensuels à trimestriels des paramètres physico- points du bassin [8,16]. Un exemple typique de profil de
chimiques – température, pH, conductivité électrique résistivité, obtenu en amont du bassin, est présenté sur la
(CE, à 25 8C) – sur une vingtaine de puits et forages. Des Fig. 2 ; ce profil montre une hétérogénéité importante
mesures ont également été effectuées sur le cumul d’une de l’altération à une échelle hectométrique. Les
dizaine d’événements pluvieux sur cinq pluviomètres de formations argileuses et/ou argilo-sableuses électrique-
l’amont du bassin, ainsi que dans les eaux de surface, en ment conductrices sont organisées en lentilles parfois
divers points du bassin (sources, affluents et cours affleurantes ; cette distribution suggère que des nappes
principal de la Donga). La station limnimétrique de perchées peuvent se mettre en place, les formations
Donga-Kolokondé (K, Fig. 1) fait l’objet depuis juin argileuses constituant des niveaux imperméables. Des
2003 d’un suivi automatique de la CE à des pas de temps cuirasses superficielles se notent par de fortes résistivités
de 30 min ; sur les piézomètres du site expérimental de électriques. Le socle, électriquement résistant, se
l’amont du bassin, des enregistrements de la CE à pas de rencontre entre 15 et 25 m sous la surface. Son modelé
temps de 30 min sont également effectués depuis avril présente de faibles ondulations sans relation avec la
2004. En parallèle à ces mesures, des analyses de la topographie. Localement, des formations massives de
chimie des ions majeurs (Ca2+, Mg2+, Na+, K+ pour les plus forte résistance (par exemple, des quartzites),
cations, HCO3 , SO42 , NO3 , Cl pour les anions) ont propices au ruissellement, affleurent (Fig. 2).

Fig. 2. Coupe de distribution des résistivités sur un versant représentatif du sous-bassin d’Ara, mai 2003 (localisée à « PG », Fig. 1). L’équivalence
entre résistivité et type de terrain est proposée sur la base de mesures de résistivité sur affleurement.
Fig. 2. Resistivity cross-section in the Ara watershed in May 2003 (located as ‘‘PG’’ in Fig. 1). The geological interpretation is based on resistivity
measurements made on different outcropping surfaces.
424 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429

respectivement de 28,3 et 13,7 % (Tableau 1). Les

mêmes conclusions peuvent être apportées pour les
sous-bassins emboı̂tés A et K (Tableau 1). En terme de
processus, l’écoulement présente une dynamique
rapide, avec des crues journalières pouvant atteindre
140 m3 s 1 (Fig. 3B). À l’échelle saisonnière, la durée
de l’écoulement apparaı̂t plus courte en 2004 qu’en
2003 (120 contre 210 jours), tributaire de l’occurrence
de pluies rapprochées.
La profondeur de la nappe est comprise entre 15 m
(saison sèche) et moins de 1 m sous le sol (saison des
pluies), en relative indépendance par rapport à la
position dans le bassin. À cette échelle, la piézométrie
de la nappe est toujours en position haute par rapport au
lit de la rivière. Cette observation est également valable
à l’échelle du sous-bassin versant d’Ara, plus densé-
ment équipé. Sur ce site, les piézomètres crépinés à 10
Fig. 3. Lames journalières précipitées (A), lames journalières écou- et 20 m montrent une piézométrie identique (nappe
lées à l’exutoire du bassin (B) et profondeur moyenne journalière de la phréatique pérenne des altérites), tandis que ceux
nappe (C) pour les années 2003 et 2004.
crépinés à 2 m captent un niveau saturé à potentiel plus
Fig. 3. Daily rainfall (A), daily runoff at the watershed outlet (B), and élevé (nappe perchée temporaire). Le maximum
mean daily water table depth (C) for years 2003 and 2004.
piézométrique de la nappe phréatique est atteint en
août–septembre, après une hausse plurimétrique des
4.2. Hydrodynamique niveaux, sans répartition particulière des amplitudes de
fluctuation sur le bassin. De même que les écoulements
Pour les deux années suivies, les pluies présentent de surface, la décrue piézométrique s’amorce dès la fin
une saisonnalité similaire, avec un pic mensuel en des précipitations régulières, pour atteindre un mini-
juillet–août (Fig. 3A). La pluviométrie présente une mum en mars–avril (Fig. 3C). Ces fluctuations piézo-
diminution du cumul de 25 % (1514 mm en 2003, métriques montrent, en comparaison interannuelle, une
1140 mm en 2004, Tableau 1). Cette différence se différence moindre que celle des débits (Fig. 3B) :
trouve accentuée (64 %) dans les écoulements de l’amplitude des chroniques varie de 3,3 à 12,5 m
surface à l’exutoire (429 mm en 2003, 156 mm en (moyenne : 7,4  2,8 m) en 2003 et de 2,4 à 10,0 m
2004), avec des coefficients d’écoulement à l’exutoire (moyenne : 5,3  2,3 m) en 2004. La nappe phréatique

Tableau 1
Bilan hydrologique pour les deux années de suivi à l’échelle du bassin versant (BV) ; la pluie (P) et le débit (Q) sont également calculés pour les deux
autres exutoires intermédiaires (A et K, cf. Fig. 1)
Table 1
Water balance for the two surveyed years at the basin scale (BV); rainfall (P) and outflow (Q) are also reported for the two intermediate outlets (A and
K, cf. Fig. 1)
P (mm) Q (mm), [Q/P] (%) R (mm), [R/P]% ETR (mm)
A K BV A K BV BV, n = 1 % BV, n = 5 % BV
2003 1634 1580 1514 480 437 429 74 370 1085
[29.4] [27.7] [28.3] [4.9] [24.4]
2004 1066 1095 1140 158 106 156 53 265 984
[14.8] [9.6] [13.7] [4.6] [23.2]
La recharge de la nappe (R) est estimée par la méthode des fluctuations piézométriques pour deux valeurs possibles de la porosité efficace de
l’aquifère (n = 1 %, n = 5 %) ; l’évapotranspiration réelle (ETR) est déduite par soustraction aux termes du bilan externe du bassin (P Q), en
négligeant le stockage interannuel dans l’aquifère (absence de tendance piézométrique) et les prélèvements anthropiques (0,2 mm an 1).
Groundwater recharge (R) is estimated by the water table fluctuation method for two possible values of the specific yield (n = 1 %, n = 5 %); effective
evapotranspiration (ETR) is deducted from the external terms of the balance (P Q), considering no groundwater storage at the interannual scale
and neglecting the human water consumption (0.2 mm yr 1).
 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429 425

Fig. 4. Suivis des fluctuations de la conductivité électrique (CE) des

eaux.
Fig. 4. Seasonal EC surveys of surface and groundwaters.

des altérites apparaı̂t donc de dynamique synchrone à

l’échelle du bassin, avec une continuité hydraulique
effective à l’échelle des versants. Fig. 5. Fluctuations typiques de la chimie des eaux mesurées à
l’exutoire au cours de la crue composite du 1er août 2003. A :
4.3. Hydrochimie concentrations (mequiv L 1), B : flux (equiv s 1).
Fig. 5. Typical changes in surface water chemistry during the com-
Les eaux échantillonnées sont peu minéralisées posite flood of 1st August 2003 at the watershed outlet. A: ionic
(hormis quelques sites pollués enrichis en NO3 et Cl), content (mequiv L 1), B: flux (equiv s 1).
avec des CE généralement inférieures à 400 mS cm 1
(Fig. 4). Les eaux sont à tendance acide (pH 5–7) et de suivies, aucune tendance à l’augmentation de la CE des
faciès chimique majoritairement bicarbonaté calco- eaux de la rivière n’est observée (Fig. 4).
magnésien, en accord avec la nature silicatée des roches. La nappe phréatique (altérites) montre des valeurs
La température des eaux de la nappe (27 à 28 8C) est nettement supérieures, en moyenne de 155  70 mS
légèrement supérieure à la température moyenne cm 1 ; plus en profondeur, les eaux des forages (socle)
interannuelle à la station de Djougou (26,4 8C). montrent des valeurs encore plus élevées (288 
Les pluies montrent des CE faibles, de l’ordre de 80 mS cm 1) ; cela témoigne d’un temps de résidence
10 mS cm 1. Ces valeurs sont logiquement inférieures à probablement plus important. Sur le site de plus fine
celles des eaux de surface (Donga et affluents) dont la investigation, à l’amont du bassin, les suivis hydro-
CE, hormis les tout premiers événements lessivants de chimiques effectués en 2004 permettent de préciser ces
la saison, reste faible en valeur absolue (<70 mS cm 1). observations. La nappe perchée (2 m) montre des CE
En valeur relative, la CE montre des fluctuations fluctuant entre 30 et 130 mS cm 1 (moyenne de
systématiques en cours de crue : à une chute de CE en 75  30 mS cm 1), alors que celles relatives aux piézo-
début d’hydrogramme, succède un recouvrement lent mètres crépinés à 10 et 20 m sont respectivement de
(quelques jours) de la CE initiale (la Fig. 5A représente l’ordre de 200 et 250 mS cm 1 (Fig. 4).
une crue composite à l’exutoire). En termes de flux En termes de minéralisation, les rapports ioniques
ioniques, la dynamique géochimique montre des montrent une grande constance dans le temps et
variations similaires, mais retardées de quelques heures l’espace. Pour les eaux de la rivière (Fig. 5A), les
par rapport au débit (Fig. 5B) ; cette diachronie peut éléments majeurs liés à l’altération (Si, Ca, Mg, Na et
s’expliquer par des contributions variables de sous- HCO3) fluctuent proportionnellement à la CE ; l’ion K
bassins aux taux d’érosion chimique différents. fait exception et montre un comportement variable,
D’amont en aval du bassin, en période d’inter-crue, vraisemblablement en raison d’une sensibilité plus
une faible augmentation de la CE est perceptible (de grande aux cycles biogéochimiques [4]. Au final, les
39  5 à 65  7 mS cm 1 pour A et BV) ; l’explication eaux issues de différentes profondeurs (subsurface,
la plus probable réside dans une augmentation de la nappe des altérites et socle profond) montrent des
longueur moyenne des versants et donc du temps de rapports ioniques peu différenciés, reflet de l’homo-
transit de l’écoulement de subsurface. À l’échelle généité d’ensemble des faciès chimiques des eaux du
saisonnière comme à l’interannuel pour les deux années bassin.
426 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429

5. Discussion avec des débits spécifiques de l’ordre de 1 m3 h 1 m 1

et des perméabilités de l’ordre de 10 5 m s 1 [16].
5.1. Processus hydrologiques Le caractère temporaire des eaux de la Donga
(absence de « débit de base ») traduit une origine
Le croisement d’informations de la géophysique, de superficielle des écoulements. Lors des pics de crue, les
l’hydrodynamique et de la géochimie permet de eaux présentent une minéralisation très faible, proche de
reconstituer un schéma cohérent des processus en jeu celle des pluies ; en période d’inter-crue, les écoulements
sur le bassin (Fig. 6). sont de faible minéralisation, de chimisme voisin de celui
La nappe phréatique, libre, est, sur tout le bassin, en des nappes perchées saisonnières, tel qu’observé à leurs
position haute par rapport au réseau hydrographique. émergences en de multiples points du bassin [16].
Elle occupe les versants en continuité hydraulique, L’écoulement à l’exutoire peut donc être considéré
comme en témoignent les mesures piézométriques. Les comme composé d’un flux rapide (ruissellement
coupes géophysiques montrent qu’elle baigne différents Hortonien sur les surfaces indurées et écoulement sur
terrains organisés en lentilles silto-argileuses ou surfaces saturées) et d’un flux retardé de subsurface
sableuses. Au cours de la saison des pluies, elle montre (vidange des nappes perchées), sans contribution
en tout point des fluctuations saisonnières, d’amplitude significative de la nappe des altérites à l’écoulement.
variable et sans organisation spatiale particulière ; la La vidange de la nappe semble donc s’effectuer
recharge apparaı̂t donc principalement de type direct, essentiellement par des transferts verticaux par évapo-
par percolation des eaux de pluie et infiltration à travers transpiration (strate arborée) et éventuellement, vers des
la zone non saturée ; en certains points du bassin fractures du socle plus en profondeur, via des brèches
cependant, l’existence de sols indurés (cuirasse) et/ou tectoniques à forte porosité [1] ; cependant, à une échelle
d’affleurements du socle peut induire un ruissellement de 10 000 km2, le haut Ouémé présente un tarissement
localisé avant infiltration. La nature silto-argileuse à annuel, ce qui milite plutôt pour des circulations dans le
sablo-argileuse du manteau altéritique milite pour un socle faibles, ou d’extension limitée.
rôle essentiellement capacitif du réservoir ; cette
interprétation a été confirmée par essais de pompage, 5.2. Bilan hydrologique

Deux méthodes sont appliquées pour contraindre le

bilan hydrologique du bassin. La première est basée sur
les fluctuations piézométriques de la nappe libre et a pour
objet une estimation de la recharge annuelle. La deuxième
est basée sur une décomposition des hydrogrammes de
crue et a pour objet de discriminer la part des écoule-
ments rapides et de subsurface dans les flux à l’exutoire.
L’estimation de la recharge par la méthode des
fluctuations piézométriques suppose une nappe libre,
rechargée par les précipitations en un temps court [7] ;
elle peut s’exprimer selon l’équation R = h  n, où R
représente la recharge (L T 1), h l’amplitude piézo-
métrique observée sur la période considérée (L T 1) et n
Fig. 6. Schéma conceptuel des flux hydriques à l’échelle d’un versant, la porosité efficace (%, adimensionnel). Pour notre cas
d’après les observations hydrologiques, géochimiques et géophysi- d’étude où la recharge s’effectue en quelques mois,
ques réalisées sur le bassin de la Donga. P : Pluie, E : évapotranspira-
l’amplitude piézométrique considérée est la valeur
tion, If : infiltration, R : recharge, RR : ruissellement rapide
(ruissellement Hortonien et sur surface saturée), SS : flux de subsur- moyenne entre le minimum de saison sèche et le
face, Fp : flux profond. Les niveaux saisonniers minimum et maximum maximum de la saison des pluies (Fig. 3C). L’estimation
de la nappe phréatique sont indiqués en tiretés. obtenue pour une valeur n considérée représente alors, du
Fig. 6. Conceptual cross-section of the water flux at a hillslope scale, fait (1) d’une probable vidange concomitante à la
according to the hydrological, geochemical, and geophysical observa- recharge et (2) d’une amplitude de fluctuation parfois
tions made in the Donga catchment. P: rainfall, E: evapotranspiration, limitée par le niveau du sol, une estimation a minima de la
If: infiltration, R: recharge, RR, runoff (Hortonian runoff and runoff on
saturated surface), SS: subsurface flow, Fp: deep groundwater flow.
recharge effective. Les prélèvements villageois, suscep-
Dotted lines represent minimum and maximum seasonal levels of the tibles d’induire un rabattement local de la nappe, ont
water table. probablement une influence minime, les fluctuations
 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429 427

mesurées dans les piézomètres expérimentaux présentant observations réalisées à la station « A » montrent des CE
des amplitudes similaires à celles obtenues dans les puits. proches de celles de l’eau de pluie durant les pics de
L’influence des flux souterrains provenant de l’amont des crue (8 < CE < 17 mS cm 1), ce qui peut justifier de
points de mesure est considérée comme négligeable, les l’assimilation de la CE des pluies (10  2 mS cm 1) au
observations réalisées sur des transects piézométriques le pôle RR. En considérant les valeurs de CE mesurées aux
long de versants montrant de faibles différences de trois stations en période d’inter-crue, une valeur de
fluctuation saisonnière (<15 %). Pour quantifier l’impact 55  6 mS cm 1 a été fixée au pôle SS ; cette valeur est
du choix de n sur la recharge estimée, une gamme de 1 à dans la gamme de celles observées pour les nappes
5 % a été retenue pour la porosité efficace (Tableau 1), en perchées (75  30 mS cm 1). Les calculs se rapportent à
accord avec les valeurs généralement admises pour les la station K (Fig. 1), équipée d’un enregistreur de CE sur
altérites de socle en Afrique de l’Ouest [9,13]. La les deux années considérées, aux valeurs de CE proches
recharge déduite est de 74 à 370 mm an 1 (2003) et de 53 de celles de l’exutoire (53  8 et 65  7 mS cm 1
à 265 mm an 1 (2004), soit de 5 à 24 % du total respectivement, en période d’inter-crue). Pour tenir
pluviométrique. Ces coefficients de recharge sont compte de l’incertitude sur les pôles, une simulation a
compatibles avec la gamme de 20 à 30 % obtenue pour été réalisée par la méthode de Monte Carlo (tirage des
des aquifères en zone de socle et sous climat soudanien valeurs des pôles dans des lois normales de moyenne et
[13]. À l’échelle du bassin, ces chiffres apparaissent écart-type précédemment donnés). Les flux de subsurface
compatibles avec des flux d’exfiltration par évapotrans- représentent 68  13 % et 83  13 % de l’écoulement
piration pendant les cinq à six mois de la saison sèche. À total, respectivement en 2003 et 2004 ; à l’échelle intra-
l’échelle des versants (Fig. 6), une estimation des flux saisonnière, aucune évolution significative des contribu-
souterrains peut être effectuée par l’application de la loi tions n’est observée. Sur cette partie du bassin versant, les
de Darcy : sur la base des paramètres de l’aquifère sommes des écoulements rapides et de subsurface
(perméabilité de 10 5 m s 1, épaisseur saturée de 10 à peuvent donc être estimées respectivement à 140 et
20 m, gradients hydrauliques de l’ordre de 2 à 3 %), les 297 mm en 2003, et à 18 et 88 mm en 2004. Pour ces deux
flux parvenant aux talwegs sont de l’ordre de années à écoulement contrasté, ces calculs suggèrent que
0,3 m3 m 1 j 1. En considérant une largeur de forêts- la baisse des écoulements n’a pas une origine unique,
galeries de 20 m, un flux d’évapotranspiration de l’ordre mais résulte d’un déficit combiné des flux de surface et de
de 15 à 20 mm j 1 serait alors suffisant pour empêcher subsurface. Ce processus s’avère donc distinct de celui
toute contribution du souterrain au réseau hydrographi- observé plus au nord au Sahel, où le déficit pluviomé-
que. Selon [10], en raison de l’effet d’oasis et de leur trique combiné au déboisement a provoqué la seule
organisation linéaire, l’ETP des forêts galeries peut être augmentation du ruissellement Hortonien sur surfaces
d’au moins un facteur 2 supérieure à l’ETP de 5 à encroûtées (par exemple, [26]) ; en zone soudanienne, un
6 mm j 1 en saison sèche, calculée selon les seuls sol plus développé favorise l’infiltration locale et les flux
paramètres atmosphériques. En première approche, de subsurface, dont la réponse au déficit pluviométrique
l’ordre de grandeur des flux souterrains apparaı̂t donc peut être différente, en amplitude, de celle du ruisselle-
compatible avec des flux d’exfiltration gouvernés par ment.
évapotranspiration ; des études complémentaires seraient
cependant nécessaires pour mieux quantifier les flux 6. Conclusion
hydriques à l’aval des versants.
En accord avec les processus identifiés, une décom- Cette étude, basée sur une approche couplée
position de l’hydrogramme en deux écoulements hydrodynamique et hydrochimique des flux hydriques,
(ruissellement, RR, et flux de subsurface, SS) a été montre que l’origine des écoulements sur le bassin
retenue. Les hypothèses classiques [2] suivantes ont été versant de la Donga est à dominante superficielle. Ce
adoptées : (1) on a invariance spatiotemporelle de la résultat précise et complète, pour le haut Ouémé, des
signature des réservoirs contributeurs ; (2) la CE études antérieures sur un sous-bassin moins étendu
représente une valeur représentative de la charge ionique (30 km2), où l’impact de la mise en culture a engendré
issue de l’altération. Pour notre cas d’étude, des variations une augmentation du ruissellement et une réduction de
temporelles du chimisme sont cependant notées en l’évapotranspiration [14], en accord avec la tendance
subsurface (Fig. 4) et K+ évolue indépendamment de la observée à l’échelle globale [25]. À terme, par une
CE (Fig. 5A) ; ces écarts restent suffisamment mineurs meilleure prise en compte des processus hydrologiques
(par exemple, R2 de 0,91 entre CE et Ca, Fig. 5A) pour et de la variabilité interannuelle des flux, une
justifier, en première approche, les hypothèses. Des amélioration des modèles hydrologiques établis pour
428 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429

le haut Ouémé [3,20,27] est attendue ; en retour, ces [5] BRGM/OBM, Carte géologique au 1/200 000e, feuille de Djou-
modèles permettront de valider la cohérence des gou–Parakou–Nikki, Office béninois des mines, Cotonou, répu-
blique du Bénin, 1984.
gammes de valeurs estimées pour les différents [6] P. Chevallier, O. Planchon, Hydrological processes in a small
paramètres du bilan hydrologique du bassin. humid savannah basin (Ivory Coast), J. Hydrol. 151 (1993)
Ces résultats constituent une première étape dans 173–191.
l’étude hydrologique de ce bassin versant. Dans le cadre [7] R.S. Crosbie, P. Binning, J.D. Kalma, A time series approach to
inferring groundwater recharge using the water table fluctuation
du programme AMMA, un suivi à plus long terme
method, Water Resour. Res. 41 (2005), doi:10.1029/
permettra de préciser encore davantage l’importance 2004WR003077.
respective du ruissellement et des écoulements de [8] M. Descloitres, M. Wubda, Y. Le Troquer, Prospections géo-
subsurface dans la variabilité interannuelle des débits. physiques sur le bassin versant d’Ara : électrique 2D et électro-
Dans une perspective plus large, cette étude représente magnétisme EM34 ; Compte- rendu de mission du 5 au 14 mai
également une contribution à une meilleure connais- 2003, Rapport technique IRD, Ouagadougou, Burkina-Faso,
2003 (15 p.).
sance de la variabilité spatiotemporelle des flux surface- [9] M. Engalenc, Méthode d’étude et de recherche de l’eau souter-
souterrain et de son impact sur l’hydrologie de l’Afrique raine des roches cristallines de l’Afrique de l’Ouest, CIEH,
tropicale [18,23]. Maisons-Alfort, France, 1978 (318 p.).
[10] FAO, Crop Evapotranspiration – Guidelines for computing crop
Remerciements water requirements. FAO Irrigation and Drainage paper, 56,
Rome, 1998. Available online at: http://www.fao.org/docrep/
X0490E/X0490E00.htm.
Une part importante de l’acquisition de données de [11] P. Faure, Carte pédologique de reconnaissance au 1:200 000e,
terrain a été assurée par S. Afouda, M. Arjounin, J.-M. feuille de Djougou, Orstom Éd., Paris, France, 1977.
Bouchez, F. Jacquin, F. Malinur et T. Ouani, techniciens [12] P. Faure, B. Volkoff, Some factors affecting regional differentia-
tion of the soils in the Republic of Benin (West Africa), Catena
de l’IRD au Bénin. Les traitements des données
32 (1998) 281–306.
pluviométriques et hydrométriques ont été coordonnés [13] C. Filippi, F. Milville, D. Thiery, Évaluation de la recharge
respectivement par L. Le Barbé et C. Peugeot (IRD). La naturelle des aquifères en climat soudano-sahélien par modéli-
direction de l’hydraulique de Djougou (MMEH, Bénin) sation hydrologique globale : application à 10 sites au Burkina-
a permis l’accès aux données et aux infrastructures Faso, J. Sci. Hydrol. 35 (1990) 29–48.
hydrauliques du bassin. Les villageois de la zone [14] S. Giertz, B. Junge, B. Diekkrüger, Assessing the effects of land
use change on soil physical properties and hydrological pro-
étudiée sont remerciés pour leur accueil et leur cesses in the sub-humid tropical environment of West Africa,
participation active à certains suivis. Ce travail a été Phys. Chem. Earth. 30 (2005) 485–496.
rendu possible par une bourse de thèse de la Côte [15] K. Hadjer, T. Klein, M. Schopp, Water consumption embedded
d’Ivoire (B. Kamagaté), par les financements institu- in its social context, north-western Benin, Phys. Chem. Earth 30
(2005) 357–364.
tionnels des programmes français ECCO-PNRH (2003–
[16] B. Kamagaté, Fonctionnement hydrologique et origine des
2006) et ORE (2002–2010), ainsi que par le soutien du écoulements sur un bassin versant en milieu tropical de socle
projet européen AMMA-EU (2005–2009). au Bénin : bassin versant de la Donga (haute vallée de l’Ouémé),
thèse, université Montpellier-2, France, 2006 (319 p.).
Références [17] L. Le Barbé, G. Alé, B. Millet, H. Texier, Y. Borel, R. Gualde,
Les ressources en eaux superficielles de la république du Bénin,
IRD Éditions, collection « Monographies hydrologiques », no 11,
[1] P. Affaton, Le bassin des Voltas (Afrique de l’Ouest) : une marge
passive, d’âge Protérozoı̈que supérieur, tectonisée au Panafricain Paris, 1993 (469 p.).
(600  50 Ma), thèse, université Aix-Marseille-3, 1987, Orstom [18] T. Lebel, T. Vischel, Climat et cycle de l’eau en zone tropicale :
Éditions, collection « Études et thèses », Paris, 1990 (499 p.). un problème d’échelle, C. R. Geoscience 337 (2005) 29–38.
[19] M. Le Lay, S. Galle, Variabilité interannuelle et intra-saison-
[2] J.-L. Boeglin, Y. Tardy, Érosion chimique et mécanique sur le
bassin amont du Niger (Guinée, Mali). Découpage en quatre nière des pluies aux échelles hydrologiques. La mousson ouest-
écoulements, C. R. Acad. Sci. Paris, Ser. IIa 325 (1997) 125–131. africaine en climat soudanien, J. Sci. Hydrol. 50 (2005)
509–524.
[3] H. Bormann, T. Fass, S. Giertz, B. Junge, B. Diekkrüger, B.
Reichert, A. Skowronek, From local hydrological process analysis [20] M. Le Lay, S. Galle, How changing rainfall regimes may affect
to regional hydrological model application in Benin: concept, the water balance: a modelling approach in West Africa, in :
results and perspectives, Phys. Chem. Earth 30 (2005) 347–356. Regional Hydrological Impacts of Climatic Change – Hydro-
climatic Variability AISH 296, 2005, 203–210.
[4] J.-J. Braun, B. Dupré, J. Viers, J.-R. Ndam Ngoupayou, J.-P.
Bedimo Bedimo, L. Sigha-Nkamdjou, R. Freydier, H. Robain, B. [21] Y. L’Hôte, G. Mahé, B. Some, J.-P. Triboulet, Analysis of a
Nyeck, J. Bodin, P. Oliva, J.-L. Boeglin, S. Stemmler, J. Berthe- Sahelian annual rainfall index from 1896 to 2000; the drought
lin, Biogeohydrodynamic in the forested humid tropical envir- continues, J. Sci. Hydrol. 47 (2002) 563–572.
[22] M.H. Loke, Electrical imaging surveys for environmental and
onment: the case study of the Nsimi small experimental
watershed (south Cameroon), Bull. Soc. geol. France 173 engineering studies. Technical note of Geotomo Software, 2000.
(2002) 347–357. Available online at: http://www.geoelectrical.com.
 B. Kamagaté et al. / C. R. Geoscience 339 (2007) 418–429 429

[23] G. Mahé, J.-C. Olivry, Évaluation des apports en eau douce à [26] L. Seguis, B. Cappelaere, G. Milesi, C. Peugeot, S. Massuel, G.
l’Atlantique depuis les côtes de l’Afrique intertropicale (1951– Favreau, Simulated impacts of climate change and land-clearing
1989), C. R. Acad. Sci. Paris, Ser. IIa 328 (1999) 621–626. on runoff from a small Sahelian catchment, Hydrol. Process. 18
[24] G. Mahé, J.-C. Olivry, R. Dessouassi, D. Orange, F. Bamba, E. (2004) 3401–3413.
Servat, Relation eaux de surface – eaux souterraines d’une [27] N. Varado, I. Braud, S. Galle, M. Le Lay, L. Seguis, B.
rivière tropicale au Mali, C. R. Acad. Sci. Paris, Ser. IIa 330 Kamagate, C. Depraetere, Multi-criteria assessment of the
(2000) 689–692. Representative Elementary Watershed approach on the
[25] A. Perrier, A. Tuzet, Le cycle de l’eau et les activités au sein de Donga catchment (Benin) using a downward approach
l’espace rural. Enjeux globaux, solutions locales et régionales, of model complexity, Hydrol. Earth Syst. Sci. 10 (2006)
C. R. Geoscience 337 (2005) 39–56. 427–442.