Application of hydrological models in a snowmelt region of Aksu River Basin

This study simulated and predicted the runoff of the Aksu River Basin, a typical river basin supplied by snowmelt in an arid mountain region, with a limited data set and few hydrological and meteorological stations. Two hydrological models, the snowmelt-runoff model （SRM） and the Danish NedbФr-AfstrФmnings rainfall-runoff model （NAM）, were used to simulate daily discharge processes in the Aksu River Basin. This study used the snow-covered area from MODIS remote sensing data as the SRM input. With the help of ArcGIS software, this study successfully derived the digital drainage network and elevation zones of the basin from digital elevation data. The simulation results showed that the SRM based on MODIS data was more accurate than NAM. This demonstrates that the application of remote sensing data to hydrological snowmelt models is a feasible and effective approach to runoff simulation and prediction in arid unguaged basins where snowmelt is a major runoff factor.

A review of integrated surface-subsurface numerical hydrological models

Hydrological modeling,leveraging mathematical formulations to represent the hydrological cycle,is a pivotal tool in representing the spatiotemporal dynamics and distribution patterns inherent in hydrology.These models serve a dual purpose:they validate theoretical robustness and applicability via observational data and project future trends,thereby bridging the understanding and prediction of natural processes.In rapid advancements in computational methodologies and the continuous evolution of observational and experimental techniques,the development of numerical hydrological models based on physicallybased surface-subsurface process coupling have accelerated.Anchored in micro-scale conservation principles and physical equations,these models employ numerical techniques to integrate surface and subsurface hydrodynamics,thus replicating the macro-scale hydrological responses of watersheds.Numerical hydrological models have emerged as a leading and predominant trend in hydrological modeling due to their explicit representation of physical processes,heightened by their spatiotemporal resolution and reliance on interdisciplinary integration.This article focuses on the theoretical foundation of surface-subsurface numerical hydrological models.It includes a comparative and analytical discussion of leading numerical hydrological models,encompassing model architecture,numerical solution strategies,spatial representation,and coupling algorithms.Additionally,this paper contrasts these models with traditional hydrological models,thereby delineating the relative merits,drawbacks,and future directions of numerical hydrological modeling.

Analysing the Potential Impact of Climate Change on the Hydrological Regime of the Upper Benue River Basin (North Cameroon)

In this study, we analyse the climate variability in the Upper Benue basin and assess its potential impact on the hydrology regime under two different greenhouse gas emission scenarios. The hydrological regime of the basin is more vulnerable to climate variability, especially precipitation and temperature. Observed hydroclimatic data (1950-2015) was analysed using a statistical approach. The potential impact of future climate change on the hydrological regime is quantified using the GR2M model and two climate models: HadGEM2-ES and MIROC5 from CMIP5 under RCP 4.5 and RCP 8.5 greenhouse gas emission scenarios. The main result shows that precipitation varies significantly according to the geographical location and time in the Upper Benue basin. The trend analysis of climatic parameters shows a decrease in annual average precipitation across the study area at a rate of -0.568 mm/year which represents about 37 mm/year over the time 1950-2015 compared to the 1961-1990 reference period. An increase of 0.7°C in mean temperature and 14% of PET are also observed according to the same reference period. The two climate models predict a warming of the basin of about 2°C for both RCP 4.5 and 8.5 scenarios and an increase in precipitation between 1% and 10% between 2015 and 2100. Similarly, the average annual flow is projected to increase by about +2% to +10% in the future for both RCP 4.5 and 8.5 scenarios between 2015 and 2100. Therefore, it is primordial to develop adaptation and mitigation measures to manage efficiently the availability of water resources.

Integrated Hydrological Modeling of the Godavari River Basin in Maharashtra Using the SWAT Model: Streamflow Simulation and Analysis

