RESUMEN El uso de datos de precipitación distribuida en modelos hidrológicos es importante para reflejar la variabilidad espacial de los procesos hidrológicos. En este estudio se ha implementado la herramienta HydroBID en la cuenca del río Guadalquivir utilizando tres productos de precipitación: i) estaciones pluviométricas; ii) GSMaP.v6_Gauge; iii) producto combinado GS. Este último se generó combinando el producto con base satelital y la red local de pluviómetros, a nivel de subcuenca. También se implementó el modelo hidrológico con los productos de precipitación y se obtuvo el caudal diario en la red hídrica, mostrando una óptima correlación de 0,99 y eficiencia de 0,96 durante el periodo 2000 - 2016 a nivel subcuenca. Los resultados simulados con el producto combinado GS presentaron mejor similitud a los caudales observados. Posteriormente, se aplicó el método de reducción de escala K-NN al modelo climático global MIROC5 RCP8.5 con las variables de precipitación y temperatura, considerando dos horizontes futuros de análisis, cercano (2031-2060) y futuro lejano (2061-2090). Los resultados muestran que la temperatura presenta patrones y comportamientos similares a los observados, con valores promedios 17 y 18°C, y variaciones de 1.5 y 2.7°C para los escenarios futuro cercano y lejano, respectivamente. En cambio, la precipitación mostró diferentes comportamientos en los tres productos empleados. Los valores de las estaciones pluviométricas y el producto GS presentan una reducción drástica al promedio observado histórico de precipitación empleado, con diferencias de hasta 230 mm/año en promedio. El producto GSMaP presentó valores de precipitación más altos, aproximadamente 80 mm/año sobre el promedio, los resultados de modelación empleando estos datos de precipitación indican que los caudales simulados presentan valores más bajos con respecto a los observados. Los caudales observados presentan un promedio máximo de 40 m3/s en marzo y un mínimo promedio de 3 m3/s en octubre. Por otro lado, los caudales simulados bajo efecto de cambio climático presentan un caudal máximo promedio de 22 m3/s en marzo y caudal mínimo de 3 m3/s en octubre, indicando una reducción de caudal en la época de lluvias.

ABSTRACT The use of distributed precipitation data in hydrological models is important to reflect the spatial variability of hydrological processes. In this study, the HydroBID tool has been implemented in the Guadalquivir River basin using three precipitation products: i) rainfall stations; ii) GSMaP.v6_Gauge; iii) combined GS product. The latter was generated by combining the satellite-based product and the local network of rain gauges, at the sub-basin level. The hydrological model was also implemented with the precipitation products and the daily flow in the water network was obtained, showing an optimum correlation of 0.99 and efficiency of 0.96 during the period 2000 - 2016 at the sub-basin level. The results simulated with the combined GS product showed better similarity to the observed flows. Subsequently, the K-NN downscaling method was applied to the MIROC5 RCP8.5 global climate model with the precipitation and temperature variables, considering two future horizons of analysis, near (2031-2060) and far future (2061-2090). The results show that temperature shows patterns and behaviors similar to those observed, with average values of 17 and 18°C, and variations of 1.5 and 2.7°C for the near and far future scenarios, respectively. On the other hand, precipitation showed different behaviors in the three products used. The values of the rain gauges and the GS product show a drastic reduction to the observed historical average of precipitation used, with differences of up to 230 mm/year on average. The GSMaP product presented higher precipitation values, approximately 80 mm/year above the average, the modeling results using this precipitation data indicate that the simulated flows present lower values with respect to the observed ones. The observed flows show a maximum average of 40 m3/s in March and a minimum average of 3 m3/s in October. On the other hand, the simulated flows under the effect of climate change show an average maximum flow of 22 m3/s in March and a minimum flow of 3 m3/s in October, indicating a reduction in flow during the rainy season.

