Cette étude concerne la partie sud-ouest du Massif du Falacro au nord de la Grèce (Macédoine orientale), un . massif de marbres limité par les plaines de Drama et de Nevrokopl et des Intrusions granitiques. Sa base imperméable est formée de schistes et gneiss. La fracturatlon représentée par quatre familles de failles 0° 20°N, 4O°-6O°N, l00°-120°N et 140°-170° a déterminé la morphologie contemporaine du massif. Elle conditionne également les écoulements souterrains et l'apparition des sources en voisinage de ces failles et le creusement des gouffres sur les failles et le creusement de la grotte Mara selon les directions 140°N et 110°N. L'évaluation des volumes d'eau infiltrée à partir des données hydroclimatiques et des volumes de l'alimentation par les gouffres prouve que le massif reçoit une alimentation supplémentaire par ailleurs. Les études hydrochlmiques, hydrodynamiques montrent que les réserves en eau dans ce massif sont faibles et que les circulations se font à deux niveaux : l'un qui est un niveau surtout fissuré dans les marbres blancs et l'autre dans le niveau des marbres graphiteux où il y a augmentation du nombre des chenaux qui aboutissent à un drain principal correspondant à l'axe des gouffres vers la source Mara. Toutes les études précédentes ont démontré l'alimentation de la source Mara par les gouffres et l'existence d'apports supplémentaires tout en permettant de différencier cette source des sources Milopotamos et Drama qui sortent au sud-est du Massif du Falacro.

La détermination des modalités de transfert de l'eau dans un aquifère et la connaissance des volumes d'eau qu'il contient , fait l'objet de nombreuses études en hydrologie. Pour obtenir ces informations, l'hydrogéologue cherche à représenter le transfert de l'eau par une fonction simple qui lui permettrait de prévoir la réponse du milieu souterrain aux stimuli d'entrée , tels que les pluies ou les pollutions, de connaître le temps de séjour moyen de l'eau dans la matrice et par suite de calculer les volumes d'eau contenus dans celle-ci . Dans cette étude, il est proposé de considérer un aquifère comme une "boite noire" dont on ne connait que les fonctions d'entrée et de sortie , et de lui appliquer les méthodes de l'analyse-système qui se sont développées en électronique. A partir de deux hypothèses , vérifiées sur un modèle physique de laboratoire, il est montré qu'un aquifère dont l'entrée est représentée par la surface du bassin hydrologique , peut être assimilée à un système linéaire dont la réponse impulsionnelle est une fonction exponentielle unlque. L'argument de cette exponentielle définit le temps de séjour moyen de l'eau dans le système et le produit du débit moyen des écoulements par ce temps moyen permet de calculer le volume d'eau contenu dans l'aquifère. L'application de cette méthode sur un bassin hydrologique situé à St Genis, Hautes Alpes, s'est avérée satisfaisante . Elle permet de définir une réponse impulsionnelle constituée de deux fonctions exponentielles représentatives de deux systèmes qui composent l'aquifère du bassin et d'en calculer les volumes utiles . L'utilisation des données isotopiques du milieu, nécessairement pondérées par les débits, apparaît comme un supplément d'information indispensable à la bonne détermination de la réponse impulsionnelle et à l'estimation des réserves de l'aquifère.