La evaluación de la recarga al acuífero del Norte de Gran Canaria mediante el balance de agua en el suelo

   LA EVALUACIÓN DE LA RECARGA AL ACUÍFERO DEL NORTE DE GRAN CANARIA MEDIANTE LE BALANCE DE AGUA EN EL SUELO Fabio De Paoli Curso 2015/2016 TUTORA: María del Carmen Cabrera Santana COTUTORA: Tatiana Cruz Fuentes Trabajo Fin de Título para la obtención del Master Universitario en Gestión Costera    El trabajo de Fin de Título (TFT) del Máster Universitario de Gestión Costera de la Universidad de Las Palmas de Gran Canaria, titulado "La evaluación de la recarga al acuífero del norte de Gran Canaria mediante le balance de agua en el suelo", ha sido realizado en el curso 2015/16 por el alumno Fabio De Paoli, con NIE Y26841668, bajo la codirección de la Dra. María del Carmen Cabrera Santana (Opto. de Física) y la Dra. Tatiana Cruz Fuentes (colaboradora externa del grupo de investigación GEOVOL), ambas vinculadas al iUNAT. El trabajo se ha llevado a cabo a partir de los datos del Proyecto de Investigación REDESAC (CLG2009/1291 0.C03). Fdo. Fabio De Paoli � Carmen Cabrera Santana Fdo. Tatiana Cruz Fuentes Las Palmas de Gran Canaria, 7 de julio de 2016 3   4  ÍNDICE  1  INDICE Agradecimientos 3 RESUMEN 5 ABSTRACT 5 1. INTRODUCCIÓN 7 1.1. Antecedentes 11 1.2. Objetivos 12 1. MATERIALES Y MÉTODOS 13 1.1. Caracterización de la zona de estudio 13 1.2. Programas informáticos utilizados 16 1.3. El código Visual-Balan 16 1.4. Datos de entrada utilizados 18 2. RESULTADOS 27 3. DISCUSION 31 4. CONCLUSIONES 33 5. BIBLIOGRAFÍA 34  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  2  AGRADECIMIENTOS  3  AGRADECIMIENTOS Quiero agradecer de manera especial y sincera a la Profesora María del Carmen Cabrera Santana y directora de este trabajo,por la orientación, el seguimiento y la supervisión continúa y sobre todo por la motivación y el apoyo recibido a lo largo de toda esta experiencia. Especial reconocimiento merece el interés mostrado por mi trabajo y las sugerencias recibidas de la Dra. Tatiana Cruz Fuentes, codirectora de este trabajo,que con su ayuda colaboró en la realización del presente trabajo. También quiero agradecer sinceramente todas aquellas personas que compartieron conmigo sus conocimientos para hacer posible la conclusión de este trabajo. Especialmente agradezco al Profesor Emilio Custodio, a D. Luis Fernando Martin, a Dña. Diana Rodríguez, a la Dra. Esmeralda Estévez, a Dña. Gema Naranjo y al Dr. Alejandro Rodríguez.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  4  RESUMEN Y ABSTRACT  5  RESUMEN La explotación de las aguas subterráneas ha sido y todavía sigue siendo realizada de forma no controlada. La utilización de los recursos hídricos subterráneos de manera intensiva amenaza a su sostenibilidad y los ecosistemas que de ellos dependen. La recarga de los acuíferos es un fenómeno natural muy complejo. Su estudio abarca muchas simplificaciones que conllevan una notable incertidumbre, por ello hay que abordar su cálculo mediante diferentes técnicas, lo más independientes posible entre sí. El conocimiento de la recarga es básico para la correcta planificación y aprovechamiento de los recursos subterráneos de una zona. En el cuadrante Nororiental de la isla de Gran Canaria las actividades agrícolas representan el mayor gasto hídrico, mientras que esta zona recibe la mayor parte de la recarga insular. La recarga en esta zona ha sido estimada ya en diferentes estudios anteriores dentro del proyecto REDESAC mediante el método del balance del ión cloruro de origen atmosférico y del balance diario de agua en el suelo. En el presente trabajo se propone una nueva evaluación de la recarga mediante el programa VISUAL BALAN, que oscila entre 12% y 26% de la Precipitación, con una media ponderada por la superficie del 19% y su comparación con las estimaciones anteriores. ABSTRACT The exploitation of groundwater has been and is still being done in an uncontrolled manner. The intense use of groundwater resources threatens its sustainability and sustainability of ecosystems that depend on it. The recharge of the aquifers is a complex natural phenomenon. The study of this phenomenon covers many simplifications that involve significant uncertainty that is the reason why the calculation has to be done using different techniques, which have to be independent from each other as much as possible. The evaluation of groundwater recharge is essential for a proper planning and utilization of groundwater resources in an area. Agricultural activities is the largest water waste activity In the Northeast quadrant of the island of Gran Canaria, and this area is where the largest insular recharge occurs. Recharge of this area has already been estimated in several previous studies within the project REDESAC by the method of chloride atmospheric ion balance and daily water soil balance. In this paper a new assessment of recharge is proposed by VISUAL BALAN program, which ranges between 12% and 26% of the rainfall, with a surface weighted average of 19% and also its comparison with previous estimations.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  12  atmosférica de cloruro y métodos químicos e isotópicos ambientales los acuíferos estudiados. Dentro del proyecto REDESAC, entre los años 2008 y 2014, el grupo GEOVOL de la Universidad de Las Palmas de Gran Canaria, dirigidos por la Dra. María del Carmen Cabrera llevo a cabo varios trabajos con el objetivo de estimar la recarga natural en la zona noreste de la isla de Gran Canaria mediante el método del balance del ión cloruro de origen atmosférico. Como resultado del muestreo de agua de lluvia entre 2008 y 2013, se estimó una recarga media multianual estimada de 25,4±18 hm3/año, o sea el 22±16% de la precipitación (Naranjo et al., 2014). Para contrastar esta estimación y la del Plan Hidrológico de Gran Canaria, que resultó el 19% de la precipitación, el mismo grupo realizó un balance diario de agua en el suelo entre los años 1980 y 2013. Para este trabajo se utilizó la hoja de cálculo Easy-Bal 3.0, obteniendo una recarga de 15±4 hm3/a, lo que supone entre un 9% y 17% de la precipitación (Naranjo et al., 2015). 1.2. Objetivos En este trabajo se va a utilizar el código Visual-Balan desarrollado por el grupo de hidrología de la Universidad de A Coruña para cuantificar y evaluar la recarga en el Norte de Gran Canaria utilizando los datos obtenidos en el proyecto REDESAC y comparar los resultados con los trabajos anteriormente citados. Ello ayudará a comprender mejor el sistema, y detectar los datos más sensibles a las variaciones que deberían ser estudiados para mejorar la planificación hidrológica de la zona. MATERIALES Y MÉTODOS  13  2. MATERIALES Y MÉTODOS 2.1. Caracterización de la zona de estudio El área de estudio ocupa las zonas N2, N3 y N4 de la zonificación hidrogeológica de Gran Canaria (PHGC, 1999) (Figura 4). Se extiende por la costa, desde las inmediaciones de Las Palmas de Gran Canaria al este, hasta el litoral que pertenece al municipio de Moya en el noroeste. Hacia el interior la zona penetra hasta las proximidades del Pico de las Nieves (1.949 m de altitud), siguiendo las divisorias hidrográficas. La zona tiene una superficie de 312 Km2. En superficie afloran los materiales más jóvenes y permeables, con suelos pocos retentivos y vegetación de moderada densidad, lo que favorece la recarga. El elevado gradiente altitudinal y la dirección predominante de los vientos alisios Isoyetas medias (mm/a) periodo 1970-2006 Figura4:Situacióndelazonadeestudioconindicacióndelascuencashidrográficasconsideradasylas isoyetasanualesmediascalculadasentre1970y2006.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  14  hacen que el clima en la zona sea muy variable con una precipitación media anual varía entre 800 mm/a a cotas altas y 115 mm/a a cotas bajas. Las temperaturas máximas y mínimas medias anuales varían desde los 14˚C en las zonas altas a 21˚C en la costa, con una temperatura media de 18˚C (Naranjo et al., 2015 a y b ). Desde un punto de vista geológico, la Isla de Gran Canariaforma parte de un archipiélago volcánico situado en la placa africana y están formada por acumulación de materiales efusivos sobre corteza oceánica como producto de la actividad de un punto caliente (Pérez-Torrado, 2008). En la zona de estudio afloran fundamentalmente lavas e ignimbritas Pliocenas y Cuaternarias y sedimentos de la Formación Detritica de Las Palmasde edad Miocena y Pliocena (Figura 5).   Figura5:GeologíasuperficialsimplificadadeGranCanaria(Pérez‐Torrado,2008). MATERIALES Y MÉTODOS  15  Según estudios hidrogeológicos anteriores de la zona (Naranjo y López-Galayo, 2006; Hernández-Quesada, 2016; Estévez, 2016), se puede deducir que la mayoría de los pozos explotan los materiales Pliocenos del Grupo Roque Nublo y los fonoliticos miocenos infrayacentes. Según el modelo hidrogeológico aceptado para las islas volcánicas como Canarias (Figura 6), la isla se caracteriza por la existencia de un acuífero único con una superficie piezométrica en forma de domo (SPA-15, 1975; Cabrera y Custodio, 2004). Dentro de este esquema, la recarga natural se genera en las zonas de cumbre y medianía por infiltración de la lluvia. El agua circula preferentemente por los materiales volcánicos más recientes y permeables hacia la costa. La descarga se produce al mar y mediante la extracción de pozos y galerías que han ido deprimiendo las descargas por manantiales (Hernández-Quesada, 2016).  Figura6:PerfilesquemáticodelaevoluciónhidrogeológicadeGranCanariaalser sometidaaexplotación(modificadadeCustodioyCabrera,2002).  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  16  2.2. Programas informáticos utilizados Durante el desarrollo del presente trabajo se han utilizado diferentes programas, en mayor o menor medida. Los datos han sido recogidos en bases de datos en Microsoft Access, mientras su manejo y la realización de cálculos específicos se han llevado a cabo con Microsoft Excel. Estos programas han sido fundamentales para el tratamiento e introducción posterior de los datos por el código Visual BALAN v.2.0. La captura, edición, análisis, tratamiento y diseño de la información geográfica ha sido realizado con las aplicaciones de ArcGis Desktop. 2.3. El código Visual BALAN El código Visual-BALAN ha sido desarrollado por el grupo de Hidrología Subterránea de la Universidad de A Coruña, cuya primera versión vio la luz en 1999(Samper et al., 1999). A partir de dicho año, el equipo de trabajo ha seguido actualizando y enriqueciendo el mismo con más opciones de forma que puede reproducir mejor la realidad. La versión utilizada en el presente trabajo es la v.2.0, cuyo manual de instrucciones puede ser descargado directamente de la web (ftp://ceres.udc.es/ITS_Caminos/Optativas/Hidrologia_Subterranea/VISUAL%20BALAN/Vi sual_Balan/manual%20Visual%20Balan%20v2.pdf) Dicho programa realiza balances de agua en el suelo edáfico, en la zona saturada y en el acuífero utilizando un entorno interactivo de entrada de datos y salida de resultados. Los términos fundamentales del balance son las entradas por precipitación y riegos, las salidas por escorrentía superficial, evapotranspiración, intercepción, flujo hipodérmico y flujo subterráneo, y la variación del contenido de humedad del suelo y de la zona no saturada y del nivel de agua en el acuífero (Figura 7). Figura 7 :EsquemaconceptualdelbalancehídricosegúnelManualdelprogramaVisualBalanv.2.0. MATERIALES Y MÉTODOS  17  El programa VISUAL BALAN V.2.0 se estructura en cuatro fases: A, B, C y D. En la Figura 8 se representa el esquema general del funcionamiento del programa. En la primera fase (A) se lee la definición propia del balance: situación geográfica de la zona, características del suelo, zona no saturada y el acuífero, y los datos climatológicos de las estaciones de medida utilizadas para una única serie de años hidrológicos. También se leen las opciones de los cálculos de las diferentes componentes y las opciones de salida de los resultados. La segunda fase (B) comprende las determinaciones previas para el cálculo de la evapotranspiración potencial (ETP), y se realiza eligiendo el método entre los propuestos, que son Thornthwaite, Blanney – Cridlle o Penman. Se introducen los valores de las diferentes características climatológicas, tales como el número teórico de horas de sol diarias de cada mes y el índice de calor mensual para el método de Thornthwaite; el porcentaje del número máximo de horas de insolación mensual respecto al total anual para el método de BlanneyCriddle; la radiación global para el método de Makkink; y la radiación extraterrestre y corrección por la altura del anemómetro para el método de Penman.En la tercera fase (C) se realiza el balance hidrometeorológico propiamente Figura8EsquemageneraldelfuncionamientodelprogramaVisualBalanv.2.0.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  18  dicho. Para cada día, el programa realiza una serie de cálculos de forma secuencial. En primer lugar, y sólo en el caso de que haya riego, calcula el volumen de agua aplicada al suelo. A continuación calcula el valor de la evapotranspiración potencial mediante uno de los cuatro métodos implementados. Posteriormente VISUAL BALAN calcula la intercepción, la precipitación nival, la fusión y escorrentía nivales, y el flujo preferente. Luego determina la escorrentía superficial, la evapotranspiración real y el flujo de salida de la zona edáfica. Una vez realizado el balance en el suelo, se efectúa el balance en la zona no saturada, calculando el flujo hipodérmico y la recarga al acuífero. Finalmente, en el balance en el acuífero se determina el caudal de salida y la variación del nivel piezométrico. En la cuarta fase (D) se presentan mediante un interfaz gráfico todos los resultados del balance para cada subcuenca y para la cuenca total, permitiendo al usuario visualizar todos los resultados del programa que incluyen: los resultados diarios, mensuales y anuales, el análisis de sensibilidad, los valores acumulados y medios en el período, los niveles medidos y calculados; los caudales medidos y calculados; las funciones de recarga y/o goteo y los resultados del ajuste automático. 2.4. Datos de entrada utilizados La entrada de datos del programa Visual BALAN se divide en dos fases principales:  DATOS BASICOS o División de la cuenca o Datos de las estaciones o Datos del suelo edáfico o Datos climatológico (calculo ETP) o Calculo de la evapotranspiración real o Método del cálculo de la escorrentía  OPCIONALES o Zona no saturada o Parámetros del acuífero o Calculo de la intercepción En la primera fase de datos básicos, es donde se concentró la mayor parte del trabajo de la investigación. Después de una esquemática descripción de la zona y de sus subcuencas, se decidió dividir la zona en las tres cuencas hidrográficas que la conforman, llamadas por la planificación hidrológica de la isla como N2 N3 y N4 (Figura 9). En la zona de estudio hay tres redes diferentes de estaciones meteorológicas: la red del Consejo Insular de Agua de Gran Canaria (CIAGC), la red de la Agencia Estatal MATERIALES Y MÉTODOS  19  de Meteorología (AEMET) y la red de la Consejería de Agricultura del Gobierno de Canarias. Para empezar el trabajo se llevó a cabo una selección de las estaciones más idóneas, atendiendo a su situación espacial y a la cantidad y calidad de los datos disponibles. Las estaciones consideradas, cuya localización espacial se muestra en la Figura 9 son las siguientes:  Estaciones CIAGC o Hoya de Gamonal 002 o Las Madres 055 o Toscón-Tamaraceite 088 o Retamilla 098 – Montañon 174 o Lomo San Pedro 199 o Campo de Golf 213  Estaciones AEMET (Agencia Estatal de Meteorología) o Bañadero C669B o Tafira C658X o C Tejeda C612F o Osorio C656V o Valleseco C665T  Estaciones Agroclimáticas (Gobierno de Canarias) o Arucas GC05 o San Mateo GC04 Figura9:Divisióndelazonaenlastressubcuencasconsideradas.Lafiguramuestratambiénlasituaciónespacialde lasestacionesmeteorológicasqueexistenenlazonaydelospozosconmedidasdenivelestáticosconsiderados.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  20  Los valores meteorológicos diarios requeridos, para cada estación, son:  Precipitaciones (mm)  Temperatura (ºC)  Concentración agua de lluvia (mm/l)  Número real de horas de sol  Velocidad media viento  Humedad relativa (%)  Indicador entre viento diario y nocturno  ETP usuario en (mm/día) Para el programa son variables imprescindibles las fechas, la precipitación y la temperatura. Para las demás variables se puede poner datos mensuales o simplemente no disponer de datos. Los datos diarios, se organizaron vinculados a la primera columna, donde aparece la serie temporal como se muestra en la Figura 10. En esta columna no puede faltar ninguna fecha, para el correcto funcionamiento del programa y debe coincidir para cada estación meteorológica, de aforos y de niveles. Figura10EjemplodehojaExceldeunaestaciónmeteorológica. MATERIALES Y MÉTODOS  21  Ninguna de las estaciones existentes es capaz de suministrar la totalidad de las variables necesarias para ejecutar el programa. Por ello, se mezclaron datos de diversas procedencias, teniendo en cuenta la cercanía, altitud y orientación de las distintas estaciones escogiendo la base de datos de REDESAC como punto de partida. De esta base se han obtenido datos de lluvia y temperatura desde el año 2008 hasta el 2013, procedentes respectivamente de las estaciones pluviométricas del Consejo Insular de Aguas y de las estaciones AEMET seleccionadas por el proyecto. A estos mismos dos organismos se les pidió la actualización de estos datos hasta el 2015. No se pudo llegar al 2016 ya que algunas de las variables consideradas no estaban disponibles para esas fechas. Se obtuvieron también los datos diarios de insolación total, velocidad y dirección media del viento y humedad relativa de las estaciones de la Agencia Estatal de Meteorología (AEMET). Estos datos fueron suministrados por la oficina AEMET presente en el campus universitario de la Universidad de Las Palmas de Gran Canaria. El año 2008 se tuvo que elegir como límite inferior de la serie temporal ya que estas estaciones empezaron a recoger datos de viento y de humedad relativa desde esta fecha. La insolación diaria, medida en porcentajes en cuatro momentos distintos del día (00:00h, 07:00h, 13:00h, 18:00h), tuvo que pasarse a horas reales de sol, considerando que mediamente un día en canarias tiene 12 horas de luz. Esta variable resulta medida exclusivamente en la estación C658X, de Tafira. Estos datos fueron extrapolados a las demás estaciones seleccionadas. Los datos de viento de la red, se miden en horas climatológicas. En cuatro distintos momentos del día (00:00h, 07:00h, 13:00h, 18:00h) se mide la velocidad media (km/h) y la dirección media (en decenas de grado) sobre 10 minutos. En las hojas Excel solo se puso la media de los cuatro valores de velocidad media obtenida diariamente. Para la humedad relativa diaria (%), se calculó el valor medio entre la humedad mínima y máxima medidas en las estaciones. Las dos estaciones Agroclimáticas del Gobierno de Canarias sirvieron para comparar y validar los datos de las demás estaciones. Los datos meteorológicos en las dos subcuencas que contienen más de una estación (N2 y N3), han sido obtenidos mediante ponderación según la superficie de influencia de cada estación, de manera que en cada subcuenca se han introducido un solo conjunto de datos. La selección de los pozos utilizados como puntos de observación de los niveles piezométricos se llevó a cabo a partir de la base de datos del proyecto REDESAC, que comprende una notable cantidad de datos de pozos en “Terrenos Particulares”. De los pozos sin explotación localizados en la zona de estudio (representados en la Figura 9) presentes en el área de estudio, se seleccionó uno para cada subcuenca. Alrededor de los pozos en explotación se genera una depresión piezométrica, que afecta a los resultados. Por ello, se seleccionaron los pozos que llevaban más tiempo abandonados, de forma que los niveles se asemejaran más a niveles estáticos. Los pozos escogidos (Figura 11) son:  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  28  Los resultados finalmente obtenidos en mm/año y hm3/año se muestran en la Tabla 10 y los porcentajes de los diferentes parámetros respecto a la Precipitación total en los años considerados para cada zona se muestran en la Tabla 11. Tabla10:Resultadosdelbalancehidráulicoobtenidosparalasdistintassubcuencasenlazonadeestudio,expresadosen mm/añoyhm3/año.P=precipitación;ETR=Evapotranspiraciónreal;R=Recarga;Int=IntercepciónyFHIP=Flujo hipodérmico. ZONA Área (km2) P ETR R Int F Hip P ETR R Int F Hip mm/año hm3/año N2 105.3 678.7 440.09 174.70 52.78 13.77 71.5 46.3 18.4 5.6 1.4 N3 151.4 541.2 404.10 92.99 44.12 6.42 81.9 61.2 14.1 6.7 1.0 N4 42.7 340.3 266.70 41.76 34.04 4.49 14.5 11.4 1.8 1.4 0.2 TOTAL 299.4 167.9 118.9 34.3 13.7 2.6 Tabla11:Resultadosdelbalancehidráulicoobtenidosparalasdistintassubcuencasenlazonadeestudio,expresadosen porcentajesrespectoalaPrecipitacióndecadasubcuencaenhm3/año.P=precipitación;ETR=Evapotranspiraciónreal; R=Recarga;Int=IntercepciónyFHIP=Flujohipodérmico. ZONA P (hm3/año) ETR (%) R (%) Int (%) F Hip (%) N2 71.47 64.84 25.74 7.78 2.03 N3 81.94 74.67 17.18 8.15 1.19 N4 14.53 78.37 12.27 10.00 1.32 Según estos resultados, la recarga al acuífero varía entre 12 y 26% según la subcuenca considerada, con un valor de recarga media ponderada según la superficie de cada subcuenca de 19%. En las Figuras 12, 13 y 14 se muestran los resultados medios anuales para las tres subcuencas obtenidos a través del interfaz del programa VISUAL BALAN de forma gráfica. RESULTADOS  29   Figura12:ResultadosobtenidosmedianteelprogramaVISUALBALANparalasubcuencaN2(enmm).  Figura13:ResultadosobtenidosmedianteelprogramaVISUALBALANparalasubcuencaN3(mm).  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  30   Figura14:ResultadosobtenidosmedianteelprogramaVISUALBALANparalasubcuencaN4(enmm). DISCUSIÓN  31  4. DISCUSIÓN El cálculo de la recarga a los acuíferos conlleva un gran número de factores que le otorgan un grado alto de incertidumbre. Esta incertidumbre proviene en primer lugar de la que tienen en sí mismos los datos de partida, desde los datos meteorológicos (precipitación, temperatura, horas de insolación, etc.) hasta los parámetros del suelo, zona no saturada y acuífero. En la zona de estudio, hay que tener en cuenta que existen solamente 8 perfiles del suelo, cuyos datos fueron extrapolados al total de la superficie de las tres subcuencas. Ello conlleva un alto grado de incertidumbre en la estimación del Punto de Marchitez Permanente (PMP), Capacidad de Campo (CC) y espesor del suelo que no ha podido ser evaluada. Además, hay que tener en cuenta que parte de la zona presenta pendientes elevadas y suelos poco desarrollados o directamente inexistentes, por lo que el espesor del suelo considerado representa una media estimada cuyo valor real no puede ser verificado. Por otro lado, los parámetros CRPG y CEPG, necesarios para el cálculo de la ETR por el método de Penman-Grindley fueron estimados por calibración porque no existen precedentes en la zona. Otra fuente de incertidumbre son los valores de conductividad hidráulica. En el caso del suelo, se estimaron valores a partir de los perfiles de que se disponía, presentando el mismo problema que los parámetros anteriores del suelo (CC y PMP). Para la zona no saturada y el acuífero, los valores considerados se obtuvieron mediante un modelo matemático (INTECSA, 1981), que no ha sido validados con datos de campo. Los resultados obtenidos son plausibles, aunque no pueden ser verificados dadas las características del acuífero: niveles piezométricos profundos que no responden rápidamente a la recarga y existencia de gran cantidad de pozos en explotación que impiden tener datos piezométricos estáticos, aun cuando existen pozos abandonados durante años con niveles piezométricos influenciados. Ello queda reflejado en la Figura 15, que representa los niveles medidos en el pozo 2250 TP (abandonado desde hace varios años), que supuestamente representa el nivel estático del acuífero. Sin embargo, se observa que los niveles están muy por debajo de los niveles calculados por el programa, lo que confirma la importancia de la explotación de las aguas subterráneas en el acuífero En los estudios anteriores se había dado mucho peso a las largas serie de datos disponible de lluvia y temperatura, contrariamente a este estudio. Así, gracias a las posibilidades del VISUAL BALAN, se tendió a organizar una base de datos temporalmente más corta (2008-2015) que pero englobara más variables, lo cual puede suponer una ventaja, aunque dependiendo de la calidad de los datos puede ser una desventaja.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  32  La estimación de la recarga en la zona de estudio según el método de la estimación de la recarga mediante el balance de cloruros de origen atmosférico es del 22% de la Precipitación (Naranjo et al., 2015a) y según un primer intento de balance de agua en el suelo varía entre el 9% y el 17% de la Precipitación (Naranjo et al., 2015b). Los resultados del presente estudio son coherentes con los obtenidos en los trabajos anteriores para la misma zona, lo cual es de esperar debido a que los datos de partida utilizados son los mismos ampliándose únicamente el tiempo considerado para los datos meteorológicos. Figura15:Nivelespiezométricoscalculadosporelprogramaymedidosenelpozo2250TP,representativo delazonaN2. CONCLUSIONES  33  5. CONCLUSIONES En el presente trabajo se ha calculado la recarga al sector nororiental de la isla de Gran Canaria mediante técnicas numéricas, aplicando el código VISUAL BALAN. Este programa permite establecer los balances hídricos en el suelo edáfico, la zona no saturada y el acuífero utilizando datos meteorológicos y parámetros de las zonas consideradas. La recarga media anual estimada para los años hidrológicos 2008/2015 para la toda la zona de estudio es de 34.3 hm3/año, lo que supone un porcentaje medio del 20% respecto a la precipitación. La zona de estudio fue dividida en tres subcuencas atendiendo a la división hidrológica del Consejo Insular de Aguas de Gran Canaria, observándose diferencias en los valores obtenidos entre las mismas. Así, los porcentajes de Recarga frente a la Precipitación varía el 12% en la zona más oriental (N4) y el 26% en la zona más occidental. (N2). Estos valores están dentro de los rangos establecidos en trabajos anteriores. Las principales dificultades encontradas durante el desarrollo del trabajo se refieren a la disponibilidad insuficiente de parámetros de partida medidos en campo para calcular la recarga que atiendan a la hetereogeneidad de los medios volcánicos. Ello hace necesario un profundo trabajo de caracterización básica del medio, mediante la toma de datos territorialmente distribuidos: ensayos de suelo en laboratorio, ensayos de bombeo en pozos, etc. Estos datos deberían ser tratados mediante tecnologías GIS de forma que representaran más fiablemente la realidad.  34  6. BIBLIOGRAFÍA Cabrera MC, Custodio E. (2004): Groundwater flow in a volcanic-sedimentary coastal aquifer: Telde area, Gran Canaria, Canary Islands, Spain. Hydrogeology Journal, 12: 305320. CoDo Curso Internacional de Hidrología Subterránea (2009): Hidrogeología. Conceptos básicos de hidrología subterránea. Ed. Fundación Curso Internacional de Hidrología Subterránea, Barcelona. 768 pp. Cruz Fuentes T. (2008): Contribución de la modelación numérica e hidroquímica al conocimiento del acuífero de La Aldea (Gran Canaria, Islas Canarias).Tesis doctoral ULPGC, 312 pp. Cruz-Fuentes, T., Naranjo, G., Cabrera, M.C., Custodio, E. (2014). Estimación de la escorrentía superficial mediante la generación automática del número de curva con un SIG. Aplicación al norte de Gran Canaria (Islas Canarias, España). II Congreso Ibérico de las Aguas Subterráneas, CIAS2014. Valencia: 217-232. Custodio E. (1997): Recarga a los acuíferos: aspectos generales sobre el proceso, la evaluación y la incertidumbre. En E. Custodio, M.R. Llamas y J. Samper (eds.): La evaluación de la recarga a los acuíferos en la planificación hidrológica: 19-40. Ed. Instituto Tecnológico GeoMinero de España/AIH-Grupo Español. Custodio, E. y Cabrera, M.C. (2002): ¿Cómo convivir con la escasez de agua? El caso de las Islas Canarias. Boletín Geológico Minero. Vol. 113, 3, 243-258. Custodio, E. y Llamas, M.R. (1983): Hidrología Subterránea. 2 Vol. Ed. Omega. 2º Edición. Estévez, E. (2016): Afección a las aguas subterráneas como consecuencia del uso de aguas regeneradas para riego en Gran Canaria. Tesis doctoral ULPGC, 140 pp. Heredia, J. y Murillo, J.M. (2002): Balance hídrico y estimación de la recarga mediante modelación numérica en pequeñas cuencas del sector sudoriental de la isla de Gran Canaria. XI Congreso Internacional de la Industria, Minería y Metalurgia, Zaragoza. 19 pp. Hernández-Quesada, M.P. (2016): Funcionamiento hidrogeológico del área de los Barrancos de Moya y Azuaje, Norte de Gran Canaria. Tesis doctoral ULPGC, 220 pp. + anejos. INTECSA (1981): Informe sobre el modelo matemático de flujo subterráneo de Gran Canaria. Excmo. Cabildo Insular de Gran Canaria. Informe interno, sin publicar. Naranjo, G. (2015a): Estimating Natural Recharge by Means of Chloride Mass Balance in a Volcanic Aquifer: Northeastern Gran Canaria (Canary Islands, Spain). Water 2015, 7, 2555-2574; doi:w10.3390/7062555. Naranjo, G.; Cabrera M.C.; Cruz-Fuentes T.; Martin L.S.; Custodio E. (2015b): Cuantificación de la recarga natural al acuífero del norte de Gran Canaria mediante el balance diario de agua en el suelo. En M.C. Cabrera et al. (eds): II Workshop: Estudio, aprovechamiento y gestión del agua en terrenos volcánicos, 247-254. BIBLIOGRAFÍA  35  Naranjo, G. y López-Galayo, O. (2006): Estudio hidrogeológico de las cuencas del Barranco de Guiniguada, Barranco de Tenoya e Interfluvio, del Sector Norte de Gran Canaria (Islas Canarias). Trabajo del 40 Curso Internacional de Hidrología Subterránea, 117 pp. PHCG (Plan Hidrológico de Gran Canaria) (1999): Boletín Oficial de Canarias 73, 08/06/1999. Decreto 82/1999. Pérez-Torrado, F.J. (2008): Geología de Gran Canaria. En: Itinerarios Geológicos por las Islas Canarias. Gran Canaria (F.J. Pérez Torrado y M.C. Cabrera, eds.). Sociedad Geológica de España, Geo-Guías, 5: 41-46. Sánchez, J. (1975): Características y distribución de los suelos en la isla de Gran Canaria. Tesis Doctoral, Universidad de La Laguna, 393 pp. Sánchez, J.; Ríos, C.; Pérez-Chacón, E. y Suárez, C. (1995): Cartografía del Potencial del Medio Natural de Gran Canaria. Memorias y Planos. Cabildo Insular de Gran Canaria, Universitat de Valencia, Universidad de Las Palmas de Gran Canaria. 165 pp+planos. Samper, J. Huguet, L.; García-Vera, M.A. y Ares, J. (1999): VISUAL BALAN: Versión 0 interactiva del programa BALAN para la elaboración de balances hidrológicos. Informe Técnico elaborado para ENRESA. Dpto. Tecnología de la Construcción. Universidad de La Coruña, 97 pp + disquete. Santana, A.I.; Cabrera, M.C.; Cruz-Fuentes, T. y Pérez-Torrado, F.J. (2014): El Ciclo del agua y las aguas subterráneas en la educación secundaria: ideas preconcebidas y propuestas de actuación. En Libro de comunicaciones XVIII Simposio sobre la enseñanza de la geología: 225-235. SPA-15, (1975): Estudio científico de los recursos de agua en las Islas Canarias (SPA/69/515). Ministerio de Obras Públicas. Dirección General de Obras Hídricas. UNESCO, Las Palmas de Gran Canaria-Madrid. MANUALES Manual del Usuario del Programa Visual Balan v. 2.0. Código interactivo para la realización de balances hidrológicos y la estimación de la recarga. ENLACES WEB PHCG (Plan Hidrológico de Gran Canaria) (2010): http://www.aguasgrancanaria.com/ciagcweb/p_noticias2.nsf/E312FE3BC4EC14B580257 687004AECFE/$File/Apartados%201-2.pdf. PHCG (Plan Hidrológico de Gran Canaria) (2015): http://www.aguasgrancanaria.com/. USGS, 2016: http://pubs.usgs.gov/gip/gw/how_a.html. Ültima consulta 14 de junio de 2016.  F. De Paoli La evaluación de la recarga del acuífero del Norte de Gran Canaria  36   37  Experiencia de la realización del Trabajo de Fin de Titulo Las actividades llevadas a cabo durante la realización del TFT han sido principalmente de búsqueda bibliográfica y de datos, análisis estadísticas entre los datos recogidos y estudio y ejecución del programa VISUAL BALAN. La realización del trabajo me permitió mejorar mis habilidades con programas como Excel y ArcGIS y de introducirme en el uso de distintas metodología de balance hídrico. Mi experiencia en el departamento de Física, como en el grupo de investigación GEOVOL ha sido muy positiva, y considero fundamental para mi futuro como persona y como profesional toda la sabiduría que el personal supo comunicarme.