1er CONGRESO IBEROAMERICANO SOBRE SEDIMENTOS Y ECOLOGA

QUERTARO, QUERTARO MXICO, 21-24 JULIO 2015

Evaluacin y anlisis de la produccin de escurrimiento y sedimentos

en una ladera con dos condiciones contrastantes en Tzintzuntzan,
Michoacn.
Rubn I. Huerto Delgadillo
Instituto Mexicano de Tecnologa del Agua
[email protected]

Alfredo Amador Garca

Universidad Michoacana de San Nicols de Hidalgo
[email protected]
Introduccin
En Mxico existen cerca de 70 lagos naturales, cuyas
extensiones varan entre 1,000 y ms de 10,000 ha, y cubren
en conjunto una superficie de 370 891 ha de nuestro pas
(Arriaga et al., 2000 citados por Carabias 2005). En la cuenca
del Lago de Ptzcuaro, Mich., el IMTA (2002) revis
distintos aspectos de la carga de sedimentos, estimando
aportes promedio de 23 628 ton/ao. De entonces a la fecha
personal de IMTA, Comisin Forestal de estado de Michoacn
(COFOM) y habitantes locales han atendido mediante la
construccin de obras de conservacin de suelo y agua
(CONSA) como son piedra acomodada, presa de gavin y
geocostales en distintas reas sumamente degradadas de la
cuenca, especialmente en el cerro Tariqueri en el municipio de
Tzintzuntzan, Mich., donde recientemente se ha incorporado
el uso de pasto vetiver, maguey y tinas ciegas.
Desde hace ms de 6 aos, se han conducido en la zona
diversos estudios y obras como parte del Programa de
Recuperacin Ambiental de la cuenca del Lago de Ptzcuaro.
El presente estudio forma parte de dichos proyectos, mismos
que han permitido con sus resultados precisar la magnitud de
las acciones para contener el problema de la carga de
sedimentos que llegan al lago por distintos puntos de la
cuenca, influyendo directamente en su eutroficacin y
proliferacin de vegetacin acutica.

Antecedentes
En las zonas ecuatoriales, las lluvias en general presentan
grandes variaciones estacionales y por consecuencia el
escurrimiento. Todos los ros del mundo no aportan
nicamente agua sino que tambin transportan sedimentos, es
decir, materiales que han sido erosionados aguas arriba en su
cuenca o cauce. De manera que aunque el proceso tome
mucho tiempo, todo embalse natural o artificial llegar a ser
llenado por los sedimentos. Por ello cuando es

importante la planeacin de obras de control de sedimentos

tanto en la cuenca como en sus cauces, as como en los
torrentes y arroyos que descargan directamente en el vaso para
asegurar una prolongada vida til de los mismos

(Campos-Aranda 2007 y 2010). El principal fenmeno que

contribuye al proceso de sedimentacin en la superficie de las
cuencas es la erosin. Se considera a esta como la cantidad
total de material erosionando en la cuenca, incluyendo
perdidas del suelo, socavaciones en cauces y aportaciones de
otras fuentes. En cambio, la aportacin de sedimentos es la
porcin de la erosin total que llegar a un punto de inters, en
el lapso de un ao, por ello se expresa en ton/ao o en m3/ao.
Monsalve (1999) seala que la discusin sobre la
determinacin o clculo de sedimentos en corrientes de agua
sigue siendo un tema abierto a la investigacin dada la
variabilidad en la prediccin de los distintos mtodos
empricos. No obstante, se admite una estrecha relacin con el
caudal o flujo que se presente en el cauce. De hecho lo
recomendado por este y otros autores (Brooks et al. 1991,
Ward y Trimble 2004 y Wanielista et al. 1997) es tener las
medidas observadas de la descarga tanto de agua (caudal)
como de sedimentos a distintas alturas del cauce. Quintanilla
et al. (2005) destacan la importancia de aforar y medir
caudales ya que estos permiten encontrar relaciones con las
cargas de sedimentos as como, en general, establecer
relaciones de este con la calidad del agua.
Para calcular caudales pequeos como los que se forman en
las laderas de las montaas, lo ms sencillo es la medicin
directa. Adems, se han desarrollado varios modelos de
aforadores para ser utilizados en distintas condiciones (Bos
1976, citado por Hudson 1997). As, el Servicio de
Conservacin de Suelos del Departamento de Agricultura de
los Estados Unidos dise un grupo de aforadores especiales
denominados aforadores H para medir los caudales con
exactitud y continuidad a partir de parcelas o de pequeas
cuencas experimentales. Son aforadores que miden caudales
escasos con exactitud, tienen buena capacidad para caudales
elevados, y no requieren estanque de amortiguacin. Adems,

