2.1 Descripción del sitio

La zona de protección ambiental PALE, decretada en el año 2011 por el gobierno del Estado de Michoacán de Ocampo; se localiza al Sureste de la capital del Estado, entre las coordenadas geográficas: 19° 39' 48.522'' y 19° 30' 31.325'' de latitud Norte y los 101° 12'24.765'' y 100° 56' 48.567'' de longitud Oeste, comparte territorio de los municipios Morelia, Charo y Madero (Figura 1). Esta ANP, es parte del Sistema Estatal de Áreas para la Conservación del Patrimonio Natural de Michoacán, administrada por la Secretaría de Medio Ambiente, Cambio Climático y Recursos Naturales, con un esquema de protección, preservación, conservación y restauración. Es una de las principales ANP del Estado de Michoacán, con 231.07 km2 de extensión. De gran importancia por representar una fuente potencial de recursos hídricos para la región, principalmente para la ciudad de Morelia (Instituto Nacional de Geografía y Estadística ^{(INEGI), 2017}).

Figura 1 Localización del ANP Pico Azul-La Escalera. 

Para el territorio del ANP, ^{García (1998}) señala el predominio de un clima templado C(w2), la temperatura media anual oscila entre 12 y 18 °C, la del mes más frío oscila entre -3 y 18 °C y la del mes más cálido está por debajo de 22 °C. Es un clima subhúmedo, la precipitación anual oscila de 200 a 1800 mm y la precipitación en el mes más seco, es de 0 a 40 mm ^{(Instituto Nacional de Estadística y Geografía (INEGI), 2017}). ^{Mas et al. (2016}) identifican en el área la existencia de dos grandes comunidades vegetales: los bosques de coníferas y de latifoliadas.

2.2 Cartografía de los paisajes

La cartografía de los paisajes f-g se realiza a partir de la clasificación tipológica a escala local. El resultado permite distinguir las unidades de paisaje a partir de la integración de los componentes naturales (clima, relieve, litología, hidrología, edafología y vegetación) e identificar en ellas las modificaciones antrópicas presentes. Ello se obtiene por el cruce del mapa de las unidades de paisajes f-g y el mapa de coberturas y uso de suelo actual. Así, se determinan, clasifican y cartografían los complejos físico-geográficos o paisajes de carácter tipológico naturales o bien los modificados por las actividades antrópicas (^{Mateo, 2011}).

La clasificación de los paisajes a nivel local (1: 50 000) se conforma por tres niveles de unidades taxonómicas: localidad, comarca y subcomarca f-g (^{Mateo 2011}). Las unidades se delimitan a partir de una taxonomía estructurada en una ordenación jerárquica y anidada, en la que los paisajes se identifican sobre la base de elementos indicadores y factores diferenciadores que determinan el comportamiento y distribución de los componentes naturales, y dan lugar a áreas relativamente homogéneas y repetibles con rasgos comunes (^{Bocco et al., 2010}).

La cartografía se realiza combinando métodos deductivos e inductivos. El primer paso es delimitar las geoformas (mesoformas del relieve). Estas unidades son originadas por la desmembración erosiva del relieve y constituyen el índice diagnóstico primario entre los establecidos para cada nivel de paisaje. Para su diferenciación se tomó la metodología propuesta por ^{García et al. (2022}), la cual consiste en identificar formas del terreno con el software SAGA para clasificar el relieve a partir de los patrones espaciales conformados por la asociación de los polígonos obtenidos semiautomáticamente (Figura 2).

Figura 2 Esquema metodológico para la delimitación y cartografía de las unidades de paisajes, según ^{García et al. (2022}). 

El procedimiento con la herramienta TPI-BLC del SAGA consiste en sintetizar la información geomorfológica del área, la selección de los radios, la interpretación de las formas del terreno y la delimitación y clasificación de las geoformas. Según los índices diagnósticos de la clasificación tipológica, las mesoformas del relieve definen las comarcas f-g. Son las primeras que se identifican y clasifican a partir de la asociación espacial de las formas del terreno, continuación, se sintetiza la información relativa a los índices diagnósticos: régimen de humedad, litología, suelos y vegetación. Las localidades-g se cartografían por asociación de las comarcas, diferenciándose por los datos de las asociaciones litológicas y la génesis del relieve. Las subcomarcas f-g, se identifican cuando se distinguen los elementos o partes de las mesoformas que definieron a las comarcas, para ello se utilizan las categorías generadas por TPI-BLC al interior de las mesoformas.

