Terrain Discontinuity Effects in the Regional Flow of a Complex Karstified Aquifer

Karst aquifers are characterized by spatial heterogeneity due to the presence of highly permeable channels and conduits in low-permeable fractured rocks (matrix block). Recent studies have reported a close relationship between surface and subsurface water in karstic regions due to the water flow through a complicated network of paths formed by fracture intersections. Subsurface flow in karstified aquifers ranges between conduit flow, in large passages with relatively high flow velocities, and diffuse flow, in the matrix block where Darcy’s law is still valid. In this paper, we present the simulation of a complex karstified aquifer system in Crete, Greece, where the presence of main faults drastically affects the regional flow. A discrete fracture approach in conjunction with an equivalent porous medium approach was adopted to simulate the mixed flow in the area of interest. The simulation results have shown that the length and the orientation of the dominant faults, primarily during the rainy season, affect the flow field.