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.

Comparative analysis of recent hydrological models and an attempt to generate new combined products for monitoring terrestrial water storage change

Hydrological models are crucial for characterizing large-scale water quantity variations and correcting GNSS reference station vertical displacements.We evaluated the robustness of multiple models,such as the Global Land Data Assimilation System (GLDAS),the Famine Early Warning System Network Land Data Assimilation System (FLDAS),the National Centers for Environmental Prediction (NCEP),and the WaterGAP Global Hydrology Model (WGHM).Inter-model and outer comparisons with Global Positioning System (GPS) coordinate time series,satellite gravity field Mascon solutions,and Global Precipitation Climatology Centre (GPCC) guide our assessment.Results confirm WGHM's 26% greater effectiveness in correcting nonlinear variations in GPS height time series compared to NCEP.In the Amazon River Basin,a 5-month lag between FLDAS,GLDAS,and satellite gravity results is observed.In eastern Asia and Australia,NCEP's Terrestrial Water Storage Changes (TWSC)-derived surface displacements correlate differently with precipitation compared to other models.Three combined hydrological models (H-VCE,H-EWM,and H-CVM) utilizing Variance Component Estimation (VCE),Entropy Weight Method (EWM),and Coefficient of Variation Method (CVM) are formulated.Correcting nonlinear variations with combined models enhances global GPS height scatter by 15%-17%.Correlation with precipitation increases by 25%-30%,and with satellite gravity,rises from 0.2 to 0.8 at maximum.The combined model eliminates time lag in the Amazon Basin TWSC analysis,exhibiting a four times higher signal-to-noise ratio than single models.H-VCE demonstrates the highest accuracy.In summary,the combined hydrological model minimizes discrepancies among individual models,significantly improving accuracy for monitoring large-scale TWSC.

Hydrological daily rainfall-runoff simulation with BTOPMC model and comparison with Xin'anjiang model

A grid-based distributed hydrological model, the Block-wise use of TOPMODEL （BTOPMC）, which was developed from the original TOPMODEL, was used for hydrological daily rainfall-runoff simulation. In the BTOPMC model, the runoff is explicitly calculated on a cell-by-cell basis, and the Muskingum-Cunge flow concentration method is used. In order to test the model＇s applicability, the BTOPMC model and the Xin＇anjiang model were applied to the simulation of a humid watershed and a semi-humid to semi-arid watershed in China. The model parameters were optimized with the Shuffle Complex Evolution （SCE-UA） method. Results show that both models can effectively simulate the daily hydrograph in humid watersheds, but that the BTOPMC model performs poorly in semi-humid to semi-arid watersheds. The excess-infiltration mechanism should be incorporated into the BTOPMC model to broaden the model＇s applicability.

Application of Hydrological Model to Simulate Rainfall-Runoff into An Khe Reservoir in the Ba River Basin, Vietnam

The BR （Ba River） basin is one of 9 main river basins in Vietnam. In the past 20 years, natural hazards such as flood andinundation have been complex and increased dramatically in both frequency and intensity in the BR basin. Recently, there have beenapproximately 198 large and small operating reservoirs which lead to increase natural hazards in the river basin. An Khe reservoir,one of big reservoirs in the upstream of the Ba river, impacts significantly on flooding in the downstream. This paper useshydrological model to simulate the flows as a basic for the safety operation of An Khe reservoir in order to prevent the downstreamfloods. The results indicate the Nash-Sutcliffe coefficient higher than 0.8 is stable and reliable parameter.

Rainfall-Runoff Modeling and Hydrological Responses to the Projected Climate Change for Upper Baro Basin, Ethiopia

This paper presents the results of Rainfall-Runoff modeling and simulation of hydrological responses under changing climate using HEC-HMS model. The basin spatial data was processed by HEC-GeoHMS and imported to HEC-HMS. The calibration and validation of the HEC-HMS model was done using the observed hydrometeorological data (1989-2018) and HEC-GeoHMS output data. The goodness-of-fit of the model was measured using three performance indices: Nash and Sutcliffe coefficient (NSE) = 0.8, Coefficient of Determination (R2) = 0.8, and Percent Difference (D) = 0.03, with values showing very good performance of the model. Finally, the optimized HEC-HMS model has been applied to simulate the hydrological responses of Upper Baro Basin to the projected climate change for mid-term (2040s) and long-term (2090s) A1B emission scenarios. The simulation results have shown a mean annual percent decrease of 3.6 and an increase of 8.1 for Baro River flow in the 2040s and 2090s scenarios, respectively, compared to the baseline period (2000s). A pronounced flow variation is rather observed on a seasonal basis, reaching a reduction of 50% in spring and an increase of 50% in autumn for both mid-term and long-term scenarios with respect to the base period. Generally, the rainfall-runoff model is developed to solve, in a complementary way, the two main problems in water resources management: the lack of gauged sites and future hydrological response to climate change data of the basin and the region in general. The study results imply that seasonal and time variation in the hydrologic cycle would most likely cause hydrologic extremes. And hence, the developed model and output data are of paramount importance for adaptive strategies and sustainable water resources development in the basin.

