General comments

The authors propose an automatic approach to infer and classify the dynamic characteristics of karst systems based on the analysis of the recession hydrograph. The approach introduces variability through multiple recession extraction methods, fitting individual recession segments and optimization approaches with various degrees of freedom for a fast and slow flow recession model. Recognition of the duality of porosity and processes in karst systems leads to a vision of a generally accepted two-reservoir conceptual model. The desire to develop a method of recession analysis that is consistent with this view is well justified. The automated and multi-angle traits of the approach are indeed essential to cope with the known biases of single and visually supervised approaches to recession analysis. In my opinion, this work will improve the robustness of the comparison of the dynamic characteristics of karst hydrological systems. Furthermore, I believe that tools for comparing hydrological systems are an area of research that will continue to develop and contribute to a better understanding of hydrological systems. For these reasons, I consider that the article is following the scope of HESS.

In general, the paper is intelligible. The results are interpretable and make sense. The introduction correctly sets the context, and the discussion provides elements that facilitate the interpretation of the results. Nevertheless, I find that these two sections do not motivate well enough the choices that have been adopted and do not give a good account of the strengths and weaknesses of the approach. As a result, the modeling philosophy and epistemological approach are a bit blurry and sometimes inconsistent (see specific issues). Clearly stating these choices would allow the discussion section to be more structured, sharpened, and provide better perspectives and research avenues. Currently, I found the latter to be relatively weak or absent. Still, I enjoyed the approach and the content, particularly the idea of varying the degrees of freedom in the model optimization. I feel that this is a reliable piece of work that I could reuse myself and will certainly be useful to other researchers. Therefore, I consider the scientific quality to be good. Since it is a paper that proposes a method, I also regret the lack of additional material that would facilitate re-implementation (input/output and/or code recipes). I would also have appreciated having access to a bit more details about the results, not especially in the paper itself but in the appendix or the supplementary materials, about the seasonality of parameters and the statistics of recession segments.

The figures are clear and good. The introduction and discussion sections are ok but could probably be improved and sharpened. But in my opinion, this is more a question of substance than form (see above and specific issues). However, I regret the long and complex sentence constructions, plus a large number of grammatical errors. I think they are below acceptable standards and could be easily removed with a good grammar/spellchecker.

Reply to General Comments: We thank reviewer 1 for her/his useful and valuable comments that will help to improve the manuscript, in particular its readability. We will carefully rework the entire text in order to improve its structure and remove typographical errors. We will improve the introduction and discussion sections by accounting more detailed on the core strengths and weakness of the recession analysis approaches. We will make the R-codes for recession analysis as well as input data (spring hydrographs) publicly available via our Github repository. We will provide further details of the results of our analysis in the appendices. All suggestions will be taken into account. According to her/his comments, we will perform the following changes.

Specific issues

As stated in my general comment, I have very few concerns about the approach itself (except perhaps SI-2), and I have a favorable opinion about the paper. Still, I am concerned about how the methodological choices are introduced and motivated and how the results are interpreted. It lacks coherence and proper motivation, which in some respects affects the clarity of the reading, might suggest that the method is inappropriate, and prevents the development of clear research perspectives. I present below the nature of my concerns and my suggestions on how they can be addressed in a motivated and coherent rationale as rooms for improvement through possible minor revisions of the introduction and discussion sections.

 SI-1: Motivation for individual recession event: seasonality.

There is a first inconsistency in the motivation for fitting individual recession events. It is said that fitting individual recession events would allow capturing the variability for better inference of the structural and dynamic traits of a system (L48-51) and would be more consistent because the parameters obtained from the entire set of recession points could not be actually representative of individual cases (L216-217). Yet, in the other section, the authors tend to interpret the mean of all parameters and indicators as representative of the karst system. In the meantime, the spread of the parameters is deemed both necessary and valuable as it reflects dynamic properties but, paradoxically, also considered as unwanted uncertainty that we shall reduce for having a set of "representative" parameters (L25-27, L450-460, L525-526). I discuss the author's suggestion for reducing uncertainties in SI-5.  For now, l continue with the value of the individual segments approach.

