Gridded precipitation products on the Hindu Kush-Himalaya: Performance and accuracy of seven precipitation products

Climate change is expected to change precipitation and temperature patterns, which will impact the hydrological regime in Asia. Most river systems in the region originate from the Hindu Kush-Himalayas, and the altered precipitation patterns pose a threat to their sustainability, making it a major concern for planners and stakeholders. Obtaining accurate data on precipitation distribution is crucial for water accounting, which poses challenge. To address this, gridded precipitation products developed from satellite imagery and modeling techniques have become a viable alternative or addition to observed rainfall. However, the accuracy of these products in the region is uncertain. In this study, we aim to evaluate and compare the seven most commonly used precipitation products for the regions to address this gap. The study evaluated seven rainfall products, namely APHRODITE, TRMM, CHIRPS, PERSIANN-CDR, CMORPH, WFDEI, and GPCC by comparing daily, dekadal, and monthly rainfall data to 168 stations data in six countries and 11 river basins in the HKH region. The analysis used four continuous statistical indicators (Pearson correlation coefficient, Bias, Root Mean Square Error, and Nash–Sutcliffe Efficiency coefficient) and two categorical indicators (Probability of Detection and False Alarm Ratio). APHRODITE consistently performed well in several basins with high r values and low RMSE values, but had positive or negative bias values in different basins. CMORPH had the lowest positive bias value in the Ganga_Brahmaputra basin, while GPCC showed the largest r value and lowest RMSE


Introduction
The Himalayas have been called the 'water tower' of Asia as more than 3 billion people are estimated to depend on the combined flows of the major rivers originating in the region [1].The pressure on the area's water resources is immense [2] and is expected to increase significantly in the future due to a combination of climate change and socioeconomic developments [1,3].The combination of steep topography, rapidly draining soils and a strongly seasonal climate delivering 80% of the rainfall during the summer monsoon in the greater part of the Himalayas, renders the natural ecosystems of the region particularly fragile and thus vulnerable to disturbance [4].Studies in the Lesser Himalayas (3,700 to 4,500 m) by [4] suggested a significant increase in rainfall quantity and intensity during the monsoon and decreased in dry season rainfall.Such changes in precipitation patterns are a big threat to the sustainability of the river systems in the region.Owing to the socioeconomic and ecological importance of these river systems, the effects of climate change are of extreme concern among the planners and stakeholders in the region.
Accurate and detailed information on spatial variation in precipitation is critical in understanding the effects of climate change on Himalayan River systems.A robust climate change analysis and adaptation planning will not be possible without long-term hydro-meteorological monitoring.However, a lack of reliable and consistent data on precipitations severely limits our scientific understanding of the climate change impacts on these rivers [5].Moreover, long-term observed data are rare in the higher altitude area of the Himalayan region, largely due to the remoteness.On top of that, systematic bias is also present through the urbanization effect on meteorological observations, and wind effect observed precipitations [6].Therefore, a better understanding of the spatial and temporal extents of Himalayan precipitation is instrumental in water resources planning in the region.
Gridded precipitation products developed from satellite imagery sources have well served these gaps and stood as a strong alternative to the in-situ stations.These products come with varying spatial and temporal resolution along with individual strengths and weaknesses [7].The use of gridded precipitation data has been common for a long time for various applications such as hydrological, disaster management, and agricultural studies.With the increasing importance of these data for applications such as hydrological modeling, assessing climate change impacts, extreme climate induced hazards, and vegetation studies, it is important assess their quality as different gridded datasets may perform differently.Li et al. (2018) found that the distribution of simulated daily discharge values agreed well with observations while using in the WRF-Hydro (v3.5.1) modeling system in the Beas Basin [8].In another study, Immerzeel et.al. (2009) found that the discharge obtained from the TRMM 3B43 matched with the observed discharge in the Himalayan river basins [9].However it is noted that the snowmelt constitutes up to 50% of the total annual discharge in the Indus catchments, (far western) and ~25% in the Tsangpo (far eastern) catchment [10], therefore the performance depends on various factors such as other input data sources like temperature, soil moisture, land cover etc. as well as the model structures and parameters.The accumulated precipitation was found in overestimation in Southern China [11] while GPM-IERG products perform satisfactorily to catch the flash flood in Yunnan China even though some products over-estimation the precipitation [12].Therefore, such products perform differently in a different region of the Himalayas and therefore separate evaluation is important.
Those products are validated for the different locations of the world however, their general consistencies are not ensured all over the world as the accuracy varies from place to place, and product by product.Same product does not perform equally all over the world [13][14][15][16][17].
Therefore, in this study we aim to evaluate the most common precipitation products in the HKH regions.This study considered the seven most commonly used precipitation products based on the literature review.The product includes Asian Precipitation Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APHRODITE) 1801_R1 [26,27], A Climate Hazard Group InfraRed Precipitation with Station Data (CHIRPS) [27,28], The climate prediction center morphing method (CMORPH) [29], The Precipitation Estimation from Remotely Sensed Information using the Artificial Neural Networks Climate Data Record (PERSIANN-CDR) [27,30].The WFDEI meteorological forcing data set has been produced using the WATCH Forcing Data (WFD) [31,32], Tropical Rainfall Measuring Mission (TRMM 3B42) [29,33,34] and Global Precipitation Climatology Centre (GPCC) [30,35,36].These products were identified as the most commonly used products in the region.We compared the identified gridded precipitation with the ground measurement obtained from the 168 stations scattered throughout the HKH region.Six statistical measures, with four continuous and two categorical statics, were considered to compare with three different time accumulation, daily, 10-day, and 30-day.
This study distinguishes itself through its notable features, including a broader spatial coverage, a comprehensive evaluation of multiple precipitation products, an analysis considering three-time accumulation durations, and localized assessments at the national and basin levels.These contributions play a significant role in advancing our understanding of gridded precipitation patterns and their reliability within the HKH region.

