Área de estudio

El área de estudio comprende bosques de clima templado del municipio Pueblo Nuevo, específicamente, en los ejidos El Brillante, La Ciudad, La Victoria y La Campana del estado de Durango, México (Figura 1). Pertenece a la provincia fisiográfica de la Sierra Madre Occidental y a la subprovincia denominada Gran Meseta y Cañones Duranguenses. Sus elevaciones varían de los 1 400 a los 2 600 msnm. El clima es semihúmedo templado o semifrío, que se torna templado o semiseco en el lado oriental de la Sierra. La temperatura media anual oscila entre 20 y 22 °C, mientras que la precipitación media anual de 800 mm a 1 200 mm. Por su situación geográfica, reúne diversas condiciones de vegetación, desde masas puras de encino hasta bosques mezclados de pino-encino (^{Ucodefo 6, 1997}).

Figura 1 Ubicación geográfica. 

Trabajo de campo

La cuantificación de los almacenes de carbono se derivó de datos dasométricos de 39 sitios cuadrangulares de 50 x 50 m de la Red de Sitios Permanentes de Investigación Forestal, establecidos por la Universidad Juárez del Estado de Durango y el Instituto Tecnológico de El Salto (^{Corral et al., 2009}). La estimación de biomasa en cada sitio de muestreo se realizó por medio de ecuaciones alométricas previamente generadas para las especies más importantes (producción de biomasa) del estado de Durango, y se aplicaron a los datos dasométricos obtenidos en los sitios de muestreo. Los resultados de biomasa se convirtieron a cantidad de carbono almacenado, mediante la multiplicación de biomasa por el factor de carbono total específico para cada especie. Dichos porcentajes se distribuyen entre valores de 0.43 y 0.51. Los modelos fueron validados a partir de los datos de 10 sitios permanentes de muestreo establecidos en el ejido La Campana, debido a que esta área tiene condiciones ecológicas y silvícolas similares a los ejidos El Brillante, La Ciudad y La Victoria.

La determinación de biomasa y contenido de carbono total en cada sitio de muestreo se hizo con la sumatoria de la cantidad de ambos por árbol. Asimismo, los resultados de dichas variables procedentes de los 39 sitios fueron extrapolados a toneladas por hectárea, con la finalidad de facilitar la generación y construcción de los modelos para predecir su estimación, a través de la percepción remota.

Datos espectrales Landsat TM (Thematic Mapper)

La imagen de satélite se basó en la página del Servicio Geológico de los Estados Unidos (^{USGS, 2013}), y se tomó el 19 de junio de 201, con una resolución espacial de 30 m. Está ortorrectificada; es decir, representa los objetos del terreno en sus verdaderas coordenadas. Este sensor opera con siete bandas espectrales, que fueron elegidas especialmente para el monitoreo de la vegetación a excepción de la banda 7, que se agregó para aplicaciones geológicas (Cuadro 1).

Cuadro 1 Características de las bandas y principales aplicaciones que componen el satélite Landsat TM. 

Bi = Bandas espectrales; IRC = Infrarrojo cercano; IRM = Infrarrojo medio; IRT = Infrarrojo termal; μm = Micras.

Los datos espectrales fueron extraídos de las coordenadas UTM centrales de los sitios de muestreo, los cuales fueron identificados por los pixeles de las imágenes. Los datos espectrales consistieron en las reflectancias obtenidas de cuatro bandas individuales de la imagen y de cuatro transformaciones matemáticas aplicadas a la reflectancia: 1) Índice de vegetación de diferencia normalizada (NDVI), compuesto por las bandas espectrales infrarrojo cercano (IRC) y rojo (R), calculado como NDVI = (IRC-R)/(IRC+R) (^{Rouse et al., 1974}); 2) NDVI41, compuesto por las bandas espectrales infrarrojo medio (IRM) y verde (V), calculado como NDVI41 = (IRM-V)/(IRM+V); 3) NDVI42, calculado como NDVI42 = (IRM-R)/(IRM+R), se compone por las bandas espectrales infrarrojo medio y rojo; 4) NDVI43, calculado como NDVI43 = (IRM-IRC)/(IRM+IRC), el índice es conocido como de estrés hídrico (^{Rock et al., 1986}), y se integra por las bandas espectrales infrarrojo medio e infrarrojo cercano. Los índices se construyeron con el fin de resaltar algunas características de la vegetación, la cual está en función, principalmente, de la clorofila, de la estructura celular y del contenido de agua de la vegetación (^{Chuvieco, 1998}; ^{Soria y Granados, 2005}; ^{Aguirre et al., 2009}); y son más sensibles en alguna parte del espectro electromagnético, además se detectan mejor en bandas específicas (^{Cohen y Goward, 2004}).

Relación entre datos de campo y datos espectrales

