El incremento del nivel del mar es un evento extremo que impactará gravemente las zonas costeras bajas, como es el caso de Tabasco. A nivel nacional la costa tabasqueña es de los sitios más vulnerables ante tal evento, por lo que es necesario estimar localmente el grado de su vulnerabilidad. Por ello, en este estudio se estimó esa vulnerabilidad mediante la aplicación de la metodología del Indice de Vulnerabilidad Costera utilizando Sistemas de Información Geográfica. Este modelo se ajustó a las condiciones locales de la costa tabasqueña, dado que los datos que exigen las variables que integran la metodología (olas, mareas, nivel del mar, pendiente costera, tasas de erosión y geomorfología) estaban disponibles para el área de estudio, y así contar con la oportunidad de encontrar datos locales del grado de vulnerabilidad de la costa. De acuerdo con los resultados encontrados, tras realizar un álgebra de mapas, los sectores con mayor vulnerabilidad se encuentran justamente frente a los sistemas lagunares costeros más importantes del estado, Carmen-Pajonal Machona y Mecoacán, este último con un mayor grado de vulnerabilidad. Es importante destacar que dicha área se encuentra dentro de la zona de influencia del puerto de Dos Bocas, donde la actividad antropogénica es muy notable. Los datos encontrados serán útiles para la toma de decisiones para el manejo del espacio costero de Tabasco.
Sea level rise is one of the most serious events that will impact low-lying lands, as is the case of most of Tabasco State. Historically, the State of Tabasco has been repeatedly impacted by extreme floods, the most recent one occurring in 2007. However, recent studies have shown that coastal erosion is the effect that most directly has impacted the Tabasco's coastline, as this has even modified soil strata; this is also related to extreme hydrometeorological events associated with environmental changes and changes in the salinity gradient off the coast. In such a situation, future changes in the coastline are almost certain. Tabasco's coastline has been recognized as one of the most vulnerable zones in the country since Mexico's first national communication to the unfccc in 1997. Therefore, it is important to evaluate the vulnerability of this zone. The purpose of this study was to estimate the vulnerability of the Tabasco's coastline by applying the coastal vulnerability index method (ivm) using a geographic information system (gis). This method has been successfully applied in several different places around the world including Canada, the United State, Spain and Indonesia. This model is suitable for the local conditions of Tabasco coast, as the input variables it requires (including waves, tides, sea level, coastal slope, erosion rates and geomorphology) are available for the study area, thus allowing the possibility of estimating the coast's vulnerability based on local data. Results from map algebra operations showed that the zones of very high or high vulnerability encompass a six-kilometer stretch around the Sánchez Magallanes community, near the del Carmen lagoon in the municipality of H. Cárdenas, Tabasco. This is due to the high-waves regime and other conditions associated to the coastal dune geomorphology as well as the unconsolidated fine sediments prevailing therein. Other high vulnerability zones are found just in front of the major lagoon system in the State of Tabasco, the Carmen-Pajonal-Machona and Mecoacan lagoons; being the last one the most vulnerable of all. It is worth pointing out that this zone is located within the influence zone of the Dos Bocas port where an intense anthropogenic activity occurs. It is also important to point out that this study is one of the first approaches to the estimation of coastal vulnerability in Mexico. The implementation and application of this model of coastal vulnerability evaluation are significant given the spatial scale of the study and that this is the first time that these variables have been collected for the state of Tabasco. Our results might inform decision-making processes on the proper management of the Tabasco coastline, thus benefiting the local communities.
Las zonas costeras están severamente amenazadas por el cambio climático (CC), y el incremento del nivel del mar (INM) representa el mayor peligro para estas regiones (IPCC, 2007). Albergan ecosistemas altamente productivos (manglares, arrecifes coralinos, lagunas costeras, marismas, estuarios, entre otros) (Burningham, 2004; Bijlsma, 1996), pero también son de las regiones más pobladas. Actualmente, las áreas costeras están convertidas en grandes urbes debido a la alta concentración económica; se estima que cerca del 50% de la población mundial habita en una franja costera de 100km (Small y Nicholls, 2003), y este crecimiento poblacional es mayor entre los 0 y 10 m de altitud (McGranahan et al., 2007). Las ciudades costeras bajas, densamente pobladas y con niveles altos de industrialización, pueden enfrentar mayores impactos ante el INM (Aris y King, 2008; McGranahan et al., 2007). La magnitud de los impactos puede ser diferente de una región a otra, ya que nivel del mar cambia de un lugar a otro por condiciones y procesos locales, como la topografía, la geología, los movimientos naturales de tierra, la climatología regional y las actividades antrópicas que contribuyen a cambios en el nivel del agua, como la extracción de fluidos del subsuelo y la oferta de sedimentos que eventualmente llegan a la costa; aun así se recomienda evaluar el impacto del INM en la medida de lo posible (Kebede et al., 2012; Nicholls y Cazenave, 2010). Los impactos potenciales del INM figuran en los planes de desarrollo de gobiernos en todo el mundo, pero la preocupación es mayor para los países en vías de desarrollo, donde ninguna o muy pocas ciudades costeras están preparadas para enfrentar el INM (Nicholls et al., 2008).
En México, la costa de Tabasco es considerada una de las zonas más vulnerables ante este evento por su ubicación en una extensa llanura inundable y su elevación menor a un metro (SEMARNAP, 1997). Tabasco es territorialmente una de las entidades más pequeñas, pero económicamente una de las más productivas, al ubicarse entre los diez estados con mayores aportaciones al producto interno bruto (PIB) nacional (INEGI, 2013). Con 200km de litoral costero repartido entre los municipios de Huimanguillo, Cárdenas, Paraíso y Frontera, la región costera del estado constituye el principal motor económico de la entidad, pues allí se realizan actividades económicas preponderantes, como la extracción de crudo, pesca comercial y artesanal y, en menor grado, agricultura y ganadería (Palma et al., 2007). La presencia de importante infraestructura, como el puerto de Dos Bocas, ha sido clave para el crecimiento y desarrollo económicos de la región y el país, impulsando al sector industrial, comercial y petrolero (Cato, 2009). En ese sentido, si la costa tabasqueña es afectada por el cambio climático, particularmente por el INM, tendría un impacto socioeconómico en la región y el país.
Tabasco se ha visto históricamente afectado por eventos extremos, principalmente por inundaciones ocurridas en los últimos 20 años (1995, 1999 y la más drástica en 2007 (Gama et al., 2008; IMTA-IPN-CP, 2008). También existe evidencia de erosión costera, por ejemplo, Hernández et al. (2008) reportaron que en los últimos 30 años el litoral ha sufrido retrocesos costeros en Sánchez Magallanes (-3 a -5 m/año) y en la desembocadura río San Pedro y San Pablo (-8 a -9m). Gornitz (1990) considera que las costas de alto riesgo ante el INM son aquellas con bajo relieve, sustrato erosionable, retroceso histórico de la línea costera y oleaje de alta energía; esas son condiciones que presenta la costa tabasqueña y por lo tanto constituye uno de los sitios más críticos ante la elevación del nivel del mar (Ortiz y Méndez, 1999; Magaña y Gay, 2002). Los promedios anuales de INM en la región Golfo de México oscilan entre 3.1±0.8mm yr-1 con un incremento al 2090 de 60 cm (IPCC, 2007), no obstante modelos recientes indican que el aumento al final de este siglo será mayor a un metro (Rahmstorf, 2012). Los cálculos realizados considerando un escenario de 1 m de elevación en el nivel del mar estiman inundaciones del 12% del territorio tabasqueño, lo que representa poco más de 3 000km2 (INE- UNAM-SEMARNAT, 2008). Ello representaría la modificación en la ecología de recursos de alto valor ecológicos (humedales, sistemas lagunares, entre otros) y socioeconómicos (recursos pesqueros, activos de producción petrolera, puertos, entre otros), que son especialmente importantes para la economía estatal y nacional. Asociado a ello, habitan en una franja de 5km poco más de 90 000 personas en 147 localidades, algunas ubicadas en la cota de 5 m (INEGI, 2010). Debido a lo anterior, en Tabasco existe gran preocupación gubernamental para atender la problemática del INM y las posibles afectaciones sobre el gran patrimonio socioambiental del estado (Hernández-Santana et al., 2008).
El IPCC (2001) define vulnerabilidad como “la incapacidad de las zonas costeras para afrontar la variabilidad climática y eventos extremos”. Este concepto enfatiza la idea de enfrentar el evento adverso con las fortalezas y debilidades inherentes al sistema costero; en este sentido, los componentes costeros locales ejercen un papel determinante en la vulnerabilidad física de la costa frente al INM (Pendleton, Thieler y Williams, 2010). El IPCC, organismo reconocido para atender todo lo relacionado con el tema de CC, ha determinado que las evaluaciones de vulnerabilidad sobre los sistemas naturales y humanos representan el primer paso para estimar los impactos del CC. Es decir, cualquier estudio de vulnerabilidad parte de la evaluación de impactos físicos en el medio natural (Klein y Nicholls, 1999). Según el Marco de Hyogo una acción clave para evaluar vulnerabilidad es desarrollar sistemas de indicadores junto a metodologías que permitan evaluar el impacto de posibles desastres en el medio (Kappes et al., 2012; EIRD-ONU, 2008).
Bajo este criterio, el presente trabajo intenta evaluar la vulnerabilidad costera del litoral tabasqueño considerando seis variables/indicadores conjuntos (geomorfología, tasa de erosión, pendiente costera, oleaje, rango mareal y nivel del mar), mediante la aplicación de la metodología de Indice de Vulnerabilidad Costera (IVC) (Ojeda et al., 2011; Thieler y Hammar 1999; Shaw et al., 1998; Gornitz, et al., 1992) a fin de identificar, a nivel local, zonas de mayor vulnerabilidad frente a este evento extremo.
Materiales y métodosÁrea de estudio y localización geográficaLa zona costera de Tabasco se ubica entre los 92° 28’ y los 94° 10’ de longitud Oeste, y los 17° 15’ y los 18°39’ de latitud Norte (Hernández- Santana et al., 2008; Figura 1), con un litoral costero de 200km repartido entre los municipios de Huimanguillo, Cárdenas, Paraíso y Centla, que en su conjunto albergan 147 localidades en una franja de 5km tierra adentro con altitudes que van de 0-5 m, con un aproximado de 92 500 habitantes (INEGI, 2010). El mar adyacente a su línea de costa es el Golfo de México, uno de los mares más complejos y productivos del mundo (Tunnell et al., 2004).
La costa tabasqueña es resultado de importantes procesos aluviales (West et al., 1969), originados por la dinámica fluvial ejercida por la cuenca Grijalva-Usumacinta (la más extensa y caudalosa de México) (CONAGUA, 2012; Hudson et al., 2005) y los ríos González y Tonalá. De acuerdo con estaciones mareográficas de la región, la marea es de tipo mixta con influencia diurna, el oleaje es moderado en dirección E-W con alturas máximas de 2 metros en condiciones meteorológicas normales; sin embargo, los vientos provenientes del norte aumentan considerablemente el oleaje con alturas de hasta 5 metros (SEMAR, s/f; Lankford, 1977). El litoral tabasqueño alberga recursos naturales de gran importancia ecológica, entre los que destaca la Reserva de la Biosfera Pantanos de Centla (RBPC), considerada la unidad de almacenamiento energético más importante de Norteamérica por su alta productividad primaria (Toledo, 2003). Otros ecosistemas de gran valor son las lagunas costeras Carmen-Pajonal-Machona y Mecoacán, las que al conjuntarse con las aguas interiores provenientes de la subcuenca del Grijalva establecen condiciones estuarinas, que permiten albergar, espacial y temporalmente, fauna de origen marino y dulceacuícola así como peces e invertebrados (Contreras et al., 2002; Sánchez y Barba, 2005). Esta mezcla de ambientes acuáticos constituyen áreas de importancia comercial de donde se extrae el 90% de la pesca del estado (Contreras y Castañeda, 2004).
Asociados a estos cuerpos lagunares costeros, en la parte posterior del cordón litoral de Tabasco, se encuentran comunidades de manglar de las especiesRhizophora mangle (mangle rojo), Avicennia germinan (mangle negro), Laguncularia racemosa (mangle blanco) y Conocarpus erectus (mangle botoncillo) (García y Palma, 1998; Villanueva et al., 2010), adaptadas a diferentes gradientes de salinidad e inundaciones recurrentes.
Implementación del modeloLa utilización de esta metodología se justifica por la aplicabilidad del IVC a cualquier entorno costero, dado que a diferencia del rango mareal los demás factores responsables del cambio en la costa son los mismos (Pendleton, Thieler y Williams, 2010). Partiendo de tal criterio, este modelo ha sido aplicado en diversos litorales costeros, tanto en latitudes tropicales como India, subtropicales como España y aún en regiones cercanas al Ártico como Estados Unidos y Canadá, cada uno con diferentes consideraciones y particularidades (Cuadro 1). El IVC está basado en variables mayormente cuantitativas; no obstante, es necesario aclarar que solo indican la vulnerabilidad relativa de la costa ante el INM. En este contexto, el concepto de vulnerabilidad adoptado, al cual se hará referencia en adelante, queda definido como el grado (o la magnitud) del daño que puede experimentar la costa tabasqueña por el INM en función de su exposición y su capacidad de adaptación.
Variables consideradas en cada uno de los países evaluados.
| Variables | USA | Canadá | USA | España |
|---|---|---|---|---|
| Gornitz et al., 1992 | Shaw et al., (1998) | Thieler y Hammar (1999) | Ojeda et al., (2009) | |
| Relieve (altitud del terreno) | ✓ | ✓ | ||
| Tipo de roca (resistencia a erosión)* | ✓ | ✓ | ||
| Forma costera* | ✓ | ✓ | ✓1 | ✓1 |
| Pendiente costera | ✓ | ✓ | ✓ | ✓2 |
| Cambio relativo del nivel del mar | ✓ | ✓ | ✓ | ✓ |
| Cambio en la línea de costa (erosión-acreción) | ✓ | ✓ | ✓ | ✓ |
| Rango mareal | ✓ | ✓ | ✓ | ✓ |
| Altura de oleaje medio significante | ✓ | ✓ | ✓ | ✓ |
Incluye tipo de roca (litología)1, sustituida por Índice Topograico 2, *variables descriptivas.
La primera fase del modelo se centra en la cuantificación de las variables que contribuyen a la evolución de la costa en el área determinada. En este caso, a diferencia de los estudios anteriormente referidos (por ej., EE.UU., Canadá), aplicamos el modelo haciendo consideraciones parecidas a las del estudio de Ojeda et al. (2009), utilizando solo seis variables que se separan en dos grupos (Cuadro 2):
- a)
Geológicas/geomorfológicas. Consideradas en función de la resistencia que ofrece cada paisaje costero de la costa frente a la erosión y la susceptibilidad a los procesos de inundación.
- b)
Físicas/hidrodinámicas. Consideradas por su contribución a los procesos de erosión e inundación.
Variables consideradas para la estimación de vulnerabilidad en Tabasco.
| Tipo de variable | Variable | Indicador |
|---|---|---|
| Geológicas Geomorfológicas | Geomorfología/Geología | Resistencia a la erosión |
| Erosión | Cambio de la línea de costa | |
| Pendiente costera | Susceptibilidad a inundación | |
| Hidrodinámicas | Cambio nivel del mar (cm) | Peligro a inundación |
| Oleaje medio significante (m) | Peligro a inundación y erosión | |
| Rango mareal (m) | Peligro a inundación y erosión |
Las variables consideradas se centran en las características físicas de la costa y exposición a diversos procesos marinos (oleaje, mareas y cambio relativo en el nivel del mar) (Ojeda et al., 2011) por lo que, en sentido práctico, el mecanismo del modelo implica contraponer las variables constituidas como forzamientos de la vulnerabilidad (variables hidrodinámicas) y la respuesta de la costa mediante sus condiciones terrestres (variables geológicas/ geomorfológicas (Cuadro 2).
Cada una de las seis variables consideradas se clasifican individualmente en intervalos de vulnerabilidad que van desde 1 (muy baja vulnerabilidad) a 5 (muy alta vulnerabilidad). Esta simple clasificación permite la integración en la siguiente ecuación que finalmente proporciona el IVC.
Cálculo de las variables
La corrida del modelo implicó un sin número de procesos de análisis geográficos con softwares especializados (Arc GIS, Quantum GIS, Arc View), dentro de los que destacan:
- 1.
La segmentación del litoral costero en rectángulos perpendiculares a la línea de costa de 500 m de ancho aproximadamente y delimitado tierra adentro por las curvas de nivel 5 m (donde la curva de nivel no apareció en las inmediaciones del litoral costero, la medida estándar fue de 5km, suponiendo que es un parteaguas entre el mar y tierra adentro, como lo indica la metodología de IVC). Este polígono se utilizó como medida estándar para proyectar los datos de cada variable sobre el litoral costero tabasqueño.
- 2.
La homogeneización de todas las variables, a fin de realizar eficientemente el análisis geográfico dentro del SIG, se realizó mediante el proceso de rasterización de todos los archivos vectoriales de cada variable, en el módulo Conversion Tools de Arc GIS. Las fuentes de las variables se muestran en el Cuadro 3.
Cuadro 3.Base de datos (geodatabase) generada para correr el ivc en Tabasco.
Variable Formato Fuente/Autor Geología/Geomorfología Vectorial Ramos et al., 2014 Erosión (m/año) Vectorial Barrios (2012), laige- ecosur (2012) Pendiente costera Raster vectorial inegi, 2010 Altura de ola significante (m) Vectorial (Appendini et al., 2013), Appendini et al., 2014 the waves generated by the Caribbean low-level jet (cllj) Tasa de nivel del mar (cm) Raster unam- 2uniatmos (2012) Rango mareal (m) Vectorial semar (2014) 1Laboratorio de Ingeniería y Procesos Costeros
Fuente: Elaboración propia. - 3.
La corrida general del modelo IVC se realizó mediante el proceso Algebra de Mapas en el módulo Spatial Analysis de Arc GIS.
Los datos utilizados para cada una de las variables, provenientes de diversas fuentes se insertaron en una base de datos espacial (geodatabase). Los insumos necesarios y el procedimiento de cálculo se resumen en el Cuadro 4.
Resumen del procedimiento de cálculo de cada variable.
| Variable | Procedimiento de cálculo |
|---|---|
| Geología Geomorfología | Clasificación de acuerdo al paisaje muestreado en cada sector (500 metros). La clasificación geomorfológica utilizada se elaboró mediante recorridos de campo y se consideraron características geológicas de la zona. |
| Erosión (m/año) | Comparación de las líneas de costa (1995-2012). La diferencia resultante entre ambos polígonos se subdivide en tramos de 500 m y se aplica la fórmula de ANIDE (2005).1 |
| Pendiente costera | Se extrajeron las curvas de nivel hasta 5km tierra adentro de un lidar propiedad del inegi (4 m de resolución), posteriormente se creó un mapa de pendientes y mediante una interpolación se proyectó a la línea de costa para asignarle valor a cada tramo costero. Los beneficios de la utilización de esta herramienta en el modelo se reflejan en la obtención de datos muy cercanos a la realidad. |
| Altura de ola significante (m) | Proyección de los puntos de oleaje sobre el litoral tabasqueño, posteriormente se calcularon promedios de altura de los meses (octubre-noviembre y diciembre), mediante una interpolación se proyectaron hacia la línea de costa. |
| Tasa de nivel del mar (cm) | Se sobrepuso el área de estudio sobre una imagen Geotiff obtenida del Servicio Mareográfico Nacional y se extrajeron los valores correspondientes a cada sector. |
| Rango mareal (m) | Se obtuvieron los datos de semar de cada una de las estaciones mareográficas. Se pasaron a una base de datos de Excel para calcular la diferencia entre la marea alta y la marea baja. |
La clasificación obtenida indica que existen sólo dos grados de vulnerabilidad: alto y muy alto, debido a que todo el frente litoral se encuentra ocupado por dunas costeras y cordones litorales de baja altitud. Las geoformas costeras mencionadas se han mantenido gracias a la entrada de sedimentos que son removidos por la hidrodinámica costera y acumulados en esta porción del litoral (West et al., 1969). Actualmente, estas geoformas están perdiendo la batalla contra la erosión marina debido en gran parte a la disminución en la oferta de sedimentos provocado por la construcción del sistema de presas en el río Grijalva, uno de los principales ríos de la cuenca que transportan las mayores cantidades de sedimentos hacia la costa (Syvitski et al., 2005).
Por otro lado, la remoción artificial de la duna costera para construir carreteras o viviendas en la costa (Brammer, 2014; Kokot y Chomnalez, 2012) han eliminado gran parte de este relieve que fungía como barrera y que, precisamente, estaba conformado por las zonas de dunas costeras, como se muestra en la Figura 2. Otro factor que ha contribuido a la desestabilización del sistema de dunas ha sido la remocion de la vegetación natural que sirve de protección a estos elementos, lo que ha contribuido a la perdida por erosion eólica. Los múltiples factores que modifican estos paisajes costeros aumentan el grado de vulnerabilidad, aun ante pequeñas variaciones del nivel del mar (Ortiz y Méndez, 1999).
Línea de costa
De acuerdo con los resultados obtenidos, uno de los sitios donde el retroceso costero es significativo corresponde al delta Grijalva-Usumacinta (Figura 3). Al igual que en los trabajos anteriores (Ortiz et al., 2010; Hernández et al., 2008), los resultados obtenidos coinciden en este retroceso de la costa. La situación es algo contradictoria si se piensa que en este sitio representa la desembocadura del sistema fluvial más caudaloso del país y por ende debería tener mayores tasas de sedimentación, lo cual se traduciría en una acreción de la costa mar adentro, sin embargo, hay que considerar los efectos de las corrientes costeras que representan un contrapeso importante para la acumulación de material proveniente de la cuenca alta.
El retroceso costero de este sitio es dinámico debido a la entrada de sedimentos, procesos geológicos internos y la corriente costera, por lo que alcanza pulsos de erosión periódica. Sin embargo, Ortiz y Méndez (1999) comentan que, debido a condiciones naturales, el delta Grijalva-Usumacinta presenta subsidencia de descenso asociado a fallas geológicas que ocasionaron desplazamientos activos, provocando mayor erosión de la playa y retroceso de la línea costera (West et al., 1969), aunque Nicholls y Cazenave (2010) y Gornitz et al., (1994) coinciden en que las actividades antropogénicas, como la extracción de fluidos (petróleo, gas, agua), incrementan la vulnerabilidad local al favorecer el hundimiento de la costa. Esto es especialmente importante, ya que la presencia de activos de producción de hidrocarburos es alta, lo que sumado a la subsidencia que mencionan Ortiz y Méndez (1999), incrementará el grado de vulnerabilidad. Por otro lado, la reducción de las fuentes de sedimentos está directamente vinculada al movimiento de línea de costa (Nicholls y Cazenave, 2010), al considerarse que anualmente las cuencas hidrológicas en condiciones naturales acarrean hacia las costas 20 millones de toneladas métricas (Milliman y Syvitski, 1992).
En este sentido, la cuenca Grijalva-Usumacinta que desemboca en este sector costero representa la mayor descarga al Golfo de México (del lado mexicano), sin embargo, el afluente del Grijalva se encuentra represado. Es claro que existe una disminución en la oferta de sedimentos debido a la modificación de esta cuenca, pero la situación puede ser aún más grave, ya que se encuentra en fase de proyecto la construcción de una presa en el afluente del Usumacinta, un afluente rodeado de zonas serranas que están sufriendo deforestación y representan fuentes eventuales de sedimentos; esto es especialmente importante para la costa tabasqueña, alimentada principalmente por los sedimentos acarreados por esta cuenca (West et al., 1969). El embalse de esta presa contribuirá a un mayor grado de vulnerabilidad si se considera que más de 100 millones de toneladas de sedimentos han sido atrapadas en diques y presas construidos en los últimos 50 años (Syvitski et al., 2005).
Pendiente costeraLa costa de Tabasco es prácticamente una planicie costera. Dada esa situación presenta pendientes inferiores al 10%. De acuerdo con los datos extraídos del LIDAR las pendientes máximas son de 9%, lo que indica que un proceso de inundación se extendería rápidamente tierra adentro de la costa tabasqueña (Carbajal, 2011). Este indicador representa directamente la fragilidad de la costa de Tabasco ante el proceso de inundación, por lo que se convierte en una variable dependiente de otras consideradas en este modelo, como la geomorfología.
OleajeLos rangos de oleaje más altos se presentan frente al sistema lagunar Carmen-Pajonal-Machona y las comunidades del Ejido Sinaloa 1ra, El Alacrán 1ra y Alacrán Sección Manatinero, en el municipio de Cárdenas; en el municipio de Paraíso las localidades de Barra de Tupilco, Guano Solo, Unión 1ra y 2da. Los sitios mencionados son los mismos que reporta Hernández et al. (2008) con altos índices de erosión; de acuerdo con ello, existe relación entre la altura de oleaje y el grado de erosión costera. Estos efectos son más fuertes durante la incidencia de fenómenos hidrometeorologicos de gran fuerza, como pasó en 1997 con los huracanes Roxanne y Opal (Galmiche-Tejeda y Solana-Villanueva, 2011).
Integrando el modelo de IVCLas mayoría de las variables que integran el IVC se relacionan entre sí mediante clasificaciones numéricas; estas, a su vez, indican un grado de vulnerabilidad que es asignado a cada sector de la costa según corresponda (Cuadro 5). Por su parte, la variable no numérica (geología/geomorfología) da un valor de vulnerabilidad clasificada cualitativamente de acuerdo con la resistencia relativa del tipo de material del paisaje costero. La tasa de erosión es significativa a lo largo del litoral, las tasas entre -2.32 y 3.31 m/año persisten a lo largo del litoral. Por su parte, la pendiente costera entre 0-2 m indica una clara vulnerabilidad de la costa ante los procesos de inundación, cubriendo un 85% del litoral. La tasa de INM rebasa los pronósticos regionales proyectados por el IPCC (2007), por lo que está clasificado como muy vulnerable.
Rangos de clasiicación de vulnerabilidad para cada variable en el área de estudio.
| Variables | Muy bajo 1 | Bajo 2 | Moderado 3 | Alto 4 | Muy alto 5 |
|---|---|---|---|---|---|
| Geomorfología | Planicie estructural baja, terraza estructural baja y dolina incipiente | Llan. aluv. alta, transic. de llan. alta-baja y plan. resid. inund., terr. monocl. baja, valle acumul., plan. conf. | Planicie marginal a la costa, llanura aluvial baja, cauce inactivo y cubeta de decantación. | Cordones de playa bajos, duna costera, depresión de turba y depresión de resurgencia. | Cordones de playa bajos con turba, llanura fluviomarina, depresión de turba salina y delta de marea. |
| Cambio de la línea de costa (m/año) | 0.12 - 12.39 | -1.44 - 0.12 | -2.32 - 1.44 | -3.31, - 2.32 | -17.98, -3.31 |
| Pendiente costera | > 9 | 4-8 | 2 - 4 | 1 - 2 | 0 - 1 |
| Cambio relativo del nivel del mar (mm/año) | -- | -- | -- | -- | > 6.65 |
| Altura de ola significante (m) | 0 - 0.38 | 0.38 - 0.74 | 0.74 - 0.83 | 0.83 - 0.89 | 0.89 - 0.94 |
| Rango mareal medio (m) | < 0.33 | -- | -- | -- | -- |
Una vez asignado el grado de vulnerabilidad a cada sector costero, todas las variables se combinaron mediante el álgebra de mapas para obtener el IVC (Figura 4). Es importante señalar que aunque el resultado final se relaciona con una cifra numérica, no debe asociarse con cambios específicos. El número resultante es una combinación de indicadores que al final permitirá identificar zonas específicas, que resultarán mayormente impactados por los efectos combinados del INM. Es claro que estos impactos se darán sobre la porción frontal del sector costero, por lo que a nivel de políticas públicas se ha avanzado en la gestión de esta zona considerándola como zona prioritaria de conservación POET (2013).
Los valores calculados se dividieron en cuatro clases, utilizando como límites los percentiles 25%, 50% y 75% para definir cuatro grados de vulnerabilidad (Ojeda et al., 2011). La utilización de los percentiles sirve para asignar un grado de vulnerabilidad relativa, sin que el valor como tal del IVC se asocie a cambios específicos de tipo físico (Sheik Mujabar y Chandrasekar, 2013). Los resultados finales del IVC indican valores entre 3.16 y 20.39. Los intervalos de vulnerabilidad total calculada por este índice se muestran en el Cuadro 6.
Las evaluaciones de vulnerabilidad representa una herramienta prioritaria para lograr el manejo sustentable de los recursos costeros al proporcionar datos esenciales para hacer una gestión adecuada de la zona costera (Sheik Mujabar y Chandrasekar, 2013; CEPAL, 2012). Una misión importante de estas evaluaciones es la de ayudar a los tomadores de decisiones en el rediseño del sistema costero (Farhan y Lim, 2011).
Al correr el modelo de IVC donde se analizaron las seis variables, la zona con el más alto valor de vulnerabilidad resultó ser el sector costero ubicado frente al complejo lagunar Mecoacán, que abarca una longitud lineal costera de 10km (Figura 5). En los últimos años esta laguna costera ha sido objeto de múltiples estudios, motivado principalmente por su cercanía con el Puerto Dos Bocas, que representa una de las zonas con mayor actividad industrial y que, posiblemente, está impactando de una u otra manera en el sistema lagunar. Medellín et al. (2013) por ejemplo, comentan que con mareas vivas la descarga del río González llega a ser nula; sin embargo, la verdadera explicación puede venir de las múltiples acciones antropogénicas que se desarrollan en la zona de influencia del puerto, además de los procesos de asolvamiento del canal del río.
Por otro lado, en la zona correspondiente al municipio de Cárdenas, a la altura de la localidad de Sánchez Magallanes, se presenta un efecto importante de vulnerabilidad, entre alta y muy alta. El sector vulnerable es de aproximadamente cinco kilómetros (Figura 6), esto es especialmente relevante dado que ahí se ubica una población importante del municipio de H. Cárdenas, trabajos recientes indican la vulnerabilidad de esta zona (Ortiz et al., 2010; Hernández et al., 2008).
Conclusiones
El modelo de medición adoptado de índice de vulnerabilidad costera fue adecuado para usarlo en la costa de Tabasco, sin embargo, la disponibilidad y accesibilidad de la información es muy limitada, por lo menos de parte de organismos, dependencias y algunas instituciones del país.
La variable que incide en menor medida dentro del modelo de índice de vulnerabilidad costera fue la marea, y las más importantes el oleaje y el nivel del mar.
De acuerdo con el índice de vulnerabilidad costera, las áreas cercanas a los sistemas lagunares Carmen-Pajonal-Machona y Mecoacán son de los sitios con mayor vulnerabilidad en la costa de Tabasco.
La vulnerabilidad muy alta y alta, que en conjunto forman seis kilómetros aproximadamente, están ubicadas por la comunidad de Sánchez Magallanes, cercana al sistema lagunar del Carmen, del municipio de H. Cárdenas, Tabasco, esto es debido a cuestiones de altos regímenes de oleaje asociado a una geomorfología de dunas costeras con material de sedimentos finos no consolidados.
Del litoral costero tabasqueño, la zona de muy alta vulnerabilidad es de unos 10km, ubicada frente al sistema lagunar Mecoacán. Es importante destacar que dicha área se encuentra inmersa en la zona de influencia del puerto de Dos Bocas, donde la actividad antropogénica es notable.
También es preciso destacar que este trabajo representa una de las primeras aproximaciones al cálculo de la vulnerabilidad local de la costa tabasqueña ante el potencial incremento del nivel del mar. Los resultados obtenidos a través de la implementación y ejecución de este modelo de índice vulnerabilidad costera son de gran alcance, dada la escala de estudio y la inclusión de variables en conjunto que por primera vez se realiza en el estado de Tabasco.