Beyond river discharge gauging: hydrologic predictions using remote sensing alone

As a result of no direct measurement data, accurate stream flow prediction in ungauged basins has remained a challenge in water resource management and flood forecasting. To meet this challenge, advanced hydrological models can be integrated with remote sensing and machine learning techniques. Hydrological models are used to understand the dynamics of streamflow patterns while remotely sensed information provides real-time information on land surface conditions as well as weather parameters thus enhancing the accuracy of these models. Using classical approaches, machine learning methods enable the analysis of complex datasets to detect interrelations or patterns that may not be obvious. Therefore, an integrated approach is crucial for improving flow predictions in ungauged basins that ultimately support effective water management and flood mitigation strategies. This can become even more complicated with climate change as shown by Morocco's high sensitivity to irregular rainfall patterns brought about by climate variations. One such place is the Ouzoud watershed, which stretches 315.6 km 2 across Central High Atlas Mountains and is characterized by frequent flooding leading to socio-economic impacts. This study aims at assessing streamflow prediction methods particularly the Quantile-Quantile Plot (QPPQ) method endorsed by USGS for its applicability in predicting streamflow at ungauged locations. The study area includes both gauged and ungauged watersheds such as the Bernat River, Lakhdar River, Ghassaf River, and Tassaouet River, with Ouzoud River being the only ungauged basin. The geological setting is characterized by carbonate rocks (limestone) throughout the catchment. Mediterranean climate mountainous conditions prevail in this region that experience significant rainfall variations. Other watersheds within the survey site are those of Tassaout and Bernat which have different geological and climatic influences on streamflow patterns. The QPPQ method entails estimating the Flow Duration Curve (FDC) for the catchment of interest, identifying matching donor sites in terms of hydrological characteristics, transferring non-excessive probabilities from these donor sites and mapping these probabilities to streamflow estimates. The performance of the QPPQ method is assessed based on several statistical criteria including Correlation Coefficient (CC), Coefficient of Determination (R 2 ), Mean Absolute Error (MAE), Root Mean Squared Error (RMSE), Nash-Sutcliffe Efficiency (NSE), and Kling-Gupta Efficiency (KGE). For instance, the CC value of 0.55 at Bernat River Sgatt station means that there exists a moderate linear relationship between predicted and observed streamflows while R 2 value at 0.36 indicates that about 36 % of the explained variation in observed flow is covered by predicted ones. Meanwhile, MAE=22.50 mirrors an average absolute prediction error while RMSE equals 55.01 in this case refers to the square root of mean squared prediction errors. Furthermore, NSE gives information concerning how far off the model's results are from the actual data with NSE=0.36 denoting a moderate fit and KGE which equals 0.34 showing good agreement in overall terms. Together, these measures expose the strengths and weaknesses of QPPQ as a streamflow prediction technique. Although the method gives some useful indications, it should be improved further by integrating it with more data sources and processes to increase the accuracy of its predictions and make them more reliable for ungauged basins. Graphical Summary of the Study Using the QPPQ Method Applied to a Gauged Basin (Bernat Basin) and an Ungauged Basin (Ouzoud Basin) • Improved QPPQ method enhances streamflow prediction in ungauged basins. • Seasonal analysis shows best QPPQ performance in winter, weakest in autumn. • Hybrid QPPQ with LSTM and RF improves prediction of extreme streamflows. • QPPQ shows strong results in Ouzoud Basin; basin-specific calibration is crucial. • Integrating remote sensing and machine learning supports water management.

