1. INTRODUCCIÓN
El manejo integral del agua es un proceso que permite mantener y mejorar el estado de los recursos naturales a través de la implementación de enfoques de gestión adaptativa y participativa de los beneficiarios. Últimamente, se ha notado que la gestión del agua se va adecuando en base a la demanda del recurso, debido al crecimiento de la población [1]. Como ejemplo se observa el Valle Central de Cochabamba, que frente a la variación espacial y temporal del ciclo del agua y el aumento en la demanda de este recurso es que se desarrolló el Plan Director de la cuenca del río Rocha. [2].
La modelación hidrológica de cuencas es una necesidad crucial en regiones con manejo integral de recursos hídricos y problemas de abastecimiento de agua [3]. Se entiende por modelación hidrológica a una herramienta para la simulación de procesos hidrológicos. Hoy en día se emplean estos modelos para el análisis de inundaciones, apoyo a la toma de decisiones, propuestas de ordenamiento territorial, proponer el diseño de obras hidráulicas y establecer protocolos de actuación ante posibles situaciones de peligro por intensas lluvias [4]. Existen diferentes modelos hidrológicos, de acuerdo a su representación espacial, pueden ser semi distribuidos o distribuidos. Los modelos distribuidos permiten considerar la variación espacial de los procesos hidrológicos como la precipitación [5].
Como ejemplo de un modelo semi distribuido se tiene el Sistema de Modelación Hidrológico del Centro de Ingeniería Hidrológica del Cuerpo de Ingenieros del Ejército de los EE. UU. (HEC-HMS por sus siglas en inglés) [6]. Este modelo utiliza métodos de precipitación-escorrentía para estimar los hidrogramas de escorrentía directa generados por las precipitaciones en una cuenca o región durante un período especificado. Es un programa libre y gratuito, muy flexible que permite al usuario la selección de diferentes métodos para el cálculo de pérdidas, hidrogramas, flujo base y propagación en cauces; posibilitando la realización de simulaciones de los procesos hidrológicos con respecto a eventos o en forma continua. [7].
Dentro de los datos de entrada necesarios en la modelación se encuentra la precipitación, la cual es una variable determinante para realizar estudios hidrológicos, teniendo como características la resolución espacial y temporal. La captura de datos de precipitación se realiza a través de diferentes métodos, ya sean directos o indirectos. Los métodos directos emplean pluviómetros, herramientas que miden la intensidad de precipitación en intervalos de horas o días. [8].
En el caso de Bolivia, el Servicio Nacional de Meteorología e Hidrología (SENAMHI) es el encargado de captar datos de precipitación y temperatura, mediante estaciones meteorológicas y pluviómetros. En todo el territorio boliviano se tiene 547 estaciones pluviométricas en operación continua [9], las cuales en su mayoría se encuentran cercanas a los ejes metropolitanos debido a que requieren mantenimiento, control y manejo [10]. Motivo por el cual no se logra abarcar todo el territorio, por lo que es necesario utilizar métodos indirectos de captación de precipitación. Los métodos indirectos emplean estimaciones basadas en el espectro electromagnético. Por ejemplo, los radares terrestres o sensores a bordo de satélites utilizan diversos algoritmos los cuales ayudan a determinar estimaciones de la precipitación, mediante la cantidad de agua precipitable. Algunos ejemplos que se pueden mencionar son: CHIRPS, GSMaP, IMERG y PERSIANN.
PERSIANN es un producto que emplea el sensor infrarrojo GridSat-B1, la precipitación estimada se obtenía a partir de la información de detección utilizando la neurología sintética: los registros de datos climáticos (PERSIANN-CDR) proporciona estimaciones de lluvia diarias con una resolución espacial de 0.25° en la banda de latitud 60 S - 60 N desde 1983 hasta el presente cercano [11]. IMERG es un algoritmo estadounidense estandarizado por el grupo GPM, usando la versión 2014 del algoritmo de generación de perfiles de Goddard (GPROF2014), luego se cuadriculan, se inter calibran con el producto de instrumento combinado GPM y se combinan en la mitad - Campos de 0,1 ° x 0,1 ° por hora [12]. CHIRPS es un producto satelital el cual obtiene datos de precipitación infrarroja con estación Climate Hazards Group, es el conjunto de datos de precipitación a nivel mundial durante más de 30 años. CHIRPS combina imágenes satelitales de 0,05° con datos de estaciones in situ para crear series temporales superpuestas de precipitaciones para analizar tendencias y monitorear sequías estacionales [13]. GSMaP, el cual es un producto de la misión “Global Precipitation Measurement” (GPM), este proporciona una tasa de lluvia global por hora, con una resolución de 0,1° x 0,1°. que proporciona observaciones de precipitación global en intervalos de tres horas. Los valores se estiman utilizando radiómetros pasivos de microondas e infrarrojos multibanda del satélite GPM Core Observatory y con la ayuda de una constelación de otros satélites. El algoritmo de recuperación de la tasa de precipitación de GPM se basa en un modelo de transferencia [14].
La aplicación de sensores satelitales permite estimar la precipitación distribuida sobre la superficie terrestre. Las estimaciones se obtienen a partir de algoritmos de interpretación de las señales de los sensores. Para reducir la incertidumbre es recomendable generar un producto combinado entre las mediciones por los pluviómetros y estimaciones de precipitación satelital. Además, que con la aplicación de un producto combinado en una modelación hidrológica se pueden evaluar los caudales simulados para determinar eventos extremos y parámetros relacionados al caudal ecológico.
El objetivo de este estudio es la aplicación de productos combinados de sensores de precipitación a bordo de satélites con datos de pluviómetros y la modelación hidrológica con análisis de parámetros del caudal ecológico del Río Rocha para el análisis de eventos extremos.
2. ZONA DE ESTUDIO
La cuenca del río Rocha tiene una extensión aproximada de 3 700 km2. Está ubicada en el Departamento de Cochabamba, en ella habitan más de 1.4 millones de personas, lo que representa el 13% de la población del país. La cuenca enfrenta muchos desafíos relacionados con la gestión del agua.
La zona de estudio se divide en tres subcuencas: Maylanco, Valle bajo y Sulty, como en la Figura 1, las cuales integran un total de 186 microcuencas y unidades hidrológicas. La subcuenca Rocha es la que tiene la mayor densidad poblacional de 840 hab/km2, la subcuenca Maylanco tiene 366 hab/km2 y la subcuenca Sulty tiene una densidad promedio de 72 hab/km 2[2].
3. METODOLOGÍA
Para la ejecución de este trabajo se siguió el esquema de la Figura 2, donde se muestras los procesos de combinación de precipitación, configuración del modelo HEC-HMS y la aplicación del paquete IHA7.
3.1. Productos de precipitación utilizados
Los productos satelitales utilizados en este estudio cuentan con diferentes resoluciones como se observa en la Tabla 1. Debido a ello y con ayuda del programa de ArcGIS, se logró generar una grilla general para todo el proceso. Para esto se fijó una grilla con resolución espacial de 0.0025° o 250 m y una escala temporal diaria.
TABLA 1 CARACTERÍSTICAS DE PRODUCTOS DE PRECIPITACIÓN
Producto | Resolución Temporal | Resolución Espacial | Metodología Aplicada | Tecnología Aplicada | Formato de Salida | Referencia |
---|---|---|---|---|---|---|
CHIRPS | Diario | 0.05° | Ajuste y Calibración CMORPH | Infrarrojo y Estaciones Terrestres | TIFF | [13] |
GSMaP | Horario y Diario | 0.1° | Ajuste y Calibración GPM | Radiómetro de Microondas e Infrarrojo | DAT | [14] |
IMERG | Diario | 0.1° | Ajuste y Calibración TRMM | Sensores de Microondas e Infrarrojo | NETCDF | [12] |
PERSIANN | Diario | 0.25° | Ajuste y Calibración CMORPH | Redes Neuronales Artificiales e Infrarrojo | TIFF | [11] |
Fuente: Elaboración Propia.
Antes del proceso iterativo de cada uno de los productos, se hizo la selección de estaciones pluviométricas, las cuales se escogieron por tres criterios: el primero fue que la ubicación de cada estación se situara dentro de los límites de la cuenca del río Rocha, el segundo fue que tengan la mayor cantidad de datos entre los períodos del 2014 al 2020 y, por último, que los meses faltantes no sobrepasen los 3 meses seguidos vacíos. Bajo esas condiciones, se obtuvieron un total de 13 estaciones que fueron empleadas en el estudio.
3.2. Rellenado de datos de Pluviómetros
Después del proceso de selección, las 13 estaciones fueron separadas en dos grupos, esto debido al límite de estaciones programadas en HEC4 es de 10 estaciones como máximo. Se logró rellenar los días faltantes en las estaciones, además, para comprobar la calidad de los datos obtenidos, se comparó con la base de datos empleados en el Balance Hídrico Superficial de Bolivia del 2018.
3.3. Proceso de combinación de productos
Una vez preparados los datos terrestres y formateados los productos satelitales, se comenzó con el proceso de combinación de productos, en el que mediante lenguajes de programación se obtuvo el error relativo y se sacó un promedio para cada subcuenca, posterior a esto se sacó el factor de ajuste para el producto. Luego se aplicó este factor al mapa de precipitación satelital por subcuenca, cinco iteraciones fueron necesarias similar al estudio de Katari [15]. En la Figura 3 se puede ver el proceso que se aplicó en la combinación de productos. Para este proceso se generaron productos con dos escalas temporales diferentes, unos a escala mensual en los que se utilizó los cuatro productos satelitales y otros a escala diaria en los que se usó los 2 mejores productos del proceso de combinación anterior.
Los indicadores escogidos fueron: Coeficiente de correlación y determinación, error absoluto medio, error cuadrático medio y eficiencia de Nash según lo sugerido por la bibliografía [16].
3.4. Modelación en HEC-HMS
Para la configuración del Modelo Hidrológico con HEC-HMS, se optó por hacer un modelo continuo aprovechando la generación de los productos combinados, el cual permitió mostrar la generación de caudales en todo el período de estudio. En la Figura 4 se puede observar esquema de cómo se realizó la modelación. Para este proceso, primeramente, se determinó los períodos con los que se harán las respectivas calibraciones y validaciones, optando por tomar los períodos de 2017-2018 para la calibración, debido a variación en sus mediciones por la actividad humana en la zona y el período de 2014 para la validación, debido a que no hubo interferencia humana en la medición de ese año.
En la Figura 5, se observa el mapa de subcuencas de nivel 8, el cual cuenta con 59 subcuencas que fueron empleadas para realizar el proceso de modelación. En el presente modelo se fijaron cada una de las subcuencas, como también uniones y tramos importantes, además se optó por colocar reservorios en el caso de la presa México.
En la configuración del modelo HEC-HMS, se utilizaron los métodos de cobertura y superficie simples, método de pérdida por déficit y constante, para la obtención de estos valores se obtuvieron datos del documento de “Modelación Estratégica de la Cuenca del Río Rocha” [17] aplicado en el Balance Hídrico Superficial Nacional del 2018, los datos que se obtuvieron de dicho documento fueron: evapotranspiración, infiltración, perdidas en la región (por consumo humano o riego) presentes en las 3 grandes cuencas, debido a esto se aplicó en cada subcuenca el factor de peso por áreas, obtenidos de los shapes de la cuenca, para los valores de cada microcuenca.
Una vez efectuados los procesos de calibración y validación, se realizó la modelación con datos de tierra de todo el período de estudio, es decir de 2014-2020, posterior a la modelación con productos de pluviómetros, se aplicaron los mismos valores y métodos para la obtención de caudales de productos combinados de CHIRPS y GSMaP. Finalmente, se realizó una comparación entre productos, teniendo como base los caudales obtenidos a partir de pluviómetros.
3.5. Aplicación del paquete IHA7
Para determinar los índices de caudal ecológico se siguió el esquema de la Figura 6. Una vez seleccionado el mejor producto combinado se inicia la aplicación del programa IHA7, para este análisis se escogieron 6 puntos importantes, los cuales fueron: Sacaba, Tiquipaya, Punata, Cliza, Entrada a La Angostura y la salida de la cuenca del Río Rocha.
El programa IHA7 que se basa en la metodología de Richter [18], analiza los caudales obtenidos y determina los posibles eventos extremos que puedan ocurrir, este programa estudia un total de 33 índices. Estos se encuentran divididos en 5 grupos: en el primero se tiene las medias de cada mes que se ingresó y la variabilidad que tienen con respecto a su media. En el segundo grupo se obtienen los mínimos y máximos diarios, cada 3 días, semanal, mensual y estacional, y las variaciones con respecto a sus valores. En el tercer grupo se marca la fecha en la que podrían pasar los eventos más extremos como lo son de mayor y menor caudal, lo que puede generar una sequía o una inundación, este resultado está en días del calendario Juliano. En el cuarto grupo se tiene el número de veces al año que puede ocurrir ciertos eventos y la duración de estos, y el valor de los umbrales que se pueda tener, cabe recordar que estos valores de umbrales máximos y mínimos están con relación a los porcentajes de los percentiles que en este caso fue 75% y 25% respectivamente. Y como último grupo se tiene la cantidad de reservas que existe para eventos extremos y las tasas de subidas y bajadas, tomando en cuenta que se está contando con eventos extremos que puedan pasar de un día a otro.
4. RESULTADOS
Como resultado del estudio se obtuvieron productos del proceso de rellenado, la combinación de productos precipitación, configuración del modelo hidrológico y análisis de índices con IHA.
4.1. Resultados proceso de rellenado
Una vez ejecutado el programa HEC4 y generados los resultados, se compararon los productos obtenidos con la base de datos utilizada en el Balance Hídrico Superficial Nacional, de los cuales se buscó las estaciones similares y se comparó los períodos de 2014-2016, que es lo último que se tiene registro de esta base. Como se ve en las Figuras 7 y 8, el coeficiente de correlación de precipitación diaria fue mayores al 80 % por lo que hay una gran similitud, por otro lado, la precipitación mensual tuvo una correlación del 100 %.

