Search

Drought is a devastating natural disaster, during which water shortage often manifests itself in the health of vegetation. Unfortunately, it is difficult to obtain high-resolution vegetation drought impact information that is spatially and temporally consistent. While remotely sensed products can provide part of this information, they often suffer from data gaps and limitations with respect to their spatial or temporal resolution. A persistent feature among remote-sensing products is the trade-off between the spatial resolution and revisit time: high temporal resolution is met with coarse spatial resolution and vice versa. Machine learning methods have been successfully applied in a wide range of remote-sensing and hydrological studies. However, global applications to resolve drought impacts on vegetation dynamics still need to be made available, as there is significant potential for such a product to aid with improved drought impact monitoring. To this end, this study predicted global vegetation dynamics based on the enhanced vegetation index (evi) and the popular Random forest (RF) regressor algorithm at 0.1°. We assessed the applicability of RF as a gap-filling and downscaling tool to generate global evi estimates that are spatially and temporally consistent. To do this, we trained an RF regressor with 0.1° evi data, using a host of features indicative of the water and energy balances experienced by vegetation, and evaluated the performance of this new product. Next, to test whether the RF is robust in terms of spatial resolution, we downscale the global evi: the model trained on 0.1° data is used to predict evi at a 0.01° resolution. The results show that the RF can capture global evi dynamics at both a 0.1° resolution (RMSE: 0.02–0.4) and at a finer 0.01° resolution (RMSE: 0.04–0.6). Overall errors were higher in the downscaled 0.01° product compared with the 0.1° product. Nevertheless, relative increases remained small, demonstrating that RF can be used to create downscaled and temporally consistent evi products. Additional error analysis revealed that errors vary spatiotemporally, with underrepresented land cover types and periods of extreme vegetation conditions having the highest errors. Finally, this model is used to produce global, spatially continuous evi products at both a 0.1 and 0.01° spatial resolution for 2003–2013 at an 8 d frequency. [ABSTRACT FROM AUTHOR]

Maintaining good surface water quality is crucial to protect ecosystem health and for safeguarding human water use activities. However, our quantitative understanding of surface water quality is mostly predicated upon observations at monitoring stations that are highly limited in space and fragmented across time. Physical models based upon pollutant emissions and subsequent routing through the hydrological network provide opportunities to overcome these shortcomings. To this end, we have developed the dynamical surface water quality model (DynQual) for simulating water temperature (Tw) and concentrations of total dissolved solids (TDS), biological oxygen demand (BOD) and fecal coliform (FC) with a daily time step and at 5 arcmin (∼ 10 km) spatial resolution. Here, we describe the main components of this new global surface water quality model and evaluate model performance against in situ water quality observations. Furthermore, we describe both the spatial patterns and temporal trends in TDS, BOD and FC concentrations for the period 1980–2019, and we also attribute the dominant contributing sectors to surface water pollution. Modelled output indicates that multi-pollutant hotspots are especially prevalent across northern India and eastern China but that surface water quality issues exist across all world regions. Trends towards water quality deterioration have been most profound in the developing world, particularly sub-Saharan Africa and South Asia. The model code is available open source (10.5281/zenodo.7932317, Jones et al., 2023), and we provide global datasets of simulated hydrology, Tw , TDS, BOD and FC at 5 arcmin resolution with a monthly time step (10.5281/zenodo.7139222, Jones et al., 2022b). These data have the potential to inform assessments in a broad range of fields, including ecological, human health and water scarcity studies. [ABSTRACT FROM AUTHOR]

The quest for hydrological hyper-resolution modelling has been on-going for more than a decade. While global hydrological models (GHMs) have seen a reduction in grid size, they have thus far never been consistently applied at a hyper-resolution (<=1 km) at the large scale. Here, we present the first application of the GHM PCR-GLOBWB at 1 km over Europe. We thoroughly evaluated simulated discharge, evaporation, soil moisture, and terrestrial water storage anomalies against long-term observations and subsequently compared results with the established 10 and 50 km resolutions of PCR-GLOBWB. Subsequently, we could assess the added value of this first hyper-resolution version of PCR-GLOBWB and assess the scale dependencies of model and forcing resolution. Eventually, these insights can help us in understanding the current challenges and opportunities from hyper-resolution models and in formulating the model and data requirements for future improvements. We found that, for most variables, epistemic uncertainty is still large, and issues with scale commensurability exist with respect to the long-term yet coarse observations used. Merely for simulated discharge, we can confidently state that model output at hyper-resolution improves over coarser resolutions due to better representation of the river network at 1 km. However, currently available observations are not yet widely available at hyper-resolution or lack a sufficiently long time series, which makes it difficult to assess the performance of the model for other variables at hyper-resolution. Here, additional model validation efforts are needed. On the model side, hyper-resolution applications require careful revisiting of model parameterization and possibly the implementation of more physical processes to be able to resemble the dynamics and spatial heterogeneity at 1 km. With this first application of PCR-GLOBWB at 1 km, we contribute to meeting the grand challenge of hyper-resolution modelling. Even though the model was only assessed at the continental scale, valuable insights could be gained which have global validity. As such, it should be seen as a modest milestone on a longer journey towards locally relevant model output. This, however, requires a community effort from domain experts, model developers, research software engineers, and data providers. [ABSTRACT FROM AUTHOR]

Drought is a devastating natural disaster, where water shortage often manifests itself in the health of vegetation. Unfortunately, it is difficult to obtain high-resolution vegetation drought impact, which is spatially and temporally consistent. While remotely sensed products can provide part of this information, they often suffer from data gaps and limitations in spatial or temporal resolutions. A persistent feature among remote sensing products is tradeoffs between spatial resolution and revisiting times, where high temporal resolution is met by coarse spatial resolution and vice verse. Machine learning methods have been successfully applied in a wide range of remote sensing and hydrological studies. However, global applications to resolve drought impacts on vegetation dynamics still need to be made available, while there is significant potential for such a product to aid improved drought impact monitoring. To this end, this study predicted global vegetation dynamics based on the Enhanced Vegetation Index (evi) and the popular Random Forest algorithm (RF) at 0.1°. We assessed the applicability of RF as a gap filling and downscaling tool to generate spatial and temporal consistent global evi estimates. To do this, we trained an RF regressor with 0.1° evi data using a host of features indicative of water and energy balances experienced by vegetation and we evaluated the performance of this new product. Next, to test whether the RF is robust in terms of spatial resolution, we downscale global evi, the model trained on 0.1° data is used to predict evi at 0.01° resolution. The results show that the RF can capture global evi dynamics at both the 0.1° (RMSE: 0.02 - 0.4) and at the finer 0.01° (RMSE: 0.04 - 0.6) resolution. Overall errors were higher in the down-scaled 0.01° compared to the 0.1° product. Yet, relative increases remained small, thus demonstrating that RF can be used to create downscaled and temporally consistent evi products. Additional error analysis reveals that errors vary spatiotemporally, with underrepresented landcover types and periods of extreme vegetation conditions having the highest errors. Finally, this model is used to produce global spatially continuous evi products at both the 0.1° and 0.01° spatial resolution for 2003-2013 at an 8-day frequency. [ABSTRACT FROM AUTHOR]