Historically, recession analysis methods applied on all points pursuit the same goal of inferring the dynamic and structural properties of the system. It is also possible to have uncertainty bounds, analytical or bootstrapped ones, for estimated parameters. So, why go for individual recession fitting if you can already depict the average and the variability of the estimated parameters? The answer is that those who have adopted this approach have learned something about the dynamics of the system. The parameters were projected on the temporal/seasonal axis, and interesting patterns related to the seasonality or the antecedent conditions of hydrological variables appeared [See 1-5]. In doing so, the individual segment approaches proved themselves useful and showed that the low dimensional models that we commonly use (Mangin is no exception) are actually underfitting recession dynamics [6].

Reply to SI-1 (Paragraph 1 and 2)

We thank the reviewer for the useful and critical reviews. We agree that the premise given for individual recession segment (IRS) analysis seems inconsistent with our suggestion for reducing parameters uncertainties. However, L48-51 and L216-217 were meant to provide justification for our methodological choice of using the IRS analysis approach. Our focus is to explore the possibility of identifying intrinsic karst aquifer dynamic with IRS analysis. Even though karst systems are very heterogeneous, yet the systems’ geometry are somewhat fixed – although dissolution of carbonate rocks would alter the geometry, but this happens on a century to millennia time scale. So we expected a systematic pattern in parameters variability, which should be linked to the system’s dynamics. We identified that more of the variability found resulted from varying length of recessions and that was why we gave recommendation (L450-460) to consider grouping recession events by length in order to identify actual variability resulting from system’s heterogeneity. Our main concern for providing such suggestion is for users (e.g. water resource managers) who ideally would rely on specific values for resource planning and management.

Therefore, I would encourage the authors to report the seasonal variation of the estimated parameters because it is a valuable insight gained through the individual segments approach and, moreover, a legitimate motivation in favor of the approach. They should capitalize on it and promote this feature as an additional strength of the approach. Furthermore, no paper, to my knowledge, reports the seasonal variability of η and ε for karst systems, same for i and K. For now, the importance of seasonality is barely mentioned despite that it is recognized as a source of variability (L435-437), and its contribution is not pictured or quantified. I understand that the paper is focused on proposing a method, but a graph about the seasonal dynamics of parameters or the classification is not expansive and could be placed in the supplementary materials.

Reply to SI-1 (Paragraph 3)

The length of recession is directly related to duration/intensity of precipitation event preceding the recession event. This means, longer recession in summer and shorter one in winter. In fact, differing length is actually influence of seasonal variation. Hence, we completely agree with the reviewer recommendation to explore the parameter variability along seasonal dimension. Therefore, we will provide additional analysis regarding parameters and seasonal variability and we will review lines L25-27, L450-460, and L525-526. Additional discoveries from this analysis will be included in the appendices as well.

 SI-2: Motivation for the Mangin model and framework.

The hypotheses behind the Mangin model are not discussed enough, and the classification framework is not sufficiently described in section 2.2. The hypothesis of the duality of porosity and processes in karst systems is well described. However, the hypothesis of the matrix flow following the linear recession equation of Maillet (L170, Eq. 3) is not. Nonlinear equations are common to describe baseflow, including in karst systems [7]. Some comments on why the Maillet approximation can be used would have been appreciated, if not theoretical, with an empirical justification of its application showing that Mangin's framework is widely used and will benefit from more robustness.

Reply to SI-2 (Paragraph 1)

We will expand the discussion on the Mangin model in section 2.2 and add a new subsection to specifically introduce the Mangin’s classification framework.

Besides, the authors paradoxically rely on the Askoy & Wittenberg (A&W) REM that explicitly states that baseflow is nonlinear (L141, Eq. 1). At first sight, it appears as a rationale problem. However, if the importance of the classification Mangin framework is more emphasized in the objective of the paper and better described in section 2.2 (or best, in its own section), using the linear model that is part of this classification framework of the authors' choice make sense. Also, even if the Mangin baseflow model is linear, extracting recession with the nonlinear method of A&W could be justifiable, because Mangin's model with both components and sufficient degrees of freedom, is in fact, a nonlinear recession model. Still, it would have been more consistent, in my opinion, to use A&W approach with the linear Maillet Eq.3, eventually, by choosing a higher CV to compensate the lack of degrees of freedom and the fear that the outcome would be too restrictive on the number of recession point.

Reply to SI-2 (Paragraph 2)