Se aplicó un análisis de correlación entre las variables forestales de interés y los datos espectrales con el fin de averiguar el comportamiento y relación entre ellos. Se empleó el procedimiento estadístico de regresión por pasos (Stepwise), con el objetivo de seleccionar un modelo de regresión lineal múltiple que estime la biomasa (Mg ha^-1) y el carbono (Mg ha^-1) con un mínimo de variables. Como indicadores estadísticos de bondad de ajuste se consideraron el coeficiente de determinación (R²), la raíz del cuadrado medio del error (RMSE) y un nivel de significancia menor a 0.05 para cada variable. Una vez desarrollados los mejores modelos de regresión, se estimó la biomasa y el carbono para los sitios de muestreo que corresponden a los datos espectrales, a partir de lo cual se generaron mapas de biomasa y carbono arbóreo aéreo para las localidades de estudio.

Validación

Con apoyo en la metodología antes descrita, se procesaron datos dasométricos y coordenadas UTM centrales de 10 sitios de muestreo, distintos a los 39 empleados para la construcción de los modelos, los cuales se ubicaron dentro del ejido La Campana; los resultados obtenidos de la aplicación de los modelos de regresión, predichos, fueron comparados con los datos de campo, observados, estimados por medio de ecuaciones alométricas y factores de carbono.

Relación entre las variables forestales y los datos espectrales

Las variables forestales biomasa y carbono mostraron una correlación negativa con los datos espectrales (Cuadro 2). De acuerdo con ^{Hall et al. (2006)} y ^{Aguirre et al. (2009)}, lo anterior se explica por la disminución del albedo en zonas con vegetación densa y cerrada. Aunque se presentó poca correlación entre la biomasa y el carbono con las reflectancias de las bandas 1, 2, 4 y el índice de vegetación NDVI42; el índice de vegetación NDVI43 tuvo la correlación más alta, con un valor de -0.62 para biomasa y -0.63 para carbono. ^{Aguirre et al. (2009)} registraron un valor de correlación de -0.80 entre el NDVI43 y el carbono arbóreo aéreo. Los resultados obtenidos corroboraron lo señalado por ^{Gong et al. (2003)} y ^{Aguirre et al. (2009)}, respecto a que la densidad de los bosques de coníferas (hoja acicular) tienen una mejor explicación por índices que son sensibles a la humedad, que por aquéllos que son sensibles a la clorofila. El signo negativo indica que la densidad forestal es inversamente proporcional al estrés hídrico de la vegetación (NDVI43). Por lo tanto, el NDVI presenta buenas correlaciones con las variables forestales de interés, debido a la reflectividad ocasionada por la clorofila de la vegetación. En concordancia con ^{Chuvieco (1998)} y ^{Soria y Granados (2005)}, cuando la vegetación es sometida a un tipo de estrés (plagas, sequía, etcétera.), su reflectancia en NDVI será inferior y, por lo tanto, un área con vegetación sana presentará un mayor contraste entre las bandas que componen este índice (infrarrojo cercano y rojo).

Cuadro 2 Coeficientes de correlación entre las variables forestales y los datos espectrales. 

B = Biomasa (Mg ha^-1); C = Carbono (Mg ha^-1); Bi = Reflectancia por banda espectral; NDVI = Índice de vegetación de diferencia normalizada.

Modelos seleccionados

Los índices de vegetación basados en la reflectancia de las bandas verde, rojo, infrarrojo cercano y medio, fueron las que describieron de manera satisfactoria el comportamiento de la biomasa y el carbono, debido a la fuerte absorción de energía causada por la presencia del vigor de la vegetación sana, la absorción de la clorofila y la humedad presente en la vegetación (^{Soria y Granados, 2005}). Los modelos seleccionados, mediante el procedimiento de regresión (Stepwise) para estimar biomasa y carbono a partir de datos espectrales se presentan en el Cuadro 3.

Cuadro 3 Modelos de regresión seleccionados. 

B= Biomasa (Mg ha^-1); C= Carbono (Mg ha^-1); B2 = Reflectancia en banda verde; B3 = Reflectancia en banda rojo; NDVI = Índice de vegetación de diferencia normalizada.

En diversos estudios se han utilizado datos espectrales para la estimación de variables forestales, que han presentado niveles de precisión (RMSE) y coeficientes de determinación (R²) muy variados. ^{Aguirre et al. (2009)} calcularon 22.24 Mg ha^-1 de RMSE para carbono en un bosque de Pinus patula Schiede ex Schltdl. & Cham. en la región central de México con datos del sensor SPOT 5 HRG; realizaron estimaciones con regresión lineal múltiple y el procedimiento no paramétrico del vecino más cercano (k-nn), cuyo valor de R² fue 0.70 para el modelo de regresión seleccionado para cuantificar carbono con un nivel de confiabilidad de 99 %.

