Assessing Near Real-Time Satellite Precipitation Products for Flood Simulations at Sub-Daily Scales in a Sparsely Gauged Watershed in Peruvian Andes

This study investigates the applicability of Satellite Precipitation Products (SPPs) in near real-time for the simulation of sub-daily runoff in the Vilcanota River basin, located in the southeastern Andes of Peru. The data from rain gauge stations are used to evaluate the quality of Integrated Multi-satellite Retrievals for GPM–Early (IMERG-E), Global Satellite Mapping of Precipitation–Near Real-Time (GSMaP-NRT), Climate Prediction Center Morphing Method (CMORPH), and HydroEstimator (HE) at the pixel-station level; and these SPPs are used as meteorological inputs for the hourly hydrological modeling. The GR4H model is calibrated with the hydrometric station of the longest record, and model simulations are also verified at one station upstream and two stations downstream of the calibration point. Comparing the sub-daily precipitation data observed, the results show that the IMERG-E product generally presents higher quality, followed by GSMaP-NRT, CMORPH, and HE. Although the SPPs present positive and negative biases, ranging from mild to moderate, they do represent the diurnal and seasonal variability of the hourly precipitation in the study area. In terms of the average of Kling-Gupta metric (KGE), the GR4H_GSMaP-NRT’ yielded the best representation of hourly discharges (0.686), followed by GR4H_IMERG-E’ (0.623), GR4H_Ensemble-Mean (0.617) and GR4H_CMORPH’ (0.606), and GR4H_HE’ (0.516). Finally, the SPPs showed a high potential for monitoring floods in the Vilcanota basin in near real-time at the operational level. The results obtained in this research are very useful for implementing flood early warning systems in the Vilcanota basin and will allow the monitoring and short-term hydrological forecasting of floods by the Peruvian National Weather and Hydrological Service.


Introduction
Floods are one of the most frequent natural disasters in Peru [1], negatively impacting the country at a socioeconomic level, affecting the physical and mental health of the population and causing severe damage to public and private property each year [2]. For example, heavy precipitation in January 2010 caused a great impact on the Andean region of Cusco. The Peruvian National Institute of Civil Defense (INDECI) reported that approximately 12,000 homes were affected. Total economic losses were estimated at USD 220 million, with 55% corresponding to infrastructure, 35% to health, and 10% to agriculture and tourism [3]. Despite the efforts made in recent decades, floods remain a recurring problem and tend to be increasingly frequent in the context of a changing climate. Therefore, flood monitoring and forecasting are key to reducing risk, improving decision-making and population response capacity, and mitigating its negative impacts [4].
Precipitation is the most common meteorological forcing of a hydrological model [5,6]. Precipitation data typically come from rain gauge stations; however, in mountainous and

