Large-scale Hydrological Models Research Articles

Large-Scale Hydrological Modeling and Decision-Making for Agricultural Water Consumption and Allocation in the Main Stem Tarim River, China

A large-scale hydrological model (MIKE HYDRO) was established for the purpose of sustainable agricultural water management in the main stem Tarim River, located in northwest China. In this arid region, agricultural water consumption and allocation management are crucial to address the conflicts among irrigation water users from upstream to downstream. The results of model calibration indicated a close correlation between simulated and observed values. Scenarios with the change on irrigation strategies and land use distributions were investigated. Irrigation scenarios revealed that the available irrigation water has significant and varying effects on the yields of different crops. Irrigation water saving could reach up to 40% in the water-saving irrigation scenario. Land use scenarios illustrated that an increase of farmland area in the lower reach gravely aggravated the water deficit, while a decrease of farmland in the upper reaches resulted in considerable benefits for all sub-catchments. A substitution of crops was also investigated, which demonstrated the potential for saving considerable amounts of irrigation water in upper and middle reaches. Overall, the results of this study provide a scientific basis for decision-making on agricultural water consumption and allocation in the study area.

Projecting future climate change effects on the extreme hydrological drought events in the Weihe River basin, China

Abstract. In this study, a framework to project the potential future climate change impacts on extreme hydrological drought events in the Weihe River basin in North China is presented. This framework includes a large-scale hydrological model driven by climate outputs from a regional climate model for historical streamflow simulations and future streamflow projections, and models for univariate drought assessment and copula-based bivariate drought analysis. It is projected by the univariate drought analysis that future climate change would lead to increased frequencies of extreme hydrological drought events with higher severity. The bivariate drought assessment using copula shows that future droughts in the same return periods as historical droughts would be potentially longer and more severe, in terms of drought duration and severity. This trend would deteriorate the hydrological drought situation in the Weihe River basin. In addition, the uncertainties associated with climate models, hydrological models, and univariate and bivariate drought analysis should be quantified in the future research to improve the reliability of this study.

Hydrological projections under climate change in the near future by RegCM4 in Southern Africa using a large-scale hydrological model

This study aims to provide model estimates of changes in hydrological elements, such as EvapoTranspiration (ET) and runoff, in Southern Africa in the near future until 2029. The climate change scenarios are projected by a high-resolution Regional Climate Model (RCM), RegCM4, which is the latest version of this model developed by the Abdus Salam International Centre for Theoretical Physics (ICTP). The hydrological projections are performed by using a large-scale hydrological model (WASMOD-D), which has been tested and customized on this region prior to this study. The results reveal that (1) the projected temperature shows an increasing tendency over Southern Africa in the near future, especially eastward of 25°E, while the precipitation changes are varying between different months and sub-regions; (2) an increase in runoff (and ET) was found in eastern part of Southern Africa, i.e. Southern Mozambique and Malawi, while a decrease was estimated across the driest region in a wide area encompassing Kalahari Desert, Namibia, southwest of South Africa and Angola; (3) the strongest climate change signals are found over humid tropical areas, i.e. north of Angola and Malawi and south of Dem Rep of Congo; and (4) large spatial and temporal variability of climate change signals is found in the near future over Southern Africa. This study presents the main results of work-package 2 (WP2) of the ‘Socioeconomic Consequences of Climate Change in Sub-equatorial Africa (SoCoCA)’ project, which is funded by the Research Council of Norway.

Simulation of soil moisture for typical plain region using the Variable Infiltration Capacity model

Abstract. Droughts have a severe impact on the development of the social economy in developed plain areas. Soil moisture is a good index, it can reasonably reflect changes in drought. In this study, Jiangsu province in the Yangtze River Plain was selected as the research region, and the VIC (Variable Infiltration Capacity) large-scale hydrological model was selected to simulate the daily soil moisture with a resolution of 0.125 × 0.125 degree from 1956 to 2009. The simulated soil moisture was verified by measured soil moisture. The results indicate that the simulated soil moisture distribution is relatively consistent for the three soil layers (0–20, 20−100 and 0–100 cm), showing an increasing trend from northwest to southeast. The simulated soil moisture anomalies agreed well with in situ observations. The simulated soil moisture data thus can be used to analyze the spatio-temporal variation of the regional water content and to provide support for drought monitoring and forecasting.