1er CONGRESO IBEROAMERICANO SOBRE SEDIMENTOS Y ECOLOGA

QUERTARO, QUERTARO MXICO, 21-24 JULIO 2015
dan paso a escurrimientos que contienen una fuerte carga de
sedimentos.
De los tres modelos de aforadores en H, el ms pequeo (HS)
puede registrar caudales de hasta 22 l/s, el tipo normal
(H) puede medir caudales de hasta 2.36 m/s y el mayor
(HL) caudales de hasta 3.32 m/s. Cada tipo se puede construir
en diversas dimensiones que se determinan por la profundidad
mxima del caudal (D); las dimensiones de fabricacin se dan
como proporciones de D, pero las proporciones de los lados
del aforador, son diferentes para cada uno de los tres tipos HS,
H y HL. Los aforadores en H pueden funcionar parcialmente
sumergidos. Se suelen prefabricar con lminas de metal y
pueden utilizarse en forma provisional empleando sacos de
arena para formar un canal de acceso o aproximacin o
tambin como instalaciones permanentes utilizando hormign
o mampostera. Los aforadores en H tienen otras dos ventajas.
El agua fluye a travs de la escotadura rpidamente de manera
que no se produce depsito de sedimentos en el aforador. De
acuerdo con Maza (2005) existe una relacin estrecha entre la
carga real de sedimentos que un ro transporta a su paso por
una seccin determinada y la produccin de material de
arrastre en la cuenca por procesos de erosin pluvial y de
remocin en masa, y que la aplicacin de ecuaciones
empricas para calcular o estimar la prdida de suelo en un
perodo dado no produce resultados precisos, pero permite
determinar ndices para comparar los potenciales erosivos en
cuencas de la misma regin.

precipitacin en el rea de estudio con la ayuda del Sistema de

Informacin Geogrfica (SIG) ArcView 3.2. Los datos diarios
faltantes de las estaciones se obtuvieron aplicando una funcin
de probabilidad gamma mixta o log-pearson inversa (CamposAranda 2007).

(a)

Metodologa
Se proces el Modelo Digital de Elevacin (MDE) de la carta
E14A22 (Ptzcuaro) escala 1:50,000 del INEGI y la extensin
DetermHydro desarrollado por Valtierra (2007), para obtener
la morfometra de las unidades de escurrimiento sobre dos
cauces en una ladera de exposicin al Noroeste del cerro
Tariaqueri, en el municipio de Tzintzuntzan dentro de la
cuenca del lago de Ptzcuaro, Mich. Para la estimacin del
escurrimiento fueron diseados y construidos dos vertedores
tipo H (Fotografa 1).
El vertedor H1 correspondi a la obra de aforo de una ladera
con obras de conservacin de suelo y agua (CONSA) mientras
que el vertedor H2 correspondi a la obra de aforo de una
ladera de similares caractersticas biofsicas y morfomtricas
adyacente, pero sin obras CONSA. Los vertedores fueron
equipados con dispositivos automatizados (level loggers de la
marca Solinst) que registraron el nivel de columna de agua
D en cada vertedor. Los datos de escurrimiento obtenidos se
extrajeron del dispositivo a una base de datos de Excel donde
fueron analizados para realizar la relacin lluvia-escurrimiento
mediante el mtodo del Soil Conservation Service (SCS) el
cual considera abstracciones de la precipitacin de cada
tormenta.
Se cont con las bases de datos de precipitacin de los aos
2004-2011 de las estaciones automticas que administra el
IMTA en las localidades de Buenavista, Chupcuaro, isla La
Pacanda y Ptzcuaro. Se determinaron valores puntuales
diarios, mensuales y anuales que se utilizaron para interpolar
(por inverso de la distancia ponderada o IDW) superficies de

