Se analizó la dinámica espacial y composición de la cobertura vegetal del Área Prioritaria para la Conservación de Pastizales (apcp) Cuchillas de la Zarca, determinando su afectación por los procesos de fragmentación. Se generó una clasificación supervisada y un análisis de fragmentación del periodo 2008 a 2011. El análisis espectral permitió definir con claridad seis clases, donde el ecosistema mayormente representado es el de pastizal natural con un 41.6% del área total. El análisis de fragmentación de las métricas del paisaje generadas muestra un proceso gradual de fragmentación. Un aumento en el número de fragmentos totales en función del tiempo confirma el proceso de fragmentación en la zona, mientras que para el tamaño medio de los parches se observa una disminución, lo que indica que los parches grandes fueron fragmentados o divididos. Aunado a lo anterior, los índices de Shannon y Simpson presentan una tendencia ascendente (0.66 a 0.89; 0.68 a 0.82), indicando el aumento de fragmentos en el tiempo. La exactitud global de las clasificaciones evaluadas en el periodo 2008-2011 muestra una alta precisión que fluctuó del 91 al 94%, mientras que el índice discreto multivariado Kappa presentó una variación de 0.90 al 0.93.
We analyzed the spatial dynamics and composition of the vegetation cover of the Priority Area for Grasslands Conservation (pagc) Cuchillas de la Zarca, determining their involvement in the fragmentation processes. A supervised classification and fragmentation analysis in the period from 2008 to 2011 was generated. Spectral analysis allowed to clearly define six classes, where the ecosystem is represented mostly by natural grassland, 41.6% of the total area. The fragmentation analysis of the generated landscape metrics shows a gradual process of fragmentation. An increase in the total number of fragments versus time confirms the fragmentation process in the area, whereas the average size of the patches shows decreasing, indicating that large patches were fragmented or divided. In addition to this, the Shannon and Simpson indicators are trending upward (0.66 to 0.89), (0.68 to 0.82), indicating the increase in time pieces. The overall accuracy of the classifications evaluated in the period 2008 - 2011 shows a high accuracy that ranged from 91 to 94%, while the index Kappa discrete multivariate presents a variation of 0.90 to 0.93.
El campo de estudio de la percepción remota ha crecido rápidamente en los últimos 15 años debido a la influencia persistente de la historia del uso del suelo y la alteración de los ecosistemas naturales, estos estudios se enfocan en las interacciones recíprocas entre patrones espaciales del terreno y los procesos ecológicos (McGarigal y Cushman, 2005; Turner, 2005). La modelización ha sido objeto de estudio de un gran número de investigadores que buscan explicar las tendencias espaciales y temporales de los atributos de la vegetación en relación con la estructura y fragmentación de los ecosistemas terrestres (Cohen y Goward, 2004).
En ese contexto, el uso de los sensores remotos, que generan imágenes satelitales, representan un método de bajo costo y permiten evaluar grandes extensiones de terrenos o zonas remotas e inaccesibles (Poso y Paananen, 1984; Wulder 1998; Wynne et al., 2000).
La aplicación de estas técnicas proporciona una base fuerte en el empleo de metodologías que permiten monitorear los patrones que intervienen en los procesos de fragmentación, mismos que están fuertemente ligados a la dinámica espacial y temporal de actividades antrópicas. Recientemente se han realizado estudios enfocados a la aplicación de métricas del paisaje con la finalidad de cuantificar la heterogeneidad a nivel ecosistema (Hoechstetter et al., 2008; McGarigal et al., 2009).
A su vez, los estudios sobre cambios progresivos en la cubierta vegetal a menudo emplean indicadores para cuantificarlos a través del tiempo, utilizando los análisis de los mapas de cobertura terrestre derivados de la clasificación de imágenes de sensores remotos a partir de dos o más periodos de tiempo (Linke et al., 2009; Velázquez et al., 2010).
Esa cuantificación produce información numé-rica denominada métricas (Subirós et al., 2006), que arrojan información sobre las características de la cubierta vegetal y los procesos ecológicos que ahí se presentan, tanto del punto de vista espacial como temporal (Burel y Baudry, 2002).
Con base en lo anterior, el ecosistema pastizal cubre entre 45 y 56millones dekm2 de la superficie terrestre (wri, 2000), y representa uno de los más extensos, mayormente productivos, diversos y amenazados biomas terrestres (Sala et al., 2000; Olson et al., 2001). El cambio de uso de suelo para ampliar la frontera agrícola, derivado de una mayor demanda de alimentos, ha resultado en la modificación de los límites geográficos de las principales regiones biogeográficas en el mundo. Esta modificación incluye al Desierto Chihuahuense, donde la fragmentación de la cubierta vegetal por acciones antropogénicas, como el pastoreo a gran escala y la conversión de áreas grandes en pequeñas parcelas, ha transformado sus frontreras (Curtin et al., 2002; Pinto, 2006; Hoth, 2009; Sánchez et al., 2009; pmarp, 2012).
El efecto principal de estas alteraciones es la fragmentación, lo cual deriva en el corto o mediano plazo en el cambio y/o la pérdida de ecosistemas y por ende de la biodiversidad (Sala et al., 2000).
En este sentido, el área de conservación y protección de pastizales de Norteamérica, región a la que pertenece Cuchillas de la Zarca, ha representado el hábitat invernal exclusivo en pastizal para mamíferos y aves reconocidas como especies focales o especies indicadoras, de importancia continental, que están disminuyendo sus poblaciones de manera sostenida anualmente (Macías et al., 2011; Martínez et al., 2011a; Macías et al., 2013) por lo que se han desarrollado esfuerzos trinacionales (Canadá, Estados Unidos y México) para la conservación de estos ecosistemas y su biodiversidad (cec y tnc, 2005; ecopad, 2007; Berlanga et al., 2010).
El objetivo de este trabajo fue realizar un análisis de fragmentación de la cubierta vegetal mediante técnicas de percepción remota en el Área Prioritaria para la Conservación de Pastizales (apcp) denominada Cuchillas de la Zarca, mediante el uso de imá-genes satelitales Landsat TM y la obtención de mé- tricas del paisaje de la zona, que permite describir los cambios que ha seguido un territorio a lo largo de un periodo de tiempo de cuatro años, donde los resultados ayudan a interpretar la situación y/o su dinámica espacial y temporal desde un punto de vista ecológico.
Área en estudioEl Área Prioritaria para la Conservación de Pastizales (apcp) denominada Cuchillas de la Zarca, forma parte del Desierto Chihuahuense en México y comprende una pequeña parte de la zona sur del estado de Chihuahua y el centro-norte del estado de Durango, cuenta con superficie de 11 588km2, y se ubica entre las coordenadas geográficas extremas 25°20’ 00” y 27°00’ 00” de latitud Norte y 105°50’ 00” y 104°25’ 00” de longitud Oeste (Figura 1).
En la zona se presentan tres tipos de clima, el primero es el semiárido-templado (BS1kw), con temperatura media anual entre 12 y 18°C, temperatura del mes más frío entre -3 y 18°C, temperatura del mes más caliente menor de 22°C; lluvias de verano del 5 al 10.2% anual, respectivamente. El segundo, árido-templado (BSokw), presenta temperaturas entre 12 y 18°C, temperatura del mes más frío entre -3 y 18°C, temperatura del mes más caliente menor de 22°C; lluvias de verano del 5 al 10.2% anual, respectivamente. El último, muy árido-semicálido (BWhw), se dan temperaturas entre 18 y 22°C, temperatura del mes más frío menor de 18°C, temperatura del mes más caliente mayor de 22°C y lluvias de verano del 5 al 10.2% anual.
La vegetación en el área en estudio se compone de matorral xerófilo y pastizal natural (Rzedowski, 1981). El estrato arbustivo está conformado por especies de los géneros Larrea,Prosopis, Acacia, Ephedra, Nolina, Opuntía, Juniperus y Quercus; el estrato herbáceo por los géneros Boutelotia, Aristida, Buchloe, Andropogon, Melinis, Muhlenbergia, Sporobolus, Heteropogon y Pleuraphis.
Se considera como región prioritaria para la conservación, porque constituye el extremo sur de los pastizales asociados a la vertiente oriental de la Sierra Madre Occidental en los estados de Chihuahua y Durango. La diversidad de origen de los suelos contribuye a la formación de los pastizales con diferente composición de especies. Existe además una alta riqueza de aves. La región se delimita hacia el oeste y el norte por subcuencas donde el parteaguas llega hasta la cota de los 2 000 msnm, que coincide con bosques bajos-abiertos de encinos, y hacia el este y sur con pastizales naturales. La región protege un gradiente entre los pastizales naturales, matorrales xerófilos y las zonas ecotonales de la Sierra Madre Occidental (piedemonte) como son los chaparrales y los bosques bajos-abiertos (Martínez et al., 2011b).
MétodosFuentes de datos y procesamiento de imágenes de satéliteSe adquirieron ocho escenas del sensor Landsat TM 5 (path/row31-41 y 31-42) con fechas de captura del 2 al 24 febrero, correspondiendo al periodo de adquisición del 2008, 2009, 2010 y 2011; las escenas fueron delimitadas usando la poligonal para el apcp Cuchillas de la Zarca definida por la Commission for Environmental Cooperation and The Nature Conservancy (cec-tnc, 2005).
Las imágenes fueron registradas usando los sistemas de referencia de la proyección Universal Transversal de Mercator (utm) ubicándose a la zona del huso horario 13N y con un datum y elipsoide World Geodesic System 1984. Las imá-genes poseen una resolución espacial de 30m. El procesamiento digital de las imágenes de satélite, implicó realizar operaciones básicas relacionadas con las correcciones y calibración de las mismas para obtener una representación de la superficie terrestre tan real como sea posible (Eastman, 2003).
Se aplicó un preprocesamiento de calibración de las escenas Landsat TM 5; al comparar escenas de di- ferentes periodos se obtienen mejores resultados si se lleva a cabo una transformación a valores de reflectancia. Las escenas fueron calibradas de acuerdo con los valores de radianza mínima y máxima (Lmax y Lmin) derivado de los valores publicados por Chander et al. (2009). Estos cálculos permiten transformar los valores originales de cada imagen en parámetros físicos, permitiendo así que sean comparables en el tiempo, aun bajo condiciones atmosféricas diferentes (Chuvieco, 2008). La radiancia espectral fue calculada usando la siguiente ecuación:
donde:
qcal: es la calibración y la radiación cuantificadas a escala en unidades de números digitales,
linλ: es la radiancia espectral en qcal=0,
lmaxλ: es la radiancia espectral en qcal=qcalmax;
qcalmin: es el valor mínimo cuantificado de pixel calibrado (correspondiente a lminλ) en Niveles Digitales,
qcalmax: es el valor máximo cuantificado de pixel calibrado (correspondiente a lmaxλ) en Niveles Digitales.
La reflectancia exoatmosférica (ρp) fue calculada usando la siguiente ecuación:
donde:
Lλ: es la radiancia espectral,
d: es la distancia en unidades astronómicas Tierra-Sol,
esunλ: es la media de irradiación solar exoatmosférica,
Φs: es el ángulo cenital solar en grados.
Las escenas calibradas fueron unidas para construir un mosaico de cada banda espectral de cada escena del área, esto permitió obtener el área en estudio empleando una extracción por máscara.
Clasificación supervisadaSe aplicó una clasificación de máxima probabilidad basada en los estadísticos de cada clase, este método asume que todas las clases en cada banda del sensor multiespectral tienen una distribución normal y calcula la probabilidad de que el pixel puede pertenecer a una clase específica basada en la información que éste posee (Richards, 1999).
donde:
i=clase,
x=n-datos dimensionales (donde n es el número de bandas),
|Σi|=determinante de la matriz de covarianza de los datos de la clase
∑i-1 = matriz inversa,
mi=vector.
Este procedimiento es de los más sofisticados, y sin duda el clasificador más utilizado en la clasificación de imágenes de teledetección.
Validación de la informaciónUn procedimiento de verificación permitió cuantificar el error de la clasificación de coberturas bajo procesos de fragmentación. El diseño de muestreo empleado, fue el aleatorio estratificado basado en el entendido de que cada clase temática registra un número determinado de puntos de muestreo de acuerdo con la superficie de cada unidad temática. Para el análisis se generaron áreas de entrenamiento, las cuales fueron utilizadas como zonas de verificación, en el entendido de que se conocía con amplia certeza el tipo de vegetación. Se empleó la información de 96 transectos distribuidos en toda la zona.
donde:
KAPPA=Índice de Kappa,
k=número de filas en la matriz,
Xii=número de observaciones en la fila i y columna i,
Xi+ ;Xi+i=total marginal para la fila i y columna i,
N=número total de observaciones.
Análisis de fragmentaciónUna vez procesadas y generadas las clases de cobertura para las áreas en estudio, se importaron al software ArcGis y se generaron las métricas del paisaje utilizando la extensión de análisis de fragmentación Patch Analysis, el cual tiene el principio de acuerdo con las metodologías sugeridas por Elkie et al. (1999).
Se obtuvieron las siguientes métricas del paisaje para determinar con precisión las áreas fragmentadas de las zonas evaluadas: Número de parches (NumP), Media del tamaño del parche (mps), Media del borde del parche (mpe), Coeficiente de varianza del tamaño del parche (PSCoV), Índice de diversidad de Simpson (sdi) e Índice de diver-sidad de Shannon (sei), donde cabe resaltar a los índices de diversidad, ya que nos dan la pauta para la determinación de proceso de fragmentación de las áreas en estudio (McGarigal y Marks, 1995).
ResultadosEl análisis espectral permitió definir con claridad seis clases (Figura 2), donde el ecosistema mayor-mente representado es el de pastizal con un 61% del área total del apcp, distribuidos entre pastizal natural con un 41.6% y pastizal inducido con un 19.3%. Posteriormente se encuentra el matorral desértico micrófilo con un 17.5% seguido de los bosques de encino-pino con 11.6% y en una menor representación en la zona se encontraron las áreas de chaparral 5.3% y los cuerpos de agua 4.6%. Pinedo (2008) en un estudio similar utilizando imágenes Landsat TM 5, determinó seis clases espectrales (agricultura de temporal, bosque pinoencino, bosque de encino-pino, bosques de encino, comunidades de pastizal y selva baja caducifolia), como las más acertadas para la representación de la vegetación presente en un predio ubicado en el municipio de Guadalupe y Calvo, Chihuahua.
De acuerdo con la Tabla 1, los datos evaluados para los diferentes periodos mostraron valores altos de exactitud, el valor relativamente más bajo lo presentó la clase de pastizal inducido para todos los años evaluados, mientras que el valor más alto lo obtuvo la clase de cuerpos de agua, disminuyendo un poco en el 2009 e incrementando para 2010 y 2011. Wang et al. (2009) evaluaron la precisión de clasificaciones de vegetación generadas del procesamiento de imágenes Landsat TM de diferentes periodos (1976, 1986 y 2002) obteniendo coeficientes de precisión Kappa de 0.87, 0.86 y 0.84. En cuanto a la exactitud global, los valores más bajos los presentaron las clasificaciones de 2008 y 2009, aumentando la precisión para las clases del 2010 y 2011. El índice discreto multivariado Kappa tuvo un comportamiento muy similar a la exactitud global siendo las clasificaciones del 2010 y 2011 con los índices más altos relativamente comparándolos con los demás. Nutini et al. (2013) evaluaron la dinámica de la cubierta vegetal de la región semiárida del Níger utilizando un conjunto de imágenes multitemporales de la serie Landsat de 2000 y 2007, emplearon métodos supervisados en donde la precisión fluctuó entre 0.80 y 0.90. La alta precisión obtenida de acuerdo con la matriz de error, se genera probablemente en las condiciones que imperan en los ecosistemas del semidésierto, asociado a la fecha (febrero) de captura de la imagen, estas superficies muestran diferente absorción de la radiación electromagnética entre clases, en este sentido los CA, por ejemplo, presentaban valores de reflectancia media en las longitudes de onda de los infrarrojos lo que los hace fácilmente identificables.
Resumen de la precisión del periodo de la clasificación de Landsat TM 5 evaluada a través de la matriz de error para el área en estudio (2008 al 2011)
| Clases | Clasificaciones | |||
|---|---|---|---|---|
| 2008 | 2009 | 2010 | 2011 | |
| CA | 0.97 | 0.93 | 0.96 | 0.97 |
| bqp | 0.90 | 0.92 | 0.95 | 0.92 |
| mdm | 0.93 | 0.91 | 0.94 | 0.96 |
| PI | 0.89 | 0.90 | 0.91 | 0.88 |
| AC | 0.93 | 0.88 | 0.92 | 0.93 |
| PN | 0.89 | 0.92 | 0.96 | 0.96 |
| Exactitud global (%) | 91 | 91 | 94 | 93 |
| kappa | 0.90 | 0.91 | 0.93 | 0.92 |
Derivado de la clasificación supervisada se obtuvieron datos de tasa de cambio basados en un periodo de cuatro años, mostrando una disminución gradual de la superficie de las áreas de chaparral con un cambio de -661 ha anuales, para una pérdida 2 644 ha (22.3%) en el periodo analizado, seguido de los bosques de encino-pino con una tasa de cambio anual de -795.6 ha (-6.7%) y por último el pastizal natural con -296.7 ha y una tasa de cambio de -0.35%. Posteriormente se identifica un aumento en la superficie ocupada por el pastizal inducido 3 523.7 ha (9.7%) y del matorral desértico micrófilo con 1 204.6 ha (3.5%). En el caso de los cuerpos de agua se encontró una disminución de -131.5 ha (1.44 %) entre el 2008 al 2010, sin embargo, en 2011 hubo una aumento de 571.33 ha (6.2%; Figura 3). De acuerdo con los datos publicados por el Servicio Meteorológico Nacional de conagua (2013), este fenómeno obedece al comportamiento de los datos mensuales de afectación por sequía en México, donde específicamente en el periodo de lluvias del 2010 hasta principios de 2011, el porcentaje de afectación por sequía en el área en estudio fue bajo, mostrando afectaciones de datos anormalmente secos y sequía moderada, lo cual permitió que los cuerpos de agua temporales se mantuvieran por un periodo de tiempo mayor. Velázquez et al. (2010) encontraron que los cambios de las coberturas del suelo sin disturbio que no evidenciaban cambio al inicio, presentaron una tasa de cambio de 13.7% en relación con las cober- turas de uso de suelo con algún disturbio, las cuales tuvieron una recuperación 6.7%. Caciano et al. (2011) mencionan que de acuerdo con las actividades humanas en el área natural protegida Sierra de Los Lobos, Guanajuato, ha disminuido un 10.5% de superficie agrícola de temporal, de manera contraria, se incrementó la superficie de matorrales y pastizal inducido con una tasa de cambio anual de 10.9%. Por otra parte, en un estudio realizado en el noroeste de Turquía se encontró una disminución del 33.5% de las áreas forestales, mientras que las áreas agrícolas tuvieron un incremento del 11.2%; los cambios se deben a la rápida urbanización y actividades agrícolas (Ikiel et al., 2013).
El análisis y la interpretación de las métricas y el índice de fragmentación que determinan la evolución del paisaje en el apcp Cuchillas de la Zarca, se describe de manera resumida en la Tabla 2.
Valores de las métricas e índices de fragmentación obtenidos en el análisis del paisaje de las clases temáticas evaluadas de 2008 a 2011
| Periodo | NumP | mps | mpe | PSCoV | sdi | sei |
|---|---|---|---|---|---|---|
| 2008 | 245.00 | 110.30 | 5 156.56 | 547.12 | 0.68 | 0.66 |
| 2009 | 210.00 | 116.78 | 6 662.67 | 545.71 | 0.72 | 0.70 |
| 2010 | 381.00 | 49.34 | 4 969.45 | 689.19 | 0.79 | 0.81 |
| 2011 | 1 300.00 | 25.45 | 3 008.78 | 1 787.25 | 0.82 | 0.89 |
Es claro el aumento en el número de fragmentos totales en función del tiempo, confirmando que existe un proceso de fragmentación o división de los grupos homogéneos de las diferentes clases o tipos de vegetación que existían en periodos pasados, se observa el aumento progresivo durante los cuatro años de estudio (Tabla 2). Plue y Cousins (2013) mencionan que en comunidades de pastizal existe una importante asociación con el número de parches, la dispersión temporal permite una prolongada presencia de numerosas especies típicas de pastizales en paisajes fragmentados.
Media del tamaño del parche.En la Tabla 2 se observa la disminución del tamaño promedio de los parches, lo cual indica que los parches grandes fueron fragmentados o divididos debido a diversas causas propias de las actividades productivas desempeñadas en el apcp (Figura 2).
Coeficiente de variación del tamaño medio de los parches.Los resultados muestran que este parámetro tiene un comportamiento similar al de número de parches, mostrando valores bajos en 2008 y con tendencia a la alta conforme se avanza en el tiempo, efectos producidos quizá por el aprovechamiento de recursos naturales, actividades productivas o por fenómenos naturales tales como la sequía que se ha acentuado en los últimos años en la zona. Se puede apreciar el aumento en el coeficiente de variación del tamaño del parche, éste es propiciado por el aumento en el número de fragmento volviendo más heterogéneo el ecosistema (Figura 4).
Media del borde del parche.
Este valor mide los valores de una misma clase (uso del suelo) y muestra el grado de fragmentación y conectividad, los valores obtenidos fluctúan de 0 al infinito. Paisajes con valores reducidos indican que se encuentran más fragmentados y aislados, y paisajes con valores altos indican que se encuentran menos fragmentados. Para el periodo evaluado existe una menor fragmentación en los primeros periodos de estudio, el valor tiende a disminuir específicamente en 2011 lo cual muestra de nuevo que es el año en el cual es más notorio el proceso de fragmentación para el apcp (Figura 5).
Índices de diversidad. De acuerdo con el análisis efectuado a las clasificaciones de las imágenes satelitales para los periodos evaluados, tanto el índice de Simpson así como el de Shannon, tienen una clara tendencia ascendente. El índice de Simpson muestra que conforme avanza el tiempo existe una mayor diversidad paisajística; por otra parte, el índice de Shannon nos da como resultado que del 2008 al 2011 está incrementándose el número de fragmentos de cada clase. En ambos casos se muestra una evolución paisajística debido al incremento de los parches que presentan un claro proceso de fragmentación en el área en estudio (Figura 6). Wu et al. (2013) en su estudio en comunidades de pastizal del noreste de China, encontraron que el número de grupos funcionales de plantas mostraron respuestas distintas debido a la fragmentación; el aumento de la diversidad era producto del aumento de la fragmentación.
De acuerdo con McGarigal et al. (2002), así como con Matteucci et al. (2004), esto nos da una indicación del grado de fragmentación, especialmente si en su estado inicial la región era relativamente homogénea en cuanto a las clases o uso del suelo.
ConclusionesEn términos generales, durante el periodo analizado, en el apcp Cuchillas de la Zarca resulta evidente el proceso de fragmentación de acuerdo con las métricas de paisaje analizadas.
En el caso de los índices de diversidad paisajística, cabe ser muy prudente en la interpretación de los resultados, pues aportan información sobre heterogeneidad y homogeneidad desde un punto de vista estrictamente cuantitativo, por lo que no se pueden extraer valoraciones de carácter cualitativo sobre los resultados obtenidos en un paisaje concreto.
Para las operaciones de mapeo y clasificación, la resolución y nivel de detalle, así como la precisión alcanzada por los sensores evaluados en este estudio, permiten trabajar grandes extensiones de terreno de forma precisa, facilitando los trabajos de cuantificación, monitoreo y comprensión de la dinámica de los ecosistemas.
Los resultados obtenidos serán de gran importancia para otros estudios que se realicen en la zona e identificar la posible relación o efecto del proceso de fragmentación que se está presentando en la zona.