Libro-Energia-Undimotriz-Cap 4

Centro Mexicano de Innovación

en Energía - Océano

Energía Undimotriz en la Costa Veracruzana, México.

Una Evaluación Geólogo-Geomorfológica de Sitios Idóneos

para el Emplazamiento de Prototipos Ingenieriles

DOI: 10.26359/EPOMEX.CEMIE0220224

4

Morfodinámica Costera:

Tendencias en las Últimas

Cuatro Décadas (1976-2017)

Daniel Morales Méndez, Emilio Saavedra Gallardo, Andrea Mancera Flores,

José Ramón Hernández Santana y Ana Patricia Méndez Linares

Instituto de Geografía, UNAM

Resumen

El presente trabajo ofrece un análisis de la morfodinámica costera veracruzana

en los últimos cuarenta y dos años (1976-2017) a nivel regional para toda la costa

veracruzana, así como en siete sitios particulares con potencial para funcionar

como enclaves de prototipos ingenieriles para el aprovechamiento de energía

undimotriz: Barra de Cazones, Palma Sola, Villa Rica, Punta Roca Partida, Playa

Hermosa, Montepío y Balzapote. La evaluación morfodinámica contempló

la interpretación de materiales aero-satelitales (1973, 1976, 1986, 1995, 2000,

2011, 2017), procesados en ArcGis 10.2, mismo software donde se calculó su

movimiento neto (NSM) y la tasa de punto final (EPR); además de lo anterior, se

157

CEMIE-Océano

la elaboraron e interpretaron perfiles de playa para determinar el estado morfodinámico

de la línea costera.

En el contexto regional, la progradación se observa extendida en la mayor parte

del litoral norte, asociado a condiciones poligenéticas, con sectores aislados de

retroceso en bloques rocosos, tanto sedimentarios de la Llanura Costera del

Golfo Norte (Cazones) o estructuras volcánicas del Cinturón Volcánico Mexicano

(Palma Sola y Villa Rica). Por otra parte, en la costa sur, los elementos acumulativos

de la Llanura Costera del Golfo Sur son interrumpidos por acantilados

volcánicos de la Sierra de los Tuxtlas (donde se localizan los 4 sitios correspondientes

a esta región), aun así, hay un importante componente acumulativo a

manera de bahías interdigitidas entre los promontorios abrasivos.

En cuanto a los sitios particulares, Cazones presenta una tasa de cambio promedio

regresiva de -0.31 m/año, Palma Sola y Villa Rica de +0.58 y +0.47 m/año, respectivamente.

De los cuatro sitios ubicados en el área de los Tuxtlas, las tasas

de cambio promedio en tres (Punta Roca Partida, Playa Hermosa y Montepío)

son menores a ±0.20 m/año, lo cual evidencia cierta estabilidad en los procesos

erosivos y acumulativos, mientras que en Balzapote sobresale el valor de +0.85

m/año, donde la colocación de infraestructura de protección costera favoreció

la acumulación de sedimentos y por ende la progradación costera.

Entre los aciertos se encuentran el cálculo automatizado, la incertidumbre a partir

de los insumos, el enfoque multiescalar y el manejo gráfico que se le puede

dar a los datos, además de proveer argumentos para estrategias de planeación

y manejo integrado del territorio costero. Aun así, por la dependencia a la calidad

de los insumos, la distancia entre transectos, varianza de los resultados y

generalización cartográfica, la evaluación a escala regional tiene un carácter

indicativo y sintético, a diferencia del interpretativo analítico de la escala local.

Palabras clave: morfodinámica costera, progradación y regresión, movimiento

neto de línea de costa, tasa de punto final, Veracruz.

158

Abstract

The present work offers an analysis of the Veracruz coastal morphodynamics

in the last forty-two years (1976-2017) at a regional level for the entire Veracruz

coast, as well as in seven particular sites with the potential to function as

enclaves of engineering prototypes for the use of wave energy: Barra de Cazones,

Palma Sola, Villa Rica, Punta Roca Partida, Playa Hermosa, Montepío and

Balzapote. The morphodynamic evaluation contemplated the interpretation of

aero-satellite materials (1973, 1976, 1986, 1995, 2000, 2011, 2017), processed in

ArcGis 10.2, the same software where its Net Shoreline Movement (NSM) and

the End Point Rate (EPR); in addition to the elaboration and interpretation of

beach profiles and the determination of the morphodynamic state.

In the regional context, progradation is observed extended in most of the north

coast, associated with polygenetic conditions, with isolated sectors of retreat in

rocky blocks, both sedimentary of the North Gulf Coastal Plain (Cazones) or volcanic

structures of the Trans-Mexican Volcanic Belt (Palma Sola and Villa Rica).

On the other hand, on the south coast, the cumulative elements of the South

Gulf Coastal Plain are interrupted by volcanic cliffs of the Sierra de los Tuxtlas

(where the 4 sites corresponding to this region are located), even so, there is an


Una Evaluación Geólogo-Geomorfológica de Sitios Idóneos para el Emplazamiento de Prototipos Ingenieriles

important component cumulative as interdigitized bays between the abrasive

promontories.

Regarding the particular sites, Cazones presents a regressive average rate of

change of -0.31 m/year, Palma sola and Villa Rica of +0.58 and +0.47 m/year,

respectively. Of the four sites located in the Tuxtlas area, the average rates of

change in three (Punta Roca Partida, Playa Hermosa and Montepío) are less

than ±0.20 m/year, which shows certain stability in the erosive and cumulative

processes. while in Balzapote the value of +0.85 m/year stands out, where the

placement of coastal protection infrastructure favored the accumulation of sediments

and therefore coastal progradation.

Among the successes are the automated calculation, the uncertainty from the

inputs, the multi-scalar approach and the graphic management that can be given

to the data, in addition to providing arguments for planning strategies and integrated

management of the coastal territory. Even so, due to the dependence on

the quality of the inputs, the distance between transects, variance of the results

and cartographic generalization, the evaluation at the regional scale has an indicative

and synthetic character, unlike the analytical interpretation of the local

scale.

Keywords: coastal morphodynamics, progradation and regression, net shoreline

movement, end point rate, Veracruz.

Conceptos básicos de morfodinámica

La morfodinámica costera se puede definir como la interacción y la adaptación

continua entre la morfología de la costa y los procesos hidrodinámicos y

eólicos que actúan sobre ella (Silva-Casarín et al., 2014), por lo tanto, se encuentra

determinada por un ciclo de construcción-destrucción, producto de

la acción de diferentes procesos ambientales clasificados en dos componentes:

el externo o forzador, que enmarca a todos aquellos agentes ajenos a

la dinámica continental del relieve, como son las mareas, las olas, el viento,

los organismos biológicos, las actividades antrópicas, etc.; y el interno, que

corresponde a todos aquellos elementos que rigen la morfodinámica terrestre

erosiva y acumulativa, es decir, los elementos condicionados por factores

como la sedimentación; estos procesos también son llamados auto-organizativos

(self-organized processes) (Flaqués et al., 2011; Ruíz-Mártinez et al., 2013;

Jaramillo et al., 2020).

La variación morfológica de la línea de costa es impulsada por variables dependientes

de las condiciones climáticas, las variaciones en el nivel medio del

mar, la energía del oleaje, las tendencias tectónicas, los fenómenos antrópicos,

entre otras (Aiello, et al., 2013; Zagórski et al., 2020). A partir de estos procesos

es posible determinar una serie de parámetros que representen, de una mejor

manera, la línea de costa en su dimensión estática y dinámica. Sin embargo, la

definición convencional que sitúa a la línea de costa como una zona de inter-

159

CEMIE-Océano

faz, presenta problemas en la práctica, si no se toma en cuenta a la dimensión

temporal, que depende de cada caso específico (Kraus y Rosati, 1997; Moore,

2000; Boak y Turner, 2005; Alesheikh et al., 2007; Aiello et al., 2013).

Como ya se mencionó, cada tipo de costa responde de forma diferencial en

dependencia de la acción de las características geológicas, los sedimentos y

los agentes geomorfológicos (Bird, 2008). Conforme a la interacción diferencial

de estos componentes, existen dos clasificaciones, la primera con base en

el cambio relativo en el nivel del mar, y la segunda, sobre la cual se enfoca el

presente estudio, en dependencia del material y su dinámica geofísica (Stanica

y Ungureanu, 2010; Gómez et al., 2014).

Para entender las variaciones en la línea de costa, es necesario además emplear

diferentes líneas de tiempo, con el fin de identificar alteraciones en la

dinámica litoral. En este sentido, se considera la exactitud de las mediciones

y su nivel de incertidumbre, la escala temporal y espacial empleadas, la metodología,

el parámetro seleccionado, la escala de estudio, la resolución de los

insumos, entre otros elementos que son fundamentales para la precisión de

las tasas de cambio resultantes (Dolan et al., 1991; Jaud et al., 2020).

En el presente capítulo, la elección de los insumos y parámetros para la determinación

de la línea de costa se elaboró en función de su disponibilidad, la

resolución de la información, la cantidad de líneas temporales, la continuidad

del parámetro en el tiempo y las características físicas de la zona. En la costa

del estado de Veracruz, debido a su extensión y a las diferentes características

físicas que presenta, se requiere de la elección de múltiples parámetros en

dos escalas de análisis, un estudio con enfoque regional para el análisis de

toda la costa y otro de carácter local para los fines del proyecto del Centro

Mexicano de Innovación en Energía del Océano (cemie-o). De este segundo

enfoque, se desprenden los estudios específicos.

Criterios metodológicos y plataformas tecnológicas

Métodos de análisis de la morfodinámica costera

Estudiar la morfodinámica costera presenta una gran complejidad, debido a

las numerosas formas de abordaje, la amplia variedad de escalas cartográficas

geomorfológicas y el modelo de análisis elegido (Pavlopoulus et al., 2009).

Esta situación está aunada a la constante necesidad de entrelazar entre sí,

condiciones morfogenéticas, morfodinámicas, morfométricas, morfoevolutivas

y morfocronológicas que, a su vez, cuentan con el potencial de ser tratadas en

su conjunto con análisis ambientales de diversas índoles (biológicos, ecológicos,

geológicos, edafológicos, sociales, por mencionar algunas direcciones).

Como presentación de un estado del arte, se dividen a continuación los criterios

metodológicos en métodos cualitativos y semi-cuantitativos, y métodos

meramente cuantitativos. Los métodos conforman diferentes etapas del análisis

morfodinámico.

160



Métodos cualitativos y semi-cuantitativos

Caracterización geomorfológica

De las alternativas de caracterización geomorfológica existentes, a continuación

se mencionan tres tipos principales y útiles para el análisis de la

morfodinámica costera, como se muestra en los estudios de Lalbiakzuali et al.

(2013), Ruíz-Martínez et al. (2013) y Doyle y Woodroffe (2018).

La caracterización morfogenética y morfocronológica, que consiste en la caracterización

del relieve en función del origen, la edad y los procesos involucrados

en la configuración actual de la superficie terrestre (Pedraza, 1996).

La caracterización morfográfica, se refiere a las características geométrico-topológicas

del relieve y su configuración ejemplificada en mapas de unidades

geomorfológicas, donde también se resalta el origen y los procesos actuales

(Serrano-García, 2017) y la morfométrica, que emplea datos de cartas topográficas

o modelos digitales de elevación para obtener modelos topográficos,

pendiente, energía del relieve, rugosidad del terreno, entre otros parámetros

métricos que aportan a la identificación de áreas con mayor susceptibilidad a

determinadas dinámicas morfológicas, como procesos de remoción en masa,

erosión o acumulación eólica, marina o fluvial (Serrano-García, 2017; Hernández-Santana

et al., 2017).

Detección de la línea de costa

La línea de costa es una de las principales características morfológicas de la

zona costera, suele usarse como referencia para los modelos numéricos que

enuncian la variación en su posición durante periodos relativamente cortos

de tiempo, lo cual es producto de la propia naturaleza altamente dinámica del

litoral (Boak y Turner, 2005). Los métodos para su definición espacial son variados,

según la disponibilidad de datos, el parámetro elegido y los objetivos

de la investigación.

Se destacan tres métodos principales para definir la línea de costa:

• El levantamiento de campo: se trazan perfiles de playa, es posible definir

la línea de costa como punto de referencia con la toma de punto gps, ya

sea mediante el establecimiento de la posición de la marea, o la zona

de rompiente. Este método de detección en campo resulta útil para el

monitoreo de procesos a muy corto plazo y el establecimiento de líneas

actuales, a escala local (Navarrete Ramírez, 2014).

• Visualización cartográfica: se refiere a la digitalización manual de la línea

de costa en un sig (posterior a la selección de un parámetro visual) en fotografías

aéreas, ortofotos o cartas topográficas (Boak y Turner, 2005); es

un método ideal para la obtención de líneas de costa antiguas de media

resolución, sin embargo, es necesario considerar el error de digitalización

al analizar tasas de cambio.

• Percepción remota: a través del uso de sensores, imágenes satelitales

y aviones, se definen los parámetros, algoritmos o índices para detectar

161

CEMIE-Océano

la línea de costa, con el uso de imágenes multiespectrales, datos de elevación,

etc. La aplicación de este método permite realizar el monitoreo

constante, a corto plazo, de la línea de costa y abarcar una mayor superficie,

sin embargo, tiene limitaciones en la resolución espacial y, por lo

tanto, en la pérdida de información, tal como lo explica Kuleli (2010).

Perfiles de Playa

La playa corresponde a todo depósito de sedimentos no consolidados, surge

principalmente como producto de la interacción entre la hidrodinámica, la

topografía marina, así como las características y disponibilidad del sedimento

(Cuevas-Jiménez y Euán-Ávila, 2009). Para reconocer patrones y explicar los

procesos de erosión-deposición y el análisis de la morfodinámica costera, se

recurre a la medición de perfiles de playa (Lynch, 2015).

Un perfil es la representación bidimensional de un corte trasversal de la topografía,

de forma general y en el caso de los perfiles de playa, todos muestran

una configuración similar, donde la porción interior del perfil muestra una

mayor inclinación, que disminuye de forma gradual hacia mar adentro, y los

cambios en la pendiente corresponden a las diferentes morfologías presentes

en la playa (Medina et al., 2001).

Se pueden enunciar dos principales metodologías para obtener perfiles de

playa. La primera, empleada desde principios del siglo xx, tal como lo explica

Posada y Buitrago (2009), corresponde al levantamiento topográfico que, en

la actualidad, consiste en un transecto que abarca desde la zona supramareal

hasta el nivel medio de la bajamar o de ser posible hasta la zona donde la altura

del oleaje se modifica por el fondo marino.

Existen distintos métodos para trazar un perfil de playa, en función de las herramientas

con las que se cuente, sin embargo, todos parten desde un punto

fijo en la parte superior o inferior de la playa y el trazo de un transecto de forma

perpendicular a la línea de costa.

• Método de Emery o de Horizonte (1961): Consiste en localizar las coordenadas

iniciales de cada perfil y trazar una línea perpendicular a la costa,

para posteriormente nivelarla al horizonte con ayuda de un distanciómetro,

estadal y clinómetro, a su vez, se establecen marcas de referencia en

cada cambio notable en la inclinación.

• Método de Andrade (2006): Es un método modificado de Emery (1961),

que a través del uso de una manguera de agua hace mediciones secuenciales

de la elevación.

• Método de Abney: Este método, similar al de Emery, emplea un nivel Abney

y un estadal para determinar la distancia y la pendiente de cada porción

de la playa.

La segunda metodología, aplicada por el desarrollo de métodos accesibles

de teledetección de alta resolución en los últimos años, es mediante el uso de

datos de elevación lidar y el trazo de perfiles perpendiculares a la línea de

162



costa, cuya precisión se relaciona directamente con la resolución del raster

(Dean et al., 2013). Sin embargo, su aplicación es limitada en la escala temporal,

dado que los cambios morfológicos son significativos y ocurren en una

temporalidad relativamente corta (Posada y Buitrago, 2009).

En general, los perfiles de playa y la caracterización morfológica y sedimentaria

de la zona costera son herramientas esenciales para su definición

física e hidrodinámica, y puede aplicarse para el monitoreo de fenómenos y

alteraciones a corto plazo, así como en la generación de modelos predictivos

de respuesta del litoral, conforme a la alteración de los patrones morfodinámicos

normales.

Con fines de establecer una línea base previa al establecimiento de los dispositivos

energéticos, se levantaron 56 perfiles, a ambos lados de cada uno de

los siete sitios seleccionados como idóneos (Anexo 1).

Métodos cuantitativos

Índices del estado morfodinámico

El estado morfodinámico es un conjunto de modelos numéricos relacionados

con la evolución a largo plazo, parten del entendimiento y monitoreo de patrones

de alteración en corto plazo, en miras de definir tendencias. Desde mediados

del siglo xx, se ha creado una serie de modelos de evolución morfológica,

como los de Wright y Short (1984), y Masselink y Short (1993), que son una

herramienta útil para evaluar los procesos y morfologías involucrados en la

dinámica costera (Vidal et al., 1995). Entre los más importantes están:

• Parámetro de Dean (W): constituye el primer índice empleado para determinar

el estado modal, refiere que tanto es reflectiva o disipativa una playa

en función de las características del sedimento y del oleaje (McLachlan

y Dorvlo, 2005).

• Índice del estado de la playa (bsI): es el producto de la multiplicación del

parámetro de Dean por el rango de marea (McLachlan y Dorvlo, 2005).

• Índice de depósitos en la playa (bdi): desarrollado por Soares (2003), es

el producto de la pendiente del frente de playa, según el tamaño de grano

del sedimento, útil en zonas de condiciones micromareales (McLachlan y

Dorvlo, 2005; Celentano y Defeo, 2006).

• Área de playa: obtenida mediante la división del rango de marea entre la

pendiente de la playa, expresa el área de interacción mar-tierra y, por lo

tanto, la zona de mecánica sedimentaria (McLachlan y Dorvlo, 2005).

• Índice de playa (bi): similar al anterior, sin embargo, al considerar el tamaño

de la partícula resulta útil para playas con rango mareales uniformes

(Celentano y Defeo, 2006).

• Carrera de marea relativa (cmr): desarrollado por Masselink y Short

(1993), relaciona la velocidad de caída del sedimento con el rango relativo

de la marea y la altura del oleaje, para definir si la playa está controlada

por el oleaje o por la marea (Vidal et al., 1995).

163

CEMIE-Océano

Donde:

H b

/ H sb

: altura significativa del oleaje (en m)

W s

: velocidad de caída del sedimento (en cm/s)

T: periodo de la ola (segundos)

Tide: rango máximo de marea primaveral (en m)

B/Slope: pendiente de la playa

Mz/Sand: tamaño medio de partículas de arena, expresado en unidades phi +1

CM: carrera de marea viva media (en m)

Tasas de cambio de la línea de costa

Se trata del reconocimiento de los patrones de avance o retroceso de la línea

de costa por los procesos de erosión, abrasión y/o deposición, por lo cual

se emplean una serie de métodos estadísticos, que toman bases litorales de

distinta temporalidad para obtener mediciones de su variación espacial y determinar

qué procesos ocurren en la zona costera, en un tiempo determinado

(Dolan et al., 1991).

Plataformas tecnológicas

En la actualidad, los análisis sobre la morfodinámica costera emplean herramientas