Study area
HKH is the largest continuous mountainous range extending over 3500km from Afghanistan in the west to Myanmar in the east and hosts all 14 mountains above 8000m heigh (Fig 1).Ten major river systems of Asia-Amu Darya, Tarim, Indus, Ganges, Brahmaputra, Irrawaddy, Yangtze, Yellow, Mekong, and Salween are originated from the higher mountainous region of HKH.The total population of HKH is more than 240 million covering more than 4.2 million km 2 [37].
Monsoon precipitation contributes higher in the Siwalik and Pir-Panjal range of lower Himalayas and continuously decreases towards the northwards into the High Himalaya.The western part of the Karakoram Himalaya is likely to have heavy snowfall during winter due to the western disturbances that emanate from upper tropospheric westerlies [38,39].Monsoon contributes above 80% of annual precipitation in the central and eastern Himalayas.Precipitation is highly variable across the Himalayans river basins with annual mean precipitation is estimated 435, 1094 and 2143mm in the Indus, the Ganges, and the Brahmaputra river basins respectively [40,41].

Precipitation data sets
Ground observations.The Global Historical Climatology Network (GHCN) is an integrated database of daily climate summaries from land surface stations across the globe [43].A total of 317 stations are available in the Hindu Kush Himalayas region, however, only 168 stations have data for more than a year which are considered for analysis in this study.As in many other places, their distribution is not even and is mostly concentrated in low-altitude areas of the Indian Himalayas.The highest elevation of the measurement station is 4880m.The data are maintained and provided by the National Climatic Data Centre (NCDC) of the National Oceanic and Atmospheric Administration (NOAA).Although data are available for a long duration (less than one year to more than 175 years), we have selected the period from 2001 to 2016 in this study, as the gridded precipitation of all products and the observed precipitation overlap during this period.
Gridded precipitation products.APHRODITE is a 0.25˚spatial resolution, daily, gridded precipitation of monsoon Asia (MA) (60.0˚E-150.0˚E,15.0˚S-55.0˚N).The dataset is available from 1998 to 2015 and generated through the interpolation of the ground observation through the modified version of the distance weighting interpolation method [44].It includes the rainfall information from the Global Telecommunication System (GTS) network, and precompiled datasets from meteorological organizations of Asian countries [23].The series of monthly and daily rain-gauge data having a period of 5 years or more were considered while constructing this product.We obtained the dataset from the Research Institute for Humanity and Nature (RIHN) and the Meteorological Research Institute of Japan Meteorological Agency (MRI/JMA).
A CHIRPS is a near-global (50˚S to 50˚N) precipitation product of 0.05˚spatial resolution available from 1981 to present [45].It was developed to support the United States Agency for International Development Famine Early Warning Systems Network (FEWS NET) and in general agricultural drought monitoring.The major input data includes (i) expected annual rainfall sequence, computed through CHPClim (ii) thermal infrared images from the Climate Prediction Center (CPC) in every 30 minutes at 4km spatial resolution and B1 IR, a 3-hourly dataset of 8km spatial resolution from National Climatic Data Center (NCDC); (iii) 3B42 dataset from NASA obtained from the Tropical Rainfall Measuring Mission (TRMM 3B42), (iv) atmospheric model rainfall fields from climate forecast system version 2, NOAA, and (v) ground measurement from multiple sources.The data product is generated into two steps first includes the precipitation estimation from infrared images and in the second the estimated precipitation is blended with stations to produce the final CHIRPS product [46].
The CMORPH is a 30-minute global precipitation data that is derived from passive microwave images obtained from geostationary satellites [47].CMORPH uses motion vectors for precipitation estimation.Time-weighted linear precipitation is used to modify precipitation features such as shape and intensity.
The PERSIANN-CDR is available from 1983 to the near future from 60˚S to 60˚N with the spatial resolution of 0.25˚on the daily basis.It uses an artificial neural network that was trained using National Centers of Environmental Prediction (NCEP) stage IV hourly precipitation data [22].The input dataset includes the B1 infrared window (IRWIN) Channel of the Gridded Satellite Data (GridSat-B), Global Precipitation Climatology Project (GPCP) V2.2.The accuracy is improved by adaptive adjustment of the network parameters using rainfall estimates from the TMI.
The WFDEI meteorological forcing data set has been produced using the WATCH Forcing Data (WFD) by using of ERA-Interim reanalysis data [48].The ERA-Interim reanalysis data improves in precipitation and wind speed data while changes in the aerosol corrections improve downward shortwave fluxes.All precipitation corrections are conducted based on the CRU option.
GPCC is the product based on quality-controlled 67000 stations around the Globe that feature a record duration of 10 years or longer [49].Source from national meteorological and hydrological services, global and regional data collections as well as WMO GIS data.This is a daily product of 0.5˚× 0.5˚is an hourly dataset.The rainfall rate is based on GPCC bias-correction, under catch corrected measures.
TRMM 3B42 V7, a peri-global precipitation product, is a 0.25˚spatial resolution available in 50˚S to 50˚N from 1998 to the present [34].The estimated precipitation was generated through the TRMM Multi-satellite Precipitation Analysis (TMPA) system that utilizes multisatellite gauge data for bias correction.The daily precipitation was computed from a 3-hourly dataset.The major input includes TRMM Precipitation Radar (PR), TRMM Microwave Imager (TMI) the advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) on Aqua SSM/I and space sensor microwave imager (SSMIS), the AMSU-B and Microwave humidity Sounder, the IR data from an international constellation of Geosynchronous Earth Orbit (GEO) satellite and GPCP precipitation gauge analysis data.
We have summarized the details of all data products in Table 1 and visually represented the selected gridded precipitation products in Fig 2 .The figure shows that the precipitation distribution and intensity vary among the different products.We have highlighted some of the identified precipitation disparities between different products using red and purple circles in the figure.

