Aqua-LAC - Vol. 6 - Nº 1 - Mar. 2014. pp.

38 - 49
PRODUCCIÓN DE SEDIMENTOS EN UNA CUENCA DE SIERRAS PAMPEANAS,
CÓRDOBA, ARGENTINA: ESTIMACIÓN PARA DISTINTOS ESCENARIOS

SEDIMENT YIELD IN A PAMPEAN RANGE BASIN (CORDOBA, ARGENTINA):

ESTIMATIONS FOR DIFFERENT SCENARIOS

M. Jimena Andreazzini1-2, Susana Degiovanni1, Pablo Spalletti3 y Martín Irigoyen3

Resumen
La producción de sedimentos (G) en una cuenca, sintetiza la compleja interrelación entre aspectos geológico-geomor-
fológicos, climáticos, bióticos y de uso del territorio, y es una de las principales variables que definen la morfodinámica de
los cursos colectores. Ante la complejidad de efectuar mediciones directas, G puede estimarse mediante técnicas de pre-
dicción, como la ecuación de Gavrilovic, que relaciona coeficientes de producción (W) y de retención (R) de sedimentos en
una cuenca. Esta metodología fue aplicada en la cuenca del río Las Cañitas (Sierras Pampeanas, Córdoba) analizando la
influencia de las variables de control y considerando hipotéticos cambios en el uso del suelo y condiciones climáticas. Las
Sierras Pampeanas de Córdoba conforman bloques submeridianos y asimétricos, de basamento cristalino precámbrico-
paleozoico inferior, cuyas alturas varían desde 2600 a 650 m.s.n.m., descendiendo hacia el Este y Sur. Particularmente,
la cuenca de estudio se desarrolla sobre la vertiente oriental de la Sierra de Comechingones (Sierras Grandes) y sector
pedemontano. Se definieron ocho ambientes geomorfológicos donde el relieve está fuertemente condicionado por las
características litológicas y topográficas, los cuales fueron además caracterizados según las variables que intervienen en
la metodología aplicada, tales como el tipo de cobertura vegetal y el uso del suelo, entre otros. El valor de producción de
sedimentos obtenido fue de 21,43 tn/ha/año. El valor promedio de R fue 0,15 indicando que del total de la producción de
sedimentos en la cuenca, el 15% pasa a transporte y el resto quedaría retenido dentro de la misma, en las distintas sub-
cuencas. Los aspectos litológico-geomorfológicos mostraron gran influencia en los valores de G obtenidos, especialmente
sobre el parámetro W. El parámetro R, dependiente de la morfometría, mostró menor variabilidad. Los incrementos en G
asociados a escenarios post-incendios forestales mostraron que la cobertura vegetal es una variable de alta sensibilidad,
y que su pérdida incide diferencialmente según el contexto litológico-geomorfológico afectado. Considerando los esce-
narios de incrementos en precipitación y temperatura proyectados para el centro-sur de Córdoba para los años 2020 y
2050, G aumentó aproximadamente un 9% y un 20%, respectivamente, en relación a las condiciones actuales.
Palabras clave: producción de sedimentos, Sierras Pampeanas, incendios forestales, cambio climático, río Las Cañitas