Fuente: Elaboración Propia.
Figura 7: Comparación entre productos rellenados de HEC4 y base de datos del Balance Hídrico Nacional escala temporal diaria.
4.2. Resultados del proceso de combinación
Para el análisis de resultados de los productos combinados se realizó una comparación entre los productos satelitales iniciales y sus productos finales, obteniendo como mejor producto a CH5M con mejores valores iniciales en su producto base CHIRPS y con buenas tendencias y coeficiente de Nash cercano a 1, por otro lado, el producto con mayor cambio fue GS5M, el cual de tener una sobreestimación considerable en su producto base GSMaP, se igualaron los valores de sus indicadores con relación al resto de productos. Por esta razón, se seleccionaron estos dos productos para realizar las iteraciones a escala temporal diaria como también para la modelación en HEC-HMS, como se ven en la Tabla 2.
Como se puede observar en la Figura 9, el cambio en el mapa precipitaciones promedio anuales del producto satelital CHIRPS en el que se estimaba tener mayor precipitación en la zona sur de la cuenca y posterior al proceso de combinado en el mapa de precipitaciones promedio anual del producto CH5M, se observó una acumulación de precipitación en la zona noreste de la cuenca, de igual forma a como se presentaba en el mapa de precipitaciones promedio anual generado por pluviómetros.