Data analysis
The gridded products were compared to the ground station data on a point-to-pixel basis [50].The timing and intensity of rainfall event and the number of rainy days is highly spatially variables in the topographically heterogeneous region of HKH Mountain.Thus, we overlaid the stations over the raster images and extracted the timeseries dataset for the respective pixels on which intersected stations lie.We wrote a python script to get the timeseries dataset of the intersected pixel in the gridded product.The time series dataset in the selected pixels were subjected to analysis for three different temporal scales i.e., daily, dekadal (10-day), and monthly  (30-day).The 10-day and 30-day data were considered for analysis when at least 90% of daily observations were available.This resulted in 13,511 dekadal totals, and 3458 monthly totals over the 168 synoptic stations for the satellite products during the 16 years of study.Evaluation statistics.A pair-wise continuous statistics were used to evaluate the performance of the gridded products in estimating the rainfall quantity and categorical statistics were used to assess the rainfall detection capabilities.We considered 6 statistical indicators with 4 continuous and two categorical indicators for the analysis.The statistical measures considered in this study for the pairwise comparison are summarized in Table 2, it includes (i) the Pearson correlation coefficient (r-squared), this gives the degree of matching on estimating the gridded products against the observed values; its range is -1 to +1 with the +1 is the perfect results, (ii) bias, the degree of estimated values over or under-estimation was measured through bias, its range is -1 to 1, with 0 is the perfect value, (iii) RMSE, the commonly used measures of differences between two variables are measured through the root mean square error (RMSE)-it measures the average magnitude of the estimated errors; lower the RMSE values, greater the central tendencies and generally smaller extreme errors; (iv) Nash-Sutcliffe Efficiency coefficient (E) measures the pattern of predicted observed time series.It varies from -1 to 1, with 1 being the perfect value.(v) The probability of detection (POD) is the proportion of the correctly detected rainfall event to the total number of forecasted events.It gives the capacity of the correctly detected capacity of rainfall estimation; the range is 0 to 1.The perfect score is 1, (vi) False alarm ratio, a statistic for the categorical forecast, is the ratio of false alarm to the total number of forecasted events.It measures the reliability of the forecast.The perfect score is 0. The summary of the statistics is depicted in Table 2.

