Introducción
Las comunidades del bosque de manglar son de vital importancia ya que desempeñan valiosas funciones tales como: dar resguardo a las costas contra el desgaste hídrico y eólico (1) (2), dar un valor natural al aporte de la biodiversidad marino-costera (3), retener nutrientes y sedimentos (4), captura de carbono (4), exportar biomasa de los suelos costeros (5), proporcionar refugio y alimento para varios organismos marinos-dulceacuícolas en varias etapas de su vida y preservar el material genético de muchos individuos (6), (7). Por consiguiente, los manglares dan la estabilidad ecológica y productiva de gran cantidad de recursos marinos en las zonas tropicales del mundo (8), (9). A pesar de ello, en la actualidad hay muy poca generación de información científica acerca del cambio de uso de suelo (10), (11), (12) caracterización de entornos o definición de los tipos de cobertura relacionados con estos ecosistemas (13).
En Costa Rica se han realizado diversos estudios de la composición florística y estructura de los manglares del Pacífico (14), (15), (16), (17), (18), así como del Caribe (19), (20), (21). También se han descrito a nivel de flora y fauna los manglares de Térraba-Sierpe (22) y del Pacífico norte (23); sin embargo, a la fecha no se han realizado estudios científicos que cuantifiquen los desechos sólidos presentes en estos ecosistemas.
Un ejemplo de lo anteriormente dicho es la ausencia de información científica del Humedal Estero Puntarenas y Manglares Asociados (HEPyMA), el cual se ubica en el Pacífico central, y se estableció como área silvestre protegida mediante el decreto ejecutivo N° 29277 MINAE del año 2001. En el año 2006 se amplió su extensión, además se declaró de interés público toda iniciativa en favor de la preservación y recuperación de sus recursos naturales (Decreto ejecutivo N° 33327 MINAE 2006), lo anterior debido a que el estero presentaba altos índices de contaminación (24) lo que motivó a que diversas instituciones y miembros de la comunidad se unieran para atender y controlar la problemática (Decreto ejecutivo N°24693 MAG 1995) (25).
Si bien dicha iniciativa contribuyó con la mejora de la calidad del agua del HEPyMA, la contaminación de sus aguas aún persiste, así como otros conflictos que han afectado al humedal desde hace muchos años, tales como invasión de la frontera agrícola, instalación de salineras y estanques para el cultivo de camarón, tala, sedimentación y modificación de la red hidrográfica (26); como consecuencia entre el periodo de 1945-2018 se perdieron 865.7 ha de manglares, con una mayor tasa de pérdida entre 1945-1979 con 22.2 ha/año; para un total de 753 ha de manglar (27).
Actualmente las amenazas más fuertes que presenta el manglar son la sedimentación en el sector interno, la agricultura y la construcción de casas (27), además de la acumulación de basura, es por ello que el objetivo de este estudio fue evaluar la estructura y el impacto humano generado en la zona interna del HEPyMA.
Materiales y métodos
Área de estudio El Humedal Estero Puntarenas y Manglares Asociados se ubica en la provincia de Puntarenas, en la costa oriental del golfo de Nicoya. Se extiende desde Chacarita, en su extremo sur, hasta la comunidad de Chomes, al norte y hacia la ciudad de Puntarenas (27). Esta zona se caracteriza por tener un periodo seco que se extiende de diciembre hasta abril, y un periodo húmedo que se extiende de mayo a noviembre, durante esos periodos la temperatura fluctúa entre los 22 y 33°C (28), (29). Por otra parte, la precipitación media anual varía entre 700 mm y 2000 mm (30). Tiene una extensión aproximada de 57,7 km2, de los cuales 15,49 km2 pertenecen al área marina protegida, se encuentra bajo la administración del Área de Conservación Pacífico Central (3), y es bordeado por el estero de Puntarenas al oeste y por terrenos de uso agrícola por el lado este.
El estudio se realizó en marzo del año 2017 en un área comprendida entre las coordenadas 9° 59´26.90 N y 84°34´14.36 W (3), (31), la cual se seleccionó previamente mediante el análisis de mapas e imágenes satelitales con el fin de determinar la extensión y facilidades de acceso, además de la representación de los conflictos presentes en la zona interna de todo el manglar debido a su cercanía a poblados rurales y a la actividad agrícola intensiva (Figura 1).
Figure 1. Location of the study area in the inner zone of the Puntarenas mangrove.
Caracterización de la vegetación
Para los muestreos realizados en abril del año 2017 se establecieron 10 parcelas (por medio del análisis de las imágenes satelitales se determinó las áreas similares en composición) de monitoreo de 10 m x 25 m que se ubicaron estratégicamente en la zona interna del manglar, en los bordes más cercanos a terrenos de uso agrícola (Criterio de experto: (32)).
En cada parcela se cuantificó el número de árboles totales y por especie, las alturas y los diámetros (DAP) a 1,30 metros de altura, basándose en los criterios de (1) y (33). Además, se establecieron mediciones de árboles ≥2,5 cm, tanto vivos como muertos, para determinar alguna relación con la actividad antropogénica
Según lo establecido por (34), para el área de estudio y cada una de las parcelas se determinaron los siguientes atributos estructurales (E1):
Para la estimación del índice de complejidad de toda el área evaluada se promediaron los valores de altura, área basal y densidad, y se usó el número más alto de especies halladas en el sitio.
Caracterización fisicoquímica de agua: Para las tomas de agua de cada parcela se enterró un tubo de PVC de 1,20 metros (en cada parcela) y en su parte inferior tenía huecos que permitió que se acumulara el agua. Con ayuda de un barreno, se introducía por el tubo de PVC, se muestrearon solo 250 ml de agua por parcela de manera igual en todas.
En cada una de las parcelas, en la época lluviosa y durante la marea baja se colectaron n=10 muestras de agua, con el fin de caracterizar los parámetros fisicoquímicos con unmultiparámetros YSI 556: los parámetros físico-químicos: pH, Temperatura (Celsius), Oxígeno (%) = porcentaje de oxígeno, UPS= Salinidad y O2 mg (mg)= concentración de oxígeno.
Caracterización de la basura: Se usaron los criterios establecidos por Cordeiro y Costa (36) para la cuantificación de la basura en n= 5 parcelas transversales a los transectos realizados para la medición estructural (Fig. 1), con un largo de 200 metros y un ancho de 50 metros. El muestreo se realizó en abril del año 2019 durante la marea baja, para ello se llevaron a cabo recorridos simultáneos, sistemáticamente predefinidos a lo largo de cada parcela, durante los cuales se elaboraron listados para la clasificación de la basura, según número y tipo de basura. Los tipos de basura fueron también predefinidos mediante observaciones previas que determinaron ocurrencia y material de basura, de esta manera se establecieron como desechos sólidos más frecuentes: material plástico, vidrio, hules/poliuretano, tela, madera y zapatos. Se consideró además la cercanía e influencia portuaria, turística y urbana que tiene este manglar desde la ciudad de Puntarenas.
Procedimiento de análisis de la información: Para caracterizar el manglar, se usó la información estructural que se procesó mediante una tabla dinámica de Microsoft Excel. Además, se determinó la transformación de los datos y se realizó un escalamiento multidimensional (MDS) y el análisis de un clúster (37), (38) y (39).
Además, para evaluar la relación entre parcelas y los datos fisicoquímicos, en los casos de tener una distribución normal, se realizaron pruebas paramétricas (Software R Core TEAM), correlación de Pearson, ANOVA y un análisis de ordenación canónica (38), (39).
Para observar la diferencia entre parcelas y entre grupos de basura establecidos, se aplicó el índice de Simpson (D), además se calculó la abundancia relativa encontrada en cada parcela en el área de 200*50 metros. Por último, se realizó una correlación de Spearman, de la densidad de basura por grupo.
Resultados
Caracterización de la estructura: Se midieron un total de 217 árboles, distribuidos en 5 especies: Avicennia bicolor Standl, Avicennia germinans (L.) L., Laguncularia racemosa (L.) C.F. Gaertn., Rhizophora racemosa G. Mey. y Pelliciera rhizophorae Planch. & Triana (Cuadro 1). De ellos 151 árboles presentaron un diámetro mayor a 10 cm, 64 individuos tuvieron un diámetro mayor a 2.5 cm y menor a 10 cm, y solamente 2 árboles mostraron un diámetro menor a 2.5 cm. Además, se encontraron 13 árboles muertos, 8 de ellos en la parcela 2; y 6 ± 1.14 individuos secos. El índice de complejidad para todo el manglar fue de 4,6100, en cuanto a cada parcela el mayor índice de complejidad lo presentó la parcela 4 con un valor de 4,6398, seguido por la parcela cinco con un valor de 2,8025; por otra parte, el menor índice de complejidad lo obtuvo la parcela dos con un valor de 0,4367. Si bien, en la parcela cuatro se pudieron encontrar cinco de las siete especies de mangle nuclear presentes en el HEPyMA, además una de las alturas promedio y área basal total más altas, el alto índice de complejidad de esta parcela está dado mayormente por las especies del género Avicennia, lo mismo ocurre en la parcela cinco dominada completamente por las especies de dicho género.
Las especies que presentaron los mayores índices de valor de importancia fueron A. germinans presente en nueve de las diez parcelas y A. bicolor presente en siete de las diez parcelas. El menor índice de valor de importancia lo obtuvo L. racemosa con un 11%, presente solamente en la parcela 4 (Figura 2).
Figure 2. Importance value index of the species present in each plot established in the Puntarenas mangrove.
Procedimiento de análisis de la información: Para el análisis estadístico para la caracterización vegetal se denominó a las parcelas externas como m01, m03, m05, m07 y m09, y las más internas como m02, m04, m06, m08 y m10. Debido a que no hubo necesidad de transformar los datos se efectuó un escalamiento multidimensional (MDS) (Figura 3), el cual determina que los datos son muy homogéneos, y que las diferencias entre grupos se debieron principalmente a los valores de riqueza y abundancia de individuos por la mayor presencia de especies en las parcelas 4 y 10.
Figure 3. Importance value index of the species present in each plot established in the Puntarenas mangrove.
Por otro lado, la figura 4 representa la riqueza de especies entre las 10 parcelas por medio del análisis de un clúster, y el dendrograma establece 4 grupos según la composición de especies que se encontró en la parte interna del manglar de Puntarenas, donde nuevamente las parcelas 4 y 10 muestran la mayor diferencia debido a que presentaron más especies (5 y 4 respectivamente).
Figure 4. Multidimensional scaling (MDS) analysis of species richness per sampled plot in the inner part of the Puntarenas mangrove swamp.
Además, se realizó un PERMANOVA, para comparar si hay diferencias entre la composición de especies en las parcelas que están más internas o externas en el área de estudio. No se determinaron diferencias en la composición de especies es igual en todo el manglar en la zona interna, dado que no se encontraron diferencias significativas (F=2,99, 1 g.l., P>0,05), con un r2 =0,27, lo que confirma la homogeneidad del manglar presente en esta área.
En el Cuadro 2, los parámetros fisicoquímicos presentaron valores similares entre las parcelas, a excepción de la salinidad cuyo menor valor lo presentó la parcela 10. Además, la parcela 01 presentó el valor más alto de temperatura y el menor porcentaje de O2 en la muestra de agua.
Posteriormente, se hizo un análisis de parámetros físicos y químicos donde no se encontraron diferencias significativas (F= 1,4978, 1,8 gl., P>0,05). Primero, se realizó un análisis canónico de tendencia (DCA) con la longitud de gradientes y como sus valores eran menores a 4, se optó por hacer un análisis de redundancia (RDA), el cual explica la riqueza de las especies en un 65%, por variables físicoquímicas, no obstante, el modelo no es significativo (P>0.05).
Cuantificación y clasificación de basura acumulada en el manglar de Puntarenas: Para la cuantificación de la basura se denominó a las parcelas como M1, M2, M3, M4 y M5. Esto significa que en la zona interna del manglar de Puntarenas se hallaron 170 objetos clasificados en 6 grupos según su composición (Cuadro 3). Se debe destacar que los plásticos y zapatos fueron los dos grupos con mayor porcentaje de presencia en el manglar con valores de 22,4±29,58 y 7,4±15,44 respectivamente.
En la parcela 1 no se encontró ninguno de los 6 grupos de basura. Por otra parte, los grupos que se encontraron con mayor cantidad en las restantes parcelas fueron el plástico y los zapatos (Figura 5).
Figure 5. Total values (%) of the 6 garbage groups, according to their composition, in the 5 sampled plots in the inner part of the Puntarenas mangrove.
Discusión
Caracterización estructural:En el HEPyMA se han identificado las siete especies nucleares de manglar reportadas en el país, pero en el sector evaluado se encontraron cinco de ellas. Asimismo, esta área presentó un índice de complejidad total mayor al reportado para otros manglares del país como Palo Seco (18) y Rincón de Osa (17), pero más bajo que el reportado para Puerto Soley (19), playa Blanca y Escondido (17), y Potrero Grande y Santa Elena (16).
Según Jiménez y Soto (15) en las secciones internas del manglar de Puntarenas, donde la influencia de los ríos es baja, los suelos presentan salinidades más altas y la vegetación muestra un desarrollo estructural reducido. Lo anterior se evidencia al comparar el valor promedio de DAP de esta investigación (16,85 cm), con el obtenido por (36) que fue de 18,41 cm, y en cuyo estudio los transectos fueron establecidos en sitios más externos del mismo manglar.
En cuanto a la mayor presencia de las especies de Avicennia spp., podría explicarse a su preferencia por suelos más consolidados, pero de granulometría fina (suelos finos o arenas (40) que suelen encontrarse en el área interna del manglar, así como a la capacidad para desarrollarse con mayor facilidad que las otras especies nucleares en suelos arenosos (41). Asimismo, en las zonas internas del manglar donde la escorrentía y el drenaje están ausentes, la especie dominante suele ser A. germinans (25), (42), (43), (44), esta especie está presente en todo el Pacífico del país, pero es más abundante en el Pacífico seco (15) región a la que pertenece el HEPyMA. Por su parte, A. bicolor se ve favorecida por su capacidad de resistir las condiciones típicas de las zonas con climas secos pudiendo formar amplios rodales en sitios donde la salinidad es menor al 50% (45), (46).
Por otro lado, la presencia de L. racemosa suele estar muy restringida en los manglares, tanto de climas lluviosos como secos (46), (47). Sin embargo, estudios recientes han determinado que la especie es recolonizadora de las áreas donde se ha generado una fuerte extracción por tala de Rhizophora mangle (46), (42). Además, al igual que las especies de Avicennia spp., L. racemosa muestra mayor capacidad para desarrollarse en granulometrías más gruesas (41).
En el caso de P. rhizophorae se considera que es más sensible a suelos con alta salinidad y suele crecer en sitios elevados cercanos a lagunas salobres o drenajes de agua dulce (48), (49). Además, en condiciones controladas de laboratorio la especie presenta baja tolerancia a altas condiciones de luz y salinidad en las que sufre mayor mortalidad, mientras que en condiciones de agua fresca y luminosidad baja-media crece mejor, lo que demuestra la poca capacidad de la especie a tolerar condiciones estresantes, contrario a lo ocurrido con A. germinans que sobrevive bien en los diversos tratamientos, aunque en condiciones de alta salinidad disminuye su biomasa (43), lo anterior podría explicar la mayor abundancia de dicha especie en la zona de estudio. Por otra parte, al igual que R. racemosa; P. rhizophorae se distribuye de mejor forma en sedimentos lodo arenosos y lodosos donde suelen encontrarse formando rodales mixtos en los que R. racemosa, ubicada en zonas más externas; sirve como barrera otorgando protección a P. rhizophorae (41).
En relación con la salinidad, se ha determinado que existen patrones similares de distribución de las especies en función del gradiente de salinidad, de manera que las especies de L. racemosa suelen ubicarse en la zona media y superior, las de R. mangle lo hacen en las zonas inferiores y las de A. germinans suelen ubicarse en las zonas superiores de los estuarios (45), (50). Asimismo, las especies del género Avicennia spp. son consideradas como indicadoras de zonas con alta salinidad por las que muestran una mayor preferencia y ocurrencia (51), y suele dominar las zonas más altas del manglar lo que demuestra su gran adaptabilidad a condiciones desfavorables (31), (33); lo anterior podría explicar la mayor abundancia de las especies de este género en este sector del manglar.
Por otra parte, en este estudio se identificó que la parcela más cercana al río presentó el valor más alto de temperatura y el menor porcentaje de O2 en la muestra de agua, esto posiblemente se deba a la poca y lenta renovación del agua, lo que conlleva que haya poco oxígeno (31).
Caracterización de la basura: El plástico es el residuo sólido encontrado con más frecuencia en las playas a nivel mundial (52), (53). En el caso del manglar de Puntarenas también fue el desecho más recurrente, junto con los zapatos debido a que ambos materiales son ampliamente usados por la población dada su alta versatilidad (54).
Asimismo, estudios recientes reportan al plástico como el residuo sólido con mayor presencia en el manglar (55), (56), (57), (58), (59), (60), (61), dentro de las implicaciones negativas que genera su presencia está la disminución de espacio para el establecimiento de las madrigueras y disminución de diversidad de cangrejos (47), lo que a su vez genera deterioro de la calidad del suelo (57) además de cambios en su estructura, disminución de la calidad del agua, afectación en el crecimiento de las plántulas y establecimiento de propágulos, degradación de la belleza escénica, entre otros (58).
Por otro lado, (55), (59) y (60) encontraron que la cantidad de desechos sólidos varía en las distintas zonas del manglar de acuerdo con la cercanía a las comunidades; así como en las diferentes épocas del año de manera que en las temporadas altas de visitación a las playas adyacentes se incrementa la presencia de residuos en el manglar (54). Adicionalmente, (60), sugieren que los manglares podrían estar funcionando como sumideros de la basura proveniente tanto del mar como de la tierra, ya que los residuos suelen quedar atrapados en las ramas, troncos, y raíces de los árboles lo que impide que se trasladen y/o regresen al mar o playas contiguas.
Además, (55) reportan que los días con clima lluvioso y ventoso previos a los muestreos pudieron favorecer la acumulación de desechos sólidos livianos, principalmente plástico; en tanto que (59, 61) indican que la abundancia de basura y microplásticos también depende del estado de alteración del manglar, intensidad de las actividades pesqueras y los tipos de canales de agua presentes en el sitio evaluado, lo anterior podría explicar la escasa presencia de residuos en este sector del manglar de Puntarenas ya que el muestreo se llevó a cabo durante la época seca y en un área lejana a asentamientos humanos.
Conclusiones
El HEPyMA evidencia que estructuralmente es un ecosistema resiliente, que tiene características ecológicas que favorecen su alta diversidad y estabilidad a pesar de las diferentes actividades antropogénicas que se desarrollan en el manglar y sus alrededores.
Es posible que la flecha arenosa que conforma la ciudad de Puntarenas esté funcionando como barrera, impidiendo que una mayor cantidad de la basura arrastrada por el mar ingrese a este sector del manglar, y por ello se acumule principalmente en la playa.
A pesar de lo anterior, hubo presencia de desechos en el manglar que afectan negativamente el funcionamiento ecológico del ecosistema, por lo que es necesario y urgente tomar medidas de conservación con el fin de evitar daños irreversibles.