UNIVERSIDAD TÉCNICA DEL NORTE

FACULTAD DE INGENIERÍA EN CIENCIAS AGROPECUARIAS Y AMBIENTALES

CARRERA DE INGENIERÍA EN RECURSOS NATURALES RENOVABLES

ARTÍCULO CIENTÍFICO

“VALIDACIÓN DE LOS MODELOS RUSLE, WEQ Y CLAVOS DE EROSIÓN EN

LA PÉRDIDA DE SUELO DE LA MICROCUENCA LA MERCED, CARCHI”

AUTORAS:
ARROYO SALAZAR LUISA GABRIELA
JARAMILLO AGUILAR JENNIFER LISBETH

DIRECTOR:
MSc. OSCAR ROSALES

ASESORES:

MSc. Gladys Yaguana

MSc. Eleonora Layana
MSc. Paul Arias

Ibarra – Ecuador

2018

Lugar de investigación: Valle del Chota, sector San Vicente de Pusir Microcuenca la Merced.
 “VALIDACIÓN DE LOS MODELOS RUSLE, WEQ Y CLAVOS DE EROSIÓN
EN LA PÉRDIDA DE SUELO DE LA MICROCUENCA LA MERCED,
CARCHI”

DATOS INFORMATIVOS

APELLIDOS: Arroyo Salazar

NOMBRES: Luisa Gabriela

C. CIUDADANÍA: 1003687769

CORREO ELECTRÓNICO: [email protected]

DIRECCIÓN: Ibarra La Esperanza, Barrio Rumipamba 26- 167

FECHA DE DEFENSA DE TRABAJO DE GRADO: 30 de enero del 2019

 “VALIDACIÓN DE LOS MODELOS RUSLE, WEQ Y CLAVOS DE EROSIÓN
EN LA gsipamba
PÉRDIDA DE SUELO DE LA MICROCUENCA LA MERCED,
CARCHI”

DATOS INFORMATIVOS

APELLIDOS: Jaramillo Aguilar

NOMBRES: Jennifer Lisbeth

C. CIUDADANÍA: 1003687769

CORREO ELECTRÓNICO: [email protected]

DIRECCIÓN: San Pablo Av. Chiriboga

FECHA DE DEFENSA DE TRABAJO DE GRADO: 30 de enero del 2019

“VALIDACIÓN DE LOS MODELOS RUSLE, WEQ Y CLAVOS DE EROSIÓN EN LA
PÉRDIDA DE SUELO DE LA MICROCUENCA LA MERCED, CARCHI”
Armando Oscar Rosales Enríquez1, Luisa Gabriela Arroyo Salazar2 y Jennifer Lisbeth Jaramillo Aguilar3
1Universidad Técnica del Norte-Instituto de Posgrado. Ave. 17 de Julio, 5-21. Ibarra, Ecuador. C.P.:199.
2Universidad Técnica del Norte-Carrera de Ingeniería en Recursos Naturales Renovables. Ave. 17 de Julio, 5-21. Ibarra,
Ecuador. C.P.: 199.
3Universidad Técnica del Norte-Carrera de Ingeniería en Recursos Naturales Renovables. Ave. 17 de Julio, 5-21. Ibarra,

Ecuador. C.P.: 199.

RESUMEN
El valle del río Chota presenta geo-formas planas, inclinadas y muy inclinadas, que sumadas a los factores de
precipitación, viento y uso inadecuado del suelo han incidido en la pérdida de este recurso. El objetivo del presente
trabajo fue evaluar la pérdida de suelos de la microcuenca La Merced localizada en la provincia del Carchi al norte
del Ecuador. Para esto, se utilizó una metodología directa de medición de la erosión del suelo denominada clavos
de erosión, la misma que fue relacionada con los modelos RUSLE y WEQ con una resolución espacial de 3 m en
los raster respecto de los parámetros de los modelos. Se analizaron los valores obtenidos en campo de la pérdida
de erosión del suelo por medio de tres cuadrantes de 10x10 m de clavos de erosión ubicados en la parte alta,
media y baja de la microcuenca, durante el período de un año y dos meses para incluir la época seca y lluviosa.
Además, se utilizó información climática: precipitación, temperatura y velocidad del viento e imágenes satelitales
(TRMM 3B42) y muestreos de Sólidos Secos Totales (SST) y Sólidos Totales Disueltos (STD) y medición de
caudales mensuales en el punto de aforo de la microcuenca. Se utilizó Sistemas de Información Geográfica (SIG)
para el análisis de la distribución espacial de las diferentes propiedades del suelo de la zona de estudio. Para el
modelo RUSLE se obtuvo valores de 494,192 t/ha/año y 159,94 t/ha/año para el modelo WEQ, y un resultado
final de 654,13 t/ha/año de pérdida de suelo y mediante la técnica de clavos de erosión y (SST) (STD), fue de
642,73 t/ha/año. El análisis de correlación de Nash entre las dos metodologías dio como resultado de la validación
un valor de r = 0,99 para el cuadrante uno y r=0,22 para el cuadrante 3 representando una categoría calificada
como Muy Bueno y Satisfactorio.