Light-absorbing Particles in Snow and Ice: Measurement and Modeling of Climatic and Hydrological impact

Light absorbing particles（LAP, e.g., black carbon, brown carbon, and dust） influence water and energy budgets of the atmosphere and snowpack in multiple ways. In addition to their effects associated with atmospheric heating by absorption of solar radiation and interactions with clouds, LAP in snow on land and ice can reduce the surface reflectance（a.k.a., surface darkening）, which is likely to accelerate the snow aging process and further reduces snow albedo and increases the speed of snowpack melt. LAP in snow and ice（LAPSI） has been identified as one of major forcings affecting climate change, e.g.in the fourth and fifth assessment reports of IPCC. However, the uncertainty level in quantifying this effect remains very high. In this review paper, we document various technical methods of measuring LAPSI and review the progress made in measuring the LAPSI in Arctic, Tibetan Plateau and other mid-latitude regions. We also report the progress in modeling the mass concentrations, albedo reduction, radiative forcing, and climatic and hydrological impact of LAPSI at global and regional scales. Finally we identify some research needs for reducing the uncertainties in the impact of LAPSI on global and regional climate and the hydrological cycle.

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.

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.

Application of Radar-Measured Rain Data in Hydrological Processes Modeling during the Intensified Observation Period of HUBEX

On the basis of Digital Elevation Model data, the raster flow vectors, watershed delineation, and spatial topological relationship are generated by the Martz and Garbrecht method for the upper area of Huangnizhuang station in the Shihe Catchment with 805 km2 of area, an intensified observation field for the HUBEX/GAME Project. Then, the Xin’anjiang Model is applied for runoff production in each grid element where rain data measured by radar at Fuyang station is utilized as the input of the hydrological model. The elements are connected by flow vectors to the outlet of the drainage catchment where runoff is routed by the Muskingum method from each grid element to the outlet according to the length between each grid and the outlet. The Nash-Sutcliffe model efficiency coefficient is 92.41% from 31 May to 3 August 1998, and 85.64%, 86.62%, 92.57%, and 83.91%, respectively for the 1st, 2nd, 3rd, and 4th flood events during the whole computational period. As compared with the case where rain-gauge data are used in simulating the hourly hydrograph at Huangnizhuang station in the Shihe Catchment, the index of model efficiency improvement is positive, ranging from 27.56% to 69.39%. This justifies the claim that radar-measured data are superior to rain-gauge data as inputs to hydrological modeling. As a result, the grid-based hydrological model provides a good platform for runoff computation when radar-measured rain data with highly spatiotemporal resolution are taken as the input of the hydrological model.

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.

Distributed hydrological models for addressing effects of spatial variability of roughness on overland flow

In this study, we investigated the origin of the overland flow roughness problem and divided the current overland flow roughness research into three types, as follows: the first type of research takes into account the effects of roughness on the volume and velocity of surface runoff, flood peaks, and the scouring capability of flows, but has not addressed the spatial variability of roughness in detail; the second type of research considers that surface roughness varies spatially with different land usage types, land-cover conditions, and different tillage forms, but lacks a quantitative study of the spatial variability; and the third type of research simply deals with the spatial variability of roughness in each grid cell or land type. We present three shortcomings of the current overland flow roughness research, including(1) the neglect of roughness in distributed hydrological models when simulating the overland flow direction and distribution,(2) the lack of consideration of spatial variability of roughness in hydrological models, and(3) the failure to distinguish the roughness formulas in different overland flow regimes. To solve these problems,distributed hydrological model research should focus on four aspects in regard to overland flow: velocity field observations, flow regime mechanisms, a basic roughness theory, and scale problems.

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.

Sensitivity Analysis and Calibration of Hydrological Modeling of the Watershed Northeast Brazil

Mathematical models of the quantity and quality of water in hydrographic basins enable simulation of a wide variety of processes, including the production of water and sediments, and the dynamics of point and nonpoint sources of pollution. These models have become increasingly complex, requiring large amounts of input data, which can increase the uncertainty of the results of simulations. For this reason, it is essential to perform calibration and validation procedures. The objective of this work was to conduct sensitivity analysis and calibration of a distributed hydrological model (SWAT) applied to the flows of water in the watershed of the Poxim River. Satisfactory performance of the model was indicated by the values obtained for the Nash-Sutcliffe efficiency coefficient (0.77), the percent bias (5.05), the root mean square error (0.48), and the ratio of the RMSE to the standard deviation of the observations (RSR) (0.49). The set of parameters identified here could be used for the simulation and evaluation of other scenarios.

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.

