INTRODUCCIÓN

En años recientes, diversas investigaciones han contribuido con la comprensión de la dinámica que afecta a los litorales en relación con los impactos que ocurren en las zonas costeras del mundo (^{Maulud et al., 2018}; ^{Silva et al., 2014}). Uno de los temas que mayor preocupación causa en dichas áreas es la variación de la línea de costa, que se ha relacionado con el incremento en el nivel del mar, y se considera un factor determinante en el deterioro de los ecosistemas costeros (^{Cazenave y Cozannet, 2014}; ^{Orejarena-Rondón et al., 2019}).

Según ^{Dinesen (2019)} y ^{Labuz (2015)}, la variación de la línea de costa guarda estrecha relación con la erosión que afecta los ecosistemas de playas arenosas del planeta, lo que influye en su equilibrio ecológico. Sus efectos no solo presionan la conservación de diversos ecosistemas litorales, sino que afecta la morfología de la costa y contribuye con el deterioro de la dinámica costera en general (^{Dugan et al., 2011}).

Precisamente, la erosión de las playas arenosas se ha incrementado por la acción del cambio climático, siendo uno de los agentes que mayores efectos ocasiona (^{Tsoukala et al., 2015}; ^{Vousdoukas et al., 2020}). Factores como el deterioro de las barreras coralinas, afectación de los organismos productores de arena, el aumento en la intensidad y periodicidad de los eventos meteorológicos extremos, pérdida progresiva e intensa de las funciones de amortiguamiento y protección costera que cumplen los manglares (^{Pennings et al., 2021}), siendo estas parte de las causas que la comunidad científica ha relacionado con la variación de la línea de costa a nivel mundial (^{Dugan et al., 2008}; ^{Cazenave y Cozannet, 2014}).

Otros agentes precursores del deterioro de la costa son el cambio de uso del suelo, la sedimentación, sobre extracción de materiales de ríos, que inciden en la dinámica sedimentaria de la que depende la conformación de playas arenosas (^{Tsoukala et al., 2015}). A estos agentes se suma la construcción de obras de alto impacto, como puertos de embarque, complejos turísticos y obras de contención (^{Dugan et al., 2008}; ²⁰¹¹; ^{Silva et al., 2014}; ^{Vargas et al., 2016}), muchos de estas costosas obras incluso se construyeron con el fin de funcionar como barreras de protección de la costa, aun cuando su efecto parece ser contrario (^{Petrişor et al., 2020}; ^{Sanjaume y Pardo-Pascual, 2005}).

En Costa Rica, al igual que otros países, se exhibe la pérdida de costa en playas arenosas del litoral Pacífico y Caribe, siendo evidente el cambio en el paisaje y los daños ocasionados a infraestructuras (^{Lizano, 2007}). Sin embargo, la investigación que se realiza en este tema es reciente y con pocos estudios que evalúen las tasas de erosión o acreción y sus relaciones con el deterioro del litoral.

En las playas del Caribe sur (CS) de Costa Rica se han observado cambios progresivos en el desplazamiento de la línea de costa en la última década que parecen agravarse con el ingreso del mar hacia el continente. En este litoral se han registrado secciones críticas donde los procesos erosivos se extienden entre el 55%-80% (^{Lizano, 2013}; ^{Piedra-Castro et al., 2013}). Esta situación se agrava con las características que mantiene la plataforma continental, la cual es clasificada por ^{Bergoeing (2017)} como angosta a lo largo del caribe de Costa Rica, hasta la Laguna Chiriquí en Panamá.

La variación de la línea de costa del CS es causada por la influencia en la dinámica sedimentaria y el consecuente ingreso del mar hacia el continente, lo que constituye una amenaza creciente. Sin embargo, la investigación en este tema en el litoral caribeño de Costa Rica es incipiente. De ahí que sean escasas las investigaciones en el CS relacionadas con este tema, por lo que es fundamental su análisis considerando la influencia que tiene sobre la dinámica sedimentaria el ingreso del mar hacia el continente. El objetivo de esta investigación fue analizar la variación espacial y multitemporal de la línea de costa en playas arenosas del Caribe sur de Costa Rica.

METODOLOGÍA

Área de estudio