Parameter identification and global sensitivity analysis of Xin'anjiang model using meta-modeling approach

Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity analysis of hydrological model is a key step in model uncertainty quantification, which can identify the dominant parameters, reduce the model calibration uncertainty, and enhance the model optimization efficiency. There are, however, some shortcomings in classical approaches, including the long duration of time and high computation cost required to quantitatively assess the sensitivity of a multiple-parameter hydrological model. For this reason, a two-step statistical evaluation framework using global techniques is presented. It is based on (1) a screening method (Morris) for qualitative ranking of parameters, and (2) a variance-based method integrated with a meta-model for quantitative sensitivity analysis, i.e., the Sobol method integrated with the response surface model (RSMSobol). First, the Morris screening method was used to qualitatively identify the parameters' sensitivity, and then ten parameters were selected to quantify the sensitivity indices. Subsequently, the RSMSobol method was used to quantify the sensitivity, i.e., the first-order and total sensitivity indices based on the response surface model (RSM) were calculated. The RSMSobol method can not only quantify the sensitivity, but also reduce the computational cost, with good accuracy compared to the classical approaches. This approach will be effective and reliable in the global sensitivity analysis of a complex large-scale distributed hydrological model.

Large-scale hydrological modeling for calculating water stress indices: implications of improved spatiotemporal resolution, surface-groundwater differentiation, and uncertainty characterization.

Physical water scarcities can be described by water stress indices. These are often determined at an annual scale and a watershed level; however, such scales mask seasonal fluctuations and spatial heterogeneity within a watershed. In order to account for this level of detail, first and foremost, water availability estimates must be improved and refined. State-of-the-art global hydrological models such as WaterGAP and UNH/GRDC have previously been unable to reliably reflect water availability at the subbasin scale. In this study, the Soil and Water Assessment Tool (SWAT) was tested as an alternative to global models, using the case study of the Mississippi watershed. While SWAT clearly outperformed the global models at the scale of a large watershed, it was judged to be unsuitable for global scale simulations due to the high calibration efforts required. The results obtained in this study show that global assessments miss out on key aspects related to upstream/downstream relations and monthly fluctuations, which are important both for the characterization of water scarcity in the Mississippi watershed and for water footprints. Especially in arid regions, where scarcity is high, these models provide unsatisfying results.

Fractional snow-covered area parameterization over complex topography

Abstract. Fractional snow-covered area (SCA) is a key parameter in large-scale hydrological, meteorological and regional climate models. Since SCA affects albedos and surface energy balance fluxes, it is especially of interest over mountainous terrain where generally a reduced SCA is observed in large grid cells. Temporal and spatial snow distributions are, however, difficult to measure over complex topography. We therefore present a parameterization of SCA based on a new subgrid parameterization for the standard deviation of snow depth over complex topography. Highly resolved snow depth data at the peak of winter were used from two distinct climatic regions, in eastern Switzerland and in the Spanish Pyrenees. Topographic scaling parameters are derived assuming Gaussian slope characteristics. We use computationally cheap terrain parameters, namely, the correlation length of subgrid topographic features and the mean squared slope. A scale dependent analysis was performed by randomly aggregating the alpine catchments in domain sizes ranging from 50 m to 3 km. For the larger domain sizes, snow depth was predominantly normally distributed. Trends between terrain parameters and standard deviation of snow depth were similar for both climatic regions, allowing one to parameterize the standard deviation of snow depth based on terrain parameters. To make the parameterization widely applicable, we introduced the mean snow depth as a climate indicator. Assuming a normal snow distribution and spatially homogeneous melt, snow-cover depletion (SCD) curves were derived for a broad range of coefficients of variations. The most accurate closed form fit resembled an existing fractional SCA parameterization. By including the subgrid parameterization for the standard deviation of snow depth, we extended the fractional SCA parameterization for topographic influences. For all domain sizes we obtained errors lower than 10% between measured and parameterized SCA.

