Η λεκάνη του ποταμού Γλαύκου εκτείνεται επί των δυτικών κλιτύων του Παναχαϊκού, στο ΒΔ/κό τμήμα του νομού Αχαΐας, NA/κά της πόλης των Πατρών καταλαμβάνοντας συνολική έκταση 101,67 Km2 και αποτελεί την κύρια πηγή κάλυψης των υδατικών αναγκών της. Χαρακτηρίζεται από μέτρια ανάπτυξη υδρογραφικού δικτύου δενδριτικού τύπου αποστράγγισης και κατά συνέπεια υψηλή διαπερατότητα των πετρωμάτων και σημαντική κατείσδυση, σε βάρος της επιφανειακής απορροής. Οι σπουδαιότεροι από τους χείμαρρους που διαρρέουν την περιοχή έρευνας είναι ο Ξηροπόταμος, ο Φίλιουρας, η Σούτα, του Βερβενίου και το Ελεκιστριάνικο, οι οποίοι συμβάλουν εντός του Γλαύκου. Το αλπικά πτυχωμένο υπόβαθρο της λεκάνης περιλαμβάνει μια ακολουθία προορογενετικών ιζημάτων (πελαγικοί ασβεστόλιθοι, ραδιολαρίτες, κ.α.) της ζώνης Πίνδου, και μια συνορογενετική κλαστική ακολουθία (φλύσχης), που αποτελείται από εναλλαγές ιλυολίθων, ψαμμιτών και σπανιότερα κροκαλοπαγών. Τα στρώματα της ζώνης Πίνδου αναδύθηκαν με την τελική φάση των πτυχώσεων στο Κάτω Ολιγόκαινο, όπου έλαβε χώρα η επώθηση της Πίνδου, υπό μορφή καλύμματος, πάνω στη ζώνη Γαβρόβου – Τριπόλεως. Τα μέτωπα των επωθήσεων μεταξύ των λεπιών και οι άξονες των πτυχώσεων έχουν διεύθυνση ΒΔ – ΝΑ με κατεύθυνση της κίνησης από τα ανατολικά προς τα δυτικά και άμεση σχέση με την πλειονότητα των καρστικών πηγών. Στην περιοχή του Πατραϊκού κόλπου επικρατεί ένα ρομβοειδές σύστημα ρηγμάτωσης με μεγάλης κλίσης και ΑΒΑ/κής και ΔΒΔ/κής διεύθυνσης κανονικά ρήγματα που δημιουργεί ένα μωσαϊκό κατατεμαχισμένων τεκτονικών μπλοκ και μια γενική ανύψωση της περιοχής που συνεχίζεται μέχρι σήμερα. Τη διάνοιξη των κοιλάδων της ευρύτερης περιοχής έρευνας, η οποία συχνά συνοδεύεται από κατολισθήσεις, την καθορίζουν 4 δέσμες νεοτεκτονικών ρηγμάτων, ΒΔ/κής, ΒΑ/κής, ΔΒΔ/κής, και ΒΒΑ/κής διεύθυνσης. Το Πλειοτεταρτογενές επικάλυμμα που εμφανίζεται στην παράκτια και λοφώδη περιοχή αποτελείται από δύο λιθοστρωματογραφικές ενότητες. Η κατώτερη συνίσταται από αργίλους και άμμους, που αποτέθηκαν σε ένα λιμναίο έως υφάλμυρο περιβάλλον ιζηματογένεσης κατά τη διάρκεια του Πλειόκαινου έως το Κατώτερο Πλειστόκαινο και η ανώτερη από Καλάβριας ηλικίας δελταϊκά και χερσαία κροκαλοπαγή. Γενικά οι προσχώσεις χαρακτηρίζονται από λιθολογική ανομοιομορφία τόσο κατά την κατακόρυφη όσο και κατά την οριζόντια διεύθυνση και το πάχος τους παρουσιάζει βαθμιαία αύξηση από τα βορειότερα και νοτιότερα όρια προς το κεντρικό τμήμα της λεκάνης, όπου και υπερβαίνει τα 200 μ. Το ετήσιο ύψος βροχόπτωσης (1071 mm) στην περιοχή έρευνας είναι αρκετά υψηλό για τα δεδομένα της χώρας μας και κυμαίνεται από 640 – 900 mm στις πεδινές περιοχές, ενώ στις ορεινές περιοχές φθάνει τα 1560 mm. Τα τελευταία χρόνια παρουσιάζεται μια εποχιακή μετακίνηση των βροχοπτώσεων από τον χειμώνα προς την άνοιξη και το φθινόπωρο, η οποία συνεπάγεται απώλεια του διαθέσιμου νερού της ενεργής κατείσδυσης, λόγω αυξημένης εξατμισοδιαπνοής. Από υδρολογικής άποψης η λεκάνη του Γλαύκου δέχεται ετησίως κατά μέσο όρο 108,89 106 m3 νερού από βροχόπτωση, εκ’ των οποίων ποσοστό 48,9 % επανέρχεται στην ατμόσφαιρα με τις διαδικασίες της εξάτμισης και διαπνοής, 28,6 % απορρέει επιφανειακά και 22,5 % κατεισδύει, κυρίως σε σημαντικούς από υδρογεωλογική άποψη σχηματισμούς (ασβεστόλιθους και τεταρτογενείς αποθέσεις). Χαρακτηριστικό υδρογεωλογικό γνώρισμα της περιοχής έρευνας αποτελούν οι ημιαυτόνομες υδρογεωλογικά καρστικές ενότητες, που αναπτύσσονται λόγω της λεπιοειδούς και ρηξιγενούς τεκτονικής στην περιοχή και διαφοροποιούν την κίνηση του υπόγειου νερού. Σχετικά με τις προσχωματικές αποθέσεις του Τεταρτογενούς ο ελεύθερος υδροφόρος ορίζοντας της ανώτερης ζώνης κάμπτεται προοδευτικά και μεταπίπτει σε επάλληλους υπό πίεση υδροφόρους ορίζοντες στην παράκτια ζώνη. Στις αρχές της δεκαετίας του ’90 στην παράκτια περιοχή, λόγω της υπεράντλησης, διαπιστώθηκαν αρνητικές τιμές της πιεζομετρικής επιφάνειας σε απόσταση έως 3 Κm από την ακτή, με συνέπεια την αλμύρινση του προσχωματικού υδροφόρου από τη διείσδυση της θάλασσας. Η μείωση των αντλήσεων στα μέσα της δεκαετίας του ’90 οδήγησαν στην άνοδο των πιεζομετρικών στάθμεων και στην επανεμφάνιση φαινομένων αρτεσιανισμού των γεωτρήσεων της παράκτιας ζώνης. Από την επεξεργασία των αντλητικών δεδομένων διαπιστώθηκε η εμφάνιση υψηλών τιμών του συντελεστή της υδραυλικής αγωγιμότητας τόσο στους προσχωματικούς υδροφόρους της παράκτιας ζώνης, όσο και στους καρστικούς σχηματισμούς. Μικρότερες ταχύτητες κίνησης του υπόγειου νερού εμφανίζουν οι υδροφόροι των αποθέσεων του Νεογενούς, αλλά και των προσχώσεων της νότιας όχθης του ποταμού Γλαύκου σε σχέση με αυτούς της βόρειας όχθης, λόγω της επικράτησης των λεπτομερών έναντι των αδρομερών υλικών. Ο σημαντικότερος παράγοντας τροφοδοσίας του προσχωσιγενή υδροφόρου της λεκάνης του Γλαύκου είναι η πλευρική τροφοδοσία κατά μήκος της κοίτης του ποταμού και εν’ συνεχεία αυτός της πλευρικής τροφοδοσίας από τους ασβεστόλιθους, ενώ λιγότερο σημαντικοί εμφανίζονται οι παράγοντες του απευθείας κατεισδύοντος νερού από τις βροχοπτώσεις και από το περίσσευμα του προς αρδευτική χρήση νερού. Με βάση το ισοζύγιο εισροών – εκροών των προσχωματικών υδροφόρων της λεκάνης του Γλαύκου, για τα υδρολογικά έτη 1999 – 2002 προέκυψε μέσο ετήσιο πλεόνασμα 6,5*106 m3 νερού. Γενικά, αν εξαιρεθεί η παράκτια ζώνη της λεκάνης του Γλαύκου, η ποιότητα του νερού των προσχωσιγενών υδροφόρων, καθώς και των νεογενών αποθέσεων, αλλά και των καρστικών σχηματισμών της λεκάνης, κρίνεται καλή. Στην πλειονότητά τους τα υπόγεια νερά, όπως και τα επιφανειακά και τα πηγαία, εμφάνιζαν τον τύπο Ca-HCO3, που περιλαμβάνει γενικά νερά με καλή τροφοδοσία και ανανέωση κυρίως κατά μήκος των αξόνων αποστράγγισης. Η χωροχρονική μετακίνηση της ζώνης υφαλμύρινσης στην παράκτια ζώνη, την περίοδο 1999 – 2002, διαπιστώθηκε εκτός από την αυξημένη κατανομή διάφορων ιόντων και με την χρησιμοποίηση των συντελεστών S.A.R., Revelle και Schoeller, αλλά και των λόγων rBr-/rCl- και rBr-/rI-. Στη διάρκεια της περιόδου 1999 – 2002, προέκυψε μια αισθητή μείωση της ρύπανσης από τη διείσδυση θαλασσινού νερού στην παράκτια περιοχή. Σημαντικό ρόλο στην επιτάχυνση του ρυθμού της αποκατάστασης των παράκτιων υδροφόρων της λεκάνης του Γλαύκου διαδραμάτισε η πληθώρα γεωτρήσεων αυτόματης ροής στην παράκτια ζώνη, οι οποίες παροχετεύουν περί τα 3,5*106 m3 ρυπασμένου νερού, σε ετήσια βάση, προς τη θάλασσα. Η θέση και η χωροχρονική μετακίνηση της ζώνης υφαλμύρινσης στην παράκτια ζώνη διαπιστώθηκε επίσης από γεωηλεκτρικές έρευνες που πραγματοποιήθηκαν με τις διατάξεις Schlumberger, Wenner – Schlumberger και Pole – Pole. Από τη γεωηλεκτρική έρευνα διαπιστώθηκε ότι η διείσδυση του υφάλμυρου νερού σπανίως εμφανίζεται με τη μορφή μετώπου, αλλά ακολουθεί επιλεκτικές οδούς μέσω των αδρομερέστερων σχηματισμών, ενώ οι λεπτομερείς αργιλικοί σχηματισμοί παίζουν τον ρόλο υδραυλικού φραγμού. Τέλος, με τη μέθοδο της σεισμικής ανάκλασης προσδιορίστηκαν η γεωμετρία και το βάθος συνάντησης του υποβάθρου στην παράκτια ζώνη. Επισημαίνεται η ανάγκη τεχνητού εμπλουτισμού του προσχωσιγενή υδροφόρου ορίζοντα με τις χειμερινές απορροές του Γλαύκου για την ταχύτερη αποκατάσταση της ποιότητας του νερού και την αύξηση του υδατικού δυναμικού της λεκάνης. Εφαρμόστηκε η μέθοδος DRASTIC για την εκτίμηση της τρωτότητας των υπόγειων νερών απέναντι σε εξωτερική ρύπανση, όπου με βάση τις τιμές του δείκτη DRASTIC κατασκευάστηκε χάρτης τρωτότητας της περιοχής. Περιοχές υψηλής τρωτότητας εντοπίζονται στο παράκτιο τμήμα του υδροφορέα, ενώ αντίθετα περιοχές μέσης τρωτότητας στο Α/κό τμήμα του υδροφορέα. Τέλος, κατασκευάστηκε μαθηματικό ομοίωμα των υδροφορέων της προσχωματικής λεκάνης του Γλαύκου με την εφαρμογή του κώδικα MODFLOW παρέχοντας έτσι ένα όσο το δυνατόν πιο αξιόπιστο πρότυπο διαχείρισης του υπόγειου υδατικού δυναμικού της περιοχή έρευνας. / The Glafkos’ river basin is extended on the west side of the Panaxaikon mountain in the northwest area of the prefecture of Achaia, southeast of the city of Patras, occupying a total extent of 101,67 Km2 and constituting the main source of Patras’ water needs. It is characterized by the mediocre growth of the hydrographical network (a dendrite type of draining) and consequently the high permeability of rocks and the important infiltration at the expense of the runoff. The most important torrents that flow through the region of research are Xiropotamos, Filiouras, Souta, Verveniou and Elekistrianiko, which contribute inside the Glafkos River. The Alpine folded basin’s basement includes a sequence of pre-orogenetic sediments (marine limestone, radiolarites, etc) of the Pindos zone, and a clastic orogenetic sequence (flysch), which is constituted by alternations of silt rocks, sandstones and rarely conglomerates. The layers of the Pindos zone emerged at the final phase of the folding tectonics in the Lower Oligocene, where the Pindos’ thrust occurred in the form of a nappe, above the Gavrovo – Tripolitsa zone. The foreheads of the thrusts between the slices as well as the axes of folds have a NW – SE direction, with the direction of movement from east to west, and have a direct relation to the majority of the karstic springs. In the region of Patraikos Gulf, what prevails is a rhomboid faulting system with high slope and ENE, WNW directed normal faults that create a mosaic of fragmentized tectonic blocks and a general rising of region that continues until today. The opening up of the valleys in the wider region of research, which is often accompanied by landslides, is determined from 4 bunches of neotectonic faults directed NW, NE, WNW and NNE. The Plio-quaternary cover that appears in the coastal and hilly region is constituted by two stratigraphical units. The lower one consists of clays and sands, which were deposited in a lacustrine to brackish environment of sedimentation during the Pliocene to the Lower Pleistocene and the upper one consists of Calabria dating from the deltaic and land conglomerates. Generally the alluvium deposits are characterized by lithological heterogeneity both at the vertical and the horizontal direction and their thickness presents a gradual increase from the northerner and southerner limits to the central department of the basin, which exceeds 200 m. The annual height of rainfall (1071 mms) in the region of research is quite high for our country’s data and fluctuates from 640 to 900 mms in flat regions, whereas in the mountainous regions it reaches 1560 mms. During the last years, what takes place is a seasonal displacement of rainfall from winter to spring and autumn, which leads to the loss of the available water of the active infiltration because of increased evaporotranspiration. From a hydrologic aspect, the Glafkos basin receives annually, on average, 108,89 106 m3 of the water from rainfall, from which 48,9 % comes back to atmosphere with the processes of evaporation and transpiration, 28,6 % rises on the surface and 22,5 % infiltrates, mainly in important, from a hydrogeological aspect, formations (limestones and quaternary depositions). The semi-independent hydrogeologically karstic units, which are developed because of the laminated and faulty tectonics in the region and differentiate the movement of groundwater, constitute a characteristic hydrogeological feature of the region of research. In relation to the alluvium deposits, the unconfined aquifer of the superior area is progressively bent and falls into successive under pressure (confined) aquifers in the coastal area. At the beginning of the ‘90s in the coastal region, because of the overpumping, negative values of the piezometric surface in a distance up to 3 Km from the coastline were realised, having as a consequence the salinity of coastal alluvium aquifer from the seawater intrusion. The reduction of pumping in the middle of the ‘90s led to the rise of piezometric levels and to the reappearance of the phenomena of artesian wells in the coastal area. Due to the elaboration of the pumping data, what was realised is the appearance of high values of the hydraulic conductivity’s factor in the alluvium aquifers of the coastal zone as well as in the karstic aquifers. The aquifers of the Neogene formations as well as the alluvium aquifers of the southern bank of the river Glafkos present a smaller speed of movement of groundwater as compared to those of the northern bank, because of the predominance of the fine-grained clastics. The most important factor of the recharge of the Glafkos’ basin alluvium aquifer is basically the costal recharge from the river bank and then follows the inductive recharge from the limestones, whereas the direct infiltrated water from the rainfall and from the surplus of the irrigation water seem to be the least important factors. Based on the balance of the inflows – outflows of the Glafkos’ basin alluvium aquifer, an average annual surplus of 6,5*106 m3 of water has resulted, during the hydrologic years 1999 –2002. Generally, the quality of water in the alluvium aquifers as well as in the aquifers of Neogene formations, but also that in the karstic aquifers is considered to be good, if the coastal area of the Glafkos’ basin is excluded. Most of the groundwater, as well as the surface and spring water, presented the geochemical type Ca-HCO3, which generally includes water of good recharge and renewal, mainly along the drainage axes. The time-space movement of the brackish zone in the coastal area, in the period 1999 –2002, was realised not only by the increased distribution of various ions but also by the utilisation of factors S. A. R., Revelle and Schoeller, as well as the ratios rBr-/rCl- and rBr-/rI-. During the period 1999 – 2002, a perceptible reduction of the pollution from the saline intrusion in the coastal region took place. The abundance of wells of automatic flow into coastal area, which drain 3,5*106 m3 of brackish water, annually, to the sea played an important role in the acceleration of the restoration of the coastal basin aquifers of the river Glafkos. The position and the time-space movement of the brackish zone in the coastal area was realised also by geoelectrical researches that took place with the dispositions Schlumberger, Wenner – Schlumberger and Pole – Pole. Geoelectrical research revealed that seawater intrusion rarely appears in the form of a forehead, but follows selective paths via the coarse-grained clastics, while the fine-grained clastic formations play the role of a hydraulic barrier. Finally, the geometry and the depth of the basement in the coastal area were determined with the application of the seismic reflection method. The need of artificial enrichment of the alluvium aquifer from the wintry flows of the Glafkos River is pointed out, for a more rapid restoration of the water quality and the increase of the Glafkos’ basin water potential. The DRASTIC method was applied for the estimation of the vulnerability of the groundwater towards exterior pollution, where based on the DRASTIC indicator’s values, a map of the vulnerability of the region of research was constructed. Areas of high vulnerability are located in the coastal part of the aquifer, while areas of medium vulnerability are located in the eastern part of the basin. Finally, a mathematic model of the Glafkos’ basin alluvium aquifers was constructed with the application of code MODFLOW, thus providing an as much as possible reliable model of the management of the Glafkos’ basin groundwater potential.

