Materiales y métodos

Área de estudio. El municipio de Alvarado es una región portuaria ubicada en la planicie costera del centro de Veracruz, entre los paralelos 18° 33' y 19° 6' N y los meridianos -96° 7' y -95° 30' O, con una superficie de 840.63 km² (Figura 1). Este municipio tiene una baja densidad poblacional (63 hab / km²) donde predominan el cultivo de la caña de azúcar y la ganadería extensiva (^{Domínguez-Machin & Silva-López 2020}). Lamentablemente, en las últimas décadas se ha generado una gran presión en sus ecosistemas, principalmente por el desarrollo de la economía del sector primario (siendo la ganadería una de las actividades más importantes) y por el desarrollo inmobiliario (sobre todo en la zona conurbada con el municipio de Veracruz). Por lo anterior, ha prevalecido el cambio de uso del suelo, con la conversión de la vegetación original a pastizales o a construcción urbana.

Figura 1 Área de estudio, municipio de Alvarado, Veracruz mostrando el uso del suelo y vegetación (Serie VII, ^{INEGI 2021}). 

El municipio cuenta con un complejo sistema lagunar provocado por la desembocadura del río Papaloapan y el río Blanco al golfo de México, formando las lagunas Camaronera, Buen País, Alvarado y Tlalixcoyan, junto a muchas otras de menor tamaño. El municipio se encuentra en una de las regiones más biodiversas de Veracruz (^{Contreras-Espinosa 2010}) y al presentar alta diversidad avifaunística, algunas zonas del municipio son parte del Área de Importancia para la Conservación de las Aves (AICA, ^{Ramírez-Bastida et al. 2000}). Además, aproximadamente el 70 % del municipio de Alvarado fue incluido como humedal de importancia internacional en el año 2004, en la Convención sobre los humedales Ramsar (^{Portilla-Ochoa et al. 2007}).

Sin embargo, la falta de planificación y regulación de las actividades productivas y extractivas en el municipio han desencadenado diversos problemas de tipo ecológico (p. ej. fragmentación del hábitat y deterioro del manglar, ^{Silva-López et al. 2012}), económico (p. ej. pesquerías de baja productividad y falta de valoración de los recursos que ofrece el manglar) y social (p. ej. conflictos organizacionales entre cooperativas y el sector gubernamental, ^{Portilla-Ochoa et al. 1998}). La vegetación original del municipio está compuesta en su mayoría por manglar (R. mangle, A. germinans y L. racemosa), y en menor proporción por vegetación acuática de hidrófitas (tular y popal), vegetación de dunas costeras, selva baja inundable y bosque de encino tropical (Figura 1).

Datos y preprocesamiento de imágenes Landsat en Google Earth Engine. Se utilizaron todas las imágenes de reflectancia de superficie disponibles de Landsat 5, 7 y 8 del año 1985 al 2021. Las imágenes se corrigieron utilizando el Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) (^{Schmidt et al. 2013}), que es generado por Google y proporcionado por el Servicio Geológico de Estados Unidos (USGS, por sus siglas en inglés). Además, se utilizó la ruta de la huella de Landsat 024 fila 047 y el sistema de referencia mundial (WRS) -2, para mapear la dinámica de la cobertura del manglar.

Las imágenes del satélite Landsat 5 provienen del sensor TM con 7 bandas, las del satélite Landsat 7 provienen del sensor ETM que añade un canal pancromático (0.5 - 0.9 µ) a las bandas ya disponibles en el TM; por último, las imágenes del satélite Landsat 8 provienen del sensor OLI con 9 bandas. Las bandas de los satélites utilizados tienen una resolución de 30 m, por lo que se homogenizaron con las bandas de los satélites Landsat 7 y 8, a través de la función de armonización propuesta por ^{Roy et al. (2016)} que se encuentra en la implementación del algoritmo LandTrendr en GEE (^{Kennedy et al. 2018}).