Hydrological modeling plays a crucial role in efficiently managing water resources and understanding the hydrologic behavior of watersheds. This study aims to simulate daily streamflow in the Godavari River Basin in Maharashtra using the Soil and Water Assessment Tool (SWAT). SWAT is a process-based hydrological model used to predict water balance components, sediment levels, and nutrient contamination. In this research, we used integrated remote sensing and GIS data, including Digital Elevation Models (DEM), land use and land cover (LULC) maps, soil maps, and observed precipitation and temperature data, as input for developing the SWAT model to assess surface runoff in this large river basin. The Godavari River Basin under study was divided into 25 sub-basins, comprising 151 hydrological response units categorized by unique land cover, soil, and slope characteristics using the SWAT model. The model was calibrated and validated against observed runoff data for two time periods: 2003-2006 and 2007-2010 respectively. Model performance was assessed using the Nash-Sutcliffe efficiency (NSE) and the coefficient of determination (R2). The results show the effectiveness of the SWAT2012 model, with R2 value of 0.84 during calibration and 0.86 during validation. NSE values also ranged from 0.84 during calibration to 0.85 during validation. These findings enhance our understanding of surface runoff dynamics in the Godavari River Basin under study and highlight the suit-ability of the SWAT model for this region.

Flood Forecasting and Warning System: A Survey of Models and Their Applications in West Africa

Flood events occurrences and frequencies in the world are of immense worry for the stability of the economy and life safety. Africa continent is the third continent the most negatively affected by the flood events after Asia and Europe. Eastern Africa is the most hit in Africa. However, Africa continent is at the early stage in term of flood forecasting models development and implementation. Very few hydrological models for flood forecasting are available and implemented in Africa for the flood mitigation. And for the majority of the cases, they need to be improved because of the time evolution. Flash flood in Bamako (Mali) has been putting both human life and the economy in jeopardy. Studying this phenomenon, as to propose applicable solutions for its alleviation in Bamako is a great concern. Therefore, it is of upmost importance to know the existing scientific works related to this situation in Mali and elsewhere. The main aim was to point out the various solutions implemented by various local and international institutions, in order to fight against the flood events. Two types of methods are used for the flood events adaptation: the structural and non-structural methods. The structural methods are essentially based on the implementation of the structures like the dams, dykes, levees, etc. The problem of these methods is that they may reduce the volume of water that will inundate the area but are not efficient for the prediction of the coming floods and cannot alert the population with any lead time in advance. The non-structural methods are the one allowing to perform the prediction with acceptable lead time. They used the hydrological rainfall-runoff models and are the widely methods used for the flood adaptation. This review is more accentuated on the various types non-structural methods and their application in African countries in general and West African countries in particular with their strengths and weaknesses. Hydrologiska Byråns Vattenbalansavdelning (HBV), Hydrologic Engineer Center Hydrologic Model System (HEC-HMS) and Soil and Water Assessment Tool (SWAT) are the hydrological models that are the most widely used in West Africa for the purpose of flood forecasting. The easily way of calibration and the weak number of input data make these models appropriate for the West Africa region where the data are scarce and often with bad quality. These models when implemented and applied, can predict the coming floods, allow the population to adapt and mitigate the flood events and reduce considerably the impacts of floods especially in terms of loss of life.

Evapotranspiration estimation methods in hydrological models

Actual evapotranspiration is a key process of hydrological cycle and a sole term that links land surface water balance and land surface energy balance. Evapotranspiration plays a key role in simulating hydrological effect of climate change, and a review of evapotranspiration estimation methods in hydrological models is of vital importance. This paper firstly summarizes the evapotranspiration estimation methods applied in hydrological models and then classifies them into the integrated converting methods and the classification gathering methods by their mechanism. Integrated converting methods are usually used in hydrological models and two differences exist among them： one is in the potential evaporation estimation methods, while the other in the function for defining relationship between potential evapora tion and actual evapotranspiration. Due to the higher information requirements of the Pen- man-Monteith method and the existing data uncertainty, simplified empirical methods for calculating potential and actual evapotranspiration are widely used in hydrological models. Different evapotranspiration calculation methods are used depending on the complexity of the hydrological model, and importance and difficulty in the selection of the most suitable evapotranspiration methods is discussed. Finally, this paper points out the prospective de velopment trends of the evapotranspiration estimating methods in hydrological modeling.

A review on evapotranspiration data assimilation based on hydrological models

