Prospective scenarios of the saltwater intrusion in an estuary under climate change context using Bayesian neural networks

In recent decades, saltwater intrusion over some low-lying coastal regions was deteriorated by rising sea-level and decreasing streamflow in the context of climate change. Though physically-based hydrodynamic models are the most detailed means to simulate salinity processes, they are commonly restricted by data insufficiency issues both in spatial resolution and temporal lasting. This motivates us to build a statistical model enable simulation and scenario analysis for coastal salinity change with limited observations. A Bayesian neural network (BNN) model is built hereby to simulate salinity. It offers more precise estimation compared with the conventional artificial neural network. Meanwhile, the model gives the uncertainty behaviors of the final salinity simulation which is not available for other methods. Future scenarios of salinity change are constructed and analyzed in different time periods on the basis of the validated BNN model. Results indicate that the water quality over lower Pearl River is degrading along with more significant uncertainties. Further analysis suggests that streamflow alteration has a more direct impact on salinity variations than the sea-level change does. The method allows a profound analysis of the potential influence on water quality degradation in coastal and low-lying regions in support of water management and adaptation toward global climate change.     

River flow modelling: comparison of performance and evaluation of uncertainty using data-driven models and conceptual hydrological model

Hydrological and statistical models are playing an increasing role in hydrological forecasting, particularly for river basins with data of different temporal scales. In this study, statistical models, e.g. artificial neural networks, adaptive network-based fuzzy inference system, genetic programming, least squares support vector machine, multiple linear regression, were developed, based on parametric optimization methods such as particle swarm optimization (PSO), genetic algorithm (GA), and data-preprocessing techniques such as wavelet decomposition (WD) for river flow modelling using daily streamflow data from four hydrological stations for a period of 1954–2009. These models were used for 1-, 3- and 5-day streamflow forecasting and the better model was used for uncertainty evaluation using bootstrap resampling method. Meanwhile, a simple conceptual hydrological model GR4J was used to evaluate parametric uncertainty based on generalized likelihood uncertainty estimation method. Results indicated that: (1) GA and PSO did not help improve the forecast performance of the model. However, the hybrid model with WD significantly improved the forecast performance; (2) the hybrid model with WD as a data preprocessing procedure can clarify hydrological effects of water reservoirs and can capture peak high/low flow changes; (3) Forecast accuracy of data-driven models is significantly influenced by the availability of streamflow data. More human interferences from the upper to the lower East River basin can help to introduce greater uncertainty in streamflow forecasts; (4) The structure of GR4J may introduce larger parametric uncertainty at the Longchuan station than at the Boluo station in the East river basin. This study provides a theoretical background for data-driven model-based streamflow forecasting and a comprehensive view about data and parametric uncertainty in data-scarce river basins.     

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.     

Can climate variability information constrain a hydrological model for an ungauged Costa Rican catchment?

Long‐term hydrological data are key to understanding catchment behaviour and for decision making within water management and planning. Given the lack of observed data in many regions worldwide, such as Central America, hydrological models are an alternative for reproducing historical streamflow series. Additional types of information—to locally observed discharge—can be used to constrain model parameter uncertainty for ungauged catchments. Given the strong influence that climatic large‐scale processes exert on streamflow variability in the Central American region, we explored the use of climate variability knowledge as process constraints to constrain the simulated discharge uncertainty for a Costa Rican catchment, assumed to be ungauged. To reduce model uncertainty, we first rejected parameter relationships that disagreed with our understanding of the system. Then, based on this reduced parameter space, we applied the climate‐based process constraints at long‐term, inter‐annual, and intra‐annual timescales. In the first step, we reduced the initial number of parameters by 52%, and then, we further reduced the number of parameters by 3% with the climate constraints. Finally, we compared the climate‐based constraints with a constraint based on global maps of low‐flow statistics. This latter constraint proved to be more restrictive than those based on climate variability (further reducing the number of parameters by 66% compared with 3%). Even so, the climate‐based constraints rejected inconsistent model simulations that were not rejected by the low‐flow statistics constraint. When taken all together, the constraints produced constrained simulation uncertainty bands, and the median simulated discharge followed the observed time series to a similar level as an optimized model. All the constraints were found useful in constraining model uncertainty for an—assumed to be—ungauged basin. This shows that our method is promising for modelling long‐term flow data for ungauged catchments on the Pacific side of Central America and that similar methods can be developed for ungauged basins in other regions where climate variability exerts a strong control on streamflow variability.     