Impacts of future climatic and land cover changes on the hydrological regime of the Madeira River basin

Brazilian strategic interest in the Madeira River basin, one of the most important of the southern Amazon tributaries, includes the development of hydropower to satisfy the country’s growing energy needs and new waterways to boost regional trade and economic development. Because of evidences that climate change impacts the hydrological regime of rivers, the aim of this study was to assess how global climate change and regional land cover change caused by deforestation could affect the river’s hydrological regime. To achieve this goal, we calibrated a large-scale hydrological model for the period from 1970–1990 and analyzed the ability of the model to simulate the present hydrological regime when climate model simulations were used as input. Climate change projections produced by climate models were used in the hydrological model to generate scenarios with and without regional land-use and land-cover changes induced by forest conversion to pasture for the period from 2011–2099. Although results show variability among models, consensus scenarios indicated a decrease in the low-flow regime. When the simulations included forest conversion to pasture, climate change impacts on low flows were reduced in the upper basin, while, in the lower basin, discharges were affected along the whole year due to the more vigorous land-use conversion in the Brazilian region of the basin.

Comparing Large-Scale Hydrological Model Predictions with Observed Streamflow in the Pacific Northwest: Effects of Climate and Groundwater*

Abstract Assessing uncertainties in hydrologic models can improve accuracy in predicting future streamflow. Here, simulated streamflows using the Variable Infiltration Capacity (VIC) model at coarse (°) and fine (°) spatial resolutions were evaluated against observed streamflows from 217 watersheds. In particular, the adequacy of VIC simulations in groundwater- versus runoff-dominated watersheds using a range of flow metrics relevant for water supply and aquatic habitat was examined. These flow metrics were 1) total annual streamflow; 2) total fall, winter, spring, and summer season streamflows; and 3) 5th, 25th, 50th, 75th, and 95th flow percentiles. The effect of climate on model performance was also evaluated by comparing the observed and simulated streamflow sensitivities to temperature and precipitation. Model performance was evaluated using four quantitative statistics: nonparametric rank correlation ρ, normalized Nash–Sutcliffe efficiency NNSE, root-mean-square error RMSE, and percent bias PBIAS. The VIC model captured the sensitivity of streamflow for temperature better than for precipitation and was in poor agreement with the corresponding temperature and precipitation sensitivities derived from observed streamflow. The model was able to capture the hydrologic behavior of the study watersheds with reasonable accuracy. Both total streamflow and flow percentiles, however, are subject to strong systematic model bias. For example, summer streamflows were underpredicted (PBIAS = −13%) in groundwater-dominated watersheds and overpredicted (PBIAS = 48%) in runoff-dominated watersheds. Similarly, the 5th flow percentile was underpredicted (PBIAS = −51%) in groundwater-dominated watersheds and overpredicted (PBIAS = 19%) in runoff-dominated watersheds. These results provide a foundation for improving model parameterization and calibration in ungauged basins.

Assimilation of satellite data to optimize large-scale hydrological model parameters: a case study for the SWOT mission

Abstract. During the last few decades, satellite measurements have been widely used to study the continental water cycle, especially in regions where in situ measurements are not readily available. The future Surface Water and Ocean Topography (SWOT) satellite mission will deliver maps of water surface elevation (WSE) with an unprecedented resolution and provide observation of rivers wider than 100 m and water surface areas greater than approximately 250 x 250 m over continental surfaces between 78° S and 78° N. This study aims to investigate the potential of SWOT data for parameter optimization for large-scale river routing models. The method consists in applying a data assimilation approach, the extended Kalman filter (EKF) algorithm, to correct the Manning roughness coefficients of the ISBA (Interactions between Soil, Biosphere, and Atmosphere)-TRIP (Total Runoff Integrating Pathways) continental hydrologic system. Parameters such as the Manning coefficient, used within such models to describe water basin characteristics, are generally derived from geomorphological relationships, which leads to significant errors at reach and large scales. The current study focuses on the Niger Basin, a transboundary river. Since the SWOT observations are not available yet and also to assess the proposed assimilation method, the study is carried out under the framework of an observing system simulation experiment (OSSE). It is assumed that modeling errors are only due to uncertainties in the Manning coefficient. The true Manning coefficients are then supposed to be known and are used to generate synthetic SWOT observations over the period 2002–2003. The impact of the assimilation system on the Niger Basin hydrological cycle is then quantified. The optimization of the Manning coefficient using the EKF (extended Kalman filter) algorithm over an 18-month period led to a significant improvement of the river water levels. The relative bias of the water level is globally improved (a 30% reduction). The relative bias of the Manning coefficient is also reduced (40% reduction) and it converges towards an optimal value. Discharge is also improved by the assimilation, but to a lesser extent than for the water levels (7%). Moreover, the method allows for a better simulation of the occurrence and intensity of flood events in the inner delta and shows skill in simulating the maxima and minima of water storage anomalies, especially in the groundwater and the aquifer reservoirs. The application of the assimilation method in the framework of an observing system simulation experiment allows evaluating the skill of the EKF algorithm to improve hydrological model parameters and to demonstrate SWOT's promising potential for global hydrology issues. However, further studies (e.g., considering multiple error sources and the difference between synthetic and real observations) are needed to achieve the evaluation of the method.

