We study the controls on open water evaporation of a large lowland reservoir in the Netherlands. To this end, we analyse the dynamics of open water evaporation at two locations, Stavoren and Trintelhaven, at the border of Lake IJssel (1100ĝ€¯km2); eddy covariance systems were installed at these locations during the summer seasons of 2019 and 2020. These measurements were used to develop data-driven models for both locations. Such a statistical model is a clean and simple approach that can provide a direct indication of (and insight into) the most relevant input parameters involved in explaining the variance in open water evaporation, without making a priori assumptions regarding the process itself. We found that a combination of wind speed and the vertical vapour pressure gradient can explain most of the variability in observed hourly open water evaporation. This is in agreement with Dalton's model, which is a well-established model often used in oceanographic studies for calculating open water evaporation. Validation of the data-driven models demonstrates that a simple model using only two variables yields satisfactory results at Stavoren, with R2 values of 0.84 and 0.78 for hourly and daily data respectively. However, the validation results for Trintelhaven fall short, with R2 values of 0.67 and 0.65 for hourly and daily data respectively. Validation of the simple models that only use routinely measured meteorological variables shows adequate performance at hourly (R2Combining double low line0.78 at Stavoren and R2Combining double low line0.51 at Trintelhaven) and daily (R2Combining double low line0.82 at Stavoren and R2Combining double low line0.87 at Trintelhaven) timescales. These results for the summer periods show that open water evaporation is not directly coupled to global radiation at the hourly or daily timescale. Rather a combination of wind speed and vertical gradient of vapour pressure is the main driver at these timescales. We would like to stress the importance of including the correct drivers of open water evaporation in the parametrization in hydrological models in order to adequately represent the role of evaporation in the surface-Atmosphere coupling of inland waterbodies.

@en

The representation of snow processes in most large-scale hydrological and climate models is known to introduce considerable uncertainty into the predictions and projections of water availability. During the critical snowmelt period, the main challenge in snow modeling is that net radiation is spatially highly variable for a patchy snow cover, resulting in large horizontal differences in temperatures and heat fluxes. When a wind blows over such a system, these differences can drive advection of sensible and latent heat from the snow-free areas to the snow patches, potentially enhancing the melt rates at the leading edge and increasing the variability of subgrid melt rates. To get more insight into these processes, we examine the melt along the upwind and downwind edges of a 50 m long snow patch in the Finseelvi catchment, Norway, and try to explain the observed behavior with idealized simulations of heat fluxes and air movement over patchy snow. The melt of the snow patch was monitored from 11 June until 15 June 2019 by making use of height maps obtained through the photogrammetric structure-from-motion principle. A vertical melt of 23 ± 2.0 cm was observed at the upwind edge over the course of the field campaign, whereas the downwind edge melted only 3 ± 0.4 cm. When comparing this with meteorological measurements, we estimate the turbulent heat fluxes to be responsible for 60 % to 80 % of the upwind melt, of which a significant part is caused by the latent heat flux. The melt at the downwind edge approximately matches the melt occurring due to net radiation. To better understand the dominant processes, we represented this behavior in idealized direct numerical simulations, which are based on the measurements on a single snow patch by and resemble a flat, patchy snow cover with typical snow patch sizes of 15, 30, and 60 m. Using these simulations, we found that the reduction of the vertical temperature gradient over the snow patch was the main cause of the reductions in vertical sensible heat flux over distance from the leading edge, independent of the typical snow patch size. Moreover, we observed that the sensible heat fluxes at the leading edge and the decay over distance were independent of snow patch size as well, which resulted in a 15 % and 25 % reduction in average snowmelt for, respectively, a doubling and quadrupling of the typical snow patch size. These findings lay out pathways to include the effect of highly variable turbulent heat fluxes based on the typical snow patch size in large-scale hydrological and climate models to improve snowmelt modeling.

@en