(b)
Fotografa N 1. Vertedor en unidad H1 (a) vista desde el interior
del canal de aproximacin ya con sensor de nivel Solinst y (b) vista
en perspectiva.

La carga de sedimentos se midi en las dos laderas

contabilizando el nmero de cubetas de 20 lt de cada vertedor,
obteniendo una muestra de 500 ml para cuantificar el peso seco
y relacionarlo con el volumen retenido y removido
semanalmente de los vertedores. Los resultados obtenidos
fueron capturados en una matriz de datos, con los cuales se
analizaron las estadsticas descriptivas generales para la
aprobacin de la hiptesis de que existen diferencias
significativas en la cantidad de escurrimiento y arrastre de
sedimentos en ambas laderas.
Para estimar la prdida de suelo por erosin, el modelo ms
utilizado para esta prediccin es la Ecuacin Universal de
Prdida de Suelo (EUPS) con base en lo sealado por SARHCP (1991) y con la que se cuantifica la erosin a travs de seis
factores (R, K, L, S, C y P) con la ayuda del SIG ArcView 3.2.

1er CONGRESO IBEROAMERICANO SOBRE SEDIMENTOS Y ECOLOGA

QUERTARO, QUERTARO MXICO, 21-24 JULIO 2015
El factor R o erosividad de la lluvia (MJ mm ha-1 h-1) se basa en
estimar la lluvia anual y aplicar un anlisis de regresin
con la ecuacin que corresponde a la regin V la cual pertenece
al estado de Michoacn, y tiene un coeficiente de correlacin de
0.94 (1):
(1)
Dnde: Y= EI30 anual (MJ mm ha-1 h-1) y X= lluvia anual (mm)
Esta ecuacin se distribuy en ArcView mediante la creacin
de imagen raster. Para estimar el factor K o erodabilidad del
suelo se obtiene a partir de la textura superficial del suelo (% de
limos + arenas muy finas, % arenas), % de materia orgnica,
estructura y permeabilidad; estos valores deben ser promedio
para los 15 o 20 cm superficiales para ser usados en un
nomograma y tablas (SARH-CP 1991). Se clasific con estos
criterios el archivo vectorial y se rasteriz con las mismas
caractersticas que la imagen del factor R.
Para obtener el valor del factor Longitud-Pendiente LS es
necesario conocer la longitud de la pendiente L; se define como
la distancia que existe desde el punto donde se origina el
escurrimiento y hasta donde la pendiente disminuye a tal grado
que ocurre la sedimentacin, o bien, donde el flujo entra a un
cauce bien definido que puede ser parte de una red de drenaje o
un canal construido tal y como el de una terraza o zanja de
desviacin. Se rasteriza con las mismas caractersticas que los
factores R y K.
Para el clculo del factor C se toman valores de suelos
agrcolas, terrenos de pastizal y forestales. Esto debido a que
sus efectos estn influenciados por muchas correlaciones
importantes como: rotacin de cultivos, labores de labranza, los
residuos vegetales se pueden remover o dejarse, eficiencia de
los cultivos para evitar la erosin. Tambin se obtienen de
tablas (SARH-CP 1991). Para su distribucin espacial se
interpretaron escenas Google Earth, se crearon los polgonos de
las cuatro coberturas en las categoras de suelo desnudo, suelo
cultivado, matorral y bosque. Los archivos KMZ se exportaron
a SHP y se rasterizaron con las mismas caractersticas que las
anteriores imgenes. Solo se consider un valor constante de
presencia de prcticas de conservacin en la unidad de
escurrimiento aforada en H1. Por ltimo se obtuvieron los
valores de erosin actual por multiplicacin de los valores de
todas las imgenes obtenidas y se extrajo la geoestadstica para
cada unidad de escurrimiento.

