the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Technical note: Analyzing river network dynamics and the active length–discharge relationship using water presence sensors
Francesca Zanetti
Nicola Durighetto
Filippo Vingiani
Gianluca Botter
Download
- Final revised paper (published on 08 Jul 2022)
- Preprint (discussion started on 29 Mar 2021)
Interactive discussion
Status: closed
-
RC1: 'Comment on hess-2021-103', Anonymous Referee #1, 12 May 2021
Review of the manuscript HESS-2021-103 “Analysing river network dynamics and active length – discharge relationship using water presence sensors” by Zanetti et al.
General comment
This is an interesting work based on a relatively new technology that has the potential to offer new insight into the dynamics of stream network in small catchments. Streams are highly dynamic systems and characterizing such “liveliness” is important to move towards a better understanding of catchments’ functioning.
Overall, this manuscript is well written, logically structured, and clearly illustrated. The introduction is solid and the results are well supported by the data. I have only some specific and minor comments that I would like the Authors to address (see below). Overall, I recommend a moderate revision but since the MODERATE grade is not available to HESS reviewers, I indicated minor review in the review form.
Specific comments
L78-83. I think that the three research questions could be more appealing than the current ones. I understand that, as this is a relatively new measurement approach, aiming at providing the reader with some methodological observation is useful for possible further application of the methods. So, I like question 1. However, I would move it as last question because, in my opinion, it’s more important to focus on processes that the method is able to describe and understand, rather than on the method itself. The second question is a bit too narrow because it implies a “yes” or “nor” and does not lead to much insight into hydrological processes. Similarly, question 3 sounds a bit too “methodological” and not so oriented towards process understanding. So, I suggest moving question 1 as third, and to rephrase question 2 and 3 and bit. Of course, the results and discussions should be reorganized to reflect these changes. Indeed, my overall impression is that this work is much technically- and methodologically-oriented and less prone to describe and understand hydrological processes, and I guess that reader of this journal are more interested in knowing how catchment works rather than know if some data can fit a certain model or no. So, I invite the Authors to consider revising the work to reflect this aspect.
L127. What are the criteria for the choice of the field deployment? Please, specify.
L145. This sounds a very short period to me. I understand that practical issues might be arisen but typically we need a longer time period to observe hydrologically processes that often highly variable in time. Si the average rainfall and streamflow in this period comparable to long-term rainfall and streamflow (or compared to the other years where observations are available, since this is aa relatively new experimental catchment). I think it’s important for the Authors to discuss this issue, explain why (if) they reckon this is a suitable spell and why, and why (if) this is a representative period for derive information on the hydrological functioning of this catchment. Moreover, the should discuss how this short period of time can potentially impact on the results.
L184. Where was the sensor placed? In the grass in a convergent zone, I guess, where water was not flowing? Please, give more information on the aspect of sensors deployment in the field.
Minor comments
L5. I suggest considering the term “customized” instead than “personalized”.
L6. I suggest removing “analysed,”.
L6. The expression “nodes’ persistency” is not clear without reading the manuscript. Please, clarify.
L37. Perhaps here the citation to Godsay and Kirchner (2014) fits well. In any case, a more recent references would be a nice complement.
L129. “heterogeneous persistencies so as to avoid redundancy in the data.” This is not clear to me, please specify.
L160. Which observations? Please specify. What does “some modelling” mean? Please, explain.
L176. Building a reliable flow rating curve, especially in mountain catchments, is a challenge. How many points were collected to build the FRC, and what was the range of streamflow values and the resulting goodness of fit measures? In other terms, is the FRC reliable to infer robust streamflow values? Please, explain.
L178-179. I understand that the Authors want to give light to the ERC project but this is not the right place. Please, remove.
L199. Did the Authors create a definition of “reliable” for their purposes? The distinction between reliable and not reliable data can be vague. Please, specify what you mean by “reliable”.
L356. I suggest considering replacing “got” with “became”.
Citation: https://doi.org/10.5194/hess-2021-103-RC1 -
AC1: 'Reply on RC1', Francesca Zanetti, 21 Jun 2021
We thank the referee for providing useful comments which will help us to improve our Manuscript. In what follows, we provide a detailed set of responses to the reviewers' comments (italic text).
General comment
This is an interesting work based on a relatively new technology that has the potential to offer new insight into the dynamics of stream network in small catchments. Streams are highly dynamic systems and characterizing such “liveliness” is important to move towards a better understanding of catchments’ functioning.
Overall, this manuscript is well written, logically structured, and clearly illustrated. The introduction is solid and the results are well supported by the data. I have only some specific and minor comments that I would like the Authors to address (see below). Overall, I recommend a moderate revision but since the MODERATE grade is not available to HESS reviewers, I indicated minor review in the review form.
We thank the referee for the positive assessment of our manuscript.
Specific comments
L78-83. I think that the three research questions could be more appealing than the current ones. I understand that, as this is a relatively new measurement approach, aiming at providing the reader with some methodological observation is useful for possible further application of the methods. So, I like question 1. However, I would move it as last question because, in my opinion, it’s more important to focus on processes that the method is able to describe and understand, rather than on the method itself. The second question is a bit too narrow because it implies a “yes” or “nor” and does not lead to much insight into hydrological processes. Similarly, question 3 sounds a bit too “methodological” and not so oriented towards process understanding. So, I suggest moving question 1 as third, and to rephrase question 2 and 3 and bit.
Thank you for the suggestion. We will rephrase the second and third research questions, and we will change their order as suggested. In particular, we plan to reformulate the research question n.2 in order to avoid it to be a closed-answered question. This could be rephrased as "how is the persistency of individual nodes of the network reflected by the statistical features of ER signals across a range of different substrates and various degree of flow intermittency?". Likewise the third question could be rephrased as follows: "how do wet length and catchment discharge co-evolve in response to a sequence of rain events, and how does their mutual relation depend on the temporal resolution of the available observations?"
Indeed, my overall impression is that this work is much technically- and methodologically-oriented and less prone to describe and understand hydrological processes, and I guess that reader of this journal are more interested in knowing how catchment works rather than know if some data can fit a certain model or no. So, I invite the Authors to consider revising the work to reflect this aspect.
Thank you for the comment. We trust our Ms is in line with the scope of HESS, as the aim and scope of the journal explicitly mentions "the study of the spatial and temporal characteristics of the global water resources (solid, liquid and vapour) and related budgets, in all compartments of the Earth system (atmosphere, oceans, estuaries, river, lakes and land masses), including water stocks, residence times, interfacial fluxes, and the pathways between various compartments". Moreover, in the aim and scope section of the journal it is explicitly claimed that "papers should contribute to the advancement of hydrological modeling, hydrological monitoring and data analysis, [...], experimental design and technology". For these reasons we believe that our work could be of interest for the readership of HESS. Furthermore, we would like to emphasize that we are aware of several highly cited, technically-oriented/methodological papers recently published by the journal (i.e. Kaplan et al., 2020 (DOI: 10.5194/hess-24-5453-2020); Kirchner 2016 part 1 (DOI: 10.5194/hess-20-279-2016) and part 2 (DOI: 10.5194/hess-20-299-2016)), an instance which makes us believe that being methodological does not represent per se a shortcoming for an HESS Manuscript. At the same time, we acknowledge the importance of deciphering the underlying physical mechanisms behind the empirical data and we have tried to provide some insight on the major hydrological processes underlying our data (lines 379-387). We are also willing to further emphasize the implications in terms of processes in the revised version of the paper (see also reply to the second referee), as per the suggestion of the referee.
L127. What are the criteria for the choice of the field deployment? Please, specify.
Thank you for the comment. Sensors were deployed along the most dynamical part of the tributaries that were observed during the field surveys, in a catchment with a size that makes the field surveys feasible with the human resources available (say, catchment area < 3 km2). Technical difficulties and the time needed to reach some locations (due to the harsh environment) were thus taken into account in the selection process. The specific location of the sensors was also chosen considering the heterogeneous substrates of the catchment, making it possible to analyze the sensors' behaviour on different type of soils, and ensuring an heterogeneous distribution of the nodes' persistency. These points will be better specified in the revised text.
L145. This sounds a very short period to me. I understand that practical issues might be arisen but typically we need a longer time period to observe hydrologically processes that often highly variable in time. Si the average rainfall and stream flow in this period comparable to long-term rainfall and stream flow (or compared to the other years where observations are available, since this is a relatively new experimental catchment). I think it’s important for the Authors to discuss this issue, explain why (if) they reckon this is a suitable spell and why, and why (if) this is a representative period for derive information on the hydrological functioning of this catchment. Moreover, the should discuss how this short period of time can potentially impact on the results.
Thanks for the comment. The length of the study period is affected by the characteristics of the study catchment, which is a high relief catchment located in the southern side of the Alps. The basin is snow-covered usually from late November until the end of spring (June). This constraints the time window that can be used to take field measurements and analyze the underlying network dynamics, as winter renews the underlying network dynamics every year. Moreover, the sensors' deployment is highly time-consuming and the set up of all the sensors require at least some months. During the 2019 the underlying climatic conditions were particularly unlucky, as the winter season came earlier than expected (early November instead of early December). At any rate, it is important to stress that though relatively short, the study period covers a wide range of climatic conditions and network configurations, as discussed below. A detailed statistical analysis was performed to analyse the representativeness of the climatic conditions observed during the period of record. Figure 1 of the supplement material compares the distribution of the daily rainfall depths (h) observed during the reference time window (Sep and Oct of 2019) and that observed in the long term (2010-2020) during the whole period within which the network is dynamical (from July 1 to November 30). The comparison shows that the frequency distribution of the rain depth for our study period matches quite well the corresponding long-term distribution, an instance which suggests that the rain regime during the analyzed period is in line with that driving the observed longer term network dynamics. The significant variations of the underlying hydro climatic conditions during the reference time window is also manifested in the statistic of the cumulated rainfall over 5 and 35 days, which have been shown to be the main drivers of the total active length in the Valfredda (Durighetto et al. 2020). Our analysis indicates that the 5-days cumulative rainfall depth (h5) during the study period ranges from 0 up to the upper quartile of the long-term distribution of h5 within the years 2010 – 2020, indicating that the event variations in rain amounts are well reproduced by the study period. Likewise, for h35 the min/max observed range during the period of record (Sep - Oct 2019) is comparable to 2σh35 (σh35 being the long-term standard deviation of h35). Therefore, while the wet summer could not be included because of the time needed to deploy the sensors and the study period spans only two months, we think it is representative of the longer term network dynamics observed in the Valfredda. All these arguments will be added to the revised Ms. Please see also our response to the same issue in the section dedicated to the second referee.
L184. Where was the sensor placed? In the grass in a convergent zone, I guess, where water was not flowing? Please, give more information on the aspect of sensors deployment in the field.
Thanks for the comment. The sensor was placed on the grass where water flow was not observed permanently, as the channel activated only during and after precipitation events. In the revised text, we will describe more comprehensively the deployment in the field of the sensors.
Minor comments
L5. I suggest considering the term “customized” instead than “personalized”.
Agreed, thanks. We will modify the text as suggested.
L6. I suggest removing “analysed,”.
Agreed, thanks.
L6. The expression “nodes’ persistency” is not clear without reading the manuscript. Please, clarify.
Thank you for the comment. We will clarify this point by saying "nodes' persistency (i.e. a proxy for the probability to observe water flowing over a given node, see e.g. Durighetto et al., 2020)".
L37. Perhaps here the citation to Godsay and Kirchner (2014) fits well. In any case, a more recent references would be a nice complement.
Thank you for the suggestion. We will add more recent references including the following: Floriancic et al., 2018 (DOI: 10.1002/hyp.13302); Godsey and Kirchner 2014 (DOI: 10.1002/hyp.10310); Jaeger et al., 2007 (DOI: 10.1007/s00267-005-0311-2); Lovill et al., 2018 (DOI: 10.1029/2017WR021903).
L129. “heterogeneous persistencies so as to avoid redundancy in the data.” This is not clear to me, please specify.
Thank you for the comment. What we mean here is that the sensors were deployed close to nodes that had a different degree of persistency in order to assess the consistency of the sensors in various field conditions (from nodes almost always dry to nodes that dries down only sporadically). We will better specify this point in the revised text.
L160. Which observations? Please specify. What does “some modelling” mean? Please, explain.
Thank you for the comment. The models referred to here are described in the paper cited in L161, Durighetto and Botter 2021 (DOI: 10.1002/hyp.14053). More information will be added in the revised text, including a brief summary about the model formulation and the main assumptions.
L176. Building a reliable flow rating curve, especially in mountain catchments, is a challenge. How many points were collected to build the FRC, and what was the range of stream flow values and the resulting goodness of fit measures? In other terms, is the FRC reliable to infer robust stream flow values? Please, explain.
Thank you for pointing that out, we agree on the comment. The rating curve was built using 7 points with discharge ranging between 9 and 300 l/s and a coefficient of determination R2 = 0.99. We will include these data in the revised text.
L178-179. I understand that the Authors want to give light to the ERC project but this is not the right place. Please, remove.
Agreed, thanks.
L199. Did the Authors create a definition of “reliable” for their purposes? The distinction between reliable and not reliable data can be vague. Please, specify what you mean by “reliable”.
Thank you for pointing that out. Non reliable data are meant to be missing data, no-data and zeros induced by malfunctioning of the sensors. We acknowledge that the term "reliable" could be considered as misleading and subjective and we will remove it from the text. In the revised version we will specify the reasons for which some data were not included in the analysis, as indicated above.
L356. I suggest considering replacing “got” with “became”.
Agreed, thanks.
-
AC1: 'Reply on RC1', Francesca Zanetti, 21 Jun 2021
-
RC2: 'Comment on hess-2021-103', Anonymous Referee #2, 13 May 2021
This paper reports on what appears to be a pilot project using electrical conductivity sensors to detect presence/absence of flow in a small Alpine basin. I appreciate the difficulty in making these kinds of field measurements, and the paper does a good job of highlighting some of the challenges involved.
However, the abstract, discussion, and conclusions make some rather confident claims, without acknowledging the severe limitations of the data and the unusual characteristics of the Valfredda study site. These issues substantially undermine the strong claims presented here. The major issues are enumerated below.
1 – I am sure that these field data were hard-won, but they span only two months (or maybe only one month – Figures 6 and B2 both refer to "the study period" but one is only about half as long as the other…?), and include only a small handful of precipitation events. It is hard to draw robust conclusions from such limited evidence. The study by Jensen et al. (2019, cited in the references) provides an illustrative contrast, with a much more extensive set of observations, and thus more robust inferences, drawn from a similar number of sensors along a similarly sized channel network (but a longer study period with more precipitation events). I will leave it to the editors to decide whether HESS wants to publish such a limited data set – speaking for myself I would have waited for a more comprehensive picture to emerge.
2 – The resulting uncertainties are very large (see figure 9), but this is not adequately accounted for in the presentation. The text (line 327) says that b varies by about 1% as the temporal resolution changes, but given that the uncertainty in b can be over 10%, it is actually unknown how stable b really is (or isn't). The text (line 331) even argues for a systematic increase in R^2 from 0.485 to 0.522, even though the uncertainty in R^2 can be over 20%, making this "systematic" increase statistically meaningless. Even these very large uncertainties may be underestimates, because the underlying data are serially correlated, meaning that (for example) few of the points in Figure 8 are statistically independent of one another. From the methods it is unclear whether this has been taken into account, as it should be.
The last main conclusion of the paper is that (lines 411ff): "The mean value of the exponent of the power law relationship between catchment discharge and total active length was found to be almost independent on the frequency of the observational data, which instead had a larger impact on the goodness of fit of the power-law model. When the frequency of the data is lower, the observed values of R2 are, on average, larger…" In view of the vast uncertainties in Fig. 9, these conclusions are reckless. Within the uncertainties, either of these trends could be strongly increasing, strongly decreasing, or zero. There is simply no robust conclusion that can be drawn from the data.
3 – The limitations of the study site are severe, particularly for analyses of network dynamics. The basic problem is that roughly 80% of the basin seems to have no surficial drainage network at all, consisting instead of talus slopes and moraines. The critical issue here – which is not acknowledged anywhere in the paper – is that this ~80% of the basin is still generating discharge (at least some of which is presumably measured at the outlet), but the accompanying network dynamics are invisible because they are occurring beneath piles of rock debris. Outside of the mapped network there appears to be roughly two square kilometers of drainage area with no surficial drainage at all.
At best, that means that any observations here cannot be compared with the rest of the network dynamics literature, in which the discharge from the whole basin is compared with the flowing stream network across the entire basin. Thus, for example, there is no way to compare Figure 8 with similar diagrams from other studies, because in this case most of the discharge appears to be generated by subsurface flow that is presumably strongly damped and lagged, suppressing the variability in Q (this may account for the sharp vertical lines in Figure 8, for example).
The manuscript doesn't confront (or even disclose) this problem anywhere, which is surprising given the abstract's mention of "the diversity of the hydrological behaviour of the study catchment" – by which the paper seems to mean only the two small drainage networks that were studied, not the other roughly three-fourths of the catchment.
It is virtually a truism in catchment studies that each site has its own idiosyncrasies, but here this particular "uniqueness of place" makes a network dynamics study particularly difficult. Why study network dynamics in a catchment where the great majority of the drainage area has no network at all? Such a site makes it particularly difficult to draw any mechanistic inferences from the observed network behavior.
The manuscript says (lines 367ff): "Network length was found to be more sensitive than discharge to small precipitation inputs: while most rain events induced visible changes in the active channel length, the catchment stream flow was sensitive only to the rain events lasting for several consecutive days (6-9/09, 13-18/10, 20-24/10) and to intense storms (more than 20-30 mm in 9-12 hours)." This is exactly the behavior that one would expect from a field site like this one, with most of the discharge being generated by relatively slow subsurface flowpaths over ~80 percent of the catchment, but with network lengths being measured on the very few surface drainages in the remaining small fraction of the catchment.
The manuscript continues (lines 375ff): "In our case study, the standard deviation of the wet length as derived from the sensors’ data is 360 m, while the standard deviation of L predicted by the power-law model based on the observed variability of the discharges is only 224 m (about 40% lower). This underestimation is induced by the poor ability of the power law model to capture the observed network dynamics produced by small precipitation inputs." It would rather seem that the problem is that *no* model could possibly capture the relationship between the network dynamics in a small part of the catchment, and the discharge generated by completely different mechanisms in the great majority of the catchment.
All of the conclusions concerning the relationship between stream length and discharge (essentially everything after line 10 in the abstract and after line 406 in the conclusions) are based on very thin data from a catchment in which discharge mostly comes from subsurface flow through rock debris (with the result that changes network length in the small fraction of the catchment with surface drainage are unsurprisingly not clearly related to the discharge, which mostly comes from the rest of the catchment). Thus all of those conclusions are based on very thin data that does not allow straightforward interpretation even in this study catchment, and cannot be extrapolated to the great majority of catchments that lack this particularly exotic geometry.
If those conclusions are excluded – as they really should be, given their weak empirical support and their inherently problematic interpretation (network lengths are not measured in the part of the catchment that generates most of the discharge) – then we have essentially a technical note outlining a new way to deploy conductivity sensors, and conveying some lessons learned from a first deployment of these sensors. That would seem to be a more appropriate way to go, rather than trying to draw strong conclusions about length/discharge relationships from such limited data and such a problematic study site.
As noted by at least another reviewer, the language would also need work (e.g., "customized" rather than "personalized"), but any revised manuscript is likely to be substantially different so I have not marked those issues in this go-around.
Citation: https://doi.org/10.5194/hess-2021-103-RC2 -
AC2: 'Reply on RC2', Francesca Zanetti, 21 Jun 2021
We thank the referee for providing useful comments which will help us to improve our Manuscript. In what follows, we provide a detailed set of responses to the reviewer' comments (italic text).
This paper reports on what appears to be a pilot project using electrical conductivity sensors to detect presence/absence of flow in a small Alpine basin. I appreciate the difficulty in making these kinds of field measurements, and the paper does a good job of highlighting some of the challenges involved.
We thank the reviewer for acknowledging the value of this work.
However, the abstract, discussion, and conclusions make some rather confident claims, without acknowledging the severe limitations of the data and the unusual characteristics of the Valfredda study site. These issues substantially undermine the strong claims presented here.
The data have some issues, which are explicitly acknowledged and discussed in the text (Section 2.2.3 and Appendix B). We also respectfully disagree on the claim that the study catchment is so unusual. Most of the headwater catchments in the Dolomitic region, and in general in the southern side of the Alps, share similar features in terms of substrate heterogeneity and possible presence of internal karst regions. Nevertheless, the revised version of the paper will better emphasize the specific geological features of the study catchment, the duration of the study period, and the impact of these factors on the main results of the paper (see below).
1 – I am sure that these field data were hard-won, but they span only two months (or maybe only one month – Figures 6 and B2 both refer to "the study period" but one is only about half as long as the other…?), and include only a small handful of precipitation events.
Thanks for noting that. The study period actually spans two months, as seen in the plots of Figure 5 and explicitly stated in the method section. Figure B2 shows only the period during which data of the sensor S39 were corrected, we will modify the caption of the figure accordingly.
It is hard to draw robust conclusions from such limited evidence. The study by Jensen et al. (2019, cited in the references) provides an illustrative contrast, with a much more extensive set of observations, and thus more robust inferences, drawn from a similar number of sensors along a similarly sized channel network (but a longer study period with more precipitation events). I will leave it to the editors to decide whether HESS wants to publish such a limited data set – speaking for myself I would have waited for a more comprehensive picture to emerge.
The length of the study period is affected by the characteristics of the study catchment which is a small high relief catchment located in the Alps. The basin is snow-covered from late November until the early summer (June). This constraints the time window that can be used to take field measurements and study the network dynamics (potentially to 4-5 months) and renews the underlying network dynamics every winter. Moreover the sensors' deployment is highly time-consuming and the set up of all the sensors require at least some months. The 2019 fall season was particularly short, as the snow season came earlier (early November instead of early December). At any rate, it is important to stress that, though relatively short, the study period covers a wide range of climatic conditions and network configurations, as discussed below. A detailed statistical analysis was performed to analyse the representativeness of the climatic conditions observed during the period of record. Figure 1 of the supplement material compares the distribution of the daily rainfall depths (h) observed during the reference time window (Sep and Oct of 2019) and that observed in the long term (2010-2020) during the whole period within which the network is dynamical (from July 1 to November, 30). The comparison shows that the frequency distribution of the rain depth for our study period matches quite well the corresponding long-term distribution, an instance which suggests that the rain regime during the analyzed period is in line with that driving the observed longer term network dynamics. Thus, we do not see any particular reason for which the results of the paper could be strongly impaired by the specific duration of the reference time window. This is especially true if one considers the diversity of the hydrological conditions observed during the period of record. More specifically, the main issues associated to the use of water presence sensors highlighted in the text would not be changed by the use of more data, and the nature of the relationship between the mean persistency of the nodes and statistical properties of the ER signal is unlikely to be significantly modified as well. Likewise, the procedure identified for reconstructing the high frequency space-time dynamics of the stream network using ER data would be the same, even if more data were considered. Also the hysteresis in the Q vs. L relationship would remain, as they are observed within individual events, and across all the events. Of course, it would have been great to have a longer study period. Based on the duration of the snow season and the time needed for the sensors deployment, we argue that the maximum possible duration of the period of record in an Alpine field site of this type could hardly exceed 3 months, after which the winter would freeze the catchment and renew the underlying network dynamics. However, repeating this experiment to get 3 months of data instead of 2 is practically unfeasible, and would imply a huge experimental effort (e.g. another deployment of all the sensors), eventually leading to results that are unlikely to be significantly different from those obtained in this study. All these arguments will be clarified in the revised version of the paper. In any case we do not see significant discrepancies between our results and the results shown by Jensen et al., 2019: in both cases there is a good consistency between ER data and the visual observations of the network, expansions/contractions of the network occur by growth of disconnections within the streams and different types of hysteresis between L and Q are observed.
2 – The resulting uncertainties are very large (see figure 9), but this is not adequately accounted for in the presentation. The text (line 327) says that b varies by about 1% as the temporal resolution changes, but given that the uncertainty in b can be over 10%, it is actually unknown how stable b really is (or isn't). The text (line 331) even argues for a systematic increase in R2 from 0.485 to 0.522, even though the uncertainty in R2 can be over 20%, making this "systematic" increase statistically meaningless.
We think the referee is highlighting the fact that the estimate of the value of b and the agreement of the power law model with the data can be dependent on the specific dates of the surveys. As per the b exponent, we believe that the pattern shown by the mean scaling exponent b is quite interesting, and that variations of the order of 10% could be considered as acceptable in the light of the simplicity of the power law model and the presence of possible measurements errors. We acknowledge, however, that the pattern of R2 is less meaningful owing to the large standard deviation, as suggested by the referee. In the revised version of the paper, we will discuss more explicitly this point about the variability of b and R2 across the different sub samples of the data. Thanks for noting the issue.
Even these very large uncertainties may be underestimates, because the underlying data are serially correlated, meaning that (for example) few of the points in Figure 8 are statistically independent of one another. From the methods it is unclear whether this has been taken into account, as it should be.
We are puzzled by this comment. Serial correlation is unavoidable in all high frequency joint data set for L and Q. Here we are exactly reproducing a fictitious sampling campaign with sporadic measurements of active length and discharge, assuming that a standard power law is applied to the sampled data, as typically done in all previous studies where the relationship between Q and L was studied. The standard deviation shown in this Figure is a measure of the diversity of model parameters as a function of the sampling dates, not a standard model uncertainty associated to a given data set - which relies on some likelihood function, with possibly correlated residuals. This standard deviation simply resembles the heterogeneity of the fitting of a power-law model done over a set of different sub-samples of the same set of data. In the revised paper we will more explicitly discuss how this analysis should be interpreted to avoid misunderstandings of these kind. Thanks for bringing this issue to our attention.
The last main conclusion of the paper is that (lines 411ff): "The mean value of the exponent of the power law relationship between catchment discharge and total active length was found to be almost independent on the frequency of the observational data, which instead had a larger impact on the goodness of fit of the power-law model. When the frequency of the data is lower, the observed values of R2 are, on average, larger…" In view of the vast uncertainties in Fig. 9, these conclusions are reckless. Within the uncertainties, either of these trends could be strongly increasing, strongly decreasing, or zero. There is simply no robust conclusion that can be drawn from the data.
Thanks for the comment. In the light of the fact that the R2 are not particularly high and the available samples of Q and L are relatively small, we deem that a range of variability of b of about 10% could be acceptable. Therefore, we maintain the point of the relative stability of b for different values of T. On the other hand, we recognize that the observed pattern in R2 might not be particularly meaningful, in view of the standard deviation of the estimate shown in Figure 8. This will be emphasized more clearly in the revised version of the paper.
3 – The limitations of the study site are severe, particularly for analyses of network dynamics.
We respectfully disagree on this claim, for the reasons detailed in this rebuttal.
The basic problem is that roughly 80% of the basin seems to have no surficial drainage network at all, consisting instead of talus slopes and moraines. The critical issue here – which is not acknowledged anywhere in the paper – is that this ~80% of the basin is still generating discharge (at least some of which is presumably measured at the outlet), but the accompanying network dynamics are invisible because they are occurring beneath piles of rock debris. Outside of the mapped network there appears to be roughly two square kilometers of drainage area with no surficial drainage at all.
The comment is certainly relevant - and we are grateful to the referee for the insightful input - but the estimate of the percentage of catchment area drained by the existing channel network indicated by the referee is largely underestimated. In fact, in the upper Valfredda catchment 1.7 km2 out of 2.6 km2 (total catchment area) directly drains through the hydrographic network (65%) according to the analyses of the drainage flow paths performed with a high resolution DTM; thus, 35% of the catchment has no surface drainage network - not 80% as guessed by the reviewer.
At best, that means that any observations here cannot be compared with the rest of the network dynamics literature, in which the discharge from the whole basin is compared with the flowing stream network across the entire basin. Thus, for example, there is no way to compare Figure 8 with similar diagrams from other studies, because in this case most of the discharge appears to be generated by subsurface flow that is presumably strongly damped and lagged, suppressing the variability in Q (this may account for the sharp vertical lines in Figure 8, for example).
Let us note that the drainage density in the upper Valfredda is comparable with that of other study catchments used to study network dynamics (e.g. South Fork of Potts Creek, Fernow, Turbolo). Nevertheless, we recognize that the presence of a karst area which originates a localized spring that releases a quite constant discharge of about 40 l/s (roughly 30% of the mean discharge) and feeds a perennial stream needs to be taken into account when analyzing the relationship between active length and discharge. This will be done in the revised version of the paper, as explained below.
The manuscript doesn't confront (or even disclose) this problem anywhere, which is surprising given the abstract's mention of "the diversity of the hydrological behaviour of the study catchment" – by which the paper seems to mean only the two small drainage networks that were studied, not the other roughly three-fourths of the catchment.
Thanks for the comment, we will add a detailed description of the impact of the portion of the catchment without drainage network on the main results of the paper (see below). Let us note that the percentage guessed by the referee is largely overestimated (as three-fourths should read one third).
It is virtually a truism in catchment studies that each site has its own idiosyncrasies, but here this particular "uniqueness of place" makes a network dynamics study particularly difficult. Why study network dynamics in a catchment where the great majority of the drainage area has no network at all? Such a site makes it particularly difficult to draw any mechanistic inferences from the observed network behavior.
The Valfredda catchment is one of the study sites of the ERC project Dynet, because it is quite representative of the headwater catchments in the Dolomitic region. It was used in this study because of the heterogeneity of the substrate, and the availability of long term data about rainfall, discharge, water quality and persistency of the nodes of the network, which lie at the basis of several results presented in the paper (Figures 3 and 4). We respectfully disagree on the claim that mechanistic inferences can hardly be made here, as indicated in the existing literature as the study by Durighetto et al. (2020); in the cited paper three different models were developed and statistically validated in order to describe the dynamics of the active drainage network length (ADNL) starting from the wet length values measured during nine surveys in the field, the atmospheric forcing and the geologic characteristics of the catchment. Both empirical data and the models results showed the influence of the antecedent precipitation on the measured/modelled dynamics of the stream network as much as the geologic features of the study catchment. All these arguments about the choice of the site will be further stressed in the revised text.
The manuscript says (lines 367ff): "Network length was found to be more sensitive than discharge to small precipitation inputs: while most rain events induced visible changes in the active channel length, the catchment stream flow was sensitive only to the rain events lasting for several consecutive days (6-9/09, 13-18/10, 20-24/10) and to intense storms (more than 20-30 mm in 9-12 hours)." This is exactly the behavior that one would expect from a field site like this one, with most of the discharge being generated by relatively slow subsurface flow paths over \~80 percent of the catchment, but with network lengths being measured on the very few surface drainages in the remaining small fraction of the catchment.
Thanks for the comment, which gives us the opportunity to clarify the issue. First of all, in the upper Valfredda catchment analyzed in this paper the karst area does not contribute to the large majority of the discharge but only to approximately 30% of it. In fact, the water that infiltrates in the region where terrain depressions and karst substrates dominate mostly exfiltrate through a permanent spring in the northern portion of the network, originating a constant discharge contribution of about 40l/s (while the average observed discharge exceeds 130 l/s). The large majority of the total discharge is instead released through the dynamical drainage network mapped in this study. Moreover, we contend that the presence of a permanent water source (Q_p) which feeds a perennial stream of length L_p is able to produce the hysteresis in the L vs. Q relationship shown by our data. Within a single intense rain event, the hysteresis depends on the fact that the active length increases faster than Q in the early stages of the event, while decreases much slower than Q in the recession. Across different events, instead, the hysteresis is generated by shifts in the response of Q and L to different types of rain events. In particular, when rain events are moderate the channel network activates but the amount of water conveyed to the outlet is limited, owing to disconnections and limited flow velocities; conversely, when rain events are more intense the same active length contributes a much larger discharge to the outlet. Anyways, we think that the point made by the referee about the role of a permanent water source for the analysis of L vs. Q relationships is quite interesting. To further investigate the impact of a perennial channel with a constant wet length (say, L_p) supplied by a permanent constant discharge Q_p on the underlying L vs. Q relationship, we have performed a numerical simulation. Therein, we assumed that the total active length is the sum of a constant active length L_p and a dynamical stream network length L_d(t), while the total discharge was assumed to be the sum of a dynamical discharge Q_d(t) and a constant discharge Q_p. The dynamical length was then assumed to be linked to the dynamical discharge through an exact power law relationship of the type L_d(t) = aQ_d(t)b. In our exercise, we studied the relationship between the total length L_T(t) = L_d(t) + L_p and the total discharge Q_T(t) = Q_d(t) + Q_p, when different values of Q_p and L_p are used. In particular, we focus on the observed values of the exponent of the power law model properly fitting the overall values of Q_T and L_T (Figures 2 to 4 of the supplement material of this rebuttal), which was compared to the b value of the dynamical power law relationship linking L_d(t) and Q_d(t). The plots shown in Figure 2 to 4 of this rebuttal indicate that the variations of the apparent b induced by the presence of a permanent discharge and a permanent active length are in the range of 10/20%. Such variations are particularly low when Q_p and L_p follow the proportion of the power law-model holding for the dynamical fractions of Q and L - as observed in our dataset, where Q_p and L_p were estimated based on field measurements. In all cases, however, the R2 doesn't change considerably, indicating that the goodness of fit of the power law model is not impacted by the presence of a positive value of Q_p and L_p. All these arguments clearly indicate that the presence of the karst area in the northern part of the catchment doesn't affect significantly the fitting of the power-law model and the extent of the observed hysteresis in the L vs. Q relationship. This also applies to the time series of L and Q studied in the Ms. All these arguments will be included in the revised text, in which we plan to isolate the impact of Q_p and L_p on the fitting and the scaling exponent and we will be able to quantitatively demonstrate that the presence of a constant discharge feeding a perennial channel does not change the goodness of fit of the power-law model and has a moderate impact on the exponent b. We thank the referee for raising this important point.
The manuscript continues (lines 375ff): "In our case study, the standard deviation of the wet length as derived from the sensors’ data is 360 m, while the standard deviation of L predicted by the power-law model based on the observed variability of the discharges is only 224 m (about 40% lower). This underestimation is induced by the poor ability of the power law model to capture the observed network dynamics produced by small precipitation inputs." It would rather seem that the problem is that *no* model could possibly capture the relationship between the network dynamics in a small part of the catchment, and the discharge generated by completely different mechanisms in the great majority of the catchment.
As pointed out above, the discharge generated by the unchanneled part of the catchment is about 40 l/s, while the mean discharge is more than three times larger (130 l/s). The reasons behind the observed hysteresis in the L vs Q relationship are described above, and pertain to differences in the velocity and connectivity of the dynamical streams across different events -- or asymmetries in the rate of change of L and Q within individual events. At any rate, we will include an in-depth analysis of how the presence of a permanent spring which feeds a perennial stream impacts the discharge vs. active length relationship, as indicated above.
All of the conclusions concerning the relationship between stream length and discharge (essentially everything after line 10 in the abstract and after line 406 in the conclusions) are based on very thin data from a catchment in which discharge mostly comes from subsurface flow through rock debris (with the result that changes network length in the small fraction of the catchment with surface drainage are unsurprisingly not clearly related to the discharge, which mostly comes from the rest of the catchment).
As pointed out above, the permanent spring through rock debris contributes to 40 l/s while the average discharge during the study period is 130 l/s. See also our previous responses on the same point.
Thus all of those conclusions are based on very thin data that does not allow straightforward interpretation even in this study catchment, and cannot be extrapolated to the great majority of catchments that lack this particularly exotic geometry.
We believe that our conclusions could be of interest for the HESS readers, and not necessarily restricted to the Valfredda site. In fact, hysteresis in the Q vs. L relationship have been previously noted also in other catchments (Zimmer et al. 2017 DOI: 10.1002/2016WR019742; Barefoot et al. 2019 DOI:10.1029/2018WR023877; Jensen et al. 2019) though the underlying reasons were not discussed in-depth. However, the potential impact of the specific features of the study catchment on the results will be better emphasized in the revised text, as indicated above.
If those conclusions are excluded – as they really should be, given their weak empirical support and their inherently problematic interpretation (network lengths are not measured in the part of the catchment that generates most of the discharge) – then we have essentially a technical note outlining a new way to deploy conductivity sensors, and conveying some lessons learned from a first deployment of these sensors. That would seem to be a more appropriate way to go, rather than trying to draw strong conclusions about length/discharge relationships from such limited data and such a problematic study site.
The potential role of the specific features of the study catchment will be better emphasized in the revised text, as per the referee's suggestions. As per the indication of transforming this paper into a technical note, we can not agree. There is a strong length constraint for HESS technical notes (a few pages) that makes this choice unsuited to this pretty long Ms. In particular, we do believe that the analysis of the joint response of Q and L to the observed precipitation forcing (Figures 7 and 8) should not be removed because they are an important outcome of our field campaign. Moreover, both the reconstruction of the spatial and temporal dynamics of the active stream network (Figure 5) and the analysis of the relationship between the statistical properties of ER signal and persistencies (Figure 4) represent important outcomes of the study, which complement our experimental data set and the discussion of the advantages and disadvantages of the sensors' deployment. Incidentally, the text quite clearly acknowledges the inherently specific nature of the findings presented in the Ms (e.g. "our results indicate that in some cases the use of a bijective L vs. Q function to infer active length changes in catchments where discharge time series are available, might lead to significant underestimation of the actual variations of the flowing channel network. In our case study, the standard deviation of the wet length as derived from the sensors’ data is 360 m, while the standard deviation of L predicted by the power-law model is [...] 40% lower)."). Nevertheless, in the revised text we will try to emphasize even more that these results might not necessarily apply to all the catchments in the word, and that further investigations are needed to determine the generalizability of the behaviour observed in the Valfredda in the light of the uniqueness of the place.
As noted by at least another reviewer, the language would also need work (e.g., "customized" rather than "personalized"), but any revised manuscript is likely to be substantially different so I have not marked those issues in this go-around.
Thanks for the comment, the text will be polished and revised (and customized will replace personalized).
-
AC2: 'Reply on RC2', Francesca Zanetti, 21 Jun 2021