sig para expresar, de forma gráfica y bidimensional, todos los procesos

que ocurren en torno a la costa (Mitra, 2011).

Uno de los principales intereses radica en la obtención de datos morfo-evolutivos

y la generación de tasas de cambio, que resulten de utilidad en planes

de gestión costera. Para esto, es necesario contar con herramientas que faciliten

su cálculo y permitan introducir la mayor cantidad de información posible

a diferentes escalas, con el fin de tener un panorama general de la dinámica

litoral, para lo cual se han creado distintos softwares orientados a la modelación

de los procesos costeros, algunos de éstos son:

164



• mike: mediante un sistema de ecuaciones no lineales, genera modelos

bidimensionales de la hidrodinámica costera, así como el

cálculo del transporte sedimentario, los cambios morfológicos en

la zona costera, y permite crear modelos predictivos de la costa a

largo plazo (Tomazin y Re, 2019).

• XBeach: es un modelo numérico de código abierto, que calcula la hidrodinámica

de la costa y la respuesta de ésta durante el embate de una

tormenta (https://www.deltares.nl/en/software/xbeach/).

• genesis: es un sistema de modelación de respuesta costera aplicado

principalmente para reconocer los patrones de comportamiento del oleaje

y del transporte sedimentario (http://www.tvrl.se/hh/genesis.htm).

• dsms/dsas: desarrollado por Thieler et al. (1994), y aplicado posteriormente

por Moore (2000), Morton et al. (2004) y recientemente por Ford

(2013), por mencionar algunos. Consiste en una extensión para ArcGis,

que otorga tasas de cambio de la línea de costa, así como regresiones lineales

simples y ponderadas, a través de análisis multi-temporales; los datos

obtenidos dependen, en gran medida, de la precisión de los insumos.

Este sistema fue empleado en el presente subproyecto para obtener las

tendencias de la morfodinámica costera a lo largo del litoral veracruzano.

Aplicación del Sistema de Análisis Digital

de la Línea Costera (Digital Shoreline Analysis System)

Metodología de aplicación del dsas

dsas es un software que se basa en el programa general integrado de triangulación

analítica (giant), empleado por la noaa y diseñado por Thieler y Danforth

(1994). Conforma una extensión sig que, mediante datos espaciales de la

línea de costa, proporciona estadísticas de cambios en su posición, datos de la

geometría del litoral y la pendiente de la playa, así como modelos predictivos

de la evolución costera, proporciona grados de error y permite la entrada de

valores de incertidumbre (Moore, 2000; Temitope y Oyedotun, 2014).

Para obtener las estadísticas de cambio en la línea de costa en Veracruz, se

utilizó la versión dsas 4.3, compatible con ArcGis10. Para poder emplear esta

herramienta fue necesario realizar un preprocesamiento de los datos.

Preprocesamiento de los datos

Requiere de la creación de una Personal Geodatabase para el almacenamiento

y administración de los datos, así mismo, todos los recursos shape se asignaron

bajo dos sistemas de coordenadas proyectadas: México itrf 92/ utm,

zonas 14 (epsg: 4487) y 15 (epsg: 4488). Además, que se agrega un campo de

fecha, uno de longitud y uno de incertidumbre.

165

CEMIE-Océano

Obtención de líneas de costa

Se realizó la digitalización de forma manual en un sig, conforme a la visibilidad

de la línea de costa, en tres distintos insumos cartográficos y se trazaron las

líneas a partir de los parámetros de Boak y Turner (2005), a una escala media

de 1:7 000.

Los insumos cartográficos utilizados fueron:

• Fotografías aéreas: en México, desde 1967, se cuenta con aerofotos de

vuelo alto (escala 1:50,000, 1:70,000 y 1:75,000) (inegi, 2000). Una de

las principales desventajas de éstas, es que proporcionan una limitada

cobertura temporal. Para el caso veracruzano, se emplearon fotografías

escaneadas y georreferenciadas con estimaciones del error cuadrático

medio (rms) <1m, para crear foto-mosaicos de distinta temporalidad (1973,

1976, 1986 y 1995), en cada uno de los sitios específicos, empleando materiales

de INEGI (1973, 1976, 1986, 1995).

• Ortofotos: Las ortofotos son fotografías aéreas rectificadas geográfica y

geométricamente, es decir, consideran las distorsiones inherentes a la

toma de la fotografía. Permiten obtener mediciones de la línea de costa

con mayor precisión, en dependencia de los datos disponibles y los datos

de elevación, aunque los parámetros de digitalización de la costa pueden

no estar totalmente limitados a la cualidad visual de la línea (Boak y Turner,

2005). Se emplearon ortofotos pancromáticas (inegi, 2000), a escala

1: 20 000 (resolución de 1.5 m) para toda la línea de costa del estado de

Veracruz.

• Imágenes satelitales: las imágenes satelitales son el recurso más reciente

y eficiente en este tipo de estudios, debido a que otorgan una amplia

cobertura, sin embargo, su aplicación está condicionada por la resolución.

Por otra parte, y similar a las ortofotos, la definición de la línea de costa no

se limita a la visualización y puede llevarse a cabo conforme algoritmos

de clasificación y otras técnicas de percepción remota (Jiménez-Orocio,

2010). En Veracruz, se usaron dos líneas obtenidas de imágenes satelitales

RapidEye (2011), con una resolución de 5 m, y Spot 6, del año 2017, a

1.5 m (Spot Image, 2017).

Incertidumbre

Dada la amplia variedad de agentes geomorfológicos presentes en las costas

de Veracruz, en conjunto con los materiales y técnicas empleadas en los estudios

de morfodinámica costera, es común la presencia de inexactitudes que

no se pueden obviar. La incertidumbre es la manifestación cuantitativa de la

fiabilidad de una medición o la precisión de una metodología (Ruiz et al., 2010).

Para determinar la incertidumbre en la medición de los cambios en la línea

de costa de Veracruz, se calculó la incertidumbre estandarizada (ut) (Miranda,

2001) con base en la modificación de dos metodologías (u): Fletcher et al.

166



(2003) y Martin-Prieto et al. (2016), que toman como base la fiabilidad de los

insumos cartográficos y la experiencia del investigador.

Donde:

S: variación estacional de la marea.

RMS: promedio del error cuadrático medio de las fotografías aéreas de la misma

temporalidad.

O: desviación estándar del error de las digitalizaciones.

Ep: tamaño de pixel.

t: variación de la marea.

Los valores de incertidumbre para la costa veracruzana se observan en la

tabla 1.

Como puede observarse, la incertidumbre no supera los ±8 m en el total

de las mediciones. Por otro lado, cada línea de costa tiene su valor de

incertidumbre asociado, que oscila entre los ±8 m para los recursos de menor

tamaño de celda y los ± 32m para recursos de menor resolución, es necesario

mencionar, que mientras se cuente con insumos de mayor resolución y datos

hidrodinámicos más precisos, será posible reducir el grado de incertidumbre

en las mediciones (Fletcher et al., 2003).

Trazo de transectos y análisis estadístico

Una vez determinado el valor de la incertidumbre total y para cada línea de

costa, es necesario trazar una línea base, a partir de la que se generan los

transectos con los que se extraen los estadísticos de tasas de cambio. Esta

baseline consiste en una línea paralela al conjunto de líneas de costa y puede

Tabla 1. Valores de incertidumbre regionales y locales.

Zona Periodo Incertidumbre total estandarizada

UT = ±m

Región norte 2000-2017 7.26

Cazones 1995-2017 7.52

Palma Sola 1995-2017 7.42

Villa Rica 1973-2017 7.82

Región sur 2000-2017 7.26

Punta Roca Partida 1976-2017 7.36

Playa Hermosa 1976-2017 7.30

Montepío 1976-2017 7.35

Balzapote 1976-2017 7.30

167

CEMIE-Océano

trazarse onshore u offshore. En la presente investigación estas líneas base

fueron trazadas a partir de la línea de costa más antigua, de 200 m para los

análisis de carácter regional de toda la costa veracruzana y menor a 50 m para

los siete sitios específicos, considerados desde el punto de vista geólogo-geotécnico

y geomorfológico, como idóneos para el emplazamiento de prototipos

ingenieriles de aprovechamiento de la energía undimotriz. Asimismo, los transectos

se dividieron según la escala de análisis; para obtener las estadísticas

regionales se generaron cada 200 m, mientras que para cada sitio especifico,

se espaciaron a 5 m para obtener mayor precisión de las mediciones.

Estadísticos

Con la aplicación del sistema de análisis dsas es posible obtener índices y

tasas de cambio, basadas en la variabilidad en la posición de la línea de costa

y la asociación de la temporalidad.

• Cambio de la línea de costa envolvente (sce, Shoreline Change Envelope):

este índice se encuentra basado únicamente en la variabilidad espacial

de la línea de costa, sin considerar un periodo de tiempo específico

y define el desplazamiento máximo entre las líneas de costa (Temitope y

Oyedontun, 2014); presenta únicamente una distancia métrica (Martin et

al., 2016) y es útil cuando se pretende generar modelos de pérdidas sedimentarias

máximas y mínimas posibles en la zona costera.

• Movimiento neto de la línea de costa (nsm, Net Shoreline Movement): a

diferencia del índice anterior, este además de considerar la distancia del

desplazamiento de la línea de costa, toma en cuenta la variable temporal,

es decir, se refiere a la distancia total de movimiento entre la línea más

antigua y la más reciente (Martín-Prieto et al., 2016). El cálculo de este

índice permite identificar los sectores históricos de progradación (valores

positivos) y de regresión o retroceso (valores negativos).

• Tasa de punto final (epr, End Point Rate): Esta tasa indica el movimiento

anual de la línea de costa, es producto de la división del movimiento neto

de la línea de costa entre la temporalidad correspondiente (Thieler et al.,

2009). Es una tasa de fácil obtención y requiere únicamente dos líneas de

costa, sin embargo, no puede mostrar tendencias cíclicas al no considerar

la información intermedia (Himmelstoss, 2009).

• Regresión lineal simple (lrr, Linear Regression): este consiste en el ajuste

de la tasa de cambio a una regresión de mínimos cuadrados, es decir,

ajusta una línea producto de la suma de los residuos al cuadrado de todos

los puntos de la costa. La tasa de regresión lineal es la pendiente de la línea,

por lo que es útil al generar modelos predictivos del comportamiento

de la evolución de la línea costera en periodos largos de tiempo. Sin embargo,

el empleo de todos los datos de la tasa de cambio es susceptible a

los efectos de datos atípicos (Himmelstoss, 2009; Temitope y Oyedotun,

2014).

168



• Regresión lineal ponderada (wlr, Weighted Linear Regression): similar al

LRR, pero en este índice, se da mayor peso a los datos “más fiables” o con

menor incertidumbre en la medición, por lo tanto, omite los valores atípicos

de tasa de cambio para trazar la línea de regresión que mejor se ajuste,

el peso en este caso, tal como lo explica Himmelstoss (2009) “se define

como una función de la varianza en la incertidumbre de la medición”,

es decir, considera el factor de incertidumbre como un parámetro para

ponderar las ecuaciones de las regresiones (Martin-Prieto et al., 2016).

Resulta útil para el análisis y predicción de la variación en la posición de la

línea de costa y generar modelos ajustados con menores rangos de error.

Es importante recordar, que las estadísticas presentan aproximaciones sobre

los cambios de la línea de costa, y las tendencias producto de las regresiones

lineales pueden verse influidas por otros factores de la morfodinámica costera,

por lo que la interpretación de los valores siempre está en función del reconocimiento

geomorfológico y dinámico de la zona de interés. Por otra parte, las

fiabilidades de los datos siempre se encuentran condicionadas a la incertidumbre,

el periodo de tiempo analizado y la calidad de los datos (Dolan et al., 1991;

Martin-Prieto et al., 2016).

Antecedentes sobre la morfodinámica

costera en México

En México, la gran mayoría de trabajos sobre morfodinámica costera han sido

realizados durante las últimas dos décadas, con excepción de estudios aislados

en la segunda mitad del siglo xx, todos a partir de distintos enfoques metodológicos,

con el objetivo de delimitar, cuantificar y clasificar la zona costera

nacional. Cronológicamente, se tienen trabajos realizados desde la década

de los 80´s, en el sector Punta Camarones-Cerritos, en las costas mazatlecas

de Sinaloa (Montaño y Peraza, 1986, 1987); a estos estudios siguen los de

Ortiz-Pérez (1988, 1992, 1994), en las costas de Campeche; posteriormente,

Geissert y Dubroeucq (1995) elaboran una caracterización y clasificación regional

del relieve costero veracruzano; Ortiz-Pérez y Benítez (1996) realizan

investigaciones de impacto ambiental en las planicies deltaicas tabasqueñas;

y finalmente, Ortiz-Pérez y Méndez-Linares (1999) evalúan desde una perspectiva

regional, escenarios de vulnerabilidad ante el ascenso del nivel medio del

mar en el golfo de México y el mar Caribe.

A partir del siglo xxi, se amplía parcialmente el panorama regional y local

de las investigaciones sobre la morfodinámica costera. Ortiz-Pérez (2005), y

Ortiz-Pérez y de la Lanza (2006), definen bajo criterios fisiográficos y geomorfológicos,

la delimitación y clasificación de las regiones costeras mexicanas;

Hernández-Santana et al. (2008), analizan la dinámica de la línea costera de

Tabasco, bajo un enfoque estatal y local; Carbajal (2010) y Botello et al. (2010),

aportan nuevos datos sobre la vulnerabilidad costera nacional; asimismo, Ji-

169

170

CEMIE-Océano

ménez-Orocio (2010), investiga la dinámica local de algunas zonas de la costa

veracruzana septentrional; González-Conchas (2011), analiza la estructura,

morfología y morfodinámica de las playas de Mazatlán, entre río Presidio y

Punta Cerritos, Sinaloa; Cruz-González (2012) en las costas de las Marismas

Nacionales, Nayarit; Barrios-Rodríguez (2012), en las costas tabasqueñas; Martínez

et al., (2012), en las costas veracruzanas; Ruíz-Martínez et al. (2013), entre

Cancún y Tulum, Quintana Roo; y Bustamante-Fernández et al. (2016), en la

costa de Punta Diamante hasta la desembocadura del río Papagayo, Guerrero.

Por último, Hernández-Santana et al. (2007, 2016), elaboran la cartografía

morfogenética del relieve costero del estado de Veracruz, bajo un enfoque

jerárquico tipológico.

Tipología de la línea de costa:

criterios de clasificación

La zona costera es el área de convergencia entre la tierra y el mar, frontera

que se define a través de la línea de costa, conforma un sistema altamente

dinámico a través de procesos de erosión y depósito, mismos que constituyen

las variaciones espacio-temporales de la línea (Kwong y Chethika, 2014).

Las variaciones en la posición de este límite permiten cuantificar y conocer

regímenes de progradación o regresión (Oyedotun, 2014). En este sentido, es

indispensable determinar su posición con la mayor precisión espacial posible,

mediante su vectorización con base en recursos cartográficos de distintas

temporalidades disponibles (Dolan et al., 1991). Para determinar la selección

del parámetro de digitalización, en función de la calidad visual de los insumos,

se consideró la clasificación de Boak y Turner (2005), la cual se complementa

con la clasificación del relieve costero para costas secundarias, acorde con

Shepard (1963). Además de esto, se tomó en cuenta la escala de trabajo y la

tendencia de los datos mareográficos, de manera que se seleccionaron siete

parámetros visuales, los cuales se adaptan a la tipología de la costa veracruzana

(figura 1). Se dividió el litoral regional en dos porciones: norte (figura 2) y sur

(figura 3); la digitalización se elaboró conforme a los parámetros mencionados

en el apartado “Obtención de líneas de costa”.

De la longitud total de la costa veracruzana (523.67 km), el 82.88 % corresponde

a zonas interpretadas como acumulativas, el 11.05 % a abrasivas y el

6.07 % a antrópicas. En ambos casos regionales es evidente la extensión de

playas con acantilados dispersos al igual que zonas antrópicas muy particulares.

Índices de playas (Perfiles geomorfológicos)

Para determinar el estado morfodinámico de las playas veracruzanas, como

consecuencia de los escasez de datos de marea y oleaje, a nivel local, se

realizó la comparación de dos índices, el parámetro de Dean (Ω) y el Índice de



Figura 1. Estado morfodinámico de los sitios específicos, conforme a la relación

entre parámetro de Dean e Índice de playa.

Figura 2. Tipología de la costa en la zona septentrional y parámetros de digitalización,

con base en Shepard (1963), y Boak y Turner (2005).

171

CEMIE-Océano

Figura 3. Tipología de costa correspondiente a la zona meridional,

con base en Shepard (1963), y Boak y Turner (2005).

Playa (bi), para resaltar las características físicas de la costa en Veracruz. Ambos

indicadores resultan de gran utilidad para playas micromareales como las

que se presentan en el golfo de México, donde el rango promedio de la marea

es de 0.30 m (López y Pares, 1998). Ambos parámetros son medidas de la capacidad

del oleaje para remover el sedimento, mientras que Dean lo hace a

través de tres elementos: el oleaje, las mareas y el tamaño del sedimento. Por

otra parte, el Índice de Playa considera también la pendiente de la playa como

factor condicionante del oleaje, puesto que, el movimiento del agua sobre el

terreno (swash) interviene con la capacidad del oleaje para remover material

(Defeo y McLachlan, 2005). A partir de esto, el estado morfodinámico se determinó

con base en seis parámetros, obtenidos a partir de 56 perfiles de playa

(Anexo 1), medidos en seis sitios a lo largo de la costa de Veracruz y con una

revisión bibliográfica para obtener los datos de marea, valores de sedimento y

periodo del oleaje (tabla 2).

Con base en el parámetro de la amplitud mareal es posible distinguir tres

principales zonas de estudio: la región norte del estado (Barra de Cazones),

presenta el menor rango de marea, con 0.19 m; la zona centro caracterizada

por los sitios de Palma Sola y Villa Rica, con 0.32 m; y al sur, la zona de los Tuxtlas,

con 0.33m. Sin embargo, a pesar de la variación en los indicadores, en

la figura 4 se aprecia que el estado morfodinamico de las playas en Veracruz

(figura 1) refiere a playas de tendencia disipativa micromareal, con pendiente

172



Tabla 2. Resultados obtenidos de índice de playa y parámetro de Dean.

Localidad Perfil T

tamaño de

grano

unidades

φ*

Barra de

Cazones

Palma Sola

Villa Rica

P

pendiente

%

RMrango

de marea

(m)

BI Hb T Ws Dean

1 3.55 6.99 0.19 -1.02 1.5 8 1.73 10.84

2 3.55 3.93 0.19 -0.76 1.5 8 1.73 10.84

3 3.55 9.91 0.19 -1.17 1.5 8 1.73 10.84

4 3.55 5.24 0.19 -0.89 1.5 8 1.73 10.84

5 3.55 5.24 0.19 -0.89 1 8 1.73 7.23

6 3.55 4.07 0.19 -0.78 1 8 1.73 7.23

7 3.55 6.99 0.19 -1.02 1.5 8 1.73 10.84

8 3.55 6.99 0.19 -1.02 1.5 8 1.73 10.84

9 3.55 18.53 0.19 -1.44 1.5 8 1.73 10.84

10 3.55 6.99 0.19 -1.02 1.5 8 1.73 10.84

11 3.55 6.99 0.19 -1.02 2 8 1.73 14.45

12 3.55 13.17 0.19 -1.29 2 8 1.73 14.45

13 3.55 5.24 0.19 -0.89 1.5 8 1.73 10.84

14 3.55 6.12 0.19 -0.96 1 8 1.73 7.23

15 3.55 4.07 0.19 -0.78 1 8 1.73 7.23

16 3.55 5.24 0.19 -0.89 1 8 1.73 7.23

17 3.55 4.07 0.19 -0.78 1 8 1.73 7.23

18 3.32 8.45 0.32 -0.90 1.3 6 1.73 12.52

19 3.32 14.34 0.32 -1.13 1.3 6 1.73 12.52

20 3.32 8.75 0.32 -0.92 1.3 6 1.73 12.52

21 3.32 8.75 0.32 -0.92 1.3 6 1.73 12.52

22 3.32 14.64 0.32 -1.14 1 6 1.73 9.63

23 3.32 15.32 0.32 -1.16 1 6 1.73 9.63

24 3.32 11.77 0.32 -1.04 0.5 6 1.73 4.82

25 3.32 8.15 0.32 -0.89 0.5 6 1.73 4.82

26 3.32 11.98 0.32 -1.05 0.4 6 1.73 3.85

27 3.32 12.86 0.32 -1.08 0.5 6 1.73 4.82

28 3.32 11.39 0.32 -1.03 0.5 6 1.73 4.82

29 3.32 6.99 0.32 -0.82 0.5 6 1.73 4.82

30 3.32 7.27 0.32 -0.84 0.5 6 1.73 4.82

31 4.00 11.39 0.32 -0.95 0.5 5 0.62 16.13

32 4.00 9.63 0.32 -0.88 0.5 5 0.62 16.13

33 4.00 3.77 0.32 -0.47 1.4 5 0.62 45.16

34 4.00 6.99 0.32 -0.74 1.4 5 0.62 45.16

35 4.00 5.24 0.32 -0.61 1.4 5 0.62 45.16

36 4.00 7.87 0.32 -0.79 1.4 5 0.62 45.16

173

CEMIE-Océano

Tabla 2. Resultados obtenidos de índice de playa y parámetro de Dean.

Localidad Perfil T

tamaño de

grano

unidades

φ*

Villa Rica

Punta Roca

Partida

Playa

Hermosa

Montepio

P

pendiente

%

RMrango

de marea

(m)

BI Hb T Ws Dean

37 4.00 7.87 0.32 -0.79 0.8 5 0.62 25.81

38 4.00 3.49 0.32 -0.44 0.8 5 0.62 25.81

39 4.00 8.75 0.32 -0.83 0.5 5 0.62 16.13

41 4.00 6.12 0.32 -0.68 0.5 5 0.62 16.13

42 3.34 3.18 0.33 -0.46 1 7.7 1.73 7.51

43 3.34 14.11 0.33 -1.11 1 7.7 1.73 7.51

44 3.34 10.78 0.33 -0.99 1 7.7 1.73 7.51

45 3.34 1.40 0.33 -0.10 1 7.7 1.73 7.51

46 3.34 5.52 0.33 -0.70 1 7.7 1.73 7.51

47 3.34 14.95 0.33 -1.13 1 7.7 1.73 7.51

48 3.34 12.50 0.33 -1.05 1 7.7 1.73 7.51

49 3.34 4.37 0.33 -0.60 2 7.7 1.73 15.01

50 3.34 9.91 0.33 -0.95 2 7.7 1.73 15.01

51 3.34 5.68 0.33 -0.71 0.5 8.42 1.73 3.43

52 3.34 11.09 0.33 -1.00 0.5 8.42 1.73 3.43

53 3.34 8.75 0.33 -0.90 2.5 8.42 1.73 17.16

54 3.34 15.84 0.33 -1.16 1.5 8.42 1.73 10.30

55 3.34 10.51 0.33 -0.98 1.5 8.42 1.73 10.30

56 3.34 10.21 0.33 -0.97 1.5 8.42 1.73 10.30

Figura 4. A) Playa Azul en Barra de Cazones, playa disipativa. B) Playa Muñecos, en Palma Sola,

