the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Regional index flood estimation at multiple durations with generalized additive models
Abstract. Estimation of flood quantiles at ungauged basins is often achieved through regression based methods. In situations where flood retention is important, e.g. floodplain management and reservoir design, flood quantile estimates are often needed at multiple durations. This poses a problem for regressionbased models as the form of the functional relationship between catchment descriptors and the response may not be constant across different durations. A particular type of regression model that is wellsuited to this situation is a generalized additive model (GAM), which allows for flexible, semiparametric modeling and visualization of the relationship between predictors and the response. However, in practice, selecting predictors for such a flexible model can be challenging, particularly given the characteristics of available catchment descriptor datasets. We employ a machine learningbased variable preselection tool which, when combined with domain knowledge, enhances the practicality of constructing GAMs. In this study, we develop a GAM for index (median) flood estimation with the primary objective of investigating durationspecific differences in how catchment descriptors influence the median flood. As the accuracy of this explainable approach is dependent on the fitted GAM being adequate, the secondary objective of our study is prediction of the median flood at ungauged locations and multiple durations, where predictive performance and reliability at ungauged locations are used as proxies for adequacy of the GAM. Predictive performance of the GAM is compared to two benchmark models: the existing loglinear model for median flood estimation in Norway and a fully datadriven machine learning model (an extreme gradient boosting tree ensemble, XGBoost). We find that the predictive accuracy and reliability of the GAM matched or exceeded that of the benchmark models at both durations studied. Within the predictor set selected for this study, we observe durationspecific differences in the relationship between the median flood and the two catchment descriptors effective lake percentage and catchment shape. Ignoring these differences results in a statistically significant decline in predictive performance. This suggests that models developed and estimated for prediction of the index flood at one duration may have reduced performance when applied directly to situations outside of that specific duration.
 Preprint
(2598 KB)  Metadata XML
 BibTeX
 EndNote
Status: closed

RC1: 'Comment on egusphere20232335', Anonymous Referee #1, 15 Nov 2023
The manuscript "Regional index flood estimation at multiple durations with generalized additive models" by Barna et al. is a comprehensive study on prediction of the median index flood on 234 stations in Norway. It aims to compare a GAM model with two benchmark models (a loglinear model and XGBoost) on two different flood duration (1h and 24h). Additionally, a variable selection method is included and models are tested in a crossvalidation approach. The manuscript is generally well written and the results are worth publishing in HESS. However, I have two concerns about the manuscript, which I think can be addressed in a revised version.
First, one of the objectives of the study is "prediction of the median flood at ungauged locations" (Line 125126). I am not quite sure if this is addressed adequately. One of the catchment descriptors is the mean annual runoff of the catchment of the SeNorge 2.0 dataset, which is based on observational data. In my opinion prediction at ungauged locations means, that there is no information about the runoff at this stations, so also no information about the mean annual runoff can be included in a prediction model. I understand that the information maybe necessary for the comparison with the RFA_2018 model, but I think it can also be beneficial to show that the GAM model performs as good without the mean annual runoff (and leave the mean annual runoff as predictor for the RFA_2018 model for simplicity). My second point concerning the prediction in ungauged locations, is the variable selection. If I understood it correctly, the variable selection is performed on the full dataset (with a cross validation scheme), and the selected variables are then used as input for the validation study (again with a crossvalidation scheme). If this is correct, the variables are selected on the full dataset, not on a subset, so the prediction error is somehow biased, as the variable selection already included information about the full dataset. If my understanding of the validation scheme is wrong, I would suggest to make this clearer in the methods section.
Second, the manuscript is a bit too long. Some parts of the methods are in the Appendix, which is fine, but makes it hard to read at some point. Two examples: (i) The description in the Introduction between line 2946 may be shortened, (ii) or the very detailed description of the reponse variable (Sect 2./2.1) can may be written more concise.
Additional minor comments:
The legend title of Figure 1c seems wrong.
Add panels in Figure 2.
In Table 1, the catchment descriptors $R_L$ is duplicated.
Line 351, what are the actual hyperparmameters that were used in the 10 final models? If it is easy possible, the information could be added.
Line 384: XGBoost is assessed only on the MAE; optimal predictors for the other four error metrics are not accessible for XGBoost when the data are assumed log normal. Why is XGBoost then used as for the variable selection procedure  and the MAPE as error metric, if this will produce unreliable result? Additionally, was it considered to alter the loss function in XGBoost for better comparison in terms of the chosen error metrics? I think this should be clarified in the main manuscript.
I think the formulas in Figure 5(a) may be wrong ($\times 100$), and as they are both partly defined in Table 3., I would suggest using the same variables.
Figure 6, the legend is not consistent with the actual plot.
Line 452  469, I think the description of partial response curves can be shortened.
Line 487, Maybe the results on regional models can be skipped in the Appendix.
Line 526, If this sentence is referring to the paragraph above, maybe change the order of the two sentences.
Citation: https://doi.org/10.5194/egusphere20232335RC1  AC2: 'Reply on RC1', Danielle Barna, 28 Feb 2024

