ARTICULO ORIGINAL

 

Efecto de las técnicas de mejora ecohidrológica del pastizal sobre el rendimiento hídrico de la microcuenca alto andina Urpay

Effect of grassland ecohydrological improvement techniques on the water performance of high andean micro watershed Urpay

 

Luis Oscanoa¹ y Enrique Flores²

¹ Ph.D. Investigador Asociado a The Mountain Institute y Laboratorio de Ecología y Utilización de Pastizales, Departamento de Producción Animal de la Facultad de Zootecnia, Universidad Nacional Agraria La Molina, Av. La Molina s/n La Molina, Lima, Casilla Lima 12, Perú, oscanoa.luis@gmail.com
² Ph.D. Profesor Principal. Laboratorio de Ecología y Utilización de Pastizales, Departamento de Producción Animal de la Facultad de Zootecnia, Universidad Nacional Agraria La Molina, Av. La Molina s/n La Molina, Lima, Casilla Lima 12, Perú, lab_pastizales@lamolina.edu.pe

 

---

Resumen

El presente trabajo fue realizado en el ecosistema de montañas en deterioro de la comunidad campesina Cordillera Blanca, Huaraz – Perú. El objetivo fue desarrollar una propuesta de estudio basada en la aplicación integrada de los índices de la condición del pastizal (C), humedad del suelo (H), caracterización ecohidrológica y mejora del rendimiento hídrico (RH) con "surcos" y "hoyos". La C fue estudiada con el método Parker (1951), la H con el diseño experimental completamente al azar y arreglo factorial (3 técnicas x 6 épocas por parámetro de estudio) y el RH con el aplicativo HEC-HMS. La investigación duró dos años, etapa en la que se estudiaron la mejora de la C y H, atributos base para el mapeo ecohidrológico (Vulnerabilidad, conservación, función hidrológica, retención de agua) y el RH. La C muy pobre no mostró cambios pero sí la H (P=0.000<0.05), siendo los "surcos" superior (14.66%) a "hoyos" (12.95%) y control (10.95%); el mapeo ecohidrológico del territorio (2 449 ha) identificó aptitudes de calidad media (49%), alta (29%) y pobre (22%); y los "surcos" perdieron menos caudal de agua (1.17 m³.s^-1) que "hoyos" (1.19 m³.s^-1) y el control (1.20 m³.s^-1). Se concluye que la caracterización ecohidrológica proporcionó adecuados indicadores para el diseño del plan de conservación, restauración y monitoreo de la microcuenca; los índices C y H son muy confiables y los "surcos" y "hoyos" son técnicas efectivas para la restauración de la función ecohidrológica del pastizal.

Palabras clave : pastos naturales, microcuenca, ecohidrología, rendimiento hídrico.

 

---

Abstract

The present work was performed in a deteriorated mountain ecosystem in a peasant community from the Cordillera Blanca, Huaraz - Peru. The objective of the experiment was to develop a study proposal based on the integrated application of pasture condition indexes (C), soil humidity (H), ecohydrological characterization (EH) and improved water performance (RH) using "furrows" and "pits". The C was studied with the Parker method (Parker, 1951), the H with a completely randomized experimental design with a factorial arrangement (3 x 6 times techniques by study parameter) and the RH with the HEC-HMS applicative. The study lasted two years, period in which the improvement of C and H were studied, which constituted base attributes for the ecohydrological mapping (Vulnerability, conservation, hydrological function, water retention) and RH. A very poor C did not show changes but H did (P = 0.000 <0.05). "Furrows" were higher (14.66%) than "pits" (12.95%) and control (10.95%). Ecohydrological mapping of the area (2 449 ha) identified medium (49%), high (29%) and poor (22%) qualities. "Furrows" lost less water flow (1.17 m³.s^-1) than "pits" (1.19 m³.s^-1) and control (1.20 m³.s^-1). It is concluded that ecohydrological characterization provided appropriate indicators for designing conservation plans, restoration and monitoring of micro watersheds; C and H indexes are of high confidence. "Furrows" and "pits" are effective techniques for the restoration of the water function of the pasture.

Key words: grassland, micro watersheds, ecohydrological, water performance.

 

---

Introducción