Los valles reconocidos en el procedimiento anterior, en su mayoría constituyen valles intramontanos, por extensión y ubicación, correspondientes a las comarcas. Para distinguir las subcomarcas de los cauces y fondos de valles de las corrientes fluviales, se recurrió a la herramienta SIG Valley Bottom Extraction Tool (V-BET), que identifica los cauces de inundación y fondos de valles de acuerdo con el principio de que la disposición de las pendientes determina los límites del fondo de las corrientes fluviales. Este procedimiento generó una capa de polígonos que se usó a modo de completamiento o corrección de los límites de los elementos hídricos generados por TPI-BLC.

Con los elementos de las mesoformas delimitados (base espacial de las subcomarcas) y su correspondiente clasificación geomorfológica, se sintetizó la información relativa a los índices diagnósticos: litología, suelo y vegetación. El régimen de humedecimiento se adquirió de las comarcas y localidades; atribuidos a cada unidad según el principio y orden de predominio. Con la finalidad de ajustar los contornos de los elementos de las mesoformas como base espacial para las subcomarcas, se toma como referencia el modelo de sombreado generado a partir del modelo digital de elevación, imágenes satelitales y perfiles topográficos.

El siguiente paso es, a partir del cruce del mapa de paisajes f-g con el mapa de cobertura y uso del suelo actualizado, generar el mapa de los paisajes antroponaturales, dichas unidades ilustran de forma general el panorama del área de estudio, en cuanto a los diferentes usos y su distribución en los paisajes del territorio.

2.3 Evaluación del grado de transformación antropogénica de la cobertura vegetal de los paisajes

El grado de transformación antropogénica que presenta la cobertura vegetal de los paisajes f-g, se evalúa con el índice de antropización de la cobertura vegetal (IACV) de ^{Shishenko (1988}). Este índice ha sido aplicado en diversas regiones de América (^{Aguirre et al., 2017}; ^{Bollo y Velasco, 2018}; ^{Cruz-Paz et al., 2018}; ^{Morales et al., 2019}; ^{Priego et al., 2004}; ^{Ramírez 2013}), bajo la Ecuación 1:

(E. 1)

Donde:

IACV: Índice de antropización de la cobertura vegetal del paisaje. ri: Grado de transformación antropogénica de los paisajes del tipo i de utilización. A _ij : Área dedicada al tipo de utilización i en el geocomplejo j. A _j : Área de la unidad de paisaje j.

Para el cálculo del índice, la base constituye el mapa paisajes antroponaturales, con el marco establecido por ^{Shishenko (1988}). La metodología propuesta por el citado autor asigna a cada uso de suelo o cobertura vegetal contenidos en un polígono de paisaje un valor ponderado o peso relativo (0.01 a 0.99) según el grado de transformación que representa esa cobertura o uso. Se parte del principio de que si existe cobertura hay menor transformación antropogénica que cuando existe un uso, así se otorga el valor más bajo de ponderación a las coberturas menos degradadas y el valor mayor para los usos que originan una mayor intensidad de la antropización de la cubertura vegetal. Estos valores de ponderación permiten dar un peso relativo al proceso de antropización, diferenciando entre la vegetación natural (primaria), la vegetación antropizada (secundaria con predominio de arbustos o herbáceas), la vegetación cultural (cultivos y pastizales), y usos intensivos como la urbanización (Cuadro 1).

Para el cálculo del indicador en cada polígono de paisaje, el valor de ponderación de cada cobertura o uso de suelo se multiplica por la superficie que ocupa, se realiza la sumatoria de los productos obtenidos y el resultado se divide entre el área total del polígono de la comarca. El valor del índice se expresa de cero a uno, donde un valor cercano a cero expresa paisajes con baja modificación antropogénica y un valor cercano a 1, paisajes con niveles altos de transformación. A partir de los resultados obtenidos se elabora un cartograma en ArcGIS por el método de rompimiento natural; que agrupa los valores similares y maximiza las diferencias entre clases (^{Jenks, 1967}). Como resultado, los valores se dividen en clases cuyos límites establecen dónde hay diferencias considerables entre los valores de los datos. Se establecen cinco categorías de antropización: muy alta, alta, media, baja y muy baja.

3.1 Paisajes físico-geográficos del área natural protegida