Para disminuir el error por la nubosidad, en la serie de tiempo larga de Landsat, se utilizó la banda “pixel_qa” que Landsat proporciona, la cual es una banda de puntuación de nubosidad que clasifica los píxeles según su nubosidad relativa. Para crear esta banda se combinan la temperatura, el brillo y el Índice Diferencial Normalizado de Nieve (NDSI, por sus siglas en inglés), con el fin de calcular una puntuación de nubosidad que va de 0 a 100 (^{Vermote et al. 2016}). Posteriormente, ordenamos las imágenes, usando la función “sort('CLOUD_COVER', true)” incorporada en GEE, para integrar al análisis las imágenes que presentaron una puntuación de nubosidad menor a 20, lo cual elimina el error por nubosidad en las imágenes Landsat. Por último, se hizo un recorte de las imágenes usando la máscara de cobertura del manglar del municipio de Alvarado obtenido del portal de la CONABIO, el cual tiene una extensión de 234 km².

Trayectoria interanual del NDVI. Se redujeron a una sola imagen anual, todas las imágenes limpias de nubes disponibles a lo largo de un año. Para esto se hizo una colección anual de imágenes de reflectancia superficial Landsat mediante la función “medoidMosaic”, del algoritmo LandTrendr en GEE (^{Kennedy et al. 2018}). Posteriormente, se obtuvo el NDVI de cada imagen anual utilizando la siguiente fórmula:

donde R _nir es la reflectancia calibrada de la parte superior de la atmósfera de la banda del infrarrojo cercano y R _r es la reflectancia calibrada de la parte superior de la atmósfera de la banda roja. Esta fórmula refleja la diferencia óptica en la vegetación entre las bandas visible e infrarroja cercana y el fondo del suelo. Los valores del NDVI varían entre -1 y 1 para diferentes tipos de cobertura terrestre (^{Running 1990}).

Por último, se creó un compuesto de imágenes anuales del NDVI, para calcular la Trayectoria interanual del NDVI en GEE. Los resultados de GEE se exportaron al software R v. 4.1.1 (^{R Core Team 2021}) y se calculó el promedio del NDVI de los pixeles de toda el área de estudio. Se calculó la tasa de incremento del NDVI por año, realizando una regresión lineal entre el NDVI como variable dependiente y el tiempo como variable independiente. Se identificaron los cambios en la dinámica de la cobertura del manglar, utilizando los años donde hubo una modificación en la trayectoria anual del NDVI en la serie de tiempo larga. Posteriormente, se calculó el promedio del NDVI por década. Finalmente, se identificaron los cambios espaciales en la cobertura del manglar, considerando su dinámica y los datos sintéticos de la serie de tiempo del NDVI de Landsat.

Reclasificación del NDVI. Con el fin de evaluar la dinámica de la cobertura del manglar en el área de estudio, se reclasificaron todas las imágenes del NDVI, usando un umbral de corte de 0.6. Este umbral fue elegido debido a su buen desempeño para identificar píxeles con vegetación (p. ej. bosques y vegetación arbolada densa) y para discriminar píxeles con espacios abiertos y sin vegetación densa (^{Wan et al. 2019}). Todos aquellos pixeles con un valor menor a 0.6 se clasificaron como “Sin manglar” y todos aquellos con un valor igual o mayor a 0.6 se clasificaron como “Con manglar”.

Usamos la fórmula (2) propuesta por ^{Puyravaud (2003)} para obtener la tasa anual de cambio forestal:

donde A ₁ y A ₂ son el área cubierta por bosque natural en el momento t ₁ y t ₂ , respectivamente, y r es la tasa anual de cambio forestal.

Para conocer el cambio en la cobertura del manglar durante el periodo estudiado, se sumó cada imagen anual del NDVI reclasificada. Si en un pixel hubo un cambio en la clasificación “Con manglar” a “Sin manglar” se le asignó el valor de 0 en ese año. Al final se obtuvo un mapa con el número de años en que cada pixel no ha cambiado (edad del manglar), con valores que van desde 0 (sin manglar) hasta 37 (manglar sin ningún cambio en el periodo estudiado).

Variables de altura del dosel y biomasa. Se utilizaron los datos de altura del dosel forestal global, con resolución espacial de 30 m (^{Potapov et al. 2021}). Esta variable fue creada mediante la integración de las mediciones de la estructura forestal, usando datos de imágenes Lidar de la Investigación de Dinámica de Ecosistemas Globales (GEDI, por sus siglas en inglés) y las LTS. El instrumento Lidar de GEDI, a bordo de la Estación Espacial Internacional, ha estado recopilando datos únicos sobre la estructura de la vegetación desde abril de 2019. El mapa de altura forestal global fue validado obteniendo buenas métricas (error cuadrático medio = 6.6 m, error absoluto medio = 4.45 m, R ² = 0.62). Este mapa, creado con imágenes ópticas de series temporales, permite el análisis histórico multi-decenal y el seguimiento operativo de la altura del bosque y su dinámica.