El Caribe sur de Costa Rica (CS) se extiende entre el centro de la provincia de Limón (N 10° 00’29"; W 83° 15’06), y la desembocadura del río Sixaola, en el cantón de Talamanca (N 09°25'40"; W 83°01'20") (Figura 1).

Figura 1 Sectores del litoral Caribe Sur, Costa Rica con detalle de las zonas de vida. 

En la sección sur del litoral Caribe, la velocidad promedio del viento es de 3 ms^-1. Los valores máximos ocurren de noviembre a marzo, y entre julio y agosto, mientras que los mínimos ocurren entre abril y mayo, con otro período de disminución que se extiende de setiembre y octubre (^{Pérez-Briceño et al., 2017}; ^{Quesada, 2016}; ^{Ricaurte-Villota et al., 2017}), de ahí que se les clasifique como de tipo bimodal.

Las corrientes marinas exhiben circulación ciclónica (movimiento antihorario) que se extiende hasta Colombia. El agua superficial sobre el borde de la plataforma continental presenta un movimiento hacia el sur-sureste y al frente de la provincia de Limón. El oleaje se dirige hacia el noreste durante todo el año con incrementos en la altura y el periodo entre olas, vinculados con la acción de los vientos alisios (^{Lizano, 2018}; ^{Ricaurte-Villota et al., 2017}).

Tasas de erosión de las playas arenosas

Se utilizó la extensión Digital Shoreline Analysis System (DSAS 4.3) para ArGis 10.3 con el fin de obtener las tasas de erosión de las playas arenosas del Caribe sur. Se unificaron las líneas de cada periodo en un solo archivo de Geodatabase siguiendo el flujo de trabajo propuesto por ^{Himmelstoss et al., (2018)}. Se utilizó como línea base la delimitación de la zona de protección de 200 m, establecida en Ley sobre la Zona Marítimo Terrestre N° 6043, de Costa Rica de 1977, que demarca un área de protección de dicha amplitud, medida desde la pleamar ordinaria.

Se generaron 1547 transectos ortogonales en intervalos de 50 m a lo largo de esta delimitación de la zona de protección. Para evaluar las variaciones en la posición de la línea de costa se consideraron las estadísticas End Point Rate (EPR), que describen los cambios en la línea costera, y cuyo valor representa la tasa de cambio. Además, se determinó el Net Shoreline Movement (NSM) que, junto con el indicador anterior, compara la variación total que experimenta la línea costera en el tiempo. Para ambos los signos negativos representan pérdida por erosión, y los signos positivos ganancia o estabilidad por acreción (^{Sarwar y Woodroffe, 2013}; ^{Himmelstoss et al., 2018}).

Los procesos erosivos y de acreción de las playas arenosas fueron clasificados según la evolución costera con base en los criterios de ^{Cifuentes-Ossa et al. (2017)}, ^{Del Río et al. (2013)} y ^{Natesan et al. (2015)}.

La comprobación del proceso de desplazamiento de línea costera EPR y el NSM se validó mediante el análisis de Clúster con una matriz de 700 puntos georreferenciados por año que se obtuvieron de los monitoreos in situ, mismos que se registran anualmente desde el año 2012 en nueve playas del CS (^{Piedra-Castro et al., 2021}).

Una vez seleccionados los puntos con erosión, se desarrolló una base de datos espacial que incluyó las características de la morfología y porcentajes de erosión de playa. Todos estos puntos se analizaron y se procesaron con el algoritmo Getis-Ord Gi* (QGIS), misma que agrupa espacialmente las entidades con valores altos o bajos para la variable de interés y permite obtener el un indicador del desplazamiento de la línea de costa de las playas arenosas. Para la extracción del clúster se utilizaron como referencia media cinco años (2017-2021) y la década completa (2012-2021).

Adicionalmente, se realizó el perfil topográfico de las playas arenosas con el fin de establecer la relación cualitativa entre ellas, la tasa de erosión y la pendiente de la costa. Para esto se recopiló la información de campo obtenida con el método de EMERY, en el cual se mide la distancia, ángulo y elevación de un segmento hasta donde existe cambio de pendiente hasta el observador, quien verifica y anota la elevación con respecto al horizonte (Figura 2).

Figura 2 Perfil topográfico de playas arenosas del Caribe sur de Costa Rica. 

RESULTADOS

Entre 1960-2023, el 70,14% de los transectos ortogonales evidenciaron tendencia general hacia la erosión, principalmente hacia el sur del litoral. Destaca el EPR obtenido para los sectores del Parque Nacional Cahuita, Puerto Viejo, Manzanillo y Gandoca, por la severidad que presentan. Otros sectores, como Estrella y Vizcaya, evidenciaron tendencia estable, con prevalencia hacia la acreción, misma que se concentró en las desembocaduras de los ríos Bananito, Estrella y Home Creek.

Las playas arenosas del CS son dinámicas y evidenciaron variación de la espacial y temporal de la línea costera para toda el área analizada. Se determinó que el balance general neto para los 63 años experimentó variación del NSM de moderados a estables, tanto por erosión como por la acreción. Este desplazamiento de la línea de costa fue crítico hacia el sureste, afectando tanto las costas con barreras naturales como aquellas en las que no existe tal protección (Figura 3).

Figura 3 Valores de tasas erosivas y de acreción promedio para los sectores evaluados en playas arenosas del Caribe sur de Costa Rica. 

El NSM, para 1960-2023, correspondió a erosión con valores promedio que alcanzaron los -6.85 m/año. En el caso del desplazamiento por acreción, este se ubicó en 7.01 m/año, siendo el sector CSB7 en el área cercana a la desembocadura del río Sixaola el que mayor aporte tuvo para este factor (Tabla 1).

Entre 1960-1976, el NSM por erosión en general para toda el área de estudio fue de -5.36 m/año, siendo la playa de Gandoca, una de las que evidenció mayor desplazamiento. Para este mismo período, el NSM por acreción fue de 2.96 m/año. Para 1998-2005 el desplazamiento de la línea costera por procesos de erosión alcanzó los -5.68 m/año, siendo Gandoca y Estrella los sectores que presentaron afectaciones importantes. En el caso de la acreción, se observaron tasas de 5.75 m/año, siendo la zona de Cocles y Estrella las que mayor contribución tuvieron por este proceso (Tabla 1).

Los signos negativos de la Tabla 1 se refieren al desplazamiento de la costa por erosión, lo cual se traduce en la reducción de la playa. Los signos positivos denotan movimiento neto por acreción. Los valores de la tabla se establecen en esos rangos obtenidos de la totalidad de los transectos ortogonales de cada sector de playa evaluado.

Los principales hallazgos en cada sector se describen a continuación:

Gandoca-Manzanillo (REGAMA). En los sectores CSB6 y CSB7, el cambio en la posición de la línea de costa entre 1960-2023 correspondió a procesos erosivos en 90% de los casos. En el período 1976-1998 y en el 2005-2011 los sectores presentaron incrementos en las áreas de erosión entre -0.61 y - 2.07 m/año, situándose como dos de los sectores con mayor pérdida de costa. En el período total, el desplazamiento de la línea por erosión se clasifica como estable, con variaciones entre -0.44 y -1.74 m/año. El NSM mostró un mínimo de acreción de 39.26 m. Solo el 5% de los 175 transectos ortogonales experimentaron tasas de acreción entre 0.01 y 0.23 m/año (Figura 4).

Figura 4 Desplazamiento de la línea de costa entre 1960-2023 por procesos erosivos y de acreción en el Refugio Mixto de Vida Silvestre Gandoca-Manzanillo, Limón Costa Rica. A) HotSpot 2012-2021; B) HotSpot 2017-2021. 