Once again, we thank the reviewer for pointing this out. It is correct that A&W extracts baseflow recession points by fitting non-linear model. However, Eq.1 (L141) becomes linear if the value of b = 1(Aksoy and Wittenberg, 2011). In fact, we set b to 1 to make sure we have that consistency with the linear approach of Mangin. We will add one or two sentences to paragraph 3 of section 2.1 to clarify this.

 SI-3: Motivation/Discussion of the REMs

In the discussion (L411-412), the authors mention that: "Overall, the adapted REMs and the introduced three POAs provide range of results that adequately represents the karst systems. This suggests that the modified REMs are well suited for application to karst spring recession analysis". I found the conclusion a little too optimistic and leaving little room for improvement. Plus, I feel that it is contradictory with the suggestion of focusing on different segment lengths to reduce uncertainty (L455-56). Why not say that REMs are unsuited then?

Three REMs are selected on the honorable basis of what previous authors have suggested (L160, Table 1). Despite my concerns about the A&W REMs (SI-2), it is quite common to refer to and implement these former REMs in recent papers [6, and 8-9 as referred by the authors]. Thus, it is justifiable. However, this is not a strong rationale. None of the authors of REMs had the primary objective of providing a method for recession extraction. They had to, but their focus was on recession analysis, not extraction. Also, they were not aware of the recent discovery inference problems associated with the extraction method (e.g., [8]).  In fact, they look at recession from the narrow and specific scope of the criteria set in Table 1. Using them broaden the scope of the analysis, but still, these are three specific angles of approach to recession extraction, and it is delicate to affirm that they adequately represent the karst systems.

My point is that a noticeable perspective for improvement would be to rely on one or more generic method that allows varying the criteria of selections (beginning/end of recession and tolerance to anomalies). It would allow inferring the statistics for ranges of criteria instead of a few predetermined ones. Coupling a generic method to the POAs could certainly improve the framework in the future.

Reply to SI-3

We thank the reviewer for the valuable suggestion. It is true that the REMs were developed to suit catchments with specific properties. We also accept that there is potential to improvement and further development of karst-specific recession extraction method. We will review the statements of L411-412 accordingly and discuss aspects of recession extraction that could benefit from further improvement. However, we would like to mention that recession selection criteria used by the REMs could still be modified by a user to suit specific requirement.

 SI-4: Degrees of freedom and data points

There is a concern about the parsimony of the approach and the limited number of points in the recession segments [6]. Instead of three models, one for each study site, more than two thousand models were fitted, in which the degrees of freedom is very close to the number of data points. So far, I believe that the results and their interpretation make some sense, and so, I am willing to believe that what is captured is meaningful.

I could have recommended relying on a nonlinear two parameters equation to add another degree of freedom to the late recession (the nonlinear exponent, SI-2), but should we? The natural question is how far we could go in terms of degrees of freedom and when equifinality will produce unreliable meaningless parameters? I understand that this is a separate question that would deserve another paper, but it is worth mentioning it.

Reply to SI-4 (Paragraph 1)

We thank the reviewer for providing a very useful and intriguing insight. In terms of recession data points, we selected recession events longer or equal to 7 days period (L284-286), so we have enough data points for the parameters optimization. The whole analysis (recession extraction and parameter optimization) is not computationally expensive; it took couple of minutes to run the analysis in R.

We also considered introducing a fourth Parameter Optimization Approach (POA) by a using nonlinear model for the slow flow recession. However, this would result in too many dimensions of REMs and POAs; also, it would be inconsistent with the Mangin classification framework we used. Nevertheless, we will be exploring this in a different study.

Also, I miss some more statistics in Table 4 (L292), such as the average number of recession points per event, which would help appreciate this issue. Also, perhaps statistics about the flow Q and dQ/dt captured by the REM methods would also be interesting.

Reply to SI-4 (Paragraph 2)

We thank you for pointing this out. We will expand Table 4 by providing the suggested statistics of the extracted recession events.

 SI-5: Uncertainties and the subsequent suggestions

I have experienced contradictory opinions about the uncertainties. On the one hand, they are considered essential because they reflect the actual dynamics of karst. On the other hand, the authors wish to reduce them and offer suggestions for doing so. Also, attention should be paid to using "reliable parameters" when we know that recession models are underfitting recession dynamics (SI-1). The term "robust parameters" seems more appropriate.

Reply to SI-5 (Paragraph 1)

