Tesis
Tesis
Tesis
Tesis Doctoral
2012
El Dr. Ignacio Morell Evangelista, Catedrático de Universidad de Hidrogeoquímica de
la Universitat Jaume I, el Dr. David Pulido Velázquez, Científico Titular de Organismo
Público de Investigación del Instituto Geológico y Minero de España y el Dr. Jorge
Mateu Mahíques, Catedrático de Universidad de Matemáticas de la Universitat Jaume
I.
Fdo. Dr. Ignacio Morell Evangelista Fdo. Dr. David Pulido Velázquez
Previamente a la defensa de la Tesis Doctoral, este trabajo ha sido evaluado por dos
censores independientes directamente relacionados con el área de investigación, Dr.
Antonio Pulido Bosch (Universidad de Almería) y Dra. María Vicenta Esteller Alberich
(Universidad Autónoma del Estado de México).
Gracias a mis directores de tesis: al Dr. Ignacio Morell, por guiarme por el mundo de
la hidrogeología, por su confianza, dedicación y ayuda, pero sobre todo por su
paciencia a la hora de dirigir este trabajo. Al Dr. David Pulido, por enseñarme a
entender y a “disfrutar” de los modelos matemáticos aplicados a la hidrogeología. Al
Dr. Jorge Mateu, por sus aportaciones matemáticas en los inicios de este trabajo, sus
buenas palabras y consejos.
A Liliana y Sandra, compañeras de Recursos Hídricos, agradecerles sus
conversaciones en el día a día y apoyo en las salidas al campo.
A los organismos y entidades que han financiado esta investigación en algún
momento: la Fundación Marina d´Or, por otorgarme una beca de investigación (2007 –
2008), y al Instituto Geológico y Minero de España (IGME), por el convenio firmado
con la Universitat Jaume I (2009).
A la Dra. Carolina Guardiola (IGME Madrid) por dedicarme su tiempo y darme
algún que otro “truco” para conseguir que los modelos corran. A la Dra. Mª Luisa
Calvache, del departamento de Geodinámica de la Universidad de Granada, por
iniciarme en el mundo de la simulación.
A todas las personas de Alcalà de Xivert que han prestado su colaboración (de
forma desinteresada) en la recopilación de datos para llevar a cabo el modelo
matemático de la Plana, a las Comunidades de Regantes de: Estopet - Torràs,
Mesquida, Bassa d´Amet, Palaba - Estopet, Racó Vale, Foia, Morito y Tossalet y de
Torreblanca (Ayuntamiento de Torreblanca). Un agradecimiento especial a Sandra por
acompañarme a todos los sitios y ayudarme en las tareas de campo. A Julián (EPTISA),
por permitirme acompañarle a tomar datos para la red de control de piezometrías y
calidad de la Confederación Hidrográfica del Júcar.
A mis amigas/o de la Universidad de Zaragoza por estar ahí siempre, las 24h del
día, los 365 días al año, desde hace más de 10 años. A mis amigas/os de Castellón, y en
especial a Anna Marzá, por interesarse por la evolución de la investigación y darme
ánimos para poder terminar la memoria.
A mi familia, por su incondicional apoyo (aunque les cueste creerlo si leen esto
es que ya lo he terminado). Y finalmente un agradecimiento especial a Aleix y Barthe,
por aguantar, por ser como son y estar siempre a mi lado.
“El agua no se agota
si se extrae del pozo
de la sabiduría humana”
Koichiro Matsuura
Resumen
En los últimos años, propiciado por la entrada en vigor de la Directiva Marco del Agua
(UE - DMA, 2000; UE - DMA 2006), están proliferando los estudios referidos a la gestión y
uso sostenible de las aguas subterráneas (IGME - UJI, 2009b; Gómez et al., 2012; Meléndez
et al., 2012). Dichos informes tienden a definir y cuantificar las salidas subterráneas al mar
de los acuíferos costeros y presentan propuestas metodológicas interesantes para su
determinación. Pero es sabido que el caudal subterráneo de salida al mar es el resultado de
todas las acciones llevadas a cabo en un acuífero; por tanto, la contención de la intrusión
marina no debería centrarse exclusivamente en determinar el caudal de salida sino en
cuantificar, y rectificar en la medida de lo posible, las acciones desencadenantes del
problema. En esta tesis se propone una nueva metodología para determinar el volumen
ecológico de remediación (VER), definido como el volumen en que hay que reducir los
bombeos en un acuífero o sector del mismo para que la recuperación piezométrica
alcanzada comporte una intrusión menor o igual a la estimada en régimen natural para el
sistema (intrusión de referencia).
In the last years, and following the application of the Water Framework Directive (UE -
DMA, 2000; UE - DMA 2006), research on water management and sustainable use of
groundwaters has increased (IGME - UJI, 2009b; Gómez et al., 2012; Meléndez et al., 2012).
Such reports tend to define and quantify the submarine groundwater discharge of coastal
aquifers and present interesting methodological frameworks to determine it. However, it is
well known that the submarine groundwater flow is the result of all the actions taken on an
aquifer. Therefore, the study of marine intrusion should not be solely focused on
determining the optimal groundwater discharge/flow, but also on quantifying, and
rectifying if possible, the actions that led to the problem. The objective of this research is to
define a new methodology to determine the ecological repair volume (ERV or VER in the
spanish version). ERV is the volume of pumping reduction that has to be reached in an
aquifer, or a sector of it, to yield a volume of intrusion lesser or equal to the estimated
intrusion for the system’s natural state (reference intrusion).
The methodological proposal for the calculation of ERV stems from a constant-density
mathematical flow model calibrated on a steady state. It consists in the iterative reduction
of pumping involved in the intrusion process, until the water table yields an intrusion
equal or lower than the estimated value for the system’s natural state (reference intrusion).
The following steps should be followed:
- PHASE 1: Modelling the system’s natural state to establish the reference intrusion.
The resulting model will be named reference model in a natural state (NS model, or
modelo RN in the spanish version).
- PHASE 2: Delimiting the areas affected by marine intrusion and trying to reduce
their surface by iteratively reducing the pumping involved in the process, until the
water table yields an intrusion equal or lower than the estimated value for the
system’s natural state. The initial model in this phase will be named as reference
model in an influenced state (IS model, or modelo RI in the spanish version).
- PHASE 3: Calculating the ERV from the IS model and the last scenario of Phase 2
(the one associated to an intrusion equal or lower than the estimated value for the
natural state).
The ecological repair volume (ERV) can be calculated for humid (ERVh), medium (ERVm),
and dry (ERVd) years. The application of this methodology in the Oropesa – Torreblanca
Plain (MASub 080.110) shows that the advance in seawater intrusion (with respects to the
natural state intrusion), in any of the three scenarios (humid, medium or dry years) is
mainly due to pumping. On average, in humid years the sector affected by a higher
intrusion is the Southern sector, and the corresponding ERV is 1 hm3/year. In medium
years, the ERVm is 1.6 hm3/year and the process affects exclusively the Southern sector. In
dry years the process affects the Southern and Central sectors; in this case, the ERVd would
be 2.8 and 0.9 hm3/year respectively.
Before the application of this methodology, a thorough analysis of the characteristics and
behaviour of the Oropesa – Torreblanca Plain (MASub 080.110) aquifer was conducted.
The evolution and present situation of the seawater intrusion of the Plana were studied
based on the information collected by the intrusion and quality control network of the
Confederación Hidrográfica del Júcar (CHJ, www.chj.es). A specific methodology has been
applied to evaluate the actions and vulnerability of the system affected by intrusion (SITE
index). The conclusion is that, by the end of the XXth century, the process affected the
whole aquifer and so the pressure due to the intrusion is classified as high.
A mathematical groundwater flow model (MODFLOW) has been used to evaluate the
actions that affect the system in the period 1973 -2009. The parameters of the model were
calibrated, and the water table evolution during the simulation was recorded. In 1993 the
intrusion affected two thirds of the coastal strip with a flow of 2.3 hm3/year. By the end of
the stated period, year 2009, the estimated intrusion affected only the Southern sector and
its volume was 0.05 hm3/year.
A feasible future scenario has been simulated (2010 - 2035). It was based on the urban
development regulations to be applied in short and medium terms in that area: the Local
Development Framework of Torreblanca (PGOU Torreblanca, 2009) and, within this
framework, the Integrated Development Plan (IDP, or PAI on the original) Doña Blanca
Golf and the IDP Marina d’Or Golf (in the municipalities of Cabanes and Oropesa). The
impact on that area’s water table due to changes in the land use was assessed. We have
obtained that there would not be water levels below the sea since 2018.
Índice
ÍNDICE GENERAL
CAPÍTULO 1. INTRODUCCIÓN .
1.1. Presentación…………………………………………………………………………...
35
1.2. Motivación y objetivos………………………………………………………………..
37
1.3. Localización geográfica………………………………………………………………
38
1.4. Rasgos geológicos básicos……………………………………………………………
40
1.5. Localización hidrogeológica…………………………………………………………
44
1.6. Antecedentes…………………………………………………………………………..
46
1.6.1. Antecedentes geológicos………………………………………………………..
46
1.6.2. Antecedentes hidrogeológicos………………………………………………….
46
1.7. Metodología básica……………………………………………………………………
48
1.8. Estructura de la memoria…………………………………………………………….
49
CAPÍTULO 2. HIDROGEOLOGÍA .
2.1. Introducción………………………………………………………………...…………. 53
2.2. La masa de agua subterránea 080.109: Maestrazgo Oriental……………………... 53
2.2.1. Caracterización hidrogeológica………………………………………………… 53
2.2.2. Balance hídrico…………………………………………………………………… 57
2.2.3. Hidroquímica general…………………………………………………………… 59
2.3. Acuíferos de borde de la Plana………………………………………………………. 60
2.3.1. Acuífero de materiales jurásicos………………………………………………. 61
2.3.2. Acuífero de materiales aptienses (Cretácico inferior)……………………….. 63
2.3.3. Acuífero en materiales del Cretácico superior………………………………… 64
2.3.4. Acuíferos en materiales del Mioceno………………………………………….. 65
2.4. La masa de agua subterránea 080.110: Plana de Oropesa - Torreblanca………... 66
2.4.1. Situación y naturaleza………………………………………………………….. 66
2.4.2. Geometría y límites…………………………………………………………….. 67
2.4.3. Parámetros hidrodinámicos…………………………………………………… 68
2.4.4. Piezometría……………………………………………………………………… 70
XVII
Índice
3.1. Antecedentes…………………………………………………………………………. 81
3.2. Evolución de la intrusión marina en la Plana de Oropesa - Torreblanca………. 83
3.3. Determinación del índice SITE……………………………………………………... 96
3.3.1. Cálculo de la superficie afectada (parámetro S)……………………………. 97
3.3.2. Cálculo de la intensidad del proceso (parámetro I)………………………... 99
3.3.3. Cálculo de la temporalidad o estacionalidad (parámetro T)……………… 101
3.3.4. Criterio de evolución (parámetro E)…………………………………………. 103
3.3.5. Confección del índice SITE…………………………………………………… 104
3.4. Conclusiones…………………………………………………………………………..107
XVIII
Índice
XIX
Índice
BIBLIOGRAFÍA 221
XX
Índice
XXI
Índice
XXII
Índice
ÍNDICE DE FIGURAS
XXIII
Índice
XXV
Índice
XXVI
Índice
XXVII
Índice
XXVIII
Índice
ÍNDICE DE TABLAS
XXIX
Índice
XXX
Índice
XXXI
Índice
XXXII
Capítulo 1 Introducción
Capítulo 1
Introducción
1.1. Presentación
1.2. Motivación y objetivos
1.3. Localización geográfica
1.4. Rasgos geológicos básicos
1.5. Localización hidrogeológica
1.6. Antecedentes
1.6.1. Antecedentes geológicos
1.6.2. Antecedentes hidrogeológicos
1.7. Metodología básica
1.8. Estructura de la memoria
33
Capítulo 1 Introducción
34
Capítulo 1 Introducción
1. INTRODUCCIÓN
1.1. Presentación
Los acuíferos costeros presentan peculiaridades hidrodinámicas e hidrogeoquímicas
que añaden complejidad a la correcta gestión de sus recursos hídricos subterráneos. En
ellos, el equilibrio agua dulce - agua salada puede verse alterado (a veces muy
notablemente) de forma natural o antrópica al producirse variaciones en las componentes
del balance.
Dado que los acuíferos costeros, por otra parte, suelen albergar núcleos importantes de
población y soportar actividades agrícolas y/o industriales, la demanda de agua puede ser
muy elevada. Si, como suele ser frecuente en el litoral mediterráneo, los recursos de aguas
superficiales son escasos o intermitentes, la garantía de suministro se basa total o
parcialmente en la explotación de aguas subterráneas. El escenario más frecuente es que se
alcance cierto grado de sobreexplotación, entendida ésta como la explotación de un
volumen de agua superior al de los recursos renovables, que provoca el progreso del agua
de mar hacia el acuífero dando lugar a la salinización más o menos severa de sus aguas.
Para evitar que el proceso de intrusión marina progrese tierra adentro, además de
otras prescripciones, es necesario permitir que las salidas de agua dulce al mar se
mantengan. En contra de la creencia de que es “agua perdida” se trata realmente de
recursos ganados, al menos a medio plazo. Para reducir el efecto de la intrusión en las áreas
más intensamente afectadas es necesario reducir los volúmenes de explotación, evitando la
existencia de niveles por debajo del nivel del mar.
Naturalmente, si la demanda no se reduce es necesario disponer de nuevos recursos
para compensar el cierre de las captaciones. La normativa vigente hace hincapié en la
procedencia del agua que debe de ser utilizada para satisfacer las necesidades hídricas de
las áreas costeras; la reutilización de las aguas residuales regeneradas procedentes de las
EDARs (estación depuradora de aguas residuales) es la opción para suplir las necesidades
agrícolas; en cambio, para la demanda urbana el agua procedente de la desalación es la
alternativa propuesta.
35
Capítulo 1 Introducción
Por otra parte, en algunos sectores mediterráneos se ha producido en los últimos años
una reducción de la explotación de aguas subterráneas, acorde a la disminución de la
demanda hídrica experimentada como consecuencia de tres factores fundamentales:
aplicación de políticas de ahorro y eficiencia en el uso del agua, reconversión de los
sistemas de riego y abandono de superficies de cultivo. Este último factor se ha acentuado
en la última década debido a la presión por la construcción de infraestructuras, a la mala
calidad del agua y al cambio de usos del suelo (de agrícola a residencial, principalmente).
Algunos de estos cambios pueden ser muy intensos, de forma que repercuten
notablemente en las condiciones hidrogeológicas locales dando lugar a situaciones nuevas
o, al menos, a nuevas tendencias en el reparto de los componentes del balance.
36
Capítulo 1 Introducción
consecución de los siguientes objetivos parciales, como son la historia evolutiva del
acuífero, su situación actual y valorar alternativas de evolución futura. Para ello se ha
desarrollado un modelo matemático de flujo que ha contribuido a mejorar el conocimiento
del sistema y valorar diferentes escenarios de gestión.
38
Capítulo 1 Introducción
Los mapas topográficos del Institut Cartogràfic Valencià 1:10000 (ICV, 2008) que
abarcan el área de estudio son: Alcossebre (594(2-2)), Alcossebre (594(2-3)), Capicorb
(594(2-4)), Torreblanca (594(1-3)), Torrenostra (594(1-4)), Torre la Sal (617(1-1)), Orpesa
Llevant (617(1-2)), Tossal del Mortorum (616(4-1)) y Orpesa (616(4-2)), así como en las hojas
topográficas 1:50.000 (ENADINSA, 1973a; ADARO 1973a y 1973b) de Alcalá de Chivert
(594(31-23)), Cuevas de Vinromá (593(30-23)) y Villafamés (616(30-24)).
La red hidrográfica es escasamente activa ya que no existen cauces de funcionamiento
permanente y ni siquiera estacional. Los dos cauces más importantes son el barranco
Chinchilla, que desemboca al norte de la población de Oropesa, y el río Cuevas o San
Miguel, que lo hace al norte de Torrenostra (figura 2), ambos son de funcionamiento
esporádico. Otros cauces de cierto interés son los barrancos de Perchets y Campello y las
ramblas de Mañes y Estopet, que sólo están activos a causa de fuertes eventos lluviosos.
En el sector central costero de la Plana se encuentra la zona húmeda denominada Prat
de Cabanes (figura 2), que está considerada como una de las mayores marismas
39
Capítulo 1 Introducción
valencianas. Fue declarada Paraje Natural en 1988, y desde ese año pertenece a la lista
ZEPA (Zona de Especial Protección para las Aves); en 1989 se incluyó en la lista RAMSAR
(Humedales de importancia internacional). Se extiende desde Torrenostra hasta Torre la
Sal, con una anchura media de 1,5 km. Es un área permanentemente inundada, en relativo
proceso de colmatación, con abundante vegetación, limitada por una barra costera de casi 8
km de longitud, unos 20 m de anchura y hasta 3 de altura, constituida por cantos muy
redondeados, bien clasificados, con tamaños que oscilan entre 1 y 20 cm.
40
Capítulo 1 Introducción
general, la cordillera está formada por una serie de macizos mesozoicos, donde se pueden
observar núcleos paleozoicos aislados, individualizados por cubetas rellenas de materiales
terciarios y cuaternarios.
Geológicamente, la provincia de Castellón se sitúa en la unión del borde más oriental de
la Rama Aragonesa con el sector meridional de la Cordillera Costero - Catalana,
configurando la transición entre ambas. Las características estructurales de la Cordillera
Costero - Catalana se aprecian en los sectores oriental y nororiental, con direcciones de
plegamiento NNE - SSO, que condicionan la actual línea de costa (Giménez, 1994). Es
indudable la influencia ibérica en la región, dirección NO - SE, pero son mucho más
visibles las directrices catalanas NNE - SSO, posiblemente por ser ligeramente posteriores y
haber cortado a las ibéricas (ENADINSA, 1973). En la figura 3 se observa la localización
geológica regional de la zona de estudio (Canerot, 1974); el círculo rojo indica la situación
de la Plana de Oropesa - Torreblanca.
41
Capítulo 1 Introducción
En esta región, Sanchis et al. (1989) y posteriormente Morell (1992), diferenciaron cinco
dominios, basados en las características descritas por Canerot en 1974 (figura 4):
-Pliegues septentrionales: Este dominio se sitúa al noroeste de la provincia de
Castellón. Destaca por los materiales jurásicos y cretácicos, que ocasionalmente aparecen
cubiertos por depósitos terciarios. Estructuralmente se caracteriza por una sucesión de
anticlinales y sinclinales.
-Zona tabular del Maestrazgo: Ocupa la parte interior central de la provincia de
Castellón, limitando al este con las Alineaciones Costeras, al norte con la región de Pliegues
Septentrionales, y al sur con el dominio Triásico del Espadà. Se caracteriza
estructuralmente por los pliegues de gran radio, conformando estructuras claramente
tabulares. Estratigráficamente predominan los materiales cretácicos, desarrollando buenas
secuencias.
-Alineaciones costeras: Abarca la zona oriental de la provincia. Corresponde a la
parte meridional de la Cordillera Costero - Catalana. Se caracteriza por una sucesión de
fosas y elevaciones alargadas, dirección NNE - SSO, generados por el enrejado de fracturas
presente en la zona, donde la complejidad va aumentando con la cercanía de la costa.
-Dominio triásico del Espadà: Conforma la parte meridional de la provincia. Limita
al norte con la zona Tabular del Maestrazgo y al este con la Llanura Costera de la Plana de
Castellón. El núcleo principal de este domino es la Sierra de Espadà, donde predominan los
materiales triásicos y afloran esporádicamente materiales paleozoicos.
-Llanuras costeras: Dentro de ella se enmarcan las tres llanuras detríticas costeras
que se asientan a lo largo de todo el litoral castellonense, que son la Plana de Vinaroz, la
Plana de Oropesa - Torreblanca y la Plana de Castellón. En ellas existe un predominio de
materiales sedimentarios continentales, en general, provenientes del transporte fluvial,
aunque también se observan materiales mixtos y marinos.
42
Capítulo 1 Introducción
La Plana de Oropesa - Torreblanca limita la norte con la Sierra de Irta, al sur con el
Macizo del Desert de les Palmes y al este por las estribaciones meridionales del Macizo del
Maestrazgo (figura 1).
La Sierra de Irta está constituida por dos ramas separadas por una falla de dirección
NE - SO. La falla principal presenta un salto de 500 m. En cambio, las dos ramas, antes
nombradas, presentan una series de fallas menores, de orientación semejante, que
conjugadas con las fallas de dirección ibérica dan lugar a la formación de pequeños bloques
(ENADINSA, 1973).
El Macizo del Desert de les Palmes incluye las Sierras del Desert de les Palmes, les
Santes, Oropesa, Borriol y Monegre. La estructura principal en este macizo es un anticlinal
de dirección ibérica, donde el núcleo viene definido por los materiales más antiguos, de
edad paleozoica, que aparecen recubiertos por una serie triásica poco potente, mediante
una discordancia angular (Canerot, 1974). Esta estructura ibérica está afectada por el
conjunto de fallas de orientación catalánide, NNE, generando un conjunto de bloques sin
hacer desaparecer la estructura principal anticlinal (Renau - Llorens, 2010).
Tectónicamente, las estribaciones meridionales del Macizo del Maestrazgo, presentan
una característica fundamental que viene definida por una densa red de fallas verticales, de
orientación catalánide, que compartimentaron la región en fosas y pilares donde,
posteriormente, las fosas sirvieron de base para la sedimentación de los materiales
terciarios (ADARO, 1973b).
43
Capítulo 1 Introducción
44
Capítulo 1 Introducción
MATERIALES
LÍMITE TOPONÍMIA EDAD
MAYORITARIOS
45
Capítulo 1 Introducción
1.6. Antecedentes
1.6.1. Antecedentes geológicos
El autor que más profundizó en el conocimiento geológico del Baix Maestrat fue
Canerot en su Tesis Doctoral (1974), la cual ha servido de base para la mayoría de los
trabajos realizados posteriormente. En 1975, EPTISA (1975a y 1975b) realizó para el
Instituto Geológico y Minero de España y para el Instituto de Reforma y Desarrollo Agrario
(IRYDA) un informe donde se recopilaron los resultados de realizar sondeos con testigo
continuo, que permitieron conocer con bastante aproximación la litoestratigrafía de la zona.
Los aspectos geomorfológicos fueron tratados por Mateu (1982), que aborda las
características del relieve del litoral. En 1985, Morell contribuyó a conocer nuevos datos
sobre la geología de la zona. A partir de este momento y hasta la actualidad no se conocen
informes significativos donde se aporten nuevos datos sobre la geología del área de
estudio.
46
Capítulo 1 Introducción
47
Capítulo 1 Introducción
B) Aspectos hidrogeoquímicos
Se ha recopilado bibliografía sobre hidrogeoquímica de la zona, con especial
énfasis en lo que refiere al proceso de intrusión marina que afecta a la misma. El estudio
de la evolución temporal y la distribución espacial de la intrusión marina se ha
realizado con los datos proporcionados por las redes de control de intrusión y calidad
de la CHJ (www.chj.es).
C) Modelización
Se ha desarrollado un modelo en diferencias finitas MODFLOW (McDonald y
Harbaugh, 1988) del acuífero pliocuaternario de la Plana de Oropesa - Torreblanca para
simular el flujo subterráneo (Visual Modflow) en la zona saturada. Se recopilaron los
datos existentes, y en los casos que fue necesario se depuraron para obtener información
fiable y representativa del área de estudio para diseñar el modelo. El modelo fue
calibrado en régimen transitorio para el periodo comprendido entre 1973 y 2009 y
verificado en régimen estacionario representativo del periodo 1998 - 2001. Se analizó la
sensibilidad de los parámetros y acciones que lo conforman.
está incluido Doña Blanca Golf) para conocer las necesidades hídricas futuras de la
zona. Con esta información, se ha definido el escenario futuro y se ha valorado los
impactos asociados al mismo simulándolo con el modelo previamente calibrado.
49
Capítulo 1 Introducción
50
Capítulo 2 Hidrogeología
Capítulo 2
Hidrogeología
2.1. Introducción
2.2. La masa de agua subterránea 080.109: Maestrazgo Oriental
2.2.1. Caracterización hidrogeológica
2.2.2. Balance hídrico
2.2.3. Hidroquímica general
2.3. Acuíferos de borde de la Plana
2.3.1. Acuífero de materiales jurásicos
2.3.2. Acuífero de materiales aptienses (Cretácico inferior)
2.3.3. Acuífero de materiales del Cretácico superior
2.3.4. Acuíferos en materiales del Mioceno
2.4. La masa de agua subterránea 080.110: Plana de Oropesa-Torreblanca
2.4.1. Situación y naturaleza
2.4.2. Geometría y límites
2.4.3. Parámetros hidrodinámicos
2.4.4. Piezometría
2.4.5. Balance hídrico
2.4.5.1. Entradas
2.4.5.2. Salidas
2.4.5.3. Balance
2.4.6. Rasgos hidroquímicos básicos
2.5. Importancia ambiental del Prat de Cabanes
2.5.1. Datos básicos
2.5.2. Caracterización paleográfica del Prat
51
Capítulo 2 Hidrogeología
52
Capítulo 2 Hidrogeología
2. HIDROGEOLOGÍA
2.1. Introducción
La Plana de Oropesa – Torreblanca es una unidad geomorfológica e hidrogeológica
definida, aunque integrada en un conjunto geológico más amplio, al que está
estrechamente vinculada. Dada su reducida extensión y los bajos valores de precipitación
media anual que se registran en este sector, la mayor parte de la recarga del acuífero tiene
lugar por transferencias laterales subterráneas. De ahí la importancia de conocer las
relaciones estructurales y las conexiones hidráulicas existentes entre el acuífero
pliocuaternario y los materiales de borde.
53
Capítulo 2 Hidrogeología
Los materiales acuíferos principales están formados por calizas dolomíticas y dolomías
del Lias - Kimmeridgiense inferior, y calizas del Kimmeridgiense superior - Barremiense,
separados por una serie margosa (130 metros) del Kimmeridgiense superior. En el sector
meridional, también conforman el acuífero las calizas del Aptiense. Los materiales que
producen la separación hidráulica del acuífero regional están constituidos por una
alternancia de calizas y margas con horizontes arenosos del Valanginiense – Bedouliense y
arenas, areniscas y arcillas del Albiense.
El acuífero que presenta mayor interés corresponde a los materiales carbonatados del
Jurásico - Cretácico inferior, y se encuentra muy plegado y fracturado posibilitando la
conexión hidráulica de sus niveles permeables. El impermeable de base está constituido
por margas y arcillas, en facies Keuper, del Triásico superior. Sobre los materiales
carbonatados mesozoicos se disponen algunos rellenos cuaternarios, aluviales y glacis; y,
en menor proporción, calizas, margas y lutitas, miocenas.
valores de transmisividad más bajos que los mencionados en el sector oriental del acuífero,
con magnitudes variables entre 100 y 1000 m2/día, incluso se han perforado algunos
sondeos negativos (Marina, 2006).
El drenaje principal del sistema tiene una dirección genérica NNO - SSE,
aproximadamente perpendicular a la línea de costa (figura 7). Llama la atención la
existencia de un notable salto hidráulico a lo largo de un umbral de dirección NNE - SSO,
aproximadamente en la alineación de las poblaciones de Albocàsser, Tírig, San Mateo y la
Jana. Este umbral no significa que exista desconexión hidráulica ya que no se han
identificado formaciones impermeables que lo justifique (Mejías et al., 2007).
Probablemente se deba a la presencia de tramos acuíferos con diferentes cargas hidráulicas,
variable según la posición relativa de los niveles en la zona de umbral, o la existencia de
una mayor permeabilidad a ambos lados del mismo (Antón - Pacheco et al., 2007). Aguas
abajo de este umbral se aprecia una notable divergencia del flujo de manera que las líneas
de flujo del sector nororiental y la mayor parte de la descarga procedente de la Sierra de
Irta se dirige hacia los manantiales del Prat de Peñíscola y Torre Badum, así como de forma
subterránea hacia el acuífero de la Plana de Vinaroz (figura 8). Por otro lado, el flujo
subterráneo del área suroccidental se orienta hacia las surgencias de Alcossebre y hacia el
acuífero de la Plana de Oropesa – Torreblanca (Antón - Pacheco et al., 2007).
Así pues, el drenaje principal del acuífero del Maestrazgo tiene lugar a través de la
Sierras de la Irta y en menor medida por alimentación lateral a los acuíferos detríticos
instalados en las planas costeras de Vinaroz y Oropesa – Torreblanca. Estos últimos
recursos podrían ser finalmente transferidos de forma difusa al mar a través de las
formaciones mesozoicas infrayacentes a dichas planas costeras (Mejías et al., 2007).
55
Capítulo 2 Hidrogeología
56
Capítulo 2 Hidrogeología
San Mateo
N
Sierra de
Valdanch
1
Peñíscola
2
Cuevas Sierra
de Irta 3
Alcalá
4
5
Mar Mediterráneo
Según la CHJ (2005, 2007), los recursos renovables son del orden de 292 hm3/año,
obtenidas agregando tres parámetros: infiltración del agua de lluvia e infiltración de ríos
(270 hm3/año), retorno de riegos (15 hm3/año) y el 80% de las entradas laterales
subterráneas (8 hm3/año).
Para calcular los recursos disponibles es necesario detraer de los recursos renovables el
volumen medioambiental necesario para garantizar el buen estado ecológico de la masa de
agua. El volumen ambiental se calcula teniendo en cuenta el volumen ecológico (descargas
subterráneas a los ríos), el 80% del volumen de las descargas subterráneas laterales a otras
masas de agua subterránea y el volumen de las descargas al mar y a las zonas húmedas. En
este caso sólo se debe considerar como volumen ambiental el 80% de las salidas laterales
subterráneas, que es del orden de 111 hm3/año. En consecuencia, los recursos disponibles
son del orden de 181 hm3/año (CHJ, 2007).
57
Capítulo 2 Hidrogeología
Los bombeos totales del sistema se estiman en 37 hm3/año, de los que 22 son para usos
agrícolas, 14 para usos urbanos y 1 para usos industriales.
Con estos datos, el índice de explotación (K) de la masa de agua, que es el cociente entre
los bombeos y los recursos disponibles, es bajo (0,2).
En la tabla 2 se sintetiza el balance hídrico realizado por CHJ (2005, 2007) utilizando el
programa Patrical.
Tabla 2. Balance hídrico del acuífero del Maestrazgo Oriental (CHJ, 2005 y 2007).
Según este balance, la variación de la reserva es de 118,9 hm3/año. Este valor fue
estimado considerándose que las salidas al mar son nulas; sin embargo, como ya se ha
mencionado en el epígrafe anterior, son bien conocidas las descargas al mar a través de la
Sierra de Irta y, muy probablemente también bajo las Planas de Oropesa – Torreblanca y
Vinaroz.
58
Capítulo 2 Hidrogeología
59
Capítulo 2 Hidrogeología
60
Capítulo 2 Hidrogeología
En general, los materiales jurásicos presentan una alta permeabilidad por fisuración y
karstificación, especialmente los tramos de dolomías del Kimmeridgiense superior –
Portlandiense y los de calizas y dolomías del Lías – Dogger. Recientemente, en 2007, se han
realizado varios ensayos de bombeo en sondeos de investigación que han arrojado valores
de transmisividad para las formaciones jurásico – cretácicas del orden de 3000 m2/día, con
caudales específicos comprendidos entre 5 y 60 l/s/m (Antón – Pacheco et al., 2007).
Los materiales jurásicos descansan sobre arcillas y margas con yeso pertenecientes al
Triásico superior, en facies Keuper, que sólo afloran puntualmente en la Sierra
Espaneguera (Mejías et al., 2006), en el sector meridional del Desierto de les Palmes sí han
61
Capítulo 2 Hidrogeología
62
Capítulo 2 Hidrogeología
Las transmisividades que caracterizan estos materiales varían desde 500 m2/día a 4
m2/día, que evidencia la heterogeneidad hidráulica del acuífero, seguramente ligada a su
funcionalidad de índole kárstica. El coeficiente de almacenamiento es del orden de 10-2
(Morell, 1985).
63
Capítulo 2 Hidrogeología
Estos materiales se sitúan en el sector central del borde de la Plana, en las inmediaciones
de Torreblanca, en el Pico Ferradura y en las Sierras de Perchets y Sernavives. También se
pueden reconocer en el sector norte, en el barranco de Estopet (figura 9).
El macizo calcáreo principal del Cretácico superior, Sierras de los Perchets y Sernavives,
constituye un importante acuífero kárstico que se extiende hacia el interior de la provincia.
La alimentación que recibe, en rasgos generales, procede de la infiltración de las
precipitaciones, pero también cabe la posibilidad de que reciba recarga lateral subterránea,
a través de sus acuíferos de borde situados al norte y noroeste, a través del pasillo de Torre
En Domenech e, incluso, desde la Sierra de Valdancha oriental. En cambio, en el
afloramiento que se observa en las inmediaciones de Torreblanca, la alimentación debe
proceder de los conglomerados miocénicos por existir continuidad hidráulica. En el sector
norte, los materiales que se observan en el barranco de Estopet, son alimentados por el
acuífero jurásico (Morell, 1989).
64
Capítulo 2 Hidrogeología
Este nivel lacustre constituye un pequeño acuífero confinado en el que se llega a captar
agua de mar en un pozo situado cerca de la desembocadura del río Chinchilla para el
abastecimiento a la planta desaladora de Marina d’Or (Morell, 2008).
Sólo en los tres sondeos mencionados y en otros sondeos de investigación realizados por
el Servicio Geológico de Obras Públicas en los años setenta (SGOP, 1977) se han
caracterizado los materiales miocenos bajo el relleno cuaternario. En otros casos, los
materiales miocenos están ausentes y el conjunto pliocuaternario reposa directa y
discordantemente sobre materiales mesozoicos.
65
Capítulo 2 Hidrogeología
66
Capítulo 2 Hidrogeología
Figura 11. Isopacas del conjunto pliocuaternario (EPTISA 1975b; Morell, 1985 )
67
Capítulo 2 Hidrogeología
El impermeable de base lo constituyen las margas y arcillas del Mioceno, que incluyen
algunas intercalaciones de conglomerados y de calizas. Estos materiales, que afloran
extensamente al oeste de la población de Torreblanca (figura 9), han sido descritos en
algunas columnas de sondeos de investigación realizados por el Servicio Geológico de
Obras Públicas en la década de los setenta. Más recientemente, dos perforaciones
realizadas en la desembocadura del río Chinchilla, al norte de Oropesa, han permitido
identificar niveles de margas, calizas y conglomerados, que se corresponden con facies
marinas, lacustres y continentales, con un espesor total superior a 140 metros. Los niveles
arenosos miocenos constituyen un acuífero de cierto interés, que es explotado en otros
puntos del litoral castellonense, como Benicasim y Burriana.
La figura 12 presenta tres cortes esquemáticos (sur, centro y norte) que ilustran sobre la
geometría del acuífero pliocuaternario, la presencia del impermeable mioceno de base y la
relación con los materiales de borde.
68
Capítulo 2 Hidrogeología
69
Capítulo 2 Hidrogeología
2.4.4. Piezometría
En la figura 13 se comparan las geometrías de las superficies piezométricas que
caracterizan a situaciones de años secos, medios y húmedos. En general, la dirección del
flujo es del NO hacia el SE, prácticamente perpendicular a la línea de costa. La diferencias
esenciales entre las tres situaciones consideradas son las lógicas variaciones en los
gradientes hidráulicos y la distinta incidencia de situaciones de flujos invertidos.
70
Capítulo 2 Hidrogeología
2.4.5.1. Entradas
Retorno de riego
71
Capítulo 2 Hidrogeología
2.4.5.2. Salidas
2.4.5.3. Balance
tener una idea del orden de magnitud de cada una de las partidas. El balance es el
siguiente, tabla 3:
Tabla 3. Resumen del balance hídrico de la Plana Oropesa - Torreblanca (IGME - UJI, 2009b)
73
Capítulo 2 Hidrogeología
Está declarado como Zona de Especial Protección para las Aves (Zona ZEPA), de
acuerdo con la Directiva 79/409/CEE, relativa a la Conservación de las Aves Silvestres.
Este espacio está también incluido en el catálogo de Important Bird Areas in Europe,
confeccionado por las siguientes asociaciones: ICBP (International Council for Birds
Preservation), IWRB (International Waterfowl and Wet Lands Research Bureau) y RSPB
(Royal Society for the Protection of Birds).
74
Capítulo 2 Hidrogeología
El Prat, es una zona húmeda, que presenta una lámina de agua permanente en algunos
lugares y estacional en otros. La inundación de determinadas zonas tiene lugar cuando el
nivel freático alcanza la cota topográfica de las zonas deprimidas. Otros aportes de agua
de escasa importancia se producen como consecuencia de precipitaciones extraordinarias.
No existe en el Prat ningún curso fluvial de régimen continuo, aunque si se encuentran
manantiales de agua dulce que contribuyen a mantener zonas inundadas durante todo el
año.
Dicho humedal se puede clasificar como una zona húmeda del ámbito costero
mediterráneo dominado por la dinámica litoral (Durán et al., 2009). El predominio de la
dinámica litoral significa que su formación se concreta en el cierre de una bahía mediante
el crecimiento de una barra litoral que la separa total o parcialmente del mar.
Según Carmona y Ruiz (2009), la restinga de cantos se extiende entre los abanicos
aluviales del río Cuevas o San Miguel (Punta de Cap i Corp) y el abanico del río Chinchilla
a lo largo de ocho kilómetros. Intercalados entre el cordón de cantos y las facies distales de
los abanicos pleistocenos aparecen diversos afloramientos rocosos de arenisca (eolianitas)
discontinuos. Se localizan principalmente al sur del Quarter Vell, en la Torre de la Sal y en
la desembocadura del río Chinchilla. En la restinga se puede distinguir un tramo
septentrional de cantos (hasta el norte del Quarter Vell), un tramo central arenoso y un
tramo meridional (hasta la Torre de la Sal) de cantos, pero con una pequeña acera arenosa
en la base de la cresta y abundantes afloramientos de eolianita. El tramo central arenoso
podría explicarse por la refracción del oleaje en el centro de la bahía.
Nivel a: en la base (4 - 3,75 metros) hay un nivel muy compacto de color gris verdoso y con
nódulos de color gris y rosado que se identifican como facies distal del abanico pleistoceno.
Nivel b: desde 3,75 a 3 metros se individualiza un nivel de pantano continental estacional.
75
Capítulo 2 Hidrogeología
Los sedimentos son de textura arcillo limosa de color gris claro, presentan moteados ocres,
nódulos de color gris, restos vegetales y no se detecta la presencia de gasterópodos. Se
efectuó una datación de 14C a 3,40 metros que dio como resultado una fecha 8360 ± 50 BP
(cal BC 7530 a 7320). Este nivel puede corresponder a humedales continentales sobre la
superficie distal de los abanicos en las fases iniciales de la transgresión marina holocena.
76
Capítulo 2 Hidrogeología
77
Capítulo 2 Hidrogeología
78
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Capítulo 3
La intrusión marina en
la Plana de Oropesa - Torreblanca
3.1. Antecedentes
3.2. Evolución de la intrusión marina en la Plana de Oropesa - Torreblanca
3.3. Determinación del índice SITE
3.3.1. Cálculo de la superficie afectada (parámetro S)
3.3.2. Cálculo de la intensidad del proceso (parámetro I)
3.3.3. Cálculo de la temporalidad o estacionalidad (parámetro T)
3.3.4. Criterio de evolución (parámetro E)
3.3.5. Confección del índice SITE
3.4. Conclusiones
79
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
80
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Basado en el análisis de los datos proporcionados por las redes de control de intrusión,
el IGME (2000) actualizó los conocimientos sobre los acuíferos costeros de las cuencas del
Júcar, Segura e Islas Baleares. El informe destaca que la práctica totalidad del acuífero de
Oropesa - Torreblanca superaba en 1999 los 400 mg/L de Cl- y más de dos tercios los 800
mg/L. Diferencia tres zonas con valores superiores a 1200 mg/L de Cl-: al norte del
barranco de Chinchilla, al suroeste del núcleo urbano de Torreblanca, y la última desde la
desembocadura del río San Miguel extendiéndose rambla arriba. En los dos primeros casos
interpreta la existencia de domos salinos mientras que en el caso de río San Miguel parece
deberse al avance horizontal de la cuña salina.
La CHJ (Fidalgo et al., 2007) efectuó un estudio para la caracterización del estado
cualitativo de las masas de agua subterránea de la franja litoral de la cuenca del Júcar.
Según este estudio, la presión en la masa de agua subterránea 080.110 (Plana de Oropesa -
Torreblanca) conforme al índice de densidad de captaciones subterránea es significativa y
según el índice de explotación de la masa de agua subterránea costera es no significativa.
Por consiguiente, la presión global por intrusión marina según estos criterios es
81
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
significativa. Pero el estudio va más allá y evalúa los impactos con influencia de la intrusión
marina, siendo:
A la vista de los datos aportados por los estudios realizados hasta fechas muy recientes
se podría resumir que la intrusión marina en el acuífero pliocuaternario de la Plana de
Oropesa - Torreblanca es patente desde hace más de cuatro décadas y que ha ido
progresando hasta afectar a la totalidad del acuífero. Sin duda, la conversión de extensas
áreas de secano tradicional (almendros principalmente) a cultivos de regadío (cítricos) ha
sido la causa principal de la sobreexplotación que ha sufrido el acuífero; en menor medida,
también ha contribuido el auge turístico experimentado hacia finales del siglo XX. Sin
embargo, existen evidencias de que el proceso debe encontrarse en situación estacionaria e
incluso en regresión en algunos sectores ya que la modernización de los sistemas de riego,
que incluye el cambio de sistema de riego de inundación a goteo, el progresivo abandono
de las tierras de cultivo por la presión urbanística y el estancamiento del turismo han
reducido notablemente los volúmenes de bombeo, especialemente en la mitad sur del
acuífero.
82
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Se han diferenciado cuatro categorías de salinización, que son las siguientes (tabla 4):
El área de salinidad baja (< 200 mg/L de cloruros) se limita a algunos sectores del borde
occidental del acuífero, en los alrededores de Torreblanca, coincidiendo con las supuestas
áreas de alimentación preferencial desde el acuífero cretácico de borde. También se
identifica un sector al norte de la Sierra de Oropesa que coincide con el área de
alimentación del entorno del Chinchilla.
83
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Figura 16. Mapas de isocontenidos de cloruros, años: 1980, 1982, 1985 y 1987
84
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Figura 17. Mapas de isocontenidos de cloruros, años: 1990, 1995, 1998 y 2000
La distribución de las áreas de las categorías de moderada y alta salinización (200 – 1000 y
1000 – 2000 mg/L, respectivamente) ocupan la mayor parte de la superficie del acuífero y
presentan sólo ligeras variaciones en las distintas fechas consideradas.
85
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Sin duda, el rasgo más llamativo es la presencia de sectores con más de 2000 mg/L de
cloruros. El sector más intensamente afectado por la salinización se localiza en el entorno de la
desembocadura del río san Miguel, al norte de la zona húmeda del Prat de Cabanes, en el
periodo 1982 – 1987, alcanzando su máxima intensidad en 1985. Un segundo sector muy
salinizado se localiza al sur de Torreblanca a partir de 1990. Su reducido tamaño y su carácter
aislado sugieren que se trata o bien de una situación de intensa explotación local (upconing) o de
la existencia de un pozo de elevada profundidad que alcanza niveles diferencialmente más
salinizados.
Figura 18. Evolución temporal de las superficies afectadas por distintos grados de
salinidad
La superficie media ocupada por el rango de menos de 200 mg/L de cloruros es del
7,5%, prácticamente restringido a áreas marginales y con variaciones poco apreciables. Las
concentraciones comprendidas entre 200 y 1000 mg/L ocupan una superficie media del
43%. Se aprecian dos situaciones distintas, una en la primera mitad del periodo
86
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
considerado (1980 - 1990) en la que la superficie ocupada es del 39,9% y la otra el periodo
1990 – 2000 en el que alcanza el 47,7%. Este incremento coincide con la disminución relativa
del rango de 1000 – 2000 mg/L y, sobre todo, con la desaparición de los dos rangos más
salinos, que alcanzan su máxima expresión precisamente en el año 1990, en el que llegan a
sumar conjuntamente casi el 50% de la superficie del acuífero. En resumen, se aprecia
claramente que en la década de los 80 la intensidad de la intrusión marina es más severa
como consecuencia de la expansión de los regadíos producida en esa época, coincidiendo
con la ocurrencia de una prolongada sequía que afectó prácticamente a todo el territorio
nacional. La notable disminución de las concentraciones medias de cloruros en el acuífero,
expresadas en la figura anterior como porcentajes de superficie afectada es consecuencia de
la progresiva reducción de los bombeos asociada a la modernización de los sistemas de
riego.
Estas variables deben tener como consecuencia que la evolución de la relación agua
dulce – agua salada haya sido diferente en cada sector del acuífero. Para comprobar si esto
es realmente así, se ha discriminado el acuífero en seis sectores (figura 19) de sur a norte
según la dirección preferencial de flujo.
87
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Sector Puntos
OT - 1 21
OT – 2 8, 9, 20, 22, 30
OT – 3 23, 25, 38, 39
OT – 4 4, 10, 15, 35, 36
OT – 5 11, 12, 13, 16, 32
OT - 6 6, 17, 18, 19, 31, 33, 34, 37
88
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
( xt − 2 / 2) + xt −1 + xt + xt +1 + ( xt + 2 / 2)
m( xt ) =
4
89
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
90
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
91
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
92
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
93
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
94
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
desde valores del orden de 300 mg/L hasta valores cercanos a 1800 mg/L. Si bien el único
punto de observación de este sector (nº 21) con registros históricos se encuentra situado
muy cerca de la costa, datos aislados (Morell, 1985; Giménez, 1994) apuntan a que el
incremento salino en esa época fue general en el sector, en el que existen numerosos pozos
pero con caudales unitarios de bombeo muy reducidos. La rápida elevación de la salinidad
es debida a que la recarga lateral en ese sector es prácticamente nula, ya que los materiales
de borde que están en contacto con el acuífero pliocuaternario son de muy baja
permeabilidad (niveles margosos del Aptiense inferior). No hay que descartar tampoco la
posibilidad de que, como ya apuntaba Morell en 1985, la circulación del agua sea muy
restringida y el tiempo de residencia muy elevado, lo que se evidencia por concentraciones
anómalamente elevadas de potasio, bromuro y estroncio. En los últimos años la salinidad
se ha reducido progresivamente, lo que se puede justificar por la disminución de los
bombeos a causa de la mala calidad del agua pero también por la fuerte presión urbanística
que soporta la zona.
En el resto del acuífero se puede admitir que existe un cierto patrón de comportamiento,
sectores 2, 3, 4, 5 y 6, que registran a grandes rasgos, incrementos de salinidad en los
periodos secos y descensos pronunciados en los periodos húmedos, lo que demuestra que
si bien el equilibrio es precario también existe una notable capacidad de recuperación; en el
caso del sector 3, además, puede jugar algún papel amortiguador de la intrusión la
presencia de la zona húmeda y en menor intensidad la explotación del acuífero. Menor
resiliencia muestran los sectores 4 y 6, en los que la alimentación lateral parece ser algo más
restringida, por lo que el comportamiento es ligeramente diferente al resto de sectores.
En consecuencia, se puede admitir que en los últimos años se aprecia una mejoría de la
salinidad del acuífero, lo que es coherente con la realidad de que los bombeos se han
reducido considerablemente.
2500
2000
Cloruros (mg/L)
1500
1000
500
0
jul-98 dic-99 abr-01 sep-02 ene-04 may-05 oct-06 feb-08 jul-09 nov-10 abr-12 ago-13
96
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
El índice SITE utiliza como datos de partida los valores de la concentración del ión
cloruro expresado en mg/L, en un número de puntos que debe ser suficiente para poder
elaborar un mapa de isocloruros con un adecuado grado de información.
El acrónimo SITE responde a los siguientes parámetros o criterios relacionados con el
proceso intrusivo: S = Superficie afectada; I = Intensidad del proceso; T = Temporalidad o
estacionalidad y E = Evolución a medio – largo plazo.
El valor de S es el cociente entre S>250 y Stotal, siendo S>250 la superficie ocupada por
contenidos de ion cloruro superiores a 250 mg/L y Stotal la superficie total del acuífero,
considerando una anchura máxima de 10 kilómetros desde la línea de costa (“superficie
afectable”).
Para cada rango de valores de S se asigna un valor, o índice numérico, y un código de
identificación, según el baremo de la tabla 6.
S Índice
Caracterización Código
(Superficie) numérico
0,00 – 0,10 Sin intrusión o Puntual 0 W
0,11 – 0,25 Localizada 1 L
0,26 – 0,50 Zonal o media 2 M
0,51 – 0,75 General o alta 3 H
0,76 – 1,00 Total o extrema 4 E
El número de puntos utilizados para la elaboración del índice SITE ha sido de 33 (tabla 7
y figura 27). Estos puntos cuentan con datos anuales de concentración de cloruros en el
periodo comprendido entre 1977 - 2001.
97
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
elabora el correspondiente mapa de isocloruros con las isolíneas de 150 mg/L, 250 mg/L,
500 mg/L, 1.000 y 1.500 mg/L (figura 27), y se determinan las superficies correspondientes
a los rangos definidos (tabla 8).
-
Nº Valor Medio Cl Matriz de Matriz de Matriz de
Código IGME Huso Sector Coorxutm Cooryutm Cota Prof_pozo
identificativo (1998-2001) cálculo S, I cálculo T cálculo E
302440005 1 31 T 256.348 4.442.074 4,34 8 412 x x x
302440009 2 31 T 256.199 4.442.436 10 6 454 x x x
302440013 3 31 T 255.610 4.441.752 18,51 30 655 x x x
302440016 4 31 T 255.963 4.446.301 39,12 43 x x
302440017 5 31 T 255.299 4.445.817 38 55 1.173 x x
302440020 6 31 T 255.742 4.446.036 33,79 51 1.302 x
302440027 7 31 T 256.390 4.446.641 33,71 38 1.524 x x x
312320003 8 31 T 266.269 4.460.442 78,44 184 440 x x x
312350001 9 31 T 259.323 4.451.798 2,09 250 503 x x
312350003 10 31 T 259.629 4.452.979 17,82 160 x x
312350004 11 31 T 263.635 4.457.468 44,72 60 903 x x
312350005 12 31 T 261.941 4.455.632 22 - 1.209 x x
312350006 13 31 T 263.371 4.456.960 38 41 355 x
312350007 14 31 T 263.429 4.456.730 35,5 42 763 x x
312350008 15 31 T 262.678 4.455.755 25,86 42 1.077 x x
312350013 16 31 T 259.691 4.453.803 28,47 34 375 x x x
312350017 17 31 T 262.118 4.453.987 4,73 12 976 x x x
312350022 18 31 T 260.955 4.453.257 10 17 921 x x x
312350024 19 31 T 264.209 4.454.088 7,21 19 1.205 x x x
312350028 20 31 T 258.567 4.451.551 28,64 36 301 x x
312350029 21 31 T 263.608 4.455.110 15,4 20 1.044 x x
312350036 22 31 T 261.302 4.456.521 35 45 213 x x
312360002 23 31 T 266.146 4.457.402 24,48 37 1.186 x x
312360006 24 31 T 264.630 4.457.346 34,81 41 1.304 x x x
312360007 25 31 T 266.587 4.454.997 12,9 16 1.175 x x x
312360009 26 31 T 266.914 4.456.002 9,68 11 1.272 x x
312360010 27 31 T 267.378 4.456.170 9,01 12 1.136 x x x
312410003 28 31 T 257.251 4.444.066 9,68 11 611 x x x
312410012 29 31 T 258.045 4.446.618 7,61 32 844 x x x
312410013 30 31 T 258.628 4.449.689 15,21 177 1.239 x
312410014 31 31 T 258.659 4.449.435 17 15 916 x x x
312410016 32 31 T 257.493 4.448.618 23,74 36 795 x x x
312410019 33 31 T 258.654 4.448.030 2,17 4 828 x x x
S >250 86,57
S= = = 0,96
S total 89,72
El resultado (0,96) indica que el acuífero está sometido a una intrusión de carácter total o
extrema, y se le asigna un valor de 4 y el código E.
98
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Figura 27. Red de control de la intrusión marina (izquierda) y mapa del valor medio de
isocloruros para el periodo 1998 -2001 (derecha)
Tabla 8. Distribución de superficies de los rangos de cloruros diferenciados (periodo 1998 - 2001)
99
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
I=
∑S i ⋅ Cl i
S total
I Índice
I Caracterización Código
(Vr=250 mg/L) numérico
<Vr < 250 Nula 0 W
Vr-2Vr 251 - 500 Baja 1 L
2Vr-4Vr 501 -1000 Media 2 M
4Vr-8Vr 1001 - 2000 Alta 3 H
>8Vr > 2.000 Extrema 4 E
Los rangos de cloruros y sus correspondientes valores medios, así como la superficie
ocupada por cada rango y los sumatorios calculados para este acuífero se presentan en la
tabla 10.
Rango Cl- Valor medio
Sup. (km2) Cl * Sup.
(mg/L) del rango
< 150 100 0 0
150 - 250 200 3,15 630
250 - 500 375 15,01 5629
500 – 1000 750 37,82 28368
1000 – 1500 1250 33,71 42138
1500 - 2000 1750 0,03 53
TOTAL 76817
Tabla 10. Distribución ponderada de la concentración de cloruros para el periodo 1998 - 2001
100
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
I=
∑S i ⋅ Cl i
=
76.817
= 856
S total 89,72
Con estos datos el valor de I queda establecido en 856, lo que indica que la salinidad del
acuífero es media y se le asigna en el índice numérico el valor 2 y el código de
identificación M.
Para el cálculo de este parámetro se dispone de dos datos anuales, con lo que se
obtienen dos series, una para los valores máximos y otra para los valores mínimos.
n
1
T= ∑ f max − f min
n x=1
donde f max corresponde a los valores máximos anuales y f min a los valores mínimos
anuales de cada año.
En este caso se seleccionaron un total de 19 puntos, que poseían dos medidas anuales de
la concentración de cloruros en una serie temporal de 20 años (tabla 12). Para cada punto se
ha calculado su valor máximo y mínimo por año y se han obtenido las correspondientes
medias máximas y mínimas anuales para el acuífero (figura 28, tabla 12).
1 n
T= ∑ f ( x) max − f ( x) min = 220
n x =1
El valor de T resultante es igual a 220, lo que indica una situación caracterizada como de
estacionalidad extrema, que corresponde al valor 4 y al código E.
101
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
T Índice
T Caracterización Código
(Vr=250 mg/L) numérico
< 0,1Vr < 25 Sin estacionalidad 0 W
0,1Vr-0,2Vr 26 – 50 Estacionalidad baja 1 L
0,2Vr-0,4Vr 51 – 100 Estacionalidad media 2 M
Estacionalidad
0,4Vr-0,8Vr 101 - 200 3 H
elevada
Estacionalidad
> 0,8Vr > 200 4 E
extrema
1400
1300
1200
1100
1000
900
V alor mínimo_medio
700 Tx
600
500
400
300
200
100
0
1977 1978 1979 1980 1981 1982 1983 1985 1987 1988 1989 1990 1992 1993 1994 1995 1996 1998 1999 2000
Años
102
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Cl- (mg/L)
Año Valor mínimo Valor máximo medio Tx
medio
1977 665 791 125
1978 630 854 224
1979 764 1048 289
1980 752 1136 383
1981 685 1049 364
1982 744 1109 365
1983 852 1051 196
1985 845 1315 470
1987 828 1190 362
1988 760 1055 295
1989 887 1074 187
1990 845 1016 171
1992 833 912 79
1993 739 893 155
1994 771 891 121
1995 858 994 135
1996 864 975 111
1998 864 1001 136
1999 849 954 105
2000 833 966 133
Media 793 1014 220
103
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
Cln
E=
Cln −1
años anteriores).
Índice
Valor de E Tipología Código
numérico
< 0,85 Recuperación importante -2 -2
0,85 – 0,97 Recuperación moderada -1 -1
0,98 – 1,02 Estable 0 0
1,03 – 1,15 Deterioro moderado 1 1
> 1,15 Deterioro importante 2 2
El parámetro E se calcula a partir del cociente entre el valor medio del contenido en
cloruros para el periodo 1998 - 2001 y para el periodo 1994 - 1997. Así, se obtiene un valor
de 1,01, que caracteriza una situación estable y se le asigna el valor numérico y el código de
0.
Cln
E= = 895/882= 1,01
Cln −1
Índice
Parámetro Rango Caracterización Código
numérico
0,00 – 0,10 Sin intrusión o puntual 0 W
0,11 – 0,25 Localizada o baja 1 L
S 0,26 – 0,50 Zonal o media 2 M
0,51 – 0,75 General o alta 3 H
0,76 – 1,00 Total o extrema 4 E
<Vr < 250 Nula 0 W
I Vr-2Vr 251 - 500 Baja 1 L
(Vr=250 2Vr-4Vr 501-1.000 Media 2 M
mg/L) 4Vr-8Vr 1.001-2.000 Alta 3 H
> 8Vr > 2.000 Extrema 4 E
< 0,1Vr < 25 Sin estacionalidad 0 W
T 0,1Vr-0,2Vr 26 - 50 Estacionalidad baja 1 L
(Vr=250 0,2Vr-0,4Vr 51-100 Estacionalidad media 2 M
mg/L) 0,4Vr-0,8Vr 101-200 Estacionalidad elevada 3 H
> 0,8Vr > 200 Estacionalidad extrema 4 E
< 0,85 Recuperación importante -2 -2
0,85 – 0,97 Recuperación moderada -1 -1
E 0,98 – 1,02 Estable 0 0
1,03 – 1,15 Deterioro moderado 1 1
> 1,15 Deterioro importante 2 2
Tabla 14. Síntesis de valores y códigos establecidos para cada uno de los criterios del índice SITE
Como resultado, el índice SITE puede ser representado mediante dos formas, una
alfanumérica, que da una información específica de cada uno de los parámetros que
intervienen en el proceso de salinización, y otra exclusivamente numérica, cuya resultante
es una única cifra. Según el fin que se quiera dar al índice, podrá ser utilizada una forma u
otra de expresión.
105
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
3S + 3I + T + E
ValorSITE =
30
El cociente 30 es el máximo valor posible del algoritmo del numerador, por lo que un
cierto acuífero totalmente salinizado tendría un valor SITE de 1 mientras que un acuífero
sin afección alguna tendría un valor SITE de 0. Dicho valor permite establecer categorías
del estado de la intrusión, o de presión debida a intrusión marina (tabla 15).
Por último el Índice SITE general es el resultado de unir el código alfanumérico y el valor
numérico del índice.
Los valores numéricos y códigos del índice SITE que caracterizan el acuífero de la Plana
de Oropesa - Torreblanca se incluyen en la tabla 16 así como la agregación de los cuatro
códigos obtenidos dan lugar al índice alfanumérico EME0. Por su parte, de la aplicación del
algoritmo de cálculo se obtiene el valor SITE de 0,73.
(3 * 4) + (3 * 2) + 4 + 0
ValorSITE = = 0,73
30
106
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
La expresión final del índice para la Plana de Oropesa - Torreblanca sería, por tanto,
EME0 (0,73), que caracteriza a un acuífero con intrusión marina total, de intensidad media
y variaciones estacionales muy significativas y en el que se ha registrado una evolución del
ión cloruro caracterizada como estable. La valoración general indica que la presión debida
a la intrusión marina en el acuífero es ALTA.
3.4. Conclusiones
El estudio de la evolución de los cloruros en la Plana de Oropesa – Torreblanca desde
1980 hasta la actualidad pone de manifiesto que la intrusión marina a lo largo de las
últimas décadas se puede dividir en cuatro fases. La primera fase, comprendida entre 1980
y 1985, se caracteriza por el avance del frente salino. En este periodo más del 50% de la
superficie de la Plana presenta concentraciones de cloruros comprendidas entre 1000 y 2000
mg/L. Durante la segunda fase (1986 – 1991) la intrusión marina llega a valores máximos,
observándose zonas con concentraciones de cloruros superiores a 2000 mg/L que en un
determinado instante llegan a representan el 20% de la superficie total. En la tercera fase
(1992 – 2000) se produce una leve recuperación del acuífero desapareciendo la superficie
afectada por cloruros de más de 2000 mg/L y la superficie afectada por cloruros de
concentraciones comprendidas entre 1000 y 2000 mg/L se reduce aproximadamente un
20%. La última fase caracterizada es la comprendida entre 2001 y la actualidad. Los datos
existentes son escasos para obtener mapas de distribución de concentraciones de cloruros
pero si se asume que dichos datos son representativos de la zona se puede observar que la
intrusión marina está en proceso de retroceso, lo que es coherente con la realidad
socioeconómica actual de la Plana de Oropesa – Torreblanca.
Para completar el estudio de la intrusión marina se ha considerado oportuno determinar
el índice SITE de la Plana de Oropesa – Torreblanca. Dicho índice va mas allá del estudio
de la distribución de los cloruros en un momento determinado e informa sobre la
intensidad, la estacionalidad y la evolución a medio plazo del proceso; además, permite
comparar diferentes acuíferos intruidos e incluso establecer prioridades de actuación en el
caso de que sean necesarias. El problema para determinar el índice actual es la falta de
datos (cloruros) en los últimos años, por lo que el índice SITE calculado es el característico
de finales de los años 90, que es el periodo más próximo a la actualidad con una matriz de
datos adecuada para determinarlo.
107
Capítulo 3 La intrusión marina en la Plana de Oropesa-Torreblanca
108
Capítulo 4 Modelo matemático de flujo de densidad constante
Capítulo 4
Modelo matemático de flujo
de densidad constante
109
Capítulo 4 Modelo matemático de flujo de densidad constante
110
Capítulo 4 Modelo matemático de flujo de densidad constante
∆h
Q = q ·S = k · ·S (1)
L
donde:
q: velocidad de Darcy o flujo especifico [LT-1];
S: sección perpendicular al flujo [L];
L: distancia entre los potenciales (∆h) [L];
∆h
: gradiente hidráulico (i);
L
k: permeabilidad o conductividad hidráulica [LT-1].
∂h ∂h ∂h
q x = − k xx ; q y = − k yy ; q z = − k zz (2)
∂x ∂y ∂z
111
Capítulo 4 Modelo matemático de flujo de densidad constante
∂ ∂h ∂ ∂h ∂ ∂h ∂h
k xx + k yy + k zz + Ws = Ss (3)
∂x ∂x ∂y ∂y ∂z ∂z ∂t
siendo:
Ws: flujo volumétrico por unidad de volumen y por unidad de tiempo.
Representa una fuente/sumidero [T-1];
Ss: coeficiente de almacenamiento específico del medio [L-1].
∂ ∂h ∂ ∂h ∂ ∂h
k xx + k yy + k zz + Ws = 0 (4)
∂x ∂x ∂y ∂y ∂z ∂z
La ecuación del flujo tiene infinitas soluciones. Para obtener una solución concreta se
requiere especificar las condiciones iniciales y las condiciones de contorno del problema.
Las condiciones iniciales fijan el estado del sistema en el instante inicial, y vienen
definidas por el nivel piezométrico en dicho instante inicial.
Las condiciones de contorno determinan la interacción del flujo con el medio que lo rodea.
Son declaraciones matemáticas, donde se especifica la variable dependiente (potencial
hidráulico) o la derivada de la variable dependiente (el flujo). Se clasifican en tres
categorías:
La discretización espacial del acuífero se lleva a cabo mediante una red rectangular
formada por paralelepípedos (3D), que es el resultado de combinar columnas y filas. El
método de diferencias finitas permite dos tipos de mallado (Anderson y Woessner, 1992):
113
Capítulo 4 Modelo matemático de flujo de densidad constante
1) Nodo centrado en los vértices de unión de las columnas y las filas. En este caso
cada nodo será compartido por las celdas contiguas, como por ejemplo en el
modelo desarrollado por Prickett and Lonnquist (1971).
Este código aproxima la ecuación del flujo (ecuación 3) en derivadas parciales mediante
diferencias finitas centrado en bloques.
Así, para una celda genérica, como la representada en la figura 29, las relaciones
matemáticas que vinculan los caudales existentes entre la celda unidad (i,j,k) y las 6 caras
que la rodean vienen plasmadas en las ecuaciones 1 y 6, donde se expresa el flujo que
discurre a través de la cara de las celdas (i,j-1,k) e (i,j,k) para un régimen estacionario y
unidireccional. En cambio las “entradas” o “salidas” del agua en la celda unidad pueden
expresarse mediante la ecuación 7.
114
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 29. Esquema del flujo subterráneo en una celda genérica (i,j,k) proveniente de
la celda adyacente (i, j-1,k) (modificado de McDonald y Harbaugh, 1988)
∆h
q = kS = C∆h (Ley de Darcy); (1)
L
(h − hi , j ,k )
(∆ci ∆vk ) × (hi , j −1,k − hi , j ,k ) ;
i , j −1, k
q 1 = kR 1 = CR 1 (6)
i , j − ,k
2
i , j − ,k
2
∆r 1
i , j − ,k
2
j−
2
donde:
q 1 : caudal que circula a través de la cara de las celdas (i,j-1,k) e (i,j,k) [L3T-1];
i , j − ,k
2
kR 1 : conductividad hidráulica a lo largo de una fila entre los nodos (i,j,k) e (i,j-1,k)
i , j − ,k
2
[LT-1];
CR 1 : conductancia a lo largo de una fila entre los nodos (i,j,k) e (i,j-1,k) [L2T-1];
i , j − ,k
2
115
Capítulo 4 Modelo matemático de flujo de densidad constante
donde:
I ext : caudal que proviene del exterior y no depende del nivel piezométrico[L3T-1];
∆hi , j , k
q 1 +q 1 +q 1 +q 1 +q 1 +q 1 + Pi , j , k hi , j , k + I i , j , k = Ssi , j , k (∆rj ∆ci ∆vk )
i, j − ,k
2
i, j + ,k
2
i − , j ,k
2
i + , j ,k
2
i, j ,k −
2
i , j ,k +
2
∆t
donde:
q 1 : caudal que circula a través de la cara de las celdas (i,j-1,k) e (i,j,k) [L3T-1];
i , j − ,k
2
q 1 : caudal que circula a través de la cara de las celdas (i,j+1,k) e (i,j,k) [L3T-1];
i , j + ,k
2
q 1 : caudal que circula a través de la cara de las celdas (i-1,j,k) e (i,j,k) [L3T-1];
i − , j ,k
2
q 1 : caudal que circula a través de la cara de las celdas (i+1,j,k) e (i,j,k) [L3T-1];
i + , j ,k
2
q 1 : caudal que circula a través de la cara de las celdas (i,j,k-1) e (i,j,k) [L3T-1];
i , j ,k −
2
q 1 : caudal que circula a través de la cara de las celdas (i,j,k+1) e (i,j,k) [L3T-1];
i , j ,k +
2
I i , j , k : caudal que proviene del exterior y no depende del nivel piezométrico [L3T-1];
116
Capítulo 4 Modelo matemático de flujo de densidad constante
∆hi , j ,k
: Aproximación en diferencias finitas de la derivada parcial del nivel
∆t
∂hi , j , k
piezométrico respecto al tiempo [LT-1];
∂t
117
Capítulo 4 Modelo matemático de flujo de densidad constante
La base del acuífero lo constituyen dos unidades bien diferenciadas. Una de carácter
impermeable compuesta por margas y arcillas del Mioceno, y otra, constituida por calizas
del Cretácico inferior (Barremiense, Aptiense y Albiense) y Cretácico superior
(Cenomaniense), que son representados en la capa 2 del modelo. Esta última presenta
carácter acuífero por lo que se ha considerado oportuno representarla como una capa
activa que recibe parte de los aportes provenientes de los acuíferos limítrofes y
posteriormente alimentará al acuífero principal (figura 30).
118
Capítulo 4 Modelo matemático de flujo de densidad constante
119
Capítulo 4 Modelo matemático de flujo de densidad constante
120
Capítulo 4 Modelo matemático de flujo de densidad constante
Discretización espacial
121
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 31. Imagen en planta y cortes transversales del modelo, filas y columnas
(Visual Modlfow)
122
Capítulo 4 Modelo matemático de flujo de densidad constante
Discretización temporal
123
Capítulo 4 Modelo matemático de flujo de densidad constante
124
Capítulo 4 Modelo matemático de flujo de densidad constante
Coeficiente de almacenamiento
125
Capítulo 4 Modelo matemático de flujo de densidad constante
126
Capítulo 4 Modelo matemático de flujo de densidad constante
Nivel
Nomenclatura
PUNTO UTM X UTM Y Piezométrico
CHJ
(msnm)
4 08.11.004 770.109.00 4.453.787.00 0,98
6 08.11.006 776.623.00 4.458.829.00 1,85
127
Capítulo 4 Modelo matemático de flujo de densidad constante
Nivel
PUNTO UTM X UTM Y Piezométrico
(msnm)
42 767.016.00 4.451.077.00 12,5
128
Capítulo 4 Modelo matemático de flujo de densidad constante
Potencial constante
Esta condición de contorno se aplica tanto en la costa como en el resto de los
límites laterales subterráneos a través de los cuales se llevan a cabo las
entradas más importantes al acuífero.
En la línea de costa se ha establecido un potencial constante de 0 msnm, que
no variará durante la simulación.
En el resto de contactos laterales subterráneos asociados a la capa 2, en las
celdas que se va a imponer esta condición de contorno, se ha relacionado el
potencial constante (que a priori es desconocido) con las precipitaciones. Para
facilitar este proceso las precipitaciones se han dividido en tres grupos y se han
diferenciado tres tipos de potencial constante.
Grupos correspondientes a las precipitaciones son:
- Grupo 1: Precipitaciones muy abundantes, mayores a 800 mm/año
- Grupo 2: Precipitaciones abundantes, de 600 a 800 mm/año
- Grupo 3: Precipitaciones normales menores a 600 mm/año.
129
Capítulo 4 Modelo matemático de flujo de densidad constante
Drenes
donde:
131
Capítulo 4 Modelo matemático de flujo de densidad constante
4.3.4. Recarga
Los dos componentes que definen la recarga son: la infiltración de las precipitaciones y
los retornos de riego.
132
Capítulo 4 Modelo matemático de flujo de densidad constante
133
Capítulo 4 Modelo matemático de flujo de densidad constante
134
Capítulo 4 Modelo matemático de flujo de densidad constante
4.3.5. Bombeos
Los bombeos se han obtenido a partir de distintas fuentes: las administraciones públicas
(Ayuntamiento de Torreblanca, 2009), las comunidades de regantes e informes técnicos
(Anónimo, 1987; Morell, 2008).
Otro dato que en algunos casos es desconocido es el caudal de extracción entre los años
setenta y ochenta. En este caso se ha optado por aumentar de forma escalonada los
bombeos hasta alcanzar el caudal máximo (conocido) en los años noventa y
posteriormente, hasta alcanzar el año 2009, se han disminuido del mismo modo el caudal
que se ha considerado (incluso hasta su desaparición en algunos casos), en función de la
superficie de regadío existente (figura 42).
135
Capítulo 4 Modelo matemático de flujo de densidad constante
136
Capítulo 4 Modelo matemático de flujo de densidad constante
BOMBEOS
40.00
35.00
30.00
Q(hm3/año)
25.00
20.00
15.00
10.00
5.00
0.00
1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010
TIEMPO(años)
BOMBEOS
137
Capítulo 4 Modelo matemático de flujo de densidad constante
138
Capítulo 4 Modelo matemático de flujo de densidad constante
Tabla 22. Pozos ubicados en la Plana de Oropesa - Torreblanca II. Pozos calculados
Tabla 23. Pozos ubicados en la Plana de Oropesa - Torreblanca III. Pozos de abastecimiento
urbano
139
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 43. Situación de los puntos de observación y sus respectivas coordenadas UTM
140
Capítulo 4 Modelo matemático de flujo de densidad constante
NO SI
CALIBRADOS CALIBRADOS
Conductividad hidráulica X
Almacenamiento X
Nivel Piezométrico inicial X
Potencial Cte. interior (TLS) X
Potencial Cte. (mar) X
Conductancia(drenes) X
Recarga X
Pozos X
141
Capítulo 4 Modelo matemático de flujo de densidad constante
En el sector Norte (figura 44) cabe destacar el punto 6, que se caracteriza por presentar
datos de observación piezométrica por debajo del nivel del mar durante la segunda mitad
del periodo de simulación; los datos reales de observación piezométrica presentan un
máximo próximo a los 4 msnm y un mínimo cercano a los 4 mbnm en un corto periodo de
tiempo. Seguramente el causante de este salto piezométrico sea la ubicación de un pozo
activo (Pou Anna) en los alrededores del mismo.
142
Capítulo 4 Modelo matemático de flujo de densidad constante
En el sector Cabanes - Costa (figura 46) se observa como el modelo también capta en
general las tendencias y el valor de las observaciones de campo.
143
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 46. Calibración de los puntos de observación del sector Cabanes - Costa
144
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 47. Calibración de los puntos de observación del sector Cabanes - Interior
145
Capítulo 4 Modelo matemático de flujo de densidad constante
146
Capítulo 4 Modelo matemático de flujo de densidad constante
opta por elegir un periodo en el cual la intrusión marina presente una cierta situación de
equilibrio y las salidas al mar sean aproximadamente constantes y suficientes para impedir
el progreso de la misma. En dicho informe se ha seleccionado el periodo 1998 - 2001 por
cumplir los requisitos mencionados, por lo que éste será el periodo utilizado para llevar a
cabo la verificación del modelo. Del periodo 1998 - 2001 se ha calculado la media de: la
recarga, de los bombeos, de los puntos de observación piezométrica y del potencial
constante, para conseguir un dato representativo y así poder ejecutar el modelo en régimen
estacionario. El resto de partidas que conforman el modelo se mantienen con los mismos
valores que el modelo ejecutado en régimen transitorio.
Se han diferenciado cuatro zonas de recarga (figura 49 y tabla 25). La zona 1 (color
morado) se caracteriza por recibir tanto infiltración de agua de lluvia como de
retorno de riego. En las celdas correspondientes a la zona 2 (verde) se infiltra
únicamente el agua procedente de las precipitaciones. Las celdas correspondientes
a los núcleos urbanos de Oropesa y Torreblanca (zona 4) se consideran
impermeabilizadas, por ello la recarga será nula. A las celdas situadas en la línea
de costa y en el límite interior también se les ha aplicado la recarga de 0 mm/año
(zona 4). En este último caso, la infiltración no debe influir en el balance final
calculado por el programa, por lo que la infiltración debe ser nula. Existen otras
zonas con un grado de urbanización menor, aproximadamente el 50 % de la celda
afectada, por lo que se les aplica una recarga que corresponde a la mitad de la
infiltración de la precipitación (zona 3).
147
Capítulo 4 Modelo matemático de flujo de densidad constante
Retorno de
Precipitación TOTAL
Zona Color riego
(mm/año)* (mm/año)
(mm/año)*
1 morado 20 66,4 86,4
2 verde 0 66,4 66,4
3 Gris claro 0 33,2 33,2
4 rojo 0 0 0
Tabla 25. Recarga del modelo estacionario.*Datos resultantes de aplicar la tasa de infiltración
correspondiente a la dotación de riego (8%) y a la precipitación(14%) (Tuñón, 2000)
148
Capítulo 4 Modelo matemático de flujo de densidad constante
Potencial constante
Bombeos
El caudal que extrae cada bombeo viene determinado por el cálculo de la media
de las extracciones en cada pozo. En la figura 50 se muestra la ubicación de los
bombeos.
149
Capítulo 4 Modelo matemático de flujo de densidad constante
Los datos piezométricos se han obtenido a partir del cálculo de la media de los
valores existentes en cada punto de observación en el periodo correspondiente. En
la tabla 26 se muestran las coordenadas y el nivel piezométrico medio de cada
punto y la figura 51 muestra la ubicación en la Plana de dichos puntos.
UTM
N.P.medio (m)
Punto X Y
6 776513 4458619 0.42
7 767215 4442935 -0.47
8 766892 4446594 0.14
9 767498 4447242 0.71
11 773157 4456761 0.73
12 772718 4454959 0.27
13 774908 4455411 0.01
14 772219 4455676 0.57
17 777109 4456269 0.53
18 777368 4457294 0.27
19 777819 4457493 0.37
21 768531 4444730 0.08
22 769151 4447330 0.19
23 769525 4450434 0.05
24 770071 4449729 0.57
25 769663 4448781 0.44
26 769655 4447728 0.09
27 769755 4447777 0.34
28 769752 4447327 0.55
29 767110 4442234 -0.96
150
Capítulo 4 Modelo matemático de flujo de densidad constante
4.5.2. Resultados
En la tabla 27 se observa la relación entre los datos piezométricos reales y los calculados
por el programa. Se considera que el modelo estacionario es válido ya que la diferencia
entre los resultados de la simulación y las observaciones está comprendida entre 1 y -1.
Este rango se ha estimado adecuado debido a que algunos de los datos introducidos en la
simulación son medias de los datos reales. En la tabla 27 se observa que las diferencias
entre niveles piezométricos observados y los calculados se sitúan dentro del rango
anteriormente definido como válido (1 y -1), el único punto que se sitúa fuera del rango es
el 23. Desde un punto de vista global, se considera que el modelo está bien calibrado. En la
figura 52 se muestran gráficamente los resultados obtenidos
151
Capítulo 4 Modelo matemático de flujo de densidad constante
Tabla 27. Calibración del modelo ejecutado en régimen estacionario. Datos reales (Obs.) vs. datos
calculados (Calc.).
152
Capítulo 4 Modelo matemático de flujo de densidad constante
Es bastante Es poco
sensible sensible
K (capa 1) X
K (capa 2) X
PARÁMETROS Ss (capa 2) X
HIDRODINÁMICOS Sy (capa 1) X
Conductancia
(drenes)
X
Recarga X
ACCIONES Bombeos X
TLS X
Tabla 28. Parámetros y acciones a los cuales se les ha realizado el análisis de sensibilidad
La recarga resultó ser poco sensible a pequeñas variaciones, por lo que se omitirán los
resultados obtenidos para la misma en este análisis de sensibilidad. Sin embargo, las
simulaciones llevadas a cabo (para valorar el escenario futuro y determinar el VER) a
puesto de manifiesto que la intrusión marina si que es sensible a los bombeos activos.
A continuación se describen detalladamente los resultados de las simulaciones llevadas
a cabo para valorar la sensibilidad del modelo al resto de parámetros y variables
contempladas.
153
Capítulo 4 Modelo matemático de flujo de densidad constante
154
Capítulo 4 Modelo matemático de flujo de densidad constante
155
Capítulo 4 Modelo matemático de flujo de densidad constante
156
Capítulo 4 Modelo matemático de flujo de densidad constante
157
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 57. Sensibilidad del punto de observación 6 a las TLS simuladas en régimen
transitorio
Figura 58. Sensibilidad del punto de observación 9 a las TLS simuladas en régimen
transitorio
158
Capítulo 4 Modelo matemático de flujo de densidad constante
Figura 59. Sensibilidad del punto de observación 12 a las TLS simuladas en régimen
transitorio.
Figura 60. Sensibilidad del punto de observación 23 a las TLS simuladas en régimen
transitorio
159
Capítulo 4 Modelo matemático de flujo de densidad constante
40.00 1200
35.00
1000
Precipitaciones (mm/año)
Salidas al mar y bombeos
30.00
800
25.00
(hm3/año)
20.00 600
15.00
400
10.00
200
5.00
0.00 0
1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012
TIEMPO(años)
160
Capítulo 4 Modelo matemático de flujo de densidad constante
16.00
14.00
12.00
10.00
Q (hm 3/año)
8.00
6.00
4.00
2.00
0.00
-2.00
1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012
Tiempo (años)
39.00
34.00
29.00
24.00
Q hm3 /año)
19.00
14.00
9.00
4.00
-1.00
1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012
Tiempo (años)
161
Capítulo 4 Modelo matemático de flujo de densidad constante
De la imagen obtenida con el Visual Modflow (figura 64), que muestra las isopiezas en
el periodo de mayor afección de la intrusión marina, se deduce que este proceso no se
reparte homogeneamente por toda la superficie estudiada, sino que es el sector Sur
(Cabanes y Oropesa) el más afectado. En este sector, por tanto las relaciones de la intrusión
marina con la recarga, los bombeos y las salidas al mar (figura 65 y 66) se identifican más
claramente.
162
Capítulo 4 Modelo matemático de flujo de densidad constante
13.00
11.00
9.00
Q(hm 3/año)
7.00
5.00
3.00
1.00
-1.00
1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012
Tiempo (año)
Figura 65. Intrusión, recarga y bombeos del sector Sur (1973 – 2009)
4.00
3.50
3.00
2.50
Q (hm3 /año)
2.00
1.50
1.00
0.50
0.00
-0.50
1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012
Tiempo (años)
163
Capítulo 4 Modelo matemático de flujo de densidad constante
164
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Capítulo 5
Volumen ecológico de remediación (VER)
de la intrusión marina
165
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
166
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Estos caudales ecológicos se refieren al mantenimiento de los ríos (UE - DMA, 2000)
pero dicha definición se podría extrapolar y aplicar a las aguas subterráneas, para
mantener el equilibrio de los acuíferos y preservar las zonas húmedas. Sin embargo, los
condicionantes que hay que tener en cuenta en cada caso son muy diferentes y la
metodología para su determinación también es distinta.
En la práctica son varios los vocablos que se utilizan indistintamente para definir el
mismo concepto: caudales ecológicos, caudales medioambientales, caudales de
mantenimiento, caudales de compensación, caudales de reserva, caudales reservados, etc.
En el caso de los acuíferos costeros, además, se utilizan otros conceptos específicos como
caudales de descarga al mar o volúmenes ecológicos. Si se analiza detenidamente el
significado de cada una de estas expresiones no es difícil apreciar diferencias conceptuales
notables que deberían ser tenidas en cuenta con mayor rigor para que, en cada caso, se
apliquen con mayor efectividad y precisión.
Un caudal circulante por un cauce podría ser considerado como ecológico, siempre que
fuese capaz de mantener el funcionamiento, composición y estructura del ecosistema
fluvial que ese cauce contiene en condiciones naturales. Es obvio que los mejores caudales
ecológicos serán los que imiten el régimen natural, ya que las biocenosis acuáticas
evolucionan de acuerdo con las pautas históricas de avenidas y estiajes. Sin embargo, tal
como se apunta en el Libro Blanco del Agua (MIMAN, 2000), la necesidad de aprovechar el
agua para los usos solicitados por la sociedad obliga a considerar como restricciones del
sistema no esos caudales óptimos, sino otros requerimientos mínimos que mantengan las
167
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
poblaciones naturales del río y sus valores ecológicos, de tal modo que no puedan
experimentar una disminución de su cuantía sin que ello implique una pérdida marcada de
los mismos. Es evidente que existe una amplia gama de caudales circulantes que son
ecológicos para un determinado cauce. Se necesita, por tanto, definir un Régimen de
Caudales Ecológicos con el fin de proteger el hábitat fluvial en todos los estados de
desarrollo de las especies acuáticas existentes en el mismo. Las distintas metodologías de
cálculo desarrolladas para la determinación de estos caudales ecológicos mínimos tienen
como objetivo la cuantificación y distribución temporal de unos caudales que permitan
atender los requerimientos ambientales dentro del cauce. Sin entrar en detallar los métodos
más utilizados, todos ellos se basan en el análisis de los regímenes de caudales históricos y
en el análisis de la variación del hábitat con los caudales circulantes.
En el Libro Blanco del Agua (MIMAN, 2000), no se hace referencia a un concepto similar
al de caudal ecológico aplicable a masas de agua subterránea. Poco después, la Directiva
Marco del Agua (UE - DMA, 2000) limita la disponibilidad de las aguas subterráneas a los
recursos naturales renovables de los acuíferos una vez detraídas sus salidas necesarias a
efectos medioambientales, que deben ser, en todo caso, establecidas en los planes de
cuenca. Estas salidas se refieren al mantenimiento de caudales mínimos permanentes con
carácter ecológico, aguas abajo de los principales embalses, al establecimiento de un criterio
de flujo genérico para el resto de la red hidrográfica, a la determinación del volumen de
agua necesaria para mantener las zona húmedas y al establecimiento de los caudales
subterráneos que deben salir al mar para evitar la intrusión marina.
En la cuenca del Júcar, por ejemplo, el documento de síntesis para el seguimiento del
Plan Hidrológico de la Cuenca del Júcar (CHJ, 2007) define los volúmenes
medioambientales de cada masa de agua subterránea como la suma de los volúmenes
ecológicos (mantenimiento de caudales en cauces fluviales), las salidas al mar y zonas
húmedas. En un documento anterior (CHJ, 2004) se consideraba también como restricción
ambiental para cada unidad hidrogeológica un valor igual al 80% de las descargas laterales
(transferencias laterales subterráneas) a las unidades hidrogeológicas contiguas, ya que
estas descargas laterales constituyen, en ocasiones, una proporción significativa del recurso
disponible en esas unidades hidrogeológicas.
En el Plan Hidrológico del Júcar (CHJ, 2007), se especifica que “… puede asimilarse a la
problemática de los caudales mínimos la necesidad de mantener en una posición estable la interfaz de
los acuíferos costeros, para lo que es imprescindible el mantenimiento de unos drenajes hacia el mar.
En la actualidad, estas salidas son suficientes en la mayor parte del frente costero salvo en puntos
concretos; sin embargo, es de prever un mayor riesgo en el futuro al intensificarse la explotación de
168
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
los recursos hidráulicos. Aunque existe una evaluación bastante ajustada de cuáles son los
volúmenes necesarios para no plantear problemas en la calidad y disponibilidad de los recursos,
debería ahondarse en su conocimiento, así como en la determinación de las líneas preferentes de
descarga”.
169
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Las salidas subterráneas al mar constituyen una partida del balance hídrico de los
acuíferos costeros. Habitualmente, esta partida se calcula directamente mediante la simple
aplicación de la fórmula de Darcy, o indirectamente mediante un balance entre las entradas
y las salidas. El cálculo directo requiere el conocimiento de la geometría del acuífero en la
franja costera, que permita determinar el espesor saturado, la conductividad hidráulica y el
gradiente hidráulico en diferentes condiciones. La estimación indirecta a través del balance
hídrico requiere abundante información sobre todas las partidas que componen el balance
y sobre las variaciones en el almacenamiento del acuífero.
Parte de este problema puede resolverse si el balance se establece para cortos periodos
de tiempo, más o menos homogéneos, pero la notable inercia de muchos de los acuíferos
costeros puede distorsionar la relación causa - efecto entre los eventos de recarga, la
explotación y las salidas al mar.
La aplicación de la ley de Darcy para estimar los caudales de salida al mar también
presenta escollos de difícil superación. En primer lugar, no todo el espesor saturado es
activo para el flujo de agua dulce sino que tiene lugar una progresiva reducción de la
sección de flujo que depende de la situación y geometría de la interfase teórica. La
170
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
171
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
La segunda fase consiste en delimitar y disminuir las áreas afectadas por la intrusión
marina. El procedimiento es el siguiente: a partir del modelo base en régimen influenciado
(modelo RI), se identifican las zonas donde existe intrusión marina y los bombeos situados
en las mismas, que son los causantes directos de la depresión del nivel piezométrico por
debajo de 0 msnm en las cercanías de la costa. El siguiente paso es disminuir un 10% los
bombeos seleccionados y volver a ejecutar el modelo. El resultado obtenido es otra
situación en la cual la superficie afectada por la intrusión marina puede coincidir o no con
la del modelo RI. A continuación se vuelve a delimitar la zona intruida y los bombeos
causantes de este nuevo escenario se vuelven a disminuir un 10% y se ejecuta otra vez el
modelo obteniendo una nueva superficie afectada. Este proceso se repite tantas veces como
sea necesario hasta que la piezometría del modelo resultante:
- En el caso de intrusión de referencia nula, esté por encima del nivel del mar
172
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
un acuífero o sector del mismo se comparan los bombeos del modelo RI con los bombeos del
último escenario realizado. La diferencia entre ambos es el volumen ecológico de remediación.
173
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Para los tres años tipo se considera que permanecerán constates en el tiempo:
- Las transferencias laterales subterráneas (TLS) se simularán mediante una
condición de contorno de tipo I (Dirichlet) y el valor introducido cambiará en
función del año calculado.
-El volumen extraído mediante captaciones es el mismo en los tres modelos base
(modelo RI) propuestos.
-El volumen reducido de los bombeos causantes de la intrusión marina se va a
sustituir por agua de otra procedencia, por lo tanto el retorno de riego no será
modificado.
174
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 1:
Como se puede observar en la figura 68 en el modelo RNh (año húmedo) no existe
intrusión marina y el flujo es perpendicular a la costa. Por tanto pretendemos limitar la
intrusión a la ausencia total de la misma (intrusión de referencia).
175
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 2:
La figura 69 muestra el modelo RIh a partir del cual se aplican los descensos de los
bombeos correspondientes para forzar el retroceso de la intrusión marina y hacerla
desaparecer. En este caso la intrusión marina únicamente afecta a la zona sur, por lo tanto
será esa zona donde se aplicarán los descensos de los bombeos.
Figura 69. Modelo RIh utilizado para el cálculo de VERh. Zona afectada por la
intrusión y bombeos involucrados
176
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
177
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
178
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
179
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Figura 73. Situación y reducción (%) de los bombeos para el cálculo del VERh
Reducción Volumen
Volumen
POZO bombeo reducido
(m3/día)
(%) (m3/día)
180
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 3:
El volumen ecológico de remediación para un año húmedo (VERh) es de 1,0 hm3/año
(tabla 32).
BOMBEOS hm3/año
FASE 1:
En ausencia de bombeos (modelo RNm), en el caso del año medio, el modelo tampoco
presenta niveles por debajo del nivel del mar (intrusión de referencia) (figura 74).
181
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 2:
En el año medio los niveles piezométricos por debajo del nivel del mar, como se puede
observar en la figura 75, afectan únicamente al sector sur. En dicha figura los pozos
señalados con un círculo amarillo son aquellos a los que se aplica la reducción del caudal
para propiciar la recuperación del acuífero.
182
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Figura 75. Modelo RIm utilizado para el cálculo de VERm. Zona afectada por la
intrusión y bombeos involucrados
183
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Figura 76. Piezometría y dirección del flujo subterráneo del modelo RIm y del
escenario 7 (VERm)
Figura 77. Situación y reducción (%) de los bombeos para el cálculo de VERm.
184
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Volumen
Volumen Reducción
POZO reducido
(m3/día) bombeo
(m3/día)
1 530 90% 477
2 530 90% 477
3 530 90% 477
4 530 80% 424
5 530 80% 424
6 530 70% 371
7 530 70% 371
8 530 70% 371
9 530 60% 318
10 530 50% 265
11 530 40% 212
12 530 40% 212
13 530 20% 106
14 530 10% 53
TOTAL 1663670
(m3/año)
FASE 3:
El volumen ecológico de remediación para un año medio (VERm) es de 1,6 hm3/año
(tabla 34).
BOMBEOS hm3/año
185
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 1:
Como se puede observar en la figura 78 en régimen natural (modelo RNs) el flujo
presenta una dirección preferencial NNO - SSE, observándose en toda la superficie el nivel
piezométrico por encima del nivel del mar, dicha situación propicia la contención de la
intrusión marina.
186
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 2:
La intrusión marina, en este caso, afecta a dos zonas de la Plana de Oropesa –
Torreblanca: zona Torreblanca y zona Sur. En la figura 79 se muestra el modelo RIs, las
zonas afectadas por la intrusión marina y los bombeos respectivos involucrados
directamente en el proceso.
Figura 79. Modelo RIs utilizado para el cálculo del VERs. Zonas afectadas y bombeos
involucrados
187
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
188
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Reducción Volumen
Volumen
POZO bombeo reducido
(m3/día)
(%) (m3/día)
25 530 30% 159
26 530 20% 106
27 530 20% 106
28 530 30% 159
29 1589 30% 477
30 4206 30% 1262
31 1655 10% 166
TOTAL 888344
(m3/año)
189
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
Reducción Volumen
Volumen
POZO bombeo reducido
(m3/día)
(%) (m3/día)
1 530 90% 477
2 530 90% 477
3 530 90% 477
4 530 90% 477
5 530 90% 477
6 530 80% 424
7 530 80% 424
8 530 80% 424
9 530 70% 371
10 530 70% 371
11 530 70% 371
12 530 60% 318
13 530 50% 265
14 530 50% 265
15 530 30% 159
16 530 40% 212
17 530 20% 106
18 530 40% 212
19 530 20% 106
20 530 40% 212
21 530 20% 106
22 1313 10% 131
23 2648 20% 530
24 2330 10% 233
TOTAL 2783110
(m3/año)
190
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
FASE 3:
El volumen ecológico de remediación para un año seco (VERs) en la Plana de Oropesa -
Torreblanca es de 3,7 hm3/año, que se distribuye en 0,9 hm3/año para la zona de
Torreblanca y 2,8 hm3/año en la zona Sur (tabla 37).
5.5. Conclusiones
Se ha propuesto una metodología para calcular el volumen ecológico de remediación
(VER) de la intrusión marina, definido como el volumen en que hay que reducir los
bombeos en un acuífero o sector del mismo para conseguir que la recuperación
piezométrica alcanzada comporte una intrusión menor o igual a la estimada en régimen
natural para el sistema (intrusión de referencia). La metodología parte de un modelo
matemático de flujo de densidad constante calibrado en régimen estacionario. Consiste en
reducir de forma iterativa los bombeos involucrados en el proceso de intrusión, hasta que
la superficie piezométrica obtenida comportara una intrusión menor o igual a la estimada
en régimen natural para el sistema (intrusión de referencia). Se ha demostrado su
aplicabilidad en un caso de estudio, la Plana de Oropesa - Torreblanca (MASub 080.110).
Se ha analizado el nivel piezométrico en régimen natural en el acuífero para 3 años tipo:
año húmedo (modelo RNh), año medio (modelo RNm) y año seco (modelo RNs). En los tres
supuestos los niveles piezométricos están situados por encima del nivel del mar y por lo
tanto la intrusión puede ser considerada despreciable en condiciones naturales.
Sin embargo la actual operación del sistema (modelo RI) si que lleva asociado un
avance de la intrusión marina. Se calculó de VER para años húmedos (VERh), medios
(VERm) y secos (VERs). En los tres casos se requiere aplicar restricciones en los bombeos en
el sector Sur.
Como era de esperar, en un año seco el stress hídrico se acentúa, por lo que el VERs
(3,7 hm3) es el mayor de los tres calculados. En este caso la reducción de los bombeos para
que los niveles piezométricos estén por encima del nivel del mar (VERs) es del orden de
191
Capítulo 5 Volumen ecológico de remediación (VER) de la intrusión marina
3,5 veces mayor que en el caso menos restrictivo, el año húmedo (1,0 hm3). Otra
característica asociada a la intrusión del año seco es que la superficie afectada se reparte
entre dos sectores: sector Sur, situado en la parte más meridional, y el sector Torreblanca,
situado en la mitad septentrional de la Plana de Oropesa - Torreblanca, entre el núcleo
urbano de Torreblanca y la costa. El VERs no se reparte equitativamente por zonas, sino
que el correspondiente a la zona Torreblanca (0,9 hm3) es una tercera parte del volumen
que hay que reducir en la zona Sur (2,8 hm3).
192
Capítulo 6 Escenario futuro
Capítulo 6
Escenario futuro
193
Capítulo 6 Escenario futuro
194
Capítulo 6 Escenario futuro
6. ESCENARIO FUTURO
6.1. Cambios de usos del suelo. Marco legal y objetivos
Debido a que los proyectos urbanísticos a corto y medio plazo en la Plana de Oropesa –
Torreblanca prevén cambios en el uso del suelo y, por tanto, unas necesidades hídricas de
la zona diferentes a las actuales. Resulta interesante valorar los impactos que dichas
actuaciones podrían producir sobre el sistema. Para ello se han llevado a cabo simulaciones
del sistema para un escenario futuro.
En el término municipal de Torreblanca está aprobado el Plan General de Ordenación
Urbana 2009 (PGOU Torreblanca, 2009), en virtud del cual el 70% de la superficie, que en la
actualidad está dedicada a la agricultura de cítricos, se ha transformado en superficie
catalogada como residencial urbanizable e industrial. En la zona costera situada al norte de
Prat de Cabanes, está proyectado el desarrollo urbanístico del Programa de Actuación
Integrado (PAI) denominado Doña Blanca Golf (figura 82).
En el interior de la Plana de Oropesa - Torreblanca, en los términos municipales de
Cabanes y Oropesa, está aprobado el PAI Marina d´Or Golf, que incluye tres campos de
golf, urbanizaciones privadas, complejos hoteleros y zonas ajardinadas asociadas. Su
extensión, una vez terminadas todas las acciones allí proyectadas, será de
aproximadamente 16 Km2 (figura 82). Una de las condiciones que el Gobierno Valenciano
impuso para la aprobación del PAI Marina D´Or Golf es que la totalidad del agua necesaria
para riego de los campos de golf y las zonas ajardinadas ha de ser agua residual
regenerada, y el agua destinada al abastecimiento urbano ha de proceder de una planta
desaladora. Por tanto a medida que se vayan ejecutando los proyectos diseñados, los
bombeos existentes irán reduciéndose llegando al cese total de las extracciones al terminar
el periodo de desarrollo del PAI.
El inicio de las obras de estos proyectos están sufriendo retraso debido a la situación
económica actual pero se espera que en un plazo razonable de tiempo puedan comenzar.
El objetivo principal de esta simulación para el escenario futuro, es conocer si los
cambios de los usos del suelo y la disminución de los bombeos puede modificar la
situación hidrodinámica actual.
195
Capítulo 6 Escenario futuro
196
Capítulo 6 Escenario futuro
Tabla 38. Periodos de construcción del campo de golf de Doña Blanca Golf, superficie urbanizada
asociada y sectores afectados por el PGOU Torreblanca, 2009. (*Superficie urbanizada relativa a las
urbanizaciones asociadas a los campos de Golf)
Campo de Superficie
Años
Golf urbanizada*
0-2 1º 0
3-5 2º 25%
6-10 3º 50%
11-15 --- 75%
16-20 --- 100%
Tabla 39. Periodos de construcción de los campos de golf de Marina d´Or Golf y superficie
urbanizada asociada. (*Superficie urbanizada relativa a las urbanizaciones asociadas a los campos de
Golf)
197
Capítulo 6 Escenario futuro
Figura 83. Sectores de actuación y Doña Blanca Golf. (PGOU Torreblanca, 2009)
198
Capítulo 6 Escenario futuro
Figura 84. Serie histórica anual de las precipitaciones. Periodo representativo de las
precipitaciones 1999 - 2009
AÑO 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009
Precipitación
441,2 572,2 342,6 581,5 826,2 616 513,8 525 712,6 607,8 472,6
mm/año
199
Capítulo 6 Escenario futuro
En la figura 85 y 87 muestran la evolución de los usos del suelo por celdas del
PAI de Marina d´Or Golf y del PGOU Torreblanca 2009, respectivamente. La
figura 86 corresponde a la leyenda de la figura 85 y la figura 88 corresponde a la
leyenda de la figura 87.
200
Capítulo 6 Escenario futuro
201
Capítulo 6 Escenario futuro
202
Capítulo 6 Escenario futuro
203
Capítulo 6 Escenario futuro
Figura 88. Leyenda correspondiente a la figura 87. Doña Blanca Golf, superficies
urbanizadas y jardines correspondientes al PGOU Torreblanca (2014 - 2035)
- Bombeos
Las necesidades hídricas de la Plana de Oropesa - Torreblanca al inicio de la
simulación del escenario de actuación (año 2010) se satisfacen totalmente con aguas
subterráneas extraídas del acuífero pliocuaternario, con un caudal anual de 14,3
hm3. Debido a que los proyectos de actuación no están distribuidos de manera
homogénea en la Plana, la evolución de las extracciones existentes al principio de la
simulación (2010) variará de la siguiente manera:
204
Capítulo 6 Escenario futuro
Tabla 41. Reducción de los bombeos afectados por el PGOU Torreblanca (2010 – 2035)
205
Capítulo 6 Escenario futuro
206
Capítulo 6 Escenario futuro
Tabla 42. Reducción de los bombeos afectados por el PAI Marina d´Or Golf (2010 – 2035)
207
Capítulo 6 Escenario futuro
208
Capítulo 6 Escenario futuro
6.4.2. Resultados
La figura 91 muestra la evolución de la piezometría a lo largo de la simulación. Se han
seleccionado las imágenes correspondientes al periodo de verano por ser el periodo de
mayor stress hídrico del año, observándose tanto la piezometría como la dirección del flujo
subterráneo (representado por flechas); las zonas en las que se produce entrada de agua
marina (intrusión) han sido señaladas con un círculo rojo. La evolución piezométrica a lo
largo de la simulación (figura 91) muestra que las isopiezas negativas, que en un principio
se observan en el sector sur, han ido cambiando a niveles positivos, hasta observarse una
piezométría por encima del nivel del mar en todo el acuífero. Para conocer con mayor
exactitud cuándo ha remitido la intrusión (niveles por debajo del nivel del mar) se han
estudiado con mayor detalle las piezometrías de los años 2025, 2026 y 2027 (figura 92),
observándose que en otoño del 2026 no existen piezometrías negativas, por lo que la
intrusión comienza el proceso de retroceso.
Figura 91. Evolución de las piezometrías a lo largo del periodo de simulación (2010 -
2035)
209
Capítulo 6 Escenario futuro
Figura 92. Evolución de las piezometrías: verano 2035, verano 2026, otoño 2026 y
verano 2027
210
Capítulo 6 Escenario futuro
Tabla 43. Balances hídricos de los años 2010 y 2035 según el modelo matemático de flujo
(Visual Modflow)
Precipitación (mm)
Analizando las figuras 91 y 92, se puede observar que la intrusión marina es un proceso
que afecta únicamente al sector sur. Con objeto de estudiar en profundidad la zona
afectada por la intrusión marina, se ha calculado el balance hídrico únicamente en esa zona
(figura 93). En la tabla 45 (balance hídrico del sector sur) se observa numéricamente la
desaparición de la intrusión marina en el escenario futuro planteado, la disminución de la
recarga y de los bombeos, debido a la impermeabilización del terreno; además, cabe
recordar que el agua que se necesita para regar los campos de golf, y las zonas ajardinadas,
etc, procede de fuentes externas al acuífero de forma que todo ello propicia el aumento de
las salidas subterráneas al mar.
211
Capítulo 6 Escenario futuro
6.5. Conclusiones
La simulación llevada a cabo entre el año 2010 y el año 2035 se caracteriza por un
drástico cambio del uso del suelo. En este periodo se construyen cuatro campos de golf (en
los términos municipales de Oropesa, Cabanes y Torreblanca), se impermeabilizan sectores
asociados al PAI Marina D´Or Golf y al PGOU Torreblanca (2009), se disminuyen los
bombeos y en Marina d´Or Golf, el agua utilizada para suplir las necesidades hídricas es de
procedencia externa al acuífero pliocuaternario. Los resultados obtenidos en la simulación
muestran como la intrusión marina está en fase de retroceso, debido a la ausencia de
piezometrías por debajo del nivel del mar, lo que indica que las acciones llevadas a cabo en
la Plana de Oropesa - Torreblanca favorecen a la recuperación del acuífero.
212
Capítulo 7 Conclusiones y propuestas
Capítulo 7
Conclusiones y propuestas
7.1. Conclusiones
7.1.1. Situación histórica y vulnerable del sistema (índice SITE)
7.1.2. Modelo matemático de flujo subterráneo
7.1.3. Escenario futuro
7.1.4. Volúmenes ecológicos de remediación (VER)
7.2. Propuestas
7.2.1. Sobre modelización
7.2.2. Sobre escenario futuro
7.2.3. Sobre la evaluación de la vulnerabilidad del sistema (índice SITE)
7.2.4. Sobre el cálculo de los volúmenes ecológicos de remediación
213
Capítulo 7 Conclusiones y propuestas
214
Capítulo 7 Conclusiones y propuestas
7. CONCLUSIONES Y PROPUESTAS
7.1. Conclusiones
Se incluyen las principales conclusiones relativas al análisis de impactos de la intrusión
y posibles estrategias para mitigarlo en la Plana de Oropesa - Torreblanca. Se analizan los
posibles impactos bajo distintas perspectivas: el estudio de la información histórica, la
vulnerabilidad del sistema aplicando la metodología SITE, la modelación matemática del
sistema y el análisis del escenario futuro. Se propone también una metología para definir
estrategias de gestión apropiadas para mitigar esos impactos. Dicha metodología se basa en
la cuantificación del volumen ecológico de remediación (VER), definido como el volumen en
que hay que reducir los bombeos en un acuífero o sector del mismo para que la
recuperación piezométrica conseguida comporte una intrusión menor o igual a la estimada
en régimen natural para el sistema (intrusión de referencia).
215
Capítulo 7 Conclusiones y propuestas
• En la simulación (1973 – 2009) se observa que entre 1984 y 1994 tiene lugar la
aparición y el avance de la intrusión marina, con dos periodos en los que ha
sido más severa: 1985 – 1986 y 1992 – 1994. En estas fechas tuvieron lugar las
dos sequías más importantes acaecidas entre 1973 - 2009. En estos periodos de
máxima intensidad de la intrusión las salidas al mar eran aproximadamente de
5,5 hm3/año y las entadas de agua marina cercanas a los 2,3 hm3/año. Se
diferencian tres zonas especialmente afectadas por el proceso: la zona sur,
situada entre el Prat de Cabanes y Oropesa del Mar en el sector Sur, con un
volumen de entradas de agua de mar de 2 hm3/año; la zona de Torreblanca,
situada en la mitad septentrional (sector Central), entre la población de
Torreblanca y la costa, con una intrusión de 0,2 hm3/año; y la zona de
Alcossebre (sector Norte), con un caudal medio de entrada de agua de mar de
0,1 hm3/año.
• A partir de 1994 el proceso ha ido remitiendo hasta que en 2009 sólo afecta al
sector más meridional del acuífero, en el que la intrusión marina ha sido
cuantificada en 0,05 hm3/año. En este sector existe un claro déficit de recarga
lateral.
216
Capítulo 7 Conclusiones y propuestas
• Los modelos que simulan el régimen natural (sin bombeos) para los años seco
(modelo RNs), medio (modelo RNm) y húmedo (modelo RNh) sirven para evaluar
la intrusión marina de referencia existente en la Plana de Oropesa – Torreblanca
cuando la influencia antrópica era mínima. Para los 3 escenarios considerados
dicha intrusión natural es despreciable.
• Los modelos que simulan el régimen influenciado con la actual operación del
sistema (con los bombeos y retornos de riego actuales) para los años secos (RIs),
medios (RIm) y húmedos (RIh), muestran que en todos los casos existen niveles
piezométricos por debajo del nivel del mar en el sector Sur del acuífero y que
sólo en los años secos se hace patente en la zona de Torreblanca.
• En el caso de la zona Sur, el VER para el año seco es 2,8 hm3/año, para el año
medio 1,6 hm3/año y para el año húmedo 1,0 hm3/año.
7.2. Propuestas
7.2.1. Sobre modelización
• Se emplearon modelos de flujo de densidad constante en diferencias finitas
para evaluar impactos en un acuífero costero. Se asume que la existencia de
217
Capítulo 7 Conclusiones y propuestas
218
Capítulo 7 Conclusiones y propuestas
219
Capítulo 7 Conclusiones y propuestas
220
Bibliografía
Bibliografía
221
Bibliografía
222
Bibliografía
Albuixech, J. (1986). Flora del Desierto de las Palmas y sus alrededores. Col.lecció Universitària.
Diputació de Castelló.
Anderson M.P.; Woessner W. W. (1992). Applied Groundwater Modeling. Simulation of Flow
and Advective Transport. Academic Press.
Anónimo (1987). Gestión de los recursos hidráulicos de Cabanes y su Ribera. Documento inédito.
Antón - Pacheco, C.; Ballesteros, B.J.; Mejías, M.; Gumiel, J.C.; Gómez, J.A.; Callejo, P.
(2007). Identificación y cartografía de surgencias costeras en el acuífero kárstico de el
Maestrazgo(Castellón) mediante técnicas de teledetección aeroportada en el infrarrojo
térmico. Boletín Geológico y Minero, 118: 649-670.
Aragonés, J. M.; Pulido, J.; Sanz, C. (1976). Estudio hidrogeológico de la Plana de Oropesa –
Torreblanca (Castellón). I Simposio Nacional de Hidrogeología. I: 470-490. Valencia.
Ayuntamiento de Torreblanca (2009). Inventario de captaciones del término municipal de
Torreblanca. Documento inédito.
Ballesteros B.J. (2008). Tipología y caracterización de la intrusión marina en acuíferos litorales.
Método SITE. (primeras bases). Informe inédito. Instituto Geológico y Minero de
España.
Ballesteros B.J; Morell, I.; García, O.; Renau - Pruñonosa, A. (2012a): Propuesta de un nuevo
índice para la caracterización de la intrusión marina. El método SITE. IV Congreso
Internacional sobre la tecnología de la intrusión de agua de mar en acuíferos
costeros. TIAC´12. Alicante.
Ballesteros B.J; Morell, I.; García, O.; Renau - Pruñonosa, A. ( 2012b ). Proposal of a new
seawater intrusion index: the SITE method. Submitted.
Cabezas, F.; Cabrera, E.; Morell, I. (2008). El agua: una cuestión de estado. Perspectiva desde la
Comunidad Valenciana. Ed. AVE. 156 pp.
Canerot, J. (1974). Recherches geologiques aux confins des Chaines Iberiques et Catalanes. Espagne.
Université Paul Sabatier. Toulouse.
Canerot, J.; Faure, Ph.; Rahal, M. (1985a). Massif du Desierto de las Palmas (Castellón). Strata, 2
(2): 85-100.
Canerot, J.; Goy A. (1985b). Le Jurassique des Iberides Orientales (Espagne). Actes du laboratoire
de Géologie Sédimentaire et Paléontologie de l´Univertité Paul Sabatier. Toulouse. Sèrie
2: Mémories. Generalites Livret d´Excursion. Vol.2.
Carmona, P.; Ruiz, J. M. (2009). Cambios geomorfológicos ambientales en el litoral de Torre la
Sal durante el Holoceno. In: Torre la Sal (Ribera de Cabanes, Castellón). Evolución del
paisaje antrópico desde la prehistoria hasta el medioevo. Coordinador E. Flors.
Monografies de Prehistòria y Arqueologia Castellonenques, nº 8.
Confederación hidrográfica del Júcar (CHJ) www.chj.es
Confederación Hidrográfica de Júcar (CHJ, 2004). Seguimiento del Plan Hidrológico de la
Cuenca del Júcar. Documento de síntesis. Confederación Hidrográfica del Júcar.
Confederación Hidrográfica del Júcar (CHJ, 2005). Informe para la Comisión Europea sobre los
artículos 5 y 6 de la Directiva Marco del Agua. Confederación Hidrográfica del
Júcar. 528 pp.
Confederación Hidrográfica del Júcar (CHJ, 2007). Balances hídricos y evaluación del riesgo
conforme al estado cuantitativo de las masas de agua subterránea de la
Confederación Hidrográfica del Júcar. Informe interno. 66 pp.
223
Bibliografía
224
Bibliografía
García, O.; I. Morell, I.; Ballesteros, B. J.; Renau - Pruñonosa, A. (2012). Determinación de la
Intrusión marina en tres acuíferos costeros mediterráneos mediante el método SITE. IV
Congreso Internacional sobre la tecnología de la intrusión de agua de mar en
acuíferos costeros. TIAC´12. Alicante. Instituto Geológico y Minero de España. Vol 1 :
283 – 295.
García - Orellana, J.; García - Solsona, E.; Ballesteros, B.; López, J.; Mejías, M.; Marina , M.
(2006). Evaluación de la descarga de un manantial costero mediante isótopos de radio: las
Fuentes de Alcossebre (Castellón). Congreso internacional sobre el Agua Subterránea en
los Países Mediterráneos (AQUA in Med-06). Málaga.
Generalitat Valenciana (1990). Apéndice a la guía de la naturaleza de la Comunidad Valenciana:
Parques y paisajes naturales. 43. Paraje natural del Desert de les Palmes. Levante - El
Mercantil Valenciano.
Giménez, E; Morell, I. (1991). Consideraciones sobre la utilización de iones minoritarios en la
caracterización de la intrusión marina. El agua en Andalucía, I: 401 – 412. Córdoba.
Giménez, E. (1994). Caracterización hidrogeoquímica de los procesos de salinización en el acuífero
detrítico costero de la Plana de Castellón (España). Tesis Doctoral. Universidad de
Granada. Granada.
Giménez, E.; Fidelibus, M.D.; Morell, I. (1995). Metodología de análisis de facies hidroquímicas
aplicada al estudio de la intrusión marina en acuíferos detríticos costeros. Aplicación a la
Plana de Oropesa (Castellón). Hidrogeología 11.
Giménez, E.; Morell, I. (1997). Hidrogeochemical analysis of salinization processes in the coastal
aquifer of Oropesa (Castellón, Spain). Enviromental Geology.
Gómez, J. D., Ballesteros, B. J., Morell, I., Renau - Pruñonosa, A. (2010). Characterization of
seawater intrusión in coastal aquifers by application of index SITE. European
Groundwater Conference. Madrid, España.
Gómez, J. D., Díaz, J. A., Meléndez, M. (2012). Cálculo de descargas ambientales al mar en
acuíferos detríticos costeros mediante modelos matemáticos de flujo. Aplicaciones en la
demarcación del Júcar. IV Congreso Internacional sobre la tecnología de la intrusión de
agua de mar en acuíferos costeros. TIAC´12. Alicante. Instituto Geológico y Minero
de España. Vol 1 : 127 - 237.
Gómez - Serrano, M. A.; Mayoral, O.; Domingo, J. (2001). Guía de la Naturaleza del litoral de
Castellón. Itinerarios para conocer su fauna, flora, paisaje e historia. Ed. Antinea. Vinaroz.
Gómez - Serrano, M. A.; Calabuig, J.D.; Mayoral, O. (1999). Vegetación litoral y cambios en el
paisaje de la provincia de Castellón. Ayuntamiento de Castellón de la Plana. Castellón
de la Plana.
González del Tánago, M.; García de Jalón, D. ( 2006a ). Attributes for assessing the
environmental quality of riparian zones. Limnetica, 25, 389-402.
González del Tánago, M.; García de Jalón, D. (2006b). Caracterización jerárquica de los ríos
españoles. Propuesta de tipología de tramos fluviales para su clasificacion atendiendo a la
directiva marco del agua. Limnetica. 25, 81-98.
González - Hidalgo J. C. (2008). Proyecto Cambio Climático: base de datos de recipitaciones,
análisis de tendencias e impactos en los sistemas naturales. Ministerio de Ciencia e
Innovación CGL2008-05112-C02-01/CLI.
González - Hidalgo J. C., Brunetti M., de Luis M. (2010) A new tool for monthly precipitation
analysis in Spain: MOPREDAS database (Monthly precipitation trends December 1945-
November 2005). International Journal of Climatology, DOI: 10.1002/joc. 2115.
225
Bibliografía
Herrera - Carvajal, S.; Villalba - Atondo, I.; Molinar - Tabares, E. ( 2007 ). Modelación
hidrogeológica y evaluación de la sustentabilidad del acuífero transfronterizo del Río San
Pedro en Sonora, México. Revista latinoamericana de Recursos Naturales, 3(1):21-30.
Ibáñez S. P. (2005). Comparación de la aplicación de distintos modelos matemáticos sobre acuíferos
costeros detríticos. Tesis Doctoral. Universidad de Granada. Granada.
Institut Cartogràfic Valencià (ICV, 2008). Hojas topográficas 1: 10000. Generalitat
Valenciana. Conselleria de Medi Ambient, Aigua, Urbanisme i Habitatge. Edició
2006-2008. Valencia.
Instituto Geológico y Minero de España (IGME, 1975). Hidrogeología de la Plana Oropesa –
Torreblanca. Proyecto de Investigación Hidrogeológica de la Cuenca Media y Baja del Júcar.
Ministerio de Industria. Madrid.
Instituto Geológico y Minero de España (IGME, 1977). Evolución de niveles piezométricos en
los Sistemas Acuíferos Cuenca Media y baja del Júcar. Servicio de publicaciones.
Ministerio de Industria. Madrid.
Instituto Geológico y Minero de España (IGME, 2000). Estado de la intrusión de agua de mar
en los acuíferos costeros españoles. Vol.II: Cuencas mediterráneas I: Segura, Júcar y Baleares.
Instituto Geológico y Minero de España. Madrid.
Instituto Geológico y Minero de España; Universitat Jaume I de Castelló (IGME – UJI,
2009a). Propuesta de indicador de estado de intrusión marina. Aplicación a las masas de agua
subterránea de la provincia de Castellón. Instituto Geológico y Minero de España y
Dirección General del Agua. Madrid.
Instituto Geológico y Minero de España; Universitat Jaume I de Castelló (IGME – UJI,
2009b). Descargas ambientales al mar en las masas de agua subterránea costeras de la
provincia de Castellón (Cuenca del Júcar). Instituto Geológico y Minero de España y
Dirección General del Agua. Madrid.
Instituto Nacional de Estadística (INE, 2011). www.ine.es
Instituto Nacional de Reforma y Desarrollo Agrario (IRYDA, 1983). Estudio geomorfológico de
la marjal de Cabanes - Torreblanca. Ministerio de Agricultura, Pesca y Alimentación.
Inédito. Madrid.
Instituto Tecnológico Geominero de España (ITGE, 1989). Estudio hidrogeológico del
Maestrazgo. Castellón. Informe interno.
Instituto Valenciano de Investigaciones Agrarias (IVIA). www.ivia.es .
Iríbar V.; Carrera J.; Custodio, E.; Medina, A. (1997). Inverse modeling of seawater intrusion in
the Llobregat delta deep aquifer. Journal of Hydrology. Vol 44: 198-226.
Júnior, R. F.; Fernández - Rubio, R.; Baquero - Úbeda, J.C.; Lorca - Fernández D.; Verdejo –
Serrano, J. (2003). Acuífero detrítico costero de Ponta da Madeira (Maranhao, Brasil).
Modelización hidrogeológica preliminar. Tecnología de la intrusión de mar en acuíferos
costeros: Países mediterráneos. IGME. Madrid. ISBN: 84-7840-470-8.
Kantorovich L. V.; Krylov V. I. (1964). Approximate Methods of Higher Analysis. Interscience
Publishers.
Marina, M. (2006). Caracterización hidrogeológica de formaciones acuíferas profundas. Aplicación
al acuífero regional jurásico de la unidad hidrogeológica 08.07 (el Maestrazgo). Boletín
Geológico y Minero. Vol. 117: 605-6011.
Martín, J. D. (1998). Geologia i mineralogia de les lutites permotriàsiques de la Serra de Vilafamés
(Castelló). Tesis de Licenciatura. Universitat Autònoma de Barcelona. Barcelona.
226
Bibliografía
Mateu, J.F. (1982). El norte del País Valenciano. Geomorfología litoral y prelitoral. Tesis Doctoral.
Universidad de Valencia. Valencia.
McDonald M.G.; Harbaugh A. W. (1984). A modular three-dimensional finite-difference ground-
water flow model: U.S. Geological Survey. Open-File. Report 83-875, 528 p.
McDonald M. G.; Harbaugh A. W. (1988). Chapter A1. A Modular Three-Dimensional Finite-
Difference Ground-Water Flow Model. Book 6. Modeling Techniques. USGS.
(http://pubs.usgs.gov/twri/twri6a1/pdf/).
Mejías, M.; García - Orellana, J.; Plata, J.L.; Marina, M.; García - Solsona E.; Ballesteros, B.;
Masqué, P.; López. J.; Fernández - Arrojo, C. (2007). Methodology of hydrogeological
characterization of deep carbonate aquifers as potential reservoirs of groundwater: Case of
study: the jurassic aquifer of El Maestrazgo (Castellón, Spain). Environmental Geology
DOI. 10.1007/S00254-007-0845-0.
Meléndez, M.; Díaz, J. A.; Gómez, J.D. (2012). Estimación de las descargas ambientales al mar y
el estado de la intrusión marina en acuíferos costeros de la zona central de la cornisa
cantábrica de España (Asturias y Cantabria). IV Congreso Internacional sobre la
tecnología de la intrusión de agua de mar en acuíferos costeros. TIAC´12. Alicante.
Instituto Geológico y Minero de España. Vol 1 : 489 - 499.
Ministerio de Medio Ambiente (MIMAM, 1996). Guía oficial de las playas de España: las Playas
de Levante, de Castellón a Murcia. Madrid.
Ministerio de Medio Ambiente (MIMAM, 2000). Libro Blanco del Agua en España. Madrid.
Ministerio de Medio Ambiente (MIMAM, 2005). Estudio inicial para la identificación y
caracterización de las masas de agua subterránea de las cuencas intercomunitarias. Dirección
General del Agua. Ministerio de Medio Ambiente.
Molina, J.A.; Casemiro, M.A.; Pertiñez, C.; Moreno, P.S. (2000). Relación suelo-vegetación en
un ecosistema litoral mediterráneo. El parque natural del Prat de Cabanes - Torreblanca
(Castellón). Fundación Dávalos-Fletcher.
Morell, I. (1985). Caracterización hidroquímica de la intrusión marina en la Plana de Oropesa-
Torreblanca (prov. Castellón). Tesis Doctoral. Universidad de Granada. Granada.
Morell, I.; Pulido - Bosch, A. (1986). Study of the freshwater - saltwater interface in the coastal
aquifer of Oropesa - Torreblanca, Castellón, Spain. Proc. 9th SWIM, pp: 61 - 72. Delft.
Morell, I.; Giménez, E.; Esteller, M.V. (1987). Ejemplo de utilización de registros verticales de
conductividad, temperatura, pH y oxigeno disuelto. Aplicación a los acuíferos de Oropesa –
Torreblanca y Moncofar (Castellón). Hidrología y Recursos Hídráulicos, XI: 171-182.
Morell, I. (1988). Intrusión marina y salinización en la Comunidad Valenciana. Jornadas del
Agua. Factor: el Desarrollo Valenciano. Valencia.
Morell, I.; Giménez, E.; Esteller, M.V. (1988). Comportamiento iónico y procesos físicoquímicos
en acuíferos detríticos costeros de las Planas de Oropesa, Castellón y Gandía. Hidrogeología
3: 21-34.
Morell I. (1989). Salinización por intrusión marina. Aplicación a los acuíferos costeros de la
provincia de Castellón. Col.lecció Universitària. . Ed. Excma. Diputación de Castellón.
ISBN 84-8695-07-3.
Morell, I. (1992). Manantiales de la provincia de Castellón. Col.leccio Universitaria. Ed. Excma.
Diputación de Castellón.
227
Bibliografía
Morell, I.; Esteller, M. V.; Cerezo, M.; Lapeña, L.; García, P. (1996a). Balance de masas de
nitratos en una parcela experimental de cítricos regados con agua subterránea y agua residual
urbana depurada. Contribuciones a la investigación y gestión del agua subterránea.
Morell & Fagundo (Eds).
Morell, I.; Tuñón, J.; Gascón, M.; Jiménez, J.R. (1996b). Caracterización hidroquímica de flujos
preferenciales en un acuífero detrítico. Contribuciones a la investigación y gestión del
agua subterránea. Morell & Fagundo (Eds). Pp:179 - 189.
Morell, I.; Padilla, F.; Alberto, J. M. (1996c). Modelling of lecahing of chloride and nitrogen
species in an experimental citrus grove. Ecological Modelling. International journal on
ecological modelling and systems ecology (87). Elsevier. Pp:235-247.
Morell, I; Bandenay, L.; Renau – Llorens, E. A.; Renau – Pruñonosa, A. (2006). Recarga
artificial de acuíferos con aguas residuals urbanas regeneradas. Ayudas a la investigación
2006. Fundación dávalos – Fetcher. ISBN: 978 – 84 – 693 – 2077 – 8.
Morell I. (2008). Actualización del inventario de puntos de agua de la Plan a de Oropesa -
Torreblanca. Informe inédito.
Morell, I. (2010). Funcionamiento hidrogeológico de la zona húmeda denominada Prat de Cabanes
(Torreblanca). Origen del agua y de la salinidad. Procesos implicados. Informe técnico.
Instituto Universitario de Plaguicidas y Aguas (IUPA). Castellón.
Morell, I.; Ballesteros, B. J.; Renau - Pruñonosa; A, Gómez, J. D. (2010). Application of a new
seawater intrusion index: SITE method. Keynote in Global Forum on Salinization and
Climate Change. FAO-Valencia University. Valencia.
Peaceman, D.W. (1977). Fundamentals of numerical reservoir simulation. Elsevier.
Pérez M. A. (2005). Modelo distribuido de simulación del ciclo hidrológico y calidad del agua,
integrado en sistemas de información geográfica, para grandes cuencas. Aportación al análisis
de presiones e impactos de la Directiva Marco del Agua. Tesis Doctoral. Universidad
Politécnica de Valencia. Valencia.
Pisinaras, V.; Petalas C.; Tsihrintzis; V. A.; Zegana, E. (2007). A groundwater flow model for
water resources Management in the Ismarida plain, North Greece. Environmental
Modeling an Assessment . Springer. Pp: 75-89.
Plan General de Ordenación Urbana de Torreblanca ( PGOU Torreblanca, 2009 ).
Clasificación del suelo e instalaciones náuticas 1:10000. www.torreblanca.es .
Plan Hidrológico de la Cuenca del Júcar (CHJ, 1997). Plan Hidrológico de la Cuenca del Júcar.
Confederación Hidrográfica del Júcar. Valencia.
Prickett, T. A.; C. G. Lonnquist (1971). Selected Digital Computer Techniques for Groundwater
Resources Evaluation. State Water Survey, Urbana, III.,Bull. 55.
Psarropoulou, E. T; KAratzas, G. P. (2012). Transient groundwater modelling with spatio -
temporally variable fluxes in complex aquifer systems: new approach in defining boundary
conditions for transient flow model. Civil Engineering and environmental systems.
Vol29(1):1 – 21.
Pulido A.; Molina, L. (2006). El agua y el medio ambiente. Instituto de Estudios Almerienses.
Colección Ciencia y Tecnología. Nº 2. pp. 159.
Reeve, A. S.; Warzocha, J.; Glaser, P. H.; Siegel, D. I. (2000). Regional ground-water flow
modeling of the Glaciar Lake Agassiz Peatlands, Minnesota. Journal of Hydrology.
Elsevier.
228
Bibliografía
Rejani, R.; Jha, M. K.; Panda, S. N.; Mull, R. (2008). Simulation modelling for efficient
groundwater management in Balasore coastal basin, India. Water Resources Management.
Vol: 22 (1): 23 - 50.
Remson, I.; Hornberger, G.M.; Molz, F.J. ( 1971). Numerical Methods in Subsurface Hydrology.
New York: Wiley-Interscience. Pp: 389.
Renau - Llorens, E. A. (2010). Elementos minoritarios y traza en la masa de agua subterránea
080.021 (Plana de Castelló). Origen y procesos asociados. Tesis Doctoral. Universitat
Jaume I. Castellón de la Plana.
Renau - Pruñonosa, A. (2008). Estudio hidrogeológico del sector meridional de la Plana de
Oropesa - Torreblanca (Castellón). Trabajo fin de posgrado. Curso Internacional de
Hidrogeología Subterránea. FCIHS y Universitat Politècnica de Catalunya.
Barcelona.
Renau - Pruñonosa, A. (2009). Cálculo de volúmenes ambientales de descarga al mar en situación
de explotación sostenible. Aplicación a la Plana de Oropesa - Torreblanca (prov. Castellón).
Trabajo de Investigación (DEA). Universitat Jaume I. Castellón de la Plana.
Sanchis, E.; Rodríguez, T.; Morell, I. (1989). La naturaleza de la Comunidad Valenciana. El
medio físico, Geología. Institución Valenciana de Estudios e Investigación. Edicions Alfons
el Magnànim.
Serrano, V.; Toca, A.; Pérez, A.; Nieto, M.; Ortiz, S.; Capilla, J. (1995). Estudio de investigación
de las descargas de aguas subterráneas al mar por la Sierra de Irta (U.H. 08.07: Maestrazgo).
Mancomunidad del Baix Maestrat, Castellón. IV Simposio Hidrogeología. Tomo XIX: 679
- 697.
Servicio Geológico de Obras Públicas (SGOP, 1975). Estudio hidrogeológico de la Plana de
Oropesa-Torreblanca para mejorar el abastecimiento a Oropesa del Mar. Ministerio de
Obras públicas y Urbanismo. Madrid.
Servicio Geológico de Obras Pública (SGOP, 1977). Estudio hidrogeológico de la Plana de
Castellón. Ministerio de Obras públicas y Urbanismo. Madrid.
Southwell R.V. (1940). Relaxation Methods in Engineering Science. Oxford University Press.
Stone, H. L. (1968). Interactive Solution of Implicit Approximations of Multidimensional Partial
Differential Equations. SIAM Journal of Numerical Analysis 5: 530-558.
Témez, J.R. (1977). Modelo matemático de transformación precipitación aportación. ASINEL.
Tuñón, J. (2000). Determinación experimental del balance hídrico del suelo y evaluación de la
contaminación asociada a las prácticas agrícolas. Tesis Doctoral. Universitat Jaume I,
Castellón de la Plana.
Unión Europea (UE - DMA, 2000). Directiva Marco del Agua, “Directiva 2000 / 60 / CE del
Parlamento Europeo y del Consejo de 23 de Octubre de 2000”. Diario Oficial de las
Comunidades Europeas de 22/12/2000. L 327/1-327/32.
Unión Europea (UE, 2006). Directiva 2006/118/CE del Parlamento Europeo y del Consejo de 12
de diciembre de 2006. Diario Oficial de la Unión Europea de 27/12/2006. L 372/19 –
372/31.
Viñals, M. J. (2004). Repercussions de l´acció antròpica de les albuferes, marjals i altres zones
humides costaneres valencianes. Anuari de l´Agrupació Borrianenca de Cultura. 15:109-
129.
229
Bibliografía
230
Anexos
Anexo 1
A1-1. Propuesta de un nuevo índice para la
caracterización de la intrusión marina.
El método SITE
A3
Anexo 1 Índice SITE
A4
PROPUESTA DE UN NUEVO ÍNDICE PARA LA CARACTERIZACIÓN
DE LA INTRUSIÓN MARINA. EL MÉTODO SITE
B.J. Ballesteros(1), I. Morell(2), O. García(1) y A. Renau-Pruñonosa(2)
(1) Unidad de Valencia. Instituto Geológico y Minero de España. c/ Cirilo Amorós, 42, entreplanta, 46009
Valencia, [email protected]. (2) Instituto Universitario de Plaguicidas y Aguas. Universitat Jaume I. Castellón,
España
RESUMEN
Conseguir el buen estado cuantitativo y cualitativo de las masas de agua subterránea es el objetivo con el que
se diseñaron las Directivas 2000/60/EC y 2006/118/EC, uno de cuyos fines es alcanzar dicha situación hacia
el año 2015. Para ello es necesario disponer de metodologías que permitan caracterizar el estado de los
distintos sistemas hidrogeológicos de una forma lo más objetiva posible. En el caso de los acuíferos costeros
los procesos de intrusión marina generan una elevada complejidad hidroquímica en sus aguas, por lo que sólo
metodologías basadas en aproximaciones simples y en criterios elementales pueden ser aplicadas de forma
universal y con resultados prácticos. Para ello, estos métodos deben cumplir dos condiciones: 1) facilidad de
elaboración a partir de datos sencillos y habitualmente disponibles 2) ser capaces de proporcionar
información realista, discriminatoria y comprensible. De acuerdo con esto el método SITE permite
caracterizar la intrusión marina sufrida por un acuífero costero con la determinación de cuatro parámetros
basados exclusivamente en el contenido en ion cloruro de sus aguas. El acrónimo SITE refleja las iniciales de
dichos parámetros: Superficie afectada, Intensidad, Temporalidad y Evolución a medio-largo plazo de la
salinización. Los resultados se concretan en la obtención de un índice numérico y de valores tipológicos,
fáciles de conceptualizar, que proporcionan información acerca de las características más representativas del
proceso intrusivo. Dicha información permite conocer de forma estandarizada la situación de las masas de
agua subterránea, y poder ser comparadas entre sí, de manera que se pueden establecer prioridades de
actuación, optimizando los recursos y orientando adecuadamente la mejora de su gestión mediante la
identificación objetiva de los problemas.
Palabras clave: Intrusión marina, acuíferos costeros, índice, gestión de aguas subterráneas
ABSTRACT
One of the main objectives of the Directives 2000/60/EC and 2006/118/EC is to achieve good quantitative
and chemical status of groundwater bodies by 2015, accordingly there is need to have available
methodologies that allow to define the condition of the different hydrogeological systems in a most objective
way. The seawater intrusion processes in coastal aquifers leads to the occurrence of a high complexity of a
groundwater quality. Consequently methodologies only based in very simple approaches and elemental
criteria can be applied in a universal pattern with practical results, therefore these methods must meet two
basic requirements: 1) easiness of elaboration from simple and readily avaible data and 2) ability to provide
realistic, discriminatory and understandable information. In accordance with that, the method SITE allows to
characterize the seawater intrusion process in a coast aquifer by the determination of four parameters only
based on chloride concentration. The SITE acronym reflects the initials of these parameters: S=Spatial
affection, I=Intensity, T=Temporal variability or seasonality and E=Evolution to medium-long term of
salinization. The key output of the method is an numerical index and typological values, easy to
conceptualize, that provide information about the most representative factors of the intrusive process. This
information will be very useful to know the real conditions of groundwater bodies in a standardized way and
to be compared with each others in order to prioritize the interventions, improving their resources and leading
suitably their management by the objective identification of their problems.
1. INTRODUCCIÓN
Figura 1. Consideración hidrostática de la intrusión marina según Ghijben-Herzberg (izda.). Diferentes zonas del equilibrio
agua dulce–agua salada por la anisotropía del acuífero y las variaciones verticales de la conductividad hidráulica (centro).
Efecto upconing (dcha.)
Las propuestas de índices de intrusión marina, con muy diferente complejidad y éxito, han sido
numerosas. En una primera aproximación, pueden ser clasificados en las siguientes categorías:
basados en la relación explotación/recarga, de calidad general, basados en el contenido en ión cloruro,
hidrogeoquímicos de elementos mayoritarios y minoritarios, mixtos a partir de datos morfológicos,
piezométricos y de calidad, basados en métodos discretizados para la determinación de la
vulnerabilidad a la contaminación y, por último, tipológicos o descriptivos. Respecto al primero de
ellos, el índice de explotación de un acuífero puede utilizarse para informar sobre la presión global
soportada por un sistema hidrogeológico. Su aparente simplicidad, sin embargo, no debe ocultar que
se basa en dos conceptos cuya determinación está sujeta a notables imprecisiones, y quela explotación
de un acuífero no es homogénea. Los índices referidos a la calidad general del agua han suelen basarse
en el contenido total en sólidos disueltos (Voudouris et al., 2004), el residuo seco y las medidas de
conductividad eléctrica. Sin embargo, los de mayor interés son los que utilizan el ión cloruro, que es
un excelente guía para conocer el porcentaje de agua salada en la mezcla debido a su carácter de soluto
no reactivo. Otros índices se basan en métodos de distribución espacial de la salinidad en el acuífero,
como el de Murillo y Navarro (1991) que proponen un Índice Particularizado de Intrusión marina
(IPI), o el porcentaje de superficie de acuífero con contenido en cloruros superior a 1.000 mg/l
(MARM, 2009).
Las relaciones iónicas equivalentes de elementos mayoritarios son también utilizadas para
evaluar la intrusión marina. Las más comúnmente utilizadas son Mg/Ca, Cl/HCO3, SO4/Cl y Na/Cl. En
cuanto a las relaciones iónicas entre iones minoritarios, pueden proporcionar información de interés
las Br/Cl, I/Cl y Sr/Cl (Morell, 1985; Jiménez, 1994, Giménez y Morell, 1997; Custodio y Alcalá,
2003; Mandilaras et al., 2007). De todos ellos, el más utilizado es el Br/Cl (o Cl/Br) que tiene su mejor
aplicación en la identificación del origen de la salinidad. Intentos de aplicar nuevos índices
hidrogeoquímicos con determinaciones analíticas complejas son los basados la relación isotópica
B10/B11 que permite discriminar orígenes de salinidad y evaluar procesos de mezcla (Giménez y
Morell, 2008; Morell y Pulido-Bosch, 2008).
Por otro lado, índices catalogados como mixtos son los de Navarro y Murillo (1997 y 2007) el
cual, a partir de las características geométricas del acuífero, nivel piezométrico y contenido en ión
cloruro, permite estimar los volúmenes de agua continental y marina. En esta misma línea, Edet y
Okereke (2001) proponen un Indice de Valoración (AI) semicuantitativo de la intrusión marina,
basado en el control de cinco parámetros: total de sólidos disueltos, densidad, concentración de
cloruros y de sodio, y relación Br/Cl. También se han utilizado algunos métodos discretizados para la
determinación de la vulnerabilidad a la contaminación, tales como los modelos GALDIT (Chachadi,
2005) y GOD (Ertekin, 2005). Por último, los índices tipológicos o descriptivos se fundamentan en la
determinación de parámetros que reflejan las principales características del proceso intrusivo. A esta
clase de índices pertenece el del tipo SITE, cuyo desarrollo metodológico se expone en este artículo.
Las primeras bases del método SITE fueron establecidas por Ballesteros (2008), más tarde los
trabajos desarrollados para su aplicación a varios acuíferos mediterráneos (IGME-DGA, 2010) y,
especialmente, los realizados para la elaboración de este artículo, permitieron terminar de perfilar los
procedimientos de estimación de los parámetros que conforman el índice. Conceptualmente, el índice
se basa en la idea de que la intrusión marina puede ser establecida desde dos puntos de vista: el
espacial y el temporal, y puede procederse a su caracterización a partir de cuatro parámetros
elementales: superficie de acuífero afectada, intensidad del proceso, carácter temporal de la afección y
evolución de la calidad a medio-largo plazo. Debido al carácter tridimensional del fenómeno intrusivo,
estos son diferentes para cada profundidad del agua en el acuífero, así como para una determinada
fecha o momento. Como se explica más adelante, para ser realmente útil el método necesita de la
asunción de unas ciertas simplificaciones conceptuales. Con estas consideraciones se ha diseñado el
método SITE, cuyo acrónimo responde a los cuatros parámetros mencionados relacionados con el
proceso intrusivo:
Por otro lado, el índice utiliza información de fácil manejo y habitualmente disponible, ya que
se trata de que sea sencillo de elaborar y fácil de entender, así como de ofrecer información elemental
que permita, por un lado, evaluar la situación general de un acuífero y, por otro, discriminar entre
acuíferos distintos. En este aspecto, y como se expone más adelante, el método utiliza exclusivamente
como dato de partida el contenido en ion cloruro del agua subterránea a lo largo del tiempo en
diferentes puntos del sistema hidrogeológico. Con esta información, el método es capaz de diferenciar
entre distintas tipologías, representadas bajo símbolos y números, que permiten definir de forma
cuantitativa y/o cualitativa la intrusión marina sufrida por un acuífero.
Establecida en apartados anteriores la complejidad del proceso, para conseguir los objetivos de
practicidad y aplicabilidad perseguidos por el método, es necesario asumir un modelo conceptual
simple, y común a todos ellos, consistente en la existencia de dos fases (agua dulce y salada) que se
mezclan en el seno de un medio homogéneo. En otro orden de cosas, dado que el fenómeno intrusivo
es de carácter tridimensional, lo más adecuado sería estructurar los datos de la salinidad del agua
(referidos a contenido en ión cloruro) según diferentes rangos de profundidad, con objeto de definir
espacialmente la cuña salina. Sin embargo, el índice propuesto trata los datos de manera
bidimensional, es decir, obvia el parámetro profundidad. Las razones que han llevado a adoptar esta
importante y sustancial simplificación pueden resumirse, básicamente, en las siguientes: la gran
anisotropía hidrodinámica de los acuíferos que complica extraordinariamente la definición de la cuña
salina, la habitual falta de información sobre las características constructivas de las perforaciones y de
los horizontes transmisivos atravesados (columna litológica), y la existencia de diferentes cargas
hidráulicas en los distintos tramos permeables que generan flujos ascendentes y descendentes en las
perforaciones. Las razones expuestas implican que es difícil, y normalmente muy costoso, poder
asignar la salinidad del agua a diferentes profundidades del acuífero y, por tanto, la obtención de datos
fiables. Además, tener en cuenta todas estas circunstancias dificultaría enormemente la utilización del
índice y restringiría su aplicación a sólo algunos pocos casos muy estudiados y con información
pormenorizada. En consecuencia, el método carecería de la cualidad de poder ser empleado de forma
universal, que es uno de sus principales objetivos. También conllevaría un notable esfuerzo, sin que en
la práctica se consiguieran resultados sustancialmente diferentes, hecho que ha podido ser constatado
en las pruebas realizadas en los acuíferos en los que se ha desarrollado el método.
De acuerdo con lo expuesto en el apartado anterior, el índice propuesto utiliza como único
parámetro la concentración de ión cloruro expresado en mg/l, siendo lo ideal disponer de dos medidas
anuales que correspondan, respectivamente, a situaciones de mayor y menor salinidad (invierno-
verano, primavera-otoño o cualquier otra específica del acuífero). El número de puntos deberá ser
suficiente como para poder elaborar un mapa de isocloruros con un adecuado grado de detalle. A partir
de los datos básicos se obtendrán las matrices de datos específicas para los diferentes parámetros
SITE, que deberán confeccionarse con los criterios más adecuados para obtener el mayor rendimiento
de la información disponible. Se deben eliminarán los datos claramente anómalos (errores de lectura,
trascripción, etc.), los que carezcan de información suficiente, así como puntos muy próximos.
Ocasionalmente, se podrán completar datos con la media de los valores inmediatamente anterior y
posterior, asumiendo que las variaciones de salinidad son graduales.
Se define el valor de referencia como el contenido en ion cloruro por encima del cual se
considera que el agua de un sistema está afectada por la intrusión marina. En principio, dicho valor
debe venir dado por el fondo geoquímico natural de este compuesto en el acuífero, sin embargo, su
determinación no es trivial y puede requerir un esfuerzo investigador considerable. En muchos casos,
este valor está relacionado con la concentración del agua de alimentación subterránea procedente de
los acuíferos limítrofes, pero pueden existir otros mecanismos que alteren este fondo, tanto naturales
como inducidos. En consecuencia, para simplificar el método se proponen dos opciones: utilizar como
valor de referencia el fondo geoquímico en el caso de que pueda ser determinado con facilidad y
seguridad, o utilizar el valor de 250 mg/l de contenido en ion cloruro, límite establecido como
referencia para el agua de consumo humano en algunos países.
Uno de los elementos de entrada del índice es la superficie total del acuífero, expresada en km2.
Sin embargo, para los de gran extensión, o con una morfología alargada perpendicular a la costa, la
aplicación del índice a toda su superficie puede distorsionar los resultados. Esto se debe a que el
proceso intrusivo queda limitado a una determinada distancia a la línea de costa, ya que hay zonas en
el interior de los acuíferos por las que no puede extenderse al existir condicionantes, tales como la
elevación de la base impermeable o las propias dimensiones y geometría del sistema, que lo impiden o
dificultan. Al igual que en el caso del valor de referencia, esta distancia puede ser determinada si se
conoce muy bien el acuífero, pero el proceso entraña cierta complejidad. Para obtener resultados
comparables, y como solución de compromiso, se lleva a cabo la corrección de superficie consistente
en limitar la superficie potencialmente afectable superficie afectable por la intrusión marina a la
localizada a menos de 10 km de la línea de costa, más allá de la cual se estima que difícilmente pueden
darse procesos intrusivos. Esta sería la considerada, a todos los efectos de la confección del índice
como la extensión del acuífero, e implica que en el cálculo de sus distintos parámetros no se tienen en
cuenta los datos externos a ella.
5. PARÁMETROS SITE
Establecidas las condiciones y datos de partida, así como las definiciones conceptuales, se
describe la determinación de los cuatro parámetros en los que se basa el índice.
Tiene por objeto cuantificar la afección espacial causada en el acuífero por la intrusión marina,
independientemente de su intensidad. Se determina como la superficie (S) del sistema, en proyección
horizontal, afectada por el fenómeno intrusivo, y viene definida por aquélla en la que las aguas
subterráneas superan el valor de referencia. Este criterio permite diferenciar cinco tipos de acuíferos
según la mayor o menor extensión relativa del proceso: sin intrusión o puntual, localizada, zonal,
general y total. El parámetro S se calcula a partir de la isolínea del valor de referencia (Vr) del mapa de
isoconcentración del ion cloruro correspondiente a la situación actual. Como se ha comentado, para
casos sin una casuística específica dicho valor de referencia será 250 mg/l. También se determinarán
las isolíneas (en mg/l) correspondientes a 0,6Vr, 2Vr, 4Vr, 6Vr y 8Vr, que servirán para el cálculo del
próximo parámetro (a partir de la isolínea de 8Vr mg/l se podrán establecer intervalos mayores). En
consecuencia, el cálculo de S se determina según:
S >Vr
S=
S total
donde Stotal es la superficie total del acuífero, que puede estar corregida (superficie afectable), y S>Vr la
superficie con un contenido en cloruros superior al valor de referencia. Para cada rango de valores de
S se asigna un valor, o índice numérico, y un código de identificación (tabla 1).
Este parámetro persigue obtener una aproximación a la salinidad media del acuífero para la
situación actual o, en su defecto, para los años más recientes disponibles. Debe entenderse que ésta no
tiene en cuenta de manera explícita variaciones espaciales ni tendencias, aunque, en cualquier caso, su
valor depende de dichas variaciones y tendencias. Este criterio permite diferenciar cinco tipos de
acuíferos según la intensidad de la intrusión: nula, baja, media, alta y extrema. Para la obtención de
este parámetro I se parte del mapa de isocontenido en cloruros para la situación actual, ya elaborado
para el cálculo del parámetro S. A partir de él se calculan las diferentes superficies comprendidas entre
isolíneas (Si) consecutivas, y a cada superficie se le asigna el valor medio de dichas isolíneas (Cli). A
la superficie con contenido inferior a 0,6Vr mg/l se le otorga, por convención, un valor medio de
0,4Vr mg/l, mientras que a la superficie que supere la isolínea más elevada se le atribuye el valor
medio entre dicha isolínea y la que correspondería a la inmediatamente superior. Finalmente, se
calcula el parámetro I según:
I=
∑S i ⋅ Cli
(mg/l)
S total
donde Si es la superficie en km2 existente entre la isolínea i y su inmediatamente superior, y Cli el
contenido medio en cloruros en mg/l asignado a dicha superficie. Para cada rango de valores de I se
asigna un valor numérico y un código de identificación (tabla 1).
Para el cálculo de este parámetro, lo oportuno es disponer de dos o más datos anuales en cada
punto y para las fechas en las que se obtenga, respectivamente, la mejor y peor calidad del agua en el
sistema. En este caso, estrictamente, no es necesaria la continuidad temporal de los datos, ya que sólo
se requiere conocer las diferencias entre el valor máximo y mínimo de un número determinado de
años. Cuando no sea posible la elaboración de las dos series (un solo dato anual por punto) el método
también aporta un cálculo de T alternativo. En el caso de disponer de dos o más campañas anuales se
procede a hallar la media de las diferencias anuales entre las concentraciones máxima y mínima de
cloruros en los puntos de control. Como se puede deducir, si existen más de dos datos al año por punto
sólo se tendrá en cuenta el valor más alto y el más bajo. Este proceso es equivalente al que se
obtendría al confeccionar dos series de salinidad, una de valores máximos anuales y otra de mínimos,
en las que la diferencia de sus medias reflejaría el comportamiento del sistema. Si el valor hallado es
pequeño se entiende que el sistema es muy inercial y presenta una gran resilencia al proceso intrusivo,
pero también una gran dificultad para conseguir recuperar la calidad de sus aguas. Por el contrario, los
valores elevados implican que el sistema tiene muy poca inercia, lo que supone una mayor facilidad
para desarrollar procesos intrusivos, pero también una mejor aptitud para poder ser recuperado. En
consecuencia, en este caso el parámetro T viene dado por la siguiente expresión:
1 n
T= ∑ ( f ( x) max − f ( x) min)
n x =1
(mg/l)
donde f (x) max corresponde los valores máximos anuales en todo el acuífero y f (x) min a los
valores mínimos anuales, ambos en mg/l. En el caso de existir una sola campaña por año, se determina
el parámetro T0, que cuantifica la desviación del contenido anual en cloruros respecto de la media de
toda la serie y viene dado por la siguiente expresión:
T0 = ⎛⎜ f ( x ) − f
⎝
)
2
(mg/l)
donde f(x) es el valor de la concentración en cloruros para cada año y f el valor medio de la serie
completa. En consecuencia, el valor del parámetro T0 está atenuado y es de menor cuantía que el
obtenido a partir de las dos series de datos de valores extremos. Por este motivo, ha sido necesario
aplicar un factor de corrección, determinado empíricamente a partir de tres acuíferos estudiados que
disponían de puntos con dos o más valores anuales de concentración en ion cloruro. En ellos se calculó
el parámetro T a partir de las dos series de datos extremos y T0 con los de una sola serie (en este caso
para las mismas campañas estacionales). La correlación obtenida entre ambos resultados fue próxima a
1,5. Por lo que: T = 1,5 × T0
Al igual que en para los parámetros anteriores, a cada rango de valores de T se le asigna un
valor numérico y un código de identificación (tabla 1).
Este parámetro (E) tiene por objeto expresar de manera simple la tendencia actual del proceso
intrusivo, que se cataloga en cinco tipologías según su comportamiento: estable, progreso moderado,
progreso importante, deterioro moderado y deterioro importante. Para su cálculo se parte de la serie
temporal de la concentración anual media de cloruros en el acuífero, obtenida a partir de la matriz
básica de datos depurados. En primer lugar, se determina la situación actual, según lo establecido en el
capítulo de definiciones conceptuales, y posteriormente la situación precedente. A continuación se
calcula el parámetro E mediante la expresión:
Cln
E=
Cln −1
siendo Cln el valor medio de la concentración de cloruros en todo el acuífero para la situación actual
(últimos cuatro años) y Cln-1 el valor medio de la concentración de cloruros para la situación
precedente (cuatro años inmediatamente anteriores). Dicho valor medio se obtiene, a su vez, a partir
del contenido medio en cloruros en cada punto para el periodo considerado. Para cada rango de
valores de E se asigna un valor numérico y un código de identificación (tabla 1), en este caso
coincidentes.
Como ya se ha expuesto, los parámetros S, I y T permiten diferenciar cinco rangos a los que se
asigna un valor numérico comprendido entre 0 y 4 (tabla 1). El valor 0 correspondería a la situación
más favorable y el valor 4 a la más desfavorable. Los códigos correspondientes a estos valores son: W
(Without), L (Low), M (Medium), H (High) y E (Extreme). Por su parte, el parámetro E define el
comportamiento de la salinidad general del sistema en los cuatro últimos años, expresada mediante un
código numérico. Como resultado, el índice SITE puede ser representado mediante dos formas, una
alfanumérica, que da información específica sobre cada uno de los parámetros que intervienen en el
proceso de salinización, y otra exclusivamente numérica, cuya resultante es una única cifra que
informa sobre la situación global del sistema. La composición del código SITE alfanumérico se realiza
mediante la agregación de los cuatro códigos correspondientes a cada parámetro. Así, un acuífero en la
situación más favorable tendría un código WWW-2 y un acuífero en la situación más desfavorable
respondería al código EEE2. Para la determinación numérica del índice se utiliza un algoritmo que
permite ponderar los cuatro parámetros que lo conforman, de acuerdo con la siguiente expresión:
3S + 3I + T + E
ValorSITE =
30
El resultado de este algoritmo, cuyo valor máximo (1,0) correspondería a la situación más
negativa, permite establecer cinco categorías del estado de la intrusión marina, o de la presión sufrida
por un acuífero a causa de este proceso (tabla 2). Finalmente, el índice SITE completo se obtiene de
unir el código y el valor SITE. Así, un acuífero cuyo índice SITE sea MML1(0,47) corresponde a un
acuífero con intrusión marina zonal, de intensidad media y variaciones estacionales poco
significativas, en el que se ha registrado un moderado incremento de la salinidad. La valoración
general arroja el resultado de que la presión debida a intrusión marina en el acuífero es MEDIA.
También, el código de temporalidad T apuntaría a que se trata de un acuífero bastante inercial, con
cierta dificultad para ser recuperado. La aplicación del método a varios acuíferos mediterráneos puede
verse en García-Menéndez et al. (2012).
6. CONCLUSIONES
El método propuesto para la caracterización de la intrusión marina en los acuíferos costeros, que
ha sido aplicado a varios acuíferos costeros mediterráneos, permite sintetizar los principales rasgos de
este tipo de procesos de una manera objetiva y fácil de elaborar, al ofrecer información esencial para
poder evaluar, por un lado, la situación general de un acuífero y, por otro, discriminar entre acuíferos
distintos. Dichas cualidades hacen al método SITE especialmente eficaz para el cumplimiento de las
Directivas 2000/60/EC y 2006/118/EC, cuyo fin último es alcanzar el buen estado cuantitativo y
cualitativo de las masas de agua subterránea. La información obtenida de esta forma se ajusta a los
requerimientos exigidos para conocer de forma estandarizada la situación de las masas de agua
subterránea, permitiendo ser comparadas entre ellas de manera que se puedan establecer prioridades de
actuación, así como las estrategias de recuperación más adecuadas, que optimicen los recursos
humanos y financieros disponibles, y orienten sobre la mejora de su gestión mediante la identificación
objetiva de los problemas.
Agradecimientos
Los autores de este artículo agradecen a Juan de Dios Gómez Gómez, Manuel Varela Sánchez y Isaac Sánchez Navarro que,
como responsables de la Actividad 7 de la Encomienda de Gestión entre la Dirección General del Agua y el Instituto
Geológico y Minero de España, permitieron la aplicación de esta metodología en dichos trabajos, lo que ayudó a avanzar en
el desarrollo de la misma. También se agradece a David Pulido Velázquez las sugerencias y comentarios realizados.
REFERENCIAS
Anexo 1
A1-2. Determinación de la intrusión marina
en tres acuíferos costeros mediterráneos
mediante el método SITE
A5
Anexo 1 Índice SITE
A6
DETERMINACIÓN DE LA INTRUSIÓN MARINA EN TRES
ACUÍFEROS COSTEROS MEDITERRÁNEOS MEDIANTE EL
MÉTODO SITE
(1) Unidad de Valencia. Instituto Geológico y Minero de España. c/ Cirilo Amorós, 42, entreplanta, 46009
Valencia, España, [email protected].
(2) Instituto Universitario de Plaguicidas y Aguas. Universitat Jaume I, Castellón, España
RESUMEN
Se aplica el índice SITE a los acuíferos de la Plana de Vinaroz, Plana de Oropesa-Torreblanca y Plana
Sagunto, localizados en la costa mediterránea española (Castellón y Valencia), para determinar la presión que
sufren a causa de la intrusión marina. Los parámetros evaluados, basados exclusivamente en la concentración
del ion cloruro de las aguas subterráneas, son la superficie afectada por la salinización (S), la intensidad de la
misma (I), la temporalidad o variación estacional (T) y la evolución a medio-largo plazo (E). El método
utilizado, que es de fácil aplicación y proporciona información realista, discriminatoria y comprensible,
permite conocer de una forma codificada y estandarizada el estado de los procesos de intrusión marina en las
masas de agua subterránea, por lo que son comparables entre sí, de manera que se pueden establecer
prioridades de actuación, optimizar los recursos y orientar adecuadamente la mejora de su gestión mediante la
identificación objetiva de los problemas. Los resultados indican que la presión por intrusión marina es baja en
la Plana de Vinaroz, alta en la Plana de Oropesa-Torreblanca y media en la Plana de Sagunto.
ABSTRACT
In order to evaluate the pressure due to seawater intrusion, SITE index is applied to three coastal detritic
aquifers (Vinaroz Plain, Oropesa-Torreblanca Plain and Sagunto Plain) located in Castellon and Valencia
(Spain). The criteria assessed, based on the groundwater chloride concentration, are the following: spatial
affection (S), intensity of salinization (I), temporal variability or seasonality (T) and evolution to medium-
long term (E). The SITE index is easy to be applied and provides realistic, discriminatory and understandable
information. This information is very useful to know the seawater intrusion process of coastal groundwater
bodies in a codified and standarized way, allowing them to be compared with each others. By the objective
allocation of the aquifers problems is possible to prioritize the different interventions, optimize the resources
and improve their management. The results show that the pressure of marine intrusion is low in Vinaroz
Plain, high in Oropesa-Torreblanca Plain and media in Sagunto Plain.
1. INTRODUCCIÓN
Alcanzar el buen estado cuantitativo y cualitativo de las masas de agua subterránea es uno de los
principales objetivos de las Directivas 2000/60/CE y 2006/118/CE, que establecen un primer horizonte
para la consecución de sus fines en el año 2015. En consecuencia, y para poder cumplir con estos
requerimientos, es necesario contar con metodologías fiables y basadas en un enfoque común,
orientadas a establecer criterios comparables para la evaluación del estado químico de las masas de
agua subterránea, determinar las tendencias significativas o sostenidas al aumento de sus componentes
y definir los puntos de partida de las inversiones de tendencia. Dentro de este contexto, los acuíferos
costeros afectados por procesos de intrusión marina muestran una casuística muy variada,
condicionada tanto por factores intrínsecos (tipo de acuífero, tamaño, parámetros hidráulicos y
condiciones hidrostáticas) como por factores extrínsecos (recarga natural e inducida y bombeos).
Como resultado de esta influencia múltiple, la caracterización del proceso intrusivo es compleja, y
sólo en aproximaciones simples y mediante el empleo de criterios básicos puede ser abordada de
forma práctica y universal. No obstante, el método elegido debe reflejar las principales características
del proceso de salinización y proporcionar información cuantitativa y cualitativa sobre la afección
espacial, la intensidad del proceso, la estacionalidad y la evolución de la intrusión. En estas premisas
se fundamenta el método SITE (Ballesteros et al., 2012), cuyas primeras bases fueron establecidas por
B.J. Ballesteros (2008) en un documento interno del IGME. Más tarde, se aplicó una primera versión
del índice a varios acuíferos costeros españoles (IGME-DGA, 2010; Morell et al., 2010), mientras que
en trabajos posteriores se perfilaron y terminaron de establecer los procedimientos de estimación de
los parámetros que conforman definitivamente el método. En este documento se presenta la aplicación
de dicho índice a tres acuíferos costeros localizados en las provincias de Castellón y Valencia (figura
1): Plana de Vinaroz, Plana de Oropesa-Torreblanca y Plana de Sagunto, afectados por procesos de
intrusión marina, en mayor o menor grado, a causa de la explotación de sus recursos hídricos.
Figura 1. Situación de los acuíferos de la Plana de Vinaroz, la Plana de Plana de Oropesa-Torreblanca y la Plana de Sagunto.
Estos tres acuíferos, de naturaleza detrítica y tipo multicapa, con superficies entre los 90 km2 (P.
de Oropesa-Torreblanca) y los 129 km2 (P. de Sagunto), están formados por conglomerados y gravas
distribuidos en una matriz limo arenosa. Su comportamiento puede asimilarse al de un acuífero libre
en el que la presencia de niveles de baja permeabilidad puede provocar la existencia de situaciones
más o menos locales de semiconfinamiento. La heterogeneidad espacial da lugar a un amplio rango de
valores de porosidad y conductividad hidráulica aunque, en general, presentan alto rendimiento en las
numerosas captaciones existentes, cuyos caudales superan en algunos puntos los 100 l/s.
Los recursos de estos sistemas hídricos están entre los 43,4 hm3/año de la Plana de Sagunto y
los 63,9 hm3/año de la Plana de Vinaroz, mientras que las extracciones por bombeo están entre los
39,0 hm3/año de la Plana de Oropesa-Torreblanca y los 55,3 hm3/año de la Plana de Sagunto (IGME-
DGA, 2010, Ballesteros et al., 2004) (tabla 1).
Tabla 1. Estado cuantitativo de la Plana de Vinaroz, la Plana de Oropesa-Torreblanca y la Plana de Sagunto para el periodo
1995-2001 (1) y 2002-03 (2).
2. FUNDAMENTOS DEL MÉTODO SITE
Como primer paso, se debe elaborar una matriz básica que recoja para cada punto sus
características esenciales (coordenadas, cota, profundidad, etc.) y las medidas de las que dispone. De
este modo, y con el fin de homogenizar la información tanto espacial como temporalmente, se facilita
el tratamiento y depuración de los datos iniciales. A partir de dicha matriz se obtendrán las matrices de
datos específicas para los diferentes parámetros SITE. El método requiere de la definición previa de
algunos conceptos utilizados para su determinación, tales como el valor de referencia que es el fondo
geoquímico del acuífero referido a ión cloruro o el contenido por encima del cual se considera que el
agua de un sistema está afectado por la intrusión marina, la superficie total del acuífero, expresada en
km2, y la superficie afectable, que es aquella para la cual se estima que pueden darse procesos
intrusivos y se establece por defecto a la situada a menos de 10 km de la línea de costa. Por otra parte,
se define la situación actual como la correspondiente a los últimos cuatro años, de manera que los
valores utilizados en los cálculos sean los resultantes de la media de los obtenidos anualmente en cada
punto a lo largo de dicho periodo, y la denominada situación precedente determinada de igual forma
pero con el periodo correspondiente a los cuatro años previos. El desarrollo del método se expone de
forma más detallada y pormenorizada en Ballesteros et al. (2012).
Para el cálculo del parámetro S se elabora el mapa de isoconcentración del ión cloruro
correspondiente a la situación actual y se define la isolínea del valor de referencia (Vr), que para los
acuíferos seleccionados se ha considerado la de 250 mg/l (valor por defecto). También se determinan
las isolíneas correspondientes a 0,6Vr (150 mg/l), 2Vr (500 mg/l), 4Vr (1.000 mg/l), 6Vr (1.500 mg/l) y
8Vr (2.000 mg/l), que servirán para el cálculo del próximo parámetro. El valor de S se calcula según:
S > 250
S=
S total
donde S>250 es la superficie ocupada por contenidos de ión cloruro superiores a 250 mg/l y Stotal la
superficie total del acuífero, ya que no ha sido necesario aplicar la corrección de superficie. Para cada
rango de valores de S se asigna un valor, o índice numérico, y un código de identificación según el
baremo de la tabla 2. Los códigos obedecen a las iniciales en lengua inglesa de Without (W), Low (L),
Medium (M), High (H) y Extreme (E).
Para la elaboración del índice SITE se han utilizado los datos disponibles de la Red de Intrusión
del IGME, operativa con datos semestrales desde 1977 hasta el 2001, procedentes de 59 puntos en la
Plana de Vinaroz, 33 puntos en la Plana de Oropesa-Torreblanca y 30 puntos en la Plana de Sagunto.
A partir de ellos se ha definido la situación más reciente, que corresponde al periodo 1998-2001. Con
los valores medios de la concentración del ión cloruro en cada punto en dicho periodo se elaboran los
correspondientes mapas de isocloruros (figura 2) y se determinan las superficies correspondientes a los
rangos definidos (tabla 3).
Parámetro Rango Caracterización Valor Código
0,00 – 0,10 Sin intrusión o puntual 0 W
0,11 – 0,25 Localizada o baja 1 L
S 0,26 – 0,50 Zonal o media 2 M
0,51 – 0,75 General o alta 3 H
0,76 – 1,00 Total o extrema 4 E
< Vr < 250 Nula 0 W
I Vr-2Vr 251 - 500 Baja 1 L
(Vr=250 mg/l) 2Vr-4Vr 501-1.000 Media 2 M
4Vr-8Vr 1.001-2.000 Alta 3 H
> 8Vr > 2.000 Extrema 4 E
< 0,1Vr < 25 Sin estacionalidad 0 W
T 0,1Vr-0,2Vr 26 - 50 Estacionalidad baja 1 L
(Vr=250 mg/l) 0,2Vr-0,4Vr 51-100 Estacionalidad media 2 M
0,4Vr-0,8Vr 101-200 Estacionalidad elevada 3 H
> 0,8Vr > 200 Estacionalidad extrema 4 E
< 0,85 Recuperación importante -2 -2
0,85 – 0,97 Recuperación moderada -1 -1
E 0,98 – 1,02 Estable 0 0
1,03 – 1,15 Deterioro moderado 1 1
> 1,15 Deterioro importante 2 2
De este modo, el valor del parámetro S obtenido para la Plana de Viñaroz es igual a 0,20, lo que
indica que el acuífero está sometido a una intrusión de carácter localizado, y se le asigna un valor de 1
y el código L. Por su parte, en la Plana de Oropesa-Torreblanca el parámetro S es igual a 0,96, lo que
corresponde a un acuífero sometido a una intrusión de carácter extremo, y se le asigna un valor de 4 y
el código E. Finalmente, en la Plana de Sagunto el parámetro S es igual a 0,59, es decir un acuífero
sometido a una intrusión de carácter general o alta, con un valor de 3 y un código M.
El parámetro Intensidad persigue obtener una aproximación a la salinidad media del acuífero,
para lo cual se parte del mapa de isocloruros correspondiente a la situación actual elaborado para la
obtención del parámetro anterior. A partir de él se calculan las superficies (Si) comprendidas entre
isolíneas consecutivas y a cada una de ellas se le asigna el valor medio de dichas isolíneas (Cli). A la
superficie de contenido inferior a 0,6Vr (mg/l) se le otorga, por convención, un valor medio de 0,4Vr
(mg/l), mientras que a la superficie que supere la isolínea más elevada se le atribuye el valor de esta
isolínea más la diferencia existente entre dicha isolínea y la que correspondería ser la inmediatamente
superior. El parámetro I resulta de la expresión:
I=
∑S i ( > 250 ) ⋅ Cl i > 250
S ( > 250 )
2
donde Si es la superficie en km entre la isolínea i y su inmediatamente superior, y Cli la concentración
media de cloruros en mg/l asignada a dicha superficie. Para cada rango de valores de I se asigna el
valor numérico y el código de identificación de la tabla 2.
Figura 2. Mapas del valor medio de isocloruros para el periodo 1998-2001 de los acuíferos de la Plana de Vinaroz
(izquierda), Plana de Oropesa-Torreblanca (centro) y Plana de Sagunto (derecha).
Rango Cl
- Supeficie (km2)
(mg/l) Plana de Vinaroz P. de Oropesa-Torreblanca Plana de Sagunto
<150 72,90 0 25,64
84,65 3,15 53,4
150-250 11,75 3,15 27,76
250 - 500 9,07 15,01 51,78
500 - 1.000 7,38 37,82 23,63
21,01 86,57 75,41
1.000-1.500 0,68 33,71 0
1.500-2.000 3,89 0,03 0
S>250 21,01 86,57 75,41
S>total 105,66 89,72 128,81
Parámetro S 0,20 0,96 0,59
Tabla 3. Superficies de los rangos de contenido en ión cloruro y valor del parámetro S en los acuíferos de la Plana de
Vinaroz, Plana de Oropesa-Torreblanca y Plana de Sagunto (periodo 1998-2001).
Con el valor de referencia considerado (250 mg/l), el valor medio de contenido en cloruros para
cada rango, la superficie correspondiente a dicho rango y la sumatoria del producto de ambos factores
para cada uno de los tres acuíferos se expone en la tabla 4. Con estos datos el valor de I queda
establecido para la Plana de Vinaroz en 248 mg/l, lo que indica que la salinidad del acuífero puede
considerarse como nula (índice numérico: 0 y código: W), para la Plana de Oropesa-Torreblanca en
856 mg/l, lo que indica que la salinidad del acuífero es media (índice numérico: 2 y código: M) y para
la Plana de Sagunto en 351 mg/l, lo que supone una salinidad del acuífero baja (índice numérico: 1 y
código: L)
Tabla 4. Rangos de cloruros, valor medio de cada rango, sumatorios Si.Cli y valor del parámetro I en la Plana de
Vinaroz, Plana de Oropesa-Torreblanca y Plana de Sagunto (periodo 1998-2001).
5. CÁLCULO DE LA TEMPORALIDAD O ESTACIONALIDAD (PARÁMETRO T)
En la Plana de Vinaroz se han seleccionado 27 puntos que poseen dos medidas anuales de la
concentración de cloruros en una serie temporal de 14 años (figura 3), para los que el valor de T ha
resultado ser de 97 mg/l, lo que indica una situación caracterizada como de estacionalidad media
(índice numérico: 2 y código: M). En la Plana de Oropesa-Torreblanca, con un total de 19 puntos con
dos datos anuales en una serie de 20 años (figura 4), el valor de T es igual a 220 mg/l, expresión de
una estacionalidad extrema (índice numérico: 4 y código: E). En cuanto a la Plana de Sagunto, el
parámetro T calculado a partir de un total de 11 puntos con datos en 13 años es igual a 60 mg/l (figura
5), lo que indica también una estacionalidad media (índice numérico: 2 y código: M).
El parámetro E expresa de manera simple la tendencia temporal del proceso intrusivo. Para
hallarlo se determinan las situaciones actual y precedente según lo establecido en las definiciones
conceptuales y, posteriormente, se calcula E mediante la expresión:
Cln
E=
Cln −1
En los tres acuíferos se ha comparado la concentración del ión cloruro en el periodo 1998-2001
con la calculada en el periodo 1994-1997 (figura 6). Los resultados arrojan valores para el acuífero de
la Plana de Vinaroz de 1,04, que caracteriza una situación de deterioro moderado (índice numérico y
código: 1), para el acuífero de la Plana de Oropesa-Torreblanca de 1,01 (índice numérico y código: 0),
que caracteriza una situación estable, y en el acuífero de la Plana de Sagunto de 1,09, que caracteriza
una situación de un moderado aumento de la concentración de cloruros (índice numérico y código: 1).
El índice SITE representado de forma alfanumérica da una información específica sobre cada
uno de los parámetros que intervienen en el proceso de salinización, y resulta de la agregación de los
cuatro códigos correspondientes a cada parámetro. Por su parte, la forma numérica refiere de manera
global la situación del acuífero, y se obtiene de la siguiente expresión:
3S + 3I + T + E
ValorSITE =
30
Cl (mg/l)
1989 408 499 91 350
-
1990 398 525 127 300
Figura 3: Valores y evolución de las concentraciones anuales del ión cloruro en el acuífero de la Plana de Vinaroz: máxima
(rojo), mínima (azul) y diferencia (verde)
Cl- (mg/l)
AÑO
Valor mínimo_medio Valor máximo_medio Tx 1400
Figura 4: Valores y evolución de las concentraciones anuales del ión cloruro en el acuífero de la Plana de Oropesa-
Torreblanca: máxima (rojo), mínima (azul) y diferencia (verde)
Figura 5: Valores y evolución de las concentraciones anuales del ión cloruro en el acuífero de la Plana de Sagunto: máxima
(rojo), mínima (azul) y diferencia (verde)
400 950
E= 1,01
375 925
E= 1,04 E= 1,09
895
350 900
335 882
323 321
325 875
Cl - (mg/l)
Cl - (mg/l)
295
300 850
275 825
250 800
225 775
200 750
Plana de V inaroz Plana de Sagunto Plana de Orope s a-Torre blanca
Figura 6: Concentración de cloruros en el periodo 1994-1997 (situación precedente) en gris claro y en el periodo 1998-
2001(situación actual) en gris oscuro en la Plana de Vinaroz, la Plana de Sagunto y la Plana de Oropesa-Torreblanca
El cociente 30 es el valor máximo del numerador, por lo que un acuífero totalmente salinizado
tendría un valor SITE de 1 mientras que un acuífero sin afección alguna tendría un valor SITE de 0, lo
que permite establecer diferentes categorías del estado de la intrusión marina (tabla 5). Por último el
Índice SITE general es el resultado de unir el código alfanumérico y el valor numérico del índice.
Los valores numéricos y códigos del índice SITE que caracterizan a los tres acuíferos se
incluyen en la tabla 6. Así, la agregación de los cuatro códigos obtenidos para la Plana de Vinaroz
arrojan un índice alfanumérico de LWM 1 y un valor SITE de 0,20, por lo que la expresión final del
índice es LWM 1(0,20). En la Plana de Oropesa-Torreblanca, el índice alfanumérico es EME 0 y de la
aplicación del algoritmo de cálculo se obtiene el valor SITE de 0,73, por tanto la expresión final del
índice es EME 0 (0,73). Por último, los valores numéricos y códigos que caracterizan el acuífero de la
Plana de Sagunto, dan lugar al índice alfanumérico HLM 1 y el valor SITE de 0,5, con una expresión
final del índice HLM 1 (0,50).
Criterio Valor Índ_num Código Caracterización Valor Índ_num Código Caracterización Valor Índ_num Código Caracterización
SITE 0,20 LWM1 Intrusión BAJA 0,73 EME0 Intrusión ALTA 0,50 HLM1 Intrusión MEDIA
Tabla 6. Valores alfanuméricos y códigos del índice SITE en los acuíferos de la Plana de Vinaroz, Plana de Oropesa-
Torreblanca y Plana de Sagunto
8. CONCLUSIONES
Los resultados obtenidos mediante la aplicación del índice SITE han permitido determinar de
una forma estandarizada la situación y el grado de afección al que se encuentran sometidos los
acuíferos analizados a causa de los procesos de intrusión marina, y han aportado una información
objetiva de los mismos desde los puntos de vista espacial y temporal. Así, la expresión del índice
obtenido para el acuífero de la Plana de Vinaroz, LWM 1 (0,20), caracteriza a un sistema con intrusión
marina local, de intensidad nula y variaciones estacionales medias, que presenta una tendencia
moderada hacia la intensificación de la salinización. El resultado final indica que la presión debida a la
intrusión marina en el sistema es BAJA. Por su parte, en el acuífero de la Plana de Oropesa-
Torreblanca el índice obtenido, EME 0 (0,73), corresponde a un acuífero con intrusión marina total, de
intensidad media y variaciones estacionales muy significativas, con una estabilización de los procesos
de salinización. La valoración general refleja que la presión debida a la intrusión marina en el acuífero
es ALTA. Finalmente, en el acuífero de la Plana de Sagunto el índice obtenido, HLM1(0,50),
caracteriza a un acuífero con intrusión marina general, de intensidad baja y variación estacional media,
en el que se detecta un moderado incremento de los procesos salinos, por lo que la valoración general
de la presión sufrida por el sistema debida a la intrusión marina es MEDIA.
Como conclusión final se puede afirmar que la aplicación del índice SITE para la
caracterización de la intrusión marina a tres acuíferos costeros mediterráneos revela que esta
metodología es capaz de sintetizar los principales rasgos de este tipo de procesos de una manera
objetiva y de fácil elaboración, ofreciendo información esencial para poder evaluar, por un lado, la
situación general de un acuífero y, por otro, discriminar entre acuíferos distintos.
REFERENCIAS
Ballesteros, B.J., Barba-Romero Muñoz, J. y Navarro Odriozola, J.O. (2004). Estudio piloto de
caracterización adicional de masas de agua subterránea aplicado a la U.H. (08.21) Plana de
Sagunto. Instituto Geológico y Minero de España.
Ballesteros, B.J. (2008): Tipología y caracterización de la intrusión marina en acuíferos litorales.
Método SITE. Primeras bases. Informe inédito. Instituto Geológico y Minero de España.
Ballesteros B.J, Morell, I., García, O. y Renau-Pruñonosa, A. (2012): Propuesta de un nuevo índice
para la caracterización de la intrusión marina. El método SITE. IV Congreso Internacional
sobre la tecnología de la intrusión de agua de mar en acuíferos costeros. TIAC´12. Alicante.
IGME-DGA. (2010): Encomienda de gestión para la realización de trabajos científico-técnicos de
apoyo a la sostenibilidad y protección de las aguas subterráneas. Actividad 7:
Establecimiento de indicadores de intrusión marina y cálculo de volúmenes ambientales al
mar. Documento interno. Instituto Geológico y Minero de España-Dirección General del
Agua.
Morell, I., Ballesteros, B.J., Renau, A, Gómez, J.D. (2010): Application of a new seawater intrusion
index: SITE method. Keynote in Global Forum on Salinization and Climate Change. FAO-
Valencia University, Valencia.
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
Anexo 2
Datos de entrada (inputs) del
modelo matemático
A7
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
A8
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
Los datos mensuales de precipitación comprendido entre 1973 y 2005 (ambos inclusive)
se han obtenido de los estudios realizados por González - Hidalgo (2008 y 2010), la estación
meteorológica utilizada en dichos estudios está ubicada en un área próxima, situada en las
cercanías del polígono industrial del Campaner, en la Asociación de Viveristas de Ágrios
(AVASA), término municipal de Alcalà de Xivert (tabla A-1).
A9
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
MESES TOTAL
AÑO
(mm)
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Sept. Octubre Nov. Dic.
1973 6,2 2,9 11,1 25,3 13 51,3 40,3 49,1 37,3 37,8 10,7 46,8 331,8
1974 1,4 9,6 119,9 114,8 47,7 16 14,8 49 18,1 18,6 2,5 3 415,4
1975 3,1 111 63,2 14 116,4 49,9 5,8 32,4 98,1 41,8 22,8 115,3 673,8
1976 2 25,9 7,8 68,3 35,7 19,5 52,1 74,9 24 55,8 4,7 118 488,7
1977 138,7 9,5 13,2 46,5 144,7 59,9 49,2 19 72,1 42,6 74,2 95,9 765,5
1978 9,3 17,5 22,8 35,1 46,5 15,6 3,3 7,4 34,6 33 50,2 21,9 297,2
1979 123,3 5,3 13,2 24,3 22 17,4 6,6 16,4 47,4 68,5 6,6 1,3 352,3
1980 32,4 90,5 12,4 49 76,4 19,2 14,8 34,7 7,5 5,2 42,9 0,9 385,9
1981 1,6 77,4 8,5 227,8 15,2 122,6 34,7 72,8 29,5 11,1 0,1 2,6 603,9
1982 148,2 68,7 118,2 43,5 50,3 1,6 30,4 72,3 109,4 221,2 26,1 4 893,9
1983 0,1 9,4 16,8 0,7 10,8 32,6 0,2 92,1 0,2 11,2 192,7 5,3 372,1
1984 0,3 50,1 18,6 36 109,8 4,1 0,1 17,4 21,3 16,9 141,4 0,4 416,4
1985 15,6 12 11,6 18,8 47,2 1,9 0,1 1,4 92,2 53,7 105,1 16,4 376
1986 0,5 5,5 15,9 23,3 10,6 0,3 21,6 4,1 107,2 150,4 77 9,2 425,6
1987 152,7 30,4 9,6 0,1 15,6 3,6 45,4 14,3 11,1 222,7 137,3 40,5 683,3
1988 98 1,9 1,8 186,9 79,7 75,4 1 2,1 29,1 40,5 177,5 0,1 694
1989 32,3 56,8 153,1 42,4 50 25,9 9,6 44,8 308,1 19,9 171,4 198 1112,3
1990 97,1 0,2 21,3 34,1 13,6 64,6 9,1 25,6 92,9 353,1 14,6 18,1 744,3
1991 62,9 182,6 76,1 88,4 24,1 18,6 14,1 1,7 49,7 39,7 79,9 187,5 825,3
1992 32,6 52,2 3,6 20,3 39,6 69,4 13,6 11,6 20,1 45,2 0,2 118,2 426,6
1993 0,3 89,4 45,8 18,2 19,5 0,1 1,9 2,2 101,6 64,6 56,2 0,1 399,9
1994 0,1 20,8 2,3 34,4 12,3 4,4 8,1 35,1 55,7 157,2 78,9 20,9 430,2
1995 1,5 15,6 2,6 26,8 7,9 6,6 12,9 23,9 61,6 48,9 31,3 124 363,6
1996 132,1 93,4 9,5 24,6 31,4 3 6,2 33,4 66,5 13,2 27,6 119,6 560,5
1997 81 0,3 0,1 78,1 13,1 29,1 15,7 34,3 40,1 14,2 24,1 78,2 408,3
1998 168,8 12,7 0,1 20,1 39 0,1 0,1 22,1 36,4 18,1 18,1 164,4 500
1999 47,8 4,3 69 75,8 6,8 10,4 25,6 5,1 82,8 84,9 28,6 0,1 441,2
2000 54,4 0,1 25,8 26,5 38,2 47,7 29,5 0,1 18,7 283,5 17 30,7 572,2
2001 24,4 21,1 15,4 27,6 42,4 2,9 18,9 1,5 54,6 36,3 56,1 41,4 342,6
2002 60,8 0,1 26,3 77,3 186,8 16,2 11,9 67,9 74,7 16,8 25,6 17,1 581,5
2003 0,1 34,3 49,4 102,4 127,3 9,3 6,8 82,1 15,4 96,3 269,5 33,3 826,2
2004 0,1 74,4 126,5 104,2 123,1 29,9 5,1 65,8 33,2 18,8 5,8 29,1 616
2005 0,1 123,6 9,3 30,6 27,2 0,1 0,1 49,6 35,8 52,3 174,1 11 513,8
2006 137,4 86,2 7,2 5,4 62,6 4 0,2 1,8 102,4 30,2 83,2 4,4 525
2007 39,8 42,4 52,2 237 19,6 11,4 0 0 94,6 121 0 94,6 712,6
2008 0 56 14 5 142 26 58 21 79 114 36,8 56 607,8
2009 56 9 49 31,6 5,6 0 9 19,4 233 10 3 47 472,6
A10
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
MESES TOTAL
AÑO
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Sept. Octubre Nov. Dic. (mm)
2010 47,8 4,3 69 75,8 6,8 10,4 25,6 5,1 82,8 84,9 28,6 0,1 441,2
2011 54,4 0,1 25,8 26,5 38,2 47,7 29,5 0,1 18,7 283,5 17 30,7 572,2
2012 24,4 21,1 15,4 27,6 42,4 2,9 18,9 1,5 54,6 36,3 56,1 41,4 342,6
2013 60,8 0,1 26,3 77,3 186,8 16,2 11,9 67,9 74,7 16,8 25,6 17,1 581,5
2014 0,1 34,3 49,4 102,4 127,3 9,3 6,8 82,1 15,4 96,3 269,5 33,3 826,2
2015 0,1 74,4 126,5 104,2 123,1 29,9 5,1 65,8 33,2 18,8 5,8 29,1 616
2016 0,1 123,6 9,3 30,6 27,2 0,1 0,1 49,6 35,8 52,3 174,1 11 513,8
2017 137,4 86,2 7,2 5,4 62,6 4 0,2 1,8 102,4 30,2 83,2 4,4 525
2018 39,8 42,4 52,2 237 19,6 11,4 0 0 94,6 121 0 94,6 712,6
2019 0 56 14 5 142 26 58 21 79 114 36,8 56 607,8
2020 56 9 49 31,6 5,6 0 9 19,4 233 10 3 47 472,6
2021 47,8 4,3 69 75,8 6,8 10,4 25,6 5,1 82,8 84,9 28,6 0,1 441,2
2022 54,4 0,1 25,8 26,5 38,2 47,7 29,5 0,1 18,7 283,5 17 30,7 572,2
2023 24,4 21,1 15,4 27,6 42,4 2,9 18,9 1,5 54,6 36,3 56,1 41,4 342,6
2024 60,8 0,1 26,3 77,3 186,8 16,2 11,9 67,9 74,7 16,8 25,6 17,1 581,5
2025 0,1 34,3 49,4 102,4 127,3 9,3 6,8 82,1 15,4 96,3 269,5 33,3 826,2
2026 0,1 74,4 126,5 104,2 123,1 29,9 5,1 65,8 33,2 18,8 5,8 29,1 616
2027 0,1 123,6 9,3 30,6 27,2 0,1 0,1 49,6 35,8 52,3 174,1 11 513,8
2028 137,4 86,2 7,2 5,4 62,6 4 0,2 1,8 102,4 30,2 83,2 4,4 525
2029 39,8 42,4 52,2 237 19,6 11,4 0 0 94,6 121 0 94,6 712,6
2030 0 56 14 5 142 26 58 21 79 114 36,8 56 607,8
2031 56 9 49 31,6 5,6 0 9 19,4 233 10 3 47 472,6
2032 47,8 4,3 69 75,8 6,8 10,4 25,6 5,1 82,8 84,9 28,6 0,1 441,2
2033 54,4 0,1 25,8 26,5 38,2 47,7 29,5 0,1 18,7 283,5 17 30,7 572,2
2034 24,4 21,1 15,4 27,6 42,4 2,9 18,9 1,5 54,6 36,3 56,1 41,4 342,6
2035 60,8 0,1 26,3 77,3 186,8 16,2 11,9 67,9 74,7 16,8 25,6 17,1 581,5
A11
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
A12
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
Figura A-1. Zonas de recarga correspondientes a Alcossebre (Alcalà de Xivert) ( 1973 – 2009)
A13
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
A14
Anexo 2 Recarga. Infiltración por precipitación (Pp) y retorno de riego (Rr) (1973 – 2009)
ZONA INFILTRACIÓN
1 Pp+Rr (1973-2009)
2 Pp+Rr (1973-2003)
3 Pp+Rr (1980-2009)
4 Pp+Rr (1976-2005)
5 Pp+Rr (1976-2007)
6 Pp+Rr (1976-2009)
7 Pp+Rr (1980-2003)
8 Pp+Rr (1980-2005)
9 Pp+Rr (1980-2007)
10 Pp+Rr (1985-2003)
11 Pp+Rr (1985-2005)
12 Pp (1973-1995)
13 Pp+Rr (1973-2007)
14 Pp+Rr (1985-2007)
15 Pp+Rr (2004-2009)
16 Pp
17 Recarga 0 mm.
18 Pp+Rr (1980-1985)+ 1/2Rr(1986-2005)
19 Pp (1973-2004)
20 Pp+ Rr (1976-1980)+ 1/2 Rr(1981-2009)
21 Pp+Rr (1976-2003)+ 1/2Rr(2004-2009)
22 Pp+ 1/2 Rr(1980-2007)
A15
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 6
4.00
3.00
Nivel piezométrico (m)
2.00
1.00
0.00
-1.00
-2.00
-3.00
-4.00
-5.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99 sep-02 may-05 feb-08 nov-10
Tiempo
Punto de observación 7
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
-2.00
-3.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99
Tiempo
Punto de observación 8
4.00
Nivel piezométrico (m)
3.00
2.00
1.00
0.00
-1.00
-2.00
-3.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99 sep-02
Tiempo
Punto de observación 9
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
-2.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99
Tiempo
A17
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 11
4.00
Nivel piezométrico (m)
3.00
2.00
1.00
0.00
-1.00
-2.00
ene-72 ene-75 ene-78 ene-81 ene-84 ene-87 ene-90 ene-93 feb-96 feb-99
Tiempo
Punto de observación 12
3.00
Nivel piezométrico (m)
2.00
1.00
0.00
-1.00
mar-71 mar-74 mar-77 mar-80 mar-83 mar-86 abr-89 abr-92 abr-95 abr-98 abr-01
Tiempo
A18
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 13
1.50
Nivel piezométrico (m)
1.00
0.50
0.00
-0.50
-1.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
Punto de observación 14
3.00
Nivel piezométrico (m)
2.00
1.00
0.00
-1.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99
Tiempo
A19
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 17
2.00
1.00
0.00
-1.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99
Tiempo
Punto de observación 18
2.00
1.00
0.00
-1.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99
Tiempo
A20
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 19
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
-2.00
jul-72 jul-75 jul-78 ago-81 ago-84 ago-87 ago-90 ago-93 ago-96 ago-99
Tiempo
Punto de observación 20
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
-2.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97
Tiempo
A21
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 21
3
Nivel piezométrico (m)
-1
-2
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
Punto de observación 22
3.00
Nivel piezométrico (m)
2.00
1.00
0.00
-1.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
tiempo
A22
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 23
3.00
Nivel piezométrico (m)
2.00
1.00
0.00
-1.00
-2.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
Punto de observación 24
2
Nivel piezométrico (m)
-1
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
A23
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 25
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
Punto de observación 26
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
A24
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 27
2
Nivel piezométrico (m)
-1
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
Punto de observación 28
2
Nivel piezométrico (m)
-1
jul-72 abr-75 ene-78 oct-80 jul-83 mar-86 dic-88 sep-91 jun-94 mar-97 dic-99
Tiempo
A25
Anexo 2 Puntos de observación. Niveles piezométricos
Punto de observación 29
2.00
Nivel piezométrico (m)
1.00
0.00
-1.00
-2.00
-3.00
-4.00
jul-72 nov-75 feb-79 may-82 sep-85 dic-88 abr-92 jul-95 nov-98 feb-02 may-05
Tiempo
A26
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
Anexo 3
Balance hídrico anual
A27
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A28
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A29
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A30
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A31
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A3-1.4. Balance hídrico anual del sector Sur de la Plana de Oropesa – Torreblanca
(1973 – 2009)
A32
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A33
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A34
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A35
Anexo 3 Balance hídrico anual del modelo matemático (1973 – 2009)
A3-2.4. Balance hídrico anual del sector Sur de la Plana de Oropesa – Torreblanca
(2010 - 2035)
A36