Por otro lado, se utilizaron los datos de las reservas de carbono sobre la superficie de los bosques en México, a una resolución de 30 m (^{Cartus et al. 2014}). Entre 2004 y 2007, la Comisión Nacional Forestal de México (CONAFOR) estableció una red de ~ 26,000 parcelas de inventario permanente en el marco de su programa de inventario nacional, el Inventario Nacional Forestal y de Suelos (INFyS). Con esa información, generaron un mapa continuo del carbono almacenado en la biomasa viva sobre el suelo, mediante un modelo entrenado con los datos de campo del INFyS. La validación usando un conjunto independiente de parcelas del INFyS resultó en un coeficiente de determinación de 0.5 con un error cuadrático medio de 14 tC / ha, en el caso de terreno plano.

Cálculo de la calidad de la cobertura del manglar. Las variables edad del bosque, altura del dosel y reservas de carbono se estandarizaron a la unidad utilizando la fórmula:

donde x es el valor real y x’ es el valor normalizado. La suma de las variables estandarizadas fue normalizada, usando la fórmula 3, para crear un índice llamado calidad de la cobertura del manglar. Este índice varía de 0 a 1, donde 1 representa manglar maduro (dosel alto y alta cantidad de biomasa) y 0 representa el extremo contrario, es decir, la carencia de manglar. Consideramos un valor del índice de 0.731 como umbral para identificar aquellos pixeles con mayor calidad de cobertura de manglar y que podrían representar áreas de alta prioridad para su conservación. Se presentan los resultados geográficamente explícitos.

Resultados

El promedio anual del NDVI manglar del municipio de Alvarado ha incrementado linealmente durante el periodo de estudio (β = 0.004, P < 0.001, R ² _aju = 0.990). Hay una caída leve en el NDVI promedio entre los años 1992 y 1993 y los años 2011 y 2012. El promedio del NDVI por década ha ido en aumento en el área de estudio (β = 0.007, P < 0.001, R ² _aju = 0.981, Tabla 1), con una tasa de incremento forestal de 1.1 % durante todo el periodo estudiado (Tabla 1). Sin embargo, existen zonas donde ha habido pérdida de cobertura del manglar recientemente, sobre todo en la zona norte (en la conurbación con el municipio de Veracruz, en la Laguna Mandinga grande y la Laguna Camaronera, Figura 2). En todo el municipio predomina la recuperación de la cobertura del manglar durante el periodo estudiado (manglares jóvenes, NDVI < 0.6), se observa en los mapas como a través del tiempo aumentan los valores de NDVI sobre todo en la Laguna de Alvarado (Figura 2).

Sources: Esri, USGS, NOAA.

Sources: Esri, Garmin, USGS, NPS.

Figura 2 Promedio del NDVI del manglar en el municipio de Alvarado por décadas. Los colores representan el NDVI, los valores del NDVI varían entre -0.5 (representando agua y suelo descubierto, color rojo) y 1 (representando manglar denso, color verde). 

Dinámica de la cobertura del manglar. Las áreas sin cambio en la cobertura del manglar son las que mayor superficie ocupan en el municipio, observándose un aumento de estas en el periodo estudiado (de 107 a 169 km²), estas áreas se encuentran principalmente en la parte sur de todo el municipio (Figura 3). Considerando la superficie, en segundo lugar, se encuentran las áreas sin manglar, observándose un decremento en esta categoría (de 87 a 47 km²), estas áreas se encuentran sobre todo en la parte norte y sur del municipio. Las áreas que tuvieron un aumento en la cobertura del manglar disminuyeron durante el periodo estudiado (de 30 a 14 km²), estas zonas se encuentran dispersas por todo el municipio, aunque se observan unas aglomeraciones en el norte y sur del municipio. Por último, las áreas que presentaron pérdida en la cobertura del manglar se redujeron durante las primeras décadas estudiadas (11 km² en la primera década Figura 3A, 4 km² en la segunda década Figura 3B), aumentaron entre los años 2000 y 2010 (10 km²Figura 3C) y disminuyeron en la última década (llegando a 3 km²Figura 3D). Estas áreas se concentran en el norte (especialmente en la conurbación con el municipio de Veracruz) y en el sur del municipio (Figura 3).