playa intermedia, se puede observar la disipación del oleaje, con más de tres olas en la zona

de barrido. Trabajo de campo mayo-junio de 2018.

174



suave, arena de grano fino <0.25 mm, y más de tres olas en la zona de barrido

(figura 4 A y B).

Acorde con los valores del parámetro de Dean, en la figura 1, es posible distinguir

tres principales tipos de playas: las intermedias (figura 4B), con valores

de 2< Ω <6, que corresponde a Playa Muñecos, en Palma Sola (perfiles 24 al

30) y Playa Montepío (perfiles 51 al 52). En estas playas se presenta la menor

altura del oleaje, mientras que el BI expresa una tendencia a playa disipativa,

manteniendo valores de pendiente bajos, similares a las presentes playas con

parámetros de Dean más altos, correspondientes a playas disipativas (<6 Ω

<20), mismas que concentran más del 70 % de las mediciones. Las características

generales de estas playas reflejan que son anchas, con arena fina y de

escasa pendiente, así mismo, en este rango es posible encontrar los valores

de BI más altos (>-0.50), que corresponden a los perfiles 33, 38 y 40, en Villa

Rica, y 45, en Punta Roca Partida, que son las playas con menor pendiente,

donde el oleaje pierde toda su energía en la zona de barrido, lo cual provoca

menos turbulencia y, por lo tanto, menor capacidad de arrastre del oleaje, en

comparación con los valores de BI más altos >-1.2, concordante con los perfiles

3, 9 y 12, en Barra de Cazones (figura 4A). Estos perfiles exhiben valores

de pendiente superiores a 5° y donde el rango de mareas es el más bajo de

todos, aunque siguen correspondiendo a playas disipativas, donde la turbulencia

presente en la zona de barrido puede ser mayor en comparación a los

sitios donde el índice de playa es menor, facilitando el arrastre sedimentario.

Por último, los valores de Dean Ω <20, corresponden a las playas más disipativas

de toda la zona y pertenecen al sitio de Villa Rica, que tiene el mayor

tamaño del sedimento y el periodo de la ola más corto (5 s), lo que propicia la

presencia de más de cuatro olas en la zona de barrido, evidencia de la predisposición

a la disipación del oleaje y, por lo tanto, a la disminución gradual de su

capacidad de arrastre. En la tabla 2, se muestran los resultados obtenidos para

cada uno de los sitios, además, se anexa toda la información de los perfiles

elaborados en campo (Anexo 1).

Morfodinámica de sitios idóneos para el emplazamiento

de prototipos ingenieriles de conversión energética

La selección de sitios idóneos para el emplazamiento de prototipos ingenieriles

de conversión de energía oceánica a eléctrica se enfocó principalmente

en la energía undimotriz. Para ello, se realizó un reconocimiento general del

litoral del golfo de México, donde a pesar del predominio de zonas netamente

acumulativas, con extensas playas bajas y plataforma continental ancha (Ortiz-Pérez,

2005), se identificaron siete sitios con mediano a alto potencial.

Para la selección de cada sitio, se tomó en cuenta la tipología de la costa,

se consideraron acantilados amplios y mayores a los 10 m de altura, todos

cercanos a pequeñas localidades, con la finalidad de brindar microgeneración

175

CEMIE-Océano

eléctrica durante la etapa experimental y, posteriormente, durante su conexión

a la red eléctrica nacional. En general, la costa de Veracruz se divide en

siete regiones costeras y en tres fueron identificados los siete sitios con alto

potencial geólogo-geotécnico y geomorfológico para el emplazamiento de los

prototipos ingenieriles (tabla 3).

La caracterización morfodinámica, según Masselink et al. (2011), depende directamente

de tres componentes, la geología, la morfología y los forzantes

que corresponden a la acción marina del oleaje, el viento y las actividades

antrópicas, los que fueron considerados para cada uno de los siete sitios.

El primer sitio se localiza al norte del estado, donde la costa se caracteriza

por una extensa llanura de carácter acumulativo fluvial y eólico-marino, playas

arenosas bajas y rectilíneas, donde destaca la presencia de un acantilado sedimentario,

compuesto de caliza margosa y coralina del Oligoceno, al norte de

Barra de Cazones (Cruz-Ortiz, 2007), que forma una terraza tectónico-abrasiva

escalonada, con alturas de 2-3, 5-7 y 10-15 m y, hacia el sur, se extiende una

barra fluvio-marina formada en el estuario del río Cazones (Moreno-Casasola,

2010). Este sitio de Barra de Cazones constituye posiblemente el único sitio

de toda la costa veracruzana, donde aparecen los niveles de terrazas abrasivas

marinas de 2-3 y 5-7 m, del Holoceno tardío y temprano respectivamente

(Hernández-Santana et al., 2007) y, a su vez, junto con el bloque Coatzacoalcos-Agua

Dulce, aledaño al límite meridional del estado, representan las dos

estructuras con ascensos neotectónicos evidentes, en el contexto acumulativo

y de frentes volcánicos de Veracruz.

En la porción central, la llanura costera se reduce debido a la presencia del

macizo volcánico de Palma Sola, lo que da lugar a la formación de salientes rocosos

intercalados con extensas playas abiertas al mar (Geissert-Keintz, 2006).

Otros dos sitios seleccionados como idóneos se corresponden con Palma Sola

y Villa Rica, que son dos promontorios acantilados rocosos, predominantemente

extrusivos y de composición básica, que alcanzan 25 m y 50 m respectivamente.

Tabla 3. Localización de sitios idóneos por regiones costeras veracruzanas.

Región

Sitio potencialmente idóneo

Tuxpan-Nautla

Barra de Cazones

Palma Sola (Miradores)

Laguna Verde- La Mancha

Tómbolo de Villa Rica

Punta Roca Partida

Playa Hermosa

Los Tuxtlas

Montepío

Balzapote

176



Hacia el sector centro-sur, la llanura costera se extiende a manera de extensas

playas y campos de dunas, que se interrumpen únicamente por la Sierra de

los Tuxtlas, un macizo ígneo, donde dominan las costas mixtas abrasivo-acumulativas

(Ortiz-Pérez, 2005). En los frentes lávicos de estas edificaciones volcánicas

se localizan los últimos cuatro sitos, que presentan los acantilados

de mayor altura, compuestos por rocas de origen volcánico extrusivo, que va

desde los más bajos en Balzapote y Playa Hermosa, con 60-80 m, a los más

altos en Punta Roca Partida y Montepío, de hasta 120 m (Hernández-Santana

et al., 2017).

Otro de los parámetros relevantes para la selección de los sitios potenciales

está directamente relacionado con el viento y el oleaje. La costa de Veracruz,

gracias a su orientación, se encuentra todo el año sometida a la incidencia perpendicular

y oblicua de los vientos provenientes del norte, que incrementan su

fuerza y velocidad durante el invierno (Ortiz-Pérez, 2005), esto no solo influye

directamente en el oleaje que mantiene una velocidad y altura constante, que

oscila entre los 0.75 y 1.5 m (cca, 2019), únicamente alterado por la presencia

de fenómenos hidrometeorológicos, cuando la altura del oleaje puede llegar

a ser superior a los 4 m (Ortiz-Pérez y De la Lanza, 2006; Moreno-Casasola,

2010), al igual que su energía.

El oleaje medio es uno de los principales indicadores de la intensidad de

los procesos costeros como la abrasión, la penetración marina e inundación,

mientras mayor sea la altura significante de la ola, mayor es su capacidad de

transporte y abrasión, de tal forma, que la energía del mismo es proporcional

al cuadrado de la altura de la onda (Gornitz et al., 1994; Kumar y Kunte, 2012;

Núñez-Gómez et al., 2017), es decir, que en Veracruz la energía del oleaje promedio

oscila entre 0.57 m y 2.25 m, y llega a ser mayor a 16 m durante alguna

tormenta tropical severa, de modo que, la capacidad máxima del oleaje se

encuentra asociada a la temporada de tormentas.

Este comportamiento no basta para determinar el efecto del oleaje, también

varía en función de la batimetría; cuando el oleaje incidente se acerca a la

línea de costa, la fricción con el fondo marino modifica su altura, orientación y

la fuerza con la que impacta, por lo que mientras mayor sea la pendiente del

fondo marino, menor será la fricción con el mismo, lo cual reduce la disipación

del oleaje, que se traduce en un impacto con mayor fuerza. Esto se puede

observar claramente en los acantilados, donde al no existir una pendiente gradual

que ayude a disminuir la velocidad y altura del oleaje, su efecto abrasivo

es mayor en la región de los Tuxtlas.

Por último, es necesario considerar las mareas, que es la variación periódica

de ascenso y descenso del nivel del mar, provocado por los efectos gravitacionales

de la luna y del sol (Lugo, 2015). Conocer adecuadamente el régimen

mareográfico resulta importante para la selección de los sitios potencialmente

idóneos para el emplazamiento de los prototipos ingenieriles, ya que un rango

mareal muy amplio implicaría un ajuste constante del equipamiento, no solo

177

CEMIE-Océano

por los continuos cambios en el régimen del oleaje y sus incrementos durante

la temporada de tormentas, sino además por el incremento del nivel de la

marea durante la ocurrencia de algún fenómeno hidrometereológico extremo.

En el caso del golfo de México, existe un régimen micromareal predominantemente

diurno, cuyo rango no supera el medio metro (López y Sierra, 1998). Los

datos del Servicio Mareográfico Nacional (smn, 2017) para cada sitio, como se

muestran en la tabla 4.

La morfodinámica costera de cada sitio es el resultado de la constante interacción

de múltiples elementos, procesos, así como fenómenos naturales y

antrópicos; el propio relieve costero es el producto del constante intercambio

de materia y energía, que determina sus rasgos morfológico-morfométricos y

evolución geomorfológica. La aplicación del sistema dsas permitió calcular las

modificaciones de la costa en los últimos cuarenta y dos años (1976-2017) e

identificar los sectores con alto dinamismo de los procesos progradativo y regresivo

de la costa, así como identificar las zonas con una estabilidad notable.

Cambios morfodinámicos de la línea costera

Zona Norte Regional

Movimiento neto de la línea de costa (nsm) 2000-2017

Para la región septentrional veracruzana se generaron 1 190 transectos, con

una longitud de 400 m, perpendiculares a la costa, que se extienden desde

una línea base trazada desde un buffer de distancia de 200 m, con respecto la

línea de costa del año 2000; la incertidumbre total estandarizada es de ±7.26

m. Los parámetros de longitud y distancia de línea base se replicaron en la

región sur.

Durante el periodo de 17 años, comprendidos entre los años 2000 y 2017,

el movimiento neto de la línea de costa muestra valores dispersos, la mayor

progradación (+99.04 m) se localiza en la planicie acumulativa de la Llanura

Tabla 4. Rangos mareográficos por sitios potencialmente idóneos.

Estación mareográfica Sitio potencialmente idóneo Rango (m)

Tuxpan Barra de Cazones 0.065

Veracruz

Palma Sola

Villa Rica

0.522

Alvarado

Punta Roca Partida

Playa Hermosa

Montepío

0.395

Balzapote

178



Costera del Golfo Norte; mientras que el mínimo (-97.17 m) se encuentra en su

contraparte sur. De manera general, predomina la progradación (603 transectos

positivos), en cambio, se detectaron 581 transectos de regresión y seis sin

cambios (figura 5). Los valores acumulativos, entre 0 y +22.11 m, se pueden

observar en ambas llanuras costeras, mientras que los negativos de -17.86 m a

0, corresponden a promontorios aislados del Cinturón Volcánico Mexicano, así

como en la transición litológica entre la provincia antes mencionada y la Sierra

Madre Oriental.

Los valores menores a -17.86 m y mayores a +22.12 m, se localizan en puntos

aislados, asociados principalmente a desembocaduras, campos de dunas y

elementos antrópicos (centros urbanos, puertos, escolleras, etc.). La tendencia

heterogénea de norte a sur se puede verificar en la figura 6, donde la distribución

de los valores de NSM no muestra una aglomeración continua (figura 6A).

Se agrega también la gráfica de la tasa de punto final (epr), que indica -8.96

m/año, como el punto de mayor retroceso, y +5.83 m/año, como el máximo

avance (figura 6B).

Figura 5. Movimiento Neto de la línea de costa (NSM), regional norte. Elaboración propia.

179

CEMIE-Océano

Figura 6. Gráficas lineales Norte. A. Movimiento neto de la línea de costa (NSM);

B. Tasa de punto final (EPR). Elaboración Propia.

En cuanto al cambio en la línea de costa por periodos, el periodo de 2000-

2011 mantiene un régimen progradativo mayor que en el periodo 2011-2017,

debido a una mayor concentración de valores negativos debajo de la media;

se puede notar que el nsm de 2000-2011 tiene un intervalo de valores más

estrecho en los cuartiles centrales (-1 a +2 m sin contar los outlayers), en comparación

con el segundo periodo, donde el intervalo con la mayor cantidad

de datos oscila de -3 a +1 m (figura 7A). Los diagramas de caja que gráfican la

tasa de punto final se muestran similares a los de nsm, de 2000 a 2011, el 50

% de los datos se encuentran en un rango de -1 a 0 m al año, mientras de 2011

a 2017, de -1.5 a +1.8 m (figura 7B).

Para los cálculos, tanto regionales como locales, el módulo dsas genera valores

nulos (null), es decir, líneas que por su orientación y las características

de la costa no calculan de manera correcta los estadísticos; dichos transectos

fueron obviados y discriminados de las evaluaciones. Los resultados fueron

estratificados con el método Natural Breaks (Smith, 2011) para facilitar su visualización

e interpretación.

Cabe agregar, que gráficamente los esquemas de epr y nsm son similares,

debido a la lógica algebraica de los índices en sí mismos (el EPR es un cociente

del nsm dividido entre los años de la información), la mayor diferencia está en

los outlayers y en la agrupación de la caja. En los casos regionales, el epr no

se mapeó por su difusa visualización, al tomar en cuenta que no se aprecia una

tendencia cartográfica, la interpretación a esta escala parte principalmente de

las gráficas estadísticas.

Zona Sur Regional

Al igual que la región norte, visualmente no se aprecia un patrón de distribución

espacial definido. Se delinearon 1 152 transectos, predominando la regresión,

con 641 transectos, 510 de avance y uno sin cambios. El mínimo (-51.05

180



m) coincide con el contacto entre la Llanura Costera del Golfo Sur y la Sierra

de los Tuxtlas; mientras que el máximo (+50.61 m) se emplaza en la curvatura

interior, al sureste (figura 8.).

Coincidente con la región septentrional, en la costa sur de Veracruz los valores

menores a -13.91 m y mayores a +10.42 m, se encuentran asociados a

Figura 7. Diagramas de caja de la zona norte. A. Movimiento neto de la línea de costa (NSM);

B. Tasa de punto final (EPR). Elaboración propia.

Figura 8. Movimiento Neto de la línea de costa (NSM), regional Sur. Elaboración propia.

181

CEMIE-Océano

condiciones morfodinámicas complejas, es decir, donde convergen regímenes

erosivos y acumulativos de varios agentes (fluviales, antrópicos, eólicos, etc.).