Fuente: Elaboración propia.
Figura 9: Mapas de precipitación promedio anual en mm del periodo de 2014-2020 de: A) Pluviómetros, B) CHIRPS y C) CH5D.
En la Figura 10 se pudo apreciar en el mapa de precipitación promedio anual del producto GSMaP una gran sobreestimación, viendo que los valores de precipitación mínimos en todo el mapa son mayores a los de los rangos de mapas generados con pluviómetros y posterior al proceso de combinación se observó una similitud entre los mapas generados con pluviómetros y el producto combinado GS5D.

Fuente: Elaboración propia.
Figura 10: Mapas de precipitación promedio anual en mm del periodo de 2014-2020 de: A) Pluviómetros, B) GSMaP y C) GS5D.
Posterior a este proceso, se seleccionó dos productos, que fueron GSMaP y CHIRPS y se realizó otro proceso de combinación, pero a escala diaria, en los que se obtuvo como mejor producto a GS5D, el cual como se ve en la Tabla 3, tuvo un Nash muy cercano a 1 y una similitud del 95% con los datos de estaciones pluviométricas.
Comparando ambos procesos de combinación a diferentes escalas temporales, se pudo observar en la Figura 11 una mayor delimitación en los datos de precipitación con el producto combinado GS5D con respeto al producto GS5M, esto debido a que el factor de ajuste en el proceso de combinación era variado para cada día a escala diaria y en el caso de la escala mensual se repetía el valor del mensual para todos los días de ese mes.
4.3. Resultados de la modelación en HEC-HMS
Como validación del modelo HEC-HMS configurado, se compararon los caudales simulados con los caudales medidos en Puente Cajón en el período de 2014, obteniendo valores de Nash mayores al 0.2 lo que se podría considerar como aceptable, además que tuvo una correlación del 70% con respecto a los caudales medidos en campo. Como se ve en la Tabla 4.
TABLA 4 INDICADORES DE COMPARACIÓN DE CAUDALES SIMULADOS Y OBSERVADOS
INDICADORES | |
---|---|
R2 | 0.5017 |
R | 0.7083 |
RMSE | 3.8756 |
MAE | 15.0203 |
NASH | 0.2515 |
Fuente: Elaboración Propia.
Como se puede observar en la Figura 12, el caudal simulado con el producto GS5D es muy similar al caudal simulado con los pluviómetros. Adicionalmente, se vio una correlación entre la variación de la precipitación y de los caudales.
Como resultados de esta modelación se pudo observar que los caudales simulados con el producto CH5D comparados con caudal simulado con los pluviómetros presentan similitud como se ve en la Figura 13.