A new parallel framework of distributed SWAT calibration

With the development of large-scale hydrologic modeling, computational efficiency is becoming more and more important. Rapid modeling and analysis are needed to deal with emergency environmental disasters. The Soil and Water Assessment Tool (SWAT) is a popular hydrologic model, which is less applied in large-scale watershed simulation because of its sequential characteristics. For improving the computational efficiency of the SWAT model, we present a new parallel processing solution for hydrologic cycle and calibration based on MPI (Message Passing Interface). We partitioned sub-basins during the processes based on a load balancing method. Then the calibration was parallelized using a master-slave scheme, in which different input parameters were allocated to different processes to run the hydrologic cycle and compute the function value. Because of the slow convergence and local optimization of the SCE-UA (Shuffled Complex Evolution-developed by University of Arizona) algorithm in SWAT calibration, a genetic algorithm (GA) is developed to optimize the calibration step. Then by dividing the default communicator into several sub-communicators, all the hydrologic cycles were parallelized in their own sub-communicators to achieve further acceleration. In this paper the results show speedups for the hydrologic cycle calculations, as well as in the optimized calibration step. In the case study, we tested the parallel hydrologic cycle by four processes, and got a speedup of 3.06. In the calibration section, after applying the GA optimization, with 10 cores, we got a speed increase of 8.0 in our GA parallel framework compared with the GA sequential calibration, which is much better than the original SWAT calibration. After the sub-communicators added, this process was speeded up even further. The study demonstrated that the GA parallel framework with multi-sub-communicators is an effective and efficient solution for the hydrologists in large scale hydrology simulations.

A hydrogeologic framework for characterizing summer streamflow sensitivity to climate warming in the Pacific Northwest, USA

Abstract. Summer streamflows in the Pacific Northwest are largely derived from melting snow and groundwater discharge. As the climate warms, diminishing snowpack and earlier snowmelt will cause reductions in summer streamflow. Most regional-scale assessments of climate change impacts on streamflow use downscaled temperature and precipitation projections from general circulation models (GCMs) coupled with large-scale hydrologic models. Here we develop and apply an analytical hydrogeologic framework for characterizing summer streamflow sensitivity to a change in the timing and magnitude of recharge in a spatially explicit fashion. In particular, we incorporate the role of deep groundwater, which large-scale hydrologic models generally fail to capture, into streamflow sensitivity assessments. We validate our analytical streamflow sensitivities against two empirical measures of sensitivity derived using historical observations of temperature, precipitation, and streamflow from 217 watersheds. In general, empirically and analytically derived streamflow sensitivity values correspond. Although the selected watersheds cover a range of hydrologic regimes (e.g., rain-dominated, mixture of rain and snow, and snow-dominated), sensitivity validation was primarily driven by the snow-dominated watersheds, which are subjected to a wider range of change in recharge timing and magnitude as a result of increased temperature. Overall, two patterns emerge from this analysis: first, areas with high streamflow sensitivity also have higher summer streamflows as compared to low-sensitivity areas. Second, the level of sensitivity and spatial extent of highly sensitive areas diminishes over time as the summer progresses. Results of this analysis point to a robust, practical, and scalable approach that can help assess risk at the landscape scale, complement the downscaling approach, be applied to any climate scenario of interest, and provide a framework to assist land and water managers in adapting to an uncertain and potentially challenging future.

