We discuss the various performance aspects of parallelizing our transient global-scale groundwater model at 30′′ resolution (30arcsec; °1/41km at the Equator) on large distributed memory parallel clusters. This model, referred to as GLOBGM, is the successor of our 5′ (5arcmin; °1/410km at the Equator) PCR-GLOBWB 2 (PCRaster Global Water Balance model) groundwater model, based on MODFLOW having two model layers. The current version of GLOBGM (v1.0) used in this study also has two model layers, is uncalibrated, and uses available 30′′ PCR-GLOBWB data. Increasing the model resolution from 5′ to 30′′ creates challenges, including increased runtime, memory usage, and data storage that exceed the capacity of a single computer. We show that our parallelization tackles these problems with relatively low parallel hardware requirements to meet the needs of users or modelers who do not have exclusive access to hundreds or thousands of nodes within a supercomputer. For our simulation, we use unstructured grids and a prototype version of MODFLOW 6 that we have parallelized using the message-passing interface. We construct independent unstructured grids with a total of 278 million active cells to cancel all redundant sea and land cells, while satisfying all necessary boundary conditions, and distribute them over three continental-scale groundwater models (168 million - Afro-Eurasia; 77 million - the Americas; 16 million - Australia) and one remaining model for the smaller islands (17 million). Each of the four groundwater models is partitioned into multiple non-overlapping submodels that are tightly coupled within the MODFLOW linear solver, where each submodel is uniquely assigned to one processor core, and associated submodel data are written in parallel during the pre-processing, using data tiles. For balancing the parallel workload in advance, we apply the widely used METIS graph partitioner in two ways: it is straightforwardly applied to all (lateral) model grid cells, and it is applied in an area-based manner to HydroBASINS catchments that are assigned to submodels for pre-sorting to a future coupling with surface water. We consider an experiment for simulating the years 1958-2015 with daily time steps and monthly input, including a 20-year spin-up, on the Dutch national supercomputer Snellius. Given that the serial simulation would require °1/44.5 months of runtime, we set a hypothetical target of a maximum of 16h of simulation runtime. We show that 12 nodes (32 cores per node; 384 cores in total) are sufficient to achieve this target, resulting in a speedup of 138 for the largest Afro-Eurasia model when using 7 nodes (224 cores) in parallel. A limited evaluation of the model output using the United States Geological Survey (USGS) National Water Information System (NWIS) head observations for the contiguous United States was conducted. This showed that increasing the resolution from 5′ to 30′′ results in a significant improvement with GLOBGM for the steady-state simulation when compared to the 5′ PCR-GLOBWB groundwater model. However, results for the transient simulation are quite similar, and there is much room for improvement. Monthly and multi-year total terrestrial water storage anomalies derived from the GLOBGM and PCR-GLOBWB models, however, compared favorably with observations from the GRACE satellite. For the next versions of GLOBGM, further improvements require a more detailed (hydro)geological schematization and better information on the locations, depths, and pumping rates of abstraction wells.

@en

Digital twins of the Earth are digital representations of the Earth system, spanning scales and domains. Their purpose is to monitor, forecast and assess the Earth system and the consequences of human interventions on the Earth system. Providing users with the capability to interact with and interrogate the system, digital twins of the Earth are decision support systems for addressing environmental challenges. By informing humans of their impact on the Earth system, digital twins aspire to promote new pathways moving forward. By answering causal queries through intervention analysis, they can enhance evidence-based policy making. Existing digital twins of the Earth are primarily technological information systems that represent the physical world. However, as the social and physical worlds are intrinsically interconnected, we argue that humans must be accounted for both within and outside digital twins of the Earth: Within twins to represent human impacts and responses that are integral to the Earth system; and outside twins to govern access and development and to guide responsible use of information acquired from twins. Incorporating human interactions in digital twins of the Earth represents a transformative frontier, promising unparalleled insights into Earth system dynamics and empower humans for action.@en

Aquitards are common hydrogeological features in the subsurface. Typically, pumping tests are used to parameterize the hydraulic conductivity of heterogeneous aquitards. However, they do not take spatial variability and uncertainty into account. Alternatively, core-scale measurements of hydraulic conductivity are used in geostatistical upscaling methods, for which their correlation lengths are needed, but this information is extremely difficult to obtain. This study investigates whether a pumping test can be used to obtain the correlation lengths needed for geostatistical upscaling and account for the uncertainty about heterogeneous aquitard conductivity. Random realizations are generated from core-scale data with varying correlation lengths and inserted into a groundwater flow model which simulates the outcome of an actual pumping test. The realizations yielded a better fit to the pumping test data than the traditional pumping test result, assuming homogeneous layers are selected. Ranges of horizontal and vertical correlation lengths that fit the pumping-test well are found. However, considerable uncertainty regarding the correlation lengths remains, which should be considered when parameterizing a regional groundwater flow model.

