Luis Latorre Albesanace en Beceite el primer día del verano de 1952. A los diez años, junto con tres niñas y dos niños más del pueblo, deja la escuela normal y se matricula de bachillerato, por libre, en el Instituto Nacional de Enseñanza Media “Joaquín Bau” de Tortosa, siéndole impartidas las clases, fuera del horario escolar, por los maestros y el cura del pueblo.
Trabaja durante 17 años en la casa Pelikan, (Pericano no) llegando a ser Jefe de Laboratorio, Investigación y Calidad, donde, aparte de sus funciones, realiza diseños de estuches escolares de plástico, muy de moda en aquel tiempo.
TOSSAL DELS TRES REIS. Muchas veces se apunta a la “tradición” como fuente inefable de que un determinado aserto o suceso es rigurosamente cierto. Tal sucede en la cita que hace Luis Latorre Albesa en su libro sobre Beceite sobre este topónimo de su lugar natal. Y ello, estaría bien como explicación, si la historia no hubiera dado ya respuesta clara y rotunda sobre el tema de losreinos y los condados en la Corona de Aragón.
Sin embargo, siempre hay quien se empeña en continuar caminos de confusión y duda bajo el tópico: “Si el pueblo lo dice y, como el pueblo nunca se equivoca, así será”.
Hay, también, otro topónimo en el Rincón de Ademuz conocido como La Cruz de los Tres Reinos. En este caso nada que objetar. Tres reinos son lindantes; Valencia, Castilla y Aragón. Sin embargo, en el caso de Beceite ni son tres reinos, ni son tres reyes. De sobras es conocido que la formación de la Corona de Aragón se realizó incorporando, repito, incorporando el Condado de Barcelona al Reino de Aragón mediante el llamado Matrimonio en Casa (en la Casa Real de Aragón).
El Reino de Valencia (cristiano) es una creación que corresponde al rey de AragónJaime I el Conquistador y se origina mediante la división de la Potestad Real. (Conquista) Esto es, el reino de Valencia tiene como madre u origen, el Reino de Aragón. Así se ve en los emblemas que usa actualmente la Generalitat Valenciana. El reino de valencia no es otorgado por quien puede crear un reino, esto es, la Iglesia Católica (El Papado), es una división de esa potestad que ya otorgó la Iglesia a los primeros Reyes de Aragón. Particularmente significativo es el viaje de Sancho Ramírez (hijo de Ramiro I) a Roma. En ese Tozal o Tossal confluyeron tres territorios: dos reinos y un condado (reino de Aragón, reino de Valencia y condado de Barcelona). Sí alguna vez subió a ese Tozal / Tossal una potestad regia fueron, en la misma persona, los reyes de Aragón y Valencia. La tercera dignidad representada, menor, también en la misma persona, sería la de conde de Barcelona. Con Felipe V desparecen los reinos de Aragón y Valencia. No desaparece el título de conde de Barcelona que ostentaba el padre del rey Juan Carlos I (don Juan, conde de Barcelona), por ser dignidad menor y no entrar en competencia con la dignidad de REAL.
Las flores también abundan en la comarca del Matarraña.
Se pueden ver amapolas, (ruelles) margaritas e incluso paetonias officinalis, muy difíciles de ver.
En los Puertos de Beceite, gracias a su altitud, aprox. 1500 metros, podemos ver flores propias de estas altitudes.
Algunas se encuentran en peligro de extinción en la zona, por lo que debemos cuidarlas y respetarlas.
Algas en el río Matarraña. En el río Algás o Algars también puede haber, y en el Ulldemó.
Sin duda, destaca en la flora de Beceite y del Matarraña, Teruel, el olivo mandarino, injerto de Tomaset en la finca de los olivares, camino del Pantano de Pena. Esta campaña 2017 será la segunda.
Nuestro colaborador de Paüls (Baix Ebre), que abre antes la boca que los ojos, Jaume Benaiges ha localizado en el archivo parroquial del pueblo evidencias documentales de una tragedia, sin sospechar que formasen parte de una serie de víctimas infantiles atribuidas a carnívoros antropófagos. Unos hechos que sucedieron en el macizo de Els Ports y que al cabo de más de 150 años son difíciles de evaluar. Se hace necesario analizar más casos de Cataluña y sobre todo de Francia, donde resuenan las muertes por unas bêtes que continúan generando debate y atención mediática.
En el Boletín de Historia y Geografía del Bajo Aragón, de los meses de marzo y abril del 1908, encontramos:
- Marín y Vidal, hablando de Peñarroya en su libro La Provincia de Teruel, pág. 147, dice: “El año 1835 apareció en los puertos de Beceite una hiena á la que las gentes bautizaron con el nombre de el lobo blanco que fué por espacio de varios años el terror de todos los pueblos próximos á la referida sierra. Esta fiera penetraba en los cementerios y desenterraba los cadáveres, y llegaba á las cercanías de las poblaciones y acometía á los niños. En Peñarroya hizo cuatro víctimas en distintas ocasiones. La última fué un muchacho de 14 años que murió entre las garras de la fiera el día 3 de Mayo de 1839. Alarmóse la población, salieron todos los cazadores en persecución de la feroz bestia y aquel mismo día fué muerta en la canal de En Pavía por D. Juan Antonio de la Torre, que era una de las personas más distinguidas de Peñarroya y el mejor tirador del país.
José de la Torre, hijo de Juan Antonio, tenía trece años cuando su padre mató a la fiera misteriosa. Más adelante dejaría un escrito donde recordaba qué pueblos fueron los afectados entre aquellos 1835 y 1839: Paüls, Arnes, Miravet, Horta de San Juan, Cretas, Torre del Compte y Peñarroya, finalizando el párrafo con un etcétera que hace plantear más interrogantes.
Menciona tres víctimas infantiles indicando las fechas y los lugares de los ataques, con los nombres y datos de los padres. Informa de una niña y un niño que tenían nueve años de edad y de un tercer niño de catorce. Una niña solo resultó herida.
Cinco ilustres cazadores protagonizarían una batida para intentar hacerse con el animal, batida que fue desde el Cap de la Vinyassa hasta el río de la Canal. El padre del cronista fue quien lo vio primero. Al salirle a cincuenta y cuatro pasos, pudo dispararle la bala que le atravesó el cuerpo.
Lo más curioso del escrito es cuando afirma que posteriormente se descubrió que aquella bestia era una hiena que se había escapado de «unos desbravadores» o domadores de circo.
Abatiendo la fiera exótica de un hipotético circo, se debería haber acabado con el problema, aunque ahora se podrían cuestionar muchos de los detalles si tuviésemos más informaciones. Pero el documento hallado recientemente en Pauls nos anima a dedicar nuestra atención a estas seis víctimas. Más aún porque en el obituario el entonces rector Ramon Valls atribuye las muertes al lobo. Algunos curas omitían las causas de defunción, excepto cuando alguien no recibía los sacramentos. Contrariamente, aquí se consignan casi todas, a pesar de que el mosén sea lacónico cuando lo escribe.
Los historiadores dicen que el origen de los registros parroquiales viene del año 1348 a raíz del alud de víctimas de la peste negra. Las ceremonias religiosas se multiplicaban y había que controlar los cobros. Aun cuando no hubiese intención de que las generaciones futuras pudiesen guardar memoria de algunos hechos, muchos sacerdotes debían intuir esta conveniencia y gracias a ellos encontramos detalles valiosos para investigar. Veamos, si no, los datos de Pauls mencionados en el encabezamiento.
El 11 de junio de 1838, en el camino que va a San Roque, arriba del Freginal de la Retora, se encuentran «las reliquias ó restos del cadáver que ha degollado y en parte devorado un lobo a las diez de la mañana de hoy…». Pertenecen a la niña de siete años Juliana Cortilla Fandos. El 7 de agosto de 1838, en la Solana d’Emperia [sic] localizan los restos «que no se comieron los lobos […] y se dijo ser el de Cosme Tarrega Celma», un zagal de trece años. Diez días después, el 17 de agosto de 1838, se entierra el cadáver de Josep Sanz Lluís. Se había encontrado muerto en los prados «de más abajo del mas de Agustí». La defunción podría ser por el mismo motivo que las otras, pero el cura no lo anota. ¿O quizá es que no ve clara la causa? La víctima parece que tiene catorce años. El 24 de agosto de 1838 dice que Tomasa Laloma Martínez «fue degollada por un lobo en la heredad de sus padres llamada la Coma y cuasi a la vista de los mismos». Tenía cartorce años. Mosén Valls esta vez debe considerar que una muerte tan dramática requiere que se anote con más precisión. El 13 de septiembre de 1838 una nueva víctima de diez años, Francesca Sanmartí Gracià, «fue degollada por un lobo en el barranco de les Ferreres».
En los meses siguientes, se puede suponer que por el éxito de la batida de Peñarroya, no habrá otros casos. En cambio, en Pauls faltaba un último desenlace fatal. El 28 de mayo de 1839, ocho días después de haber matado la supuesta hiena del Matarraña, entierran a Benet Gabriel Viña, de siete años, por «mordedura de un lobo». No se indica dónde pasó, si bien se deduce que estaría herido gravemente de antes, y es el único caso con «Misa de Cuerpo presente» y que «recibió el Santo Sacramento de la Extremaunción». La Extremaunción entonces era la unción in extremis, administrada cuando la muerte se veía venir. A los otros los habían depositado directamente en el cementerio. Se trataba solo de restos y a menudo se lee «no recibió Sacramento alguno».
En todos los casos del pueblo, los escritos no hacen concesiones; no se llama «lobo blanco» a la fiera, ni «hiena», ni «bestia diabólica», a pesar de que podríamos intuir que en aquellos territorios se sintiese pánico. En circunstancias así, únicamente se entendería que si los niños continuaban arriesgándose por los alrededores, era porque debían considerarse imprescindibles en las labores del campo o haciendo de pastores. Constatación que ahora veríamos como de irresponsabilidad de los padres, a menos que entonces una desgracia de estas aún se asumiese como inevitable. Actualmente también sabemos que circular en coche puede tener riesgos y no nos privamos de ello.
Para poder ocupar espacios productivos amenazados por los depredadores, la lucha de nuestros ancestros, que aún depredaban más, tuvo que ser cruel y dolorosa durante generaciones, pero los humanos no solo queríamos sobrevivir. Lo arrasamos todo hasta provocar que muchas especies desapareciesen. Una de las afectadas fue el lobo, y en la península Ibérica hace solo casi un siglo que se extinguió de los territorios mediterráneos y de bastante más allá. Sin memoria de las épocas en que esta fiera se consideraba habitual aquí, hemos perdido infinidad de pautas sobre cómo nos debíamos comportar para prevenirnos y defendernos, aunque más a menudo actuasen movidos por la indignación por los daños, o bien se resignasen a la fatalidad y no hiciesen nada. Existía todo un mundo de estrategias o de infortunios, así como de actitudes a la hora de enfrentarse a los depredadores para subsistir o de arriesgarse para disponer de más recursos. Muchas veces
la gente se dejaba la piel porque la intención de enriquecerse les hacía exponerse más.
Ya en siglo xx el rastro de las desgracias se olvidaría, enterrado con los recuerdos de los más viejos, que se asilaban o bien quedaban inmersos en los suburbios industriales de muchas ciudades. El desarrollo desbocado pasaba por alto la consideración ética del entorno, aunque el contacto con la naturaleza y sus valores empezaba a añorarse. Otra visión sustituiría la realidad dura y vivida. Faltaba poco para que se propiciase una cierta nostalgia, quizá bucólica. El lobo corría un peligro inminente de extinción irreversible en buena parte de Europa y los movimientos conservacionistas querían que se diese una nueva oportunidad para cohabitar con la fiera más legendaria del hemisferio norte. Se implantarían acciones de recuperación y la finalidad sería la de salvar Canis lupus. Pese a las treguas humanas, hoy aún no es posible sacarlo de muchas listas de especies amenazadas, pero determinados núcleos de varios países se han recuperado y se expanden.
Para favorecer estos intentos, con buena fe o desconocimiento histórico y situándose en uno de los extremos ideológicos, aparecían los que negaban a ultranza cualquier indicio de víctimas humanas atribuidas al lobo. Datos que también se van encontrando en Cataluña, aunque de manera exigua quizá porque se ha investigado poco.
En esta reflexión tenemos que preguntarnos: ¿Una hiena podría vivir durante cinco años en Els Ports? ¿Sería albino el animal? ¿Existió el circo de que habla la crónica? ¿Por qué la gente, que entonces conocía bien al lobo, lo llamaba «el lobo blanco»? ¿O es que se trataba de un lobo albino singular o de pelaje claro por ser viejo? ¿No sería un mastín cruzado con lobo? ¿O un perro asilvestrado?
La diagnosis no se puede aventurar sin más pruebas. Debemos advertir que la impunidad de que se gozaba en todas partes cargando al lobo con las culpas fue real y documentada, hasta que llegó la criminología moderna. Ahora obtendríamos muestras de pelo, saliva o excrementos de cada depredador para analizar el ADN y dictaminar quiénes son los auténticos responsables.
Sobre estos casos de Els Ports, hemos querido consultar a Jean-Marc Moriceau,profesor de la Universidad de Caen, en la Baja Normandía, y experto en la temática de víctimas humanas atribuidas a las fieras. Sus libros han dado pie al debate al haber hecho públicos muchos registros galos semejantes a los de Pauls. Según él, «se trata de un lobo muy peligroso, de comportamiento inusual».
El profesor Moriceau explica con detalle su postura: «La designación de “hiena” es genérica para calificar cualquier animal salvaje que ataca al hombre y que se distingue por su comportamiento, en lugar de los depredadores del territorio. En la Francia del siglo xviii (e incluso en 1814) en cuatro casos analizados las mentalidades populares evocaban a la hiena para subrayar bien al agresor y aislarlo de la conducta de sus congéneres: en 1756 en el Lionés, en 1764-1770 en Gavaldá, en 1788 en el Esterèu y en 1815-1818 en el Nivernés. En el caso de Cataluña, excepto que se encontrase el origen exacto de la supuesta hiena y la identidad de su propietario (que no parece ser el caso) nos atenemos a los rumores populares que reflejan un choque psicológico por contraste con la conducta habitual de los animales depredadores de la región, a saber, los lobos. En resumen, salvo nuevas pruebas, creo que esta supuesta hiena es un lobo devorador de hombres.»
En Cataluña había habido acciones oficiales, con contenidos significativos porque dan claves de cómo algunos lobos podían acostumbrarse a comer carne humana. En octubre de 1816 el capitán general ordena batidas advirtiendo que hay «personas mordidas por lobos» de los mismos alrededores de Barcelona. Decían que, rabiosos, habían devorado a una chica, matado a dos hombres y herido a otro, cuando si los lobos hubiesen estado avanzadamente rabiosos no habrían podido comer nada. También decretaban cacerías porque «aun quando los lobos no sean rabiosos parece que han contraído el hábito de cebarse en carne humana […] tal vez por la proporción que han tenido en esta última guerra».
Por lo que respecta a episodios igualmente documentados, hay uno de El Vallès Oriental que daba a conocer el médico Bonaventura Sendil en 1825:
En Vilamajor cuatro muertos entre niños y niñas. Cogidos tres, dos niños y una niña, pero salvados, aunque gravemente heridos. En Monseny un niño muerto y devorado. En Santa Susana uno muerto y devorado. En Santa Susana uno muerto y otro gravemente herido, pero defendido. En San Esteban de Palau Tordera, uno muerto y la mayor parte devorado, y tres gravemente heridos pero defendidos. En Palau Tordera una niña muerta, y devorada la mayor parte. En Campins una niña muerta, y un niño cogido, pero defendido sin ser herido. En Collsabadell una niña muerta, y una cogida y gravemente herida. En Vilalba Saserra una niña cogida y defendida sin ser herida.
Acaba diciendo que el último ataque se produjo el día 18 de marzo contra una niña, «devorada enteramente, a excepción de un pedazo de brazo y un poquito de cabeza. Advierto todos estos daños y puedo dar justificación siempre que sea necesario».
Las personas no somos presas habituales del lobo, que prefiere herbívoros. Pero tampoco se pueden descartar sus capacidades como depredador oportunista, inteligente y carroñero saneador del medio. De primavera a verano puede arriesgarse cuando no encuentra presas para criar a los lobeznos. No podemos pasar por alto que los meses con más personas afectadas coinciden mucho con los de los estudios estadísticos de Moriceau. Es cuando había más presencia humana en los campos, los sembrados podían servir de escondrijo y permitir a los lobos observar a las posibles presas vulnerables. Sin el ganado, que durante medio año solía alejarse por la trashumancia, era difícil la supervivencia de las camadas que campaban cerca de la gente. Más aún si habían criado y contando con el recurso del ganado, enfermo o sano y de la carroña que generaba.
No es extraño, por tanto, que los depredadores se sintiesen atraídos por los pueblos donde había otro tipo de ganado, desperdicios y olores de alimentos o que se decidiesen finalmente por los humanos, dada la desorganización social y los muertos debidos a las guerras. Quedan archivos de poblaciones aquí mencionadas donde hacer investigación y obtener nuevos datos. Serían útiles para conocer más detalles de interés sobre estos hechos que ahora nos parecen insólitos. Bibliografía Bevanger, K. y L. L. Olden, 2002. La peur du loup (en chapurriau, la po del llop). Recueil d’attaques de loups sur des humains. Norsk Institutt for Naturfonskning. Trondheim, Noruega. http://www2.nina.no/lcie_new/pdf/634993231241795116_Linnell%20NINA%20OP%20La%20Peur%20du%20loup.pdf Burnier, F., 2001. «Comment peut on vivre avec des loups sauvages?». En chapurriau : Com se pot viure en los llops salvaches ? http://controverses.mines-paristech.fr/public/promo03/C03B18/liens/Commentpeut-onvivreavecdesloupssauvages.2001-11-116.PDF Bulletin des medecins suisses, 11: 549-554.
Costa i Savoia, E., 2012. Paüls. Arola Editors. Tarragona. Curto i Homedes, A. y A. Martínez Silvestre, 1999. «La presència del llop a l’antic terme de Tortosa durant la baixa edat mitjana».
ALBERT CURTO HOMEDES, ALBERT MARTÍNEZ SILVESTRE PDF aquí https://www.raco.cat/index.php/ActaHistorica/article/view/188798/254596 Acta Historica et Archaeologia Medieaevalia, 20-21: 455-475.
Manent, A., 1995.
El llop i els humans. Passat i present a Catalunya. Arola. Tarragona.
Massip, J. M. y J. Maroto (eds.), 2011. Llops i humans a Catalunya. Del passat al present. Centre d’Estudis Comarcals de Banyoles. Banyoles.
Moriceau, J. M., 2007. Histoire du mechant loup. 3000 attaques sur l’homme en France. xve-xxe. siècle. Fayard. París.
Moriceau, J. M., 2011.
Para la correcta estimación de la recarga por precipitación a los acuíferos tiene especial relevancia tener en cuenta la generación de escorrentía superficial. No considerarla en el cálculo de los recursos hídricos subterráneos puede suponer una sobreestimación de los mismos. En el sistema acuífero del Baix Ebre, en el Sur de Cataluña, es preciso evaluar la escorrentía superficial y de la zona vadosa que se producen en las formaciones carbonatadas karstificadas del macizo de los Ports de Beceit con el objetivo de realizar una estimación más aproximada de los recursos transferidos desde ese macizo a la Plana de La Galera. A partir del modelo conceptual hidrogeológicose realiza la estimación de la escorrentía superficial media anual, que incluye la de los acuíferos colgados temporales de los Ports de Beseit, tanto en su vertiente NW, cuenca del río Matarraña, como en su vertiente SE, Plana de La Galera. Para ello se analizan los hidrogramas de caudales de los ríos y sus afluentes, los hidrogramas de llenado y vaciado del embalse de Ulldecona y el balance de agua en el suelo realizado con el código Visual Balan en la cuenca receptora del embalse de Ulldecona.
La escorrentía superficial en los Ports se ha estimado en 105±20 mm·a−1, lo cual supone el 20–30% de la recarga media anual calculada por los métodos de balance de agua en el suelo y deposición atmosférica de ión cloruro, en torno a 350–500 mm·a−1, cuya mayor parte es transferida lateralmente a la Plana de La Galera.
Escorrentía superficial; recursos hídricos; karst; Baix Ebre; Ports de Beseit; embalse de Ulldecona; análisis de hidrogramas.
recibido el 8 de mayo de 2015 / Aceptado el 23 de octubre de 2015 / Publicado online el 18 de abril de 2016
Citation / Cómo citar este artículo: Espinosa Martínez, S. & Custodio, E. (2016). Estimación de la escorrentía superficial para el cálculo de la recarga a los acuíferos del macizo kárstico de los Ports de Beseit (Tarragona, España) combinando balance de agua en el suelo y análisis de hidrogramas de caudales. Estudios Geológicos 72(1): e045. http://dx.doi.org/10.3989/egeol.42132.374.
En cuencas de clima árido-semiárido y sub-húmedo mediterráneas, la disminución de la incertidumbre en la estimación de la recarga debido a la escasez de recursos hídricos que suele caracterizarlas tiene especial importancia. En dichas áreas, la demanda de caudal para el abastecimiento poblacional y para el desarrollo de actividades agrícolas es notable y por ello un buen conocimiento de sus recursos hídricos, en especial de los subterráneos, juega un papel clave en su desarrollo socio-económico (Sahuquillo et al., 2009; Custodio, 2011; Molinero et al., 2011).
En estas áreas las actividades humanas se desarrollan frecuentemente en zonas llanas y los recursos de agua disponibles proceden de acuíferos recargados desde sistemas montañosos adyacentes. Su estudio ha tomado especial relevancia en las últimas dos décadas en áreas áridas y semiáridas de diversos lugares del mundo, incluyendo el oeste Norteamericano (Manning, 2002; Nelson & Mayo, 2014; Welch & Allen, 2014), Nuevo México (Wilson & Guan, 2004), Sudamérica (Varni & Custodio, 2013), sur de Australia (Guanet al., 2010), sur peninsular español (Alcaláet al., 2006; Alcaláet al. 2011), área mediterránea continental de la Cordillera Prelitoral Catalana (Lambán, 1998), sistemas insulares españoles (Cardoso, 1997; Cruzet al., 2011; Cabreraet al., 2013) y el Baix Ebre en el noreste peninsular español aquí considerado (Espinosa, 2014; Espinosaet al., 2015). Un aspecto importante que afecta a la validez de las evaluaciones de la recarga en áreas montañosas es la estimación de la escorrentía superficial y subsuperficial media anual, incluyendo la producida en el medio no saturado (zona vadosa) debido a la formación de acuíferos colgados de funcionamiento temporal (Beven, 1989; Hardie et al., 2012; O’Brien et al., 2013; Custodio et al., 1997).
La importancia del estudio de la escorrentía superficial ya ha sido establecida en macizos montañosos de clima atlántico europeo en los que la escorrentía supone del 35 al 55% de la precipitación media anual (Samper et al., 2015) y en sistemas insulares áridos (Cruz-Fuentes, 2008). Estas estimaciones parten de un modelo conceptual de funcionamiento hidrogeológico, basado en los datos y observaciones disponibles y en lo posible validado mediante estudios y consideraciones hidrogeoquímicas e isotópicas ambientales (Allison, 1988; Gaspariniet al., 1990; Allisonet al., 1994; Wood & Sanford, 1995; Horneroet al., 2013).
Este trabajo tiene como objetivo la estimación de la escorrentía superficial media anual para contribuir a acotar la incertidumbre en la estimación de los recursos hídricos subterráneos del sistema kárstico de los Ports de Beseit (en adelante referido también como Els Ports).
Metodología
Para la estimación de la escorrentía media anual en los Ports de Beseit se ha procedido a la recuperación de las series históricas de caudales diarios de las bases de datos de la Confederación Hidrográfica del Júcar (CHJ) en lo referente al embalse de Ulldecona (vertiente sureste de los Ports) y de la Confederación Hidrográfica Ebro (CHE) en cuanto a los cauces afluentes al río Matarraña con nacimiento en los Ports (ríos Figuerales, Pena, Ulldemó) y del río Algars, como apoyo a la interpretación del comportamiento hídrico del macizo de los Ports de Beseit (Figura 2), por estar enclavados en el macizo y presentar características geomorfológicas en cuanto a pendientes, escarpes, superficies de afloramiento de roca, karstificación y vegetación que se pueden considerar como representativas del conjunto territorial.
A partir de los datos de caudales diarios recopilados de trabajos previos (Tourís, 1986; CHE, 1991; CHE, 2009; Badiella, 2009) y de los registros temporales de las estaciones de aforo, se realizó un estudio preliminar general de las condiciones hidrométricas de contorno y de los límites de la zona de estudio (Espinosa, 2014; Espinosa et al., 2014a y b) para determinar la idoneidad de dichos datos para la aplicación de los métodos de estimación de la escorrentía superficial.
Se dispone de las series pluviométricas de las estaciones meteorológicas de Fredes y del Parque Natural (PN) Els Ports, esta última la más representativa del área considerada. Como la estación PN Els Ports entró en funcionamiento en enero de 2007, se ha tenido que extender la serie de datos a partir de los de la estación de Mas de Barberans, que está óptimamente ubicada en el centro de la cuenca de la Plana de La Galera (Figura 4), realizado una corrección por altitud (Espinosa, 2014). Con los hidrogramas de aforo para el periodo 1994–2009 se ha realizado un estudio de la distribución de años húmedos y secos, en función de la desviación de la pluviometría acumulada anual respecto a la pluviometría media anual, existiendo coherencia entre los datos registrados por las estaciones de aforo y los eventos de lluvia registrados por las estaciones meteorológicas. En la serie se aprecia a grandes rasgos una secuencia de periodos lustrales de años secos y años húmedos. Las investigaciones se han desarrollado en una época de tendencia seca.
A partir de la representación gráfica de las series históricas de caudales, tanto del embalse de Ulldecona como de los afluentes de los ríos Matarraña y Algars, se ha determinado de forma aproximada la escorrentía media anual, considerando las superficies de las cuencas receptoras correspondientes.
Para la interpretación de los hidrogramas se definen tres horizontes: a) suelo superior y roca aflorante, con escasa cubierta vegetal y donde tiene lugar la evapotranspiración y lo que se puede denominar escorrentía superficial y subsuperficicial o rápida; b) tramo de gran espesor no saturado o zona vadosa en cuyo seno se pueden formar acuíferos colgados semipermanentes y ocasionales sobre intercalaciones litológicas de muy baja permeabilidad y cuyo agotamiento produce una escorrentía lenta que se suma a la escorrentía superficial; y c) zona saturada donde tiene lugar el flujo subterráneo. La designación de suelo superior debe considerarse como el suelo que pueda existir, en general pobre y de escasa profundidad, más la roca alterada y fracturada somera, que en algunos lugares puede llegar a ser un exokarst.
Paralelamente se ha realizado el estudio del comportamiento hidrológico de la cuenca receptora del embalse de Ulldecona, (ojo de corteza en castellano) aplicando técnicas hidroquímicase isotópicas ambientales, que han sido básicas para formular el modelo conceptual de funcionamiento hidrológico que soporta la interpretación del balance de agua en el suelo con el código Visual Balan en dicha cuenca. Para el balance de agua en el suelo se ha partido de los parámetros hidráulicos del suelo superior, zona no saturada y acuífero utilizados en la resolución del balance de agua en el suelo general (Espinosa, 2014). El balance de agua para la cuencareceptora del embalse de Ulldecona se ha calibrado con los datos de aforo del embalse de Ulldecona, que consisten en los valores diarios de la reserva de agua en el embalse y los caudales de salida, en el periodo de octubre de 1994 a septiembre de 2009.
Área de estudio
El área de estudio se encuentra localizada en el Baix Ebre, noreste de la Península Ibérica, entre las provincias de Tarragona, Castellón y Teruel. Els Ports actúa de límite provincial y de divisoria de aguas superficiales (Figura 1). Se trata de un gran macizo montañoso, por lo que es importante trabajar a escala regional para poder considerar su contribución a la recarga de los acuíferos del llano adyacente de la Plana de La Galera.
Figura 1.—Localización de la zona de estudio, Baix Ebre, S de Cataluña, NE de la Península Ibérica (izquierda) y Embalse de Ulldecona en el sistema Ports de Beseit, con los límites de la zona de estudio (derecha).
Se ha partido del modelo hidrológico conceptual general de la zona definido en la Tesis Doctoral del primer autor (Espinosa, 2014), habiendo sido necesario avanzar previamente en el conocimiento del funcionamiento hidrodinámico local a partir de la aplicación de técnicas hidroquímicas e isotópicas ambientales en el entorno del embalse de Ulldecona.
En la vertiente SE de los Ports de Beseit, parte de la escorrentía superficial producida en Els Ports es transferida a través de los materiales de piedemonte hacia la Plana de La Galera, infiltrándose en el acuífero a medida que avanza en esa dirección, habiéndose estimado su valor entre 95 y 115 mm·a−1 con el balance general de agua en el suelo, tal y como se describe en Espinosa (2014). Estos resultados proporcionan una idea de la magnitud de la escorrentía superficial, pero carecen de la robustez necesaria al ser muy dependientes de los parámetros utilizados en los balances, aunque se han intentado calibrados a partir de los niveles piezométricos. Además, no han podido realizarse mediciones in situ del caudal de escorrentía superficial debido a que su existencia en el Baix Ebre queda limitada a los eventos de precipitación extraordinarios, los cuales son ocasionales, con periodicidad groseramente anual y de duración menor a una semana. Ha resultado imposible coordinar logísticamente alguna posible campaña de aforo y muestreo durante alguno de estos eventos.
Tres factores principales hicieron plantear la hipótesis de la importancia de la escorrentía superficial en el sistema Ports de Beseit-Plana de La Galera: (1) ocurrencia de eventos de precipitación estacionales de fuerte intensidad, con escurrimiento superficial, los cuales se observaron en campo y son bien conocidos por la población local, (2) existencia de un embalse en la zona sur de los Ports de Beseit con un comportamiento de llenado y vaciado ligado a estos eventos de precipitación y (3) los resultados previos obtenidos en el balance de agua en el suelo con el código Visual Balan, que muestran que la escorrentía superficial media anual en la zona de montaña, puede suponer aproximadamente un 20–25% de la recarga media anual (Espinosa, 2014; Espinosa et al., 2015).
El área estudiada se caracteriza por un fuerte contraste de relieves, entre 2 y 300 m en la Plana y hasta cerca de 1500 m de los Ports. La precipitación media es de 550 mm·a−1 en la Plana de la Galera (condiciones de semi-aridez), aproximadamente 1100 mm·a−1 en la zona alta de los Ports de Beseit (clima mediterráneo sub-húmedo) y de 690 mm·a−1 en la cuenca del río Matarraña. Se trata de lluvias estacionales que pueden llegar a ser intensas. Se pueden diferenciar 4 grupos principales de usos de suelo y cobertera vegetal: coníferas y matorral bajo en la zona de los Ports y cultivo de secano y mixto (con regadío de apoyo) en la Plana y cuenca del río Matarraña.
Contexto geológico regional
El área pertenece a la rama NE de la Cordillera Ibérica en su tránsito a la Cordillera Prelitoral Catalana. Se pueden distinguir tres zonas: los Ports de Beseit, que son la cabecera de la red hidrográfica, con sus dos vertientes, la cuenca del Matarraña limitando con la depresión del Ebro al noroeste de los Ports y la Plana de La Galera al sudeste de los Ports.
En los Ports de Beseit predominan los materiales mesozoicos, marcando un tránsito entre los dominios de la Cordillera Ibérica y las Cordilleras Costero-Catalanas. Estos materiales se apilan en escamas cabalgantes y pliegues con direcciones dominantes SW-NE en la zona oriental y dirección NW-SE en la zona occidental y vergencia dominante al NW, donde los niveles arcillosos del Muschekalk medio y del Keuper actúan como planos de despegue y su existencia da lugar a manantiales que generan el caudal de base en la cabecera de los ríos Matarraña y Algars. La estructura geológica facilita la presencia repetitiva de afloramientos permeables de edad mesozoica, en ocasiones con alto grado de karstificación evidenciado por la existencia de lapiaces y en especial de cuevas y simas, las que están más desarrolladas en la cara NW del macizo. Se dispone de una cartografía de los “avencs” (galerías) de la potente zona no saturada explorados por el Club de Espeleología de Tortosa y realizada por ellos mismos (CHE, 2008).
Los materiales dominantes, notablemente fallados y fracturados, son las calizas y dolomías mesozoicas, además de conglomerados terciarios asociados al frente de cabalgamiento. Las calizas y dolomías mesozoicas también afloran formando las sierras de Godall y Montsià al SE de la zona de estudio (Figura 2).
Figura 2.—Mapa geológico de la zona de estudio según la Confederación Hidrográfica del Ebro (http://hidrologia.geoslab.com/HydroGeoEbro/). Corte geológico A–B (Beseit-Mar Mediterráneo), según Bayó et al., (1992).
Los acuíferos presentes en los Ports de Beseit corresponden a una serie fundamentalmente dolomítico-calcárea que abarca el Triásico, Jurásico y Cretácico (Mesozoico). Aunque existen niveles margosos de baja permeabilidad intercalados, especialmente en las series del Triásico y del Dogger (Jurásico), la compleja fracturación tectónica pone en contacto entre sí los diferentes niveles permeables, de forma que cabe considerar al conjunto como un sólo acuífero kárstico de elevada heterogeneidad, aunque con la posible existencia de acuíferos colgados de funcionamiento temporal. La karstificacióntiene un desarrollo muy diferente en función de su localización en la serie estratigráfica. Los materiales mesozoicos karstificados que afloran en los Ports continúan en profundidad hacia la Plana de La Galera, lo que da lugar a la posibilidad de un flujo de agua subterránea perpendicular a las estructuras desde los Ports a la Plana de La Galera.
Los materiales más abundantes que afloran en la depresión del Ebro son de edad cenozoica (Terciario y Cuaternario). Son materiales terrígenos (conglomerados, areniscas y margas), calizas lacustres y evaporitas (yesos). Los cauces principales (ríos Matarraña y Algars) se encajan en su tramo medio-bajo en depósitos aluviales que afloran en el lecho del río a lo largo de todo el recorrido (CHE, 2008).
La Plana de la Galera forma parte de una fosa tectónica (graben) delimitada por un conjunto de fallas subverticales que separan la fosa de la Plana de La Galera, rellena de materiales de edad terciaria a cuaternaria con predominancia de materiales plio-cuaternarios, de los bloques levantados (al SE horst del Montsià y Godall y al NE-SW de los Ports de Beseit) donde afloran materiales de edad mesozoica. El río Ebro circula por una falla de desgarre que desplaza su margen izquierdo en dirección NW respecto a su margen derecho; esta falla actúa de límite geológico e hidrogeológico (Tourís, 1986; Badiella, 2009).
Condiciones hidrológicas en los contornos
Los principales cauces presentes en la vertiente NW de los Ports de Beseit, ordenados de SW a NE, son los ríos Tastavins, Pena, Matarraña y Ulldemó, que constituyen la cabecera del río Matarraña, y los ríos Algars y Canaleta. Sus cursos están adaptados a los cambios orográficos, que a su vez están condicionados por la geología y la tectónica, con barrancos de mayor pendiente en cabecera y que se suavizan hacia la depresión del Ebro. La orientación de los cauces es de S a N en el sector occidental y cambia a SE-NW en el sector oriental.
En la vertiente SE, la red fluvial de los Ports de Beseit está compuesta por una serie de barrancos que recogen tanto la descarga de las surgencias de los niveles acuíferos colgados como la escorrentía superficial producida en los Ports, desembocando en el margen derecho del río Ebro. Los barrancos de mayor entidad son los barrancos de Xalamera, Paüls, Conca, Lloret, Sant Antoni y La Galera (Figura 3).
Figura 3.—Ubicación de subcuencas de la vertiente Matarraña y embalse de Ulldecona, así como de los puntos de aforo y estaciones meteorológicas utilizadas. Se resaltan los límites de las cuencas principales y el área semiendorreica notablemente karstificada de la Mola de Catí.
Existen dos embalses construidos en la década de 1960, el embalse de Pena en la vertiente W y el embalse de Ulldecona en la vertiente E, de 18 hm3 y 11 hm3 de capacidad respectivamente. En ambos casos el agua embalsada se destina al abastecimiento local y riego.
Resultados
Caracterización de la cuenca receptora del embalse de Ulldecona para la definición del modelo conceptual de funcionamiento
El embalse de Ulldecona, emplazado en el sureste de los Ports de Beseit (Figura 3), capta la cabecera del río de La Sénia y regula su curso aguas arriba de la pedanía de El Castell. Sus reservas se destinan al regadío en La Sènia y Ulldecona (ver Figuras 1 y 3).
La cuenca receptora, de 126 km2, consiste en dos subcuencas distintas en su ubicación y orientación: el barranco de La Fou hacia el NNE y el barranco de Corruscars hacia el SW (Figura 3) y desarrollo predominantemente según la de los Ports. La subcuenca del barranco de la Fou tiene mayor longitud y profundidad que la del barranco de Corruscars. Además, en esta última subcuenca la densidad de pinar es mayor que en la del barranco de La Fou. Estas características geomorfológicas parecen reflejarse en la química elemental del agua del embalse muestreada en el barranco de Corruscars (estación JUA802 en la Figura 3), no mostrada aquí, a través de una mayor concentración en iones mayoritarios respecto a la química del agua del embalse muestreada a la salida del mismo (estación JUA801 en la Figura 3), procedente en su mayoría del barranco de La Fou.
El embalse puede llenarse hasta verter por coronación con eventos de lluvia de dos o tres días de duración que totalicen unos 100 mm. Cuando esto ocurre, el embalse mantiene el nivel máximo de la lámina de agua durante aproximadamente 1 a 1,5 meses, con una evacuación diaria en el frente de presa de unos 0,5 m3·s−1 de agua. Por tanto, el llenado del embalse se produce tanto por la escorrentía superficial rápida como por la que se deriva de acuíferos colgados temporales y de cierta descarga subterránea. Para comprobarlo se procedió a la caracterización hidroquímica del agua embalsada en comparación con la composición química del agua de lluvia (punto de control TMM P3-B, Figura 4) y con la de agua subterránea de surgencias y pozos de los alrededores del embalse (Figura 4), no mostrada aquí.
Figura 4.—Localización de los puntos de control de calidad del agua del embalse de Ulldecona de la red de puntos de aforo de la Confederación Hidrográfica del Júcar (Mapa base: Google Earth-Institut Cartogràfic de Catalunya).
Se trata de aguas bicarbonatadas cálcico-magnésicas con concentración de sulfato variable. Se pueden diferenciar 3 grupos a partir de los cuales se plantea la hipótesis del modelo conceptual de funcionamiento: a) grupo de menor concentración que corresponde la lámina de agua más superficial del embalse y las surgencias Font del Teix y Font de La Fou que descargan al embalse; su concentración es comparable a la del agua de lluvia, lo que indica que la descarga ha sido rápida, con escasa interacción con el terreno b) grupo con concentración media que corresponde a la Font dels Rossegadors y a las muestras del frente de presa localizadas al pie de la misma; indica que el agua embalsada procede tanto en la escorrentía superficial como del agotamiento de los acuíferos colgados de la zona vadosa como de la propia lluvia y c) grupo con mayor concentración de iones en general y en especial de sulfatos que corresponde a las muestras de agua de la Font de Sant Pere y piezómetro surgente aguas abajo del embalse, que es surgente la mayor parte del año (Figura 4); mayor tiempo de circulación por el terreno.
Análisis de hidrogramas
Estudio de las series temporales hidrométricas diarias del embalse de Ulldecona
Se dispone de las series temporales diarias de las reservas y de las salidas del embalse de Ulldecona para un periodo de 15 años. Por tanto es posible calcular las entradas diarias al embalse a partir de:
E=ΔV+S+Ev+f
siendo E=entradas, ΔV=variación de volumen, S=salidas, Ev=evaporación, f=fugas del embalse.
La serie de entradas al embalse (E) ha sido considerada como la escorrentía total de la cuenca, que es la suma de la escorrentía superficial y la de acuíferos colgados temporales. Las variaciones de volumen se han calculado a partir de la serie diaria de reservas calculadas en el embalse. Se supone que las fugas del embalse por la presa y su entorno son nulas ya que en la década de 1970 se realizaron inyecciones de cemento de la cerrada y no se produce drenaje significativo por sus galerías.
La evaporación media anual de agua del embalse (superficie máxima del embalse de 0,47 km2) calculada para el periodo 1994–2009 (Espinosa, 2014) representa el 5% de las entradas anuales, por lo que no se considerará en el cálculo de la serie temporal diaria de entradas ya que su valor es similar a los errores de medida y estimación del resto de variables.
La Figura 5 muestra la relación lineal ajustada de la serie diaria de entradas (escorrentía, E), con la precipitación anual (P) en la estación meteorológica PN Els Ports. Esta relación es E=0,23 (P-728) para valores anuales en mm. El umbral de precipitación (728 mm·a−1) por debajo del cual no se produce escorrentía tiene asociada una notable incertidumbre. Esta incertidumbre es atribuible a que la ocurrencia de escorrentía a lo largo del año depende de las condiciones previas de la cuenca, tanto meteorológicas como del estado de humedad del suelo. A escala de detalle, según observaciones en campo, los eventos concentrados de precipitación superiores a 100 mm en un día suelen dar lugar a escorrentía superficial, pero también en función de las condiciones previas de humedad en que se encuentre el terreno.
Figura 5.—Relación entre entradas anuales (mm·a−1) del embalse de Ulldecona y precipitación anual (mm·a−1) de la estación PN Els Ports (izquierda) y de los valores acumulados (mm) (derecha).
En el hidrograma se diferencian 2 fases en el agotamiento de los eventos de caudal de entrada (fase 1 y 2 en la Figura 6), una con pendiente mayor que la otra. Los tiempos de semi-agotamiento (tiempo para que el caudal se reduzca a la mitad) se han calculado en unos 6 días en promedio para la fase 1 y en unos 35 días para la fase 2. Se plantea la hipótesis de que la fase 1 representa la escorrentía superficial y la subsuperficial que se produce en la parte superior del terreno y la fase 2 la escorrentía procedente del agotamiento de los acuíferos colgados temporales.
Figura 6.—Cálculo del agotamiento en los eventos de caudales de entrada del embalse de Ulldecona. Figura superior: hidrograma de caudales (no se han representado los caudales mayores que 14 m3·s−1) y selección de los tramos de ampliación; los valores de <0,2 m3·s−1 no tienen significación porque están dentro del error de estimación. Figura inferior: detalle de los agotamientos de caudal de entrada al embalse de Ulldecona con ordenadas en escala logarítmica.
De la serie diaria de caudales se ha calculado la aportación total, que es de 103 mm·a−1 para las entradas y de 100 mm·a−1 para salidas.
A partir de la representación gráfica de las entradas diarias al embalse las aportaciones medias anuales de escorrentía superficial del suelo superior se han evaluado en 30 mm·a−1 y en 73 mm·a−1 para la zona vadosa, como el área bajo la curva descrita por el hidrograma (Figura 7). El volumen de agua de escorrentía de la zona vadosa se calcula como la suma de las áreas delimitadas por el dibujo entre los agotamientos intermedios de cada evento de entrada (fase 2) y el volumen de agua de escorrentía del suelo superior (fase 1).
Figura 7.—Entradas al embalse de Ulldecona en el periodo octubre 1994-septiembre 2009 e interpretación del hidrograma. La parte sombreada en azul es lo se que supone que corresponde a la escorrentía de la zona vadosa, apreciado por ajuste manual (Qe = caudal de entrada en hm3·d−1).
Análisis de hidrogramas de los afluentes del río Matarraña y Algars
Como apoyo a la interpretación del comportamiento hídrico del macizo de los Ports de Beseit frente a eventos de lluvia se han estudiado los hidrogramas de las series temporales de caudal para los ríos con nacimiento en los Ports localizados en la vertiente noroeste del macizo (Figura 3).
A partir de la representación gráfica y tras separar la escorrentía del suelo superior de la de la zona vadosa (sombreados de color rojo y azul respectivamente en la Figura 8), se han obtenido los valores de aportación total, aportación por escorrentía superficial y aportación por escorrentía de la zona vadosa (Tabla 1), considerando las superficies de las respectivas cuencas aguas arriba del punto de aforo correspondiente (Figura 3).
Figura 8.—Caudal del río Ulldemó en el periodo octubre 1994-septiembre 2009 e interpretación del hidrograma. El sombreado en rojo es lo que se identifica como escorrentía de los acuíferos colgados temporales en la zona vadosa y el en azul es la escorrentía del medio subterráneo.
Tabla 1.—Aportaciones medias anuales para los afluentes del río Matarraña y para el río Algars. La situación espacial de la estaciones se muestra en la Figura 3
PERIODO
ID
ESTACIÓN
Años
Sup. km2
Aportaciones medias anuales (mm·a−1)
AT
AS
ASS
AL
1994–2010
9113
Figuerales
16
55
2,25
0,25
2,00
1994–2010
9110
Pena (Beceite)
16
49
52,49
4,25
48,24
1994–2007
9052
Matarraña (Ulldemó)
13
48
176,63
15,46
81,7
80,00
1994–2010
9153
Algars (Horta S. Joan)
16
115
180,4
19,20
160,94
*A: aportación (T: total, S: escorrentía superficial del terreno, SS: escorrentía zona vadosa, L: flujo subterráneo, ID: identificación).
El río Ulldemó es el que reproduce mejor los eventos de precipitación con generación de escorrentía similar a la producida en el embalse de Ulldecona. Los restantes difieren ya que el río Figuerales es de pequeña entidad y no nace propiamente en el macizo, el río Pena se alimenta de manantiales permanentes al pie de los Ports y el río Algars nace en el área notablemente karstificada de la Mola del Catí (Figura 3) y además recibe a lo largo de su recorrido aportes permanentes que incrementan considerablemente su caudal.
Balance de agua en el suelo con el código Visual Balan en la cuenca receptora del embalse de Ulldecona
La obtención de la escorrentía media anual en la cuenca receptora del embalse de Ulldecona se ha realizado mediante el balance de agua en el suelo con el código Visual Balan v2.0 (Samper, 1997/2005). Se ha partido del balance de agua en el suelo superior de la zona de estudio (Espinosa, 2014), utilizando el método del Número de Curva (US-SCS, 1972; 1985; Ferrer et al., 1995) para calcular la escorrentía superficial de los eventos de lluvia y se han utilizado los mismos parámetros de la cubierta vegetal y suelo superior, dado que las características son similares. Dichos parámetros representan un suelo poco desarrollado con roca aflorante y pendientes mayores al 1% (Tabla 2, Espinosa, 2014). Además, se ha considerado que existe interceptación puesto que hay zonas con pinares relativamente densos.
Tabla 2.—Parámetros hidráulicos finales introducidos en el código Visual Balan para el balance de agua en el suelo en la cuenca del Embalse de Ulldecona
DATOS
Años
15
HIDROMETEO
EMET
PN Els Ports
Punto Control
Aforo salida Embalse de Ulldecona
área cuenca receptora (km2)
126
CUBIERTA VEGETAL
Interceptación-Cubierta vegetal
Bosque de pinares
Método de cálculo: Singh
Capacidad de almacenamiento (mm)
1,27
Coeficiente de interceptación
0,2
SUELO EDÁFICO
Porosidad total
0,3
Humedad inicial
0,04
Conductividad hidráulica (m·d−1)
0.1
Reserva útil máx. (mm)
20
ETP
datos usuario
ETR por el método exponencial
Umbral déficit hídrico, CEME (mm)
10
CKRD
0,1
Escorrentía superficial (Número de curva)
Lluvia mínima del aguacero (mm)
1,5
Número de curva
50
ZONA NO SATURADA
Resolución por Método Explícito
Coef. agotam. flujo hipodérmico (d−1)
0.099
Tiempo semi-agotamiento (d)
7
Contenido inicial de agua (mm)
10
Cond hidráulica vertical (m/d)
0,1
Coeficiente de percolación (d−1)
0,231
Tiempo semi-agotamiento (d)
3
ACUÍFERO
Coeficiente de agotamiento (α) (d−1)
0,0036
Tiempo semi-agotamiento (d)
192,5
Coeficiente almacenamiento (porosidad drenable)
0,025
Nivel inicial (m snm)
43,63
Nivel de descarga (m snm)
30
Para el balance hidrometeorológico se ha utilizado la serie diaria de datos de precipitación (mm), temperatura media (°C), velocidad del viento (km·h−1), humedad relativa (%) y evapotranspiración potencial (mm) de la estación meteorológica PN Els Ports.
El balance de agua en el suelo se calcula teniendo en cuenta la infiltración de agua de lluvia, que es la precipitada menos la interceptación y la escorrentía superficial, descontando la evapotranspiración real calculada por el método exponencial (Samper 1997/2005), con un umbral de déficit hídrico de 10 mm (valor por debajo del cual la evapotranspiración real decrece respecto a la potencial). El excedente más el agua de precipitación que pasa por el suelo como flujo preferente por fisuras atraviesa descendentemente la base del suelo constituye la recarga en tránsito diferida que percola verticalmente por el medio no saturado y que en su recorrido se retrasa a través de variación del almacenamiento y en parte sale lateralmente dando lugar a la escorrentí diferida o lenta. El cálculo se realiza considerando que la salida de agua del almacenamiento en el suelo decrece exponencialmente (método modificado según Samper 1997/2005). Los resultados del balance de agua en el suelo en la cuenca del embalse de Ulldecona se muestran en la Tabla 3 (periodo de estudio 1994–2009). La media anual de la escorrentía superficial y de la zona vadosa es de 107 mm·a−1, suponiendo el 20 % de la recarga producida, en torno a 500 mm·a−1 (Figura 9).
Figura 9.—Hidrograma de precipitación versus escorrentía superficial y de la zona vadosa producida en los Ports (mm·d−1), derivado de la escorrentía calculada por el balance de agua en el suelo con el código Visual Balan.
Tabla 3. — Resultados medios del balance de agua en el suelo, en mm·a−1, para la cuenca receptora del Embalse de Ulldecona en el periodo octubre de 1994 a septiembre de 2009
Precipitación
Interceptación
Escorrentía superficial
Escorrentía zona vadosa
ETR
Recarga
Caudal total
1087
98
51
56
370
515
619
Para la calibración del modelo del balance de agua se ha utilizado la serie histórica de caudales diarios de entrada al embalse de Ulldecona. Los caudales totales calculados a partir del balance de agua en el suelo son considerablemente mayores que los medidos debido a que son la suma de los caudales de escorrentía superficial, escorrentía en zona vadosa y flujo subterráneo en el punto de aforo considerado (Figura 10). Según el modelo conceptual de funcionamiento considerado una parte del embalse de Ulldecona está emplazada por encima del nivel freático del acuífero regional y, el caudal aforado corresponde sólo a las entradas por escorrentía superficial y de la zona vadosa y una pequeña fracción de la recarga, ya que la casi totalidad de la recarga es transferida lateralmente a la Plana de La Galera; lo que está favorecido por la disposición alargada transversal del embalse. De esta forma, los caudales de entrada acumulados calculados y medidos (2438 m3·s−1 y 2240 m3·s−1 respectivamente) de escorrentía del suelo superior y de los acuíferos colgados temporales para el periodo 1994–2009 no difieren de forma significativa.
Figura 10.—Evolución temporal de los caudales diarios medidos y calculados (m3·s−1) y de la precipitación (mm·d−1). No se han representado los caudales >50 m3·s−1.
Para determinar de la bondad del ajuste entre los caudales de entrada del embalse con la escorrentía superficial más la vadosa calculada por el modelo se ha aplicado el Criterio D de Schultz (Cabrera, 2012), el cual representa la desviación de los caudales calculados respecto de los medidos:
donde Qcal,i es el valor simulado, Qi es el valor medido, Qmax es el caudal máximo observado en el periodo de estudio y n el número de datos de la serie. D puede variar entre 0 y 18, considerándose que el ajuste es muy bueno cuando D queda en el intervalo 0 a 3. El valor obtenido es D=0,008.
Se ha realizado un análisis de sensibilidad de los diferentes caudales que componen el caudal total y de los otros componentes del balance frente a los parámetros que presentan mayor incertidumbre y que mayor influencia pueden tener en su variación (Figura 11). También se ha obtenido el coeficiente de sensibilidad (CS) equivalente al factor de ponderación en la estimación de la varianza del parámetro analizado, que es la pendiente de la curva de sensibilidad en el entorno del valor central calculado.
Figura 11.—Sensibilidad de los componentes del balance de agua frente a los parámetros de mayor influencia: (a) espesor equivalente del suelo; la escorrentía superficial (del suelo superior) no está afectada por el espesor equivalente del suelo según el modo de cálculo realizado, (b) número de curva. (c), coeficiente de flujo preferente, CRKD y (d) coeficiente de agotamiento de la zona vadosa.
La variación de la reserva máxima en el suelo se representa como espesor equivalente del suelo (Figura 11), ya que el programa no permite su consideración directa. La reserva útil máxima en mm equivale a 40 veces el espesor del suelo en m para una diferencia volumétrica entre la capacidad de campo y el punto de marchitez de 0,04.
El análisis de sensibilidad muestra que la escorrentía superficial sólo es sensible de forma importante al número de curva con un coeficiente de sensibilidad (CS) para valores de Número de Curva entre 45 y 55 es de 4 mm·a−1 por unidad de número de curva. El número de curva, definido por las características de cobertura de suelo (uso y pendiente), determina en gran parte el comportamiento hidrológico de la cuenca en relación con la capacidad de producir escorrentía en función de la magnitud del evento de precipitación que la genera. El coeficiente de flujo preferente y el coeficiente de agotamiento de la zona vadosa son otros dos parámetros que afectan sensiblemente al caudal de agua drenado de la zona no saturada y a la recarga.
Discusión y conclusiones
El modelo conceptual de funcionamiento hidrodinámico de la cuenca receptora del embalse de Ulldecona parte de la hipótesis de que el llenado del embalse es mayoritariamente debido a la aportación de dos tipos de flujo, uno superficial con un agotamiento rápido y otro de escorrentía de la zona vadosa, más lento, que corresponde al agotamiento de niveles acuíferos colgados temporales (Figura 6).
Dicha hipótesis se apoya en el estudio hidroquímico (Espinosa, 2014), en el cual se ha comparado la concentración de calcio, magnesio y cloruro de las surgencias situadas aguas arriba del embalse de Ulldecona, del agua de la zona superficial, media y profunda del embalse y de piezómetros y surgencias aguas abajo del frente de presa (Figura 12). En dicho estudio se aprecia una notable influencia del agua de lluvia en el agua muestreada en el embalse y en las surgencias localizadas aguas arriba del mismo, corroborando el rápido agotamiento de estos manantiales y de las laderas de la cuenca que vierten al mismo. También se refleja la contribución del terreno en las muestras del frente de presa y de manantial existente a pie de presa, con mayor concentración que las muestras del agua superficial del embalse pero conservando la marca del agua de lluvia. Esta composición hizo considerar una posible mezcla de aguas debida a la conexión hidráulica entre el embalse y el acuífero regional, siendo el flujo subterráneo el que estuviera contribuyendo a la recarga del embalse, pero el estudio de la comparación de la evolución de los hidrogramas de llenado y vaciado del embalse con las oscilaciones de niveles de piezómetros situados en el entorno (no aportado aquí) no refleja esta posibilidad.
Figura 12.—Esquema del funcionamiento hidrodinámico de la cuenca receptora del embalse de Ulldecona y localización de los puntos de muestreo de agua superficial y subterránea. R: recarga en los Ports de Beseit, 1: Font del Teix, 2: Font de la Fou, 6: Piezómetro, 7: Font dels Rossegadors, 8: Font de Sant Pere; A: muestra de la zona media del embalse, B: muestra de la zona profunda del embalse y D: muestra de agua superficial del embalse.
Con todo lo anterior puede definirse que el modelo de funcionamiento de la cuenca receptora del embalse de Ulldecona consiste en escorrentía rápida en el suelo superior y escorrentía más lenta en la zona vadosa como consecuencia del agotamiento de niveles acuíferos colgados temporales. Ambos flujos alimentan al embalse, pero gran parte de la recarga se transfiere a otras cuencas o a aguas abajo, ya que la mayor parte del embalse está situado por encima del nivel piezométrico del acuífero regional.
Este modelo de funcionamiento de escala local puede extrapolarse al resto de los Ports de Beseit ya que aunque el grado de karstificación en todo el macizo sea variable a nivel de cuenca de barranco, en general las características geomorfológicas e hidrodinámicas son similares (Espinosa et al., 2015). Según este modelo, en la vertiente SE de los Ports, esta escorrentía superficial y sub-superficial es transferida lateralmente por los materiales de piedemonte a la cuenca adyacente de la Plana de La Galera. En muy raras ocasiones alcanzan al río Ebro.
Una vez obtenidos los valores de escorrentía media anual con el estudio de los hidrogramas de llenado-vaciado del embalse de Ulldecona y caudales de agua superficial de la vertiente NW del macizo y con el balance de agua en el suelo en la cuenca receptora del embalse de Ulldecona, ha sido posible comparar los resultados obtenidos con el balance general de agua en el suelo (Espinosa, 2014). Dichos valores de escorrentía varían entre 95 y 115 mm·a−1 (Tabla 4). A continuación se hace una pequeña reflexión sobre la validez de dichos resultados en el sistema de los Ports de Beseit.
Tabla 4. — Valores de escorrentía media anual obtenida por los diferentes métodos de estimación. S: superficie en km2, A: aportación en mm·a−1, AT: aportación total, ASS: aportación escorrentía de la zona vadosa, AS: aportación superficial, AL: aportación de escorrentía lenta (flujo subterráneo)
Aportación media anual mm a−1
PERIODO
ID
ESTACIÓN
Años
S km2
AT
ASS
AS
AL
HIDROGRAMAS
1994–2010
9113
Figuerales
16
55
2
2
1994–2010
9110
Pena (Beceite)
16
49
52
48
4
1994–2007
9052
Matarraña (Ulldemó)
13
48
177
81
15
80
1994–2010
9153
Algars (Horta SJ)
16
115
180
161
19
1994–2009
Embalse Ulldecona
Entradas
15
126
103
73
30
Salidas
15
126
100
VISUAL BALAN
1994–2009
Embalse Ulldecona
15
126
107
56
51
1994–2011
Balance general de agua en el suelo (Espinosa, 2014)
SZ4
17
175
103
SZ5
17
138
95–113
Las aportaciones de la cuenca receptora del embalse de Ulldecona calculadas a partir del hidrograma de volumen de entradas al embalse son de 103 mm·a−1. La diferencia con las salidas es de 3 mm·a−1, lo que supone un error del 3% que es menor que el de las propias medidas hidrométricas y está dentro de los errores de estimación atribuibles a la evaporación en el embalse y posibles pequeñas infiltraciones en el vaso del embalse no identificadas. Es posible que existan errores en la medición o incluso correcciones introducidas en las series, de las cuales se desconoce la metodología aplicada.
Con el balance de agua en el suelo se ha obtenido como resultado de la aportación 107 mm·a−1, que es sensible al parámetro Número de Curva, con un coeficiente de sensibilidad de 4 mm a−1 por unidad de número de curva. El número de curva se ajusta en función de las características del suelo, por lo que un estudio más en detalle en que se distingan zonas en función de la cubierta vegetal podría reducir de la incertidumbre de la escorrentía.
Los resultados de los hidrogramas tienen una incertidumbre asociada considerable, propia del método, y a pesar de ello hay coincidencia de magnitudes con los balances de agua.
Las conversaciones con especialistas de la zona y las observaciones de campo apuntan a que la escorrentía superficial se produce sólo cuando se tienen determinadas condiciones del estado del suelo previas al evento de lluvia. Por ejemplo, se produce escorrentía si ha llovido durante algunos días consecutivos y la lluvia acumulada en dichos días supera los 100 mm,, como se observa en los meses de invierno de los años 1995, 1996, 1997, 2002, 2003, 2004 y 2006 (ver Figura 9). Si no ha llovido durante un dilatado periodo de tiempo y se produce un evento de precipitación superior a los 100 mm·d−1 también se produce escorrentía, como se observa en octubre de 2000 y en diciembre de 2007. Además, analizando las series hidrométricas de escorrentía y precipitación diaria, la ocurrencia de escorrentía en función de los eventos de precipitación tiene lugar en momentos concretos, los cuales responderían al comportamiento observado en campo (Figura 9).
Para la estimación de los recursos hídricos del sistema Ports de Beseit-Plana de La Galera (Espinosa, 2014) se planteó la hipótesis de que la desviación estándar de la escorrentía superficial podría suponer el 20% de su valor. Por tanto, la magnitud de la escorrentía superficial media anual que tiene lugar en los Ports de Beseit se evalúa como 105±20 mm·a−1. Este valor podría definirse con mayor exactitud realizando un estudio de la incertidumbre asociada a cada una de las variables que intervienen en la estimación de la escorrentía superficial mediante la estimación de sus varianzas.
Se ha comprobado que existen dos tipos de flujo que llenan el embalse, uno superficial con agotamiento rápido, con un tiempo de semi-agotamiento medio de unos 6 días, y otro de escorrentía de la zona vadosa más lento con un tiempo de semi-agotamiento de unos 35 días, que corresponde a la descarga de los niveles acuíferos temporales colgados en el medio no saturado. Este modelo de funcionamiento de escala local es extrapolable al resto de los Ports de Beseit ya que, aunque el grado de karstificación en todo el macizo sea variable, en general las características geomorfológicas e hidrodinámicas son similares (Espinosa et al., 2015).
La escorrentía media anual estimada en la cuenca receptora del embalse de Ulldecona es comparable y de la misma magnitud cuando se calcula a partir de los caudales de entrada al mismo (103 mm·a−1) o por balance de agua en el suelo con el código Visual Balan (107 mm·a−1), siendo todos ellos a su vez coherentes con el rango de escorrentía media anual del balance de agua en el suelo con el código Visual Balan aplicado en toda la superficie de los Ports de Beseit (95 a 115 mm·a−1).
El análisis de sensibilidad de las variables del balance de agua en el suelo en la cuenca receptora del embalse de Ulldecona muestra que la escorrentía depende especialmente del número de curva seleccionado con un CS de 4 mm·a−1 por unidad de NC. La realización de un análisis de las varianzas de las variables que intervienen en el cálculo de la escorrentía superficial contribuiría a la disminución de su incertidumbre.
Como el valor de la escorrentía superficial media anual se evalúa en 105±20 mm·a−1, la escorrentía superficial representa entre un 20–30% de la recarga calculada a partir del balance general de agua en el suelo y del balance de la deposición atmosférica de ión cloruro, en torno a 500±60 mm·a−1 y a 350±150 mm·a−1 respectivamente (Espinosa 2014). La importancia de tener en cuenta o no la escorrentía superficial en la estimación de recursos de una cuenca de estas características queda patente ya que afecta significativamente a los resultados de recarga por la precipitación, sobreestimándolos, sobre todo en la aplicación del balance de la deposición atmosférica de ión cloruro (Espinosa, 2014).
AGRADECIMIENTOS
La investigación se enmarca en el desarrollo de la beca de FPI (CGL2009-12910-C03-01) del primer autor, dentro del proyecto REDESAC (MICINN CGLl2009-12910-C03-01). Los autores agradecen la disponibilidad de todas las personas que han hecho posible el desarrollo de los trabajos de campo, tanto en instalación de captadores de agua de lluvia como en localización y muestreo de puntos de interés. También se agradece a los laboratorios de análisis químico elemental e isotópico ambiental la eficiencia de su trabajo y aportación de ideas, en especial a Santiago Balagué, químico profesional y a Albert Solé de la Universidad de Barcelona. Diversos expertos en la zona han proporcionado notables datos e información, en especial Teresa Carceller, Javier San Román y Miguel Ángel García Vera de la Confederación Hidrográfica del Ebro, Javier Lambán y Antonio Azcón de la Unidad del IGME en Zaragoza y Josep Torrens, hidrogeólogo de Tarragona. Se agradece también la ayuda para acceder a los datos a la Confederación Hidrográfica del Ebro y de la Confederación Hidrográfica del Júcar.
