Water balance modeling over variable time scales based on the Budyko framework – Model development and testing

Partitioning of precipitation into evapotranspiration and runoff is controlled by climate and catchment characteristics. The degree of control exerted by these factors varies with the spatial and temporal scales of processes modeled. The Budyko framework or the “limits” concept was used to model water balance at four temporal scales (mean annual, annual, monthly and daily). The method represents a top-down approach to hydrologic modeling and is expected to achieve parsimony of model parameters. Daily precipitation, potential evapotranspiration, and streamflow from 265 catchments in Australia were used. On a mean annual basis, the index of dryness defined as the ratio of potential evapotranspiration to precipitation was confirmed to be a dominant factor in determining the water balance with one model parameter. Analysis of the data, however, suggested increased model complexity is necessary on finer time scale such as monthly. In response, the Budyko framework for mean annual water balance was extended to include additional factors and this resulted in a parsimonious lumped conceptual model on shorter-time scale. The model was calibrated and tested against measured streamflow at variable time scales and showed promising results. The strengths of the model are consistent water balance relationships across different time scales, and model parsimony and robustness. As result, the model has the potential to be used to predict streamflow for ungauged catchments.     

Young runoff fractions control streamwater age and solute concentration dynamics

We introduce a new representation of coupled solute and water age dynamics at the catchment scale, which shows how the contributions of young runoff waters can be directly referenced to observed water quality patterns. The methodology stems from recent trends in hydrologic transport that acknowledge the dynamic nature of streamflow age and explores the use of water age fractions as an alternative to the mean age. The approach uses a travel time‐based transport model to compute the fractions of streamflow that are younger than some thresholds (e.g., younger than a few weeks) and compares them to observed solute concentration patterns. The method is here validated with data from the Hubbard Brook Experimental Forest during spring 2008, where we show that the presence of water younger than roughly 2 weeks, tracked using a hydrologic transport model and deuterium measurements, mimics the variation in dissolved silicon concentrations. Our approach suggests that an age–discharge relationship can be coupled to classic concentration–discharge relationship, to identify the links between transport timescales and solute concentration. Our results highlight that the younger streamflow components can be crucial for determining water quality variations and for characterizing the dominant hydrologic transport dynamics.     

Hydrological analysis of the Upper Tiber River Basin,Central Italy: a watershed modelling approach

The quantification of the various components of hydrological processes in a watershed remains a challenging topic as the hydrological system is altered by internal and external drivers. Watershed models have become essential tools to understand the behaviour of a catchment under dynamic processes. In this study, a physically based watershed model called Soil Water Assessment Tool was used to understand the hydrologic behaviour of the Upper Tiber River Basin, Central Italy. The model was successfully calibrated and validated using observed weather and flow data for the period of 1963–1970 and 1971–1978, respectively. Eighteen parameters were evaluated, and the model showed high relative sensitivity to groundwater flow parameters than the surface flow parameters. An analysis of annual hydrological water balance was performed for the entire upper Tiber watershed and selected subbasins. The overall behaviour of the watershed was represented by three categories of parameters governing surface flow, subsurface flow and whole basin response. The base flow contribution has shown that 60% of the streamflow is from shallow aquifer in the subbasins. The model evaluation statistics that evaluate the agreement between the simulated and the observed streamflow at the outlet of a watershed and other three different subbasins has shown a coefficient of determination (R²) from 0.68 to 0.81 and a Nash–Sutcliffe efficiency (E_NS) between 0.51 and 0.8 for the validation period. The components of the hydrologic cycle showed variation for dry and wet periods within the watershed for the same parameter sets. On the basis of the calibrated parameters, the model can be used for the prediction of the impact of climate and land use changes and water resources planning and management. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of riverbed incision on the hydrology of the Vietnamese Mekong Delta