RC2: 'Comment on egusphere20232335', Anonymous Referee #2, 13 Jan 2024
Report on the manuscript untitled:
Regional index flood estimation at multiple durations with generalized additive models
The topic of the manuscript is appealing and at first sight it appears interesting. However, there are high number of major issues based on which I cannot recommend publication of this manuscript (such as validation of the employed approach is missing, recent and related literature is ignored, confusing terminology). I’m providing in the following major comments only:
 The title is confusing and misleading in different aspects, such as:  different durations could be understood simultaneously or as a variable  the index flood model is based on quantile while the paper treats only the median (fixed quantile order)
 why not using directly and fully Machine Learning? The compatibility between a statistical model GAM and variable selection method based on ML should be discussed. Especially, later in the paper, there is a formal method to select variables for GAM (implemented in mgcv package). In the same idea, (line 43) I’m wondering if the authors are using a modified version of IIS. Hence, this method should be checked and validated before for this choice/context. The question is about the compatibility of this change.
 around lines 4045: this text is ambiguous and not well justified/ motivated. It is based on a unique old paper (see next comments for recent papers). This is part is crucial and motivates the study. Hence, the paper motivation and foundation are questionable.
 Indeed it is more informative to include the duration in the modeling. However, to deal with the duration, it is now appropriate to consider a multivariate framework involving the duration as a variable and simultaneously with other variable like the peak and/volume. The multivariate regional framework, index flood model, is already developed (e.g. Requena et al. 2016, J. of Hydrology; Azam et al. 2018, Water).
 why this and only these values (1h and 24h)?
 Around line 60: Dealing with nonlinearity is not only through transformation but directly using nonlinear approaches (see e.g. Ouali et al. 2017, J. Advances in Modeling Earth Systems; Cannon 2018, Stochastic environmental research and risk).
 line 108: not sure about this statement, especially no refs given. As far as I know, variable selection is not the strength of ML. I don’t know what is reported in Guisan et al. 2002, but it may be not up to date (given the fast development of ML).
 line 115: I'm surprised to see that such an important topic is treated only in the hydrological framework. It is questionable to heavily rely on this.
 Last line page 4: This assumption is either strong or in contradiction with the problematic to be treated in the paper.
 It is important to provide an equation for “median annual maximum flood” to be explicit and avoid confusion.
 Some parts of the methodology should be in the results section (section 4.2 and from line 250).
 Equation 6: something is missing or wrong. The righthand side does not depend on i (so the summation is over what?).
 Using the term permutation test could be misleading since this is a generic term on how to obtain pvalue.
Citation: https://doi.org/10.5194/egusphere20232335RC2  AC1: 'Reply on RC2', Danielle Barna, 28 Feb 2024
Status: closed