Los paisajes del ANP se cartografiaron a escala 1: 50 000 (Figura 3), quedaron conformados por 2 390 polígonos, distribuidos al interior de dos grandes regiones f-g descritas por ^{Bollo et al. (2019}): la Provincia f-g Eje Neovolcánico Transmexicano, subprovincia Mil Cumbres, en la cual se identificaron cuatro localidades, 18 comarcas y 95 subcomarcas f-g; y en la Provincia Sierra Madre del Sur, subprovincia Depresión del Balsas, en la cual se identificó una localidad, seis comarcas y 19 subcomarcas. El Apéndice 1 muestra la leyenda resumida del mapa de paisajes.

Figura 3 Paisajes físico-geográficos del ANP Pico Azul-La Escalera. 

El mapa de paisajes antroponaturales, que es el resultado final de este procedimiento y la base para la evaluación de la antropización de la cobertura vegetal de los paisajes, se obtuvo al cruzar el mapa de paisajes f-g, con el mapa base de vegetación y uso de suelo, elaborado a partir del procesamiento de la imagen Landsat LC08 L1TP 027046 20210226 01 T1, del 26 de febrero del año 2021, apoyado en la clasificación de la vegetación realizada en el territorio por ^{Mas et al. (2016}). Como resultado se obtuvo que la cobertura principal, constituida por los bosques de pino-encino, presenta un 43.02 % de su superficie considerada como vegetación primaria, un 10.07 % por vegetación secundaria arbustiva y el 46.91 % por vegetación secundaria herbácea. En la cobertura de bosques de pinos, el 57.18 % es vegetación primaria, el 9.35 % es vegetación secundaria arbustiva y el 33.46 % bosques con vegetación secundaria herbácea. En la cobertura de bosques de encino, el 28.38 % es vegetación primaria, el 24.06 % vegetación secundaria arbustiva y el 47.54 % bosques con vegetación secundaria herbácea. La cobertura de selva baja caducifolia ha desaparecido por completo como vegetación primaria, el 61.61 % es vegetación secundaria arbustiva y el 38.43 % restante, vegetación secundaria herbácea.

3.2 Antropización de los paisajes del área natural protegida

Acorde con el cálculo del IACV, se obtuvieron las siguientes definiciones y rangos del grado de antropización de la cobertura vegetal de los paisajes mediante el método de rompimiento natural en ArcGIS (Cuadro 2).

En el ANP, los paisajes mejor conservados en la categoría muy bajo ocupan tan solo el 8.6 % del territorio; mientras que los paisajes en la categoría de bajo IACV son los más extendidos, en el 42.13 %. Ello muestra que el ANP presenta un grado de conservación de la cobertura vegetal aceptable, con el 50 % de sus paisajes en buen estado de conservación. Los paisajes modificados, con categorías medias y altas del índice ocupan la superficie restante, con excepción de un 2.5 % del ANP con categoría muy alta, que presenta fuertes modificaciones antropogénicas (Cuadro 3).

En el Cuadro 4 se muestra la diferenciación del índice por localidades y comarcas f-g del ANP.

En el Cuadro 4 se puede observar la distribución de las categorías por localidades de paisajes f-g. En las localidades f-g I y II, las de mayor extensión, predomina la categoría baja del índice, que representa el 20 % de la superficie del ANP en cada una de ellas, esta categoría de conservación indica baja modificación (Cuadro 2). La distribución se puede observar en la Figura 4 para las 24 comarcas presentadas en el mapa de paisajes y calculado en los 136 polígonos que las integran.

Figura 4 Distribución de las categorías de IACV en las comarcas físico-geográficas del ANP. 

La categoría muy baja del IACV tiene superficies similares en las localidades f-g I a la IV, presente en tan solo 6 comarcas y 24 polígonos, de pequeña extensión, estas mismas deben ser observadas por ser las de muy poca modificación. La categoría media del índice solo se presenta en las localidades f-g II y V, en 10 comarcas y 24 polígonos con superficies extensas en su mayoría, representa territorios donde la modificación comienza a manifestarse de manera importante, es una categoría que merece atención en el sentido de tratar de detener los procesos que provocan la antropización. La categoría alta sigue la misma distribución que la categoría media, en las localidades f-g II y V con extensiones similares, se presenta en 9 localidades y 26 polígonos, en ellas las modificaciones ya son intensas, necesitan de medidas de restauración de la vegetación original, son paisajes donde se realizan las principales actividades económicas del ANP (Cuadro 2). Los paisajes muy modificados, que se corresponden con la categoría muy alta del índice, se presentan solamente en la localidad f-g I, en 1 comarca y un polígono, que se encuentra urbanizado y dedicado a actividades agropecuarias.

