Next Article in Journal
Recent Advances in Biofiltration for PPCP Removal from Water
Next Article in Special Issue
Evaluation of Future Changes in Climate Extremes over Southeast Asia Using Downscaled CMIP6 GCM Projections
Previous Article in Journal
Under the Strong Influence of Human Activities: The Patterns and Controlling Factors of River Water Chemistry Changes—A Case Study of the Lower Yellow River
Previous Article in Special Issue
Risk of Natural Hazards Caused by Extreme Precipitation in Poland in 1951–2020
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving Flood Forecasting Skill by Combining Ensemble Precipitation Forecasts and Multiple Hydrological Models in a Mountainous Basin

1
China Meteorological Administration Basin Heavy Rainfall Key Laboratory/Hubei Key Laboratory for Heavy Rain Monitoring and Warning Research, Institute of Heavy Rain, China Meteorological Administration, Wuhan 430205, China
2
State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China
3
Three Gorges National Climatological Observatory, Yichang 443099, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(13), 1887; https://doi.org/10.3390/w16131887
Submission received: 11 May 2024 / Revised: 20 June 2024 / Accepted: 20 June 2024 / Published: 1 July 2024
(This article belongs to the Special Issue Analysis of Extreme Precipitation Under Climate Change)

Abstract

:
Ensemble precipitation forecasts (EPFs) derived from single numerical weather predictions (NWPs) often miss extreme events, and individual hydrological models (HMs) often fail to accurately capture all types of flows, including flood peaks. To address these shortcomings, this study introduced four “EPF + HM” schemes for ensemble flood forecasting (EFF) by combining two EPFs and two HMs. A generator-based post-processing (GPP) method was applied to correct biases and under-dispersion within the raw EPF data. The effectiveness of these schemes in delivering high-quality flood forecasts was assessed using both deterministic and probabilistic metrics. The results indicate that, once post-processed by GPP, all proposed schemes show improvements in both deterministic and probabilistic performances, with skillful flood forecasts for 1–7 lead days. The deterioration in forecast performance with extended lead times is also lessened. Notably, the results indicate that uncertainty within hydrological models has a more pronounced impact on capturing flood peaks than uncertainty in precipitation inputs. This study recommends combining individual EPF with multiple hydrological models for reliable flood forecasting. In conclusion, effective flood forecasting necessitates employing post-processing techniques to correct EPFs and accounting for the uncertainty inherent in hydrological models, rather than relying solely on the uncertainty of the input data.

1. Introduction

Floods are the most frequent natural disasters worldwide, posing significant threats to human lives and socio-economic development [1]. Research by Alfieri, Bisselink [2] indicates that intensifying hydrological cycles across global watersheds are causing large-scale flooding events, often exceeding the capacity of existing flood protection infrastructure, a trend likely exacerbated by climate change. Since the late 1980s, flood management strategies have shifted from reliance on structural measures to embracing non-structural approaches [3,4]. Among these, flood forecasting and early warning systems are considered the most effective and cost-efficient methods for mitigating flood risks and minimizing their impacts [5,6,7]. These systems provide critical information regarding water flow, levels, and the timing of peak floods. Timely warnings enable authorities and affected populations to implement effective response measures, underscoring the need to improve forecast accuracy and extend lead times.
Flood forecasting can be categorized into deterministic and probabilistic forecasts. Deterministic forecasts provide a specific prediction, representing one possible outcome from a range of possibilities. They typically do not consider inherent uncertainties and often have shorter lead times [8,9,10]. Despite being derived from optimal initial conditions, deterministic forecasts may underestimate or overestimate actual outcomes. To address these limitations, the integration of ensemble precipitation forecasts (EPFs) into flood forecasting has gained traction in hydrological and meteorological research [11,12,13]. EPFs are generated by running multiple simulations from numerical weather prediction (NWP) models with slightly varied initial conditions or different model physics. By incorporating EPFs from NWPs into hydrological models, ensemble hydrological forecasts can be produced, offering a range of possible future states rather than a single deterministic prediction. Studies have shown that EPFs can significantly improve the accuracy and reliability of flood forecasts by accounting for the uncertainty in precipitation inputs [14,15,16,17].
Uncertainties in ensemble hydrological forecasts primarily originate from input data and the hydrological simulation process [18]. The uncertainty in hydrological simulations stems mainly from boundary conditions and the structural and parametric characteristics of the employed models [19,20]. This paper focuses on uncertainties arising from these two sources. To quantify uncertainties in input data, utilizing multiple EPFs from various weather forecasting centers—a “grand ensemble” (GE)—has proven effective [21]. A GE assimilates a wider array of uncertainty sources in the inputs, creating a more robust forecast probability density function (PDF) [22]. Conversely, a single EPF may address uncertainties specific to NWP models, potentially overlooking extreme events [23]. Since 2006, the THORPEX Interactive Grand Global Ensemble (TIGGE), part of the THORPEX scientific program launched by the World Meteorological Organization, has gathered ensemble forecasts from over 10 global NWP centers [24]. The study by Pappenberger, Bartholmes [23] demonstrates that forecasts from a GE capture flood peaks more effectively than those based on a single EPF, providing more reliable predictions.
A single hydrological model, however, does not encompass all potential output errors nor does it capture every aspect of the flood hydrological process [25]. Therefore, employing multiple hydrological models with varied structures helps quantify model structure uncertainty while enhancing forecasting capability and accuracy [26,27]. Recently, ensemble forecasting integrating multiple hydrological models has gained prominence [28,29,30]. Addressing uncertainties in both input data and hydrological model structures can be effectively managed by combining various EPF and multiple hydrological models (“EPF + HM”). This approach harnesses the strengths of different models while mitigating their limitations, thus enhancing the precision and predictability of hydrological forecasts [31,32,33].
Additionally, systematic bias and insufficient dispersion in meteorological ensemble forecasts, coupled with low spatial and temporal resolution when integrated with hydrological models and flawed spatio-temporal correlation structures between variables, necessitate enhancements through statistical post-processing techniques [34,35]. Over the past two decades, a variety of statistical post-processing methods have been developed, including linear bias correction [36], rank histogram techniques [37], reliability diagram techniques [38], Bayesian model averaging [39], kernel density estimation [40], logistic regression [41], non-homogeneous Gaussian regression [42], and so on. Research has demonstrated that applying these post-processing techniques to raw meteorological ensemble forecasts significantly improves the accuracy and reliability of hydrological predictions [43,44,45].
In light of these advances, this study hypothesizes that integrating statistical post-processing methods with GE and multiple hydrological models (MHM) will significantly enhance flood forecasting accuracy and reliability. To test this hypothesis, we propose the following research objectives. (1) Evaluate the effectiveness of the statistical post-processing method: assess how post-processing improves EPF performance in flood forecasting by correcting biases and addressing insufficient dispersion. (2) Investigate the utility of a GE: determine if a GE, which aggregates EPFs from multiple numerical weather prediction centers, provides more accurate flood peak predictions than a single EPF. (3) Assess the potential of MHM: evaluate whether combining different hydrological models enhances flood forecasting reliability compared to a single hydrological model (SHM), particularly in capturing event variability and uncertainty. (4) Examine the benefits of integrating GE with MHM: explore whether the combined use of GE and MHM improves the overall quality of flood forecasting, leveraging the strengths of both approaches for more precise and dependable predictions. The primary objective of this research is to identify the most effective combination of EPF and hydrological models (EPF+HM) for improving flood forecasting. The novelty of this study lies in developing and implementing a comprehensive framework that integrates GE and MHM for ensemble flood forecasting (EFF). By applying this innovative framework to the Qingjiang River Basin in the mountainous central-western region of China—characterized by high population density and vulnerability to frequent rainstorms and floods—we aim to address challenges such as systematic biases, insufficient dispersion, and spatio-temporal discrepancies in traditional models. This research represents the first attempt to apply multi-model techniques for ensemble flood forecasting in this region. The outcomes are expected to contribute to developing more effective flood forecasts, particularly in regions with similar topographical and climatic challenges.