RC1: 'Comment on egusphere20232335', Anonymous Referee #1, 15 Nov 2023
The manuscript "Regional index flood estimation at multiple durations with generalized additive models" by Barna et al. is a comprehensive study on prediction of the median index flood on 234 stations in Norway. It aims to compare a GAM model with two benchmark models (a loglinear model and XGBoost) on two different flood duration (1h and 24h). Additionally, a variable selection method is included and models are tested in a crossvalidation approach. The manuscript is generally well written and the results are worth publishing in HESS. However, I have two concerns about the manuscript, which I think can be addressed in a revised version.
First, one of the objectives of the study is "prediction of the median flood at ungauged locations" (Line 125126). I am not quite sure if this is addressed adequately. One of the catchment descriptors is the mean annual runoff of the catchment of the SeNorge 2.0 dataset, which is based on observational data. In my opinion prediction at ungauged locations means, that there is no information about the runoff at this stations, so also no information about the mean annual runoff can be included in a prediction model. I understand that the information maybe necessary for the comparison with the RFA_2018 model, but I think it can also be beneficial to show that the GAM model performs as good without the mean annual runoff (and leave the mean annual runoff as predictor for the RFA_2018 model for simplicity). My second point concerning the prediction in ungauged locations, is the variable selection. If I understood it correctly, the variable selection is performed on the full dataset (with a cross validation scheme), and the selected variables are then used as input for the validation study (again with a crossvalidation scheme). If this is correct, the variables are selected on the full dataset, not on a subset, so the prediction error is somehow biased, as the variable selection already included information about the full dataset. If my understanding of the validation scheme is wrong, I would suggest to make this clearer in the methods section.
Second, the manuscript is a bit too long. Some parts of the methods are in the Appendix, which is fine, but makes it hard to read at some point. Two examples: (i) The description in the Introduction between line 2946 may be shortened, (ii) or the very detailed description of the reponse variable (Sect 2./2.1) can may be written more concise.
Additional minor comments:
The legend title of Figure 1c seems wrong.
Add panels in Figure 2.
In Table 1, the catchment descriptors $R_L$ is duplicated.
Line 351, what are the actual hyperparmameters that were used in the 10 final models? If it is easy possible, the information could be added.
Line 384: XGBoost is assessed only on the MAE; optimal predictors for the other four error metrics are not accessible for XGBoost when the data are assumed log normal. Why is XGBoost then used as for the variable selection procedure  and the MAPE as error metric, if this will produce unreliable result? Additionally, was it considered to alter the loss function in XGBoost for better comparison in terms of the chosen error metrics? I think this should be clarified in the main manuscript.
I think the formulas in Figure 5(a) may be wrong ($\times 100$), and as they are both partly defined in Table 3., I would suggest using the same variables.
Figure 6, the legend is not consistent with the actual plot.
Line 452  469, I think the description of partial response curves can be shortened.
Line 487, Maybe the results on regional models can be skipped in the Appendix.
Line 526, If this sentence is referring to the paragraph above, maybe change the order of the two sentences.
Citation: https://doi.org/10.5194/egusphere20232335RC1  AC2: 'Reply on RC1', Danielle Barna, 28 Feb 2024

RC2: 'Comment on egusphere20232335', Anonymous Referee #2, 13 Jan 2024
Report on the manuscript untitled:
Regional index flood estimation at multiple durations with generalized additive models
The topic of the manuscript is appealing and at first sight it appears interesting. However, there are high number of major issues based on which I cannot recommend publication of this manuscript (such as validation of the employed approach is missing, recent and related literature is ignored, confusing terminology). I’m providing in the following major comments only:
 The title is confusing and misleading in different aspects, such as:  different durations could be understood simultaneously or as a variable  the index flood model is based on quantile while the paper treats only the median (fixed quantile order)
 why not using directly and fully Machine Learning? The compatibility between a statistical model GAM and variable selection method based on ML should be discussed. Especially, later in the paper, there is a formal method to select variables for GAM (implemented in mgcv package). In the same idea, (line 43) I’m wondering if the authors are using a modified version of IIS. Hence, this method should be checked and validated before for this choice/context. The question is about the compatibility of this change.
 around lines 4045: this text is ambiguous and not well justified/ motivated. It is based on a unique old paper (see next comments for recent papers). This is part is crucial and motivates the study. Hence, the paper motivation and foundation are questionable.
 Indeed it is more informative to include the duration in the modeling. However, to deal with the duration, it is now appropriate to consider a multivariate framework involving the duration as a variable and simultaneously with other variable like the peak and/volume. The multivariate regional framework, index flood model, is already developed (e.g. Requena et al. 2016, J. of Hydrology; Azam et al. 2018, Water).
 why this and only these values (1h and 24h)?
 Around line 60: Dealing with nonlinearity is not only through transformation but directly using nonlinear approaches (see e.g. Ouali et al. 2017, J. Advances in Modeling Earth Systems; Cannon 2018, Stochastic environmental research and risk).
 line 108: not sure about this statement, especially no refs given. As far as I know, variable selection is not the strength of ML. I don’t know what is reported in Guisan et al. 2002, but it may be not up to date (given the fast development of ML).
 line 115: I'm surprised to see that such an important topic is treated only in the hydrological framework. It is questionable to heavily rely on this.
 Last line page 4: This assumption is either strong or in contradiction with the problematic to be treated in the paper.
 It is important to provide an equation for “median annual maximum flood” to be explicit and avoid confusion.
 Some parts of the methodology should be in the results section (section 4.2 and from line 250).
 Equation 6: something is missing or wrong. The righthand side does not depend on i (so the summation is over what?).
 Using the term permutation test could be misleading since this is a generic term on how to obtain pvalue.
Citation: https://doi.org/10.5194/egusphere20232335RC2  AC1: 'Reply on RC2', Danielle Barna, 28 Feb 2024
Viewed
HTML  XML  Total  BibTeX  EndNote  

1,228  117  31  1,376  25  24 
 HTML: 1,228
 PDF: 117
 XML: 31
 Total: 1,376
 BibTeX: 25
 EndNote: 24
Viewed (geographical distribution)
Country  #  Views  % 

Total:  0 
HTML:  0 
PDF:  0 
XML:  0 
 1