The hydrogeomorphology of the Vietnamese Mekong Delta (VMD) has been significantly altered by natural and anthropogenic drivers. In this study, the spatiotemporal changes of the flow regime were examined by analysing the long-term daily, monthly, annual and extreme discharges and water levels from 1980 to 2018, supported by further investigation of the long-term annual sediment load (from the 1960s to 2015), river bathymetric data (in 1998, 2014 and 2017) and daily salinity concentration (from the 1990s to 2015) using various statistical methods and a coupled numerical model. Then, the effects of riverbed incision on the hydrology were investigated. The results show that the dry season discharge (i.e., in March–June) of the Tien River increased by up to 23% from the predam period (1980–1992) to the postdam period (1993–2018) but that the dry season water level at My Thuan decreased by up to −46%. The annual mean and monthly water levels in June at Tan Chau and in January and June–October at My Thuan in the Tien River decreased statistically, even though the respective discharges increased significantly. These decreased water levels instead of the increased discharges were attributed to the accelerated riverbed incision upstream from My Thuan, which increased by more than three times, from a mean rate of −0.16 m/year (−16.7 Mm³/year) in 1998–2014 to −0.5 m/year (−52.5 Mm³/year) in 2014–2017. This accelerated riverbed incision was likely caused by the reduction in the sediment load of the VMD (from 166.7 Mt/year in the predam period to 57.6 Mt/year in the postdam period) and increase in sand mining (from 3.9 Mm³ in 2012 to 13.43 Mm³ in 2018). Collectively, the decreased dry season water level in the Tien River is likely one of the main causes of the enhanced salinity intrusion.     

Streamflow forecasting in remote areas

Forecasts of runoff volumes are required in order to maximize the utility of water-supply sources. In remote areas where hydrologic and land-use data are sparse, forecast models are needed; such models should be conceptually rational so they can be transferred to remote watersheds where data are sparse. A series of models were calibrated for a large watershed in India. A spatially-distributed seasonally-varying model having a structure similar to the rational method was found to provide precise, unbiased estimates of 10-day streamflow volumes. The model was tested on a watershed that was not used for calibration, with the results indicating a high correlation between the observed and measured streamflow. Thus, the model should provide good estimates of streamflow volumes on other ungaged watersheds.     

Sensitivity of streamflow simulation in the Delaware River Basin to forecasted land‐cover change for 2030 and 2060

In order to simulate the potential effect of forecasted land‐cover change on streamflow and water availability, there has to be confidence that the hydrologic model used is sensitive to small changes in land cover (<10%) and that this land‐cover change exceeds the inherent uncertainty in forecasted conditions. To investigate this, a 26‐year streamflow record was simulated for 33 basins (54–928 km²) in the Delaware River Basin using three dates of land cover: the 2011 National Land‐Cover Dataset (Homer, Fry, & Barnes, 2012 ), 2030 land‐cover conditions representing median values from 101 equally‐likely forecasts, and 2060 land‐cover conditions corresponding to the same iterations used to represent 2030. Streamflow was simulated using a process‐based hydrologic model that includes both pervious and impervious methods as parameterized by three land‐cover‐based hydrologic response units (HRUs)—forested, agricultural, and developed land. Small, but significant differences in streamflow magnitude, variability, and seasonality were seen among the three time periods—2011, 2030, and 2060. Temporal differences were discernible from the range of conditions simulated with 101 equally likely forecasts for 2030. Development was co‐located with the most frequent landscape components, as characterized by topographic wetness index, resulting in a change in hydrology for each HRU, highlighting that knowing the location of disturbance is key to understanding potential streamflow changes. These results show that streamflow simulation using regional calibration that incorporates land‐cover‐based HRUs can be sensitive to relatively small changes in land‐cover and that temporal trends resulting from land‐cover change can be isolated in order to evaluate other changes that might affect water resources.     

Changes in streamflow regimes and their responses to different soil and water conservation measures in the Loess Plateau watersheds,China

Investigating the changes in streamflow regimes in response to various influencing factors contributes to our understanding of the mechanisms of hydrological processes in different watersheds and to water resource management strategies. This study examined streamflow regime changes by applying the indicators of hydrologic alteration method and eco-flow metrics to daily runoff data (1965–2016) from the Sandu, Hulu and Dali Rivers on the Chinese Loess Plateau, and then determined their responses to terracing, afforestation and damming. The Budyko water balance equation and the double mass curve method were used to separate the impacts of climate change and human activities on the mean discharge changes. The results showed that the terraced and dammed watersheds exhibited significant decreases in annual runoff. All hydrologic metrics indicated that the highest degree of hydrologic alteration was in the Sandu River watershed (terraced), where the monthly and extreme flows reduced significantly. In contrast, the annual eco-deficit increased significantly, indicating the highest reduction in streamflow among the three watersheds. The regulation of dams and reservoirs in the Dali River watershed has altered the flow regime, and obvious decreases in the maximum flow and slight increases in the minimum flow and baseflow indices were observed. In the Hulu River watershed (afforested), the monthly flow and extreme flows decreased slightly and were categorized as low-degree alteration, indicating that the long-term delayed effects of afforestation on hydrological processes. The magnitude of the eco-flow metrics varied with the alteration of annual precipitation. Climate change contributed 67.47% to the runoff reduction in the Hulu River watershed, while human activities played predominant roles in reducing runoff in the Sandu and Dali River watersheds. The findings revealed distinct patterns and causes of streamflow regime alteration due to different conservation measures, emphasizing the need to optimize the spatial allocation of measures to control soil erosion and utilize water resources on the Loess Plateau.     