Impact of climate change on hydrological extremes in the Yangtze River Basin, China

The recent (1970–1999) and future (2070–2099) climates under the SRES A1B scenario, simulated by the regional climate model RegCM4.0 driven with lateral boundary conditions from the ECHAM5 general circulation model, are utilized to force a large-scale hydrological model for assessing the hydrological response to climate changes in the Yangtze River Basin, China. The variable infiltration capacity model (VIC) is utilized to simulate various hydrological components for examining the changes in streamflow at various locations throughout the Yangtze River Basin. In the end of the twenty-first century, most of the Yangtze River Basin stands out as “hotspots” of climate change in China, with an annual temperature increase of approximately 3.5 °C, an increase of annual precipitation in North and a decrease in South. Runoff in the upper reach of Yangtze River is projected to increase throughout the year in the future, especially in spring when the increase will be approximately 30 %. Runoff from the catchments in the northern part of Yangtze River will increase by approximately 10 %, whereas that in the southern part will decrease, especially in the dry season, following precipitation changes. The frequency of extreme floods at three mainstream stations (Cuntan, Yichang, and Datong) is projected to increase significantly. The original extreme floods with return periods of 50, 20, and 10 years will change into floods with return periods of no more than 20, 10, and 5 years. The projected increase in extreme floods will have significant impacts on water resources management and flood control systems in the Yangtze River Basin.

A modelling framework to project future climate change impacts on streamflow variability and extremes in the West River, China

Abstract. In this study, a hydrological modelling framework was introduced to assess the climate change impacts on future river flow in the West River basin, China, especially on streamflow variability and extremes. The modelling framework includes a delta-change method with the quantile-mapping technique to construct future climate forcings on the basis of observed meteorological data and the downscaled climate model outputs. This method is able to retain the signals of extreme weather events, as projected by climate models, in the constructed future forcing scenarios. Fed with the historical and future forcing data, a large-scale hydrologic model (the Variable Infiltration Capacity model, VIC) was executed for streamflow simulations and projections at daily time scales. A bootstrapping resample approach was used as an indirect alternative to test the equality of means, standard deviations and the coefficients of variation for the baseline and future streamflow time series, and to assess the future changes in flood return levels. The West River basin case study confirms that the introduced modelling framework is an efficient effective tool to quantify streamflow variability and extremes in response to future climate change.

The potential of coordinated reservoir operation for flood mitigation in large basins – A case study on the Bavarian Danube using coupled hydrological–hydrodynamic models

Summary The coordinated operation of reservoirs in large-scale river basins has great potential to improve flood mitigation. However, this requires large scale hydrological models to translate the effect of reservoir operation to downstream points of interest, in a quality sufficient for the iterative development of optimized operation strategies. And, of course, it requires reservoirs large enough to make a noticeable impact. In this paper, we present and discuss several methods dealing with these prerequisites for reservoir operation using the example of three major floods in the Bavarian Danube basin (45,000 km2) and nine reservoirs therein: We start by presenting an approach for multi-criteria evaluation of model performance during floods, including aspects of local sensitivity to simulation quality. Then we investigate the potential of joint hydrologic-2d-hydrodynamic modeling to improve model performance. Based on this, we evaluate upper limits of reservoir impact under idealized conditions (perfect knowledge of future rainfall) with two methods: Detailed simulations and statistical analysis of the reservoirs’ specific retention volume. Finally, we investigate to what degree reservoir operation strategies optimized for local (downstream vicinity to the reservoir) and regional (at the Danube) points of interest are compatible. With respect to model evaluation, we found that the consideration of local sensitivities to simulation quality added valuable information not included in the other evaluation criteria (Nash–Sutcliffe efficiency and Peak timing). With respect to the second question, adding hydrodynamic models to the model chain did, contrary to our expectations, not improve simulations, despite the fact that under idealized conditions (using observed instead of simulated lateral inflow) the hydrodynamic models clearly outperformed the routing schemes of the hydrological models. Apparently, the advantages of hydrodynamic models could not be fully exploited when fed by output from hydrological models afflicted with systematic errors in volume and timing. This effect could potentially be reduced by joint calibration of the hydrological–hydrodynamic model chain. Finally, based on the combination of the simulation-based and statistical impact assessment, we identified one reservoir potentially useful for coordinated, regional flood mitigation for the Danube. While this finding is specific to our test basin, the more interesting and generally valid finding is that operation strategies optimized for local and regional flood mitigation are not necessarily mutually exclusive, sometimes they are identical, sometimes they can, due to temporal offsets, be pursued simultaneously.

