Contribución de las técnicas hidroquímicas, isotópicas e hidrodinámicas a la caracterización de acuíferos carbonatados y su relación con la red de flujo superficial. Aplicación al acuífero de Alcadozo (cuenca del Segura)

Abstract

[SPA] El incremento de las demandas de agua en regiones áridas y semiáridas y los escenarios de cambio climático que proyectan un paulatino descenso pluviométrico en nuestra latitud, aumentan el interés en estudiar adecuadamente todos los acuíferos existentes en una zona, incluso aquellos poco explotados y conocidos, y a evaluar los recursos que pueden aportar al conjunto del sistema hídrico. En este contexto, esta tesis doctoral ha tenido como objetivo general demostrar que el uso de múltiples técnicas de caracterización e investigación hidrogeológica a distintas escalas espaciales, y la integración de la información que aportan todas ellas, es una metodología contundente para generar un modelo conceptual sobre el funcionamiento de acuíferos carbonatados de grandes dimensiones de los cuales existe poca información. Otro objetivo principal ha sido avanzar en la integración de algunos métodos de cuantificación de la recarga con el fin de mejorar la confiabilidad de los valores que proporcionan, de la variabilidad espacial de los mismos y de la información sobre los procesos que la controlan. El área de estudio es el acuífero de Alcadozo, que está localizado en el centro de la provincia de Albacete y es parte del dominio Prebético Externo de las Cordilleras Béticas. Con una superficie de 509 km2, el acuífero es parte de la MASb Alcadozo, constituida por materiales del Jurásico (Dogger y Lías), Cretácico Inferior y Paleógeno. Las calizas y dolomías jurásicas constituyen las principales formaciones permeables sobre el impermeable de base regional, que es el Triásico, y forman el acuífero estudiado. Tectónicamente es una zona de pliegues y escamas. El clima es mediterráneo semiárido, salvo en su sector occidental, donde hay una marcada influencia continental. El relieve es accidentado (cotas entre 944 y 1577 m s.n.m.) y presenta un fuerte gradiente pluviométrico (600 mm en el O a 360 mm en el E) y termométrico (12 ºC en el O a 17 ºC en el E). Su principal eje de drenaje superficial es el río Mundo, un río de montaña mediterránea calcárea que es afluente del río Segura. La explotación de agua subterránea es mínima en la mayor parte de la MASb y sólo es algo relevante en la zona oriental, donde hay agricultura de regadío. El objetivo general de la tesis se ha desglosado en varios objetivos particulares: definir el esquema regional del flujo de agua subterránea; cuantificar la recarga y conocer su distribución a escala regional, así como los rangos de incertidumbre de las estimaciones; caracterizar hidroquímica e isotópicamente las aguas subterráneas y superficiales, identificar las fuentes principales de solutos y deducir los procesos hidrogeoquímicos más relevantes a escala regional; analizar el impacto que tuvo la perforación del túnel de Talave sobre la red de flujo de agua subterránea y conocer la relación hidrodinámica actual con el macizo; elaborar un modelo hidrogeológico conceptual sobre el funcionamiento del acuífero y validarlo contrastándolo con modelación numérica del flujo de agua subterránea. Para conseguir estos objetivos, la metodología de investigación ha consistido en la aplicación de diferentes técnicas de estudio (hidroquímicas, isotópicas, hidrodinámicas, estadísticas y de modelación numérica) a un conjunto de datos empíricos procedentes, principalmente, de trabajos de campo realizados específicamente para esta tesis (2010 a 2014), en segundo lugar, de la tesis de maestría del autor (2008 a 2009) y, en menor medida, de estudios e informes previos (del Instituto Geológico y Minero de España, la Confederación Hidrográfica del Segura, Confederación Hidrográfica del Tajo y el Sindicato Central de Regantes del acueducto Tajo-Segura). Para conocer la configuración de la red de flujo a escala regional y la ubicación y magnitud de las zonas de recarga y descarga se han usado datos hidrodinámicos (niveles y caudales) medidos en diferentes manantiales, sondeos y tramos del río Mundo, así como en los piezómetros y drenes relacionados con el túnel de Talave. Esos datos también han permitido caracterizar hidrodinámicamente los dos manantiales principales del acuífero, Ayna y Liétor. Se ha realizado un trabajo específico en el interior del túnel de Talave para localiza, muestrear y estimar las principales filtraciones de agua subterránea. Se ha estudiado la relación río Mundo-acuífero desde su origen hasta el embalse de Talave mediante balances diferenciales de agua, de masa de soluto (cloruro) y de actividad de radón. La recarga de agua subterránea por lluvia ha sido estimada en condiciones estacionarias a escala subregional (cuenca hidrogeológica de los manantiales de Ayna y Liétor) y regional (acuífero) mediante los métodos de balance de agua en el suelo (BAS, modelado con el código Visual Balan) y balance de cloruro atmosférico (BMC) respectivamente. Se han calculado los intervalos de confianza (incertidumbre) asociados a los valores calculados y la consistencia de ambos métodos. Para caracterizar hidroquímicamente las aguas superficiales y subterráneas se han tomado unas 150 muestras de agua de ríos, manantiales, pozos y sondeos y se han usado algunos datos históricos de la literatura. También se han muestreado las principales surgencias interiores del túnel de Talave (una sola vez). Para conocer la composición química del agua de lluvia y el aporte de solutos al agua subterránea por vía atmosférica, entre los años 2011 y 2014 se ha muestreado y analizado química e isotópicamente el agua de lluvia en 4 estaciones que cubren razonablemente toda la superficie del área de estudio. En todas las muestras de agua de cualquier procedencia se han analizado los componentes mayoritarios; en muchas muestras se ha analizado Br y en unas pocas Sr. Se ha evaluado la consistencia de los análisis y se ha identificado la existencia de una incertidumbre analítica significativa, aunque aceptable a los fines de la mayor parte de las elaboraciones realizadas. La identificación del impacto de la actividad agrícola se ha apoyado analizando en algunas muestras la relación 87Sr/86Sr. Los datos se han estudiado con gráficos, análisis bivariante y multivariante (Análisis de Componentes Principales, ACP). Para identificar las zonas de recarga y la procedencia de las líneas de flujo, en todas las muestras de agua de lluvia y en un amplio conjunto de muestras de agua subterránea se han analizado los isótopos estables de la molécula de agua (18O y 2H). En un grupo reducido de muestras se ha medido la actividad de tritio (3H) para obtener información sobre los tiempos de permanencia del agua subterránea y la posible existencia de mezcla de líneas de flujo. La composición isotópica del agua de lluvia de la zona se ha contrastado con datos de diferentes estaciones nacionales de la red REVIP-AEMET-CEDEX. Con los datos hidrodinámicos y de caudal se estableció un modelo hidrogeológico conceptual inicial de la MASb Alcadozo que permitió conocer la configuración regional del flujo, las principales formaciones permeables y su disposición en la estructura geológica. Los trabajos iniciales también aportaron algunos valores de recarga, descarga natural (concentrada y difusa), explotación de agua subterránea y coeficientes de agotamiento de los manantiales principales. Esta información muestra un acuífero con una elevada capacidad de regulación y con un balance hídrico equilibrado. A esta información se sumó después la recolectada y generada durante el estudio del túnel de Talave, que aumentó notablemente la información disponible para el diseño del modelo conceptual. El análisis de la evolución hidrodinámica temporal de las descargas de agua subterránea al túnel muestra que la relación acuífero-túnel está ya en equilibrio, y que las oscilaciones temporales dependen principalmente de la recarga estacional. De ese análisis se ha estimado un descenso del nivel piezométrico de entre 69 y 104 m en un amplio sector del acuífero debido a la construcción del túnel. La cuantificación de la recarga a escala regional mediante BMC muestra una notable variabilidad espacial en las tasas de recarga, entre 20 y 243 mm/año, y coeficientes de variación entre 0,16 y 0,38. A escala subregional, los resultados obtenidos mediante BAS para las cuencas de Liétor y Ayna (calibrados con caudales de manantiales) han sido de 35 y 50 mm/año respectivamente, con coeficientes de variación de 0,49 y 0,36. Los resultados obtenidos a ambas escalas espaciales son consistentes entre sí. Se ha confirmado un sistema en régimen natural. Las aguas de lluvia son de dos tipos químicos, HCO3-Mg y HCO3SO4-Mg. Presentan concentraciones pequeñas de los componentes principales del aerosol marino Cl y Na, cuyo único origen es el marino. Existe una gradiente positivo de la concentración de Cl en el aporte atmosférico de NO y N a E y SE. Los contenidos de SO4 muestran una o más fuentes de sulfato, adicionales al aerosol marino, en el aporte atmosférico. Las fuentes potenciales son: la resuspensión de polvo litológico (afloramientos del Triásico), el aporte de polvo sahariano y de sulfuros de las emisiones del complejo industrial de Puertollano (a unos 170 km al O). Las aguas subterráneas tienen mineralizaciones bajas y medias (valores de CE entre 300 y 900 S/cm, pero más frecuentemente entre 400 y 600 S/cm) y son de tipo HCO3-CaMg y HCO3-MgCa, en coherencia con la litología dominante. En general, no presentan variaciones temporales significativas en su composición química. Una vez infiltrada en el terreno, la composición química del agua evoluciona por contacto con los minerales más abundantes del acuífero, calcita y dolomita. En zonas muy concretas del acuífero las aguas tienen la marca litológica de los sedimentos Triásicos, que afloran o subafloran puntualmente en el centro, y regionalmente en el O y SE. En el sector oriental del acuífero los excedentes de riego son también una fuente de agua y de solutos, incorporando al acuífero NO3. La composición química de siete muestras de agua tomadas en el interior del túnel de Talave muestra también el impacto de la recarga de excedentes de riego en esa zona, ya que todas tienen más NO3 del atribuible al fondo natural (entre 14 y 30 mg/L). Las concentraciones de Sr y los valores de la relación 87Sr/86Sr en un conjunto de muestras seleccionadas son consistentes con este modelo conceptual: en la mayor parte del acuífero las aguas tienen la huella química e isotópica de los carbonatos (Jurásico), en puntos localizados tienen la huella de las evaporitas del Triásico y en muestras del sector oriental muestran una mezcla de las huellas de los carbonatos y los fertilizantes. El ACP ha confirmado los principales procesos responsables de la composición química del agua subterránea y de su distribución regional. La composición química de las aguas del río Mundo sugiere notable interacción río-acuífero y transferencia de la marca química del agua subterránea al río. Los valores de δ18O y δ2H en el agua de lluvia son consistentes, en general, con los de la línea meteórica media mundial (LMM; δ2H=8δ18O+10), si bien una parte de las muestras está enriquecida en distinto grado por evaporación, lo que podría suponer valores originales algo más ligeros de los medidos. Los valores del exceso de deuterio indicarían que las precipitaciones más abundantes son de origen atlántico, pero el enriquecimiento evaporativo podría enmascarar una mayor contribución de lluvias de origen mediterráneo. Las aguas subterráneas medidas tienen rangos amplios, entre -6,60 y -8,47 ‰ para el δ18O y entre -45,8 y -56,14 ‰ para el δ2H. Esto indica que la recarga se produce a distintas altitudes. El gradiente isotópico altitudinal estimado para el δ18O con datos de manantiales y pozos someros es aproximadamente de -0,45 ‰/100 m, mayor que en otras zonas de estudio cercanas. Las lluvias tomadas durante cuatro años en tres estaciones con un rango de cotas de casi 600 m no muestran gradiente isotópico alguno, pero el conjunto de las aguas subterráneas analizadas sí muestra la impronta de un gradiente isotópico en la recarga que es consistente con el arriba mencionado. Esto sugiere que durante el periodo de observación han dominado las lluvias convectivas en vez de las frontales u orográficas, pero que a largo plazo dominan las segundas. Considerando válido dicho gradiente isotópico, la mayoría de las aguas subterráneas estudiadas, cuyas cotas de muestreo oscilan entre algo menos de 600 y algo más de 1400 m s.n.m., se habrían recargado en las zonas más altas del sistema, por encima de los 1200 m s.n.m. Las muestras de las filtraciones al túnel de Talave tienen dos rangos de cotas de recarga, uno por encima de los 1200 m s.n.m., consistente con los relieves más altos situados a unos 20 km al O del túnel, y otro algo más bajo, consistente con los relieves existentes en el entorno del túnel. Todas las muestras de agua subterránea medidas tienen contenidos de 3H que indican tiempos de residencia entre relativamente cortos y contemporáneos (entre 0,4 UT ± 0,2 UT y 5,4 ± 0,6 UT). El conjunto de los valores sugiere la existencia de dos partes relativamente diferenciadas en la red de flujo, una con flujos más rápidos y cortos, situados en zonas más altas, y otra con flujos de mayor extensión espacial y tiempo de residencia, situada en zonas más bajas. En cualquier caso, en el segundo grupo es factible que existan mezclas de líneas de flujo de distinto tiempo de permanencia. Las aguas filtradas al túnel son mezcla, en distintas proporciones, de líneas de flujo locales y subregionales. La moderada mineralización general de las aguas subterráneas, la aparición de NO3 en pozos agrícolas profundos y en las filtraciones al túnel de Talave y la presencia de tritio contemporáneo en muchas muestras indican que las aguas subterráneas estudiadas tienen tiempos de permanencia en el sistema relativamente cortos. Dado que los espesores de zona no saturada son notables (10 m al O a >120 m al E), eso implica que la infiltración del agua de lluvia hasta la zona saturada es rápida e importante. Además, los isótopos y la química muestran una mayor relevancia de las mezclas verticales de la prevista inicialmente, singularmente en el sector oriental del acuífero. Aunque a escala regional domina la componente horizontal del flujo, la escala espacial de las líneas de flujo es más bien subregional, y localmente dominan incluso los flujos locales. Este modelo conceptual sobre el conjunto del acuífero podría referirse a la parte superior del mismo, bien por ser esta la zona más permeable y conductora del conjunto (donde tendría lugar el flujo activo y la renovación), o bien por ser donde se ubican los manantiales y quizás también la mayoría de los pozos y sondeos muestreados. Confirmar esta hipótesis queda pendiente para futuros trabajos. La modelación numérica del flujo se ha realizado con el código MODFLOW. El modelo ha reproducido razonablemente el modelo conceptual, si bien es mejorable, especialmente en la definición de las condiciones de contorno. Una primera modelación en régimen permanente de las condiciones actuales proporciona niveles que se ajustan bastante bien a los observados, aunque algunos valores sugieren sectores con funcionamiento más complejo. Los balances zonales han sido coherentes con el balance hídrico de la MASb y reproducen adecuadamente los caudales observados (R2=0,97). En un segundo modelo (también en régimen natural) se han simulado las condiciones previas a la construcción del túnel de Talave. Las isopiezas resultantes son coherentes con el previsible esquema de flujo regional, y los caudales de los principales manantiales muestran algunas diferencias significativas con los valores actuales. El análisis tentativo del impacto hidrogeológico del túnel muestra una buena relación entre los valores calculados y los valores históricos disponibles, salvo para el manantial de Liétor. Comparando ambos modelos, el descenso estimado de la superficie piezométrica en el sector oriental oscila entre 30 m y 190 m.