Globally, the impacts of climate change can vary across different regions. This study uses a probability framework to evaluate recent historical (1976–2016) and near-future projected (until 2049) climate change across Europe using Climate Research Unit and ensemble climate model datasets (under RCPs 2.6 and 8.5). A historical assessment shows that although the east and west of the domain experienced the largest and smallest increase in temperature, changes in precipitation are not as smooth as temperature. Results indicate that the maximum changes between distributions of the variables (temperature and precipitation) mainly occur at extreme percentiles (e.g., 10% and 90%). A group analysis of four subregions of Europe, namely east (G1), north (G2), west/south (G3), and UK/Ireland (G4), shows that G1 and G4 are expected to have the largest increase in temperature and precipitation extremes, respectively. Although maximum increases in temperature in G3 and G4 occur at larger percentiles, G1 and G2 experience maximum increases at both large and small percentiles. The maximum increase of precipitation over the study domain, however, occurs mainly at larger extremes. To quantify changes in the distribution of projection (2020–2049) relative to the historical reference (1976–2005), two measures are defined, namely a change in occurrences (KS statistic) and intensities at different quantiles (Δ). Results confirm that the temperature distribution tends to shift to higher temperatures. Changes in distribution and extremes of precipitation are spatially variable. Furthermore, extremes are expected to be intensified under RCP 8.5. The quantile analysis and defined measures reveal changes in the entire probability distribution, reflecting possible climate changes in the future.

@en

Developments such as climate change and a growing demand for drinking water threaten the sustainability of drinking water supply worldwide. To deal with this threat, adaptation of drinking water supply systems is imperative, not only on a global and national scale but particularly on a local scale. This investigation sought to establish characteristics that describe the sustainability of local drinking water supply. The hypothesis of this research was that sustainability characteristics depend on the context that is analysed, and therefore, a variety of cases must be analysed to reach a better understanding of the sustainability of drinking water supply in the Netherlands. Therefore, three divergent cases on drinking water supply in the Netherlands were analysed. One case related to a short-term development (2018 summer drought), and two concerned long-term phenomena (changes in water quality and growth in drinking water demand). We used an integrated systems approach, describing the local drinking water supply system in terms of hydrological, technical, and socio-economic characteristics that determine the sustainability of a local drinking water supply system. To gain a perspective on the case study findings that are broader than the Dutch context, the sustainability aspects identified were paired with global aspects concerning sustainable drinking water supply. This resulted in the following set of hydrological, technical, and socio-economic sustainability characteristics: (1) water quality, water resource availability, and impact of drinking water abstraction; (2) reliability and resilience of the technical system and energy use and environmental impact; (3) drinking water availability, water governance, and land and water use. Elaboration of these sustainability characteristics and criteria into a sustainability assessment can provide information on the challenges and trade-offs inherent in the sustainable development and management of a local drinking water supply system.@en

Hydrological models are being applied for impact assessment across a wide range of resolutions. In this study, we quantify the effect of model resolution on the simulated hydrological response in five mesoscale basins in the Swiss Alps using the distributed hydrological model Spatial Processes in Hydrology (SPHY). We introduce a new metric to compare a range of values resulting from a distributed model with a single value: the density-weighted distance (DWD). Model simulations are performed at two different spatial resolutions, matching common practices in hydrology: 500m×500m matching regional-scale models, and 40km×40km matching global-scale modeling. We investigate both the intra-basin response in seasonal streamflow and evapotranspiration from the high-resolution model and the difference induced by the two different spatial resolutions, with a focus on four seasonal extremes, selected based on temperature and precipitation. Results from the high-resolution model show that the intra-basin response covers a surprisingly large range of anomalies and show that it is not uncommon to have both extreme positive and negative flux anomalies occurring simultaneously within a catchment. The intra-basin response was grouped by land cover, where different dominant runoff-generating processes are driving the differences between these groups. The low-resolution model failed to capture the diverse and contrasting response from the high-resolution model, since neither the complex topography nor land cover classes were properly represented. DWD values show that, locally, the hydrological response simulated with a high-resolution model can be a lot more extreme than a low-resolution model might indicate, which has important implications for global or continental scale assessments carried out at coarse grids of 0:5°×0:5° or 0:25°× 0:25° resolution.