Resultados
Las principales caractersticas morfomtricas de las laderas
estudiadas se presentan en la tabla 1. Los registros de caudal
en los vertedores se ajustan al modelo polinmico emprico
(R2=0.99) en el que x es la altura registrada en el sensor de
nivel e y el caudal en litros por segundo (lps)(2).
y = 0.0005x3 + 0.0609x2 + 0.1027x

(2)

Los registros de los sensores de nivel cada 5 minutos

corresponden a la altura de la cabeza de vertedor (x) que se
utiliza para el clculo del caudal (y) con el modelo sealado.

Los mayores picos de descarga se obtuvieron como resultado

de las tormentas del 24 y 29 de julio de 2011 respectivamente,
con gasto prximos a los 200 lps en el vertedor H1. Para esas
mismas tormentas, el hidrograma observado con los registros
del vertedor H2 los gastos son de aproximadamente 20 lps.
De las 11 estaciones meteorolgicas con pluvimetros
automatizados que dispone el IMTA en la cuenca solo fueron
de utilidad los 4 ms cercanos al rea de estudio (Tabla 2), de
los cuales solo una mnima proporcin de informacin faltante
se complet con el mtodo descrito.
Tabla N1. Morfometra bsica de las unidades de escurrimiento
aforadas en el rea de estudio.

Rasgo
rea (ha)
Permetro (m)
Pendiente minima (%)
Pendiente mxima (%)
Pendiente promedio (%)
Elevacin mnima (m)
Elevacin mxima (m)
Elevacin Promedio (m)
Elevacin (des. est.)
Longitud mxima (km)
Longitud equivalente (km)
Long. Mxima (m)
ndice de forma
Coeficiente de compacidad
Relacin de elongacin
Relacin de bifurcacin
Densidad de drenaje (km/km2)
Densidad de corriente (corr./km2)
Orden de cuenca
Long. Mxima (m)

H1
(con obras
CONSA)
10.3
1,548.4
0
91.3
30.8
2,120.0
2,320.0
2,239.4
42.2
0.77
0.60
682.5
0.22
1.36
0.53
1.50
9.42
58.01
2.0
682.5

H2
(sin obras
CONSA)
18.8
2,388.6
0
83.3
31.1
2,159.1
2,440.0
2,320.2
57.4
1.18
1.00
1,031.25
0.18
1.55
0.48
0.75
8.44
84.81
3.0
1,031.2

Tabla N2. Coordenadas geogrficas (fraccin decimal de grado)

y UTM (datum ITRF92) de localizacin de las estaciones
meteorolgicas y/o pluvimetros automatizados y proporcin de
registros pluviomtricos en 2011.
Estacin
Buena Vista
Chupcuaro
Pacanda
Ptzcuaro

LONG
101.540
101.579
101.649
101.630

LAT
19.556
19.674
19.608
19.547

X
233 496.8
229 601.0
222 102.0
224 001.7

Y
2 164 283.4
2 177 456.6
2 170 260.7
2 163 418.8

%
93.3
82.8
88.6
83.1

Con los valores observados de escurrimiento y carga de

sedimentos, as como la extraccin geoestadstica de las
imgenes raster resultado de interpolacin de los valores
puntuales diarios de las 4 estaciones automticas
mencionadas, se obtuvieron los valores promedio diario de
carga de sedimentos en la unidad de escurrimiento aforada por
el vertedor H1 (con obras CONSA) y la marcha anual de todos
los parmetros medidos y estimados que se muestran en la
Fotografa N1. De igual modo que lo concerniente a la unidad
de escurrimiento aforada por el vertedor H2 (sin obras

1er CONGRESO IBEROAMERICANO SOBRE SEDIMENTOS Y ECOLOGA

QUERTARO, QUERTARO MXICO, 21-24 JULIO 2015
CONSA) en que se muestra su marcha anual en 2011 en la
Ilustracin N 2.

jul

4.0

pp (mm)

3.0

2.0

6
8

1.0

10

0.0
Sed (kg)