Accurate estimation of evapotranspiration(ET),especially at the regional scale,is an extensively investigated topic in the field of water science. The ability to obtain a continuous time series of highly precise ET values is necessary for improving our knowledge of fundamental hydrological processes and for addressing various problems regarding the use of water. This objective can be achieved by means of ET data assimilation based on hydrological modeling. In this paper,a comprehensive review of ET data assimilation based on hydrological modeling is provided. The difficulties and bottlenecks of using ET,being a non-state variable,to construct data assimilation relationships are elaborated upon,with a discussion and analysis of the feasibility of assimilating ET into various hydrological models. Based on this,a new easy-to-operate ET assimilation scheme that includes a water circulation physical mechanism is proposed. The scheme was developed with an improved data assimilation system that uses a distributed time-variant gain model(DTVGM),and the ET-soil humidity nonlinear time response relationship of this model. Moreover,the ET mechanism in the DTVGM was improved to perfect the ET data assimilation system. The new scheme may provide the best spatial and temporal characteristics for hydrological states,and may be referenced for accurate estimation of regional evapotranspiration.

Separating impacts of vegetation change and climate variability on streamflow using hydrological models together with vegetation data

Vegetation information is seldom considered in lumped conceptual rainfall-runoff models.This paper uses two modified rainfall-runoff models,the Xinanjiang-ET and SIMHYD-ET models in which vegetation leaf area index is incorporated,to investigate impacts of vegetation change and climate variability on streamflow in a Southern Australian catchment,the Crawford River experimental catchment,where Tasmanian blue gum plantations were introduced gradually from 1998 till 2005.The Xinanjiang-ET and SIMHYD-ET models incorporate remotely-sensed leaf area index(LAI) data obtained from the Advanced Very High Resolution Radiometer(AVHRR) on board NOAA polar orbiting satellites.Compared to the original versions,the Xinanjiang-ET and SIMHYD-ET models show marginal improvements in runoff simulations in the pre-plantation period(1882-1997).The calibrated Xinanjaing-ET and SIMHYD-ET models are then used to simulate plantation impact on streamflow in the post-plantation period.The total change in streamflow between the pre-plantation and post-plantation periods is 32.4 mm/a.The modelling results from the two models show that plantation reduces streamflow by 20.5 mm/a,and climate variability reduces streamflow by 11.9 mm/a.These results suggest that increase in plantations can reduce streamflow substantially,even more than climate variability.

Analysis of Hydrological Simulation Models Using the Parameter Combinatorial Diagram

The aim of this paper is to present graphically the behaviour of a simulation model to the varying parameters and to establish the suitability of this representation as a valid tool for the analysis of the same parameters. In this paper, we define parameter combinatorial diagram as the joint graphical representation of all box plots related to the adjustment between real and simulated data, by setting and/or changing the parameters of the simulation model. To do this, we start with a box plot representing the values of an objective adjustment function, achieving these results when varying all the parameters of the simulation model, Then we draw the box plot when setting all the parameters of the model, for example, using the median or average. Later, we get all the box plots when carrying out simulations combining fixed or variable values of the model parameters. Finally, all box plots obtained are represented neatly in a single graph. It is intended that the new parameter combinatorial diagram is used to examine and analyze simulation models useful in practice. This paper presents combinatorial diagrams of different examples of application as in the case of hydrologic models of one, two, three, and five parameters.

Modelling Hydrological Consequences of Climate Change—Progress and Challenges

The simulation of hydrological consequences of climate change has received increasing attention from the hydrology and land-surface modelling communities. There have been many studies of climate-change effects on hydrology and water resources which usually consist of three steps： （1） use of general circulation models （GCMs） to provide future global climate scenarios under the effect of increasing greenhouse gases, （2） use of downscaling techniques （both nested regional climate models, RCMs, and statistical methods） for ＂downscaling＂ the GCM output to the scales compatible with hydrological models, and （3） use of hydrologic models to simulate the effects of climate change on hydrological regimes at various scales. Great progress has been achieved in all three steps during the past few years, however, large uncertainties still exist in every stage of such study. This paper first reviews the present achievements in this field and then discusses the challenges for future studies of the hydrological impacts of climate change.

Multi-Model Approach for Assessing the Influence of Calibration Criteria on the Water Balance in Ouémé Basin

Hydrological models are very useful tools for evaluating water resources, and the hydroclimatic hazards associated with the water cycle. However, their calibration and validation require the use of performance criteria which choice is not straightforward. This paper aims to evaluate the influence of the performance criteria on water balance components and water extremes using two global rainfall-runoff models (HBV and GR4J) over the Ouémé watershed at the Bonou and Savè outlets. Three (3) Efficacy criteria (Nash, coefficient of determination, and KGE) were considered for calibration and validation. The results show that the Nash criterion provides a good assessment of the simulation of the different parts of the hydrograph. KGE is better for simulating peak flows and water balance elements than other efficiency criteria. This study could serve as a basis for the choice of performance criteria in hydrological modelling.