@en

It is generally acknowledged in the environmental sciences that the choice of a computational model impacts the research results. In this study of a flood and drought event in the Swiss Thur basin, we show that modeling decisions during the model configuration, beyond the model choice, also impact the model results. In our carefully designed experiment we investigated four modeling decisions in ten nested basins: the spatial resolution of the model, the spatial representation of the forcing data, the calibration period, and the performance metric. The flood characteristics were mainly affected by the performance metric, whereas the drought characteristics were mainly affected by the calibration period. The results could be related to the processes that triggered the particular events studied. The impact of the modeling decisions on the simulations did, however, vary among the investigated sub-basins. In spite of the limitations of this study, our findings have important implications for the understanding and quantification of uncertainty in any hydrological or even environmental model. Modeling decisions during model configuration introduce subjectivity from the modeler. Multiple working hypotheses during model configuration can provide insights on the impact of such subjective modeling decisions.

@en

Hydrologic projections are of vital socio-economic importance. However, they are also prone to uncertainty. In order to establish a meaningful range of storylines to support water managers in decision making, we need to reveal the relevant sources of uncertainty. Here, we systematically and extensively investigate uncertainty in hydrologic projections for 605 basins throughout the contiguous US. We show that in the majority of the basins, the sign of change in average annual runoff and discharge timing for the period 2070-2100 compared to 1985-2008 differs among combinations of climate models, hydrologic models, and parameters. Mapping the results revealed that different sources of uncertainty dominate in different regions. Hydrologic model induced uncertainty in the sign of change in mean runoff was related to snow processes and aridity, whereas uncertainty in both mean runoff and discharge timing induced by the climate models was related to disagreement among the models regarding the change in precipitation. Overall, disagreement on the sign of change was more widespread for the mean runoff than for the discharge timing. The results demonstrate the need to define a wide range of quantitative hydrologic storylines, including parameter, hydrologic model, and climate model forcing uncertainty, to support water resource planning.

@en

The 6.5-km² Hupsel Brook catchment has been used as an example rural lowland area to understand rainfall–runoff processes, land–atmosphere interactions, and solute transport and to investigate how they are affected by changes in land use and water management. Meteorological and hydrological variables have been measured nearly continuously since 1964, including the 1976 drought and 2010 flood. In addition, more than 2200 water quality samples have been analyzed since the 1980s, with dedicated field campaigns focused on soil physics, evapotranspiration, and rainfall measurement. Novel insights based on these observations include the conclusion that shallow groundwater tables result in a coupled saturated–unsaturated zone and sustain evaporation in dry periods. Partitioning of rainwater between various flow paths is storage driven and therefore catchment wetness determines, together with groundwater–surface water interaction, the response of runoff to rainfall as well as solute transport.

@en

Climate change and other future developments can influence the availability of groundwater resources for drinking water. The uncertainty of the projected impact is a challenge given the urgency to decide on adaptation measures to secure the drinking water supply. Improved understanding on how climate change affects the groundwater system is necessary to develop adaptation strategies. AZURE is used, a detailed, well-calibrated hydrological model to study the projected impact of climate change scenarios on the large Veluwe aquifer in the Netherlands. The Veluwe area is an important source of drinking water. However, some existing groundwater extractions in the area affect nearby groundwater-dependent ecosystems. Redistribution of the licensed extraction volumes of these sites is considered to reduce the impact on these ecosystems. The projected impact of climate change and redistribution to groundwater levels is studied. The research shows that in a slowly responding large aquifer the projected climate change may cause rising groundwater levels despite the projected increase in summer dryness. The results indicate that this impact may exceed the impact of redistribution of extraction volumes. In addition, it is shown that the combined effect strongly depends on local conditions, thus highlighting the need for high-resolution modelling.

@en