Palabras claves: Erosión, Suelo, Clavos de Erosión, RUSLE, WEQ, Validación

PALABRAS CLAVE: oso andino, conflicto, peligro de extinción, ganado vacuno, hábitat natural

ABSTRACT
n
The Chota Valley presents plain, tilted and very tilted geological formations. In addition with precipitation and wind
factors, and inadequate use of the soil influenced the wearing of this natural resource. The aim of the present study
was to evaluate the soil degradation in the La Merced micro-watershed, located in the Carchi porovince, at the
North of Ecuador. In order to achieve this, a direct erosion measurement technique was used, called erosion nails.
This technique was related to RUSLE and WEQ models, with an spatial resolution of 3 m of the parameters of the
models. The obtained values of soil erosion were analyzed in three quadrants of erosion nails with 10 x 10 m each.
These were situated in high, medium, and low altitude zones of the micro-watershed during a year and two months
long time period, for including the dry and wet seasons. Furthermore, multiple climate data was considered such
as: precipitation, temperature, wind flow velocity, satellite imagery, total dry solids (Sólidos Secos Totales, SST)
and total disolved solids (Sólidos Totales Disueltos, STD). Geographical Information Systems (GIS) were employed
for the spatial distribution analysis of the zone of study. For the RUSLE and WEQ models, an amount of 494,192
and 159,94 ton/ha/year respectively were obtained, with a final result of ton/ha/year of soil loss. By the other hand,
applying the erosion nails technique and SST & STD, the obtained result was 606.38 ton/ha/year of aparent
erosion. The Nash correlation analysis between both methodologies yielded r = 0.99 for the one quadrant and r=
0,22 for the three quadrant representing a category qualified as Very Good and satisfactory.

Keywords: Erosion, Soil, Erosion Nails, RUSLE, WEQ, Validation

INTRODUCCIÓN Hasta el presente, se han desarrollado una
variedad de metodologías para el estudio de
El suelo es un recurso natural no renovable por
erosión, algunas están basadas en observación
tal razón su degradación afecta severamente a
directa en campo y modelación de factores,
las funciones naturales que son base de la vida, conocidos como métodos directos e indirectos.
al igual que a la regulación del ciclo hidrológico y
Los métodos directos son estudios técnicos, que
a la calidad ambiental, generando daño severo a miden directamente los valores de erosión, como
todo un ecosistema, destruyendo el habita de la
por ejemplo técnica de clavos de erosión
parte superficial de la corteza terrestre (Toy,
(Morgan, 1994). Mientras que los métodos
Foster y Renard, 1999 ; Park, Jeon, Jung y Choi,
indirectos son modelos matemáticos (RUSLE Y
2011). Los procesos de degradación física de los
WEQ), que estiman las pérdidas de suelo a
suelos se manifiestan a través de la erosión
través de variables climatológicas y
hídrica y eólica que se considera comúnmente
características propias de cada tipo de suelo
como una de las principales causas de la pérdida
(Pizarro y Curtiño, 1997; Vega y Febles, 2005).
de suelo (Ostovari, Ghorbani-Dashtaki, Bahrami,
Naderi & Melo, 2017). El concepto de erosión del La presente investigación determinó la tasa de
suelo también puede referirse a la erosión erosión por acción hídrica y eólica de la
antrópica, que es de desarrollo rápido en microcuenca La Merced y se clasifico según los
comparación con la erosión natural o geológica rangos de clasificación de pérdida de suelo
de evolución muy lenta (Mutchler, 1988). propuesta por la FAO (1981).
Según la FAO (2008), las consecuencias de la
degradación del suelo incluyen: disminución de
MÉTODO
la productividad agrícola, migración humana y de
especies animales, inseguridad alimentaria, El área de estudio se localiza en la parte norte
daños a los ecosistemas básicos y pérdida de del Ecuador en la provincia del Carchi, Cantón
biodiversidad. Bolívar; la Microcuenca La Merced posee una
superficie de 337,77 hectáreas (Figura 1).
En el valle del río Chota son evidentes los
Pourrut (1995), define como un clima Ecuatorial
procesos de erosión de origen antrópico y natural
Mesotérmico con temperaturas medias anuales
(precipitación y viento) que ocasionan el
que fluctúan entre 12-20 °C con mínima variación
desgaste del horizonte superficial de los suelos y
entre los meses ecológicamente secos y lluvioso;
la pérdida de fertilidad de este recurso. La
precipitaciones anuales menores a 500 mm.
microcuenca La Merced es una pequeña cuenca
afluente de la cuenca del río Chota
Un indicador del proceso erosivo de la
microcuenca constituye la textura gruesa del
suelo (franco-arenosa) que influye en la
susceptibilidad de pérdida por erosión hídrica y
eólica en la acumulación del material en los
cauces de quebradas situados en las cotas más
bajas del cauce principal de la microcuenca. El
material depositado proviene de la parte alta y
está constituido por partículas de arena, limo y
arcilla (Agrocalidad, 2014). Con el transcurso del
tiempo, los sistemas agropecuarios se han
deteriorado debido al mal manejo de los recursos
suelo y agua, generando acumulación de
sedimentos en las vías de tercer orden que Figura 1. Ubicación
conecta los sectores de Pusir Grande y Tumbatu. La Merced presenta una topografía muy variada
entre plana, ondulada y ligeramente ondulada
que influye en el proceso erosivo; por este motivo utilizó la RUSLE propuesta por Renard et al.,
se optó por delimitar los drenajes con mayor (1991), ver ecuación 1.
representación en el área para poder tener un
Ec 1. Erosión Hidríca
mejor análisis de erosión hídrica, y se definió
como 7 microcuencas (Figura 2). A= R*K*L*S*C
Donde: A es la pérdida de suelo por erosión
hídrica en t/ha/año, R es el factor de erosividad
en (MJ.mm)/(ha.h.año), K es el factor de
erodabilidad del suelo en (t.ha.h)(ha.Mj.mm), LS
es una combinación de la longitud y grado de la
pendiente (adimensional), C es el factor
cobertura del suelo (adimensional).