Calibration of hydrological models in ungauged basins based on satellite radar altimetry observations of river water level

This paper proposes a new orientation to address the problem of hydrological model calibration in ungauged basin. Satellite radar altimetric observations of river water level at basin outlet are used to calibrate the model, as a surrogate of streamflow data. To shift the calibration objective, the hydrological model is coupled with a hydraulic model describing the relation between streamflow and water stage. The methodology is illustrated by a case study in the Upper Mississippi Basin using TOPEX/Poseidon (T/P) satellite data. The generalized likelihood uncertainty estimation (GLUE) is employed for model calibration and uncertainty analysis. We found that even without any streamflow information for regulating model behavior, the calibrated hydrological model can make fairly reasonable streamflow estimation. In order to illustrate the degree of additional uncertainty associated with shifting calibration objective and identifying its sources, the posterior distributions of hydrological parameters derived from calibration based on T/P data, streamflow data and T/P data with fixed hydraulic parameters are compared. The results show that the main source is the model parameter uncertainty. And the contribution of remote sensing data uncertainty is minor. Furthermore, the influence of removing high error satellite observations on streamflow estimation is also examined. Under the precondition of sufficient temporal coverage of calibration data, such data screening can eliminate some unrealistic parameter sets from the behavioral group. The study contributes to improve streamflow estimation in ungauged basin and evaluate the value of remote sensing in hydrological modeling. Copyright © 2011 John Wiley & Sons, Ltd.     

Distinguishing streamflow trends caused by changes in climate,forest cover,and permafrost in a large watershed in northeastern China

Understanding how rivers respond to changes in land cover, climate, and subsurface conditions is critical for sustainably managing water resources and ecosystems. In this study, long‐term hydrologic, climate, and satellite data (1973–2012) from the Upper Tahe River watershed (2359 km²) in the Da Hinggan Mountains of northeast China were analysed to quantify the relative hydrologic effects of climate variability (system input) and the combined influences of forest cover change and permafrost thaw (system characteristics) on average annual streamflow (system response) using 2 methods: the sensitivity‐based method and the Kendall–Theil robust line method. The study period was subdivided into a forest harvesting period (1973–1987), a forest stability period (1988–2001), and a forest recovery period (2002–2012). The results indicated that the combined effects of forest harvesting and permafrost thaw on streamflow (+ 47.0 mm) from the forest harvesting period to the forest stability period was approximately twice as large as the effect associated with climate variability (+20.2 mm). Similarly, from the forest stability period to the forest recovery period, the decrease in average annual streamflow attributed to the combined effects of forest recovery and permafrost thaw (?38.0 mm) was much greater than the decrease due to climate variability (?22.2 mm). A simple method was used to separate the distinct impacts of forest cover change and permafrost thaw, but distinguishing these influences is difficult due to changes in surface and subsurface hydrologic connectivity associated with permafrost thaw. The results highlight the need to consider multiple streamflow drivers in future watershed and aquatic ecosystem management. Due to the ecological and hydrological susceptibility to disturbances in the Da Hinggan Mountains, forest harvesting will likely negatively impact ecohydrological processes in this region, and the effects of forest species transition in the forest recovery process should be further investigated.     

Hydrologic indicators of hydraulic conditions that drive flow–biota relationships

Abstract

This article tests the association between streamflow alteration and the alteration of ecologically significant hydraulic environments. There has been a recent shift in environmental flow assessments to develop rapid desktop-based approaches that are applicable in a regional context. Streamflow statistics (e.g. minimum monthly flow) are often chosen to predict the impact of streamflow alteration on aquatic ecosystems. The assumption that the flow–biota relationship will be obscured by the effect of how streamflow interacts with channel morphology is often acknowledged, but not quantified. In this study, streamflow statistics are derived for 19 reaches in four river systems in Victoria, Australia. Hydraulic metrics were used to quantify ecologically significant surface flow conditions (Froude number) and the area of bench inundation, shallow and deep water. Multivariate analysis was used to investigate the correlation between streamflow statistics altered with regulation and the hydraulic metrics. It was found that streamflow statistics have a weak correlation to surface flow condition and the area of shallow water under natural streamflow conditions. The results show that hydrologic statistics have limited utility in quantifying changes in hydraulic environments. A similar magnitude of flow alteration can produce diverse hydraulic results. The confounding influence of channel morphology prevents streamflow statistics being an adequate surrogate for the assessment of hydraulic alteration. Modelling flow–biota relationships in a regional context is limited by the inadequacy of streamflow statistics to model ecologically significant hydraulic function. Improving knowledge of ecohydraulically significant hydrologic statistics will improve the effectiveness of environmental flow planning to sustain instream habitat conditions. A probabilistic approach is required to enable a risk-based approach to desktop generalization of flow–biota relations.