rea del
vertedor

sep nov

Q (m3)

escurrimiento (m3) y carga de

sedimentos (kg)

ene mar may

pp(mm)

Ilustracin N1. Fluctuacin mensual de los valores promedio

diarios de precipitacin, escurrimiento y carga de sedimentos en
la Unidad de Escurrimiento H1 (con obras CONSA).

ene mar may jul

sep nov

4.0

pp (mm)

3.0

4
2.0
6
1.0

8
10

escurrimiento (m3) y carga de

sedimentos (kg)

0.0
Sed (kg)

Q (m3)

Tabla N3. Valores observados de escurrimiento y emisin de

sedimentos, as como estimados de precipitacin en las unidades
de estudio aforadas en los vertedores.

H1(con obras
CONSA)
H2 (sin obras
CONSA)

A
(km2)

pp
(mm)

SEDIMENTOS
Q anual
(m3)

ton/ao

ton/ha/ao

0.10

777.1

50.04

18.02

1.74

0.18

725.2

28.27

187.13

9.92

Discusin y Conclusiones
Los datos recopilados hasta ahora muestran que la unidad de
escurrimiento del aforador H2 sin obras de conservacin,
teniendo casi el doble de superficie que la unidad de
escurrimiento del aforador H1, muestra casi diez veces menos
escurrimiento y casi tres veces ms emisin de sedimentos que
su unidad apareada adyacente.
Los valores anuales y diarios de emisin de sedimentos en la
unidad sujeta a obras CONSA (H1), destacan por ser bajos, lo
que refleja el xito -al menos temporal- de las obras realizadas
en la zona. En contraste, el vertedor H2 se consider como
aforo de ejemplo por su relativamente buena condicin de
cobertura de uso de suelo y vegetacin y carece de obras
CONSA, pero muestra una mayor retencin de agua y por otro
lado, muestra sobre todo un elevado valor promedio de
emisin de sedimentos al lago.
La distribucin espacial de los factores de la EUPS y los
valores resultantes de erosin en ambas unidades de
escurrimiento no muestran una diferencia notable y si en
cambio para el caso de unidad del aforador H2 el valor es
inferior a lo observado en cuanto a transporte de sedimentos.

pp(mm)

Agradecimientos
Ilustracin N2. Fluctuacin mensual de los valores promedio
diarios de precipitacin, escurrimiento y carga de sedimentos en
la Unidad de Escurrimiento H2 (sin obras CONSA).

Corresponde un valor anual promedio de 0.05 kg de emisin

diaria de sedimentos para el vertedor H1 y 0.11 kg para el
vertedor H2. En contraste, el vertedor H2 que carece de obras
CONSA y muestra valores ms elevados de emisin de
sedimentos al lago.
Por otro lado, siguiendo el clculo distribuido de la erosin y
algebra de mapas generados para este fin conforme al modelo
de la EUPS se tiene los valores 2.63 y 2.69 ton/ha/ao
respectivamente para las unidades de escurrimiento aforadas
por los vertedores H1 y H2.
A partir del anlisis anterior, finalmente, en la tabla 3 se
presentan los resultados de este estudio mediante un
comparativo de valores de precipitacin, escurrimiento y carga
de sedimentos en las dos microcuencas para el ao 2011.

Deseamos expresar nuestro agradecimiento a la Fundacin

Gonzalo Ro Arronte por el apoyo econmico y tcnico para la
realizacin de este proyecto, as como a la Comisin Forestal
del Estado de Michoacn, y a los bilogos lvaro San Agustn
y Juan Carlos Alvarez por su valiosa participacin.

1er CONGRESO IBEROAMERICANO SOBRE SEDIMENTOS Y ECOLOGA

