Se seleccionaron 30 parcelas de muestreo para las mediciones de biomasa y estructura del bosque de manglar por cada zona, a excepción de Potosí que era un área más pequeña con solo ocho parcelas. Cada parcela media 20 m x 10 m (200 m²). Las parcelas por cada sitio fueron seleccionados al azar, en función de la heterogeneidad de la cobertura del bosque de manglar, para ello, se hizo uso de imágenes satelitales de Sentinel 2B con tamaño de píxel de 10 m x 10 m, las que fueron preprocesadas con el software QGIS 3.2 (acceso libre). Los análisis estadísticos y gráficos fueron realizados con el software RStudio 3.5.0 (acceso libre).

Para conocer la estructura del bosque de manglar, en cada parcela se identificaron las especies, la densidad poblacional, la regeneración natural, daños antrópicos y daños naturales. Se realizaron medición del diámetro a la altura del pecho (DAP) y altura total para cada especie.

Para determinar la biomasa aérea (AGB) y biomasa subterránea (BGB), se utilizaron ecuaciones alométricas (Komiyama et al., 2005) ; tomando en cuenta los datos morfometricos para cada especie y el factor de densidad de cada especie (Bindu et al., 2018) .

Biomasa aérea

AGB = 0.251 * p * DBH ^2.46

Donde: AGB: Biomasa aérea, 0.251: Constante, : Densidad específica para cada tipo de especies, DBH: Diametro a la Altura del Pecho.

Biomasa subterránea (raíces)

BGB = 0.199 * p ^0.899 * DBH ^2.22

Donde: BGB: Biomasa subterránea, 0.199: Constante, : Densidad específica para cada tipo de especies, DBH: Diametro a la altura del pecho.

La suma de la AGB y BGB corresponde a la biomasa total de la parcela, el cual fue usado para calcular la biomasa en toda la cobertura boscosa del ecosistema de manglar para cada sitio de muestreo.

Se descargaron imágenes satelitales con ayuda del sensor Sentinel 2B (imágenes corregidas atmosféricamente) de la página web https://scihub.copernicus.eu/ de preferencia con un tamaño de píxel de 15 m x 15 m. Las imágenes fueron procesadas con el software QGIS 3.2 (acceso libre).

Una vez descargadas las imágenes se procedió a corregir vacios (gaps), los cuales son causados por una falla en el instrumento que corrige las líneas escaneadas, provocando vacíos en las distintas imágenes adquiridas del satélite, que se rellenan usando el Phase 2 Gap- Fill Algorithm, propuesto y utilizado por el servicio geológico de los Estados Unidos (USGS), basado en un ajuste lineal local del histograma (Rocha, 1999).

Con las imágenes corregidas se calculó el índice de vegetación de diferencia normalizada (NDVI) para estimar la vigorosidad de la vegetación (Rouse et al., 1974). El NDVI está basado en la absorción del espectro rojo y el infrarrojo cercano.

NDVI = (NIR - R) / (NIR + R)

Donde: NDVI: Índice de vegetación de diferencia normalizada, NIR: Espectro de infrarrojo cercano, R: Espectro rojo.

A partir de los datos de la biomasa aérea (AGB), biomasa subterránea (BGB) y los valores del NDVI para las parcelas seleccionadas, se realizó una regresión exponencial, a partir de la cual, se calculó la biomasa total para cada zona de muestreo, tomando en cuenta la formula de Myeong et al. (2006).

AGB = a * e ^(NDVI*b)

Donde a y b son constantes que se obtuvieron de la regresión entre el valor del píxel del NDVI y el AGB (muestreo de campo) de las parcelas seleccionadas. El valor de la parcela se ajustará al tamaño del píxel.

Las muestras de suelo se recolectaron en 10 parcelas seleccionadas al azar de las 30 parcelas del muestreo de estructura del bosque de manglar. La profundidad de muestreo fue de 20 cm, se depositaron en un recipiente limpio y seco, mezclándo hasta obtener una muestra homogénea de aproximadamente 1 kg y transferida a bolsas plásticas Ziploc®.