Editor Z.W. Kundzewicz; Guest editor M. Acreman

Citation Turner, M. and Stewardson, M., 2014. Hydrologic indicators of hydraulic conditions that drive flow–biota relationships. Hydrological Sciences Journal, 59 (3–4), 659–672.     

Evaporation in Mediterranean conditions: Estimations based on isotopic approaches at the watershed scale

In Mediterranean regions, the marked climatic seasonality and uneven precipitation distribution complicate the application of isotope mass balances to obtain meaningful basin-wide annual average evaporation rates. In the present study, a mass balance approach carried out on the Tavignanu River watershed in Corsica (France), showed unrealistic evaporation rate estimates: 10% for 2017–2018 and 1% for 2018–2019. This suggests that not only does evaporation alter the seasonal isotopic composition in the river, but that there is complex variability of the dominant water reservoirs contributing to the streamflow. Therefore, we propose a modified mass balance approach, including monthly quantifications of different water sources contributing to the river discharge. This allows the discrimination of isotopic variation occurring by evaporation from that originating by mixing processes. By applying this modified approach, we estimated evaporation rates on the Tavignanu River watershed that were in good agreement with results obtained by hydrological modelling: 40% for 2017–2018 and 46% for 2018–2019, respectively. An uncertainty analysis showed that evaporation rates obtained with the modified isotopic approach are close to those obtained with the non-modified approach. Therefore, we recommend using this modified isotope mass balance approach to estimate evaporation rates in such regions as the Mediterranean with high seasonality in hydrological processes.     

How much complexity is needed to simulate watershed streamflow and water quality? A test combining time series and hydrological models

Modelled hydrologic processes are represented in a set of numerical equations; the complexity of which can be measured by the total number of variables needed. A single dominant hydrologic process could control the hydrologic response of a watershed, and so the identification of the corresponding dominant variable(s) would aid in identifying a parsimonious model and in collecting more reliable data. By accounting for both model complexity and serial correlation in the variables, a model is used to identify the dominant variables for representing watershed scale streamflow, sediment transport and phosphorus yields. Long‐term water quantity and quality data were used to show that rainfall and non‐linear soil water storage were the dominant variables for weekly streamflow, suspended sediment and particulate phosphorus. Model accuracy did not consistently improve when other statistically significant variables were included. The results suggest that improved model performance may not justify the added model complexity. As such, identification of dominant variables would be the priority for developing parsimonious hydrologic models, especially at watershed scales. Copyright © 2013 John Wiley & Sons, Ltd.     

A stepwise-cluster forecasting approach for monthly streamflows based on climate teleconnections

In this study, a stepwise cluster forecasting (SCF) framework is proposed for monthly streamflow prediction in Xiangxi River, China. The developed SCF method can capture discrete and nonlinear relationships between explanatory and response variables. Cluster trees are generated through the SCF method to reflect complex relationships between independent (i.e. explanatory) and dependent (i.e. response) variables in the hydrologic system without determining specific linear/nonlinear functions. The developed SCF method is applied for monthly streamflow prediction in Xiangxi River based on the local meteorological records as well as some climate index. Comparison among SCF, multiple linear regression, generalized regression neural network, and least square support vector machine methods would be conducted. The results indicate that the SCF method would produce good predictions in both training and testing periods. Besides, the inherent probabilistic characteristics of the SCF predictions are further analyzed. The results obtained by SCF can presented as intervals, formulated by the minimum and maximum predictions as well as the 5 and 95 % percentile values of the predictions, which can reflect the variations in streamflow forecasts. Therefore, the developed SCF method can be applied for monthly streamflow prediction in various watersheds with complicated hydrologic processes.     

Global sensitivity analysis of hydrologic processes in major snow‐dominated mountainous river basins in Colorado

