La estimación de las pérdidas de suelo por erosión es un elemento de suma utilidad para la planificación y toma de decisiones a diferentes niveles. A nivel predial permite comparar objetivamente diferentes alternativas de uso y manejo y seleccionar la que, ofreciendo el nivel de conservación deseado, cumpla con los objetivos de producción y sea más simple de llevar a la práctica. Para la elaboración de la presente guía se trabajó con los datos proporcionados por nuestro amigo Geomatico Josué Díaz, así como en el videoturial, el cual puede ser consultado en el siguiente enlace: https://www.youtube.com/watch?v=gFKHDa0dHfo
T.I. Hugo Jiménez Reyes
[email protected]
Introducción El modelo de erosión USLE (Universal Soil Equation) es uno de los más extendidos en el mundo, permitiendo calcular la erosión potencial de un territorio y su caracterización en estados erosivos de diferentes niveles. Las pérdidas de suelo según este modelo, se calculan mediante la siguiente expresión:
A=R⋅K⋅L⋅S⋅C⋅P Siendo:
A las pérdidas de suelo por unidad de superficie (t/ha・ano) El factor R se define como el índice de erosividad de las precipitaciones. El factor K se define como el índice de erosionabilidad del suelo El factor L·S mide la erosionabilidad debida a la forma del terreno, combina de forma conjunta el factor L, longitud de ladera, con el factor S, pendiente del terreno El factor C se define como el factor de cobertura vegetal y muestra la relación de las pérdidas de suelo que se producen bajo unas condiciones determinadas de tipo de vegetación, cobertura vegetal (%) y técnicas de cultivo El factor P se denomina factor de prácticas de conservación, es considerado independiente del factor C, en aquellos casos en que los cultivos se realicen según cultivo a nivel, cultivo por fajas y cultivo en terrazas. Este factor también varía con la pendiente
Para conocer el significado de los valores de perdida de suelo obtenidos se utiliza la Clasificación Provisional para la Evaluación de la Degradación de los Suelos, confeccionada por la F.A.O., P.N.U.M.A. y U.N.E.S.C.O. y que se muestra a continuación: Pérdidas de suelo A (t/ha*año) <10 10-50 50-200 >200
Grado de erosión hídrica Ninguna o ligera Moderada Alta Muy Alta
Otro factor para los suelos es el llamado “Valor T” que significa “Pérdida tolerable de suelos”. No es directamente utilizado en la RUSLE, pero se usa en conjunto con ella para planificación de la conservación. La pérdida tolerable de suelos (T) es la máxima cantidad de pérdida de suelo en Ton/ha/año, que se puede tolerar y todavía permitir una alta productividad de cultivos de manera económica y sostenible. La producción de sedimentos no debe de confundirse con la erosión y los términos no son, por lo tanto, intercambiables. La producción de sedimentos es la cantidad de suelos erodado que se transporta a un punto en la cuenca que está lejos del origen del desprendimiento de partículas de suelo. En una cuenca, la producción de sedimento incluye la erosión de las pendientes, canales y remoción en masa menos el sedimento que se deposita después de que se erosiona, pero antes de que alcance el punto de interés (figura 1). La USLE y la RUSLE no estiman la producción de sedimento.
Figura 1. Perfil esquemático de la pendiente para aplicaciones de la RUSLE en erosión en surcos y entre surcos. El factor R R es el factor de erosividad de las lluvias. Es la suma de los valores promedio anual (EI) en un año normal de lluvia. El índice de erosión es una medida de la fuerza erosiva de la lluvia para una lluvia determinada. Cuando los otros factores son constantes, las pérdidas debidas a la lluvia son directamente proporcionales al producto de la energía cinética total de la tormenta (E) por la máxima intensidad de lluvia con 30 minutos de duración (I). Las tormentas de menos de 0.5” (12.5 mm) no se incluyen en los cálculos de erosividad ya que estas tormentas generalmente agregan muy poco al valor total de R. El factor R representa la EI de la tormenta promedio y por lo tanto entre más largo sea el registro con el que se determina mejor (se considera como mínimo un registro de 8 años). El factor R es una indicación de las características más importantes de la erosividad: la cantidad de lluvia y el pico de intensidad sostenido en un período determinado (30 minutos). Es de hacer notar que para hacer el análisis anual, solamente se escogerán aquellas tormentas que estén separadas más de 6 horas entre una y otra, para considerarlas independientes. Los índices de erosividad varían con la localidad y sus dimensionales son en MJ.mm.ha-1.hr-1.año-1. El factor K K es el factor de erodabilidad del suelo que representa tanto la susceptibilidad del suelo a la erosión como la tasa de escorrentía, medida bajo las condiciones de una parcela estándar unitaria. Los suelos con alto contenido de arcilla tienen bajos valores de K, de cerca de 0.05 a 0.15, por su resistencia al desprendimiento de las partículas de suelo. Suelos con texturas gruesas, como los suelos arenosos, tienen bajos valores de cerca de 0.05 a 0.2, ya que la baja escorrentía, aunque estos suelos son fácilmente erodables. Suelos con texturas medias como los francos, tienen un valor moderado de K, de cerca de 0.25 a 0.4 ya que son moderadamente susceptibles al desprendimiento y producen una escorrentía moderada. Suelos que contienen un alto contenido de limo son los más erodables. Son fácilmente desprendibles, tienden a formar una costra y producir altas tasas de escorrentía. Los valores de K de estos suelos tienden a ser mayores que 0.4.
La materia orgánica reduce la erodabilidad porque reduce la susceptibilidad del suelo al desprendimiento e incrementa la infiltración, lo cual reduce la escorrentía y por lo tanto la erosión. La adición o acumulación de materia orgánica a través del manejo como la incorporación de abonos verdes o estiércol se representa en el factor C, más que en el factor K. La extrapolación del factor K en el nomógrafo más allá de un contenido de 4% no se recomienda o se permite en la RUSLE. En la RUSLE, el factor K considera que el suelo como un todo y el factor Kf considera solamente la fracción fina o sea el material de menos de 2 mm de diámetro. Para la mayor parte de los suelos Kf = K. La estructura del suelo afecta, tanto la susceptibilidad al desprendimiento como la infiltración. permeabilidad del perfil del suelo afecta el factor K ya que afecta la escorrentía.
La
Aunque el factor K se seleccionó para representar un suelo en su condición natural, el manejo pasado o el mal uso de un suelo por cultivos intensivos pueden incrementar la erodabilidad del suelo. El factor K puede necesitar ser incrementado si el subsuelo está expuesto o donde la materia orgánica ha sido agotada, la estructura del suelo ha sido destruida o la compactación del suelo ha reducido la permeabilidad. Un especialista en suelos puede asistirlo para hacer esta interpretación. El factor de erodabilidad del suelo se calcula con la ecuación del nomograma de Wischmeier citado por Mannaerts (1999)
Donde: K = Factor de erodabilidad del suelo (t*acre*hr/100*acre.pie*tf*pulg) OM = Materia orgánica [%] S = Código de la estructura del suelo P = Código de permeabilidad M = Producto de las fracciones del tamaño de las partículas primarias ó (% limo + % arena muy fina)*(100 - % arcilla) A continuación se presenta la tabla con los datos de agua del suelo para las clases principales de textura de suelo: (Mannaerts, 1999) Textura clase
Permeabilidad
Conductividad saturada
SCS Grupo Hidrológico
código
(mm/hr)
de suelo
6
<1
D
5
1-2
C-D
4
2-5
C
3
5-10
B
2
10-60
A
1
>60
A
Arcilla, franco arcilloso Arcillo arenoso, franco arcilloso limoso Franco arcillo arenoso, franco arcilloso Franco limoso, franco Areno francoso, franco arenoso Arena
Códigos de permeabilidad y estructura del suelo en función de su textura
En la figura siguiente se muestra el nomograma para calcular el factor K. Debemos ser cuidadosos de verificar las unidades que arroja el nomograma pues los mismos deben multiplicarse por 0.1317 para que queden en unidades del sistema internacional (SI). Las dimensionales de K en el SI son: T*ha*h*ha-1*MJ-1*mm-1
Nomograma para el cálculo del factor K. (Unidades del Sistema Inglés, US).
Nomograma para la determinación de la textura del suelo.
Si no se cuenta con información de suelos, se puede calcular el valor de K basados en el tipo de roca formadora del suelo (litofacies) de acuerdo a la siguiente tabla:
Recuerde siempre multiplicar por el factor de conversión a unidades SI.
El factor LS L es el factor de longitud de la pendiente que representa el efecto de la longitud de la pendiente en la erosión (adimensional). Es la razón entre la pérdida de suelo en la longitud de pendiente en campo con respecto a la longitud de una parcela de 22.1 m en el mismo tipo de suelo y gradiente. La longitud de la pendiente es la distancia desde el origen de la escorrentía superficial a lo largo de su vía de flujo a la localidad donde se concentra el flujo o donde se deposita el suelo desprendido. Afortunadamente, los valores de pérdida de suelo no son especialmente sensitivos a la longitud de la pendiente y las diferencias en longitud de + o – 10% no son importantes en la mayoría de pendientes, especialmente en paisajes planos. Las longitudes de pendiente se determinan mejor visitando el sitio, caminando las vías de flujo y haciendo medidas directamente en el campo. Obtenga L a través de la medición perpendicular a la curva que origina el flujo superficial hasta donde se inicia la deposición o donde la escorrentía entra a un canal bien definido. Los mapas de curvas a nivel mayores a 0.50 m deben ser utilizados de manera cautelosa para determinar longitudes de pendiente. Los valores de longitud de pendiente son a menudo muy grandes cuando se utilizan mapas de curvas a nivel. Las longitudes de pendiente, generalmente no exceden los 120 m. Longitudes de pendiente mayores a 300 m no se deben utilizar en la RUSLE debido a la confiabilidad de la ecuación cuando se aplican estas longitudes es cuestionable y el flujo comienza a concentrarse en la mayoría de los paisajes antes de que ocurran estas longitudes de pendiente. S es el grado de inclinación de la pendiente. Representa los efectos de la inclinación en la erosión. La pérdida de suelo se incrementa más rápidamente con pendientes muy escarpadas en comparación con el incremento en longitud de la pendiente. Se define como el cociente de la pérdida de suelo en el campo entre la que ocurre en una parcela de 9% de inclinación bajo idénticas condiciones. La relación entre la pérdida del suelo y la gradientes se ve influenciada por la densidad de la cobertura vegetal y el tamaño de las partículas del suelo. Generalmente, el factor L y S se consideran juntos. Los factores LS juntos son iguales a 1 bajo condiciones de la parcela unitaria de 22.1 m y 9% de inclinación. Los valores de L y S son relativos y representan que tan erodable es una longitud e inclinación particulares en relación a los 22.1 m de largo y 9% de inclinación de la parcela unitaria. Por lo tanto, algunos valores de LS son menores que 1 y otros mayores que 1. Los cultivos en hilera o en contorno no afectan los valores de LS. En la figura siguiente se ilustran los diferentes tipos típicos de longitud de pendiente (Dissmeyer y Foster, 1980). La pendiente A: Si el bosque sin disturbar en la parte de arriba no produce escorrentía superficial, la pendiente superior comienza en el borde del suelo del bosque sin disturbar y se extiende aguas abajo hasta el surco o el canal, si la escorrentía es concentrada por el canal. La pendiente B: Punto de origen de escorrentía a un canal, si la escorrentía es concentrada por el canal. Pendiente C: Desde el canal a un punto de concentración de flujo. Pendiente D: Punto de origen de la escorrentía a un camino que concentra la escorrentía. Pendiente E: Desde el camino a una planicie de inundación donde ocurre la deposición. Pendiente F: En la cima de una colina, desde el punto de origen de la escorrentía a la planicie de inundación donde ocurre la deposición. Pendiente G: Del punto de origen de la escorrentía a una depresión ligera donde se concentra la escorrentía.
El efecto de la topografía en la erosión de suelos en la RUSLE tiene dos componentes, como se apuntó arriba: el factor de longitud de la pendiente (L) y el factor de inclinación de la pendiente (S).
El factor L: Donde λ es la longitud de la pendiente (m), m es el exponente de la longitud de la pendiente y β es el ángulo de la pendiente.1. La longitud de la pendiente se define como la distancia horizontal desde donde se origina el flujo superficial a el punto donde comienza la deposición o donde la escorrentía fluye a un canal definido.
El factor L con el área de drenaje aportadora (Desmet & Govers, 1996)
donde A(i,j)[m] es el área aportadora unitaria a la entrada de un pixel (celda), D es el tamaño del pixel y x es el factor de corrección de forma.
El factor S: El ángulo β se toma como el ángulo medio a todos los subgrids en la dirección de mayor pendiente (McCool et al, 1987,1989).
El factor C C es el factor de manejo de cultivo. El factor C se utiliza para reflejar el efecto del cultivo y las prácticas de manejo en las tasas de erosión. Es el factor que se utiliza más a menudo para comparar los impactos relativos de las opciones de manejo en los planes de conservación. El facto C indica cómo el plan de conservación afectará la pérdida promedio anual y cómo la pérdida potencial de suelo será distribuida en el tiempo durante las actividades de construcción, rotaciones de cultivo u otros esquemas de manejo. El factor C se basa en el concepto de la desviación en relación a un estándar, en este caso un área limpia de cultivos y en barbecho continuo. La tasa de pérdida de suelo (SLR) es entonces un estimado de la tasa de pérdida de suelo en las condiciones actuales en relación a aquella que se experimentaría bajo las condiciones de referencia. “C” representa los efectos de las plantas, la cobertura del suelo, la biomasa del suelo y las actividades disturbadoras del suelo en los procesos de erosión. La RUSLE usa un método del subfactor para calcular las tasas de pérdida del suelo, que son las tasas de pérdida en un momento dado, en una secuencia de manejo de cobertura, en relación a la pérdida de suelo en la parcela unitaria. Las tasas de pérdida del suelo varían con el tiempo, así como con el dosel o estrato, cobertura del suelo, rugosidad, biomasa del suelo y el grado de cambio en la compactación. El valor del factor “C” es una tasa promedio ponderado con respecto al factor R durante el año. Los subfactores utilizados para calcular la razón de pérdida de suelo son cobertura (dosel), cobertura de la superficie y rugosidad, uso anterior del suelo y humedad antecedente del suelo. La cobertura de la superficie es el material que está en contacto con el suelo y que intercepta las gotas de lluvia y detiene o hace más lenta la escorrentía. El porcentaje total de superficie cubierta es la característica utilizada por la RUSLE para calcular cómo la cobertura de la superficie afecta la erosión. La cobertura de la superficie incluye toda la cobertura presente incluyendo fragmentos de roca, vegetación viva, criptógamas y residuos de plantas. El tamaño mínimo requerido de material para ser tomado en cuenta como residuo de la superficie, es que sea de tamaño lo suficientemente grande para no ser removido por la escorrentía.
La RUSLE toma en cuenta la rugosidad de la superficie en el cálculo del valor C. La rugosidad del suelo retiene el agua en las depresiones y reduce la erosividad del impacto de las gotas y del flujo de agua. Si la deposición es lo suficientemente profunda, ocurre mucha deposición en ellos. Con el tiempo, la rugosidad desaparece conforme las depresiones se llenan con sedimento. La rugosidad se reduce en la USLE en función de la lluvia acumulada después de una operación de cultivo. La rugosidad también indica el grado de obstrucción y la probabilidad que la superficie se sellará, aumentando la escorrentía y la erodabilidad del suelo. Ya que se toma en cuenta la rugosidad en los sistemas de siembra, es una de las razones por las que el factor C son menores en la RUSLE. Si un factor C de 0.15 representa el sistema de manejo de un determinado cultivo, esto significa que la erosión será reducida a 15% de la cantidad que hubiera ocurrido bajo condiciones de barbecho contínuo. 1
Cuando se aplica esta fórmula en el Raster Calculator de ArcGIS se debe tomar en cuenta que el ángulo deberá ser convertido a radianes (1 grado sexagesimal = 0.01745 radianes)
A continuación se presentan los valores de C para diferentes cultivos y prácticas.
Otra tabla citada por Jung et al (sf) donde se da una combinación entre el factor C y P, se da a continuación:
Tipo de cobertura Área urbana Suelo desnudo Bosque denso Bosque ralo Bosque mezclado con cultivos (Agroforestería) Cultivos Cultivos de inundación Pastura densa Pastura rala Pasturas con cultivos Humedales Cuerpos de agua
Factor C 0.1 0.35 0.001 0.01 0.1 0.5 0.1 0.08 0.2 0.25 0.05 0.01
Factor P 1 1 1 1 0.8 0.5 0.5 1 1 0.8 1 1
El factor P P es el factor de prácticas de control de erosión. El factor P de la RUSLE refleja el impacto que las prácticas de control tienen sobre la tasa de erosión. Se conceptualiza como al razón entre la pérdida de suelo de un cultivo en contornos y siembra en fajas en comparación con cultivos con surcos rectos pendiente arriba o abajo. Como en el caso de otros factores, el factor P diferencia entre cultivos y pasturas permanentes. Ambas opciones permiten el terraceo o la siembra en contornos, pero los cultivos contienen la opción de la siembra en franjas, mientras que las pasturas permanentes contienen otro tipo de alteración mecánica. Para los propósitos de este factor, la opción de pastura permanente se basa en que las operaciones de control se hacen con poca frecuencia, mientras que en los cultivos son parte de las prácticas anuales de manejo. A continuación se lista una tabla con los valores de P promedio más utilizados.
Objetivos del ejercicio Espacializar los factores componenentes de la ERPUS. Elaborar el mapa de pérdida de suelo por erosión calculado a través de la ERPUS. Contexto del Ejercicio Los datos de este ejercicio corresponden a la microcuenca Arroyo Grande ubicada en el municipio de Jarabacoa, La Vega, Republica Dominicana.
Microcuenca Arroyo Grande Uso de suelo Edafología Estaciones meteorológicas con datos de precipitación media mensual, precipitación total promedio anual, índice de Fournier y factor R. Modelo de elevación Digital Aster (30 mts de resolución)
Abrir el programa ArcGIS-ArcMap y cargar el ejercicio Buscar en el Escritorio el icono y hacer doble click. Si no aparece el icono de ArcMap en el Escritorio puedes encontrarlo en: Inicio-Programas-ArcGIS-ArcMap.
1.
Cargar el documento de mapa llamado Rusle.mxd que se encuentra en el directorio C:\SIG\RUSLE.
Se abre un mapa que contiene las capas descritas anteriormente.
A continuación se describen los ejercicios, siguiendo la secuencia descrita en la metodología 1. Llenado de huecos del DEM (Fill the sinks) Esta función llena los huecos presentes en la malla (o grid). Si celdas con mayor elevación rodean a una celda, el agua es atrapada en esa celda y no puede fluir. La función Fill modifica el valor de la elevación para eliminar este problema. Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-Hydrology-Fill Abra el Arc Toolbox y escoja la opción Spatial Analyst-Tools-Hydrology-Fill. Esta herramienta permite llenar todos los mínimos o depresiones locales que pueda tener el modelo de elevación digital.
Llene la caja de diálogo con los siguientes datos:
Nombre al archivo de salida fill y guárdelo en el directorio de trabajo. Presione el botón OK.
2. Pendiente Ahora calcularemos la pendiente en grados a partir del DEM. Elija la opción del menú de ArcToolbox-Spatial-Analyst Tools-Surface-Slope
Llene la caja de diálogo con los siguientes datos: Nombre a la capa de datos resultante slope y guárdela en el directorio de trabajo.
3. Asignación de Direcciones de flujo (Flow Direction) Esta función calcula la dirección de flujo para una malla dada. Los valores en las celdas de la malla de dirección de flujo indican la dirección del descenso máximo para cada celda. Utilizando las herramientas de hidrología del Arc Toolbox calcularemos ahora el grid de dirección de flujo. Escoja la opción de Arc Toolbox-Spatial Analyst Tools-Hydrology-Flow Direction. Esta herramienta crea un grid de dirección de flujo.
Llene la caja de diálogo con los siguientes datos:
Nombre al archivo de salida FlowDir y guárdelo en el directorio de trabajo. Haga clic en OK.
4. Acumulación de flujo
Esta función calcula la malla de acumulación de flujo. Esta malla contiene el número acumulado de celdas localizadas aguas arriba de la celda. Para cada celda de la malla de entrada. Escoja la opción de Arc Toolbox-Spatial Analyst Tools-Hydrology-FlowAcumulation. Esta herramienta crea un grid de acumulación de flujo. Llene la caja de diálogo con los siguientes datos: Nombre al archivo de salida FlowAcc y guárdelo en el directorio de trabajo. Haga clic en OK.
5. Elaboración del mapa del factor K (erodabilidad del suelo) Ahora se elaborará el mapa del factor de erodabilidad del suelo tomando en cuenta los siguientes factores: Textura del suelo (fracciones) Permeabilidad Estructura Haga click derecho sobre la capa Tipo de suelos y escoja la opción Properties. Haga click en la pestaña Fields y desmarque todas las columnas a excepción de las columnas de Area, Textura, Drenaje, Arcilla_, Arena_gr, Arena_fin y Limo. Haga clic en Aceptar. Ahora exportaremos la tabla de suelos para trabajarla en Excel y así poder hacer el cálculo del factor K. Haga click derecho sobre el nombre de la capa Tipo de Suelo y abra la tabla haciendo click en la opción Open Attribute Table.
Cuando se abra la tabla haga click en el botón Options y luego en la opción Export
Guarde la tabla con el nombre factor_K.txt (text file).
Haga click en Save y luego en OK. No agregue la tabla al mapa todavía.
Abra Excel, navegue al directorio de trabajo y cargue el archivo factor_k.txt. Recuerde que el Tipo de archivo es Text Files.
En la siguiente ventana haga clic en Next.
En la siguiente ventana escoja Semicolon (punto y coma) como el separador de columnas. Haga clic en Finish.
Ahora calcule en una columna aparte la suma de las columnas Arena_fin + Limo_. Cree una columna para él % de Materia Orgánica (MO), el código de estructura y el de permeabilidad. Llene la tabla con los siguientes datos (asumiremos 3% de MO y estructura 3 para todos los datos). Haciendo uso del nomograma provisto calcule el valor de K (para ganar tiempo ya se le proveen en la tabla anterior). Recuerde que estos valores son proveidos en unidades del sistema Inglés y por lo tanto deberá convertirlos a valores del SI multiplicando por el factor 0.1317 (FACT_K_ (SI)).
TEXTURA
DRENAJE
ARCILLA
ARENA_GR
ARENA_FIN
LIMO
LIMO_ARENFIN
MO
ESTRUC
PERMEAB
0
0
0
0
0
3
3
0
35
10
30
25
55
3
3
0 35
0 10
0 30
0 25
0 55
3 3
25
5
55
15
70
0 0 38
0 0 7
0 0 35
0 0 20
0 0 0 0 35
0 0 0 0 10
0 0 0 0 30
25
5
Modarado
25
Modrado
55 0 0 0
6 Franco arcillosa
Moderado a bajo
Franco arcillosa Franco arcillo arenosa
Moderado a bajo Bueno a excesivo
Franco arcillosa
Franco arcillosa Franco arcillo arenosa Franco arcillo arenosa Arcillosa
Moderado
Moderado a bajo Moderado
FACT_K
FACT_K_(SI)
4
0.38
0.050046
3 3
0 4
0.38
0.050046
3
3
2
0.52
0.068484
0 0 55
3 3 3
3 3 3
0 0 3
0.31
0.040827
0 0 0 0 25
0 0 0 0 55
3 3 3 3 3
3 3 3 3 3
0 0 0 0 0
0.38
0.050046
55
15
70
3
3
4
0.6
0.07902
5
55
15
70
3
3
3
0.56
0.073752
5 0 0 0
30 0 0 0
10 0 0 0
40 0 0 0
3 3 3 3
3 3 3 3
3 3 3 0
0.4
0.05268
Ahora transferiremos los datos a la tabla del shape de Tipos de suelo. Abra la tabla del shape Tipos de suelo. Haga click en el botón Options y seleccione la opción Add Field para agregar una nueva columna a la tabla.
Nombre a la nueva columna factor_k y escriba una Precision de 8 (número de dígitos máximo que puede tener el número) y una Scale de 6 (el número de dígitos a la derecha del punto decimal). Luego haga click en OK.
En el menú Editor haga click en Start Editing.
Si aparece esta caja de diálogo, seleccione como el directorio para editar el resaltado en la siguiente figura y luego haga clic en OK.
Copie los valores de K a la tabla como se muestra abajo.
Haga click de nuevo en el menú Editor y luego haga clic en Stop Editing y guarde los cambios hechos en la tabla. Luego cierre la tabla. Haga clic en Ok.
Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-InterpolationIDW
Llene la caja de diálogo con los siguientes datos:
Nombre a la capa de datos resultante Factor_K y guárdela en el directorio de trabajo. Haga click en OK.
6.
Ahora procederemos a calcular el valor de “F” en la fórmula para calcular el factor L:
Elija la opción del menú de ArcToolbox Spatial Analyst Tools-Map Algebra-Raster Calculator
Ingrese la siguiente fórmula: ((Sin("slope" * 0.01745) / 0.0896) / (3 * Power(Sin("slope" * 0.01745 ),0.8 ) + 0.56)) , haciendo click en las teclas de la Calculadora ráster (no use las teclas del teclado) llene los campos, presione Ok. Se agregará un grid llamado factor_f a la tabla de contenidos.
Hacemos doble clic en la capa factor_f
Abra nuevamente el Raster Calculator e introduzca la fórmula "factor_f" / 1 + "factor_f") para calcular el exponente m de la fórmula, llene los campos como se muestra en la figura y presione Ok.
7. Factor L (Factor de longitud de la pendiente) Ahora procederemos a calcular el valor de “L” Para calcular el factor L:
Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-Map AlgebraRaster Calculator Abra Raster Calculator e ingrese la siguiente fórmula para calcular el valor de L. (Power(("flow_acc" + 625),("factor_m" + 1)) - Power("flow_acc",("factor_m" + 1))) / (Power(25,("factor_m" + 2)) * Power(22.13,"factor_m"))
Elija un destino donde se guardara el resultado y nómbrelo factor L. En esta fórmula se utiliza el grid de acumulación de flujo ([flowacc]), la dimensión del pixel (30m x 30m) y el exponente m.
8. Factor S Ahora se calculará el valor del factor S. Ya que la fórmula varía de acuerdo a dos condiciones, se utilizará el comando CON que es la función Condicional de ArcGIS (puede consultar la ayuda del programa si desea profundizar en la sintaxis del comando). La sintaxis simplificada del comando es la siguiente: Con(
, , ) Para calcular el factor S: Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-Map Algebra-Raster Calculator
Abra Raster Calculator e ingrese la siguiente fórmula para calcular el valor de S. Con((Tan("slope" * 0.01745) < 0.09),(10.8 * Sin("slope" * 0.01745) + 0.03),(16.8 * Sin("slope" * 0.01745) - 0.5))
Elija un destino donde se guardara el resultado y nómbrelo factor S.
9. A continuación calcularemos el Factor LS combinado Para calcular el factor LS: Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-Map Algebra-Raster Calculator Abra Raster Calculator e ingrese la siguiente fórmula para calcular el valor de LS. "factor_l" * "factor_s" Elija un destino donde se guardara el resultado y nómbrelo factor LS.
10. Elaboración del mapa del factor R (erosividad de las lluvias) Para la elaboración del mapa del factor R se hará una interpolación de los valores de R calculados para cada una de las estaciones existentes en y alrededor del área de estudio.
Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-InterpolationSpline Llene la caja de diálogo con los siguientes datos:
Nombre a la capa de datos resultante Factor_R y guárdela en el directorio de trabajo.
Si desea profundizar en los métodos de interpolación puede consultar la página: http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=An_overview_of_the_Interpolation_tools
11. Ahora vamos a crear el grid que contendrá el mapa del factor C. Abra el ArcToolBox y seleccione la opción Conversion Tools|To Raster|Feature to Raster
Llene la caja de diálogo que se abre, de la siguiente manera: Nombre a la capa de datos resultante Factor_C y guárdela en el directorio de trabajo. Haga click en OK.
12. Elaboración del mapa del factor A (Pérdida estimada promedio de suelo en Ton/ha/año) Para calcular el factor A: Elija la opción del menú de ArcToolbox-Spatial Analyst Tools-Map Algebra-Raster Calculator
Abra Raster Calculator e ingrese la siguiente fórmula para calcular el valor de A. "factor_r" * "factor_k" * "factor_ls" * "factor_c"
Elija un destino donde se guardara el resultado y nómbrelo Pérdida de suelo.
Bibliografía Abbar Mushtak, T. 2003a. Application of GIS to estimate soil erosion using RUSLE. J. Geo-Spatial Inform. Sci., 1: 34-37. Almoza, Y. et al. 2007. Metodología para la estimación del factor erosividad de las precipitaciones en el modelo (RUSLE). Zonas Áridas, 11(1). 150-158. Barrios R., A.G. Distribución Espacial del factor LS (RUSLE) usando procedimientos SIG compatibles con IDRISI. Aplicación en una cuenca andina. 2000. Revista Forestal Venezolana 44(1): 57-63. Calvo. C. 1993. Dinámica, uso apropiado y sostenible de la tierra en la cuenca del Río Turrialba. Tesis Mag. Sc. Turrialba, CR, CATIE. 186 p. Colotti, E. 2004. Aplicabilidad de los datos de lluvia horaria en el cálculo de la erosividad. Venezuela. Fondo Editorial de Humanidades y Educación, Universidad Central de Venezuela. 182 p. Desmet, P.J.; Govers, G. 1996a. A GIS-procedure for the automated calculation of the USLE LS-factor on topographically complex landscape units. J. Soil and Water Conservation, 51: 427-433. Di Leo, C.M. et al. 1999. Erosividad de las precipitaciones en Tandil, Provincia de Buenos Aires. Ciencia del Suelo. 17 (2): 58-61 FAO. Food and Agriculture Organization. 1992. Erosión de suelos en América Latina. Consultado: 3/11/2008. Disponible en línea: http://www.fao.org/docrep/T2351S/T2351S00.HTM#contents INETER. Instituto Nicaragüense de Estudios Territoriales. 2005. Erosión Hídrica. Mapas de Amenazas. Recomendaciones técnicas para su elaboración. 52 p. Mannaerts, C. 1999. Factores de erosión. Módulo 11: Degradación de suelos. ITC. Curso de Postgrado en levantamiento de recursos hídricos. Notas de clase y ejercicios. CLAS. Cochabamba, Bolivia. McCool, D.K. et al. 1987. Revised slope steepness factor for the Universal Soil Loss Equation. Transactions of ASAE 30 (5): 1387-1396.
McCool, D.K. et al. 1989. Revised slope length factor for the Universal Soil Loss Equation. Transactions of ASAE 32: 1571-1576. Molnar, D. K., Julien, P.J. (1998) Estimation o f u p l a n d e r o s i o n u s i n g G I S . Computers & Geosciences, 24, 183-192. Mora, I. 1987. Evaluación de la pérdida de suelo mediante la Ecuación Universal (EUPS): Aplicación para definir acciones de manejo en la cuenca del río Pejibaye, Vertiente Atlántica, Costa Rica. Tesis Mag. Sc. Turrialba, CR, CATIE. 104 p. Moore, I., Burch, G. (1986). Physical basis of length-slope factor in the universal soil loss equation. Soil Science Society of America Journal, 50, 1294-1298 Oñate V., F. 2004. Metodología para la evaluación del riesgo de erosión hídrica en zonas áridas y su aplicación en el manejo y protección de proyectos hidráulicos. Revista Electrónica de la REDLACH, 1(1). Patriche, C.V. et al. Aspects regarding soil erosion spatial modeling using the USLE/RUSLE within GIS. 2006. Geographia Technica, No.2. p 87-96. Pérez A., J.D.; Mesa S., O.J. 2002. Estimación del factor de erosividad de la lluvia en Colombia. In: Simposio Latinoamericano de Control de Erosión. s.p. Renard, K.G. et al. 1997. Predicting Soil Erosion by Water: “A Guide to Conservation Planning with the Revised Universal Soil Loss Equation (RUSLE)”. USDA, Agricultural Research Service, Agricultural Handbook, No. 703. 384 p. Vahrson, W.G. 1990. El potencial erosivo de la lluvia en Costa Rica. Agronomía Costarricense, 14(1): 15-24. Vahrson, W.G. 1991. Taller de Erosión de Suelos: Resultados, comentarios y recomendaciones. Agronomía Costarricense, 15(1/2): 197-203. Wischmeier, W.H; Smith, D.D. 1978. Predicting rainfall erosion losses. USDA, Agriculture Research Service, Agricultural Handbook, No. 537.