Study Area
The basin of the Vilcanota River delimited at Intihuatana km105 stream gauge station ( Figure 1) lies in the Cusco Region, southeast of Peru, and it is part of the Andean region of the Amazon basin [27]. It has a drainage area of 9594 km 2 and a mean annual discharge of 130 m 3 /s. Its relief is predominantly rugged, characteristic of the Cordillera Vilcanota, with a mean slope and elevation of 15.6 • and 4176 masl, respectively. Annual precipitation ranges from 800 to 1000 mm. The Vilcanota basin is located in a transition zone with Remote Sens. 2021, 13, 826 3 of 18 tropical, sub-, and extratropical climates, a short wet season (November to March), and a long dry period during the rest of the year [28]. The basin continually suffers from floods that severely harm people's health, private property, and result in large economic losses to the local community and the Cusco region [29]. Therefore, there is an urgent need to implement a monitoring and flood forecasting system in the basin in quasi-real time. However, the low density and short recording times of the current network of rain gauge stations in the basin makes this task difficult [28]. Therefore, the use of estimated hourly satellite precipitation products with fine spatiotemporal resolution is a promising alternative for short-term flood simulations.
Remote Sens. 2021, 13, x FOR PEER REVIEW 3 of 18 of 130 m 3 /s. Its relief is predominantly rugged, characteristic of the Cordillera Vilcanota, with a mean slope and elevation of 15.6° and 4176 masl, respectively. Annual precipitation ranges from 800 to 1000 mm. The Vilcanota basin is located in a transition zone with tropical, sub-, and extratropical climates, a short wet season (November to March), and a long dry period during the rest of the year [28]. The basin continually suffers from floods that severely harm people's health, private property, and result in large economic losses to the local community and the Cusco region [29]. Therefore, there is an urgent need to implement a monitoring and flood forecasting system in the basin in quasi-real time. However, the low density and short recording times of the current network of rain gauge stations in the basin makes this task difficult [28]. Therefore, the use of estimated hourly satellite precipitation products with fine spatiotemporal resolution is a promising alternative for short-term flood simulations. In this work, a representative study domain was defined for the extraction, processing, and bias correction of the hourly SPPs, between the range of latitudes south 12.9º-14.8º, and longitudes west 70.6º-72.7º. Additionally, 11 sub-basins and river sections were delineated for semi-distributed hydrological modeling using the GR4H model (Figure 1), considering the location of the existing hydrometric stations in the basin. The details of the main physiographic characteristics of each sub-basin are shown in Table 1.

Ground-Based data
Hourly data were collected from 15 rain gauge stations and four hydrometric stations from the automated network of the Peruvian National Weather and Hydrological Service (SENAMHI) for the period from 1 January 2016 to 31 March 2020. Table 2 summarizes the selected stations. The local time is five hours less than the Coordinated Universal Time . Of the total rain gauge stations in the study domain, six are within the basin, between 1774 masl (INH) and 3920 masl (SAC). Data coverage in the study period ranges In this work, a representative study domain was defined for the extraction, processing, and bias correction of the hourly SPPs, between the range of latitudes south 12.9º-14.8º, and longitudes west 70.6º-72.7º. Additionally, 11 sub-basins and river sections were delineated for semi-distributed hydrological modeling using the GR4H model (Figure 1), considering the location of the existing hydrometric stations in the basin. The details of the main physiographic characteristics of each sub-basin are shown in Table 1.

Ground-Based Data
Hourly data were collected from 15 rain gauge stations and four hydrometric stations from the automated network of the Peruvian National Weather and Hydrological Service (SENAMHI) for the period from 1 January 2016 to 31 March 2020. Table 2 summarizes the selected stations. The local time is five hours less than the Coordinated Universal Time (UTC-5). Of the total rain gauge stations in the study domain, six are within the basin, between 1774 masl (INH) and 3920 masl (SAC). Data coverage in the study period ranges between 63.8% (INH) and 99.5% (SIC) at rain gauge stations and between 12% (INT) and 36.5% (CHI) at hydrometric stations, except for the PIS station, which has the longest record of hourly streamflows in the basin (98.3%). The distribution of rain gauge and Remote Sens. 2021, 13, 826 4 of 18 hydrometric stations in the study domain is shown in Figure 1. There is a greater density of rain gauge stations towards the northwest of the study domain and a low density towards the southeast. Hydrometric stations are distributed heterogeneously throughout the entire basin, with the middle basin (PIS) having the longest data record. Therefore, of the four stations collected, only one was selected for calibrating the GR4H model, and the remaining three were selected for its validation.

Satellite Precipitation Products (SPPs)
This study seeks to evaluate the feasibility of implementing a near real-time flood monitoring system at the operational level so that only gridded products with latencies less than 12 h (relative to UTC-5) were selected. Data from the products IMERG-E [30], GSMaP-NRT [31], CMORPH [32], and HE [33] were downloaded and processed for the period from 1 January 2016, to 31 March 2020. The selected SPPs have hourly temporal resolutions and spatial resolutions that vary between 0.1º (~10 km) and 0.05º (~5 km). The details of the characteristics of the selected SPPs are presented in Table 3.

