Search

Abstract. This paper describes the rationale and the protocol of the first component of the third simulation round of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3a, http://www.isimip.org, last access: 2 November 2023) and the associated set of climate-related and direct human forcing data (CRF and DHF, respectively). The observation-based climate-related forcings for the first time include high-resolution observational climate forcings derived by orographic downscaling, monthly to hourly coastal water levels, and wind fields associated with historical tropical cyclones. The DHFs include land use patterns, population densities, information about water and agricultural management, and fishing intensities. The ISIMIP3a impact model simulations driven by these observation-based climate-related and direct human forcings are designed to test to what degree the impact models can explain observed changes in natural and human systems. In a second set of ISIMIP3a experiments the participating impact models are forced by the same DHFs but a counterfactual set of atmospheric forcings and coastal water levels where observed trends have been removed. These experiments are designed to allow for the attribution of observed changes in natural, human, and managed systems to climate change, rising CH4 and CO2 concentrations, and sea level rise according to the definition of the Working Group II contribution to the IPCC AR6.

Water - Global Assessment and Prognosis (WaterGAP) is a modelling approach for quantifying water resources and water use for all land areas of the Earth that has served science and society since 1996. In this paper, the refinements, new algorithms and new data of the most recent model version v2.2e are described, together with a thorough evaluation of simulated water use, streamflow and total water storage anomaly against observation data. WaterGAP v2.2e improves the handling of inland sinks and now excludes not only large but also small man-made reservoirs when simulating naturalized conditions. The reservoir and non-irrigation water use data were updated. In addition, the model was calibrated against an updated and extended dataset of streamflow observations at 1509 gauging stations. The model can now be started using pre-scribed water storages and other conditions, which facilitates data assimilation as well as near real-time monitoring and forecast simulations. For specific applications, the model can consider the output of a glacier model, approximate the effect of rising CO2 concentrations on evapotranspiration or calculate the water temperature in rivers. In the paper, the publicly available standard model output is described and caveats of the model version are provided alongside the description of the model setup in the ISIMIP3 framework. [ABSTRACT FROM AUTHOR]

Although river discharge simulations from global hydrological models have undergone extensive validation, there has been less validation of reservoir operations, primarily because of limited observational data. However, recent advancements in satellite remote sensing technology have facilitated the collection of valuable data regarding water surface area and elevation, thereby providing the ability to validate reservoir storage. In this study, we sought to establish a methodology for validation and intercomparison of reservoir storage within global hydrological model simulations using satellite-derived data. Accordingly, we chose two satellite-derived reservoir operation products, DAHITI and GRSAD, to create monthly time series storage data for seven reservoirs in the contiguous United States (CONUS), with access to long-term ground truth data (the total catchment area accounts for about 9% of CONUS). We assessed two global hydrological models that participated in the Inter Sectoral Model Intercomparison Project (ISIMIP) Phase 3 project, H08 and WaterGAP2, with three distinct forcing datasets: GSWP3-W5E5 (GW), CR20v3-W5E5 (CW), and CR20v3-ERA5 (CE). The results indicated that WaterGAP2 generally outperforms H08; the CW forcing dataset demonstrated superior results compared with GW and CE; the DAHITI showed better consistency with ground observations than GRSAD if temporal coverage is sufficient. Overall, our study emphasizes the potential uses of satellite remote sensing data in reservoir operations validation and underscores the importance of normalization and decomposition techniques for improved validation efficacy. The results highlight the relative performances of different hydrological models and forcing datasets, yielding insights concerning future advancements in reservoir simulation and operational 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 riverine ecosystems. However, at the global scale, the assessment of EFs is associated with significant uncertainty. Here, we present a novel method to determine EFs by Environmental Flow Envelopes (EFE), which is an envelope of variability bounded by discharge limits within which riverine ecosystems are not seriously compromised. The EFE is defined globally in approximately 4,400 sub-basins at monthly time resolution, considering also the methodological uncertainties related with global EF studies. In addition to a lower bound of discharge, the EFE introduces an upper bound of discharge, identifying areas where streamflow has increased substantially. Further, instead of only showing whether EFs are violated, as commonly done, we quantify, for the first time, the frequency, severity, and trends of EFE violations, which can be considered as potential threats to riverine ecosystems. We use pre-industrial (1801-1860) quasi-natural discharge and a suite of hydrological EFR methods and global hydrological models to estimate EFE, applying data from the ISIMIP 2b ensemble. We then compare the EFEs to recent past (1976-2005) discharge to assess the violations of the EFE. We found that the EFE violations most commonly manifest themselves by insufficient streamflow during the low flow season, with less violations during intermediate flow season, and only few violations during high flow season. These violations are widespread: discharge in half of the sub-basins of the world has violated the EFE during more than 5 % of the months between 1976 and 2005. The trends in EFE violations have mainly been increasing during the past decades and will likely remain problematic with projected increases in anthropogenic water use and hydroclimatic changes. Indications of excessive streamflow through EFE upper bound violations are relatively scarce and spatially distributed, although signs of increasing trends can be identified and potentially attributed to climate change. While the EFE provides a quick and globally robust way of determining environmental flow allocations at the sub-basin scale, local fine-tuning is necessary for practical applications and further research on the coupling between quantitative discharge and riverine ecosystem responses is required. [ABSTRACT FROM AUTHOR]