En Manzanillo, en el sector CSB5, la tendencia al igual que la anterior fue hacia la erosión de la playa. El 67% de los transectos evidenció tasas de EPR entre -2.07 y -0.01 m/año. Los valores más críticos de erosión se observaron en 2005-2011. En el resto de los periodos el comportamiento presentó poca variación con respecto a la acreción. Para el período total entre 1960-2023 la tasa de acreción alcanzó un máximo de 2.4 m/año, mientras que el NSM fue de 42.64 m en sectores cercanos a las desembocaduras de los ríos (Figura 4).

Mediante el análisis de puntos calientes en el REGAMA, se evidencia que el sector noroeste del CSB7 experimenta un proceso erosivo constante y repetitivo que se ha documentado efectivamente desde el 2012. En los otros dos transectos CSB6 y CSB5, se documentó este comportamiento, aunque en menor magnitud que en el anterior (Figura 4 A y 4 B).

Cocles-Puerto Viejo. De los 151 transectos ortogonales de CSB4 y CSB3, el 90% de ellos evidenciaron procesos erosivos en los períodos 1960-1976 y 2005-2011, con valores por encima de -0.06 m/año. En el período total 1960-2023, entre el 67%-83% de estos transectos denotaron erosión según las tasas del EPR, con variación entre el -0.2 m/año y -0.4 m/año, lo que clasifica este desplazamiento como estable. En contraste con lo anterior, la acreción se clasificó como moderada con valores entre 0.01 y 3.51 m/año entre 2011-2023, en sitios puntuales de ambos sectores. El NSM fue de 39.16 m para los lugares donde ocurre este tipo de desplazamiento de la línea de costa (Figura 5).