Fuente: Elaboración propia.
Figura 12: Caudales simulados con datos de pluviómetros y GS5D y hietogramas de precipitación combinada GS5D a la salida.

Fuente: Elaboración propia.
Figura 13: Caudales simulados con datos de pluviómetros y CH5D y hietogramas de precipitación combinada CH5D.
Los coeficientes de determinación, correlación y Nash con GS5M y CH5M se muestran en la Tabla 5. El producto de GS5D mostró resultados muy favorables con un coeficiente de Nash en un rango de 0.6 a 0.8, similar a los reportado [19], como también coeficientes de correlación y determinación.
4.4. Resultados de índices de alteración hidrológica
Una vez aplicado IHA7 se obtuvieron resultados, que se aprecian en la Tabla 6, donde se observa que en las zonas de Tiquipaya y Sacaba se tiene bajo caudal, no obstante, podían tener eventos extremos sorpresivos, mientras que en las zonas de Punata, Cliza y cercanías a La Angostura el caudal fue más estable y proporcional en cuanto a su incremento en crecidas y bajadas de caudal.
TABLA 6 DATOS SOBRE EL ESTUDIO GENERADOS POR EL IHA 7.1
Período del análisis: 2014-2021 (7 años) | SACABA | PUNATA | ANGOSTURA | CLIZA | TIQUIPAYA | SALIDA RIO ROCHA |
---|---|---|---|---|---|---|
Factor de normalización | 1 | 1 | 1 | 1 | 1 | 1 |
Caudal medio anual | 0.12 | 3.91 | 5.94 | 0.90 | 0.02 | 8.02 |
Caudal medio no normalizado | 0.12 | 3.91 | 5.94 | 0.90 | 0.02 | 8.02 |
Coeficiente Var. Anual | 7.10 | 3.88 | 3.85 | 4.79 | 28.38 | 3.94 |
Predictibilidad del caudal | 0.90 | 0.74 | 0.74 | 0.8 | 1 | 0.74 |
Constancia/predictibilidad | 0.88 | 0.73 | 0.73 | 0.79 | 1 | 0.72 |
% de inundaciones en períodos de 60 días | 1 | 0.26 | 0.25 | 0.67 | 0 | 0.27 |
Período libre de inundaciones | 50 | 8 | 8 | 11 | 315 | 8 |
Fuente: Elaboración Propia.
En Tiquipaya, como se puede apreciar en la Figura 14, se presentan caudales con bajas magnitudes; pero ante cualquier incremento de la precipitación este incrementa el caudal. Esto se debe a las pendientes inclinadas que existe en las microcuencas del Parque Nacional Tunari, donde precipitaciones intensas generaron eventos extremos con daños registrados, tal como señala el evento en febrero de 2018 el cual coincide con lo ocurrido en el diario Opinión [20], en el que se reportaron desbordes que generaron daño a varias viviendas debido al incremento del agua en los ríos de la zona.
En el caso de la zona del embalse de La Angostura se pudo apreciar en la Figura 15, existieron temporadas de lluvia que estuvieron por debajo del promedio, además que generaron diversos problemas de sequías, como se indicó en el periódico Los Tiempos en Diciembre de 2017 [21], Enero de 2018 [22] y Marzo de 2019 [23], en los que coincidieron que los niveles de la presa estaban bajos pese a las lluvias que se dieron en los días anteriores.
Se analizaron de igual forma los caudales en la salida del río Rocha para poder proponer un caudal ecológico que ayude a la zona. Como se observa en la Tabla 7, se debería mantener un 25% del caudal debido a que en el parámetro mínimo del 10% el caudal ecológico no sería recordable.
TABLA 7 PARÁMETROS PORCENTUALES DEL ANÁLISIS DEL CAUDAL ECOLÓGICO EN IHA7
MES | 10% | 25% | 50% | 75% | 90% | (75-25) /50 |
---|---|---|---|---|---|---|
Octubre | 0 | 0 | 0 | 0.020 | 0.091 | 0 |
Noviembre | 0 | 0 | 0.071 | 0.310 | 1.929 | 4.364 |
Diciembre | 0 | 0.030 | 2.021 | 4.915 | 8.766 | 2.417 |
Enero | 0 | 0.772 | 3.681 | 18.250 | 37.570 | 4.747 |
Febrero | 0 | 1.494 | 7.515 | 19.490 | 38.540 | 2.395 |
Marzo | 0 | 0 | 0.365 | 0.798 | 6.462 | 2.189 |
Abril | 0 | 0 | 0 | 0 | 0.039 | 0 |
Mayo | 0 | 0 | 0 | 0 | 0 | 0 |
Junio | 0 | 0 | 0 | 0 | 0 | 0 |
Julio | 0 | 0 | 0 | 0 | 0 | 0 |
Agosto | 0 | 0 | 0 | 0 | 0 | 0 |
Septiembre | 0 | 0 | 0 | 0 | 0 | 0 |
Fuente: Software IHA7.
5. CONCLUSIONES
En este estudio se han generado productos combinados de precipitación en la cuenca del rio Rocha usando pluviómetros locales con datos de sensores de precipitación a bordo de satélites de los productos CHIRPS, GSMAP, IMERG y PERSIANN. Se ha configurado un modelo en HEC-HMS en modo continuo con 59 unidades hidrológicas para obtener el caudal en puntos clave en la cuenca del río Rocha. Se logró evaluar el caudal simulado con las mediciones en el Puente Cajón. Posteriormente se utilizó el programa IHA7 para obtener parámetros de caudal ecológico.
Durante el proceso de rellenado para la preparación de datos, se obtuvieron resultados satisfactorios con el programa HEC4, ya que se compararon los resultados obtenidos con la base de datos generada para el Balance Hídrico Superficial Nacional y se logró observar un comportamiento similar entre los productos de precipitación estudiados. Con correlaciones a escala temporal diaria de más del 80 % y a escala temporal mensual del 100%.
En la combinación de productos de precipitación a escala mensual, luego de generar los cuatro productos combinados y compararlos mediante indicadores escogidos como el coeficiente de correlación y el coeficiente de Nash, se observa que el mejor producto es CH5M, ya que se obtuvo una eficiencia de Nash de 0.99. Durante el proceso de combinación a escala diaria, entre los productos CHIRPS y GSMAP se encontraron mejores resultados con el producto GS5D, ya que se obtuvo una eficiencia de Nash de 0.83 y una similitud del 95%. Además, se observa una mayor similitud entre los productos combinados a escala temporal diaria con las estaciones pluviométricas debido a que el factor de ajuste es aplicado cada día del mes, a diferencia de los mensuales en los que se aplica el mismo factor de un mes a todos los días de ese mes.
Obtenidos los datos de precipitación, se procedió a realizar el modelo hidrológico de la cuenca en el programa HEC-HMS, donde se obtuvo valores de eficiencia de Nash de más de 0.2, lo cual en la teoría de los criterios de los coeficientes de Nash se lo considera aceptable y una correlación del 70% con respecto a los datos medidos en Puente Cajón. En la modelación hidrológica de los productos combinados GS5D y CH5D, se obtuvo como mejor producto para simular a GS5D con una eficiencia de 0.73 y una correlación de 86 % al generado con datos de precipitación de estaciones terrestres.
Por último, se utilizó el programa IHA7 empleando los caudales simulados para el análisis de caudal ecológico en la cuenca. Empleando los eventos extremos se pudo determinar los percentiles de caudal que puedan afectar a la población, esto mediante la comparación de eventos reales ocurridos. Los productos de precipitación diaria permiten un análisis de eventos extremos en las microcuencas del Parque Tunari, las cuales son propensas a inundaciones y deslizamientos. Por otro lado, se observó que el caudal ecológico recomendado a la salida de la cuenca del río Rocha sería un 25 % del caudal promedio anual en el lugar, esto debido a que en el 10 % los caudales tienden a 0 por lo que no se mantendría el caudal ecológico. Se recomienda utilizar el esquema propuesto para evaluar otras cuencas en condiciones similares.