En cuencas hidrográficas, se han ensayado modelos hidrológicos para estimar la descarga de agua y erosión del suelo. Los parámetros utilizados para el análisis del ordenamiento territorial fueron: biomasa de cultivos (Pieri et al., 2007), precipitación-escorrentía (Amaya et al., 2009), erosión del suelo (Nearing et al., 2011), cobertura vegetal (Di Luzio et al., 2002) y uso de la tierra (Albhaisi et al., 2013). En este sentido, Nearing et al. (2011) desarrollaron el Range Hidrologic Erotion Model (RHEM) con el fin de modelar la erosión del pastizal basado en informaciones de tierras de cultivo (escorrentía y erosión) ellos consideraron de que las tasas de pérdida de suelo son uno de los pocos indicadores para evaluar la salud del pastizal y la eficiencia práctica de su conservación.

Uno de los problemas preponderantes en la gestión hídrica de ecosistemas de montañas es el sobrepastoreo del pastizal. En esta, las especies vegetales del pastizal afectan diferenciadamente el desarrollo de los parámetros composición florística, biomasa aérea y humedad del suelo del ecosistema pastizal y las correlaciones entre estas son cambiantes con la época del año (Eviner, 2004); de la misma forma, el estado fenológico de las plantas varía con la especie, el suelo y el clima (Tovar & Oscanoa, 2002). Por problemas de variabilidad del clima, el modelamiento requiere considerar un enfoque integral para el diseño de planes de control de los servicios ecosistémicos porque la duración y la intensidad de pastoreo incrementa la escorrentía y erosión del suelo, amenazas constantes para la desertificación y la condición de los pastizales (Zhang et al., 2012; Yalew, 2014).

Desde el punto de vista de la ciencia del pastizal, Florez & Bryant (1989) definieron como "condición" al diagnóstico de los cambios que ocurren con respecto a la comunidad clímax (situación de máxima relación trófica e equilibrio del ecosistema). Pyke et al. (2002) recomendaron incorporar para estudios de suelos, cubiertos de pastos naturales el enfoque ecológico, involucrando para ello atributos de estabilidad del sitio, función hidrológica e integridad biótica de la comunidad vegetal que anteriores metodologías no lo consideraron. Y desde la ciencia del suelo se propone que en los próximos estudios de la condición ecológica del pastizal se integre al modelo de evaluación el factor humedad del suelo que muy pocas veces se tuvieron en cuenta (Eekeren et al., 2010; Oscanoa & Flores, 2016). Snymana & Preez (2005) encontraron que los cambios en los niveles de humedad del suelo varían en concordancia al grado de degradación del suelo y siendo la lluvia la única fuente de recarga de agua en el pastizal estas muestran influencia significativa en la cantidad de humedad del suelo (Salve & Allen-Diaz, 2001).

El objetivo de este trabajo fue desarrollar una propuesta basada en la aplicación integrada de los principales atributos de la ciencia del pastizal (condición - C) y suelo (humedad - H), caracterización ecohidrológica y simulación de mejora del rendimiento hídrico (RH) de la microcuenca con las técnicas de "hoyos" y "surcos" de infiltración. Las preguntas de investigación fueron ¿cuál es el nivel de influencia de los atributos de la condición del pastizal y humedad del suelo?, ¿son los submodelos de vulnerabilidad, conservación, función hidrológica y retención del agua adecuados para el estudio ecohidrológico del pastizal?, ¿cuál es el nivel de mejora del rendimiento hídrico con las técnicas "hoyos" y "surcos"?

Materiales y métodos

Diseño del estudio y participantes

Llevamos a cabo este trabajo en pastos naturales alto andinos (4 100 msnm) entre mayo 2011 y diciembre 2014 en terrenos de la comunidad campesina Cordillera Blanca, Sector Canrey Chico, Provincia de Recuay, Huaraz, Perú (Figura 1). La delimitación de la subcuenca Urpay se realizó participativamente con ganaderos de la comunidad y el apoyo de imágenes satelitales Landsat con sensor TM del año 2005 aportada por el Ministerio de Educación del Perú (MINEDU). Holdridge (1987) clasificó a esta zona como Páramo muy Húmedo-Subalpino Tropical (Pmh- SaT) con temperatura promedio de 12 °C, precipitación promedio de 700 mm por año y suelo franco arenoso con tipo de pastizal Calamagrostis macrophylla y Stipa brachyphylla de condición muy pobre.

 

Población del estudio