Hydrological forecasts are important for operational water management and near-future planning, even more so in light of the increased occurrences of extreme events such as floods and droughts. Having a forecasting framework, which is flexible in terms of input forcings and forecasting locations (local, regional, or national) that can deliver this information in fast and computational efficient manner, is critical. In this study, the suitability of a hybrid forecasting framework, combining data-driven approaches and seasonal (re)forecasting information from dynamical models, to predict hydrological variables was explored. Target variables include discharge and surface water levels for various stations at a national scale, with the Netherlands as the focus. Five different machine learning (ML) models, ranging from simple to more complex and trained on historical observations of discharge, precipitation, evaporation, and seawater levels, were run with seasonal (re)forecast data, including the European Flood Awareness System (EFAS) and ECMWF seasonal forecast system (SEAS5), of these driver variables in a hindcast setting. The results were evaluated using the evaluation metrics, i.e. anomaly correlation coefficient (ACC), continuous ranked probability (skill) score (CRPS and CRPSS), and Brier skill score (BSS), in comparison to a climatological reference hindcast. Aggregating the results of all stations and ML models revealed that the hindcasting framework outperformed the climatological reference forecasts by roughly 60 % for discharge predictions (80 % for surface water level predictions). Skilful prediction for the first lead month, independently of the initialization month, can be made for discharge. The skill extends up to 2–3 months for spring months due to snowmelt dynamic captured in the training phase of the model. Surface water level hindcasts showed similar skill and skilful lead times. While the different ML models showed differences in performance during a testing and training phase using historical observations, running the ML framework in a hindcast setting showed only minor differences between the models, which is attributed to the uncertainty in seasonal forecasts. However, despite being trained on historical observations, the hybrid framework used in this study shows similar skilful predictions to previous large-scale forecasting systems. With our study, we show that a hybrid framework is able to bring location-specific skilful seasonal forecast information with global seasonal forecast inputs. At the same time, our hybrid approach is flexible and fast, and as such, a hybrid framework could be adapted to make it even more interesting to water managers and their needs, for instance, as part of a fast model-predictive control framework. [ABSTRACT FROM AUTHOR]

Wildfires are recurrent natural hazards that affect terrestrial ecosystems, the carbon cycle, climate and society. They are typically hard to predict, as their exact location and occurrence are driven by a variety of factors. Identifying a selection of dominant controls can ultimately improve predictions and projections of wildfires in both the current and a future climate. Data-driven models are suitable for identification of dominant factors of complex and partly unknown processes and can both help improve process-based models and work as independent models. In this study, we applied a data-driven machine learning approach to identify dominant hydrometeorological factors determining fire occurrence over Fennoscandia and produced spatiotemporally resolved fire danger probability maps. A random forest learner was applied to predict fire danger probabilities over space and time, using a monthly (2001–2019) satellite-based fire occurrence dataset at a 0.25 ∘ spatial grid as the target variable. The final data-driven model slightly outperformed the established Canadian Forest Fire Weather Index (FWI) used for comparison. Half of the 30 potential predictors included in the study were automatically selected for the model. Shallow volumetric soil water anomaly stood out as the dominant predictor, followed by predictors related to temperature and deep volumetric soil water. Using a local fire occurrence record for Norway as target data in a separate analysis, the test set performance increased considerably. This demonstrates the potential of developing reliable data-driven models for regions with a high-quality fire occurrence record and the limitation of using satellite-based fire occurrence data in regions subject to small fires not identified by satellites. We conclude that data-driven fire danger probability models are promising, both as a tool to identify the dominant predictors and for fire danger probability mapping. The derived relationships between wildfires and the selected predictors can further be used to assess potential changes in fire danger probability under different (future) climate scenarios. [ABSTRACT FROM AUTHOR]