Studies of the global sea-level budget (SLB) and the global ocean-mass budget (OMB) are essential to assess the reliability of our knowledge of sea-level change and its contributions. Here we present datasets for times series of the SLB and OMB elements developed in the framework of ESA's Climate Change Initiative. We use these datasets to assess the SLB and the OMB simultaneously, utilising a consistent framework of uncertainty characterisation. The time series, given at monthly sampling, include global mean sea-level (GMSL) anomalies from satellite altimetry; the global mean steric component from Argo drifter data with incorporation of sea surface temperature data; the ocean mass component from Gravity Recovery and Climate Experiment (GRACE) satellite gravimetry; the contribution from global glacier mass changes assessed by a global glacier model; the contribution from Greenland Ice Sheet and Antarctic Ice Sheet mass changes, assessed from satellite radar altimetry and from GRACE; and the contribution from land water storage anomalies assessed by the WaterGAP global hydrological model. Over the period Jan 1993–Dec 2016 (P1, covered by the satellite altimetry records), the mean rate (linear trend) of GMSL is 3.05 ± 0.24 mm yr−1. The steric component is 1.15 ± 0.12 mm yr−1 (38 % of the GMSL trend) and the mass component is 1.75 ± 0.12 mm yr−1 (57 %). The mass component includes 0.64 ± 0.03 mm yr−1 (21 % of the GMSL trend) from glaciers outside Greenland and Antarctica, 0.60 ± 0.04 mm yr−1 (20 %) from Greenland, 0.19 ± 0.04 mm yr−1 (6 %) from Antarctica, and 0.32 ± 0.10 mm yr−1 (10 %) from changes of land water storage. In the period Jan 2003–Aug 2016 (P2, covered by GRACE and the Argo drifter system), GMSL rise is higher than in P1 at 3.64 ± 0.26 mm yr−1. This is due to an increase of the mass contributions (now about 2.22 ± 0.15 mm yr−1, 61 % of the GMSL trend), with the largest increase contributed from Greenland. The SLB of linear trends is closed for P1 and P2, that is, the GMSL trend agrees with the sum of the steric and mass components within their combined uncertainties. The OMB budget, which can be evaluated only for P2, is also closed, that is, the GRACE-based ocean-mass trend agrees with the sum of assessed mass contributions within uncertainties. Combined uncertainties (1-sigma) of the elements involved in the budgets are between 0.26 and 0.40 mm yr−1, about 10 % of GMSL rise. Interannual variations that overlie the long-term trends are coherently represented by the elements of the SLB and the OMB. Even at the level of monthly anomalies the budgets are closed within uncertainties, while also indicating possible origins of remaining misclosures. [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]

WaterGAP is a global hydrological model that quantifies human use of groundwater and surface water as well as water flows and water storage and thus water resources on all land areas of the Earth. Since 1996, it has served to assess water resources and water stress both historically and in the future, in particular under climate change. It has improved our understanding of continental water storage variations, with a focus on overexploitation and depletion of water resources. In this paper, we describe the most recent model version WaterGAP 2.2d, including the water use models, the linking model that computes net abstractions from groundwater and surface water and the WaterGAP Global Hydrology Model WGHM. Standard model output variables that are freely available at a data repository are explained. In addition, the most requested model outputs, total water storage anomalies, streamflow and water use, are evaluated against observation data. Finally, we show examples of assessments of the global freshwater system that can be done with WaterGAP2.2d model output. [ABSTRACT FROM AUTHOR]