Streamflow prediction in ungauged basins (PUB) is challenging, and Long Short‐Term Memory (LSTM) is widely used to for such predictions, owing to its excellent migration performance. Traditional LSTM forced by meteorological data and catchment attribute data barely highlight the optimum data integration strategy for LSTM and its migration from data‐rich basins to ungauged ones. In this study, we experimented with 1,897 global catchments and found that LSTM‐corrected Global Hydrological Models (GHMs) outperformed uncorrected GHMs, improving the median Nash‐Sutcliff efficiency (NSE) from 0.03 to 0.66. Notably, there was a large gap between traditional LSTM modeling in ungauged basins and autoregressive modeling in data‐rich basins, and GHM‐forced LSTM were an effective way to close this gap in ungauged basins. The spatial heterogeneity of the performance of GHM‐forced LSTM was mainly influenced by three metrics (dryness, the leaf area index and latitude), which described the hydrological similarity among catchments. Weaker hydrological similarity among continental catchments results in larger variability in GHM‐forced LSTM, with the best performance in Siberia (NSE, 0.54) and the worst in North America (NSE, 0.10). However, the migration performance of GHM‐forced LSTM was significantly improved (NSE, 0.63) in ungauged basins when hydrological similarity was considered. This study stressed the advantages of GHM‐forced LSTM and due significance should be attached to hydrological similarities among catchments to improve hydrological prediction in ungauged catchments.

Streamflow prediction is a vital public service that helps to establish flash-flood early warning systems or assess the impact of projected climate change on water management. However, the availability of streamflow observations limits the utilization of the state-of-the-art streamflow prediction techniques to the basins where hydrometric gauging stations exist. Since the most river basins in the world are ungauged, the development of the specialized techniques for the reliable streamflow prediction in ungauged basins (PUB) is of crucial importance. In recent years, the emerging field of deep learning provides a myriad of new models that can breathe new life into the stagnating PUB methods. In the presented study, we benchmark the streamflow prediction efficiency of Long Short-Term Memory (LSTM) networks against the standard technique of GR4J hydrological model parameters regionalization (HMREG) at 200 basins in Northwest Russia. Results show that the LSTM-based regional hydrological model significantly outperforms the HMREG scheme in terms of median Nash-Sutcliffe efficiency (NSE), which is 0.73 and 0.61 for LSTM and HMREG, respectively. Moreover, LSTM demonstrates the comparable median NSE with that for basin-scale calibration of GR4J (0.75). Therefore, this study underlines the high utilization potential of deep learning for the PUB by demonstrating the new state-of-the-art performance in this field.

A dynamic water quality model, HYPE, was applied to a large, data-sparse region to study whether reliable information on water quantity and water quality could be obtained for both gauged and ungauged waterbodies. The model (called S-HYPE) was set up for all of Sweden (∼450 000 km2), divided into sub-basins with an average area of 28 km2. Readily available national databases were used for physiographic data, emissions and agricultural practices, fixed values for representative years were used. Daily precipitation and temperature were used as the dynamic forcing of the model. Model evaluation was based on data from several hundred monitoring sites, of which approximately 90% had not been used in calibration on a daily scale. Results were evaluated using the Nash-Sutcliffe efficiency (NSE), correlation and relative errors: 92% of the spatial variation was explained for specific water discharge, and 88% and 59% for total nitrogen and total phosphorus concentrations, respectively. Day-to-day variations were modelled with satisfactory results for water discharge and the seasonal variation of nitrogen concentrations was also generally well captured. In 20 large, unregulated rivers the median NSE for water discharge was 0.84, and the corresponding number for 76 partly-regulated river basins was 0.52. In small basins, the NSE was typically above 0.6. These major achievements relative to previous similar experiments were ascribed to the step-wise calibration process using representative gauged basins and the use of a modelling concept, whereby coefficients are linked to physiographic variables rather than to specific sites. Editor D. Koutsoyiannis Citation Strömqvist, J., Arheimer, B., Dahné, J., Donnelly, C. and Lindström, G., 2012. Water and nutrient predictions in ungauged basins: set-up and evaluation of a model at the national scale. Hydrological Sciences Journal, 57 (2), 229–247.

Calibrating a hydrological model in ungauged small river basins of the northeastern Tibetan Plateau based on near-infrared images

