text new page (beta)
Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
Introducción
La erosión es un problema de degradación del suelo que implica impactos ambientales, económicos y sociales. Las cárcavas son formas de erosión avanzada, donde las partículas de suelo desprendidas afectan la calidad del agua, sedimentan cauces y colmatan embalses aguas abajo (Pathak et al., 2005).
Las cárcavas son estrechas incisiones que se forman cuando se concentra el escurrimiento en dirección de la pendiente del terreno y generalmente transportan agua durante y después de un evento de lluvia (Poesen, 1993; Poesen et al., 2003). Las cárcavas se clasifican en efímeras y permanentes, las primeras son aquellas que se borran con el laboreo agrícola mientras que las segundas son difíciles de eliminar con implementos agrícolas.
Las covariables que intervienen en la formación cárcavas, se dividen en dos tipos: covariables antrópicas y físicas. Entre las covariables antrópicas destaca el sobrepastoreo, producto de una elevada carga animal en los ecosistemas (Nyssen et al., 2002); la agricultura sin control de escurrimientos y su expansión en pendientes pronunciadas; y la deforestación por su impacto en la disminución de la cubierta vegetal (Nyssen et al., 2002; Pathak et al., 2005). De acuerdo con Moore et al. (1991), Poesen et al. (2003) y Valcárcel et al. (2003), entre las covariables físicas que influyen en la erosión en cárcavas, se tienen las geológicas, climáticas, edáficas, topográficas, hidrológicas, de uso de suelo, vegetación y manejo de los terrenos. Para las covariables geológicas destaca el material parental (Moore et al., 1991); entre climáticas están aquellas ligadas a la distribución, intensidad y magnitud de los eventos extremos de lluvia (Poesen y Hooke, 1997; Nyssen et al., 2002); en las edáficas, las relacionadas con la distribución espacial de las propiedades del suelo (Lal, 1992; Bull y Kirkby, 1997); topográficas, aquellas vinculadas con la velocidad y concentración del escurrimiento como la forma del terreno, el índice topográfico de humedad y de posición topográfica, la longitud y grado de la pendiente (Kreznor et al., 1989; Zinck y Valenzuela, 1990; Moore et al., 1991; Verstappen y Van Zuidam, 1991; Stolt et al.,1993; Pathak et al., 2005); y las hidrológicas, que dependen del tamaño y forma de la cuenca (Meijerink, 1988; Bocco, 1991).
Entre las herramientas y procedimientos que se disponen para la detención de cárcavas están los medios ópticos, a través de sensores remotos, y los métodos de clasificación digital, a través de sistemas de información geográfica (SIG) y el uso de técnicas de análisis multicriterio. Los sensores remotos son considerados como una herramienta de bajo costo para colectar datos de uso de suelo, densidad de vegetación e identificar áreas con problemas de erosión (Chuvieco, 2002; Bagyaraj et al., 2013). Por otra parte, los SIG son una herramienta que facilita el análisis multicriterio (Basildo-Martín y López, 1998) de diversas capas de información ambiental (covariables) para un tema geográfico de interés. Este análisis, a partir de una mezcla lineal de probabilidades ponderadas, permite identif icar, según la aptitud de cada cobertura temática, la interdependencia entre estas covariables (Bagyaraj et al., 2013). Los análisis multicritero, por ejemplo, son aplicables para analizar las covariables que influyen en la erosión (Tamene et al., 2006).
La clasificación supervisada es una técnica usada en investigaciones ambientales para agrupar información (pixeles) proveniente de diversas fuentes; por ejemplo, en estudios de usos de suelo y vegetación, sitios con riesgo a deslizamientos, áreas con problemas de erosión hídrica, determinación de potencial de recarga de acuíferos, entre otros. La clasificación supervisada de la imagen resultante del análisis de componentes principales, consiste en categorizar estadísticamente los datos de campo según las clases observadas en el área de estudio. El método parte de definir las categorías, observar la distribución de los sitios de muestreo y verificar los resultados (Chuvieco, 1990). Para clasificar la información, que describe un proceso, es necesario conocer el área de trabajo para definir correctamente los polígonos de entrenamiento y los valores estadísticos (media, rango, varianza, desviación estándar, covarianza, etc.) entre covariables.
De acuerdo con Guerrero-Arenas et al. (2010), los ecosistemas de la Mixteca alta sufrieron, desde el Pleistoceno tardío hasta el Holoceno, una serie de transformaciones geológicas y orográficas. Después de la conquista española, el área ha experimentado deforestación, sobrepastoreo con ganado caprino (Mendoza-García, 2002) y expansión de la frontera agrícola, lo que ha provocado amplias superficies de erosión en cárcavas. En vista del estado avanzado de erosión que experimentan los depósitos eólicos y aluviales, de origen volcánico, de la Mixteca, surgió la inquietud de determinar las covariables ambientales que influyen en la formación de cárcavas e identificar áreas susceptibles al avance de este tipo de degradación. El objetivo de la investigación fue identificar áreas susceptibles a la erosión por cárcavas, a través de clasificación supervisada de componentes principales, según las covariables ambientales con mayor influencia en su formación.
Materiales y Métodos
La zona de estudio se localiza al noroeste del estado de Oaxaca y cubre el 39.3% de la Mixteca Oaxaqueña, región conocida a nivel mundial por su importancia biológica y cultural. El área de trabajo abarca una superficie de 6167.5 km2, se ubica en los distritos político-administrativos de Coixtlahuaca (24.28%), Juxtlahuaca (5.93%), Nochixtlán (20.29%), Teposcolula (19.53%) y Tlaxiaco (29.97%), específicamente entre las coordenadas geográficas 16° 55’ 32.46” y 18° 7’ 2.74” N y 97° 58’ 13.35” y 97° 01’ 7.42” O (Figura 1).
La zona de estudio se ubica en la provincia fisiográfica Sierra Madre del Sur (INEGI, 2010) cuyo origen es el choque de la placa tectónica de Cocos con la de Norte América, generando el levantamiento de la sierra y la gran complejidad geológica de sus subprovincias. Así, el área de análisis es cubierta por las subprovincias Mixteca Alta (68.53%), Sierras Centrales de Oaxaca (27.64%) y Cordillera Costera del Sur (3.83%).
De acuerdo con SGM (2000a y 2001), en orden cronológico, en el área de estudio las rocas más antiguas son del Precámbrico y cubren el 4.3% de la superficie. Durante el Paleozoico se originaron litologías de conglomerados poligénicos de lutita-areniscas, areaniscas-lutitas y caliza-lutitas, que cubren el 0.6% del área de estudio. El 31.2% del área analizada es de origen mesozoico y se compone de materiales como caliza, calizas con dolomía, calizas-lutitas, marga-limolita, lutita-caliza, arenisca-lutita, arenisca y toba andesítica. Durante el periodo reciente (era cenozoica), se formó el 63.8% de la litología del área de estudio, con materiales producto de erupciones volcánicas y arrastres aluviales proveniente de toba andesítica-andesita, toba andesítica-arenisca y conglomerados poligénico-arenisca. Estos materiales volcánicos se han acumulado en las partes bajas producto de la erosión producida en las partes más altas.
Según el SGM (2000b y 2001), en el área de estudio se presentan las siguientes formaciones estratigráficas: Tamazulapam (39.72%), Teposcolula (14.76%), Tlaxiaco (5.6%), Chimeco-Sabinal (4.36%), Tecocoyunca (2.18%), Tilaltongo (1.73%), Zapotitlán (0.44%), Matzitzi (0.35%), Ixtaltepec (0.21%), Huajuapan (0.12%), Tecomazuchitl (0.08%), Tiñú (0.03%), Yucunama (0.03%), Yododeñe (0.02%) y sin formación no identif icad (30.37%). La parte norte del área de estudio que limita con las formaciones Teposcolula y Tamazulapam son depósitos que se acumularon en una cuenca sedimentaria formada al lado un anticlinal (Santamaría-Díaz et al., 2008). La formación Tamazulapam, controlada principalmente por la falla del mismo nombre, es donde se alinea la mayor concentración de cárcavas del área de estudio.
El 74.1% del área de estudio presenta un clima templado con lluvias en verano y precipitaciones que varían de 506 a 1136 mm y temperaturas entre 8.1 y 18 °C; el 17.06% del área corresponde a un clima seco con precipitaciones de 316 a 651 mm y temperaturas entre 11.1 y 25.8 °C; y el 8.83% de la superficie analizada presenta un clima semicálido con precipitaciones de 640 a 1136 mm y temperaturas entre 18 y 19.6 °C (INEGI, 2008). Los suelos presentes en el área de trabajo, según la serie II del INEGI (2006), se agrupan en leptosoles (38.18%), luvisoles (25.13%), phaeozem (15.60%), regosoles (11.31%), vertisoles (5.44%), cambisoles (3.66%), gleysoles (0.42%), fluvisoles (0.11%) y el 0.15% es área urbana.
De acuerdo con uso de suelo y vegetación, realizado por WWF (2014), la zona de estudio cuenta con 15 tipos de cubiertas. Los pastizales ocupan el primer lugar (20.99%), con importantes áreas con problema de cárcavas. El segundo lugar corresponde al bosque de encino (19.11%), seguido por bosque de pino encino y encino pino (16.60%), uso agrícola (12.92%), matorral (12.12%), y otras cubiertas con una proporción menor al 10%.
La metodología (Figura 2) inició con una revisión bibliográfica y recopilación de información para integrar el sistema de información geográfica, conformado por el modelo digital de elevaciones (INEGI, 2012), distribución de suelos (INEGI, 2006), geología (SGM, 2000a y 2001), clima (INEGI, 2008), fisiografía (INEGI, 2010), vegetación y uso de suelo (WWF, 2014), población (INEGI, 2014 y 2015), imágenes de satélite y recorridos de campo. La estación receptora del sistema ERMEX proporcionó tres imágenes del satélite Spot 7 (dos del 29 de enero del 2015 y una del 4 de noviembre de 2015) y tres del Spot 6 (del 16 de febrero de 2015). El recorrido de campo se realizó entre el 8 y el 10 de diciembre de 2015 del cual se identificaron 503 puntos de entrenamiento. Se utilizaron imágenes Google Earth Pro 6.0 para generar áreas de entrenamiento y verificar los resultados obtenidos, ya que la resolución de las imágenes Spot multiespectrales son de 6 m y de las primeras es de 1.2 m.
Figura 2: Diagrama metodológico para obtener las covariables que influyen en la formación de cárcavas.
La identificación de covariables ambientales, que influyen en la formación de cárcavas, se realizó con base en la revisión bibliográfica y las correlaciones entre las covariables, de este modo, las covariables seleccionadas para el análisis fueron diez: 1) pendiente del terreno (PEN), acorde a Burrough y McDonnell (1998), se obtuvo con el comando 3D Analyst Tools>Raster Surface>Slope de ArcGis 10; 2) curvatura perfil o vertical (CV) y 3) curvatura plana u horizontal (CH), según Moore et al. (1991), ambas se generaron con el comando 3D Analyst Tools>Raster Surface>Curvatuure de ArcGis 10; 4) longitud de la pendiente (LP), de acuerdo con Moore et al. (1991), se calculó con el comando Terrain Analysis >Hydrology>Slope Lenght de Qgis 2.8.1; 5) índice de Posición topográf ica (IPT), como indica Weiss (2001), se determinó con la extensión para ArcGis Land Facet Corridor Designer; 6) índice topográf ico de humedad (ITH), se calculó en ArcGis 10, basado en lo propuesto por Beven y Kirkby (1979); 7) distancia a picos (DP), corresponde a la distancia euclidiana medida a partir de las crestas generadas con IPT y para su cálculo se usó la herramienta Spatial Analyst Tools>Distance> Euclidean Distance de ArcGis 10; 8) distancia a la falla Tamazulapam (DF), se obtuvo con la herramienta Spatial Analyst Tools>Distance> Euclidean Distance; 9) índice de vegetación de diferencia normalizada (NDVI), como indicador de la densidad de la cobertura vegetal, se obtuvo de imágenes Spot corregidas por radianza y reflectancia a través del software Qgis 2.8.1 con el método al tope de la atmósfera; y 10) densidad de población (D), se calculó con base en el último censo de población de INEGI (2010) con los comandos Spatial Analyst Tools>Density> Kernel Density de Arc Gis 10.
Las covariables utilizadas fueron escaladas a valores de 0 a 256 con el propósito de homogenizar la información, y someterlas al análisis de componentes principales (ACP). Con gráfica de densidad se detectó que de las diez capas que representaban a las covariables (o los 10 componentes principales), con solo seis componentes principales se pudo explicar el proceso de formación de cárcavas, representando el 96.7% de la información. Estos seis componentes se usaron para identificar, con base en las áreas de entrenamiento y el método de máxima verosimilitud, la probabilidad de presencia de cárcavas en cada pixel del área de estudio. De acuerdo con lo que proponen Rossiter (2002), la validación de resultados se realizó con 1.1 puntos de inspección por km2. Con el comando Data Management Tools>Feature class> Create Random Points, de Arc Gis 10, se establecieron 289 puntos al azar y en cada uno de ellos se establecieron tres áreas concéntricas (con radios de 50, 100 y 150 m) para observar en cada círculo la proporción (el área que abarca) de cobertura de cárcavas (Figura 3).
Resultados y Discusión
Relación entre covariables y la presencia de cárcavas
Con base en 502 áreas de entrenamiento empleadas en este trabajo, a continuación, se indica como ocurre la presencia de cárcavas según la magnitud de cada covariable ambiental analizada.
Pendiente (PEN). Las cárcavas del área de estudio se presentan en pendientes entre 1 y 77%, pero concentrándose con mayor frecuencia en el rango 2.5 a 6.5% (29.16%) y en el intervalo 7.5 a 22.5% (51.57%). Las pendientes más bajas se asocian a atrincheramientos en las inmediaciones de valles intermontanos y las más altas a cárcavas ubicadas en lomeríos compuestos de materiales residuales. Para esta covariable, Palacio (1990) encontró que en el centro de México las cárcavas se desarrollan en pendientes de 0 a 18% y más usualmente entre 2 y 8%. Aun cuando en pendientes superiores existe mayor riesgo de erosión remontante, es menos común su presencia ya que generalmente estos terrenos son someros y cubiertos con vegetación.
Curvatura horizontal (CH). Las cárcavas en la zona se localizan principalmente entre -4.75 y 5.65 (adim.) y se concentran con mayor frecuencia en los rangos -1.75 y -0.45 (32.7%) y entre 0.05 y 1.75 (64.5%). En general, el rango negativo corresponde a flujos convergentes asociados a topografías cercanas a valles donde se generan corrientes con mayor acumulación de caudal y potencial erosivo, por otro lado, las curvaturas positivas se vinculan a flujos divergentes, como ocurre en topografías cónicas (volcánicas) con pendientes más pronunciadas.
Curvatura vertical o de perfil (CV). Las cárcavas en el área de estudio se presentan en curvaturas verticales entre -3.4 a 3.45; principalmente se concentran en el rango -0.75 a -0.05 (27.15%) y 0.05 a1.05 (68.49%). El rango negativo indica que la velocidad de la corriente se incrementa conforme desciende, lo que genera flujos con mayor potencial para causar erosión en materiales residuales consolidados en los lechos de las cárcavas. Por otra parte, los valores positivos se relacionan con cárcavas localizadas en zonas donde la sedimentación es el proceso dominante, es decir, estas se localizan en las inmediaciones de los valles aluviales.
Longitud de la pendiente (LP). En el área de estudio las cárcavas se presentan en laderas con longitudes entre 0 y 1300 m, pero principalmente se concentran en longitudes entre 10 y 30 m (22.18%) y en longitudes de 35 y 110 m (33.52%). En general se observa que el rango más bajo corresponde a áreas de depósitos aluviales, mientras las longitudes mayores a terrenos de ladera baja.
Índice de posición topográf ica (IPT). Se encontró que la mayoría de las cárcavas se ubican en IPT>1 y IPT<1, es decir, en crestas (49.69%) y en valles intermontanos (42.22%), respectivamente. Esto indica que las cárcavas ocurren principalmente en laderas con presencia de materiales residuales de toba volcánica y en las inmediaciones de los valles intermontanos donde hay presencia de materiales aluviales y producción de escurrimientos.
Índice topográfico de humedad (ITH). se encontró que las cárcavas se distribuyen entre 0 y 14 m2 m-1 °; pero ocurren principalmente entre 0.0 y 2 m2 m-1 ° (42.43%), que se relacionan con acumulaciones medias de humedad, pendiente moderada y materiales residuales de toba volcánica; también estas ocurren entre 2 y 7 m2 m-1 ° (56.64%) que corresponden a áreas con las mayores acumulaciones de agua, pendientes bajas y materiales aluviales, características que se conjugan en las inmediaciones de los valles intermontanos.
Distancia a picos (DP). Las cárcavas del área de estudio se distribuyen a una distancia, a los picos topográficos más cercanos, de 0 y 1300 m, pero se concentran a distancias entre 100 y 400 m (46.91%) y entre 500 y 900 m (29.66%). Las cárcavas más cercanas a picos se registran en materiales volcánicos residuales, que circundan antiguas estructuras volcánicas, y las más alejadas a depósitos aluviales perimetrales a los valles intermontanos.
Distancia a fallas (DF). En la zona las cárcavas se ubican a una distancia entre 0 y 12 000 m de la falla Tamazulapam, pero la mayor frecuencia ocurre entre 0 y 4000 m (49.87%), seguido del rango 4000 y 10 000 m (45.87%). Así, las cárcavas más cercanas a la falla están asociadas a sitios con mayor actividad volcánica y las más alejadas a la formación de valles aluviales. En este sentido, Irigaray y Chacón (2002) mencionan que la distancia a fallas es una covariable que contribuye a la inestabilidad de las masas de suelo.
Densidad de vegetación (DV). Las cárcavas en el área de estudio presentan valores entre -0.54 y 0.29 concentrándose estas en el rango entre -0.16 y 0.2 (99.02%), es decir se localizan en suelos desnudos y vegetación baja (López-Pérez et al., 2015). Esto concuerda con los resultados de Alatorre y Beguería (2009), que mencionan que las cárcavas se ubican principalmente en zonas marginales donde es difícil que se establezca la vegetación en pendientes superiores al 15%.
Densidad de población (D). en el área analizada se encontró que las cárcavas se ubican en sitios con densidades entre 0 a 170 hab km-2, distribuidas principalmente en los rangos 5 a 30 hab km-2 (60.91%) y de 30 a 70 hab km-2 (31.40%). Las densidades bajas de población corresponden a cárcavas localizadas en sitios de difícil acceso y escaso potencial agrícola mientras que las cárcavas localizadas en lugares con densidades medias corresponden a entornos con mayor actividad agrícola, como ocurre en las inmediaciones de los valles aluviales.
En la zona de estudio se identificaron dos tipos de cárcavas, según su origen biofísico, las primeras localizadas en zonas cercanas a fallas geológicas con presencia de actividad volcánica. Sobre todo, este tipo de cárcavas se presentan al pie de las laderas que circundan antiguas estructuras volcánicas donde se presentan cauces divergentes, el agua alcanza velocidades críticas conforme desciende la estructura volcánica y ocurren cambios de uso de suelo; condiciones similares a las encontradas por Bean et al. (2017) en las laderas compuestas por materiales volcánicos blandos. En estas cárcavas la acumulación y velocidad del flujo son suficientes para causar erosión en los materiales residuales de toba volcánica localizados en lechos y paredes de las cárcavas; materiales similares a los observados por Bocco et al., 1991. (Dardis et al., 1988; Bocco et al., 1991; Oostwoud W. et al., 1999; Prat et al., 2003 ; Bean et al., 2017).
El segundo tipo de cárcavas se observó en zonas donde la sedimentación ha sido el proceso dominante, es decir, en las inmediaciones de los valles aluviales intermontanos; similar a lo observado por Dardis et al. (1988) y Oostwoud Wijdenes et al. (1999) en el suroeste de España y Sudáfrica, respectivamente, donde encontraron este tipo de erosión en las colinas cercanas a los valles. En este segundo tipo de cárcavas, aunque la pendiente del terreno es menor que las anteriores, los escurrimientos deben su poder erosivo a una mayor concentración de los mismos, por la convergencia de los cauces y mayores superficies de sus microcuencas, que produce el atrincheramiento de los materiales aluviales depositados con anterioridad.
Covariables que influyen en la formación de cárcavas
En la matriz de covarianzas, resultado del ACP para las diez covariables analizadas, se observó en la diagonal (Cuadro 1) que la distancia a la falla Tamazulapam tiene mayor variabilidad (55.8%), seguida de la distancia a picos (12.3%), índice topográfico de humedad (11.1%), densidad de vegetación (9.0%), pendiente (5.0%), índice de posición topográfica (2.8%), densidad de población (1.5%), longitud de la pendiente (1.3%), curvatura vertical (0.6%) y curvatura plana (0.5%). Esto implica que las covariables DF, DP, ITH, DV explican el 88.2% de la variabilidad de la información, es decir, son las que más aportarán a los resultados.
Cuadro 1: Matriz de covarianzas del análisis de componentes principales.
| Covariable | CH | CV | D | DF | DP | ITH | LP | COV | PEN | ITP |
|---|---|---|---|---|---|---|---|---|---|---|
| CH (Adim.) | 15.14 | 11.87 | -0.01 | 0.00 | -2.79 | -23.19 | -6.12 | 3.02 | -5.24 | 12.52 |
| CV (Adim.) | 11.87 | 17.68 | 0.01 | -0.38 | -0.99 | -26.81 | -8.00 | 2.28 | -1.48 | 8.96 |
| D (hab km -2) | -0.01 | 0.01 | 42.86 | 15.45 | -1.11 | 5.62 | -1.34 | 2.51 | -5.07 | -0.37 |
| DF (m) | 0.00 | -0.38 | 15.45 | 1590.44 | -10.34 | -17.74 | 1.70 | 44.64 | 27.83 | 0.60 |
| DP (m) | -2.79 | -0.99 | -1.11 | -10.34 | 349.72 | 32.37 | 14.53 | -22.28 | -13.56 | -40.53 |
| ITH (m 2 m -1 °) | -23.19 | -26.81 | 5.62 | -17.74 | 32.37 | 317.26 | 56.79 | -9.97 | -20.56 | -45.25 |
| LP (m) | -6.12 | -8.00 | -1.34 | 1.70 | 14.53 | 56.79 | 37.04 | -1.64 | 1.22 | -14.97 |
| COV | 3.02 | 2.28 | 2.51 | 44.64 | -22.28 | -9.97 | -1.64 | 256.60 | -0.10 | 9.51 |
| PEN (%) | -5.24 | -1.48 | -5.07 | 27.83 | -13.56 | -20.56 | 1.22 | -0.10 | 143.48 | -13.82 |
| IPT (Adim.) | 12.52 | 8.96 | -0.37 | 0.60 | -40.53 | -45.25 | -14.97 | 9.51 | -13.82 | 79.31 |
CH = curvatura plana u horizontal; CV = curvatura perfil o vertical; D = densidad de población; DF = distancia a la falla Tamazulapam; DP = distancia a picos; ITH = índice topográfico de humedad; LP = longitud de la pendiente; COV = cobertura vegetal; PEN = pendiente; IPT = índice de posición topográfica.
Al observar la matriz de correlaciones (Cuadro 2) destaca que la curvatura vertical y horizontal presentan alta correlación entre ellas (0.73) por lo que se puede prescindir de una de las dos sin alterar significativamente los resultados. Entre longitud de pendiente e índice topográfico de humedad se tiene una correlación intermedia (0.52), debido a que el índice de humedad considera la forma y longitud de la pendiente para modelar la distribución de la humedad en el terreno.
Cuadro 2: Matriz de correlaciones del análisis de componentes principales.
CH = curvatura plana u horizontal; CV = curvatura perfil o vertical; D = densidad de población; DF = distancia a la falla Tamazulapam-+}; DP = distancia a picos; ITH = índice topográfico de humedad; LP = longitud de la pendiente; COV = cobertura vegetal; PEN = pendiente; IPT = índice de posición topográfica.
Los vectores propios, generados por el ACP (Cuadro 3), mostró que en el componente uno (CP1) la covariable con mayor peso es la distancia a la falla de Tamazulapam, en los CP2 y CP3 distancia a picos y el índice topográfico de humedad, CP4 NDVI, CP5 pendiente, CP6 índice de posición topográfica, CP7 densidad de población, CP8 longitud de la pendiente, CP9 curvatura vertical, plana y longitud de la pendiente y en CP10 el mayor impacto lo tienen la curvatura plana y la curvatura vertical. El análisis estadístico reflejo que las variables ambientales detonantes de la formación de cárcavas son de tipo geológico, topográfico, de densidad de vegetación y antropogénico, de aquí que el área por naturaleza es frágil, al interaccionar con el incremento de población, la demanda de los recursos naturales se eleva lo que a su vez puede generar un mal manejo, traduciéndose como una región condenada a presentar una alta degradación física, es decir erosión con cárcavas.
Cuadro 3: Matriz de vectores propios del análisis de componentes principales.
CH = curvatura plana u horizontal; CV = curvatura perfil o vertical; D = densidad de población; DF = distancia a la falla Tamazulapam; DP = distancia a picos; ITH = índice topográfico de humedad; LP = longitud de la pendiente; COV = cobertura vegetal; PEN = pendiente; IPT = índice de posición topográfica.
Los autovalores para cada componente principal, generados por el análisis, proporcionan el porcentaje de variación que explica cada uno de ellos. Estos valores son útiles para seleccionar el número de componentes principales con la mayor concentración de datos originales. La manera más común de seleccionar el número de componentes idóneos es a través de una gráfica de densidad donde se muestre el porcentaje acumulado de autovalores por componente principal (Figura 4). En este caso, para determinar la probabilidad de presencia de cárcavas, se seleccionaron los primeros seis componentes principales que explican el 96.7% de la información contenida en las diez covariables analizadas.
Con base en la matriz de valores propios y la densidad de datos, que aporta cada componente principal, se observó que los aspectos ambientales que dominaron cada componente fueron: CP1 tipo geológico, CP2 geohidrológico, CP3 hidrogeológico, CP4 cubierta vegetal, CP5 topográfico y CP6 topohidrológico.
En la zona de estudio la presencia de cárcavas se asocia principalmente a la distancia a una falla geológica de importancia regional, seguido de la distancia a conos volcánicos, la concentración de escurrimientos conforme a la ubicación topográfica sobre el cauce, la pendiente del lecho del cauce y la posición topográfica relativa de los sitios. Esta última asociada a cárcavas en terreno de pie de monte, sobre suelos residuales de toba volcánica, y a atrincheramientos localizados en las inmediaciones de los valles intermontanos, conformados por antiguos depósitos aluviales.
Probabilidad de formación de cárcavas en el área de estudio
Al clasificar la imagen multicapa, con los seis componentes principales seleccionados, se encontró que el 90.17% del área de estudio no presenta cárcavas, mientras que el 9.83% de la superficie total (606.54 km2) tiene una probabilidad de existencia de cárcavas mayor 0.5%. Asimismo, el 4.23% de la superficie (261.11 km2) tiene una probabilidad superior de 50% de que exista esta forma de erosión (Figura 5, Cuadro 4).
Cuadro 4: Probabilidad de existencia de cárcavas
Probabilidad de existencia |
Área |
|
|---|---|---|
km2 |
% |
|
0 |
5560.92 |
90.17 |
0.5 |
5.68 |
0.09 |
1 |
5.07 |
0.08 |
2.5 |
18.75 |
0.30 |
5 |
39.16 |
0.64 |
10 |
90.01 |
1.46 |
25 |
186.76 |
3.03 |
50 |
146.95 |
2.38 |
75 |
75.7 |
1.23 |
90 |
27.56 |
0.45 |
95 |
6.4 |
0.10 |
97.5 |
2.6 |
0.04 |
99 |
1.23 |
0.02 |
99.5 |
0.35 |
0.01 |
100 |
0.32 |
0.01 |
Total general |
6167.46 |
100.00 |
Al contabilizar la presencia o ausencia de cárcavas en 289 puntos de verificación, con una probabilidad de existencia superior al 50%, se encontraron cárcavas en 278 puntos, lo que produce una confiabilidad global del mapa generado de 96.2%. Al considerar el 20% de afectación en los buffer se identificaron 251 sitios con cárcavas y 38 sin ellas, por lo que la certeza de predicción disminuyó a 86.9%.
La Figura 6 ilustra el comportamiento con relación a la superficie porcentual que abarca cada uno de los buffer con respecto al número de estos con afectación, la tendencia es clara, entre mayor sea el radio de los buffer la predicción de sitios con cárcavas disminuye, tomando en cuenta un porcentaje de afectación de 30‑40%, el buffer de 50 m predice una confiabilidad de 83.4%, el de 100 m de 81.7% y el de 150 m de 77.9%.
Al examinar el método empleado se dedujo que produce buenos resultados ya que conforme se reducía el radio de análisis, respecto a un punto con cárcavas observadas, se incrementaba la predicción de los resultados. En tanto que las covariables ambientales distancia a fallas, distancia a picos, índice topográfico de humedad, densidad de vegetación, pendiente e índice de posición topográfica, son una alternativa para identificar las áreas susceptibles a la formación de cárcavas de forma precisa y rápida.
Conclusiones
En el área de estudio las cárcavas se presentan en suelos residuales de toba al pie de monte de antiguos aparatos volcánicos y en depósitos aluviales que circundan los valles intermontanos.
A nivel regional la distancia a la falla de Tamazulapam y la presencia de aparatos volcánicos condicionaron la formación de cárcavas, seguido de factores geomorfológicos que controlan procesos de erosión y sedimentación, lo que se infirió a través del índice topográfico de humedad, pendiente del terreno, densidad de vegetación e índice de posición topográfica.
Las covariables ambientales, distancia a fallas, distancia a picos, índice topográfico de humedad, densidad de vegetación, pendiente e índice de posición topográfica, son una alternativa para identificar las áreas susceptibles a la formación de cárcavas en sitios conformados por suelos residuales de origen volcánico y en depósitos aluviales localizados en la periferia de los valles intermontanos.
El método de clasificación supervisada, usando ACP, permite identificar áreas con presencia de cárcavas con una precisión aceptable (96.2%) y es útil para determinar el potencial de un sitio para desarrollar este tipo de erosión.