Por su parte, ^{Labrecque et al. (2006)} determinaron valores de RMSE de 59, 59 y 58 (Mg ha^-1) para bosques de coníferas, caducifolios y mixtos, respectivamente, para mapear biomasa forestal en una región piloto situada al oeste de Canadá, a través de imágenes Landsat TM. ^{Escandón et al. (1999)} usaron el mismo sensor, pero con un enfoque diferente, y evaluaron dos métodos para estimar biomasa arbórea. El primer método consistió en una clasificación supervisada multiespectral con seis bandas; utilizaron tipos de vegetación identificados a partir de la composición de la biomasa de los géneros dominantes y de la altura promedio estimada del dosel. En el segundo método trabajaron con índices de vegetación de diferencia normalizada (NDVI) de las bandas TM4/TM3; TM4/TM5 y TM4/TM7; posteriormente, aplicaron un modelo de regresión que relaciona la biomasa promedio con los valores o niveles digitales (ND) con los NDVI. Resultó que el modelo exponencial fue el de mejor ajuste para los tres NDVI; sus valores para los NDVI fueron TM4/TM3 (R²= 0.61), TM4/TM5 (R²= 0.67) y TM4/TM7 (R²= 0.68). Comparativamente, los citados por ^{Aguirre et al. (2009)} son más similares a los del presente estudio y se consideran como aceptables para este tipo de investigaciones.

Las diferencias identificadas respecto a las experiencias previamente citadas pueden deberse a los satélites, por efecto de la resolución espacial que ofrecen el LANDSAT y el SPOT (^{Aguirre et al., 2009}) y a la resolución espectral (^{Soria y Granados, 2005}).

Comparación de las estimaciones biomasa-carbono con los modelos seleccionados

En la Figura 2 se muestra la relación de los valores de biomasa y de carbono observados en campo y los predichos con los modelos de regresión desarrollados, como una forma viable de ser utilizados para un bosque. La variable empleada fue el índice de estrés hídrico (NDVI43) debido a que fue el que resultó altamente significativo (p≤0.0001) y presentó mayor correlación con las variables de interés. Se advierte una relación proporcional inversa, ya que al disminuir el valor del índice de estrés hídrico (NDVI43), aumentan los de biomasa y carbono, lo que refleja el potencial de captura de carbono que hay en este tipo de vegetación.

Figura 2 Valores de biomasa (Mg ha^-1 ) y carbono (Mg ha^-1 ) observados y predichos construido con el NDVI43. 

Aplicación y validación de los modelos

Los modelos seleccionados fueron aplicados a las reflectancias de las bandas e índices de vegetación correspondientes. Al extrapolarlos al área de estudio, se obtuvo la cantidad de biomasa y el contenido de carbono arbóreo aéreo.

En el Cuadro 4 se consigna la comparación de las estimaciones de biomasa y carbono arbóreo aéreo total en campo y las derivadas de los modelos de regresión generados.

Cuadro 4 Biomasa (Mg ha^-1) y contenido de carbono (Mg ha¹), estimados en campo vs. modelos generados. 

Las valoraciones anteriores se pueden mejorar, si se cuenta con datos de sensores remotos con una resolución espectral mejorada (^{Thenkabail et al., 2004}). Al respecto, los resultados obtenidos coinciden con las conclusiones de ^{Fernández et al. (2005)}, quienes afirman que, a pesar de ciertas limitaciones, las imágenes de esta procedencia son válidas para la estimación de parámetros forestales a una escala regional.

Generación de mapas de biomasa y carbono arbóreo aéreo

Con base en el análisis del valor máximo de biomasa y carbono de los sitios de muestreo, se determinaron cifras para los pixeles de las imágenes de 400 Mg ha^-1 para la biomasa y de 200 Mg ha^-1 para el contenido de carbono. Para facilitar la interpretación de ambos mapas, los datos se agruparon en clases e intervalos (Cuadro 5).

Cuadro 5 Clases e intervalos de valores utilizados para la generación de los mapas. 

Las figuras 3 y 4 exhiben los mapas generados de la biomasa para los ejidos El Brillante y La Campana y contenido de carbono total para los predios de los ejidos La Ciudad y La Victoria.

Los diferentes colores representan los valores de biomasa y contenido de carbono arbóreo aéreo en toneladas por hectárea para los ejidos evaluados en los mapas generados para las regiones forestales. Por ejemplo, al tomar como parámetro el mapa de carbono para el ejido La Victoria, cada color se refiere al contenido de carbono en las áreas boscosas con un valor máximo de 200 toneladas.

Asimismo, este tipo de mapas puede ser usado para la atención de las causas de la pérdida de bosques y del carbono forestal, al fortalecer el manejo de los bosques y la conservación de su biodiversidad.

Figura 3 Mapas de biomasa total (Mg ha^-1 ). 

Figura 4 Mapas de carbono arbóreo aéreo total (Mg ha^-1 ).