Las gráficas de nsm y epr no muestran una tendencia de norte a sur, sobresale

de manera muy general una agrupación central de valores de retroceso, apreciable

en la región asociada a la Sierra de los Tuxtlas (figura 9 A y B).

La dinámica por periodos, a largo plazo en dos lapsos de avance neto, marca

tendencias de retroceso en ambos periodos, de 2000 a 2011, donde el 50 %

de los datos se encuentra entre -10 y +8.5 m, y de 2011 a 2017, entre -2.8 a +1.9

m para epr (figura 10 A y B).

Es evidente la problemática de aplicar el módulo en un contexto regional,

principalmente por la distancia entre transectos, la generalización cartográfica

Figura 9. Gráficas lineales Sur. A. Movimiento neto de la línea de costa (NSM);


182

Figura 10. Diagramas de caja zona sur. A. Movimiento neto de la línea de costa (NSM);




que implica y la ambigua interpretación de los modelos estadísticos. Por lo anterior

y ligado a los objetivos del cemie-o, se realizó el análisis en cada uno de

los siete sitios con potencial para el aprovechamiento undimotriz, los cuales se

desarrollan a continuación.

Barra de Cazones

Movimiento neto de la línea de costa (NSM) 1995-2017

En Barra de Cazones se generaron 971 transectos de 50 m, espaciados a 5

m entre sí, a partir de una línea base a 30 m de distancia, extendida desde el

insumo más antiguo. La incertidumbre total estandarizada fue de ±7.52 m.

Se calcularon 644 transectos con valores negativos, 603 progradativos y 6

sin cambios. Dentro del ámbito predominante de retroceso, el rango menor a

-27.80 m coincide con la desembocadura del río Cazones en su margen norte,

el mínimo (-40.49 m) se localiza justo en el contacto de la porción de influencia

fluvio-marina y el escarpe marino-abrasivo adyacente. Al norte del rasgo abrasivo

mencionado, se localiza una región acumulativa en playa Chaparrales, la

cual es una franja de avance con un máximo de +26.16 m (figura 11).

Resaltan, en el bloque acantilado al sur a la desembocadura del rio Cazones,

valores de retroceso de 0 a -15.13 m, intercalados con algunos transectos de

avance en bloques basculados o con incidencia de procesos de remoción en

masa (figura 12).

Para los casos de sitios específicos, se eligieron los perfiles de playa más

cercanos al desplazamiento máximo y mínimo, en el caso de que exista alguno.

En playa Chaparrales, donde se localiza el perfil más cercano al máximo

(figura 13, perfil A), se puede apreciar una morfología acumulativa homogénea,

mientras que cercano al mínimo (figura 14, perfil B), en primer plano, se observa

la zona donde se trazó el perfil, desde la rompiente de las olas hasta una terraza

fluviomarina erosionada y, en segundo plano, se puede observar el inicio

del promontorio abrasivo.

Tasa de punto final (epr) 1995-2017.

Se dividió la costa de cada sitio específico, según su tipología y los valores de

nsm obtenidos, para facilitar la interpretación por periodos. Para la porción A,

en los tres periodos predomina el retroceso, especialmente de 1995 a 2000,

donde se encuentra el mínimo de ese periodo (-7.24 m/año). Para la porción

B, se ubica de playa Chaparrales como acumulativa, adyacente a Punta Pulpo

hasta el bloque abrasivo; aquí se encuentra el máximo avance del periodo

1995-2000 (+5.5 m/año) y del periodo 2000-2011 (+1.9 m/año) (figura 15).

En el bloque calizo abrasivo (C), donde el retroceso es evidente, se encuentra

el mayor retroceso -4.57 m/año en el periodo 2000-2011; mientras que el

retroceso más importante del periodo 2011-2017, se localiza en la zona D (-7.26

m/año), justamente en la transición entre la playa del sur y el promontorio abrasivo.

Por último, en la zona E, se encuentra la progradación más importante del

periodo 2011-2017 (+5.29 m/año) (figura 15).

183

CEMIE-Océano

Figura 11. Movimiento neto de la línea de costa (NSM): Cazones. Elaboración propia.

Morfodinámica general

De norte a sur, las gráficas lineales de nsm y epr (figura 16 A y B) confirman la

clasificación elaborada para los mapas de cambios anuales. Para ambas gráficas

lo más llamativo es que la zona de retroceso más importante e identificable

se encuentra en la barra adyacente al río Cazones, en su margen norte, a

pesar de configurar una plataforma arenosa, la acción fluvial en combinación

con los agentes costeros y eólicos pueden sugerir el desarrollo primario de

una llanura abrasiva.

Por periodos temporales, en el periodo 1995-2000, el movimiento neto concentra

el 50% de los datos en el intervalo de -7.5 a +6 m de nsm y -1.4 a +1.2

184



Figura 12. Acantilado sedimentario, al sur de la desembocadura del río Cazones. Se observa

el bloque ligeramente basculado y depósitos producto de procesos de remoción en masa.

Trabajo de campo mayo-junio de 2017.

Figura 13, perfil A. Ficha litoral con perfil de playa en Chaparrales, al norte de Punta Pulpo.

Elaboración propia. Trabajo de campo mayo a junio de 2018.

185

CEMIE-Océano

Figura 14, perfil B. Ficha litoral con perfil de Playa Sur, al sur del río Cazones.


m/año de epr; en el periodo 2000-2011, el intervalo con más transectos coincidentes

es de -16 a +4 m de nsm y de -1.8 a +0.2 m/año; por último, el periodo

2011-2017, se ubica entre -10 a +10 m de movimiento neto y de nsm entre -1.8 a

+1.8 m/año. En los dos primeros periodos el régimen predominante fue regresivo,

mientras que de 2011 a 2017, se muestra una distribución normal de los

datos (figura 17 A y B).

Palma Sola


Para el sitio de Palma Sola, el modelo estadístico generó 1 030 transectos con

una longitud de 150 m, perpendiculares a la costa y con una incertidumbre

total estandarizada de ±7.42 m. Las tendencias generales para el nsm, en la

temporalidad 1995-2017 (figura 18) mostraron una progradación de la línea de

costa (con un porcentaje de 56.35 %), especialmente en la porción norte de

playa Palma Sola, donde los valores máximos llegaron hasta los +51.59 m, con

un equivalente de +1.16 m/año. Esta zona se caracteriza por un desarrollo moderado

de dunas, con una altura aproximada de 3 m y un ancho de 25 m (figura

19, perfil A). Los valores más significativos de progradación (+ 27.6 a +51.29) se

localizaron a lo largo de playa Palma Sola, desde el sur de La Loma, hasta aproximadamente

500 m antes de la desembocadura de río Palma Sola. En dicha

186



Figura 15. Tasa de Punto final (EPR), Cazones. Elaboración propia.

Figura 16. Gráficas lineales de Barra de Cazones. A. Movimiento neto de la línea de costa (NSM);


187

CEMIE-Océano

Figura 17. Diagramas de caja de Barra de Cazones. A. Movimiento neto de la línea de costa (NSM);


Figura 18. Movimiento neto de la línea de costa (NSM), Palma Sola. Elaboración propia.

188



Figura 19, perfil A. Ficha litoral con perfil de playa Palma Sola, al sur de La Loma.


playa, se localiza un campo de dunas parabólicas y embrionarias, que llegan a

alcanzar una altura de cuatro metros.

El proceso de regresión más importante se localiza en el promontorio rocoso

La Loma, situado al norte del área, con un valor de -28.29 m, equivalente a

-1.28 m/año. Esta zona se distingue por la existencia de acantilados abrasivos,

en donde se dan procesos gravitacionales importantes, principalmente de

desprendimiento de rocas.

Además de La Loma, existen procesos de regresión significativos en la desembocadura

del río Palma Sola y en el acantilado Miradores. En el río, la acción

fluvial transporta una cantidad importante de sedimentos, dando como resultado

una regresión en el frente de la desembocadura, con valores de -28.29

a -11.47 m; este proceso se extiende al norte, aproximadamente 500 m, donde

alcanza valores de -11.46 a -5.85 m.

Miradores, al igual que La Loma, consiste en una serie de promontorios rocosos

con acantilados de pendientes muy pronunciadas, donde predominan

los procesos gravitacionales. En estos frentes rocosos los valores de regresión

189

CEMIE-Océano

también oscilan entre los -28.29 a los -5.85m. Al sur de Miradores se localiza

Playa Andrea, en donde ocurrieron las menores tasas de regresión -5.84 a 0

m.

Tasa de punto final (epr) 1995-2017

Para la tasa de punto final, se dividió el área en cinco zonas, de acuerdo con

el proceso más significativo con el que coinciden, en su mayoría, al menos dos

temporalidades. En la primera zona (A), se observa que existe un predominio

de regresión, que alcanza los valores máximos entre -4.74 a -1.92 m/año para

el periodo 1995-2000. Esta zona corresponde a una parte del promontorio

rocoso de La Loma; aunque a lo largo de estos acantilados hay porciones de

progradación, sin embargo, éstos no representan la mayoría de los datos (figura

20).

Figura 20. Tasa de punto final (EPR), Palma Sola. Elaboración propia.

190



La zona B es la que mayor discrepancia presenta en los tres periodos. Esta

zona se extiende aproximadamente desde la porción central de La Loma, hasta

el centro de playa Palma Sola. El comportamiento del promontorio rocoso

es similar en las tres temporalidades, sin embargo, discrepa en la playa. Para

el periodo 2000-2011, se ubica el máximo valor de progradación dentro de

este sector, con un valor de +5.51 m/año; por otra parte, el periodo 2011-2017

muestra una regresión importante respecto a los anteriores años.

El sector C, se delimita por el área de influencia del río Palma Sola. Dentro de

esta zona se localizaron tres máximos de valores, tanto de regresión como de

progradación, siendo así el sector con mayor dinámica. Para la temporalidad

2011-2017 existe un marcado contraste entre la porción regresiva (-3.96 m/año)

y la progradativa (+4.02 m/año). Por otro lado, a pesar de que el periodo 2000-

2011 es en su mayoría progradativo, presenta el valor con más regresión (-2.26

m/año). La alta dinámica de este sector se debe a los procesos fluviales, que

son un importante factor de transporte de sedimentos.

En la zona D existe una relativa estabilidad entre los procesos de acumulación

y retroceso en los primeros periodos, no obstante, en la temporalidad

1995-2000, se encuentran los extremos (-4.74 m/año y +3.65 m/año) en el

acantilado Miradores. En los años 2011-2017, se observa que hay un predominio

del proceso acumulativo, a excepción del acantilado Miradores.

La última zona (E), se caracteriza por una dinámica baja, respecto a las demás

zonas. Este sector se localiza desde el extremo sur del acantilado Miradores

hasta la porción sur de playa Andrea, donde para el periodo 1995-2017 ocurre

la máxima progradación costera, en la porción adyacente al acantilado de Miradores,

sin embargo, de forma general esta porción de la costa muestra una

relativa estabilidad con una tasa de cambio promedio de menos de 1 m/año

durante los tres periodos.


En las gráficas lineales se muestran las zonas con mayores datos de progradación

o regresión (de norte a sur). Se observa que, en general, los sectores con

avance o acreción corresponden a playas y campos de dunas, donde hay una

tipología acumulativa del relieve costero; por el contrario, las áreas con mayor

retroceso pertenecen a los promontorios rocosos, los cuáles se encuentran

dentro de la tipología abrasiva. Como se mencionó anteriormente, existe una

mayor dinámica de progradación en Palma Sola, la cual se extiende principalmente

en la playa homónima (figura 21 A y B).

El análisis de los boxplots muestra las principales tendencias en nsm y epr

para cada periodo (figura 22 A y B). Se observa que las temporalidades 1995-

2000 y 2000-2011 del nsm, son claramente progradativas, no obstante, difieren

en los outlayers y la amplitud de los rangos. El primer periodo tiene valores

atípicos de -24 a +18 m, mientras que el segundo carece de ellos, pero cuenta

con una mayor amplitud de valores positivos, además de representar al de

191

CEMIE-Océano

Figura 21. Gráficas lineales, Palma Sola. A. Movimiento neto de la línea de costa (NSM);


Figura 22. Diagramas de caja, Palma Sola. A. Movimiento neto de la línea de costa (NSM);


mayor dinámica entre los tres. Por otro lado, el periodo 2011-2017 muestra una

tendencia negativa, con valores atípicos de ±24 m.

El epr muestra una dinámica muy similar al nsm, ya que los primeros dos

periodos tienen una tendencia positiva y el tercero una negativa. En el periodo

1995-2000, se observa un mayor equilibrio de los datos, con valores atípicos

de -4.8 a +3.8 m/año. La segunda temporalidad también muestra una mayor

amplitud de los valores positivos y una mayor dinámica. En el tercer periodo,

la amplitud de datos negativos y positivos es similar, y los valores atípicos de

esta temporalidad alcanzan los ±4 m/año.

192



Tómbolo de Villa Rica


Para el análisis morfodinámico en el tómbolo de Villa Rica, se tomaron cinco

líneas de costa, en un periodo de 44 años (1973-2017), obteniéndose 677 transectos

de 150 m y cada 5 m, con una incertidumbre estandarizada de ±7.82 m.

La dinámica de la línea costera mantiene un comportamiento similar a la región

central de la costa veracruzana, con dominio de la regresión en el frente acantilado

y rocoso del tómbolo, donde el movimiento oscila entre los -17.03 m y los

-1.32 m, que equivale a un retroceso promedio de -0.38 m/año. La regresión

máxima se presenta en la porción adyacente al promontorio rocoso, donde se

encuentra el contacto entre la costa acumulativa y la abrasiva. En esta zona

el viento incrementa su capacidad de arrastre, como consecuencia del direccionamiento

del flujo, ocasionado por la barrera que forma el acantilado, de

tal forma, que favorece la acumulación de material y la formación de dunas

parabólicas y transversales en dirección ne-sw, sobre la barrera de arena que

forma el tómbolo (figura 23).

Figura 23. Movimiento neto de la línea de costa (NSM) Villa Rica. Elaboración propia.

193

CEMIE-Océano

En las porciones acumulativas de la costa existen dos regímenes diferentes,

el primero corresponde a la zona acumulativa de playa Barra Limón, al norte

del tómbolo, que forma una costa recta abierta al mar, que facilita el arrastre

marino y eólico; se puede reconocer una marcada tendencia erosiva por la

presencia de áreas de deflación y abrasión marina, que exponen una plataforma

de abrasión (figura 24). Aquí la variación en la línea de costa fluctúa

entre los -16.39 m y los +38.32 m, con un retroceso máximo de -0.53 m/año, en

el extremo norte, cerca de la desembocadura del río Limón. La progradación

máxima de +1.01 m/año, se localiza en el contacto con el acantilado, donde se

encuentra la transición morfogenética del tipo de costa.

Por otra parte, al sur del tómbolo, en Playa Villa Rica, la costa adquiere forma

curvilínea y el acantilado forma una barrera para la acción erosiva del viento y

del oleaje. En esta área, el movimiento neto de la línea de costa es de +51.57

m, en promedio, con una progradación máxima de +90.27 m, equivalente a

+2.05 m/año, y una progradación mínima de +0.99 m. En el límite sur, se aprecia

un cambio en la tendencia positiva y encontramos zonas regresivas de

hasta -11.95 m, por la presencia de la boca de la laguna El Llano.

Esta diferencia en la dinámica se aprecia también en la morfología de la costa,

mientras que al norte es más evidente la acción erosivo-acumulativa del

viento, que promueve la formación de campos de dunas parabólicas, transversales

y embrionarias de más de 75 m de ancho y 10 m de elevación (figura 25,

perfil B).

194

Figura 24. Plataforma de abrasión expuesta en la zona acumulativa, al norte

del tómbolo de Villa Rica. Trabajo de campo mayo a junio de 2018.



Figura 25, perfil B. Ficha litoral con perfil de playa Barra Limón, al norte de Punta Las Quebradas.


Al sur, la presencia de la zona habitacional y la laguna El Llano han limitado

el desarrollo dunar, mientras que el ancho de la playa es superior a los 30 m

(figura 26, perfil A); solo es posible encontrar una cadena de dunas de no más

de 3 m, sobre las que se distribuyen construcciones habitacionales de carácter

rural y palapas.


En el análisis de la tendencia de la línea de costa para el tómbolo de Villa Rica,

se consideraron cinco líneas de costa de distinta temporalidad, que constituyen

cuatro periodos y que explican la evolución de la morfodinámica costera,

como se observa en la figura 27. En 44 años se mantiene la tendencia

morfológica de la playa dividida en tres zonas principales: en la primera, es

evidente un cambio en la tendencia durante las últimas cuatro décadas, que

obedece a los procesos dominantes en esta porción de la costa y a la disminución

en los aportes fluviales del río el Limón. De 1973 a 1986, el movimiento

de la línea de costa es relativamente homogéneo, con variaciones que oscilan

hasta los +6.09 m/año. El desplazamiento máximo se aprecia en la zona transicional

con el acantilado semicircular, que corresponde a la zona B, que es el

sector que corresponde a la costa abrasiva. En esta zona de transición ocurren

desplazamientos que van de +5.41 m/año, para el periodo 1973-1986, hasta

-7.03 m/año, durante el periodo de 1986-2000.

195

CEMIE-Océano

Figura 26, perfil A. Ficha litoral con perfil de playa Villa Rica al sur de Punta Las Quebradas.

Elaboración propia. Trabajo de campo mayo-junio 2018.

En contraste, el comportamiento de la línea de costa abrasiva refleja desplazamientos

no mayores a -1.41 m/año, esto por la resistencia del frente rocoso

a la abrasión marina. En comparación, en la porción A se advierte un cambio

brusco en la dinámica para el periodo 1986-2000, con la aparición de una

tendencia regresiva, con retrocesos de hasta -2.2 m/año en la cercanía a la

desembocadura del rio El Limón, que alcanza su máximo retroceso de -3.17 m/

año, durante el periodo 2011-2017.

Para la costa acumulativa, situada al sur del tómbolo, se aprecia que el cambio

a la tendencia progradativa es algo relativamente reciente, solo en el último

periodo de 2011 a 2017, cuando ocurre un desplazamiento máximo de

+10.05 m/año. Sin embargo, en otros periodos anteriores la dinámica costera

presentó una tendencia regresiva, con un máximo de -1.97 m/año, entre 1973-

1986, sobre todo en la boca de la laguna El Llano.

Durante el período comprendido entre los años 1986 y 2000, domina la tendencia

progradativa, con variaciones de hasta +2.33 m/año; mientras que en el

periodo del 2000 al 2011, el retroceso máximo fue de -2 m/año.


Se aprecia el dominio de la progradación, en general, en el sitio del tómbolo

de Villa Rica con dos secciones acumulativas, al norte y al sur del tómbolo, en

196



Figura 27. Tasa de punto final (EPR), tómbolo de Villa Rica. Elaboración propia.

la porción norte encontramos una pequeña zona en regresión menor a -0.5m/

año, que se explica por la presencia de la desembocadura del río Limón, donde

los aportes fluviales favorecen la fluctuación de la línea de costa. Hacia el

sur, la progradación incrementa conforme se acerca al acantilado (hasta +1 m/

año), seguido de un abrupto cambio de la dinámica por la presencia del acantilado

rocoso, donde los valores negativos, indicadores de la regresión, no

