texto en 
Inglés (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
El Parque Nacional Cumbres de Monterrey (PNCM) forma parte del sistema montañoso Sierra Madre Oriental; posee un amplio gradiente altitudinal, así como una diversidad de exposición a la insolación e influencia de masas de aire húmedo que favorecen un extraordinario mosaico de ecosistemas y diversidad biológica; por ello, se le considera el área de mayor relevancia y valor ecológico del estado de Nuevo León, principalmente, para el Área Metropolitana de Monterrey, ya que suministra más de 40 % del agua que consume la población (Cantú-Ayala et al., 2013; Rovalo-Merino et al., 2013).
En el PNCM se distribuyen 3 729 especies de flora y fauna, de las cuales 141 están en algún estatus de riesgo: mamíferos (11), aves (32), reptiles (36), anfibios (4), peces (8), plantas (42) y hongos (8) (Semarnat, 2010; González-Saldivar et al., 2013; Naturalista, 2021). Sus ecosistemas capturan 6 113 920 Mg C, lo que representa una contribución muy importante a la mitigación de los efectos de las emisiones de gases de efecto invernadero (Jiménez-Pérez et al., 2013).
Al respecto, las comunidades vegetales con las cifras más altas son los bosques de pino-encino, con un promedio de 82 Mg C ha-1 (Rodríguez-Laguna et al., 2009); encino-pino, con 70 Mg C ha-1 (Aguirre-Calderón y Jiménez-Pérez, 2011); bosques de pino, 62 Mg C ha-1 (Pimienta de la Torre et al., 2007); matorral submontano, con 12 Mg C ha-1 (Montaño et al., 2016); y el pastizal, 4 Mg C ha-1 (Yerena et al., 2014).
El Parque Nacional Cumbres de Monterrey es una de las áreas más visitadas en Nuevo León, debido a sus atractivos naturales idóneos para el esparcimiento y la recreación (Correa-Sandoval y Mayén, 2013). Estas actividades pueden convertirse en un factor significativo para conservar el ambiente, siempre y cuando se consideren los factores ecológicos, sociales y económicos de la región (Díaz, 2010). Sin embargo, hay prácticas como el uso de vehículos todo terreno que se consideran como una de las causas de mayor disturbio ecológico e hidrológico a los cauces de ríos y arroyos (Brenner, 2006; Menchaca y Alvarado, 2011); la cual se ha incrementado en la zona.
La situación antes descrita ejerce una fuerte presión sobre los recursos naturales disponibles, con la consecuente degradación de los ecosistemas vulnerables (Rodríguez y Acevedo, 2015); además, existe crecimiento de las comunidades que habitan en el Parque y falta de áreas idóneas para los cultivos agrícolas (Aragón-Palacios, 2013).
El análisis multitemporal es de gran importancia para el monitoreo de la dinámica en el cambio del uso del suelo y de la vegetación, ya que permite comprender los factores que detonan el proceso de sucesión ecológica, que a su vez, inciden en la recuperación de las funciones vitales de los ecosistemas que han sido expuestos tanto a impactos naturales, como a efectos de los disturbios antrópicos (Giri et al., 2007; Menchaca y Alvarado, 2011). Asimismo, el uso de imágenes satelitales de alta resolución es una importante alternativa para evaluar variables estructurales y ecofisiológicas a escala regional y a través del tiempo, para analizar la evolución de los ecosistemas forestales posterior a las alteraciones en su estructura y composición (Tirpak y Giuliano, 2010).
En el PNCM se han desarrollado pocas investigaciones en las que se evalúe la dinámica de sus comunidades vegetales; la mayoría de los estudios se han enfocado en el crecimiento de las especies de pinos (González, 2019) y en la reconstrucción climática mediante dendrocronología (Luna, 2020); sin una perspectiva a nivel de paisaje que permita comprender el efecto que ha tenido la exclusión de las actividades productivas en la restauración pasiva de las comunidades vegetales o en la degradación de los ecosistemas posincendio.
El objetivo de la presente investigación fue realizar una clasificación del uso del suelo y vegetación en tres años: 2000, 2010 y 2018. Como hipótesis se planteó que la superficie de las comunidades vegetales se recupera de manera natural conforme pasa el tiempo, por la exclusión de actividades productivas en el PNCM.
Materiales y Métodos
El PNCM se ubica en el estado de Nuevo León, entre las coordenadas geográficas 26°31´00" de latitud norte, 100°17´20" longitud oeste (Figura 1) y colinda con Coahuila. Forma parte de la cuenca hidrográfica RH24 Río Bravo-Conchos, situada entre los límites de los municipios: Allende, García, Montemorelos, Monterrey, Rayones, Santa Catarina, Santiago y San Pedro Garza García, Nuevo León; con una superficie total de 177 395.95 ha (DOF, 2000).
Los climas dominantes en el PNCM son Semicálido húmedo (A)C(w1), (A)C(w2); Templado subhúmedo C(w1); Árido cálido BSohw y Semiárido cálido BS1hw. La precipitación pluvial varía de 344 mm a 983 mm al año y la temperatura promedio anual es de 9 a 24 °C (Cuervo-Robayo et al., 2014). Los suelos son de tipo semiárido, asociados con vegetación desértica; la mayoría son poco profundos, de texturas gruesas y en ocasiones presentan subsuelos duros o poco permeables; los dominantes son Leptosol (78.88 %) Pheozem (17.68 %), Fluvisol (1.45 %), Luvisol (0.98 %), Calcisol (0.55 %), Regosol (0.31 %) y Vertisol (0.15 %) (Inegi, 2019).
Las aguas superficiales de la zona de estudio son drenadas, principalmente, por las subcuencas del río Santa Catarina y el río Ramos; pequeñas porciones del PNCM corresponden a otras subcuencas: por la parte noroeste drena a la subcuenca del río Pesquería; la región sureste, al río Pilón; y por su región central, hacia el noreste, el río San Juan (SIATL, 2020).
Procesamiento y generación de ortomosaicos
Con la finalidad de identificar la resolución óptima en procesos de clasificación supervisada para determinar las tasas de cambio de uso de suelo y vegetación, se realizó un comparativo entre las imágenes satelitales de mayor uso en este tipo de análisis: Landsat 8 OLI (Figura 2A) y Sentinel 2A (Figura 2B), en comparación con ortofotos (Figura 2C), imágenes satelitales Birdseye (Figura 2D) y de Airbus Defence and Space.
Figura 2 Comparativo de imagen satelital Landsat 8 OLI (A), Sentinel 2A (B), Ortofotos (C) e imágenes de alta resolución (D).
De acuerdo con Earth Observing System (2020) las imágenes Landsat 8 (30 m/píxel; bandas 7, 5, 4) y Sentinel 2 (10 m/pixel; bandas 8,4,3) se clasifican como de resolución media; mientras que, las ortofotos (1.5 m/pixel), las imágenes satelitales Birdseye (0.28 m/píxel) y las Airbus Defence and Space (1.14 m/pixel) se consideran de alta resolución, porque que tienen una resolución por de debajo de los 5 m/pixel de la clasificación de Earth Observing System.
El procesamiento de las imágenes satelitales de alta resolución se basa en la georreferenciación, corrección radiométrica, normalización radiométrica de la serie y porcentaje de nubosidad (Tirpak y Giuliano, 2010); por ello, la disponibilidad de imágenes con estas características se limitó a los años 2000, 2010 y 2018.
Las imágenes satelitales de media resolución se obtuvieron de la plataforma Earth Explorer del Servicio Geológico de los Estados Unidos de América (USGS, 2020), las ortofotos del servidor Espacios y Datos de México (Inegi, 2020) y las imágenes Birdseye y Airbus Defence and Space del software SASPlanet (SASPlanet, 2020).
Clasificación supervisada
Se realizó una clasificación supervisada, a partir de la generación de tres ortomosaicos, compuestos por 27 ortofotos (año 2000; Figura 3A), 578 imágenes de Birdseye (año 2010; Figura 3B), 144 de Airbus Defence and Space (año 2018; Figura 3C); las cuales se compararon con imágenes Landsat (año 2018; Figura 3D) y Sentinel (2018).
Para la detección de cambios se utilizó el software QGIS 2.18.25 “Las Palmas” de código abierto (QGIS Development Team, 2020).
La corrección atmosférica a las imágenes de cada periodo se hizo recortándolas y sometiéndolas a un proceso de clasificación no supervisada con el módulo K-means analysis, que agrupa los valores de celda en clases con el método de análisis de conglomerados de datos multivariados (Jumb et al., 2014; Rashmi et al., 2016); posteriormente, se transformaron los archivos de formato raster a vectorial, para una clasificación supervisada (Figura 4).
La clasificación supervisada se llevó a cabo mediante puntos de control, que consistieron en 148 sitios distribuidos de manera sistemática en áreas con alto grado de confusión por la reflectancia de las imágenes, la exposición, ruido y nubosidad.
Se obtuvo información relacionada con los diferentes usos del suelo, cobertura vegetal y tipo de vegetación, la cual se cotejó con la desarrollada por Conafor (2013) e Inegi (2017), y se generó una clasificación de diferentes usos del suelo (agricultura, asentamientos humanos, caminos, líneas de transmisión eléctrica y ríos), así como cobertura forestal (Bosque de Ayarín, Encino, Encino-Pino, Pino, Pino-Encino, Matorral Desértico Micrófilo, Matorral Desértico Rosetófilo, Matorral Submontano, Pastizal y Sin Vegetación) distribuidas en el PNCM (Figura 5).
La concordancia y precisión de los resultados de clasificación de las imágenes satelitales de alta resolución se calculó en el módulo r.kappa en GRASS 7.6.0 (QGIS Development Team, 2020), el cual genera una matriz de error y determina el Coeficiente Kappa de Cohen.
Determinación de pérdidas y ganancias
La determinación de la pérdida o ganancia en el proceso de restauración de los diferentes tipos de vegetación, se obtuvo mediante la construcción de matrices de transición y tasas de cambio para los años considerados en el estudio, mediante la ecuación desarrollada por la FAO (1996) y adaptada por Palacio-Prieto et al. (2004).
Donde:
Para describir los cambios de uso del suelo y vegetación, se generó una tabulación cruzada del tiempo 1 y del tiempo 2, lo cual permite obtener una matriz de cambio. La diagonal muestra la cantidad total del paisaje estable entre una fecha y otra, mientras que fuera de la diagonal se ubican las transiciones de las clases entre el Año 1 y el Año 2. Un valor positivo de tasa de cambio indica ganancia de superficie restaurada, mientras que a un valor negativo le corresponde una pérdida de cobertura forestal.
Resultados y Discusión
A partir del proceso de clasificación supervisada se derivaron 15 clases de uso de suelo y vegetación; en las categorías de asentamientos humanos, caminos, líneas de transmisión eléctrica y ríos no se observó un cambio significativo; por ello, el análisis se realizó a partir de las 11 clasificaciones principales (Cuadro 1).
Cuadro 1 Clasificación del uso del suelo y vegetación en el PNCM.
| Año | 2000 | 2010 | 2018 | |||
|---|---|---|---|---|---|---|
| Vegetación | ha | % | ha | % | ha | % |
| A | 1 550.00 | 0.87 | 1 701.33 | 0.96 | 1 865.57 | 1.05 |
| BA | 3 600.69 | 2.03 | 3 696.70 | 2.08 | 3 760.52 | 2.12 |
| BQ | 20 119.76 | 11.34 | 20 531.17 | 11.57 | 21 035.22 | 11.86 |
| BQP | 15 148.80 | 8.54 | 15 518.38 | 8.75 | 15 915.20 | 8.97 |
| BP | 24 957.86 | 14.07 | 25 999.67 | 14.66 | 26 866.54 | 15.14 |
| BPQ | 26 276.07 | 14.81 | 26 468.60 | 14.92 | 27 336.31 | 15.41 |
| MDM | 219.14 | 0.12 | 217.51 | 0.12 | 197.79 | 0.11 |
| MDR | 13 664.03 | 7.70 | 12 463.80 | 7.03 | 11 949.78 | 6.74 |
| MS | 37 452.30 | 21.11 | 37 820.18 | 21.32 | 39 751.12 | 22.41 |
| P | 1 152.91 | 0.65 | 1 030.07 | 0.58 | 968.54 | 0.55 |
| SV | 25 860.23 | 14.58 | 24 554.38 | 13.84 | 20 355.19 | 11.47 |
A = Agricultura; BA = Bosque de Ayarín; BQ = Bosque de Encino; BQP = Bosque de Encino-Pino; BP = Bosque de Pino; BPQ = Bosque de Pino-Encino; MDM = Matorral Desértico Micrófilo; MDR = Matorral Desértico Rosetófilo; MS = Matorral Submontano; P = Pastizal; SV = Sin Vegetación.
El análisis comparativo entre imágenes Landsat 8 OLI y Sentinel 2ª, con respecto a las de alta resolución demostró que existe una sobreestimación superior a 24 % para superficies con cobertura forestal, debido a que algunas imágenes presentaron desplazamiento entre las coordenadas de las bandas, lo cual puede producir un error sistemático relacionado con la sobreposición en el mismo path, pero diferente row, en combinación con el grado de nubosidad en las imágenes (Cristóbal et al., 2004; Astola et al., 2019).
La clasificación de las imágenes de alta resolución presentó valores promedios de índice de Kappa de 0.83, que se consideran dentro de un nivel de precisión muy buena, y fueron superiores a lo citado por Mendes et al. (2015), quienes obtuvieron un valor de 0.58 para imágenes del satélite Geoeye-1.
El PNCM abarca una superficie de 177 395.87 ha, de las cuales para el año 2018 (Cuadro 1), 81.25 % correspondió a cubierta forestal, 14.58 % sin vegetación aparente, 3.18 % a ríos, 0.70 % asentamientos humanos, 0.19 % caminos y 0.10 % a líneas de transmisión eléctrica. Entre los principales tipos de vegetación identificados destacaron: matorral submontano (39 751.12 ha), pino-encino (27 336.31 ha), pino (26 866.54 ha), encino (21 035.22 ha), encino-pino (15 915.20 ha), matorral desértico rosetófilo (11 949.78 ha), bosque de ayarín (3 760.52 ha) y matorral desértico micrófilo (197.79 ha).
En general se aprecia un proceso de recuperación de la cobertura forestal en el PNCM; principalmente, en bosques de pino, pino-encino, encino-pino, encino y ayarín como resultado de la capacidad de resiliencia de los ecosistemas ante disturbios naturales o antropogénicos (Lloret, 2012; Mora-Donjuán y Alanís-Rodríguez, 2016), así como por los esfuerzos de restauración que han realizado las comunidades, organizaciones no gubernamentales e instituciones públicas ( Rovalo-Merino et al., 2013).
Sin embargo, en ecosistemas como el matorral desértico micrófilo, matorral desértico rosetófilo y pastizales continúa la pérdida de cobertura, debido a los procesos de cambio de uso del suelo con fines agrícolas y ganaderos, así como por el aumento de los asentamientos humanos (Figura 6).
Detección de cambio de uso del suelo
La recuperación de la vegetación en el PNCM es positiva, como se aprecia durante el periodo 2000-2010, en el que se registró una recuperación de 13.03 % en la cobertura forestal, con respecto al año 2000; principalmente, en áreas de bosques de pino, ayarín, encino-pino, encino, pino-encino y en el matorral submontano. Sin embargo, también hubo una pérdida de 20.18 % en los pastizales, matorral desértico rosetófilo y matorral desértico micrófilo; en cuanto a la agricultura se determinó un aumento de 9.76 % y las superficies sin cobertura aparente presentaron una disminución de 5.05 % (Cuadro 2).
Cuadro 2 Cambios en la cobertura del suelo de 2000 a 2018 en el PNCM.
| Período | ∆ 2000 - 2010 | ∆ 2010 - 2018 | ∆ 2000 - 2018 | |||
|---|---|---|---|---|---|---|
| Vegetación | ha | % | ha | % | ha | % |
| Agricultura | 151.33 | 9.76 | 164.24 | 9.65 | 315.57 | 20.36 |
| Bosque de Ayarín | 96.01 | 2.67 | 63.82 | 1.73 | 159.83 | 4.44 |
| Bosque de Pino | 1 041.81 | 4.17 | 866.87 | 3.33 | 1 908.68 | 7.65 |
| Bosque de Pino-Encino | 192.53 | 0.73 | 867.71 | 3.28 | 1 060.24 | 4.04 |
| Bosque de Encino | 411.41 | 2.04 | 504.05 | 2.46 | 915.46 | 4.55 |
| Bosque de Encino-Pino | 369.58 | 2.44 | 396.82 | 2.56 | 766.40 | 5.06 |
| Matorral Desértico Micrófilo | -1.63 | -0.75 | -19.72 | -9.07 | -21.35 | -9.74 |
| Matorral Desértico Rosetófilo | -1 200.23 | -8.78 | -514.02 | -4.12 | -1 714.25 | -12.55 |
| Matorral Submontano | 367.88 | 0.98 | 1 930.94 | 5.11 | 2 298.82 | 6.14 |
| Pastizal | -122.84 | -10.65 | -61.53 | -5.97 | -184.37 | -15.99 |
| Sin Vegetación | -1305.85 | -5.05 | -4 199.19 | -17.10 | -5 505.04 | -21.29 |
∆ = Incremento o perdida de cobertura.
Estos resultados presentan una tendencia similar a los obtenidos en diversos estudios en áreas protegidas o de exclusión, como lo citado por Flórez-Yepes et al. (2017), quienes estimaron una recuperación de 30.44 % para bosques de minas de oro ubicadas en Manizales, Colombia; así como, lo documentado por Sanhouse-Garcia et al. (2016) en Sinaloa, México, cuyos bosques de pino-encino y pino aumentaron 38.98 % y 29.56 %, respectivamente.
En relación a la agricultura se observaron aumentos de 9.76 % (2000-2010), 9.65 % (2010-2018) y 20.36 % (2000-2018); valores similares a los consignados por Kumar et al. (2020) en un análisis de cambio de cobertura en los márgenes del río Ganges en el distrito de Haridwar, India, para los años 1996, 2003, 2010 y 2017 en donde la agricultura tuvo un incremento de 17.32 %; pero inferior al obtenido por Martin et al. (2019) en el corredor de vida silvestre Kwakuchinja de Tanzania, sitio en el cual la agricultura presentó un incremento de 35.6 %.
Respecto a las investigaciones captura de carbono realizadas en diversos tipos de vegetación (Pimienta de la Torre et al., 2007; Rodríguez et al., 2009; Aguirre-Calderón y Jiménez-Pérez, 2011; Montaño et al., 2016), la recuperación de ecosistemas forestales en el PNCM genera una adicionalidad de 86 939.68 Mg C en bosques de pino-encino, 53 648.00 Mg C en encino-pino, 118 338.16 Mg C en pino y 27 585.84 Mg C en el matorral submontano.
La tasa de cambio progresiva de la agricultura y negativa de los matorrales desértico micrófilo, rosetófilo y del pastizal (Figura 7) pone en riesgo la integridad de los ecosistemas, así como la provisión de bienes y servicios como la recarga de cuerpos de agua que abastecen el Área Metropolitana de Monterrey (López e Ixtacuy, 2017).
A = Agricultura; BA = Bosque de Ayarín; BQ = Bosque de Encino; BQP = Bosque de Encino-Pino; BP = Bosque de Pino, BPQ = Bosque de Pino-Encino; MDM = Matorral Desértico Micrófilo; MDR = Matorral Desértico Rosetófilo, MS = Matorral Submontano; P = Pastizal; SV = Sin Vegetación.
Figura 7 Dinámica de la tasa de cambio del suelo y vegetación.
Uno de los principales factores que continúa afectando a los ecosistemas son los incendios forestales, debido a que en el PNCM se han afectado a más de 9 085.09 ha, durante el período de 2000 a 2011 (Conanp, 2011) y 2 493.68 ha, de 2012 a 2019 (Conafor, 2020).
Los incendios forestales están estrechamente relacionados con periodos previos de excesiva precipitación, como los registrados en el 2008 cuando se dañaron 4 249.31 ha por la acumulación de combustible que resultó de la precipitación (1 600 mm) del 2005; y los del 2011 que abarcaron una extensión de 5 037.51 ha (Yerena et al., 2013), incendios que se presentaron después de una precipitación de 1 915 mm ocurrida en 2010 (Figura 8).
Conclusiones
El análisis multitemporal indica un proceso gradual de recuperación de la cobertura forestal, principalmente, bosques de pino, pino-encino, encino-pino, encino y ayarín; como respuesta de la capacidad de resiliencia y la priorización de estrategias de restauración ecológica en dichos ecosistemas. Una situación contraria se observa para los matorrales desértico rosetófilo, micrófilo y en el pastizal que presentan una pérdida continua, como efecto del cambio de uso del suelo para fines agrícolas.
Las tasas de cambio de uso del suelo y vegetación del presente estudio indican que existe una dinámica en el proceso sucesional de los ecosistemas del PNCM, con una tendencia similar a lo citado en áreas naturales protegidas o de exclusión, en donde existe una recuperación gradual de la cobertura forestal en ciertos tipos de vegetación y pérdida de ecosistemas de menor interés económico, los cuales son sustituidos por el avance de la agricultura y de los asentamientos humanos.
La exclusión de actividades productivas en áreas naturales protegidas no garantiza la conservación de la biodiversidad ni el funcionamiento óptimo de los servicios ambientales, debido a que existen factores externos que ponen en riesgo la integridad de los ecosistemas, como la presencia recurrente de incendios forestales que merman la capacidad de resiliencia y aumentan la vulnerabilidad de los servicios ecosistémicos.
El uso de las imágenes de alta resolución mejora en gran medida la interpretación de los diferentes escenarios de cambio del uso del suelo y vegetación, en comparación con las imágenes Landsat y Sentinel, lo cual permite tener un panorama real de la pérdida o ganancia de la cobertura forestal en áreas naturales protegidas.
