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.

SimET: An open-source tool for estimating crop evapotranspiration and soil water balance for plants with multiple growth cycles

Accurate estimation of crop evapotranspiration(ETc) and soil water balance, which is vital for optimizing water management strategy in crop production, can be performed by simulation. But existing software has many deficiencies, including complex operation, limited scalability, lack of batch processing, and a single ETc model. Here we present simET, an open-source software package written in the R programming language. Many concepts involved in crop ETc simulation are condensed into functions in the package. It includes three widely used crop ETc models built on these functions: the single-crop coefficient,double-crop coefficient, and Shuttleworth–Wallace models, along with tools for preparing model data and comparing estimates. SimET supports ETc simulation in crops with repeated growth cycles such as alfalfa, a perennial forage crop that is cut multiple times annually.

Estimation of water balance in the source region of the Yellow River based on GRACE satellite data

Water storage has important significance for understanding water cycles of global and local domains and for monitoring climate and environmental changes. As a key variable in hydrology, water storage change represents the sum of precipitation, evaporation, surface runoff, soil water and groundwater exchanges. Water storage change data during the period of 2003-2008 for the source region of the Yellow River were collected from Gravity Recovery and Climate Experiment （GRACE） satellite data. The monthly actual evaporation was estimated according to the water balance equation. The simulated actual evaporation was significantly consistent and correlative with not only the observed pan （20 cm） data, but also the simulated results of the version 2 of Simple Biosphere model. The average annual evaporation of the Tangnaihai Basin was 506.4 mm, where evaporation in spring, summer, autumn and winter was 130.9 mm, 275.2 mm, 74.3 mm and 26.1 mm, and accounted for 25.8%, 54.3%, 14.7% and 5.2% of the average annual evaporation, respectively, The precipitation increased slightly and the actual evaporation showed an obvious decrease. The water storage change of the source region of the Yellow River displayed an increase of 0.51 mm per month from 2003 to 2008, which indicated that the storage capacity has significantly increased, probably caused by the degradation of permafrost and the increase of the thickness of active layers. The decline of actual evaporation and the increase of water storage capacity resulted in the increase of river runoff.

Understanding the impact of mountain landscapes on water balance in the upper Heihe River watershed in northwestern China

Estimating the impact of mountain landscape on hydrology or water balance is essential for the sus- tainable development strategies of water resources. Specifically, understanding how the change of each landscape influences hydrological components will greatly improve the predictability of hydrological responses to mountain landscape changes and thus can help the government make sounder decisions. In the paper, we used the VIC （Variable Infiltration Capacity） model to conduct hydrological modeling in the upper Heihe River watershed, along with a frozen-soil module and a glacier melting module to improve the simulation. The improved model performed satisfactorily. We concluded that there are differences in the runoff generation of mountain landscape both in space and time. About 50% of the total runoff at the catchment outlet were generated in mid-mountain zone （2,900-4,000 m asl）, and water was mainly consumed in low mountain region （1,700-2,900 m asl） because of the higher requirements of trees and grasses. The runoff coefficient was 0.37 in the upper Heihe River watershed. Barren landscape produced the largest runoff yields （52.46% of the total runoff） in the upper Heihe River watershed, fol- lowed by grassland （34.15%）, shrub （9.02%）, glacier （3.57%）, and forest （0.49%）. In order to simulate the impact of landscape change on hydrological components, three landscape change scenarios were designed in the study. Scenario 1, 2 and 3 were to convert all shady slope landscapes at 2,000-3,300 m, 2,000-3,700 m, and 2,000-4,000 m asl respectively to forest lands, with forest coverage rate increased to 12.4%, 28.5% and 42.0%, respectively. The runoff at the catchment outlet correspondingly declined by 3.5%, 13.1% and 24.2% under the three scenarios. The forest landscape is very important in water conservation as it reduced the flood peak and increased the base flow. The mountains as ＂water towers＂ play important roles in water resources generation and the impact of mountain landscapes on hydrology is significant.

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.

Hydrological Impact Assessment of Climate Change on Lake Tana’s Water Balance, Ethiopia

