A TimeVarying Distributed Unit Hydrograph considering soil moisture content
 ^{1}School of Civil and Hydraulic Engineering, Huazhong University of Science and Technology, Wuhan, 430074, China
 ^{2}Hubei Key Laboratory of Digital Valley Science and Technology, Wuhan 430074, China
 ^{3}PowerChina Huadong Engineering Corporation Limited, Hangzhou 310014, China
 ^{4}Meizhou Hydrological Bureau, Guangdong Province, Meizhou 514000, China
 ^{1}School of Civil and Hydraulic Engineering, Huazhong University of Science and Technology, Wuhan, 430074, China
 ^{2}Hubei Key Laboratory of Digital Valley Science and Technology, Wuhan 430074, China
 ^{3}PowerChina Huadong Engineering Corporation Limited, Hangzhou 310014, China
 ^{4}Meizhou Hydrological Bureau, Guangdong Province, Meizhou 514000, China
Abstract. The distributed unit hydrograph (DUH) method has been widely used for flood routing simulation, because it can well characterize the underlying surface characteristics and various rainfall intensities. The core of the DUH is the calculation of flow velocity. However, the current velocity formula assumed a global equilibrium of the watershed and ignored the impact of timevarying soil moisture content on flow velocity, which leads to a larger flow velocity value. The goal of this study is to identify a soil moisture content factor, which was derived based on the water storage capacity curve, to explore the responses of DUH to soil moisture content in unsaturated areas. Thus, an improved distributed unit hydrograph based on timevarying soil moisture content was proposed in this paper. The proposed method considered the impact of both the timevarying rainfall intensity and soil moisture content on the flow velocity, and the watershed is assumed not to be equilibrium but vary with the soil moisture. The Qin River Basin was selected as a case study, and results of the timevarying distributed unit hydrograph (TDUH) and current DUH methods were used as comparisons with that of proposed method. Influence mechanism of timevarying soil moisture content on the flow velocity and flood forecasts were explored. Results show that the proposed method performs the best among the three methods. The shape and duration of the unit hydrograph can be mainly related to the soil moisture content at initial stage of a storm. When the watershed is approximately saturated, the grid flow velocity is majorly dominated by the excess rainfall.
Bin Yi et al.
Status: closed

RC1: 'Comment on hess2021470', Anonymous Referee #1, 21 Oct 2021
The manuscript proposes a Unit Hydrograph estimating travel times including also the timevarying rainfall intensity and soil moisture content information.
The topic is surely interesting although is well crystallized in literature and in practical hydrology.
Unfortunately, the manuscript has some drawbacks that do not allow me to suggest its publication. I see practical, technical, and theoretical issues to be addressed.
Firstly, the text (language, structure, figure captions, typos) should be significantly improved since presently, it does not allow a full understanding of described methods.
It is not clear (although the title is clear) if the manuscript proposes a IUH model, that is a rainfallrunoff method, or a simplified routing approach. Reading the title and the introduction it seems the first option, looking the case study the second one.
As recalled by the authors in the Introduction the IUH approach assumes some hypotheses (i.e. linear system, time invariant, rainfall spatially homogeneous) that of course are far from the watersheds real behavior, however this is an accepted compromise in challenging hydrological studies (i.e. ungauged basins). Â In my opinion there is a contradiction in trying to make nonlinear a linear approach, maybe it is better to use a non linear approach or an other approach. If a timevarying rainfall and spatial distribution information are introduced, the nature of IUH is lost and we do not know any more what we are implementing. In my experience I also tried to force the IUH concept, however I limited it to the estimation of the hillslope velocity cell by cell indeed further adaptation (i.e. empirical estimation of channel velocity) would have been in contrast to the IUH theoretical definition (Grimaldi et al. 2010; 2012).
Most importantly, it should be clarified the practical context on which we are referring. Personally, I consider the WFIUH approach particularly brilliant in small and ungauged basin application since it optimizes the topographic information and since the IUH approximations are coherent with by the basin dimension. In other contexts, maybe, other approaches should be preferred (semidistributed or fully distributed models).
The present manuscript does not clarify these aspects. Â It provides a case study with a large basin, eliminating the hillslopes (including an area threshold of 1 km2), assuming calibration, increasing the number of parameters: without a clear context it disorients the reader. I would expect to see nine case studies (each subbasin) in order to evaluate the expected improvements of IUH given by the soil moisture content, in ungauged contest. If the practical aim is different (large gauged basins) the comparison should be performed with other models underlying the added value of the proposed approach.
Grimaldi, S., Petroselli, A., Alonso, G., Nardi, F. Flow time estimation with spatially variable hillslope velocity in ungauged basins (2010) Advances in Water Resources, 33 (10), pp. 12161223.
Grimaldi, S., Petroselli, A., Nardi, F. A parsimonious geomorphological unit hydrograph for rainfallrunoff modelling in small ungauged basins (2012) Hydrological Sciences Journal, 57 (1), pp. 7383.Â
 AC1: 'Reply on RC1', Bin Yi, 10 Nov 2021

