Máster en Técnicas de Gestión del Medio Ambiente y el Territorio Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor Carlos Cabezas Rabadán Tutor: Dr. Josep Eliseu Pardo Pascual 2014 – 2015
Apr 07, 2017
Máster en Técnicas de Gestión del Medio Ambiente y el
Territorio
Análisis de la línea de costa y su
relación con los parámetros
morfológicos en playas de la Safor
Carlos Cabezas Rabadán
Tutor: Dr. Josep Eliseu Pardo Pascual
2014 – 2015
Resumen
Las playas son espacios naturales frágiles de gran importancia socioeconómica. El ser
humano condiciona con sus acciones el comportamiento de la costa, la cual se encuentra
en muchos lugares amenazada por procesos erosivos que conducen al estrechamiento e
incluso desaparición de las playas, afectando a la actividad turística y de ocio que en ellas
se desarrolla.
La metodología seguida para obtener la posición de la línea de costa a partir de imágenes
Landsat cuenta con un elevado grado de precisión y automatización, permitiendo analizar
posteriormente con ayuda de herramientas SIG la evolución seguida por las playas,
esencial para poder llevar a cabo una gestión integrada y adecuada de las zonas costeras.
Al mismo tiempo se intentó comprobar la existencia de relaciones entre la variación de la
línea de costa a lo largo de un año y las características morfológicas de la playa.
Se determinó la existencia de relaciones entre la variación de la costa, la pendiente de la
playa y la textura de sus sedimentos. La obtención de la posición de la línea de costa
permite por tanto no solo conocer los patrones evolutivos de las playas, sino también
estimar algunas de sus características morfológicas como la textura y la pendiente.
Palabras Clave: Extracción de la línea de costa, evolución de la línea de costa, procesos
costeros, gestión de playas, herramientas SIG, imágenes Landsat
Abstract
Beaches are fragile natural areas with great socio-economic importance. Human acts
define the evolution of the coast, which in several cases is endangered by erosion
processes. The processes lead to diminishment or even disappearing of the beaches,
affecting tourism and leisure activities taking place in the given area.
The methodology used to obtain the position of the coastline from Landsat images works
with high level of accuracy and automatization allowing later analysis of the evolution of
the beach with help of GIS tools, which is essential for consistent and adequate coastal
zone management.
Simultaneously, we wanted to verify the existence of relations between annual changes
of the coastline and the morphologic characteristics of the beach.
The conclusion was that there is an existing relation between changes in the shoreline
position, slope of the beach and the texture of its sediments. Obtaining the position of the
coastline therefore allows us not only to study the evolutionary patterns of the beaches,
but also to estimate some of their morphological characteristics such as texture and slope.
Keywords: Coastline extraction, shoreline evolution, coastal processes, beach
management, GIS tools, Landsat images.
Índice 1. Introducción .............................................................................................................. 1
1.1 Importancia de las playas ................................................................................... 1
1.2 Utilidad del análisis de la línea de costa ............................................................ 3
1.3 Alternativas existentes ....................................................................................... 4
1.4 Extracción de la línea de costa con precisión subpíxel a partir de imágenes
Landsat TM .................................................................................................................. 5
2. Objetivos ................................................................................................................... 6
3. Zona de estudio ......................................................................................................... 7
3.1 Segmentos de playa estudiados .......................................................................... 8
Zona 4 - Gandia ........................................................................................................ 9
Zona 5 - Piles .......................................................................................................... 10
Zona 6 - Oliva ......................................................................................................... 10
3.2 Factores que influyen en el régimen sedimentario del litoral ............................... 11
3.2.1 Agentes naturales ...................................................................................... 11
3.2.2 Actuaciones antrópicas ............................................................................. 15
4. Metodología ............................................................................................................ 21
4.1 Realización de las tareas del trabajo ................................................................ 21
4.2 Software empleado .......................................................................................... 21
4.3 Extracción y análisis de la línea de costa ......................................................... 22
4.3.1 Extracción de los puntos de la línea de costa a partir de imágenes Landsat.
…………………………………………………………………………...22
4.3.2 Obtención de las líneas de costa. Paso de puntos a líneas ........................ 25
4.3.3 Influencia del número de líneas registrado en cada año en el análisis ..... 28
4.3.4 Preparación del análisis de las líneas de costa en DSAS. Construcción de la
línea base y de los transectos .................................................................................. 30
4.3.5 Filtrado y eliminación de errores ............................................................. 33
4.3.6 Cálculo de resultados estadísticos mediante DSAS ................................. 34
4.2 Granulometrías del sedimento de las playas......................................................... 37
4.2.1 Muestreo ........................................................................................................ 38
4.2.2 Procesado de las muestras ............................................................................. 38
4.2.3 Estudio de los resultados ............................................................................... 40
4.3 Cálculo de las pendientes de las playas ................................................................ 41
5. Resultados ............................................................................................................... 44
5.1 Evolución de las líneas de costa ...................................................................... 44
5.1.1 Cambios interanuales. Evolución de las líneas de costa a largo plazo ..... 44
5.1.2 Cambios intranuales. Variación de las líneas de costa durante un año. ... 54
5.2 Sedimentología ..................................................................................................... 57
5.3 Cálculo de la pendiente de las playas ................................................................... 62
6. Discusión ................................................................................................................ 64
6.1 Análisis de la evolución de las líneas de costa ..................................................... 64
Zona 4 - Gandia ...................................................................................................... 64
Zona 5 - Piles .......................................................................................................... 65
Zona 6 - Oliva ......................................................................................................... 66
Análisis conjunto .................................................................................................... 67
6.2 Análisis sedimentológico ...................................................................................... 71
6.3 Análisis de la relación entre las desviaciones estándar anuales de las líneas de costa
y las características morfológicas de las playas.......................................................... 72
6.3.1 Relación entre la pendiente de las playas y la variación intranual de la línea
de costa …………………………………………………………………………...73
6.3.2 Relación entre la textura de las playas y la variación intranual de la línea de
costa …………………………………………………………………………...74
7. Conclusiones y recomendaciones ........................................................................... 77
Bibliografía ..................................................................................................................... 79
Anexos ............................................................................................................................ 83
Anexo I – Datos referentes a las zonas 1, 2 y 3 .......................................................... 83
1 - Zona de estudio ................................................................................................. 84
2 – Actuaciones y movimientos de arena ............................................................... 86
3 – Evolución de la línea de costa .......................................................................... 87
4 – Sedimentología ................................................................................................. 94
5 – Pendiente del frente de playa ............................................................................ 95
Anexo II – Información adicional .............................................................................. 96
Anexo III - Parámetros estadísticos DSAS............................................................... 100
Bibliografía de los anexos ........................................................................................ 102
Agradecimientos ........................................................................................................... 103
Currículum Vitae .......................................................................................................... 104
Listado de acrónimos
DGPS - Differential Global Positioning System
DSAS - Digital Shoreline Analysis System
EPR - End Point Rate
GIZC – Gestión Integral de Zonas Costeras
GPS - Global Positioning System
LIDAR - Light Detection And Ranging
LMS - Least Median of Squares
LRR - Linear Regression Rate
MDE – Modelo digital de elevaciones
NSM - Net Shoreline Movement (NSM),
PNOA - Plan Nacional de Ortofotografía Aérea
SCE - Shoreline Change Envelope
SIG o GIS – Sistema de Información Geográfica
1
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
1. Introducción
1.1 Importancia de las playas
Las playas son espacios con una gran importancia a nivel global asociada al turismo,
debido al impacto económico que generan al constituir un foco para el desarrollo de una
actividad turística fundamentada en el atractivo de las playas y la existencia de un clima
suave (Yepes y Medina, 2005). En este sentido, España constituye según la Organización
Mundial del Turismo uno de los destinos turísticos más populares (OMT, 2012),
alcanzando los 57,7 millones de visitantes internacionales, motivados en gran medida por
las ofertas de turismo de masas de “sol y playa” (Yepes y Medina, 2005; Obiol y Pitarch,
2011).
En zonas como el litoral de Valencia resulta evidente la importancia que tiene el
aprovechamiento de las playas por parte de este sector, hasta llegar a ser el producto
económico más importante de la región (Yepes y Medina, 2005), generando el 13% del
PIB (Cámara de comercio de Valencia, 2011).
El turismo debe ser considerado como una actividad económica, asociada a zonas costeras
que gozan de buen clima, que genera presión sobre el medio, donde las playas, espacios
naturales de gran fragilidad, constituyen un recurso limitante (Yepes y Medina, 2005).
La actividad turística también ha dado lugar a modificaciones importantes a nivel
socioeconómico al desplazar a otras actividades primarias tradicionales, y a nivel
territorial al motivar la expansión urbanística que el litoral valenciano ha sufrido desde
los años 60 (Obiol y Pitarch, 2011; Obiol, 2003; Brocal et al., 2001, generando un modelo
territorial y turístico donde el cuidado del entorno constituye una preocupación secundaria
(Obiol y Pitarch, 2011).
En este sentido se pueden destacar actuaciones de crecimiento urbano que han tenido un
efecto directo sobre el frente costero degradándolo y alterando la interacción en el sistema
playa-duna, como por ejemplo la construcción de edificaciones sobre cordones dunares,
inmovilizando o eliminando la arena que en ellos se encontraba, y reduciendo la
disponibilidad de sedimento en las playas, lo que en muchos casos se traduce en patrones
erosivos que llegan a poner en riesgo el mantenimiento de las playas.
2
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
No obstante, las alteraciones antrópicas incluyen una gran diversidad de acciones tales
como la disminución de entrada de sedimentos al sistema costero por actuaciones sobre
la red fluvial, así como la construcción de diferentes tipos de obras duras en recintos
portuarios, desembocaduras de ríos, e incluso con el objetivo de conservar playas, lo cual
se traduce en la mayoría de las ocasiones en un freno al transporte longitudinal de los
sedimentos, agravando el problema erosivo. Todo esto conduce a un escenario en el que
la mayoría de las playas del Golfo de Valencia sufren procesos de erosión o acreción
inducidos por construcciones humanas o una gestión del espacio litoral inadecuadas
(Serra y Medina, 1996). En numerosas playas encontramos situaciones graves de
inestabilidad, con importantes costes ambientales, económicos y sociales.
La fragilidad de estas áreas viene dada no únicamente por la presión que el ser humano
aplica sobre ellas, sino también por la vulnerabilidad que presentan ante un escenario de
ascenso del nivel del mar (IPCC Working Group II, 2014), constituyendo una
preocupación global ante este tipo de amenazas.
Dada la gran importancia y la fragilidad de los espacios costeros, se vuelve indispensable
alcanzar una gestión adecuada, que considere a estas zonas como un sistema
multidimensional, donde se producen interacciones entre sistemas naturales (como la
biota, los sedimentos y el agua), socio-culturales y de gestión (James, 2000). La gestión
integrada de las zonas costeras (GIZC) aparece como una solución viable en la búsqueda
de beneficios económicos, ambientales y sociales, donde la gestión se debería realizar
siguiendo un enfoque ecosistémico, adaptado al manejo de playas (Ariza et al., 2008). La
GIZC pretende alcanzar la sostenibilidad de las zonas costeras mediante un proceso
continuo de gestión responsable, que toma como punto de partida la interacción entre los
diferentes actores principales de éstas y los avances científicos.
Para la implementación de una gestión integrada adecuada en la costa se necesita disponer
de información y datos suficientes, ya que es habitual encontrarse problemas de falta de
los mismos al plantear sistemas de gestión (Osorio et al., 2007). Por tanto, resulta
necesario realizar una adecuada recogida de datos (Ariza et al., 2008; DGC, 2008) físicos,
naturales y socioeconómicos, que permitan llevar a cabo una caracterización de este
entorno.
El análisis de líneas de costa se presenta en este sentido como una alternativa interesante
y capaz de responder al objetivo buscado.
3
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
1.2 Utilidad del análisis de la línea de costa
La morfodinámica costera considera, de forma conjunta y con interacciones entre sí, la
morfología de la costa y la dinámica de los fluidos, las corrientes litorales, el oleaje y el
transporte de sedimentos (Voigt, 1998). La topografía del fondo juega un papel
fundamental en el movimiento del agua, responsable del transporte de sedimentos, que a
su vez modifica la morfología del fondo para, finalmente, condicionar los procesos
hidrodinámicos (Ojeda, E., 2008). La evolución de la zona costera responde por tanto a
la interacción de todos estos elementos, siendo por tanto objeto de diferentes procesos de
gran dinamismo, a diferente escala temporal y espacial.
Estos cambios pueden ser considerados desde diferentes escalas espaciales y temporales
(Ojeda, E., 2008), que van desde cambios en la tendencia sedimentaria general de toda
una célula costera a la formación de pequeños ripples que pueden cambiar en minutos
(Becker et al., 2007), pasando por cambios estacionales en el perfil de la playa a lo largo
de un año (Komar, 1998), o a eventos de tormenta que pueden generar cambios bruscos
en el perfil de playa en periodos de tiempo cortos (Shepard, 1950).
Los esfuerzos por encontrar técnicas que permitan analizar estos cambios pasan por el
desarrollo de nuevos métodos que sean capaces de integrar los grandes volúmenes de
información disponible, haciendo que resulten significativos a nivel geomorfológico, y
útiles para la posterior gestión del litoral. Además, deben permitir un nivel de
automatización lo más elevado posible, pudiendo explotar al máximo los registros
existentes en la actualidad (Pardo et al., 2008).
En la búsqueda de indicadores del estado de la costa que permitan llevar a cabo un
seguimiento de su estado y evolución a lo largo del tiempo, el estudio de la posición de
la línea de costa parece ser una gran opción debido a la posibilidad de detectar cambios
ante variaciones a corto plazo, como las asociadas a tormentas (Ojeda, E., 2008), como
también al permitir analizar cambios de tendencias a nivel general durante periodos de
tiempo mayores, constituyendo un elemento de gran utilidad a llevar a cabo labores de
gestión y planificación de la zona costera (Ojeda, J. et al., 2013). Este estudio considera
los cambios registrados dentro de una escala temporal corta, recogiendo los datos
existentes en un periodo de 30 años, considerando procesos que puedan generar una
tendencia erosiva o acumulativa en los diferentes sectores del litoral.
4
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
1.3 Alternativas existentes
Existen diferentes técnicas disponibles que permiten obtener información referente a la
posición de la línea de costa, para posteriormente estudiar los cambios sufridos por la
misma.
En primer lugar, como instrumento de mayor recorrido histórico, se encuentra la
cartografía. Aquí se engloban tanto los mapas topográficos como las cartas náuticas. En
ambos casos la información disponible con calidad suficiente para ser utilizada para
detectar cambios en la línea de costa es inexistente hasta finales del siglo XIX, lo que
unido a una baja frecuencia de muestreos, separados en el tiempo, hace que únicamente
resulten útiles para detectar tendencias generales muy marcadas (Ojeda, E., 2000).
La fotografía aérea, de menor recorrido histórico, constituye instrumento más utilizado
para detectar cambios y patrones de la línea de costa. El empleo de estos documentos,
existentes en España desde los años 40 (y con mayor frecuencia desde los 70), permite
definir la línea de costa mediante fotointerpretación. La línea de contacto entre la duna
costera (foredune) y la playa alta se puede emplear como un buen indicador de los
cambios a largo plazo. En cualquier caso, al igual que el empleo de cartografía únicamente
permite identificar cambios significativos (Ojeda, E., 2000).
A escala de mayor detalle se encuentra el procedimiento clásico de emplear estaciones
totales, para la realización de perfiles transversales entre los que realizar interpolaciones.
Como sustitutivo de levantamientos topográficos con estación total aparece el DGPS, un
sistema de posicionamiento global diferencial para medir de forma continua, basado en
el empleo de GPS fijo en una base georreferenciada junto a un GPS móvil que se desplaza.
También con un gran nivel de detalle, la videomonitorización cuenta con la ventaja de
permitir el análisis de cambios al capturar diferentes imágenes de la misma zona,
permitiendo el análisis de fenómenos como los temporales. Esta técnica toma fotografías
de la costa desde un punto elevado, para posteriormente reproyectarlas de forma oblicua
y obtener imágenes en planta. Sin embargo, estas técnicas que aportan un mayor nivel de
detalle requieren mucho tiempo y dinero, a la vez que no permiten una comparación
retrospectiva, quedando relegadas al análisis de zonas muy localizadas (Ojeda, E., 2000).
La altimetría láser aerotransportada, conocida como LIDAR, es una técnica de
popularidad creciente, basada en la combinación de un vector aéreo equipado con GPS y
5
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
sistema de navegación inercial junto a un altímetro láser de barrido. Es una técnica que
permite una toma de datos eficiente, rápida y gran precisión, (Pardo et al., 2014) aporta
información de carácter continuo. Su principal limitación va ligada a su corta existencia,
de forma que no tiene posibilidad retrospectiva, por lo que su aplicación únicamente
puede ir ligada al estudio de cambios muy recientes o de cambios en el futuro, contando
asimismo con una baja frecuencia de toma de datos (Pardo et al., 2014).
Finalmente nos encontramos con la teledetección espacial, fuente de los datos empleados
en este trabajo. Esta técnica ha ofrecido un gran volumen de datos de las zonas litorales
desde mediados del siglo pasado, con una periodicidad alta, aunque su baja resolución
espacial limita su empleo para calcular cambios en la línea de costa (Ojeda, E., 2000).
Ninguna de las técnicas anteriormente citadas responde adecuadamente al interés por
realizar un análisis de la evolución de la línea de costa a largo plazo, y que pueda ser
aplicado de forma generalizada a toda la zona costera a un coste bajo.
1.4 Extracción de la línea de costa con precisión subpíxel a partir
de imágenes Landsat TM
A partir del año 2008 aparece la posibilidad de tener acceso libre a la información de los
datos geográficos almacenados por el United States Geological Survey (USGS), lo que
da lugar a que centros de investigación trabajen en el desarrollo de diferentes aplicaciones
para aprovechar esta información.
En este trabajo se pretende seguir el método propuesto por Pardo et al. (2012), que emplea
las imágenes Landsat TM disponibles desde marzo de 1984 a la actualidad, junto con
técnicas de análisis digital de imagen, para alcanzar una resolución de 5 m en la
determinación de la posición de la línea de costa. El método utilizado, que se expone con
mayor detalle en el apartado metodológico, permite llevar a cabo una determinación de
la posición de la línea de costa con una gran precisión, a nivel subpíxel, y un elevado
grado de automatización, permitiendo su aplicación a escala global, contando con
información retrospectiva con la que realizar comparaciones, abriendo la puerta a un
mayor conocimiento de la variación de la línea de costa, y por lo tanto a una gestión más
eficiente de las playas (Pardo et al., 2012).
El empleo de este método sobre líneas de costas micromareales permite por lo tanto llegar
a obtener información referente a la variabilidad registrada a nivel intranual, ya que se
6
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
dispone de un gran volumen de información sobre la línea de costa (Fig 1.1), conociendo
su posición en diferentes momentos del año, lo que permite que pequeños cambios sí
puedan ser registrados. Además consigue cuantificar patrones erosivos o acumulativos a
medio plazo, al contar con información desde 1984 hasta la actualidad Pardo et al. (2012).
Fig 1.1. Evolución de la línea de costa en la playa de Gandia Nord. Las líneas de 1984 a 2014 se muestran en una escala de marrones (más antiguas), verdes y azules (más recientes). Se observa una situación de clara acreción, con las líneas de
costa más recientes (azul oscuro) localizadas mar adentro. Se puede apreciar el gran número de líneas analizado y la complejidad del análisis.
2. Objetivos
En este trabajo se pretende determinar la evolución seguida por la línea de costa en
diferentes segmentos de playas valencianas, aplicando una metodología que permite, a
partir de imágenes satelitales Landsat TM y ETM+, la extracción automatizada de líneas
de costa a nivel sub-píxel.
A partir de las líneas extraídas se analiza la variabilidad de la costa para intentar establecer
una relación con la pendiente y la granulometría de las playas, así como determinar
patrones y tendencias evolutivas. Se evalúa la metodología seguida y su utilidad para
determinar la evolución y las tendencias seguidas por la línea de costa en playas arenosas
tanto a lo largo de un año como a largo plazo, ya que esta metodología constituye una
alternativa más económica y sencilla para conocer las características de las playas.
La hipótesis planteada indica que, en un régimen natural, las playas que sufren el efecto
de la erosión muestran un estrechamiento y un aumento de su pendiente, lo que genera
una mayor variabilidad intranual, a la vez que presentan una granulometría del sedimento
más gruesa. No obstante, esta relación puede verse alterada por las diferentes actuaciones
antrópicas realizadas sobre la costa, tales como realimentaciones de arena y
7
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
construcciones, por cambios naturales del nivel del mar, así como por errores en la
metodología de la obtención de las líneas de costa.
3. Zona de estudio
La zona estudiada se engloba en una unidad morfodinámica localizada en el Mediterráneo
occidental, enclavada entre la desembocadura del río Ebro, al norte, y el cabo de San
Antonio, al sur (Serra y Medina, 1996). El litoral valenciano está constituido por 454 km
de costa, de los cuales 265 son playas, (MOPU, 1970 (citado en Yepes y Medina, 2005))
la mayor parte arenosas, aunque más al sur presenta playas de grava de manera
intermitente. El litoral valenciano, donde se localizan los segmentos de playa estudiados,
se encuentra caracterizado por presentar una gran diversidad de formas acumulativas. La
formación de barras litorales ha dado lugar a las playas, pero también a diferentes zonas
de marjal y lagunas costeras, como la Albufera de Valencia, así como a abanicos
aluviales, llanuras de inundación y cordones dunares (Sanjaume, 1985; MAGRAMA,
2007).
Todas las líneas de costa estudiadas aparecen asociadas a sectores de arena, encontrando
tanto segmentos de playa urbanizados y transformados por la presencia humana, como
tramos en los que todavía aparecen cordones dunares (Fig 3.1). En este sentido, de norte
a sur encontramos, como segmentos dunares de mayor entidad, las zonas de la Devesa
del Saler, con un campo dunar que supera las 850 ha (Sanjaume, 1985), la zona de la
playa de l’Auir entre Xeraco y Gandia, el extremo sur del frente costero de Piles, y la
zona sur del municipio de Oliva. La presencia de estos cordones podría jugar un papel
clave a la hora de disminuir la acción del proceso erosivo frente a lugares donde no existen
o han sido destruidos (Pardo y Sanjaume, 2001).
En cuanto a la playa sumergida, muestra importantes diferencias respecto a las pendientes
y a la presencia de barras sumergidas entre las diferentes zonas. Pardo et al. (2014)
analizaron los perfiles elaborados por la autoridad costera de Valencia (Dirección General
de Costas de Valencia) entre la costa y la isobata de 5m de profundidad, para así poder
obtener las pendientes. Se determinó que las pendientes suaves aparecían habitualmente
ligadas a lugares acumulativos, como por ejemplo en las playas de la Patacona o Gandia,
al norte de los puertos de Valencia y Gandia. Por el contrario, al sur de estos puntos las
playas tienen pendientes más pronunciadas y patrones erosivos.
8
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
A nivel morfosedimentario, los tramos estudiados se localizan dentro de una célula
sedimentaria que engloba todo el sector meridional del golfo de Valencia, en el que se
produce un transporte longitudinal N-S de sedimentos (Brocal et al., 2001). Las playas
localizadas al sur del puerto de Valencia han sufrido durante las últimas décadas una
tendencia erosiva, con una situación deficitaria que afecta a sectores muy amplios del
litoral (Pardo y Sanjaume, 2001), motivada principalmente por la interrupción a cargo del
puerto de la deriva litoral, así como por el descenso de los aportes de nuevos materiales
desde los ríos.
3.1 Segmentos de playa estudiados
La información disponible sobre la localización de la línea de costa aparece dividida en 6
segmentos diferentes de costa, todos ellos de sustrato arenoso. Los tres segmentos de
playa localizados más al norte, en las playas de la Patacona, la Malvarrosa, els Ferros, la
Garrofera, l’Alcatí y la Punta aparecen localizados en el municipio de Valencia, a
excepción del primero de ellos, que en parte se localiza en Alboraia. Los otros tres
segmentos de línea de costa considerados se encuentran en los municipios de Gandia,
Piles y Oliva, todos ellos en la comarca de la Safor.
Tabla 3.1. Zonas de estudio de la línea de costa
Zona Municipio Playa Longitud (m)
1 Alboraia - Valencia La Patacona y la Malvarrosa 2840
2 Valencia Els Ferros y la Garrofera 1710
3 Valencia L’Alcatí y la Punta 1800
4 Gandia l’Auir, Gandia Nord 3630
5 Piles Piles 1710
6 Oliva Terranova, Pau Pi, l'Aigua Blanca 2990
En este apartado se describen las zonas o segmentos costeros Nº 4, 5 y 6, objeto de este
trabajo, y cuyos resultados aparecen plasmados en la presente memoria, en el apartado
correspondiente. No obstante es importante recordar que se ha trabajado de forma
conjunta con las zonas Nº 1, 2 y 3, localizadas más al norte (Tabla 3.1), cuyos resultados,
obtenidos por Jesús Soriano en su trabajo de fin de máster, se introducirán en el apartado
de discusión, para así dotar al análisis de una mayor riqueza de datos, y mayores
posibilidades de entender los fenómenos que tienen lugar.
9
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 3.1. Localización de las zonas de línea de costa estudiadas (azul). Encontramos la zona 4 (izquierda) en las playas de l’Auir i Gandia Nord, y las zonas 5 y 6 (derecha) en el municipio de Piles y Oliva. Los frentes costeros se dividen entre los
que conservan cordones dunares (verde) y los antropizados (rosa).
Zona 4 - Gandia
Esta zona de estudio, de 3630 metros de longitud, incluye la parte sur de la playa de l’Auir
y la parte norte de la playa de Gandia Nord (Fig 3.2). L’Auir es la playa localizada en el
extremo norte del municipio de Gandia, cuenta con un cordón dunar que se extiende desde
el río Vaca, al norte, y constituye un espacio litoral sin urbanizar. Por el contrario, la parte
norte de este segmento incluye la playa de Gandia Nord, urbanizada, con un frente costero
totalmente artificial, y en cuyo extremo sur aparece el puerto de la localidad. Ambas
playas presentan anchuras muy considerables, llegando a superarse los 130 metros junto
a la instalación portuaria.
Fig 3.2. Playa de l’Auir y su cordón dunar (Izq.) Autor: Carlos Cabezas. Playa de Gandia Nord (drcha..) Fuente: http://www.magrama.gob.es/es/costas/servicios/guia-playas/
10
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Zona 5 - Piles
Este sector presenta las líneas de costa de la playa de Piles, que ocupa todo el frente
costero del municipio, urbanizado en su mayor parte, excepto en el extremo sur, que
conserva una pequeña formación dunar (Fig 3.3). Se trata de una playa relativamente
estrecha, en la que en ocasiones el paseo marítimo es alcanzado por las olas durante
temporales.
Fig 3.3. Imágenes de la playa de Piles: su tramo con paseo marítimo (izq.) y la zona sur que presenta cordón dunar (drcha.). Autor: Carlos Cabezas.
Zona 6 - Oliva
El último segmento de línea de costa estudiado se sitúa en el municipio de Oliva, al sur
del puerto, en las playas de Terranova, Pau Pi y l'Aigua Blanca (Fig 3.4). La parte norte
de este segmento, en la playa de Terranova, presenta un frente costero urbanizado,
mientras que más hacia el sur aparecen tramos sin urbanizar, o bien en los que la
urbanización se ha llevado a cabo a una mayor distancia del mar, conservándose el cordón
dunar en muchos casos.
Fig. 3.4. Playas de Pau Pi y l’Aigua Blanca, ambas en la zona 6. Autor: Carlos Cabezas.
11
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
3.2 Factores que influyen en el régimen sedimentario del litoral
La presencia de sedimentos en el sistema litoral juega un papel fundamental a la hora de
entender la evolución de las playas. Así pues, los fenómenos erosivos van en muchos
casos ligados a la falta de sedimento en el sistema, desapareciendo la arena presente en
las playas de forma lenta pero progresiva sin que sea reemplazada.
Esto puede deberse tanto a la falta de entrada de nuevos sedimentos al sistema litoral,
como a cambios en su distribución lo que, a su vez, puede tener un origen natural, o bien
deberse a actuaciones llevadas a cabo por el ser humano sobre la costa.
3.2.1 Agentes naturales
Entre los agentes de origen natural que influyen en la introducción, distribución y
retención de sedimentos en el sistema litoral encontramos el oleaje, las mareas, las
corrientes marinas, el nivel del mar, el viento, el clima, y la geomorfología de la zona.
-Oleaje
Las olas son oscilaciones regulares de la superficie del mar, motivadas generalmente por
la acción del viento, y que constituyen el principal agente responsable del modelado de
la costa, erosionando o acumulando los sedimentos presentes a lo largo del litoral. La
costa valenciana presenta un oleaje de periodos cortos y poca altura, en el que las brisas
marinas juegan un papel importante regulando los periodos y la altura de las olas (Pardo,
1991). La altura de ola para esta zona de costa se sitúa en torno a los 80 cm, con periodos
de retorno cercanos a los 6 s (www.puertos.es). Las situaciones con mayores olas suelen
ir asociadas a periodos de tormenta, con una situación de bajas presiones motivada por
westerlies afectando a las latitudes medias (Pardo et al., 2014). En esta zona la presencia
de tormentas es habitual durante el otoño, apareciendo asociadas a la presencia de áreas
de baja presión en el golfo de Lyon o en el mar de Alborán (Pardo, 1991).
Se debe destacar en este apartado la incidencia en noviembre de 2001 de un fuerte
temporal sobre el Mediterráneo occidental (Ramis et al., 2002; Obiol 2003; Julià, 2003;
www.puertos.es) que afectó a la zona estudiada, siendo uno de los más severos de los que
se tiene constancia, y que generó importantes procesos erosivos en las playas al movilizar
gran cantidad de sedimento. Asimismo, diferentes tormentas sucedieron de forma
consecutiva durante los siguientes meses hasta la primavera del 2002 (Pardo et al., 2014).
12
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
En la zona del golfo de Valencia el oleaje tiene un origen mayoritariamente del este,
aunque para la zona de la Safor (Fig. 3.5) sigue una tendencia muy marcada, actuando
principalmente del NE (www.puertos.es), lo cual está relacionado con los vientos
predominantes en la zona, de gregal (NE). El oleaje no ha sufrido variaciones
significativas entre los años 1984 y 2014 (www.puertos.es) por lo que su influencia en
los resultados del análisis evolutivo de la línea de costa no es significativa.
Fig. 3.5. Altura significativa del oleaje medido en medias y máximos mensuales (izq.) y rosa de oleaje (drcha.), ambos medidos para el periodo 1984-2014 en el punto SIMAR 2084108. Fuente: www.puertos.es.
-Corrientes
El patrón de corrientes viene determinado por la circulación existente a gran escala en el
mar Mediterráneo, constante debido a la distribución de masas de agua en la región
occidental, en equilibrio dinámico, lo que se traduce en una circulación termohalina
ciclónica con dirección N-S (Font et al., 1986). Vientos locales y descargas continentales
generan alteraciones que dan lugar a patrones de mayor variabilidad en la zona de la
plataforma continental.
El golfo de Valencia presenta una deriva bien definida en dirección norte-sur entre
Benicarló y Oliva (Pardo y Sanjaume, 2001), incluyendo a todas las playas que son objeto
de este estudio. La diferente orientación de la costa, con una inclinación respecto a las
olas variable, genera cambios en la capacidad de transporte a lo largo del litoral (Brocal
et al., 2001). Al norte del puerto de Valencia el transporte era alto (Pardo et al., 2014),
mientras que en las playas situadas justo al sur del puerto de Valencia, en la zona del
Saler, el transporte longitudinal era muy intenso en el pasado, llegando a alcanzar valores
superiores a los 100.000 m3/año (MAGRAMA, 2007). No obstante, el cambio en la
configuración de la costa, con la realineación de la misma, ha generado que el transporte
disminuya hasta valores cercanos a los 50-70.000 m3/año (MAGRAMA, 2007).
13
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Siguiendo hacia el sur encontramos el Cap de Cullera, que no constituye un obstáculo al
paso de sedimentos, continuando el mismo transporte hacia el sur. A partir de la
desembocadura del río Xúquer el transporte es menor (Pardo et al., 2014) y se encuentra
un patrón erosivo asociado a tasas de transporte estimadas en 60-80.000 m3/año, con un
descenso gradual desde la zona de Xeraco hasta el puerto de Gandia, donde se registran
de 20 a 30.000 m3/año. El puerto actúa como obstáculo, generando acumulación al norte
del mismo, mientras que en las playas del sur existen patrones erosivos importantes. Más
hacia el sur el transporte experimenta un crecimiento que va desde los 35.000 m3/año en
la playa de Miramar hasta unos 40-50.000 m3/año en el entorno del Puerto de Oliva, al
sur del cual vuelve a registrarse un nuevo descenso (MAGRAMA, 2007) hasta ser
prácticamente nulo (Pardo et al., 2014)
-Mareas y nivel del mar
El nivel del mar puede influir al analizar la evolución de la línea de costa de diferentes
formas: bien a través de variaciones temporales, a consecuencia de la marea astronómica
o debido a cambios atmosféricos, o bien asociado a variaciones globales del nivel marino.
La oscilación media de la marea en el Mediterráneo es escasa, habitualmente inferior a
los 20 cm, por lo que juega un papel menor en el modelado costero de este espacio,
pudiendo no obstante alcanzarse valores bastante mayores de manera puntual.
La variación estacional del nivel del mar es otro factor a tener en cuenta. La costa
valenciana posee una variación anual máxima media de 0,37m (Pardo et al., 2011),
ocasionando que se puedan registrar movimientos en la posición de la orilla asociados al
nivel del mar de hasta 5,7m (Pardo et al., 2012).
Por otra parte, el ascenso del nivel del mar a nivel global, también registrado en el
Mediterráneo, parece ser un factor a tener cada vez más en cuenta (IPCC Working Group
II, 2014), y podría resultar esencial a la hora de considerar la evolución de la línea de
costa y el avance o retroceso de nuestras playas. Respecto a este ascenso existen diferentes
opiniones a la hora de cuantificar la variación sufrida a lo largo del siglo XX. Tal y como
apunta Church (2001), los valores de ascenso anual estimados por los diferentes
organismos y autores son poco clarificadores, pues cada modelo empleado pronostica un
patrón de elevación del nivel del mar diferente. Según el IPCC, durante el siglo XX el
nivel medio del mar a nivel global ascendió a un ritmo cercano a los 1,7 mm al año, siendo
mayor este ascenso a partir de la década de los 90, y alcanzando un ritmo de unos 3
mm/año, bastante mayor al promedio observado en la primera parte del siglo. Esta
14
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
variación en la tendencia concuerda con los datos mostrados por el Instituto Español de
Oceanografía (Vargas et al., 2010), que cifra el ascenso para el Mediterráneo español en
1,5 mm al año, hasta un total de 15 a 20 cm a lo largo del siglo, siendo este incremento
mayor a partir de los años 90.
-Viento
La acción del viento condiciona el modelaje de la costa de forma directa al realizarse el
transporte eólico de partículas de sedimento mientras que, indirectamente, lo hace al ser
el generador principal del oleaje. En playas de arena la acción del viento y de las olas
puede generar cambios importantes con bastante rapidez, por lo que se vuelve necesario
considerar su intensidad, frecuencia y dirección. La naturaleza micromareal del
Mediterráneo no hace sino acrecentar la importancia del viento.
Fig. 3.6. Velocidad del viento medido en medias y máximos mensuales (izq.) y rosa de los vientos (drcha.), ambos medidos para el periodo 1984-2014 en el punto SIMAR 2084108. Fuente: www.puertos.es
Con este objetivo se muestra de forma gráfica la frecuencia y la rosa de vientos (Fig.3.6)
registrada en esta zona para el periodo 1984-2014. Esta información procede de un punto
SIMAR, localizado frente a la costa de la Safor, en las coordenadas 39ºN, 0ºW
(www.puertos.es).
Se puede decir respeto a la dirección de los vientos que el del NE o gregal, y el poniente
(W) son los más habituales y que alcanzan mayores intensidades. No obstante existe una
alternancia estacional entre los vientos de levante y poniente, siendo estos últimos más
habituales durante el invierno. Asimismo se puede afirmar que no han existido cambios
de tendencia significativos en cuanto a la velocidad del viento se refiere a lo largo del
periodo estudiado.
15
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
3.2.2 Actuaciones antrópicas
Las actuaciones llevadas a cabo por el ser humano sobre los espacios costeros pueden
tener diferentes efectos a nivel temporal, espacial y de magnitud dependiendo de su
naturaleza.
Actuaciones sobre el medio litoral
Así pues, por una parte se pueden encontrar cambios rápidos y delimitados en el espacio,
habitualmente relacionados con modificaciones puntuales del transporte de sedimentos a
lo largo del sistema costero debido al establecimiento de obstáculos artificiales al
transporte longitudinal, o bien con actuaciones directas de redistribución, retirada y
retención del sedimento.
En la costa valenciana la mayoría de los desequilibrios sedimentarios tienen su origen en
la construcción de obras artificiales (Pardo y Sanjaume, 2001). La alteración del
transporte longitudinal tiene lugar cuando obras antrópicas, como puertos o diques, se
localizan en la zona litoral, de forma más o menos perpendicular a la costa, disminuyendo
la deriva natural del sedimento y retener parte del mismo. En la costa Valenciana, donde
tiene lugar un transporte norte-sur, este tipo de obstáculos generan la sobreacumulación
de sedimento en su parte norte, con un progresivo avance de la línea de costa, que genera
una rotura de las olas cada vez más lejos, llegando a la orilla con menos energía,
retroalimentando el proceso, mientras que dejan desprovistas de sedimento a las playas
del sur, en las que se producen procesos erosivos y de retroceso de la línea de costa. Este
retranqueo se debe a que las arenas retenidas en el obstáculo ya no pueden reponer a las
desplazadas por la erosión natural aguas abajo, generando una situación deficitaria
crónica y aumentando progresivamente la pendiente de la playa, en un proceso inverso al
mencionado anteriormente.
En nuestra zona de estudio encontramos como principales actuaciones antrópicas con
consecuencias sobre el transporte sedimentario las sucesivas ampliaciones del puerto de
Valencia (Fig 3.7), el cual actúa como un elemento modificador muy importante de la
deriva litoral sobre la zona estudiada, al obstaculizar el paso de gran cantidad de
sedimentos, con consecuencias para las playas del entorno (Serra, 1998), así como el
puerto de Gandia y el puerto deportivo de Oliva. También aparecen diferentes espigones
asociados a las golas y desembocaduras de ríos, como en el caso del Xúquer y el Túria
16
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig.3.7. Ampliaciones del puerto de Valencia, años 1980, 2006 y 2015. Fuente: http://googlemaps.es
No obstante, la acción humana también juega un papel clave al alterar directamente la
distribución natural de los sedimentos mediante la realización de dragados y vertidos de
arena, de elevado coste (Obiol, 2003), con el objetivo de aportar arena a playas en
situación de retroceso, permitiendo que el uso de la playa se pueda seguir llevando a cabo.
La tabla 3.2, recopila las actuaciones llevadas a cabo sobre las playas estudiadas (Fig.
3.8, 3.9), siendo de gran importancia pues pueden jugar un papel clave en la evolución
de la línea de costa, desvirtuando los patrones que las playas estaban experimentando
antes de las actuaciones.
Tal y como se puede observar, las playas de Gandia y Piles han sufrido diferentes
movimientos de arena en las últimas décadas, actuando en muchos casos la playa de Piles
como receptora de arena extraída de las playas de Gandia Nord y l’Auir.
17
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 3.2. Actuaciones de dragado y vertido realizadas en la zona de estudio Fuente:
elaboración a partir de datos de la Dirección General de Costas (DGC).
Descripción de la actuación
Lugar de préstamo o dragado
Playa de vertido
Volumen (m3)
X Y Fecha inicio
Fecha finalización
Vertido de arena (Regenerar playas)
Fondo submarino Sur de la
Safor 60000 - - 1993 1993
Movimiento de arena Aparcamientos Subterráneos Gandia
Gandia Nord
49627 745468 4321836 2001 2001
Movimiento de arena (Emergencia)
Norte del dique de levante Gandia
Piles 6000 750078 4315480 2003 2003
Reacomodación de arenas (Mantenimiento)
Playa de Piles Piles 3200 750094 4315464 2006/04 2007/11
Reacomodación de arenas (Mantenimiento)
Paseo Marítimo de Piles Piles 7300 750110 4315441 2006/04 2007/11
Aportación de arenas (Mantenimiento)
Excavación fuera DPM-T Piles 14370 750157 4315388 2006/04 2007/11
Aportación de Arena (Mantenimiento)
Excavación fuera DPM-T Piles 5000 750170 4315377 2006/04 2007/11
Reacomodación de arenas en la playa
Playa de Gandia Nord Gandia Nord
23366 744747 4323149 2007/07 2007
Reacomodación de arenas (Emergencia)
Playa de Gandia Piles 46542 750128 4315418 2008/01 2008/06
Reacomodación de arenas en la playa
Playa de Gandia Bellreguard 50419 744916 4322876 2009/02 2009/04
Reacomodación de arenas en la playa
Playa de Gandia Miramar 49669 744707 4323193 2009/02 2009/04
Reacomodación de arenas
Playa de Gandia Piles 50789 750143 4315401 2009/03 2009/05
Dragado del Puerto de Gandia
Bocana Puerto de Gandia
Piles 38000 750064 4315495 2010/05 2010/06
Aportación de Arena (Mantenimiento)
Excavación fuera DPM-T Piles 500 750197 4315346 2011/06 2014/06
Aportación de Arena (Mantenimiento)
Excavación fuera DPM-T Piles 3440 750185 4315360 2013/02 2013/02
18
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 3.8. Actuaciones de dragado y vertido realizadas frente al municipio de Gandia. Fuente: elaboración a partir de datos de la Dirección General de Costas <www.magrama.gob.es/es/costas/servicios/guia-playas>
19
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 3.9. Actuaciones de dragado y vertido realizadas frente al municipio de Piles. Fuente: elaboración a partir de datos de la Dirección General de Costas <www.magrama.gob.es/es/costas/servicios/guia-playas>
-Retención del sedimento-eliminación cordón dunar
El ser humano también puede afectar a la distribución del sedimento costero por el mero
hecho de inmovilizar las arenas existentes, construyendo sobre las playas o destruyendo
los cordones dunares, y actuando sobre el correcto funcionamiento del sistema playa-
duna.
Un ejemplo de ello lo encontramos al sur del puerto de Valencia, en la zona del Saler, en
la que durante los años 70 se planteó e inició un gran proyecto urbanizador, en el que se
removieron las dunas delanteras, localizando un paseo marítimo en la parte norte. Estas
actuaciones contribuyeron a que las playas de la zona entrasen en un proceso recesivo
(Fernández, 2002). Algunos años después, en un intento de revertir la situación, se
destruyó el paseo marítimo y se fue procediendo a recuperar el cordón dunar, extrayendo
la arena original que había sido depositada en las restingas, así como importando arena
20
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
procedente de playa de la Malvarrosa y del dragado de fondos marinos. Más al sur, y ya
directamente en la zona estudiada se tiene constancia de actuaciones de regeneración
artificial del cordón dunar en Oliva, en la platja de l’Aigua Blanca (Obiol, 2003).
No obstante, pese a que el caso de las dunas de l’Albufera es el más llamativo, en la zona
estudiada encontramos otros ejemplos, como el frente costero del municipio de Piles,
donde se puede observar cómo se optó por localizar el paseo marítimo sobre las dunas
que existían en la zona (Pardo y Sanjaume, 2001), lo que genera la inmovilización de una
gran cantidad de sedimento, unido a un estrechamiento de la playa, lo que plantea
consecuencias ligadas a un aumento de la erosión, ya que se impide la libre adecuación
del perfil playero a las condiciones energéticas del mar, generando efectos semejantes a
las escolleras aunque más lentos, provocando un retroceso de la línea de costa a través de
la migración de materiales mar adentro, algo que ocurre en las playas de Piles y
Bellreguard, justo al norte de la misma (Pardo y Sanjaume, 2001).
Alteración del régimen de entrada de sedimento al sistema
El ser humano también puede actuar sobre el régimen de entrada de sedimento al sistema
litoral, ligado directamente a actuaciones sobre el sistema fluvial, con trampas y
obstáculos que impiden el paso del sedimento, así como a causa de la alteración de la
cuenca de drenaje. De este modo, las actuaciones antrópicas pueden derivar en cambios
de tendencia a largo plazo, con efectos a mayor escala que pueden afectar a la totalidad
de la célula sedimentaria costera.
En las costas valencianas, el principal aporte de sedimento tiene su origen en el sistema
fluvial (Sanjaume, 1985). Los ríos Túria y Xúquer constituyen las principales fuentes de
sedimento a este sector de costa (Brocal et al., 2001) y la regulación que se ha llevado a
cabo de forma generalizada en el sistema fluvial a lo largo del siglo XX mediante presas
y pantanos ha alterado de una forma importante la entrada de sedimentos al litoral (Pardo
y Sanjaume, 2001), al atraparlos e impedir su desplazamiento, a la vez que al regular el
caudal de los ríos y disminuir su capacidad de transporte, destacando en este sentido la
regularización del Xúquer en su cuenca baja mediante diferentes obras transversales.
21
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
4. Metodología
El presente apartado se ha dividido en 3 grandes secciones: análisis de las líneas de costa,
granulometrías del sedimento y cálculo de las pendientes de las playas.
4.1 Realización de las tareas del trabajo
El presente trabajo, que de forma global trabaja con 6 zonas del litoral, se ha realizado de
forma paralela al de Jesús Soriano González, quien en su trabajo expone los resultados
obtenidos para los 3 sectores de costa localizados más al norte (identificados como Zonas
1, 2 y 3), en playas de los municipios de Valencia y Alboraya, mientras que este trabajo
presenta los resultados asociados a las líneas de costa localizadas más al sur, en la comarca
de la Safor, que serán identificados como Zona 4 “Gandia”, Zona 5 “Piles” y Zona 6
“Oliva”. Para el desarrollo de ambos proyectos se han realizado algunas tareas de forma
conjunta, pasando a realizar el análisis de los datos de forma individual.
De este modo, este trabajo muestra únicamente los resultados de las zonas 4, 5 y 6, pese
a que en el apartado de discusión se consideran las 6 zonas y se emplean la totalidad de
los resultados obtenidos (los resultados de las zonas 1, 2 y 3 aparecen recopilados en los
anexos), ya que esto ayuda a aumentar la calidad de ambos trabajos, enriquecer el estudio
en su conjunto y ayudar a comprobar la hipótesis formulada inicialmente.
Se ha optado por realizar los proyectos de este modo ya que se contaba con una gran
cantidad de datos iniciales, que unidos al importante trabajo de procesado y tareas
informáticas requeridas por los mismos, hacían más idónea una división de tareas, para
así reducir el tiempo destinado al procesado de los datos.
4.2 Software empleado
Para poder desarrollar el presente trabajo ha sido necesario emplear una serie de
programas al llevar a cabo el tratamiento de los datos:
Programas GIS. Diferentes sistemas de información gráfica han sido utilizados para poder
procesar, analizar, almacenar y gestionar los datos. También han sido empleados para
elaborar mapas que permitan visualizar de forma gráfica los resultados. Se ha trabajado
con dos programas diferentes, optando por emplear uno u otro dependiendo de la tarea
concreta que se pretendía realizar.
22
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
ArcMap. Este programa, desarrollado por la empresa ESRI, ha sido utilizado en su
versión 10.1. En él se ha añadido la extensión DSAS que se describe posteriormente
QGIS. La versión Desktop Wien de este programa, de software libre, ha sido empleada
para algunas de las tareas o pasos llevados a cabo, al considerarlo más liviano o presentar
algunas funciones no disponibles en ArcMap.
DSAS. La versión 4.3 del Digital Shoreline Analysis System, desarrollada por el USGS,
es un software libre que se ejecuta como una extensión de ArcMap, permitiendo el cálculo
de la velocidad de cambio de la línea de costa a partir de diferentes estadísticos. (Thieler
et al., 2009).
GRADISTAT. Este programa, desarrollado por Blott y Pye (2001), permite agilizar el
análisis de las estadísticas granulométricas del sedimento de las playas, lo que resulta útil
de forma posterior a la realización de granulometrías mediante tamizado. El programa
permite calcular parámetros como la media, clasificación y la asimetría.
4.3 Extracción y análisis de la línea de costa
El análisis planteado en este trabajo pasa por la obtención de unos datos de partida, que
suministran una localización puntual de la línea de costa en diferentes fechas. A partir de
estos datos se obtienen las diferentes líneas de costa, que permiten finalmente llevar a
cabo el cálculo de diferentes parámetros estadísticos con la herramienta DSAS. Los
resultados obtenidos posibilitan comprender la variación sufrida por las líneas de costa
en las franjas de litoral estudiadas, tanto en análisis intranuales que permitan observar la
variación a lo largo de un año, como en un análisis interanual, que dé a conocer la
evolución sufrida en el sector a medio y largo plazo. El periodo de estudio considerado
es de 30 años, entre 1984 y 2014.
4.3.1 Extracción de los puntos de la línea de costa a partir de imágenes
Landsat.
El estudio de los cambios y la evolución en la línea de costa parte de los registros
puntuales de la localización de esta línea en diferentes momentos del tiempo. En el
presente trabajo, los datos puntuales han sido obtenidos siguiendo la metodología
propuesta por Pardo et al. (2012) para la extracción automática de las líneas a nivel
subpíxel mediante la aplicación de un algoritmo sobre imágenes multiespectrales de
satélites Landsat, que permite una alta precisión, cercana a los 5 metros.
23
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
La metodología propuesta toma las imágenes de la base de datos de USGS
(http://earthexplorer.usgs.gov/) (NASA, 2006). Se incluyen imágenes tomadas entre los
años 1984 y 2014, adquiridas a través de los sensores TM en el caso del Landsat 5, y
ETM+ en el caso del Landsat 7.
El proceso, partiendo de las imágenes satelitales iniciales, atraviesa diferentes etapas
hasta la obtención de la línea de costa final. Para conseguirlo, se vuelve necesario realizar
una extracción grosera de la línea de costa en cada una de las imágenes, para a
continuación realizar correcciones a nivel radiométrico y geométrico, que permiten
disminuir el error del proceso.
Georreferenciación a nivel subpíxel
En primer lugar se parte de las imágenes suministradas por la NASA, las cuales ya
cuentan con una buena georreferenciación (mejor a 0,4 píxeles según la NASA, (2006)),
pero dado los objetivos perseguidos se vuelve necesario alcanzar el nivel subpíxel en estas
imágenes. Esto se consigue mediante un algoritmo basado en una transformación discreta
de Fourier, y la aplicación de una correlación cruzada fundamentada en la correlación de
Pearson (Guizar-Sicarios et al., 2008), para determinar el desplazamiento entre las
imágenes que deben presentar la misma resolución. En el proceso se usa, como referencia,
un mosaico de ortofotos del PNOA del espectro visible e infrarrojo con resolución
0.5m/píxel que son degradadas hasta 30m/píxel para así presentar la misma resolución
que las imágenes obtenidas del Landsat. La correlación permite conocer el
desplazamiento existente en el eje x e y, y permite mejorar sustancialmente la
georreferenciación de las coordenadas de las líneas que se extraerán, haciendo que el error
máximo derivado de la georreferenciación pase a estar entre los 2 y 3 metros, cuando
originalmente era de 15 metros.
Extracción automática de la línea de costa a nivel subpíxel
El método seguido se basa en la identificación de la interfaz tierra-agua a partir de las
diferencias espectrales existentes entre el agua y la tierra en el espectro infrarrojo de las
imágenes Landsat, en este caso sobre la banda 5 del Landsat 5 y 7, y la banda 6 del Landsat
8.
Se inicia el proceso con una aproximación desde la escala de píxel mediante una
binarización de la imagen de partida, llevando a cabo una diferenciación entre la
24
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
superficie de tierra y mar mediante una umbralización automática del histograma en la
banda del infrarrojo.
Sobre esta primera línea se pasa a realizar la detección de la posición a nivel subpíxel,
con un remuestreo para trabajar sobre píxeles de tamaño menor. Para lograrlo, en primer
lugar se transforma la interfaz tierra-agua definida anteriormente de forma aproximada
mediante morfología matemática, con la modelización un polinomio de 5º grado, que se
resuelve partiendo de los valores digitales de un vecindario de 7 x 7 píxeles en bandas
infrarrojas. En las bandas del infrarrojo cercano y medio (bandas 4,5 y 7 del Landsat 5 y
7) el agua presenta un valor mínimo y el continente un valor mucho más alto. La
localización de la costa se establece en el punto (definido sobre la superficie matemática)
en que se produce el máximo gradiente entre el continente y el mar. Como cada píxel
aparece incluido en varios vecindarios, lo que genera varias soluciones posibles, siendo
el promedio de ellas la línea definitiva, localizada sobre los puntos en que con mayor
probabilidad se encuentra la interfaz agua-tierra. El análisis se repite de forma sucesiva
a lo largo de la línea definida inicialmente.
Sobre el resultado obtenido resulta necesario aplicar una corrección radiométrica pues en
las imágenes de partida las diferentes superficies llevan asociadas diferentes respuestas
espectrales, que afectan negativamente al modelizado, pudiendo generar desplazamientos
de la línea obtenida a nivel subpíxel. Se actúa determinando mediante una modelización
la relación existente entre el nivel digital del vecindario y la adecuación del ajuste de la
función matemática, y por tanto, el desplazamiento sufrido. De un ajuste por función
polinomial de los valores de tierra se deduce la corrección radiométrica a aplicar sobre
todos los puntos obtenidos
En cuanto al error obtenido, Pardo et al. (2012) plantea una evaluación de los resultados,
mediante la comparación con líneas de costa fijas a lo largo del periodo de estudio. Con
ese propósito se empleaban zonas de escollera, sobre las que se comparaban las líneas de
costa obtenidas a partir de imágenes Landsat con una línea de referencia trazada a partir
de una ortoimagen. La evaluación de errores consistía en el cálculo de las distancias o
desplazamientos respecto a la línea de referencia.
Aplicando la metodología explicada, se obtienen para los segmentos de costa objeto de
estudio un total de 281.618 puntos, de los cuales 151.620 aparecen en la zona 4, en la
zona 5 encontramos 35.689, y en la zona 6 se localizan 94.309 (Tabla 4.1)
25
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 4.1. Longitud de los sectores y número de puntos registrados
Zona Longitud (m)
Número de puntos
4 3630 151620
5 1710 35689
6 2990 94309
Los puntos presentan información de su localización y fecha, localizados cada 7.5 m a lo
largo de la costa, y que constituirán las líneas de costa que se estudian en el presente
trabajo. El flujo de trabajo planteado por Pardo et al. (2012) aparece esquematizado en
la Fig 4.1.
Fig. 4.1. Metodología para la extracción de líneas de costa a partir de imágenes satelitales Landsat (Pardo et al., 2012).
4.3.2 Obtención de las líneas de costa. Paso de puntos a líneas
Mediante el procedimiento descrito anteriormente se obtienen los puntos de localización
de la línea de costa a lo largo de diferentes fechas y para diferentes segmentos del litoral
valenciano. Para llevar a cabo el procesado de los datos hasta llegar a obtener líneas de
costa, se emplean tanto el software ArcMap 10.1 (ESRI), como QGIS Desktop 2.8.1,
utilizados de forma alternativa dependiendo de las tareas y operaciones que se van
desarrollando en cada momento.
El conjunto de puntos del que partimos con la información contenida anteriormente
aparece en un único fichero shape (.shp) operable desde sistemas de información
geográfica. No obstante, para poder llevar a cabo análisis referentes a la evolución de la
26
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
línea de costa se debe contar con un shape de líneas, sobre las cuales sí es posible
identificar cambios, avances y retrocesos entre las diferentes fechas en que han sido
tomadas las imágenes Landsat, así como operar con la herramienta DSAS. Para poder
alcanzar nuestro objetivo en primer lugar, tras reproyectar los puntos al sistema de
referencia ETRS89 y a la proyección UTM huso 30 (N), optamos por separar en
diferentes ficheros shape los puntos de cada una de las 6 zonas diferentes, lo que permitirá
agilizar las tareas en el programa SIG.
En cada una de estas 6 zonas pasamos a seleccionar los puntos referidos a una misma
fecha y, al mismo tiempo, a transformarlos en líneas. Para llevar a cabo esta tarea nos
ayudamos del programa QGIS, que posibilita agilizar el proceso al llevarlo a cabo en un
solo paso y de manera automática, sin tener que realizar una extracción manual de los
puntos asociados a cada una de las fechas. De esta manera, empleando la herramienta
“Points2One” conseguimos obtener un fichero shape de polilíneas, en el que se generan
líneas no cerradas formadas cada una de ellas a partir de los puntos con un campo común,
en este caso el campo de la fecha (Fig. 4.2)
Fig. 4.2. Extensión Points2One de QGIS (izq.). Líneas obtenidas de la unión de los puntos en la Zona 5 (Drcha.).
Posteriormente buscamos que cada una de las líneas pase a estar en un shape diferente,
algo que QGIS permite hacer de forma prácticamente automática a través de la
herramienta “Split vectorial layer”. De este modo obtenemos, para los 6 tramos costeros
diferentes, una serie de líneas de costa, cada una de ellas perteneciente a una fecha
diferente, entre los años 1984 y 2014.
27
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.3. Herramienta “Split vectorial layer” de QGIS.
En este punto debemos destacar que el número de líneas de costa de que disponemos en
cada uno de los tramos, o lo que es lo mismo, el número de fechas para las que tenemos
información de la línea de costa, difiere de unos tramos a otros (Tabla 4.2). De este modo,
la distribución de líneas en cada uno de los años no es homogénea a lo largo de los años
del periodo considerado, existiendo años en que únicamente se han podido obtener 1 ó 2
líneas de costa, o incluso ninguna. Ello se debe a que el segmento 4 (correspondiente al
norte de Gandia) toma imágenes de cuatro escenas de los dos satélites funcionales en el
periodo estudiado (Landsat 5 y 7), dos de ellos centrados directamente en la costa
valenciana y otros dos centrados en las de Baleares, pero que también alcanzan este sector.
Más al sur, sin embargo, solo se dispone de dos escenas, lo que explica el menor número
de líneas. Por otra parte, es importante señalar que en las imágenes Landsat 7, debido al
fallo que el 31 de Mayo de 2003 sufrió el corrector de líneas de escaneado (SLC), las
imágenes llegan con espacios vacíos. Las áreas con registro de imagen sí pueden ser
analizadas (Almonacid et al., 2013) pero obviamente de las zonas sin información no se
puede extraer la línea de costa. Se trata, en todo caso, de segmentos cortos sin
información.
28
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 4.2. Líneas disponibles en cada uno de los segmentos de costa analizados
Segmento 4 Segmento 5 Segmento 6
Año Nº líneas Nº líneas Nº líneas
2014 24 9 9
2013 17 6 6
2012 6
2011 22 13 13
2010 21 11 11
2009 27 12 12
2008 10 4 5
2007 31 12 12
2006 17 6 6
2005 10 2 2
2004 10 7 7
2003 13 6 6
2002 13 7 7
2001 12 8 8
2000 12 6 6
1999 3 1 1
1990 2 1 1
1987 12 7 7
1986 5
1985 3 3 2
1984 8 4 4
Total 278 125 125
4.3.3 Influencia del número de líneas registrado en cada año en el análisis
Tal y como se ha señalado anteriormente, el número de líneas de costa disponible difiere
enormemente de unos años a otros. Debido a ello, se analiza la posible relación existente
entre el número de líneas que se analiza en cada año y su repercusión sobre el resto de
resultados, ya que una alta variabilidad anual podría ir asociada a un número de líneas
disponibles muy bajo en ese año, mientras que en años con mayor disponibilidad de datos
esto podría traducirse en variabilidades menores.
29
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Para comprobarlo se cruzan los datos del número de líneas con las desviaciones estándar
obtenidas para cada uno de los años (Fig 4.4), observando que no existe una relación
marcada entre ambos factores. La desviación estándar no varía al trabajar con más líneas,
posiblemente debido a que la presencia de estos valores extremos se compensa con un
mayor número de líneas con valores cercanos a la media. No obstante, se debe señalar
que para mejorar la fiabilidad de los análisis intranuales, directamente se ha optado por
descartar los años con un número de líneas inferior a 7.
Fig. 4.4. Desviación estándar frente al número de líneas
2
4
6
8
10
0 5 10 15 20 25 30 35Nº líneas
Zona 4 - Desviación estándar - nº líneas
DesvEst
2
4
6
8
10
0 5 10 15 20 25 30 35Nº líneas
Zona 5 - Desviación estándar - nº líneas
DesvEst
2
4
6
8
10
0 5 10 15 20 25 30 35Nº líneas
Zona 6 - Desviación estándar - nº líneas
DesvEst
30
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
4.3.4 Preparación del análisis de las líneas de costa en DSAS. Construcción
de la línea base y de los transectos
El análisis de la evolución de la línea de costa en la franja litoral de los tres tramos
estudiados se realiza en base a los resultados obtenidos en la aplicación DSAS (Digital
Analysis Shoreline System), desarrollada por el United States Geology Service (Thieler
et al. 2009), y utilizada como complemento del software SIG ArcMap (ESRI). Esta
herramienta se basa en el empleo de las líneas de costa referenciadas espacial y
temporalmente, para a partir de ellas calcular parámetros de variabilidad que permitan
comprender la dinámica evolutiva de la costa tanto intranual como interanualmente dentro
del periodo de estudio. Para un correcto desarrollo del análisis de las líneas de costa en la
aplicación DSAS se debe tener en cuenta, tal y como señalan Thieler et al. (2009), los
siguientes elementos básicos:
a) Geodatabase
Todos los datos y ficheros empleados en DSAS deben gestionarse desde una geodatabase
personal, que actúa como un fichero en el que se almacenan los resultados generados, a
la vez que alberga todos los datos de entrada. Se crea una geodatabase para cada una de
las zonas de estudio, en las que se inserta un shape con las líneas de costa de esa zona y
la línea base generada a partir de ellas. La geodatabase también albergará el fichero con
los transectos una vez se genere desde DSAS.
Para insertar los shapes con las líneas de cada año en la geodatabase, en primer lugar se
crea una nueva geodatabase personal desde el menú ArcCatalog. A continuación se
cargan las capas en el proyecto y, posteriormente, se insertan en la geodatabase como
clases de entidad.
b) Líneas de costa
Se incluyen en una “clase de entidad” todas líneas de costa que vayan a ser objeto del
análisis, agrupándolas por ejemplo las líneas de un mismo año para hacer el análisis
intranual, a la vez que se importa el fichero a la geodatabase personal de la zona
Para el correcto funcionamiento de DSAS resulta indispensable, por una parte, que las
tablas de atributos de los ficheros de entrada sigan la estructura definida por la aplicación,
incluyendo en un nuevo campo el formato de fecha (DATE_) solicitado por el programa
(mes/día/año) (Fig. 4.5). Al mismo tiempo, debemos incluir en las tablas de datos una
31
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
columna que haga referencia a la incertidumbre de los datos, que en este caso, tal y como
se ha indicado en el apartado de obtención de los datos, es de 5m. Este es el error medio
cuadrático estimado asociado a la metodología de obtención de las imágenes seguida,
estimada en trabajos anteriores de Pardo et al., (2012).
Fig. 4.5. Requerimientos de las líneas de costa. Fuente: Thieler et al., 2005.
c) Línea base, baseline o línea de referencia.
Para realizar el análisis de las líneas de costa resulta necesario definir una línea localizada
paralelamente a ellas, y que se emplea como referencia para medir la evolución o cambios
registrados entre las diferentes líneas de costa. Al mismo tiempo, esta línea se toma como
base a la hora de generar los transectos. En este caso se opta por definir la línea base tierra
adentro, lo cual se consigue mediante la realización de un buffer (se opta por 150m) sobre
el conjunto de las líneas de costa. El resultado se recorta, dejando únicamente el
fragmento de línea localizado tierra adentro, y que finalmente constituirá la baseline (Fig.
4.6).
Fig. 4.6. Buffer sobre el conjunto de líneas del sector 5 para obtener la línea de referencia (Izq.) y conjunto de líneas de
costa del sector 4 (azul), la línea de referencia (verde) y los transectos (rojo) (drcha.).
32
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
d) Transectos.
Con ayuda de DSAS definimos una serie de transectos perpendiculares a las líneas de
costa, distribuidos de forma equidistante a largo de la baseline, y que permiten medir las
diferencias de distancia registradas entre las intersecciones con las diferentes líneas de
costa y la línea base (Fig. 4.6). Resulta necesario definir su localización, en este caso “on
shore” o tierra adentro, de manera que los transectos serán creados en la dirección opuesta
para así poder intersectar a las diferentes líneas de costa.
Se ha optado por distribuir los transectos con una separación de 100m entre ellos,
asegurándonos de que los transectos, de 500m de longitud, tienen una longitud suficiente
para intersectar todas las líneas de costa. Este fichero, al igual que las líneas de costa y la
baseline, también se inserta en la carpeta de la geodatabase.
Fig. 4.7. Generación de los transectos perpendiculares a la baseline, y preparación del análisis.
Así se obtiene un conjunto de líneas de costa referenciadas espacial y temporalmente,
localizadas en paralelo a la baseline, y cruzadas por los diferentes transectos (Fig. 4.6)
Cuando ya se dispone de todos los datos indispensables se pasa a ejecutar DSAS.
33
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.8. La distancia desde la baseline y la fecha de las líneas de costa permite a DSAS realizar el cálculo de los estadísticos de las tasas de variación. Fuente: Thieler et al., 2005.
4.3.5 Filtrado y eliminación de errores
A la hora de proceder en nuestro análisis es importante destacar que hemos realizado una
comprobación y filtrado de los datos. Una vez obtenidas las líneas de costa pasamos a
realizar una comprobación de los valores obtenidos, llevando a cabo una tarea de filtrado
y de eliminación de valores anómalos y erróneos. En primer lugar detectamos que en las
zonas con desembocaduras y “golas” algunas de las líneas de costa se desvían siguiendo
los cursos de agua que las originan, por lo que en caso de coincidir estas zonas con alguno
de los transectos se opta por eliminarlo y descartar los resultados que podría llevar
asociados.
También se comprueba que a lo largo de los tramos de costa estudiados no aparezcan
zonas incompletas en cuanto al número de líneas se refiere, faltando datos de algunas de
las fechas, ya que ello alteraría el resultado de los análisis, haciendo necesario descartar
esas zonas.
Por otra parte, es importante destacar que en algunos tramos, como al sur del tramo
localizado más al sur, se ha observado que la línea de costa describía formas claramente
antinaturales, localizándose los puntos alternamente en la interfaz tierra/agua y mar
adentro. Este error, posiblemente asociado a la rotura de las olas, ha hecho que se opte
por no considerar el segmento de playa donde esto ocurre de forma generalizada (Fig.
4.9).
34
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.9. Error registrado en la detección de la línea de costa en el sector 6.
4.3.6 Cálculo de resultados estadísticos mediante DSAS
Para realizar los cálculos en DSAS se procede a cargar las líneas de costa pertenecientes
a un mismo sector, tanto incluyéndolas a todas ellas en un shape conjuntamente para
poder estudiar la evolución del sector, como incluyendo las pertenecientes a un mismo
año para analizar los cambios intranuales. Asimismo se cargan la línea de referencia y los
transectos. Una vez cargada toda la información necesaria, en los menús de DSAS se
seleccionan estos ficheros necesarios para el análisis, especificando en el caso del shape
de la línea de costa el campo que contiene la fecha (DATE_) y la incertidumbre.
Finalmente, en el menú de cálculo de estadísticos de cambio seleccionamos los módulos
Shoreline Change Envelope (SCE), Net Shoreline Movement (NSM), End Point Rate
(EPR), Linear Regression Rate (LRR) y Least Median of Squares (LMS) (Fig. 4.10). El
programa calcula los resultados estadísticos, que aparecen en formato de tabla (.dbf)
dentro de la geodatabase personal en la que se está trabajando.
Fig. 4.10. Selección de los estadísticos calculados en el análisis mediante DSAS.
Se generan dos tablas de resultados para cada grupo de líneas de costa analizadas (Fig.
4.11). La primera de ellas, conocida con el nombre “intersect” recoge la información
referente a los puntos en que los transectos intersectan con las líneas de costa
pertenecientes a diferentes fechas. Este fichero incluye información de la distancia a la
línea de referencia a la que se encuentra cada una de las líneas en cada uno de los
transectos, incluyendo las coordenadas en que se produce la intersección.
35
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.11. Resultados del análisis de las líneas de costa realizado en DSAS. La primera tabla (“Intersect”) muestra los puntos de intersección entre las líneas de costa y los transectos, mientras que la segunda (“Rates”) muestra un resumen
de los parámetros estadísticos.
Por otra parte, la segunda de las tablas de atributos, llamada “rates”, presenta los
resultados de los diferentes parámetros estadísticos calculados por la aplicación. En el
análisis con DSAS se emplean los estadísticos descritos por Thieler et al. (2005)
siguientes, por considerar que son los que mejor se ajustan a los objetivos del estudio y
aportan una información de mayor interés:
– Net Shoreline Movement (NSM): Este estadístico muestra la distancia recorrida entre
la línea de costa en la fecha más antigua y la más reciente, es decir, su movimiento neto
(Fig. 4.12). Puede resultar de utilidad para definir una tendencia global.
Fig. 4.12. Ejemplo del procedimiento de cálculo del estadístico NSM por parte de DSAS. Fuente: Thieler et al., 2005.
- Line Regresion Rate (LRR): Este estadístico refleja el ajuste de una recta de regresión
por mínimos cuadrados a todos los puntos del litoral de un transecto en particular, siendo
la pendiente de la línea la que refleja la tasa de cambio. Se obtiene del cruce de las
posiciones de las líneas de costa frente a las diferentes fechas, respetando su cronología,
36
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
y realizando una regresión lineal simple de los puntos (Fig. 4.13). El valor de LRR suele
ser un complemento útil del estadístico NSM, confirmando su tendencia. Los resultados
positivos van asociados a tendencias acumulativas, mientras que los negativos se suelen
asociar a zonas con una tendencia erosiva.
– Shoreline Change Envelope (SCE): Representa el cambio total sufrido por la línea de
costa, como una variación máxima independientemente de los años en que se produce.
– End Point Rate (EPR): Es una tasa, fruto de dividir la distancia entre la línea de costa
más reciente y la más antigua por el tiempo transcurrido entre los años estudiados.
Fig. 4.13. Ejemplo del procedimiento de cálculo del estadístico LRR por parte de DSAS. Fuente: Thieler et al., 2005.
Pese a que el análisis con DSAS permite obtener otros parámetros estadísticos como ECI
(Confidence of End Point Rate), LSE (Standard Error of the Estimate), y R2, la mayoría
de ellos reflejan fenómenos similares a SCE, EPR, LRR y NSM, descritos anteriormente,
o bien describen su ajuste incertidumbre. De este modo, los cambios en la línea de costa
pueden explicarse a través de estos últimos 4 parámetros, por lo que se decide no entrar a
considerar el resto de parámetros al considerarlo innecesario para nuestros objetivos.
Los resultados obtenidos en formato de tabla pasan a ser analizados en un procesador de
cálculo tipo Excel, donde se calcula el promedio de los diferentes parámetros estadísticos
para cada uno de los años en que se disponen líneas de costa, así como para el total del
periodo 1984-2014.
37
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Estos resultados ayudan a entender cuál es la evolución sufrida por la costa tanto a lo
largo de un año como a más largo plazo, así como determinar si existen patrones erosivos
o acumulativos en los diferentes sectores estudiados. La información recogida en los
parámetros estadísticos que resulta de mayor interés se representa en forma de gráficas y
mapas qué ayudan a interpretar visualmente la tendencia sufrida a lo largo de toda la
costa. A nivel de gráficos se estudia la variación de NSM y LRR a lo largo de todo el
periodo estudiado, mientras que se han elaborado mapas con los valores de NSM y LRR
en los diferentes transectos,
De forma complementaria se analiza la variación de la línea de costa para el conjunto de
cada zona estudiada, mediante el cálculo de la distancia de la línea de costa a la baseline
en el conjunto de transectos de cada zona. Esto se realiza a partir de la tabla “intersect”
de resultados de DSAS, tomando como referencia la distancia existente a línea de costa
durante la primera fecha de la que se tienen datos, a partir de la cual se determina si la
variación sufrida ha sido positiva o negativa, alejándose o acercándose a la baseline.
Otro parámetro que es posible obtener a partir de las líneas de costa es la variabilidad
existente entre las líneas pertenecientes a un mismo año, medida como la desviación
estándar de los cambios en la distancia existente entre cada una de las líneas y la línea
base, en la intersección con los diferentes transectos. Estas desviaciones estándar se
obtienen con la ayuda del cálculo de estadísticas que el mismo ArcMap posee a partir de
las tablas “intersect”, para todos los años del periodo de estudio utilizando la herramienta
“summarize”.
4.2 Granulometrías del sedimento de las playas
El estudio sedimentológico del proyecto aparece centrado en la realización de análisis
granulométricos del sedimento presente en las diferentes zonas de litoral que son objeto
de estudio, dada la relación existente entre la morfodinámica de las playas y sus
características sedimentológicas. Mediante estos análisis se pretende estudiar la
frecuencia con que aparecen granos de diferentes dimensiones en las playas estudiadas.
La granulometría puede aportar información sobre los procesos de erosión y transporte
que tienen lugar en el entorno, así como sobre el origen de los sedimentos.
En este caso se busca analizar los materiales de las playas, arenosos en todos los casos.
Para alcanzar este objetivo, el método propuesto es la separación de los diferentes
tamaños con la ayuda de una columna de tamices.
38
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
4.2.1 Muestreo
Se ha llevado a cabo una toma de un total de 36 muestras (18 en los sectores 4, 5 y 6)
mediante diferentes salidas de campo llevadas a cabo durante el mes de abril de 2015, en
diferentes transectos de las playas, localizados cerca de los extremos de cada uno de los
3 sectores estudiados, y tomando la arena en tres puntos de cada uno de ellos, de forma
perpendicular a la línea de costa: en la zona de rompiente, en el estrán, y al pie de duna
(o en caso de que la playa no presente, en la zona más alejada de la orilla) (Tabla 4.3, Fig.
4.14). La localización de los puntos de muestreo también obedece a la elección llevada a
cabo por Sanjaume (1985), intentando un acercamiento a sus puntos de muestreo para así
poder realizar comparaciones entre los resultados de ambos estudios. Se ha optado por
realizar únicamente dos transectos en cada uno de los sectores estudiados, ya que éstos
tienen una longitud relativamente corta, sin elementos como desembocaduras,
formaciones rocosas o construcciones antrópicas que a priori podrían generar un cambio
en la naturaleza de los sedimentos a lo largo del segmento.
Las muestras se han tomado siempre a una profundidad cercana a los 10-20 cm de forma
manual, retirando previamente el material superficial (Sanjaume, 1985).
Tabla 4.3. Puntos de muestreo sobre el rompiente de la playa.
Transecto Coordenadas
Zona Transecto Playa X Y
4 Norte L'Auir (Gandia) 743920 4324125
4 Sur Gandia Nord 745435 4321755
5 Norte Piles Norte 749690 4315590
5 Sur Piles Sur 750530 4314625
6 Norte Terranova (Oliva) 752575 4312400
6 Sur L'aigua blanca (Oliva) 754290 4310725
4.2.2 Procesado de las muestras
El procesado las muestras se realizó en el laboratorio de la Facultad de Geografía e
Historia de la Universidad de Valencia. En primer lugar se llevó a cabo el secado de las
muestras, mediante la deposición de la arena sobre papel de filtro, en cubetas de aluminio,
exponiéndolas a temperatura ambiente y al sol, hasta eliminar completamente la humedad
(Fig. 4.15).
39
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.14. Los transectos se localizan perpendicularmente a la línea de costa a partir de los puntos indicados en los mapas
Fig. 4.15. Toma de muestras de sedimento (Izq.) y secado de las mismas en el laboratorio (drcha.). Autor: Carlos Cabezas.
Fig. 4.16. Torre de tamices. Autor: Carlos Cabezas.
40
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Una vez secas, 100g de cada una de las muestras han sido tamizados durante 15 minutos
mediante una columna de tamices (Fig 4.16), dispuestos con una luz de malla cada vez
menor, permitiendo una separación de las distintas fracciones de sedimento y conocer en
qué porcentaje se encuentran. A partir de los diferentes diámetros de luz de malla de los
tamices, expresados en milímetros, se calcula el valor en la escala de Φ (Krumbein,
1934).
Los tamices empleados cuentan con una luz de malla que va desde 1400 µm (-0,5 Φ) de
la parte superior hasta los 63 µm (4 Φ) de la inferior (Tabla 4.4). Una vez realizado el
tamizado se procede a pesar los sedimentos que han quedado retenidos por los diferentes
tamices. Considerando estas cantidades retenidas, junto la cantidad de sedimento exacta
que se ha introducido en la columna (algo superior a 100 g), se calculan los porcentajes
de sedimento que se encuentran de cada tamaño.
Tabla 4.4. Luz de malla de los tamices expresada en µm y Φ.
Diámetro (µm) Φ
1.400 -0,5
1.000 0,00
900 0,15
710 0,5
400 1,32
224 2,16
180 2,5
125 3
90 3,5
63 4
4.2.3 Estudio de los resultados
Los resultados se han representado mediante histogramas de frecuencia y curvas de
distribución acumuladas. El programa GRADISTAT (Blott y Pye, 2001) ha sido
empleado en esta etapa del proceso, agilizando el proceso de obtención de resultados. La
nomenclatura empleada para interpretar los parámetros es la de Folk y Ward (1957),
obteniendo diferentes índices granulométricos de interés al analizar cambios en el
sedimento (Friedman y Sanders, 1978).
41
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Los parámetros considerados en el análisis del sedimento principalmente el valor medio
del tamaño del sedimento (2), el grado de clasificación, la curtosis (3), y el grado de
asimetría (4), ya que son comúnmente utilizados para describir la textura y forma del
grano, pudiendo ayudar a determinar cambios y evolución del sedimento (Batman y
Dougherty, 1997).
(2) (3) (4)
Fórmulas para el cálculo de: tamaño medio del sedimento (2), asimetría (3) y curtosis (4), donde P es el diámetro de grano en el percentil x (Blott y Pye, 2001).
Se presta especial atención a la media del tamaño de las partículas en las diferentes zonas
de la playa, ya que es este el parámetro que se considerará al realizar un análisis conjunto
con las pendientes y las variaciones en la línea de costa intranualmente. Se ha optado por
elegir este parámetro ya que representa con bastante eficacia el tamaño de sedimentos que
se encuentra en cada una de las zonas.
4.3 Cálculo de las pendientes de las playas
El análisis del comportamiento de la línea de costa en el presente estudio aparece
complementado con datos referentes a las pendientes de las playas estudiadas, poniendo
especial atención sobre el frente de playa, es decir, sobre los primeros metros, la zona de
rompiente o swash.
Para ello partimos del modelo digital de elevaciones (MDE) suministrado por el Institut
Cartogràfic Valencià (terrasit.gva.es/es/descargas). Este modelo, en formato raster y de
1m de precisión, fue obtenido en el año 2009 a partir de datos LIDAR, y muestra la
elevación de cada uno de los píxeles.
Dado que los ficheros cubren una superficie muy grande de terreno, su manejo en el
programa SIG es lento y poco eficiente, más aun considerando que tan sólo resulta
interesante para este propósito la información contenida en una pequeña franja en la
interfaz tierra/agua. De este modo, se opta por recortar los ficheros, realizando
previamente un buffer de 150m en torno a las líneas de costa, que será empleado de molde
en un “Clip” o recorte posterior del MDE (Fig. 4.17).
42
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.17. Buffer sobre el MDE de 150m en torno a las líneas de costa.
Debido a que los diferentes sectores presentan longitudes no muy elevadas, se considera
innecesario realizar una división en diferentes tramos más pequeños. La reducción de
tamaño de ficheros realizada permite, ahora sí, trabajar con fluidez suficiente en el
programa GIS.
En este punto se pasa a calcular la pendiente a lo largo de la superficie de las playas
mediante la herramienta “Slope”.
Fig. 4.18. Curvas de nivel obtenidas cada 10 cm a partir del MDE. Permiten determinar posteriormente la línea de costa.
Asimismo, para poder delimitar con mayor exactitud la extensión de las playas se plantea
una extracción de las curvas de nivel. Para ello, se utiliza en el software SIG la
herramienta “Contour”, generando curvas de niveles cada 10 cm a partir de los datos de
altitud del MDE (Fig. 4.18).
Las curvas de nivel obtenidas nos ayudan a trazar la línea de la orilla y, al mismo tiempo
y con ayuda de una ortofoto, delimitamos toda la extensión de la playa, hasta zonas con
vegetación, cordones dunares o construcciones. Al delimitar la playa en sus límites con
el mar y tierra adentro, conociendo los valores de altitud en ambas zonas, es posible llevar
a cabo un cálculo de la pendiente total de la misma. No obstante, resulta especialmente
43
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
útil conocer la pendiente existente en los primeros metros de playa, en el llamado
beachface o frente de playa. Estudios previos realizan el cálculo de la pendiente
considerando los primeros 5 m de playa (Pardo et al., 2014), por lo que también se toma
esta distancia al considerarla como representativa. Para poder considerar la pendiente
existente en estos primeros metros se vuelve necesario generar una línea localizada a 5 m
de distancia, lo que se realiza mediante un buffer a partir de la línea de costa determinada
anteriormente, de cuyo resultado únicamente se emplea la línea paralela generada tierra
adentro (Fig. 4.17).
Fig. 4.19. Definición de la totalidad de la superficie de la playa (Izq.). Línea de costa y buffer de 5m a su alrededor, que permite obtener la pendiente del beach face.
Una vez delimitada la superficie total de la playa y el beach face mediante líneas, se toman
los transectos generados en DSAS, equidistantes y perpendiculares a la línea de costa,
para determinar los puntos en que se tomarán datos de altitud. Para facilitar la obtención
de la información se genera una capa de puntos, en la que éstos se localizan sobre las
intersecciones de las capas de líneas de transectos, línea de costa, beach face y beach
back. Una vez localizados los puntos, la herramienta “Extract by point” permite recopilar
los datos de altitud de forma rápida.
A partir de los datos de altitud se resta la altitud de la línea de costa a la de las líneas de
beach face y beach back. De los resultados, por trigonometría se obtienen las pendientes
a lo largo de la playa (1), expresadas como tangente de α, en las intersecciones con los
diferentes transectos, tanto para para los primeros 5m de la playa como para la totalidad
de la misma, es decir, hasta el pie de duna.
tan α = 𝐻𝑡−𝐻𝑚
5 (1) siendo Ht la altura 5m tierra adentro o al pie de duna, y Hm la altura en la línea de costa
44
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 4.20. Puntos a 5m (verde) y sobre la línea de costa (morado), cuya altitud se utiliza para calcular la pendiente de la
playa. Los transectos generados por DSAS (amarillo), localizados cada 100m, sirven de guía.
5. Resultados
Los resultados se muestran agrupados en tres apartados: en primer lugar se muestra el
análisis realizado sobre los cambios registrados por la línea de costa, seguidamente se
muestra el resultado del análisis sedimentológico, para por último presentar los resultados
del estudio de las pendientes de las playas.
5.1 Evolución de las líneas de costa
En primer lugar se analizan los cambios sufridos por la línea de costa en las tres zonas
estudiadas, tanto a nivel interanual considerando la totalidad del periodo 1984 – 2014,
como a nivel intranual, a partir de los valores estadísticos obtenidos con la herramienta
DSAS.
5.1.1 Cambios interanuales. Evolución de las líneas de costa a largo plazo
Se ha optado por presentar los valores obtenidos a lo durante el periodo 1984-2014 de
NSM y LRR en los diferentes transectos de la zona 4 (Gandia), 5 (Piles) y 6 (Oliva) (Fig.
5.1 – 5.6). En el caso del NSM se diferencian por colores los lugares en que el movimiento
neto de la línea de costa presenta valores positivos (acumulación) y negativos (erosión),
señalando asimismo los puntos en que el NSM presenta valores máximos y mínimos de
erosión y de acreción, mientras que en el caso de LRR, se distingue entre los valores,
positivos, negativos, máximos y mínimos del estadístico.
45
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
También se ha plasmado la evolución de los estadísticos NSM a lo largo del periodo
1984-2014, mostrando sus valores anuales fruto de analizar la variación de las líneas
existente en cada uno de los años (Fig 1-3 del Anexo III).
Por otra parte se muestra la variación media de la línea de costa en cada una de las de las
zonas, considerando de forma conjunta la posición en todos los transectos. Para ello se
calcula el cambio de posición de la línea de costa a lo largo de las fechas en que se dispone
de datos, tomando como posición inicial de referencia la de la primera de las fechas con
datos, en 1984. Estas figuras ayudan a comprender cómo ha tenido lugar la erosión o
acreción y qué tendencia han seguido a lo largo del tiempo, en lugar de mostrar
únicamente un resultado para el periodo global (Fig. 5.7).
Realizando un análisis por zonas, se determina que el comportamiento seguido es el
siguiente:
Zona 4 - Gandia
Para la zona número 4, que incluye las playas de Gandia Nord y l’Auir, se cuenta con un
total de 37 transectos, optando por eliminar el localizado más al norte al presentar
anomalías en los datos asociadas a la desembocadura del río Vaca. Atendiendo a los
valores registrados de NSM se puede decir que la zona ha sufrido dos tendencias
enfrentadas, de forma que mientras que la parte norte del segmento estudiado ha sufrido
globalmente un proceso erosivo, la parte sur, más cercana al puerto de Gandia, ha gozado
de un proceso acumulativo, que ha alcanzado máximos de 14 m de acreción en la parte
central de la playa de Gandia Nord. La comparativa con los valores de LRR, positivos en
la parte sur del segmento y superiores a los 0,25 m/año, con los máximos en los transectos
más cercanos al puerto, refuerza la idea de la existencia de acreción en esta zona, mientras
que los valores de LRR negativos, aunque de una magnitud inferior a los 0,25 m/año en
la mayoría de los casos, obtenidos en la parte norte, complementan la idea de que los
segmentos localizados más hacia el norte, en la playa de l’Auir, sí sufren de cierto
retroceso de la línea de costa.
La figura 5.7 considera la evolución fecha a fecha de la línea de costa, considerando toda
la zona en su conjunto. Se observa que las oscilaciones son muy marcadas entre las
posiciones de la costa tomadas a lo largo de un año, destacando que se llegan a alcanzar
retrocesos puntuales, asociados al periodo entre mitad del otoño y mitad de la primavera.
46
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Según la evolución registrada, la zona presenta un patrón de cierta acreción, aunque éste
no es muy acentuado si se la considera en conjunto.
Fig 5.1. Valores de NSM en la zona 4 considerando el periodo 1984 – 2014.
47
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 5.2. Valores de LRR (m/año) en la zona 4 considerando el periodo 1984 – 2014.
48
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Zona 5 – Piles
Sobre la playa de Piles se ha realizado el análisis a partir de 17 transectos, optando por
descartar el localizado más al sur, influenciado por la desembocadura de un barranco que
desvirtuaba los resultados obtenidos.
Se trata de una playa oscilante a nivel longitudinal, lo cual se traduce en la inexistencia
de una tendencia evolutiva clara de la playa, alternando transectos con valores de NSM
positivos y negativos (Fig. 5.3). Al mismo tiempo, el valor de LRR es negativo en todos
los transectos, registrando la mayoría de ellos valores asociados a un retroceso mayor a
los 0,25 m/año (Fig 5.4).
Analizando la evolución de la línea de costa para la zona en conjunto a lo largo del periodo
de estudio (Fig 5.7) se observa una estabilidad en la línea de costa desde el primero de los
datos, en 1984, hasta la década del 2000, con un marcado retroceso en el año 2002,
seguido de estabilidad hasta el año 2009, a lo largo del cual se registró una acreción en la
playa para, posteriormente, volver a posiciones anteriores asociadas a un cierto retroceso,
aparentemente interrumpido de nuevo durante el año 2013 donde temporalmente se
registró una progradación. Analizando el periodo en conjunto parece haberse registrado
una tendencia erosiva moderada en esta zona desde los años 80.
Zona 6 – Oliva
Respecto a la zona 6, que engloba diferentes playas al sur del puerto de Oliva, se ha optado
por eliminar el grupo de transectos localizados más al sur, disponiendo para el análisis
final de los 30 transectos localizados más al norte.
A nivel general parece existir una tendencia acumulativa, aunque no muy marcada, pues
la mayoría de transectos presentan valores de NSM positivos, con magnitudes mayores
que los negativos (Fig. 5.5), siendo el máximo valor de acreción de 16m (localizado en
la parte norte de la zona), mientras que el de máxima erosión de tan sólo 4m.
No obstante, los valores de LRR son negativos en casi todos los transectos estudiados,
apareciendo una tendencia marcadamente diferente entre la parte norte del segmento, en
cuyo extremo aparecen los únicos transectos con LRR positivo, y los LRR negativos
tienen valores pequeños, mientras que en la parte sur la totalidad de los transectos presenta
valores de LRR más marcadamente negativos, que irían asociados a tendencias erosivas
(Fig 5.6).
49
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Analizando la evolución de la costa en promedio, se observa una acreción durante los
años 2000-2001 y 2006-2007, para posteriormente a lo largo de 2010 y 2011 sufrir un
retroceso, y aun aumento a lo largo del 2013 (Fig. 5.7).
Fig 5.3. Valores de NSM en la zona 5 considerando el periodo 1984 – 2014.
50
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 5.4. Valores de LRR (m/año) en la zona 5 considerando el periodo 1984 – 2014.
51
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 5.5. Valores de NSM en la zona 6 considerando el periodo 1984 – 2014.
52
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 5.6. Valores de LRR (m/año) en la zona 6 considerando el periodo 1984 – 2014.
53
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 5.7. Variación media de la línea de costa (m) en cada una de las zonas a lo largo de las fechas en que se dispone de
datos, en el periodo 1984-2014. Las fechas se muestran en formato (mes/día/año).
54
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
5.1.2 Cambios intranuales. Variación de las líneas de costa durante un año.
A nivel intranual, los cambios que sufre la línea de costa a lo largo de un año han sido
analizados estadísticamente mediante el programa DSAS, obteniéndose los parámetros
estadísticos de desviación estándar y SCE, los cuales aparecen representados en las Fig.
5.8 y 5.9. Se han considerado únicamente los años en los que se disponía de más de 7
líneas de costa, con el objetivo de aumentar la fiabilidad.
Los valores estadísticos intranuales de SCE y NSM, para las tres zonas estudiadas
aparecen recogidos en el Anexo III, Tabla 1 y Fig. 1-3.
Evolución de SCE y la desviación estándar
Resulta importante estudiar la magnitud de la variación de las líneas de costa a nivel
intranual, es decir, a lo largo del periodo de un año.
Esto es posible mediante el análisis los valores obtenidos de SCE, que indica la distancia
máxima entre las líneas de costa registradas en un mismo año, así como en la desviación
estándar de la posición de las líneas de costa (Fig. 5.8).
Al comparar los valores de la desviación estándar entre las diferentes playas se puede
observar cómo éstos son netamente inferiores para la zona 5 (Piles) que para las otras dos
playas. De éstas, la zona 4 (Gandia) presenta una desviación estándar mayor que la 6
(Oliva).
Algo similar ocurre con los cambios registrados por el SCE, siendo la zona 5 la que menor
SCE registra, seguida por la zona 6 y, finalmente, la zona 4, con valores claramente más
altos que en los otros 2 tramos de playa (Fig. 5.9).
55
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 5.8. Valores promedio de la desviación estándar en las Zona 4, 5 y 6 en el periodo 1984-2014.
Fig 5.9. Valores promedio de SCE (m) en las Zona 4, 5 y 6 en el periodo 1984-2014.
Dado que posteriormente se pretende establecer una relación entre estos datos y la
pendiente de la costa, se ha optado por tomar el valor promedio de la desviación estándar
entre los años 2007 y 2011 (Tabla 5.1) ya que son los más cercanos en el tiempo a la fecha
de toma de los datos de altitud por parte del LIDAR, en el 2009, que se han empleado
para el cálculo de las pendientes.
2
3
4
5
6
7
8
9
10
1982 1986 1990 1994 1998 2002 2006 2010 2014
(m)
Año
Desviación estándar en zonas 4, 5, 6
Zona 4 Zona 5 Zona 6
10
15
20
25
30
35
40
1982 1986 1990 1994 1998 2002 2006 2010 2014
(m)
Año
SCE en zonas 4, 5, 6
Zona 4 Zona 5 Zona 6
56
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 5.1. Desviación estándar (en m) intranual de cada una de las zonas estudiadas durante el
periodo 1984-2014.
Zona
Año 4 - Gandia 5 - Piles 6 - Oliva
2014 5,32 3,79 4,34
2013 9,22
2012
2011 6,51 4,77 5,44
2010 7,29 7,50 7,01
2009 6,84 4,70
2008 7,55
2007 7,50 6,77 7,72
2006 7,39
2005 8,36
2004 7,73 6,35 8,08
2003 8,58
2002 8,36 5,42 4,61
2001 8,68 6,77 7,21
2000 7,56
1987 8,74 4,65 4,48
Total 7,71 564 6,11
2007-2011 7,14 5,93 6,73
Se puede observar cómo la desviación estándar promedio de la zona 4, en las playas de
Gandia, es la más alta de todas las zonas estudiadas (7,14 m), seguida por la zona 6 (6.73
m) que incluye las playas de Oliva, mientras que la zona 5 es la que menor desviación
presenta entre sus líneas de costa a lo largo de un año (5,93 m).
Se puede comprobar que, pese a que existen cambios marcados en la desviación registrada
entre los diferentes años, la relación de desviaciones considerando la totalidad de los años
de que se dispone datos es similar a la existente al considerar únicamente el periodo 2007
- 2011.
57
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
5.2 Sedimentología
El análisis del sedimento se ha realizado mediante análisis granulométricos, cuyos
resultados aparecen plasmados en las figuras 5.10 – 5.15, así como en la Tabla 5.2.
Dada su importancia y utilidad posterior, en la tabla 5.2 se muestran los valores promedio
de la media del tamaño sedimento en las zonas de rompiente y estrán de cada una de las
zonas estudiadas.
Las figuras muestran la distribución del tamaño de las partículas en la zona de rompiente,
estrán, y pie de duna, en los transectos de muestreo llevados a cabo en la parte norte y sur
de las zonas estudiadas. Se incluyen los porcentajes acumulados que representan cada una
de las clases granulométricas (en µm). También se han elaborado gráficas que muestran los
porcentajes de cada clase granulométrica sin acumular (Anexo II, Fig. 3 y 4).
A nivel general, encontramos materiales más gruesos en la zona de rompiente que en el
estrán y que a pie de duna, a excepción de la zona 6 (Oliva), donde el grosor de los
sedimentos a pie de duna supera a los del estrán. En todos los casos se puede hablar de
arena fina o, en algunos puntos de la rompiente, de arena media, considerando que son
arenas bien clasificadas.
Fig. 5.10. Distribución del sedimento según el tamaño de grano en el transecto norte de la zona 4, en la playa de l’Auir.
0
20
40
60
80
100
0 200 400 600 800 1000 1200 1400
Po
rce
nta
je a
cu
mu
lad
o (
%)
Diámetro de partícula (µm)
Zona 4 norte - L'Auir
Rompiente
Estrán
Pie de duna
58
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 5.11. Distribución del sedimento según el tamaño de grano en el transecto sur de la zona 4, en la playa de Gandia Nord.
Fig. 5.12. Distribución del sedimento según el tamaño de grano en el transecto norte de la zona 5, en la playa de Piles
Fig. 5.13. Distribución del sedimento según el tamaño de grano en el transecto sur de la zona 5, en la playa de Piles.
0
20
40
60
80
100
0 200 400 600 800 1000 1200 1400
Po
rce
nta
je a
cu
mu
lad
o (
%)
Diámetro de partícula (µm)
Zona 4 sur - Gandia Nord
Rompiente
Estrán
Pie de duna
0
20
40
60
80
100
0 200 400 600 800 1000 1200 1400
Po
rce
nta
je a
cu
mu
lad
o (
%)
Diámetro de partícula (µm)
Zona 5 norte - Piles norte
Rompiente
Estrán
Pie de duna
0
20
40
60
80
100
0 200 400 600 800 1000 1200 1400
Po
rce
nta
je a
cu
mu
lad
o (
%)
Diámetro de partícula (µm)
Zona 5 norte - Piles sur
Rompiente
Estrán
Pie de duna
59
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 5.14. Distribución del sedimento según el tamaño de grano en el transecto norte de la zona 6, en la playa de Terranova (Oliva).
Fig. 5.15. Distribución del sedimento según el tamaño de grano en el transecto sur de la zona 6, en la playa de l’aigua blanca (Oliva).
A continuación se realiza el análisis siguiendo los transectos de norte a sur. En primer
lugar se muestran los resultados del transecto norte de la zona 4, en la playa de l’Auir.
(Fig. 5.10). El valor medio del tamaño de sedimento es de 247.3 µm en la zona de
rompiente, 243.6 µm en el estrán y 212.8 µm a pie de duna de modo que, tal y como
cabría esperar, los sedimentos son más finos a medida que nos alejamos del agua.
Respecto a la curtosis, el sedimento de rompiente y estrán es platicúrtico, es decir,
presenta un reducido grado de concentración alrededor de los valores centrales, mientras
que en el punto de muestreo localizado a pie de duna la arena es mesocúrtica, con un
grado de concentración medio alrededor de los valores centrales de la variable.
0
20
40
60
80
100
0 200 400 600 800 1000 1200 1400
Po
rce
nta
je a
cu
mu
lad
o (
%)
Diámetro de partícula (µm)
Zona 6 norte - Terranova (Oliva)
Rompiente
Estrán
Pie de duna
0
20
40
60
80
100
0 200 400 600 800 1000 1200 1400
Po
rce
nta
je a
cu
mu
lad
o (
%)
Diámetro de partícula (µm)
Zona 6 sur - L'aigua blanca (Oliva)
Rompiente
Estrán
Pie de duna
60
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
El siguiente transecto, al sur de la zona 4, se localiza en la playa de Gandia (Fig. 5.11).
En su zona de rompiente aparecen sedimentos con una media de 207.1 µm, en el estrán
de 182.7 µm, y 174.7 µm en el extremo de la playa, junto al paseo marítimo. Los puntos
de rompiente y junto al paseo marítimo presentan un sedimento mesocúrtico, mientras
que el estrán es leptocúrtico, con un elevado grado de concentración de sedimentos
alrededor de los valores centrales. El sedimento a pie de duna se considera simétrico,
mientras que en la zona de rompiente y el estrán los valores son negativos.
El transecto al norte de la zona 5 (Fig. 5.12) muestra unos valores de tamaño medio de
294.7 µm para la rompiente, 220.7 µm para el estrán, y 221.1 µm para el pie de duna. Los
sedimentos en estrán y junto al paseo marítimo son leptocúrticos, y presentan valores de
asimetría negativos, mientras que en el rompiente son simétricos y mesocúrticos.
Al sur de la zona 5, en la zona de Piles que presenta cordón dunar (Fig. 5.13), los
sedimentos de la rompiente presentan una media de 242.2 µm, 235.7 µm en el estrán, y
222.7 µm a pie de duna. Al igual que en el otro transecto realizado en la misma zona, la
asimetría en estrán y pie de duna es negativa, mientras que en la zona de rompiente existe
simetría. Aquí rompiente y estrán presentan sedimentos platicúrticos, mientras que a pie
de duna son mesocúrticos.
El transecto al norte de la zona 6 (Fig. 5.14), localizado en la playa de Terranova, en
Oliva, muestra valores medios de 251.0, 187.9 y 212.1 µm respectivamente para
rompiente, estrán y pie de duna, siendo por tanto éstos últimos mayores que los del estrán.
Existe simetría en rompiente y estrán, mientras que a pie de duna los valores son
negativos. En cuanto a la curtosis, estrán y pie de duna presentan sedimentos
leptocúrticos, mientras que la zona de rompiente los presenta platicúrticos.
Por último, el transecto norte de la zona 6, en la playa de l’Aigua Blanca (Fig. 5.15),
presenta una media de 251.4, 189.6 y 243.0 µm en rompiente, estrán y pie de duna. Los
sedimentos presentan simetría excepto a pie de una, donde ésta presenta valores
negativos, y tanto ésta zona como la rompiente presentan una curtosis de distribución
platicúrtica, mientras que en el estrán es leptocúrtica.
61
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 5.2. Resumen de la media del sedimento promediada para cada una de las zonas
Zona Media en la rompiente (µm) Media en el estrán (µm)
Zona 4 227.20 213.19
Zona 5 268.48 228.19
Zona 6 251.20 188.79
Tabla 5.3. Resumen del análisis sedimentológico de las zonas estudiadas, de las muestras
tomadas en el la parte norte (1) y sur (2) en la zona de rompiente (A), estrán (B) y pie de duna
(C),donde se muestran los valores de media (mean), clasificación, asimetría (skewness) y curtosis
(kurtosis), obtenidos según el método de Folk y Ward (1957).
Zona 4 - Gandia
Transecto Norte - Ahuir Sur - Gandia Nord
(µm) 4_1A 4_1B 4_1C 4_2A 4_2B 4_2C
MEAN 247.3 243.6 212.8 207.1 182.7 174.7
SORTING 1.397 1.391 1.395 1.397 1.329 1.284
SKEWNESS 0.000 -0.042 0.157 0.176 0.158 0.096
KURTOSIS 0.855 0.858 0.964 0.977 1.128 1.092
Zona 5 - Piles
Transecto Norte Sur
(µm) 5_1A 5_1B 5_1C 5_2A 5_2B 5_2C
MEAN 294.7 220.7 221.1 242.2 235.7 222.7
SORTING 1.348 1.395 1.372 1.360 1.367 1.360
SKEWNESS 0.079 0.143 0.182 0.095 0.141 0.207
KURTOSIS 1.135 0.949 1.040 0.883 0.899 1.059
Zona 6 - Oliva
Transecto Norte - Terranova Sur - L'Aigua blanca
(µm) 6_1A 6_1B 6_1C 6_2A 6_2B 6_2C
MEAN 251.0 187.9 212.1 251.4 189.6 243.0
SORTING 1.359 1.303 1.365 1.350 1.308 1.347
SKEWNESS -0.001 0.032 0.173 0.017 0.051 0.135
KURTOSIS 0.864 1.302 1.223 0.860 1.318 0.886
62
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
5.3 Cálculo de la pendiente de las playas
En este apartado se muestra la información de las pendientes de las diferentes zonas de
playa estudiadas, obtenidas en los diferentes transectos definidos mediante el programa
DSAS, a partir de datos LIDAR de agosto de 2009 de las altitudes de la línea de costa y
de los diferentes puntos de la playa.
De este modo se obtuvieron, de norte a sur a lo largo de los diferentes segmentos de
costa, las altitudes de la línea de costa, pie de duna y de la línea que delimita el frente de
playa, localizada 5 m tierra adentro paralela a la línea de costa.
A partir de los datos de las altitudes han obtenido las pendientes de la totalidad de la playa
y del frente de playa, representados en las Fig. 5.16 - 5.18.
Fig. 5.16. Pendiente en la zona 4 entre la línea de costa y la línea de beach face, 5m tierra adentro (azul), y de la totalidad de la playa (naranja).
No obstante, la pendiente del frente de playa se considera mucha más representativa del
comportamiento morfodinámico de la playa, y que guarda una relación más estrecha con
la posición de la línea de costa. Esto ha hecho que se opte por emplear en el análisis
posterior estos valores de pendiente y no los obtenidos a partir de las altitudes a pie de
duna. Al mismo tiempo, la pendiente de la totalidad de la playa parece verse en gran
medida determinada por su anchura, lo que a su vez depende del grado de antropización,
existiendo diferencias entre las zonas con presencia de cordón dunar o paseo marítimo,
además de las realización de tareas de limpieza y redistribución de la arena.
0.00
0.02
0.04
0.06
0.08
0.10
0.12
0.14
0.16
0.18
0 500 1000 1500 2000 2500 3000 3500 4000
Pen
die
nte
(tg
án
gulo
)
Distancia al norte (m)
Pendiente Zona 4
Beach face 5m
Pie de duna
63
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 5.17. Pendiente en la zona 5 entre la línea de costa y la línea de beach face, 5m tierra adentro (azul), y de la totalidad de la playa (naranja).
Fig. 5.18. Pendiente en la zona 6 entre la línea de costa y la línea de beach face, 5m tierra adentro (azul), y de la totalidad de la playa (naranja).
Con el objetivo de intentar establecer una relación entre las pendientes de las playas y la
variación de la línea de costa que sufren anualmente, se necesita trabajar con un único
valor por playa. De este modo, en cada una de las zonas se ha obtenido el valor promedio
de la pendiente en el frente de playa (Tabla 5.4).
Tabla 5.4. Pendiente media de las zonas 4, 5 y 6 en el año 2009, expresada como tangente del
ángulo entre la altitud de la línea de costa y la altitud a 5m de la misma.
Zona Pendiente media 2009 (tangente α)
4 0.070
5 0.106
6 0.050
0.00
0.02
0.04
0.06
0.08
0.10
0.12
0.14
0.16
0.18
0 200 400 600 800 1000 1200 1400 1600
Pen
die
nte
(tg
án
gulo
)
Distancia al norte (m)
Pendiente Zona 5
Beach face 5m
Pie de duna
0.00
0.02
0.04
0.06
0.08
0.10
0.12
0.14
0.16
0.18
0 500 1000 1500 2000 2500 3000 3500
Pen
die
nte
(tg
án
gulo
)
Pendiente Zona 6
Beach face 5m
Pie de duna
64
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Se puede observar fácilmente como los valores de pendiente del frente de playa en la
zona 5 (Piles), con un promedio de 0,106, son superiores a los registrados en las otras dos
zonas. A continuación encontramos la zona 4 (Gandia), que presenta un valor promedio
de 0,070, mientras que en la zona 6 (Oliva) éste es de 0,050.
6. Discusión
6.1 Análisis de la evolución de las líneas de costa
Las 3 zonas de costa estudiadas presentan un comportamiento a largo plazo que puede ser
analizado mediante una interpretación de los resultados mostrados en el presente trabajo.
Este patrón evolutivo se encuentra influenciado tanto por fenómenos que condicionan a
la totalidad del sistema litoral o a la célula sedimentaria en que se encuentran, como por
elementos y alteraciones que condicionan el comportamiento del medio litoral de forma
puntual.
Zona 4 - Gandia
En primer lugar, la zona 4, constituida por las playas de l’Auir (mitad norte) y Gandia
Nord (mitad sur), parece mostrar un comportamiento diferente entre su parte norte y su
parte sur. En la mitad norte parece tener lugar un patrón erosivo débil, con un retroceso
medio que en la mayoría de los casos no supera los 0,25 m/año (Fig. 5.2). Por su parte, la
mitad sur presenta una tendencia acumulativa de hasta 0,56 m/año, con avances de la línea
de costa máximos de 14 m entre 1984 y 2014 (Fig. 5.1).
Considerando la línea de costa de toda la zona en conjunto, el análisis del cambio de
posición a lo largo de las fechas de que se dispone de información (Fig 5.7) muestra que
la zona presenta un patrón de cierta acreción. La línea de costa sufrió un proceso
progradante registrado a partir del año 2000, con un importante retroceso al final del 2001,
del que se fue recuperando en los años siguientes. En el año 2009 la línea sufrió un
retroceso muy marcado, apareciendo en varios momentos del año hasta 20 m desplazada
tierra adentro respecto a su posición original en 1984, produciéndose su recuperación
durante la primavera de 2011 y 2013. El año 2014, el último de la serie, muestra una línea
localizada en una posición similar a la registrada al inicio del estudio en 1984.
Esta zona de costa se encuentra condicionada a nivel antrópico tanto por la presencia del
puerto de Gandia a poca distancia hacia el sur (Fig. 3.1), que actúa como freno al
65
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
transporte longitudinal de sedimentos N-S que en este sector se produce, así como por los
diferentes movimientos de arena que se han realizado en la zona (Tabla 3.2, Fig. 3.8). Las
diferentes actuaciones de extracción y vertido de arenas llevadas a cabo en esta zona
pueden haber condicionado el comportamiento de la línea de costa. Destacan las
operaciones de dragado realizadas en la playa de l’Auir, hasta en 4 ocasiones, varias de
ellas destinadas a alimentar las playas del sur de la Safor (Fig. 3.9), como Bellreguard y
Piles.
Otro factor condicionante del comportamiento de esta zona es la presencia de un cordón
dunar que ocupa la totalidad de la playa de l’Auir, en la mitad norte de este sector (Fig.
3.1) que puede haber actuado como reservorio de arena para las playas de la zona, a la
vez que de factor favorable para la progradación registrada con el paso de los años en la
playa de Gandia Nord.
Al mismo tiempo, se debe destacar que ésta zona costera ha seguido tradicionalmente un
patrón progradante (Sanjaume y Pardo, 2008), en la que la creación y ampliación del
puerto de Gandia no ha hecho más que actuar como un obstáculo al transporte
longitudinal, generando acumulación justo al norte del puerto, generando una anchura de
playa que alcanza los 130 m en algunos puntos.
Zona 5 - Piles
La zona 5, delimitada básicamente por la playa de Piles, presenta un patrón marcadamente
erosivo, que se traduce en unas tasas regresivas superiores los 0,25 m/año en la mayoría
de transectos, con máximos de hasta 0,43 metros/año (Fig. 5.4). No obstante,
considerando la variación registrada entre la primera y la última fecha de las que se
dispone de datos, la posición de la línea de costa se ha desplazado incrementando la playa
hasta 10 m en puntos al sur de la zona, en la que en la que encontramos una pequeña
formación dunar (Fig. 3.1). No obstante, los valores de LRR, así como el análisis de la
posición de las líneas de costa y de ortofotos históricas, lleva a pensar que éstas
acumulaciones puntuales no responden a un patrón acumulativo en esta zona, si no que
se deben a oscilaciones temporales de la línea de costa a nivel longitudinal, resultando
evidente que se trata de una zona de marcada tendencia erosiva.
Analizando la variación de la línea de costa de toda la zona en conjunto, a lo largo de las
diferentes fechas de que se dispone de información de la posición (Fig 5.7), se puede
observar cómo la posición de la línea de costa se mantuvo más o menos constante desde
66
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
1984 hasta el año 1990, observando un importante cambio a partir del año 2000, cuando
se aprecia una progradación que dura hasta el año 2002, a partir del cual la costa parece
sufrir un retroceso que se mantiene hasta el año 2009, lo que se explica por los importantes
aportes artificiales realizados durante la década de los noventa. A partir de este momento
la playa parece recuperarse parcialmente, para de nuevo volver al patrón recesivo hasta
el año 2013, en el que durante un periodo de un año la playa sufrió una progradación
importante, dejando la posición final cercana a la registrada en 1984.
Esta zona de costa, localizada al sur del puerto de Gandia, aparece enormemente
condicionada por la presencia del mismo, ya que tras sus sucesivas ampliaciones impide
el paso de materiales hacia las playas del sur (Sanjaume y Pardo, 2008). La playa
también parece estar muy determinada por la existencia de un paseo marítimo construido
sobre el cordón dunar que existía anteriormente en esa zona (Pardo y Sanjaume, 2001),
lo que ha generado la inmovilización de una gran cantidad de arena, eliminando el
almacén natural o reservorio de la misma. Además, la construcción de un paseo, una
estructura sólida, a una distancia reducida de la línea de costa, no hace más que
incrementar la acción erosiva.
El extremo sur de esta zona presenta una pequeña formación dunar, frente a la cual la
playa presenta una anchura notablemente superior a la zona que cuenta con paseo
marítimo, así como una pendiente menor considerando el conjunto de la playa (Fig. 5.16)
debido a que la línea de costa ha retrocedido, pero la playa ha podido encontrar el
equilibrio gracias a que el frente costero en esta zona no has sido inmovilizado, y a que
el cordón dunar actúa como suministro de arena en casos de temporal, permitiendo el
mantenimiento de una anchura mayor.
Zona 6 - Oliva
La zona 6, localizada a lo largo de las playas al sur del puerto de Oliva, presenta una
situación erosiva, bastante elevada en su parte sur, en la playa de l’Aigua Blanca, con la
totalidad de los transectos mostrando valores de LRR asociados a retrocesos mayores a
0,25 m/año, y un máximo de 0,43 m/año. Por su parte, el extremo norte de la zona 6, en
la playa de Terranova, cuenta con los únicos dos transectos de LRR positivo, con el
máximo acumulativo únicamente en 0,12 metros por año (Fig. 5.6). Atendiendo a los
valores de NSM considerando el periodo de 1984 a 2014, la zona presenta una alternancia
67
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
de transectos con erosión y con acreción, aunque predomina la erosión en la parte sur,
mientas que la acreción es más visible en la parte norte (Fig. 5.5).
Analizando la evolución de la costa en promedio (Fig. 5.7), se observa una acreción
durante los años 2000-2001 y 2006-2007, un retroceso a lo largo del 2010, una
recuperación a lo largo del 2011, y una importante acreción durante el 2013, presentando
en 2014 valores cercanos a los iniciales. Se observa un comportamiento similar al sufrido
en la zona 5.
Considerando el conjunto de datos analizados sobre esta zona, se puede determinar que
la tendencia seguida no es muy marcada, aunque parece existir un proceso ligeramente
erosivo, especialmente en la parte sur de la zona, posiblemente motivado por una falta de
transporte de arena asociado a la deriva longitudinal, la cual ya es prácticamente nula en
esta zona (Pardo et al., 2014).
Análisis conjunto
Analizando la posición media de la línea de costa para el conjunto de las tres zonas (Fig.
5.7), se puede observar que las oscilaciones a lo largo de un año son muy marcadas
(MAGRAMA, 2007) destacando que se llegan a alcanzar los mayores retrocesos
habitualmente entre mitad del otoño y mitad de la primavera.
Resulta llamativo el comportamiento similar de las diferentes zonas en algunos momentos
de la serie histórica. En primer lugar se observa la acreción experimentada por las 3 zonas
en la década del 2000 respecto a los datos de los años 80 y 90. Esta progradación sufrida
aparentemente por la costa, que se ve reflejada al comparar los datos del año 2000 con
los de 1990, podría deberse al vertido en 1993 de 600.000 m3 de arenas al sur del puerto
de Gandia, tomadas de fondos submarinos, con el objetivo de regenerar las playas
localizadas al sur del puerto (Tabla 3.2), una actuación de una gran magnitud que generó
un cambio profundo en la configuración morfológica de ese sector. La gran alimentación
que entre 1996 y 1999 se realizó en las playas del Saler, al sur del puerto de Valencia, al
norte de estas tres zonas (Anexo I, Tabla I) también puede haber ayudado a dar lugar a
esta acreción. En el caso de Gandia, la acreción iría asociada a la acumulación generada
por la presencia del puerto.
También resulta significativo el transporte durante el año 2009 de un gran volumen de
arena desde la zona 4, en concreto desde la playa de l’Auir, hacia las playas del sur de la
Safor, entre las que se encuentra la de Piles (zona 5). Esta época coincide con retrocesos
68
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
evidentes registrados en la zona 4 (Fig. 5.7), lo cual podría suponer una evidencia de la
estrecha relación que hay entre las operaciones de dragado y el retroceso de la costa.
Asimismo, los vertidos llevados a cabo sobre la playa de Gandia Nord coinciden en su
localización con los transectos donde el parámetro NSM registra los mayores avances y
retrocesos, lo cual volvería a reforzar la idea de que existe una relación directa entre las
actuaciones humanas de este tipo y cambios en la posición de la línea de costa.
Por su parte, en la zona 5 se han realizado vertidos de arena de forma repetida en la década
de los años 2000 y 2010 (Tabla 3.2) con el objetivo de frenar el retroceso de la línea de
costa, ya que la playa, estrecha, ha llegado a desaparecer en determinados temporales
(Fig. 6.1). Estas actuaciones parecen tener una repercusión directa sobre la posición de la
línea de costa, destacando los grandes vertidos llevados a cabo durante la primavera de
2009 sobre los municipios de Bellreguard y Miramar, al norte de esta zona, y sobre el
propio municipio de Piles, precedidas por otra actuación importante a principios de 2008
también sobre la playa de Piles, que presumiblemente juegan un papel esencial en la
recuperación sufrida por la playa a lo largo de ese año. Este proceso de acreción a lo largo
del 2009 es opuesto al registrado en la zona 4, donde tal y como se ha comentado la playa
sufrió un importante retroceso en el mismo periodo de tiempo.
Asimismo, a finales del 2001 y principios del 2002 en todas las zonas se produce un
retroceso muy importante, del cual las playas parecen tardar en recuperarse. Este hecho
lo asociamos con la presencia de un fuerte temporal en noviembre de 2001 sobre el
Mediterráneo occidental (Ramis et al. 2002; Julià-Ventura, 2003; Obiol, 2003;
www.puertos.es), y que junto con una serie de sucesivas tormentas que tuvieron lugar
hasta la primavera del 2002 (Pardo et al., 2014), provocó una gran movilización de
sedimento e impulsó la erosión de la costa. Tal y como indican Pardo et al. (2014) estas
zonas de costa, con orientación noroeste – sudeste sufrieron la acción de las tormentas
sucesivas, lo que generó un retroceso progresivo que alcanzó su máximo en mayo de 2002
(Fig 5.7). Este fenómeno también resulta apreciable en las zonas 1, 2 y 3 (Fig. 12, Anexo
I).
De forma similar, aunque no tan acentuada, en octubre del año 2007, a finales del año
2009 y principios de 2010, y en octubre de 2010, se registran importantes retrocesos de
la línea de costa en las zonas 4, 5 y 6, algunos de los cuales, como el de octubre de 2010,
también con repercusión en las zonas 1, 2 y 3 (Anexo I, Fig 12). Estos retrocesos, al igual
que ocurre durante el temporal de noviembre de 2001, van asociados con periodos en los
69
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
que se registran niveles muy elevados del nivel del mar (Anexo II, Fig II), lo que da lugar
a que la posición de la línea de costa se localice temporalmente tierra adentro.
En cuanto a la evolución a largo plazo, el litoral de la Safor, una zona históricamente
acumulativa parece haber sufrido un cambio de tendencia, tal y como apuntan (Sanjaume
y Pardo, 2008). Las zonas 4 y 6 deberían presentar claros patrones acumulativos, algo que
sin embargo no ocurre, ya que por ejemplo la zona de Oliva, localizada al final de la célula
sedimentaria, muestra una tendencia de ligera erosión, lo que podría actuar como una
evidencia de una falta de sedimento en el sistema litoral. Esta suposición aparece
reforzada por estudios sobre el cordón dunar, como en el caso de Oliva (Sanjaume y
Pardo, 2008) con retrocesos y pérdidas significativas en el volumen de sedimentos.
No obstante, el ser humano actúa sobre el medio litoral acentuando los procesos erosivos
o, por el contrario, enmascarándolos, reduciendo o anulando las tendencias negativas
(Sanjaume y Pardo, 2008). Tal es el caso de la zona 5 (Piles), donde la erosión no alcanza
la magnitud esperada y frena su tendencia en diferentes momentos gracias a los vertidos
llevados a cabo por el ser humano.
Las actuaciones humanas pueden suponer, como en el caso de los puertos de Gandia y
Valencia, un freno al transporte longitudinal, generando acumulaciones de sedimento en
su parte norte y erosión en las playas del sur, algo que se manifiesta en las zonas 5 (Piles)
y 2 (Anexo I, Fig. 8, 9 y 12) Estas actuaciones pueden tener repercusiones en zonas muy
alejadas (Sanjaume y Pardo, 2008), y ser las causantes de cambios de tendencia o de
ausencia de sedimento en zonas distantes como en el caso de la zona 6, en Oliva.
Finalmente se debe resaltar los efectos que pueden tener actuaciones llevadas a cabo fuera
del medio litoral, como en la red fluvial, frenando los aportes de sedimento al medio
litoral e impidiendo su entrada en el sistema (Sanjaume y Pardo, 2008).
De este modo, encontramos una dinámica costera y unos patrones evolutivos que en gran
parte no son el resultado de una tendencia natural, si no que se encuentran determinados
por las actuaciones humanas. Estas actuaciones en muchos casos responden a decisiones
políticas y económicas, cuyo coste económico, ecológico y social debería ser evaluado
(Brocal et al., 2001).
Las principales consecuencias son la alteración de las tendencias naturales de las playas,
en muchos casos generando estrechamientos que pueden conllevar finalmente la
desaparición de las playas.
70
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
La erosión de la playa conduce a una reducción de la anchura, limitando la superficie
disponible para los usuarios, que a su vez se traduce en una disminución de su confort y
disfrute de la playa (Yepes y Medina, 2007). El estrechamiento o, en el peor de los casos,
desaparición de las playas, va a tener importantes consecuencias negativas sobre el sector
turístico, el cual es básico en la economía de la zona.
Fig. 6.1. Desaparición de la playa de Piles en algunos puntos. 13/11/2004 Fuente: Google Earth, 2015.
Asimismo, si el patrón recesivo natural se mantiene, nuevas actuaciones antrópicas de
movilización de arenas van a ser necesarias, con el consiguiente coste que ello conlleva,
ya que de lo contrario las playas más estrechas como la de Piles corren el riesgo de
desaparecer, algo que ya ha ocurrido en alguna ocasión (Fig. 6.1).
Frente a los procesos erosivos, resulta de vital importancia el mantenimiento de dunas así
como evitar la inmovilización de los frentes costeros. La arena presente en los cordones
dunares, como al norte de Gandia, y al sur de Piles y Oliva, puede contrarrestar o paliar
la pérdida de aporte sedimentario en estas zonas. Parece existir una mayor afección del
proceso erosivo en lugares donde no existen cordones (Pardo y Sanjaume, 2001),
generando diferentes magnitudes de cambio pese a contar con un tipo de sedimentos
similar (Pardo et al., 2014), En momentos críticos para la playa, como los temporales, el
mar erosiona las dunas, que liberan sedimento al sistema costero (Pye y Blott, 2008)
disminuyendo la acción de las olas y ayudando a mantener la anchura de playa (Davidson,
2010).
71
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
6.2 Análisis sedimentológico
En este apartado se discuten los resultados sedimentológicos obtenidos para las zonas 4,
5 y 6 (Apartado 5.2) realizando una comparativa con los resultados obtenidos sobre las
mismas zonas por Sanjaume (1985).
Con el propósito de realizar esta comparativa, la tabla 6.1 muestra los resultados
obtenidos por Sanjaume (1985), de norte a sur, en los análisis sedimentológicos llevados
a cabo en puntos localizados en las zonas 4, 5 y 6. La localización está disponible en el
Anexo II, Fig. 5. Los valores originales han sido transformados a µm para facilitar su
interpretación.
Tabla 6.1. Análisis sedimentológico llevado a cabo por Sanjaume (1985) en las zonas estudiadas
Sector Muestra Media (µm) Asimetría Curtosis Zona
4 -
Gan
dia
195 237.66 -0.039 1.465 Rompiente
196 214.34 -0.213 0.813 Rompiente
197 191.44 0.046 0.889 Estrán
5 -
Pile
s
211 262.61 0.166 1.541 Rompiente
212 231.33 -0.185 1.322 Estrán
213 227.19 -0.227 1.231 Estrán
214 236.19 -0.163 1.35 Estrán
6 -
Oliv
a
219 204.76 -0.109 0.736 Rompiente
220 198.2 0.035 0.743 Estrán
221 194.66 0.100 0.766 Estrán
Comparando los datos de la tabla anterior con la tabla 5.3, la playa de l’Auir, al norte de
la zona 4, parece haber sufrido un aumento en el tamaño de sedimento en la zona de
rompiente (actualmente 247.3 µm) mientras que la playa de Gandia Nord, al sur, ha
registrado una disminución en el tamaño medio, tanto en estrán (207.1 µm) como en
rompiente (182.7 µm), frente a los valores de principios de los años 80.
En cuanto a la zona 5, en su extremo norte el sedimento de rompiente se ha vuelto más
grueso (294.7 µm). Para el resto de la playa se dispone de datos del estrán, siendo el
tamaño actual (235.7 µm) similar al que existía en el pasado.
Por su parte, en la zona de Oliva resulta llamativo el incremento de grosor en el sedimento
de rompiente en la zona norte (251.0 µm), mientras que para la totalidad de la playa el
sedimento del estrán es ligeramente más fino.
72
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Analizando los datos históricos y actuales, parecen confirmarse algunos patrones. En
primer lugar, se puede interpretar que resulta más común encontrar sedimentos finos en
playas como Gandia Nord (Tabla 5.3) y en la Malva-Rosa en su tramo más cercano al
puerto (Anexo I, Tabla 3), donde predomina la acumulación, en este caso ligada a la
presencia de obstáculos al transporte longitudinal, por lo tanto con condiciones más
favorables para la sedimentación.
Por el contrario, procesos erosivos van asociados a un hidrodinamismo más alto y a un
incremento en el tamaño de los sedimentos, especialmente en la zona de rompiente. Esto
ocurre en las zonas 5 (Piles) y 2. Las playas de l’Auir y Oliva, tal y como se ha indicado
en el apartado 6.1, aparentemente presentan ligeros patrones erosivos, algo que parece
estar apoyado por aumentos en el tamaño del sedimento en la rompiente de ambas playas.
El incremento del tamaño de rompiente en Oliva resulta especialmente llamativo, y
estaría asociado a un progresivo lavado de los materiales más finos, los cuales no están
siendo reemplazados por otros nuevos, constituyendo un síntoma de la falta de sedimento
en el sistema, al tratarse del tramo final del mismo. Este hecho confirmaría la suposición
de que se está produciendo un cambio de tendencia en la evolución seguida por la costa.
El análisis de los datos actuales también muestra cómo dentro de las distintas zonas, el
tamaño de sedimento parece ser mayor en la parte norte frente a la parte sur (Tabla 5.3;
Anexo I, tabla 3), posiblemente debido a la movilización que se produce del sedimento
fino con mayor facilidad hacia el sur, hasta encontrar obstáculos al transporte.
Por otra parte, resulta necesario destacar la gran cantidad de actuaciones antrópicas que
se han llevado a cabo en estos tramos de costa, con importantes volúmenes de vertido de
sedimento durante las últimas décadas. La existencia de alimentación artificial y otras
actuaciones pueden hacer que los resultados se muestren como menos coherentes. Debido
a la artificialización de los aportes de sedimentos (Sanjaume y Pardo, 2008) esto puede
traducirse en una mayor uniformidad de los sedimentos, lo que en el playa de Piles podría
estar relacionado con la presencia de un sedimento platicúrtico.
6.3 Análisis de la relación entre las desviaciones estándar anuales de
las líneas de costa y las características morfológicas de las playas
El comportamiento de las playas parece estar determinado por diferentes factores, entre
los que se encuentra la variación registrada por la línea de costa a lo largo de un año, y
73
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
características morfológicas como la pendiente y la textura, los cuales están relacionados
entre sí (Heitor y Gama, 2009), y pueden determinar la respuesta de la playa ante
fenómenos como los temporales (Haerens et al., 2012; Pardo et al., 2014; Reyes et al.,
1999). El conocimiento de estas relaciones permitiría estimar el comportamiento de un
segmento costero y determinar sus características morfológicas a partir de la variación de
la línea de costa, la cual es posible conocer siguiendo la metodología descrita en este
trabajo, sin necesidad de llevar a cabo análisis sedimentológicos o de pendiente.
6.3.1 Relación entre la pendiente de las playas y la variación intranual de
la línea de costa
Parece existir una relación inversa entre la pendiente media de las playas y la variación
registrada por la línea de costa, de forma que la pendiente condiciona de algún modo la
magnitud de cambios de posición de la línea de costa Pardo et al., 2014.
Con el objetivo de intentar relacionar la pendiente de las playas estudiadas con la
variación registrada por la línea de costa, y así poder responder a la hipótesis inicial del
estudio, se enfrenta la pendiente media registrada en 2009 en las 6 zonas de que se dispone
de datos (Tabla 5.4) con la desviación estándar media de la línea de costa registrada entre
los años 2007 y 2011 (Tabla 5.1).
Se observa la tendencia lógica esperable, apareciendo una relación inversa entre la
pendiente y la desviación promedio de la línea de costa, de forma que las zonas con una
mayor pendiente sufren menores variaciones intranuales de la línea de costa (Tabla 6.2,
Fig. 6.2).
Tabla 6.2.Variación de la línea de costa y pendiente media en las diferentes zonas.
Zona Desviación estándar media (2007-2011) Pendiente media (2009)
1 – La Patacona y Malva-Rosa 7,068 0,044
2 - Els Ferros y la Garrofera 5,702 0,087
3 – L’Alcatí y la Punta 5,036 0,067
4 -Gandia 7,136 0,070
5 - Piles 5,934 0,106
6 - Oliva 6,725 0,047
74
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig 6.2. Pendientes registrada en el año 2009 en las 6 zonas de costa, expresada como tan α, frente a la desviación
estándar del periodo 2007-2011, en m.
El coeficiente de determinación obtenido presenta un valor bajo (r2=0,24), pero se debe
tener en cuenta el hecho de que los datos de pendiente, muy cambiantes, han sido en una
única fecha (agosto de 2009), lo que condiciona el ajuste del resultado.
Anteriores trabajos, como el de Pardo et al. (2014), presentan resultados que concuerdan
con los expuestos en este trabajo, asociando las pendientes suaves a tramos de costa
acumulativos, como los de las zonas 1 y 4, localizados al norte de barreras transversales
como son los puertos Valencia y Gandia respectivamente, en los cuales podemos observar
que la variación de la línea de costa a lo largo del año es máxima (Fig. 6.2).
En cambio, las zonas de costa con pendientes más pronunciadas van asociadas a lugares
con patrones erosivos, como ocurre en las zonas 2 y 5, las cuales a su vez muestran una
variación intranual de la línea de costa menor. Se trata en ambos casos de playas
localizadas al sur de los puertos Valencia y Gandia, importantes obstáculos al transporte
de sedimentos.
6.3.2 Relación entre la textura de las playas y la variación intranual de la
línea de costa
También parece posible establecer una relación entre el tipo de sedimento y la variación
intranual media de la línea de costa en cada una de las zonas estudiadas. Con este objetivo
se consideran los valores de tamaño medio del sedimento, tanto de rompiente como del
estrán (Tabla 6.3), los cuales se enfrentan a las desviación estándar media del periodo
2007-2011 (Fig. 6.3 y 6.4).
1
2
3 4
5
6
y = -0.0139x + 0.1571R² = 0.2441
0.000
0.020
0.040
0.060
0.080
0.100
0.120
4.5 5.0 5.5 6.0 6.5 7.0 7.5
Pen
die
nte
med
ia 2
00
9 (
tg α
)
Desviación estándar media (m)
Variación línea de costa - Pendiente
75
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 6.3. Variación de la línea de costa y tamaño del sedimento en las diferentes zonas.
Zona Desviación estándar media
(2007-2011) Media rompiente (µm) Media estrán (µm)
1 – La Patacona y Malva-Rosa 7,068 190,4 195,4
2 - Els Ferros y la Garrofera 5,702 430,7 250,0
3 – L’Alcatí y la Punta 5,036 269,7 265,1
4 -Gandia 7,136 227,2 213,2
5 - Piles 5,934 268,5 228,2
6 - Oliva 6,725 251,2 188,8
Fig. 6.3. Variación de la línea de costa frente al tamaño de sedimento en la zona de rompiente.
Fig. 6.4. Variación de la línea de costa frente al tamaño de sedimento en el estrán.
1
2
3
4
56
y = -56.172x + 624.96R² = 0.3267
0
50
100
150
200
250
300
350
400
450
500
4.5 5.0 5.5 6.0 6.5 7.0 7.5
Tam
año
med
io s
edim
ento
en
ro
mp
ien
te
(µm
)
Desviación estándar (m)
Variación de la línea de costa - Tamaño del sedimento de rompiente
1
23
45
6
y = -32.41x + 426.53R² = 0.819
0
50
100
150
200
250
300
4.5 5.0 5.5 6.0 6.5 7.0 7.5
Tam
año
med
io s
edim
ento
en
est
rán
(µ
m)
Desviación estándar (m)
Variación de la línea de costa - Tamaño del sedimento del estrán
76
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
A partir de este análisis parece evidente que las playas con materiales más finos presentan
una mayor variabilidad intranual de la línea de costa, cumpliendo la tendencia esperable.
De este modo, las zonas 1 y 4, con tamaño de sedimento menor en el estrán y rompiente,
son las zonas con mayor variabilidad de la línea de costa a lo largo de un año. Por el
contrario, las zonas 2, 3 y 5, que son las que menor desviación estándar registran, son las
que presentan un tamaño de sedimento más grueso (Fig. 6.3).
El análisis de la textura del sedimento frente a la variación de la línea de costa muestra
un coeficiente de determinación (r2=0,33) más elevado que en el caso de la pendiente, de
lo que se deduce que el tamaño del sedimento es un parámetro más estable que la
pendiente de la playa cuyo valor, tal y como se ha señalado anteriormente, además de ser
bastante cambiante, responde a una toma de datos puntual. No obstante, si se atiende al
tamaño del sedimento del estrán, el ajuste es mucho mejor (r2=0,82) (Fig. 6.4),
posiblemente al presentar esta zona una mayor estabilidad, eliminándose del análisis
alteraciones como la presencia de piedras en la zona de rompiente que pueden distorsionar
el análisis.
Atendiendo a la relaciones que existe entre los parámetros analizados, parece posible
afirmar que la variabilidad de las líneas de costa, medida a través de la desviación estándar
media, actúa como un buen indicador de algunas características de las playas, como la
textura y la pendiente, las cuales son mucho más complejas y caras de obtener.
El análisis puede estar desdibujado tanto por errores inherentes a la metodología seguida,
como por variaciones registradas en la morfología y la textura de las playas, asociados
tanto a causas naturales que hayan tenido lugar previamente a la toma de datos, como a
diferentes actuaciones humanas, que en conjunto pueden contribuir a que las relaciones
entre parámetros no se den con la magnitud esperada.
El análisis confirma la hipótesis inicial al encontrarse relaciones entre las diferentes
características morfológicas de las playas y su comportamiento. De este modo, playas con
sedimento más fino presentan variabilidades mayores de la posición de la línea de costa
a lo largo del año, llevando asociadas pendientes menores. Estas situaciones suelen darse
en playas con tendencia acumulativa. Por el contrario, playas con sedimento más grueso
llevan asociadas mayores pendientes y variabilidades menores de la posición de la línea
de costa, llevando asociados patrones más erosivos.
77
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
7. Conclusiones y recomendaciones
Las playas, espacios naturales de gran importancia para la sociedad en el Golfo de
Valencia, sufren en la actualidad procesos de erosión o acreción inducidos por el ser
humano, a causa de una gestión inadecuada del espacio litoral. Para alcanzar una gestión
adecuada, resulta necesario disponer de información y datos suficientes.
En la búsqueda de indicadores del estado y evolución del litoral, el análisis de las líneas
de costa se presenta como una alternativa útil, al permitir detectar tanto variaciones a
corto plazo como tendencias generales, favoreciendo las labores de gestión y
planificación del medio costero. En este sentido, las imágenes Landsat TM junto al
empleo de técnicas de análisis digital permiten determinar la posición de la línea de costa
con gran precisión.
La herramienta DSAS junto con el empleo de software GIS posibilitan el análisis y
procesado de gran cantidad de información referente a la posición de la línea de costa,
determinando variaciones y tendencias evolutivas. El elevado grado de automatización
del proceso y la existencia de información retrospectiva permiten su aplicación a escala
global para realizar comparaciones y observar tendencias evolutivas de la costa a largo
plazo.
La evolución de la costa está condicionada por factores naturales, como el oleaje, las
corrientes, las mareas y nivel del mar, y el viento, además de fenómenos como las
tormentas, con gran capacidad para movilizar sedimento y desplazar la posición de la
línea de costa a corto plazo. No obstante, el ser humano condiciona enormemente los
patrones evolutivos del litoral mediante sus actuaciones, al construir barreras al transporte
longitudinal de sedimentos como los puertos, al realizar extracciones y vertidos de arenas,
y al modificar el régimen hídrico y actuar sobre la red fluvial, imposibilitando la entrada
de sedimento al sistema litoral.
El litoral de la Safor, una zona históricamente acumulativa, parece haber sufrido un
cambio de tendencia en sus playas, lo que podría suponer una evidencia de la falta de
sedimento en el sistema litoral. De las tres zonas consideradas, únicamente la playa de
Gandia Nord, al norte del puerto, presenta una tendencia acumulativa clara, con una
disminución del tamaño de los sedimentos y una pendiente muy baja. En cambio, la playa
de Piles presenta una clara tendencia erosiva, que se traduce en un aumento de su
78
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
pendiente, mayor tamaño de sedimento, y estrechamiento de la playa, que hace necesarias
operaciones de vertido de arenas.
Los procesos erosivos que algunas playas registran son paliados mediante actuaciones de
gran coste económico. La tendencia seguida por algunas playas conduce a situaciones de
riesgo, en las que sin la continua intervención humana podrían de sufrir estrechamientos
importantes o incluso desaparecer. Frente a los procesos erosivos, el mantenimiento de
las dunas, como reservorios de arena, y una correcta planificación y ordenación territorial
resultan básicos.
El estudio del sedimento permite identificar tendencias acumulativas o erosivas en las
playas, así como la deriva longitudinal y la existencia de actuaciones antrópicas, estando
el tamaño y tipo de sedimento condicionados por la realización de actuaciones de vertido.
Los análisis granulométricos y de pendientes llevados a cabo han sido útiles para probar
la hipótesis formulada en el trabajo, aunque para mejorar la fiabilidad sería necesario
trabajar con un volumen de datos mayor, con más puntos de muestreo o réplicas en el
caso de las granulometrías. En cuanto al cálculo de las pendientes, sería recomendable
contar con datos referentes a otras fechas, así como calcular las pendientes considerando
anchuras de frente de playa diferentes.
Al relacionar el comportamiento y la variabilidad de la posición de la línea de costa a lo
largo de un año, la textura y la pendiente de la playa, se ha demostrado que sí existen
estrechas relaciones entre el comportamiento de las playas y sus características
morfológicas. De este modo, playas en las que se producen situaciones de acumulación
y, por lo tanto, aumentan su anchura, presentan menos pendiente, mayor variación de la
línea de costa a lo largo del año, y sedimento más fino, y suelen localizarse al norte de
obstáculos al transporte longitudinal de sedimentos.
También se ha determinado que la variabilidad de la línea de costa, medida a través de la
desviación estándar media, parece actuar como un indicador de otras características de la
playa, como la textura y la pendiente, lo cual resulta de gran interés a la hora de obtener
información sobre las características morfológicas de la playa, ya que éstas son mucho
más complejas y caras de obtener.
La metodología planteada para obtener la variación de la línea de costa resulta
especialmente útil, ya que no sólo permite determinar la evolución de las playas, sino que
79
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
también permite estimar características morfológicas de las playas como la textura y la
pendiente.
Bibliografía
Almonacid, J., Pardo, J. E., Ruiz. L. (2013). Detección automática de la línea de costa con
precisión subpíxel en imágenes Landsat 7 con error de bandeado. XV Congreso de la Asociación
Española de Teledetección INTA, Torrejón de Ardoz (Madrid). 22-24 octubre, 2013.
Ariza, E., Jiménez, J. A., y Sardá, R. (2008). A critical assessment of beach management on the
Catalan coast. Ocean & Coastal Management, 51(2), 141-160.
Batman, S. y Dougherty, E. (1997). Size distributions for multivariate morphological
granulometries: texture classification and statistical properties. Optical Engineering, 36(5), 1518-
1529.
Becker, J. M., Firing, Y. L., Aucan, J., Holman, R., Merrifield, M., y Pawlak, G. (2007). Video‐
based observations of nearshore sand ripples and ripple migration. Journal of Geophysical
Research: Oceans (1978–2012), 112(C1).
Blott, S. Pye, K. (2001). GRADISTAT: a grain size distribution and statistics package for the
analysis of unconsolidated sediments. Earth surface processes and Landforms, 26(11), 1237-
1248.
Brocal, R., López García, M. J., y Pardo, J. E. (2001). Cambios en la línea de costa mediante
fotografía aérea e imágenes IRS-PAN en el litoral valenciano: sector Cullera-Tavernes (1956-
1999). Teledetección. Medio ambiente y cambio global, 225-228.
Cámara oficial de Comercio, Industria y Navegación de Valencia (2011). La Economía de la
Comunidad Valenciana. Disponible en <http://www.camaravalencia.com/es-
ES/informacion/economica/estadisticas_economicas/Documents/La_economia_de_la_CV_2011
.pdf> [Consultado mayo 2015]
Church, J. A. (2001). How fast are sea levels rising? Science, 294(5543), 802-803
Davidson, R. (2010). Introduction to coastal processes and geomorphology. Cambridge
University Press.
DGC, 2008. Directrices sobre actuaciones en playas. Secretaría General para el Territorio y la
Biodiversidad, Dirección General de Costas, Ministerio de Medio Ambiente de España. (41 pp.).
Disponible en: <www.magrama.es/es/costas/publicaciones/directrices_sobre_playas_tcm7-
153279.pdf> [Consultado mayo 2015]
DGC. Informes de la Dirección General de Costas sobre vertidos y movimientos de arena.
Dirección General de Costas, Ministerio de Medio Ambiente de España
80
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fernández, A. (2002). Ordenación del frente litoral de la Albufera: sector Dehesa de El Saler,
Valencia. Ciudades: Revista del Instituto Universitario de Urbanística de la Universidad de
Valladolid, (7), 145-149.
Folk, R., y Ward, W. (1957). Brazos River bar: a study in the significance of grain size
parameters. Journal of Sedimentary Research, 27(1).
Font J, Julià A, Rovira J, Salat J, Sánchez-Pardo J, (1986). Circulación marina en la plataforma
continental del Ebro determinada a partir de la distribución de masas de agua y los
microcontaminantes orgánicos en el sedimento. Actas Geológicas Hispánicas 21-22, 483-489.
Friedman, G. y Sanders, J. (1978). Principles of sedimentology (Vol. 8). New York: Wiley.
Google Maps <http://googlemaps.es> [Consultado abril 2015]
Google Earth (2015). Imagen del 13/11/2004 de la zona de Piles. Escala indeterminada.
[Consultado agosto 2015].
Guizar-Sicairos, M., Thurman, S. T., Fienup, J. R. 2008. Efficient subpixel image registration
algorithms. Optics Letters, 33(2), 156-158.
Haerens, P., Bolle, A., Trouw, K., Houthuys, R., 2012. Definition of storm thresholds for
significant morphological change of the sandy beaches along the Belgian coastline.
Geomorphology 143–144, 104–117.
ICV (Institut Cartogràfic Valencià). MDE-LIDAR (2009). <www.icv.gva.es> [Consultada en
mayo de 2015]
IPCC Working Group II, 2014. Climate Change 2014: Impacts, Adaptation, and Vulnerability.
Intergovernmental Panel on Climate Change.
James, R. J. (2000). From beaches to beach environments: linking the ecology, human-use and
management of beaches in Australia. Ocean & Coastal Management, 43(6), 495-514.
Julià, T. (2003). El temporal de mar de noviembre de 2001 y sus efectos en el litoral catalán.
Universitat Politècnica de Catalunya. Dept. d'Enginyeria Hidràulica, Marítima i Ambiental.
Komar, P. (1976). Beach processes and sedimentation. New Jersey: Prentice Hall.
Krumbein, W. (1934). Size frequency distributions of sediments. Journal of sedimentary
Research, 4(2).
MAGRAMA (2007). Estudio ecocartográfico del litoral de las provincias de Alicante y
Valencia. Memoria divulgativa. Ministerio de Medio Ambiente de España.
MAGRAMA. Guía de playas. Ministerio de Medio Ambiente de España.
<http://www.magrama.gob.es/es/costas/servicios/guia-playas> [Consultado abril 2015]
MOPU (Ministerio de obras públicas y urbanismo). (1970). Playas. Modelos, tipos y sugerencias
para su ordenación. Dirección general de Puertos y Costas. Madrid, 64p.
NASA, 2006. Landsat 7 science data user's handbook.
<http://landsathandbook.gsfc.nasa.gov/pdfs/Landsat7_Handbook.pdf> [Consultado mayo 2015]
Obiol, E. (2003). La regeneración de playas como factor clave del avance del turismo
valenciano. Cuadernos de geografía, (73), 121-146.
81
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Obiol, E. y Pitarch, M. (2011). El litoral turístico valenciano: intereses y controversias en un
territorio tensionado por el residencialismo. Boletín de la Asociación de Geógrafos Españoles,
(56), 177-200.
Ojeda, J. (2000). Métodos para el cálculo de la erosión costera. Revisión, tendencias y
propuesta. Boletín de la Asociación de Geógrafos Españoles, (30), 103-118.
Ojeda. E. (2008): “Shoreline and nearshore bar morphodynamics of beaches affected by artificial
nourishment”. Tesis doctoral. Universidad Politécnica de Cataluña
Ojeda. J, Díaz, M., Prieto, A. y Álvarez. J. (2013). Línea de costa y sistemas de información
geográfica: Modelo de datos para la caracterización y cálculo de indicadores en la costa andaluza.
Investigaciones geográficas, 60, 37-52.
OMT, Organización Mundial del Turismo (2012). UNWTO Annual Report.Diponible en:
< www2.unwto.org/es/annual-reports> [Consultado abril 2015]
Osorio. A, Pérez. J, Ortíz. C y Medina. R (2007). Técnicas basadas en imágenes de video para
cuantificar variables ambientales en zonas costeras. Avances en Recursos Hidráulicos,
(16), 51-64.
Pardo, J.E. (1991). La erosión antrópica en el litoral valenciano. Conselleria de’Obres
Públicques, Urbanisme i Transports, GVA. (240 págs).
Pardo, J. E., y Sanjaume, E. (2001). Análisis multiescalar de la evolución costera. Cuadernos de
geografía, (69), 95-125.
Pardo, J. E., Ruiz, L. A., Palomar, J. M., Calafa., Colmenarez, G. R., Almonacid, J. y Gracia, G.
(2008). Teledetección, GPS y Lidar: Nuevas técnicas e análisis y evolución de la línea de costa y
de los espacios playa-duna. Actas de las Jornadas Técnicas “Las nuevas técnicas de información
geográfica al servicio de la gestión de zonas costeras: Análisis de la evolución de playas y dunas,
62-68.
Pardo, J. E., Palomar, J. M., García, L. y Garrigues, P. (2011). Determinación de la tendencia
evolutiva de un segmento de playa basándose en múltiples levantamientos tridimensionales. En
Montoya I., Rodríguez, I., y Sánchez, M. J. (Eds.), Avances en Geomorfología Litoral. Actas de
las VI Jornadas de Geomorfología Litoral (27–30). Madrid: Universidad Rey Juan Carlos.
Pardo, J. E., Almonacid, J., Ruiz, L.A., y Palomar, J. (2012). Automatic extraction of shorelines
from Landsat TM and ETM+ multi-temporal images with subpixel precision. Remote Sensing of
Environment, 123, 1-11.
Pardo, J. E., Almonacid, J., Ruiz, L. A., Palomar, J., y Rodrigo, R. (2014). Evaluation of storm
impact on sandy beaches of the Gulf of Valencia using Landsat imagery series. Geomorphology,
214, 388-401.
Pye, K., Blott, S.J., 2008. Decadal-scale variation in dune erosion and accretion rates: an
investigation of the significance of changing storm tide frequency and magnitude on the Sefton
coast UK. Geomorphology 102(3), 652–666.
Puertos del Estado. Datos históricos. < www.puertos.es> [Consultado en agosto 2015]
82
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Ramis, C., Alonso, S., Romero, R., y Homar, V. (2002). Análisis preliminar del temporal del 10
al 12 de Noviembre de 2001 en Baleares Preliminar analysis of the 10-12 November 2001 storm
in the Balearics. III Asamblea Hispano-Portuguesa de Geodesia y Geofísica (4-7 Febrero).
Reyes, J.L., Martins, J.T., Benavente, J., Ferreira, Ó., Gracia, F.J., Alveirinho-Dias, J.M., López-
Aguayo, G., 1999. Gulf of Cádiz beaches: a comparative response to storm events. Bol Inst. Esp.
Oceanogr. 15, 221–228.
Sanjaume, E. (1985). Las costas valencianas. Valencia, Universitat de València, 505 p.
Sanjaume, E. y Pardo, J. E. (2008). Cambios de tendencias recientes en la evolución costera del
golfo de Valencia: análisis espaciales y sedimentológicos. Actas de las Jornadas Técnicas “Las
nuevas técnicas de información geográfica al servicio de la gestión de zonas costeras: Análisis
de la evolución de playas y dunas.
Serra, J. (1998): Impacto de las obras de abrigo del puerto de Valencia sobre las playas del sur,
Littoral'98, 361-369
Serra, J. y Medina, J. R. (1996). Beach monitoring program of Valencia (Spain). Coastal
Engineering Proceedings, 1(25).
Shepard. F (1950): “Beach cycles in southern California”, U.S. Army Corps of Engineers, Beach
Erosion Board, Technical Memo, 20.
Thieler, E. R., Himmelstoss, E. A., Zichichi, J. L. y Ergul, A. (2009). The Digital Shoreline
Analysis System (DSAS) Version 4.0-An ArcGIS Extension for Calculating Shoreline Change
(No. 2008-1278). US Geological Survey.
Vargas, M., García, M.C., Moya, F., Tel, E., Parrilla, G., Plaza, F., Lavín, A., García M.J.
(2010). Cambio climático en el Mediterráneo español. Madrid: Instituto Español de
Oceanografía.
Voigt B. (1998). Glossary of Coastal Terminology. Washington State Department of Ecology,
Coastal Monitoring and Analysis Program, 98-105.
Yepes, V. y Medina, J. R. (2005). Land use tourism models in Spanish coastal areas. A case study
of the Valencia region. Journal of Coastal Research, 83-88.
Yepes, V. y Medina, J. R. (2007). Gestión de playas encajadas de uso intensivo. Libro de
Resúmenes de las IX Jornadas Españolas de Ingeniería de Costas y Puertos, 175-176.
83
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Anexos
Anexo I – Datos referentes a las zonas 1, 2 y 3
El presente anexo recoge información perteneciente a las zonas de costa Nº 1, 2 y 3, la
cual aparece recopilada en el TFM de Jesús Soriano. Esta información ha sido utilizada
en el trabajo para elaborar un análisis global más coherente, al poder disponer de una
mayor cantidad de datos.
Las zonas de estudio aparecen localizadas en la Fig 1., incluyendo información referente
a la antropización del frente costero. A continuación, las Fig. 2, 3 y 4 muestran las 6
playas que de norte a sur constituyen los 6 segmentos de costa estudiados.
La tabla 1, junto a la figura 5, muestran las actuaciones de préstamo y vertido de
sedimento llevadas a cabo en las playas estudiadas desde la década de los 90 hasta la
actualidad.
Las Fig. 6 – 11 muestran los resultados estadísticos para los parámetros NSM y LRR en
las 3 zonas estudiadas, empleados para estudiar la evolución seguida por los diferentes
segmentos de playa entre 1984 y 2014. Asimismo, esta evolución también se estudia a
partir de la variación media de la línea de costa de cada una de las zonas a lo largo de las
fechas en que se dispone de datos (Fig. 12). La variación intranual de la línea de costa
aparece reflejada mediante la desviación estándar media de la línea de costa en la Fig. 13,
así como mediante la variación intranual del SCE (Fig. 14).
En cuanto a los resultados de las granulometrías, estos aparecen reflejados en las tablas 2
y 3. La tabla 2 muestra el valor promedio del tamaño medio del sedimento, obtenido en
cada una de las zonas, mientras que la tabla 3 resume el análisis sedimentológico de las
muestras tomadas en los transectos norte y sur de las zonas estudiadas.
Finalmente, la pendiente del frente de playa en las diferentes zonas aparece recogida en
la tabla 4.
84
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
1 - Zona de estudio
Fig. 1 Zonas de estudio 1, 2 y 3. Perpendicularmente a la costa se muestran los transectos
estudiados en las playas. Paralelamente a la costa se muestran la naturaleza del frente costero,
antropizado (en rojo) o con presencia dunar (amarillo)..
85
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 2. Playas de la zona 1: La Patacona (Izq.) y Malva-Rosa (Drcha.). Fuente:
www.magrama.gob.es/es/costas/servicios/guia-playas y www.eltiempo.es.
Fig. 3. Playas de la zona 2: Els Ferros (Izq.) y la Garrofera (Drcha.). Fuente:
www.valencia.es
Fig. 4. Playas de L’Alcatí (Izq.) y la Punta (Drcha.). Fuente: www.valencia.es y
parador.es/es/paradores/parador-de-el-saler
86
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
2 – Actuaciones y movimientos de arena
Tabla 1. Préstamos y Vertidos de sedimento en las playas estudiadas desde 1980-2010.
Fuente: elaboración a partir de datos de la Dirección General de Costas (DGC).
Fig. 5. Dragados y vertidos realizados en las zonas estudiadas durante el período 1984 – 2014.
año comienzo fin DESCRIPCIÓN PRÉSTAMO VERTIDO TIPO VOLUMEN (m3)
1980 1980 1980 Movimiento de arena N puerto de Valencia Restauración Saler Arena 500000
1982 1982 1991 Movimiento de arena Zona interna Restauración Saler Arena -
1983 12/20/1983 8/4/1984 Contrucción dique exento Cantera Playa Pinedo escollera 57796
1988 1988 1998 Movimiento de arena Malladas, laguna racó de l'olla, canteras Restauración Saler Arena -
1996 12/18/1996 4/25/1999 Movimiento de arena Playa Saler (parte Trasera) Playa Saler (Parte activa) Arena 249146
1996 12/18/1996 4/25/1999 Movimiento de arena Playa de la Malvarrosa Playa Saler Arena 210915
1996 12/18/1996 4/25/1999 Movimiento de arena Playa de la Malvarrosa Playa del Saler (Formación núcleo duna) Arena 81014
1997 1997 2000 Movimiento de arena Malladas y N de puerto de Valencia Restauración Saler Arena -
1998 dic-98 Movimiento de arena Gola del Perellonet Playa de Pinedo Arena 29900
1998 dic-98 Movimiento de arena Gola de Pujol Playa de Pinedo Arena 11100
1998 10/30/1998 12/30/1998 Contrucción espigón Cantera Playa de Pinedo escollera 2860.92
1999 1999 1999 Movimiento de arena Gola del Perellonet Playa de Pinedo Arena 27000
1999 1999 1999 Movimiento de arena Gola del Perellonet Playa de Pinedo Arena 8500
1999 1999 1999 Movimiento de arena Gola de Pujol Playa de Pinedo Arena 6500
2001 2001 2007 Extracción arenaMalladas Saler y sobre antiguo Paseo
marítimoRestauración Saler Arena 86183
2006 abr-06 nov-07 Aportación de Arena Excavación fuera DPM-T Pinedo-arbre de Gos Arena 35000
2009 dic-09 nov-10Reacomodación de arenas
en la playaGola del Perellonet Playa de L'abre del Gos Arena 15244
2009 dic-09 nov-10Reacomodación de arenas
en la playaGola de Pujol Playa de L'abre del Gos Arena 5350
2010 2010 2010 Movimiento de arena Cantera Restauración Saler Arena 121431
2010 2010 2010
Crecimiento del espigón Sur
de la desembocadura del
nuevo cauce del río Turia
(Plan Sur). Eliminación de
los espigones preexistentes.
- - - -
2010 2010 2010 Movimiento arenas Golas del Pujol y del Perellonet Restauración Saler - -
2012 nov-12 dic-12Reacomodación de arenas
en la playaGola del Perellonet Playa la Garrofera (El Saler) Arena 13000
2014 dic_2014 dic_2014 Formación dunar Playa Malvarrosa (Valencia) Playa del Saler (Valencia) Arena 20000
87
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
3 – Evolución de la línea de costa
Fig. 6. Net shoreline movement de 1984-2014 en la zona 1.
88
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 7. Linear Regresion Rate de 1984-2014 en la zona 1.
89
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 8. NSM de 1984-2014 en las playas de Els Ferros y La Garrofera, Zona 2
90
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 9. LRR de 1984.2014 en las playas de Els Ferros y La Garrofera. Zona 2
91
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 10. NSM de las playas de L’Alcatí y La Punta, 1984-2014. Zona 3
92
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 11. LRR en las playas de L’Alcatí y La Punta, 1984-2014. Zona 3
93
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 12. Variación media de la línea de costa (m) en cada una de las zonas a lo largo de las
fechas en que se dispone de datos, en el periodo 1984-2014. Las fechas se muestran en
formato (mes/día/año).
94
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 13. Desviación estándar media intranual de la línea de costa en las 3 zonas (1984-2014).
Fig. 14. SCE intranual en las 3 zonas estudiadas (1984-2014).
4 – Sedimentología
Tabla 2. Media del tamaño de sedimento promediada para cada zona.
Zona Media en la rompiente (µm) Media en el estrán (µm)
Zona 1 190.4 195.4
Zona 2 430.7 250.0
Zona 3 269.7 265.1
95
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Tabla 3. Resumen del análisis sedimentológico de las zonas estudiadas, de las muestras
tomadas en el la parte norte (1) y sur (2) en la rompiente (A), estrán (B) y pie de duna (C).
5 – Pendiente del frente de playa
Tabla 4. Pendientes de 2009 en las zonas estudiadas. Se expresan como la tangente del
ángulo.
96
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Anexo II – Información adicional
El presente anexo recoge información adicional a la mostrada en la memoria del TFM,
por no considerarse estrictamente esencial, aunque en ocasiones se hace referencia a estos
datos.
Se recoge información referente a la variación del nivel del mar en la zona estudiada,
mostrada en las Fig. 1 y 2.
Fig. 1. Nivel del mar según el mareógrafo de Gandia en el periodo 2008-2014. Fuente:
www.puertos.es
Fig. 2. Nivel del mar según el mareógrafo de Valencia para el periodo 1984-2014. Fuente:
www.puertos.es
Por otra parte, a nivel sedimentológico se muestran los resultados de los análisis, en este
caso con la distribución de tamaño de las partículas en porcentaje (Fig. 3 y 4).
97
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 3. Distribución de tamaño de las partículas de sedimento en los puntos de muestreo
98
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig.4. Distribución de tamaño de las partículas de sedimento en los puntos de muestreo
99
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Finalmente se muestra el mapa con la localización de los puntos de muestreo utilizados por
Sanjaume (1985) en los análisis de sedimento (Fig. 5).
Fig.5. Localización de los puntos de muestreo en los análisis de Sanjaume (1985).
100
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Anexo III - Parámetros estadísticos DSAS
Este anexo recoge la información adicional obtenida mediante DSAS referente a la
evolución y cambios sufridos por la costa. Esta información no ha sido empleada para
realizar los análisis plasmados en la memoria, o bien ya aparece en la memoria expresada
de otra manera, por lo que se consideraba redundante incluirla.
Las Fig. 1 – 3 muestran la evolución de los estadísticos NSM a lo largo del periodo 1984-
2014, mostrando sus valores anuales fruto de analizar la variación de las líneas existente
en cada uno de los años. Por su parte, la tabla 1 recoge los valores de NSM y SCE
intanualmente a lo largo del periodo estudiado.
Fig. 1. Evolución del estadístico NSM entre 1984 y 2014 en la zona 4.
Fig. 2. Evolución del estadístico NSM entre 1984 – 2014 en la Zona 5.
-30
-20
-10
0
10
20
30
20
14
20
13
20
12
20
11
20
10
20
09
20
08
20
07
20
06
20
05
20
04
20
03
20
02
20
01
20
00
19
87
19
86
19
85
19
84
NSM (m) Zona 4 1984-2014
NSM
-25
-20
-15
-10
-5
0
5
10
15
20
14
20
13
20
11
20
10
20
09
20
08
20
07
20
06
20
04
20
03
20
02
20
01
20
00
19
87
19
85
19
84
NSM (m) Zona 5 1984-2014
NSM
101
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Fig. 3. Evolución del estadístico NSM entre 1984 – 2014 en la zona 6.
Tabla 1. Estadísticos asociados a cambios anuales de las zonas 4, 5 y 6.
Zona 4 Zona 5 Zona 6 AÑO SCE NSM SCE NSM SCE NSM
2014 22.08 -0.55 11.8 8.45 13.4 9.66
2013 37.04 -3.21 17.57 9.1 21.38 11.26
2012 27.59 -7.82 2011 24.7 10.02 15.68 4.85 17.84 7.91
2010 28.02 -3.5 24.27 -18.85 23.85 -18.28
2009 26.91 2.31 16.14 7.09 12.93 0.46
2008 23.34 -8.98 8.42 -0.23 7.82 -3.96
2007 32.81 -0.12 23.82 -8.17 26.09 -16.61
2006 25.41 3.38 20.25 -9.09 2005 29.01 -8.54 5.42 5.03
2004 24.07 -5.52 17.32 0.56 22.27 -6.95
2003 32.27 -15.8 21.81 0.68 22.36 0.04
2002 29.6 16.93 15.3 8.02 12.78 6.24
2001 29.79 -23.92 21.28 -15.51 22.21 -16.75
2000 24.89 -14.63 11.25 -4.86 1987 31.28 24.75 13.22 6.33 12.6 3.07
1986 15.22 -6.88 1985 10.2 9.41 7.35 5.52 7.69 6.43
1984 18.85 -13.89 9.05 -3.83 8.81 -5.79
-20
-15
-10
-5
0
5
10
15
20
14
20
13
20
11
20
10
20
09
20
08
20
07
20
06
20
05
20
04
20
03
20
02
20
01
19
87
19
85
19
84
NSM (m) Zona 6 1984-2014
NSM
102
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Bibliografía de los anexos
Ayuntamiento de Valencia. Disponible en: <www.valencia.es> [Consultado en agosto
2015]
El tiempo.es. Disponible en: <www.eltiempo.es> [Consultado en agosto 2015]
MAGRAMA. Ministerio de Agricultura, Alimentación y Medio Ambiente. Guía de
playas. Disponible en: <www.magrama.gob.es/es/costas/servicios/guia-playas>
[Consultado en agosto 2015]
MAGRAMA. Ministerio de Agricultura, Alimentación y Medio Ambiente.Dirección
General de Sostenibilidad de la Costa. Informes sobre actuaciones de movimientos de
sedimento, dragados y vertidos.
Parador de El Saler. Disponible en: <parador.es/es/paradores/parador-de-el-saler>
Consultado en agosto 2015]
Puertos del Estado. Información histórica. <www.puertos.es> [Consultado en mayo 2015]
103
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Agradecimientos
Me gustaría aprovechar para dar las gracias a las personas que me han ayudado con este
trabajo. A Jesús, por ser quien me ha acompañado a lo largo de todo el proyecto, con
quien he tropezado en los mismos obstáculos y con quien a la fuerza me he tenido que
entender para buscar soluciones. A Josep, mi tutor, por guiarme a lo largo de este trabajo,
y a Francisca, por sus consejos. A León, por echarnos un cable en el laboratorio. A Borja
y Carlos, por sus clases magistrales de GIS. A Adri, por acompañarme a tomar muestras
de arena. A Anna, por sus traducciones. A mis padres, y a los compañeros de piso Mati y
Fani, por aguantar mis cabreos cuando DSAS se ponía cabezón. Gracias a todos.
104
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
Currículum Vitae
Carlos Cabezas Rabadán
>Educación y formación
-Graduado en Ciencias ambientales (UPV). Nota: 9.2 [2009-2013]
-Máster en Evaluación y Seguimiento Ambiental de Ecosistemas Marinos y Costeros
(UPV). Nota: 9.3 [2013-2014]
-Máster en Técnicas para la Gestión del Medio Ambiente y del Territorio (UV)
Cursándolo en la actualidad. Nota media: 8.77 [2014-2015]
>Idiomas
-Español y Valenciano -Inglés: B2 (Intermedio-alto)
>Proyectos
-Propuesta para la adecuación de una vía verde con contenidos divulgativos
entre Gandia y l'Orxa. [2013]
Propuestas para la adecuación de una vía verde siguiendo el trazado de un antiguo ferrocarril, incluyendo elementos para la divulgación ambiental. Trabajo de Final de Grado. Calificación: 9, Excelente.
-Análisis de la relación existente entre las características de las playas y el perfil,
preferencias y percepción de sus usuarios [2014]
Trabajo de investigación que busca determinar la existencia de relaciones entre los perfiles de usuario, tipo de playa al que acuden, y motivación y percepción al visitarla. TFM. Calificación: 10, Mención especial.
>Experiencia
-Prácticas en el departamento de urbanismo del ayuntamiento de Gandía. Programa de prácticas
de la Diputación de Valencia. [Julio - Agosto 2015]
-Prácticas en l’Institut Cartogràfic Valencià (ICV). Análisis espacial para el estudio técnico del
Catálogo de playas de la Comunidad Valenciana, desarrollando una metodología para clasificar
las playas de acuerdo a su naturalidad. [Junio 2014]
-Beca de colaboración del Ministerio de Educación y Ciencia. Colaboración con el Departamento
de Ingeniería hidráulica y medio ambiente (UPV) en el proyecto “Influencia de factores
ambientales en la distribución espacio temporal de la macrofauna bentónica de playas arenosas.” [Noviembre 2013 – Julio 2014]
-Prácticas en empresa. NORAY ingeniería y medio ambiente. [Febrero – Junio 2013]
-Diseño y puesta en práctica de los talleres de la semana de la ciencia en la Escuela Politécnica
Superior de Gandia [Noviembre 2012]
105
Análisis de la línea de costa y su relación con los parámetros morfológicos en playas de la Safor
>Cursos independientes
-Deporte, turismo y medio ambiente (Universidad Politécnica de Valencia). [Abril 2013]
-Biodiversidad en los ecosistemas mediterráneos: su conservación como deber científico y social
(Universidad de Valencia). [Julio 2011]
>Reconocimientos y premios
-Beca a la excelencia académica de la Generalitat Valenciana
-Premio Fin de Grado 2012-13 de la Escuela Politécnica Superior de Gandia
-Premio a los Alumnos de Cuarto Curso con Mejor Expediente de la UPV
-Premio a los Alumnos de Tercer Curso con Mejor Expediente de la UPV
-Premio a los Alumnos de Segundo Curso con Mejor Expediente de la UPV
-Premio a los Alumnos de Primer Curso con Mejor Expediente de la UPV