Optimal design of rain gauge network in the Middle Yarra River catchment,Australia   总被引：1，自引：0，他引：1       下载免费PDF全文

Sajal Kumar Adhikary  Abdullah Gokhan Yilmaz  Nitin Muttil 《水文研究》2015,29(11):2582-2599

Rainfall data are a fundamental input for effective planning, designing and operating of water resources projects. A well‐designed rain gauge network is capable of providing accurate estimates of necessary areal average and/or point rainfall estimates at any desired ungauged location in a catchment. Increasing network density with additional rain gauge stations has been the main underlying criterion in the past to reduce error and uncertainty in rainfall estimates. However, installing and operation of additional stations in a network involves large cost and manpower. Hence, the objective of this study is to design an optimal rain gauge network in the Middle Yarra River catchment in Victoria, Australia. The optimal positioning of additional stations as well as optimally relocating of existing redundant stations using the kriging‐based geostatistical approach was undertaken in this study. Reduction of kriging error was considered as an indicator for optimal spatial positioning of the stations. Daily rainfall records of 1997 (an El Niño year) and 2010 (a La Niña year) were used for the analysis. Ordinary kriging was applied for rainfall data interpolation to estimate the kriging error for the network. The results indicate that significant reduction in the kriging error can be achieved by the optimal spatial positioning of the additional as well as redundant stations. Thus, the obtained optimal rain gauge network is expected to be appropriate for providing high quality rainfall estimates over the catchment. The concept proposed in this study for optimal rain gauge network design through combined use of additional and redundant stations together is equally applicable to any other catchment. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.  相似文献   

Rainfall network design using kriging and entropy   总被引：4，自引：0，他引：4  

Yen‐Chang Chen  Chiang Wei  Hui‐Chung Yeh 《水文研究》2008,22(3):340-346

The spatial distribution of rainfall is related to meteorological and topographical factors. An understanding of the weather and topography is required to select the locations of the rain gauge stations in the catchment to obtain the optimum information. In theory, a well‐designed rainfall network can accurately represent and provide the needed information of rainfall in the catchment. However, the available rainfall data are rarely adequate in the mountainous area of Taiwan. In order to provide enough rainfall data to assure the success of water projects, the rainfall network based on the existing rain gauge stations has to be redesigned. A method composed of kriging and entropy that can determine the optimum number and spatial distribution of rain gauge stations in catchments is proposed. Kriging as an interpolator, which performs linear averaging to reconstruct the rainfall over the catchment on the basis of the observed rainfall, is used to compute the spatial variations of rainfall. Thus, the rainfall data at the locations of the candidate rain gauge stations can be reconstructed. The information entropy reveals the rainfall information of the each rain gauge station in the catchment. By calculating the joint entropy and the transmitted information, the candidate rain gauge stations are prioritized. In addition, the saturation of rainfall information can be used to add or remove the rain gauge stations. Thus, the optimum spatial distribution and the minimum number of rain gauge stations in the network can be determined. The catchment of the Shimen Reservoir in Taiwan is used to illustrate the method. The result shows that only seven rain gauge stations are needed to provide the necessary information. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Rain‐gauge network evaluation and augmentation using geostatistics     

Ke‐Sheng Cheng  Yun‐Ching Lin  Jun‐Jih Liou 《水文研究》2008,22(14):2554-2564

Rain‐gauge networks are often used to provide estimates of area average rainfall or point rainfalls at ungauged locations. The level of accuracy a network can achieve depends on the total number and locations of gauges in the network. A geostatistical approach for evaluation and augmentation of an existing rain‐gauge network is proposed in this study. Through variogram analysis, hourly rainfalls are shown to have higher spatial variability than annual rainfalls, with hourly Mei‐Yu rainfalls having the highest spatial variability. A criterion using ordinary kriging variance is proposed to assess the accuracy of rainfall estimation using the acceptance probability defined as the probability that estimation error falls within a desired range. Based on the criterion, the percentage of the total area with acceptable accuracy A_p under certain network configuration can be calculated. A sequential algorithm is also proposed to prioritize rain‐gauges of the existing network, identify the base network, and relocate non‐base gauges. Percentage of the total area with acceptable accuracy is mostly contributed by the base network. In contrast, non‐base gauges provide little contribution to A_p and are subject to removal or relocation. Using a case study in northern Taiwan, the proposed approach demonstrates that the identified base network which comprises of approximately two‐thirds of the total rain‐gauges can achieve almost the same level of performance (expressed in terms of percentage of the total area with acceptable accuracy) as the complete network for hourly Mei‐Yu rainfall estimation. The percentage of area with acceptable accuracy can be raised from 56% to 88% using an augmented network. A threshold value for the percentage of area with acceptable accuracy is also recommended to help determine the number of non‐base gauges which need to be relocated. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods     

