Reply on RC2

This study investigates the potential predictability of streamflow and soil moisture in the Alabama–Coosa–Tallapoosa (ACT) river basin in the southeastern United States. The study employs the state-of-the-art National Water Model (NWM) and compares the effects of initial soil moisture condition with those of seasonal climate anomalies on streamflow and soil moisture forecast skills. We have designed and implemented seasonal streamflow forecast ensemble experiments following the methodology suggested by Dirmeyer et al. The study also compares the soil moisture variability in the NWM with in situ measurements and remote sensing data from the Soil Moisture Active and Passive (SMAP) satellite. The NWM skillfully simulates the observed streamflow in the ACT basin. The soil moisture variability is 46% smaller in the NWM compared with the SMAP data, mainly due to a weaker amplitude of the seasonal cycle. This study finds that initial soil moisture condition is a major source of predictability for the seasonal streamflow forecast. The contribution of the initial soil moisture condition is comparable or even higher than that of seasonal climate anomaly effects in dry seasons. In the boreal summer season, the initial soil moisture condition contributes to 65% and 48% improvements in the seasonal streamflow and soil moisture forecast skills, respectively. This study attributes a greater improvement in the streamflow forecast skill to the lag effects between the soil moisture and streamflow anomalies. The results of this study can inform the development and improvement of the operational streamflow forecasting system.

Abstract. Skillful seasonal streamflow forecasts obtained from climate and land surface conditions could significantly improve water and energy management. Since climate forecasts are updated on a monthly basis, we evaluate the potential in developing operational monthly streamflow forecasts on a continuous basis throughout the year. Further, basins in the rainfall–runoff regime critically depend on the forecasted precipitation in the upcoming months as opposed to snowmelt regimes where initial hydrological conditions (IHC) play a critical role. The goal of this study is to quantify the role of updated monthly precipitation forecasts and IHC in forecasting 6-month lead monthly streamflow and soil moisture for a rainfall–runoff mechanism dominated basin – Apalachicola River at Chattahoochee, FL. The Variable Infiltration Capacity (VIC) land surface model is implemented with two forcings: (a) updated monthly precipitation forecasts from ECHAM4.5 Atmospheric General Circulation Model (AGCM) forced with sea surface temperature forecasts and (b) daily climatological ensembles. The difference in skill between the above two quantifies the improvements that could be attainable using the AGCM forecasts. Monthly retrospective streamflow forecasts are developed from 1981 to 2010 and streamflow forecasts estimated from the VIC model are also compared with those predicted by using the principal component regression (PCR) model. The mean square error (MSE) in predicting monthly streamflows, using the VIC model, are compared with the MSE of streamflow climatology under ENSO (El Niño Southern Oscilation) conditions as well as under normal years. Results indicate that VIC forecasts obtained using ECHAM4.5 are significantly better than VIC forecasts obtained using climatological ensembles and PCR models over 2–6 month lead time during winter and spring seasons in capturing streamflow variability and reduced mean square errors. However, at 1-month lead time, streamflow utilizing the climatological forcing scheme outperformed ECHAM4.5 based streamflow forecasts during winter and spring, indicating a dominant role of IHCs up to a 1-month lead time. During ENSO years, streamflow forecasts exhibit better skill even up to a six-month lead time. Comparisons of the seasonal soil moisture forecasts, developed using ECHAM4.5 forcings, with seasonal streamflows also show significant skill, up to a 6-month lead time, in the four seasons.