GIS Application in Urban Flood Risk Analysis: Midar as a Case Study

The significance of this study lies in its exploration of the advanced applications of Geographic Information Systems (GIS) in assessing urban flood risks, with a specific focus on Midar, Morocco. This research is pivotal as it showcases that GIS technology is not just a tool for mapping, but a critical component in urban planning and emergency management strategies. By meticulously identifying and mapping flood-prone areas in Midar, the study provides invaluable insights into the potential vulnerabilities of urban landscapes to flooding. Moreover, this research demonstrates the practical utility of GIS in mitigating material losses, a significant concern in flood-prone urban areas. The proactive approach proposed in this study, centered around the use of GIS, aims to safeguard Midar’s population and infrastructure from the devastating impacts of floods. This approach serves as a model for other urban areas facing similar challenges, highlighting the indispensable role of GIS in disaster preparedness and response. Overall, the study underscores the transformative potential of GIS in enhancing urban resilience, making it a crucial tool in the fight against natural disasters like floods.

Simulations of a Hydrological Model as Coupled to a Regional Climate Model

Considering a detailed hydrologic model in the land surface scheme helps to improve the simulation of regional hydro-climatology. A hydrologic model, which includes spatial heterogeneities in precipitation and infiltration, is constructed and incorporated into the land surface scheme BATS. Via the coupled-model (i.e., a regional climate model) simulations, the following major conclusions are obtained: the simulation of surface hydrology is sensitive to the inclusion of heterogeneities in precipitation and infiltration; the runoff ratio is increased after considering the infiltration heterogeneity, a result which is more consistent with the observations of surface moisture balance over humid areas; the introduction of the parameterization of infiltration heterogeneity can have a greater influence on the regional hydro-climatology than the precipitation heterogeneity; and the consideration of the impermeable fraction for the region reveals some features that are closer to the trend of aridification over northern China.

Hydrological characteristics and changes in the Nu-Salween River basin revealed with model-based reconstructed data

The Nu-Salween River(NSR),the longest free-flow river in Southeast Asia,plays an irreplaceable role in social development and ecological protection.The lower NSR region is particularly valuable as it is inhabited by approximately 6.7 million people.The basin has limited hydraulic conservancy infrastructure and insufficient ability to cope with climate change risks.Studying the hydrological characteristics and changes in the basin provides the scientific basis for rational protection and development of the basin.However,owing to the limitation of observation data,previous studies have focused on the local area and neglected the study of the lower reaches,which is not enough to reflect the spatial characteristics of the entire basin.In this study,the ECMWF 5th generation reanalysis data(ERA5)and Multi-Source Weighted-Ensemble Precipitation(MSWEP)were applied to develop a geomorphology-based hydrological model(GBHM)for reconstructing hydrological datasets(i.e.GBHM-ERA5 and GBHM-MSWEP).The reconstructed datasets covering the complete basin were verified against the gauge observation and compared with other commonly used streamflow products,including Global Flood Awareness System v2.1,GloFAS-Reanalysis dataset v3.0,and linear optimal runoff aggregate(LORA).The comparison results revealed that GBHM-ERA5 is significantly better than the other four datasets and provides a good reproduction of the hydrological characteristics and trends of the NSR.Detailed analysis of GBHM-ERA5 revealed that:(1)A multi-year mean surface runoff represented 39%of precipitation over the basin during 1980–2018,which had low surface runoff in the upstream,while areas around the Three Parallel Rivers Area and the estuary had abundant surface runoff.(2)The surface runoff and discharge coefficient of variations in spring were larger than those in other seasons,and the inter-annual variation in the downstream was smaller than that in the upstream and midstream regions.(3)More than 70%of the basin areas showed a decreasing trend in the surface runoff,except for parts of Nagqu,south of Shan State in Myanmar,and Thailand,where surface runoff has an increasing trend.(4)The downstream discharge has dropped significantly at a rate of approximately 680 million cubic metresper year,and the decline rate is greater than that of upstream and midstream,especially in summer.This study provides a data basis for subsequent studies in the NSR basin and further elucidates the impact of climate change on the basin,which is beneficial to river planning and promotes international cooperation on the water-and eco-security of the basin.