Figura 5 Desplazamiento de la línea de costa entre 1960-2023 por procesos erosivos y de acreción en Puerto Viejo, Limón Costa Rica. A) HotSpot 2012-2021; B) HotSpot 2017-2021. 

En ambos sectores CSB4 y CSB3, el análisis HotSpot mostró que las tasas erosivas que se mantienen a lo largo del tiempo incluso alternan en las áreas cercanas a Puerto Viejo, y en menor magnitud hacia el sureste, en dirección a Cocles (Figura 5 A y 5 B).

Parque Nacional Cahuita-Estrella. En el sector CSB2, entre el 80% y el 85% de los transectos ortogonales correspondieron a tasas negativas de EPR, para los períodos 2005-2011 y entre 2011-2023 con valores promedio de -0.65 m/año, catalogando la erosión como estable. Las áreas críticas por pérdida de playa se concentraron en punta Cahuita y en cerca del estero Perezoso. En el sector CSB1, los mayores desplazamientos por EPR negativo se presentaron entre 1960-1976 y 1998-2005, con valores del 74% al 82% del total de transectos respectivamente. En cuanto al NSM, este alcanzó los -0.93 m/año para el período total 1960-2023, observando leve tendencia de los datos hacia el desplazamiento de la línea costera por erosión en ambos sectores (Figura 6).

Figura 6 Desplazamiento de la línea de costa entre 1960-2023 por procesos erosivos y de acreción en el Parque Nacional Cahuita, Limón Costa Rica. A) HotSpot 2012-2021; B) HotSpot 2017-2021. 

Los puntos críticos de erosión de costa en el Parque Nacional Cahuita se mantienen espacialmente en las mismas áreas, siendo el sector de Puerto Vargas donde se observa mayor afectación por este tipo de eventos (Figura 6 A y 6 B)

Estrella-Limón. En el tramo litoral CSA3 de los 131 transectos ortogonales, entre el 68% y 83% obtuvieron valores EPR erosivos para los períodos 1998-2005 y 2011-2023, respectivamente, con tasas que oscilaron entre -0.7 m/año a -0.8 m/año. Las áreas que presentaron acreción se ubicaron en las desembocaduras de los ríos Estrella y Bananito principalmente. En el caso de los sectores CSA2 y CSA1, los EPR de tipo erosivo se presentaron en el 69% y el 74% de los transectos en los períodos de análisis que se extienden desde el 2005 y hasta el 2023, con valores que oscilaron entre estables a moderadas de -3.2 m/año en las inmediaciones del río Bananito (Figura 7).

Figura 7 Desplazamiento de la línea de costa entre 1960-2023 por procesos erosivos y de acreción en Limón, Costa Rica. A) HotSpot 2012-2021; B) HotSpot 2017-2021. 

Se destaca el posible efecto que tiene la Zona de Convergencia Intertropical en la variación de los perfiles de costa, ya que se ha demostrado que esta puede incidir de manera directa en las interacciones océano-atmósfera y tierra-atmósfera a escala local (^{Waliser y Jiang 2015}), por lo que se convierte en un factor adicional que puede potencialmente afectar las playas arenosas, al favorecer la formación de eventos meteorológicos intensos y más frecuentes, los cuales golpean con mayor fuerza estas áreas de la costa provocando variaciones en el perfil topográfico.

Al igual que ocurre con las playas de Gandoca y Estrella, las que se ubican en estos tramos CSA3, CSA2 y CSA1 mantienen fondos no consolidados y blandos, sin formaciones geológicas que favorezcan la reducción de la energía de ingreso del oleaje, de ahí que su afectación sigue patrones puntuales en los que ocurre tasas de EPR elevadas, siendo focalizada la erosión (Figura 7).