Society’s increasing demand for water and the need for its long-term management have motivated efforts toward improving seasonal streamflow forecasts. Currently, seasonal climate forecasts are routinely issued in meteorological centers around the world, generating information for decision-making and seasonal streamflow forecasting (SSF) studies that are becoming more frequent. Seasonal streamflow forecast skill derives from land surface initial conditions and atmospheric boundary conditions that mostly depend on large-scale climate phenomena (such as ENSO). Thus, seasonal rainfall predictions produced by dynamic climate models that represent ocean-atmosphere interactions may have a positive impact on streamflow forecasts. In South America, seasonal streamflow forecasts are essential for the hydropower sector, which is responsible for ~65% of the electric energy produced in countries such as Brazil. In this work, we assessed seasonal streamflow forecasts over South America based on a continental-scale application of a hydrologic-hydrodynamic model and precipitation forecasts from the ECMWF's fifth generation seasonal forecast system (SEAS5). Seasonal streamflow forecasts (SEAS5-SF) were evaluated against a reference model run and forecast skill was estimated relative to the Ensemble Streamflow Prediction (ESP) method. The bias correction of SEAS5 predicted precipitation improved the performance of the seasonal streamflow forecasts, frequently turning negative skill results into near null to positive skill. Results indicate that the ESP remains a hard-to-beat method for seasonal streamflow forecasting in South America. SEAS5-SF skill was found to be dependent on initialization month, season, basin and forecast lead time, with greater skill on the initialization month lead time. Rivers where the forecast skill is higher were Amazon, Araguaia, Tocantins and Paraná. Acknowledgments: This work presents part of the results obtained during the project granted by the Brazilian Agency of Electrical Energy (ANEEL) under its Research and Development program Project PD 6491-0503/2018 – “Previsão Hidroclimática com Abrangência no Sistema Interligado Nacional de Energia Elétrica” developed by the Paraná State electric company (COPEL GeT), the Meteorological System of Paraná (SIMEPAR) and the RHAMA Consulting company. The Hydraulic Research Institute (IPH) from the Federal University of Rio Grande do Sul (UFRGS) contribute to part of the project through an agreement with the RHAMA company (IAP-001313).

Use of an Ensemble Re-ordering Method for disaggregation of seasonal categorical rainfall forecasts into conditioned ensembles of daily rainfall for hydrological forecasting

<p>Streamflow information for the months ahead is of great value to existing decision-making practices, particularly to those affected by the vagaries of the climate and who would benefit from better understanding and managing climate-related risks. Despite the large effort, there is still limited knowledge of the key drivers controlling the quality of the seasonal streamflow forecasts. In this investigation, we show that the seasonal streamflow predictability can be clustered, and hence regionalised, based on a priori knowledge of local hydro-climatic conditions. To reach these conclusions we analyse the seasonal forecasts of streamflow volumes across about 35400 basins in Europe, which vary in terms of climatology, scale and hydrological regime. We then link the forecast quality to various descriptors including physiography, hydro-climatic characteristics and meteorological biases. This allows the identification of the key drivers along a strong hydro-climatic gradient. Results show that, as expected, the seasonal streamflow predictability varies geographically and seasonally with acceptable values for the first lead months. In addition, the predictability deteriorates with increasing lead months particularly in the winter months. Nevertheless, we show that the forecast quality is well correlated to a set of drivers, which vary depending on the initialization month. The forecast quality of seasonal streamflow volumes is strongly dependent on the basin’s hydrological regime, with quickly reacting basins (of low river memory) showing limited predictability. On the contrary, snow and/or baseflow dominated regions with long recessions (and hence high river memory) show high streamflow predictability. Finally, climatology and precipitation biases are also strongly related to streamflow predictability, highlighting the importance of developing robust bias-adjustment methods.</p>

We provide a framework for careful analysis of the different methodological choices we make when constructing multimodel ensemble seasonal forecasts of hydroclimatic variables. Specifically, we focus on three common modeling decisions: (i) number of models, (ii) multimodel combination approach, and (iii) lead time for prediction. The analysis scheme includes a multimodel ensemble forecasting algorithm based on nonparametric regression, a set of alternatives for the options previously pointed, and a selection of probabilistic verification methods for ensemble forecast evaluation. The usefulness of this framework is tested through an example application aimed to generate spring/summer streamflow forecasts at multiple locations in Central Chile. Results demonstrate the high impact that subjectivity in decision‐making may have on the quality of ensemble seasonal hydroclimatic forecasts. In particular, we note that the probabilistic verification criteria may lead to different choices regarding the number of models or the multimodel combination method. We also illustrate how this objective analysis scheme may lead to results that are extremely relevant for the case study presented here, such as skillful seasonal streamflow predictions for very dry conditions.