Στην παρούσα εργασία έγινε προσπάθεια να διερευνηθούν οι υδρογεωλογικές συνθήκες των καρστικών υδροφόρων οριζόντων του κεντρικού τμήματος του νομού Αχαΐας και να οριοθετηθούν οι κύριες υδρογεωλογικές λεκάνες ή ενότητες της περιοχής. Στα πλαίσια της έρευνας αυτής εκτελέσθηκαν οι παρακάτω εργασίες: α) Γεωλογική αναγνώριση της περιοχής. β) Ανάλυση του υδρογραφικού δικτύου και του μορφολογικού αναγλύφου της περιοχής έρευνας. γ) Χάραξη γεωλογικών τομών που καλύπτουν όλη την περιοχή μελέτης με σκοπό τη διερεύνηση του υδρογεωλογικού καθεστώτος στην περιοχή. δ) Αναγνώριση των κυριοτέρων πηγών στην ύπαιθρο και σχεδιασμός σκαριφημάτων για κάθε μία από αυτές. ε) Χάραξη υδρογεωλογικών λεκανών για την περιοχή μελέτης που αποστραγγίζονται από τις αντίστοιχες πηγές. ζ) Επεξεργασία των υδροχημικών δεδομένων στα νερά των πηγών. Σύμφωνα με τα αποτελέσματα που προέκυψαν από την παρούσα έρευνα φαίνεται ότι η έντονη τεκτονική καταπόνηση που έχουν δεχθεί οι γεωλογικοί σχηματισμοί κατά το παρελθόν παίζει σημαντικό ρόλο στη διαμόρφωση των υδρογεωλογικών λεκανών της περιοχής, καθώς τα ρήγματα ελέγχουν την εξάπλωσή τους. Ακόμη, τα παράλληλα τεκτονικά λέπια δημιουργούν υδρογεωλογικές λεκάνες διεύθυνσης περίπου Β-Ν. Τοπικό επίπεδο βάσης αποτελούν πολλές φορές οι πυριτικές παρενστρώσεις των ανωκρητιδικών ασβεστολίθων, αφού συντελούν στο σχηματισμό μικρής σημασίας υδροφόρων οριζόντων. Στα μέτωπα των επωθήσεων παρουσιάζεται μεγαλύτερο ποσοστό καρστοποίησης των ασβεστολίθων (π.χ. Καλούσι, Πλατανόβρυση) λόγω μεγαλύτερης καταπόνησης των πετρωμάτων στις περιοχές αυτές. Τα διαγράμματα Piper δίνουν γενικά τύπο νερού Ca – HCO3 που δείχνει νερά που φιλοξενούνται σε καθαρά ασβεστολιθικούς υδροφόρους. Τέλος, η πλήρης ή μερική υδρομάστευση των πηγών συντελεί στην αλλαγή του υδραυλικού καθεστώτος των πηγών, με αποτέλεσμα μερικές από αυτές να μην λειτουργούν πλέον, ενώ η ρύπανση του νερού των πηγών σε ορισμένες θέσεις οφείλεται καθαρά σε ανθρωπογενείς δραστηριότητες. / This paper presents the results of the hydrogeological research in the karst aquifers conducted in the central part of the Achaia province, as well as the limitation of the hydrogeological basins. During the research, the tasks that were performed were the following: a) Geological research of the whole area. b) Geomorphological and hydrographical analysis of the research area. c) Side elevations in order to define the hydrogeological conditions in the exploration area. d) Location of the springs and drawing side elevations for each one of them. e) Drawing of the hydrogeological basins that are dewatered of the springs. f) Elaboration of the hydrochemical data of the spring waters. According to the results that came up at the end of this research, it seems that the tectonic pressure to the rocks and the sediments of the area has infected the development of the hydrogeological basins because the faults usually control their diffusion. Furthermore, the parallel tectonic flakes form hydrogeological basins with a North-South orientation. Local base-level of the karst aquifers are usually the clay intrusions inside the limestone layers, that form aquifers of small importance. In the areas near the thrust, there is bigger solution in the karst aquifers (e.g. Platanovrisi, Kalousi etc.), because of the amount of pressure that affected those areas. The Piper diagrams suggest a water type of Ca-HCO3 that shows water from limestone aquifers. Finally, the use of the spring waters in many cases is so big that it affects the level of the underground water, while the contamination of the water seems to be imputed to the stock farming in some places.

Deux expériences ont été menées sur la plate-forme expérimentale "SCERES" afin d'évaluer les concentrations et les flux de vapeurs de TCE dans SCERES en présence de deux dalles de béton fissurées installées, l'une après l'autre, à la surface de SCERES. Cet aquifère poreux est un milieu hétérogène de grande échelle (25 x 12 x 3 m3). Les résultats ont montré que le panache de vapeur de TCE couvre la plupart du bassin au bout de 3 semaines depuis la création de la zone source de TCE dans le sous sol. L'hétérogénéité du site SCERES a engendrée une distribution verticale non uniforme de la concentration de vapeurs de TCE. La simulation du panache de vapeur dans SCERES a été effectuée au moyen du code de calcul multiphasique "SIMUSCOPP". La présence sur SCERES de la dalle de béton, un milieu peu perméable et peu diffusif, a constitué une "barrière" en vue du transfert de vapeurs de TCE vers l'interface dalle/atmosphère. Afin de mieux quantifier le flux de vapeurs à travers la dalle de béton, une étude de coefficient de diffusion et de perméabilité des deux dalles a été réalisée. Un mouvement vertical ascendant du toit de la nappe a généré un fort gradient de pression motrice de l'air du sol. Ceci a engendré une forte augmentation des flux de vapeurs à l'interface sol/atmosphère. La quantification de ces flux de vapeurs a été effectuée à l'aide d'une solution semi analytique basée sur la loi de Fick et la loi de Darcy en tenant compte à la fois de l'effet de gradient de pression motrice et l'effet de densité de vapeurs sur le transfert de vapeurs vers la surface du sol. L'intrusion de vapeurs de TCE dans le bâtiment modèle, installé sur la dalle de béton, a été générée par une mise en dépression dans ce dernier. Ce qui a fait augmenter la concentration de vapeurs de TCE sous la dalle ainsi dans le bâtiment. La simulation de l'intrusion de vapeurs dans l'air intérieur de bâtiment a été réalisée par l'intermédiaire du code de calcul multiphysics "COMSOL", avec lequel nous avons démontré l'évolution de la concentration de vapeurs obtenues expérimentalement dans le bâtiment et qui dépend directement de la variation spatio-temporelle du flux massique à travers la dalle. / Two experiments were conducted on the experimental platform "SCERES" to assess the TCE vapour concentrations and fluxes in SCERES with two concrete slabs installed, one after the other, on the ground surface. This artificial aquifer is a large scale (25 x 12 x 3 m3) heterogeneous porous medium. The results showed that the TCE vapour plume covers most of the basin 3 weeks after the creation of the TCE source area in unsaturated zone. The heterogeneity of SCERES has generated a non uniform vertical distribution of the TCE vapour concentration. Simulation of vapour plume in SCERES was carried out by the multiphase code "SIMUSCOPP".The presence in SCERES of a low permeability and low diffusive medium compared to the sand in the basin,as a concrete slab, constituted a "barrier" for the transfer of TCE vapour to the interface concrete slab / atmosphere. To better quantify the TOE fluxes through the concrete slab, a study of diffusion coefficient and permeability of both concrete slabs was done. An upward vertical movement of the water table has generated a strong soil air pressure driving gradient, which led to a strong increase in the TCE vapour concentrations near the surface which has increased the vapour fluxes at the interface soil / atmosphere. Quantification of vapour fluxes at the interfaces soil / atmosphere and concrete slab / atmosphere was performed using a semi analytical approach based onFick's and Darcy's laws by taking into account both the effect of the driving pressure gradient and the effect of density vapour on the vapour transfer towards the soil surface.The intrusion of TCE vapours into the model building installed on the concrete slab was generated by creating a vacuum. The results showed that, during the TCE vapour suction from the model building, the concentration of TCE vapours under the concrete slab and in the building increases. Simulation of vapour intrusion into indoor air was done by the computational Multiphysics code "COMSOL", allowing simulation of the evolution of the vapour concentration obtained experimentally in the building. Il was shown that they depend directly on the spatial-temporal variation of the mass flux through the slab.