Posteriormente las muestras fueron secadas sobre un papel limpio y seco, por un periodo de 48 horas. Una vez seca se tamizaron y se molieron con un mortero y luego se pasaron por un tamiz de 2 mm de diámetro. Posteriormente se empacaron en bolsas plásticas para su posterior análisis químico del carbono almacenado en el suelo, a través de la técnica volumétrica en el laboratorio de química del suelo de la Universidad Nacional Autónoma de Nicaragua, León (UNAN-León).

Para el cálculo del carbono azul almacenado en el bosque de mangle, se usó el método estandarizado internacionalmente para las mediciones del almacenamiento de carbono azul en los ecosistemas, incluidos los manglares (Howard et al., 2014). Es así que, para crear el mapa de carbono se multiplicó el total de la biomasa (aérea y subterránea) por el factor de conversión 0.475 (47.5 % de biomasa) (Intergovernmental Panel on Climate Change (IPCC), 2013). Al resultado se sumó el valor del carbono total en el suelo, para obtener el valor total de reservorio de carbono en los bosques de manglar para cada zona de estudio.

Avicennia germinans y Rhizophora spp presentaron los valores más alto de biomasa por hectáreas, siendo muy heterogénea a lo largo de los seis sitios de muestreo (Figura 2). La isla Juan Venado registró os mayores valores de biomasa con la especie de A. germinans, contrario al resto de muestreo, característico de lugares altamente afectados por las acciones antrópicas (González Quiroz, 2016).

Figura 2. Figura 2. Biomasa aérea por especies según zona de muestreo

Potosí tiene los valores promedio más alto de biomasa por hectárea (Mg ha^-1), superando las 200 Mg ha^-1. El valor más bajo de biomasa registrado fue en Puerto Sandino (Cuadro 2).

Se identificaron 20 810.6 ha de cubierta de manglar, con distintos niveles de salud. El área identificada coincide con datos reportados por (Carvajal y Soto Valdivia, 2006), quienes indican que la zona del Pacífico Norte de Nicaragua tiene unas 19 700 ha cubiertas por manglar. Corinto con sus 5 976 ha, representó el 28.7 % del manglar del Pacífico Norte de Nicaragua, seguidos de Padre Ramos con 22.7 % y Aserradores con 20.2 % (Cuadro 2).

El valor más alto de biomasa por hectárea se registró en la zona de Potosí, coincidiendo con las mayores alturas y diámetros registrados (esta zona pertenece al área protegida Estero Real). Corinto y Padre Ramos registraron los mayores valores de biomasa en total, debido a su mayor cobertura de mangle. De las tres áreas protegidas (Isla Juan Venado, Padre Ramos y Estero Real), la Reserva Natural Isla Juan Venado presentó los valores más bajo de biomasa con 155.9 Mg ha^-1. El promedio de biomasa fue de 182.2±123.3 Mg ha^-1, valores similares a los reportados por (Donato et al., 2011; Intergovernmental Panel on Climate Change (IPCC), 2013; Simard et al., 2006). Estos valores suelen observarse en ecosistemas afectados por factores antrópicos (Velázquez-Pérez et al., 2019). Los resultados son inferiores a las medias globales (166.48±84.97 Mg ha^-1; solo biomasa area) de los bosques tropicales (Hu et al., 2016).

Cuadro 2. Cuadro 2. Biomasa acumulada en seis ecosistemas de manglar en el Pacífico Norte de Nicaragua

Zona de muestreo	Número de píxel	Área (ha)	Biomasa aérea (Mg Píxel-1)	Biomasa subterránea (Mg Píxel-1)	Biomasa (Mg ha-1)	Total de biomasa
Puerto Sandino	267 041	2 670.4	0.94±0.98	0.56±0.53	149.5	399 275.2
Isla Juan Venado	233 065	2 330.7	0.97±0.80	0.58±0.43	155.9	363 299.7
Corinto	597 661	5 976.6	1.14±0.77	0.68±0.41	181.6	1 085 543.0
Aserradores	420 751	4 207.5	1.13±0.65	0.68±0.35	180.6	759 971.5
Padre Ramos	472 272	4 722.7	1.29±0.96	0.75±0.50	204.6	966 211.8
Potosí	90 265	902.7	1.40±1.04	0.81±0.46	220.9	199 399.4
	2 081 055	20 810.6	x̅=1.15±0.87	x̅=0.68±0.46	x̅=182.2	3 773 700.6