-Factor de erosividad (R)

El factor R determina la agresividad o erosividad
de la lluvia, define la capacidad potencial de las
precipitaciones para provocar erosión en el suelo
durante un periodo de lluvia determinado
(Renard et al., 1997).
Figura 2. Microcuenca la Merced Se realizó mediante dos tipos de datos:
Estaciones meteorológicas: Se utilizaron datos
Para evaluar la pérdida de suelo en la de precipitación mensual de todos los 15 años de
microcuenca La Merced se utilizó las ecuaciones registros meteorológico de cinco estaciones
RUSLE Y WEQ. Para la obtención de los aledañas a La Merced (Tabla 1) y se utilizó la
parámetros utilizados en las ecuaciones se ecuación 3. Los datos de precipitación se
empleó Sistemas de Información Geográfica interpolaron mediante el método de Kriging el
(SIG), Instituto Geográfico Militar (IGM) e Instituto cual es un proceso geoestadístico que cálculo un
Nacional de Meteorología e Hidrología (INAMHI). promedio ponderado de las observaciones
Posteriormente se tomaron datos en campo muéstrales y hace predicción en sitios no
durante el periodo de 14 meses desde Mayo muestreados dentro del área de estudio.
2017 a Junio 2018, medidos mediante tres
parcelas situadas en el área de estudio, mediante Tabla 1. Estaciones Meteorológicas
la técnica de Clavos de Erosión. Además, en esta CÓD ESTACION COOR COOR ALTITUD
X Y msnm
fase de campo se realizaron análisis de
laboratorio para determinar los sedimentos secos M104 Mira 830900 10061205 2430
y totales del recurso hídrico del cauce principal. M314 Ambuqui 832801 10046983 1833
Finalmente, se validaron entre los datos M597 Piman chico 828105 10048661 1850
obtenidos en campo en el periodo Mayo 2017- M603 Salinas 819850 10055320 1730
Abril 2018; y los valores de erosión calculados M0086 San Vicente 829454 10054693 1891
con las ecuaciones del modelo RUSLE de Pusir

Determinación de la erosión hídrica y eólica. Imágenes Satelitales: Para la validación del

modelo RUSLE, fue necesario descargarse las
Ecuación Universal de Pérdida de Suelo imágenes satelitales del servidor
Revisada (RUSLE) https://mirador.gsfc.nasa.gov específicamente
Para la determinación de la erosión hídrica se del satélite Tropical Rainfall Measuring Mission
laminar y en surcos de terrenos agrícolas se (TRMM 3B42) ver figura 3, las cuales fueron
procesadas en ArcMap 10.4; finalmente se La ecuación 5 fue desarrollada por Moore y
extrajeron los valores con la herramienta de Burch (1986); y se basada en la acumulación del
extract by points a los puntos de ubicación flujo de caudal y la inclinación dela pendiente
geográfica de las estaciones meteorológicas y se (longitud de la pendiente y pendiente inclinada)
utilizó la ecuación 3 de Fournier (1960).
Ec 5. Longitud del Terreno (LS)
LS= ((flow accumulation ×cell size) /22,13) ^0,4
×(sinslope/0,0896)^1,3

Donde LS es longitud del terreno, Flow

accumulation es el número de celdas que
contribuyen al flujo de una celda dada, Cell size
es la longitud del tamaño de las celdas, 22.13 es
la constante, Sinslope es el seno de la pendiente
en radianes.
Figura 3. Imágenes del satélite Tropical Rainfall
Measuring Mission (TRMM 3B42)
-Factor Cobertura Vegetal ( C )