RESUMEN La cuenca del río Rocha es una de las cuencas más importantes del departamento de Cochabamba pues en ella alberga aproximadamente 76% de la población del departamento, El río Rocha es un ecosistema urbano pues se encuentra situado dentro de la ciudad. Debido a su accidentado relieve topográfico se genera una gran cantidad de caudal en la época de lluvia. Por otra parte, en la época seca este río presenta un pequeño caudal proveniente de las descargas líquidas, residuales domésticas y/o industriales. El objetivo de este trabajo fue obtener series continuas de caudales en 3 secciones del río Rocha a un nivel sub horario para el periodo 15-dic-2017 a 27-abr-2023 esto realizando la recopilación, procesamiento y donde fue necesario las correcciones necesarias a los niveles medidos para posteriormente convertirlos en caudales mediante curvas H-Q. Una vez obtenidos los caudales estos junto a la precipitación fueron usados para el cálculo de los coeficientes de escurrimiento de las 3 subcuencas de estudio (el Abra, puente Cajón y puente Pico de Loro) obteniendo valores de 4.5%, 5..4% y 19% respectivamente , también se presenta el resultado de un análisis de extremos empírico para un periodo de retorno de 5 años con caudales medidos cada 15 min obteniendo caudales de hasta 16 m3/seg para la estación El Abra, 52 m3/seg para la estación Puente Cajón y 64 m3/seg para la estación Pico de Loro.

ABSTRACT The Rocha River basin is one of the most important basins in the department of Cochabamba because it is home to approximately 76% of the department's population. The Rocha River is an urban ecosystem because it is located within the city. Due to its rugged topographic relief, it generates a large amount of water during the rainy season. On the other hand, during the dry season, this river has a small flow from liquid discharges, domestic and/or industrial wastewater. The objective of this work was to obtain continuous series of flow rates in 3 sections of the Rocha River at a sub-hourly level for the period 15-Dec-2017 to 27-Apr-2023 by collecting, processing and, where necessary, making the necessary corrections to the measured levels to later convert them into flow rates by means of H-Q curves. Once the flow rates were obtained, these together with the precipitation were used to calculate the runoff coefficients of the 3 sub-basins of study (El Abra, Cajón bridge and Pico de Loro bridge), obtaining values of 4.5%, 5.4% and 19% respectively. The results of an empirical extremes analysis for a return period of 5 years with flows measured every 15 min are also presented, obtaining flows of up to 16 m3/sec for the El Abra station, 52 m3/sec for the Puente Cajón station and 64 m3/sec for the Pico de Loro station.

RESUMEN Durante el año hidrológico comprendido entre agosto de 2022 y julio de 2023, se evidenció una notable escasez de agua en el departamento de Cochabamba y en varias regiones del país, lo que motivó la necesidad de evaluar el comportamiento de la precipitación y su tendencia. Para este propósito, se utilizaron series anuales de totales de precipitación con la mayor cantidad de datos posibles, provenientes de puntos de medición en el departamento de Cochabamba. El análisis realizado mostró que la precipitación está disminuyendo, principalmente en el sur del departamento, en las provincias de Mizque, Campero y el sur de Carrasco. Este fenómeno no fue detectado en la vertiente norte de la cordillera del Tunari, donde, por el contrario, se presentó un incremento en la cantidad de lluvia precipitada. A su vez, se detectó que la precipitación del año hidrológico 2022-2023 presenta un déficit de entre el 31% y el 56% de la precipitación anual de largo período. Con el propósito de incorporar en el estudio los efectos del El Niño y La Niña, el análisis de la variabilidad, la oscilación y la tendencia en los ciclos de estos fenómenos, y su teleconexión con la precipitación, revela que los eventos de El Niño podrían acentuarse en el futuro, lo cual incide negativamente en la disponibilidad de agua proveniente de la lluvia en nuestra región. Además, nos encontramos en la meseta de un ciclo de sequía que podría durar algunos años más. Por lo tanto, es indispensable continuar con el análisis de los datos meteorológicos para apoyar la toma de decisiones anticipadas, con el fin de reducir los impactos del déficit de precipitación que pueda suscitarse en el futuro.

ABSTRACT During the hydrological year from August 2022 to July 2023, a notable water shortage was observed in the department of Cochabamba and in several regions of the country, which prompted the need to assess the behavior of precipitation and its trend. For this purpose, annual series of total precipitation with the most data possible were used, sourced from measurement points in the department of Cochabamba. The analysis revealed that precipitation is decreasing, primarily in the southern part of the department, in the provinces of Mizque, Campero, and the southern part of Carrasco. This phenomenon was not observed in the northern slope of the Tunari mountain range, where, on the contrary, an increase in the amount of rainfall was recorded. Additionally, it was found that the precipitation for the 2022-2023 hydrological year presents a deficit ranging from 31% to 56% of the long-term annual precipitation. To incorporate the effects of El Niño and La Niña into the study, the analysis of variability, oscillation, and trends in these cycles and their teleconnection with precipitation reveals that El Niño events may become more pronounced in the future, which negatively impacts the availability of rainwater in our region. Furthermore, we are currently in the midst of a drought cycle that could last for several more years. Therefore, it is essential to continue analyzing meteorological data to support decision-making in advance, in order to reduce the impacts of potential precipitation deficits in the future.