The aim of this study is to evaluate the hydrological impacts of climate change on the water balance of Lake Tana in Ethiopia. Impact assessments are by downscaled General Circulation Model (GCM) output and hydrological modeling. For A2 and B2 emission scenarios, precipitation, maximum and minimum temperature estimates from the HadCM3 GCM were used. GCM output was downscaled using the Statistical DownScaling Model (SDSM 4.2). Impact analyses were applied for three future time periods: early, mid and late 21st century. Over-lake evaporation is estimated by Hardgrave’s method, and over-lake precipitation is estimated by inverse distance weighing interpolation, whereas inflows from gauged and ungauged catchments are simulated by the HBV hydrological model. Findings indicate increases in maximum and minimum temperature on annual base for both emission scenarios. The projection of mean annual over lake precipitation for both A2 and B2 emission scenarios shows increasing pattern for 21st century in comparison to the baseline period. The increase of mean annual precipitation for A2 emission scenario is 9% (112 mm/year), 10% (125 mm/year) and 11% (137 mm/year) for the three future periods respectively. B2 emission scenario mean annual precipitation shows increase by 9% (111 mm/year), 10% (122 mm/year) and 10% (130 mm/year) respectively for the three future periods. Findings indicate consistent increases of lake storage for all three future periods for both A2 and B2 emission scenarios.

High Elevation Energy and Water Balance:the Roles of Surface Albedo and Temperature

Observation and modeling of the coupled energy and water balance is the key to understand hydrospheric and cryospheric processes at high elevation.The paper summarizes the progress to address this aspect in relation with different earth system elements,from glaciers to wetlands.The energy budget of two glaciers,i.e.Xiao Dongkemadi and Parlung No.4,was studied by means of extended field measurements and a distributed model of the coupled energy and mass balance was developed and evaluated.The need for accurate characterization of surface albedo was further documented for the entire Qinghai Tibet Plateau by numerical experiments with Weather Research and Forecast(WRF)on the sensitivity of the atmospheric boundary layer to the parameterization of land surface processes.A new approach to the calibration of a coupled distributed watershed model of the energy and water balance was demonstrated by a case study on the Heihe River Basin in northwestern China.The assimilation of land surface temperature did lead to the retrieval of critical soil and vegetation properties as the soil permeability and the canopy resistance to the exchange of vapour and carbon dioxide.The retrievals of actual Evapo-Transpiration(ET)were generated by the ETMonitor system and evaluated against eddy covariance measurements at sites spread throughout Asia.As regards glacier response to climate variability,the combined findings based on satellite data and model experiments showed that the spatial variability of surface albedo and temperature is significant and controls both glacier mass balance and flow.Experiments with both atmospheric and hydrosphere-cryosphere models documented the need and advantages of using accurate retrievals of land surface albedo to capture lan-atmosphere interactions at high elevation.

On the Downscaling of Meteorological Fields Using Recurrent Networks for Modelling the Water Balance in a Meso-Scale Catchment Area of Saxony,Germany

In this study, recurrent networks to downscale meteorological fields of the ERA-40 re-analysis dataset with focus on the meso-scale water balance were investigated. Therefore two types of recurrent neural networks were used. The first approach is a coupling between a recurrent neural network and a distributed watershed model and the second a nonlinear autoregressive with exogenous inputs (NARX) network, which directly predicted the component of the water balance. The approaches were deployed for a meso-scale catchment area in the Free State of Saxony, Germany. The results show that the coupled approach did not perform as well as the NARX network. But the meteorological output of the coupled approach already reaches an adequate quality. However the coupled model generates as input for the watershed model insufficient daily precipitation sums and not enough wet days were predicted. Hence the long-term annual cycle of the water balance could not be preserved with acceptable quality in contrary to the NARX approach. The residual storage change term indicates physical restrictions of the plausibility of the neural networks, whereas the physically based correlations among?the components of the water balance were preserved more accurately by the coupled approach.

Impacts of forestation on the annual and seasonal water balance of a tropical catchment under climate change