Statistical Evaluation of SPPs against Rain Gauges
Because a pixel-to-pixel comparison requires a high density of stations to reduce the uncertainty of the interpolations, added to the different resolutions of the SPPs, in this study, the SPPs were evaluated using a point-to-pixel approach. The pixels of each SPP where the rain gauge stations are located in the study domain were extracted. SPP performance was evaluated through the calculation of statistics ( Table 4) that include the relative error (RE), coefficient of correlation (R), root mean square error (RMSE), and mean absolute error (MAE) using the hourly data of all stations for each month of the year. Table 4. Statistical metrics and their corresponding equations used for evaluating the meteorological performance of SPPs.

Statistical Metric Unit Equation Optimal Value
Relative Error (RE) - Note: n, number of samples; G i , precipitation from rain gauges; S i , precipitation estimates from satellite products.
Likewise, the representation of the diurnal cycle of precipitation in the rainy season (November to March) was evaluated for each station point. For this, only the days with daily precipitation greater than 15 mm were selected, and the coefficient of correlation was calculated for each station point and SPPs. This will allow evaluating not only the magnitude of the SPPs but also their diurnal variability.

Bias Correction of SPPs
In this work, rain gauge observations and the SPPs was merged using a simple bias adjustment correction [34]. For this, in each time interval (t), the pixel corresponding to the station point is extracted, and the bias to the station point (Ybias) is calculated as: then Ybias is interpolated by inverse distance weighting (IDW), using a weighting factor of 2, to generate a bias map (Mbias). Finally, the corrected bias is added to the SPPs (SPP ), as shown below:

Model Description and Setup
The GR4H model [23], which is a variant of the GR4J model of [35] with an hourly time step, is used to simulate the runoff in each of the 11 sub-basins of the Vilcanota River ( Table 1). The structure of the GR4H model is shown in Figure 2. The contribution of precipitation after interception (Pn) is divided into direct runoff (Pn−Ps) and infiltration (Ps) in the soil moisture reservoir of maximum storage capacity equal to X1. The percolation (Perc) from the soil moisture reservoir is mixed with the direct runoff to form the total runoff (Pr). Then, the total runoff is partitioned into two routing processes using a rate defined by parameter X2. Ninety percent of the runoff is routed by the combination of a unit hydrograph (UH1) and a nonlinear reservoir of maximum capacity X3. The other 10% is routed by another unit hydrograph (UH2). The base period of UH2 (2X4) is double that of UH1 (X4).  The GR4H model [23], which is a variant of the GR4J model of [35] with an hourly ime step, is used to simulate the runoff in each of the 11 sub-basins of the Vilcanota River Table 1). The structure of the GR4H model is shown in Figure 2. The contribution of preipitation after interception (Pn) is divided into direct runoff (Pn−Ps) and infiltration (Ps) n the soil moisture reservoir of maximum storage capacity equal to X1. The percolation Perc) from the soil moisture reservoir is mixed with the direct runoff to form the total unoff (Pr). Then, the total runoff is partitioned into two routing processes using a rate efined by parameter X2. Ninety percent of the runoff is routed by the combination of a nit hydrograph (UH1) and a nonlinear reservoir of maximum capacity X3. The other 10% s routed by another unit hydrograph (UH2). The base period of UH2 (2X4) is double that f UH1 (X4). To implement a flood forecast model in each sub-basin for operational purposes in he SENAMHI and given that the GR4H model is an aggregate model, the Muskingumunge method [37] was used to transfer the volume of runoff generated in each sub-basin semi-distributed model). In this way, the hourly flow was simulated for 11 sub-basins nd river sections (Table 1). Additionally, the GR4H model requires the input of hourly recipitation (PHR) and potential evapotranspiration (ETPHR) data for each sub-basin. The To implement a flood forecast model in each sub-basin for operational purposes in the SENAMHI and given that the GR4H model is an aggregate model, the Muskingum-Cunge method [37] was used to transfer the volume of runoff generated in each sub-basin (semi-distributed model). In this way, the hourly flow was simulated for 11 sub-basins and river sections (Table 1). Additionally, the GR4H model requires the input of hourly precipitation (P HR ) and potential evapotranspiration (ETP HR ) data for each sub-basin. The mean P HR of each sub-basin was calculated for each time interval as the mean of the bias-corrected SPP grids, while the ETP HR was estimated by the Hargreaves-Samani method [38] from the disaggregation of the daily climatic data (1981-2016) of the mean air temperature from the Peruvian Interpolated data of SENAMHI's Climatological and hydrological Observations (PISCO) product. The PISCO product is a gridded database of precipitation and temperature that coverage throughout the Peruvian territory for the period 1981-2016. The temperature product is generated by applying geostatistical techniques to combining air temperature data from MODIS images and observations from weather stations nationwide and is available online: http://iridl.ldeo.columbia.edu/ SOURCES/.SENAMHI/.HSR/.PISCO/ (accessed on 7 April 2019). Hydrological modeling was performed using the "airGR" package [39] in R language.