Sources: Esri, USGS, NOAA.

Figura 3 Dinámica de la cobertura del manglar en el municipio de Alvarado entre la década de 1980 y 1990 (A), 1990 y 2000 (B), 2000 y 2010 (C) y 2010 y 2020 (D). El color rojo representa pérdida en la cobertura de manglar, el azul ganancia en la cobertura del manglar y el verde representa la permanencia en la cobertura de manglar. 

El 50 % de la cobertura del manglar del área de estudio (118 km²) está conformada por zonas en donde dicha cobertura no ha sido perturbada desde el inicio del periodo estudiado, acumulando entre 25 y 36 años sin cambio (color verde Figura 4). Sin embargo, el 36 % de la cobertura del manglar del municipio de Alvarado (85 km²) se compone de cobertura que ha sido perturbada recientemente (entre 0 y 12 años, color rojo Figura 4). El resto del municipio se compone por cobertura del manglar con una edad entre 13 y 24 años sin perturbación (14 %, 32 km², color amarillo Figura 4).

Sources: Esri, USGS, NOAA.

Figura 4 Edad de la cobertura del manglar en el municipio de Alvarado, Veracruz. 

Por otro lado, la mayoría de la cobertura del manglar en el municipio se encuentra ocupada por manglar que tiene entre 0 y 10 m de altura (53 % del área total), seguida por manglar que tiene entre 11 y 20 m de altura (37 %). En el municipio, aún existen zonas de manglar con una altura del dosel entre 21 y 29 m (10 %,) ubicadas principalmente al sur del municipio (Figura 5A).

Sources: Esri, USGS, NOAA.

Figura 5 A) Altura del dosel y B) Biomasa aérea del manglar en el municipio de Alvarado. 

La mayoría del manglar en el municipio se encuentra ocupado por zonas que contienen entre 1 y 20 tC / ha (aproximadamente 189 km²), seguido por áreas que contienen entre 21 y 40 tC /ha (45 km²). Mientras que, solo el 0.2 % del área de estudio (0.4 km²) almacena una alta cantidad de biomasa que representa entre 41 y 58 tC / ha (Figura 5B), estas zonas se ubican en el sur del municipio.

Con el índice de calidad de la cobertura del manglar, se lograron identificar áreas que cuentan con manglar maduro, dosel alto y gran cantidad de biomasa, estas se ubican en el sur del municipio. En general, el municipio tiene una cobertura del manglar con baja calidad, la cual está constituida por árboles en desarrollo de bajo dosel y poca biomasa (Figura 6A). Por último, los sitios prioritarios para la conservación de la cobertura del manglar en el municipio se encuentran ubicados al sur, representando un 3 % de la superficie total del municipio y 9 % del manglar (22 km², Figura 6B). Estos resultados pueden ser visualizados de manera interactiva en rpubs.com/jgomezd/veg_alv.

Sources: Esri, USGS, NOAA.

Figura 6 A) Índice de calidad de la cobertura de manglar utilizando edad del manglar, altura del dosel y biomasa; y B) Sitios prioritarios para la conservación del manglar en el municipio de Alvarado, Veracruz. 

Discusión

Dinámica de la cobertura del manglar. Las variaciones interanuales del NDVI dentro del área de estudio son evidentes. El 50 % de la cobertura del manglar en el municipio es joven (entre 0 y 24 años), como consecuencia de la transición forestal. Aunque múltiples factores pueden explicar el proceso, en particular el abandono de las tierras agrícolas, por las altas tasas de emigración a Estados Unidos de América, podría ser uno de los más importantes (^{Tuirán 2002}, ^{López-Barrera et al. 2014}). Asimismo, dentro del municipio de Alvarado se ha registrado que un considerable porcentaje de agricultores ha abandonado la actividad agrícola, debido a la implementación de los acuerdos de libre comercio del TLCAN durante la última década (^{López-Barrera et al. 2014}). Existe una relación entre los procesos de abandono, crecimiento de bosque secundario y recuperación del manglar, que se refleja en la disminución de la tasa de pérdida de manglar. De manera similar, la recuperación de la vegetación se ha observado principalmente en bosques templados y de mayor altitud, en el centro del estado de Veracruz (^{Gómez-Díaz et al. 2018}), por fenómenos sociales similares a los que ocurren en el área de estudio.