Ec 3. Índice Modificado De Fournier Fue necesario adquirir una ortofoto del Cantón
Bolívar del proyecto SIGTIERRAS del año 2011
IMF= p^2/P con resolución espacial de 0,29 cm, que
mediante la técnica de teledetección se procedió
Donde p^2 es la precipitación del mes más
a dibujar a detalle los polígonos de uso y
lluvioso y P promedio anual de precipitación.
cobertura vegetal, y se obtuvo ocho tipos de uso
de suelos; y se asignaron valores de las Tablas 8
-Factor Erodabilidad (K) y 9 de los estudios realizados por Wischmeier y
Smith (1978), y posteriormente se transformó de
El factor K se describe como la susceptibilidad, formato shaperfile a raster.
vulnerabilidad y pérdida principalmente por la
acción cinética de la lluvia. Y se obtuvo mediante
la ecuación 4 y 5 de Wischmeier y Smith (1978).
Ecuación de Pérdida de Suelo por Erosión
Eólica
Ec 4. Erodabilidad
Para el cálculo de la erosión eólica se utilizó la
K=(2,1 M^(1.14 ) 〖(10〗^(-4))(12-MO)+3.25 Ecuación 6 WEQ desarrollada por Woodruff y
Siddoway (1965).
(s-2)+2.5 (p-3))/100
Ec 6. Erosión Eólica
Ec 5. Tamaño de partículas
E= (I*K*C*L*V)
M=(% limo+% arena muy fina)(100-%arcilla)
Donde E es la pérdida de suelo por erosión eólica
en t/ha/año, I es el factor de erodabilidad del
Donde K es el factor de erodabilidad del suelo
suelo (t.ha.h)(ha.Mj.mm), K es el factor de
(t.ha.h) / (ha.Mj.mm), MO = Materia orgánica, S
rugosidad del surco del suelo (adimensional), C
es el código de la estructura, p es el código de
es el factor climático (adimensional), L es el
permeabilidad y M es el producto de las
factor longitud del terreno en la dirección
fracciones del tamaño de las partículas
prevaleciente de los vientos (adimensional) y V =
es el factor equivalente de cobertura de
vegetación medida en (kg/ha).
-Factor de Longitud del Terreno (LS)
- Factor Erodabilidad (I) parámetros de vegetación, altura y cobertura del
follaje (área foliar) son esenciales para evaluar la
En el Índice de erodabilidad del suelo se utilizó la
erosión eólica (Wischmeier y Smith, 1978). El
capa del factor K obtenida en la ecuación RUSLE
factor vegetal se determina por la cantidad de
cobertura vegetal y los residuos del cultivo en
-Factor Climático ( C ) kg/ha. El área de estudio se encuentra
clasificada en la estratificación vegetal como
Está representado por la velocidad del viento, Bosque Seco Andino (BSA) y FRA, (2015) cita
humedad superficial del suelo en un lugar a Evaluación Nacional Forestal, (2013) el cual en
determinado; Panebianco y Buschiazzo (2008) y su estudio se calculó la biomasa total para BSA
Buschiazzo y Zobeck (2008) desarrollaron la con rangos de 1200-1700 mm de precipitación y
ecuación 7. con más de 5 meses secos, características
Ec 7. Factor Climático geográficas y climáticas de la microcuenca,
otorgando un valor de 95.82 t/ha, dato de
C= 3,86 [U^3/((P⁄2,54)/(1,8 T+ 32))^(10/9) ] entrada que fue convertido a kg/ha, este valor se
asignó al mapa de cobertura vegetal obtenido del
factor C del modelo RUSLE, asignando un valor
Donde, C es el factor climático, U es la velocidad
de biomasa en la categoría de vegetación
media del viento y T es la temperatura media.
xerofítica para posteriormente multiplicar con el
área, y transformar el shapefile a formato raster.
-Factor Rugosidad (K)
El índice de rugosidad superficial fue
Clavos de erosión
determinado por la influencia de crestas y surcos
creados por labores de labranza, estos tipos de Ubicación de cuadrantes
rugosidad reduce la velocidad del viento,
actuando como trampas para las partículas La metodología implantada en este estudio fue
transportadas por acción eólica Zingg y observar áreas evidentes de desgaste de la capa
Wooddruff (1951) y se representa de la siguiente superficial de la corteza terrestre, causado por
ecuación 8: factores antrópicos y naturales, se instaló tres
cuadrantes de 10 por 10 m (100 m2)
Ec 8. Factor de Rugosidad georreferenciados con coordenadas UTM como
lo muestra la figura 1. Posteriormente se
K=4 〖 H〗^2 / I procedió a tomar muestras de los tres cuadrantes
de suelo con tubos de PVC (diámetro de 3 cm)
Donde, H es la altura de los surcos en mm e I es con la técnica adaptada del barreno para ser
la distancia de los surcos en mm. trasladados al Laboratorio de la Agencia
Ecuatoriana de Aseguramiento de Calidad del
Agro (AGROCAILIDAD).
-Factor Longitud en Dirección de los Vientos
(L)
Para la realización la ecuación WEQ se Caracterización de los cuadrantes
consideró la longitud del campo en dirección La metodología consistió en medir y cuantificar el
prevaleciente de la erosión eólica según proceso de erosión hídrica por factores naturales
(Woodruff y Siddoway, 1965). y antrópicos, a través de 3 cuadrantes con la
técnica de clavos de erosión, correspondiente a
un método experimental sencillo, directo y de
-Factor Cobertura de Vegetación (V) bajo costo, por lo cual es considerado el más
La erosión eólica afecta a la superficie del suelo preciso en cuanto a la estimación de pérdida de
sin cobertura vegetal por lo que es necesario los suelo sobre un terreno (Pizarro y Cutiño, 1997).
El material utilizado fue clavos de 6 pulgadas,
colocados a 1 metro de distancia entre clavo a Toma de muestras para análisis de sólidos
clavo, alcanzando un total de 96 clavos. totales, disueltos
Medición con la técnica clavos de erosión Se estableció tres criterios para el punto de
muestreo accesibilidad, representatividad y
El trabajo de investigación fue instalado en el
seguridad, según el Instituto Ecuatoriano De
mes de Abril 2017 y posteriormente se realizó la
Normalización (2016). Para la metodología
primera toma de datos en Mayo 2017 hasta
empleada en la investigación se tomó en cuenta
Junio del 218 realizado durante 1 año 2 meses,
la norma INEN 1105. La metodología en cuestión
en periodos mensuales permitiendo obtener un
se realizó “in situ” permitió cuantificar el nivel de
promedio de pérdida se suelo. El registro de
partículas de suelo que se encuentran en
mediciones consistió en utilizar un calibrador pie
suspensión presentes por acción de arrastre
de rey digital que permitió una medición exacta
natural del recurso agua hacia el punto de aforo
de pérdida de sustrato, para lo cual se colocó el
de la microcuenca La Merced en el periodo de
calibrador en la rodela o arandela hasta la parte
julio 2017- junio2018.
superior del clavo obteniendo el valor correcto de
pérdida de sedimentos.
RESULTADOS Y DISCUSIÓN
Cuantificación de la erosión hídrica La microcuenca presenta una aparente
superficial susceptibilidad ante los factores de erosión
hídrica y eólica registrando pérdidas
Para obtener resultados de erosión hídrica se
significativas de suelo, esto se comprobó
procedió a cuantificar los datos obtenidos en
mediante el análisis espacial de cada factor de
campo, mediante la técnica clavos de erosión
los modelos RUSLE y WEQ, en los cuales se
para lo cual se realizó el cálculo del cociente
utilizaron datos de precipitación del INAMHI de
entre la sumatoria de las mediciones del total de
un periodo de 15 años y se obtuvo el valor de
clavos que presentaron erosión y considerando
297,39 t/ha/año (Figura 4) ; y se puede observar
el valor cero a los clavos que mostraron
que en la figura 5 con los datos de precipitación
sedimentación y aquellos que no presentaron de las imágenes satelitales TRMM 3B43 (mayo
variación.
2017- abril 2018) se obtuvo un valor de 451,31
Para obtener la cuantificación del suelo t/ha/año, valores que se ubican dentro del rango
erosionado se tomó en cuenta las unidades t/ha de clasificación Muy Fuerte de pérdida de suelo
para el periodo de estudio, se transformó las propuesta por la FAO (1981).
medidas de espesor o altura de suelo erosionado
de milímetros a metros y se multiplicó por el valor
de la densidad aparente y área de cada
cuadrante en estudio Da (t/m3) Hudson (1976).
La expresión matemática que determina estos
procesos, planteada por Pizarro y Cuitiño (1997)
y posteriormente redefinida por Pizarro y Cuitiño
(2002), se presenta a la ecuación 9.
Ec 9. Suelo sedimentado o erosionado
X (t / ha)= Y*Da*A
Dónde : X es el suelo sedimentado o erosionado
(t ⁄ ha), Y es la altura del suelo sedimentado (m),
Da es densidad aparente del suelo (t /m^3 ) y A
es el área de cada cuadrante (m2).
Figura 4. Pérdida de suelo con datos de INAMHI
La microcuenca cuenta con el 75% de suelos de En el modelo WEQ de erosión eólica Eltaif y
categoría taxonómica de Vitrandic Haplustols Gharaibeh (2011) y Nyamtseren, (2016) exponen
característicos por ser secos en todas partes de que el clima es el principal factor erosivo. Sobre
los horizontes del suelo y con una infiltración la estructura del suelo, se menciona que los
media a alta según USDA (2006), haciéndolos agregados de este recurso pueden ser fácilmente
más propenso a tener escorrentía superficial transportados por el viento, cuando se disminuye
como explica López (2017). En un estudio la humedad en la superficie, especialmente en
realizado por Jaramillo (2015) de la subcuenca suelos de transición de cultivos a barbechos; la
del río Portoviejo un 0,21% de áreas con erosión microcuenca cuenta con 14,7% de suelos
fuerte lo que explica que se debe a sus fuertes desprotegidos de cobertura vegetal dejando
pendientes y superficies que se encuentran con expuestos a los procesos erosivos eólicos e
actividades agropecuarias sin un manejo hídricos. Las pendientes fuertemente
adecuado del recurso suelo pronunciadas tienen cierta influencia en la
erosión eólica, actuando como corredores de
vientos la principal dirección de viento fue NW
(Noroeste); se observó que las áreas con una
llanura inclinada y relieve llano rodeadas de
montañas se distinguen como zonas más
afectadas por la erosión eólica.
En la tabla 2 se muestran los datos obtenidos en
con los modelos de pérdida de suelo por erosión
hídrica (RUSLE) y eólica (WEQ).
Tabla 2. Total de pérdida de suelo de la Merced
RUSLE WEQ Pérdida
total de
t/ha/año t/ha/año suelo