RESUMEN El cálculo de los caudales máximos por medio de la generación estocástica de tormentas elípticas dentro de la cuenca Misicuni se realizó con los softwares especializados IT-LluviaNH y IT-Inundación del Banco Interamericano de Desarrollo desarrollados por la plataforma CAPRA, obteniendo las curvas PADF características y 300 escenarios de tormentas estocásticas evaluados en los modelos hidrológicos e hidrodinámicos resultando en hidrogramas a nivel de presa y de subcuencas para diferentes periodos de retorno. Se determinaron también caudales máximos por medio de metodologías convencionales relacionando las curvas IDF y tormentas de diseño. La comparación gráfica y numérica de los hidrogramas a la salida de la cuenca permitió verificar la variación entre las metodologías aplicadas y el impacto en el aliviadero de la presa a través del tránsito de avenidas sobre el embalse. Los caudales máximos estocásticos resultaron inferiores a los determinados por las metodologías convencionales. Esto se debe a la no uniformidad de las intensidades de precipitación sobre el área de la cuenca. Sin embargo, contemplando la simultaneidad de escenarios de evento extremo independientes en cada subcuenca (aplicando porcentajes muy bajos a la probabilidad de excedencia) se tienen caudales máximos similares entre las metodologías. Asimismo, se calculó la proyección de los caudales con periodo de retorno de 2, 5, 10, 20, 50 y 100 años a periodos de retorno de 1000 y 10000 años para ambas metodologías. El impacto de los caudales estocásticos máximos laminados no representa un riesgo para el aliviadero de la presa al ser menores al caudal de diseño.

ABSTRACT The calculation of peak flows by means of stochastic generation of elliptical storms within the Misicuni basin was performed with the specialized software IT-LluviaNH and IT-Flood of the Inter-American Development Bank developed by the CAPRA platform, obtaining the characteristic PADF curves and 300 stochastic storm scenarios evaluated in the hydrological and hydrodynamic models resulting in hydrographs at the dam and sub-basin levels for different return periods. Peak flows were also determined by conventional methodologies relating to IDF curves and design storms. The graphical and numerical comparison of the hydrographs at the basin outlet made it possible to verify the variation between the methodologies applied and the impact on the dam spillway through the transit of floods over the reservoir. The stochastic peak flows were lower than those determined by conventional methodologies. This is due to the non-uniformity of precipitation intensities over the basin area. However, considering the simultaneity of independent extreme event scenarios in each sub-basin (applying very low percentages to the probability of exceedance), the maximum flows are similar between methodologies. Likewise, the projection of flows with return periods of 2, 5, 10, 20, 20, 50 and 100 years to return periods of 1000 and 10000 years was calculated for both methodologies. The impact of the maximum stochastic flows does not represent a risk for the dam spillway since they are lower than the design flow.

RESUMEN Los embalses del Valle Central de Tarija y su sedimentación constante ponen en riesgo el suministro de agua y consiguientemente el desarrollo de la economía departamental. El objetivo del análisis realizado es plantear un nuevo paradigma conceptual de uso del agua embalsada a largo plazo, y la necesidad de analizar e implementar la gestión de los sedimentos en el Valle Central de Tarija. Se propone una metodología que permite caracterizar la dinámica sedimentaria de embalses con escasa información y que están ubicados en zona de montaña, considerando la inexistencia de información hidrométrica y sedimentológica. La metodología consiste en la aplicación de ecuaciones empíricas que permiten determinar la capacidad de retención del embalse, la densidad aparente de los sedimentos, la degradación de la cuenca, verificar la vida de diseño y estimar la vida útil del embalse para finalmente categorizar el problema generado por los sedimentos mediante un indicador de colmatación. Se emplearon drones para la generación de la topografía de las riberas y levantamiento batimétrico con el objetivo de actualizar la superficie topo batimétrica. La metodología disminuye los tiempos de trabajo de campo y simplifica la logística requerida, generando productos de buena precisión y mayor detalle, y de bajo costo operacional para el monitoreo de la sedimentación. A través de la evaluación de tales mediciones en combinación con otros datos correlacionados se ha podido formular una ecuación de producción de sedimentos regional, destacando la necesidad de replicar la experiencia en otras cuencas de Bolivia.