Model Calibration, Validation, and Verification
The calibration of the model considers the hourly recording at PIS stations from 1 January 2016, to 31 March 2018. The first 300 values were considered as part of the warmup period to reduce the uncertainty associated with the initial conditions of the model. The parameters of the GR4H model (X1, X2, X3 and X4) were automatically calibrated using the Shuffled Complex Evolutionary algorithm [40], considering the Kling-Gupta efficiency criterion [41] as an objective function. The mean absolute relative error (MARE), percentage bias (PBIAS), and root mean square error (RMSE) were used to evaluate the model performance in terms of low flows representation, model bias, and high flows errors, respectively. The summary of the statistics used to evaluate the performance of the hydrological model is detailed in Table 5.
The validation consisted of evaluating the outputs of the GR4H model at the PIS station for the period from 1 April 2018, to 31 March 2020, using the parameters obtained in the calibration. The verification process consisted of evaluating the model performance at the remaining hydrometric stations. For this, the hourly flow records available between 1 August 2018, and 31 March 2020, at stations SAL (upstream of PIS) and CHI and INT (downstream of PIS) were used. Likewise, this study evaluates the performance of the GR4H model by using the ensemble mean of simulations generated by IMERG-E', GSMaP-NRT', CMORPH', and HE'. In this way, the aim is to evaluate the representativeness of the mean of the simulations commonly used in a flow monitoring system at the operational level. Table 5. Statistical metrics and their corresponding equations used for evaluating the hydrological performance of SPPs.

Statistical Metric Unit Equation Optimal Value
Kling-Gupta efficiency (KGE) - Note: n, number of samples; O i , observed streamflow; X i , simulated streamflow. Figure 3 shows the seasonal variation in the metrics calculated to evaluate the meteorological performance of the SPPs using all the stations within the study domain. Among all of them, the RE (Figure 3a) shows that the differences between the SPPs are more pronounced. Specifically, IMERG-E tends to underestimate hourly precipitation during much of the year, reaching underestimates of up to −53% during the rainy season (November to March). In contrast, GSMaP-NRT predominantly overestimates up to 42.5% for the same period. The CMORPH and HE products show moderate differences between the dry period (April to October), but they coincide in underestimations close to −20% in the rainy period. The R values (Figure 3b) are less than 0.30, as expected for an hourly time scale; however, fit improves during the wet period. In terms of RMSE and MAE (Figure 3c,d), the GSMaP-NRT product presents the greatest bias of up to 1.137 and 0.43 mm/h, respectively, between February and March.