Combining the SWAT model with sequential uncertainty fitting algorithm for streamflow prediction and uncertainty analysis for the Lake Dianchi Basin,China

Streams play an important role in linking the land with lakes. Nutrients released from agricultural or urban sources flow via streams to lakes, causing water quality deterioration and eutrophication. Therefore, accurate simulation of streamflow is helpful for water quality improvement in lake basins. Lake Dianchi has been listed in the ‘Three Important Lakes Restoration Act’ in China, and the degradation of its water quality has been of great concern since the 1980s. To assist environmental decision making, it is important to assess and predict hydrological processes at the basin scale. This study evaluated the performance of the soil and water assessment tool (SWAT) and the feasibility of using this model as a decision support tool for predicting streamflow in the Lake Dianchi Basin. The model was calibrated and validated using monthly observed streamflow values at three flow stations within the Lake Dianchi Basin through application of the sequential uncertainty fitting algorithm (SUFI‐2). The results of the autocalibration method for calibrating and the prediction uncertainty from different sources were also examined. Together, the p‐factor (the percentage of measured data bracketed by 95% prediction of uncertainty, or 95PPU) and the r‐factor (the average thickness of the 95PPU band divided by the standard deviation of the measured data) indicated the strength of the calibration and uncertainty analysis. The results showed that the SUFI‐2 algorithm performed better than the autocalibration method. Comparison of the SUFI‐2 algorithm and autocalibration results showed that some snowmelt factors were sensitive to model output upstream at the Panlongjiang flow station. The 95PPU captured more than 70% of the observed streamflow at the three flow stations. The corresponding p‐factors and r‐factors suggested that some flow stations had relatively large uncertainty, especially in the prediction of some peak flows. Although uncertainty existed, statistical criteria including R² and Nash–Sutcliffe efficiency were reasonably determined. The model produced a useful result and can be used for further applications. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of single and multi-site calibration techniques on hydrological model performance,parameter estimation and predictive uncertainty: a case study in the Logone catchment,Lake Chad basin

Understanding hydrological processes at catchment scale through the use of hydrological model parameters is essential for enhancing water resource management. Given the difficulty of using lump parameters to calibrate distributed catchment hydrological models in spatially heterogeneous catchments, a multiple calibration technique was adopted to enhance model calibration in this study. Different calibration techniques were used to calibrate the Soil and Water Assessment Tool (SWAT) model at different locations along the Logone river channel. These were: single-site calibration (SSC); sequential calibration (SC); and simultaneous multi-site calibration (SMSC). Results indicate that it is possible to reveal differences in hydrological behavior between the upstream and downstream parts of the catchment using different parameter values. Using all calibration techniques, model performance indicators were mostly above the minimum threshold of 0.60 and 0.65 for Nash Sutcliff Efficiency (NSE) and coefficient of determination (R ²) respectively, at both daily and monthly time-steps. Model uncertainty analysis showed that more than 60% of observed streamflow values were bracketed within the 95% prediction uncertainty (95PPU) band after calibration and validation. Furthermore, results indicated that the SC technique out-performed the other two methods (SSC and SMSC). It was also observed that although the SMSC technique uses streamflow data from all gauging stations during calibration and validation, thereby taking into account the catchment spatial variability, the choice of each calibration method will depend on the application and spatial scale of implementation of the modelling results in the catchment.     

Amplification of hydrological model uncertainties in projected climate simulations of the Upper Indus Basin: Does it matter where the water is coming from?

There is a growing appreciation of the uncertainties in the estimation of snow-melt and glacier-melt as a result of climate change in high elevation catchments. Through a detailed examination of three hydrological models in two catchments, and interpretation of results from previous studies, we observed that many variations in estimated streamflow could be explained by the selection of a best parameter set from the possible good model parameters. The importance of understanding changing glacial dynamics is critically important for our study areas in the Upper Indus Basin where Pakistan's policymakers are planning infrastructure to meet the future energy and water needs of hundreds of millions of people downstream. Yet, the effect of climate on glacial runoff and climate on snowmelt runoff is poorly understood. With the HBV model, for example, we estimated glacial melt as between 56% and 89% for the Hunza catchment. When rainfall was a scaled parameter, the models estimated glacial melt as between 20% and 100% of streamflow. These parameter sets produced wildly different projections of future climate for RCP8.5 scenarios in 2046–2075 compared to 1976–2005. Assuming no glacial shrinkage, for one climate projection, we found that the choice among good parameter sets resulted in projected values of future streamflow across a range from +54% to +125%. Parameter selection was the most significant source of uncertainty in the glaciated catchment and amplified climate model uncertainty, whereas climate model choice was more important in the rainfall dominated catchment. Although the study focuses on Pakistan, the overall conclusions are instructive for other similar regions in the world. We suggest that modellers of glaciated catchments should present results from at least the book-ends: models with low sensitivity to ice-melt and models with high sensitivity to ice-melt. This would reduce confusion among decision makers when they are faced with similar contrasting results.     