Haifa Feki  Mohamed Slimani  Christophe Cudennec 《水文科学杂志》2013,58(7):1294-1314

Abstract

This paper compares the performance of three geostatistical algorithms, which integrate elevation as an auxiliary variable: kriging with external drift (KED); kriging combined with regression, called regression kriging (RK) or kriging after detrending; and co-kriging (CK). These three methods differ by the way by in which the secondary information is introduced into the prediction procedure. They are applied to improve the prediction of the monthly average rainfall observations measured at 106 climatic stations in Tunisia over an area of 164 150 km² using the elevation as the auxiliary variable. The experimental sample semivariograms, residual semivariograms and cross-variograms are constructed and fitted to estimate the rainfall levels and the estimation variance at the nodes of a square grid of 20 km?×?20 km resolution and to develop corresponding contour maps. Contour diagrams for KED and RK were similar and exhibited a pattern corresponding more closely to local topographic features when (a) the network is sparse and (b) the rainfall–elevation correlation is poor, while CK showed a smooth zonal pattern. Smaller prediction variances are obtained for the RK algorithm. The cross-validation showed that the RMSE obtained for CK gave better results than for KED or RK.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Feki, H., Slimani, M., and Cudennec, C., 2012. Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods. Hydrological Sciences Journal, 57 (7), 1294–1314.  相似文献   

Temporal and spatial rainfall analysis across a humid tropical catchment     

P. Campling  A. Gobin  J. Feyen 《水文研究》2001,15(3):359-375

Temporal and spatial rainfall patterns were analysed to describe the distribution of daily rainfall across a medium‐sized (379km²) tropical catchment. Investigations were carried out to assess whether a climatological variogram model was appropriate for mapping rainfall taking into consideration the changing rainfall characteristics through the wet season. Exploratory, frequency and moving average analyses of 30 years' daily precipitation data were used to describe the reliability and structure of the rainfall regime. Four phases in the wet season were distinguished, with the peak period (mid‐August to mid‐September) representing the wettest period. A low‐cost rain gauge network of 36 plastic gauges with overflow reservoirs was installed and monitored to obtain spatially distributed rainfall data. Geostatistical techniques were used to develop global and wet season phase climatological variograms. The unscaled climatological variograms were cross‐validated and compared using a range of rainfall events. Ordinary Kriging was used as the interpolation method. The global climatological variogram performed better, and was used to optimize the number and location of rain gauges in the network. The research showed that although distinct wet season phases could be established based on the temporal analysis of daily rainfall characteristics, the interpolation of daily rainfall across a medium‐sized catchment based on spatial analysis was better served by using the global rather than the wet season phase climatological variogram model. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan          下载免费PDF全文

Meng‐Jung Tsai  Robert J. Abrahart  Nick J. Mount  Fi‐John Chang 《水文研究》2014,28(3):1055-1070

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Interpolation processes using multivariate geostatistics for mapping of climatological precipitation mean in the Sannio Mountains (southern Italy)     

Nazzareno Diodato  Michele Ceccarelli 《地球表面变化过程与地形》2005,30(3):259-268

The spatial variability of precipitation has often been a topic of research, since accurate modelling of precipitation is a crucial condition for obtaining reliable results in hydrology and geomorphology. In mountainous areas, the sparsity of the measurement networks makes an accurate and reliable spatialization of rainfall amounts at the local scale difficult. The purpose of this paper is to show how the use of a digital elevation model can improve interpolation processes at the subregional scale for mapping the mean annual and monthly precipitation from rainfall observations (40 years) recorded in a region of 1400 km² in southern Italy. Besides linear regression of precipitation against elevation, two methods of interpolation are applied: inverse squared distance and ordinary cokriging. Cross‐validation indicates that the inverse distance interpolation, which ignores the information on elevation, yields the largest prediction errors. Smaller prediction errors are produced by linear regression and ordinary cokriging. However, the results seem to favour the multivariate geostatistical method including auxiliary information (related to elevation). We conclude that ordinary cokriging is a very flexible and robust interpolation method because it can take into account several properties of the landscape; it should therefore be applicable in other mountainous regions, especially where precipitation is an important geomorphological factor. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