Lowlands are vulnerable to flooding due to their mild topography in often densely populated areas with high social and economic value. Moreover, multiple physical processes coincide in lowland areas, such as those involved in river-sea interactions and in merging rivers at confluences. Simultaneous occurrence of such processes can result in amplifying or attenuating effects on water levels. Our aim is to understand the mechanisms behind simultaneous occurrence of discharge waves in a river and its lowland tributaries. Here, we introduce a new way of analyzing lowland discharge and water level dynamics, by tracing individual flood waves based on dynamic time warping. We take the confluence of the Meuse River (∼ 33000km2) with the joining tributaries of the Dommel and Aa rivers as an example, especially because the January 1995 flood at this confluence was the result of the simultaneous occurrence of discharge peaks in the main stream and the tributaries and because independent observations of water levels and discharge are available for a longer period. The analysis shows that the exact timing of the arrival of discharge peaks is of little relevance because of the long duration of the average discharge wave compared to typical time lags between peaks. The discharge waves last on average 9 days, whereas the lag time between discharge peaks in the main river and the tributaries is typically 3 days. This results in backwaters that can rise up to 1.5 m over a distance of 4 km from the confluence. Thus, local measures to reduce the impact of flooding around the confluence should account for the long duration of flood peaks in the main system.@en

We discuss two definitions of reproducibility, and question if both definitions are required to be met in computational hydrological studies.@en

Wetlands are important reservoirs of water, carbon and biodiversity. They are typical landscapes of lowland regions that have high potential for water retention. However, the hydrology of these wetlands in tropical regions is often studied in isolation from the processes taking place at the catchment scale. Our main objective is to study the hydrological dynamics of one of the largest tropical rainforest regions on an island using a combination of satellite remote sensing and novel observations from dedicated field campaigns. This contribution offers a comprehensive analysis of the hydrological dynamics of two neighbouring poorly gauged tropical basins; the Kapuas basin (98g700gkm2) in West Kalimantan and the Mahakam basin (77g100gkm2) in East Kalimantan, Indonesia. Both basins are characterised by vast areas of inland lowlands. Hereby, we put specific emphasis on key hydrological variables and indicators such as discharge and flood extent. The hydroclimatological data described herein were obtained during fieldwork campaigns carried out in the Kapuas over the period 2013-2015 and in the Mahakam over the period 2008-2010. Additionally, we used the Tropical Rainfall Measuring Mission (TRMM) rainfall estimates over the period 1998-2015 to analyse the distribution of rainfall and the influence of El-Niño - Southern Oscillation. Flood occurrence maps were obtained from the analysis of the Phase Array type L-band Synthetic Aperture Radar (PALSAR) images from 2007 to 2010. Drought events were derived from time series of simulated groundwater recharge using time series of TRMM rainfall estimates, potential evapotranspiration estimates and the threshold level approach. The Kapuas and the Mahakam lake regions are vast reservoirs of water of about 1000 and 1500gkm2 that can store as much as 3 and 6.5 billiongm3 of water, respectively. These storage capacity values can be doubled considering the area of flooding under vegetation cover. Discharge time series show that backwater effects are highly influential in the wetland regions, which can be partly explained by inundation dynamics shown by flood occurrence maps obtained from PALSAR images. In contrast to their nature as wetlands, both lowland areas have frequent periods with low soil moisture conditions and low groundwater recharge. The Mahakam wetland area regularly exhibits low groundwater recharge, which may lead to prolonged drought events that can last up to 13 months. It appears that the Mahakam lowland is more vulnerable to hydrological drought, leading to more frequent fire occurrences than in the Kapuas basin.@en