Las comarcas f-g con muy baja modificación en el ANP (Figura 5), son la comarca 4, de la localidad f-g I, en los piedemontes de basaltos, con suelo andosol húmico, y coberturas de bosques primarios y en menor medida secundario con predominio de arbustos, de pino encino, de encinos y mesófilo de montaña, tienen la superficie más extensa de paisajes poco modificados; así como las comarcas 6 y 7, en conos y domos volcánicos, sobre andesita y brecha volcánica andesítica, con suelos acrisol órtico, andosol húmico y coberturas primarias y algunas secundarias de bosques de pino encino, de encinos y mesófilo de montaña secundario con predominio de arbustos. Las comarcas de esta categoría se caracterizan por fuertes pendientes, suelos de baja calidad para actividades agropecuarias, ausencia o muy poca red fluvial de carácter temporal y carencias de viales de acceso, lo que se refleja en su baja modificación antrópica.

Figura 5 Comarcas con muy baja modificación (6 y 7). 

Las comarcas 1,3 y 5 de la localidad I presentan, junto a las comarcas 6, 7, 8, 9 y 13 de la localidad f-g II, las más importantes superficies de paisajes con baja modificación (Figura 6). Las primeras son parte de un cono volcánico, domo y piedemontes con pendientes de (0-30o) sobre andesita y brecha volcánica andesítica, con suelos andosol húmico y luvisol crómico, y coberturas de bosques primarios de pino encino, de encino y mesófilo de montaña secundario con predominio de arbustos y herbáceas; en estas comarcas, se presentan algunas áreas de cultivos perennes y de agricultura de temporal.

En las comarcas 6, 7, 8, 9 y 13 de la localidad f-g II, que conforman el cono y piedemonte con pendientes de (0-30o), sobre rocas ígneas extrusivas ácidas con suelos andosol húmico y acrisol órtico, y coberturas de bosques primarios y secundarios arbustivos de pino encino y de encino, se localizan algunas áreas de cultivos perennes; y en la depresión vulcano-erosiva, sobre andesita y brecha volcánica andesítica, con suelo luvisol crómico, andosol húmico y acrisol órtico, predominan los bosques primarios de pino encino, de encinos y mesófilos de montaña y bosques secundarios de estos tipos con predominio de arbustos.

Figura 6 Comarcas con baja modificación antrópica (1, 3 y 5). 

Las comarcas de esta categoría se caracterizan por el predominio de fuertes pendientes, suelos de calidad baja a media para actividades agropecuarias, ausencia o muy poca red fluvial de carácter permanente y muy pocas vialidades de acceso, lo que se refleja en su baja modificación antrópica.

Las comarcas con categoría media del índice muestran modificaciones y transformaciones variables (Figura 7), las mayores superficies las encontramos en los polígonos de las comarcas 7, 8, 10 y 13, de la localidad f-g II, ello indica un proceso de ocupación por diferentes actividades agropecuarias y la construcción de viviendas rurales que conlleva al inicio y extensión del proceso de antropización.

Figura 7 Comarcas con categoría media del índice de modificación (7, 8, 10 y 13). 

La categoría media del índice de modificación la encontramos también en la localidad f-g V, donde el proceso de antropización ya se ha extendido al 50 % de la superficie de esta. Está representada por polígonos de las comarcas 19, 21, 22, 23 y 24, conformadas por paisajes de valles fluvio-denudativos-erosivos, encajados, y cauces con corrientes fluviales perennes e intermitentes, en andesita y brecha volcánica andesítica, con suelos luvisol crómico, feozem háplico y coberturas secundarias con predominio de arbustivas y herbáceas de selva baja caducifolia, bosques de pino encino y encino; son notables en estas comarcas las áreas de bosques convertidas en pastizales inducidos y, en menor medida, las convertidas al cultivo del aguacate. Esta categoría es característica en la localidad f-g II en los paisajes contiguos a la localidad V, desde los cuales se han extendido la ocupación del territorio a través del desmonte de los bosques para la creación de pastizales inducidos y la siembra de aguacate.

Las comarcas f-g en esta categoría media se caracterizan, generalmente, por pendientes medias (0-30o), suelos de calidad media para activades agropecuarias, una red fluvial con predominio de corrientes permanentes y poca presencia de viales de acceso, lo que se refleja en una categoría medio de antropización; las principales acciones modificadoras son la tala de los bosques para su aprovechamiento, la introducción de cultivos permanentes de aguacate y la creación de pastizales inducidos.