A new parameter parsimonious rainfall–run‐off model, the Distance Distribution Dynamics (DDD) model, is used to simulate hydrological time series at ungauged sites in the Lygne basin in Norway. The model parameters were estimated as functions of catchment characteristics determined by geographical information system. The multiple regression equations relating catchment characteristics and model parameters were trained from 84 calibrated catchments located all over Norway, and all model parameters showed significant correlations with catchment characteristics. The significant correlation coefficients (with p‐value < 0.05) ranged from 0.22 to 0.55. The suitability of DDD for predictions in ungauged basins was tested for 17 catchments not used to estimate the multiple regression equations. For ten of the 17 catchments, deviations in Nash–Sutcliffe efficiency (NSE) criteria between the calibrated and regionalised model were less than 0.1, and for two calibrated catchments within the Lygne basin, the deviations were less than 0.08. The median NSE for the regionalized DDD for the 17 catchments for two time series was 0.66 and 0.72. Deviations in NSE between calibrated and regionalised models are well explained by the deviations between calibrated and regressed parameters describing spatial snow distribution and snowmelt respectively. The quality of the simulated run‐off series for the ungauged sites in the Lygne basin was assessed by comparing flow indices describing high, medium and low flow estimated from observed run‐off at the 17 catchments and for the simulated run‐off series. The indices estimated for the simulated series were generally well within the ranges defined by the 17 observed series. Copyright © 2014 John Wiley & Sons, Ltd.

Quantifying continental-scale river discharge is essential to understanding the terrestrial water cycle but is susceptible to errors caused by a lack of observations and the limitations of hydrodynamic modeling. Data assimilation (DA) methods are increasingly used to estimate river discharge in combination with emerging river-related remote sensing products (e.g., water surface elevation [WSE], water surface slope, river width, and flood extent). However, directly comparing simulated WSE to satellite altimetry data remains challenging (e.g., because of large biases between simulations and observations or uncertainties in parameters), and large errors can be introduced when satellite observations are assimilated into hydrodynamic models. In this study we performed direct, anomaly, and normalized value assimilation experiments to investigate the capacity of DA to improve river discharge within the current limitations of hydrodynamic modeling. We performed hydrological DA using a physically-based empirical localization method applied to the Amazon Basin. We used satellite altimetry data from ENVISAT, Jason 1, and Jason 2. Direct DA was the baseline assimilation method and was subject to errors due to biases in the simulated WSE. To overcome these errors, we used anomaly DA as an alternative to direct DA. We found that the modeled and observed WSE distributions differed considerably (e.g., differences in amplitude, seasonal flow variation, and a skewed distribution due to limitations of the hydrodynamic models). Therefore, normalized value DA was performed to improve discharge estimation. River discharge estimates were improved at 24 %, 38 %, and 62 % of stream gauges in the direct, anomaly, and normalized value assimilations relative to simulations without DA. Normalized value assimilation performed best for estimating river discharge given the current limitations of hydrodynamic models. Most gauges within the river reaches covered by satellite observations accurately estimated river discharge, with Nash-Sutcliffe efficiency (NSE) > 0.6. The amplitudes of WSE variation were improved in the normalized DA experiment. Furthermore, in the Amazon Basin, normalized assimilation (median NSE = 0.50) improved river discharge estimation compared to open-loop simulation with the global hydrodynamic model (median NSE = 0.42). River discharge estimation using direct DA methods was improved by 7 % with calibration of river bathymetry based on NSE. The direct DA approach outperformed the other DA approaches when runoff was considerably biased, but anomaly DA performed best when the river bathymetry was erroneous. The uncertainties in hydrodynamic modeling (e.g., river bottom elevation, river width, simplified floodplain dynamics, and the rectangular cross-section assumption) should be improved to fully realize the advantages of river discharge DA through the assimilation of satellite altimetry. This study contributes to the development of a global river discharge reanalysis product that is consistent spatially and temporally.