2. Study Area and Data

The Qingjiang River Basin (108°35′–111°35′ E, 22°33′–30°50′ N), situated in the southwestern part of Hubei Province, China, flows from west to east through ten counties and cities and joins the Yangtze River at Yidu. The main channel stretches approximately 423 km, encompassing a controlled drainage area of about 16,700 square km. The basin experiences a mid-subtropical monsoon climate, with an average annual temperature of about 14 °C and an annual precipitation of about 1400 mm, predominantly concentrated in the summer months. The terrain of the basin, primarily consisting of mountainous regions that account for over 80% of its area, slopes from west to east, forming a landscape of high mountains and deep valleys. This region, located in western Hubei, is known for its severe rainstorms. Influenced by the summer rainstorms, summer runoff accounts for about 40% of the annual runoff, with peak floods typically occurring in June and July. Given the significant impact of human activities, particularly dam construction, on the middle and lower reaches of the basin, this study focuses on the upper Qingjiang River Basin. This area is less affected by human interventions, with the Shuibuya Hydropower Station serving as the upstream catchment area (Figure 1).
This study incorporates meteorological and hydrological observation data, as well as EPFs from NWP centers. The observation dataset includes daily temperature, precipitation, potential evaporation, and runoff data spanning from 2014 to 2020. Meteorological data were sourced from 175 automatic meteorological stations distributed across the study area, and hydrological data were collected from the Shuibuya Hydropower Station (Figure 1). The EPFs utilized in this research were obtained from two highly regarded NWP centers within the TIGGE database, specifically ECMWF and NCEP. Detailed information about these models is provided in Table 1.

3. Methodology

This study employs two EPFs and two hydrological models to conduct EFF, aiming to effectively manage the uncertainties in the flood forecasting process. These uncertainties originate from the input data and are propagated through various hydrological models, leading to uncertain outputs. Moreover, a statistical post-processing method is applied to correct biases and inadequate dispersion within raw EPFs. Finally, multiple evaluation metrics are utilized to assess the ensemble hydrological forecasts, encompassing both the raw and post-processed forecasts.

3.1. Post-Processing Method

In this study, the generator-based post-processing (GPP) method [36] was employed to correct EPFs. The GPP method is specifically designed for statistical post-process meteorological ensemble forecasts and has proven effective in previous research [45,46]. This method integrates meteorological forecast information into a weather generator, allowing for the regeneration of ensemble forecast outcomes that are consistent with both the ensemble meteorological forecasts and historical observational characteristics. The statistical distribution parameters for forecast variables within the weather generator are estimated from historical observational data. Firstly, the probability distribution function (PDF) for precipitation parameters is defined as follows:
y | x 1 , , x M ~ g ( y | x 1 , , x k )
where y represents the precipitation variable, x 1 , , x k corresponds to an ensemble of precipitation forecasts composed of K members. The PDF characterizes a mixed discrete/continuous distribution where the probability of zero precipitation is positive, and the probability of positive precipitation follows a continuous skewed distribution. The precipitation mixed distribution model proposed by [47] is as follows:
g ( y | f k ) = P ( y = 0 | f k ) I ( y = 0 ) + P ( y > 0 | f k ) g k ( y | f k ) I ( y > 0 )
where g ( y | f k ) denotes the probability distribution given the member forecasts; I() is the indicator function, which equals 1 if the expression inside the parenthesis holds true, and, otherwise, it equals 0; P ( y = 0 | f k ) and g ( y > 0 | f k ) represent the probabilities of no precipitation and precipitation, respectively; and g k ( y | f k ) signifies a two-parameter gamma distribution. When applying the GPP method to correct precipitation, the probability density functions for precipitation in different seasons or of different intensities are independently calibrated to match observational data. Subsequently, using the forecast information from the original EPF, samples are redrawn from these calibrated probability density functions to generate the post-processed EPFs. For a detailed methodology, refer to Chen, Brissette [36].

3.2. Hydrological Models