L’Atlas saharien sud oranais où les Monts des Ksour forment une région montagneuse constituée de chaînonsallongés or ientés SO–NE. Les dépressions crétacées sont représentées par des synclinaux à fond plat surplombés par desanticlinaux pr incipalement jurassiques. Les formations gréseuses crétacées appelées « Continental Intercalaire » de la régiond’Aïn Séfra renferment une impor tante ressource en eau destinée à l’alimentation en eau potable et à l’ir r igation. L’étudehydrogéologique et l’identi fication des aquifères de la région montrent l’impor tante interconnexion ent re l’aquifère jurassiqueet le CI. La recharge naturelle des aquifères sous un climat semi-aride est le terme du bilan hydrique le plus délicat àdéterminer . La méthodologie proposée est une approche par interpolations spatiales et temporelles des données deprécipitations journalières et de températures moyennes mensuelles mesurées à une station, durant trois années repères(sèche, moyenne, humide). Elle est très « pratique » dans des r égions où les données climatiques sont limi tées et discontinuesdans l’espace et dans le temps. Son application nécessite des connaissances sur la nature du sol, la surface et l’altitude desimpluviums. Cette étude nous a permis d’estimer , à par t ir des gradients pluviométr ique et thermique, la lame d’eau tombéesur les affleurements jurassiques à une altitude moyenne de 1600 m, et les valeurs de l’évapotranspiration dans tout le bassinétudié. La pluie efficace calculée, et le débit d’alimentation de l’aquifère, par son impluvium calcaréo-gréseux en altitude etgréseux dans les dépressions, montrent que l’aquifère du Continental Intercalaire est alimenté pr incipalement par lesaffleurements du Jurassique surplombant les dépressions gréseuses, durant les années pluvieuses où les précipitationsdépassent 300 mm par an. Les caractér istiques chimiques des eaux ont fait l’objet de différentes analyses statistiques afin detenter d’expliquer l’or igine d’une minéralisation assez élevée (ent re 2000 et 4000 mg/ l). Dans une première approche,l’emploi de l’analyse en composantes pr incipales (ACP) met en évidence trois fami lles d’eau cor respondant aux t rois niveauxaquifères exploités. Dans une seconde étape, l’analyse factor ielle discr iminante croisée (AFD) caractér ise l’impor tance et ladominance des paramètres responsables de cette minéralisation, en relation avec le cadre géostructural du systèmehydrogéologique des gouttières synclinales d’Aïn Séfra ; les contacts anormaux entre les formations gréseuses du ContinentalIntercalaire et les structures diapir iques argilo-gypso-salines du Tr ias forment les pr incipales zones de contamination. Aprèsavoi r estimé la recharge efficace (qui var ie entre 0 pour une année sèche et 684 mm/ an une année humide) et identi fié lesécoulements souter rains et pour prévoir le compor tement hydrodynamique des gouttières synclinales d’Aïn Séfra, i l estpr imordial de réaliser un modèle mathématique. Nous n’avons pu réaliser que le modèle de l’unité hydrogéologique de GaretEl Kheil, vu l’indisponibilité des données sur les deux aut res unités (El Hanjir -Aïn Séfra, Aïn Séfra-Tiout). Nous avons choisi uncalcul par différences finies à l’aide d’un code simple et pratique « ASMWIN ». Afin de déterminer l’impact du changementclimatique et la surexploitation, deux scénarios d’exploitation ont été proposés : une recharge réduite de moitié (-50% : 810 6m3/ an au lieu de 16 10 6) et un débit d’exploitation doublé (+100% : 2.72 106m3/ an au lieu de 1.36 10 6). Ils nous ontrévélé que la nappe est plus vulnérable aux changements climatiques qu’à la surexploitation, ce qui nous a permis de déduireque la nappe est encore en sous-exploitée. / The Saharan Atlas Mountains, South of Oran includes the ‘Ksour Ridge’, a mountainous region compr ising SW-NEranges. Cretaceous depressions are represented by flat-bottomed synclines overhung by anticlines, mainly Jurassic. Thesandstone formations, so-called Continental Intercalary, of the Ain Sefra region contain an unconfined significant waterresource for dr inking water and ir r igation. The hydrogeological study and identification of aquifers prove the impor tantinterconnection between the Continental Intercalary and the Jurassic aquifers. In a recharge area character ized by a semiaridor an arid climate, therefore a low annual rainfall, effective rainfall is the hardiest item to estimate, as it was demonstrated byclassical monthly hydrological balance computations. The methodology proposed in this work relies on spatial and temporalinterpolations of scarce climatic data. This methodology can be used to determine the recharge flow to aquifers of these zones.The Ain Sefra’s syncline includes the sandstone aquifers of Continental Intercalary and of the Plio-Quaternary layers. Threeyears, with contrasting annual rainfalls, are chosen in this application: 1983 as dry year , 2006 as average year , and 2008 ashumid one. This application requires information about the soil nature, thickness and porosity, the sur face and the altitude ofthe recharge areas. The rainfall and the recharge flow to aquifers indicate that the pr incipal recharge areas are the Jurassicoutcrops, mainly dur ing the humid years, when the annual rainfall exceeds 300 mm. The chemical character istics of waterhave been submitted to different statistical analyses. As a fi rst approach, a pr incipal component analysis (PCA) displayed threeclusters of water cor responding to the three aquifer layers which are exploited. Discr iminant factor analysis (DFA) was usedto character ize the impor tance and to rank the parameters responsible for this high mineralization (between 2000 and 4000mg/ l), namely: system character istics (lithology), structural geology (fault and Tr iassic diapirs). After estimating the effectiverecharge (ranging between 0 for a dry year, and 684 mm/ year for a humid one) and identifying groundwater flow, and inorder to predict the hydrodynamic behavior of Ain Sefra’s syncline, it is essential to implement a mathematical model. Wewere only able to establish the model of the Garet El Kheil hydrogeological unit, given the data unavailability on the other twounits (El Hanjir -Ain Sefra, Ain Sefra-Tiout). We chose a finite difference method in permanent mode to improve knowledge ofthe hydraulic character istics of aquifer , and to evaluate a complete hydrological balance, using the simple and convenient"ASMWIN" code. Exploitation of the model made it possible to determine the impact of the pumped yields, and the weathertrends on the piezometr ic evolution of the aquifer . Two operating scenar ios have been proposed: halving recharge (-50%: 810 6 m 3 / year instead of 16 10 6) and doubling the well’s yield (+100 %: 2.72 10 6 m 3 / year , instead of 1.36 10 6). Weconclude that the groundwater resource is more vulnerable to the climate change than intensive use of aquifers, meaning thatthe aqui fer is underexploited at present.