A Drought Monitoring and Forecasting System for Sub-Sahara African Water Resources and Food Security

Drought is one of the leading impediments to development in Africa. Much of the continent is dependent on rain-fed agriculture, which makes it particularly susceptible to climate variability. Monitoring drought and providing timely seasonal forecasts are essential for integrated drought risk reduction. Current approaches in developing regions have generally been limited, however, in part because of unreliable monitoring networks. Operational seasonal climate forecasts are also deficient and often reliant on statistical regressions, which are unable to provide detailed information relevant for drought assessment. However, the wealth of data from satellites and recent advancements in large-scale hydrological modeling and seasonal climate model predictions have enabled the development of state-of-the-art monitoring and prediction systems that can help address many of the problems inherent to developing regions. An experimental drought monitoring and forecast system for sub-Saharan Africa is described that is...

Identification of changes in hydrological drought characteristics from a multi-GCM driven ensemble constrained by observed discharge

• Global hydrological (discharge) drought is identified in a multi-model analysis. • Models were constrained by observed discharge for robust projections. • Both low flows and drought characteristics are needed to quantify future change. • Intensification of future hydrological drought is projected in most basins. Drought severity and related socio-economic impacts are expected to increase due to climate change. To better adapt to these impacts, more knowledge on changes in future hydrological drought characteristics (e.g. frequency, duration) is needed rather than only knowledge on changes in meteorological or soil moisture drought characteristics. In this study, effects of climate change on droughts in several river basins across the globe were investigated. Downscaled and bias-corrected data from three General Circulation Models (GCMs) for the A2 emission scenario were used as forcing for large-scale models. Results from five large-scale hydrological models (GHMs) run within the EU-WATCH project were used to identify low flows and hydrological drought characteristics in the control period (1971–2000) and the future period (2071–2100). Low flows were defined by the monthly 20th percentile from discharge (Q20). The variable threshold level method was applied to determine hydrological drought characteristics. The climatology of normalized Q20 from model results for the control period was compared with the climatology of normalized Q20 from observed discharge of the Global Runoff Data Centre. An observation-constrained selection of model combinations (GHM and GCM) was made based on this comparison. Prior to the assessment of future change, the selected model combinations were evaluated against observations in the period 2001–2010 for a number of river basins. The majority of the combinations (82%) that performed sufficiently in the control period, also performed sufficiently in the period 2001–2010. With the selected model combinations, future changes in drought for each river basin were identified. In cold climates, model combinations projected a regime shift and increase in low flows between the control period and future period. Arid climates were found to become even drier in the future by all model combinations. Agreement between the combinations on future low flows was low in humid climates. Changes in hydrological drought characteristics relative to the control period did not correspond to changes in low flows in all river basins. In most basins (around 65%), drought duration and deficit were projected to increase by the majority of the selected model combinations, while a decrease in low flows was projected in less basins (around 51%). Even if low discharge (monthly Q20) was not projected to decrease for each month, droughts became more severe, for example in some basins in cold climates. This is partly caused by the use of the threshold of the control period to determine drought events in the future, which led to unintended droughts in terms of expected impacts. It is important to consider both low discharge and hydrological drought characteristics to anticipate on changes in droughts for implementation of correct adaptation measures to safeguard future water resources.

Comparison of simulated water, nitrate, and bromide transport using a Hooghoudt-based and a dynamic drainage model