ABSTRACT The reservoirs of the Central Valley of Tarija and their constant sedimentation pose a risk to the water supply and, consequently, the development of the regional economy. The goal of the analysis is to propose a new conceptual paradigm for the long-term use of stored water, as well as the need to analyze and implement sediment management in the Central Valley of Tarija. A methodology is proposed that allows for the characterization of the sedimentary dynamics of reservoirs with limited information and located in mountainous areas, considering the lack of hydrometric and sedimentological data. This methodology involves the application of empirical equations to determine the retention capacity of the reservoir, the apparent density of the sediments, the degradation of the watershed, verify the design life, and estimate the useful life of the reservoir, ultimately categorizing the sediment-related problem using a siltation indicator. Additionally, drones were employed to generate the topography of the banks and perform bathymetric surveys with the goal of updating the topo bathymetric surface. This methodology reduces fieldwork time and simplifies the required logistics, generating products with good precision and greater detail, at a low operational cost for sedimentation monitoring. Through the evaluation of these measurements, combined with other correlated data, a regional sediment production equation has been formulated, highlighting the need to replicate this experience in other watersheds in Bolivia.

RESUMEN El Laboratorio de Hidráulica de la Universidad Mayor de San Simón en Bolivia (LHUMSS) ha iniciado el estudio de la actividad morfológica del río Pilcomayo mediante la ejecución de trabajos de modelación física y numérica. La modelación numérica se ejecuta mediante la aplicación del modelo computacional Delft3D. En este proceso, el modelo computacional ha sido calibrado para la reproducción de las variables hidrodinámicas. El río Pilcomayo es un río morfológicamente activo, con cambios constantes en la forma del lecho debido al flujo y a las elevadas tasas de transporte de sedimentos en suspensión. El presente trabajo de investigación detalla la continuación en los trabajos de modelación numérica sobre el río. Se presenta la calibración sedimentológica, donde se reproduce el transporte de sedimentos observado en la Estación Villamontes, calibración refrendada con el respectivo análisis estadístico para la cuantificación de los errores residuales entre los valores observados y calculados por el modelo obteniéndose resultados calificados como “muy buenos” y “buenos” según la literatura existente. A partir de la calibración se ha efectuado la modelación de escenarios morfológicos de “largo plazo” y “corto plazo”, donde los resultados muestran la actividad morfológica de la zona de estudio.

ABSTRACT The Hydraulic Laboratory of the Universidad Mayor de San Simón in Bolivia (LHUMSS) has begun studying the morphological activity of the Pilcomayo River through the execution of physical and numerical modeling work. Numerical modeling is carried out using the Delft3D computational model. In this process, the computational model has been calibrated for the reproduction of hydrodynamic variables. The Pilcomayo River is morphologically active, with constant changes in the bed shape due to flow and high rates of suspended sediment transport. The present research work details the continuation of numerical modeling efforts on the river. The sedimentological calibration is presented, where the sediment transport observed at the Villamontes Station is reproduced, with the calibration supported by the respective statistical analysis for quantifying the residual errors between the observed and calculated values by the model, obtaining results rated as "very good" and "good" according to existing literature. Based on the calibration, the modeling of "long-term" and "short-term" morphological scenarios has been carried out, where the results show the morphological activity of the study area.

RESUMEN El propósito del presente estudio es conocer la eficiencia hidráulica de las obras de captación tipo Coanda en tres ríos de montaña del subandino de Bolivia, dichas obras forman parte de la Central Hidroeléctrica San José 1 y están ubicadas en el departamento de Cochabamba, provincia Chapare del Municipio de Villa Tunari. Para este estudio se comparó los caudales captados aforados frente a los caudales captados simulados por el programa Numerical Model of Coanda Screen Performance. Se llevaron campañas de aforo en los tres ríos durante el periodo marzo del 2018 a febrero del 2022. El número de las campañas y su frecuencia fue limitado por el advenimiento de la Pandemia COVID 19. Para la medición de los caudales se empleó el método de vadeo y utilizando correntómetros, los cuales se aforaron en diferentes puntos de cada obra de captación. Los resultados de la presente investigación permiten conocer valores del performance global para la pantalla Coanda en el río Ronquito del 80%, para Ronco del 74% y Antara del 80%. En atención a ello, las pantallas Coanda se constituyen en una buena alternativa de obras de captación de caudales en ríos de montaña, más aún si ellas reducen la necesidad a la inclusión de obras adicionales como desarenadores.

