GEOLOGIJA 25/2, 309—325 (1982), Ljubljana
UDK 551.44:551.466(497.12)=z40
Transmission des ondes de marée à travers l'aquifère côtier
de Kras
Razširjanje valov plimovanja skozi obalni vodonosnik Krasa
Primož Krivic
Geološki zavod, 61000 Ljubljana, Parmova 33
Résumé
L'étude de la propagation d'ondes de marée marine dans l'aquifère
côtier de Kras a permis de caractériser le comportement hydrodynamique
général de l'aquifère et d'évaluer certains paramètres caractéristiques
du milieu, tels que la diffusivité. En étudiant la propagation de ces ondes,
nous avons examiné plusieurs modèles hydrogéologiques. Le modèle qui
correspond le mieux aux conditions géologiques connues est celui d'une
nappe captive avec un effet d'étanchéité partielle du toit de la zone
drainante. Ce modèle, ainsi que celui d'une nappe libre avec un effet
d'égoixttement, permettent la meilleure reconstitution des marées ob-
servées.
Kratka vsebina
S študijem razširjanja valov plime in oseke v obalnem vodonosniku
Krasa smo lahko ugotovili splošne hidrodinamične lastnosti vodonosnika
in določili nekatere značilne parametre, kot na primer difuzivnost. Pri
analizi periodičnih nihanj gladine podtalnice smo uporabili več hidro-
geoloških modelov. Znanim geološkim razmeram najbolje ustreza model
vodonosnika pod pritiskom z delno vodotesnostjo krovnine drenažne cone.
Ta model omogoča tudi najboljšo rekonstrukcijo v vodonosniku opazo-
vanih nihanj, enako kot tudi model vodonosnika s prosto gladino in vpli-
vom izcejanja.
Avertissement. Cet article constitue la deuxième partie d'un mémoire sous titre:
«Étude hydrodynamique d'un aquifère karstique» présenté par l'auteur en vue d'une
soutenance de thèse de Docteur-ingénieur à l'Université Montpellier II, France, année
1981. La suite paraîtra dans le fascicule suivant de la GEOLOGIJA.
Pojasnilo. Članek predstavlja drugi del doktorske disertacije z naslovom: »Studija
hidrodinamike kraškega vodonosnika«, ki jo je avtor zagovarjal leta 1981 na Univerzi
Montpellier II v Franciji. Nadaljevanje bo izšlo v naslednjem zvezku revije GEOLO-
GIJA.
310	Primož Krivic
1. Objectif de l'étude
Lorsque, sur un littoral les nappes aquifères se déversent en mer, elles sont
ainsi en contact avec les eaux libres océaniques dont les fluctuations de
niveaux, conséquences des marées, se répercutent sous forme d'ondes de dé-
formation, sur les surfaces piézométriques.
L'étude de la propagation de ces ondes peut permettre, dans des conditions
données, de caractériser le comportement hydrodynamique de l'aquifère et
d'évaluer certains paramètres caractéristiques du milieu, tel que la diffusività.
Ce problème sera abordé dans ce chapitre, à propos d'observations effec-
tuées sur des forages des karsts côtiers. En effet, l'aquifère karstique est in-
fluencé, tout au moins en partie, comme nous allons le voir, par les marées
du Golfe de Trst (Trieste), malgré la faible amplitude de celles-ci.
2. Relations hydrodynamiques entre nappe aquifère et plan d'eau libre
Rappels théoriques
Considérons le système formé par un aquifère semi-infini, d'épaisseur cons-
tante en communication avec la mer sur toute sa hauteur.
La fonction h = h (x, t) qui définit la surface piézométrique est solution de
l'équation générale de diffusività hydraulique (J. B-oussinesq, 1877; J.
Ferris, 1951):
S est le coefficient d'emmagasinement
T la transmissivité
h la charge hydraulique
La résolution de l'équation fondamentale est obtenue par deux approches
différentes, selon la forme de l'évolution dans le temps du plan d'eau libre,
—	soit que cette évolution est de forme quelconque,
—	soit qu'elle puisse être assimilée à une sinusoïde.
2.1. Evolution de forme quelconque
La résolution de l'équation différentielle ci-dessus dépend de la complexité
des conditions aux limites.
Le cas d'un modèle simple est défini par les hypothèses suivantes:
—	couche aquifère homogène, isotrope et de même épaisseur dans toute la
zone influencée par la marée,
—	limite rectiligne entre plan d'eau et aquifère,
—	aquifère captif semi-infini,
—	libéralisation de l'eau instantanée,
—	pas d'autre influence que celle de la marée,
—	écoulement laminaire.
Transmission des ondes de marée à travers l'aquifère côtier de Kras	311
Fig. 1. Propagation d'ondes piézométriques dans l'aquifère côtier
SI. 1. Razširjanje piezometričnih valov v obalnem vodonosniku
Fig. 2. Propagation d'ondes piézométriques dans l'aquifère côtier
SI. 2. Razširjanje piezometričnih valov v obalnem vodonosniku
312
Primož Krivic
Fig. 3. Vue générale de la vallée de Brestovica. A l'avant
plan, au centre, le piézomètre B-4, permettant l'enregistre-
ment des marées dans l'aquifère, est situé dans une doline.
Au delà, la mer Adriatique et la bouche de la rivière Soča
à l'arrière plan
SI. 3. Dolina Brestovice. V ospredju, v srednjem delu slike
je v manjši vrtači piezometer B-4, ki omogoča meritve plimo-
vanja v vodonosniku, v ozadju pa je Jadransko morje in ustje
reke Soče
On considère la fonction h (x, t) satisfaisant aux conditions aux limites
suivantes:
h (x, t) = 0 pour t < 0
h (0, t) = ho pour t> 0
Après une remontée brusque ho sur une nappe initialement horizontale, ou
après un abaissement brusque ho de la cote du plan d'eau libre, on obtient,
pour les charges h à la distance x de cette limite, l'expression suivante:
qui représente la résolution de l'équation fondamentale obtenue par double
transformation de CARSON-LAPLACE (M. Bonnet etJ. Schneider.
1969).
Transmission des ondes de marée à travers l'aquifère côtier de Kras
313
Fig. 4. Piézomètre naturel Drča jama. Exploré par les
plongeurs jusqu' à la profondeur de 15 m sous le niveau de la
mer. Les fluctuations de niveau piézométrique sont enregi-
strés dans le forage B-5 qui recoupe le chenal immergé de
l'aven Drča jama à la profondeur de 66,7 m à 72,5 m
SI. 4. Naravni piezometer Drča jama. Potapljači so raziskali
podvodni rov do globine 15 m pod gladino morja. Nihanja
gladine podtalnice beleži limnigraf na vrtini B-5, ki je pre-
vrtala podvodni kanal Drče jame v globini od 66,7 m do
72,5 m
La fonction erfc (^^fonction erreur complémentaire>>) est définie par l'inté-
grale suivante:
En appliquant le théorème de superposition des états transitoires, la solution
élémentaire précédente peut être étendue aux cas plus complexes.
La solution générale h (x, t) devient alors:
314	Primož Krivic
La méthode de transformation de CARSON-LAPLACE a été appliquée par
différents auteurs pour le calcul de diffusi vité, notamment en milieu poreux:
M. Bonnet etJ. Schneider (1969), R. D eg allier (1969), P. Peau-
decerf (1973), mais aussi en milieu fissuré: J. C. G arni er et C. L.
Clarion (1967), M. Razack, C. Drogue, C. Romariz et C. Al-
meida (1980).
2.2. Evolution de forme sinusoïdale
2.2.1. Nappe captive
La propagation d'ondes piézométriques dans l'aquifère côtier peut être
représentée par la résolution de l'équation différentielle de base de l'hydro-
dynamique souterraine, avec les conditions aux limites suivantes (J. B o u s -
s i nés q, 1877; J. Ferris, 1951):
ho étant la demi amplitude et to la période de fluctuation du plan d'eau libre.
En faisant les mêmes hypothèses que dans le chapitre précédent, on obtient
la solution suivante pour le mouvement de la surface piézométrique dans une
nappe captive (modèle BOUSSINESQ):
d'où la demi amplitude hx des fluctuations de surface piézométrique à la
distance x du rivage:
et la vitesse V de transmission de l'onde:
La diffusivi té D == T/S peut être aisément calculée,
— soit à partir de l'amortissement de l'amplitude:
— soit à partir du déphasage:
Cette démarche, très classique, à partir de l'équation de BOUSSINESQ,
a été reprise et développée par J. Ferris (1951), puis largement utilisée
pour l'interprétation des fluctuations périodiques dans les nappes alluviales
(G. Trup in -1969-, R. D eg allier -1969-, E. De Cazenove -1971-,
Transmission des ondes de marée à travers l'aquifère côtier de Kras	315
R. Degallier -1978-) ainsi que dans les aquifères karstiques (L. M o u 1 -
lard, B. Mijatović, R. Kareh etB. Masaad -1965-, J. C. Gar-
nier et C. L. Clarion-1967-, M. Razack, C. Drogue, C. Roma-
riz et C. Almeida -1980-).
2.2.2. Nappe semi-captive
Supposons maintenant que l'aquifère étudié (perméabilité horizontale K,
épaisseur M) soit en communication à travers son éponte semi-perméable (per-
méabilité verticale K', épaisseur M') avec une nappe à niveau d'eau constant.
C'est le schéma classique de la drainance.
D'après E. De Cazenove (1971), la grandeur caractéristique B (distance
de drainance):
est liée à un paramètre ô différent de zéro, caractérisant l'intensité de la
drainance, par:
L'équation de diffusivité et sa solution s'écrivent:
Les valeurs de l'amortissement et du déphasage permettent de déduire
xo/(5 d'une part, et xo ô d'autre part, d'où xq, <5 et B. A partir de ces valeurs
la diffusivité D peut être calculée aisément.
2.2.3. Nappe libre
Dans la mesure où les variations périodiques du niveau h sont faibles par
rapport à l'épaisseur de la zone saturée de l'aquifère, il est admis d'appliquer
les équations valables pour les nappes captives (E. De Cazenove, 1971).
2.2.4. Nappe libre avec un effet d'égouttement
Dans le cas d'un aquifère libre, avec un effet d'égouttement non négligeable,
en communication directe avec la mer (la méthode d'interprétation présentée
ci-après suppose maintenues toutes les autres hypothèses du § 2.1.), l'équation
de la diffusivité s'écrit (par analogie avec le phénomène de drainance):
316
Primož Krivic
Fig. 5. Aven Drča jama. On aperçoit, derrière l'échelle limni-
métrique, le plan d'eau rabattu par le pompage d'essai dans le
forage B-5
SI. 5. Drča jama. Za vodomerno letvijo lahko opazimo vodno
gladino sifonskega jezera, ki se je znižala zaradi črpalnega
poizkusa na vrtini B-5
ß qui est un paramètre différent de zéro intervenant dans le facteur d'égout-
tement B':
A partir de l'amortissement et du déphasage, on obtient les valeurs de xo,
ß et B'. La diffusivité T/S peut donc être calculée aisément.
Transmission des ondes de marée à travers l'aquifère côtier de Kras	317
3. Conditions expérimentales
3.1. Données géologiques et hydrogéologiques
Rappelons que l'épaisseur de la série carbonatée aquifère est très grande,
car ni les forages de reconnaissance ni la prospection géophysique (sondages
électriques) n'ont pu découvrir le substratum imperméable (F. Drobne,
P. Krivic, D, Ravnik et U. Premru, 1977). Il n'y a pas non plus
de recouvrement sauf des dépôts d'argiles de décalcification dans les dolines,
mais qui ne font pas partie de la zone saturée de l'aquifère. Il est donc plausible
d'envisager l'hypothèse d'une nappe libre avec des variations piézométriques
négligeables par rapport à l'épaisseur de la zone saturée.
Mais il faut noter qu'il y a eu probablement, au moins sur la zone proche
du littoral, un développement karstique important sous le niveau actuel de la
mer, du fait des régressions marines du Tertiaire et du Quaternaire.
Effectivement, les plongeurs spéléologues connaissent de nombreux chenaux
karstiques à des profondeurs sous la mer de l'ordre de 30 à 50 m. De la sorte,
on peut envisager l'existence d'une zonation plus ou moins horizontale de la
conductivité hydraulique du karst, et des circulations actuelles en conduits
localement en charge.
L'écoulement n'est évidemment pas totalement en charge, car les calcaires
affleurent largement et sont parcourus de fractures ouvertes subverticales.
3.2. Marée dans le Golfe de Trst (Trieste)
L'Adriatique est une mer intérieure, à forme allongée, dirigée du NW vers
le SE. Dans la partie septentrionale, les profondeurs sont modérées, mais
l'amplitude de marées est plus grande que celle d'autres régions méditerranéen-
nes (F. Mosetti et M. Pic otti, 1960).
Cette amplitude augmente, peu à peu, vers le Nord jusqu'au Golfe de Trieste
où elle atteint un maxima de 0,6 et même de 0,9 m. Il s'agit d'une marée
semi-diurne avec une période de 12 h 25 m et (ou) d'une marée mixte avec
deux maxima, dont l'un est prédominant. L'allure très complexe de la marée
observée (fig. 1) est dûe notamment à la forme particulière de l'Adriatique et
également à l'influence des vents dominants (bora, siroco) qui agissent sur le
niveau marin en l'abaissant ou en le surélevant selon leur direction.
3.3. Observations disponibles
Les premiers renseignements sur les fluctuations de la nappe datent d'avril
1977, époque à laquelle nous avions installé des limnigraphes permettant un
enregistrement continu des niveaux piézométriques dans les forages.
Ensuite, dans le cadre des travaux de reconnaissance de l'aquifère côtier,
des enregistrements (voir fig. 2) ont été réalisés dans les piézomètres B-4 (fig. 3)
et B-5 (fig. 4, 5), distants de 4000 m de la côte, ainsi que sur le marégraphe
situé dans le Golfe de Trst (Trieste).
Durant la période d'observations du 28. 4. 1977 au 12. 12. 1978, le phénomè-
ne de marée est apparu quatorze fois avec des amplitudes sur les forages de
2 à 10 cm. La période ininterrompue des fluctuations du niveau piézométrique
318
Primož Krivic
la plus longue a duré du 25. 10. 1978 au 12. 12. 1978, soit 49 jours. Cette période
correspond à un étiage d'automne prolongé.
Pour l'application pratique des modèles d'interprétation présentés dans le
paragraphe 2, nous avons utilisé les enregistrements réalisés dans le forage
B-4 qui sont d'une très bonne qualité.
4. Propagation de la marée dans l'aquifère karstique
4.1. Hypothèse d'une nappe captive
4.1.1. Application de la transformation de CARSON-LAPLACE
Nous avons d'abord appliqué la méthode utilisant la double transformation
de CARSON-LAPLACE. Les calculs effectués sur les données du piézomètre
B-4 distant de 4000 m de la côte fournissent une diffusivité: D = 193 mVs.
A l'aide de cette valeur de T/S et des conditions aux limites (marée marine
observée), il est possible de déterminer la courbe théorique des fluctuations du
niveau piézométrique dans le B-4. Les calculs sont effectués à l'aide du
programme informatique HYGP1 (P. Peaudecerf, 1973).
Les fluctuations piézométriques observées et calcules sont reportées sur la
figure 6. La courbe théorique ne coïncide pas avec la courbe observée. C'est
donc que l'hypothèse d'une nappe captive ne subissant aucune autre influence
Fig. 6. Reconstitution des fluctuations sur le forage B-4 par application de la trans-
formation de Carson-Laplace
SI. 6. Rekonstrukcija nihanj gladine podtalnice v vrtini B-4 po metodi Carson-Laplace
Transmission des ondes de marée à travers l'aquifère côtier de Kras
319
Fig. 7. Reconstitution de la marée dans le Golfe de Trst
SI. 7. Rekonstrukcija plimovanja morja v Tržaškem zalivu
Fig. 8. Reconstitution de la marée marine
SI. 8. Rekonstrukcija plimovanja morja
320	Primož Krivic
que celle de la marée n'est pas verifiée et qu'il faut faire intervenir d'autres
facteurs; soit le phénomène de débit retardé et d'égouttement, soit l'effet
d'étanchéité partielle du toit de la zone drainante, comme nous le verrons plus
loin.
4.1.2. Nappe captive (évolution sinusoïdale)
La diffusivité est calculée à l'aide de l'amortissement de l'amplitude et du
déphasage. La marée est très irrégulière dans le Golfe de Trieste. Sa forme
précise est obtenue en analysant les enregistrements avec des intervalles de
temps de 6 à 12 heures (fig. 7). Pour chaque intervalle sont recherchées l'ampli-
tude et la période apparentes (fig. 8). Ces valeurs permettent de calculer la
grandeur xo qui contient l'inconnue T/S que l'on cherche à déterminer.
L'interprétation consiste ensuite à reproduire les fluctuations dans le forage
B-4. Les résultats obtenus sont présentés dans le tableau 1 et les fig. 9 et 10.
Il faut noter que les paramètres calculés d'après le déphasage sont systé-
matiquement plus élevés que ceux obtenus à l'aide de l'amortissement, les
moyennes étant une fois et demi plus fortes.
4.2. Hypothèse d'une nappe libre avec un effet d'égouttement
Dans ce modèle, on envisage l'hypothèse de libération progressive de l'eau
dûe au drainage de la hauter dénoyée. Ce phénomène d'égouttement se mani-
feste notamment dans les zones à minces fissures où la libération de l'eau est
retardée dans une certaine mesure.
Comme précédemment, les rapports d'amplitudes et les déphasages ont été
reproduits par le modèle. La courbe de réponse théorique a été ainsi déterminée.
Elle coïncide bien avec la courbe réellement observée (fig. 11). Les valeurs
obtenues sont présentées dans le tableau 2.
Fig. 9. Reconstitution des fluctuations au forage B-4 (nappe captive, modèle Bous-
sinesq)
SI. 9. Rekonstrukcija nihanj v vrtini B-4 po Boussinesqovem modelu za vodonosnik
pod pritiskom
Transmission des ondes de marée à travers l'aquifère côtier de Kras	321
Fig. 10. Reconstitution des fluctuations au forage B-4 (nappe captive, modèle Bous-
sinesq)
SI. 10. Rekonstrukcija nihanj v vrtini B-4 po Boussinesqovem modelu za vodonosnik
pod pritiskom
Fig. 11. Reconstitution des fluctuations au forage B-4 (nappe libre avec un effet
d'égouttement, marée sinusoïdale)
SI. 11. Rekonstrukcija nihanj v vrtini B-4 za vodonosnik s prosto gladino in izceja-
njem ter sinusoidalnim plimovanjem
Ces valeurs sont plus cohérentes que les précédentes et les écarts à leur
moyenne sont inférieurs à 10 "/o. On obtient ainsi des valeurs de B' et D moins
dispersées que dans les calculs antérieurs.
4.3. Hypothèse d'une nappe captive avec un effet d'étanchéité partielle du toit
de la zone drainante
L'aquifère a été décrit comme étant de forme d'un niveau très conducteur
— et capacitif — surmonté par une zone relativement moins karstifiée.
8 — Geologija 25/2
322
Primož Krivic
Tableau 1. Résultats d'interprétation d'après la méthode de Bous-
sinesq (nappe captive)
Tabela 1. Rezultati interpretacije po Boussinesqovi metodi za vodo-
nosnik pod pritiskom
Aussi, peut-on retenir par analogie, un modèle de nappe captive présentant
une étanchéité partielle à son toit. Avec:
X = x/xq: distance réduite;
Un facteur d'étanchéité B" élevé correspond à un écoulement plutôt en
charge. Les ondes piézométriques sont alors peu modifiées. Inversement, si B"
est faible, l'influence de l'effet d'étanchéité partielle devient considérable.
Transmission des ondes de marée à travers l'aquifère côtier de Kras
323
Tableau 2. Résultats d'interprétation d'après la méthode tenant
compte de l'effet d'égouttement (nappe libre)
Tabela 2. Rezultati interpretacije po metodi, ki upošteva izcejanje
pri vodonosnikih s prosto gladino
La dimension de B" est celle d'une longueur: L
S": coefficient d'emmagasinement retardé
1/a: indice de retard
La solution générale h(x, t) est la même que dans le cas d'une nappe libre
avec un effet d'égouttement; les calculs sont identiques.
5. Conclusion
En étudiant la propagation d'ondes piézométriques dans l'aquifère côtier,
sous l'effet des marées marines, nous avons examiné plusieurs modèles hydro-
géologiques. Notre démarche a été de chercher par tâtonnements une super-
position entre les courbes expérimentales et les figures dérivées des lois théo-
riques. Seul le cas d'écoulements parallèles a été traité, et les paramètres
hydrodynamiques ont été considérés comme constants dans l'espace et dans
le temps.
Le modèle qui correspond le mieux aux conditions géologiques connues est
celui d'une nappe captive avec un effet d'étanchéité partielle du toit de la zone
324	Primož Krivic
drainante. Ce modèle, ainsi que celui d'une nappe libre avec un effet d'égout-
tement, permettent la meilleure reconstitution des marées observées. Le choix
entre ces deux modèles reste pourtant difficile car il est bien certain que la
nappe n'est ni parfaitement captive ni tout à fait libre.
Nous donnons dans la suite quelques justifications de notre choix, tout
d'abord sur la différence entre la célérité de propagation d'ondes piézométriques
et la vitesse propre des particules d*eau. Il y a évidemment une différence
importante entre la célérité d'ondes qui se propagent dans l'aquifère étudié et
la vitesse effective de l'eau. Des célérités élevées indiqueraient un écoulement
en charge. Le traçage effectué à partir de la perte de la Notranjska Reka (F.
B i d o v e c , 1965) a donné le temps de passage de l'eau entre cette perte et
les sources de Timav de 16 jours, d'où la vitesse apparente de 90 m/h. En ce
qui concerne la célérité, nous avons montré (P. Krivic, 1981) que le temps
de passage d'ondes de crue sur le même chemin n'est que de deux jours et la
célérité de 730 m/h.
A signaler aussi que l'effet d'étanchéité partielle n'est pas très important
dans l'aquifère étudié; ß étant voisin de 1 (/? = 1Д41). L'écoulement se fait
donc surtout en charge.
Les essais par pompage dans le forage B-4 ont fourni une transmissivité
apparente de 10~^ mVs. Si l'on estime la porosité efficace d'une nappe libre
à 1 "/a, la diffusivité serait donc de l'ordre de 10 mVs. Pour une nappe captive
avec un coefficient d'emmagasinement 10~® < S < 10"^ la diffusivité D = T/S
serait de 100 mVs à 1000 m^/s.
D'après le modèle, la valeur de diffusivité de l'ordre de 200 m^s pourrait
être retenue pour la zone drainante dans laquelle se situent le forage B-4
et les sources côtières. Cette valeur est bien dans la fourchette caractérisant le
secteur étudié.	—
Bibliographie
Bonnet, M. & Schneider, J. 1969. Un exemple d' application de la mé-
thode de Carson-Laplace pour l'étude des relations entre nappe et rivière. Analyse des
variations du niveau du Chari sur la nappe de la depression Tchadienne. Mém.
BRGM. no 76, pp. 53—89.
Boussinesq, J. 1877. Essai sur la théorie des eaux courantes. Mém. Acad.
Sci. Paris, Imprimerie Nationale, Paris, 252 pp.
De Cazenove, E. 1971. Ondes phréatiques sinusoïdales. La Houille Blanche,
no 7, pp. 601—615.
Degallier, R. 1969. Interprétation des variations naturelles du niveau des
nappes. Bull. BRGM. 2ème série, section III. no 2, pp. 7—56.
Degallier, R. 1978. Propagation de fluctuations périodiques dans une nappe
semi-captive non illimitée. Application à la nappe de Moulineaux (Rouen). Rapport
BRGM 78 SGN 319 HYD, 14 p., Orléans.
Drobne, F., Krivic, P., Ravnik, D. & Premru, U. 1977. Študija
o podzemeljski vodi v Tržaško-komenskem krasu, I. faza, leto 1976/77. Rapport GZL,
Ljubljana.
Ferris, J. 1951, Cyclic fluctuations of water level as a basis for determining
aquifer transmissibility. Assoc. Int. Hydrol. Sci., Assem. Gén. Bruxelles, t II. 148 pp.
Garnier, J. C. & Clarion, C. L. 1967.Transmission des ondes de crue
du Verdón à travers le massif d'appui rive droite du futur barrage de Sainte-Croix.
La Houille Blanche, no 1, pp. 61—71.
Transmission des ondes de marée à travers l'aquifère côtier de Kras	325
Krivic, P. 1981. Étude hydrodynamique d'un aquifère karstique côtier: le
Kras de Slovénie, Yougoslavie. Thèse Doct. Ing. Univ. Montpellier. 108 ., Montpellier.
Krivic, P. & Drobne, F. 1980. Hidrogeološke raziskave Tržaško-komen-
skega krasa. Zbornik 6. Jug. simp, hidrogeol. inž. geol., 1, pp. 233—239. Portorož.
Mosetti, F. & Picotti, M. 1960. Actualité sur les recherches de thalas-
sographie dans la mer Adriatique. Boll. Soc. Adr. Sc., 1, pp. 73—77.
Moullard, L., Mijatović, B., Kareh, R. & Massaad,B. 1965.
Exploitation d'une nappe karstique captive à exutoires sous-marins. Problèmes posés
et solution adoptée. Côte libanaise. Actes Coll. Dubrovnik. AIHS-UNESCO, 1, pp.
237—250.
Peaudecerf, P. 1973. Relations entre une nappe alluviale et un cours d'eau
de surface. Les programmes HYGP1-HYGP2. Bull. BRGM. 2ème série, section III.
no 2, pp. 77—85.
Razack, M., Drogue, C., Romariz, C. & Almeida, C. 1980. Étude
de l'effet de marée océanique sur un aquifère carbonate côtier (Miocène de l'Algarve
— Portugal). Journ. HydroL, 45, pp. 57—69.
Todd, D. K. 1964. Ground water Hydrology. John Wiley & Sons, 336 p.. New
York.
T r u p i n , G. 1969. Étude des caractéristiques hydrauliques d'un aquifère captif
par l'observation des effets de pression dus à la marée océanique. Mém. BRGM. no 76.
pp. 359—378.