The influence of parameter distribution uncertainty on hydrological and sediment modeling: a case study of SWAT model applied to the Daning watershed of the Three Gorges Reservoir Region, China

Parameter uncertainty involved in hydrological and sediment modeling often refers to the parameter dispersion and the sensitivity of the parameter. However, a limitation of the previous studies lies in that the assignment of range and specification of probability distribution for each parameter is usually difficult and subjective. Therefore, there is great uncertainty in the process of parameter calibration and model prediction. In this study, the impact of probability parameter distribution on hydrological and sediment modeling was evaluated using a semi-distributed model—the Soil and Water Assessment Tool (SWAT) and Monte Carlo method (MC)—in the Daning River watershed of the Three Gorges Reservoir Region (TGRA), China. The classic types of parameter distribution such as uniform, normal and logarithmic normal distribution were involved in this study. Based on results, parameter probability distribution showed a diverse degree of influence on the hydrological and sediment prediction, such as the sampling size, the width of 95% confidence interval (CI), the ranking of the parameter related to uncertainty, as well as the sensitivity of the parameter on model output. It can be further inferred that model parameters presented greater uncertainty in certain regions of the primitive parameter range and parameter samples densely obtained from these regions would lead to a wider 95 CI, resulting in a more doubtful prediction. This study suggested the value of the optimized value obtained by the parameter calibration process could may also be of vital importance in selecting the probability distribution function (PDF). Such cases, where parameter value corresponds to the watershed characteristic, can be used to provide a more credible distribution for both hydrological and sediment modeling.     

Assessing the relative importance of parameter and forcing uncertainty and their interactions in conceptual hydrological model simulations

Predictions of river flow dynamics provide vital information for many aspects of water management including water resource planning, climate adaptation, and flood and drought assessments. Many of the subjective choices that modellers make including model and criteria selection can have a significant impact on the magnitude and distribution of the output uncertainty. Hydrological modellers are tasked with understanding and minimising the uncertainty surrounding streamflow predictions before communicating the overall uncertainty to decision makers. Parameter uncertainty in conceptual rainfall-runoff models has been widely investigated, and model structural uncertainty and forcing data have been receiving increasing attention. This study aimed to assess uncertainties in streamflow predictions due to forcing data and the identification of behavioural parameter sets in 31 Irish catchments. By combining stochastic rainfall ensembles and multiple parameter sets for three conceptual rainfall-runoff models, an analysis of variance model was used to decompose the total uncertainty in streamflow simulations into contributions from (i) forcing data, (ii) identification of model parameters and (iii) interactions between the two. The analysis illustrates that, for our subjective choices, hydrological model selection had a greater contribution to overall uncertainty, while performance criteria selection influenced the relative intra-annual uncertainties in streamflow predictions. Uncertainties in streamflow predictions due to the method of determining parameters were relatively lower for wetter catchments, and more evenly distributed throughout the year when the Nash-Sutcliffe Efficiency of logarithmic values of flow (lnNSE) was the evaluation criterion.     

Evaluating the contribution of multi-model combination to streamflow hindcasting by empirical and conceptual models

The contribution of multi-model combination to daily streamflow hindcasting was evaluated through the HBV (Hydrologiska Byråns Vattenbalansavdelning) and RNN (recurrent neural networks) models with 100 ensemble members generated with different initial conditions for both. In the calibration phase, the analysis showed that the HBV and RNN models with 20 members have better accuracy and require less calibration time. The combination of two models, however, did not provide significant improvements when 80 more members were added in the combination. In the validation phase, the results indicated that both HBV and RNN models with 20 members not only accurately produce reliable and stable streamflow hindcasting, but also effectively simulate the timing and the value of peak flows. From the consistency of calibration and validation results, the study provides an important contribution, namely, that ensemble size is not sensitive to the type of hydrological model in terms of streamflow hindcasting.     

Effects of conditional parameterization on performance of rainfall‐runoff model regarding hydrologic non‐stationarity