Abstract
In a drainage basin, the sediment yield (G) synthesizes the complex relationships between geological, geomorphological,
climatic, biotic and land uses features, being one of key variables that define the dynamics of the main stream. Consider-
ing the complexity for direct measures, G can be estimated using prediction methods such as Gavrilovic`s equation, which
link production (W) and retention (R) sediment coefficients. This procedure was applied in Las Cañitas river basin (Com-
echingones Mountains). The influence of control variables was analyzed and hypothetical land uses and climatic condi-
tions changes were taken into account. The Córdoba Pampean Mountains are formed by lower Precambrian- Paleozoic
bedrock blocks with an asymmetrical, sub meridian arrangement and heights descending from 2600-650 m.a.s.l. towards
the East and South. The studied basin is developed on the eastern slope of Comechingones Mountains and its piedmont
area. Eight geomorphological environments characterized by a relief strongly controlled by lithological and topographic
features were defined. Each unit was characterized considering the variables involved in Gavrilovic methodology such as
vegetation cover and land use, among others. The obtained sediment yield value was 21.43 tons/ha/yr. The R average
value was 0.15, indicating that 15% of the sediment yield is transported while the rest would be retained in the different
sub-basins. Lithological-geomorphological features show great influence on the obtained G values, especially on W. The R
parameter, dependent from basin morphometry, showed less variability. G values increments, associated with post-wildfire
scenarios, showed that vegetation cover is a highly sensitive variable whose loss influences G in different ways according
to the lithological-geomorphological setting. Considering the projected climatic scenarios (2020 and 2050) for the South-
Central Córdoba region, with higher temperature and precipitation than present conditions, G increased approximately 9%
and 20%, respectively.
Keywords: yield sediments, Pampean Ranges, wildfire, landuse, climatic change, Las Cañitas river

1
Departamento de Geología, Facultad de Ciencias Exactas, Físico-Químicas y Naturales, Universidad Nacional de Río Cuarto.
Ruta 36 – Km 601, Río Cuarto, Córdoba, Argentina. CP: X5804BYA. Tel: 0358-4676229. E-mail: mandreazzini@exa.unrc.edu.ar
2
CONICET. Consejo Nacional de Investigaciones Científicas y Técnicas
3
Laboratorio de Hidráulica, Instituto Nacional del Agua (INA), Ezeiza.

Recibido: 28/08/2013
Aceptado: 25/03/2014

38 Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014

 Producción de sedimentos en una cuenca de sierras pampeanas, Córdoba, Argentina: estimación para distintos escenarios
INTRODUCCIÓN geológicos, climáticos y de uso del territotio, con
resultados en general aceptables. En numerosas
La producción de sedimentos sintetiza la compleja
cuencas de montaña ha podido ser validada a partir
interrelación entre los aspectos geológico-geomor-
del contraste de los resultados con mediciones de
fológicos, climáticos, bióticos y de uso del territorio
caudales sólidos en suspensión y la comparación de
de una cuenca y a su vez es una de las variables
batimetrías en embalses. Originalmente, este méto-
principales que definen la morfodinámica de los cur-
do fue desarrollado para ser aplicado en cuencas de
sos colectores. Su conocimiento adquiere alto valor,
montaña en el Sur y Sudeste de Yugoslavia, y luego
en tanto permite efectuar diagnósticos y pronósticos
fue utilizado en Suiza, Italia (Bazzoffi, 1985), Croacia
del funcionamiento de los sistemas fluviales, tanto
(Globevnik et al. 2003), Grecia (Emmanouloudis et
en condiciones naturales como intervenidas; sin em-
al. 2003), entre otros. Más recientemente diferentes
bargo, su cuantificación es compleja.
autores utilizan esta metodología en conjunto con
Para definir propuestas de manejo de cuencas
la Ecuación Universal de Pérdida de Suelo (USLE)
fluviales, ya sean de ordenamiento territorial o
(Panagos et al. 2014), o para contrastar los resul-
intervenciones directas en los cursos, resulta
tados alcanzados con otros modelos, tales como el
imprescindible conocer no sólo el volumen sino el
PSIAC (Pacific Southwest Inter Agency Committee)
comportamiento del material sólido generado en una
(Vezzoli et al. 2013), obteniendo un buen nivel de
cuenca, que es transportado luego por el sistema
aproximación. Tosic et al. (2012) emplea la ecua-
fluvial.
ción de Gavrilovic para comparar producción de
La cantidad de sedimentos entregados por una sedimentos en relación a cambios importantes en
cuenca hidrológica en la unidad de tiempo se define el uso del territorio, como los ocurridos en la situa-
como producción de sedimentos. No todo el mate- ción post-guerra en Bosnia y Herzegovina. En otros
rial movilizado por procesos erosivos en una cuenca estudios, como los realizados por Abadi & Ahmadi
alcanza la red de drenaje, por lo que se considera (2011) en la cuenca Kasilian (Irán), se observa que
como producción de sedimentos al volumen que se la metodología de Gavrilovic sobrestima los datos
transporta hacia los colectores principales. En cuen- medidos, aunque los autores sugieren que es po-
cas montañosas dicha producción se puede asociar sible mejorar los resultados utilizando los Sistemas
fundamentalmente a dos fuentes, una relacionada de Información Geográfica para discretizar los pará-
a la erosión superficial generalizada del sustrato, y metros intervinientes. En Argentina se ha empleado
otra a los procesos de remoción en masa en sitios satisfactoriamente en cuencas del Noroeste (Brea
puntuales. La cantidad de sedimentos generados et al. 1999; Olmos et al. 2008; COBINABE, 2010;
en una determinada cuenca hídrica y, consecuente- Brea y Spalletti, 2010) y en la Patagonia (Spalletti
mente, el grado de erosión de la misma, pueden eva- y Brea, 2006). Beltramone et al. (2005) presentan
luarse en forma indirecta a través de relevamientos estimaciones realizadas por este método para cuen-
topo-batimétricos de los vasos de lagos y embalses cas serranas de Córdoba. Diversos trabajos han uti-
o mediante la utilización de mediciones sistemáticas lizado esta metodología en combinación con técni-
de caudales sólidos en los cursos de agua (COBI- cas de SIG, demostrando un gran potencial para la
NABE, 2010). Sin embargo, es un problema común estimación de la erosión de suelos en las cuencas
de muchas cuencas la inexistencia de estaciones de (Globevnik et al. 2003; Amini et al. 2010; Olmos et
aforos sólidos, de caracterizaciones granulométri- al. 2008; COBINABE, 2010; Brea y Spalletti, 2010;
cas del material transportado en suspensión y por Tangestani, 2006, entre otros).
arrastre de fondo, y de datos de sedimentación en
Específicamente, para los cursos que drenan las
espejos de agua. Resulta necesario entonces utili-
Sierras Pampeanas de Córdoba, en especial su
zar técnicas de predicción para determinar órdenes
extremo sur, aún son pocos los trabajos que se re-
de magnitud de la producción de sedimentos. En
fieren a la producción y/o transporte de sedimentos
general, las formulaciones analizan la interacción de
y sus principales características (textura, composi-
factores tales como clima, litología, cobertura vege-
ción, etc.) entre los que pueden citarse Corral et al.
tal, uso del suelo y topografía.
(2009), Ulla (2008) y Degiovanni et al. (2010), entre
A pesar de la gran cantidad de metodologías exis- otros. Considerando los aspectos morfodinámicos e
tentes, la valoración cuantitativa precisa de la pro- intervenciones realizadas y previstas en la red de
ducción de sedimentos, es imposible de asegurar drenaje de esta región, el conocimiento de todas es-
debido a la complejidad del proceso erosivo, a la tas variables y sus proyecciones futuras, resulta de
cantidad de factores que influyen en el mismo y a la primordial importancia para mejorar la gestión de los
escala de los fenómenos actuantes. Sin embargo, a recursos hídricos.
pesar de las dificultades, pueden ajustarse ciertas
La mayoría de los cursos del sur de Córdoba nacen
metodologías con una precisión aceptable, que per-
en la Sierra de Comechingones (Sierras Pampea-
mita cumplir los objetivos perseguidos (UNESCO,
nas) y se integran a la cuenca del Plata, ya sea natu-
2010).
ral o artificialmente (Degiovanni et al., 2005, Blarasin
Entre ellas, la metodología de Gavrilovic (1972, et al., 2005). Especialmente en la última centuria, en
1988) ha sido empleada en diferentes contextos respuesta a cambios climáticos, de uso del territorio

Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014 39

M. Jimena Andreazzini, Susana Degiovanni, Pablo Spalletti y Martín Irigoyen
e intervenciones ingenieriles, la red de drenaje en su La cuenca del río Las Cañitas forma parte de la
conjunto se encuentra desajustada, dominando los cuenca alta del río Cuarto (provincia de Córdoba,
procesos de incisión y erosión retrocedente, en las Argentina), y drena una superficie de 267 Km2 en el
cuencas medias y bajas, con una alta tasa de trans- faldeo oriental de la Sierra de Comechingones (Fig.
porte de sedimentos hacia las áreas de descarga, 1), entre los 1800 y los 650 m.s.n.m. Esta vertiente
lo que promueve nuevas intervenciones (Blarasin et está fuertemente disectada por la red de drenaje,
al., 1993; Cantero et al., 1998; Degiovanni, 2008; que preserva relictos de superficies de erosión en
Degiovanni et al., 2009; 2014). los sectores cumbrales (“pampas de altura”). En
En las cuencas altas se registran alteraciones vin- el área serrana dominan rocas metamórficas poco
culadas a sobrepastoreo, obras hidráulicas y, espe- alteradas y, subordinadamente, rocas miloníticas y
cialmente, incendios forestales. Esta problemática, graníticas. En el pedemonte, valles serranos y pam-
muy común y recurrente en las sierras de Córdoba, pas de altura, se presentan sedimentos cuaternarios
tiene alto impacto en los caudales líquidos y sólidos aluviales, coluviales y eólicos.
de los cursos de agua. La mayoría de los años se Desde el punto de vista geomorfológico se recono-
registran numerosos focos de incendio (150-200), cieron ocho ambientes principales (Fig. 1), los cua-
especialmente entre los meses de julio y noviembre, les son caracterizados considerando además aque-
que involucran importantes áreas de bosques y pas- llas variables que intervienen en la metodología
tizales naturales. En los años 2008, 2010 y 2013, se aplicada, como la vegetación y el suelo, entre otras.
afectaron 61.541,65 has, 10.305 has (Secretaría de Las pampas de altura (I-A) presentan bajas pendien-
Ambiente y Desarrollo Sustentable, 2009 y 2011) y tes, afloramientos rocosos parcialmente cubiertos
más de 50.000 has, respectivamente. por loess, suelos de alto desarrollo, erosión hídrica
Por otra parte, considerando los cambios en las localizada, canales de lecho rocoso poco profundos,
variables climáticas (precipitación y temperatura) escasos materiales aluviales y alta cobertura de
proyectados para la región centro-Sur de Córdoba pastizales. La superficie metamórfica central (I-B)
(AIACC, 2006) también es esperable que surjan mo- exhibe un 40 a 60% de superficie rocosa, pendien-
dificaciones en la producción de sedimentos, por lo tes moderadas-altas, cursos muy incididos de lecho
que resulta de interés evaluar también estos esce- rocoso y aluvión subordinado, procesos gravitato-
narios hipotéticos. rios e hídricos muy localizados asociados a depósi-
tos coluviales, parcialmente edafizados y cubiertos
Así, el objetivo de este trabajo fue estimar la tasa de
de pastizales. El ambiente granítico (I-C) posee alto
producción actual de sedimentos de una subcuenca
porcentaje de afloramientos de roca alterada (areni-
representativa de la problemática descripta, analizar
zación) y erosión activa, cursos con mayor material
la influencia de las variables de control y efectuar
aluvional, y cubierta de pastizales discontinua. El
estimaciones de producción considerando cambios
sector serrano oriental y bloque pedemontano (I-D
en el uso del suelo y en las condiciones climáticas.
y E) presenta menor pendiente y áreas de aflora-
Se seleccionó la subcuenca del río Las Cañitas, que
mientos, mayor agradación en valles, mayor cober-
integra el sistema del río Cuarto, uno de los más in-
tura arbórea-arbustiva y pastizales, mejor desarrollo
tervenidos de la región, especialmente en su cuenca
edáfico, y cárcavas locales en ejes de valles. La
media debido a que presenta extracción de áridos,
planicie de agradación pedemontana (IIA, B, C) esta
rectificaciones y obras transversales, y en su cuenca
incidida por los ríos principales, que desarrollaron
baja por la existencia de canalizaciones.
cursos meandriformes hasta su semiconfinamiento
en secuencias cohesivas. Los procesos de erosión
Cuenca del río las cañitas: características gene- hídrica son moderados a leves.
rales El clima de la región es de tipo Templado Subhú-
La cuenca de estudio se desarrolla sobre la vertien- medo, disminuyendo las precipitaciones medias
te oriental de la Sierra de Comechingones (Sierras anuales (PMA) y la temperatura media anual (TMA)
Grandes) y sector pedemontano adyacente. Estas desde la zona pedemontana (Estación Las Tapias,
serranías forman parte de las Sierras Pampeanas de 728 m.s.n.m., con 910.9 mm y 14 ºC, Fig. 1) hacia
Córdoba, constituidas por bloques submeridianos de las pampas de altura (Estación La Cumbre, 1500
basamento cristalino precámbrico-paleozoico infe- m.s.n.m., con 760.8 mm y 11.9 ºC, Fig. 1). Las llu-
rior, que conforman cordones montañosos de perfil vias son estacionales, concentrándose el 80% en el
transversal asimétrico. Su borde occidental es más periodo primavera-verano.
abrupto y está representado por la escarpa de falla Respecto al uso del suelo, el sector serrano se des-
principal, mientras que el oriental constituye la pen- tina a ganadería extensiva, con mayor sobrecarga
diente estructural de menor pendiente. Sus alturas de ganado vacuno en las pampas de altura. En el
varían desde 2600 a 650 m.s.n.m., descendiendo sector pedemontano predomina la agricultura.
hacia el Este y Sur.

40 Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014

 Producción de sedimentos en una cuenca de sierras pampeanas, Córdoba, Argentina: estimación para distintos escenarios

Figura1. Ubicación de la cuenca Las Cañitas y mapa geomorfológico.

Se sitúan además las estaciones meteorológicas La Cumbre y Las Tapias.

METODOLOGÍA APLICADA el autor de la metodología y representan respectiva-

mente, el grado de protección del suelo dado por la
En este trabajo se aplicó el método de Gavrilovic vegetación y la intervención antrópica, el grado de
(1972, 1988), que considera el volumen de sedimen- resistencia a la erosión del suelo considerando sus
tos “G” (en m3/año) producido por erosión y transpor- características geolitológicas, y el estado erosivo de
tado a la sección de cierre de una cuenca, como el la cuenca. I es el gradiente de la pendiente superfi-
producto de la producción media anual de sedimento cial (en %).
por erosión superficial “W” y el coeficiente de reten-
ción de sedimentos “R”. Por su parte, el coeficiente R indica la relación entre
el volumen de sedimento que efectivamente pasa
La expresión para determinar el volumen promedio por la sección de cierre de la cuenca y el volumen
anual erosionado de sedimento es: total de material producido por erosión superficial. La
ecuación para su cálculo según Zemlijc (1971) es la
W = T · h · π · Z3/2 · F (m3/año) (1) siguiente:

Donde: R = [O · D1/2 · (L+Li)]/(L+10) · F (4)

T = coeficiente de temperatura, que se obtiene de:

T = [(t/10) + 0.1]1/2 (2) Donde O es el perímetro de la cuenca (km), D la

diferencia de nivel media en la cuenca (cota media
menos cota mínima, en km), Li la longitud total de los
Siendo t la temperatura promedio anual (°C), h la afluentes fluviales laterales (km), L la longitud de la
precipitación media anual (mm/año) y F la superficie cuenca por el talweg del cauce principal (km) y F la
de la cuenca (km²). Z es el coeficiente de erosión, superficie de la cuenca (km²).
cuya expresión es: Toda la información fue trabajada en ambiente SIG,
obteniendo el valor de “G” en 47 subcuencas.
Z = X · Y · (φ + I)1/2 (3) El valor del coeficiente X se asignó a partir de un
mapa de cobertura vegetal y uso del suelo, elabora-
Los valores de los coeficientes X (entre 0.1 y 1), Y do desde una imagen Landsat 5 de fecha 31/12/08,
(entre 0.2 y 2) y φ (entre 0.1 y 1) son propuestos por y tomando como base el trabajo de Suárez (1993) e

Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014 41

M. Jimena Andreazzini, Susana Degiovanni, Pablo Spalletti y Martín Irigoyen
imágenes Google Earth. Se realizó una clasificación Centro Hadley; y escenarios de emisión A2 (mayor
supervisada por el método de máxima verosimilitud, desarrollo económico) y B2 (mayor sustentabilidad
obteniendo 9 clases. ambiental). Para el presente trabajo se tomaron so-
Para definir la distribución espacial de los paráme- lamente los resultados obtenidos con el primer MCG
tros Y y φ, se utilizó la información litológica y geo- (denominados escenarios E-A2 y E-B2), ya que fue-
morfológica, la de las características edáficas y la de ron los que proyectaron los máximos incrementos
los procesos activos observados en la cuenca, ba- en precipitación, y sólo se evaluaron cambios en la
sadas en datos de campo, trabajos antecedentes e producción de sedimentos modificando las variables
información cartográfica disponible. climáticas en la metodología de Gavrilovic.
El mapa de pendientes (en %) fue obtenido a par- Para el caso del escenario E-A2, la PMA de la serie
tir de un modelo de elevación digital de la banda X Las Tapias-Los Chañares y de la Estación La Cum-
del Shuttle Radar Topographic Mission (25 m de re- bre fue incrementada en un 6.57% y la TMA en 0.07
solución), disponible gratuitamente en EOWEB-NG ºC para 2020, mientras que para 2050, la PMA se
(2013). incrementó en un 16.94% y la TMA en 0.64 ºC.
La distribución de las lluvias y temperatura en fun- En el escenario E-B2, la PMA de las series utiliza-
ción de la cota se obtuvo por regresión, utilizando das fue incrementada en un 6.10% y la TMA en 0.19
los datos de las estaciones Las Tapias y La Cumbre ºC para el escenario 2020. Para el 2050, la TMA se
(Fig. 1). incrementó en un 13.14% y la TMA en 0.55 ºC. Si
bien el registro de la Estación La Cumbre es de corta
Desde el modelo de elevación digital y empleando
duración, se lo utilizó ya que son datos medidos en la
herramientas SIG, se obtuvo la red de drenaje y un
cuenca, de modo que los resultados obtenidos sur-
total de 47 subcuencas, en general de orden 3. A
gen de considerar estos datos como representativos
partir de ello se calcularon las áreas y perímetros de
de un periodo de mayor longitud (el evaluado como
las subcuencas y la longitud de los cursos principa-
escenario de base).
les y tributarios.
Se calculó el valor de W para las distintas subcuen-
RESULTADOS
cas, y finalmente se obtuvo el volumen de sedimento
G (m3/año) para cada una de ellas y para la cuenca Producción de sedimentos para las condiciones
completa. actuales
Luego de obtenida la producción de sedimentos para La Figura 2 presenta el mapa de clases de vegeta-
las condiciones climáticas y de uso del suelo actual, ción (incluyendo las áreas correspondientes a aflora-
se realizaron cálculos para compararla con la pro- mientos rocosos y cauces fluviales) y la Figura 3 el
ducción en escenarios hipotéticos de: 1) situación mapa del coeficiente X para el área de estudio. Se
post-incendios forestales, y 2) incremento en la pre- asignaron los menores valores de X a las unidades
cipitación y temperatura media anual. Bosque (que es de tipo alto y abierto en este caso,
En el primer caso y teniendo en cuenta la hetero- con vegetación arbustiva y herbácea en la base),
geneidad espacial de las variables que controlan los Pradera (pastizales con más de 70% de cobertu-
procesos erosivos, se consideró conveniente evaluar ra), Mallines (cobertura del 90% de herbáceas, do-
las modificaciones en la producción total de sedi- minando pastizales), y Estepa Arbustiva (cobertura
mentos, si el incendio afectara a diferentes unidades herbácea con monte bajo). Valores intermedios se
geomorfológicas. Se asumió un incendio hipotético asignaron a la unidad Estepa Herbácea, que posee
de 550 has en los ambientes de relieves graníticos (I- un 60-70 % de cobertura herbácea y 30-40 % de
C), la unidad metamórfica central (I-B), y las pampas afloramientos rocosos, y a las Áreas Cultivadas, don-
de altura (I-A), siendo el área quemada igual dentro de la cobertura dada por los cultivos genera menos
de cada ambiente. Fueron modificados únicamente protección que en los casos anteriores. Los valores
los valores del coeficiente X del método de Gavrilovic mayores se indicaron para los sectores donde domi-
en los tres casos analizados, y dentro de esas áreas nan los afloramientos rocosos (60% o más) y para
no se incluyeron los afloramientos rocosos. los cauces fluviales, siendo en ambos casos la co-
En el segundo caso se empleó información sobre bertura vegetal mínima o ausente.
escenarios de cambio climático proyectados para Para el coeficiente Y (Fig. 4), se asignaron los valo-
la región centro-sur de la provincia de Córdoba, res más altos a los suelos arenoso-limosos con bajo
disponibles en AIACC (2006). En dicho estudio se desarrollo edáfico del sector extraserrano, y a los
construyeron escenarios climáticos base a partir de desarrollados en la faja fluvial del río Las Cañitas. A
datos del periodo 1961-1990 de las variables tem- los suelos de las pampas de altura, si bien son loés-
peratura y precipitación, y escenarios de cambio sicos, se les otorgó un valor menor dado su mayor
climático para el 2020 y 2050. Estos autores utiliza- contenido en arcillas y materia orgánica, su fuerte
ron las salidas de tres Modelos Climáticos Globales estructura y su mayor espesor. Valores intermedios
(MCG): EH4TR98 del European Center for Medium se seleccionaron para los suelos de escaso desarro-
Range Weather Forecast, GFDLTR90 del Geophy- llo, sobre materiales heterogéneos, presentes en la-
sical Fluids Dynamics Laboratory, y HAD3TR00 del deras de variada pendiente dentro del ambiente gra-

42 Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014

 Producción de sedimentos en una cuenca de sierras pampeanas, Córdoba, Argentina: estimación para distintos escenarios
nítico y metamórfico. Finalmente, los más resistentes mientos rocosos (con una mayor resistencia en las
se asociaron a los suelos de los mallines (resistencia metamorfitas, dado que la roca se encuentra con
debido a la saturación), sedimentos aluviales de gra- muy bajo grado de alteración).
nulometría muy gruesa de los cauces y a los aflora-

Figura 2. Mapa con las clases de vegetación definidas para la cuenca del río Las Cañitas.

Figura 3. Mapa del coeficiente X de la metodología de Gavrilovic para la cuenca Las Cañitas.

Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014 43

M. Jimena Andreazzini, Susana Degiovanni, Pablo Spalletti y Martín Irigoyen
En el mapa para el coeficiente φ (Fig. 5) las unida- comprobarse comparando su extensión y desarrollo
des con menor presencia de rasgos erosivos son en fotografías aéreas del año 1970 e imágenes sate-
las vinculadas al ambiente metamórfico, donde los litales actuales. En el ambiente granítico, los rasgos
afloramientos se encuentran prácticamente inaltera- erosivos son más importantes, y es posible observar
dos, y en las laderas de los valles se observa mí- el grado de alteración que presenta la roca expuesta,
nima erosión, asociada a pequeñas cárcavas que que lleva a la arenización de la misma. Se observan
en general se encuentran estabilizadas. Desde las rasgos de erosión laminar y rills, además de cárca-
laderas de mayor pendiente llegan al canal aportes vas en los materiales sueltos que conforman las la-
por procesos gravitatorios muy localizados. Valores deras.
bajos también fueron asignados a los sectores de Por otra parte, el mapa de pendientes para la cuen-
Bosque y Estepa Arbustiva, asociados al sector más ca se muestra en la Figura 6a. La pampa de altura,
oriental de la sierra, donde localmente se desarrollan el sector serrano oriental y el extraserrano muestran
cárcavas en los ejes y laderas de los valles con relle- los menores valores, mientras que el sector serrano
no sedimentario. En el sector extraserrano, también central los más elevados, en especial la unidad me-
los procesos de erosión hídrica en surcos y cárcavas tamórfica, con abruptas laderas asociadas al marca-
están muy localizados, especialmente en los bordes do proceso de incisión fluvial, controlado por estruc-
de terrazas fluviales, por lo que el coeficiente asig- turas de dirección NO-SE y ENE-OSO.
nado fue intermedio a bajo. Entre las unidades con
Por su parte, la Figura 6b muestra el modelo de
mayor evidencia de procesos de erosión se incluye-
elevación del terreno utilizado y las subcuencas y
ron algunos sectores en las pampas de altura, donde
cursos fluviales principales, a partir de los cuales se
se han desarrollado cárcavas de cientos de metros
calcularon los parámetros para determinar el coefi-
de longitud, aunque en la actualidad se encuentran
ciente “R”.
parcialmente vegetadas y estabilizadas, como puede

Figura 4. Mapa del coeficiente Y de la metodología de Gavrilovic para la cuenca Las Cañitas.

El valor de producción de sedimentos obtenido, su- distintas subcuencas. La Figura 7 y Tabla 1 muestran
mando la producción de las 47 subcuencas fue de la distribución espacial de la producción de sedimen-
216220.41 m3/año, lo que da un total de 21.43 tn/ha/ tos en los diferentes ambientes geomorfológicos de
año. Según esta metodología, el valor promedio de R la cuenca.
dio 0.15, indicando que del total de la producción de Como puede observarse hay una marcada variación
sedimentos en la cuenca, el 15% pasa a transporte y de la producción de sedimentos, con valores máxi-
el resto quedaría retenido dentro de la misma, en las

44 Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014

 Producción de sedimentos en una cuenca de sierras pampeanas, Córdoba, Argentina: estimación para distintos escenarios
mos en el ambiente granítico, controlado fundamen- y una cubierta de pastizales más discontinua. Este
talmente por los coeficientes Y y φ, que evidencian ambiente con sólo el 4.7 % del área total aporta el
el mayor grado de meteorización de las rocas gra- 8.4 % de G, al igual que otras unidades que triplican
níticas, la presencia de procesos erosivos activos el área.

Figura 5. Mapa del coeficiente φ de la metodología de Gavrilovic para la cuenca Las Cañitas.

Figura 6. a- Mapa de pendientes y b- Modelo de elevación digital, red de drenaje principal

y subcuencas definidas para la cuenca Las Cañitas.

Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014 45

M. Jimena Andreazzini, Susana Degiovanni, Pablo Spalletti y Martín Irigoyen
Las unidades con menor producción son las del sector finos, de moderado desarrollo y uso agrícola (II-C).
serrano oriental (I-D y E) y pampas de altura (I-A) dado Finalmente, la unidad I-B con una producción interme-
sus menores pendientes, desarrollo edáfico y cobertu- dia, aporta casi el 67% de G debido a su mayor super-
ra vegetal. Se destaca la baja producción de la unidad ficie, y es la que define el aporte medio de la cuenca
I-A, asociada a suelos con alto contenido en materia (21.4 tn/ha/año). En este caso la distribución interna
orgánica y arcillas (y saturados en los mallines), e im- de G es muy variable, y los más altos valores se aso-
portante cobertura de pastizales. La zona extraserrana cian a las laderas de moderada a alta pendiente de los
muestra una producción intermedia, donde los mayo- valles incididos, cubiertas de derrubios con suelos de
res valores se asocian a suelos limosos-arenosos muy alta infiltración y poco profundos.

Figura 7. Mapa de la tasa anual promedio de producción de sedimentos

por erosión superficial cuenca Las Cañitas (m3/año).

Tabla 1. Cuenca Las Cañitas. Producción de sedimentos en condiciones actuales

(G) y post-incendio (Gi=promedio en cada unidad x total cuenca, Git=total cuenca x
promedio en cada unidad). Se incluye el valor promedio de R de cada unidad.
Unidad R G % cambio
Área Gi Git
Geomorfológica promedio tn/ha/año - % Gi vs. G
Total cuenca 100 0.16 21.4-100 --- --- ---
I-A 13.2 0.16 9.8-6 15 22 53.1
I-B 57.3 0.18 25-66.9 26 22.4 6
IC 4.7 0.2 38.9-8.4 66 22.8 69.9
I-D+I-E 13.3 0.1 14.3-8.8 19 22.2 32.9
II-A, ByC 11.7 0.1 18-9.8 --- --- ---

PRODUCCIÓN DE SEDIMENTOS PARA ESCENA- los cambios en la producción total de la cuenca (Git),
RIOS HIPOTÉTICOS al considerar incendios en diferentes unidades, sur-
gen incrementos mayores, especialmente en I-A y
Incendios forestales
I-C, demostrando la alta sensibilidad de estos am-
En la Tabla 1 se presentan los valores de producción bientes a la pérdida de vegetación.
de sedimentos para una situación post-incendios fo- En el caso de la pampa de altura, con un 13% del
restales en la cuenca. Si bien no son significativos área quemada, la producción se incrementa en más

46 Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014

 Producción de sedimentos en una cuenca de sierras pampeanas, Córdoba, Argentina: estimación para distintos escenarios
del 50%. La pérdida de los pastizales de alta co- alterar la protección del suelo dada por la cobertura
bertura deja expuestos suelos desarrollados sobre vegetal, ya que los otros coeficientes de la metodolo-
materiales loessoides, los cuales a pesar de la baja gía poseen valores intermedios a altos.
pendiente que posee esta unidad, muestran alta sen- Como ejemplo, en el año 2011 y 2013 ocurrieron in-
sibilidad a los procesos erosivos. Si bien estos suelos cendios en la cuenca Las Cañitas (Fig. 8a), y el con-
poseen buena estructura y alto contenido en materia tenido de sólidos en suspensión de los arroyos se in-
orgánica, debe tenerse en cuenta que los incendios crementó durante el periodo de lluvias siguiente (Fig.
también pueden afectar las propiedades físicas del 8b). El caudal sólido medido en el río Las Tapa, a la
suelo (agregación-estructura), convirtiéndose en salida del área serrana, para caudales del orden de
muchas oportunidades en suelos hidrofóbicos, con 0,7-0,8 m3/seg varió de 26,9 a 92 gr/seg en situacio-
lo que se potencia la tasa de erosión. Para el caso nes normales y con posterioridad a los incendios del
de la unidad granítica, con un 44% de área afecta- 2011, respectivamente. Para el caso del río Las Ta-
da, el valor de G aumenta un 70%, demostrando que pias, la concentración de sólidos en suspensión varió
esta unidad posee un alto potencial de producción al de 0-0,07 gr/L en condiciones normales a más de 5
gr/L para crecidas post-incendios 2013 (Fig. 8b).

Figura 8. a) Vista de un sector de la subcuenca del río Las Tapias afectado por los incendios del 2013.
b) Incremento de sólidos en suspensión en el río Las Tapias a la salida del sector serrano,
en una crecida posterior a los incendios.

Cambio climático Tabla 2. Producción de sedimentos (tn/ha/año)

cuenca Las Cañitas en los 4 escenarios climáticos
Estos escenarios también mostraron incrementos en
proyectados.
la producción de sedimentos (Tabla 2). Como puede
observarse, la producción para el 2020 aumenta de % cambio
manera similar en ambos escenarios, en algo más Unidad E-A2 E-A2 E-B2 E-B2 E-A2 2050
G
del 9% respecto a la condición actual. Aquí no se Geomorfológica 2020 2050 2020 2050 respecto
aG
han marcado diferencias apreciables ya que la va-
riación en la temperatura en el escenario A2 es muy I-A 9.8 10.4 11.7 10.5 11.5 19.4
pequeña. Para el 2050 los valores se incrementan
I-B 25 27.3 30.6 27.4 29.8 22.4
aproximadamente en un 20% respecto a la produc-
ción total actual de la cuenca y se observan más di- IC 38.9 42.8 46.2 42.6 46.4 18.8
ferencias entre escenarios, correspondiendo la ma-
I-D+I-E 14.3 15.8 17.7 15.8 17.2 23.8
yor producción al E-A2.
Para comparar los cambios en la producción dentro II-A, ByC 18 19.9 22.4 19.9 21.7 24.4
de las distintas unidades geomorfológicas, a modo Total cuenca 21.4 23.4 26.2 23.4 25.5
de ejemplo se muestran los porcentajes de cambio
del escenario E-A2 para 2050 respecto a las condi- % cambio de
ciones actuales (Tabla 2), observándose que las ma- la producción --- 9.3 22.4 9.4 19.2
yores diferencias son para el sector serrano oriental total
y extraserrano. Tal como lo muestra la ecuación (1),
al no modificar ninguno de los demás factores del CONCLUSIONES
método, las mayores producciones se asocian a las
mayores precipitaciones y temperaturas, en este Los resultados obtenidos indican que los aspectos
caso, en el sector oriental de la cuenca. litológico-geomorfológicos tienen una marcada in-
fluencia en la producción de sedimentos en la cuen-

Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014 47

M. Jimena Andreazzini, Susana Degiovanni, Pablo Spalletti y Martín Irigoyen
ca, especialmente sobre el parámetro W. El paráme- del rio Cuarto. Pcia. de Córdoba. Actas del XVI Con-
tro R, dependiente de la morfometría, mostró menor greso Geológico Argentino. La Plata. 581-585.
variabilidad. Además se observó que la distribución Blarasin, M., Cantú, M., Degiovanni, S., Becker, A.
areal de las lluvias y temperaturas tiene menor in- & C. Eric, 1993. Estudio geoambiental de un área
cidencia para explicar variaciones en dicha produc- de llanura loéssica sometida a alta presión de uso
ción. agrícola: el caso de General Deheza, Córdoba. XII
Con relación a los escenarios evaluados, globalmen- Congreso Geológico Argentino, Mendoza. Tomo VI:
te la producción de sedimentos se incrementa para 274-283.
los escenarios climáticos proyectados, aunque sería
Blarasin M., Degiovanni S., Cabrera A. & Villegas
conveniente analizar también las variaciones de X y
M. 2005. Aguas superficiales y subterráneas en el
φ ante tales cambios, ya que las características de la
Sur de Córdoba: Una perspectiva geoambiental. 1ª
cobertura vegetal y uso del suelo como así también
ed.- Universidad nacional de Río Cuarto. Río Cuarto,
la presencia de rasgos erosivos en la cuenca posi-
Argentina, 346 p.
blemente, se verán afectados.
Por otra parte, los incrementos en la producción de Brea, D. & P. Spalletti. 2010. Generación y transpor-
sedimentos asociados a incendios forestales mos- te de sedimentos en cuencas de ríos de montaña.
traron que en general la cobertura vegetal es una Documento Técnico Nº 22, 52-74. En Brea, J.D. &
variable de alta sensibilidad, y que su pérdida incide F. Balocchi (Eds.). Procesos de erosión - sedimenta-
diferencialmente según el contexto litológico-geo- ción en cauces y cuencas. Volumen 1. Programa Hi-
morfológico afectado. drológico Internacional de la UNESCO para América
Latina y el Caribe.
Si bien los resultados obtenidos son preliminares,
demuestran la importancia de discretizar el análisis Cantero, A., M. Cantú, J. Cisneros, J. Cantero, M.
de la producción de sedimentos en una cuenca de- Blarasin, V. Becerra, J. González, A. Degioanni, H.
terminada, para lograr una mayor representatividad Gil, J. de Prada, S. Degiovanni, C. Cholaky, M. Ville-
de las variaciones de los parámetros que controlan gas, A. Cabrera & C. Eric, 1998. Las tierras y aguas
la producción y el almacenamiento de sedimentos. del Sur de Córdoba: propuestas para un manejo sus-
tentable. Ed. UNRC. ISBN: 950-665-101-9.
Finalmente, a pesar del carácter estimativo de la
metodología empleada en el presente trabajo, los COBINABE - Comisión Binacional para el desarrollo
valores asignados a los diferentes coeficientes eva- de la Alta Cuenca del Río Bermejo y el Río Grande
luados en la metodología de Gavrilovic surgen de un de Tarija-. 2010. Generación y Transporte de Sedi-
detallado tratamiento de la información disponible, mentos en la Cuenca Binacional del Río Bermejo.
apoyado en un importante relevamiento de campo, Caracterización y análisis de los procesos intervi-
por lo cual se considera que los resultados obtenidos nientes. 230 p.
son consistentes y que este análisis podría ser ex- Corral, M., Baldissone, M., Farías, H., Rodriguez, A.
tendido a otras cuencas de características similares & F. López. 2009. Balance sedimentológico simplifi-
en Sierras Pampeanas. cado para el río Cuarto, Córdoba. Actas RIOS 2009:
IV Simposio Regional sobre Hidráulica de Ríos, Sal-
REFERENCIAS BIBLIOGRÁFICAS ta, Argentina.
Degiovanni, S. 2008. Análisis geoambiental del com-
Aalto, R., Dunne, T. & J.L. Guyot 2006. Geomorphic
portamiento de los sistemas fluviales del Sur de
controls on Andean denudation rates. The Journal of
Córdoba, en especial del Aº Achiras-del Gato, como
Geology (114): 85-99.Abadi, l.Z. & H. Ahmadi, 2011.
base para su gestión sustentable. Tesis Doctoral,
Comparison of EPM and geomorphology methods UNRC. Inédita. 597 pág.- Río Cuarto-Córdoba.
for erosion and sediment yield assessment in Ka-
silian Watershed, Mazandaran Province, Iran. DE- Degiovanni, S., N. Doffo & M. Villegas, 2009. Reju-
SERT 16: 103-109. venecimiento de la red de drenaje del Sur de Cór-
doba como consecuencia del cambio climático du-
Amini, S., Rafiei, B., Khodabakhsh, S. & M. Heyda- rante la última centuria. En: Geomorfología y Cambio
ri, 2010. Estimation of erosion and sediment yield of Climático- Capítulo 6. Sayago y Collantes (Eds.), IN-
Ekbatan Dam drainage basin with EPM, using GISS. GEMA (Universidad Nacional de Tucumán)- MAGNA
Iranian Journal of Earth Sciences (2): 173-180. Ediciones. 20 pág.
AIACC - Assessments of Impacts and Adaptations to Degiovanni, S., Andreazzini, J., Villegas, M., Eric,
Climate Change -. 2006. Vulnerability and Adapta- C. & Doffo, N. 2010. Caracterización de la carga en
tion to Climate Variability and Change: The Case of suspensión de los arroyos del sur de Córdoba. Im-
Farmers in Mexico and Argentina. Final Report Pro- plicancias morfodinámicas y ambientales, 185-192.
ject Nº LA29. 173 p. En Varni, M., Entraigas, I. & L. Vives (Eds.). Hacia
Beltramone, C., Barbeito, O., Ambrosino, S. & R. la gestión integral de los recursos hídricos en zonas
Ugarte. 2005. Determinación del área fuente y pro- de llanura. I Congreso Internacional de Hidrología de
ducción de sedimento anual para la cuenca superior Llanuras. Azul, Buenos Aires, Argentina, 1301 p.

48 Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014

 Producción de sedimentos en una cuenca de sierras pampeanas, Córdoba, Argentina: estimación para distintos escenarios
Degiovanni, S., Echevarría, K., Andreazzini, J. & N. Secretaría de Ambiente y Desarrollo Sustentable,
Doffo, 2014. Los sistemas fluviales del sur de Cór- Argentina. 2011. Estadística de incendios forestales
doba: alteraciones históricas, procesos de ajuste y 2010, 01/09/13, www.ambiente.gov.ar/archivos/web/
consecuencias ambientales. Actas CD, XIX Congre- PNEF/file/Incendios%202010/Documento%20com-
so Geológico Argentino. Córdoba. 2 p. pleto.PDF.
EOWEB- Earth Observation Center-. German Re- Suarez, S.A. 1993. Mapeo fisionómico-estructural de
mote Sensing Data Center, Alemania, 01/09/2013. la vegetación natural del Dpto. Río Cuarto. III Achi-
http://eoweb.dlr.de:8080/, e-mail: helpdesk-dfd@dlr. ras y Las Cañitas. Informe final Beca de Iniciación
de Científica. Universidad Nacional de Río Cuarto. Río
Cuarto, Argentina.
Gavrilovic, S. 1972. Inzenjering o bujicnim tokovima i
eroziji. Izgradnja, specijalno izdanje, 261-264. Beo- Tangestani,H. 2006. Comparison EPM and PSIAC
grad, Serbia and Montenegro. Model in Erosion and Sediment Yield Assessment
with GIS tools in Afzar Sub-Catchments, Fars Prov-
Gavrilovic, Z.,1988. The use of an empirical method
ince, Iran, Elsevier Magazine, Vol. 27 (5): 585-597.
(erosion potential method) forcalculating sediment
production and transportation in unstudied or torren- Tosic, R., Dragicevic, S., Zlatic, M., Todosijevic, M. &
tialstreams. In: International conference of river re- S. Kostadinov, 2012. The impact of socio-demogra-
gime, 18–20 May. John Wileyand Sons, Wallingford, phic changes on land use and soil erosion (Case stu-
Chichester, pp. 411–422. dy: Ukrina river catchment). Географски разгледи
(46):69-78.
INA, 2007. Proyecto “Estudios de erosión, transporte
y sedimentación en la cuenca alta del Río Pilco- Ulla, J.P. 2008. Variación de parámetros sedimento-
mayo”: determinación de la producción de sedimen- lógicos y su relación con la descarga e intervencio-
tos por erosión superficial en la cuenca alta del Río nes de cauce en la cuenca media alta del río Cho-
Pilcomayo. INA (Instituto Nacional del Agua, Argen- cancharava, Córdoba. Tesis de Licenciatura. Dpto.
tina). de Geología, Universidad Nacional de Río Cuarto,
Río Cuarto, Argentina. 159 p.
Olmos, L.A., Spalletti, P., Farías, H.D. & M.T. Pilán
2008. Verificación de modelos de producción de UNESCO. 2010. Procesos de erosión - sedimenta-
sedimentos en cuencas del NO Argentino. Revista ción en cauces y cuencas. Volumen 1. Programa Hi-
de Ingeniería del Agua, Vol. 15, N°2. drológico Internacional de la UNESCO para América
Latina y el Caribe. Documento Técnico Nº 22. Brea,
Panagos, P., Karydas, C., Ballabio, C. e I. Gitas,
J.D. & F. Balocchi (Eds.), 138 p.
2014. Seasonal monitoring of soil erosion at regional
scale: An application ofthe G2 model in Crete focus- Vezzoli, G., Ghielmi, G., Mondaca, G., Resentini, A.,
ing on agricultural land uses. International Journal of Villarroel, E.K., Padoan, M. & P. Gentile 2013. Quan-
Applied Earth Observation and Geoinformation 27 tifying modern erosion rates and river-sediment con-
(2014) 147–155 tamination in the Bolivian Andes. Journal of South
American Earth Sciences 45: 42-55.
Secretaría de Ambiente y Desarrollo Sustentable,
Argentina. 2009. Estadística de incendios forestales Zemljic M. 1971. Calcul du debit solide. Evaluation
2008. 01/09/13, www.ambiente.gob.ar/archivos/web/ de la vegetation comme un des facteurs antierosif.
PNEF/file/Incendios%202008/incend-2008-%20cor- International Symposium Interpreevent, Villach, Aus-
reg3/Documento%20completo.pdf tria.

Aqua-LAC - Vol. 6 - Nº. 1 - Mar. 2014 49