The paper analyzes the national DK-model hydrological information and prediction (HIP) system and HIP portal viewed as a ‘digital twin’ and how the introduction of real-time dynamic updating of the DK-model HIP simulations can make room for plug-in submodels with real-time boundary conditions made available from an HIP portal. The possible feedback to a national real-time risk knowledge base during extreme events (flooding and drought) is also discussed. Under climate change conditions, Denmark is likely to experience more rain in winter, more evapotranspiration in summer, intensified cloudbursts, drought, and sea level rise. These challenges were addressed as part of the Joint Governmental Digitalization Strategy 2016–2020 for better use and sharing of public data about the terrain, water, and climate to support climate adaptation, water management, and disaster risk reduction. This initiative included the development of a new web-based data portal (HIP portal) developed by the Danish Agency for Data Supply and Infrastructure (SDFI). GEUS delivered 5 terabytes of hydrological model data to the portal, with robust calibration methods and hybrid machine learning (ML) being key parts of the deliverables. This paper discusses the challenges and potentials of further developing the HIP digital twin with ‘plug-in digital twins’ for local river basins, including feedback to the national level.

Accurate streamflow estimation is essential for effective water resource management and adapting to extreme events in the face of changing climate conditions. Hydrological models have been the conventional approach for streamflow interpolation and extrapolation in time and space for the past few decades. However, their large-scale applications have encountered challenges, including issues related to efficiency, complex parameterization, and constrained performance. Deep learning methods, such as long short-term memory (LSTM) networks, have emerged as a promising and efficient approach for large-scale streamflow estimation. In this study, we have conducted a series of experiments to identify optimal hybrid modeling schemes to consolidate physically based models with LSTM aimed at enhancing streamflow estimation in Denmark. The results show that the hybrid modeling schemes outperformed the Danish National Water Resources Model (DKM) in both gauged and ungauged basins. While the standalone LSTM rainfall–runoff model outperformed DKM in many basins, it faced challenges when predicting the streamflow in groundwater-dependent catchments. A serial hybrid modeling scheme (LSTM-q), which used DKM outputs and climate forcings as dynamic inputs for LSTM training, demonstrated higher performance. LSTM-q improved the mean Nash–Sutcliffe efficiency (NSE) by 0.22 in gauged basins and 0.12 in ungauged basins compared to DKM. Similar accuracy improvements were achieved with alternative hybrid schemes, i.e., by predicting the residuals between DKM-simulated streamflow and observations using LSTM. Moreover, the developed hybrid models enhanced the accuracy of extreme events, which encourages the integration of hybrid models within an operational forecasting framework. This study highlights the advantages of synergizing existing physically based hydrological models (PBMs) with LSTM models, and the proposed hybrid schemes hold the potential to achieve high-quality large-scale streamflow estimations.

Monitoring the variation of rivers and lakes is of great importance. Satellite radar altimetry is a promising technology to do this on a regional to global scale. Satellite radar altimetry data has been used successfully to observe water levels in lakes and (large) rivers, and has also been combined with hydrologic/hydrodynamic models. Except CryoSat-2, all radar altimetry missions have been operated in conventional low resolution mode with a short repeat orbit (35 days or less). CryoSat-2, carrying a Synthetic Aperture Radar (SAR) altimeter, has a 369-day repeat and a drifting ground track pattern and provides new opportunities for hydrologic research. The narrow inter-track distance (7.5 km at the equator) makes it possible to monitor many lakes and rivers and SAR mode provides a finer along-track resolution, higher return power and speckle reduction through multi-looks. However, CryoSat-2 challenges conventional ways of dealing with satellite inland water altimetry data because virtual station time series cannot be directly derived for rivers. We review the CryoSat-2 mission characteristics, data products, and its use and perspectives for inland water applications. We discuss all the important steps in the workflow for hydrologic analysis with CryoSat-2, and conclude with a discussion of promising future research directions.

Availability of in situ river monitoring data, especially of data shared across boundaries, is decreasing, despite growing challenges for water resource management across the entire globe. This is especially valid for the case study of this work, the Brahmaputra Basin in South Asia. Commonly, satellite altimeters are used in various ways to provide information about such river basins. Most missions provide virtual station time series of water levels at locations where their repeat orbits cross rivers. CryoSat-2 is equipped with a new type of altimeter, providing estimates of the actual ground location seen in the reflected signal. It also uses a drifting orbit, challenging conventional ways of processing altimetry data to river water levels and their incorporation in hydrologic–hydrodynamic models. However, CryoSat-2 altimetry data provides an unprecedentedly high spatial resolution. This paper suggests a procedure to (i) filter CryoSat-2 observations over rivers to extract water-level profiles along the river, and (ii) use this information in combination with a hydrologic–hydrodynamic model to fit the simulated water levels with an accuracy that cannot be reached using information from globally available digital elevation models (DEMs) such as from the Shuttle Radar Topography Mission (SRTM) only. The filtering was done based on dynamic river masks extracted from Landsat imagery, providing spatial and temporal resolutions high enough to map the braided river channels and their dynamic morphology. This allowed extraction of river water levels over previously unmonitored narrow stretches of the river. In the Assam Valley section of the Brahmaputra River, CryoSat-2 data and Envisat virtual station data were combined to calibrate cross sections in a 1-D hydrodynamic model of the river. The hydrologic–hydrodynamic model setup and calibration are almost exclusively based on openly available remote sensing data and other global data sources, ensuring transferability of the developed methods. They provide an opportunity to achieve forecasts of both discharge and water levels in a poorly gauged river system.

There is an urgent demand for assessments of climate change impacts on the hydrological cycle at high spatial resolutions. In particular, the impacts on shallow groundwater levels, which can lead to both flooding and drought, have major implications for agriculture, adaptation, and urban planning. Predicting such hydrological impacts is typically performed using physically based hydrological models (HMs). However, such models are computationally expensive, especially at high spatial resolutions. This study is based on the Danish national groundwater model, set up as a distributed, integrated surface–subsurface model at a 500 m horizontal resolution. Recently, a version at a higher resolution of 100 m was created, amongst others, to better represent the uppermost groundwater table and to meet end-user demands for water management and climate adaptation. However, the increase in resolution of the hydrological model also increases computational bottleneck. To evaluate climate change impacts, a large ensemble of climate models was run with the 500 m hydrological model, while performing the same ensemble run with the 100 m resolution nationwide model was deemed infeasible. The desired outputs at the 100 m resolution were produced by developing a novel, hybrid downscaling method based on machine learning (ML). Hydrological models for five subcatchments, covering around 9 % of Denmark and selected to represent a range of hydrogeological settings, were run at 100 m resolutions with forcings from a reduced ensemble of climate models. Random forest (RF) algorithms were established using the simulated climate change impacts (future – present) on water table depth at 100 m resolution from those submodels as training data. The trained downscaling algorithms were then applied to create nationwide maps of climate-change-induced impacts on the shallow groundwater table at 100 m resolutions. These downscaled maps were successfully validated against results from a validation submodel at a 100 m resolution excluded from training the algorithms, and compared to the impact signals from the 500 m HM across Denmark. The suggested downscaling algorithm also opens for the spatial downscaling of other model outputs. It has the potential for further applications where, for example, computational limitations inhibit running distributed HMs at fine resolutions.