La unidad de análisis del comportamiento hídrico del pastizal fue la microcuenca Urpay, ubicada a una altitud de 3 450 - 4 800 msnm y con punto de descarga tributaria en la coordenada UTM 230800 y 8926312 (Figura 1). En la cabecera del territorio nace el rio Urpay generado por la caída de lluvias. En la Figura 2 se presenta las precipitaciones registradas por las estaciones pluviométricas Pocco I y Pocco II, implementadas durante el periodo de estudio. El caudal base fue sostenido en el tiempo por pequeñas lagunas, puquiales, bofedales y filtraciones que proveyeron agua al rio, las que finalmente como tributario desembocaron en el rio Santa. La microcuenca fue considerada representativa del sistema de cabeceras de cuenca de la cordillera blanca y proporcionó información confiable para la simulación del comportamiento hidrológico del ecosistema pastizal alto andino. Es importante este trabajo debido a que la influencia de los atributos del estado ecohidrológico sobre la función productiva de la microcuenca se encuentra inmersa dentro un sistema dinámico de alta complejidad.

 

Tratamientos Experimentales

Se ensayaron 18 tratamientos producto de la combinación de tres técnicas y seis épocas de muestreo. El monitoreo detallado de los tratamientos se realizó bajo el modelo sistémico aleatorio en parcelas ("surcos", "hoyos" y uso comunal) en las épocas de fin de lluvias (may-11), plena seca (ago-11), inicio de lluvias (oct-11), plena lluvia (feb-13) y plena seca (set-13) consideradas estaciones claves para el manejo ganadero. Estos tratamientos fueron estudiados en parcelas cercadas 25 x 25 m, donde los surcos (15 cm de ancho por 10 cm de profundidad) fueron construidos a intervalos de dos metros y los hoyos (7.5 cm diámetro por 10 cm de profundidad) a intervalos de un metro, en referencia al área de pastoreo libre comunal considerada como testigo (Figura 1).

Parámetros Evaluados

La condición del pastizal de las parcelas se evaluó al inicio y final del experimento con la metodología propuesta por Oscanoa & Flores (2016) a través de los índices de especies deseables (D en %), índice forrajero (IF = % de especies deseables más poco deseables), índice de suelo desnudo más pavimento de erosión más roca (BRP en %) y el índice de vigor (V en %) de las especies vegetales claves (Parker, 1951). La integración matemática de los índices parciales de condición se realizó con la Fórmula:

Puntaje = 0.5D + 0.2IF + 0.2COB + 0.1V
Donde:
Índice de cobertura (COB) = 100 - BRP

El puntaje total se comparó con la escala puntaje – condición (Tabla 1) propuesta por Florez et al. (1991).

 

Las muestras de suelo para el análisis de la humedad del suelo fueron obtenidos con cilindros de muestreo de 7.62 cm de altura y 7.62 cm de diámetro a 10 y 20 cm de profundidad del suelo. Con una comba se introdujo el cilindro en el suelo, se uniformizó la superficie de los extremos del cilindro con cuchillo plano, se sellaron con papel de hornear y las muestras codificadas fueron trasladadas al laboratorio (Oscanoa & Flores, 2016). En el laboratorio, las muestras frescas fueron pesadas, oreadas y secadas en estufa a 105 °C por 24 horas y luego las muestras secas fueron pesadas en balanza de precisión y los resultados expresados en porcentaje (Blake et al., 1965).

La caracterización ecohidrológica se realizó con el objeto de integrar los parámetros ecológicos e hidrológicos del pastizal para la restauración de las funciones ambientales del ecosistema; en razón del cual se generó los submodelos de vulnerabilidad, estado de conservación, función hidrológica y capacidad de retención de agua por el suelo (Pyke et al., 2002; Oscanoa, 2016); el primero evaluó el daño provocado por el hombre, en razón del cual se consideraron los indicadores de fisionomía vegetal, pendiente y precipitación; el segundo submodelo midió el efecto del sobre pastoreo sobre el estado de la comunidad vegetal, a través de la condición del pastizal, modalidad de tenencia de tierras y práctica local de conservación de suelos; en el tercer submodelo se analizaron el efecto del comportamiento biofísico del agua de lluvias sobre la producción de biomasa, humedad del suelo, producción de mantillo y temperatura del suelo (Oscanoa & Flores, 2016); y en el último submodelo se estudiaron la capacidad de retener el agua de lluvias en contra de la fuerza de gravedad a través de la erosión del suelo, tasa de infiltración y análisis de textura del suelo.

El modelamiento ecohidrológico de la subcuenca se realizó con el software ENVI 5.0 y ArcGis 10.0; con el primero se hizo la fotointerpretación del área de estudio y con el segundo el desarrollo de los mapas temáticos y el modelamiento final. Se procesó las características ecohidrológicas de la subcuenca tomando como indicadores el ponderado de ponderados de los submodelos temáticos de vulnerabilidad, conservación, función hidrológica y retención del agua. El proceso determinó las características hídricas de las unidades zonificadas que luego fueron sometidos a simulación de mejora de las áreas con potencial ecohidrológico de la subcuenca.