Understanding and simulating of three-dimensional subsurface hydrological partitioning in an alpine mountainous area, China

Critical zone(CZ)plays a vital role in sustaining biodiversity and humanity.However,flux quantification within CZ,particularly in terms of subsurface hydrological partitioning,remains a significant challenge.This study focused on quantifying subsurface hydrological partitioning,specifically in an alpine mountainous area,and highlighted the important role of lateral flow during this process.Precipitation was usually classified as two parts into the soil:increased soil water content(SWC)and lateral flow out of the soil pit.It was found that 65%–88%precipitation contributed to lateral flow.The second common partitioning class showed an increase in SWC caused by both precipitation and lateral flow into the soil pit.In this case,lateral flow contributed to the SWC increase ranging from 43%to 74%,which was notably larger than the SWC increase caused by precipitation.On alpine meadows,lateral flow from the soil pit occurred when the shallow soil was wetter than the field capacity.This result highlighted the need for three-dimensional simulation between soil layers in Earth system models(ESMs).During evapotranspiration process,significant differences were observed in the classification of subsurface hydrological partitioning among different vegetation types.Due to tangled and aggregated fine roots in the surface soil on alpine meadows,the majority of subsurface responses involved lateral flow,which provided 98%–100%of evapotranspiration(ET).On grassland,there was a high probability(0.87),which ET was entirely provided by lateral flow.The main reason for underestimating transpiration through soil water dynamics in previous research was the neglect of lateral root water uptake.Furthermore,there was a probability of 0.12,which ET was entirely provided by SWC decrease on grassland.In this case,there was a high probability(0.98)that soil water responses only occurred at layer 2(10–20 cm),because grass roots mainly distributed in this soil layer,and grasses often used their deep roots for water uptake during ET.To improve the estimation of soil water dynamics and ET,we established a random forest(RF)model to simulate lateral flow and then corrected the community land model(CLM).RF model demonstrated good performance and led to significant improvements in CLM simulation.These findings enhance our understanding of subsurface hydrological partitioning and emphasize the importance of considering lateral flow in ESMs and hydrological research.

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.

Developing a Framework to Measure Watershed Sustainability by Using Hydrological/Water Quality Model

A framework is built, wherein hydrological/water quality model is used to measure watershed sustainability. For this framework, watershed sustainability has been defined and quantified by defining social, environmental and biodiversity indicators. By providing weightage to these indicators, a “River Basin Sustainability Index” is built. The watershed sustainability is then calculated based on the concepts of reliability, resilience and vulnerability. The framework is then applied to a case study, where, based on watershed management principles, four land use scenarios are created in GIS. The Soil and Water Assessment Tool (SWAT) is used as a hydrology/water quality model. Based on the results the land uses are ranked for sustainability and policy implications have been discussed. This results show that landuse (both type and location) impact watershed sustainability. The existing land use is weak in environmental sustainability. Also, riparian zones play a critical role in watershed sustainability, although beyond certain width their contribution is not significant.

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.

Isotope hydrogeochemical models for assessing the hydrological processes in a part of the largest continental flood basalts province of India

Continental Flood Basalts(CFB)occupy one fourth of the world’s land area.Hence,it is important to discern the hydrological processes in this complex hydrogeological setup for the sustainable water resources development.A model assisted isotope,geochemical,geospatial and geophysical study was conducted to understand the monsoonal characteristics,recharge processes,renewability and geochemical evolution in one of the largest continental flood basalt provinces of India.HYSPLIT modelling and stable isotopes were used to assess the monsoonal characteristics.Rayleigh distillation model were used to understand the climatic conditions at the time of groundwater recharge.Lumped parameter models(LPM)were employed to quantify the mean transit time(MTT)of groundwater.Statistical and geochemical models were adopted to understand the geochemical evolution along the groundwater flow path.A geophysical model was used to understand the geometry of the aquifer.The back trajectory analysis confirms the isotopic finding that precipitation in this region is caused by orographic uplifting of air masses originating from the Arabian Sea.Stable isotopic data of groundwater showed its meteoric origin and two recharge processes were discerned;(i)quick and direct recharge by precipitation through fractured and weathered basalt,(ii)low infiltration through the clayey black cotton soil and subjected to evaporation prior to the recharge.Tritium data showed that the groundwater is a renewable source and have shorter transit times(from present day to<30 years).The hydrogeochemical study indicated multiple sources/processes such as:the minerals dissolution,silicate weathering,ion exchange,anthropogenic influences etc.control the chemistry of the groundwater.Based on the geo-electrical resistivity survey,the potential zones(weathered and fractured)were delineated for the groundwater development.Thus,the study highlights the usefulness of model assisted isotopic hydrogeochemical techniques for understanding the recharge and geochemical processes in a basaltic aquifer system.

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.