Abstract. Non-parametric Bayesian networks (NPBNs) are graphical tools for statistical inference widely used for reliability analysis and risk assessment and present several advantages, such as the embedded uncertainty quantification and limited computational time for the inference process. However, their implementation in hydrological studies is still scarce. Hence, to increase our understanding of their applicability and extend their use in hydrology, we explore the potential of NPBNs to reproduce catchment-scale hydrological dynamics. Long-term data from 240 river catchments with contrasting climates across the United States from the Catchment Attributes and Meteorology for Large-sample Studies (CAMELS) data set will be used as actual means to test the utility of NPBNs as descriptive models and to evaluate them as predictive models for maximum daily river discharge in any given month. We analyse the performance of three networks, one unsaturated (hereafter UN-1), one saturated (hereafter SN-1), both defined only by hydro-meteorological variables and their bivariate correlations, and one saturated network (hereafter SN-C), consisting of the SN-1 network and including physical catchments' attributes. The results indicate that the UN-1 network is suitable for catchments with a positive dependence between precipitation and discharge, while the SN-1 network can also reproduce discharge in catchments with negative dependence. The latter can reproduce statistical characteristics of discharge (tested via the Kolmogorov–Smirnov statistic) and have a Nash–Sutcliffe efficiency (NSE) ≥0.5 in ∼40 % of the catchments analysed, receiving precipitation mainly in winter and located in energy-limited regions at low to moderate elevation. Further, the SN-C network, based on similarity of the catchments, can reproduce discharge statistics in ∼10 % of the catchments analysed. We show that once a NPBN is defined, it is straightforward to infer discharge and to extend the network itself with additional variables, i.e. going from the SN-1 network to the SN-C network. However, the results also suggest considerable challenges in defining a suitable NPBN, particularly for predictions in ungauged basins. These are mainly due to the discrepancies in the timescale of the different physical processes generating discharge, the presence of a “memory” in the system, and the Gaussian-copula assumption used for modelling multivariate dependence.

Large-sample hydrology datasets have become increasingly popular in recent years, providing hydro-meteorological time series and catchment attributes for thousands of catchments worldwide. However, the role of such catchment attributes in informing model regionalization, and particularly the effect of their level of spatial detail on prediction in ungauged basins (PUB), remains poorly explored. This study addresses this gap by examining whether catchment attributes derived from geological maps of varying levels of detail improve model regionalization and, in turn, PUB, for both a bucket-type model and a data-driven Long Short-Term Memory (LSTM) model across 130 catchments in two independent basins: the Moselle (27 100 km²) and the Garonne (13 730 km²). We conducted five modeling experiments: a benchmark without geological information and four geology-informed configurations with increasing levels of detail (random, global, continental, and regional). A fold-based space–time cross-evaluation strategy was used to assess model performance on both time periods and catchments unseen during calibration. Performance was evaluated using a modified Nash–Sutcliffe Efficiency (NSE) and a set of streamflow signatures describing flow variability, storage, and regime behavior. Across both basins and model types, benchmark experiments yielded the lowest space–time performance, followed by the random experiment and the global geology experiment, while the experiments using continental and regional geology consistently resulted in higher NSE values. Improvements were strongest for the bucket-type model, with the most detailed geological attributes leading to consistent gains in median performance and robustness. Differences among experiments were more pronounced for streamflow signatures. For the bucket-type model, only the experiments adopting the continental or regional geology reproduced observed signatures with Spearman’s correlations exceeding 0.60, whereas the LSTM model already showed reasonable skill in the benchmark case but still benefited systematically from increasing geological detail. Together, these results demonstrate that incorporating detailed geological information can enhance streamflow representation and model transferability in PUB applications, and that the level of geological detail is a critical, yet often overlooked, factor in large-sample hydrology and regionalization studies.