El grado de influencia de las prácticas de conservación de suelos ("hoyos", "surcos") sobre la variable caudal fueron estimados con el software HEC- HMS; de modo que los atributos fueron analizados detalladamente con la finalidad de caracterizar realísticamente la tipología ecohidrológica y sus implicancias sobre el rendimiento hídrico de la microcuenca. Los parámetros considerados fueron: (1) datos de la superficie, tasa de pérdida determinado con el método curva número (Balvanshi & Tiwari, 2014), conversión de lluvia en escurrimiento con el método del Unit Hydrography (SCS) y el flujo base constante fue el valor más bajo de la época seca registrados por los vertederos de la estación Pocco (Gil, 2014), (2) la curva número ponderada fue estimada con el ArcGis 10.0, la impermeabilidad considerada fue del 5% y el tiempo de concentración (Tc) fue calculado con la ecuación Kirpich (Vélez & Botero, 2011; Taghvaye et al., 2017):

Tc = 0.000323(0.77L/0.385S)
Dónde:

L es la longitud del cauce, S es la semisuma de la cota máxima y mínima dividido entre L. Y para el parámetro tiempo de retardo el valor fue considerado como el 60% del tiempo de concentración.

Diseño Experimental

Se empleó el Diseño Completamente al Azar con arreglo factorial con dos factores A (dos técnicas de mejora y testigo) y B (seis meses de muestreo), donde los tratamientos considerados fueron la interacción de las técnicas con los meses de muestreo, cuyas combinaciones hicieron un total de 18 tratamientos. El Diseño Factorial fue no balanceado (desigual cantidad de observaciones por casilla) y fue completo, ya que se realizaron todas las posibles combinaciones entre los valores de las variables (Montgomery, 1991). Los promedios de los parámetros estudiados fueron comparados con el test de comparaciones múltiples Duncan (p<0.05).

Modelo Estadístico:

Yijk = µ + αi +τj + ατij + eijk
i = 1, 2 ,3 ; j = 1, 2 ,3 ,4,5,6 ; k = 1, 2 ....,nij
Donde:

Yijk es la respuesta obtenida (rendimiento de parámetros hidrológicos del pastizal) en la unidad experimental correspondiente a la k-ésima repetición a la cual se le aplicó el i-ésimo nivel del factor A, con el j-ésimo nivel del factor B. µ es el efecto de la media general. αi es el efecto del i-ésimo nivel del factor técnica de mejora (A). τj es el efecto del j-ésimo nivel del factor época de muestreo (B). ατij es el efecto de la interacción entre el i-ésimo nivel del factor A con el j- ésimo nivel del factor B. eijk es la variación del error asociado con las ijk unidades.

Resultados y discusión

Humedad del suelo y condición del pastizal

En dos años de estudio la condición del pastizal muy pobre para ovinos, vacunos y vicuñas no lograron mejora con ninguna de las técnicas aplicadas ("hoyos" y "surcos"); en cambio la humedad del suelo de "surcos" mostró superioridad significativa en los factores técnica, mes y su interacción (R²=0.812; p<0.05). Al comparar el factor mes (noviembre 2012 y febrero 2013) con los otros meses estudiados, estas mostraron diferencias estadísticas significativas en todas sus combinaciones (Tabla 2).

 

Estos datos sugieren que los atributos que evalúan el pastizal en deterioro, basado en la composición florística, están sujeto a una lenta recuperación del vigor de las plantas de importancia ganadera al aplicarse técnicas de mejora; sin embargo, el factor humedad del suelo se muestra más sensible a cambios que ocurren en el suelo; en ese sentido, los "surcos" mostraron superioridad sostenida en todos los meses estudiados y alcanzó el pico más alto durante el mes de febrero 2013 (Figura 3). Los resultados obtenidos expresan que los atributos ensayados son de alta complementariedad para la evaluación de la biomasa aérea del pastizal.

 

Deutsch et al. (2010) observaron que la humedad del suelo bajo sequía severa limita la producción de los pastizales templados-fríos, los cambios en los niveles de humedad del suelo variaron en concordancia al grado de degradación del suelo (Snymana & Preez, 2005), siendo la lluvia la única fuente de recarga de agua; la textura, cubierta vegetal, y elevación influyen en el contenido de humedad del suelo (Salve & Allen-Diaz, 2001). La dinámica de la humedad del suelo mostraron tendencias similares en los perfiles del suelo de bosques y del pastizal con secado más intenso en el bosque durante la temporada de crecimiento, posiblemente por la mayor utilización de agua en las capas por los árboles (Hagyó et al., 2006). De manera que la función de los tipos de vegetación es afectada por la distribución de agua disponible en el suelo y la composición florística del pastizal (Breshears & Barnes, 1999).