Background:This study aims to assess the effects of a forestation program and climate change on the annual and seasonal water balance of the Bogowonto catchment(597 km^(2))in Java,Indonesia.The catchment study is rare example in Indonesia where forestation has been applied at the catchment level.However,since the forestation program has been initiated,evaluations of the program only focus on the planting area targets,while the environmental success e.g.,impacts on the hydrological processes have never been assessed.This study used a calibrated Soil and Water Assessment Tool(SWAT)model to diagnose the isolated and combined effects of forestation and climate change on five water balance components,namely streamflow(Q),evapotranspiration(ET),surface runoff(Q_(s)),lateral flow(Q_(l))and base flow(Q_(b)).Results:The results show that from 2006 to 2019,forest cover has increased from 2.7% to 12.8% of the total area,while in the same period there was an increase in the mean annual and seasonal temperature,rainfall,and streamflow.Results of SWAT simulations show that changes in the mean annual and seasonal water balance under the forestation only scenario were relatively minor,while changes were more pronounced under the climate change only scenario.Based on the combined impacts scenario,it was observed that the effects of a larger forest area on the water balance were smaller than the effects of climate change.Conclusions:Although we found that forestation program has minor impacts compared to that of climate change on the hydrological processes in the Bogowonto catchment,seasonally,forestation activity has decreased the streamflow and surface runoff during the wet season which may reduce the risk of moderate floods.However,much attention should be paid to the way how forestation may result in severe drought events during the dry season.Finally,we urge the importance of accounting for the positive and negative effects in future forestation programs.

Investigation of Water Balance and Map Crisis Preparation in the Recent Two Decades of Kalacho Plain,Dehdasht,Iran

Considering the water restrictions,it is important to check the water balance in each area.These restrictions are even more important in arid and semi-arid regions.The purpose of this study was to study the water balance in the Kalacho plain of Kohgiluyeh city.Kalacho plain is about 15 km east of Dehdasht city.In recent years,due to the excessive extraction of lowland wells,the quality of groundwater has also declined steadily.Therefore,hydrogeological assessment and aquifer management seems necessary.Groundwater management requires understanding and functioning of the aquifer under natural conditions(firstly)and then predicting the effects of harvesting or feeding.Undoubtedly,understanding the actual behaviors of a natural system requires some research for each particular area.The average annual temperature and precipitation in the meteorological stations is about 18.1°C and 394 mm,respectively.The water crisis map of this plain has been prepared based on the data of 45 km2 area during the 1991-1992 blue water year.Groundwater Balance Range of the Kalacho Plain is selected based on available statistics and distribution of piezometers,wells and agricultural areas.According to the long-term hydrograph,the plain faces an average loss of 0.15 m and a reservoir deficit of 2.5 million m3.Also during the 20-year period,a total of 24 m of drop and 68 million cubic meters of reservoir deficits were created in the plain.By calculating the effective parameters in the general water balance equation(inputs and outputs),the plain was finally divided into five critical areas.Zones 1,2 and 3 are the most critical lowland areas in the vicinity of the hydraulic connection with the Gachsaran Fm.,and areas 4 and 5 are less critical because of the adjacent Asmari Formation,which feed on this karstic aquifer.Overall,this negative balance indicates an increasing drop in groundwater level and its reservoir deficit.

Assessing the Influence of Global Climate and Anthropogenic Activities on the Water Balance of an Andean Lake

Tropical regions along the Andean Cordillera face an uncertain future as mountain lakes and snow peaks exhibit receding trends associated with factors such as climatic precursors and local anthropogenic activities. Tota, the largest mountain lake in the Colombian Andes, exemplifies the role played by these factors on the lake's hydrologic evolution. A monthly water balance in Tota Lake was performed using available hydrological information from 1958 to 2007 to address interannual and multiannual level fluctuations associated with human activities and climatic precursors. The balance shows that net water uses fluctuated around 2 m3/s during this period with a pattern that, although constrained during years of severe decline in lake levels, is able to explain most of the multiannual decaying trend of 1.5 cm/year in the last 50 years. The lake’s naturalized levels were used to determine the influence of climate precursors on the lake evolution. Using Multichannel Singular Spectrum Analysis (M-SSA), significant five-year ENSO and 20-year PDO related quasi-oscillations were detected, explaining 54% of the variance associated with mean annual naturalized level fluctuations. ENSO is markedly in-phase with lake levels, with critical declines associated with low precipitation and high evaporation rates during El Ni?o years, whereas the PDO signal exhibits a phase opposition with lake levels, with low naturalized levels during a positive PDO phase and high levels during a negative PDO phase (an important result to consider given the current cooling trend of the PDO signal).