En concordancia con lo anterior, el análisis de clúster evidenció que los sectores CSA1, CSA2 y CSA3 mantienen áreas críticas por erosión, mismas que fueron identificadas por el EPR. Al igual que se desprende del resultado del DSAS, el HotSpot detectó las áreas de las desembocaduras de los ríos Estrella, Bananito, Estero Negro y Vizcaya como puntos fríos o de baja significancia. Esta condición prevalece debido al aporte y deposición de los sedimentos que ingresan en el sistema sedimentario de playas arenosas de origen terrígeno que se forman a la salida de estos ríos (Figura 7 A y 7 B).

Con respecto a los perfiles topográficos de las playas arenosas del CS, se destaca la variación que mantienen en la pendiente, la cual ha ido incrementando su ángulo, siendo Gandoca y Puerto Vargas los que presentaron la mayor reducción en el ancho e inclinación. Estas playas se mantienen una tendencia de tipo reflexiva (^{McLachlan y Defeo, 2017}).

En Vizcaya, al igual que en el resto del tramo CSA1, CSA2 y CSA3, se observa el mismo comportamiento de playa tipo semi-reflexiva, con inclinación reducida y con anchos que superan los 40 m (Figura 8).

Figura 8 Perfiles topográficos de playas arenosas del Caribe Sur de Costa Rica 2016-2019; A) Gandoca; B) Puerto Vargas; C) Estrella; D) Vizcaya. 

DISCUSIÓN

A nivel mundial, la variación de la línea costera se ha vinculado a los efectos que son provocados y multiplicados por la acción del cambio climático, como resultado del incremento en el nivel medio del mar (^{Tsoukala et al., 2015}; ^{Vousdoukas et al., 2020}). Una de las consecuencias fuertes y que mayor cambio provoca en las playas arenosas es la erosión de la costa (Dinesen, 2018; ^{Labuz, 2015}). Precisamente, el retroceso de la costa prevé valores mundiales de pérdida alrededor del 80% de las zonas costeras del mundo en escenarios RCP 4.5 y RCP 8.5 (^{Dugan et al., 2008}; ^{Labuz, 2015}; ^{López et al., 2017}; ^{Vousdoukas et al., 2020}). ^{Sytnik, Del Río, Greggio y Bonetti (2018)} destacan que la información cuantitativa sobre la posición de la costa en distintos períodos de tiempo es vital para sustentar acciones puntuales de la gestión costera, incluyendo el ingreso del mar hacia el continente y la colocación de defensas en la zona costera.

Durante el periodo comprendido entre 1960-2023, el CS ha experimentado desplazamiento de la línea costera por erosión, que se ha intensificado en algunas de las playas arenosas que conforman este espacio. Respecto al escenario a nivel mundial, existen factores en el CS que pueden vincularse a estos impactos, por ejemplo, el intenso arrastre de sedimentos (^{López et al., 2017}; ^{Ordaz et al., 2016}; ^{Ricaurte-Villota et al., 2017}), deterioro generalizado de los productores de arena, reducción del efecto barrera de arrecifes coralinos, poca o nula presencia de reservas de arena, dunas incipientes, incremento de fenómenos climáticos extremos como el de ENOS (^{Mentaschi et al., 2018}), afectación del sistema sedimentario por efecto de infraestructura humana (^{Cazenave y Cozannet, 2014}; ^{Dungan et al., 2011}; ^{Silva et al., 2014}).

Entre 1960 y 2023, más del 70% de la línea de costa del CS evidenció cambios severos de posición principalmente relacionado con la erosión. Se obtuvieron valores del EPR por cada playa del CS, y se logró determinar cómo el NSM que ocurre es debido a erosión por encima de los eventos de deposición. Estas variaciones espacio-temporales de línea costera destacan los datos obtenidos para el Parque Nacional Cahuita, Puerto Viejo, Manzanillo y Gandoca. Los eventos de variación de la línea coinciden con otros documentados en países como Cuba, Colombia, Puerto Rico, Haití y República Dominicana, donde el desplazamiento de la línea costera afecta ecosistemas e infraestructuras humanas (^{Pacheco et al., 2015}; ^{Paneque et al., 2010}; ^{Rangel-Buitrago et al., 2015}; ^{Rangel y Posada, 2013}; ^{Silva et al., 2014}).