RC2: 'Comment on hess2021470', Jasper Vrugt, 26 Oct 2021
Review of "A TimeVarying Distributed Unit Hydrograph considering soil moisture content"Â
Summary: In this paper the authors propose an extension of the distributed unit hydrograph (DUH) method for routing of water in distributed hydrologic models. The default DUB approach takes into consideration the topography of the land surface in computation of surface runoff and river inflow but does not account for spatial variability of the soil moisture capacity within the catchment. The authors present a simple extension of the flow velocity equation which is thought to account for a varying moisture content.
Evaluation: I am not an expert on routing methods, in fact, never taken a course in surface hydrology, but I do know something about math, hydrology, and modeling. The authors presented in this paper are interesting, yet, I believe the paper warrants a major revision before it can be judged to making a significant contribution to hydrology and water resources. Specifically, the paper is not very well written, syntax and grammar need major improvements. I believe this is a first requirement before the paper is ready for detailed review. What is more, the methodology needs a much better physical underpinning citing properly past publications on this topic, presenting relevant units of variables, addressing sensitivity of results to variables such as gamma. Furthermore, the case study / demonstration of the methodology is not particularly convincing. I highlight my main comments below  not in a particular order of importance.Â
1. The authors should reference previous work on routing/modeling approaches. For example, Line 269, Eq. 10 does not provide a reference, whereas this equation is simply the Pareto distribution function, used by Moore (1985) to describe the spatial variability of the soil moisture storage capacity in the watershed. This is just one example  this comment applies to many equations used by the authors; the continuity equation, Manning's equation, etc.
2. The derivation of the TDUH method needs significant improvement. Not well written  the derivation has been given in other manuscripts so either present it briefly and clearly, otherwise, I recommend the authors refer to the original sources as the present derivation brings up more questions than it answers. I am particularly bothered with the way the equations and variables are presented. It reads as a collection of some equations with some variables.Â Also, I am left wondering whether the equations presented are derived for the first time by the authors or whether they have been presented in the literature years ago? For example, Eqs. (6), (7), (8) and (9).Â Â Â
3. Authors should give units of variables they use. This will make it easier for readers to digest the material, and for students to reproduce/implement the approach the authors have presented.
4. Line 208209: Variable "m" is introduced, but is not used in Equation (1) or (2). Introduce variables when they are used and not ahead of time, unless this makes sense to do.
5. Line 261  302: The extension the authors propose, essentially applies the ideas of Moore (1985) to the unit hydrograph. Authors should do a much better job connecting what they do to the literature.
6. Equation (15)  is this equation a simple extension of Eq. (14) with "w_{t}" raised to gamma ? or is this equation from Bhattacharya et al. (2012) and/or Bunster et al. (2019)?Â
7. The authors state that their routing method takes account of soil moisture content, but soil moisture content does not appear in any of the equations. Instead, they use the Pareto distribution function to express the spatial variability of the soil moisture capacity. They assume that this capcity represents temporal variations in soil moisture content. Furthermore, the representative volume of the soil moisture content is unknown. Are we considering the topsoil moisture content, or the moisture content of the first 50 or 100 cm of the profile? I guess I am looking for a better physical underpinning of the presented method.Â
8. Line 261  268: I read this paragraph several times, and it is still not clear to me. What is B_{t}? Is "A_{t}" the same A as used in the continuity equation? Why not use symbol theta for moisture content? Why do I need to compute the ratio of A and A+B?Â
9. Why does WM' use the prime symbol?Â
10. Line 261  268: We have A for moisture content, B for the maximum soil moisture storage, W for current soil moisture storage. Why using so many different variables for essentially the same thing? Also, what is the unit of storage and moisture content? Are they similar, or different (as they should be). Is soil moisture storage not simply equal to soil moisture content x depth of profile. Why not use theta for soil moisture content and S for soil moisture storage, where S = L*theta, where L is the depth of the profile?
11.Â Line 269  289: I have a hard time following the different steps of the methodology. I think the authors unnecessarily confuse readers  the methodology can be presented in much easier to understand language  and in a much more coherent style.Â
12. Equation (10) presents how we compute alpha, but then in a next equation, alpha is a function of time. Please make clear in your entire derivation which variables are constant (scalars), which ones vary as function of time/space (scalars) and, if necessary, which are vectors and/or matrices. WM' varies as function of time? Otherwise alpha is constant.Â
13. The exponent b of the Pareto distribution function. Is this constant, or varies throughout your watershed?Â
14. Line 283: Mention that w_t varies between 0 and 1, thus, $w_{t} \in (0,1]$ (in latex). I assume that w_t cannot be zero as soil can never be entirely be depleted from water.Â
15. Equation 13: The denominator may need further explanation. Either solve analytically the integral of Eq. (12) and substitute this in Eq. 13, or do the analytic integration explicitly in Eq. 13.Â
16. The denominator of Eq. 13  last step  is wrongly formulated. Do we do 1  (b/(b+1))*1  alpha_t^{1/b} or Â 1  (b/(b+1)) * ( 1  alpha_t^{1/b} ); If the first then remove the 1 to get 1  (b/(b+1))  alpha_t^{1/b}, etc. The present formulation is unclear.
17. Line 296: The variable 'gamma', is this constant for the entire watershed, or varies per subcatchment or grid cell or ? I recommend that for each parameter the text explains how this parameter is treated, besides its units, a description of what the parameter represents, etc.
18. The authors use an aggregated objective function. Why use a single aggregate objective function? I would advise analyzing the performance metrics seperately. The SCE_UA method can do three separate trials  each using a different objective. Then you can compare the results of the proposed method against existing unit hydrograph routing methods proposed in the literature. If so desired, you can even consider the aggregate objective function  but then as fourth option.Â
19. What is the definition of flood peak? Need more information to compute the first objective function; are you looking at the peaks of the measured discharge record? and then use the exact same indices of the simulated record to compute the objective function? Or are you getting the indices of the peaks from the simulated record? and then use the corresponding measured values to compute the OF? This difference in implementation may seem insignificant, but can lead to widely different results.Â Â
20. Why define the peak error as bias? Why is this preferred over a standard squared residual metric? L2norm versus L1 norm. Same question for the timing error  and see above comment as well for this 2nd metric.Â
21. Why not use a metric such as the sum of squared residuals to compare the measured and simulated discharge records of the different routing methods?Â Â
22. I doubt that the improvement of the new routing method proposed by the authors is related to incorporating what the authors believe to be soil moisture content. What is key to the performance of the new routing method is what is done to the parameter gamma. The authors articulate what they have done with gamma on Lines 365  373. The value of gamma will determine the results of the new method; hence, is the performance improvement related to the gamma parameter, simply as this provides additional flexibility to routing?Â
23. The different routing methods amount to a model selection problem  and proper techniques such as information criteria or the marginal likelihood should be used to compare the different routing methods. I strongly doubt that what is shown in this paper is the result of soil moisture. The methodology and statistical analysis should be much improved to inspire confidence in this conclusion. As it stands right now, there are many other reasons so as to why the new method outperforms the existing routing methods. One of which is the parameter gamma. If nothing else, the analysis should show the sensitivity of the new routing method to the choice of gamma. The same should be done for parameter k in the existing formulation.Â
24. Some section names confuse the reader; for example, section 4.2 is labeled "Derivation of TDUH considering timevarying soil moisture content", but the derivation has already been presented. In 4.2, the authors simply apply their method to a case study. Unless "derivation" has another meaning and refers to the computation of the TDUH.Â
I leave it with this for now as further review will reiterate similar points. I very much appreciate the efforts of the authors in trying to improve the description of the unit hydrograph for distributed hydrologic modeling. Yet, the present paper needs a major revisition before it can be judged to making a significant contribution to the field and warrant publication in HESS/HESSD. As it stands right now, the methodology needs to be much better embedded into the literature and cite relevant papers when using existing equations, etc. and improve considerably the physical underpinning of the presented method. The authors should also revisit their calibration method  and provide compelling information about the sensitivity of their findings to the choice of gamma and k. Furthermore, the presentation and writing need considerable improvement.Â
I hope my comments are useful to improve the presentation and description of the methodology, including its application in a case study.
 AC2: 'Reply on RC2', Bin Yi, 03 Dec 2021