Hydrological behaviour and water balance analysis for Xitiaoxi catchment of Taihu Basin

With the rapid social and economic development of the Taihu region, Taihu Lake now faces an increasingly severe eutrophication problem. Pollution from surrounding catchments contributes greatly to the eutrophication of water bodies in the region. Investigation of surface flow and associated mass transport for the Xitiaoxi catchment is of a significant degree of importance as the Xitiaoxi catchment is one of the major catchments within the Taihu region. A SWAT-based distributed hydrological model was established for the Xitiaoxi catchment. The model was calibrated and verified using hydrometeorological data from 1988 to 2001. The results indicate that the modeled daily and annual stream flow match the observed data both in the calibration period and the verification period, with a linear regression coefficient R2 and a coefficient e for modeled daily stream flow greater than 0.8 at Hengtangcun and Fanjiacun gauge stations. The results show that the runoff process in the Xitiaoxi catchment is affected both by rainfall and human activities （e.g., reservoirs and polder areas）. Moreover, the human activities weaken flood peaks more noticeably during rainstorms. The Water balance analysis reveals the percentages of precipitation made up by surface flow, evapotranspiration, groundwater recharge＇ and the change of soil storage, all of which are considered useful to the further understanding of the hydrological processes in the Xitiaoxi catchment. This study provides a good base for further studies in mass transport modeling and comparison of modeling results from similar hydrological models.

Effect of calibration data series length on performance and optimal parameters of hydrological model

In order to assess the effects of calibration data series length on the performance and optimal parameter values of a hydrological model in ungauged or data-limited catchments （data are non-continuous and fragmental in some catchments）, we used non-continuous calibration periods for more independent streamflow data for SIMHYD （simple hydrology） model calibration. Nash-Sutcliffe efficiency and percentage water balance error were used as performance measures. The particle swarm optimization （PSO） method was used to calibrate the rainfall-runoff models. Different lengths of data series ranging from one year to ten years, randomly sampled, were used to study the impact of calibration data series length. Fifty-five relatively unimpaired catchments located all over Australia with daily precipitation, potential evapotranspiration, and streamflow data were tested to obtain more general conclusions. The results show that longer calibration data series do not necessarily result in better model performance. In general, eight years of data are sufficient to obtain steady estimates of model performance and parameters for the SIMHYD model. It is also shown that most humid catchments require fewer calibration data to obtain a good performance and stable parameter values. The model performs better in humid and semi-humid catchments than in arid catchments. Our results may have useful and interesting implications for the efficiency of using limited observation data for hydrological model calibration in different climates.

Spatial distribution characteristics and mechanism of nonhydrological time-variable gravity in China's Mainland

The purpose of this study is to explore nonhydrological mass transfer in China's Mainland.For this purpose,gravity recovery and climate experiment(GRACE)data were obtained to study the spatial distribution of time variant gravity signals in China's Mainland.Then,from auxiliary hydrological data processed according to the current hydrological model,a new more comprehensive hydrological model of China's Mainland was constructed.Finally,the time variant signals of this new hydrological model were removed from the time variant gravity field computed from GRACE data,thus obtaining a description of the nonhydrological mass transfer of China's Mainland.The physical sources and mechanisms of the resulting mass transfer are then discussed.The improved,more realistic,hydrological model used here was created by selecting the hydrological components with the best correlations in existing hydrological models,by use of correlation calculation,analysis,and comparison.This improved model includes water in soils and deeper strata,in the vegetation canopy,in lakes,snow,and glaciers,and in other water components(mainly reservoir storage,swamps,and rivers).The spatial distribution of the transfer signals due to nonhydrological mass in China's Mainland was obtained by subtracting the combined hydrological model from the GRACE time-variable gravity field.The results show that the nonhydrological signals in China's Mainland collected in GRACE data were mainly positive signals,and were distributed in the Bohai Rim and the northern and eastern parts of the Tibetan Plateau.The above nonhydrological mass transfer signals have been studied further and are discussed.The results show that the nonhydrological mass migration signals in the Bohai Rim region originate primarily from sea level change and marine sediment accumulation.The mass accumulation from Indian plate collision in the Tibetan Plateau appears to be the main reason for the increase in the residual gravity field in that region.

Hydrological Response to Environment Change in Himalayan Watersheds: Assessment from Integrated Modeling Approach

