Los acuíferos costeros mediterráneos, en su mayoría, presentan un grado de explotación muy elevado, que ha dado lugar a la presencia de procesos de salinización por intrusión marina. Si bien el proceso de intrusión marina puede ser hidrodinámicamente muy complejo, puede asumirse que debe existir relación directa entre los caudales de descarga subterránea al mar y la intrusión marina. 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). La metodología propuesta para el cálculo del VER 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 sea menor o igual a la intrusión estimada en régimen natural para el sistema (intrusión de referencia). Para su aplicación se propone seguir los siguientes pasos: FASE 1: Simular la situación del sistema en régimen natural para conocer la intrusión de referencia. Al modelo resultante se le denominará modelo de referencia en régimen natural (modelo RN). FASE 2: Delimitar las áreas afectadas por la intrusión marina e intentar disminuir su superficie mediante la reducción iterativa de los bombeos involucrados en el proceso, hasta que la superficie piezométrica obtenida presente una intrusión igual o menor a la estimada en régimen natural. El modelo de partida en esta fase se denominará modelo base en régimen influenciado (modelo RI). FASE 3: Cálculo del VER a partir del modelo RI y del último escenario contemplado en la Fase 2 (el que se asocia a una intrusión igual o menor a la estimada en régimen natural) El volumen ecológico de remediación (VER) se puede calcular para años húmedos (VERh), medios (VERm) y secos (VERs). La aplicación de esta metodología en la Plana de Oropesa – Torreblanca (MASub 080.110) pone de manifiesto que el avance en la intrusión marina (respecto a la existente en régimen natural), en cualquiera de los tres casos (año húmedo, medio o seco) es debida principalmente a la acción de los bombeos. En término medio, en los años húmedos el sector afectado por el avance de la intrusión es el Sur y el VERh correspondiente es de 1 hm3/año. En los años medios, el VERm es de 1,6 hm3/año y el proceso afecta únicamente al sector Sur. En los años secos el proceso afecta al sector Sur y al sector Central; en este caso el VERs sería de 2,8 y 0,9 hm3/año para cada sector respectivamente. Previamente a la aplicación de la metodología propuesta a la Plana de Oropesa - Torreblanca (MASub 080.110), ha sido necesario profundizar en el conocimiento de las características del acuífero y su funcionamiento. En cuanto a la intrusión marina de la Plana, la evolución del proceso y su situación actual se han estudiado a partir de la información recogida en las redes de control de intrusión y calidad de la Confederación Hidrográfica del Júcar (CHJ, www.chj.es). Se ha aplicado una metodología para valorar los impactos y la vulnerabilidad del sistema frente a la intrusión (índice SITE), concluyéndose que, a finales del siglo XX, el proceso afectaba a la totalidad del acuífero y la presión debida a la intrusión se clasifica como alta. El desarrollo de un modelo matemático de flujo subterráneo (MODFLOW), para el periodo 1973 – 2009, ha permitido evaluar las acciones que afectan al sistema, calibrar los parámetros y acciones que lo conforman y observar la evolución de la piezometría a lo largo de la simulación. Al final del periodo considerado, año 2009, se estima que la intrusión, que durante el año 1993 llega a afectar a dos tercios de la franja costera con un caudal de entrada de 2,3 hm3/año, se concentra únicamente en el sector Sur y se cuantifica en 0,05 hm3/año. Se ha simulado un escenario futuro plausible (2010 – 2035), definido a partir de las ordenaciones urbanísticas proyectadas a corto y medio plazo en la zona: el Plan General de Ordenación Urbana 2009 (PGOU Torreblanca, 2009) dentro del cual se desarrolla el Programa de Actuación Integrada (PAI) Doña Blanca Golf y el PAI Marina d´Or Golf (en los términos municipales de Cabanes y Oropesa). Esto ha permitido valorar el impacto de los cambios del uso del suelo en la piezometría de la zona, observándose que a partir del año 2018 los niveles piezométricos por debajo del nivel del mar son inexistentes.