Our motivation for using IRS analysis is to capture the karst system dynamics reflected by spring discharge (Jachens et al., 2020; Kovács et al., 2005). As discussed in our reply to specific comment 1 (Paragraph 1&2), we found the variability of parameters estimated with the IRS analysis to be too large and we believe this can be reduced in a systematic way or better quantify. We address issues relating to our suggestion for reducing the uncertainty in our reply to subsequent paragraphs where the reviewer further discusses them. We will replace “reliable parameters” with “robust parameters” in the manuscript.

The authors refer to suggestions in the abstract and conclusion, but nothing is said explicitly about them in these sections. The suggestions are two, not very well evidenced and hidden in lines L455-460 : (i) focusing on segments of different recession lengths, (ii) using the master recession curve (MRC). I found none of them to be very relevant nor developed. With one (i), it is basically said that REMs are inappropriate (SI-3), and, by focusing on one type of segment length, some processes are dismissed, and you reduce uncertainty by occulting the natural variability that you aim at capturing. Perhaps, presenting the results per season would be more appropriate (SI-1) to reduce the uncertainty of sensitivity to initial conditions.

Reply to SI-5 (Paragraph 2)

As discussed in our reply to Paragraph 3 of specific comment 1, our suggestion to group recession events by length is encompassed by considering seasonal variability. We agree that exploring the parameters variability with respect to seasonality will provide better approach to explain/account for the parameters uncertainty. Therefore, we will implement this change as we already mention in our reply to specific comment 1.

Regarding the use of MRC (ii), I am not sure I understand what the authors mean. I understand that since the MRC is an average behavior, fitting this average will produce a single estimate of the parameters, which will naturally reduce the uncertainty. Well, this is not reducing, it is still occulting, which is why I believe the authors ultimately reject their own suggestion in the last sentence (L459-60). Note that the MRC approaches with uncertainty quantification do exist (e.g.  [10]), so what is said in L99-101 and L459-60 is not correct. It is entirely possible to consider an approach to MRC that takes variability into account. That said, I would not recommend it. I don't see the point of building an MRC, which is a delicate process when you have a model that you can calibrate directly on the whole hydrograph. I understand that some of the references cited do, but it doesn't make much sense to me. The MRC, in my opinion, is an empirical method that is applied when one wants to abstain from model fitting because of their too restrictive and inappropriate hypotheses on the dynamic (see the introduction of [11]). The authors' approach is hypothetical-deductive and relies on intelligible parameters to classify the dynamics. The references to the empirical MRC are, in my opinion, inappropriate and confusing.

Reply to SI-5 (Paragraph 3)

Thank you for pointing out that there are MRC approaches with uncertainty quantification. However, we do not know of any MRC analysis approach for multi-segment (quick and slow flow) recession event that quantifies uncertainty. Having mention this, we will review the statements made in L99-101 and L459-60. Our suggestion of MRC in L457-458 read thus “Another way of coping with this problem is to consider master recession curve analysis which is often criticize for its inability to adequately represent storage variability (Gregor and Malík, 2012; Kovacs, 2021; Kresic and Bonacci, 2010)”. By this statement, we do not intend to mislead the reader that MRC approach is better. We only provide this suggestion as an alternative approach if one does not want to be worried with issues of parameters variability arising from IRS analysis. However, we understand the concerns of the reviewer and we will rework L457-458 to reflect our main intention for providing the suggestion.

In L459-460, we quote “However, since per event analysis is useful for better understanding of the system’s dynamic, defining a systematic approach to quantify parameter uncertainty will help to increase the confidence of individual recession segment analysis”. With improved modeling approaches and computational ability available nowadays, we believe a systematic uncertainty quantification can be incorporated into IRS analysis. We are presently working on a separate study that integrate parameter uncertainty with IRS and if possible, we will be glad to share this with the reviewer in the future.

I think better suggestions would most likely come from criticizing the approach, the Mangin model, and the classification framework. If one projects data from a complex, high dimensional reality onto a lower-dimensional model and classification framework, there will be uncertainties that are both natural variability and projection artifacts. If they are too large, the message is that the model or framework is not helpful, and another projection must be found. Mangin developed his framework without taking uncertainties into account. Now that the authors have done so, relevant questions arise: should the dispersion of the parameter distribution, which reflects the sensitivity of the watershed, be taken into account in the watershed classification framework? How can this be done? Furthermore, the i indicator was found to be a poor discriminant for the classification of the three study sites. Is this a poor indicator that should be removed from a classification framework, or can we expect that with other study sites, i may show more distinct and clear patterns? Should a nonlinear recession be considered in place of the Maillet equation?  These are just a few interesting questions that arise when the authors allow themselves to critique the analytical framework.