Billions of people rely on groundwater as an accessible source for drinking water and irrigation, especially in times of drought. Its importance will likely increase with a changing climate. It is still unclear, however, how climate change will impact groundwater systems globally and thus the availability of this vital resource. This study investigates uncertainties in groundwater recharge projections using a multi-model ensemble of eight global hydrological models (GHMs) that are driven by the bias-adjusted output of four global circulation models (GCMs). Preindustrial and current groundwater recharge values are compared with recharge for different global warming (GW) levels as a result of three representative concentration pathways (RCPs). Results suggest that the uncertainty range is extensive, and projections with confidence can only be made for specific regions of the world. In some regions, reversals of groundwater recharge trends can be observed with global warming. On average, a consistent median increase of groundwater recharge in northern Europe of 19 % and a decrease of 10 % in the Amazon at 3 °C GW compared to preindustrial levels are simulated. In the Mediterranean, a 2 °C GW leads to a reduction of GWR of 38 %. Because most GHMs do not include CO2 driven vegetation processes, we investigate how, including the effect of evolving CO2 concentrations into the calculation of future groundwater recharge impacts the results. In some regions, the inclusion of these processes leads to differences in groundwater recharge changes of up to 100 mm year-1. Overall, models that include CO2 driven vegetation processes simulate less severe decreases of groundwater recharge and in some regions even increases instead of decreases. In regions where GCMs predict decreases in precipitation, and groundwater availability is most important, the model agreement among GHMs with dynamic vegetation is lowest in contrast to GHMs without, which show a high agreement. [ABSTRACT FROM AUTHOR]

Continental water mass change affects ocean mass change (OMC). Assessing the net contribution, however, remains a challenge. We present an integrated version of the WaterGAP global hydrological model that is able to simulate total continental water storage anomalies (TWSA) over the global continental area (except Greenland and Antarctica) consistently by integrating the output from the global glacier model of Marzeion et al. (2012) as an input to WaterGAP. Monthly time series of global mean TWSA obtained with an ensemble of four variants of the integrated model, corresponding to different precipitation input and irrigation water use assumptions, were validated against an ensemble of four TWSA solutions based on GRACE satellite gravimetry over January 2003 to August 2016. The overall fit to GRACE, measured by the Nash-Sutcliffe efficiency (NSE) coefficient, was found to be 0.87. By decomposing the original TWSA signal into its seasonal, linear trend and inter-annual components, we find that the seasonal amplitude and phase are very well reproduced (NSE=0.88), the linear trend is overestimated by 30-50% (NSE=0.65) and inter-annual variability is captured to a certain extent (NSE=0.57) by the integrated model. During the period 1948-2016, we find that continents lost 34-41mm of sea level equivalent (SLE) to the oceans, with global glacier mass loss accounting for 81% of the cumulated mass loss and glacier-free land water storage anomalies (LWSA) accounting for the remaining 19%. Over 1948-2016, the mass gain on land from impoundment of water in man-made reservoirs, equivalent to 8mm SLE, was offset by the mass loss from water abstractions, amounting to 15-21mm SLE and reflecting a cumulated groundwater depletion of 13-19mm SLE. Climate-driven LWSA are highly sensitive to precipitation input and correlate with El Niño Southern Oscillation multi-year modulations. Significant uncertainty remains in trends of modelled LWSA, which are highly sensitive to simulation of irrigation water use and man-made reservoirs. [ABSTRACT FROM AUTHOR]

To improve the understanding of trends in extreme flows related to flood events at the global scale, historical and future changes of annual maximum streamflow are investigated, using a comprehensive streamflow archive and six global hydrological models. The models' capacity to characterise trends in annual maximum streamflow is evaluated across 3,666 river gauge locations over the period from 1971 to 2005, focusing on four aspects of trends over continental and global scale: (i) mean, (ii) standard deviation, (iii) percentage of locations showing significant trends and (iv) spatial pattern. Compared to observed trends, simulated trends driven by observed climate forcing generally have a higher mean, lower spread, and a similar percentage of locations showing significant trends. Models show a moderate capacity to simulate spatial patterns of historical trends, with approximately only 12-25 % of the spatial variance of observed trends across all gauge stations accounted for by the simulations. Interestingly, there are significant differences between trends simulated by GHMs forced with historical climate and forced by bias corrected climate model output during the historical period, suggesting the important role of the stochastic natural (decadal, inter-annual) climate variability. Significant differences were found in simulated flood trend results when averaged only at gauged locations compared to when averaging across all simulated grid cells, highlighting the potential for bias toward well-observed regions in the state-of-understanding of changes in floods. Future climate projections (simulated under RCP2.6 and RCP6.0 greenhouse gas concentration scenario) suggest a potentially high level of change in individual regions, with up to 35 % of cells showing a statistically significant trend (increase or decrease) and greater changes indicated for the higher concentration pathway. Importantly, the observed streamflow database under-samples the percentage of high-risk locations under RCP6.0 greenhouse gas concentration scenario by more than an order of magnitude (0.9 % compared to 11.7 %). This finding indicates a highly uncertain future for both flood-prone communities and decision makers in the context of climate change. [ABSTRACT FROM AUTHOR]