llegan a superar los -0.60 m/año, evidencia de la resistencia y la estabilidad

del acantilado. Al sur del tómbolo, en playa Villa Rica, la progradación costera

supera los +2 m/año (figura 28 A y B).

Con la interpretación de los datos agrupados (figura 29), la tendencia general

progradativa de la costa se mantiene para casi todos los periodos, si bien

del 2000 al 2011 existe un ligero cambio hacia la regresión costera, con una

modificación promedio de -0.08 m/año. Este comportamiento heterogéneo de

la costa se explica por tres razones: la primera, por el cierre de la boca de la laguna

El Llano, debido al emplazamiento de un gasoducto, en 1974, que propi-

197

CEMIE-Océano

Figura 28. Gráficas lineales. A. Movimiento neto de la línea de costa (NSM);


Figura 29. Diagramas de caja, tómbolo de Villa Rica. A. Movimiento neto de la línea de costa (NSM);


ció el cese del intercambio sedimentario entre la laguna y la costa; la segunda,

por el incremento de la mancha habitacional; y la tercera, por el frecuente embate

de huracanes, como el huracán Dean, en el año 2007. En contraste, para

el periodo 2011-2017, es evidente el regreso a la tendencia progradativa con un

avance promedio de +0.72 m/año; es también en este año, cuando se presenta

la mayor cantidad de datos atípicos, que pueden explicarse por el cierre de la

boca de la laguna y el depósito de material en las áreas adyacentes a ésta.

Punta Roca Partida


Para el sitio de Punta Roca Partida, se generaron 795 transectos separados

cada 5 m y con una línea base variable, de acuerdo con la tipología, de en-

198



tre 15 y 20 m; la incertidumbre estandarizada fue de ±7.36 m. En el periodo

1976-2017, la tendencia predominante (figura 30) es la progradación (56.35

%), con un valor máximo de +23.35 m, equivalente a +0.57 m/año. Dicha cifra

se localiza en la porción sur de la desembocadura del Arroyo de Liza. Este

valor es el único positivo inmediatamente al sur de la desembocadura, donde

todos los demás son regresivos; esto se debe a la dinámica fluvial, ya que el

arroyo deposita una gran cantidad de sedimentos en un banco de arena una

vez que desemboca al mar. Las zonas con mayor progradación, además de la

mencionada anteriormente, se intercalan en porciones donde existen cadenas

de dunas en zonas acumulativas, o en morfologías cóncavas del promontorio

rocoso (figura 31). Este avance en las porciones abrasivas se puede interpretar

como un proceso de basculamiento; los mayores avances en dichas zonas son

de +2.23 a +23.34 m.

Por otra parte, el mayor retroceso en Punta Roca Partida alcanza los -27.56

m, equivalente a -0.67 m/año. Este valor se localiza en la porción sur del área,

en la cadena de dunas embrionarias (figura 32, perfil A). Las zonas de regre-

Figura 30. Movimiento neto de la línea de costa (NSM), Punta Roca Partida. Elaboración propia.

199

CEMIE-Océano

sión se encuentran dispersas sobre el promontorio y algunas partes de la

playa. Estos retrocesos son consecuencia de procesos eólicos y de remoción

en masa, principalmente desprendimiento de rocas.

Figura 31. Acantilado volcánico de Punta Roca Partida de hasta 65 m de altura y depósitos al pie del

acantilado, producto de los procesos de remoción en masa. Vista desde el norte, desde Playa Ensenada.

Trabajo de campo mayo a junio de 2018.

Figura 32, perfil A. Ficha litoral con perfil de playa en Punta Roca Partida, al sur del acantilado

homónimo. Elaboración propia. Trabajo de campo mayo a junio de 2018.

200




Para la evaluación del EPR se dividió al área de estudio en tres sectores, de

acuerdo con su dinámica (figura 33). La zona A muestra una transición progresiva

hacia el retroceso conforme cambia la temporalidad. Este sector se localiza

sobre el promontorio rocoso, en la porción oeste. El periodo 1976-2000, se

caracterizó por el valor de mayor progradación (+0.51 m/año), aunque en general

este sector es progradativo en esta temporalidad. Para el siguiente periodo

hay un cambio de dinámica, ya que se encuentran porciones de regresión, y

para la temporalidad 2011-2017, la dinámica es predominantemente regresiva.

El sector B está localizado en la porción oriental de los acantilados. El periodo

1976-2000 es totalmente progradativo y para la siguiente temporalidad

se localiza el mayor valor de progradación (+2.82 m/año), sin embargo, la dinámica

cambia, ya que en el último periodo se encuentra el valor de mayor

retroceso en esta sección (-4.40 m/año).

La zona C, debido a que es un sector acumulativo, es altamente dinámico.

Esto se puede observar en los tres periodos, ya que cada uno cuenta con valores

extremos. En el primer periodo se encuentra el de mayor retroceso (-1.56

m/año), en el segundo también está la cifra con más regresión (-2.72 m/año) y,

finalmente, el tercer periodo muestra un cambio importante, debido a que es

totalmente progradativo y su valor máximo es de +6.57 m/año.


Las gráficas lineales muestran las zonas con retroceso o progradación, de este

a oeste. Se advierte que no hay un patrón totalmente definido en cuanto a la

predominancia de un proceso u otro. Esto se debe a que tanto la zona abrasiva,

como la acumulativa, son altamente dinámicas y sus valores de avance o

retroceso se encuentran distribuidos de manera dispersa (figura 34).

Figura 33. Tasa de punto final (EPR), Punta Roca Partida. Elaboración propia.

201

CEMIE-Océano

Los boxplots de Punta Roca Partida muestran que, para los primeros dos

periodos, tanto en nsm, como en epr, hay una ligera tendencia al retroceso.

En el periodo 1976-2000 existen valores atípicos que alcanzan los +12 m, -38

m, en nsm, y +0.5, -1.6 m/año, en epr. La temporalidad 2000-2011 muestra una

distribución similar, tanto en valores negativos como en valores positivos. Por

último, en el periodo 2011-2017 hay un cambio de dinámica hacia la progradación,

presentando una distribución de rangos similares (figura 35).



Figura 35. Diagramas de caja, Punta Roca Partida. A. Movimiento neto de la línea de costa (NSM);


202



Playa Hermosa


Se trazaron transectos separados por 5 m, la longitud de línea base fue 15 m

para costas acantiladas y 20 m para porciones acumulativas desde el insumo

más antiguo (1976), la incertidumbre estandarizada se definió en ±7.3 metros.

El módulo dsas obtuvo 301 transectos, 128 con valores positivos, 19 correspondientes

a valores en 0 y 154 negativos. Ante el predominio de valores de retroceso,

el valor más importante es de -22.43 m, en la bahía ubicada al oeste del

bloque rocoso central. De manera interdigitizada al promontorio mencionado,

se logran apreciar bahías en las que se presenta el máximo de progradación

(+33.83 metros) (figura 36).

Se elaboraron dos perfiles en la bahía localizada al sur del bloque acantilado

(figura 37, perfil A), que muestra una pendiente continua suave, propia de un

entorno acumulativo aislado, con un aporte bajo de material marino derivado

de una rompiente en derrame (figura 37).

Figura 36. Movimiento neto de la línea (NSM), playa Hermosa. Elaboración propia.

203

CEMIE-Océano

Figura 37, perfil A. Ficha litoral con perfil de playa en Playa Hermosa, al sur de Punta Lagarto.


Tasa de punto Final (epr) 1976-2017

La costa de playa Hermosa se dividió en tres sectores: A y C, bahías de amplitud

media con predominio de acumulación litoral; y B, que es un promontorio

abrasivo con bahías intercaladas entre los acantilados. En A y C no se registró

ningún máximo ni mínimo. Para el primer periodo (1976-2000), el mínimo (-1.23

m) se localizó en el contacto del bloque rocoso y A, y el máximo (+1.72 m) en

una de las bahías internas del bloque. En el periodo 2000-2011, el mayor avance

se encontró en la porción acumulativa entre A y B, el mayor retroceso en

el acantilado entre B y C. Curiosamente, para 2011-2017 estas ubicaciones se

invierten, con el máximo retroceso (-3.3 m) entre A y B, mientras que el mayor

avance (+5.59 m) se localiza entre B y C (figura 38).


La gráfica lineal de oeste a este, muestra dos patrones que resaltan, uno de

regresión al oeste, que se ubica en la transición de las regiones A y B, en el

epr, y el central, correspondiente al bloque volcánico, en el que se aprecia la

intercalación de valores en las estructuras disyuntivas, que propiciaron la formación

de bahías (figura 39 A y B).

En el análisis por periodos, se aprecia un predominio regresivo en los dos

primeros lapsos de tiempo, el que cambia a progradativo de 2011 a 2017. En el

204



primer periodo, el intervalo de valores con la mayor cantidad de datos fue de

-6 m a +4 m, en nsm, y -0.20 m a +0.15 m de epr; de 2000 a 2011, -8 m a +4 m

de movimiento neto y -0.80 m a +0.30 m/año; y en el último periodo, de -5 m a

+9 m en todo el periodo y -0.95 m a +1.40 m de epr (figura 40 A y B).

Montepío


La morfodinámica en Montepío, entre 1976 y 2017, mantiene un comportamiento

similar a la zona costera de Los Tuxtlas, en Veracruz, predominando principalmente

la costa abrasiva, que abarca un 64 % de acantilado rocoso, con

Figura 38. Tasa de punto final (EPR), playa Hermosa. Elaboración propia.



205

CEMIE-Océano

un retroceso promedio de -11.17 m y una regresión máxima de -0.86 m/año en

las zonas de debilidad, definidas por fracturas en el promontorio rocoso. Por

otra parte, la sección acumulativa es producto de los aportes fluviales del río

Máquina, donde se presentan los mayores valores de la variación de la costa,

con un avance máximo de +36.68 m y un retroceso de -58.34 m, es decir, una

tasa de cambio de -3.43 m/año en la zona de la desembocadura sobre la que

se encuentra la localidad de Montepío (figura 41).

Figura 40. Diagramas de caja, Playa Hermosa. A. Movimiento neto de la línea de costa (NSM);


Figura 41. Movimiento neto de la línea de costa (NSM), Montepío. Elaboración propia.

206



El alto dinamismo en este sector de la costa, la presencia de una zona habitacional

y el uso turístico de la playa han promovido la modificación de esta zona

con el fin de disminuir el desplazamiento de la costa, como se observa en la

figura 43 A. Los pobladores han realizado trabajos de protección costera con

el objetivo de promover la sedimentación y el incremento de la playa, dado

que en el año 2018 la porción acumulativa era estrecha, de no más de 20 m

de ancho, y completamente antropizada con una morfología convexa (figura

43, perfil B).

Figura 42. Perfil de playa B de Montepío, en la desembocadura del río Máquina,


Figura 43. Área acumulativa de la costa. A. Modificación de la desembocadura del río Máquina.

B. Escarpe de playa al norte de la desembocadura. Fotografías del trabajo

de campo, mayo a junio de 2018.

207

CEMIE-Océano

Esta alteración de los patrones normales de deposición en la desembocadura

ha provocado una disminución en la disponibilidad de sedimentos en porciones

aledañas de la costa, donde si bien para el 2017 mantenían la máxima

progradación, para el 2018 se visualiza un escarpe en la playa, que supera el

metro de alto (figura 43B y figura 44, perfil A), cuyo desarrollo ha sido favorecido

por la acción del oleaje, del viento y de los fenómenos hidrometeorológicos

durante los meses de verano e invierno.

Hacia el norte, la costa acumulativa incrementa su longitud, y la porción de la

playa y las dunas superan los 30 m de ancho. En esta sección se encuentran

dos escarpes como evidencia del arrastre de sedimentos por el oleaje y la

marea. De igual manera, la configuración de la costa mantiene una morfología

rectilínea abierta al mar, lo que ha contribuido al desarrollo de dunas con más

de tres metros de altura, poco comunes en la región de Los Tuxtlas, sobre las

que se ha desarrollado vegetación arbustiva densa, que estabiliza los sedimentos.


La costa de Montepío se puede dividir en tres zonas, acorde con la tendencia

morfodinámica general que manifiestan; al norte y al sur se desarrolla la costa

abrasiva y los acantilados constituyen una barrera que promueve el depósito

Figura 44. Perfil de playa A de Montepío, al norte de la desembocadura del río Máquina,

Elaboración propia. Trabajo de campo, mayo a junio de 2018.

208



de material y el desarrollo de una costa acumulativa, según el sentido de la

corriente de deriva litoral. La primera zona acantilada (A) en la sección septentrional,

está formada por un pequeño bloque rocoso de no más de 150 m de

longitud, que mantiene una tendencia homogénea durante los tres periodos

de tiempo, con un retroceso promedio de -2.99 m, es una zona relativamente

estable en comparación con el área B, que corresponde al sector acumulativo

(figura 45). En esta zona se emplazan la mayoría de los valores extremos,

fundamentalmente cercanos a la desembocadura con valores de -10.14 m/año,

para el periodo de 2011 al 2017, y el comportamiento de la costa es heterogéneo.

Durante el periodo comprendido entre los años 1976 al 2000, existe una tendencia

progradativa con un desplazamiento promedio de +0.26 m/año, y desplazamientos

máximos que no superan los ± 2.5m/año. Este panorama cambia

para el periodo entre 2000 y 2011, cuando la regresión máxima supera los -4

m/año, en la desembocadura del río Máquina, donde en el periodo anterior, se

presentó la máxima progradación. Para el último periodo, del 2011 al 2017, la

tendencia regresiva se mantiene con valores de hasta -10.14 m/año.

La sección C, que corresponde a la mayor extensión de la costa abrasiva,

mantiene un desplazamiento promedio de -0.46 m/año. Cabe destacar, que la

ausencia de datos para el año 2011, sobre todo hacia el extremo sur, dificulta el

análisis de la evolución del acantilado (figura 45).


En las gráficas, lo primero que se puede apreciar son los valores extremos en

la zona de la desembocadura del río Máquina, con cambios abruptos en la

morfodinámica costera, que pasa de una progradación que supera los +25 m,

Figura 45. Tasa de punto final (EPR), Montepío. Elaboración propia.

209

CEMIE-Océano

a un retroceso de cerca de -60 m. En contraste, hacia el sur del sitio, prevalece

un dominio de procesos abrasivos asociados al acantilado rocoso de 120 m

de altura, mientras que la sección de playa no tiene un comportamiento únicamente

acumulativo, sino que se mezclan zonas de acumulación con zonas de

regresión, lo que indica un mayor dinamismo en esta área en comparación con

la zona del acantilado (figura 46 A y B).

Sin embargo, el estado morfodinámico en Montepío mantiene un comportamiento

similar al de otros sitios, como Barra de Cazones, con una morfología

costera de playa disipativa, motivada por el oleaje que alcanza una altura

máxima de 2.5 m. Esta situación podría alterarse al inicio de la temporada de

lluvias y el desarrollo de fenómenos hidrometeorológicos, que incrementan la

velocidad del viento y la intensidad del oleaje.

El análisis de los datos agrupados para cada periodo (figura 47) refleja claramente

la tendencia regresiva en toda la costa, principalmente por la presencia

del acantilado, que abarca la mayor parte del sitio. El promedio del movimiento

de la línea de costa fluctúa entre los -0.034 m/año para el periodo de 1976 y

2000, siendo la menor variabilidad de la costa, con tasas de cambio promedio,

menores a ±1m/año.

Este comportamiento se explica por el amplio rango temporal de 24 años,

en comparación con los periodos subsecuentes de 11 y 6 años, lo que resalta

la importancia no solo de la escala espacial dentro del análisis de la morfodinámica

costera, sino también de la temporal. En contraste, para el periodo

del 2011 al 2017, la tendencia negativa se incrementa, con una tasa de cambio

promedio de -0.57 m/año.

De manera idéntica al análisis de las gráficas lineales, resaltan los datos atípicos

en los gráficos de cajas, que simplemente destacan el papel que juegan

las desembocaduras de los ríos y la interacción entre dos ambientes geomorfológicos

diferentes para el análisis de la variación de la línea de costa.

Figura 46. Gráficas lineales, Montepío. A. Movimiento neto de la línea de costa (NSM);


210



Figura 47. Diagramas de caja Montepío. A) Movimiento neto de la línea de costa (NSM);

B) Tasa de punto final (EPR). Elaboración propia.

Balzapote


Todos los sitios identificados como idóneos para el sur veracruzano se caracterizan

por la formación y consolidación volcánica de su relieve, bajo la influencia

compleja marina, eólica y fluvial, con sectores alternados de playas en

el macizo rocoso, sobre todo en correspondencia con estructuras geológicas

disyuntivas. El sitio Balzapote sobresale por el carácter antrópico de su costa,

donde se construyó un puerto protegido, con escolleras y una pequeña bahía,

a manera de “bolsillo acumulativo”, en el periodo 1976-2000 (figura 48). Este

cambio brusco en la geometría costera determinó el trazado de 323 transectos,

53 transectos regresivos, uno sin cambios y 269 progradativos, delineados

desde una línea base a 20 m hacia la línea más antigua y todos espaciados

a cinco metros.

Los valores progradativos mayores a +57 m, se clasificaron como antrópicos,

debido su distribución en la bahía del puerto artificial, y las escolleras con su

adyacencia abrasiva hacia el este. El máximo valor está centrado en la bahía,

con +136.93 m. En las zonas distales a los elementos construidos, tanto al este

como al oeste, se ubican zonas regresivas, con valores de hasta -14.92 m,

mientras en las proximales de la influencia antrópica, se observan transectos

progradativos con un máximo de +57.29 m (figura 49).


Para evaluar la morfodinámica, se dividió la costa en tres sectores: A, playa occidental

a la escollera; B, zona antropizada; y C, acantilado volcánico oriental.

En el periodo 1976-2000, época de la construcción del puerto, el máximo se lo-

211

CEMIE-Océano

Figura 48. Evolución de la línea de costa del sitio Balzapote, propiciada

por la construcción de escolleras. A. 1976; B. 2000; C. 2011; D. 2017.

212



Figura 49. Movimiento neto de la línea de costa (NSM), Balzapote. Elaboración propia.

calizó en la escollera occidental (+5.24 m/año) y el mínimo en el acantilado del

sector C (-1.31 m/año). Del año 2000 al 2011, la mayor progradación se observa

en la bahía antrópica (+4.15 m/año), y la regresión más notable (-2.74 m/año) en

la región abrasiva C. Para el último periodo, la mayor regresión (-2.57 m/año)

se concentró entre las regiones A y B, adyacente a la escollera occidental, y la

mayor acreción (+5.77 m/año) en uno de los “bolsillos acumulativos” intercalados

en el promontorio abrasivo C (figura 50).


Las gráficas lineales del sitio Balzapote, tanto de EPR como NSM, muestran el

predominio progradativo, con sectores de retroceso muy débiles, correspondientes

a los acantilados rocosos orientales (figura 51).

213

CEMIE-Océano

Figura 50. Tasa de punto final (EPR), Balzapote. Elaboración propia.



En cuanto a los periodos analizados, el periodo 1976-2000, considerado

como de antropización, registró valores progradativos entre +2 y +48 m, con

tasas promedio anual entre +0.1 y +2 m/año. Para los subsecuentes periodos