Questioning it is not shooting oneself in the foot. I think the Mangin model and framework is a recognized and intelligible way of framing the dynamics of the karst system, and the authors have succeeded in providing a more robust way of doing this. The authors could stress this success while recognizing that they have also highlighted potential limits. Even so, the proposed automated framework could also be used more widely to assess the relevance of the Mangin model and classification framework in the future to propose relevant improvements. Finally, their approach is also transferable to other models and classification frameworks that will be developed in the future.

Reply to SI-5 (Paragraph 4 and 5)

We thank the reviewer for the interesting and intriguing questions raised about the Mangin classification framework. In fact few authors (e.g. El-Hakim and Bakalowicz, 2007) have also identified some limitations of the Mangin framework and modified the classification scheme. Almost in every studies where Mangin classification has been used, estimation of K and i were based on few recession events, usually the longest summer event(s). So the robustness of the classification framework has not really been tested with a more dynamic sample of recession events as we have done in this study. Therefore, with respect to the suggestion of the reviewer, we will provide additional discussion to highlight the questions raised about the classification and possibility for improvements.

 SI-7: Reproducibility

This issue is not related to the introduction or the discussion section but to the lack of supplementary materials. The authors propose a technique. This technique is itself motivated by the idea of reducing user error and bias. Although the mathematics remains clear and straightforward, providing input/output files, code recipes, or examples will significantly help future users ensure that the technique is correctly applied and, therefore, be more in phase the paper's purpose.

Reply to SI-7

As previously mentioned in our reply to the general comments, we will provide the R-code and input datasets with usage instructions via our GitHub repository.

 Comment on the title

The title is straightforward. However, "efficient" and "accurate" are not the best term.

Efficient can mean many things, for instance, "computationally efficient" or well organized. Accurate is not the best given that the point is to depict and leverage the parameters' variability. The fact that the method is automated is stressed in the introduction and is probably more appropriate. The central classification task is missing in the title. I would suggest, for instance: "Karst spring recession analysis and classification: an automated method considering both fast and slow flow components."

Reply to comment on the title

We thank the reviewer for the suggestion. We considered the suggestions made by the reviewer. However, we will like to retain “efficient” because the methods are fast and computationally inexpensive. Therefore, we will change the title to “Karst spring recession analysis and classification: efficient, automated methods for both fast and slow components”.

References

Aksoy, H. and Wittenberg, H.: Nonlinear baseflow recession analysis in watersheds with intermittent streamflow, Hydrol. Sci. J., 56(2), 226–237, doi:10.1080/02626667.2011.553614, 2011.

El-Hakim, M. and Bakalowicz, M.: Significance and origin of very large regulating power of some karst aquifers in the Middle East. Implication on karst aquifer classification, J. Hydrol., 333(2–4), 329–339, doi:10.1016/j.jhydrol.2006.09.003, 2007.

Gregor, M. and Malík, P.: Construction of master recession curve using genetic algorithms, J. Hydrol. Hydromechanics, 60(1), 3–15, doi:10.2478/v10098-012-0001-8, 2012.

Jachens, E. R., Rupp, D. E., Roques, C. and Selker, J. S.: Recession analysis revisited: impacts of climate on parameter estimation, Hydrol. Earth Syst. Sci, 24, 1159–1170, doi:10.5194/hess-24-1159-2020, 2020.

Jeannin, P.-Y. and Sauter, M.: Analysis of karst hydrodynamic behaviour using global approaches: a review., 1998.

Kovacs, A.: Quantitative classification of carbonate aquifers based on hydrodynamic behaviour, , (February), doi:10.1007/s10040-020-02285-w, 2021.

Kovács, A., Perrochet, P., Király, L. and Jeannin, P. Y.: A quantitative method for the characterisation of karst aquifers based on spring hydrograph analysis, J. Hydrol., 303(1–4), 152–164, doi:10.1016/j.jhydrol.2004.08.023, 2005.

Kresic, N. and Bonacci, O.: Spring discharge hydrograph., 2010.

Technical issues