Las comarcas con un alto índice de modificación antrópica se distribuyen fundamentalmente al suroeste del territorio (Figura 8), son representativos los polígonos de las comarcas 6, 7, 8, 11 y 13 de la localidad f-g II y la comarca 14 y 15 de la localidad III, que van desde el cono y domo volcánico hasta los piedemontes, con pendientes de planas a escarpadas (0-45o), formadas por andesita-brecha volcánica andesítica, basaltos, rocas ígneas extrusivas ácidas, limolitas y conglomerados, con suelos acrisol órtico, andosol húmico y ranker, y coberturas de bosques secundarios con predominio de arbustos y herbazales, de bosque mesófilo de montaña, de pino encino y de encinos; extensas áreas de pastizales inducidos fundamentalmente por la tala de los bosques de pino encino y áreas dedicadas a la agricultura de temporal.

Figura 8 Comarcas con alto índice de modificación antrópica (6 y 11). 

Al suroeste del territorio con un alto índice de modificación antrópica, se encuentran polígonos de las comarcas 19, 20, 21 y 22 de la localidad f-g V, formadas por montañas vulcano-denudativas con pendientes de planas a fuertemente inclinadas (0-30o) y valles fluvio-denudativos, erosivos encajados, en forma de V, con laderas de ligeramente a muy fuertemente inclinadas (5-45o), y cauces con corrientes fluviales permanentes e intermitentes, constituidos por andesita-brecha volcánica andesítica y conglomerados, con suelos luvisol crómico, feozem háplico y regosol éutrico, coberturas de bosques secundarios con predominio de herbáceas de pino encino, de encinos y selva baja caducifolia; así como extensas áreas de pastizales inducidos y de áreas dedicadas a la agricultura de temporal y en menor medida permanente.

Las comarcas de esta categoría se caracterizan, de forma general, por pendientes medias, suelos de buena calidad para activades agropecuarias, una red fluvial con predominio de corrientes fluviales permanentes, presencia de asentamientos y viales de acceso, lo que se refleja en un alto índice de antropización inducido por la eliminación de los bosques de pino encino y encino fundamentalmente, para áreas dedicadas a la agricultura de temporal y cultivos perennes. Las comarcas con un muy alto índice de modificación antrópica se localizan al oeste del ANP (Figura 9), corresponden a la totalidad del polígono de la comarca 2 de la localidad f-g I, caracterizada por superficies planas a ligeramente inclinadas (0-5o), con remanentes de coberturas de bosques secundarios de pino encino con predominio de herbáceas y el predominio de extensas áreas dedicadas a la agricultura de temporal, permanente y pastizales. Se encuentran 7 polígonos de la comarca 22 en la localidad f-g V, formada por valles fluvio-denudativos, erosivos encajados, en forma de V, con laderas de ligeramente a muy fuertemente inclinadas (5-45o), y cauces con corrientes fluviales intermitentes, constituidos por andesita-brecha volcánica andesítica, con suelos luvisol crómico, feozem háplico y regosol éutrico, y coberturas de remanentes de bosques secundarios con predominio de herbáceas de pino encino y selva baja caducifolia y el predominio de pastizales inducidos.

En su mayoría son comarcas con pendientes planas (0-15o), donde el bosque natural ha sido fuertemente afectado y en parte remplazado por pastizales inducidos y áreas agrícolas de cultivos de temporal y perennes, fundamentalmente aguacate; presentan poblaciones importantes, con alta accesibilidad y presencia de cursos de aguas de carácter temporal y permanente.

Figura 9 Comarcas muy fuertemente modificadas (2). 

La intensidad de los procesos de antropización, expresados a través de las categorías del IACV, muestran una relativa correspondencia espacial con la composición geomorfológica de los paisajes, relación más evidente en las comarcas f-g que se corresponden a las geoformas con menores pendientes o cercanas a valles y cauces fluviales, con suelos aptos para el cultivo; en estas geoformas es mayor la accesibilidad, por lo que se presenta mayor número de vialidades que facilitan los procesos de antropogenización; hay también una correspondencia con la accesibilidad a las unidades de paisaje, lo que facilita un mayor desarrollo de infraestructuras viales que posibilitan la apropiación del territorio y la tala de sus bosques. Coincide que la mayoría de las comarcas con fuertes pendientes y suelos poco fértiles presentan poca a muy poca antropización, fenómeno lógico y frecuente en los territorios montañosos del país.