Maintaining good surface water quality is crucial to protect ecosystem health and for safeguarding human water use activities. Yet, our quantitative understanding of surface water quality is mostly predicated upon observations at monitoring stations that are highly limited in space and fragmented across time. Physically-based models, based upon pollutant emissions and subsequent routing through the hydrological network, provide opportunities to overcome these shortcomings. To this end, we have developed the dynamical surface water quality model (DynQual) for simulating water temperature (Tw) and concentrations of total dissolved solids (TDS), biological oxygen demand (BOD) and fecal coliform (FC) with a daily timestep and at 5 arc-minute (~10km) spatial resolution. Here, we describe the main components of this new global surface water quality model and evaluate model performance against in-situ water quality observations. Furthermore, we describe both the spatial patterns and temporal trends in TDS, BOD and FC concentrations for the period 1980-2019, also attributing the dominant contributing sectors. The model code is available open-source (https://github.com/UU-Hydro/DYNQUAL) and we provide global datasets of simulated hydrology, Tw, TDS, BOD and FC at 5 arc-minute resolution with a monthly timestep (https://doi.org/10.5281/zenodo.7139222). This data has potential to inform assessments in a broad range of fields, including ecological, human health and water scarcity studies. [ABSTRACT FROM AUTHOR]

Human actions and climate change have drastically altered river flows across the world, resulting in adverse effects on riverine ecosystems. Environmental flows (EFs) have emerged as a prominent tool for safeguarding the riverine ecosystems, but at the global scale, the assessment of EFs is associated with high uncertainty related to the hydrological data and EF methods employed. Here, we present a novel, in-depth global EF assessment using environmental flow envelopes (EFEs). Sub-basin-specific EFEs are determined for approximately 4400 sub-basins at a monthly time resolution, and their derivation considers the methodological uncertainties related to global-scale EF studies. In addition to a lower bound of discharge based on existing EF methods, we introduce an upper bound of discharge in the EFE. This upper bound enables areas to be identified where streamflow has substantially increased above natural levels. Further, instead of only showing whether EFs are violated over a time period, we quantify, for the first time, the frequency, severity, and trends of EFE violations during the recent historical period. Discharge was derived from global hydrological model outputs from the ISIMIP 2b ensemble. We use pre-industrial (1801–1860) quasi-natural discharge together with a suite of hydrological EF methods to estimate the EFEs. We then compare the EFEs with recent historical (1976–2005) discharge to assess the violations of the EFE. These violations most commonly manifest as insufficient streamflow during the low-flow season, with fewer violations during the intermediate-flow season, and only a few violations during the high-flow season. The EFE violations are widespread and occur in half of the sub-basins of the world during more than 5 % of the months between 1976 and 2005, which is double compared with the pre-industrial period. The trends in EFE violations have mainly been increasing, which will likely continue in the future with the projected hydroclimatic changes and increases in anthropogenic water use. Indications of increased upper extreme streamflow through EFE upper bound violations are relatively scarce and dispersed. Although local fine-tuning is necessary for practical applications, and further research on the coupling between quantitative discharge and riverine ecosystem responses at the global scale is required, the EFEs provide a quick and globally robust way of determining environmental flow allocations at the sub-basin scale to inform global research and policies on water resources management. [ABSTRACT FROM AUTHOR]

The quest for hydrological hyper-resolution modelling is already on-going for more than a decade. While global hydrological models (GHMs) have seen a reduction in grid size, thus far they never have been consistently applied at hyper-resolution (<= 1km) at the large scale. Here, we present the first application of the GHM PCR-GLOBWB at 1 km over Europe. We thoroughly evaluated simulated discharge, evaporation, soil moisture, and terrestrial water storage anomalies, and subsequently compared results with the 'established' 10 km 50 km resolutions of PCR-GLOBWB. Subsequently, we could assess the added value of this first hyperresolution version of PCR-GLOBWB as well as understand model and data requirements for future improvements. found that for most variables epistemic uncertainty is still large. Merely for simulated discharge we can confidently state that model output at hyper-resolution improves over coarser resolutions. This first large-scale hyper-resolution modelling attempt shows that applying a GHM consistently is by now feasible with improved data availability and computer power. Also, simulated discharge improves due to better representation of the river network at 1 km. However, currently available observations are not yet widely available at hyper-resolution or lack sufficiently long timeseries, which makes it difficult to assess the performance of the model for other variables at hyper resolution. At the model side, hyper-resolution applications require improved parameterization implementation of physical processes to be able to resemble the dynamics and spatial heterogeneity at 1 km. With this first application of PCR-GLOBWB at 1 km, we contribute to meeting the 'grand challenge' of hyperresolution modelling. As such, it should be seen as a modest milestone on a longer journey towards locally relevant model output which requires a community effort from both model developers and data providers. [ABSTRACT FROM AUTHOR]

Hydrological forecasts are important for operational water management and near future planning, even more so in light of increased occurrences of extreme events such as floods and droughts. Having a flexible forecasting framework that can deliver this information in fast and computational efficient manner is critical. In this study, the suitability of a hybrid forecasting framework, combining data-driven approaches and seasonal (re)forecasting information to predict hydrological variables was explored. Target variables include discharge and surface water levels for various stations at national scale with the Netherlands as focus. Five different ML models, ranging from simple to more complex and trained on historical observations of discharge, precipitation, evaporation and sea water levels, were run with seasonal (re)forecast data (EFAS and SEAS5) of these driver variables in a hindcast setting. The results were evaluated using the evaluation metrics Anomaly Correlation Coefficient (ACC), Continuous Ranked Probability (Skill) Score (CRPS and CRPSS), and Brier Skill Score (BSS) in comparison to a climatological reference hindcast. Aggregating results of all stations and ML models revealed that the hindcasting framework outperformed the climatological reference forecasts by roughly 60% for discharge predictions (80% for surface water level predictions). Skilful prediction for the first lead month, independently of initialization month, can be made for discharge. The skill extends up to 2-3 months for spring months due to snow melt dynamics captured in the training phase of the model. Surface water levels hindcasts showed similar skill and skilful lead times. While the different ML models showed differences in performance during a testing and training phase using historical observations, running the ML framework in a hindcast setting showed only minor differences between the models, which is attributed to the uncertainty in seasonal forecasts. However, despite being trained on historical observations, the hybrid framework used in this study shows similar skilful predictions as previous large scale forecasting systems. With our study we show that a hybrid framework is able to bring location specific skilful seasonal forecast information with global seasonal forecast inputs. At the same time our hybrid approach is flexible and fast, and as such a hybrid framework could be adapted to make it even more interesting to water managers and their needs, for instance a part of a fast model-predictive control framework. [ABSTRACT FROM AUTHOR]

Wildfires are recurrent natural hazards that affect terrestrial ecosystems, the carbon cycle, climate and society. They are typically hard to predict, as their exact location and occurrence are driven by a variety of factors. Identifying a selection of dominant controls can ultimately improve predictions and projections of wildfires in both the current and a future climate. In this study, we applied a data-driven machine learning approach to identify dominant hydrometeorological factors determining fire occurrence over Fennoscandia, and produced spatiotemporally resolved fire danger probability maps. A random forest learner was applied to predict fire danger probabilities over space and time, using a monthly 2001-2019 satellite-based fire occurrence dataset at a 0.25° spatial grid as the target variable. The final data-driven model slightly outperformed the established Canadian fire weather index (FWI) used for comparison. Half of the 30 potential predictors included in the study were automatically selected for the model. Shallow volumetric soil water anomaly stood out as the dominant predictor, followed by predictors related to temperature and deep volumetric soil water. Using a local fire occurrence record for Norway as target data in a separate analysis, the test set performance increased considerably. This improvement shows the potential of developing reliable data-driven prediction models for regions with a high quality fire occurrence record, and the limitation of using satellite-based fire occurrence data in regions subject to small fires not picked up by satellites. We conclude that data-driven fire prediction models are promising, both as a tool to identify the dominant predictors and for fire danger probability mapping. The derived relationships between wildfires and its compound predictors can further be used to assess potential changes in fire danger probability under future climate scenarios. [ABSTRACT FROM AUTHOR]

Increasing population, economic growth and changes in diet have dramatically increased the demand for food and water over the last decades. To meet increasing demands, irrigated agriculture has expanded into semi-arid areas with limited precipitation and surface water availability. This has greatly intensified the dependence of irrigated crops on groundwater withdrawal and caused a steady increase in groundwater withdrawal and groundwater depletion. One of the effects of groundwater pumping is the reduction in streamflow through capture of groundwater recharge, with detrimental effects on aquatic ecosystems. The degree to which groundwater withdrawal affects streamflow or groundwater storage depends on the nature of the groundwater–surface water interaction (GWSI). So far, analytical solutions that have been derived to calculate the impact of groundwater on streamflow depletion involve single wells and streams and do not allow the GWSI to shift from connected to disconnected, i.e. from a situation with two-way interaction to one with a one-way interaction between groundwater and surface water. Including this shift and also analysing the effects of many wells requires numerical groundwater models that are expensive to set up. Here, we introduce an analytical framework based on a simple lumped conceptual model that allows us to estimate to what extent groundwater withdrawal affects groundwater heads and streamflow at regional scales. It accounts for a shift in GWSI, calculates at which critical withdrawal rate such a shift is expected, and when it is likely to occur after withdrawal commences. It also provides estimates of streamflow depletion and which part of the groundwater withdrawal comes out of groundwater storage and which parts from a reduction in streamflow. After a local sensitivity analysis, the framework is combined with parameters and inputs from a global hydrological model and subsequently used to provide global maps of critical withdrawal rates and timing, the areas where current withdrawal exceeds critical limits and maps of groundwater and streamflow depletion rates that result from groundwater withdrawal. The resulting global depletion rates are compared with estimates from in situ observations and regional and global groundwater models and satellites. Pairing of the analytical framework with more complex global hydrological models presents a screening tool for fast first-order assessments of regional-scale groundwater sustainability and for supporting hydro-economic models that require simple relationships between groundwater withdrawal rates and the evolution of pumping costs and environmental externalities. [ABSTRACT FROM AUTHOR]

Global water models (GWMs) simulate the terrestrial water cycle on the global scale and are used to assess the impacts of climate change on freshwater systems. GWMs are developed within different modelling frameworks and consider different underlying hydrological processes, leading to varied model structures. Furthermore, the equations used to describe various processes take different forms and are generally accessible only from within the individual model codes. These factors have hindered a holistic and detailed understanding of how different models operate, yet such an understanding is crucial for explaining the results of model evaluation studies, understanding inter-model differences in their simulations, and identifying areas for future model development. This study provides a comprehensive overview of how 16 state-of-the-art GWMs are designed. We analyse water storage compartments, water flows, and human water use sectors included in models that provide simulations for the Inter-Sectoral Impact Model Intercomparison Project phase 2b (ISIMIP2b). We develop a standard writing style for the model equations to enhance model intercomparison, improvement, and communication. In this study, WaterGAP2 used the highest number of water storage compartments, 11, and CWatM used 10 compartments. Six models used six compartments, while four models (DBH, JULES-W1, Mac-PDM.20, and VIC) used the lowest number, three compartments. WaterGAP2 simulates five human water use sectors, while four models (CLM4.5, CLM5.0, LPJmL, and MPI-HM) simulate only water for the irrigation sector. We conclude that, even though hydrological processes are often based on similar equations for various processes, in the end these equations have been adjusted or models have used different values for specific parameters or specific variables. The similarities and differences found among the models analysed in this study are expected to enable us to reduce the uncertainty in multi-model ensembles, improve existing hydrological processes, and integrate new processes. [ABSTRACT FROM AUTHOR]

Soil moisture is highly variable in space and time, and deficits (i.e., droughts) play an important role in modulating crop yields. Limited hydroclimate and yield data, however, hamper drought impact monitoring and assessment at the farm field scale. This study demonstrates the potential of using field-scale soil moisture simulations to support high-resolution agricultural yield prediction and drought monitoring at the smallholder farm field scale. We present a multiscale modeling approach that combines HydroBlocks – a physically based hyper-resolution land surface model (LSM) – with machine learning. We used HydroBlocks to simulate root zone soil moisture and soil temperature in Zambia at 3 h 30 m resolution. These simulations, along with remotely sensed vegetation indices, meteorological data, and descriptors of the physical landscape (related to topography, land cover, and soils) were combined with district-level maize data to train a random forest (RF) model to predict maize yields at district and field scales (250 m). Our model predicted yields with an average testing coefficient of determination (R2) of 0.57 and mean absolute error (MAE) of 310 kgha-1 using year-based cross-validation. Our predicted maize losses due to the 2015–2016 El Niño drought agreed well with losses reported by the Food and Agriculture Organization (FAO). Our results reveal that soil moisture is the strongest and most reliable predictor of maize yield, driving its spatial and temporal variability. Soil moisture was also a more effective indicator of drought impacts on crops than precipitation, soil and air temperatures, and remotely sensed normalized difference vegetation index (NDVI)-based drought indices. This study demonstrates how field-scale modeling can help bridge the spatial-scale gap between drought monitoring and agricultural impacts. [ABSTRACT FROM AUTHOR]

Global water models (GWMs) simulate the terrestrial water cycle, on the global scale, and are used to assess the impacts of climate change on freshwater systems. GWMs are developed within different modeling frameworks and consider different underlying hydrological processes, leading to varied model structures. Furthermore, the equations used to describe various processes take different forms and are generally accessible only from within the individual model codes. These factors have hindered a holistic and detailed understanding of how different models operate, yet such an understanding is crucial for explaining the results of model evaluation studies, understanding inter-model differences in their simulations, and identifying areas for future model development. This study provides a comprehensive overview of how state-of-the-art GWMs are designed. We analyze water storage compartments, water flows, and human water use sectors included in 16 GWMs that provide simulations for the Inter-Sectoral Impact Model Intercomparison Project phase 2b (ISIMIP2b). We develop a standard writing style for the model equations to further enhance model improvement, intercomparison, and communication. In this study, WaterGAP2 used the highest number of water storage compartments, 11, and CWatM used 10 compartments. Seven models used six compartments, while three models (JULES-W1, Mac-PDM.20, and VIC) used the lowest number, three compartments. WaterGAP2 simulates five human water use sectors, while four models (CLM4.5, CLM5.0, LPJmL, and MPI-HM) simulate only water used by humans for the irrigation sector. We conclude that even though hydrologic processes are often based on similar equations, in the end, these equations have been adjusted or have used different values for specific parameters or specific variables. Our results highlight that the predictive uncertainty of GWMs can be reduced through improvements of the existing hydrologic processes, implementation of new processes in the models, and high-quality input data. [ABSTRACT FROM AUTHOR]

In eastern Africa droughts can cause crop failure and lead to food insecurity. With increasing temperatures, there is an a priori assumption that droughts are becoming more severe. However, the link between droughts and climate change is not sufficiently understood. Here we investigate trends in long-term agricultural drought and the influence of increasing temperatures and precipitation deficits. Using a combination of models and observational datasets, we studied trends, spanning the period from 1900 (to approximate pre-industrial conditions) to 2018, for six regions in eastern Africa in four drought-related annually averaged variables: soil moisture, precipitation, temperature, and evaporative demand (E0). In standardized soil moisture data, we found no discernible trends. The strongest influence on soil moisture variability was from precipitation, especially in the drier or water-limited study regions; temperature and E0 did not demonstrate strong relations to soil moisture. However, the error margins on precipitation trend estimates are large and no clear trend is evident, whereas significant positive trends were observed in local temperatures. The trends in E0 are predominantly positive, but we do not find strong relations between E0 and soil moisture trends. Nevertheless, the E0 trend results can still be of interest for irrigation purposes because it is E0 that determines the maximum evaporation rate. We conclude that until now the impact of increasing local temperatures on agricultural drought in eastern Africa is limited and we recommend that any soil moisture analysis be supplemented by an analysis of precipitation deficit. [ABSTRACT FROM AUTHOR]

Increasing population, economic growth and changes in diet have dramatically increased the demand for food and water over the last decades. To meet increasing demands, irrigated agriculture has expanded into semi-arid areas with limited precipitation and surface water availability. This has greatly intensified the dependence of irrigated crops on groundwater withdrawal and caused a steady increase of non-renewable groundwater use, i.e. groundwater taken out of aquifer storage that will not be replenished in human time scales. One of the effects of groundwater pumping is the reduction in streamflow through capture of groundwater recharge, with detrimental effects on aquatic ecosystems. The degree to which groundwater withdrawal affects streamflow or groundwater storage depends on the nature of the groundwater-surface water interaction (GWSI). So far, analytical solutions that have been derived to calculate the impact of groundwater on streamflow depletion involve single wells and streams and do not allow the GWSI to shift from connected to disconnected, i.e. from a situation with two-way interaction to one with a one-way interaction between groundwater and surface water. Including this shift and also analyse the effects of many wells, requires numerical groundwater models that are expensive to setup. Here, we introduce a simple conceptual analytical framework that allows to estimate to what extent groundwater withdrawal affects groundwater heads and streamflow. It allows for a shift in GWSI, calculates at which critical withdrawal rate such a shift is expected and when it is likely to occur after withdrawal commences. It also provides estimates of streamflow depletion and which part of the groundwater withdrawal comes out of groundwater storage and which parts from a reduction in streamflow. After a local sensitivity analysis, the framework is used to provide global maps of critical withdrawal rates and timing, the areas where current withdrawal exceeds critical limits, and maps of groundwater depletion and streamflow depletion rates that result from groundwater withdrawal. The resulting global depletion rates are similar to those obtained from global hydrological models and satellites. The analytical framework is particularly useful for performing first-order sensitivity studies and for supporting hydroeconomic models that require simple relationships between groundwater withdrawal rates and the evolution of pumping costs and environmental externalities. [ABSTRACT FROM AUTHOR]

Soil moisture is highly variable in space, and its deficits (i.e. droughts) plays an important role in modulating crop yields and its variability across landscapes. Limited hydroclimate and yield data, however, hampers drought impact monitoring and assessment at the farmer field-scale. This study demonstrates the potential of field-scale soil moisture simulations to advance high-resolution agricultural yield prediction and drought monitoring at the smallholder farm field-scale. We present a multi-scale modeling approach that combines HydroBlocks, a physically-based hyper-resolution Land Surface Model (LSM), and machine learning. We applied HydroBlocks to simulate root zone soil moisture and soil temperature in Zambia at 3-hourly 30-m resolution. These simulations along with remotely sensed vegetation indices, meteorological conditions, and data describing the physical properties of the landscape (topography, land cover, soil properties) were combined with district-level maize data to train a random forest model (RF) to predict maize yields at the district- and field-scale (250-m) levels. Our model predicted yields with a coefficient of variation (R²) of 0.61, Mean Absolute Error (MAE) of 349 kg ha-1, and mean normalized error of 22 %. We captured maize losses due to the 2015/2016 El Niño drought at similar levels to losses reported by the Food and Agriculture Organization (FAO). Our results revealed that soil moisture is the strongest and most reliable predictor of maize yield, driving its spatial and temporal variability. Consequently, soil moisture was also the most effective indicator of drought impacts in crops when compared with precipitation, soil and air temperatures, and remotely-sensed NDVI-based drought indices. By combining field-scale root zone soil moisture estimates with observed maize yield data, this research demonstrates how field-scale modeling can help bridge the spatial scale discontinuity gap between drought monitoring and agricultural impacts. [ABSTRACT FROM AUTHOR]

In August 2017 Bangladesh faced one of its worst river flooding events in recent history. This paper presents, for the first time, an attribution of this precipitation-induced flooding to anthropogenic climate change from a combined meteorological and hydrological perspective. Experiments were conducted with three observational datasets and two climate models to estimate changes in the extreme 10-day precipitation event frequency over the Brahmaputra basin up to the present and, additionally, an outlook to 2 ∘ C warming since pre-industrial times. The precipitation fields were then used as meteorological input for four different hydrological models to estimate the corresponding changes in river discharge, allowing for comparison between approaches and for the robustness of the attribution results to be assessed. In all three observational precipitation datasets the climate change trends for extreme precipitation similar to that observed in August 2017 are not significant, however in two out of three series, the sign of this insignificant trend is positive. One climate model ensemble shows a significant positive influence of anthropogenic climate change, whereas the other large ensemble model simulates a cancellation between the increase due to greenhouse gases (GHGs) and a decrease due to sulfate aerosols. Considering discharge rather than precipitation, the hydrological models show that attribution of the change in discharge towards higher values is somewhat less uncertain than in precipitation, but the 95 % confidence intervals still encompass no change in risk. Extending the analysis to the future, all models project an increase in probability of extreme events at 2 ∘ C global heating since pre-industrial times, becoming more than 1.7 times more likely for high 10-day precipitation and being more likely by a factor of about 1.5 for discharge. Our best estimate on the trend in flooding events similar to the Brahmaputra event of August 2017 is derived by synthesizing the observational and model results: we find the change in risk to be greater than 1 and of a similar order of magnitude (between 1 and 2) for both the meteorological and hydrological approach. This study shows that, for precipitation-induced flooding events, investigating changes in precipitation is useful, either as an alternative when hydrological models are not available or as an additional measure to confirm qualitative conclusions. Besides this, it highlights the importance of using multiple models in attribution studies, particularly where the climate change signal is not strong relative to natural variability or is confounded by other factors such as aerosols. [ABSTRACT FROM AUTHOR]

In eastern Africa droughts can cause crop failure and lead to food insecurity. With increasing temperatures, there is an a priori assumption that droughts are becoming more severe, however, the link between droughts and climate change is not sufficiently understood. In the current study we focus on agricultural drought and the influence of high temperatures and precipitation deficits on this. Using a combination of models and observational datasets, we studied trends in six regions in eastern Africa in four drought-related annually averaged variables – soil moisture, precipitation, temperature and, as a measure of evaporative demand, potential evapotranspiration (PET). In standardized soil moisture data, we find no discernible trends. Precipitation was found to have a stronger influence on soil moisture variability than temperature or PET, especially in the drier, or water-limited, study regions. The error margins on precipitation-trend estimates are however large and no clear trend is evident. We find significant positive trends in local temperatures. However, the influence of these on soil moisture annual trends appears limited as evaporation is water limited. The trends in PET are predominantly positive, but we do not find strong relations between PET and soil moisture trends. Nevertheless, the PET-trend results can still be of interest for irrigation purposes as it is PET that determines the maximum evaporation rate. We conclude that, until now, the impact of increasing local temperatures on agricultural drought in eastern Africa is limited and recommend that any soil moisture analysis be supplemented by analysis of precipitation deficit. [ABSTRACT FROM AUTHOR]

In August 2017 Bangladesh faced one of its worst river flooding events in recent history. This paper presents for the first time an attribution of this precipitation-induced flooding from a combined meteorological and hydrological perspective. Experiments were conducted with three observational data sets and two climate models to estimate changes in extreme 10-day precipitation event frequency over the Brahmaputra basin. The precipitation fields were then used as meteorological input for four different hydrological models to estimate the corresponding changes in river discharge, allowing for comparison between approaches and for the robustness of the attribution results to be assessed. In all three observational precipitation data sets the climate change trends for extreme precipitation similar to observed in August 2017 are not significant, however in two out of three series, the sign of this insignificant trend is positive. One climate model shows a significant positive influence of anthropogenic climate change, whereas the other simulates a cancellation between the increase due to greenhouse gases and a decrease due to sulphate aerosols. Considering discharge rather than precipitation, the hydrological models show that attribution of the change in discharge towards higher values is somewhat less uncertain than for precipitation, but the 95 % confidence interval still encompasses no change in risk. For the future, all models project an increase in probability of extreme events at 2 °C global heating since pre-industrial times, becoming more than 1.7 times more likely for high 10-day precipitation, and about a factor 1.5 more likely for discharge. Our best estimate on the trend in flooding events similar to the Brahmaputra event of August 2017 is derived by synthesizing the observational and model results: We find the change in risk to be greater than one and of similar order of magnitude (between 1 and 2) for both the meteorological and hydrological approach. This study shows that, for precipitation-induced flooding events, investigating changes in precipitation is useful, either as an alternative when hydrological models are not available, or as an additional measure to confirm qualitative conclusions. Besides, it highlights the importance of using multiple models in attribution studies, particularly where the climate change signal is not strong relative to natural variability or is confounded by other factors such as aerosols. [ABSTRACT FROM AUTHOR]

We present PCR-GLOBWB 2, a global hydrology and water resources model. Compared to previous versions of PCR-GLOBWB, this version fully integrates water use. Sector-specific water demand, groundwater and surface water withdrawal, water consumption, and return flows are dynamically calculated at every time step and interact directly with the simulated hydrology. PCR-GLOBWB 2 has been fully rewritten in Python and PCRaster Python and has a modular structure, allowing easier replacement, maintenance, and development of model components. PCR-GLOBWB 2 has been implemented at 5 arcmin resolution, but a version parameterized at 30 arcmin resolution is also available. Both versions are available as open-source codes on https://github. com/UU-Hydro/PCR-GLOBWB_model (Sutanudjaja et al., 2017a). PCR-GLOBWB 2 has its own routines for groundwater dynamics and surface water routing. These relatively simple routines can alternatively be replaced by dynamically coupling PCR-GLOBWB 2 to a global two-layer groundwater model and 1-D-2-D hydrodynamic models. Here, we describe the main components of the model, compare results of the 30 and 5 arcmin versions, and evaluate their model performance using Global Runoff Data Centre discharge data. Results show that model performance of the 5 arcmin version is notably better than that of the 30 arcmin version. Furthermore, we compare simulated time series of total water storage (TWS) of the 5 arcmin model with those observed with GRACE, showing similar negative trends in areas of prevalent groundwater depletion. Also, we find that simulated total water withdrawal matches reasonably well with reported water withdrawal from AQUASTAT, while water withdrawal by source and sector provide mixed results. [ABSTRACT FROM AUTHOR]

There is growing evidence that climate change will alter water availability in Europe. Here, we investigate how hydrological low flows are affected under different levels of future global warming (i.e. 1.5, 2, and 3K with respect to the pre-industrial period) in rivers with a contributing area of more than 1000 km². The analysis is based on a multi-model ensemble of 45 hydrological simulations based on three representative concentration pathways (RCP2.6, RCP6.0, RCP8.5), five Coupled Model Intercomparison Project Phase 5 (CMIP5) general circulation models (GCMs: GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, NorESM1-M) and three state-of-theart hydrological models (HMs: mHM, Noah-MP, and PCRGLOBWB). High-resolution model results are available at a spatial resolution of 5 km across the pan-European domain at a daily temporal resolution. Low river flow is described as the percentile of daily streamflow that is exceeded 90% of the time. It is determined separately for each GCM/HM combination and warming scenario. The results show that the low-flow change signal amplifies with increasing warming levels. Low flows decrease in the Mediterranean region, while they increase in the Alpine and Northern regions. In the Mediterranean, the level of warming amplifies the signal from -12% under 1.5 K, compared to the baseline period 1971-2000, to -35% under global warming of 3 K, largely due to the projected decreases in annual precipitation. In contrast, the signal is amplified from C22 (1.5 K) to C45% (3 K) in the Alpine region due to changes in snow accumulation. The changes in low flows are significant for regions with relatively large change signals and under higher levels of warming. However, it is not possible to distinguish climate-induced differences in low flows between 1.5 and 2K warming because of (1) the large inter-annual variability which prevents distinguishing statistical estimates of periodaveraged changes for a given GCM/HM combination, and (2) the uncertainty in the multi-model ensemble expressed by the signal-to-noise ratio. The contribution by the GCMs to the uncertainty in the model results is generally higher than the one by the HMs. However, the uncertainty due to HMs cannot be neglected. In the Alpine, Northern, and Mediterranean regions, the uncertainty contribution by the HMs is partly higher than those by the GCMs due to different representations of processes such as snow, soil moisture and evapotranspiration. Based on the analysis results, it is recommended (1) to use multiple HMs in climate impact studies and (2) to embrace uncertainty information on the multimodel ensemble as well as its single members in the adaptation process. [ABSTRACT FROM AUTHOR]

Closing the terrestrial water budget is necessary to provide consistent estimates of budget components for understanding water resources and changes over time. Given the lack of in situ observations of budget components at anything but local scale, merging information from multiple data sources (e.g., in situ observation, satellite remote sensing, land surface model, and reanalysis) through data assimilation techniques that optimize the estimation of fluxes is a promising approach. Conditioned on the current limited data availability, a systematic method is developed to optimally combine multiple available data sources for precipitation (P), evapotranspiration (ET), runoff (R), and the total water storage change (TWSC) at 0.5° spatial resolution globally and to obtain water budget closure (i.e., to enforce P-ET-R-TWSC =0) through a constrained Kalman filter (CKF) data assimilation technique under the assumption that the deviation from the ensemble mean of all data sources for the same budget variable is used as a proxy of the uncertainty in individual water budget variables. The resulting long-term (1984-2010), monthly 0.5° resolution global terrestrial water cycle Climate Data Record (CDR) data set is developed under the auspices of the National Aeronautics and Space Administration (NASA) Earth System Data Records (ESDRs) program. This data set serves to bridge the gap between sparsely gauged regions and the regions with sufficient in situ observations in investigating the temporal and spatial variability in the terrestrial hydrology at multiple scales. The CDR created in this study is validated against in situ measurements like river discharge from the Global Runoff Data Centre (GRDC) and the United States Geological Survey (USGS), and ET from FLUXNET. The data set is shown to be reliable and can serve the scientific community in understanding historical climate variability in water cycle fluxes and stores, benchmarking the current climate, and validating models. [ABSTRACT FROM AUTHOR]

Glaciers are essential hydrological reservoirs, storing and releasing water at various timescales. Short-term variability in glacier melt is one of the causes of streamflow droughts, here defined as deficiencies from the flow regime. Streamflow droughts in glacierised catchments have a wide range of interlinked causing factors related to precipitation and temperature on short and long timescales. Climate change affects glacier storage capacity, with resulting consequences for discharge regimes and streamflow drought. Future projections of streamflow drought in glacierised basins can, however, strongly depend on the modelling strategies and analysis approaches applied. Here, we examine the effect of different approaches, concerning the glacier modelling and the drought threshold, on the characterisation of streamflow droughts in glacierised catchments. Streamflow is simulated with the Hydrologiska Byråns Vattenbalansavdelning (HBV-light) model for two case study catchments, the Nigardsbreen catchment in Norway and theWolverine catchment in Alaska, and two future climate change scenarios (RCP4.5 and RCP8.5). Two types of glacier modelling are applied, a constant and dynamic glacier area conceptualisation. Streamflow droughts are identified with the variable threshold level method and their characteristics are compared between two periods, a historical (1975-2004) and future (2071-2100) period. Two existing threshold approaches to define future droughts are employed: (1) the threshold from the historical period; (2) a transient threshold approach, whereby the threshold adapts every year in the future to the changing regimes. Results show that drought characteristics differ among the combinations of glacier area modelling and thresholds. The historical threshold combined with a dynamic glacier area projects extreme increases in drought severity in the future, caused by the regime shift due to a reduction in glacier area. The historical threshold combined with a constant glacier area results in a drastic decrease of the number of droughts. The drought characteristics between future and historical periods are more similar when the transient threshold is used, for both glacier area conceptualisations. With the transient threshold, factors causing future droughts can be analysed. This study revealed the different effects of methodological choices on future streamflow drought projections and it highlights how the options can be used to analyse different aspects of future droughts: the transient threshold for analysing future drought processes, the historical threshold to assess changes between periods, the constant glacier area to analyse the effect of short-term climate variability on droughts and the dynamic glacier area to model more realistic future discharges under climate change. [ABSTRACT FROM AUTHOR]

Closing the terrestrial water budget is necessary to provide consistent estimates of budget components for understanding water resources and changes over time. Given the lack of in situ observations of budget components at anything but local scale, merging information from multiple data sources (e.g., in situ observation, satellite remote sensing, land surface model, and reanalysis) through data assimilation techniques that optimize the estimation of fluxes is a promising approach. Conditioned on the current limited data availability, a systematic method is developed to optimally combine multiple available data sources for precipitation (P), evapotranspiration (ET), runoff (R), and the total water storage change (TWSC) at 0.5° spatial resolution globally and to obtain water budget closure (i.e., to enforce P -ET-R-TWSC =0) through a constrained Kalman filter (CKF) data assimilation technique under the assumption that the deviation from the ensemble mean of all data sources for the same budget variable is used as a proxy of the uncertainty in individual water budget variables. The resulting long-term (1984-2010), monthly 0.5° resolution global terrestrial water cycle Climate Data Record (CDR) data set is developed under the auspices of the National Aeronautics and Space Administration (NASA) Earth System Data Records (ESDRs) program. This data set serves to bridge the gap between sparsely gauged regions and the regions with sufficient in situ observations in investigating the temporal and spatial variability in the terrestrial hydrology at multiple scales. The CDR created in this study is validated against in situ measurements like river discharge from the Global Runoff Data Centre (GRDC) and the United States Geological Survey (USGS), and ET from FLUXNET. The data set is shown to be reliable and can serve the scientific community in understanding historical climate variability in water cycle fluxes and stores, benchmarking the current climate, and validating models. [ABSTRACT FROM AUTHOR]

We present PCR-GLOBWB 2, a global hydrology and water resources model. Compared to previous versions of PCR-GLOBWB, this version fully integrates water use. Sector-specific water demand, groundwater and surface water withdrawal, water consumption and return flows are dynamically calculated at every time step and interact directly with the simulated hydrology. PCR-GLOBWB 2 has been fully rewritten in Python and PCRaster-Python and has a modular structure, allowing easier replacement, maintenance, and development of model components. PCR-GLOBWB 2 has been implemented at 5 arc-minute resolution, but a version parameterized at 30 arc-minute resolution is also available. Both versions are available as open source codes on https://github.com/UU-Hydro/PCR-GLOBWB_model. PCR-GLOBWB 2 has its own routines for groundwater dynamics and surface water routing. These relatively simple routines can alternatively be replaced by dynamically coupling PCR-GLOBWB 2 to a global two-layer groundwater model and 1D-2D-hydrodynamic models, respectively. Here, we describe the main components of the model, compare results of the 30 arcminute and the 5 arc-minute versions and evaluate their model performance using GRDC discharge data. Results show that model performance of the 5 arc-minute version is notably better than that of the 30 arc-minute version. Furthermore, we compare simulated time series of total water storage (TWS) of the 5 arc-minute model with those observed with GRACE, showing similar negative trends in areas of prevalent groundwater depletion. Also, we find that simulated water withdrawal, by source and sector, matches reasonably well with reported water withdrawal from AQUASTAT. [ABSTRACT FROM AUTHOR]

Drought is an abnormal and prolonged deficit in available water. Possible drought impacts are crop losses, famine, fatalities, power blackouts and degraded ecosystems. These severe socio-economic and environmental impacts show the need to carefully monitor drought conditions using a suitable index. Our objective is to provide an intercomparison of frequently used physical drought indices to show to which degree they are interchangeable for monitoring drought in precipitation, soil moisture, groundwater and streamflow. Physical indices are commonly introduced to predict drought impacts, because appropriate drought impact indices are still missing. Correlations (R) between frequently used indices for different drought types were calculated at the global scale. We have made the index timeseries available to the community for future studies. Precipitation drought indices show low to intermediate correlations (ranging from R = 0.1 to 0.75), soil moisture drought indices show an even lower similarity (R = 0.25). Indices for streamflow drought show the highest correlation (R = 0.5 to 0.95). Additionally, meteorological drought indices do not capture the soil moisture drought correctly (R = 0.0 to 0.6) nor streamflow drought (R = 0.0 to 0.7). These findings have implications for drought monitoring systems: (i) for each drought type, a different index should carefully be identified; (ii) drought indices that are designed to monitor the same drought type show large discrepancies in their anomalies and hence drought detection; (iii) there is no single superior physical drought index that is capable of accurately capturing the diverse set of drought impacts in all parts of the hydrological cycle. [ABSTRACT FROM AUTHOR]

There is growing evidence that climate change will alter water availability in Europe. Here, we investigate how hydrological low flows are affected under different levels of future global warming (i.e., 1.5, 2 and 3 K). The analysis is based on a multi-model ensemble of 45 hydrological simulations based on three RCPs (rcp2p6, rcp6p0, rcp8p5), five CMIP5 GCMs (GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, NorESM1-M) and three state-of-the-art hydrological models (HMs: mHM, Noah-MP, and PCR-GLOBWB). High resolution model results are available at the unprecedented spatial resolution of 5 km across the pan-European domain at daily temporal resolution. Low river flow is described as the percentile of daily streamflow that is exceeded 90 % of the time. It is determined separately for each GCM/HM combinations and the warming scenarios. The results show that the change signal amplifies with increasing warming levels. Low flows decrease in the Mediterranean while they increase in the Alpine and Northern regions. In the Mediterranean, the level of warming amplifies the signal from -12 % under 1.5 K to -35 % under 3 K global warming largely due to the projected decreases in annual precipitation. In contrast, the signal is amplified from +22 % (1.5 K) to +45 % (3 K) in the Alpine region because of the reduced snow melt contribution. The changes in low flows are significant for regions with relatively large change signals and under higher levels of warming. Nevertheless, it is not possible to distinguish climate induced differences in low flows between 1.5 and 2 K warming because of the large variability inherent in the multi-model ensemble. The contribution by the GCMs to the uncertainty in the model results is generally higher than the one by the HMs. However, the uncertainty due to HMs cannot be neglected. In the Alpine and Northern region as well as the Mediterranean, the uncertainty contribution by the HMs is partly higher than those by the GCMs due to different representations of processes such as snow, soil moisture and evapotranspiration. [ABSTRACT FROM AUTHOR]

Over recent decades, the global population has been rapidly increasing and human activities have altered terrestrial water fluxes to an unprecedented extent. The phenomenal growth of the human footprint has significantly modified hydrological processes in various ways (e.g. irrigation, artificial dams, and water diversion) and at various scales (from a watershed to the globe). During the early 1990s, awareness of the potential for increased water scarcity led to the first detailed global water resource assessments. Shortly thereafter, in order to analyse the human perturbation on terrestrial water resources, the first generation of largescale hydrological models (LHMs) was produced. However, at this early stage few models considered the interaction between terrestrial water fluxes and human activities, including water use and reservoir regulation, and even fewer models distinguished water use from surface water and groundwater resources. Since the early 2000s, a growing number of LHMs have incorporated human impacts on the hydrological cycle, yet the representation of human activities in hydrological models remains challenging. In this paper we provide a synthesis of progress in the development and application of human impact modelling in LHMs. We highlight a number of key challenges and discuss possible improvements in order to better represent the human-water interface in hydrological models. [ABSTRACT FROM AUTHOR]

Snow is an important component of water storage in the Himalayas. Previous snowmelt studies in the Himalayas have predominantly relied on remotely sensed snow cover. However, snow cover data provide no direct information on the actual amount of water stored in a snowpack, i.e., the snow water equivalent (SWE). Therefore, in this study remotely sensed snow cover was combined with in situ observations and a modified version of the seNorge snow model to estimate (climate sensitivity of) SWE and snowmelt runoff in the Langtang catchment in Nepal. Snow cover data from Landsat 8 and the MOD10A2 snow cover product were validated with in situ snow cover observations provided by surface temperature and snow depth measurements resulting in classification accuracies of 85.7 and 83.1% respectively. Optimal model parameter values were obtained through data assimilation of MOD10A2 snow maps and snow depth measurements using an ensemble Kalman filter (EnKF). Independent validations of simulated snow depth and snow cover with observations show improvement after data assimilation compared to simulations without data assimilation. The approach of modeling snow depth in a Kalman filter framework allows for data-constrained estimation of snow depth rather than snow cover alone, and this has great potential for future studies in complex terrain, especially in the Himalayas. Climate sensitivity tests with the optimized snow model revealed that snowmelt runoff increases in winter and the early melt season (December to May) and decreases during the late melt season (June to September) as a result of the earlier onset of snowmelt due to increasing temperature. At high elevation a decrease in SWE due to higher air temperature is (partly) compensated by an increase in precipitation, which emphasizes the need for accurate predictions on the changes in the spatial distribution of precipitation along with changes in temperature. [ABSTRACT FROM AUTHOR]

Closing the terrestrial water budget is necessary to providing consistent estimates of budget components for understanding water resources and changes over time. Given the lack of in-situ observations of budget components at anything but local scale, merging information from multiple data sources (e.g. in-situ observation, satellite remote sensing, land surface model and reanalysis) through data assimilation techniques that optimize the estimation of fluxes is a promising approach. In this study, a systematic method is developed to optimally combine multiple available data sources for precipitation (P), evapotranspiration (ET), runoff (R) and the total water storage change (TWSC) at 0.5° spatial resolution globally and to obtain water budget closure (i.e. to enforce P - ET - R - TWSC = 0) through a Constrained Kalman Filter (CKF) data assimilation technique. The resulting long-term (1984-2010), monthly, 0.5° resolution global terrestrial water cycle Climate Data Record (CDR) dataset is developed under the auspices of the National Aeronautics and Space Administration (NASA) Earth System Data Records (ESDRs) program. This dataset serves to bridge the gap between sparsely gauged regions and the regions with sufficient in-situ observations in investigating the temporal and spatial variability in the terrestrial hydrology at multiple scales. The CDR created in this study is validated against in-situ measurements like river discharge from the Global Runoff Data Centre (GRDC) and the United States Geological Survey (USGS) and ET from FLUXNET. The dataset is shown to be reliable and can serve the scientific community in understanding historical climate variability in water cycle fluxes and stores, benchmarking the current climate, and validating models. [ABSTRACT FROM AUTHOR]

Glaciers are essential hydrological reservoirs, storing and releasing water at various time scales. Short-term variability in glacier melt is one of the causes of streamflow droughts, defined as below normal water availabilities. Streamflow droughts in glacierised catchments have a wide range of interlinked causing factors related to precipitation and temperature on short and long time scales. Climate change affects glacier storage capacity, with resulting consequences for discharge regimes and drought. Future projections of streamflow drought in glacierised basins can, however, strongly depend on the modelling strategies and analysis approaches applied. Here, we examine the effect of different approaches, concerning the glacier modelling and the drought threshold, on the characterisation of streamflow droughts in glacierised catchments. Streamflow is simulated with the HBV-light model for two case study catchments, the Nigardsbreen catchment in Norway and the Wolverine catchment in Alaska, and two future climate change scenarios (RCP4.5 and RCP8.5). Two types of glacier modelling are applied, a constant and dynamical glacier area conceptualisation. Streamflow droughts are identified with the variable threshold level method and their characteristics are compared between two periods, a historical (1975-2004) and future (2071-2100) period. Two existing threshold approaches to define future droughts are employed, (1) the threshold from the historical period and (2) a transient threshold approach, whereby the threshold adapts every year in the future to the changing regimes. Results show that drought characteristics differ among the combinations of glacier area modelling and thresholds. The historical threshold combined with a dynamical glacier area projects extreme increases in drought severity in the future, caused by the regime shift due to a reduction in glacier area. The historical threshold combined with a constant glacier area results in a drastic decrease of the number of droughts. The drought characteristics between future and historic periods are more similar when the transient threshold is used, for both glacier dynamics conceptualisations. With the transient threshold causing factors of future droughts, can be analysed. This study revealed the different effects of methodological choices on future streamflow drought projections and it highlights how the options can be used to analyse different aspects of future droughts: the transient threshold for analysing future drought processes, the historical threshold to assess changes between periods, the constant glacier area to analyse the effect of short term climate variability on droughts and the dynamical glacier area to model realistic future discharges under climate change. [ABSTRACT FROM AUTHOR]

Land surface and hydrologic models (LSM/HM) are used at diverse spatial resolutions ranging from 1-10 km in catchment-scale applications to over 50 km in global-scale applications. Application of the same model structure at different spatial scales requires that the model estimates similar fluxes independent of the model resolution and fulfills a flux-matching condition across scales. An analysis of state-of-the-art LSMs and HMs reveals that most do not have consistent and realistic parameter fields for land surface geophysical properties. Multiple experiments with the mHM, Noah-MP, PCR-GLOBWB and WaterGAP models are conducted to demonstrate the pitfalls of poor parameterization practices currently used in most operational models, which are insufficient to satisfy the flux-matching condition. These examples demonstrate that J. Dooge's 1982 statement on the unsolved problem of parameterization in these models remains true. We provide a short review of existing parameter regionalization techniques and discuss a method for obtaining seamless hydrological predictions of water fluxes and states across multiple spatial resolutions. The multiscale parameter regionalization (MPR) technique is a practical and robust method that provides consistent (seamless) parameter and flux fields across scales. A general model protocol is presented to describe how MPR can be applied to a specific model, with an example of this application using the PCR-GLOBWB model. Applying MPR to PCR-GLOBWB substantially improves the flux-matching condition. Estimation of evapotranspiration without MPR at 5 arcmin and 30 arcmin spatial resolutions for the Rhine river basin results in a difference of approximately 29 %. Applying MPR reduce this difference to 9 %. For total soil water, the differences without and with MPR are 25 % and 7 %, respectively. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In the past century, the number and the total storage capacity of large dams and reservoirs have consistently increased. Reservoirs are typically built to satisfy a wide diversity of uses from domestic to industrial purposes. They have (deliberately or not) altered the frequency, the magnitude and the spatial distribution of hydrological extremes, i.e. floods and droughts. Moreover, as their scope can be the mitigation of droughts and/or floods, the coexistence of competing objectives in the management strategy may lead to unintended consequences, as seen for example in the 2011 flooding in Brisbane (Australia), and in the 2018 flooding in Kerala (India). In this work, we investigate the interplay between hydrological extremes and society with a focus on changes in reservoir operation rules. We explore how drought and flood can influence the human impacts on hydrological extremes, and investigate upstream-downstream cascade effects. To this end, the large-scale hydrological model PCR-GLOBWB is implemented to simulate the hydrological response of European catchments. The model also simulate, for the first time, the dynamic interaction between hydrological extremes and changes in reservoirs' operation rules. This work is a first attempt to unravel the feedback mechanisms between society and hydrological extremes throughout the time at large scale. Results provide insights that can be interesting for the interpretation of the role of reservoirs in a wide variety of case studies. [ABSTRACT FROM AUTHOR]

To support the decision-making process at seasonal time scales, hydrological forecasts with a high temporal and spatial resolution are required to provide the level of information needed by water managers. Following the 2018 Northwest European drought, seasonal forecasting became one of the major focus points of European water managers.As part of the EDgE (End-to-end Demonstrator for improved decision making in the water sector in Europe) project commissioned by ECMWF (C3S), we have created a unique dataset of hydrological seasonal forecasts derived from four global climate models (CanCM4, FLOR-B01, ECMF, LFPW) in combination with four global hydrological models (GHMs: PCR-GLOBWB, VIC, mHM, Noah-MP). The forecasts are provided at a daily temporal resolution and 5-km spatial resolution and are bias corrected against E-OBS meteorological observations.Results show that skillful discharge forecasts can be made throughout Europe 3 months in advance, with predictability up to 6 months for Northern Europe due to the impact of snow. Drought forecasts have high skill over Europe due to the large spatial and temporal extent of drought events that is picked up by the seasonal forecasts. The model forecast skill is model dependent, where some GHMs show high skill for high flows, whereas other GHMs have a superior performance with regard to seasonal forecasting of low flows. Finally, the results indicate that the initial model conditions can dominate the forecast skill for multiple months.Even though multi-model seasonal forecast systems like EDgE show a great potential, the communication of the skills of these multi-model systems remains challenging. For drought mitigation, water managers indicated that long lead times are required to deploy low cost, low regret mitigation strategies. Flood mitigation was not helped by the use of seasonal forecast, because the costs of preventive actions are high and the high uncertainty in the forecast made it unsuitable to support decision making. In the end, water managers indicated that they are helped by simplified metrics of uncertainty rather than getting the full uncertainty information for the forecast system.Overall, we conclude that multi-model systems have an advantage as they are better capable of capturing the uncertainty of a forecast, but are challenging when it comes to real-time decision making and operational applications. [ABSTRACT FROM AUTHOR]