En los "hoyos", al final del experimento, hay un declive pronunciado del contenido de humedad del suelo, el cual podría deberse a que estas han perdido su función plena debido a que los hoyos se han rellenado nuevamente de suelo; este efecto es similar a la investigación comunicado por Al-Seekh et al. (2009) quienes observaron que la profundidad del suelo, precipitación y pisoteo del ganado son los principales factores para la pobre productividad de los pastizales. Zhang et al. (2010) determinaron que las tasas basales de respiración y temperatura varían con la gradiente de humedad, siendo la temperatura del suelo y contenido de agua los factores ambientales claves para la variación en la respiración del suelo; por lo que el manejo de los niveles de la humedad del suelo son parámetros importantes para el rendimiento de la función hidrológica del pastizal. En general, los resultados demostraron que las técnicas de mejora de suelos influenciaron en el incremento de la humedad del suelo y sustenta la confianza para la recuperación de las especies vegetales nativas de importancia ecohidrológica.

Modelo ecohidrológico de la microcuenca

El modelamiento determinó que la microcuenca Urpay ostenta sitios con vocación para la reserva de agua de lluvias. Proporcionalmente, de las 2 445 ha del territorio, la extensión de bofedales y humedales representó el 29% de la extensión, las zonas de "aptitud media" fue de 49% y el 22% de la superficie restante mostraron signos de franco proceso de deterioro (Tabla 3). Superficialmente, estas significaron que las áreas se agruparon en media (1 211 ha), baja (427 ha) y muy baja (117 ha), todas ellas con solo recarga de agua de lluvias y cubiertos por pastos naturales de puna y pajonales (Figura 4).

 

 

La ciencia ecohidrológica de montañas busca integrar los parámetros ecológicos e hidricos de los pastos naturales altoandinos con el propósito de mejorar las funciones ambientales dañados. En ese sentido dentro la Tabla 1 de valoración del modelo ecohidrológico resaltan por nivel de importancia las temáticas de vulnerabilidad e hidrología del pastizal (Tabla 3). La primera expresa el mayor peso del efecto de vulnerabilidad para el ecosistema pastizal (35%) ocacionados por la acción directa del hombre y la segunda el grado (30%) de influencia del deterioro autoregulatorio de la caída de lluvias en la recarga natural del agua del rio Urpay (Figura 4).

Estos datos confirmaron que la aptitud ecohidrológica de la microcuenca se encontró íntimamente vinculada al tipo de uso de los pastos naturales, fisiografía y clima (inundación, friaje y sequia). Con respecto a los parámetros en estudio se encontró que el proceso ponderativo de los sub modelos fueron representativos, pues el modelo consolidado representó realísticamente la aptitud ecohidrológica de la microcuenca; de manera que las zonas calificadas como "medio" se muestran parcelas adecuados para la introducción de técnicas de mejora y la aptitud alta como área natural para la conservación hidrológica. A medida que se mejora el estudio geoespacial del ecosistema se fortalece la calidad de información para el desarrollo del plan de uso de los recursos naturales en armonía con el estado de las fuentes de agua.

Los gestores de recursos naturales requieren generar modelos explícitos de respuesta hidrológica a cambios de conductores climáticos en un escenario donde las condiciones de humedad futura son más difíciles de proyectar que la temperatura (Flint et al., 2013). El problema principal para el modelamiento es la poca disponibilidad de información hidroclimática (Amaya et al., 2009); por lo que el estado de pastoreo de terrenos tendría cuantitativamente ser evaluados por el estudio edafo-agrostológico (Florez & Bryant, 1989) y la humedad del suelo (Snymana & Preez, 2005; Salve & Allen-Diaz, 2001); cualitativamente a través de los atributos de estabilidad del sitio, función hidrológica e integridad biótica de la comunidad vegetal (Pyke et al., 2002).

La existencia de factores de información por lo cual cada parámetro ecohidrológico pudo expresar su peso de influencia en el estudio de recarga de agua son índices dependientes del sistema hidrológico complejo que va más allá de la precipitación, tales como fuentes de deshielo del glaciar, filtraciones, lagunas, bofedales, y humedales. Los submodelos y modelos permitieron además ponderar adecuadamente el conocimiento local con propuestas metodológicas capaces de procesar integralmente las informaciones vinculadas a acuerdos comunales, testimonios de usuarios y percepciones locales.