Overall, the reading is understandable and straightforward as one feels the fear of not being understood correctly. However, as a result, the reading is sometimes wordy (lots of conjunctions) and redundant, with punctuation errors, past/present tenses getting mixed up, plural/singular mistakes, many missing articles, and other grammar errors. I recommend a few more iterations to correct, simplify and shorten the sentences. Accordingly, I have not reported all the errors and have gradually focused on the minor problems of interpretation that were not part of the SI sections.

Reply to Technical issues

We thank the reviewer for pointing out these mistakes/errors. We will diligently review the manuscript and rework all complex sentences with complex structure.

• Abstract

L21: Mixing past and present. Missing word: … parameters [that ?] reasonably represent …

Reply: This sentence will reviewed.

L22: The average KGE seems by far higher than 0.7 (see figure 3)

Reply: The average KGE is actually >0.9, we will replace 0.7 with the appropriate value.

L25-27: On the recommendation for reducing uncertainties, see SP5 (Same remarks for

the last sentence of the conclusion L525-526).

Reply: We will make the appropriate revision to L25-27.

• Introduction

L40: missing word: time series [analysis?]

Reply: We will include the missing word.

L43: a lumped attributes for karst systems

Reply: We will review this sentence.

L49: describes

Reply: We will review this sentence.

L69: I don't like the term "manually". It gives an impression of a "pen and paper" approach. I guess that recession fitting is still computationally assisted. I would replace manually with a more precise word such as visually supervised or validated by the user. If changed, run a find and replace to spot the word elsewhere in the paper.

Reply: We will replace “manually” to “visually supervised” and anywhere else it appears on the paper.

L81: we can modify them

Reply: We will review this sentence.

L85-99: the two first bullet items are redundant with the first sentence. I suggest removing the bullet list and keep the objective and method summarized in one paragraph.

Reply: We will review this sentences as suggested by the reviewer.

L99-101: the? Master recession curve, the analysis. However, I recommend removing this last sentence. The motivation for fitting individual recession should be (and is already succinctly) discussed before so the last paragraph is entirely dedicated to what is shown in the paper.

Reply: We will remove this sentence as suggested by the reviewer.

• Data and Methods

L104: "generic", I would remove. These methods seem rather specific to me. They are called specific later on (L155).

Reply: We will modify the text accordingly.

L116: Do you apply the recession fitting on the smoothed points?

Reply: We fitted the recession model to the observed Q (non-smoothed points)

L120: "baseflow (non-influenced)" use the consistent term slow flow defined above (L108-110). Otherwise, it is confusing.

Reply: We will modify text and use consistent terminology.

L123: a more strict

Reply: We will modify the text accordingly.

L128: Stoelzlz --> Stoelzle

Reply: We will modify the text accordingly.

L133: the Askoy and Wittenberg iterative procedure is not described. How iterative? By truncating progressively the head of the segment?

Reply: If coefficient of variation (CV) for fitting of a flow recession is larger than 0.10, then the first day of the recession period is excluded and the remaining part is fitted again. If CV is still higher than 0.10, one more day is excluded and so forth until CV is lower or equal to 0.10.

L147: rule based and exclusion criteria --> rules and exclusion criteria

Reply: We will modify the text accordingly.

L148: I would use "aim at eliminating …" instead of "ensure". By adding degrees of freedom, you acknowledge that they actually don't ensure anything.

Reply: We will modify the text accordingly.

L170-176-184: symbols can be simplified by using Q + subscript more often instead of Q, q, phi, psy

Reply: We will consider changing the symbols used in equations.

L173: is referenced without being defined.

Reply: We do not understand what this comment means.

L188: requires + remove comma

Reply: We will modify the text accordingly.

L189: Eq.5 --> Eq.4

Reply: We will modify the text accordingly.

L192-193. The last sentence comes from nowhere. Put this last sentence below, after introducing K. Alternatively, you can start a new section 2.3 called "Karst classification".

Reply: We will add a new subsection for karst system classification.

L197: what five classes? I would definitely provide a more detailed description. Eventually say that the classification is illustrated later on with Figure 5.

Reply: We will elaborate this under the new “karst system classification” subsection.

L203: I don't understand the definition of i with respect to Eq. 8. Should "the value of the quickflow recession" be "the proportion of remaining quickflow…"?

Reply: Infiltration delay, i, defines the retardation effect between infiltration and output (discharge). The value of t in equation is set to 2 (Jeannin and Sauter, 1998) for comparing rate of infiltration through saturated zone. We will give more explanation under the new subsection “Karst system classification”.