The performance of watershed models in simulating stream discharge depends on the adequate representation of important watershed processes. In snow‐dominated systems, snow, surface and subsurface hydrologic processes comprise a complex network of nonlinear interactions that influence the magnitude and timing of discharge. This study aims to identify critical processes and interactions that control discharge hydrographs in five major mountainous snow‐dominated river basins in Colorado, USA. A comprehensive watershed model (Soil and Water Assessment Tool) and a variance‐based global sensitivity analysis technique (Fourier Amplitude Sensitivity Test) were used in conjunction to identify critical models parameters and processes that they represent. Average monthly streamflow and streamflow root mean square error over a period of 20 years were used as two separate objective functions in this analysis. Examination of the sensitivity of monthly streamflow revealed the influence of parameters on flow volume, whereas the sensitivity of streamflow root mean square error also exposed the influence of parameters on the timing of the hydrographs. A stability analysis was performed to investigate the computational requirements for a robust sensitivity analysis. Results show that streamflow volume is mostly influenced by shallow subsurface processes, whereas interactions between groundwater and snow processes were the key in the timing of streamflows. A large majority of important parameters were common among all study watersheds, which underlies the prospect for regionalization of process‐based hydrologic modelling in headwater river basins in Colorado. Copyright © 2013 John Wiley & Sons, Ltd.     

Remote sensing and ground‐based measurements of evapotranspiration in an extreme cold Patagonian desert

Accurate estimates of seasonal evapotranspiration (ET) at different temporal and spatial scales are essential for understanding the biological and environmental determinants of ecosystem water balance in arid regions and the patterns of water utilization by the vegetation. For this purpose, remote sensing ET estimates of a Patagonian desert in Southern Argentina were verified with field measurements of soil evaporation and plant transpiration using an open top chamber. Root distribution and seasonal variation in soil volumetric water content were also analysed. There was a high correlation between remote sensing and field measurements of ecosystem water fluxes. A substantial amount of the annual ET occurred in spring and early summer (73.4 mm) using winter rain stored in the soil profile and resulting in water content depletion of the upper soil layers. A smaller amount of annual ET was derived from few rainfall events occurring during the mid or late summer (41.4 mm). According to remote sensing, the 92.9% of the mean annual precipitation returns to the atmosphere by transpiration or evaporation from the bare soil and by canopy interception. Only 7.1% infiltrates to soil layers deeper than 200 cm contributing to the water table recharge. Fourier time series analysis, cross‐correlation methods and multiple linear regression models were used to analyse 11 years of remote sensing data to assess determinants of water fluxes. A linear model predicts well the variables that drive complex ecosystem processes such as ET. Leaf area index and air temperature were not linearly correlated to ET because of the multiple interaction among variables resulting in time lags with ET variations and thus these two variables were not included in the linear model. Soil water content, the fraction of photosynthetic active radiation and precipitation explained 86% of the ET monthly variations. The high volumetric water content and the small seasonal variations at 200‐cm depth were probably the result of little water uptake from deeper soil horizons by roots with low hydraulic conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Spatial assessment of hydrologic alteration across the Pearl River Delta,China, and possible underlying causes

The alterations of the water level across the Pearl River Delta (PRD) were investigated using a ‘range of variability approach’ (RVA) based on monthly water level datasets extracted from 17 gauging stations. A mapping method was used to illustrate the spatial patterns in the degrees of alteration of water levels. The results indicated that more stations showing moderate and high alterations in monthly mean maximum and minimum water levels when compared with monthly maximum and minimum water levels. River channels characterized by higher alterations of water levels were observed mainly in the regions north of 22° 30′N. Alterations of water levels across the PRD were a consequence of various influencing factors. However, changed hypsography due to extensive and intensive human activities, particularly the large‐scale dredging and excavation of the river sand, may be taken as one of the major causes for the substantial hydrologic alteration. This study indicated that the river channels characterized by altered water levels are mostly those characterized by highly and moderately intensive sand dredging. The changed ratio of the streamflow between Makou and Sanshui stations, the major upstream flow control stations, also influenced the water level alterations of the Pearl River delta. The results of this study will be of great significance in water resources management and better human mitigation of the natural hazards due to the altered water level under the changing environment. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatiotemporal dynamics of water sources in a mountain river basin inferred through δ2H and δ18O of water