[1] Seasonal streamflow forecasts are valuable for planning and allocation of water resources. In Australia, the Bureau of Meteorology employs a statistical method to forecast seasonal streamflows. The method uses predictors that are related to catchment wetness at the start of a forecast period and to climate during the forecast period. For the latter, a predictor is selected among a number of lagged climate indices as candidates to give the “best” model in terms of model performance in cross validation. This study investigates two strategies for further improvement in seasonal streamflow forecasts. The first is to combine, through Bayesian model averaging, multiple candidate models with different lagged climate indices as predictors, to take advantage of different predictive strengths of the multiple models. The second strategy is to introduce additional candidate models, using rainfall and sea surface temperature predictions from a global climate model as predictors. This is to take advantage of the direct simulations of various dynamic processes. The results show that combining forecasts from multiple statistical models generally yields more skillful forecasts than using only the best model and appears to moderate the worst forecast errors. The use of rainfall predictions from the dynamical climate model marginally improves the streamflow forecasts when viewed over all the study catchments and seasons, but the use of sea surface temperature predictions provide little additional benefit.

Droughts, prolonged heatwaves, heavy rainfall and multiple large-scale floods - recent years have shown that global climate change requires more sustainable and timely water management at all levels. In particular, for the optimized use of water resources for irrigation or hydropower generation, it is essential to know their likely availability in the coming months anywhere in the world. This sub-seasonal to seasonal time range is covered by seasonal forecasting systems such as SEAS5 developed by the European Center for Medium-Range Weather Forecasts (ECMWF). These systems have the potential to provide important data for improving water management practices. However, without a bias correction, the data deviate strongly from the actual data. We have shown for several regions of the world that the Bias Correction and Spatial Disaggregation (BCSD) method can significantly improve predictive capability. By storing fixed cumulative distribution functions (CDFs) and parallelization, we were able to extend the system from the regional to the global level, i.e. to produce BCSD predictions for the entire globe and present this version at EGU24.Our next step is to evaluate the resulting seasonal forecasts in terms of their predictive quality. This evaluation is carried out using several measures of quality, including the Continuous Ranked Probability Skill Score (CRPSS) and the Brier Skill Score (BSS). To illustrate this, the strengths and weaknesses of the bias-corrected seasonal forecasts are highlighted using two regions. The focus is on the Sahel region, which has a lower forecast quality despite its high social relevance, and on the Lake Victoria catchment area, for which a high forecast quality is achieved. The aim is to achieve as precise an assessment as possible of the global forecast quality, which allows a realistic assessment of the forecasts, particularly in regions with strong fluctuations in water availability. By providing this bias-corrected forecast data in near real time together with an analysis of its quality, better estimates will be available for direct use by water managers or as input for subsequent modeling processes.