Better parameterization of a hydrological model can lead to improved streamflow prediction. This is particularly important for seasonal streamflow forecasting with the use of hydrological modelling. Considering the possible effects of hydrologic non‐stationarity, this paper examined ten parameterization schemes at 12 catchments located in three different climatic zones in east Australia. These schemes are grouped into four categories according to the period when the data are used for model calibration, i.e. calibration using data: (1) from a fixed period in the historical records; (2) from different lengths of historical records prior to prediction year; (3) from different climatic analogue years in the past; and (4) data from the individual months. Parameterization schemes were evaluated according to model efficiency in both the calibration and verification period. The results show that the calibration skill changes with the different historic periods when data are used at all catchments. Comparison of model performance between the calibration schemes indicates that it is worth calibrating the model with the use of data from each individual month for the purpose of seasonal streamflow forecasting. For the catchments in the winter‐dominant rainfall region of south‐east Australia, a more significant shift in rainfall‐runoff relationships at different periods was found. For those catchments, model calibration with the use of 20 years of data prior to the prediction year leads to a more consistent performance. Copyright © 2011 John Wiley & Sons, Ltd.     

Quantification of parametric uncertainty of ANN models with GLUE method for different streamflow dynamics

This study describes the parametric uncertainty of artificial neural networks (ANNs) by employing the generalized likelihood uncertainty estimation (GLUE) method. The ANNs are used to forecast daily streamflow for three sub-basins of the Rhine Basin (East Alpine, Main, and Mosel) having different hydrological and climatological characteristics. We have obtained prior parameter distributions from 5000 ANNs in the training period to capture the parametric uncertainty and subsequently 125,000 correlated parameter sets were generated. These parameter sets were used to quantify the uncertainty in the forecasted streamflow in the testing period using three uncertainty measures: percentage of coverage, average relative length, and average asymmetry degree. The results indicated that the highest uncertainty was obtained for the Mosel sub-basin and the lowest for the East Alpine sub-basin mainly due to hydro-climatic differences between these basins. The prediction results and uncertainty estimates of the proposed methodology were compared to the direct ensemble and bootstrap methods. The GLUE method successfully captured the observed discharges with the generated prediction intervals, especially the peak flows. It was also illustrated that uncertainty bands are sensitive to the selection of the threshold value for the Nash–Sutcliffe efficiency measure used in the GLUE method by employing the Wilcoxon–Mann–Whitney test.     

How can the uncertainty in the natural inflow regime propagate into the assessment of water resource systems?

The Canadian Rocky Mountain headwaters support the water resource systems of the Canadian Prairies. Significant variations in natural headwater contributions have been observed due to warming climate. Projecting future natural headwater flows under climate change effects, however, has large uncertainty. First, there are difficulties in climate modeling and downscaling in alpine regions. Second, streamflow modeling in mountainous areas is extremely challenging. There is therefore a need to understand the effects of uncertainty in the natural inflow regime, and in particular how this translates into uncertainty in representing the state and the outflow of water resource systems. Considering the Oldman River basin in Alberta, Canada, we synthesized different inflow regimes based on site/inter-site properties of the historical inflow regime. The water resources system was then conditioned on the synthesized inflow regimes to identify the mechanisms of error propagation from the headwater streamflows to the water allocations. The results show that the response of the water resource system to the uncertainty in the generated inflow regime depends on the system state, flow condition and the component of interest. Generally, the response of the reservoirs to the uncertainty in the estimated inflow regime is more significant in dry years, in particular during low flow conditions. The response at the system outlet is rather different, as the propagation of the headwater uncertainty is more significant during high flow conditions. Also, similar inflow estimates in terms of error and uncertainty may result in different error and uncertainty estimates in the simulated outflows; therefore, lower bias and uncertainty in estimating the regional inflow regime does not necessarily mean lower bias and uncertainty in simulating the streamflow at the outlet of the system. Our results provide improved understanding of uncertainty propagation through complex water resource systems, but also portray the need for better climate and hydrological modeling in the Rocky Mountains for improved water management in the Canadian Prairies, particularly in the face of uncertain climate futures. This will be crucial if the natural headwater inflows decline and/or the system faces drought conditions.     

Uncertainty in high and low flows due to model structure and parameter errors