se presenta una distribución de la mayoría de los datos bajo una tendencia

progradativa, pero conforme a una distribución normal; entre 2000 y 2011, los

valores de nsm se agruparon principalmente entre -6 y +12 m, y de epr entre

-0.8 y +1.2 m/año. El último periodo ubica el 50 % de los datos entre 0 y 10 m

de nsm, así como -0.40 a +1.6 m/año (figura 52).

Resumen estadístico

La costa veracruzana septentrional muestra una tendencia progradativa con

un desplazamiento promedio total de +1.72 m (tabla 5). Analizando cada uno de

los sitios del sector septentrional, Barra de Cazones es el único que mantiene

una tendencia regresiva, con una tasa de cambio promedio de -0.31 m/año,

que se asocia a la convergencia de los procesos fluviales y marinos.

En los otros dos sitios norteños - Palma Sola y tómbolo de Villa Rica -, el comportamiento

y la distribución progradativa coincide con la tendencia general

214



Figura 52. Diagramas de caja zona norte. A. Movimiento neto de la línea de costa (NSM);


Zonas de estudio

Tabla 5. Resumen estadístico general..

Desplazamiento total

promedio (m)

Tasa de cambio

promedio (m/año)

Zona norte regional +1.72 +0.09

Barra de Cazones -6.64 -0.31

Palma Sola +12.76 +0.58

Tómbolo de Villa Rica +20.46 +0.47

Zona sur regional -1.95 -0.10

Punta Roca Partida +0.19 +0.01

Playa Hermosa -0.09 0.00

Montepío -7.57 -0.19

Balzapote +35.03 +0.85

de toda la región septentrional, matizada por los aportes sedimentarios fluviales.

Esta situación contrasta con la región meridional de la costa veracruzana,

donde se distingue una tendencia regresiva, con una variación promedio -1.95

m. Los cuatro sitios del área de Los Tuxtlas muestran tasas de cambio promedio,

menores a ±0.20 m/año (Punta Roca Partida, Playa Hermosa y Montepío),

lo cual evidencia cierta estabilidad en los procesos erosivos y acumulativos,

mientras que en el sitio Balzapote, el emplazamiento de infraestructuras de

protección costera como escolleras, favoreció la acumulación de sedimentos

y, por ende, la progradación costera.

215

CEMIE-Océano

Posibles implicaciones de los cambios

de la línea de costa

En el contexto regional, en la costa situada al norte de la ciudad de Veracruz, la

progradación ocupa la mayor parte del litoral, principalmente asociada a condiciones

geológicas y geomórficas poligenéticas (marinas, fluviales, lacuno-palustres,

antrópicas). A su vez, se proyectan sectores aislados de regresión, en

bloques tectónicos de substrato sedimentario consolidado de la Llanura Costera

del Golfo Norte, como el sitio Barra de Cazones o frentes lávicos de las

estructuras volcánicas de Palma Sola y Villa Rica.

Para la región meridional, extendida hacia el sur de la ciudad de Veracruz,

las llanuras y costas acumulativas de la Llanura Costera del Golfo Sur, se interrumpen

por los acantilados volcánicos de la Sierra de los Tuxtlas, donde

se localizan los cuatro sitios potencialmente idóneos para el emplazamiento

de prototipos ingenieriles de aprovechamiento de la energía undimotriz. No

obstante, se presentan geoformas acumulativas, a manera de bahías interdigitadas,

entre los promontorios abrasivos acantilados.

Si bien la tipología morfogenética general de toda la costa veracruzana coincide

con ciertas tendencias morfodinámicas, es en su interpretación local, donde

se toma un sentido más concreto de las morfologías del relieve costero y

sus procesos formadores primarios y los activos actuales.

En Barra de Cazones, el sector regresivo ubicado al norte de la desembocadura

del río Cazones, es un ejemplo de una plataforma morfogenéticamente

acumulativa, con indicios primarios de abrasión, en un contexto interactivo de

procesos fluviales, eólicos y marinos. Por lo tanto, las zonas en los dos rangos

de nsm, cercanos a cero (-5.13 m a 0, y 0 a +4.74 m) se pueden reconocer como

las porciones más estables e idóneas para la instalación prototipos ingenieriles.

Estas zonas son Punta Pulpo, norte de la playa Chaparrales y la porción

norte del bloque acantilado central.

Palma Sola representa uno de los ejemplos de costa poligenética, donde

los valores de mayor regresión ocurren en el relieve abrasivo, acelerado por

procesos de remoción, y los de mayor progradación en el relieve acumulativo

fluvial y eólico-marino, a excepción de la desembocadura del río homónimo.

Uno de los puntos con mayor embate del oleaje es La Loma, ubicado al norte

de la zona costera de Palma Sola. A pesar de que este promontorio rocoso es

predominantemente regresivo, en la porción central se extiende una geoforma

ligeramente cóncava, donde se encuentran categorías de menor de regresión

(-5.84 m a 0 m), e incluso, de progradación (+0.01 m a +27.59 m). Por lo tanto,

este sitio pudiera considerarse para el enclave de prototipos que aprovechen

la energía undimotriz.

En el tómbolo de Villa Rica existe una relativa estabilidad de la línea de costa,

sobre todo en la zona del frente del acantilado rocoso, con tasas de cambio

inferiores a 1 m/año, en comparación con los dos sectores acumulativos. En las

216



áreas adyacentes al acantilado, se advierte mayor dinamismo, con tasas de

cambio de -4.55 m/año hasta +3.67 m/año, en las cercanías al acantilado y al

frente a la zona habitacional, donde en el último año la costa ha retrocedido

más de 40 m (Arcos, 2020). Esta situación determinó la colocación de una escollera

paralela a la costa, de una longitud aproximada de 100 m (figura 53 B),

con el objetivo de detener el arrastre de sedimentos por el oleaje y las corrientes

litorales (figura 53 A), y disminuir el impacto socioeconómico en la población

de Villa Rica, por la pérdida del espacio, principalmente de uso turístico.

En el sitio Punta Roca Partida existe una dinámica controlada, en gran parte,

por el promontorio volcánico, ubicado en toda su porción septentrional y

central. La morfodinámica tiende a un equilibrio relativo entre la progradación

(56.35%) y la regresión (43.65%), debido a que los procesos de remoción son

compensados por los de basculamiento. Dada la energía de oleaje en este

macizo volcánico, Punta Roca Partida representa uno de los sitios con mayor

potencial para el aprovechamiento de la energía undimotriz.

El incremento de las tasas anuales de cambio en playa Hermosa, durante el

período 2011-2017, si bien no exceden -1 m/año, puede corresponder a eventos

hidrometeorológicos extremos, entre ellos seis huracanes de categorías 1 y 2

de la escala Saffir-Simpson, y 15 tormentas tropicales. Aun así, es uno de los

sitios más estables, tanto en su sector acantilado como en las bahías acumulativas

interiores.

Los acantilados más desarrollados y vigorosos de la costa veracruzana se

localizan en el sitio de Montepío. Esta costa abrasiva acantilada mantiene tasas

promedio de cambio por debajo de 1 m/año, tendencia que se mantuvo en los

distintos periodos de análisis. El cerro Borrego presentó un desplazamiento

mínimo de la costa, con una variación promedio de -0.46 m/año para los tres

periodos de análisis. Por su estabilidad y accesibilidad este sitio es uno de

los idóneos para la instalación de los prototipos ingenieriles en la costa. En

contraste con esta estabilidad, en la zona de la desembocadura se presenta el

Figura 53. Cambios antrópicos en el litoral del tómbolo de Villa Rica: A. Imagen obtenida

por un vehículo aéreo no tripulado; B. Vista en planta (Fuente: CIPAES MX, 2021).

217

CEMIE-Océano

mayor dinamismo, producto de la interacción fluvio-marina, donde la tasa de

cambio puede superar los -10 m/año.

Por último, en el sitio Balzapote la acción antrópica determinó las tendencias

progradativas máximas entre las escolleras que delimitan el puerto, aunado al

transporte sedimentario fluvial local. El impacto del componente humano y la

tendencia irregular de los valores de cambio en los sectores acumulativos y

en el bloque abrasivo, adyacentes a la bahía, vuelven evidente la necesidad

de monitorear con estudios geomorfológicos multitemporales, que aporten a

la comprensión de sistemas complejos con variaciones en cortos lapsos de

tiempo.

Conclusiones

La morfodinámica costera, mediante el cálculo de las tasas regresivas y progradativas,

ofrece una mirada multitemporal de los sistemas litorales, a partir

del cambio en la línea de costa. Entre los aciertos del módulo dsas se encuentra

el cálculo automatizado, la consideración de la incertidumbre a partir de

los insumos, la representación espacial en entornos locales y el fácil manejo

gráfico de los datos, además de proveer argumentos cuantitativos para las

estrategias de planeación y manejo integrado del territorio costero. Las dos

principales limitantes son claramente la dependencia a la calidad de los insumos

y al manejo regional, más que la distancia entre los transectos, la varianza

de los resultados y la ambigua cartografía. Esto determina la necesidad de los

estudios locales, funcionando como una guía, más que como un instrumento

certero. La evaluación a escala regional tiene un carácter indicativo y sintético,

sobre el interpretativo y analítico de la escala local.

El objetivo de generar una línea base para efectos de estudios litorales está

satisfecho, pues además de ofrecer un instrumento metodológico para evaluar

la morfodinámica costera, estos indicadores en su conjunto con interpretaciones

globales del sistema terrestre pueden ser aplicados en aras de la instalación

de prototipos ingenieriles de aprovechamiento energético undimotriz.

Por otro lado, es importante relacionar estos aspectos físicos con las

dinámicas de espacios socialmente construidos, con miras a la generación de

políticas públicas municipales de preservación y protección de ecosistemas

costeros, de planes de desarrollo, de manejo integrado del paisaje, de intervenciones

locales, de iniciativa participativa y demás proyectos proyectos interseccionales

e intersectoriales que permitan a las comunidades costeras,

tanto en el golfo de México como en todo el litoral nacional, interactuar y construir

su medio ambiente en vías de mejorar las realidades que cada contexto

específico involucra.

218



Referencias

Aiello, A., Canora, F., Pasquariello, G., Spilotro, G., 2013. Shoreline Variations an coastal dynamics:

A space-time data analysis of the Jonian litoral, Italy. Esturaine, Coastal and Shelf

Science, 129: 124-135. DOI:10.1016/j.ecss.2013.06.012.

Alesheikh, A., Ghorbanali, A., Nouri, N., 2007. Coastline change detection using remote

sensing. International Journal of Environmental Science and Technology, 4(1): 61-66.

DOI:10.1007/BF03325962.

Andrade, F. y Ferreira, M. A., 2006. A simple method of measuring beach profiles. Journal of

Coastal Research, 22(4): 995-999. DOI: 10.2112/04-0387.1.

Arcos, A., 2020. Ve en qué lugar de Veracruz se quedan sin playa y sin trabajo. El Sol de

Córdoba. Diario de Xalapa. 2020/09/20. Disponible en línea https://www.elsoldecordoba.

com.mx/local/ve-en-que-lugar-de-veracruz-se-quedan-sin-playa-y-sin-trabajo-pandemiamedio-ambiente-desempleo-villa-rica-5784807.html

consultado el 3 de marzo del 2021.

Barrios-Rodríguez, M. L., 2012. Variaciones morfodinámicas de la costa deltaica de Tabasco.

Tesis de Licenciatura. Universidad Nacional Autónoma de México, Ciudad de México.

Bird, E., 2008. Coastal Geomorphology: an introduction. University of Melbourne, Australia,

Ed. Wiley, 2nd. edition, 436 p.

Boak, E. y Turner, I., 2005. Shoreline definition and detection: A review. Journal of Coastal

Research, 21(4):688-703. DOI:10.2112/03-0071.1.

Botello, A., Villanueva, S., Gutiérrez, J., Rojas, J. L. (eds.), 2010. Vulnerabilidad de las zonas

costeras mexicanas ante el cambio climático. SEMARNAT-INE, UNAM-ICMYL, Universidad

Autónoma de Campeche, 514 pp.

Bustamante-Fernández, J. E., Méndez-Linares, A. P., Hernández-Santana, J. R., 2016. Variaciones

morfodinámicas hiperanuales (1951-2007) de la línea costera entre Punta Diamante y

Río Papagayo, Acapulco, México. Boletín de la asociación de Geógrafos españoles, 71:301-

322. DOI:10.21138/bage.2284.

Carbajal-Domínguez, J. A., 2010. Zonas costeras bajas en el Golfo de México ante el incremento

del nivel del mar, pp. 359-380. En Alfonso Vázquez Botello, Susana Villanueva-Fragoso,

Jorge Gutiérrez, y José Luis Rojas Galaviz (eds.) Vulnerabilidad de las zonas costeras mexicanas

ante el cambio climático. SEMARNAT-INE, UNAM-ICMYL, Universidad Autónoma de

Campeche, 514 pp.

CCA, Centro de Ciencias de la Atmósfera, 2019. Altura del oleaje en el Golfo de México. Universidad

Nacional Autónoma de México. Disponible en http://grupoioa.atmosfera.unam.

mx/pronosticos/index.php/oleaje/gom/gom-alt

Celentano, E. y Defeo, O., 2006. Habitat harshness and morphodynamics: life history traits of

the mole crab Emerita brasiliensis in Uruguayan sandy beaches. Marine Biology, 149(6),

1453-1461. DOI:10.1007/s00227-006-0309-1.

Centro Integral de Planeación y Análisis Espacial (CIPAES MX), 2021. Fotografías de campo

en planta y obtenidas por un VANT, Villa Rica, Veracruz. Disponible en https://www.facebook.com/CipaesMex.

Cruz-González, M., 2012. Variaciones estacionales en la línea de costa entre el Canal de Cuautla

y el estero San Cristóbal, Nayarit. Tesis de Licenciatura. Universidad Nacional Autónoma

de México, Ciudad de México.

Cruz-Ortiz, N, L., 2007. Petrografía y geoquímica de arenas de playa de México: Implicaciones

de procedencia y ambientes tectónicos. Tesis de Ingeniería. Instituto de Ciencias Básicas

e Ingeniería, Centro de investigaciones de Ciencia de la Tierra, Universidad Autónoma del

Estado de Hidalgo, 144 pp.

219

CEMIE-Océano

Cuevas-Jiménez, A. y Euán-Ávila, J., 2009. Morphodynamics of carbonate beaches in the Yucatán

Peninsula. Ciencias Marinas, 35(3), 307-320. DOI: 10.1007/s00227-006-0309-1.

Dean, R., Walton, T., Rosati, J. y Absalonsen, L., 2013. Beach Erosion: Causes and Stabilization.

Coastal Hazards, 6:319-365. DOI: 10.1007/978-94-007-5234-4_13.

Defeo, O. y McLachlan, A., 2005. Patterns, processes and regulatory mechanisms in sandy

beach macrofauna: a multi-scale analysis. Marine Ecology Progress Series, 295:1-20.

DOI:10.3354/meps295001.

Dolan, R., Fenster, M.S., Holme, S. J., 1991. Temporal analysis of shoreline recession and accretion.

Journal of Coastal Research, 7(3):723-744. https://www.jstor.org/stable/4297888.

Doyle, T. B. y Woodroffe, C. D., 2018. The application of LiDAR to investigate foredune

morphology and vegetation. Geomorphology, 303:106-121. DOI: 10.1016/j.geomorph.2017.11.005.

Emery, K. O., 1961. A simple method of measuring beach profiles. Limnology and Oceanography,

6(1):90-93. DOI:10.4319/lo.1961.6.1.0090.

Flaqués, A., Calvete D., Van den Berg N., Ribas F., Fernández A. y Caballeria M., 2011. Self

organized morphological patterns in coastal dynamics. Recerca i innovació a l’escola de

camins:1-4.

Fletcher, C., Rooney, J., Berbee, M., Lim, S., Richmond, B., 2003. Mapping shoreline change

using digital orthophotogrammetry on Maui, Hawaii. Journal of Coastal Research, 38: 106-

124. https://www.jstor.org/stable/25736602.

Ford, M., 2013. Shoreline changes interpreted from multi-temporal aerial photographs and

high resolution satellite images:Wotje Atoll, Marshall Islands. Remote Sensing of Environment,

135:130-140. DOI: 10.1016/j.rse.2013.03.027.

Geissert, D. y Dubroeucq, D.,1995. Influencia de la geomorfología en la evolución de los suelos

de dunas costeras en Veracruz, México. Investigaciones Geográficas, 3:37-51.

Geissert-Keintz, D., 2006. El entorno Físico: Procesos y Cambios. En Patricia Moreno-Casasola

(ed.) Entornos Veracruzanos: La costa de la Mancha. Xalapa de Enríquez: INECOL, 576 pp.

Gómez, C., Alberto, J., Arce, G., 2014. Morfología litoral: Agentes, procesos y formas resultantes.

Resúmenes. Facultad de Humanidades. UNNE Revista geográfica digital, 11(21):53.

DOI:10.30972/geo.11212118.

González-Conchas, M. A., 2011. Morfodinámica de la línea costera comprendida entre la

desembocadura del Río Presidio y Punta Cerritos, Mazatlán, Sinaloa, México, durante la

segunda mitad del siglo XX y el presente. Tesis de Maestría, Posgrado en Geografía de la

UNAM, Ciudad de México.

Gornitz, V. M., Daniels, R. C., White, T. W. y Birdwell, K. R., 1994. The development of a coastal

risk assessment database: vulnerability to sea-level rise in the US Southeast. Journal of

Coastal Research, Special Issue 12:327-338.

Hernández Santana, J. R., A. P., Méndez Linares, M., Figueroa Mah-Eng, 2007. Caracterización

del relieve plio-cuaternario del entorno costero del Estado de Veracruz, México. Cuaternario

y Geomorfología, 21(3-4):113-131. http://tierra.rediris.es/CuaternarioyGeomorfologia/images/vol21_3_4/21Hernandez.pdf

Hernández-Santana, J. R., Ortiz-Pérez, M. A., Méndez-Linares, A. P. y Gama-Campillo, L. 2008.

Morfodinámica de la línea de costa del estado de Tabasco, México: Tendencias desde la

segunda mitad del siglo XX hasta el presente. Investigaciones Geográficas, 65:7-21.

Hernández-Santana, J. R., Méndez-Linares, A. P., López-Portillo, J. y Preciado-López, J. C.

2016. Coastal geomorphological cartography of Veracruz State, México. Journal of Maps,

12(2):316-323. DOI:10.1080/17445647.2015.1016128.

220



Hernández-Santana, J. R., J. L., Pérez-Damián, F., Rosete-Vérges, M., Villalobos-Delgado, A. P.,

Méndez-Linares, E., Navarro-Salas, 2017. Clasificación geomorfométrica del relieve mexicano:

una aproximación morfográfica por densidad de curvas de nivel y la energía del

relieve. Investigaciones Geográficas, 94:1-15. DOI: 10.14350/rig.57019.

Himmelstoss, E. A., 2009. DSAS 4.0 installation instructions and user guide. Thieler, ER, Himmelstoss,

EA, Zichichi, JL, and Ergul, Ayhan, 2008-1278.

Instituto Nacional de Estadística, Geografía e Informática (INEGI), 1973. Fotografías aéreas, a