Reduction of uncertainty in hydrologic predictions is important for water resources development and management, especially in ungauged basins. The objective of this study was to predict rainfall and runoff in an ungauged basin located in a mountainous region of northern Thailand by using a weather radar to estimate rainfall and a modeling approach to simulate rainfall-runoff processes. The calibrated Ze-R relationship (Ze=18R1.45) was used to estimate rainfall over a gauged basin and an assumed ungauged basin. The runoff simulations in the gauged and ungauged basins were conducted using the quasi-distributed hydrological model HEC-HMS. The simulation results in the gauged basin, for both calibration and verification, were in agreement with the observed data. Model parameters in the assumed ungauged basin were estimated using transposition and regionalization techniques. The predicted direct runoff hydrograph in the assumed ungauged basin using radar measured rainfall and the regionalization technique matches well with the measurements at the basin outlet. Based on the study results, the radar rainfall estimates and regionalization of the hydrologic model parameters are promising alternatives for improved hydrologic predictions in ungauged basins.

Prediction in Ungauged Basins (PUB) is an important task for water resources planning and management and remains a fundamental challenge for the hydrological community. In recent years, geostatistical methods have proven valuable for estimating hydrological variables in ungauged catchments. However, four major problems restrict the development of geostatistical methods. We established a new information diffusion model based on genetic algorithm (GIDM) for spatial interpolating of runoff in the ungauged basins. Genetic algorithms (GA) are used to generate high-quality solutions to optimization and search problems. So, using GA, the parameter of optimal window width can be obtained. To test our new method, seven experiments for the annual runoff interpolation based on GIDM at 17 stations on the mainstream and tributaries of the Yellow River are carried out and compared with the inverse distance weighting (IDW) method, Cokriging (COK) method, and conventional IDMs using the same sparse observed data. The seven experiments all show that the GIDM method can solve four problems of the previous geostatistical methods to some extent and obtains best accuracy among four different models. The key problems of the PUB research are the lack of observation data and the difficulties in information extraction. So the GIDM is a new and useful tool to solve the Prediction in Ungauged Basins (PUB) problem and to improve the water management.

The focus in the search for more reliable predictions in ungauged basins (PUB) has generally been on reducing uncertainty in watershed models (mainly their parameters). More recently, however, we seem to remember that the ultimate objective is not to define the parameters of a specific model but to understand the watershed: What behavior do we expect the ungauged watershed to exhibit? And what behavior should not occur in a particular ungauged watershed? The answers to these questions actually provide additional information that can be assimilated in watershed models for uncertainty reduction in PUB. This extension to hydrologic modeling approaches provides a quantitative link between watershed modeling and statistical hydrology as well as process hydrology that has to be explored. We witness a convergence of approaches—Bayesian, set theoretic, and optimization based—toward utilizing this link. The result is an opportunity for the (quantitative) dialog between modelers, statistical hydrologists, and experimentalists. We close our discussion of this development by presenting new and exciting research questions that we now have to address.

Study RegionThe study area consists of 44 catchments across Ireland. Study FocusWe regionalize two hydrological models (GR4J and GR6J) to produce continuous discharge simulations and compare performance in simulating high, median and low flow conditions with other established approaches to prediction in ungauged basins and a simple benchmark of using the median parameter set across all catchments. These include K-nearest neighbor (KNN) and statistical methods for predicting flow quantiles using catchment characteristics. Different objective functions were selected for different parts of flow regime and the success of different methods for regionalizing hydrological model parameters; including multiple linear regression (MLR), non-linear regression (NL) and random forests (RF) were evaluated. New Hydrological insights for the RegionAll regionalization approaches perform well for average flow conditions. The GR4J model regionalized using RF performs best for simulating high flows, though all regionalized models underestimate the median annual flood. GR6J regionalized using RF performs best for low flows. While KNN and statistical approaches that directly leverage physical catchment descriptors provide comparable median performances across catchments, the spread in relative error across our sample is reduced using regionalized hydrological models. Our results highlight that the choice of hydrological model, objective functions for optimization and approach to linking model parameters and physical catchment descriptors significantly influence the success of regionalization for low and high flows.

Modelling ungauged catchments using the catchment runoff response similarity

A computer vision-based approach to fusing spatiotemporal data for hydrological modeling