[ENG] The increase of water demands in arid and semi-arid areas and the climate change scenarios that project a gradual decrease in rainfall in our latitude, augment the need to adequately study all aquifers in an area, even those littleexploited and known, and evaluate the resources they can contribute to the whole water resources system. In this context, the general objective of this thesis is to demonstrate that the use of multiple characterization techniques and hydrogeological research at different spatial scales, and the integration of the information they all contribute, is a conclusive methodology to stablish solid conceptual models on large carbonate aquifers of which there is little information. Another main objective is to advance in the integration of recharge estimation methods in order to improve the reliability of the values they provide, their spatial variability, and the information on the processes involved. The study area is the Alcadozo aquifer, which is located in the centre of the Albacete province and is part of the External Prebetic domain of the Betic Range. With an area of 509 km2, the aquifer is part of the Alcadozo Groundwater Body (GWB), which is formed by Jurassic (Dogger and Lias), Lower Cretaceous and Paleogene materials. The Jurassic limestones and dolomites constitute the main permeable formations, which overlay lowpermeability Triassic materials and form the aquifer studied. The area is characterized by a complex geologic structure consisting of a succession of folds and thrusts. The climate is semi-arid Mediterranean except in the western sector, where there is a marked continental influence. The relief is rugged (altitude between 944 and 1577 m a.s.l.) and display strong rainfall (600 mm to 360 mm from W to E) and temperature gradients (12 ºC in the W to 17 ºC in the E). The Mundo River is the main surface drainage axis. It is a calcareous Mediterranean mountain river and a tributary of the Segura River. Groundwater is almost unexploited in most of the GWB except in the eastern area, where there is irrigated agriculture. The general objective of the thesis has been split into several specific objectives: to define the regional groundwater flow; to estimate recharge rates and their spatial distribution at regional scale, as well as the uncertainty ranges of the estimates; to characterize groundwater and surface water hydrochemical and isotopic composition; to identify the main solute sources and deduce the most relevant hydrogeochemical processes at regional scale; to analyse the hydrodynamic impact that drilling of the Talave tunnel had on the groundwater flow and to know the current hydrodynamic relationship with the rock massif; to develop a conceptual hydrogeological model on the functioning of the aquifer and to validate it with groundwater flow numerical modeling. The methodology used to achieve these objectives consisted in the application Spring and river flow rates and borehole and wells piezometric levels have been used to know the groundwater flow network pattern at regional scale, the magnitude of recharge and the location of recharge and discharge areas. Those data also allowed understanding the hydrodynamic behaviour of the two man outlets of the aquifer, the Ayna and Liétor springs. A particular survey was carried out in the Talave tunnel to locate and sample the main groundwater leaks inside the tunnel. The relationship between the Mundo River and the aquifer has been studied by means of water balances, chloride mass balances and radon activity balances. Interannual recharge to the aquifer has been estimated for the Ayna and Liétor springs catchments by modelling the Soil Water Balance (SWB), and for the whole aquifer surface by Atmospheric Chloride Mass Balance (ACMB). The confidence intervals of the recharge values have also been estimated, as well as the consistency of both methods. To perform the hydrochemical and isotopic studies, some 150 water samples has been collected from rivers, springs, wells, boreholes and the main leaks inside the Talave tunnel. To study the supply of atmospheric solutes to the aquifer, rain water has been collected between years 2011 and 2014 in four rain stations. Major solute components have been analysed in all the water samples; Br has been measured in many samples and Sr has been measured in a selected set of samples. The internal consistency of chemical analyses has been evaluated and a significant uncertainty has been found, but it is acceptable for the purpose of most of the elaborations performed. The study of the impact of agricultural activities in groundwater chemistry has been supported with the measurement of the 87Sr/86Sr ratio in a selected set of samples. The chemical data have been studied by means of graphics, bivariate and multivariate analysis (Principal Component Analysis, PCA). To identify the main recharge areas the two stable isotopes of water, 18O and 2H, have been measured in all the rain water samples and most of the groundwater samples. To get information on groundwater residence time and the possible existence of mixing conditions, tritium activity has been measured in a set of selected samples. The isotopic composition of rain water in the area has been contrasted by comparison with the composition of rain water in several stations of the REVIP/AEMET/CEDEX network. The first piezometric and spring flow rate data allowed to stablish a preliminary conceptual model of the Alcadozo GWB, which informed about the regional groundwater flow pattern, the main permeable geological formations and their location within the geological structure. The first works also provided some values on recharge and discharge rates, groundwater exploitation, and depletion coefficients of the main springs. This information showed that the aquifer has a notable regulation capacity and its balance is equilibrated. The information available was later augmented with the many reports and data collected and generated during the study of the Talave tunnel, which contributed notably to the review and final design of the aquifer’s conceptual model. The study of the temporal evolution of the tunnel discharge shows that the tunnel is already in equilibrium with the aquifer. They analysis suggest that the piezometric drawdown produced by the drilling could be between 69 and of different study techniques (hydrochemical, isotopic, hydrodynamic, statistical and numerical modeling) to empirical data mainly generated from field works carried out specifically for this thesis (2010 to 2014), and also taken from the Master Thesis of the author (2008 to 2009) and from previous studies and reports of different organizations (Geological Survey of Spain, Segura River Water Authority, Tagus River Water Authority and Tagus-Segura Aqueduct Irrigators Authority). 104 m, and that it affected to a wide area in the central and eastern sector of the aquifer. Recharge quantification at aquifer scale by CMB shows a notable spatial variability between 20 and 243 mm/yr, with variation coefficients between 0.6 and 0.38. At subregional scale, modeling the SWB for the Ayna and Liétor springs catchments provided recharge values of 35 and 50 mm/yr, respectively, with variation coefficients of 0.49 and 0.36. The results obtained at the two spatial scales are consistent. This confirmed that the aquifer is under natural conditions. Rain water is of the Mg-HCO3 and Mg-HCO3SO4 HCO3-Mg y HCO3SO4-Mg chemical types. They have small contents of the two main components of marine aerosol, Na and Cl. Seawater is their only origin. A positive gradient of the Cl content in the atmospheric supply has been observed from NO and N to SE and E. The SO4 contents point to one or more sources in addition to the marine one. The potential sources of SO4 are: resuspension of lithological dust, supply of Sahara dust, and supply of sulphide particles from the Puertollano industrial area (some 170 km to the O of the study site). Groundwater has low to medium mineralization values (EC between 300 and 900 S/cm, but most frequently between 400 and 600 S/cm) and shows two very close chemical types, CaMg-HCO3 and MgCa-HCO3, which are consistent with the main lithology. Quite locally some spring and borehole samples have the chemical signature of the Triassic sediments, which crop out locally in the centre of the area and to the W and SE. In the eastern sector of the aquifer, including the leaks to the Talave tunnel, groundwater has agricultural NO3. In most of the area groundwater composition does not show temporal changes. The Sr contents and the 87Sr/86Sr values supports the conceptual model of the hydrochemistry origin: in the major part of the aquifer groundwater has the chemical signature of the Jurassic carbonates; in a few localized sites it has the signature of the Triassic evaporates, and in the eastern sector it shows the signature of a mixture between carbonates and fertilizers. The results of the PCA confirmed the main processes responsible of groundwater chemical composition and its regional distribution. The chemical composition of the Mundo River suggests a notable interaction with the aquifer, with produces a transfer of the aquifer composition to the river water. The δ18O y δ2H values of rainwater are consistent, in general, with the mean meteoric water line, although some samples are enriched by evaporation and their original values could be somewhat lighter than measured. The deuterium excess values point to a dominance of the Atlantic-origin precipitations, though the evaporative enrichment could mask a more relevant than observed Mediterranean contribution. Groundwater shows wide isotopic ranges, between -6.60 and -8.47 ‰ for δ18O and between -45.8 and -56.14 ‰ for δ2H, which points to recharge at different heights. The δ18O altitudinal isotopic gradient estimated with spring and handdug wells is around 0.45 ‰/100 m, larger than in other close areas. Rain samples collected during four years in three stations with an altitude range of almost 600 m do not show any isotopic gradient. However, the groundwater samples show the existence of an isotopic gradient in recharge water similar to the one mentioned above. This is suggesting that convective precipitations could have dominated during the study period while frontal precipitations dominate in the long term. Assuming valid the above mentioned gradient, most of the groundwater samples taken between 600 and 1400 m a.s.l. would have recharged above 1200 m a.s.l. The samples of the Talave tunnel show two recharge zones, one above 1200 m a.s.l,, compatible with an area located 19 km to the O of the tunnel, and another one lower, consistent with the topographic levels around the tunnel. The measured tritium contents (between 0.4 UT ± 0.2 UT and 5.4 ± 0.6 UT) indicate short residence times, with some samples compatible with contemporary recharge. The set of values belong to two zones of the groundwater flow network: the samples of the highest springs and wells are from a zone with rapid and short flow lines; the samples of the lowest springs and wells show somewhat longer flow lines and residence times. The second group could represent mixtures of groundwater lines with different residence times. Groundwater leaking to the Talave tunnel is mixtures of local and subregional flow lines. The general moderate mineralization of groundwater, the presence of agricultural NO3 in deep wells of the E part of the aquifer and in the leaks to the Talave tunnel, and the presence of contemporary tritium in most of the samples indicate that the studied groundwater have rather short residence times. Given that the thickness of the unsaturated zone is relevant (from 10 m to the O to >120 m to the E), this implies that rainwater infiltration to the saturated zone is a relevant and fast process. Moreover, the isotopic and chemical composition a most relevant role to vertical mixtures than expected, especially in the eastern part of the aquifer. This means that even though at regional scale the main component of groundwater flow is the horizontal one, the magnitude of the flows is subregional rather than regional, and local flows can dominate in particular zones. Given the large thicknesses of the Jurassic formations, this conceptual model could be related to the upper part of the aquifer and not to the whole aquifer. This could be explained in two ways: by the existence of a more permeable and conductive zone in the upper part, where the active flow and renovation would take place; by the fact many wells are not fully penetrants and many springs have small catchments. However, to check these hypothesis remain pending for a future study. The modeling of groundwater flow has been performed with MODFLOW. The model has reproduced reasonably the main features of the conceptual model, but it has to be improved. A first model of the present day situation in steady state conditions provided piezometric levels close to the observed ones, though in some areas the values suggest a more complex behaviour. The water balances performed by zones agree quite well with the balances estimated by other methods, and the modelled springs flow rates approach the observed ones (R2=0.97). A second model simulated the situation before the drilling of the Talave tunnel. The piezometric levels modelled are consistent with the expected regional flow pattern, and the springs flow rates are larger than the present day ones. However, those flow rates are consistent with the few historical data available, with the exception of the Liétor spring. A tentative analysis of the possible impact of the drilling shows that the piezometric drawdown produced by the tunnel in the oriental part of the aquifer could have been between 30 and 190 m.