In recent years, alternating drought and extreme precipitation events have highlighted the need for subseasonal to seasonal forecasts of the terrestrial water cycle. In particular, predictions of the impacts of dry and wet extremes on subsurface water resources are crucial to provide stakeholders in agriculture, forestry, the water sector, and other fields with information supporting the sustainable use of these resources.In this context, we release an experimental Water Resources Bulletin (https://adapter-projekt.de/bulletin/index.html) four times per year, offering probabilistic forecasts of the total subsurface water storage (TSS) anomaly at a 0.6 km resolution, from the surface down to 60 m depth, for the upcoming seven months across Germany. These seasonal forecasts are generated using the integrated, physics-based hydrological model ParFlow/CLM, forced by 50 ensemble members of the SEAS5 seasonal forecast from the European Centre for Medium-Range Weather Forecasts (ECMWF).To evaluate our forecasts, we evaluated six 7-months probabilistic forecasts covering the vegetation period (March to September) for the years 2018 to 2023 with a reference long-term historical time series based on the same ParFlow/CLM setup. The forecast skill was assessed by comparing these seasonal forecasts to a climatology-based 10-member pseudo-forecast over the 2013–2023 period (using the leave-one-out method), extracted from the reference time series.The monthly Continuous Ranked Probability Skill Score (CRPSS), which evaluates the ensemble distribution based on daily TSS data, indicates that the probabilistic forecast outperforms the climatology-based pseudo-forecast in most regions, except in 2018 and, to a lesser extent, in 2020 and 2022. This can be attributed to an under-representation of extremely dry members in the ensemble, combined with the memory effect of the initial conditions at increasing soil depths. For example, while March 2018 started with a slightly above-average TSS and experienced a strong meteorological drought leading to an agricultural drought, the initial TSS anomaly in March 2019 was already negative, with a less pronounced precipitation deficit during the vegetation period. This resulted in a much higher forecast skill, because of the memory effect accurately simulated with the physics-based model. Notably, the forecast skill only slightly decreases with increasing lead time, both for precipitation and TSS.The analysis of the Relative Operating Characteristic Skill Score (ROCSS) for the lower quintile of the TSS distribution assesses whether negative TSS anomalies (i.e., droughts) are adequately represented within the probabilistic forecast ensemble. The results are consistent with those of the CRPSS, showing lower skill in 2018. Nevertheless, the ROCSS analysis overall indicates moderate to high skill for the probabilistic forecast, while the climatology-based pseudo-forecast demonstrates no skill. This confirms that the dry conditions experienced in central Europe in recent years were captured within the probabilistic forecast, underlining the added value of these forecasts and their usefulness in the experimental Water Resources Bulletin.

This study asks the question of whether GCMs are ready to be operationalized for streamflow forecasting in South Asian river basins, and if so, at what temporal scales and for which water management decisions are they likely to be relevant? The authors focused on the Ganges, Brahmaputra, and Meghna basins for which there is a gridded hydrologic model calibrated for the 2002–10 period. The North American Multimodel Ensemble (NMME) suite of eight GCM hindcasts was applied to generate precipitation forecasts for each month of the 1982–2012 (30 year) period at up to 6 months of lead time, which were then downscaled according to the bias-corrected statistical downscaling (BCSD) procedure to daily time steps. A global retrospective forcing dataset was used for this downscaling procedure. The study clearly revealed that a regionally consistent forcing for BCSD, which is currently unavailable for the region, is one of the primary conditions to realize reasonable skill in streamflow forecasting. In terms of relative RMSE (normalized by reference flow obtained from the global retrospective forcings used in downscaling), streamflow forecast uncertainty (RMSE) was found to be 38%–50% at monthly scale and 22%–35% at seasonal (3 monthly) scale. The Ganges River (regulated) experienced higher uncertainty than the Brahmaputra River (unregulated). In terms of anomaly correlation coefficient (ACC), the streamflow forecasting at seasonal (3 monthly) scale was found to have less uncertainty (>0.3) than at monthly scale (<0.25). The forecast skill in the Brahmaputra basin showed more improvement when the time horizon was aggregated from monthly to seasonal than the Ganges basin. Finally, the skill assessment for the individual seasons revealed that the flow forecasting using NMME data had less uncertainty during monsoon season (July–September) in the Brahmaputra basin and in postmonsoon season (October–December) in the Ganges basin. Overall, the study indicated that GCMs can have value for management decisions only at seasonal or annual water balance applications at best if appropriate historical forcings are used in downscaling. The take-home message of this study is that GCMs are not yet ready for prime-time operationalization for a wide variety of multiscale water management decisions for the Ganges and Brahmaputra River basins.

Despite recent advances in operational streamflow forecasting systems, anticipating and forecasting droughts and associated low-flow conditions remain major challenges in hydrology, with substantial impacts on water-dependent sectors such as agriculture and industry. Enhancing sub-seasonal to seasonal streamflow forecasting is therefore critical for improving water resources management. This study investigates the performance of a Handoff forecast Long Short-Term Memory (LSTM) architecture (Nearing et al., 2024) for probabilistic streamflow forecasting at lead times extending up to six months, with a particular emphasis on low-flow conditions.The Handoff forecast LSTM is trained regionally on a subset of 292 basins from the Catchment Attributes and MEteorology for Large-sample Studies - FRance dataset (CAMELS-FR) (Delaigue et al., 2025), after excluding basins affected by unreliable low-flow measurements. Model training relies on basin-averaged hydro-meteorological reanalysis data provided by CAMELS-FR. Evaluation of the model is conducted using ensemble streamflow forecasts generated from historical scenarios and meteorological ensemble predictions from the SEAS5 model from the European Center for Medium-Range Weather Forecasts (ECMWF) (Johnson et al., 2019)The generated ensemble streamflow forecasts are evaluated using a set of probabilistic metrics such as the Continuous Ranked Probability Score (CRPS), the Brier Score, the Area under the ROC curve, and the Talagrand diagram, and using the natural streamflow climatology as a reference. In addition, a sensitivity analysis of static catchment attributes is performed to assess their relative contribution to model performance and to better understand the drivers of predictability across basins.Delaigue, O., Guimarães, G. M., Brigode, P., Génot, B., Perrin, C., Soubeyroux, J.-M., Janet, B., Addor, N., & Andréassian, V. (2025). CAMELS-FR dataset: a large-sample hydroclimatic dataset for France to explore hydrological diversity and support model benchmarking. Earth System Science Data, 17(4), 1461–1479. https://doi.org/10.5194/essd-17-1461-2025Johnson, S. J., Stockdale, T. N., Ferranti, L., Balmaseda, M. A., Molteni, F., Magnusson, L., Tietsche, S., Decremer, D., Weisheimer, A., Balsamo, G., Keeley, S. P. E., Mogensen, K., Zuo, H., & Monge-Sanz, B. M. (2019). SEAS5: The new ECMWF seasonal forecast system. Geoscientific Model Development, 12(3), 1087–1117. https://doi.org/10.5194/gmd-12-1087-2019Nearing, G., Cohen, D., Dube, V., Gauch, M., Gilon, O., Harrigan, S., Hassidim, A., Klotz, D., Kratzert, F., Metzger, A., Nevo, S., Pappenberger, F., Prudhomme, C., Shalev, G., Shenzis, S., Tekalign, T. Y., Weitzner, D., & Matias, Y. (2024). Global prediction of extreme floods in ungauged watersheds. Nature, 627(8004), 559–563. https://doi.org/10.1038/s41586-024-07145-

Climate change is expected to make droughts more frequent and severe all over the world. South Korea is no exception and is already suffering from extreme droughts such as one that prolonged from 2013 to 2015 and caused nation‐wide damage. To mitigate drought damages, better management of existing water infrastructure is essential. A promising opportunity to improve operational decisions is to make use of seasonal weather forecasts that are provided by general circulation models and can be downscaled to the catchment scale. This study hence assesses the skill of seasonal forecasts over 20 catchments in South Korea, where the largest reservoirs are located. Datasets from four weather forecasting centres (ECMWF, UK Met Office, Météo France and DWD) have been evaluated over the period 2011–2020, and their skill quantified using the Continuous Ranked Probability Skill Score (CRPSS). We analyse how skill varies across the seasons and years, and if it can be linked to catchments characteristics. In doing so, we develop a methodology and a Python package to implement it, which is freely available for future applications to other regions. As for the study case, our results showed that among the four forecasting centres, ECMWF's forecasts were the most skilful in South Korea. In particular, seasonal forecasts outperform the climatology for 2 months of lead time and are more skilful during the wet season and in dry years. Linear bias correction is found to be useful to correct systematic seasonal biases, whereas we found no significant correlation between the catchment characteristics and forecast skill. We also investigated the possibility of anticipating dry years from seasonal forecasts and/or ENSO indices but found no significant link.

We investigate the ability of hydrological multimodel ensemble predictions to enhance the skill of streamflow forecasts at short‐ to medium‐range timescales. To generate the multimodel ensembles, we implement a new statistical postprocessor, namely, quantile regression‐Bayesian model averaging (QR‐BMA). Quantile regression‐Bayesian model averaging uses quantile regression to bias correct the ensemble streamflow forecasts from the individual models and Bayesian model averaging to optimally combine their probability density functions. Additionally, we use an information‐theoretic measure, namely, conditional mutual information, to quantify the skill enhancements from the multimodel forecasts. We generate ensemble streamflow forecasts at lead times from 1 to 7 days using three hydrological models: (i) Antecedent Precipitation Index‐Continuous, (ii) Hydrology Laboratory‐Research Distributed Hydrologic Model, and (iii) Weather Research and Forecasting Hydrological modeling system. As forcing to the hydrological models, we use weather ensemble forecasts from the National Centers for Environmental Prediction 11‐member Global Ensemble Forecast System Reforecast version 2. The forecasting experiments are performed for four nested basins of the North Branch Susquehanna River, USA. We find that after bias correcting the streamflow forecasts from each model, their skill performance becomes comparable. We find that the multimodel ensemble forecasts have higher skill than the best single‐model forecasts. Furthermore, the skill enhancements obtained by the multimodel ensemble forecasts are found to be dominated by model diversity, rather than by increased ensemble size alone. This result, obtained using conditional mutual information, indicates that each hydrological model contributes additional information to enhance forecast skill. Overall, our results highlight benefits of hydrological multimodel forecasting for improving streamflow predictions.

Forecast verification is important across scientific disciplines, as it provides a framework for evaluating the performance of a forecasting system. In the atmospheric sciences, probabilistic skill scores are often used for verification, as they provide a way of ranking the performance of different probabilistic forecasts unambiguously. In order to be useful, a skill score must be proper: it must encourage honesty in the forecaster and reward forecasts that are reliable and have good resolution. A new score, the error‐spread score (ES), is proposed, which is particularly suitable for evaluation of ensemble forecasts. It is formulated with respect to the moments of the forecast. The ES is confirmed to be a proper score and is therefore sensitive to both resolution and reliability. The ES is tested on forecasts made using the Lorenz ′96 system and found to be useful for summarizing the skill of the forecasts. The European Centre for Medium‐Range Weather Forecasts (ECMWF) ensemble prediction system (EPS) is evaluated using the ES. Its performance is compared with a perfect statistical probabilistic forecast: the ECMWF high‐resolution deterministic forecast dressed with the observed error distribution. This generates a forecast that is perfectly reliable if considered over all time, but does not vary from day to day with the predictability of the atmospheric flow. The ES distinguishes between the dynamically reliable EPS forecasts and the statically reliable dressed deterministic forecasts. Other skill scores are tested and found to be comparatively insensitive to this desirable forecast quality. The ES is used to evaluate seasonal range ensemble forecasts made with the ECMWF System 4. The ensemble forecasts are found to be skilful when compared with climatological or persistence forecasts, though this skill is dependent on the region and time of year.

Accurate and reliable seasonal climate forecasts are frequently sought by climate‐sensitive sectors to support decision‐making under climate variability and change. Temperature trend is discernible globally over the past decades, but seasonal forecasts produced by a global climate model (GCM) generally underestimate such trend. Current statistical methods used for calibrating seasonal climate forecasts mostly do not explicitly account for climate trends. Consequently, the calibrated forecasts also fail to capture the observed trend. Solving this problem can enhance user confidence in seasonal climate forecasts. In this study, we extend the capability of the Bayesian joint probability (BJP) modelling approach for statistical calibration of seasonal climate forecasts. A trend component is introduced into the BJP algorithm for embedding the observed trend into calibrated ensemble forecasts. We apply the new model (named BJP‐t) to three test stations in Australia. Seasonal forecasts of daily maximum temperatures from the SEAS5 model, operated by the European Centre for Medium‐Range Weather Forecasts (ECMWF), are calibrated and evaluated. The BJP‐t calibrated ensemble forecasts can reproduce the observed trend, when the raw ensemble forecasts and the BJP calibrated ensemble forecasts both fail to do so. The BJP‐t calibration leads to more skilful, more reliable and sharper forecasts than the BJP calibration.