El marco de transición forestal considera el abandono de tierras como una oportunidad para la conservación de la biodiversidad (^{Rudel et al. 2005}). Según ^{Hosonuma et al. (2012)}, México se encuentra en la Fase 3 del modelo (transición tardía), donde por el momento no se experimenta una renovación forestal neta. Una gran proporción de los bosques, en los países que se encuentran en la última fase de transición, siguen degradados (^{Hosonuma et al. 2012}), lo que ocurre en Alvarado.

Algunos impulsores de la degradación del manglar que se registran en el área son la extracción de madera, la tala y los incendios descontrolados, los cuales son los principales en América Latina y México (^{Hosonuma et al. 2012}, ^{Kolb et al. 2013}, ^{Ellis et al. 2017}). Observamos que en el área aún se presentan procesos de deforestación en el manglar, debido a la expansión de la ganadería con la consecuente creación de pastizales. Además, ^{Gómez-Díaz et al. (2018)} encontraron que las tierras bajas del centro de Veracruz tienen la mayor proporción de bosque convertido. Lo cual se corrobora en nuestra área de estudio, siendo estas zonas las más utilizadas para la ganadería, el desarrollo inmobiliario y la agricultura.

Validación y comparación con otros estudios. Los resultados de la pérdida y ganancia de la cobertura del manglar en Alvarado son precisos, cuando se comparan con imágenes históricas de alta resolución. Estos resultados demuestran que las imágenes de teledetección a largo plazo pueden detectar la dinámica del uso del suelo, en áreas con intensa actividad humana. ^{Gómez-Díaz et al. (2018)} analizaron el cambio en la cobertura forestal en el centro de Veracruz, en el período de 1993 al 2014 utilizando imágenes Landsat, encontrando que el área cubierta por bosques remanentes disminuyó un 57 %. También mostraron que la tasa de cambio anual, de la cubierta forestal neta de 1993 al 2000, fue de -0.44 %, desde entonces, la cubierta forestal aumentó a una tasa del 0.11 % entre los años 2000 y 2014. Estos resultados son similares a lo observado en nuestra tendencia anual, donde hubo un decremento en el NDVI promedio del municipio entre los años 1992 y 1993, mostrando una gradual recuperación desde entonces.

De acuerdo con la ^{CONABIO (2016)}, en un análisis realizado entre el 2010 y 2015, existen zonas de pérdida y ganancia de manglar al sur y al norte de la laguna de Alvarado y la laguna Pajarillos. Estos patrones son semejantes a los encontrados en este estudio, donde incluso el índice de calidad de la cobertura del manglar muestra que esta área ha sufrido una perturbación. Asimismo, ^{Versañez-Vences (2018)} documentó patrones similares en la misma zona de estudio durante el periodo de 1985 al 2015.

También se ha demostrado la eficacia de utilizar series temporales largas con imágenes Landsat y GEE en humedales costeros de China (2,798 imágenes ETM+/OLI) (^{Wang et al. 2020}). La dinámica interanual de dichos humedales, durante las últimas tres décadas, siguió un patrón similar a lo encontrado en el municipio de Alvarado: un período estable durante los años 1986 -1992, un periodo creciente durante 1993 - 2001 y un período decreciente durante 2002 - 2016. Al igual que en nuestro estudio, las actividades humanas desempeñaron un papel importante en el cambio de área, la erosión y la expansión de las planicies intermareales en China. Los impactos antrópicos comenzaron a pequeña escala durante 1986 y 1990. Sin embargo, en décadas recientes el aumento en la población y la presión para cubrir la demanda de recursos, han disminuido la cobertura de las planicies intermareales. Este parece ser un patrón común al expandirse el crecimiento urbano en zonas planas y costeras (^{Cissell et al. 2018}, ^{Pinheiro 2020}, ^{Wang et al. 2020}).