Rendimiento hídrico de la microcuenca

La caracterización ecohidrológica de la microcuenca proporciono información razonable al modelo HEC-HMS para la estimación de descarga de agua de lluvias al implementarse las técnicas de mejora de suelos y calcular el efecto de mejora de la capacidad de retención de agua superficial de la microcuenca. En promedio la técnica "surcos" (1.17m³.s-¹) mostró un menor nivel de descarga en el punto final de recorrido del rio Urpay que los "hoyos" (1.19 m³.s-¹) y el control (1.20 m³.s-¹) (Tabla 4).

 

Cuantitativamente, las curvas de tendencia expresaron que la técnica "surcos" conserva mejor el volumen de agua al final de la época de lluvias (Figura 5). De manera que la introducción de técnicas de mejoramiento de suelo en el área de estudio ha mejorado la retención de agua superficial provenientes de las precipitaciones que no se infiltra ni se evapora en forma diferenciada durante la época de lluvias, los efectos coincidieron plenamente con el estado fenológico de floración de la comunidad vegetal y el nivel de humedad de suelo de la microcuenca (Figura 5).

 

Los resultados demostraron que si es posible estimar la descarga de agua superficial dentro el complejo sistema del pastizal. Yan et al. (2013) observaron semejantes resultados al presente ensayo y señalaron que los principales cambios que afectaron el caudal son los cambios en las tierras de cultivo, bosques y áreas urbanas, y los cambios en los pastizales no ejercieron influencia significativa en el caudal y producción de sedimentos. En contraposición a estos resultados, Stehr et al. (2010) manifestaron que la sustitución de los pastizales con plantaciones de especies introducidas causó una reducción de la media anual de descarga en un 10%. Finalmente, Nearing et al. (2011) indicaron la necesidad imperiosa de contar con un modelo para predecir la erosión del suelo de pastoreo porque la que existen se desarrollaron para estudiar tierras de cultivo con características de procesos hidrológicos diferentes al pastizal debido a la heterogeneidad de las propiedades y uso del suelo.

Conclusiones

El estudio basado en la aplicación integrada de los principales atributos de la ciencia del pastizal (condición - C) y suelo (humedad - H) encontró que los índices son muy confiables y precisa el proceso de análisis de los cambios en la composición aérea y subsuelo de la comunidad vegetal a corto y mediano plazo. La caracterización ecohidrológica de la microcuenca Urpay logró una adecuada delimitación del territorio por aptitud y cada respuesta fue fortalecida con el análisis individualizada de los mapas temáticos de vulnerabilidad, conservación, función hidrológica y retención de agua que interactuaron en la regulación de descarga del agua al rio. En los sitios ecohidrológicos, el modelo simulado encontró mejora del rendimiento hídrico con la aplicación de técnicas de conservación de suelos ("hoyos" y "surcos") y evidenciaron que si es posible mejorar, a mínimo costo económico y de tiempo, la función ecohidrológica de los servicios ambientales, planificación participativa del territorio y diseño de programas de gestión racional de microcuencas hidrográficas basadas en el conocimiento local y científico. Las limitantes de datos de precipitación detectados por los investigadores sugieren la necesidad de desarrollar pluviómetros de tecnología intermedia, confiables y de fácil manejo por el productor.

Agradecimientos

Al Comité de Agropecuaria Local "Alli Pastu" de Canrey Chico - Huaraz, a The Mountain Institute, al Profesor Juancarlos Cruz y al Laboratorio Utilización de Pastizales de la Universidad Nacional Agraria La Molina por el apoyo de análisis de muestras; a The Collaborative Crop Research Program/The MckNight Foundation por la beca doctoral y al Ph.D. Jorge Recharte por los apoyos y comentarios al contenido del presente artículo.

Literatura citada

Albhaisi M., Brendonck L. & Batelaan O. 2013. Predicted impacts of land use change on groundwater recharge of the upper Berg catchment, South Africa. Water SA, 39(2): 211-220. DOI: http://dx.doi.org/10.4314/wsa.v39i2.4Links ]doi.org/10.4314/wsa.v39i2.4" target="_blank">.

Al-seekh S.H., Mohammad A.G. & Amor Y.A. 2009. Effect of Grazing on Soil Properties at Southern Part of West Bank Rangeland. Hebron University Research Journal. Vol.(4), No.(1), pp.(35 – 53).         [ Links ]