Status: closed

RC1: 'Comment on hess2021470', Anonymous Referee #1, 21 Oct 2021
The manuscript proposes a Unit Hydrograph estimating travel times including also the timevarying rainfall intensity and soil moisture content information.
The topic is surely interesting although is well crystallized in literature and in practical hydrology.
Unfortunately, the manuscript has some drawbacks that do not allow me to suggest its publication. I see practical, technical, and theoretical issues to be addressed.
Firstly, the text (language, structure, figure captions, typos) should be significantly improved since presently, it does not allow a full understanding of described methods.
It is not clear (although the title is clear) if the manuscript proposes a IUH model, that is a rainfallrunoff method, or a simplified routing approach. Reading the title and the introduction it seems the first option, looking the case study the second one.
As recalled by the authors in the Introduction the IUH approach assumes some hypotheses (i.e. linear system, time invariant, rainfall spatially homogeneous) that of course are far from the watersheds real behavior, however this is an accepted compromise in challenging hydrological studies (i.e. ungauged basins). Â In my opinion there is a contradiction in trying to make nonlinear a linear approach, maybe it is better to use a non linear approach or an other approach. If a timevarying rainfall and spatial distribution information are introduced, the nature of IUH is lost and we do not know any more what we are implementing. In my experience I also tried to force the IUH concept, however I limited it to the estimation of the hillslope velocity cell by cell indeed further adaptation (i.e. empirical estimation of channel velocity) would have been in contrast to the IUH theoretical definition (Grimaldi et al. 2010; 2012).
Most importantly, it should be clarified the practical context on which we are referring. Personally, I consider the WFIUH approach particularly brilliant in small and ungauged basin application since it optimizes the topographic information and since the IUH approximations are coherent with by the basin dimension. In other contexts, maybe, other approaches should be preferred (semidistributed or fully distributed models).
The present manuscript does not clarify these aspects. Â It provides a case study with a large basin, eliminating the hillslopes (including an area threshold of 1 km2), assuming calibration, increasing the number of parameters: without a clear context it disorients the reader. I would expect to see nine case studies (each subbasin) in order to evaluate the expected improvements of IUH given by the soil moisture content, in ungauged contest. If the practical aim is different (large gauged basins) the comparison should be performed with other models underlying the added value of the proposed approach.
Grimaldi, S., Petroselli, A., Alonso, G., Nardi, F. Flow time estimation with spatially variable hillslope velocity in ungauged basins (2010) Advances in Water Resources, 33 (10), pp. 12161223.
Grimaldi, S., Petroselli, A., Nardi, F. A parsimonious geomorphological unit hydrograph for rainfallrunoff modelling in small ungauged basins (2012) Hydrological Sciences Journal, 57 (1), pp. 7383.Â
 AC1: 'Reply on RC1', Bin Yi, 10 Nov 2021