Adicionalmente, las playas arenosas del CS han experimentado cambios constantes de la pendiente y del ancho. Siendo evidente en Puerto Vargas, Manzanillo, Puerto Viejo y Gandoca, la variación y la temporalidad corta que ocurre son indicadores efectivos del avance de la erosión de la costa y de la vulnerabilidad a la que se enfrentan. Al respecto, ^{Aizam (2021)}, concluye que ambos factores juegan un rol importante en la mitigación del efecto que tiene la entrada de oleajes de alta energía y la consecuente pérdida de grandes volúmenes de sedimentos.

En ese mismo lapso, la tasa EPR de acreción varió principalmente en áreas cercanas a las desembocaduras de los ríos Estrella, Bananito, Estero Negro y Vizcaya. Entre 1960 y 1998, hubo una reducción del desplazamiento de línea de costa de 30% aproximadamente, que afectó en sectores como Manzanillo (CSB5), Puerto Vargas (CSB2), Estrella (CSA3) y Banano (CSA2). Dicha reducción se relaciona con el posible efecto que tuvo el terremoto del año 1991, que modificó la topografía costera y expuso la plataforma marina entre 0.3 m y 1.85 m (^{Quesada, 2016}). Este afloramiento de sustrato rocoso conformó una barrera natural que mitigó el impacto del nivel de mar en sectores puntuales al sur de la provincia (^{Amador et al., 1994}; ^{Denyer et al., 1994}). Debido a lo anterior, se logró evidenciar como se redujo la erosión en playas, principalmente en aquellas donde emergió la plataforma rocosa, con la consecuente progradación en áreas cercanas a las desembocaduras de ríos.

Otro factor importante de estas playas tiene que ver con el origen de la partícula y su tamaño, pues determinan la vulnerabilidad que tienen éstas a ser erosionadas, por lo cual son fundamentales en la comprensión del transporte de sedimentos (^{Cortés et al., 1998}; ^{Salazar et al., 2004}). Puerto Vargas, Cocles, Punta Uva y Manzanillo tienen sedimentos provenientes de la desintegración de organismos marinos con estructuras duras de carbonato de calcio y de conformación fina en su mayoría que, al ser arrastrados por la acción del oleaje se pueden acumular en las cuencas submarinas de la barrera arrecifal a profundidades mayores a los 30 m, en donde ocurren procesos de litificación que conlleva a pérdida de arena (^{López et al., 2017}; ^{Medvedev y Juanes, 1981}).

En las playas del CS, el tamaño en general de las partículas muestra una tendencia a la pérdida de sedimentos finos de la duna (^{Piedra et al., 2021}), lo cual limita la capacidad de asimilación que tiene el sistema sedimentario ante el arrastre ocasionado por el oleaje de corto período y alta energía. Según ^{Nguyen et al. (2016)}, al variar el tamaño de la partícula a tamaños mayores de 1 mm, se incrementa la vulnerabilidad a la erosión, esto como consecuencia de la excavación que provoca la movilización de partículas gruesas ubicadas en la postplaya y punto de inflexión. En este proceso juega un rol de importancia la dirección del oleaje, el cual se ha determinado que sigue un rumbo noreste, y con una relación proporcional entre la altura y el período que se incrementa considerablemente en los últimos y primeros meses del año debido a la acción de los vientos alisios (^{Lizano, 2007}).

Gandoca-Manzanillo (REGAMA). En playas como Gandoca (CSB7) y Manzanillo (CSB5), denotan un patrón de tipo semi-reflexivo, en donde el ingreso de olas rompe directo en la barrera de vegetación ubicada en la duna, lo que contribuye con la extracción de importante cantidad de sedimento, el cual posiblemente se deposita en zonas profundas del lecho marino (^{López et al., 2017}; ^{Ordaz et al., 2016}; ^{Ricaurte-Villota et al., 2017}). Este aspecto es respaldado por ^{Bergoeing (2017)}, quien destaca que las isobatas que ingresan en dirección noroeste-sureste y hasta la laguna de Chiriquí, en Panamá, se acercan progresivamente al continente, lo que provoca que la plataforma continental sea angosta en dichos sectores, lo que avorece que el sedimento recorra menor distancia para ser depositado en áreas profundas.

