08 CREANDO UN MODELO DE EROSI N H DRICA

CREANDO UN MODELO DE EROSIÓN HÍDRICA
Ulises Padilla-García

1.- Para este ejercicio se requiere tener los shapes de:
Altura (curvas de nivel)
Precipitación (isoyetas)
Suelos (edafología)
Vegetación (uso de duelo y vegatación)
MDE (modelo digital de elvaciones)

- Hay que activar la herramienta 3D Analyst
Tools – Customize – 3D Analyst y Spatial Analyst

2.- Para obtener el MDE es necesario contarcon un shape de curvas de nivel (altura) y:
Caja de herramientas – Index – Topo to raster – Imput (tema con las curvas de nivel) – Feature Layer
(tema) – Field (campo con datos de elevación) – Type (Contour) – Ouput (nombre del raster MDE) –
ouput cell size(tamaño de la celda 20) -OK

3.- Para obtener la Erosión
Factor R

Potencial y Actual hay que obtener el factor R, K, LS y C.

Las ecuacionesusadas para la determinación de R (ver Mapa) fueron obtenidas del estudio de Cortés
et al, (1990), para el estado de Querétaro aplican dos regiones con sus respectivas ecuasiones:
Región noreste:
Región suroeste:

R = 1.2078P + 0.002276P2
R = 3.4880P + 0.000188P2

Donde R es el índice de erosión o erosividad promedio anual en MJ.mm.ha-1.h-1, y P es la
precipitación media anual en mm.

4.- Deltema de precipitación (isoyetas) obtenemos el Factor R, para esto primero hay que calcular la
precipitación media para cada región (en tabla – add field – “pre-prom” – agregar la información
correspondiente).

Ahora se convierte a grid el shape de precipitación:
Caja de herramientas – Index – polygon to raster – Imput (precipitacion) – Value (elegir el campo
sobre el cual se realizara el raste“pre_pro”) –Ouput (nombrar el archivo, cuidando que no sea muy
grande “precip”) – Cell assigment (elegir el método para asignar la celda – generalmente debe de ser
“Maximun area”) – Prioritary (none) – Cell size (20) – ok.

El raster de “precip” representa la precipitación promedio anual de la ecuación, por lo que es necesario
aplicar la fórmula con este shape.
Obteniendo el factor R: de las opciones deEspatial Analys, selecciona “Raster Calculator” para
poder realizar la siguiente operación con el raster de precip: R = 3.4880P + 0.000188P2
En resumen debes tener algo como esto:

Nota: Es muy importante respetar el número de paréntesis para que las operaciones se realicen adecuadamente.

Recuerda que los grid (raster) generados en la calculadora son temporales, esto significa que si lo
eliminasde la vista se pierde. Para guardar el grid, selecciónalo (botón derecho –export data) y en este
caso guárdalo como factor_r.

5.- Factor K
Para obtener el grid del factor K, vamos a utilizar los valores del factor K que dependen del tipo de
suelo y su textura definidos por la FAO.

o

Para generar este raster hay que abrir la tabla y agregar un campo
nuevo llamado k (que sea de tipo doble, conprecisión y escala de
10), al cual se le agregará el valor de K a cada tipo de suelo según
corresponda en la tabla.

Ejemplo del calculo del factor K: la columna llamada “DFIN_CLAVEPRE y DEFIN_CVEP_1”
contiene el tipo y subtipo principal de suelo; en ocasiones hay presencia de más de un tipo de suelo en
este caso también hay que considerarlo “DFN_CVESEC y DFN_CVES_1” y sumarlo al valor
principal.En edafología se representan a través de la formula edáfica general para cada región, en ella se
mencionan los tipos de suelo presentes separados por un +, seguido por una / y un número que señala
la textura (1 gruesa, 2 media y 3 fina). Por lo que el valor de K es igual al promedio de todos los tipos
de suelo presentes en cada unidad edáfica:
Vp / 3, su valor de K sería 0.026 / 3 = 0.0086
L + Hh/ 2, su valor de K sería 0.040 + 0.020 / 2 = 0.03

Una vez que calculaste todos los factores K para cada unidad de suelo, guarda los cambios y convierte
el shape a grid (utiliza el campo del factork para este procedimiento), dale nombre de factor_k.

Grid Factor K

6.- Factor C
Para obtener el factor C, simplemente debes añadir una columna en el shape de vegetación llamado
factor C (ejercicio...