^{Pinheiro (2020)} cuantificó los cambios en la superficie cubierta por los manglares, en los años 2005, 2016 y 2019, a través de imágenes de satélite Landsat y analizó el impacto ambiental de las actividades antropogénicas. Este autor determinó que los bosques de manglar se redujeron aproximadamente en un 18 % entre 2005 - 2016, siendo sustituidos en su mayoría por cubiertas “Urbanas”, mientras que entre 2016 y 2019 hubo un incremento del 2.08 % de estos ecosistemas. Él atribuyó como causa principal de decrecimiento, las actividades antrópicas asociadas a proyectos turísticos (^{Pinheiro 2020}). En el caso del municipio de Alvarado, la reducción en la cobertura del manglar se debe al aumento en las actividades antropogénicas tales como la expansión urbana.

^{Cissell et al. (2018)} realizaron la clasificación no supervisada de datos satelitales Landsat, de 1999 y 2016, para mapear y cuantificar los cambios en la cubierta forestal de manglares, a lo largo de la costa de Campeche. Ellos encontraron una disminución del 14 % en el área de manglares y un aumento del 19 % en la extensión urbana de la ciudad de Campeche durante este periodo. La disminución de los manglares se concentró en las franjas costeras. Gran parte de la degradación se debió a la alteración de los flujos de agua de los canales de los manglares, debido a los proyectos de construcción de carreteras, así como al aumento de la presión indirecta de la población en expansión de Campeche. El corredor costero de manglares de Campeche está protegido como Reserva de la Biosfera por la UNESCO, sin embargo, su proximidad al centro urbano y el crecimiento de la ciudad lo hace vulnerable a los impulsores directos e indirectos de la degradación de los manglares, lo cual también ocurre en el municipio de Alvarado, donde la pérdida de cobertura del manglar se concentró en las zonas conurbadas.

Implementación de series de tiempo largas en GEE. La composición promedio del NDVI integrada con los datos de imagen disponibles en un año determinado, puede simplificar el procesamiento de datos, lo que resulta en un algoritmo más eficiente. Aunque el procedimiento de máscara de nube se utilizó antes del análisis del NDVI, algunas nubes residuales y neblina en las imágenes podrían afectar la precisión de la teledetección. Sin embargo, como demostramos en este estudio, reducir las imágenes disponibles de un año a una sola imagen anual puede mejorar la precisión y reducir el impacto de las nubes (^{Huang et al. 2017}). Otros factores que podrían afectar los cambios reportados en el NDVI son la variabilidad en los valores de reflectancia, debido a la heterogeneidad del paisaje y por la influencia de la estacionalidad climática. Esto se pudo controlar debido a que la zona de estudio es homogénea en cuanto al tipo de vegetación (manglar) y la cobertura dominante (pastizal). Además, la mayoría de las imágenes utilizadas en este estudio fueron de la estación seca (debido a su baja cobertura de nubes). También cabe mencionar que en este estudio solo utilizamos el NDVI, el cual es recomendado en los análisis de detección de cambios de tendencia a largo plazo (^{Forkel et al. 2013}). En un futuro próximo se podría considerar una comparación entre diferentes índices.

La computación en la nube, con los datos de Landsat archivados en GEE, tiene muchas ventajas en el mapeo a gran escala y de series de tiempo largas, como el monitoreo de la tala de vegetación leñosa y la detección de la extensión de los asentamientos humanos (^{Huang et al. 2017}). Asimismo, es fácil integrar series de tiempo del NDVI estacionales para mejorar la precisión de la clasificación de la vegetación, disminuyendo el error entre las sombras de las nubes y el agua, además de reducir la confusión entre los diferentes tipos de vegetación.

A medida que se integren más algoritmos de aprendizaje automático y datos de teledetección a la plataforma de GEE y al catálogo de datos, se espera que la extracción de información de teledetección se simplifique aún más. Mientras tanto, la falta de funcionalidad requiere la descarga de datos para el procesamiento local adicional y la carga de datos a GEE. Las limitaciones de ancho de banda pueden hacer que este proceso sea lento y difícil de manejar. La capacidad para materializar o leer productos de datos intermedios, hacia o desde el almacenamiento basado en la nube, podría aumentar la capacidad de usar otras herramientas informáticas basadas en la nube (por ejemplo, Compute Engine) en lugar de Earth Engine, disminuyendo la necesidad de procesamiento local y la transferencia de datos intensiva en ancho de banda.