Changes of Terrestrial Water Storage in the Yellow River Basin Under Global Warming

The increasing temperature in the Yellow River Basin has led to a rapid rise in the melting level height,at a rate of 5.98 m yr^(-1)during the cold season,which further contributes to the transition from snowfall to rainfall patterns.Between 1979 and 2020,there has been a decrease in snowfall in the Yellow River Basin at a rate of-3.03 mm dec^(-1),while rainfall has been increasing at a rate of 1.00 mm dec^(-1).Consequently,the snowfall-to-rainfall ratio(SRR)has decreased.Snowfall directly replenishes terrestrial water storage(TWS)in solid form until it melts,while rainfall is rapidly lost through runoff and evaporation,in addition to infiltrating underground or remaining on the surface.Therefore,the decreasing SRR accelerates the depletion of water resources.According to the surface water balance equation,the reduction in precipitation and runoff,along with an increase in evaporation,results in a decrease in TWS during the cold season within the Yellow River Basin.In addition to climate change,human activities,considering the region's dense population and extensive agricultural land,also accelerate the decline of TWS.Notably,irrigation accounts for the largest proportion of water withdrawals in the Yellow River Basin(71.8%)and primarily occurs during the warm season(especially from June to August).The impact of human activities and climate change on the water cycle requires further in-depth research.

Transformation of surface water and groundwater and water balance in the agricultural irrigation area of the Manas River Basin,China

Calculation of the water balance is very important to relieve the pressure on water resources in arid agricultural irrigation areas.This research focused on irrigation water balance calculations in the Manas River Basin of the southern margin of the Junggar Basin of China,and aimed to analyze the groundwater level dynamic trend and response characteristics of the basin water cycle under water-saving irrigation measures.The surface water and groundwater coupling model of MIKE 11-Visual MODFLOW was used to simulate rainfall runoff in mountainous areas,and quantify the contribution of water balance components in the plain irrigation area.Convergence of the delayed watershed in the mountain area was obvious,and when the river runoff exceeded 200 m^(3)/s,the error in the runoff simulation was large.The water balance in the plain agricultural irrigation area was in a negative equilibrium state,and the difference between recharge and discharge was−2.65 billion m^(3).The evapotranspiration was 24.49 billion m^(3) under drip irrigation,accounting for only approximately 51%of the total discharge.The lateral discharge of the unsaturated and saturated aquifers was 15.38 billion m^(3),accounting for approximately 32%of the total discharge.The main reason for the groundwater decline in the irrigation area was closely related to the extraction of groundwater,because the amount of recharge and discharge in the natural state was approximately identical.The MIKE 11-Visual MODFLOW model produced accurate results,and the research method provided a new exploration technique to quantify the effect of water supply mode on the groundwater table.The model is suitable for the management of water resources in arid areas.

A Setup for a Scenario-Driven Water Balance at Landscape Scale—Assessment with AKWA-M~-Embedded in a Model Framework for Land-Use Planning’s Decision Support in Mountainous Southwest China

Land-cover changes cause a loss of natural vegetation in many parts of the world. In the Xishuangbanna （西双版纳） district （Yunnan （云南） Province）, rubber plantations replace tropical rainforests covering already an area of about 10% of the study area （2007）. There, land-use allocation is mostly driven by economic considerations. Thus, local planning authorities need decision support for land-use planning issues, which integrate socio-economic and ecological aspects. Within the NabanFrame, an agro-economic, ecological and social model was applied, which, altogether, interacted with a land allocation model via defined interfaces. Effects on the water cycle, ecological conditions as well as socio-economic should be considered by integrating the spatially distributed rainfali-runoff and water balance model AKWA-M in the model setup.