In mountainous river basins of the Pacific Northwest, climate models predict that winter warming will result in increased precipitation falling as rain and decreased snowpack. A detailed understanding of the spatial and temporal dynamics of water sources across river networks will help illuminate climate change impacts on river flow regimes. Because the stable isotopic composition of precipitation varies geographically, variation in surface water isotope ratios indicates the volume-weighted integration of upstream source water. We measured the stable isotope ratios of surface water samples collected in the Snoqualmie River basin in western Washington over June and September 2017 and the 2018 water year. We used ordinary least squares regression and geostatistical Spatial Stream Network models to relate surface water isotope ratios to mean watershed elevation (MWE) across seasons. Geologic and discharge data was integrated with water isotopes to create a conceptual model of streamflow generation for the Snoqualmie River. We found that surface water stable isotope ratios were lowest in the spring and highest in the dry, Mediterranean summer, but related strongly to MWE throughout the year. Low isotope ratios in spring reflect the input of snowmelt into high elevation tributaries. High summer isotope ratios suggest that groundwater is sourced from low elevation areas and recharged by winter precipitation. Overall, our results suggest that baseflow in the Snoqualmie River may be relatively resilient to predicted warming and subsequent changes to snowpack in the Pacific Northwest.     

Using stable isotopes to estimate young water fractions in a heavily regulated,tropical lowland river basin

The young water fraction of streamflow (Fyw), an important hydrological variable, has been calculated for the first time, for a monsoon-fed coastal catchment in northern Vietnam. Oxygen stable isotopes (δ¹⁸O) from six river sites in the Day River Basin (DRB) were analysed monthly, between January 2015 and December 2018. River δ¹⁸O signatures showed sine wave variability, reflecting the amount effect and tropical (dry-rainy) seasonality of the region. The δ¹⁸O composition of precipitation ranged from −12.67 to +1.68‰, with a mean value of −5.14‰, and in-streamflow signatures ranged from −11.63 to −1.37‰ with a mean of −5.02‰. Fractions of young water (Fyw) were calculated from the unweighted and flow-weighted δ¹⁸O composition of samples. Unweighted Fyw ranged between 29 ± 8% and 82 ± 21% with a mean value of 51 ± 19%, and was not significantly different from flow-weighted Fyw (range between 33 ± 25% and 92 ± 73%, mean 52 ± 36%). Both unweighted and flow-weighted Fyw were highest in the middle of stream and lowest in downstream sites, capturing the impacts of landuse changes, hydrology and human activities in the catchment. Our calculations imply that more than a half of rainwater reaches the DRB river mainstream within the first 3 months. The Fyw is much higher than the global average (of one-third) and insensitive to discharge due to the combination of a humid catchment with high rainfall, low storage capacity, flat landscape and an intensive drainage system in the DRB. Also the low discharge sensitivity of Fyw in the DRB implies that the regional hydrology is severely altered by humans.     

Streamflow response to the rapid retreat of a lake‐calving glacier

There has been increasing attention over the last decade to the potential effects of glacier retreat on downstream discharge and aquatic habitat. This study focused on streamflow variability downstream of Bridge Glacier in the southern Coast Mountains of BC between 1979 and 2014, prior to and during a period in which the glacier experienced enhanced calving and rapid retreat across a lake‐filled basin. Here we combined empirical trend detection and a conceptual‐parametric hydrological model to address the following hypotheses: (1) streamflow trends in late summer and early autumn should reflect the opposing influences of climatic warming (which would tend to increase unit‐area meltwater production) and the reduction in glacier area (which would tend to reduce the total volume of meltwater generated), and (2) winter streamflow should increase because of displacement of lake water as ice flows past the grounding line and calves into the lake basin. In relation to the first hypothesis, we found no significant trends in monthly discharge during summer. However, applying regression analysis to account for air temperature and precipitation variations, weak but statistically significant negative trends were detected for August and melt season discharge. The HBV‐EC model was applied using time‐varying glacier cover, as derived from Landsat imagery. Relative to simulations based on constant glacier extent, model results indicated that glacier recession caused a decline in mean monthly streamflow of 9% in August and 11% in September. These declines in late‐summer streamflow are consistent with the results from our empirical analysis. The second hypothesis is supported by the finding of positive trends for December, January, and February discharge. Despite the modelled declines in late‐summer mean monthly streamflow, recorded discharge data exhibited neither positive nor negative trends during the melt season, suggesting that Bridge Glacier may currently be at or close to the point of peak water. Further analysis of the impact of lake‐terminating glaciers on downstream discharge is needed to refine the peak water model. Copyright © 2016 John Wiley & Sons, Ltd.