@en

Aquitards are common hydrogeological features and their hydraulic conductivity is an important property for various groundwater management issues. Predicting their hydraulic conductivity proves challenging, given its dependence on numerous variables. In this study, the dominant factors for predicting aquitard hydraulic conductivity are identified. To this end, a random forest model is trained on a dataset consisting of more than 1000 hydraulic conductivity measurements of core-scale sediment samples from a wide range of stratigraphic units and depths in the Netherlands. The dataset contains textural properties, such as the grain size distribution and porosity, as well as structural data, such as location, sampling depth, stratigraphical unit, lithofacies, organic carbon content, carbonate content and groundwater chloride concentration. Results show that clay fraction, stratigraphic unit, depth, lithofacies and x-coordinate are the most important features for predicting the hydraulic conductivity. Here, x-coordinate is presumably a proxy for distance from marine influence. Using a more detailed grain size distribution or using derived parameters such as the grain size percentiles does not improve the model any further. Our findings indicate that structural properties play a significant role in predicting aquitard conductivity, as they serve as indicators of processes such as compaction and soft-sediment deformation. The model is furthermore an effective method to estimate hydraulic conductivity for sediment samples without conducting costly and time-consuming hydraulic conductivity measurements.

@en

Fresh groundwater reserves, being of vital importance for more than a billion of people living in the coastal zone, are being threatened by saltwater intrusion due to anthropogenic activities and climate change. High resolution three-dimensional (3D), variable-density (VD), groundwater flow and salt transport (FT) numerical models are increasingly being used to support water managers and decision makers in their strategic planning and measures for dealing with the problem of fresh water shortages. However, these computer models typically require long runtimes and large memory usage, making them impractical to use without parallelization. Here, we parallelize SEAWAT, and show that with our parallelization 3D-VD-FT modeling is now feasible for a wide range of hydrogeologists, since a) speedups of more than two orders of magnitude can be obtained as illustrated in this paper, and b) large 3D-VD-FT models are feasible with memory requirements far exceeding single machine memory.

@en

Worldwide, billions of people rely on fresh groundwater reserves for their domestic, agricultural and industrial water use. Extreme droughts and excessive groundwater pumping put pressure on water authorities in maintaining sustainable water usage. High-resolution integrated models are valuable assets in supporting them. The Netherlands Hydrological Instrument (NHI) provides the Dutch water authorities with open source modeling software and data. However, NHI integrated groundwater models often require long run times and large memory usage, therefore strongly limiting their application. As a solution, we present a distributed memory parallelization, focusing on the National Hydrological Model. Depending on the level of integration, we show that significant speedups can be obtained up to two orders of magnitude. As far as we know, this is the first reported integrated groundwater parallelization of an operational hydrological model used for national-scale integrated water management and policy making. The parallel model code and data are freely available.

@en

Whether a coastal area is suitable for beach nourishments and can induce a growth in fresh groundwater resources depends on the appropriateness of the intended site for beach nourishments, and the attainable growth in fresh groundwater resources. In this study we presume that all eroding sandy beaches are suitable for large beach nourishments, and focus on the impact of these nourishments on fresh groundwater in various coastal settings. The growth in fresh groundwater resources – as a consequence of the construction of a beach nourishment – was quantified with 2-D variable-density groundwater models, for a global range in geological parameters and hydrological processes. Our simulation results suggest that large beach nourishments will likely lead to a (temporary) increase of fresh groundwater resources in most settings. However, for a substantial growth in fresh groundwater, the coastal site should receive sufficient groundwater recharge, consist of sediment with a low to medium hydraulic conductivity, and be subject to a limited number of land-surface inundations. Our global analysis shows that 17% of shorelines may consist of erosive sandy beaches, and of these sites 50% have a high potential suitability. This shows a considerable potential worldwide to combine coastal protection with an increase in fresh groundwater resources.@en