Esta particularidad de la batimetría propicia variaciones de la profundidad en distancias cortas (^{Brenes-Rodríguez y Benavides-Morera, 2015}), favoreciendo la acción de oleaje sobre el arrastre de sedimentos que ocurre en distintas secciones de la playa (^{López et al., 2017}; ^{Ricaurte-Villota et al., 2017}). Esta condición favorece el desplazamiento de la línea de costa y pérdida de territorio continental por inundación de los ecosistemas continentales, y el posible cambio de la estructura y composición de la barrera de vegetación de la duna y posduna, de ahí que en lugares como Gandoca (CSB7, CSB6), Manzanillo (CSB5), Cocles (CBS4) y Puerto Vargas (CBS2), es recurrente la caída de árboles y palmeras de altura superior a los 10 m. Esta sería una consecuencia directa del bajo aporte de sedimentos fluviales, la influencia de la variabilidad climática que afecta el mar Caribe, los procesos tectónicos y el aumento progresivo del nivel del mar que inciden en la movilidad de la línea de costa (^{Bergoeing, 2017}; ^{Quesada-Román y Pérez-Briceño, 2019}).

Los sectores de las playas arenosas de Gandoca (CSB7, CSB6) y Manzanillo (CSB5), son parte de los bloques dinámicos y con mayores desplazamientos documentados mediante recolección de datos in situ en la última década, siendo la sección del noreste la que denota valores altos para el HotSpot erosivo. Destaca en estas secciones la existencia de dos orígenes distintos del sedimento y la presencia de barreras coralinas, que pueden restringir el arrastre de material particulado fino, aunque su efecto de protección ha sido paulatinamente superado dada la intensidad de la erosión.

Tanto Gandoca como Manzanillo exhiben un patrón de rompimiento de olas en las áreas cercanas a la zona infralitoral (ZIL), donde se ubican los bancos de arena, los cuales tienen un rol fundamental en la mitigación de los efectos del oleaje (^{Andreeva et al., 2021}), sin embargo, estos son afectados por la intensidad con que son impactados e influenciados por las corrientes de resaca, las cuales posiblemente son las responsables de desplazar el sedimento de estas barreras hacia sitios de profundidad en la plataforma continental (Castelle y Masselink, 2022), de ahí que el desplazamiento de la costa de este tramo litoral entre 1960-2023 mantuvo una tendencia general hacia la erosión.

Dicho patrón de desplazamiento de la costa es evidenciado en las continuas variaciones del perfil, mismos que cambian el ángulo de inclinación en pocos metros, conduciendo a la playa hacia una de tipo semi-reflexivo. Al respecto, ^{Van Rijn, (2009)} señala que la variación en ambos factores afecta de manera significativa el desplazamiento de la costa, siendo erosionada con facilidad y periodos cortos de tiempo, principalmente por la acción de rompimiento que tienen olas grandes en la barra, lo que afecta el flujo de retorno del sedimento desde el lecho marino (^{Andreeva et al., 2021}). De igual manera, al aumentar el nivel del mar en el CS, la zona de rompiente de las olas posiblemente sufre variaciones y es desplazada hacia el continente, favoreciendo los procesos erosivos por encima de la acreción, lo que acrecienta el desequilibrio en el ingreso de los sedimentos a la playa.

Cocles-Puerto Viejo. En contraste con lo anterior, el desplazamiento de la línea costera en el tramo Cocles (CSB4) y Puerto Viejo (CSB3) presentó poca o nula variación en el tiempo. La tasa de erosión se clasifica como estable con baja fluctuación, no así la acreción, que es importante en las zonas cercanas a las desembocaduras. La función de barrera y aporte sedimentario son parte de las características llamativas, pues su comportamiento de mitigación es efectivo al compararlo con las áreas anteriores.

Existe una concordancia entre los resultados obtenidos por DSAS y el análisis HotSpot para los segmentos CSB4 y CSB3, ya que, aunque los procesos erosivos son menos intensos, se documentó la permanencia de áreas con desplazamiento de costa por este factor, principalmente en las áreas cercanas a Puerto Viejo, y en menor magnitud hacia el sureste, en dirección a Cocles. Al igual que Manzanillo, en estos sectores existen barreras coralinas y fondos duros que contribuyen en la mitigación, reducción de la erosión y retención del sedimento, que favorece la conformación de bancos de arena en la ZIL que actúan como una barrera natural que mitiga el efecto del oleaje alta intensidad.