La nappe superficielle de la plaine d’Avignon, située entre le Rhône et la Durance, présente les caractéristiques d’une nappe péri urbaine, avec un secteur agricole étendu, mais dont la surface se réduit rapidement sous l’effet de l’extension urbaine et du développement économique. Peu profonde et proche de la surface, cette nappe est directement sous l’influence des eaux de surface. En particulier, l’eau de la Durance participe fortement à la recharge soit directement soit indirectement via les nombreux canaux d’irrigation et donc localement de manière décalée par rapport au régime hydrologique méditerranéen. L’occupation des sols, et en particulier l’urbanisation, joue donc un rôle important sur le fonctionnement hydrodynamique. En particulier le recul des terres agricoles et la diminution parallèle des volumes d’irrigation tendent à réduire la recharge de la nappe.Pour essayer d’anticiper les conséquences de cette évolution, la présente étude vise à proposer un modèle de fonctionnement hydrogéologique de la zone intégrant le rôle du sol et de son usage. Il fallait pour cela compléter les connaissances sur le fonctionnement de la nappe préciser les conditions aux limites et établir un bilan en eau au pas de temps mensuel. Ce travail repose sur une démarche multi-approches utilisant, en particulier, la complémentarité entre le suivi piézométrique et la technique du traçage naturel de l'eau (traceurs chimiques et isotopiques). Les différentes approches ont permis de localiser et de quantifier les trois principales composantes de la recharge : la Durance ou la recharge latérale, l’irrigation et la pluie ou la recharge verticale. Cette dernière est contrôlée surtout par la nature et l’usage du sol. Très clairement, la majeure partie de cette recharge se trouve dans des zones encore irriguées qui apportent un volume moyen annuel de 19 Mm3. Avec une moyenne annuelle de 135 mm, la recharge par la pluie joue un rôle secondaire. Son impact sur la nappe est visible en dehors des zones irriguées en particulier en amont de la plaine sur les versants des collines qui bordent la plaine à l’Est. Cet apport de la pluie n’est plus visible dans la partie médiane, fortement diluée par les eaux d’irrigation. Plus à l’ouest en zone urbaine, l’effet de la pluie est réduit et discontinu à cause des surfaces imperméabilisées ; il subit en outre une dilution par l’apport des eaux de la Durance. L’importance de la recharge par la Durance a été confirmée par une modélisation hydrodynamique et du transfert du signal isotopique. Elle constitue l’entrée principale du système avec un volume moyen de 70 Mm3/an. Les résultats de différents scénarios du modèle numérique valident les ordres de grandeur de différentes composantes et les caractéristiques de l’écoulement (direction, vitesse…) et montrent que dans le cas d’un abandon de l’irrigation, la baisse piézométrique serait de l’ordre de 3 m sous les périmètres irrigués. La nappe retrouverait alors les fluctuations piézométriques naturelles sous la dépendance de la pluie d’une part et de la Durance d’autre part. En fonction de son régime futur, l’alimentation par la Durance pourrait augmenter et compenser une partie de la baisse vers l’intérieur de la plaine. / The groundwater system of the alluvial plain of Avignon, located between the Rhone and the Durance rivers, is characterized by an expanded agricultural area, whose surface is rapidly shrinking due to urban sprawl and economic development. This shallow aquifer is directly under the influence of surface water. In particular, the Durance water is heavily involved in the recharge process either directly or indirectly via the numerous irrigation canals and therefore with a time shift according to the Mediterranean hydrological regime. Land use, particularly urbanization, plays an important role on the hydrodynamic occurrence. In particular, the loss of agricultural land and accordingly the reduction in irrigation volumes tend to reduce groundwater recharge. To anticipate the consequences of these developments, this study aimed at proposing a hydrogeological model of the area including the role of land use. First of all, This required a better understanding of the aquifer system, namely specify boundary conditions and establish a monthly water balance An experimental investigation was carried out using different approaches including water level monitoring and environmental tracing of water (chemical and isotopic tracers). This work enabled us locating and quantifying the three main components of recharge: Durance or lateral recharge, irrigation and rainfall or vertical recharge. This recharge proved to be controlled primarily by the soil nature and land use. Clearly, much of this recharge occur in irrigated areas which provide an average annual volume of 19 millions m3. With an annual average of 135 mm, recharge by the rainfall plays a secondary role. Its impact on the groundwater was observed out of the irrigated lands especially, on the eastern boundaries of the plain. This contribution of rainfall was no longer detected downwards, where irrigation water was getting predominant in the recharge process. In urban areas down west, the effect of rainfall drops owing to the increase of impervious surfaces and the rising of Durance influence.. The importance of the recharge from the Durance was confirmed by hydrodynamic and isotope transfer modeling. It is the main input to the system with an average volume of 70 millions m3 / year.