escala 1:70000. México, D. F.

Instituto Nacional de Estadística, Geografía e Informática (INEGI), 1976. Fotografías aéreas a

escalas 1:50,000 y 1:70,000. México, D.F.

Instituto Nacional de Estadística, Geografía e Informática (INEGI),1986. Fotografías aéreas, a

escala 1:75000. México, D. F.

Instituto Nacional de Estadística, Geografía e Informática (INEGI),1995. Fotografías aéreas, a

escala 1:75000. Aguascalientes, México.

Instituto Nacional de Estadística, Geografía e Informática (INEGI), 2000. Ortofotos digitales, a

escala 1:20 000. Aguascalientes, México.

Jaud, M., Bertin, S., Beauverger, M., Augereau, E., & Delacourt, C., 2020. RTK GNSS-Assisted

Terrestrial SfM Photogrammetry without GCP: Application to Coastal Morphodynamics

Monitoring. Remote Sensing, 12(11): 1889. DOI:10.3390/rs12111889.

Jaramillo, C., González, M., Medina, R. y Turki, I. 2020. An equilibrium-based shoreline rotation

model. Coastal Engineering, 163:17 p. DOI: 10.1016/j.coastaleng.2020.103789.

Jiménez Orocio, O. A., 2010. Análisis espacial-temporal de la erosión y acreción de las costas

del Estado de Veracruz. Tesis de maestría en ciencias. Instituto de Ecología A. C., Xalapa,

México.

Kraus, N. C. y Rosati, J. D., 1997. Interpretation of shoreline-position data for coastal engineering

analysis. Coastal Engineering Technical Note CETN-II-39, U.S. Army Engineer Waterways

Experiment Station, Vicksburg, 13 pp. Retrieved from https://erdc-library.erdc.

dren.mil/jspui/bitstream/11681/2118/1/CETN-II-39.pdf.

Kuleli, T. 2010. Quantitative analysis of shoreline changes at the Mediterranean Coast in

Turkey. Environmental Monitoring Assessment, 167, 387-397. DOI:10.1007/s10661-009-

1057-8.

Kumar, A. A. y Kunte, P. D., 2012. Coastal vulnerability assessment for Chennai, east coast of

India using geospatial techniques. Natural Hazards, 64(1), 853-872. DOI:10.1007/s11069-

012-0276-4.

Kwong, F.A. y Chethika, W.D.G., 2014. Impact of Coastal Land Use Change on Shoreline

Dynamics in Yunlin County, Taiwan. Environments, 1:124-136. DOI:10.3390/environments1020124.

Lalbiakzuali, C., Suresh-Gandhi, M. y Ram-Mohan, V., 2013. Geomorphology and evolution

of coastal landforms between Mimisal and Rajamadam, Palk Strait, east coast of India.

Current Science, 105(12):1730-1735. https://www.jstor.org/stable/i24099524.

López, B. M. y Pares-Sierra, A., 1998. Circulación del Golfo de México inducida por mareas,

viento y la corriente de Yucatán. Ciencias Marinas, 24(1):65-93. DOI:10.7773/cm.

v24i1.740.

Lugo-Hubp, J., 2015. Diccionario geomorfológico. Geografía para el siglo XXI, UNAM, Instituto

de Geografía, Textos Universitarios, 7:479.

221

CEMIE-Océano

Lynch, M. G. A., 2015. Propuesta de un índice de vulnerabilidad costera aplicado a tres sectores

