Χρήση Γεωγραφικών Συστημάτων Πληροφοριών (G.I.S.) για επίλυση προβλημάτων περιβαλλοντικής υδρογεωλογίας στο νομό Κορινθίας

Abstract

The present doctoral thesis examines the potential applications of Geographic Information Systems (GIS) in the scientific field of environmental hydrogeology. The research area for the application of these methods is the northeastern part of the prefecture of Corinth which includes from east to the west the Xerias, Rachianis, Zapantis and Asopos river basins. The wider study area includes both prealpine and postalpine geological formations of terrestrial, lagoonal and marine origin. The dominant geomorphological and tectonic structure in the wider study area is the Corinthian Gulf which is a wide rift zone with high rates of displacement between the two parts of the Peloponnese and Central Greece, as a result of the tensile tectonics that has affected and still affects the region after the end of the Alpine Orogenic episode in the middle Miocene. The neotectonic of the study area is characterized by intense tensile stresses responsible for the creation of the Corinthian trench and numerous normal and oblique faults, parallel and perpendicular to the axis of the trench. A large proportion of these faults are active and responsible for the genesis of medium and large magnitude earthquakes during historical times. The relief of the survey area is characterized by three large plane sections, the coastal area between Corinth and Kiato known as Vochas plane, the area between Agios Basilios and Spathovouni villages and the plane area around the historical city of Nemea. The remainder of the study area can be described as hilly or semi - mountainous with moderate slopes and rounded mountain peaks, with the exception of the SW and SE edge of the area where the predominance of carbonate formations of the prealpine basement creates mountainous terrain with steep slopes Slope deep ravines and sharp mountain peaks. The morphometric analysis shows that watersheds in the study area are expected to have hydrographs of “slow runoff” and low flooding potential. The morphometric characteristics of Zapantis, Ellison and Seliandros basins make them more susceptible to “flash flood” events. From the hydrogeological point of view, a number of ground water bodies (aquifers) and practically impermeable formations are developing in the study area. To investigate the movement of groundwater in aquifers, six different spatial interpolation methods were implemented and it was found that the interpolation method, considerably affects the final formation of the piezometric surface. The criteria for the evaluation of each method were firstly whether the resulting surface does not alter the physical properties a piezometric surface must have (small hydraulic gradient, agreement with the topography, etc.) and secondly a number of statistical error indicators quantify the agreement of the resulting surface with the existing values at the measuring points. Combining the above mentioned observations, we conclude that the most appropriate method for creating groundwater level surfaces, with the process of spatial interpolation, for cases of regional areas with alternating aquifer and practically non- aquifer formations and general heterogeneity in hydraulic characteristics, is the method of Ordinary Cokriging with prior removal of data trends. The annual variation of ground water level, as it was shown by the comparison of ground water level surface of periods 5/2003 and 10/2004 in alluvial aquifers, except the coastal portion of Vocha aquifer, we observe an increase in ground water levels varying from 1 to 10 meters while within the conglomerate aquifers and the confined aquifer that develops within the marl formation we observe a slight decline ranging from 1 to 10 meters. This variation can be attributed to the fact that shallow alluvial aquifers are recharged faster than the deeper confined aquifers and are influenced by the water balance of only the current hydrological year while the deeper confined aquifers are affected by more than one hydrologic year. The corresponding situation for the periods 10/2004 and 10/2005 is almost inverse because a rising in ground water level is observed both in the deeper confined aquifers and in the coastal part of the alluvial aquifer of Vocha while drawdown is observed in other alluvial aquifers. The repeated differentiation between aluvial and deeper confined aquifers suggests that they are responding to recharge with a different time lag. Based mainly on hydrolithological classification of geological formations but also on their spatial distribution as it is represented in the hydrogeological map but additionally on the basis of information from borehole lithological sections and also the spatial distribution of the groundwater level we end in the delineation of “groundwater bodies” or “aquifer systems” for the study area. According to the above mentioned evidence five porous and two karstic aquifer systems were defined. Both karstic aquifer systems are developed partly in the study area with their bigger part located outside of the study area. For this reason, they were not investigated thoroughly in this research. The feature that was examined in these systems was the way they interact with other aquifer systems and how they affect the hydrological balance of the watersheds of the study area. Three of the five porous aquifer systems (Vocha, Agios Basilios and Nemea) developed entirely within the study area and they include mostly shallow alluvial aquifers which form unified or multilevel aquifers with deeper confined aquifers of the Plio - Pleistocene formations. The other two porous aquifer systems (Aidonia - Petri and Kryoneri - Baltos) develope within the Pleistocene conglomerates (Calabrian conglomerates) and despite the fact that they are partially included in the study area they were investigated in detail because they play a key role in water balance of the study area within which they are discharged through high capacity contact springs in places where with come in contact with the underlying impermeable Pliocene marls. The most prolific aquifers in the study area are the conglomerate aquifers while aquifers of Vocha, Nemea and Agios Basilios - Spathovouni have a highly variable productivity and have average to low yields. Finally it was found that the aquifer developed in the marl formation has low productivity and small practical importance. The extensive development of cultivations in the study area bring significant pressures on water resources in both the quantitative part due to increased irrigation needs and the qualitative part because of diffuse pollution on water resources resulting from agricultural practices (fertilizers, pesticides etc.). The various aquifers in the study area have significant hydrochemical variations which are found both in the distribution they have in Piper diagram and from the various relationships between hydrochemical parameters such as the relationship between alkalinity and calcium or alkalinity and concentrations of REE and Eh-pH, but also in the statistical indicators of hydrochemical parameters such as alkalinity, conductivity, sulphate, nitrate and heavy metals concentrations. These variations are due to the different composition of the aquifer materials, the difference in residence time of groundwater within the aquifer, the redox conditions and other factors. These variations are normalized in cases where there is hydraulic connection between different aquifers. As it was shown by the spatial distribution of pH values, acidic conditions prevailing almost in the entire coastal alluvial aquifer system of Vocha with the most acidic conditions located in the areas of Velo and Assos - Lecheo while alkaline conditions prevailing in the water systems of conglomerates and also in the deeper confined aquifer that develops in the Isthmia - Examillia - Ancient Corinth area. Correspondingly the spatial distribution of redox potential reveal that oxidizing conditions prevail in most of the study area except for the part Vrahati -Kokoni - Evangelistria in Vocha alluvial aquifer, the Pasio - Laliotis – Sikyon area and large parts of the aquifer developed in the formation of marls. By studying the spatial distribution of nitrogen compounds in groundwater, we come to the conclusion that there is a widespread contamination of groundwater by nitrate, mainly in the Vocha coastal aquifer system except the area between Vrahati - Kokoni - Neratza and also throughout the whole area between Examilia - Solomos - Spathovouni - Agios Basilios both in the upper unconfined and deeper confined aquifers. As shown by the spatial distribution of the parameters associated with seawater intrusion, there seems to be a propagation of the seawater intrusion front, to the axis between Assos and Spathovouni even at great distance from the coast. This conclusion is confirmed by the fact that, deep boreholes in the Spathovouni area have been abandoned due to salinization. Similarly there seems to be a retreat of the seawater intrusion front near Isthmia-Examillia area. From the simultaneous observation and study of the spatial distribution of chloride, ground water level and the hydrochemical types of groundwater, it is obvious that in the study area, sea water intrusion represent a very serious problem mainly along the axis indicated above. In the rest of the coastal zone, although hydrogeological conditions do not differs seawater intrusion is quite limited. This fact may be attributed both on good recharge of coastal aquifers and on the relative uplift of the practically impermeable marl formation which acts as a natural barrier to seawater intrusion. The alluvial aquifers of the study area are showing the greatest concentration of heavy metals. In some cases the concentration values are so high that they can be attributed to contamination of groundwater from point sources or to the influence by hydrothermal fluids. For the better understanding of hydrogeological and hydrochemical processes in aquifers in the study area, the spatial distributions of rare earth elements (REE), as well as their evolution along groundwater flow lines were investigated in response to evolution of other hydrochemical parameters such as concentrations of main elements and redox conditions. The aquifers developed in marls and marine terraces show the higher values of REE elements, which is an indication that in fact these aquifers are in full hydraulic connection or even compose a single confined aquifer. The conglomerates and alluvial aquifers of Vocha and Nemea show generally lower REE concentration values while the aquifer system of Agios Basilios-Spathovouni generally have high REE concentration values which also indicates that the deeper confined aquifer of the system, is hydraulically connected or developed within the Plio - Pleistocene marls. The cumulative concentration of rare earth elements seems to increase with the decrease in redox potential as well as the decrease in phi Acidic and reducing conditions appear to favor the preservation of dissolved rare earth elements in groundwater. According to the spatial distribution of δ18Ο values and according to the relationship between δ18Ο values and the altitude of the sampling point, the aquifer systems of the conglomerates appears to be recharged directly from the area where they have surface development, or possibly laterally recharged with groundwater of the adjacent karst aquifers initially recharged in areas of higher altitude. The remainder of the study area appears to be recharged from the areas where the formations that host the aquifer have surface development, excluding the areas of Sikyon-Velo and Lecheo that seems to be recharged with water from areas with higher altitudes. High tritium values were observed in parts of the coastal alluvial aquifer of Vocha, but also in the porous aquifer system of Nemea, particularly around Koutsi while no big discrepancies between δ18Ο and altitude of sampling points occurred in the same sites. High rates of tritium can be attributed mainly to problematic recharge and discharge of aquifers, resulting in long residence times of groundwater in the aquifer. Under the framework of application of advanced methods using GIS, simulation models and spatial multi-criteria analysis models were used and developed, for solving problems related to hydrological and hydrochemical processes at river basin scale, the vulnerability of aquifers and the placement of new water wells. The river basin simulation model that was applied is called “SWAT”, acronym for Soil and Water Assessment Tool, which is a fully distributed physically based model that may be applied in ArcGIS environment. Application of SWAT model in the study area was conducted for the four main river basins (Asopos, Zapantis, Rachiani and Xerias) for the period 1991-2001 with annual and monthly time step. After the initial implementation, the model was manually calibrated using time series of runoff measurements from a station on the River Asopos, namely the bridge between Nemea and Petrio. The results of the model include distributions of parameters such as the direct and total runoff, recharge of the phreatic aquifer and the quantity of nitrogen and nitrate that reaches the surface water and water of aquifer recharge, on a sub basin or basin scale. The investigation of aquifer vulnerability is carried out by the optimization of the standard DRASTIC method and various modifications and transformations are proposed, based on statistical parameters of the distribution of nitrate concentration as an indicator of aquifer pollution. Eventually, an integrated method, that can successfully determine the specific vulnerability or the pollution risk of an aquifer under intense environmental pressure, is provided. Besides the above optimization process of the DRASTIC model, two other hybrid methodologies based on the principles of DRASTIC method were developed. These two methodologies are based on two methods of multivariate spatial statistics and more specifically the method of " multiple logistic regression » and the method of “weighted evidence”, in order to determine the specific vulnerability (pollution risk) of aquifers. The results of the three methods are compared in terms of success on the coincidence of high values of pollution risk with high nitrate concentrations in groundwater. By creating a multi-criteria analysis model using fuzzy logic methods and hydrogeological, environmental and socioeconomic criteria, the distribution of the suitability of each location of study the area for drilling new water wells was accomplished. From the distribution of suitability it appears that large parts of the study area are totally unsuitable for drilling new water wells. This is due to the large number of existing water wells and the large extend of the marl formation which holds very low productivity aquifers. The choice of the method of fuzzy logic and the combination of criteria used, manage to focus in fully appropriate areas for the drilling of new water wells. These areas are then further classified based on the combination of the above mentioned criteria. The methodology of multi-criteria analysis using fuzzy logic has great flexibility and adaptability. By altering the type and coefficients of fuzzy functions and choosing a different combination of operators the final result can be tailored to the priorities of research. The best positions in the study area are located in the formation of conglomerates, west of Nemea plane area and around villages Kryoneri and Souli as well as the upstream part of the alluvial aquifer of Vocha and within recent river deposits nearby all the rivers of the study area