ABSTRACT The purpose of this study is to know the hydraulic efficiency of the Coanda type intake works in three mountain rivers of the sub-Andean region of Bolivia, these works are part of the San José 1 Hydroelectric Power Plant and are located in the department of Cochabamba, Chapare province of the Municipality of Villa Tunari. For this study, the gauged captured flows were compared with the captured flows simulated by the Numerical Model of Coanda Screen Performance program. Gauging campaigns were carried out in the three rivers during the period March 2018 to February 2022. The number of campaigns and their frequency was limited by the advent of Pandemic COVID 19. For the measurement of flow rates, the wading method was used and using current meters, which were gauged at different points of each catchment. The results of the present investigation show global performance values of 80% for the Coanda screen in the Ronquito river, 74% for Ronco and 80% for Antara. In view of this, the Coanda screens are a good alternative for catchment works in mountain rivers, especially if they reduce the need to include additional works such as sand traps.

RESUMEN En el valle central de Cochabamba, el aumento de la población y la creciente demanda de agua para consumo humano, industrial y riego ha provocado un incremento en la cantidad de pozos perforados, elevando el riesgo de la contaminación de las aguas subterráneas. Por esta razón, se ha visto necesario identificar la vulnerabilidad a la contaminación de acuíferos en una zona productiva del valle central de Cochabamba, para el cual se ha utilizado la herramienta DRASTIC. El método considera la combinación de siete factores hidrogeológicos utilizando herramientas SIG para los cuales se ha alimentado con información hidrogeológica local y actualizada. Para obtener el índice de vulnerabilidad se han analizado cuatro escenarios en base a los pesos asignados a cada factor. Los resultados indican que los factores con mayor afectación son la profundidad del nivel estático, la geología del acuífero, el impacto en la zona vadosa, la recarga y el tipo de suelo. La zona norte presenta vulnerabilidad Media Alta al ser zona de recarga con pozos profundos, la zona central indica vulnerabilidad Alta a Muy Alta por la presencia de pozos surgentes y actividades agrícolas, y la zona sur presenta de forma general vulnerabilidad Baja a Media Baja al ser una zona de descarga, urbanizada, acuíferos confinados y niveles de agua profundos. Por tanto, se recomienda dar un mayor énfasis a la zona norte y zona central, validar con estudios hidrogeoquímicos, elaborar planes de gestión, al igual que emplear el esquema propuesto en otras zonas de gran interés hidrogeológico.

ABSTRACT In the central valley of Cochabamba, the increase in population and the growing demand for water for human consumption, industry and irrigation has led to an increase in the number of wells drilled, raising the risk of groundwater contamination. For this reason, it has been necessary to identify the vulnerability to aquifer contamination in a productive zone of the central valley of Cochabamba, for which the DRASTIC tool has been used. The method considers the combination of seven hydrogeological factors using GIS tools for which it has been fed with local and updated hydrogeological information. To obtain the vulnerability index, four scenarios were analyzed based on the weights assigned to each factor. The results indicate that the factors with the greatest impact are the depth of the static level, the geology of the aquifer, the impact on the vadose zone, recharge and soil type. The northern zone presents Medium High vulnerability because it is a recharge zone with deep wells, the central zone indicates High to Very High vulnerability due to the presence of emerging wells and agricultural activities, and the southern zone generally presents Low to Medium Low vulnerability because it is a discharge zone, urbanized, confined aquifers and deep-water levels. Therefore, it is recommended that greater emphasis be given to the northern and central zones, validate with hydrogeochemical studies, develop management plans, and use the proposed scheme in other zones of great hydrogeological interest.

RESUMEN En la región metropolitana de Cochabamba, el crecimiento urbano hacia áreas de recarga acuífera ha aumentado el estrés hídrico en las aguas subterráneas, que abastecen aproximadamente el 65% del consumo total de agua. Por ello, se llevó a cabo un estudio hidrogeológico en una zona de abanicos aluviales en el Valle Central de Cochabamba. Se estableció una red de monitoreo de niveles freáticos en pozos de suministro de agua potable y riego en los municipios de Tiquipaya y Colcapirhua, realizando 26 mediciones entre abril y septiembre de 2021. Se desarrolló un modelo hidrogeológico simplificado con Visual MODFLOW Flex, calibrado con alta correlación (0.94) y un RMS de 17 m. En la validación, la correlación fue de 0.63, con un RMS normalizado del 78%. Se identificó la dirección principal del flujo subterráneo del norte a sur. Durante el monitoreo, se observó un descenso significativo del nivel freático, especialmente en los pozos OG-5 y OG-7, mientras que en OG-28 y OG-10 se registró un aumento de 3 a 5 metros desde 2019, indicando recarga asociada a los abanicos aluviales, principalmente el río Chijlawiri. Se recomienda emplear el modelo en estudios posteriores, mejorar la gestión de acuíferos y realizar un monitoreo continuo con datos mensuales para una mejor planificación hídrica.

