Groundwater vulnerability and delineation of protection zones in the discharge area of a karstic aquifer—application in Agyia’s karst system (Crete, Greece)Groundwater vulnerability and delineation of protection zones in the discharge area of a karstic aquifer—application in Agyia’s karst system (Crete, Greece) Peer-Reviewed Journal Publication Δημοσίευση σε Περιοδικό με Κριτές 2025-02-212023enThis work represents a contribution to the protection techniques of karst aquifers against groundwater pollution. The paper sets out the methodology being introduced for the protection of the karstic system that gives rise to five (5) major groups of springs and supplies fourteen (14) pumping wells near Agyia Chania (Crete, Greece). Starting from a geological and hydrogeological survey of the area, the work presents a vulnerability assessment of the karstic aquifer based on the application of three index-based methods (EPIK, PRESK and DRISTPI). The protection zones for the discharge area of the aquifer were delineated through an integrated geomorphological approach and groundwater flow modeling. At first, the risk of polluting substances migration from ground surface to groundwater was considered based on the spatial distribution of vulnerability. Following this, the vulnerability was evaluated in the saturated zone, where the attenuation mechanisms of contaminants were reducing due to the raised flow velocity. The groundwater flow and contaminant transport processes was considered using the MODFLOW code. Next, the data from the vulnerability mapping and the groundwater flow simulation were merged into an integrated assessment to delimit the protection zones for the water abstraction points. The vulnerability assessment outlines zones of high vulnerability in the SE part of the area, far away from the discharge zone of the aquifer and the water abstraction points. These zones are associated with an intensive infiltration process via carbonate formations. Protection Zone I was delineated 20 m around the water abstraction points, and it should be excluded from any anthropogenic activity. Protection Zone II coves part of the very high and high vulnerability zones defined by the DRISTPI method (located upwards of the water abstraction points), as well as an area downwards of springs and wells, where the flow path lines which demonstrate the subsurface travelling time of 50 days are projected to the ground surface. Protection Zone III extends outside Zone Ι and Zone ΙΙ, up to the limits of the hydrogeological or hydrological basin, whichever is larger. It includes the entire capture zone (i.e., the surface and underground catchment area) that feeds the water abstraction points. In this manner the protection zones include the entire contributing area to water abstraction points, not just the ground surface recharge zone.http://creativecommons.org/licenses/by/4.0/Water152Steiakakis_et_al_Water_15(2)_2023.pdfChania [Greece]Library of TUC2025-02-21application/pdf4.7 MBfree Steiakakis Emmanouil Στειακακης Εμμανουηλ Vavadakis Dionysios Βαβαδακης Διονυσιος Mourkakou Ourania Μουρκακου Ουρανια MDPI Karst aquifer Vulnerability Index-based methods EPIK PRESK DRISTPI Groundwater flow modeling 50-days isochronous curves Protection zones