Evaluation and comparison of interpolated gauge rainfall data and gridded rainfall data in Florida,USA     

Meijing Zhang  Conrado de Leon 《水文科学杂志》2018,63(4):561-582

A variety of spatially continuous rainfall products are available but little evaluation of their accuracy has been published for areas with high spatial variability in rainfall. Five gridded rainfall products (PRISM, RTMA, and the interpolated Florida Automated Weather Network, FAWN, rainfall layers based on three interpolated methods) were assessed for Florida State. Point-to-pixel and pixel-to-pixel comparisons were performed to compare the five products. On average, the PRISM and RTMA products resulted in a better fit with the daily FAWN rainfall datasets, while FAWN-based interpolated products resulted in a better fit with the monthly FAWN rainfall datasets based on point-to-pixel analysis. Inverse distance weighting and ordinary kriging methods performed slightly better than the thin plate spline method in predicting daily rainfall. In general, monthly and seasonal rainfall amounts from PRISM and RTMA products were higher and lower, respectively, than reference rainfall amounts from FAWN gauge stations and FAWN-based interpolated products.  相似文献   

Spatial and temporal variability of the hydrological response in a small Mediterranean research catchment (Vallcebre,Eastern Pyrenees)     

J. Latron  M. Soler  P. Llorens  F. Gallart 《水文研究》2008,22(6):775-787

This paper analyses the spatial and temporal variability of the hydrological response in a small Mediterranean catchment (Cal Rodó). The first part of the analysis focuses on the rainfall–runoff relationship at seasonal and monthly scale, using an 8‐year data set. Then, using storm‐flow volume and coefficient, the temporal variability of the rainfall–runoff relationship and its relationship with several hydrological variables are analysed at the event scale from hydrographs observed over a 3‐year period. Finally, the spatial non‐linearity of the hydrological response is examined by comparing the Cal Rodó hydrological response with the Can Vila sub‐catchment response at the event scale. Results show that, on a seasonal and monthly scale, there is no simple relationship between rainfall and runoff depths, and that evapotranspiration is a factor that introduced some non‐linearity in the rainfall–runoff relationship. The analysis of monthly values also reveals the existence of a threshold in the relationship between rainfall and runoff depths, denoting a more contrasted hydrological response than the one usually observed in humid catchments. At the event scale, the storm‐flow coefficient has a clear seasonal pattern with an alternance between a wet period, when the catchment is hydrologically responsive, and a dry summer period, when the catchment is much less reactive to any rainfall. The relationship between the storm‐flow coefficient and rainfall depth, rainfall maximum intensity and base‐flow shows that observed correlations are the same as those observed for humid conditions, even if correlation coefficients are notably lower. Comparison with the Can Vila sub‐catchment highlights the spatial heterogeneity of the rainfall‐runoff relationship at the small catchment scale. Although interpretation in terms of runoff processes remains delicate, heterogeneities between the two catchments seem to be related to changes in the ratio between infiltration excess and saturation processes in runoff formation. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Spatial prediction of river channel topography by kriging   总被引：2，自引：0，他引：2  

Carl J. Legleiter  Phaedon C. Kyriakidis 《地球表面变化过程与地形》2008,33(6):841-867