L208: move in the first paragraph.

Reply: We will modify the text accordingly.

L211: The figure is very clear but, perhaps, an illustration of the effect of η and ε would be nice as well.

Reply: We will consider the reviewer suggestion.

L214-215: I will add some references about the inference of recession parameters from the MRC. Jachens et al. seems important enough to be placed in the introduction (see SI-1).

Reply: We will include the reference in the introduction.

L221: "which is in reality untrue": I would not speak of reality or truthfulness. I would rather say that REMs are possibly not identifying the most relevant or optimal t_i.

Reply: We will modify the text accordingly.

L231: for clarity, "a predefined value" -->"the REMs predefined values".

Reply: We will modify the text accordingly.

L259: I would briefly precise why KGE is chosen (over Nash and Sutcliffe Efficiency, for instance).

Reply: We will provide justification for choosing KGE to determine model efficiency.

L261-2: The Mangin classification is actually not described in section 2.2, only K and i. See comment above (L192-193). I would remove the last two sentences if the classification is correctly described in the appropriate section.

Reply: as previously mentioned, we will add a new subsection for karst system classification and text will be modified according.

L262: This is used

Reply: We will modify the text accordingly.

• Test springs and data

L279 / Table 3 / Temporal res. : Since a reader first looks at the figures and tables in an article, specifying the time resolution as "hourly" and "sub-hourly" may mislead the reader into thinking that the methods have been applied at the time resolutions mentioned.

Reply: This is the temporal resolution of original dataset. All datasets were aggregated to daily time step for the analysis (L273-75).

• Results

L292/Table 4: I miss some important statistics about the length of recession events, Q, and -dQ/dt. I believe that they could help in the interpretation of the results. If not, consider presenting these statistics in the supplementary materials.

Reply: We will additional statistics of extracted recession events in appendices/supplement materials.

L312: The modification … produces … interactions. To simplify …

Reply: We will modify the text accordingly.

L314: Remove comma + a defined category

Reply: We will modify the text accordingly.

L316: The results show

Reply: We will modify the text accordingly.

L317: The estimation or estimations

Reply: We will modify the text accordingly.

L320-324: several missing articles. A much …

Reply: We will add missing articles in these lines.

L326: The infiltration rate, the curve concavity, …

Reply: We will modify the text accordingly.

L329: "vice-versa" is confusing.

Reply: We will review the sentence for better understanding.

L331-335: I was frustrated about the patterns not being explained, but they are in the discussion section. Therefore, I suggest referring to the discussion section.

Reply: We will modify the text accordingly.

L337: The estimation .. the curve concavity …

Reply: We will modify the text accordingly.

L347: What about Q0? It is an important parameter to evaluate aquifer recharge and to perform hydrograph separation. I would recommend adding it or putting it in the supplementary materials. Besides, you mention the Q0 compensation effect later on (L466).

Reply: We will provide a plot of Q0 vs REMs and POAs for all springs in the appendices.

L357: “The grouped mean” is confusing. Grouped how? Are the boxes min-max ranges?

Reply: Based on the Mangin classification framework, we derived per-event classification with each pair of REM and POA. Then we estimated the mean (average of all event classes) classification for each REM-POA pair. The colour-filled boxes are drawn over the area covered by the mean classes derived with all pair of REM and POA for each spring.

L417-460: I miss a discussion about the role of spatio-temporal variability of rainfall and the full set of other factors affecting recession sensitivity to initial conditions

Reply: As mentioned in our reply to specific comment 1, we will provide additional analysis on the influence of seasonality on parameters variability and relevant discussion (as highlighted here by the reviewer) will be provided.

L431: I think that less permissive REMs will tend to include the early recession, hence, increase the recession coefficient alpha if the Maillet equation is used alone because the recession rate would be, on average, more pronounced when high flows are considered. However, I am not sure that this is true with the entire Mangin model that separates high from low recession flows, and I do not personally see in Figure 4 that "permissive recession extraction tends to produce higher estimates of slow flow recession coefficients."

Reply: Permissive extraction methods (Vogel and Brutsaert) produced slightly higher value of alpha for Saivu and Qachquoch spring (Figure 4). Although this is not seen in Lehnbachquellen spring which is located in a snow-dominated catchment. However, we believe that our additional analysis on seasonality influence on parameter variability will give us a clearer result about this.