In the current human-modified world, or Anthropocene, the state of water stores and fluxes has become dependent on human as well as natural processes. Water deficits (or droughts) are the result of a complex interaction between meteorological anomalies, land surface processes, and human inflows, outflows, and storage changes. Our current inability to adequately analyse and manage drought in many places points to gaps in our understanding and to inadequate data and tools. The Anthropocene requires a new framework for drought definitions and research. Drought definitions need to be revisited to explicitly include human processes driving and modifying soil moisture drought and hydrological drought development. We give recommendations for robust drought definitions to clarify timescales of drought and prevent confusion with related terms such as water scarcity and overexploitation. Additionally, our understanding and analysis of drought need to move from single driver to multiple drivers and from uni-directional to multi-directional. We identify research gaps and propose analysis approaches on (1) drivers, (2) modifiers, (3) impacts, (4) feedbacks, and (5) changing the baseline of drought in the Anthropocene. The most pressing research questions are related to the attribution of drought to its causes, to linking drought impacts to drought characteristics, and to societal adaptation and responses to drought. Example questions include i) What are the dominant drivers of drought in different parts of the world? (ii) How do human modifications of drought enhance or alleviate drought severity? (iii) How do impacts of drought depend on the physical characteristics of drought vs. the vulnerability of people or the environment? (iv) To what extent are physical and human drought processes coupled, and can feedback loops be identified and altered to lessen or mitigate drought? (v) How should we adapt our drought analysis to accommodate changes in the normal situation (i.e. what are considered normal or reference conditions) over time? Answering these questions requires exploration of qualitative and quantitative data as well as mixed modelling approaches. The challenges related to drought research and management in the Anthropocene are not unique to drought, but do require urgent attention. We give recommendations drawn from the fields of flood research, ecology, water management, and water resources studies. The framework presented here provides a holistic view on drought in the Anthropocene, which will help improve management strategies for mitigating the severity and reducing the impacts of droughts in future.

@en

The transfer of parameter sets over different temporal and spatial resolutions is common practice in many large-domain hydrological modelling studies. The degree to which parameters are transferable across temporal and spatial resolutions is an indicator of how well spatial and temporal variability is represented in the models. A large degree of transferability may well indicate a poor representation of such variability in the employed models. To investigate parameter transferability over resolution in time and space we have set up a study in which the Variable Infiltration Capacity (VIC) model for the Thur basin in Switzerland was run with four different spatial resolutions (1 km × 1 km, 5 km × 5 km, 10 km × 10 km, lumped) and evaluated for three relevant temporal resolutions (hour, day, month), both applied with uniform and distributed forcing. The model was run 3150 times using the Hierarchical Latin Hypercube Sample and the best 1 % of the runs was selected as behavioural. The overlap in behavioural sets for different spatial and temporal resolutions was used as an indicator of parameter transferability. A key result from this study is that the overlap in parameter sets for different spatial resolutions was much larger than for different temporal resolutions, also when the forcing was applied in a distributed fashion. This result suggests that it is easier to transfer parameters across different spatial resolutions than across different temporal resolutions. However, the result also indicates a substantial underestimation in the spatial variability represented in the hydrological simulations, suggesting that the high spatial transferability may occur because the current generation of large-domain models has an inadequate representation of spatial variability and hydrologic connectivity. The results presented in this paper provide a strong motivation to further investigate and substantially improve the representation of spatial and temporal variability in large-domain hydrological models.

@en

A meta-analysis on 192 peer-reviewed articles reporting on applications of the variable infiltration capacity (VIC) model in a distributed way reveals that the spatial resolution at which the model is applied has increased over the years, while the calibration and validation time interval has remained unchanged. We argue that the calibration and validation time interval should keep pace with the increase in spatial resolution in order to resolve the processes that are relevant at the applied spatial resolution. We identified six time concepts in hydrological models, which all impact the model results and conclusions. Process-based model evaluation is particularly relevant when models are applied at hyper-resolution, where stakeholders expect credible results both at a high spatial and temporal resolution.

@en

Water management in lowland areas has focused on facilitating rapid drainage over the past centuries leading to increased flows during wet periods, but also to lower groundwater tables during dry periods. This has impacted society e.g. by increased need for irrigation and subsidence of soils. Water managers in lowland areas therefore try to focus on smart water management, wherein water is retained as much as possible. This study focuses on regional effects of delaying discharge peaks for water retention, which could alter the occurrence of discharge peaks at confluences so that the discharge peaks coincide. There is currently a lack of knowledge on processes and effects of simultaneous occurrence of discharge peaks in both hydrological and hydraulic research areas. Here, we investigate the simultaneous occurrence of discharge peaks for the River Meuse and a lowland tributary, by making use of the Dynamic Time Warping method. The outcome will help water authorities with making informed decisions to reduce flood risk and helps preparing them for floodings like the unexpected flooding of January 1995.@en

