Refine
Document Type
- Article (8)
- Doctoral Thesis (1)
Language
- English (9)
Has Fulltext
- yes (9)
Is part of the Bibliography
- no (9)
Keywords
- Open Science (1)
- global (1)
- groundwater (1)
- hydrology (1)
- models (1)
- reproducibility (1)
- software (1)
Institute
Groundwater is the largest source of accessible freshwater with its dynamics having significantly changed due to human withdrawals, and being projected to continue to as a result of climate change. The pumping of groundwater has led to lowered water tables, decreased base flow, and depletion.
Global hydrological models (GHMs) are used to simulate the global freshwater cycle, assessing impacts of changes in climate and human freshwater use. Currently, groundwater is commonly represented by a bucket-like linear storage component in these models. Bucket models, however, cannot provide information on the location of the groundwater table. Due to this limitation, they can only simulate groundwater discharge to surface water bodies but not recharge from surface water to groundwater and calculate no lateral and vertical groundwater flow whatsoever among grid cells. For instance this may lead to an underestimation of groundwater resources in semiarid areas, where groundwater is often replenished by surface water. In order to overcome these limitations it is necessary to replace the linear groundwater model in GHMs with a hydraulic head gradient-based groundwater flow model
This thesis presents the newly developed global groundwater model G3M and its coupling to the GHM WaterGAP spanning over 70,000 lines of newly developed code. Development and validation of the modeling software are discussed along with numerical challenges. Based on the newly developed software, a global natural equilibrium groundwater model is presented showing better agreements with observations than previous models. Groundwater discharge to rivers is found to be the most dominant flow component globally, compared to flows to other surface water bodies and lateral flows. Furthermore, first global maps of the distribution of gaining and losing surface water bodies are displayed.
For the purpose of determining the uncertainty in model outcomes a sensitivity study is conducted with an innovative approach through applying a global sensitivity analysis for a computationally complex model. First global maps of spatially distributed parameter sensitivities are presented. The results at hand indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, even though parameter sensitivities do vary regionally.
A high resolution model of New Zealand is developed to further understand the involved uncertainties connected to the spatial resolution of the global model. This thesis finds that a new understanding is necessary how these models can be evaluated and that a simple increase in spatial resolution is not improving the model performance when compared to observations.
Alongside the assessment of the natural equilibrium, the concept of a fully coupled transient model as integrated storage component replacing the former model in the hydrological model WaterGAP is discussed. First results reveal that the model shows reasonable response to seasonal variability although it contains persistent head trends leading to global overestimates of water table depth due to an incomplete coupling. Nonetheless, WaterGAP-G3M is already able to show plausible long term storage trends for areas that are known to be affected by groundwater depletion. In comparison with two established regional models in the Central Valley the coupled model shows a highly promising simulation of storage declines.
Global‐scale gradient‐based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady‐state hydraulic heads simulated by two versions of the global groundwater model G3M, at spatial resolutions of 5 arc‐minutes (9 km) and 30 arc‐seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low‐resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global‐scale groundwater models will not improve the simulation of the global freshwater system.
When assessing global water resources with hydrological models, it is essential to know about methodological uncertainties. The values of simulated water balance components may vary due to different spatial and temporal aggregations, reference periods, and applied climate forcings, as well as due to the consideration of human water use, or the lack thereof. We analyzed these variations over the period 1901–2010 by forcing the global hydrological model WaterGAP 2.2 (ISIMIP2a) with five state-of-the-art climate data sets, including a homogenized version of the concatenated WFD/WFDEI data set. Absolute values and temporal variations of global water balance components are strongly affected by the uncertainty in the climate forcing, and no temporal trends of the global water balance components are detected for the four homogeneous climate forcings considered (except for human water abstractions). The calibration of WaterGAP against observed long-term average river discharge Q significantly reduces the impact of climate forcing uncertainty on estimated Q and renewable water resources. For the homogeneous forcings, Q of the calibrated and non-calibrated regions of the globe varies by 1.6 and 18.5 %, respectively, for 1971–2000. On the continental scale, most differences for long-term average precipitation P and Q estimates occur in Africa and, due to snow undercatch of rain gauges, also in the data-rich continents Europe and North America. Variations of Q at the grid-cell scale are large, except in a few grid cells upstream and downstream of calibration stations, with an average variation of 37 and 74 % among the four homogeneous forcings in calibrated and non-calibrated regions, respectively. Considering only the forcings GSWP3 and WFDEI_hom, i.e., excluding the forcing without undercatch correction (PGFv2.1) and the one with a much lower shortwave downward radiation SWD than the others (WFD), Q variations are reduced to 16 and 31 % in calibrated and non-calibrated regions, respectively. These simulation results support the need for extended Q measurements and data sharing for better constraining global water balance assessments. Over the 20th century, the human footprint on natural water resources has become larger. For 11–18% of the global land area, the change of Q between 1941–1970 and 1971–2000 was driven more strongly by change of human water use including dam construction than by change in precipitation, while this was true for only 9–13 % of the land area from 1911–1940 to 1941–1970.
When assessing global water resources with hydrological models, it is essential to know the methodological uncertainties in the water resources estimates. The study presented here quantifies effects of the uncertainty in the spatial and temporal patterns of meteorological variables on water balance components at the global, continental and grid cell scale by forcing the global hydrological model WaterGAP 2.2 (ISI-MIP 2.1) with five state-of-the-art climate forcing input data-sets. While global precipitation over land during 1971–2000 varies between 103 500 and 111 000 km3 yr−1, global river discharge varies between 39 200 and 42 200 km3 yr−1. Temporal trends of global wa- ter balance components are strongly affected by the uncertainty in the climate forcing (except human water abstractions), and there is a need for temporal homogenization of climate forcings (in particular WFD/WFDEI). On about 10–20 % of the global land area, change of river discharge between two consecutive 30 year periods was driven more strongly by changes of human water use including dam construction than by changes in precipitation. This number increases towards the end of the 20th century due to intensified human water use and dam construction. The calibration approach of WaterGAP against observed long-term average river discharge reduces the impact of climate forcing uncertainty on estimated river discharge significantly. Different homgeneous climate forcings lead to a variation of Q of only 1.6 % for the 54 % of global land area that are constrained by discharge observations, while estimated renewable water resources in the remaining uncalibrated regions vary by 18.5 %. Uncertainties are especially high in Southeast Asia where Global Runoff Data Centre (GRDC) data availability is very sparse. By sharing already available discharge data, or installing new streamflow gauging stations in such regions, water balance uncertainties could be reduced which would lead to an improved assessment of the world’s water resources.
The assessment of water balance components using global hydrological models is subject to climate forcing uncertainty as well as to an increasing intensity of human water use within the 20th century. The uncertainty of five state-of-the-art climate forcings and the resulting range of cell runoff that is simulated by the global hydrological model WaterGAP is presented. On the global land surface, about 62 % of precipitation evapotranspires, whereas 38 % discharges into oceans and inland sinks. During 1971–2000, evapotranspiration due to human water use amounted to almost 1 % of precipitation, while this anthropogenic water flow increased by a factor of approximately 5 between 1901 and 2010. Deviation of estimated global discharge from the ensemble mean due to climate forcing uncertainty is approximately 4 %. Precipitation uncertainty is the most important reason for the uncertainty of discharge and evapotranspiration, followed by shortwave downward radiation. At continental levels, deviations of water balance components due to uncertain climate forcing are higher, with the highest discharge deviations occurring for river discharge in Africa (−6 to 11 % from the ensemble mean). Uncertain climate forcings also affect the estimation of irrigation water use and thus the estimated human impact of river discharge. The uncertainty range of global irrigation water consumption amounts to approximately 50 % of the global sum of water consumption in the other water use sector.
In global hydrological models, groundwater (GW) is typically represented by a bucket-like linear groundwater reservoir. Reservoir models, however, (1) can only simulate GW discharge to surface water (SW) bodies but not recharge from SW to GW, (2) provide no information on the location of the GW table, and (3) assume that there is no GW flow among grid cells. This may lead, for example, to an underestimation of groundwater resources in semiarid areas where GW is often replenished by SW or to an underestimation of evapotranspiration where the GW table is close to the land surface. To overcome these limitations, it is necessary to replace the reservoir model in global hydrological models with a hydraulic head gradient-based GW flow model.
We present G3M, a new global gradient-based GW model with a spatial resolution of 5′ (arcminutes), which is to be integrated into the 0.5∘ WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables in-memory coupling to WGHM while keeping overall runtime relatively low, which allows sensitivity analyses, calibration, and data assimilation. This paper presents the G3M concept and model design decisions that are specific to the large grid size required for a global-scale model. Model results under steady-state naturalized conditions, i.e., neglecting GW abstractions, are shown. Simulated hydraulic heads show better agreement to observations around the world compared to the model output of de Graaf et al. (2015). Locations of simulated SW recharge to GW are found, as is expected, in dry and mountainous regions but areal extent of SW recharge may be underestimated. Globally, GW discharge to rivers is by far the dominant flow component such that lateral GW flows only become a large fraction of total diffuse and focused recharge in the case of losing rivers, some mountainous areas, and some areas with very low GW recharge. A strong sensitivity of simulated hydraulic heads to the spatial resolution of the model and the related choice of the water table elevation of surface water bodies was found. We suggest to investigate how global-scale groundwater modeling at 5′ spatial resolution can benefit from more highly resolved land surface elevation data.
In global hydrological models, groundwater storages and flows are generally simulated by linear reservoir models. Recently, the first global gradient-based groundwater models were developed in order to improve the representation of groundwater-surface water interactions, capillary rise, lateral flows and human water use impacts. However, the reliability of model outputs is limited by a lack of data as well as model assumptions required due to the necessarily coarse spatial resolution. The impact of data quality is presented by showing the sensitivity of a groundwater model to changes in the only available global hydraulic conductivity data-set. To better understand the sensitivity of model output to uncertain spatially distributed parameter inputs, we present the first application of a global sensitivity method for a global-scale groundwater model using nearly 2000 steady-state model runs of the global gradient-based groundwater model G3M. By applying the Morris method in a novel domain decomposition approach that identifies global hydrological response units, spatially distributed parameter sensitivities are determined for a computationally expensive model. Results indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, though parameter sensitivities vary regionally. For large areas of the globe, rivers are simulated to be either losing or gaining, depending on the parameter combination, indicating a high uncertainty of simulating the direction of flow between the two compartments. Mountainous and dry regions show a high variance in simulated head due to numerical difficulties of the model, limiting the reliability of computed sensitivities in these regions. This instability is likely caused by the uncertainty in surface water body elevation. We conclude that maps of spatially distributed sensitivities can help to understand complex behaviour of models that incorporate data with varying spatial uncertainties. The findings support the selection of possible calibration parameters and help to anticipate challenges for a transient coupling of the model.
To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of Abstract. To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of capillary rise on evapotranspiration by global hydrological models (GHMs), it is necessary to replace the bucket-like linear GW reservoir model typical for hydrological models with a fully integrated gradient-based GW flow model. Linear reservoir models can only simulate GW discharge to SW bodies, provide no information on the location of the GW table and assume that there is no GW flow among grid cells. A gradient-based GW model simulates not only GW storage but also hydraulic head, which together with information on SW table elevation enables the quantification of water flows from GW to SW and vice versa. In addition, hydraulic heads are the basis for calculating lateral GW flow among grid cells and capillary rise.
G³M is a new global gradient-based GW model with a spatial resolution of 5' that will replace the current linear GW reservoir in the 0.5° WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables inmemory coupling to WGHM while keeping overall runtime relatively low, allowing sensitivity analyses and data assimilation. This paper presents the G³M concept and specific model design decisions together with results under steady-state naturalized conditions, i.e. neglecting GW abstractions. Cell-specific conductances of river beds, which govern GW-SW interaction, were determined based on the 30'' steady-state water table computed by Fan et al. (2013). Together with an appropriate choice for the effective elevation of the SW table within each grid cell, this enables a reasonable simulation of drainage from GW to SW such that, in contrast to the GW model of de Graaf et al. (2015, 2017), no additional drainage based on externally provided values for GW storage above the floodplain is required in G³M. Comparison of simulated hydraulic heads to observations around the world shows better agreement than de Graaf et al. (2015). In addition, G³M output is compared to the output of two established macro-scale models for the Central Valley, California, and the continental United States, respectively. As expected, depth to GW table is highest in mountainous and lowest in flat regions. A first analysis of losing and gaining rivers and lakes/wetlands indicates that GW discharge to rivers is by far the dominant flow, draining diffuse GW recharge, such that lateral flows only become a large fraction of total diffuse and focused recharge in case of losing rivers and some areas with very low GW recharge. G³M does not represent losing rivers in some dry regions. This study presents the first steps towards replacing the linear GW reservoir model in a GHM while improving on recent efforts, demonstrating the feasibility of the approach and the robustness of the newly developed framework.