As precipitation is a fundamental component of the global hydrological cycle that governs water resource distribution, the understanding of its temporal and spatial behavior is of great interest, and exact estimates of it are crucial in multiple lines of research. Meteorological data provide input for hydroclimatic models and predictions, which generally lack complete series. Many studies have addressed techniques to fill gaps in precipitation series at annual and monthly scales, but few have provided results at a daily scale due to the complexity of orographic characteristics and in some cases the non-linearity of precipitation. The objective of this study was to assess different methods of filling gaps in daily precipitation data using regression model (RM) and machine learning (ML) techniques. RM included linear regression (LRM) and multiple regression (MRM) algorithms, while ML included multiple regression algorithms (ML-MRM), K-nearest neighbors (ML-KNN), gradient boosting trees (ML-GBT), and random forest (ML-RF). This study covered the Malas, Omas, and Cañete River (MOC) watersheds, which are located on the Pacific Slope of central Peru, and a nineteen-year period of records (2001–2019). To assess model performance, different statistical metrics were applied. The results showed that the optimized machine learning (OML) models presented the least variability in estimation errors and the best approximation of the actual data from the study zone. In addition, this investigation shows that ML interprets and analyzes non-linear relationships between rain gauges at a daily scale and can be used as an efficient method of filling gaps in daily precipitation series.

The Penman-Monteith (PM) equation is considered the most accurate method for estimating reference evapotranspiration (ETo); however, its application requires a large amount of data that is not always available. This study aimed to regionalize the coefficients of the Hargreaves-Samani (HS) equation to estimate ETo in high-altitude areas, specifically the Peruvian Altiplano (PA). The methodology included (1) evaluation of the original HS equation, (2) calibration and validation of the empirical coefficient ($CH$) and empirical exponent ($EH$) at each weather station, and (3) regionalization of the calibrated coefficients using a multiple linear regression approach. The results showed that the original HS equation had NSE values ranging from −0.57 to 0.87, PBIAS from −18.60% to 12.70%, MAE from 0.16 to 0.65 mm/d, and RMSE from 0.20 to 0.67 mm/d. After calibrating $CH$ and $EH$, performance improved significantly, achieving validation values of NSE ranging from 0.67 to 0.94, PBIAS from −0.55% to 1.37%, MAE from 0.01 to 0.05 mm/d, and RMSE from 0.13 to 0.21 mm/d. Finally, the regionalization of 0.859 and 0.744, respectively. These results indicate that the HS equation, with calibrated and regionalized coefficients, is a viable alternative for estimating ETo in regions with limited meteorological data.