Η παρούσα διδακτορική διατριβή εξετάζει τις πιθανές εφαρμογές των Γεωγραφικών Συστημάτων Πληροφοριών (ΓΣΠ) στο επιστημονικό πεδίο της περιβαλλοντικής υδρογεωλογίας. Ως περιοχή έρευνας για την εφαρμογή των μεθόδων αυτών οριοθετείται το ΒΑ τμήμα του νομού Κορινθίας που περιλαμβάνει της λεκάνες απορροής Ξεριά, Ράχιανης, Ζαπάντη και Ασωπού. Η γεωλογική δομή της ευρύτερης περιοχής έρευνας περιλαμβάνει σχηματισμούς τόσο του Προνεογενούς υποβάθρου όσο και μεταλπικούς χερσαίους, λιμνοθαλάσσιους και θαλάσσιους σχηματισμούς. Το κυρίαρχο από γεωμορφολογικής και γεωτεκτονικής άποψης στοιχείο στην ευρύτερη περιοχή έρευνας αποτελεί ο Κορινθιακός κόλπος ο οποίος είναι ένα ευρύ τεκτονικό βύθισμα με υψηλούς ρυθμούς απομάκρυνσης, μεταξύ των δυο τεμαχών, της Πελοποννήσου και της στερεάς Ελλάδας, αποτέλεσμα της εφελκυστικής τεκτονικής που επέδρασε και συνεχίζει να επιδρά στην ευρύτερη περιοχή μετά το τέλος των αλπικών πτυχώσεων κατά το Μέσο Μειόκαινο. Το προνεογενές υπόβαθρο της περιοχής έρευνας δομείτε από τις γεωτεκτονικές ζώνες της Υποπελαγονικής, της ζώνης Γαβρόβου-Τριπόλεως και της ζώνης Πίνδου οι οποίες εμφανίζονται σε μικρά τμήματα στο ΝΑ και ΝΔ άκρο της περιοχής έρευνας. Τα ιζήματα του Τεταρτογενούς αποτελούνται από ποταμοχερσαία αλουβιακά ριπίδια και παράκτια λιμνοθαλάσσια κροκαλοπαγή στα νότια, λιμνοθαλάσσια δελταϊκά κροκαλοπαγή τύπου Gilbert και λιμνοθαλάσσιες μάργες στο κέντρο και τέλος λιμνοθαλάσσιες μάργες και θαλάσσιες αναβαθμίδες από μάργες, ψαμμίτες και κροκαλοπαγή στο βορειότερο παράκτιο τμήμα. Η νεοτεκτονική της περιοχής έρευνας χαρακτηρίζεται από τις έντονες εφελκυστικές τάσεις που είναι υπεύθυνες για την δημιουργία της Κορινθιακής τάφρου και των πολυάριθμων κανονικών και πλαγιοκανονικών ρηγμάτων παράλληλα και κάθετα προς τον άξονα της τάφρου. Μεγάλο ποσοστό των ρηγμάτων αυτών είναι ενεργά και υπεύθυνα για την γένεση σεισμών μεσαίου και μεγάλου μεγέθους κατά τους ιστορικούς χρόνους. Το ανάγλυφο της περιοχής έρευνας χαρακτηρίζεται από τρεις μεγάλες πεδινές ενότητες, την παράκτια περιοχή μεταξύ Κορίνθου και Κιάτου γνωστή και ως κάμπος της Βόχας, την περιοχή Αγίου Βασιλείου – Σπαθοβουνίου και την πεδινή περιοχή Νεμέας. Το υπόλοιπο της περιοχής έρευνας μπορεί να χαρακτηρισθεί ως λοφώδες-ημιορεινό με μέτριες κλίσεις πρανών και ομαλές απολήξεις κορυφογραμμών με εξαίρεση το ΝΔ και ΝΑ άκρο της περιοχής όπου η επικράτηση των ανθρακικών σχηματισμών του προνεογενούς υποβάθρου δημιουργεί ορεινό ανάγλυφο με μεγάλες κλίσεις πρανών βαθιές χαραδρώσεις και οξύληκτες κορυφογραμμές. Η μορφομετρική ανάλυση καταδεικνύει ότι οι λεκάνες της περιοχής έρευνας αναμένεται να παρουσιάζουν υδρογραφήματα «βραδείας απορροής» και χαμηλό δυναμικό εκδήλωσης πλημμυρών. Οι λεκάνες των Ζαπάντη, Ελλισώνα και Σέλιανδρου είναι αυτές που ο συνδυασμός μορφομετρικών χαρακτηριστικών καθιστούν πιο επιδεκτικές στην εκδήλωση «ταχείας απορροής» και «αστραπιαίων» πλημμυρικών επεισοδίων. Από υδρογεωλογικής άποψης στην περιοχή αναπτύσσονται μια σειρά από υδροφόρα συστήματα καθώς και πρακτικά στεγανοί σχηματισμοί. Για την διερεύνηση της κίνησης του υπόγειου νερού στα υδροφόρα συστήματα αυτά εφαρμόστηκαν έξι διαφορετικές μέθοδοι χωρικής παρεμβολής και διαπιστώθηκε πώς η επιλογή μεθόδου, επηρεάζει καθοριστικά το τελικό αποτέλεσμα. Τα κριτήρια για την αξιολόγηση κάθε μεθόδου ήταν αφενός το κατά πόσο η επιφάνεια που προκύπτει δεν αλλοιώνει τις φυσικές ιδιότητες που πρέπει να έχει μια πιεζομετρική επιφάνεια (μικρές υδραυλικές κλίσεις, συμφωνία με την τοπογραφία κ.τ.λ.) και αφετέρου μια σειρά από στατιστικούς δείκτες σφάλματος που εκφράζουν ποσοτικά την συμφωνία της παραγόμενης επιφάνειας με τις υφιστάμενες τιμές στα σημεία μέτρησης. Συνδυάζοντας τις ανωτέρω παρατηρήσεις καταλήγουμε στο συμπέρασμα ότι, η πλέον ενδεικνυόμενη μέθοδος για την δημιουργία κατανομής στάθμης με την διαδικασία της χωρικής παρεμβολής, για περιπτώσεις ευρύτερων περιοχών με εναλλαγές υδροφόρων και πρακτικά μη υδροφόρων σχηματισμών και γενικά ανομοιογένεια στα υδραυλικά χαρακτηριστικά είναι η μέθοδος Ordinary Cokriging με απομάκρυνση των τάσεων (trend remove). Σχετικά με την υπερετήσια μεταβολή της στάθμης του υπόγειου νερού, όπως προκύπτει από την σύγκριση των επιφανειών στάθμης μεταξύ των περιόδων 5/2003 και 10/2004 στους προσχωματικούς υδροφόρους, εκτός του παράκτιου τμήματος του υδροφόρου της Βόχας, παρατηρούμε μια άνοδο του επιπέδου στάθμης που κυμαίνεται από 1 έως και 10 μέτρα ενώ στους υδροφόρους των κροκαλοπαγών αλλά και του υπό πίεση ορίζοντα που αναπτύσσεται εντός του σχηματισμού των μαργών παρατηρούμε μια μικρή πτώση που κυμαίνεται από 1 έως και 10 μέτρα. Η διαφοροποίηση αυτή μπορεί να αποδοθεί στο ότι οι αβαθείς προσχωματικοί υδροφόροι έχουν ταχύτερη επαναπλήρωση σε σχέση με τους βαθύτερους ελεύθερους ή υπό πίεση ορίζοντες και επηρεάζονται από το υδρολογικό ισοζύγιο μόνο του τρέχοντος υδρολογικού έτους ενώ οι βαθύτεροι υδροφόροι ορίζοντες επηρεάζονται από παραπάνω από ένα υδρολογικά έτη. Η αντίστοιχη εικόνα για τις περιόδους 10/2004 και 10/2005 είναι σχεδόν αντεστραμμένη με την άνοδο της στάθμης να παρατηρείται τόσο στους βαθύτερους υδροφόρους όσο και στο παράκτιο τμήμα του προσχωματικού υδροφόρου της Βόχας, ενώ πτώση στάθμης παρατηρείται στους υπόλοιπους προσχωματικούς υδροφόρους. Η επανάληψη της διαφοροποίησης μεταξύ προσχωματικών και βαθύτερων υδροφόρων ενισχύει την άποψη ότι επηρεάζονται με διαφορετική χρονική υστέρηση από το υδρολογικό ισοζύγιο. Με βάση κυρίως την υδρολιθολογική κατάταξη των σχηματισμών αλλά και την χωρική τους κατανομή όπως αποτυπώνεται στον υδρογεωλογικό χάρτη αλλά επιπρόσθετα με βάση πληροφορίες από τις τομές των γεωτρήσεων που απογράφτηκαν αλλά και από την κατανομή της στάθμης του υπόγειου νερού καταλήγουμε στην οριοθέτηση των υδρογεωλογικών ενοτήτων για την περιοχή έρευνας. Σύμφωνα με τα προαναφερθέντα στοιχεία οριοθετούνται πέντε κοκκώδη υδροφόρα συστήματα, και δύο καρστικά υδροφόρα συστήματα στα νότια όρια της περιοχής έρευνας. Και τα δύο καρστικά υδροφόρα συστήματα αναπτύσσονται τμηματικά στην περιοχή έρευνας με το μεγαλύτερο τμήμα τους να βρίσκεται εκτός περιοχής. Για το λόγο αυτό δεν ερευνήθηκαν διεξοδικά στο πλαίσιο της παρούσας διατριβής. Το χαρακτηριστικό που εξετάστηκε σε αυτά τα συστήματα ήταν ο τρόπος με τον οποίο αλληλεπιδρούν με τα υπόλοιπα υδροφόρα συστήματα αλλά και ο τρόπος που επηρεάζουν τα υδρολογικά ισοζύγια των λεκανών απορροής της περιοχής έρευνας. Από τα πέντε κοκκώδη υδροφόρα συστήματα τα τρία (Βόχας, Αγίου Βασιλείου και Νεμέας) αναπτύσσονται εξ ολοκλήρου στην περιοχή έρευνας και περιλαμβάνουν κυρίως αβαθείς προσχωματικούς υδροφόρους οι οποίοι αναπτύσσουν ενιαία ή και επάλληλη υδροφορία με βαθύτερους υδροφόρους ορίζοντες των Πλειο-Πλειστοκαινικών ιζημάτων. Τα υπόλοιπα δύο κοκκώδη συστήματα (Αηδονίων-Πετρίου και Κρυονερίου-Βάλτου) αναπτύσσονται εντός των κροκαλοπαγών του Πλειστοκαίνου (Καλάβρια κροκαλοπαγή) και παρά το γεγονός ότι περιλαμβάνονται τμηματικά στην περιοχή έρευνας ερευνήθηκαν διεξοδικά διότι παίζουν καθοριστικό ρόλο στην υδροοικονομία της περιοχής έρευνας εντός της οποίας εκφορτίζονται μέσω υψηλής δυναμικότητας πηγών επαφής με τις υποκείμενες μάργες του Πλειοκαίνου. Οι αποδοτικότεροι υδροφόροι στην περιοχή έρευνας είναι εκείνοι των κροκαλοπαγών ενώ οι υδροφόροι που αναπτύσσονται στα υδροφόρα συστήματα Βόχας, Νεμέας και Αγίου Βασιλείου-Σπαθοβουνίου παρουσιάζουν κυμαινόμενη απόδοση με γενικά μέσες έως χαμηλές αποδόσεις. Τέλος διαπιστώνεται ότι η υδροφορία που αναπτύσσεται εντός του σχηματισμού των μαργών είναι χαμηλής απόδοσης χωρίς ιδιαίτερη πρακτική σημασία. Η ανάπτυξη των καλλιεργειών στην περιοχή έρευνας επιφέρει σημαντικές πιέσεις στους υδατικούς πόρους τόσο στο ποσοτικό κομμάτι λόγω των αυξημένων αρδευτικών αναγκών όσο και στο ποιοτικό κομμάτι λόγω της διάχυτης μόλυνσης που καταλήγει στους υδατικούς πόρους από τις καλλιεργητικές πρακτικές (λιπάσματα, φυτοφάρμακα). Οι διάφοροι υδροφόροι της περιοχής έρευνας παρουσιάζουν σημαντικές υδροχημικές διαφοροποιήσεις οι οποίες εντοπίζονται τόσο από την κατανομή τους στο διάγραμμα Piper όσο και από τους διάφορους συσχετισμούς μεταξύ των υδροχημικών παραμέτρων όπως για παράδειγμα η σχέση αλκαλικότητας – ασβεστίου, αλκαλικότητας-συγκέντρωσης σπάνιων γαιών και Eh–pH, αλλά και στους στατιστικούς δείκτες υδροχημικών παραμέτρων όπως η αλκαλικότητα, η αγωγιμότητα, τα θειικά και νιτρικά ιόντα και η συγκέντρωση βαρέων μετάλλων. Οι διαφοροποιήσεις αυτές οφείλονται στην διαφορετική σύσταση του σκελετού των υδροφόρων, στην διαφορά του χρόνου παραμονής του υπόγειου νερού εντός των υδροφόρων, των οξειδοαναγωγικών συνθηκών και άλλων παραγόντων οι οποίες βέβαια ομαλοποιούνται στις περιπτώσεις που υφίσταται υδραυλική επικοινωνία μεταξύ των υδροφόρων. Όπως προκύπτει από την μελέτη της κατανομής του pH, όξινες συνθήκες επικρατούν σχεδόν στο σύνολο του παράκτιου προσχωματικού υδροφόρου συστήματος της Βόχας με επίκεντρο τις περιοχές Βέλου και Άσσου-Λεχαίου, ενώ αλκαλικές συνθήκες επικρατούν στα υδροφόρα συστήματα των κροκαλοπαγών αλλά και στην περιοχή Ισθμίων-Εξαμιλλίων-Αρχαίας Κορίνθου στον βαθύτερο υδροφόρο ορίζοντα που αναπτύσσεται στην περιοχή. Αντίστοιχα από την κατανομή του δυναμικού οξειδοαναγωγής προκύπτει ότι οξειδωτικές συνθήκες επικρατούν στο μεγαλύτερο τμήμα της περιοχής έρευνας με εξαίρεση το τμήμα Βραχατίου-Κοκκωνίου-Ευαγγελίστριας στον προσχωματικό υδροφόρο της Βόχας το τμήμα Πάσιου-Λαλιώτη-Σικυώνος καθώς και μεγάλα τμήματα των υδροφόρων που αναπτύσσονται στον σχηματισμό των μαργών. Από την μελέτη των κατανομών των ενώσεων αζώτου στο υπόγειο νερό, διαπιστώνεται μια ευρεία ρύπανση του υπόγειου νερού από νιτρικά ιόντα κυρίως στο παράκτιο υδροφόρο σύστημα της Βόχας με εξαίρεση την περιοχή μεταξύ Βραχατίου-Κοκκωνίου-Νεράτζας αλλά και σε ολόκληρη την ζώνη μεταξύ Εξαμιλίων-Σολωμού-Σπαθοβουνίου-Αγίου Βασιλείου τόσο στους ελεύθερους όσο και στους βαθύτερους υπό πίεση υδροφόρους ορίζοντες. Όπως φαίνεται από την κατανομή των παραμέτρων που σχετίζονται με υφαλμύρινση, επέκταση του μετώπου υφαλμύρινσης φαίνεται να έχουμε στον άξονα μεταξύ Άσσου και Σπαθοβουνίου και μάλιστα σε μεγάλη απόσταση από την ακτογραμμή. Η παρατήρηση αυτή ενισχύεται από το γεγονός ότι γεωτρήσεις μεγάλου βάθους στην περιοχή του Σπαθοβουνίου έχουν εγκαταλειφθεί λόγω υφαλμύρινσης. Αντίστοιχα υποχώρηση του μετώπου υφαλμύρινσης, φαίνεται να έχουμε στην περιοχή Ισθμίων-Εξαμιλλίων. Όπως λοιπόν διαπιστώνουμε, από την ταυτόχρονη παρατήρηση και μελέτη της κατανομής των χλωριόντων, της πιεζομετρίας και των υδροχημικών τύπων του υπόγειου νερού, στην περιοχή έρευνας εντοπίζονται έντονα προβλήματα υφαλμύρινσης, τα οποία αναπτύσσονται κυρίως κατά μήκος άξονα που προαναφέρθηκε. Στην υπόλοιπη παράκτια ζώνη αν και οι συνθήκες άντλησης δεν διαφοροποιούνται έχουμε περιορισμένες ενδείξεις υφαλμύρινσης γεγονός που μπορεί να αποδοθεί τόσο στην καλή τροφοδοσία και επαναπλήρωση των υδροφόρων όσο και στην σχετική ανύψωση του λεπτομερούς σχηματισμού των μαργών ο οποίος λειτουργεί ως φυσικό φράγμα προς την διείσδυση του θαλασσινού νερού. Οι προσχωματικοί υδροφόροι της περιοχής έρευνας παρουσιάζουν την μεγαλύτερη συγκέντρωση βαρέων μετάλλων ενώ σε ορισμένες περιπτώσεις οι τιμές είναι τόσο υψηλές που παραπέμπουν σε μόλυνση του υδροφόρου από σημειακή πηγές ή σε επίδραση από υδροθερμικά ρευστά. Στο πλαίσιο διερεύνησης των υδρογεωλογικών και υδροχημικών διεργασιών στους υδροφόρους της περιοχής έρευνας, μελετήθηκαν τόσο οι κατανομές των σπάνιων γαιών, όσο και οι μεταβολές των συγκεντρώσεών τους κατά μήκος γραμμών ροής του υπόγειου νερού, σε συνάρτηση με τις μεταβολές των υπόλοιπων υδροχημικών παραμέτρων, όπως οι συγκεντρώσεις των κύριων στοιχείων και οι οξειδοαναγωγικές συνθήκες. Υψηλότερες τιμές σπάνιων γαιών παρουσιάζουν οι υδροφόροι που αναπτύσσονται εντός του σχηματισμού των μαργών και των θαλάσσιων αναβαθμίδων, γεγονός που αποτελεί ένδειξη ότι ουσιαστικά πρόκειται για υδροφόρους σε πλήρη υδραυλική επικοινωνία ή ακόμα και για ένα ενιαίο υδροφόρο τουλάχιστον στην βαθύτερη υπό πίεση υδροφορία. Οι υδροφόροι των κροκαλοπαγών και οι προσχωματικοί υδροφόροι Βόχας και Νεμέας παρουσιάζουν γενικά χαμηλότερες τιμές ενώ στο υδροφόρο σύστημα Αγίου Βασιλείου-Σπαθοβουνίου έχουμε υψηλές γενικά τιμές γεγονός που ομοίως αποτελεί ένδειξη ότι η βαθύτερη υδροφορία του συστήματος αυτού είναι σε υδραυλική επικοινωνία ή αναπτύσσεται εντός των υδροφόρων στρωμάτων που εντοπίζονται στις μάργες του Πλειο-Πλειστοκαίνου. Η αθροιστική συγκέντρωση των σπάνιων γαιών φαίνεται να αυξάνεται με την μείωση του δυναμικού οξειδοαναγωγής όπως επίσης και με την μείωση του pH. Όξινες και αναγωγικές συνθήκες φαίνεται ότι ευνοούν την εν διαλύσει παραμονή των σπάνιων γαιών στο υπόγειο νερό Βάση της κατανομής δ18Ο και σύμφωνα με την σχέση της τιμής του δ18Ο με το υψόμετρο της περιοχής προέλευσης, τα υδροφόρα συστήματα των κροκαλοπαγών φαίνεται να τροφοδοτούνται απευθείας από την περιοχή στην οποία αναπτύσσονται επιφανειακά ή ενδεχομένως και πλευρικά με υπόγειο νερό των καρστικών υδροφόρων που προέρχεται από περιοχές με υψηλότερο υψόμετρο. Το υπόλοιπο της περιοχής έρευνας φαίνεται να τροφοδοτείτε αποκλειστικά από τις περιοχές στις οποίες οι σχηματισμοί που φιλοξενούν την υδροφορία αναπτύσσονται επιφανειακά, με εξαίρεση τις περιοχές Σικυώνος-Βέλου και Λεχαίου που φαίνεται να τροφοδοτούνται με νερά από περιοχές με υψηλότερα υψόμετρα. Υψηλές τιμές τριτίου παρατηρούνται σε τμήματα του παράκτιου προσχωματικού υδροφόρου της Βόχας αλλά και στο κοκκώδες υδροφόρο σύστημα Νεμέας και κυρίως στην περιοχή Κουτσίου χωρίς να παρατηρούνται μεγάλες αποκλίσεις μεταξύ δ18Ο και υψομέτρου δειγματοληψίας στις ίδιες θέσεις. Οι υψηλές τιμές τριτίου μπορούν να αποδοθούν κυρίως στην προβληματική τροφοδοσία και εκφόρτιση του υδροφόρου που οδηγεί σε μεγάλους χρόνους παραμονής του υπόγειου νερού στον υδροφόρο. Στο πλαίσιο εφαρμογής προηγμένων μεθόδων με την χρήση ΓΠΣ εφαρμόσθηκαν και δημιουργήθηκαν μοντέλα προσομοίωσης και μοντέλα πολυκριτηριακής χωρικής ανάλυσης για την επίλυση προβλημάτων που αφορούν τις υδρολογικές και υδροχημικές διεργασίες σε κλίμακα λεκάνης απορροής, την τρωτότητα των υδροφόρων οριζόντων και την χωροθέτηση νέων υδρομαστευτικών έργων. Το μοντέλο προσομοίωσης λεκάνης απορροής που εφαρμόσθηκε είναι το SWAT, ακρωνύμιο των λέξεων Soil and Water Assessment Tool, το οποίο είναι ένα πλήρως κατανεμημένο (fully distributed), φυσικής βάσεως (physically based) μοντέλο που εκτελείται σε περιβάλλον ArcGIS. Η εφαρμογή του μοντέλου SWAT στην περιοχή έρευνας πραγματοποιήθηκε για τις τέσσερις κύριες υδρολογικές λεκάνες (Ασωπός, Ζαπάντης, Ράχιανη, Ξεριάς) και για την περίοδο 1991-2001 με ετήσιο και μηνιαίο χρονικό βήμα. Μετά την αρχική εφαρμογή του μοντέλου ακολούθησε χειροκίνητη ρύθμιση του μοντέλου με την χρήση χρονοσειράς μετρήσεων απορροής από σταθμό μέτρησης στον ποταμό Ασωπό και συγκεκριμένα στη γέφυρα Νεμέας-Πετρίου. Στα αποτελέσματα του μοντέλου περιλαμβάνονται οι κατανομές ανά υδρολογική υπολεκάνη και ανά υδρολογική λεκάνη παραμέτρων όπως η άμεση και η ολική επιφανειακή απορροή, η τροφοδοσία του φρεάτιου υδροφόρου ορίζοντα και η ποσότητα αζώτου και νιτρικών που φτάνει στο επιφανειακό νερό και στο νερό τροφοδοσίας του υδροφόρου. Η μελέτη της τρωτότητας των υδροφόρων πραγματοποιεί¬ται με την βελτιστοποίηση της πρότυπης μεθόδου DRASTIC και προτείνονται διάφορες τροποποιήσεις και μετασχηματισ¬μοί, που βασίζονται στις στατιστικές παραμέτρους της κατανομής της συγκέντρωσης νιτρικών ως δείκτη μόλυνσης. Τελικά επιτυγχάνεται η ανάπτυξη μιας ολοκληρωμένης μεθόδου, η οποία θα μπορεί να προβλέπει επιτυχώς την ειδική (specific) τρωτότητα ή το δυναμικό μόλυνσης (Pollution Risk) ενός υδροφόρου ορίζοντα, ο οποίος βρίσκεται κάτω από έντονη περιβαλλοντική πίεση. Εκτός της παραπάνω διαδικασίας βελτιστοποίησης αναπτύσσονται και δύο άλλες μεθοδολογίες οι οποίες βασίζονται στις αρχές της μεθόδου DRASTIC αλλά χρησιμοποιούν δύο μεθόδους χωρικής στατιστικής και πιο συγκεκριμένα την μέθοδο της «πολλαπλής λογιστικής παλινδρόμησης» (multiple logistic regression) και την μέθοδο «βαρύτητας των στοιχείων» (weighted evidence), για τον προσδιορισμό της ειδικής τρωτότητας των υδροφόρων. Τα αποτελέσματα των τριών μεθόδων συγκρίνονται με κριτήριο επιτυχίας την σύμπτωση των υψηλών τιμών δυναμικού μόλυνσης με τις υψηλές συγκεντρώσεις νιτρικών στο υπόγειο νερό. Με την δημιουργία ενός μοντέλου πολυκριτηριακής ανάλυσης με την χρήση μεθόδων ασαφούς λογικής και κριτήρια υδρογεωλογικά, περιβαλλοντικά και οικονομοτεχνικά προσδιορίζεται η κατανομή της καταλληλότητας κάθε θέσης της περιοχής έρευνας για την ανόρυξη νέων υδρομαστευτικών έργων. Από την κατανομή της καταλληλότητας φαίνεται ότι μεγάλες εκτάσεις της περιοχής έρευνας είναι εντελώς ακατάλληλες για ανόρυξη νέων υδρομαστευτικών έργων, γεγονός που οφείλεται κυρίως στην ύπαρξη πολυάριθμων υδρομαστευτικών έργων που καθιστούν αδύνατη την ανόρυξη νέων έργων στην περιοχή γειτονίας τους αλλά και στην ευρεία εξάπλωση του σχηματισμού των μαργών εντός του οποίου δεν αναπτύσσεται αξιόλογη υδροφορία. Η επιλογή της μεθόδου της ασαφούς λογικής και του συνδυασμού τελεστών που χρησιμοποιήθηκαν αποκλείουν εντελώς τις περιοχές αυτές και καταφέρνουν να επικεντρώσουν τα αποτελέσματα σε θέσεις κατάλληλες από κάθε άποψη, βαθμονομώντας περαιτέρω τις θέσεις αυτές με βάση τη συγκέντρωση των βέλτιστων συνθηκών. Η μεθοδολογία της πολυκριτηριακής ανάλυσης με την χρήση ασαφούς λογικής διαθέτει μεγάλη ευελιξία και προσαρμοστικότητα. Με αλλαγές στο είδος και τους δείκτες των συναρτήσεων συμμετοχής αλλά και επιλογή διαφορετικού συνδυασμού τελεστών μπορεί το τελικό αποτέλεσμα να διαμορφωθεί ανάλογα με τις προτεραιότητες της έρευνας. Εάν για παράδειγμα η έρευνα επικεντρώνεται σε έργα υδρομάστευσης με σκοπό την ύδρευση μπορούσε να αυξήσουμε την επίδραση των κριτηρίων ποιότητας και να μειώσουμε την επίδραση των κριτηρίων απόδοσης δεδομένου ότι για την ύδρευση απαιτούνται πολύ καλά ποιοτικά χαρακτηριστικά και σχετικά μικρές αποδόσεις. Οι βέλτιστες θέσεις στην περιοχή έρευνας τοποθετούνται στους σχηματισμούς των κροκαλοπαγών δυτικά της πεδινής περιοχής της Νεμέας και πέριξ των οικισμών Κρυονέρι και Σούλι καθώς επίσης στο ανάντη τμήμα του προσχωματικού υδροφόρου της Βόχας και πλησίον της κοίτης όλων των ποταμών της περιοχής έρευνας.