Human activities, as well as climate change, have had increasing impacts on natural hydrological systems, particularly streamflow. However, quantitative assessments of these impacts are lacking on large scales. In this study, we use the simulations from six global hydrological models driven by three meteorological forcings to investigate direct human impact (DHI) and climate change impact on streamflow in China. Results show that, in the sub-periods of 1971-1990 and 1991-2010, one-fifth to one-third of mean annual streamflow (MAF) reduced due to DHI in northern basins and much smaller (<4%) MAF reduced in southern basins. From 1971-1990 to 1991-2010, total MAF changes range from -13% to 10% across basins, wherein the relative contributions of DHI change and climate change show distinct spatial patterns. DHI change caused decreases in MAF in 70% of river segments, but climate change dominated the total MAF changes in 88% of river segments of China. In most northern basins, climate change results in changes of -9% to 18% of MAF, while DHI change results in decreases of 2% to 8% in MAF. In contrast with the impacts of climate change that may increase or decrease streamflow, DHI change almost always contributes to decreases in MAF over time, wherein water withdrawals are supposed to be the major impact on streamflow. This quantitative assessment can be a reference for attribution of streamflow changes at large scales despite uncertainty remains. We highlight the significant DHI in northern basins and the necessity to modulate DHI through improved water management towards a better adaptation to future climate change. [ABSTRACT FROM AUTHOR]

This paper describes the IMPACT Global Hydrological Model (IGHM), a component of the International Model for Policy Analysis of Agricultural Commodities and Trade (IMPACT) integrated modeling system. IMPACT has been developed in the early 1990s to identify and analyze long-term challenges and opportunities for food, agriculture, and natural resources at global and regional scales and builds on a series of previous food demand and supply projections models developed at the International Food Policy Research Institute since the early 1980s. The IGHM has been developed to assess water availability and variability as drivers of water use and irrigated crop production in IMPACT. It adopts a saturation runoff generation scheme and uses a linear groundwater reservoir to simulate base flow in 0.5º latitude by 0.5º longitude grid cells over the global land surface excluding Antarctica. The IGHM has four cell-specific calibration parameters, which are determined through maximizing the Kling-Gupta efficiency (KGE) with a genetic algorithm at the grid cell level, using gridded natural runoff series generated by the WaterGAP Global Hydrological Model (WGHM). During the calibration and validation periods, globally, the majority of grid cells attain KGE values greater than 0.50. As a meta-model of the more computationally expensive WGHM, IGHM transfers the climate-hydrology dynamics provided by WGHM into the integrated IMPACT model at a lower computational cost and enables coupling hydrology and other related processes considered in IMPACT which are important for analyzing long-term water and food security under a range of environmental and socioeconomic changes. [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]

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 many hydrological models, the amount of evapotranspired water is calculated using the potential evapotranspiration (PET) approach. The main driver of several PET approaches is net radiation, whose downward components are usually obtained from meteorological input data, whereas the upward components are calculated by the model itself. Thus, uncertainties can be large due to both the input data and model assumptions. In this study, we compare the radiation components of the WaterGAP Global Hydrology Model, driven by two meteorological input datasets and two radiation setups from ERA-Interim reanalysis. We assess the performance with respect to monthly observations provided by the Baseline Surface Radiation Network (BSRN) and the Global Energy Balance Archive (GEBA). The assessment is done for the global land area and specifically for energy/water limited regions. The results indicate that there is no optimal radiation input throughout the model variants, but standard meteorological input datasets perform better than those directly obtained by ERA-Interim reanalysis for the key variable net radiation. The low number of observations for some radiation components, as well as the scale mismatch between station observations and 0.5° × 0.5° grid cell size, limits the assessment. [ABSTRACT FROM AUTHOR]