nueva página del texto (beta)
Español (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
La erosión de suelos en México, calculada con el método Ecuación Universal de Pérdida de Suelos (Universal Soil Loss Equation [USLE]) (Wischmeier y Smith, 1960), varía entre 492 y 1296 t/ha/año (Maass, 1998). La clasificación de la Organización de las Naciones Unidas para la Agricultura y la Alimentación (FAO, por sus siglas en inglés), el Programa de las Naciones Unidas para el Medio Ambiente (PNUMA) y la Organización de las Naciones Unidas para la Educación, la Ciencia y la Cultura (UNESCO, por sus siglas en inglés) para el USLE reúne las tasas superiores a 250 t/ha/año en un mismo rango como el valor máximo (erosión extrema), y de este modo se han registrado en los trabajos realizados en México. El mapa nacional de erosión potencial de México (Montes et al., 2011) establece más de 64% del territorio mexicano en este rango. Pablo (2014) calcula tasas de erosión con el método USLE en esta misma cuenca, con valores medios para cada atributo, propuestos por diferentes autores, para México, que establecen sobre la formación Yanhuitlán los valores máximos de erosión (superiores a 250 t/ha/año).
La erosión es un proceso sumamente complejo de analizar y ha sido estudiada (principalmente en suelos) con diferentes metodologías y modelos desde la década de 1950 hasta hoy, como la USLE (Wischmeier y Smith, 1960), el Chemical, Runoff, and Erosion from Agricultural Management Systems - CREAMS- (Knisel, 1980) o el Water erosion Prediction Project -WEPP- (Nearing et al., 1989) aplicados principalmente en terrenos agrícolas para analizar la erosión laminar y en surcos; no obstante, éstos no se pueden aplicar sobre superficies sin suelo. El caso de la erosión en cárcavas tampoco presenta modelos matemáticos de estimación y tiene que medirse mediante la utilización de métodos directos sobre el terreno. En México, los primeros trabajos de estudios de erosión acelerada en cárcavas con el uso de Sistemas de Información Geográfica (SIG), como herramienta de cartografía digital, se realizaron en los años noventa (Bocco, 1991, Palacio y Vázquez, 1990). Es también en estos años cuando se llevó a cabo el primer trabajo de fotogrametría digital para análisis de erosión en cárcavas en dos dimensiones con el uso de un globo de helio, como vehículo aéreo no tripulado (Unmaned Aerial Vehicle [UAV]), equipado con una cámara para evaluar la dinámica regresiva de cárcavas en varias fechas, con lo cual se obtuvieron resultados satisfactorios (Palacio, 1990).
La erosión hídrica de una formación litológica depende de factores exógenos (clima y cubierta vegetal) y factores endógenos (densidad aparente, permeabilidad-saturación, dureza, pendiente). El agente externo fundamental es la precipitación, la cual tiene una fuerza erosiva diferencial según su intensidad y duración (Palacio y Vázquez 1990); los factores internos proporcionan la erodabilidad de cada formación. Uno de los intereses de esta investigación es mostrar cómo la fotogrametría digital puede resolver el problema del cálculo de erosión en zonas carentes de suelo, donde la USLE no se puede aplicar. Esto es de suma importancia para los estudios de erosión en esta zona de cárcavas, ya que es la que aporta mayores tasas de sedimentos.
La fotogrametría digital ha avanzado mucho desde entonces, tanto en software como en hardware, en obtención, procesamiento y producción de productos derivados de la fotografía digital, como los Modelos Digitales de Superficie (MDS) de alta resolución, pero, sin duda, el uso de drones para este fin ha supuesto toda una revolución, no sólo porque hace posible la adquisición ilimitada de fotografía aérea en tiempo y forma a bajo costo, sino también el desarrollo de la tercera dimensión en los MDS, lo que permite cálculos de volumen en los análisis del terreno (Stöcker, Eltner y Karrasch, 2015; Gómez-Gutierrez et al., 2014; Frankl et al., 2015; Eltner, 2013). Estos MDS de alta resolución se producen a partir de nubes de puntos con información XYZ generada por los softwares fotogramétricos. Dichas nubes de puntos se obtienen a partir de estereoscopia con las fotos aéreas, a diferencia de las que se generan con Light Detection and Ranging o Laser Imaging Detection and Ranging (Lidar) a través de láser. Ambas tecnologías tienen sus limitaciones a la hora de trabajar en zonas con vegetación. Si bien el láser ofrece mejores resultados al atravesar el follaje, ninguno de los dos ofrece un resultado satisfactorio para nuestro objetivo por la falta de densidad en las nubes de puntos bajo el follaje, por ello este trabajo se realiza en zonas sin vegetación. La ventaja de la fotogrametría digital de alta resolución con respecto a Lidar en el resto de superficies es el menor costo económico en la obtención de datos.
El uso de drones para estos fines es algo relativamente nuevo. Las primeras imágenes aéreas de vehículos no tripulados con radiocontrol para fines científicos fueron tomadas en 1979 (Przybilla y Wester-Ebbinghaus, 1979), pero no será hasta el siglo xxi cuando se desarrollen los que hoy conocemos como drones y comiencen los primeros trabajos científicos. En 2004 se publicaron los tres primeros artículos relacionados con mapeo fotogramétrico con drones en el isprs Journal of Photogrammetry and Remote Sensing, y en 2008 empezó a ser tendencia con 21 artículos (Remondino et al., 2011). A partir de este año se amplió el uso de drones para múltiples fines científicos. La aplicación para estudios de geomorfometría y erosión fue creada a comienzos de la segunda década de este siglo (Eltner et al., 2016).
SFM (STRUCTURE FOR MOTION)
La fotogrametría clásica se basa en el principio de estereoscopía a través de imágenes aéreas con más de 60% de sobreposición, con el que se generan modelos digitales del terreno. No obstante, para poder establecer estudios de cambio superficial multitemporal, los métodos fotogramétricos tradicionales requerían muchos ajustes técnicos de orientación de las imágenes, tanto en la toma como en el procesamiento de datos, los cuales no llegan a ofrecer precisión y confianza en los resultados (Eltner, 2016; Stöcker et al , 2015). La técnica Structure for Motion (SfM) (Ullman, 1979) permite solventar estos ajustes de orientación de manera automatizada y confiable en los actuales softwares fotogramétricos. Con el desarrollo tecnológico de los drones, el geoposicionamiento de cada fotografía ya es automático, no sólo en coordenadas XYZ, sino también en (valores de rotación del eje del aparato que pueden afectar cada fotografía), lo que asegura la correcta geometría de los modelos por medio de los puntos coincidentes en las fotos. El SfM permite, por tanto, posicionar un modelo sobre otro con el conocimiento de que la geometría de ambos ya ha sido rectificada por el software fotogramétrico. El uso de dicha técnica es de vital importancia en esta investigación porque elimina el error geométrico en los modelos 3D.
OBJETIVOS
Los objetivos de este trabajo son analizar la manera en que la erosión modela diferentes litologías y calcular la tasa de pérdida en volumen (m3/area) y masa (t/area) en tres litologías diferentes en el valle de Yanhuitlán, Oaxaca. Para ello se levantó un modelo digital del terreno de alta resolución (10 cm) de una ladera con 200 m de longitud con similares condiciones de pendiente, orientación, clima y cubierta vegetal, sobre la cual se analizó la morfometría del terreno antes y después de las lluvias, con lo que se pudieron calcular, al mismo tiempo, los cambios en volumen y la masa perdida.
ÁREA DE ESTUDIO
Para poner a prueba los métodos planteados en esta investigación, se eligió un sector del valle de Yanhuitlán, situado en la Mixteca Alta, una región geográficamente diversa en el oeste de Oaxaca y áreas de los estados vecinos de Puebla y Guerrero, al sureste de México. Esta es una zona de contacto entre la Sierra Madre del Sur y la Sierra Madre Oriental, con una topografía de montañas medias y amplios valles con laderas fuertemente erosionadas (Palacio et al., 2016). No obstante, la erosión no sigue patrones similares en estas laderas, con diferencias en profundidad de disección, densidad de disección y patrón de la red de escurrimientos, lo que produce un modelado erosivo diferente en cada formación litológica. Debido a estos contrastes y a la intensa dinámica geomorfológica, este escenario es ideal para los objetivos de esta investigación.
El clima de la cuenca de Yanhuitlán es templado subhúmedo (Cw1), de acuerdo con la clasificación climática de Köppen, modificada por García (2004) para México, con una temperatura media anual de 15 °C y una precipitación media anual de 633 mm (según datos de la Estación Climatológica del Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias [INIFAP] en Santo Domingo Yanhuitlán), con lluvias concentradas en verano y un periodo de canícula marcado de julio a agosto. La vegetación es, en su mayor parte, secundaria, debido al sobrepastoreo y/o al cambio de uso del suelo por agricultura y ganadería. Predomina el chaparral y el pastizal inducido con algunos relictos de bosque de táscate en la parte baja e intermedia de la cuenca y bosque de encino-pino perturbado en las partes altas de la cuenca con algunos relictos primarios. No obstante, este intenso uso ha causado la pérdida total de suelo y vegetación en muchas zonas de la cuenca, donde hoy día se produce una intensa erosión de la roca madre, objeto de este estudio.
La geología del valle de Yanhuitlán comprende principalmente tres litologías de diferente naturaleza, cuyas edades están comprendidas entre el Eoceno y el Cuaternario (Figura 2). Por orden de antigüedad son las siguientes.
La formación Yanhuitán es una sucesión de finas capas rítmicas arcillo-limosas. Los sedimentos fueron depositados en ambientes lacustres durante el Eoceno, cuya edad estimada es de 40.5 ± 1.7 Ma (Martiny et al., 2000) con base en fechamientos de las intrusiones magmáticas de la formación Yanhuitlán. Es la formación litológica que se espera más deleznable, debido a su escasa consolidación de origen. El tipo de erosión de esta formación es en cárcavas. El test de dureza (USDA, 2012) practicado en campo reportó categoría 1 (0.60-1.25 MPa); puede romperse con el dedo.
La toba Llano de Lobos cubre la formación Yanhuitlán. Está conformada por secuencias de tobas riodacíticas y andesíticas con intercalaciones de calizas en una matriz arenosa. Su edad se estima entre 40.5 ± 1.7 Ma y 28.9 ± 0.6 Ma, por relaciones estratigráficas (Santamaría-Díaz et al., 2008). El tipo de erosión de esta formación es laminar, con una dureza categoría 3 (5-12.5 MPa); los fragmentos pueden ser rotos con la mano (USDA, 2012).
La andesita Yucudaac, de origen volcánico, sobreyace discordantemente en la toba Llano de Lobos y tiene una edad estimada inferior a 28.9 ± 0.6 Ma (Ferrusquía-Villafranca, 1976). Se trata de una serie de derrames lávicos andesíticos de color gris oscuro o negro, con un alto grado de intemperismo a través de sus fisuras. El tipo de erosión de esta formación es en surcos o regueros, con una dureza categoría 5 Hard Rock (50-100 MPa); se necesitan varios golpes para romper y puede ser rayada con una aguja (USDA, 2012).
MATERIALES Y MÉTODOS
Levantamiento fotogramétrico con dron
1) Con un drone Phantom 3 se realizó un vuelo fotogramétrico tipo POI (Point of Interest), que consiste en elegir un punto central sobre el cual volar, con una altura y radio determinados, para realizar fotografías aéreas oblicuas con una sobreposición suficiente para poder generar la nube de puntos a partir de estereoscopía.
En este caso se realizaron dos vuelos similares a 120 m de altura y 150 m de radio del punto central, con lo que se obtuvo un total de 68 fotos con más de 60% de sobreposición (Figura 3).
2) Una vez realizado el trabajo en campo, se procedió a la descarga e inclusión de las imágenes al software fotogramétrico Pix4D. Éste reconoce los datos técnicos de la cámara que soporta el dron, tanto la geolocalización de cada foto como las variaciones en producidas por el movimiento de la aeronave, que vienen incluidas en los metadatos de las fotos. Con este software se obtuvo la nube de puntos en formato LAS (Lidar) a partir de la estereoscopía digital obtenida de los diferentes ángulos de visión de un mismo punto en cada fotografía aérea.
3) Una vez obtenida la nube de puntos en formato, LAS se incluyó en el software Global Mapper 16, con el que se generaron los MDS en formato raster a resolución de 10 cm/pixel. Estos MDS se generan a través de la triangulación de estos puntos con el método de triangulación de redes irregulares (TIN, por sus siglas en inglés), donde se elige la opción minimum digital surface que toma el valor mínimo de los puntos en la nube.
4) Al MDS se le generó un análisis hídrico superficial en formato shape con la herramienta Watershed en Global Mapper 16, tanto de los escurrimientos como de su área de captación, los cuales fueron analizados por litología en parcelas de 50 x 50 m establecidas en zonas sin vegetación (Figura 4).
Fuente: elaboración propia.
Figura 4 Escurrimientos, áreas de captación y zonas de muestreo: andesita (1), toba (2) y arcillas (3).
5) Posteriormente se incluyeron las coberturas de escorrentía superficial en el software ArcGIS 10.2, con el que se analizaron diversos parámetros morfométricos: longitud de los escurrimientos o cauces, área de captación de los escurrimientos, disección vertical y densidad de disección, así como los órdenes de drenaje que se forman por litología. Los cálculos de estos parámetros se realizaron por parcelas de 50 x 50 m y por m2.
Dureza y densidad de las rocas
La dureza de las rocas fue determinada en campo con la metodología y terminología del United States Departament of Agriculture (USDA, 2012) en referencia al comportamiento por el golpeo de un martillo geológico. Los valores están referidos en megapascales (MPa).
Se determinó la densidad aparente de las rocas con objeto de cuantificar la masa real que se pierde por erosión con base en el volumen calculado. Se aplicaron dos métodos diferentes de laboratorio según la dureza de cada tipo litológico. En el caso de rocas blandas, se introdujo un cilindro metálico de 100 cm3 de volumen directamente en la superficie, luego se secó el material extraído en un horno a 100° durante 24 horas y se obtuvo su masa en seco; la densidad se determinó a partir de la masa seca y el volumen del cilindro. En el caso de la andesita, al no ser posible introducir el cilindro por su dureza, se calculó mediante la inmersión en agua de una muestra seca de roca previamente recubierta con barniz para evitar la saturación de los poros, observando el volumen desplazado por la roca al introducirla en una probeta milimetrada.
La literatura presenta valores de densidad de este tipo de rocas entre 2.07-2.49 g/cm3 para las andesitas, 1-1.3 g/cm3 para las tobas y 0.8-1.6 g/cm3 para las arcillas redepositadas (mudstones) (Paniukov, 1981; USDA, 2012).
Ecuación Universal de Pérdida de Suelo (USLE en inglés )
Los valores para la USLE fueron calculados con la metodología de Wischmeier y Smith (1960). Los valores de erodabilidad fueron estimados en campo mediante la descripción de perfiles de suelo y calculados en una supuesta ladera de 50 m (igual a las de las parcelas estudiadas en cárcavas). Se realizaron cálculos con pendiente de 10° y de 20° en cultivo de maíz con arado en primavera y surcos transversales a la pendiente (como se practica en la zona) y un valor R de erosividad de la lluvia basado en un total anual de 633 mm y calculado con base en el Índice Modificado de Fournier (IFM) (Velasco y Cortés, 2009).
Pérdida de material por erosión
El cálculo de erosión por litología se realizó con la comparación de los MDS de dos fechas (antes y después del periodo de lluvias) para determinar el volumen de material que se pierde por cada litología, las cuales fueron el 16 de mayo y el 26 de septiembre de 2016. La precipitación medida en estas fechas fue de 523 mm (Estación Climatológica del INIFAP en Santo Domingo Yanhuitlán).
Para poder comparar ambos MDS, se geoposicionaron exactamente sobre puntos de control terrestres tomados con GPS de alta precisión en zonas donde no pueden ser modificados por la erosión en estos lapsos de tiempo, como estructuras viales, construcciones, etcétera (Figura 5). Se establecieron tres Puntos de Control Terrestre (GCP, por sus siglas en inglés) y se calculó el error diferencial con la estación geodésica más cercana, en este caso, la de Nochixtlán, a 29 km, para geolocalizarlos también con precisión en el espacio (geoposición absoluta).
La comparación de MDS de las dos fechas se realizó con álgebra de mapas en un Sistema de Información Geográfica (ArcGIS 10.2). La ecuación 1 consistió en restar el modelo del día 25 de septiembre al del 16 de mayo, de manera que se generó un nuevo MDS que mostró cambios en la altura (Z) con valores negativos en las zonas donde se produjo la erosión.
Ecuación 1: MDS CAMBIO = MDS1 (25sep) - MDS2 (16mayo)
Al conocer el tamaño de pixel (10 x 10 cm = 100 cm2) podemos establecer el cambio en volumen por pixel al sumarle la diferencia en Z (resultado de la ecuación 1). Para ello se multiplicó por100 el MDS CAMBIO (ecuación 2), y se obtuvo el volumen de cambio en cm3 (antes de realizar esta operación, se pasó la información de MDS CAMBIO de metros a centímetros).
Ecuación 2: MDS CAMBIO x 100 = MDS VOLUMEN
De esta manera, al sumar el valor de todos los pixeles con valor negativo en el MDS VOLUMEN, obtenemos la pérdida de suelo en volumen de cada unidad litológica en un área de 2500 m2. Con este dato pudimos calcular tasas de erosión en m3/m2 y m3/ha. Las parcelas tienen una pendiente de 20%, orientación SE y una precipitación de 532 mm. Una vez calculada la pérdida en volumen, se puede determinar la pérdida de erosión en masa de cada formación superficial mediante la fórmula expresada en la ecuación 3.
Ecuación 3: Masa = Densidad aparente * Volumen
Cálculo de error entre modelos
Para calcular el error en las mediciones, se establecieron seis puntos de revisión en zonas de “no cambio”: tres puntos usados para el corregistro de los modelos (gcp1, gcp2 y gcp3) y tres independientes en las zonas circundantes al modelo (chp1, chp2 y chp3) para calcular el error en la geometría y posición de los modelos en Z, ya que X y Y quedan resueltas en el software con los puntos GCP a través del geoposicionamiento de cada fotografía. La distribución de los puntos dependió de la ubicación de estructuras fijas (construcciones de mampostería), y se persiguió siempre la mejor relación geométrica por triangulación (Figura 5). El error, por tanto, se calculó con la diferencia en el valor Z de estos seis puntos, de efectivo carácter fijo, entre el MDS1 y el MDS2.
El error en Z entre modelos es debido a las diferencias geométricas en la triangulación de las nubes de puntos, ya que al extraer la información altimétrica de los modelos con los puntos GPS, se obtiene el valor que genera la triangulación en ese punto, que no coincide exactamente sobre el punto GPS. Los errores de estos puntos oscilan entre 1 y 5 cm (Tabla 1), por tanto, se tomará el valor máximo de error altimétrico (5 cm) para calcular el error máximo en volumen. Para disminuir el error es necesario realizar un vuelo más próximo al área de interés que elimine estas variaciones espaciales en la nube de puntos de manera que la triangulación sea más suavizada (Tabla 1).
Tabla 1 Error entre MDS1 y MDS2.
| PUNTO | LONGITUD | LATITUD | ALTURA MDS1 | ALTURA MDS2 | ERROR (mds2-mds1) |
|---|---|---|---|---|---|
| ch1 | 17.580 | -97.391 | 2471.287 | 2471.276 | -0.011 |
| ch2 | 17.580 | -97.392 | 2486.790 | 2486.758 | -0.032 |
| ch3 | 17.581 | -97.391 | 2413.255 | 2413.244 | -0.011 |
| GCP1 | 17.580 | -97.391 | 2454.489 | 2454.434 | -0.055 |
| GCP2 | 17.580 | -97.390 | 2441.897 | 2441.881 | -0.016 |
| GCP3 | 17.581 | -97.390 | 2434.385 | 2434.339 | -0.046 |
RESULTADOS
Morfometría
La longitud de los escurrimientos y su área de captación indica la forma en que se modela una superficie. La disección vertical muestra qué tan profundo se erosiona, y la densidad de disección, cuántos escurrimientos se forman en un área dada. Algunos de estos parámetros tienen relación entre sí, por ejemplo: a mayor área de captación, menor tamaño de escurrimientos. A mayor longitud de escurrimientos, menores órdenes de drenaje.
Los resultados obtenidos de las mediciones de los parámetros morfométricos de densidad de disección, disección vertical, longitud de escurrimiento y su área de captación, muestran resultados diferenciales en cada litología por los siguientes atributos: dureza, densidad aparente y consolidación (Tabla 2).
Tabla 2 Parámetros morfométricos por litología.
| Parámetros morfométricos (50 m2) | Andesita | Toba | Formación Yanhuitlán |
|---|---|---|---|
| Áreas de captación | |||
| N° microcuencas de drenaje | 2065 | 1813 | 2254 |
| Área máxima de drenaje (m2) | 8.06 | 9.36 | 6 |
| Área media de drenaje (m2) | 1.10 | 1.24 | 1.01 |
| Escurrimientos | |||
| N° escurrimientos | 2266 | 1939 | 2499 |
| Longitud máxima de escurrimiento | 8.46 | 12.15 | 6.90 |
| Longitud media de escurrimiento | 1.10 | 1.44 | 1.08 |
| Disección vertical (m/m2) | |||
| Máxima | 1.96 | 2.29 | 3.04 |
| Media | 0.33 | 0.44 | 0.86 |
| Densidad de disección | |||
| Escurrimientos (50 m2) | 2413 | 2091 | 2652 |
| Escurrimientos/m2 | 0.9652 | 0.8364 | 1.0608 |
| Longitud máxima/m2 | 3.90 | 3.60 | 5.01 |
| Longitud media | 1.10 | 1.23 | 1.25 |
Área de captación de los escurrimientos
El área de captación (microcuenca) de cada cauce es mayor en la toba, con una media de 1.24 m2, seguida de la andesita, con 1.10 m2, y la formación Yanhuitlán, con 1.01 m2, este último debido a la baja consistencia del material que se desprende lateralmente por las paredes de las cárcavas, lo cual favorece más la disección vertical. Es por ello que en la formación Yanhuitlán es donde existen más cuencas de captación (2254 microcuencas/2500 m2), debido a su menor tamaño, frente a 2065 de la toba y 1813 de la andesita.
Longitud y patrones de los escurrimientos
La longitud de los escurrimientos indica la capacidad de infiltración de la formación litológica: cuanto más elongadas, menor infiltración. La toba tiene una expresión en el paisaje más elongada, con una longitud media de 1.44 m; la andesita y la formación Yanhuitlán tienen valores similares de 1.10 m y 1.08 m respectivamente, lo cual coincide con los datos de los parámetros areales. Los patrones de los escurrimientos reflejan una mayor concentración en las arcillas (dentrítico), erosión laminar en la toba (patrones paralelos) y patrón de erosión de trellis en la andesita, con escurrimientos pequeños y enrejados (Figura 6).
Disección vertical
La disección vertical de las formaciones superficiales en esta ladera responde a la dureza del material geológico y su estructura, cuanto más deleznable es el material, más incisión a través de sus canales de escorrentía. En este caso, la formación Yanhuitlán (arcillas poco consolidadas) es la más deleznable, tiene una profundidad media de disección de 0.86 m y llegan a producirse profundidades de hasta 3 m. La toba Llano de Lobos (cenizas litificadas) responde con una profundidad de escorrentía media de 0.44 m, mientras que la andesita yucudaac (ígnea extrusiva), la más dura, responde con una disección media de 0.33 m (Figura 7).
Densidad de disección y órdenes de drenaje
La densidad de disección se refiere a la longitud de cauces por unidad de área. Los resultados del parámetro a esta resolución no reflejan con tanta verosimilitud las diferencias que se aprecian visualmente en cada litología, dado que el detalle detecta cualquier diferencia centrimétrica en la superficie, por lo que se obtuvieron todos los escurrimientos que se pueden generar, esto homogeneiza los resultados. Para conseguir mayores diferencias en los valores y reflejarlas, es necesario disminuir la resolución.
No obstante, los datos siguen la lógica que se espera con las arcillas de la formación Yanhuitlán que presenten una mayor densidad de disección tanto en número de escurrimientos/m2 (1.06) como por la suma de la longitud/m2 (1.25 m/m2), seguida de la andesita (0.96 m/m2), por su mayor susceptibilidad de meteorización, y por la toba, con 0.83 m/m2.
Los órdenes de drenaje aportan resultados más interesantes que refuerzan las diferencias en los tipos de escurrimientos. El porcentaje en tipo de orden de escurrimiento por litología confiere mayores escurrimientos de primer orden en la toba, más escurrimientos de segundo orden en la andesita y más de tercer orden en la formación Yanhuitlán. Este dato refuerza la información del tipo de modelado que sufre cada litología, donde una mayor cantidad de drenaje de primer orden representa un tipo de erosión más laminar en la toba y mayor concentración de escurrimientos de tercer orden en la formación Yanhuitlán, donde se encuentran las cárcavas y la red de drenaje está más ramificada (Tabla 3).
Pérdida de suelo por erosión en volumen y masa
En 2016, la Estación Climatológica de Yanhuitlán registró un total de precipitación de 696 mm, que fue superior al promedio histórico anual de 633 mm (Uniatmos, 2013). La precipitación registrada en la temporada de lluvias (mayo-septiembre) fue de 532 mm (84% del total registrado en el año). Las tasas de erosión efectivamente medidas fueron las correspondientes a la precipitación entre el 16 de mayo y el 25 de septiembre de 2016. En la Figura 8 podemos apreciar visualmente el cambio en valores de Z dentro de cada parcela, con gradación de tonos rojizos donde se han perdido sedimentos y gradación de tonos verdes donde se han ganado; en amarillo aparece la zona donde no se han registrado cambios de erosión ni sedimentación.
Fuente: elaboración propia.
Figura 8 Dinámica de sedimentos y tasa de erosión por pixel: andesita (1), toba (2) y arcillas (3).
Los pixeles con erosión coinciden con los de mayor disección vertical de los parámetros morfométricos. La mayoría de pixeles con valores negativos por erosión se encuentran entre 0 y 3 cm de altura (Figura 8), la suma de estos pixeles negativos multiplicada por el área del pixel arrojó el volumen de material erosionado. Éste presenta diferencias importantes por litología.
Los resultados se presentan en rangos de valores, que representan la medición en la diferencia volumétrica de los modelos (dato mayor); el dato menor es el resultado de restarle el error máximo (5 cm). Al resultado de cambio entre modelos (MDScambio) se le restará en cada pixel el error máximo, lo cual implica generar en este intervalo valores que no representen error alguno. Por parcela de análisis en el periodo de lluvias (16 de mayo al 25 septiembre de 2016), los valores de pérdida estuvieron entre 95-123 m3/2500 m2 en la formación Yanhuitlán, 25-32 m3/2500 m2 en la toba y 9-12 m3/2500 m2 en la andesita. Esto, traducido a tasa de erosión por m2, equivale a 0.038-0.0492 m3/m2, 0.010-0.013 m3/m2 y 0.004-0.005 m3/m2. Si lo calculamos en tasa por hectárea, representa una pérdida entre 380-492, 100-128 y 37.48-48 m3/ha, respectivamente, con una precipitación de 532 mm. Si se interpolan estos datos con la tasa anual de precipitación de 696 mm, en 2016 quedan en 500-641 m3/ha en la formación Yanhuitlán, 130-167 m3/ha en la toba Llano de Lobos y 52-67 m3/ha en la andesita (Tabla 4).
Tabla 4 Comparativa de erosión en volumen y masa por litología.
| PARÁMETROS MORFOMÉTRICOS (50X50m) 10 cm resolución | Andesita | Toba | F. Yanhuitlán YYanhuitlán |
|---|---|---|---|
| Áreas (microcuencas) m 2 | |||
| N° microcuencas de drenaje | 2065 | 1813 | 2254 |
| Área máxima de drenaje | 8.06 | 9.36 | 6 |
| Área media de drenaje | 1.10 | 1.24 | 1.01 |
| Lineal (m) | |||
| N° escurrimientos | 2266 | 1939 | 2499 |
| Longitud máxima de escurrimiento | 8.46 | 12.15 | 6.90 |
| Longitud media de escurrimiento | 1.10 | 1.44 | 1.08 |
| Disección vertical (m/m 2 ) | |||
| Máxima | 1.96 | 2.29 | 3.04 |
| Media | 0.33 | 0.44 | 0.86 |
| Densidad de disección (escurrimientos/m 2 ) | |||
| líneas drenaje | 2413 | 2091 | 2652 |
| líneas de drenaje/m2 | 0,9652 | 0,8364 | 1,0608 |
| longitud máxima/m2 | 3.90 | 3.60 | 5.01 |
| longitud media | 1.10 | 1.23 | 1.25 |
Para calcular la pérdida en valores de masa, toneladas/hectárea (t/ha), se aplicaron los valores de densidad aparente superficial de cada litología al resultado obtenido del volumen, según la fórmula de cálculo de masa. Las densidades que obtuvimos en campo y laboratorio fueron las siguientes:
AY= 2.45 g/cm3
TLL= 1.10 g/cm3
YAN=1.40 g/cm3
Por tanto, los resultados de pérdida de suelo en masa en el periodo medido resultan en 532-689 t/ha en la formación Yanhuitlán, 110-141 t/ha en la toba Llano de Lobos y 92-118 t/ha/año en la andesita yucudaac. Las tasas anuales de producción de sedimento son de 701-898 t/ha/año, 143-183 t/ha/año y 128-164 t/ha/año, respectivamente. La Figura 9 muestra gráficamente las diferencias marcadas por litología.
Discusión: cálculo de erosión en suelo y zonas carentes de suelo
Los datos de erosión en este trabajo fueron calculados en zoas donde la cubierta edáfica se ha perdido debido a la intensa erosión, por lo cual los valores calculados sólo se podrían extrapolar en zonas similares. Los modelos de predicción de la erosión, como la USLE, sólo pueden ser aplicados en laderas sujetas a erosión pluvial y laminar. El método aquí empleado representa un complemento ideal para ajustar los cálculos de producción de sedimentos en cuencas en las zonas carentes de suelo, las cuales representan altos porcentajes.
El caso más relevante en el establecimiento de esta diferenciación de tasas de erosión por zonas es la formación Yanhuitlán, donde se generan las cárcavas cuando el suelo desaparece. Si observamos los valores de erosión de la formación Yanhuitlán en zonas cultivadas, con datos tomados en campo y los resultados de pérdida de suelo sobre cárcavas calculados con fotogrametría digital, los valores medidos en cárcavas superan sustancialmente (485 a 637 t/ha/año) a los obtenidos en la usle para campos cultivados (Tabla 5). Esto implica que las tasas de erosión calculadas en esta litología en zonas de cárcavas casi pueden llegar a duplicarse.
USLE= R x K x LS x C x P
Donde:
R: |
lluvia: 633 mm = 226.7 |
K: |
erodabilidad: 0.52 |
LS: |
pendiente y longitud de pendiente 15.27 |
C: |
tipo de cultivo/manejo: maíz (0.40) x arado en primavera (0.90) = 0.40 x 0.90 = 0.36 |
P: |
Prácticas de conservación de suelo: en contra a la pendiente (0.75). |
Tabla 5 Tasas de erosión en la formación Yanhuitlán calculadas por medio de la USLE (Wischmeier, 1960) para campos cultivados con maíz y por medio de fotogrametría digital con drones para zonas con cárcavas. Elaboración Propia.
| EROSION EN FORMACIÓN YANHUITLÁN | Tn/Ha |
|---|---|
| Cultivado con maiz y arado horizontal (USLE) Pend=10° | 138.2 |
| Cultivado con maiz y arado horizontal (USLE) Pend=20° | 485.4 |
| Cárcavas (DRON) Pend=20° | 637.3-819.3 |
CONCLUSIONES
Las imágenes aéreas de alta resolución obtenidas desde drones o UAV y la fotogrametría digital son herramientas útiles para el análisis morfométrico del relieve a gran detalle y su uso se puede extender a muchos campos de la investigación, entre ellos, el modelado del paisaje a través de la erosión.
Los efectos de la diferente erodabilidad en las tres litologías de la cuenca estudiada producen modelados erosivos característicos en cada una, expresados en diferencias de densidad de disección, profundidad de disección, microcuencas de captación y longitudes de los escurrimientos, así como en las tasas de erosión. Esto puede apreciarse a simple vista en cárcavas en la formación Yanhuitlán, por su alto contenido de arcilla con alta desagregación, pequeños surcos producidos por el alto grado de intemperismo en la andesita yucudaac o erosión laminar en la toba Llano de Lobos, debido a su poca infiltración.
El detalle del levantamiento fotogramétrico con drones permite analizar el diferente modelado de la erosión en cada formación litológica.
Las arcillas de la formación Yanhuitlán son el material más deleznable (erosionable) por su escasa consolidación, con mayores índices de densidad de disección y disección vertical. Por ello también tiene el mayor número de cuencas de captación con mínimo tamaño.
La toba es la que mayor superficie de captación genera para los escurrimientos.
La andesita, al ser la roca más dura, tiene la menor disección vertical.
Con levantamientos multitemporales por medio de fotogrametría digital con drones es posible cuantificar la erosión real en volumen y masa de una unidad litológica y estimar tasas de erosión y producción de sedimentos en zonas carentes de suelo y/o sujetas a erosión en cárcavas.
Las arcillas de la formación Yanhuitlán se caracterizan como las de mayor pérdida de volumen y masa, con una gran diferencia respecto de las demás litologías. En el periodo analizado (mayo a septiembre de 2016) se determinó una tasa de erosión de 532-689 t/ha. Si se extrapola este dato a todo el año con base en el total de lluvia anual, entonces resulta una tasa de 701-898 t/ha/año.
La toba Llano de Lobos tuvo una tasa de pérdida entre 110-141 t/ha en el periodo analizado y una tasa anual estimada de 143-183 t/ha/año.
La andesita es la formación que menos pérdida de material tiene, con una tasa de pérdida entre 98-112 t/ha en el periodo analizado y una tasa anual estimada de 128-164 t/ha/año.
Este método resuelve el problema del cálculo de erosión en zonas carentes de suelo, donde los métodos de estimación no pueden determinarlo, sobre todo en zonas de cárcavas, que son las que mayores tasas de erosión generan. Por tanto, este método contribuye a estudios de erosión en cuencas más integrales, precisos y realistas, e incluye el análisis de la evolución de laderas, válido para estudios de análisis de riesgo.
De igual forma, se precisa mejorar la metodología de vuelo para obtener mejores resultados en la nube de puntos que permita minimizar el error de los MDS y así obtener tasas más precisas y no rangos.