Land use changes such as deforestation,increase in cropping or grazing areas and built-up land, likely modify the water balance and land surface behavior in the Himalayan watersheds.An integrated approach of hydrological and hydraulic modeling was adopted for comparative analysis of hydrological pattern in three Himalayan watersheds i.e.Khanpur,Rawal and Simly situated in the Northern territory of Pakistan.The rainfall-runoff model SWAT- Soil and water assessment tool and Hydro CAD were calibrated for the selected watersheds.The correlation analysis of the precipitation data of two climate stations i.e.Murree and Islamabad, with the discharge data of three rivers was utilized to select best suitable input precipitation data for Hydro CAD rainfall-runoff modeling.The peak flood hydrograph were generated using Hydro CAD runoff to optimize the basin parameters like CN, runoff volume, peak flows of the three watersheds.The hydrological response of the Rawal watershed was studied as a case study to different scenarios of land use change using SWAT model.The scenario of high deforestation indicated a decline of about 6.3% in the groundwater recharge tostream while increase of 7.1% in the surface runoff has been observed under the scenario of growth in urbanization in the recent decades.The integrated modeling approach proved helpful in investigating the hydrological behavior under changing environment at watershed level in the Himalayan region.

Optimization of the Hydrological Model Using Multi-objective Particle Swarm Optimization Algorithm

An application of multi-objective particle swarm optimization （MOPSO） algorithm for optimization of the hydrological model （HYMOD） is presented in this paper. MOPSO algorithm is used to find non-dominated solutions with two objectives： high flow Nash-Sutcliffe efficiency and low flow Nash- Sutcliffe efficiency. The two sets＇ coverage rate and Pareto front spacing metric are two criterions to analyze the performance of the algorithms. MOPSO algorithm surpasses multi-objective shuffled complex evolution metcopolis （MOSCEM_UA） algorithr~, in terms of the two sets＇ coverage rate. But when we come to Pareto front spacing rate, the non-dominated solutions of MOSCEM_ UA algorithm are better-distributed than that of MOPSO algorithm when the iteration is set to 40 000. In addition, there are obvious conflicts between the two objectives. But a compromise solution can be acquired by adopting the MOPSO algorithm.

Projection of future streamflow of the Hunza River Basin,Karakoram Range(Pakistan)using HBV hydrological model

Hydrologiska Byrans Vattenbalansavdeling(HBV) Light model was used to evaluate the performance of the model in response to climate change in the snowy and glaciated catchment area of Hunza River Basin. The study aimed to understand the temporal variation of streamflow of Hunza River and its contribution to Indus River System(IRS). HBV model performed fairly well both during calibration(R2=0.87, Reff=0.85, PBIAS=-0.36) and validation(R2=0.86, Reff=0.83, PBIAS=-13.58) periods on daily time scale in the Hunza River Basin. Model performed better on monthly time scale with slightly underestimated low flows period during bothcalibration(R2=0.94, Reff=0.88, PBIAS=0.47) and validation(R2=0.92, Reff=0.85, PBIAS=15.83) periods. Simulated streamflow analysis from 1995-2010 unveiled that the average percentage contribution of snow, rain and glacier melt to the streamflow of Hunza River is about 16.5%, 19.4% and 64% respectively. In addition, the HBV-Light model performance was also evaluated for prediction of future streamflow in the Hunza River using future projected data of three General Circulation Model(GCMs) i.e. BCC-CSM1.1, CanESM2, and MIROCESM under RCP2.6, 4.5 and 8.5 and predictions were made over three time periods, 2010-2039, 2040-2069 and 2070-2099, using 1980-2010 as the control period. Overall projected climate results reveal that temperature and precipitation are the most sensitiveparameters to the streamflow of Hunza River. MIROC-ESM predicted the highest increase in the future streamflow of the Hunza River due to increase in temperature and precipitation under RCP4.5 and 8.5 scenarios from 2010-2099 while predicted slight increase in the streamflow under RCP2.6 during the start and end of the 21 th century. However, BCCCSM1.1 predicted decrease in the streamflow under RCP8.5 due to decrease in temperature and precipitation from 2010-2099. However, Can ESM2 predicted 22%-88% increase in the streamflow under RCP4.5 from 2010-2099. The results of this study could be useful for decision making and effective future strategic plans for water management and their sustainability in the region.