RC2: 'Comment on hess2021470', Jasper Vrugt, 26 Oct 2021
Review of "A TimeVarying Distributed Unit Hydrograph considering soil moisture content"Â
Summary: In this paper the authors propose an extension of the distributed unit hydrograph (DUH) method for routing of water in distributed hydrologic models. The default DUB approach takes into consideration the topography of the land surface in computation of surface runoff and river inflow but does not account for spatial variability of the soil moisture capacity within the catchment. The authors present a simple extension of the flow velocity equation which is thought to account for a varying moisture content.
Evaluation: I am not an expert on routing methods, in fact, never taken a course in surface hydrology, but I do know something about math, hydrology, and modeling. The authors presented in this paper are interesting, yet, I believe the paper warrants a major revision before it can be judged to making a significant contribution to hydrology and water resources. Specifically, the paper is not very well written, syntax and grammar need major improvements. I believe this is a first requirement before the paper is ready for detailed review. What is more, the methodology needs a much better physical underpinning citing properly past publications on this topic, presenting relevant units of variables, addressing sensitivity of results to variables such as gamma. Furthermore, the case study / demonstration of the methodology is not particularly convincing. I highlight my main comments below  not in a particular order of importance.Â
1. The authors should reference previous work on routing/modeling approaches. For example, Line 269, Eq. 10 does not provide a reference, whereas this equation is simply the Pareto distribution function, used by Moore (1985) to describe the spatial variability of the soil moisture storage capacity in the watershed. This is just one example  this comment applies to many equations used by the authors; the continuity equation, Manning's equation, etc.
2. The derivation of the TDUH method needs significant improvement. Not well written  the derivation has been given in other manuscripts so either present it briefly and clearly, otherwise, I recommend the authors refer to the original sources as the present derivation brings up more questions than it answers. I am particularly bothered with the way the equations and variables are presented. It reads as a collection of some equations with some variables.Â Also, I am left wondering whether the equations presented are derived for the first time by the authors or whether they have been presented in the literature years ago? For example, Eqs. (6), (7), (8) and (9).Â Â Â
3. Authors should give units of variables they use. This will make it easier for readers to digest the material, and for students to reproduce/implement the approach the authors have presented.
4. Line 208209: Variable "m" is introduced, but is not used in Equation (1) or (2). Introduce variables when they are used and not ahead of time, unless this makes sense to do.
5. Line 261  302: The extension the authors propose, essentially applies the ideas of Moore (1985) to the unit hydrograph. Authors should do a much better job connecting what they do to the literature.
6. Equation (15)  is this equation a simple extension of Eq. (14) with "w_{t}" raised to gamma ? or is this equation from Bhattacharya et al. (2012) and/or Bunster et al. (2019)?Â
7. The authors state that their routing method takes account of soil moisture content, but soil moisture content does not appear in any of the equations. Instead, they use the Pareto distribution function to express the spatial variability of the soil moisture capacity. They assume that this capcity represents temporal variations in soil moisture content. Furthermore, the representative volume of the soil moisture content is unknown. Are we considering the topsoil moisture content, or the moisture content of the first 50 or 100 cm of the profile? I guess I am looking for a better physical underpinning of the presented method.Â
8. Line 261  268: I read this paragraph several times, and it is still not clear to me. What is B_{t}? Is "A_{t}" the same A as used in the continuity equation? Why not use symbol theta for moisture content? Why do I need to compute the ratio of A and A+B?Â
9. Why does WM' use the prime symbol?Â
10. Line 261  268: We have A for moisture content, B for the maximum soil moisture storage, W for current soil moisture storage. Why using so many different variables for essentially the same thing? Also, what is the unit of storage and moisture content? Are they similar, or different (as they should be). Is soil moisture storage not simply equal to soil moisture content x depth of profile. Why not use theta for soil moisture content and S for soil moisture storage, where S = L*theta, where L is the depth of the profile?
11.Â Line 269  289: I have a hard time following the different steps of the methodology. I think the authors unnecessarily confuse readers  the methodology can be presented in much easier to understand language  and in a much more coherent style.Â
12. Equation (10) presents how we compute alpha, but then in a next equation, alpha is a function of time. Please make clear in your entire derivation which variables are constant (scalars), which ones vary as function of time/space (scalars) and, if necessary, which are vectors and/or matrices. WM' varies as function of time? Otherwise alpha is constant.Â
13. The exponent b of the Pareto distribution function. Is this constant, or varies throughout your watershed?Â
14. Line 283: Mention that w_t varies between 0 and 1, thus, $w_{t} \in (0,1]$ (in latex). I assume that w_t cannot be zero as soil can never be entirely be depleted from water.Â
15. Equation 13: The denominator may need further explanation. Either solve analytically the integral of Eq. (12) and substitute this in Eq. 13, or do the analytic integration explicitly in Eq. 13.Â
16. The denominator of Eq. 13  last step  is wrongly formulated. Do we do 1  (b/(b+1))*1  alpha_t^{1/b} or Â 1  (b/(b+1)) * ( 1  alpha_t^{1/b} ); If the first then remove the 1 to get 1  (b/(b+1))  alpha_t^{1/b}, etc. The present formulation is unclear.
17. Line 296: The variable 'gamma', is this constant for the entire watershed, or varies per subcatchment or grid cell or ? I recommend that for each parameter the text explains how this parameter is treated, besides its units, a description of what the parameter represents, etc.
18. The authors use an aggregated objective function. Why use a single aggregate objective function? I would advise analyzing the performance metrics seperately. The SCE_UA method can do three separate trials  each using a different objective. Then you can compare the results of the proposed method against existing unit hydrograph routing methods proposed in the literature. If so desired, you can even consider the aggregate objective function  but then as fourth option.Â
19. What is the definition of flood peak? Need more information to compute the first objective function; are you looking at the peaks of the measured discharge record? and then use the exact same indices of the simulated record to compute the objective function? Or are you getting the indices of the peaks from the simulated record? and then use the corresponding measured values to compute the OF? This difference in implementation may seem insignificant, but can lead to widely different results.Â Â
20. Why define the peak error as bias? Why is this preferred over a standard squared residual metric? L2norm versus L1 norm. Same question for the timing error  and see above comment as well for this 2nd metric.Â
21. Why not use a metric such as the sum of squared residuals to compare the measured and simulated discharge records of the different routing methods?Â Â
22. I doubt that the improvement of the new routing method proposed by the authors is related to incorporating what the authors believe to be soil moisture content. What is key to the performance of the new routing method is what is done to the parameter gamma. The authors articulate what they have done with gamma on Lines 365  373. The value of gamma will determine the results of the new method; hence, is the performance improvement related to the gamma parameter, simply as this provides additional flexibility to routing?Â
23. The different routing methods amount to a model selection problem  and proper techniques such as information criteria or the marginal likelihood should be used to compare the different routing methods. I strongly doubt that what is shown in this paper is the result of soil moisture. The methodology and statistical analysis should be much improved to inspire confidence in this conclusion. As it stands right now, there are many other reasons so as to why the new method outperforms the existing routing methods. One of which is the parameter gamma. If nothing else, the analysis should show the sensitivity of the new routing method to the choice of gamma. The same should be done for parameter k in the existing formulation.Â
24. Some section names confuse the reader; for example, section 4.2 is labeled "Derivation of TDUH considering timevarying soil moisture content", but the derivation has already been presented. In 4.2, the authors simply apply their method to a case study. Unless "derivation" has another meaning and refers to the computation of the TDUH.Â
I leave it with this for now as further review will reiterate similar points. I very much appreciate the efforts of the authors in trying to improve the description of the unit hydrograph for distributed hydrologic modeling. Yet, the present paper needs a major revisition before it can be judged to making a significant contribution to the field and warrant publication in HESS/HESSD. As it stands right now, the methodology needs to be much better embedded into the literature and cite relevant papers when using existing equations, etc. and improve considerably the physical underpinning of the presented method. The authors should also revisit their calibration method  and provide compelling information about the sensitivity of their findings to the choice of gamma and k. Furthermore, the presentation and writing need considerable improvement.Â
I hope my comments are useful to improve the presentation and description of the methodology, including its application in a case study.
 AC2: 'Reply on RC2', Bin Yi, 03 Dec 2021
Bin Yi et al.
Bin Yi et al.
Viewed
HTML  XML  Total  BibTeX  EndNote  

586  137  16  739  10  9 
 HTML: 586
 PDF: 137
 XML: 16
 Total: 739
 BibTeX: 10
 EndNote: 9
Viewed (geographical distribution)
Country  #  Views  % 

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