Figura 5. Pérdida de suelo con datos de Microcuenca 297,39 159,94 457,33

La Merced (Datos
imágenes satelitales TRMM 3B43 INAMHI)
El mapa de erosión eólica de la microcuenca La Microcuenca
Merced se obtuvo aplicando la técnica de algebra La Merced 494,192 159,94 654,132
de mapas de todos los factores individuales que (Datos TRMM)
la ecuación WEQ propone. En la figura 6 se
muestra la tasa de la erosión eólica en la
Pérdida de suelo por erosión hídrica mediante
microcuenca con valores entre 0 a 159,92
la técnica de clavos
t/ha/año y con un valor medio de 12,18 t/ha/año.
Se obtuvo valores de 642,73 t/ha/año mediante
la técnica de clavos de erosión y análisis de
sólidos totales y disueltos (figura 7).

Figura 6. Pérdida de suelo por erosión eólica Figura 7. Pérdida de suelo por erosión hídrica
 señala que el estudio elaborado obtuvo una
La técnica de clavos de erosión instalada en la erosión neta de 31,2 t/ha y sedimentación de
microcuenca fue de gran importancia para el 13,5 t/ha, sin duda en el ensayo realizado en la
estudio en cuestión, debido a que se calculó el Microcuenca La Merced en los tres puntos
volumen de suelo erosionado en los tres puntos establecidos por cuadrantes se determinó la
establecidos para la investigación. Por su parte dinámica del suelo en su respuesta a la acción
Moreno (2007) menciona que al estimar la de la precipitación y el escurrimiento superficial.
erosión laminar se puede evidenciar y captar
De igual manera, para proceder a realizar la
valores acertados con la técnica de clavos de
validación de los datos obtenidos en campo se
erosión; tal es el caso que al cabo de 1 año y 2
procedió a la metodología de Zapata (2015),
meses del periodo 2017 y 2018 se estableció tres
menciona que para una evaluación precisa de
cuadrantes en áreas idóneas para el estudio,
variables es necesario utilizar criterios entre los
resultando valores totales en campo de 642,73
cuales se incluye el coeficiente de correlación de
t/ha de erosión aparente.
Nash (1970); por lo que al comparar con las
evidencias obtenidas en campo de la presente
La FAO (1981) señala que los niveles de erosión
investigación se determinó que los valores se
tanto hídrica como eólica, por lo que señala que
encuentran en porcentajes excelentes y muy
los niveles de erosión son clasificados según el
buenos; según la categorización de Nash.
grado de pérdida de suelo, con valores que
sobrepasen o sean mayores que 200 t/ha/año, se El empleo de un metodología para obtención de
encuentran en un nivel de erosión muy fuerte, sin sólidos totales disueltos y sólidos secos totales
duda los valores obtenidos en el estudio se lo asocia con el estudio realizado por Vega y
realizado en el periodo 2017 - 2018 indica que la Febles (2005), menciona que la tasa de erosión
microcuenca La Merced se encuentra en un se asume a la pérdida de suelo directamente
rango de 642,73 t/ha/año ubicándose en la proporcional a la cantidad de moléculas de suelo
escala de erosión hídrica Muy Fuerte. perdido al arrastre y transporte de sedimentos
hacia cuerpos de agua; la microcuenca la Merced
se obtuvo un gasto total de sólidos secos y
Validación de la ecuación RUSLE empleando
los datos de medición de clavos de erosión y sólidos suspendidos totales de 642,73 t/ha/año
por lo que cabe recalcar que los factores tanto
datos de análisis de sólidos totales disueltos
y sólidos secos totales. antrópicos como naturales son fundamentales
para la el proceso de sedimentación y transporte.
La validación del modelo RUSLE y Clavos de
CONCLUCSIONES
erosión; se obtuvo para el cuadrante uno valores
de r2 = 0.99 ubicándose en una categoría de En la determinación de erosión hídrica,
excelente. El cuadrante dos presento un r2 = - mediante el programa ArcMap 10.4 y el
1,11 ubicándose en una categoría de insuficiente modelo RUSLE utilizado en el presente
representando que el método de clavos de estudio, se determinó valores de 494,19
erosión es más recomendable y finalmente para t/ha/año correspondiente a valores
el cuadrante tres se obtuvo un valor de r2= 0,22 categorizados según la FAO para un nivel de
categorizando como valor satisfactorio para el erosión muy fuerte; para el cálculo de
modelo según Nash (1970). erosión eólica mediante el modelo WEQ se
La técnica de clavos en la presente investigación obtuvo valores de 159,94 t/ha/año.
y metodología experimental, es de gran
influencia ya que se considera una técnica En la microcuenca La Merced, la erosión
idónea para la determinación de erosión neta y eólica es un grave problema por sus
suelo movilizado; por lo que menciona Pizarro et características de zonas seca y árida; siendo
al., (2008) en su estudio ejecutado en suelos más alta la afectación en los periodos de
desnudos de la pre-cordillera andina y valle transición de cultivo a barbecho.
central de la región de chile; con esta referencia
El cálculo de erosión con los modelos
RUSLE Y WEQ mediante la utilización de BIBLIOGRAFÍA
cartografía permitió obtener valores Arnoldus, H., M., J. (1977). Methodology used to
significativos de erosión hídrica y eólica, determine the maximum potential average
con un total de 654,13 t/ha/año. Este valor se annual soil loss due to sheet and rill erosion in
relaciona con los datos obtenidos en campo Morocco. FAO Soils Bulletin , (34), 39–51.
utilizando la técnica de clavos de erosión,
análisis de sólidos totales y disueltos cuyo AGROCALIDAD (Agencia Ecuatoriana de
dato fue de 642,73 t/ha/año. Aseguramiento de Calidad de Agro). (2014).
Guía de Buenas Prácticas Agrícolas. vol. 03. pág.
(14-68).
La erosión total estimada para la
microcuenca La Merced, mediante la técnica Eltaif, N. I., y Gharaibeh, M. A. (2011). Aplicación
de clavos de erosión, fue de 606,39 t/ha, este de un modelo matemático para predecir y
valor indica que el área en estudio se reducción de la erosión eólica en tierras áridas no
encuentra bajo un proceso acelerado de protegidas; Revista Chapingo Serie Ciencias
erosión hídrica. Forestales y del Ambiente, (17), 195-206.
Evaluación Nacional Forestal. (2013). Evaluación
La medición de erosión hídrica, tanto por la Nacional Forestal-Resultados. MAE.
técnica de clavos de erosión como aplicando Recuperado de:
el modelo RUSLE, generó valores mayores http://suia.ambiente.gob.ec/documents/10179/1
a 200 t/ha/año que corresponden a un rango 85860/Evaluaci%C3%B3n+Nacional+Forestal_
de pérdida de suelo severa según la FAO NREFD+1.pdf/955aaa38-34b6-4b4d-9278-
(1981). 8fe915df893f