Cuadro 3. Cuadro 3. Carbono acumulado en seis ecosistemas de manglar en el Pacífico Norte de Nicaragua

Zona de muestreo	Carbono Azul (Mg C)	Carbono en suelo (Mg C)	Carbono Azul Total (Mg C)	Carbono Azul (Mg C ha-1)	Carbono en Suelo (Mg C ha-1)
Puerto Sandino	189 655.7	155 708.2	345 363.9	71.0	58.3
Isla Juan Venado	172 567.4	146 752.1	319 319.5	74.0	63.0
Corinto	515 632.9	425 329.6	940 962.5	86.3	71.2
Aserradores	360 986.5	301 963.4	662 949.9	85.8	71.8
Padre Ramos	458 950.6	359 873.0	818 823.6	97.2	76.2
Potosí	94 714.7	72 479.0	167 193.7	104.9	80.3
	1 792 507.8	1 462 105.3	3 254 613.1	x̅=86.5	x̅=70.1

Figura 3. Figura 3. Distribución de biomasa en seis sistemas de manglar del Pacífico Norte de Nicaragua.

El sumidero de carbono azul para el Pacífico Norte de Nicaragua es de 1 792 507.8 Mg C para la parte aérea y subterránea de los árboles en las 20 810.6 ha de manglar estimadas. Esto representa 86.5±10.7 Mg C ha^-1 en promedio de almacenamiento para los seis sistemas de manglar evaluados (Cuadro 3). Los resultados son similares a los reportados por (Bautista-Olivas et al., 2018) en el bosque El Sargento en Sonora (México) en una zona bien conservadas (54.1-158.4 Mg C ha^-1), así mismo ligeramente superiores (36-65.5 Mg C ha^-1) a zonas altamente afectada por factores antrópicos localizados en la Bahía del Tóbari, Sonora. El nivel de impactos antrópicos sobre los ecosistemas de manglar tiene efectos sobre el almacenamiento de carbono (Blanco-Libreros et al., 2015).

En solo 20 cm de profundidad de suelo, se contabilizaron 1 462 105.3 Mg C almacenado en el suelo de las seis zonas de estudio (Cuadro 3), evidenciando la importancia de los manglares como sumideros de carbono (Howard et al., 2014). Los 20 cm representan tan solo el 23.1 % del total de carbono almacena en 100 cm de profundidad (Bhomia et al., 2016). Los valores de carbono en suelo están relacionados con la biomasa, siendo Padre Ramos y Potosí los que presentaron los mayores registros en promedio con 76.2 y 80.3 Mg C ha^-1. Sumado los dos componentes de almacenamiento de carbono [Carbono Azul (Mg C) y Carbono en Suelo (Mg C)], el bosque de manglar del Pacífico Norte de Nicaragua almacena 3 254 613.1 Mg C en total (aéreo, subterráneo y suelo), con los mayores registros para Corinto con 940 962.5 Mg C. En promedio se observaron 156.7±16.5 Mg C ha^-1 de carbono azul y carbono en el suelo para el Pacífico Norte de Nicaragua, superiores a los 122±16 Mg C ha^-1 reportados por (Doughty et al., 2016) en Florida, USA.

De acuerdo con (Zeng et al., 2021) Nicaragua tiene un potencial de financiamiento de 3 900±600 ha, aportando un potencial de mitigación de cambio climático de 51 000±14 000 t CO₂ por año^-1, y con un valor presente neto de 1 760 000±645 000 USD por año^-1. El total de carbono almacenado anualmente en los bosques de manglar del Pacífico Norte de Nicaragua podría ser parte del mercado de almacenamiento de carbono y además contribuir a la conservación a través de los fondos de este mecanismo. Los resultados evidencian que la cantidad de carbono almacenado podría aumentar en correspondencia con las estrategias de conservación. Es necesario garantizar mayor salud de bosque, restaurar y mantener las coberturas de manglar para conservar los importante e irremplazables servicios ecosistémicos del manglar, en especial el almacenamiento de carbono.

Figura 4. Figura 4. Carbono almacenado en los seis sistemas de manglar del Pacífico Norte de Nicaragua