ABSTRACT In the metropolitan region of Cochabamba, urban growth towards areas of aquifer recharge has increased water stress in groundwater, which supplies approximately 65% of total water consumption. Therefore, a hydrogeological study was carried out in an alluvial fan zone in the Central Valley of Cochabamba. A water table monitoring network was established in drinking water supply and irrigation wells in the municipalities of Tiquipaya and Colcapirhua, taking 26 measurements between April and September 2021. A simplified hydrogeological model was developed with Visual MODFLOW Flex, calibrated with a high correlation (0.94) and a 17 m RMS. In the validation, the correlation was 0.63, with a normalized RMS of 78%. The main direction of subsurface flow was identified from north to south. During monitoring, a significant decrease in the water table was observed, especially in wells OG-5 and OG-7, while in OG-28 and OG-10 an increase of 3 to 5 meters was recorded since 2019, indicating recharge associated with alluvial fans, mainly the Chijlawiri River. It is recommended to use the model in subsequent studies, improve aquifer management and carry out continuous monitoring with monthly data for better water planning.

RESUMEN La cuenca del Alto Paraguay, tiene ecosistemas de gran importancia debido a su rica biodiversidad y su papel crucial en la regulación de los recursos hídricos en la región. Sin embargo, en los últimos años se ha observado un creciente interés en comprender los cambios que están afectando a esta zona. Aunque existen investigaciones previas sobre el cambio climático a nivel global, se ha reconocido la necesidad de realizar un análisis detallado de la variabilidad y tendencia climática específicamente en la cuenca del Alto Paraguay. La cuenca abarca una superficie aproximada de 600.000 km² entre Bolivia, Brasil y Paraguay. Este estudio analiza la variabilidad y tendencia climática de las estaciones hidrometeorológicas del Río Alto Paraguay, con el fin de proporcionar información fundamental para la toma de decisiones en la planificación, diseño y operación del Puerto Busch. De esa manera, garantizar su sostenibilidad y resiliencia frente a los desafíos climáticos presentes y futuros, utilizando el método de Perturbación de Cuantiles. Este método requiere dos series temporales: una de las series se toma como serie de referencia o línea base, mientras que la otra se selecciona como una subserie. Las subseries son subperiodos seleccionados de estas series temporales completas y representan los periodos de interés con una longitud 12 años en el presente estudio. Con la aplicación de este método se obtendrá las anomalías climáticas para cada estación analizada. Los resultados sugirieron que existe una fuerte variabilidad temporal, donde se exploró la posible teleconexión con los patrones del Índice de Oscilación del Sur IOS en eventos de lluvia. Se encontraron patrones donde las estaciones están iniciando un ciclo de periodo seco.

ABSTRACT The Alto Paraguay Basin has ecosystems of great importance due to its rich biodiversity and its crucial role in regulating water resources in the region. However, in recent years, there has been a growing interest in understanding the changes affecting this area. Although there is previous research on climate change at a global level, the need for a detailed analysis of climate variability and trends specifically in the Alto Paraguay Basin has been recognized, which covers an area of approximately 600,000 km² across Bolivia, Brazil, and Paraguay. This study analyzes the variability and climate trend of the hydrometeorological stations of the Alto Paraguay River, aiming to provide essential information for decision-making in the planning, design, and operation of the Puerto Busch, ensuring its sustainability and resilience in the face of current and future climate challenges, using the Quantile Perturbation method. This method requires two time series: one series is taken as a reference or baseline series, while the other is selected as a subseries. The subseries are subperiods selected from these complete time series and represent the periods of interest, with a length of 12 years in this study. The application of this method will provide the climatic anomalies for each station analyzed. The results suggested a strong temporal variability, where the possible teleconnection with the Southern Oscillation Index (SOI) patterns in rainfall events was explored. Patterns were found where the stations are entering a dry period cycle.