El método de los clavos de erosión es FAO. (1981). EL ESTADO MUNDIAL DE LA

adecuado para cuantificar la erosión hídrica, AGRICULTURA Y LA ALIMENTACIÓN ROMA
producida en un sector determinado 1982; Análisis mundial: situación en los países
conociendo también la cantidad de menos adelantados Cómo aliviar la pobreza
sedimentos producidos en un área de rural. Colección FAO: Agricultura. 14, 204.
estudio, por tal razón es una técnica apta para
establecer medidas de manejo y FAO. (2008). FAO Sala de Prensa. Recuperado
conservación de suelos. el 2 de 06 de 2015, de
La validación de los datos simulados http://www.fao.org/newsroom/es/news/2008/100
(RUSLE) y los datos observados (Clavos de 0874/
erosión) de los tres cuadrantes; se obtuvo FRA (Evaluación de Recursos Forestales)
que el cuadrante uno se encuentra en la (2015). Paquete de informe sobre bosques 2015.
categoría de excelente con un r2 = 0.99. El Recuperado de: http://www.fao.org/3/a-
cuadrante 2 presento un r2 = -1,11 az203s.pdf
representa que el método de clavos de
erosión es más recomendable; hay que tener
en cuenta que la pendiente de cuadrante 2 es Hudson, N. W. (1976). Efficiency of Water and
la más alta de 64,77; y el modelo asigna Fertilizer Use in Semi-Arid Regions; Joint
valores enteros a los clavos que presentan Fao/Iaea Division of Atomic Energy in Food and
sedimentación que en campo se asignó al Agriculture. London-England.
valor de cero. El cuadrante tres obtuvo un
INEN (Instituto Ecuatoriano De
r2= 0,22 categorizando como valor
Normalización).(2016). Censo Agropecuarío.
satisfactorio para el modelo.
Recuperado de:
http://www.ecuadorencifras.gob.ec/estadisticas- Panebianco, J. E., y Buschiazzo, D. E. (2008).
agropecuarias-2/ Erosion Predictions with the Wind Erosion
Equation (Weq) Using Different Climatic Factors.
Jaramillo, J. J. (2015). Estudio del Riesgo por
InterScience. Argentina. Doi:10.1002/ldr.813
Erosión Hídrica del Suelo Utilizando el Modelo
U.S.L.E, Mediante Herramientas S.I.G, Aplicado Pizarro, T. Cutiño, M. (1997). Evaluación
en la Subcuenca Río Portoviejo, Provincia de cuantitativa de la erosión hídrica superficial en
Manabí (Tesis de Maestría). Universidad Central suelos desnudos de la pre-cordillera andina y
del Ecuador. Quito-Ecuador valle central de la VII Región. VI Jornadas del
CONAPHU-CHILW, 6,16.
López, J. (2017). Erosión Hídrica en la Unidad
Hidrográfica del Río Ambato, en un Escenario de Pizarro, T., y Cuitiño, M., H. (2002). Método de
Cambio Climático, Mediante el Modelo Evaluación de la Erosión Hídrica Superficial en
USLE/RUSLE (Tesis pregrado).Universidad Suelo Desnudo en Chile. Cuad.Soc.Esp.Ciencia
Nacional de Chimborazo, Riobamba- Ecuador. Forestal, (13) ,165-170.
Meyer, L. D., y Wischmeier , W. H. (1969). Panebianco, J. E., y Buschiazzo, D. E. (2008).
Mathematical Simulation of the Process of Soil Erosion Predictions with the Wind Erosion
Erosion by Water.Transactions of the ASAE, Equation (Weq) Using Different Climatic Factors.
12(6): 754-758,762.e. InterScience. Argentina. Doi:10.1002/ldr.813
Pourrunt, P. (1995) El Agua en el Ecuador.
Moore, I. D., y Burch, G. J. (1986). Base física del
Instituto Nacional de Hidrología y Meteorología
factor de pendiente longitudinal en la ecuación
(INAMHI), (7). Quito, Ecuador: Corporación
universal de pérdida de suelo. Soil Science
Editora Nacional y Colegio de Geógrafos del
Society of America Journal, 50(5):1294–1298
Ecuador.
Moreno-Brotons, J. (2007). Erosión eólica e
Renard, K. G., Foster, G. A., Weesies, G. A.,
hídrica en estériles de minería en el campo de
McCool, D. K y Yoder,D.,C .(1997). Predicting
Cartagena. PFC Ciencias Ambientales. Murcia.
soil erosion by water: a guide to conservation
Morgan, R. 1., (997). Erosión y conservación de planning with the Revised Universal Soil Loss
suelos. Madrid, España: Editorial Mundi-Prensa. Equation. Agricultural Research Service (USDA-
ARS) Handbook No. 703. US Department of
Mutchler, K .(1988). Soil erosion research
Agriculture: Washington, DC.
methods. Arliton Vriginia. Estados Unidos.
Toy T, J., Foster G. R., y Renard K.G. (1999).
Nash, J E; Sutcliffe, J V (1970). “Previsión del
RUSLE for mining, construction and reclamation
flujo del río a través de modelos conceptuales,
lands J. Soil Water Conservancv, (54),462-467.
parte I - Una discusión de principios”. Revista de
Hidrología. 10 (3): 282 - 290. USDA (Departamento de Agricultura de los
Estados Unidos Servicio de Conservación de
Nyamtseren, M., Jamsran, T., Doljin, D., y Sodov,
Recursos Naturales). (2006). Claves para la
K. (2016). Evaluación espacial de la erosión
Taxonomía de Suelos. Decima Edició.Texaco-
eólica del suelo usando el enfoque WEQ en México.
Mongolia. Revista de Ciencias Geográficas, 26
(4), 473-483. Vega, M., y Febles, J. (2005).La Investigación De
Suelos Erosionados: Métodos E Índices De
Ostovari, Y., Ghorbani-Dashtaki, S., Bahrami, H.,
Diagnóstico. Minería y Geología, 21 (2), 1-18.
Naderi, M., y Melo, J.A.(2017). Soil loss
Holguín, Cuba.
prediction by an integrated system using RUSLE,
GIS and remote sensing in semi-arid región. Wischmeier, W. H. y Smith, D. D. (1978).
GEODRS. 1-44.doi: Predicting rainfall erosion losses a gudie to
10.1016/j.geodrs.2017.06.003 conservation planning. Agriculture Handbook No.
537. USDA/Science and Education
Administration, Washington, DC
Woodruff, N. P., y Siddoway, F. H. (1965). A wind
erosion equation. Soil Science Society of
America Proceedings. Recuperado de:
https://www.ars.usda.gov/ARSUserFiles/302005
25/897%20A%20wind%20erosion%20equation.
pdf
Zingg, A.W., Woodruff, N. P., 1951. Calibration of
a Portable Tunnel for Simple Determination of
Roughness and Drag on Field Surfaces. Agron.
J. 43 (4), 191–
193