de la Región Lima (Doctoral dissertation, Pontificia Universidad Católica del Perú,

Facultad de Letras y Ciencias Humanas. Mención: Geografía y Medio Ambiente.

Martín-Prieto, J. A., Roig-Munar, F. X., Rodríguez-Perea, A., Mir-Gual, M., Pons-Buades, G. X.,

Gelabert-Ferrer, B., 2016. La erosión histórica de la playa de Sa Rápita (S. Mallorca). Investigaciones

Geográficas, 66:135-154. DOI:10.14198/ INGEO2016.66.08.

Martínez, M. L., Vázquez, G., López-Portillo, J., Psuty, N. P., García-Franzo, J. G., Silveria, T. M.,

Rodríguez-Revelo, N. A., 2012. Dinámica de un paisaje complejo de la costa de Veracruz.

Investigación Integral, 4(1):151-160.

Masselink, G. y Short, A. D., 1993. The effect of tide range on beach morphodynamics and

morphology: a conceptual beach model. Journal of Coastal Research, 9(3):785-800.

Masselink, G., Aagaard, T. y Kroon, A., 2011. Destruction of intertidal bar morphology during

a summer storm surge event: Example of positive morphodynamic feedback. Journal of

Coastal Research, Special Issue 64:105-109.

McLachlan, A. y Dorvlo, A., 2005. Global patterns in sandy beach macrobenthic communities.

Journal of Coastal Research, 214:674-687.

Medina, R., Bernabeu, A. M., Vidal, C. y Gonzalez, M., 2001. Relationship between beach

morphodynamics and equilibrium profiles. In Coastal Engineering 2000, pp. 2589-2601.

Miranda, J., 2001. Evaluación de la incertidumbre en datos experimentales. Instituto de Física,

UNAM, Ciudad de México, 42 pp.

Mitra, D., 2011. Remote Sensing and GIS for Coastal Zone Management: Indian Experience.

Geoinformatics in Applied Geomorphology, 63.

Montaño-Ley, Y. y Peraza-Vizcarra, R., 1986. Aplicación de un modelo de tipo predictivo para

cuantificar el transporte litoral de arenas en dos playas de la costa Sur de Sinaloa. Anales

del Instituto de Ciencias del Mar y Limnología, UNAM. Disponible en http://biblioweb.tic.

unam.mx/cienciasdelmar/instituto/1986-1/articulo206.html, consultado el 7 de enero

de 2017.

Montaño-Ley, Y. y Gutiérrez-Estrada, M., 1987. Control de perfiles de playa en el área de

Mazatlán, Sinaloa, México. Anales del Instituto de Ciencias del Mar y Limnología, UNAM,

14(1):1–11.

Moreno-Casasola, P., 2010. Veracruz, Mar de arena. VII Construcción y destrucción de un

paisaje costero. Xalapa de Enríquez: INECOL, 216-259.

Moore, L., 2000. Shoreline mapping techniques. Journal of Coastal Research, 16:111-124.

http://www.jstor.org/stable/4300016.

Morton, R. A., Miller, T. L., y Moore, L. J., 2004. National assessment of shoreline change:

Part 1: Historical shoreline changes and associated coastal land loss along the U.S. Gulf of

Mexico: U.S. Geological Survey Open-file Report 2004-1043, 45 p.

Navarrete-Ramírez, S. M., 2014. Protocolo Indicador Variación línea de costa: perfiles de

playa. Indicadores de monitoreo biológico del Subsistema de Áreas Marinas Protegidas

(SAMP). Invemar, GEF y PNUD. Serie de Publicaciones Generales del Invemar, 73:1-36.

Núñez-Gómez, J, C., Ramos-Reyes R., Barba-Macías E., Espinoza-Tenorio A. y Gama-Campillo

L. M., 2017. Índice de vulnerabilidad costera del litoral tabasqueño, México. Investigaciones

Geográficas, 91:70-85. DOI:10.14350/rig.50172.

Ortiz-Pérez, M. A., 1988. Evidencias de cambios geomorfológlcos del sistema litoral mediante

el análisis de imágenes aéreas, Memorias, Ecología y Conservación del delta de los ríos

Usumacinta y Grijalva, Instituto Nacional de Investigaciones sobre Recursos Bióticos, División

Regional Tabasco, Gobierno del Estado de Tabasco, pp. 43-54.

222



Ortiz-Pérez, M. A., 1992, Retroceso reciente de la línea de costa del frente deltaico del río San

Pedro, Campeche, Tabasco. Investigaciones Geográficas, 25:7-23.

Ortiz-Pérez, M. A., 1994. Repercusiones del ascenso del nivel del mar en costas bajas de planicies

deltaicas. Geografía y Desarrollo, 2(11):69-72.

Ortiz-Pérez, M. A. y Benítez, J., 1996. Elementos teóricos para el entendimiento de los problemas

de impacto ambiental en las planicies delticas: la región de Tabasco y Campeche, pp.

483-503. En Alfonso Botello, José Luis Rojas Galaviz, Jorge Benítez y David Zárate Lomelí

(eds.) Golfo de México, contaminación e impacto ambiental: diagnóstico y tendencias.

EPOMEX, Serie científica 5, Universidad Autónoma de Campeche.

Ortiz-Pérez, M. A., Méndez-Linares, A. P., 1999. Escenarios de vulnerabilidad por ascenso del

nivel del mar en la costa mexicana del Golfo de México y el Mar Caribe. Investigaciones

Geográficas, 39:68-81.

Ortiz-Pérez, M. A., 2005. Características físicas de las costas: base para su regionalización. El

caso de la costa veracruzana, pp. 81-95. En Patricia Moreno-Casasola, Elisa Peresbarbosa

Rojas y Ana Cecilia Travieso Bello (eds.) Manejo Costero Integral: el enfoque municipal.

Instituto de Ecología, Xalapa, 1266 pp.

Ortiz -Pérez, M. A. y De La Lanza-Espino, G., 2006. Diferenciación del espacio costero de México:

Un inventario regional. Geografía para el siglo XXI. Instituto de Geografía, UNAM, Serie

de Textos Universitarios, 3, 143 p.

Oyedotun, T., D., 2014. Shoreline Geometry: DSAS as a tool for historical trend analysis. British

Society for Geomorphology. In Geomorphological techniques, Chap. 3, Sec. 2, pp. 2-12.

Pavlopoulos, K., Evelpidou, N. y Vassilopoulos, A., 2009. Mapping geomorphological environments.

Springer Science & Business Media.

Pedraza Gilsanz, J. D., 1996. Geomorfología: principios, métodos y aplicaciones (No. 551.4

PED).

Posada, B. O. y Buitrago, N. R., 2009. Metodología para el levantamiento de perfiles de playa:

Manual del programa Iberoamericano de Ciencia y Tecnología para el Desarrollo.

RAPIDEYEAG, 2011. Imágenes satelitales Rapideye. Brandenburg an der Havel, Alemania.

Ruíz, A., Tejeda Martínez, A., Miranda Alonso, S. y Flores Zamudio, R.H. 2010. Climatología.

En G. Benítez Badillo y C. Welsh (coords.). Atlas del patrimonio natural, histórico y cultural

de Veracruz. Tomo I. Patrimonio natural. Comisión del Estado de Veracruz para la Conmemoración

de la Independencia Nacional y de la Revolución Mexicana. Xalapa de Enríquez,

pp. 65-84.

Ruiz-Martínez, G., Silva-Casarín, R., Posada-Venegas, G., 2013. Comparación morfodinámica

de la costa noroeste del estado de Quintana Roo, México. Tecnología y Ciencias del Agua,

1(3):47-65.

Servicio Mareográfico Nacional (SMN), 2017. Datos de las estaciones mareográficas de Veracruz.

UNAM, Instituto de Geofísica. Recuperado el 12 de junio de 2019 de http://www.

mareografico.unam.mx/portal/index.php?page=Estaciones.

Serrano García, M. A., 2017. El mapa geomorfológico. Una cartografía geológica peculiar y

útil. Memorias de la Real Sociedad Española de Historia Natural, 2(14):127-153.

Shepard, F. P., 1963. Submarine geology. Harper & Row, New York, 557 p.

Silva-Casarín R., Villatoro-Lacoutre M. M., Ramo-Durón F. J., Pedroza-Paez D., Ortíz-Pérez M.

A, Mendoza-Baldwin E. G., Delgadillo-Calzadilla M. A., Escudero-Castillo M. C., Félix-Delgado

A. y Cid., 2014. Caracterización de la zona costera y planteamiento de elementos

técnicos para la elaboración de criterios de regulación y manejo sustentable. Instituto de

Ingeniería, UNAM, 125 pp.

223

CEMIE-Océano

Spot Image, 2017. Imágenes satelitales Spot 6/7. Sistema de Información Geográfica, S. A. de

C. V., Ciudad de México, 7 imágenes.

Soares, A. G., 2003. Sandy beach morphodynamics and macrobenthic communities in temperate,

subtropical and tropical region-a macroecological approach. PhD thesis, University

of Port Elizabeth.

Stanica, A., Dan, S. y Ungureanu, V. G. 2010. Coastal changes at the Sulina mouth of the Danube

River as a result of human activities. Marine Pollution Bulletin, 55(10-12):555-563.

DOI:10.1016/j.marpolbul.2007.09.015.

Temitope, D. y Oyedotun, T., 2014. Shoreline geometry: DSAS as a tool for historical trend

analysis. Geomorphological Techniques, 3(2.2):12.

Thieler, E, R., Danforth, W. W., 1994. Historical Shoreline mapping (II): Appplication of the

Digital Shoreline mapping an Analysis System (DSMS/DSAS) to Shoreline change mapping

in Puerto Rico. Journal of Coastal Research, 10(3):600-620.

Thieler, E.R., Himmelstoss, E.A., Zichichi, J.L., & Ergul, A. 2009. The digital shoreline analysis

system (DSAS) Version 4.0 – An ArcGis Extension for calculating shoreline change. Open-File

Report. U. S. Geological Survey Report No. 2008-1278. http://woodshole.er.usgs.gov/

projectpages/dsas/version4/

Tomazin, N. y Re, M., 2019. Herramientas de modelación para la gestión costera de la provincia

de Buenos Aires. INA, imfia. Technologies for the design of a regional strategic plan

for the coastal management and adaptation to Climate Change in the Province of Buenos

Aires, 40 pp.

Vidal, C., Losada, M. Á., Medina, R. y Losada, Í., 1995. Modelos de morfodinámica de playas.

Ingeniería del agua, 2(1 Extraordinario):55-74.

Wright, L. D. y Short, A. D., 1984. Morphodynamic variability of surf zones and beaches: a

synthesis. Marine geology, 56(1-4):93-118. DOI:10.1016/0025-3227(84)90008-2.

Zagorski, P., Jarosz, K. y Superson, J., 2020. Integrated assessment of shoreline change along

the Calypsostranda (Svalbard) from remote sensing, field survey and GIS. Marine Geodesy,

43(5), 433-471. DOI:10.1080/01490419.2020.1715516.

224



Anexo 1

Estación: 1 Coordenadas: 20°42'18''N, 97°11'43.5''W 20°42'16.4''N,97°11'43.6''W

Localidad:

Cazones,

Playa Sur

Fecha: 01/06/2018

Escala vertical 1:200

Escala horizontal 1:250

Altitud inicial 0.0 msnm Az 263

Altitud Final 3.43 Elementos hidrodinámicos

Altura máxima 3.43 Rompiente En derrame

Altura total 3.43

1.5 m

Ancho del perfil 14.32 Altura del oleaje

Elementos morfológicos Altura (m) Longitud (m)

Dunas 1.08 4.48

Berma 1.41 9.84

Escarpe ------------ ------------

Otros: ------------ ------------

INSERTAR IMAGEN DEL PERFIL

Tipo de Vegetación Inducida Elementos antrópicos

Ríos o cuerpos de agua ------------ Infraestructura ------------

Caracteristicas del

sedimento Fino-grueso Asentamiento ------------

Uso de suelo Agropecuario

225

CEMIE-Océano

Es tación: 2

Coordenadas: 20°42'12.8''N, 97°11'42.7''W 20°42'14.5''N,97°11'3.8''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

263

Elementos hidrodinámicos

3.09

3.09 Rompiente En derrame

3.09

34.34

Altura del

oleaje

1.5 m

Elementos morfológicos

Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

------------

------------

------------

Longitud (m)

------------

------------

------------

------------

Tipo de Vegetación

Ríos o cuerpos de agua

Caracteristicas del sedimento

Inducida

------------

Fino-grueso

Localidad: Cazones, Playa Sur Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

------------

Elementos antrópicos

------------

Agrícola

226



Es tación: 3


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

241

2.42 Elementos hidrodinámicos


2.42 Altura del

1.5 m

45.7 oleaje


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

------------

------------

------------

Longitud (m)

------------

------------

------------

------------

Tipo de Vegetación



Inducida

------------

Fino-medio

Localidad: Cazones, Playa Sur Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

Cabañas


Localidad rural

Habitacional y turístico

227

CEMIE-Océano

Es tación: 4

Coordenadas: 20°41'0.5''N, 97°11'9.5''W 20°41'00''N,97°11'10.5''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

242

3.28 3.28 Elementos hidrodinámicos

Rompiente En voluta/derrame

3.28 Altura del

33.5 oleaje

1.5 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

1.66

1.62

------------

------------

Longitud (m)

3.3

30.2

------------

------------

Tipo de Vegetación



De duna

------------

Fino-grueso

Localidad: Cazones, Playa Azul Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

Cerca


------------

------------

228



Es tación: 5


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

242



2.14 Altura del

1 m

28.7 oleaje


Dunas

Berma

Escarpe

Otros:

Altura (m)

0.67

1.47

------------

------------

Longitud (m)

3

25.7

------------

------------

Tipo de Vegetación



De duna e inducida

Río Cazones

Medio-fino





Infraestructura

Asentamiento

Uso de suelo

Cerca


------------

------------

229

CEMIE-Océano

Es tación:6


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

233


2.39


2.39

17.9

Altura del

oleaje

1 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

0.48

1.91

------------

------------

Longitud (m)

1.8

16.1

------------

------------

Tipo de Vegetación



De duna, con perturbación

------------

Medio





Infraestructura

Asentamiento

Uso de suelo

Cerca


------------

------------

230



Es tación: 7

Coordenadas: 20°41'19.4''N, 97°11'20.5''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

245



2.05 Altura del

1.5 m

28.4 oleaje


Dunas

Berma

Escarpe

Otros:

Altura (m)

0.44

1.61

------------

------------

Longitud (m)

2.5

25.9

------------

------------

Tipo de Vegetación



Inducida

------------

Medio





Infraestructura

Asentamiento

Uso de suelo

Casas de tablas o de cemento


Construcción rústica


231

CEMIE-Océano

Es tación: 8


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

246



1.78

42.99

Altura del

oleaje

1.5 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

1.78

------------

------------

Longitud (m)

------------

22.99

------------

------------

Tipo de Vegetación



Inducida

------------

Fino-medio

Localidad: Playa Azul Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

Habitacional


Turístico rural

Residencial turístico

232



Es tación: 9

Coordenadas: 20°44.4'49.7''N, 97°11'54.1''W 20°44.4'49.9''N,97°11'54.6''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

232


2.59

2.59 Rompiente En voluta

2.59

21.08

Altura del

oleaje

1.5 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

1.24

------------

------------

Longitud (m)

------------

4.43

------------

------------

Tipo de Vegetación



Matorral de zona costera

------------

Grueso-medio

Localidad: Chaparrales Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

------------


------------

------------

233

CEMIE-Océano

Es tación: 10

Coordenadas: 20°44'56.9''N, 97°11'58''W 20°44'56.5''N,97°11'58.9''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

142


3.1


3.1

31.4

Altura del

oleaje

1.5 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

0.78

1.42

0.32

------------

Longitud (m)

5.28

12.77

7.55

------------

Tipo de Vegetación



De duna

Río Cazones

Grueso-medio

Localidad: Chaparrales Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

------------


------------

------------

234



Estación: 11 Coordenadas: 20°45'7.8''N, 97°12'4.2''W 20°45'7.4''N,97°12'4.9''W

Localidad:

Chaparrales

Fecha: 01/06/2018






Altura total 1.53

2 m



Dunas ------------ ------------

Berma 1.31 9.85

Escarpe 0.54 5.03

Otros: ------------ ------------


Tipo de Vegetación De duna, sin perturbación Elementos antrópicos


Caracteristicas del

sedimento Medio-fino Asentamiento ------------

Uso de suelo ------------

235

CEMIE-Océano

Estación: 12 Coordenadas: 20°45'17.4''N, 97°12'9.3''W 20°45'17.2''N,97°12'10''W

Localidad:

Chaparrales

Fecha: 01/06/2018






Altura total 2.95

2 m



Dunas 2.12 7.53

Berma 0.83 16.04

Escarpe ------------ ------------

Otros: ------------ ------------


Tipo de Vegetación De duna, sin perturbación Elementos antrópicos


Caracteristicas del

sedimento Grueso-medio Asentamiento ------------

Uso de suelo ------------

236



Es tación: 13


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

249


1.07


1.07

17.47

Altura del

oleaje

1.5 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

1.07

------------

------------

Longitud (m)

------------

17.47

------------

------------

Tipo de Vegetación



Inducida

------------

Grueso-fino

Localidad: Chaparralesr Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

Caseríos


Localidad rural

Vivienda

237

CEMIE-Océano

Es tación: 14


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

315



1.01 Altura del

1 m

36.3 oleaje


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

0.64

------------

------------

Longitud (m)

------------

8.87

------------

------------

Tipo de Vegetación



De duna

Río Cazones

Fino

Localidad: Barra de Cazones Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

------------


------------

Pecuario

238



Es tación: 15


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

284


2.04


2.04

42.18

Altura del

oleaje

1 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

0.79

------------

------------

Longitud (m)

------------

10.13

------------

------------

Tipo de Vegetación



De duna

Río Cazones

Fino





Infraestructura

Asentamiento

Uso de suelo

------------


------------

Agropecuario

239

CEMIE-Océano

Es tación: 16


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

276



1.89

41.38

Altura del

oleaje

1 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

0.7

------------

------------

Longitud (m)

------------

6.38

------------

------------

Tipo de Vegetación



------------

Río Cazones

Fino

Localidad:Barra de Cazones Fecha: 01/06/2018




Infraestructura

Asentamiento

Uso de suelo

------------


------------

Agropecuario

240



Es tación: 17

Coordenadas: 20°43'45''N, 97°11'56.1''W 20°43'45.2''N,97°11'58.3''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

263



1.8

66.04

Altura del

oleaje

1 m


Dunas

Berma

Escarpe

Otros:

Altura (m)

------------

0.85

------------

------------

Longitud (m)

------------

9.7

------------

------------

Tipo de Vegetación



------------

Río Cazones

Fino





Infraestructura

Asentamiento

Uso de suelo

------------


------------

Agropecuario

241

242

CEMIE-Océano



243

244

CEMIE-Océano



245

246

CEMIE-Océano



247

248

CEMIE-Océano



249

CEMIE-Océano

Estación: 26


Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

246



14.11 74.45 Altura del

oleaje 0.4m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura

6.25

4.33

1.09

------------

------------

Longitud

17.58

9.15

21.06

------------

------------

Arbustiva

------------

Fino

Localidad: Playa Muñecos

Infraestructura

Asentamiento

Uso de suelo

Fecha: 05/06/2018




Caminos y casas abandonadas


------------

------------

250



Estación: 27 Coordenadas: 19°45'18.8''N, 96°24'53.3''W

Localidad: Playa Muñecos

Fecha: 05/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

242



15.64

60.4

Altura del

oleaje

0.5 m


Dunas

Berma

Escarpe

Otros:

Altura

0.42

1.08

3.26

7.31

0.88

------------

------------

3.52 INSERTAR IMAGEN DEL PERFIL

Longitud

5.62

7.12

1.8

27.04

------------

------------

251

252

CEMIE-Océano



Estación: 29 Coordenadas: 19°45'03.3''N, 96°24'40.5''W 19°45'02.5''N,96°24'42.1''W Localidad: Playa Muñecos

Fecha: 05/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

236 INSERTAR IMAGEN DEL PERFIL



10.46

100.07

Altura del

oleaje

10.5 m

Eelementos morfológicos

Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura

2.07

2.99

6.32

1.14

------------

------------

Longitud

8.35

18.27

17.17

23.94

------------

------------

De duna, sin perturbación

------------

Fino

Infraestructura

Asentamiento

------------


------------

253

CEMIE-Océano

Estación: 30

Coordenadas: 19°45'03.3''N, 96°24'40.5''W 19°45'02.5''N,96°24'42.1''W Localidad: Playa Muñecos

Fecha: 05/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Longitud del perfil


Dunas

Berma

Escarpe

Otros:

0.0 msnm Az

241



5.69

53.24

Altura

1.34

2.16

3.01

------------

------------

Altura del

oleaje

Longitud

6.38

5.7

38.96

------------

------------


Tipo de Vegetación



Arbustiva

------------

Fino

Infraestructura

Asentamiento

Uso de suelo

------------


------------

------------

254



Estación: 31

Coordenadas: 19°41'4.8''N,96°23'53.4''W;19°41'8.1''N,96°23'58.9''W Localidad: Villa Rica (Nte.)

Fecha: 05/06/2018

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

204

10.88msnm Elementos hidrodinámicos

10.88m Rompiente En voluta

10.88m

100.51

Altura del

oleaje

0.5m

Escala Vertical 1:200

Escala Horizontal 1:500


Dunas

Berma

Tipo de Vegetación



Altura m

0.78

1.03

6.48

2.59

Longitud m

9.65

18.77

19.33

51.12

Vegetación de duna sin perturbación

Al norte desembocadura del río Limón

Fino

Infraestructura

Asentamiento

Uso de suelo

--------------------------------------------


--------------------------------------------

--------------------------------------------

255

CEMIE-Océano

Estación: 32

Coordenadas: 19°41'4.8''N, 96°23'53.4''O 19°41'3.4''N,96°23'56.7''W

Localidad: Villa Rica, Norte. Fecha: 05/06/2018

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

243°



9.83

108.17

Altura del

oleaje

0.5m


Escala Vertical 1:200


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación


Caracteristicas del

sedimento

Altura (m)

2.03

3.67

1.58

1.33

2.69

--------

--------

Longitud (m)

6.08

11.91

12.08

16.65

32.33

--------

--------

De duna sin perturbación

----------------

Fino de origen fluvial

Infraestructura

Asentamiento

Uso de suelo

----------------


----------------

----------------

256



Estación: 33

Coordenadas: 19°38'26.5''N,96°23'36.1''W; 19°38'25''N,96°23'38.8''W

Localidad: Villa Rica, Sur.

Fecha: 06/06/2018


Escala Veritical 1:200

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

239°



6.04

99.09

Altura del

oleaje

1.40m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura (m)

2.3

4.14

2.02

-----------

-----------

Longitud (m)

28.48

22.25

44.53

-----------

-----------

De dunas con perturbación de casuarina Elementos antrópicos

Laguna El Farallón

Infraestructura

------------------------------

Fino de origen terrestre (fluvial)

Asentamiento

------------------------------

Uso de suelo

------------------------------

257

CEMIE-Océano

Estación: 34

Coordenadas: 19°38'35.7''N,96°23'41.2''W;19°38'345''N,96°23'41.2''W Localidad: Villa Rica (Sur)

Fecha 06/06/2018

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

251


6.45m Rompiente En derrame

6.45m

102.2m

Altura del

oleaje

1.40m




Dunas

Berma

Tipo de Vegetación



Altura m

1.68

0.96

1.11

2.67

2.99

Longitud m

10.38

6.71

11.76

14.53

45.6

Vegetación natural de duna perturbada con

cerca viva de casuarina

Laguna El Farallón

Fino

Infraestructura

Asentamiento

Uso de suelo


---------------------------

---------------------------

---------------------------

258



Estación: 35

Coordenadas: 19°38'43.3''N,96°23'45''W;19°38'42.6''N,96°23'47.5''W Localidad: Villa Rica (Sur)

Fecha: 06/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

250



3.64m

76.18m

Altura del

oleaje

1.40m


Dunas

Berma

Tipo de Vegetación



Altura m

2

0.54

2

1.51

Longitud m

17.86

6.81

23.33

25

Natural de duna perturbada con casuarina

Laguna El Farallón

Fino

Infraestructura

Asentamiento

Uso de suelo

----------------------------------------


----------------------------------------

----------------------------------------

259

CEMIE-Océano

Estación: 36

Coordenadas: 19°38'56''N,96°23'48.8''W;19°38'55.5''N,96°23'51.6''W Localidad: Villa Rica (Sur)

Fecha: 06/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

255



5.81m

85.67m

Altura del

oleaje

1.40m


Dunas

Berma

Tipo de Vegetación



Altura m

2.61

0.97

3.08

2.54

Longitud m

18.82

13.56

10.75

37.76

Natural de duna

Laguna El Llano

Fino

Infraestructura

Asentamiento

Uso de suelo

----------------------------------


----------------------------------

----------------------------------

260



Estación: 37

Coordenadas: 19°39'4.2''N,96°23'50.1''W;19°39'03.8''N,96°23'53.7''W Localidad: Villa Rica

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

265m

6msnm Elementos hidrodinámicos

6m Rompiente En derrame

6m

108.6

Altura del

oleaje

0.80m


Dunas

Berma

Otros

Tipo de Vegetación



Altura m

3.46

4.18

2.42

----------

Longitud m

31.25

14.59

34.85

----------

Natural de duna

Laguna El Llano

Fino

Infraestructura

Asentamiento

Uso de suelo

Fecha: 06/06/2018



---------------------------------------


---------------------------------------

---------------------------------------

261

262

CEMIE-Océano



Estación: 39

19°39'39.6''N,96°23'55''W;19°39'39.5''N,96°23'56.4''W Localidad: Villa Rica

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

257



2.57m

43.3m

Altura del

oleaje

0.5m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura m

2.37

0.9

0.26

---------------

Longitud m

15.4

27.9

0

---------------

Natural de duna sin perturbación

Laguna El Llano

Fino

Infraestructura

Asentamiento

Uso de suelo

Fecha: 06/06/2018



----------------------------------------------------


----------------------------------------------------

----------------------------------------------------

263

CEMIE-Océano

Estación: 40

Coordenadas: : 19°39'56.8''N, 96°23'53.8''W 19°39'58.2''N,96°23'56.4''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

275

3.21 msnm Elementos hidrodinámicos

3.21 m Rompiente En derrame

3.21 m

89.39 m

Altura del

oleaje

< 0.5 m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura (m)

1.7

0.9

------------

------------

Longitud (m)

9.4

11.09

------------

------------


Boca de la laguna El Llano

Fino-medio

Localidad: Villa Rica, Sur. Fecha: 06/06/2018




Infraestructura

Asentamiento

Uso de suelo

Turística (palapas)


Localidad rural


264



Estación: 41

Coordenadas: 19°40'14.9''N,96°23'55''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az (°) 273

2.44 msnm Elementos hidrodinámicos

2.44 m Rompiente En derrame

2.44 m

40.97 m

Altura del

oleaje < 0.5 m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura (m)

0.67

1.02

-------------

-------------

Longitud (m)

3.27

8.9

-------------

-------------


---------------------

Fino-medio

Localidad: Villa Rica

Fecha: 06/06/2018




Infraestructura

Asentamiento

Uso de suelo

Turística (palapas)


Localidad rural


265

CEMIE-Océano

Estación Coordenadas Localidad

Fecha: 08-06-2018

42 18°42'14.5''N, 95°11'46.7''W Punta Roca Partida

18°42'12.7''N, 95°11'46.1''W

Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil 71.08

0.0 msnm Az

179


Rompiente En voluta

Altura del

oleaje

1


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura Longitud tipo

------- ------- embrionarias y estabilizadas

1.09 3.12 Subhorizontal

------- ------- no

arbustiva

no

arenoso

Infraestructura

Asentamiento

Uso de suelo


no

no

------------------------------

266



Estación Coordenadas: 18°42' 14.8''N, 95°11'48.6''W Localidad

Fecha: 08-06-2018

43 18°42'13.3''N,95°11'48.6''W Punta Roca Partida

Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil

0.0 msnm Az

168


Rompiente En voluta

46.77

Altura del

oleaje

1


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación




3.98 8.38 embrionarias y estabilizadas


----- ----- no

----- -----

arbustiva

no

arenoso

Infraestructura

Asentamiento

Uso de suelo


no

no

267

CEMIE-Océano

Estación Coordenadas 18°42'15.4'N, 96°11'53.6''W Localidad

Fecha: 08-06-2018


Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



0.0 msnm Az

200


Rompiente en voluta

37.68

Altura del

oleaje




----------- ----------- no

----------- -----------

arbustiva

no

arenoso

1

Infraestructura

Asentamiento

Uso de suelo



no

no

pecuaria

268



Estación

Coordenadas 18°42'14.5''N, 95°11'46.7''W Localidad

Fecha: 08-06-2018


Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



0.0 msnm Az

197



34.76

Altura del

oleaje


1 14.93 embrionarias y estabilizadas


no

arbustiva

no

arenoso

1

Infraestructura

Asentamiento

Uso de suelo



no

no

pecuario

269

CEMIE-Océano

Estación Coordenadas 18°41'53.8''N, 95°10'54.5''W Localidad

Fecha: 08-06-2018

46 18°41'52.5''N, 95°10'56''W Punta Roca Partida

Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



0.0 msnm Az

226



54.74

Altura del

oleaje




no

arbustiva

no

arenoso

1

Infraestructura

Asentamiento

Uso de suelo


no


no

270



Estación

Coordenadas 18°41'45.9''N, 95°10'48.1''W Localidad: Punta Roca Patida

Fecha: 08-06-2019

47 18°41'451''N, 95°10'48.9''W

Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



0.0 msnm Az

218



30.24

Altura del

oleaje


2.74 0.9 embrionarias y


no

arbustiva

no

arenoso

1

Infraestructura

Asentamiento

Uso de suelo


no


no

271

CEMIE-Océano

Estación Coordenadas 18°41'41.4''N,95°10'43.4''W Localidad: Punta Roca Partida

Fecha: 08-06-2018

48 18°41'40.7''N,95°10'44.3''W

Altitud inicial

Altitud Final

Altura máxima

Longitud del perfil


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



0.0 msnm Az

212



42.41

Altura del

oleaje


1.27 5.59 embrionarias y


no

arbustiva

no

arenoso

1

Infraestructura

Asentamiento

Uso de suelo



no

no

272




Localidad: Playa Hermosa

Fecha: 08/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az




2.6

46.98

Altura del

oleaje

2 m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura Longitud

0.95 3.18

0.92 15.45

------------ ------------

------------ ------------


------------

Medio-fino

Infraestructura

Asentamiento

Uso de suelo

Construcciones de cemento


------------

------------

273

CEMIE-Océano


Localidad: Playa Hermosa

Fecha: 08/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az




3.43

38.81

Altura del

oleaje

2 m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura Longitud

1.58 9.2

1.11 12.99

------------ ------------

------------ ------------


------------

Medio-fino

Infraestructura

Asentamiento

Uso de suelo

Construcciones de cemento


------------

------------

274



Estación: 51

Coordenadas: 18°38'41.3''N,95°05'51.3''W;18°38'40.3''N,95°05'51.3''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

168

1.3m Elementos hidrodinámicos


1.3m

17.28m

Altura del

oleaje

0.5m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura m

-----------

1.30

-----------

-----------

Longitud m

----------

17.28

---------

---------

Sin vegetación

Desembocadura modificada del Rio Cold y

Máquina

Fino-medio

Localidad: Montpío

Fecha: 08/06/2018

Escala vertical

1:200


Infraestructura

Turistica (Palapas)

Asentamiento

Uso de suelo

Localidad rural

Habitaciónal y turístico

275

CEMIE-Océano

Estación: 52

Coordenadas:18°38'40.8''N,95°05'56.6''W;18°38'40.1''N,95°05'56.5''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

178



1.87m

20.36m

Altura del

oleaje

0.5m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura m

--------------

2.55

-----------

-----------

Longitud m

----------

11.64

-----------

-----------

Sin vegetación

Desemboradura modificada del rio Máquina y

Cold

Fino-medio

Localidad:Montepío

Fecha: 08/06/2018


Infraestructura

Turística (Palapas)

Asentamiento

Uso de suelo

Localidad rural

Habitacional y turística

276



Estación: 53

Coordenadas:18°38'46.9''Nn,95°6'21.8''W;18°38'45.5''N,95°6'22.4''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

187



9.17m

33.84m

Altura del

oleaje

≥2.5m


Dunas

Berma

Escarpe

Tipo de Vegetación



Altura m

4.21

1.25

0.72

0.34

Longitud m

11.42

10.06

0

0

Vegetacion dunar (matorral 1.5m-4m)

Al sur presencia de río

Fino-medio

Localidad: Montepío

Fecha: 08/06/2018



Infraestructura

Asentamiento

Uso de suelo

--------------------------------


-------------------------------

--------------------------------

277

CEMIE-Océano

Estación: 54

Coordenadas: 18°38'50.4''N,95°6'31.1''W;18°38'49.2''N,95°6'31.4''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

201



8.67m

44.59

Altura del

oleaje

≥1.5m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura m

4.96

1.63

0.23

--------------

Longitud m

12.93

31.66

0

--------------

Vegetación dunar arbustiva (0.5m-3m)

---------------------------

fino-medio


Fecha: 08/06/2018



Infraestructura

Asentamiento

Uso de suelo

-------------------------


-----------------------------

------------------------------

278



Estación: 55

Coordenadas:18°38'54.2''N,95°6'40.7''W;18°38'52.4''N,95°6'41.7''W Localidad: Montepío

Fecha : 08/06/2018



Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

204



5.60m

66.85

Altura del

oleaje

≥1.5m


Dunas

Berma

Escarpe

Otros:

Tipo de Vegetación



Altura m

3.49

2.18

0.52

----------------

Longitud m

36.56

30.29

0

----------------

Vegetacion de duna arbustiva (0.5-3m)

-------------------------

fino-medio

Infraestructura

Asentamiento

Uso de suelo

------------------------------------------


------------------------------------------

------------------------------------------

279

CEMIE-Océano

Estación:56

Coordenadas: 18°38'55.5''N,95°6'44.2''W;18°38'53.8''N,95°6'45''W

Altitud inicial

Altitud Final

Altura máxima

Altura total

Ancho del perfil

0.0 msnm Az

202



4.54m

71.67m

Altura del

oleaje

≥1.5m


Dunas

Berma

Escarpe

Tipo de Vegetación



Altura m

6.94

1.93

0.33

0.78

Longitud m

42.13

20.7

2.2

2

Vegetacion de duna arbustiva (0.5-3m)

----------------------------------

fino-medio


Fecha: 08/06/2018



Infraestructura

Asentamiento

Uso de suelo

camino rural atrás de la duna


Localidad rural

--------------------------------------------------

280

Libro-Energia-Undimotriz-Cap 4

Create successful ePaper yourself

Delete template?

Save as template?