This paper aims to investigate the uncertainty in simulated extreme low and high flows originating from hydrological model structure and parameters. To this end, three different rainfall-runoff models, namely GR4J, HBV and Xinanjiang, are applied to two subbasins of Qiantang River basin, eastern China. The Generalised Likelihood Uncertainty Estimation approach is used for estimating the uncertainty of the three models due to parameter values, henceforth referred as parameter uncertainty. Uncertainty in simulated extreme flows is evaluated by means of the annual maximum discharge and mean annual 7-day minimum discharge. The results show that although the models have good performance for the daily flows, the uncertainty in the extreme flows could not be neglected. The uncertainty originating from parameters is larger than uncertainty due to model structure. The parameter uncertainty of the extreme flows increases with the observed discharge. The parameter uncertainty in both the extreme high flows and the extreme low flows is the largest for the HBV model and the smallest for the Xinanjiang model. It is noted that the extreme low flows are mostly underestimated by all models with optimum parameter sets for both subbasins. The largest underestimation is from Xinanjiang model. Therefore it is not reliable enough to use only one set of the parameters to make the prediction and carrying out the uncertainty study in the extreme discharge simulation could give an overall picture for the planners.     

Catchment water storage variation with elevation

One of the most important functions of catchments is the storage of water. Catchment storage buffers meteorological extremes and interannual streamflow variability, controls the partitioning between evaporation and runoff, and influences transit times of water. Hydrogeological data to estimate storage are usually scarce and seldom available for a larger set of catchments. This study focused on storage in prealpine and alpine catchments, using a set of 21 Swiss catchments comprising different elevation ranges. Catchment storage comparisons depend on storage definitions. This study defines different types of storage including definitions of dynamic and mobile catchment storage. We then estimated dynamic storage using four methods, water balance analysis, streamflow recession analysis, calibration of a bucket‐type hydrological model Hydrologiska Byråns Vattenbalansavdelning model (HBV), and calibration of a transfer function hydrograph separation model using stable isotope observations. The HBV model allowed quantifying the contributions of snow, soil and groundwater storages compared to the dynamic catchment storage. With the transfer function hydrograph separation model both dynamic and mobile storage was estimated. Dynamic storage of one catchment estimated by the four methods differed up to one order of magnitude. Nevertheless, the storage estimates ranked similarly among the 21 catchments. The largest dynamic and mobile storage estimates were found in high‐elevation catchments. Besides snow, groundwater contributed considerably to this larger storage. Generally, we found that with increasing elevation the relative contribution to the dynamic catchment storage increased for snow, decreased for soil, but remained similar for groundwater storage.     

How informative are stream level observations in different geographic regions?

Simple runoff models with a low number of model parameters are generally able to simulate catchment runoff reasonably well, but they rely on model calibration, which makes their use in ungauged basins challenging. In a previous study it has been shown that a limited number of streamflow measurements can be quite informative for constraining runoff models. In practice, however, instead of performing such repeated flow measurements, it might be easier to install a stream level logger. Here, a dataset of 600+ gauged basins in the USA was used to study how well models perform when only stream level data, rather than streamflow data, are available. A runoff model (the HBV model) was calibrated assuming that only stream level observations were available, and the simulations were evaluated on the full observed streamflow record. The results indicate that stream level data alone can already provide surprisingly good model simulation results in humid catchments, whereas in arid catchments some form of quantitative information (e.g. a streamflow observation or a regional average value) is needed to obtain good results. These results are encouraging for hydrological observations in data scarce regions as level observations are much easier to obtain than streamflow measurements. Based on runoff modelling, it might even be possible to derive streamflow time series from the level data obtained from loggers, satellites or community‐based approaches. Copyright © 2016 John Wiley & Sons, Ltd.     

Parameter identification and uncertainty analysis for simulating streamflow in a river basin of Eastern India