Amaya G., Restrepo-Tamayo C., Vélez M.V., Vélez J.I. & Álvarez O.D. 2009. Modelación del comportamiento hidrológico de tres cuencas en el Urabá Antioqueño – Colombia. Avances en Recursos Hidraúlicos,19: 21-38. Medellín – Colombia.         [ Links ]

Balvanshi A. & Tiwari H.L. 2014. A Comprehensive review of runoff estimation by the curve number method. International Journal of Innovative Research in Science, Engineering and Technology, Vol. 3, Issue 11.         [ Links ]

Blake G.R. & Hartge K.H. 1965. Bulk Density. In: Methods of Soil Analysis, Part 1. Physical and Mineralogical Properties, Including statistics of Measurement and Sampling. Agronomy: Monograph 9: 374-390. ASA-SSSA.         [ Links ]

Breshears D.D. & Barnes F.J. 1999. Interrelationships between plant functional types and soil moisture heterogeneity for semiarid landscapes within the grassland/forest continuum: a unified conceptual model. Landscape Ecology 14: 465–478, 1999.         [ Links ]

Deutsch E.S., Bork E.W. & Willms W.D. 2010. Soil moisture and plant growth responses to litter and defoliation impacts in Parkland grasslands. Agriculture, Ecosystems and Environment 135: 1-9.         [ Links ]

Di Luzio M., Srinivasan R. & Arnold J.G. 2002. Integration of watershed tools and SWAT model into BASINS. JAWRA, 38(4): 1127-1141. DOI:https://doi.org/10.1111/j.1752-1688.2002.tb05551.x.         [ Links ]

Eekeren N.V., Boer H., Hanegraaf M., Bokhorst J., Nierop D., Bloem J., Schouten T., Goede R. & Brussaard L. 2010. Ecosystem services in grassland associated with biotic and abiotic soil parameters. Soil Biology & Biochemistry, 42(9): 1491-1504. DOI: 10.1016/j.soilbio.2010.05.016Links ]doi.org/10.1016/j.soilbio.2010.05.016" target="_blank">.

Eviner V. 2004. Plant traits that influence ecosystem  processes vary independently among species. Ecology, 85(8): 2215–2229. DOI: https://doi.org/10.1890/03-0405Links ]org/10.1890/03-0405" target="_blank">.

Flint E.L., Flint L.A., Thorne H.J. & Boynton R. 2013. Fine- scale hydrologic modeling for regional landscape applications: the California Basin Characterization Model development and performance. Ecological Processes 2013.         [ Links ]

Florez A. & Bryant E. 1989. Manual de Pastos y Forrajes. Programa colaborativo de apoyo a la investigación en rumiantes menores. INIIA, Universidad de California. Lima – Perú         [ Links ].

Florez M.A., Malpartida I.E. & San Martin F. 1991. Manual de forrajes para zonas áridas y semiáridas andinas.         [ Links ]

Gil R.J. 2014. Monitoreo Hidrológico en ecosistemas de la puna húmeda bajo el método de Cuencas Pareadas: C.C. Cordillera Blanca, Canrey Chico, Recuay, Ancash, Perú. Proyecto Puna-agua. Conv. Instituto de Montaña/ The MckNight Foundation. Documento en prensa.         [ Links ]

Hagyó A., Rajkai K. & Nagy Z. 2006. Effect of forest and grassland vegetation on soil hydrology in Mátra Mountains (Hungary). Biologia, Bratislava, 61/Suppl. 19:S261—S265.         [ Links ]

Holdridge L.R. 1987. Ecología basada en las zonas de vida. Instituto Interamericano de Cooperación para la Agricultura (IICA). San José de Costa Rica.         [ Links ]

Montgomery D.C. 1991. Diseño y análisis de experimentos. Arizona State University , Estados Unidos de América. Grupo Editorial Iberoamérica. Pag: 589.         [ Links ]

Nearing M.A., Wei H., Stone J.J., Pierson F.B., Spaeth K.E., Weltz M.A., Flanagan D.C. & Hernandez M. 2011. A rangeland hydrology and erosion model. Transactions of the ASABE, 54(3): 901-908. American Society of Agricultural and Biological Engineers. DOI:10.13031/2013.37115Links ]org/10.13031/2013.37115" target="_blank">.

Oscanoa L. 2016. Influencia de las prácticas de conservación de suelos y mejora de la función hidrológica de praderas naturales altoandinas. Tesis para optar el grado Doctoris Philosophiae en Ciencia Animal. Universidad Nacional Agraria La Molina. Lima, Perú         [ Links ].