Topographic information is fundamental to geomorphic inquiry, and spatial prediction of bed elevation from irregular survey data is an important component of many reach‐scale studies. Kriging is a geostatistical technique for obtaining these predictions along with measures of their reliability, and this paper outlines a specialized framework intended for application to river channels. Our modular approach includes an algorithm for transforming the coordinates of data and prediction locations to a channel‐centered coordinate system, several different methods of representing the trend component of topographic variation and search strategies that incorporate geomorphic information to determine which survey data are used to make a prediction at a specific location. For example, a relationship between curvature and the lateral position of maximum depth can be used to include cross‐sectional asymmetry in a two‐dimensional trend surface model, and topographic breaklines can be used to restrict which data are retained in a local neighborhood around each prediction location. Using survey data from a restored gravel‐bed river, we demonstrate how transformation to the channel‐centered coordinate system facilitates interpretation of the variogram, a statistical model of reach‐scale spatial structure used in kriging, and how the choice of a trend model affects the variogram of the residuals from that trend. Similarly, we show how decomposing kriging predictions into their trend and residual components can yield useful information on channel morphology. Cross‐validation analyses involving different data configurations and kriging variants indicate that kriging is quite robust and that survey density is the primary control on the accuracy of bed elevation predictions. The root mean‐square error of these predictions is directly proportional to the spacing between surveyed cross‐sections, even in a reconfigured channel with a relatively simple morphology; sophisticated methods of spatial prediction are no substitute for field data. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

A non-parametric automatic blending methodology to estimate rainfall fields from rain gauge and radar data   总被引：1，自引：0，他引：1  

Carlos A. Velasco-Forero  Daniel Sempere-Torres  Eduardo F. Cassiraga  J. Jaime Gómez-Hernández 《Advances in water resources》2009

Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.  相似文献   

Effect of spatial variation characteristics on contouring of design storm depth     

Ke‐Sheng Cheng  Chiang Wei  Yen‐Ben Cheng  Hui‐Chung Yeh 《水文研究》2003,17(9):1755-1769

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

Geostatistical radar–raingauge merging: A novel method for the quantification of rain estimation accuracy     

《Advances in water resources》2014

Compared to other estimation techniques, one advantage of geostatistical techniques is that they provide an index of the estimation accuracy of the variable of interest with the kriging estimation standard deviation (ESD). In the context of radar–raingauge quantitative precipitation estimation (QPE), we address in this article the question of how the kriging ESD can be transformed into a local spread of error by using the dependency of radar errors to the rain amount analyzed in previous work. The proposed approach is implemented for the most significant rain events observed in 2008 in the Cévennes-Vivarais region, France, by considering both the kriging with external drift (KED) and the ordinary kriging (OK) methods. A two-step procedure is implemented for estimating the rain estimation accuracy: (i) first kriging normalized ESDs are computed by using normalized variograms (sill equal to 1) to account for the observation system configuration and the spatial structure of the variable of interest (rainfall amount, residuals to the drift); (ii) based on the assumption of a linear relationship between the standard deviation and the mean of the variable of interest, a denormalization of the kriging ESDs is performed globally for a given rain event by using a cross-validation procedure. Despite the fact that the KED normalized ESDs are usually greater than the OK ones (due to an additional constraint in the kriging system and a weaker spatial structure of the residuals to the drift), the KED denormalized ESDs are generally smaller the OK ones, a result consistent with the better performance observed for the KED technique. The evolution of the mean and the standard deviation of the rainfall-scaled ESDs over a range of spatial (5–300 km²) and temporal (1–6 h) scales demonstrates that there is clear added value of the radar with respect to the raingauge network for the shortest scales, which are those of interest for flash-flood prediction in the considered region.  相似文献   

Continuous daily hydrograph simulation using duration curves of a precipitation index     

V. Yu. Smakhtin  B. Masse 《水文研究》2000,14(6):1083-1100

The paper describes a parsimonious approach for generating continuous daily stream‐flow time‐series from observed daily rainfall data in a catchment. The key characteristic in the method is a duration curve. It is used to convert the daily rainfall information from source rain gauges into a continuous daily hydrograph at the destination river site. For each source rain gauge a time‐series of rainfall related ‘current precipitation index’ is generated and its duration curve is established. The current precipitation index reflects the current catchment wetness and is defined as a continuous function of precipitation, which accumulates on rainy days and exponentially decays during the periods of no rainfall. The process of rainfall‐to‐runoff conversion is based on the assumption that daily current precipitation index values at rainfall site(s) in a catchment and the destination site's daily flows correspond to similar probabilities on their respective duration curves. The method is tested in several small catchments in South Africa. The method is designed primarily for application at ungauged sites in data‐poor regions where the use of more complex and information consuming techniques of data generation may not be justified. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