Nuestros resultados resaltan la importancia de la teledetección para un monitoreo ambiental consistente y verdadero, en comparación con las estadísticas publicadas oficialmente. Las imágenes de mayor resolución, como los datos de Sentinel-2 y ZY-3, pueden reducir la confusión dentro de la vegetación y mejorar la precisión de la clasificación (^{Huang et al. 2017}). Por ejemplo, la mayor parte de la ganancia en cobertura del manglar se detectó al sur del municipio, precisamente en áreas donde se han desarrollado programas de reforestación (^{Lucio-Palacio et al. 2014}).

Implicaciones para la conservación. El manglar de este municipio cumple múltiples beneficios; como la estabilización de los ecosistemas de humedales, la prevención de la erosión y la protección a zonas pobladas adyacentes (absorbiendo los impactos de las marejadas ciclónicas durante fenómenos meteorológicos extremos, como en los huracanes). Por lo que de acuerdo con ^{Pérez-Vega & Ortiz-Pérez (2002)}, la pérdida de la cobertura del manglar incrementa los gastos máximos hidrométricos, en las subcuencas del margen derecho del río Papaloapan. Además, de acuerdo con ^{Brena-Zepeda et al. (2012)}, la pérdida de manglar maduro representa una amenaza a la sustentabilidad y a la biodiversidad del área estudiada. Lamentablemente, la superficie destinada para agricultura y ganadería ha ocupado terrenos que antes fueron áreas con manglares, propiciando la transformación de la vegetación original a pastizal inducido (^{Brena-Zepeda et al. 2012}).

La mayoría de los sitios prioritarios para la conservación del manglar y la biodiversidad se encuentran fragmentados, sobre todo en la zona norte del municipio, mientras que en la zona sur se encuentra un continuo de cobertura del manglar con buena conectividad. Estos sitios se han conservado principalmente por ser de difícil acceso, la dificultad para establecer pastizales o ganado (por las características ecológicas de este tipo de vegetación) y por la protección de las especies de manglar mediante las leyes mexicanas (^{Silva-López et al. 2022}). Lo anterior, refuerza la importancia de preservar estos últimos relictos de manglar original con buena calidad mediante el decreto de un Área Natural Protegida (ANP).

El 70 % de la extensión del municipio de Alvarado se encuentra protegido por el Convenio de Ramsar, sin embargo, la vegetación en muchas zonas se encuentra altamente transformada. El humedal de Alvarado, al ser un sitio Ramsar, tiene reconocimiento internacional, lamentablemente, la conservación de esta área no se encuentra regulada en la práctica por la leyes nacionales, estatales o municipales, por lo que existe un vacío legal que ha permitido el cambio de uso de suelo, con la consiguiente deforestación y fragmentación de la vegetación original. Por ejemplo, de acuerdo con ^{Vázquez-Lule (2008)}, la laguna de Alvarado tuvo un cambio en la cobertura del manglar pasando de 195 km² en 1973 a 138 km² en 2006 con una tasa de deforestación anual del 1.04 %. Por lo anterior, es importante proteger mediante otro mecanismo (ANP de nivel federal o estatal) las áreas con cobertura de manglar más extensas que aún quedan en el municipio, debido a su importancia para la conservación de la diversidad de flora y fauna, en particular, porque están amenazadas por la expansión urbana y agrícola.

Muchos de los sitios prioritarios para la conservación tienen una valiosa integración de atributos ecológicos, como el aislamiento producido por el cordón de dunas costeras que forman una pequeña cuenca en la que quedan inmersos los manglares, el encontrarse en una zona de transición con las selvas bajas de la región de Los Tuxtlas, la vegetación acuática característica de la cuenca del Papaloapan, y la riqueza de especies de su fauna de peces, anfibios, reptiles, aves y mamíferos, conformando alrededor de 400 taxones (^{León et al. 2007}). Lamentablemente, nuestros análisis muestran que esta zona ha sido altamente impactada por el cambio de uso de suelo, por lo que urgen programas de restauración y conservación en esta área.