Uncertainty is inherent in modelling studies. However, the quantification of uncertainties associated with a model is a challenging task, and hence, such studies are somewhat limited. As distributed or semi‐distributed hydrological models are being increasingly used these days to simulate hydrological processes, it is vital that these models should be equipped with robust calibration and uncertainty analysis techniques. The goal of the present study was to calibrate and validate the Soil and Water Assessment Tool (SWAT) model for simulating streamflow in a river basin of Eastern India, and to evaluate the performance of salient optimization techniques in quantifying uncertainties. The SWAT model for the study basin was developed and calibrated using Parameter Solution (ParaSol), Sequential Uncertainty Fitting Algorithm (SUFI‐2) and Generalized Likelihood Uncertainty Estimation (GLUE) optimization techniques. The daily observed streamflow data from 1998 to 2003 were used for model calibration, and those for 2004–2005 were used for model validation. Modelling results indicated that all the three techniques invariably yield better results for the monthly time step than for the daily time step during both calibration and validation. The model performances for the daily streamflow simulation using ParaSol and SUFI‐2 during calibration are reasonably good with a Nash–Sutcliffe efficiency and mean absolute error (MAE) of 0.88 and 9.70 m³/s for ParaSol, and 0.86 and 10.07 m³/s for SUFI‐2, respectively. The simulation results of GLUE revealed that the model simulates daily streamflow during calibration with the highest accuracy in the case of GLUE (R² = 0.88, MAE = 9.56 m³/s and root mean square error = 19.70 m³/s). The results of uncertainty analyses by SUFI‐2 and GLUE were compared in terms of parameter uncertainty. It was found that SUFI‐2 is capable of estimating uncertainties in complex hydrological models like SWAT, but it warrants sound knowledge of the parameters and their effects on the model output. On the other hand, GLUE predicts more reliable uncertainty ranges (R‐factor = 0.52 for daily calibration and 0.48 for validation) compared to SUFI‐2 (R‐factor = 0.59 for daily calibration and 0.55 for validation), though it is computationally demanding. Although both SUFI‐2 and GLUE appear to be promising techniques for the uncertainty analysis of modelling results, more and more studies in this direction are required under varying agro‐climatic conditions for assessing their generic capability. Copyright © 2015 John Wiley & Sons, Ltd.     

Sensitivity analysis and uncertainty assessment in water budgets simulated by the variable infiltration capacity model for Canadian subarctic watersheds

In this study, we evaluate uncertainties propagated through different climate data sets in seasonal and annual hydrological simulations over 10 subarctic watersheds of northern Manitoba, Canada, using the variable infiltration capacity (VIC) model. Further, we perform a comprehensive sensitivity and uncertainty analysis of the VIC model using a robust and state-of-the-art approach. The VIC model simulations utilize the recently developed variogram analysis of response surfaces (VARS) technique that requires in this application more than 6,000 model simulations for a 30-year (1981–2010) study period. The method seeks parameter sensitivity, identifies influential parameters, and showcases streamflow sensitivity to parameter uncertainty at seasonal and annual timescales. Results suggest that the Ensemble VIC simulations match observed streamflow closest, whereas global reanalysis products yield high flows (0.5–3.0 mm day⁻¹) against observations and an overestimation (10–60%) in seasonal and annual water balance terms. VIC parameters exhibit seasonal importance in VARS, and the choice of input data and performance metrics substantially affect sensitivity analysis. Uncertainty propagation due to input forcing selection in each water balance term (i.e., total runoff, soil moisture, and evapotranspiration) is examined separately to show both time and space dimensionality in available forcing data at seasonal and annual timescales. Reliable input forcing, the most influential model parameters, and the uncertainty envelope in streamflow prediction are presented for the VIC model. These results, along with some specific recommendations, are expected to assist the broader VIC modelling community and other users of VARS and land surface schemes, to enhance their modelling applications.     

Hydrogeological response to climate change in alpine hillslopes

Climate change threatens water resources in snowmelt‐dependent regions by altering the fraction of snow and rain and spurring an earlier snowmelt season. The bulk of hydrological research has focused on forecasting response in streamflow volumes and timing to a shrinking snowpack; however, the degree to which subsurface storage offsets the loss of snow storage in various alpine geologic settings, i.e. the hydrogeologic buffering capacity, is still largely unknown. We address this research need by assessing the affects of climate change on storage and runoff generation for two distinct hydrogeologic settings present in alpine systems: a low storage granitic and a greater storage volcanic hillslope. We use a physically based integrated hydrologic model fully coupled to a land surface model to run a base scenario and then three progressive warming scenarios, and account for the shifts in each component of the water budget. For hillslopes with greater water retention, the larger storage volcanic hillslope buffered streamflow volumes and timing, but at the cost of greater reductions in groundwater storage relative to the low storage granite hillslope. We found that the results were highly sensitive to the unsaturated zone retention parameters, which in the case of alpine systems can be a mix of matrix or fracture flow. The presence of fractures and thus less retention in the unsaturated zone significantly decreased the reduction in recharge and runoff for the volcanic hillslope in climate warming scenarios. This approach highlights the importance of incorporating physically based subsurface flow in to alpine hydrology models, and our findings provide ways forward to arrive at a conceptual model that is both consistent with geology and hydrologic principles. Copyright © 2016 John Wiley & Sons, Ltd.