Raingauge network optimization in a tropical urban area by coupling cross-validation with the geostatistical technique     

Mohd Zaharifudin Muhamad Ali 《水文科学杂志》2018,63(3):474-491

An adequate and reliable raingauge network is essential for observing rainfall data in hydrology and water resource applications. A raingauge network developed for a catchment area is commonly extended periodically to increase data accuracy. Due to financial constraints, the network is reviewed for the optimal number of stations. A new optimization approach is developed in this study by coupling a cross-validation technique with a geostatistical method for raingauge network optimization to prioritize raingauge stations. The spatial interpolation error of the spatial rainfall distribution, measured as the root mean square error (E_rms) optimization criterion is applied to a raingauge network in a tropical urban area. The results indicate that this method can successfully optimize the number of rainfall stations in an existing raingauge network, as the stations are prioritized based on their importance in the network.  相似文献   

Towards flash-flood prediction in the dry Dead Sea region utilizing radar rainfall information     

Efrat Morin  Yael JacobyShilo Navon  Erez Bet-Halachmi 《Advances in water resources》2009

Flash-flood warning models can save lives and protect various kinds of infrastructure. In dry climate regions, rainfall is highly variable and can be of high-intensity. Since rain gauge networks in such areas are sparse, rainfall information derived from weather radar systems can provide useful input for flash-flood models. This paper presents a flash-flood warning model which utilizes radar rainfall data and applies it to two catchments that drain into the dry Dead Sea region. Radar-based quantitative precipitation estimates (QPEs) were derived using a rain gauge adjustment approach, either on a daily basis (allowing the adjustment factor to change over time, assuming available real-time gauge data) or using a constant factor value (derived from rain gauge data) over the entire period of the analysis. The QPEs served as input for a continuous hydrological model that represents the main hydrological processes in the region, namely infiltration, flow routing and transmission losses. The infiltration function is applied in a distributed mode while the routing and transmission loss functions are applied in a lumped mode. Model parameters were found by calibration based on the 5 years of data for one of the catchments. Validation was performed for a subsequent 5-year period for the same catchment and then for an entire 10-year record for the second catchment. The probability of detection and false alarm rates for the validation cases were reasonable. Probabilistic flash-flood prediction is presented applying Monte Carlo simulations with an uncertainty range for the QPEs and model parameters. With low probability thresholds, one can maintain more than 70% detection with no more than 30% false alarms. The study demonstrates that a flash-flood warning model is feasible for catchments in the area studied.  相似文献   

The impact of different gridding methods on catchment geomorphology and soil erosion over long timescales using a landscape evolution model     

G. R. Hancock 《地球表面变化过程与地形》2006,31(8):1035-1050

  相似文献   

Evaluation of satellite-based rainfall products for hydrological modelling in Morocco     

Yves Tramblay  Vera Thiemig  Alain Dezetter  Lahoucine Hanich 《水文科学杂志》2013,58(14):2509-2519

ABSTRACT

Several satellite-based precipitation estimates are becoming available at a global scale, providing new possibilities for water resources modelling, particularly in data-sparse regions and developing countries. This work provides a first validation of five different satellite-based precipitation products (TRMM-3B42 v6 and v7, RFE 2.0, PERSIANN-CDR, CMORPH1.0 version 0.x) in the 1785 km² Makhazine catchment (Morocco). Precipitation products are first compared against ground observations. Ten raingauges and four different interpolation methods (inverse distance, nearest neighbour, ordinary kriging and residual kriging with altitude) were used to compute a set of interpolated precipitation reference fields. Second, a parsimonious conceptual hydrological model is considered, with a simulation approach based on the random generation of model parameters drawn from existing parameter set libraries, to compare the different precipitation inputs. The results indicate that (1) all four interpolation methods, except the nearest neighbour approach, give similar and valid precipitation estimates at the catchment scale; (2) among the different satellite-based precipitation estimates verified, the TRMM-3B42 v7 product is the closest to observed precipitation, and (3) despite poor performance at the daily time step when used in the hydrological model, TRMM-3B42 v7 estimates are found adequate to reproduce monthly dynamics of discharge in the catchment. The results provide valuable perspectives for water resources modelling of data-scarce catchments with satellite-based rainfall data in this region.