En cuanto a los perfiles y el ancho, se destaca que ambos sectores tienen características similares. Estos factores antes mencionados mitigan levemente las variaciones en estas condiciones. Estas características, como las que ocurren en el tramo CSB4 y CSB3, corresponden con playas de tipo semi reflexivo o intermedio, en las cuales el oleaje rompe directamente en la ZIL y en condiciones de tormenta ingresan hasta la barrera de vegetación, lo cual puede variar la granulometría y pendiente de la playa.

Parque Nacional Cahuita-Estrella. Los tramos de costa comprendidos entre Parque Nacional Cahuita (CBS2) y Estrella (CBS1) evidenciaron variación continua con fluctuaciones mínimas en las tasas de acreción, ocasionado importantes procesos erosivos, y responsables del desplazamiento de la línea costera. Las zonas críticas corresponden con punta Cahuita, y Perezoso, donde las altas tasas de EPR sugieren que existe un claro déficit entre el volumen de sedimentos que se mueve en el transporte litoral por deriva neta respecto a los que están disponibles en la duna (^{Rangel-Buitrago et al., 2015}).

Mediante el análisis de puntos calientes se logró determinar que los sectores CBS2 y CBS1 muestran comportamiento similar al descrito por las tasas EPR para el desplazamiento de la línea costera por erosión, coincidencias que aparecen principalmente en las zonas cercanas al estero Perezoso, ubicado al noroeste del área protegida. En el caso de CBS1, las zonas cercanas al río Estrella mostraron baja significancia en el clúster, lo cual puede tener relación con el flujo de sedimentos de origen continental que llegan a la desembocadura y que son depositados en la playa por acción del oleaje. Según ^{Dipper (2022)}, el transporte litoral de sedimentos por deriva neta juega un rol de importancia para la nutrición y distribución horizontal de los materiales que se depositan, esto al facilitar el reemplazo del sedimento que se mueve en semanas, meses o años de una a otra playa, favorece de esta manera el intercambio entre ellas (^{Patsch y Griggs, 2006}).

Estrella-Limón. Este tramo de la costa Estrella-Limón (CSA1, CSA2 y CSA3) presenta variación de las tasas EPR por erosión en los que hay alternancia con la acreción. En general, corresponde a uno de los sectores estables y con variación equitativa entre ambos procesos. Los escarpes en las playas arenosas aparecen durante periodos cortos que se extienden de unas semanas a meses, sin que esto pueda ser utilizado como un indicador de la reducción en el avance del mar hacia el continente.

Adicionalmente, en este lugar puede existir una contribución significativa de los bancos de arena para reducir los oleajes intensos, que son la única protección que existe, de ahí que la morfología de las playas de este tramo costero exhibe reducida exposición a cambios bruscos en la posición de línea costera, variación de los perfiles o del oleaje incidente. Al respecto ^{Taveira-Pinto et al., (2011)}, señalan que factores como las olas y las corrientes de la zona infralitoral son los responsables directos del transporte de sedimentos y por tanto, los que producen variaciones evidentes en la morfología. De hecho, los sectores CSA1, CSA2 y CSA3 parecen mantener un balance positivo de ingreso de sedimentos, los cuales son resultado de la alimentación fluvial y terrígena en la cual intervienen los ríos Vizcaya, Banano y Bananito.

CONCLUSIÓN

En general, durante los últimos 63 años, las playas arenosas del Caribe sur experimentaron procesos erosivos fuertes en aproximadamente el 70% de la costa, caracterizados por cambios abruptos en el perfil y ancho de las mismas. La variación multitemporal de la línea costera ha ido intensificando su desplazamiento en dirección hacia el continente, posiblemente como resultado de una menor resiliencia ante los eventos climáticos, al incremento en el nivel medio del mar y al deterioro generalizado de las reservas de sedimento como la duna y los bancos de arena. La situación crítica que atraviesan estos ecosistemas de transición mar-continente provocará en el mediano plazo impactos que incrementarán la vulnerabilidad y afectación de ecosistemas costeros y comunidades humanas presentes en esta franja litoral.