Oscanoa L. & Flores E. 2016. Influencia de técnicas de mejora de suelos sobre la función hídrica de pastos naturales altoandinos. Ecol. Apl. Vol. 15 N° 2, pp. 91-99.         [ Links ]

Parker K.W. 1951. A method for measuring trend and range condition on National Forest Ranges. USDA. Forest Service.         [ Links ]

Pieri L., Bittelli M., Wu Q.J., Dun S., Flanagan C.D., Rossi P.P., Ventura F. & Salvatorelli F. 2007. Using the Water Erosion Prediction Project (WEPP) model to simulate field-observed runoff and erosion in the Apennines mountain range, Italy. Journal of Hydrology, 336: 84–97.         [ Links ]

Pyke D.A., Herrick J.E., Shaver P. & Pellant M. 2002. Rangeland health attributes and indicators for qualitative. J. Range Manage, 55(6): 584-597. DOI:10.2458/azu_jrm_v55i6_pykeLinks ]org/10.2458/azu_jrm_v55i6_pyke" target="_blank">.

Salve R. & Allen-Diaz B. 2001. Variations in soil moisture content in a rangeland catchment. J. Range Manage. 54(1): 44–51. DOI: 10.2307/4003526Links ]org/10.2307/4003526" target="_blank">.

Snymana H. & du Preez D. 2005. Rangeland degradation in a semi-arid South Africa—II: influence on soil quality. Journal of Arid Environments, 60(3): 483–507. DOI: https://doi.org/10.1016/j.jaridenv.2004.06.005Links ]org/10.1016/j.jaridenv.2004.06.005" target="_blank">.

Stehr A., Aguayo M., Link O., Parra O., Romero F. & Alcayaga H. 2010. Modelling the hydrologic response of a mesoscale Andean watershed to changes in land use patterns for environmental planning. Hydrol. Earth Syst. Sci., 14, 1963–1977.         [ Links ]

Taghvaye S.E., Nohegar A., Malekian A., Hoseini M. & Holisaz A. 2017. Estimating time of concentration in large watersheds. Paddy Water Environ 15:123–132.         [ Links ]

Tovar O. & Oscanoa L. 2002. Guía para la identificación de pastos naturales alto andinos de mayor importancia ganadera. Instituto de Montaña. Proyecto FOCAL Obc. Primera edición. Huaraz, Perú         [ Links ].

Vélez U.J. & Botero G.A. 2011. Estimación del tiempo de concentración y tiempo de rezago en la cuenca experimental urbana de la quebrada San Luis, Manizales. Dyna, 78(165): 58•71. Medellín.

Yalew G.S., Tobias P., Schweitzer C., Liersch S., van der Kwast J., Mul M.L., van Griensven A. & van der Zaag P. 2014. Dynamic Feedback between Land Use and Hydrology for Ecosystem Services Assessment. In: Ames D.P., Quinn N.W.T. & Rizzoli A.E. (Eds.). Proceedings of the 7th International Congress on Environmental Modelling and Software (iEMSs), June 15-19, San Diego, California, USA. URL: https://scholarsarchive.byu.edu/iemssconference/2014/Links ]byu.edu/iemssconference/2014/" target="_blank">.

Yan B., Fang N.F., Zhang P.C. & Shi Z.H. 2013. Impacts of land use change on watershed streamflow and sediment yield: An assessment using hydrologic modelling and partial least squares regression. Journal of Hydrology, 484: 26–37. DOI:https://doi.org/10.1016/j.jhydrol.2013.01.008Links ]org/10.1016/j.jhydrol.2013.01.008" target="_blank">.

Zhang L.H., Chen Y.N., Zhao R.F. & Li W.H. 2010. Significance of temperature and soil water content on soil respiration in three desert ecosystems in Northwest China. Journal of Arid Environments, 74(10): 1200 - 1211. DOI:https://doi.org/10.1016/j.jaridenv.2010.05.031Links ]org/10.1016/j.jaridenv.2010.05.031" target="_blank">.

Zhang Y., Hernandez M., Anson E., Nearing M.A., Wei M.A., Stone J.J. & Heilman P. 2012. Modeling climate change effects on runoff and soil erosion in southeastern Arizona rangelands and implications for mitigation with conservation practices. Journal of Soil and Water Conservation, 67(5): 390-405. DOI:10.2489/jswc.67.5.390        [ Links ]

 

Presentado: 16/0/2018

Aceptado: 16/01/2019