This paper is the outcome of a community initiative to identify major unsolved scientific problems in hydrology motivated by a need for stronger harmonisation of research efforts. The procedure involved a public consultation through online media, followed by two workshops through which a large number of potential science questions were collated, prioritised, and synthesised. In spite of the diversity of the participants (230 scientists in total), the process revealed much about community priorities and the state of our science: a preference for continuity in research questions rather than radical departures or redirections from past and current work. Questions remain focused on the process-based understanding of hydrological variability and causality at all space and time scales. Increased attention to environmental change drives a new emphasis on understanding how change propagates across interfaces within the hydrological system and across disciplinary boundaries. In particular, the expansion of the human footprint raises a new set of questions related to human interactions with nature and water cycle feedbacks in the context of complex water management problems. We hope that this reflection and synthesis of the 23 unsolved problems in hydrology will help guide research efforts for some years to come.@en

The Geological Survey of the Netherlands (TNO-GSN) provides public information on the subsurface (http://www.dinoloket.nl). Among this information are (hydro)geological models with hydraulic parametrization. GeoTOP is a voxel model of the upper 50 meters with voxels of 100m x 100m x 0.5m. The basis of the hydraulic parametrisation is a large collection of laboratory measurements of the hydraulic conductivity determined for samples from undisturbed cores, which TNO-GSN collects systematically throughout the Netherlands. The values are upscaled from laboratory to voxel scale for each combination of geological unit and lithoclass. The hydraulic parametrization of GeoTOP for the province of Zeeland (TNO-GSN, 2016) resulted in very high resistances of the confining layer covering the main regional aquifer. The high resistance was predicted for a large area, which made it feasible to validate it. The magnitiude of the resistance is important for the salinization risk for the shallow groundwater and the increase of this risk due to sea level rise. @en

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.

@en

The storage and dynamics of water in the root zone control many important hydrological processes such as saturation excess overland flow, interflow, recharge, capillary rise, soil evaporation and transpiration. Consequently, the effect of the root zone parameterization on the runoff prediction of global hydrological or land surface models (hereafter, GHMs) can be quite large. Many studies have found that climate is the most important factor in explaining rooting depths of natural vegetation, and thus outweighs the importance of soil and vegetation type. Most GHMs, however, use only soil maps and look-up tables for root zone properties to parameterize the root zone. This is problematic, as the observations of the root zone may indirectly have some climate properties, but is still likely to contribute to errors as the data is scarce, often incomplete (i.e. not the entire vertical root profile), at the wrong scale (i.e. individual plant vs. ecosystem value vs. grid cell value), and it is uncertain which part of root zone is hydrologically active. Recent studies developed different optimization methods to a priori determine the root zone storage or rooting depth at the global scale. These studies also suggested that their values could potentially be directly exported to other GHMs. In this study we test which root zone storage and rooting depth would be optimal for the global hydrological model PCR-GLOBWB. We want to answer the following questions: • How much do the “best fit” root zone storages and rooting depths differ from each other with different model settings in PCR-GLOBWB? • Are calibrated root zone storage capacities and rooting depths for PCR-GLOBWB similar to root zone storage capacities and rooting depths obtained by other global studies? The optimal root zone parameters follow from taking the parameters that best reproduce the evaporative flux from a state-of-the-art ensemble of satellite remote sensing based gridded evaporation products. We find that the global spatial patterns of maximum root zone storage capacity are quite independent of other parameters in PCR-GLOBWB. These patterns are similar to the patterns found in recent global studies, but the absolute values are not. We believe that the future for GHMs lies in dynamically determining the root zone storage based on logical – climate based – optimization concepts that are tailor-made for each GHM independently. @en

For a large beach nourishment called the Sand Engine – constructed in 2011 at the Dutch coast – we have examined the impact of coastal forcing (i.e. natural processes that drive coastal hydro- and morphodynamics) and groundwater
recharge on the growth of a fresh groundwater lens between 2011 and 2016. Measurements of the morphological change and the tidal dynamics at the study site were incorporated in a calibrated three-dimensional and variabledensity
groundwater model of the study area. Simulations with this model showed that the detailed incorporation of both the local hydro- and morphodynamics and the actual recharge rate can result in a reliable reconstruction of the growth in fresh groundwater resources. In contrast, the neglect of tidal dynamics, land-surface inundations, and morphological changes in model simulations can result in considerable overestimations of the volume of fresh groundwater. In particular, wave runup and coinciding coastal erosion during storm surges limit the growth in fresh groundwater resources in dynamic coastal environments, and should be considered at potential nourishment sites to delineate the area that is vulnerable to salinization.@en

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.

@en

Recent research in large-scale hydroclimatic variability is surveyed, focusing on five topics: (i) variability in general, (ii) droughts, (iii) floods, (iv) land-atmosphere coupling, and (v) hydroclimatic prediction. Each surveyed topic is supplemented by illustrative examples of recent research, as presented at a 2016 symposium honoring the career of Professor Eric Wood. Taken together, the recent literature and the illustrative examples clearly show that current research into hydroclimatic variability is strong, vibrant, and multifaceted.

@en

Water related catastrophes such as floods are putting more and more people at risk. Moreover this has a large economic impact as well. For example, in 2011 a flood in Bangkok wiped out a large number of harddrive manufacturing plants, leading to a global shortage and increase in price for a two year period. We have a decent grasp on forecasting the weather, especially on the short to medium term (a few days to a week). We have no such grasp for flood forecasting, especially not on the global scale.

@en

The diversity in hydrologic models has historically led to great controversy on the "correct" approach to process-based hydrologic modeling, with debates centered on the adequacy of process parameterizations, data limitations and uncertainty, and computational constraints on model analysis. In this paper, we revisit key modeling challenges on requirements to (1) define suitable model equations, (2) define adequate model parameters, and (3) cope with limitations in computing power. We outline the historical modeling challenges, provide examples of modeling advances that address these challenges, and define outstanding research needs. We illustrate how modeling advances have been made by groups using models of different type and complexity, and we argue for the need to more effectively use our diversity of modeling approaches in order to advance our collective quest for physically realistic hydrologic models.

@en

This study investigates the added value of operational radar with respect to rain gauges in obtaining high-resolution daily rainfall fields as required in distributed hydrological modelling. To this end we combine data from the Netherlands operational national rain gauge network (330 gauges nation wide) with our own experimental network (30 gauges within 225 km2). Based on 74 selected rainfall events (March-October 2004) the spatial variability of daily rainfall is investigated at three spatial extents: small (225 km 2), medium (10,000 km2) and large (82,875 km2). From this analysis it is shown that semivariograms show no clear dependence on season. Predictions of point rainfall are performed for all three extents using three different geostatistical methods: (i) ordinary kriging (OK; rain gauge data only), (ii) kriging with external drift (KED) and (iii) ordinary colocated cokriging (OCCK), the latter two using both rain gauge data and range-corrected daily radar composites - a standard operational radar product from the Royal Netherlands Meteorological Institute (KNMI). Focus here is on automatic prediction. For small extent rain gauge data alone perform better than radar while for larger extents with lower gauge densities radar performs overall better than rain gauge data alone (OK). Methods using both radar and rain gauge data (KED and OCCK) prove to be more accurate than using either rain gauge data alone (OK) or radar, in particular for larger extents. The added value of radar is positively related to the correlation between radar and rain gauge data. Using a pooled semivariogram is almost as good as using event-based semivariograms, which is convenient if the prediction is to be automated. An interesting result is that the pooled semivariograms perform better in terms of estimating the prediction error (kriging variance) especially for the small- and medium extent, where the number of data points to estimate semivariograms is small and event-based semivariograms are rather unstable.@en

This study investigates the added value of operational radar with respect to rain gauges in obtaining high-resolution daily rainfall fields as required in distributed hydrological modeling. To this end data from the Netherlands operational national rain gauge network (330 gauges nationwide) is combined with an experimental network (30 gauges within 225 km2). Based on 74 selected rainfall events (March-October 2004) the spatial variability of daily rainfall is investigated at three spatial extents: small (225 km2), medium (10 000 km2), and large (82 875 km2). From this analysis it is shown that semivariograms show no clear dependence on season. Predictions of point rainfall are performed for all three extents using three different geostatistical methods: (i) ordinary kriging (OK; rain gauge data only), (ii) kriging with external drift (KED), and (iii) ordinary collocated cokriging (OCCK), with the latter two using both rain gauge data and range-corrected daily radar composites - a standard operational radar product from the Royal Netherlands Meteorological Institute (KNMI). The focus here is on automatic prediction. For the small extent, rain gauge data alone perform better than radar, while for larger extents with lower gauge densities, radar performs overall better than rain gauge data alone (OK). Methods using both radar and rain gauge data (KED and OCCK) prove to be more accurate than using either rain gauge data alone (OK) or radar, in particular, for larger extents. The added value of radar is positively related to the correlation between radar and rain gauge data. Using a pooled semivariogram is almost as good as using event-based semivariograms, which is convenient if the prediction is to be automated. An interesting result is that the pooled semivariograms perform better in terms of estimating the prediction error (kriging variance) especially for the small and medium extent, where the number of data points to estimate semivariograms is small and event-based semivariograms are rather unstable.@en