Validation of SPP against Rain Gauges
Additionally, Figure 4 shows the correlation (R) map of the diurnal cycle of the observed precipitation versus the SPP estimates for daily precipitation above 15 mm/d. The spatial distribution of R shows good positive fits (R ≥ 0.8) at station points located predominantly in the southwest portion of the study domain. The IMERG-E product (Figure 4a Figure 3 shows the seasonal variation in the metrics calculated to evaluate the meteorological performance of the SPPs using all the stations within the study domain. Among all of them, the RE (Figure 3a) shows that the differences between the SPPs are more pronounced. Specifically, IMERG-E tends to underestimate hourly precipitation during much of the year, reaching underestimates of up to −53% during the rainy season (November to March). In contrast, GSMaP-NRT predominantly overestimates up to 42.5% for the same period. The CMORPH and HE products show moderate differences between the dry period (April to October), but they coincide in underestimations close to −20% in the rainy period. The R values (Figure 3b) are less than 0.30, as expected for an hourly time scale; however, fit improves during the wet period. In terms of RMSE and MAE (Figure 3c and 3d), the GSMaP-NRT product presents the greatest bias of up to 1.137 and 0.43 mm/h, respectively, between February and March.   for PIS station illustrate the high variability of the hourly streamflows in the calibration and validation periods, with the exception of December 2019 to March 2020, where the simulations overestimate the observed streamflows. This overestimation could be attributed to the uncertainty of the latest PIS station's rating curve as a result of abrupt changes in the river section due to recent floods that occurred in the basin on December 2019 and January 2020, however this does not reduce model performance in the total period.

Evaluation of SPP's Hydrological Performance Using GR4H Model
The GR4H model was calibrated and validated using bias-corrected SPPs (IMERG-E', GSMaP-NRT', CMORPH', and HE') to further evaluate its hydrological utility in quasi real-time. Figure 5 shows the hourly hydrographs simulated for the PIS station ( Figure 1 and Table 2) using the four different hourly precipitation forcing inputs for the period from 1 January 2016, to 31 March 2020, and the hydrograph corresponding to the mean of the set of simulated flows is also shown (Figure 5e,j). The figure illustrates the feasibility of using near real-time SPPs to generate sub-daily runoff in a basin in the Andean region of Peru with scarce hydrometeorological information. The simulated hydrographs for PIS station illustrate the high variability of the hourly streamflows in the calibration and validation periods, with the exception of December 2019 to March 2020, where the simulations overestimate the observed streamflows. This overestimation could be attributed to the uncertainty of the latest PIS station's rating curve as a result of abrupt changes in the river section due to recent floods that occurred in the basin on December 2019 and January 2020, however this does not reduce model performance in the total period.  Figure 7 shows the hourly, daily, and monthly dispersion between the discharges observed at the PIS station and the simulated discharges using the four SPPs with bias    shows that the bias-corrected GR4H model forced by SPPs can largely reproduce the hourly discharge at hydrometric stations upstream and downstream of the PIS station (calibration). The figure shows the hourly discharges observed for the period from 1 August 2018 to 31 March 2020 (with missing data), and the simulated hourly discharges with four SPPs in near real-time and the hydrograph corresponding to the mean of the simulation datasets.       During the verification process, the performance statistics of the GR4H model are improved (Figure 9) when using GSMaP-NRT' as its precipitation forcing data, followed by IMERG-E', CMORPH' and HE', with the exception of the MARE index, where the simulations with CMORPH' yield better values. On the other hand, comparing the simulations at the three hydrometric verification stations, the values of KGE (Figure 9a)  During the verification process, the performance statistics of the GR4H model are improved (Figure 9) when using GSMaP-NRT' as its precipitation forcing data, followed by IMERG-E', CMORPH' and HE', with the exception of the MARE index, where the sim-  (Figure 9d).
In Figure 9, the statistical metrics using the ensemble mean approach (black bars) shows better results than HE' simulations and it has a performance comparable to other SPPs. For instance, in the SAL station, the ensemble means approach has values of KGE (0.644) and PBIAS (−26.4%) higher than CMORPH' (KGE = 0.534 and PBIAS of −30.5%).

Discussion
In this study, Satellite Precipitation Products (SPP) of fine spatial (0.05º to 0.10º) and temporal (hourly) resolution are used, and a conceptual hydrological model is used to simulate sub-daily streamflows in an Andean basin with scarce hydrometeorological information located in the southeastern region of Peru.
The precipitation data of the SPP versions in near real-time and rain gauge observations in the study domain were compared (Figure 3). This study found that of the four SPPs evaluated, IMERG-E (0.10º) has the lowest bias in hourly precipitation during the rainy period (November to March), followed by CMORPH (0.08º), HE (0.05º), and GSMaP-NRT (0.10º). This verifies that a finer spatial resolution and shorter latency product does not necessarily guarantee a better representation of precipitation. It was also found that In Figure 9, the statistical metrics using the ensemble mean approach (black bars) shows better results than HE' simulations and it has a performance comparable to other SPPs. For instance, in the SAL station, the ensemble means approach has values of KGE (0.644) and PBIAS (−26.4%) higher than CMORPH' (KGE = 0.534 and PBIAS of −30.5%).

Discussion
In this study, Satellite Precipitation Products (SPP) of fine spatial (0.05º to 0.10º) and temporal (hourly) resolution are used, and a conceptual hydrological model is used to simulate sub-daily streamflows in an Andean basin with scarce hydrometeorological information located in the southeastern region of Peru.
The precipitation data of the SPP versions in near real-time and rain gauge observations in the study domain were compared (Figure 3). This study found that of the four SPPs evaluated, IMERG-E (0.10º) has the lowest bias in hourly precipitation during the rainy period (November to March), followed by CMORPH (0.08º), HE (0.05º), and GSMaP-NRT (0.10º). This verifies that a finer spatial resolution and shorter latency product does not necessarily guarantee a better representation of precipitation. It was also found that the CMORPH, HE, and IMERG-E products tend to underestimate hourly precipitation throughout the year, while GSMaP-NRT tends to overestimate it (Figure 3a). Likewise, it is shown that IMERG-E has the best capacity to represent the cycle of daily precipitation above 15 mm/d (Figure 4a), followed again by CMORPH, HE, and GSMaP-NRT. The limited number of rain gauges (seven stations) in a medium-sized Andean basin (9594 km 2 ) makes it difficult to verify and validate the spatial patterns of hourly precipitation and would increase the uncertainty in the bias-corrected SPP. Therefore, the correction applied in this work also considered the stations located outside the basin; however, in the near future, this work would incorporate evaluating the effect of different bias-correcting methods on reducing the uncertainty of the SPPs.
According to [43], in addition to the uncertainties associated with the meteorological forcing, the performance of the hydrological model can be affected by sources of uncertainty related to the calibration of parameters, the structure of the model, and the streamflow measurements (uncertainty in rating curves). This work found that running the GR4H model provides satisfactory streamflow simulations at three hydrometric stations (PIS, SAL, and CHI) in the Vilcanota basin, as shown in Figures 5 and 8.This good performance is due to reducing the uncertainty of the SPPs during the bias correction but can also be associated with the compensation effect in the calibration mentioned in [44]. The latter refers to the fact that model parameters could compensate SPP errors during the calibration, in order to reproduce the hydrograph observed, and losing its predictable capacity in other sub-basins. In our work, this effect is only observed at the INT station since at the SAL and CHI stations, located upstream and downstream of the PIS calibration point, respectively, as shown in Figure 9. It should be noted that although our results provide information on the suitability of the SPPs for sub-daily hydrological modeling, the hydrological performance of these products is based only on the conceptual GR4H semi-distributed model, so the results could be different from a more complex hydrological model and with a greater spatial distribution of the meteorological forcings.
Finally, the objective of this study is to investigate the hydrological utility of SPPs in near real-time for the continuous and operational monitoring of sub-daily streamflows in the Vilcanota basin, so SPPs with the greatest near real-time post-corrections are not evaluated, such as IMERG-L, IMERG-F, GSMaP-MVK, etc. However, these products can be analyzed in subsequent studies. Despite this, we consider that the present study demonstrates the feasibility of the use of precipitation products estimated by satellites in the simulation of sub-daily streamflows in a data scarce basin and provides valuable information for the monitoring of floods at the operational level, as shown in Figure 10, where the four SPPs are used in near real-time to map the streamflows in 11 sub-basins (Table 1) during flood events on 6-7 February 2020. Figure 10. Example of simulated discharges in 11 Vilcanota's river streams for the 22:00 and 02:00 hours from 6 and 7 February 2020, respectively; using four SPPs as meteorological forcing data of the GR4H model.

Conclusions
There are few studies in the Andean region of Peru that focus on the near real-time hydrometeorological evaluation of SPPs and at sub-daily scales. In this work, we evaluated four SPPs (IMERG-E, GSMaP-NRT, CMORPH, and HE) and their hydrological applications in the Vilcanota River basin for a period of just over three years (2016-2020) at the hourly scale. We applied statistical metrics to evaluate the consistency between the SPPs and rain gauge observations, and we applied the GR4H model to analyze the precipitation-runoff transformation for the purposes of monitoring and forecasting sub-daily streamflows.
The statistical evaluation at the pixel-point level shows that the IMERG-E product generally presents the highest quality, followed by GSMaP-NRT, CMORPH, and HE. Although the magnitude of precipitation has positive and negative biases during the rainy season, the products are able to represent the diurnal and seasonal variability of the hourly precipitation at most of the stations in the study domain, including the stations within the basin.
The SPPs are capable of simulating the hydrological regime in the Vilcanota basin. The GSMaP-NRT'-forced simulations yielded better performance statistics in the simulated hourly streamflow rates, followed by IMERG-E' and CMORPH'. HE' was least able to represent runoff in the basin despite having the shortest latency. Operationally, and considering the latency of the SPPs with the best hydrological performance, hydrological monitoring of floods is feasible using the GSMaP-NRT' and IMERG-E' products with a 5 h delay. Figure 10. Example of simulated discharges in 11 Vilcanota's river streams for the 22:00 and 02:00 hours from 6 and 7 February 2020, respectively; using four SPPs as meteorological forcing data of the GR4H model.

Conclusions
There are few studies in the Andean region of Peru that focus on the near real-time hydrometeorological evaluation of SPPs and at sub-daily scales. In this work, we evaluated four SPPs (IMERG-E, GSMaP-NRT, CMORPH, and HE) and their hydrological applications in the Vilcanota River basin for a period of just over three years (2016-2020) at the hourly scale. We applied statistical metrics to evaluate the consistency between the SPPs and rain gauge observations, and we applied the GR4H model to analyze the precipitation-runoff transformation for the purposes of monitoring and forecasting sub-daily streamflows.
The statistical evaluation at the pixel-point level shows that the IMERG-E product generally presents the highest quality, followed by GSMaP-NRT, CMORPH, and HE. Although the magnitude of precipitation has positive and negative biases during the rainy season, the products are able to represent the diurnal and seasonal variability of the hourly precipitation at most of the stations in the study domain, including the stations within the basin.
The SPPs are capable of simulating the hydrological regime in the Vilcanota basin. The GSMaP-NRT'-forced simulations yielded better performance statistics in the simulated hourly streamflow rates, followed by IMERG-E' and CMORPH'. HE' was least able to represent runoff in the basin despite having the shortest latency. Operationally, and considering the latency of the SPPs with the best hydrological performance, hydrological monitoring of floods is feasible using the GSMaP-NRT' and IMERG-E' products with a 5 h delay.
Near real-time SPPs have shown high potential for flood monitoring in basins in the absence of a wide network of rain gauges, which facilitates the implementation of a hydrological early warning system. However, a greater number of new rain gauge stations in the basin would provide more information for a better evaluation of the ability of SPPs to capture the spatial variability of precipitation in the basin and the generation of runoff in sub-basins. Therefore, future studies should analyze the uncertainty of the SPP bias correction methods on the hydrological performance of the model.  Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not yet publicly available but are being organized into a larger dataset for publication in a public database.