Hydrological models are utilized for various applications. Previous studies comparing lumped and distributed models in outlet streamflow simulation have found minimal differences in performance between the two model types [48,49]. Given the scope of this study and the prevalent use of lumped models for hydrological simulation and forecasting, lumped models were selected. Distributed models were excluded due to their complexity and data demands. Two lumped hydrological models were used in this study: the Xinanjiang model (XAJ) [50] and the mod`ele du G’enie Rural `a 4 param`etres (GR4J) [51]. These two models are widely used for streamflow simulation and have proven effective in various hydrological contexts [52,53,54].
The XAJ proposed by [50] is a lumped hydrological model designed for humid and semi-humid regions. This model conceptualizes the catchment as a soil box, consisting of upper, lower, and deep soil layers. The XAJ model’s computational process consists of four modules (Figure 2): (1) Calculation of catchment evapotranspiration and soil moisture, using a three-layer evapotranspiration model—rainfall first replenishes the upper layer, then the lower and deep layers; evapotranspiration initially depletes the upper layer’s soil moisture, followed sequentially by the lower and deep layers. (2) Runoff generation, which utilizes a saturation excess runoff model where a storage capacity curve controls the distribution of precipitation throughout the catchment to calculate total runoff. (3) Water source partitioning, using a free water storage capacity curve to separate the total runoff into surface runoff, subsurface flow, and underground runoff. (4) Catchment routing, where surface runoff typically uses the unit hydrograph method, while interflow and groundwater flow use the linear reservoir method. The XAJ model is characterized by 15 parameters: 4 related to evaporation, 2 for runoff generation, and 9 for runoff routing. These parameters (Table 2) require calibration using observed streamflow to accurately reflect the hydrological characteristics of a specific watershed. Further details can be found in the study by Zhao [50].
Developed in 2003 by French scholar Perrin et al., based on the GR3J model [55], the GR4J model has a simple structure with only four free parameters (Table 2) and uses two nonlinear reservoirs for production and routing calculations (Figure 3): (1) The production reservoir first determines effective precipitation Pn and residual evaporation capacity En based on catchment precipitation P and evaporation capacity E, then calculates the precipitation Ps and replenishes the production reservoir and the evapotranspiration Es. (2) The routing reservoir employs time segment unit hydrographs for routing calculations. Different runoff components have different routing times; thus, rainfall is split into two parts: 90% of the rainfall is routed using unit hydrograph UH1, needing adjustment through the routing reservoir, and the remaining 10% flows directly to the catchment outlet section, routed using unit hydrograph UH2. More details are available from Perrin, Michel [51].
The basic inputs for both models include catchment-average precipitation, temperature, and potential evapotranspiration, with outputs being the runoff at the catchment outlet. Catchment-average precipitation is calculated using the Thiessen Polygon method. The model calibration period spans from 2014 to 2017, with a validation period that spans from 2018 to 2020. During calibration, the Nash–Sutcliffe Efficiency coefficient (NSE) and the volumetric error (RE) are used as objective functions, with optimization performed using the SCE-UA algorithm. The formulas for NSE and RE are as follows:
N S E = 1 i = 1 n ( Q i o b s Q i s i m ) i = 1 n ( Q i o b s Q i o b s ¯ )
R E ( % ) = i = 1 n ( Q i o b s Q i s i m ) i = 1 n ( Q i o b s ) × 100
where Q i o b s represents the observed daily runoff, Q i s i m represents the simulated daily runoff, and Q ¯ o b s denotes the mean of the observed runoff over multiple days. The NSE value ranges from −∞ to 1, with 1 being the optimal value. The RE value ranges from 0 to ∞, with 0 being the optimal value.

3.3. Generation and Evaluation of Ensemble Flood Forecasting

The EFF is generated based on the following “EPF + HM” schemes: (A) ECMWF-driven XAJ, denoted as ECMWF + XAJ; (B) EPFs (ECMWF and NCEP) driving XAJ, referred to as GE + XAJ; (C) ECMWF driving multiple models (XAJ and GR4J), denoted as ECMWF + MHM; and (D) EPFs driving multiple hydrological models (XAJ and GR4J), known as GE + MHM. Scheme (D) represents the multi-ensemble, multi-model combination [27]. These four schemes generate ensemble members numbering 50, 70, 100, and 140, respectively. Additionally, ensemble forecast products corrected by the GPP method are used for flood ensemble forecasting in these four schemes. The results are evaluated and compared before and after post-processing.
This paper evaluates the performance of EFF using the ensemble mean and the 90% confidence interval (the difference between the 95th and 5th percentiles). The ensemble mean reflects the central tendency, while the 90% confidence interval quantifies the range, variability, and uncertainty within the ensemble. This dual assessment allows for a comprehensive evaluation of the ensemble’s accuracy in predicting flood events and in assessing the associated uncertainty. Both deterministic and probabilistic indicators are used. The deterministic indicators include the correlation coefficient (R) and the Relative Mean Error (RME), which measure the average deviation between the forecast mean and actual observations. Probabilistic indicators include the Continuous Ranked Probability Skill Score (CRPSS) and quantile histograms. To evaluate EFF uncertainty, three indicators are employed: Average Relative Interval Length (ARIL), the percentage of observations enclosed by the confidence interval (PCI), as well as the percentage of observations enclosed by the Unit Confidence Interval (PUCI) [56].
(1) R measures the correlation between forecasted and observed values, serving as an indicator of forecast reliability. It is calculated as follows:
R = i = 1 n ( O i O ¯ ) ( F i F ¯ ) i = 1 n ( O i O ¯ ) 2 i = 1 n ( F i F ¯ ) 2
where Oi and Fi represent observed and ensemble forecasted precipitation, respectively, while O ¯ and F ¯ are the daily mean values of observed and forecasted precipitation over an n-day time series. R values range from [0, 1], with values closer to 1 indicating more reliable forecasts.
(2) RME reflects the relative deviation between forecast values and the observed mean, assessing the accuracy of the forecasts. It is expressed as follows:
R M E = i = 1 n ( F i O ¯ ) i = 1 n Q i × 100 %
RME values range from [−100%, 100%], with values closer to 0 indicating smaller errors between the ensemble mean and observations.
(3) CRPSS evaluates the performance of the quantitative precipitation forecast (QPF) ensemble distribution relative to the climatological distribution. It is defined for a variable n at time t, estimated from the ensemble forecast probability density function as follows:
C R P S n t = + F n , t ( y ) I y X 2 d y
The continuous skill score, CRPSS, is defined as follows:
C R P S S n = 1 C R P S n C R P S n , r e f
where C R P S n , r e f is the CRPS obtained using a reference forecast. Higher CRPSS values approaching 1 indicate better ensemble forecast performance.
(4) Quantile histograms are used to check if the ensemble forecast values and observed values come from the same probability distribution. A flat, horizontal line shape indicates moderate dispersion within the ensemble; a concave shape with high sides and a low middle suggests under-dispersion; a convex shape with low sides and a high middle indicates over-dispersion. A flat histogram shape is a necessary condition for reliable ensemble forecasts. Additionally, the uniformity index of the quantile histogram is expressed as follows:
Δ = i = 1 m + 1 F i 1 m + 1
as Δ approaches 0, the flatness of the histogram shape improves, indicating more reliable ensemble forecasts.
(5) ARIL assesses the resolution of the predictive distributions, defined over the entire duration of the simulation as follows:
A R I L = 1 n i = 1 n Q u p p e r , i s i m Q l o w e r , i s i m Q i o b s
where n is the number of days in the observed record, Q u p p e r , i s i m and Q l o w e r , i s i m are the upper and lower simulated discharges of the 90% confidence interval, respectively, and Q i o b s is the observed discharge.
PCI quantifies the proportion of observed data falling within a specified range of simulated intervals, defined as follows:
P C I = 1 n i = 1 n Q i n , i o b s
where Q i n , i o b s represents the number of observed discharges within the 90% confidence interval. A smaller ARIL indicates a narrower uncertainty interval, while a higher PCI suggests a greater reliability of the interval.
PUCI, derived from ARIL and PCI, is defined as follows:
P U C I = ( 1 A b s ( P C I 0.9 ) ) A R I L
Higher PUCI values generally indicate lower uncertainty in the 90% confidence interval of discharge.

4. Results

4.1. Calibration and Validation of Two Hydrological Models

Based on daily meteorological and hydrological observations within the study area, the periods from 2014 to 2017 were designated as the calibration period, and from 2018 to 2020 as the validation period. During these periods, the XAJ and GR4J hydrological models were calibrated and validated. Figure 4 displays the observed and simulated hydrological process graphs for both XAJ and GR4J hydrological models during the calibration and validation periods. The graphs indicate differences in the models’ abilities to capture flood peaks, which can be attributed to their distinct structural characteristics. However, both models effectively captured flood peaks during the calibration period.
Moreover, during the calibration and validation periods, the NSE values for the XAJ model are 0.86 and 0.88, respectively, and the RE values are 8.58% and 7.26%, respectively. For the GR4J model, the NSE values are 0.92 and 0.91, and the RE values are −0.15% and 8.26%, respectively. Compared to earlier studies [57,58], both the XAJ and GR4J models demonstrate good applicability, as indicated by NSE and RE values, which suggests that the models are well calibrated and reliable for further flood forecasting research. The calibrated parameters for both hydrological models are detailed in Table 2.

4.2. Evaluation of Ensemble Flood Forecasting Performance for Four Schemes

Figure 5 displays the R between the ensemble mean and observed runoff, generated by the four schemes, both before and after correction with the GPP method, over a forecast horizon of 1–7 days. As shown in the figure, the R values for the 1 lead day are relatively high, at approximately 0.7 before correction and about 0.8 after correction. As the forecast horizon extends, the R values gradually decrease. Scheme (C), which focuses solely on the uncertainty of the hydrological model, consistently shows relatively higher correlations both before and after correction. In contrast, schemes (B) and (D), which consider the uncertainty of the forecasting models, exhibit a rapid decline in R values after the third day of the forecast period. Overall, the R values in the right graph are higher than those in the left graph, indicating that runoff forecasts derived by the corrected EPF exhibit a higher degree of correlation with observed values. In conclusion, these results suggest that considering the uncertainty of hydrological models in runoff forecasting, coupled with effective post-processing methods, can enhance the performance of runoff forecasts.
Figure 6 illustrates the RME between the ensemble mean and observed runoff forecasted using the four schemes, both before and after correction with the GPP method. The plots indicate a general tendency of the ensemble mean to underestimate observed values, with schemes (A) and (C) showing a relatively lower degree of underestimation, about −40%, before correction. In contrast, schemes (B) and (D) exhibit a more significant underestimation of observed runoff, with RME ranging from −80% to −40%, and this underestimation intensifies as the forecast horizon extends. After correction, the underestimation in all schemes is markedly improved, reduced to within −20%, demonstrating that effective post-processing techniques can significantly enhance the accuracy of runoff forecasts.
Figure 7 presents the CRPSS of the forecasted runoff ensembles obtained from the four schemes, before and after correction with the GPP method. This evaluation measures the probabilistic skills of the ensemble forecasts. Using the classifications established by Harrigan, Prudhomme [59] and Bennett, Wang [60], the flood forecasting skill levels are defined as follows: very high if CRPSS is between 0.75 and 1; high if CRPSS is between 0.5 and 0.75; moderate if CRPSS is between 0.25 and 0.5; low if CRPSS is between 0 and 0.25; no skill if CRPSS equals 0; and negative skill if CRPSS is less than 0. Accordingly, it is evident that before correction, schemes (A) and (C) show moderate skills for 1–7 lead days. However, schemes (B) and (D) show low skills once the forecast period exceeds 3 days. After applying the GPP correction, the CRPSS for all schemes improves. All four schemes show high skills for the first two lead days and moderate skills for 3–7 lead days. What is more, scheme (A) exhibits the lowest CRPSS on 1 lead day, and the others show slightly higher CRPSS values.
Figure 8 shows the quantile histograms of the forecasted runoff ensembles obtained from the four schemes, both before and after correction with the GPP method. The ranks of schemes (A) through (D) are 50, 70, 100, and 140, respectively. Initially, the quantile histograms of the uncorrected EPF display a U-shape, with a large concentration of forecast values at the lowest and highest ranks, indicating a large concentration of forecast values at the lowest and highest ranks. This pattern indicates under-dispersion, reflecting an underestimation of forecast uncertainty. After applying GPP correction, the histograms become flatter, indicating a more uniform distribution. Specifically, the uniformity indices for the quantile histograms of schemes (A) through (D) decreased from 1.17 to 0.42, 1.11 to 0.26, 1.13 to 0.35, and 1.09 to 0.31, respectively. This demonstrates that the GPP method effectively mitigates the issue of under-dispersion in ensemble forecasts. Additionally, compared to scheme (A), the other schemes have slightly lower uniformity indices, which indicates that considering the uncertainty in both forecasting models and hydrological models can further enhance the probabilistic skills of ensemble forecasts, to a certain extent.
In addition, Table 3 presents three uncertainty metrics for evaluating 90% confidence intervals of discharge derived from different schemes. It can be determined that (1) generally, both ARIL and PCI values derived by raw EPFs are smaller than those derived by GPP post-processed EPFs, making it difficult to determine whether post-processing is useful or not; (2) considering the PUCI, those derived by GPP post-processed EPFs have larger values than those derived by raw EPFs, indicating that their results are more reliable and exhibit less uncertainty; and (3) the PUCI values of schemes (A) and (C) are generally larger than those of schemes (B) and (D), suggesting that former approaches yield more robust and reliable forecasts.

4.3. Application to a Typical Flood Event

On 17 July 2020, the Qingjiang River Basin experienced a significant torrential rain event, leading to a rapid rise in river levels. On 18 July, the flow rate at the Shuibuya station in the Qingjiang Basin reached a peak of 6767 m3/s, marking an extreme flood event. Figure 9 presents the EFFs for this flood event simulated by the four “EPF + HM” schemes on 1 lead day. Figure 9a,b shows the raw EPF simulation results and the post-processed EPF results using the GPP method, respectively. The schemes (A)~(D) generated ensemble member counts of 50, 70, 100, and 140, respectively. The figures display the ensemble mean and the 90% confidence interval to evaluate the EFFs’ performance in capturing the flood event.
From Figure 9a, it is evident that the EFFs generated using the raw EPF exhibit biases and under-dispersion. The ensemble mean and 90% confidence interval in scheme (A) underestimate the observed flood peak. Compared to scheme (A), scheme (B) shows no improvement in the ensemble mean and 90% confidence interval, while scheme (C) performs slightly better in capturing the observed peak event, indicating that GE does not necessarily enhance the accuracy of the EFFs. Schemes (C) and (D), compared to schemes (A) and (B), have a wider range of forecast ensembles, suggesting that incorporating hydrological uncertainty can improve the quality of the EFFs. Thus, including multiple hydrological models (MHM) results in a broader range of ensemble forecasts and reduces bias compared to using only single hydrological models (SHMs).
From Figure 9b, the flood forecasts conducted with the GPP-corrected EPF show improvements in the range of the forecast ensemble. The ensemble mean in scheme (A) is still underestimated, but the forecast range can better cover the observed flood. The performance of scheme (B) is similar to that of scheme (A), indicating that when using a single hydrological model, GE does not demonstrate an improvement over a single EPF. Schemes (C) and (D) outperform other schemes in terms of the ensemble mean, providing higher estimates for the flood peak. These results suggest that while post-processing of the EPF has been applied, using GE does not effectively enhance the performance of the EFFs, whereas using MHM better captures the flood event.

5. Discussion

In this study, EFFs were generated by integrating both raw and post-processed EPFs, i.e., ECMWF and NCEP, into calibrated hydrological models, i.e., XAJ and GR4J. The flood forecasts were produced from four ‘NWP + HM’ schemes, i.e., (A) ECMWF + XAJ, (B) GE + XAJ, (C) ECMWF + MHM, and (D) GE + MHM. The performance of both the raw and post-processed flood forecasts was evaluated using deterministic indicators such as R and RME, alongside probabilistic indicators like CRPSS and quantile histograms, as well as three uncertainty evaluation indicators—ARIL, PCI, and PUCI.
The analysis reveals that forecasts driven by raw EPFs yield satisfactory correlation and relative error, but are limited in probabilistic metrics such as CRPSS and quantile histograms. Results demonstrate that after GPP correction, all schemes show marked improvements in both deterministic and probabilistic terms, with high deterministic skills and moderate probabilistic skills over 1–7 lead days. This enhancement also mitigated the degradation of forecast quality over lead days, validating the effectiveness of the GPP method, consistent with the findings of Zhang, Chen [46]. These findings underscore the necessity of employing post-processing methods to enhance the reliability of flood forecasts.
A comparative analysis between schemes (A) and (B) shows that using a single hydrological model, flood forecasts based on GE do not show a significant improvement over those using a single EPF. Comparisons between schemes (A) and (C) reveal that when using a single EPF, forecasts based on MHM offer a broader ensemble range, a smaller ensemble mean error, and an increased accuracy of the ensemble mean. This improvement in flood forecast quality is credited to the integration of hydrological model uncertainty. Comparing scheme (A) with scheme (D) shows that although using multiple forecasting models and hydrological models yields forecasts with higher dispersion, they exhibit lower correlation, larger biases, and lower CRPSS. This suggests that while it is reasonable to use GE to quantify input uncertainty, compared to a single EPS, GE may not necessarily enhance flood forecasting performance, aligning with the findings of Teja, Manikanta [31]. Therefore, it is recommended to combine a single EPF with MHM (i.e., scheme (C)) to achieve reliable flood forecasting. Moreover, results indicate that compared to input uncertainty, the uncertainty in hydrological models can more effectively enhance the predictability of flood forecasts. Therefore, focusing on hydrological model uncertainty rather than input data uncertainty could decrease computational time and burden in flood forecasting.

6. Conclusions

This study assesses the efficacy of EFFs in accurately predicting streamflow and flood events using calibrated hydrological models driven by TIGGE ensemble forecasts. Two lumped models and two EPFs were employed to address uncertainties in hydrological model structures and input uncertainties, respectively. Four forecasting schemes combining different hydrological models and EPFs were analyzed to determine the most effective setup for delivering high-quality flood forecasts. The GPP method was used to correct bias and under-dispersion within raw EPFs. The performance of EFFs was assessed using both deterministic and probabilistic metrics.
The results underscore the significance of applying post-processing techniques like GPP to correct EPFs to improve the reliability of flood forecasting. The application of the GPP has proven effective in maintaining reliable forecast skills for 1–7 lead days and reducing quality degradation over extended periods.
A critical finding from this study is the limited improvement in flood forecasts when using a GE with a single hydrological model, which did not substantially outperform forecasts derived from a single EPF. However, the integration of MHM with a single EPF, as seen in scheme (C), markedly enhanced forecast accuracy by offering a wider range of ensemble forecasts, lower bias, and improved correlation. This suggests that uncertainty in hydrological models plays a more critical role in enhancing forecast predictability than input data uncertainty alone. Thus, addressing hydrological model uncertainty is crucial for optimizing flood forecasting systems.
Nonetheless, this study has limitations, as the GPP post-processing method corrects only temporal errors and does not consider spatial corrections. Scholars suggest that addressing errors from the spatial distribution of precipitation could further improve runoff forecast accuracy [61,62]. Moreover, given the complex terrain of the study area and lack of observational data, this study only used two lumped hydrological models. However, semi-distributed or distributed models, which better handle spatial variability in rainfall, might provide better hydrological forecasting outcomes [63]. Exploring how different hydrological models perform under varied geographic and climatic conditions in future research could provide valuable insights.
In conclusion, this research advocates for the strategic combination of a single EPF with multiple hydrological models to achieve more reliable flood forecasting. This approach not only harnesses the strengths of post-processing techniques but also aligns with the practical needs of reducing computational demands in flood forecasting practices.

Author Contributions

Conceptualization, Y.X.; Methodology, Y.X.; Software, T.S.; Validation, Z.Y.; Formal Analysis, H.Q.; Investigation, Y.X. and H.Q.; Resources, T.P.; Data Curation, Z.Y. and T.S.; Writing—Original Draft, Y.X.; Writing—Review and Editing, Y.X. and T.P.; Project Administration, T.P.; Funding Acquisition, T.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work is jointly supported by the Yangtze River Water Science Joint Research Fund (U2340201), the Natural Science Foundation of Hubei Province (Joint Fund) (2023AFD094 and 2022CFD129), the Open Grants of the State Key Laboratory of Severe Weather (2024LASW-B04), the Project of Yangtze River Basin Meteorological Opening Fund (CJLY2022Y06), and the Basic Research Fund of WHIHR (WHIHRKYYW2024).

Data Availability Statement

The meteorological observations are publicly archived and can be accessed through the China Meteorological Administration. Streamflow observations can be obtained from the Shuibuya Hydropower Station. Public datasets from ECMWF and NCEP utilized in this study are available at https://apps.ecmwf.int/datasets/data/tigge/levtype=sfc/type=cf/ (accessed on 28 November 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ICHARM Report. Global Trends in Water Related Disasters: An Insight for Policymakers; International Centre for Water Hazard and Risk Management (UNESCO): Tsukuba, Japan, 2009; Available online: http://www.icharm.pwri.go.jp (accessed on 25 March 2024).
  2. Alfieri, L.; Bisselink, B.; Dottori, F.; Naumann, G.; de Roo, A.; Salamon, P.; Wyser, K.; Feyen, L. Global projections of river flood risk in a warmer world. Earth’s Future 2017, 5, 171–182. [Google Scholar] [CrossRef]
  3. WMO. Manual on Flood Forecasting and Warning; WMO No. 1072; World Meteorological Organization: Geneva, Switzerland, 2011. [Google Scholar]
  4. Difrancesco, K.; Tullos, D. Flexibility in Water Resources Management: Review of Concepts and Development of Assessment Measures for Flood Management Systems. JAWRA J. Am. Water Resour. Assoc. 2014, 50, 1527–1539. [Google Scholar] [CrossRef]
  5. Hopson, T.M.; Webster, P.J. A 1–10-day ensemble forecasting scheme for the major river basins of Bangladesh: Forecasting severe floods of 2003–2007. J. Hydrometeorol. 2010, 11, 618–641. [Google Scholar] [CrossRef]
  6. Alfieri, L.; Burek, P.; Dutra, E.; Krzeminski, B.; Muraro, D.; Thielen, J.; Pappenberger, F. GloFAS–global ensemble streamflow forecasting and flood early warning. Hydrol. Earth Syst. Sci. 2013, 17, 1161–1175. [Google Scholar] [CrossRef]
  7. Thiemig, V.; Bisselink, B.; Pappenberger, F.; Thielen, J. A pan-African medium-range ensemble flood forecast system. Hydrol. Earth Syst. Sci. 2015, 19, 3365–3385. [Google Scholar] [CrossRef]
  8. Zhao, P.; Wang, Q.J.; Wu, W.; Yang, Q. Which precipitation forecasts to use? Deterministic versus coarser-resolution ensemble NWP models. Q. J. R. Meteorol. Soc. 2021, 147, 900–913. [Google Scholar] [CrossRef]
  9. Toth, Z.; Talagrand, O.; Candille, G.; Zhu, Y. Probability and ensemble forecasts. Forecast Verif. Pract. Guide Atmos. Sci. 2003, 137, 163. [Google Scholar]
  10. Boucher, M.A.; Anctil, F.; Perreault, L.; Tremblay, D. A comparison between ensemble and deterministic hydrological forecasts in an operational context. Adv. Geosci. 2011, 29, 85–94. [Google Scholar] [CrossRef]
  11. Yucel, I.; Onen, A.; Yilmaz, K.; Gochis, D. Calibration and evaluation of a flood forecasting system: Utility of numerical weather prediction model, data assimilation and satellite-based rainfall. J. Hydrol. 2015, 523, 49–66. [Google Scholar] [CrossRef]
  12. Rodwell, M.; Palmer, T. Using numerical weather prediction to assess climate models. Q. J. R. Meteorol. Soc. J. Atmos. Sci. Appl. Meteorol. Phys. Oceanogr. 2007, 133, 129–146. [Google Scholar] [CrossRef]
  13. Bauer, P.; Thorpe, A.; Brunet, G. The quiet revolution of numerical weather prediction. Nature 2015, 525, 47–55. [Google Scholar] [CrossRef] [PubMed]
  14. Ming, X.; Liang, Q.; Xia, X.; Li, D.; Fowler, H.J. Real-time flood forecasting based on a high-performance 2-D hydrodynamic model and numerical weather predictions. Water Resour. Res. 2020, 56, e2019WR025583. [Google Scholar] [CrossRef]
  15. Xiang, Y.; Peng, T.; Gao, Q.; Shen, T.; Qi, H. Evaluation of TIGGE precipitation forecast and its applicability in streamflow predictions over a Mountain River Basin, China. Water 2022, 14, 2432. [Google Scholar] [CrossRef]
  16. Alfieri, L.; Pappenberger, F.; Wetterhall, F.; Haiden, T.; Richardson, D.; Salamon, P. Evaluation of ensemble streamflow predictions in Europe. J. Hydrol. 2014, 517, 913–922. [Google Scholar] [CrossRef]
  17. Jiang, X.; Zhang, L.; Liang, Z.; Fu, X.; Wang, J.; Xu, J.; Zhang, Y.; Zhong, Q. Study of early flood warning based on postprocessed predicted precipitation and Xinanjiang model. Weather Clim. Extrem. 2023, 42, 100611. [Google Scholar] [CrossRef]
  18. Krzysztofowicz, R. Bayesian theory of probabilistic forecasting via deterministic hydrologic model. Water Resour. Res. 1999, 35, 2739–2750. [Google Scholar] [CrossRef]
  19. Pappenberger, F.; Matgen, P.; Beven, K.J.; Henry, J.-B.; Pfister, L. Influence of uncertain boundary conditions and model structure on flood inundation predictions. Adv. Water Resour. 2006, 29, 1430–1449. [Google Scholar] [CrossRef]
  20. Ajami, N.K.; Duan, Q.; Sorooshian, S. An integrated hydrologic Bayesian multimodel combination framework: Confronting input, parameter, and model structural uncertainty in hydrologic prediction. Water Resour. Res. 2007, 45, 208–214. [Google Scholar] [CrossRef]
  21. Qu, B.; Zhang, X.; Pappenberger, F.; Zhang, T.; Fang, Y. Multi-model grand ensemble hydrologic forecasting in the Fu river basin using Bayesian model averaging. Water 2017, 9, 74. [Google Scholar] [CrossRef]
  22. Cuo, L.; Pagano, T.C.; Wang, Q. A review of quantitative precipitation forecasts and their use in short-to medium-range streamflow forecasting. J. Hydrometeorol. 2011, 12, 713–728. [Google Scholar] [CrossRef]
  23. Pappenberger, F.; Bartholmes, J.; Thielen, J.; Cloke, H.L.; Buizza, R.; de Roo, A. New dimensions in early flood warning across the globe using grand-ensemble weather predictions. Geophys. Res. Lett. 2008, 35, 1956–1964. [Google Scholar] [CrossRef]
  24. Swinbank, R.; Kyouda, M.; Buchanan, P.; Froude, L.; Hamill, T.M.; Hewson, T.D.; Keller, J.H.; Matsueda, M.; Methven, J.; Pappenberger, F. The TIGGE project and its achievements. Bull. Am. Meteorol. Soc. 2016, 97, 49–67. [Google Scholar] [CrossRef]
  25. Ajami, N.K.; Duan, Q.; Gao, X.; Sorooshian, S. Multimodel combination techniques for analysis of hydrological simulations: Application to distributed model intercomparison project results. J. Hydrometeorol. 2006, 7, 755–768. [Google Scholar] [CrossRef]
  26. Dion, P.; Martel, J.-L.; Arsenault, R. Hydrological ensemble forecasting using a multi-model framework. J. Hydrol. 2021, 600, 126537. [Google Scholar] [CrossRef]
  27. Troin, M.; Arsenault, R.; Wood, A.W.; Brissette, F.; Martel, J.-L. Generating Ensemble Streamflow Forecasts: A Review of Methods and Approaches over the Past 40 Years. Water Resour. Res. 2021, 57, e2020WR028392. [Google Scholar] [CrossRef]
  28. Li, Z.; Yu, J.; Xu, X.; Sun, W.; Pang, B.; Yue, J. Multi-model ensemble hydrological simulation using a BP Neural Network for the upper Yalongjiang River Basin, China. Proc. Int. Assoc. Hydrol. Sci. 2018, 379, 335–341. [Google Scholar] [CrossRef]
  29. Sharma, S.; Siddique, R.; Reed, S.; Ahnert, P.; Mejia, A. Hydrological model diversity enhances streamflow forecast skill at short-to medium-range timescales. Water Resour. Res. 2019, 55, 1510–1530. [Google Scholar] [CrossRef]
  30. Xu, J.; Anctil, F.; Boucher, M.-A. Hydrological post-processing of streamflow forecasts issued from multimodel ensemble prediction systems. J. Hydrol. 2019, 578, 124002. [Google Scholar] [CrossRef]
  31. Teja, K.N.; Manikanta, V.; Das, J.; Umamahesh, N. Enhancing the predictability of flood forecasts by combining Numerical Weather Prediction ensembles with multiple hydrological models. J. Hydrol. 2023, 625, 130176. [Google Scholar] [CrossRef]
  32. Wanders, N.; Thober, S.; Kumar, R.; Pan, M.; Sheffield, J.; Samaniego, L.; Wood, E.F. Development and evaluation of a pan-European multimodel seasonal hydrological forecasting system. J. Hydrometeorol. 2019, 20, 99–115. [Google Scholar] [CrossRef]
  33. Han, S.; Coulibaly, P. Probabilistic Flood Forecasting Using Hydrologic Uncertainty Processor with Ensemble Weather Forecasts. J. Hydrometeorol. 2019, 20, 1379–1398. [Google Scholar] [CrossRef]
  34. Scheuerer, M.; Hamill, T.M. Statistical postprocessing of ensemble precipitation forecasts by fitting censored, shifted gamma distributions. Mon. Weather Rev. 2015, 143, 4578–4596. [Google Scholar] [CrossRef]
  35. Wilks, D.S. Comparison of ensemble-MOS methods in the Lorenz’96 setting. Meteorol. Appl. 2006, 13, 243–256. [Google Scholar] [CrossRef]
  36. Chen, J.; Brissette, F.P.; Li, Z. Postprocessing of ensemble weather forecasts using a stochastic weather generator. Mon. Weather Rev. 2014, 142, 1106–1124. [Google Scholar] [CrossRef]
  37. Hamill, T.M.; Colucci, S.J. Verification of Eta–RSM short-range ensemble forecasts. Mon. Weather Rev. 1997, 125, 1312–1327. [Google Scholar] [CrossRef]
  38. Atger, F. Spatial and interannual variability of the reliability of ensemble-based probabilistic forecasts: Consequences for calibration. Mon. Weather Rev. 2003, 131, 1509–1523. [Google Scholar] [CrossRef]
  39. Raftery, A.E.; Gneiting, T.; Balabdaoui, F.; Polakowski, M. Using Bayesian model averaging to calibrate forecast ensembles. Mon. Weather Rev. 2005, 133, 1155–1174. [Google Scholar] [CrossRef]
  40. Roulston, M.S.; Smith, L.A. Combining dynamical and statistical ensembles. Tellus A Dyn. Meteorol. Oceanogr. 2003, 55, 16–30. [Google Scholar] [CrossRef]
  41. Wilks, D.S. Extending logistic regression to provide full-probability-distribution MOS forecasts. Meteorol. Appl. J. Forecast. Pract. Appl. Train. Tech. Model. 2009, 16, 361–368. [Google Scholar] [CrossRef]
  42. Gneiting, T.; Raftery, A.E.; Westveld, A.H.; Goldman, T. Calibrated probabilistic forecasting using ensemble model output statistics and minimum CRPS estimation. Mon. Weather Rev. 2005, 133, 1098–1118. [Google Scholar] [CrossRef]
  43. Cloke, H.L.; Pappenberger, F.; van Andel, S.J.; Schaake, J.; Thielen, J.; Ramos, M.-H. Hydrological Ensemble Prediction Systems; Wiley-Blackwell: Chichester, UK, 2013. [Google Scholar]
  44. Javanshiri, Z.; Fathi, M.; Mohammadi, S.A. Comparison of the BMA and EMOS statistical methods for probabilistic quantitative precipitation forecasting. Meteorol. Appl. 2021, 28, e1974. [Google Scholar] [CrossRef]
  45. Xiang, Y.; Liu, Y.; Zou, X.; Peng, T.; Yin, Z.; Ren, Y. Post-Processing Ensemble Precipitation Forecasts and Their Applications in Summer Streamflow Prediction over a Mountain River Basin. Atmosphere 2023, 14, 1645. [Google Scholar] [CrossRef]
  46. Zhang, J.; Chen, J.; Li, X.; Chen, H.; Xie, P.; Li, W. Combining postprocessed ensemble weather forecasts and multiple hydrological models for ensemble streamflow predictions. J. Hydrol. Eng. 2020, 25, 04019060. [Google Scholar] [CrossRef]
  47. Sloughter, J.M.L.; Raftery, A.E.; Gneiting, T.; Fraley, C. Probabilistic quantitative precipitation forecasting using Bayesian model averaging. Mon. Weather Rev. 2007, 135, 3209–3220. [Google Scholar] [CrossRef]
  48. Kumar, A.; Singh, R.; Jena, P.P.; Chatterjee, C.; Mishra, A. Identification of the best multi-model combination for simulating river discharge. J. Hydrol. 2015, 525, 313–325. [Google Scholar] [CrossRef]
  49. Vansteenkiste, T.; Tavakoli, M.; Ntegeka, V.; De Smedt, F.; Batelaan, O.; Pereira, F.; Willems, P. Intercomparison of hydrological model structures and calibration approaches in climate scenario impact projections. J. Hydrol. 2014, 519, 743–755. [Google Scholar] [CrossRef]
  50. Zhao, R.-J. The Xinanjiang model applied in China. J. Hydrol. 1992, 135, 371–381. [Google Scholar]
  51. Perrin, C.; Michel, C.; Andréassian, V. Improvement of a parsimonious model for streamflow simulation. J. Hydrol. 2003, 279, 275–289. [Google Scholar] [CrossRef]
  52. Arsenault, R.; Gatien, P.; Renaud, B.; Brissette, F.; Martel, J.-L. A comparative analysis of 9 multi-model averaging approaches in hydrological continuous streamflow simulation. J. Hydrol. 2015, 529, 754–767. [Google Scholar] [CrossRef]
  53. Mathevet, T.; Gupta, H.; Perrin, C.; Andréassian, V.; Le Moine, N. Assessing the performance and robustness of two conceptual rainfall-runoff models on a worldwide sample of watersheds. J. Hydrol. 2020, 585, 124698. [Google Scholar] [CrossRef]
  54. Xiang, Y.; Chen, J.; Li, L.; Peng, T.; Yin, Z. Evaluation of eight global precipitation datasets in hydrological modeling. Remote Sens. 2021, 13, 2831. [Google Scholar] [CrossRef]
  55. Nascimento, N.D.O.; Yang, X.; Makhlouf, Z.; Michel, C. GR3J: A daily watershed model with three free parameters. Hydrol. Sci. J. 1999, 44, 263–277. [Google Scholar]
  56. Li, L.; Xu, C.Y.; Xia, J.; Engeland, K.; Reggiani, P. Uncertainty estimates by Bayesian method with likelihood of AR (1) plus Normal model and AR (1) plus Multi-Normal model in different time-scales hydrological models. J. Hydrol. 2011, 406, 54–65. [Google Scholar] [CrossRef]
  57. Li, X.-Q.; Chen, J.; Xu, C.-Y.; Li, L.; Chen, H. Performance of post-processed methods in hydrological predictions evaluated by deterministic and probabilistic criteria. Water Resour. Manag. 2019, 33, 3289–3302. [Google Scholar] [CrossRef]
  58. Shu, Z.; Zhang, J.; Jin, J.; Wang, L.; Wang, G.; Wang, J.; Sun, Z.; Liu, J.; Liu, Y.; He, R. Evaluation and application of quantitative precipitation forecast products for mainland China based on TIGGE multimodel data. J. Hydrometeorol. 2021, 22, 1199–1219. [Google Scholar] [CrossRef]
  59. Harrigan, S.; Prudhomme, C.; Parry, S.; Smith, K.; Tanguy, M. Benchmarking ensemble streamflow prediction skill in the UK. Hydrol. Earth Syst. Sci. 2018, 22, 2023–2039. [Google Scholar] [CrossRef]
  60. Bennett, J.C.; Wang, Q.J.; Robertson, D.E.; Schepen, A.; Li, M.; Michael, K. Assessment of an ensemble seasonal streamflow forecasting system for Australia. Hydrol. Earth Syst. Sci. 2017, 21, 6007–6030. [Google Scholar] [CrossRef]
  61. Chen, J.; Li, X.; Xu, C.-Y.; Zhang, X.J.; Xiong, L.; Guo, Q. Postprocessing Ensemble Weather Forecasts for Introducing Multisite and Multivariable Correlations Using Rank Shuffle and Copula Theory. Mon. Weather Rev. 2022, 150, 551–565. [Google Scholar] [CrossRef]
  62. Wang, L.; Bao, H. Ensemble flood forecasting based on ensemble NWP and the GMKHM distributed hydrological model. MATEC Web Conf. 2018, 246, 01108. [Google Scholar] [CrossRef]
  63. Das, T.; Bárdossy, A.; Zehe, E.; He, Y. Comparison of conceptual model performance using different representations of spatial variability. J. Hydrol. 2008, 356, 106–118. [Google Scholar] [CrossRef]
Figure 1. The Qingjiang River Basin.
Figure 1. The Qingjiang River Basin.
Water 16 01887 g001
Figure 2. Structure of the XAJ model (Zhao, 1992 [50]).
Figure 2. Structure of the XAJ model (Zhao, 1992 [50]).
Water 16 01887 g002
Figure 3. Structure of the GR4J model (Perrin, 2003 [51]).
Figure 3. Structure of the GR4J model (Perrin, 2003 [51]).
Water 16 01887 g003
Figure 4. Calibration and validation hydrographs generated by the XAJ and GR4J hydrological models, along with observed streamflow at the Shuibuya station.
Figure 4. Calibration and validation hydrographs generated by the XAJ and GR4J hydrological models, along with observed streamflow at the Shuibuya station.
Water 16 01887 g004
Figure 5. Correlation (R) plots for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 4 forecast schemes for 1–7 lead days.
Figure 5. Correlation (R) plots for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 4 forecast schemes for 1–7 lead days.
Water 16 01887 g005
Figure 6. Relative Mean Error (RME) plots for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 4 forecast schemes for 1–7 lead days.
Figure 6. Relative Mean Error (RME) plots for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 4 forecast schemes for 1–7 lead days.
Water 16 01887 g006
Figure 7. CRPSS plots for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 4 forecast schemes for 1–7 lead days.
Figure 7. CRPSS plots for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 4 forecast schemes for 1–7 lead days.
Water 16 01887 g007
Figure 8. Rank Histograms for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 1 lead day.
Figure 8. Rank Histograms for flood forecasts driven by RAW and GPP post-processed precipitation forecasts for 1 lead day.
Water 16 01887 g008
Figure 9. Hydrograph of ensemble forecasts for a flooding event driven by RAW (left) and GPP post-processed (right) precipitation forecasts at 1 lead day.
Figure 9. Hydrograph of ensemble forecasts for a flooding event driven by RAW (left) and GPP post-processed (right) precipitation forecasts at 1 lead day.
Water 16 01887 g009
Table 1. Overview of the datasets used in this study.
Table 1. Overview of the datasets used in this study.
DatasetsSourcesEnsemble Members
(Perturbed)
Temporal ResolutionSpatial
Resolution
Temporal Span
Observed meteorological data175 meteorological stations-DailyStation2014–2020
Ensemble precipitation forecasts (EPFs)ECMWF (European Centre for Medium-Range Weather Forecasts)507 lead days
at 6 h
0.5 grid2014–2020
NCEP (National Centers for Environmental Prediction)20
Observed
discharge data
Shuibuya hydropower station-DailyStation 2014–2020
Table 2. Calibrated parameters in both XAJ and GR4J models.
Table 2. Calibrated parameters in both XAJ and GR4J models.
ParametersDescriptionRangeCalibrated Value
XAJWMAreal tension water capacity (mm)Humid regions: 100~160134.16
WUMFraction of upper water in WM0~0.60.26
WLMFraction of lower water in WM0.4~10.65
KEvaporation coefficient0~0.60.55
CCoefficient of evapotranspiration in the lower soil layer0.2~10.62
IMImpermeable coeflicient0~0.20.07
BTension water distribution index0~0.020.01
EXFree water distribution index1~1.51.39
SMAreal free water capacity (mm)0~6059.98
KIFraction of free water to interflow0~10.50
KGFraction of free water to groundwater0~1, KI + KG < 10.50
CIInterflow recession coefficient0.8~10.80
CGGroundwater recession coefficient0.8~10.90
NNumber of reservoirs of the Nash model0~50.27
KStorage constant of the Nash model0~53.04
GR4JX1Production reservoir: storage of rainfall on the soil surface (mm)1~750187.07
X2Groundwater exchange coefficient: a
function of groundwater exchange
that influences the routing reservoir
−10~102.91
X3Routing storage: amount of water that can be stored in soil porous (mm)1~40068.84
X4Time peak: the time when the ordinate
peak of flood hydrograph is created
0.5~101.33
Table 3. The Average Relative Interval Length (ARIL), the percentage of observations enclosed by the confidence interval (PCI), and the percentage of observations enclosed by the Unit Confidence Interval (PUCI) values for the 90% confidence interval of discharge for different schemes.
Table 3. The Average Relative Interval Length (ARIL), the percentage of observations enclosed by the confidence interval (PCI), and the percentage of observations enclosed by the Unit Confidence Interval (PUCI) values for the 90% confidence interval of discharge for different schemes.
SchemesRawGPP
ARILPCIPUCIARILPCIPUCI
(A) CMWF + XAJ1.3040.3310.3302.6150.7640.331
(B) GE + XAJ1.8770.3600.2452.9350.8250.315
(C) ECMWF + MHM1.7650.4500.3122.6710.8160.343
(D) GE + MHM2.3170.4700.2462.9800.8630.323
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xiang, Y.; Peng, T.; Qi, H.; Yin, Z.; Shen, T. Improving Flood Forecasting Skill by Combining Ensemble Precipitation Forecasts and Multiple Hydrological Models in a Mountainous Basin. Water 2024, 16, 1887. https://doi.org/10.3390/w16131887

AMA Style

Xiang Y, Peng T, Qi H, Yin Z, Shen T. Improving Flood Forecasting Skill by Combining Ensemble Precipitation Forecasts and Multiple Hydrological Models in a Mountainous Basin. Water. 2024; 16(13):1887. https://doi.org/10.3390/w16131887

Chicago/Turabian Style

Xiang, Yiheng, Tao Peng, Haixia Qi, Zhiyuan Yin, and Tieyuan Shen. 2024. "Improving Flood Forecasting Skill by Combining Ensemble Precipitation Forecasts and Multiple Hydrological Models in a Mountainous Basin" Water 16, no. 13: 1887. https://doi.org/10.3390/w16131887

APA Style

Xiang, Y., Peng, T., Qi, H., Yin, Z., & Shen, T. (2024). Improving Flood Forecasting Skill by Combining Ensemble Precipitation Forecasts and Multiple Hydrological Models in a Mountainous Basin. Water, 16(13), 1887. https://doi.org/10.3390/w16131887

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop