analisis de la linea de costa y su relacion con los parametros morfologicos en playas de la safor...

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


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


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


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


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


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


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


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


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


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


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


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).

http://www.puertos.es/

13


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


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


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


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


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)


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


Playa de Gandia Bellreguard 50419 744916 4322876 2009/02 2009/04


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





18


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


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


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


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


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


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


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


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


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


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


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


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


DesvEst

2

4

6

8

10

0 5 10 15 20 25 30 35Nº líneas


DesvEst

30


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


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


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


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


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


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


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


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


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


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


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


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).

http://www.google.es/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&cad=rja&uact=8&ved=0CCIQFjAA&url=http%3A%2F%2Fwww.icv.gva.es%2F&ei=459bVYvdKcu3UYGEgagI&usg=AFQjCNHOrQQHfeNRb_7ibRwqpoNqhfpbTw&bvm=bv.93756505,d.d24


42


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


ú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


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


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


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


Fig 5.2. Valores de LRR (m/año) en la zona 4 considerando el periodo 1984 – 2014.

48


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


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).


50



51



52



53


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


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


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


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


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


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 (

%)


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 (

%)


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 (

%)


Zona 5 norte - Piles sur

Rompiente

Estrán

Pie de duna

59


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 (

%)


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 (

%)


Zona 6 sur - L'aigua blanca (Oliva)

Rompiente

Estrán

Pie de duna

60


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


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


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




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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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

http://www.cyclopaedia.es/link?u=http%3A%2F%2Fcgat.webs.upv.es%2Fbigfiles%2FAlmonacid-Caballer_et_al_2013_AET.pdf&b=Subpixel

http://www.cyclopaedia.es/link?u=http%3A%2F%2Fcgat.webs.upv.es%2Fbigfiles%2FAlmonacid-Caballer_et_al_2013_AET.pdf&b=Subpixel

80


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


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]

http://www.puertos.es/

82


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


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


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


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


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


3 – Evolución de la línea de costa

Fig. 6. Net shoreline movement de 1984-2014 en la zona 1.

88


Fig. 7. Linear Regresion Rate de 1984-2014 en la zona 1.

89


Fig. 8. NSM de 1984-2014 en las playas de Els Ferros y La Garrofera, Zona 2

90


Fig. 9. LRR de 1984.2014 en las playas de Els Ferros y La Garrofera. Zona 2

91


Fig. 10. NSM de las playas de L’Alcatí y La Punta, 1984-2014. Zona 3

92


Fig. 11. LRR en las playas de L’Alcatí y La Punta, 1984-2014. Zona 3

93


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


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


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


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


Fig. 3. Distribución de tamaño de las partículas de sedimento en los puntos de muestreo

98


Fig.4. Distribución de tamaño de las partículas de sedimento en los puntos de muestreo

99


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


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


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


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


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


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]

https://www.linkedin.com/vsearch/p?title=Beca+de+colaboraci%C3%B3n&trk=prof-exp-title

https://www.linkedin.com/vsearch/p?company=Departamento+de+Ingenier%C3%ADa+hidr%C3%A1ulica+y+medio+ambiente+%28UPV%29&trk=prof-exp-company-name

https://www.linkedin.com/vsearch/p?company=Departamento+de+Ingenier%C3%ADa+hidr%C3%A1ulica+y+medio+ambiente+%28UPV%29&trk=prof-exp-company-name

https://www.linkedin.com/vsearch/p?title=Pr%C3%A1cticas+en+empresa&trk=prof-exp-title

https://www.linkedin.com/vsearch/p?company=NORAY+ingenier%C3%ADa+y+medio+ambiente&trk=prof-exp-company-name

https://www.linkedin.com/vsearch/p?title=Dise%C3%B1o+y+realizaci%C3%B3n+de+los+talleres+de+la+semana+de+la+ciencia&trk=prof-exp-title

105


>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

analisis de la linea de costa y su relacion con los parametros morfologicos en playas de la safor...

Technology