L'importance des enjeux liés à la qualité et à la connaissance de l'eau souterraine nécessite de caractériser le fonctionnement des nappes d’eaux souterraines. En hydrogéologie, les méthodes qui permettent de réaliser le modèle conceptuel d’un aquifère reposent généralement sur l’observation et la caractérisation du milieu, à partir notamment de la mise en place de forages et la réalisation de pompages d’essai. Cependant, ces méthodes présentent des limites pour caractériser la variabilité spatiale des aquifères hétérogènes. Dans cette thèse, nous nous sommes intéressés au développement d’approches expérimentales qui utilisent des informations issues de traçage thermique ou salin pour la caractérisation des paramètres hydrauliques des milieux hétérogènes poreux. Ce choix expérimental a été justifié 1) par les difficultés de prise en compte des hétérogénéités hydrauliques dans les milieux poreux naturels et 2) pour une meilleure compréhension des phénomènes mis en jeu en situation contrôlée. Dans les travaux réalisés, nous avons tout d’abord proposé d’utiliser des mesures thermiques ponctuelles, puis obtenues dans l'infrarouge thermique pour estimer les paramètres hydrauliques en milieux poreux hétérogènes. Cependant, les méthodes de traçage thermiques étant limitées dans les zones faiblement perméables, nous avons fait évoluer notre approche en la combinant avec une méthode géophysique (Potentiel Spontané). Cette méthode, sensible aux écoulements des fluides et aux processus thermoélectrique et électrochimique, nous a permis de collecter des données temporelles intéressantes pour surveiller des variations thermiques ou saline provoquées dans les milieux. Les mesures obtenues ont été par la suite utilisées dans des algorithmes d’inversion pour estimer les distributions spatiales des propriétés hydrauliques. Les résultats obtenus nous ont permis de montrer l’efficacité de ces approches pour la caractérisation des milieux hétérogènes et par conséquent pour la modélisation des écoulements des fluides dans ces milieux. / The importance of issues related to the quality of groundwater requires characterizing the operation of groundwater aquifers. Hydrogeological methods used to carry out conceptual model of an aquifer are generally based on the observation and characterization of the medium, by using boreholes and carrying out of pumping test operations. However, these methods present some limitations in characterizing the spatial variability of heterogeneous aquifers. In this thesis, we develop experimental approaches that use information from thermal or saline tracing for the characterization of hydraulic parameters in heterogeneous porous media. The choice of using experimental approaches was justified 1) by difficulties of taking into account hydraulic heterogeneities in natural porous media and 2) for a better understanding of the phenomena involved in a controlled situation. In this work, we first proposed the use of punctual thermal measurements, then measurements obtained by using thermal infrared to estimate hydraulic parameters in heterogeneous porous media. However, as thermal tracing methods are limited in low permeability areas, we improved our approach by combining it with a geophysical method (Spontaneous Potential). This method, sensitive to fluid flows and thermoelectric and electrochemical processes, has allowed us to collect interesting temporal data to monitor thermal or saline variations in the media. The measurements obtained were then used in inversion algorithms to estimate the spatial distributions of hydraulic properties. The results obtained allowed to demonstrate the effectiveness of these approaches for the characterization of heterogeneous media and consequently for the modelling of fluid flows.