Climatic water balance dynamics over the last five decades in Romania＇s most arid region, Dobrogea

The analysis of a territory＇s climatic water balance dynamics against the background of climate change is a key component for increasing water resource management efficiency. The present study aims to analyse climatic water balance （CWB） dynamics in Romania＇s most arid region, Dobrogea, located in the southeast. The study covers the 1961-2009 period, and is based on annual and seasonal CWB values （mm）, provided by nine weather stations located throughout the region. The study, based on statistical and GIS techniques, is divided into two main stages, both carried out at annual and seasonal scales - trend analysis using the Mann-Kendall test, the Sen＇s slope method, and CWB value distribution type analysis. In order to identify the probabilistic types of distributions four mathematical models were identified-Pearson, Gamma, Chi-Squared and Wakeby, statistically verified with the P-P Plot, Q-Q Plot and Probabillity Difference Graph （PDG） curve tests. Thus, in terms of trends, the results showed a deficit increase especially at the northern stations, mainly for annual values （with a peak in the northeast, where CWB rates reached -3.2 mm/yr）. While general CWB declines occurred in winter, spring and summer, apparent decrease rates were found in the northern region （highest negative rates-summer, northwest, -1.4 mm/yr）. Autumn is an exception, due to overall increase rates which peaked in the southwest （2.3 mm/yr）. However, the entire trend analysis indicated a general lack of statistical significance. The distribution type histogram analysis showed that, annually and seasonally, deficit values are generally dominant （more noticeable in the northern region）, except for the winter season, mainly characterized by surplus intervals. Thus, the results suggest a climatic water deficit increase over the last five decades especially in northern Dobrogea, which signals the need for a spatial prioritization targeting a more efficient water resource management, necessary first and foremost for increasing regional agricultural system productivity.

Improvement in estimation of soil water deficit by integrating airborne imagery data into a soil water balance model

In this study,an approach that integrates airborne imagery data as inputs was used to improve the estimation of soil water deficit(SWD)for maize and sunflower grown under full and deficit irrigation treatments.The proposed model was applied to optimize the maximum total available soil water(TAWr)by minimizing the difference between a water stress coefficient ks and crop water stress index(1-CWSI).The optimal value of maximum TAWr was then used to calibrate a soil water balance model which in turn updated the estimation of soil water deficit.The estimates of SWD in the soil profile of both irrigated maize and sunflower fields were evaluated with the crop root zone SWD derived from neutron probe measurements and the FAO-56 SWD procedure.The results indicated a good agreement between the estimated SWD from the proposed approach and measured SWD for both maize and sunflower.The statistical analyses indicated that the maximum TAWr estimated from CWSI significantly improved the estimates of SWD,which reduced the mean absolute error(MAE)and root mean square error(RMSE)by 40%and 44%for maize and 22%for sunflower,compared with the FAO-56 model.The proposed procedure works better for crops under deficit irrigation condition.With the availability of higher spatial and temporal resolution airborne imagery during the growing season,the optimization procedure can be further improved.

Spatio-temporal variation of water conservation and its impact factors on the southern slope of Qilian Mountains