For large-scale hydrological modeling, the accuracy of the models used is a trade-off with the computational requirements. The models that perform well on the daily/meter scale may not perform well when applied at the yearly/kilometer scale. We compare two models of water flow and nitrate and bromide transport in a tile drained soil. The first model is based on a 2-D grid with an explicit drain node, here called the Dynamic Drainage Model (DDM). The second and less computationally expensive model is based on an 1-D vertical discretization where the horizontal flow is included as a sink term based on the Hooghoudt theory, here called the Hooghoudt Drainage Model (HDM). Both are based on Finite Volume Method solutions to Richard's equation and to the advection-dispersion equation (ADE), and embedded within the Daisy agroecological model, which includes the nitrogen cycle. The two models are run with 10 years of weather data and three different lower-boundary conditions. Losses of water, nitrogen, and bromide to both drain pipes and deep percolation/leaching are compared between the models, at daily and yearly time scales. In no case do we find the discrepancy large enough to warrant a rejection of the use of the faster HDM instead of DDM. For the daily time scale, we find in general a higher Nash-Sutcliffe efficiency coefficient for water (0.98–1.00) than for nitrate (0.97–1.00), and the lowest for bromide (0.95–1.00). The results are explained with a low concentration gradient along the water flow pathway toward the drain.

Data-driven scale extrapolation: estimating yearly discharge for a large region by small sub-basins

Abstract. Large-scale hydrological models and land surface models are so far the only tools for assessing current and future water resources. Those models estimate discharge with large uncertainties, due to the complex interaction between climate and hydrology, the limited availability and quality of data, as well as model uncertainties. A new purely data-driven scale-extrapolation method to estimate discharge for a large region solely from selected small sub-basins, which are typically 1–2 orders of magnitude smaller than the large region, is proposed. Those small sub-basins contain sufficient information, not only on climate and land surface, but also on hydrological characteristics for the large basin. In the Baltic Sea drainage basin, best discharge estimation for the gauged area was achieved with sub-basins that cover 5% of the gauged area. There exist multiple sets of sub-basins whose climate and hydrology resemble those of the gauged area equally well. Those multiple sets estimate annual discharge for the gauged area consistently well with 6 % average error. The scale-extrapolation method is completely data-driven; therefore it does not force any modelling error into the prediction. The multiple predictions are expected to bracket the inherent variations and uncertainties of the climate and hydrology of the basin.

Towards an improved land surface scheme for prairie landscapes

The prairie region of Canada and the United States is characterized by millions of small depressions of glacial origin called prairie potholes. The transfer of surface runoff in this landscape is mainly through a “fill and spill” mechanism among neighboring potholes. While non-contributing areas, that is small internally drained basins, are common on this landscape, during wet periods these areas can become hydrologically connected to larger regional drainage systems. Accurate prediction of prairie surface runoff generation and streamflow thus requires realistic representation of the dynamic threshold-mediated nature of these contributing areas. This paper presents a new prairie surface runoff generation algorithm for land surface schemes and large scale hydrological models that conceptualizes a hydrologic unit as a combination of variable and interacting storage elements. The proposed surface runoff generation algorithm uses a probability density function to represent the spatial variation of pothole storages and assumes a unique relationship between storage and the fractional contributing area for runoff (and hence amount of direct runoff generated) within a grid cell. In this paper the parameters that define this relationship are obtained by calibration against streamflow. The model was compared to an existing hydrology-land surface scheme (HLSS) applied to a typical Canadian prairie catchment, the Assiniboine River. The existing configuration is based on the Canadian Land Surface Scheme (CLASS) and WATROF (a physically-based overland and interflow scheme). The new configuration consists of CLASS coupled with the new PDMROF model. Results showed that the proposed surface runoff generation algorithm performed better at simulating streamflow, and appears to capture the dynamic nature of contributing areas in an effective and parsimonious manner. A pilot evaluation based on 1 m LiDAR data from a small (10 km2) experimental area suggests that the shape of the modeled storage-contributing area relationship is broadly consistent with that inferred from terrain analysis, under certain simplifying assumptions. The direct identification of storage–runoff parameters from terrain analysis is an outstanding challenge, and a promising area for future research.