Près de 25% de la population mondiale est alimentée par de l’eau en provenance d’aquifères karstiques. La compréhension et la protection de ces derniers apparait donc comme essentielle dans le cadre d’une augmentation des besoins en eau potable. De plus, une contamination des forages d'alimentation en eau potable par une eau turbide peut s'avérer fortement dommageable car entrainant une possible contamination des populations desservies. Dans le cas de la Normandie, des coupures régulières son nécessaires afin de préserver la santé des habitants. La modélisation et la prédiction des augmentations de turbidité apparaissent comme un travail difficile du fait des nombreux phénomènes et paramètres régissant la turbidité ainsi que la non-linéarité de la réponse entre les précipitations et la turbidité. Peu de modèles à l'heure actuelle ont été proposés pour représenter la relation liant la turbidité avec les précipitations. C'est ainsi, en s'intéressant au forage AEP d'Yport responsable de l'alimentation en eau potable de la ville du Havre, que nous proposons une application des réseaux de neurones pour la prévision de la turbidité. Durant les travaux de thèse, nous avons mis en avant la nécessité d'effectuer des campagnes d'échantillonnages des produits phytosanitaires afin de permettre l'identification des éventuels proxies des produits phytosanitaires tel que la turbidité, les précipitations ou bien la conductivité. Par la suite, les travaux effectués dans cette thèse nous ont permis (i) de monter que les modèles par réseaux de neurones permettent de prévoir à 12h et 24h les variations de turbidité, (ii) de tester plusieurs voies d'amélioration de ces modèles, (iii) d'intégrer l'analyse multirésolution aux modèles par réseaux de neurones et pour finir (iiii) d'identifier un semi proxy des contaminations en produits phytosanitaires. / Approximately 25% of the world's population is supplied by water from karstic aquifers. The understanding and protection of these appears to be essential in the context of drinking water needs increasing. In addition, contamination of drinking water by turbid water can be highly damaging by resulting in possible contamination of the served populations. In the case of Normandy, regular drinking water cut-off are necessary to preserve the health of the inhabitants. The modeling and prediction of turbidity event appears as a challenging work because of the number of phenomenon and parameters involves in turbidity variation as well as the non-linearity of the link between rainfall and turbidity. Actually, few models have been proposed to represent the relationship between turbidity and rainfall. In this context, by focusing on Yport's pumping well which is responsible for Half of Le Havre city drinking water supply, we propose an application of neural networks for turbidity prediction. During this thesis work, we emphasized the need to carry out sampling campaigns for phytosanitary products to enable the identification of possible phytosanitary product proxies such as turbidity, rainfall or conductivity. Subsequently, the work carried out in this thesis enabled us to (i) designed neural network models allow to predict at 12h and 24h the turbidity variations, (ii) test several ways to improve these models, (iii) integrate multiresolution analysis into neural networks models and finally (iiii) identify a semi proxy for phytosanitary product contamination.

Depuis 2006, l’arrêt des pompages d’exhaure dans le Bassin Houiller Lorrain (France) a conduit à l’ennoyage des travaux miniers abandonnés, avec pour conséquence la mise en place d’un nouvel équilibre hydrodynamique régional. De récentes recherches portant sur l’exploitation de la chaleur des réservoirs ennoyés ont suscité de nouvelles interrogations, auxquelles nous nous sommes proposé de répondre. Notre travail avait pour objectif de chercher à comprendre le comportement thermo-hydrodynamique de l’eau de mine au sein d’un système en cours d’ennoyage ou récemment ennoyé. Dans un premier temps, les contextes géographique, géologique et hydrogéologique du Bassin Houiller Lorrain ont fait l’objet d’une synthèse, et une zone d’étude a été choisie. Dans un second temps, des profils de température et de conductivité électrique, complétés par des jaugeages, ont été mesurés dans d’anciens puits de mine du Bassin Houiller Lorrain, offrant une meilleure compréhension du comportement de l’eau à l’échelle d’un ouvrage profond. À partir de l’analyse de ces données, un modèle thermo-hydrodynamique et des simulations numériques ont pu être réalisés à cette échelle. Les résultats permettent d’expliquer les phénomènes observés. Leur stabilité est également étudiée. Dans un troisième temps, un modèle spatialisé maillé a été construit pour aborder la problématique du comportement hydrodynamique d’un réservoir minier entier. La montée de l’eau observée a été correctement reproduite ; le modèle peut par ailleurs être utilisé de façon prédictive pour la période suivant l’ennoyage. Plusieurs outils ont été testés, améliorés ou développés afin de faciliter l’étude des réservoirs ennoyés, notamment concernant l’homogénéisation tridimensionnelle des conductivités hydrauliques et le couplage d’un modèle spatialisé maillé avec un réseau de drains-conduits / Since 2006, cessation of dewatering in Lorraine Coal Basin (France) led to the flooding of abandoned mines, resulting in a new hydrodynamical balance in the area. Recent researches concerning geothermal exploitation of flooded reservoirs raised new questions, which we propose to answer. Our work aimed to understand the thermo-hydrodynamical behaviour of mine water in a flooding or flooded system. Firstly, we synthetized the geographical, geological and hydrogeological contexts of the Lorraine Coal Basin, and we chose a specific area for our studies. Secondly, temperature and electric conductivity log profiles were measured in old pits of the Lorraine Coal Basin, giving a better understanding of the water behaviour at a deep mineshaft scale. We were able to build a thermo-hydrodynamical model and simulate water behaviour at this scale. Flow regime stability is also studied. Thirdly, a hydrodynamical spatialized meshed model was realized to study the hydrodynamical behaviour of a mine reservoir as a whole. Observed water-table rise was correctly reproduced : moreover, the model can be used in a predictive way after the flooding. Several tools were tested, improved or developed to ease the study of flooded reservoirs, as three-dimensional upscaling of hydraulic conductivities and a coupled spatialized meshed model with a pipe network

