C. R.

Geoscience 336 (2004) 1099–1109

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

Comportement d’une source karstique soumise à une sécheresse

prolongée : la source Bittit (Maroc)
Fouad Amraoui a , Moumtaz Razack b,∗ , Lhoussaine Bouchaou c
a Laboratoire d’hydrogéologie, faculté des sciences Ain Chock, université Hassan-II, BP 5366, Maarif, Casablanca, Maroc
b Département d’hydrogéologie, UMR 6532, université de Poitiers, 40, av. du Recteur-Pineau, 86022 Poitiers cedex, France
c Laboratory de géologie appliquée et de géo-environnement, faculté des sciences Ibn-Zohr, BP 28/5, 80000 Agadir, Maroc

Reçu le 19 mai 2003 ; accepté après révision le 29 mars 2004

Présenté par Ghislain de Marsily

Résumé
Le comportement d’une émergence karstique (source Bittit, Moyen Atlas marocain) soumise à une sécheresse prolongée
(> 20 ans), est analysé par des méthodes d’étude des hydrogrammes et des méthodes corrélatoires et spectrales. Deux séries
de débit, correspondant à une période pluviométrique normale et à une période de sécheresse, sont analysées. Les résultats
montrent que le système karstique est inertiel, avec des réserves très importantes. Cependant malgré le grand pouvoir régulateur
des réserves, la longue période de sécheresse a eu pour conséquence une diminution significative des débits, du fait d’apports
moins diversifiés et d’une réduction des réserves. Pour citer cet article : F. Amraoui et al., C. R. Geoscience 336 (2004).
 2004 Académie des sciences. Publié par Elsevier SAS. Tous droits réservés.
Abstract
Behaviour of a karstic spring subjected to a long drought period: Bittit spring (Morocco). The behaviour of a karstic
spring (Bittit spring, Middle Atlas, Morocco) subjected to a prolonged drought (> 20 yr), is analysed using spring hydrographs
analysis methods and correlation and spectral methods. Two series of discharges, corresponding to a normal rainfall period and
a drought period, are analysed. The results show that the karstic system is inertial, with very significant storage. However, in
spite of the great regulating capacity of the storage, the long drought period significantly reduced the spring discharges, due to
less diversified inputs, and to the karstic system storage depletion. To cite this article: F. Amraoui et al., C. R. Geoscience 336
(2004).
 2004 Académie des sciences. Publié par Elsevier SAS. Tous droits réservés.

Mots-clés : source karstique ; Moyen Atlas ; Maroc ; sécheresse ; hydrogrammes ; analyses corrélatoire et spectrale

Keywords: karstic spring; Middle Atlas; Morocco; drought; hydrographs; correlation and spectral analysis

* Auteur correspondant.
Adresse e-mail : [email protected] (M. Razack).

1631-0713/$ – see front matter  2004 Académie des sciences. Publié par Elsevier SAS. Tous droits réservés.
doi:10.1016/j.crte.2004.03.016
1100 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109