QUERTARO, QUERTARO MXICO, 21-24 JULIO 2015
Referencias:
Amador G.A. 2000. Simulacin dinmica del impacto de
la agricultura en la cuenca de Ptzcuaro, Mich.
Tesis Maestra. UMSNH. 96 pp.
Brooks, K.N., Folliott, P.F. Gregersen, H.M and J.L.
Thames. 1991.
Hydrology and the
management of watersheds. Iowa State
University Press. Ames, Iowa. 392 pp.
Campos- Aranda, D.F. 2007. Estimacin y
aprovechamiento del escurrimiento. UASL.
Mxico 440 pp.
Campos, A.D.F. 2010. Introduccin a la Hidrologa
Urbana. UASL. Mxico. 220 pp.
Carabias, J, R. Landa, J. Collado y P. Martnez. 2005.
Agua, Medio Ambiente y Sociedad. UNAM.
Mxico. 219 pp.
Hudson N.H. 1997. Medicin sobre el Terreno de la
Erosin del Suelo y de la Escorrenta.
(Boletn de Suelos de la FAO - 68). FAO Organizacin de las Naciones Unidas para la
Agricultura y la Alimentacin. Roma. Versin
html:
http://www.fao.org/docrep/T0848S/T0848S00
.
IMTA 2002. Estimacin del impacto de las descargas
del dren Zurumtaro en el Lago de Ptzcuaro
y alternativas de tratamiento. Informe
Tcnico. 215 pp.
Instituto Nacional de Estadstica y Geografa (INEGI).
Modelo Digital de Terreno. E14A22.
Disponible
en:
http://www.inegi.org.mx/geo/contenidos/datosr
elieve/continental/default.aspx
Maza A.J.A. 2005. Velocidad media de inicio de inicio de
erosin en pilas, velocidad media critica de
arrastre y velocidad media de equilibrio de
erosin. In: Rivera et al. (eds.) La medicin de
sedimentos en Mxico. IMTA-UJAT. Mxico.
319 pp.
Monsalve, S. G. 1999. Hidrologa en la Ingeniera. 2.
Ed. Alfaomega. Bogot 383 pp.
Quintanilla C.A., Rivera T.F., Barajas F.J., Sanchez
R.P.A y Soto C.G. 2005. Procesamiento de
datos de sedimentos. In: Rivera et al. (eds.) La
medicin de sedimentos en Mxico. IMTAUJAT. Mxico. 319 pp.
Secretara de Agricultura y Recursos HidrulicosColegio de Posgraduados. SARH-CP. 1991.
Manual de prediccin de prdida de suelo por
erosin. 3. Ed. Direccin General de
Conservacin del Suelo y Agua. Colegio de
Postgraduados. Mxico, D.F.
San Agustn, J. A. 2012. Estimacin de la carga de
sedimentos en suspensin del dren
Tzurumutaro al lago de Ptzcuaro,
Michoacn, Mx. Tesis de Licenciatura.
Facultad de Biologa-UMSNH. Morelia,
Michoacn. 91 pp.
Tapia, V. L.M., Larios, G.A., Oropeza, M.J.L., Velsquez
V.M y A. Amador G. 2007. Simulacin de

lluvia para parametrizar el modelo WEPP en

dos sistemas de manejo agrcola de ladera. In
Snchez-Brito C. et al. (ed.) Avances de
Investigacin en Agricultura Sostenible IV.
Bases Metodolgicas para el Manejo Integral
de Cuencas Hidrolgicas. Libro Tcnico Nm.
7 INIFAP-CIRPAC. Uruapan, Mxico.566 p.
Tapia, V. L.M., Tiscareo L. M. Amador G. A.,
Fernndez, O. D. y Bez, A.D. 2001.
Simulacin de prdida de suelo en la cuenca
del lago de Ptzcuaro a partir de sistemas de
manejo agrcola con percepcin remota.
Revista de Ingeniera Hidrulica en Mxico.
XVI (4).
Valtierra, J.G. 2007. Desarrollo de una herramienta
computacional
para
el
diagnstico
hidrolgico de cuencas. Tesis Maestra. UAQ.
148pp
Wanielista, M., Kersten, R. and R. Eaglin. 1997.
Hydrology Water Quantity and Quality
Control. John Wiley & Sons. 2nd ed. USA.
Ward. A.D. and Trimble S.W. 2004. Environmental
Hydrology. Lewis Pub. 475 pp.