Les eaux souterraines constituent une réserve d’eau potable non négligeable, leur alimentation se fait en majeure partie par les précipitations, appelée recharge des nappes. Du fait de leur grande importance, la compréhension du fonctionnement de ces ressources en eau est plus que jamais indispensable. Celle-ci passe par l’élaboration de modèles mathématiques. Ces outils nous offrent une meilleure appréhension et une bonne prévision des phénomènes physiques. Les systèmes hydrogéologiques sont généralement très complexes et caractérisés par des dynamiques hydriques très variables dans le temps et dans l’espace. Cette complexité a attiré l’attention de nombreux hydrogéologues et un grand nombre de modèles très sophistiqués a été développé afin de décrire ces systèmes avec précision. Cependant, la prise en compte de la recharge de ces réservoirs reste toujours un défi dans la modélisation hydrogéologique. En effet, le plus souvent, les modèles hydrogéologiques simulent l’écoulement dans la nappe tout en considérant la recharge comme une constante sur le domaine et indépendante du système. De plus, elle est souvent calculée de façon simplifiée. Or, la recharge traduisant la quantité des précipitations atteignant les nappes est une composante hydrologique complexe et variable car elle interagit avec les nappes et dépend des conditions climatiques, du couvert végétal et du transfert de l’eau dans le sol. Ce présent travail vise à intégrer cette recharge variable et complexe aux modèles hydrogéologiques. À cet effet, un modèle couplé a été développé. Une première partie de ce modèle permet de calculer la recharge des nappes en modélisant les interactions précipitations-sol et l’hydrodynamique dans le sol. Cette modélisation a été effectuée en utilisant des modèles conceptuels simples basés sur des lois empiriques (Gardénia, Nash) et des modèles physiques résolvant l’équation de Richards. La recharge ainsi calculée est intégrée à la deuxième partie du modèle simulant l’hydrodynamique dans les nappes décrite par l’équation de diffusivité. Des méthodes numériques précises et robustes ont été utilisées pour résoudre les équations du modèle mathématique : les éléments finis non conformes ont été utilisés pour résoudre l’équation de diffusivité et l’équation de Richards est résolue sous sa forme mixte par une méthode itérative en temps. En somme, ce modèle couplé permet de décrire les variations de niveaux de nappe à partir des données météorologiques connaissant les paramètres caractéristiques de cet aquifère. [...] / Groundwater is the main available water resource for many countries; they are mainly replenished by water from precipitation, called groundwater recharge. Due to its great importance, management of groundwater resources is more essential than ever, and is achieved through mathematical models which offer us a better understanding of physical phenomena as well as their prediction. Hydrogeological systems are generally complex thus characterized by a highly variable dynamic over time and space. These complexities have attracted the attention of many hydrogeologists and many sophisticated models that can handle these issues and describe these systems accurately were developed. Unfortunately, modeling groundwater recharge is still a challenge in groundwater resource management. Generally, groundwater models are used to simulate aquifers flow without a good estimation of recharge and its spatial-temporal distribution. As groundwater recharge rates show spatial-temporal variability due to climatic conditions, land use, and hydrogeological heterogeneity, these methods have limitations in dealing with these characteristics. To overcome these limitations, a coupled model which simulates flow in the unsaturated zone and recharge as well as groundwater flow was developed. The flow in the unsaturated zone is solved either with resolution of Richards equation or with empirical models while the diffusivity equation governs flow in the saturated zone. Robust numerical methods were used to solve these equations: we apply non-conforming finite element to solve the diffusivity equation and we used an accurate and efficient method for solving the Richards equation. [...]

L’aquifère molassique du Bas-Dauphiné est situé le long de la vallée du Rhône, dans le Sud-Est de la France. Cet aquifère d’une superficie proche de3000 km², et d’une épaisseur moyenne de 400m renferme une eau d’excellente qualité, utilisée par de nombreuses collectivités pour l’alimentation en eau potable. Cependant certains secteurs montrent une forte vulnérabilité de la nappe aux activités agricoles. L’utilisation d’outils géochimiques etisotopiques a permis de préciser le fonctionnement hydrodynamique de l’aquifère. Dans un premier temps, nous avons montré que la stratification des écoulements décrite par De La Vaissière (2006) sur la partie drômoise de l’aquifère s’étend au secteur isérois. Les eaux les plus profondes ont des vitesses de circulation de l’ordre du mètre par an alors que les flux superficiels ont des vitesses de circulations d’une centaine de mètres par an. D’autre part, le marquage des nappes superficielles et des rivières par des teneurs faibles en tritium (de 3 à 4 UT) et forte en magnésium (jusqu’à 18 mg/L)indique un apport d’eaux anciennes, issues de l’aquifère molassique vers ces eaux superficielles. La définition de deux pôles d’eaux et l’application d’une équation de mélange couplés à la réalisation de bilans hydrogéologique a permis d’appréhender les volumes échangés. Il apparaît finalement que les réservoirs d’eaux superficiels constituent l’exutoire principal de l’aquifère molassique. L’utilisation des éléments traces a mis en avant le rôleessentiel du temps de séjour des eaux dans l’aquifère ainsi que des conditions d’oxydo-réduction dans l’acquisition de la minéralisation. L’évaluation de la qualité naturelle des eaux de la nappe a montré l’impact des activités agricoles sur l’aquifère, avec prés de 80% des échantillons ayant une teneur en nitrates supérieure à la concentration naturelle supposée. L’étude des teneurs en pesticides conforte ce constat. De plus l’étude de l’évolution des concentrations en polluants montre une dégradation de la ressource. / The molassic basin of Bas-Dauphiné is located in south-eastern France, in the Rhône valley. With an averagethickness of about 400 m, and a surface area of about 2900 km², this aquifer is an important groundwaterresource for freshwater supply and agriculture. However, this resource is also vulnerable and is impacted byhuman activities. The use of geochemical and isotopic analyses made it possible to understand the hydrodynamicsurvey of the aquifer. Firstly, we showed the stratification of the groundwater previously described in thesouthern part of the molassic aquifer could be extended to the north. The groundwater flow velocity is about onemeter / year for the deepest flow and around 100-200 meter /year for the shallowest flow. Secondly, superficialaquifers and rivers are marked by low tritium activities (3 to 4 UT) and high magnesium concentrations (until 18mg/L), which indicate ancient molassic water contribution. The definition of two water types and the applicationof a mixing equation combined with hydrogeological balances lead to an estimation of the contribution of thedeep aquifer to the shallow aquifer. It is finally established that surface water (aquifers and rivers) constitute themain outlet of the molassic aquifer. The use of trace elements shows the importance of the groundwaterresidence time and of redox conditions in the water mineralization. The assessment of baseline quality shows theimpact of agricultural activities upon the aquifer, with nearly 80 % of samples showing higher nitrateconcentration than the maximal natural concentration. This is confirmed by the study of pesticidesconcentrations. Furthermore, the evolution of pollutants concentrations points out a deterioration of the resource.