A meta-analysis on 192 peer-reviewed articles reporting applications of the Variable Infiltration Capacity (VIC) model in a distributed way reveals that the spatial resolution at which the model is applied has increased over the years, while the calibration and validation time interval has remained unchanged. We argue that the calibration and validation time interval should keep pace with the increase in spatial resolution in order to resolve the processes that are relevant at the applied spatial resolution. We identified six time concepts in hydrological models, which all impact the model results and conclusions. Process-based model evaluation is particularly relevant when models are applied at hyper-resolution, where stakeholders expect credible results both at a high spatial and temporal resolution.@en

This paper presents a hybrid local-global sensitivity analysis method termed the Distributed Evaluation of Local Sensitivity Analysis (DELSA), which is used here to identify important and unimportant parameters and evaluate how model parameter importance changes as parameter values change. DELSA uses derivative-based "local" methods to obtain the distribution of parameter sensitivity across the parameter space, which promotes consideration of sensitivity analysis results in the context of simulated dynamics. This work presents DELSA, discusses how it relates to existing methods, and uses two hydrologic test cases to compare its performance with the popular global, variance-based Sobol' method. The first test case is a simple nonlinear reservoir model with two parameters. The second test case involves five alternative "bucket-style" hydrologic models with up to 14 parameters applied to a medium-sized catchment (200 km²) in the Belgian Ardennes. Results show that in both examples, Sobol' and DELSA identify similar important and unimportant parameters, with DELSA enabling more detailed insight at much lower computational cost. For example, in the real-world problem the time delay in runoff is the most important parameter in all models, but DELSA shows that for about 20% of parameter sets it is not important at all and alternative mechanisms and parameters dominate. Moreover, the time delay was identified as important in regions producing poor model fits, whereas other parameters were identified as more important in regions of the parameter space producing better model fits. The ability to understand how parameter importance varies through parameter space is critical to inform decisions about, for example, additional data collection and model development. The ability to perform such analyses with modest computational requirements provides exciting opportunities to evaluate complicated models as well as many alternative models. Key Points Sobol' and DELSA methods showed that they yield similar sensitivity results DELSA provides more information at lower calculation costs than Sobol' DELSA is recommended as effective and efficient method for parameter sensitivity

@en

Subsurface flow and storage dynamics at hillslope scale are difficult to ascertain, often in part due to a lack of sufficient high-resolution measurements and an incomplete understanding of boundary conditions, soil properties, and other environmental aspects. A continuous and extreme rainfall experiment on an artificial hillslope at Biosphere 2's Landscape Evolution Observatory (LEO) resulted in saturation excess overland flow and gully erosion in the convergent hillslope area. An array of 496 soil moisture sensors revealed a two-step saturation process. First, the downward movement of the wetting front brought soils to a relatively constant but still unsaturated moisture content. Second, soils were brought to saturated conditions from below in response to rising water tables. Convergent areas responded faster than upslope areas, due to contributions from lateral subsurface flow driven by the topography of the bottom boundary, which is comparable to impermeable bedrock in natural environments. This led to the formation of a groundwater ridge in the convergent area, triggering saturation excess runoff generation. This unique experiment demonstrates, at very high spatial and temporal resolution, the role of convergence on subsurface storage and flow dynamics. The results bring into question the representation of saturation excess overland flow in conceptual rainfall-runoff models and land-surface models, since flow is gravity-driven in many of these models and upper layers cannot become saturated from below. The results also provide a baseline to study the role of the coevolution of ecological and hydrological processes in determining landscape water dynamics during future experiments in LEO.

@en