The ecology of Qilian Mountains has been seriously threatened by uncontrolled grazing and wasteland reclamation. This study examined the ecological changes on the southern slope of Qilian Mountains in China from the perspective of water conservation by classifying different clusters of water conservation functional areas to efficiently use limited human resources to tackle the water conservation protection problem. In this study, we used Integrate Valuation of Ecosystem Services and Tradeoffs(InVEST) model to estimate water conservation and analyzed the factors that influence the function. The results of this study include:(1) from 2000 to 2015, the water conservation of the southern slope of Qilian Mountains generally showed an increasing trend, and the total water conservation in 2015 increased by 42.18% compared with that in 2000.(2) Rainfall, fractional vegetation cover(FVC), and evapotranspiration have the most significant influence on the water conservation of the study area. Among them, water conservation is positively correlated with rainfall and FVC(P<0.05) and negatively correlated with evapotranspiration(P<0.05).(3) The importance level of water conservation functional areas gradually increases from northwest to southeast, and the region surrounding Menyuan Hui Autonomous County in the southeast of the southern slope of Qilian Mountains is the core water conservation functional area. And(4) the study area was divided into five clusters(Cluster Ⅰ–Cluster Ⅴ) of water conservation, with the areas of Clusters Ⅰ through Ⅴ accounting for 0.58%, 13.74%, 41.23%, 32.43%, and 12.01% of the whole study area, respectively.

Effective groundwater level recovery from mining reduction: Case study of Baoding and Shijiazhuang Plain area

The effective recovery of water level is a crucial measure of the success of comprehensive groundwater over-exploitation management actions in North China.However,traditional evaluation method do not directly capture the relationship between mining and other equilibrium elements.This study presents an innovative evaluation method to assess the water level recovery resulting from mining reduction based on the relationship between variation in exploitation and recharge.Firstly,the recharge variability of source and sink terms for both the base year and evaluation year is calculated and the coefficient of recharge variationβis introduced,which is then used to calculate the effective mining reduction and solve the water level recovery value caused by the effective mining reduction,and finally the water level recovery contribution by mining reduction is calculated by combining with the actual volume of mining reduction in the evaluation area.This research focuses on Baoding and Shijiazhuang Plain area,which share similar hydrogeological conditions but vary in groundwater exploitation and utilization.As the effect of groundwater level recovery with mining reduction was evaluated in these two areas as case study.In 2018,the results showed an effective water level recovery of 0.17 m and 0.13 m in the shallow groundwater of Shijiazhuang and Baoding Plain areas,respectively.The contributions of recovery from mining reduction were 76%and 57.98%for these two areas,respectively.It was notable that the water level recovery was most prominent in the foothill plain regions.From the evaluation results,it is evident that water level recovery depends not only on the intensity of groundwater mining reduction,but also on its effectiveness.The value of water level recovery alone cannot accurately indicate the intensity of mining reduction,as recharge variation significantly influences water level changes.Therefore,in practice,it is crucial to comprehensively assess the impact of mining reduction on water level recovery by combining the coefficient of recharge variation with the contribution of water level recovery from mining reduction.This integrated approach provide a more reasonable and scientifically supported basis,offering essential data support for groundwater management and conservation.To improve the accuracy and reliability of evaluation results,future work will focus on the standardizing and normalizing raw data processing.

Effect of glaciers on the annual catchment water balance within Budyko framework

The effects of catchment characteristics and climate variables on water partitioning into evapotranspiratio n and runoff can be evaluated using the Budyko framework.However,the influence of glaciers on catchment characteristics within the framework has yet been adequately investigated.Here we extend the Budyko framework and apply the elasticity method to examine the effects of glaciers on runoff between 2001 and 2010 in 25 upstream catchments of the Tarim River Basin in western China.The consideration of glacier mass balance and glacier fraction improves the performance of the Budyko framework,especially for the catchments with a high glacier fraction.We found that the catchment characteristic parameterwas strongly affected by glacier fraction,and it changes from 1.15 to 2.09 when glacier fraction decreases from 0.4191 to 0.0005.This also reflects the change in water-energy partitioning that eventually effects on evapotranspiration and runoff.We further assessed the average runoff responses to changes in precipitation,potential evapotranspiration,glacier mass balance,and glacier fraction in the 25 catchments.Although the runoff appears most sensitive to precipitation in average,its sensitivity to glacier mass balance and glacier area in fact rises from-0.07%to 0.17% and about 0-0.54%,respectively,when the glacier fraction increases from 0.0005 to 0.4191,further demonstrating the increasing influence of glaciers when the fraction becomes larger.After all,the inclusion of glacier factors in the Budyko framework allows us to understand more about the impacts and contributions of glaciers to runoff at a catchment scale.