Abridged English version calculated over the period 1975–2001) and a reduced
dispersion (standard deviation = 186 l s−1 ; CV =
1. Introduction 12.3%). The spring thus presents a relatively stable
discharge, indicating a good inter-annual regulation
The object of this article is, from a general point by the karstic system. Fig. 3 illustrates the discharge
of view, to characterize the hydrodynamic behaviour evolution, which shows a regular decrease since the
of large karstic systems according to their storage beginning of the drought. The discharge time series of
level. The study site is represented by the Bittit spring, the Bittit spring, used in this work, cover respectively
located in the Middle Atlas, in Morocco. The reason the periods 1975–1977 and 1995–1998.
for this choice is the large size of this system and its
complexity, the existence of a long series of rainfall 2. Methods
deficit and the availability of more than 70 yr of
data. A better understanding of the evolution of such Three methods have been used in this work to
large karstic systems, taking into account their storage analyse the influence of the long drought period on
level, is particularly important within the context the spring hydrodynamic behaviour. These methods
of climatic changes, evolving in the studied region are the discharge distribution analysis, the recession
towards aridity. curves analysis and correlation and spectral analysis.
The plain of Sais is one of the great agricultural Other methods, like the use of natural tracers [4,7] will
regions of Morocco. Its water resources are provided be applied shortly.
mainly by the karstic aquifer of the Liasic carbon- The first method deals with the distribution of the
ate plateau of Middle Atlas (Fig. 1). This plateau is discharge cumulative frequencies [11] and allows us
drained by a group of springs of which the most signif- to check (or not) the unity of the karstic systems.
icant is the Bittit spring (average annual discharge = Thus, one can highlight sinks (or sources) out of
1.3 m3 s−1 ) [1–3,5]. This spring is currently exploited (or into) the studied karstic system. These features
to partly supply the town of Meknès with drinking wa- appear as discontinuities on the frequency curve. The
ter. The Middle Atlas area has undergone since two fitted probability distribution function is Gaussian.
decades a regular rainfall deficit (Fig. 2A) with signif- In practice, one plots the discharge frequencies on a
icant consequences on water resources (Fig. 2B). Gaussian functional diagram. The observation on the
The Middle Atlas plateau is made of limestone diagram of several segments displays the existence of
of Liasic age. Several springs are located along the several flow regimes. These changes of flow regimes
contact between the carbonate plateau and the plain can have internal origins related to hydrodynamic
of Sais. The Middle Atlas plateau is the recharge area modifications within the karstic system (storage and/or
of the karstic system, which is drained by the springs. release of water; overflows), or external origins related
Fig. 2A represents the rainfall variations compared to relationships with nearby systems resulting in water
to the annual average rainfall, from 1934 to 2000. It contributions or losses.
confirms the irregularity of the rains and highlights The recession curves analysis [11] allows us to
a significant reduction since 1980. The total rainfall characterize the behaviour of two subsystems: the un-
deficit since 1980 at the El Hajeb station is equivalent saturated zone or infiltration subsystem and the sat-
to three years of average rainfall. urated zone or drowned karst subsystem. The flood
The hydrogeologic system of Bittit has a very discharge shows how the infiltration zone is drained,
significant storage capacity, allowing storage of a whereas the base flow discharge characterizes the
significant recharge from rainfall and snow melting, draining of the drowned karst. The interest of this
without causing brutal floods. The spring discharge is method is to define the characteristics of each subsys-
indeed quite constant all the year around. The presence tem using some simple parameters.
of sandy dolomites at the base of the Liasic aquifer Correlation and spectral methods are inspired by
suggests the existence of an aquifer with interstitial signal processing methods. Their application to karstic
porosity and relatively slow flow. The Bittit spring is systems proved efficient to identify the organization
characterized by a high average discharge (1510 l s−1 and the dynamics of karstic systems [6,8,9,12,13].
 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109 1101

The correlogram highlights the dependence between In the high-discharge domain, a discontinuity with
each other of successive events for increasing time a slope increase (segments 2 and 3, period 1975–1977)
intervals. Analytically, it shows the evolution of the means water outlets outside the system and can sug-
autocorrelation coefficients rk . If the correlogram gest the probable existence of overflows (for example,
decreases quickly, it indicates events independent a drain above the principal spring). In a dry period,
from each other. On the other hand, a slow decrease one observes a discontinuity with reduction in slope
characterizes events dependent on each other. The between segments 2 and 3. This characteristic means
value rk = 0.2 is a remarkable value, as the events the existence of a different water contribution to the
can be considered independent when rk < 0.2. The system and can represent a temporary retention of a
spectral analysis corresponds to a change from a time certain volume of water in additional structures at the
domain to a frequency domain by a Fourier transform time of the flood and its release at the beginning of the
of the correlogram. The spectral function S(f ) of recession. The discontinuity between segments 3 and 4
density represents the decomposition of the time of period 1975–1977 suggests the contribution of an
series variance in the frequency domain. Interpretation external water volume only at the time of the flood
of this function consists in identifying the various periods (modification of the basin limits, snow melt-
peaks representing periodic phenomena, which helps ing, rain). During the period 1975–1977, the spring
to characterize the system. The weaker the cut-off discharge was higher than at present and it was dis-
frequency (the frequency such as S(f ) is almost null), tributed in four segments. Recharge of the spring was
the more inertial the system. The S(f ) function also more constant and probably more diversified before
allows the calculation of the regulation time of the the years 1980 than during the current drought period.
system, which represents the influence of the duration Indeed, if the groundwater level drops because of a re-
of the input signal [Treg = S(F = 0)/2]. The longer duction of the contributions, secondary drains become
the regulation time, the more inertial the system. unsaturated and the recharge area is reduced. This can
be seen on the frequency diagram by the reduction of
the number of flow regimes (three segments instead of
3. Results and discussion
four). The modifications of the recharge area consid-
ered above are only assumptions, which could not be
Frequency distributions of the discharge were ana- checked as no tracing tests [3] were performed in the
lysed for two periods. The first corresponds to two system.
hydrological cycles with normal rainfall 09/1975– The choice of the various floods that were used
08/1977 and the second with two hydrological cycles for the recession analysis was made in such a way
located in a drought period 09/1995–08/1997. The that they represent the maximum amplitudes of each
diagrams (Fig. 4) show three (period 1975–1977) hydrological cycle and they precede long periods of
to four (period 1995–1997) segments. This feature low water level (Fig. 5). The recession coefficients
reveals that discharge evolution is different between (Table 2) are low (1.5 × 10−4 d−1 ) and seem similar
normal and dry periods. Table 1 summarizes the from one recession to another, thus translating the
results. The changes in flow regimes noted for a difficulty for the system to drain stored water and
low discharge are explained by the passage from a consequently its facility to create reserves. Dynamic
flood flow (segment 2, Fig. 4) to a recession flow volumes are large and can be estimated at nearly
(segment 1, Fig. 4). The flow regime related to the first a billion cube metres. It is significant to note the
segment corresponds to a Maillet distribution [10]. variation of the dynamic reserve between 1977 and
The steeper slope of segment 1 of period 75/77 shows 1996. Indeed, in 19 years, the system lost 179 Mm3 ,
that in a normal period, the recession discharge is less that is an average annual volume of 9.4 Mm3 . This
dispersed. Segments 2, 3, 4 (period 1975–1977) and 2, loss is to be connected directly with the drought
3 (period 1995–1997) represent the flood discharge. that prevails in the area since 1981. The recession
The larger number of segments in the normal period is generally spread out: 1 month to 45 days for the
indicates that the origin of the water contributions is studied spring. This expresses a delayed response to
more diversified than in a dry period. the input signal represented by the rain.
1102 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109

The Bittit spring is the outlet of a wide and complex in the input-output relations of the system, expressing
hydrogeologic system. The water volume sustaining the non-steadiness of the system. Thus, the drought
recession is very significant (907 Mm3 in 1996), which phenomenon became a feature that should at present
implies a huge reservoir and places this aquifer among be taken into account for a sustainable management of
weakly karstified domains [11]. Sandy dolomites are water resources, particularly in semi-arid zones with
porous systems with very high regulation capacity. irregular climates.
Finally, the flood volume (VI = 0.5 Mm3 ) is much less
important compared to the dynamic volume Vr .
Spectral and correlation analysis is carried out 1. Introduction
with a daily time step and a 120-day truncation.
Fig. 6 shows the discharge simple correlograms during L’objet de cet article est, d’un point de vue général,
the two periods. One notes a rather slow decrease, de caractériser le comportement hydrodynamique des
indicating a significant memory effect. The system grands systèmes karstiques en fonction de leur état
modulates the input signal, it has a memory effect in de remplissage. Le site d’étude est représenté par
relation with the regulating role of the reserves. During la source de Bittit, localisée dans le Moyen Atlas
the period 1975–1977, the decrease is much slower marocain. Ce choix est lié à la grande taille de ce
(rk = 0.2 for k = 117 days) than during the period système et à sa complexité, à l’existence d’une longue
1995–1998 (rk = 0.2 for k = 37 days). This highlights série de pluies déficitaires et à la disponibilité de
the effect of the groundwater reserve reduction on the données sur plus de 70 ans.
damping of the input signal. This result also highlights La plaine du Sais est l’une des grandes régions agri-
the vertical heterogeneity of the karstic system and coles du Maroc. Ses ressources en eau sont fournies en
shows the significant influence of its storage level on grande partie par l’aquifère karstique du Lias du pla-
input–output relations. The spectral density function teau carbonaté du Moyen Atlas (Fig. 1). Ce causse se
(Fig. 7) shows narrow cut-off frequencies (0.11 d−1 ), vidange par un groupe de sources dont la plus impor-
which confirms that the system is very inertial. The tante est celle de Bittit (Qmoy annuel = 1,3 m3 s−1 ) [1–
calculated regulation times are of 53 days over the 3,5]. Cette dernière est actuellement exploitée pour ali-
period 1975–1977 and 37 days over the period 1995– menter en partie la ville de Meknès en eau potable.
1998. This result shows a significant reduction in the Le Moyen Atlas marocain connaît, depuis une ving-
inertia of the system during the drought period due to taine d’années, un déficit pluviométrique chronique
the reserve decrease. (Fig. 2A) dont les conséquences sont sensibles sur les
ressources en eau (Fig. 2B). La région a été ainsi sou-
4. Conclusion mise, au cours de ces dernières décades, à un régime
pluviométrique irrégulier.
The various hydrodynamic methods used in this Mieux comprendre l’évolution de ces grands sys-
study show that the Bittit hydrogeologic system is very tèmes karstiques, compte tenu de leur état de remplis-
capacitive: significant storage, sustained discharge sage, revêt une importance particulière au moment où
all year round, relatively slow flow because of the l’on parle de plus en plus de changements climatiques,
presence of a sandy dolomite base. However, in spite allant, dans la région étudiée, vers l’aridification.
of the great storage capacity of this karstic aquifer,
the spring discharge shows a regular and significant
lowering tendency due to a long drought period 2. Région d’étude
that considerably reduced the aquifer recharge. The
analyses of two discharge series allowed comparing Le Causse moyen-atlasique est constitué essentiel-
the spring differential behaviour between a normal lement de formations calcaires d’âge Liasique, qui re-
period and a dry period. This shows that the water posent sur le Trias argileux et basaltique. Une dizaine
deficit induced a reduction of the recharge area and de sources sont localisées le long du contact entre
of the reserve (estimated at 9.4 Million m3 yr−1 on le plateau carbonaté et la plaine du Sais. Le causse
average), during the two last decades, and a change moyen-atlasique constitue le secteur de recharge du
 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109 1103

Fig. 1. Carte hydrogéologique de la bordure septentrionale du causse du Moyen Atlas (extrait de la carte géologique au 1/100 000 d’El Hajeb).
1, Formations mio-plio-quaternaires ; 2, affleurements carbonatés liasiques ; 3, dépressions karstiques à remplissage argileux ; 4, source
karstique ; 5, gouffres ou avens ; 6, faille ou flexure ; 7, cours d’eau ; 8, ville.
Fig. 1. Hydrogeological map of the northern border of the Middle Atlas Plateau (from the geological map of El Hajeb at 1/100 000).
1, Mio-Plio-Quaternary formations; 2, Liasic formations; 3, karstic depressions with argillaceous deposits; 4, karstic spring; 5, sink-holes;
6, fault; 7, River; 8, Town.

système karstique, qui est vidangé par les sources. Le faille qui aurait guidé son installation. Les forages
régime des précipitations est du type méditerranéen. de reconnaissance creusés au voisinage de la source
Le mois de décembre est en général le plus arrosé montrent que le Lias est peu épais à cet endroit et que
(83 mm en moyenne à El Hajeb), alors que les mois le Trias marque une remontée, jouant ainsi un rôle de
de juillet et d’août sont les plus secs (4 à 5 mm à barrière, en empêchant l’écoulement de l’eau vers la
El Hajeb). Des averses orageuses concentrées et bru- plaine.
tales peuvent déclencher l’écoulement dans les vallées Le système hydrogéologique de Bittit jouit d’une
sèches. capacité de stockage très importante, permettant d’em-
La Fig. 2A représente les écarts pluviométriques magasiner une recharge importante provenant des
par rapport à la pluie moyenne annuelle (612 mm) pluies et des fontes de neige, sans provoquer de crues
entre 1934 et 2000. Elle confirme l’irrégularité des brutales. Le débit de la source est, en effet, bien sou-
pluies et met en évidence une diminution importante tenu durant toute l’année. La présence de dolomies sa-
des apports depuis 1980. Le déficit cumulé des pluies bleuses à la base de l’aquifère liasique plaide en faveur
depuis 1981 à El Hajeb équivaut à trois années de l’existence d’un aquifère à porosité d’interstices et
moyennes. à écoulement relativement lent. On aurait donc, dans
L’émergence Bittit jaillit au pied de la flexure l’ensemble, un système ayant un pouvoir de régula-
de Sidi Lmir, au sein d’une surface d’aplanissement tion qui permet d’amortir les pics de recharge et de
recouverte par les basaltes quaternaires qui jouent le restituer l’eau avec modération.
rôle d’un toit relativement imperméable. La source La source Bittit est caractérisée par un débit moyen
sort sous les basaltes en profitant probablement d’une élevé (1510 l s−1 calculé sur la période 1975–2001) et
1104 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109

Fig. 3. Évolution du débit annuel de la source Bittit (A) et pluie

annuelle à El Hajeb (B) (1976–2000).
Fig. 3. Annual discharge evolution of the Bittit spring (A) and
annual rainfall at El Hajeb station (B) (1976/2000).

3. Méthodologie

3.1. Étude des débits classés

Fig. 2. Déficit pluviométrique et conséquences sur les ressources

en eau. (A) Écarts pluviométriques relatifs ( %) par rapport à la La méthode des débits classés [11] consiste à
moyenne à El Hajeb (1934–2001). (B) Suivi piézométrique sur le étudier la distribution des fréquences cumulées des
forage 290/15 dans la plaine du Sais, montrant une baisse régulière débits observés à la source. Elle permet de caractériser
du niveau de la nappe consécutive à la sécheresse. la structure de l’hydrogramme unitaire. Cette méthode
Fig. 2. Rainfall deficit and its effects on water resources. (A) se révèle être un outil intéressant pour vérifier ou
Relative deviations (%) from average rainfall at El Hajeb station
non l’unité des systèmes karstiques. Ainsi, on peut
(1934–2001). (B) Groundwater monitoring at borehole 290/15 in
the plain of Sais showing a regular decrease of the water level due mettre en évidence une fuite ou un apport d’eau vers le
to the drought. système karstique étudié. Ceci se traduit sur la courbe
par des discontinuités.
La loi de probabilité ajustée est la loi normale.
En pratique, on reporte les fréquences cumulées de
une dispersion réduite (écart type = 186 l s−1 ; CV = débits sur un diagramme fonctionnel. L’observation
12,3 %). L’exutoire présente donc un écoulement re- sur la courbe de plusieurs segments de droites tra-
lativement stable, témoignant d’une bonne régulation duit l’existence de plusieurs lois d’écoulement. Ces
interannuelle par le système aquifère karstique. La changements de loi d’écoulement peuvent avoir des
Fig. 3 illustre l’évolution du débit de la source, qui origines diverses, liées soit à des modifications hy-
marque une tendance baissière régulière depuis le dé- drodynamiques internes au système karstique, comme
but de la sécheresse. Les séries chronologiques des dé- le stockage–déstockage des réserves ou l’entrée en
bits de la source Bittit, utilisées dans ce travail, recou- fonctionnement de trop-pleins, soit à des interrelations
vrent respectivement les périodes 1975–1977 et 1995– avec un système voisin, qui se traduisent par des ap-
1998. ports ou des fuites.
 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109 1105

3.2. Étude des courbes de récession est inertiel. L’inertie du système karstique est liée à
son organisation de stockage (mise en réserve) et de
L’étude des courbes de récession [11] permet de drainage.
caractériser le comportement des deux sous-systèmes
qui constituent le système karstique : la zone non sa-
turée, ou sous-système « infiltration », et la zone sa- 4. Résultats et discussion
turée ou sous-système « karst noyé ». L’écoulement
de décrue est représentatif de la manière dont s’ef- 4.1. Analyse des débits classés
fectue la vidange de la zone d’infiltration, alors que
l’écoulement de tarissement caractérise la vidange du Deux courbes des débits classés de la source Bit-
karst noyé. L’intérêt de cette méthode est de définir tit ont été construites. La première correspond aux
les caractéristiques de chaque sous-système à l’aide de deux cycles hydrologiques à pluviométrie normale
quelques paramètres simples. 09/1975–08/1977 et la seconde aux deux cycles hy-
drologiques situés en période de sécheresse 09/1995–
3.3. Analyses corrélatoires et spectrales 08/1997.
Les graphiques des débits classés de l’émergence
Les analyses corrélatoires et spectrales s’inspirent Bittit (Fig. 4) sont constitués de trois (période 1975–
directement des méthodes de traitement du signal. 1977) à quatre (période 1995–1997) segments de
Leur application à l’étude des systèmes karstiques droite. Cette organisation révèle que l’évolution des
s’est avérée tout à fait efficace pour identifier l’orga- débits est différente entre période normale et période
nisation et la dynamique des systèmes karstiques [6,8, sèche. Il est important de tenter de comprendre les réa-
9,12,13]. lités hydrodynamiques sous-jacentes à ces différentes
Le corrélogramme met en évidence la dépendance, lois. L’interprétation ci-dessous distingue le compor-
les uns avec les autres, des évènements successifs pour tement de l’émergence par rapport à son débit moyen
des intervalles de temps croissants. Analytiquement, pour les deux périodes étudiées. Le Tableau 1 résume
il montre l’évolution du coefficient rk d’autocorréla- les résultats obtenus.
tion. Si le corrélogramme décroît très vite, il indique Les changements des lois d’écoulement relevés au
des événements indépendants les uns des autres. En niveau des faibles débits s’expliquent par le passage
revanche, une décroissance lente caractérise des évé- d’un écoulement de décrue (segment 2, Fig. 4) à un
nements dépendants les uns des autres. La valeur rk = écoulement de tarissement (segment 1, Fig. 4). La
0,2 est une valeur remarquable, car des événements loi d’écoulement relative à ce premier segment de
peuvent être considérés indépendants quand rk < 0,2. droite correspond à une loi de Maillet [10]. La pente
L’analyse spectrale correspond à un changement plus importante du segment 1 de la période 1975–
d’un domaine de temps à un domaine de fréquence 1977 montre qu’en période normale, les débits de
par une transformée de Fourrier du corrélogramme. tarissement sont moins dispersés.
La fonction spectrale S(f ) de densité représente la Les segments 2–4 (période 1975–1977) et 2, 3 (pé-
décomposition de la variance de la série chronologique riode 1995–1997) représentent les débits de décrue. Le
dans le domaine des fréquences. Les phénomènes nombre de segments, plus grand en période normale,
périodiques apparaissent sous forme de pics dans le est une indication d’origines des apports, plus diversi-
graphe de S(f ). Quand ces pics caractérisent le signal fiées qu’en période sèche.
de sortie, ceci montre comment le signal d’entrée Dans le domaine des débits élevés, une disconti-
a été modifié par le système. Plus la fréquence de nuité avec augmentation de pente (segments 2 et 3,
coupure (la fréquence telle que S(f ) est presque période 1975–1977), signifie une fuite d’eau vers l’ex-
nulle) est faible, plus le système est inertiel. La térieur du système et peut traduire l’existence probable
fonction S(f ) permet également le calcul du temps de « trop-pleins » (par exemple un drain au-dessus de
de régulation du système, qui représente la durée la source principale). En période sèche, on observe
d’influence du signal d’entrée [Treg = S(f = 0)/2]. une discontinuité avec réduction de pente entre les
Plus le temps de régulation est long, plus le système segments 2 et 3. Cette particularité signifie un apport
1106 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109

Fig. 4. Distribution des débits en période normale (A : 1975–1977) et en période sèche (B : 1995–1997).
Fig. 4. Discharges distribution during a normal period (A: 1975–1977) and a drought period (B: 95/97).

Tableau 1
Analyse de la distribution des débits (l s−1 )
Table 1
Analysis of the discharges (l s−1 ) distribution
Débit moyen (l s−1 ) / 1er intervalle / 2e intervalle / 3e intervalle / 4e intervalle /
Average discharge (l s−1 ) 1st interval 2nd interval 3rd interval 4th interval
Chronique / Series 1975–1977 1737 < 1640 1640–1840 1840–1920 > 1920
Chronique / Series 1995–1997 1400 < 1400 1400–1515 > 1515 –

d’eau dans le système et peut traduire la mise en ré- faibles des limites des intervalles de débits. Les mo-
serve temporaire d’un certain volume d’eau dans des difications d’impluvium envisagées ci-dessus sont des
structures annexes lors des crues et sa libération dès le hypothèses, qui n’ont pu être vérifiées faute d’essais
début de la décrue. de traçage dans le système [3].
La discontinuité entre les segments 3 et 4 de
la période 1975–1977 traduit l’apport d’un volume 4.2. Analyse des courbes de récession
d’eau extérieur, uniquement lors des périodes de crue
(modification des limites du bassin, pluie, fonte de Le choix des différentes crues qui ont servi à l’ana-
neige). lyse des courbes de récession a été fait de telle fa-
Durant la période 1975–1977, les débits de l’émer- çon qu’elles représentent les amplitudes maximales
gence étaient plus importants qu’actuellement et ils de chaque cycle hydrologique et précèdent de longues
s’organisaient en quatre segments de droite. L’alimen- périodes d’étiage (Fig. 5). Les coefficients de tarisse-
tation de la source était plus soutenue et probablement ment (Tableau 2) sont faibles (1,5 × 10−4 j−1 ) et sem-
plus diversifiée avant les années 1980 que durant la blent reconductibles d’un tarissement à l’autre, tradui-
période de sécheresse actuelle. En effet, si le niveau sant ainsi la vidange lente des réserves et, par consé-
de l’aquifère baisse à cause d’une réduction des ap- quent, la tendance au stockage. Les volumes dyna-
ports, des drains secondaires se retrouvent dénoyés et miques sont considérables et peuvent être estimés à
l’aire d’alimentation s’en retrouve réduite. Ceci se tra- près d’un milliard de mètres cube. Il est important
duit sur le graphique des débits classés par la réduc- de noter la variation de la réserve dynamique de la
tion des lois d’écoulement (trois segments au lieu de source entre 1977 et 1996. En effet, en l’espace de
quatre) et par un changement vers des valeurs plus 19 ans, le système a perdu 179 millions de mètres
 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109 1107

Fig. 5. Courbes de récession de la source Bittit en période normale (A) et en période sèche (B).
Fig. 5. Recession curves of the Bittit spring in a normal period (A) and in a drought period (B). ‘décrue’: flood recession; ‘tarissement’: baseflow.

Tableau 2 cube, soit un volume moyen annuel de 9,4 Mm3 . Cette

Analyse des courbes de récession perte est à mettre en relation directe avec la séche-
Table 2 resse qui sévit dans la région depuis 1981. L’absence
Analysis of recession curves de forage sur le causse ne permet pas de visualiser
Date début de décrue / 22/05/77 19/03/96 la réaction de la nappe à cette sécheresse. L’aquifère
Recession starting time du Lias se prolonge sous la plaine du Sais, vers le
Durée de la décrue ti (jour) / 49 34 nord (Fig. 1), où il devient captif sous les formations
Flood recession (d)
mio-plio-quaternaires. Le suivi piézométrique dans le
Durée du tarissement (jour) / > 53 147
Baseflow recession (d) forage 290/22 (X = 497,5 ; Y = 353,8) depuis fé-
Débit de début Q0 de la décrue (l s−1 ) / 2072 1630 vrier 1968 (Fig. 2B) montre une baisse régulière de
Flood starting discharge Q0 (l s−1 ) la charge de la nappe depuis 1980 (−2,7 m an−1 en
Débit de base QRo du tarissement (l s−1 ) / 1898 1477 moyenne).
Baseflow discharge QRo (l s−1 ) La décrue est généralement étalée : 30 à 45 j pour
q0 (l s−1 ) = Q0 − QRo 174 153
Débit QRo à la fin de la décrue (l s−1 ) / 1886 1470
la source étudiée. Ceci traduit une réponse retardée
Discharge QRo at the end of the flood (l s−1 ) au signal d’entrée représenté par la pluie. Les durées
Coefficient de tarissement α(10−3 ) / 0,15 0,14 des décrues peuvent varier d’une année hydrologique à
Recession coefficient α(10−3 ) une autre, en fonction de l’importance de l’impluvium
Volume dynamique VR (Mm3 ) / 1086,3 907,2 et de sa répartition temporelle.
Dynamic volume VR (Mm3 ) La source Bittit offre un système hydrogéologique
Volume de décrue Vi (Mm3 ) / 0,7 0,5
étendu et complexe. Le volume d’eau soutenant le
Flood volume Vi (Mm3 )
tarissement est très important (907 Mm3 en 1996),
1108 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109

Fig. 6. Corrélogrammes simples des débits de la source Bittit Fig. 7. Spectres de densité de variance des débits de la source
(m = 120) (septembre 1975–août 1977 : n = 731 et décembre Bittit (m = 120) (septembre 1975–août 1977 : n = 731 et décembre
1995–février 1998 : n = 770) ; l’abscisse est graduée en jours. 1995–février 1998 : n = 770) ; l’abscisse est graduée en jours−1 .
Fig. 6. Simple correlograms of the Bittit spring discharges Fig. 7. Spectral density functions of the Bittit spring discharges
(m = 120 d) (September 1975–August 1977: n = 731 and Decem- (m = 120 d) (September 1975–August 1977: n = 731 and Decem-
ber 1995–February 1998: n = 770); the horizontal axis is in days. ber 1995–February 1998: n = 770); the horizontal axis is in days−1 .

ce qui implique un immense réservoir et place cet sensible de l’inertie du système durant la période de
aquifère parmi les domaines peu karstifiés [12]. Les sécheresse, due à la baisse des réserves.
dolomies sableuses constitueraient un système poreux Les paramètres obtenus en analysant les débits
à pouvoir régulateur très élevé. Enfin, le volume de de la source montrent un effet mémoire grand à
décrue (Vi = 0,5 Mm3 ) est insignifiant par rapport au très grand (37 à 117 j), une bande spectrale étroite
volume dynamique Vr . (0,11 j−1 ), et un temps de régulation de 37 à 53 j. Ces
données peuvent être comparées à celles des aquifères
4.3. Analyse corrélatoire et spectrale carbonatés étudiés par ailleurs [6,11]. La source Bittit
peut être classée parmi les aquifères à effet mémoire
L’analyse est réalisée à court terme, avec un pas important de type Fontestorbes.
de temps journalier et une troncature de 120 j. La
Fig. 6 montre les corrélogrammes simples des débits
de la source durant les deux périodes d’observation. 5. Conclusion
On note une décroissance assez lente, témoignant d’un
effet mémoire important. Le système module le si- Les différentes méthodes hydrodynamiques utili-
gnal d’entrée ; il a un effet mémoire qui est en rela- sées dans cette étude montrent que le système hydro-
tion avec le rôle régulateur des réserves. Durant la pé- géologique de Bittit est très capacitif : stockage impor-
riode 1975–1977, la décroissance est beaucoup plus tant, débit régulé et soutenu durant toute l’année, écou-
lente (rk = 0,2 pour k = 117 j) qu’au cours de la pé- lement relativement lent à cause de la présence à la
riode 1995–1998 (rk = 0,2 pour k = 37 j). Ceci traduit base d’une assise de dolomies sableuses, karstification
l’effet de la diminution de la réserve de la nappe sur peu développée. Cependant, malgré la grande capacité
l’amortissement du signal d’entrée. Ce résultat sou- de réserve de cet aquifère karstique, le débit à l’exu-
ligne l’hétérogénéité verticale des systèmes karstiques toire marque une tendance régulière et significative à
et montre l’influence significative de leur état de rem- la baisse, due à une longue période de sécheresse, qui
plissage sur les relations entrées–sorties. réduit considérablement l’alimentation de l’aquifère.
Les spectres de densité de variance de la source L’analyse de deux chroniques de débits de l’émer-
(Fig. 7) montrent des fréquences de coupure ou bandes gence a permis de comparer le comportement diffé-
spectrales étroites (0,11 j−1 ), et confirment qu’il s’agit rentiel de la source entre une période normale et une
d’un système inertiel. Les temps de régulation calculés période sèche. Il en ressort que le déficit hydrique a
sont de 53 j pour la période 1975/77 et 37 j pour la produit une réduction de l’aire d’alimentation, une di-
période 1995–1998. Ce résultat montre une réduction minution des réserves estimée à 9,4 millions de mètres
 F. Amraoui et al. / C. R. Geoscience 336 (2004) 1099–1109 1109

cube par an en moyenne durant les deux dernières dé- [2] F. Amraoui, L’aquifère karstique du Moyen Atlas tabulaire :
cennies et un changement dans les relations entrée– problème de la turbidité affectant les sources du complexe
sortie du système, dû à l’influence significative de son Ribaa-Bittit (Maroc), Les Cahiers de la Recherche, Éditions
du Rectorat de Casablanca (Maroc), Série A, Sciences et
état de remplissage, mettant en évidence son caractère Techniques 5 (2003) 81–92.
non stationnaire. [3] F. Amraoui, M. Razack, L. Bouchaou, Turbidity dynamics in
Les méthodes mises en œuvre dans ce travail (étude karstic systems. Example of Ribaa and Bittit springs in the
des débits classés et des courbes de récession, analyses Middle Atlas (Morocco), Hydrol. Sci. J. 48 (6) (2003) 971–
corrélatoires et spectrales) ont permis de mieux com- 984.
[4] C. Batiot, C. Emblanch, B. Blavoux, Carbone organique to-
prendre l’évolution d’un grand système karstique tel tal (COT) et magnésium (Mg2+ ) : deux traceurs complémen-
que celui de Bittit. D’autres méthodes, notamment les taires du temps de séjour dans l’aquifère karstique, C. R. Geos-
traceurs naturels [4,7], peuvent contribuer à cet objec- cience 335 (2) (2003) 205–210.
tif. Leur application fera l’objet d’un prochain travail. [5] A. Benjbara, Observations sur les sources de la plaine de
Ainsi, le phénomène de sécheresse est devenu un Fès–Meknès, Rev. Eau et Développement, in : Direction de
l’Hydraulique (Maroc), 1987, pp. 31–37.
paramètre dont il faut tenir compte pour une gestion [6] L. Bouchaou, P. Chauve, J. Mudry, J. Mania, Y. Hsissou,
durable de la ressource, notamment en zone semi-aride Structure et fonctionnement d’un hydrosystème karstique de
à climat irrégulier, où les changements climatiques montagne sous climat semi-aride : cas de l’Atlas de Beni-
tendent vers une aridification. Mellal (Maroc), J. Afr. Earth Sci. 25 (2) (1997) 225–236.
[7] A.M. Delaunay, M. Calleja, J.-P. Faillat, S. Pistre, P. Cour,
Traçage pollinique naturel dans un hydrogéosystème karstique
thermo-minéral, C. R. Acad. Sci., Ser. IIa 331 (7) (2000) 467–
Remerciements 472.
[8] D. Labat, R. Ababou, A. Mangin, Analyse multirésolution
Ce travail a bénéficié du soutien financier de croisée de pluies et débits de sources karstiques, C. R.
l’Action intégrée France–Maroc AI 221/STU/00. Les Geoscience 334 (8) (2002) 551–556.
[9] M. Larocque, A. Mangin, M. Razack, O. Banton, Contribution
auteurs tiennent à remercier chaleureusement le Pr. of correlation and spectral analyses to the regional study of
G. de Marsily et le second réviseur anonyme, pour a large karst aquifer (Charente, France), J. Hydrol. 205 (3–4)
l’ensemble de leurs remarques et suggestions, qui ont (1998) 217–231.
permis d’améliorer cet article. [10] E. Maillet, Essais d’hydraulique et fluviale, Hermann, Paris,
1905, 218 p.
[11] A. Mangin, Étude des débits classés d’exutoires karstiques
portant sur un cycle hydrologique, Ann. Spéléol. 26 (2) (1971)
Références 283–329.
[12] A. Mangin, Pour une meilleure connaissance des systèmes
[1] F. Amraoui, L’aquifère karstique du Moyen Atlas tabulaire : hydrologiques à partir des analyses corrélatoires et spectrales,
problème de la turbidité affectant les sources du complexe J. Hydrol. 67 (1984) 25–43.
Ribaa-Bittit (Maroc), Proc. 7th Conf. On Limestones Hyd. and [13] A. Padilla, A. Pulido-Bosch, Study of hydrographs of karstic
Fissured Media, Besançon, France, Sci. Tech. Environ. 113 aquifers by means of correlation and cross-spectral analysis,
(2001) 45–48 (Mém. H.S.). J. Hydrol. 168 (1995) 73–89.