Ground observation availability and coverage
The availability of ground observation stations is very low in the Hindu Kush Himalayan regions in comparison to their counterparts in Europe and America.Table 3 shows the statistics of available ground observation stations globally based on GHCN [43].It is very high in North America, followed by Australia, Europe, South America, Asia, and Africa.There are 317 stations in the Hindu Kush Himalayas; however, 47% of the total stations were maintained for less than a year.Only the remaining 168 stations depicted in Fig 3 were subjected to analysis.

Continuous Statistics
Pearson correlation coefficient r ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi P n i¼1 p ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi P n i¼1 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi Most of the stations in the HKH were concentrated in the Southeast (Bangladesh) and the southern part of the Indian Himalayas.Tibetan plateau, the upstream part of the Yangtze River, has sparse stations but they are uniform.The central part of the Tibetan Plateau has very few or no stations.Most of the stations are located at an altitude less than 3000m whereas 57 stations are located at an altitude less than 1000m, 54 stations are in between 1000 and 2000m, 29 in between 2000 and 3000m while only 20 in between 3000 and 4000m and only eight stations lie above 4000m, where 4880m is the highest station.The details of the stations are presented in Fig 3.
Statistical analysis was performed for seven daily gridded precipitation products with observed precipitation (GHCN) collected from the ground stations at a point to pixel basis.The statistical analysis of seven gridded precipitation products overall six HKH countries and all eleven river basins are presented in the following subsections.
Results based on the country.This section provides an overview of the findings from a study that compared the performance of seven different gridded precipitation products in six countries in HKH regions.Fig 4 illustrates the statistical matrices for seven gridded precipitation products over six countries.Daily, dekadal and monthly dispersion are presented in the first, second and third columns respectively.The countrywide best-performing products are summarized in Table 4.
In the daily dataset, CHIRPS showed the highest r value (0.04) with less RMSE (13.62 mm) and 1.36 mm positive bias value among all seven gridded precipitation products in Afghanistan.In China, APHRODITE product exhibited the largest r value (0.51) and smallest RMSE value (4.87 mm) with 0.63 mm of positive bias value.Likewise in India, again APHRODITE product exhibited the largest r value (0.05) and less RMSE value (13.67 mm) with 0.05 mm of positive bias value.In Myanmar, GPCC displayed r value (0.11) and smallest RMSE value (12.05 mm) with 4.04 mm of negative bias value.In Nepal, APHRODITE product showed the highest r value (0.39) with less RMSE (11.05 mm) and 1.47 mm of negative bias value, whereas in Pakistan, GPCC displayed the largest r value (0.16) and less RMSE value (11.43 mm) with 1.46 mm of positive bias value.
The dekadal evaluation based on country data, CHIRPS showed the highest r value (0.16) with less RMSE (59.77 mm) and 15.72 mm of positive bias value among all seven gridded precipitation products in Afghanistan.In China, APHRODITE product exhibited the largest r value (0.76) and smallest RMSE value (18.87 mm) with 6.29 of positive bias value.Likewise in India, again APHRODITE product exhibited r value (0.15) and less RMSE value (73.62 mm) with the lowest 3.55 mm of negative bias value.In Nepal, APHRODITE product showed the highest r value (0.79) with less RMSE (48.52 mm) and 16.12 mm negative bias value, whereas in Pakistan, GPCC displayed the largest r value (0.46) and less RMSE value (31.14 mm) with 6.42 mm positive bias value.The monthly evaluation based on country data, CHIRPS showed the highest r value (0.29) with the lowest RMSE (121.5 mm) and 41.41 mm of positive bias value among all seven gridded precipitation products in Afghanistan.In China, APHRODITE product exhibited the largest r value (0.   5 summarizes the best-performing products in each basin. The daily evaluation based on basin data, APHRODITE showed the highest r value (0.13) with less RMSE (16.85 mm) and 0.35 mm of positive bias value among all seven gridded precipitation products in the Brahmaputra basin.In the Ganga basin, APHRODITE product exhibited the largest r value (0.09) and less RMSE value (14.99 mm) with 0.87 mm of positive bias value.Likewise in Ganga_Bhramaputra basin, CMORPH product exhibited the largest r value (0.12) and lowest RMSE value (12.4 mm) with 3.66 mm of positive bias value.In the Huang basin, APHRODITE displayed the highest r value (0.42) and smallest RMSE value (5.62 mm) with 0.98 mm positive bias value.In the Irrawaddy basin, again APHRODITE product showed the highest r value (0.19) with less RMSE (13.82 mm) and 0.18 mm of positive bias value.In the same way, in the Mekong basin, APHRODITE displayed the largest r value (0.73) and lowest RMSE value (3.71 mm) with no bias value.In Salween basin, APHRODITE product exhibited the largest r value (0.76) and lowest RMSE value (3.27 mm) with 0.7 mm of positive bias value.Likewise, in Sindha basin, APHRODITE product exhibited largest r value (0.05) and less RMSE value (11.6 mm) with 0.85 mm of positive bias value.In Tarim basin, APHRO-DITE displayed highest r value (0.4) and smallest RMSE value (2.75 mm) with 0.3 mm of positive bias value.In Tibet_interior basin, again APHRODITE product showed the highest r value (0.46) with lowest RMSE (5.62 mm) and 1.31 mm of positive bias value.In the same way, in Yangtze basin, again APHRODITE displayed largest r value (0.5) and lowest RMSE value (5.8 mm) with 0.67 mm of positive bias value.
The dekadal evaluation based on basin data, APHRODITE showed the highest r value (0.22) with less RMSE (89.65 mm) and 4.6 mm of negative bias value among all seven gridded precipitation products in Bhramaputra basin.In Ganga basin, APHRODITE product exhibited the largest r value (0.26) and less RMSE value (82.37 mm) with 8.32 mm of positive bias value, whereas, in Ganga_Bhramaputra basin, CMORPH product exhibited the smallest r value (0.29) but lowest RMSE value (58.75 mm) with 26.45 mm (low) of positive bias value.In the Huang basin, APHRODITE displayed the highest r value (0.62) and smallest RMSE value (23.3 mm) with 9.79 mm positive bias value.In the Irrawaddy basin, again APHRODITE product showed the high r value (0.45) with the lowest RMSE (68.7 mm) and 2.8 mm of positive bias value.In the same way, in the Mekong basin, APHRODITE displayed the largest r value (0.85) and lowest RMSE value (16.29 mm) with 0.35 mm of positive bias value.In Salween basin, APHRODITE product exhibited the largest r value (0.86) and lowest RMSE value (17.23 mm) with 7.69 mm of positive bias value.In Sindha basin, the GPCC product exhibited the largest r value (0.39) and lowest RMSE value (39.08 mm) with 8.54 mm positive bias value.In the Tarim basin, APHRODITE displayed the highest r value (0.69) and smallest RMSE value (9.35 mm) with 2.7 mm of positive bias value.In Tibet_interior basin, again APHRODITE product showed the highest r value (0.63) with the lowest RMSE (24.92 mm) and 12.6 mm of positive bias value.In the same way, in the Yangtze basin, again APHRODITE displayed the largest r value (0.81) and lowest RMSE value (19.79 mm) with 6.38 mm of positive bias value.
The monthly evaluation based on basin data, APHRODITE showed the high r value (0.29) with RMSE (201.51) and 27.62 (low) of negative bias value among all seven gridded precipitation products in the Brahmaputra basin.In Ganga basin, TRMM product exhibited the largest r value (0.31) and RMSE value (179.93) with 5.14 (low) positive bias value.Likewise in In summary, APHRODITE consistently showed the highest r values and lowest RMSE values across several basins, including the Brahmaputra, Ganga, Huang, Irrawaddy, Mekong, Salween, Tarim, Tibet_interior, and Yangtze basins.The product exhibited positive bias values in most basins, except for the dekadal evaluation in the Bhramaputra basin where it showed negative bias.CMORPH had the lowest positive bias value in the Ganga_Brahmaputra basin, and GPCC showed the largest r value and lowest RMSE value in the Sindha basin.

Discussion
The gridded precipitation datasets used in this study are typically based on observations from rain gauges, radar, and satellite remote sensing, which are interpolated onto a regular grid to create a spatially continuous dataset.APHRODITE has consistently demonstrated superior performance in most of the river basins and countries within the HKH region, when compared to the seven other products considered.Similar results were obtained in a comparative study with CHIRPS and PERSIANN-CDR in the Tibetan Plateau and its surroundings [26].In another study, APHRODITE outperformed the other five precipitation products (CPC-RFE, GSMaP, TRMM-3B42, TRMM-2B31) in the five small river basins in the Himalayas [51].However, all products were found to be quite reliable when evaluated on a monthly basis, as they demonstrated high correlation coefficients and low RMSE.A similar observation was made in an analysis of the Northwest Himalayan Region [36].Overall, APHRODITE consistently outperforms other products by better tracking above-threshold precipitation events, in which other products tend to underestimate the precipitation quantity [51].
The source of error in different products, obtained from various methodology, might be different, including the spatial and temporal resolution of the dataset, the interpolation method used, and the quality and density of the input observations.In the higher altitude regions of the Himalayas, the quality and density of input observations are reported to be relatively low and mostly very sparse [35,40], making it challenging to accurately represent precipitation patterns and amounts.Additionally, the high-altitude regions in the Himalayas are characterized by steep slopes, deep valleys, and complex atmospheric processes [52], which further complicates precipitation measurement and estimation.Furthermore, the rain shadow regions in the Himalayas, such as the regions to the north of the Great Himalayas, receive less precipitation due to the orographic effect [51], which can be difficult to accurately capture in gridded precipitation products.The higher temporal/spatial resolution datasets (hourly or few km scale) may have lower accuracy compared to datasets with lower temporal/spatial resolution (longer time of several km scales) due to the higher variability and uncertainty in precipitation measurements at this resolution.In the Himalayas, where precipitation is highly variable in space and time, accurately representing precipitation patterns and amounts can be particularly challenging.However, if gridded precipitation datasets are carefully validated and calibrated against observed streamflow data, they can still provide useful information for hydrological modeling in the Himalayas.
By considering a wider geographic area, evaluating multiple products, analyzing different time scales, and providing localized insights, our study contributes greatly to enhancing knowledge in this field.Ultimately, these findings strengthen our understanding of gridded precipitation in the HKH region.However, a large number of precipitation products, including ERA5 and ERA5-Land, Multi-Source Weighted-Ensemble Precipitation (MSWEP) v2.2, Modern-Era Retrospective Analysis for Research and Applications 2 (MERRA-2), and Integrated Multi-Satellite Retrievals (IMERG) for Global Precipitation Mission (GPM), near-realtime precipitation estimations based on machine learning algorithms, such as FY4QPE-MSA and PECA-FY4A etc. among others, are available for use in this region and have gained popularity in various applications.Therefore, we encourage future studies to assess the quality of these products.Furthermore, it would be worthwhile to study how different precipitation products perform in various topographical settings, as we have not taken this into account in this study.

Conclusion and recommendation
The performance of gridded precipitation products varies across countries and basins.After considering statistical parameters, it was found that CHIRPS is the better product for Afghanistan, APHRODITE for China, India, and Nepal, and GPCC for Pakistan, for all daily, dekadal, and monthly analyses.However, the performance of the products is not consistent across different basins.APHRODITE performs better than other products in most basins, except for the Ganga_Bhramaputra basin where CMORPH had favorable results in daily analysis.In dekadal analysis, APHRODITE outperformed other products in all basins except Ganga_Bhramaputra (CMORPH) and Sindha (GPCC).When it comes to monthly analysis, APHRODITE performed better in most basins except for Ganga (TRMM), Ganga_Bhramaputra (CMORPH), Irrawaddy (CHIRPS), and Sindha (CDR).In general, APHRODITE and TRMM performed best in all the indicators for all stations in the region, while CHIRPS performed the worst.
Due to the highly heterogeneous topography, gridded precipitation may not capture the timing and spatial variation within a short distance.As the gridded precipitation is an aggregation of a spatial extent of a single-pixel, it may moderate the results, and the timing could also have a significant impact on it.Given the highly heterogeneous topography and potential impact of timing and spatial variation on the accuracy of gridded precipitation products, it is recommended that daily estimates be used with caution in the HKH region.Nonetheless, some products' dekadal and monthly estimates can be used with confidence.It is advisable for researchers and professionals to evaluate the products for their specific study area and objectives before incorporating these gridded products into their analysis.

Fig 1 .
Fig 1. Map of the Hindu Kush Himalayan region, displaying the major river network and locations of precipitation stations that were taken into account for analysis.Source base layer: [42].https://doi.org/10.1371/journal.pwat.0000145.g001

Fig 2 .
Fig 2. Illustration of the spatial variation of precipitation as obtained from seven different gridded products; source base layer: [42].https://doi.org/10.1371/journal.pwat.0000145.g002 8) and less RMSE value (44.3 mm) with 18.49 mm positive bias value.Likewise in India, again APHRODITE product exhibited r value (0.22) and less RMSE value (165.04 mm) with the lowest 18.23 mm of negative bias value.In Nepal, APHRODITE product showed the highest r value (0.84) with less RMSE (140.61 mm) and 74.22 mm negative bias value, whereas in Pakistan, GPCC displayed the largest r value (0.54) and the lowest RMSE value (65.36 mm) with 17.75 mm positive bias value.In overall, the performance of different gridded precipitation products in estimation of precipitation in different countries varies.CHIRPS performed well in Afghanistan with high r values and low RMSE, but had positive bias values.APHRODITE performed well in China and India with high r values and less RMSE, but had positive or negative bias values in different evaluations.GPCC performed well in Myanmar and Pakistan with high r values and low RMSE, but had negative or positive bias values in different evaluations.APHRODITE performed consistently well in Nepal with high r values and less RMSE, but had negative bias values in some evaluations.Overall, different products had different strengths and weaknesses in predicting precipitation in different countries and types of evaluations.

Table 4 . Rankings of the top-performing data sources in various counties and different time intervals of aggregation.
Ganga_Bhramaputra basin, CMORPH product exhibited r value (0.45) and the lowest RMSE value (141.32) with 80.47 (low) positive bias value.In Huang basin, APHRODITE displayed the highest r value (0.66) and smallest RMSE value (53.29) with 29.08 of positive bias value.In Irrawaddy basin, CHIRPS product showed the highest r value (0.61) with the lowest RMSE (155.93) and 7.33 positive bias value.In Mekong basin, APHRODITE displayed largest r value (0.86) and RMSE value (41.33) with 2.13 of bias value.In Salween basin, APHRODITE

Fig 5. A summary of statistical metrics for seven gridded precipitation products and observed data, categorized by basin. https
://doi.org/10.1371/journal.pwat.0000145.g005product exhibited largest r value (0.88) and RMSE value (45.36) with 22.92 positive bias values.Likewise, in Sindha basin, CDR product exhibited the largest r value (0.22) and lowest RMSE value (98.5) with 6.63 of positive bias value.In Tarim basin, APHRODITE displayed the highest r value (0.74) and smallest RMSE value (22.78) with 8.03 of positive bias value.In Tibet_interior basin, again APHRODITE product showed the highest r value (0.68) with RMSE (58.32) and 37.25 of positive bias value.In the same way, in Yangtze basin, again APHRODITE displayed the largest r value (0.86) and RMSE value (44.74) with 18.64 of positive bias value.