L441-442: what does “under” or “over-estimated” mean since there is no reference?

Reply: We will change “under” to “low” and “over” to “high”

L443: higher infiltration rates

Reply: We will modify the text accordingly.

L486: a recent study … categorises

Reply: We will modify the text accordingly.

L487: calculated

Reply: We will modify the text accordingly.

L491: methods involve

Reply: We will modify the text accordingly.

L658: Table A1 in the supplementary materials is unused and unexplained.

Reply: We will replace Table A1 and provide explanation for the new one.

Reply to comments by anonymous referee (Reviewer 2) on the manuscript “Karst spring recession curve analysis: efficient, accurate methods for both fast and slow flow components” by Tunde Olarinoye et al., Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hess-2021-249-RC1, 2021

General comments

Thank you for the opportunity to review this manuscript " Karst spring recession curve analysis: efficient, accurate methods for both fast and slow flow components." The authors have presented an intriguing combination of techniques to separate and quantify fast and slow flow components of a spring hydrograph recession.  This is an exciting and valuable contribution to the literature and, while the methodology appears technically sound, there are a few overarching issues that I feel need to be addressed prior to final publication.

Reply to general comments

We thank reviewer 2 for her/his useful and constructive comments. The following suggestions made by the reviewer will be taken into account and following changes will be made accordingly.

First, the use of a snow-dominated system violates much of the underlying assumptions of recession curve analysis. Namely, that the recession behavior is controlled by infiltration and flow path properties within the rock itself. In a snowmelt system, this is not often the case - snowpack often acts as a separate storage component in the system that is not controlled by aquifer properties but instead by a combination of sublimation, melt, and refreeze processes. To continue utilizing this site, the authors should be clear up front about this potential violation of assumptions and utilize this in the discussion of the results.

Reply: We thank the reviewer for pointing out issues about other external controls in a snow-dominated catchment. We will clearly highlight this issue and a potential violation of assumptions in section 3 (Test springs and data).

Similarly, it would be valuable to discuss the potential implications of seasonality of precipitation (and soil moisture) on infiltration capacity in the mediterranean system. Either more information on when recession curves are in relation to seasonal precipitation or discussion of this impact on the fit of the results would be valuable.

Reply: We thank the reviewer for her/his suggestion. Similar to our reply to first reviewer’s specific comments 1 and 5, we will provide additional analysis and discussion on the effects of seasonality on the variability of recession parameters

Generally, the document could use some minor grammatical review to correct some minor issues throughout.

Reply: We will carefully review the entire text to correct all grammatical and spelling errors.

General comments

Thank you for the opportunity to review this manuscript " Karst spring recession curve analysis: efficient, accurate methods for both fast and slow flow components." The authors have presented an intriguing combination of techniques to separate and quantify fast and slow flow components of a spring hydrograph recession.  This is an exciting and valuable contribution to the literature and, while the methodology appears technically sound, there are a few overarching issues that I feel need to be addressed prior to final publication.

Reply to general comments

We thank reviewer 2 for her/his useful and constructive comments. The following suggestions made by the reviewer will be taken into account and following changes will be made accordingly.

First, the use of a snow-dominated system violates much of the underlying assumptions of recession curve analysis. Namely, that the recession behavior is controlled by infiltration and flow path properties within the rock itself. In a snowmelt system, this is not often the case - snowpack often acts as a separate storage component in the system that is not controlled by aquifer properties but instead by a combination of sublimation, melt, and refreeze processes. To continue utilizing this site, the authors should be clear up front about this potential violation of assumptions and utilize this in the discussion of the results.

Reply: We thank the reviewer for pointing out issues about other external controls in a snow-dominated catchment. We will clearly highlight this issue and a potential violation of assumptions in section 3 (Test springs and data).

Similarly, it would be valuable to discuss the potential implications of seasonality of precipitation (and soil moisture) on infiltration capacity in the mediterranean system. Either more information on when recession curves are in relation to seasonal precipitation or discussion of this impact on the fit of the results would be valuable.

Reply: We thank the reviewer for her/his suggestion. Similar to our reply to first reviewer’s specific comments 1 and 5, we will provide additional analysis and discussion on the effects of seasonality on the variability of recession parameters

Generally, the document could use some minor grammatical review to correct some minor issues throughout.

Reply: We will carefully review the entire text to correct all grammatical and spelling errors.